Next Article in Journal
Novel Fractional Models Compatible with Real World Problems
Previous Article in Journal
Intrinsic Metric Formulas on Some Self-Similar Sets via the Code Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation

1
Laboratoire Lmdan, Département de Mathématiques de la Décision, Faculté des Sciences Economiques et Gestion, Université Cheikh Anta Diop de Dakar, BP 5683 Dakar Fann, Senegal
2
Centre de Recherche Economique Appliquées, Laboratoire Ingénierie Financiére et Economique (LIFE), Faculté des Sciences Économiques et de Gestion, Université Cheikh Anta Diop de Dakar, BP 5683 Dakar Fann, Senegal
*
Author to whom correspondence should be addressed.
Fractal Fract. 2019, 3(2), 14; https://doi.org/10.3390/fractalfract3020014
Submission received: 5 March 2019 / Revised: 25 March 2019 / Accepted: 26 March 2019 / Published: 27 March 2019

Abstract

:
In this paper, the approximate solutions of the fractional diffusion equations described by the fractional derivative operator were investigated. The homotopy perturbation Laplace transform method of getting the approximate solution was proposed. The Caputo generalized fractional derivative was used. The effects of the orders α and ρ in the diffusion processes was addressed. The graphical representations of the approximate solutions of the fractional diffusion equation and the fractional diffusion-reaction equation both described by the Caputo generalized fractional derivative were provided.

1. Introduction

Fractional calculus has many applications in physics and has attracted several researchers. The popularity of fractional calculus is due to the variety of the fractional derivatives operators. We can cite: the Caputo fractional derivative and the Riemann–Liouville fractional derivative [1,2], the Atangana–Baleanu fractional derivative [3,4,5,6], the Caputo–Fabrizio fractional derivative [7], and many others. The Atangana–Baleanu fractional derivative and the Caputo–Fabrizio fractional derivative were recently used in modeling physical phenomena [5,7,8,9,10,11]. The Atangana–Baleanu fractional derivative and the Caputo–Fabrizio fractional derivative are great compromises in modeling real world phenomena. The most relevant application of fractional calculus in physics is the fractional diffusion equation. In Reference [1], Podlubny introduced the diffusion equation using the fractional derivative operators. In Reference [12], Santoz studied the Fokker–Planck equation using the Atangana–Baleanu fractional derivative and the Caputo–Fabrizio fractional derivative. In Reference [13], Santoz used the fractional Prabhakar Derivative for studying the diffusion equation. In Reference [14], Santoz proposed a new approach of the random walk with the fractional diffusion equation. In Reference [7], Hristov gave a complete review of the fractional models and suggests a novel method for solving fractional diffusion equations.
The problem of solving differential equations is one of the most important field in fractional calculus. Many methods are used for getting the analytical solutions [15], the numerical solutions [16], and the approximate solutions. We have the Fourier sine transform [15,17]. We have the classical Laplace transform. We have the standard power series. We have the homotopy method [18,19] and others. Liao introduced the homotopy method in Reference [20]. The technique is favorite a favorite today for getting the approximate solutions of the fractional differential equations described by the fractional derivatives. Several results and applications of the homotopy method exist. In Reference [21], Das et al. solved the fractional diffusion equation described by the Caputo fractional derivative. They discussed the numerical results. They analyzed the convergence of the homotopy method. In Reference [18], Das et al. investigated getting the approximate solution of the fractional diffusion-reaction equation described by the Caputo fractional derivative. The authors used the optimal homotopy analysis method. They gave the graphical representations of the diffusion processes to illustrate the results. In Reference [22], Delgado et al. used the Laplace homotopy analysis for getting the approximate solutions of the linear partial differential equations described by the Caputo–Fabrizio fractional derivative and the Atangana–Baleanu fractional derivative. The main advancement of Delagado et al.’s work is the use of the Laplace transform in the application of the homotopy method. In Reference [23], Guo et al. presented the fractional variational homotopy perturbation iteration method. They applied their method for getting the solution of the fractional diffusion equation. They used the Riemann–Liouville fractional derivative in their studies. In Reference [24], Kumar et al. used the optimal homotopy analysis method for solving the nonlinear fractional diffusion equation described by the Caputo fractional derivative. In Reference [25], Yan used a modified perturbation homotopy method for getting the approximate solution of the fractional nonlinear heat transfer equation involving with the Caputo fractional derivative. Many other investigations exist relative to the homotopy methods, as in Abuasad [26] and Rajeev [27].
Other methods exist for solving the fractional diffusion equations. In Reference [28], Jassim et al. used the local fractional Laplace variational iteration method for solving the fractional diffusion equation and the fractional wave equation described by the local fractional derivative operators. In Reference [29], Khader et al. used the Legendre approximation for getting the approximate solution of the fractional diffusion-reaction equation. They considered the Caputo fractional derivative in their studies. In Reference [30], Hristov proposed the integral balance methods as the heat balance integral method and the double integral method for solving the fractional sub-diffusion equations.
Recently, the generalized fractional derivative was introduced in the literature by Udita [31]. Udita’s definition is the generalization of the Caputo and the Riemann–Liouville fractional derivative operators (see more investigations in [32,33,34]). In 2018, Fahd et al. introduced the Laplace transforms of the generalized fractional derivatives [35]. In our paper, we investigated the homotopy perturbation ρ -Laplace transform method for solving the fractional diffusion equation and the fractional diffusion-reaction equation described by the Caputo generalized fractional derivative. We used the Riemann–Liouville fractional integral in the classical homotopy perturbation method for solving the fractional differential equations. We had many iterations in the homotopy perturbation method. In all iterations, the Riemann–Liouville fractional integral was used for solving the fractional differential equations. The main difficulty in fractional calculus is to calculate the Riemann–Liouville fractional integral of a given function explicitly. It is not trivial and impossible in several cases. The novelty of the homotopy perturbation ρ -Laplace transform method is that, in all the iterations in our works, we replaced the Riemann–Liouville integrator by the ρ -Laplace transform for solving the fractional differential equations. This novel method was applied to solve the fractional diffusion equation and the fractional diffusion-reaction equation. We solved a classical open problem: getting the analytical solutions or approximates solutions of the diffusion equations. Note, many methods exist for solving the diffusion equations, such as the Fourier sine transform [36], the integral balance methods [30], the powers series, and many others. The application of the techniques cited above depends on the used boundary conditions. The methods have some limitations. In many cases, they cannot be applied. We proved our method is useful of getting the approximate solutions of the fractional diffusion equation and the fractional diffusion-reaction equation. The method takes into account the given boundary condition. The homotopy perturbation method has a history. The homotopy method was introduced by Liao [20]. Later, Shader et al. [37] proposed a modified version of the homotopy perturbation method. The authors in Reference [37] introduced a new interpretation of the parameter expansion in the homotopy perturbation method. Odibat et al. [38] applied the modified homotopy method for solving the fractional quadratic Riccati differential equation. Kumar et al. [39] applied the homotopy perturbation method for solving the fractional Black–Scholes European option pricing by using the Riemann–Liouville integral operator. See also Yavuz et al. [19]. Darzi et al. [40] proposed a new version of the homotopy perturbation method. The authors in Reference [40] introduced the fractional power series in the homotopy perturbation method. In our method, we proposed the ρ -Laplace transform for solving the fractional differential equations.
In Section 2, we recall the preliminary definitions and the functions for future use in the paper. In Section 3, we describe the homotopy Laplace transform method for solving the fractional differential equation represented by the Caputo generalized fractional derivative. In Section 4, we use the homotopy perturbation Laplace transform method for solving the generalized fractional diffusion equation. We analyzed the impact of the orders α and ρ in the diffusion processes. In Section 5, we use the homotopy perturbation Laplace transform method for solving the generalized fractional diffusion-reaction equation. In Section 6, we give the concluding remarks.

2. Preliminary Defintions and Lemmas

In this section, the generalized fractional integral, the left generalized fractional derivative, the ρ -Laplace transform, and the Mittag–Leffler function have been recalled.
The generalized fractional integral is recalled in the following definition. Let’s the function g : [ 0 , + [ R . As in [35], the generalized integral of order α , ρ > 0 of the function g is defined in the form:
I α , ρ g ( t ) = 1 Γ ( α ) 0 t t ρ s ρ ρ α 1 g ( s ) d s s 1 ρ ,
where Γ ( ) is the Gamma function, for a.e t > 0 , and 0 < α < 1 .
The generalized fractional derivative is recalled in the following definition. Let the function g : [ 0 , + [ R . As in Reference [35], the left generalized fractional derivative of order α , ρ > 0 of the function g is defined in the form:
D α , ρ g ( t ) = I 1 α , ρ g ( t ) = 1 Γ ( 1 α ) d d t 0 t t ρ s ρ ρ α g ( s ) d s s 1 ρ ,
where Γ ( ) is the Gamma function, for a.e t > 0 , and 0 < α < 1 .
The Caputo generalized fractional derivative is recalled in the following definition. Let the function g : [ 0 , + [ R . As in Reference [35], the Caputo generalized fractional derivative of order α , ρ > 0 of the function g is defined in the form:
D c α , ρ g ( t ) = 1 Γ ( 1 α ) 0 t t ρ s ρ ρ α γ n g ( s ) d s s 1 ρ ,
where Γ ( ) is the Gamma function, for a.e t > 0 , γ = t 1 ρ d d t and 0 < α < 1 .
For the rest of this paper, the Laplace transform will be used to help us in solving the fractional differential equations. The ρ -Laplace transform was recently introduced in the literature [35]. The ρ -Laplace transform of the Caputo generalized fractional derivative is expressed in the following form:
L ρ D c α , ρ g ( t ) = s α L ρ g ( t ) k = 0 n 1 s α k 1 γ k g ( 0 ) .
Besides, the ρ -Laplace transform of a given function g is described in the form:
L ρ g ( t ) ( s ) = 0 e s t ρ ρ g ( t ) d t t 1 ρ .
As in Reference [41], the Mittag–Leffler function with two parameters is defined as follows:
E α , β z = k = 0 z k Γ ( α k + β ) ,
where α > 0 , β R and z C . The classical exponential function is obtained with α = β = 1 .

3. Homotopy Perturbation and Laplace Transform Method

In this section, we describe the homotopy perturbation method and ρ -Laplace transform. We will use this to get the solutions of the fractional diffusion and the fractional diffusion-reaction equation, both represented by the Caputo generalized fractional derivative. The homotopy perturbation method for solving the differential equation was first proposed by Liao in [20]. Here, we combined both the homotopy perturbation method and the ρ -Laplace transform. Let the fractional differential equation described by the Caputo generalized fractional derivative given by:
D c α , ρ u ( x , t ) + L u ( x , t ) + N u ( x , t ) = g ( x , t ) ,
with the initial boundary condition defined as u ( x , 0 ) = f ( x ) , where L represents the linear operator which contains the integer or non-integer derivative operators, N denotes the nonlinear operator, and the function g is considered as the input term.
The homotopy perturbation is defined as the following construction [20]:
D c α , ρ u ( x , t ) + p L u ( x , t ) + N u ( x , t ) g ( x , t ) = 0 ,
where p 0 , 1 represents the homotopy perturbation parameter. It is straightforward to notice, when the parameter p = 0 , we recover the generalized fractional differential equation defined by:
D c α , ρ u ( x , t ) = 0 ,
with the initial boundary condition defined as u ( x , 0 ) = f ( x ) . Note that when the parameter p = 1 , we obtained the fractional differential equation described by the Caputo generalized fractional derivative defined in Equation (7). We used the parameter p to expand the solution of the fractional differential Equation (7), in the following form:
u ( x , t ) = u 0 ( x , t ) + p u 1 ( x , t ) + p 2 u 2 ( x , t ) + p 3 u 3 ( x , t ) + .
The next step is to put Equation (10) into Equation (7). The functions u 0 , u 1 , u 2 , u 3 , … are the solutions of the fractional differential equations described by the Caputo generalized fractional derivative:
p 0 : D c α , ρ u 0 ( x , t ) = 0 ; p 1 : D c α , ρ u 1 ( x , t ) = L u 0 ( x , t ) m 1 ( u 0 ( x , t ) ) + g ( x , t ) ; p 2 : D c α , ρ u 2 ( x , t ) = L u 1 ( x , t ) m 2 ( u 0 ( x , t ) , u 1 ( x , t ) ) ; p 3 : D c α , ρ u 3 ( x , t ) = L u 2 ( x , t ) m 3 ( u 0 ( x , t ) , u 1 ( x , t ) , u 2 ( x , t ) ) : ,
where the functions m 1 , m 2 , m 2 , … satisfy the following condition:
m ( u 0 ( x , t ) + p u 1 ( x , t ) + p 2 u 2 ( x , t ) + p 3 u 3 ( x , t ) + ) = m 1 ( ( u 0 ( x , t ) ) + p m 2 ( u 0 ( x , t ) , u 1 ( x , t ) ) + p 2 m 3 ( u 0 ( x , t ) , u 1 ( x , t ) , u 2 ( x , t ) ) + .
At each step, we applied the ρ -Laplace transform of getting the functions u 0 , u 1 , u 2 , u 3 , …. The boundary conditions of the fractional differential equation (Equation (10)) are given respectively by:
u 1 ( x , 0 ) = u 2 ( x , 0 ) = u 3 ( x , 0 ) = = 0 .

4. Approximate Solution of the Fractional Diffusion Equation

This section addresses the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative. We used both the homotopy perturbation method and the ρ -Laplace transform of the Caputo generalized fractional derivative to get the approximate solution of the fractional diffusion equation. The fractional diffusion equation under consideration is expressed as the following form:
D t α , ρ u ( x , t ) = 2 u ( x , t ) x 2 ,
with the initial boundary condition defined as for all x > 0 ,
u ( x , 0 ) = e x 1 .
We notice the Taylor expansion of the boundary condition at order 1, and we recover the classical boundary condition. That is u ( x , 0 ) = x . In other words, Equation (14) is the general boundary condition of the fractional diffusion equation defined by Equation (13).
In the first step of the resolution, we solve the fractional differential equation described by the Caputo generalized fractional derivative, expressed as:
D α , ρ u 0 ( x , t ) = 0 ,
with the initial boundary condition defined by u 0 ( x , t ) = e x 1 . When we apply the ρ -Laplace transform to both sides of Equation (15), we obtain:
s α u 0 ( x , s ) s α 1 u 0 ( x , 0 ) = 0 s α u 0 ( x , s ) = s α 1 u 0 ( x , 0 ) s α u 0 ( x , s ) = s α 1 e x 1 u 0 ( x , s ) = e x 1 s .
When we apply the inverse of ρ -Laplace transform to both sides of Equation (16), we obtain the analytical solution of Equation (15),
u 0 ( x , t ) = e x 1 .
In the second step, we solve the fractional differential equation described by the Caputo fractional derivative expressed as:
D α , ρ u 1 ( x , t ) = 2 u 0 ( x , t ) x 2 = e x ,
with the initial boundary condition defined by u 1 ( x , t ) = 0 . We apply the ρ -Laplace transform to both sides of Equation (18), and we obtain:
s α u 1 ( x , s ) s α 1 u 1 ( x , 0 ) = e x s s α u 1 ( x , s ) = e x s u 1 ( x , s ) = e x s 1 + α .
When we apply the inverse of ρ -Laplace transform to both sides of Equation (19), we obtain the analytical solution of equation:
u 1 ( x , t ) = e x Γ 1 + α t ρ ρ α .
In the third step, we solve the fractional differential equation described by the Caputo fractional derivative expressed as:
D α , ρ u 2 ( x , t ) = 2 u 1 ( x , t ) x 2 = e x Γ 1 + α t ρ ρ α ,
with initial boundary condition defined by u 2 ( x , t ) = 0 . We apply the ρ -Laplace transform to both sides of Equation (19), and we obtain:
s α u 2 ( x , s ) s α 1 u 2 ( x , 0 ) = e x s 1 + α s α u 2 ( x , s ) = e x s 1 + α u 2 ( x , s ) = e x s 1 + 2 α ,
When we apply the inverse of ρ -Laplace transform to both sides of Equation (19), we obtain the analytical solution of Equation (18),
u 2 ( x , t ) = e x Γ 1 + 2 α t ρ ρ 2 α .
We adopt the same reasoning in the other steps. Finally, the approximate solution of the fractional diffusion in Equation described by the Caputo generalized fractional derivative Equation (13), under the condition Equation (14), is given by:
u ( x , t ) = lim p 1 i = 0 n p i u i ( x , t ) u ( x , t ) = u 0 ( x , t ) + u 1 ( x , t ) + u 2 ( x , t ) + u ( x , t ) = e x 1 + e x Γ 1 + α t ρ ρ α + e x Γ 1 + 2 α t ρ ρ 2 α + u ( x , t ) = e x 1 + 1 Γ 1 + α t ρ ρ α + 1 Γ 1 + 2 α t ρ ρ 2 α + . . . 1 u ( x , t ) = e x E α t ρ ρ α 1 .
The approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative Equation (13) is given by:
u ( x , t ) = e x E α t ρ ρ α 1 .
We can notice that when α = ρ = 1 , we recover the analytical solution of the classical diffusion equation defined by:
t u ( x , t ) = 2 u ( x , t ) x 2 ,
with initial boundary conditions defined as for all x > 0 :
u ( x , 0 ) = e x 1 ,
given by:
u ( x , t ) = e x + t 1 .
In Figure 1, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative in the space coordinates and the time, when α = 0.5 and ρ = 1 .
In Figure 2, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative in the space coordinates and the time, when α = 1 and ρ = 0.5 .
In Figure 3, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative, when α 1 , ρ = 1 and t = 0.5 .
We observed all the curves increase and the order of the curves follow the increase of the order α , (see the direction depicted by the arrow). In general, when we fixed the order ρ = 1 , we noticed an acceleration effect generated by the orders α 1 .
In Figure 4, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative, when α 1 , ρ = 1 and t = 0.5 . We observed all the curves increase and the order of the curves follow the increase of the order α , so we have the same behaviors as in Figure 3. But in general, when we fixed the order ρ = 1 , we noticed a retardation effect generated by the orders α 1 .
In Figure 5, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative, when ρ 1 , α = 1 and t = 0.5 .
We observed all the curves increase and the order of the curves follow the increase of the order α (see the direction depicted by the arrow). In general, regarding the curve obtained with α = 1 and ρ = 1 , we notice an acceleration effect generated by the orders ρ 1 .
In Figure 6, we depict the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative, when ρ 1 , α = 1 and t = 0.5 . We observed all the curves increase and the order of the curves follow the increase of the orders α . In general, when we fixed the order α = 1 , we noticed a retardation effect generated by the order ρ 1 .

5. Approximate Solutions of the Fractional Diffusion-Reaction Equation

This section addresses the approximate solution of the fractional diffusion-reaction equation described by the Caputo generalized fractional derivative. We used both the homotopy perturbation method and the Laplace transform of the Caputo generalized fractional derivative of getting the solution of the fractional diffusion-reaction equation. The fractional diffusion equation under consideration is expressed as:
D t α , ρ u ( x , t ) = 2 u ( x , t ) x 2 + u ( 1 u ) ,
with initial boundary condition defined as for all x > 0 :
u ( x , 0 ) = 1 / 2 .
In the first step of the resolution, we solve the fractional differential equation described by the Caputo generalized fractional derivative expressed as:
D α , ρ u 0 ( x , t ) = 0 ,
with the initial boundary condition defined by u 0 ( x , t ) = 1 / 2 . We apply the ρ -Laplace transform to both sides of Equation (31), and we obtain:
s α u 0 ( x , s ) s α 1 u 0 ( x , 0 ) = 0 s α u 0 ( x , s ) = s α 1 u 0 ( x , 0 ) s α u 0 ( x , s ) = 1 2 s α 1 u 0 ( x , s ) = 1 2 s
We apply the inverse of ρ -Laplace transform to both sides of Equation (32), and we obtain the analytical solution of Equation (31),
u 0 ( x , t ) = 1 2 .
In the second step, we solve the fractional differential equation described by the Caputo fractional derivative expressed as:
D α , ρ u 1 ( x , t ) = 2 u 0 ( x , t ) x 2 + u 0 ( 1 u 0 ) = 1 4 ,
with initial boundary condition defined by u 1 ( x , t ) = 0 . We apply the ρ -Laplace transform to both sides of Equation (34), and we obtain:
s α u 1 ( x , s ) s α 1 u 1 ( x , 0 ) = 1 4 s s α u 1 ( x , s ) = 1 4 s u 1 ( x , s ) = 1 4 s 1 + α .
We apply the inverse of ρ -Laplace transform to both sides of Equation (35), and we obtain the analytical solution of Equation (34),
u 1 ( x , t ) = 1 4 Γ 1 + α t ρ ρ α .
In the third step, we solve the fractional differential equation described by the Caputo generalized fractional derivative expressed as:
D α , ρ u 2 ( x , t ) = 2 u 1 ( x , t ) x 2 + u 1 ( 1 u 1 ) = 1 4 Γ 1 + α t ρ ρ α 1 4 Γ 1 + α 2 t ρ ρ 2 α ,
with the initial boundary condition defined by u 2 ( x , t ) = 0 . We apply the ρ -Laplace transform to both sides of Equation (37), we obtain
s α u 2 ( x , s ) s α 1 u 2 ( x , 0 ) = 1 4 s 1 + α Γ 1 + 2 α 4 Γ 1 + α 2 1 s 1 + 2 α s α u 2 ( x , s ) = 1 4 s 1 + α Γ 1 + 2 α 4 Γ 1 + α 2 1 s 1 + 2 α u 2 ( x , s ) = 1 4 s 1 + 2 α Γ 1 + 2 α 4 Γ 1 + α 2 1 s 1 + 3 α
We apply the inverse of ρ -Laplace transform to both sides of Equation (38), and we obtain the analytical solution of Equation (37),
u 2 ( x , t ) = 1 4 Γ 1 + 2 α t ρ ρ 2 α Γ 1 + 2 α 4 Γ 1 + α 2 1 Γ 1 + 3 α t ρ ρ 3 α .
We adopt the same reasoning in the other steps. Finally, the approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative Equation (29), is given by:
u ( x , t ) = lim p 1 i = 0 n p i u i ( x , t ) u ( x , t ) = u 0 ( x , t ) + u 1 ( x , t ) + u 2 ( x , t ) + u ( x , t ) = 1 4 + 1 4 E α t ρ ρ α Γ 1 + 2 α 4 Γ 1 + α 2 1 Γ 1 + 3 α t ρ ρ 3 α .
The approximate solution of the fractional diffusion equation described by the Caputo generalized fractional derivative Equation (29) is given by:
u ( x , t ) = 1 4 + 1 4 E α t ρ ρ α h ( t ) ,
where
h ( t ) = Γ 1 + 2 α 4 Γ 1 + α 2 1 Γ 1 + 3 α t ρ ρ 3 α .
We can notice when α = ρ = 1 , we recover the analytical solution of the considered diffusion equation given by
u ( x , t ) = 1 4 + 1 4 e t 1 48 t 3 .
In Figure 7, we depict the curve of the analytical solution of the diffusion-reaction equation obtained when α = ρ = 1 . We observed the curve increase in time.

6. Convergence Analysis of the Method

In this section, we analyze the stability and the convergence of the homotopy perturbation ρ -Laplace transform method. The solution of the fractional diffusion equation (Equation (13)) and the fractional diffusion-reaction equation (Equation (29)) converge when the series defined in Equation (24) and (40) converge. To study the convergence of the homotopy perturbation ρ -Laplace transform method, we defined the error function expressed as:
e ( x , t ) = u e x a c t u a p p r o ,
where the function u e x a c t represents the exact solution and u a p p r o represents the approximate solution of the fractional diffusion equation. For the illustration of the convergence, we consider the order α = ρ = 1 . The exact solution is given by u e x a c t = e x + t 1 . The approximate solution is obtained with Equation (24), after three iterations. The results are consigned in Table 1.
We noted the results were in good agreement. Thus, the convergence of the homotopy perturbation ρ -Laplace transform method converged.

7. Conclusion

In this paper, the approximate solutions of the fractional diffusion equation and the fractional diffusion-reaction equation, both described by the Caputo fractional derivative, were provided, using the homotopy perturbation method. The impact of the order α and ρ in the diffusion processes was addressed. In general, the orders α and ρ had retardation or an acceleration effect in the diffusion processes. The behavior of the solutions of the fractional diffusion equation and the fractional diffusion-reaction equations both described by the Caputo fractional derivative in coordinate space and time were provided graphically. The problem of getting the analytical solution or the approximate solution of the fractional nonlinear sub-diffusion equation has been an open problem in fractional calculus. The problem is not trivial. Applying our method of combining both the homotopy perturbation method and the ρ -Laplace transform to get the approximate solution seemed adequate. We investigated using this method for solving this class fractional diffusion equations. In forthcoming works, some investigations will be done in this direction.

Author Contributions

All authors contributed equally and significantly in writing this article. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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: New York, NY, USA, 1998. [Google Scholar]
  2. Sene, N. Mittag–Leffler input stability of fractional differential equations and its applications. Discret. Contin. Dyn. Syst. Ser. S 2019. [Google Scholar] [CrossRef]
  3. Atangana, A.; Gomez-Aguilar, J.F. A new derivative with normal distribution kernel: Theory, methods and applications. Phys. A Stat. Mech. Appl. 2017, 476, 1–14. [Google Scholar] [CrossRef]
  4. Atangana, A.; Gomez-Aguilar, J.F. Decolonisation of fractional calculus rules: Breaking commutativity and associativity to capture more natural phenomena. Eur. Phys. J. Plus 2019, 133, 1–22. [Google Scholar] [CrossRef]
  5. Escamilla, A.C.; Gómez-Aguilar, J.F.; Baleanu, D.; Córdova-Fraga, T.; Escobar-Jiménez, R.F.; Olivares-Peregrino, V.H.; Qurashi, M.M.A. Bateman–Feshbach Tikochinsky and Caldirola–Kanai Oscillators with New Fractional Differentiation. Entropy 2017, 19, 1–21. [Google Scholar]
  6. Gomez-Aguilar, J.F.; Atangana, A. New insight in fractional differentiation: Power, exponential decay and Mittag–Leffler laws and applications. Eur. Phys. J. Plus 2018, 132, 13. [Google Scholar] [CrossRef]
  7. Hristov, J. Derivatives with non-singular kernels from the Caputo–Fabrizio definition and beyond: Appraising analysis with emphasis on diffusion models. Front. Fract. Calc. 2017, 1, 270–342. [Google Scholar]
  8. Delgado, V.F.M.; Taneco-Hernández, M.A.; Gómez-Aguilar, J.F. On the solutions of fractional order of evolution equations. Eur. Phys. J. Plus 2017, 132, 47. [Google Scholar] [CrossRef]
  9. Hristov, J. Response functions in linear viscoelastic constitutive equations and related fractional operators. Math. Model. Nat. Phenom. 2019, 14, 305. [Google Scholar] [CrossRef]
  10. Gomez-Aguilar, J.F.; Torres, L.; Yépez-Martínez, H.; Baleanu, D.; Reyes, J.M.; Sosa, I.O. Fractional Liénard type model of a pipeline within the fractional derivative without singular kernel. Adv. Differ. Equ. 2016, 2016, 173. [Google Scholar] [CrossRef]
  11. Yepez-Martıneza, H.; Gomez-Aguilar, J.F.; Sosa, I.O.; Reyes, J.M.; Torres-Jimenez, J. The Feng’s first integral method applied to the nonlinear mKdV space-time fractional partial differential equation. Revista Mexicana de Fısica 2016, 62, 310–316. [Google Scholar]
  12. Santos, M.D.; Gomez, I.S. A fractional Fokker–Planck equation for non-singular kernel operators. J. Stat. Mech. Theory Exp. 2018, 2018, 123205. [Google Scholar] [CrossRef]
  13. Santos, M.D. Fractional Prabhakar Derivative in Diffusion Equation with Non-Static Stochastic Resetting. Physics 2019, 1, 5. [Google Scholar] [CrossRef]
  14. Santos, M.D. Non-Gaussian Distributions to Random Walk in the Context of Memory Kernels. Fractal Fract. 2018, 2, 20. [Google Scholar] [CrossRef]
  15. Sene, N. Analytical solutions of Hristov diffusion equations with non-singular fractional derivatives. Chaos 2019, 29, 023112. [Google Scholar] [CrossRef] [PubMed]
  16. Owolabi, K.M.; Atangana, A. Robustness of fractional difference schemes via the Caputo subdiffusion-reaction equations. Chaos Solitons Fractals 2018, 111, 119–127. [Google Scholar] [CrossRef]
  17. Sene, N. Stokes’ first problem for heated flat plate with Atangana–Baleanu fractional derivative. Chaos Solitons Fractals 2018, 117, 68–75. [Google Scholar] [CrossRef]
  18. Das, S.; Vishal, K.; Gupta, P.K. Approximate analytical solution of diffusion equation with fractional time derivative using optimal homotopy analysis method. Surv. Math. Appl. 2013, 8, 35–49. [Google Scholar]
  19. Yavuz, M.; Özdemir, N. A different approach to the European option pricing model with new fractional operator. Math. Model. Nat. Phenom. 2018, 13, 12. [Google Scholar] [CrossRef]
  20. Liao, S.J. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. 2004, 147, 499–513. [Google Scholar] [CrossRef]
  21. Das, S.; Vishal, K.; Gupta, P.K.; Saha Ray, S. Homotopy analysis method for solving fractional diffusion equation. Int. J. Appl. Math. Mech. 2011, 7, 28–37. [Google Scholar]
  22. Delgado, V.F.M.; Aguilar, J.F.G.; Martínez, H.Y.; Baleanu, D.; Jimenez, R.F.E.; Peregrino, V.H.O. Laplace homotopy analysis method for solving linear partial differential equations using a fractional derivative with and without kernel singular. Adv. Differ. Equ. 2016, 2016, 164. [Google Scholar] [CrossRef]
  23. Guo, S.; Mei, L.; Li, Y. Fractional variational homotopy perturbation iteration method and its application to a fractional diffusion equation. Appl. Math. Comput. 2013, 219, 5909–5917. [Google Scholar] [CrossRef]
  24. Vishala, K.; Das, S. Solution of the Nonlinear Fractional Diffusion Equation with Absorbent Term and External Force Using Optimal Homotopy-Analysis Method. Zeitschrift für Naturforschung A 2012, 67a, 203–209. [Google Scholar] [CrossRef]
  25. Yan, L.M. Modified Homotopy Perturbation Method Coupled with Laplace Transform for Fractional Heat Transfer and Porous Media Equations. Therm. Sci. 2013, 17, 1409–1414. [Google Scholar] [CrossRef]
  26. Abuasad, S.; Hashim, I. Homotopy Decomposition Method for Solving Higher-Order Time-Fractional Diffusion Equation via Modified Beta Derivative. Sains Malaysiana 2018, 47, 2899–2905. [Google Scholar] [CrossRef]
  27. Kushwaha, M.S. Homotopy perturbation method for a limit case Stefan problem governed by fractional diffusion equation. Appl. Math. Model. 2013, 37, 3589–3599. [Google Scholar]
  28. Jassim, H.K.; Ünlü, C.; Moshokoa, S.P.; Khalique, C.M. Local Fractional Laplace Variational Iteration Method for Solving Diffusion and Wave Equations on Cantor Sets within Local Fractional Operators. Math. Probl. Eng. 2015, 2015, 309870. [Google Scholar] [CrossRef]
  29. Khader, M.M.; Sweilam, N.H.; Mahdy, A.M.S. An Efficient Numerical Method for Solving the Fractional Diffusion Equation. J. Appl. Math. Bioinf. 2011, 1, 1–12. [Google Scholar]
  30. Hristov, J. Approximate solutions to fractional subdiffusion equations. Eur. Phys. J. Spec. Top. 2011, 193, 229–243. [Google Scholar] [CrossRef]
  31. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef]
  32. Gambo, Y.Y.; Ameen, R.; Jarad, F.; Abdeljawad, T. Existence and uniqueness of solutions to fractional differential equations in the frame of generalized Caputo fractional derivatives. Adv. Differ. Equ. 2018, 2018, 134. [Google Scholar] [CrossRef]
  33. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their Caputo modification. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. [Google Scholar] [CrossRef]
  34. Jarad, F.; Gurlu, E.U.; Abdeljawad, T.; Baleanu, D. On a new class of fractional operators. Adv. Differ. Equ. 2017, 2017, 247. [Google Scholar] [CrossRef]
  35. Fahd, J.; Abdeljawad, T. A modified Laplace transform for certain generalized fractional operators. Results Nonlinear Anal. 2018, 2, 88–98. [Google Scholar]
  36. Sene, N. Solution of fractional diffusion equations and Cattaneo–Hristov model. Int. J. Anal. Appl. 2019, 17, 191–207. [Google Scholar]
  37. Shaher, M.; Odibat, Z. Homotopy perturbation method for nonlinear partial differential equations of fractional order. Phys. Lett. A 2007, 365, 345–350. [Google Scholar]
  38. Odibat, Z.; Momani, S. Modified homotopy perturbation method: Application to quadratic Riccati differential equation of fractional order. Chaos Solitons Fractals 2008, 36, 167–174. [Google Scholar] [CrossRef]
  39. Kumar, S.; Kumar, D.; Singh, J. Numerical computation of fractional Black–Scholes equation arising in financial market. Egypt. J. Basic Appl. Sci. 2014, 1, 177–183. [Google Scholar] [CrossRef]
  40. Darzi, R.; Agheli, B. Homotopy analysis of transient magneto-bio-fluid dynamics of micropolar squeeze film in a porous medium: A model for magneto-bio-rheological lubrication. Math. Morav. 2018, 22, 93–105. [Google Scholar] [CrossRef]
  41. Sene, N. Exponential form for Lyapunov function and stability analysis of the fractional differential equations. J. Math. Comput. Sci. 2018, 18, 388–397. [Google Scholar] [CrossRef]
Figure 1. Approximate solutions of the fractional diffusion equation, α = 0.5 and ρ = 1 .
Figure 1. Approximate solutions of the fractional diffusion equation, α = 0.5 and ρ = 1 .
Fractalfract 03 00014 g001
Figure 2. Approximate solutions of the fractional diffusion equation, α = 1 , ρ = 0.5 and t = 0.5 .
Figure 2. Approximate solutions of the fractional diffusion equation, α = 1 , ρ = 0.5 and t = 0.5 .
Fractalfract 03 00014 g002
Figure 3. Approximate solutions of the fractional diffusion equation, α 1 , ρ = 1 and t = 0.5 .
Figure 3. Approximate solutions of the fractional diffusion equation, α 1 , ρ = 1 and t = 0.5 .
Fractalfract 03 00014 g003
Figure 4. Approximate solutions of the fractional diffusion equation, α 1 , ρ = 1 and t = 0.5 .
Figure 4. Approximate solutions of the fractional diffusion equation, α 1 , ρ = 1 and t = 0.5 .
Fractalfract 03 00014 g004
Figure 5. Approximate solutions of the fractional diffusion equation, α = 1 , ρ 1 and t = 0.5 .
Figure 5. Approximate solutions of the fractional diffusion equation, α = 1 , ρ 1 and t = 0.5 .
Fractalfract 03 00014 g005
Figure 6. Approximate solutions of the fractional diffusion equation, α = 1 , ρ 1 and t = 0.5 .
Figure 6. Approximate solutions of the fractional diffusion equation, α = 1 , ρ 1 and t = 0.5 .
Fractalfract 03 00014 g006
Figure 7. Approximate solutions of the fractional diffusion-reaction equation.
Figure 7. Approximate solutions of the fractional diffusion-reaction equation.
Fractalfract 03 00014 g007
Table 1. Numerical values of the error function.
Table 1. Numerical values of the error function.
x00.10.20.30.4
t0.10.10.10.10.1
u e x a c t 0.105170.221400.349850.491820.64872
u a p p r o 0.105000.221210.349650.491590.64846
e = u e x a c t u a p p r o 1.7 × 10 4 1.9 × 10 4 2 × 10 4 2.3 × 10 4 2.6 × 10 4

Share and Cite

MDPI and ACS Style

Sene, N.; Fall, A.N. Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation. Fractal Fract. 2019, 3, 14. https://doi.org/10.3390/fractalfract3020014

AMA Style

Sene N, Fall AN. Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation. Fractal and Fractional. 2019; 3(2):14. https://doi.org/10.3390/fractalfract3020014

Chicago/Turabian Style

Sene, Ndolane, and Aliou Niang Fall. 2019. "Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation" Fractal and Fractional 3, no. 2: 14. https://doi.org/10.3390/fractalfract3020014

APA Style

Sene, N., & Fall, A. N. (2019). Homotopy Perturbation ρ-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation. Fractal and Fractional, 3(2), 14. https://doi.org/10.3390/fractalfract3020014

Article Metrics

Back to TopTop