Forward invariant set preservation in discrete dynamical systems and numerical schemes for ODEs: application in biosciences
Tóm tắt
We present two results on the analysis of discrete dynamical systems and finite difference discretizations of continuous dynamical systems, which preserve their dynamics and essential properties. The first result provides a sufficient condition for forward invariance of a set under discrete dynamical systems of specific type, namely time-reversible ones. The condition involves only the boundary of the set. It is a discrete analog of the widely used tangent condition for continuous systems (viz. the vector field points either inwards or is tangent to the boundary of the set). The second result is nonstandard finite difference (NSFD) scheme for dynamical systems defined by systems of ordinary differential equations. The NSFD scheme preserves the hyperbolic equilibria of the continuous system as well as their stability. Further, the scheme is time reversible and, through the first result, inherits from the continuous model the forward invariance of the domain. We show that the scheme is of second order, thereby solving a pending problem on the construction of higher-order nonstandard schemes without spurious solutions. It is shown that the new scheme applies directly for mass action-based models of biological and chemical processes. The application of these results, including some numerical simulations for invariant sets, is exemplified on a general Susceptible-Infective-Recovered/Removed (SIR)-type epidemiological model, which may have arbitrary large number of infective or recovered/removed compartments.
Tài liệu tham khảo
Anguelov, R., Berge, T., Chapwanya, M., Djoko, J.K., Kama, P., Lubuma, J.M.-S., Terefe, Y.: Nonstandard finite difference method revisited and application to the Ebola virus disease dynamics transmission. J. Differ. Equ. Appl. 26(6), 818–854 (2020)
Anguelov, R., Dumont, Y., Lubuma, J.M.-S., Shillor, M.: Dynamically consistent nonstandard finite difference schemes for epidemiological models. J. Comput. Appl. Math. 255, 161–182 (2014)
Anguelov, R., Kama, P., Lubuma, J.M.-S.: On nonstandard finite difference models of reaction-diffusion equations. J. Comput. Appl. Math. 175, 11–29 (2005)
Anguelov, R., Lubuma, J.M.-S.: Contributions to the mathematics of the nonstandard finite difference method and applications. Numer. Methods Partial Differ. Equ. 17, 518–543 (2001)
Anguelov, R., Lubuma, J.M.-S.: A second-order nonstandard finite difference scheme and application to model of biological and chemical processes. In: Gumel, A.B. (ed.) Mathematical and Computational Modeling of Phenomena Arising in Population Biology and Nonlinear Oscillations. AMS Contemporary Mathematics (in press)
Anguelov, R., Lubuma, J.M.-S., Shillor, M.: Topological dynamic consistency of nonstandard finite difference schemes for dynamical systems. J. Differ. Equ. Appl. 17(12), 1769–1791 (2011)
Bony, J.M.: Principe du maximum, inégalité de Harnack et unicité du problème de Cauchy pour les opérateurs elliptiques dégénérés. Ann. Inst. Fourier 19(1), 277–304 (1969)
Brezis, H.: On a characterization of flow-invariant sets. Commun. Pure Appl. Math. XXIII, 261–263 (1970)
Castillo-Chavez, C., Hethcote, H.W., Andreasen, V., Levin, S.A., Liu, W.M.: Epidemiological models with age structure, proportionate mixing, and cross-immunity. J. Math. Biol. 27, 233–258 (1989)
Crandall, M.G.: A generalization of Peano’s existence theorem and flow invariance. Proc. Am. Math. Soc. 36(1), 151–155 (1972)
Gumel, A.B., Connel McCluskey, C., Watmough, J.: An sveir model for assessing potential impact of an imperfect anti-SARS vaccine. Math. Biosci. Eng. 3(3), 485–512 (2006)
Hairer, E., Lubich, C., Wanner, G.: Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations. Springer, Berlin (2006)
Hartman, P.: On invariant sets and on a theorem of Wazewski. Proc. Am. Math. Soc. 32(2), 511–520 (1972)
Hethcote, H.W.: The mathematics of infectious disease. SIAM Rev. 42, 599–653 (2000)
Kojouharov, H.V., Roy, S., Gupta, M., Alalhareth, F., Slezak, J.M.: A second-order modified nonstandard theta method for autonomous differential equations. Appl. Math. Lett. 112, 106775 (2021)
Lubuma, J.M.-S., Roux, A.: An improved theta method for systems of ordinary differential equations. J. Differ. Equ. Appl. 9, 1023–1035 (2003)
Luenberger, D.G.: Optimization by Vector Space Methods. Wiley, New York (1969)
Mickens, R.E.: Nonstandard Finite Difference Models of Differential Equations. World Scientific, Singapore (1994)
Mickens, R.E.: Nonstandard Finite Difference Schemes: Methodology and Applications. World Scientific, Singapore (2021)
Moya, E.D., Rodrigues, D.S., Pietrus, A., Severo, A.M.: A mathematical model for HIV/AIDS under pre-exposure and post-exposure prophylaxis. Biomathematics 11, 2208319 (2022). https://doi.org/10.55630/j.biomath.2022.08.3191/28
Nagumo, M.: Über die Lage der Integralkurven gewöhnlicher Differentialgleichungen. In: Proceedings of the Physico-Mathematical Society of Japan, 3rd Series, vol. 24, pp. 551–559 (1942)
Otunuga, O.M., Ogunsolu, M.O.: Qualitative analysis of a stochastic SEITR epidemic model with multiple stages of infection and treatment. Infect. Dis. Model. 5, 61–90 (2020)
Redheffer, R.M.: The theorems of Bony and Brezis on flow-invariant sets. Am. Math. Mon. 79(7), 740–747 (1972)
Smith, H.L.: Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems. Am. Math. Soc., Providence (1995)
Stuart, A.M., Humphries, A.R.: Dynamical Systems and Numerical Analysis. Cambridge University Press, New York (1998)
van den Driessche, P., Watmough, J.: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 180, 29–48 (2002)
Walter, W.: Ordinary Differential Equations. Springer, New York (1998)
Yamazaki, K., Wang, X.: Global well-posedness and asymptotic behaviour of solutions to a reaction-convective-diffusion cholera epidemic model. Discrete Contin. Dyn. Syst., Ser. B 21, 1297–1316 (2016)