Next Article in Journal
Analysis and Optimal Control of φ-Hilfer Fractional Semilinear Equations Involving Nonlocal Impulsive Conditions
Previous Article in Journal
Deep Learning Framework to Detect Ischemic Stroke Lesion in Brain MRI Slices of Flair/DW/T1 Modalities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetry Analysis, Exact Solutions and Conservation Laws of a Benjamin–Bona–Mahony–Burgers Equation in 2+1-Dimensions

by
María S. Bruzón
,
Tamara M. Garrido-Letrán
and
Rafael de la Rosa
*
Departamento de Matemáticas, Facultad de Ciencias, Universidad de Cádiz, 11510 Puerto Real, Cádiz, Spain
*
Author to whom correspondence should be addressed.
M.S. Bruzón sadly passed away prior to the revision of this paper. This is one of her last works.
Symmetry 2021, 13(11), 2083; https://doi.org/10.3390/sym13112083
Submission received: 2 August 2021 / Revised: 2 October 2021 / Accepted: 26 October 2021 / Published: 3 November 2021
(This article belongs to the Section Physics)

Abstract

:
The Benjamin–Bona–Mahony equation describes the unidirectional propagation of small-amplitude long waves on the surface of water in a channel. In this paper, we consider a family of generalized Benjamin–Bona–Mahony–Burgers equations depending on three arbitrary constants and an arbitrary function G ( u ) . We study this family from the standpoint of the theory of symmetry reductions of partial differential equations. Firstly, we obtain the Lie point symmetries admitted by the considered family. Moreover, taking into account the admitted point symmetries, we perform symmetry reductions. In particular, for G ( u ) 0 , we construct an optimal system of one-dimensional subalgebras for each maximal Lie algebra and deduce the corresponding (1+1)-dimensional nonlinear third-order partial differential equations. Then, we apply Kudryashov’s method to look for exact solutions of the nonlinear differential equation. We also determine line soliton solutions of the family of equations in a particular case. Lastly, through the multipliers method, we have constructed low-order conservation laws admitted by the family of equations.

1. Introduction

The Benjamin–Bona–Mahony equation (BBM), or regularised long wave (RLW) equation,
u t + u x + u u x u x x t = 0 ,
was first proposed by Benjamin et al. [1] as an alternative mathematical model to the Korteweg–de Vries for modelling long wave motions in nonlinear dispersive systems. The authors stressed that both models are applicable at the same level of approximation, however, from a computational mathematics point of view, the BBM equation has some advantages over the KdV equation. Generalized forms of the BBM equation have been widely studied. In this paper, we consider a generalized family of Benjamin–Bona–Mahony–Burgers (gBBMB) equations given by
u t γ ( u t x x + u t y y ) α ( u x x + u y y ) + β ( u x + u y ) = ( F u u x + F u u y ) + f ,
where u = u ( t , x , y ) is an analytic function of the time coordinate t and the spatial coordinates x and y, γ 0 , α , β and f are arbitrary constants, whereas F ( u ) is a nonlinear function. The gBBMB Equation (2) was considered in [2]. Here, the authors considered a new class of polynomial functions equipped with a parameter to approximate the solutions of the gBBMB Equation (2).
Regarding the physical interpretation, the original BBM equation was proposed as a model of long waves in channel flows where nonlinear dispersion is incorporated. Its solutions approximate solutions of the two-dimensional Euler equations. Hence, for the two-dimensional model considered here (2), the variables t , x , y keep the same physical interpretations as in Euler equations [3,4]. For two-dimensional water waves in a channel with a flat bottom, the independent variable t is interpreted as the elapsed time and x , y determine the position, the horizontal and vertical coordinates along the channel, respectively.
Currently, nonlinear equations involving parameters and arbitrary functions, which arise in diverse fields as mathematical biology and physics, nonlinear dynamics and plasma physics, have attracted the attention of numerous researchers. Nevertheless, the analysis of nonlinear equations involving arbitrary functions is often difficult and laborious. In particular, for dealing with the determination of exact solutions, several direct methods have been elaborated, among them Kudryashov method [5,6,7,8], tanh-sech method [9,10], Painlevé analysis [11,12], Adomian decomposition method [13,14,15,16] and other special methods [17,18]. However, these methods just work for a limited kind of equations.
Lie symmetries of a partial differential equation (PDE) are transformations acting on the space of independent and dependent variables which transforms the PDE solution space into itself. The analysis of Lie symmetries is one of the most effective algorithms to analyse PDE equations including the construction of invariant solutions, construction of mappings between equivalent equations of the same family, finding invertible mappings of nonlinear PDEs to linear PDEs or finding conservation laws. Among these applications, we highlight the construction of invariant solutions, i.e., interesting special classes of solutions which are invariant under a Lie group of point transformations. In the case of a PDE, the invariance under a Lie group of point transformations allows one to obtain, constructively, similarity solutions (invariant solutions), which are invariant under a subgroup of the full Lie symmetry group admitted by the PDE. Similarity solutions arise as solutions of a reduced system of differential equations (DEs) with fewer number of independent variables. For further information on Lie symmetries and their applications one can refer, as example, to [19,20,21,22,23] and references therein.
In many natural processes, such as physical or chemical ones, conservation laws emerge. Considering an isolated physical system, these laws characterise physical properties that do not change over time. Regarding PDEs, not all conservation laws have a physical interpretation. However, when a PDE has a large number of conservation laws it is usually an indicator of its integrability. Conservation laws are also studied for their applicability in numerical methods for PDEs since they are useful to investigate the existence, uniqueness and stability of solutions. Moreover, exact solutions can be constructed by using conservation laws.
With respect to the concept of conservation laws, different results can be found in the the recent literature. In [24], Ibragimov proved a theorem to find conservation laws which do not require the existence of a classical Lagrangian. This theorem is based on the concept of adjoint equation for nonlinear DEs [25,26,27]. Nontrivial conservation laws have been determined by using Ibragimov’s conservation law theorem, see, for instance, Refs. [28,29,30]. Anco and Bluman presented a direct algorithmic method using adjoint-symmetries for finding all local conservation laws of a given DE system [31,32,33,34]. Moreover, in [35], this general method was further developed and reviewed in detail. In [36,37], symmetry properties of conservation laws of PDEs are analysed by using the multiplier method. In particular, in [36], it was proved that Ibragimov’s formula [24] for determining conservation laws is equivalent to a standard formula for the action of an infinitesimal symmetry on a conservation law multiplier. The same author recently showed [38] that Ibragimov’s method can lead to trivial conservation laws. Most importantly, the formula does not necessarily yield all non-trivial conservation laws unless the symmetry action on the set of these conservation laws is transitive, property which cannot be verified until all conservation laws have been determined. It is worth highlighting the symmetry multi-reduction method [39], which is a generalization of the double reduction method [40,41,42]. The method proposed in [39] gives an explicit algorithm for PDEs with n 2 independent variables admitting a symmetry algebra whose dimension is at least n 1 that allows us to find all symmetry-invariant conservation laws, which will reduce to first integrals for the ODE that describe the symmetry-invariant solutions of the PDE.
In [43], the following generalized BBMB equations with power functions were studied
u t + b u x + a ( u m ) x + ( u n ) t x x + c ( u k ) x x = 0 .
Lie symmetries were used to reduce the equation into ordinary differential equations (ODEs) in order to obtain new solutions. Furthermore, conservation laws were obtained for special values of the parameters. Here we similarly analyse Equation (2). In fact, considering n = k = 1 and changing u m to F ( u ) in Equation (3), the one-dimensional version of Equation (2) arises.
The aim of this paper is to analyse Equation (2) from the viewpoint of Lie symmetries, conservation laws and analytical solutions. First, in Section 2, we determine a complete classification of the point symmetries admitted by PDE (2) depending on the arbitrary parameters γ , α , β and f and the arbitrary function F ( u ) . In Section 3, we use the one-dimensional point symmetry groups admitted by Equation (2) to reduce the PDE into ODEs. Then, in Section 4, exact solutions are obtained by using Kudryashov method. In Section 5, by using the multipliers method, a complete classification of the low-order conservation laws of PDE (6) has been achieved. In Section 6, we obtain line soliton solutions u ( t , x , y ) = h ( x + μ y λ t ) , where μ and λ represent the direction and the speed of propagation of the line soliton. Taking into account the line soliton formulation into PDE (2), a nonlinear third-order ODE for h is obtained. This equation is reduced to a second-order ODE through the corresponding differential invariants. Moreover, the second-order ODE can be easily integrated that leads to its complete quadrature. Undoing the change of variables we obtain a first-order separable ODE which can be integrated to determine the analytical expression of the line soliton solutions. Finally, in Section 7, we present the conclusions.

2. Lie Point Symmetries

To determine the Lie point symmetries of PDE (2), we consider a one-parameter group of infinitesimal transformations in ( t , x , y , u ) given by
t ˜ = t + ε τ ( t , x , y , u ) + O ( ε 2 ) , x ˜ = x + ε ξ ( t , x , y , u ) + O ( ε 2 ) , y ˜ = y + ε ϕ ( t , x , y , u ) + O ( ε 2 ) , u ˜ = u + ε η ( t , x , y , u ) + O ( ε 2 ) ,
where ε is the group parameter. An admitted infinitesimal generator is of the form
X = τ ( t , x , y , u ) t + ξ ( t , x , y , u ) x + ϕ ( t , x , y , u ) y + η ( t , x , y , u ) u .
We recall that a group of transformations (4) is a Lie point symmetry of PDEEquation (2) if and only if the solution space of PDE Equation (2) is invariant under the action of the one-parameter group of point transformations (4). Lie’s fundamental theorems prove that Lie groups of transformations are entirely characterized by their infinitesimal generators. An infinitesimal formulation of Lie groups of transformations is crucial since it replaces nonlinear conditions of invariance of a given DE by equivalent, albeit far simpler, linear conditions which determine the infinitesimal generators of the group. Furthermore, for a given DE, the essential applications of Lie groups only need knowledge of the admitted infinitesimal generators. For further information one can refer, for example, to references [19,20,21,23]. For the sake of simplicity, it is helpful to introduce the function G ( u ) = β F u . Therefore, Equation (2) becomes
u t γ ( u t x x + u t y y ) α ( u x x + u y y ) + G ( u x + u y ) f = 0 .
Point symmetries are determined by applying the symmetry invariance criterion [19,20,21,23], which requires that
X ( 3 ) u t γ ( u t x x + u t y y ) α ( u x x + u y y ) + G ( u x + u y ) f = 0 ,
when Equation (6) holds. Here, X ( 3 ) represents the third prolongation of the infinitesimal generator (5), which is given by
X ( 3 ) = X + η i 1 ( 1 ) u i 1 + η i 1 i 2 ( 2 ) u i 1 i 2 + η i 1 i 2 i 3 ( 3 ) u i 1 i 2 i 3 ,
with coefficients
η i 1 ( 1 ) = D i 1 ( η ) u t D i 1 ( τ ) u x D i 1 ( ξ ) u y D i 1 ( ϕ ) , η i 1 i 2 ( 2 ) = D i 2 ( η i 1 ( 1 ) ) u i 1 t D i 2 ( τ ) u i 1 x D i 2 ( ξ ) u i 1 y D i 2 ( ϕ ) , η i 1 i 2 i 3 ( 3 ) = D i 3 ( η i 1 i 2 ( 2 ) ) u i 1 i 2 t D i 3 ( τ ) u i 1 i 2 x D i 3 ( ξ ) u i 1 i 2 y D i 3 ( ϕ ) ,
with D the total derivative operator, u i = u x i , i = 1 , 2 , 3 with x 1 = t , x 2 = x and x 3 = y , and i j = 1 , 2 , 3 for j = 1 , 2 , 3 . The invariance condition (7) splits with respect to the differential consequences of u, yielding a set of 64 determining equations. By simplifying this system, we obtain τ = τ ( t ) , ξ = ξ ( x , y ) , ϕ = ϕ ( x , y ) , η = η ( t , x , y , u ) , and the parameters α , γ and f along with the arbitrary function G ( u ) are related by the following conditions:
η u u = 0 , ξ y + ϕ x = 0 , ξ x ϕ y = 0 , γ η t u + α τ t = 0 , 2 ξ x γ η x x u + η y y u = 0 , ξ x x + ξ y y 2 η x u = 0 , ϕ x x + ϕ y y 2 η y u = 0 , 2 γ η t x u η G u G τ t + ξ x ξ y = 0 , 2 γ η t y u η G u G τ t + ϕ y ϕ x = 0 , γ η t x x + η t y y + α η x x + η y y G η x + η y + f τ t + 2 ξ x η u η t = 0 .
The determining system has been derived and solved with the aid of Maple commands “rifsimp” and “pdsolve”. Moreover, it should be noted that the gBBMB Equation (6) are preserved under the equivalence transformation given by
u ˜ u + k ,
with k constant. So we obtain the result:
Theorem 1.
The classification of point symmetries admitted by the gBBMB Equation (6) depending on the arbitrary constants α, γ and f, and the arbitrary function G ( u ) is given by the following cases:
(i) 
For arbitrary α, γ, f and G ( u ) the admitted point symmetries are generated by:
X 1 = t , t i m e - t r a n s l a t i o n . X 2 = x , s p a c e - t r a n s l a t i o n . X 3 = y , s p a c e - t r a n s l a t i o n .
(ii) 
Extra point symmetries of the gBBMB Equation (6) are admitted in the following cases:
(a) 
If α = 0 and G ( u ) = g 1 exp ( m u )
  • For f 0
    X 4 = exp ( f m t ) t + f u , d i l a t i o n .
  • For f = 0
    X 5 = m t t u , d i l a t i o n   a n d   s h i f t .
(b) 
If α = 0 and G ( u ) = g 1 u
X 6 = t t + u u , s c a l i n g .
(c) 
If α = f = 0 and G ( u ) = g 1 u n
X 7 = n t t u u , s c a l i n g .
(d) 
If G ( u ) = g 1
X ϕ = ϕ u , X ψ = u ψ u ,
where ϕ ( t , x , y ) and ψ ( t , x , y ) satisfies, respectively,
ϕ t γ ( ϕ t x x + ϕ t y y ) α ( ϕ x x + ϕ y y ) + g 1 ( ϕ x + ϕ y ) = 0 ,
and
ψ t γ ( ψ t x x + ψ t y y ) α ( ψ x x + ψ y y ) + g 1 ( ϕ x + ϕ y ) f = 0 .
In the above, g 1 0 , m 0 , n 0 are arbitrary constants.

3. Symmetry Reductions

In this section, we will restrict our attention to those cases where G ( u ) 0 . By using the point symmetries admitted by PDE (6), we can reduce PDE (6) to an equation with fewer number of independent variables. Each point symmetry of PDE (6) leads to a characteristic system
d t τ = d x ξ = d y ϕ = d u η .
Solving the characteristic equations, one obtains similarity variables z and r, and similarity solutions w ( z , r ) . The substitution of these variables into PDE (6) leads to third-order nonlinear PDEs for w ( z , r ) . In general, it is not always possible to determine all the possible group-invariant solutions, since the Lie symmetry group of a given PDE can contain an infinite number of Lie subgroups. The aim is to classify all the feasible group-invariant solutions into classes in a way that solutions belonging to the same class are equivalent, i.e., these solutions are related through an element of the Lie symmetry group; vice versa, solutions belonging to different classes are not equivalent and therefore there exists no element of the Lie symmetry group that maps one solution into the other. To address this problem, we will find an optimal system of one-dimensional subalgebras [23,44]. For that purpose, it is very useful to determine the most general symmetry Lie algebra that PDE (6) admits depending on the arbitrary function G ( u ) and the arbitrary parameters α , γ and f. The following theorem shows a basis of generators for each maximal Lie algebra admitted by PDE (6).
Theorem 2.
The maximal Lie algebras for PDE (6), with G ( u ) 0 , along with their non-zero commutator structure are given by:
1.
For arbitrary α, γ, f and G ( u )
A 1 = s p a n ( X 1 , X 2 , X 3 ) .
2.
If α = 0 , f 0 and G ( u ) = g 1 exp ( m u )
A 2 = s p a n ( X 1 , X 2 , X 3 , X 4 ) , X 1 , X 4 = f m X 4 .
3.
If α = f = 0 and G ( u ) = g 1 exp ( m u )
A 3 = s p a n ( X 1 , X 2 , X 3 , X 5 ) , X 1 , X 5 = m X 1 .
4.
If α = 0 and G ( u ) = g 1 u
A 4 = s p a n ( X 1 , X 2 , X 3 , X 6 ) , X 1 , X 6 = X 1 .
5.
If α = f = 0 and G ( u ) = g 1 u n
A 5 = s p a n ( X 1 , X 2 , X 3 , X 7 ) , X 1 , X 7 = n X 1 .
Theorem 3.
For the gBBMB Equation (6), with G ( u ) 0 , an optimal system of one-dimensional subalgebras for each maximal Lie algebra is given by:
1.
For arbitrary α, γ, f and G ( u )
X 1 + ν X 2 + μ X 3 , X 2 , X 3 .
2.
If α = 0 , f 0 and G ( u ) = g 1 exp ( m u )
X 1 + ν X 2 + μ X 3 , X 4 + λ X 2 + σ X 3 , X 2 , X 3 .
3.
If α = f = 0 and G ( u ) = g 1 exp ( m u )
X 1 + ν X 2 + μ X 3 , X 5 + λ X 2 + σ X 3 , X 2 , X 3 .
4.
If α = 0 and G ( u ) = g 1 u
X 1 + ν X 2 + μ X 3 , X 6 + λ X 2 + σ X 3 , X 2 , X 3 .
5.
If α = f = 0 and G ( u ) = g 1 u n
X 1 + ν X 2 + μ X 3 , X 7 + λ X 2 + σ X 3 , X 2 , X 3 .
In the above, ν, μ, λ and σ are arbitrary constants.
Now, taking into account the optimal system of one-dimensional subalgebras for each maximal Lie algebra given in Theorem 3, we will consider some one-dimensional reductions which allow us to transform PDE (6) into PDEs in 1+1-dimensions. Moreover, since these PDEs admit symmetry groups, we can reduce the number of independent variables again to obtain third-order nonlinear ODEs.

3.1. Reduction under X 1 + ν X 2 + μ X 3

Let us start considering the symmetry generator X 1 + ν X 2 + μ X 3 , where ν and μ are constants. Using this symmetry, we reduce (6) to a PDE with two independent variables. The symmetry gives the invariants
z = x ν t , r = y μ t u = w ( z , r ) .
Using these invariants, Equation (6) reduces to
γ ν w z z z + w r r z + γ μ w r r r + w r z z α w z z + w r r + G w z + w r ν w z μ w r f = 0 ,
which admits the symmetries
Z 1 = r , Z 2 = z .
The symmetry Z 1 + δ Z 2 , with δ constant, provides the invariants
q = z δ r , w = h ( q ) ,
and therefore PDE (11) transforms to the third-order nonlinear ODE
γ δ 2 + 1 ν δ μ h α δ 2 + 1 h + G 1 δ + δ μ ν h f = 0 .
Taking h ( q ) = V ( q ) , Equation (14) can be written as follows
V = 1 γ δ 2 + 1 ν δ μ α δ 2 + 1 V H 1 δ + δ μ ν V + f ,
where H ( ω ( V ) ) = G ( V d q ) = G ( h ) .

3.2. Reduction under X 4 + λ X 2 + σ X 3

Now, we consider the symmetry X 4 + λ X 2 + σ X 3 . This symmetry produces the invariants
z = x λ f m e f m t , r = y σ f m e f m t , u = f t + w ( z , r ) ,
where w ( z , r ) verifies
γ λ w z z z + w r r z + γ σ w r r r + w r z z + g 1 e m w w z + w r λ w z σ w r = 0 .
PDE (17) admits symmetries (12). Taking into account invariants (13), PDE (17) is transformed into the third-order nonlinear ODE
γ δ 2 + 1 λ δ σ h + g 1 e m h 1 δ + δ σ λ h = 0 .

3.3. Reduction under X 5 + λ X 2 + σ X 3

Consider the symmetry X 5 + λ X 2 + σ X 3 . This symmetry yields the invariants
z = x λ m log t , r = y σ m log t , u = w ( z , r ) 1 m log t ,
where w ( z , r ) satisfies
γ λ w z z z + w r r z + γ σ w r r r + w r z z + g 1 m e m w w z + w r λ w z σ w r 1 = 0 .
This equation admits symmetries (12). From Z 1 + δ Z 2 one obtains the invariants (13). Taking into account invariants (13), PDE (20) is transformed into the third-order nonlinear ODE
γ δ 2 + 1 λ δ σ h + g 1 m e m h 1 δ + δ σ λ h 1 = 0 .

3.4. Reduction under X 6 + λ X 2 + σ X 3

Consider the symmetry X 6 + λ X 2 + σ X 3 . This symmetry produces the invariants
z = x λ log t , r = y σ log t , u = t w ( z , r ) ,
where w ( z , r ) must satisfy
γ λ w z z z + w r r z + γ σ w r r r + w r z z γ w z z + w r r + g 1 w w z + w r λ w z σ w r + w f = 0 .
PDE (23) admits symmetries (12). By considering invariants (13), PDE (23) is transformed into the third-order nonlinear ODE
γ δ 2 + 1 λ δ σ h γ δ 2 + 1 h + g 1 h 1 δ + δ σ λ h + h f = 0 .

3.5. Reduction under X 7 + λ X 2 + σ X 3

Consider the symmetry X 7 + λ X 2 + σ X 3 . This symmetry yields the invariants
z = x λ n log t , r = y σ n log t , u = t 1 / n w ( z , r ) ,
where w ( z , r ) satisfies
γ λ w z z z + w r r z + γ σ w r r r + w r z z + γ w z z + w r r + g 1 n w n w z + w r λ w z σ w r w = 0 .
This equation admits symmetries (12). Taking into account invariants (13), PDE (26) is transformed into the third-order nonlinear ODE
γ δ 2 + 1 λ δ σ h + γ δ 2 + 1 h + g 1 n h n 1 δ + δ σ λ h h = 0 .

4. Exact Solutions via Kudryashov’s Method

In this section, we determine the function G for which Equation (6) admits solutions which are obtained by employing the Kudryashov’s method [8]. Kudryashov [6] called the simplest equation to any nonlinear ordinary differential equation of lesser order than the original equation with a known general solution.
  • The Riccati equation was the first example of the simplest equation
V + V 2 p V q = 0
If V ( q ) is a solution of Equation (28), then the equation
V = 2 V 3 3 p V 2 + ( p 2 2 q ) V + p q
has special solutions that are expressed via the general solution of Equation (28). It was proved by differentiating Equation (28) with respect to z and substituting V from Equation (28) into expressions obtained.
Equation (15), with H = γ δ 2 + 1 δ μ ν 1 δ ( 2 V 2 3 p V ) and α = 0 , can be written in the same form that Equation (29) by considering
p q = f γ δ 2 + 1 ν δ μ , p 2 2 q = 1 γ δ 2 + 1 .
Thus, Equation (15) has special solutions that are expressed via the general solution of Riccati Equation (28).
In order to apply the Kudryashov’s method to the nonlinear PDE in 1+1-dimensions (11), the first step is to reduce the nonlinear PDE into a nonlinear ODE, which we have already done using Lie symmetries in the previous section. Thus, we consider the ODE (15), which can be written in the form
V = a 1 H V + a 2 V + a 3 V + a 4 ,
where a 1 = δ 1 γ δ 2 + 1 ν δ μ , a 2 = 1 γ δ 2 + 1 , a 3 = α γ ν δ μ and a 4 = f γ δ 2 + 1 ν δ μ . We suppose that the solution of ODE (30) can be expressed in terms of a polynomial of the form
V = n = 0 N A n Y n ,
where Y = Y ( q ) satisfies the first-order nonlinear ODE
Y ( q ) = Y 2 ( q ) Y ( q ) ,
A n , n = 0 , , N , are constants to be determined later, A N 0 . We note that the solution of (32) is
1 1 + exp ( q ) .
Considering the homogeneous balance between V and H ( V ) V in (30), we obtain:
  • If H V = V 2 then N = 2 .
  • If H V = V 3 then N = 1 .
We consider H V = V 2 and N = 2 then we can write (31) as
V = A 0 + A 1 Y ( q ) + A 2 Y 2 ( q ) ,
with A 2 0 . In the following we determine A n , n = 0 , , 2 . We substitute V , V 2 , V and V expressions in Equation (30). Equating each coefficient of Y n , n = 0 , , 2 , to zero, yields a set of simultaneous algebraic equations for A n .
A 0 2 a 1 + A 0 a 2 + a 4 = 0 , A 1 ( 2 A 0 A 1 + A 2 + a 3 1 ) = 0 , 2 A 0 A 1 A 2 + A 1 A 1 2 + A 1 a 3 3 A 1 + A 2 A 2 + 2 A 2 , a 3 4 A 2 = 0 , 2 ( A 1 A 1 A 2 A 1 + A 2 a 3 5 A 2 ) = 0 , A 2 ( A 1 A 2 6 ) = 0 .
Solving system (35) for a 1 = 36 + a 2 2 4 a 4 and a 3 = 5 , we obtain the set of solutions:
A 2 = 24 a 4 36 + a 2 2 , A 1 = 48 a 4 36 + a 2 2 , A 0 = 2 ( 150 25 a 2 ) a 4 25 ( 36 + a 2 2 ) .
Consequently, the solution of Equation (30) with a 1 = 36 + a 2 2 4 a 4 and a 3 = 5 is
V ( q ) = 2 a 4 ( a 2 cosh ( q ) + a 2 6 sinh ( q ) + 6 ) a 2 2 36 ( cosh ( q ) + 1 ) .
  • The second example, presented by Kudryashov, was the Jacobi elliptic function equation
( V ) 2 V 4 a V 3 b V 2 c V d = 0 .
If V ( q ) is a solution of Equation (37) then the equation
V = 2 V 3 + 3 2 a V 2 + b V + 1 2 c ,
has special solutions that are expressed via the general solution of Equation (37). It was proved on similar lines given above.
Equation (15), with H ( V ) = γ δ 2 + 1 δ μ ν 1 δ ( 2 V 2 + 3 2 a V ) and α = 0 , can be written in the same form that of Equation (38) by considering
c = 2 f γ δ 2 + 1 ν δ μ , b = 1 γ δ 2 + 1 .
Thus, Equation (15) has special solutions that are expressed via the general solution of Jacobi Equation (37).
If V 1 , V 2 , V 3 and V 4 are the roots of
V 4 + a V 3 + b V 2 + c V + d = 0 ,
then Equation (37) with b and c given in (39) is
( V ) 2 = ( V V 1 ) ( V V 2 ) ( V V 3 ) ( V V 4 ) .
By a transformation [6], its solution could be written in terms of the Jacobi elliptic function sn(mq, k), where sn is the elliptic sine function.
  • The Weierstrass elliptic function equation was the third example of the simplest equation
( V ) 2 4 V 3 a V 2 b V c = 0 .
If V ( q ) is a solution of Equation (41), then the equation
V = 6 V 2 + a V + 1 2 b ,
has special solutions that are expressed via the general solution of Equation (41). It was proved on similar lines given above.
Equation (15), with H ( V ) = 6 γ δ 2 + 1 δ μ ν 1 δ V and α = 0 , can be written in the same form that Equation (42) by considering
b = 2 f γ δ 2 + 1 ν δ μ , a = 1 γ δ 2 + 1 .
Thus, Equation (15) has special solutions that are expressed via the general solution of Weierstrass Equation (41).

5. Conservation Laws

We construct conservation laws for the ( 2 + 1 ) -dimensional generalized BBMB Equation (6) by employing the multiplier method [23,34]. First, we determine low-order multipliers
Q = Q ( t , x , y , u )
admitted by the Equation (6). Recall that multipliers Q ( t , x , y , u ) for the equation under study (6) are obtained from the following equation known as the determining equation
δ δ u u t γ ( u t x x + u t y y ) α ( u x x + u y y ) + G ( u x + u y ) f Q = 0 ,
which holds off of the set of solutions of Equation (6) and where δ δ u is the Euler operator with respect to the variable u [20].
Splitting Equation (44) with respect to u and its derivatives, we obtain a determining system depending on the arbitrary parameters α , γ , f and the arbitrary function G ( u ) which Equation (6) involves. We solve the determining system by using “rifsimp” and “pdsolve” commands in Maple, and we determine the following classification for the low-order multipliers (43). We obtain the following result.
Theorem 4.
The low-order multipliers of differential order zero (43) admitted by the (2+1)-dimensional generalized BBMB Equation (6) depending on the arbitrary constants α , γ , f and an arbitrary function G ( u ) are given in the cases:
(i) 
For arbitrary f and G ( u ) , and α 0 , γ verifying c c γ α > 0 , with c 0 any fixed real constant, the admitted multipliers are:
Q 1 = exp c t + 1 2 c c γ α ( x y ) ,
Q 2 = exp c t 1 2 c c γ α ( x y ) .
(ii) 
For arbitrary f and G ( u ) , and α 0 , γ verifying c c γ α < 0 , with c 0 any fixed real constant, the admitted multipliers are:
Q 3 = e c t sin 1 2 c c γ α ( x y ) ,
Q 4 = e c t cos 1 2 c c γ α ( x y ) .
(iii) 
For α = 0 , γ > 0 , arbitrary f and G ( u ) , the admitted multipliers are:
Q 5 = f t + u ,
Q 6 = f 1 ( y x ) ,
Q 7 = f 2 ( t ) exp x y 2 γ ,
Q 8 = f 3 ( t ) exp x y 2 γ .
(iv) 
For α = 0 , γ < 0 , arbitrary f and G ( u ) , the admitted multipliers are Q 5 , Q 6 and
Q 9 = f 4 ( t ) sin x y 2 γ ,
Q 10 = f 5 ( t ) cos x y 2 γ .
In the above, f 1 ( y x ) , f 2 ( t ) , f 3 ( t ) , f 4 ( t ) and f 5 ( t ) are arbitrary functions of their arguments.
A local conservation law for the (2+1)-dimensional generalized BBMB Equation (6) is a divergence expression of the form
D t T ( t , x , y , u , u t , u x , ) + D x X ( t , x , y , u , u t , u x , ) + D y Y ( t , x , y , u , u t , u x , ) = 0 ,
that holds for the solutions u ( t , x , y ) of Equation (6), and where D t , D x and D y denote the total derivative operators with respect to t , x and y, respectively. T is a function called conserved density, and X , Y are functions called spatial fluxes. All these functions depend on t , x , y , u and derivatives of u.
Using the multipliers obtained previously (Theorem 4), we can calculate the conservation laws by integrating the following equation, called the characteristic equation,
D t T + D x X + D y Y = u t γ ( u t x x + u t y y ) α ( u x x + u y y ) + G ( u x + u y ) f Q ,
or via an homotopy formula [33]. We obtain the following result.
Theorem 5.
Low-order conservation laws (55) admitted by the (2+1)-dimensional generalized BBMB Equation (6) depending on the arbitrary constants α , γ and f and an arbitrary function G ( u ) are given by:
(i) 
For arbitrary f and G ( u ) , and α 0 , γ verifying c c γ α > 0 , with c 0 any fixed real constant, the admitted conservation laws are:
T 1 = α u c γ α exp c t + 1 2 c c γ α ( x y ) , X 1 = ( 2 γ c f + γ c u t + 2 α f 2 c ( c γ α ) α u x γ u t x + 2 2 α c c γ α u + G ( u ) d u ) exp c t + 1 2 c c γ α ( x y ) , Y 1 = c 2 ( c γ α ) γ u t + γ u t y + α u y + 2 2 α c c γ α u G ( u ) d u ) exp c t + 1 2 c c γ α ( x y ) .
T 2 = α u c γ α exp c t 1 2 c c γ α ( x y ) , X 2 = ( 2 γ c f γ c u t 2 α f 2 c ( c γ α ) α u x γ u t x 2 2 α c c γ α u + G ( u ) d u ) exp c t 1 2 c c γ α ( x y ) , Y 2 = c 2 ( c γ α ) γ u t γ u t y α u y + 2 2 α c c γ α u + G ( u ) d u ) exp c t 1 2 c c γ α ( x y ) .
(ii) 
For arbitrary f and G ( u ) , and α 0 , γ verifying c c γ α < 0 , with c 0 any fixed real constant, the admitted conservation laws are:
T 3 = α u e c t c γ α sin 1 2 c c γ α ( x y ) , X 3 = e c t α u x + γ u t x G ( u ) d u sin 1 2 c c γ α ( x y ) + e c t c ( c γ α ) α c u 2 α f + γ f c 2 γ c 2 2 u t cos 1 2 c c γ α ( x y ) , Y 3 = e c t 1 2 c c γ α α u + γ u t cos 1 2 c c γ α ( x y ) + e c t γ u t y + α u y G ( u ) d u sin 1 2 c c γ α ( x y ) .
T 4 = α u e c t c γ α cos 1 2 c c γ α ( x y ) , X 4 = e c t α u x + γ u t x G ( u ) d u cos 1 2 c c γ α ( x y ) + e c t c ( c γ α ) α c u 2 α f + γ f c 2 γ c 2 2 u t sin 1 2 c c γ α ( x y ) , Y 4 = e c t 1 2 c c γ α α u + γ u t sin 1 2 c c γ α ( x y ) + e c t γ u t y + α u y G ( u ) d u cos 1 2 c c γ α ( x y ) .
(iii) 
For α = 0 , γ > 0 , arbitrary f and G ( u ) , the admitted conservation laws are:
T 5 = γ u x 2 + u y 2 2 + f t u 2 2 , X 5 = γ f t u u t x G f t u d u , Y 5 = γ f t u u t y + G f t u d u .
T 6 = u f 1 ( y x ) 2 γ u d 2 d ( y x ) 2 f 1 ( y x ) , X 6 = γ u t d d ( y x ) f 1 ( y x ) f f 1 ( y x ) d x γ u t x f 1 ( y x ) + f 1 ( y x ) G ( u ) d u , Y 6 = γ u t d d ( y x ) f 1 ( y x ) γ u t y f 1 ( y x ) + f 1 ( y x ) G ( u ) d u .
T 7 = f exp x y 2 γ f 2 ( t ) d t , X 7 = f 2 ( t ) γ 2 u t γ u t x + G ( u ) d u exp x y 2 γ , Y 7 = f 2 ( t ) γ 2 u t γ u t y + G ( u ) d u exp x y 2 γ .
T 8 = f exp x y 2 γ f 3 ( t ) d t , X 8 = f 3 ( t ) γ 2 u t γ u t x + G ( u ) d u exp x y 2 γ , Y 8 = f 3 ( t ) γ 2 u t γ u t y + G ( u ) d u exp x y 2 γ .
(iv) 
For α = 0 , γ < 0 , arbitrary f and G ( u ) , the admitted conservation laws are:
T 9 = u x x + u y y γ u cos x y 2 γ , X 9 = γ f 4 t ( t ) α f 4 ( t ) u x sin x y 2 γ + γ 2 f 4 t ( t ) + α 2 γ u cos x y 2 γ + f 4 G d u sin x y 2 γ + f 2 γ f 4 ( t ) cos x y 2 γ , Y 9 = γ f 4 t ( t ) + α f 4 ( t ) u y sin x y 2 γ γ 2 f 4 t + α 2 γ f 4 u cos x y 2 γ + f 4 G d u sin x y 2 γ .
T 10 = u x x + u y y γ u sin x y 2 γ , X 10 = γ f 3 t ( t ) α f 3 ( t ) u x cos x y 2 γ γ 2 f 3 t ( t ) + α 2 γ u sin x y 2 γ + f 4 G d u cos x y 2 γ + f 2 γ f 4 ( t ) sin x y 2 γ , Y 10 = γ f 3 t ( t ) + α f 3 ( t ) u y cos x y 2 γ + γ 2 f 3 t + α 2 γ f 3 u sin x y 2 γ + f 3 G d u cos x y 2 γ .
In the above, f 1 ( y x ) , f 2 ( t ) , f 3 ( t ) , f 4 ( t ) and f 5 ( t ) are arbitrary functions of their arguments.

6. Line Soliton Solution

Consider G ( u ) = g 1 u n + g 2 , α = f = 0 , where n 2 , 1 , 0 , g 1 0 and g 2 are arbitrary constants, and the two-dimensional subalgebra spanned by
X = X 1 + λ X 2 and Y = X 3 μ X 2 ,
which satisfies X , Y = 0 . By applying the abelian subalgebra (67), the (2+1)-dimensional generalized BBMB Equation (6) can be reduced into a nonlinear third-order ODE through the use of two independent invariants z and h satisfying
X z = 0 , X h = 0 , Y z = 0 , Y h = 0 .
From (68), one obtains
z = x + μ y λ t , h = u .
The group invariant solution
u = h ( x + μ y λ t ) ,
is a line travelling wave. The amplitude h of a line travelling wave is invariant under translations in the perpendicular direction. The solution (70) depends on two parameters, μ = tan ϕ shows the direction of the wave propagation, i.e., the inclination of the line travelling wave in the ( x , y ) -plane, with ϕ the angle from the positive y-axis in counterclockwise direction, and where c = λ 1 + μ 2 is the speed of the line wave. We are interested in line travelling waves whose amplitude exhibits exponential asymptotic decay for large | z | . These line travelling waves are known as line solitons. The group-invariant solution (70) transforms Equation (6) into the nonlinear ODE given by
γ λ ( 1 + μ 2 ) h + ( 1 + μ ) ( g 1 h n + g 2 ) λ h = 0 .
For arbitrary n, λ and μ , ODE (71) only admits the obvious point symmetry (invariance under translations in z)
V = z .
The differential invariants corresponding to the infinitesimal generator V (72) are
ω = h and χ = h .
After substituting invariants (73) into the third-order ODE (71) one finds that χ ( ω ) verifies the second-order ODE
γ λ ( 1 + μ 2 ) χ χ + χ 2 + ( 1 + μ ) ( g 1 ω n + g 2 ) λ = 0 .
Fortunately, Equation (74) can be readily integrated to yield
γ λ ( 1 + μ 2 ) χ χ + ( 1 + μ ) g 1 n + 1 ω n + g 2 λ ω k 1 = 0 ,
with k 1 arbitrary constant. Integrating again Equation (75) with respect to ω , we obtain the general solution of Equation (74), which is given by
χ 2 = 1 γ λ ( 1 + μ 2 ) k 2 + 2 k 1 ω + λ g 2 ( 1 + μ ) ω 2 2 ( 1 + μ ) g 1 ( n + 1 ) ( n + 2 ) ω n + 2 .
We recall that we look for a localized smooth solution such that h , h , h 0 as z ± . This implies k 1 = k 2 = 0 . Undoing change of variables (73), one obtains
h = ± h 1 γ λ ( 1 + μ 2 ) λ ( 1 + μ ) g 2 2 ( 1 + μ ) g 1 ( n + 1 ) ( n + 2 ) h n .
Equation (77) is a separable first-order ODE. Suppose μ 1 , after separating and integrating, one obtains
h ( z ) = ( n + 1 ) ( n + 2 ) λ ( 1 + μ ) g 2 2 g 1 ( 1 + μ ) sech 2 n 2 λ ( 1 + μ ) g 2 γ λ ( 1 + μ 2 ) ( z + z 0 ) 1 / n ,
where z 0 is an arbitrary constant. Note that sech 2 ( ϕ ) is smooth in ϕ and vanishes as ϕ ± .
Thus, for
λ ( 1 + μ ) g 2 γ λ > 0
we have Equation (78) is a smooth solution in z which asymptotically decays to 0 as z ± . Finally, undoing change of variable (69) one obtains a line soliton solution for PDE (6) given by
u ( t , x , y ) = ( n + 1 ) ( n + 2 ) λ ( 1 + μ ) g 2 2 g 1 ( 1 + μ ) sech 2 n 2 λ ( 1 + μ ) g 2 γ λ ( 1 + μ 2 ) ( x + μ y λ t + z 0 ) 1 / n .
In Figure 1, it is represented the line soliton solution (79) considering μ = 1 , γ = 1 , λ = 2 , n = 1 , g 1 = 3 , g 2 = 1 and z 0 = 0 , for some fixed values of t. Here, it can be seen how the amplitude of the line soliton presents exponential asymptotic decay for large x + μ y λ t .

7. Conclusions

In this paper, we have classified the point symmetries admitted by the generalized family of BBMB equations (6). We have determined the most general symmetry Lie algebra along with its non-zero commutator structure that PDE (6) admits depending on three arbitrary parameters α , γ and f, and the arbitrary function G ( u ) . We have also constructed an optimal system of one-dimensional subalgebras for each maximal Lie algebra admitted by the considered equation. By using the subalgebras of this optimal system, we determine similarity variables and similarity solutions, which allow us to transform the generalized (2+1)-dimensional BBMB (6) into (1+1)-dimensional nonlinear third-order PDEs. Moreover, taking into account the symmetries of the reduced equations, we transform the considered family to nonlinear third-order ODEs. We have explored several exact solutions for this family. On the other hand, we have determined a classification of admitted low-order conservation laws of PDE (6) depending on the arbitrary parameters and the arbitrary function. Finally, line soliton solutions have been obtained in a particular case.

Author Contributions

Conceptualization, M.S.B., T.M.G.-L. and R.d.l.R.; methodology, M.S.B., T.M.G.-L. and R.d.l.R.; software, M.S.B., T.M.G.-L. and R.d.l.R.; validation, M.S.B., T.M.G.-L. and R.d.l.R.; formal analysis, M.S.B., T.M.G.-L. and R.d.l.R.; investigation, M.S.B., T.M.G.-L. and R.d.l.R.; writing—original draft preparation, M.S.B., T.M.G.-L. and R.d.l.R.; writing—review and editing, T.M.G.-L. and R.d.l.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This article has no additional data.

Acknowledgments

The authors acknowledge the financial support from Junta de Andalucía group FQM-201. The authors warmly thank the referees for their valuable comments and recommending changes that significantly improved this paper. In memory of María de los Santos Bruzón Gallego: thank you for dedicating your time and effort to care us and help us. You will always be our role model. May Maruchi rest in peace.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Benjamin, T.B.; Bona, J.L.; Mahony, J.J. Model Equations for Long Waves in Nonlinear Dispersive Systems. Phil. Trans. R. Soc. A 1972, 272, 47–78. [Google Scholar]
  2. Hajishafieiha, J.; Abbasba, S. A new class of polynomial functions for approximate solution of generalized Benjamin-Bona-Mahony-Burgers (gBBMB) equations. Appl. Math. Comput. 2020, 367, 124765. [Google Scholar] [CrossRef]
  3. Alazman, A.A.; Albert, J.P.; Bona, J.L.; Chen, M.; Wu, J. Comparisons between the BBM equation and a Boussinesq system. Adv. Differ. Equ. 2006, 11, 121–166. [Google Scholar]
  4. Bona, J.L.; Pritchard, W.G.; Scott, L.R. A Comparison of Solutions of Two Model Equations for Long Waves; Wisconsin Univ-Madison Mathematics Research Center: Madison, WI, USA, 1983. [Google Scholar]
  5. Kudryashov, N.A. Analytical Theory of Nonlinear Differential Equations; Institute of Computer Investigations: Moskow/Igevsk, Russia, 2004. [Google Scholar]
  6. Kudryashov, N.A. Simplest equation method to look for exact solutions of nonlinear differential equations. Chaos Solitons Fractals 2005, 24, 1217–1231. [Google Scholar] [CrossRef] [Green Version]
  7. Kudryashov, N.A.; Loguinova, N.B. Extended simplest equation method for nonlinear differential equations. Appl. Math. Comput. 2008, 205, 396–402. [Google Scholar] [CrossRef]
  8. Kudryashov, N.A. One method for finding exact solutions of nonlinear differential equations. Commun. Nonlinear Sci. Numer Simulat 2012, 17, 2248–2253. [Google Scholar] [CrossRef] [Green Version]
  9. Gómez, C.A.; Salas, A.H.; Frias, B.A. New periodic and soliton solutions for the Generalized BBM and Burgers-BBM equations. Appl. Math. Comput. 2010, 217, 1430–1434. [Google Scholar] [CrossRef]
  10. Seadawy, A.R.; Sayed, A. Traveling wave solutions of the Benjamin-Bona-Mahony water wave equations. Abstr. Appl. Anal. 2014, 926838. [Google Scholar] [CrossRef] [Green Version]
  11. Paliathanasis, A. Lie symmetries and singularity analysis for generalized shallow-water equations. Int. J. Nonlinear Sci. Numer. Simul. 2020, 21, 739–747. [Google Scholar] [CrossRef]
  12. Yuping, Z.; Jing, L.; Guangmei, W. Lax pair, auto-Bäcklund transformation and conservation law for a generalized variable-coefficient KdV equation with external-force term. Appl. Math. Lett. 2015, 45, 58–63. [Google Scholar]
  13. Adomian, G. Nonlinear Stochastic Operator Equations; Kluwer Academic Publishers: Boston, MA, USA, 1986. [Google Scholar]
  14. Adomian, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers: Boston, MA, USA, 1994. [Google Scholar]
  15. Dhaigude, D.B.; Birajdar, G.A.; Nikam, V.R. Adomain decomposition method for fractional Benjamin-Bona-Mahony-Burger’s equations. Int. J. Appl. Math. Mech. 2012, 8, 42–51. [Google Scholar]
  16. Zeidan, D.; Chau, C.K.; Lu, T.T.; Zheng, W.Q. Mathematical studies of the solution of Burgers’ equations by Adomian decomposition method. Math. Meth. Appl. Sci. 2019, 43, 2171–2188. [Google Scholar] [CrossRef]
  17. Lu, D.; Seadawy, A.R.; Iqbal, M. Construction of new solitary wave solutions of generalized Zakharov-Kuznetsov-Benjamin-Bona-Mahony and simplified modified form of Camassa-Holm equations. Open Phys. 2018, 16, 896–909. [Google Scholar] [CrossRef]
  18. Simbanefayi, I.; Khalique, C.M. Travelling wave solutions and conservation laws for the Korteweg-de Vries-Benjamin-Bona-Mahony equation. Results Phys. 2018, 8, 57–63. [Google Scholar] [CrossRef]
  19. Bluman, G.W.; Anco, S.C. Symmetry and Integration Methods for Differential Equations; Springer: New York, NY, USA, 2002. [Google Scholar]
  20. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  21. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  22. Clarkson, P.A. Nonclassical symmetry reductions of the Boussinesq equation. Chaos Solitons Fractals 1995, 5, 2261–2301. [Google Scholar] [CrossRef]
  23. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  24. Ibragimov, N.H. A new conservation theorem. J. Math. Anal. Appl. 2007, 333, 311–328. [Google Scholar] [CrossRef] [Green Version]
  25. Ibragimov, N.H. Quasi-self-adjoint differential equations. Arch. ALGA 2007, 4, 55–60. [Google Scholar]
  26. Ibragimov, N.H. Nonlinear self-adjointness and conservation laws. J. Phys. A Math. Theor. 2011, 44, 432002. [Google Scholar] [CrossRef]
  27. Galiakberova, L.R.; Ibragimov, N.H. Nonlinear self-adjointness of the Krichever-Novikov equation. Commun. Nonlinear Sci. Numer. Simulat. 2014, 19, 361–363. [Google Scholar] [CrossRef]
  28. Bruzón, M.S.; Gandarias, M.L.; Ibragimov, N.H. Self-adjoint sub-classes of generalized thin film equations. J. Math. Anal. Appl. 2009, 357, 307–313. [Google Scholar] [CrossRef]
  29. Ibragimov, N.H.; Torrisi, M.; Tracina, R. Quasi self-adjoint nonlinear wave equations. J. Phys. A Math. Theor. 2010, 43, 442001. [Google Scholar] [CrossRef] [Green Version]
  30. Ibragimov N., H.; Torrisi, M.; Tracina, R. Self-adjointness and conservation laws of a generalized Burgers equation. J. Phys. A Math. Theor. 2011, 44, 145201. [Google Scholar] [CrossRef] [Green Version]
  31. Anco, S.C. Conservation laws of scaling-invariant field equations. J. Phys. A Math. Gen. 2003, 36, 8623–8638. [Google Scholar] [CrossRef] [Green Version]
  32. Anco, S.C.; Bluman, G.W. Direct Construction of Conservation Laws from Field Equations. Phys. Rev. Lett. 1997, 78, 2869–2873. [Google Scholar] [CrossRef]
  33. Anco, S.C.; Bluman, G.W. 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] [CrossRef] [Green Version]
  34. Anco, S.C.; Bluman, G.W. Direct construction method for conservation laws of partial differential equations. Part II: General treatment. Eur. J. Appl. Math. 2002, 13, 567–585. [Google Scholar] [CrossRef] [Green Version]
  35. Anco, S.C. Generalization of Noether’s theorem in modern form to non-variational partial differential equations. In Recent progress and Modern Challenges in Applied Mathematics, Modeling and Computational Science; Fields Institute Communications 79; Springer: New York, NY, USA, 2017. [Google Scholar]
  36. Anco, S.C. Symmetry properties of conservation laws. Internat. J. Modern Phys. B 2016, 30, 1640003. [Google Scholar] [CrossRef] [Green Version]
  37. Anco, S.C.; Kara, A. Symmetry-invariance conservation laws of partial differential equations. Eur. J. Appl. Math. 2018, 29, 78–117. [Google Scholar]
  38. Anco, S.C. On the incompleteness of Ibragimov’s conservation law theorem and its equivalence to a standard formula using symmetries and adjoint-symmetries. Symmetry 2017, 9, 33. [Google Scholar] [CrossRef] [Green Version]
  39. Anco, S.C.; Gandarias, M.L. Symmetry multi-reduction method for partial differential equations with conservation laws. Commun. Nonlinear Sci. Numer. Simulat. 2020, 91, 105349. [Google Scholar] [CrossRef]
  40. Kara, A.H.; Mahomed, F.M. Relationship between symmetries and conservation laws. Int. J. Theor. Phys. 2000, 39, 23–40. [Google Scholar] [CrossRef]
  41. Sjöberg, A. Double reduction of PDEs from the association of symmetries with conservation laws with applications. Appl. Math. Comput. 2007, 184, 608–616. [Google Scholar] [CrossRef]
  42. Sjöberg, A. On double reduction from symmetries and conservation laws. Nonlinear Anal. Real World Appl. 2009, 10, 3472–3477. [Google Scholar] [CrossRef]
  43. Bruzón, M.S.; Garrido, T.M.; de la Rosa, R. Conservation laws and exact solutions of a Generalized Benjamin–Bona–Mahony–Burgers equation. Chaos Solitons Fractals 2016, 89, 578–583. [Google Scholar] [CrossRef]
  44. Patera, J.; Winternitz, P. Subalgebras of real three- and four-dimensional Lie algebras. J. Math. Phys. 1977, 18, 1449. [Google Scholar] [CrossRef]
Figure 1. Line soliton solutions (79) for some fixed values of t, with μ = 1 , γ = 1 , λ = 2 , n = 1 , g 1 = 3 , g 2 = 1 and z 0 = 0 : [a] t = 0 (blue), [b] t = 5 (green), [c] t = 10 (purple).
Figure 1. Line soliton solutions (79) for some fixed values of t, with μ = 1 , γ = 1 , λ = 2 , n = 1 , g 1 = 3 , g 2 = 1 and z 0 = 0 : [a] t = 0 (blue), [b] t = 5 (green), [c] t = 10 (purple).
Symmetry 13 02083 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bruzón, M.S.; Garrido-Letrán, T.M.; de la Rosa, R. Symmetry Analysis, Exact Solutions and Conservation Laws of a Benjamin–Bona–Mahony–Burgers Equation in 2+1-Dimensions. Symmetry 2021, 13, 2083. https://doi.org/10.3390/sym13112083

AMA Style

Bruzón MS, Garrido-Letrán TM, de la Rosa R. Symmetry Analysis, Exact Solutions and Conservation Laws of a Benjamin–Bona–Mahony–Burgers Equation in 2+1-Dimensions. Symmetry. 2021; 13(11):2083. https://doi.org/10.3390/sym13112083

Chicago/Turabian Style

Bruzón, María S., Tamara M. Garrido-Letrán, and Rafael de la Rosa. 2021. "Symmetry Analysis, Exact Solutions and Conservation Laws of a Benjamin–Bona–Mahony–Burgers Equation in 2+1-Dimensions" Symmetry 13, no. 11: 2083. https://doi.org/10.3390/sym13112083

APA Style

Bruzón, M. S., Garrido-Letrán, T. M., & de la Rosa, R. (2021). Symmetry Analysis, Exact Solutions and Conservation Laws of a Benjamin–Bona–Mahony–Burgers Equation in 2+1-Dimensions. Symmetry, 13(11), 2083. https://doi.org/10.3390/sym13112083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop