Positivity preserving operator splitting nonstandard finite difference methods for SEIR reaction diffusion model

Abstract In this work, we will introduce two novel positivity preserving operator splitting nonstandard finite difference (NSFD) schemes for the numerical solution of SEIR reaction diffusion epidemic model. In epidemic model of infection diseases, positivity is an important property of the continuous system because negative value of a subpopulation is meaningless. The proposed operator splitting NSFD schemes are dynamically consistent with the solution of the continuous model. First scheme is conditionally stable while second operator splitting scheme is unconditionally stable. The stability of the diffusive SEIR model is also verified numerically with the help of Routh-Hurwitz stability condition. Bifurcation value of transmission coefficient is also carried out with and without diffusion. The proposed operator splitting NSFD schemes are compared with the well-known operator splitting finite difference (FD) schemes.

[1]  Herbert W. Hethcote,et al.  Optimal ages of vaccination for measles , 1988 .

[2]  M. Samsuzzoha.  A study on numerical solutions of epidemic models , 2012 .

[3]  Francisco M. Fernández,et al.  On some approximate methods for nonlinear models , 2009, Appl. Math. Comput..

[4]  Ronald E. Mickens.  A Nonstandard Finite Difference Scheme for an Advection-Reaction Equation , 2004 .

[5]  R. Mickens.  Nonstandard Finite Difference Models of Differential Equations , 1993 .

[6]  D. Edwards,et al.  Lead-Acid Battery Model Under Discharge With a Fast Splitting Method , 2011, IEEE Transactions on Energy Conversion.

[7]  M. O. Ahmad,et al.  Numerical Modeling of SEIR Measles Dynamics with Diffusion , 2018 .

[8]  Hanquan Wang,et al.  Numerical studies on the split-step finite difference method for nonlinear Schrödinger equations , 2005, Appl. Math. Comput..

[9]  M. Rafiq.  NUMERICAL MODELING OF INFECTIOUS DISEASES DYNAMICS , 2017 .

[10]  R. Harwood.  Operator Splitting Method and Applications for Semilinear Parabolic Partial Differential Equations , 2011 .

[11]  Kazeem Oare Okosun,et al.  Global stability analysis and control of leptospirosis , 2016 .

[12]  Afia Naheed.  Effect of Treatment on Transmission Dynamics of SARS Epidemic , 2016 .

[13]  Wei Zhang,et al.  Hidden hyperchaos and electronic circuit application in a 5D self-exciting homopolar disc dynamo. , 2017, Chaos.

[14]  C. Paquet,et al.  Measles vaccine effectiveness in standard and early immunization strategies, Niger, 1995. , 1998, The Pediatric infectious disease journal.

[15]  Settapat Chinviriyasit,et al.  Numerical modelling of an SIR epidemic model with diffusion , 2010, Appl. Math. Comput..

[16]  Julien Clinton Sprott,et al.  Detecting Hidden Chaotic Regions and Complex Dynamics in the Self-Exciting Homopolar Disc Dynamo , 2017, Int. J. Bifurc. Chaos.

[17]  T. Fujimoto,et al.  TWO CHARACTERIZATIONS OF INVERSE-POSITIVE MATRICES: THE HAWKINS-SIMON CONDITION AND THE LE CHATELIER-BRAUN PRINCIPLE ∗ , 2004 .

[18]  Vadim Zharnitsky.  Averaging for split-step scheme , 2003 .

[19]  Edward H. Twizell,et al.  An unconditionally convergent discretization of the SEIR model , 2002, Math. Comput. Simul..

[20]  J. Mendlewicz,et al.  PRIMARY ENURESIS AND AFFECTIVE ILLNESS , 1974 .

[21]  M. Rafiq,et al.  Numerical modeling of susceptible latent breaking-out quarantine computer virus epidemic dynamics , 2018, Heliyon.

[22]  Edward H. Twizell,et al.  One-dimensional measles dynamics , 2004, Appl. Math. Comput..

[23]  I SALK.  Epidemic disease. , 1963, Indian medical journal.

[24]  Ronald E. Mickens,et al.  A nonstandard finite difference scheme for a Fisher PDE having nonlinear diffusion , 2003 .

[25]  Matjaz Perc,et al.  Bifurcation analysis of two disc dynamos with viscous friction and multiple time delays , 2019, Appl. Math. Comput..

[26]  M. Yao,et al.  Study of hidden attractors, multiple limit cycles from Hopf bifurcation and boundedness of motion in the generalized hyperchaotic Rabinovich system , 2015 .

[27]  E. A. Bakare,et al.  Modelling and Simulation of the Dynamics of the Transmission of Measles , 2012 .

[28]  Edward H. Twizell,et al.  Chaos-free numerical solutions of reaction-diffusion equations , 1990, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[29]  Kalyan Manna,et al.  Global stability and a non-standard finite difference scheme for a diffusion driven HBV model with capsids , 2015 .

[30]  M. Rafiq,et al.  Numerical modeling of three dimensional Brusselator reaction diffusion system , 2019, AIP Advances.

[31]  Manmohan Singh,et al.  Modelling of tumor cells regression in response to chemotherapeutic treatment , 2017 .

[32]  Manmohan Singh,et al.  Predator-prey model with prey-taxis and diffusion , 2007, Math. Comput. Model..

[33]  Manmohan Singh,et al.  Numerical study of SARS epidemic model with the inclusion of diffusion in the system , 2014, Applied Mathematics and Computation.

[34]  Afia Naheed.  A study of spatio-temporal spread of infectious disease: SARS , 2016 .

[35]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[36]  Xiaohua Ding,et al.  A non-standard finite difference method for a hepatitis B virus infection model with spatial diffusion , 2014 .