Analysis and Nonstandard Numerical Design of a Discrete Three-Dimensional Hepatitis B Epidemic Model

Abstract

In this work, we numerically investigate a three-dimensional nonlinear reaction-diffusion susceptible-infected-recovered hepatitis B epidemic model. To that end, the stability and bifurcation analyses of the mathematical model are rigorously discussed using the Routh–Hurwitz condition. Numerically, an efficient structure-preserving nonstandard finite-difference time-splitting method is proposed to approximate the solutions of the hepatitis B model. The dynamical consistency of the splitting method is verified mathematically and graphically. Moreover, we perform a mathematical study of the stability of the proposed scheme. The properties of consistency, stability and convergence of our technique are thoroughly analyzed in this work. Some comparisons are provided against existing standard techniques in order to validate the efficacy of our scheme. Our computational results show a superior performance of the present approach when compared against existing methods available in the literature.