A dynamically consistent computational method to solve numerically a mathematical model of polio propagation with spatial diffusion. (May 2022)

• Record Type: Journal Article
• Authors: Ahmed, Nauman; Macías-Díaz, Jorge E.; Shahid, Naveed; Raza, Ali; Rafiq, Muhammad
• Is Part Of: Computer methods and programs in biomedicine. Volume 218(2022)
• Publication Date: 2022-05
• Subjects: Diffusive mathematical model -- Basic reproductive number -- Global stability analysis -- Positivity-preserving scheme -- Neumann stability analysis; 65M06 -- 39A14 -- 35L53 -- 92D25; Medicine -- Computer programs -- Periodicals; Biology -- Computer programs -- Periodicals; Computers -- Periodicals; Medicine -- Periodicals; Médecine -- Logiciels -- Périodiques; Biologie -- Logiciels -- Périodiques; Biology -- Computer programs; Medicine -- Computer programs; Periodicals; Electronic journals; 610.28
• DOI: 10.1016/j.cmpb.2022.106709
• Languages: English
• ISSNs: 0169-2607

Abstract

Highlights: A diffusive deterministic model for the propagation of poliovirus is proposed. The existence of positive and bounded solutions is proved mathematically. The stability analysis of the equilibrium solutions is performed. A positivity-preserving numerical model to solve this system is introduced. Simulations confirm the validity of the analytical and computational results. Abstract: Background and objective: In this work, a mathematical model based on differential equations is proposed to describe the propagation of polio in a human population. The motivating system is a compartmental nonlinear model which is based on the use of ordinary differential equations and four compartments, namely, susceptible, exposed, infected and vaccinated individuals. Methods: In this manuscript, the mathematical model is extended in order to account for spatial diffusion in one dimension. Nonnegative initial conditions are used, and we impose homogeneous Neumann conditions at the boundary. We determine analytically the disease-free and the endemic equilibria of the system along with the basic reproductive number. Results: We establish thoroughly the nonnegativity and the boundedness of the solutions of this problem, and the stability analysis of the equilibrium solutions is carried out rigorously. In order to confirm the validity of these results, we propose an implicit and linear finite-difference method to approximate the solutions of the continuous model. Conclusions: The numerical model is stable in the sense of von Neumann, it yields consistent approximations to the exact solutions of the differential problem, and that it is capable of preserving unconditionally the positivity of the approximations. For illustration purposes, we provide some computer simulations that confirm some theoretical results derived in the present manuscript.