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

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.


Introduction
The present research work is dedicated to studying the analytical solution of fractionalorder 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.
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].

Preliminaries
Definition 1. The following transformation off (η) is called natural transformation and is expressed as [29,30] N where the transformation parameters are s and u.

Definition 2.
The following transformation off (η) is called inverse natural transformation and is expressed as 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 NTf (η) is thef (η) and is defined as 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 a where h 1 * l 1 represents the convolution of h 1 and l 1 .
where Γ is a function defined by Definition 5. The Caputo fractional derivative operator with order γ is defined as [33]

Idea of FNTM
The general fractional-order PDEs are given as 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 Applying the natural transformation to Equation (1), we get Using the natural transform's differentiation property, we get The following infinite series represents the NTDM solution υ(ξ, η) and Adomian polynomials as We get replacement Equations (5) and (6) in Equation (4).
Applying the natural transformation's linearity, We can generally write Equations (9) and (10) implementing the inverse natural transformation
The Equation (25) correction function is provided by The first term being Then we get for j = 0 The following terms are The series forms a solution to Problems (4.2), for example, Then γ = 1; the integer NDM results as 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.

Problem
Consider fractional-order three-dimensional parabolic equation: with initial condition with boundary conditions Concerning the natural transformation of (32), we get using the inverse natural transform.
The Equation (35) correction function is provided by The first term being Then we get for j = 0 The following terms are The series forms a solution to Problems (4.3), for example, Then γ = 1; the integer NDM results as 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.

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.