Fully-Discrete Lyapunov Consistent Discretizations for Parabolic Reaction-Diffusion Equations with r Species

Abstract

Reaction-diffusion equations model various biological, physical, sociological, and environmental phenomena. Often, numerical simulations are used to understand and discover the dynamics of such systems. Following the extension of the nonlinear Lyapunov theory applied to some class of reaction-diffusion partial differential equations (PDEs), we develop the first fully discrete Lyapunov discretizations that are consistent with the stability properties of the continuous parabolic reaction-diffusion models. The proposed framework provides a systematic procedure to develop fully discrete schemes of arbitrary order in space and time for solving a broad class of equations equipped with a Lyapunov functional. The new schemes are applied to solve systems of PDEs, which arise in epidemiology and oncolytic M1 virotherapy. The new computational framework provides physically consistent and accurate results without exhibiting scheme-dependent instabilities and converging to unphysical solutions. The proposed approach represents a capstone for developing efficient, robust, and predictive technologies for simulating complex phenomena.

Keywords

Reaction-diffusion / Lyapunov functional / Fully discrete Lyapunov consistent discretizations / Epidemiology / Tumor growth / Oncolytic M1 virotherapy / 93D05 / 65L06 / 35K57 / 65N12 / 65N35 / 92D30

References

[1] Abdelmalek S, Bendoukha S. Global asymptotic stability for a SEI reaction-diffusion model of infectious diseases with immigration. Int. J. Biomath., 2018, 11(03): 1850044

[2] Al Jahdali R, Dalcin L, Parsani M. On the performance of relaxation and adaptive explicit Runge-Kutta schemes for high-order compressible flow simulations. J. Comput. Phys., 2022, 464 111333

[3] Al Jahdali, R., Fernández, D., Dalcin, L., Sayyari, M., Markowich, P., Parsani, M.: Reproducibility repository for “Fully discrete Lyapunov consistent discretizations of any order for parabolic reaction-diffusion equations with r species”. (2023). https://github.com/aanslab-papers/2023-fully-discrete-lyapunov-consistent-discretizations

[4] Anosov DV, Aranson SK, Arnold VI, Bronshtein IU, Il’yashenko YS, Grines VZ. Ordinary Differential Equations and Smooth Dynamical Systems, 1997, Berlin, Heidelberg, Springer

[5] Arnold DN. An interior penalty finite element method with discontinuous elements. SIAM J. Numer. Anal., 1982, 194742-760

[6] Biazar J, Mehrlatifan MB. A compact finite difference scheme for reaction-convection-diffusion equation. Chiang Mai J. Sci., 2018, 45(3): 1559-1568

[7] Brezzi, F., Ushiki, S., Fujii, H.: “Real” and “ghost” bifurcation dynamics in difference schemes for ODEs. In: Numerical Methods for Bifurcation Problems, pp. 79–104. Springer, Basel (1984)

[8] Bueno-Orovio A, Perez-Garcia VM, Fenton FH. Spectral methods for partial differential equations in irregular domains: the spectral smoothed boundary method. SIAM J. Sci. Comput., 2006, 28(3): 886-900

[9] Buonomo B, Rionero S. On the Lyapunov stability for SIRS epidemic models with general nonlinear incidence rate. Appl. Math. Comput., 2010, 217(8): 4010-4016

[10] Burton, T.A.: Stability by Fixed Point Theory for Functional Differential Equations. Courier Corporation, Chicago (2013)

[11] Butcher, J.C.: Numerical Methods for Ordinary Differential Equations. John Wiley & Sons Ltd (2008)

[12] Cai Y, Chi D, Liu W, Wang W. Stationary patterns of a cross-diffusion epidemic model. Abstr. Appl. Anal., 2013

[13] Carpenter, M.H., Gottlieb, D., Abarbanel, S.: Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes. J. Comput. Phys. 111(2), 220–236 (1994). https://doi.org/10.1006/jcph.1994.1057

[14] Carpenter MH, Fisher TC, Nielsen EJ, Frankel SH. Entropy stable spectral collocation schemes for the Navier-Stokes equations: discontinuous interfaces. SIAM J. Sci. Comput., 2014, 36(5): 835-867

[15] Carpenter MH, Nordström J, Gottlieb D. A stable and conservative interface treatment of arbitrary spatial accuracy. J. Comput. Phys., 1999, 148(2): 341-365

[16] Chernyshenko AY, Olshanskii MA. An adaptive octree finite element method for PDEs posed on surfaces. Comput. Methods Appl. Mech. Eng., 2015, 291: 146-172

[17] Cockburn B, Shu C-W. The local discontinuous Galerkin method for time-dependent convection-diffusion systems. SIAM J. Numer. Anal., 1998, 35(6): 2440-2463

[18] Curtain R, Zwart H. Introduction to Infinite-Dimensional Systems Theory: a State-Space Approach, 2020, New York, Springer

[19] David, A.: The simulations driving the world’s response to COVID-19: how epidemiologists rushed to model the coronavirus pandemic. Nature 580(7803), 316–318 (2020)

[20] Del Rey Fernández, D.C., Carpenter, M.H., Dalcin, L., Zampini, S., Parsani, M.: Entropy stable h/p\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$h/p$$\end{document}-nonconforming discretization with the summation-by-parts property for the compressible Euler and Navier-Stokes equations. SN Partial Differential Equations and Applications 1, 9 (2020). https://doi.org/10.1007/s42985-020-00009-z

[21] Del Rey Fernández DC, Hicken JE, Zingg DW. Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations. Computers & Fluids, 2014, 95(22): 171-196

[22] Elaiw AM, Hobiny AD, Al Agha AD. Global dynamics of reaction-diffusion oncolytic M1 virotherapy with immune response. Appl. Math. Comput., 2020, 367 124758

[23] Gatenby RA, Gawlinski ET. A reaction-diffusion model of cancer invasion. Can. Res., 1996, 56(24): 5745-5753

[24] Greenberg JM, Ta Tsien L. The effect of boundary damping for the quasilinear wave equation. J. Differential Equations, 1984, 52(1): 66-75

[25] Griffiths D, Sweby P, Yee HC. On spurious asymptotic numerical solutions of explicit Runge-Kutta methods. IMA J. Numer. Anal., 1992, 12(3): 319-338

[26] Gustafsson, B., Kreiss, H.O., Oliger, J.: Time-Dependent Problems and Difference Methods. Pure and Applied Mathematics: a Wiley Series of Texts, Monographs and Tracts. John Wiley & Sons Ltd (2013)

[27] Haddad, W.M., Chellaboina, V.: Nonlinear Dynamical Systems and Control: a Lyapunov-based Approach. Princeton University Press, Princeton, New Jersey (2008)

[28] Hairer E, Iserles A, Sanz-Serna JM. Equilibria of Runge-Kutta methods. Numer. Math., 1990, 58(1): 243-254

[29] Haushofer J, Metcalf CJE. Which interventions work best in a pandemic?. Science, 2020, 368(6495): 1063-1065

[30] Heidari, M., Ghovatmand, M., Skandari, M.N., Baleanu, D.: Numerical solution of reaction-diffusion equations with convergence analysis. Journal of Nonlinear Mathematical Physics, 30, 384–399 (2022). https://doi.org/10.1007/s44198-022-00086-1

[31] Iserles A. Stability and dynamics of numerical methods for nonlinear ordinary differential equations. IMA J. Numer. Anal., 1990, 10: 1-30

[32] Ketcheson DI. Relaxation Runge-Kutta methods: conservation and stability for inner-product norms. SIAM J. Numer. Anal., 2019, 57(6): 2850-2870

[33] Khalil, H.K.: Nonlinear Systems, 3rd edn. Prentice-Hall, Upper Saddle River, NJ (2002)

[34] Kolev MK, Koleva MN, Vulkov LG. An unconditional positivity-preserving difference scheme for models of cancer migration and invasion. Mathematics, 2022, 10(1): 131

[35] Komornik V. Exact Controllability and Stabilization: the Multiplier Method, 1994, Chichester, John and Wiley

[36] Kreiss, H.-O., Scherer, G.: Finite element and finite difference methods for hyperbolic partial differential equations. In: Boor, C. (ed.) Mathematical Aspects of Finite Elements in Partial Differential Equations. Academic Press (1974)

[37] Krstic M, Guo B-Z, Balogh A, Smyshlyaev A. Output-feedback stabilization of an unstable wave equation. Automatica, 2008, 44163-74

[38] Krstic M, Smyshlyaev A. Backstepping boundary control for first-order hyperbolic PDEs and application to systems with actuator and sensor delays. Systems & Control Letters, 2008, 57(9): 750-758

[39] LaSalle, J.P., Artstein, Z.: The Stability of Dynamical Systems. CBMS-NSF Regional Conference Series in Applied Mathematics. Society for Industrial and Applied Mathematics, Philadelphia (1976)

[40] LaSalle, J.P., Lefschetz, S.: Stability by Lyapunov’s Direct Method: with Applications. Mathematics in Science and Engineering: a Series of Monographs and Textbooks. Academic Press, New York (1961)

[41] Liu, Y., Cai, J., Liu, W., Lin, Y., Guo, L., Liu, X., Qin, Z., Xu, C., Zhang, Y., Su, X., Deng, K., Yan, G., Liang, J.: Intravenous injection of the oncolytic virus M1 awakens antitumor T cells and overcomes resistance to checkpoint blockade. Cell Death & Disease 11(12), 1062 (2020). https://doi.org/10.1038/s41419-020-03285-0

[42] Lyapunov, A.M.: The general problem of the stability of motion (In Russian). Ph.D. thesis, University of Kharkov (1892)

[43] Lyapunov AM. The general problem of the stability of motion. Int. J. Control, 1992, 55(3): 531-534

[44] Macías-Díaz J. On the numerical and structural properties of a logarithmic scheme for diffusion-reaction equations. Appl. Numer. Math., 2019, 140: 104-114

[45] Macías-Díaz JE, Puri A. An explicit positivity-preserving finite-difference scheme for the classical Fisher-Kolmogorov-Petrovsky-Piscounov equation. Appl. Math. Comput., 2012, 218(9): 5829-5837

[46] Mañosa V. Periodic travelling waves in nonlinear reaction-diffusion equations via multiple Hopf bifurcation. Chaos, Solitons & Fractals, 2003, 18(2): 241-257

[47] Marelli G, Howells A, Lemoine NR, Wang Y. Oncolytic viral therapy and the immune system: a double-edged sword against cancer. Front. Immunol., 2018, 9: 866

[48] Mattsson K. Boundary procedures for summation-by-parts operators. J. Sci. Comput., 2003, 18(1): 133-153

[49] Mattsson K, Ham F, Iaccarino G. Stable boundary treatment for the wave equation on second-order form. J. Sci. Comput., 2009, 41(3): 366-383

[50] Newell AC. Finite amplitude instabilities of partial difference equations. SIAM J. Appl. Math., 1977, 33(1): 133-160

[51] Nordström J, Carpenter MH. Boundary and interface conditions for high-order finite-difference methods applied to the Euler and Navier-Stokes equations. J. Comput. Phys., 1999, 148(2): 621-645

[52] Nordström J, Carpenter MH. High-order finite-difference methods, multidimensional linear problems, and curvilinear coordinates. J. Comput. Phys., 2001, 173(1): 149-174

[53] Parsani, M., Boukharfane, R., Nolasco, I.R., Del Rey Fernández, D.C., Zampini, S., Hadri, B., Dalcin, L.: High-order accurate entropy-stable discontinuous collocated Galerkin methods with the summation-by-parts property for compressible CFD frameworks: scalable ssdc algorithms and flow solver. J. Comput. Phys. 424, 109844 (2021). https://doi.org/10.1016/j.jcp.2020.109844

[54] Parsani M, Carpenter MH, Fisher TC, Nielsen EJ. Entropy stable staggered grid discontinuous spectral collocation methods of any order for the compressible Navier-Stokes equations. SIAM J. Sci. Comput., 2016, 38(5): 3129-3162

[55] Parsani M, Carpenter MH, Nielsen EJ. Entropy stable discontinuous interfaces coupling for the three-dimensional compressible Navier-Stokes equations. J. Comput. Phys., 2015, 290: 132-138

[56] Ranocha, H., Sayyari, M., Dalcin, L., Parsani, M., Ketcheson, D.I.: Relaxation Runge-Kutta methods: fully-discrete explicit entropy-stable schemes for the Euler and Navier-Stokes equations. SIAM J. Sci. Comput. 42(2), A585–C68 (2020). https://doi.org/10.1137/19M1263480

[57] Reitz RD. A study of numerical methods for reaction-diffusion equations. SIAM J. Sci. Stat. Comput., 1981, 2(1): 95-106

[58] Roy CJ. Review of code and solution verification procedures for computational simulation. J. Comput. Phys., 2005, 205(1): 131-156

[59] Sigdel RP, McCluskey CC. Global stability for an SEI model of infectious disease with immigration. Appl. Math. Comput., 2014, 243: 684-689

[60] Smyshlyaev A, Cerpa E, Krstic M. Boundary stabilization of a 1-D wave equation with in-domain antidamping. SIAM J. Control. Optim., 2010, 48(6): 4014-4031

[61] Svärd M, Carpenter MH, Nordström J. A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions. J. Comput. Phys., 2007, 225(1): 1020-1038

[62] Svärd M, Nordström J. Review of summation-by-parts schemes for initial-boundary-value-problems. J. Comput. Phys., 2014, 268(1): 17-38

[63] Thielscher, A., Antunes, A., Saturnino, G.B.: Field modeling for transcranial magnetic stimulation: a useful tool to understand the physiological effects of TMS? In: 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), pp. 222–225 (2015). https://doi.org/10.1109/EMBC.2015.7318340

[64] Thomas PD, Lombard CK. Geometric conservation law and its application to flow computations on moving grids. AIAA J., 1979, 17(10): 1030-1037

[65] Toubaei S, Garshasbi M, Jalalvand M. A numerical treatment of a reaction-diffusion model of spatial pattern in the embryo. Computational Methods for Differential Equations, 2016, 4(2): 116-127

[66] Wang W, Cai Y, Wu M, Wang K, Li Z. Complex dynamics of a reaction-diffusion epidemic model. Nonlinear Anal. Real World Appl., 2012, 13(5): 2240-2258

[67] Wang W, Zhao X-Q. Basic reproduction numbers for reaction-diffusion epidemic models. SIAM J. Appl. Dyn. Syst., 2012, 11(4): 1652-1673

[68] Xiao X, Wang K, Feng X. A lifted local Galerkin method for solving the reaction-diffusion equations on implicit surfaces. Comput. Phys. Commun., 2018, 231: 107-113

[69] Zhang, J., Liu, Y., Tan, J., Zhang, Y., Wong, C.-W., Lin, Z., Liu, X., Sander, M., Yang, X., Liang, L., Song, D., Dan, J., Zhou, Y., Cai, J., Lin, Y., Liang, J., Hu, J., Yan, G., Zhu, W.: Necroptotic virotherapy of oncolytic alphavirus M1 cooperated with Doxorubicin displays promising therapeutic efficacy in TNBC. Oncogene 40(29), 4783–4795 (2021). https://doi.org/10.1038/s41388-021-01869-4

Funding

Global Collaborative Research, King Abdullah University of Science and Technology(BAS/1/1663-01-01)

RIGHTS & PERMISSIONS

The Author(s)