Next Article in Journal
Revisiting the Schrödinger–Dirac Equation
Previous Article in Journal
Computational Simulation and Parametric Analysis of the Effectiveness of Ternary Nano-composites in Improving Magneto-Micropolar Liquid Heat Transport Performance
Previous Article in Special Issue
Modified Exp-Function Method to Find Exact Solutions of Microtubules Nonlinear Dynamics Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Analysis of Fractional-Order Fokker–Planck Equation

1
Department of Mathematical Sciences College of Sciences, Princess Nourah Bint Abdulrahman University, Riyadh 11546, Saudi Arabia
2
Department of Mathematics, Abdul Wali khan University, Mardan 23200, Pakistan
3
Mathematics Department, Khalifa University, Abu Dhabi 127788, United Arab Emirates
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 430; https://doi.org/10.3390/sym15020430
Submission received: 2 January 2023 / Revised: 13 January 2023 / Accepted: 21 January 2023 / Published: 6 February 2023
(This article belongs to the Special Issue Differential/Difference Equations and Its Application)

Abstract

:
The importance of partial differential equations in physics, mathematics and engineering cannot be emphasized enough. Partial differential equations are used to represent physical processes, which are then solved analytically or numerically to examine the dynamical behaviour of the system. The new iterative approach and the Homotopy perturbation method are used in this article to solve the fractional order Fokker–Planck equation numerically. The Caputo sense is used to characterize the fractional derivatives. The suggested approach’s accuracy and applicability are demonstrated using illustrations. The proposed method’s accuracy is expressed in terms of absolute error. The proposed methods are found to be in good agreement with the exact solution of the problems using graphs and tables. The results acquired using the given approaches are also obtained at various fractional orders of the derivative. It is observed from the graphs and tables that fractional order solutions converge to an integer solution when the fractional orders approach the integer-order of the problems. The tabular and graphical view for the given problems is obtained through Maple. The presented approaches can be applied to existing non-linear fractional partial differential equations due to their accurate, simple and straightforward implementation.

1. Introduction

The theory of fractional calculus (FC) has received a lot of interest in recent years because of its applications to complex systems. The simulation of significant-world issues employing fractional order derivatives gives higher accuracy than modeling involving integer-order derivatives, according to fractional derivative principles. FC refers to the background and non-local dispersed effects of any physical system in phenomena such as wave motion analysis, solitary waves, phase turbulence in reaction–diffusion schemes [1,2,3,4], chaotic drifting waves induced by photon collision [5], wrinkled flame front propagation [6], time fractional-coupled mKdV equation [7,8,9], fractional order wave equations [10] and fractional space-time diffusion equation [11,12,13]. Fractional Differential Equations (FDEs) have received a lot of interest from mathematicians because they enable fractional modeling of various natural processes [14,15,16]. As a result, the use of FDEs to represent many physical systems and processes has increased, such as coloured noise [17], economics [18], earthquake oscillation [19] and bioengineering [20]. Rheology [21], control theory [22], visco-elastic materials [23], signal processing [24], polymers [25], damping method [26] and so on are some of the additional applications [27,28,29,30]. The fractional differential equation is an effective tool for modeling nonlinear phenomena in scientific and engineering models. In applied mathematics and engineering, partial differential equations, especially nonlinear ones, have been utilized to model a vast array of scientific phenomena [31,32,33]. Parallel to their work in the physical sciences, researchers were able to identify and model a vast array of relevant and real-world physical problems using fractional order partial differential equations (FPDEs). It has long been asserted how crucial it is to establish estimates for them using numerical or analytical techniques. Consequently, symmetry analysis is a useful method for comprehending partial differential equations, particularly when examining equations generated from accounting-related mathematical concepts [34,35,36]. Even though symmetry is the cornerstone of nature, the majority of natural observations lack symmetry. A sophisticated strategy for concealing symmetry is the appearance of unexpected symmetry-breaking events. Two types of symmetry exist: finite and infinitesimal. Both discrete and continuous finite symmetries exist. Space is a continuous transformation, but parity and temporal inversion are discrete natural symmetries. Patterns have always intrigued mathematicians. In the seventeenth century, classification of spatial and planar patterns gained significant traction. Unfortunately, exact solution of fractional nonlinear differential equations has proven to be quite difficult.
The Fokker–Planck equation (FPE) was developed by Adriaan Fokker and Max Planck to describe the time evolution of the probability density function of a particle’s position and velocity and it is one of the most extensively used statistical physics equations [37]. FPE appears in a variety of natural science domains; Brownian motion [38] and the diffusion model of chemical reactions [39] are now widely used in physics, chemistry, engineering and biology in various modified forms. The FPE first appears in kinetic theory [40], where it represents the behavior of one-particle distribution function of a dilute gas with long-range collisions, such as a Coulomb gas. Some applications of this type of equation can be found in research by He and Wu [41], Jumarie [42], Kamitani and Matsuba [43], Xu et al. [44] and Zak [45].
Among these applications, we considered Fokker–Planck equations of fractional order with the general form
ξ ϑ β ( ϕ , ϑ ) = L ξ ϕ ( ϕ , ϑ ) + ξ ϕ ϕ ( ϕ , ϑ ) + N ξ ϕ ϕ ( ϕ , ϑ ) , ϕ , ϑ > 0 , β ( 0 , 1 ] ,
with initial source
ξ ( ϕ , 0 ) = ζ ( ϕ ) .
The function ( ϕ , ϑ ) is assumed to be a causal function of time and space, i.e., vanishing for ϕ < 0 and ϑ < 0 , β is the parameter describing the order of the fractional time and space derivative. Fokker and Planck proposed the Fokker–Planck equation (Equation (1)) to describe brownian motion of particles [46]. The Fokker–Plank equation, which explains solute transport, depicts the change in probability of a random function in time and space. PDEs of both time and space fractional order are used to describe a variety of phenomena, including wave propagation, continuous random walks, charge carrier transport in amorphous semiconductors, anomalous diffusion, ribosome mobility along mRNA and pattern generation in polymeric networks [47]. In the current study [48], we use both the innovative iterative approach offered by Gejji and Jafari [49] and the homotopy perturbation transform method proposed by Madani et al. [50], Khan and Wu. Daftardar-Gejji and Jafari introduced a new iterative method for finding numerical solutions to nonlinear functional equations in 2006 [51]. Many nonlinear differential equations of integer and fractional order [52] and fractional boundary value problems have been solved using the iterative method. In a simple manner, the second strategy combines the Elzaki transformation, the homotopy perturbation method and He’s polynomials. He [53,54] invented the homotopy perturbation technique (HPM), which is a series expansion approach for solving nonlinear partial differential equations. To ensure convergence of approximation series over a certain interval of physical parameters, the HPM employs a so-called convergence-control parameter.
The rest of this work is arranged in the following manner. The Abel–Riemann fractional derivative, Caputo fractional derivative and Elzaki transform are all defined in Section 1. The new iterative transform method for solving fractional partial differential equations is described in Section 2. The Homotopy perturbation transform is described in Section 3 for solving fractional partial differential equations. In Section 4, we show five examples of how the approaches can be used to solve Fokker–Planck equations. Section 5 presents the conclusion.

2. Basic Definitions

2.1. Definition

The fractional derivative D β in the Abel–Riemann sense having order β is given as [55]
D β μ ( φ ) = d κ d φ κ μ ( φ ) , β = κ , 1 Γ ( κ β ) d d φ κ 0 φ μ ( φ ) ( φ ϕ ) β κ + 1 d ϕ , κ 1 < β < κ ,
where κ Z + , β R + and
D β μ ( φ ) = 1 Γ ( β ) 0 φ ( φ ϕ ) β 1 μ ( ϕ ) d ϕ , 0 < β 1 .

2.2. Definition

The fractional integration operator κ ϕ in Abel–Riemann sense is defined as [55]
κ β μ ( φ ) = 1 Γ ( β ) 0 φ ( φ ϕ ) β 1 μ ( φ ) d φ , φ > 0 , β > 0 .
having properties:
κ β φ κ = Γ ( κ + 1 ) Γ ( κ + β + 1 ) φ κ + ϕ , D β φ κ = Γ ( κ + 1 ) Γ ( κ β + 1 ) φ κ ϕ .

2.3. Definition

The Caputo derivative D β of fractional order β is given as [55]
C D β μ ( φ ) = 1 Γ ( κ β ) 0 φ μ κ ( ϕ ) ( φ ϕ ) β κ + 1 d ϕ , κ 1 < β < κ , d κ d φ κ μ ( φ ) , κ = β .
with the properties
κ φ β D φ β g ( φ ) = g ( φ ) k = 0 m g k ( 0 + ) φ k k ! , f o r φ > 0 , a n d κ 1 < β κ , κ N . D φ β κ φ β g ( φ ) = g ( φ ) .

2.4. Definition

The Caputo operator in terms of Elzaki transform is [55]:
E [ D φ β g ( φ ) ] = s β E [ g ( φ ) ] k = 0 κ 1 s 2 β + k g ( k ) ( 0 ) , w h e r e κ 1 < β < κ .

3. Methodology of NITM

Consider fractional order PDE of the form
D ϑ β ξ ( ϕ , ϑ ) + N ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) = h ( ϕ , ϑ ) , n N , n 1 < β n ,
subjected to initial condition
ξ k ( ϕ , 0 ) = g k ( ϕ ) , k = 0 , 1 , 2 , . . . , n 1 ,
where N and M represents linear and non-linear terms.
By employing Elzaki transform for Equation (4), we obtain
E [ D ϑ β ξ ( ϕ , ϑ ) ] + E [ N ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) ] = E [ h ( ϕ , ϑ ) ] .
By employing the Elzaki differentiation property
E [ ξ ( ϕ , ϑ ) ] = k = 0 m s 2 β + k u ( k ) ( ϕ , 0 ) + s β E [ h ( ϕ , ϑ ) ] s β E [ N ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) ] ,
By applying the inverse Elzaki transform to Equation (7),
ξ ( ϕ , ϑ ) = E 1 k = 0 m s 2 β + k u k ( ϕ , 0 ) + s β E [ h ( ϕ , ϑ ) ] E 1 s β E [ N ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) ] .
By means of the iterative technique, we have
ξ ( ϕ , ϑ ) = m = 0 ξ m ( ϕ , ϑ ) .
N m = 0 ξ m ( ϕ , ϑ ) = m = 0 N ξ m ( ϕ , ϑ ) ,
the non-linear term N is determined as
N m = 0 ξ m ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + N k = 0 m ξ k ( ϕ , ϑ ) M k = 0 m ξ k ( ϕ , ϑ ) .
On putting (9)–(11) into Equation (8), we have
m = 0 ξ m ( ϕ , ϑ ) = E 1 s β k = 0 m s 2 ϕ + k u k ( ϕ , 0 ) + E [ h ( ϕ , ϑ ) ] E 1 s β E N k = 0 m ξ k ( ϕ , ϑ ) M k = 0 m ξ k ( ϕ , ϑ ) .
Hence, the iteration formula is defined as
ξ 0 ( ϕ , ϑ ) = E 1 s β k = 0 m s 2 ϕ + k u k ( ϕ , 0 ) + s β E ( g ( ϕ , ϑ ) ) ,
ξ 1 ( ϕ , ϑ ) = E 1 s β E [ N [ ξ 0 ( ϕ , ϑ ) ] + M [ ξ 0 ( ϕ , ϑ ) ] ,
ξ m + 1 ( ϕ , ϑ ) = E 1 s β E N k = 0 m ξ k ( ϕ , ϑ ) M k = 0 m ξ k ( ϕ , ϑ ) , m 1 .
Thus, the solution for the m-term in the series is obtained by means of Equations (4) and (5)
ξ ( ϕ , ϑ ) ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + . . . . . , + ξ m ( ϕ , ϑ ) , m = 1 , 2 , .

4. Methodology of HPTM

Consider the fractional order PDE of the form
D ϑ β ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) + N ξ ( ϕ , ϑ ) = h ( ϕ , ϑ ) , ϑ > 0 , 0 < β 1 , ξ ( ϕ , 0 ) = g ( ϕ ) , ϕ .
By employing the Elzaki transform in Equation (17), we obtain
E [ D ϑ β ξ ( ϕ , ϑ ) + M ξ ( ϕ , ϑ ) + N ξ ( ϕ , ϑ ) ] = E [ h ( ϕ , ϑ ) ] , ϑ > 0 , 0 < β 1 , ξ ( ϕ , ϑ ) = s 2 g ( ϕ ) + s β E [ h ( ϕ , ϑ ) ] s β E [ M ξ ( ϕ , ϑ ) + N ξ ( ϕ , ϑ ) ] .
By applying inverse Elzaki transform, we obtain
ξ ( ϕ , ϑ ) = F ( ϕ , ϑ ) E 1 s β E { M ξ ( ϕ , ϑ ) + N ξ ( ϕ , ϑ ) } ,
where
F ( ϕ , ϑ ) = E 1 s 2 g ( ϕ ) + s β E [ h ( ϕ , ϑ ) ] = g ( ϕ ) + E 1 s β E [ h ( ϕ , ϑ ) ] .
For parameter p, the perturbation technique is determined as
ξ ( ϕ , ϑ ) = k = 0 p k ξ k ( ϕ , ϑ ) ,
here p is the perturbation parameter and p [ 0 , 1 ] .
The nonlinear components are defined as
N ξ ( ϕ , ϑ ) = k = 0 p k H k ( ξ k ) ,
where He`s polynomials are represented by H n with ξ 0 , ξ 1 , ξ 2 , . . . , ξ n , and are given as
H n ( ξ 0 , ξ 1 , , ξ n ) = 1 β ( n + 1 ) D p k N k = 0 p k ξ k p = 0 ,
where D p k = k p k .
On putting Equations (22) and (23) into Equation (19), we obtain
k = 0 p k ξ k ( ϕ , ϑ ) = F ( ϕ , ϑ ) p × E 1 s β E { M k = 0 p k ξ k ( ϕ , ϑ ) + k = 0 p k H k ( ξ k ) } .
By comparing both sides of the coefficient of p, we have
p 0 : ξ 0 ( ϕ , ϑ ) = F ( ϕ , ϑ ) , p 1 : ξ 1 ( ϕ , ϑ ) = E 1 s β E ( M ξ 0 ( ϕ , ϑ ) + H 0 ( ξ ) ) , p 2 : ξ 2 ( ϕ , ϑ ) = E 1 s β E ( M ξ 1 ( ϕ , ϑ ) + H 1 ( ξ ) ) , p k : ξ k ( ϕ , ϑ ) = E 1 s β E ( M ξ k 1 ( ϕ , ϑ ) + H k 1 ( ξ ) ) , k > 0 , k N .
The terms ξ k ( ϕ , ϑ ) are easily computable giving convergent series. On taking p 1 ,
ξ ( ϕ , ϑ ) = lim M k = 1 M ξ k ( ϕ , ϑ ) .

4.1. Example

Let us consider the time-fractional Fokker–Planck equation as
β ϑ β ξ ( ϕ , ϑ ) + ϕ ϕ 6 ξ ( ϕ , ϑ ) 2 ϕ 2 ϕ 2 12 ξ ( ϕ , ϑ ) = 0 , ϕ , ϑ > 0 , β ( 0 , 1 ] ,
subjected to initial condition
ξ ( ϕ , 0 ) = ϕ 2
for special value β = 1 ; the exact solution is
ξ ( ϕ , ϑ ) = ϕ 2 exp ϑ 2
By employing the Elzaki transform in Equation (27), we have
E [ ξ ( ϕ , ϑ ) ] = s 2 ( ϕ 2 ) + s β E ϕ ϕ 6 ξ ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 12 ξ ( ϕ , ϑ ) ,
By applying the inverse Elzaki transform, we obtain
ξ ( ϕ , ϑ ) = e ϕ + E 1 s β E ϕ ϕ 6 ξ ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 12 ξ ( ϕ , ϑ ) .
Hence, by implementing NITM, we obtain
ξ 0 ( ϕ , ϑ ) = ϕ 2 , ξ 1 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ 6 ξ 0 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 12 ξ 0 ( ϕ , ϑ ) = ϕ 2 ϑ β 2 Γ ( β + 1 ) ,
ξ 2 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ 6 ξ 1 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 12 ξ 1 ( ϕ , ϑ ) = ϕ 2 ( ϑ β ) 2 4 Γ ( 2 β + 1 ) ,
ξ 3 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ 6 ξ 2 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 12 ξ 2 ( ϕ , ϑ ) = ϕ 2 ( ϑ β ) 3 8 Γ ( 3 β + 1 ) ,
Thus, we obtain solution in series form as
ξ ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + ξ 3 ( ϕ , ϑ ) + .
So we obtain
ξ ( ϕ , ϑ ) = ϕ 2 1 + ϑ β 2 Γ ( β + 1 ) + ϑ 2 β 4 Γ ( 2 β + 1 ) + ϑ 3 β 8 Γ ( 3 β + 1 ) + .
Now, by implementing HPETM, we obtain
n = 0 p n w n ( ϕ , ϑ ) = ( ϕ 2 ) + p { E 1 ( s β E [ n = 0 p n H n ( w ) ] ) } .
By comparing both sides of the coefficient of p, we obtain:
p 0 : w 0 ( ϕ , ϑ ) = ϕ 2 , p 1 : w 1 ( ϕ , ϑ ) = [ E 1 { s β E ( H 0 ( w ) ) } ] = ϕ 2 ϑ β 2 Γ ( β + 1 ) , p 2 : w 2 ( ϕ , ϑ ) = [ E 1 { s β E ( H 1 ( w ) ) } ] = ϕ 2 ( ϑ β ) 2 4 Γ ( 2 β + 1 ) , p 3 : w 3 ( ϕ , ϑ ) = [ E 1 { s β E ( H 2 ( w ) ) } ] = ϕ 2 ( ϑ β ) 3 8 Γ ( 3 β + 1 ) ,
Thus, we obtain the solution in series form in terms of HPM as
ξ ( ϕ , ϑ ) = n = 0 p n w n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ 2 1 + ϑ β 2 Γ ( β + 1 ) + ϑ 2 β 4 Γ ( 2 β + 1 ) + ϑ 3 β 8 Γ ( 3 β + 1 ) + .
The graphs in Figure 1 depict how the exact and suggested techniques solved the problem when β = 1 .  Figure 1 depicts our method’s solution at various fractional orders of β = 1 , 0.75 , 0.50 , 0.25 inside the domain of 0 ϕ , ϑ 1 , while Figure 1 depicts the solution for problem 1 at ϑ = 0.5 and 0 ϕ 1 , respectively. Additionally, Table 1 compares the proposed method results in terms of absolute error at various fractional orders.

4.2. Example

Let us consider the time-fractional Fokker–Planck equation as
ϑ β ξ ( ϕ , ϑ ) + ϕ ϕ ξ ( ϕ , ϑ ) 2 ϕ 2 ϕ 2 2 ξ ( ϕ , ϑ ) = 0 , ϕ , ϑ > 0 , β ( 0 , 1 ] ,
subjected to initial condition
ξ ( ϕ , 0 ) = ϕ ,
for special value β = 1 , the exact solution is
ξ ( ϕ , ϑ ) = ϕ exp ϑ .
By employing Elzaki transform in Equation (38), we have
E [ ξ ( ϕ , ϕ , ϑ ) ] = s 2 ( ϕ ) + s β E ϕ ϕ ξ ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 2 ξ ( ϕ , ϑ ) ,
By applying inverse Elzaki transform, we obtain
ξ ( ϕ , ϕ , ϑ ) = ( ϕ ) + E 1 s β E ϕ ϕ ξ ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 2 ξ ( ϕ , ϑ ) .
Hence, by implementing NITM, we obtain
ξ 0 ( ϕ , ϑ ) = ϕ , ξ 1 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ ξ 0 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 2 ξ 0 ( ϕ , ϑ ) = ϕ ϑ β Γ ( β + 1 ) ,
ξ 2 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ ξ 1 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 2 ξ 1 ( ϕ , ϑ ) = ϕ ( ϑ β ) 2 Γ ( 2 β + 1 ) ,
ξ 3 ( ϕ , ϑ ) = E 1 s β E ϕ ϕ ξ 2 ( ϕ , ϑ ) + 2 ϕ 2 ϕ 2 2 ξ 2 ( ϕ , ϑ ) = ϕ ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
Thus, we obtain the solution in series form as
ξ ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + ξ 3 ( ϕ , ϑ ) + ξ n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ + ϕ ϑ β Γ ( β + 1 ) + ϕ ( ϑ β ) 2 Γ ( 2 β + 1 ) + ϕ ( ϑ β ) 3 Γ ( 3 β + 1 ) + .
Now, by implementing HPETM, we obtain
n = 0 p n w n ( ϕ , ϑ ) = ( e ϕ ) + p { E 1 ( s β E [ n = 0 p n H n ( w ) ] ) } .
By comparing both sides of the coefficient of p, we obtain:
p 0 : w 0 ( ϕ , ϑ ) = ϕ , p 1 : w 1 ( ϕ , ϑ ) = [ E 1 { s β E ( H 0 ( w ) ) } ] = ϕ ϑ β Γ ( β + 1 ) , p 2 : w 2 ( ϕ , ϑ ) = [ E 1 { s β E ( H 1 ( w ) ) } ] = ϕ ( ϑ β ) 2 Γ ( 2 β + 1 ) , p 3 : w 3 ( ϕ , ϑ ) = [ E 1 { s β E ( H 2 ( w ) ) } ] = ϕ ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
Thus, we obtain solution in series form in terms of HPM as
ξ ( ϕ , ϑ ) = n = 0 p n w n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ + ϕ ϑ β Γ ( β + 1 ) + ϕ ( ϑ β ) 2 Γ ( 2 β + 1 ) + ϕ ( ϑ β ) 3 Γ ( 3 β + 1 ) + .
The graphs in Figure 2 depict how the exact and suggested techniques solved the problem when β = 1 .  Figure 2 depicts our method’s solution at various fractional orders of β = 1 , 0.75 , 0.50 , 0.25 inside the domain of 0 ϕ , ϑ 5 , while Figure 2 depicts the solution for problem 2 at ϑ = 0.5 and 0 ϕ 5 , respectively. Additionally, Table 2 compares the proposed method results in terms of absolute error at various fractional orders.

4.3. Example

Let us consider the time-fractional Fokker–Planck equation as
ϑ β ξ ( ϕ , ϑ ) + ϕ 4 ϕ ξ 2 ( ϕ , ϑ ) ϕ ϕ 3 ξ ( ϕ , ϑ ) 2 ϕ 2 ξ 2 ( ϕ , ϑ ) = 0 , ϕ , ϑ > 0 , β ( 0 , 1 ] ,
subjected to initial condition
ξ ( ϕ , 0 ) = ϕ 2 .
for special value β = 1 ; the exact solution is
ξ ( ϕ , ϑ ) = ϕ 2 exp ϑ .
By employing the Elzaki transform in Equation (49), we have
E [ ξ ( ϕ , ϕ , ϑ ) ] = s 2 ( ϕ 2 ) + s β E ϕ ϕ 3 ξ ( ϕ , ϑ ) + 2 ϕ 2 ξ 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 2 ( ϕ , ϑ ) ,
By applying the inverse Elzaki transform, we obtain
ξ ( ϕ , ϕ , ϑ ) = ϕ 2 + E 1 s β E ϕ ϕ 3 ξ ( ϕ , ϑ ) + 2 ϕ 2 ξ 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 2 ( ϕ , ϑ ) .
Hence, by implementing NITM, we obtain
ξ 0 ( ϕ , ϑ ) = ϕ 2 , ξ 1 ( ϕ , ϑ ) = E 1 [ s β E { ϕ ϕ 3 ξ 0 ( ϕ , ϑ ) + 2 ϕ 2 ξ 0 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 0 2 ( ϕ , ϑ ) } ] = ϕ 2 ϑ β Γ ( β + 1 ) ,
ξ 2 ( ϕ , ϑ ) = E 1 [ s β E { ϕ ϕ 3 ξ 1 ( ϕ , ϑ ) + 2 ϕ 2 ξ 1 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 1 2 ( ϕ , ϑ ) } ] = ϕ 2 ( ϑ β ) 2 Γ ( 2 β + 1 ) ,
ξ 3 ( ϕ , ϑ ) = E 1 [ s β E { ϕ ϕ 3 ξ 2 ( ϕ , ϑ ) + 2 ϕ 2 ξ 2 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 2 2 ( ϕ , ϑ ) } ] = ϕ 2 ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
Thus, we obtain the solution in series form as
ξ ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + ξ 3 ( ϕ , ϑ ) + ξ n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ 2 + ϕ 2 ϑ β Γ ( β + 1 ) + ϕ 2 ( ϑ β ) 2 Γ ( 2 β + 1 ) + ϕ 2 ( ϑ β ) 3 Γ ( 3 β + 1 ) + .
Now, by implementing HPETM, we obtain
n = 0 p n w n ( ϕ , ϑ ) = ( e ϕ ) + p { E 1 ( s β E [ n = 0 p n H n ( w ) ] ) } .
The non-linear terms are represented by the polynomials H n ( w ) . The elements of He’s polynomials, for example, are obtained using the recursive relationship H n ( w ) = 2 ϕ 2 ξ 2 ( ϕ , ϑ ) ϕ 4 ϕ ξ 2 ( ϕ , ϑ ) , n N . The following approximation is achieved by equating the equivalent power coefficient of p on both sides:
p 0 : w 0 ( ϕ , ϑ ) = cos ( ϕ ) , p 1 : w 1 ( ϕ , ϑ ) = [ E 1 { s β E ( H 0 ( w ) ) } ] = ϕ 2 ϑ β Γ ( β + 1 ) , p 2 : w 2 ( ϕ , ϑ ) = [ E 1 { s β E ( H 1 ( w ) ) } ] = ϕ 2 ( ϑ β ) 2 Γ ( 2 β + 1 ) , p 3 : w 3 ( ϕ , ϑ ) = [ E 1 { s β E ( H 2 ( w ) ) } ] = ϕ 2 ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
The solution in series form by means of HPM is given as
ξ ( ϕ , ϑ ) = n = 0 p n w n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ 2 + ϕ 2 ϑ β Γ ( β + 1 ) + ϕ 2 ( ϑ β ) 2 Γ ( 2 β + 1 ) + ϕ 2 ( ϑ β ) 3 Γ ( 3 β + 1 ) .
The graphs in Figure 3 depict how the exact and suggested techniques solved the problem when β = 1 .  Figure 3 depicts our method’s solution at various fractional orders of β = 1 , 0.75 , 0.50 , 0.25 inside the domain of 0 ϕ , ϑ 5 , while Figure 3 depicts the solution for problem 3 at ϑ = 0.5 and 0 ϕ 10 , respectively.

4.4. Example

Let us consider the time-fractional Fokker–Planck equation as
ϑ β ξ ( ϕ , ϑ ) ϕ ξ ( ϕ , ϑ ) 2 ϕ 2 ξ ( ϕ , ϑ ) = 0 , ϑ > 0 , β ( 0 , 1 ] ,
subjected to initial condition
ξ ( ϕ , 0 ) = ϕ ,
for special value β = 1 ; the exact solution is
ξ ( ϕ , ϑ ) = ϕ + ϑ .
By employing the Elzaki transform in Equation (60), we have
E [ ξ ( ϕ , ϕ , ϑ ) ] = s 2 ( ϕ ) + s β E ϕ ξ ( ϕ , ϑ ) + 2 ϕ 2 ξ ( ϕ , ϑ ) ,
By applying the inverse Elzaki transform, we obtain
ξ ( ϕ , ϕ , ϑ ) = ( ϕ ) + E 1 s β E ϕ ξ ( ϕ , ϑ ) + 2 ϕ 2 ξ ( ϕ , ϑ ) .
Hence, by implementing NITM, we obtain
ξ 0 ( ϕ , ϑ ) = ϕ , ξ 1 ( ϕ , ϑ ) = E 1 s β E ϕ ξ 0 ( ϕ , ϑ ) + 2 ϕ 2 ξ 0 ( ϕ , ϑ ) = ϑ β Γ ( β + 1 ) ,
ξ 2 ( ϕ , ϑ ) = E 1 s β E ϕ ξ 1 ( ϕ , ϑ ) + 2 ϕ 2 ξ 1 ( ϕ , ϑ ) = 0 ,
ξ 3 ( ϕ , ϑ ) = E 1 s β E ϕ ξ 2 ( ϕ , ϑ ) + 2 ϕ 2 ξ 2 ( ϕ , ϑ ) = 0 ,
Thus, we obtain the solution in series form as
ξ ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + ξ 3 ( ϕ , ϑ ) + ξ n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ + ϑ β Γ ( β + 1 ) .
Now, by implementing HPETM, we obtain
n = 0 p n w n ( ϕ , ϑ ) = ( e ϕ ) + p { E 1 ( s β E [ n = 0 p n H n ( w ) ] ) } .
By comparing both sides of the coefficient of p, we obtain:
p 0 : w 0 ( ϕ , ϑ ) = ϕ , p 1 : w 1 ( ϕ , ϑ ) = [ E 1 { s β E ( H 0 ( w ) ) } ] = ϑ β Γ ( β + 1 ) , p 2 : w 2 ( ϕ , ϑ ) = [ E 1 { s β E ( H 1 ( w ) ) } ] = 0 , p 3 : w 3 ( ϕ , ϑ ) = [ E 1 { s β E ( H 2 ( w ) ) } ] = 0 ,
Thus, we obtain the solution in series form in terms of HPM as
ξ ( ϕ , ϑ ) = n = 0 p n w n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ϕ + ϑ β Γ ( β + 1 ) .
The graphs in Figure 4 depict how the exact and suggested techniques solved the problem when β = 1 .  Figure 4 depicts our method’s solution at various fractional orders of β = 1 , 0.75 , 0.50 , 0.25 inside the domain of 0 ϕ , ϑ 5 , while Figure 4 depicts the solution for problem 4 at ϑ = 0.5 and 0 ϕ 5 , respectively.

4.5. Example

Let us consider the time-fractional Fokker–Planck equation as
β ϑ β ξ ( ϕ , ϑ ) ( 1 ϕ ) ϕ ξ ( ϕ , ϑ ) ( e ϑ ϕ 2 ) 2 ϕ 2 ξ ( ϕ , ϑ ) = 0 , ϑ > 0 , β ( 0 , 1 ] ,
subjected to initial condition
ξ ( ϕ , 0 ) = 1 + ϕ ,
for special value β = 1 ; the exact solution is
ξ ( ϕ , ϑ ) = exp ϑ ( 1 + ϕ ) .
By employing the Elzaki transform in Equation (71), we have
E [ ξ ( ϕ , ϕ , ϑ ) ] = s 2 ( 1 + ϕ ) + s β E ( 1 ϕ ) ϕ ξ ( ϕ , ϑ ) + ( e ϑ ϕ 2 ) 2 ϕ 2 ξ ( ϕ , ϑ ) ,
By applying the inverse Elzaki transform, we obtain
ξ ( ϕ , ϕ , ϑ ) = ( 1 + ϕ ) + E 1 s β E ( 1 ϕ ) ϕ ξ ( ϕ , ϑ ) + ( e ϑ ϕ 2 ) 2 ϕ 2 ξ ( ϕ , ϑ ) .
Hence, by implementing NITM, we obtain
ξ 0 ( ϕ , ϑ ) = 1 + ϕ , ξ 1 ( ϕ , ϑ ) = E 1 s β E ( 1 ϕ ) ϕ ξ 0 ( ϕ , ϑ ) + ( e ϑ ϕ 2 ) 2 ϕ 2 ξ 0 ( ϕ , ϑ ) = ( 1 + ϕ ) ϑ β Γ ( β + 1 ) ,
ξ 2 ( ϕ , ϑ ) = E 1 s β E ( 1 ϕ ) ϕ ξ 1 ( ϕ , ϑ ) + ( e ϑ ϕ 2 ) 2 ϕ 2 ξ 1 ( ϕ , ϑ ) = ( 1 + ϕ ) ( ϑ β ) 2 Γ ( 2 β + 1 ) ,
ξ 3 ( ϕ , ϑ ) = E 1 s β E ( 1 ϕ ) ϕ ξ 2 ( ϕ , ϑ ) + ( e ϑ ϕ 2 ) 2 ϕ 2 ξ 2 ( ϕ , ϑ ) = ( 1 + ϕ ) ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
Thus, we obtain the solution in series form as
ξ ( ϕ , ϑ ) = ξ 0 ( ϕ , ϑ ) + ξ 1 ( ϕ , ϑ ) + ξ 2 ( ϕ , ϑ ) + ξ 3 ( ϕ , ϑ ) + ξ n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ( 1 + ϕ ) + ( 1 + ϕ ) ϑ β Γ ( β + 1 ) + ( 1 + ϕ ) ( ϑ β ) 2 Γ ( 2 β + 1 ) + ( 1 + ϕ ) ( ϑ β ) 3 Γ ( 3 β + 1 ) + .
Now, by implementing HPETM, we obtain
n = 0 p n w n ( ϕ , ϑ ) = ( e ϕ ) + p { E 1 ( s β E [ n = 0 p n H n ( w ) ] ) } .
By comparing both sides of the coefficient of p, we obtain:
p 0 : w 0 ( ϕ , ϑ ) = 1 + ϕ , p 1 : w 1 ( ϕ , ϑ ) = [ E 1 { s β E ( H 0 ( w ) ) } ] = ( 1 + ϕ ) ϑ β Γ ( β + 1 ) , p 2 : w 2 ( ϕ , ϑ ) = [ E 1 { s β E ( H 1 ( w ) ) } ] = ( 1 + ϕ ) ( ϑ β ) 2 Γ ( 2 β + 1 ) , p 3 : w 3 ( ϕ , ϑ ) = [ E 1 { s β E ( H 2 ( w ) ) } ] = ( 1 + ϕ ) ( ϑ β ) 3 Γ ( 3 β + 1 ) ,
Thus, we obtain the solution in series form in terms of HPM as
ξ ( ϕ , ϑ ) = n = 0 p n w n ( ϕ , ϑ ) .
So we obtain
ξ ( ϕ , ϑ ) = ( 1 + ϕ ) + ( 1 + ϕ ) ϑ β Γ ( β + 1 ) + ( 1 + ϕ ) ( ϑ β ) 2 Γ ( 2 β + 1 ) + ( 1 + ϕ ) ( ϑ β ) 3 Γ ( 3 β + 1 ) + .
The graphs in Figure 5 depict how the exact and suggested techniques solved the problem when β = 1 .  Figure 5 depicts our method’s solution at various fractional orders of β = 1 , 0.75 , 0.50 , 0.25 inside the domain of 0 ϕ , ϑ 5 , while Figure 5 depicts the solution for problem 5 at ϑ = 0.5 and 0 ϕ 10 , respectively.

5. Conclusions

To solve the space and time-fractional Fokker–Planck equation, the new iterative approach and the homotopy perturbation method are used in this article. The two methods are particularly powerful and efficient in finding analytical and numerical solutions for a wide range of space-time fractional partial differential equations. Without employing linearization, perturbation or limiting assumptions, they give results in terms of convergent series with easily computed components. The study demonstrates that the two methodologies need less computational effort than previous methods while providing quantitatively accurate results. In all examples, the excellent agreement of numerical findings between the two approaches is also evident and notable. Finally, the proposed approaches are more efficient and solve the complexity of calculating fractional order PDE solutions, which occurs frequently in science and engineering.

Author Contributions

Methodology, F.M. and A.K.; Software, R.S. and A.A.; Investigation, F.M.; Writing—original draft, F.M. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

The author, Fatemah Mofarreh, expresses her gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R27), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. Thanks to editor, reviewers and to Khalifa University, Abu Dhabi, United Arab Emirates, for financial support in this research.

Data Availability Statement

The numerical data used to support the findings of this study are included within the article.

Acknowledgments

The author, Fatemah Mofarreh, expresses her gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R27), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. Thanks to editor, reviewers and to Khalifa University, Abu Dhabi, United Arab Emirates, for financial support in this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yépez-Martínez, H.; Gómez-Aguilar, F.; Sosa, I.O.; Reyes, J.M.; Torres-Jiménez, J. The Fengs first integral method applied to the nonlinear mKdV space-time fractional partial differential equation. Rev. Mex. Fis. 2016, 62, 310–316. [Google Scholar]
  2. Baleanu, D.; Machado, J.A.; Luo, A.C. Fractional Dynamics and Control; Springer Science and Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  3. Magin Richard, L. Fractional Calculus in Bioengineering; Begell House Redding: Chicago, IL, USA, 2006. [Google Scholar]
  4. Ellahi, R.; Alamri, S.Z.; Basit, A.; Majeed, A. Effects of MHD and slip on heat transfer boundary layer flow over a moving plate based on specific entropy generation. J. Taibah Univ. Sci. 2018, 12, 476–482. [Google Scholar] [CrossRef]
  5. Machado, J.A.T.; Silva, M.F.; Barbosa, R.S.; Jesus, I.S.; Reis, C.M.; Marcos, M.G.; Galhano, A.F. Some applications of fractional calculus in engineering. Math. Probl. Eng. 2010, 2010, 639801. [Google Scholar]
  6. Li, C.; Qian, D.; Chen, Y. On Riemann-Liouville and Caputo derivatives. Discret. Dyn. Nat. Soc. 2011, 2011, 562494. [Google Scholar] [CrossRef]
  7. Kirane, M.; Abdeljabbar, A. Nonexistence of Global Solutions of Systems of Time Fractional Differential equations posed on the Heisenberg group. Math. Methods Appl. Sci. 2022, 45, 7336–7345. [Google Scholar] [CrossRef]
  8. Al-Sawalha, M.M.; Khan, A.; Ababneh, O.Y.; Botmart, T. Fractional view analysis of Kersten-Krasil’shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Math. 2022, 7, 18334–18359. [Google Scholar] [CrossRef]
  9. Rahman, Z.; Abdeljabbar, A.; Harun-Or-Roshid; Ali, M.Z. Novel Precise Solitary Wave Solutions of Two Time Fractional Nonlinear Evolution Models via the MSE Scheme. Fractal Fract. 2022, 6, 444. [Google Scholar] [CrossRef]
  10. Cuahutenango-Barro, B.; Taneco-Hernández, M.A.; Gómez-Aguilar, J.F. On the solutions of fractional-time wave equation with memory effect involving operators with regular kernel. Chaos Solitons Fractals 2018, 115, 283–299. [Google Scholar] [CrossRef]
  11. Gómez-Aguilar, J.F.; Miranda-Hernández, M.; López-López, M.G.; Alvarado-Martínez, V.M.; Baleanu, D. Modeling and simulation of the fractional space-time diffusion equation. Commun. Nonlinear Sci. Numer. Simul. 2016, 30, 115–127. [Google Scholar] [CrossRef]
  12. Eftekhari, T.; Rashidinia, J. A new operational vector approach for time-fractional subdiffusion equations of distributed order based on hybrid functions. Math. Methods Appl. Sci. 2022, 46, 388–407. [Google Scholar] [CrossRef]
  13. Sunthrayuth, P.; Alyousef, H.A.; El-Tantawy, S.A.; Khan, A.; Wyal, N. Solving Fractional-Order Diffusion Equations in a Plasma and Fluids via a Novel Transform. J. Funct. Spaces 2022, 2022, 1899130. [Google Scholar] [CrossRef]
  14. Wang, L.; Liu, G.; Xue, J.; Wong, K. Channel Prediction Using Ordinary Differential Equations for MIMO systems. IEEE Trans. Veh. Technol. 2022, 1–9. [Google Scholar] [CrossRef]
  15. Katsikis, V.N.; Mourtas, S.D.; Simos, T.E. Zeroing Neural Network for Pseudoinversionof an Arbitrary Time-Varying Matrix Based on Singular Value Decomposition. Mathematics 2022, 10, 1208. [Google Scholar] [CrossRef]
  16. Chen, H.; Li, S. Multi-Sensor Fusion by CWT-PARAFAC-IPSO-SVM for Intelligent Mechanical Fault Diagnosis. Sensors 2022, 22, 3647. [Google Scholar] [CrossRef] [PubMed]
  17. Xu, Y.; Li, Y.; Liu, D.; Jia, W.; Huang, H. Responses of Duffing oscillator with fractional damping and random phase. Nonlinear Dyn. 2013, 74, 745–753. [Google Scholar] [CrossRef]
  18. Caputo, M. Linear models of dissipation whose Q is almost frequency independent II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  19. Ford, N.J.; Simpson, A.C. The numerical solution of fractional differential equations: Speed versus accuracy. Numer. Algorithms 2001, 26, 333–346. [Google Scholar] [CrossRef]
  20. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  21. Liu, L.; Zhang, S.; Zhang, L.; Pan, G.; Yu, J. Multi-UUV Maneuvering Counter-Game for Dynamic Target Scenario Based on Fractional-Order Recurrent Neural Network. IEEE Trans. Cybern. 2022, 1–14. [Google Scholar] [CrossRef]
  22. Shah, N.A.; Dassios, I.; Chung, J.D. A decomposition method for a fractional order multi-dimensional telegraph equation via the Elzaki transform. Symmetry 2021, 13, 8. [Google Scholar] [CrossRef]
  23. Saadeh, R.; Qazza, A.; Burqan, A. A new integral transform: ARA transform and its properties and applications. Symmetry 2020, 12, 925. [Google Scholar] [CrossRef]
  24. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Fract. Differ. Appl. 2015, 2, 731–785. [Google Scholar]
  25. Yang, X.J.; Baleanu, D.; Srivastava, H.M. Local Fractional Integral Transforms and Their Applications; Academic Press: Cambridge, MA, USA, 2015. [Google Scholar]
  26. Losada, J.; Nieto, J.J. Properties of the new fractional derivative without singular kernel. Fract. Differ. Appl. 2015, 2, 87–92. [Google Scholar]
  27. Zidan, A.M.; Khan, A.; Shah, R.; Alaoui, M.K.; Weera, W. Evaluation of time-fractional Fisher’s equations with the help of analytical methods. AIMS Math. 2022, 7, 18746–18766. [Google Scholar] [CrossRef]
  28. Kbiri Alaoui, M.; Nonlaopon, K.; Zidan, A.M.; Khan, A.; Shah, R. Analytical investigation of fractional order cahn-hilliard and gardner equations using two novel techniques. Mathematics 2022, 10, 1643. [Google Scholar] [CrossRef]
  29. Areshi, M.; Khan, A.; Shah, R.; Nonlaopon, K. Analytical investigation of fractional order Newell-Whitehead-Segel equations via a novel transform. AIMS Math. 2022, 7, 6936–6958. [Google Scholar] [CrossRef]
  30. Alyobi, S.; Shah, R.; Khan, A.; Shah, N.A.; Nonlaopon, K. Fractional Analysis of Nonlinear Boussinesq Equation under Atangana-Baleanu-Caputo Operator. Symmetry 2022, 14, 2417. [Google Scholar] [CrossRef]
  31. Fan, X.; Wei, G.; Lin, X.; Wang, X.; Si, Z.; Zhang, X.; Zhao, W. Reversible Switching of Interlayer Exchange Coupling through Atomically Thin VO2 via Electronic State Modulation. Matter 2020, 2, 1582–1593. [Google Scholar] [CrossRef]
  32. Meng, F.; Pang, A.; Dong, X.; Han, C.; Sha, X.; Jing, N.; Na, J. H-infinity Optimal Performance Design of an Unstable Plant under Bode Integral Constraint. Complexity 2018, 2018, 1–10. [Google Scholar] [CrossRef]
  33. Meng, F.; Wang, D.; Yang, P.; Xie, G.; Cutberto, R.; Romero-Melendez, C. Application of Sum of Squares Method in Nonlinear H Control for Satellite Attitude Maneuvers. Complexity 2019, 2019, 1–10. [Google Scholar] [CrossRef]
  34. Jin, H.; Wang, Z. Boundedness, blowup and critical mass phenomenon in competing chemotaxis. J. Differ. Equ. 2016, 260, 162–196. [Google Scholar] [CrossRef]
  35. Liu, P.; Shi, J.; Wang, Z.-A. Pattern formation of the attraction-repulsion Keller-Segel system. Discret. Contin. Dyn. Syst. B 2013, 18, 2597–2625. [Google Scholar] [CrossRef]
  36. He, H.M.; Peng, J.G.; Li, H.Y. Iterative approximation of fixed point problems and variational inequality problems on Hadamard manifolds. UPB Bull. Ser. A 2022, 84, 25–36. [Google Scholar]
  37. Risken, H. The FokkerPlanck Equation: Method of Solution and Applications; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  38. Chandresekhar, S. Stochastic problems in physics and astronomy. Rev. Mod. Phys. 1943, 15, 1–89. [Google Scholar] [CrossRef]
  39. Kramers, H.A. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica 1940, 7, 284–304. [Google Scholar] [CrossRef]
  40. Fokker, A. The median energy of rotating electrical dipoles in radiation fields. Annalen Der Physik 1914, 43, 810–820. [Google Scholar] [CrossRef]
  41. He, J.H.; Wu, X.H. Construction of solitary solution and compaction-like solution by variational iteration method. Chaos Solitons Fractals 2006, 29, 108–113. [Google Scholar] [CrossRef]
  42. Jumarie, G. Fractional Brownian motions via random walk in the complex plane and via fractional derivative, comparison and further results on their Fokker–Planck equations. Chaos Solitons Fractals 2004, 22, 907–925. [Google Scholar] [CrossRef]
  43. Kamitani, Y.; Matsuba, I. Self-similar characteristics of neural networks based on Fokker–Planck equation. Chaos Solitons Fractals 2004, 20, 329–335. [Google Scholar] [CrossRef]
  44. Xu, Y.; Ren, F.Y.; Liang, J.R.; Qiu, W.Y. Stretched Gaussian asymptotic behavior for fractional Fokker–Planck equation on fractal structure in external force fields. Chaos Solitons Fractals 2004, 20, 581–586. [Google Scholar] [CrossRef]
  45. Zak, M. Expectation-based intelligent control. Chaos Solitons Fractals 2006, 28, 616–626. [Google Scholar] [CrossRef]
  46. Risken, H. Fokker–Planck Equation; Springer: Berlin/Heidelberg, Germany, 1996; pp. 63–95. [Google Scholar]
  47. Herau, F. Short and long time behavior of the Fokker–Planck equation in a confining potential and applications. J. Funct. Anal. 2007, 244, 95–118. [Google Scholar] [CrossRef]
  48. Daftardar-Gejji, V.; Bhalekar, V.S. Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method. Comput. Math. Appl. 2010, 59, 1801–1809. [Google Scholar] [CrossRef]
  49. Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. 2006, 316, 753–763. [Google Scholar] [CrossRef]
  50. Bhalekar, S.; Daftardar-Gejji, V. New iterative method: Application to partial differential equations. Appl. Math. Comput. 2008, 203, 778–783. [Google Scholar] [CrossRef]
  51. Jafari, H. Iterative Methods for Solving System of Fractional Differential Equations. Ph.D. Thesis, Pune University, Pune, India, 2006. [Google Scholar]
  52. Bhalekar, S.; Daftardar-Gejji, V. Solving evolution equations using a new iterative method. Numer. Methods Partial. Differ. Equ. 2010, 26, 906–916. [Google Scholar] [CrossRef]
  53. He, J.-H. The homotopy perturbation method nonlinear oscillators with discontinuities. Appl. Math. Comput. 2004, 151, 287–292. [Google Scholar] [CrossRef]
  54. He, H. Homotopy perturbation technique. Comput. Methods Appl. Mech. Eng. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  55. Nonlaopon, K.; Alsharif, A.M.; Zidan, A.M.; Khan, A.; Hamed, Y.S.; Shah, R. Numerical investigation of fractional order Swift-Hohenberg equations via a novel transform. Symmetry 2021, 13, 1263. [Google Scholar] [CrossRef]
Figure 1. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 1.
Figure 1. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 1.
Symmetry 15 00430 g001
Figure 2. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 2.
Figure 2. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 2.
Symmetry 15 00430 g002
Figure 3. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 3.
Figure 3. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 3.
Symmetry 15 00430 g003
Figure 4. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 4.
Figure 4. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 4.
Symmetry 15 00430 g004
Figure 5. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 5.
Figure 5. Graphical layout of exact solution, proposed method solution, and 3D and 2D behavior at various fractional orders of example 5.
Symmetry 15 00430 g005
Table 1. Comparative analysis of NITM and HPTM solution of example 1.
Table 1. Comparative analysis of NITM and HPTM solution of example 1.
η ξ | Exact NITM | | Exact NITM | | Exact HPTM | | Exact HPTM |
β = 0 . 6 β = 1 β = 0 . 8 β = 1
0.58.48379150 × 10 03 1.5700000 × 10 08 4.09610800 × 10 04 1.5700000 × 10 08
11.69675830 × 10 02 6.3000000 × 10 08 1.63844300 × 10 03 6.3000000 × 10 08
1.52.54513740 × 10 02 1.4100000 × 10 07 3.68649700 × 10 03 1.4100000 × 10 07
23.39351660 × 10 02 2.5200000 × 10 07 6.55377200 × 10 03 2.5200000 × 10 07
2.54.24189580 × 10 02 3.9300000 × 10 07 1.02402690 × 10 02 3.9300000 × 10 07
0.00135.09027490 × 10 02 5.6700000 × 10 07 1.47459870 × 10 02 5.6700000 × 10 07
3.55.93865400 × 10 02 7.7000000 × 10 07 2.00709300 × 10 02 7.7000000 × 10 07
46.78703320 × 10 02 1.0100000 × 10 06 2.62150900 × 10 02 1.0100000 × 10 06
4.57.63541240 × 10 02 1.2700000 × 10 06 3.31784700 × 10 02 1.2700000 × 10 06
58.48379150 × 10 02 1.5700000 × 10 06 4.09610800 × 10 02 1.5700000 × 10 06
0.51.27076980 × 10 02 6.2500000 × 10 08 6.81188800 × 10 04 6.2500000 × 10 08
12.54153960 × 10 02 2.5000000 × 10 07 2.72475500 × 10 03 2.5000000 × 10 07
1.53.81230940 × 10 02 5.6300000 × 10 07 6.13069900 × 10 03 5.6300000 × 10 07
25.08307920 × 10 02 1.0000000 × 10 06 1.08990200 × 10 02 1.0000000 × 10 06
2.56.35384900 × 10 02 1.5630000 × 10 06 1.70297190 × 10 02 1.5630000 × 10 06
0.00237.62461880 × 10 02 2.2500000 × 10 06 2.45227950 × 10 02 2.2500000 × 10 06
3.58.89538860 × 10 02 3.0600000 × 10 06 3.33782500 × 10 02 3.0600000 × 10 06
41.01661584 × 10 01 4.0000000 × 10 06 4.35960800 × 10 02 4.0000000 × 10 06
4.51.14369282 × 10 01 5.0600000 × 10 06 5.51762900 × 10 02 5.0600000 × 10 06
51.2707698 × 10 01 6.2500000 × 10 06 6.81188800 × 10 02 6.2500000 × 10 06
Table 2. Comparative analysis of NITM and HPTM solution of example 2.
Table 2. Comparative analysis of NITM and HPTM solution of example 2.
η ξ | Exact NITM | | Exact NITM | | Exact HPTM | | Exact HPTM |
β = 0 . 6 β = 1 β = 0 . 8 β = 1
0.54.36677680 × 10 03 3.5000000 × 10 09 1.64248450 × 10 03 3.5000000 × 10 09
11.74671070 × 10 02 7.0000000 × 10 09 3.28496900 × 10 03 7.0000000 × 10 09
1.53.93009910 × 10 02 1.0000000 × 10 08 4.92745400 × 10 03 1.0000000 × 10 08
26.98684280 × 10 02 1.4000000 × 10 08 6.56993800 × 10 03 1.4000000 × 10 08
2.51.09169419 × 10 02 1.8000000 × 10 08 8.21242200 × 10 03 1.8000000 × 10 08
0.00131.57203963 × 10 02 2.1000000 × 10 08 9.85490700 × 10 03 2.1000000 × 10 08
3.52.13972060 × 10 01 2.4000000 × 10 08 1.14973920 × 10 02 2.4000000 × 10 08
42.79473710 × 10 01 2.8000000 × 10 08 1.31398760 × 10 02 2.8000000 × 10 08
4.53.53708920 × 10 01 3.2000000 × 10 08 1.47823600 × 10 02 3.2000000 × 10 08
54.36677680 × 10 01 3.5000000 × 10 08 1.64248450 × 10 02 3.5000000 × 10 08
0.56.12090700 × 10 03 1.7000000 × 10 08 2.73690750 × 10 03 1.7000000 × 10 08
12.44836280 × 10 02 3.4000000 × 10 08 5.47381500 × 10 03 3.4000000 × 10 08
1.55.50881630 × 10 02 5.1000000 × 10 08 8.21072200 × 10 03 5.1000000 × 10 08
29.79345120 × 10 02 6.8000000 × 10 08 1.09476300 × 10 02 6.8000000 × 10 08
2.51.53022675 × 10 02 8.5000000 × 10 08 1.36845380 × 10 02 8.5000000 × 10 08
0.00232.20352652 × 10 01 1.0200000 × 10 07 1.64214450 × 10 02 1.0200000 × 10 07
3.52.99924450 × 10 01 1.1900000 × 10 07 1.91583520 × 10 02 1.1900000 × 10 07
43.91738050 × 10 01 1.3600000 × 10 07 2.18952600 × 10 02 1.3600000 × 10 07
4.54.95793470 × 10 01 1.5300000 × 10 07 2.46321680 × 10 02 1.5300000 × 10 07
56.12090700 × 10 01 1.7000000 × 10 07 2.73690750 × 10 02 1.7000000 × 10 07
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

Mofarreh, F.; Khan, A.; Shah, R.; Abdeljabbar, A. A Comparative Analysis of Fractional-Order Fokker–Planck Equation. Symmetry 2023, 15, 430. https://doi.org/10.3390/sym15020430

AMA Style

Mofarreh F, Khan A, Shah R, Abdeljabbar A. A Comparative Analysis of Fractional-Order Fokker–Planck Equation. Symmetry. 2023; 15(2):430. https://doi.org/10.3390/sym15020430

Chicago/Turabian Style

Mofarreh, Fatemah, Asfandyar Khan, Rasool Shah, and Alrazi Abdeljabbar. 2023. "A Comparative Analysis of Fractional-Order Fokker–Planck Equation" Symmetry 15, no. 2: 430. https://doi.org/10.3390/sym15020430

APA Style

Mofarreh, F., Khan, A., Shah, R., & Abdeljabbar, A. (2023). A Comparative Analysis of Fractional-Order Fokker–Planck Equation. Symmetry, 15(2), 430. https://doi.org/10.3390/sym15020430

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