Next Article in Journal
Charge Asymmetry in Top Quark Pair Production
Previous Article in Journal
The Response of Nanobeams with Temperature-Dependent Properties Using State-Space Method via Modified Couple Stress Theory
Previous Article in Special Issue
Conservation Laws and Travelling Wave Solutions for Double Dispersion Equations in (1+1) and (2+1) Dimensions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lie Symmetries and Low-Order Conservation Laws of a Family of Zakharov-Kuznetsov Equations in 2 + 1 Dimensions

Department of Mathematics, Universidad de Cádiz, Puerto Real, 11510 Cádiz, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(8), 1277; https://doi.org/10.3390/sym12081277
Submission received: 8 June 2020 / Revised: 10 July 2020 / Accepted: 14 July 2020 / Published: 2 August 2020
(This article belongs to the Special Issue Conservation Laws and Symmetries of Differential Equations)

Abstract

:
In this work, we study a generalised ( 2 + 1 ) equation of the Zakharov–Kuznetsov (ZK) ( m , n , k ) equation involving three arbitrary functions. From the point of view of the Lie symmetry theory, we have derived all Lie symmetries of this equation depending on the arbitrary functions. Line soliton solutions have also been obtained. Moreover, we study the low-order conservation laws by applying the multiplier method. This family of equations is rich in Lie symmetries and conservation laws. Finally, when the equation is expressed in potential form, it admits a variational structure in the case when two of the arbitrary functions are linear. In addition, the corresponding Hamiltonian formulation is presented.

1. Introduction

In the context of plasma physics, the Zakharov–Kuznetsov equation (ZK) arises to describe ion-sound waves propagating along the magnetic field [1]. Its formal derivation involves the Euler–Poisson system for uniformly magnetised plasmas [2]. The two-dimensional ZK equation has the form
u t + u x u + x Δ u = 0 , u = u ( t , x , y ) .
The ZK equation describes the behaviour of weakly nonlinear ion-acoustic waves in plasma comprising cold ions and hot isothermal electrons in a uniform magnetic field. It generalised the well-known Korteweg de Vries (KdV) equation, it is not completely integrable and has a Hamiltonian.
There are several papers in which different generalisations of Equation (1), in two and three dimensions, have been studied. We shall now proceed to show some well-known generalised ZK equations in ( 2 + 1 ) -dimensions:
  • The modified Zakharov–Kuznetsov equation [3,4,5]
    u t + a u n / 2 u x + u x x + u y y x = 0 ,
    where a 0 and n 1 are arbitrary constants.
  • The ZK ( m , n , k )
    u t + a ( u m ) x + b ( u n ) x x x + c ( u k ) y y x = 0 ,
    where a , b , c are arbitrary constants while m , n and k are positive integers [6,7,8].
  • ZK ( m , n , k ) equation with generalised evolution and time-dependent coefficients
    ( u l ) t + a ( t ) ( u m ) x + b ( t ) ( u n ) x x x + c ( t ) ( u k ) y y x = α ( t ) u l ,
    where a ( t ) , b ( t ) , c ( t ) , α ( t ) are time-dependent coefficients and m , n , k , l are integers [9,10].
In the present paper, we are interested in studying a generalised ZK equation involving arbitrary functions. The simplest generalization would be
u t + u x f ( u ) + x Δ g ( u ) = 0 .
In the previous Equation (5), x Δ g ( u ) represents the dispersion term, with the same dispersion effect for directions x and y. However, in order to get a more general classification, we consider the case in which the dispersion effect could be different in the two directions.
Therefore, we study the following ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov equation involving three arbitrary functions (gZK)
u t + f ( u ) x + g ( u ) x x x + h ( u ) y y x = 0 ,
where f ( u ) , g ( u ) and h ( u ) are nonzero arbitrary functions.
When functions f , g and h are p-power nonlinearities, it was proved that the equation admits new solitary pattern, solitary wave and singular solitary wave solutions. Moreover, the authors proved a theorem on the convergence of the homotopy analysis method to solve this equation [11].
Due to the physical background of the equation, the objectives of this work are twofold. First, we seek point symmetries of Equation (6). Symmetries of a partial differential equation (PDE) leave invariant the whole space of solutions of the equation and, in that way, symmetries can be used to obtain reductions and exact group-invariant solutions. The Lie method determines all the Lie symmetries that a given PDE admits [12,13,14,15]. Moreover, for different cases of the ZK Equation (2)–(4), the Lie method has been shown as a useful tool to get exact solutions, including soliton solutions, cnoidal waves and travelling wave solutions [16,17,18,19]. In Section 3, line soliton solutions have been obtained.
Second, we aim to study the conservation laws of the gZK Equation (6) with physical interest. Conservation laws provide basic conserved quantities for all solutions, such as, mass, energy and so on. By using the multiplier method [20,21,22], we find all local conservation laws admitted by gZK Equation (6). Conservation laws for some special cases of the previous Equations (2)–(4) can be found in [23,24,25]. Consequently, the study of conservation laws of Equation (6) is also motivated to determine special cases for the arbitrary functions, f , g and h, with extra conservation laws.
Finally, making use of the fact that an equation admits a variational structure if and only if the Frechet derivative of the equation is self-adjoint (i.e., the Helmholtz conditions hold) [15,22], we present the case for the arbitrary functions f, g and h when the potential form of the gZK Equation (6) possesses a variational structure and we write the associated Hamiltonian formulation for the gZK Equation (6). Next, we determine the variational symmetries of the potential equation admitting the variational structure.

2. Lie Symmetries

We apply the Lie method to the two-dimensional gZK Equation (6). Therefore, we consider a one-parameter Lie group of point transformations acting on independent and dependent variables
t ˜ = t + ε τ ( t , x , y , u ) + O ( ε 2 ) , x ˜ = x + ε ξ 1 ( t , x , y , u ) + O ( ε 2 ) , y ˜ = y + ε ξ 2 ( t , x , y , u ) + O ( ε 2 ) , u ˜ = u + ε η ( t , x , y , u ) + O ( ε 2 ) ,
where ε is the group parameter. A Lie point symmetry for Equation (6) is a transformation (7) that leaves (6) invariant. From (7) one can obtain the associated vector field which is given by
X = τ ( t , x , y , u ) t + ξ 1 ( t , x , y , u ) x + ξ 2 ( t , x , y , u ) y + η ( t , x , y , u ) u .
A generator (8) is a point symmetry of gZK Equation (6) if
X ( 3 ) u t + f ( u ) x + g ( u ) x x x + h ( u ) y y x = 0
when Equation (6) holds. Here, X ( 3 ) is the third prolongation of the vector field (8) defined 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 ( ξ 1 ) u y D i 1 ( ξ 2 ) , η 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 ( ξ 1 ) u i 1 y D i 2 ( ξ 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 ( ξ 1 ) u i 1 i 2 y D i 3 ( ξ 2 ) ,
where D is 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 .
Invariance condition (9) can be expanded and split with respect to the derivatives of u obtaining an overdetermined linear system of equations for the infinitesimals τ ( t , x , y , u ) , ξ 1 ( t , x , y , u ) , ξ 2 ( t , x , y , u ) , η ( t , x , y , u ) along with the arbitrary functions f ( u ) , g ( u ) and h ( u ) . Thus we have the following result:
Theorem 1.
The classification of point symmetries admitted by the ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov Equation (6) is given by the following cases:
(i) 
For arbitrary f ( u ) , g ( u ) and h ( 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) 
Additional point symmetries are admitted by the ( 2 + 1 ) -dimensional generalised ZK Equation (6) in the following cases:
  • For arbitrary g ( u ) , h ( u ) , and f ( u ) = f 1 u + f 2 , there is an extra generator:
    X 4 = 3 t t + 2 f 1 t + x x + y y , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
  • For f ( u ) = f 1 ( u + a ) m + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ( u + a ) q + h 2 , the additional symmetry is
    X 5 = ( 3 m n 2 ) t t + 2 f 2 ( m 1 ) t + ( m n ) x x + ( m q ) y y 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    Moreover:
    For m = 1 and f 2 = 0 , generator X 4 is also admitted.
    For m = 2 , n = q = 1 , a = 0 , the following generator is admitted
    X 6 = 2 f 1 t x + u , G a l i l e a n   b o o s t
    For n = 1 , m = q = 1 3 and f 1 h 1 > 0 , two additional generators are admitted
    X 7 = 2 cos 2 f 1 h 1 y 1 y + 3 f 1 h 1 sin 2 f 1 h 1 y ( u + a ) u , X 8 = sin 2 f 1 h 1 y y 3 f 1 h 1 2 cos 2 f 1 h 1 y 1 ( u + a ) u .
    For n = 1 , m = q = 1 3 and f 1 h 1 < 0 , two additional generators are admitted
    X 9 = exp 2 f 1 h 1 y y 3 f 1 h 1 ( u + a ) u , X 10 = exp 2 f 1 h 1 y y + 3 f 1 h 1 ( u + a ) u .
  • For f ( u ) = f 1 ( u + a ) m + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 , and h ( u ) = h 1 ( u + a ) q + h 2 ,
    X 5 n = 0 = ( 3 m 2 ) t t + 2 f 2 ( m 1 ) t + m x x + ( m q ) y y 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    Moreover, for m = 1 and f 2 = 0 , X 4 is also admitted.
  • For f ( u ) = f 1 ( u + a ) m + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 5 q = 0 = ( 3 m n 2 ) t t + 2 f 2 ( m 1 ) t + ( m n ) x x + m y y 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    As before, for m = 1 and f 2 = 0 , X 4 is admitted.
  • For f ( u ) = f 1 ( u + a ) m + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 5 n = q = 0 = ( 3 m 2 ) t t + 2 f 2 ( m 1 ) t + m x x + m y y 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    Moreover, for m = 1 and f 2 = 0 , X 4 is admitted.
  • For f ( u ) = f 1 ln ( u + a ) + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ( u + a ) q + h 2
    X 5 m = 0 = ( n + 2 ) t t + 2 f 2 t + n x x + q y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
  • For f ( u ) = f 1 ln ( u + a ) + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 and h ( u ) = h 1 ( u + a ) q + h 2
    X 5 m = n = 0 = 2 t t + 2 f 2 t x + q y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
  • For f ( u ) = f 1 ln ( u + a ) + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 5 m = q = 0 = ( n + 2 ) t t + 2 f 2 t + n x x + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
  • For f ( u ) = f 1 ln ( u + a ) + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 5 m = n = q = 0 = t t + f 2 t x + ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 ( u + a ) ln ( u + a ) 1 + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ( u + a ) q + h 2
    X 11 = ( n 1 ) t t + 2 f 1 t + ( n 1 ) x x + ( q 1 ) y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 ( u + a ) ln ( u + a ) 1 + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 and h ( u ) = h 1 ( u + a ) q + h 2
    X 11 n = 0 = t t + 2 f 1 t x x + ( q 1 ) y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 ( u + a ) ln ( u + a ) 1 + f 2 u + f 3 , g ( u ) = g 1 ( u + a ) n + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 11 q = 0 = ( n 1 ) t t + 2 f 1 t + ( n 1 ) x x y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 ( u + a ) ln ( u + a ) 1 + f 2 u + f 3 , g ( u ) = g 1 ln ( u + a ) + g 2 and h ( u ) = h 1 ln ( u + a ) + h 2
    X 11 n = q = 0 = t t + 2 f 1 t x x y y + 2 ( u + a ) u , s h i f t   c o m b i n e d   w i t h   a s c a l i n g   a n d   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 e m u + f 2 u + f 3 , g ( u ) = g 1 e n u + g 2 and h ( u ) = h 1 e q u + h 2
    X 12 = ( 3 m n ) t t + 2 m f 2 t + ( m n ) x x + ( m q ) y y 2 u , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 e m u + f 2 u + f 3 , g ( u ) = g 1 u + g 2 and h ( u ) = h 1 e q u + h 2
    X 12 n = 0 = 3 m t t + m 2 f 2 t + x x + ( m q ) y y 2 u , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 e m u + f 2 u + f 3 , g ( u ) = g 1 e n u + g 2 and h ( u ) = h 1 u + h 2
    X 12 q = 0 = ( 3 m n ) t t + 2 m f 2 t + ( m n ) x x + m y y 2 u , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    For f ( u ) = f 1 e m u + f 2 u + f 3 , g ( u ) = g 1 u + g 2 and h ( u ) = h 1 u + h 2
    X 12 n = q = 0 = 3 m t t + m 2 f 2 t + x x + m y y 2 u , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    In the above, f 1 0 , f 2 , f 3 , g 1 0 , g 2 , h 1 0 , h 2 , a, m 0 , n 0 and q 0 are arbitrary constants.

    3. Line Soliton Solution

    A line soliton solution
    u ( t , x , y ) = w ( x + μ y λ t ) ,
    is a travelling wave solution of the form u ( t , x , y ) = w K · x λ t with x = ( x , y ) and K = ( 1 , μ ) . The solution depends on two parameters, where μ represents the direction of propagation of the line soliton, i.e., the inclination of the line soliton in the ( x , y ) -plane is μ = tan α , with α the angle from the positive y-axis in counterclockwise direction; whereas c = λ K = λ 1 + μ 2 represents the speed of the wave.
    In this section, we will determine line soliton solutions for the ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov Equation (6). Taking into account the line soliton formulation (10) into Equation (6) we obtain the nonlinear third-order ODE
    μ 2 h ( w ) + g ( w ) w ( z ) + 3 μ 2 h ( w ) + g ( w ) w ( z ) w ( z ) + μ 2 h ( w ) + g ( w ) ( w ( z ) ) 3 + f ( w ) λ w ( z ) = 0 ,
    where z = x + μ y λ t . Integrating Equation (11) with respect to z we obtain
    μ 2 h ( w ) + g ( w ) w ( z ) + μ 2 h ( w ) + g ( w ) ( w ( z ) ) 2 + f ( w ) λ w = 0 ,
    omitting the constant of integration. We recall that we are interested in a solitary wave solution, therefore we impose w, w , w 0 as z ± . On the other hand, if one considers f ( w ) = λ w , g ( w ) = μ 2 h ( w ) , Equation (12) is solved. Thus, it is obvious that
    u ( t , x , y ) = a sech m b ( x + μ y λ t ) ,
    with m, a and b nonzero arbitrary constants, is a smooth solution of Equation (6) for f ( u ) = λ u , g ( u ) = μ 2 h ( u ) satisfying the desired asymptotic conditions.
    In [9], the authors proved the existence of solitary wave solutions for Equation (6) with f ( u ) , g ( u ) and h ( u ) power functions. As far as we are concerned, line soliton solutions given by (13) have not been previously obtained.

    4. Conservation Laws

    We apply the direct method of the multipliers in order to obtain conservation laws of the ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov equation with three arbitrary functions (6).
    In two dimensions, a local conservation law is a divergence expression of the form
    D t T + D x X + D y Y = 0
    that holds for the whole set of solutions u ( t , x , y ) .
    T is called conserved density, and X , Y are spatial fluxes. The three are functions depending on t , x , y , u and spatial derivatives of u because, by using the equation, time derivatives can be vanished. A conserved current is an expression ( T , X , Y ) .
    Every non-trivial local conservation law (14) for Equation (6) can be expressed in its characteristic form
    D t T + D x X + D y Y = u t + f ( u ) x + g ( u ) x x x + h ( u ) y y x Q ,
    where
    Q = δ T δ u
    is the so called multiplier. Solving the following system of determining equations
    δ δ h u t + f ( u ) x + g ( u ) x x x + h ( u ) y y x Q = 0 ,
    that holds off the set of solutions of Equation (6), all multipliers Q can be found.
    All non-trivial conservation laws arise from multipliers (16). Moreover, as conservation laws of physical importance come from low-order multipliers, we have considered multipliers of the form Q ( t , x , y , u , u x , u y , u x x , u y y ) . Similarly to Lie symmetries, here the determining Equation (17) for low-order multipliers splits with respect to u t , u x x x , u x y y and their differential consequences, yielding an overdetermined system of equations for Q and the arbitrary functions which Equation (6) involves.
    Once the multipliers are found, the corresponding non-trivial local conservation laws can be obtained by integrating the characteristic Equation (15).
    Theorem 2.
    All low-order conservation laws of the ( 2 + 1 ) -dimensional gZK Equation (6) are
    (i) 
    For arbitrary f ( u ) , g ( u ) and h ( u ) , there is one multiplier Q 1 = F ( u ) of which the conservation law is
    T 1 = F y u , X 1 = F y g u u u x 2 + d 2 d y 2 F y h u + F y g u u x x + f u , Y 1 = F y h u u x y + F y h u u u x u y d d y F y h u u x .
    (ii) 
    Additional conservation laws are admitted in the following cases:
    • For arbitrary h ( u ) , f ( u ) = f 1 h u + f 2 u + f 3 , and g ( u ) = g 1 u + g 2 , there are seven extra multipliers.
      For the multiplier Q 2 = e f 1 y ,
      T 2 = e f 1 y u , X 2 = e f 1 y h u f 1 u y + f 1 h u + f 2 u + g 1 u x x , Y 2 = e f 1 y h u u u x u y + h u u x y .
      For the multiplier Q 3 = e f 1 y ,
      T 3 = e f 1 y u , X 3 = h u f 1 u y + f 1 h u + f 2 u + g 1 u x x e f 1 y , Y 3 = h u u u x u y + h u u x y e f 1 y .
      For the multiplier Q 4 = e f 1 y g 1 f 2 sin f 2 x g 1 ,
      T 4 = e f 1 y g 1 u f 2 sin f 2 x g 1 , X 4 = e f 1 y f 2 g 1 f 1 h u u y g 1 h u f 1 g 1 3 / 2 u x x sin f 2 x g 1 + cos f 2 x g 1 f 2 g 1 u x , Y 4 = e f 1 y f 2 g 1 h u u u x u y + h u u x y sin f 2 x g 1 + cos f 2 x g 1 h u f 2 f 1 .
      For the multiplier Q 5 = g 1 f 2 e f 1 y sin f 2 x g 1 ,
      T 5 = g 1 e f 1 y u f 2 sin f 2 x g 1 , X 5 = e f 1 y f 2 g 1 f 1 h u u y g 1 h u f 1 g 1 3 / 2 u x x sin f 2 x g 1 + cos f 2 x g 1 f 2 g 1 u x , Y 5 = e f 1 y f 2 g 1 h u u u x u y + h u u x y sin f 2 x g 1 + cos f 2 x g 1 h u f 2 f 1 .
      For the multiplier Q 6 = g 1 f 2 e f 1 y cos f 2 x g 1 ,
      T 6 = g 1 e f 1 y u f 2 cos f 2 x g 1 , X 6 = e f 1 y f 2 g 1 f 1 h u u y + g 1 h u f 1 + g 1 3 / 2 u x x cos f 2 x g 1 + g 1 sin f 2 x g 1 u x f 2 , Y 6 = e f 1 y f 2 g 1 h u u u x u y + h u u x y cos f 2 x g 1 + sin f 2 x g 1 f 1 h u f 2 .
      For the multiplier Q 7 = e f 1 y g 1 f 2 cos f 2 x g 1 ,
      T 7 = e f 1 y g 1 u f 2 cos f 2 x g 1 , X 7 = e f 1 y f 2 g 1 f 1 h u u y + g 1 h u f 1 + g 1 3 / 2 u x x cos f 2 x g 1 + g 1 sin f 2 x g 1 u x f 2 , Y 7 = e f 1 y f 2 g 1 h u u u x u y + h u u x y cos f 2 x g 1 + sin f 2 x g 1 f 1 h u f 2 .
      For the multiplier Q 8 = e f 1 y 2 + 1 t f 2 x e f 1 y f 2 ,
      T 8 = 1 f 2 u e f 1 y + e f 1 y t f 2 x , X 8 = 1 f 2 h u u y t f 2 x f 1 + f 1 t f 2 x h u + f 2 2 t u + g 1 t u x x x u f 2 g 1 x u x x u x e f 1 y + e f 1 y h u u y t f 2 x f 1 + f 1 t f 2 x h u + f 2 2 t u + g 1 t u x x x u f 2 g 1 x u x x u x , Y 8 = 1 f 2 f 1 h u + h u u u x u y + h u u x y t f 2 x e f 1 y + a 1 h u + h u u u x u y + h u u x y t f 2 x e f 1 y .
    • For arbitrary f ( u ) , h ( u ) , and g ( u ) = g 1 + g 2 h u , there is an additional multiplier Q 9 = h u , whose associated conservation law is
      T 9 = h u d u , X 9 = h u f u d u + h u h u u g 2 u x 2 + 1 2 g 2 u x 2 u y 2 h u 2 + h u h u g 2 u x x , Y 9 = h u h u u u x u y + h u u x y .
    • For arbitrary h ( u ) , f ( u ) = f 1 h u + f 2 u + f 3 , and g ( u ) = g 1 + g 2 h u , the multiplier is
      Q 10 = 4 g 2 k 1 e k 1 y 4 k 1 2 + f 1 e 1 / 4 g 2 f 1 f 2 t g 2 k 1 e g 2 k 1 x g 2 e 1 / 4 g 2 f 1 x g 2 k 1 e 1 / 4 y f 1 k 1 e g 2 k 1 f 2 t g 2 1 ,
      and its associated conservation law is
      T 10 = 4 g 2 k 1 u 4 k 1 2 + f 1 e 1 / 4 y 4 k 1 2 + f 1 g 2 + 4 k 1 2 + f 1 t f 2 x g 2 _ c 1 , X 10 = 16 16 _ c 1 3 + 4 f 1 k 1 k 1 h u u x k 1 2 + f 1 / 4 g 2 k 1 2 u x 2 g 2 h u u + k 1 u y k 1 2 k 1 g 2 u x x 1 / 4 u y f 1 h u + k 1 2 f 1 / 4 2 h u k 1 2 f 2 u ) g 2 e 1 / 4 y 4 k 1 2 + f 1 g 2 + 4 k 1 2 + f 1 t f 2 x g 2 k 1 ) , Y 10 = 4 4 k 1 3 + f 1 k 1 k 1 2 h u u u x u y + h u u x y g 2 + h u k 1 4 1 / 16 f 1 2 e 1 / 4 y 4 k 1 2 + f 1 g 2 + 4 k 1 2 + f 1 t f 2 x g 2 k 1 ) .
    • For f ( u ) = f 3 + f 2 u + f 1 e h 2 u , g ( u ) = g 2 + g 1 e h 2 u , and h ( u ) = h 1 e h 2 u + h 3 , there are three multipliers.
      For the multiplier Q 11 = e h 2 u ,
      T 11 = e h 2 u h 2 , X 11 = 1 2 e h 2 u b 1 u x 2 h 1 u y 2 h 2 2 + 2 g 1 h 2 u x x + f 1 h 2 e u h 2 + 2 f 2 h 2 , Y 11 = e 2 h 2 u h 1 h 2 h 2 u x u y + u x y .
      For the multiplier Q 12 = t f 1 h 2 e h 2 u + t f 2 x f 1 h 2 ,
      T 12 = t e h 2 u f 1 + u t f 2 x f 1 h 2 , X 12 = 1 2 f 1 h 2 t h 2 f 1 g 1 u x 2 + h 1 u y 2 h 2 2 2 g 1 h 2 u x x f 1 e 2 h 2 u + 2 g 1 u x 2 t f 2 x h 2 2 + 2 g 1 f 2 t u x x x u x x + u x h 2 + 4 f 1 t f 2 x 2 e h 2 u + 2 f 2 u t f 2 x , Y 12 = h 1 f 1 e h 2 u t f 1 c 2 e h 2 u + t f 2 x h 2 u x u y + u x y .
      For the multiplier
      Q 13 = 4 g 1 h 1 k 1 e k 1 y g 1 h 1 4 k 1 2 h 1 + f 1 e k 1 g 1 h 1 x g 1 e 1 / 4 g 1 h 1 f 1 f 2 t g 1 h 1 k 1 e 1 / 4 g 1 h 1 f 1 x g 1 h 1 k 1 e 1 / 4 y f 1 h 1 k 1 e k 1 g 1 h 1 f 2 t g 1 1 , T 13 = 4 g 1 h 1 k 1 u g 1 h 1 4 k 1 2 h 1 + f 1 e 1 / 4 4 k 1 2 h 1 + f 1 t f 2 x g 1 h 1 y g 1 4 k 1 2 h 1 + f 1 g 1 h 1 k 1 , X 13 = 16 1 g 1 h 1 16 h 1 k 1 3 + 4 f 1 k 1 g 1 k 1 k 1 2 h 1 + f 1 / 4 e h 2 u u x h 2 g 1 h 1 + h 1 2 k 1 4 h 1 2 h 2 k 1 3 u y + h 1 g 1 h 2 u x x + g 1 h 2 2 u x 2 + f 1 / 2 k 1 2 + 1 / 4 f 1 h 1 h 2 k 1 u y 1 / 16 f 1 2 e h 2 u + f 2 h 1 k 1 2 u e 1 / 4 4 k 1 2 h 1 + f 1 t f 2 x g 1 h 1 y g 1 4 k 1 2 h 1 + f 1 g 1 h 1 k 1 , Y 13 = 16 h 1 2 k 1 4 1 / 16 f 1 2 g 1 h 1 + k 1 2 g 1 h 1 2 h 2 h 2 u x u y + u x y g 1 h 1 16 h 1 k 1 3 + 4 f 1 k 1 e 1 / 4 4 t f 2 x k 1 2 h 1 + f 1 / 4 g 1 h 1 + 4 g 1 k 1 2 y h 1 + k 1 h 1 h 2 u 1 / 4 y f 1 g 1 h 1 k 1 .
    • For arbitrary f ( u ) , g ( u ) = g 1 u + g 2 , and h ( u ) = h 1 u + h 2 , there are two additional multipliers.
      For the multiplier Q 14 = u ,
      T 14 = 1 / 2 u 2 , X 14 = g 1 u u x x 1 / 2 h 1 u y 2 1 / 2 g 1 u x 2 + u d d u f u d u , Y 14 = h 1 u u x y .
      For the multiplier Q 15 = u y y + g 1 u x x h 1 + f u h 1 ,
      T 15 = 1 / 2 g 1 u x 2 h 1 u y 2 + 2 f u d u h 1 , X 15 = 1 / 2 g 1 2 u x x 2 + 2 g 1 h 1 u x x u y , y + h 1 2 u y y 2 + 2 g 1 f u u x x + 2 h 1 f u u y y + 2 g 1 u x + f u 2 h 1 , Y 15 = u x u y .
    • For f ( u ) = 1 / 2 f 1 u 2 + f 2 u + a 3 , g ( u ) = g 1 u + g 2 , and h ( u ) = h 1 u + h 2 , there are two extra multipliers.
      For the multiplier Q 16 = f 1 u f 2 t + x ,
      T 16 = 1 / 2 f 1 u + 2 f 2 t 2 x u , X 16 = 1 / 6 2 f 1 2 u 3 6 f 1 f 2 u 2 + 6 u x x g 1 6 h 1 u y y f 1 6 f 2 2 u + 3 g 1 u x 2 3 h 1 u y 2 f 1 6 f 2 g 1 u x x + h 1 u y y t + 1 / 2 f 1 x u 2 + f 2 x u + 1 / 6 6 x u x x 6 u x g 1 + h 1 x u y y , Y 16 = h 1 f 1 t u x 1 u y .
      For the multiplier Q 17 = u y y + g 1 u x x h 1 + 1 / 2 u 2 f 1 h 1 ,
      T 17 = 1 / 6 f 1 u 3 3 g 1 u x 2 3 h 1 u y 2 h 1 , X 17 = 1 24 h 1 12 g 1 2 u x x 2 + 24 h 1 u x x u y y + 12 f 1 u 2 u x x + 12 u x f 2 u x + 2 u t g 1 + 12 h 1 2 u y y 2 + 12 f 1 u 2 u y y 12 f 2 u y 2 h 1 + 3 f 1 2 u 4 + 4 f 1 f 2 u 3 , Y 17 = f 2 u x + u t u y .
    In the above, f 1 0 , f 2 , f 3 , g 1 0 , g 2 , h 1 0 , h 2 0 , h 3 0 and k 1 0 are arbitrary constants.

    5. Variational Structure and Hamiltonian Formulation

    The ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov Equation (6) can be written in potential form
    v t x + f ( v x ) x + g ( v x ) x x x + h ( v x ) y y x = 0 ,
    in terms of the potential
    u = v x .
    It is straightforward to prove that the potential Equation (35) admits a local Lagrangian structure
    δ L δ v = 0
    if and only if g ( v x ) = g 1 v x + g 2 and h ( v x ) = h 1 v x + h 2 , where δ L / δ v is the variational derivative with respect to the variable v
    δ δ v = v D t v t D x v x D y v y + D t 2 v t t + D x 2 v x x + D y 2 v y y + D t D x v t x + D t D y v t y + D x D y v x y + ,
    and where the Lagrangian is given by
    L = 1 2 g 1 v x x 2 + 1 2 h 1 v x y 2 1 2 v x v y f ( v x ) d v x
    in terms of the variable v.
    For g ( u ) = g 1 u + g 2 and h ( u ) = h 1 u + h 2 , the variational structure yields a Hamiltonian formulation for the gZK Equation (6), given by
    u t = D x δ H δ u
    where D x is the Hamiltonian operator and
    H = R 2 f ( u ) d u 1 2 g 1 u x 2 1 2 h 1 u y 2 d x d y
    is the Hamiltonian density.

    Variational Symmetries of the gZK Potential Equation

    The ( 2 + 1 ) -dimensional gZK potential Equation (35) with g ( v x ) = g 1 v x + g 2 and h ( v x ) = h 1 v x + h 2 , given by
    v t x + f ( v x ) x + g 1 v x x x x + h 1 v y y x x = 0 ,
    possesses a local Lagrangian structure (37) in terms of the Lagrangian (39), which we use to find the variational symmetries of Equation (42).
    A generator
    Y = τ ( t , x , y , v ) t + ξ 1 ( t , x , y , v ) x + ξ 2 ( t , x , y , v ) y + η ( t , x , y , v ) v ,
    is a Lie point symmetry of the gZK potential Equation (42) if
    Y ( 4 ) v t x + f ( v x ) x + g 1 v x x x x + h 1 v y y x x = 0
    when Equation (42) holds. Here Y ( 4 ) is the fourth prolongation of (43).
    The point symmetry (43) has an equivalent characteristic form given by
    Y ^ = P v , P = η τ v t ξ 1 v x ξ 2 v y
    in terms of the symmetry characteristic P.
    A symmetry (45) will be a variational symmetry if and only if it leaves invariant the Lagrangian density L up to a total divergence,
    Y ^ L = D t Ψ t + D x Ψ x + D y Ψ y ,
    with Ψ t , Ψ x , Ψ y depending on t, x, y, v, and derivatives of v. The invariance condition is equivalent to the equation:
    δ δ v P δ L δ v = 0 ,
    where δ δ v is the variational derivative (38), in terms of the Lagrangian functional L and the characteristic P [15,22]. In order to determine the variational symmetries of the ZK potential Equation (42), we verify which of the symmetries (43) satisfy the variational symmetry condition (47). We are interested in the cases when f ( v x ) 0 , since the Equation (42) is then nonlinear. Thus, we obtain the following results.
    Theorem 3.
    The classification of point symmetries admitted by the ( 2 + 1 ) -dimensional generalised Zakharov–Kuznetsov potential Equation (42) is given by the following cases:
    (i) 
    For arbitrary f ( v x ) , g 1 and h 1 , the admitted point symmetries are generated by:
    Y 1 = t , t i m e t r a n s l a t i o n . Y 2 = x , s p a c e t r a n s l a t i o n . Y 3 = y , s p a c e t r a n s l a t i o n . Y 4 = F ( t , y ) v .
    (ii) 
    Additional point symmetries are admitted by the (2+1)-dimensional generalised ZK potential Equation (6) in the following cases:
    • For f ( v x ) = f 1 ( v x + a ) m + f 2 v x + f 3 ,
      Y 5 = ( m 1 ) 3 t t + ( m 1 ) 2 f 2 t + x x + ( m 1 ) y y + ( ( m 3 ) v 2 a x ) v , s c a l i n g   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    • For f ( v x ) = 1 6 f 1 v x 3 + 1 2 f 2 v x 2 + f 3 v x + f 4 ,
      Y 6 = 3 f 1 t t f 1 ( 2 f 3 t + x ) t f 2 2 x f 1 y y + f 2 x v , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    • For f ( v x ) = 1 2 f 1 v x 2 + f 2 v x + f 3 ,
      Y 7 = f 1 t x + ( x f 2 t ) v , G a l i l e a n   b o o s t . Y 8 = 3 t t + 2 f 2 t + x x + y y v v , s c a l i n g   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    • For f ( v x ) = f 1 v x 3 + f 2 v x + f 3 ,
      Y 9 = 3 t t + 2 f 2 t + x x + y y , d i l a t i o n   c o m b i n e d   w i t h   a   G a l i l e a n   b o o s t .
    Theorem 4.
    The ( 2 + 1 ) -dimensional gZK potential Equation (42) admits the variational point symmetries spanned by Y 1 , Y 2 , Y 3 , Y 4 and Y 7 . Symmetry Y 5 is only variational for m = 7 3 and symmetry Y 6 is only variational for f 1 = 0 , which is equivalent to Y 7 . Finally, Y 8 and Y 9 are not variational symmetries.
    Noether’s theorem states that there is a one-to-one correspondence between conservation laws and variational symmetries for any equation admitting a variational structure [15]. The correspondence is equivalent to the condition:
    Q = P ,
    where Q is the conservation law multiplier and P is the variational symmetry characteristic. Therefore, for each of the variational symmetries of the gZK potential Equation (42) in Theorem 4, there is a corresponding conservation law associated to the multiplier Q = P .
    We focus our attention on the variational symmetries and conservation laws of the gZK potential Equation (42) that correspond to local symmetries and local conservation laws of the gZK Equation (6). We observe that the space-translation symmetry Y 2 of the potential Equation (42) correspond to the multiplier Q = v x of the potential Equation (42), that yields the multiplier Q 14 = u and conservation law (31) of gZK Equation (6). We also observe that the Galilean boost symmetry Y 7 of the potential Equation (42) corresponds to the multiplier Q = x f 2 t f 1 t v x of the potential Equation (42), that yields the multiplier Q 16 = x f 2 t f 1 t u and conservation law (33) of gZK Equation (6).

    6. Conclusions

    In this paper, we have considered a generalised ZK equation in ( 2 + 1 ) -dimensions depending on three arbitrary functions (6). We have performed a classification of point symmetries admitted by Equation (6) depending on the arbitrary functions f ( u ) , g ( u ) and h ( u ) . We have shown that if f ( u ) , g ( u ) and h ( u ) satisfy certain conditions, Equation (6) admits line soliton solutions. To our knowledge, these solutions have not appeared previously in the literature. We have constructed multipliers of the generalised ZK Equation (6) having dependence on dependent variables, independent variables, and derivatives of dependent variables up to second order. After computing multipliers, we have derived the corresponding low-order local conservation laws. Finally, we have presented the potential form of the gZK Equation (6) and obtained the Lagrangian structure in the case when g and h are linear functions. In addition, we have presented the associated Hamiltonian formulation. Moreover, we have determined the variational symmetries of the potential equation admitting a Lagrangian structure (42).

    Author Contributions

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

    Funding

    This research received no external funding.

    Acknowledgments

    We warmly thank the referees for their valuable comments and suggestions. The support of Junta de Andalucía group FQM-201 is gratefully acknowledged.

    Conflicts of Interest

    The authors declare no conflict of interest.

    Abbreviations

    The following abbreviations are used in this manuscript:
    ZKZakharov–Kuznetsov equation
    KdVKorteweg de Vries
    gZKGeneralised Zakharov–Kuznetsov equation with three arbitrary functions
    PDEPartial differential equation

    References

    1. Zakharov, V.E.; Kuznetsov, E.A. Three-dimensional solitons. Zhurnal Eksp. Teoret. Fiz 1974, 66, 594–597. [Google Scholar]
    2. Lannes, D.; Linares, F.; Saut, J.C. The Cauchy Problem for the Euler–Poisson System and Derivation of the Zakharov–Kuznetsov Equation. In Studies in Phase Space Analysis with Applications to PDEs; Cicognani, M., Colombini, F., Del Santo, D., Eds.; Springer: New York, NY, USA, 2013; pp. 181–213. [Google Scholar]
    3. Munro, S.; Parkes, E. The derivation of a modified Zakharov—Kuznetsov equation and the stability of its solutions. J. Plasma Phys. 1999, 62, 305–317. [Google Scholar] [CrossRef]
    4. Schamel, H. A modified Korteweg-de Vries equation for ion acoustic waves due to resonant electrons. J. Plasma Phys. 1973, 9, 377–387. [Google Scholar] [CrossRef]
    5. Wazwaz, A.M. The extended tanh method for the Zakharov–Kuznetsov (ZK) equation, the modified ZK equation, and its generalized forms. Commun. Nonlinear Sci. 2008, 13, 1039–1047. [Google Scholar] [CrossRef]
    6. Ma, H.C.; Yu, Y.D.; Ge, D.J. The auxiliary equation method for solving the Zakharov–Kuznetsov (ZK) equation. Comput. Math. Appl. 2009, 58, 2523–2527. [Google Scholar] [CrossRef] [Green Version]
    7. Hong-Cai, M.; Yao-Dong, Y.; Dong-Jie, G. New Exact Travelling Wave Solutions for Zakharov–Kuznetsov Equation. Commun. Theor. Phys. 2009, 51, 609–612. [Google Scholar] [CrossRef]
    8. Zayed, E.; Al-Nowehy, A. Exact solutions of the Biswas-Milovic equation, the ZK(m,n,k) equation and the K(m,n) equation using the generalized Kudryashov method. Open Phys. 2016, 14, 129–139. [Google Scholar] [CrossRef]
    9. Triki, H.; Wazwaz, A.M. A one-soliton solution of the ZK(m,n,k) equation with generalized evolution and time-dependent coefficients. Nonlinear Anal. Real World Appl. 2011, 12, 2822–2825. [Google Scholar] [CrossRef]
    10. Li, W.; Zhao, Y.M. Exact Solutions of the ZK(M,N,K) Equation with Generalized Evolution and Time-Dependent Coefficients. J. Inf. Comput. Sci. 2012, 7, 97–110. [Google Scholar]
    11. Fallahzadeh, A.; Fariborzi, M.A. A note on the convergence of the Zakharov-Kuznetsov equation by homotopy analysis method. J. Linear Topol. Algebra 2014, 3, 7–13. [Google Scholar]
    12. Bluman, G.W.; Cheviakov, A.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2009. [Google Scholar]
    13. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
    14. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations, Symmetries, Exact Solutions, and Conservation Laws; CRC Press: Boca Raton, FL, USA, 1993; Volume 1. [Google Scholar]
    15. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer-Verlag: New York, NY, USA, 1986. [Google Scholar]
    16. Johnpillai, A.; Kara, A.; Biswas, A. Symmetry Solutions and Reductions of a Class of Generalized (2+1)-dimensional Zakharov—Kuznetsov Equation. Int. J. Nonlinear Sci. Numer. Simul. 2012, 12, 45–50. [Google Scholar] [CrossRef]
    17. Khalique, C.M.; Adem, K.R. Exact solutions of the (2+1)-dimensional Zakharov–Kuznetsov modified equal width equation using Lie group analysis. Math. Comput. Model. 2011, 54, 184–189. [Google Scholar] [CrossRef]
    18. Moleleki, L.D.; Johnpillai, A.G.; Khalique, C.M. Symmetry Reductions and Exact Solutions of a Variable Coefficient (2+1)-Zakharov-Kuznetsov Equation. Math. Comput. Appl. 2012, 17, 132–139. [Google Scholar] [CrossRef] [Green Version]
    19. Li, C.; Zhang, J. Lie Symmetry Analysis and Exact Solutions of Generalized Fractional Zakharov-Kuznetsov Equations. Symmetry 2019, 11, 601. [Google Scholar] [CrossRef] [Green Version]
    20. 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]
    21. 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]
    22. 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; Melnik, R., Makarov, R., Belair, J., Eds.; Springer: New York, NY, USA, 2016; pp. 119–182. [Google Scholar]
    23. Mothibi, D.M.; Khalique, C.M. Conservation Laws and Exact Solutions of a Generalized Zakharov–Kuznetsov Equation. Symmetry 2015, 7, 949–961. [Google Scholar] [CrossRef]
    24. Tracinà, R. On the nonlinear self-adjointness of the Zakharov–Kuznetsov equation. Commun. Nonlinear Sci. 2014, 19, 377–382. [Google Scholar] [CrossRef]
    25. Kara, A.H. An Analysis of the Symmetries and Conservation Laws of the Class of Zakharov-Kuznetsov Equations. Math. Comput. Appl. 2010, 4, 658–664. [Google Scholar] [CrossRef] [Green Version]

    Share and Cite

    MDPI and ACS Style

    Bruzón, M.S.; Garrido, T.M.; Recio, E.; de la Rosa, R. Lie Symmetries and Low-Order Conservation Laws of a Family of Zakharov-Kuznetsov Equations in 2 + 1 Dimensions. Symmetry 2020, 12, 1277. https://doi.org/10.3390/sym12081277

    AMA Style

    Bruzón MS, Garrido TM, Recio E, de la Rosa R. Lie Symmetries and Low-Order Conservation Laws of a Family of Zakharov-Kuznetsov Equations in 2 + 1 Dimensions. Symmetry. 2020; 12(8):1277. https://doi.org/10.3390/sym12081277

    Chicago/Turabian Style

    Bruzón, María S., Tamara M. Garrido, Elena Recio, and Rafael de la Rosa. 2020. "Lie Symmetries and Low-Order Conservation Laws of a Family of Zakharov-Kuznetsov Equations in 2 + 1 Dimensions" Symmetry 12, no. 8: 1277. https://doi.org/10.3390/sym12081277

    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