# American Institute of Mathematical Sciences

July  2019, 24(7): 3439-3451. doi: 10.3934/dcdsb.2018327

## Complex dynamics in a discrete-time size-structured chemostat model with inhibitory kinetics

 a. School of Mathematics and Computational Science, Xiangtan University, Xiangtan, 411105, Hunan, China b. School of Computer Science and Network Security, Dongguan University of Technology, Dongguan, 523808, Guangdong, China c. College of Finance and Statistics, Hunan University, Changsha, 410079, Hunan, China d. Department of Mathematics and Statistics, University of New Brunswick, Fredericton, NB, Canada

* Corresponding author: Lin Wang

Received  May 2018 Revised  August 2018 Published  January 2019

Fund Project: The work of DZ was partially supported by the National Natural Science Foundation of China (No. 11501193) and the China Post Doctorial Fund (No. 2015M582335). LW was partially supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (NSERC)

An inhibitory uptake function is incorporated into the discrete, size-structured nonlinear chemostat model developed by Arino et al. (Journal of Mathematical Biology, 45(2002)). Different from the model with a monotonically increasing uptake function, we show that the inhibitory kinetics can induce very complex dynamics including stable equilibria, cycles and chaos (via the period-doubling cascade). In particular, when the nutrient concentration in the input feed to the chemostat $S^0$ is larger than the upper break-even concentration value $\mu$, the model exhibits three types of bistability allowing a stable equilibrium to coexist with another stable equilibrium, or a stable cycle or a chaotic attractor.

Citation: Dan Zhang, Xiaochun Cai, Lin Wang. Complex dynamics in a discrete-time size-structured chemostat model with inhibitory kinetics. Discrete & Continuous Dynamical Systems - B, 2019, 24 (7) : 3439-3451. doi: 10.3934/dcdsb.2018327
##### References:

show all references

##### References:
Bifurcation diagram of the limiting system (7). Here $f(s) = \frac{as}{1+0.1s+0.01s^2}, E = 0.1, S^{0} = 100, U_{0} = 70$ with $a\in[0.049, 0.059]$ Thus $aS^{0}>1$ and $S^{0}> \mu$. A cascade of period-doublings to chaos occurs as $a$ increases
Left: a numerical solution of the limiting system (7) with $U_0 = 60$; Right: a numerical solution of system (11) with $(U_0, S_0) = (60,6)$. Here $f(s) = \frac{0.05s}{1+0.1s+0.01s^2}$, $E = 0.1$, $S^{0} = 100$
Numerical solutions of system (11) with $f(s) = \frac{0.05s}{1+0.1s+0.01s^2}$, $E = 0.1$, $S^{0} = 100$. Left: $(U_t,S_t)\to E_0 = (0, S^0)$ as $t\to\infty$, initial condition $(U_0, S_0) = (10,6)$ was used; Right: $(U_t, S_t) \to E_1 = (S^0-\lambda, \lambda)$ as $t\to\infty$, initial condition was $(U_0, S_0) = (80,6)$
Numerical solutions of system (11) with $f(s) = \frac{0.054s}{1+0.1s+0.01s^2}$, $E = 0.1$, $S^{0} = 100$. Left: $(U_t,S_t)\to (0, S^0)$ as $t\to\infty$, initial condition $(U_0, S_0) = (10,6)$ was used; Right: $(U_t, S_t)$ approaches a stable $2-cycle$, initial condition was $(U_0, S_0) = (80,6)$
Numerical solutions of system (11) with $f(s) = \frac{0.059s}{1+0.1s+0.01s^2}$, $E = 0.1$, $S^{0} = 100$. Left: $(U_t,S_t)\to (0, S^0)$ as $t\to\infty$, initial condition $(U_0, S_0) = (10,6)$ was used; Right: $(U_t, S_t)$ approaches a chaotic attractor, initial condition was $(U_0, S_0) = (80,6)$
 [1] H. L. Smith, X. Q. Zhao. Competitive exclusion in a discrete-time, size-structured chemostat model. Discrete & Continuous Dynamical Systems - B, 2001, 1 (2) : 183-191. doi: 10.3934/dcdsb.2001.1.183 [2] Josselin Garnier, George Papanicolaou, Tzu-Wei Yang. Mean field model for collective motion bistability. Discrete & Continuous Dynamical Systems - B, 2019, 24 (2) : 851-879. doi: 10.3934/dcdsb.2018210 [3] Nabil T. Fadai, Michael J. Ward, Juncheng Wei. A time-delay in the activator kinetics enhances the stability of a spike solution to the gierer-meinhardt model. Discrete & Continuous Dynamical Systems - B, 2018, 23 (4) : 1431-1458. doi: 10.3934/dcdsb.2018158 [4] Sze-Bi Hsu, Cheng-Che Li. A discrete-delayed model with plasmid-bearing, plasmid-free competition in a chemostat. Discrete & Continuous Dynamical Systems - B, 2005, 5 (3) : 699-718. doi: 10.3934/dcdsb.2005.5.699 [5] Zhipeng Qiu, Jun Yu, Yun Zou. The asymptotic behavior of a chemostat model. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 721-727. doi: 10.3934/dcdsb.2004.4.721 [6] Eric A. Carlen, Maria C. Carvalho, Jonathan Le Roux, Michael Loss, Cédric Villani. Entropy and chaos in the Kac model. Kinetic & Related Models, 2010, 3 (1) : 85-122. doi: 10.3934/krm.2010.3.85 [7] Kaijen Cheng, Kenneth Palmer. Chaos in a model for masting. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 1917-1932. doi: 10.3934/dcdsb.2015.20.1917 [8] Jianquan Li, Zuren Feng, Juan Zhang, Jie Lou. A competition model of the chemostat with an external inhibitor. Mathematical Biosciences & Engineering, 2006, 3 (1) : 111-123. doi: 10.3934/mbe.2006.3.111 [9] Carole Guillevin, Rémy Guillevin, Alain Miranville, Angélique Perrillat-Mercerot. Analysis of a mathematical model for brain lactate kinetics. Mathematical Biosciences & Engineering, 2018, 15 (5) : 1225-1242. doi: 10.3934/mbe.2018056 [10] María J. Cáceres, Ricarda Schneider. Blow-up, steady states and long time behaviour of excitatory-inhibitory nonlinear neuron models. Kinetic & Related Models, 2017, 10 (3) : 587-612. doi: 10.3934/krm.2017024 [11] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [12] Masaki Sekiguchi, Emiko Ishiwata, Yukihiko Nakata. Dynamics of an ultra-discrete SIR epidemic model with time delay. Mathematical Biosciences & Engineering, 2018, 15 (3) : 653-666. doi: 10.3934/mbe.2018029 [13] Eduardo Liz. A new flexible discrete-time model for stable populations. Discrete & Continuous Dynamical Systems - B, 2018, 23 (6) : 2487-2498. doi: 10.3934/dcdsb.2018066 [14] Edoardo Beretta, Fortunata Solimano, Yanbin Tang. Analysis of a chemostat model for bacteria and virulent bacteriophage. Discrete & Continuous Dynamical Systems - B, 2002, 2 (4) : 495-520. doi: 10.3934/dcdsb.2002.2.495 [15] Yangjin Kim, Khalid Boushaba. An enzyme kinetics model of tumor dormancy, regulation of secondary metastases. Discrete & Continuous Dynamical Systems - S, 2011, 4 (6) : 1465-1498. doi: 10.3934/dcdss.2011.4.1465 [16] Jose S. Cánovas, Tönu Puu, Manuel Ruiz Marín. Detecting chaos in a duopoly model via symbolic dynamics. Discrete & Continuous Dynamical Systems - B, 2010, 13 (2) : 269-278. doi: 10.3934/dcdsb.2010.13.269 [17] Wenbo Cheng, Wanbiao Ma, Songbai Guo. A class of virus dynamic model with inhibitory effect on the growth of uninfected T cells caused by infected T cells and its stability analysis. Communications on Pure & Applied Analysis, 2016, 15 (3) : 795-806. doi: 10.3934/cpaa.2016.15.795 [18] Sie Long Kek, Mohd Ismail Abd Aziz, Kok Lay Teo, Rohanin Ahmad. An iterative algorithm based on model-reality differences for discrete-time nonlinear stochastic optimal control problems. Numerical Algebra, Control & Optimization, 2013, 3 (1) : 109-125. doi: 10.3934/naco.2013.3.109 [19] Rui Xu, M.A.J. Chaplain, F.A. Davidson. Periodic solutions of a discrete nonautonomous Lotka-Volterra predator-prey model with time delays. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 823-831. doi: 10.3934/dcdsb.2004.4.823 [20] Yun Kang. Permanence of a general discrete-time two-species-interaction model with nonlinear per-capita growth rates. Discrete & Continuous Dynamical Systems - B, 2013, 18 (8) : 2123-2142. doi: 10.3934/dcdsb.2013.18.2123

2018 Impact Factor: 1.008