Next Article in Journal
Modeling Double Stochastic Opinion Dynamics with Fractional Inflow of New Opinions
Previous Article in Journal
Enhancing Medical Image Quality Using Fractional Order Denoising Integrated with Transfer Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis Study of Time-Fractional Nonlinear Modified Kawahara Equation Based on the Homotopy Perturbation Sadik Transform

1
Institute of Computing Science and Technology, Guangzhou University, Guangzhou 510006, China
2
Wuhan Huamao Automation Co., Ltd., Wuhan 430074, China
3
Department of Computer Engineering and Information, College of Engineering in Wadi Alddawasir, Prince Sattam Bin Abdulaziz University, Wadi Alddawasir 11991, Saudi Arabia
4
Department of Industrial Engineering, Apadana Institute of Higher Education, Shiraz 7187985443, Iran
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(9), 512; https://doi.org/10.3390/fractalfract8090512
Submission received: 19 July 2024 / Revised: 24 August 2024 / Accepted: 26 August 2024 / Published: 29 August 2024

Abstract

:
In this manuscript, we survey a numerical algorithm based on the combination of the homotopy perturbation method and the Sadik transform for solving the time-fractional nonlinear modified shallow water waves (called Kawahara equation) within the frame of the Caputo–Prabhakar (CP) operator. The nonlinear terms are handled with the assistance of the homotopy polynomials. The stability analysis of the implemented method is studied by using S-stable mapping and the Banach contraction principle. Also, we use the fixed-point method to determine the existence and uniqueness of solutions in the given suggested model. Finally, some numerical simulations are illustrated to display the accuracy and efficiency of the present numerical method. Moreover, numerical behaviors are captured to validate the reliability and efficiency of the scheme.

1. Introduction

The integral and derivatives of fractional order play a vital role and are basic in many branches of engineering and physical sciences [1,2,3,4,5]. The differential equations containing fractional operators have substantial applications in different scopes of applied sciences [6,7], engineering [8,9], and physics [10,11,12,13,14,15]. The Kawahara differential models of fractional order are important nonlinear models that play a vital role in mathematical and engineering sciences. These partial differential models describe physical phenomena such as water transfer and shallow water waves. Due to the important applications of these types of nonlinear equations in physical sciences [16], these equations are considered in this work. This paper studies the fractional expansion of the modified Kawahara model by applying the Caputo–Prabhakar fractional derivative in the sense of the three-parameter Mittag–Leffler function. The Kawahara model is one of the several influential models and is important in ocean engineering and physics. To study the water wave equations in the long-wave regime, it was introduced and developed by Kawahara as a type of modest value of Weber numbers (WNs) near 1 3 and the surface tension [17]. Moreover, ref. [18] and Amick and Kirchgässner [19] derived the stimulating results associated with the model of water wave problems. The Kawahara equation and its modified forms have recently been discussed by many authors [20,21]. In this manuscript, we consider the following time-fractional modified equal-width wave equation [22,23,24]:
  C D t μ z ( x , t ) + z 2 ( x , t ) z x ( x , t ) + c z x x x ( x , t ) + d z x x x x x = 0 , z ( x , 0 ) = z 0 ( x ) ,
in which ( μ ) > 0 , and c and d are non-zero real constants. Here, the symbol   C D t μ is Caputo–Prabhakar fractional operator of order μ and is defined by the following formula:
  C D t μ z ( x , t ) = I 1 μ , η d d t z ( x , t ) ,
in which I 1 μ , η is the Prabhakar fractional operator and given as [25]
( I μ , η z ) ( x , t ) = 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) ρ ) z ( x , s ) d s , σ , ρ C ,
in which E λ , μ η is the Mittag–Leffler (ML) function and for z C , defined by [25,26]
E λ , μ η ( z ) = n = 0 Γ ( η + n ) Γ ( η ) n ! Γ ( λ n + μ ) z n , ( λ ) , ( μ ) > 0 , η > 0 .
In general, for all the fractional differential equations (FDEs), the process of finding the analytical solution is not easy to compute. Hence, researchers are drawn towards obtaining an approximate or numerical solution of these types of equations, which they find to be valuable. To gain the approximate solutions of Equation (1), we express a substantial and reliable method based on a combination of the Sadik transform (ST) with the homotopy perturbation (HP) technique. Prabhakar [25] showed a novel derivative with the ML kernel as an extension of Riemann–Liouville and Caputo derivatives and displayed the essence of generalizing the models associated with probability and mathematical physics. The main reason for choosing the ML function in this manuscript is related to its main applications in various models, such as heterogeneous [27], Havriliak–Negami [28,29], viscoelastic [30], and stochastic [31] models studied in physical sciences [32,33,34].
One of the numerical methods for solving the FDEs is the Sadik transform. Recently, Shaikh [35,36] discussed a modern integral representation called the Sadik transform (ST). In [35], Shaikh found some properties of this transform, like the existing theorem and duality theorem of the Sadik transform. Moreover, the author illustrated that Sumudu, Laplace, and Laplace–Carson transforms are particular cases of ST. Soon after, Shaikh in [37] presented a shifting theorem for ST and its properties for the derivative of functions. Also, for the dynamical system, the transfer function is obtained using control theory via ST. The integral transform technique is an effective tool and useful for solving fractional differential equations. However, it is difficult to find the exact solution for all kinds of FDEs using the integral transform method (see [38,39]). In this connection, the authors of [40] studied some basic notions and properties of ST, and then employed them to investigate nonlinear models. The integral transformation method with different structures has been performed to solve numerous types of differential equations in various scientific fields, for example, in the Adomian decomposition method [41,42,43], homotopy-based schemes [44,45,46], the iterative perturbation method [47], numerical methods based on the Petrov–Galerkin [48], the Petrov–Galerkin finite method [49], the variational iteration method [50], and the methods studied in the references [22,51,52,53,54,55].
With this aim, the rest of this paper is presented as follows. Section 2 contains a summary of the definitions and mathematical preliminaries of fractional calculus. We describe the suggested method for solving Equation (1) in Section 3. The convergence, stability, existence, and uniqueness of the suggested method are studied in Section 4. To verify and demonstrate the proposed approach, we solve the proposed model in Section 5. Finally, some conclusions are expressed in Section 6.

2. Important Preliminaries

This part studies the theorems and lemmas which will be applied in the following sections.
Definition 1
([56]). Let ψ ( t ) be a piecewise continuous function under the interval [ 0 , B ] such that | ψ ( t ) | C 1 exp ( a t ) , where B > 0 , a > 0 , and C 1 > 0 are the real constants. Then, the Sadik transform for ψ ( t ) is given by
T ( ψ ( t ) ) = 1 υ β 0 + e t υ α ψ 1 ( t ) d t .
Lemma 1
([4]). Under the assumptions of Equation (4), we have
( I μ , η t ν 1 ) ( x ) = Γ ( ν ) t μ + ν 1 E λ , μ + ν η ( σ t λ ) .
Lemma 2
([56]). Suppose that Φ α , β υ is an ST of ψ 1 ( t ) and Ψ α , β υ is an ST of ψ 2 ( t ) . Then, ST of ( ψ 1 * ψ 2 ) ( t ) is displayed by
T [ ( ψ 1 * ψ 2 ) ( t ) ] = υ β Φ α , β υ . Ψ α , β υ ,
in which the symbol * shows the convolution operator. Here, we consider Φ α , β υ = T ( ψ 1 ( t ) ) = 1 υ β 0 + e t υ α ψ 1 ( t ) d t and Ψ α , β υ = T ( ψ 2 ( t ) ) = 1 υ β 0 + e t υ α ψ 2 ( t ) d t .
Lemma 3.
Under the assumptions of Lemma 1, we have
T t μ 1 E λ , μ η ( σ t λ ) ; υ = υ α μ β ( 1 σ υ λ α ) η ,
in which T is the Sadik transform of function.
Proof. 
With the help of the definition of ST and using Equation (4), we obtain
T ( t μ 1 E λ , μ η ( σ t λ ) ; υ ) = 1 υ β 0 + e t υ α t μ 1 E λ , μ η ( σ t λ ) d t = 1 υ β Q ( λ , μ , η ) 0 + e t υ α t λ n + μ 1 d t .
in which Q ( λ , μ , η ) = n = 0 Γ ( η + n ) σ n Γ ( λ n + μ ) Γ ( η ) n ! . By the use of the gamma function, the above equations can be expressed as follows:
T ( t μ 1 E λ , μ η ( σ t λ ) ; υ ) = 1 υ α μ + β n = 0 η n ( σ υ λ α ) n n ! = υ α μ β ( 1 σ υ λ α ) η .
In this case, Equation (7) is obtained. □
Theorem 1
([56]). Suppose Φ α , β υ = T ( ψ ( i ) ( t ) ) and ψ ( i ) ( t ) , i N { 0 } are continuous. Then,
T [ ψ i ( t ) ] = υ i α Φ α , β υ k = 0 i 1 υ k α β ψ i k 1 ( 0 ) .
Also, the integration of φ ( i ) ( t ) is achieved by using the ST displayed as follows:
T [ 0 t ψ ( t ) ( i ) ] = 1 υ α Φ α , β υ .
Theorem 2
The ST of fractional derivative defined in (2) for i = 1 is calculated by
T   C D t μ z ( x , t ) = υ α μ ( 1 σ υ λ α ) η Z ( υ , α , β ) υ α μ ( α + β ) ( 1 σ υ λ α ) η z ( x , 0 + ) ,
in which Z ( υ , α , β ) = T ( z ( x , t ) ) .
Proof. 
By using Equations (2), (6), (7) and (10), we have
T   C D t μ z ( x , t ) = 1 υ β 0 e t υ α   C D t μ z ( x , t ) d t = 1 υ β 0 e t υ α 0 t ( t s ) μ E λ , 1 μ η ( σ ( t s ) λ ) d d s z ( x , s ) d s d t = 1 υ β 0 s e t υ α ( t s ) μ E λ , 1 μ η ( σ ( t s ) λ ) d d s z ( x , s ) d t d s = 1 υ β 0 e s υ α d d s z ( x , s ) 0 e y υ α y μ E λ , 1 μ η ( σ y λ ) d y d s = υ β T d d t z ( x , t ) × T t μ E λ , 1 μ η ( σ t λ ) = υ α Z ( υ , α , β ) υ β z ( x , 0 + ) υ α ( μ 1 ) ( 1 σ υ λ α ) η = υ α μ ( 1 σ υ λ α ) η Z ( υ , α , β ) υ α μ ( α + β ) ( 1 σ υ λ α ) η z ( x , 0 + ) .
The proof is completed. □
Lemma 4
([12]). Let S : X X be a mapping and ( X , . ) be a Banach space. Then, for w , v X
S w S v σ 1 w S w + η 1 w v ,
in which σ 1 0 and 0 η 1 1 . Moreover, there is a fixed point of S and it is Picard’s S -stable.

3. The Projected Scheme

This part demonstrates the HP method and the ST for obtaining the following equation solutions:
  C D t μ z ( x , t ) + M z ( x , t ) + V z ( x , t ) = w ( x , t ) , z ( x , 0 ) = z 0 ( x ) ,
in which M z ( x , t ) is a linear part of Equation (14), V z ( x , t ) is the nonlinear part of Equation (14), and w ( x , t ) is a known function. By applying the ST and Equation (12), one can obtain
υ α μ ( 1 σ υ λ α ) η Z ( υ , α , β ) υ α μ ( α + β ) ( 1 σ υ λ α ) η z 0 ( x ) + T M z ( x , t ) + V z ( x , t ) = T ( w ( x , t ) ) .
By using the Sadik inverse transform, we obtain
z ( x , t ) = T 1 υ α μ ( 1 σ υ λ α ) η T ( w ( x , t ) ) + 1 υ α + β z 0 ( x ) T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M z ( x , t ) + V z ( x , t ) ] .
To solve Equation (16), we consider the set of solutions by
z = n = 0 β n z n ,
where β n is the unknown coefficient and is calculated from Equation (16). Now, we consider the nonlinear part V z by
V z = n = 0 β n F n ( z ) ,
in which F n ( z ( x , t ) ) = 1 n ! n β n N i = 0 ϱ i z i β = 0 is the Adomian function that is defined in [57]. By inputting Equations (17) and (18) into Equation (16), we obtain
n = 0 β n z n = W T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M n = 0 β n z n + V n = 0 β n z n = W T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M n = 0 β n z n + n = 0 β n F n ( z ( x , t ) )
in which W = T 1 1 υ α μ ( 1 σ υ λ α ) η T ( w ) + 1 υ α + β z 0 ( x ) . For n = 1 , 2 , , we have
β 0 : z 0 = W ( x , t ) , β 1 : z 1 = T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M z 0 + F 0 ( z ) ; t , β 2 : z 2 = T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M z 1 + F 1 ( z ) ; t ,   β n + 1 : z n + 1 = T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T M z n + F n ( z ) ; t .
Therefore, the solution of Equation (14) is obtained as
z = n = 0 z n .

4. Existence and Uniqueness of Solutions with Caputo–Prabhakar Derivative

In this section, we used the fixed-point theorem to show the existence and uniqueness of solutions for the time-fractional modified Kawahara equation. To this aim, by applying the Prabhakar fractional integral given in Equation (3) on both sides of Equation (1), we obtain
z ( x , t ) z ( x , 0 ) = [ 0 t ( t s ) μ 1 E λ , μ η ( ω ( t s ) λ ) ( z 2 ( x , s ) z x ( x , s ) + c z x x x ( x , s ) + d z x x x x x ( x , s ) ) ] d s .
For simplicity, we assume that
ζ ( z , x , t ) = z 2 ( x , t ) z x ( x , t ) + c z x x x ( x , t ) + d z x x x x x ( x , t ) .
Then, we have
z ( x , t ) z 0 = 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z , x , s ) d s .
Now, we show that the function ζ ( Z , x , t ) admits the Lipschitz condition when ζ ( Z , x , t ) is bounded. For this aim, we assume that the continuous functions U and V are upper bounded; thus, we have
ζ ( U , x , t ) ζ ( V , x , t ) = ( U 2 U x V 2 V x ) + c ( U x x x V x x x ) + d ( U x x x x V x x x x ) U 2 U x V 2 V x + c U x x x V x x x + d U x x x x x V x x x x x = 1 3 x ( U 3 V 3 ) + a x 3 ( U V ) + d x 5 ( U V ) = 1 3 x ( U V ) ( U 2 + U V + V 2 ) + c x 3 ( U V ) + d x 5 ( U V ) c L 1 3 ( l 1 2 + l 1 l 2 + l 2 2 ) U V + c L 2 U V + d L 3 U V = c L 1 3 ( l 1 2 + l 1 l 2 + l 2 2 ) + + c L 2 + d L 3 U V = ϱ U V ,
in which U l 1 , V l 2 and α = c L 1 3 ( l 1 2 + l 1 l 2 + l 2 2 ) + + c L 2 + d L 3 . Then, the Lipschitz condition for ζ ( Z , x , t ) is satisfied and if 0 α < 1 , it is a contraction. From Equation (23), we use the recursive relation
z n = 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z n 1 , x , s ) d s ,
with initial condition z 0 = z ( x , 0 ) . Furthermore, we consider the following differences:
Ψ n + 1 ( x , t ) = z n + 1 z n 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ( ζ ( z n , x , τ ) ζ ( z n 1 , x , s ) ) d s ,
by taking the norm on both sides of Equation (25), one can obtain
Ψ n + 1 ( x , t ) = z n + 1 z n 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ( ζ ( z n , x , s ) ζ ( z n 1 , x , s ) ) d τ 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z n , x , s ) ζ ( z n 1 , x , s ) d s .
Applying Equation (27), we show the existence of a solution for the projected equation.
Theorem 3.
The time-fractional modified Kawahara equation with Caputo–Prabhakar derivative (1) has a solution, provided that the following holds true:
M ϱ < 1 .
Proof. 
Let the function Φ n ( x , t ) = Z n + 1 Z + Z ( x , 0 ) . Then, by applying Equations (24) and (27) and Lemma 1, we obtain
Φ n ( x , t ) = z n + 1 z + z ( x , 0 ) 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ( ζ ( z n , x , s ) ζ ( z , x , s ) ) d s 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z n , x , s ) ζ ( z , x , s ) d s = M ϱ z n z ( M ϱ ) n z z 1 ,
in which Lemma 2.1 in [6] ensures the existence of E λ , μ + 1 η ( σ t λ ) . Then, we can put M = E λ , μ + 1 η ( σ t λ ) . Thus, Equation (29) implies that the functions Φ n ( x , t ) 0 when n for M ϱ < 1 , which further gives that lim n z n + 1 = z . Consequently, solutions of the time-fractional modified Kawahara Equation (1) exist. □
Theorem 4.
Equation (1) has a unique solution when true:
1 M ϱ 0 .
Proof. 
For the uniqueness solution of the time-fractional modified Kawahara equation with Caputo–Prabhakar derivative (1), we use the contrary path for the proof. Suppose that u 1 is the other solution of (1); then, we have
z z 1 = 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z z 1 , x , s ) d s .
Therefore, we use the norm, for which we thus have
z z 1 = 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z z 1 , x , s ) d s 0 t ( t s ) μ 1 E λ , μ η ( σ ( t s ) λ ) ζ ( z z 1 , x , s ) d s M ϱ z z 1 ,
which implies
z z 1 1 M ϱ 0 .
Thus, Equation (33) implies that the functions z z 1 0 when 1 M ϱ 0 . Consequently, z = z 1 . Thus, the solution of the time-fractional modified Kawahara equation with Caputo–Prabhakar derivative (1) is unique. □
Theorem 5.
Suppose that z and z n are the exact and numerical solutions of Equation (14), such that z n + 1 ϖ z n , ϖ ( 0 , 1 ) , n N . Then, the series given by (21) converges.
Proof. 
From z n + 1 ϖ z n , we obtain
z 1 ϖ z 0 = ϖ z 0 ( x ) , z 2 ϖ z 1 ϖ ϖ z 0 ( x ) , z n + 1 ϖ z n ϖ n + 1 z 0 ( x ) .
Then,
i = n + 1 z i ( x , t ) i = n + 1 ϖ i z 0 ( x ) = z 0 ( x ) i = n + 1 ϖ i .
Therefore,
z z n + 1 = i = n + 1 z i i = n + 1 z i i = n + 1 ϖ i z 0 ( x ) = z 0 ( x ) i = n + 1 ϖ i ϖ n + 1 1 ϖ z 0 ( x ) .
Because ϖ ( 0 , 1 ) , then z z n 0 as n . □

Stability Analysis

In this part, we study the stability of the projected scheme. To achieve the stability of the method proposed in Equation (1), we define a self-mapping T as follows:
T ( z i ) = z i + 1 = T 1 1 υ α + β z 0 ( x ) T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c z x x x + d z x x x x x + z 2 ( x , t ) z x ] = z i T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c z x x x + d z x x x x x + z 2 z x ] ,
where T ( 1 ) = 1 υ α + β . To obtain the required result, we consider the following for ( i , j ) N × N :
T ( z i ) T ( z j ) = z i z j [ T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c ( z i ) x x x ( z j ) x x x + d ( z i ) x x x x x ( z j ) x x x x x + z i 2 ( z i ) x z j 2 ( z j ) x ] ] = z i ( x , t ) z j ( x , t ) [ T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c 3 3 x 3 z i u j + d 5 x 5 z i z j + 1 3 x z i 3 z j 3 ] ] .
On employing the norm on Equation (38) and after simplifying, one can obtain
T ( z i ) T ( z j ) z i z j + [ T 1 [ ( υ α μ ( 1 σ υ λ α ) η ) × T c 3 3 x 3 z i z j + d 5 x 5 z i z j + 1 3 x z i 3 z j 3 ] ] .
We let the continuous function z i be upper bounded; then,
T ( z i ( x , t ) ) T ( z j ( x , t ) ) c 3 β H 1 ( λ , μ , η ) + d β 1 H 2 ( λ , μ , η ) + 1 3 β 2 ( γ 1 2 + γ 1 γ 2 + γ 2 2 ) H 3 ( ρ , μ , γ ) z i ( x , t ) z j ( x , t ) ,
in which z i ( x , t ) γ 1 , z j ( x , t ) γ 2 , β = 3 x 3 , β 1 = 5 x 5 , β 2 = x and H k ( λ , μ , η ) , k = 1 , 2 , 3 are functions retrieved from T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) T . . The given operator in Equation (40) is a contraction mapping if the following relation is satisfied:
c 3 β H 1 ( λ , μ , η ) + d β 1 H 2 ( λ , μ , η ) + 1 3 β 2 ( γ 1 2 + γ 1 γ 2 + γ 2 2 ) H 3 ( ρ , μ , γ ) < 1 .
Thus, we can state that the self-mapping S has a unique fixed point. Also, S proves Lemma 3 with
σ 1 = 0 , η 1 = c 3 β H 1 ( λ , μ , η ) + d β 1 H 2 ( λ , μ , η ) + 1 3 β 2 ( γ 1 2 + γ 1 γ 2 + γ 2 2 ) H 3 ( λ , μ , η ) .
Therefore, S is in agreement with Picard’s S -stable. Hence the numerical method is stable.

5. Computational Results

To illustrate the numerical application of the hired scheme to investigate its efficiency and applicability, we consider the numerical results and solve them by applying the present method. In our implementation, the proposed method was carried out by applying Matlab. The absolute error function z z n is obtained by the L 2 n o r m , where z n is the numerical solution and z is the exact solution computed in the present method. To show the computational results, we consider the time-fractional Equation (1) with z 0 = 3 c 10 d s e c h 2 λ x , where λ = c 20 d . The exact solution of this equation for the special case η = 0 , μ = 1 , c > 0 and d < 0 is z ( x , t ) = 3 c 10 d s e c h 2 λ ( x c t ) , where c = 25 d 4 c 2 25 d . By decomposing the presented method on Equation (1), we have
n = 0 β n z n ( x , t ) = 3 c 10 d s e c h 2 λ x T 1 [ ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c ( n = 0 β n z n ( x , t ) ) x x x ( x , t ) + d ( n = 0 β n z n ( x , t ) ) x x x x x + n = 0 β n F n ( z ( x , t ) ) ; t ] ,
in which F n ( z ) is the nonlinear part defined by
n = 0 β n F n ( z ) = z 2 z x .
The polynomials F n ( z ) for some of the factors are calculated as
F 0 ( z ) = z 0 2 ( z 0 ) x , F 1 ( z ) = z 0 2 ( z 1 ) x + 2 z 0 z 1 ( z 0 ) x ,
We can compute p in Equation (44) as
β 0 : z 0 ( x , t ) = 3 c 10 d s e c h 2 λ x , β 1 : z 1 = T 1 ( 1 υ α μ ( 1 σ υ λ α ) η ) × T c ( z 0 ) x x x + d ( z 0 ( x , t ) ) x x x x x + F 0 ( z ) ; t = [ 24 λ 3 c 2 10 d s e c h 2 λ x tanh 3 λ x + 24 λ 3 c 2 10 d s e c h 4 λ x tanh λ x + 12 λ 3 c 2 10 d s e c h 2 λ x tanh λ x 96 λ 5 c d 10 d s e c h 2 λ x tanh 5 λ x + 192 λ 5 c d 10 d s e c h 4 λ x tanh 3 λ x + 168 λ 5 c d 10 d s e c h 2 λ x tanh 3 λ x 168 λ 5 c d 10 d s e c h 2 λ x tanh λ x + 384 λ 5 c d 10 d s e c h 5 λ x tanh 3 λ x 336 λ 5 c d 10 d s e c h 6 λ x tanh λ x + 27 c 3 5 d 10 d s e c h 6 λ x tanh λ x ] t μ 1 E λ , μ η ( σ t λ ) ,
Finally, the series numerical solution is calculated by
z = n = 0 z n = 3 c 10 d s e c h 2 λ x + [ 24 λ 3 c 2 10 d s e c h 2 λ x tanh 3 λ x + 24 λ 3 c 2 10 d s e c h 4 λ x tanh λ x + 12 λ 3 c 2 10 b s e c h 2 λ x tanh λ x 96 λ 5 c d 10 b s e c h 2 λ x tanh 5 λ x + 192 λ 5 c d 10 b s e c h 4 λ x tanh 3 λ x + 168 λ 5 c d 10 b s e c h 2 λ x tanh 3 λ x 168 λ 5 c d 10 b s e c h 2 λ x tanh λ x + 384 λ 5 c d 10 b s e c h 5 λ x tanh 3 λ x 336 λ 5 c d 10 b s e c h 6 λ x tanh λ x + 27 c 3 5 d 10 d s e c h 6 λ x tanh λ x ] t μ α E λ α , 1 μ α η ( σ t λ α ) + .
Now, we study a comparative examination of the absolute error function results for classical Equation (1), for the special case γ = 0 , μ = 1 . Table 1 presents the computational results given by the application of our proposed method and the other one obtained by the authors of [58,59] for the special case γ = 0 , μ = 1 . The computational results for different values of μ are demonstrated in Figure 1 and Figure 2 by selecting μ = 0.75 , 0.85 , 0.95 , and λ = 0.5 , η = 0.75 , σ = 1 . In particular, the computational results for the solutions corresponding to λ = 0.5 , η = 0.75 , σ = 1 at disparate values of μ with t = 0.5 are presented in Figure 3. Figure 4 shows the graphs of the absolute error function for the case η = 0 , μ = 1 . According to the established results, the achieved method is very accurate and effective in all cases.

6. Conclusions

In the present investigation, we studied the modified Kawahara equation associated with the Caputo–Prabhakar fractional operator. With the help of a unified computational technique with the homotopy perturbation algorithm and Sadik transform, the iterative solution was attained. To show and examine the stability, existence, and uniqueness of the solutions of the given suggested model, we used the fixed-point theory as well as the theory of S -stable mapping. Some stimulating numerical results were obtained at different values of the parameters μ , ρ , γ , and ω . The attained results illustrate that the Caputo–Prabhakar fractional operator can be applied to demonstrate real-world problems.

Author Contributions

Investigation, Methodology and Writing—original draft, Z.C., S.K., J.S. and M.H.D.; Software and Conceptualization, Z.C., S.K., J.S. and M.H.D.; Writing—review and editing, Z.C., S.K., J.S. and M.H.D.; visualization and supervision, Z.C., S.K., J.S. and M.H.D. All authors have read and agreed to the published version of the manuscript.

Funding

The work is supported by the Chinese National Natural Science Foundation under Grant 62272112, Central Science and Technology Development Project led by Hubei Province (2023EGA001), and the National key R&D Program of China (Grant 2019YFA0706338402).

Data Availability Statement

Data sharing is not applicable to this article, as no datasets were generated or analyzed during the current study.

Conflicts of Interest

Author Zhihua Chen was employed by the Wuhan Huamao Automation Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Amiataei, A.; Derakhshan, M. Stability analysis and numerical simulation of non-steady partial differential model in the human pulmonary capillaries using finite differences technique. Hacet. J. Math. Stat. 2023, 52, 1658–1676. [Google Scholar] [CrossRef]
  2. Diethelm, K.; Baleanu, D.; Scalas, E. Fractional Calculus: Models and Numerical Methods; World Scientific: Singapore, 2012. [Google Scholar]
  3. Goswami, A.; Singh, J.; Kumar, D. An efficient analytical approach for fractional equal width equations describing hydro-magnetic waves in cold plasma. Phys. A Stat. Mech. Appl. 2019, 524, 563–575. [Google Scholar] [CrossRef]
  4. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Generalized Mittag-Leffler function and generalized fractional calculus operators. Integral Transform. Spec. Funct. 2004, 15, 31–49. [Google Scholar] [CrossRef]
  5. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Academic Press: Cambridge, MA, USA; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  6. Farid, G.; Pecaric, J.; Tomovski, Z. Opial-type inequalities for fractional integral operator involving Mittag-Leffler function. Fract. Differ. Calc. 2015, 5, 93–106. [Google Scholar] [CrossRef]
  7. Kako, M. Nonlinear wave modulation in cold magnetized plasmas. J. Phys. Soc. Jpn. 1972, 33, 1678–1687. [Google Scholar] [CrossRef]
  8. Xiong, T.; Lv, Y.; Yi, W. Nonlinear Vibration and Control of Underwater Supercavitating Vehicles. IEEE Access 2018, 6, 62503–62513. [Google Scholar] [CrossRef]
  9. Yépez-Martínez, H.; Gómez-Aguilar, J.F. A new modified definition of Caputo-Fabrizio fractional-order derivative and their applications to the Multi Step Homotopy Analysis Method (MHAM). J. Comput. Appl. Math. 2019, 346, 247–260. [Google Scholar] [CrossRef]
  10. Kosari, S.; Xu, P.; Shafi, J.; Derakhshan, M.H. An efficient hybrid numerical approach for solving two-dimensional fractional cable model involving time-fractional operator of distributed order with error analysis. Numer. Algorithms 2024, 1–20. [Google Scholar] [CrossRef]
  11. Lu, D.; Seadawy, A.R.; Arshad, M. Bright-dark solitary wave and elliptic function solutions of unstable nonlinear Schrödinger equation and their applications. Opt. Quantum Electron. 2018, 50, 23. [Google Scholar] [CrossRef]
  12. Singh, J.; Kumar, D.; Qurashi, M.A.; Baleanu, D. Analysis of a new fractional model for damped Bergers’ equation. Open Phys. 2017, 5, 35–41. [Google Scholar] [CrossRef]
  13. Shchepetil’inikov, V.A.; Starodubtsev, P.A.; Traskovskii, V.I.; Senchenko, A.G.; Alifanov, R.N. Method of the finding the nonlinear phenomenas in ocean from moving undersea object and their theoretical explanation. J. Sib. Fed. Univ. Math. Phys. 2010, 3, 267–275. [Google Scholar]
  14. Derakhshan, M.H.; Kumar, P.; Salahshour, S. A high-order space-time spectral method for the distributed-order time-fractional telegraph equation. Int. J. Dyn. Control 2024, 12, 2778–2794. [Google Scholar] [CrossRef]
  15. Garra, R.; Garrappa, R. The Prabhakar or three parameter Mittag-Leffler function: Theory and application. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 314–329. [Google Scholar] [CrossRef]
  16. Garra, R.; Gorenflo, R.; Polito, F.; Tomovski, Ž. Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef]
  17. Kawahara, T. Oscillatory solitary waves in dispersive media. J. Phys. Soc. Jpn. 1972, 33, 260–264. [Google Scholar] [CrossRef]
  18. Amick, C.J.; McLeod, J.B. A singular perturbation problem in water waves. Stab. Appl. Anal. Cont. Media 1992, 1, 127–148. [Google Scholar]
  19. Amick, C.J.; Kirchgässner, K.A. Theory of solitary water-waves in the presence of surface tension. Arch. Ration. Mech. Anal. 1989, 105, 1–49. [Google Scholar] [CrossRef]
  20. Kurulay, M. Approximate analytic solutions of the modified Kawahara equation with homotopy analysis method. Adv. Differ. Equ. 2012, 2012, 178. [Google Scholar] [CrossRef]
  21. Kumar, D.; Singh, J.; Baleanu, D. Modified Kawahara equation within a fractional derivative with non-singular kernel. Therm. Sci. 2018, 22, 789–796. [Google Scholar] [CrossRef]
  22. Arora, R.; Siddiqui, M.J.; Singh, V.P. Solution of modified equal width wave equation, its variant and non-homogeneous Burgers’ equation by RDT method. Am. J. Comput. Appl. Math. 2011, 1, 53–56. [Google Scholar] [CrossRef]
  23. Gardner, L.R.; Gardner, G.A. Solitary waves of the equal width wave equation. J. Comput. Phys. 1992, 101, 218–223. [Google Scholar] [CrossRef]
  24. Raslan, K.R. Collocation method using quartic B-spline for the equal width (EW) equation. Appl. Math. Comput. 2005, 168, 795–805. [Google Scholar] [CrossRef]
  25. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  26. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  27. Miskinis, P. The Havriliak-Negami susceptibility as a nonlinear and nonlocal process. Phys. Scr. 2009, 2009, 014019. [Google Scholar] [CrossRef]
  28. Garrappa, R.; Mainardi, F.; Guido, M. Models of dielectric relaxation based on completely monotone functions. Fract. Calc. Appl. Anal. 2016, 19, 1105–1160. [Google Scholar] [CrossRef]
  29. Stanislavsky, A.; Weron, K. Atypical case of the dielectric relaxation responses and its fractional kinetic equation. Fract. Calc. Appl. Anal. 2016, 19, 212–228. [Google Scholar] [CrossRef]
  30. Giusti, A.; Colombaro, I. Prabhakar-like fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2018, 6, 38–143. [Google Scholar] [CrossRef]
  31. D’Ovidio, M.; Polito, F. Fractional Diffusion-Telegraph Equations and Their Associated Stochastic Solutions. Theory Probab. Appl. 2018, 62, 552–574. [Google Scholar] [CrossRef]
  32. An, J.; Hese, E.V.; Baes, M. Phase-space consistency of stellar dynamical models determined by separable augmented densities. Mon. Not. R. Astron. Soc. 2012, 422, 652–664. [Google Scholar] [CrossRef]
  33. Beghin, L.; Orsingher, E. Fractional Poisson processes and related planar random motions. Electron. J. Probab. 2009, 14, 1790–1826. [Google Scholar] [CrossRef]
  34. Górska, K.; Horzela, A.; Bratek, L.; Penson, K.A.; Dattoli, G. The probability density function for the Havriliak-Negami relaxation. arXiv 2016, arXiv:1611.06433. [Google Scholar]
  35. Shaikh, S.L. Introducing a new integral transform Sadik transform. Am. Int. J. Res. Sci. Technol. Eng. Math. 2018, 22, 100–102. [Google Scholar]
  36. Shaikh, S.L. Some applications of the new integral transform for Partial differential Equations. Math. J. Interdiscip. Sci. 2018, 7, 45–49. [Google Scholar] [CrossRef]
  37. Shaikh, S.L. Sadik transform in control theory. Int. J. Innov. Sci. Res. Technol. 2018, 3, 1–3. [Google Scholar]
  38. Jumarie, G. Laplace’s transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville derivative. Appl. Math. Lett. 2009, 22, 1659–1664. [Google Scholar] [CrossRef]
  39. Kexue, L.; Jigen, P. Laplace transform and fractional differential equations. Appl. Math. Lett. 2011, 24, 2019–2023. [Google Scholar] [CrossRef]
  40. Abhale, N.S.; Pawar, S.S. Fundamental properties of Sadik transform and it’s applications. J. Appl. Math. Comput. 2019, 6, 995–999. [Google Scholar]
  41. Lu, T.T.; Zheng, W.Q. Adomian decomposition method for first order PDEs with unprescribed data. Alex. Eng. J. 2020, 60, 2563–2572. [Google Scholar] [CrossRef]
  42. Li, W.; Pang, Y. Application of Adomian decomposition method to nonlinear systems. Adv. Differ. Equ. 2020, 2020, 67. [Google Scholar] [CrossRef]
  43. Zeidan, D.; Chau, C.K.; Lu, T.T.; Zheng, W.Q. Mathematical studies of the solution of Burgers’ equations by Adomian decomposition method. Math. Methods Appl. Sci. 2020, 43, 2171–2188. [Google Scholar] [CrossRef]
  44. Pasha, S.A.; Nawaz, Y.; Arif, M.S. The modified homotopy perturbation method with an auxiliary term for the nonlinear oscillator with discontinuity. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1363–1373. [Google Scholar] [CrossRef]
  45. Veeresha, P.; Prakasha, D.G.; Kumar, D. An efficient technique for nonlinear time-fractional Klein-Fock-Gordon equation. Appl. Math. Comput. 2020, 364, 124637. [Google Scholar] [CrossRef]
  46. Yu, D.N.; He, J.H.; Garcıa, A.G. Homotopy perturbation method with an auxiliary parameter for nonlinear oscillators. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1540–1554. [Google Scholar] [CrossRef]
  47. Khalid, M.; Khan, F.S.; Sultana, M. A highly accurate numerical method for solving nonlinear time-fractional differential difference equation. Math. Methods Appl. Sci. 2019, 44, 8243–8253. [Google Scholar] [CrossRef]
  48. Bhowmik, S.K.; Karakoc, S.B. Numerical solutions of the generalized equal width wave equation using the Petrov-Galerkin method. Appl. Anal. 2019, 21, 714–734. [Google Scholar] [CrossRef]
  49. GaziKarakoc, S.B.; Ali, K.K. Analytical and computational approaches on solitary wave solutions of the generalized equal width equation. Appl. Math. Comput. 2020, 71, 124933. [Google Scholar] [CrossRef]
  50. Anjum, N.; He, J.H. Laplace transform: Making the variational iteration method easier. Appl. Math. Lett. 2019, 1, 134–138. [Google Scholar] [CrossRef]
  51. Alomari, A.K.; Syam, M.I.; Anakira, N.R.; Jameel, A.F. Homotopy Sumudu transform method for solving applications in physics. Results Phys. 2020, 18, 103265. [Google Scholar] [CrossRef]
  52. Eshaghi, J.; Kazem, S.; Adibi, H. The local discontinuous Galerkin method for 2D nonlinear time-fractional advection-diffusion equations. Eng. Comput. 2019, 35, 1317–1332. [Google Scholar] [CrossRef]
  53. Goyal, M.; Prakash, A.; Gupta, S. An Efficient Perturbation Sumudu Transform Technique for the Time-Fractional Vibration Equation with a Memory Dependent Fractional Derivative in Liouville–Caputo Sense. Int. J. Appl. Comput. Math. 2021, 7, 156. [Google Scholar] [CrossRef]
  54. Prakash, A.; Goyal, M.; Gupta, S. Fractional variational iteration method for solving time-fractional Newell-Whitehead-Segel equation. Nonlinear Eng. 2019, 8, 164–171. [Google Scholar] [CrossRef]
  55. Prakash, A.; Kaur, H. A New Numerical Method for a Fractional Model of Non-Linear Zakharov–Kuznetsov Equations via Sumudu Transform. In Methods of Mathematical Modelling: Fractional Differential Equations; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  56. Redhwan, S.S.; Shaikh, S.L.; Abdo, M.S. Some properties of Sadik transform and its applications of fractional-order dynamical systems in control theory. Adv. Theory Nonlinear Anal. Its Appl. 2019, 4, 51–66. [Google Scholar] [CrossRef]
  57. Ghorbani, A. Beyond Adomian polynomials: He polynomials. Chaos Solitons Fractals 2009, 39, 1486–1492. [Google Scholar] [CrossRef]
  58. Bhatter, S.; Mathur, A.; Kumar, D.; Nisar, K.S.; Singh, J. Fractional modified Kawahara equation with Mittag-Leffler law. Chaos Solitons Fractals 2020, 131, 109508. [Google Scholar] [CrossRef]
  59. Jin, L. Application of variational iteration method and homotopy perturbation method to the modified Kawahara equation. Math. Comput. Model. 2009, 49, 573–578. [Google Scholar] [CrossRef]
Figure 1. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with c = 1 , d = 1 .
Figure 1. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with c = 1 , d = 1 .
Fractalfract 08 00512 g001
Figure 2. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with c = 1 , d = 1 .
Figure 2. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with c = 1 , d = 1 .
Fractalfract 08 00512 g002
Figure 3. The absolute error function graph for η = 0 and μ = 1 when t = 0.5 with c = 1 , d = 1 .
Figure 3. The absolute error function graph for η = 0 and μ = 1 when t = 0.5 with c = 1 , d = 1 .
Fractalfract 08 00512 g003
Figure 4. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with t = 0.5 , c = 1 , d = 1 .
Figure 4. The approximate solution graph corresponding to λ = 0.5 , η = 0.75 , σ = 1 at different values of μ with t = 0.5 , c = 1 , d = 1 .
Fractalfract 08 00512 g004
Table 1. The exact solution values when μ = 1 and comparison between the proposed method with the methods VIM [59], HPM [59], and HATM [58] corresponding to ρ = 0.5 , γ = 0.75 , ω = 1 and μ = 0.99 , α = 0.5 , a = 0.001 , b = 1 .
Table 1. The exact solution values when μ = 1 and comparison between the proposed method with the methods VIM [59], HPM [59], and HATM [58] corresponding to ρ = 0.5 , γ = 0.75 , ω = 1 and μ = 0.99 , α = 0.5 , a = 0.001 , b = 1 .
xtExact SolutionVIM [59]HPM [59]HATM [58]
5 0.2 9.474889415 × 10 7 9.474984315 × 10 4 9.474984315 × 10 4 9.474984315 × 10 4
2.5 0.4 9.483773375 × 10 6 9.483868961 × 10 4 9.483868961 × 10 4 9.483868965 × 10 4
0 0.6 9.486831272 × 10 7 9.486832980 × 10 4 9.486832980 × 10 4 9.486832980 × 10 4
2.5 0.8 9.484055589 × 10 6 9.483868961 × 10 4 9.483868961 × 10 4 9.483868965 × 10 4
5 1.0 9.475453144 × 10 6 9.474984315 × 10 4 9.474984315 × 10 4 9.474984315 × 10 4
xtOur Method
5 0.2 9.745606079 × 10 7
2.5 0.4 1.526682172 × 10 6
0 0.6 9.720373521 × 10 7
2.5 0.8 1.076279876 × 10 6
5 1.0 2.221689661 × 10 6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Z.; Kosari, S.; Shafi, J.; Derakhshan, M.H. Stability Analysis Study of Time-Fractional Nonlinear Modified Kawahara Equation Based on the Homotopy Perturbation Sadik Transform. Fractal Fract. 2024, 8, 512. https://doi.org/10.3390/fractalfract8090512

AMA Style

Chen Z, Kosari S, Shafi J, Derakhshan MH. Stability Analysis Study of Time-Fractional Nonlinear Modified Kawahara Equation Based on the Homotopy Perturbation Sadik Transform. Fractal and Fractional. 2024; 8(9):512. https://doi.org/10.3390/fractalfract8090512

Chicago/Turabian Style

Chen, Zhihua, Saeed Kosari, Jana Shafi, and Mohammad Hossein Derakhshan. 2024. "Stability Analysis Study of Time-Fractional Nonlinear Modified Kawahara Equation Based on the Homotopy Perturbation Sadik Transform" Fractal and Fractional 8, no. 9: 512. https://doi.org/10.3390/fractalfract8090512

APA Style

Chen, Z., Kosari, S., Shafi, J., & Derakhshan, M. H. (2024). Stability Analysis Study of Time-Fractional Nonlinear Modified Kawahara Equation Based on the Homotopy Perturbation Sadik Transform. Fractal and Fractional, 8(9), 512. https://doi.org/10.3390/fractalfract8090512

Article Metrics

Back to TopTop