# New Solitons and Multishock Wave Structures for the Conformable Space Fractional Burger and Time Fractional Sharma-Tasso-Olver Models

Authors: Selina Akter; Harun-Or-Roshid; N. F. M. Noor
Journal: Advances in Mathematical Physics
 (2022)

---

## Abstract

This present paper studies the conformable space fractional Burgers, and the time fractional Sharma-Tasso-Olver models; both are highly important for nonlinear diffusive waves in fluid dynamics, sound waves in a viscous medium, and flow in field soils, as well as in gas and plasma dynamics. To retrieve explicit solutions of the fractional differential models, we propose an integral scheme, namely, Modified Kudryashov method. We obtain periodic, solitary, mixed periodic-soliton, and polynomial solutions through the approach. In particular, we exhibit topological kink-dark bell wave, topological kink, singular kink, bright bell, peakon solitons, and periodic shape waves to apply suitable values on parameters for both distinct models. The impact of fractionality on the wave shape and its deformation is analyzed and discussed graphically. We also investigate multishock wave’s solutions of both models and analyzed the effect of each existing parameters involved in the obtained solutions. To visualize the real characters of the solitary solutions, the graphical elucidation in 3D and 2D profiles are plotted. In computational effort and realization, it is emphasized that the proposed scheme is friendly useful, highly effective, and a powerful mathematical tool to extract exact solitary wave solutions for the differential models, as well as fractional differential models.


---

## Body

## 1. Introduction

Many natural phenomena in the fields of engineering; ocean engineering, fluid flows, mathematical biology, signal processing, optical communications systems, and electromagnetic theories are modeled in terms of fractional derivatives [1–11]. They have attracted high curiosities among the researchers owing to their occurrences in the modeling of complex physical phenomena. Therefore, solitary wave solutions for the relevant nonlinear problems are fundamental issues. A variety of powerful procedures have been explored to evaluate the solitary solutions of the various differential models with and without fractional derivatives in literature, such as Weierstrass elliptic function method [2], Laplace perturbation method [3], extended sinh-Gordon equation expansion method [4, 5], Darboux transformation [6], modified residual power series method [7], Bernoulli polynomial approaches [8], exp-expansion method [9], Adomian decomposition method [10], homotopy perturbation method [11], generalized Kudryashov method [12–14], (G′/G)-expansion method [15], and Hirota bilinear methods [16–20].Nowadays, numerous efforts have been invested in execution of the above methods to derive exact soliton solutions for nonlinear differential models with space-time fractional. Moreover, fractional differential models have apparent huge devotion due to proper explanation of many nonlinear events [21, 22]. Event symmetry analysis, analytical solutions, and conservation laws of both total and fractional nonlinear models [22–24] have attracted more attention day by day. The space time fractional Burger (sFB) [8–10] and the time fractional Sharma-Tasso-Olver (tFSTO) [9, 11] models are among the most important fractional models to acquire deep understandings in real nonlinear phenomena in the relevant fields. Both the sFB and tFSTO models are initially converted from the generalized Burger model with some parameter evaluations which have many applications in stochastic problem, water and shock waves, gas dynamics, acoustic science, thermal conduction, etcetera.Motivated by the importance of the fractional models of sFB and tFSTO and their solutions, present investigation shed the light on the new approach for solving the space time fractional Burger (sFB) and the time fractional Sharma-Tasso-Olver (tFSTO) models by introducing Modified Kudryashov method (MKM). Moreover, we went to derive multishock soliton solutions of the models. Hence, the structure of this paper is as follows: In Section2, we incorporate definition of conformable fractional derivative and some of its properties. In Section 3, we describe details on the proposed Modified Kudryashov method (MKM) and its applications. In Section 4, we derive new plentiful analytic solutions of the sFBM and the tFSTOM models. Next, we present different types of shapes of the gained wave solutions in 3D and 2D graphics, respectively, in Section 5. Finally, the conclusions of this research are deduced in the last part.


## 2. Conformable Fractional Derivative and Its Properties

In this section, some important properties of the conformable fractional derivatives according to fractional calculus theory in [3–5] are revisited:Let us consider a functionϕ:0,∞⟶R where the conformable fractional derivative of ϕ with order δ is defined as ∂δϕ/∂tδ=limε⟶0+ϕt+εt1−δ−ϕt/ε, t>0 and 0<δ≤1.Hence, few key properties of the fractional derivative are given as:(1)∂δ∂tδlϕ+mψ=l∂δ∂tδϕ+m∂δ∂tδψ,∀l,m∈R,∂δ∂tδtβ=βtβ−δ,∀β∈R,and
(2)∂δ∂tδλ=0,λ=const.,∂δ∂tδϕ∘ψt=t1−δϕ′ψtψ′t.


## 3. The Proposed Modified Kudryashov Method (MKM)

Here with onwards, we will extend the Kudryashov method by using a new suitable integration technique so that our modified method of Kudryashov can further be applied to any nonlinear evolution equations. Let us consider a general form of any nonlinear model as(3)QU,DxαDtβU,DxαDxαU,DxαU,DtβU,⋯⋯..=0,0<α,β<1.Note thatQ is a polynomial of U=Ux,t and its various partial fractional derivatives Dα include the nonlinear terms.Step-01. To convert the nonlinear partial differential Eq. (5) into ordinary differential equation (ODE), we first use the traveling variable,
(4)Ux,t=Uξ,ξ=kxαΓ1+α−ωxβΓ1+β,where k, ω are nonzero constants. Applying the chain rule
(5)DxαU=ρxdUdξDxαξ,DtβU=ρtdUdξDtβξ,where ρx,ρt are sigma indexes [3] which can be written as ρx,ρt=Kconst..Now, by inserting Eq. (4) along with (5) into Eq. (3), we can reach an ordinary differential equation as follows:
(6)Q1Uξ,Uξξξ,Uξξ,UUξ,Uξξξ,..⋯⋯=0.Next, we will integrate Eq. (6) as many times as possible while keeping every integral constant to be zero.Step-2. Assuming the solution of Eq. (6) takes the general form of a polynomial function Fξ, we can write
(7)Uξ=∑i=0NaiFξi∑j=0MbjFξj,where ai,bj,ω,k,i=0,⋯⋯,m,j=0,⋯⋯,n are constants to be evaluated later. Here, aN≠0,bM≠0 whereas N and Mare positive integers, which can be found from the principle of balancing technique. Next, we have a condition of Fξ that satisfies the general form of the Riccati equation,
(8)F′ξ=A+BFξ+CF2ξ.Step-3. By substituting Eq. (7) into Eq. (6) and using the ODE (8), we can execute a polynomial for Fξ by collecting all same order of Fξ together. Equating coefficients of Fξ yields a set of algebraic systems for ai,bj,ω,k,i=0,⋯⋯,m,j=0,⋯⋯,n. Solving the unknowns and setting them into Eq. (7), we can attain the analytic and exact traveling wave solutions of the nonlinear model in Eq. (3).It is further pointed out that diverse types of solutions are possible for (8):Case-1: WhenB=0,C=1, then (8) has the solution
(9)Fξ=−−Atanh−Aξ,A<0−−Acoth−Aξ,A<0,(10)Fξ=−1ξ,A=0,(11)Fξ=AtanAξ,A>0−AcotAξ,A>0.Case-2: WhenA=0,B=−1,C=1, then (8) has the solution
(12)Fξ=11+heξ.Case-3: WhenC=0, then (8) has the solution
(13)Fξ=heBξ−AB.In general, there may have more solutions of (8), but that will be applicable when we can solve any nonlinear models without particular values of the constants A,B,C.Remark. When we consider the particular valuesA=0,B=−1,C=1 in the auxiliary equation of the proposed method, we reach to the generalized Kudryashov method [12, 13] and when b0=1,bj=0 in the trail solution (7) of our method, then we reach to the extended Kudryashov method [14]. Thus, the methods [12–14] are the special cases of this newly proposed MKM.


## 4. Applications to the Models

In these parts, we present the applications of the MKM for obtaining exact traveling wave solutions to the space fractional Burgers and the time fractional Sharma-Tasso-Olver models. Additionally, we present multishock wave solutions by a different approach of the models.


### 4.1. The sFBM and tFSTOM Models

The sFBM model given by Eq. (14) is a nonlinear model of diffusive waves which found many applications in fluid dynamics, one-dimensional soundwaves in a viscous medium, gas dynamics, shock waves in a viscous medium, magneto-hydrodynamic waves in a medium with finite electrical conductivity, turbulence, and in plasma dynamics applied to unsaturated flow in field soils [8–10]. This model has the usual form [8–10]:
(14)νt+ννx−kνxx+qDxδν=0,x,t>0,0<δ≤1,ν0,t=0,νx0,t=1t−π22kt2,where k,q are constants and δ is an aspect relating to the order of the fractional derivative.On the other hand, the tFSTOM model has real-life applications in the fields of relativistic quantum atom dynamics, the relativistic energy evaluation, and the nonlinear electromagnetic phenomena. The tFSTOM model has the usual form [9, 11]:
(15)Dtδν+3aνx2+3aν2νx+3aννxx+aνxxx=0,t>0,0<δ<1,νx,0=−2b0tan2b02x,where a,b0 are constants and δ is an aspect relating to the order of the fractional derivative.


### 4.2. The Application of the MKM to the sFBM

In this section, we explore a reliable treatment to the sFBM with the help of the MKM.Primarily, we will use the traveling wave variable,νx,t=νξ,ξ=lxδ/Γ1+δ−wt, to renovate Eq. (14) to the following ordinary differential equation:
(16)−wν′+lνν′−kl2ν″+qlν′=0.Integrating (16) one time gives us
(17)p+lq−wν+lν2/2−kl2ν′=0,where p is an integral constant.Now we compute the balance number of Eq. (17) from the linear terms ν″ and ν3 gives n=m+1. For m=1, we have n=2. So, the trial solution of Eq. (7) takes the following form,
(18)νξ=l0+l1Fξ+l2Fξ2m0+m1Fξ.We now differentiate Eq. (18) with respect to ξ along with Eq. (8), and then inserting ν,ν′ intoEq. (17), which gives an equation. Equating the coefficients of Flξ from the required equation equal to zero yields the following set of constraints:Set-1:(19)p=2k2l3CAm12+Bm0m1−Bm12+Cm02−2Cm0m1+Cm12m12,w=lq+2Ckl−Bkl−2Cklm0/m1,l0=2Cklm0m1−m0/m1,l1=2Cklm1,l2=2Cklm1.Set-2:(20)p=2Ck2l3m1Am1−Bm022m0+m1A2Bm14+A2Cm14−2AB2m0m13+4BC2m04m0−m1+2ABCm0m124m02−m0m2−1+4AC2m03m12+2B2m02m1Bm0m1+Bm12−4Cm02+B2Cm02m13Am12−Bm0m1+2Cm024,w=lklABm12Am12−2Bm0+2A2Cm13m0+m1−2ABCm0m12m0+2m1+m02m14AC2+B3m1+2B2Cm02m12−2m02+qm12A2m12+B2m02−2Am0m12Bm1−2Cm0−4Cm02Bm1−CAm12−Bm0m1+2Cm022,l0=2klm1Am1−Bm02Bm0m1+Cm0m1−Am12−2Cm02Am12−Bm0m1+2Cm022,l1=2Cklm13Am1−Bm0Am12−Bm0m1+2Cm022,l2=2Cklm12Bm0−Am1Am12−Bm0m1+2Cm02.Set-3:(21)p=B8BC3kll2m022A+B+B3l23l2−4Cklm0−322C3k2l2m02+BCkll2m0R128C6k2lm04,w=8C3klqm02−4B2Ckll2m0+B2l228C3km02,l0=±Rm0klC3,l1=Bl224C2klm0,m1=2Cm0B,where
(22)R=−18B4C3klm024AC−B2+B3l22Cklm0−l2±2Cklm024AC−B24C4klm024AC−B2+B3l24Cklm0−l2.Setting the constraints of Set-1 into (18) along with the results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we attain the following solutions:Case-1S1: ForB=0,C=1, the gained solutions are:
(23)ν1x,t=2klm0m1−m0/m1−2klm1−Atanh−Aξ−2Aklm1tanh2−Aξm0−m1−Atanh−Aξ,A<0,(24)ν2x,t=2klm0m1−m0/m1−2klm1−Acoth−Aξ−2Aklm1coth2−Aξm0−m1−Acoth−Aξ,A<0,(25)ν3x,t=2klm0m1−m0ξ2−2klm12ξ+2klm12m1m0ξ2−m1ξ,(26)ν4x,t=2klm0m1−m0/m1+2klm1AtanAξ+2Aklm1tan2Aξm0+m1AtanAξ,A>0,(27)ν5x,t=2klm0m1−m0/m1−2klm1AcotAξ+2Aklm1cot2Aξm0−m1AcotAξ,A>0,where ξ=lxδ/Γ1+δ−lq+2kl−2klm0/m1t.Case-2S1: ForA=0,B=−1,C=1, the gained solutions are:
(28)ν6x,t=2klm1+2klm11+heξ+2klm0m1−m0/m11+heξ2m11+heξ+m01+heξ2,where ξ=lxδ/Γ1+δ−lq+2kl+kl−2klm0/m1t.Case-3S1: The third case gives the invalid results asC=0, gives l2=0.Setting the constraints of Set-2 into (18) along with results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we further attain the following solutions:Case-1S2: ForB=0,C=1, the gained solutions are:
(29)ν7x,t=2A2klm13m0m1−Am12−2m02−2A2klm15−Atanh−Aξ+2A2klm13Am12+2m02tanh2−AξAm12+2m022m0−m1−Atanh−Aξ,A<0,ν8x,t=2A2klm13m0m1−Am12−2m02−2A2klm15−Acoth−Aξ+2A2klm13Am12+2m02coth2−AξAm12+2m022m0−m1−Acoth−Aξ,A<0,ν9x,t=2A2klm13m0m1−Am12−2m02+2A2klm15AtanAξ−2A2klm13Am12+2m02tan2AξAm12+2m022m0+m1AtanAξ,A>0,ν10x,t=2A2klm13m0m1−Am12−2m02−2A2klm15AcotAξ−2A2klm13Am12+2m02cot2AξAm12+2m022m0m1AcotAξ,A>0,where ξ=lxδ/Γ1+δ−l2A2klm13m0+m1+4Aklm02m1+4qA2m14+Am02m12+m04t/Am12+2m022.Case-2S2: ForA=0,B=−1,C=1, the gained solutions are:
(30)ν11x,t=−2klm0m12m0m1+2m02+2klm02m131+heξ−4klm1m021+heξ2m0m1+2m022m11+heξ+m01+heξ2,where ξ=lxδ/Γ1+δ−lklm024m02+m12+qm024m02+4qm0m1+m12t/m0m1+2m022.Case-3S2: The third case gives the invalid results asC=0, gives l2=0.Setting the constraints of Set-3 into (18) along with results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we attain the following solutions:Case-1S3: ForB=0,C=1, Case-1 and Case-3 give the undefined statements.Case-2S3: ForA=0,B=−1,C=1, the gained solutions are:
(31)ν12x,t=4klm0l2−l221+heξ±4R1+heξ24klm021+heξheξ−1,where ξ=lxδ/Γ1+δ−8kqlm02−4kll2m0+l22/8km02t and R=l22/8−k2l2m02−1/2kll2m0orl22/8.Case-3S3: The third case gives the invalid results asC=0, gives l2=0.


### 4.3. The Multishock Waves to the Time sFBM

In this section, we now explore a reliable technique to achieve multishock wave of the fractional model. To do so, we consider the simple Burger equation(32)νt−2ννx−νxx=0.Taking traveling variableνx,t=Yξ,ξ=lx−wt, Eq. (32) leads to ODE as follows:
(33)Yξx,t=Y−1lY2.Equation (33) has a multisoliton solution [21].
(34)Y=∑j=1nljeljx−wjt1+∑j=1neljx−wjt.Let us consider the trial solution of Eq. (17) is
(35)v=a0+a1Y.asbalancenumberofBurgermodelisunityKeeping arbitrary constantp=0 in Eq. (17), differentiating Eq. (35) together with Eq. (33) into Eq. (17), leads an algebraic system and solving, we reach to multisoliton shock wave solution
(36)v=2k∑j=1nlj−2k∑j=1nδlj/Γ1+δeljxδ/Γ1+δ−wjt1+∑j=1neljxδ/Γ1+δ−wjt,where
(37)wj=∑j=1nklj2+qlj.


### 4.4. The Application of the MKM to the Time FSTOM

In this section, we now explore a reliable treatment to the time FSTOM with the help of the MKM.Firstly, we will use the traveling wave variableνx,t=νξ,ξ=x−ltδ/Γ1+δ, to renovate Eq. (4) to the following ordinary differential equation:
(38)−lν′+3aν′2+3aν2ν′+3aνν″+aν‴=0.Integrating Eq. (38) one time gives us
(39)p−lν+3aνν′+aν3+aν″=0.Now we compute the balance number of (39) from the linear terms ν″ and ν3 gives n=m+1.Form=1, we have n=2. So, the trial solution of Eq. (7) takes the following form,
(40)νξ=l0+l1Fξ+l2Fξ2m0+m1Fξ.We now differentiate Eq. (40) with respect to ξ along with (10), and then inserting v,v′,v″ into Eq. (39), it gives an equation. Equating the coefficients of Flξ from the required equation equal to zero yields the following findings:It is very difficult to solve the obtained system of equations in general case even in computational software. Thus, we choose the particular cases (Case-1 and Case-2) to get the results as follows:Case-1: ForB=0,C=1, the system of equation leads to the solutions:Set-1:p=0,l=−4aA,l0=0,l2=l1,m0=−l1/2,m1=−l1/2.Set-2:p=0,l=−aA,l0=0,l2=l1,m0=−l1,m1=−l1.Setting the constraints of Set-1 into (40) along with the results of auxiliary equation from Case-1 of our method, we attain the following solutions:Case-1S1: ForB=0,C=1, the gained solutions are:
(41)ν1x,t=2−Atanh−Ax+4aAtδ/Γ1+δ+2Atanh2−Ax+4aAtδ/Γ1+δ1−−Atanh−Ax+4aAtδ/Γ1+δ,A<0,(42)ν2x,t=2−Acoth−Ax+4aAtδ/Γ1+δ+2Acoth2−Ax+4aAtδ/Γ1+δ1−−Acoth−Ax+4aAtδ/Γ1+δ,A<0,(43)ν3x,t=−2AtanAx+4aAtδ/Γ1+δ−2Atan2Ax+4aAtδ/Γ1+δ1+AtanAx+4aAtδ/Γ1+δ,A>0,(44)ν4x,t=2AcotAx+4aAtδ/Γ1+δ−2Acot2Ax+4aAtδ/Γ1+δ1−AcotAx+4aAtδ/Γ1+δ,A>0.Setting the constraints of Set-2 into (40) along with results of auxiliary equation from Case-1 of our method, we attain to the following solutions:Case-1S2: ForB=0,C=1, the solutions are:
(45)ν5x,t=−Atanh−Ax+aAtδ/Γ1+δ+Atanh2−Ax+aAtδ/Γ1+δ1−−Atanh−Ax+aAtδ/Γ1+δ,A<0,ν6x,t=−Acoth−Ax+aAtδ/Γ1+δ+Acoth2−Ax+aAtδ/Γ1+δ1−−Acoth−Ax+aAtδ/Γ1+δ,A<0,ν7x,t=−AtanAx+aAtδ/Γ1+δ−Atan2Ax+aAtδ/Γ1+δ1+AtanAx+aAtδ/Γ1+δ,A>0,ν8x,t=AcotAx+aAtδ/Γ1+δ−Acot2Ax+aAtδ/Γ1+δ1−AcotAx+aAtδ/Γ1+δ,A>0.Case-2: ForA=0,B=−1,C=1, the system of equations leads to the solutions:Set-1:(46)p=1m122m0m1+m123am02−m0m1−6am02+12am0m1−7am12m12+8am0m1−m02m1+am12+10m0m1−8m02,l=3am0m1−m02+6am02−12am0m1+7am12m12,l0=m02−m0m1m1,l1=−m1,l2=−m1.Set-2:p=0,l=a,l0=3m1/2,l1=−2m1,l2=−2m1,m0=3m1/2.Set-3:p=−12a,l=13a,l0=m1/2,l1=−m1,l2=−m1,m0=−m1/2.Set-4:p=0,l=a,l0=m1/2,l1=−m1,l2=−m1,m0=m1/2.Set-5:(47)p=0,l=−a275−3/2±17/23+1815−3/2±17/22+2412−3/2±17/23−200422−3/2±17/2+78,l0=−m13−3/2±17/2+1111−3/2±17/2+39,l1=−32±172m1,l2=−2m1−3/2±17/2+3,m0=−m1/2.Setting the constraints of Set-1 to Set-5 into (40) along with the results of auxiliary equation from Case-2 yields the solutions:
(48)ν9x,t=m02−m0m11+heξ2−m122+heξm11+heξm01+heξ+m1,(49)ξ=x−3am0m1−m02+6am02−12am0m1+7am12tδm12Γ1+δ,(50)ν10x,t=−8−4hexpx−atδ/Γ1+δ+1+hexpx−atδ/Γ1+δ21+hexpx−atδ/Γ1+δ5+3hexpx−atδ/Γ1+δ,(51)ν11x,t=−4−2hexpx+12atδ/Γ1+δ+1+hexpx+12atδ/Γ1+δ21−h2exp2x+24atδ/Γ1+δ,(52)ν12x,t=−4−2hexpx−atδ/Γ1+δ+1+hexpx−atδ/Γ1+δ21+hexpx−atδ/Γ1+δ3+hexpx−atδ/Γ1+δ,(53)ν13x,t=−6L+22/11L+391+hexpx−ltδ/Γ1+δ2+2L1+hexpx−ltδ/Γ1+δ−4/L+31−1+hexpx−ltδ/Γ1+δ2,where l=−a275L3+1815L2+2412L−2004/22L+78 and L=−3/2±17/2.


### 4.5. The Multishock Waves to the FSTOM

In this section, we now explore a reliable technique [25] to achieve multishock wave of the fractional FSTOM. To do so, we consider the simple Burger equation (32) and Eq. (33). Then, we consider the wave transformation
(54)νx,t=νξ,ξ=lx−wtδΓ1+δ.Using Eq. (54), Eq. (4) leads to
(55)−wv′+3al2v′2+3alv2v′+3al2vv′′+al3v′′′=0.Integrating once with putting integral constant as zero, Eq. (55) provides
(56)–wv+3al2vν′+alν3+al3v′′=0.Let us consider the trial solution of Eq. (56) is
(57)v=a0+a1Y.asbalancenumberofthemodelisunityDifferentiating Eq. (57) together with Eq. (33) into Eq. (56), leads an algebraic system and solving, we reach to mutlishock wave solution.
(58)v=−∑j=1nlj+2∑j=1nljeljx−wtδ/Γ1+δ1+∑j=1neljx−wtδ/Γ1+δ,where
(59)wj=∑j=1nalj3.


## 4.1. The sFBM and tFSTOM Models

The sFBM model given by Eq. (14) is a nonlinear model of diffusive waves which found many applications in fluid dynamics, one-dimensional soundwaves in a viscous medium, gas dynamics, shock waves in a viscous medium, magneto-hydrodynamic waves in a medium with finite electrical conductivity, turbulence, and in plasma dynamics applied to unsaturated flow in field soils [8–10]. This model has the usual form [8–10]:
(14)νt+ννx−kνxx+qDxδν=0,x,t>0,0<δ≤1,ν0,t=0,νx0,t=1t−π22kt2,where k,q are constants and δ is an aspect relating to the order of the fractional derivative.On the other hand, the tFSTOM model has real-life applications in the fields of relativistic quantum atom dynamics, the relativistic energy evaluation, and the nonlinear electromagnetic phenomena. The tFSTOM model has the usual form [9, 11]:
(15)Dtδν+3aνx2+3aν2νx+3aννxx+aνxxx=0,t>0,0<δ<1,νx,0=−2b0tan2b02x,where a,b0 are constants and δ is an aspect relating to the order of the fractional derivative.


## 4.2. The Application of the MKM to the sFBM

In this section, we explore a reliable treatment to the sFBM with the help of the MKM.Primarily, we will use the traveling wave variable,νx,t=νξ,ξ=lxδ/Γ1+δ−wt, to renovate Eq. (14) to the following ordinary differential equation:
(16)−wν′+lνν′−kl2ν″+qlν′=0.Integrating (16) one time gives us
(17)p+lq−wν+lν2/2−kl2ν′=0,where p is an integral constant.Now we compute the balance number of Eq. (17) from the linear terms ν″ and ν3 gives n=m+1. For m=1, we have n=2. So, the trial solution of Eq. (7) takes the following form,
(18)νξ=l0+l1Fξ+l2Fξ2m0+m1Fξ.We now differentiate Eq. (18) with respect to ξ along with Eq. (8), and then inserting ν,ν′ intoEq. (17), which gives an equation. Equating the coefficients of Flξ from the required equation equal to zero yields the following set of constraints:Set-1:(19)p=2k2l3CAm12+Bm0m1−Bm12+Cm02−2Cm0m1+Cm12m12,w=lq+2Ckl−Bkl−2Cklm0/m1,l0=2Cklm0m1−m0/m1,l1=2Cklm1,l2=2Cklm1.Set-2:(20)p=2Ck2l3m1Am1−Bm022m0+m1A2Bm14+A2Cm14−2AB2m0m13+4BC2m04m0−m1+2ABCm0m124m02−m0m2−1+4AC2m03m12+2B2m02m1Bm0m1+Bm12−4Cm02+B2Cm02m13Am12−Bm0m1+2Cm024,w=lklABm12Am12−2Bm0+2A2Cm13m0+m1−2ABCm0m12m0+2m1+m02m14AC2+B3m1+2B2Cm02m12−2m02+qm12A2m12+B2m02−2Am0m12Bm1−2Cm0−4Cm02Bm1−CAm12−Bm0m1+2Cm022,l0=2klm1Am1−Bm02Bm0m1+Cm0m1−Am12−2Cm02Am12−Bm0m1+2Cm022,l1=2Cklm13Am1−Bm0Am12−Bm0m1+2Cm022,l2=2Cklm12Bm0−Am1Am12−Bm0m1+2Cm02.Set-3:(21)p=B8BC3kll2m022A+B+B3l23l2−4Cklm0−322C3k2l2m02+BCkll2m0R128C6k2lm04,w=8C3klqm02−4B2Ckll2m0+B2l228C3km02,l0=±Rm0klC3,l1=Bl224C2klm0,m1=2Cm0B,where
(22)R=−18B4C3klm024AC−B2+B3l22Cklm0−l2±2Cklm024AC−B24C4klm024AC−B2+B3l24Cklm0−l2.Setting the constraints of Set-1 into (18) along with the results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we attain the following solutions:Case-1S1: ForB=0,C=1, the gained solutions are:
(23)ν1x,t=2klm0m1−m0/m1−2klm1−Atanh−Aξ−2Aklm1tanh2−Aξm0−m1−Atanh−Aξ,A<0,(24)ν2x,t=2klm0m1−m0/m1−2klm1−Acoth−Aξ−2Aklm1coth2−Aξm0−m1−Acoth−Aξ,A<0,(25)ν3x,t=2klm0m1−m0ξ2−2klm12ξ+2klm12m1m0ξ2−m1ξ,(26)ν4x,t=2klm0m1−m0/m1+2klm1AtanAξ+2Aklm1tan2Aξm0+m1AtanAξ,A>0,(27)ν5x,t=2klm0m1−m0/m1−2klm1AcotAξ+2Aklm1cot2Aξm0−m1AcotAξ,A>0,where ξ=lxδ/Γ1+δ−lq+2kl−2klm0/m1t.Case-2S1: ForA=0,B=−1,C=1, the gained solutions are:
(28)ν6x,t=2klm1+2klm11+heξ+2klm0m1−m0/m11+heξ2m11+heξ+m01+heξ2,where ξ=lxδ/Γ1+δ−lq+2kl+kl−2klm0/m1t.Case-3S1: The third case gives the invalid results asC=0, gives l2=0.Setting the constraints of Set-2 into (18) along with results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we further attain the following solutions:Case-1S2: ForB=0,C=1, the gained solutions are:
(29)ν7x,t=2A2klm13m0m1−Am12−2m02−2A2klm15−Atanh−Aξ+2A2klm13Am12+2m02tanh2−AξAm12+2m022m0−m1−Atanh−Aξ,A<0,ν8x,t=2A2klm13m0m1−Am12−2m02−2A2klm15−Acoth−Aξ+2A2klm13Am12+2m02coth2−AξAm12+2m022m0−m1−Acoth−Aξ,A<0,ν9x,t=2A2klm13m0m1−Am12−2m02+2A2klm15AtanAξ−2A2klm13Am12+2m02tan2AξAm12+2m022m0+m1AtanAξ,A>0,ν10x,t=2A2klm13m0m1−Am12−2m02−2A2klm15AcotAξ−2A2klm13Am12+2m02cot2AξAm12+2m022m0m1AcotAξ,A>0,where ξ=lxδ/Γ1+δ−l2A2klm13m0+m1+4Aklm02m1+4qA2m14+Am02m12+m04t/Am12+2m022.Case-2S2: ForA=0,B=−1,C=1, the gained solutions are:
(30)ν11x,t=−2klm0m12m0m1+2m02+2klm02m131+heξ−4klm1m021+heξ2m0m1+2m022m11+heξ+m01+heξ2,where ξ=lxδ/Γ1+δ−lklm024m02+m12+qm024m02+4qm0m1+m12t/m0m1+2m022.Case-3S2: The third case gives the invalid results asC=0, gives l2=0.Setting the constraints of Set-3 into (18) along with results of auxiliary equation from Case-1, Case-2, and Case-3 of our method, we attain the following solutions:Case-1S3: ForB=0,C=1, Case-1 and Case-3 give the undefined statements.Case-2S3: ForA=0,B=−1,C=1, the gained solutions are:
(31)ν12x,t=4klm0l2−l221+heξ±4R1+heξ24klm021+heξheξ−1,where ξ=lxδ/Γ1+δ−8kqlm02−4kll2m0+l22/8km02t and R=l22/8−k2l2m02−1/2kll2m0orl22/8.Case-3S3: The third case gives the invalid results asC=0, gives l2=0.


## 4.3. The Multishock Waves to the Time sFBM

In this section, we now explore a reliable technique to achieve multishock wave of the fractional model. To do so, we consider the simple Burger equation(32)νt−2ννx−νxx=0.Taking traveling variableνx,t=Yξ,ξ=lx−wt, Eq. (32) leads to ODE as follows:
(33)Yξx,t=Y−1lY2.Equation (33) has a multisoliton solution [21].
(34)Y=∑j=1nljeljx−wjt1+∑j=1neljx−wjt.Let us consider the trial solution of Eq. (17) is
(35)v=a0+a1Y.asbalancenumberofBurgermodelisunityKeeping arbitrary constantp=0 in Eq. (17), differentiating Eq. (35) together with Eq. (33) into Eq. (17), leads an algebraic system and solving, we reach to multisoliton shock wave solution
(36)v=2k∑j=1nlj−2k∑j=1nδlj/Γ1+δeljxδ/Γ1+δ−wjt1+∑j=1neljxδ/Γ1+δ−wjt,where
(37)wj=∑j=1nklj2+qlj.


## 4.4. The Application of the MKM to the Time FSTOM

In this section, we now explore a reliable treatment to the time FSTOM with the help of the MKM.Firstly, we will use the traveling wave variableνx,t=νξ,ξ=x−ltδ/Γ1+δ, to renovate Eq. (4) to the following ordinary differential equation:
(38)−lν′+3aν′2+3aν2ν′+3aνν″+aν‴=0.Integrating Eq. (38) one time gives us
(39)p−lν+3aνν′+aν3+aν″=0.Now we compute the balance number of (39) from the linear terms ν″ and ν3 gives n=m+1.Form=1, we have n=2. So, the trial solution of Eq. (7) takes the following form,
(40)νξ=l0+l1Fξ+l2Fξ2m0+m1Fξ.We now differentiate Eq. (40) with respect to ξ along with (10), and then inserting v,v′,v″ into Eq. (39), it gives an equation. Equating the coefficients of Flξ from the required equation equal to zero yields the following findings:It is very difficult to solve the obtained system of equations in general case even in computational software. Thus, we choose the particular cases (Case-1 and Case-2) to get the results as follows:Case-1: ForB=0,C=1, the system of equation leads to the solutions:Set-1:p=0,l=−4aA,l0=0,l2=l1,m0=−l1/2,m1=−l1/2.Set-2:p=0,l=−aA,l0=0,l2=l1,m0=−l1,m1=−l1.Setting the constraints of Set-1 into (40) along with the results of auxiliary equation from Case-1 of our method, we attain the following solutions:Case-1S1: ForB=0,C=1, the gained solutions are:
(41)ν1x,t=2−Atanh−Ax+4aAtδ/Γ1+δ+2Atanh2−Ax+4aAtδ/Γ1+δ1−−Atanh−Ax+4aAtδ/Γ1+δ,A<0,(42)ν2x,t=2−Acoth−Ax+4aAtδ/Γ1+δ+2Acoth2−Ax+4aAtδ/Γ1+δ1−−Acoth−Ax+4aAtδ/Γ1+δ,A<0,(43)ν3x,t=−2AtanAx+4aAtδ/Γ1+δ−2Atan2Ax+4aAtδ/Γ1+δ1+AtanAx+4aAtδ/Γ1+δ,A>0,(44)ν4x,t=2AcotAx+4aAtδ/Γ1+δ−2Acot2Ax+4aAtδ/Γ1+δ1−AcotAx+4aAtδ/Γ1+δ,A>0.Setting the constraints of Set-2 into (40) along with results of auxiliary equation from Case-1 of our method, we attain to the following solutions:Case-1S2: ForB=0,C=1, the solutions are:
(45)ν5x,t=−Atanh−Ax+aAtδ/Γ1+δ+Atanh2−Ax+aAtδ/Γ1+δ1−−Atanh−Ax+aAtδ/Γ1+δ,A<0,ν6x,t=−Acoth−Ax+aAtδ/Γ1+δ+Acoth2−Ax+aAtδ/Γ1+δ1−−Acoth−Ax+aAtδ/Γ1+δ,A<0,ν7x,t=−AtanAx+aAtδ/Γ1+δ−Atan2Ax+aAtδ/Γ1+δ1+AtanAx+aAtδ/Γ1+δ,A>0,ν8x,t=AcotAx+aAtδ/Γ1+δ−Acot2Ax+aAtδ/Γ1+δ1−AcotAx+aAtδ/Γ1+δ,A>0.Case-2: ForA=0,B=−1,C=1, the system of equations leads to the solutions:Set-1:(46)p=1m122m0m1+m123am02−m0m1−6am02+12am0m1−7am12m12+8am0m1−m02m1+am12+10m0m1−8m02,l=3am0m1−m02+6am02−12am0m1+7am12m12,l0=m02−m0m1m1,l1=−m1,l2=−m1.Set-2:p=0,l=a,l0=3m1/2,l1=−2m1,l2=−2m1,m0=3m1/2.Set-3:p=−12a,l=13a,l0=m1/2,l1=−m1,l2=−m1,m0=−m1/2.Set-4:p=0,l=a,l0=m1/2,l1=−m1,l2=−m1,m0=m1/2.Set-5:(47)p=0,l=−a275−3/2±17/23+1815−3/2±17/22+2412−3/2±17/23−200422−3/2±17/2+78,l0=−m13−3/2±17/2+1111−3/2±17/2+39,l1=−32±172m1,l2=−2m1−3/2±17/2+3,m0=−m1/2.Setting the constraints of Set-1 to Set-5 into (40) along with the results of auxiliary equation from Case-2 yields the solutions:
(48)ν9x,t=m02−m0m11+heξ2−m122+heξm11+heξm01+heξ+m1,(49)ξ=x−3am0m1−m02+6am02−12am0m1+7am12tδm12Γ1+δ,(50)ν10x,t=−8−4hexpx−atδ/Γ1+δ+1+hexpx−atδ/Γ1+δ21+hexpx−atδ/Γ1+δ5+3hexpx−atδ/Γ1+δ,(51)ν11x,t=−4−2hexpx+12atδ/Γ1+δ+1+hexpx+12atδ/Γ1+δ21−h2exp2x+24atδ/Γ1+δ,(52)ν12x,t=−4−2hexpx−atδ/Γ1+δ+1+hexpx−atδ/Γ1+δ21+hexpx−atδ/Γ1+δ3+hexpx−atδ/Γ1+δ,(53)ν13x,t=−6L+22/11L+391+hexpx−ltδ/Γ1+δ2+2L1+hexpx−ltδ/Γ1+δ−4/L+31−1+hexpx−ltδ/Γ1+δ2,where l=−a275L3+1815L2+2412L−2004/22L+78 and L=−3/2±17/2.


## 4.5. The Multishock Waves to the FSTOM

In this section, we now explore a reliable technique [25] to achieve multishock wave of the fractional FSTOM. To do so, we consider the simple Burger equation (32) and Eq. (33). Then, we consider the wave transformation
(54)νx,t=νξ,ξ=lx−wtδΓ1+δ.Using Eq. (54), Eq. (4) leads to
(55)−wv′+3al2v′2+3alv2v′+3al2vv′′+al3v′′′=0.Integrating once with putting integral constant as zero, Eq. (55) provides
(56)–wv+3al2vν′+alν3+al3v′′=0.Let us consider the trial solution of Eq. (56) is
(57)v=a0+a1Y.asbalancenumberofthemodelisunityDifferentiating Eq. (57) together with Eq. (33) into Eq. (56), leads an algebraic system and solving, we reach to mutlishock wave solution.
(58)v=−∑j=1nlj+2∑j=1nljeljx−wtδ/Γ1+δ1+∑j=1neljx−wtδ/Γ1+δ,where
(59)wj=∑j=1nalj3.


## 5. Numerical Results and Discussions

Numerical illustrations of wave profiles are essential keys for exploration, analysis, and as a tool to communicate the achieved results of the problems coherently. In order to demonstrate the effect of waves in the corresponding application fields, we require definitive knowledge of securing the use of the properties in daily life. Hence, in this section, we physically explain and discuss the features of the solutions of the sFBM and tFSTOM via the MKM in 3D graphics. We discuss the impact of fractionality and parameters on the wave solutions and observe the potential changes of wave properties.


### 5.1. The Physical Descriptions for the Solutions of the sFBM via the MKM

We present here dynamical properties such as topological kink-dark bell and singular kink, bright bell solitons, and periodic waves of the sFBM. We obtain twelve solutions of the sFBM, among them the solutionsν1x,t and ν7x,t exhibit solitonic characters while the solutions ν2x,t and ν8x,t exhibit similar solitonic nature with singularities. Figure 1 represents topological kink-dark bell waves of ν1x,t. It is signified that the effect of order of differential fractionality is observed on the wave surface which gives the real situation as really most of the natural waves proceed in a curvy pattern. In Figure 1, one can observe that the wave surface is curvier for the lower fractional order (see Figure 1(a) and corresponding contour plot (d)), and curvy pattern gradually reduces to straight line surface as the order of fractionality increases (see Figure 1(b) and corresponding contour plot (e)), and ultimately being straight surface for complete order of the derivative, i.e., δ=1 (see Figure 1(c) and corresponding contour plot (f)). The solution ν3x,t comes from the rational solutions as a topological singular kink soliton (see Figure 2). However, the solutions ν4x,t and ν9x,t are periodic waves while the solutions ν5x,t and ν10x,t are similar periodic waves with singularities. Figure 3 represents periodic waves of ν4x,t. Besides, the solutions ν6x,t,ν11x,t, and ν12x,t come in terms of exponential functions that exhibit bright bell solitary waves as depicted in Figure 4 of ν11x,t. The impact of fractionality is similarly observed in all wave solutions like explanation of Figure 1.Figure 1
Topological kink-dark bell wave form1=1,m0=2, A=−1,l=1.5,k=q=1 of (23): (a, b) curvy surface for fractional order, (c) straight surface for complete order, (d, e, f) corresponding contour plots of beyond.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1(d)
δ=0.25(e)
δ=0.25(f)
δ=1Figure 2
Topological singular kink wave forl=k=q=m0=m1=1 of (25): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=0.5(b)
δ=1Figure 3
Periodic wave solution form1=1,m0=2, A=1,l=1.5,k=q=1 of (26): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=0.5(b)
δ=1Figure 4
Bright bell wave forl=k=q=h=m0=m1=1 of (30): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=3/4(b)
δ=1


### 5.2. The Physical Descriptions for the Multishock Wave Solutions of the sFBM

Shock waves and its multistages jumping process (multitimes shocks) are very important in the fields of shock waves in a viscous medium, magneto-hydrodynamic waves in a medium with finite electrical conductivity, turbulence, and in plasma dynamics. We present here overtaking collision of multishock waves to the sFBM and its different scattering way, distances from each other, amplitudes, and height of each shock are showed for different changes of parametera in Figures 5–7. The effect of dispersive parameters k, the coefficient of fractional derivative q, and fractional order on overtaking collisions of multishocks are illustrated graphically.Figure 5
Single shock wave forl1=1: (a) q=1,δ=0.9, (b) k=3,δ=0.9, K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)Figure 6
Two shock waves forl1=1,l2=0.5: (a) q=1,δ=0.9, (b) k=3,δ=0.9K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)Figure 7
Three shock waves forl1=1,l2=0.5,l3=1.5: (a) q=1,δ=0.9, (b) k=3,δ=0.9, K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)


### 5.3. The Physical Descriptions for the Solutions of the tFSTOM via the MKM

We present here dynamical properties such as topological kink and singular kink, bright and bark bell solitons, periodic waves, and peakon soliton of the tFSTOM. We obtain thirteen solutions of the tFSTOM, and among them, the solutionsν1x,t,ν5x,t exhibit solitonic character while ν2x,t,ν6x,t represent similar solitonic nature with singularities. Figure 8 represents topological kink waves specified by ν1x,t. It is signified that the effect of order of differential fractionality is observed on the wave surface which gives the real situation as really most of the natural waves proceed in a curvy pattern. In Figure 8, one can observe that the wave surface is curvy for the lower fractional order (see Figure 8(a) and corresponding contour plot (d)) and curvy pattern gradually reduces to straight line surface as the order of fractionality increases (see Figure 8(b) and corresponding contour plot (e)), and ultimately being straight surface for complete order of the derivative, i.e., δ=1 (see Figure 8(c) and corresponding contour plot (f)). The solutions ν3x,t and ν7x,t are periodic waves while the solutions ν4x,t,ν8x,t represent similar periodic waves with singularities. Figure 9 represents periodic waves of ν3x,t. The other solutions of ν9x,t,ν10x,t,ν11x,t,ν12x,t, and ν13x,t come in an exponential fashion and represent different solitonic characters for dissimilar values of the existing parameters. Figures 10–12, demonstrate topological kink, singular kink, and dark bell wave solutions by ν9x,t, ν11x,t, and ν13x,t, respectively. The nature of the solution ν10x,t,ν12x,t is similar to the result ν9x,t and the nature of the solution deformed into a single peak wave when the fractionality tend to unity (see Figures 10(b) and 10(c)). Besides this, Figure 11 specified by ν11x,t exhibits multipeakon soliton only for lower fractionality δ⟶0<0.25 (see Figure 11(a)), but it reduces to singular kink wave with multipeak as increases of the fractionality δ⟶1;δ∈0,1 (see Figures 11(b) and 11(c)). Figure 12 represents dark bell wave for the parametric value L=−3/2+17/2, but it is just being a single peak wave depicted by Figure 13 for L=−3/2−17/2. It is observed for Figure 13 that at the time of interactions or nearer, the waves’ height decreases exponentially with utmost deeper peak. The impact of fractionality is similarly observed in all wave solutions like explanation of Figure 8.Figure 8
Topological kink wave forA=−1,a=0.5 of (41): (a, b) curvy surface for fractional order, (c) straight surface for complete order, and (d, e, f) corresponding contour plots of (a, b, c).
(a)
δ=0.25(b)
δ=0.5(c)
δ=1(d)
δ=0.25(e)
δ=0.5(f)
δ=1Figure 9
Periodic wave solution forA=a=1 of (43): (a) for δ=0.25, (b) for δ=0.5, (c) for δ=1: (a, b) curvy surface for fractional order, (c) straight surface for complete order.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1Figure 10
Kink wave forA=a=m1=1,δ=0.5 of (48): (a) m0=1, (b) m0=2.5, (c) m0=3: (a, b) curvy surface for fractional order, (c) straight surface for complete order and deform into a single peakon wave.
(a)
m0=1(b)
m0=2.5(c)
m0=3Figure 11
Multipeakon soliton fora=h=1 of (51): (a) multipeakon soliton, (b) singular kink with curvy surface for fractional order, (c) singular kink with multipeakon straight surface for complete order of fractionality.
(a)
δ=0.25(b)
δ=0.75(c)
δ=1Figure 12
Dark bell wave fora=h=1 of (53): (a, b) curvy surface for fractional order, (c) straight surface for complete order of fractionality.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1Figure 13
Peakon wave fora=h=1 of (53): L=−3/2−17/2: (a) curvy surface for fractional order, (b) comparatively straight surface for fractionality δ=0.5.
(a)
δ=0.25(b)
δ=0.5


### 5.4. The Physical Descriptions for the Multishock Wave Solutions of the tFSTOM

We present here overtaking collision of multishock waves to the tFSTOM and its different scattering way, amplitudes, and height of each shock are showed for different changes of parametera in Figure 14. From the view, one interested researcher and engineer might realize the effect of parameters, which parameter, and how the value can proceed quickly to overtake the other waves. Through the images, we can effectively analyze the effect of fractional power and parameters on the wave solutions, and observe the potential changes of characteristics of waves.Figure 14
(a) Single shock wave forδ=0.9,l1=0.5 at t=40. (b) Two shock waves for δ=0.9,l1=0.5,l2=1 at t=10. (c) Five shock waves for δ=0.9,l1=0.5,l2=1,l3=1.5,l4=2,l5=2.5 at x=40.
(a)(b)(c)


## 5.1. The Physical Descriptions for the Solutions of the sFBM via the MKM

We present here dynamical properties such as topological kink-dark bell and singular kink, bright bell solitons, and periodic waves of the sFBM. We obtain twelve solutions of the sFBM, among them the solutionsν1x,t and ν7x,t exhibit solitonic characters while the solutions ν2x,t and ν8x,t exhibit similar solitonic nature with singularities. Figure 1 represents topological kink-dark bell waves of ν1x,t. It is signified that the effect of order of differential fractionality is observed on the wave surface which gives the real situation as really most of the natural waves proceed in a curvy pattern. In Figure 1, one can observe that the wave surface is curvier for the lower fractional order (see Figure 1(a) and corresponding contour plot (d)), and curvy pattern gradually reduces to straight line surface as the order of fractionality increases (see Figure 1(b) and corresponding contour plot (e)), and ultimately being straight surface for complete order of the derivative, i.e., δ=1 (see Figure 1(c) and corresponding contour plot (f)). The solution ν3x,t comes from the rational solutions as a topological singular kink soliton (see Figure 2). However, the solutions ν4x,t and ν9x,t are periodic waves while the solutions ν5x,t and ν10x,t are similar periodic waves with singularities. Figure 3 represents periodic waves of ν4x,t. Besides, the solutions ν6x,t,ν11x,t, and ν12x,t come in terms of exponential functions that exhibit bright bell solitary waves as depicted in Figure 4 of ν11x,t. The impact of fractionality is similarly observed in all wave solutions like explanation of Figure 1.Figure 1
Topological kink-dark bell wave form1=1,m0=2, A=−1,l=1.5,k=q=1 of (23): (a, b) curvy surface for fractional order, (c) straight surface for complete order, (d, e, f) corresponding contour plots of beyond.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1(d)
δ=0.25(e)
δ=0.25(f)
δ=1Figure 2
Topological singular kink wave forl=k=q=m0=m1=1 of (25): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=0.5(b)
δ=1Figure 3
Periodic wave solution form1=1,m0=2, A=1,l=1.5,k=q=1 of (26): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=0.5(b)
δ=1Figure 4
Bright bell wave forl=k=q=h=m0=m1=1 of (30): (a) curvy surface for fractional order, (b) straight surface for complete order.
(a)
δ=3/4(b)
δ=1


## 5.2. The Physical Descriptions for the Multishock Wave Solutions of the sFBM

Shock waves and its multistages jumping process (multitimes shocks) are very important in the fields of shock waves in a viscous medium, magneto-hydrodynamic waves in a medium with finite electrical conductivity, turbulence, and in plasma dynamics. We present here overtaking collision of multishock waves to the sFBM and its different scattering way, distances from each other, amplitudes, and height of each shock are showed for different changes of parametera in Figures 5–7. The effect of dispersive parameters k, the coefficient of fractional derivative q, and fractional order on overtaking collisions of multishocks are illustrated graphically.Figure 5
Single shock wave forl1=1: (a) q=1,δ=0.9, (b) k=3,δ=0.9, K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)Figure 6
Two shock waves forl1=1,l2=0.5: (a) q=1,δ=0.9, (b) k=3,δ=0.9K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)Figure 7
Three shock waves forl1=1,l2=0.5,l3=1.5: (a) q=1,δ=0.9, (b) k=3,δ=0.9, K=3, q=1, and (c) q=1,k=3 at x=40.
(a)(b)(c)


## 5.3. The Physical Descriptions for the Solutions of the tFSTOM via the MKM

We present here dynamical properties such as topological kink and singular kink, bright and bark bell solitons, periodic waves, and peakon soliton of the tFSTOM. We obtain thirteen solutions of the tFSTOM, and among them, the solutionsν1x,t,ν5x,t exhibit solitonic character while ν2x,t,ν6x,t represent similar solitonic nature with singularities. Figure 8 represents topological kink waves specified by ν1x,t. It is signified that the effect of order of differential fractionality is observed on the wave surface which gives the real situation as really most of the natural waves proceed in a curvy pattern. In Figure 8, one can observe that the wave surface is curvy for the lower fractional order (see Figure 8(a) and corresponding contour plot (d)) and curvy pattern gradually reduces to straight line surface as the order of fractionality increases (see Figure 8(b) and corresponding contour plot (e)), and ultimately being straight surface for complete order of the derivative, i.e., δ=1 (see Figure 8(c) and corresponding contour plot (f)). The solutions ν3x,t and ν7x,t are periodic waves while the solutions ν4x,t,ν8x,t represent similar periodic waves with singularities. Figure 9 represents periodic waves of ν3x,t. The other solutions of ν9x,t,ν10x,t,ν11x,t,ν12x,t, and ν13x,t come in an exponential fashion and represent different solitonic characters for dissimilar values of the existing parameters. Figures 10–12, demonstrate topological kink, singular kink, and dark bell wave solutions by ν9x,t, ν11x,t, and ν13x,t, respectively. The nature of the solution ν10x,t,ν12x,t is similar to the result ν9x,t and the nature of the solution deformed into a single peak wave when the fractionality tend to unity (see Figures 10(b) and 10(c)). Besides this, Figure 11 specified by ν11x,t exhibits multipeakon soliton only for lower fractionality δ⟶0<0.25 (see Figure 11(a)), but it reduces to singular kink wave with multipeak as increases of the fractionality δ⟶1;δ∈0,1 (see Figures 11(b) and 11(c)). Figure 12 represents dark bell wave for the parametric value L=−3/2+17/2, but it is just being a single peak wave depicted by Figure 13 for L=−3/2−17/2. It is observed for Figure 13 that at the time of interactions or nearer, the waves’ height decreases exponentially with utmost deeper peak. The impact of fractionality is similarly observed in all wave solutions like explanation of Figure 8.Figure 8
Topological kink wave forA=−1,a=0.5 of (41): (a, b) curvy surface for fractional order, (c) straight surface for complete order, and (d, e, f) corresponding contour plots of (a, b, c).
(a)
δ=0.25(b)
δ=0.5(c)
δ=1(d)
δ=0.25(e)
δ=0.5(f)
δ=1Figure 9
Periodic wave solution forA=a=1 of (43): (a) for δ=0.25, (b) for δ=0.5, (c) for δ=1: (a, b) curvy surface for fractional order, (c) straight surface for complete order.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1Figure 10
Kink wave forA=a=m1=1,δ=0.5 of (48): (a) m0=1, (b) m0=2.5, (c) m0=3: (a, b) curvy surface for fractional order, (c) straight surface for complete order and deform into a single peakon wave.
(a)
m0=1(b)
m0=2.5(c)
m0=3Figure 11
Multipeakon soliton fora=h=1 of (51): (a) multipeakon soliton, (b) singular kink with curvy surface for fractional order, (c) singular kink with multipeakon straight surface for complete order of fractionality.
(a)
δ=0.25(b)
δ=0.75(c)
δ=1Figure 12
Dark bell wave fora=h=1 of (53): (a, b) curvy surface for fractional order, (c) straight surface for complete order of fractionality.
(a)
δ=0.25(b)
δ=0.5(c)
δ=1Figure 13
Peakon wave fora=h=1 of (53): L=−3/2−17/2: (a) curvy surface for fractional order, (b) comparatively straight surface for fractionality δ=0.5.
(a)
δ=0.25(b)
δ=0.5


## 5.4. The Physical Descriptions for the Multishock Wave Solutions of the tFSTOM

We present here overtaking collision of multishock waves to the tFSTOM and its different scattering way, amplitudes, and height of each shock are showed for different changes of parametera in Figure 14. From the view, one interested researcher and engineer might realize the effect of parameters, which parameter, and how the value can proceed quickly to overtake the other waves. Through the images, we can effectively analyze the effect of fractional power and parameters on the wave solutions, and observe the potential changes of characteristics of waves.Figure 14
(a) Single shock wave forδ=0.9,l1=0.5 at t=40. (b) Two shock waves for δ=0.9,l1=0.5,l2=1 at t=10. (c) Five shock waves for δ=0.9,l1=0.5,l2=1,l3=1.5,l4=2,l5=2.5 at x=40.
(a)(b)(c)


## 6. Conclusions

We have successfully achieved the solitary wave solutions of a class of the conformable sFB and the tFSTO models via the proposed resourceful integral MKM. Periodic, topological kink, singular kink, peakon soliton, and periodic wave soliton solutions are derived by the straight forward approach. The retrieved explicit solutions are presented in a more general form of the fractional differential models. The impact of fractionality on the wave shape and its deformation is analyzed and discussed graphically. Beside this, we achieved multishock wave solutions to the both fractional models and showed the parametric effect on the changes of shock waves. To visualize the real properties of the solutions, the graphical elucidation in 3D and 2D profiles of the gained solutions are presented for suitable parametric selection. Moreover, the generalized Kudryashov [12, 13] and the extended Kudryashov methods [14] are the special cases of this proposed MKM. From our computational effort and achievement, we emphasize that the proposed scheme is very simple, highly effective, and a powerful mathematical tool to extract exact solitary wave solutions of both differential and fractional differential models. The achieved solutions are realistic and appropriate to use and they penetrate the fission and fusion occurrences in solitons, fluid dynamics, sound waves in a viscous medium, gas dynamics, and plasma dynamics.



---

*Source: 7096486-2022-05-23.xml*