Numerical analysis of diffusive susceptible-infected-recovered epidemic model in three space dimension

Abstract

In this article, numerical solution of three dimensional susceptible-infected-recovered (SIR) reaction-diffusion epidemic system is furnished with a time efficient operator splitting nonstandard finite difference (OS-NSFD) method. We perform the comparison of proposed OS-NSFD method with popular forward Euler explicit finite difference (FD) method and time efficient backward Euler operator splitting finite difference (OS-FD) implicit method. The proposed OS-NSFD method is implicit in nature but computationally efficient as compared to forward Euler explicit (FD) scheme. The numerical stability and bifurcation value of transmission coefficient for SIR reaction-diffusion epidemic system is also investigated with the aid of Routh-Hurwitz method. At the end, we give two numerical experiments and simulation. In first experiment, all the numerical schemes are compared with the help of simulations. In second experiment we show the simulations of proposed NSFD technique at different values of parameters. Also we discuss the importance of transmission rate to control the spread of disease with the help of simulations. (C) 2019 Elsevier Ltd. All rights reserved.