Stochastic Systems and Control: Theory and ApplicationsView this Special Issue
Research Article | Open Access
Threshold Dynamics of a Stochastic Chemostat Model with Two Nutrients and One Microorganism
A new stochastic chemostat model with two substitutable nutrients and one microorganism is proposed and investigated. Firstly, for the corresponding deterministic model, the threshold for extinction and permanence of the microorganism is obtained by analyzing the stability of the equilibria. Then, for the stochastic model, the threshold of the stochastic chemostat for extinction and permanence of the microorganism is explored. Difference of the threshold of the deterministic model and the stochastic model shows that a large stochastic disturbance can affect the persistence of the microorganism and is harmful to the cultivation of the microorganism. To illustrate this phenomenon, we give some computer simulations with different intensity of stochastic noise disturbance.
Chemostat is commonly used to describe the dynamics of a microbial population in a continuous bioreactor in which microorganisms grow on a substrate and has attracted great interest of many scholars [1–8], since it was first introduced by Monod . A single simple species chemostat model with Michaelis-Menten-Monod functional response was proposed by  as follows:where is the concentration of the nutrient, is the concentration of the organism, is the dilution (or washout) rate, is the maximal growth rate, is the Michaelis-Menten (or half-saturation) constant with units of concentration, and is a “yield” constant reflecting the conversion of nutrient to organism.
However, experimental results have indicated that the microorganisms depend on a variety of nutrition substances such as carbon, nitrogen, energy, growth factors, inorganic salts, and water. Then the model of microorganisms species growth in the chemostat on two nutrients is considered by [10–14]. A model of single-species growth in the chemostat on two substitutable resources with Michaelis-Menten-Monod functional response was proposed by  as follows:
However, it is now well known that stochastic noise is widely present in biological systems and so on [15–33] and microorganisms are inevitably influenced by some random factors in the process of cultivation. To better understand the dynamic behavior of the chemostat, a host of scholars proposed a slice of stochastic chemostat models and studied the effect of the random noise on the dynamic behavior of the stochastic models. As an example, Imhof and Walcher  proposed a stochastic chemostat model for a single microorganism species consuming a single nutrient. They found that random effects may lead to extinction in scenarios where the deterministic model predicts persistence. Recently, Xu and Yuan  established a stochastic chemostat model in which the maximal growth rate is influenced by the white noise in environment as follows:They got an analogue break-even concentration involving the white noise which can determine the exclusion and persistence of the microorganism. And more stochastic chemostat models can be found in [36–39].
Motivated by the papers mentioned above, in this paper, we further consider a model of single-species growth in the chemostat on two supplementary resources with Michaelis-Menten-Monod functional response and environmental noise. We assume that the maximal growth rate is perturbed by white noises so that where is a standard Brownian motion with intensity Then the resultant model takes the following form:Our main objective in the rest of this paper is to investigate the threshold dynamics of stochastic chemostat model (5) and explore the conditions under which microorganisms will die out or exist.
In this section, we will give some notations, definitions, and lemmas which will be used for analyzing our main results. To this end, throughout this paper, we let be a complete probability space with a filtration satisfying the usual conditions: it is increasing and right continuous while contains all -null sets; we use to represent a scalar Brownian motion defined on the complete probability space ; also let . If for an integrable function on , define Then we have the following.
Definition 1. For system (5),(i)the microorganism is said to be extinctive if ,(ii)the microorganism is said to be permanent in mean if there exists a positive constant such that .
Then, one can show the following lemmas.
Lemma 3. Letting be a solution of system (5) with initial value , then
3. Dynamics of Deterministic System (2)
In this section, we will focus on the deterministic system (2). It is easy to see that the equilibria point of (2) satisfyand, obviously, model (2) has a microorganism extinction equilibrium Let be the coexistence equilibrium of model (2), which satisfieswhere Then we have that Denotewhere ,
Obviously, If , we have If , we have Thus, equation has one positive root at least, and
From the third equation of (12), we have Then we have the following theorem.
Theorem 4. If and , then system (2) has unique positive equilibrium
Regarding the stability of these equilibria, we have the following theorem.
Theorem 5. Then for system (2), one has the following. (i)If , microorganism extinction equilibrium is locally stable; if it is unstable.(ii)If and , the coexistence equilibrium is locally stable.
Proof. Linearizing the system at the equilibrium gives the Jacobianwhere The characteristic equation giveswhere Obviously, we have, at , Then we have and thus if , all the eigenvalues of (23) have negative real part; then, by the stability theory, is stable.
And, at , we have here is used. Then all the eigenvalues of (23) have negative real part; thus, by the stability theory, the diseases equilibrium is stable as long as it exists.
4. Dynamics of Stochastic System (5)
In this section, we explore the conditions leading to the extinction of the two infectious diseases. Denote where is introduced in (16). Then we have the following.
Theorem 6. For system (5), if one of the following holds, (i), , and ,(ii), , and ,(iii), , and ,(iv), , and ,then the microorganism of system (5) goes to extinction almost surely. Moreover, almost surely.
Proof. Let be a solution of system (5) with initial value . Applying Itô’s formula to system (5) results inwhere ,
Integrating both sides of (30) from to giveswhere known as the local continuous martingale, and . Obviously, we need to estimate the maximum value of
Let us consider quadratic functionIt is easy to verify that when , reaches its maximum value at ; and when , achieve its maximum value at Then, in (31), we have four cases to be discussed, depending on whether or , which are as follows: Case 1: , ; Case 2: , ; Case 3: , ; and Case 4: ,
For Case 1, since , , then achieve the maximum value Then we can easily see from (31) thatDividing both sides of (34) by , we haveand, by Lemma 3, we have Then, taking the limit superior on both sides of (35) leads to which implies , and here is used.
Case 2. ,
In this case, we can easily see from (31) thatDividing both sides of (38) by , we haveand, by Lemma 3, we have Then, taking the limit superior on both sides of (38) leads to which implies
The same discussion can be used in Case 3; here we omit it.
Next, we consider Case 4: , From (31), we haveDividing both sides of (42) by , we haveand, by Lemma 3, we have Then, taking the limit superior on both sides of (43) leads to which implies
Next, we prove the last conclusion. Given , since , we have for large enough. By the first equation of system (5), we have Then when we haveOn the other hand from the proof of Lemma 2, we have Let . Then one hasFrom (47) and (49), we have almost surely.
By employing the method similar above, it then follows that almost surely. This completes the proof of Theorem 6.
4.2. Permanence in Mean
Theorem 7. If , then the microorganism is permanent in mean; moreover, satisfies where
Proof. Integrating from to and dividing by on both sides of system (5) yieldThen one can getApplying Itô’s formula givesIntegrating from 0 to and dividing by on both sides of (55) yields where , Noticing that then we haveIf , we can get By inequality (59), we haveBy Lemma 3, we get that According to Lemma 2, one sees that , , and , and then one has and Thus taking the inferior limit of both sides of (60) yieldsAnd if , we can get where ,
Inequality (62) can be rewritten asTaking the inferior limit of both sides of (63) yieldsLet , and we get from (61) and (64)This completes the proof of Theorem 7.
Remark 8. Theorems 6 and 7 show that the condition for the microorganism to go to extinction or permanence depends on the intensity of the noise disturbances completely. And small noise disturbances will be beneficial to the cultivation of the microorganism; conversely, large white noise disturbance is harmful to the cultivation of the microorganism.
5. Conclusion and Numerical Simulation
This paper proposes and investigates a new stochastic chemostat model with two substitutable nutrients and one microorganism. Then main objective in this paper is to investigate the threshold dynamics of stochastic chemostat model (5) and explore the conditions which can determine the extinction and permanence of the microorganism using two substitutable nutrients. Firstly, for the corresponding deterministic model, the threshold for extinction or existence of the microorganism is obtained by analyzing the stability of the equilibria. Then the threshold of the stochastic chemostat for the extinction and the permanence in mean of the microorganism is explored. The results show that there exists a significant difference between the threshold of the deterministic system and the stochastic system, which makes the persistent microorganism of a deterministic system become extinct due to large stochastic disturbance. That is, large stochastic disturbance is harmful to the cultivation of the microorganism. It is worth mentioning that this paper is a promotion of the work of Xu and Yuan .
Next, using the Euler Maruyama (EM) method , we give some numerical simulation to illustrate the extinction and persistence of the microorganism in stochastic system and corresponding deterministic system for comparison.
Firstly, we begin from a deterministic system; the basic parameters are set as , , , , , , and Direct calculation shows that , , and Then according to Theorems 4 and 5, the deterministic system has a unique stable positive equilibrium , which is locally stable and the deterministic system is permanent (see Figure 1).
Next, we consider the influence of stochastic disturbance on the above deterministic system. According to Theorem 6, different parameters are chosen to give insights into the reasonability of the results stated in Theorem 6.
We choose different value of parameters and and discuss below five different cases.
(d) , ,
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work is supported by Shandong Provincial Natural Science Foundation (no. ZR2015AQ001), the National Natural Science Foundation of China (no. 11371230), and Research Funds for Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources by Shandong Province and SDUST (2014TDJH102).
- D. Herbert, R. Elsworth, and R. C. Telling, “The continuous culture of bacteria; a theoretical and experimental study,” Journal of General Microbiology, vol. 14, no. 3, pp. 601–622, 1956.
- H. L. Smith and P. Waltman, The Theory of the Chemostat: Dynamics of Microbial Competition, Cambridge Studies in Mathematical Biology, Cambridge University Press, Cambridge, UK, 1995.
- T. Zhang, T. Zhang, and X. Meng, “Stability analysis of a chemostat model with maintenance energy,” Applied Mathematics Letters. An International Journal of Rapid Publication, vol. 68, pp. 1–7, 2017.
- X. Meng, L. Wang, and T. Zhang, “Global dynamics analysis of a nonlinear impulsive stochastic chemostat system in a polluted environment,” Journal of Applied Analysis and Computation, vol. 6, no. 3, pp. 865–875, 2016.
- W. Wang, W. Ma, and H. Yan, “Global dynamics of modeling flocculation of microorganism,” Applied Sciences (Switzerland), vol. 6, no. 8, article no. 221, 2016.
- X. Meng, Q. Gao, and Z. Li, “The effects of delayed growth response on the dynamic behaviors of the Monod type chemostat model with impulsive input nutrient concentration,” Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal, vol. 11, no. 5, pp. 4476–4486, 2010.
- J. Gao, B. Shen, E. Feng, and Z. Xiu, “Modelling and optimal control for an impulsive dynamical system in microbial fed-batch culture,” Computational and Applied Mathematics, vol. 32, no. 2, pp. 275–290, 2013.
- T. Zhang, W. Ma, and X. Meng, “Global dynamics of a delayed chemostat model with harvest by impulsive flocculant input,” Advances in Difference Equations, Paper No. 115, 17 pages, 2017.