Next Article in Journal
Pion Productions with Isospin-Dependent In-Medium Cross Sections
Next Article in Special Issue
Two Analytical Schemes for the Optical Soliton Solution of the (2 + 1) Hirota–Maccari System Observed in Single-Mode Fibers
Previous Article in Journal
Ionospheric TEC Forecasting over an Indian Low Latitude Location Using Long Short-Term Memory (LSTM) Deep Learning Network
Previous Article in Special Issue
Quiescent Optical Solitons with Kudryashov’s Generalized Quintuple-Power and Nonlocal Nonlinearity Having Nonlinear Chromatic Dispersion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Soliton Solutions of Generalized Fifth Order Time-Fractional KdV Models: Laplace Transform with Homotopy Perturbation Algorithm

1
Department of Sciences and Humanities, National University of Computer and Emerging Sciences, Lahore 54770, Pakistan
2
Faculty of Technical Physics, Information Technology and Applied Mathematics, Lodz University of Technology, 90-924 Lodz, Poland
3
Faculty of Mechanical Engineering, Lodz University of Technology, 90-924 Lodz, Poland
*
Author to whom correspondence should be addressed.
Universe 2022, 8(11), 563; https://doi.org/10.3390/universe8110563
Submission received: 28 September 2022 / Revised: 22 October 2022 / Accepted: 25 October 2022 / Published: 27 October 2022
(This article belongs to the Special Issue Research on Optical Soliton Perturbation)

Abstract

:
The main purpose of this research is to propose a new methodology to observe a class of time-fractional generalized fifth-order Korteweg–de Vries equations. Laplace transform along with a homotopy perturbation algorithm is utilized for the solution and analysis purpose in the current study. This extended technique provides improved and convergent series solutions through symbolic computation. The proposed methodology is applied to time-fractional Sawada–Kotera, Ito, Lax’s, and Kaup–Kupershmidt models, which are induced from a generalized fifth-order KdV equation. For validity purposes, obtained and existing results at integral orders are compared. Convergence analysis was also performed by computing solutions and errors at different values in a fractional domain. Dynamic behavior of the fractional parameter is also studied graphically. Simulations affirm the dominance of the proposed algorithm in terms of accuracy and fewer computations as compared to other available schemes for fractional KdVs. Hence, the projected algorithm can be utilized for more advanced fractional models in physics and engineering.

1. Introduction

Modeling and analysis of physical phenomena is an essential step in scientific work. There are many important mathematical models for capturing different situations in various fields. Among many, Korteweg–de Vries models are essential due to their ability of capturing physical situations, such as the motion of long waves in shallow water with weak non-linearities. Initially, this model was derived as 1D small amplitude and long surface gravity waves. Flood analysis and ocean flow analysis are a few uses of shallow water equations. Scott Russell in 1834 [1], Boussinesq and Rayleigh in 1870 [2,3,4], and Korteweg and De Vries in 1895 [5] played a major role in the discovery of the KdV model. These models have many uses in collision-free hydromagnetic waves, stratified internal waves, shock wave formation [6], fluid and quantum mechanics [7], and in biology [8].
A generalized fifth-order time-fractional KdV equation is usually in the following form of non-linear partial differential equation(PDE):
β J t β + a J 2 J x + b J x 2 J x 2 + c J 3 J x 3 + 5 J x 5 = 0 ,
where J is the wave function with x and t as the space and time variables, respectively. In addition, β represents the fractional order parameter. The model in Equation (1) has various physical applications in acoustic magnetic propagation in plasma, incompressible and inviscid fluids, gravitational field, etc. It typically consists of one linear dispersive term and three non-linear terms. The linear dispersive term has a momentous role in the balancing of non-linearity and dispersion effects of soliton behavior [9]; a, b, and c in Equation (1) are real constant parameters. Different values of these parameters give different versions of time-fractional KdV equations. For instance, a = 45, b = 15, and c = 15 give a time-fractional, Sawada–Kotera equation [10]. Fixing a, b, and c as 2, 6, and 3, respectively, leads to a time-fractional Ito equation [11]. In addition, when a = 30, b = 20, and c = 10, it gives a time-fractional Lax’s KdV equation [12,13]. Similarly, a = 45, b = −15, and c = −15 leads to a time-fractional Kaup–Kupershmidt equation [14]. Many researchers tried to find exact solutions of these equations in [15,16,17,18].
Since fractional calculus deals with the calculation of non-integer order derivatives and integrals, it appeals to a wide range of audiences while dealing with fractional ordinary and partial differential equations. Different numerical and analytical methods, including the homotopy analysis method [19], the residual power series method [20], the variation iteration method [21], the Adomian decomposition method [22], the Bernstein polynomials method [23], homotopy perturbation method [24], etc., are used in the literature while dealing with fractional problems. Many modifications of existing techniques are also applied to fractional KdV equations [25,26,27].
The homotopy perturbation method (HPM), introduced by He [28], provides a convenient technique to find analytical solutions of linear and non-linear differential equations in both fractional and integer form. To reduce errors and increase its reliability, many modifications of HPM [29,30,31,32] have been adapted in the literature. The Laplace homotopy perturbation method (LHPM), which combines HPM with Laplace transform, is an efficient and convenient modification of HPM for solving fractional and integral differential equations [33]. Results obtained from LHPM are highly accurate and straightforward without imposing any restrictions on the concerned model. In this paper, we have extended LHPM to generalized fifth-order highly non-linear time-fractional KdV models. In the rest of the manuscript, Section 2 and Section 3 consist of preliminaries and the general concept of LHPM for fractional KdV equations. Convergence and error estimation are presented in Section 4. Applications and solutions of different KdV models are explained in Section 5. Discussion and conclusion are presented in Section 6 and Section 7, respectively.

2. Preliminaries

Definition 1
([34]). The Laplace transform L of the Riemann–Liouville time-fractional integral I t β on a function J ( v , t ) is described as:
L [ I t β J ( v , t ) ] = s β L [ J ( v , t ) ] , f 1 < β f .
Definition 2
([35]). The Laplace transform L of Caputo’s time-fractional derivative D t β on a function J ( v , t ) is described as:
L [ D t β J ( v , t ) ] = s β L [ J ( v , t ) ] j = 0 f 1 s β j 1 J ( f ) ( v , 0 ) , f 1 < β f .
where J ( f ) ( v , 0 ) represents the initial conditions.
Lemma 1
([36]). Let J be a function, then for a positive constant c, the stability result states
J ( u ) J ( v ) c u v , u , v ϵ R n .

3. Fundamental Concept of the Laplace Homotopy Perturbation Method for Fifth-Order Time-Fractional KdV Models

Let us consider a general non-linear, fifth-order, time-fractional KdV equation as:
β J ( x , t ) t β + a J 2 ( x , t ) J ( x , t ) x + b J ( x , t ) x 2 J ( x , t ) x 2 + c J ( x , t ) 3 J ( x , t ) x 3 + 5 J ( x , t ) x 5 = 0 , x ϵ Ω , t > 0 , f 1 < β f ,
with initial conditions
J ( f ) ( x , 0 ) = I f , f = 0 , 1 , 2 , ,
where J ( x , t ) is an unknown function with β t β as its fractional derivative; a, b, and c are constant parameters that give distant versions of time fractional KdV equations for different values.
The first step of the Laplace transform algorithm is applied on both sides of (5), which gives
L { β J ( x , t ) t β + a J 2 ( x , t ) J ( x , t ) x + b J ( x , t ) x 2 J ( x , t ) x 2 + c J ( x , t ) 3 J ( x , t ) x 3 + 5 J ( x , t ) x 5 } = 0 ,
Using (3) gives
L [ J ( x , t ) ] 1 s β j = 0 f 1 s β j 1 J ( f ) ( x , 0 ) + 1 s β L { a J 2 ( x , t ) J ( x , t ) x + b J ( x , t ) x 2 J ( x , t ) x 2 + c J ( x , t ) 3 J ( x , t ) x 3 + 5 J ( x , t ) x 5 } = 0 ,
Now, we can construct a homotopy
H = ( 1 p ) ( L { J ( x , t ; p ) } J 0 ( x , t ) ) + p ( L { J ( x , t ; p ) } 1 s β j = 0 f 1 s β j 1 J ( f ) ( x , 0 ) + 1 s β L a J 2 ( x , t ) J ( x , t ) x + b J ( x , t ) x 2 J ( x , t ) x 2 + c J ( x , t ) 3 J ( x , t ) x 3 + 5 J ( x , t ) x 5 ) ,
where J 0 ( x , t ) is the initial guess, and it satisfies the given conditions.
Expanding J ( x , t ) in Taylor series with regard to p gives:
J ( x , t ; p ) = m = 1 p m J m ,
Substituting (10) into (9) and then comparing the coefficients of p, we have distinct order problems.
The problem at the first order is:
L { J 1 ( x , t ) } + J 0 ( x , t ) 1 s β j = 0 f 1 s β j 1 J ( f ) ( x , 0 ) + 1 s β L { a J 0 2 ( x , t ) J 0 ( x , t ) x + b J 0 ( x , t ) x 2 J 0 ( x , t ) x 2 + c J 0 ( x , t ) 3 J 0 ( x , t ) x 3 + 5 J 0 ( x , t ) x 5 } = 0 ,
Employing inverse Laplace transform gives
J 1 ( x , t ) = L 1 J 0 ( x , t ) + 1 s β j = 0 f 1 s β j 1 J ( f ) ( x , 0 ) L 1 { 1 s β L { a J 0 2 ( x , t ) J 0 ( x , t ) x + b J 0 ( x , t ) x 2 J 0 ( x , t ) x 2 + c J 0 ( x , t ) 3 J 0 ( x , t ) x 3 + 5 J 0 ( x , t ) x 5 } } ,
The the problem at the kth order is:
L { J k ( x , t ) } + 1 s β L { a J k 1 2 ( x , t ) J k 1 ( x , t ) x + b J k 1 ( x , t ) x 2 J k 1 ( x , t ) x 2 + c J k 1 ( x , t ) 3 J k 1 ( x , t ) x 3 + 5 J k 1 ( x , t ) x 5 } = 0 ,
By taking Inverse Laplace transform, we have
J k ( x , t ) = L 1 { 1 s β L { a J k 1 2 ( x , t ) J k 1 ( x , t ) x + b J k 1 ( x , t ) x 2 J k 1 ( x , t ) x 2 + c J k 1 ( x , t ) 3 J k 1 ( x , t ) x 3 + 5 J k 1 ( x , t ) x 5 } } ,
The approximate solution of the general fifth-order KdV equation is
J ˜ = J 0 ( x , t ) + J 1 ( x , t ) + J 2 ( x , t ) + J 3 ( x , t ) + ,
Residual error is observed by substituting (15) in a given time fractional KdV Equation (5) as:
Res = β J ˜ t β + a J ˜ 2 J ˜ x + b J ˜ x 2 J ˜ x 2 + c J ˜ ( x , t ) 3 J ˜ x 3 + 5 J ˜ x 5 .

4. Convergence and Error Estimation of LHPM for Fractional KdV Equation

4.1. Convergence

Theorem 1.
Consider a Banach space ( B [0,T], . ) and suppose J n ( x , t ) and J ( x , t ) are defined in it. Then, for a constant ζ where 0 < ζ < 1, the series solution in (15) converges to the solution of fractional KdV (5).
Proof. 
Let { a n } be the sequence of partial sums of (15). We have to prove that a n ( x , t ) is a Cauchy sequence in ( B [0,T], . ). Consider
a n + 1 ( x , t ) a n ( x , t ) = J n + 1 ( x , t ) ζ J n ( x , t ) ζ 2 J n 1 ( x , t ) ζ n + 1 J 0 ( x , t ) ,
Now, for partial sums a n and a m where n,m ϵ N and n m , using triangle inequality, we obtain
a n a m = a n ( x , t ) a n 1 ( x , t ) + a n 1 ( x , t ) a n 2 ( x , t ) + + a m + 1 ( x , t ) a m ( x , t ) a n ( x , t ) a n 1 ( x , t ) + a n 1 ( x , t ) a n 2 ( x , t ) + + a m + 1 ( x , t ) a m ( x , t ) ,
Using (17), we have
a n a m ζ n J 0 ( x , t ) + ζ n 1 J 0 ( x , t ) + + ζ m + 1 J 0 ( x , t ) ( ζ n + ζ n 1 + + ζ m + 1 ) J 0 ( x , t ) ζ m + 1 ( ζ n m 1 + ζ n m 2 + + ζ + 1 ) J 0 ( x , t ) ζ m + 1 1 ζ n m 1 ζ J 0 ( x , t ) ,
Since 0 < ζ < 1, therefore, 1 ζ n m < 1. Thus, we have
a n a m ζ m + 1 1 ζ m a x | J 0 ( x , t ) | , t ϵ [ 0 , T ] ,
Since J 0 is bounded, so
lim n , m a n ( x , t ) a m ( x , t ) = 0 .
Hence, we have proved that a n ( x , t ) is a Cauchy sequence in Banach space ( B [0,T], . ). Thus, the series solution in (15) converges to the solution of (5). □

4.2. Error Estimation

Theorem 2.
Consider the time fractional KdV Equation (5), then the maximum absolute truncation error of its solution (15) is
| J ( x , t ) j = 0 m J j ( x , t ) | ζ m + 1 1 ζ J 0 ( x , t ) .
Proof. 
From (19) we have,
J ( x , t ) a m ζ m + 1 1 ζ n m 1 ζ J 0 ( x , t ) ,
Since 0 < ζ < 1, therefore 1 ζ n m < 1. Thus we have
| J ( x , t ) j = 0 m J j ( x , t ) | ζ m + 1 1 ζ J 0 ( x , t ) .

5. Solutions of Time-Fractional KdV Models Using Laplace Homotopy Perturbation Method

Example 1.
Consider the time fractional, Sawada–Kotera equation
β J t β = 45 J 2 J x 15 J x 2 J x 2 15 J 3 J x 3 5 J x 5 , 0 < β 1 , t > 0 ,
associated with initial condition
J ( x , 0 ) = 2 k 2 sec h 2 ( k x ) ,
where k 0. The exact solution of (25) is:
J ( x , t ) = 2 k 2 sec h 2 ( k ( x 16 k 4 t ) ) .
Solution: Taking Laplace transform of (25) and then by using (3), we have
s β L [ J ( x , t ) ] s β 1 2 k 2 sec h 2 ( k x ) L 45 J 2 J x 15 J x 2 J x 2 15 J 3 J x 3 5 J x 5 = 0 ,
Using (9), homotopy can be constructed as follows:
H = ( 1 p ) ( L { J ( x , t ) } J 0 ( x , t ) ) + p ( L { J ( x , t ) } 1 s 2 k 2 sec h 2 ( k x ) 1 s β L 45 J 2 J x 15 J x 2 J x 2 15 J 3 J x 3 5 J x 5 ) ,
where J 0 ( x , t ) is the initial guess which satisfies the initial condition (26). For the current problem, the following initial guess is taken
J 0 ( x , t ) = 2 k 2 sec h 2 ( k x ) ,
Expanding J ( x , t ) in Taylor’s series with regard to p, and then comparing the coefficients of identical powers of p, we have.
First-order problem:
L { J 1 ( x , t ) } + J 0 ( x , t ) 1 s 2 k 2 sec h 2 ( k x ) 1 s β L { 45 J 0 2 J 0 x 15 J 0 x 2 J 0 x 2 15 J 0 3 J 0 x 3 5 J 0 x 5 } = 0 , J 1 ( x , 0 ) = 0 ,
Taking inverse Laplace transform leads to a solution at the first order:
J 1 ( x , t ) = t k 7 ( 64 sec h 6 ( x k ) tanh ( x k ) + 128 sec h 4 ( x k ) tanh 3 ( x k ) + 64 sec h 2 ( x k ) tanh 5 ( x k ) ) ,
Second-order problem:
L { J 2 ( x , t ) } 1 s β L 45 J 1 2 J 1 x 15 J 1 x 2 J 1 x 2 15 J 1 3 J 1 x 3 5 J 1 x 5 = 0 , J 2 ( x , 0 ) = 0 ,
Second-order solution:
J 2 ( x , t ) = 512 t 2 k 12 sec h 2 ( x , t ) ( sec h 2 ( x , t ) 2 tanh 2 ( x , t ) ) ( sec h 2 ( x , t ) + tanh 2 ( x , t ) ) 4 ,
Third-order problem:
L { J 3 ( x , t ) } 1 s β L 45 J 2 2 J 2 x 15 J 2 x 2 J 2 x 2 15 J 2 3 J 2 x 3 5 J 2 x 5 = 0 , J 3 ( x , 0 ) = 0 ,
Third order-solution:
J 3 ( x , t ) = 32768 3 t 3 k 17 sec h 2 ( x , t ) tanh ( x , t ) ( 2 sec h 2 ( x , t ) tanh 2 ( x , t ) ) ( sec h 2 ( x , t ) + tanh 2 ( x , t ) ) 6 ,
Fourth-order problem:
L { J 4 ( x , t ) } 1 s β L 45 J 3 2 J 3 x 15 J 3 x 2 J 3 x 2 15 J 3 3 J 3 x 3 5 J 3 x 5 = 0 , J 4 ( x , 0 ) = 0 ,
Fourth-order solution:
J 4 ( x , t ) = 131072 3 t 4 k 22 sec h 2 ( x , t ) ( sec h 2 ( x , t ) + tanh 2 ( x , t ) ) 8 ( 2 sec h 4 ( x , t ) 11 sec h 2 ( x , t ) tanh 2 ( x , t ) + 2 tanh 4 ( x , t ) ) ,
The rest of values for J i ( x , t ) with i ≥ 5 can be computed in a similar way. The approximate solution of (25) can be captured by
J ˜ = i = 0 4 J i ( x , t )
Residual error is obtained by substituting (39) in (25) as follows:
Res J = β J ˜ t β + 45 J ˜ 2 J ˜ x + 15 J ˜ x 2 J ˜ x 2 + 15 J ˜ 3 J ˜ x 3 + 5 J ˜ x 5 .
Example 2.
Consider the following time-fractional Ito model
β J t β = 2 J 2 J x 6 J x 2 J x 2 3 J 3 J x 3 5 J x 5 , 0 < β 1 , t > 0 ,
along with the initial condition
J ( x , 0 ) = 20 k 2 30 k 2 tanh 2 ( k x ) ,
where k 0. The exact solution of (41) is
J ( x , t ) = 20 k 2 30 k 2 tanh 2 ( k ( x 96 k 4 t ) ) .
Solution: Applying the Laplace transform on (41) and using (3) give the following:
s β L [ J ( x , t ) ] s β 1 20 k 2 30 k 2 tanh 2 ( k x ) L 2 J 2 J x 6 J x 2 J x 2 3 J 3 J x 3 5 J x 5 = 0 ,
Utilizing (9) homotopy for the given problem is
H = ( 1 p ) ( L { J ( x , t ) } J 0 ( x , t ) ) + p ( L { J ( x , t ) } 1 s 20 k 2 30 k 2 tanh 2 ( k x ) 1 s β L 2 J 2 J x 6 J x 2 J x 2 3 J 3 J x 3 5 J x 5 ) ,
where
J 0 ( x , t ) = 20 k 2 30 k 2 tanh 2 ( k x ) ,
By replacing (10) in (41) and then comparing coefficients of like power of p, leads to the following:
First-order problem:
L { J 1 ( x , t ) } + J 0 ( x , t ) 1 s 20 k 2 30 k 2 tanh 2 ( k x ) 1 s β L { 2 J 0 2 J 0 x 6 J 0 x 2 J 0 x 2 3 J 0 3 J 0 x 3 5 J 0 x 5 } = 0 , J 1 ( x , 0 ) = 0 ,
Second-order problem:
L { J 2 ( x , t ) } 1 s β L 2 J 1 2 J 1 x 6 J 1 x 2 J 1 x 2 3 J 1 3 J 1 x 3 5 J 1 x 5 = 0 , J 2 ( x , 0 ) = 0 ,
Third-order problem:
L { J 3 ( x , t ) } 1 s β L 2 J 2 2 J 2 x 6 J 2 x 2 J 2 x 2 3 J 2 3 J 2 x 3 5 J 2 x 5 = 0 , J 3 ( x , 0 ) = 0 ,
Fourth-order problem:
L { J 4 ( x , t ) } 1 s β L 2 J 3 2 J 3 x 6 J 3 x 2 J 3 x 2 3 J 3 3 J 3 x 3 5 J 3 x 5 = 0 , J 4 ( x , 0 ) = 0 ,
Continuing this way, higher order problems and solutions can be formulated. The approximate solution of (41) is
J ˜ = J 0 ( x , t ) + J 1 ( x , t ) + J 2 ( x , t ) + J 3 ( x , t ) + J 4 ( x , t ) +
Residual error is derived by substituting the obtained approximate solution in (41)
Res J = β J ˜ t β + 2 J ˜ 2 J ˜ x + 6 J ˜ x 2 J ˜ x 2 + 3 J ˜ 3 J ˜ x 3 + 5 J ˜ x 5 .
Example 3.
Consider the following time-fractional Lax’s KdV model
β J t β = 30 J 2 J x 20 J x 2 J x 2 10 J 3 J x 3 5 J x 5 , 0 < β 1 , t > 0 ,
with initial condition
J ( x , 0 ) = 2 k 2 ( 2 3 tanh 2 ( k x ) ) ,
where k 0 is an arbitrary constant. The exact solution of (53) is
J ( x , t ) = 2 k 2 ( 2 3 tanh 2 ( k ( x 56 k 4 t ) ) ) .
Solution:
Following the similar steps mapped out in Section 3, we have
First-order problem:
L { J 1 ( x , t ) } + J 0 ( x , t ) 1 s 2 k 2 ( 2 3 tanh 2 ( k x ) ) 1 s β L { 30 J 0 2 J 0 x 20 J 0 x 2 J 0 x 2 10 J 0 3 J 0 x 3 5 J 0 x 5 } = 0 , J 1 ( x , 0 ) = 0 ,
Second-order problem:
L { J 2 ( x , t ) } 1 s β L 30 J 1 2 J 1 x 20 J 1 x 2 J 1 x 2 10 J 1 3 J 1 x 3 5 J 1 x 5 = 0 , J 2 ( x , 0 ) = 0 ,
Third-order problem:
L { J 3 ( x , t ) } 1 s β L 30 J 2 2 J 2 x 20 J 2 x 2 J 2 x 2 10 J 2 3 J 2 x 3 5 J 2 x 5 = 0 , J 3 ( x , 0 ) = 0 ,
Fourth-order problem:
L { J 4 ( x , t ) } 1 s β L 30 J 3 2 J 3 x 20 J 3 x 2 J 3 x 2 10 J 3 3 J 3 x 3 5 J 3 x 5 = 0 , J 4 ( x , 0 ) = 0 ,
The inverse Laplace transform leads to the approximate solution J ˜ ( x , t ) . The residual error of the current problem is
Res J = β J ˜ t β + 30 J ˜ 2 J ˜ x + 20 J ˜ x 2 J ˜ x 2 + 10 J ˜ 3 J ˜ x 3 + 5 J ˜ x 5 .
Results related to Example 3 are in Table 1 and Table 2 and Figure 1, Figure 2 and Figure 3.
Example 4.
Consider the following time-fractional Kaup–Kupershmidt model
β J t β = 45 J 2 J x + 15 ρ J x 2 J x 2 + 15 J 3 J x 3 5 J x 5 , 0 < β 1 , t > 0 ,
associated with initial condition
J ( x , 0 ) = 1 4 w 2 λ 2 sec h 2 ( w λ x 2 ) + w 2 λ 2 12 ,
where λ and w 0 are arbitrary constants. The exact solution of (61) is:
J ( x , t ) = 1 4 w 2 λ 2 sec h 2 ( w 5 λ t β ( 8 λ 2 μ + 16 μ 2 + λ 4 ) 32 β + w λ x 2 ) + w 2 λ 2 12 ,
Solution:
Using basic steps given in Section 3, we obtain the following
First-order problem:
L { J 1 ( x , t ) } + J 0 ( x , t ) 1 s 1 4 w 2 λ 2 sec h 2 ( w λ x 2 ) + w 2 λ 2 12 1 s β L { 45 J 0 2 J 0 x + 15 ρ J 0 x 2 J 0 x 2 + 15 J 0 3 J 0 x 3 5 J 0 x 5 } = 0 , J 1 ( x , 0 ) = 0 ,
Second-order problem:
L { J 2 ( x , t ) } 1 s β L 45 J 1 2 J 1 x + 15 ρ J 1 x 2 J 1 x 2 + 15 J 1 3 J 1 x 3 5 J 1 x 5 = 0 , J 2 ( x , 0 ) = 0 ,
Third-order problem:
L { J 3 ( x , t ) } 1 s β L 45 J 2 2 J 2 x + 15 ρ J 2 x 2 J 2 x 2 + 15 J 2 3 J 2 x 3 5 J 2 x 5 = 0 , J 3 ( x , 0 ) = 0 ,
Fourth-order problem:
L { J 4 ( x , t ) } 1 s β L 45 J 3 2 J 3 x + 15 ρ J 3 x 2 J 3 x 2 + 15 J 3 3 J 3 x 3 5 J 3 x 5 = 0 , J 4 ( x , 0 ) = 0 ,
The inverse Laplace transform leads to the approximate solution J ˜ ( x , t ) . The residual error of the concerned problem is
Res J = β J ˜ t β + 45 J ˜ 2 J ˜ x 15 ρ J ˜ x 2 J ˜ x 2 15 J ˜ 3 J ˜ x 3 + 5 J ˜ x 5 .

6. Discussion

In this manuscript, the Laplace transform with homotopy perturbation method is proposed for the solution of generalized fifth-order time-fractional KdV models. The proposed method is tested against different time-fractional KdV models, including Sawada–Kotera, Ito, Lax, and Kaup–Kupershmidt KdV equations. These equations belong to the prominent form of fractional KdV family hierarchy. The main focus of this paper is to propose a new method for the solution and analysis of KdV equations in a fractional environment. Graphical and numerical comparisons of the obtained results were also made to provide strong evidence for the usage of the proposed technique. Obtained results are compared with existing ones from the literature.
The proposed method is firstly applied to the Sawada–Kotera model in Example 1, and the results are shown in Table 3 and Table 4. Table 3 presents a comparison of LHPM with exact and fractional conformable residual power series algorithm (FCRPSA) results. These results showed the efficiency of LHPM over FCRPSA. In Table 4, error analysis has been performed in the fractional domain by finding residual errors at different values of β . The results are consistent through out the fractional domain. It is also observed that increasing the value of the fractional parameter decreases the error. Figure 4 illustrates the three-dimensional plot of approximate solutions and errors at β = 1. The 3D solution curve (part a) depicts that at x and t = 0, the crest of the wave is highest, and as the distance and time increased, the altitude of the wave started to decrease. Figure 5 represents the effect of fractional parameter β on the water waves at x = 10. It is observed that for the time range between 0 and 1, increasing β decreases the water wave profile. Elevation and motion of the surface of water waves in two-dimensional form at a fixed time for various β can be seen in Figure 6. Near x = 0, the altitude of the wave was highest, but it started to smoothen out as the distance increased.
For the time-fractional Ito equation in Example 2, Table 5 is the comparison of LHPM with the modified Adomian decomposition method (mADM) and FCRPSA error at β = 1 and k = 0.01, whereas Table 6 displays the residual errors at different values of β . In both tables, it can be seen that LHPM is a reliable and powerful technique. Figure 7 is the LHPM solution and error graphs in 3D form, which displays the decrease in wave altitude with the increase of distance and time. The effect of different values of fractional parameter on the water waves surface at distance x = 10 can be seen in Figure 8, which shows that increasing β decreases the water wave level for the interval 0 < t < 1. A 2D plot to evaluate the motion of the wave surface is displayed in Figure 9 for k = 1, t = 3, and β = 0.2, 0.4, 0.8, and 1.0. It can be seen that nearby x = 0, the wave is at its peak but starts declining as the distance becomes larger.
In Example 3, a time-fractional Lax’s KdV equation is solved numerically and analytically by LHPM. A comparison of FCRPSA, mADM, and LHPM absolute errors at β = 1 in Table 1 shows that LHPM is more productive than mADM and FCRPSA. Residual errors for various β values are also shown in Table 2. A 3D LHPM solution plot and an error plot (Figure 1) are also displayed. From Figure 2 at k = 0.1 and x = 10, it is observed that increasing the value of the fractional parameter for the fifth-order Lax’s KdV model decreases the waves level. The effect of various β values on surface waves for fixed t = 3 is shown in 2D form in Figure 3. Analysis reveals that the behaviour of escalation of a wave is similar to that of Example 2. Table 7, Table 8 and Table 9 provide a comparison of LHPM errors with the optical optimal homotopy asymptotic method (OHAM) errors of the Kaup–Kupershmidt equation (Example 4) for β = 0.5, 0.75, and 1, respectively. Observation showed that LHPM is more reliable than OHAM. A 3D solution graph and error graph of the Kaup–Kupershmidt equation can be seen in Figure 10. For different β values, surface waves level in a 2D form for t = 3 is displayed in Figure 11 which depicts a sinusoidal cycle between x = −0.3 and 0.3. Moreover, wave level is declining at a greater value of x .

7. Conclusions

This paper is focused on the analysis of generalized fifth-order time-fractional KdV models through the Laplace transform along with homotopy perturbation. For checking the validity and efficiency of the proposed method, it is applied to the Sawada–Kotera, Ito, Lax, and Kaup–Kupershmidt KdV models in fractional sense, and residual errors are computed for different values of fractional parameter in the fractional domain. LHPM solutions are obtained without the imposition of any restrictions on the structure of models. The obtained approximate solutions and errors are illustrated in three-dimensional plots for reader convenience. Plots against different values of fractional parameter on water surface level are also displayed to provide better understanding of the models. A comparison of proposed and existing techniques affirm the efficiency and accuracy of LHPM over other methods. Hence, LHPM is helpful in managing complex non-linear, fractional, higher order KdV equations with improved accuracy.

Author Contributions

Conceptualization, M.Q. and M.B.R.; methodology, E.A.; software, J.A.; validation, J.A., M.Q. and M.B.R.; formal analysis, J.A.; investigation, E.A.; resources, M.Q.; data curation, M.Q.; writing—original draft preparation, E.A.; writing—review and editing, J.A., M.Q. and M.B.R.; visualization, M.Q.; supervision, M.B.R., M.Q.; project administration, J.A.; funding acquisition, J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Russell, J.S. Report on Waves. In Proceedings of the 14th Meeting of the British Association for the Advancement of Science, York, UK, September 1844. [Google Scholar]
  2. Boussinesq, J. Theorie de l’intumescence liquide appelee onde solitaire ou de translation se propageant dans un canal rectangulaire. CR Acad. Sci. Paris 1871, 72, 755–759. [Google Scholar]
  3. Boussinesq, J.V. Theorie generale des mouvements qui sont propages dans un canal rectangulaire horizontal. CR Acad. Sci. Paris 1871, 73, 256–260. [Google Scholar]
  4. Boussinesq, J. Theorie des ondes et des remous qui se propagent le long d’un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond. J. Math. Pures Appl. 1872, 338, 55–108. [Google Scholar]
  5. Korteweg, D.J.; de Vries, G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1895, 39, 422–443. [Google Scholar] [CrossRef]
  6. Al-Smadi, M.; Arqub, O.A.; Gaith, M. Numerical simulation of telegraph and Cattaneo fractional-type models using adaptive reproducing kernel framework. Math. Methods Appl. Sci. 2020, 40, 8472–8489. [Google Scholar] [CrossRef]
  7. Al-Smadi, M.; Arqub, O.A.; Momani, S. Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense. Phys. Scr. 2020, 95, 075218. [Google Scholar] [CrossRef]
  8. Kumar, S.; Kumar, A.; Samet, B.; Dutta, H. A study on fractional host-parasitoid population dynamical model to describe insect species. Numer. Methods Partial. Differ. Equ. 2020, 37, 1673–1692. [Google Scholar] [CrossRef]
  9. Ito, M. An Extension of Nonlinear Evolution Equations of the K-dV (mK-dV) Type to Higher Orders. J. Phys. Soc. Jpn. 1980, 49, 771–778. [Google Scholar] [CrossRef]
  10. Iyiola, O.S. A numerical study of ito equation and Sawada–Kotera equation both of time-fractional type. Adv. Math. Sci. J. 2013, 2, 71–79. [Google Scholar]
  11. Wang, J.; Xu, T.Z.; Wang, G.W. Numerical algorithm for time-fractional Sawada–Kotera equation and Ito equation with Bernstein polynomials. Appl. Math. Comput. 2018, 338, 1–11. [Google Scholar] [CrossRef]
  12. Lee, C. Some remarks on the fifth-order KdV equations. J. Math. Anal. Appl. 2015, 425, 281–294. [Google Scholar] [CrossRef]
  13. Bakodah, H.O. Modified Adomain Decomposition Method for the Generalized Fifth Order KdV Equations. Am. J. Comput. Math. 2013, 3, 53–58. [Google Scholar] [CrossRef] [Green Version]
  14. Arqub, O.A.; Al-Smadi, M.; Almusawa, H.; Baleanu, D.; Hayat, T.; Alhodaly, M.; Osman, M. A novel analytical algorithm for generalized fifth-order time-fractional nonlinear evolution equations with conformable time derivative arising in shallow water waves. Alex. Eng. J. 2022, 61, 5753–5769. [Google Scholar] [CrossRef]
  15. Guo, Y.; Li, D.; Wang, J. The new exact solutions of the Fifth-Order Sawada–Kotera equation using three wave method. Appl. Math. Lett. 2019, 94, 232–237. [Google Scholar] [CrossRef]
  16. Liu, T. Exact Solutions to Time-Fractional Fifth Order KdV Equation by Trial Equation Method Based on Symmetry. Symmetry 2019, 11, 742. [Google Scholar] [CrossRef] [Green Version]
  17. Park, C.; Nuruddeen, R.I.; Ali, K.K.; Muhammad, L.; Osman, M.S.; Baleanu, D. Novel hyperbolic and exponential ansatz methods to the fractional fifth-order Korteweg–de Vries equations. Adv. Differ. Equ. 2020, 2020, 1–12. [Google Scholar] [CrossRef]
  18. Singh, J.; Gupta, A.; Baleanu, D. On the analysis of an analytical approach for fractional Caudrey-Dodd-Gibbon equations. Alex. Eng. J. 2022, 61, 5073–5082. [Google Scholar] [CrossRef]
  19. Saratha, S.R.; Bagyalakshmi, M.; Krishnan, G.S.S. Fractional generalised homotopy analysis method for solving nonlinear fractional differential equations. Comput. Appl. Math. 2020, 39, 1–32. [Google Scholar] [CrossRef]
  20. Ismail, G.M.; Abdl-Rahim, H.R.; Ahmad, H.; Chu, Y.M. Fractional residual power series method for the analytical and approximate studies of fractional physical phenomena. Open Phys. 2020, 18, 799–805. [Google Scholar] [CrossRef]
  21. Ziane, D.; Cherif, M.H. Variational iteration transform method for fractional differential equations. J. Interdiscip. Math. 2018, 21, 185–199. [Google Scholar] [CrossRef]
  22. Guo, P. The Adomian Decomposition Method for a Type of Fractional Differential Equations. J. Appl. Math. Phys. 2019, 7, 2459–2466. [Google Scholar] [CrossRef] [Green Version]
  23. Song, H.; Yi, M.; Huang, J.; Pan, Y. Bernstein polynomials method for a class of generalized variable order fractional differential equations. IAENG Int. J. Appl. Math. 2016, 46, 437–444. [Google Scholar]
  24. Nadeem, M.; He, J.H. The homotopy perturbation method for fractional differential equations: Part 2, two-scale transform. Int. J. Numer. Methods Heat Fluid Flow 2021, 32, 559–567. [Google Scholar] [CrossRef]
  25. Veeresha, P.; Prakasha, D.G.; Singh, J. Solution for fractional forced KdV equation using fractional natural decomposition method. AIMS Math. 2020, 5, 798–810. [Google Scholar] [CrossRef]
  26. Xing, Z.; Wen, L.; Wang, W. An efficient difference scheme for time-fractional KdV equation. Comput. Appl. Math. 2021, 40. [Google Scholar] [CrossRef]
  27. Alderremy, A.A.; Aly, S.; Fayyaz, R.; Khan, A.; Shah, R.; Wyal, N. The Analysis of Fractional-Order Nonlinear Systems of Third Order KdV and Burgers Equations via a Novel Transform. Complexity 2022, 2022, 4935809. [Google Scholar] [CrossRef]
  28. He, J.H. Homotopy perturbation technique. Comput. Methods Appl. Mech. Eng. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  29. He, J.H.; El-Dib, Y.O. The enhanced homotopy perturbation method for axial vibration of strings. Facta Univ. Ser. Mech. Eng. 2021, 19, 735. [Google Scholar] [CrossRef]
  30. Anjum, N.; He, J.H.; Ain, Q.T.; Tian, D. He’s modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system. Facta Univ. Ser. Mech. Eng. 2021, 19, 601. [Google Scholar] [CrossRef]
  31. Thabet, H.; Kendre, S. Modified least squares homotopy perturbation method for solving fractional partial differential equations. Malaya J. Mat. 2018, 6, 420–427. [Google Scholar] [CrossRef] [Green Version]
  32. Habib, S.; Batool, A.; Islam, A.; Nadeem, M.; Gepreel, K.A.; He, J.H. Study of nonlinear Hirota-Satsuma coupled kdv and coupled mkdv system with time fractional derivative. Fractals 2021, 29, 2150108. [Google Scholar] [CrossRef]
  33. Johnston, S.J.; Jafari, H.; Moshokoa, S.P.; Ariyan, V.M.; Baleanu, D. Laplace homotopy perturbation method for Burgers equation with space- and time-fractional order. Open Phys. 2016, 14, 247–252. [Google Scholar] [CrossRef]
  34. Yin, X.B.; Kumar, S.; Kumar, D. A modified homotopy analysis method for solution of fractional wave equations. Adv. Mech. Eng. 2015, 7, 168781401562033. [Google Scholar] [CrossRef] [Green Version]
  35. Bohner, M.; Tun, O.; Tun, C. Qualitative analysis of caputo fractional integro-differential equations with constant delays. Comput. Appl. Math. 2021, 40, 1–17. [Google Scholar] [CrossRef]
  36. Tun, C.; Tun, O. On the stability, integrability and boundedness analyses of systems of integro-differential equations with time-delay retardation. Rev. Real Acad. Cienc. Exactas Fis. Nat. Serie A Mat. 2021, 115. [Google Scholar] [CrossRef]
  37. Gupta, A.K.; Ray, S.S. The comparison of two reliable methods for accurate solution of time-fractional Kaup–Kupershmidt equation arising in capillary gravity waves. Math. Methods Appl. Sci. 2015, 39, 583–592. [Google Scholar] [CrossRef]
Figure 1. 3D graphical illustration of LHPM solution (a) and error (b) in Example 3 when β = 1 and k = 0.15.
Figure 1. 3D graphical illustration of LHPM solution (a) and error (b) in Example 3 when β = 1 and k = 0.15.
Universe 08 00563 g001
Figure 2. Effect of fractional parameter β on time-fractional Lax’s KdV equation (Example 3) when k = 0.1 and x = 10.
Figure 2. Effect of fractional parameter β on time-fractional Lax’s KdV equation (Example 3) when k = 0.1 and x = 10.
Universe 08 00563 g002
Figure 3. 2D surface formation at different values of β in Example 3 when k = 1, t = 3.
Figure 3. 2D surface formation at different values of β in Example 3 when k = 1, t = 3.
Universe 08 00563 g003
Figure 4. 3D graphical illustration of LHPM solution (a) and error (b) in Example 1 when β = 1 and k = 0.15.
Figure 4. 3D graphical illustration of LHPM solution (a) and error (b) in Example 1 when β = 1 and k = 0.15.
Universe 08 00563 g004
Figure 5. Effect of fractional parameter β on time-fractional Sawada–Kotera equation (Example 1) when k = 0.1 and x = 10.
Figure 5. Effect of fractional parameter β on time-fractional Sawada–Kotera equation (Example 1) when k = 0.1 and x = 10.
Universe 08 00563 g005
Figure 6. 2D surface formation at different values of β in Example 1 when k = 1, t = 3.
Figure 6. 2D surface formation at different values of β in Example 1 when k = 1, t = 3.
Universe 08 00563 g006
Figure 7. 3D graphical illustration of LHPM solution (a) and error (b) in Example 2 when β = 1 and k = 0.15.
Figure 7. 3D graphical illustration of LHPM solution (a) and error (b) in Example 2 when β = 1 and k = 0.15.
Universe 08 00563 g007
Figure 8. Effect of fractional parameter β on a time-fractional Ito equation (Example 2) when k = 0.1 and x = 10.
Figure 8. Effect of fractional parameter β on a time-fractional Ito equation (Example 2) when k = 0.1 and x = 10.
Universe 08 00563 g008
Figure 9. 2D surface formation at different values of β in Example 2 when k = 1 , t = 3 .
Figure 9. 2D surface formation at different values of β in Example 2 when k = 1 , t = 3 .
Universe 08 00563 g009
Figure 10. 3D graphical illustration of LHPM solution (a) and error (b) in Example 4 when β = 1, λ = 0.1, w = 1, ρ = 2.5, and μ = 0.
Figure 10. 3D graphical illustration of LHPM solution (a) and error (b) in Example 4 when β = 1, λ = 0.1, w = 1, ρ = 2.5, and μ = 0.
Universe 08 00563 g010
Figure 11. 2D surface formation at different values of β in Example 4 when λ = 1, w = 1, ρ = 2.5, μ = 1 and t = 3.
Figure 11. 2D surface formation at different values of β in Example 4 when λ = 1, w = 1, ρ = 2.5, μ = 1 and t = 3.
Universe 08 00563 g011
Table 1. Comparison of LHPM errors with FCRPSA and mADM errors in Example 3 when β = 1 and k = 0.01.
Table 1. Comparison of LHPM errors with FCRPSA and mADM errors in Example 3 when β = 1 and k = 0.01.
t x Exact SolutionLHPM SolutionLHPM ErrorFCRPSA Error [14]mADM Error [13]
20.000399760.000399765.42 × 10 20 1.34 × 10 17 2.30 × 10 13
40.000399040.0003990401.30 × 10 17 4.60 × 10 13
0.860.000397840.0003978401.27 × 10 17 6.19 × 10 13
80.000396170.0003961701.21 × 10 17 9.21 × 10 13
100.000394040.0003940401.14 × 10 17 1.15 × 10 12
20.000399760.000399765.42 × 10 20 5.17 × 10 16 1.42 × 10 12
40.000399040.000399045.42 × 10 20 5.08 × 10 16 2.84 × 10 12
560.000397840.000397845.42 × 10 20 4.94 × 10 16 4.26 × 10 12
80.000396170.000396175.42 × 10 20 4.73 × 10 16 5.68 × 10 12
100.000394040.0003940404.49 × 10 16 7.10 × 10 12
Table 2. Error analysis of LHPM in fractional domain for time-fractional Lax’s KdV model (Example 3) when k = 0.14 and x = 4 .
Table 2. Error analysis of LHPM in fractional domain for time-fractional Lax’s KdV model (Example 3) when k = 0.14 and x = 4 .
t β = 0.2 β = 0.4 β = 0.6 β = 0.8
0.1 3.15 × 10 8 1.35 × 10 9 3.43 × 10 11 4.87 × 10 13
0.3 9.45 × 10 8 1.22 × 10 8 9.27 × 10 10 3.94 × 10 11
0.5 1.57 × 10 7 3.39 × 10 8 4.29 × 10 9 3.04 × 10 10
0.7 2.20 × 10 7 6.65 × 10 8 1.17 × 10 8 1.16 × 10 9
0.9 2.83 × 10 7 1.10 × 10 7 2.50 × 10 8 3.19 × 10 9
Table 3. Comparison of LHPM and FCRPSA errors in Example 1 when β = 1 and k = 0.01.
Table 3. Comparison of LHPM and FCRPSA errors in Example 1 when β = 1 and k = 0.01.
t x Exact SolutionLHPM SolutionLHPM ErrorFCRPSA Error [14]
20.000199920.000199921.38 × 10 20 3.46 × 10 18
40.000199680.000199683.13 × 10 21 3.46 × 10 18
0.160.000199280.000199282.57 × 10 20 1.04 × 10 17
80.000198720.000198721.38 × 10 20 1.38 × 10 17
100.000198010.000198014.48 × 10 20 1.73 × 10 17
20.000199920.000199926.59 × 10 20 2.70 × 10 16
40.000199680.000199683.86 × 10 20 8.22 × 10 16
0.560.000199280.000199282.05 × 10 20 1.35 × 10 15
80.000198720.000198721.18 × 10 20 1.88 × 10 15
100.000198010.000198016.15 × 10 20 2.39 × 10 15
20.000199920.000199923.71 × 10 20 1.59 × 10 15
40.000199680.000199682.63 × 10 20 4.75 × 10 15
0.960.000199280.000199281.54 × 10 20 7.90 × 10 15
80.000198720.000198724.39 × 10 20 1.09 × 10 14
100.000198010.000198015.10 × 10 20 1.40 × 10 14
Table 4. Error analysis of LHPM in fractional domain for time-fractional Sawada–Kotera model (Example 1) when k = 0.15 and x = 1 .
Table 4. Error analysis of LHPM in fractional domain for time-fractional Sawada–Kotera model (Example 1) when k = 0.15 and x = 1 .
t β = 0.2 β = 0.4 β = 0.6 β = 0.8
0.1 9.50 × 10 9 4.09 × 10 10 1.03 × 10 11 1.47 × 10 13
0.3 2.85 × 10 8 3.68 × 10 9 2.79 × 10 10 1.19 × 10 11
0.5 4.75 × 10 8 1.02 × 10 8 1.29 × 10 9 9.19 × 10 11
0.7 6.65 × 10 8 2.00 × 10 8 3.55 × 10 9 3.53 × 10 10
0.9 8.55 × 10 8 3.31 × 10 8 7.55 × 10 9 9.64 × 10 10
Table 5. Comparison of LHPM errors with FCRPSA and mADM in Example 2 when β = 1 and k = 0.01.
Table 5. Comparison of LHPM errors with FCRPSA and mADM in Example 2 when β = 1 and k = 0.01.
t x Exact SolutionLHPM SolutionLHPM ErrorFCRPSA Error [14]mADM Error [13]
20.00199880.00199882.11 × 10 22 2.77 × 10 17 1.41 × 10 16
40.00199520.00199528.47 × 10 22 3.88 × 10 16 5.63 × 10 16
0.260.00198920.00198925.08 × 10 21 4.13 × 10 16 1.26 × 10 15
80.00198080.00198083.38 × 10 21 1.04 × 10 15 2.25 × 10 15
100.00197020.001970202.31 × 10 15 3.52 × 10 15
20.00199880.0019988001.41 × 10 16
40.00199520.00199528.47 × 10 22 3.60 × 10 16 5.64 × 10 16
0.660.00198920.00198925.08 × 10 21 4.05 × 10 16 1.27 × 10 15
80.00198080.00198083.38 × 10 21 9.90 × 10 16 2.25 × 10 15
100.00197020.00197023.38 × 10 21 2.20 × 10 15 3.52 × 10 15
20.00199880.001998802.77 × 10 17 1.41 × 10 16
40.00199520.00199521.69 × 10 21 3.33 × 10 16 5.66 × 10 16
1.060.00198920.00198921.69 × 10 21 3.94 × 10 16 1.27 × 10 15
80.00198080.00198083.38 × 10 21 8.99 × 10 16 2.26 × 10 15
100.00197020.00197023.38 × 10 21 2.00 × 10 15 3.53 × 10 15
Table 6. Error analysis of LHPM in a fractional domain for a time-fractional Ito model (Example 2) when k = 0.14 and x = 1 .
Table 6. Error analysis of LHPM in a fractional domain for a time-fractional Ito model (Example 2) when k = 0.14 and x = 1 .
t β = 0.2 β = 0.4 β = 0.6 β = 0.8
0.1 2.16 × 10 7 9.34 × 10 9 2.36 × 10 10 3.35 × 10 12
0.3 6.50 × 10 7 8.41 × 10 8 6.38 × 10 9 2.71 × 10 10
0.5 1.08 × 10 6 2.33 × 10 7 2.95 × 10 8 2.09 × 10 9
0.7 1.51 × 10 6 4.58 × 10 7 8.11 × 10 8 8.05 × 10 9
0.9 1.95 × 10 6 7.57 × 10 7 1.72 × 10 7 2.20 × 10 8
Table 7. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt model (Example 4) at different t when β = 0.5, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
Table 7. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt model (Example 4) at different t when β = 0.5, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
x t = 0.3 t = 0.5 t = 0.8
OHAM [37]LHPMOHAM [37]LHPMOHAM [37]LHPM
0.2 7.23 × 10 6 4.22 × 10 9 7.21 × 10 6 5.41 × 10 9 7.18 × 10 6 6.78 × 10 9
0.4 5.79 × 10 5 8.55 × 10 9 5.78 × 10 5 1.10 × 10 8 5.78 × 10 5 1.38 × 10 8
0.6 1.53 × 10 4 1.28 × 10 8 1.53 × 10 4 1.65 × 10 8 1.53 × 10 4 2.08 × 10 8
0.8 2.86 × 10 4 1.71 × 10 8 2.86 × 10 4 2.20 × 10 8 2.85 × 10 4 2.78 × 10 8
Table 8. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt equation (Example 4) at different t when β = 0.75, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
Table 8. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt equation (Example 4) at different t when β = 0.75, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
x t = 0.3 t = 0.5 t = 0.8
OHAM [37]LHPMOHAM [37]LHPMOHAM [37]LHPM
0.2 7.22 × 10 6 3.05 × 10 9 7.17 × 10 6 4.44 × 10 9 7.12 × 10 6 6.26 × 10 9
0.4 5.79 × 10 5 6.14 × 10 9 5.78 × 10 5 8.98 × 10 9 5.76 × 10 5 1.27 × 10 8
0.6 1.53 × 10 4 9.21 × 10 9 1.53 × 10 4 1.34 × 10 8 1.52 × 10 4 1.91 × 10 8
0.8 2.86 × 10 4 1.22 × 10 8 2.85 × 10 4 1.79 × 10 8 2.85 × 10 4 2.55 × 10 8
Table 9. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt model (Example 4) at different t when β = 1, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
Table 9. Comparison of LHPM and OHAM errors in the Kaup–Kupershmidt model (Example 4) at different t when β = 1, λ = 0.1, w = 1, ρ = 2.5 and μ = 0.
x t = 0.3 t = 0.5 t = 0.8
OHAM [37]LHPMOHAM [37]LHPMOHAM [37]LHPM
0.2 7.21 × 10 6 2.09 × 10 9 7.14 × 10 6 3.46 × 10 9 7.03 × 10 6 5.05 × 10 9
0.4 5.79 × 10 5 4.19 × 10 9 5.77 × 10 5 6.97 × 10 9 5.74 × 10 5 1.11 × 10 8
0.6 1.53 × 10 4 6.28 × 10 9 1.52 × 10 4 1.04 × 10 8 1.52 × 10 4 1.66 × 10 8
0.8 2.86 × 10 4 8.36 × 10 9 2.85 × 10 4 1.39 × 10 8 2.85 × 10 4 2.22 × 10 8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qayyum, M.; Ahmad, E.; Riaz, M.B.; Awrejcewicz, J. Improved Soliton Solutions of Generalized Fifth Order Time-Fractional KdV Models: Laplace Transform with Homotopy Perturbation Algorithm. Universe 2022, 8, 563. https://doi.org/10.3390/universe8110563

AMA Style

Qayyum M, Ahmad E, Riaz MB, Awrejcewicz J. Improved Soliton Solutions of Generalized Fifth Order Time-Fractional KdV Models: Laplace Transform with Homotopy Perturbation Algorithm. Universe. 2022; 8(11):563. https://doi.org/10.3390/universe8110563

Chicago/Turabian Style

Qayyum, Mubashir, Efaza Ahmad, Muhammad Bilal Riaz, and Jan Awrejcewicz. 2022. "Improved Soliton Solutions of Generalized Fifth Order Time-Fractional KdV Models: Laplace Transform with Homotopy Perturbation Algorithm" Universe 8, no. 11: 563. https://doi.org/10.3390/universe8110563

APA Style

Qayyum, M., Ahmad, E., Riaz, M. B., & Awrejcewicz, J. (2022). Improved Soliton Solutions of Generalized Fifth Order Time-Fractional KdV Models: Laplace Transform with Homotopy Perturbation Algorithm. Universe, 8(11), 563. https://doi.org/10.3390/universe8110563

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