Next Article in Journal
Special Issue Editorial “Symmetric Distributions, Moments and Applications”
Next Article in Special Issue
A New Approach to Approximate Solutions of Single Time-Delayed Stochastic Integral Equations via Orthogonal Functions
Previous Article in Journal
Performance Analysis on the Small-Scale Reusable Launch Vehicle
Previous Article in Special Issue
Extended Approach to the Asymptotic Behavior and Symmetric Solutions of Advanced Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conservation Laws and Travelling Wave Solutions for a Negative-Order KdV-CBS Equation in 3+1 Dimensions

1
Department of Mathematics, University of Cadiz, 11510 Puerto Real, Spain
2
Department of Mathematics, Quaid-e-Azam Campus, University of the Punjab, Lahore 54000, Pakistan
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(9), 1861; https://doi.org/10.3390/sym14091861
Submission received: 18 July 2022 / Revised: 18 August 2022 / Accepted: 19 August 2022 / Published: 6 September 2022
(This article belongs to the Special Issue Differential/Difference Equations and Its Application)

Abstract

:
In this paper, we study a new negative-order KdV-CBS equation in ( 3 + 1 ) dimensions which is a combination of the Korteweg-de Vries (KdV) equation and Calogero–Bogoyavlenskii–Schiff (CBS) equation. Firstly, we determine the Lie point symmetries of the equation and conservation laws by using the multiplier method. The conservation laws will be used to obtain a triple reduction to a second order ordinary differential equation (ODE), which lead to line travelling waves and soliton solutions. Such solitons are obtained via the modified form of simple equation method and are displayed through three-dimensional plots at specific parameter values to lend physical meaning to nonlinear phenomena. It illustrates that these solutions might be extremely beneficial in understanding physical phenomena in a variety of applied mathematics areas.

1. Introduction

The analysis of higher-dimensional nonlinear systems, particularly integrable systems, has exploded in popularity in recent years. Solitary waves theory finds applications in a wide range of scientific domains, including telecommunications, transport phenomena, ocean waves, quantum mechanics, plasma physics, nonlinear fibre optics, and many more. In latest years, research into generating higher-dimensional integrable equations has gotten greater attention, with various integrable models established in the setting of (2+1) and (3+1)-dimensional equations [1,2]. The domain of integrable equations is crucial to investigate as it clarifies the true nature of nonlinearity in science disciplines and reveals its scientific nature.
A variety of chemical, biological and physical phenomena are modeled via nonlinear partial differential equations (PDEs) which contribute a key role in nonlinear science [3]. It provides an abundance of physical data and a deeper understanding of the problem’s physical characteristics, resulting in additional applications [4,5,6,7]. Many tracks have been established for physical issues in recent years in order to come up with exact solutions for nonlinear PDEs using modern computer technologies. There are a number of effective approaches, for instance, the modified extended tanh-function technique [8], the extended modified auxiliary equation mapping method [9,10], the Riccati–Bernoulli Sub-ODE technique [11], the auxiliary equation technique [12], the ( G / G ) -expansion technique [13], homotopy perturbation technique [14], the unified technique [15], the generalized unified technique [16], modified simple equation technique [17], the extended trial function technique [18], the homogeneous balance technique [19], the extended direct algebraic technique [20], the modified extended direct algebraic technique [21] and the collective variable technique [22].
Among the most powerful strategies for constructing nonlinear PDEs solutions is Lie group transformations theory. The investigation of the invariance property of a particular differential equation under a continuous group of transformations is the core principle behind this theory. The application of Lie symmetry transformation technique to PDEs, results in reductions and invariant solutions. Invariant solutions are frequently employed to investigate analytical features. Conservation laws are well acknowledged to participate actively in the solution of a PDE. Not all PDEs obeying conservation laws have physical interpretations, but they are crucial in investigating the integrability of PDEs.
In a wide range of scientific and technical domains, the classical Lie symmetry theory is frequently used [23,24]. This idea was originally suggested by a mathematician named Sophus Lie in the early nineteenth century [25,26]. It has significantly increased applications in nonlinear PDEs and has proven beneficial in various fields of differential equations [27]. The fundamental purpose underlying Lie group theory is to use the invariance requirement of nonlinear PDEs to achieve the similarity variable and reduction equation, and subsequently to derive similarity and solitary wave solutions [28]. These solutions explain a fundamental and significant physical phenomenon. The Chen–Lee–Liu equation, Sawada–Kotera equation, Boussinesq equation, and many more models have recently been studied using the Lie symmetry technique [29,30,31,32]. In 2017, Wazwaz investigated the (2+1)-dimensional modified KdV-Calogero–Bogoyavlenskii–Schiff equation to discover abundant solutions with various physical features [2].
We consider a new negative-order KdV-CBS model in ( 3 + 1 ) -dimensions given by
u x t + u x x x y + 4 u x u x y + 2 u x x u y + λ u x x + μ u x y + ν u x z = 0 ,
where λ , μ and ν are unspecified coefficients. Equation (1) is the combination of the Korteweg–De Vries (KdV) equation and Calogero–Bogoyavlenskii–Schiff (CBS) equation. It should go without saying that for ν = 0 and μ = 0 , Equation (1) will be simplified to the negative-order KdV equation. Although, for ν = 0 and λ = 0 , Equation (1) will be simplified to the negative-order CBS equation. In addition, the investigated model passes Painlevé test with no restrictions on the compatibility criteria or the variables used in the equation. The considered equation has been recently modeled by Wazwaz [33].
The aim of this paper is to determine the conservation laws, symmetries, line travelling waves and line soliton solutions admitted by the negative-order KDV-CBS equation in (3+1) dimensions. First, in Section 2 all Lie point symmetries are derived. In Section 3 the conservation laws arising from the low order multipliers are obtained. In Section 4 travelling wave reduction and first integrals are constructed. The line travelling waves u = U ( x + b y + c z a t ) , where a , b and c are arbitrary constants, and ξ is a travelling wave transformation are considered for the negative-order KDV-CBS equation in (3+1) dimensions. The resulting fourth-order nonlinear differential equation for U is directly reduced to a first-order differential equation by using the multi-reduction method [34] to the invariant under translation conservation laws. These differential equation is then integrated to get the explicit form of the line travelling wave solutions. For some particular cases of the parameters the general solution is obtained in terms of the Weierstrass P function and the Weierstrass ζ function.
This paper then concentrates in Section 5 on the use of the extended simple equation method (ESEM) to extract solitons from the negative-order KdV-CBS equation in three dimensions. The ESEM method is a powerful mathematical technique for finding exact solutions to nonlinear partial differential equations (NLPDEs). Soliton solutions in three dimensions are established. Section 6 concludes with a few closing remarks.

2. Lie Point Symmetries

A one-parameter Lie group of transformations on ( t , x , y , z , u ) introduced by a vector field leaves the invariant solution space of the equation specified by Equation (1).
Each admitted point symmetry may be exploited to minimize the number of independent variables in Equation (1). For example, we might convert PDEs to ODEs. There may be symmetries of these ODEs that allow us to reduce the order of the equation and whose solutions appear to be related to invariant solutions u ( t , x , y , z ) of Equation (1).
We consider a one-parameter Lie group of infinitesimals transformations acting on independent and dependent variables of the new ( 3 + 1 ) -dimensional negative-order KdV-CBS (nKdV-nCBS) Equation (1) given by
t ˜ = t + ε τ ( t , x , y , z , u ) + O ( ε 2 ) , x ˜ = x + ε ξ 1 ( t , x , y , z , u ) + O ( ε 2 ) , y ˜ = y + ε ξ 2 ( t , x , y , z , u ) + O ( ε 2 ) , z ˜ = y + ε ξ 3 ( t , x , y , z , u ) + O ( ε 2 ) , u ˜ = u + ε η ( t , x , y , z , u ) + O ( ε 2 ) ,
where ε is the group parameter and the accompanying vector field looks like
X = τ ( t , x , y , z , u ) t + ξ 1 ( t , x , y , z , u ) x + ξ 2 ( t , x , y , z , u ) y + ξ 3 ( t , x , y , z , u ) z + η ( t , x , y , z , u ) u .
The resulting transformation group will be a point symmetry if and only if
p r ( 4 ) X u x t + u x x x y + 4 u x u x y + 2 u x x u y + λ u x x + μ u x y + ν u x z | ε = 0 ,
where p r ( 4 ) X is the fourth prolongation of the vector field Equation (2), and ε denotes the solution space of Equation (1). This determining equation Equation (3) splits with respect to derivatives of u, yielding an over-determined linear system of 168 equations for the infinitesimals τ ( t , x , y , u ) , ξ ( t , x , y , u ) , ϕ ( t , x , y , u ) and η ( t , x , y , u ) called determining system. The following are the outcomes of solving the determining system:
τ = f 2 ( z ν t ) t 2 + f 3 ( z ν t ) + f 4 ( z ν t ) , ξ 1 = 1 2 f 1 ( z ν t ) x + 1 2 f 2 ( z ν t ) t x + 1 2 f 4 ( z ν t ) x + F 8 ( t , z ) , ξ 2 = f 1 ( z ν t ) y + f 2 ( z ν t ) t y + f 6 ( z ν t ) + f 7 ( z ν t ) t , ξ 3 = f 2 ( z ν t ) ν t 2 + f 4 ( z ν t ) ν t + f 5 ( z ν t ) , η = 1 2 f 2 ( z ν t ) t u 1 2 f 4 ( z ν t ) u + 1 2 f 1 ( z ν t ) u + 1 4 f 2 ( z ν t ) x y 1 4 f 2 ( z ν t ) μ x t 1 4 f 4 ( z ν t ) μ x + 1 4 f 1 ( z ν t ) μ x + 1 4 f 7 ( z ν t ) x 3 4 f 2 ( z ν t ) λ y t + 1 2 F 8 t y + 1 2 F 8 z ν y 1 4 f 4 ( z ν t ) λ y 1 4 f 1 ( z ν t ) λ y + F 11 ( t , z ) ,
where f i ( z ν t ) ( i = 1 , 2 , , 7 ) , F 8 ( t , z ) and F 11 ( t , z ) are arbitrary functions.
Theorem 1.
(i) The point symmetries admitted by the ( 3 + 1 ) -dimensional negative-order KdV-CBS Equation (1) [33] are generated by:
X 1 = t t + x x + ( 2 μ t x ) y + ν t z u u , X 2 = t t , X 3 = z , X 4 = 4 t y + x u , X 5 = y , X F = F ( t , z ) x + ( 1 2 F t y + 1 2 F z ν y ) u , X G = G ( t , z ) u .
(ii) These symmetries comprise a five-dimensional algebra
X 1 , X 2 , X 3 , X 4 , X 5 .
(iii) These symmetries comprise two infinite-dimensional generators
X F = F ( t , z ) x + ( 1 2 F t y + 1 2 F z ν y ) u , X G = G ( t , z ) u .
Their commutator is given by
[ X 1 , X 2 ] = 2 μ X 5 ν X 3 ,
[ X 1 , X 4 ] = 2 X 4 ,
[ X 1 , X 5 ] = X 5 ,
[ X 2 , X 4 ] = 4 X 5 .

3. Conservation Laws

A local conservation law of a scalar PDE
G ( t , x , y , z , u , u t , u x , u y , u z , ) = 0 ,
for u ( t , x , y , z ) is the equation of continuity
D t T + D x Φ x + D y Φ y + D z Φ z = 0 ,
keeping on the space E of solutions of the PDE, where T is represents the conserved density and Φ = ( Φ x , Φ y , Φ z ) is the spatial flux vector, which are functions of t , x , y , z , u , and derivatives of u. While, ( T , Φ ) indicates the conserved current.
Every non-trivial conservation law of the PDE G = 0 derives from a multiplier, and there establish one-to-one correspondence between non-trivial conserved currents ( T , Φ ) | E modulo trivial ones and non-zero multipliers Q | E , with Q G = D t T + D x Φ x + D y Φ y + + D z Φ z retaining as an identity. Here, Q is a function of t , x , y , z , u , with u derivatives, such that Q | E is non-singular. Several explicit approaches may be used to produce a conserved current ( T , Φ ) | E for each solution Q.
For the negative-order KdV-CBS equation in ( 3 + 1 ) -dimensions (1), conservation laws have the characteristic form
D t T + D x Φ x + D y Φ y + D z Φ z = ( u x t + u x x x y + 4 u x u x y + 2 u x x u y + λ u x x + μ u x y + ν u x z ) Q .
Taking into account the following form of low-order multipliers:
Q ( x , z , u , u t , u x , u y , u z , u x x , u x y , u x x x , u x x y ) ,
all the previous form low-order multipliers are found by requiring that the divergence condition
E u ( ( u x t + u x x x y + 4 u x u x y + 2 u x x u y + λ u x x + μ u x y + ν u x z ) Q ) = 0 ,
is satisfied, where E u represents the Euler operator with respect to u [35,36,37]. The divergence condition for the multipliers (12) splits with respect to the derivatives of u ( t , x , y , z ) leading to an overdetermined system of 476 equations for Q with λ 0 , μ 0 , ν 0 , which can be easily solved. This yield the following proposition
Proposition 1.
The low-order multipliers admitted by the ( 3 + 1 ) -dimensional KdV-CBS Equation (1), with λ 0 , μ 0 , ν 0 , are given by
Q 1 = F ( z ) , Q 2 = u t , Q 3 = u x , Q 4 = u y , Q 5 = u z , Q 6 = x 4 z u y ν , Q 7 = 3 u x 2 + u x x x .
These multipliers generate non-trivial low-order conservation laws, which are listed below:
Theorem 2.
The conservation laws for the ( 3 + 1 ) -dimensional KdV-CBS Equation (1), with λ 0 , μ 0 , ν 0 are given by:
T 1 = F ( z ) u x , X 1 = ν F ( z ) u + F ( z ) ( λ + 2 u y ) u x , Y 1 = F ( z ) ( μ u x + u x 2 + u x x x ) , Z 1 = ν F ( z ) u x .
T 2 = 1 2 ( λ + 2 u y ) u x 2 1 2 ( μ u y + ν u z ) u x + 1 2 ( u x y u x x ) , X 2 = u t 2 2 + u t 2 ( ( 2 λ + 4 u y ) u x + ν u y + μ u z + 2 u x x y ) u t 1 2 ( u x y u t x ) , Y 2 = 1 2 u t u x ( μ + 2 u x ) 1 2 ( u t x u x x ) , Z 2 = 1 2 ν u t u x .
T 3 = 1 2 u x 2 , X 3 = 1 2 u x ( λ + 2 u y ) u x + 2 u x x y , Y 3 = 1 2 u x x 2 + 1 2 μ u x 2 + u x 3 , Z 3 = 1 2 ν u x 2 .
T 4 = 1 2 u x u y , X 4 = 1 2 ( μ + 4 u x ) u y 2 + 1 2 ( 2 λ u x + ν u z + u t + 2 u x x y ) 1 2 u x y 2 , Y 4 = 1 2 ( λ u x + ν u z + u t ) u x , Z 4 = 1 2 ν u x u y .
T 5 = 1 2 u x u z , X 5 = 1 2 ν u z 2 + 1 2 ( ( 2 λ + 4 u y ) u x + μ u y + 2 u x x y + u t ) u z 1 2 u x z u x y , Y 5 = 1 2 ( μ + 2 u x ) u x u z 1 2 u x x u x z , Z 5 = 1 2 ( λ + 2 u y ) u x 2 1 2 ( u t + μ u y ) u x + 1 2 u x x u x y .
T 6 = 1 ν ( ν x 2 z u y ) u x , X 6 = ν x u z + ( 2 x u x 2 z u z ) u y + x λ u x + x u x x y λ u u x y + 2 z ν ( ( μ 4 u x ) u y 2 + ( 2 λ u x u t 2 u x x y u y + u x y 2 ) , Y 6 = ( μ x + x u x + 2 z u z ) + 2 z ν ( u t + λ u x ) u x , Z 6 = 2 z u x u y ν u .
T 7 = 1 2 u x x 2 + u x 3 , X 7 = 1 2 ( λ + 2 u y ) u x x 2 + 1 2 ( 2 μ u x y + 2 ν u x z 4 u x u x y + 2 u t x ) u x x + ( λ + 2 u y ) u x 3 , Y 7 = 3 u x 2 u x x x + 1 2 u x x x 2 1 2 μ u x x 2 + u x u x x 2 + μ u x 3 + 5 2 u x 4 , Z 7 = 1 2 ν u x x 2 + ν u x 3 .

4. Travelling Wave Reduction and First Integrals

It is well-known that the most popular application of symmetry reduction is the reduction to ODE.
A line travelling wave takes the following form:
u ( t , x , y , z ) = U ( ξ ) , ξ = x + b y + c z a t ,
where a , b and c are arbitrary constants, and ξ is a travelling wave transformation [38].
By plugging (22) into Equation (1) gives a nonlinear fourth-order ODE
b U + ( 6 b U + λ + b μ + c ν a ) U = 0 .
Since this ODE arises from a symmetry reduction under translations of equation (1), the corresponding conservations laws of the equation invariant under translations will similarly reduce to a first integral of the ODE [39,40,41]. Furthermore, all first integrals that arise from symmetry invariant conservation laws can be found directly, by using the symmetry, through the general multi-reduction method introduced in [34]. This reduction yields two first integrals. The resulting first integrals of ODE (23) are derived from the corresponding symmetry invariant multipliers
Q 1 = u ξ , Q 2 = 1 ,
and are given by
Ψ 1 = 2 b U U + b U 2 4 b U 3 + ( a b μ c ν λ ) U 2 = C 1 , Ψ 2 = b U + 3 b U 2 + ( λ + b μ + c ν a ) U = C 2 .
Eliminating U yields a second-order ODE
( U ) 2 + 2 ( U ) 3 + ( μ + c ν a λ b ) ( U ) 2 2 C 2 b U + 2 C 1 b = 0 .
Consequently we have a direct triple reduction from the fourth-order PDE (1) to a second-order ODE (23).
Setting U = V yields a first-order ODE
( V ) 2 + 2 V 3 + ( μ + c ν a λ b ) V 2 2 C 2 b V + 2 C 1 b = 0 .
For a = b μ + c ν + λ , the general solution can be given in terms of the Weierstrass P function, i.e., the general solution is
V ( ξ ) = P ( 1 2 ξ , g 2 , g 3 ) ,
with g 2 = 4 C 2 b and g 3 = 4 C 1 b . Consequently, integrating once with respect to ξ yields
U ( ξ ) = I 2 ζ ( I 2 2 ξ , 4 C 2 b , 4 C 1 b ) .
where P ( ξ ; g 2 , g 3 ) is the Weierstrass elliptic function, general solution of
P 2 = 4 P 3 g 2 P g 3
with g 2 , g 3 arbitrary constants (cf. ref. [42] ), and ζ ( ξ ; g 2 , g 3 ) is the Wierstrass zeta function defined by ζ ( ξ ) = ξ P ( s ) d s , (cf. ref. [43], p. 641).
Setting C 1 = C 2 = 0 , α = 2 β and β = 1 4 b ( b μ c ν + a λ ) in (26) an exact solution is
V ( ξ ) = α sech ( ( β ) ξ ) 2 .
Consequently an exact solution for (23) is
U ( ξ ) = 2 β tanh ( ( β ) ξ ) ,
yielding the line-kink solution
u ( t , x , y , z ) = 2 β tanh ( ( β ) ( x + b y + c z a t ) ) ,
where a , b and c are arbitrary constants, for the ( 3 + 1 ) -dimensional negative-order KdV-CBS Equation (1).

5. Extraction of Solitons from a Negative-Order KdV-CBS Equation

In this section, we find the optical solitary wave solutions for Equation (26) using the extended simple equation method (ESEM). Using this method, the initial solution takes the following form [44]:
U ( ξ ) = i = N N B i φ i ( ξ ) ,
where B i ( i = N , N + 1 , , 1 , 0 , 1 , , N 1 , N ) are the unknown coefficients to be found and N is the positive integer which can be calculated using balance principle on Equation (26). By applying balance principle, we get N = 1 and therefore, Equation (27) takes the form
U ( ξ ) = B 1 φ ( ξ ) + B 0 + B 1 φ ( ξ ) ,
which satisfies the ansatz equation given by
φ ( ξ ) = b 0 + b 1 φ ( ξ ) + b 2 φ 2 ( ξ ) ,
where b 0 , b 1 and b 2 are arbitrary constants.
If we take b 1 = 0 in Equation (29), then the ansatz Equation (29) converts to Riccati equation, and we have
φ ( ξ ) = b 0 b 2 b 2 t a n b 0 b 2 ( ξ + ξ 0 ) , b 0 b 2 > 0 ,
φ ( ξ ) = b 0 b 2 b 2 t a n h b 0 b 2 ξ m ln ( ξ 0 ) 2 , ξ 0 > 0 , b 0 b 2 < 0 , m = ± 1 .
If we take b 0 = 0 in Equation (29), then the ansatz Equation (29) becomes Bernoulli equation, and we have
φ ( ξ ) = b 1 e b 1 ξ b 1 ξ 0 b 2 e b 1 ξ , b 1 > 0 ,
φ ( ξ ) = b 1 e b 1 ξ b 1 ξ 0 + b 2 e b 1 ξ , b 1 < 0 .
If we take b 0 = b 1 = 0 in Equation (29), then the ansatz Equation (29) converts to sparable equation, and in this case we have
φ ( ξ ) = 1 b 2 ξ + ξ 0 , b 2 0 .
Thus, following is the general solution to ansatz Equation (29):
φ ( ξ ) = b 1 4 b 0 b 2 b 1 2 t a n 4 b 0 b 2 b 1 2 2 ( ξ + ξ 0 ) 2 b 2 , 4 b 0 b 2 > b 1 2 , b 2 > 0 ,
φ ( ξ ) = b 1 + 4 b 0 b 2 b 1 2 t a n 4 b 0 b 2 b 1 2 2 ( ξ + ξ 0 ) 2 b 2 , 4 b 0 b 2 > b 1 2 , b 2 < 0 ,
where τ 0 represents the constant of integration. Now, after putting Equation (28) along with Equation (29) into Equation (26) generates a system of algebraic equations and further doing some algebra calculations, we get the following cases:
  • Case 1. When b 1 = 0 ,
A 1 = a λ c v b μ 8 b b 2 , A 1 = 2 b 2 , b 0 = λ + c v + b μ a 16 b 2 b .
By using parameters defined in Equation (37) along with the corresponding ansatz equation yields the following solution
U 1 ( ξ ) = a + λ + c v + b μ + 8 A 0 tan b 0 b 2 ( ξ + ξ 0 ) b 0 b 2 b 16 b 2 tan b 0 b 2 ( ξ + ξ 0 ) 2 b 0 b 8 b tan b 0 b 2 ( ξ + ξ 0 ) b 0 b 2 ,
U 2 ( ξ ) = 2 b 0 b 2 + A 0 b 0 b 2 tanh m ln ( ξ 0 ) 2 b 0 b 2 ξ + 2 b 2 b 0 tanh m ln ( ξ 0 ) 2 b 0 b 2 ξ 2 b 0 b 2 tanh b 0 b 2 ξ + m ln ( ξ 0 ) 2 ,
where b 0 = λ + c v + b μ a 16 b 2 b , ξ = x + b y + c z a t and m = ± 1 (Figure 1).
  • Case 2. When b 0 = 0 ,
A 1 = 0 , A 1 = 2 b 2 , b 1 = a b μ λ c v b .
By using paramors defined in Equation (40) along with the corresponding ansatz equation yields the following solution
U 3 ( ξ ) = A 0 a b μ λ c v b ξ 0 A 0 b 2 A 1 a b μ λ c v b e a b μ λ c v b ξ a b μ λ c v b ξ 0 b 2 e a b μ λ c v b ξ , b 0 ;
where A 1 = 2 b 2 and ξ = x + b y + c z a t (Figure 2).
  • Case 3. When b 0 = b 1 = 0 ,
A 1 = a λ c v b μ 6 b 2 b , A 1 = 2 b 2 .
By using paramors defined in Equation (42) along with the corresponding ansatz equation yields the following solution
U 4 ( ξ ) = A 1 b 2 2 ξ 2 2 A 1 b 2 ξ ξ 0 + A 1 ξ 0 2 A 0 b 2 ξ + A 0 ξ 0 + A 1 ξ 0 b 2 ξ ,
where A 1 = a λ c v b μ 6 b 2 b , A 1 = 2 b 2 and ξ = x + b y + c z a t (Figure 3).
  • Case 4.
( i ) : A 1 = b b 1 2 + λ + c v + b μ a 2 b 2 b , A 1 = 0 , b 0 = b b 1 2 + λ + c v + b μ a 4 b 2 b ; ( ii ) : A 1 = 0 , A 1 = 2 b 2 , b 0 = b b 1 2 + λ + c v + b μ a 4 b 2 b , b , b 2 0 .
By using paramors defined in Equation (44) along with the corresponding anstaz equation yields the following solution
U 5 ( ξ ) = 2 A 1 b 2 A 0 b 1 + A 0 tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b b 1 + tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b , 4 b 0 b 2 > b 1 2 , b 2 > 0 ,
U 6 ( ξ ) = 2 A 1 b 2 + A 0 b 1 + A 0 tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b b 1 + tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b , 4 b 0 b 2 > b 1 2 , b 2 < 0 ,
and
U 7 ( ξ ) = 2 A 0 b 2 A 1 b 1 + A 1 tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b 2 b 2 , b 2 > 0 ,
U 8 ( ξ ) = 2 A 0 b 2 + A 1 b 1 + A 1 tan 1 2 λ + c v + b μ a b ( ξ + ξ 0 ) λ + c v + b μ a b 2 b 2 , b 2 < 0 ,
where A 1 = a b b 1 2 λ c v b μ 2 b 2 b , A 1 = 2 b 2 and ξ = x + b y + c z a t (Figure 4 and Figure 5).

6. Concluding Remarks

The present work has obtained several new results for the ( 3 + 1 ) negative-order KdV-CBS equation. Firstly, Lie point symmetries and low-order conservation laws have been derived. From the translation symmetries, travelling wave reductions are obtained. Moreover, from the invariant conservation laws under translations two first integrals are obtained for the travelling wave ODE. It happens that these two first integrals are functionally independent and a triple reduction has been derived, yielding a first-order ODE. For some particular cases of the parameters, the general solution is obtained in terms of the Weierstrass function and a line-kink solution is obtained for the ( 3 + 1 ) negative-order KdV-CBS Equation (1). These explicit invariant solutions can be expected to be important in understanding the asymptotics of more general solutions of (1).
Secondly, the extended simple equation approach has been used to obtain the exact solutions for the negative-order KdV-CBS equation, which is a powerful method to solve nonlinear PDEs. These solutions include solitary wave solutions, specifically singular and dark solitons. In addition, numerical solutions are represented in 3D plots as shown in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 to provide physical meaning of nonlinear phenomena. These solutions could be highly useful in comprehending physical phenomena in several disciplines of applied mathematics.

Author Contributions

All authors contributed equally to the study conception, All authors have read and agreed to the published version of the manuscript.

Funding

MG is supported by FQM-201 from Junta de Andalucia and Cádiz University Plan propio.

Data Availability Statement

All the data included within the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wazwaz, A.M. A new (2+ 1)-dimensional Korteweg-de Vries equation and its extension to a new (3+ 1)-dimensional Kadomtsev-Petviashvili equation. Phys. Scr. 2011, 84, 035010. [Google Scholar] [CrossRef]
  2. Wazwaz, A.M. Abundant solutions of various physical features for the (2+ 1)-dimensional modified KdV-Calogero–Bogoyavlenskii–Schiff equation. Nonlinear Dyn. 2017, 89, 1727–1732. [Google Scholar] [CrossRef]
  3. El-Kalaawy, O.H.; Ibrahim, R.S. Solitary Wave Solution of the Two-Dimensional Regularized Long-Wave and Davey-Stewartson Equations in Fluids and Plasmas. Appl. Math. 2012, 3, 833–843. [Google Scholar] [CrossRef]
  4. Al-Amr, M.O. New applications of reduced differential transform method. Alex. Eng. J. 2014, 53, 243–247. [Google Scholar] [CrossRef]
  5. Helal, M.A.; Seadawy, A.R. Variational method for the derivative nonlinear Schrödinger equation with computational applications. Phys. Scr. 2009, 80, 350–360. [Google Scholar] [CrossRef]
  6. Salman, F.; Raza, N.; Basendwah, G.A.; Jaradat, M.M. Optical solitons and qualitative analysis of nonlinear Schrodinger equation in the presence of self steepening and self frequency shift. Results Phys. 2022, 39, 105753. [Google Scholar] [CrossRef]
  7. Jhangeer, A.; Raza, N.; Rezazadeh, H.; Seadawy, A. Nonlinear self-adjointness, conserved quantities, bifurcation analysis and travelling wave solutions of a family of long-wave unstable lubrication model. Pramana 2020, 94, 1–9. [Google Scholar] [CrossRef]
  8. Taghizadeh, N.; Mirzazadeh, M. The modified extended tanh method with the Riccati equation for solving nonlinear partial differential equations. Math. Aeterna 2012, 2, 145–153. [Google Scholar]
  9. Cheemaa, N.; Chen, S.; Seadawy, A.R. Propagation of isolated waves of coupled nonlinear (2+ 1)-dimensional Maccari system in plasma physics. Results Phys. 2020, 17, 102987. [Google Scholar] [CrossRef]
  10. Lu, D.; Seadawy, A. R.; Iqbal, M. Mathematical methods via construction of traveling and solitary wave solutions of three coupled system of nonlinear partial differential equations and their applications. Results Phys. 2018, 11, 1161–1171. [Google Scholar] [CrossRef]
  11. Yang, X.F.; Deng, Z.C.; Wei, Y. A Riccati-Bernoulli sub-ODE method for nonlinear partial differential equations and its application. Adv. Differ. Equ. 2015, 2015, 1–17. [Google Scholar] [CrossRef] [Green Version]
  12. Jiong, S. Auxiliary equation method for solving nonlinear partial differential equations. Phys. Lett. A 2003, 309, 387–396. [Google Scholar]
  13. Zhang, H. New application of the (G′/G)-expansion method. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 3220–3225. [Google Scholar] [CrossRef]
  14. Liao, S. Comparison between the homotopy analysis method and homotopy perturbation method. Appl. Math. Comput. 2005, 169, 1186–1194. [Google Scholar] [CrossRef]
  15. Raza, N.; Rafiq, M.H.; Kaplan, M.; Kumar, S.; Chu, Y. The unified method for abundant soliton solutions of local time fractional nonlinear evolution equations. Results Phys. 2021, 22, 103979. [Google Scholar] [CrossRef]
  16. Wazwaz, A.M.; Osman, M.S. Analyzing the combined multi-waves polynomial solutions in a two-layer-liquid medium. Comput. Math. Appl. 2018, 76, 276–283. [Google Scholar] [CrossRef]
  17. Javid, A.; Raza, N. Singular and dark optical solitons to the well posed Lakshamanan-Porsezian-Daniel model. Optik 2018, 171, 120–129. [Google Scholar] [CrossRef]
  18. Ekici, M.; Sonmezoglu, A.; Zhou, Q.; Biswas, A.; Ullah, M.Z.; Asma, M.; Moshokoa, S.P.; Belic, M. Optical solitons in DWDM system by extended trial equation method. Optik 2017, 141, 157–167. [Google Scholar] [CrossRef]
  19. Wang, M.L.; Zhou, Y.B.; Li, Z.B. Application of homogeneous balance method to exact solutions of nonlinear equations in mathematical physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  20. Vahidi, J.; Zabihi, A.; Rezazadeh, H.; Ansari, R. New extended direct algebraic method for the resonant nonlinear Schrödinger equation with Kerr law nonlinearity. Optik 2021, 227, 165936. [Google Scholar] [CrossRef]
  21. Arshad, M.; Seadawy, A.R.; Lu, D. Optical soliton solutions of the generalized higher-order nonlinear Schrödinger equations and their applications. Opt. Quantum Electron. 2018, 5, 1–6. [Google Scholar] [CrossRef]
  22. Alhussain, Z.A.; Raza, N. New Optical Solitons with Variational Principle and Collective Variable Strategy for Cold Bosons in Zig-Zag Optical Lattices. J. Math. 2022, 2022. [Google Scholar] [CrossRef]
  23. Kumar, M.; Kumar, R.; Kumar, A. On similarity solutions of Zabolotskaya-Khokhlov equation. Comput. Math. Appl. 2014, 68, 454–463. [Google Scholar] [CrossRef]
  24. Ray, S.S. Lie symmetry analysis and reduction for exact solution of (2+ 1)-dimensional Bogoyavlensky-Konopelchenko equation by geometric approach. Mod. Phys. Lett. B 2018, 32, 1850127. [Google Scholar] [CrossRef]
  25. Vinita; Ray, S.S. Optical soliton group invariant solutions by optimal system of Lie subalgebra with conservation laws of the resonance nonlinear Schrödinger equation. Mod. Phys. Lett. B 2020, 34, 2050402. [Google Scholar] [CrossRef]
  26. Wazwaz, A.M. Exact solutions of compact and noncompact structures for the KP-BBM equation. Appl. Math. Comput. 2005, 169, 700–712. [Google Scholar] [CrossRef]
  27. Olver, P.J. Application of Lie Groups to Differential Equations; Springer: NewYork, NY, USA, 1993. [Google Scholar]
  28. Kumar, S.; Kumar, D.; Wazwaz, A.M. Group invariant solutions of (3+ 1)-dimensional generalized B-type Kadomstsev Petviashvili equation using optimal system of Lie subalgebra. Phys. Scr. 2019, 94, 065204. [Google Scholar] [CrossRef]
  29. Bai, Y.; Bilige, S.; Chaolu, T. Potential symmetries, one-dimensional optimal system and invariant solutions of the coupled Burgers? equations. J. Appl. Math. Phys. 2018, 6, 1825. [Google Scholar] [CrossRef]
  30. Kara, A.H.; Biswas, A.; Zhou, Q.; Moraru, L.; Moshokoa, S.P.; Belic, M. Conservation laws for optical solitons with Chen?Lee?Liu equation. Optik 2018, 174, 195–198. [Google Scholar] [CrossRef]
  31. Jafari, H.; Kadkhoda, N.; Azadi, M.; Yaghobi, M. Group classification of the time-fractional Kaup-Kupershmidt equation. Sci. Iran. 2017, 24, 302–307. [Google Scholar] [CrossRef]
  32. Kadkhoda, N.; Jafari, H.; Moremedi, G.M.; Baleanu, D.; Qaenat, I.R.; Babolsar, I.R. Optimal system and symmetry reduction of the (1+ 1) dimensional Sawada–Kotera equation. Int. J. Pure Appl. Math. 2016, 103, 215–226. [Google Scholar] [CrossRef]
  33. Wazwaz, A.M. Two new Painlevé integrable KdV-Calogero–Bogoyavlenskii–Schiff (KdV-CBS) equation and new negative-order KdV-CBS equation. Nonlinear Dyn. 2021, 104, 4311–4315. [Google Scholar] [CrossRef]
  34. Anco, S.C.; Gandarias, M. Symmetry multi-reduction method for partial differential equations with conservation laws. Commun. Nonlin. Sci. Numer. Simul. 2020, 91, 105349. [Google Scholar] [CrossRef]
  35. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: Berlin/Heidelberg, Germany, 1986. [Google Scholar]
  36. Bluman, G.W.; Cheviakov, A.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2009. [Google Scholar]
  37. 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; Springer: New York, NY, USA, 2017; pp. 119–182. [Google Scholar]
  38. Anco, S.C.; Gandarias, M.L.; Recio, E. Line-solitons, line-shocks, and conservation laws of a universal KP-like equation in 2+1 dimensions. J. Math. Anal. Appl. 2021, 504, 125319. [Google Scholar] [CrossRef]
  39. 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]
  40. Sjöberg, A. On double reduction from symmetries and conservation laws. Nonlin. Analysis Real World Appl. 2009, 10, 3472–3477. [Google Scholar] [CrossRef]
  41. Bokhari, A.H.; Dweik, A.Y.; Zaman, F.D.; Kara, A.H.; Mahomed, F.M. Generalization of the double reduction theory. Nonlin. Analysis Real World Appl. 2010, 11, 3763–3769. [Google Scholar] [CrossRef]
  42. Whittaker, E.E.; Watson, G.M. Modern Analysis, 4th ed.; Cambridge University Press: Cambridge, UK, 1927. [Google Scholar]
  43. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Tables; Dover: New York, NY, USA, 1965. [Google Scholar]
  44. Lu, D.; Seadawy, A.; Arshad, M. Applications of extended simple equation method on unstable nonlinear Schrödinger equations. Optik 2017, 140, 136–144. [Google Scholar] [CrossRef]
Figure 1. Graphical representation of Equations (38) and (39) with suitable parameters b = 0.2 , c = 0.1 , v = 0.4 , λ = 0.1 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , m = 1 , y = 0 , z = 0 : (a) a = 0.5 , b 2 = 0.2 (b) a = 0.5 , b 2 = 0.2 .
Figure 1. Graphical representation of Equations (38) and (39) with suitable parameters b = 0.2 , c = 0.1 , v = 0.4 , λ = 0.1 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , m = 1 , y = 0 , z = 0 : (a) a = 0.5 , b 2 = 0.2 (b) a = 0.5 , b 2 = 0.2 .
Symmetry 14 01861 g001
Figure 2. Graphical representation of Equation (41) with suitable parameters b = 0.1 , c = 0.3 , v = 1.2 , λ = 0.1 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , m = 1 , y = 0 , z = 0 , a = 0.4 , b 2 = 0.2 .
Figure 2. Graphical representation of Equation (41) with suitable parameters b = 0.1 , c = 0.3 , v = 1.2 , λ = 0.1 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , m = 1 , y = 0 , z = 0 , a = 0.4 , b 2 = 0.2 .
Symmetry 14 01861 g002
Figure 3. Graphical representation of Equation (43) with suitable parameters a = 0.5 , b = 0.1 , c = 0.1 , v = 0.2 , λ = 0.4 , μ = 1 , A 0 = 0.6 , ξ 0 = 1 , b 2 = 1 , m = 1 , y = 0 , z = .
Figure 3. Graphical representation of Equation (43) with suitable parameters a = 0.5 , b = 0.1 , c = 0.1 , v = 0.2 , λ = 0.4 , μ = 1 , A 0 = 0.6 , ξ 0 = 1 , b 2 = 1 , m = 1 , y = 0 , z = .
Symmetry 14 01861 g003
Figure 4. Graphical representation of Equations (45) and (46) with suitable parameters a = 1.5 , b = 0.2 , c = 0.2 , v = 0.3 , λ = 0.4 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , b 1 = 0.2 , m = 1 , y = 0 , z = 0 : (a) b 2 = 0.2 (b) b 2 = 0.2 .
Figure 4. Graphical representation of Equations (45) and (46) with suitable parameters a = 1.5 , b = 0.2 , c = 0.2 , v = 0.3 , λ = 0.4 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , b 1 = 0.2 , m = 1 , y = 0 , z = 0 : (a) b 2 = 0.2 (b) b 2 = 0.2 .
Symmetry 14 01861 g004
Figure 5. Graphical representation of Equations (47) and (48) with suitable parameters a = 1.5 , b = 0.2 , c = 0.2 , v = 0.3 , λ = 0.4 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , b 1 = 0.2 , m = 1 , y = 0 , z = 0 : (a) b 2 = 0.2 (b) b 2 = 0.2 .
Figure 5. Graphical representation of Equations (47) and (48) with suitable parameters a = 1.5 , b = 0.2 , c = 0.2 , v = 0.3 , λ = 0.4 , μ = 1 , A 0 = 0.1 , ξ 0 = 1 , b 1 = 0.2 , m = 1 , y = 0 , z = 0 : (a) b 2 = 0.2 (b) b 2 = 0.2 .
Symmetry 14 01861 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gandarias, M.L.; Raza, N. Conservation Laws and Travelling Wave Solutions for a Negative-Order KdV-CBS Equation in 3+1 Dimensions. Symmetry 2022, 14, 1861. https://doi.org/10.3390/sym14091861

AMA Style

Gandarias ML, Raza N. Conservation Laws and Travelling Wave Solutions for a Negative-Order KdV-CBS Equation in 3+1 Dimensions. Symmetry. 2022; 14(9):1861. https://doi.org/10.3390/sym14091861

Chicago/Turabian Style

Gandarias, Maria Luz, and Nauman Raza. 2022. "Conservation Laws and Travelling Wave Solutions for a Negative-Order KdV-CBS Equation in 3+1 Dimensions" Symmetry 14, no. 9: 1861. https://doi.org/10.3390/sym14091861

APA Style

Gandarias, M. L., & Raza, N. (2022). Conservation Laws and Travelling Wave Solutions for a Negative-Order KdV-CBS Equation in 3+1 Dimensions. Symmetry, 14(9), 1861. https://doi.org/10.3390/sym14091861

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