
Citation: Tomás Caraballo, Renato Colucci, Javier López-de-la-Cruz, Alain Rapaport. Study of the chemostat model with non-monotonic growth under random disturbances on the removal rate[J]. Mathematical Biosciences and Engineering, 2020, 17(6): 7480-7501. doi: 10.3934/mbe.2020382
[1] | Tewfik Sari, Frederic Mazenc . Global dynamics of the chemostat with different removal rates and variable yields. Mathematical Biosciences and Engineering, 2011, 8(3): 827-840. doi: 10.3934/mbe.2011.8.827 |
[2] | Marion Weedermann . Analysis of a model for the effects of an external toxin on anaerobic digestion. Mathematical Biosciences and Engineering, 2012, 9(2): 445-459. doi: 10.3934/mbe.2012.9.445 |
[3] | Alain Rapaport, Jérôme Harmand . Biological control of the chemostat with nonmonotonic response and different removal rates. Mathematical Biosciences and Engineering, 2008, 5(3): 539-547. doi: 10.3934/mbe.2008.5.539 |
[4] | Alejandro Rincón, Fredy E. Hoyos, Gloria Restrepo . Global stability of a continuous bioreactor model under persistent variation of the dilution rate. Mathematical Biosciences and Engineering, 2023, 20(2): 3396-3424. doi: 10.3934/mbe.2023160 |
[5] | Zhihao Ke, Chaoqun Xu . Structure analysis of the attracting sets for plankton models driven by bounded noises. Mathematical Biosciences and Engineering, 2023, 20(4): 6400-6421. doi: 10.3934/mbe.2023277 |
[6] | Virginia Giorno, Serena Spina . On the return process with refractoriness for a non-homogeneous Ornstein-Uhlenbeck neuronal model. Mathematical Biosciences and Engineering, 2014, 11(2): 285-302. doi: 10.3934/mbe.2014.11.285 |
[7] | Ihab Haidar, Alain Rapaport, Frédéric Gérard . Effects of spatial structure and diffusion on the performances of the chemostat. Mathematical Biosciences and Engineering, 2011, 8(4): 953-971. doi: 10.3934/mbe.2011.8.953 |
[8] | Meng Gao, Xiaohui Ai . A stochastic Gilpin-Ayala mutualism model driven by mean-reverting OU process with Lévy jumps. Mathematical Biosciences and Engineering, 2024, 21(3): 4117-4141. doi: 10.3934/mbe.2024182 |
[9] | Jean-Jacques Kengwoung-Keumo . Dynamics of two phytoplankton populations under predation. Mathematical Biosciences and Engineering, 2014, 11(6): 1319-1336. doi: 10.3934/mbe.2014.11.1319 |
[10] | Gonzalo Robledo . Feedback stabilization for a chemostat with delayed output. Mathematical Biosciences and Engineering, 2009, 6(3): 629-647. doi: 10.3934/mbe.2009.6.629 |
The chemostat refers to a laboratory device used for the growth of micro-organisms in a culture environment [1,2], that has been regarded as an idealization of the nature to study microbial ecosystems in stationary stage [3]. It turned out to be an important investigation field due to a large number of applications, especially in waste water treatment [4,5] but also in ecological and environmental sciences (see [6,7,8,9,10]).
It is worth mentioning that the chemostat has been subject to a large number of scientific publications and books, not only in biology and ecology but also in mathematics. Indeed there exists a specific research area about the so-called "theory of the chemostat" [11,12] where many researchers have been involved in the last years. This interest has been strengthened by the fact that the chemostat device can be mathematically modeled in a simple way which reproduces quite faithfully real bio-processes.
Let us recall quickly in what consists the original chemostat device. It is composed of three tanks, the feed bottle, the culture vessel and the collection vessel, which are interconnected by pumps (see Figure 1). The substrate is stored in the feed bottle and provided with a given flow rate to the culture vessel, where interactions between nutrient and microbial biomass take place. The media from the culture vessel is also withdrawn towards the collection vessel with the same flow rate, to keep the volume in the culture vessel constant.
The mathematical model of the chemostat is given by the following dynamical system
dsdt=D(sin−s)−μ(s)x, | (1.1) |
dxdt=−Dx+μ(s)x, | (1.2) |
where s=s(t) and x=x(t) denote the concentration of substrate and species, respectively, sin is the input concentration of nutrient, D is the removal rate (also called input flow rate), and μ(⋅) is the specific growth function describing the kinetics of the nutrient consumption by the bacterial species. Here, we assume that the yield coefficient of the conversion of the substrate into biomass is equal to 1 (that is always possible to impose by a change of the unit of the biomass concentration). More specifically, we consider in the present work the Haldane growth function
μ(s)=μ0sks+s+s2ki,s≥0, | (1.3) |
where ks is the affinity constant and ki a parameter modeling the growth inhibition under large concentrations of substrate. Since it will be useful later, let us define the number
sm:=argmaxμ(s)=√kiks. | (1.4) |
Many works have been dedicated to this classical deterministic model of the chemostat (see for instance [11,12]) but most of time the removal rate D is kept constant, although it is well-known that in practice it is frequently subject to disturbances (see for instance [13] where some chronicles of time varying removal rates are depicted).
Motivated by this fact, it has been proposed to model the perturbations on the input flow rate in the chemostat model (1.1)–(1.2) by D+Φ(z∗(θtω)), where z∗(θtω) denotes the Ornstein-Uhlenbeck process (introduced in more detail in Section 2) and Φ is a bounded function defined as
Φ(z)=2dπarctan(z), | (1.5) |
where d>0 (see, for instance, [14,15] for other possible functions). Indeed, any odd function Φ which is increasing and such that limz→+∞Φ(z)=d can be considered.
In this way, once practitioners provide an interval [Dl,Dr]⊂R, 0<Dl<D<Dr<∞, (typically obtained from observations) we can define d=Dr−D=D−Dl and then the perturbed input flow is bounded for every time and any realization of the noise, i.e.,
Dl≤D+Φ(z∗(θtω))≤Dr∀t∈R. | (1.6) |
Then, the random chemostat model becomes
dsdt=(D+Φ(z∗(θtω)))(sin−s)−μ(s)x, | (1.7) |
dxdt=(D+Φ(z∗(θtω)))x+μ(s)x. | (1.8) |
One may wonder the reason why we consider this way of modeling bounded random fluctuations to perturb the input flow D in the model (1.1)–(1.2), instead of considering other stochastic process, such as for instance the well-known standard Wiener process. Indeed, this way has been typically used to model real perturbations and provides several advantages from both mathematical and biological points of view, see for instance [15].
On the one hand, this way of modeling noise fits in a loyal way the bounded variations of the input flow rate observed in real life. On the contrary, the Wiener process is unbounded with probability one which leads to arbitrary large (possibly negative) values of the corresponding perturbed input flow rate, which is not realistic at all from the biological point of view. We refer readers to [15,16] where the authors explain the relevant drawbacks found when perturbing the input flow in the classical deterministic chemostat (1.1)–(1.2) with a Wiener process (where μ is a Monod growth function).
Moreover, the approach proposed in this paper allows to prove the persistence of the bacterial species (under some conditions on the growth function), as it is observed by practitioners on very long time periods despite variations of the input flow rate. This is not the case when considering the Wiener process where persistence cannot be ensured (see [16] and [17] where the Wiener process is used to model disturbances on the input flow and environmental perturbations in the classical deterministic chemostat).
In industrial setup, large concentrations of the input substrate sin can be observed and it is also well-known that bacterial species may suffer from growth inhibition under very large concentrations s. The non-monotonic growth function (1.3) precisely models this fact (see [18]). Differently to the classical case, where the growth function μ is assumed to be increasing (as this is the case for the Monod function), the dynamics of the deterministic chemostat model with such non-monotonic growth function may exhibit a bi-stability for certain values of the dilution rate D (see for instance [11]). Depending on the initial condition, the state of the system converges asymptotically to the wash-out of the biomass (which is not a desirable state) or to a positive equilibrium. This kind of instability is observed in practice and present an issue in industrial applications because it requires a good monitoring of the system to detect if the state belongs the attraction basin of the wash-out equilibrium [19,20,21,22]. Most of the time, practitioners prefer to size the process to avoid such a behavior, i.e. such that the system admits an unique globally stable equilibrium (see also [23,24]). This mathematically amounts to have the following condition
D<min(μ(sin),μ(sm)) | (1.9) |
(this result is recalled later). The purpose of the present work is to study the behavior of the random dynamics when the constant removal rate is replaced by D+Φ(z∗(θtω)) where D satisfies condition (1.9).
The realizations of this variable may satisfy or not condition (1.9) at some times t. If not, one may wonder if this could lead the biomass to extinction. This question is of primer importance for the practitioners for the good health of the bio-process. In other words, for a nominal removal rate that satisfies condition (1.9), is the persistence of the biomass always guaranteed, even when the realizations of the noise provide effective values of the removal rate that do not satisfy this condition? Precise definitions of persistence in the framework of the chemostat will be given later.
Let us underline that stochastic modeling of the chemostat has received a great attention in the literature, considering different kinds of demographic noise [25,26,27,28,29,30], but few works have dealt with noise on the input, and much less in the case of considering Haldane consumption function, whereas this is quite natural for an open system as the chemostat, which is often the main source of fluctuations.
The paper is organized in the following way: in Section 2 we provide preliminaries and classical results about the deterministic chemostat model (1.1)–(1.2). In Section 3 we study the properties of the solutions of the random chemostat model with Haldane consumption kinetics (1.7)–(1.8). Then, in Section 4 we give conditions to ensure both uniform weak and strong persistence of the species. In Section 5, several simulations are presented to support the theoretical study. Finally, we present some conclusions in Section 6.
In this section we recall briefly some results that are useful in this paper. For the sake of clarity, we split this preliminary section in three different parts: the first one concerns classical results of the deterministic chemostat model (1.1)–(1.2), the second one is about the Ornstein-Uhlenbeck process and the third one recalls the definitions of persistence.
The next proposition recalls the classical results about the chemostat model (1.1)–(1.2) when μ is a non-monotonic function. We refer readers to [11,12] for proofs and more details.
Proposition 2.1. Assume that there exists ˆs∈(0,sin) such that the function μ is increasing on (0,ˆs) and decreasing on (ˆs,sin). Define the break-even concentrations λ−(D), λ+(D) as follows
λ−(D)=min{s∈[0,ˆs];μ(s)≥D},D∈[0,μ(ˆs)]λ+(D)=max{s∈[ˆs,sin];μ(s)≥D},D∈[μ(sin),μ(ˆs)] |
1. If D>μ(ˆs), the system (1.1)–(1.2) possesses a unique equilibrium E0:=(0,sin), which is globally asymptotically stable on R2+.
2. If D<μ(sin), the system (1.1)–(1.2) admits a unique positive equilibrium E−:=(sin−λ−(D),λ−(D)) which is globally asymptotically stable on R⋆+×R+.
3. If D∈[μ(sin),μ(ˆs)], the system (1.1)–(1.2) presents a bi-stability between E− and E0. From any initial condition in R⋆+×R+ excepted on a set of null measure, the solution converges asymptotically to E− or E0.
Remark 1. In practice, only the second case is desirable because it guarantees that in any situation the wash-out of the biomass is avoided.
Remark 2. For the Haldane expression (1.3), one has explicit expressions of the functions λ±
λ±(D)=(μ0−D)ki±√(μ0−D)2k2i−4D2kski2D. |
Let us recall that the concept of break-even concentrations has been revisited in the context of stochastic models of the chemostat [31,32] but we will not need it here. We keep the classical deterministic definition.
In the rest of the paper, we shall consider that we are in conditions of Proposition 2.1, that is μ non-monotonic on the interval [0,sin] (otherwise the analysis is similar to monotonic growth function and cannot present bi-stability). Throughout the paper we shall consider the following hypothesis
Assumption 1. There exists sm∈(0,sin) such that the function μ is increasing on (0,sm) and decreasing on (sm,sin).
We present here briefly the Ornstein-Uhlenbeck (OU) process. For more details we refer readers to [13,33,34].
The Ornstein-Uhlenbeck (OU) process is a stationary mean-reverting Gaussian stochastic process defined as
(t,ω)↦z(t,ω):=z∗(θtω)=−βν0∫−∞eβsθtω(s)ds,for allt∈R,ω∈Ω,β,ν>0, | (2.1) |
where ω denotes a standard Wiener process in a probability space (Ω,F,P), β is the mean reversion constant representing the strength with which the process is attracted by the mean, ν>0 is the volatility constant describing the variation or the size of the noise and θt denotes the usual Wiener shift flow given by
θtω(⋅)=ω(⋅+t)−ω(t),t∈R. |
We note that the OU process (2.1) can be obtained as the stationary solution of the Langevin equation
dz+βzdt=νdω. | (2.2) |
Typically, the OU process (2.1) can model the position of a particle by taking into account its friction in a fluid (which is the main difference with the typical standard Wiener process). Indeed, it can be considered as a generalization of the standard Wiener process and provides a link between the standard Wiener process (β=0, ν=1) and no noise at all (β=1, ν=0).
From now on we consider β and ν fixed and z∗(θtω) the OU process defined above.
We recall in the next proposition some of its properties.
Proposition 2.2 (See [35,36]). There exists a θt-invariant set ˜Ω∈F of Ω of full measure such that for ω∈˜Ω and β,ν>0, we have
(i) the random variable |z∗(ω)| is tempered with respect to {θt}t∈R, i.e., for a.e. ω∈˜Ω,
limt→∞e−ηtsupt∈R|z∗(θ−tω)|=0,for all η>0. |
(ii) this mapping is a stationary solution of (2.2) with continuous trajectories
(t,ω)→z∗(θtω)=−βν0∫−∞eβs(θtω)(s)ds; |
(iii) for any ω∈˜Ω one has:
limt→±∞|z∗(θtω)|t=0;limt→±∞1t∫t0z∗(θsω)ds=0; |
limt→±∞1t∫t0|z∗(θsω)|ds=E[|z∗|]<∞. |
We recall here the definitions of (uniform) persistence (see for instance [37]) that we consider in the present work.
Definition 2.1. The dynamics of (1.7)–(1.8) is strongly persistent if there exists ε>0 such that for any initial condition (s(0),x(0)) in R+×R⋆+ and any realization of the D+Φ(z∗(θ(⋅)ω)), the solution verifies
lim inft→+∞x(t)>ε. |
Definition 2.2. The dynamics of (1.7)–(1.8) is weakly persistent if there exists ε>0 such that for any initial condition (s(0),x(0)) in R+×R⋆+ and any realization of D+Φ(z∗(θ(⋅)ω)), the solution verifies
lim supt→+∞x(t)>ε. |
In this section we study the random chemostat model (1.7)–(1.8) presented in the introduction. We prove the existence and uniqueness of a global positive solution and provide results about the existence of absorbing and attracting sets which, in addition, are deterministic (i.e. that do not depend on the realization of the noise). In addition, we derive first conditions under which extinction of species cannot be avoided.
In the sequel, we denote X:={(s,x)∈R2:s,x≥0} for the positive cone.
Theorem 3.1. For any initial condition u(0)∈X, system (1.7)–(1.8) possesses a unique global solution u(t;0,ω,u(0))=(s(t;0,ω,u(0)),x(t;0,ω,u(0)))∈C1([0,+∞);X). In addition, it remains in X for every t>0.
Proof. Let us first write system (1.7)–(1.8) as
dudt=L(θtω)u+F(u,θtω), |
where
u=(sx),L(θtω)=(−(D+Φ(z∗(θtω)))00−(D+Φ(z∗(θtω)))) |
and
F(u,θtω)=((D+Φ(z∗(θtω)))sin−μ0sks+s+s2kixμ0sks+s+s2kix). |
F(⋅,θtω)∈C1(X×[0,+∞);X) whence F is locally Lipschitz respect to u∈X. Therefore, for each realization of the noise, du/dt=L(θtω)u+F(u,θtω) is an non-autonomous differential equation with a right member Lipschitz with respect to u and continuous with respect to t (recall that z∗(θtω) is continuous respect to t and Φ is also continuous). Therefore, the solution of the Cauchy problem admits an unique local solution of system (1.7)–(1.8) from the theory of ordinary differential equations.
Let us first show that both x and s remain in X for any u0∈X. To this end, we notice first that x(t)=0 solves (1.7). By uniqueness of the Cauchy problem, we deduce that any other solution is such that x(t)≠0 for any t. Notice also that one has, thanks to (1.6),
dsdt|s=0=(D+Φ(z∗(θtω)))sin>0, |
which proves that the axis s=0 is repulsive in X. This demonstrates the positiveness of the solution of system (1.7)–(1.8).
Now we prove that both the substrate and the microorganisms concentrations remain bounded for every time. To this end, define v(t):=s(t)+x(t). Then the variable v is solution of the following differential equation
dvdt=(D+Φ(z∗(θtω)))(sin−v) |
and, thanks to (1.6),
dvdt≤Drsin−Dlv, |
By comparison of solutions of scalar ODEs (see [38]) we obtain
v(t;0,ω,v(0))≤v(0)e−Dlt+DrDlsin(1−e−Dlt),t≥0. |
Then, v is forward bounded and since v=s+x≥0, we deduce that both the solution s and x are also bounded for positive time. The unique solution of system (1.7)–(1.8) is thus defined for any t≥0.
Now, we prove the existence of deterministic attracting and absorbing sets for the solutions of the random chemostat model (1.7)–(1.8).
Theorem 3.2. The system (1.7)–(1.8) possesses a deterministic (forward) attracting set
A:={(s,x)∈X:s+x=sin}. | (3.1) |
Proof. Define the variable q(t):=s(t)−sin+x(t). Then q satisfies the following differential equation
dqdt=−(D+Φ(z∗(θtω)))q |
whose solution is given by
q(t;0,ω,q(0))=q(0)e−∫t0(D+Φ(z∗(θsω)))ds. | (3.2) |
From (1.6), we have
q(0)e−Drt≤q(t;0,ω,q(0))≤q(0)e−Dlt |
and taking limit when t goes to infinity we have
limt→+∞q(t;0,ω,q(0))=0, |
which means that, for any ε>0, there exists T(ω,ε)>0 such that
q(t;0,ω,q(0))∈[−ε,ε],t≥T(ω,ε). |
This proves that
Bε:={(s,x)∈X:sin−ε≤s+x≤sin+ε}. |
is a (forward) deterministic absorbing set for system (1.7)–(1.8). As ε>0 is arbitrary, the proof is done.
Remark 3. Let us underline that the attracting set (3.1) obtained in Theorem 3.2 does not depend on the event ω∈Ω. This is another particularity that we obtain when considering this way of modeling random bounded realizations.
Our aim now is to provide first conditions on the parameters of system (1.7)–(1.8) under which extinction of the species occurs.
Theorem 3.3. Assume that the inequality
Dl>μ(sm) | (3.3) |
is fulfilled, where sm is defined in (1.4). Then, the singleton
Ae:={(sin,0)}⊂A |
is a (forward) attracting set.
Proof.
From the equation describing the dynamics of the species (1.8) and (1.6) one has
dxdt=−(D+Φ(z∗(θtω)))x+μ(s)x≤−(Dl−μ(sm))x. |
and by comparison of solutions of scalar ordinary differential equations [38], we have
x(t;0,ω,x(0))≤x(0)e−(Dl−μ(sm))t. | (3.4) |
Taking limit when t goes to infinity in (3.4), we obtain from (3.3) and Proposition 2.2 (ⅲ),
limt→+∞x(t;0,ω,x(0))=0. |
In this section our aim is to provide conditions to ensure persistence of the species in the sense defined in Section 2.3.
We shall distinguish in the following weak and strong persistence, as defined for instance in [37]. Before stating a result concerning the weak persistence of species, let us first provide a property linked to the choice of the function Φ.
Proposition 4.1. Let be Φ(z)=2dπarctan(z) defined as in (1.5). Then
limt→+∞1t∫t0Φ(z∗(θsω))ds=0,a.s. inΩ. | (4.1) |
Proof. Since Φ(z)=2dπarctan(z), then we have
∫Ω|2dπarctan(z∗(ω))|dP(ω)≤d|Ω|=d |
whence Φ(z∗(⋅))∈L1(Ω,F,P). Then, since P is invariant by θ (see [34,35]), from the Birkhoff ergodic theorem it yields that
limt→+∞1t∫t0Φ(z∗(θsω))ds=E[Φ(z∗(ω))],a.s. inΩ. |
Hence, it is enough to show that E[Φ(z∗(ω))]=0. In fact, we have
E[Φ(z∗(ω))]=∫RΦ(x)fOU(x)dx=0, |
where fOU denotes the density function of the random variable z∗(⋅), which is an even function since it is gaussian, and Φ is an odd function.
Remark 4. From Proposition 4.1 we notice that we could model random perturbations as in this paper, i.e., by means of Φ(z∗(θtω)), and the ergodic property (4.1) remains true as long as Φ is an odd measurable function.
Theorem 4.1. Assume that the inequality
μ(sin)>D. | (4.2) |
is fulfilled. Then, the random chemostat model (1.7)–(1.8) is weakly persistent, that is there exits ε>0 such that for any initial condition x(0)∈X, any realization satisfies
lim supt→+∞x(t;0,ω,x(0))≥ε. |
Proof. Let ε>0 be such that μ(s)>D for all s∈[sin−ε,sin+ε] and define
η:=min{μ(s)−D:s∈[sin−ε,sin+ε]}>0. |
From Theorem 3.2, we know that s(t;0,ω,s(0)) converges asymptotically to the set [0,sin]. Therefore, there exists T(ω,ε)>0 such that s(t;0,ω,s(0))<sin+ε for any t>T(ω,ε).
Consider now, for t>T(ω,ε), the sets
U(t):={τ∈[T(ω,ε),t]:s(τ;0,ω,s(0))<sin−ε},V(t):={τ∈[T(ω,ε),t]:s(τ;0,ω,s(0))∈[sin−ε,sin+ε]} |
and the functions in [0,1]
u(t):=measU(t)t−T(ω,ε),v(t):=measV(t)t−T(ω,ε)=1−u(t). |
Then, from (1.8), one can write the following inequality
˙x(t)≥{−(D+Φ(z∗(θtω)))x(t),t∈U(t)(η−Φ(z∗(θtω)))x(t),t∈V(t) |
and by integration between T(ω,ε) and t>T(ω,ε) one obtain
x(t)≥x(T(ω,ε))e(t−T(ω,ε))[u(t)(−D)+(1−u(t))η−1t−T(ω,ε)∫tT(ω,ε)Φ(z∗(θrω))dr],t>T(ω,ε). | (4.3) |
Assume that one has
limt→+∞u(t)=0. |
From Proposition 4.1, one has
limt→+∞e−1t−T(ω,ε)∫tT(ω,ε)Φ(z∗(θrω))dr=1 |
and then (4.3) gives
limt→+∞x(t)=+∞ |
which is a contradiction since x is bounded. We deduce thus that the sets U(⋅) are necessarily such that
measU(t)→+∞ when t→+∞ |
that is
meas({t:s(t)<sin−ε})=+∞ |
which shows that one has
lim inft→+∞s(t)≤sin−ε |
or equivalently
lim supt→+∞x(t)≥ε>0 |
since s(t)+x(t) converges to sin for any realization (see the proof of Theorem 3.2).
Remark 5. Theorem 4.1 proves the weak (uniform) persistence of species as long as D<μ(sin) is fulfilled, the same condition that guarantees the persistence in the deterministic case (see Proposition 2.1). However, let us underline that we do not impose the upper bound Dr of the variations of the removal rate to fulfill this inequality. This means that one could have realizations of the disturbances such that the effective value of the removal rate is above μ(sin) on large periods of time making the species to be arbitrary closed to the extinction but it will always persist.
In the next result we consider a stronger condition on the upper bound Dr of the removal rate which ensures the strong persistence of the species.
Theorem 4.2. Assume that the inequality
Dr<μ(sin) | (4.4) |
is fulfilled. Then, the random chemostat model (1.7)–(1.8) is strongly persistent and the set
Ap:={(x,y)∈X:s+x=sin,x≥sin−λ−(Dr)}⊂A |
is (forward) attracting, where the function λ− is defined in Proposition 2.1.
Proof. Take ε>0 such that μ(s)>Dr for any s∈[sin−ε,sin+ε] and posit
η:=min{μ(s)−Dr:s∈[sin−ε,sin+ε]}>0. |
From now on, we will s(t), x(t) and q(t) instead of s(t;0,ω,s(0)), x(t;0,ω,x(0)) and q(t;0,ω,q(0)) to make the readability easier even though we recall that every state variable depends on the noise.
From Theorem 3.2, we know that the variable q(t)=s(t)−sin+x(t) converges asymptotically to zero for any realization, and thus any solution s(t) of the random model (1.7)–(1.8) converges to the set [0,sin]. Consequently, there exists T(ω,ε)>0 such that
s(t)<sin+ε,q(t)>−κ,∀t>T(ω,ε), |
where
κ:=μ(sin−ε/2)−Drμ(sin−ε/2)(ε/2)>0. |
If s(t)∈[sin−ε,sin+ε] for any t>T(ω,ε), then one has, from (1.8), ˙x(t)>ηx(t) for any t>T(ω,ε), which implies that x is unbounded, thus a contradiction. We deduce that there exists a finite time ˉT(ω)≥T(ω,ε) such that s(ˉT(ω))≤sin−ε/2. On another hand, from (1.7), s(t) can be written as the solution of the non-autonomous dynamics
dσdt=F(t,σ):=(D+Φ(z∗(θtω))−μ(σ))(sin−σ)−μ(σ)q(t). | (4.5) |
Note that one has
F(t,sin−ε/2)≤(Dr−μ(sin−ε/2))(ε/2)−μ(sin−ε/2)q(t) |
and then
F(t,sin−ε/2)≤(Dr−μ(sin−ε/2))(ε/2)+μ(sin−ε/2)κ=0,t>T(ω,ε) |
(from the definition of κ). We deduce that the set [0,sin−ε/2] is forward invariant for the semi-flow {˙σ=F(t,σ),t>T(ω,ε)}. Therefore, one has
s(t)≤sin−ε/2,t>ˉT(ω). |
Then, once identified σ with s, from (4.5)
dsdt≤(Dr−μ(s))(sin−s)−μ(s)q(t),t>ˉT(ω) |
and from the comparison of solutions of scalar ordinary differential equations [38], one has the inequality s(t)≤s+(t) for any t>ˉT(ω), where s+(t) is solution of the Cauchy problem
ds+dt=(Dr−μ(s+))(sin−s+)−μ(s+)q(t),s+(ˉT(ω))=s(ˉT(ω)). |
Note that the solution s+(t) belongs to the interval [0,sin] for any t≥ˉT(ω) (and is thus bounded) and that its dynamics is asymptotic autonomous with limiting dynamics
ds†dt=(Dr−μ(s†))(sin−s†) | (4.6) |
Under assumption Dr<μ(sin), one has necessarily λ−(Dr)<sin and the property
(Dr−μ(σ))(λ−(Dr)−σ)>0,∀σ∈[0,sin]∖{λ−(Dr)} |
is fulfilled. One finally obtains that any solution of (4.6) in [0,sin] is such that s†(t)→λ−(Dr) when t→+∞. From the theory of asymptotically autonomous dynamical systems [39], one concludes that s+(t) converges also to λ−(Dr) when t→+∞, which proves that one has
lim supt→+∞s(t)≤λ−(Dr) |
or equivalently
lim inft→+∞x(t)≥sin−λ−(Dr)>0 |
since s(t)+x(t) converges to sin for any realization. This demonstrates the strong persistence of the random dynamics (1.7)–(1.8) with the explicit lower bound sin−λ−(Dr).
Now we consider the last situation when μ(sin)<D<μ(sm) which corresponds to the bi-stability in the deterministic case. In the random framework, one cannot guarantee persistence nor wash-out of the biomass. As we shall see later on simulations, the asymptotic behavior of the solutions depends on the initial condition and the realization of the noise. However we show that an upper bound on the biomass can be provided.
Theorem 4.3. Assume x(0)>0. Then, the upper bound
lim supt→+∞x(t;0,ω,x(0))≤sin−λ−(Dl) |
for any solution of the random dynamics (1.7)–(1.8) as long as Dl<μ(sm).
Proof. Let us remark that in this proof we will write x(t), s(t) and q(t), or simply x, s and q, instead of x(t;0,ω,x(0)), s(t;0,ω,s(0)) and q(t;0,ω,q(0)) for the sake of simplicity.
Consider x(0)>0. Then, from (1.8) one has
dxdt=(μ(s)−(D+Φ(z∗(θtω))))x=(μ(sin−x−q)−(D+Φ(z∗(θtω))))x:=F(t,x), | (4.7) |
where q=sin−x−s.
We know that limt→+∞q(t)=0 then, for every ε>0 there exists T(ω,ε)>0 such that |q(t)|<ε for every t>T(ω,ε). Thus x(t)<sin+ε for every t>T(ω,ε).
Moreover, for all x∈(sin−λ−(Dl)+ε,sin+ε), one has sin−x−q∈(−q−ε,λ−(Dl)−q−ε). In fact, since sin−x−q=s≥0 (see Theorem 3.1), then we have sin−x−q∈(0,λ−(Dl)−q−ε)⊂(0,λ−(Dl)) whence F(t,x)<0 since μ(0,λ−(Dl))<Dl and D+Φ(z∗(θtω))≥Dl.
Thus,
lim supt→+∞x(t)≤sin−λ−(Dl)+ε |
for every ε>0.
In this section we present several numerical simulations to support our theoretical results in three different cases: extinction, strong persistence and weak persistence of species.
Two different figures are presented in each case. The first figure (with two panels) concerns the dynamics of the substrate (top) and the species (bottom) in the random chemostat model (1.7)–(1.8) for different realizations of the noise (in colored continuous lines) along with the case without noise i.e. the deterministic chemostat model (in blue dashed lines). In addition, in each panel we display a little box with a zoom to illustrate better the dynamics around the attracting set. The second figure offers information about the consumption function of the consumer species and those parameters involved in conditions to have extinction, weak persistence and strong persistence of species.
We first present the case of extinction in Figure 2. We consider s(0)=14, x(0)=5 as initial conditions and set sin=14, μ0=4, ks=7.5, ki=4, D=1.4 and the interval given by practitioners from observations is [Dl,Dr]=[1.15,1.65], then d=0.25. In addition, ν=1 and β=1. We can observe in Figure 3 that Dl>μ(sm) holds true and then, as proved in Theorem 3.3, the species extinguishes (the attracting set in this case is {(0,sin)}).
This is indeed quite intuitive as condition (3.3) basically means that the nominal value D of the removal rate is too large large compared to the growth kinetics (recall that Dl>μ(sm) implies D>μ(sm)), even though we start far away from the washout (sin,0).
Now we present situations where weak and strong persistence is ensured. To compare, we start with the case of weak persistence by considering s(0)=16.7, x(0)=0.01 as initial conditions and set sin=16.7, μ0=7, ks=7, ki=7. In addition, β=1 and ν=1. Once fixed these parameters, the removal rate D will be the set to allow conditions (4.4) and (4.2) be true or false.
In Figure 4 we set D=1.7 and the interval [Dl,Dr]=[1.45,1.95] then d=0.25. In this case D<μ(sin) is fulfilled (observe Figure 5) then we have weak persistence of species, see Theorem 4.1. Although we start very closed to the washout (sin,0), we observe persistence of species.
In Figure 6 we set D=1.4 and the interval [Dl,Dr]=[1.15,1.65] then d=0.25. In this case Dr<μ(sin) holds true (see Figure 10) and then strong persistence of species is obtained from Theorem 4.1.
Moreover, we would like to notice that OU process proves again to be a powerful tool when modeling real (bounded) noises. In addition, it allows us to observe clearly the difference between weak and strong persistence.
Finally we present some numerical simulations to observe that both extinction and persistence can be obtained when Dl<μ(sin)<D is fulfilled, as in the deterministic case. However, the issue here depends also on the realization of the noise.
We have considered the chemostat model (1.1)–(1.2) with Haldane consumption kinetics, under bounded perturbations on the input flow rate, motivated by real cases in industrial setup and biotechnology. To this end, we use a bounded saturated function of the Ornstein-Uhlenbeck process that allows us to ensure the noise to be bounded in realistic intervals.
We prove in Theorem 3.1 existence and uniqueness of global positive solution of the corresponding random chemostat (1.7)–(1.8) by means of standard results from the theory of ODEs and thanks to properties of the Ornstein-Uhlenbeck process. Then, in Theorem 3.2 we establish the existence of an absorbing and attracting set which has the nice property to be deterministic, i.e. that does not depend on the realization of the noise.
We then focused on the long-time behavior of the random dynamics inside this attracting set. To this end, we first proved in Theorem 3.3 that extinction of species cannot be avoided as long as Dl>μ(sm) whatever is the input concentration sin. On the opposite, we proved in Theorem 4.1 the weakly uniform persistence of species when D<μ(sin), which means that species can be temporarily arbitrary closed to the extinction but still persist even when having random disturbances in the input flow. The condition D<μ(sin) ensures persistence in the deterministic case but the effective removal rate does not necessarily fulfill this condition depending on the realizations. Finally, we prove the strong persistence of the species in Theorem 4.2 under the stronger condition Dr<μ(sin). In this case, we provide an explicit lower bound for the asymptotic concentration of the species, a useful information for practitioners.
In addition, we support the theoretical results with several numerical simulations which depict the possible behaviors of the random dynamics. Moreover, this allows us to illustrate the difference between weak and strong persistence: once fixed every parameter we change the removal rate D to let conditions μ(sin)>D and μ(sin)>Dr be true or not. Condition to have weak persistence (μ(sin)>D) does not involve the noise and then the random realizations of the species can fluctuate closed to zero and we still ensure persistence whereas condition to have strong persistence (μ(sin)>Dr) only fulfills for realizations of the perturbed removal rate inside a small enough interval.
To conclude, we would like to stress that the Ornstein-Uhlenbeck process has proved once again to be relevant tool to model random disturbances in a biological framework. This stochastic process fits in a quite loyal way the bounded fluctuations that are observed in practice, and justify the (weak or strong) persistence of the biomass despite the possible realizations of noise, as also observed in practice.
This work has been partially supported by Proyecto de Excelencia P12-FQM-1492 from Junta de Andalucía, Project PGC2018-096540-B-I00 from Ministerio de Ciencia, Innovación y Universidades (Spanish government) and Project US-1254251 from Consejería de Economía y Conocimiento (Junta de Andalucía).
The authors declare there is no conflict of interests.
[1] | J. Monod, La technique de culture continue: Théorie et applications, Annales de l'Institute Pasteur, 79 (1950), 390-410. |
[2] |
A. Novick, L. Szilard, Experiments with the chemostat on spontaneous mutations of bacteria, Proc. Natl. Acad. Sci., 36 (1950), 708-719. doi: 10.1073/pnas.36.12.708
![]() |
[3] |
H. W. Jannasch, Steady state and the chemostat in ecology, Limnol. Oceanogr., 19 (1974), 716- 720. doi: 10.4319/lo.1974.19.4.0716
![]() |
[4] |
G. D'Ans, P. Kokotovic, D. Gottlieb, A nonlinear regulator problem for a model of biological waste treatment, IEEE Trans. Autom. Control., 16 (1971), 341-347. doi: 10.1109/TAC.1971.1099745
![]() |
[5] | J. W. M. La Rivière, Microbial ecology of liquid waste treatment, in Advances in Microbial Ecology, 1 (1977), 215-259, Springer US. |
[6] | J. Barlow, F. de Noyelles, B. Peterson, J. Peterson, W. Schaffner, "Continuous flow nutrient bioassays with natural phytoplankton populations". G. Glass (Editor): Bioassay Techniques and Environmental Chemistry, John Wiley & Sons Ltd., 1973. |
[7] |
I. F. Creed, D. M. McKnight, B. A. Pellerin, M. B. Green, B. A. Bergamaschi, G. R. Aiken, et al., The river as a chemostat: fresh perspectives on dissolved organic matter flowing down the river continuum, Can. J. Fish. Aquatic Sci., 72 (2015), 1272-1285. doi: 10.1139/cjfas-2014-0400
![]() |
[8] | S. Jorgensen, B. Fath, Fundamentals of Ecological Modelling Applications in Environmental Management and Research. Elsevier, 2011. |
[9] |
J. Kalff, R. Knoechel, Phytoplankton and their dynamics in oligotrophic and eutrophic lakes, Annu. Rev. Ecol. Syst., 9 (1978), 475-495. doi: 10.1146/annurev.es.09.110178.002355
![]() |
[10] |
E. Rurangwa, M. C. J. Verdegem, Microorganisms in recirculating aquaculture systems and their management, Rev. Aquac., 7 (2015), 117-130. doi: 10.1111/raq.12057
![]() |
[11] | J. Harmand, C. Lobry, A. Rapaport, T. Sari, The Chemostat: Mathematical Theory of Microorganisms Cultures. Wiley, Chemical Engineering Series, John Wiley & Sons, Inc., 2017. |
[12] | H. L. Smith, P. Waltman, The theory of the chemostat: dynamics of microbial competition. Cambridge University Press, 1995. |
[13] | T. Caraballo, M. J. Garrido-Atienza, J. López-de-la-Cruz, A. Rapaport, Modeling and analysis of random and stochastic input flows in the chemostat model, Discret. Continuous Dyn. Syst. Ser. B, 24 (2018), 3591-3614. |
[14] | Y. Asai, P. Kloeden, Numerical schemes for random odes via stochastic differential equations, Commun. Appl. Analysis, 17 (2013), 521-528. |
[15] |
T. Caraballo, X. Han, P. E. Kloeden, Chemostats with random inputs and wall growth, Math. Methods Appl. Sci., 38 (2015), 3538-3550. doi: 10.1002/mma.3437
![]() |
[16] | T. Caraballo, M. J. Garrido-Atienza, J. López-de-la-Cruz, Some Aspects Concerning the Dynamics of Stochastic Chemostats, vol. 69, ch. 11, pp. 227-246. Springer International Publishing, 2016. |
[17] |
T. Caraballo, M. J. Garrido-Atienza, J. López-de-la-Cruz, Dynamics of some stochastic chemostat models with multiplicative noise, Commun. Pure Appl. Analysis, 16 (2017), 1893-1914. doi: 10.3934/cpaa.2017092
![]() |
[18] |
J. Andrews, A mathematical model for the continuous culture of microorganisms utilizing inhibitory substrates, Biotechnol. Bioeng., 10 (1968), 707-723. doi: 10.1002/bit.260100602
![]() |
[19] | G. Bastin, D. Dochain, On-line estimation and adaptive control of bioreactors. Elsevier, 1990. |
[20] |
A. Rapaport, J. Harmand, Robust regulation of a class of partially observed nonlinear continuous bioreactors, J. Process. Control., 12 (2002), 291-302. doi: 10.1016/S0959-1524(01)00029-4
![]() |
[21] |
B. Satishkumar, M. Chidambaram, Control of unstable bioreactor using fuzzy tuned PI controller, Bioprocess Eng., 20, (1999), 127-132. doi: 10.1007/s004490050570
![]() |
[22] |
A. Schaum, J. Alvarez, T. Lopez-Arenas, Saturated PI control of continuous bioreactors with haldane kinetics, Chem. Eng. Sci., 68 (2012), 520-529. doi: 10.1016/j.ces.2011.10.006
![]() |
[23] | A. Rapaport, I. Haidar, J. Harmand, Global dynamics of the buffered chemostat for a general class of response functions, J. Math. Biol., 71 (2014), 69-98. |
[24] |
A. Rapaport, J. Harmand, Biological control of the chemostat with nonmonotonic response and different removal rates, Math. Biosci. Eng., 5 (2008), 539-547. doi: 10.3934/mbe.2008.5.539
![]() |
[25] |
B. Cloez, C. Fritsch, Gaussian approximations for chemostat models in finite and infinite dimensions, J. Math. Biol., 75 (2017), 805-843. doi: 10.1007/s00285-017-1097-6
![]() |
[26] |
P. Collet, S. Martínez, S. Méléard, J. S. Martín, Stochastic models for a chemostat and long-time behavior, Adv. Appl. Probab., 45 (2013), 822-836. doi: 10.1017/S0001867800006595
![]() |
[27] |
C. Fritsch, J. Harmand, F. Campillo, A modeling approach of the chemostat, Ecol. Model., 299 (2015), 1-13. doi: 10.1016/j.ecolmodel.2014.11.021
![]() |
[28] |
J. Grasman, M. D. Gee, O. A. V. Herwaarden, Breakdown of a chemostat exposed to stochastic noise, J. Eng. Math., 53 (2005), 291-300. doi: 10.1007/s10665-005-9004-3
![]() |
[29] |
L. Imhof, S. Walcher, Exclusion and persistence in deterministic and stochastic chemostat models, J. Differ. Equations, 217 (2005), 26-53. doi: 10.1016/j.jde.2005.06.017
![]() |
[30] |
G. Stephanopoulos, R. Aris, A. Fredrickson, A stochastic analysis of the growth of competing microbial populations in a continuous biochemical reactor, Math. Biosci., 45 (1979), 99-135. doi: 10.1016/0025-5564(79)90098-1
![]() |
[31] |
C. Xu, S. Yuan, An analogue of break-even concentration in a simple stochastic chemostat model, Appl. Math. Lett., 48 (2015), 62-68. doi: 10.1016/j.aml.2015.03.012
![]() |
[32] |
D. Zhao, S. Yuan, Critical result on the break-even concentration in a single-species stochastic chemostat model, J. Math. Anal. Appl., 434 (2016), 1336-1345. doi: 10.1016/j.jmaa.2015.09.070
![]() |
[33] |
T. Caraballo, R. Colucci, J. López-de-la-Cruz, A. Rapaport, A way to model stochastic perturba-tions in population dynamics models with bounded realizations, Commun. Nonlinear Sci. Numer. Simul., 77 (2019), 239-257. doi: 10.1016/j.cnsns.2019.04.019
![]() |
[34] | T. Caraballo, X. Han, Applied Nonautonomous and Random Dynamical Systems, Applied Dynamical Systems. Springer International Publishing, 2016. |
[35] | L. Arnold, Random Dynamical Systems. Springer Berlin Heidelberg, 1998. |
[36] |
T. Caraballo, P. E. Kloeden, B. Schmalfuss, Exponentially stable stationary solutions for stochastic evolution equations and their perturbation, Appl. Math. Optim., 50 (2004), 183-207. doi: 10.1007/s00245-004-0802-1
![]() |
[37] |
H. I. Freedman, P. Moson, Persistence definitions and their connections, Proc. Am. Math. Soc., 109 (1990), 1025-1033. doi: 10.1090/S0002-9939-1990-1012928-6
![]() |
[38] | W. Walter, Ordinary Differential Equations. Springer New York, 1998. |
[39] | H. R. Thieme, Convergence results and a poincare-bendixson trichotomy for asymptotically autonomous differential equations, J. Math. Biol., 30 (1992), 755-763. |
1. | Tomás Caraballo, Javier López-de-la-Cruz, Bounded random fluctuations on the input flow in chemostat models with wall growth and non-monotonic kinetics, 2021, 6, 2473-6988, 4025, 10.3934/math.2021239 | |
2. | Tomás Caraballo, Javier López-de-la-Cruz, Survey on chemostat models with bounded random input flow, 2021, 1, 2767-8946, 52, 10.3934/mmc.2021005 | |
3. | Mateusz Czyżniewski, Rafał Łangowski, A robust sliding mode observer for non-linear uncertain biochemical systems, 2022, 123, 00190578, 25, 10.1016/j.isatra.2021.05.040 | |
4. | Tomás Caraballo, Javier López-de-la-Cruz, Verónica Caraballo-Romero, Non-autonomous chemostat models with non-monotonic growth, 2022, 0020-7160, 1, 10.1080/00207160.2022.2148467 | |
5. | Tomás Caraballo, José A. Langa, Alexandre N. Carvalho, Alexandre N. Oliveira-Sousa, Continuity and topological structural stability for nonautonomous random attractors, 2022, 22, 0219-4937, 10.1142/S021949372240024X | |
6. | Tomás Caraballo, Javier López-de-la-Cruz, Alain Rapaport, Study of the dynamics of two chemostats connected by Fickian diffusion with bounded random fluctuations, 2022, 22, 0219-4937, 10.1142/S0219493722400020 | |
7. | Yu Mu, Zuxiong Li, Bifurcation dynamics of a delayed chemostat system with spatial diffusion, 2023, 205, 03784754, 186, 10.1016/j.matcom.2022.09.022 | |
8. | Alejandro Rincón, Fredy E. Hoyos, Gloria Restrepo, Global stability of a continuous bioreactor model under persistent variation of the dilution rate, 2022, 20, 1551-0018, 3396, 10.3934/mbe.2023160 | |
9. | Tomás Caraballo, Javier López‐de‐la‐Cruz, Verónica Caraballo‐Romero, Chemostat models with Monod and Haldane consumption functions and random environmental fluctuations, 2023, 0170-4214, 10.1002/mma.9061 | |
10. | Huijian Zhu, Yuming Peng, Yiyang Li, Caibin Zeng, Forward dynamics and memory effect in a fractional order chemostat minimal model with non-monotonic growth, 2023, 0, 1937-1632, 0, 10.3934/dcdss.2023019 | |
11. | Hongcui Wang, Chaoqun Xu, Existence and Internal Structure of the Deterministic Attracting Set for a Random Ant Colonies Model, 2024, 23, 0219-4775, 10.1142/S0219477524500068 | |
12. | Javier López-de-la-Cruz, Alexandre N. Oliveira-Sousa, SIR models with vital dynamics, reinfection, and randomness to investigate the spread of infectious diseases, 2025, 140, 10075704, 108359, 10.1016/j.cnsns.2024.108359 | |
13. | Huijian Zhu, Lijie Li, Weiquan Pan, Extinction and strong persistence in the Beddington–DeAngelis predator–prey random model, 2023, 46, 0170-4214, 19351, 10.1002/mma.9629 | |
14. | Tomás Caraballo, Javier López-de-la-Cruz, Verónica Caraballo-Romero, Effects of real random perturbations on Monod and Haldane consumption functions in the chemostat model, 2024, 218, 03784754, 482, 10.1016/j.matcom.2023.11.035 |