Abstract
In this paper, two structure-preserving nonstandard finite difference (NSFD) operator splitting schemes are designed for the solution of reaction diffusion epidemic models. The proposed schemes preserve all the essential properties possessed by the continuous systems. These schemes are applied on a diffusive SEIQV epidemic model with a saturated incidence rate to validate the results. Furthermore, the stability of the continuous system is proved, and the bifurcation value is evaluated. A comparison is also made with the existing operator splitting numerical scheme. Simulations are also performed for numerical experiments.