Fifth China Wen
Fifth China Wen
Research Article
Invariant Solutions and Conservation Laws of the
(2 + 1)-Dimensional Boussinesq Equation
Copyright © 2014 Wenjuan Rui et al. This is an open access article distributed under the Creative Commons Attribution License,
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Invariant solutions and conservation laws of the (2 + 1)-dimensional Boussinesq equation are studied. The Lie symmetry approach is
used to obtain the invariant solutions. Conservation laws for the underlying equation are derived by utilizing the new conservation
theorem and the partial Lagrange approach.
obtained many explicit exact solutions of (1) by using the One can write the Lie-Bäcklund operator (7) in characteristic
new generalized transformation in homogeneous balance form as follows:
method. More new double periodic and multiple soliton 𝜕 𝜕
solutions are obtained for the generalized (2 + 1)-dimensional 𝑋 = 𝜉𝑖 𝐷𝑖 + 𝑊𝛼 𝛼 + ∑𝐷𝑖1 𝑖2 ⋅⋅⋅𝑖𝑠 (𝑊𝛼 ) 𝛼 . (10)
𝜕𝑢 𝑠≥1 𝜕𝑢 𝑖1 𝑖2 ⋅⋅⋅𝑖𝑠
Boussinesq equation [20].
The main goal of this paper is to use the Lie symmetry The Noether operators associated with a Lie-Bäcklund sym-
method [21, 22] to obtain the invariant solutions. In addition metry operator 𝑋 are given by
to this, conservation laws will be derived for (1) by using the
𝛿 𝛿
new conservation theorem [23, 24] and the partial Lagrange 𝑁𝑖 = 𝜉𝑖 + 𝑊𝛼 + ∑𝐷 (𝑊𝛼 ) 𝛼 ,
approach [25, 26]. 𝛿𝑢𝑖𝛼 𝑠≥1 𝑖1 𝑖2 ⋅⋅⋅𝑖𝑠 𝛿𝑢𝑖𝑖1 𝑖2 ⋅⋅⋅𝑖𝑠 (11)
The outline of this paper is as follows. In Section 2,
we present symmetry group analysis and group-invariant 𝑖 = 1, 2, . . . , 𝑛,
solutions of (1). In Section 3, the conservation laws for (1) are where the Euler-Lagrange operators with respect to deriva-
established. Finally, some conclusions are given in Section 4. tives of 𝑢𝛼 are obtained from (5) by replacing 𝑢𝛼 by the
corresponding derivatives. For example,
2. Method of Lie Symmetries 𝛿 𝜕 𝜕
= + ∑(−1)𝑠 𝐷𝑗1 ⋅ ⋅ ⋅ 𝐷𝑗𝑠 𝛼 ,
2.1. Preliminaries. In this section we briefly present the 𝛿𝑢𝑖𝛼 𝜕𝑢𝑖𝛼 𝑠≥1 𝜕𝑢𝑖𝑗1 𝑗2 ⋅⋅⋅𝑗𝑠 (12)
notation and pertinent results used in this paper [21, 22].
𝑖 = 1, . . . , 𝑛, 𝛼 = 1, . . . , 𝑚,
Consider a 𝑘th-order system of PDEs of 𝑛 independent
variables 𝑥 = (𝑥1 , 𝑥2 , . . . , 𝑥𝑛 ) and 𝑚 dependent variables 𝑢 = and the Euler-Lagrange, Lie-Bäcklund, and Noether opera-
(𝑢1 , 𝑢2 , . . . , 𝑢𝑚 ): tors are connected by the operator identity:
𝐸𝛼 (𝑥, 𝑢, 𝑢(1) , . . . , 𝑢(𝑘) ) = 0, 𝛼 = 1, . . . , 𝑚, 𝛿
(3) 𝑋 + 𝐷𝑖 (𝜉𝑖 ) = 𝑊𝛼 + 𝐷𝑖 𝑁𝑖 . (13)
𝛿𝑢𝛼
where 𝑢(1) , 𝑢(2) , . . . , 𝑢(𝑘) denote the collections of all first,
The 𝑛-tuple vector 𝑇 = (𝑇1 , 𝑇2 , . . . , 𝑇𝑛 ), 𝑇𝑗 ∈ A, 𝑗 = 1, . . . , 𝑛
second,. . ., 𝑘th-order partial derivatives; that is, 𝑢𝑖𝛼 =
𝐷𝑖 (𝑢𝛼 ), 𝑢𝑖𝑗𝛼 = 𝐷𝑗 𝐷𝑖 (𝑢𝛼 ), . . ., respectively, with the total is a conserved vector of (3) if 𝑇𝑖 satisfies
derivative operator with respect to 𝑥𝑖 given by 𝐷𝑖 𝑇𝑖 |(3) = 0. (14)
where the additional coefficients are determined uniquely by 𝜉𝑥𝑥 = 𝜂𝑥𝑥 = 𝜂𝑥𝑢 = 𝜂𝑢𝑢 = 0,
the prolongation formulae: 𝜏𝑡 − 𝜙𝑦 = 0, 𝜏𝑡 − 2𝜉𝑥 = 0, 𝜙𝑡 − 𝜏𝑦 = 0,
𝜁𝑖𝛼 = 𝐷𝑖 (𝑊 ) + 𝛼
𝜉𝑗 𝑢𝑖𝑗𝛼 , 𝜏𝑡𝑡 + 2𝜂𝑢𝑡 + 𝜏𝑦𝑦 = 0, 𝜙𝑡𝑡 − 2𝜂𝑢𝑡 + 𝜙𝑦𝑦 = 0, (17)
(8)
𝜁𝑖𝛼1 ⋅⋅⋅𝑖𝑠 = 𝐷𝑖1 ⋅ ⋅ ⋅ 𝐷𝑖𝑠 (𝑊𝛼 ) + 𝜉𝑗 𝑢𝑗𝑖
𝛼
1 ⋅⋅⋅𝑗𝑠
, 𝑠 > 1, 𝜂𝑥𝑥 + 𝜂𝑦𝑦 − 𝜂𝑡𝑡 + 𝜂𝑥𝑥𝑥𝑥 = 0,
Solving the above system (17) we reach that the symmetry of 3.1. Variational Method for a System and Its Adjoint. The
(1) is spanned by the five vector fields: system of adjoint equations to the system of 𝑘th-order
differential equation (3) is defined by [23]
𝜕 𝜕 𝜕
𝑋1 = , 𝑋2 = , 𝑋3 = ,
𝜕𝑥 𝜕𝑦 𝜕𝑡 𝐸𝛼∗ (𝑥, 𝑢, V, . . . , 𝑢(𝑘) , V(𝑘) ) = 0, 𝛼 = 1, . . . , 𝑚, (27)
𝜕 𝜕
𝑋4 = 𝑦 +𝑡 , (18) where
𝜕𝑡 𝜕𝑦
𝛿 (]𝛽 𝐸𝛽 )
𝜕 𝜕 𝜕 𝜕 𝐸𝛼∗ (𝑥, 𝑢, V, . . . , 𝑢(𝑘) , V(𝑘) ) = ,
𝑋5 = 𝑥 + 2𝑦 + 2𝑡 + (2𝑢 + 1) . 𝛿𝑢𝛼 (28)
𝜕𝑥 𝜕𝑦 𝜕𝑡 𝜕𝑢
𝛼 = 1, . . . , 𝑚, ] = ] (𝑥) ,
2.3. Group-Invariant Solutions. In this section we present the
reduction forms of (1) by using symmetry group method. and ] = (]1 , ]2 , . . . , ]𝑚 ) are new dependent variables.
To do this, particular linear combinations of infinitesimals We recall here the following results as given in Ibragimov
are considered and their corresponding invariants are deter- [23].
mined. For example, for the symmetry operator 𝑋 = 𝑋1 + A system of (3) is said to be self-adjoint if the substitution
𝑋2 + 𝑋3 , we can compute the invariants by integrating the of ] = 𝑢 into the system of adjoint equation (27) yields the
characteristic equations: same system (3).
Assume the system of (3) admits the symmetry generator:
𝑑𝑥 𝑑𝑦 𝑑𝑡
= = . (19)
1 1 1 𝜕 𝜕
𝑋 = 𝜉𝑖 𝑖
+ 𝜂𝛼 𝛼 , (29)
The corresponding invariants are 𝑞 = 𝑥 − 𝑦, 𝑝 = 𝑡 − 𝑦, 𝜔 = 𝜕𝑥 𝜕𝑢
𝑢. Now treating 𝑝, 𝑞 as the new independent variables and
and then the system of adjoint equation (27) admits the
𝜔 as the new dependent variable, we obtain the ordinary
operator:
differential equation:
2 𝜕 𝜕 𝜕
(2𝜔 (𝑞) + 1) 𝜔 (𝑞) + 2(𝜔 (𝑞)) + 𝜔(4) (𝑞) = 0. (20) 𝑌 = 𝜉𝑖 + 𝜂𝛼 𝛼 + 𝜂∗𝛼 𝛼 ,
𝜕𝑥𝑖 𝜕𝑢 𝜕] (30)
Integrating (20) twice with respect to 𝑞 we obtain 𝜂∗𝛼 = − [𝜆𝛼𝛽 ]𝛽 𝛼 𝑖
+ ] 𝐷𝑖 (𝜉 )] ,
2
𝜔 (𝑞) + 𝜔 (𝑞) + 𝜔 (𝑞) = 𝑐0 𝑞 + 𝑐1 , 𝑐0 , 𝑐1 ∈ 𝑅. (21)
where the operator (30) is an extension of (29) to the variable
Taking 𝑐0 = 0 and 𝜔 (𝑞) = 𝜓(𝜔) we get ]𝛼 and the 𝜆𝛼𝛽 are obtainable from
If there exist a function 𝐿 = 𝐿(𝑥, 𝑢, 𝑢(1) , . . . , 𝑢(𝑙) ) ∈ A, 𝑙 ≤ 𝑘 where V = V(𝑥, 𝑦, 𝑡) is a new dependent variable and (40) gives
and nonzero functions 𝑓𝛼𝛽 ∈ A such that (34) can be written
V𝑡𝑡 − 2𝑢𝑥𝑥 V − 4𝑢𝑥 V𝑥 − (2𝑢 + 1) V𝑥𝑥 − V𝑦𝑦 − V𝑥𝑥𝑥𝑥 = 0. (41)
as 𝛿𝐿/𝛿𝑢𝛼 = 𝑓𝛼𝛽 𝐸𝛽1 , provided 𝐸𝛽1 ≠ 0, then 𝐿 is known
as a partial Lagrangian of (34); otherwise it is the standard By using Theorem 1, we obtain the following Lagrangian for
Lagrangian. The differential equations of the form the system of (38) and (40):
𝛿𝐿
= 𝑓𝛼𝛽 𝐸𝛽1 (35) 𝐿 = V (𝑢𝑡𝑡 − (2𝑢 + 1) 𝑢𝑥𝑥 − 2𝑢𝑥2 − 𝑢𝑦𝑦 − 𝑢𝑥𝑥𝑥𝑥 ) . (42)
𝛿𝑢𝛼
are called a system of partial Euler-Lagrange equations.
The operator 𝑋 in (10) is a partial Noether operator
corresponding to a partial Lagrangian 𝐿 ∈ A of the system (1) We first consider the Lie point symmetry generator
(35) if it can be determined from 𝑋1 = 𝜕/𝜕𝑥, and we have 𝑊 = −𝑢𝑥 . Hence using (32),
we obtain the following components of the conserved
𝛿𝐿 vector 𝑇:
𝑋 (𝐿) + 𝐿𝐷𝑖 (𝜉𝑖 ) = 𝑊𝛼 + 𝐷𝑖 (𝐵𝑖 ) , (36)
𝛿𝑢𝛼
𝑇𝑥1 = (2𝑢 + 2V − 1) 𝑢𝑥 V𝑥 − 𝑢𝑥𝑥 V𝑥𝑥 − 𝑢𝑥 V𝑥𝑥𝑥 ,
1 2 𝑛 𝑖
for some vector 𝐵 = (𝐵 , 𝐵 , . . . , 𝐵 ), 𝐵 ∈ A. Here 𝑊 =
(𝑊1 , 𝑊2 , . . . , 𝑊𝑚 ), 𝑊𝛼 ∈ A is the characteristic of 𝑋. 𝑇𝑡1 = − 𝑢𝑥𝑡 V + 𝑢𝑥 V𝑡 , (43)
3.3. Conservation Law. We now construct conservation law (3) Now the Lie point symmetry generator 𝑋3 = 𝜕/𝜕𝑡 has
of (1) using the two approaches outlined above. the Lie characteristic function 𝑊 = −𝑢𝑡 . Hence we
can obtain the conserved vector whose components
Application of the New Conservation Theorem. The (2 + 1)- are
dimensional Boussinesq equation is given by
𝑇𝑥3 = [2𝑢𝑥 𝑢𝑡 + (2𝑢 + 1) 𝑢𝑥𝑡 + 𝑢𝑥𝑥𝑥𝑡 ] V − (2𝑢 + 1) 𝑢𝑡 V𝑥
𝐸𝛼 = 𝑢𝑡𝑡 − (2𝑢 + 1) 𝑢𝑥𝑥 − 2𝑢𝑥2 − 𝑢𝑦𝑦 − 𝑢𝑥𝑥𝑥𝑥 = 0. (38) − 𝑢𝑥𝑡 V𝑥𝑥 − 𝑢𝑡 V𝑥𝑥𝑥 ,
We recall that (38) admits the following five Lie point
symmetry generators: 𝑇𝑡3 = − [2𝑢𝑥2 + (2𝑢 + 1) 𝑢𝑥𝑥 + 𝑢𝑦𝑦 ] V + 𝑢𝑡 V𝑡 − VV𝑥𝑥𝑥𝑥 ,
𝜕 𝜕 𝑇𝑦3 = − 𝑢𝑦𝑡 V − 𝑢𝑡 V𝑦 .
𝑋1 = , 𝑋2 = ,
𝜕𝑥 𝜕𝑦 (45)
𝜕 𝜕 𝜕
𝑋3 = , 𝑋4 = 𝑦 +𝑡 , (39) (4) The Lie point symmetry generator 𝑋4 = 𝑦(𝜕/𝜕𝑡) +
𝜕𝑡 𝜕𝑡 𝜕𝑦 𝑡(𝜕/𝜕𝑦) has the Lie characteristic function 𝑊 =
𝜕 𝜕 𝜕 𝜕 −𝑡𝑢𝑦 − 𝑦𝑢𝑡 . Hence using (32), the components of the
𝑋5 = 𝑥 + 2𝑦 + 2𝑡 + (2𝑢 + 1) . conserved vector are given by
𝜕𝑥 𝜕𝑦 𝜕𝑡 𝜕𝑢
The adjoint equation of (38), by invoking (27), is 𝑇𝑥4 = [2𝑢𝑥 (𝑡𝑢𝑦 + 𝑦𝑢𝑡 ) + (2𝑢 + 1) (𝑡𝑢𝑥𝑦 + 𝑦𝑢𝑥𝑡 )
𝛿 [V (𝑢𝑡𝑡 − (2𝑢 + 1) 𝑢𝑥𝑥 − 2𝑢𝑥2 − 𝑢𝑦𝑦 − 𝑢𝑥𝑥𝑥𝑥 )] + 𝑡𝑢𝑥𝑥𝑥𝑦 + 𝑦𝑢𝑥𝑥𝑥𝑡 ] V − (2𝑢 + 1) (𝑡𝑢𝑦 + 𝑦𝑢𝑡 ) V𝑥
𝐸𝛼∗ = = 0,
𝛿𝑢
(40) − (𝑡𝑢𝑥𝑦 + 𝑦𝑢𝑥𝑡 ) V𝑥𝑥 − (𝑡𝑢𝑦 + 𝑦𝑢𝑡 ) V𝑥𝑥𝑥 ,
Abstract and Applied Analysis 5
𝑇𝑡4 = − [𝑢𝑦 + 2𝑦𝑢𝑥2 + (2𝑢 + 1) 𝑦𝑢𝑥𝑥 + 𝑦𝑢𝑦𝑦 to determine the Noether-type operators. If we substitute (48)
and (50) into the partial Noether-type operators determining
+ 𝑡𝑢𝑦𝑡 + 𝑦𝑢𝑥𝑥𝑥𝑥 ] V + (𝑡𝑢𝑦 + 𝑦𝑢𝑡 ) V𝑡 , (36), we obtain
𝜉 = 𝜙 = 𝜏 = 𝜂𝑢 = 0, 𝜂 = 𝐶 (𝑥, 𝑦, 𝑡) ,
𝑇𝑦4 = [𝑢𝑡 − 2𝑡𝑢𝑥2 + 𝑡𝑢𝑡𝑡 − (2𝑢 + 1) 𝑡𝑢𝑥𝑥 + 𝑦𝑢𝑦𝑡 − 𝑡𝑢𝑥𝑥𝑥𝑥 ] V
𝐵1 = − 𝐶𝑡 𝑢 + 𝛼 (𝑥, 𝑦, 𝑡) ,
− (𝑡𝑢𝑦 + 𝑦𝑢𝑡 ) V𝑦 . (51)
(46) 𝐵2 = − (𝑢2 + 𝑢) 𝐶𝑥 + 𝛽 (𝑥, 𝑦, 𝑡) ,
1 1 1 1 2
𝐿 = 𝑢𝑡2 − (𝑢 + ) 𝑢𝑥2 − 𝑢𝑦2 − 𝑢𝑥𝑥 . (48) References
2 2 2 2
[1] M. J. Ablowitz and P. A. Clarkson, Solitons, Nonlinear Evolution
The Euler-Lagrange-type equation is Equations and Inverse Scattering, Cambridge University Press,
New York, NY, USA, 1991.
𝛿𝐿
= −𝑢𝑡𝑡 + (2𝑢 + 1) 𝑢𝑥𝑥 + 𝑢𝑥2 + 𝑢𝑦𝑦 + 𝑢𝑥𝑥𝑥𝑥 . (49) [2] V. B. Matveev and M. A. Salle, Darboux Transformations and
𝛿𝑢 Solitons, Springer, Berlin, Germany, 1991.
So (1) can be written as [3] Y. Chen, Z. Y. Yan, and H. Q. Zhang, “Exact solutions for a
family of variable-coefficient reaction—duffing equations via
𝛿𝐿 the Backlund transformation,” Theoretical and Mathematical
= −𝑢𝑥2 , (50) Physics, vol. 132, no. 1, pp. 970–975, 2002.
𝛿𝑢
6 Abstract and Applied Analysis
[4] A. Wazwaz, “The Hirota’s direct method for multiple-soliton [22] P. J. Olver, Applications of Lie Groups to Differential Equations,
solutions for three model equations of shallow water waves,” vol. 107 of Graduate Texts in Mathematics, Springer, Berlin,
Applied Mathematics and Computation, vol. 201, no. 1-2, pp. Germany, 2nd edition, 1993.
489–503, 2008. [23] N. H. Ibragimov, “A new conservation theorem,” Journal of
[5] E. Fan and H. Zhang, “A note on the homogeneous balance Mathematical Analysis and Applications, vol. 333, no. 1, pp. 311–
method,” Physics Letters A, vol. 246, no. 5, pp. 403–406, 1998. 328, 2007.
[6] E. Fan, “Two new applications of the homogeneous balance [24] N. H. Ibragimov and T. Kolsrud, “Lagrangian approach to evo-
method,” Physics Letters A, vol. 265, no. 5-6, pp. 353–357, 2000. lution equations: symmetries and conservation laws,” Nonlinear
[7] S. Y. Lou, H. Y. Ruan, D. F. Chen, and W. Z. Chen, “Similarity Dynamics, vol. 36, no. 1, pp. 29–40, 2004.
reductions of the KP equation by a direct method,” Journal of [25] A. H. Kara and F. M. Mahomed, “Relationship between symme-
Physics A: Mathematical and General, vol. 24, no. 7, pp. 1455– tries and conservation laws,” International Journal of Theoretical
1467, 1991. Physics, vol. 39, no. 1, pp. 23–40, 2000.
[8] S. Lou, X. Tang, and J. Lin, “Similarity and conditional similarity [26] A. H. Kara and F. M. Mahomed, “Noether-type symmetries and
reductions of a (2 + 1)-dimensional KdV equation via a direct conservation laws via partial Lagrangians,” Nonlinear Dynam-
method,” Journal of Mathematical Physics, vol. 41, no. 12, pp. ics, vol. 45, no. 3-4, pp. 367–383, 2006.
8286–8303, 2000.
[9] A. Wazwaz, “The tanh method for traveling wave solutions of
nonlinear equations,” Applied Mathematics and Computation,
vol. 154, no. 3, pp. 713–723, 2004.
[10] E. Yusufoglu and A. Bekir, “Solitons and periodic solutions of
coupled nonlinear evolution equations by using the sine-cosine
method,” International Journal of Computer Mathematics, vol.
83, no. 12, pp. 915–924, 2006.
[11] P. E. Hydon, Symmetry Methods for Differential Equations: A
Beginners Guide, Cambridge University Press, Cambridge, UK,
2000.
[12] F. M. Mahomed, “Symmetry group classification of ordinary
differential equations: survey of some results,” Mathematical
Methods in the Applied Sciences, vol. 30, no. 16, pp. 1995–2012,
2007.
[13] R. Hirota, “Direct methods in soliton theory,” in Solitons, R. K.
Bullogh and P. J. Caudrey, Eds., Springer, Berlin, Germany, 1980.
[14] R. Hirota, “Exact envelope-soliton solutions of a nonlinear wave
equation,” Journal of Mathematical Physics, vol. 14, pp. 805–809,
1973.
[15] E. V. Krishnan, S. Kumar, and A. Biswas, “Solitons and
other nonlinear waves of the Boussinesq equation,” Nonlinear
Dynamics, vol. 70, no. 2, pp. 1213–1221, 2012.
[16] H. W. Yang, B. S. Yin, and Y. L. Shi, “Forced dissipative
Boussinesq equation for solitary waves excited by unstable
topography,” Nonlinear Dynamics, vol. 70, no. 2, pp. 1389–1396,
2012.
[17] S. M. El-Sayed and D. Kaya, “The decomposition method
for solving (2+1)-dimensional Boussinesq equation and (3+1)-
dimensional KP equation,” Applied Mathematics and Computa-
tion, vol. 157, no. 2, pp. 523–534, 2004.
[18] M. Senthilvelan, “On the extended applications of homoge-
neous balance method,” Applied Mathematics and Computation,
vol. 123, no. 3, pp. 381–388, 2001.
[19] Y. Chen, Z. Yan, and H. Zhang, “New explicit solitary wave
solutions for (2 + 1)-dimensional Boussinesq equation and (3 +
1)-dimensional KP equation,” Physics Letters A, vol. 307, no. 2-3,
pp. 107–113, 2003.
[20] C. Huai-Tang and Z. Hong-Qing, “New double periodic and
multiple soliton solutions of the generalized (2+1)-dimensional
Boussinesq equation,” Chaos, Solitons and Fractals, vol. 20, no.
4, pp. 765–769, 2004.
[21] G. W. Bluman and S. Kumei, Symmetries and Differential
Equations, vol. 81 of Applied Mathematical Sciences, Springer,
New York, NY, USA, 1989.
Advances in Advances in Journal of Journal of
Operations Research
Hindawi Publishing Corporation
Decision Sciences
Hindawi Publishing Corporation
Applied Mathematics
Hindawi Publishing Corporation
Algebra
Hindawi Publishing Corporation
Probability and Statistics
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
International
Journal of Journal of
Mathematics and
Mathematical
Discrete Mathematics
Sciences