Next Article in Journal
Multiple Novels and Accurate Traveling Wave and Numerical Solutions of the (2+1) Dimensional Fisher-Kolmogorov- Petrovskii-Piskunov Equation
Next Article in Special Issue
Lagrangian Formulation, Conservation Laws, Travelling Wave Solutions: A Generalized Benney-Luke Equation
Previous Article in Journal
Research on Vibration Suppression Method Based on Coaxial Stacking Measurement
Previous Article in Special Issue
Generalized Camassa–Holm Equations: Symmetry, Conservation Laws and Regular Pulse and Front Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Solutions and Conserved Vectors of the Two-Dimensional Generalized Shallow Water Wave Equation

by
Chaudry Masood Khalique
1,2,*,† and
Karabo Plaatjie
1,†
1
International Institute for Symmetry Analysis and Mathematical Modelling, Department of Mathematical Sciences, Mafikeng Campus, North-West University, Private Bag X 2046, Mmabatho 2735, South Africa
2
Department of Mathematics and Informatics, Azerbaijan University, Jeyhun Hajibeyli Str., 71, Baku AZ1007, Azerbaijan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(12), 1439; https://doi.org/10.3390/math9121439
Submission received: 26 March 2021 / Revised: 1 May 2021 / Accepted: 10 May 2021 / Published: 21 June 2021

Abstract

:
In this article, we investigate a two-dimensional generalized shallow water wave equation. Lie symmetries of the equation are computed first and then used to perform symmetry reductions. By utilizing the three translation symmetries of the equation, a fourth-order ordinary differential equation is obtained and solved in terms of an incomplete elliptic integral. Moreover, with the aid of Kudryashov’s approach, more closed-form solutions are constructed. In addition, energy and linear momentum conservation laws for the underlying equation are computed by engaging the multiplier approach as well as Noether’s theorem.

1. Introduction

It is widely known that most physical phenomena of the natural world are governed by nonlinear partial differential equations (NPDEs). Such equations arise in several physical and other problems such as in the study of electromagnetism, electrodynamics, fluid flow, elasticity, propagation of heat or sound, quantum mechanics, meteorology, and oceanography, just to mention a few. Some of these models include the modified Zakharov–Kuznetsov model [1] that recounts the ion-acoustic drift solitary waves existing in a magnetoplasma with electron–positron ions, which are found in a primordial universe. A generalized system of three-dimensional modified Kadomtsev–Petviashvili–Burgers-type equation, which describes ion-acoustic and dust-magneto-acoustic waves in one of cosmic or laboratory dusty plasmas, was studied in [2]. Moreover, in [3], the vector bright solitons of the Fokas–Lenells system, which models the femtosecond pulses in an optical fibre, was investigated. In addition, the Boussinesq–Burgers-type system of equations, which delineates shallow water waves appearing close to lakes or ocean beaches, was studied in [4]. The list continues; see also [5,6,7,8,9,10,11,12,13,14,15].
To obtain a good understanding of physical phenomena that are governed by NPDEs, it is imperative to look for their closed-form solutions. However, this is an arduous exercise, and for several years, researchers and scientists have been working on developing methods for finding closed-form solutions.
A good number of methods for obtaining closed-form solutions of NPDEs have been proposed. We have, amongst others, the following various special methods for finding closed-form solutions of NPDEs: the inverse scattering transform technique [16], the bifurcation method [17], the simplest equation technique [18], the extended simplest equation technique [19], Kudryashov’s technique [20], Hirota’s technique [21], Bäcklund transformation [22], Darboux transformation [23], the homogeneous balance method [24], ( G / G ) —expansion technique [25], Lie symmetry method [26,27,28], etc.
It was during the halfway point of the nineteenth century that Sophus Lie, a Norwegian mathematician, endowed the study of continuous transformation groups, presently called Lie groups. These groups have had an intense impact on pure and applied mathematics in addition to engineering, physics, and other applied sciences. It provides the most effective and powerful techniques for obtaining closed-form solutions to NPDEs. For example, see [26,27,28,29,30,31,32,33]. It is worth noting here that the notion of a Galois group had influenced Lie’s work on differential equations (DEs).
In the study of DEs, conservation laws play a crucial part. Conservation laws can be used to check whether a given differential equation is integrable. The accuracy of numerical solution methods can be checked by invoking conservation laws. They can be utilized in constructing solutions of partial differential equations (PDEs) by reducing their order. Conservation laws refer to certain physical properties that do not change subsequently within an isolated physical system. For example, conservation of energy, linear momentum, angular momentum, and electric charge. Recently, exact solutions were obtained for some PDEs using conservation laws. For differential equations that arise from the variational principle, there exists a connection between the symmetries of differential equations and conservation laws. This fact was established by Emmy Noether, a German mathematician, in 1918 and is stated in Noether’s theorem. For more details, see for example [34,35,36,37,38,39,40,41,42,43,44,45] and the references therein.
In this work, we study the generalized two-dimensional generalized shallow water wave equation, given by [46]
u t x + α u y u x x + 2 α u x u x y + β u x y + ν u x x x y = 0 ,
where α , β , and ν are nonzero constants. Here, t , x , and y represent, respectively, time and two scaled space coordinates whereas u ( t , x , y ) represents the height of the Riemann wave.
In [46], the authors presented painlévé analysis, nonlocal symmetries, and exact solutions of a variable coefficients form of (1), whereas in [47], the binary Bell polynomials technique was employed to a variable coefficient form of Equation (1). We note that, in the study of fluids, this equation decomposes to well-known equations when taking special cases of constants α , β , and ν . For example, when α = 2 , β = 0 , and ν = 1 , it becomes a ( 2 + 1 ) -dimensional breaking soliton [48]:
u t x 2 u y u x x 4 u x u x y u x x x y = 0 ,
which is one of the models used in describing the two-dimensional interplay of a Riemann wave distributed along the y-axis with a long wave generated along the x-axis. Wazwaz [49] obtained solutions of this equation using the modified form of Hirota’s bilinear method. The authors in [48] presented solutions of Equation (2) using the symmetry method. In the case where α = 2 , β = 0 , and ν = 1 , Equation (1) transforms to a ( 2 + 1 )-dimensional Calogero–Bogoyavlenskii–Schiff [50]
u t x + 2 u y u x x + 4 u x u x y + u x x x y = 0 ,
which describes the exchange of a long generating wave in the x direction with a Riemann generating wave in the y direction. Moreover, the researcher indicated that Equation (3) possesses soliton as well as N-soliton solutions, which are smooth in one coordinate. In [51], Wazwaz invoked Hirota’s bilinear approach to secure solutions of (3) and further utilized the tanh-coth method to gain its travelling wave solutions in [52]. The authors in [53] derived a Bäcklund transformation alongside a Lax pair of the form in the Calogero–Bogoyavlenskii–Schiff Equation (3) via singular manifold method. Furthermore, in [54], the authors invoked a tanh function as well as improved the ( G / G ) —expansion methods to construct travelling wave solutions of (3). It was also revealed that some analytic solutions secured via these methods are analogous. However, a generalized version of (3), namely
u t x + a u x u x y + b u y u x x + u x x x y = 0 ,
has been studied by many researchers by simply replacing 4 by a and 2 by b in (3). For instance, the ( G / G ) —expansion technique was invoked to achieve some travelling wave solutions of (4) in [55]. Moreover, in [56], the authors obtained breather-type and analytic soliton solutions of (4) by utilizing Hirota’s bilinear method and the extended homoclinic test technique. Using the idea of a tanh-coth approach, some closed-form solutions of (4) were obtained in [57]. Furthermore, sine-cosine and complex techniques were engaged to secure, respectively, analytic and meromorphic exact solutions of (4) in [58,59].
The formation of this paper is as follows: In Section 2, first, we employ Lie’s theory to reduce (1) to an ordinary differential equation (ODE). Closed-form solutions are then derived by an integration of the ODE and by invoking Kudryashov’s method. Thereafter, in Section 3, we construct conserved densities and spatial fluxes for (1) by invoking the multiplier method and Noether’s theorem. Finally, we present concluding remarks in Section 4.

2. Exact Solutions

In this section, we derive exact solutions of the two-dimensional generalized shallow water wave Equation (1). We begin by deriving symmetries of (1) and thereafter use them to obtain an ODE. We then employ the direct integration procedure together with Kudryashov’s method to construct closed-form solutions of (1).

2.1. Symmetries and Symmetry Reductions

Consider the symmetry group of (1) brought about by
U = T t + X x + Y y + η u ,
where all T, X, Y, and η are functions of variables t, x, y, and u. Using the fourth prolongation pr ( 4 ) U in (1) and with the help of Maple, we obtain an overdetermined system of linear homogeneous PDEs:
T u = 0 , T x = 0 , T y = 0 , X y = 0 , Y u = 0 , Y x = 0 , Y t t = 0 , T t t t = 0 , X u = 0 , Y y y = 0 , 2 η u + T t Y y = 0 , 2 α η x + β T t β Y y Y t = 0 , α η y + X t = 0 , 2 Y t y T t t = 0 , 2 X x T t + Y y = 0 .
The above linear PDEs may be easily solved and yield the following eight Lie point symmetries:
U 1 = t , U 2 = x , U 3 = y , U 4 = u , U 5 = 2 α t y + x u , U 6 = 2 α t t + α x x β x + α u u , U 7 = α x x 2 α y y α u + β x u , U 8 = 2 α t 2 t + α t x x + 2 α t y y β t x x y + α t u u .
Consider the symmetry U = U 1 + a U 2 + b U 3 , where a and b are constants. Using this symmetry U, we reduce (1) to a PDE with one less independent variable, that is, two independent variables. The symmetry U gives three invariants:
f = x a t , g = y b t , θ = u .
Using these invariants, Equation (1) reduces to
ν θ f f f g + β θ f g + α θ g θ f f + 2 α θ f θ f g b θ f g a θ f f = 0 ,
which has
R 1 = f , R 2 = g , R 3 = α H ( g ) g + a H ( g ) θ , R 4 = α f f α g g + b f β f α θ θ
as its symmetries. The symmetry R = R 1 + c R 2 , with c a constant, provides two invariants q = g c f and U = θ . These invariants transform (7) into a nonlinear ordinary differential equation (NLODE)
k 1 U ( q ) k 2 U ( q ) U ( q ) + k 3 U ( q ) = 0 ,
where k 1 = b a c 2 β , k 2 = 3 α c , k 3 = ν c 2 and q = a c b t c x + y .

2.2. Solution Using Direct Integration

In this subsection, we derive a solution of the two-dimensional generalized shallow water Equation (1) by direct integration of the ODE (8). Taking U ( q ) = V ( q ) , Equation (8) becomes
k 1 V k 2 V V + k 3 V = 0 .
Twice integration of (9) with respect to q gives
1 2 k 1 V 2 1 6 k 2 V 3 + 1 2 k 3 V 2 + r 1 V + r 2 = 0
with r 1 and r 2 constants. Equation (10) can be rewritten as
V 2 = k 2 3 k 1 V 3 k 3 k 1 V 2 2 r 1 k 1 V 2 r 2 k 1 .
Assume that n 1 , n 2 , and n 3 are roots of
V 3 3 k 3 k 2 V 2 6 r 1 k 2 V 6 r 2 k 2 = 0
with n 1 > n 2 > n 3 . Equation (11) now becomes
V 2 = k 2 3 k 1 ( V n 1 ) ( V n 2 ) ( V n 3 )
and its solution could be written in terms of the Jacobi elliptic function [60,61,62]
V ( q ) = n 2 + ( n 1 n 2 ) cn 2 k 2 ( n 1 n 3 ) 12 k 1 q , R 2 , R 2 = n 1 n 2 n 1 n 3 ,
where cn is the elliptic cosine function. Integrating Equation (12) with respect to q and returning to the original variables, we accomplish the solution of (1) as
u ( t , x , y ) = 12 k 1 n 1 n 2 2 k 2 ( n 1 n 3 ) R 8 EllipticE sn k 2 ( n 1 n 3 ) 12 k 1 q , R 2 , R 2 + n 2 n 1 n 2 1 R 4 R 4 q + K ,
with q = a c b t c x + y , K being a constant, and EllipticE [ g , k ] being the incomplete elliptic integral given as [60]
EllipticE [ g , k ] = 0 g 1 k 2 s 2 1 s 2 d s .
The wave profile of the periodic solution (13) for parametric values k 1 = 70 , k 2 = 10 , a = 4 , b = 0.2 , c = 0.6 , K = 1 , n 1 = 100 , n 2 = 50.05 , and n 3 = 60 at t = 14 , can be seen in Figure 1.

2.3. Solution via Kudryashov’s Method

In this subsection, we present the solution of (1), which is obtained by employing Kudryashov’s method [20]. This method is one of the most popular approaches for determining closed-form solutions of NPDEs. The first step is to reduce the NPDE (1) to a nonlinear ODE, which we already performed using the Lie symmetries in the previous section. Thus, we work with the ODE (8). We suppose that a solution of (8) can be expressed as
U ( q ) = n = 0 N A n Y n ( q ) ,
where Y ( q ) satisfies the first-order nonlinear ODE
Y ( q ) = Y 2 ( q ) Y ( q ) .
We note that the solution of (16) is
Y ( q ) = 1 1 + exp q .
For Equation (8), the balancing procedure yields N = 1 . Thus, from (15), we have
U ( q ) = A 0 + A 1 Y ( q ) .
Now substituting (18) into (8) and using (16), we obtain
2 b A 1 c Y 3 ( q ) 2 a A 1 c 2 Y 3 ( q ) + 3 a A 1 c 2 Y 2 ( q ) 3 b A 1 c Y 2 ( q ) A 1 c 2 Y ( q ) a + b A 1 c Y ( q ) + 6 α A 1 2 c 2 Y 5 ( q ) 15 α A 1 2 c 2 Y 4 ( q ) + 12 α A 1 2 c 2 Y 3 ( q ) 3 α A 1 2 c 2 Y 2 ( q ) 2 β A 1 c Y 3 ( q ) + 3 β A 1 c Y 2 ( q ) β A 1 c Y q 24 ν A 1 c 3 Y 5 ( q ) + 60 ν A 1 c 3 Y 4 ( q ) 50 ν A 1 c 3 Y 3 ( q ) + 15 ν A 1 c 3 Y 2 ( q ) ν A 1 c 3 Y ( q ) = 0 .
Separating on the powers of Y ( q ) gives the algebraic equations for the coefficients A 0 and A 1 as
α A 1 2 c 2 4 ν A 1 c 3 = 0 , b c A 1 ν A 1 c 3 a c 2 A 1 β A 1 c = 0 , β A 1 c α A 1 2 c 2 + 5 ν A 1 c 3 + a c 2 A 1 b c A 1 = 0 , 6 α A 1 2 c 2 25 ν A 1 c 3 a c 2 A 1 + b c A 1 β A 1 c = 0 .
The solution of these equations is
A 0 = A 0 , A 1 = 4 ν c α , a = b β c 2 ν c .
Thus, the solution of (1) reads
u ( t , x , y ) = A 0 + 4 ν c α 1 + exp β c 2 ν t c x + y .
The wave profile of the solution (19) for parametric values A 0 = 0 , α = 1 , β = 0.05 , c = 0.01 , ν = 1 , t = 0 is presented in Figure 2.
Invariant solution using U 5
The Lie point symmetry U 5 produces the invariants I 1 = t , I 2 = x , and I 3 = u x y / ( 2 α t ) . The use of these invariants furnishes us with the group invariant solution of (1) given by
u ( t , x , y ) = x y 2 α t β 2 α x C x + h ( t ) ,
where C is a constant and h is a function of t.
Invariant solution using U 6
The point symmetry U 6 of (1) yields the group invariant solution:
u ( t , x , y ) = 1 t G ( f , g ) β 2 α x , f = x t , g = y .
The substitution of this expression for u into (1) gives the PDE
2 ν G f f f g + 4 α G f G f g + 2 α G g G f f f G f f 2 G f = 0 ,
for which the symmetries are
Γ 1 = g , Γ 2 = G , Γ 3 = 2 α f + g g , Γ 4 = f f 2 g g G G .
The symmetry Γ 3 gives the invariant solution G ( f , g ) = U ( z ) + f g / ( 2 α ) , where z = g . This expression for G satisfies the Equation (22) for arbitrary function U. Thus, the group-invariant solution under U 6 of (1) is
u ( t , x , y ) = x y 2 α t β 2 α x U ( y ) t ,
where U ( y ) is an arbitrary function of y.
Invariant solution using U 7
Here, we use the symmetry operator U 7 to obtain the invariant solution associated with this symmetry. The usual procedure gives
u ( t , x , y ) = 1 x U ( ξ ) β 2 α x , ξ = x 2 y .
Substituting this expression of u into Equation (1) gives the nonlinear ODE
α U 2 + 6 α ξ U U 2 α U U + 3 ν U + 4 ν ξ 2 U + 12 ν ξ U = 0 .
Invariant solution using U 8
Finally, the use of symmetry U 8 gives the invariant solution
u ( t , x , y ) = 1 t Φ ( f , g ) + x y 2 α t β x 2 α ,
with f = x / t 1 / 2 and g = y / t , which when substituted in (1) gives the reduced PDE
α Φ f f Φ g + 2 α Φ f Φ f g + ν ϕ f f f g = 0 .
The above PDE (28) has three Lie symmetries
Γ 1 = f + F 1 ( g ) g , Γ 2 = F 2 ( g ) g + Φ , Γ 3 = f f + F 3 ( g ) g Φ Φ .
Now utilizing Γ 1 with F 1 ( g ) = 1 , we perform a reduction of Equation (28). We obtain the invariant solution Φ ( f , g ) = U ( z ) , where z = g f . This reduces Equation (28) to an ODE
ν U 3 α U U = 0 ,
for which the solution is
U ( z ) = 2 8 / 3 ν α 4 / 3 ζ α 2 / 3 2 4 / 3 ν z + C 0 ; s 1 , s 2 , s 1 = 2 5 / 2 α 3 C 0 , s 2 = 2 C 1 ,
where C 0 , C 1 are arbitrary constants and ζ is the Weierstrass zeta function [60]. Thus, the solution of (1) is
u ( t , x , y ) = 2 8 / 3 ν α 4 / 3 t ζ α 2 / 3 2 4 / 3 ν y t x t + C 0 ; s 1 , s 2 + x y 2 α t β x 2 α .
For the operator Γ 2 with F 2 ( g ) = 1 , we obtain the invariant solution Φ ( f , g ) = U ( f ) + g , which reduces (28) to the ODE U f f = 0 , for which the solution is U = C 1 f + C 2 , where C 1 , C 2 are constants. Reverting to the original variables, we attain the solution of (1) given by
u ( t , x , y ) = C 1 x t + C 2 t + y t 3 / 2 + x y 2 α t β x 2 α .
The dynamic behaviour of the solution (32), for parametric values α = 10 , β = 5 , C 1 = 0.3 , C 2 = 0.01 , and x = 0.5 , is sketched in Figure 3.

3. Conservation Laws

We now construct conservation laws for the two-dimensional generalized constant-coefficient shallow water Equation (1). To achieve this task, we utilize two methods; the first one is the multiplier method, and the second is classical Noether’s theorem.

3.1. Conservation Laws Utilizing the Multiplier Method

We determine the conserved vectors of Equation (1) by employing the multiplier method [26]. Here, we seek first-order multipliers Q = Q ( t , x , y , u , u t , u x , u y ) , which are established from
δ δ u Q u t x + α u y u x x + 2 α u x u x y + β u x y + ν u x x x y = 0 ,
where δ / δ u is the Euler operator [37]
δ δ u = u + s = 1 ( 1 ) s D i 1 D i s u i 1 i s
and for our purpose
δ δ u = u D x u x D y u y + D t D x u t x + D x 2 u x x + D x D y u x y + D x 3 D y u x x x y ,
with D t , D x , and D y being the total derivative operators given by
D t = t + u t u + u t t u t + u t x u x + u t y u y + . D x = x + u x u + u t x u t + u x x u x + u x y u y + . D y = y + u y u + u t y u t + u x y u x + u y y u y + .
Expanding (33) and separating on the appropriate derivatives of u, we obtain twenty-six multiplier-determining equations:
Q t x + α u x Q x y + β Q x y = 0 , Q t u + α Q x y = 0 , Q t u t + 4 Q u = 0 , Q x u t = 0 , Q x y = 0 , Q x u x Q u = 0 , Q t u x 2 α u y Q u + α Q y = 0 , Q t u y 2 α u x Q u + 2 α Q x 2 β Q u = 0 , Q x x = 0 , Q u x = 0 , Q y y = 0 , Q u y = 0 , Q y u t = 0 , Q y u x = 0 , Q y u y 2 Q u = 0 , Q u u = 0 , Q u u t = 0 , Q u u x = 0 , Q u u y = 0 , Q u t u t = 0 , Q u t u x = 0 , Q u t u y = 0 , Q u x u x = 0 , Q u x u y = 0 , Q u y u y = 0 , Q x u y = 0 ,
for which the solution is
Q = y F ( t ) + G ( t ) α u x F ( t ) 2 K 5 α t 2 u t K 5 α t x u x 2 K 5 α t y u y K 5 α t u + 2 K 6 β t u y 2 K 1 α t u y K 5 β t x + 4 K 6 t u t + K 6 x u x + 2 K 6 y u y + K 5 x y + K 6 u + K 2 u t + K 3 u x + K 4 u y + K 1 x ,
where F and G are arbitrary functions of t and K i , where i = 1 , , 6 , are constants of integration. The conserved quantities of Equation (1) are constructed using the divergence identity
D t T t + D x T x + D y T y = Q u t x + α u y u x x + 2 α u x u x y + β u x y + ν u x x x y ,
where T t represents the conserved density and T x and T y are spatial fluxes. Thus, after some reckoning, we obtain low-order conservation laws corresponding to the eight multipliers and these are given below.
Case 1. For Q 1 = x 2 α t u y , we have the conserved vector
T 1 t = 1 2 α t u u x y 1 2 α t u x u y + 1 2 x u x 1 2 u , T 1 x = 2 3 α 2 t u u x u y y + 2 3 α 2 t u u x y u y 4 3 α 2 t u x u y 2 + 3 4 α ν t u u x x y y + 1 2 α β t u u y y 1 2 α ν t u x u x y y + 1 4 α ν t u x x u y y α ν t u x x y u y + 1 2 α ν t u x y 2 1 2 α β t u y 2 + 1 2 α t u u t y 1 2 α t u t u y + α x u x u y + 1 2 β x u y + 3 4 ν x u x x y 1 2 ν u x y + 1 2 x u t , T 1 y = 2 3 α 2 t u u x u x y 2 3 α 2 t u u x x u y 2 3 α 2 t u x 2 u y 3 4 α ν t u u x x x y 1 2 α β t u u x y 1 4 α ν t u x u x x y 1 2 α β t u x u y + 1 4 α ν t u x x u x y 1 4 α ν t u x x x u y α t u u t x + 1 2 α x u x 2 1 2 α u u x + 1 2 β x u x 1 2 β u + 1 4 ν x u x x x 1 4 ν u x x ;
Case 2. For Q 2 = u t , we have the conserved vector
T 2 t = 2 3 α u u x u x y + 1 3 α u u x x u y + 1 2 ν u u x x x y + 1 2 β u u x y + 1 4 u u t x + 1 4 u t u x , T 2 x = 3 8 ν u u t x x y + 3 8 ν u t u x x y 1 4 ν u t x u x y + 1 8 ν u t x x u y + 1 4 ν u t x y u x 1 8 ν u t y u x x 1 3 α u u t x u y 1 3 α u u t y u x 1 4 β u u t y 1 4 u u t t + 2 3 α u t u y u x + 1 4 β u t u y + 1 4 u t 2 , T 2 y = 1 3 α u t u x 2 1 3 α u u t x u x + 1 4 β u t u x + 1 8 ν u t x x u x 1 8 ν u u t x x x + 1 8 ν u t u x x x 1 8 ν u t x u x x 1 4 β u u t x ;
Case 3. For Q 3 = u x , we have the conserved vector
T 3 t = 1 4 u x 2 1 4 u u x x , T 3 x = 2 3 α u x 2 u y + 1 3 α u u x u x y + 5 8 ν u x u x x y + 1 4 β u x u y + 1 4 u t u x + 1 8 ν u u x x x y 3 8 ν u x x u x y + 1 8 ν u x x x u y + 1 4 β u u x y + 1 4 u u t x , T 3 y = 1 3 α u x 3 + 1 4 β u x 2 1 3 α u u x u x x + 1 4 ν u x u x x x 1 8 ν u u x x x x 1 4 β u u x x 1 8 ν u x x 2 ;
Case 4. For Q 4 = u y , we have the conserved vector
T 4 t = 1 4 u x u y 1 4 u u x y , T 4 x = 2 3 α u x u y 2 + 1 4 β u y 2 1 3 α u u x y u y + 1 2 ν u x x y u y + 1 4 u t u y 3 8 ν u u x x y y + 1 4 ν u x u x y y 1 8 ν u x x u y y 1 4 ν u x y 2 1 3 α u u x u y y 1 4 β u u y y 1 4 u u t y , T 4 y = 1 3 α u u x u x y + 1 3 α u u x x u y + 3 8 ν u u x x x y + 1 4 β u u x y + 1 2 u u t x + 1 3 α u x 2 u y + 1 8 ν u x u x x y + 1 4 β u x u y 1 8 ν u x x u x y + 1 8 ν u x x x u y ;
Case 5. For Q 5 = 2 α t 2 u t α t x u x 2 α t y u y α t u β t x + x y , we have the conserved vector
T 5 t = 4 3 α 2 t 2 u u x u x y 2 3 α 2 t 2 u u y u x x α β t 2 u u x y α ν t 2 u u x x x y 1 2 α t 2 u u t x 1 2 α t 2 u t u x + 1 4 α t x u u x x + 1 2 α t y u u x y 1 4 α t x u x 2 1 2 α t y u x u y + 1 4 α t u u x 1 2 β t x u x + 1 2 β t u + 1 2 x y u x 1 2 u y , T 5 x = 1 2 α t 2 u u t t 1 2 β 2 t x u y + α t u u t 1 2 β t x u t + 1 2 β ν t u x y 1 4 α x u u x + 3 4 ν x y u x x y + 1 2 β x y u y 1 3 α 2 t x u u x u x y + 2 3 α 2 t y u u x u y y + 2 3 α 2 t y u u x y u y 1 8 α ν t x u u x x x y + 3 4 α ν t y u u x x y y 1 4 α β t x u u x y + 1 2 α β t y u u y y 5 8 α ν t x u x u x x y 1 2 α ν t y u x u x y y 5 4 α β t x u x u y + 3 8 α ν t x u x x u x y + 1 4 α ν t y u x x u y y 1 8 α ν t x u x x x u y α ν t y u x x y u y 4 3 α 2 t 2 u t u x u y + 3 4 α t 2 ν u u t x x y + 1 2 α t 2 ν u t x u x y + 1 4 α t 2 ν u t y u x x 1 4 α t 2 ν u t x x u y 1 2 α t 2 ν u t x y u x 3 4 α t 2 ν u x x y u t + 1 2 α β t 2 u u t y 1 2 α β t 2 u t u y 2 3 α 2 t x u x 2 u y 4 3 α 2 t y u x u y 2 + α 2 t u u x u y + 1 2 α ν t y u x y 2 1 2 α β t y u y 2 1 4 α t x u u t x + 1 2 α t y u u t y + 3 2 α ν t u u x x y + α β t u u y 1 4 α t x u t u x 1 2 α t y u t u y 1 2 α ν t u x u x y 3 4 β ν t x u x x y + α x y u x u y + 2 3 α 2 t 2 u u t x u y + 2 3 α 2 t 2 u u t y u x + 1 4 u 2 α + 1 2 ν u x 1 2 α t 2 u t 2 1 4 ν x u x x 1 2 ν y u x y + 1 2 u t x y , T 5 y = 1 3 α 2 t x u u x u x x 2 3 α 2 t y u u x u x y 2 3 α 2 t y u u x x u y + 1 4 α β t x u u x x + 1 8 α ν t x u u x x x x 3 4 α ν t y u u x x x y 1 2 α β t y u u x y 1 4 α ν t x u x u x x x 1 4 α ν t y u x u x x y 1 2 α β t y u x u y + 1 4 α ν t y u x x u x y 1 4 α ν t y u x x x u y + 2 3 α 2 t 2 u u t x u x 1 4 α ν t 2 u x u t x x 1 2 α β t 2 u x u t + 1 4 α t 2 ν u u t x x x + 1 4 α t 2 ν u x x u t x 1 4 α t 2 ν u x x x u t + 1 2 α β u t 2 u t x 2 3 α 2 t y u x 2 u y 3 4 α β t x u x 2 + 1 8 α ν t x u x x 2 α t y u u t x + 3 4 α β t u u x + 3 8 α ν t u u x x x 1 8 α ν t u x u x x 1 4 β ν t x u x x x 1 4 y ν u x x 1 2 β y u + 1 2 β 2 t u 2 3 α 2 t 2 u t u x 2 + 1 2 α x y u x 2 1 2 α y u u x + 1 2 y β x u x + 1 4 y ν u x x x x 1 3 α 2 t x u x 3 + 1 3 α 2 t u u x 2 1 2 β 2 t x u x + 1 4 β ν t u x x ;
Case 6. For Q 6 = 2 β t u y + 2 y u y + 4 t u t + x u x + u , we have the conserved vector
T 6 t = 8 3 u u x u x y α t + 4 3 u u x x u y α t + 2 u u x x x y ν t + 3 2 u u x y β t + u u t x t 1 4 u u x x x 1 2 u u x y y 1 4 u u x + 1 2 u x u y β t + u x u t t + 1 4 u x 2 x + 1 2 u x u y y , T 6 x = 1 2 β y u y 2 + 1 4 u t u x x + 1 2 u t u y y + 1 2 β 2 t u y 2 u u t t t 1 2 ν y u x y 2 3 2 ν u u x x y + 1 2 ν u x u x y + 1 4 x u u t x 1 2 y u u t y β u u y + u t 2 t u u t + 4 3 α β t u x u y 2 4 3 α t u u t x u y 4 3 α t u u t y u x 3 4 β ν t u u x x y y + 8 3 α t u t u x u y + 1 2 β ν t u x u x y y 1 4 β ν t u x x u y y + β ν t u x x y u y + 1 3 α x u u x u x y 2 3 α y u u x u y y 2 3 α y u u x y u y 2 3 α β t u u x u y y 2 3 α β t u u x y u y + 2 3 α x u x 2 u y + 4 3 α y u x u y 2 + 1 4 β x u x u y 1 2 β 2 t u u y y 1 2 β ν t u x y 2 3 2 ν t u u t x x y 3 2 β t u u t y + 3 2 ν t u t u x x y + 3 2 β t u t u y ν t u t x u x y + 1 2 ν t u t x x u y + ν t u t x y u x 1 2 ν t u t y u x x + 1 8 ν x u u x x x y 3 4 ν y u u x x y y + 5 8 ν x u x u x x y + 1 2 ν y u x u x y y 3 8 ν x u x x u x y 1 4 ν y u x x u y y + 1 8 ν x u x x x u y + ν y u x x y u y α u u x u y + 1 4 β x u u x y 1 2 β y u u y y , T 6 y = 2 3 α y u u y u x x + 2 3 α β t u x 2 u y 4 3 α t u u t x u x 1 3 α x u u x u x x + 2 3 α y u u x u x y + 1 4 β ν t u x u x x y + 3 4 β ν t u u x x x y 1 4 β ν t u x x u x y + 1 4 β ν t u x x x u y + 2 3 α β t u u y u x x + 2 3 α β t u u x u x y + 1 4 ν y u x x x u y + 1 2 u β 2 t u x y 1 4 u β u x x x + 1 2 u β u x y y + 4 3 α t u t u x 2 + 2 3 α y u x 2 u y + 1 2 β 2 t u x u y + β t u t u x + 1 2 ν t u t x x u x + 1 4 ν x u x u x x x + 1 4 ν y u x u x x y + 1 2 β y u x u y 1 2 ν t u u t x x x 1 8 ν x u u x x x x + 3 4 ν y u u x x x y + 1 2 ν t u t u x x x 1 2 ν t u t x u x x 1 4 ν y u x x u x y + u y u t x 1 3 α u u x 2 + 1 4 β x u x 2 1 4 β u u x + 1 8 ν u x u x x 1 8 ν x u x x 2 3 8 ν u u x x x + 1 3 α x u x 3 ;
Case 7. For Q 7 = y F ( t ) α u x F ( t ) , we have the conserved vector
T 7 t = 1 2 y u x F ( t ) + 1 4 α u u x x 1 4 α u x 2 F ( t ) , T 7 x = 1 2 y u F ( t ) + α y u x u y 1 4 α u u x + 3 4 ν y u x x y + 1 2 β y u y 1 2 u β + 1 2 u t y 1 4 ν u x x F ( t ) 1 3 α 2 u u x u x y + 2 3 α 2 u x 2 u y + 5 8 α ν u x u x x y + 1 4 α β u x u y + 1 4 α u t u x + 1 4 α β u u x y + 1 8 α ν u u x x x y + 1 4 α u u t x 3 8 α ν u x x u x y + 1 8 α ν u x x x u y F ( t ) T 7 y = 1 2 α y u x 2 + 1 2 β y u x + 1 4 ν y u x x x x F ( t ) + 1 3 α 3 u u x u x x 1 3 α 2 u x 3 1 4 α β u x 2 + 1 8 α ν u x x 2 1 4 α ν u x u x x x + 1 8 α ν u u x x x x + 1 4 α β u u x x F ( t ) ;
Case 8. For Q 8 = G ( t ) , we have the conserved vector
T 8 t = 1 2 u x G ( t ) , T 8 x = 1 2 u G ( t ) + α u x u y + 1 2 β u y + 3 4 ν u x x y + 1 2 u t G ( t ) , T 8 y = 1 2 α u x 2 + 1 2 β u x + 1 2 ν u x x x G ( t ) .

3.2. Conservation Laws Utilizing Noether’s Theorem

We now construct conservation laws for Equation (1). We utilize the classical Noether’s theorem [34,43] to achieve this task.
Equation (1) has a second-order Lagrangian
L = 1 2 ν u x x u x y u t u x α u x 2 u y β u x u y ,
since δ L / δ u = 0 on (1). Here, δ / δ u is the Euler operator. Thus, we can determine the Noether symmetries of (1) by employing the Lagrangian (37) on the determining equation
: X [ 2 ] L + L D t ( τ ) + D x ( ξ ) + D y ( ϕ ) D t ( G t ) D x ( G x ) D y ( G y ) = 0 ,
with gauge functions G t , G x , and G y . Expanding (38) and solving the resultant linear PDEs yield the following Noether symmetries and gauge functions:
X 1 = t , G t = 0 , G x = 0 , G y = 0 , X 2 = y , G t = 0 , G x = 0 , G y = 0 , X 3 = g ( t ) u , G t = 0 , G x = 1 2 u g ( t ) , G y = 0 , X 4 = 2 α t y + x u , G t = 1 2 u , G x = 0 , G y = 1 2 β u , X 5 = α f ( t ) x + y f ( t ) u , G t = 0 , G x = 1 2 u y f ( t ) 1 2 β u f ( t ) , G y = 0 , X 6 = 4 α t t + α x x + 2 α y y ( β x + α u ) u , G t = 1 2 β u , G x = 0 , G y = 1 2 β 2 u .
Thus, the conserved vectors corresponding to the above obtained Noether point symmetries, using [41]
T k = L ξ k + η α u x j α ξ j L u x k α l = 1 k D x l L u x l x k α + l = k n ζ l α u x l x j α ξ j L u x k x l α
are given by, respectively,
T 1 t = 1 2 α u x 2 u y 1 2 β u x u y + 1 2 ν u x x u x y , T 1 x = α u t u x u y + 3 4 ν u t u x x y 1 4 ν u x x u t y 1 2 ν u t x u x y + 1 2 β u t u y + 1 2 u t 2 , T 1 y = 1 2 α u t u x 2 + 1 2 β u t u x + 1 4 ν u t u x x x 1 4 ν u x x u t x ;
T 2 t = 1 2 u x u y , T 2 x = 1 2 u t u y + α u x u y 2 + 3 4 ν u y u x x y 1 2 ν u x y 2 1 4 ν u x x u y y + 1 2 β u y 2 , T 2 y = 1 4 ν u x x u x y + 1 4 ν u x x x u y 1 2 u t u x ; T 3 t = 1 2 u x g ( t ) , T 3 x = α u x u y 3 4 ν u x x y 1 2 β u y 1 2 u t g ( t ) + 1 2 u g ( t ) , T 3 y = 1 2 α u x 2 1 2 β u x 1 4 ν u x x x g ( t ) ; T 4 t = 1 2 u 1 2 x u x + α t u x u y , T 4 x = 1 2 ν u x y 1 2 x u t α x u x u y 1 2 β x u y 1 2 ν x u x x y + α t u t u y + 2 α 2 t u x u y 2 + α ν t u y u x x y α ν t u x y 2 α ν t u x x u y y + α β t u y 2 , T 4 y = 1 2 β u + α ν t u y u x x x 1 2 ν x u x x x 1 2 β x u x 1 2 α x u x 2 α t u t u x + α ν t u x x u x y ; T 5 t = 1 2 α u x 2 f ( t ) 1 2 y u x f ( t ) , T 5 x = 1 2 α 2 u x 2 u y + 3 4 α ν u x u x x y 1 4 α ν u x x u x y f ( t ) + 1 2 y u t α y u x u y 3 4 ν y u x x y + 1 4 ν u x x 1 2 β y u y + 1 2 β u f ( t ) + 1 2 y u f ( t ) , T 5 y = 1 2 α 2 u x 3 + 1 2 α β u x 2 + 1 4 α ν u x x x u x 1 4 α ν u x x 2 f ( t ) 1 2 α y u x 2 + 1 2 β y u x + 1 4 ν y u x x x f ( t ) ; T 6 t = 1 2 x α u x 2 2 t α 2 u y u x 2 + 1 2 x β u x + 1 2 α u u x + y α u y u x 2 t α β u y u x + 2 t α ν u x y u x x 1 2 β u , T 6 x = 1 2 x u y u x 2 α 2 + 2 y u y 2 u x α 2 + u u y u x α 2 + 4 t u y u x u t α 2 + y β u y 2 α y ν u x y 2 α + 2 t u t 2 α + 1 2 β u u y α + x β u y u x α ν u x u x y α 3 4 ν u y u x x α 1 2 y ν u y y u x x α 1 4 x ν u x y u x x α + 3 4 ν u u x x y α + 3 2 y ν u y u x x y α + 3 4 x ν u x u x x y α + 1 2 u u t α + y u y u t α + 2 t β u y u t α + 3 t ν u x x y u t α t ν u x x u t y α 2 t ν u x y u t x α + 1 2 x β 2 u y 1 2 β ν u x y + 3 4 x β ν u x x y + 1 2 x β u t ,
T 6 y = 1 2 x α 2 u x 3 + x α β u x 2 + 1 2 α 2 u u x 2 + 2 t α 2 u t u x 2 + 1 2 x β 2 u x + 1 2 α β u u x 1 2 α ν u x x u x + 1 4 x α ν u x x x u x y α u t u x + 2 t α β u t u x 1 4 x α ν u x x 2 1 4 β ν u x x + 1 2 y α ν u x y u x x + 1 4 x β ν u x x x + 1 4 α ν u u x x x + 1 2 y α ν u y u x x x + t α ν u x x x u t t α ν u x x u t x 1 2 β 2 u .
It should be noted that we derived eight conservation laws of (1) by invoking the multiplier method whereas Noether’s theorem yielded six conservation laws. However, we see that only one conservation law is common when using both methods. The conservation laws obtained with the Lagrange multiplier method could be obtained by means of Noether’s theorem using contact transformations. Thus, we conclude that there are benefits and constraints of the different methods used when deriving the conservation laws of nonlinear partial differential equations.

4. Concluding Remarks

In this paper, we investigated the generalized two-dimensional generalized shallow water wave Equation (1). We determined the Lie point symmetries of this equation, performed symmetry reductions, and reduced it to an NLODE (8). This NLODE (8) was integrated directly and its solution was obtained in the form of an incomplete elliptic integral. Moreover, Kudrayshov’s method was employed to obtain the solution of the NLODE (8). These solutions were presented graphically. Finally, we used the multiplier method and Noether’s theorem to derive the conservation laws of (1). The multiplier method yielded eight multipliers, which gave eight local conservation laws for Equation (1), whereas Noether’s theorem provided us with six local conservation laws, which included energy and linear momentum conservation laws.

Author Contributions

Conceptualization, C.M.K. and K.P.; Methodology, C.M.K.; Software, K.P.; Validation, C.M.K.; Writing—original draft, K.P.; Writing—review & editing, C.M.K. 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

Not applicable.

Acknowledgments

The authors thank OD Adeyemo for fruitful discussions and North-West University, Mafikeng campus, for their continued support. Additionally, the authors thank the reviewers for their positive suggestions, which helped to improve the paper enormously.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, X.X.; Tian, B.; Qu, Q.X.; Yuan, Y.Q.; Zhao, X.H. Lie group analysis, solitons, self-adjointness and conservation laws of the modified Zakharov-Kuznetsov equation in an electron-positron-ion magnetoplasma. Chaos Solitons Fractals 2020, 134, 109709. [Google Scholar] [CrossRef]
  2. Gao, X.Y. Mathematical view with observational/experimental consideration on certain (2 + 1)-dimensional waves in the cosmic/laboratory dusty plasmas. Appl. Math. Lett. 2019, 91, 165–172. [Google Scholar] [CrossRef]
  3. Zhang, C.R.; Tian, B.; Qu, Q.X.; Liu, L.; Tian, H.Y. Vector bright solitons and their interactions of the couple Fokas–Lenells system in a birefringent optical fiber, Z. Angew. Math. Phys. 2020, 71, 1–19. [Google Scholar]
  4. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Water-wave symbolic computation for the Earth, Enceladus and Titan: The higher-order Boussinesq-Burgers system, auto-and non-auto-Bäcklund transformations. Appl. Math. Lett. 2020, 104, 106170. [Google Scholar] [CrossRef]
  5. Benoudina, N.; Zhang, Y.; Khalique, C.M. Lie symmetry analysis, optimal system, new solitary wave solutions and conservation laws of the Pavlov equation. Commun. Nonlinear Sci. Numer. Simul. 2021, 94, 105560. [Google Scholar] [CrossRef]
  6. Gandarias, M.L.; Duran, M.R.; Khalique, C.M. Conservation laws and travelling wave solutions for double dispersion equations in (1 + 1) and (2 + 1) dimensions. Symmetry 2020, 12, 950. [Google Scholar] [CrossRef]
  7. Shafiq, A.; Khalique, C.M. Lie group analysis of upper convected Maxwell fluid flow along stretching surface. Alex. Eng. J. 2020, 59, 2533–2541. [Google Scholar] [CrossRef]
  8. Rosa, M.; Chuliána, S.; Gandarias, M.L.; Traciná, R. Application of Lie point symmetries to the resolution of an interface problem in a generalized Fisher equation. Physica D 2020, 405, 132411. [Google Scholar] [CrossRef]
  9. Wazwaz, A.M. Exact soliton and kink solutions for new (3 + 1)-dimensional nonlinear modified equations of wave propagation. Open Eng. 2017, 7, 169–174. [Google Scholar] [CrossRef]
  10. Wazwaz, A.M.; El-Tantawy, S. A new (3 + 1)-dimensional generalized Kadomtsev-Petviashvili equation. Nonlinear Dyn. 2016, 84, 1107–1112. [Google Scholar] [CrossRef]
  11. Wazwaz, A.M. Multiple-soliton solutions for a (3 + 1)-dimensional generalized KP equation. Commun. Nonlinear. Sci. Numer. Simul. 2012, 17, 491–495. [Google Scholar] [CrossRef]
  12. Wazwaz, A.M. The tanh and sine-cosine method for compact and noncompact solutions of nonlinear Klein-Gordon equation. Appl. Math. Comput. 2005, 167, 1179–1195. [Google Scholar] [CrossRef]
  13. Yang, J.Y.; Ma, W.X.; Khalique, C.M. Determining lump solutions for a combined soliton equation in (2 + 1)-dimensions. Eur. Phys. J. Plus 2020, 135, 494. [Google Scholar] [CrossRef]
  14. Yildirim, Y.; Yasar, E. An extended Korteweg–de Vries equation: Multi-soliton solutions and conservation laws. Nonlinear Dyn. 2017, 90, 1571–1579. [Google Scholar] [CrossRef]
  15. Mhlanga, I.E.; Khalique, C.M. A study of a generalized Benney–Luke equation with time-dependent coefficients. Nonlinear Dyn. 2017, 90, 1535–1544. [Google Scholar] [CrossRef]
  16. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  17. Zhang, L.; Khalique, C.M. Classification and bifurcation of a class of second-order ODEs and its application to nonlinear PDEs. Discret. Contin. Dyn.-Syst. 2018, 11, 777–790. [Google Scholar] [CrossRef]
  18. 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]
  19. Kudryashov, N.A.; Loguinova, N.B. Extended simplest equation method for nonlinear differential equations. Appl. Math. Comput. 2008, 205, 396–402. [Google Scholar] [CrossRef]
  20. Kudryashov, N.A. One method for finding exact solutions of nonlinear differential equations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 2248–2253. [Google Scholar] [CrossRef] [Green Version]
  21. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  22. Gu, C.H. Soliton Theory and Its Application; Zhejiang Science and Technology Press: Hangzhou, China, 1990. [Google Scholar]
  23. Matveev, V.B.; Salle, M.A. Darboux Transformations and Solitons; Springer: New York, NY, USA, 1991. [Google Scholar]
  24. Wang, M.; Zhou, Y.; Li, Z. Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  25. Wang, M.; Li, X.; Zhang, J. The (G/G)—Expansion method and travelling wave solutions for linear evolution equations in mathematical physics. Phys. Lett. A 2005, 24, 1257–1268. [Google Scholar]
  26. Olver, P.J. Applications of Lie Groups to Differential Equations, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1993. [Google Scholar]
  27. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1994; Volume 1. [Google Scholar]
  28. Ibragimov, N.H. Elementary Lie Group Analysis and Ordinary Differential Equations; John Wiley & Sons: Chichester, NY, USA, 1999. [Google Scholar]
  29. Khalique, C.M.; Adem, K.R. Explicit solutions and conservation laws of a (2 + 1)-dimensional KP-Joseph-Egri equation with power law nonlinearity. J. Appl. Nonlinear Dyn. 2018, 7, 1–9. [Google Scholar] [CrossRef]
  30. Khalique, C.M.; Adeyemo, O.D.I. Simbanefayi, On optimal system, exact solutions and conservation laws of the modified equal-width equation. J. Appl. Math. Nonlinear Sci. 2018, 3, 409–418. [Google Scholar] [CrossRef] [Green Version]
  31. Khalique, C.M.; Moleleki, L.D. A (3 + 1)-dimensional generalized BKP-Boussinesq equation: Lie group approach. Results Phys. 2019, 13, 102239. [Google Scholar] [CrossRef]
  32. Khalique, C.M.; Plaatjie, K.; Simbanefayi, I. Exact solutions of equal-width equation and its conservation laws. Open Phys. 2019, 17, 505–511. [Google Scholar] [CrossRef]
  33. Motsepa, T.; Khalique, C.M.; Gandarias, M.L. Symmetry analysis and conservation laws of the Zoomeron equation. Symmetry 2017, 9, 27. [Google Scholar] [CrossRef] [Green Version]
  34. Noether, E. Invariante variationsprobleme, Nachr. v. d. Ges. d. Wiss. zu Göttingen. Math.-Phys. Klasse 1918, 2, 235–257. [Google Scholar]
  35. 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]
  36. Leveque, R.J. Numerical Methods for Conservation Laws, 2nd ed.; Birkhäuser-Verlag: Basel, Switzerland, 1992. [Google Scholar]
  37. Ibragimov, N.H. A new conservation theorem. J. Math. Anal. Appl. 2007, 333, 311–328. [Google Scholar] [CrossRef] [Green Version]
  38. Naz, R.; Mahomed, F.M.; Mason, D.P. Comparison of different approaches to conservation laws for some partial differential equations in fluid mechanics. Appl. Math. Comput. 2008, 205, 212–230. [Google Scholar] [CrossRef]
  39. Sjöberg, A. On double reductions from symmetries and conservation laws. Nonlinear Anal. Real World Appl. 2009, 10, 3472–3477. [Google Scholar] [CrossRef]
  40. Yasar, E.; Özer, T. On symmetries, conservations laws and similarity solutions of foam drainage equation, Internat. J. Non-Linear Mech. 2011, 46, 357–362. [Google Scholar] [CrossRef]
  41. Sarlet, W. Comment on ‘Conservation laws of higher order nonlinear PDEs and the variational conservation laws in the class with mixed derivatives’. J. Phys. A Math. Theor. 2010, 43, 458001. [Google Scholar] [CrossRef]
  42. Motsepa, T.; Abudiab, M.; Khalique, C.M. A Study of an extended generalized (2 + 1)-dimensional Jaulent-Miodek equation. Int. J. Nonlin. Sci. Num. 2018, 19, 391–395. [Google Scholar] [CrossRef]
  43. 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; Springer: New York, NY, USA, 2017; Volume 79, pp. 119–182. [Google Scholar]
  44. Khalique, C.M.; Abdallah, S.A. Coupled Burgers equations governing polydispersive sedimentation; a Lie symmetry approach. Results Phys. 2020, 16, 102967. [Google Scholar] [CrossRef]
  45. Bruzón, M.S.; Gandarias, M.L. Traveling wave solutions of the K(m, n) equation with generalized evolution. Math. Meth. Appl. Sci. 2018, 41, 5851–5857. [Google Scholar] [CrossRef]
  46. Xin, X.; Zhang, L.; Xia, Y.; Liu, H. Nonlocal symmetries and exact solutions of the (2 + 1)-dimensional generalized variable coefficient shallow water wave equation. Appl. Math. Lett. 2019, 94, 112–119. [Google Scholar] [CrossRef] [Green Version]
  47. Lan, Z.; Gao, Y.; Yang, J.; Su, C.; Zuo, D. Solitons, bäcklund transformation, lax pair, and infinitely many conservation law for a (2 + 1)-dimensional generalised variable-coefficient shallow water wave equation. Naturforsch 2016, 71, 69–79. [Google Scholar] [CrossRef]
  48. Yildirim, Y.; Yasar, E. A (2 + 1) dimensional breaking soliton equation: Solutions and conservation laws, Chaos, Solitons and Fractals. Chaos Solitons Fractals 2018, 107, 146–155. [Google Scholar] [CrossRef]
  49. Wazwaz, A.M. Integrable (2 + 1)-dimensional and (3 + 1)-dimensional breaking soliton equations. Phys. Scr. 2010, 81, 035005. [Google Scholar] [CrossRef]
  50. Wazwaz, A.M. A variety of completely integrable Calogero-Bogoyavlenskii-Schiff equations with time-dependent coefficients. Int. J. Numer. Method H 2020. [Google Scholar] [CrossRef]
  51. Wazwaz, A.M. Multiple-soliton solutions for the Calogero-Bogoyavlenskii-Schiff, Jimbo-Miwa and YTSF equations. Appl. Math. Comput. 2008, 203, 592–597. [Google Scholar] [CrossRef]
  52. Wazwaz, A.M. New solutions of distinct physical structures to high-dimensional nonlinear evolution equations. Appl. Math. Comput. 2008, 196, 363–370. [Google Scholar] [CrossRef]
  53. Salem, S.; Kassem, M.; Mabrouk, S.M. Similarity Solution of (2 + 1)-dimensional Calogero-Bogoyavlenskii-Schiff equation Lax pair. Am. J. Appl. Math. 2019, 7, 137–144. [Google Scholar] [CrossRef] [Green Version]
  54. Shakeel, M.; Mohyud-Din, S.T. Improved (G/G)—Expansion and extended tanh methods for (2 + 1)-dimensional Calogero-Bogoyavlenskii-Schiff equation. Alex. Eng. J. 2015, 54, 27–33. [Google Scholar] [CrossRef] [Green Version]
  55. Najafi, M.; Najafi, M.; Arbabi, S. New application of (G/G)–Expansion method for generalized (2 + 1)-dimensional nonlinear evolution equations. J. Eng. Math. 2013, 5, 595–603. [Google Scholar] [CrossRef] [Green Version]
  56. Darvishi, M.T.; Najafi, M.; Najafi, M. New application of EHTA for the generalized (2 + 1)-dimensional nonlinear evolution equations. Int. J. Math. Comput. Sci. 2010, 6, 132–138. [Google Scholar]
  57. Najafi, M.; Najafi, M.; Arbabi, S. New exact solutions for the generalized (2 + 1)-dimensional nonlinear evolution equations by tanh-coth method. Int. J. Modern Theor. Phys. 2013, 2, 79–85. [Google Scholar]
  58. Najafi, M.; Arbabi, S.; Najafi, M. New application of sine-cosine method for the generalized (2 + 1)-dimensional nonlinear evolution equations. Int. J. Adv. Math. Sci. 2013, 1, 45–49. [Google Scholar] [CrossRef] [Green Version]
  59. Aminakbari, N.; Gu, Y.; Yuan, W. Meromorphic exact solutions of the (2 + 1)-dimensional generalized Calogero-Bogoyavlenskii-Schiff equation. Open Math. J. 2020, 18, 1342–1351. [Google Scholar] [CrossRef]
  60. Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1972. [Google Scholar]
  61. Kudryashov, N.A. Analytical Theory of Nonlinear Differential Equations; Institute of Computer Investigations: Moskow, Russia, 2004. [Google Scholar]
  62. Motsepa, T.; Khalique, C.M. Cnoidal and snoidal waves solutions and conservation laws of a generalized (2 + 1)-dimensional KdV equation. In Proceedings of the 14th Regional Conference on Mathematical Physics, Islamabad, Pakistan, 9–14 November 2015; World Scientific Publishing Co. Pte Ltd.: Singapore, 2018; pp. 253–263. [Google Scholar]
Figure 1. The 3D and 2D solution profiles of (13).
Figure 1. The 3D and 2D solution profiles of (13).
Mathematics 09 01439 g001
Figure 2. The 3D and 2D solution profile of (19).
Figure 2. The 3D and 2D solution profile of (19).
Mathematics 09 01439 g002
Figure 3. The 3D and 2D solution profiles of (32).
Figure 3. The 3D and 2D solution profiles of (32).
Mathematics 09 01439 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khalique, C.M.; Plaatjie, K. Exact Solutions and Conserved Vectors of the Two-Dimensional Generalized Shallow Water Wave Equation. Mathematics 2021, 9, 1439. https://doi.org/10.3390/math9121439

AMA Style

Khalique CM, Plaatjie K. Exact Solutions and Conserved Vectors of the Two-Dimensional Generalized Shallow Water Wave Equation. Mathematics. 2021; 9(12):1439. https://doi.org/10.3390/math9121439

Chicago/Turabian Style

Khalique, Chaudry Masood, and Karabo Plaatjie. 2021. "Exact Solutions and Conserved Vectors of the Two-Dimensional Generalized Shallow Water Wave Equation" Mathematics 9, no. 12: 1439. https://doi.org/10.3390/math9121439

APA Style

Khalique, C. M., & Plaatjie, K. (2021). Exact Solutions and Conserved Vectors of the Two-Dimensional Generalized Shallow Water Wave Equation. Mathematics, 9(12), 1439. https://doi.org/10.3390/math9121439

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