A new eco-epidemiological model with diffusion and square root functional response

Document Type : Original Paper

Authors

Department of Mathematical sciences, Yazd University, Yazd, Iran

Abstract

The aim of this study is to investigate an eco-epidemiological model with a bilinear incidence rate and a square root functional response. First, we examine the equilibrium and stability points of the system for various parameter values. The main challenge in population models is finding a numerical method for approximating non-negative solutions. Some numerical methods, such as the Euler method, are inefficient as they sometimes fail to produce non-negative solutions. Non-negative approximations obtained from non-standard finite difference methods are also conditional. In this paper, we propose a numerical method that provides unconditional and acceptable solutions. We then discuss the compatibility of the proposed numerical method. Finally, we compare the efficiency of the proposed method with the Euler and non standard methods using numerical simulations. We also investigate the effect of the hunting behavior of the prey species on the population density and analyze the dynamic model for some numerical values of the parameters involved in the problem.

Keywords

Main Subjects


[2] Ahmed N., Wei Z., Baleanu D., Rafiq M. and Rehman M. A., Spatio-temporal numerical modeling of reaction-diffusion measles epidemic system, Chaos: An Interdisciplinary Journal of Nonlinear Science, 29(10) (2019), 103101.
[3] Ahmed N., SS T., Rafiq M., Rehman M. A., Ali M., and Ahmad M. O., Positivity preserving operator splitting nonstandard finite difference methods for SEIR reaction diffusion model, Open Mathematics, 17(1) (2019), 313-330.
[4] Ahmed N., Wei Z., Baleanu D., Rafiq M., Rehman M. A., Spatio-temporal numerical modeling of
reaction-diffusion measles epidemic system, Chaos: An Interdisciplinary Journal of Nonlinear Science, 29(10) (2019), 103101.
[5] Ajraldi V., Pittavino M. and Venturino E., Modeling herd behavior in population systems, Nonlinear Analysis: Real World Applications, 12(4) (2011), 2319-2338.
[6] Banerjee, M., Kooi, B. W., Venturino, E., An ecoepidemic model with prey herd behavior and predator feeding saturation response on both healthy and diseased prey, Mathematical Modelling of Natural Phenomena, 12(2) (2017), 133-161.
[7] Braza P. A., Predator–prey dynamics with square root functional responses, Nonlinear Analysis: Real World Applications, 13(4) (2012), 1837-1843.
[8] Bate A. M. and Hilker F. M., Preytaxis and travelling waves in an eco-epidemiological model, Bulletin of mathematical biology, 81 (2019) 995-1030.
[9] Baurmann M., Gross T. and Feudel U., Instabilities in spatially extended predator–prey systems:
Spatio-temporal patterns in the neighborhood of Turing–Hopf bifurcations, Journal of Theoretical
Biology, 245(2) (2007), 220-229.
[10] Bendahmane M. and Langlais M., A reaction-diffusion system with cross-diffusion modeling the spread of an epidemic disease, Journal of Evolution Equations, 10 (2010), 883-904.
[11] Chen-Charpentier B. M. and Kojouharov H. V., An unconditionally positivity preserving scheme for advection–diffusion reaction equations, Mathematical and computer modelling, 57(9-10) (2013),2177-2185.
[12] Gimmelli G., Kooi B. W. and Venturino E., Ecoepidemic models with prey group defense and feeding saturation, Ecological Complexity, 22 (2015), 50-58.
[13] Greenhalgh D., Khan Q. J. and Al-Kharousi F. A., Eco-epidemiological model with fatal disease in the prey, Nonlinear Analysis: Real World Application, 53 (2020), 103072.
[14] Haque M., A detailed study of the Beddington–DeAngelis predator–prey model, Mathematical Biosciences, 234(1) (2011), 1-16.
[15] Hillen T. and Painter K. J., A user’s guide to PDE models for chemotaxis, Journal of mathematical biology, 58(1-2) (2009), 183.
[16] Holling C. S., The components of predation as revealed by a study of small-mammal predation of the European Pine Sawfly1, The Canadian Entomologist, 91(5) (1959), 293-320.
[17] Holling C. S., Some characteristics of simple types of predation and parasitism1, The canadian entomologist, 91(7) (1959), 385-398.
[18] Kim S., Antwi-Fordjour K., Prey group defense to predator aggregated induced fear, The European Physical Journal Plus, 137(6) (2022), 1-17.
[19] Hsu S. B., Hwang T. W. and Kuang Y., Global analysis of the Michaelis–Menten-type ratio-dependent predator-prey system, Journal of mathematical biology, 42 (2001), 489-506.
[20] Kiss K. and Tóth J., n-Dimensional ratio-dependent predator-prey systems with memory, Differential Equations and Dynamical Systems, 17(1-2) (2009), 17-35.
[21] Kooij R. E. and Zegeling A., A predator–prey model with Ivlev’s functional response Journal of Mathematical Analysis and Applications, 198(2) (1996), 473-489.
[22] Lin C. M. and Ho C. P., Local and global stability for a predator–prey model of modified Leslie–Gower and Holling-type II with time-delay, Tunghai Sci, 8 (2006), 33-61.
[23] Lotka A. J., Elements of mathematical biology, Dover Publications, (1956) [24] Liu H. and Cheng H., Dynamic analysis of a prey–predator model with state-dependent control strategy and square root response function, Advances in Difference Equations, 1 (2018), 1-13.
[25] Mickens R.E., Exact solutions to a finite‐difference model of a nonlinear reaction‐advection equation: Implications for numerical analysis, Numerical Methods for Partial Differential Equations, 5(1989), 313-325.
[26] Mickens R. E., A best finite‐difference scheme for the fisher equation, Numerical Methods for Partial Differential Equations, 10(5) (1994), 581-585.
[27] Pal D., Santra P. and Mahapatra G. S., Predator–prey dynamical behavior and stability analysis with square root functional response, International Journal of Applied and Computational Mathematics, (3)(2017), 1833-1845.
[28] Pal A. K and Samanta G. P., Stability analysis of an eco-epidemiological model incorporating a prey refuge, Nonlinear Analysis: Modelling and Control, 15 (2010), 473-491.
[29] Rahman M. S. and Chakravarty S., A predator-prey model with disease in prey, Modelling and Control, 18 (2013), 191-209.
[30] Salman S. M., Yousef A. M. and Elsadany A. A., Stability, bifurcation analysis and chaos control of a discrete predator-prey system with square root functional response, Chaos, Solitons and Fractals, 93(2016), 20-31.
[31] Sarwardi S., Haque M. and Mandal P. K., Persistence and global stability of Bazykin predator–
prey model with Beddington–DeAngelis response function, Communications in Nonlinear Science and Numerical Simulation, 19(1) (2014), 189-209.
[32] Sapoukhina N., Tyutyunov Y. and Arditi R., The role of prey taxis in biological control: a spatial
theoretical model, The American Naturalist, 162(1) (2003), 61-76.
[33] Smith G. D., Numerical solution of partial differential equations: finite difference methods, Oxford university press, (1985)
[34] Sieber M., Malchow H. and Hilker F. M., Disease-induced modification of prey competition in ecoepidemiological models, MoEcological Complexity, 18 (2014), 74-82.
[35] Tan Y., Cai Y., Yao R., Hu M. and Wang W., Complex dynamics in an eco-epidemiological model with the cost of anti-predator behaviors, Nonlinear Dynamics, (2022), 1-15.
[36] Twizell E. H., Wang Y., and Price W. G., Chaos-free numerical solutions of reaction-diffusion equations, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 430(1880)(1990), 541-576.
[37] Wang J. and Nie H., Invasion dynamics of a predator-prey system in closed advective environments, Journal of Differential quations, 318 (2022), 298-322.