# Invariant Solutions and Conservation Laws of the (2 + 1)-Dimensional Boussinesq Equation

Authors: Wenjuan Rui; Peiyi Zhao; Yufeng Zhang
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

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.


---

## Body

## 1. Introduction

In recent years, searching for explicit solutions of nonlinear evolution equations (NEEs) has attracted the attention of many mathematicians and physicists. Particularly, various effective methods have been used to explore different kinds of solutions of NEEs, such as the inverse scattering method [1], the Darboux transformation [2] and the Bäcklund transformation [3], the Hirota method [4], the homogeneous balance method [5, 6], the similarity reduced method [7, 8], the tanh method [9], and the sine-cosine method [10]. But up to now a unified method that can be used to deal with all types of NEEs has not been discovered. Among the above mentioned methods, the Lie symmetry method is one of the most effective methods to determine solutions of differential equations. In the past decades, there have been considerable developments in symmetry methods for differential equations [11, 12].In the present paper we consider the (2 + 1)-dimensional Boussinesq equation:(1)
u
t
t
-
u
x
x
-
u
y
y
-
(
u
2
)
x
x
-
u
x
x
x
x
=
0
.
A well-known soliton equation is the Boussinesq equation:
(2)
u
t
t
-
u
x
x
-
3
(
u
2
)
x
x
-
u
x
x
x
x
=
0
.
Equation (2) was introduced by Boussinesq to describe motions of long waves in shallow water [13, 14]. It also appears in a wide variety of physical systems such as nonlinear lattice waves, iron sound waves in plasma, and vibrations in a nonlinear string. For the transonic speed perturbations, by neglecting the interaction of waves moving in the opposite directions, the Boussinesq equation (2) can be reduced to the KdV equation. Equation (2) itself is also a dimensional reduction of the KP equation in the moving frame. Moreover, the Boussinesq equation (2) is completely integrable and admits inverse scattering. Owing to its profound importance and nice mathematical properties, a great deal of research work has been invested in recent years for the study of the Boussinesq equation. Krishnan et al. [15] studied the dynamics of shallow water waves that are governed by the Boussinesq equations. Yang et al. [16] obtained solutions of homogeneous and inhomogeneous dissipative Boussinesq equation by using the modified Jacobi elliptic function expansion method and the pseudospectral method. However, there has also been a growing interest in the study of NEEs in multidimensions, especially in (1 + 2) and (1 + 3) dimensions. To find some exact soliton solutions in higher dimensions is much more difficult than in 1 + 1 dimensions. Recently, El-Sayed and Kaya [17] used the decomposition method to obtain the exact solitary-wave solutions of (1). Senthilvelan [18] obtained the travelling wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation by homogeneous balance method and explored certain new solutions of the equations. Chen et al. [19] obtained many explicit exact solutions of (1) by using the new generalized transformation in homogeneous balance method. More new double periodic and multiple soliton solutions are obtained for the generalized (2 + 1)-dimensional Boussinesq equation [20].The main goal of this paper is to use the Lie symmetry method [21, 22] to obtain the invariant solutions. In addition to this, conservation laws will be derived for (1) by using the new conservation theorem [23, 24] and the partial Lagrange approach [25, 26].The outline of this paper is as follows. In Section2, we present symmetry group analysis and group-invariant solutions of (1). In Section 3, the conservation laws for (1) are established. Finally, some conclusions are given in Section 4.


## 2. Method of Lie Symmetries

### 2.1. Preliminaries

In this section we briefly present the notation and pertinent results used in this paper [21, 22].Consider akth-order system of PDEs of n independent variables x
=
(
x
1
,
x
2
,
…
,
x
n
) and m dependent variables u
=
(
u
1
,
u
2
,
…
,
u
m
):
(3)
E
α
(
x
,
u
,
u
(
1
)
,
…
,
u
(
k
)
)
=
0
,
α
=
1
,
…
,
m
,
where u
(
1
)
,
u
(
2
)
,
…
,
u
(
k
) denote the collections of all first, second,…, kth-order partial derivatives; that is, u
i
α
=
D
i
(
u
α
)
,
u
i
j
α
=
D
j
D
i
(
u
α
)
,
…, respectively, with the total derivative operator with respect to x
i given by
(4)
D
i
=
∂
∂
x
i
+
u
i
α
∂
∂
u
α
+
u
i
j
α
∂
∂
u
j
α
+
⋯
,
i
=
1
,
…
,
n
,
where the summation convention is used whenever appropriate.The Euler-Lagrange operator, for eachα, is given by
(5)
δ
δ
u
α
=
∂
∂
u
α
+
∑
s
≥
1
‍
(
-
1
)
s
D
i
1
⋯
D
i
s
∂
∂
u
i
1
i
2
⋯
i
s
α
,
and the Lie-Bäcklund operator is
(6)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
,
ξ
i
,
η
α
∈
A
,
where A is the space of differential functions. The operator (6) is an abbreviated form of infinite formal sum:
(7)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
+
∑
s
≥
1
‍
ζ
i
1
i
2
⋯
i
s
α
∂
∂
u
i
1
i
2
⋯
i
s
α
,
where the additional coefficients are determined uniquely by the prolongation formulae:
(8)
ζ
i
α
=
D
i
(
W
α
)
+
ξ
j
u
i
j
α
,
ζ
i
1
·
·
·
i
s
α
=
D
i
1
⋯
D
i
s
(
W
α
)
+
ξ
j
u
j
i
1
⋯
j
s
α
,
s
>
1
,
in which W
α is the Lie characteristic function:
(9)
W
α
=
η
α
-
ξ
j
u
j
α
.
One can write the Lie-Bäcklund operator (7) in characteristic form as follows:
(10)
X
=
ξ
i
D
i
+
W
α
∂
∂
u
α
+
∑
s
≥
1
‍
D
i
1
i
2
⋯
i
s
(
W
α
)
∂
∂
u
i
1
i
2
⋯
i
s
α
.
The Noether operators associated with a Lie-Bäcklund symmetry operator X are given by
(11)
N
i
=
ξ
i
+
W
α
δ
δ
u
i
α
+
∑
s
≥
1
‍
D
i
1
i
2
⋯
i
s
(
W
α
)
δ
δ
u
i
i
1
i
2
⋯
i
s
α
,
i
=
1,2
,
…
,
n
,
where the Euler-Lagrange operators with respect to derivatives of u
α are obtained from (5) by replacing u
α by the corresponding derivatives. For example,
(12)
δ
δ
u
i
α
=
∂
∂
u
i
α
+
∑
s
≥
1
‍
(
-
1
)
s
D
j
1
⋯
D
j
s
∂
∂
u
i
j
1
j
2
⋯
j
s
α
,
i
=
1
,
…
,
n
,
α
=
1
,
…
,
m
,
and the Euler-Lagrange, Lie-Bäcklund, and Noether operators are connected by the operator identity:
(13)
X
+
D
i
(
ξ
i
)
=
W
α
δ
δ
u
α
+
D
i
N
i
.
The n-tuple vector T
=
(
T
1
,
T
2
,
…
,
T
n
)
,
T
j
∈
A
,
j
=
1
,
…
,
n is a conserved vector of (3) if T
i satisfies
(14)
D
i
T
i
|
(
3
)
=
0
.
Equation (14) defines a local conservation law of system (3).


### 2.2. Lie Point Symmetries

Now, we consider the following Lie group of transformations with independent variablesx
,
y
,
t and dependent variable u:
(15)
x
*
=
x
*
(
x
,
y
,
t
,
ε
)
,
y
*
=
y
*
(
x
,
y
,
t
,
ε
)
,
t
*
=
t
*
(
x
,
y
,
t
,
ε
)
,
u
*
=
u
*
(
x
,
y
,
t
,
ε
)
,
where ε is the group parameter. The infinitesimal generator for the Lie group transformations can be expressed in the following form:
(16)
X
=
ξ
∂
∂
x
+
ϕ
∂
∂
y
+
τ
∂
∂
t
+
η
∂
∂
u
.
Applying the fourth prolongation p
r
(
4
)
X to (1), we obtain the following determining equations:
(17)
ξ
y
=
ξ
t
=
ξ
u
=
ϕ
x
=
ϕ
u
=
τ
x
=
τ
u
=
0
,
ξ
x
x
=
η
x
x
=
η
x
u
=
η
u
u
=
0
,
τ
t
-
ϕ
y
=
0
,
τ
t
-
2
ξ
x
=
0
,
ϕ
t
-
τ
y
=
0
,
τ
t
t
+
2
η
u
t
+
τ
y
y
=
0
,
ϕ
t
t
-
2
η
u
t
+
ϕ
y
y
=
0
,
η
x
x
+
η
y
y
-
η
t
t
+
η
x
x
x
x
=
0
,
η
-
ξ
x
+
τ
t
+
2
u
(
τ
t
-
ξ
x
)
=
0
,
η
u
+
2
τ
t
-
2
ξ
x
=
0
,
4
η
x
+
ξ
t
t
+
2
η
x
u
-
ξ
y
y
=
0
.
Solving the above system (17) we reach that the symmetry of (1) is spanned by the five vector fields:
(18)
X
1
=
∂
∂
x
,
X
2
=
∂
∂
y
,
X
3
=
∂
∂
t
,
X
4
=
y
∂
∂
t
+
t
∂
∂
y
,
X
5
=
x
∂
∂
x
+
2
y
∂
∂
y
+
2
t
∂
∂
t
+
(
2
u
+
1
)
∂
∂
u
.


### 2.3. Group-Invariant Solutions

In this section we present the reduction forms of (1) by using symmetry group method. To do this, particular linear combinations of infinitesimals are considered and their corresponding invariants are determined. For example, for the symmetry operator X
=
X
1
+
X
2
+
X
3, we can compute the invariants by integrating the characteristic equations:
(19)
d
x
1
=
d
y
1
=
d
t
1
.
The corresponding invariants are q
=
x
-
y
,
p
=
t
-
y
,
ω
=
u. Now treating p
,
q as the new independent variables and ω as the new dependent variable, we obtain the ordinary differential equation:
(20)
(
2
ω
(
q
)
+
1
)
ω
′′
(
q
)
+
2
(
ω
′
(
q
)
)
2
+
ω
(
4
)
(
q
)
=
0
.
Integrating (20) twice with respect to q we obtain
(21)
ω
′′
(
q
)
+
ω
2
(
q
)
+
ω
(
q
)
=
c
0
q
+
c
1
,
c
0
,
c
1
∈
R
.
Taking c
0
=
0 and ω
′
(
q
)
=
ψ
(
ω
) we get
(22)
ψ
′
(
ω
)
ψ
(
ω
)
+
ω
2
+
ω
=
c
1
.
The solution of (22) is
(23)
ψ
(
ω
)
=
±
1
3
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
,
where c
2 is a constant. Notice that ω
′
(
q
)
=
ψ
(
ω
); we have
(24)
ω
′
(
q
)
=
±
1
3
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
,
so we get
(25)
q
=
±
3
∫
‍
1
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
d
ω
.
Recall q
=
x
-
y
,
ω
=
u; we get the solution of (1):
(26)
x
-
y
=
±
3
∫
‍
1
-
6
u
3
+
18
c
1
u
-
9
u
2
+
9
c
2
d
u
.


## 2.1. Preliminaries

In this section we briefly present the notation and pertinent results used in this paper [21, 22].Consider akth-order system of PDEs of n independent variables x
=
(
x
1
,
x
2
,
…
,
x
n
) and m dependent variables u
=
(
u
1
,
u
2
,
…
,
u
m
):
(3)
E
α
(
x
,
u
,
u
(
1
)
,
…
,
u
(
k
)
)
=
0
,
α
=
1
,
…
,
m
,
where u
(
1
)
,
u
(
2
)
,
…
,
u
(
k
) denote the collections of all first, second,…, kth-order partial derivatives; that is, u
i
α
=
D
i
(
u
α
)
,
u
i
j
α
=
D
j
D
i
(
u
α
)
,
…, respectively, with the total derivative operator with respect to x
i given by
(4)
D
i
=
∂
∂
x
i
+
u
i
α
∂
∂
u
α
+
u
i
j
α
∂
∂
u
j
α
+
⋯
,
i
=
1
,
…
,
n
,
where the summation convention is used whenever appropriate.The Euler-Lagrange operator, for eachα, is given by
(5)
δ
δ
u
α
=
∂
∂
u
α
+
∑
s
≥
1
‍
(
-
1
)
s
D
i
1
⋯
D
i
s
∂
∂
u
i
1
i
2
⋯
i
s
α
,
and the Lie-Bäcklund operator is
(6)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
,
ξ
i
,
η
α
∈
A
,
where A is the space of differential functions. The operator (6) is an abbreviated form of infinite formal sum:
(7)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
+
∑
s
≥
1
‍
ζ
i
1
i
2
⋯
i
s
α
∂
∂
u
i
1
i
2
⋯
i
s
α
,
where the additional coefficients are determined uniquely by the prolongation formulae:
(8)
ζ
i
α
=
D
i
(
W
α
)
+
ξ
j
u
i
j
α
,
ζ
i
1
·
·
·
i
s
α
=
D
i
1
⋯
D
i
s
(
W
α
)
+
ξ
j
u
j
i
1
⋯
j
s
α
,
s
>
1
,
in which W
α is the Lie characteristic function:
(9)
W
α
=
η
α
-
ξ
j
u
j
α
.
One can write the Lie-Bäcklund operator (7) in characteristic form as follows:
(10)
X
=
ξ
i
D
i
+
W
α
∂
∂
u
α
+
∑
s
≥
1
‍
D
i
1
i
2
⋯
i
s
(
W
α
)
∂
∂
u
i
1
i
2
⋯
i
s
α
.
The Noether operators associated with a Lie-Bäcklund symmetry operator X are given by
(11)
N
i
=
ξ
i
+
W
α
δ
δ
u
i
α
+
∑
s
≥
1
‍
D
i
1
i
2
⋯
i
s
(
W
α
)
δ
δ
u
i
i
1
i
2
⋯
i
s
α
,
i
=
1,2
,
…
,
n
,
where the Euler-Lagrange operators with respect to derivatives of u
α are obtained from (5) by replacing u
α by the corresponding derivatives. For example,
(12)
δ
δ
u
i
α
=
∂
∂
u
i
α
+
∑
s
≥
1
‍
(
-
1
)
s
D
j
1
⋯
D
j
s
∂
∂
u
i
j
1
j
2
⋯
j
s
α
,
i
=
1
,
…
,
n
,
α
=
1
,
…
,
m
,
and the Euler-Lagrange, Lie-Bäcklund, and Noether operators are connected by the operator identity:
(13)
X
+
D
i
(
ξ
i
)
=
W
α
δ
δ
u
α
+
D
i
N
i
.
The n-tuple vector T
=
(
T
1
,
T
2
,
…
,
T
n
)
,
T
j
∈
A
,
j
=
1
,
…
,
n is a conserved vector of (3) if T
i satisfies
(14)
D
i
T
i
|
(
3
)
=
0
.
Equation (14) defines a local conservation law of system (3).


## 2.2. Lie Point Symmetries

Now, we consider the following Lie group of transformations with independent variablesx
,
y
,
t and dependent variable u:
(15)
x
*
=
x
*
(
x
,
y
,
t
,
ε
)
,
y
*
=
y
*
(
x
,
y
,
t
,
ε
)
,
t
*
=
t
*
(
x
,
y
,
t
,
ε
)
,
u
*
=
u
*
(
x
,
y
,
t
,
ε
)
,
where ε is the group parameter. The infinitesimal generator for the Lie group transformations can be expressed in the following form:
(16)
X
=
ξ
∂
∂
x
+
ϕ
∂
∂
y
+
τ
∂
∂
t
+
η
∂
∂
u
.
Applying the fourth prolongation p
r
(
4
)
X to (1), we obtain the following determining equations:
(17)
ξ
y
=
ξ
t
=
ξ
u
=
ϕ
x
=
ϕ
u
=
τ
x
=
τ
u
=
0
,
ξ
x
x
=
η
x
x
=
η
x
u
=
η
u
u
=
0
,
τ
t
-
ϕ
y
=
0
,
τ
t
-
2
ξ
x
=
0
,
ϕ
t
-
τ
y
=
0
,
τ
t
t
+
2
η
u
t
+
τ
y
y
=
0
,
ϕ
t
t
-
2
η
u
t
+
ϕ
y
y
=
0
,
η
x
x
+
η
y
y
-
η
t
t
+
η
x
x
x
x
=
0
,
η
-
ξ
x
+
τ
t
+
2
u
(
τ
t
-
ξ
x
)
=
0
,
η
u
+
2
τ
t
-
2
ξ
x
=
0
,
4
η
x
+
ξ
t
t
+
2
η
x
u
-
ξ
y
y
=
0
.
Solving the above system (17) we reach that the symmetry of (1) is spanned by the five vector fields:
(18)
X
1
=
∂
∂
x
,
X
2
=
∂
∂
y
,
X
3
=
∂
∂
t
,
X
4
=
y
∂
∂
t
+
t
∂
∂
y
,
X
5
=
x
∂
∂
x
+
2
y
∂
∂
y
+
2
t
∂
∂
t
+
(
2
u
+
1
)
∂
∂
u
.


## 2.3. Group-Invariant Solutions

In this section we present the reduction forms of (1) by using symmetry group method. To do this, particular linear combinations of infinitesimals are considered and their corresponding invariants are determined. For example, for the symmetry operator X
=
X
1
+
X
2
+
X
3, we can compute the invariants by integrating the characteristic equations:
(19)
d
x
1
=
d
y
1
=
d
t
1
.
The corresponding invariants are q
=
x
-
y
,
p
=
t
-
y
,
ω
=
u. Now treating p
,
q as the new independent variables and ω as the new dependent variable, we obtain the ordinary differential equation:
(20)
(
2
ω
(
q
)
+
1
)
ω
′′
(
q
)
+
2
(
ω
′
(
q
)
)
2
+
ω
(
4
)
(
q
)
=
0
.
Integrating (20) twice with respect to q we obtain
(21)
ω
′′
(
q
)
+
ω
2
(
q
)
+
ω
(
q
)
=
c
0
q
+
c
1
,
c
0
,
c
1
∈
R
.
Taking c
0
=
0 and ω
′
(
q
)
=
ψ
(
ω
) we get
(22)
ψ
′
(
ω
)
ψ
(
ω
)
+
ω
2
+
ω
=
c
1
.
The solution of (22) is
(23)
ψ
(
ω
)
=
±
1
3
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
,
where c
2 is a constant. Notice that ω
′
(
q
)
=
ψ
(
ω
); we have
(24)
ω
′
(
q
)
=
±
1
3
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
,
so we get
(25)
q
=
±
3
∫
‍
1
-
6
ω
3
+
18
c
1
ω
-
9
ω
2
+
9
c
2
d
ω
.
Recall q
=
x
-
y
,
ω
=
u; we get the solution of (1):
(26)
x
-
y
=
±
3
∫
‍
1
-
6
u
3
+
18
c
1
u
-
9
u
2
+
9
c
2
d
u
.


## 3. Conservation Law

Here we briefly present the pertinent results of the two variational methods we utilize below.


### 3.1. Variational Method for a System and Its Adjoint

The system of adjoint equations to the system ofkth-order differential equation (3) is defined by [23]
(27)
E
α
*
(
x
,
u
,
v
,
…
,
u
(
k
)
,
v
(
k
)
)
=
0
,
α
=
1
,
…
,
m
,
where
(28)
E
α
*
(
x
,
u
,
v
,
…
,
u
(
k
)
,
v
(
k
)
)
=
δ
(
ν
β
E
β
)
δ
u
α
,
α
=
1
,
…
,
m
,
ν
=
ν
(
x
)
,
and ν
=
(
ν
1
,
ν
2
,
…
,
ν
m
) are new dependent variables.We recall here the following results as given in Ibragimov [23].A system of (3) is said to be self-adjoint if the substitution of ν
=
u into the system of adjoint equation (27) yields the same system (3).Assume the system of (3) admits the symmetry generator:
(29)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
,and then the system of adjoint equation (27) admits the operator:
(30)
Y
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
+
η
*
α
∂
∂
ν
α
,
η
*
α
=
-
[
λ
β
α
ν
β
+
ν
α
D
i
(
ξ
i
)
]
,
where the operator (30) is an extension of (29) to the variable ν
α and the λ
β
α are obtainable from
(31)
X
(
E
α
)
=
λ
α
β
E
β
.Theorem 1 (see [23]).
Every Lie point, Lie-Bäcklund, and nonlocal symmetry (29) admitted by the system of (3) gives rise to a conservation law for the system consisting of (3) and the adjoint equation (27), where the components T
i of the conserved vector T
=
(
T
1
,
…
,
T
n
) are determined by
(32)
T
i
=
ξ
i
L
+
W
α
δ
L
δ
u
i
α
+
∑
s
≥
1
‍
D
i
1
⋯
D
i
s
(
W
α
)
δ
L
δ
u
i
i
1
i
2
⋯
i
s
α
,
i
=
1
,
…
,
n
,
with Lagrangian given by
(33)
L
=
ν
α
E
α
(
x
,
u
,
…
,
u
(
k
)
)
.


### 3.2. Partial Noether Approach for a System of PDEs

The following results are due to Kara and Mahomed [26] that are based on the partial Lagrangian approach to construct conservation laws for a system of PDEs.Suppose that the system of (3) is written as
(34)
E
α
=
E
α
0
+
E
α
1
=
0
,
α
=
1
,
…
,
m
.
If there exist a function L
=
L
(
x
,
u
,
u
(
1
)
,
…
,
u
(
l
)
)
∈
A
,
l
≤
k and nonzero functions f
α
β
∈
A such that (34) can be written as δ
L
/
δ
u
α
=
f
α
β
E
β
1, provided E
β
1
≠
0, then L is known as a partial Lagrangian of (34); otherwise it is the standard Lagrangian. The differential equations of the form
(35)
δ
L
δ
u
α
=
f
α
β
E
β
1
are called a system of partial Euler-Lagrange equations.The operatorX in (10) is a partial Noether operator corresponding to a partial Lagrangian L
∈
A of the system (35) if it can be determined from
(36)
X
(
L
)
+
L
D
i
(
ξ
i
)
=
W
α
δ
L
δ
u
α
+
D
i
(
B
i
)
,
for some vector B
=
(
B
1
,
B
2
,
…
,
B
n
)
,
B
i
∈
A. Here W
=
(
W
1
,
W
2
,
…
,
W
m
)
,
W
α
∈
A is the characteristic of X.Theorem 2 (see [26]).
If the operatorX as in (10) is a partial Noether operator of a partial Lagrangian L corresponding to a partial Eule-Lagrange system of the form (35), the components T
i of the conserved vector T of (3) or (35) can be constructed by the following formula:
(37)
T
i
=
B
i
-
ξ
i
L
-
W
α
δ
L
δ
u
i
α
-
∑
s
≥
1
‍
D
i
1
⋯
D
i
s
(
W
α
)
δ
L
δ
u
i
i
1
i
2
·
·
·
i
s
α
,
i
=
1
,
…
,
n
,
where the characteristic W
=
(
W
1
,
W
2
,
…
,
W
m
)
,
W
α
∈
A of X is also the characteristic of the conservation law D
i
T
i
=
0 of (3).


### 3.3. Conservation Law

We now construct conservation law of (1) using the two approaches outlined above.Application of the New Conservation Theorem. The (2 + 1)-dimensional Boussinesq equation is given by
(38)
E
α
=
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
=
0
.
We recall that (38) admits the following five Lie point symmetry generators:
(39)
X
1
=
∂
∂
x
,
X
2
=
∂
∂
y
,
X
3
=
∂
∂
t
,
X
4
=
y
∂
∂
t
+
t
∂
∂
y
,
X
5
=
x
∂
∂
x
+
2
y
∂
∂
y
+
2
t
∂
∂
t
+
(
2
u
+
1
)
∂
∂
u
.
The adjoint equation of (38), by invoking (27), is
(40)
E
α
*
=
δ
[
v
(
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
)
]
δ
u
=
0
,
where v
=
v
(
x
,
y
,
t
) is a new dependent variable and (40) gives
(41)
v
t
t
-
2
u
x
x
v
-
4
u
x
v
x
-
(
2
u
+
1
)
v
x
x
-
v
y
y
-
v
x
x
x
x
=
0
.
By using Theorem 1, we obtain the following Lagrangian for the system of (38) and (40):
(42)
L
=
v
(
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
)
.
(1)
We first consider the Lie point symmetry generatorX
1
=
∂
/
∂
x, and we have W
=
-
u
x. Hence using (32), we obtain the following components of the conserved vector T:
(43)
T
x
1
=
(
2
u
+
2
v
-
1
)
u
x
v
x
-
u
x
x
v
x
x
-
u
x
v
x
x
x
,
T
t
1
=
-
u
x
t
v
+
u
x
v
t
,
T
y
1
=
u
x
y
v
-
u
x
v
y
.
(2)
The Lie point symmetry generatorX
2
=
∂
/
∂
y has the Lie characteristic function W
=
-
u
y. Thus by using (32), the components of the conserved vector are given by
(44)
T
x
2
=
[
2
u
x
u
y
+
(
2
u
+
1
)
u
x
y
+
u
y
x
x
x
]
v
-
(
2
u
u
y
+
u
y
)
v
x
-
u
x
y
v
x
x
-
u
y
v
x
x
x
,
T
t
2
=
-
u
y
t
v
+
u
y
v
t
,
T
y
2
=
[
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
]
v
-
u
y
v
y
-
v
v
x
x
x
x
.
(3)
Now the Lie point symmetry generatorX
3
=
∂
/
∂
t has the Lie characteristic function W
=
-
u
t. Hence we can obtain the conserved vector whose components are
(45)
T
x
3
=
[
2
u
x
u
t
+
(
2
u
+
1
)
u
x
t
+
u
x
x
x
t
]
v
-
(
2
u
+
1
)
u
t
v
x
-
u
x
t
v
x
x
-
u
t
v
x
x
x
,
T
t
3
=
-
[
2
u
x
2
+
(
2
u
+
1
)
u
x
x
+
u
y
y
]
v
+
u
t
v
t
-
v
v
x
x
x
x
,
T
y
3
=
-
u
y
t
v
-
u
t
v
y
.
(4)
The Lie point symmetry generatorX
4
=
y
(
∂
/
∂
t
)
+
t
(
∂
/
∂
y
) has the Lie characteristic function W
=
-
t
u
y
-
y
u
t. Hence using (32), the components of the conserved vector are given by
(46)
T
x
4
=
[
2
u
x
(
t
u
y
+
y
u
t
)
+
(
2
u
+
1
)
(
t
u
x
y
+
y
u
x
t
)
+
t
u
x
x
x
y
+
y
u
x
x
x
t
]
v
-
(
2
u
+
1
)
(
t
u
y
+
y
u
t
)
v
x
-
(
t
u
x
y
+
y
u
x
t
)
v
x
x
-
(
t
u
y
+
y
u
t
)
v
x
x
x
,
T
t
4
=
-
[
u
y
+
2
y
u
x
2
+
(
2
u
+
1
)
y
u
x
x
+
y
u
y
y
+
t
u
y
t
+
y
u
x
x
x
x
]
v
+
(
t
u
y
+
y
u
t
)
v
t
,
T
y
4
=
[
u
t
-
2
t
u
x
2
+
t
u
t
t
-
(
2
u
+
1
)
t
u
x
x
+
y
u
y
t
-
t
u
x
x
x
x
]
v
-
(
t
u
y
+
y
u
t
)
v
y
.
(5)
Finally, we consider the Lie point symmetry generatorX
5
=
x
(
∂
/
∂
x
)   +  2
y
(
∂
/
∂
y
)
+
2
t
(
∂
/
∂
t
)
+
(
2
u
+
1
)
(
∂
/
∂
u
) and the Lie characteristic function W
=
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t. Hence we can obtain the conserved vector whose components are
(47)
T
x
5
=
[
x
u
t
t
-
u
y
y
+
2
u
x
(
2
y
u
y
+
2
t
u
t
-
2
u
-
1
)
+
(
4
u
+
2
)
(
y
u
x
y
+
t
u
x
t
)
-
u
x
x
x
-
2
y
u
x
x
x
y
-
2
t
u
x
x
x
t
]
v
+
(
2
u
+
1
)
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
v
x
-
(
x
u
x
x
+
2
y
u
x
y
+
2
t
u
x
t
)
v
x
x
+
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
v
x
x
x
,
T
t
5
=
v
[
-
2
(
2
u
+
1
)
t
u
x
x
-
4
t
u
x
2
-
2
t
u
y
y
-
2
t
u
x
x
x
x
-
x
u
x
t
-
2
y
u
y
t
]
-
v
t
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
,
T
y
5
=
v
[
2
y
u
t
t
-
2
y
(
2
u
+
1
)
u
x
x
-
4
y
u
x
2
-
2
y
u
x
x
x
x
+
x
u
x
y
+
2
t
u
y
t
]
+
v
y
(
2
u
+
1
-
x
u
x
-
2
y
u
y
-
2
t
u
t
)
.Remark 3.
The conserved vector contains the arbitrary solutionv of the adjoint equation (40) and hence one can obtain an infinite number of conservation laws.Application of the Partial Lagrange Method. Consider the partial Lagrange approach given by Kara and Mahomed [26]. A partial Lagrange for (1) is
(48)
L
=
1
2
u
t
2
-
(
u
+
1
2
)
u
x
2
-
1
2
u
y
2
-
1
2
u
x
x
2
.
The Euler-Lagrange-type equation is
(49)
δ
L
δ
u
=
-
u
t
t
+
(
2
u
+
1
)
u
x
x
+
u
x
2
+
u
y
y
+
u
x
x
x
x
.
So (1) can be written as
(50)
δ
L
δ
u
=
-
u
x
2
,
to determine the Noether-type operators. If we substitute (48) and (50) into the partial Noether-type operators determining (36), we obtain
(51)
ξ
=
ϕ
=
τ
=
η
u
=
0
,
η
=
C
(
x
,
y
,
t
)
,
B
1
=
-
C
t
u
+
α
(
x
,
y
,
t
)
,
B
2
=
-
(
u
2
+
u
)
C
x
+
β
(
x
,
y
,
t
)
,
B
3
=
-
C
y
u
+
γ
(
x
,
y
,
t
)
,
subject to the condition
(52)
B
t
1
+
B
x
2
+
B
y
3
=
0
.
The conserved components, using (37), are
(53)
T
1
=
C
t
u
-
C
u
t
+
α
,
T
2
=
-
(
u
2
+
u
)
C
x
+
C
(
2
u
+
1
)
u
x
+
C
x
u
x
x
+
β
,
T
3
=
-
C
y
u
+
C
u
y
+
γ
.


## 3.1. Variational Method for a System and Its Adjoint

The system of adjoint equations to the system ofkth-order differential equation (3) is defined by [23]
(27)
E
α
*
(
x
,
u
,
v
,
…
,
u
(
k
)
,
v
(
k
)
)
=
0
,
α
=
1
,
…
,
m
,
where
(28)
E
α
*
(
x
,
u
,
v
,
…
,
u
(
k
)
,
v
(
k
)
)
=
δ
(
ν
β
E
β
)
δ
u
α
,
α
=
1
,
…
,
m
,
ν
=
ν
(
x
)
,
and ν
=
(
ν
1
,
ν
2
,
…
,
ν
m
) are new dependent variables.We recall here the following results as given in Ibragimov [23].A system of (3) is said to be self-adjoint if the substitution of ν
=
u into the system of adjoint equation (27) yields the same system (3).Assume the system of (3) admits the symmetry generator:
(29)
X
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
,and then the system of adjoint equation (27) admits the operator:
(30)
Y
=
ξ
i
∂
∂
x
i
+
η
α
∂
∂
u
α
+
η
*
α
∂
∂
ν
α
,
η
*
α
=
-
[
λ
β
α
ν
β
+
ν
α
D
i
(
ξ
i
)
]
,
where the operator (30) is an extension of (29) to the variable ν
α and the λ
β
α are obtainable from
(31)
X
(
E
α
)
=
λ
α
β
E
β
.Theorem 1 (see [23]).
Every Lie point, Lie-Bäcklund, and nonlocal symmetry (29) admitted by the system of (3) gives rise to a conservation law for the system consisting of (3) and the adjoint equation (27), where the components T
i of the conserved vector T
=
(
T
1
,
…
,
T
n
) are determined by
(32)
T
i
=
ξ
i
L
+
W
α
δ
L
δ
u
i
α
+
∑
s
≥
1
‍
D
i
1
⋯
D
i
s
(
W
α
)
δ
L
δ
u
i
i
1
i
2
⋯
i
s
α
,
i
=
1
,
…
,
n
,
with Lagrangian given by
(33)
L
=
ν
α
E
α
(
x
,
u
,
…
,
u
(
k
)
)
.


## 3.2. Partial Noether Approach for a System of PDEs

The following results are due to Kara and Mahomed [26] that are based on the partial Lagrangian approach to construct conservation laws for a system of PDEs.Suppose that the system of (3) is written as
(34)
E
α
=
E
α
0
+
E
α
1
=
0
,
α
=
1
,
…
,
m
.
If there exist a function L
=
L
(
x
,
u
,
u
(
1
)
,
…
,
u
(
l
)
)
∈
A
,
l
≤
k and nonzero functions f
α
β
∈
A such that (34) can be written as δ
L
/
δ
u
α
=
f
α
β
E
β
1, provided E
β
1
≠
0, then L is known as a partial Lagrangian of (34); otherwise it is the standard Lagrangian. The differential equations of the form
(35)
δ
L
δ
u
α
=
f
α
β
E
β
1
are called a system of partial Euler-Lagrange equations.The operatorX in (10) is a partial Noether operator corresponding to a partial Lagrangian L
∈
A of the system (35) if it can be determined from
(36)
X
(
L
)
+
L
D
i
(
ξ
i
)
=
W
α
δ
L
δ
u
α
+
D
i
(
B
i
)
,
for some vector B
=
(
B
1
,
B
2
,
…
,
B
n
)
,
B
i
∈
A. Here W
=
(
W
1
,
W
2
,
…
,
W
m
)
,
W
α
∈
A is the characteristic of X.Theorem 2 (see [26]).
If the operatorX as in (10) is a partial Noether operator of a partial Lagrangian L corresponding to a partial Eule-Lagrange system of the form (35), the components T
i of the conserved vector T of (3) or (35) can be constructed by the following formula:
(37)
T
i
=
B
i
-
ξ
i
L
-
W
α
δ
L
δ
u
i
α
-
∑
s
≥
1
‍
D
i
1
⋯
D
i
s
(
W
α
)
δ
L
δ
u
i
i
1
i
2
·
·
·
i
s
α
,
i
=
1
,
…
,
n
,
where the characteristic W
=
(
W
1
,
W
2
,
…
,
W
m
)
,
W
α
∈
A of X is also the characteristic of the conservation law D
i
T
i
=
0 of (3).


## 3.3. Conservation Law

We now construct conservation law of (1) using the two approaches outlined above.Application of the New Conservation Theorem. The (2 + 1)-dimensional Boussinesq equation is given by
(38)
E
α
=
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
=
0
.
We recall that (38) admits the following five Lie point symmetry generators:
(39)
X
1
=
∂
∂
x
,
X
2
=
∂
∂
y
,
X
3
=
∂
∂
t
,
X
4
=
y
∂
∂
t
+
t
∂
∂
y
,
X
5
=
x
∂
∂
x
+
2
y
∂
∂
y
+
2
t
∂
∂
t
+
(
2
u
+
1
)
∂
∂
u
.
The adjoint equation of (38), by invoking (27), is
(40)
E
α
*
=
δ
[
v
(
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
)
]
δ
u
=
0
,
where v
=
v
(
x
,
y
,
t
) is a new dependent variable and (40) gives
(41)
v
t
t
-
2
u
x
x
v
-
4
u
x
v
x
-
(
2
u
+
1
)
v
x
x
-
v
y
y
-
v
x
x
x
x
=
0
.
By using Theorem 1, we obtain the following Lagrangian for the system of (38) and (40):
(42)
L
=
v
(
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
-
u
y
y
-
u
x
x
x
x
)
.
(1)
We first consider the Lie point symmetry generatorX
1
=
∂
/
∂
x, and we have W
=
-
u
x. Hence using (32), we obtain the following components of the conserved vector T:
(43)
T
x
1
=
(
2
u
+
2
v
-
1
)
u
x
v
x
-
u
x
x
v
x
x
-
u
x
v
x
x
x
,
T
t
1
=
-
u
x
t
v
+
u
x
v
t
,
T
y
1
=
u
x
y
v
-
u
x
v
y
.
(2)
The Lie point symmetry generatorX
2
=
∂
/
∂
y has the Lie characteristic function W
=
-
u
y. Thus by using (32), the components of the conserved vector are given by
(44)
T
x
2
=
[
2
u
x
u
y
+
(
2
u
+
1
)
u
x
y
+
u
y
x
x
x
]
v
-
(
2
u
u
y
+
u
y
)
v
x
-
u
x
y
v
x
x
-
u
y
v
x
x
x
,
T
t
2
=
-
u
y
t
v
+
u
y
v
t
,
T
y
2
=
[
u
t
t
-
(
2
u
+
1
)
u
x
x
-
2
u
x
2
]
v
-
u
y
v
y
-
v
v
x
x
x
x
.
(3)
Now the Lie point symmetry generatorX
3
=
∂
/
∂
t has the Lie characteristic function W
=
-
u
t. Hence we can obtain the conserved vector whose components are
(45)
T
x
3
=
[
2
u
x
u
t
+
(
2
u
+
1
)
u
x
t
+
u
x
x
x
t
]
v
-
(
2
u
+
1
)
u
t
v
x
-
u
x
t
v
x
x
-
u
t
v
x
x
x
,
T
t
3
=
-
[
2
u
x
2
+
(
2
u
+
1
)
u
x
x
+
u
y
y
]
v
+
u
t
v
t
-
v
v
x
x
x
x
,
T
y
3
=
-
u
y
t
v
-
u
t
v
y
.
(4)
The Lie point symmetry generatorX
4
=
y
(
∂
/
∂
t
)
+
t
(
∂
/
∂
y
) has the Lie characteristic function W
=
-
t
u
y
-
y
u
t. Hence using (32), the components of the conserved vector are given by
(46)
T
x
4
=
[
2
u
x
(
t
u
y
+
y
u
t
)
+
(
2
u
+
1
)
(
t
u
x
y
+
y
u
x
t
)
+
t
u
x
x
x
y
+
y
u
x
x
x
t
]
v
-
(
2
u
+
1
)
(
t
u
y
+
y
u
t
)
v
x
-
(
t
u
x
y
+
y
u
x
t
)
v
x
x
-
(
t
u
y
+
y
u
t
)
v
x
x
x
,
T
t
4
=
-
[
u
y
+
2
y
u
x
2
+
(
2
u
+
1
)
y
u
x
x
+
y
u
y
y
+
t
u
y
t
+
y
u
x
x
x
x
]
v
+
(
t
u
y
+
y
u
t
)
v
t
,
T
y
4
=
[
u
t
-
2
t
u
x
2
+
t
u
t
t
-
(
2
u
+
1
)
t
u
x
x
+
y
u
y
t
-
t
u
x
x
x
x
]
v
-
(
t
u
y
+
y
u
t
)
v
y
.
(5)
Finally, we consider the Lie point symmetry generatorX
5
=
x
(
∂
/
∂
x
)   +  2
y
(
∂
/
∂
y
)
+
2
t
(
∂
/
∂
t
)
+
(
2
u
+
1
)
(
∂
/
∂
u
) and the Lie characteristic function W
=
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t. Hence we can obtain the conserved vector whose components are
(47)
T
x
5
=
[
x
u
t
t
-
u
y
y
+
2
u
x
(
2
y
u
y
+
2
t
u
t
-
2
u
-
1
)
+
(
4
u
+
2
)
(
y
u
x
y
+
t
u
x
t
)
-
u
x
x
x
-
2
y
u
x
x
x
y
-
2
t
u
x
x
x
t
]
v
+
(
2
u
+
1
)
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
v
x
-
(
x
u
x
x
+
2
y
u
x
y
+
2
t
u
x
t
)
v
x
x
+
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
v
x
x
x
,
T
t
5
=
v
[
-
2
(
2
u
+
1
)
t
u
x
x
-
4
t
u
x
2
-
2
t
u
y
y
-
2
t
u
x
x
x
x
-
x
u
x
t
-
2
y
u
y
t
]
-
v
t
[
(
2
u
+
1
)
-
x
u
x
-
2
y
u
y
-
2
t
u
t
]
,
T
y
5
=
v
[
2
y
u
t
t
-
2
y
(
2
u
+
1
)
u
x
x
-
4
y
u
x
2
-
2
y
u
x
x
x
x
+
x
u
x
y
+
2
t
u
y
t
]
+
v
y
(
2
u
+
1
-
x
u
x
-
2
y
u
y
-
2
t
u
t
)
.Remark 3.
The conserved vector contains the arbitrary solutionv of the adjoint equation (40) and hence one can obtain an infinite number of conservation laws.Application of the Partial Lagrange Method. Consider the partial Lagrange approach given by Kara and Mahomed [26]. A partial Lagrange for (1) is
(48)
L
=
1
2
u
t
2
-
(
u
+
1
2
)
u
x
2
-
1
2
u
y
2
-
1
2
u
x
x
2
.
The Euler-Lagrange-type equation is
(49)
δ
L
δ
u
=
-
u
t
t
+
(
2
u
+
1
)
u
x
x
+
u
x
2
+
u
y
y
+
u
x
x
x
x
.
So (1) can be written as
(50)
δ
L
δ
u
=
-
u
x
2
,
to determine the Noether-type operators. If we substitute (48) and (50) into the partial Noether-type operators determining (36), we obtain
(51)
ξ
=
ϕ
=
τ
=
η
u
=
0
,
η
=
C
(
x
,
y
,
t
)
,
B
1
=
-
C
t
u
+
α
(
x
,
y
,
t
)
,
B
2
=
-
(
u
2
+
u
)
C
x
+
β
(
x
,
y
,
t
)
,
B
3
=
-
C
y
u
+
γ
(
x
,
y
,
t
)
,
subject to the condition
(52)
B
t
1
+
B
x
2
+
B
y
3
=
0
.
The conserved components, using (37), are
(53)
T
1
=
C
t
u
-
C
u
t
+
α
,
T
2
=
-
(
u
2
+
u
)
C
x
+
C
(
2
u
+
1
)
u
x
+
C
x
u
x
x
+
β
,
T
3
=
-
C
y
u
+
C
u
y
+
γ
.


## 4. Conclusions

In this paper, we investigate Lie point symmetries, similarity reduction, invariant solutions, and conservation laws of the (2 + 1)-dimensional Boussinesq equation. One of the most important applications of theory of Lie groups is to obtain the conservation laws of differential equations. It is well known that the famous Noether theorem establishes a connection between symmetries and conservation laws of differential equations provided that the equations are Euler-Lagrange equations. However, (1) does not admit Lagrangians. We obtain conservation laws of (1) by utilizing the new conservation theorem and the partial Lagrange approach.



---

*Source: 840405-2014-08-07.xml*