Next Article in Journal
The Existence of rG Family and tG Family, and Their Geometric Invariants
Next Article in Special Issue
The Application of Accurate Exponential Solution of a Differential Equation in Optimizing Stability Control of One Class of Chaotic System
Previous Article in Journal
A New Geometric Metric in the Shape and Size Space of Curves in R n
Previous Article in Special Issue
One-Stage Multiple Comparisons with the Control for Exponential Median Lifetimes under Heteroscedasticity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science

by
Chaudry Masood Khalique
1,2,3,* and
Oke Davies Adeyemo
1
1
International Institute for Symmetry Analysis and Mathematical Modelling, Department of Mathematical Sciences, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho 2735, South Africa
2
College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, Shandong, China
3
Department of Mathematics and Informatics, Azerbaijan University, Jeyhun Hajibeyli str., 71, Baku AZ1007, Azerbaijan
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(10), 1692; https://doi.org/10.3390/math8101692
Submission received: 21 August 2020 / Revised: 19 September 2020 / Accepted: 21 September 2020 / Published: 1 October 2020
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
In this article, we examine a (3+1)-dimensional generalized breaking soliton equation which is highly applicable in the fields of engineering and nonlinear sciences. Closed-form solutions in the form of Jacobi elliptic functions of the underlying equation are derived by the method of Lie symmetry reductions together with direct integration. Moreover, the ( G / G ) -expansion technique is engaged, which consequently guarantees closed-form solutions of the equation structured in the form of trigonometric and hyperbolic functions. In addition, we secure a power series analytical solution of the underlying equation. Finally, we construct local conserved vectors of the aforementioned equation by employing two approaches: the general multiplier method and Ibragimov’s theorem.

1. Introduction

It is an indisputable fact that nonlinear evolutionary equations (NLEEs) have provided a great deal of assistance in modeling several real-world problems consisting of diffusion, dispersion and convection that possess nonlinear effects. Achievements in solving these NLEEs through analytical approaches or numerical techniques guarantee a better understanding of the physical phenomena which form the background for the models [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23].
Solitons are defined as localized wave disturbances which are propagated with their shape neither changing nor even spreading out. Soliton notion was first established by Zabusky and Kruskal in the year 1965 [8]. Solitons play crucial roles in engineering and nonlinear sciences applications and furnish us with the sagacity of the pertinent phenomena in nonlinear science. Here, we give an example of a soliton-related equation.
The popularly celebrated Korteweg-de Vrie (KdV) equation ascribed to nonlinear dispersive partial differential equation (NLDPDE),
u t 6 u u x + u x x x = 0 ,
is a mathematical model that describes the special waves that are referred to as solitons on shallow water surfaces [8]. The KdV equation possesses numerous inter-connectivities with physical problems that comprise long embedded waves submerged in a density-partitioned ocean alongside waves in shallow water characterized by nonlinear weak reinstating forces, acoustic waves on a crystal lattice and plasma-accompanying ion-acoustic waves.
Besides the inverse scattering transform [9], which was developed to find solutions to the KdV Equation (1), many other essential approaches were introduced by researchers. Lately, many mathematicians and physicists have established effective and efficient techniques for constructing viable analytic solutions to NLEEs, namely, the Painlévé expansion technique [10], the mapping method and extended mapping technique [11,12], the Bäcklund transformation [13], the rational expansion method [14], the F -expansion technique [15], the tanh and sine-cosine method [16], the extended simplest equation method [17], Hirota’s technique [18], Lie symmetry analysis [19,20], the bifurcation technique [21], the ( G / G ) -expansion method [22], the Darboux transformation [23], the sine-Gordon equation expansion technique [24] and Kudryashov’s method [25], just to mention a few.
The existent physical models are basically related via NLEEs depending on space variables in three dimensions and one other time variable. Furthermore, the paramount goal of the mathematicians has moved away from nonlinear models involving two-dimensions or three-dimensions [26,27] to non-linear models with dimensions higher than three [28,29,30].
Here, in this work our target is to consider the generalized (3+1)-dimensional breaking soliton Equation (3D-gBSe) [31]
Δ u x t + α u x ( u x y + u x z ) + β u x x ( u y + u z ) + γ ( u x x x y + u x x x z ) = 0
with parameters α , β and γ which are nonzero real-valued constants. If γ = 1 , β = 2 , α = 4 and u = u ( x , y , t ) , Equation (2) converts to two-dimensional breaking soliton equation [32]:
u x t 4 u x u x y 2 u x x u y u x x x y = 0 .
A good number of different kinds of results have been obtained in the literature for breaking solitons in two dimensions [33,34,35] that describe the interconnections between the long wave and the Riemann wave, which propagate respectively along x-axis and y-axis.
It is pertinent to discuss briefly some special cases of the equation under study. If γ = 1 , β = 2 , α = 4 and u = u ( x , y , t ) , then Equation (2) reduces to a potential form of two-dimensional Calogero–Bogoyavlenskii–Schiff equation [36]. With constants in Equation (2) taken as γ = 1 , β = 2 , α = 4 and the dependent variable u = u ( x , y , t ) , Darvishi et al. [37] obtained the solution of Equation (2) by employing a three-wave method. Al-Amr [38] constructed closed-form solutions of another version of (2) when γ = 1 , β = 4 , α = 4 with u independent of t via modified simple equation method. Pallavi et al. [31] found some exact solutions of 3D-gBSe (2) by engaging the extended exp ( ϕ ( ξ ) ) -expansion method. In addition, the authors in [39] engaged the symmetry analysis to procure the residual symmetry of a special case of Equation (2) when γ = 1 , β = 2 and α = 4 . They further proved that the (3+1)-dimensional breaking soliton equation is integrable in the sense of its consistent possession of Riccati expansion.
Symmetry analysis is one of the most systematic techniques with which to find closed-form solutions of differential equations. It was pioneered, towards the end of the nineteenth century by a Norwegian mathematician, Sophus Lie (1842–1899), who realized that the ad hoc methods for solving differential equations could be unified. This technique has found applications in many areas of mathematically-based sciences [19,20].
Conservation laws are basic laws of nature that do not change in an isolated system and have a wide range of applications in physics and in other fields, such as engineering, applied sciences and so on. They can be used in the study of numerical techniques and as well as to determine whether the solution of a differential equation exists and is unique [40,41,42,43,44,45].
In this work, we catalog our article in the following way. In Section 2, we carry out symmetry analysis of 3D-gBSe (2) and perform symmetry reduction of the equation. We find closed-form solutions in the form of elliptic functions by direct integration of the reduced ordinary differential equation. Moreover, the ( G / G ) -expansion method and the power series method are employed to obtain more solutions of Equation (2). Furthermore, in Section 3, we present conserved vectors of the underlying equation by implementing the general multiplier method and Ibragimov’s theorem. In the end, concluding remarks are given.

2. Solutions of the 3D-gBSe (2)

This section is concerned with the construction of closed-form solutions of 3D-gBSe (2) which are obtained by using Lie symmetry analysis, the ( G / G ) -expansion method and the power series approach.

2.1. Lie Symmetry Analysis of Equation (2)

Here we calculate the Lie point symmetries of (2) and eventually utilize them to generate exact solutions for the Equation (2).

2.1.1. Lie Point Symmetries of Equation (2)

We take into consideration a one-parameter Lie group of infinitesimal transformations furnished with a parameter ( ϵ ) that acts on the independent and dependent variables of the equation, thereby keeping equation invariant. The infinitesimal transformations are
x ˜ = x + ϵ ξ 1 ( x , y , z , t , u ) + O ( ϵ 2 ) , y ˜ = y + ϵ ξ 2 ( x , y , z , t , u ) + O ( ϵ 2 ) , z ˜ = z + ϵ ξ 3 ( x , y , z , t , u ) + O ( ϵ 2 ) , t ˜ = t + ϵ ξ 4 ( x , y , z , t , u ) + O ( ϵ 2 ) , u ˜ = u + ϵ ϕ ( x , y , z , t , u ) + O ( ϵ 2 ) ,
wherein the parameter ϵ 1 (a small expansion term) is the group parameter; and ξ 1 , ξ 2 , ξ 3 , ξ 4 and ϕ are infinitesimal coefficients. Thus, related infinitesimal Lie algebra is spanned by the vector field
W = ξ 1 ( x , y , z , t , u ) x + ξ 2 ( x , y , z , t , u ) y + ξ 3 ( x , y , z , t , u ) z + ξ 4 ( x , y , z , t , u ) t + ϕ ( x , y , z , t , u ) u .
The vector field W secures a symmetry of 3D-gBSe (2) provided
pr ( 4 ) W ( Δ ) | Δ = 0 = 0 .
Here pr ( 4 ) W denotes the fourth prolongation of W , which is defines as [19]
pr ( 4 ) W = W + ϕ x u x + ϕ y u y + ϕ z u z + ϕ t u t + ϕ x x u x x + ϕ x y u x y + ϕ x z u x z + ϕ x t u x t + ϕ x x x y u x x x y + ϕ x x x z u x x x z ,
with ϕ x , ϕ y , ϕ z , ϕ t , ϕ x t , ϕ x x , ϕ x y , ϕ x z , ϕ x x x y and ϕ x x x z given by
ϕ x = D x ( ϕ ) u x D x ( ξ 1 ) u y D x ( ξ 2 ) u z D x ( ξ 3 ) u t D x ( ξ 4 ) , ϕ y = D y ( ϕ ) u x D y ( ξ 1 ) u y D y ( ξ 2 ) u z D y ( ξ 3 ) u t D y ( ξ 4 ) , ϕ z = D z ( ϕ ) u x D z ( ξ 1 ) u y D z ( ξ 2 ) u z D z ( ξ 3 ) u t D z ( ξ 4 ) , ϕ t = D t ( ϕ ) u x D t ( ξ 1 ) u y D t ( ξ 2 ) u z D t ( ξ 3 ) u t D t ( ξ 4 ) , ϕ x t = D x ( ϕ t ) u x t D x ( ξ 1 ) u y t D x ( ξ 2 ) u z t D x ( ξ 3 ) u t t D x ( ξ 4 ) , ϕ x x = D x ( ϕ x ) u x x D x ( ξ 1 ) u x y D x ( ξ 2 ) u x z D x ( ξ 3 ) u x t D x ( ξ 4 ) , ϕ x y = D x ( ϕ y ) u x y D x ( ξ 1 ) u y y D x ( ξ 2 ) u z y D x ( ξ 3 ) u y t D x ( ξ 4 ) , ϕ x z = D x ( ϕ z ) u x z D x ( ξ 1 ) u z y D x ( ξ 2 ) u z z D x ( ξ 3 ) u z t D x ( ξ 4 ) , ϕ x x x y = D x ( ϕ x x y ) u x x x y D x ( ξ 1 ) u x x y y D x ( ξ 2 ) u z y x x D x ( ξ 3 ) u x x y t D x ( ξ 4 ) , ϕ x x x z = D x ( ϕ x x z ) u x x x z D x ( ξ 1 ) u z y x x D x ( ξ 2 ) u z z x x D x ( ξ 3 ) u z x x t D x ( ξ 4 ) ,
and D x , D y , D z and D t are total derivatives with respect to independent variables x, y, z and t respectively [19]. For illustrative purposes, it can be indicated that
D i = x i + u i u + u i j u j + , i , j = 1 , 2 , 3 , 4
where ( x 1 , x 2 , x 3 , x 4 ) = ( x , y , z , t ) . Thus, by expanding Equation (3), one gets an overdetermined system of linear PDEs:
ξ u 1 = 0 , ξ u 3 = 0 , ξ u 2 = 0 , ξ x 3 = 0 , ξ u 4 = 0 , ξ x 4 = 0 , ξ y 4 + ξ z 4 = 0 , ξ t 2 α ϕ x = 0 , ξ x 2 = 0 , ϕ t x = 0 , ϕ u t = 0 , ϕ x x = 0 , ϕ x y + ϕ x z = 0 , ϕ u x = 0 , ϕ u y + ϕ u z = 0 , ξ y y 3 + 2 ξ y z 3 + ξ z z 3 = 0 , ϕ y y + ϕ z z + 2 ϕ y z = 0 , ϕ u u = 0 ξ t 4 + 2 ϕ u ξ z 3 ξ y 3 = 0 , ξ t 1 β ( ϕ y + ϕ z ) = 0 , ξ t 3 α ϕ x = 0 , ξ x 1 + ϕ u = 0 , ξ y 1 + ξ z 1 = 0 , ξ y 2 + ξ z 2 ξ z 3 ξ y 3 = 0 .
The solutions of the system give infinitesimal generators of one-parameter Lie group of symmetries for Equation (2) as
ξ 1 = x F 2 ( z y ) + β F 4 ( t , z y ) + F 8 ( z y ) , ξ 2 = α t F 1 ( z y ) + y F 6 ( z y ) + F 7 ( z y ) , ξ 3 = α t F 1 ( z y ) + F 5 ( z y ) + y F 6 ( z y ) , ξ 4 = 2 t F 2 ( z y ) + t F 6 ( z y ) + F 9 ( z y ) ϕ = x F 1 ( z y ) + u F 2 ( z y ) + F 3 ( t , z y ) + y F t 4 ( t , z y ) ,
where F i , i = 1 , 2 , , 9 are arbitrary functions of their arguments. Hence, the Lie algebra of infinitesimal symmetries of Equation (2) is spanned by the vector fields
W 1 = F 8 ( z y ) x , W 2 = F 7 ( z y ) y , W 3 = F 5 ( z y ) z , W 4 = F 9 ( z y ) t , W 5 = F 3 ( t , z y ) u , W 6 = β F 4 ( t , z y ) x + y F t 4 ( t , z y ) u , W 7 = F 2 ( z y ) x x + 2 t t u u , W 8 = F 6 ( z y ) y y + y z + t t , W 9 = F 1 ( z y ) α t y + α t z + x u .
If we take F 1 , , F 9 to be equal to 1, then 3D-gBSe (2) admits an eight-dimensional Lie algebra L 8 spanned by the vectors
X 1 = x , X 2 = y , X 3 = z , X 4 = t , X 5 = u , X 6 = 2 t t + x x u u , X 7 = y y + y z + t t , X 8 = α t y + α t z + x u .

2.1.2. Symmetry Reduction of 3D-gBSe (2)

Here we perform symmetry reduction of the 3D-gBSe (2). In the first place, we use the symmetry X = X 1 + ν X 2 + X 3 + X 4 with nonzero constant ν to reduce the 3D-gBSe (2) to a PDE containing three independent variables. As a result, solving the associated Lagrange system for X gives rise to four invariants, viz.,
f = y ν t , g = t x , h = y ν z , θ = u .
Now regarding θ as the new dependent variable, and f, g and h as new independent variables, the 3D-gBSe (2) then transforms into
ν θ g f θ g g + α ( θ g θ g f + θ g θ g h ν θ g θ g h ) + β ( θ f θ g g + θ h θ g g ν θ h θ g g ) + γ ( ν θ g g g h θ g g g f θ g g g h ) = 0 ,
which is a nonlinear PDE in three independent variables. By exploiting the symmetries of Equation (9), we transform Equation (9) into a PDE with two independent variables. Equation (9) possesses eight point symmetries, which include the three translation symmetries Σ 1 = / f , Σ 2 = / g and Σ 3 = / h . The symmetry Σ = Σ 1 + Σ 2 + Σ 3 has three invariants r = f g , s = h g and θ = H , and using these invariants, Equation (9) transforms to
H r r 2 H r s H s s ν H r r ν H r s + α ( H r + H s ) [ ( H r r + 2 H r s + H s s ) ν ( H r s + H s s ) ] + β ( H r r + 2 H r s + H s s ) H r + H s ν H s + γ [ ( H r r r r + 6 H r r s s + 4 H r r r s + 4 H r s s s + H s s s s ) ν H r r r s + 3 H r s s s + 3 H r r s s + H s s s s ] = 0 ,
which is a NLPE in two independent variables. Equation (10) has four point symmetries that include Ω 1 = / r , Ω 2 = / s . The symmetry Ω = Ω 1 + ω Ω 2 , with ω 0 has two invariants ζ = s ω r a n d H = G , which leads to the group-invariant solution H = G ( ζ ) , where G satisfies the fourth-order nonlinear ordinary differential equation (NLODE)
( ω 2 + 2 ω + ν ω ν ω 2 1 ) G α β ν ( ( ω 1 ) ( ω 2 2 ω + 1 ) ( 1 ω ν ) ) G G + γ ( 7 ω 2 + ν ω 3 4 ω 3 4 ω + 3 ν ω 3 ν ω 2 ν + 1 ) G = 0 .
For simplicity, we write Equation (11) as
P G ( ζ ) Q G ( ζ ) G ( ζ ) + R G ( ζ ) = 0 ,
where P = ω ( ω ( 1 ν ) + ν + 2 ) 1 , Q = α β ν ( ( ω 1 ) ( ω 2 2 ω + 1 ) ( 1 ω ν ) ) , R = γ ( ω 2 ( ω ( ν 4 ) 3 ν + 7 ) + ω ( 3 ν 4 ) ν + 1 ) and ζ = ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z .

2.1.3. Solution of Equation (2) by Direct Integration of Equation (12)

Integrating Equation (12) with respect to ζ once gets a third-order ODE
P G 1 2 Q ( G ) 2 + R G + C 1 = 0 ,
where C 1 is a constant of integration. By multiplying Equation (13) by G , and integrating and simplifying the resulting equation, we have the second-order NLODE
1 2 P ( G ) 2 1 6 Q ( G ) 3 + 1 2 R ( G ) 2 + C 1 G + C 2 = 0
with C 2 an integration constant. Equation (14) can be rewritten as
( G ) 2 = Q 3 R ( G ) 3 P R ( G ) 2 2 C 1 R G 2 C 2 R .
Suppose Θ = G . Equation (15) becomes
Θ 2 = Q 3 R Θ 3 P R Θ 2 2 C 1 R Θ 2 C 2 R .
We consider the cubic equation
Θ 3 3 P Q Θ 2 6 C 1 Q Θ 6 C 2 Q = 0
and assume that the roots of the equation are r 1 , r 2 and r 3 such that r 3 r 2 r 1 . Then Equation (16) becomes
Θ 2 = Q 3 R ( Θ r 1 ) ( Θ r 2 ) ( Θ r 3 )
and the solution to Equation (16) can be expressed in terms of Jacobi elliptic function [46,47,48] as
Θ ( ζ ) = r 2 + ( r 1 r 2 ) c n 2 Q ( r 1 r 2 ) 12 R ζ , S 2 , S 2 = r 1 r 2 r 1 r 3
with cn denoting the cosine elliptic function. By integrating Equation (19) and then returning to the original variables ( x , y , z , t ) , the 3D-gBSe (2) possesses a periodic solution
u ( x , y , z , t ) = 12 R ( r 1 r 2 ) 2 Q ( r 1 r 3 ) S 8 EllipticE s n Q ( r 1 r 3 ) 12 R ζ , S 2 , S 2 + r 2 ( r 1 r 2 ) 1 S 4 S 4 ζ + C 3 ,
where sn represents the sine elliptic function, ζ = ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z , C 3 is an integration constant and EllipticE [ q , v ] denotes the incomplete elliptic integral given by [48]
EllipticE [ q , v ] = 0 v 1 v 2 n 2 1 n 2 d n .
To view the dynamics of solution (20) graphically, we sketched its graphs in Figure 1, Figure 2 and Figure 3, for certain values of the parameters and for t = 10 and z = 1 t = 11 and z = 0 and t = 25 and z = 9 , respectively.

2.2. Exact Solutions of Equation (2) by ( G / G ) -Expansion Method

This subsection engages the ( G / G ) -expansion method [22] to get exact solutions of the 3D-gBSe (2). We consider a solution of the form
G ( ζ ) = j = 0 M B j H ( ζ ) H ( ζ ) j ,
where H ( ζ ) satisfies
H ( ζ ) + λ H ( ζ ) + μ H ( ζ ) = 0
with λ and μ taken as constants. In this case, the constant M is generated by balancing the nonlinear term with the highest order as well as the highest order derivatives emerging in Equation (2). Here B 0 , , B M , are parameters to be determined.
Engagement of the balancing procedure for Equation (11) gives M = 1 , and consequently, the solution of Equation (11) assumes the form
G ( ζ ) = B 0 + B 1 H ( ζ ) H ( ζ ) .
Substituting the value of G ( ζ ) from Equation (23) into Equation (11) and using Equation (22) leads to an algebraic equation in B 0 and B 1 , which splits over various powers of H ( ζ ) to give a system of algebraic equations:
α λ μ 2 ν ω 2 B 1 2 2 α λ μ 2 ν ω B 1 2 + α λ μ 2 ν B 1 2 + α λ μ 2 ω 3 B 1 2 3 α λ μ 2 ω 2 B 1 2 + 3 α λ μ 2 ω B 1 2 α λ μ 2 B 1 2 + β λ μ 2 ν ω 2 B 1 2 2 β λ μ 2 ν ω B 1 2 + β λ μ 2 ν B 1 2 + β λ μ 2 ω 3 B 1 2 3 β λ μ 2 ω 2 B 1 2 + 3 β λ μ 2 ω B 1 2 β λ μ 2 B 1 2 + γ λ 3 μ ν ω 3 B 1 3 γ λ 3 μ ν ω 2 B 1 + 3 γ λ 3 μ ν ω B 1 γ λ 3 μ ν B 1 + γ λ 3 μ ω 4 B 1 4 γ λ 3 μ ω 3 B 1 + 6 γ λ 3 μ ω 2 B 1 4 γ λ 3 μ ω B 1 + γ λ 3 μ B 1 + 8 γ λ μ 2 ν ω 3 B 1 24 γ λ μ 2 ν ω 2 B 1 + 24 γ λ μ 2 ν ω B 1 8 γ λ μ 2 ν B 1 + 8 γ λ μ 2 ω 4 B 1 32 γ λ μ 2 ω 3 B 1 + 48 γ λ μ 2 ω 2 B 1 32 γ λ μ 2 ω B 1 + 8 γ λ μ 2 B 1 λ μ ν ω 2 B 1 + λ μ ν ω B 1 λ μ ω 2 B 1 + 2 λ μ ω B 1 λ μ B 1 = 0 , 2 α λ 2 μ ν ω 2 B 1 2 4 α λ 2 μ ν ω B 1 2 + 2 α λ 2 μ ν B 1 2 + 2 α λ 2 μ ω 3 B 1 2 6 α λ 2 μ ω 2 B 1 2 + 6 α λ 2 μ ω B 1 2 2 α λ 2 μ B 1 2 + 2 α μ 2 ν ω 2 B 1 2 4 α μ 2 ν ω B 1 2 + 2 α μ 2 ν B 1 2 + 2 α μ 2 ω 3 B 1 2 6 α μ 2 ω 2 B 1 2 + 6 α μ 2 ω B 1 2 2 α μ 2 B 1 2 + 2 β λ 2 μ ν ω 2 B 1 2 4 β λ 2 μ ν ω B 1 2 + 2 β λ 2 μ ν B 1 2 + 2 β λ 2 μ ω 3 B 1 2 6 β λ 2 μ ω 2 B 1 2 + 6 β λ 2 μ ω B 1 2 2 β λ 2 μ B 1 2 + 2 β μ 2 ν ω 2 B 1 2 4 β μ 2 ν ω B 1 2 + 2 β μ 2 ν B 1 2 + 2 β μ 2 ω 3 B 1 2 6 β μ 2 ω 2 B 1 2 + 6 β μ 2 ω B 1 2 2 β μ 2 B 1 2 + γ λ 4 ν ω 3 B 1 3 γ λ 4 ν ω 2 B 1 + 3 γ λ 4 ν ω B 1 γ λ 4 ν B 1 + γ λ 4 ω 4 B 1 4 γ λ 4 ω 3 B 1 + 6 γ λ 4 ω 2 B 1 4 γ λ 4 ω B 1 + γ λ 4 B 1 + 22 γ λ 2 μ ν ω 3 B 1 66 γ λ 2 μ ν ω 2 B 1 + 66 γ λ 2 μ ν ω B 1 22 γ λ 2 μ ν B 1 + 22 γ λ 2 μ ω 4 B 1 88 γ λ 2 μ ω 3 B 1 + 132 γ λ 2 μ ω 2 B 1 88 γ λ 2 μ ω B 1 + 22 γ λ 2 μ B 1 + 16 γ μ 2 ν ω 3 B 1 48 γ μ 2 ν ω 2 B 1 + 48 γ μ 2 ν ω B 1 16 γ μ 2 ν B 1 + 16 γ μ 2 ω 4 B 1 64 γ μ 2 ω 3 B 1 + 96 γ μ 2 ω 2 B 1 64 γ μ 2 ω B 1 + 16 γ μ 2 B 1 λ 2 ν ω 2 B 1 + λ 2 ν ω B 1 λ 2 ω 2 B 1 + 2 λ 2 ω B 1 λ 2 B 1 2 μ ν ω 2 B 1 + 2 μ ν ω B 1 2 μ ω 2 B 1 + 4 μ ω B 1 2 μ B 1 = 0 , α λ 3 ν ω 2 B 1 2 2 α λ 3 ν ω B 1 2 + α λ 3 ν B 1 2 + α 3 ω 3 B 1 2 3 α λ 3 ω 2 B 1 2 + 3 α λ 3 ω B 1 2 α λ 3 B 1 2 + 6 α λ μ ν ω 2 B 1 2 12 α λ μ ν ω B 1 2 + 6 α λ μ ν B 1 2 + 6 α λ μ ω 3 B 1 2 18 α λ μ ω 2 B 1 2 + 18 α λ μ ω B 1 2 6 α λ μ B 1 2 + β λ 3 ν ω 2 B 1 2 2 β λ 3 ν ω B 1 2 + β λ 3 ν B 1 2 + β λ 3 ω 3 B 1 2 3 β λ 3 ω 2 B 1 2 + 3 β λ 3 ω B 1 2 β λ 3 B 1 2 + 6 β λ μ ν ω 2 B 1 2 12 β λ μ ν ω B 1 2 + 6 β λ μ ν B 1 2 + 6 β λ μ ω 3 B 1 2 18 β λ μ ω 2 B 1 2 + 18 β λ μ ω B 1 2 6 β λ μ B 1 2 + 15 γ λ 3 ν ω 3 B 1 45 γ λ 3 ν ω 2 B 1 + 45 γ λ 3 ν ω B 1 15 γ λ 3 ν B 1 + 15 γ λ 3 ω 4 B 1 60 γ λ 3 ω 3 B 1 + 90 γ λ 3 ω 2 B 1 60 γ λ 3 ω B 1 + 15 γ λ 3 B 1 + 60 γ λ μ ν ω 3 B 1 180 γ λ μ ν ω 2 B 1 + 180 γ λ μ ν ω B 1 60 γ λ μ ν B 1 + 60 γ λ μ ω 4 B 1 240 γ λ μ ω 3 B 1 + 360 γ λ μ ω 2 B 1 240 γ λ μ ω B 1 + 60 γ λ μ B 1 3 λ ν ω 2 B 1 + 3 λ ν ω B 1 3 λ ω 2 B 1 + 6 λ ω B 1 3 λ B 1 = 0 , 4 α λ 2 ν ω 2 B 1 2 8 α λ 2 ν ω B 1 2 + 4 α λ 2 ν B 1 2 + 4 α λ 2 ω 3 B 1 2 12 α λ 2 ω 2 B 1 2 + 12 α λ 2 ω B 1 2 4 α λ 2 B 1 2 + 4 α μ ν ω 2 B 1 2 8 α μ ν ω B 1 2 + 4 α μ ν B 1 2 + 4 α μ ω 3 B 1 2 12 α μ ω 2 B 1 2 + 12 α μ ω B 1 2 4 α μ B 1 2 + 4 β λ 2 ν ω 2 B 1 2 8 β λ 2 ν ω B 1 2 + 4 β λ 2 ν B 1 2 + 4 β λ 2 ω 3 B 1 2 12 β λ 2 ω 2 B 1 2 + 12 β λ 2 ω B 1 2 4 β λ 2 B 1 2 + 4 β μ ν ω 2 B 1 2 8 β μ ν ω B 1 2 + 4 β μ ν B 1 2 + 4 β μ ω 3 B 1 2 12 β μ ω 2 B 1 2 + 12 β μ ω B 1 2 4 β μ B 1 2 + 50 γ λ 2 ν ω 3 B 1 150 γ λ 2 ν ω 2 B 1 + 150 γ λ 2 ν ω B 1 50 γ λ 2 ν B 1 + 50 γ λ 2 ω 4 B 1 200 γ λ 2 ω 3 B 1 + 300 γ λ 2 ω 2 B 1 200 γ λ 2 ω B 1 + 50 γ λ 2 B 1 + 40 γ μ ν ω 3 B 1 120 γ μ ν ω 2 B 1 + 120 γ μ ν ω B 1 40 γ μ ν B 1 + 40 γ μ ω 4 B 1 160 γ μ ω 3 B 1 + 240 γ μ ω 2 B 1 160 γ μ ω B 1 + 40 γ μ B 1 2 ν ω 2 B 1 + 2 ν ω B 1 2 ω 2 B 1 + 4 ω B 1 2 B 1 = 0 , 5 α λ ν ω 2 B 1 2 10 α λ ν ω B 1 2 + 5 α λ ν B 1 2 + 5 α λ ω 3 B 1 2 15 α λ ω 2 B 1 2 + 15 α λ ω B 1 2 5 α λ B 1 2 + 5 β λ ν ω 2 B 1 2 10 β λ ν ω B 1 2 + 5 β λ ν B 1 2 + 5 β λ ω 3 B 1 2 15 β λ ω 2 B 1 2 + 15 β λ ω B 1 2 5 β λ B 1 2 + 60 γ λ ν ω 3 B 1 180 γ λ ν ω 2 B 1 + 180 γ λ ν ω B 1 60 γ λ ν B 1 + 60 γ λ ω 4 B 1 240 γ λ ω 3 B 1 + 360 γ λ ω 2 B 1 240 γ λ ω B 1 + 60 γ λ B 1 = 0 , 2 α ν ω 2 B 1 2 4 α ν ω B 1 2 + 2 α ν B 1 2 + 2 α ω 3 B 1 2 6 α ω 2 B 1 2 + 6 α ω B 1 2 2 α B 1 2 + 2 β ν ω 2 B 1 2 4 β ν ω B 1 2 + 2 β ν B 1 2 + 2 β ω 3 B 1 2 6 β ω 2 B 1 2 + 6 β ω B 1 2 2 β B 1 2 + 24 γ ν ω 3 B 1 72 γ ν ω 2 B 1 + 72 γ ν ω B 1 24 γ ν B 1 + 24 γ ω 4 B 1 96 γ ω 3 B 1 + 144 γ ω 2 B 1 96 γ ω B 1 + 24 γ B 1 = 0 .
The solution of the above system of equations via Mathematica is given by
B 0 = B 0 , B 1 = 12 γ ( ω 1 ) α + β and γ = ν ω + ω 1 ( ω 1 ) 2 λ 2 4 μ ( ν + ω 1 ) .
Thus, we have two types of solutions of the 3D-gBSe (2) according to the sign of λ 2 4 μ :
When λ 2 4 μ > 0 , we gain the hyperbolic function solution
u ( x , y , z , t ) = B 0 + B 1 λ 2 + Δ 1 C 1 sinh Δ 1 ζ + C 2 cosh Δ 1 ζ C 1 cosh Δ 1 ζ + C 2 sinh Δ 1 ζ ,
where ζ = ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z , Δ 1 = 1 2 λ 2 4 μ together with C 1 and C 2 arbitrary constants.
When λ 2 4 μ < 0 , we achieve the trigonometric function solution
u ( x , y , z , t ) = B 0 + B 1 λ 2 + Δ 2 C 1 sin Δ 2 ζ + C 2 cos Δ 2 ζ C 1 cos Δ 2 ζ + C 2 sin Δ 2 ζ ,
where ζ = ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ; Δ 2 = 1 2 4 μ λ 2 ; and C 1 and C 2 are arbitrary constants.
In a bid to see the behavior of solution (24), we sketched graphs of solution (24) for various parameter values in Figure 4, Figure 5 and Figure 6 t = 1 and z = 0 ; t = 2 and z = 3.5 ; and t = 0.2 and z = 20 , respectively.
Likewise, in Figure 7, Figure 8 and Figure 9 we present the profiles of solution (25) for various parameter values, for t = 2.1 and z = 3 ; t = 6 and z = 25 ; and t = 4 and z = 30 , respectively.

2.3. The Jacobi Elliptic Function Solutions of Equation (2)

In this subsection, by using the auxiliary ordinary differential equations, we obtain the Jacobi function solutions to the 3D-gBSe (2) with the help of symbolic computation system Maple. Recall that the cosine-amplitude function, cn ( ζ ; m ) , and the sine amplitude function, sn ( ζ ; m ) , satisfy, respectively [48], the ordinary differential equations:
H ( ζ ) + 2 m H ( ζ ) 3 + ( 1 2 m ) H ( ζ ) = 0 , H ( ζ ) 2 m H ( ζ ) 3 + ( m + 1 ) H ( ζ ) = 0 .
We now engage the ODEs (26) instead of the ODE (22) in the ( G / G ) -expansion technique. As a consequence, following the procedure of the previous subsection we secure new closed-form solutions of the 3D-gBSe (2), which are periodic solutions in terms of the Jacobi elliptic functions; that is,
u ( x , y , z , t ) = B 0 B 1 dn ( ζ ; m ) sn ( ζ ; m ) cn ( ζ ; m ) ,
with parameters
B 0 = B 0 , B 1 = 12 γ ( ω 1 ) α + β , γ = 1 ω ν ω 4 ( ω 1 ) 2 ( ω + ν 1 ) , m = 4 γ ν 4 γ + ν ω + 12 γ ω 8 γ ν ω 12 γ ω 2 + 4 γ ν ω 2 + 4 γ ω 3 + ω 1 8 γ ( ω 1 ) 2 ( ω + ν 1 )
and
u ( x , y , z , t ) = B 0 + B 1 cn ( ζ ; m ) dn ( ζ ; m ) sn ( ζ ; m ) ,
where we have
B 0 = B 0 , B 1 = 12 γ ( ω 1 ) α + β , γ = 1 ω ν ω 4 ( ω 1 ) 2 ( ω + ν 1 ) , m = 4 γ 4 γ ν ν ω 12 γ ω + 8 γ ν ω + 12 γ ω 2 4 γ ν ω 2 4 γ ω 3 ω + 1 4 γ ( ω 1 ) 2 ( ω + ν 1 ) .
The profiles of solutions (27) and (28), for different values of the parameters and t and x, are shown in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, respectively.

2.4. Solution of Equation (2) Using the Power Series Method

Here we shall present the analytic solution to ordinary differential Equation (12) by utilizing the power series method [49,50,51]. In this case we procure a solution of Equation (12) in a power series which assumes the structure
G ( ζ ) = p = 0 c p ζ p ,
where c 0 = G ( 0 ) 0 . Furthermore
G ( ζ ) = p = 1 p c p ζ p 1 , G ( ζ ) = p = 2 p ( p 1 ) c p ζ p 2 , G ( ζ ) = p = 4 p ( p 1 ) ( p 2 ) ( p 3 ) c p ζ p 4 .
By inserting the values of G , G and G into Equation (12), the equation is transformed to
P p = 0 ( p + 1 ) ( p + 2 ) c p + 2 ζ p Q p = 0 ( p + 1 ) c p + 1 ζ p p = 0 ( p + 1 ) ( p + 2 ) c p + 2 ζ p + R p = 0 ( p + 1 ) ( p + 2 ) ( p + 3 ) ( p + 4 ) c p + 4 ζ p = 0 ,
which simplifies further to
P 2 c 2 + p = 1 ( p + 1 ) ( p + 2 ) c p + 2 ζ p Q ( 2 c 1 c 2 p = 1 k = 0 p ( k + 1 ) ( p + 1 k ) ( p + 2 k ) c k + 1 c p + 2 k ζ p ) + R 24 c 4 + p = 1 ( p + 1 ) ( p + 2 ) ( p + 3 ) ( p + 4 ) c p + 4 ζ p = 0 .
Comparing coefficients in Equation (31), for p = 0 , one obtains
c 4 = c 1 c 2 Q c 2 P 12 R .
In a more general way, if we consider p 1 , we get a recursive relation
c p + 4 = 1 R ( p + 1 ) ( p + 2 ) ( p + 3 ) ( p + 4 ) { Q k = 0 p ( k + 1 ) ( p + 1 k ) ( p + 2 k ) c k + 1 c p + 2 k P ( p + 1 ) ( p + 2 ) c p + 2 } .
One can calculate all the coefficients c p ( p 4 ) from Equations (32) and recursive Formula (33). Successive determination of other terms can as well be made from Equation (32) and Equation (33) in a unique way. Furthermore, it can be ascertained that power series Equation (29) converges along side the coefficients c p s stated in Equation (32) and (33), which will not be considered here see [50]. Therefore, this power series solution can be considered as a closed-form solution. Thus,
G ( ζ ) = c 0 + c 1 ζ + c 2 ζ 2 + c 3 ζ 3 + c 4 ζ 4 + p = 1 c p + 4 ζ p + 4 , = c 0 + c 1 ζ + c 2 ζ 2 + c 3 ζ 3 + c 1 c 2 Q c 2 P 12 R ζ 4 + p = 1 1 R ( p + 1 ) ( p + 2 ) ( p + 3 ) ( p + 4 ) { Q k = 0 p ( k + 1 ) ( p + 1 k ) ( p + 2 k ) c k + 1 c p + 2 k P ( p + 1 ) ( p + 2 ) c p + 2 } ζ p + 4 .
Therefore, the power series solution of the 3D-gBSe (2) is
u ( x , y , z , t ) = c 0 + c 1 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) + c 2 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 2 + c 3 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 3 + c 4 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 4 + p = 1 c p + 4 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) p + 4 = c 0 + c 1 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) + c 2 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 2 + c 3 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 3 + Q c 1 c 2 P c 2 12 R ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 4 + p = 1 1 R ( p + 1 ) ( p + 2 ) ( p + 3 ) ( p + 4 ) { Q k = 0 p ( k + 1 ) ( p + 1 k ) ( p + 2 k ) c k + 1 c p + 2 k P ( p + 1 ) ( p + 2 ) c p + 2 } ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) p + 4 ,
where c i , ( i = 0 , 1 , 2 , 3 , 4 ) are constants. Other relevant coefficients c p ( p 4 ) can be then obtained from Equations (32) and (33). In physical applications, we can approximately give the series solution as
u ( x , y , z , t ) = c 0 + c 1 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) + c 2 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 2 + c 3 ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 3 + Q c 1 c 2 P c 2 12 R ( ( ω ( ν + 1 ) 1 ) t ( ω 1 ) x + ( 1 ω ) y ν z ) 4 + ,
which is the power series solution of the 3D-gBSe (2).
To demonstrate the behaviour of power series solution (34) pictorially, we present its profile in Figure 16, Figure 17 and Figure 18 for certain values of the parameters, for t = 1.6 and z = 1.1 ; t = 1 and z = 8 ; and t = 5 and z = 4 , respectively.
Remark 1.
We note here that the assumption that c 0 = G ( 0 ) 0 in the power series solution (29) is germane for our argument; otherwise, we would not be able to secure exact analytic solution for the underlying equation by using this method.

3. Conserved Vectors of 3D-gBSe (2)

This section provides the conserved vectors for the 3D-gBSe (2). In this study, we utilize the general multiplier technique [19] and the Ibragimov’s theorem [44] to derive the conserved vectors.

3.1. Multiplier Technique

The multiplier approach happens to be one of the most sturdy and robust in deriving the conserved vectors of a differential equation [19,52,53,54,55,56]. We seek the third-order multipliers represented by Λ and so write down its determining equation [56]
δ δ u [ Λ ( u x t + α u x ( u x y + u x z ) + β u x x ( u y + u z ) + γ ( u x x x y + u x x x z ) ) ] = 0
with the Euler operator
δ δ u = D x u x D y u y D z u z + D x 2 u x x + D x D y u x y + D x D z u x z + D x D t u x t + D x 3 D y u x x x y + D x 3 D z u x x x z .
The expansion of Equation (35), after splitting and simplifying the result with regard to the various derivatives of ( u ) involved, yields equations which when solved give the value of Equation Λ as determining
Λ = F 1 ( t , z y ) + u x F 2 ( t , z y ) + 1 β α y F t 2 ( t , z y ) + u x x x F 3 ( z y ) + α + β 2 γ u x 2 F 3 ( z y ) ,
where F 1 , F 2 and F 3 are arbitrary functions. The conserved quantities are then obtained by invoking the divergence identity D t T t + D x T x + D y T y + D z T z = Δ Λ with T t as conserved density and ( T x , T y , T z ) being spatial fluxes [40]. Thus, after some calculations, we obtain the following three local conserved vectors of (2):
Case 1 . Λ 1 = F 1 ( t , z y ) . C 1 t = 1 2 F 1 ( t , z y ) u x , C 1 x = 1 4 { [ ( α 2 β ) ( u x y + u x z ) u + α ( u y + u z ) u x + 2 β ( u y + u z ) u x + 3 γ ( u x x y + u x x z ) + 2 u t ] F 1 ( t , z y ) 2 F t 1 ( t , z y ) u } , C 1 y = 1 4 α ( u x 2 u u x x ) + 2 β u u x x + γ u x x x F 1 ( t , z y ) , C 1 z = 1 4 α ( u x 2 u u x x ) + 2 β u u x x + γ u x x x F 1 ( t , z y ) . Case 2 . Λ 2 = 1 β α y F t 2 ( t , z y ) + u x F 2 ( t , z y ) . C 2 t = 1 4 ( β α ) ( α β ) ( u u x x u x 2 ) F 2 ( t , z y ) + 2 F t 2 ( t , z y ) y u x , C 2 x = 1 24 ( β α ) { 9 ( α β ) [ 8 9 β u x u ( u x y + u x z ) u x ( u y + u z ) 8 9 α u x ( 1 2 ( u y + u z ) u x + u ( u x y + u x z ) ) 5 3 u x γ ( u x x y + u x x z ) 2 5 u t + 1 3 γ ( 3 u x x u x y + 3 u x x u x z u ( u x x x y + u x x x z ) u x x x ( u y + u z ) ) 2 3 u u t x ] F 2 ( t , z y ) + ( ( ( 6 u + ( 12 u y + 12 u z ) y ) u x 12 u y ( u x y + u x z ) ) β + 6 y ( ( u y + u z ) u x + u ( u x y + u x z ) ) α + ( ( 18 u x x y + 18 u x x z ) y 6 u x x ) γ + 12 y u t ) F t 2 ( t , z y ) 12 F t t 2 ( t , z y ) u y } , C 2 y = 1 24 ( β α ) { 3 [ α 8 3 u u x u x x 4 3 u x 3 8 3 β u u x u x x + γ ( u u x x x x 2 u x u x x x + u x x 2 ) ] ( β α ) F 2 ( t , z y ) + 6 y ( ( u x 2 u u x x ) α + 2 β u u x x + γ u x x x ) F t 2 ( t , z y ) } , C 2 z = 1 24 ( β α ) { 3 [ α 8 3 u u x u x x 4 3 u x 3 8 3 β u u x u x x + γ ( u u x x x x 2 u x u x x x + u x x 2 ) ] ( β α ) F 2 ( t , z y ) + 6 y ( ( u x 2 u u x x ) α + 2 β u u x x + γ u x x x ) F t 2 ( t , z y ) } . Case 3 . Λ 3 = 1 2 γ ( α + β ) u x 2 + 2 γ u x x x F 3 ( z y ) . C 3 t = 1 4 γ 1 3 u x 3 ( α + β ) + 2 3 ( α + β ) u u x x γ u x x x + γ u u x x x x F 3 ( z y ) , C 3 x = 1 8 γ { F 3 ( z y ) [ γ 2 ( ( 2 u x x x x y + 2 u x x x z ) u ( u x x x x x y + u x x x x x z ) 3 u x x u x x x y 3 u x x u x x x z 3 u x x x ( u x x y + u x x z ) + 2 u x y u x x x x + 2 u x z u x x x x ) γ 2 u y u x x x x x γ 2 u z u x x x x x + γ [ 7 3 β 1 7 α ( u x x y + u x x z ) u x 2 2 α u u x x x y 2 α u u x x x z + u x { α 2 3 u x y + 2 3 u x z u x x 2 u x x x ( u y + u z ) + β [ u 2 3 u x x x y + 2 3 u x x x z + u x x 10 3 u x y + 10 3 u x z 2 3 u x x x ( u y + u z ) ] + 4 u t x x } + α { u [ 4 3 u x x u x x y 4 3 u x x u x x z 2 u x x x ( u x y + u x z ) ] 2 3 u x x 2 ( u y + u z ) } + β { u [ 4 3 u x x u x x y + 4 3 u x x u x x z 2 3 u x x x ( u x y + u x z ) ] 10 3 u x x 2 ( u y + u z ) } 2 u u t x x x 2 u t u x x x 4 u t x u x x ] + ( α + β ) u x { u x 2 β + 1 2 α ( u y + u z ) + u x ( 3 2 α u ( u x y + u x z ) + β u ( u x y + u x z ) 2 3 u t ) 4 3 u u t x } ] } , C 3 y = 1 12 γ { F 3 ( z y ) [ γ 2 3 2 u u x x x x x x + 3 2 u x x u x x x x 3 2 u x u x x x x x 3 2 u x x x 2 + γ 3 2 u x 2 u x x x β + 7 3 α + u u x x x x ( β + 3 α ) u x u u x x u x x x ( β 5 α ) 3 2 1 2 α u x 2 + u u x x β 3 2 α ( α + β ) u x 2 ] } , C 3 z = 1 12 γ { F 3 ( z y ) [ γ 2 3 2 u u x x x x x x + 3 2 u x x u x x x x 3 2 u x u x x x x x 3 2 u x x x 2 + γ 3 2 u x 2 u x x x β + 7 3 α + u u x x x x ( β + 3 α ) u x u u x x u x x x ( β 5 α ) 3 2 1 2 α u x 2 + u u x x β 3 2 α ( α + β ) u x 2 ] } .

3.2. Ibragimov Approach

This approach suggests that every infinitesimal generator is ascribed uniquely to a conserved quantity. We now compute the conserved vectors of 3D-gBSe (2) by exploiting a new conservation theorem by Ibragimov [43,44]. The reader is addressed to the given references for a detailed understanding of the approach. Thus, we have the adjoint of (2) [44] as
H * δ δ u ( v ( u x t + α u x ( u x y + u x z ) + β u x x ( u y + u z ) + γ ( u x x x y + u x x x z ) ) ) = 0
with the Euler operator defined in Equation (36). The expansion of Equation (38) gives
H * v t x + α u x x v y β u x x v y + β u y v x x + α u x v x y + v x z + 2 β v x u x y + u x z + α u x x v z β u x x v z + β u z v x x + γ v x x x y + γ v x x x z = 0 .
The formal Lagrangian of 3D-gBSe (2) alongside its adjoint given in Equation (39) is expressed as
L = v u x t + α u x u x y + β u x x u y + α u x u x z + β u x x u z + γ u x x x y + γ u x x x z .
We recall here that 3D-gBSe (2) admits nine Lie point symmetries stated in Equation (7). Hence, we construct the associated conserved vectors for the stated Lagrangian by making use of the formula [44]
C i = ξ i L + W α L u i α D j L u i j α + D j D k L u i j k α + + D j ( W α ) L u i j α D k L u i j k α + + D j D k ( W α ) L u i j k + ,
with the Lie characteristic function W α expressed as W α = η α ξ j u j α , for α = 1 , 2 and j = 1 , 2 , 3 , 4 . Thus, by the reckoning of the relation in Equation (41), we achieve the conserved vectors [55] ( T i , X i , Y i , Z i ) of 3D-gBSe (2) for the respective symmetries in Equation (7) as
T 1 = 1 2 u x v x F 1 ( z y ) 1 2 u x x v F 1 ( z y ) , X 1 = β u x u x z v F 1 ( z y ) + β u x u x y v F 1 ( z y ) + 1 4 γ u x x x z v F 1 ( z y ) + 1 4 γ u x x x y v F 1 ( z y ) + 1 2 u t x v F 1 ( z y ) + 1 2 v t u x F 1 ( z y ) + 1 2 α u x 2 v z F 1 ( z y ) + 1 2 α u x 2 v y F 1 ( z y ) + β u x u z v x F 1 ( z y ) + β u x u y v x F 1 ( z y ) + 3 4 γ u x v x x z F 1 ( z y ) + 3 4 γ u x v x x y F 1 ( z y ) 1 2 γ u x x v x z F 1 ( z y ) 1 2 γ u x x v x y F 1 ( z y ) 1 4 γ v x x u x z F 1 ( z y ) 1 4 γ v x x u x y F 1 ( z y ) + 1 2 γ v x u x x z F 1 ( z y ) + 1 2 γ v x u x x y F 1 ( z y ) + 1 4 γ u x x x v z F 1 ( z y ) + 1 4 γ u x x x v y F 1 ( z y ) , Y 1 = 1 2 α u x 2 v x F 1 ( z y ) β u x x u x v F 1 ( z y ) 1 4 γ u x x x x v F 1 ( z y ) + 1 4 γ u x v x x x F 1 ( z y ) 1 4 γ u x x v x x F 1 ( z y ) + 1 4 γ u x x x v x F 1 ( z y ) , Z 1 = 1 4 γ u x v x x x F 1 ( z y ) β u x x u x v F 1 ( z y ) 1 4 γ u x x x x v F 1 ( z y ) + 1 2 α u x 2 v x F 1 ( z y ) 1 4 γ u x x v x x F 1 ( z y ) + 1 4 γ u x x x v x F 1 ( z y ) ; T 2 = 1 2 u y v x F 2 ( z y ) 1 2 u x y v F 2 ( z y ) , X 2 = β v x u y 2 F 2 ( z y ) + 1 2 α v z u x u y F 2 ( z y ) + 1 2 α v y u x u y F 2 ( z y ) + β u z v x u y F 2 ( z y ) 1 2 α v u x z u y F 2 ( z y ) + β v u x z u y F 2 ( z y ) 1 2 α v u x y u y F 2 ( z y ) + 3 4 γ v x x z u y F 2 ( z y ) + 3 4 γ v x x y u y F 2 ( z y ) + 1 2 v t u y F 2 ( z y ) 1 2 α v u y z u x F 2 ( z y ) 1 2 α v u y y u x F 2 ( z y ) β v u z u x y F 2 ( z y ) 1 2 γ v x z u x y F 2 ( z y ) 1 2 γ u x y v x y F 2 ( z y ) + 1 2 γ v x u x y z F 2 ( z y ) + 1 2 γ v x u x y y F 2 ( z y ) 1 4 γ u y z v x x F 2 ( z y ) 1 4 γ u y y v x x F 2 ( z y ) + 1 4 γ v z u x x y F 2 ( z y ) + 1 4 γ v y u x x y F 2 ( z y ) 3 4 γ v u x x y z F 2 ( z y ) 3 4 γ v u x x y y F 2 ( z y ) 1 2 v u t y F 2 ( z y ) , Y 2 = 1 2 α u x x u y v F 2 ( z y ) + α u x u x z v F 2 ( z y ) + 1 2 α u x u x y v F 2 ( z y ) + β u x x u z v F 2 ( z y ) + γ u x x x z v F 2 ( z y ) + 3 4 γ u x x x y v F 2 ( z y ) + u t x v F 2 ( z y ) + 1 2 α u x u y v x F 2 ( z y ) + 1 4 γ v x u x x y F 2 ( z y ) 1 4 γ v x x u x y F 2 ( z y ) + 1 4 γ u y v x x x F 2 ( z y ) , Z 2 = 1 2 α u x x u y v F 2 ( z y ) 1 2 α u x u x y v F 2 ( z y ) β u x x u y v F 2 ( z y ) 1 4 γ u x x x y v F 2 ( z y ) + 1 2 α u x u y v x F 2 ( z y ) + 1 4 γ v x u x x y F 2 ( z y ) 1 4 γ v x x u x y F 2 ( z y ) + 1 4 γ u y v x x x F 2 ( z y ) ; T 3 = 1 2 u z v x F 3 ( z y ) 1 2 u x z v F 3 ( z y ) , X 3 = 1 2 α v z u x u z F 3 ( z y ) + β v x u z 2 F 3 ( z y ) + 1 2 α v y u x u z F 3 ( z y ) + β u y v x u z F 3 ( z y ) 1 2 α v u x z u z F 3 ( z y ) 1 2 α v u x y u z F 3 ( z y ) + β v u x y u z F 3 ( z y ) + 3 4 γ v x x z u z F 3 ( z y ) + 3 4 γ v x x y u z F 3 ( z y ) + 1 2 v t u z F 3 ( z y ) 1 2 α v u z z u x F 3 ( z y ) 1 2 α v u y z u x F 3 ( z y ) β v u y u x z F 3 ( z y ) 1 2 γ u x z v x z F 3 ( z y ) + 1 2 γ v x u x z z F 3 ( z y ) 1 2 γ u x z v x y F 3 ( z y ) + 1 2 γ v x u x y z F 3 ( z y ) 1 4 γ u z z v x x F 3 ( z y ) 1 4 γ u y z v x x F 3 ( z y ) + 1 4 γ v z u x x z F 3 ( z y ) + 1 4 γ v y u x x z F 3 ( z y ) 3 4 γ v u x x z z F 3 ( z y ) 3 4 γ v u x x y z F 3 ( z y ) 1 2 v u t z F 3 ( z y ) , Y 3 = 1 2 α u x x u z v F 3 ( z y ) 1 2 α u x u x z v F 3 ( z y ) β u x x u z v F 3 ( z y ) 1 4 γ u x x x z v F 3 ( z y ) + 1 2 α u x u z v x F 3 ( z y ) + 1 4 γ v x u x x z F 3 ( z y ) 1 4 γ v x x u x z F 3 ( z y ) + 1 4 γ u z v x x x F 3 ( z y ) , Z 3 = 1 2 α u x u x z v F 3 ( z y ) + α u x u x y v F 3 ( z y ) + 1 2 α u x x u z v F 3 ( z y ) + β u x x u y v F 3 ( z y ) + 3 4 γ u x x x z v F 3 ( z y ) + γ u x x x y v F 3 ( z y ) + u t x v F 3 ( z y ) + 1 2 α u x u z v x F 3 ( z y ) + 1 4 γ v x u x x z F 3 ( z y ) 1 4 γ v x x u x z F 3 ( z y ) + 1 4 γ u z v x x x F 3 ( z y ) ; T 4 = α u x u x z v F 4 ( z y ) + α u x u x y v F 4 ( z y ) + β u x x u z v F 4 ( z y ) + β u x x u y v F 4 ( z y ) + γ u x x x z v F 4 ( z y ) + γ u x x x y v F 4 ( z y ) + 1 2 u t x v F 4 ( z y ) + 1 2 u t v x F 4 ( z y ) , X 4 = 1 2 α v z u x u t F 4 ( z y ) + 1 2 α v y u x u t F 4 ( z y ) + β u z v x u t F 4 ( z y ) + β u y v x u t F 4 ( z y ) 1 2 α v u x z u t F 4 ( z y ) + β v u x z u t F 4 ( z y ) 1 2 α v u x y u t F 4 ( z y ) + β v u x y u t F 4 ( z y ) + 3 4 γ v x x z u t F 4 ( z y ) + 3 4 γ v x x y u t F 4 ( z y ) + 1 2 v t u t F 4 ( z y ) 1 2 α v u x u t z F 4 ( z y ) 1 4 γ v x x u t z F 4 ( z y ) 1 2 α v u x u t y F 4 ( z y ) 1 4 γ v x x u t y F 4 ( z y ) β v u z u t x F 4 ( z y ) β v u y u t x F 4 ( z y ) 1 2 γ v x z u t x F 4 ( z y ) 1 2 γ v x y u t x F 4 ( z y ) + 1 2 γ v x u t x z F 4 ( z y ) + 1 2 γ v x u t x y F 4 ( z y ) + 1 4 γ v z u t x x F 4 ( z y ) + 1 4 γ v y u t x x F 4 ( z y ) 3 4 γ v u t x x z F 4 ( z y ) 3 4 γ v u t x x y F 4 ( z y ) 1 2 v u t t F 4 ( z y ) , Y 4 = 1 2 α u t u x x v F 4 ( z y ) 1 2 α u x u t x v F 4 ( z y ) β u t u x x v F 4 ( z y ) 1 4 γ u t x x x v F 4 ( z y ) + 1 2 α u t u x v x F 4 ( z y ) + 1 4 γ u t v x x x F 4 ( z y ) 1 4 γ v x x u t x F 4 ( z y ) + 1 4 γ v x u t x x F 4 ( z y ) , Z 4 = 1 2 α u t u x x v F 4 ( z y ) 1 2 α u x u t x v F 4 ( z y ) β u t u x x v F 4 ( z y ) 1 4 γ u t x x x v F 4 ( z y ) + 1 2 α u t u x v x F 4 ( z y ) + 1 4 γ u t v x x x F 4 ( z y ) 1 4 γ v x x u t x F 4 ( z y ) + 1 4 γ v x u t x x F 4 ( z y ) ; T 5 = 1 2 v x F 5 ( t , z y ) , X 5 = 1 2 α u x v z F 5 ( t , z y ) 1 2 α u x v y F 5 ( t , z y ) + 1 2 α u x z v F 5 ( t , z y ) + 1 2 α u x y v F 5 ( t , z y ) β u z v x F 5 ( t , z y ) β u y v x F 5 ( t , z y ) β u x z v F 5 ( t , z y ) β u x y v F 5 ( t , z y ) 3 4 γ v x x z F 5 ( t , z y ) 3 4 γ v x x y F 5 ( t , z y ) + 1 2 v F t 5 ( t , z y ) 1 2 v t F 5 ( t , z y ) , Y 5 = 1 2 α u x v x F 5 ( t , z y ) 1 2 α u x x v F 5 ( t , z y ) + β u x x v F 5 ( t , z y ) 1 4 γ v x x x F 5 ( t , z y ) , Z 5 = 1 2 α u x v x F 5 ( t , z y ) 1 2 α u x x v F 5 ( t , z y ) + β u x x v F 5 ( t , z y ) 1 4 γ v x x x F 5 ( t , z y ) ; T 6 = 1 2 β u x v x F 6 ( t , z y ) 1 2 β u x x v F 6 ( t , z y ) 1 2 y v x F t 6 ( t , z y ) , X 6 = β 2 u z u x v x F 6 ( t , z y ) + β 2 u y u x v x F 6 ( t , z y ) + β 2 u x u x z v F 6 ( t , z y ) + β 2 u x u x y v F 6 ( t , z y ) + 1 2 α β v z u x 2 F 6 ( t , z y ) + 1 2 α β v y u x 2 F 6 ( t , z y ) 1 2 β u x v F t 6 β y u z v x F t 6 ( t , z y ) β y u y v x F t 6 ( t , z y ) β y v u x z F t 6 ( t , z y ) β y v u x y F t 6 ( t , z y ) 1 2 γ v x z u x x β F 6 ( t , z y ) 1 2 γ β v x y u x x F 6 ( t , z y ) 1 4 γ β u x z v x x F 6 ( t , z y ) 1 4 γ β u x y v x x F 6 ( t , z y ) + 1 2 γ β v x u x x z F 6 ( t , z y ) + 3 4 γ β u x v x x z F 6 ( t , z y ) + 1 2 γ β v x u x x y F 6 ( t , z y ) + 3 4 γ β u x v x x y F 6 ( t , z y ) + 1 4 γ β v z u x x x F 6 ( t , z y ) + 1 4 γ β v y u x x x F 6 ( t , z y ) + 1 4 γ β v u x x x z F 6 ( t , z y ) + 1 4 γ β v u x x x y F 6 ( t , z y ) + 1 2 β u x v t F 6 ( t , z y ) + 1 2 β v u t x F 6 ( t , z y ) + 1 2 y v F t t 6 ( t , z y ) + 1 2 α u x v F t 6 ( t , z y ) 1 2 α y v z u x F t 6 ( t , z y ) 1 2 α y v y u x F t 6 ( t , z y ) + 1 2 α y u x z v F t 6 ( t , z y ) + 1 2 α y u x y v F t 6 ( t , z y ) + 1 4 γ v x x F t 6 ( t , z y ) 3 4 γ y v x x z F t 6 ( t , z y ) 3 4 γ y v x x y F t 6 ( t , z y ) 1 2 y v t F t 6 ( t , z y ) + 1 2 y v 6 F t t 5 ( t , z y ) , Y 6 = 1 2 α β u x 2 v x F 6 ( t , z y ) β 2 u x u x x v F 6 ( t , z y ) 1 4 β γ u x x v x x F 6 ( t , z y ) + 1 4 β γ u x x x v x F 6 ( t , z y ) + 1 4 β γ u x v x x x F 6 ( t , z y ) 1 4 β γ u x x x x v F 6 ( t , z y ) 1 2 α y u x v x F t 6 ( t , z y ) 1 2 α y u x x v F t 6 ( t , z y ) + β y u x x v F t 6 ( t , z y ) 1 4 γ y v x x x F t 6 ( t , z y ) , Z 6 = 1 2 α β u x 2 v x F 6 ( t , z y ) β 2 u x u x x v F 6 ( t , z y ) 1 4 β γ u x x v x x F 6 ( t , z y ) + 1 4 β γ u x x x v x F 6 ( t , z y ) + 1 4 β γ u x v x x x F 6 ( t , z y ) 1 4 β γ u x x x x v F 6 ( t , z y ) 1 2 α y u x v x F t 6 ( t , z y ) 1 2 α y u x x v F t 6 ( t , z y ) + β y u x x v F t 6 ( t , z y ) 1 4 γ y v x x x F t 6 ( t , z y ) ; T 7 = 2 α t u x u x z v F 7 ( z y ) + 2 α t u x u x y v F 7 ( z y ) + 2 β t u x x u z v F 7 ( z y ) + 2 β t u x x u y v F 7 ( z y ) + 2 γ t u x x x z v F 7 ( z y ) + 2 γ t u x x x y v F 7 ( z y ) u x v F 7 ( z y ) + 1 2 v x u F 7 ( z y ) 1 2 x u x x v F 7 ( z y ) + t u t x v F 7 ( z y ) + t u t v x F 7 ( z y ) + 1 2 x u x v x F 7 ( z y ) , X 7 = 1 2 α x v z u x 2 F 7 ( z y ) + 1 2 α x v y u x 2 F 7 ( z y ) 1 2 α v u z u x F 7 ( z y ) 2 β v u z u x F 7 ( z y ) + 1 2 α u v z u x F 7 ( z y ) 1 2 α v u y u x F 7 ( z y ) 2 β v u y u x F 7 ( z y ) + 1 2 α u v y u x F 7 ( z y ) + β x u z v x u x F 7 ( z y ) + β x u y v x u x F 7 ( z y ) + β x v u x z u x F 7 ( z y ) γ v x z u x F 7 ( z y ) + β x v u x y u x F 7 ( z y ) γ v x y u x F 7 ( z y ) + 3 4 γ x v x x z u x F 7 ( z y ) + 3 4 γ x v x x y u x F 7 ( z y ) + α t v z u t u x F 7 ( z y ) + α t v y u t u x F 7 ( z y ) + 1 2 x v t u x F 7 ( z y ) α t v u t z u x F 7 ( z y ) α t v u t y u x F 7 ( z y ) + β u u z v x F 7 ( z y ) + β u u y v x F 7 ( z y ) 1 2 α u v u x z F 7 ( z y ) + β u v u x z F 7 ( z y ) + γ v x u x z F 7 ( z y ) 1 2 α u v u x y F 7 ( z y ) + β u v u x y F 7 ( z y ) + γ v x u x y F 7 ( z y ) + 3 4 γ v z u x x F 7 ( z y ) + 3 4 γ v y u x x F 7 ( z y ) 1 2 γ x v x z u x x F 7 ( z y ) 1 2 γ x v x y u x x F 7 ( z y ) 1 4 γ u z v x x F 7 ( z y ) 1 4 γ u y v x x F 7 ( z y ) 1 4 γ x u x z v x x F 7 ( z y ) 1 4 γ x u x y v x x F 7 ( z y ) 9 4 γ v u x x z F 7 ( z y ) + 1 2 γ x v x u x x z F 7 ( z y ) + 3 4 γ u v x x z F 7 ( z y ) 9 4 γ v u x x y F 7 ( z y ) + 1 2 γ x v x u x x y F 7 ( z y ) + 3 4 γ u v x x y F 7 ( z y ) + 1 4 γ x v z u x x x F 7 ( z y ) + 1 4 γ x v y u x x x F 7 ( z y ) + 1 4 γ x v u x x x z F 7 ( z y ) + 1 4 γ x v u x x x y F 7 ( z y ) 3 2 v u t F 7 ( z y ) + 2 β t u z v x u t F 7 ( z y ) + 2 β t u y v x u t F 7 ( z y ) α t v u x z u t F 7 ( z y ) + 2 β t v u x z u t F 7 ( z y ) α t v u x y u t F 7 ( z y ) + 2 β t v u x y u t F 7 ( z y ) + 3 2 γ t v x x z u t F 7 ( z y ) + 3 2 γ t v x x y u t F 7 ( z y ) + 1 2 u v t F 7 ( z y ) + t u t v t F 7 ( z y ) 1 2 γ t v x x u t z F 7 ( z y ) 1 2 γ t v x x u t y F 7 ( z y ) + 1 2 x v u t x F 7 ( z y ) 2 β t v u z u t x F 7 ( z y ) 2 β t v u y u t x F 7 ( z y ) γ t v x z u t x F 7 ( z y ) γ t v x y u t x F 7 ( z y ) + γ t v x u t x z F 7 ( z y ) + γ t v x u t x y F 7 ( z y ) + 1 2 γ t v z u t x x F 7 ( z y ) + 1 2 γ t v y u t x x F 7 ( z y ) 3 2 γ t v u t x x z F 7 ( z y ) 3 2 γ t v u t x x y F 7 ( z y ) t v u t t F 7 ( z y ) , Y 7 = 1 2 α u x v x u F 7 ( z y ) α u x 2 v F 7 ( z y ) α t u x u t x v F 7 ( z y ) + 1 2 α u x x u v F 7 ( z y ) + α t u t u x x v F 7 ( z y ) β x u x x u x v F 7 ( z y ) β u x x u v F 7 ( z y ) 2 β t u t u x x v F 7 ( z y ) γ u x x x v F 7 ( z y ) + 1 4 γ v x x x u F 7 ( z y ) 1 4 γ x u x x x x v F 7 ( z y ) 1 2 γ t u t x x x v F 7 ( z y ) + α t u t u x v x F 7 ( z y ) + 1 2 γ t u t v x x x F 7 ( z y ) 1 2 γ t v x x u t x F 7 ( z y ) + 1 2 γ t v x u t x x F 7 ( z y ) + 1 2 α x u x 2 v x F 7 ( z y ) 1 2 γ u x v x x F 7 ( z y ) + 1 4 γ x u x v x x x F 7 ( z y ) + 3 4 γ u x x v x F 7 ( z y ) 1 4 γ x u x x v x x F 7 ( z y ) + 1 4 γ x u x x x v x F 7 ( z y ) , Z 7 = 1 2 α u x v x u F 7 ( z y ) α u x 2 v F 7 ( z y ) α t u x u t x v F 7 ( z y ) + 1 2 α u x x u v F 7 ( z y ) + α t u t u x x v F 7 ( z y ) β x u x x u x v F 7 ( z y ) β u x x u v F 7 ( z y ) 2 β t u t u x x v F 7 ( z y ) γ u x x x v F 7 ( z y ) + 1 4 γ v x x x u F 7 ( z y ) 1 4 γ x u x x x x v F 7 ( z y ) 1 2 γ t u t x x x v F 7 ( z y ) + α t u t u x v x F 7 ( z y ) + 1 2 γ t u t v x x x F 7 ( z y ) 1 2 γ t v x x u t x F 7 ( z y ) + 1 2 γ t v x u t x x F 7 ( z y ) + 1 2 α x u x 2 v x F 7 ( z y ) 1 2 γ u x v x x F 7 ( z y ) + 1 4 γ x u x v x x x F 7 ( z y ) + 3 4 γ u x x v x F 7 ( z y ) 1 4 γ x u x x v x x F 7 ( z y ) + 1 4 γ x u x x x v x F 7 ( z y ) ; T 8 = α t u x u x z v F 8 ( z y ) + α t u x u x y v F 8 ( z y ) + β t u x x u z v F 8 ( z y ) + β t u x x u y v F 8 ( z y ) + γ t u x x x z v F 8 ( z y ) + γ t u x x x y v F 8 ( z y ) 1 2 y u x z v F 8 ( z y ) 1 2 y u x y v F 8 ( z y ) + 1 2 t u t x v F 8 ( z y ) + 1 2 t u t v x F 8 ( z y ) + 1 2 y u z v x F 8 ( z y ) + 1 2 y u y v x F 8 ( z y ) , X 8 = β y v x u z 2 F 8 ( z y ) 1 2 α v u x u z F 8 ( z y ) + 1 2 α y v z u x u z F 8 ( z y ) + 1 2 α y v y u x u z F 8 ( z y ) + 2 β y u y v x u z F 8 ( z y ) 1 2 α y v u x z u z F 8 ( z y ) 1 2 α y v u x y u z F 8 ( z y ) 1 4 γ v x x u z F 8 ( z y ) + 3 4 γ y v x x z u z F 8 ( z y ) + 3 4 γ y v x x y u z F 8 ( z y ) + β t v x u t u z F 8 ( z y ) + 1 2 y v t u z F 8 ( z y ) β t v u t x u z F 8 ( z y ) 1 2 α y v u z z u x F 8 ( z y ) 1 2 α v u y u x F 8 ( z y ) + 1 2 α y v z u y u x F 8 ( z y ) + 1 2 α y u y v y u x F 8 ( z y ) α y v u y z u x F 8 ( z y ) 1 2 α y v u y y u x F 8 ( z y ) + β y u y 2 v x F 8 ( z y ) 1 2 α y v u y u x z F 8 ( z y ) + 1 2 γ v x u x z F 8 ( z y ) 1 2 γ y u x z v x z F 8 ( z y ) + 1 2 γ y v x u x z z F 8 ( z y ) 1 2 α y v u y u x y F 8 ( z y ) + 1 2 γ v x u x y F 8 ( z y ) 1 2 γ y v x z u x y F 8 ( z y ) 1 2 γ y u x z v x y F 8 ( z y ) 1 2 γ y u x y v x y F 8 ( z y ) + γ y v x u x y z F 8 ( z y ) + 1 2 γ y v x u x y y F 8 ( z y ) 1 4 γ y u z z v x x F 8 ( z y ) 1 4 γ u y v x x F 8 ( z y ) 1 2 γ y u y z v x x F 8 ( z y ) 1 4 γ y u y y v x x F 8 ( z y ) 3 4 γ v u x x z F 8 ( z y ) + 1 4 γ y v z u x x z F 8 ( z y ) + 1 4 γ y v y u x x z F 8 ( z y ) + 3 4 γ y u y v x x z F 8 ( z y ) 3 4 γ y v u x x z z F 8 ( z y ) 3 4 γ v u x x y F 8 ( z y ) + 1 4 γ y v z u x x y F 8 ( z y ) + 1 4 γ y v y u x x y F 8 ( z y ) + 3 4 γ y u y v x x y F 8 ( z y ) 3 2 γ y v u x x y z F 8 ( z y ) 3 4 γ y v u x x y y F 8 ( z y ) 1 2 v u t F 8 ( z y ) + 1 2 α t v z u x u t F 8 ( z y ) + 1 2 α t v y u x u t F 8 ( z y ) + β t u y v x u t F 8 ( z y ) 1 2 α t v u x z u t F 8 ( z y ) + β t v u x z u t F 8 ( z y ) 1 2 α t v u x y u t F 8 ( z y ) + β t v u x y u t F 8 ( z y ) + 3 4 γ t v x x z u t F 8 ( z y ) + 3 4 γ t v x x y u t F 8 ( z y ) + 1 2 y u y v t F 8 ( z y ) + 1 2 t u t v t F 8 ( z y ) 1 2 y v u t z F 8 ( z y ) 1 2 α t v u x u t z F 8 ( z y ) 1 4 γ t v x x u t z F 8 ( z y ) 1 2 y v u t y F 8 ( z y ) 1 2 α t v u x u t y F 8 ( z y ) 1 4 γ t v x x u t y F 8 ( z y ) β t v u y u t x F 8 ( z y ) 1 2 γ t v x z u t x F 8 ( z y ) 1 2 γ t v x y u t x F 8 ( z y ) + 1 2 γ t v x u t x z F 8 ( z y ) + 1 2 γ t v x u t x y F 8 ( z y ) + 1 4 γ t v z u t x x F 8 ( z y ) + 1 4 γ t v y u t x x F 8 ( z y ) 3 4 γ t v u t x x z F 8 ( z y ) 3 4 γ t v u t x x y F 8 ( z y ) 1 2 t v u t t F 8 ( z y ) , Y 8 = 1 2 α y u z u x v x F 8 ( z y ) + 1 2 α y u y u x v x F 8 ( z y ) + 1 4 γ y u x x z v x F 8 ( z y ) + 1 4 γ y u x x y v x F 8 ( z y ) + 1 2 α t u x u t v x F 8 ( z y ) + 1 4 γ t u t x x v x F 8 ( z y ) + 1 2 α y v u x u x z F 8 ( z y ) + 1 2 α y v u x u x y F 8 ( z y ) + 1 2 α y v u z u x x F 8 ( z y ) + 1 2 α y v u y u x x F 8 ( z y ) 1 4 γ y u x z v x x F 8 ( z y ) 1 4 γ y u x y v x x F 8 ( z y ) + 1 4 γ y u z v x x x F 8 ( z y ) + 1 4 γ y u y v x x x F 8 ( z y ) + 3 4 γ y v u x x x z F 8 ( z y ) + 3 4 γ y v u x x x y F 8 ( z y ) + 1 2 α t v u x x u t F 8 ( z y ) β t v u x x u t F 8 ( z y ) + 1 4 γ t v x x x u t F 8 ( z y ) + y v u t x F 8 ( z y ) 1 2 α t v u x u t x F 8 ( z y ) 1 4 γ t v x x u t x F 8 ( z y ) 1 4 γ t v u t x x x F 8 ( z y ) , Z 8 = 1 2 α y u z u x v x F 8 ( z y ) + 1 2 α y u y u x v x F 8 ( z y ) + 1 4 γ y u x x z v x F 8 ( z y ) + 1 4 γ y u x x y v x F 8 ( z y ) + 1 2 α t u x u t v x F 8 ( z y ) + 1 4 γ t u t x x v x F 8 ( z y ) + 1 2 α y v u x u x z F 8 ( z y ) + 1 2 α y v u x u x y F 8 ( z y ) + 1 2 α y v u z u x x F 8 ( z y ) + 1 2 α y v u y u x x F 8 ( z y ) 1 4 γ y u x z v x x F 8 ( z y ) 1 4 γ y u x y v x x F 8 ( z y ) + 1 4 γ y u z v x x x F 8 ( z y ) + 1 4 γ y u y v x x x F 8 ( z y ) + 3 4 γ y v u x x x z F 8 ( z y ) + 3 4 γ y v u x x x y F 8 ( z y ) + 1 2 α t v u x x u t F 8 ( z y ) β t v u x x u t F 8 ( z y ) + 1 4 γ t v x x x u t F 8 ( z y ) + y v u t x F 8 ( z y ) 1 2 α t v u x u t x F 8 ( z y ) 1 4 γ t v x x u t x F 8 ( z y ) 1 4 γ t v u t x x x F 8 ( z y ) ; T 9 = 1 2 α t u x z v F 9 ( z y ) 1 2 α t u x y v F 9 ( z y ) + 1 2 v F 9 ( z y ) + 1 2 α t u z v x F 9 ( z y ) + 1 2 α t u y v x F 9 ( z y ) 1 2 x v x F 9 ( z y ) , X 9 = 1 2 α 2 t u z v z u x F 9 ( z y ) 1 2 α 2 t v u z z u x F 9 ( z y ) + 1 2 α 2 t v z u y u x F 9 ( z y ) + 1 2 α 2 t u z v y u x F 9 ( z y ) + 1 2 α 2 t u y v y u x F 9 ( z y ) α 2 t v u y z u x F 9 ( z y ) α 2 1 2 t v u y y u x F 9 ( z y ) 1 2 α 2 t v u z u x z F 9 ( z y ) 1 2 α 2 t v u y u x z F 9 ( z y ) 1 2 α 2 t v u z u x y F 9 ( z y ) 1 2 α 2 t v u y u x y F 9 ( z y ) 1 2 α v u z F 9 ( z y ) 1 2 α v u y F 9 ( z y ) 1 2 α x v z u x F 9 ( z y ) 1 2 α x v y u x F 9 ( z y ) + α β t u z 2 v x F 9 ( z y ) + α β t u y 2 v x F 9 ( z y ) + 2 α β t u z u y v x F 9 ( z y ) + 1 2 α x v u x z F 9 ( z y ) 1 2 α γ t u x z v x z F 9 ( z y ) + 1 2 α γ t v x u x z z F 9 ( z y ) + 1 2 α x v u x y F 9 ( z y ) 1 2 α γ t v x z u x y F 9 ( z y ) 1 2 α γ t u x z v x y F 9 ( z y ) 1 2 α γ t u x y v x y F 9 ( z y ) + α γ t v x u x y z F 9 ( z y ) + 1 2 α γ t v x u x y y F 9 ( z y ) 1 4 α γ t u z z v x x F 9 ( z y ) 1 2 α γ t u y z v x x F 9 ( z y ) 1 4 α γ t u y y v x x F 9 ( z y ) + 1 4 α γ t v z u x x z F 9 ( z y ) + 1 4 α γ t v y u x x z F 9 ( z y ) + 3 4 α γ t u z v x x z F 9 ( z y ) + 3 4 α γ t u y v x x z F 9 ( z y ) 3 4 α γ t v u x x z z F 9 ( z y ) + 1 4 α γ t v z u x x y F 9 ( z y ) + 1 4 α γ t v y u x x y F 9 ( z y ) + 3 4 α γ t u z v x x y F 9 ( z y ) + 3 4 α γ t u y v x x y F 9 ( z y ) 3 2 α γ t v u x x y z F 9 ( z y ) 3 4 α γ t v u x x y y F 9 ( z y ) + 1 2 α t u z v t F 9 ( z y ) + 1 2 α t u y v t F 9 ( z y ) 1 2 α t v u t z F 9 ( z y ) 1 2 α t v u t y F 9 ( z y ) + β v u z F 9 ( z y ) + β v u y F 9 ( z y ) β x u z v x F 9 ( z y ) β x u y v x F 9 ( z y ) β x v u x z F 9 ( z y ) + 1 2 γ v x z F 9 ( z y ) β x v u x y F 9 ( z y ) + 1 2 γ v x y F 9 ( z y ) 3 4 γ x v x x z F 9 ( z y ) 3 4 γ x v x x y F 9 ( z y ) 1 2 x v t F 9 ( z y ) , Y 9 = 1 2 α 2 t u x u x z v F 9 ( z y ) + 1 2 α 2 t u x u x y v F 9 ( z y ) + 1 2 α 2 t u x x u z v F 9 ( z y ) + 1 2 α 2 t u x x u y v F 9 ( z y ) + 3 4 α γ t u x x x z v F 9 ( z y ) + 3 4 α γ t u x x x y v F 9 ( z y ) + 1 2 α u x v F 9 ( z y ) 1 2 α x u x x v F 9 ( z y ) + α t u t x v F 9 ( z y ) + β x u x x v F 9 ( z y ) + 1 2 α 2 t u x u z v x F 9 ( z y ) + 1 2 α 2 t u x u y v x F 9 ( z y ) 1 4 α γ t v x x u x z F 9 ( z y ) 1 4 α γ t v x x u x y F 9 ( z y ) + 1 4 α γ t v x u x x z F 9 ( z y ) + 1 4 α γ t v x u x x y F 9 ( z y ) + 1 4 α γ t u z v x x x F 9 ( z y ) + 1 4 α γ t u y v x x x F 9 ( z y ) 1 2 α x u x v x F 9 ( z y ) + 1 4 γ v x x F 9 ( z y ) 1 4 γ x v x x x F 9 ( z y ) , Z 9 = 1 2 α 2 t u x u x z v F 9 ( z y ) + 1 2 α 2 t u x u x y v F 9 ( z y ) + 1 2 α 2 t u x x u z v F 9 ( z y ) + 1 2 α 2 t u x x u y v F 9 ( z y ) + 3 4 α γ t u x x x z v F 9 ( z y ) + 3 4 α γ t u x x x y v F 9 ( z y ) + 1 2 α u x v F 9 ( z y ) 1 2 α x u x x v F 9 ( z y ) + α t u t x v F 9 ( z y ) + β x u x x v F 9 ( z y ) + 1 2 α 2 t u x u z v x F 9 ( z y ) + 1 2 α 2 t u x u y v x F 9 ( z y ) 1 4 α γ t v x x u x z F 9 ( z y ) 1 4 α γ t v x x u x y F 9 ( z y ) + 1 4 α γ t v x u x x z F 9 ( z y ) + 1 4 α γ t v x u x x y F 9 ( z y ) + 1 4 α γ t u z v x x x F 9 ( z y ) + 1 4 α γ t u y v x x x F 9 ( z y ) 1 2 α x u x v x F 9 ( z y ) + 1 4 γ v x x F 9 ( z y ) 1 4 γ x v x x x F 9 ( z y ) .
Remark 2.
We notice that the conserved vectors obtained using multiplier as well as Ibragimov’s theorem contain arbitrary functions. Moreover, this approach also secures nine conserved quantities corresponding to the nine Lie symmetries obtained. All these attest to the fact that the 3D-gBSe (2) possesses conserved vectors that are infinitely many in number. Theoretically speaking, the obtained conserved vectors could be used to find exact solutions of the3D-gBSe (2). However, this task is not so simple and could be possibly done in the future work.

4. Discussion of the Graphical Representations of Solutions Obtained for 3D-gBSe (2)

The dynamics of solutions in the model of the generalized (3+1)-dimensional breaking soliton equation are described graphically. In a bid to gain better knowledge regarding the physical properties of the secured results, the derived solutions were depicted by choosing different values of parameters. For instance, the 3D-plot, density plot and 2D-plot of the periodic soliton solution (20) are represented with different values of parameters in Figure 1, where we assign values to the constants; in Figure 1, we let ω = 0.2 , ν = 0.6 , r 1 = 100 , r 2 = 50.05 , r 3 = 20.05 , Q = 10 , R = 70 and C 3 = 1 ; in Figure 1 we made ω = 0.2 , ν = 0.6 , r 1 = 90 , r 2 = 40.05 , r 3 = 0.05 , Q = 10 , R = 70 and C 3 = 1 ; and in Figure 3 we assumed ω = 0.2 , ν = 0.6 , r 1 = 100 , r 2 = 50.05 , r 3 = 80 , Q = 10 , R = 70 and C 3 = 1 ,. In addition, we varied the values of t and z such that they were allocated values t = 10 , z = 1 , t = 11 , and z = 0 ; and t = 25 and z = 9 for the Figures respectively. It is observable that the wave deflection changes as time t reduces from 10 to 25 .
Furthermore, the singular soliton solutions (24) obtained in the form of hyperbolic and trigonometric functions were plotted. 3D, density and 2D plots of (24) were done by setting positive values to the constants; in Figure 4 we supposed λ = 0.08 , ω = 2.02 , ν = 2.02 , α = 30 , μ = 10.125 , β = 20.2 , C 1 = 0.5 and C 2 = 2.3 ; in Figure 5 we let λ = 0.08 , ω = 2.02 , ν = 0.02 , α = 30 , μ = 10.125 , β = 20.2 , C 1 = 1 and C 2 = 2 ; and in Figure 6 we assumed λ = 0.0008 , ω = 0.02 , ν = 0.02 , α = 300 , μ = 0.125 , β = 200.2 , C 1 = 1 and C 2 = 0.8 . More so, we randomized t and z in which the figures respectively have t = 1 , z = 0 , t = 2 and z = 3.5 ; and t = 0.2 and z = 20 . Figure 7 shows the 3D, density and 2D-plots of solution (25) with a singularity when t = 2.1 and z = 3 ; Figure 8 reflects the same with a singularity when t = 6 and z = 25 ; and Figure 9 exhibits (25) with a singularity at t = 4 , z = 30 . Other parameters involved were allotted values: λ = 0.02 , ω = 2.02 , ν = 0.02 , μ = 0.005 , α = 0.025 , β = 1 , C 1 = 0 and C 2 = 2 ; in Figure 7, λ = 0.02 , ω = 2.02 , ν = 0.02 and μ = 0.005 ; in Figure 8, α = 0.025 , β = 1 , C 1 = 0.1 and C 2 = 2 ; and in Figure 9, λ = 0.02 , ω = 2.02 , ν = 0.02 , μ = 0.005 , α = 0.025 , β = 1 , C 1 = 0.1 and C 2 = 2 .
In addition, we made graphical demonstrations of the soliton solution (27) with singularities, by using 3D-plots, density plots and 2D-plots. In this case, we allocated values to the constants; for instance, in Figure 10 we let ω = 0.48 , ν = 0 , α = 0.2 , β = 4 , B 0 = 0 , t = 1 and x = 2 ; in Figure 11 we put ω = 0.51 , ν = 0 , α = 0.5 , β = 4 , B 0 = 0 , t = 0 and x = 3 . In Figure 12 we assumed ω = 0.51 , ν = 0 , α = 0.5 , β = 4 , B 0 = 0 , t = 1 and x = 4 .
Moreover, Figure 13 displays 3D, density and 2D-plots pictorial representation of soliton solution (28) with singularities when ω = 0.48 , ν = 0 , α = 0.2 , β = 4 , B 0 = 0 , t = 1 and x = 4 . Figure 14 also reveals the 3D-plot, density and 2D-plots graphic display of (28) with singularities when ω = 0.48 , ν = 0 , α = 0.2 , β = 4 , B 0 = 0 , t = 0 and x = 4 . Besides, with ω = 0.6 , ν = 0 , α = 0.2 , β = 4 , B 0 = 0 , t = 6 and x = 7 , we show 3D, density and 2D-plots of Figure 15.
Finally, we show the solitary wave profile of series solution (34) with 3D-plot, density plot and 2D-plot in Figure 16 when ω = 1.6 , ν = 0.3 , c 0 = 0 , c 1 = 0.2 , c 2 = 7.825 , c 3 = 0.125 , P = 1.8 , Q = 83 , R = 25 , t = 1.6 and z = 1.1 . Figure 17 exhibits the 3D, density and 2D plots of solution (34) when ω = 1.6 , ν = 0.4 , c 0 = 0 , c 1 = 0.2 , c 2 = 10.825 , c 3 = 0.125 , P = 10.8 , Q = 83 , R = 25 , t = 1 and z = 8 . When ω = 0.1 , ν = 0.4 , c 0 = 0 , c 1 = 0.2 , c 2 = 8.825 , c 3 = 0.125 , P = 10.8 , Q = 83 , R = 25 , t = 5 and z = 4 Figure 18 shows 3D-plot, density and 2D-plots of series solution (34).
We can see that the changes in the parameters cause clear and obvious effects on the dynamics of the solitons and the series solutions. By making a comparison of the results secured in this work and those of the previous studies, it was revealed that the solutions gained here are completely different.

5. Conclusions

In this article, we investigated a (3+1)-dimensional generalized breaking soliton Equation (2). We constructed exact solutions for Equation (2) via the implementation of Lie symmetry reductions, direct integration, ( G / G ) -expansion and power series solution methods. Solutions obtained were in the form of trigonometric, hyperbolic functions and Jacobi elliptic functions. Moreover, the graphical representations of the solutions were shown with various values of the parameters and t, x and z to give a view of the wave dynamical appearances of the solutions. Lastly, we constructed the conserved vectors of Equation (2) through the engagement of the general multiplier method and Ibragimov’s theorem. The results of this paper could be of interest to scientists working in the fields of engineering and nonlinear sciences.

Author Contributions

Conceptualization, O.D.A. and C.M.K.; methodology, C.M.K.; software, O.D.A.; validation, C.M.K.; writing—original draft preparation, O.D.A.; writing—review and editing, C.M.K.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank North-West University for its unalloyed sustenance and sincerely 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. Wazwaz, A.M. Multiple optical kink solutions for new Painlevé integrable (3+1)-dimensional sine-Gordon equations with constant and time-dependent coefficients. Optik 2020, in press. [Google Scholar]
  2. Wazwaz, A.M. Painlevé analysis for new (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equations with constant and time-dependent equations. Int. Numer. Methods Heat Fluid Flow 2019, 30, 996–1008. [Google Scholar] [CrossRef]
  3. Motsepa, T.; Khalique, C.M. Closed-form solutions and conserved vectors of the (3+1)-dimensional negative-order KdV equation. Adv. Math. Model. Appl. 2020, 5, 7–18. [Google Scholar]
  4. Wazwaz, A.M.; Xu, G.Q. Integrability aspects and localized wave solutions for a new (4+1)-dimensional Boiti-Leon-Manna-Pempinelli equation. Nonlinear Dyn. 2019, 98, 1379–1390. [Google Scholar]
  5. Gandarias, M.L.; Rosa, M.; 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]
  6. 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]
  7. 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]
  8. Zabusky, N.J.; Kruskal, M.D. Interaction of solitons in a collisionless plasma and recurrence of initial states. Phys. Rev. Lett. 1965, 15, 240–243. [Google Scholar] [CrossRef] [Green Version]
  9. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  10. Weiss, J.; Tabor, M.; Carnevale, G. The Painleve property and a partial differential equations with an essential singularity. Phys. Lett. A 1985, 109, 205–208. [Google Scholar]
  11. Zheng, C.L.; Fang, J.P. New exact solutions and fractional patterns of generalized Broer-Kaup system via a mapping approach. Chaos Soliton Fract. 2006, 27, 1321–1327. [Google Scholar] [CrossRef]
  12. Wen, X.Y. Construction of new exact rational form non-travelling wave solutions to the (2+1)-dimensional generalized BroerKaup system. Appl. Math. Comput. 2010, 217, 1367–1375. [Google Scholar] [CrossRef]
  13. Gu, C.H. Soliton Theory and Its Application; Zhejiang Science and Technology Press: Hangzhou, China, 1990. [Google Scholar]
  14. Zeng, X.; Wang, D.S. A generalized extended rational expansion method and its application to (1+1)-dimensional dispersive long wave equation. Appl. Math. Comput. 2009, 212, 296–304. [Google Scholar] [CrossRef]
  15. Zhou, Y.; Wang, M.; Wang, Y. Periodic wave solutions to a coupled KdV equations with variable coefficients. Phys. Lett. A 2003, 308, 31–36. [Google Scholar] [CrossRef]
  16. 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]
  17. Kudryashov, N.A.; Loguinova, N.B. Extended simplest equation method for nonlinear differential equations. Appl. Math. Comput. 2008, 205, 396–402. [Google Scholar] [CrossRef]
  18. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  19. Olver, P.J. Applications of Lie Groups to Differential Equations, 2nd ed.; Springer: Berlin, Germany, 1993. [Google Scholar]
  20. Ibragimov, N.H. Elementary Lie Group Analysis and Ordinary Differential Equations; John Wiley & Sons: Chichester, NY, USA, 1999. [Google Scholar]
  21. 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. Ser. S 2018, 11, 777–790. [Google Scholar] [CrossRef]
  22. 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]
  23. Matveev, V.B.; Salle, M.A. Darboux Transformations and Solitons; Springer: New York, NY, USA, 1991. [Google Scholar]
  24. Chen, Y.; Yan, Z. New exact solutions of (2+1)-dimensional Gardner equation via the new sine-Gordon equation expansion method. Chaos Solitons Fractals 2005, 26, 399–406. [Google Scholar] [CrossRef]
  25. 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]
  26. Chargui, Y.; Trabelsi, A.; Chetouani, L. Exact solution of the (1+1)-dimensional Dirac equation with vector and scalar linear potentials in the presence of a minimal length. Phys. Lett. A 2010, 374, 531–534. [Google Scholar]
  27. 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]
  28. Wazwaz, A.M. Multiple-soliton solutions for extended (3+1)-dimensional Jimbo-Miwa equations. Appl. Math. Lett. 2017, 64, 21–26. [Google Scholar] [CrossRef]
  29. Kaur, L.; Wazwaz, A.M. Dynamical analysis of lump solutions for (3+1)-dimensional generalized KP-Boussinesq equation and its dimensionally reduced equations. Phys. Scr. 2018, 93, 075203. [Google Scholar] [CrossRef]
  30. Ali, M.R.; Ma, W.X. New exact solutions of nonlinear (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equation. Adv. Math. Phys. 2019, 2019, 7. [Google Scholar] [CrossRef] [Green Version]
  31. Verma, P.; Kaur, L. Extended exp ( ϕ ( ξ ) ) -expansion method for generalized breaking soliton equation. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2020; Volume 2214. [Google Scholar]
  32. Yildirim, Y.; Yasar, E. A (2+1)-dimensional breaking soliton equation: Solutions and conservation laws. Chaos Solitons Fract. 2018, 107, 146–155. [Google Scholar] [CrossRef]
  33. Mei, J.; Zhang, H. New types of exact solutions for a breaking soliton equation. Chaos Solitons Fract. 2004, 20, 771–777. [Google Scholar] [CrossRef]
  34. Yan-Ze, P.; Krishnan, E.V. Two classes of new exact solutions to (2+1)-dimensional breaking soliton equation. Commun. Theor. Phys. 2005, 44, 807. [Google Scholar]
  35. Wazwaz, A.M. Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders. Pramana 2016, 87, 68. [Google Scholar] [CrossRef]
  36. Wazwaz, A.M. The (2+1) and (3+1)-dimensional CBS equations: Multiple soliton solutions and multiple singular soliton solutions. Z. Naturforsch A 2010, 65, 173–181. [Google Scholar] [CrossRef]
  37. Darvishi, M.T.; Najafi, M.; Najafi, M. New exact solutions for the (3+1)-dimensional breaking soliton equation. Int. J. Math. 2010, 6, 134–137. [Google Scholar]
  38. Al-Amr, M.O. Exact solutions of the generalized (2+1)-dimensional nonlinear evolution equations via the modified simple equation method. Comput. Math. Appl. 2015, 69, 390–397. [Google Scholar] [CrossRef]
  39. Liu, X.Z.; Yu, J.; Lou, Z.M. Residual symmetry, CRE integrability and interaction solutions of the (3+1)-dimensional breaking soliton equation. Phys. Scr. 2018, 93, 085201. [Google Scholar] [CrossRef]
  40. 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]
  41. Leveque, R.J. Numerical Methods for Conservation Laws, 2nd ed.; Birkhäuser-Verlag: Basel, Switzerland, 1992. [Google Scholar]
  42. 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]
  43. Ibragimov, N.H. Integrating factors, adjoint equations and Lagrangians. J. Math. Anal. Appl. 2006, 318, 742–757. [Google Scholar] [CrossRef] [Green Version]
  44. Ibragimov, N.H. A new conservation theorem. J. Math. Anal. Appl. 2007, 333, 311–328. [Google Scholar] [CrossRef] [Green Version]
  45. 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]
  46. Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1972. [Google Scholar]
  47. Kudryashov, N.A. Analytical Theory of Nonlinear Differential Equations; Institute of Computer Investigations: Moskow-Igevsk, Russia, 2004. [Google Scholar]
  48. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products, 7th ed.; Academic Press: New York, NY, USA, 2007. [Google Scholar]
  49. Wang, G.; Liu, X.; Zhang, Y. Symmetry reduction, exact solutions and conservation laws of a new fifth-order nonlinear integrable equation. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 2313–2320. [Google Scholar] [CrossRef]
  50. Liu, H.; Li, J. Lie symmetry analysis and exact solutions for the short pulse equation. Nonlinear Anal. 2009, 21, 26–33. [Google Scholar] [CrossRef]
  51. Liu, H.; Li, J.; Zhang, Q.X. Lie symmetry analysis and exact solutions for general Burger’s equation. J. Comput. Appl. Math. 2009, 228, 1–9. [Google Scholar] [CrossRef] [Green Version]
  52. Bruzón, M.S.; Gandarias, M.L. Travelling wave solutions of the K(m, n) equation with generalized evolution. Math. Meth. Appl. Sci. 2018, 41, 5851–5857. [Google Scholar] [CrossRef]
  53. Cheviakov, A.F. Computation of fluxes of conservation laws. J. Eng. Math. 2010, 66, 153–173. [Google Scholar] [CrossRef] [Green Version]
  54. 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; World Scientific Publishing Co Pte Ltd.: Singapore, 2018; pp. 253–263. [Google Scholar]
  55. Khalique, C.M.; Moleleki, L.D. A (3+1)-dimensional generalized BKP-Boussinesq equation: Lie group approach. Results Phys. 2019, 13, 2211–3797. [Google Scholar] [CrossRef]
  56. Khalique, C.M.; Adeyemo, O.D. A study of (3+1)-dimensional generalized Korteweg-de Vries-Zakharov-Kuznetsov equation via Lie symmetry approach. Results Phys. 2020, 18, 103197. [Google Scholar] [CrossRef]
Figure 1. Evolution of periodic solution (20) at t = 10 and z = 1 .
Figure 1. Evolution of periodic solution (20) at t = 10 and z = 1 .
Mathematics 08 01692 g001
Figure 2. Evolution of periodic solution (20) at t = 11 and z = 0 .
Figure 2. Evolution of periodic solution (20) at t = 11 and z = 0 .
Mathematics 08 01692 g002
Figure 3. Evolution of periodic solution (20) at t = 25 and z = 9 .
Figure 3. Evolution of periodic solution (20) at t = 25 and z = 9 .
Mathematics 08 01692 g003
Figure 4. Evolution of traveling wave of singular soliton solution (24) at t = 1 and z = 0 .
Figure 4. Evolution of traveling wave of singular soliton solution (24) at t = 1 and z = 0 .
Mathematics 08 01692 g004
Figure 5. Evolution of traveling wave of singular soliton solution (24) at t = 2 and z = 3.5 .
Figure 5. Evolution of traveling wave of singular soliton solution (24) at t = 2 and z = 3.5 .
Mathematics 08 01692 g005
Figure 6. Evolution of traveling wave of singular soliton solution (24) at t = 0.2 and z = 20 .
Figure 6. Evolution of traveling wave of singular soliton solution (24) at t = 0.2 and z = 20 .
Mathematics 08 01692 g006
Figure 7. Evolution of traveling wave of singular soliton solution (25) at t = 2.1 and z = 3 .
Figure 7. Evolution of traveling wave of singular soliton solution (25) at t = 2.1 and z = 3 .
Mathematics 08 01692 g007
Figure 8. Evolution of traveling wave of singular soliton solution (25) at t = 6 and z = 25 .
Figure 8. Evolution of traveling wave of singular soliton solution (25) at t = 6 and z = 25 .
Mathematics 08 01692 g008
Figure 9. Evolution of traveling wave of singular soliton solution (25) at t = 4 and z = 30 .
Figure 9. Evolution of traveling wave of singular soliton solution (25) at t = 4 and z = 30 .
Mathematics 08 01692 g009
Figure 10. Profile of traveling wave solution (27) with parameters t = 1 and x = 2 .
Figure 10. Profile of traveling wave solution (27) with parameters t = 1 and x = 2 .
Mathematics 08 01692 g010
Figure 11. Profile of traveling wave solution (27) with parameters t = 0 and x = 3 .
Figure 11. Profile of traveling wave solution (27) with parameters t = 0 and x = 3 .
Mathematics 08 01692 g011
Figure 12. Profile of traveling wave solution (27) with parameters t = 1 and x = 4 .
Figure 12. Profile of traveling wave solution (27) with parameters t = 1 and x = 4 .
Mathematics 08 01692 g012
Figure 13. Profile of traveling wave solution (28) when t = 1 and x = 4 .
Figure 13. Profile of traveling wave solution (28) when t = 1 and x = 4 .
Mathematics 08 01692 g013
Figure 14. Profile of traveling wave solution (28) when t = 0 and x = 4 .
Figure 14. Profile of traveling wave solution (28) when t = 0 and x = 4 .
Mathematics 08 01692 g014
Figure 15. Profile of traveling wave solution (28) when t = 6 and x = 7 .
Figure 15. Profile of traveling wave solution (28) when t = 6 and x = 7 .
Mathematics 08 01692 g015
Figure 16. Evolution of traveling wave of series solution (34) at t = 1.6 and z = 1.1 .
Figure 16. Evolution of traveling wave of series solution (34) at t = 1.6 and z = 1.1 .
Mathematics 08 01692 g016
Figure 17. Evolution of traveling wave of series solution (34) at t = 1 and z = 8 .
Figure 17. Evolution of traveling wave of series solution (34) at t = 1 and z = 8 .
Mathematics 08 01692 g017
Figure 18. Evolution of traveling wave of series solution (34) at t = 5 and z = 4 .
Figure 18. Evolution of traveling wave of series solution (34) at t = 5 and z = 4 .
Mathematics 08 01692 g018

Share and Cite

MDPI and ACS Style

Masood Khalique, C.; Davies Adeyemo, O. Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science. Mathematics 2020, 8, 1692. https://doi.org/10.3390/math8101692

AMA Style

Masood Khalique C, Davies Adeyemo O. Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science. Mathematics. 2020; 8(10):1692. https://doi.org/10.3390/math8101692

Chicago/Turabian Style

Masood Khalique, Chaudry, and Oke Davies Adeyemo. 2020. "Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science" Mathematics 8, no. 10: 1692. https://doi.org/10.3390/math8101692

APA Style

Masood Khalique, C., & Davies Adeyemo, O. (2020). Closed-Form Solutions and Conserved Vectors of a Generalized (3+1)-Dimensional Breaking Soliton Equation of Engineering and Nonlinear Science. Mathematics, 8(10), 1692. https://doi.org/10.3390/math8101692

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