Strongly damped wave equations are of a great current interest, see Refs. [6, 8–10, 13–15, 18, 20, 24, 25] and references therein. In  an initial-boundary value problem for a nonlinear wave equation in one space dimension with a nonlinear damping term |u|m–1ut and a source term of the form |u|p–1u with m, p > 1 has been considered. The authors analyzed the global well posedness of the initial-boundary value problem corresponding to
It is well known [10, 18, 20, 21] that when the damping term |u|m–1ut is absent from the equation, then the source term |u|p–1u yields the solution to blow up in finite time. However, the interaction between the damping and source terms is often difficult to analyze. In  a new method of investigating the so-called quasilinear strongly-damped wave equations in bounded 3D domains has been presented. This method allows us to establish the existence and uniqueness of energy solutions in some cases for the growth exponent of the non-linearity and f. In  for the strongly damped nonlinear wave equation with Dirichlet boundary conditions, where both the non-linearities f and g exhibit a critical growth, while h is a time-independent forcing term the existence of an exponential attractor of optimal regularity is proven.
It is known that conservation laws play a significant role in the solution process of an equation or a system of differential equations. Although not all of the conservation laws of partial differential equations (PDEs) may have physical interpretation they are essential in studying the integrability of the PDEs.
In [2–5], Anco and Bluman gave a general algorithmic method to find all conservations laws for evolution equations Many recent papers using this method have been published in [1, 11, 12]. After Ibragimov’s results several papers appeared concerned with self-adjointness and its applications to PDEs [22, 23]. In this work, we derive conservation laws by using the direct method of the multipliers [2–5]. In this paper, we consider the strongly damped wave equations of the form (2) where f(u) and g(u) are arbitrary non-zero functions and a > 0 is a given dissipation parameter. We provide a complete classification of all point symmetries and all low-order conservation laws admitted by this equation (2).
2 Classical symmetries
Lie classical method is based on the determination of the symmetry group of a differential equation, i.e., the largest group of transformations acting on dependent and independent variables of the equation so that maps solutions of the equation into other solutions.
In order to apply Lie classical method to equation (2) we consider the one-parameter Lie group of infinitesimal transformations in (x, t, u) given by (3) (4) (5) where ∊ is the group parameter. The symmetry group of equation (2) will be given by the set of vector fields of the form (6)
with J = (j1, . . ., jk), 1 ≤ jk ≤ 2 and 1 ≤ k ≤ 3, and u(3) denotes the sets of partial derivatives up to third order . Applying (7) we obtain a set of determining equations for the infinitesimals ξ(x, t, u), τ(x, t, u) and η(x, t, u). By solving the determining system we obtain the following results.
The point symmetries admitted by a strongly damped wave equation (2) in the general case are generated by
If f(u) = α and g(u) = βu + γ, with α, β and γ arbitrary constants, besides v1 and v2 we obtain 2 new generators: where q(x, t) must satisfy equation (2).
If with γ is an arbitrary constant, besides v1, v2 and v4, we obtain the following generator:
3 Symmetry reductions and exact solutions
In this section, we use the Lie point symmetries admitted by (2) to obtain group-invariant solutions, then by using the simplest equation method due to Kudriashov we derive some exact solutions for (2).
In order to find travelling wave solutions for equation (2) we consider that f(u) and g(u) are arbitrary functions, by using the generator X = λ∂x + ∂t we obtain the similarity variables and similarity solution (8) and the following reduced ODE (9)
In  Kudryashov introduced the simplest equation method in order to look for exact solutions of nonlinear differential equations. In  it was pointed out that many of the so called “new travelling wave solutions” could be derived from the solutions of a simple nonlinear ordinary differential equation. In this paper looking for travelling waves solutions of (2) and observing that (2) can be written in conserved form for f(u) linear and g(u) quadratic we set and we search for exact solutions of (9) with
We are now considering as simplest equation (10)
The general solution of equation (10) with a1 ≠ 0 can be written in terms of the Jacobi elliptic functions sn(kz, p).
The general solution of equation (10) with a1 = 0 can be written in terms of the Jacobi elliptic function sn2(kz, p) or via de Weierstrass elliptic function. We now look for exact solutions h(z) of (11) such that h(z) satisfies (10). By substituting and into (11) we obtain that the following conditions must be satisfied (12)
We obtain the following:
As a particular case for equation (11) admits the solution h = sech2(z) and (2) admit a soliton solution (13)
In Figure 1, we plot solution u(x, t) (13), which describes a soliton.
For f(u) = α and g(u) = ßu+γ, taking by using generator
We get the invariant solution (14) where h(z) must satisfy (15)
Finally, for by using generator v5 the corresponding invariant solution is written in the form (16)
where h(z) must satisfy (17)
Integrating once with respect to z we obtain (18)
where c0 is a constant of integration. For this reduction, ODE (18) is a second order equation whose solution is
where c1, c2 are arbitrary constants and E represents the error function. Therefore, undoing transformation (16), we obtain a solution of equation (2).
4 Conservation laws
A local conservation law of equation (2) is a continuity equation (19) that holds for the whole set of solutions u(x, t), where the conserved density C1 and the spatial flux C2 are functions of x, t, u, and derivatives of u. Here Dt, Dx denote total derivatives with respect to t and x respectively. The pair of expressions (C1, C2) is called a conserved current.
Two local conservation laws are considered to be equivalent  if they differ by a trivial conservation law C1 = DXΘ, C2 = –DtΘ, where C1 and C2 are evaluated on the set of solutions of equation (2), and Θ is some function of x, t, u, and derivatives of u.
We begin by observing that equation (2) has a Cauchy-Kovaleskaya form. Consequently, the results in Ref.[3, 4] show that all non-trivial conservation laws arise from multipliers. Specifically, when we move off of the set of solutions of equation (2), every non-trivial local conservation law (19) is equivalent to one that can be expressed in the characteristic form (20)
where Λ(x, t, u, ux, ut, ...) is a multiplier, and where (C̃1, C̃2) differs from (C1, C2) by a trivial conserved current. On the set of solutions u(x, t) of equation (2), the characteristic form (20) reduces to the conservation law (19).
In general, a function Λ(x, t, u, ux, ut, ...) is a multiplier if it is non-singular on the set of solutions u(x, t) of equation (2), and if its product with equation (2) is a divergence expression with respect to t and x. There is a one-to-one correspondence between non-trivial multipliers and non-trivial conservation laws in characteristic form.
The determining equation to obtain all multipliers is (21)
This equation must hold off of the set of solutions of equation (2). Once the multipliers are found, the corresponding non-trivial conservation laws are obtained either by using a homotopy formula Ref. [2–4] or by integrating the characteristic equation (20) Ref. .
We will now find all multipliers and obtain corresponding non-trivial conservation laws. The determining equation (21) splits with respect to the variables utt, uttt, uttx, utxx, uxxx, utttx, uttxx, utxxx, uxxxx. This yields a linear determining system for Λ(x, t, u, ut, ux, utx, uxx) which can be solved by the same algorithmic method used to solve the determining equation for infinitesimal symmetries. By using Maple we solve this determining system subject to the classification conditions. After solving the determining system we do not obtain any multiplier for general f(u) and g(u)
All special cases for which additional multipliers (2) are admitted consist of: (22a) (22b) (22c) (23a) (23b) (24a) (24b) (24c) (24d) (24e)
For each solution Q, a corresponding conserved density and flux can be derived by integration of the divergence identity (20) [5, 7]. The admitted conservation laws in each case are given by: (25a) (25b) (26a) (26b) (27a) (27b) (28a) (28b) (29a) (29b) (30a) (30b) (31a) (31b)
In this paper, we have considered a strongly damped wave equations with two arbitrary non-zero functions f(u) and g(u) and a positive constant a. We have provided a complete classification of all point symmetries and all low-order conservation laws admitted by this equation. Equation (9) when f(u) is a linear function and g(u) is a quadratic function has abundant exact solutions that can be expressed in terms of the Jacobi elliptic functions and consequently the equation (9) has abundant exact solutions that can be expressed in terms of trigonometric and hyperbolic functions. Hence (2) has a plenty of periodic waves, solitary waves, etc. These solutions are derived from the solutions of a simple nonlinear ordinary differential equation.
Adem K.R., Khalique C.M., Symmetry analysis and conservation laws of a generalized two-dimensional nonlinear KP-MEW equation, Mathematical Problems in Engineering, 2015, 2015, Article ID 805763. Web of ScienceGoogle Scholar
Anco S.C., Bluman G., Direct construction method for conservation laws of partial differential equations Part I: Examples of conservation law classifications, Eur. J. Appl. Math., 2002, 13, 545–566. Google Scholar
Anco S.C., Bluman G., Direct construction method for conservation laws for partial differential equations Part II: General treatment, Euro. Jnl of Applied Mathematics, 2002, 41, 567–585. Google Scholar
Anco S.C., Generalization of Noether’s theorem in modern form to non-variational partial differential equations, To appear in Fields Institute Communications: Recent progress and Modern Challenges in Applied Mathematics, Modeling and Computational Science. arXiv: mathph/1605.08734 Google Scholar
Bluman G.W., Cheviakov A., Anco S.C., Applications of symmetry methods to partial differential equations, New York: Springer, 2009. Google Scholar
Gandarias M.L., Khalique CM., Symmetries solutions and conservation laws of a class of nonlinear dispersive wave equations, Communications in Nonlinear Science and Numerical Simulation, 2016, 32, 114–121. Web of ScienceCrossrefGoogle Scholar
Kalantarov V.K., Attractors for some nonlinear problems of mathematical physics, Zap. Nauchn. Sem. Leningrad. Otdel. Mat. Inst. Steklov (LOMI), 1986, 152, 50–54. Google Scholar
Olver P., Applications of Lie groups to differential equations, Springer-Verlag, 1993. Google Scholar
Tracinà R., Bruzón M.S., Gandarias M.L., On the nonlinear self-adjointness of a class of fourth-order evolution equations, Applied Mathematics and Computation, 2016, 275, 299–304. CrossrefWeb of ScienceGoogle Scholar
Tracinà R., Freire I. L., Torrisi M., Nonlinear self-adjointness of a class of third order nonlinear dispersive equations, Communications in Nonlinear Science and Numerical Simulation, 2016, 32, 225–233. CrossrefWeb of ScienceGoogle Scholar
Yang M., Sun C., Dynamics of strongly damped wave equations in locally uniform spaces: Attractors and asymptotic regularity, Trans. Am. Math. Soc., 2009, 361(2), 1069–1101. Google Scholar
About the article
Published Online: 2017-05-04
Citation Information: Open Physics, Volume 15, Issue 1, Pages 300–305, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2017-0033.
© 2017 M. Luz Gandarias et al.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. BY-NC-ND 3.0