Next Article in Journal
Selecting the Fintech Strategy for Supply Chain Finance: A Hybrid Decision Approach for Banks
Next Article in Special Issue
Highly Dispersive Optical Solitons in Fiber Bragg Gratings with Kerr Law of Nonlinear Refractive Index
Previous Article in Journal
A Model of a Universal Neural Computer with Hysteresis Dynamics for Avionics Problems
Previous Article in Special Issue
Diverse Forms of Breathers and Rogue Wave Solutions for the Complex Cubic Quintic Ginzburg Landau Equation with Intrapulse Raman Scattering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bifurcation Theory, Lie Group-Invariant Solutions of Subalgebras and Conservation Laws of a Generalized (2+1)-Dimensional BK Equation Type II in Plasma Physics and Fluid Mechanics

by
Oke Davies Adeyemo
1,
Lijun Zhang
1,2 and
Chaudry Masood Khalique
1,*
1
Department of Mathematical Sciences, International Institute for Symmetry Analysis and Mathematical Modelling, Mafikeng Campus, North-West University, Private Bag X 2046, Mmabatho 2735, South Africa
2
College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(14), 2391; https://doi.org/10.3390/math10142391
Submission received: 21 May 2022 / Revised: 29 June 2022 / Accepted: 29 June 2022 / Published: 7 July 2022

Abstract

:
The nonlinear phenomena in numbers are modelled in a wide range of fields such as chemical physics, ocean physics, optical fibres, plasma physics, fluid dynamics, solid-state physics, biological physics and marine engineering. This research article systematically investigates a (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation. We achieve a five-dimensional Lie algebra of the equation through Lie group analysis. This, in turn, affords us the opportunity to compute an optimal system of fourteen-dimensional Lie subalgebras related to the underlying equation. As a consequence, the various subalgebras are engaged in performing symmetry reductions of the equation leading to many solvable nonlinear ordinary differential equations. Thus, we secure different types of solitary wave solutions including periodic (Weierstrass and elliptic integral), topological kink and anti-kink, complex, trigonometry and hyperbolic functions. Moreover, we utilize the bifurcation theory of dynamical systems to obtain diverse nontrivial travelling wave solutions consisting of both bounded as well as unbounded solution-types to the equation under consideration. Consequently, we generate solutions that are algebraic, periodic, constant and trigonometric in nature. The various results gained in the study are further analyzed through numerical simulation. Finally, we achieve conservation laws of the equation under study by engaging the standard multiplier method with the inclusion of the homotopy integral formula related to the obtained multipliers. In addition, more conserved currents of the equation are secured through Noether’s theorem.

1. Introduction

Fluid mechanics is a branch of physics concerning the mechanics of fluids such as liquids, gases, and plasmas and the forces on them. Applications of fluid mechanics are found in a wide range of disciplines which include civil, chemical, mechanical as well as biomedical engineering, geophysics, oceanography, astrophysics, biology and meteorology [1,2,3,4,5]. Nonlinear partial differential equations (NLPDE) in the fields of mathematics and physics play numerous important roles in theoretical sciences. They are the most fundamental models essential for studying nonlinear phenomena. Such phenomena occur in oceanography, the aerospace industry, meteorology, nonlinear mechanics, biology, population ecology, plasma physics and fluid mechanics, to mention a few. In [1] the authors studied a generalized advection–diffusion equation which is a nonlinear partial differential equation in fluid mechanics, characterizing the motion of a buoyancy propelled plume in a bent-on absorptive medium. Moreover, in [2], a generalized Korteweg–de Vries–Zakharov–Kuznetsov equation was studied. This equation delineates mixtures of warm adiabatic fluid, hot isothermal as well as cold immobile background species applicable in fluid dynamics. Furthermore, the authors of [3] considered an NLPDE where they explored the important inclined magneto-hydrodynamic flow of an upper-convected Maxwell liquid through a leaky stretched plate. In addition, the heat transfer phenomenon was studied with the heat generation and absorption effect. Plasmas considered as ‘the most abundant form of ordinary matter in the universe’ have been observed to be associated with stars which extend to the rarefied intracluster medium and possibly the intergalactic regions [4]. For instance, the authors of [4], for various types of the cosmic dusty plasmas, considered an observationally/experimentally-supported (3+1)-dimensional generalized variable-coefficient Kadomtsev–Petviashvili (KP)-Burgers-type equation. This equation could depict the dust–magneto–acoustic, dust–acoustic, magneto–acoustic, positron–acoustic, ion–acoustic, ion, electron–acoustic, quantum–dust–ion–acoustic or dust–ion–acoustic waves in one of the cosmic/laboratory dusty plasmas. The reader can access more examples in [5,6,7,8,9,10,11,12].
Observation has shown that nonlinear partial differential equations appear to model diverse physical systems, such as found in water wave theory, condensed matters, nonlinear mechanics, the aerospace industry, plasma physics, nonlinear optics lattice dynamics and so on [13,14,15,16,17,18,19]. In order to really understand these physical phenomena, it is of immense importance to secure results for differential equations (DEs) that control these aforementioned phenomena. Moreover, the research on nonlinear travelling waves (periodic, solitary, kink together with anti-kink), as well as the integrability of diverse significant nonlinear partial differential equations in the likes of the KdV equation [20], sine-Gordon equation [21] and nonlinear Schrödinger equation [22] possess vast practical values. All these involved exact solutions afford us the opportunity of being given information that aids sound understanding of the mechanism involved in the complicated physical phenomena, as well as dynamical procedures that are modelled via these nonlinear evolution equations [23].
However, no general and systematic theory was available to be applied to NLPDEs so that their closed-form solutions can be obtained. Nonetheless, in recent times mathematicians and physicists have evolved effective techniques to achieve viable analytical solutions to NLPDEs, such as inverse scattering transform [13], Bäcklund transformation [24], F-expansion technique [25], extended simplest equation approach [26], Lie symmetry analysis [27,28,29,30,31], the G G —expansion technique [32], Darboux transformation [33], sine-Gordon equation expansion technique [34] as well as the Kudryashov approach [35], modified extended direct algebraic approach [36,37], the sine-cosine method [11], Hirota’s bilinear technique [38], the exp-function expansion technique [12], and the auxiliary ordinary differential equation approach [10]; the list continues.
Furthermore, in recent years, the bifurcation technique [39] among other techniques has been used for obtaining both bounded and unbounded solutions of NLPDE. This technique allows for the extensive study of the dynamical performance of the analytic travelling wave solutions as well as their phase portrait analysis via the engagement of the theory of dynamical systems. In [40] Jiang et al. investigated the dynamical behaviour of points of equilibrium together with the bifurcations of phase portraits involved in the travelling wave results for the CH- γ equation. In addition, Saha [41] also exhibited the existence of smooth alongside non-smooth travelling wave solutions of generalized KP-MEW equations by the exploitation of the bifurcation theory of planar dynamical system. Das et al. [42,43] equally examined the existence together with stability analysis of the dispersive solution of the KP-BBM as well as KP equations with the prevalence of dispersion consequence.
A two-dimensional generalization of the well-recognized Korteweg–de Vries equation yields the Bogoyavlensky–Konopelchenko equation [44]:
p t + α p x x x + β p x x y + 6 α p p x + 4 β p p y + 4 β p x x 1 p y = 0 ,
with constant coefficients α and β , where x 1 = d x . Inserting x 1 p = u into Equation (1), one attains the equivalent structure of (1) as [45]:
u t x + 6 α u x u x x + 4 β u x u x y + 4 β u y u x x + α u x x x x + β u x x x y = 0 .
In [45] with u y = v x in (2), the authors integrated the result once to obtain a system of NLPDE. Further, they utilized the Lie group theoretic approach to obtain solutions to the system of equations. Added to that is the fact that they engaged the method to secure conservation laws of the equations. Besides, the authors employed a new concept of nonlinear self-adjointness of differential equations in conjunction with formal Lagrangian for constructing nonlocal conservation laws of the system. In [46], Triki et al. investigated the Bogoyavlensky–Konopelchenko Equation (2) and secured some shock wave solutions to the equation. In addition, various applications of (2) were highlighted in [45,46]. This established version describes an interconnection of a long wave propagation directed towards the x-axis together with a Riemann wave propagation directed towards the y-axis [47]. Some authors examined (2) with 4 β replaced by 3 β and secured the solution of the resultant model. For instance, a Darboux transformation as well as some travelling wave solutions were given in [48] for Equation (2). We note that the replacement earlier mentioned presents Equation (2) as a special case of the KdV model in [49]. In addition to that, a few particular properties of the equation have also been explored.
Chen et al. [50] contemplated the NLPDE called (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation stated as:
v t x + α ( 6 v x v x x + v x x x x ) + β ( v x x x y + 3 v x v x y + 3 v x x v y ) + γ 1 v x x + γ 2 v x y + γ 3 v y y = 0 ,
which exists in plasma physics and fluid mechanics with α , β , γ 1 , γ 2 , γ 3 , nonzero real valued constants and v = v ( t , x , y ) . The authors got the Lump-type solutions together with lump solutions of (3) with the employment of symbolic computation given in Hirota bilinear form [51] as:
( D t D x + α D x 4 + β D x 3 D y + σ D x D y + γ D x 2 + ν D y 2 ) f · f = 0 ,
achieved under the transformations:
u = 12 α β 1 ln f x x , v = 12 α β 1 ln f x ,
with nonzero real constants σ , γ and ν , where f is an analytic function depending on x, y and t, D x , D y and D t are regarded as the bilinear derivative operators given by [38,51], which they used in constructing new closed-form and explicit solutions that include two-wave alongside polynomial solutions for the equation. In addition, the lump-type solution found comprises eleven parameters together with six independent parameters (arbitrary), as well as non-zero conditions. Not only that, lump solutions were achieved by considering a particular class of parameters, the motion track of which is also theoretically and graphically delineated. In the same vein, lower-order lump solution of (3) has been presented [52]. The authors of [53] confirmed in their work the existence of diverse wave structures for (3) delineating nonlinear waves in applied sciences. In this regard, on the basis of Hirota’s bilinear structure and diverse test schemes, various kinds of exact solutions, comprising breather-wave, double soliton, rational, cross-kink, mixed-type, as well as interaction solutions to the equation, were formally extracted.
Moreover, in [54], the authors considered a version of (3) in the form:
u t x + k 1 u x x x x + k 2 u x x x y + 2 k 1 k 3 k 2 u x u x x + k 3 u x u x y + k 3 u x x u y + γ 1 u x x + γ 2 u x y + γ 3 u y y = 0 ,
with real function u = u ( x , y , t ) with scaled time variable t as well as scaled space variables x , y and real constants k 1 , k 2 , k 3 , γ 1 , γ 2 , γ 3 . They went ahead to examine the equation which applies in fluid mechanics and plasma physics by utilizing the Lie symmetry technique to obtain symmetries of the equation. Besides, the G G —expansion technique, polynomial expansion as well as power series expansion methods were adopted to achieve some solutions of the equation by the authors.
In this article, we investigate the (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation ((2+1)-D genBKe), a version of (3) structured as:
Δ u t x + α ( 6 u x u x x + u x x x x ) + β u x x x y + 3 ( ρ u x u x y + δ u x x u y ) + γ u x x + σ u x y + ν u y y = 0 ,
applicable in plasma physics and fluid mechanics with α , β , σ , γ , ν , ρ and δ as nonzero real valued constants. In the study, we carry out explicit solutions of the (2+1)-D genBKe (4) to achieve its abundant closed-form and travelling wave solutions. Thus, we catalogue the article in the subsequent format. Section 2, presents the Lie group analysis of Equation (4) where the obtained generators are adopted in computing its optimal system of Lie subalgebras. In addition, each Lie subalgebra is explored to reduce (4) and obtain solutions of the underlying equation. In Section 3, we adopt the bifurcation theory of the dynamical system to secure some nontrivial travelling wave solutions of the under-study equation. Numerical simulations of the secured solutions are conducted for further analysis and discussion in Section 4. Furthermore, Section 5 furnishes the conservation laws of (2+1)-D genBKe to be constructed via the standard multiplier technique with the use of the homotopy formula. In addition, we engage Noether’s theorem to gain more conserved vectors of (4) with ρ = 2 δ . Shortly after, we present the concluding remarks.

2. Lie Symmetry Analysis

This section first presents the algorithm for the computation of the Lie point symmetries of (2+1)-D genBKe (4) together with its differential generators. Thereafter, we engage them to calculate the optimal system of Lie subalgebras and utilize them to generate exact solutions for (4).

2.1. Lie Point Symmetries

Here in this subsection, we contemplate the one-parameter Lie group of infinitesimal transformations
t ˜ t + ε ξ 1 ( t , x , y , u ) + O ( ε 2 ) , x ˜ x + ε ξ 2 ( t , x , y , u ) + O ( ε 2 ) , y ˜ y + ε ξ 3 ( t , x , y , u ) + O ( ε 2 ) , u ˜ u + ε η ( t , x , y , u ) + O ( ε 2 ) ,
with ε standing for the parameter of the group alongside ξ 1 , ξ 2 , ξ 3 , η serving as the infinitesimals of the transformations depending on t, x, y, and u. Thus utilizing ε (one-parameter), Lie group of infinitesimal transformation in compliance with invariant conditions [55,56], solution space ( t , x , y , u ) of (2+1)-D genBKe (4) stays invariant and can also transform into another space ( t ˜ , x ˜ , y ˜ , u ˜ ) .
In accordance with the technique for deciding the infinitesimal generators of nonlinear differential equations (NLDE), we shall secure the infinitesimal generator of (4). Symmetry group of (2+1)-D genBKe (4) will be found by exploring vector field:
X = ξ 1 ( t , x , y , u ) t + ξ 2 ( t , x , y , u ) x + ξ 3 ( t , x , y , u ) y + η ( t , x , y , u ) u ,
where ξ i , i = 1 , 2 , 3 such that ξ s and η are functions depending on t, x, y alongside u. We recall that (5) is a symmetry of (2+1)-D genBKe (4) if invariance condition,
p r ( 4 ) X Δ | Δ = 0 = 0 ,
holds. Here p r ( 4 ) X denotes the fourth prolongation of ( X ) [29] defined by:
p r ( 4 ) X = X + ζ t u t + ζ x u x + ζ y u y + ζ t x u t x + ζ x x u x x + ζ x y u x y + ζ y y u y y + ζ x x x x u x x x x + ζ x x x y u x x x y ,
with the ζ t , ζ x , ζ y , ζ t x , ζ x x , ζ x y , ζ y y , ζ x x x x and ζ x x x y , given as:
ζ t = D t ( η ) u t D t ( ξ 1 ) u x D t ( ξ 2 ) u y D t ( ξ 3 ) , ζ x = D x ( η ) u t D x ( ξ 1 ) u x D x ( ξ 2 ) u y D x ( ξ 3 ) , ζ y = D y ( η ) u t D y ( ξ 1 ) u x D y ( ξ 2 ) u y D y ( ξ 3 ) , ζ t x = D x ( ζ t ) u t t D x ( ξ 1 ) u t x D x ( ξ 2 ) u t y D x ( ξ 3 ) , ζ x x = D x ( ζ x ) u t x D x ( ξ 1 ) u x x D x ( ξ 2 ) u x y D x ( ξ 3 ) , ζ x y = D x ( ζ y ) u t y D x ( ξ 1 ) u y x D x ( ξ 2 ) u y y D x ( ξ 3 ) , ζ y y = D y ( ζ y ) u t y D y ( ξ 1 ) u x y D y ( ξ 2 ) u y y D y ( ξ 3 ) , ζ x x x x = D x ( ζ x x x ) u x x x t D x ( ξ 1 ) u x x x x D x ( ξ 2 ) u x x x y D x ( ξ 3 ) , ζ x x x y = D x ( ζ x x y ) u x x y t D x ( ξ 1 ) u x x x y D x ( ξ 2 ) u x x y y D x ( ξ 3 ) ,
and the total derivatives D t , D x as well as D y defined as:
D t = t + u t u + u t t u t + u t x u x + , D x = x + u x u + u x x u x + u t x u t + , D y = y + u y u + u y y u y + u y t u t + .
Writing out the expanded form of determining Equation (6) and splitting it over the various derivatives of u, we get twenty-two overdetermined systems of linear partial differential equations:
ξ u 2 = 0 , ξ u 3 = 0 , ξ u 1 = 0 , η u u = 0 , ξ y 1 = 0 , ξ x 1 = 0 , ξ x 3 = 0 , η x u = 0 , η y u ξ x y 2 = 0 , η y u 3 ξ x y 2 = 0 , η u + ξ x 2 = 0 , ξ y 3 3 ξ x 2 = 0 , η u + ξ x 2 = 0 , α ξ x 2 + β ξ y 2 α ξ y 3 = 0 , 2 α η u δ ξ y 2 ρ ξ y 2 + 2 α ξ y 2 = 0 , ξ x x 2 = 0 , η t u + 6 α η x x + 3 ρ η x y + σ η y u ξ t x 2 σ ξ x y 2 ν ξ y y 2 + 4 α η x x x u + 3 β η x x y u = 0 , 6 α η x + 3 δ η y ξ t 2 + γ ξ x 2 σ ξ y 2 + γ ξ y 3 + 6 α η x x u + 3 β η x y u = 0 , 2 ξ x 2 ξ t 1 + ξ y 3 = 0 , 3 δ η x x + 2 ν η y u ν ξ y y 3 + β η x x x u = 0 , η t x + γ η x x + σ η x y + ν η y y + α η x x x x + β η x x x y = 0 , 3 ρ η x + 2 σ ξ x 2 2 ν ξ y 2 ξ t 3 + 3 β η x x u = 0 .
Solving the system of linear PDEs via symbolic software MathLie, one procures ξ 1 , ξ 2 , ξ 3 and η given as:
ξ 1 = c 1 , ξ 2 = f 1 ( t ) , ξ 3 = c 2 + c 3 t , η = 1 3 δ ρ δ c 3 x 2 α c 3 y + 3 δ ρ f 2 ( t ) + ρ y f 1 ( t ) .
If we define arbitrary functions f 1 ( t ) and f 2 ( t ) as f 1 ( t ) = c 4 and f 2 ( t ) = c 5 , where c 4 and c 5 are arbitrary constants, thus with the aid of (5), the solution purveys vectors:
X 1 = x , X 2 = y , X 3 = t , X 4 = u , X 5 = t y + x 3 ρ 2 α 3 δ ρ y u .
Theorem 1.
(2+1)-D genBK Equation (4) admits a five dimensional Lie algebra L 5 spanned by the vectors X 1 , , X 5 .
The associated group transformations for X 1 , , X 5 are
G 1 : ( t ˜ , x ˜ , y ˜ , u ˜ ) ( t , x + ε 1 , y , u ) , G 2 : ( t ˜ , x ˜ , y ˜ , u ˜ ) ( t , x , y + ε 2 , u ) , G 3 : ( t ˜ , x ˜ , y ˜ , u ˜ ) ( t + ε 3 , x , y , u ) , G 4 : ( t ˜ , x ˜ , y ˜ , u ˜ ) ( t , x , y , u + ε 4 ) , G 5 : ( t ˜ , x ˜ , y ˜ , u ˜ ) t , x , y + ε 5 t , u + ε 5 3 ρ 2 α ε 5 3 δ ρ y α ε 5 2 3 δ ρ t ,
with ε 1 , , ε 5 representing real numbers. We realize that G 1 portrays the x-translation, G 2 the y-translation and G 3 the t-translation.
Theorem 2.
If u = f ( t , x , y ) is a solution of the (2+1)-D genBKe (4), then so are the functions presented as:
G 1 ( ε 1 ) : u ( t , x , y ) = f ( t , x ε 1 , y ) , G 2 ( ε 2 ) : u ( t , x , y ) = f ( t , x , y ε 2 ) , G 3 ( ε 3 ) : u ( t , x , y ) = f ( t ε 3 , x , y ) , G 4 ( ε 4 ) : u ( t , x , y ) = f ( t , x , y ) + ε 4 , G 5 ( ε 5 ) : u ( t , x , y ) = f ( t , x , y ε 5 t ) ε 5 3 ρ + 2 α ε 5 3 δ ρ y + α ε 5 2 3 δ ρ t .

2.2. Optimal System of One-Dimensional Subalgebras

It is revealed that it is unfeasible to list all possible group-invariant solutions. As a result, the situation necessitates an effective, systematic and efficient means of classifying these solutions. The moment this is achieved, the optimal system of group-invariant solutions is then formed. Ibragimov et al. [57] invoke a robust approach that depends on the commutator table in achieving the one-dimensional subalgebras optimal system. In consequence, we give the commutator table (table of Lie brackets) of (4) associated with (8) in Table 1, that is
We state here that apparently { X 1 , X 2 , X 3 , X 4 , X 5 } is closed under the Lie bracket. Besides, we express an arbitrary operator X L 5 as:
X = l 1 X 1 + l 2 X 2 + l 3 X 3 + l 4 X 4 + l 5 X 5 .
In a bid to secure the linear transformations related to vector l = ( l 1 , l 2 , l 3 , l 4 , l 5 ) , we have the generator defined as:
E i = c i j k l j l k , i = 1 , 2 , 3 , 4 , 5 ,
with c i j k given for the relation [ X i , X j ] = c i j k X k . On taking cognizance of Equation (10) alongside Table 1, generators E 1 , E 2 , E 3 , E 4 , E 5 are presented as:
E 1 = δ l 5 l 4 , E 2 = 2 α l 5 l 4 , E 3 = 3 δ ρ l 5 l 2 , E 4 = 0 , E 5 = 2 α l 2 l 4 δ l 1 l 4 3 δ ρ l 3 l 2 .
In association with E 1 , E 2 , E 3 , E 4 and E 5 , we give the Lie equations possessing parameters a 1 , a 2 , a 3 , a 4 and a 5 having the initial criteria l ˜ | a i = 0 = l , i = 1 , 5 , as
d l ˜ 1 d a 1 = 0 , d l ˜ 2 d a 1 = 0 , d l ˜ 3 d a 1 = 0 , d l ˜ 4 d a 1 = δ l ˜ 5 , d l ˜ 5 d a 1 = 0 , d l ˜ 1 d a 2 = 0 , d l ˜ 2 d a 2 = 0 , d l ˜ 3 d a 2 = 0 , d l ˜ 4 d a 2 = 2 α l ˜ 5 , d l ˜ 5 d a 2 = 0 , d l ˜ 1 d a 3 = 0 , d l ˜ 2 d a 3 = 3 δ ρ l ˜ 5 , d l ˜ 3 d a 3 = 0 , d l ˜ 4 d a 3 = 0 , d l ˜ 5 d a 3 = 0 , d l ˜ 1 d a 4 = 0 , d l ˜ 2 d a 4 = 0 , d l ˜ 3 d a 4 = 0 , d l ˜ 4 d a 4 = 0 , d l ˜ 5 d a 4 = 0 , d l ˜ 1 d a 5 = 0 , d l ˜ 2 d a 5 = 3 δ ρ l ˜ 3 , d l ˜ 3 d a 5 = 0 , d l ˜ 4 d a 5 = δ l ˜ 1 + 2 α l ˜ 2 , d l ˜ 5 d a 5 = 0 .
Consequently, we give the transformations involved in the solution of Equations (11) as
T 1 : l ˜ 1 = l 1 , l ˜ 2 = l 2 , l ˜ 3 = l 3 , l ˜ 4 = l 4 + δ a 1 l 5 , l ˜ 5 = l 5 , T 2 : l ˜ 1 = l 1 , l ˜ 2 = l 2 , l ˜ 3 = l 3 , l ˜ 4 = l 4 2 α a 2 l 5 , l ˜ 5 = l 5 , T 3 : l ˜ 1 = l 1 , l ˜ 2 = l 2 + 3 δ ρ a 3 l 5 , l ˜ 3 = l 3 , l ˜ 4 = l 4 , l ˜ 5 = l 5 , T 4 : l ˜ 1 = l 1 , l ˜ 2 = l 2 , l ˜ 3 = l 3 , l ˜ 4 = l 4 , l ˜ 5 = l 5 , T 5 : l ˜ 1 = l 1 , l ˜ 2 = l 2 3 δ ρ a 5 l 3 , l ˜ 3 = l 3 , l ˜ 4 = l 4 3 α δ ρ a 5 2 l 3 + 2 α a 5 l 2 δ a 5 l 1 , l ˜ 5 = l 5 .

Optimal Classification

We observe the fact that the transformations T i , i = 1 , , 5 actually map vector X L 5 presented by (9) to vector X ˜ L 5 expressed via the relation:
X ˜ = l ˜ 1 X 1 + l ˜ 2 X 2 + l ˜ 3 X 3 + l ˜ 4 X 4 + l ˜ 5 X 5 .
The technique involved in the construction of optimal system in this process demands the simplification of general vector structured as:
l = ( l 1 , l 2 , l 3 , l 4 , l 5 ) ,
by engaging transformations T 1 , T 2 , T 3 , T 4 , T 5 . We are captivated to seek for simplest representative of each class of alike vectors of (12) by inserting these representatives in (9) and so, we gain one-dimensional subalgebras optimal system of (2+1)-D genBKe (4). Thus, we structured the classifications into two different cases.
  • Case 1.  l 5 0
  • 1.1.  l 1 = 0 ,
We contemplate transformation T 3 by taking a 3 = l 2 3 δ ρ l 5 , we can then make l ˜ 2 = 0 . Thus vector (12) reduces to the structure:
l = ( 0 , 0 , l 3 , l 4 , l 5 ) .
Moreover, if we take a 1 = l 4 δ l 5 from T 1 which makes l ˜ 4 = 0 , then we further reduce vector (13) to:
l = ( 0 , 0 , l 3 , 0 , l 5 ) .
Evidently, since (14) cannot be further reduced, without loss of generality, we assume that l 3 = 1 and l 5 = ± 1 . Therefore, we have the optimal representative:
X 3 ± X 5 .
Next, we contemplate the case of l 3 0 and first consider the resultant subalgebra when l 2 0 .
  • 1.1.1.  l 3 0 ,
  • 1.1.1.1.  l 2 0 ,
By taking a 2 = l 4 2 α l 5 from transformation T 1 , we can make l ˜ 4 = 0 . Now, since l 1 = 0 and l 2 = l 3 = l 5 0 , then vector (12) becomes:
l = ( 0 , l 2 , l 3 , 0 , l 5 ) .
If we suppose that l 2 = 1 and l 3 = l 5 = ± 1 , then we have the representative
X 2 ± X 3 ± X 5 .
Remark 1.
We notice here that for the case of l 2 = 0 , we achieve an optimal representative earlier obtained and consequently contribute no additional subalgebra to the optimal system.
  • 1.1.2.  l 3 = 0 .
We take, in this case, a 3 = l 2 3 δ ρ l 5 from T 3 , so that we make l ˜ 2 = 0 . In addition, by considering a 5 = l 4 2 α l 2 δ l 2 in T 5 , thereby making l ˜ 4 = 0 , we secure vector:
l = ( 0 , 0 , 0 , 0 , l 5 )
and so we have the optimal representative:
X 5 .
  • 1.1.2.1.  l 4 0 .
By taking a 5 = l 2 3 δ ρ l 3 from T 5 , we have the reduced form of vector (12) as
l = ( 0 , 0 , 0 , l 4 , l 5 ) ,
which can not be simplified further and so we gain the representative:
X 4 ± X 5 .
Now, we contemplate some subcases when l 1 0 with a view to obtaining all possible optimal representatives.
  • 1.2.  l 1 0 ,
  • 1.2.1.  l 4 = 0 ,
  • 1.2.1.1.  l 3 0 ,
By making a 3 = l 2 3 δ ρ l 5 in transformation T 3 which occasions the possibility of making l ˜ 2 = 0 , we have the vector:
l = ( l 1 , 0 , l 3 , 0 , l 5 ) ,
which we can not further streamline and so we gain the optimal representative:
X 1 ± X 3 ± X 5 .
  • 1.2.1.2.  l 3 = 0 .
By taking in transformation T 5 , a 5 = l 4 2 α l 2 δ l 2 and a 5 = l 2 3 δ ρ l 3 , we have the vector:
l = ( 0 , 0 , 0 , l 4 , l 5 ) ,
which can not be simplified further and so we gain the representative:
X 1 ± X 5 .
Next, we consider the case of l 4 0 and then take into account the resultant subalgebra when l 3 = 0 .
  • 1.2.2.  l 4 0 ,
  • 1.2.2.1.  l 3 = 0 ,
By taking a 3 = l 2 3 δ ρ l 5 in transformation T 3 , we make l ˜ 2 = 0 and so we have vector:
l = ( l 1 , 0 , 0 , l 4 , l 5 ) ,
which gives rise to the optimal representative:
X 1 ± X 4 ± X 5 .
We reveal here that remark (1) absolutely applies to the case of l 4 = 0 and l 3 0 .
  • Case 2.  l 5 = 0 .
In this second part of the process, we contemplate the structure of vector (12) as:
l = ( l 1 , l 2 , l 3 , l 4 , 0 ) .
Finally, we consider the case of l 4 0 and then take into account the optimal representatives when l 1 = 0 .
  • 2.1.  l 4 0 ,
  • 2.1.1.  l 1 = 0 .
By contemplating the parameter a 5 = l 2 3 δ ρ l 3 in transformation T 5 , one can definitely make l ˜ 2 = 0 and so, we have the reduced form of vector (22) to be given as:
l = ( 0 , 0 , l 3 , l 4 , 0 ) ,
which consequently yields the optimal representative:
X 3 ± X 4 .
  • 2.1.2.  l 1 0 .
Conversely, if we consider l 1 0 with l 3 = 0 , using T 3 where a 3 = l 2 3 δ ρ l 5 , occasions vector (22) giving us:
l = ( l 1 , 0 , 0 , l 4 , 0 )
and so we gain the subalgebra
X 1 ± X 4 .
  • 2.2.  l 4 = 0 .
By taking l 2 0 and also considering the converse ( l 2 = 0 ) with the use of T 5 where a 5 = l 2 3 δ ρ l 3 , we gain the respective subalgebras:
X 1 ± X 2 ± X 3 , X 1 ± X 3 .
  • 2.2.1.  l 3 = 0 ,
If we take the parameter a 5 = l 2 3 δ ρ l 3 in transformation T 5 , that is l ˜ 2 = 0 , one gets:
X 1 .
Finally, if we take l 1 = 0 with l 2 0 and in addition contemplate a case of l 3 0 with l 1 = 0 , we get in the respective situations:
X 2 , X 3 .
Conclusively, by gathering the operators secured (that is, (15)–(21), (23)–(25) and (27)), we arrive at a theorem, which is:
Theorem 3.
The subsequent operators provide an optimal system of one-dimensional subalgebras of the Lie algebra which is spanned by vectors X 1 , X 2 , X 3 , X 4 , X 5 of (2+1)-D genBKe (4):
X 1 , X 2 , X 3 , X 5 , X 3 ± X 5 , X 4 ± X 5 , X 1 ± X 5 , X 3 ± X 4 , X 1 ± X 3 , X 1 ± X 4 , X 2 ± X 3 ± X 5 , X 1 ± X 3 ± X 5 , X 1 ± X 4 ± X 5 , X 1 ± X 2 ± X 3 .

2.3. Group-Invariants and Some Exact Solutions

This subsection presents group-invariant solutions of (2+1)-D genBKe (4) by exploring results presented in Theorem 3. Thus, furnishing some exact solutions of (4). Therefore, we utilize the Lagrangian system given as [27,29]:
d t ξ 1 ( t , x , y , u ) = d x ξ 2 ( t , x , y , u ) = d y ξ 3 ( t , x , y , u ) = d u η ( t , x , y , u ) ,
to secure the group-invariant solutions related to the vector fields.

2.3.1. Optimal Subalgebra X 1

The characteristic equation corresponding to optimal subalgebra X 1 = / x is
d t 0 = d x 1 = d y 0 = d u 0 .
On solving system (28), one gains invariants alongside their group-invariant as:
T = t , Y = y , where u ( t , x , y ) = G ( T , Y ) .
Therefore, by using the functions and variables from (29) in (4), we obtain:
G Y Y = 0 ,
which gives a solution in terms of T and Y but by back-substitution, we have
u ( t , x , y ) = f 1 ( t ) y + f 2 ( t ) .
Arbitrary functions f 1 and f 2 are depending on t in (30), a solution of (4).

2.3.2. Optimal Subalgebra X 2

The group-invariant associated with optimal subalgebra X 2 = / y is calculated as:
u ( t , x , y ) = G ( T , X ) , with T = t , X = x .
On utilizing the obtained group-invariant, (2+1)-D genBKe (4) is transformed to:
G T X + 6 α G X G X X + α G X X X X + γ G X X = 0 .
As a consequence, we gain a logarithmic-hyperbolic function solution in this regard as:
G ( T , X ) = 2 A 2 tanh A 1 T + A 2 X + A 0 + A 2 ln tanh A 1 T + A 2 X + A 0 1 tanh A 1 T + A 2 X + A 0 + 1 + 4 3 A 2 2 X γ 6 α X A 1 6 α A 2 X + f ( T ) d T ,
where A 0 , A 1 as well as A 2 are arbitrary constants. Therefore, on retrograding to the basic variables, one achieves a solution of (2+1)-D genBKe (4) in this case as:
u ( t , x , y ) = 2 A 2 tanh A 1 t + A 2 x + A 0 + A 2 ln tanh A 1 t + A 2 x + A 0 1 tanh A 1 t + A 2 x + A 0 + 1 + 4 3 A 2 2 x γ 6 α x A 1 6 α A 2 x + f ( t ) d t .
Further investigation of PDE (32) reveals that it has four Lie point symmetries,
R 1 = T + F 1 ( T ) G , R 1 = X + F 2 ( T ) G , R 3 = T X + 1 6 α X + F 3 ( T ) G , R 4 = T T + 1 3 X X + F 4 ( T ) γ 9 α X 1 3 G G .
We contemplate some special cases of the generators obtained. Letting F 1 ( T ) = 1 , we have solution of R 1 as G ( T , X ) = T + ϕ ( r ) , r = X , that further reduces (4) to:
γ ϕ ( r ) + 6 α ϕ ( r ) ϕ ( r ) + α ϕ ( r ) = 0 ,
whose result furnishes a trigonometric function solution of (2+1)-D genBKe (4) as:
u ( t , x , y ) = t γ α tan γ 4 α x ± α C 0 + C 1 .
C 0 and C 1 are integration constants. Moreover, taking F 2 ( T ) = 1 , we have G ( T , X ) = X + ϕ ( r ) , r = T , which gives a trivial solution. Besides, for F 1 ( T ) = F 2 ( T ) = 0 , we consider a linear combination Q = c 0 R 1 + c 1 R 2 whose solution is G ( T , X ) = ϕ ( r ) , r = c 0 X c 1 T . Utilizing the gained outcome, we reduce Equation (4) to:
γ c 0 ϕ ( r ) c 1 ϕ ( r ) + 6 α c 0 2 ϕ ( r ) ϕ ( r ) + α c 0 3 ϕ ( r ) = 0 .
On solving nonlinear ordinary differential equation (NODE) (35), we secure:
u ( t , x , y ) = C 1 c 1 γ c 0 α c 0 tanh 1 2 c 0 3 / 2 c 1 γ c 0 α c 0 3 / 2 α C 0 ( c 0 x c 1 t ) ,
which is an hyperbolic solution of (4) with C 0 and C 1 , integration constants. In addition, taking F 3 ( T ) = 0 , we have outcome G ( T , X ) = X 2 / 12 α T + ϕ ( r ) , r = T , which gives no solution of interest. Besides, for F 4 ( T ) = 0 , we have the result G ( T , X ) = T 1 / 3 ϕ ( r ) γ X / 6 α , r = X T 1 / 3 which eventually transforms (4) to:
18 α ϕ ( r ) ϕ ( r ) + 3 α ϕ ( r ) r ϕ ( r ) 2 ϕ ( r ) = 0 .

2.3.3. Optimal Subalgebra X 3

Lie optimal subalgebra X 3 = / t reduces (2+1)-D genBKe (4) to the PDE
σ G X Y + γ G X X + ν G Y Y + 6 α G X G X X + 3 ρ G X G X Y + 3 δ G Y G X X + α G X X X X + β G X X X Y = 0
through the group-invariant alongside its invariants calculated and presented as
u ( t , x , y ) = G ( X , Y ) , whereas X = x , Y = y .
Consequently, we secure a solution of (37) with respect to X and Y but by back-substitution, we find a steady-state hyperbolic solution of (4) in this regard as:
u ( t , x , y ) = Ω 0 ρ + Ω 0 δ + 4 α ν δ σ ρ σ 4 A 1 2 β δ 4 A 1 2 β ρ cosh Ω 1 2 ν 1 × { 4 Ω 0 A 1 β sinh Ω 1 2 ν 16 A 1 3 β 2 sinh Ω 1 2 ν 4 A 1 2 A 2 β δ cosh Ω 1 2 ν 4 A 1 2 A 2 β ρ cosh Ω 1 2 ν + 8 A 1 α ν sinh Ω 1 2 ν 4 A 1 β σ sinh Ω 1 2 ν + Ω 0 A 2 δ cosh Ω 1 2 ν + Ω 0 A 2 ρ cosh Ω 1 2 ν + 4 A 2 α ν cosh Ω 1 2 ν A 2 δ σ cosh Ω 1 2 ν A 2 ρ σ cosh Ω 1 2 ν } ,
where Ω 0 = 16 A 1 4 β 2 16 α ν A 1 2 + 8 σ β A 1 2 4 γ ν + σ 2 , Ω 1 = Ω 0 A 1 y 4 A 1 3 β y + 2 A 1 ν x A 1 σ y + 2 A 0 ν , where A 0 and A 1 are arbitrary constants of solution. On performing the Lie symmetry analysis on (37), we obtain translation symmetries
R 1 = X , R 2 = Y , R 3 = G .
We contemplate the linear combination of the three generators as Q = c 0 / X + c 1 / Y + c 2 / G . Therefore, Q furnishes the solution G ( X , Y ) = c 2 / c 0 X + ϕ ( r ) , where r = Y c 1 / c 0 X . Engaging the function and its variables, we reduce (4) to:
α c 1 4 ϕ ( 4 ) ( r ) β c 1 3 c 0 ϕ ( 4 ) ( r ) + 6 α c 1 2 c 2 c 0 ϕ ( r ) + γ c 1 2 c 0 2 ϕ ( r ) + c 0 4 ν ϕ ( r ) 3 c 1 c 2 c 0 2 ρ ϕ ( r ) c 1 c 0 3 σ ϕ ( r ) 6 α c 1 3 c 0 ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 δ ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 ρ ϕ ( r ) ϕ ( r ) = 0 .
On solving the fourth-order NODE (39), we achieve the trigonometric function:
u ( t , x , y ) = ± 1 c 0 Δ 0 c 1 2 ( 2 α c 1 c 0 ( δ + ρ ) ) { 2 Δ 0 i Δ 1 tan [ Δ 0 2 Δ 1 ( α c 1 4 A 1 β c 0 c 1 3 A 1 i ( c 0 y c 1 x ) Δ 1 c 0 ) ] } + c 2 c 0 x + A 2 ,
where Δ 0 = σ c 0 2 c 1 ν c 0 3 6 α c 1 2 c 2 + c 0 c 1 ( 3 ρ c 2 γ c 1 ) , Δ 1 = c 1 3 ( α c 1 β c 0 ) with constant of integrations A 1 and A 2 . We observe that the obtained result presented in (40) is a steady-state complex trigonometric function solution of (4).

2.3.4. Optimal Subalgebra X 3 + a X 5 , a { 1 , 1 }

The group-invariant related to subalgebra X 3 + a X 5 is calculated and presented as:
u ( t , x , y ) = G ( X , Y ) + 2 a 2 α 9 δ ρ t 3 + a 3 ρ x 2 a α 3 δ ρ y , where X = x , Y = y 1 2 a t 2 .
Invoking the function given in (41) along with the variables, we transform (4) to:
a G X Y + σ G X Y + γ G X X + ν G Y Y a t G X Y + 3 ρ G X G X Y + 6 α G X G X X + 3 δ G Y G X X + α G X X X X + β G X X X Y = 0 .
On applying the Lie theoretic approach on (42), we achieve three generators:
R 1 = X , R 2 = Y , R 3 = G .
Now, the similarity solution of R 1 = / X purveys G ( X , Y ) = ϕ ( r ) , with r = Y . Thus using the function reduces (4) to differential equation ϕ ( r ) = 0 whose solution is:
ϕ ( r ) = A 0 r + A 1 ,
where A 0 and A 1 are integration constants. On retrograding to the basic variables,
u ( t , x , y ) = 2 a 2 α 9 δ ρ t 3 + a 3 ρ x 2 a α 3 δ ρ y + A 0 y 1 2 a t 2 + A 1 .
Next, we gain the solution related to generator R 2 as G ( X , Y ) = ϕ ( r ) , with r = X . In consequence, we reduce Equation (4) to a fourth-order NODE expressed as:
γ ϕ ( r ) + 6 α ϕ ( r ) ϕ ( r ) + α ϕ ( r ) = 0 .
Thus, on solving the NODE and reverting to the fundamental variables, one obtains:
u ( t , x , y ) = 2 a 2 α 9 δ ρ t 3 + a 3 ρ x 2 a α 3 δ ρ y γ α tan γ 4 α x ± α A 1 + A 2 ,
with A 1 and A 2 , integration constants. On contemplating the combination of R 1 and R 2 as Q = c 0 R 1 + c 1 R 2 . In consequence, Q furnishes the solution G ( X , Y ) = ϕ ( r ) , where r = Y c 1 / c 0 X . Imploring the function and its variables transforms (4) to:
a c 1 c 0 3 t ϕ ( r ) a c 1 c 0 3 ϕ ( r ) + α c 1 4 ϕ ( 4 ) ( r ) β c 1 3 c 0 ϕ ( 4 ) ( r ) + γ c 1 2 c 0 2 ϕ ( r ) + c 0 4 ν ϕ ( r ) c 1 c 0 3 σ ϕ ( r ) 6 α c 1 3 c 0 ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 δ ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 ρ ϕ ( r ) ϕ ( r ) = 0 .
On solving NODE (45), we secure a complex tan-hyperbolic solution of (4) as:
u ( t , x , y ) = 2 a 2 α 9 δ ρ t 3 + a 3 ρ x 2 a α 3 δ ρ y ± 2 i Ω 1 c 1 2 ( 2 α c 1 c 0 ( δ + ρ ) ) tanh { Ω 2 2 c 1 3 ( α c 1 β c 0 ) × α c 1 4 A 1 β c 0 c 1 3 A 1 i c 0 ( c 1 3 ( α c 1 β c 0 ) ) y a 2 t 2 c 1 c 0 x } + A 2 ,
where Ω 1 = c 1 3 ( α c 1 β c 0 ) ( ν c 0 2 + γ c 1 2 + c 0 c 1 ( a ( t 1 ) σ ) ) , A 1 together with A 2 constant of integration and Ω 2 = c 0 ( ν c 0 2 + γ c 1 2 + c 0 c 1 ( a ( t 1 ) σ ) ) . Furthermore, we contemplate the combinations of all the symmetries as Q = c 0 R 1 + c 1 R 2 + c 2 R 3 . Hence, Q produces the solution G ( X , Y ) = c 2 / c 0 X + ϕ ( r ) , where r = Y c 1 / c 0 X . On utilizing function G ( X , Y ) as well as its variables, we reduce (4) to NODE
a c 1 c 0 3 t ϕ ( r ) a c 1 c 0 3 ϕ ( r ) + α c 1 4 ϕ ( 4 ) ( r ) β c 1 3 c 0 ϕ ( 4 ) ( r ) + 6 α c 1 2 c 2 c 0 ϕ ( r ) + γ c 1 2 c 0 2 ϕ ( r ) + c 0 4 ν ϕ ( r ) 3 c 1 c 2 c 0 2 ρ ϕ ( r ) c 1 c 0 3 σ ϕ ( r ) 6 α c 1 3 c 0 ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 δ ϕ ( r ) ϕ ( r ) + 3 c 1 2 c 0 2 ρ ϕ ( r ) ϕ ( r ) = 0 .
The solution of (47) gives us complex trigonometric function satisfying (4) as:
u ( t , x , y ) = 2 a 2 α 9 δ ρ t 3 + a 3 ρ x 2 a α 3 δ ρ y + c 2 c 0 x ± 2 Ω 3 i c 1 3 ( α c 1 β c 0 ) c 0 Ω 3 c 1 2 ( 2 α c 1 c 0 ( ρ + δ ) ) × tan { Ω 3 2 c 1 3 ( α c 1 β c 0 ) [ β c 0 c 1 3 A 1 α c 1 4 A 1 c 0 c 1 3 ( α c 1 β c 0 ) × y a 2 t 2 c 1 c 0 x ] } + A 2 ,
where Ω 3 = ν c 0 3 + c 0 c 1 2 ( a ( t 1 ) σ ) + 6 α c 1 2 c 2 + c 0 c 1 ( γ c 1 3 ρ c 2 ) with A 1 and A 2 representing the integration constants of the solution.

2.3.5. Optimal Subalgebra X 2 + a X 3 + b X 5 , a , b { 1 , 1 }

We reduce (4) via X 2 + a X 3 + b X 5 to a NLPDE with dependent variables X , Y as:
3 a γ ρ G X X + 3 a σ ρ G X Y + 3 a ν ρ G Y Y 3 ρ G X Y + 9 a ρ 2 G X G X Y + 18 a α ρ G X G X X + 9 a δ ρ G Y G X X + 3 a α ρ G X X X X + 3 a β ρ G X X X Y + b = 0 ,
by utilizing the invariants with their group-invariant expressed via the function
X = x , Y = 1 2 a 2 a y b t 2 2 t , where we calculated the group - invariant as u ( t , x , y ) = G ( X , Y ) + 2 b 2 α 9 a 2 δ ρ t 3 + b α 3 a 2 δ ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t .
On applying Lie symmetry algorithm to Equation (49), we achieve three generators
R 1 = X , R 2 = Y , R 3 = G .
Similarity solution to R 1 = / X yields G ( X , Y ) = ϕ ( r ) , where r = Y . Therefore using the function reduces (4) to the linear ordinary differential equation (LODE)
3 a ρ ν ϕ ( r ) + b = 0 .
The solution to the LODE is ϕ ( r ) = b r 2 / 6 a ν ρ + A 1 r + A 2 , where A 1 and A 2 are integration constants. Hence, solution to (2+1)-D genBKe (4) in this regard is:
u ( t , x , y ) = 2 b 2 α 9 a 2 δ ρ t 3 + b α 3 a 2 δ ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t b 24 a 3 ν ρ 2 a y b t 2 2 t 2 + A 1 2 a 2 a y b t 2 2 t + A 2 .
In the same vein, generator R 2 furnishes G ( X , Y ) = ϕ ( r ) , r = X , so (4) becomes:
3 a α ρ ϕ ( r ) + 3 a γ ρ ϕ ( r ) + 18 a α ρ ϕ ( r ) ϕ ( r ) + b = 0 .
No solution of (52) can be secured. However, considering a special case of the equation with b = 0 , one achieves a trigonometric solution of (4) in this regard as
u ( t , x , y ) = 2 b 2 α 9 a 2 δ ρ t 3 + b α 3 a 2 δ ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t γ α tan γ 4 α x α A 1 + A 2 ,
which is actually an algebraic-trigonometric solution of (2+1)-D genBKe (4). Further, imploring generators R 1 and R 2 , we obtain solution function G ( X , Y ) = ϕ ( r ) , r = Y c 1 / c 0 X . On applying the function in Equation (4) changes it to NODE
3 a α c 1 4 ρ ϕ ( 4 ) ( r ) 3 a β c 1 3 c 0 ρ ϕ ( 4 ) ( r ) + 3 a γ c 1 2 c 0 2 ρ ϕ ( r ) + 3 a c 0 4 ν ρ ϕ ( r ) 3 a c 1 c 0 3 ρ σ ϕ ( r ) 18 a α c 1 3 c 0 ρ ϕ ( r ) ϕ ( r ) + 9 a c 1 2 c 0 2 δ ρ ϕ ( r ) ϕ ( r ) + 9 a c 1 2 c 0 2 ρ 2 ϕ ( r ) ϕ ( r ) + b c 0 4 + 3 c 1 c 0 3 ρ ϕ ( r ) = 0 .
We let b = 0 to gain an elliptic solution of (54) and give it a simple representation:
α 0 ϕ ( r ) + α 1 ϕ ( r ) ϕ ( r ) + α 2 ϕ ( 4 ) ( r ) = 0
where α 0 = 3 a γ c 1 2 c 0 2 ρ + 3 a c 0 4 ν ρ 3 a c 1 c 0 3 ρ σ + 3 c 1 c 0 3 ρ , α 1 = 18 a α c 0 c 1 3 ρ + 9 a c 0 2 c 1 2 δ ρ + 9 a c 0 2 c 1 2 ρ 2 , α 2 = 3 a α c 1 4 ρ 3 a β c 0 c 1 3 ρ . Integrating (55) twice with ϕ ( r ) = Θ ( r ) gives
Θ ( r ) 2 = α 1 3 α 2 Θ ( r ) 3 α 0 α 2 Θ ( r ) 2 2 A 0 α 2 Θ ( r ) 2 A 1 α 2 ,
where A 0 and A 1 are integration constants. We engage the transformation,
Θ ( r ) = 12 α 2 α 1 ( r ) α 0 α 1 .
Thus, we reckon Equation (56) as NODE with Weierstrass elliptic function [58,59]
( r ) 2 4 ( r ) 3 + g 1 ( r ) + g 2 = 0 ,
with the involved Weierstrass elliptic invariants g 1 and g 2 expressed as:
g 1 = 1 12 α 2 2 α 0 2 2 α 1 A 0 , and g 2 = 1 216 α 2 3 α 0 3 + 3 α 1 α 1 A 1 α 0 A 0 .
Contemplating (57) alongside (58) and reverting to the basic variables yields:
u ( t , x , y ) = 2 b 2 α 9 a 2 δ ρ t 3 + b α 3 a 2 δ ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t α 0 2 a α 1 2 a y b t 2 2 t + α 0 c 1 α 1 c 0 x + 12 α 2 α 1 ζ { 1 2 a 2 a y b t 2 2 t c 1 c 0 x ; 1 12 α 2 2 α 0 2 2 α 1 A 0 , 1 216 α 2 3 α 0 3 + 3 α 1 α 1 A 1 α 0 A 0 } .
Next, we consider the combination of obtained symmetries as Q = c 0 / X + c 1 / Y + c 2 / G . Consequently, Q gives the function G ( X , Y ) = c 2 / c 0 X + ϕ ( r ) , where r = Y c 1 / c 0 X . Invoking the function and its variables, we reduce (4) to:
3 a α c 1 4 ρ ϕ ( 4 ) ( r ) 3 a β c 1 3 c 0 ρ ϕ ( 4 ) ( r ) + 18 a α c 1 2 c 2 c 0 ρ ϕ ( r ) + 3 a γ c 1 2 c 0 2 ρ ϕ ( r ) + 3 a c 0 4 ν ρ ϕ ( r ) 9 a c 1 c 2 c 0 2 ρ 2 ϕ ( r ) 3 a c 1 c 0 3 ρ σ ϕ ( r ) 18 a α c 1 3 c 0 ρ ϕ ( r ) ϕ ( r ) + 9 a c 1 2 c 0 2 δ ρ ϕ ( r ) ϕ ( r ) + 9 a c 1 2 c 0 2 ρ 2 ϕ ( r ) ϕ ( r ) + b c 0 4 + 3 c 1 c 0 3 ρ ϕ ( r ) = 0 .
Just as earlier demonstrated, we present simplified structure of (61) with b = 0 as:
α 5 ϕ ( 4 ) ( r ) + 6 α 4 ϕ ( r ) ϕ ( r ) α 3 ϕ ( r ) = 0 ,
where α 3 = 9 a c 1 c 2 c 0 2 ρ 2 + 3 a c 1 c 0 3 ρ σ 18 a α c 1 2 c 2 c 0 ρ 3 a γ c 1 2 c 0 2 ρ 3 c 1 c 0 3 ρ 3 a c 0 4 ν ρ , α 4 = 3 a α c 1 3 c 0 ρ + 3 / 2 a c 1 2 c 0 2 δ ρ + 3 / 2 a c 1 2 c 0 2 ρ 2 , α 5 = 3 a α c 1 4 ρ 3 a β c 1 3 c 0 ρ . On Integrating (62)
α 5 ϕ ( r ) + 3 α 4 ϕ ( r ) 2 α 3 ϕ ( r ) + K 0 = 0 ,
with integration constant K 0 . On engaging the representations expressed as:
ϕ ( r ) = α 5 α 4 Θ ( r ) , λ = α 3 α 5 , K 1 = K 0 α 4 α 5 2 ,
Equation (63) then becomes the second order nonlinear differential equation:
Θ ( r ) + 3 Θ ( r ) 2 λ Θ ( r ) + K 1 = 0
Equation (65) multiplied by Θ ( r ) and integrating the outcome furnishes,
Θ ( r ) 2 = ( 2 Θ ( r ) 3 λ Θ ( r ) 2 + 2 K 1 Θ ( r ) + 2 K 2 ) ,
with integration constant K 2 . Suppose that the algebraic equation Θ ( r ) 3 1 2 λ Θ ( r ) 2 + K 1 Θ ( r ) + K 2 = 0 possesses roots ϑ 1 , ϑ 2 , ϑ 3 with the property ϑ 1 > ϑ 2 > ϑ 3 , then
Θ ( r ) 2 = 2 ( Θ ( r ) ϑ 1 ) ( Θ ( r ) ϑ 2 ) ( Θ ( r ) ϑ 3 ) .
Equation (66) possess a highly famous solution expressed with regards to Jacobi elliptic function (cn) [58,60] which we present in the structure,
Θ ( r ) = ϑ 2 + ( ϑ 1 ϑ 2 ) cn 2 ϑ 1 ϑ 3 2 r | Δ 2 , where Δ 2 = ϑ 1 ϑ 2 ϑ 1 ϑ 3 .
Reckoning (67) as well as (64) and retrograding to the basic variables gives:
u ( t , x , y ) = 2 b 2 α 9 a 2 δ ρ t 3 + b α 3 a 2 δ ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t + α 5 α 4 { r ( ϑ 2 + ϑ 1 ( Δ 2 1 ) ) Δ 2 + 2 ϑ 1 ϑ 2 dn ϑ 1 ϑ 3 2 r | Δ 2 E am ϑ 1 ϑ 3 2 r | Δ 2 | Δ 2 ϑ 1 ϑ 3 Δ 2 dn ϑ 1 ϑ 3 2 r | Δ 2 2 } + c 2 c 0 x ,
with E representing elliptic integral of the second kind while ‘am’ and ‘dn’ are respectively amplitude and delta elliptic functions. Besides, we notice that in relation (67) and (68) some limits of Jacobi elliptic functions cn and dn exist which give rise to some other functions such as hyperbolic and trigonometric. For instance, lim Δ 2 0 cn r | Δ 2 = cos ( r ) , lim Δ 2 0 dn r | Δ 2 = 1 , lim Δ 2 1 cn r | Δ 2 = sec h ( r ) and lim Δ 2 1 dn r | Δ 2 = sec h ( r ) , whereas r = 1 / 2 a 2 a y b t 2 2 t c 1 / c 0 x .

2.3.6. Optimal Subalgebra X 1 + a X 3 + b X 5 , a , b { 1 , 1 }

Lie optimal subalgebra X 1 + a X 3 + b X 5 produces similarity transformation variables,
X = 1 a a x t , Y = 1 2 a 2 a y b t 2 , whereas the group - invariant is sec ured as u ( t , x , y ) = G ( X , Y ) + 2 b 2 α 9 a 2 δ ρ t 3 b 6 a 2 ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t .
Engaging the found similarity variables reduces (2+1)-D genBKe (4) to an NLPDE
3 a γ ρ G X X + 3 a σ ρ G X Y + 3 a ν ρ G Y Y 3 ρ G X X + 9 a ρ 2 G X G X Y + 18 a α ρ G X G X X + 9 a δ ρ G Y G X X + 3 a α ρ G X X X X + 3 a β ρ G X X X Y + b = 0 ,
The Lie theoretic approach used in studying Equation (69) yields its symmetries as:
R 1 = X , R 2 = Y , R 3 = G .
On following the usual process solution to R 1 = / X secures G ( X , Y ) = ϕ ( r ) , with r = Y . Subsequently utilizing the function obtained reduces (4) to the LODE,
3 a ρ ν ϕ ( r ) + b = 0 .
On solving the linear ordinary differential Equation (70), we obtain a solution of (4) as:
u ( t , x , y ) = A 0 2 a 2 a y b t 2 + 2 b 2 α 9 a 2 δ ρ t 3 b 6 a 2 ρ t 2 + b 3 a ρ x 2 b α 3 a δ ρ y t b 24 a 3 ν ρ 2 a y b t 2 2 t 2 + A 1 ,
with integration constants A 0 and A 1 . In addition R 2 gives the solution G ( X , Y ) = ϕ ( r ) , with r = X . On engaging the function secured, we reduce (4) to the LODE,
3 a α ρ ϕ ( r ) + 3 a γ ρ ϕ ( r ) + 18 a α ρ ϕ ( r ) ϕ ( r ) 3 ρ ϕ ( r ) + b = 0 .
In a bid to secure a solution of (4) in this instance, we let b = 0 and, as a consequence:
u ( t , x , y ) = 2 α b 2 t 3 9 a 2 δ ρ b t 2 6 a 2 ρ + b x 3 a ρ 2 α b y 3 a δ ρ ( a γ 1 ) a α ( 1 a γ ) a α ( 1 a γ ) × tanh 1 2 1 a γ a α a x t a ± a α C 1 + C 2 ,
which is an algebraic–hyperbolic solution of (2+1)-D genBKe (4) with integration constants C 1 and C 2 . On following the usual procedure, R 1 and R 2 linearly combined yields the solution G ( X , Y ) = ϕ ( r ) , r = c 0 Y c 1 X and these transform (4) to:
3 a α c 1 4 ρ ϕ ( 4 ) ( r ) 3 a β c 0 c 1 3 ρ ϕ ( 4 ) ( r ) + 3 a γ c 1 2 ρ ϕ ( r ) + 3 a c 0 2 ν ρ ϕ ( r ) 3 a c 0 c 1 ρ σ ϕ ( r ) + b 18 a α c 1 3 ρ ϕ ( r ) ϕ ( r ) + 9 a c 0 c 1 2 δ ρ ϕ ( r ) ϕ ( r ) + 9 a c 0 c 1 2 ρ 2 ϕ ( r ) ϕ ( r ) 3 c 1 2 ρ ϕ ( r ) = 0 .
Now, having observed that no solution of (73) can be secured in its current state, we take a special case b = 0 of the equation. We present in an easier way (73) as:
β 2 ϕ ( 4 ) ( r ) β 1 ϕ ( r ) ϕ ( r ) β 0 ϕ ( r ) = 0 ,
where β 0 = 3 a c 1 c 0 ρ σ + 3 c 1 2 ρ 3 a γ c 1 2 ρ 3 a c 0 2 ν ρ , β 1 = 18 a α c 1 3 ρ 9 a c 0 c 1 2 δ ρ 9 a c 0 c 1 2 ρ 2 , β 2 = 3 a α c 1 4 ρ 3 a β c 0 c 1 3 ρ . We let ϕ ( r ) = Θ ( r ) in (74) and integrating the equation gives
2 β 2 Θ ( r ) β 1 Θ ( r ) 2 2 β 0 Θ ( r ) = 2 C 0 ,
where C 0 is the integration constant. On taking the multiplication of (75) and Θ ( r ) and subsequently integrating the resulting NODE, one then achieves:
Θ ( r ) 2 = β 1 3 β 2 Θ ( r ) 3 + β 0 β 2 Θ ( r ) 2 + 2 C 0 β 2 Θ ( r ) + 2 C 1 β 2 ,
with integration constant C 1 . We get a Weierstrass elliptic solution [61] of (4) via:
Θ ( r ) = W ( r ) β 0 β 1 ,
which is the transformation needed in this regard to reduce (76) to elliptic function,
W ξ 2 = 4 W 3 g 2 W g 3 , where ξ = β 1 12 β 2 r .
That is, a Weierstrass elliptic function with elliptic invariants g 1 and g 2 secured as:
g 1 = 24 C 0 β 1 12 β 0 2 β 1 2 , and g 2 = 8 β 0 3 β 1 3 24 β 0 C 0 β 1 2 + 24 C 1 β 1 .
On reckoning (77), we possess the solution of (76) with regards to Θ ( r ) as:
Θ ( r ) = β 1 12 β 2 ( r r 0 ) ; 24 C 0 β 1 12 β 0 2 β 1 2 , 8 β 0 3 β 1 3 24 C 0 β 0 β 1 2 + 24 C 1 β 1 β 0 β 1 .
On reverting to the basic variables, one achieves the solution of Equation (4) as:
u ( t , x , y ) = 2 α b 2 t 3 9 a 2 δ ρ b t 2 6 a 2 ρ + b x 3 a ρ 2 α b y 3 a δ ρ t 2 3 β 2 β 1 ζ { 1 2 β 1 3 β 2 × c 0 2 a 2 a y b t 2 c 1 a ( a x t ) r 0 ; 24 C 0 β 1 12 β 0 2 β 1 2 , 8 β 0 3 β 1 3 24 C 0 β 0 β 1 2 + 24 C 1 β 1 } β 0 β 1 c 0 2 a 2 a y b t 2 c 1 a ( a x t ) r 0 ,
which is a Weierstrass elliptic solution of (4) where r 0 is an arbitrary constant. Next, we contemplate the combination of the three found symmetries as performed earlier, and secure G ( X , Y ) = c 2 X + c 0 ϕ ( r ) , with r = c 0 Y c 1 X which transform (4) to:
3 a α c 0 c 1 4 ρ ϕ ( 4 ) ( r ) 3 a β c 0 2 c 1 3 ρ ϕ ( 4 ) ( r ) + 18 a α c 0 c 2 c 1 2 ρ ϕ ( r ) + 3 a γ c 0 c 1 2 ρ ϕ ( r ) + b 9 a c 0 2 c 2 c 1 ρ 2 ϕ ( r ) 3 a c 0 2 c 1 ρ σ ϕ ( r ) 18 a α c 0 2 c 1 3 ρ ϕ ( r ) ϕ ( r ) + 9 a c 0 3 c 1 2 δ ρ ϕ ( r ) ϕ ( r ) + 9 a c 0 3 c 1 2 ρ 2 ϕ ( r ) ϕ ( r ) 3 c 0 c 1 2 ρ ϕ ( r ) + 3 a c 0 3 ν ρ ϕ ( r ) = 0 .
In order to gain more general solution of (4) in this regard, we let b = 0 and so:
β 3 ϕ ( 4 ) ( r ) + 12 β 4 ϕ ( r ) ϕ ( r ) β 5 ϕ ( r ) = 0 ,
where β 3 = 3 a α c 0 c 1 4 ρ 3 a β c 0 2 c 1 3 ρ , β 4 = 3 / 4 a c 0 3 c 1 2 δ ρ + 3 / 4 a c 0 3 c 1 2 ρ 2 3 / 2 a α c 0 2 c 1 3 ρ , β 5 = 9 a c 0 2 c 2 c 1 ρ 2 18 a α c 0 c 2 c 1 2 ρ 3 a γ c 0 c 1 2 ρ + 3 a c 0 2 c 1 ρ σ + 3 c 0 c 1 2 ρ 3 a c 0 3 ν ρ . On the integration of Equation (81) and invoking the representation ϕ ( r ) = β 3 / 2 β 4 Θ ( r ) , we obtain:
Θ ( r ) + 3 Θ ( r ) 2 ω Θ ( r ) + A 1 = 0 ,
where ω = β 5 / β 3 with A 1 = 2 β 4 A 0 / β 3 2 , A 0 and A 1 being integration constants. Next, we multiply (82) by Θ ( r ) and integrate the result with regards to r and secure
Θ ( r ) 2 + 2 Θ ( r ) 3 ω Θ ( r ) 2 + 2 A 1 Θ ( r ) + 2 A 2 = 0 .
Thus, (83) occasions a well notable Jacobi elliptic cosine function solution [61] with cubic polynomial roots θ 3 < θ 2 < θ 1 and besides, parameter 0 Ω 0 2 1 . In consequence, we recover u ( t , x , y ) , the solution of Equation (4) in this instance as:
u ( t , x , y ) = 2 α b 2 t 3 9 a 2 δ ρ + b x 3 a ρ t b t 2 6 a 2 ρ 2 α b y 3 a δ ρ t + c 2 a a x t + θ 2 r + c 0 β 3 2 β 4 2 θ 1 θ 2 sn θ 1 θ 3 2 r | Ω 0 2 cos 1 dn θ 1 θ 3 2 r | Ω 0 2 | Ω 0 2 θ 1 θ 3 1 dn θ 1 θ 3 2 r | Ω 0 2 2 ,
where Ω 0 2 = ( θ 1 θ 2 ) / ( θ 1 θ 3 ) and r = c 0 / 2 a ( 2 a y b t 2 ) c 1 / a ( a x t ) . Moreover, the Jacobi sine elliptic function sn possesses the property that as Ω 0 2 0 , we have sn ( r ) sin ( r ) and as Ω 0 2 1 , we also obtain sn ( r ) tanh ( r ) .

2.3.7. Optimal Subalgebra X 1 + a X 2 + b X 3 , a , b { 1 , 1 }

The Lagrangian system related to X 1 + a X 2 + b X 3 solves to give group-invariant
u ( t , x , y ) = G ( X , Y ) , where X = x t / b , Y = y a t / b .
On using the function alongside other expressions from (85) in (4), we have:
b γ G X X + b σ G X Y + b ν G Y Y a G X Y G X X + 3 b ρ G X G X Y + 6 b α G X G X X + 3 b δ G Y G X X + b α G X X X X + b β G X X X Y = 0 .
As a consequence, we secure the solution of (86) with respect to X and Y but reverting to the fundamental variables gives a solution of (2+1)-D genBKe (4) as:
u ( t , x , y ) = 4 i β Ω 0 + α β + b ( 2 α ν β σ 4 β 2 A 1 2 ) Ω 0 ( δ + ρ ) + Ω 1 { A 1 sec h [ 1 2 b 2 ν ( a 2 A 1 t + Ω 0 A 1 ( a t b y ) + b 2 ( A 1 ( σ + 4 β A 1 2 ) y ν ( 2 A 0 + 2 A 1 x ) ) b A 1 ( a y + t ( a ( σ + 4 β A 1 2 ) 2 ν ) ) ) ] } { A 2 ( Ω 0 ( δ + ρ ) + a ( δ + ρ ) + b ( 4 α ν δ σ ρ σ 4 β A 1 2 ( δ + ρ ) ) ) + sec h [ 1 2 b 2 ν ( a 2 A 1 t + Ω 0 A 1 ( a t b y ) + b 2 4 A 1 y 1 4 σ + β A 1 2 ν ( 2 A 1 x + i π + 2 A 0 ) 4 A 1 b [ 1 4 a y + t a β A 1 2 + 1 4 σ 1 2 ν ] ) ] } 1 ,
where Ω 0 = ν ( 4 b 4 b 2 γ 16 α b 2 A 1 2 ) + ( a b σ 4 b β A 1 2 ) , Ω 1 = a ( δ + ρ ) + b ( 4 α ν δ σ ρ σ 4 β ( δ + ρ ) A 1 2 ) with constants A 0 and A 1 arbitrary. Function (87) is a complex bright soliton solution of (4). Furthermore, investigation revealed that Equation (86) possesses three Lie point symmetries which are given as
R 1 = X , R 2 = Y , R 3 = G .
Linearly combining the symmetries furnishes the function G ( X , Y ) = c 2 X + c 0 ϕ ( r ) , with r = c 0 Y c 1 X . Thus, on engaging the function, we further reduce (4) to:
b β c 0 c 1 3 ϕ ( 4 ) ( r ) a c 0 c 1 ϕ ( r ) α b c 1 4 ϕ ( 4 ) ( r ) 6 α b c 2 c 1 2 ϕ ( r ) b γ c 1 2 ϕ ( r ) b c 0 2 ν ϕ ( r ) + 3 b c 0 c 2 c 1 ρ ϕ ( r ) + b c 0 c 1 σ ϕ ( r ) + 6 α b c 0 c 1 3 ϕ ( r ) ϕ ( r ) 3 b c 0 2 c 1 2 δ ϕ ( r ) ϕ ( r ) 3 b c 0 2 c 1 2 ρ ϕ ( r ) ϕ ( r ) + c 1 2 ϕ ( r ) = 0 .
Therefore, we present Equation (88) in a lesser structure as:
α 1 ϕ ( r ) α 2 ϕ ( r ) ϕ ( r ) + α 3 ϕ ( 4 ) ( r ) = 0 ,
α 1 = a c 0 c 1 + c 1 2 6 α b c 2 c 1 2 b γ c 1 2 b c 0 2 ν + 3 b c 0 c 2 c 1 ρ + b c 0 c 1 σ , α 2 = 3 b c 0 2 c 1 2 δ + 3 b c 0 2 c 1 2 ρ 6 α b c 0 c 1 3 , α 3 = b β c 0 c 1 3 α b c 1 4 . We set ϕ ( r ) = Θ ( r ) in (89) and by integrating the resulting NODE repeatedly two times, we secure a first order NODE presented as:
Θ ( r ) 2 = α 2 3 α 3 Θ ( r ) 3 α 1 α 3 Θ ( r ) 2 2 C 0 α 3 Θ ( r ) 2 C 1 α 3 ,
with constants of integration C 0 and C 1 . On contemplating the cubic polynomial α 2 3 α 3 Θ ( r ) 3 α 1 α 3 Θ ( r ) 2 2 C 0 α 3 Θ ( r ) 2 C 1 α 3 = 0 , whose real roots are a 2 < a 1 < a 0 , we have
Θ r 2 = α 2 3 α 3 ( Θ a 0 ) ( Θ a 1 ) ( Θ a 2 ) ,
with real roots a 0 , a 1 as well as a 2 satisfying algebraic relations expressed as:
a 0 a 1 + a 0 a 2 + a 1 a 2 = 2 C 0 α 3 , a 0 a 1 a 2 = 2 C 1 α 3 , a 0 + a 1 + a 2 = α 1 α 3
According to [62], we express a primitive solution of (4) via the elliptic function,
u ( t , x , y ) = c 2 x c 2 b t + c 0 { 12 α 3 ( a 0 a 1 ) 2 α 2 ( a 0 a 2 ) Δ 0 8 { EllipticE [ s n ( α 2 ( a 0 a 2 ) 12 α 3 ( r r 0 ) , Δ 0 2 ) , Δ 0 2 ] } + a 1 ( a 0 a 1 ) 1 Δ 0 4 Δ 0 4 ( r r 0 ) + C 2 } ,
with r = c 0 ( y a t / b ) c 1 ( x t / b ) , r 0 and C 2 arbitrary constants. Besides, parameter Δ 0 2 and incomplete elliptic integral EllipticE [ m ; z ] are accordingly expressed as:
Δ 0 2 = a 0 a 1 a 0 a 2 and EllipticE [ m ; z ] = 0 m 1 z 2 w 2 1 w 2 d w .

3. Travelling Wave Solutions

We examine the travelling wave solutions of the (2+1)-D genBKe (4). Generally speaking, travelling wave solutions of a partial differential equation emanates as special group-invariant solutions wherein the considered group is translational with respect to space of independent variables.
Here in this study, we engage linear combination of the translation operators X 1 , X 2 and X 3 , namely X = ρ X 1 + ε X 2 + μ X 3 with constant values σ and ε . Following the usual Lie symmetry procedure, we utilize X to reduce (4) to fourth-order NODE,
A ψ ( z ) B ψ ( z ) ψ ( z ) + C ψ ( z ) = 0 ,
via the travelling wave z = p x + q y + r t where p = ε , q = μ c ρ , r = ε c and so A = p ( r + σ q + γ p ) + ν q 2 , B = 6 p 2 ( α p + β q ) and C = p 3 ( α p + β q ) .
Integrating (91) just once supplies a third-order ODE,
A ψ 1 2 B ψ 2 + C ψ + C 1 = 0 ,
where C 1 is regarded as an integration constant. Multiplying Equation (92) by ψ , integrating once as well as simplifying the resulting equation, we have the second-order nonlinear ODE
1 2 A ( ψ ) 2 1 6 B ( ψ ) 3 + 1 2 C ( ψ ) 2 + C 1 ψ + C 2 = 0 ,
where C 2 is an integration constant. Equation (93) can be rewritten as
( ψ ) 2 = B 3 C ( ψ ) 3 A C ( ψ ) 2 2 C 1 C ψ 2 C 2 C .
Suppose Ψ = ψ , Equation (94) becomes:
Ψ 2 = B 3 C Ψ 3 A C Ψ 2 2 C 1 C Ψ 2 C 2 C .

3.1. Bifurcation and Explicit Solutions

Here we use the bifurcation theory method [39,63,64] of dynamical systems to obtain some nontrivial solutions of (95), which is the reduced form of (91).
Suppose from Equation (95) we say:
P 3 ( Ψ ) = B 3 C Ψ 3 A C Ψ 2 2 C 1 C Ψ 2 C 2 C .
We can deduce from Equation (94) that:
d 2 Ψ d z 2 = B 2 C Ψ 2 A C Ψ C 1 C .
Let Ψ = w , then (95) is equivalent to planar dynamical system,
d Ψ d z = w , d w d Ψ = B 2 C Ψ 2 A C Ψ C 1 C ,
which invariably possesses the first integral H ( Ψ , w ) calculated as:
H ( Ψ , w ) = w 2 2 B 6 C Ψ 3 + A 2 C Ψ 2 + C 1 C Ψ = h ,
where h is the constant of integration and function H ( Ψ , w ) is Hamiltonian.
It is obvious to see that Hamiltonian H ( Ψ , w ) = h = C 2 C corresponds to Equation (96). As a result, we observe that the dynamical system behaviors of ordinary differential Equation (95) from the orbits of the above system (98) relates to H ( Ψ , w ) = C 2 C . Apparently, phase orbits given via the vector field relative to system (98) decides all the results that can be gained for (96).
An investigation of bifurcation of the planar dynamical system (98) secures diverse kinds of solutions of (96) contemplated under various coefficient conditions. Thus, the dynamical character and closed-form solutions of ODE (96) are generated.
We first study the equilibrium points of the system (98) to attain the dynamical action of the system. Evidently, the roots of P 3 ( Ψ ) = 0 are regarded as the abscissas of the points of equilibrium included in the system (98). Moreover, we suppose that Ψ e is one of the roots of P 3 ( Ψ ) = 0 , meaning that, ( Ψ e , 0 ) stands as an equilibrium point of system (98). By the reason of theory of planar dynamical systems [63,64], the matrix needs to be studied.
D f ( Ψ 0 , 0 ) = 0 1 P 3 ( Ψ e ) 0 ,
where
P 3 ( Ψ e ) = 2 B C Ψ 2 A C
of the linearized system of (98) exists at a point ( Ψ e , 0 ) . The point of equilibrium ( Ψ , 0 ) is a center which has a punctured neighborhood wherein any solution procured is taken as a periodic orbit; if d e t ( D f ( Ψ e , 0 ) ) = P 3 ( Ψ e ) > 0 . It is said to be a saddle point if d e t ( D f ( Ψ e , 0 ) ) = P 3 ( Ψ e ) < 0 . Nevertheless, we call it a cusp point if d e t ( D f ( Ψ e , 0 ) ) = P 3 ( Ψ e ) = 0 . It is needed to equally investigate boundary curves related to the centers as well as the orbits that serve as a connector between the saddle points or cusp points which the Hamiltonian H ( Ψ , w ) = h determines in order to obtain the phase portraits other than the equilibriums. Evidently, system (98) possesses neither equilibrium point nor a cusp when A 2 + 2 B C 1 C 2 0 , hence system (98) has no trivial nontrivial bounded solutions. Nonetheless, (98) has two equilibrium points when A 2 + 2 B C 1 C 2 > 0 . Let
Ψ e ± = 1 B A ± A 2 + 2 B C 1
H ( Ψ e ± , 0 ) = h ± = 1 3 B 2 C ( A 2 + 2 B C 1 ) [ A ± A 2 + 2 B C 1 ] + A B C 1 ,
then ( Ψ e + , 0 ) is a saddle point, ( Ψ e , 0 ) is also a center and h + > h .
When we have h + > h > h , Hamiltonian H ( Ψ , w ) = h defines a family of periodic orbits present around the center given as ( u e , 0 ) which is confined by the boundary curves defined by function H ( Ψ , w ) = h + . Notwithstanding, H ( Ψ , w ) = h + explains a homoclinic orbit that passes through the saddle point ( Ψ e + , 0 ) .
We now consider some cases of (96) and obtain the following solutions.
Case (1.) Equation (96) possesses a bounded solution which approaches Ψ e + as z goes to infinity:
Ψ ( z ) = 1 B ( A + A 2 + 2 B C 1 ) 3 A 2 + 2 B C 1 sech 2 1 2 A 2 + 2 B C 1 C 2 4 ( z z 0 ) ,
where z 0 is an arbitrary constant. Integrating (100) and returning to the original variables secures a nontrivial solitary wave solution of (4) in this regard as:
u ( t , x , y ) = 1 B { A 2 + 2 B C 1 + A z z 0 6 A 2 + 2 B C 1 A 2 + 2 B C 1 C 2 4 × tanh 1 2 A 2 + 2 B C 1 C 2 4 z z 0 } ,
with z = p x + q y + r t and z 0 arbitrary constant. We also have a constant solution
Ψ ( z ) = 1 B A + A 2 + 2 B C 1
as well as an unbounded solution:
Ψ ( z ) = 1 B ( A + A 2 + 2 B C 1 ) + 3 A 2 + 2 B C 1 c s c h 2 1 2 A 2 + 2 B C 1 C 2 4 ( z z 0 ) .
Integrating (102) and retrograding to the original variables, we secure an unbounded solution of (2+1)-dimensional gBK (4) as:
u ( t , x , y ) = 1 B { A 2 + 2 B C 1 + A z z 0 6 A 2 + 2 B C 1 A 2 + 2 B C 1 C 2 4 × coth 1 2 A 2 + 2 B C 1 C 2 4 z z 0 } ,
where z = p x + q y + r t and z 0 is an arbitrary constant.
Case (2.) Since B 3 C > 0 , then for any arbitrary real constant
Φ A 2 A 2 + 2 B C 1 B , A A 2 + 2 B C 1 B ,
Ψ ( z ) = Φ 1 2 3 Φ 3 A B + 3 Φ 2 + 6 A B Φ + 9 A 2 B 2 + 24 C 1 B sn 2 ( Ω + ( z z 0 ) , k + ) ,
where
Ω + = 2 4 B C Φ + A C + B 3 C 3 Φ 2 + 6 A B Φ + 9 A 2 B 2 + 24 C 1 B and
k + = 2 3 Φ 2 6 A B Φ 6 C 1 B 3 Φ + 3 A B + 3 Φ 2 + 6 A B Φ + 9 A 2 B 2 + 24 C 1 B .
The integration of (104) secures a bounded nontrivial solution of (4) as:
u ( t , x , y ) = k z + P 0 6 3 A B P 1 Q + E am R | R 1 sn 2 ( S | S 1 ) S 2 dn Q 1 | Q 2 1 2 3 Q 0 sn 2 ( R ) P 1 3 A B
P 0 = 1 B 3 A 3 B k + 3 B 3 A 2 + 2 B k A B 2 k 2 + 8 B C 1 B 2 C , P 1 = 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 Q = z 3 A + B 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 3 k C 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 12 2 k ( B k 2 A ) 2 C 1 B , R = z 3 A + B 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 3 k C 2 6 | 2 3 Q 0 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 S = z 3 A + B 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 3 k C 2 6 , R 1 = 2 3 k ( B k 2 A ) 2 C 1 B 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 , S 1 = 2 3 k ( B k 2 A ) 2 C 1 B 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 , Q 1 = z 3 A + B 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 3 k C 2 6 , S 2 = 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 2 3 k ( B k 2 A ) 2 C 1 B , Q 2 = 2 3 k ( B k 2 A ) 2 C 1 B 3 A B 3 k + 9 A 2 + 6 B k A 3 B 2 k 2 + 24 B C 1 B 2 ,
where E [ am ( R | R 1 ) ] is an elliptic integral of the second kind sn ( S | S 1 ) , am ( R | R 1 ) and dn ( Q 1 | Q 2 ) denotes accordingly elliptic sine, amplitude as well as delta functions. In addition to that, variable z = p x + q y + r t with arbitrary constant z 0 is taken as zero.
Case (3.) Equation (96) possesses no nontrivial bounded solutions. However, at the instance when 2 C 2 C = 2 h , we have an unbounded solution that is expressed as
Ψ ( z ) = 1 B ( A A 2 + 2 B C 1 ) + 3 A 2 + 2 B C 1 sec 2 1 2 A 2 + 2 B C 1 C 2 4 ( z z 0 )
and a constant solution also given in this case as:
Ψ ( z ) = 1 B A A 2 + 2 B C 1 .
Integrating (105) with regards to variable z z 0 , one achieves:
u ( t , x , y ) = 1 B V 0 tan 1 2 ( z z 0 ) A 2 + 2 B C 1 C 2 4 + ( z z 0 ) A A 2 + 2 B C 1 ,
where we have V 0 = 6 A 2 + 2 B C 1 A 2 + 2 B C 1 C 2 4 , z = p x + q y + r t and z 0 as an arbitrary constant.
We note from the dynamical system earlier stated that we can deduce the fact that:
d w d Ψ = B 0 Ψ 2 + B 1 Ψ + B 2 ,
where B 0 = B 2 C , B 1 = A C and B 2 = C 1 C . In clear terms, we can suggest that phase orbits given by the vector fields of dynamical system (98) determined the collection of all the solutions of (97). Thus, we state here that bounded solutions of (97) relates to the bounded phase orbits that system (98) has which will have to be investigated. Along the orbit connected with H ( Ψ , w ) = h , we have:
d Ψ d z 2 = 2 B 0 3 Ψ 3 + B 1 Ψ 2 + 2 B 2 Ψ + 2 h .
As a consequence, the general formula associated with the solutions of (97) can as well be given viz;
0 Ψ d Ψ ( 2 B 0 Ψ 3 / 3 ) + B 1 Ψ 2 + 2 B 2 Ψ + 2 h = ± z 0 z d z .
Nonetheless, it may be laborious to know the properties as well as the shapes of (109) that are actually decided by the parameters B 0 , B 1 , B 2 and h. Obviously, the abscissas possessed by equilibrium points of dynamical system (98) are zeros of B 0 Ψ 2 + B 1 Ψ + B 2 = 0 . Clearly, the system (98) has no bounded orbits when B 1 2 4 B 0 B 2 < 0 . We suppose that B 1 2 4 B 0 B 2 > 0 in order for us to examine the bounded orbits owned by system (98). We designate Ψ ± = ( B 1 ± B 1 2 4 B 0 B 2 ) / 2 B 0 , and as such we have E + ( Ψ + , 0 ) alongside E ( Ψ , 0 ) which represent two equilibrium points of system (98). As expounded by the theory of planar dynamical system, we realize that E is a center and also E + is a saddle point. We indicate here that h ± = H ( Ψ ± , 0 ) , and, by doing a careful computation, we achieve:
h ± = 1 12 B 0 2 B 1 2 4 B 0 B 2 B 1 ± B 1 2 4 B 0 B 2 + 2 B 0 B 1 B 2 .
Evidently, h < h < h + and we have it that H ( Ψ , w ) = h + correlates to homoclinic orbits. Moreover, H ( Ψ , w ) = h relates to the center E and then H ( Ψ , w ) = h , where h + < h < h is related to a class of closed orbits that surround center E which are encompassed by a homoclinic orbit. Meaning that (109) defines bounded solutions if and only if the condition given as h + h < h holds. Precisely, (109) explains a family of periodic solutions whenever h + < h < h .
When h = h + , Equation (109) explains a bounded solution that tends towards Ψ + as z goes to infinity. In fact,
2 B 0 3 Ψ 3 + B 1 Ψ 2 + 2 B 2 Ψ + 2 h + = 2 B 0 3 ( Ψ Ψ + ) 2 ( Ψ Ψ 0 ) ,
with Ψ 0 = ( B 1 + 2 B 1 2 4 B 0 B 2 ) / 2 B 0 . In consequence (109) can be reduced to:
Ψ 0 Ψ d Ψ ( Ψ Ψ + ) B 0 ( Ψ Ψ 0 ) = 2 3 ( z z 0 ) ,
from which we can get the exact solution in the structure of a secant hyperbolic
Ψ = Ψ + ( Ψ + Ψ 0 ) sech 2 B 0 ( Ψ + Ψ 0 ) 6 ( z z 0 ) ,
where z = p x + q y + r t and z 0 is an arbitrary constant. By further simplification, Equation (111) becomes:
Ψ = Ψ + 3 B 1 2 4 B 0 B 2 2 B 0 sech 2 1 2 B 1 2 4 B 0 B 2 1 / 4 ( z z 0 ) ,
and this is regarded as an exact bounded solution of (97).
Therefore, we consider the lemma stated as follows.
Lemma 1.
The general second-order ODE (97) has bounded solutions if and only if B 1 2 4 B 0 B 2 > 0 . The bounded solutions can be expressed as (109) in an implicit form. In fact, provided h < h < h + , (109) defines a family of bounded periodic solutions and h = h + defines a bounded solution which approaches Ψ + as z goes to infinity and can be expressed explicitly as (112), where Ψ + = ( B 1 + B 1 2 4 B 0 B 2 ) / 2 B 0 and h ± is defined by (110).

Bounded Travelling Wave Solutions to the Generalized (2+1)-Dimensional Bogoyavlensky–Konopelchenko Equation

According to analysis and results in the above subsection, it is evident that (97) possesses only two kinds of bounded solutions, amidst of which one is found out to be a family of periodic solutions whereas another is discovered to be a family of solutions which approaches a fixed number as z tends to infinity. It is noteworthy to assert here that what we are targeting is to study the bounded travelling wave solutions associated with (2+1)-D genBKe (4) which are determined via Ψ = d ψ / d z , and Ψ satisfies (97). So we have to investigate how we can get the bounded solution of (97).
Visibly, ψ ( z ) = z 0 z Ψ ( z ) d z whereas Ψ ( z ) can implicitly be expressed as stated in (109). By virtue of the geometry meaning of the integral as well as the properties of the solutions of (97), we get the travelling wave solutions to the (2+1)-D genBKe (4). In order to achieve the bounded solutions needed, we choose the integral constant C 1 to be zero that implies B 2 = 0 in (112) and as such
ψ ( z ) = C 1 3 | B 1 | B 0 tanh | B 1 | 2 ( z z 0 ) ;
which means
ψ ( z ) = C 1 + 16 p 3 ( r + σ q + γ p ) p 2 ( α p + β q ) tanh 1 2 ( r + σ q + γ p ) p 2 ( α p + β q ) ( z z 0 ) ,
that is, the family of analytic bounded kink traveling wave solutions to the (2+1)-D genBKe (4), with z = p x + q y + r t and z 0 alongside C 1 regarded as arbitrary constants.
Nonetheless, we may not be able to achieve bounded solutions from the family of periodic solutions of (97). We can easily see that if Ψ ( z ) is a periodic solutions of (97), in the same vein, ψ ( z ) = z 0 z Ψ ( z ) d z is bounded if and only if 0 T Ψ ( z ) d z = 0 , where T represents the period of the function Ψ ( z ) . Recall that the period of the function Ψ ( z ) which is given by (109) with h < h < h + is dependent continuously on the parameters, B 0 , B 1 , B 2 and h. So 0 T Ψ ( z ) d z continuously depends on the parameters, B 0 , B 1 , B 2 and h as well. Suppose we have it that V ( B 0 , B 1 , B 2 , h ) = 0 T Ψ ( z ) d z ; as a consequence, V ( B 0 , B 1 , B 2 , h ) is defined as a continuous function of B 0 , B 1 , B 2 and h. The prove to showcase the existence of the root of V ( B 0 , B 1 , B 2 , h ) = 0 to furnish us with the idea of the existence of the bounded periodic travelling wave solutions to (2+1)-D genBKe (4) is given in [65].
Theorem 4.
The generalized (2+1)-dimensional Bogoyavlensky–Konopelchenko equation possesses two types of bounded travelling wave solutions given as:
(1) 
The generalized (2+1)-dimensional Bogoyavlensky–Konopelchenko equation has a family of analytic bounded kink travelling wave solutions:
u ( t , x , y ) = C 1 + 16 p 3 ( r + σ q + γ p ) p 2 ( α p + β q ) tanh 1 2 ( r + σ q + γ p ) p 2 ( α p + β q ) ( p x + q y + r t z 0 ) ,
where z 0 and C 1 are two arbitrary constants;
(2) 
The generalized (2+1)-dimensional Bogoyavlensky–Konopelchenko equation possesses at least two families of bounded periodic travelling wave solutions which are determined implicitly by (109) and
u ( z ) = z 0 z Ψ ( z ) d z ,
where z = p x + q y + r t and z 0 is an arbitrary constant.

4. Dynamical Wave Behaviour and Analysis of Solutions

The physical phenomena of those secured closed-form solutions can be captured more clearly via graphical evaluation. The obtained solutions of the (2+1)-D genBKe equation comprises kink and anti-kink waves, periodic solitons waves, multi-soliton waves, singular solitons, as well as mixed dark–bright waves of different dynamical structures. Those secure solutions contain several sets of arbitrary constants and functions, which consequently exhibit diverse dynamical structures of multiple solitons through their numerical simulations. We present the structure of the dynamical behaviour of the waves in 3D, 2D and density plots with the aid of Maple software. The singular periodic wave structure in Figure 1 depicts the dynamics of solitary wave solution (34) where we utilize the parameters values γ = 100 , α = 1 , C 0 = 1 , C 1 with variables y = 0 and 1 t , x 1 . Figure 2 represents topological kink soliton solution (36) in 3D, 2D and density plots where we engage values γ = 1 , α = 4 , C 0 = 1 , C 1 = 10 , c 0 = 1 , c 1 = 100 where y = 0 , 10 t 10 and 4 x 4 . Now, for (30), we contemplate a few different choices of arbitrary functions f 1 ( t ) and f 2 ( t ) and for the fact that the solution contains variable y, we consider another function of y as g ( y ) . Therefore, since the solution is a function of t and y, we first consider f ( t ) = 3 sech 4 ( t ) , f ( t ) = ( f 1 ( t ) , f 2 ( t ) ) and g ( y ) = cos ( y ) sin ( y ) , using Maple software, we further illustrate the solution in Figure 3 with the range π t π and 2 π y 3 π where we have x = 0 . Hence, the numerical simulation reveals a doubly-periodic interaction between two-solitons with different amplitudes. Further, we choose f ( t ) = 3 sech 4 ( t ) and g ( y ) = ( 2 tanh ( y ) + cos ( y ) ) in Figure 4 where we have variables x = 0 as well as t and y in the range π t π and 2 π y 3 π . This then exhibits periodic interaction between solitons at varying amplitude and frequency along y t -axis. Moreover, on selecting f ( t ) = 3 sec h ( t ) and g ( y ) = ( 2 tanh 2 ( y ) + sin ( y ) ) , we plot Figure 5 where π t π , 2 π y 3 π and x = 0 . This occasions periodic interaction between solitons travelling at different amplitude but moving in the same direction. In Figure 6 we choose f ( t ) = 3 sec h ( t ) Si ( t ) and g ( y ) = sin ( y ) along with 3 π t 3 π and 2 π y 4 π . We can see in the figure three soliton interactions. These include a kink with t-axis periodic and y-axis periodic, which is clearly revealed in the propagation of the amplitude. Meanwhile, selection of f ( t ) = 3 sec h ( t ) and g ( y ) = ( 2 cn ( t , y ) + sin ( y ) ) with x = 0 , 3 π t π and 2 π y 3 π furnishes doubly-periodic and 1-soliton interactions as portrayed in Figure 7. The interaction depicts an upsurge of wave propagating at varying amplitude, travelling at different velocity and time intervals. Moreover, we can see in Figure 8 a periodic interaction existing between two-solitons with opposite amplitude and propagating at a uniform frequency. This is achieved by allocating f ( t ) = 3 sec h ( t ) and g ( y ) = 3 t cos ( y ) where x = 0 , π t π and 2 π y 6 π . Besides, Figure 9 exhibits wave dynamical behaviour surfacing from a collision between a kink and a soliton solution purveyed by assigning f ( t ) = 4 sec h ( t ) and g ( y ) = t tanh ( y ) with x = 0 , π t π and π y 4 π . Finally on wave interactions, we assign functions f ( t ) = 40 sec h ( t ) and g ( y ) = 20 t sech 2 ( y ) in Figure 10 where x = 0 , π t π and π y 4 π . The resultant effect of the soliton collisions gives a two-soliton wave propagating with opposite amplitude along y t -axis.
Next, the kink solution (72) is depicted with Figure 11 with dissimilar constant values a = 1 , b = 1 , C 1 = 1 , C 2 = 2 , α = 2 , δ = 1 , γ = 1 , ρ = 1 at y = 1 and 10 t , x 10 . The various dynamical behaviour of periodic solution (84) is exhibited in Figure 12, Figure 13 and Figure 14 using parameter values a = 1 , b = 1 , c 0 = 1 , c 1 = 1 , c 2 = 1 , α = 1 , β 3 = 1 , β 4 = 1 , δ = 1 , ρ = 1 , θ 1 = 9 , θ 2 = 1 , θ 3 = 1 , Ω 0 2 = 0.09 at t = 2 and 2 x , y 2 , a = 1 , b = 1 , c 0 = 1 , c 1 = 1 , c 2 = 1 , α = 1 , β 3 = 1 , β 4 = 1 , δ = 1 , ρ = 1 , θ 1 = 9 , θ 2 = 1 , θ 3 = 1 , Ω 0 2 = 0.09 at t = 5 and 2 x , y 2 as well as a = 1 , b = 1 , c 0 = 1 , c 1 = 1 , c 2 = 1 , α = 1 , β 3 = 1 , β 4 = 1 , δ = 1 , ρ = 1 , θ 1 = 9 , θ 2 = 1 , θ 3 = 1 , Ω 0 2 = 0.09 at t = 2 and 2 x , y 2 accordingly. Moreover, the motion character of solution are further depicted in Figure 15 and Figure 16 respectively via values a = 1 , b = 1 , c 0 = 1 , c 1 = 1 , c 2 = 1 , α = 1 , β 3 = 1 , β 4 = 1 , δ = 1 , ρ = 1 , θ 1 = 40 , θ 2 = 2 , θ 3 = 5 , Ω 0 2 = 0.26 at t = 2 and 2 x , y 2 alongside a = 1 , b = 1 , c 0 = 1 , c 1 = 1 , c 2 = 1 , α = 1 , β 3 = 1 , β 4 = 1 , δ = 1 , ρ = 1 , θ 1 = 50 , θ 2 = 5 , θ 3 = 5 , Ω 0 2 = 0.26 at t = 3 and 2 x , y 2 . The Weierstrass elliptic function solution (60) is represented graphically in Figure 17 with unalike parametric values a = 1 , b = 1 , c 0 = 1 , c 1 = 2 , α = 2 , α 0 = 1 , α 1 = 1 , α 2 = 2 , δ = 1 , ρ = 1 , A 0 = 1 , A 1 = 2 where y = 1 and 10 x , y 10 . This wave depiction reveals a multi-soliton wave structure which is a significant wave in nonlinear science and engineering.
Further, we depict the elliptic integral solution (68) in Figure 18, Figure 19, Figure 20 and Figure 21. This is achieved by invoking dissimilar constant values a = 1 , b = 1 , c 0 = 1 , c 2 = 1 , α = 2 , α 4 = 1 , α 5 = 1 , δ = 1 , ρ = 1 , ϑ 1 = 3 , ϑ 2 = 2 , ϑ 3 = 1 , Δ 2 = 0.09 at t = 2 and 1 x , y 1 , a = 1 , b = 1 , c 0 = 1 , c 2 = 1 , α = 2 , α 4 = 1 , α 5 = 5 , δ = 1 , ρ = 1 , ϑ 1 = 3 , ϑ 2 = 2 , ϑ 3 = 1 , Δ 2 = 0.09 when t = 1 and 1 x , y 1 , a = 1 , b = 1 , c 0 = 1 , c 2 = 1 , α = 2 , α 4 = 1 , α 5 = 5 , δ = 1 , ρ = 1 , ϑ 1 = 3 , ϑ 2 = 2 , ϑ 3 = 1 , Δ 2 = 0.09 at t = 1 and 1 x , y 1 as well as a = 1 , b = 1 , c 0 = 1 , c 2 = 0 , α = 1 , α 4 = 1 , α 5 = 1 , δ = 1 , ρ = 1 , ϑ 1 = 3 , ϑ 2 = 2 , ϑ 3 = 1 , Δ 2 = 0.08 at t = 0 and 1 x , y 1 respectively. We notice that the dynamical wave behaviour of elliptic integral solution (68) reveals a mixed dark and bright soliton wave profile which is akin to hyperbolic secant and hyperbolic tangent functions. It is known that the elliptic solution disintegrates to elementary hyperbolic functions by taking some special limits. These functions comprise secant hyperbolic and tangent hyperbolic. It will be recalled that these two constitute bell and anti-bell shapes respectively. As a consequence, this asserted relationship and the interconnections between elliptic solutions and the involved functions are conspicuously revealed in Figure 18, Figure 19, Figure 20 and Figure 21.
The various nontrivial solitary wave solutions obtained from bifurcation analysis of (2+1)-D genBKe (4) in this study, to actually view their dynamical character, numerical simulation of the involved parameters are performed using Mathematica 11.3. Therefore, we reveal the nontrivial bounded solution (101) via 3D, 2D and density plots in Figure 22 with varying parameter values r = 0.2 , p = 0.1 , q = 0.3 , A = 0.05 , B = 5 , C = 1.02 , C 1 = 8 with t = 0.4 and 6 x , y 6 . The solution (103) is portrayed in Figure 23 using unalike values r = 0.2 , p = 0.1 , q = 0.3 , A = 0.05 , B = 7 , C = 1.05 , C 1 = 9 with t = 0.7 and 8 x , y 8 . Moreover, unbounded solution (106) is represented in Figure 24 through 3D, 2D as well as the density plot with constant values r = 0.2 , p = 0.1 , q = 0.3 , A = 0.5 , B = 5 , C = 1 , C 1 = 4 with t = 0.2 and 10 x , y 10 . We further exhibit the travelling wave solution (113) in Figure 25, Figure 26, Figure 27 and Figure 28 using dissimilar values of parameters respectively given as: r = 0.5 , p = 1 , q = 1 , α = 5 , β = 200 , σ = 90 , γ = 100 , C 1 = 4 with t = 2 and 10 x , y 10 ; r = 0.1 , p = 1 , q = 1 , α = 50 , β = 200 , σ = 90 , γ = 100 , C 1 = 0 with x = 3 and 10 t , y 10 ; r = 0.1 , p = 1 , q = 1 , α = 50 , β = 200 , σ = 90 , γ = 100 , C 1 = 0 with x = 3 and 10 t , y 10 ; r = 0.1 , p = 1 , q = 1 , α = 50 , β = 200 , σ = 90 , γ = 100 , C 1 = 0 with y = 5 and 10 t , x 10 .
Significant observations
Figure 17 portrays a localized wave structure of multi-solitons of Equation (4). The dynamical structure appears due to the balance between nonlinearity and the dispersion term. Figure 18, Figure 19, Figure 20 and Figure 21 depicts the coexistence between bright and dark solitons with various wave structures. It is eminent that bright soliton profiles are identified with hyperbolic secant functions. The bright soliton solution usually assumes a bell-shaped figure and also propagates in an undistorted manner without any variation in shape for arbitrarily long distances. Nevertheless, dark soliton solutions which usually exhibit anti-bell wave structures, configured also as topological optical solitons, are characterized by hyperbolic tangent functions. Moreover, important to note is the fact that Equation (56) which can be seen in various cases of symmetry reductions via optimal subalgebras in this study is reminiscent of the ordinary differential equation (ODE) achieved in the quintessential work conducted by Korteweg along with De Vries in [18]. In addition to that, this ODE is interconnected with long waves which propagate along a rectangular canal. Moreover, ODE (56) delineates stationary waves and by imposing some certain constraints for example having the fluid undisturbed at infinity, Korteweg and De Vries secured negative and positive solitary waves alongside cnoidal wave solutions [18,66].

5. Conservation Laws

This section reveals the constructed conservation laws for (2+1)-D genBKe (4) by the engagement of the multipier method [67] along with the well-known Noether’s theorem [68].

5.1. Conserved Vectors via Homotopy Formula

It is germane to state that the multiplier technique is advantageous in the sense that it works for any PDE either with or without variational principle [6,28,67]. In other words, the multiplier method does not require the availability of variational principle before the conserved vectors of a given PDE is obtained. To derive the conserved vectors of (2+1)-D genBKe (4), we first determine the second-order multipliers via the criteria,
δ δ u ( Λ Δ ) = 0 ,
with Λ = Λ ( t , x , y , u , u t , u x , u y , u x x , u x y ) and the Euler operator δ / δ u expressed as:
δ δ u = u D t u t D x u x D y u y + D t D x u t x + D x D y u x y + D x 2 u x x + D y 2 u y y + D x 4 u x x x x + D x 3 D y u x x x y .
On expanding Equation (114) and using the standard Lie theory algorithm, one achieves:
Λ y y = 0 , Λ y u x = 0 , Λ u x u x = 0 , Λ x = 0 , Λ u = 0 , Λ u t = 0 , Λ u y = 0 , Λ u x x = 0 , Λ u x y = 0 ,
which can be solved without much tedious process thereby giving the value of Λ as
Λ ( t , x , y , u , u t , u x , u y , u x x , u x y ) = f 1 ( t ) y 3 ( ρ δ ) u x f 1 ( t ) + C 1 u x + f 2 ( t ) ,
with arbitrary functions f 1 ( t ) and f 2 ( t ) dependent on t. Meanwhile, the homotopy integral formula [69] for the multiplier can be expressed as:
T = 0 1 u Δ Λ u t u = u ( λ ) D x Δ Λ u t x u = u ( λ ) d λ , X = 0 1 u Δ Λ u x u = u ( λ ) D x Δ Λ u x x u = u ( λ ) + D x 2 Δ Λ u x x x u = u ( λ ) D x 3 Δ Λ u x x x x u = u ( λ ) + u t Δ Λ u t x u = u ( λ ) u y D x Δ Λ u x x y u = u ( λ ) D x 2 Δ Λ u x x x y u = u ( λ ) + u x Δ Λ u x x u = u ( λ ) D x Δ Λ u x x x u = u ( λ ) + u x y Δ Λ u x x y u = u ( λ ) D x Δ Λ u x x x y u = u ( λ ) + u x x Δ Λ u x x x u = u ( λ ) + u x x x Δ Λ u x x x x u = u ( λ ) + u x y D x Δ Λ u x x x y u = u ( λ ) + u x x y Δ Λ u x x x y u = u ( λ ) d λ , Y = 0 1 u Δ Λ u y u = u ( λ ) D y Δ Λ u y y u = u ( λ ) D x Δ Λ u x y u = u ( λ ) D x 3 Δ Λ u x x x y u = u ( λ ) + u y Δ Λ u y y u = u ( λ ) + u x Δ Λ u x y u = u ( λ ) + u x x x Δ Λ u x x x y u = u ( λ ) d λ .
As a consequence, the three multipliers Λ 1 = u x , Λ 2 = f 1 ( t ) y 3 ( ρ δ ) u x f 1 ( t ) and Λ 3 = f 2 ( t ) from (115) secures the conservation laws, accordingly as:
T 1 = 1 4 u x 2 1 4 u u x x , X 1 = 1 8 β u u x x x y + 1 4 σ u u x y + 1 2 ν u u y y + α u x u x x x + 5 8 β u x u x x y + 1 4 σ u x u y 3 8 β u x x u x y + 1 8 β u x x x u y + δ u x 2 u y + 1 2 ρ u x 2 u y + 1 2 γ u x 2 + 1 4 u t u x 1 2 α u x x 2 + 1 4 u u t x + 2 α u x 3 + ρ u u x u x y δ u u x u x y , Y 1 = δ u u x u x x ρ u u x u x x + 1 2 ρ u x 3 + 1 4 σ u x 2 1 8 β u x x 2 1 4 σ u u x x + 1 2 ν u x u y + 1 4 β u x u x x x 1 8 β u u x x x x 1 2 ν u u x y ; T 2 = 3 4 ρ u u x x f 1 ( t ) 3 4 δ f 1 ( t ) u u x x + 3 4 δ u x 2 f 1 ( t ) 3 4 ρ u x 2 f 1 ( t ) + 1 2 y u x f 1 ( t ) , X 2 = 3 4 δ u t u x f 1 ( t ) 3 4 ρ u t u x f 1 ( t ) + 3 α u x 2 y f 1 ( t ) + 3 2 γ δ u x 2 f 1 ( t ) + γ y u x f 1 ( t ) + α y u x x x f 1 ( t ) + 3 4 β y u x x y f 1 ( t ) + 1 2 σ y u y f 1 ( t ) + 6 α δ u x 3 f 1 ( t ) 6 α ρ u x 3 f 1 ( t ) + 3 δ 2 u x 2 u y f 1 ( t ) 3 2 ρ 2 u x 2 u y f 1 ( t ) + 3 4 δ u u t x f 1 ( t ) 3 4 ρ u u t x f 1 ( t ) + 3 2 α ρ u x x 2 f 1 ( t ) 3 4 δ u u x f 1 ( t ) 3 2 γ ρ u x 2 f 1 ( t ) 3 2 α δ u x x 2 f 1 ( t ) 3 4 ρ σ u x u y f 1 ( t ) 3 ρ 2 u u x u x y f 1 ( t ) 3 2 δ ρ u x 2 u y f 1 ( t ) 15 8 β ρ f 1 ( t ) u x u x x y + 3 4 δ σ u x u y f 1 ( t ) 3 δ 2 u u x u x y f 1 ( t ) 3 4 ρ σ u u x y f 1 ( t ) + 3 4 δ σ u u x y f 1 ( t ) 3 8 β ρ u u x x x y f 1 ( t ) + 3 8 β δ u u x x x y f 1 ( t ) + 3 4 ρ y u x u y f 1 ( t ) + 6 δ ρ u u x u x y f 1 ( t ) 1 4 β u x x f 1 ( t ) 1 2 u y f 1 ( t ) 1 2 σ u f 1 ( t ) + 1 2 y u t f 1 ( t ) 3 2 ν ρ u u y y f 1 ( t ) + 3 α δ u x u x x x f 1 ( t ) 3 α ρ u x u x x x f 1 ( t ) + 15 8 β δ u x u x x y f 1 ( t ) + 3 2 δ ν f 1 ( t ) u u y y 3 2 δ y u u x y f 1 ( t ) + 3 4 ρ y u u x y f 1 ( t ) + 3 2 δ y u x u y f 1 ( t ) 9 8 β δ u x x u x y f 1 ( t ) + 9 8 β ρ u x x u x y f 1 ( t ) + 3 8 β δ u x x x u y f 1 ( t ) 3 8 β ρ u x x x u y f 1 ( t ) , Y 2 = 3 2 ν ρ u u x y f 1 ( t ) 3 2 δ ν u u x y f 1 ( t ) + 3 8 β ρ u u x x x x f 1 ( t ) + 3 4 ρ σ u u x x f 1 ( t ) 3 8 β δ u u x x x x f 1 ( t ) 3 4 δ σ u u x x f 1 ( t ) + 3 2 δ y u u x x f 1 ( t ) 3 4 ρ y u u x x f 1 ( t ) + 3 ρ 2 u u x u x x f 1 ( t ) 3 2 ν ρ u x u y f 1 ( t ) + 3 δ 2 u u x u x x f 1 ( t ) 6 δ ρ u u x u x x f 1 ( t ) + 3 4 β δ u x u x x x f 1 ( t ) 3 4 β ρ u x u x x x f 1 ( t ) + 3 2 δ ν u x u y f 1 ( t ) + 3 2 δ ρ u x 3 f 1 ( t ) + 3 4 ρ y u x 2 f 1 ( t ) + 3 4 δ σ u x 2 f 1 ( t ) 3 4 ρ σ u x 2 f 1 ( t ) 3 8 β δ u x x 2 f 1 ( t ) + 3 8 β ρ u x x 2 f 1 ( t ) + 1 2 σ y u x f 1 ( t ) + 1 4 β y u x x x f 1 ( t ) + ν y u y f 1 ( t ) 3 2 ρ 2 u x 3 f 1 ( t ) ν u f 1 ( t ) ; T 3 = 1 2 u x f 2 ( t ) , X 3 = 3 4 β u x x y f 2 ( t ) + 3 α u x 2 f 2 ( t ) + γ u x f 2 ( t ) + 1 2 σ u y f 2 ( t ) + α u x x x f 2 ( t ) 3 2 δ u u x y f 2 ( t ) + 3 4 ρ u u x y f 2 ( t ) + 3 2 δ u x u y f 2 ( t ) + 3 4 ρ u x u y f 2 ( t ) + 1 2 u t f 2 ( t ) 1 2 u f 2 ( t ) , Y 3 = 1 2 σ u x f 2 ( t ) + 1 4 β u x x x f 2 ( t ) + ν u y f 2 ( t ) + 3 2 δ u u x x f 2 ( t ) 3 4 ρ u u x x f 2 ( t ) + 3 4 ρ u x 2 f 2 ( t ) .

5.2. Conserved Vectors via Noether Theorem

This subsection furnishes the Noether theorem [68,69] to achieve the conserved currents of the (2+1)-D genBKe (4) with ρ = 2 δ . Consequently, Equation (4) admits a Lagrangian Lagrangian ( L ) whose equivalent minimal differential order is given as:
L = 1 2 β u x x u x y 1 2 u t u x α u x 3 + 1 2 α u x x 2 1 2 γ u x 2 3 2 δ u x 2 u y 1 2 σ u x u y 1 2 ν u y 2 ,
which can easily be ascertained by inspection. Thus we arrive at a Lemma:
Lemma 2.
The (2+1)-D genBKe (4) forms the Euler–Lagrange equation with the functional
J ( v ) = 0 0 0 L ( t , x , y , u t , u x , u y , u x x , u x y ) d t d x d y ,
where the conforming function of Lagrange L is as given in (117).
We achieve variational symmetry P by employing symmetry invariance condition expressed as:
p r ( 2 ) P L + L [ D t ( ξ 1 ) + D x ( ξ 2 ) + D y ( ξ 3 ) ] = D t ( B t ) + D x ( B x ) + D y ( B y ) ,
with the gauge functions B t , B x and B y depending on ( t , x , y , u ) . In addition, the second prolongation p r ( 2 ) P of P can be recovered by the relation:
p r ( 2 ) P = P + ζ t u t + ζ x u x + ζ y u y + ζ x x u x x + ζ x y u x y ,
with the variable coefficients as defined in (7) and P = ξ 1 / x + ξ 2 / y + ξ 3 / t + η / u . Separating the monomials from the expansion of (118) secures the presented system of linear partial differential equations. They are:
ξ x 1 = 0 , B u t + 2 ξ u 1 = 0 , ξ u 1 + B u t = 0 , ξ t 1 + ξ y 3 B t t + 2 η u 3 ξ x 2 = 0 , ξ u 1 = 0 , ξ x 2 = 0 , η x = 0 , ξ u 2 = 0 , ξ u 3 = 0 , ξ u u 1 = 0 , ξ u u 2 = 0 , ξ u 1 + B u t = 0 , η u u 2 ξ x u 2 = 0 , 2 η x u ξ x x 2 = 0 , ξ u 1 = 0 , ξ x 1 = 0 , ξ u 3 = 0 , ξ u 3 = 0 , ξ x 3 = 0 , ξ x u 1 = 0 , B u t + ξ u 1 = 0 , ξ x u 1 = 0 , ξ x x 1 = 0 , ξ u u 3 = 0 , ξ x u 3 = 0 , ξ x u 3 = 0 , ξ x x 3 = 0 , η x x = 0 , B u t + 2 ξ u 1 = 0 , ξ u 3 = 0 , ξ u 1 + B u t = 0 , 4 α ξ u 3 + 5 β ξ u 2 = 0 , 2 α ξ u 3 + 3 δ ξ u 2 = 0 , 4 α ξ x u 1 + β ξ y u 1 = 0 , 2 α ξ u u 3 + β ξ u u 2 = 0 , 2 α η x x + β η x y = 0 , β ξ x y 1 + 2 α ξ x x 1 = 0 , B x x + B y y = 0 , 2 α η u u β ξ u y 2 4 α ξ x u 2 = 0 , 2 β η x u 2 β ξ x y 3 2 α ξ x x 3 = 0 , 6 δ ξ x 3 + σ ξ u 3 + ν ξ u 2 = 0 , σ ξ x 1 + 2 ν ξ y 1 + ξ x 3 = 0 , σ η x + 2 ν η y + 2 B u y = 0 , β η u u 4 α ξ x u 3 β ξ u y 3 β ξ x u 2 = 0 , β η u y β ξ x y 2 + 4 α η x u 2 α ξ x x 2 = 0 , σ ξ u 1 + 6 δ ξ x 1 + σ B u t + ξ u 3 = 0 , η t + 2 γ η x + σ η y + 2 B u x = 0 , 6 α ξ x 1 + 3 δ ξ y 1 + γ ξ u 1 + γ B u t + ξ u 2 = 0 , B t t ξ y 3 + 2 γ ξ x 1 + σ ξ y 1 2 η u = 0 , ξ t 2 γ ξ t 1 γ ξ y 3 + γ B t t 2 γ η u + γ ξ x 2 + σ ξ y 2 6 α η x 3 δ η y = 0 γ ξ u 2 2 α ξ t 1 2 α ξ y 3 + 2 α B t t 6 α η u + 4 α ξ x 2 + 3 δ ξ y 2 = 0 , 6 α ξ x 3 9 δ η u + 3 δ ξ x 2 + γ ξ u 3 + σ ξ u 2 + 3 δ B t t 3 δ ξ t 1 = 0 , 2 γ ξ x 3 6 δ η x 2 σ η u + 2 ν ξ y 2 σ ξ t 1 + σ B t t + ξ t 3 = 0 , 2 β η u 3 β ξ x 2 4 α ξ x 3 + β ξ t 1 + β ξ y 3 β B t t = 0 , σ ξ x 3 ν ξ x 2 + ν B t t ν ξ t 1 2 ν η u + ν ξ y 3 = 0 .
We achieve the solution of the system with regards to ξ 1 , ξ 2 , ξ 3 and η as
ξ 1 = c 1 t + c 2 , ξ 3 = 2 3 c 1 y 4 α ν 9 δ c 1 t + 1 3 c 1 σ t + c 3 , B t = 2 3 c 1 t + F 3 ( x , y ) , ξ 2 = 1 3 c 1 x + 2 α 9 δ c 1 y + F 1 ( t ) , η = 1 27 δ 2 ( 2 σ α c 1 + 6 γ c 1 δ 9 δ F 1 ( t ) ) y + F 2 ( t ) , B x = 1 54 δ 2 ( 6 γ c 1 δ σ 2 c 1 α σ 2 9 δ σ F 1 ( t ) 9 δ y F 1 ( t ) 27 δ 2 F 2 ( t ) ) u + G ( t , x , y ) , B y = 1 27 δ 2 ν ( 6 γ c 1 δ 2 σ α c 1 9 δ F 1 ( t ) ) u G x ( t , x , y ) d y + F 4 ( t , x ) .
Functions F 1 ( t ) , F 2 ( t ) , F 3 ( x , y ) , F 4 ( x , t ) , and G ( t , x , y ) in the solution are arbitrary so are constants c 1 , c 2 and c 3 . Thus, we have the five Noether symmetries together with their respective gauge functions as:
P 1 = t , B t = 0 , B x = 0 , B y = 0 , P 2 = y , B t = 0 , B x = 0 , B y = 0 , P 3 = t t + 1 3 x + 2 α 9 δ y x + 2 3 y 4 α ν 9 δ t + 1 3 σ t y 1 27 δ 2 6 γ δ 2 σ α u , B t = 2 3 t , B x = σ 54 δ 2 6 γ δ 2 α σ u , B y = ν 27 δ 2 6 γ δ 2 α σ u , P F 1 = F 1 ( t ) x + 1 3 δ y F 1 ( t ) u , B t = 0 , B x = 9 σ 54 δ F 1 ( t ) + 9 54 δ y F 1 ( t ) u , B y = ν 3 δ F 1 ( t ) u , P F 2 = F 2 ( t ) u , B t = 0 , B x = 1 2 F 2 ( t ) u , B y = 0 .
We invoke the relation [70]:
T k = L τ k + ( ξ α ψ x j α τ j ) L ψ x k α l = 1 k D x l L ψ x l x k α + l = k n ( η l α ψ x l x j α τ j ) L ψ x k x l α ,
to secure the conserved vectors for the six Noether symmetries respectively as:
T 1 t = 1 2 α u x x 2 α u x 3 1 2 γ u x 2 + 1 2 β u x x u x y 3 2 δ u x 2 u y 1 2 σ u x u y 1 2 ν u y 2 , T 1 x = 3 α u t u x 2 + α u t u x x x α u x x u t x + γ u t u x + 3 4 β u t u x x y 1 4 β u x x u t y 1 2 β u t x u x y + 3 δ u t u x u y + 1 2 σ u t u y + 1 2 u t 2 , T 1 y = 1 4 β u t u x x x 1 4 β u x x u t x + 3 2 δ u t u x 2 + 1 2 σ u t u x + ν u t u y ; T 2 t = 1 2 u x u y , T 2 x = 1 2 u t u y + 3 α u x 2 u y + α u x x x u y α u x x u x y + 3 4 β u y u x x y 1 2 β u x y 2 1 4 β u x x u y y + γ u x u y + 3 δ u x u y 2 + 1 2 σ u y 2 , T 2 y = 1 2 α u x x 2 1 2 u t u x α u x 3 1 2 γ u x 2 + 1 4 β u x x u x y + 1 4 β u x x x u y + 1 2 ν u y 2 ; T 3 t = 1 6 x u x 2 t α u x 3 1 2 γ t u x 2 3 2 δ t u y u x 2 + 1 9 δ α y u x 2 1 27 δ 2 α σ y u x + 1 3 y u y u x 2 9 δ α ν t u y u x 1 3 t σ u y u x + 1 9 δ y γ u x 1 2 ν t u y 2 + 1 2 α t u x x 2 + 1 2 β t u x y u x x , T 3 x = 2 3 x α u x 3 + 4 9 δ α 2 y u x 3 + 1 6 γ x u x 2 2 9 δ 2 α 2 σ y u x 2 + 7 3 α y u y u x 2 + 1 2 δ x u y u x 2 4 3 δ α 2 ν t u y u x 2 + α σ t u y u x 2 + 3 α t u t u x 2 + 7 9 δ α γ y u x 2 + 2 δ y u y 2 u x 4 3 α ν t u y 2 u x + δ σ t u y 2 u x 2 27 δ 2 α γ σ y u x + 4 3 γ y u y u x 4 9 δ α γ ν t u y u x + 1 3 γ σ t u y u x 2 9 δ α σ y u y u x 1 6 β u x y u x 1 3 α u x x u x 1 18 δ α β u x x u x + 1 4 β x u x x y u x + 1 6 δ α β y u x x y u x + 1 3 x α u x x x u x + 2 9 δ α 2 y u x x x u x + γ t u t u x + 3 δ t u y u t u x + 2 9 δ γ 2 y u x + 1 6 σ 2 t u y 2 1 6 ν x u y 2 1 9 δ α ν y u y 2 + 1 3 σ y u y 2 2 9 δ α ν σ t u y 2 1 3 β y u x y 2 + 2 9 δ α β ν t u x y 2 1 6 β σ t u x y 2 1 6 α x u x x 2 1 9 δ α 2 y u x x 2 1 27 δ 2 α σ 2 y u y + 1 9 δ γ σ y u y + 1 54 δ 2 α β σ u x x 1 6 β u y u x x 1 6 y β u y y u x x + 1 9 δ α β ν t u y y u x x 1 12 β σ t u y y u x x + 1 2 t u t 2 2 3 α y u x y u x x 1 12 β x u x y u x x + 4 9 δ α 2 ν t u x y u x x 1 3 α σ t u x y u x x 1 18 δ α β y u x y u x x 1 18 δ β γ u x x 1 18 δ 2 α β σ y u x x y + 1 2 β y u y u x x y 1 3 δ α β ν t u y u x x y + 1 4 β σ t u y u x x y + 1 6 δ β γ y u x x y 2 27 δ 2 α 2 σ y u x x x + 2 3 α y u y u x x x 4 9 δ α 2 ν t u y u x x x + 1 3 α σ t u y u x x x + 2 9 δ α γ y u x x x 1 27 δ 2 α σ y u t + 1 3 y u y u t 2 9 δ α ν t u y u t + 2 3 σ t u y u t + 3 4 β t u x x y u t + α t u x x x u t + 1 9 δ γ y u t 1 4 β t u x x u t y 1 2 β t u x y u t x α t u x x u t x 1 9 δ γ σ u + 1 27 δ 2 α σ 2 u , T 3 y = 1 2 δ x u x 3 1 3 α y u x 3 + 4 9 δ α 2 ν t u x 3 1 3 α σ t u x 3 + 2 9 δ α γ ν t u x 2 + 1 6 σ x u x 2 1 6 γ σ t u x 2 + 3 2 δ t u t u x 2 1 27 δ 2 α σ 2 y u x + 1 9 δ γ σ y u x + 1 3 ν x u y u x + 2 9 δ α ν y u y u x 1 12 β u x x u x + 1 12 β x u x x x u x + 1 18 δ α β y u x x x u x 1 3 y u t u x + 2 9 δ α ν t u t u x + 1 3 σ t u t u x 2 9 δ α ν 2 t u y 2 + 1 3 ν y u y 2 + 1 6 ν σ t u y 2 + 1 3 α y u x x 2 1 12 β x u x x 2 2 9 δ α 2 ν t u x x 2 + 1 6 α σ t u x x 2 1 18 δ α β y u x x 2 + 2 9 δ γ ν y u y 2 27 δ 2 α ν σ y u y + 1 6 β y u x y u x x 1 9 δ α β ν t u x y u x x + 1 12 β σ t u x y u x x 1 54 δ 2 α β σ y u x x x + 1 6 β y u y u x x x 1 9 δ α β ν t u y u x x x + 1 12 β σ t u y u x x x + 1 18 δ β γ y u x x x + ν t u y u t + 1 4 β t u x x x u t 1 4 β t u x x u t x 2 9 δ γ ν u + 2 27 δ 2 α ν σ u ; T F 1 t = 1 2 u x 2 F 1 ( t ) 1 6 δ y u x F 1 ( t ) , T F 1 x = 2 α u x 3 F 1 ( t ) + α u x x x u x F 1 ( t ) 1 2 α u x x 2 F 1 ( t ) + 1 2 γ u x 2 F 1 ( t ) + 3 4 β u x u x x y F 1 ( t ) 1 4 β u x x u x y F 1 ( t ) + 3 2 δ u x 2 u y F 1 ( t ) 1 6 δ y u t F 1 ( t ) 1 2 ν u y 2 F 1 ( t ) + 1 12 δ β u x x F 1 ( t ) 1 δ α y u x 2 F 1 ( t ) 1 3 δ α y u x x x F 1 ( t ) 1 4 δ β y u x x y F 1 ( t ) 1 3 δ γ y u x F 1 ( t ) y u x u y F 1 ( t ) 1 6 δ σ y u y F 1 ( t ) + 1 6 δ σ u F 1 ( t ) + 1 6 δ y u F 1 ( t ) , T F 1 y = 1 4 β u x x x u x F 1 ( t ) 1 4 β u x x 2 F 1 ( t ) + 3 2 δ u x 3 F 1 ( t ) + 1 2 σ u x 2 F 1 ( t ) + ν u x u y F 1 ( t ) 1 12 δ β y u x x x F 1 ( t ) 1 6 δ σ y u x F 1 ( t ) 1 2 y u x 2 F 1 ( t ) 1 3 δ ν y u y F 1 ( t ) + 1 3 δ ν u F 1 ( t ) ; T F 2 t = 1 2 u x F 2 ( t ) , T F 2 x = 3 α u x 2 F 2 ( t ) α u x x x F 2 ( t ) γ u x F 2 ( t ) 3 4 β u x x y F 2 ( t ) 3 δ u x u y F 2 ( t ) 1 2 σ u y F 2 ( t ) 1 2 u t F 2 ( t ) + 1 2 u F 2 ( t ) , T F 2 y = 1 4 β u x x x F 2 ( t ) 3 2 δ u x 2 F 2 ( t ) 1 2 σ u x F 2 ( t ) ν u y F 2 ( t ) .

6. Particular Notes on the Conservation Laws

In the latter part of our investigation in this study, local conservation laws, which have an important place in the use of linearization techniques, numerical schemes as well as stability analysis of solutions were achieved. It is well understood that conservation laws are the key ingredients in a bid to deduce the physical aspects of the underlying model. Some well known conserved quantities in physics are the conservation of mass (or matter), energy (power), momentum (linear or angular) as well as Hamiltonian. For instance, the conservation of energy is a consequence of the time invariance of physical systems. In this regard, added to the fact already known that the prevalence of functions in the conserved quantities reveals that the model under consideration has a limitless number of conservation laws, T 1 , X 1 and Y 1 correspond to conservation of momentum.

7. Conclusions

This paper presents a study carried out on the (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko Equation (4). Lie group analysis is invoked to obtain solutions to the equation via the corresponding optimal system of Lie subalgebras in one dimension where various members of the system are engaged to perform the reductions of (4). As a result of the action, diverse solitary wave solutions were achieved and these include elliptic integrals, trigonometric, Weierstrass, complex, topological kink and anti-kink functions. Moreover, on adopting the bifurcation theory of dynamical systems, we obtained nontrivial bounded and unbounded travelling wave solutions of (4) comprising algebraic, rational, periodic, hyperbolic as well as trigonometric functions. Numerical simulations of the various results gained are performed, analyzed and discussed. Further to that, we derived conservation laws of the equation by engaging the multiplier technique and Noether’s theorem where we secured various local conserved vectors. In addition to the diverse advantages and merits of the achieved solutions in this study in various fields of science and engineering, the conservation laws investigated are also of importance. In classical physics, we have these laws consisting of the conservation of energy, and linear as well as angular momentum. Conserved quantities are crucial to our comprehension of the physical world which are seen to be basic laws of nature. Thus, they possess a wide range of applications in physics, and in other diverse fields of study, for instance, chemistry and engineering to mention a few. Some of these applications have been given earlier. Therefore, our results can be utilized for experimental and applied purposes for further studies in various areas of research in science, technology and engineering.

Author Contributions

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

Funding

This research received no external funding.

Acknowledgments

O. D. Adeyemo thanks Shandong University of Science and Technology, People’s Republic of China for the three months financial support given during his visit. It shall not be forgotten.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Abbreviations

Abbreviations

Abbreviations

ODEsOrdinary Differential equations
LODEsLinear Ordinary Differential equations
NODEsNonlinear Ordinary Differential equations
PDEsPartial differential equations
NLDEsNonlinear differential equations
NLPDEsNonlinear partial differential equations
LIPDEsLinear partial differential equations
KdVKortweg-de Vries
KPKadomtsev–Petviashvili
KP-MEWKadomtsev–Petviashvili-Modified Equal Width equation
KP-BBMKadomtsov-Petviashivilli-Benjamin-Bona-Mahony
(2+1)-D genBKe(2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation
2DTwo-dimensional
3DThree-dimensional

References

  1. Adeyemo, O.D.; Motsepa, T.; Khalique, C.M. A study of the generalized nonlinear advection-diffusion equation arising in engineering sciences. Alex. Eng. J. 2022, 61, 185–194. [Google Scholar] [CrossRef]
  2. 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]
  3. Shafiq, A.; Khalique, C.M. Lie group analysis of upper convected Maxwell fluid flow along stretching surface. Alex. Eng. J. 2020, 59, 2533–2541. [Google Scholar] [CrossRef]
  4. Gao, X.Y.; Guo, Y.J.; Shan, W.R. Cosmic dusty plasmas via a (3+1)-dimensional generalized variable-coefficient Kadomtsev-Petviashvili-Burgers-type equation: Auto-Bäcklund transformations, solitons and similarity reductions plus observational/experimental supports. Waves Random Complex Media 2021, 1–21. [Google Scholar] [CrossRef]
  5. Khalique, C.M.; Adeyemo, O.D. Langrangian formulation and solitary wave solutions of a generalized Zakharov-Kuznetsov equation with dual power-law nonlinearity in physical sciences and engineering. J. Ocean Eng. Sci. 2021. [Google Scholar] [CrossRef]
  6. Khalique, C.M.; Adeyemo, O.D. Closed-form solutions and conserved vectors of a generalized (3+1)-dimensional breaking soliton equation of engineering and nonlinear science. Mathematics 2020, 8, 1692. [Google Scholar] [CrossRef]
  7. Adeyemo, O.D. Applications of cnoidal and snoidal wave solutions via an optimal system of subalgebras for a generalized extended (2+1)-D quantum Zakharov-Kuznetsov equation with power-law nonlinearity in oceanography and ocean engineering. J. Ocean Eng. Sci. 2022, in press. [Google Scholar] [CrossRef]
  8. Adeyemo, O.D.; Khalique, C.M. Dynamics of soliton waves of group-invariant solutions through optimal system of an extended KP-like equation in higher dimensions with applications in medical sciences and mathematical physics. J. Geom. Phys. 2022, 177, 104502. [Google Scholar] [CrossRef]
  9. Adeyemo, O.D.; Khalique, C.M. Dynamical soliton wave structures of one-dimensional Lie subalgebras via group-invariant solutions of a higher-dimensional soliton equation in ocean physics and mechatronics engineering. Commun. Appl. Math. Comput. 2022, in press. [Google Scholar]
  10. Jiong, S. Auxiliary equation method for solving nonlinear partial differential equations. Phys. Lett. A 2003, 309, 387–396. [Google Scholar]
  11. 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]
  12. He, J.H.; Wu, X.H. Exp-function method for nonlinear wave equations. Chaos Soliton Fract. 2006, 30, 70. [Google Scholar] [CrossRef]
  13. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  14. Adomian, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1994. [Google Scholar]
  15. Agrawal, G.P. Nonlinear Fiber Optics; Academic Press: New York, NY, USA, 1995. [Google Scholar]
  16. Shukla, P.K.; Mamun, A.A. Introduction to Dusty Plasma Physics; Institute of Physics Publishing: Bristol, UK, 2002. [Google Scholar]
  17. Kochanov, M.B.; Kudryashov, N.A.; Sinel’shchikov, D.I. Non-linear waves on shallow water under an ice cover, higher order expansions. J. Appl. Math. Mech. 2013, 77, 25–32. [Google Scholar] [CrossRef]
  18. Korteweg, D.J.; Vries, G.D. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Lond. Edinb. Dubl. Phil. Mag. 1895, 39, 422–443. [Google Scholar] [CrossRef]
  19. El-Tantawy, S.A.; Moslem, W.M. Nonlinear structures of the Korteweg-de Vries and modified Korteweg-de Vries equations in non-Maxwellian electron-positron-ion plasma: Solitons collision and rogue waves. Phys. Plasmas 2014, 21, 052112. [Google Scholar] [CrossRef]
  20. Gardner, C.S.; Greene, J.M.; Kruskal, M.D.; Miura, R.M. Method for solving the Korteweg-de Vries equation. Phys. Rev. Lett. 1967, 19, 1095–1097. [Google Scholar] [CrossRef]
  21. Yu, J.; Lou, S.Y. Deformation and (3+1)-dimensional integrable model. Sci. China Ser. A 2000, 43, 655–660. [Google Scholar] [CrossRef]
  22. Lou, S.Y. Searching for higher dimensional integrable models from lower ones via Painlevé analysis. Phys. Lett. 1998, 80, 5027–5031. [Google Scholar] [CrossRef]
  23. El-Wakii, S.A.; Abdou, M.A.; Elhanbaly, A. New solitons and periodic wave solutions for nonlinear evolution equations. Phys. Lett. A 2006, 353, 40–47. [Google Scholar] [CrossRef]
  24. Gu, C.H. Soliton Theory and Its Application; Zhejiang Science and Technology Press: Hangzhou, China, 1990. [Google Scholar]
  25. 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]
  26. Kudryashov, N.A.; Loguinova, N.B. Extended simplest equation method for nonlinear differential equations. Appl. Math. Comput. 2008, 205, 396–402. [Google Scholar] [CrossRef]
  27. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  28. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  29. Olver, P.J. Applications of Lie Groups to Differential Equations, Graduate Texts in Mathematics, 2nd ed.; Springer: Berlin, Germany, 1993; Volume 107. [Google Scholar]
  30. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1994; Volume 1. [Google Scholar]
  31. Ibragimov, N.H. Elementary Lie Group Analysis and Ordinary Differential Equations; John Wiley & Sons: Chichester, NY, USA, 1999. [Google Scholar]
  32. 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]
  33. Matveev, V.B.; Salle, M.A. Darboux Transformations and Solitons; Springer: New York, NY, USA, 1991. [Google Scholar]
  34. Chen, Y.; Yan, Z. New exact solutions of (2+1)–dimensional Gardner equation via the new sine-Gordon equation expansion method. Chaos Solitons Fract. 2005, 26, 399–406. [Google Scholar] [CrossRef]
  35. Kudryashov, N.A. Simplest equation method to look for exact solutions of nonlinear differential equations. Chaos Solitons Fract. 2005, 24, 1217–1231. [Google Scholar] [CrossRef] [Green Version]
  36. Hossain, A.K.M.K.S.; Akbar, M.A. Traveling wave solutions of nonlinear evolution equations via modified simple equation method. Int. J. Appl. Math. Theor. Phys. 2017, 3, 20–25. [Google Scholar] [CrossRef] [Green Version]
  37. Seadawy, A.R. Fractional solitary wave solutions of the nonlinear higher-order extended KdV equation in a stratified shear flow: Part I. Comput. Math. Appl. 2015, 70, 345–352. [Google Scholar] [CrossRef]
  38. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  39. 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]
  40. Jiang, B.; Liu, Y.; Zhang, J. Bifurcations and some new travelling wave solutions for the CH-γ equation. App. Math. Comput. 2014, 228, 220–233. [Google Scholar] [CrossRef]
  41. Saha, A. Bifurcation of travelling wave solutions for the generalized KP-MEW eequations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 3539–3551. [Google Scholar] [CrossRef]
  42. Ganguly, A.; Das, A. Explicit solutions and stability analysis of the (2+1)-dimensional KP-BBM equation with dispersion effect. Commun. Nonlinear Sci. Numer. Simul. 2015, 25, 102–117. [Google Scholar] [CrossRef]
  43. Das, A.; Ganguly, A. Existence and stability of dispersive solutions to the Kadomtsev–Petviashvili equation in the presence of dispersion effect. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 326–339. [Google Scholar] [CrossRef]
  44. Xin, X.; Liu, X.; Zhang, L. Explicit solutions of the Bogoyavlensky-Konoplechenko equation. Appl. Math. Comput. 2010, 215, 3669–3673. [Google Scholar] [CrossRef]
  45. Ray, S.S. On conservation laws by Lie symmetry analysis for (2+1)-dimensional Bogoyavlensky–Konopelchenko equation in wave propagation. Comput. Math. Appl. 2017, 74, 1158–1165. [Google Scholar] [CrossRef]
  46. Triki, H.; Jovanoski, Z.; Biswa, A. Shock wave solutions to the Bogoyavlensky–Konopelchenko equation. Indian J. Phys. 2014, 88, 71–74. [Google Scholar] [CrossRef]
  47. Konopelchenko, B.G. Solitons in Multidimensions-Inverse Spectral Transfrom Method; World Scientific: River Edge, NJ, USA, 1993. [Google Scholar]
  48. Prabhakar, M.V.; Bhate, H. Exact solutions of the Bogoyavlensky-Konoplechenco equation. Lett. Math. Phys. 2003, 64, 1–6. [Google Scholar] [CrossRef]
  49. Bogoyavlenski, O.I. Overturning solitons in new two-dimensional integrable equations. Math. USSR-Izv. 1990, 34, 245–259. [Google Scholar] [CrossRef]
  50. Chen, S.T.; Ma, W.X. Exact solutions to a generalized Bogoyavlensky–Konopelchenko equation via maple symbolic computations. Complexity 2019, 2019, 8787460. [Google Scholar] [CrossRef]
  51. Li, Q.; Chaolu, T.; Wang, Y. Lump-type solutions and lump solutions for the (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation. Comp. Math. Appl. 2019, 77, 2077–2085. [Google Scholar] [CrossRef]
  52. Chen, S.T.; Ma, W.X. Lump solution to a generalized Bogoyavlensky–Konopelchenko equation. Front. Math. China 2018, 13, 525–534. [Google Scholar] [CrossRef]
  53. Pouyanmehr, R.; Hosseini, K.; Ansari, R. Different wave structures to the (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation. Int. J. Appl. Comput. Math. 2019, 5, 1–12. [Google Scholar] [CrossRef]
  54. Liu, F.Y.; Gao, Y.T.; Yu, X.; Li, L.Q.; Ding, C.C. Lie group analysis and analytic solutions for a (2+1)-dimensional generalized Bogoyavlensky–Konopelchenko equation in fluid mechanics and plasma physics. Eur. Phys. J. Plus 2021, 136, 1–14. [Google Scholar] [CrossRef]
  55. Bluman, G.W.; Cole, J.D. Similarity Methods for Differential Equations; Springer: New York, NY, USA, 1974. [Google Scholar]
  56. Olver, P.J. Application of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  57. Grigoriev, Y.N.; Kovalev, V.F.; Meleshko, S.V.; Ibragimov, N.H. Symmetries of Integro-Differential Equations: With Applications in Mechanics and Plasma Physics; Springer: New York, NY, USA, 2010. [Google Scholar]
  58. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic Press: New York, NY, USA, 2007. [Google Scholar]
  59. Akhiezer, N.I. Elements of The Theory of Elliptic Functions; American Mathematical Society: Providence, RI, USA, 1990. [Google Scholar]
  60. Kudryashov, N.A. Analytical Theory of Nonlinear Differential Equations; Igevsk Institute of Computer Investigations: Moskow, Russia, 2004. [Google Scholar]
  61. Kudryashov, N.A. First integrals and general solution of the Fokas-Lenells equation. Optik 2019, 195, 163135. [Google Scholar] [CrossRef]
  62. Billingham, J.; King, A.C. Wave Motion; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  63. Chow, S.N.; Hale, J.K. Method of Bifurcation Theory; Springer: New York, NY, USA, 1981. [Google Scholar]
  64. Guckenheimer, J.; Holmes, P. Dynamical Systems and Bifurcations of Vector Fields; Springer: New York, NY, USA, 1983. [Google Scholar]
  65. Zhang, L.; Khalique, C.M. Traveling wave solutions and infinite-dimensional linear spaces of multiwave solutions to Jimbo-Miwa equation. Abstr. Appl. Anal. 2014, 2014, 963852. [Google Scholar] [CrossRef]
  66. Drazin, P.G.; Johnson, R.S. Solitons: An Introduction; Cambridge University Press: Cambridge, UK, 1989. [Google Scholar]
  67. Anco, S.C.; Bluman, G.W. Direct construction method for conservation laws of partial differential equations. Part I: Examples of conservation law classifications. Eur. J. Appl. Math. 2002, 13, 545–566. [Google Scholar] [CrossRef] [Green Version]
  68. Noether, E. Invariante variationsprobleme. Nachrichten Ges. Wiss. Göttingen Math.-Phys. Kl. 1918, 2, 235–257. [Google Scholar]
  69. Anco, S.C. Generalization of Noethers Theorem in modern Form to Nonvariational Partial Differential Equations. In Recent Progress and Modern Challenges in Applied Mathematics, Modeling and Computational Science; Melnik, R., Makarov, R., Beglair, J., Eds.; Fields Institute Communications; Springer: New York, NY, USA, 2017; Volume 79. [Google Scholar]
  70. 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]
Figure 1. Solitary wave depiction of singular periodic solution (34) at y = 0 .
Figure 1. Solitary wave depiction of singular periodic solution (34) at y = 0 .
Mathematics 10 02391 g001
Figure 2. Solitary wave depiction of topological anti-kink soliton (36) at y = 0 .
Figure 2. Solitary wave depiction of topological anti-kink soliton (36) at y = 0 .
Mathematics 10 02391 g002
Figure 3. Wave depiction of soliton interaction with variant amplitudes at x = 0 .
Figure 3. Wave depiction of soliton interaction with variant amplitudes at x = 0 .
Mathematics 10 02391 g003
Figure 4. Wave profile depiction of soliton interaction with different amplitudes, frequency and also propagating along the same direction when variable x = 0 .
Figure 4. Wave profile depiction of soliton interaction with different amplitudes, frequency and also propagating along the same direction when variable x = 0 .
Mathematics 10 02391 g004
Figure 5. Wave profile depiction of soliton interaction with varying amplitudes but acting and propagating along the same direction where we have variable x = 0 .
Figure 5. Wave profile depiction of soliton interaction with varying amplitudes but acting and propagating along the same direction where we have variable x = 0 .
Mathematics 10 02391 g005
Figure 6. Wave profile depiction of soliton interaction with variant amplitudes and frequency with the wave propagation taking place at different level when x = 0 .
Figure 6. Wave profile depiction of soliton interaction with variant amplitudes and frequency with the wave propagation taking place at different level when x = 0 .
Mathematics 10 02391 g006
Figure 7. Wave profile depiction of soliton interaction with varying amplitudes, frequency and also propagating at different time intervals when variable x = 0 .
Figure 7. Wave profile depiction of soliton interaction with varying amplitudes, frequency and also propagating at different time intervals when variable x = 0 .
Mathematics 10 02391 g007
Figure 8. Wave profile depiction of soliton interaction with variant amplitudes and frequency with the propagation in the opposite directions when we have variable x = 0 .
Figure 8. Wave profile depiction of soliton interaction with variant amplitudes and frequency with the propagation in the opposite directions when we have variable x = 0 .
Mathematics 10 02391 g008
Figure 9. Wave depiction of soliton interaction at different amplitude with x = 0 .
Figure 9. Wave depiction of soliton interaction at different amplitude with x = 0 .
Mathematics 10 02391 g009
Figure 10. Wave profile depiction of soliton interaction at varying amplitude and propagating at a constant velocity and also moving in different directions when x = 0 .
Figure 10. Wave profile depiction of soliton interaction at varying amplitude and propagating at a constant velocity and also moving in different directions when x = 0 .
Mathematics 10 02391 g010
Figure 11. Solitary wave depiction of hyperbolic function solution (72) at y = 1 .
Figure 11. Solitary wave depiction of hyperbolic function solution (72) at y = 1 .
Mathematics 10 02391 g011
Figure 12. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Figure 12. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Mathematics 10 02391 g012
Figure 13. Solitary wave profile depiction of elliptic solution (84) at t = 5 .
Figure 13. Solitary wave profile depiction of elliptic solution (84) at t = 5 .
Mathematics 10 02391 g013
Figure 14. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Figure 14. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Mathematics 10 02391 g014
Figure 15. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Figure 15. Solitary wave profile depiction of elliptic solution (84) at t = 2 .
Mathematics 10 02391 g015
Figure 16. Solitary wave profile depiction of elliptic solution (84) at t = 3 .
Figure 16. Solitary wave profile depiction of elliptic solution (84) at t = 3 .
Mathematics 10 02391 g016
Figure 17. Solitary wave depiction of Weierstrass elliptic solution (60) at y = 1 .
Figure 17. Solitary wave depiction of Weierstrass elliptic solution (60) at y = 1 .
Mathematics 10 02391 g017
Figure 18. Solitary wave depiction of elliptic integral solution (68) at t = 2 .
Figure 18. Solitary wave depiction of elliptic integral solution (68) at t = 2 .
Mathematics 10 02391 g018
Figure 19. Solitary wave depiction of elliptic integral solution (68) at t = 1 .
Figure 19. Solitary wave depiction of elliptic integral solution (68) at t = 1 .
Mathematics 10 02391 g019
Figure 20. Solitary wave depiction of elliptic integral solution (68) at t = 1 .
Figure 20. Solitary wave depiction of elliptic integral solution (68) at t = 1 .
Mathematics 10 02391 g020
Figure 21. Solitary wave depiction of elliptic integral solution (68) at t = 0 .
Figure 21. Solitary wave depiction of elliptic integral solution (68) at t = 0 .
Mathematics 10 02391 g021
Figure 22. Wave profile depiction of nontrivial bounded solution (101) at t = 0.4 .
Figure 22. Wave profile depiction of nontrivial bounded solution (101) at t = 0.4 .
Mathematics 10 02391 g022
Figure 23. Wave profile depiction of nontrivial unbounded solution (103) at t = 0.7 .
Figure 23. Wave profile depiction of nontrivial unbounded solution (103) at t = 0.7 .
Mathematics 10 02391 g023
Figure 24. Wave profile depiction of nontrivial unbounded solution (106) at t = 0.2 .
Figure 24. Wave profile depiction of nontrivial unbounded solution (106) at t = 0.2 .
Mathematics 10 02391 g024
Figure 25. Tavelling wave profile depiction of nontrivial solution (113) at t = 2 .
Figure 25. Tavelling wave profile depiction of nontrivial solution (113) at t = 2 .
Mathematics 10 02391 g025
Figure 26. Tavelling wave profile depiction of nontrivial solution (113) at x = 3 .
Figure 26. Tavelling wave profile depiction of nontrivial solution (113) at x = 3 .
Mathematics 10 02391 g026
Figure 27. Tavelling wave profile depiction of nontrivial solution (113) at x = 3 .
Figure 27. Tavelling wave profile depiction of nontrivial solution (113) at x = 3 .
Mathematics 10 02391 g027
Figure 28. Tavelling wave profile depiction of nontrivial solution (113) at y = 5 .
Figure 28. Tavelling wave profile depiction of nontrivial solution (113) at y = 5 .
Mathematics 10 02391 g028
Table 1. Lie brackets.
Table 1. Lie brackets.
[ X i , X j ] X 1 X 2 X 3 X 4 X 5
X 1 0000 δ X 4
X 2 0000 2 α X 4
X 3 0000 3 δ ρ X 2
X 4 00000
X 5 δ X 4 2 α X 4 3 δ ρ X 2 00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adeyemo, O.D.; Zhang, L.; Khalique, C.M. Bifurcation Theory, Lie Group-Invariant Solutions of Subalgebras and Conservation Laws of a Generalized (2+1)-Dimensional BK Equation Type II in Plasma Physics and Fluid Mechanics. Mathematics 2022, 10, 2391. https://doi.org/10.3390/math10142391

AMA Style

Adeyemo OD, Zhang L, Khalique CM. Bifurcation Theory, Lie Group-Invariant Solutions of Subalgebras and Conservation Laws of a Generalized (2+1)-Dimensional BK Equation Type II in Plasma Physics and Fluid Mechanics. Mathematics. 2022; 10(14):2391. https://doi.org/10.3390/math10142391

Chicago/Turabian Style

Adeyemo, Oke Davies, Lijun Zhang, and Chaudry Masood Khalique. 2022. "Bifurcation Theory, Lie Group-Invariant Solutions of Subalgebras and Conservation Laws of a Generalized (2+1)-Dimensional BK Equation Type II in Plasma Physics and Fluid Mechanics" Mathematics 10, no. 14: 2391. https://doi.org/10.3390/math10142391

APA Style

Adeyemo, O. D., Zhang, L., & Khalique, C. M. (2022). Bifurcation Theory, Lie Group-Invariant Solutions of Subalgebras and Conservation Laws of a Generalized (2+1)-Dimensional BK Equation Type II in Plasma Physics and Fluid Mechanics. Mathematics, 10(14), 2391. https://doi.org/10.3390/math10142391

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