Applications on Double ARA–Sumudu Transform in Solving Fractional Partial Differential Equations

: In this article, we apply the double ARA–Sumudu transformation (DARA-ST) to the nonlocal fractional Caputo derivative operator. We achieve interesting results and implement them to solve certain classes of fractional partial differential equations (FPDEs). Several physical applications are discussed and analyzed, such as telegraph, Klein–Gordon and Fokker–Planck equations. The new technique with DARA-ST is efﬁcient and accurate in examining exact solutions of FPDEs. In order to show the applicability of the presented method, some numerical examples and ﬁgures are illustrated. A symmetry analysis is used to verify the results.


Introduction
The fractional calculus generalizes the operations of differentiation and integration to noninteger orders. The fractional calculus has become an important tool for the study of some physical phenomena, engineering and science, such as electromagnetics, viscoelasticity, fluid mechanics, electrochemistry, biological population models, optics and signal processing. Furthermore, fractional calculus processes have become one of the most useful approaches in a variety of applied sciences to deal with certain properties of (long) memory effects. There are many definitions of fractional derivatives, such as Riemann-Liouville, Caputo, Caputo-Fabrizio, Atangana-Baleanu, conformable and the generalized fractional derivative [1][2][3][4][5][6].
The method of double integral transforms is a hot topic in recent research, and it basically depends on applying a single transformation twice on functions of two variables or applying two different transformations on the same function. This new approach is a powerful tool for solving PDEs. Although double integral transformations, their properties and theorems are recent studies, they have attracted the interest of many mathematicians.
where f (x, t) is a continuous function of two variables x > 0 and t > 0.
In this article, we implement DARA-ST to solve families of FPDEs of the form x, t ≥ 0, n − 1 < α ≤ n , m − 1 < β ≤ m and m, n ∈ N, (1) with the initial conditions (ICs) ∂ j f (x, 0) ∂t j = g j (x), j = 0, 1, · · · , m − 1 (2) and the boundary conditions (BCs) where A, B and C are real constants, D α x and D β t are the fractional Caputo's derivatives with respect to x and t, respectively, L is a linear operator and z(x, t) is the source function.
The main motivation of the present study is to expand the applications of DARA-ST by using it to solve FPDEs. We show the efficiency of the proposed method by applying the DARA-ST to several interesting applications to obtain the exact solutions and analyze the results. The novelty of this work arises from the establishment of a new simple formula for solving PDEs of fractional orders. The simplicity and applicability of this new formula is illustrated by handling some applications, where we use the new approach to solve some important FPDEs.
This article is organized as follows: in the next two sections, we present some basic definitions and theorems related to our work. A new algorithm for solving families of FPDEs using DARA-ST is presented in Section 4. Several examples are given in Section 5 to demonstrate the proposed technique. We illustrate the numerical evaluations of the results in Section 6. Finally, our results are discussed in Section 7.

Sumudu and ARA Transformations
In this section, we introduce the definition of Sumudu and ARA transforms with some properties. Definition 1 [33]. Sumudu transform of the function f (x) is defined as Definition 2 [34]. ARA transform of order n of the function f (x) is given by and the ARA transform of the function f (x) of order one is defined as For simplicity, let us denote Theorem 1 [33]. (The sufficient condition for the existence of Sumudu transform).
If the function f (x) is a piecewise continuous in every finite interval 0 ≤ x ≤ α and satisfies then Sumudu transform exists for all 1 u > α.
Theorem 2 [34]. (The sufficient condition for the existence of ARA transform).
If the function f (x) is a piecewise continuous in every finite interval 0 ≤ x ≤ β and satisfies then ARA transform exists for all s > β.
Definition 5. The Mittag-Leffler function is defined by The single ARA transform of The single Sumudu transform of t β−1 E α,β (λt α ) takes the value    Table 2 illustrates the values of DARA-ST for some basic functions. Table 2. DARA-ST for some functions [35].
The proof is completed.

Algorithm of DARA-ST Method
In this section, we present the technique of using DARA-ST to solve families of FPDEs. In order to achieve our goal, we have to calculate DARA-ST for the nonlocal Caputo fractional derivative in the following lemma.
from the definition of the convolution, we have Using the convolution property of ARA transform in Table 1, we obtain Applying the derivative property of ARA transform in Table 1, we obtain After simple computations, we obtain from the definition of the convolutions, we have Using the convolution property of Sumudu transform in Table 1, we obtain Applying the derivative property of Sumudu transform in Table 1, we obtain After simple computations, we obtain .

Solving FPDEs by DARA-ST
In this section, we apply DARA-ST to obtain solutions of some FPDEs. We consider the initial boundary value problems (1)- (3). To obtain the solution by the new approach, we apply DARA-ST on both sides of Equation (1), to obtain Furthermore, we apply the single ARA transform to the ICs (3), and the single Sumudu transform to the BCs (2), to obtain Simplifying Equation (4), and substituting the values in Equations (5) and (6), we have Running the inverse DARA-ST, G −1 x S −1 t on both sides of Equation (7), we obtain which is the solution of the target problem.

Illustrative Examples
In this section, we introduce some famous PDEs in mathematical physics such as Reaction-diffusion, advection-diffusion, telegraph equation, wave equation, Klein-Gordon and Fokker-Planck, we apply the new double transform on these equations and use it to obtain the solution of these problems and we implement the obtained formula in Equation (8) to solve FPDEs, to handle these problems using the new approach. The main goal here is to illustrate the applicability and ease of use of the new double transform.

Fractional Advection-Diffusion Equation
Consider the fractional advection-diffusion equation with the IC and the BCs Applying the single ARA transform on ( ) in Equation (17), we obtain Applying the single Sumudu transform on ℎ ( ) and ℎ ( ) in Equation (18)

1−2u β
in Equation (19), we obtain the solution of (20) as follows with the IC and the BCs  (19), we obtain the solution of (20) as follows

Fractional Telegraph Equation
Consider the fractional telegraph equation Applying the single ARA transform on ( ) and ( )in Equation (24), we obtain

Fractional Telegraph Equation
Consider the fractional telegraph equation with the ICs f (x, 0) = g 0 (x), f t (x, 0) = g 1 (x) (24) and the BCs Applying the single ARA transform on g 0 (x) and g 1 (x) in Equation (24), we obtain Applying the single Sumudu transform on h 0 (t) and h 1 (t) in Equation (25), we obtain , Substituting L[ f (x, t)] = C 0 f (x, t) + C 1 f t (x, t) = 0, z(x, t) = 0, n = m = 2, and the functions G 0 (s), G 1 (s), H 0 (u), H 1 (u) in the general formula in Equation (8) and after simple computations, we obtain Example 3. Consider the telegraph equation with the ICs and the BCs 1+u and H 1 (u) = 1 1+u in Equation (26), we obtain the solution of (27) as follows  (26), we obtain the solution of (27) as follows

Fractional Wave Equation
Consider the fractional wave equation with the ICs and the BCs

Fractional Fokker-Planck Equation
Consider the fractional Fokker-Planck equation with the IC and the BCs Applying the single ARA transform on ( ) in Equation (43), we obtain

Fractional Fokker-Planck Equation
Consider the fractional Fokker-Planck equation with the IC and the BCs Applying the single ARA transform on g 0 (x) in Equation (43), we obtain Applying the single Sumudu transform on h 0 (t) and h 1 (t) in Equation (44), we obtain , , z(x, t) = 0, n = 2, m =1, and the functions G 0 (s), H 0 (u), H 1 (u) in the general formula in Equation (8) and after simple computations, we obtain Example 5. Consider the fractional Fokker-Planck equation with the IC and the BCs Solution. Putting α = 2, G 0 (s) = 1 s , H 0 (u) = u β , H 1 (u) = 1 in Equation (45), we obtain the solution of (46) as follows with the IC

Numerical Simulations
In this section, we illustrate the numerical evaluation of the solutions obtained by solving the given FPDEs. We also discuss the numerical behavior of the results when solving FPDEs, then we compare it with the solution of the equation of integer order.
The solutions of Examples 1, 2 and 5 are simply computed when β = 1. We examine the numerical solutions of different values of β = 0.95, 0.85 and 0.75. As a result, we notice that, with choosing different values of β, the obtained fractional solutions are in coordination with the closed form of the solution when β = 1, as illustrated in Figure 6.
Moreover, it clearly implies that as β approaches 1, the obtained solutions of the FPDEs approach the exact solutions obtained in the integer case.

Numerical Simulations
In this section, we illustrate the numerical evaluation of the solutions obtained by solving the given FPDEs. We also discuss the numerical behavior of the results when solving FPDEs, then we compare it with the solution of the equation of integer order.
The solutions of Examples 1, 2 and 5 are simply computed when = 1. We examine the numerical solutions of different values of = 0.95, 0.85 and 0.75. As a result, we notice that, with choosing different values of , the obtained fractional solutions are in coordination with the closed form of the solution when = 1, as illustrated in Figure 6. Moreover, it clearly implies that as approaches 1, the obtained solutions of the FPDEs approach the exact solutions obtained in the integer case.  Moreover, we mention that as and approach the close integer orders, the obtained solutions of the FPDEs approach the exact solutions in the integer case.

Conclusions
In this research, DARA-ST is applied to the Caputo fractional derivative to obtain a new interesting formula, that is implemented to solve families of FPDEs. We have presented a new method to obtain exact solutions of these equations. We show the reliability and efficiency of the proposed method by presenting some interesting physical applications. In the future, we will pair DARA-ST with some iteration methods to solve nonlinear FPDEs, such as nonlinear telegraph equation, nonlinear wave equation, nonlinear Klein-Gordon and nonlinear Fokker-Planck. In addition, researchers can use new definitions of FC such as the generalized fractional derivative and others to search and obtain new results on transformations.   Moreover, we mention that as α and β approach the close integer orders, the obtained solutions of the FPDEs approach the exact solutions in the integer case.

Conclusions
In this research, DARA-ST is applied to the Caputo fractional derivative to obtain a new interesting formula, that is implemented to solve families of FPDEs. We have presented a new method to obtain exact solutions of these equations. We show the reliability and efficiency of the proposed method by presenting some interesting physical applications. In the future, we will pair DARA-ST with some iteration methods to solve nonlinear FPDEs, such as nonlinear telegraph equation, nonlinear wave equation, nonlinear Klein-Gordon and nonlinear Fokker-Planck. In addition, researchers can use new definitions of FC such as the generalized fractional derivative and others to search and obtain new results on transformations.