Next Article in Journal
A General Design Method of Cam Profile Based on Cubic Splines and Dynamic Model: Case Study of a Gravity-Driven Tricycle
Previous Article in Journal
Are Citizens Credentialist or Post-Nationalists? A Fuzzy-Eco Apostle Model Applied to National Identity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximate Solution of Nonlinear Time-Fractional PDEs by Laplace Residual Power Series Method

1
Department of Mathematical Sciences, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
2
Department of Mathematics, Faculty of Arts and Science, Amman Arab University, Amman 11953, Jordan
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(12), 1980; https://doi.org/10.3390/math10121980
Submission received: 9 May 2022 / Revised: 31 May 2022 / Accepted: 6 June 2022 / Published: 8 June 2022

Abstract

:
Most physical phenomena are formulated in the form of non-linear fractional partial differential equations to better understand the complexity of these phenomena. This article introduces a recent attractive analytic-numeric approach to investigate the approximate solutions for nonlinear time fractional partial differential equations by means of coupling the Laplace transform operator and the fractional Taylor’s formula. The validity and the applicability of the used method are illustrated via solving nonlinear time-fractional Kolmogorov and Rosenau–Hyman models with appropriate initial data. The approximate series solutions for both models are produced in a rapid convergence McLaurin series based upon the limit of the concept with fewer computations and more accuracy. Graphs in two and three dimensions are drawn to detect the effect of time-Caputo fractional derivatives on the behavior of the obtained results to the aforementioned models. Comparative results point out a more accurate approximation of the proposed method compared with existing methods such as the variational iteration method and the homotopy perturbation method. The obtained outcomes revealed that the proposed approach is a simple, applicable, and convenient scheme for solving and understanding a variety of non-linear physical models.

1. Introduction

The subject of fractional calculus is not new. It is a generalization of classical calculus that deals with the ordinary differentiation and integration of arbitrary order. It goes back to Leibniz in a letter to L’Hospital in the late seventeenth century. The main idea of fractional calculus is that natural phenomena modeling is not in integer operators; it is in fractional operators. So, the fractional calculus focuses on behaviors that cannot be modeled by traditional theory [1,2,3,4,5,6,7]. In the past, a lot of prominent contributions were made to the subject of the theory and applications of fractional partial differential equations (FPDEs). These equations are more effectively used to analyze and describe several phenomena in various fields such as mechanical systems, dynamical systems, control theory, mixed convection flows, heat transfer, unification of diffusion, image processing, and wave propagation phenomenon [8,9,10,11,12,13,14,15]. Nevertheless, no method gives an explicit solution for FPDEs due to the intricacies of the fractional calculus that includes these equations. Recently, numerical techniques such as the Adomian decomposition method (ADM), variational iteration method (VIM), reproducing kernel method (RKM), Laplace variational iteration method (LVIM), Laplace Adomian decomposition method (LADM), and residual power series method (RPSM) are used widely to find approximate solutions of many nonlinear fractional differential equations that do not have exact analytic solutions. For more information regarding the methods and numerical techniques for solving fractional differential equations [16,17,18,19,20,21,22]. On the other hand, RPSM has been widely used to find out the solutions to linear and nonlinear issues of fractional differential, and it is used to find out the solution for the system of FPDEs [22]. Additionally, it is used for well-known partial differential equations of fractional order, such as fractional Newell–Whitehead–Segel equation [23], time-fractional Fokker–Planck equations [24], fractional Kundu–Eckhaus and massive Thirring models [25], coupled fractional resonant Schrödinger equation [26], and fuzzy fractional IVPs [27,28,29,30,31,32]. The proposed algorithm is straightforward, accurate, and powerful and creates a series of solutions for different models that occur in applied mathematics without terms of perturbation, discretization, and linearization.
Creating approximation solutions for nonlinear time FPDEs using the aforementioned numeric-analytic methods and others is a significant matter for scholars. Thus, there has become an insistent requirement for efficient semi-analytic methods to construct precise solutions for both linear and nonlinear fractional models. Motivated by this, the primary contribution of this work is to create accurate approximate solutions in a closed-form series for a certain class of nonlinear time FPDEs in light of the time-Caputo fractional derivative sense via extending the application of the Laplace RPSM. This method is proposed and proved by El-Ajou [31] to investigate the exact solitary solutions for a class of nonlinear time-FPDEs. It depends basically on treating the main problem in Laplace space with the help of RPSM, where the unknown coefficients could be found via the concept limit, unlike the RPSM which uses the fractional derivatives in each step to find these coefficients [33]. The proposed method has been successfully employed to produce exact and precise approximate solutions by involving fast convergent power series for emerging realism models in physical phenomena due to its features, which are that it is easy, straightforward, handles directly to various kinds of initial conditions, needs no to linearization or restrictive assumptions, does not need major computational requirements and is performed with less time and more accuracy. More applications, analysis, and advanced techniques used to process and solve linear and non-linear real-life models are found in the references [34,35,36,37,38,39,40,41,42,43,44,45,46,47].
The structure of the article is arranged as follows. In Section 2, essential definitions, properties, and theorems about fractional calculus, Laplace transform, and Laplace fractional expansion (FE) are shown. The methodology of Laplace RPSM for solving nonlinear time-FPDEs is investigated in Section 3. In Section 4, two initial value problems (IVPs) of fractional-order Kolmogorov equation and Rosenau–Hyman equation are solved to show the applicability and accuracy of our approach. Finally, Section 5 is devoted to the conclusions.

2. Basic Concepts and Notations

In this section, we review the essential definitions and theorems of fractional derivatives in the sense of Caputo. Additionally, we revise the primary definitions and theorems related to Laplace transform which will be used mainly in the next section.
Definition 1
(See Ref. [3]). For 𝒶 + , the Riemann–Liouville fractional integral operator for a real-valued function 𝒲 ( x , t ) is denoted by 𝓘 t 𝒶 and defined as:
𝓘 t 𝒶 W ( x , t ) = { 1 Γ ( 𝒶 ) 0 t W ( x , η ) ( t η ) 1 𝒶 d η ,     0 η < t ,   𝒶 > 0 , W ( x , t ) , 𝒶 = 0 .
Definition 2
(See Ref. [3]). The time fractional derivative of order 𝒶 > 0 , for the function W ( x , t ) in the Caputo case is denoted by D t 𝒶 , and defined as:
D t 𝒶 W ( x , t ) = { 𝓘 t n α ( D t n W ( x , t ) ) , 0 < n 1 < 𝒶 n ,   D t n W ( x , t ) , 𝒶 = n ,
where D t n = n t n , and n .
Consequently, for n 1 < 𝒶 n ,   β > 1 and t 0 , the operators D t 𝒶 and 𝓘 t α satisfy the following properties:
  • D t 𝒶 c = 0 , c .
  • D t 𝒶 t β = Γ ( β + 1 ) Γ ( β + 1 𝒶 ) t β 𝒶 .
  • D t 𝒶 𝓘 t 𝒶 W ( x , t ) = W ( x , t ) .
  • 𝓘 t 𝒶 D t 𝒶 W ( x , t ) = W ( x , t ) j = 0 n 1 D t j ( x , 0 + ) t j j !   , for W C n [ a , b ] , n 1 < 𝒶 n ,   n   a n d   a , b .
Definition 3
(See Ref. [31]). Let W ( x , t ) is a piecewise continuous function on I × [ 0 , ) and of exponential order δ . Then, the Laplace transformation of the function W ( x , t ) is denoted and defined as follows:
𝝕 ( x , s ) = [ W ( x , t ) ] 0 e s t W ( x , t ) d t , s > δ ,
whereas the inverse Laplace transformation of the function W ( x , s ) is defined as follows:
W ( x , t ) = 1 [ 𝝕 ( x , s ) ] c i c + i e s t 𝝕 ( x , s ) d s ,       c = R e ( s ) > δ 0 ,
where δ 0 lies in the right half plane of the absolute convergence of the Laplace integral.
Lemma 1
(See Ref. [31]). Let W ( x , t ) and φ ( x , t ) are piecewise continuous functions on I × [ 0 , ) and of exponential order δ 1 and δ 2 , respectively, where δ 1 < δ 2 . Suppose that 𝝕 ( x , s ) = [ W ( x , t ) ] ,   Φ ( x , s ) = [ φ ( x , t ) ] , and a ,   b are constants. Then, the following properties are satisfied:
  • [ a W ( x , t ) + b φ ( x , t ) ] = a 𝝕 ( x , s ) + b Φ ( x , s ) , x I ,   s > δ 1 .
  • 1 [ a 𝝕 ( x , s ) + b Φ ( x , s ) ] = a W ( x , t ) + b φ ( x , t ) , x I ,   t 0 .
  • [ e a t W ( x , t ) ] = 𝝕 ( x , s a ) , x I ,   s > a + δ 1 .
  • lim s s 𝝕 ( x , s ) = W ( x , 0 ) ,   x I .
Lemma 2
(See Ref. [31]). Let W ( x , t ) be a piecewise continuous function on I × [ 0 , ) and of exponential order δ , and 𝝕 ( x , s ) = [ W ( x , t ) ] . Then,
  • [ 𝓘 W ( x , t ) ] = s 𝒶 1 𝝕 ( x , s ) ,   𝒶 > 0
  • [ D t 𝒶 W ( x , t )   ] = s 𝒶 𝝕 ( x , s ) k = 0 n 1 s 𝒶 k 1 D t k W ( x , 0 ) , n 1 < 𝒶 n
  • [ D t j 𝒶 W ( x , t )   ] = s j 𝒶 𝝕 ( x , s ) k = 0 j 1 s ( j k ) 𝒶 1 D t k 𝒶 W ( x , 0 ) , 0 < 𝒶 1 , where  D t j 𝒶 = D t 𝒶 · D t 𝒶 D t 𝒶 ( j t i m e s )
Theorem 1
(See Ref. [31]). Let W ( x , t ) be a piecewise continuous function on I × [ 0 , ) and of exponential order δ . Suppose that the function 𝝕 ( x , s ) = [ W ( x , t ) ] has the following fractional expansion:
𝝕 ( x , s ) = n = 0 h n ( x ) s n α + 1 ,   x I ,   s > δ ,   0 < 𝒶 1 .
Then, h n ( x ) = D t n 𝒶 W ( x , 0 ) .
Remark 1.
The inverseLaplace transformation  W ( x , t ) = 1 [ 𝝕 ( x , s ) ] in Theorem 1 is in the following form:
W ( x , t ) = n = 0 D t n 𝒶 W ( x , 0 ) t n 𝒶 Γ ( n 𝒶 + 1 ) ,   0 <   α 1 , t > 0 .
Theorem 2
(See Ref. [31]). Let W ( x , t ) be a piecewise continuous function on I × [ 0 , ) and of exponential order δ and 𝝕 ( x , s ) = [ W ( x , t ) ] can be represented as the fractional expansion in Theorem 1. If | s [ D t ( n + 1 ) a W ( x , t )   ] | M ( x ) , on I × ( δ , γ ] where 0 < 𝒶 1 , then the remainder R n ( x , s ) of the FE in Theorem 1 satisfies the following inequality:
| R n ( x , s ) | M ( x ) s 1 + ( n + 1 ) 𝒶 ,   x I ,   δ < s γ .
Theorem 3.
If 𝒶 ( 0 , 1 ) , W k + 1 ( x , t ) 𝒶 W k ( x , t ) gives k N and 0 < t < T < 1 , then the series of numerical solutions converges to an exact solution [35].
Proof. 
We notice that   0 < t < T < 1 ,
W ( x , t ) W k ( x , t ) = m = k + 1 W m ( x , t ) m = k + 1 W m ( x , t ) ( η ) m = k + 1 C m = 𝒶 k + 1 1 𝒶 ( η ) 0   a s   k .
 □

3. Methodology of Laplace RPSM

In this section, we clarify the principle of the Laplace RPSM algorithm to solve nonlinear time fractional PDEs. Our strategy to use the proposed scheme depends on coupling the Laplace transform operator and fractional RPSM. More specifically, consider the following initial value problem of nonlinear time fractional PDEs:
{ D t 𝒶 W ( x , t ) = N x [ W ( x , t ) ] , W ( x , 0 ) = f ( x ) ,
where N x is a nonlinear operator relative to x of degree r , x I , t 0 , D t 𝒶 refers to 𝒶 -th Caputo fractional derivative for 𝒶 ( 0 , 1 ] , and W ( x , t ) is an unknown function to be determined.
To construct the approximate solution of (1) by using the Laplace RPSM, one can perform the following procedure:
Step A: Apply the Laplace transform on both sides of (1), and utilizing the initial data of (1), as well as depending on Lemma 2, part (2), we get:
𝝕 ( x , s ) = f ( x ) s 1 s 𝒶 { N x [ W ( x , t ) ] } , where   𝝕 ( x , s ) = [ W ( x , t ) ] ( s ) ,   s > δ .
Step B: According to Theorem 1, we assume that the approximate solution of the Laplace Equation (2) takes the following fractional expansion:
𝝕 ( x , s ) = f ( x ) s + n = 1 h n ( x ) s n 𝒶 + 1 ,   x I ,   s > δ 0 ,
and the k-th Laplace series solution takes the following form:
𝝕 k ( x , s ) = f ( x ) s + n = 1 k h n ( x ) s n 𝒶 + 1 ,   x I ,   s > δ 0 .
Step C: We define the k-th Laplace fractional residual function of (2) as
( R e s 𝝕 k ( x , s ) ) = 𝝕 k ( x , s ) f ( x ) s + 1 s 𝒶 { N x [ W ( x , t ) ] } ,
and the Laplace residual function of (2) are defined as:
lim k ( R e s 𝝕 k ( x , s ) ) = ( R e s 𝝕 ( x , s ) ) = 𝝕 ( x , s ) f ( x ) s + 1 s 𝒶 { N x [ W ( x , t ) ] } .
As in [31,32,33], some useful facts of Laplace residual function which are essential in finding the approximate solution are listed as follows:
  • lim k ( R e s 𝝕 k ( x , s ) ) = ( R e s 𝝕 ( x , s ) ) , for x I ,   s > δ 0 .
  • ( R e s 𝝕 ( x , s ) ) = 0 , for x I ,   s > δ 0 .
  • lim s s k 𝒶 + 1 ( R e s 𝝕 k ( x , s ) ) = 0 , for x I ,   s > δ 0 , and k = 1 , 2 , 3 ,
Step D: Substitute the k-th Laplace series solution (4) into the k-th Laplace fractional residual function of (5).
Step E: The unknown coefficients h k ( x ) , for k = 1 , 2 , 3 , , could be founded by solving the system lim s s k 𝒶 + 1 ( R e s 𝝕 k ( x , s ) ) = 0 . Then, we collect the obtained coefficients in terms of fractional expansion series (4) 𝝕 k ( x , s ) .
Step F: Running the inverse Laplace transform operator on both sides of the obtained Laplace series solution to get the approximate solution W k ( x , t ) , of the main Equation (1).

4. Numerical Examples

In this section, the superiority, efficiency, and applicability of the Laplace RPSM are demonstrated by testing two non-linear time fraction IVPs. It is worth mention here that all numerical computations and symbolic have been carried out using MATHEMATICA 12 software package.
Example 1.
Consider the following nonlinear time-fractional Kolmogorov IVP:
{ D t 𝒶 W ( x , t ) ( x + 1 ) D x W ( x , t ) x 2 e t D x 2 W ( x , t ) = 0 , W ( x , 0 ) = x + 1 ,
where 0 < 𝒶 1 , and ( x , t ) [ 0 , 1 ] × . The exact solutions for standard case 𝒶 = 1 , is given as W ( x , t ) = ( x + 1 ) e t .
By applying the Laplace transform operator on the both sides of time-fractional Kolmogorov equation of (7) and using part 2 in Lemma 2 and the initial data-space of ( 7 ) , we get the following Laplace fractional equation:
𝝕 ( x , s ) = x + 1 s + x + 1 s 𝒶 D x 𝝕 ( x , s ) + x 2 s 𝒶 { 1 { 1 s 1 } D x 2 1 { 𝝕 } } ,
where 𝝕 ( x , s ) = [ W ( x , t ) ] .
Considering the last discussion, the k-th Laplace series, 𝝕 k ( x , s ) for (8) is expressed as the form:
𝝕 k ( x , s ) = x + 1 s + n = 1 k h n ( x ) s n 𝒶 + 1 .
In addition, we define the k-th Laplace residual function of (8) as
( R e s 𝝕 k ( x , s ) ) = n = 1 k h n ( x ) s n 𝒶 + 1 x + 1 s 𝒶 D x ( x + 1 s + n = 1 k h n ( x ) s n 𝒶 + 1 ) x 2 s 𝒶 { 1 { 1 s 1 } D x 2 1 { x + 1 s + n = 1 k h n ( x ) s n 𝒶 + 1 } } .
For k = 1 , in (10), we get the 1-st Laplace residual function as
( R e s 𝝕 1 ( x , s ) ) = h 1 ( x ) s 𝒶 + 1 x + 1 s 𝒶 D x ( x + 1 s + h 1 ( x ) s 𝒶 + 1 ) x 2 s 𝒶 { 1 { 1 s 1 } D x 2 1 { x + 1 s + h 1 ( x ) s 𝒶 + 1 } } = h 1 ( x ) s 𝒶 + 1 x + 1 s 𝒶 ( 1 s + h 1 ( x ) s 𝒶 + 1 ) x 2 h 1 ( x ) ( s 1 ) 2 𝒶 + 1 s 𝒶 .
Next, multiply both sides of (11) by s 𝒶 + 1 , and then solve the system l i m s s 𝒶 + 1 ( R e s 𝝕 1 ( x , s ) ) = 0 , and we get: h 1 ( x ) = x + 1 . So, the 1-st Laplace series solution of (8) could be written as: 𝝕 1 ( x , s ) = x + 1 s + x + 1 s 𝒶 + 1 .
For k = 2 , in (10), then the 2-nd Laplace residual function can be expressed as
( R e s 𝝕 2 ( x , s ) ) = x + 1 s 𝒶 + 1 + h 2 ( x ) s 2 𝒶 + 1 x + 1 s 𝒶 D x ( x + 1 s + x + 1 s 𝒶 + 1 + h 2 ( x ) s 2 𝒶 + 1 ) x 2 s 𝒶 { 1 { 1 s 1 } D x 2 1 { x + 1 s + x + 1 s 𝒶 + 1 + h 2 ( x ) s 2 𝒶 + 1 } } = x + 1 s 𝒶 + 1 + h 2 ( x ) s 2 𝒶 + 1 x + 1 s 𝒶 ( 1 s + 1 s 𝒶 + 1 + h 2 ( x ) s 2 𝒶 + 1 ) x 2 h 2 ( x ) ( s 1 ) 2 𝒶 + 1 s 𝒶
After that, we multiply the result of Equation (12) by the factor s 2 𝒶 + 1 to get the following equation:
s 2 𝒶 + 1 ( R e s 𝝕 2 ( x , s ) ) = h 2 ( x ) x 1 ( x + 1 ) h 2 ( x ) s 𝒶 x 2 h 2 ( x ) ( s 1 ) 2 𝒶 + 1 s 1 𝒶 .
By solving l i m s s 2 𝒶 + 1 ( R e s 𝝕 2 ( x , s ) ) = 0 , yields that:   h 2 ( x ) = x + 1 . So, the 2-nd Laplace series solution of (8) could be written as: 𝝕 2 ( x , s ) = x + 1 s + x + 1 s 𝒶 + 1 + x + 1 s 2 𝒶 + 1 .
Similarly, for k = 3 , we have
s 3 𝒶 + 1 ( R e s 𝝕 3 ( x , s ) ) = s 3 𝒶 + 1 ( x + 1 s 𝒶 + 1 + x + 1 s 2 𝒶 + 1 + h 3 ( x ) s 3 𝒶 + 1 x + 1 s 𝒶 D x ( x + 1 s + x + 1 s 𝒶 + 1 + x + 1 s 2 𝒶 + 1 + h 3 ( x ) s 3 𝒶 + 1 ) x 2 s 𝒶 { 1 { 1 s 1 } D x 2 1 { x + 1 s + x + 1 s 𝒶 + 1 + x + 1 s 2 𝒶 + 1 + h 3 ( x ) s 3 𝒶 + 1 } } ) = h 3 ( x ) ( x + 1 ) ( x + 1 ) h 2 ( x ) s 𝒶 x 2 h 3 ( x ) ( s 1 ) 3 a + 1 s 1 2 a ,
and by solving l i m s s 3 𝒶 + 1 ( R e s 𝝕 3 ( x , s ) ) = 0 , one can obtain that h 3 ( x ) = x + 1 . So, the 3 r d Laplace series solution of (8) could be written as: 𝝕 3 ( x , s ) = x + 1 s + x + 1 s 𝒶 + 1 + x + 1 s 2 𝒶 + 1 + x + 1 s 3 𝒶 + 1 .
Processing the previous steps for an arbitrary k , and using the fact l i m s s k 𝒶 + 1 ( R e s 𝝕 k ( x , s ) ) = 0 , one can obtain that: h k ( x ) = x + 1 ,   f o r   k = 4 , 5 , . Thus, the k-th Laplace series solution of (8) could be reformulated by the following fractional expansion:
𝝕 k ( x , s ) = ( 1 s + 1 s 𝒶 + 1 + 1 s 2 𝒶 + 1 + 1 s 3 𝒶 + 1 + + 1 s k 𝒶 + 1 ) = n = 0 k x + 1 s n 𝒶 + 1 .
Lastly, we apply the inverse Laplace transform for the obtained expansion (15) to conclude that the k-th approximate solution of the nonlinear time-fractional Kolmogorov IVP (7) has the form:
W k ( x , t ) = ( x + 1 ) n = 0 k t n 𝒶 Γ ( n 𝒶 + 1 ) ,
when k and by substituting 𝒶 = 1 in (16), we get the Maclaurin series expansion of the closed form W ( x , t ) = ( x + 1 ) e t , which is fully in agreement with the exact solution.
Numerical results of the 10-th approximate solutions for the nonlinear time-fractional Kolmogorov IVP (7) are computed and summarized in Table 1 at fixed values of the variable x , and some selected grid points in [ 0 , 1 ] with step size of 0.25 , and different values of fractional order 𝒶 ’s such that 𝒶 { 1 , 0.95 , 0.85 , 0.75 ,   0.65 } . From the table, it can be found that the present method provides us with an accurate approximate solution, which is in good agreement with each other for all values of t in [ 0 , 1 ] , especially when approaching the initial values. Further, numerical comparisons are performed in Table 2 to validate the accuracy of our approach by establishing the recurrence errors | W 8 ( x , t ) W 7 ( x , t ) | for the obtained approximate solution of IVP (7) at various values of 𝒶 .
Figure 1 shows the graphs of the exact and the tenth approximate curves solutions for the nonlinear time-fractional Kolmogorov IVP (7) at various 𝒶 values. Obviously, one can see that the obtained approximate solutions for different values of fractional order simulate the solution for the classical case. Additionally, the exact and approximate solutions match at 𝒶 = 1 , and this confirms the effectiveness and performance of our approach. While Figure 2 demonstrates the comparison of the geometric behavior between the exact solution and the obtained 10-th approximate solution to the nonlinear time-fractional Kolmogorov IVP (7) at various 𝒶 values for ( x , t ) [ 0 , 1 ] 2 . From these 3D surface plots, we see that the solution behaviors for different Caputo fractional derivatives on their domain are in close agreement with each other, particularly for classical derivative. Moreover, the total calculation cost comparison of the given figures in Example 1 is reported in Table 3.
Example 2.
Consider the following nonlinear time-fractional Rosenau–Hyman IVP:
{ D t 𝒶 W W D x 3 W W D x W 3 D x W D x 2 W = 0 , W ( x , 0 ) = 8 c 3 cos 2 x 4 ,
where  0 < 𝒶 1 , and  ( x , t ) [ 0 , 1 ] × . The exact solutions for standard case  𝒶 = 1 , is  W ( x , t ) = 8 c 3 cos 2 ( x c t 4 ) .
According to the Laplace RPSM, we firstly operate the Laplace transform to nonlinear time-fractional Rosenau–Hyman of (17), and using the initial data-space of (17), we get:
𝝕 ( x , s ) = 8 c 3 s cos 2 x 4 + 1 s 𝒶 { 1 { 𝝕 } D x 3 1 { 𝝕 } } + 1 s 𝒶 { 1 { 𝝕 } D x 1 { 𝝕 } } + 3 s 𝒶 { D x 1 { 𝝕 } D x 2 1 { 𝝕 } } .
The k-th proposed Laplace series solution of the Laplace Equation (18) will be in the form:
𝝕 k ( x , s ) = 8 c 3 s cos 2 x 4 + n = 1 k h n ( x ) s n 𝒶 + 1 .
Therefore, the k-th Laplace residual function of (18) can be defined as
( R e s 𝝕 k ( x , s ) ) = 𝝕 k ( x , s ) + 8 c 3 s cos 2 x 4 1 s 𝒶 { 1 { 𝝕 k } D x 3 1 { 𝝕 k } } 1 s 𝒶 { 1 { 𝝕 k } D x 1 { 𝝕 k } } 3 s 𝒶 { D x 1 { 𝝕 k } D x 2 1 { 𝝕 k } } .
To define h 1 ( x ) , we consider k = 1 in (20)
( R e s 𝝕 1 ( x , s ) ) = 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 + 8 c 3 s cos 2 x 4 1 s 𝒶 { 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } D x 3 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } } 1 s 𝒶 { 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } D x 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } } 3 s 𝒶 { D x 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } D x 2 1 { 8 c 3 s cos 2 x 4 + h 1 ( x ) s 𝒶 + 1 } } .
Now, we multiply both sides of (21) by s 𝒶 + 1 to get:
s 𝒶 + 1 ( R e s 𝝕 1 ( x , s ) ) = 8 3 c 2 s i n ( x 4 ) c o s 3 ( x 4 ) + 32 3 c 2 sin ( x 4 ) c o s ( x 4 ) ( 1 8 s i n 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) + 3 c 8 s 𝒶 c o s 2 ( x 4 ) h 1 ( x ) + 8 c h 1 ( x ) s 𝒶 ( 1 8 sin 2 ( x 4 ) 1 8 c o s 2 ( x 4 ) ) 4 c s 𝒶 sin ( x 4 ) cos ( x 4 ) h 1 ( x ) + 8 c h 1 ( 3 ) ( x ) 3 s 𝒶 c o s 2 ( x 4 ) c h 1 ( x ) s 𝒶 s i n ( x 4 ) cos ( x 4 ) 3 Γ ( 2 𝒶 + 1 ) h 1 ( x ) h 1 ( x ) Γ 2 ( 𝒶 + 1 ) s 2 𝒶 Γ ( 2 𝒶 + 1 ) h 1 ( x ) h 1 ( x ) Γ 2 ( 𝒶 + 1 ) s 2 𝒶 Γ ( 2 𝒶 + 1 ) h 1 ( x ) h 1 ( 3 ) ( x ) Γ 2 ( 𝒶 + 1 ) s 2 𝒶 + h 1 ( x ) .
Next, by solving lim s s 𝒶 + 1 ( R e s 𝝕 1 ( x , s ) ) = 0 , and after some algebra simplification, one can get h 1 ( x ) = 2 c 2 3 sin x 2 . So, the 1-st Laplace series solution of (17) could be written as:
𝝕 1 ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2
To find h 2 ( x ) , we consider that k = 2 , in the Laplace residual Equation (20), and by multiplying the obtained equation by the factor s 2 𝒶 + 1 , we get
s 2 𝒶 + 1 ( R e s 𝝕 2 ( x , s ) ) = h 2 ( x ) + 64 3 c 2 cos ( x 4 ) sin ( x 4 ) ( 1 8 sin 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) + 2 3 c 3 cos ( x 4 ) sin ( x 4 ) sin ( x 4 ) + c sin ( x 4 ) cos ( x 4 ) Γ ( 5 𝒶 + 2 ) h 2 ( x ) h 2 ( 3 ) ( x ) 3 Γ ( 2 𝒶 + 1 ) Γ ( 3 𝒶 + 2 ) s 4 𝒶 + 1 4 c sin ( x 4 ) cos ( x 4 ) Γ ( 5 𝒶 + 2 ) h 2 ( x ) h 2 ( x ) 3 Γ ( 2 𝒶 + 1 ) Γ ( 3 𝒶 + 2 ) s 4 𝒶 + 1 2 c 3 sin ( x 4 ) cos ( x 4 ) sin ( x 2 ) Γ ( 4 𝒶 + 2 ) h 2 ( 3 ) ( x ) 9 Γ ( 𝒶 + 1 ) Γ ( 3 𝒶 + 2 ) s 3 𝒶 + 1 + 8 c 3 sin ( x 4 ) cos ( x 4 ) sin ( x 2 ) Γ ( 4 𝒶 + 2 ) h 2 ( x ) 9 Γ ( 𝒶 + 1 ) Γ ( 3 𝒶 + 2 ) s 3 𝒶 + 1 8 9 c 2 cos 3 ( x 4 ) sin ( x 4 ) h 2 ( 3 ) ( x ) s 2 𝒶 + 1 + 32 9 c 2 cos 3 ( x 4 ) sin ( x 4 ) h 2 ( x ) s 2 𝒶 + 1 2 3 c 2 sin ( x 2 ) s 𝒶 + 8 3 c 2 cos 3 ( x 4 ) sin ( x 4 ) s 𝒶 + 32 3 c 2 cos ( x 4 ) sin ( x 4 ) ( 1 8 sin 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) s 𝒶 c sin ( x 4 ) cos ( x 4 ) h 2 ( x ) s 𝒶 + 8 c ( 1 8 sin 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) h 2 ( x ) s 𝒶 4 c sin ( x 4 ) cos ( x 4 ) h 2 ( x ) s 𝒶 + 2 2 𝒶 + 5 c 2 sin ( x 4 ) cos ( x 4 ) sin 2 ( x 4 ) ( 1 8 sin 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) Γ ( 𝒶 + 1 2 ) 3 π Γ ( 𝒶 + 1 ) s 𝒶 4 c sin ( x 4 ) cos ( x 4 ) Γ ( 3 𝒶 + 1 ) h 2 ( x ) Γ ( 𝒶 + 1 ) Γ ( 2 𝒶 + 1 ) s 2 𝒶 + 8 c Γ ( 3 𝒶 + 1 ) ( 1 8 sin 2 ( x 4 ) 1 8 cos 2 ( x 4 ) ) h 2 ( x ) Γ ( 𝒶 + 1 ) Γ ( 2 𝒶 + 1 ) s 2 𝒶 3 4 2 𝒶 Γ ( 2 𝒶 + 1 2 ) h 2 ( x ) h 2 ( x ) π Γ ( 2 𝒶 + 1 ) s 2 𝒶 .
Then, by solving lim s s 2 𝒶 + 1 ( R e s 𝝕 2 ( x , s ) ) = 0 , this gives h 2 ( x ) = c 3 3 cos x 2 . So, the 2-nd Laplace series solution of (18) could be written as:
𝝕 2 ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 .
To construct the 3-rd Laplace series solution of (18) we should substitute 𝝕 3 ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 , into the 3-rd Laplace residual function of (20), such that:
( R e s 𝝕 3 ( x , s ) ) = 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 1 s 𝒶 { 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } D x 3 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } } 1 s 𝒶 { 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } D x 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } } 3 s 𝒶 { D x 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } D x 2 1 { 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + h 3 ( x ) s 3 𝒶 + 1 } } .
Following that, one can obtain h 3 ( x ) , via looking the solution of lim s s 3 𝒶 + 1 ( R e s 𝝕 3 ( x , s ) ) = 0 , to conclude that h 3 ( x ) = c 4 6 sin x 2 , and thus the 3-rd Laplace series solution of (18) could be expressed as:
𝝕 3 ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + c 4 6 s 3 𝒶 + 1 sin x 2 .
Continue in the similar manner for k = 4 , and by utilizing the result lim s s 4 𝒶 + 1 ( R e s 𝝕 4 ( x , s ) ) = 0 , then the fourth unknown function will be h 4 ( x ) = c 5 12 cos x 2 , and hence the 4-th Laplace series solution of (18) is given by the following expansion:
𝝕 4 ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + c 4 6 s 3 𝒶 + 1 sin x 2 c 5 12 s 4 𝒶 + 1 cos x 2 .
By repeating the previous algorithm for arbitrary k , and using MATHEMATICA Software Package 12, we can find out the unknown coefficient functions h k ( x ) , in the fractional expansion (19) has the following general forms:
h k ( x ) = {   1 3 c k + 1 2 k 2 sin x 2   :               k = 1 , 5 , 9 ,     1 3 c k + 1 2 k 2 sin x 2 :     k = 3 , 7 , 11 ,     1 3 c k + 1 2 k 2 cos x 2 :   k = 2 , 6 , 10 ,   1 3 c k + 1 2 k 2 cos x 2 :   k = 4 , 8 , 12 ,  
Therefore, the k-th Laplace series solution of (18) is given by the following expansion:
𝝕 k ( x , s ) = 8 c 3 s cos 2 x 4 2 c 2 3 s 𝒶 + 1 sin x 2 + c 3 3 s 2 𝒶 + 1 cos x 2 + c 4 6 s 3 𝒶 + 1 sin x 2 c 5 12 s 4 𝒶 + 1 cos x 2 + .
If we replace the term ( 8 c 3 s cos 2 x 4 ) by the term ( 4 c 3 s 4 c 3 s cos x 2 ) in (27) and perform some algebra iterations, then we get the following fractional expansion:
𝝕 k ( x , s ) = 4 c 3 ( 1 + sin x 2 ( c 2 1 s 𝒶 + 1 ( c 2 ) 3 1 s 3 𝒶 + 1 + ) ) 4 c 3 cos x 2 ( 1 ( c 2 ) 2 1 s 2 𝒶 + 1 + ( c 2 ) 4 1 s 4 𝒶 + 1 + ) .
Correspondingly, the k-th Laplace series solution of (28) could be expressed in the following finite series terms:
  𝝕 k ( x , s ) = 4 c 3 ( 1 + sin x 2 ( n = 0 k ( 1 ) n ( c 2 ) 2 n + 1 1 s ( 2 n + 1 ) 𝒶 + 1 ) ) 4 c 3 cos x 2 ( n = 0 k 1 ( 1 ) n ( c 2 ) 2 n 1 s 2 n 𝒶 + 1 ) ,     for   k = 1 , 3 , 5 , , 𝝕 k ( x , s ) = 4 c 3 ( 1 + sin x 2 ( n = 0 k 1 ( 1 ) n ( c 2 ) 2 n + 1 1 s ( 2 n + 1 ) 𝒶 + 1 ) ) 4 c 3 cos x 2 ( n = 0 k ( 1 ) n ( c 2 ) 2 n 1 s 2 n 𝒶 + 1 ) ,     for       k = 2 , 4 , 6 ,
Lastly, we operate the inverse Laplace transform to both sides of (29) to get the k-th approximate solution of the nonlinear time-fractional Rosenau–Hyman IVP (17) as:
W k ( x , t ) = 1 { 𝝕 k ( x , s ) } = 4 c 3 ( 1 + sin x 2 ( n = 0 k ( 1 ) n ( c t 𝒶 2 ) 2 n + 1 1 Γ ( ( 2 n + 1 ) 𝒶 + 1 ) ) ) 4 c 3 cos x 2 ( n = 0 k 1 ( 1 ) n ( c t 𝒶 2 ) 2 n 1 Γ ( 2 n 𝒶 + 1 ) )   for   k = 1 , 3 , 5 , , W k ( x , t ) = 1 { 𝝕 k ( x , s ) } = 4 c 3 ( 1 + sin x 2 ( n = 0 k 1 ( 1 ) n ( c t 𝒶 2 ) 2 n + 1 1 Γ ( ( 2 n + 1 ) 𝒶 + 1 ) ) ) 4 c 3 cos x 2 ( n = 0 k ( 1 ) n ( c t 𝒶 2 ) 2 n 1 Γ ( 2 n 𝒶 + 1 ) )   for   k = 2 , 4 , 6 ,
Consequently, the approximate solution of nonlinear time-fractional Rosenau–Hyman IVP (17) is given by:
W ( x , t ) = lim k W k ( x , t ) = 4 c 3 ( 1 + sin x 2 ( n = 0 ( 1 ) n ( c t 𝒶 2 ) 2 n + 1 1 Γ ( ( 2 n + 1 ) 𝒶 + 1 ) ) ) 4 c 3 cos x 2 ( n = 0 ( 1 ) n ( c t 𝒶 2 ) 2 n 1 Γ ( 2 n 𝒶 + 1 ) ) .
Particularly, if 𝒶 = 1 , in (31), then the general form of the approximate solution of (17) can be written as:
W ( x , t ) = 4 c 3 ( 1 + sin x 2 ( n = 0 ( 1 ) n ( c t 2 ) 2 n + 1 1 ( 2 n + 1 ) ! ) ) 4 c 3 cos x 2 ( n = 0 ( 1 ) n ( c t 2 ) 2 n 1 ( 2 n ) ! ) .
Which is fully in agreement with the Maclaurin series expansion of the exact solution W ( x , t ) = 8 c 3 cos 2 ( x c t 4 ) .
The efficiency and accuracy of the Laplace RPSM are demonstrated by calculating the absolute errors | W ( x , t ) W 7 ( x , t ) | , for standard case 𝒶 = 1 , at fixed values of the spatial coordinate variable x , and some chosen grid points of t , in [ 0 , 1 ] , as shown in Table 4. As we can see from the table, the obtained approximate solution coincides with the exact solution, by using only the seventh terms of the approximate solution.
For the purpose of numerical comparisons, Table 5 shows the absolute errors of the 5-th approximate solution to the nonlinear time-fractional Rosenau–Hyman IVP (17) by Laplace RPSM, Variational Iteration Method (VIM), and Homotopy Perturbation Method (HPM) [34] at standard case 𝒶 = 1 , for fixed value of x and some chosen mesh points of t . The numerical simulation given in Table 5 reveals that the absolute errors obtained by Laplace RPSM are smaller than other errors, and this emphasizes that the Laplace RPSM is more accurate in finding the exact solution of the nonlinear time-fractional Rosenau–Hyman IVP (17).
On the other hand, to show the effect of the fractional derivative to nonlinear time-fractional Rosenau–Hyman IVP (17), the graphs of the exact and 7-th approximate solutions for different values of 𝒶 is established in Figure 3. Further, the geometric behavior of the exact and 7-th approximate solutions are plotted in 3D surface plots for t [ 0 , 1 ] , and x [ 0 , 2 π ] , at various values of values of 𝒶 , as shown in Figure 4. This shows that from these figures the obtained approximate solution converges continuously to the standard-case 𝒶 = 1 as 𝒶 moves over ( 0 , 1 ) , and that the graphs of the behavior of the obtained 7-th approximate solution are consistent with each other, especially when considering the standard derivative.

5. Conclusions

In this article, the approximate analytical solution is constructed and analyzed for nonlinear time-fractional Kolmogorov, and Rosenau–Hyman equations with suitable initial conditions utilizing the Laplace RPSM under time-Caputo differentiability. The present approach is a modification of the fractional RPSM via coupling it to the Laplace transform operator. The benefit of utilizing the Laplace RPSM is that it gives more accurate convergence McLaurin series and needs a small size of computation without involving the discretization, perturbation, or any other physical restrictive conditions. Two well-known physical applications are tested to validate the applicability and superiority of the proposed method. The obtained approximate solutions are discussed via graphics and numeric simulation. The obtained results are compared with other well-known existing methods in the literature. Therefore, the results confirm that the Laplace RPSM is a straightforward and convenient tool to deal with the range of various non-linear time fractional-PDEs that arise in engineering and science problems. In future studies, the Laplace RPSM can be extended to find exact and approximate solutions for systems of FPDEs. Consequently, the application of the Laplace RPSM can be extended to handle physical models and dynamical models.

Author Contributions

Conceptualization, H.A. and M.A.; methodology, H.A.; software, H.A.; validation, M.A., A.I. and M.D.; writing—original draft preparation, H.A.; writing—review and editing, A.I. and M.D.; supervision, A.I. and M.D.; funding acquisition, A.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universiti Kebangsaan Malaysia (Project Code: DIP-2020-001).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mainardi, F.; Raberto, M.; Gorenflo, R.; Scalas, E. Fractional calculus and continuous-time finance II: The waiting-time distribution. Phys. A Stat. Mech. Its Appl. 2000, 287, 468–481. [Google Scholar] [CrossRef] [Green Version]
  2. He, J.H. Some applications of nonlinear fractional differential equations and their approximations. Bull. Sci. Technol. 1999, 15, 86–90. [Google Scholar]
  3. Caputo, M. Linear models of dissipation whose Q is almost frequency independent—II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  4. Oldham, K.B.; Spanier, J. The Fractional Calculus. In Integrations and Differentiations of Arbitrary Order; Descartes Press: Cambridge, MA, USA, 1974. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  6. Alabedalhadi, M.; Al-Smadi, M.; Abu Arqub, O.; Baleanu, D.; Momani, S. Structure of optical soliton solution for nonlinear resonant space-time Schrödinger equation in conformable sense with full nonlinearity term. Physica Scripta. 2020, 95, 105215. [Google Scholar] [CrossRef]
  7. Al-Smadi, M.; Abu Arqub, O. Computational algorithm for solving fredholm time-fractional partial integrodifferential equations of Dirichlet functions type with error estimates. Appl. Math. Comput. 2019, 342, 280–294. [Google Scholar] [CrossRef]
  8. Al-Smadi, M.; Abu Arqub, O.; Hadid, S. An attractive analytical technique for coupled system of fractional partial differential equations in shallow water waves with conformable derivative. Commun. Theor. Phys. 2020, 72, 085001. [Google Scholar] [CrossRef]
  9. Jleli, M.; Kumar, S.; Kumar, R.; Samet, B. Analytical approach for time fractional wave equations in the sense of Yang-Abdel-Aty-Cattani via the homotopy perturbation transform method. Alex. Eng. J. 2020, 59, 2859–2863. [Google Scholar] [CrossRef]
  10. Hasan, S.; Al-Smadi, M.; Dutta, H.; Momani, S.; Hadid, S. Multi-step reproducing kernel algorithm for solving Caputo–Fabrizio fractional stiff models arising in electric circuits. Soft Comput. 2022, 26, 3713–3727. [Google Scholar] [CrossRef]
  11. Hasan, S.; El-Ajou, A.; Hadid, S.; Al-Smadi, M.; Momani, S. Atangana-Baleanu fractional framework of reproducing kernel technique in solving fractional population dynamics system. Chaos Solitons Fractals 2020, 133, 109624. [Google Scholar] [CrossRef]
  12. Al-Smadi, M.; Abu Arqub, O.; Momani, S. A computational method for two-point boundary value problems of fourth-order mixed integrodifferential equations. Math. Probl. Eng. 2013, 2013, 832074. [Google Scholar] [CrossRef]
  13. Al-Smadi, M.; Dutta, H.; Hasan, S.; Momani, S. On numerical approximation of Atangana-Baleanu-Caputo fractional integro-differential equations under uncertainty in Hilbert Space. Math. Model. Nat. Phenom. 2021, 16, 41. [Google Scholar] [CrossRef]
  14. Al-Smadi, M.; Abu Arqub, O.; Zeidan, D. Fuzzy fractional differential equations under the Mittag-Leffler kernel differential operator of the ABC approach: Theorems and applications. Chaos Solitons Fractals 2021, 146, 110891. [Google Scholar] [CrossRef]
  15. Al-Smadi, M.; Abu Arqub, O.; Gaith, M. Numerical simulation of telegraph and Cattaneo fractional-type models using adaptive reproducing kernel framework. Math. Methods Appl. Sci. 2021, 44, 8472–8489. [Google Scholar] [CrossRef]
  16. Wu, L.; Xie, L.d.; Zhang, J.f. Adomian decomposition method for nonlinear differential-difference equations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 12–18. [Google Scholar] [CrossRef]
  17. Odibat, Z.; Momani, S. Application of variational iteration method to nonlinear differential equations of fractional order. Int. J. Nonlinear Sci. Numer. Simul. 2006, 7, 27–34. [Google Scholar] [CrossRef]
  18. Xu, H.; Liao, S.-J.; You, X.-C. Analysis of nonlinear fractional partial differential equations with the homotopy analysis method. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 1152–1156. [Google Scholar] [CrossRef]
  19. Hasan, S.; Al-Smadi, M.; El-Ajou, A.; Momani, S.; Hadid, S.; Al-Zhour, Z. Numerical approach in the Hilbert space to solve a fuzzy Atangana-Baleanu fractional hybrid system. Chaos Solitons Fractals 2021, 143, 110506. [Google Scholar] [CrossRef]
  20. Al-Smadi, M.; Freihat, A.; Khalil, H.; Momani, S.; Ali Khan, R. Numerical multistep approach for solving fractional partial differential equations. Int. J. Comput. Methods 2017, 14, 1750029. [Google Scholar] [CrossRef]
  21. Kumar, S. A new fractional modeling arising in engineering sciences and its analytical approximate solution. Alex. Eng. J. 2013, 52, 813–819. [Google Scholar] [CrossRef] [Green Version]
  22. Aljarrah, H.; Alaroud, M.; Ishak, A.; Darus, M. Adaptation of Residual-Error Series Algorithm to Handle Fractional System of Partial Differential Equations. Mathematics 2021, 9, 2868. [Google Scholar] [CrossRef]
  23. Saadeh, R.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, U.K.S. Application of fractional residual power series algorithm to solve Newell–Whitehead–Segel equation of fractional order. Symmetry 2019, 11, 1431. [Google Scholar] [CrossRef] [Green Version]
  24. Freihet, A.; Hasan, S.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, U.K.S. Toward computational algorithm for time-fractional Fokker–Planck models. Adv. Mech. Eng. 2019, 11, 1–11. [Google Scholar] [CrossRef]
  25. Al-Smadi, M.; Abu Arqub, O.; Hadid, S. Approximate solutions of nonlinear fractional Kundu-Eckhaus and coupled fractional massive Thirring equations emerging in quantum field theory using conformable residual power series method. Phys. Scr. 2020, 95, 105205. [Google Scholar] [CrossRef]
  26. Al-Smadi, M.; Abu Arqub, O.; 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]
  27. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, U.K.S. An analytical numerical method for solving fuzzy fractional Volterra integro-differential equations. Symmetry 2019, 11, 205. [Google Scholar] [CrossRef] [Green Version]
  28. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, U.K.S. Computational optimization of residual power series algorithm for certain classes of fuzzy fractional differential equations. Int. J. Differ. Equ. 2018, 2018, 8686502. [Google Scholar] [CrossRef]
  29. Alaroud, M.; Ahmad, R.R.; Din, U.K.S. An efficient analytical-numerical technique for handling model of fuzzy differential equations of fractional-order. Filomat 2019, 33, 617–632. [Google Scholar] [CrossRef] [Green Version]
  30. Bataineh, M.; Alaroud, M.; Al-Omari, S.; Agarwal, P. Series Representations for Uncertain Fractional IVPs in the Fuzzy Conformable Fractional Sense. Entropy 2019, 23, 1646. [Google Scholar] [CrossRef]
  31. El-Ajou, A. Adapting the Laplace transform to create solitary solutions for the nonlinear time-fractional dispersive PDEs via a new approach. Eur. Phys. J. Plu. 2021, 136, 229. [Google Scholar] [CrossRef]
  32. Alaroud, M. Application of Laplace residual power series method for approximate solutions of fractional IVP’s. Alex. Eng. J. 2022, 61, 1585–1595. [Google Scholar] [CrossRef]
  33. Alaroud, M.; Tahat, N.; Al-Omari, S.; Suthar, D.L.; Gulyaz-Ozyurt, S. An Attractive Approach Associated with Transform Functions for Solving Certain Fractional Swift-Hohenberg Equation. J. Funct. Spaces 2021, 2021, 3230272. [Google Scholar] [CrossRef]
  34. Yulita Molliq, R.; Noorani, M.S.M. Solving the fractional Rosenau-Hyman equation via variational iteration method and homotopy perturbation method. Int. J. Differ. Equ. 2012, 2012, 472030. [Google Scholar] [CrossRef] [Green Version]
  35. Korpinar, Z.; Inc, M.; Hınçal, E.; Baleanu, D. Residual power series algorithm for fractional cancer tumor models. Alex. Eng. J. 2020, 59, 1405–1412. [Google Scholar] [CrossRef]
  36. Nikan, O.; Avazzadeh, Z. Numerical simulation of fractional evolution model arising in viscoelastic mechanics. Appl. Num. Math. 2021, 169, 303–320. [Google Scholar] [CrossRef]
  37. Ali, A.; Abbas, M.; Akram, T. New group iterative schemes for solving the two-dimensional anomalous fractional sub-diffusion equation. J. Math. Comp. Sci. 2021, 22, 119–127. [Google Scholar] [CrossRef]
  38. Almatroud, A.; Ababneh, O.; Al-sawalha, M.M. Modify adaptive combined synchronization of fractional order chaotic systems with fully unknown parameters. J. Math. Comp. Sci. 2020, 21, 99–112. [Google Scholar] [CrossRef]
  39. Akram, T.; Abbas, M.; Ali, A. A numerical study on time fractional Fisher equation using an extended cubic B-spline approximation. J. Math. Comp. Sci. 2021, 22, 85–96. [Google Scholar] [CrossRef]
  40. Qu, H.; Rahman, M.; Ahmad, S.; Riaz, M.B.; Ibrahim, M.; Saeed, T. Investigation of fractional order bacteria dependent disease with the effects of different contact rates. Chaos Sol. Fract. 2022, 159, 112169. [Google Scholar] [CrossRef]
  41. El-Deeb, A.A.; Baleanu, D. Some new dynamic Gronwall–Bellman–Pachpatte type inequalities with delay on time scales and certain applications. J. Inequal. Appl. 2022, 2022, 45. [Google Scholar] [CrossRef]
  42. Kumar, P.; Baleanu, D.; Erturk, V.S.; Inc, M.; Govindaraj, V. A delayed plant disease model with Caputo fractional derivatives. Adv. Cont. Discr. Mod. 2022, 2022, 11. [Google Scholar] [CrossRef]
  43. Wang, G.; Wazwaz, A. On the modified Gardner type equation and its time fractional form. Chaos Solit. Fract. 2022, 155, 111694. [Google Scholar] [CrossRef]
  44. Wang, G. A new (3 + 1)-dimensional Schrödinger equation: Derivation, soliton solutions and conservation laws. Nonlin. Dyn. 2021, 104, 1595–1602. [Google Scholar] [CrossRef]
  45. Wang, G. Symmetry Analysis, Analytical Solutions And Conservation Laws Of A Generalized Kdv–Burgers–Kuramoto Equation And Its Fractional Version. Fractals 2021, 29, 2150101. [Google Scholar] [CrossRef]
  46. Wang, G. A novel (3 + 1)-dimensional sine-Gorden and a sinh-Gorden equation: Derivation, symmetries and conservation laws. Appl. Math. Lett. 2021, 113, 106768. [Google Scholar] [CrossRef]
  47. Wang, G. A (2 + 1)-dimensional sine-Gordon and sinh-Gordon equations with symmetries and kink wave solutions. Nucl. Phys. B 2020, 953, 114956. [Google Scholar] [CrossRef]
Figure 1. (a) Plots of Exact W ( x , t ) and W 10 ( x , t ) at x = 0 and with various 𝒶 values of IVP (7). (b) Plots of Exact W ( x , t ) and W 10 ( x , t ) at t = 0 and with various 𝒶 values of IVP (7).
Figure 1. (a) Plots of Exact W ( x , t ) and W 10 ( x , t ) at x = 0 and with various 𝒶 values of IVP (7). (b) Plots of Exact W ( x , t ) and W 10 ( x , t ) at t = 0 and with various 𝒶 values of IVP (7).
Mathematics 10 01980 g001
Figure 2. 3D Surface Plots of Exact solution W ( x , t ) , and the 10-th approximate solution W 10 ( x , t ) , for IVP (7), with t [ 0 , 1 ] , and x [ 0 , 1 ] , at various values of 𝒶 .
Figure 2. 3D Surface Plots of Exact solution W ( x , t ) , and the 10-th approximate solution W 10 ( x , t ) , for IVP (7), with t [ 0 , 1 ] , and x [ 0 , 1 ] , at various values of 𝒶 .
Mathematics 10 01980 g002aMathematics 10 01980 g002b
Figure 3. Profile the 7-th approximate solutions W 7 ( x , t ) , at various values of 𝒶 , for the nonlinear time-fractional Rosenau–Hyman IVP (17): (a) x = π ,   c = 4 ; (b) x = π 2 ,   c = 4 .
Figure 3. Profile the 7-th approximate solutions W 7 ( x , t ) , at various values of 𝒶 , for the nonlinear time-fractional Rosenau–Hyman IVP (17): (a) x = π ,   c = 4 ; (b) x = π 2 ,   c = 4 .
Mathematics 10 01980 g003
Figure 4. 3D Surface Plots of Exact solution W ( x , t ) , and the 7-th approximate solution W 7 ( x , t ) , for IVP (17), with t [ 0 , 1 ] , and x [ 0 , 2 π ] , at various values of 𝒶 .
Figure 4. 3D Surface Plots of Exact solution W ( x , t ) , and the 7-th approximate solution W 7 ( x , t ) , for IVP (17), with t [ 0 , 1 ] , and x [ 0 , 2 π ] , at various values of 𝒶 .
Mathematics 10 01980 g004
Table 1. Results of the 10-th approximate solution at different values of 𝒶 for Example 1.
Table 1. Results of the 10-th approximate solution at different values of 𝒶 for Example 1.
x t i 𝒶 = 1 𝒶 = 0.95   𝒶 = 0.85   𝒶 = 0.75   𝒶 = 0.65  
0 0.0 1.00 1.00 1.00 1.00 1.00
0.25 1.2840254167 1.3168989383 1.3960863710 1.4990355063 1.6365268759
0.50 1.6487212707 1.7072557012 1.8456231911 2.0217199431 2.2527959051
0.75 2.1170000155 2.2042026866 2.4091674423 2.6683879567 3.0068075604
1.0 2.7182818011 2.8399806687 3.1254929139 3.4858483992 3.9554385524
0.5 0.0 1.00 1.00 1.00 1.00 1.00
0.25 1.9260381250 1.9753484075 2.0941295565 2.2485532594 2.4547903138
0.50 2.4730819060 2.5608835519 2.7684347867 3.0325799147 3.3791938577
0.75 3.1755000232 3.3063040300 3.6137511635 4.0025819351 4.5102113406
1.0 4.0774227017 4.2599710030 4.6882393708 5.2287725988 5.9331578286
Table 2. The recurrence errors | W 8 ( x , t ) W 7 ( x , t ) | of the tenth approximate solution for Example 1.
Table 2. The recurrence errors | W 8 ( x , t ) W 7 ( x , t ) | of the tenth approximate solution for Example 1.
t i 𝒶 = 0.75   𝒶 = 0.85   𝒶 = 0.95   𝒶 = 1.00  
0.16 2.3301689 × 10 8 1.1468189 × 10 9 5.1702642 × 10 11 1.0651924 × 10 11
0.32 1.4913081 × 10 6 1.2779060 × 10 7 1.0030869 × 10 8 2.7269635 × 10 9
0.48 1.6986931 × 10 5 2.0133509 × 10 6 2.1859095 × 10 7 6.9889089 × 10 8
0.64 9.5443717 × 10 5 1.4239767 × 10 5 1.9461063 × 10 6 6.9810262 × 10 7
0.80 3.6408889 × 10 4 6.4936813 × 10 5 1.0609201 × 10 5 4.1610159 × 10 6
0.80 1.0871636 × 10 3 2.2434864 × 10 4 4.2409205 × 10 4 1.7891607 × 10 5
Table 3. Total computational cost for the obtained figures in Example 1.
Table 3. Total computational cost for the obtained figures in Example 1.
IDImage Size
(KB)
Maximum Memory
(MB)
Tracing Time
(s)
Total Cost
( MB × S )
Figure 1a 11.5 48 4.4 211.2
b 11.2 48 4.2 201.6
Figure 2a 30.6 49 4.5 220.5
b 34.7 50 4.4 220.0
c 37.4 54 5.0 270.0
d 36.2 51 4.5 229.5
Table 4. Numerical results for Example 2 at 𝒶 = 1 , and c = 2 5 , with different values of t .
Table 4. Numerical results for Example 2 at 𝒶 = 1 , and c = 2 5 , with different values of t .
x i   t i   W ( x , t ) W 7 ( x , t ) | W ( x , t ) W 7 ( x , t ) |  
0 0.15 1.0664266846661268 1.0664266846661268 0.0
0.30 1.0657069546321090 1.0657069546321067 2.220446049250313 × 10 15
0.45 1.0645081242730636 1.0645081242730066 5.706546346573305 × 10 14
0.60 1.0628312724553954 1.0628312724548266 5.688782778179302 × 10 13
0.75 1.0606779082325560 1.0606779082291666 3.389288849575678 × 10 12
0.90 1.0580499694869983 1.0580499694724268 1.457145515360025 × 10 11
π 6 0.15 0.1387518027223361 0.1387518027223363 2.220446049250313 × 10 16
0.30 0.1496922557134731 0.1496922557134746 1.498801083243961 × 10 15
0.45 0.1609779597794714 0.1609779597795144 4.293787547737793 × 10 14
0.60 0.1725987585484345 0.1725987585488619 4.274081089050696 × 10 13
0.75 0.1845441940858514 0.1845441940883893 2.537886567566261 × 10 12
0.90 0.1968035163060297 0.1968035163169100 1.088035217478022 × 10 11
Table 5. Numerical comparison of 5-th approximate solution IVP (17), at 𝒶 = 1 , x = π , and c = 0.5 .
Table 5. Numerical comparison of 5-th approximate solution IVP (17), at 𝒶 = 1 , x = π , and c = 0.5 .
t i   Laplace RPSMVIMHPM
0.1 1.11022 × 10 15 5.0000 × 10 10 1.0000 × 10 10
0.2 1.03251 × 10 13 5.0000 × 10 10 1.7360 × 10 9
0.3 1.76559 × 10 12 5.0000 × 10 10 1.3182 × 10 8
0.4 1.32255 × 10 11 1.0000 × 10 10 5.5542 × 10 8
0.5 6.30599 × 10 11 4.0000 × 10 10 1.6948 × 10 7
0.6 2.25934 × 10 10 7.0000 × 10 10 4.2165 × 10 7
0.7 6.64599 × 10 10 1.2000 × 10 9 9.1117 × 10 7
0.8 1.69218 × 10 9 2.1000 × 10 9 1.7761 × 10 6
0.9 3.85878 × 10 9 4.0000 × 10 9 3.1998 × 10 6
1.0 8.06643 × 10 9 8.6000 × 10 9 5.4173 × 10 6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aljarrah, H.; Alaroud, M.; Ishak, A.; Darus, M. Approximate Solution of Nonlinear Time-Fractional PDEs by Laplace Residual Power Series Method. Mathematics 2022, 10, 1980. https://doi.org/10.3390/math10121980

AMA Style

Aljarrah H, Alaroud M, Ishak A, Darus M. Approximate Solution of Nonlinear Time-Fractional PDEs by Laplace Residual Power Series Method. Mathematics. 2022; 10(12):1980. https://doi.org/10.3390/math10121980

Chicago/Turabian Style

Aljarrah, Hussam, Mohammad Alaroud, Anuar Ishak, and Maslina Darus. 2022. "Approximate Solution of Nonlinear Time-Fractional PDEs by Laplace Residual Power Series Method" Mathematics 10, no. 12: 1980. https://doi.org/10.3390/math10121980

APA Style

Aljarrah, H., Alaroud, M., Ishak, A., & Darus, M. (2022). Approximate Solution of Nonlinear Time-Fractional PDEs by Laplace Residual Power Series Method. Mathematics, 10(12), 1980. https://doi.org/10.3390/math10121980

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