Global numerical analysis of an improved IMEX numerical scheme for a reaction diffusion SIS model in advective environments

This paper presents globally numerical properties of a new numerical scheme for a reaction-diffusion advection susceptible-infected-susceptible (SIS) model. A new numerical treatment technique is introduced in spatial discretization of advection-diffusion equation, which enables the numerical solutions to preserve the stability and positivity with less stepsize restrictions. The convergence, biological significance and globally stability of numerical solutions is explored in the paper. A threshold value, named by numerical basic reproduction number and denoted by ������& UDelta;������ 0 , is introduced in the numerical stability analysis of the model. It is proved the numerical disease free equilibrium (DFE) is globally asymptotically stable if ������& UDelta;���0 ���<1and unstable if ������& UDelta;���0 ���>1. It is shown the numerical basic number ������& UDelta;������ 0 replicates the asymptotic behaviors of the basic reproduction number ������0 for the model. Some numerical experiments are given in the end to confirm the conclusions.