A Reliable Technique for Solving Fractional Partial Differential Equation

: The development of numeric-analytic solutions and the construction of fractional-order mathematical models for practical issues are of the greatest importance in a variety of applied mathematics, physics, and engineering problems. The Laplace residual-power-series method (LRPSM), a new and dependable technique for resolving fractional partial differential equations, is introduced in this study. The residual-power-series method (RPSM), a well-known technique, and the Laplace transform (LT) are elegantly combined in the suggested technique. This innovative approach computes the fractional derivative in the Caputo sense. The proposed method for handling fractional partial differential equations is provided in detail, along with its implementation. The novel approach yields a series solution to fractional partial differential equations. To validate the simplicity, effectiveness, and viability of the suggested technique, the provided model is tested and simulated. A numerical and graphical description of the effects of the fractional order γ on approximating the solutions is provided. Comparative results show that the suggested method approximates more precisely than current methods such as the natural homotopy perturbation method. The study showed that the aforementioned method is straightforward, trustworthy, and suitable for analysing non-linear engineering and physical issues.


Introduction
It has been noted that fractional-order α derivatives, often in the range between 0 and 1, are a helpful tool for describing a variety of events [1]. To obtain the desired order of a fractional derivative, the Caputo fractional derivative computes an ordinary derivative first, then a fractional integral. In reverse order, the Riemann-Liouville fractional derivative is calculated. As a result, although the Riemann-Liouville fractional derivative permits initial conditions in terms of fractional integrals and their derivatives, and the Caputo fractional derivative only enables the involvement of conventional initial and boundary conditions. These two operators coincide under the homogeneous initial condition assumption. Anyone who has studied fundamental calculus is familiar with the differentiation operator D = d dς . Furthermore, provided that n is a positive integer, the nth derivative of u, denoted by D γ u(ς) = d n u(ς) dς n , is well defined for suitable functions u. L'Hôpital asked Leibnitz in 1695 what significance could be assigned to D n u(ς) if n were a fraction. However, it was not until 1884 that the theory of generalized operators reached a stage in its evolution that was appropriate for the modern mathematician to use as a starting point. By that time, the theory had been expanded to include operators D v , where m could be real or complex, rational or irrational, and positive or negative [2]. Although derivatives of arbitrary order were discussed by Leibniz, Euler, Laplace, Lacroix, and Fourier, Niels Henrik Abel used fractional operations for the first time in 1823. In order to solve the tautochrone paradox, Abel used fractional calculus [2]. Liouville is to be credited with making what is likely the first sincere attempt to define a fractional derivative logically. He wrote nine articles on the topic between 1832 and 1837, with the latest one in the field appearing in 1855. In recent years, it has been discovered that the fractional calculus is very effective in describing a wide variety of physical phenomena, including damping laws and diffusion processes [2][3][4][5]. Kilbas and Trujillo [1], Caputo [6], Debanth [7], Jafari and Seifi [8], Kemple and Beyer [9], Oldham and Spanier [10], Momani and Shawagfeh [11], and others provide some fundamental works on various elements of the fractional calculus [12][13][14][15][16].
Over the past forty years, fundamental research and advancements on fractional derivatives and differential equations have been made. Traditional differential equations with non-local and genetic significance in material characteristics are generalized as fractional-order differential equations. Fractional partial differential equations are increasingly used in the creation of non-linear models and the analysis of dynamical systems. The theory of fractional-order calculus has been connected to real-world projects and used to examine and investigate a variety of phenomena, such as chaos theory [17], financial models [18], a noisy environment [19], optics [20], and others [21][22][23][24]. The characteristics of non-linear issues that occur in nature are largely described by the solutions of fractional differential equations. Since it is challenging to find an accurate solution for fractional differential equations representing non-linear phenomena, many analytical and numerical techniques are employed [25][26][27][28].
In recent years, scholars and researchers have paid close attention to both the numerical and analytical solutions of PDE systems. For resolving fractional FPDEs such as this, numerous numerical and analytical algorithms have been used, including the first integral method [29], the Elzaki transform decomposition method [30,31], the double Laplace transform method [32], the homotopy perturbation transform method [33,34], the conformable fractional Laplace transform method [35], the Yang transform decomposition method [36,37], the generalized two-dimensional differential transform method [38], the Fourier transform [39], He's variational iteration method [40], the fractional complex differential transformation method [41], and the fractional variational iteration method [42].
The power-series method (PSM), which results in a closed-form solution of known functions, is well proven as an efficient method for solving linear ordinary-partial differential equations. In the case of nonlinear problems, it is impossible to obtain a closed-form solution, and finding out the series coefficients is a highly challenging task. A modified version of the PSM that treats the coefficients as transformed functions that follow a set of rules and are determined in recurrence relations is introduced to address the aforementioned limitations of the standard PSM. The differential transform method is the name of this improvement (DTM). Different kinds of integro-differential equations and linear-nonlinear equations have both been solved using it. Another advancement is the establishment of the residual-power-series method (RPSM) through the differentiation of the nth ordered coefficient of the PSM's nth partial sum of the PSM (n-1)-times.
It was necessary to increase the use of the power-series method to deal with fractional difficulties during the modification of the ordinary derivative to a fractional derivative because it is more general. Many significant models that arise in various branches of science and engineering are constructed and solved analytically using the fractional DTM and fractional RPS methodologies. By including the Laplace transform (LT) into the RPSM's technique, we hope to improve its accuracy in this work. This RPSM promotion is known as the Laplace residual-power-series method (LRPSM). Solving the FPDEs introduces the building of this innovative approach. Accuracy to the necessary level has been attained. The suggested technique has a very easy and uncomplicated process. The findings indicate that, in comparison to other analytical procedures, the current method has the appropriate accuracy.
The framework of the study is detailed as follows. First, we use key FC theory ideas and findings in Section 2. Additionally, several original findings that provide the basis for the innovative technique in Section 2 are provided. The solutions to time-fractional PDEs are then determined in Section 3 using the LRPSM. Some of the problems in Section 4 are solved using LRPSM. A brief conclusion ends Section 6.

Preliminaries
Here, we provide some definitions in terms of Caputo and Riemann-Liouville, along with the Laplace transform theorem . Definition 1. The fractional derivative in terms of Caputo is stated as [6,43] where m ∈ N and J γ η represents a fractional integral in terms of Riemann-Liouville (RL) as (2) Definition 2. The LT is stated as [43] u with inverse LT as Proof. For proof, see [44].

Remark 1.
On taking the inverse LT of Equation (5) as provided in [43]: This corresponds to the fractional Taylor's formula described in [45]. The convergence of the FPS in Theorem (1) is explained and proven by the following theorem.

General Methodology of LRPSM
with initial source By employing LT to (7), As from the fact that C D a and by utilizing the initial condition (8), we have We may express the transformed function U(ς, µ) in the following manner: The kth-truncated series of (11) can be expressed as As provided in [46], from the definition of the Laplace residual function We provide several features that emerge in the common residual power series approach [46]: . . We will now solve the system below recursively in order to define the coefficient functions The next step is to take the inverse LT of U k (ς, µ) to obtain the kth approximation u k (ς, η).

Numerical Examples
Here, we solve three problems to show the accuracy of the proposed method.

Problem
Assume the fractional partial differential equation of the following form: with the initial condition By employing LT to Equation (1) and by using Equation (9), we obtain The kth-truncated series can be expressed as so the kth-LRFs are The kth-truncated series Equation (17) will be substituted into the kth-truncated residual function Equation (18) to yield f k (ς, , µ). The resulting equation, µ kγ+1 , will then be multiplied, and the relation lim µ→∞ (µ kγ+1 L t Res u,k (ς, , µ)) = 0, k = 1, 2, 3, · · · . Several terms are as and so on.
By employing inverse LT, we obtain On putting γ = 1, we have
By employing inverse LT, we obtain On putting γ = 1, we have

Problem
Assume the fractional partial differential equation of the following form: with the initial condition By employing LT in Equation (31) and by using Equation (32), we obtain The kth-truncated series can be expressed as so the kth-LRFs are The kth-truncated series Equation (34) will be substituted into the kth-truncated residual function Equation (35) to yield f k (ς, , µ). The resulting equation, µ kγ+1 , will then be multiplied, and the relation lim µ→∞ (µ kγ+1 L t Res u,k (ς, , µ)) = 0, k = 1, 2, 3, · · · . Several terms are as and so on.
By employing inverse LT, we obtain On putting γ = 1, we have

Results and Discussion
The numerical analysis between exact and approximative solutions, as shown in Tables 1-3, has been investigated in detail and with more precision in this study. The correctness and simplicity of the suggested method are demonstrated by computing the numerical values at different fractional orders. Tables 1-3 display the numerical comparison of the accurate and approximative solutions, demonstrating that the series solution soon converges to a small value. As a result, adding more terms for an approximate solution increases the accuracy of the analytical result. Figure 1 shows how the accurate and suggested approaches behave, as well as the characteristics of the approximative solution. For a better understanding of the problem's characteristics, we also provide the suggested method solution at various fractional orders in Figures 2 and 3. Figure 4 calculates the solution to problem 2 using the suggested and actual method. Figure 5 displays the graphical representations for γ = 0.8, 0.6. The behaviour of problem 2 in 2D and 3D for different fractional orders is shown in Figure 6. Similarly, Figure 7 presents the actual and suggested methods' solutions at γ = 1, whereas Figures 8 and 9 present the proposed approach solution at different fractional orders. Based on the tables and graphs, we came to the conclusion that the proposed technique solution was in good agreement with the precise solution.

Conclusions
In this study, the Laplace residual-power-series method (LRPSM), a powerful new technique for solving fractional partial differential equations, is developed by successfully combining the residual-power-series method (RPSM) and the Laplace transform. The new technique provides a series solution with elegant computational terms that quickly converges to an exact or approximate solution. The fractional derivative is handled in the Caputo sense in this novel analytical technique. With the help of the new analytical technique, fractional partial differential equations are successfully solved precisely. The obtained results via our technique are compatible with the results obtained by the natural homotopy perturbation method. In order to understand the behavior of the provided problems, solutions at different fractional orders are taken and are shown with the help of graphs and tables, which confirm that we get closer to the exact solution as the order of γ goes from fractional-order towards the integercorder. It is clearly proven that the new reliable approach is both straightforward and highly accurate. The Laplace residualpower-series method is a powerful and reliable technique for handling fractional partial differential equations that are both linear and nonlinear.