Next Article in Journal
Causal Information Rate
Next Article in Special Issue
Trapezoidal-Type Inequalities for Strongly Convex and Quasi-Convex Functions via Post-Quantum Calculus
Previous Article in Journal
Dirichlet Polynomials and Entropy
Previous Article in Special Issue
Some Parameterized Quantum Midpoint and Quantum Trapezoid Type Inequalities for Convex Functions with Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Technique, Based on Natural Transform to Solve Fractional-Order Parabolic Equations

1
Department of Mathematics, Texas A & M University-Kingsville, Kingsville, TX 78363, USA
2
Department of Mathematical Sciences College of Sciences, Princess Nourah Bint Abdulrahman University, Riyadh 11546, Saudi Arabia
3
Department of Mathematics, Abdul Wali khan University, Mardan 23200, Pakistan
4
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
*
Author to whom correspondence should be addressed.
Submission received: 24 July 2021 / Revised: 17 August 2021 / Accepted: 17 August 2021 / Published: 21 August 2021
(This article belongs to the Special Issue Advanced Numerical Methods for Differential Equations)

Abstract

:
This research article is dedicated to solving fractional-order parabolic equations using an innovative analytical technique. The Adomian decomposition method is well supported by natural transform to establish closed form solutions for targeted problems. The procedure is simple, attractive and is preferred over other methods because it provides a closed form solution for the given problems. The solution graphs are plotted for both integer and fractional-order, which shows that the obtained results are in good contact with the exact solution of the problems. It is also observed that the solution of fractional-order problems are convergent to the solution of integer-order problem. In conclusion, the current technique is an accurate and straightforward approximate method that can be applied to solve other fractional-order partial differential equations.

1. Introduction

The present research work is dedicated to studying the analytical solution of fractional-order parabolic equations. In the literature, it is well recognized that a broad range of problems in physics, engineering, nuclear physics, and mathematics can be defined as unique boundary and initial value problems. A homogeneous beam’s transverse vibrations are controlled by fractional single fourth-order parabolic partial differential equations (PDEs). Such problem types occur in viscoelastic and inelastic flow mathematical modeling, layer deflection theories, and beam deformation [1,2,3,4,5,6,7,8]. Analyses of these problems have drawn the attention of several physicists and mathematicians.
The time fractional parabolic PDEs with variable coefficient:
γ υ η γ + κ ( ξ , ϕ , ψ ) 4 υ ξ 4 + 1 ϕ μ ( ξ , ϕ , ψ ) 4 υ ϕ 4 + 1 ψ ρ ( ξ , ϕ , ψ ) 4 υ ψ 4 = g ( ξ , ϕ , ψ , η ) , 1 < γ 2 , η 0 ,
where κ ( ξ , ϕ , ψ ) , μ ( ξ , ϕ , ψ ) and ρ ( ξ , ϕ , ψ ) are positive, with initial condition
υ ( ξ , ϕ , ψ , η ) = f 0 ( ξ , ϕ , ψ ) , υ η ( ξ , ϕ , ψ , η ) = k 0 ( ξ , ϕ , ψ ) ,
with boundary conditions
υ ( a , ϕ , ψ , η ) = h 0 ( ϕ , ψ , η ) , υ ( b , ϕ , ψ , η ) = h 1 ( ϕ , ψ , η ) , υ ( ξ , a , ψ , η ) = g 0 ( ϕ , ψ , η ) , υ ( ξ , b , ψ , η ) = g 1 ( ϕ , ψ , η ) , υ ( ξ , ϕ , a , η ) = k 0 ( ϕ , ψ , η ) , υ ( ξ , ϕ , b , η ) = k 1 ( ϕ , ψ , η ) , υ ξ ξ ( a , ϕ , ψ , η ) = h 0 ¯ ( ϕ , ψ , η ) , υ ξ ξ ( b , ϕ , ψ , η ) = h 1 ¯ ( ϕ , ψ , η ) , υ ϕ ϕ ( ξ , a , ψ , η ) = g 0 ¯ ( ϕ , ψ , η ) , υ ϕ ϕ ( ξ , b , ψ , η ) = g 1 ¯ ( ϕ , ψ , η ) , υ ψ ψ ( ξ , ϕ , a , η ) = k 0 ¯ ( ϕ , ψ , η ) , υ ψ ψ ( ξ , ϕ , b , η ) = k 1 ¯ ( ϕ , ψ , η ) .
for which h , g , k , h , g , k are continuous variables and differs between 0 and 1, which is the beam’s flexural stiffness ratio [1] in its volume per unit mass, as mentioned in [1,2,3,4,5,6,7]. After being used for the first time in underwater acoustics, the parabolic equation has undergone extensive development, including improvements in accuracy and implementation in the time domain. With the introduction of the wide-angle parabolic equation, the phase errors of parabolic equation solutions, which approximate the solution of the wave equation, were greatly reduced. While various generalisations of the wide-angle parabolic equation have been considered, the parabolic equation’s aperture limitation has remained a source of concern. The time-domain parabolic equation enables one to calculate pulse propagation without using Fourier synthesis. The time-domain parabolic equation has been expanded to account for interface conditions, nonlinear propagation, density variations and sediment attenuation, as well as wide-angle diffraction and sediment dispersion. Many researchers [9,10] have attempted to study the analytical solutions of parabolic equations of the fourth order. Different techniques have been suggested recently, such as the B-spline method [11], the decomposition method [12], the implicit scheme [13] and the Spline method [14], to analyze the solution of the partial differential fourth-order parabolic equation. Biazar and Ghazvini [15] used He’s iterative technique for the solution of parabolic PDEs. The modified version of this method was introduced in [16] to solve singular fourth-order parabolic PDEs. The fourth-order parabolic PDE analytical solution was examined in [17]. The modified Laplace variational iteration technique was discussed by [18] to solve singular fourth-order parabolic PDEs.
Rawashdeh and Maitama developed a new method, which was named the natural transform decomposition method (NTDM) in 2014, to handle linear and non-linear PDEs and ODEs that occur in several applications of mathematical engineering and physics [19]. NDM is a combination of NTM [20] ADM [21]. The suggested method provides a series from a solution which converges quickly to an exact solution in a closed form, see Belgacem and Silambarasan [22]. The NTDM provides analytical results of fractional-order heat and wave problems [23]. The NTDM provides an analytical solution by using fractional-order delay PDEs [24]. Different linear and non-linear PDEs and ODEs, such as fractional diffusion equations, are solve by using NTDM [25], fractional non-linear systems of PDEs [26], fractional telegraph equation [27], and time-fractional coupled Burger equations [28].

2. Preliminaries

Definition 1.
The following transformation of f ¯ ( η ) is called natural transformation and is expressed as [29,30]
N + [ f ¯ ( η ) ] = Q 1 ( s , u ) = 1 u 0 e s η u f ¯ ( η ) d η ; s , u > 0 ,
where the transformation parameters are s and u.
Definition 2.
The following transformation of f ¯ ( η ) is called inverse natural transformation and is expressed as
N [ Q 1 ( s , u ) ] = f ¯ ( η ) = 1 2 π i p 1 i p 1 + i e s η u Q 1 ( s , u ) d s ,
where s and u denote the transformation factors and s = p 1 in the complex plane s = ξ 1 + i ϕ 1 is taken as an integral.
Definition 3.
The nth derivative in term of NT
The nth derivative in term of NT f ¯ ( η ) is the f ¯ ( η ) and is defined as
N [ f ¯ ( η ) ] = Q ( s , u ) = s u Q 1 ( s , u ) k = 0 1 s ( k + 1 ) u k f ¯ k ( 0 ) , 1 .
Theorem 1.
If H 1 ( s , u ) , L 1 ( s , u ) is the natural transformation of the corresponding functions h 1 ( η ) and l 1 ( η ) are both set to A, then H 1 ( s , u ) is the natural transformation
N [ h 1 l 1 ] = u H 1 ( s , u ) L 1 ( s , u ) ,
a where h 1 l 1 represents the convolution of h 1 and l 1 .
Definition 4.
The Riemann–Liouville fractional-order integral [31,32]
I ξ γ f ¯ ( ξ ) = f ¯ ( ξ ) if γ = 0 , 1 Γ ( γ ) 0 ξ ( ξ υ ) γ 1 f ¯ ( υ ) d υ if γ > 0 ,
where Γ is a function defined by
Γ ( Γ ) = 0 e ξ ξ ω 1 d ξ ω C .
Definition 5.
The Caputo fractional derivative operator with order γ is defined as [33]
D γ f ¯ ( ξ ) = γ f ¯ ( ξ ) η γ = I γ γ f ¯ ( ξ ) η γ , if 1 < γ , N . γ f ¯ ( ξ ) η γ .
where N , ξ > 0 , f ¯ C η , η 1 .
Definition 6.
Natural transform of D η γ u ( η ) by means of Caputo–Fabrizio is defined as
N [ D η γ ] = 1 1 γ + γ ( u s ) N [ υ ( η ) ] 1 s υ ( 0 ) .

3. Idea of FNTM

The general fractional-order PDEs are given as
C F D γ υ ( ξ , η ) + L υ ( ξ , η ) + N υ ( ξ , η ) = q ( ξ , η ) , ξ , η 0 , m 1 < γ < m ,
where D γ = γ η γ represents the fractional derivative in term of Caputo sense. Moreover, L is the linear and N is the non-linear term in Equation (1).
The initial condition is
υ ( ξ , 0 ) = k ( ξ ) , 0 < γ 1 , η > 0 .
Applying the natural transformation to Equation (1), we get
N + D γ υ ( ξ , η ) + N + L υ ( ξ , η ) + N υ ( ξ , η ) = N + q ( ξ , η ) ,
Using the natural transform’s differentiation property, we get
N + υ ( ξ , η ) = 1 s υ ( ξ , 0 ) + u ( s γ ( s u ) ) s 4 N + q ( ξ , η ) u ( s γ ( s u ) ) s 4 N + L υ ( ξ , η ) + N υ ( ξ , η ) .
Now υ ( ξ , 0 ) = k ( ξ )
N + υ ( ξ , η ) = k ( ξ ) s + u ( s γ ( s u ) ) s 4 N + q ( ξ , η ) u ( s γ ( s u ) ) s 4 N + L υ ( ξ , η ) + N υ ( ξ , η ) .
The following infinite series represents the NTDM solution υ ( ξ , η )
υ ( ξ , η ) = j = 0 υ j ( ξ , η ) ,
and Adomian polynomials as
N υ ( ξ , η ) = j = 0 A j ,
A j = 1 j ! d j d λ j N j = 0 ( λ j υ j ) λ = 0 , j = 0 , 1 , 2
We get replacement Equations (5) and (6) in Equation (4).
N + j = 0 υ j ( ξ , η ) = k ( ξ ) s + u ( s γ ( s u ) ) s 4 N + q ( ξ , η ) u ( s γ ( s u ) ) s 4 N + L j = 0 υ j ( ξ , η ) + j = 0 A j .
Applying the natural transformation’s linearity,
N + υ 0 ( ξ , η ) = k ( ξ ) s + u ( s γ ( s u ) ) s 4 N + q ( ξ , η ) , ( 8 a )
N + υ 1 ( ξ , η ) = u ( s γ ( s u ) ) s 4 N + L υ 0 ( ξ , η ) + A 0 .
We can generally write
N + υ j + 1 ( ξ , η ) = u ( s γ ( s u ) ) s 4 N + L υ j ( ξ , η ) + A j , j 1 .
Equations (9) and (10) implementing the inverse natural transformation
υ 0 ( ξ , η ) = k ( ξ ) + N u ( s γ ( s u ) ) s 4 N + q ( ξ , η ) ,
υ j + 1 ( ξ , η ) = N u ( s γ ( s u ) ) s 4 N + L υ j ( ξ , η ) + A j .

4. Numerical Implementation

4.1. Problem

Consider fractional-order one-dimensional parabolic equation:
γ + 1 υ η γ + 1 + 1 ξ + ξ 4 120 4 υ ξ 4 = 0 , 0 < γ 1 , η 0 ,
with initial condition
υ ( ξ , 0 ) = 0 , υ η ( ξ , 0 ) = 1 + ξ 5 120 ,
with boundary conditions
υ ( 1 2 , η ) = 1 + ( 1 2 ) 5 120 sin ( η ) , υ ( 1 , η ) = 121 120 sin ( η ) , 2 υ ξ 2 1 2 , η = 1 6 1 2 3 sin ( η ) , 2 υ ξ 2 ( 1 , η ) = 1 6 sin ( η ) .
Concerning the natural transformation of (12), we get
υ ( ξ , s , u ) = 1 s ( 0 ) + u s 2 1 + ξ 5 120 u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ ξ 4 .
Using the inverse natural transformation,
υ ( ξ , η ) = N u s 2 1 + ξ 5 120 u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ ξ 4 ,
υ ( ξ , η ) = 1 + ξ 5 120 η N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ ξ 4 .
The Equation (15) correction function is provided by
= 0 υ + 1 ( ξ , η ) = 1 + ξ 5 120 η N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 = 0 4 υ ξ 4 ,
The first term
υ 0 ( ξ , η ) = 1 + ξ 5 120 η ,
Then we got
υ + 1 ( ξ , η ) = N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 = 0 4 υ ξ 4 ,
for j = 0
υ 1 ( ξ , η ) = N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ 0 ξ 4 , υ 1 ( ξ , η ) = 1 + ξ 5 120 η 2 3 ! ( 3 3 γ + γ η ) .
The following terms are
υ 2 ( ξ , η ) = N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ 1 ξ 4 = 1 + ξ 5 120 η 4 5 ! ( γ η + 5 5 γ ) , υ 3 ( ξ , η ) = N u ( s γ ( s u ) ) s 4 N + 1 ξ + ξ 4 120 4 υ 2 ξ 4 = 1 + ξ 5 120 η 6 7 ! ( γ η + 7 7 γ ) ,
The series form the solution of Problems (4.1), such as:
υ ( ξ , η ) = υ 0 ( ξ , η ) + υ 1 ( ξ , η ) + υ 2 ( ξ , η ) + υ 3 ( ξ , η ) + υ 4 ( ξ , η ) .
υ ( ξ , η ) = 1 + ξ 5 120 η η 2 3 ! ( 3 3 γ + γ η ) + η 4 5 ! ( γ η + 5 5 γ ) η 6 7 ! ( γ η + 7 7 γ ) + ,
When γ = 1 , the integer NDM solution is
υ ( ξ , η ) = 1 + ξ 5 120 η η 3 3 ! + η 5 5 ! η 7 7 ! + η 9 9 ! .
The exact solution is
υ ( ξ , η ) = 1 + ξ 5 120 sin ( η ) .
Figure 1, show that the exact and analytical solution graph of Problem 4.1. In Figure 2, the obtained solutions of Problem 4.1 are plotted at various fractional orders of the derivatives; it is confirmed that the exact and derived results are in close contact with each other. Thus the proposed method provided an accurate solution for Problem 4.1.

4.2. Problem

Consider fractional-order two-dimensional parabolic equation:
γ + 1 υ η γ + 1 + 2 1 ξ 2 + ξ 4 6 ! 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ ϕ 4 = 0 , 0 < γ 1 , η 0 ,
with initial condition
υ ( ξ , ϕ , 0 ) = 0 , υ η ( ξ , ϕ , 0 ) = 2 + ξ 6 6 ! + ϕ 6 6 ! ,
with boundary conditions
υ ( 1 2 , ϕ , η ) = 2 + ( 1 2 ) 6 6 ! + ϕ 6 6 ! sin ( η ) , υ 1 2 , ϕ , η = 2 + ( 1 ) 6 6 ! + ϕ 6 6 ! sin ( η ) , υ ξ ξ ( 1 2 , ϕ , η ) = ( 1 2 ) 4 4 ! sin ( η ) , υ ξ ξ ( 1 2 , ϕ , η ) = 1 24 sin ( η ) , υ ϕ ϕ ( ξ , 1 2 , η ) = ( 1 2 ) 4 4 ! sin ( η ) , υ ϕ ϕ ( ξ , 1 2 , η ) = 1 24 sin ( η ) .
Concerning the natural transformation of (22), we get
υ ( ξ , ϕ , s , u ) = 1 s ( 0 ) + u s 2 2 + ξ 6 6 ! + ϕ 6 6 ! u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ ϕ 4 ,
using inverse natural transformation.
υ ( ξ , ϕ , η ) = N u s 2 2 + ξ 6 6 ! + ϕ 6 6 ! u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ ϕ 4 ,
υ ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ ϕ 4
The Equation (25) correction function is provided by
= 0 υ + 1 ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! = 0 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! = 0 4 υ ϕ 4 ,
The first term being
υ 0 ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η ,
Then we get
υ + 1 ( ξ , ϕ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! = 0 4 υ ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! = 0 4 υ ϕ 4 ,
for j = 0
υ 1 ( ξ , ϕ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ 0 ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ 0 ϕ 4 , υ 1 ( ξ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η 2 3 ! ( 3 3 γ + γ η ) .
The following terms are
υ 2 ( ξ , ϕ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ 1 ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ 1 ϕ 4 , υ 2 ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η 4 5 ! ( γ η + 5 5 γ ) υ 3 ( ξ , ϕ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 1 ξ 2 + ξ 4 6 ! 4 υ 2 ξ 4 + 2 1 ϕ 2 + ϕ 4 6 ! 4 υ 2 ϕ 4 , υ 3 ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η 6 7 ! ( γ η + 7 7 γ ) ,
The series forms a solution to Problems (4.2), for example,
υ ( ξ , ϕ , η ) = υ 0 ( ξ , ϕ , η ) + υ 1 ( ξ , ϕ , η ) + υ 2 ( ξ , ϕ , η ) + υ 3 ( ξ , ϕ , η ) + υ 4 ( ξ , ϕ , η ) .
υ ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η η 2 3 ! ( 3 3 γ + γ η ) + η 4 5 ! ( γ η + 5 5 γ ) η 6 7 ! ( γ η + 7 7 γ ) + .
Then γ = 1 ; the integer NDM results as
υ ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! η η 3 3 ! + η 5 5 ! η 7 7 ! + η 9 9 ! ,
The exact solution is
υ ( ξ , ϕ , η ) = 2 + ξ 6 6 ! + ϕ 6 6 ! sin ( η ) .
Figure 3 shows the exact and analytical solution grpah of Problem 4.2. In Figure 4, the obtained solutions of Problem 4.2 are plotted at various fractional orders of the derivatives; it is confirmed that the exact and derived results are in close contact with each other. Thus the proposed method provided an accurate solution for Problem 4.2.

4.3. Problem

Consider fractional-order three-dimensional parabolic equation:
γ + 1 υ η γ + 1 + 2 ϕ + ψ 2 cos ξ 1 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ ψ 4 = 0 , 0 < γ 1 , η 0 ,
with initial condition
υ ( ξ , ϕ , ψ , 0 ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) , υ η ( ξ , ϕ , ψ , 0 ) = cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( ξ + ϕ + ψ ) ,
with boundary conditions
υ ( 0 , ϕ , ψ , η ) = 1 + ϕ + ψ cos ( ϕ ) cos ( ψ ) e η , υ ( π 3 , ϕ , ψ , η ) = 2 π 3 6 + ϕ + ψ cos ( ϕ ) cos ( ψ ) e η , υ ( ξ , 0 , ψ , η ) = 1 + ξ + ψ cos ( ξ ) cos ( ψ ) e η , υ ( ξ , π 3 , ψ , η ) = 2 π 3 6 + ξ + ψ cos ( ξ ) cos ( ψ ) e η , υ ( ξ , ϕ , 0 , η ) = 1 + ξ + ϕ cos ( ξ ) cos ( ϕ ) e η , υ ( ξ , ϕ , π 3 , η ) = 2 π 3 6 + ξ + ϕ cos ( ξ ) cos ( ϕ ) e η , υ ξ ( 0 , ϕ , ψ , η ) = υ ϕ ( ξ , 0 , ψ , η ) = υ ψ ( ξ , ϕ , 0 , η ) = e η , υ ξ ( π 3 , ϕ , ψ , η ) = υ ϕ ( ξ , π 3 , ψ , η ) = υ ψ ( ξ , ϕ , π 3 , η ) = 3 + 2 2 e η .
Concerning the natural transformation of (32), we get
υ ( ξ , ϕ , ψ , s , u ) = 1 s ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) + u s 2 cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( ξ + ϕ + ψ ) u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ ψ 4 ,
using the inverse natural transform.
υ ( ξ , ϕ , ψ , η ) = N 1 s ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) + u s 2 cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( ξ + ϕ + ψ ) N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ ψ 4 ,
υ ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( 1 η ) N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ ψ 4 ,
The Equation (35) correction function is provided by
= 0 υ + 1 ( ξ , ϕ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( 1 η ) N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 = 0 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 = 0 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 = 0 4 υ ψ 4 ,
The first term being
υ 0 ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( 1 η ) ,
Then we get
υ + 1 ( ξ , ϕ , ψ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 = 0 4 υ ξ 4 + 2 ξ + ψ 2 cos ϕ 1 = 0 4 υ ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 = 0 4 υ ψ 4 ,
for j = 0
υ 1 ( ξ , ϕ , ψ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ 0 ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ 0 ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ 0 ψ 4 , υ 1 ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) ( 1 γ + γ η ) η 2 3 ! ( 3 3 γ + γ η ) .
The following terms are
υ 2 ( ξ , ϕ , ψ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ 1 ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ 1 ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ 1 ψ 4 , υ 2 ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) η 3 4 ! ( 4 4 γ + γ η ) η 4 5 ! ( γ η + 5 5 γ ) , υ 3 ( ξ , ϕ , ψ , η ) = N u ( s γ ( s u ) ) s 4 N + 2 ϕ + ψ 2 cos ξ 1 4 υ 2 ξ 4 + 2 ξ + ψ 2 cos ϕ 1 4 υ 2 ϕ 4 + 2 ϕ + ξ 2 cos ψ 1 4 υ 2 ψ 4 , υ 3 ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) η 5 6 ! ( 6 6 γ + γ η ) η 6 7 ! ( γ η + 7 7 γ ) .
The series forms a solution to Problems (4.3), for example,
υ ( ξ , ϕ , ψ , η ) = υ 0 ( ξ , ϕ , ψ , η ) + υ 1 ( ξ , ϕ , ψ , η ) + υ 2 ( ξ , ϕ , ψ , η ) + υ 3 ( ξ , ϕ , ψ , η ) + υ 4 ( ξ , ϕ , ψ , η ) .
υ ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) { 1 η + ( 1 γ + γ η ) η 2 3 ! ( 3 3 γ + γ η ) + η 3 4 ! ( 4 4 γ + γ η ) η 4 5 ! ( γ η + 5 5 γ ) + η 5 6 ! ( 6 6 γ + γ η ) η 6 7 ! ( γ η + 7 7 γ ) } .
Then γ = 1 ; the integer NDM results as
υ ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) { 1 η + η 2 2 ! η 3 3 ! + η 4 4 ! η 5 5 ! + η 6 6 ! η 7 7 ! } .
The exact solution is
υ ( ξ , ϕ , ψ , η ) = ξ + ϕ + ψ cos ( ξ ) + cos ( ϕ ) + cos ( ψ ) e η .
Figure 5, show that the exact and analytical solution graph of Problem 4.3. In Figure 6, the obtained solutions of Problem 4.3 are plotted at various fractional orders of the derivatives; it is confirmed that the exact and derived results are in close contact with each other. Thus, the proposed method provided an accurate solution for Problem 4.3.

5. Conclusions

In the present article, an efficient analytical technique is used to solve fractional-order parabolic equations. The present method is the combinations of two well-known methods, namely the natural transform and Adomian decomposition method. The natural transform is applied to the given problem, which makes it easier. After this, we implemented the Adomian decomposition method and then the inverse natural transform to get the closed form analytical solutions for the given problems. The proposed method requires a small number of calculation to attain closed form solutions and is therefore considered to be one of the best analytical techniques to solve fractional-order partial differential equations.

Author Contributions

Data curation, R.P.A.; Funding acquisition, K.N.; Investigation, F.M.; Methodology, R.S.; Resources, W.L. and K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.

Conflicts of Interest

The authors have no conflict of interest.

References

  1. Khaliq, A.Q.M.; Twizell, E.H. A family of second order methods for variable coefficient fourth order parabolic partial differential equations. Int. J. Comput. Math. 1987, 23, 63–76. [Google Scholar] [CrossRef]
  2. Gorman, D.J. Free Vibration Analysisi of Beams and Shafts(Book); Research supported by the National Research Council of Canada; Wiley-Interscience: New York, NY, USA, 1975; 395p. [Google Scholar]
  3. Andrade, C.; McKee, S. High accuracy ADI methods for fourth order parabolic equations with variable coefficients. J. Comput. Appl. Math. 1977, 3, 11–14. [Google Scholar] [CrossRef] [Green Version]
  4. Conte, S.D. A stable implicit finite difference approximation to a fourth order parabolic equation. J. ACM 1957, 4, 18–23. [Google Scholar] [CrossRef]
  5. Royster, W.C.; Conte, S.D. Convergence of finite difference solutions to a solution of the equation of the vibrating rod. Proc. Am. Math. Soc. 1956, 7, 742–749. [Google Scholar] [CrossRef]
  6. Evans, D.J. A stable explicit method for the finite-difference solution of a fourth-order parabolic partial differential equation. Comput. J. 1965, 8, 280–287. [Google Scholar] [CrossRef] [Green Version]
  7. Evans, D.J.; Yousif, W.S. A note on solving the fourth order parabolic equation by the AGE method. Int. J. Comput. Math. 1991, 40, 93–97. [Google Scholar] [CrossRef]
  8. Wazwaz, A.-M. On the solution of the fourth order parabolic equation by the decomposition method. Int. J. Comput. Math. 1995, 57, 213–217. [Google Scholar] [CrossRef]
  9. Liao, W.; Zhu, J.; Khaliq, A.Q.M. An efficient high-order algorithm for solving systems of reaction-diffusion equations. Numer. Methods Part. Differ. Equ. 2002, 18, 340–354. [Google Scholar] [CrossRef]
  10. Jain, M.K.; Iyengar, S.R.K.; Lone, A.G. Higher order difference formulas for a fourth order parabolic partial differential equation. Int. J. Numer. Methods Eng. 1976, 10, 1357–1367. [Google Scholar] [CrossRef]
  11. Caglar, H.; Caglar, N. Fifth-degree B-spline solution for a fourth-order parabolic partial differential equations. Appl. Math. Comput. 2008, 201, 597–603. [Google Scholar] [CrossRef]
  12. Wazwaz, A.-M. Analytic treatment for variable coefficient fourth-order parabolic partial differential equations. Appl. Math. Comput. 2001, 123, 219–227. [Google Scholar] [CrossRef]
  13. Rashidinia, J.; Mohammadi, R. Sextic spline solution of variable coefficient fourth-order parabolic equations. Int. J. Comput. Math. 2010, 87, 3443–3454. [Google Scholar] [CrossRef]
  14. Aziz, T.; Khan, A.; Rashidinia, J. Spline methods for the solution of fourth-order parabolic partial differential equations. Appl. Math. Comput. 2005, 167, 153–166. [Google Scholar] [CrossRef]
  15. Biazar, J.; Ghazvini, H. Hes variational iteration method for fourth-order parabolic equations. Comput. Math. Appl. 2007, 54, 1047–1054. [Google Scholar] [CrossRef] [Green Version]
  16. Noor, M.A.; Noor, K.I.; Mohyud-Din, S.T. Modified variational iteration technique for solving singular fourth-order parabolic partial differential equations. Nonlinear Anal. Theory Methods Appl. 2009, 71, e630–e640. [Google Scholar] [CrossRef]
  17. Dehghan, M.; Manafian, J. The solution of the variable coefficients fourth-order parabolic partial differential equations by the homotopy perturbation method. Z. Naturforschung A 2009, 64, 420–430. [Google Scholar] [CrossRef]
  18. Nadeem, M.; Li, F.; Ahmad, H. Modified Laplace variational iteration method for solving fourth-order parabolic partial differential equation with variable coefficients. Comput. Math. Appl. 2019, 78, 2052–2062. [Google Scholar] [CrossRef]
  19. Rawashdeh, M.; Maitama, S. Finding exact solutions of nonlinear PDEs using the natural decomposition method. Math. Methods Appl. Sci. 2017, 40, 223–236. [Google Scholar] [CrossRef]
  20. Baskonus, H.M.; Bulut, H.; Pandir, Y. The natural transform decomposition method for linear and nonlinear partial differential equations. Math. Eng. Sci. Aerosp. (MESA) 2014, 5, 111–126. [Google Scholar]
  21. Adomian, G. A new approach to nonlinear partial differential equations. J. Math. Anal. Appl. 1984, 102, 420–434. [Google Scholar] [CrossRef] [Green Version]
  22. Belgacem, F.B.M.; Silambarasan, R. Theory of natural transform. J. MESA 2012, 3, 99–124. [Google Scholar]
  23. Khan, H.; Shah, R.; Kumam, P.; Arif, M. Analytical Solutions of Fractional-Order Heat and Wave Equations by the Natural Transform Decomposition Method. Entropy 2019, 21, 597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Shah, R.; Khan, H.; Kumam, P.; Arif, M.; Baleanu, D. Natural Transform Decomposition Method for Solving Fractional-Order Partial Differential Equations with Proportional Delay. Mathematics 2019, 7, 532. [Google Scholar] [CrossRef] [Green Version]
  25. Shah, R.; Khan, H.; Mustafa, S.; Kumam, P.; Arif, M. Analytical Solutions of Fractional-Order Diffusion Equations by Natural Transform Decomposition Method. Entropy 2019, 21, 557. [Google Scholar] [CrossRef] [Green Version]
  26. Rawashdeh, M.S.; Al-Jammal, H. New approximate solutions to fractional nonlinear systems of partial differential equations using the FNDM. Adv. Differ. Equ. 2016, 1, 235. [Google Scholar] [CrossRef] [Green Version]
  27. Eltayeb, H.; Abdalla, Y.T.; Bachar, I.; Khabir, M.H. Fractional telegraph equation and its solution by natural transform decomposition method. Symmetry 2019, 11, 334. [Google Scholar] [CrossRef] [Green Version]
  28. Prakasha, D.G.; Veeresha, P.; Rawashdeh, M.S. Numerical solution for (2+1)-dimensional time-fractional coupled Burger equations using fractional natural decomposition method. Math. Methods Appl. Sci. 2019, 42, 3409–3427. [Google Scholar] [CrossRef]
  29. Belgacem, F.B.M.; Silambarasan, R. Advances in the natural transform. AIP Conf. Proc. 2012, 1493, 106–110. [Google Scholar]
  30. Khan, Z.H.; Khan, W.A. N-transform properties and applications. NUST J. Eng. Sci. 2008, 1, 127–133. [Google Scholar]
  31. Hilfer, R. Applications of Fractional Calculus in Physics; World Sci. Publishing: River Edge, NJ, USA, 2000. [Google Scholar]
  32. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998; Volume 198. [Google Scholar]
  33. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations, 1st ed.; Wiley: Hoboken, NJ, USA, 1993. [Google Scholar]
Figure 1. Exact and NTDM solution for γ = 1 of Problem 4.1.
Figure 1. Exact and NTDM solution for γ = 1 of Problem 4.1.
Entropy 23 01086 g001
Figure 2. For different values of γ of Problem 4.1.
Figure 2. For different values of γ of Problem 4.1.
Entropy 23 01086 g002
Figure 3. The exact and NTDM solution for γ = 1 of Problem 4.2.
Figure 3. The exact and NTDM solution for γ = 1 of Problem 4.2.
Entropy 23 01086 g003
Figure 4. For different value of γ of Problem 4.2.
Figure 4. For different value of γ of Problem 4.2.
Entropy 23 01086 g004
Figure 5. The exact and NTDM solution for γ = 1 of Problem 4.3.
Figure 5. The exact and NTDM solution for γ = 1 of Problem 4.3.
Entropy 23 01086 g005
Figure 6. For different values of γ of Problem 4.3.
Figure 6. For different values of γ of Problem 4.3.
Entropy 23 01086 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Agarwal, R.P.; Mofarreh, F.; Shah, R.; Luangboon, W.; Nonlaopon, K. An Analytical Technique, Based on Natural Transform to Solve Fractional-Order Parabolic Equations. Entropy 2021, 23, 1086. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081086

AMA Style

Agarwal RP, Mofarreh F, Shah R, Luangboon W, Nonlaopon K. An Analytical Technique, Based on Natural Transform to Solve Fractional-Order Parabolic Equations. Entropy. 2021; 23(8):1086. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081086

Chicago/Turabian Style

Agarwal, Ravi P., Fatemah Mofarreh, Rasool Shah, Waewta Luangboon, and Kamsing Nonlaopon. 2021. "An Analytical Technique, Based on Natural Transform to Solve Fractional-Order Parabolic Equations" Entropy 23, no. 8: 1086. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081086

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