# American Institute of Mathematical Sciences

February  2021, 26(2): 1129-1148. doi: 10.3934/dcdsb.2020156

## Global analysis of a model of competition in the chemostat with internal inhibitor

 a. TAG_SUP Université Ibn Khaldoun, Tiaret 14000, Algérie b. TAG_SUP Ecole Normale Supérieure, Mostaganem 27000, Algérie c. TAG_SUP Université Djillali Liabès, LDM, Sidi Bel Abbès 22000, Algérie

* Corresponding author: Mohamed Dellal

Received  August 2019 Revised  January 2020 Published  February 2021 Early access  May 2020

Fund Project: The authors are supported by the "Directorate General for Scientific Research and Technological Development, Ministry of Higher Education and Scientific Research, Algeria" and the Euro-Mediterranean research network TREASURE (http://www.inra.fr/treasure)

A model of two microbial species in a chemostat competing for a single resource in the presence of an internal inhibitor is considered. The model is a four-dimensional system of ordinary differential equations. Using general growth rate functions of the species, we give a complete analysis for the existence and local stability of all steady states. We describe the behavior of the system with respect to the operating parameters represented by the dilution rate and the input concentrations of the substrate. The operating diagram has the operating parameters as its coordinates and the various regions defined in it correspond to qualitatively different asymptotic behavior: washout, competitive exclusion of one species, coexistence of the species, bistability, multiplicity of positive steady states. This bifurcation diagram which determines the effect of the operating parameters, is very useful to understand the model from both the mathematical and biological points of view, and is often constructed in the mathematical and biological literature.

Citation: Mohamed Dellal, Bachir Bar. Global analysis of a model of competition in the chemostat with internal inhibitor. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1129-1148. doi: 10.3934/dcdsb.2020156
##### References:

show all references

##### References:
Projection of equilibria $E_0$, $E_1$, $E_2$ and $E_c$ in the plane $(S,p)$ and existence and stability conditions of these points. Where stable equilibrium points are denoted by filled circles and unstable equilibrium points are denoted by empty circles
The existence of multiple positive equilibria when the condition (20) $\rm(a)$: is not satisfied; $\rm(b)$: is satisfied
Case $n = 3$: Local stability of $E_1$, $E_2$, $E_c^2$ and instability of $E_c^1$, $E_c^2$, $E_0$
are exactly the same except that the values of $\alpha_i$ have been inverted">Figure 5.  The operating diagram of the system (1) where $\mu_i$ are given by (2) and the curves $\Gamma_i^c$ do not intersect. $\rm(a)$: The occurrence of the bistability region $\mathcal{J}_7$, where the coexistence region $\mathcal{J}_6$ does not exist. $\rm(b)$: The occurrence of the coexistence region $\mathcal{J}_6$. The biological parameters used to construct Figs. Figs. 5(a, b) are exactly the same except that the values of $\alpha_i$ have been inverted
are exactly the same except that the values of $\alpha_i$ have been inverted">Figure 6.  The operating diagram of the system (1) where $\mu_i$ are given by (2) and the curves $\Gamma_i^c$ intersect. The pictures show the occurrence of the bistability region $\mathcal{J}_8$ of $E_2$ and $E_c^i$, The biological parameters used to construct Figs. Figs. 6(a, b) are exactly the same except that the values of $\alpha_i$ have been inverted
Illustrative graph of $\Gamma_i^c$, $i = 1,2$ defined by (33), showing the relative positions of the roots $S_c^1(\!D\!)$ and $S_c^2(\!D\!)$ of (39) with respect to the root $S_i(\!D\!,\!S^0\!)$ of (38). Region Ⅰ: $S_i<S_c^1<S_c^2$; region Ⅱ: $S_c^1<S_c^2<S_i$; region Ⅲ : $S_c^1<S_i<S_c^2$
. (a): Bistability of $E_1$ and $E_2$ when $(S^0,D) = (0.1,0.9)\in\mathcal{J}_7^a$. (b): Bistability of $E_c^1$ and $E_2$ when $(S^0,D) = (0.05,1.15)\in\mathcal{J}_8^a$. (c): Global stability of $E_c^1$ when $(S^0,D) = (0.02,1.15)\in\mathcal{J}_6^a$">Figure 7.  The trajectories of the reduced model (28) where $(S^0,D)$ are chosen in regions of Fig. 6(a). (a): Bistability of $E_1$ and $E_2$ when $(S^0,D) = (0.1,0.9)\in\mathcal{J}_7^a$. (b): Bistability of $E_c^1$ and $E_2$ when $(S^0,D) = (0.05,1.15)\in\mathcal{J}_8^a$. (c): Global stability of $E_c^1$ when $(S^0,D) = (0.02,1.15)\in\mathcal{J}_6^a$
Existence and local asymptotic stability of equilibria of system (1)
 Equilibria Existence Local exponential stability $E_0$ Always $\lambda_1>S^0$ and $\lambda_2>S^0$ $E_1$ $\lambda_1F_2(S_1)$ $E_2$ $\lambda_2F_1(S_2)$ $E_{c}$ (23) has a solution $(\alpha_1\gamma_1-\alpha_2\gamma_2)[F'_2(S_c)-F'_1(S_c)]>0$
 Equilibria Existence Local exponential stability $E_0$ Always $\lambda_1>S^0$ and $\lambda_2>S^0$ $E_1$ $\lambda_1F_2(S_1)$ $E_2$ $\lambda_2F_1(S_2)$ $E_{c}$ (23) has a solution $(\alpha_1\gamma_1-\alpha_2\gamma_2)[F'_2(S_c)-F'_1(S_c)]>0$
Parameter values used in Section 5 where $\mu_i$ are given by (2)
 Parameters $m_1$ $m_2$ $a_1$ $a_2$ $K_1$ $K_2$ $\alpha_1$ $\alpha_2$ Figures Units $h^{-1}$ $h^{-1}$ $gl^{-1}$ $gl^{-1}$ $gl^{-1}$ $gl^{-1}$ Case (a) 1.0 2.0 0.01 0.04 0.01 0.006 0.1 4.0 5(a) Case (b) 1.0 2.0 0.01 0.04 0.01 0.006 4.0 0.1 5(b) Case (c) 2.0 9.0 0.006 0.04 0.005 0.001 0.005 0.4 6(a), 7 Case (d) 2.0 9.0 0.006 0.04 0.005 0.001 0.4 0.005 6(b)
 Parameters $m_1$ $m_2$ $a_1$ $a_2$ $K_1$ $K_2$ $\alpha_1$ $\alpha_2$ Figures Units $h^{-1}$ $h^{-1}$ $gl^{-1}$ $gl^{-1}$ $gl^{-1}$ $gl^{-1}$ Case (a) 1.0 2.0 0.01 0.04 0.01 0.006 0.1 4.0 5(a) Case (b) 1.0 2.0 0.01 0.04 0.01 0.006 4.0 0.1 5(b) Case (c) 2.0 9.0 0.006 0.04 0.005 0.001 0.005 0.4 6(a), 7 Case (d) 2.0 9.0 0.006 0.04 0.005 0.001 0.4 0.005 6(b)
Existence and stability of equilibria in the regions of the operating diagrams of Fig. 5, when the curves $\Gamma_i^c$ do not intersect. The letter S (resp. U) means stable (resp. unstable) and empty if that equilibrium does not exist
 Region The relative positions of $S_i$ and $S_c^i$ $E_0$ $E_1$ $E_2$ $E_{c}^1$ $E_{c}^2$ $(S^0,D)\in\mathcal J_1$ $S_1$ and $S_2$ do not exist S $(S^0,D)\in\mathcal J_2$ $S_1$ does not exist U S $(S^0,D)\in\mathcal J_3$ $S_c^1  Region The relative positions of$ S_i $and$ S_c^i  E_0  E_1  E_2  E_{c}^1  E_{c}^2  (S^0,D)\in\mathcal J_1  S_1 $and$ S_2 $do not exist S$ (S^0,D)\in\mathcal J_2  S_1 $does not exist U S$ (S^0,D)\in\mathcal J_3  S_c^1
Existence and stability of equilibria in the regions of the operating diagrams of Fig. 6, when the curves $\Gamma_i^c$ intersect
 Region The relative positions of $S_i$ and $S_c^i$ $E_0$ $E_1$ $E_2$ $E_{c}^1$ $E_{c}^2$ $(S^0,D)\in\mathcal J_1$ $S_1$ and $S_2$ do not exist S $(S^0,D)\in\mathcal J_2$ $S_1$ does not exist U S $(S^0,D)\in\mathcal J_3$ $S_i  Region The relative positions of$ S_i $and$ S_c^i  E_0  E_1  E_2  E_{c}^1  E_{c}^2  (S^0,D)\in\mathcal J_1  S_1 $and$ S_2 $do not exist S$ (S^0,D)\in\mathcal J_2  S_1 $does not exist U S$ (S^0,D)\in\mathcal J_3  S_i
Parameter values of $S_i$ and $S_c^i$ used in Fig.7
 $(S^0,D)$ Regions $S_1$ $S_2$ $S_c^1$ $S_c^2$ Figures $(0.1,0.9)$ $\mathcal{J}_7^a$ 0.006 0.085 0.005 0.064 7(a) $(0.05,1.15)$ $\mathcal{J}_8^a$ 0.009 0.042 0.012 0.025 7(b) $(0.02,1.15)$ $\mathcal{J}_6^a$ 0.008 0.017 0.012 0.025 7(c)
 $(S^0,D)$ Regions $S_1$ $S_2$ $S_c^1$ $S_c^2$ Figures $(0.1,0.9)$ $\mathcal{J}_7^a$ 0.006 0.085 0.005 0.064 7(a) $(0.05,1.15)$ $\mathcal{J}_8^a$ 0.009 0.042 0.012 0.025 7(b) $(0.02,1.15)$ $\mathcal{J}_6^a$ 0.008 0.017 0.012 0.025 7(c)
 [1] Bachir Bar, Tewfik Sari. The operating diagram for a model of competition in a chemostat with an external lethal inhibitor. Discrete & Continuous Dynamical Systems - B, 2020, 25 (6) : 2093-2120. doi: 10.3934/dcdsb.2019203 [2] 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 [3] Georg Hetzer, Wenxian Shen. Two species competition with an inhibitor involved. Discrete & Continuous Dynamical Systems, 2005, 12 (1) : 39-57. doi: 10.3934/dcds.2005.12.39 [4] Jifa Jiang, Fensidi Tang. The complete classification on a model of two species competition with an inhibitor. Discrete & Continuous Dynamical Systems, 2008, 20 (3) : 659-672. doi: 10.3934/dcds.2008.20.659 [5] Hua Nie, Wenhao Xie, Jianhua Wu. Uniqueness of positive steady state solutions to the unstirred chemostat model with external inhibitor. Communications on Pure & Applied Analysis, 2013, 12 (3) : 1279-1297. doi: 10.3934/cpaa.2013.12.1279 [6] Nahla Abdellatif, Radhouane Fekih-Salem, Tewfik Sari. Competition for a single resource and coexistence of several species in the chemostat. Mathematical Biosciences & Engineering, 2016, 13 (4) : 631-652. doi: 10.3934/mbe.2016012 [7] Hua Nie, Yuan Lou, Jianhua Wu. Competition between two similar species in the unstirred chemostat. Discrete & Continuous Dynamical Systems - B, 2016, 21 (2) : 621-639. doi: 10.3934/dcdsb.2016.21.621 [8] Hua Nie, Feng-Bin Wang. Competition for one nutrient with recycling and allelopathy in an unstirred chemostat. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2129-2155. doi: 10.3934/dcdsb.2015.20.2129 [9] Mohamed Dellal, Bachir Bar, Mustapha Lakrib. A competition model in the chemostat with allelopathy and substrate inhibition. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021120 [10] Frédéric Mazenc, Michael Malisoff, Patrick D. Leenheer. On the stability of periodic solutions in the perturbed chemostat. Mathematical Biosciences & Engineering, 2007, 4 (2) : 319-338. doi: 10.3934/mbe.2007.4.319 [11] Hua Nie, Sze-bi Hsu, Jianhua Wu. A competition model with dynamically allocated toxin production in the unstirred chemostat. Communications on Pure & Applied Analysis, 2017, 16 (4) : 1373-1404. doi: 10.3934/cpaa.2017066 [12] Xiaoqing He, Sze-Bi Hsu, Feng-Bin Wang. A periodic-parabolic Droop model for two species competition in an unstirred chemostat. Discrete & Continuous Dynamical Systems, 2020, 40 (7) : 4427-4451. doi: 10.3934/dcds.2020185 [13] 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 [14] Frederic Mazenc, Gonzalo Robledo, Michael Malisoff. Stability and robustness analysis for a multispecies chemostat model with delays in the growth rates and uncertainties. Discrete & Continuous Dynamical Systems - B, 2018, 23 (4) : 1851-1872. doi: 10.3934/dcdsb.2018098 [15] Zhiqi Lu. Global stability for a chemostat-type model with delayed nutrient recycling. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 663-670. doi: 10.3934/dcdsb.2004.4.663 [16] E. Cabral Balreira, Saber Elaydi, Rafael Luís. Local stability implies global stability for the planar Ricker competition model. Discrete & Continuous Dynamical Systems - B, 2014, 19 (2) : 323-351. doi: 10.3934/dcdsb.2014.19.323 [17] Pan Zheng. Asymptotic stability in a chemotaxis-competition system with indirect signal production. Discrete & Continuous Dynamical Systems, 2021, 41 (3) : 1207-1223. doi: 10.3934/dcds.2020315 [18] Shanshan Chen, Jianshe Yu. Stability and bifurcation on predator-prey systems with nonlocal prey competition. Discrete & Continuous Dynamical Systems, 2018, 38 (1) : 43-62. doi: 10.3934/dcds.2018002 [19] Satoshi Kosugi, Yoshihisa Morita, Shoji Yotsutani. A complete bifurcation diagram of the Ginzburg-Landau equation with periodic boundary conditions. Communications on Pure & Applied Analysis, 2005, 4 (3) : 665-682. doi: 10.3934/cpaa.2005.4.665 [20] Masaaki Mizukami. Boundedness and asymptotic stability in a two-species chemotaxis-competition model with signal-dependent sensitivity. Discrete & Continuous Dynamical Systems - B, 2017, 22 (6) : 2301-2319. doi: 10.3934/dcdsb.2017097

2020 Impact Factor: 1.327