Abstract
In this manuscript, a Bazykin–Berezovskaya model with diffusion by strong Allee effects is studied. Neumann boundary conditions are used to see the positive solution of a diffusion system. Local stability analyses are discussed for all the equilibrium points. The analysis of stability for the proposed scheme is also given. Implicit finite difference schemes like: Euler, Crank–Nicolson (CN) and non-standard finite difference (NSFD) are used to verify the simulation by numerically. A comparison reveals that NSFD method is unconditionally stable for any temporal step-size.
1. Introduction
In biomathematics, the Predator-Pray (PP) model is making an impact to discuss and understand the dynamics of populations. In 1925 and 1926, Alferd J. Lotka [Citation1] and Vito Volterra [Citation2] initially developed the model for a polynomial second-degree differential equation known as the Lotka-Volterra (LV) system. Lotka developed and studied the application of insectivore animal and plant family. On the other hand, the same model discussed by Volterra to explain the percentage of caught predictory fish. Volterra considered the system of linear ordinary differential equation by taking the PP model (1) (1) where and ζ are positive numbers. Later on, the generalized form of LV system for the multiple dimensions was (2) (2) Many researchers extended this model for a variety of applications especially in two dimensions, i.e. one predator and other one is pray species [Citation3–7]. The theoretical aspect of PP model is very interesting for the development of population dynamics. In recent years, the authors focused on the phenomenon of the Allee effect. The correlation between the density (population size) and rate of growth is discussed on the Allee effect. Depending upon the nature of population growth, there are two types of the Allee effect, i.e. Strong Allee effect and Weak Allee effect. A weak Allee effect means, the population growth rate has reduced as per capita growth rate. Moreover, a population size rate becomes negative in case of a strong Allee effect under zero density. This complete procedure was developed by Allee in Ref. [Citation8]. The Allee effect can be written as, (3) (3) where δ represents growth rate. The Allee effect depends on the value of d, if d>0, it behaves as a strong Allee effect and its becomes weak if . In literature, there are different works that exist to analyse the outcomes of Allee effect on different population models and conclude that it can have more impact on the dynamical system. The stabilizing or destabilizing of a dynamical system by replacing the stability of some singularities in Allee effects, and the time taken to reach the stable solution of that dynamical system to take much longer.
In 1988, Bazykin and Berezovskaya introduced the new PP model which is known as Bazykin–Berezovskaya (BB) model [Citation9]. This model related to the strong Allee effect. Voorn et al. [Citation10] considered the following Bazykin-Berezovskaya model with strong Allee effects; (4) (4) The model is in dimensionless and X, Y are the sizes of prey and predator population, respectively. The parameter K is the carrying capacity, L is the strong Allee effect threshold, C is the feeding efficiency, M is the predator mortality rate of Lotka-Voltrra model [Citation10,Citation11]. Bifurcation review in (Equation4(4) (4) ) with weak Allee effects discussed by Saad and Boubaker [Citation12]. Further, the same authors studied the fractional order qualitative analysis. The noise-induced eradication in (Equation4(4) (4) ) investigated by Bashkirtseva and Ryashk [Citation13]. In engineering, ecology, biology, and physics, the discrete-time population (DTP) models are applicable [Citation14–21]. In addition, Din [Citation22], Pal et al. [Citation23], and Kartal et al. [Citation24] investigated that (DTP) models may be more complex than continuous time models. Recently different authors published in planar systems [Citation25–29].
The authors of Ref. [Citation30] have used the fractional derivative of polynomials of the first kind in the Caputo sense to obtain a spectral solution for delay differential equations with the help of an explicit Chebyshev tau method. The authors of [Citation31] has used the shifted Gegenbauer-Gauss collocation method for solving fractional differential equations with delay. Youssri et al. [Citation32] has used the shifted Chebeyshev polynomials of the third kind in the Caputo sense to obtain fractional differential equation and fractional delay differential equations. In another article, Youssri et al. [Citation33] has applied generalized Lucas polynomial sequence treatment to obtain the solution of fractional pantograph differential equation. The authors of the papers [Citation34–42] considered various kinds of ordinary and functional differential equations of higher order and Volterra integro-differential equations. They investigated the stability and some other qualitative analysis of solutions of that equations using suitable tools and techniques.
2. Development of diffusion model
Over here, we consider 1-D coupled BB population model: (5) (5) where is a bounded domain in with sufficiently smooth boundary . Here and symbolize the sizes of prey and predator respectively at time , and a point . The parameters and are diffusion rates for each X and Y. Certainly due to biological reasons X and Y have to be non-negative functions. and are also non-negative. We denote as ordinary derivative with respect to τ.
2.1. Steady states of model without diffusion
The system (Equation4(4) (4) ) which is given in the absence of diffusion has four equilibrium points and these are and .
2.2. Diffusionless case
In this case, we have .
For equilibrium point , the Jacobian matrix is: (6) (6) This means the two eigenvalues are and . Since all the parameters of the system are positive, so and . Thus according to the RH criterion, if all the eigenvalues are negative or have negative real parts, then the model is locally asymptotically stable at that equilibrium point. So the model under consideration is locally asymptotically stable at .
For equilibrium point , the Jacobian matrix is: (7) (7) This means the two eigenvalues are and . Since all the parameters of the system are positive, so if K<L and if L<M. Thus according to the RH criterion, if all the eigenvalues are negative or have negative real parts, then the model is locally asymptotically stable at that equilibrium point. So the model is locally asymptotically stable at if K<L<M.
For equilibrium point , the Jacobian matrix is: (8) (8) This means the two eigenvalues are and . Since the all the parameters of the system are positive, so if L<K and if K<M. Thus according to the RH criterion, if all the eigenvalues are negative or have negative real parts, then the model is locally asymptotically stable at that equilibrium point. So the model is locally asymptotically stable at if L<K<M.
For equilibrium point , the Jacobian matrix is: (9) (9) The characteristic equation of it is (10) (10) It can be written as (11) (11) where and . Then according to the RH criterion if and , then the system is locally asymptotically stable at that equilibrium point. It could be possible only if 2M−K−L>0 or , M>L and K>M. In other words, it can be written as and M>L otherwise this model is unstable.
2.3. Diffusion model with steady-states
We have to discuss the basic properties of the Bazykin–Berezovskaya population system of non-homogeneous steady states (12) (12) with the suitable conditions .
Definition 2.1
If a constant solution is stable in the absence of diffusion but becomes unstable when diffusion is present, the constant solution is said to be turing unstable.
Proposition 2.1
[Citation43]
Suppose that and . It follows that
(a) | If (13) (13) with on and , then (14) (14) | ||||
(b) | If (15) (15) with on and , then (16) (16) |
3. Numerical investigation
We used the following numerical analysis based on backward Euler (BE) method, CN method and NSFD method to investigate the solutions and their stability.
3.1. The BB model
In this section, the system (Equation5(5) (5) ) will discretized: (17) (17) Using FD for Equation (Equation17(17) (17) ), firstly are shared into with space and time step sizes like, and . Then the mesh nodes are and , where and .
We can write and as the FD approximations of and . First- and second-order temporal and spatial derivative difference formulas are: (18) (18) (19) (19) (20) (20)
3.2. BE method
Applying the values of and in the first equation of (Equation17(17) (17) ) and the values of and in the second equation of (Equation17(17) (17) ), we have (21) (21) (22) (22) After some simple calculation, we have (23) (23) (24) (24) where and .
This technique is stable for Equations (Equation23(23) (23) ) and (Equation24(24) (24) ).
3.3. CN method
Applying the values of , and in the first equation of (Equation4(4) (4) ) and the values of , and in Equation (Equation5(5) (5) ), we have (25) (25) (26) (26) After some simple calculation, we have (27) (27) (28) (28) where and .
This technique is unconditionally stable for Equations (Equation17(17) (17) ).
3.4. Proposed NSFD method
We will design the NSFD method for system (Equation17(17) (17) ). Applying the values of , , , in Equation (Equation4(4) (4) ), we have (29) (29) (30) (30) or (31) (31) (32) (32) After some calculation, we have the following explicit NSFD scheme (33) (33) where and .
This approach is completely stable for Equations (Equation17(17) (17) ).
3.4.1. Analysis of stability NSFD scheme
To obtain the steadiness norms, Von-Neumann steadiness scheme is used here. After linearizing (Equation31(31) (31) ) and then applying the Von-Neumann steadiness norms, we get (34) (34) (35) (35) (36) (36) (37) (37) Similarly, (38) (38) From (Equation37(37) (37) ) and (Equation38(38) (38) ), it is clear that the proposed nonstandard FD scheme for system (Equation17(17) (17) ) is completely stable.
3.4.2. Consistency of NSFD method
To observe the uniformity of NSFD technique, Taylor series is used for , and are given below: (39) (39) (40) (40) (41) (41) (42) (42) Putting the values of , and in above equation, we have (43) (43) By substituting the values and then taking , the above equation gives (44) (44) Similarly, Taylor series is used for , and are given below: (45) (45) (46) (46) (47) (47) Putting the values of , and in above equation, we have (48) (48) Putting and taking , the above equation simplifies to the equation given below (49) (49)
4. Result and discussion
The initial distribution of the number of Prey and Predator are given in Figure . 3D mesh simulations and 2D solutions for prey and predator using backward Euler scheme are given in Figures – . Simulations 2, 4, 6 and 8 converge to all the four equilibrium points verifying the theoretical part, whereas in the simulations 3, 5, 7 and 9 one or both simulations fails to converge the true equilibrium points when increasing the stepsize . When using Crank–Nicolson scheme, 3D mesh plots and 2D solutions for and are given in Figures – . Simulations 10,12,14, and 16 converge to true equilibrium points, whereas in Figures , , and , one or both populations fails to converge to true equilibrium points when taking temporal stepsize higher than the stepsize taken for 10, 12 14 and 16. When using NSFD scheme for the system (Equation17(17) (17) ), the simulations for both populations converge to the stable solutions even for (Figures –). All the simulations are done on Intel Corei5, 7th generation Dell Laptop.
5. Conclusions
In this manuscript, we have proposed a BB model by strong Allee effects. The main objective is to investigate the effects of diffusion term in the dynamics of the BB model. We have developed three finite difference schemes to solve system (Equation4(4) (4) ). We used Taylor's theorem to check our proposed nonstandard FD scheme is accurate in first order and consistent in spatiotemporal. The analysis of nonstandard FD is done by using the Von-Neumann stability. The graphical portraits exhibits the fact that the proposed nonstandard FD is stable for both small or large temporal step-size. On the other hand, both well-known methods fail to preserve the positivity property and show divergence on different step size . Hence the proposed scheme is more effective, accurate, reliable and consistent for under observed model. Finally, the authors have concluded that NSFD scheme has an edge on the other two schemes.
Acknowledgments
Acknowledgements
All authors contributed equally and significantly in writing this article. All authors read and approved the final manuscript.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Lotka AJ. Elements of physical biology. Baltimore (MD): Williams and Wilkins; 1925.
- Volterra V. Variazioni e fluttuazioni del numero d"Individui in specie animali conviventi. Memoire del R. Comitato talaassografico italiano, men CXXXI; 1927.
- Lou Y, Ni W-M. Diffusion, self-diffusion and cross-diffusion. J Differ Equ. 1996;131:79–131.
- Zafar ZUA, Shah Z, Ali N, et al. Numerical study and stability of the Lengyel Epstein chemical model with diffusion. Adv Differ Equ. 2020;2020:427.
- Zafar ZUA, Ali N, Inc M, et al. Mathematical modeling of corona virus (COVID-19) and stability analysis. Comput Methods Biomech Biomed Eng. 2022:1–20. DOI:10.1080/10255842.2022.2109020.
- Zafar ZUA, Hussain MT, Inc M, et al. Fractional order dynamics of human papillomavirus. Res Phys. 2022;34:Article ID 105281.
- Zafar ZUA, Ali N, Baleanu D. Dynamics and numerical investigations of a fractional-order model of toxoplasmosis in the population of human and cats. Chaos Solitons Fract. 2021;151:Article ID 111261.
- Allee WC, Bowen E. Studies in animal aggregations: mass protection against colloidal silver among goldfishes. J Exp Zool. 1932;61:185–207.
- Bazykin AD. Nonlinear dynamics of interacting populations. Singapore: World Scientific; 1998.
- Akrami MH. Dynamical behaviours of Bazykin-Berezovskayamodel with fractional-order and its discretization. Comput Methods Differ Equ. 2021;9(4):1013–1027. DOI:10.22034/cmde.2020.30802.1460.
- Ben Saad A, Boubaker O. On bifurcation analysis of the predator-prey BB-model with weak allee effect. In: 16th Int. Conf. Sciences and Techniques of Automatic Control and Computer Engineering (STA). IEEE; 2015. p. 19–23.
- Bashkirtseva l, Ryashko L. Noise-induced extinction in Bazykin-Berezovskaya population model. Eur Phys J B. 2016;89:165.
- Ben Saad A, Boubaker O. A new fractional-order predator-prey system with Allee effect. In: Azar, et al., editor. Fractional order control and synchronization of chaotic systems. Studies in computational intelligence. Vol. 688. Springer; 2017. p. 857–877.
- He Z, Lai X. Bifurcation and chaotic behavior of a discrete-time predator-prey system. Nonlinear Anal RWA. 2011;12:403–417.
- Jing Z, Yang J. Bifurcation and chaos in discretetime predator-prey system. Chaos Soliton Fract. 2006;27:259–277.
- Liu X, Xiao D. Complex dynamic behaviors of a discrete-time predator-prey system. Chaos Soliton Fract. 2007;32:80–94.
- Li B, He Z. Bifurcations and chaos in a twodimensional discrete Hindmarsh-Rose model. Nonlinear Dyn. 2014;76(1):697–715.
- Yuan LG, Yang QG. Bifurcation, invariant curve and hybrid control in a discrete-time predator-prey system. Appl Math Model. 2015;39(8):2345–2362.
- Cheng L, Cao H. Bifurcation analysis of a discretetime ratio-dependent predator-prey model with Allee Effect. Commun Nonlinear Sci Numer Simulat. 2016;38:288–302.
- Kartal S. Dynamics of a plant-herbivore model with differential-difference equations. Cogent Math. 2016;3(1):Article ID 1136198.
- Chen B, Chen J. Complex dynamic behaviors of a discrete predator-prey model with stage structure and harvesting. Int J Biomath. 2017;10(1):Article ID 1750013.
- Pal. S. K. Sasmal S, Pal N. Chaos control in a discrete-time predator-prey model with weak Allee effect. Int J Biomath. 2018;11(7):Article ID 1850089.
- Kartal S, Fuat G. Global behaviour of a predator-prey like model with piecewise constant arguments. J Biol Dyn. 2015;9:159–171.
- Din Q. Controlling chaos in a discrete-time prey. predator model with Allee effects. Int J Dyn Control. 2018;6(2):858–872.
- Elabbasy EM, Elsadany AA, Zhang Y. Bifurcation analysis and chaos in a discrete reduced Lorenz system. Appl Math Comput. 2014;228:184–194.
- Atabaigi A, Akrami MH. Dynamics and bifurcations of a hostparasite model. Int J Biomath. 2017;10(6):Article ID 1750089.
- Salman SM, Yousef AM, Elsadany AA. Stability. bifurcation analysis and chaos control of a discrete predator-prey system with square root functional response. Chaos Soliton Fract. 2016;93:20–31.
- Isik S. A study of stability and bifurcation analysis in discrete-time predatorprey system involving the Allee effect. Int J Biomath. 2019;12(1):Article ID 1950011.
- Din Q. Complexity and chaos control in a discretetime prey-predator model. Commun Nonlinear Sci Numer Simul. 2017;49:113–134.
- Abd-Elhameed WM, Machado JAT, Youssri YH. Hypergeometric fractional derivative formula of shifted Chebyshev polynomials: tau algorithm for a type of fractional delay differential equations. Int J Nonlinear Sci Numer Simul. 2022;23(7-8):1253–1268. DOI:10.1515/ijnsns-2020-0124.
- Hafez RM, Youssri YH. Shifed Gegenbaur-Gauss collocation method for solving fractional-differential equations with proportional delays. Kragujevac J Math. 2022;46(6):981–996.
- Youssri YH, Abd-Elhameed WM, Ahmad HM. New fractional derivative expression of the shifted third kind Chebyshev polynomials: application to a type of nonlinear fractional pantograph differential equations. J Funct Spaces. 2022;2022:Article ID 3966135. DOI:10.1155/2022/3966135.
- Youssri YH, Abd-Elhameed WM, Mohamed AS, et al. Generalized lucas polynomial sequence treatment of fractional pantograph differential equation. Int J Appl Comput Math. 2021;7(27):1253–1268. DOI:10.1007/s40819-021-00958-y.
- Tunç C. Some stability and boundedness conditions for non-autonomous differential equations with deviating arguments. Electron J Qual Theory Differ Equ. 2010;1:12 pp. DOI:10.14232/ejqtde.2010.1.1.
- Tunç C. Stability and bounded of solutions to non-autonomous delay differential equations of third order. Nonlinear Dyn. 2010;62:945–953. https://doi.org/10.1007/s11071-010-9776-5.
- Tunç C. Stability to vector Liénard equation with constant deviating argument. Nonlinear Dyn. 2013;73:1245–1251. DOI:10.1007/s11071-012-0704-8.
- Tunc C. Qualitative properties in nonlinear Volterra integro-differential equations with delay. J Taibah Univ Sci. 2017;1(2):309–314. DOI:10.1016/j.jtusci.2015.12.009.
- Tunc C, Dinç Y. Qualitative properties of certain non-linear differential systems of second order. J Taibah Univ Sci. 2017;11(2):359–366. DOI:10.1016/j.jtusci.2016.05.002.
- Tunç C, Tunç O. On the asymptotic stability of solutions of stochastic differential delay equations of second order. J Taibah Univ Sci 2019;13(1):875–882. DOI:10.1080/16583655.2019.1595359.
- Tunç C, Tunç O. A note on certain qualitative properties of a second order linear differential system. Appl Math Inf Sci. 2015;9(2):953–956. DOI:10.12785/amis/090245.
- Tunç C, Tunç, O. On the boundedness and integration of non-oscillatory solutions of certain linear differential equations of second order. J Adv Res. 2016;7(1):165–168. DOI:10.1016/j.jare.2015.04.005.
- Tunç C, Tunç O. A note on the stability and boundedness of solutions to non-linear differential systems of second order. J Assoc Arab Univ Basic Appl Sci. 2017;24:169–175. DOI:10.1016/j.jaubas.2016.12.004.
- Lou Y, Ni WM. Diffusion, self-diffusion and cross-diffusion. J Differ Equ. 1996;131:79–131.