Highly Efficient Energy Stable Schemes for Multi-dimensional Space Fractional Reaction-Diffusion Models

Abstract

Several important physical phenomena in engineering and scientific fields are modeled by non-local fractional models. The development of numerical schemes becomes crucial due to the non-availability of the exact solutions to such models. However, the numerical approximation of these models is challenging and imposes several computational constraints. In this paper, we have devised two highly efficient, energy-stable numerical schemes to solve space fractional reaction-diffusion models. Spatial discretization is performed using a fourth-order matrix transform technique having the advantage of straightforward extension to two and higher spatial dimensions. The time-stepping schemes are developed using an exponential time differencing approach based on a third-order real-pole restricted Padé approximation and a third-order Padé(1,2) approximation. Computationally efficient versions of the schemes are constructed using a splitting technique. Algorithms based on these schemes are constructed, allowing easy coding, and implemented to perform several numerical experiments on problems of practical interest such as the enzyme kinetics equation, Fisher’s equation, and the Allen-Cahn equation. The proposed schemes allow the accurate and efficient simulation of these dynamical models. Solution profiles are plotted to demonstrate the effectiveness of these schemes. Convergence results are computed to validate the accuracy, and central processing unit time is recorded to show the computational efficiency.

Keywords

Reaction-diffusion models / Matrix transfer technique (MTT) / Exponential time differencing (ETD) / Fractional Laplacian / Numerical methods / 65N35 / 65M22 / 65M12 / 35K11 / 35K57

References

[1] Adam, N., Franke, F., Aland, S.: A simple parallel solution method for the Navier-Stokes Cahn-Hilliard equations. Mathematics 8(8), p.1224 (2020). https://doi.org/10.3390/math8081224

[2] Allen, S.M., Cahn, J.W.: Mechanisms of phase transformations within the miscibility gap of Fe-rich Fe-Al alloys. 24(5), 425–437 (1976)

[3] Allen, S.M., Cahn, J.W.: A microscopic theory for antiphase boundary motion and its application to antipahse domain coarsening. Acta Metall. 27, 1085–1095 (1979)

[4] Alt R. A-stable one-step methods with step-size control for stiff systems of ordinary differential equations. J. Comput. Appl. Math., 1978, 4(1): 29-35

[5] Alzahrani SS, Khaliq A, Biala TA, Furati KM. Fourth-order time stepping methods with matrix transfer technique for space-fractional reaction-diffusion equations. Appl. Numer. Math., 2019, 146: 123-144

[6] Baeumer B, Benson DA, Meerschaert MM, Wheatcraft SW. Subordinated advection-dispersion equation for contaminant transport. Water Resour. Res., 2001, 37: 1543-1550

[7] Bhatt, H., Joshi, J., Argyros, I.: Fourier spectral high-order time-stepping method for numerical simulation of the multi-dimensional Allen-Cahn equations. Symmetry 13(245), 1–15 (2021)

[8] Bu, L., Wu, J.H., Mei, L.Q., Wang, Y.: Second-order linear adaptive time-stepping schemes for the fractional Allen-Cahn equation. Computers & Mathematics with Applications 145, 260–274 (2023)

[9] Bueno-Orovio A, Kay D, Burrage K. Fourier spectral methods for fractional-in-space reaction-diffusion equations. BIT Numer. Math., 2014, 54: 937-954

[10] Çelik C, Duman M. Crank-Nicolson method for the fractional diffusion equation with the Riesz fractional derivative. J. Comput. Phys., 2012, 231: 1743-1750

[11] Certaine, J.: The solution of ordinary differential equations with large time constants. In: Mathematical Methods for Digital Computers, pp. 128–132. Wiley, New York (1960)

[12] Cox, S.M., Matthews, P.C.: Exponential time differencing for stiff systems. J. Comput. Phys. 176, 430–455 (2002)

[13] Dahlquist GG. A special stability problem for linear multistep methods. BIT Numer. Math., 1963, 3(1): 27-43

[14] Ding HF. The construction of an optimal fourth-order fractional-compact-type numerical differential formula of the Riesz derivative and its application. Commun. Nonlinear Sci. Numer. Simul., 2023, 23 107272

[15] Ding HF, Li CP. High-order numerical algorithm and error analysis for the two-dimensional nonlinear spatial fractional complex Ginzburg-Landau equation. Commun. Nonlinear Sci. Numer. Simul., 2023, 120 107160

[16] Ding HF, Li CP, Chen YQ. High-order algorithms for Riesz derivative and their applications (I). Abstr. Appl. Anal., 2014, 2014: 1-17

[17] Ding HF, Li CP, Chen YQ. High-order algorithms for Riesz derivative and their applications (II). J. Comput. Phys., 2015, 293: 219-237

[18] Ding HF, Zhang YX. New numerical methods for the Riesz space fractional partial differential equations. Comput. Math. Appl., 2012, 63: 1135-1146

[19] Ding HF, Zhang YX, He WS, Yang XY. A new numerical method for the Riesz space fractional diffusion equation. Advances. Mater. Res., 2011, 213: 393-396

[20] Du Q, Liu C, Wang X. A phase field approach in the numerical study of the elastic bending energy for vesicle membranes. J. Comput. Phys., 2004, 198: 450-468

[21] Ehle, B.L.: High order A-stable methods for the numerical solution of systems of D.E.’s. BIT Comput. Sci. Numer. Math. 8(4), 276–278 (1968)

[22] Feng XB, Prohl A. Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows. Numer. Math., 2003, 94: 33-65

[23] Fisher RA. The wave of advance of advantageous genes. Ann. Eugen., 1937, 7: 355-369

[24] Furati KM, Yousuf M, Khaliq A. Fourth-order methods for space fractional reaction-diffusion equations with non-smooth data. Int. J. Comput. Math., 2018, 95: 1240-1256

[25] Gallopoulos E, Saad Y. On the parallel solution of parabolic equations. National Aeronautics and Space Administration, 1989

[26] Guo Z, Lin P, Lowengrub J, Wise SM. Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier-Stokes-Cahn-Hilliard system: primitive variable and projection-type schemes. Computational Methods Applied. Mech. Eng., 2017, 326: 144-74

[27] Hairer E, Wanner G. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2010, Springer, Berlin, Springer Series in Computational Mathematics

[28] Hilfer R. Applications of Fractional Calculus in Physics, 2000, River Edge, NJ, World Scientific Publishing Company

[29] Hochbruck M, Ostermann A. Exponential Runge-Kutta methods for parabolic problems. Appl. Numer. Math., 2005, 53: 323-339

[30] Ilic M, Liu F, Turner I, Anh V. Numerical approximation of a fractional-in-space diffusion equation (II) with nonhomogeneous boundary conditions. Fractional Calculus and Applied Analysis, 2006, 9(2): 333-349

[31] Iyiola OS, Asante-Asamani EO, Furati KM, Khaliq AQM, Wade BA. Efficient time discretization scheme for nonlinear space fractional reaction-diffusion equations. Int. J. Comput. Math., 2018, 95(6/7): 1274-1291

[32] Iyiola OS, Wade BA. Exponential integrator methods for systems of non-linear space-fractional models with super-diffusion processes in pattern formation. Comput. Math. Appl., 2018, 75(10): 3719-3736

[33] Khaliq A, Liang X, Furati KM. A fourth-order implicit-explicit scheme for the space fractional nonlinear Schrödinger equations. Numer. Algorithms, 2017, 75: 147-172

[34] Lee, C., Park, J., Kwak, S., Kim, S., Choi, Y., Ham, S., Kim, J.: An adaptive time-stepping algorithm for the Allen-Cahn equation. Journal of Function Spaces 2022, Article ID 2731593 (2022)

[35] Li C, Chen A. Numerical methods for fractional partial differential equations. Int. J. Comput. Math., 2018, 95(6/7): 1048-1099

[36] Meerschaert MM, Tadjeran C. Finite difference approximations for fractional advection-dispersion flow equations. J. Comput. Appl. Math., 2004, 172: 65-77

[37] Norsett SP. Restricted Padé approximation to the exponential function. SIAM Journal of Numerical Analysis, 1978, 15(5): 1008-1029

[38] Norsett SP, Wolfbrandt A. Attainable order of rational approximations to the exponential function with only real poles. BIT Numer. Math., 1977, 17: 200-208

[39] Ortigueira MD. Riesz potential operators and inverses via fractional centered derivatives. Int. J. Math. Math. Sci., 2006, 2006: 48391

[40] Rahman M, Mahmood A, Younis M. Improved and more feasible numerical methods for Riesz space fractional partial differential equations. Appl. Math. Comput., 2014, 237: 264-273

[41] Rida SZ, Yahya AA, Zidan NA, Bakry HM. Fractional order of mathematical systems for some biochemical applications. J. Fract. Calc. Appl., 2014, 5(3S): 1-17

[42] Salkuyeh DK. On the finite difference approximation to the convection diffusion equation. Appl. Math. Comput., 2006, 179: 79-86

[43] Scalas, E.: The application of continuous-time random walks in finance and economics. PhysicaA 362(2), 225–239 (2006)

[44] Singh J. Analysis of fractional blood alcohol model with composite fractional derivative. Chaos, Solitons Fractals, 2020, 140: 110-127

[45] Thomée V. Galerkin Finite Element Methods for Parabolic Problems, 2006, Springer, Berlin, Springer Series in Computational Mathematics25

[46] Voss DA, Khaliq A. A linearly implicit predictor-corrector method for reaction-diffusion equations. Comput. Math. Appl., 1999, 38(11): 207-216

[47] Yang Q, Liu F, Turner I. Numerical methods for fractional partial differential equations with Riesz space fractional derivatives. Appl. Math. Model., 2010, 34(1): 200-218

[48] Yousuf M, Furati KM, Khaliq A. High-order time-stepping methods for two-dimensional Riesz fractional nonlinear reaction-diffusion equations. Comput. Math. Appl., 2020, 80: 204-226

[49] Yue P, Zhou C, Feng JJ, Ollivier-Gooch CF, Hu HH. Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing. J. Comput. Phys., 2006, 219: 47-67

[50] Zhai SY, Gui DW, Zhao JP, Feng XL. High accuracy spectral method for the space-fractional diffusion equations. J. Math. Study, 2014, 47(3): 274-286

[51] Zheng, M., Jin, Z., Liu, F., Anh, V.: Matrix transfer technique for anomalous diffusion equation involving fractional Laplacian. Appl. Numer. Math. 172, 242–258 (2022)

Rights & permissions

Shanghai University