Table of Contents Author Guidelines Submit a Manuscript
Volume 2019, Article ID 4681205, 18 pages
Research Article

Dynamical Analysis of a Stochastic Multispecies Turbidostat Model

1Department of Mathematics, Hubei University for Nationalities, Enshi, Hubei 445000, China
2Department of Mathematics, City University of Hong Kong, Hong Kong, China
3School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China

Correspondence should be addressed to Zuxiong Li; moc.621@7240xzil

Received 31 July 2018; Accepted 4 December 2018; Published 10 January 2019

Academic Editor: Marcelo Messias

Copyright © 2019 Yu Mu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A stochastic turbidostat system in which the dilution rate is subject to white noise is investigated in this paper. First of all, sufficient conditions of the competitive exclusion among microorganisms are obtained by employing the techniques of stochastic analysis. Furthermore, the results demonstrate that the competition among microorganisms and stochastic disturbance will affect the dynamical behaviors of microorganisms. Finally, the theoretical results obtained in this contribution are illustrated by numerical simulations.

1. Introduction

The chemostat and turbidostat, two types of devices for continuous cultivation of microorganism, have been utilized to analyze population dynamics. Novick et al. [1] first proposed the mechanism of chemostat in 1950, and then an increasing number of researchers have devoted themselves to investigate chemostat systems [27]. However, there exist some drawbacks in chemostat model with a constant dilution rate, such as the waste of substrate and higher viscosity caused by the mass transfer efficiency. Models with the dilution rate related to the state of microorganism, which can be called turbidostat (see [8]), can overcome the above drawbacks and be helpful to design optimal strategies to improve the substrate utilization. Flegr [9] investigated turbidostat system with two species from numerical simulation and De Leenheer et al. [4] analyzed the system theoretically. Li [10] systematically investigated the turbidostat model and established sufficient conditions of coexistence of two species. Subsequently, many researchers investigated turbidostat and chemostat systems. Many valuable and interesting results were obtained [1116].

According to May [17], the parameters of the system, such as birth rate, death rate, and the input concentration of nutrient, are inevitably disturbed by environmental factors. In recent years, more and more stochastic ecological models were chosen to describe the dynamics of populations. Grasmanet et al. [18] established a stochastic chemostat model with three trophic levels. Using singular perturbation methods, they obtained the expected breakdown time and analyzed the influence of stochastic factors on dynamic behaviors of the system in detail. Considering the stochastic factors in chemostat, Campillo et al. [19] proposed stochastic systems with different population scales, and they investigated these systems and derived the field of validity for different population scales. Collet et al. [5] pointed out that the population size in the device was determined not only by the concentration of the nutrient but also by the stochastic birth and death rates. By analyzing the long time behaviors of microbes, they derived the conditions for the global existence of the solutions of the system and the existence of quasi-stationary distribution. Meng et al. [6] constructed an impulsive stochastic chemostat model and investigated the extinction and permanence of the microorganisms. They also pointed out that small perturbation from white noise could cause the extinction of microorganisms. Xu et al. [7] investigated a chemostat model containing telegraph noise with the help of Markov chain and derived the break-even concentration which determines the persistence or extinction of microorganisms. One can also find a large number of stochastic models on cultivation of microorganisms [14, 15, 2023] as well as in other fields [2428].

Competition is all around environment because of the limitation of natural resource. Ghoul et al. [29] and Kaye et al. [30] pointed out that microbes expressed many competitive behaviors with their neighbours or plants for scarce nutrients and limited spaces. Many practical results were obtained by investigating the corresponding competitive models. Butler et al. [16] established a system containing two competitors and a growth-limiting nutrient. They derived conditions for uniform persistence and local stability of the food web and also obtained the conditions for predator-mediated coexistence. Wolkowicz et al. [31] proposed a competitive chemostat model with distributed delay to describe the process of nutrient consumption. They proved that there existed only one survivor under any conditions and also pointed out that the theoretical results in their paper were valid for all systems with monotone growth response functions. Liu et al. [32] established a stochastic competitive model and obtained sufficient conditions of extinction and persistence (including weak persistence and strong persistence). They stated that only one species survived under certain stochastic noise perturbation. Xu et al. [33] analyzed a competitive chemostat model and derived the critical value of the noise. Zhao et al. [34] discussed on a stochastic competition system and derived sufficient conditions of persistence and stationary distribution for each population. Actually, a great number of stochastic competition models have been studied [3540].

In this paper, considering the dilution rate of microorganisms related to the feedback control and the influence of stochastic factors from the environment, we establish the following stochastic turbidostat system based on Zhang et al. [21] who constructed a stochastic chemostat model and obtained the conditions of competitive exclusion.where and represent the concentration of nutrient and microorganisms at the time , separately. expresses the input concentration of nutrition, is Brownian motion defined on a complete probability space () and stands for the density of the white noise, ( and are constants) denotes the dilution rate of the turbidostat and ( is constant) is the sum of the dilution rate and death rate of the microorganisms in the system. is the Holling 2 functional response function ( and are constants).

This paper is organized as follows. In Section 2, the existence and uniqueness of the positive solution of system (1) are verified. Section 3 demonstrates the main results. A discussion is given and a numerical example is offered to verify our theoretical results in Section 4.

2. Existence and Uniqueness of the Positive Solution

We in this section demonstrate the existence and uniqueness of global positive solution of system (1).

Theorem 1. If , then stochastic system (1) has a unique positive solution almost surely with initial value .

Proof. For , letAccording to the It formula, we havewith initial value . It is obvious that and are continuous and differentiable functions. Hence (3) satisfies local Lipschitz condition, which means that there exists a unique local solution for . Then is the unique positive local solution for system (1) with the initial value . Next, we prove that the positive solution of system (1) is global; that is, . We first define the filtration as and choose an enough large such that . Then for all , we define a stopping time byNext, we further show . According to the definition of the stopping time, Therefore, is increasing as . Set . Then . Now, it is sufficient for us to show . Arguing by contradiction, we assume that . Then there exist constants and such that . Therefore, there is a constant such that for all .
Define a function bywhere is a positive constant. It is obvious that . Hence . By the It formula, the following equation can be derivedwhere Through calculations, the following inequality can be obtained Defining , we have which yields Integrating from to and taking expectation in both sides of inequality (11), one can obtainLet for all , then . Based on the definition of the stopping time, we have or for all . Hence there is a positive constant such that . Furthermore, Setting , we have , which is contradictory with . Hence . This completes the proof.

Remark 2. If we take as a threshold value, then the microorganisms in system (1) will survive as and be extinct as .

3. The Principle of Competitive Exclusion

In this section, we prove that stochastic turbidostat system (1) satisfies the principle of competitive exclusion. For simplicity, we first define

Lemma 3. The solution of system (1) satisfies with any initial value . Here, is a positive constant satisfying

Proof. Let and where is a nonnegative constant decided later. By the It formula, we derive We can further obtain Denote where Let be a constant. Then For , we have where Define a function by and which gives Integrating from to and taking expectation for (25), one has Then it is easy to obtain which leads to is a continuous function, so there is a constant such that According to (24), for small enough , , integrating from to and taking expectation, we can get where here, Applying Burkholder-Davis-Gundy inequality [41], we have Hence, we have Particularly, one can choose such that which yields For any , the following inequality can be obtained by Chebyshev’s inequality [41], Using Borel-Cantelli Lemma [41], for all and sufficiently large , we have Therefore, there is a , for all and ; the above equation holds. In addition, for all , when and , which means When , where is determined by (22). Therefore, This completes the proof.

Lemma 4. If , then the solution of system (1) with initial value satisfies where is the break-even concentration of and is the equilibrium concentration of .

Proof. Define It follows from , Lemma 3, and Burkholder-Davis-Gundy inequality [41] that By Doob’s martingale inequality [41], for arbitrary , it follows that By Borel-Cantelli Lemma [41], for all and sufficient large , we derive Thus, there exists a positive constant such that for all and , which implies Taking superior limit leads to Sending , we have Then there exist a constant and a set for small enough such that . Thus for , we have which means Together with we derive and Applying the same methods, we can prove that This completes the proof.

Definition 5. Stochastic equation has a solution weakly converging to the distribution . Here is a probability measure of such that . Particularly, has density , where is a normal constant and

Theorem 6. is the solution of system (1) with any initial value . Moreover, if then we have where , , is the break-even concentration of , is determined later by (74), is the boundary equilibrium of the corresponding deterministic model, and is defined in Definition 5.

Proof. The deterministic system corresponding to (1) has an equilibrium under the condition of , where and satisfy Define function by Then by the It formula, we get where According to (66), we have It follows from some calculations that Set</