# Dynamical Analysis of Nonlocalized Wave Solutions of2+1-Dimensional CBS and RLW Equation with the Impact of Fractionality and Free Parameters

Authors: Nur Hasan Mahmud Shahen; Foyjonnesa; Md. Rafiqul Islam; Ahmet Bekir; M.M. Rahman
Journal: Advances in Mathematical Physics
 (2022)

---

## Abstract

This study retrieves some new exact solutions to the2+1-dimensional Calogero-Bogoyavlenskii-Schilf (CSB) equation and regularized long wave (RLW) equation in the context of nonlinear traveling wave phenomena. In this regard, the advanced exp−φξ-expansion method is imposed to the 2+1-dimensional CBS and RLW equation, and consequently, rogue, kink, singular kink, periodic, singular, and multiple soliton solutions are exhibited in terms of trigonometric, hyperbolic, and rational function solutions. To enucleate the underlying nonlocalized traveling wave features, accomplished exact solutions are established by making their dynamic behavior of the exact solutions exhibited in three-dimensional (3D) and two-dimensional (2D) combined chart with the help of computational software Maple 18. All of our accomplished solutions are claimed to be new in the sense of conformable derivative, chosen a unique fractional type wave transformation, dynamical behavior of fractionally and free variable, and the imposed method on our preferred equations.


---

## Body

## 1. Introduction

Most of the physical problems are now formatted into mathematical schemes through partial differential equations (PDEs). In recent years, calculus of fractional nonlinear differential condition is one of the rising concern of nonlinear dynamics as well as plays an eventual role in real-life problems such as plasma physics [1], quantum mechanics [2], fluid mechanics [3], optical fibers [4], and other fields of applied mathematics and engineering. Calculus of fractional derivative is the generality of the traditional order of integration and differentiation and a conventional technique for demonstrating complex physical behavior, especially of the nonlinear mathematical science of engineering and physics. In the recent past, nonlinear fractional partial differential equations (FPDEs) have been inventing a potential platform for the researchers to interpret the tangible phenomena. As a result, some significant, concise, and original methods have been explored and explain to exact closed solutions of nonlinear FPDEs, videlicet, namely, the Hirota bilinear method [5, 6], modified extended tanh-function method [7], G′/G-expansion method [8], fractional subequation method [9], modified trial equation method [10], advanced exp−φξ method [11], tan−φξ-expansion method [12], and improved Kudryashov method [13]. Recently, some researchers, like, Ferdous et al. [14], attained exact wave solutions to the extended Zakharov–Kuzetsov equation by implementing the generalized exp−φξ-expansion method. Alike, adopting the exp−φξ-expansion modified technique, Khater [15] has established three different types of nonlocalized traveling wave solutions to the generalized Hirota-Satsuma (HS) couple KdV system. Encouraged by ongoing research on related topics, we extract the recent and more common precise solitary wave solutions to some nonlinear FPDEs, described earlier by proposing the advanced method of exp−φξ-expansion, which can be considered the generalization of the generalized form of exp−φξ-expansion [16].Our current study mainly explores the dynamical changes of conformable time-fractional2+1-dimensional CSB equation [17] and regularized long wave (RLW) equation by employing the new progressive exp−φξ-expansion technique. Recent past, Rahman et al. [18] imposed this technique on some nonfractional NLEEs, but they did not show any productive concept to conformable nonlinear time fractional sense of NPEEs. We decisively believe that our stated advanced exp−φξ modified method will be frolicked a significant role in investigating the traveling wave arrangements of nonlinear science and mathematical physics. Here it is significant to note that our stated fractional equations are in a conformable sense so we have converted our equations in fractional form.Consider the following generalized2+1-dimensional CBS conditions [19–22].
(1)ut+ϕuuy=0,ϕu=∂x2+au+bux∂x−1,or equivalently
(2)ut+uxxy+auuy+bvx∂x−1vy=0,where ∂x−1=∫fdx and a and b are two parameters. If we put dimensional exchange, ∂x=∂y, then Equation (2) changes into a standard (KdV) equation. From there, Equation (2) can be written in the potential time fractional formation of 2+1-dimensional CBS equation [17] with a conformable sense.
(3)uxNtθu+4uxuxy+2uxxuy+uxxy=0,t>0,x,y∈R,where the coefficient of u and its higher derivative term is called the free variable and θ is a time derivative fractional-order parameter with an interval 0<θ≤1.Bogoyavlenskii and Schilf firstly originate this CBS equation in numerous ways [23]. For full possible indispensable conditions, three effective procedures, to be specific the Backlund change strategy, the reverse Scaterin Scheme, and Hirota’s bilinear approach [24–26], were minutely connected to get the precise arrangements of these equations. The CBS equation is resolved in numerous aspects of water wave mechanics, ocean science, nonlinear science, engineering mathematics, etc., under specific authoritative techniques such as the modified simple equation scheme, the Sin-Gordan expansion technique, the unified method, and the G′/G-expansion technique. [27–30]. To the best of our knowledge, our proposed advanced- exp−φξ expansion technique is not utilized yet to find the exact solutions of CBS equation with the concept of conformable derivative. Now we consider the time fractional regularized long wave (RLW) equation [31].
(4)Ntθu+ux+αu2x−uxxNtθu=0,t>0,α,β∈R,where0<θ≤1.The RLW equation is a particular class of nonlinear conformable fractional equations that generates fruitful models for predicting the water wave and physical phenomena [32]. Many authoritative approaches have been the subject of investigating the RLW equation, such as the Fourier leap-frog method [33], finite difference method [34], Hirota direct method approach [35], and modified extended tanh method [36]. Whereas so far we know, with the sense of conformable derivative, there is no prolific results have not yet been established about our declared technique to simplify the RLW equation. The exact wave solutions of RLW equation we attained are completely new in respect of our mentioned method.Recently, Ferdous et al. [14] proposed a generalized exp−φξ-expansion method by choosing the auxiliary ODE of the form φ′ξ−λexpφξ−μexp−φξ=r. In our mentioned method, we take the auxiliary nonlinear ODE of the form φ′ξ−λexpφξ−μexp−φξ=0 by setting r=0. Here λ and μ are real type parameters, and reasonably, our said supplementary form affords much better complete solutions to the FPDEs as well as has a more comprehensive operational physical explanation than the findings of Ferdous et al. [14]. The important favor of our specified technique over the other approaches [15, 18] is that it provides some straight and succinct form of new exact wave solution, as well as it is very useful, well-organized, and responsive applicable in presenting of exact and explicit wave solutions to NPDEs, arises in science, mathematical physics, and engineering.We have employed successfully theexp−φξ method on CBS and RLW equations to treasure the meticulous traveling wave solutions which have not appeared beforehand in publications.We believe that our discoveries might be used to computer simulations of water waves in shallow water, to the study of the propagation of small amplitude long waves on the surface of shallow water, and to the advancement of scientific knowledge. Theoretical physicists, applied mathematicians, and ocean engineers who are interested in the precise solutions to moving rogue waves should all find this paper to be of interest. When this project is successfully completed, it may be possible to forecast the sea level of a coastal area during natural catastrophes. Low pressure systems like cyclones are the main culprits in natural disasters.The paper’s structure is as follows: the conformable derivative’s narration and its elementary possessions are given in Section2. The advanced exp−φξ expansion scheme has been outlined in Section 3. In Section 4, we apply this approach to the time-fractional 2+1-dimensional equation of CBS and RLW. The findings and the debate are discussed in Section 5, and the conclusions are provided in Section 6.


## 2. Inceptions and Procedures

### 2.1. Definition and Some Features of Conformable Fractional Derivative

Khalil et al. [37] primarily exposed the sense of conformable derivative with a limit operator.Definition 1.
Iff:0,∝⟶R, then the conformable derivative with the fractional sense of f order δ is defined as
(5)Ttδft=limε⟶0ft+εt1−δ−ftεforallt>0,0<δ≤1.Later, in the fractional-order method, Abdeljawad [38] also discussed chain law, Gronwall’s inequality, exponential functions, definite and infinite component integration, Taylor power series expansions, and Laplace transform conformable derivative. The difficulty of exiting the modified Riemann-Liouville derivative description [39] will efficiently resolve the concept of a conformable fractional-order derivative.Theorem 2.
Letω∈0,1 and f=ft,g=gt be ω-conformably differentiable at a point t>0; then,
(i)
Ttωcf+dg=cTtωf+dTtωg,forallc,d∈R(ii)
Ttωtγ=γtγ−ω,forallγ∈R(iii)
Ttωfg=gTtωf+fTtωg,(iv)
Ttωf/g=gTtωf−fTtωg/g2
Moreover, if the functionf is differentiable, then Ttωft=t1−ωdf/dt.Theorem 3.
Considerf:0,ω⟶R be a real type function such that f is differentiable and ω-conformable derivable. Also, assume that g be a derivable function well-defined in the range of f. Then, we have Ttωfogt=t1−ωgtω−1g′tTtωftt=gt, where prime indicates the simple derivatives with respect to the point t.In this research, we have careful about the preferred equation with the concept of conformable-based derivative. In circumstance of basic calculus, numerous functions do not expand Taylor’s power-order illustrations on specific points but in the technique of conformable-order derivative, they do consume the existence. CD accomplishes good in the product and chain rule while complicated plans appear in logic of basic fractional geometry. The CD of a constant type function is correspondent to zero wherever it is not the topic for Riemann derivative of fractional order. Mittag Leffler functions show a noteworthy authoritative in fractional-order calculus as interpretation to exponential function where the fractional-order form of exponential type function of the formft=etα/α seems in circumstance of CD.


## 2.1. Definition and Some Features of Conformable Fractional Derivative

Khalil et al. [37] primarily exposed the sense of conformable derivative with a limit operator.Definition 1.
Iff:0,∝⟶R, then the conformable derivative with the fractional sense of f order δ is defined as
(5)Ttδft=limε⟶0ft+εt1−δ−ftεforallt>0,0<δ≤1.Later, in the fractional-order method, Abdeljawad [38] also discussed chain law, Gronwall’s inequality, exponential functions, definite and infinite component integration, Taylor power series expansions, and Laplace transform conformable derivative. The difficulty of exiting the modified Riemann-Liouville derivative description [39] will efficiently resolve the concept of a conformable fractional-order derivative.Theorem 2.
Letω∈0,1 and f=ft,g=gt be ω-conformably differentiable at a point t>0; then,
(i)
Ttωcf+dg=cTtωf+dTtωg,forallc,d∈R(ii)
Ttωtγ=γtγ−ω,forallγ∈R(iii)
Ttωfg=gTtωf+fTtωg,(iv)
Ttωf/g=gTtωf−fTtωg/g2
Moreover, if the functionf is differentiable, then Ttωft=t1−ωdf/dt.Theorem 3.
Considerf:0,ω⟶R be a real type function such that f is differentiable and ω-conformable derivable. Also, assume that g be a derivable function well-defined in the range of f. Then, we have Ttωfogt=t1−ωgtω−1g′tTtωftt=gt, where prime indicates the simple derivatives with respect to the point t.In this research, we have careful about the preferred equation with the concept of conformable-based derivative. In circumstance of basic calculus, numerous functions do not expand Taylor’s power-order illustrations on specific points but in the technique of conformable-order derivative, they do consume the existence. CD accomplishes good in the product and chain rule while complicated plans appear in logic of basic fractional geometry. The CD of a constant type function is correspondent to zero wherever it is not the topic for Riemann derivative of fractional order. Mittag Leffler functions show a noteworthy authoritative in fractional-order calculus as interpretation to exponential function where the fractional-order form of exponential type function of the formft=etα/α seems in circumstance of CD.


## 3. Overview of the Method

In this division, we have considered our proposed advancedexp−φξ-expansion scheme stepwise in brief. Assume a nonlinear time-fractional NPD equation in subsequent form,
(6)RΠ,Πx,TtθΠ,Πxx,Ttt2θΠ,Πxxx,⋯⋯=0,where Π=Πx,t is defined as an unknown type function and R denotes the polynomial of Π. It is a diverse type of PD, in through the nonlinear segment and the uppermost order of differential are involved.Step 1.
We deliberate a traveling variable for transferring the mentioned equation to nondimensionality. We alter all self-governing variable into a single variable, by the following way.(7)Πx,t=uξ,ξ=kxηη±Vtθθ.
By employing mentioned variable, Equation (7) allows us in dropping Equation (6) in ODE for Πx,t=uξ in the type
(8)P.⋯,u‴,u″,u′,u,=0.Step 2.
Let us consider that a polynomial can begin the solution of OD Equation (8) in the
(9)exp−φξasu=∑i=0NAiexp−φξi,AN≠0.
HereN is the positive integer, which can be integrated by consistency of the highest order of derivatives to the highest order of nonlinear portions, performed in Equation (8).
Also, the differential ofφξ satisfies ODE in the subsequent form
(10)φ′ξ−μexp−φξ−λexpφξ=0,and then, the function solutions of ODE Equation (8) are of the form.Case 1.
The function solutions are hyperbolic (whenλμ<0):
(11)φξ=lnλ−μtanh−λμξ+C,φξ=lnλ−μcoth−λμξ+C.Case 2.
The function solutions are trigonometric (whenλμ>0):
(12)φξ=lnλμtanλμξ+C,φξ=ln−λμcotλμξ+C.Case 3.
Whenμ>0 and λ=0,
(13)φξ=ln1−μξ+C.Case 4.
Whenμ=0 and λ∈R,
(14)φξ=lnλξ+C,
whereC is defined as constant, λμ<0, or λμ>0 is conditional on the symbol of μ.Step 3.
By submitting Equation (9) into Equation (8) and finally consuming Equation (6), assemble all like term of order of exp−qφξ, q=0,±1,±2,±3,..⋯ prearranged; then, we achieve a polynomial form of exp−qφξ and connecting each coefficient of the gotten polynomial corresponding to zero, yields system of algebraic equations (SAE).Step 4.
Assume the constants’ determination be attained as one or more results by determining the mathematical conditions in phase 3. Putting the constant designs laterally with the measures for Equation (7), from the evaluation equation, we can attain up-to-date, extensive, and exhaustive dynamic wave propagation (6).


## 4. Application of the Suggested Method

### 4.1. Application of2+1 CBS Equation

Consider the nonlinear time-fractional CBS equation in a conformable sense as follows (see, for example, [17]):
(15)uxNtθu+4uxuxy+2uxxuy+uxxxy=0,t>0andx,y∈R.Hereθ is a fractional parameter in a conformable sense with the interval 0<θ≤1. Here Ntθu is conformable time fractional derivative of order θ. ux and uy are dispersive terms, where
(16)ux,y,t=uξ,ξ=x+y−wtθθ.Now using this voyaging wave variableξ=x+y−wtθ/θ into Equation (15) and integrating with respect to ξ, we get the following transformation form of ordinary differential equations with the value u=uξ.
(17)−wu′+α+β2u′2+u‴=0,where primes indicate the differentiating with respect to ξ.Now with the virtue of the homogeneous balancing method ofu‴ and u′2, we get the value of N as 1. As a result, Eq.(9) will term into the following form:
(18)uξ=A0+A1exp−φξ.Hereafter differentiating Equation (18) regarding ξ and putting the value u,u′,u‴ into Equation (17), we obtain the same polynomial set.Finally, equating the coefficientse−iϕξ equal to zero where i=0,±1,±2,±3,⋯⋯, we get some system of equations as follows.
(19)12A12βμ2+12A12αμ2+wA1μ−2A1λμ2=0,αλμA12+βλμA12−8λ2μA1+λwA1=0,−6A1λ3+12A12βλ2+12A12αλ2=0.Now solving the system of equations, we get one set of solutions.Set 1(20)w=−4λμ,A0=A0,A1=12λα+β.Case 1.
Whenλμ<0, the following solutions are obtained for hyperbolic functions:
Family 1(21)u1x,y,t=A0+12λcosh−λμξcosh−λμC+sinh−λμξsinh−λμCα+β−λ/μsinh−λμξcosh−λμC+cosh−λμξsinh−λμC,where, ξ=x+y+4λμwtθ/θ.
Family 2(22)u2x,y,t=A0+12λsinh−λμξcosh−λμC+cosh−λμξsinh−λμCα+β−λ/μcosh−λμξcosh−λμC+sinh−λμξsinh−λμC,where ξ=x+y+4λμtθ/θ.Case 2.
The function solutions are trigonometric whenλμ>0.
Family 1(23)u3x,y,t=A0+12λα+βλ/μtanλμC+tanλμξ−12λtanλμCtanλμξα+βλ/μtanλμC+tanλμξ,where ξ=x+y+4λμtθ/θ.
Family 2(24)u4x,y,t=A0−12λcotλμCα+βλ/μ−1+cotλμCcotλμξ−12λcotλμξα+βλ/μ−1+cotλμCcotλμξ,where ξ=x+y+4λμtθ/θ.Case 3.
Whenμ>0 and λ=0, the solution is not sufficient, so we can neglect this case.Case 4.
u
5
x
,
y
,
t
=
A
0
+
12
/
α
+
β
ξ
+
C, where, ξ=x+y and C is an arbitrary constant.


### 4.2. Application of RLW Equation

In this section, we applied our mentioned method to the time-fractional RLW equation as follows [31].
(25)Ntθu+ux+αu2x−βuxxNtθu=0,t>0andx,y∈R.Hereθ is the fractional constant with the interval 0<θ≤1, where
(26)ux,t=uξ,ξ=mx+ntθθ.By utilizing this wave variable of Equation (26) into Equation (25) and integrating with respect to ξ, we get this ordinary differential equation form.
(27)n+mu+αu2−βm2n2u″=0,where u′ represents the differentiating of u with respect to ξ.With the homogeneous balancing of the highest order nonlinear termu2 and the highest linear term u″, we find the value of N2. As an outcome, ?added value=" Eq.(9)"?>can be written in the following form.
(28)uξ=A0+A1−φξ+A2−φξ2.Hereafter, we differentiate Equation (28) regarding ξ and putting the required value u,u2,u″ in Eq.(17).Therefore, we finally get some polynomials and equate the coefficientse−iϕξ equal to zero, where i=0,±1,±2,±3,⋯⋯, and we get some system of equations as follows.
(29)−2βm2n2A2μ2+αA02+mA0+nA0=0,−2βλm2μn2A1+2αA0A1+mA1+nA1=0,−8βλm2μn2A2+2αA0A2+αA12+mA2+nA2=0,−2βλ2m2n2A1+2αA1A2=0,−6βλ2m2n2A2+αA22=0.Now by solving this system of equations, we find two sets of solutions as below.Set 1(30)m=±18±1+16βnλμn2+1/βλμn2α,m=n,A0=12±1/8±1+16βnλμn2+1/βλμn2+nα,A1=0,A2=32±1/8±1+16βn3λμ+1/βλμn2+nλμα.Set 2(31)m=±18±1+−16βnλμn2+1/βλμn2α,A0=−32±1/8±1+−16βn3μλ+1/βλμn2+nα,A2=−32±1/8±1+−16βnλn2+1/βλμn2+nλμα.Case 1.
whenλμ<0, the following solutions are obtained for hyperbolic functions:
Set 1.
Family 1(32)u6,7x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1−λ/μtanh−λμξ+C−32±1/8±1+16βλμn3+1/βλμn2+nαtanh−λμξ+C2,where ξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(33)ξ=±18±1+16βλμn3+1βλμn2x+ntθθ.ξ=±18±1+16βλμn3+1βλμn2x+ntθθ,u8,9x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1−λ/μcoth−λμξ+C−32±1/8±1+16βλμn3+1/βλμn2+nαcoth−λμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Set 2
Family 1(34)u10,11x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1−λ/μtanh−λμξ+C+32±1/8±1+−16βλμn3+1/βλμn2+nαtanh−λμξ+C2,where ξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(35)u12,13x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1−λ/μcoth−λμξ+C+32±1/8±1+−16βλμn3+1/βλμn2+nαcoth−λμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.Case 2.
The function solutions are trigonometric whenλμ>0.
Set 1
Family 1(36)u14,15x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1λ/μtanλμξ+C+32±1/8±1+16βλμn3+1/βλμn2+nαtanλμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(37)u16,17x,t=12±1/8±1+16βλμn3+1/βλμn2+nα−A1λ/μcotλμξ+C+32±1/8±1+16βλμn3+1/βλμn2+nαcotλμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Set 2
Family 1(38)u18,19x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1λ/μtanλμξ+C−32±1/8±1+−16βλμn3+1/βλμn2+nαtanλμξ+C2,
whereξ=±1/8±1+−16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(39)u20,21x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα−A1λ/μcotλμξ+C−32±1/8±1+−16βλμn3+1/βλμn2+nαcotλμξ+C2,
whereξ=±1/8±1+−16βλμn3+1/βλμn2x+ntθ/θ.


## 4.1. Application of2+1 CBS Equation

Consider the nonlinear time-fractional CBS equation in a conformable sense as follows (see, for example, [17]):
(15)uxNtθu+4uxuxy+2uxxuy+uxxxy=0,t>0andx,y∈R.Hereθ is a fractional parameter in a conformable sense with the interval 0<θ≤1. Here Ntθu is conformable time fractional derivative of order θ. ux and uy are dispersive terms, where
(16)ux,y,t=uξ,ξ=x+y−wtθθ.Now using this voyaging wave variableξ=x+y−wtθ/θ into Equation (15) and integrating with respect to ξ, we get the following transformation form of ordinary differential equations with the value u=uξ.
(17)−wu′+α+β2u′2+u‴=0,where primes indicate the differentiating with respect to ξ.Now with the virtue of the homogeneous balancing method ofu‴ and u′2, we get the value of N as 1. As a result, Eq.(9) will term into the following form:
(18)uξ=A0+A1exp−φξ.Hereafter differentiating Equation (18) regarding ξ and putting the value u,u′,u‴ into Equation (17), we obtain the same polynomial set.Finally, equating the coefficientse−iϕξ equal to zero where i=0,±1,±2,±3,⋯⋯, we get some system of equations as follows.
(19)12A12βμ2+12A12αμ2+wA1μ−2A1λμ2=0,αλμA12+βλμA12−8λ2μA1+λwA1=0,−6A1λ3+12A12βλ2+12A12αλ2=0.Now solving the system of equations, we get one set of solutions.Set 1(20)w=−4λμ,A0=A0,A1=12λα+β.Case 1.
Whenλμ<0, the following solutions are obtained for hyperbolic functions:
Family 1(21)u1x,y,t=A0+12λcosh−λμξcosh−λμC+sinh−λμξsinh−λμCα+β−λ/μsinh−λμξcosh−λμC+cosh−λμξsinh−λμC,where, ξ=x+y+4λμwtθ/θ.
Family 2(22)u2x,y,t=A0+12λsinh−λμξcosh−λμC+cosh−λμξsinh−λμCα+β−λ/μcosh−λμξcosh−λμC+sinh−λμξsinh−λμC,where ξ=x+y+4λμtθ/θ.Case 2.
The function solutions are trigonometric whenλμ>0.
Family 1(23)u3x,y,t=A0+12λα+βλ/μtanλμC+tanλμξ−12λtanλμCtanλμξα+βλ/μtanλμC+tanλμξ,where ξ=x+y+4λμtθ/θ.
Family 2(24)u4x,y,t=A0−12λcotλμCα+βλ/μ−1+cotλμCcotλμξ−12λcotλμξα+βλ/μ−1+cotλμCcotλμξ,where ξ=x+y+4λμtθ/θ.Case 3.
Whenμ>0 and λ=0, the solution is not sufficient, so we can neglect this case.Case 4.
u
5
x
,
y
,
t
=
A
0
+
12
/
α
+
β
ξ
+
C, where, ξ=x+y and C is an arbitrary constant.


## 4.2. Application of RLW Equation

In this section, we applied our mentioned method to the time-fractional RLW equation as follows [31].
(25)Ntθu+ux+αu2x−βuxxNtθu=0,t>0andx,y∈R.Hereθ is the fractional constant with the interval 0<θ≤1, where
(26)ux,t=uξ,ξ=mx+ntθθ.By utilizing this wave variable of Equation (26) into Equation (25) and integrating with respect to ξ, we get this ordinary differential equation form.
(27)n+mu+αu2−βm2n2u″=0,where u′ represents the differentiating of u with respect to ξ.With the homogeneous balancing of the highest order nonlinear termu2 and the highest linear term u″, we find the value of N2. As an outcome, ?added value=" Eq.(9)"?>can be written in the following form.
(28)uξ=A0+A1−φξ+A2−φξ2.Hereafter, we differentiate Equation (28) regarding ξ and putting the required value u,u2,u″ in Eq.(17).Therefore, we finally get some polynomials and equate the coefficientse−iϕξ equal to zero, where i=0,±1,±2,±3,⋯⋯, and we get some system of equations as follows.
(29)−2βm2n2A2μ2+αA02+mA0+nA0=0,−2βλm2μn2A1+2αA0A1+mA1+nA1=0,−8βλm2μn2A2+2αA0A2+αA12+mA2+nA2=0,−2βλ2m2n2A1+2αA1A2=0,−6βλ2m2n2A2+αA22=0.Now by solving this system of equations, we find two sets of solutions as below.Set 1(30)m=±18±1+16βnλμn2+1/βλμn2α,m=n,A0=12±1/8±1+16βnλμn2+1/βλμn2+nα,A1=0,A2=32±1/8±1+16βn3λμ+1/βλμn2+nλμα.Set 2(31)m=±18±1+−16βnλμn2+1/βλμn2α,A0=−32±1/8±1+−16βn3μλ+1/βλμn2+nα,A2=−32±1/8±1+−16βnλn2+1/βλμn2+nλμα.Case 1.
whenλμ<0, the following solutions are obtained for hyperbolic functions:
Set 1.
Family 1(32)u6,7x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1−λ/μtanh−λμξ+C−32±1/8±1+16βλμn3+1/βλμn2+nαtanh−λμξ+C2,where ξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(33)ξ=±18±1+16βλμn3+1βλμn2x+ntθθ.ξ=±18±1+16βλμn3+1βλμn2x+ntθθ,u8,9x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1−λ/μcoth−λμξ+C−32±1/8±1+16βλμn3+1/βλμn2+nαcoth−λμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Set 2
Family 1(34)u10,11x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1−λ/μtanh−λμξ+C+32±1/8±1+−16βλμn3+1/βλμn2+nαtanh−λμξ+C2,where ξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(35)u12,13x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1−λ/μcoth−λμξ+C+32±1/8±1+−16βλμn3+1/βλμn2+nαcoth−λμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.Case 2.
The function solutions are trigonometric whenλμ>0.
Set 1
Family 1(36)u14,15x,t=12±1/8±1+16βλμn3+1/βλμn2+nα+A1λ/μtanλμξ+C+32±1/8±1+16βλμn3+1/βλμn2+nαtanλμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(37)u16,17x,t=12±1/8±1+16βλμn3+1/βλμn2+nα−A1λ/μcotλμξ+C+32±1/8±1+16βλμn3+1/βλμn2+nαcotλμξ+C2,
whereξ=±1/8±1+16βλμn3+1/βλμn2x+ntθ/θ.
Set 2
Family 1(38)u18,19x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα+A1λ/μtanλμξ+C−32±1/8±1+−16βλμn3+1/βλμn2+nαtanλμξ+C2,
whereξ=±1/8±1+−16βλμn3+1/βλμn2x+ntθ/θ.
Family 2(39)u20,21x,t=−32±1/8±1+−16βλμn3+1/βλμn2+nα−A1λ/μcotλμξ+C−32±1/8±1+−16βλμn3+1/βλμn2+nαcotλμξ+C2,
whereξ=±1/8±1+−16βλμn3+1/βλμn2x+ntθ/θ.


## 5. Results and Discussion

Physical and graphical representation of newly established solutions of the2+1-dimensional CBS and RLW equation will be well-defined in this segment. Pictorial representation is the absolute means of envisioning the complete tangible sense of real-life concerns. We also showed the graphical depiction of our resultant solutions exhausting the computational Maple features to designate the fractional-order derivative by selecting suitable fractional values of θ. We have assigned specific high costs to the unidentified parameters to envision the consequent solutions’ natural features. The established equations are presented in Figures 1–13, which have been designed.Figure 1
Dynamical behavior of function solution ofu1x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 2
Dynamical behavior of function solution ofu2x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 3
Dynamical behavior of function solution ofu3x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 4
Dynamical behavior of function solution ofu4x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 5
Dynamical behavior of function solution ofu5x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 6
Dynamical behavior of function solution ofu6,7x,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 7
Dynamical behavior of function solution ofu8,9x,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 8
Dynamical behavior of function solution ofu10,11x,y,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 9
Dynamical behavior of function solution ofu12,13x,t. (a, b) 3D plot and (c) 2D combined stripe plot for x=1.
(a)
(b)
(c)Figure 10
Dynamical behavior of function solution ofu14,15x,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)Figure 11
Dynamical behavior of function solution ofu16,17x,t. (a, b) 3D plot and (c) 2D combined stripe plot for x=1.
(a)
(b)
(c)Figure 12
Dynamical behavior of function solution ofu18,19x,y,t. (a, b) 3D plot view representation and (c) 2D combined plot for x=1, respectively.
(a)
(b)
(c)Figure 13
Dynamical behavior of function solution ofu20,21x,t. (a, b) 3D plot and (c) 2D combined plot for x=1.
(a)
(b)
(c)


### 5.1. Physical Explanation

This subsection will focus on the physical elucidation of the2+1-dimensional CBS and RLW equation exact solutions by utilizing the advanced exp−φξ-expansion scheme. The attained solutions u1x,y,t, u2x,y,t, , u8,9x,t, u10,11x,t, and u12,13x,t are hyperbolic function solutions, u3x,y,t, u4x,y,t, u14,15x,t, u6,7x,tu16,17x,t, u18,19x,t, and u20,21x,t are trigonometric function solutions, and u5x,t is an algebraic solution. Figures 1–5 are presented from CBS equation. Figures 6–13 are from the RLW equation. In every set of figures, we have changed the free parameter of our preferred equation and fractional-order derivative θ=0.25to1 with other suitable parameters to show the dynamical behavior of attained solutions. Figure 1 represents the singular soliton shape u1x,y,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−2, C=−0.5, A=1, α=1, and β=1; with the change of fractional parameter, θ=0.25to1, the solution shape changes its amplitude and size only, but there have no changes for the variety of free parameters α and β. Figure 2 represents the bright soliton solution shape u2x,y,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−2, C=−0.5, A=1, α=1, and β=1. With the change of fractional parameter, θ=0.25to1, the solution shape turns into the dark soliton shape. Here for the evolution of free parameters α and β, the solution shape changes its amplitude and size only. Figures 3 and 4 represent the periodic singular solution shape u3x,y,t and u4x,y,t for the fractional value θ=0.25 with the suitable parameters. With the variety of fractional parameters θ=0.25to1 and free parameters, the solution has no dynamical change but only its amplitude and size. Figure 5 represents the singular soliton solution shape u5x,y,t for the fractional value θ=0.25 with the suitable parameters λ=2, μ=0, C=−0.5, A=1, α=1, and β=1. With the variety of fractional parameters θ=0.25to1 and free parameters, the solution has no dynamical change. Figures 6–9 belong to the set-1 solution set, and Figures 10–13 belong to the set-2 solution set of the RLW equation. Figure 6 represents the two soliton solution shape of u6x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the two-soliton shape turns into the one-soliton shape and has no change with free parameters. Figure 7 represents the one-soliton solution shape u7x,t for the fractional value θ=0.25 with the right parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. Shape moves to W-shaped solution shape for the transformation of fractional-order derivative θ=0.25to1. Herewith, with the change of free parameters α and β, the solution shape changes its amplitude and size only. Figures 8 and 9 represent the singular soliton shape of u14,15x,t and u16,17x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the singular soliton shape turns into the two-soliton shape and has no change with free parameters. Figure 10 represents the rogue wave soliton shape u10,11x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameter, θ=0.25to1, the rogue wave turns into the singular soliton shape. Figure 11 represents the rogue wave soliton shape u12,13x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameter, θ=0.25to1, the rogue wave changes its amplitude and size only. Figures 12 and 13 represent the singular soliton shape of u18,19x,t and u20,21x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the singular soliton shape turns into the two-soliton shape.


### 5.2. Graphical Representation

The graphical decorations of our attained solutions of the2+1-dimensional CBS and RLW equation have exposed in this subsection. Our preferred equations are resolved, and attained solutions are extracted from the trigonometric and hyperbolic function by computational tool Maple 21. For detailed physical clarity, all precise solution charts are exhibited in 3D as well as 2D combined stripe charts.


## 5.1. Physical Explanation

This subsection will focus on the physical elucidation of the2+1-dimensional CBS and RLW equation exact solutions by utilizing the advanced exp−φξ-expansion scheme. The attained solutions u1x,y,t, u2x,y,t, , u8,9x,t, u10,11x,t, and u12,13x,t are hyperbolic function solutions, u3x,y,t, u4x,y,t, u14,15x,t, u6,7x,tu16,17x,t, u18,19x,t, and u20,21x,t are trigonometric function solutions, and u5x,t is an algebraic solution. Figures 1–5 are presented from CBS equation. Figures 6–13 are from the RLW equation. In every set of figures, we have changed the free parameter of our preferred equation and fractional-order derivative θ=0.25to1 with other suitable parameters to show the dynamical behavior of attained solutions. Figure 1 represents the singular soliton shape u1x,y,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−2, C=−0.5, A=1, α=1, and β=1; with the change of fractional parameter, θ=0.25to1, the solution shape changes its amplitude and size only, but there have no changes for the variety of free parameters α and β. Figure 2 represents the bright soliton solution shape u2x,y,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−2, C=−0.5, A=1, α=1, and β=1. With the change of fractional parameter, θ=0.25to1, the solution shape turns into the dark soliton shape. Here for the evolution of free parameters α and β, the solution shape changes its amplitude and size only. Figures 3 and 4 represent the periodic singular solution shape u3x,y,t and u4x,y,t for the fractional value θ=0.25 with the suitable parameters. With the variety of fractional parameters θ=0.25to1 and free parameters, the solution has no dynamical change but only its amplitude and size. Figure 5 represents the singular soliton solution shape u5x,y,t for the fractional value θ=0.25 with the suitable parameters λ=2, μ=0, C=−0.5, A=1, α=1, and β=1. With the variety of fractional parameters θ=0.25to1 and free parameters, the solution has no dynamical change. Figures 6–9 belong to the set-1 solution set, and Figures 10–13 belong to the set-2 solution set of the RLW equation. Figure 6 represents the two soliton solution shape of u6x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the two-soliton shape turns into the one-soliton shape and has no change with free parameters. Figure 7 represents the one-soliton solution shape u7x,t for the fractional value θ=0.25 with the right parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. Shape moves to W-shaped solution shape for the transformation of fractional-order derivative θ=0.25to1. Herewith, with the change of free parameters α and β, the solution shape changes its amplitude and size only. Figures 8 and 9 represent the singular soliton shape of u14,15x,t and u16,17x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=−0.5, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the singular soliton shape turns into the two-soliton shape and has no change with free parameters. Figure 10 represents the rogue wave soliton shape u10,11x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameter, θ=0.25to1, the rogue wave turns into the singular soliton shape. Figure 11 represents the rogue wave soliton shape u12,13x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameter, θ=0.25to1, the rogue wave changes its amplitude and size only. Figures 12 and 13 represent the singular soliton shape of u18,19x,t and u20,21x,t for the fractional value θ=0.25 with the suitable parameters λ=3, μ=−0.2, C=1, n=−0.4, α=6, and β=10. With the variety of fractional parameters, θ=0.25to1, the singular soliton shape turns into the two-soliton shape.


## 5.2. Graphical Representation

The graphical decorations of our attained solutions of the2+1-dimensional CBS and RLW equation have exposed in this subsection. Our preferred equations are resolved, and attained solutions are extracted from the trigonometric and hyperbolic function by computational tool Maple 21. For detailed physical clarity, all precise solution charts are exhibited in 3D as well as 2D combined stripe charts.


## 6. Conclusions

In the above study, exact wave solutions with several types of wave constructions for time fractional2+1-dimensional CBS and RLW equation have been accumulated through the advanced exp−φξ expansion technique. The dissimilar wave structures characterized the attained solutions’ dynamical behavior in this system. Preferred fractional and free parameters have significant insinuations, such as putting the dissimilar magnitude of the fractional and free parameters from a separate function of meaningful solutions can be originated exclusively and variation of suitable parameters has a significant impact on the wave. Moreover, we claim that our obtained traveling solutions are new in the sense of conformable derivative and could be operative in the research of nonlinear physical phenomena. It is very clear to decide that our promoted technique is effective, reliable, and friendly applicable and delivers sufficient well-matched explanations to NLFEEs arise in engineering, applied mathematics, nonlinear dynamics, and mathematical physics.



---

*Source: 3031117-2022-09-09.xml*