Applications and Utilization of Laplace Transform in Science and Engineering ProblemsView this Special Issue
A Hybrid Interpolation Method for Fractional PDEs and Its Applications to Fractional Diffusion and Buckmaster Equations
This study presents a novel numerical method to solve PDEs with the fractional Caputo operator. In this method, we apply the Newton interpolation numerical scheme in Laplace space, and then, the solution is returned to real space through the inverse Laplace transform. The Newton polynomial provides good results as compared to the Lagrangian polynomial, which is used to construct the Adams–Bashforth method. This procedure is used to solve fractional Buckmaster and diffusion equations. Finally, a few numerical simulations are presented, ensuring that this strategy is highly stable and quickly converges to an exact solution.
Recently, partial differential equations (PDEs) have been successfully used to handle many real problems in engineering and mathematical physics such as viscoelastic materials, seismic analysis, and viscous damping [1–3]. Over the past few decades, fractional calculus has been used for better analysis of mathematical models [4, 5]. Fractional operators have been implemented to study nonlocal behaviour of the physical model. Khan et al. investigated hidden attractors in four-dimensional dynamical systems with power law kernel . Khan et al. analyzed the imperfect testing infection disease model using the fractional operator . Some other applications of fractional calculus are listed in [8, 9]. Also, fractional PDEs have attracted the scientists due to their wide range of applications in various areas of science. For instance, Saifullah et al. investigated the Klein–Gordon equation with power law kernel . Gulalai et al. analyzed nonlinear modified KdV equation under Atangana–Baleanu–Caputo derivative . Fractional PDEs have memory terms; they are not the same as classical PDEs, and evaluating FPDEs numerically is much difficult than solving PDEs. So, it is essential to establish effective numerical techniques to deal with PFDEs. There are several analytical and numerical methods for solving fractional PDEs. Akram et al. studied applied extended cubic B-spline collocation scheme for time fractional subdiffusion equation, cable equation, and Klein–Gordan equation [12–14]. Ahmad et al. used Yang transfrom coupled with the homotopy perturbation method (HPM) to study fractional KdV and Fornberg–Whitham equations . The Legendre wavelets method has been utilized to investigate parabolic fractional PDEs . A biorthogonal Hermite cubic spline Galerkin method has been used to solve fractional Riccati equation . Iqbal et al. applied the Elzaki transformation and the new iteration technique to study fractional-order Kaup–Kupershmidt equation . Shah et al. utilized the Haar wavelet collocation method to analyze the fractional-order COVID-19 model . Many other approaches have been used for computing solution of fractional PDEs [20, 21].
The Adams–Bashforth method has been recognized as a powerful numerical scheme for solving nonlinear partial differential equations with integer and noninteger derivatives. Gnitchogna and Atangana proposed a new numerical scheme which can be seen as an extension of the Adams–Bashforth one for partial differential equations using Lagrangian interpolation . There is a need to expand Newton interpolation numerical schemes to PFDEs because of their accuracy and efficiency in integer-order and fractional derivatives . To handle PFDE with Caputo derivative, we develop a numerical scheme that combines the Newton interpolation numerical scheme and Laplace transform. By removing one variable and converting a PFDE to an FDE using the Laplace transform, the newly generated FDE can be easily solved in Laplace space, and the result may then be returned to real space using the inverse transform. Finally, some examples are provided to assess the efficacy of the developed method.
The study is organized as follows: Section 2 provides the basic concepts of fractional calculus. Section 3 deals with numerical approximation of fractional PDEs using Newton polynomial. Sections 4 and 5 are the applications of the suggested scheme. Section 6 gives the numerical simulations of the proposed results. Last, the conclusion is provided in Section 7.
2. Basic Concepts of Caputo Operator
Here, we recall important definitions of Laplace transform and fractional derivatives .
Definition 1. The Caputo-type fractional derivative is given aswhere .
Definition 2. The Laplace transform (LT) of the function is
Definition 3. The LT of iswhere .
3. Numerical Approximations
Consider a general fractional nonlinear PDE aswhere is a nonlinear operator and is a linear operator. Applying LT to (4), we get
Let and . Then, (6) becomes
Using the Caputo fractional integral operator, we get
For , one gets
Next, we approximate the above equation aswhere . So, the kernel might be approximated with in via Newton’s polynomial. So,
Substituting equation (11) in (10), we get
Using the linearity property of integral, we have
Now, we evaluate the integrals in (13) as
By substituting these integrals values into (13), we get
Rearranging the terms and simplifying, we have
The inverse Laplace transform is used to get the solution as
Furthermore, we can discretize equation (16) as follows:where . Equation (17) gives the numerical approximation of the fractional PDEs using Newton polynomial.
In the following sections, this novel method will be applied to solve a fraction diffusion problem and a fractional Buckmaster equation.
4. Applications to Fractional Diffusion Equation
Diffusion equation is parabolic PDE having many applications in real world phenomena. In physics, it represents the macroscopic behaviour of multiple microparticles in Brownian motion, which is caused by the particles’ random motions and collisions. Here, we evaluate the accuracy and effectiveness of the newly proposed numerical method by applying it to a fractional diffusion model. Consider
Using Laplace transform, we get
Let and .
Equation (19) can be rewritten as
We already established that the solution to equation (22) in the Laplace space is given in equation (15). So,
Let us define the expression to simplify the formulas aswhere , .
The inverse Laplace transform is used to get the solution as
Discretizing in the space variable, we get
Letting and , equation (24) becomes
Thus, equation (25) is the required numerical results for fractional diffusion equation using Newton polynomial.
5. Application to the Fractional Buckmaster Equation
The Buckmaster equation is a nonlinear PDE that was developed in 1977 by John D. Buckmaster. The surface of a thin sheet of viscous liquid is represented by this equation. In this section, we provide another application of the proposed scheme to show the its validity. We consider the fractional Buckmaster equation as
Taking the Laplace transform of (29), we get
Let and . Now, consider
We already established that the solution to (28) in Laplace space is represented by (15). So,
Taking the inverse Laplace transform, we getwhere , .
Discretizing in the space variable and substituting , we get
Thus, equation (31) gives the numerical solution of the proposed fractional Buckmaster equation.
6. Numerical Simulations
In this section, we demonstrate numerical solutions to diffusion and Buckmaster equations for various fractional orders. Figures 1–6 show the simulation of the numerical solution of the fractional diffusion equation for fractional order . At , the solution of the proposed problem converges to the exact solution. The numerical simulations and their associated contour plot are presented to show the efficiency and stability of the numerical scheme, which can be observed from Figures 1–6. The numerical simulation of the obtained numerical solution of the Buckmaster problem is shown in Figures 7–12 for fractional order . The last two figures reveal that the solution directly converges to the exact solution of the proposed problem. Thus, the efficiency and applicability of the suggested hybrid interpolation method is verified in these two applications.
In this study, we have devised a hybrid Newton polynomial interpolation approach to solve fractional PDEs. We have used the Caputo fractional operator. This technique combines Newton interpolation with the Laplace transform. We have derived a general numerical scheme for solution of fractional PDE using Laplace transform and Newton interpolation. The approach was verified by solving two PDEs, including fractional diffusion and Buckmaster equations. All the results have been validated through numerical simulations. In future, this method can be used to solve more difficult nonlinear fractional PDEs.
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
S. Khaliq, A. Ullah, S. Ahmad, A. Akgül, A. Yusuf, and T. A. Sulaiman, “Some novel analytical solutions of a new extented (2 + 1)-dimensional Boussinesq equation using a novel method,” J. Ocean Eng. Sci., 2022, in Press.View at: Google Scholar
S. Y. Lou and X. Hao, “Linear superposition in the general heavenly equation,” Physics Letters A, vol. 443, Article ID 128203, 2022.View at: Publisher Site | Google Scholar
R. Wu, X. Shen, and J. Zhao, “Conforming finite element methods for two-dimensional linearly elastic shallow shell and clamped plate models,” Applied Mathematics and Computation, vol. 430, Article ID 127259, 2022.View at: Publisher Site | Google Scholar
S. Ahmad, A. Ullah, M. Partohaghighi, S. Saifullah, A. Akgül, and F. Jarad, “Oscillatory and complex behaviour of Caputo-Fabrizio fractional order HIV-1 infection model,” AIMS Mathematics, vol. 7, no. 3, pp. 4778–4792, 2021.View at: Publisher Site | Google Scholar
A. Khan, H. M. Alshehri, T. Abdeljawad, Q. M. Al-Mdallal, and H. Khan, “Stability analysis of fractional nabla difference COVID-19 model,” Results in Physics, vol. 22, Article ID 103888, 2021.View at: Publisher Site | Google Scholar
Z. A. Khan, J. Khan, S. Saifullah, and A. Ali, “Dynamics of hidden attractors in four-dimensional dynamical systems with power law,” Journal of Function Spaces, vol. 2022, Article ID 3675076, 22 pages, 2022.View at: Publisher Site | Google Scholar
Z. A. Khan, A. Khan, T. Abdeljawad, and H. Khan, Computational Analysis of Fractional Order Imperfect Testing Infection Disease Model, Fractals, NewYork, NY, USA, 2022.View at: Publisher Site
A. Devi, A. Kumar, T. Abdeljawad, and A. Khan, “Stability analysis of solutions and existence theory of fractional Lagevin equation,” Alexandria Engineering Journal, vol. 60, no. 4, pp. 3641–3647, 2021.View at: Publisher Site | Google Scholar
P. Kumam, A. Seadawy, G. Zaman, and Z. Shah, “Study of mathematical model of Hepatitis B under Caputo-Fabrizo derivative,” AIMS Mathematics, vol. 6, no. 1, pp. 195–209, 2021.View at: Publisher Site | Google Scholar
S. Saifullah, A. Ali, M. Irfan, and K. Shah, “Time-fractional Klein–Gordon equation with solitary/shock waves solutions,” Mathematical Problems in Engineering, vol. 2021, Article ID 6858592, 15 pages, 2021.View at: Publisher Site | Google Scholar
S. Gulalai, F. A. Ahmad, F. A. Rihan, A. U. Rihan, Q. M. Al-Mdallal, and A. Ali, “Nonlinear analysis of a nonlinear modified KdV equation under Atangana Baleanu Caputo derivative,” AIMS Mathematics, vol. 7, no. 5, pp. 7847–7865, 2022.View at: Publisher Site | Google Scholar
T. Akram, M. Abbas, and A. I. Ismail, “An extended cubic B-spline collocation scheme for time fractional sub-diffusion equation,” in Innovation, Penang, Malaysia, December 2018.View at: Publisher Site | Google Scholar
T. Akram, M. Abbas, and A. I. Ismail, “Numerical solution of fractional cable equation via extended cubic B-spline,” in Proceedings of tand Analytics Conference & Exhibition (IACE 2019), Kedah, Malaysia, March 2019.View at: Publisher Site | Google Scholar
T. Akram, M. Abbas, M. B. Riaz, A. I. Ismail, and N. M. Ali, “Development and analysis of new approximation of extended cubic B-spline to the non-linear time fractional Klein-Gordon equation,” Fractals, vol. 28, no. 08, Article ID 2040039, 2020.View at: Publisher Site | Google Scholar
S. Ahmad, A. Ullah, A. Akgül, and F. Jarad, “A hybrid analytical technique for solving nonlinear fractional order PDEs of power law kernel: application to KdV and Fornberg-Witham equations,” AIMS Mathematics, vol. 7, no. 5, pp. 9389–9404, 2022.View at: Publisher Site | Google Scholar
S. S. Ray and A. K. Gupta, “Numerical solution of fractional partial differential equation of parabolic type with Dirichlet boundary conditions using two-dimensional Legendre wavelets method,” Journal of Computational and Nonlinear Dynamics, vol. 11, no. 1, Article ID 011012, 2016.View at: Publisher Site | Google Scholar
H. Bin Jebreen and I. Dassios, “A biorthogonal hermite cubic spline Galerkin method for solving fractional riccati equation,” Mathematics, vol. 10, no. 9, 1461 pages, 2022.View at: Publisher Site | Google Scholar
N. Iqbal, H. Yasmin, A. Rezaiguia, J. Kafle, A. O. Almatroud, and T. S. Hassan, “Analysis of the fractional-order kaup–kupershmidt equation via novel transforms,” Journal of Mathematics, vol. 2021, Article ID 2567927, 13 pages, 2021.View at: Publisher Site | Google Scholar
K. Shah, Z. A. Khan, A. Ali, R. Amin, H. Khan, and A. Khan, “Haar wavelet collocation approach for the solution of fractional order COVID-19 model using Caputo derivative,” Alexandria Engineering Journal, vol. 59, no. 5, pp. 3221–3231, 2020.View at: Publisher Site | Google Scholar
N. A. Obeidat and D. E. Bentil, “Novel technique to investigate the convergence analysis of the tempered fractional natural transform method applied to diffusion equations,” Journal of Ocean Engineering and Science, 2022.View at: Publisher Site | Google Scholar
R. P. Agarwal, F. Mofarreh, R. Shah, W. Luangboon, and K. Nonlaopon, “An analytical technique, based on natural transform to solve fractional-order parabolic equations,” Entropy, vol. 23, no. 8, 1086 pages, 2021.View at: Publisher Site | Google Scholar
R. Gnitchogna and A. Atangana, “New two step Laplace Adam-Bashforth method for integer a noninteger order partial differential equations,” Numerical Methods for Partial Differential Equations, vol. 34, no. 5, pp. 1739–1758, 2017.View at: Publisher Site | Google Scholar
A. Atangana and S. I. Araz, New Numerical Scheme with Newton Polynomial Theory, Methods, and Applications, Elsevier, Amsterdam, Netherlands, 2021.
D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional calculus: models and numerical methods, World Scientific, Singapore, 2012.