A dynamical study on stochastic reaction diffusion epidemic model with nonlinear incidence rate

Eur Phys J Plus. 2023;138(4):350. doi: 10.1140/epjp/s13360-023-03936-z. Epub 2023 Apr 22.

Abstract

The current study deals with the stochastic reaction-diffusion epidemic model numerically with two proposed schemes. Such models have many applications in the disease dynamics of wildlife, human life, and others. During the last decade, it is observed that the epidemic models cannot predict the accurate behavior of infectious diseases. The empirical data gained about the spread of the disease shows non-deterministic behavior. It is a strong challenge for researchers to consider stochastic epidemic models. The effect of the stochastic process is analyzed. So, the SIR epidemic model is considered under the influence of the stochastic process. The time noise term is taken as the stochastic source. The coefficient of the stochastic term is a Borel function, and it is used to control the random behavior in the solutions. The proposed stochastic backward Euler scheme and the proposed stochastic implicit finite difference scheme (IFDS) are used for the numerical solution of the underlying model. Both schemes are consistent in the mean square sense. The stability of the schemes is proven with Von-Neumann criteria and schemes are unconditionally stable. The proposed stochastic backward Euler scheme converges toward a disease-free equilibrium and does not converge toward an endemic equilibrium but also possesses negative behavior. The proposed stochastic IFD scheme converges toward disease-free equilibrium and endemic equilibrium. This scheme also preserves positivity. The graphical behavior of the stochastic SIR model is much similar to the classical SIR epidemic model when noise strength approaches zero. The three-dimensional plots of the susceptible and infected individuals are drawn for two cases of endemic equilibrium and disease-free equilibriums. The efficacy of the proposed scheme is shown in the graphical behavior of the test problem for the various values of the parameters.