1. Introduction
The Benjamin–Bona–Mahony equation
was, for the first time, studied by Benjamin
et al. [
1]. It is also known as the regularized long-wave equation and is applicable to shallow water waves and to the study of drift waves in plasma or the Rossby waves in rotating fluids. See, for example, [
2] and references therein. It should be noted, however, that the BBM equation was also derived by Peregrine [
3].
In the last two decades, various versions of the Benjamin–Bona–Mahony equation have been investigated in the literature. A general form of the Benjamin–Bona–Mahony equation is [
2]
where
α,
β and
δ are constants with the nonlinear and dispersion coefficients
β ≠ 0 and
δ < 0, respectively. For different values of the constants
α,
β and
δ,
Equation (2), results in various types of nonlinear equations which are very useful in the study of various physical phenomena. See [
2] and references therein. However, in certain cases, the physical situations of the problem dictate us to consider nonlinear equations with variable coefficients [
2,
4–
8]. Recently, in [
2] the variable coefficients version of the Benjamin–Bona–Mahony equation
was investigated and some exact solutions were obtained using the classical Lie group method [
9].
In this study, we consider the
Equation (3), but with
δ being an arbitrary constant, namely
where
α(
t) and
β(
t) are arbitrary functions of
t. The objective of the study is to classify the Noether operators and to construct conservation laws for the
Equation (4).
Conservation laws are mathematical expressions of the physical laws, such as conservation of energy, mass, momentum and so on. They play a very important role in the solution and reduction of partial differential equations. Conservation laws have been widely used to investigate the existence, uniqueness and stability of solutions of nonlinear partial differential equations. This can be seen in the references [
10–
12]. They have also been employed in the development and use of numerical methods (see for example, [
13,
14]). Recently, conserved vectors associated with Lie point symmetries have been used to find exact solutions of some partial differential equations [
15] and systems of ordinary equations [
16]. Noether theorem [
17] gives us an elegant way to construct conservation laws provided a Lagrangian is known for an Euler–Lagrange equation. Thus, the knowledge of a Lagrangian is essential in this case.
It is worth mentioning that no Lagrangian exists for
Equation (4) and as a result one can not invoke Noether theorem. However, an interesting approach is employed to construct conservation laws for
Equation (4). It should be noted that the present approach, which we use here, fails to construct the conservation laws for
Equation (3),
i.e., when
δ is a function of
t.
The paper is organized as follows. In Section 2 we briefly give the fundamental notations and relations concerning the Noether symmetry approach. Section 3 obtains the Noether operators and the corresponding conservation laws for the
Equation (4). In Section 4, a conservation law is used to obtain exact solution of a special case of the Benjamin–Bona–Mahony equation. Concluding remarks are presented in Section 5.
Consider the Benjamin–Bona–Mahony
Equation (4). Here, we note that
Equation (4) does not admit a Lagrangian. Nevertheless, we can transform
Equation (4) into a variational form by setting
u =
zx. Thus,
Equation (4) becomes a fourth-order equation, namely
It can easily be verified that a second-order Lagrangian of the
Equation (15) is
The substitution of
L from
Equation (16) into
Equation (12) and splitting with respect to the derivatives of
z, yields an overdetermined system of linear PDEs
After some lengthy calculations, the solution of the above system yields
Case 1. α(t) and β(t) arbitrary but not of the form contained in Cases 2–4.
In this case, we obtain two Noether point symmetries. These are given below together with their corresponding gauge functions:
Invoking Theorem 1, and reverting back to the original variables, the two nontrivial conserved vectors associated with these two Noether point symmetries are, respectively,
and
Here it can be seen that the above conserved vectors do not satisfy the divergence condition, viz.,
DiTi|(4) = 0, as some excessive terms emerge that require some further analysis. By making a slight adjustment to these terms, it can be shown that these terms can be absorbed into the divergence condition. For,
hence
We now redefine the conserved vectors in the parenthesis as:
Thus, the modified conserved vectors
and
satisfy the divergence condition. We observe that the conserved vectors
Equations (38) and
(39) are local conserved vectors. Likewise, we can then construct the second pair of the conserved quantities associated with
and
as:
Note that the conserved quantities
Equations (40) and
(41) are nonlocal conserved vectors, and since
d(
t) is an arbitrary function of
t, one obtains infinitely many nonlocal conserved vectors of
Equation (4). A special case of
Equations (40) and
(41), when
d(
t) = 1, is
Case 2. α(t) = γ, β(t) = λ, where γ and λ are non-zero constants.
This case provides us with three Noether symmetry generators, namely,
X1, X2 given by the operators
Equations (30) and
(31) and
X3 given by
Again, the above conserved flow fails to satisfy the divergence criterion. Thus, by inheriting the same procedure as in Case 1, the nontrivial conserved flows associated with
X1, X2 and
X3 are
and
respectively. The conserved vector
Equations (47) and
(48) is a local conserved vector whereas the conserved vectors
Equations (49) and
(52) are nonlocal conserved quantities. It should be noted that one can use
Equations (49) and
(50) to construct infinitely many nonlocal conserved vectors.
Case 3. α(t) = γeσt, β(t) = λeσt, where γ, σ and λ are nonzero constants.
Here we get three Noether point symmetries operators, viz.,
X1, X2 given by the operators
Equations (30) and
(31) and
X3 given by
The application of Theorem 1, with
X3, gives
and as before, the modified conserved vectors are given by
Thus, the nontrivial conserved flows associated with
X1,
X2 and
X3 are, respectively,
Case 4. α(t) = γ/a(t), β(t) = λ/a(t), where γ, λ are constants, with γ, λ ≠ 0 and a(t) an arbitrary function of t.
In this case, we obtain three Noether point symmetry generators, viz.,
X1, X2 given by
Equations (30) and
(31) and
X3 given by
Following the above procedure, the nontrivial local and nonlocal conserved quantities corresponding to
X1, X2 and
X3, in this case, are
respectively.
Remark 1. Remark: It should be noted that since the Lagrangian Equation (16) is invariant under the spatial translation symmetry, this will give rise to the linear momentum conservation laws. Moreover, if α(
t) =
γ,
β(
t) =
λ, where γ and λ are non-zero constants, then the corresponding Lagrangian Equation (16) is also invariant under the time translation symmetry and thus the linear momentum and energy are both conserved. 4. Exact Solution of Equation (4) for a Special Case Using Conservation Laws
First we recall a definition and theorem from [
15] that we utilize in this Section.
Definition 2. Suppose that X is a symmetry of Equation (10) and T a conserved vector of Equation (10). Then if X and T satisfythen X is said to be associated with T. Define a nonlocal variable
v by
Tt =
vx, Tx = −
vt. Then using the similarity variables
r,
s,
w with the generator
, we have
Tr =
vs,
Ts = −
vr, and the conservation law is then rewritten as
where
Theorem 2. An n-th order partial differential equation with two independent variables, which admits a symmetry X that is associated with a conserved vector T, is reduced to an ordinary differential equation of order n − 1
; namely, Tr =
k, where Tr is defined by Equation (72) for solutions invariant under X. We now use the above definition and theorem to obtain an exact solution of one special case of
Equation (4) by making use of its conservation law. We consider Case 2 of Section 3.
Recall that the
Equation (4) with
α(
t) =
γ,
β(
t) =
λ, where
γ and
λ are non-zero constants, admits (among others)
and possesses the conservation law with conserved vector
that is associated with both
X1 and
X2. We now define the combination of
X1 and
X2 by
X =
ρX1+
X2. Thus, the canonical coordinates of
X are given by
where
u =
u(
r) is the invariant solution under
X if
u satisfies
Equation (4) with
α(
t) =
γ,
β(
t) =
λ. Employing
Equation (72), the
r-component of
T given in
Equation (75) is
By letting
ur =
p(
u), we have ... Then
Equation (78) reduces to the first order ordinary differential equation
The integration of
Equation (79) leads to the four parameter family of solutions
of
Equation (4) invariant under
X =
ρX1 +
X2.