Tomas Caraballo

Renato Colucci

Javier Lopez

La Cruz

Tomás Caraballo

Javier López-De-La-Cruz

Alain Rapaport

Study of the chemostat model with non-monotonic growth under random disturbances on the removal rate

Keywords: chemostat model, non-monotonic growth, bounded noise, Ornstein-Uhlenbeck, absorbing set

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

The chemostat refers to a laboratory device used for the growth of micro-organisms in a culture environment [START_REF] Monod | La technique de culture continue: Théorie et applications[END_REF][START_REF] Novick | Experiments with the chemostat on spontaneous mutations of bacteria[END_REF], that has been regarded as an idealization of the nature to study microbial ecosystems in stationary stage [START_REF] Jannasch | Steady state and the chemostat in ecology[END_REF]. It turned out to be an important investigation field due to a large number of applications, especially in waste water treatment [START_REF] D'ans | A nonlinear regulator problem for a model of biological waste treatment[END_REF][START_REF] La Rivière | Microbial ecology of liquid waste treatment[END_REF] but also in ecological and environmental sciences (see [START_REF] Barlow | Continuous flow nutrient bioassays with natural phytoplankton populations[END_REF][START_REF] Creed | The river as a chemostat: fresh perspectives on dissolved organic matter flowing down the river continuum[END_REF][START_REF] Jorgensen | Fundamentals of Ecological Modelling Applications in Environmental Management and Research[END_REF][START_REF] Kalff | Phytoplankton and their dynamics in oligotrophic and eutrophic lakes[END_REF][START_REF] Rurangwa | Microorganisms in recirculating aquaculture systems and their management[END_REF]).

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" [START_REF] Harmand | The Chemostat: Mathematical Theory of Micro-organisms Cultures[END_REF][START_REF] Smith | The theory of the chemostat: dynamics of microbial competition[END_REF] 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.

dx dt = -Dx + µ(s)x, (1.2)
where s = s(t) and x = x(t) denote the concentration of substrate and species, respectively, s in 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) = µ 0 s

s + k s + s 2 k i , s ≥ 0, (1.3)
where k s is the affinity constant and k i a parameter modeling the growth inhibition under large submitted to Mathematical Biosciences and Engineering concentrations of substrate. Since it will be useful later, let us define the number s m := arg max µ(s) = k i k s .

(1.4)

Many works have been dedicated to this classical deterministic model of the chemostat (see for instance [START_REF] Harmand | The Chemostat: Mathematical Theory of Micro-organisms Cultures[END_REF][START_REF] Smith | The theory of the chemostat: dynamics of microbial competition[END_REF]) 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 [START_REF] Caraballo | Modeling and analysis of random and stochastic input flows in the chemostat model[END_REF] where some chronicles of time varying removal rates are depicted).

Motivated by this fact, it has been proposed a modeling of 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, [START_REF] Asai | Numerical schemes for random odes via stochastic differential equations[END_REF][START_REF] Caraballo | Chemostats with random inputs and wall growth[END_REF] for other possible functions).

In this way, once practitioners provide us with an interval [D l , D r ] ⊂ R, 0 < D l < D < D r < ∞, (typically obtained from observations) we can define d = D r -D = D -D l and then the perturbed input flow is bounded for every time and any realization of the noise, i.e.,

D l ≤ D + Φ(z * (θ t ω)) ≤ D r ∀t ∈ R. (1.6)
Then, the random chemostat model writes as follows

ds dt = (D + Φ(z * (θ t ω)))(s in -s) -µ(s)x, (1.7
)

dx dt = (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 the mathematical and biological points of view, see for instance [START_REF] Caraballo | Chemostats with random inputs and wall growth[END_REF].

On the one hand, this way of modeling noise fits in a loyal way the bounded variations of the input flow r ate o bserved i n r eal l ife. O n t he c ontrary, t he W iener p rocess i s u nbounded with probability one which leads to arbitrary large (possibly negative) values of the corresponding perturbed input flow r ate, w hich i s n ot r ealistic a t a ll f rom t he b iological p oint o f v iew. We refer readers to [START_REF] Caraballo | López-de-la-Cruz[END_REF][START_REF] Caraballo | Chemostats with random inputs and wall growth[END_REF] 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 r ate. T his i s n ot t he c ase when considering the Wiener process where persistence cannot be ensured (see [START_REF] Caraballo | López-de-la-Cruz[END_REF] and [START_REF] Caraballo | Dynamics of some stochastic chemostat models with multiplicative noise[END_REF] 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 s in 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 [START_REF] Andrews | A mathematical model for the continuous culture of microorganisms utilizing inhibitory substrates[END_REF]). 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 [START_REF] Harmand | The Chemostat: Mathematical Theory of Micro-organisms Cultures[END_REF]). 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 a to positive equilibrium. This kind of instability is observed in practice and present an issue in industrial applications because its requires a good monitoring of the system to detect if the state belongs the attraction basin of the wash-out equilibrium [START_REF] Bastin | On-line estimation and adaptive control of bioreactors[END_REF][START_REF] Rapaport | Robust regulation of a class of partially observed nonlinear continuous bioreactors[END_REF][START_REF] Satishkumar | Control of unstable bioreactor using fuzzy tuned PI controller[END_REF][START_REF] Schaum | Saturated PI control of continuous bioreactors with haldane kinetics[END_REF]. 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 [START_REF] Rapaport | Global dynamics of the buffered chemostat for a general class of response functions[END_REF][START_REF] Rapaport | Biological control of the chemostat with nonmonotonic response and different removal rates[END_REF]). This mathematically amounts to have the following condition

D < min(µ(s in ), µ(s m )) (1.9)
(this result is recalled later one). The purpose of the present work is to study the behavior of the random dynamics when the constant removal rate is replaced by a bounded random process D+Φ(z * (θ t ω)), which is supposed to be continuous, where D satisfies the condition (1.9).

The realizations of this variable may satisfy or not the 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 c ondition ( 1.9), i s t he p ersistence o f t he b iomass always guaranteed, even when the realizations of the noise provide effective values of the removal rate that do not satisfy this condition? Precise definitions o f p ersistence i n t he f ramework o f the chemostat will be given later one.

Let us underline that stochastic modeling of the chemostat has received a great attention in the literature, considering different kinds of demographic noise [START_REF] Cloez | Gaussian approximations for chemostat models in finite and infinite dimensions[END_REF][START_REF] Collet | Stochastic models for a chemostat and long-time behavior[END_REF][START_REF] Fritsch | A modeling approach of the chemostat[END_REF][START_REF] Grasman | Breakdown of a chemostat exposed to stochastic noise[END_REF][START_REF] Imhof | Exclusion and persistence in deterministic and stochastic chemostat models[END_REF][START_REF] Stephanopoulos | A stochastic analysis of the growth of competing microbial populations in a continuous biochemical reactor[END_REF], 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.

Preliminaries

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 deterministic chemostat model

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 [START_REF] Harmand | The Chemostat: Mathematical Theory of Micro-organisms Cultures[END_REF][START_REF] Smith | The theory of the chemostat: dynamics of microbial competition[END_REF] for proofs and more details. Proposition 2.1. Assume that there exists ŝ ∈ (0, s in ) such that the function µ is increasing on (0, ŝ) and decreasing on (ŝ, s in ). Define the break-even concentrations λ -(D), λ + (D) as follows

λ -(D) = min{s ∈ [0, ŝ] ; µ(s) ≥ D}, D ∈ [0, µ(ŝ)] λ + (D) = max{s ∈ [ŝ, s in ] ; µ(s) ≥ D}, D ∈ [µ(s in ), µ(ŝ)] 1. If D > µ(ŝ), the system (1.1)-(1.2) possesses an unique equilibrium E 0 := (0, s in ), which is globally asymptotically stable on R 2 + . 2. If D < µ(s in ), the system (1.1)-(1.2) admits an unique positive equilibrium E -:= (s in - λ -(D), λ -(D)) which is globally asymptotically stable on R + × R + . 3. If D ∈ [µ(s in ), µ(ŝ)], the system (1.1)-(1.
2) presents a bi-stability between E -and E 0 . From any initial condition in R + × R + excepted on a set of null measure, the solution converges asymptotically to E -or E 0 .

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)k i ± (µ 0 -D) 2 k 2 i -4D 2 k s k i 2D .
Let us recall that the concept of break-even concentrations has been revisited in the context of stochastic models of the chemostat [START_REF] Xu | An analogue of break-even concentration in a simple stochastic chemostat model[END_REF][START_REF] Zhao | Critical result on the break-even concentration in a single-species stochastic chemostat model[END_REF] 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, s in ] (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 s m ∈ (0, s in ) such that the function µ is increasing on (0, s m ) and decreasing on (s m , s in ).

The Ornstein-Uhlenbeck process

We present here briefly the Ornstein-Uhlenbeck (OU) process. For more details we refer readers to [START_REF] Caraballo | A way to model stochastic perturbations in population dynamics models with bounded realizations[END_REF][START_REF] Caraballo | Modeling and analysis of random and stochastic input flows in the chemostat model[END_REF][START_REF] Caraballo | Applied Nonautonomous and Random Dynamical Systems[END_REF].

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 all t ∈ 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 [2, 12]

). 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. ω ∈ Ω,

lim t→∞ e -ηt sup t∈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:

lim t→±∞ |z * (θ t ω)| t = 0; lim t→±∞ 1 t t 0 z * (θ s ω)ds = 0; lim t→±∞ 1 t t 0 |z * (θ s ω)| ds = E [|z * |] < ∞.

Persistence in the chemostat

We recall here the definitions of (uniform) persistence (see for instance [START_REF] Freedman | Persistence definitions and their connections[END_REF]) that we consider in the present work .

Properties of the solutions of the random dynamics

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) ∈ R 2 : 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 )) ∈ C 1 ([0, +∞); X ). In addition, it remains in X for every t > 0.

Proof. Let us first write system (1.7)-(1.8) as

du dt = L(θ t ω)u + F (u, θ t ω), where u = s x , L(θ t ω) = -(D + Φ(z * (θ t ω))) 0 0 -(D + Φ(z * (θ t ω)))
and

F (u, θ t ω) =          (D + Φ(z * (θ t ω)))s in -µ 0 s s + k s + s 2 k i x µ 0 s s + k s + s 2 k i x          . F (•, θ t ω) ∈ C 1 (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 u 0 ∈ 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),

ds dt s=0 = (D + Φ(z * (θ t ω)))s in > 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

dv dt = (D + Φ(z * (θ t ω)))(s in -v) ≤ D r s in -D l v
and one, thanks to (1.6),

dv dt ≤ D r s in -D l v,
By comparison of solutions of scalar ODEs (see [START_REF] Walter | Ordinary Differential Equations[END_REF]) we obtain

v(t; 0, ω, v 0 ) ≤ v 0 e -D l t + D r D l s in 1 -e -D l t , 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 = s in } . (3.1)
Proof. Define the variable q(t) := s(t) -s in + x(t). Then q satisfies the following differential equation

dq dt = -(D +Φ(z * (θ t ω)))q
whose solution is given by q(t; 0, ω, q(0)) = q(0)e - From (1.6), we have q(0)e -Dr t ≤ q(t; 0, ω, q(0)) ≤ q(0)e -D l t and taking limit when t goes to infinity we have lim t→+∞ 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 : s in -ε ≤ s + x ≤ s in + ε} .
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 get 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.

dx dt = -(D +Φ(z * (θ t ω)))x + µ(s)x (3.4) ≤ -(D l -µ(s m ))x.
and by comparison of solutions of scalar ordinary di erential equations [START_REF] Walter | Ordinary Differential Equations[END_REF], we have x(t; 0, ω, x(0)) ≤ x(0)e -(D l -µ(sm))t .

Taking limit when t goes to infinity in (3.4), we get from (3.3) and Proposition 2.2 (iii), lim t→+∞ x(t; 0, ω, x(0)) = 0.

submitted to Mathematical Biosciences and Engineering

Conditions for persistence

In this section our aim is to provide conditions to ensure persistence of the species in the sense defined in Section 2.3.

Before stating a result concerning the weak persistence of species, let us provide now the following property. Proof. SinceΦ(z) = -2d π arctan(z), then we have

Φ - 2d π arctan(z * (ω)) dP(ω) ≤ d|Φ| = d whenceΦ (z * (ω)) ∈ L 1 (Φ, F, P).
Then, since P is invariant by θ t ω (see ??), from the Bircho ergodic theorem it yields that lim t→+∞

1 t t 0 Φ(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)f OU (x)dx = 0,
where f OU denotes the density function of the Ornstein-Uhlenbeck process, which is an odd function since it is gaussian, and Φ is an even 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 even function. 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 sup t→+∞ x(t; 0, ω, x 0 ) ≥ ε.

Proof. Let ε > 0 be such that µ(s) > D for all s ∈ [s in -ε, s in + ε] and define

η := min{µ(s) -D : s ∈ [s in -ε, s in + ε]} > 0.

submitted to Mathematical Biosciences and Engineering

From Theorem 3.2, we know that s(t; 0, ω, s 0 ) onverges asymptotically to the set [0, s in ]. Therefore, there exists T (ω, ε) > 0 such that s(t; 0, ω, s 0 ) < s in + ε for any t > T (ω, ε).

Consider now, for t > T (ω, ε), the sets

U (t) := {τ ∈ [T (ω, ε), t] : s(τ ; 0, ω, s 0 ) < s in -ε} , V (t) := {τ ∈ [T (ω, ε), t] : s(τ ; 0, ω, s 0 ) ∈ [s in -ε, s in + ε]} and the functions in [0, 1] u(t) := meas U (t) t -T (ω, ε) , v(t) := meas V (t) t -T (ω, ε) = 1 -u(t).
Then, from equation (1.8), one can write the following inequality

ẋ(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 get Remark 5. Theorem 4.1 proves the weak (uniformly) persistence of species as long as D < µ(s in ) 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 D r 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 µ(s in ) on large periods of time making the species to be arbitrary closed to the extinction but it will always persist.

x(t) ≥ x(T (ω, ε))e (t-T (ω,ε)) u(t)(-D)+(1-u(t))η- 1 t-T (ω,ε) t T (ω,ε) Φ(z * (θrω))dr , t > T (ω, ε) (4.
In the next result we consider a stronger condition on the upper bound D r of the removal rate which ensures the strong persistence of the species.

Theorem 4.2. Assume that the inequality

D r < µ(s in ) (4.4)
is fulfilled. Then, the random chemostat model (1.7)-(1.8) is strongly persistent and the set

A p := {(x, y) ∈ X : s + x = s in , x ≥ s in -λ -(D r )} ⊂ A
is (forward) attracting, with (where the function λ -is defined in Proposition 2.1).

Proof. Take ε > 0 such that µ(s) > D r for any s ∈ [s in -ε, s in + ε] and posit

η := min{µ(s) -D r : s ∈ [s in -ε, s in + ε]} > 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)-s in +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, s in ]. Consequently, there exists T (ω, ε) > 0 such that

s(t) < s in + ε, q(t) > -κ, ∀t > T (ω, ε), where κ := µ(s in -ε/2) -D r µ(s in -ε/2) (ε/2) > 0. If s(t) ∈ [s in -ε, s in + ε] for any t > T (ω, ε)
, then one has, from equation (1.8)), ẋ(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 (ω)) ≤ s in -ε/2. On another hand, from equation (1.7), s(t) can be written as the solution of the non-autonomous dynamics

ds dt = F (t, s) := (D + Φ(z * (θ t ω)) -µ(s))(s in -s) -µ(s)q(t). ( 4.5)
Note that one has

F (t, s in -ε/2) ≤ (D r -µ(s in -ε/2))(ε/2) -µ(s in -ε/2)q(t)
submitted to Mathematical Biosciences and Engineering and then

F (t, s in -ε/2) ≤ (D r -µ(s in -ε/2))(ε/2) + µ(s in -ε/2)κ = 0, t > T (ω, ε)
(from the definition of κ). We deduce that the set [0, s in -ε/2] is forward invariant for the semi-flow { ṡ = F (t, s), t > T (ω, ε)}. Therefore, one has

s(t) ≤ s in -ε/2, t > T (ω).
Then, one can write from equation (4.5)

ds dt ≤ (D r -µ(s))(s in -s) -µ(s)q(t), t > T (ω)
and from the comparison of solutions of scalar ordinary differential equations [START_REF] Walter | Ordinary Differential Equations[END_REF], one has the inequality s(t) ≤ s + (t) for any t > T (ω), where s + (t) is solution of the Cauchy problem

ds + dt = (D r -µ(s + ))(s in -s + ) -µ(s + )q(t), s + (T (ω)) = s(T (ω)).
Note that the solution s + (t) belongs to the interval [0, s in ] for any t ≥ T (ω) (and is thus bounded) and that its dynamics is asymptotic autonomous with limiting dynamics

ds † dt = (D r -µ(s † ))(s in -s † ) (4.6)
Under assumption D r < µ(s in ), one has necessarily λ -(D r ) < s in and the property

(D r -µ(σ))(λ -(D r ) -σ) > 0, ∀σ ∈ [0, s in ] \ {λ -(D r )}
is fulfilled. One finally obtains that the any solution of (4.6) in [0, s in ] is such that s † (t) → λ -(D r ) when t → +∞. From the theory of asymptotically autonomous dynamical systems [START_REF] Thieme | Convergence results and a poincare-bendixson trichotomy for asymptotically autonomous differential equations[END_REF], one concludes that s + (t) converges also to λ -(D r ) when t → +∞, which proves that one has lim sup

t→+∞ s(t) ≤ λ -(D r ) or equivalently lim inf t→+∞ x(t) ≥ s in -λ -(D r ) > 0
since s(t) + x(t) converges to s in for any realization. This demonstrates the strong persistence of the random dynamics (1.7)-(1.8) with the explicit lower bound s in -λ -(D r ). Now we consider the last situation when µ(s in ) < D < µ(s m ) which corresponds to the bistability 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. 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

dx dt = (µ(s) -(D + Φ(z * (θ t ω))))x = (µ(s in -x -q) -(D + Φ(z * (θ t ω))))x := F (t, x), ( 4.7)
where q = s in -x -s.

We know that lim t→+∞ q(t) = 0 then, for every ε > 0 there exists T (ω, ε) > 0 such that |q(t)| < ε for every t > T (ω, ε). Thus x(t) < s in + ε for every t > T (ω, ε).

Moreover, for all x ∈ (s in -λ -(D l )+ε, s in +ε), one has s in -x-q ∈ (-q -ε, λ -(D l )-q -ε). In fact, since s in -x -q = s ≥ 0 (see Theorem 3.1), then we have

s in -x -q ∈ (0, λ -(D l ) -q -ε) ⊂ (0, λ -(D l )) whence F (t, x) < 0 since µ(0, λ -(D l )) < D l and D + Φ(z * (θ t ω)) ≥ D l . Thus, lim sup t→+∞ x(t) ≤ s in -λ -(D l ) + ε
for every ε > 0.

Numerical simulations

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 a re p resented i n e ach c ase. T he fi rst fig ure (wi th two pan els) 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 o ffers i nformation a bout t he c onsumption f unction 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 s in = 14, µ 0 = 4, k s = 7.5, k i = 4, D = 1.4 and the interval given by practitioners from observations is [D l , D r ] = [1.15, 1.65], then d = 0.25. In addition, ν = 1 and = 1. We can observe in Figure 3 that D l > µ(s m ) holds true and then, as proved in Theorem 3.3, the species extinguishes (the attracting set in this case is {(0, s in )}). 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 D l < µ(s in ) < D is fulfilled, as in the deterministic case. However, the issue here depends also on the realization of the noise.

Conclusion

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 D l > µ(s m ) whatever is the input concentration s in . On the opposite, we proved in Theorem 4.1 the weakly uniformly persistence of species when D < µ(s in ), 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 < µ(s in ) ensures persistence in the deterministic case but the effective removal rate does not necessarily fulfills this condition depending on the realizations. Finally, we prove the strong persistence of the species in Theorem 4.2 under the stronger condition D r < µ(s in ). 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 depicts 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 µ(s in ) > D and µ(s in ) > D r be true or not. Condition to have weak persistence (µ(s in ) > 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 (µ(s in ) > D r ) 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.

Figure 1 .

 1 Figure 1. The chemostat

Definition 2 . 1 .Definition 2 . 2 .

 2122 The dynamics (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 inf t→+∞ x(t) > . The dynamics (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 sup t→+∞ x(t) > .

t 0 (

 0 D+Φ(z * (θsω)))ds .(3.2)

Theorem 3 . 3 .

 33 Assume that the inequality(3.3) D l > µ(s m )is fulfilled, w here s m i s d efined in(1 .4). Th en, th e singletonA e := {(s in , 0)} ⊂ A is a (forward) attracting set.Proof. From the equation describing the dynamics of the species (1.8) and (1.6) one has

Proposition 4 . 1 .

 41 Let beΦ(z) = -defined as in(1.5). Then Φ (z * (θ s ω))ds = 0, a.s. inΦ .(4.1)

Theorem 4 . 1 .

 41 Assume that the inequality µ(s in ) > D. (4.2)

Theorem 4 . 3 .

 43 Assume x 0 > 0. Then, the upper boundlim sup t→+∞ x(t; 0, ω, x 0 ) ≤ s in -λ -(D l )for any solution of the random dynamics (1.7)-(1.8) as long as D l < µ(s m ).

Figure 2 .

 2 Figure 2. Extinction of species. Values of the parameters: s in = 14, µ 0 = 4, k s = 7, k i = 5, D = 1.4, D l = 1.15, D r = 1.65, d = 0.25 , β = 1, ν = 1 and the initial values s 0 = 14, x 0 = 5

Figure 3 .

 3 Figure 3. Consumption function of consumer species in case of extinction. Values of the parameters: s in = 14, s m = 5.4772, D = 1.4, D l = 1.65, D l = 1.15, µ(s in ) = 0.7943, µ(s m ) = 10699.

Figure 4 .Figure 5 .

 45 Figure 4. Persistence of species in weak sense. Values of the parameters: s in = 16.7, µ 0 = 7, k s = 7, k i = 7, D = 1.7, D l = 1.45, D r = 1.95, d = 0.25, β = 1, ν = 1 and the initial values s 0 = 16.7, x 0 = 0.01

Figure 6 .Figure 7 .

 67 Figure 6. Persistence of species in strong sense. Values of the parameters: s in = 16.7, µ 0 = 7, k s = 7, k i = 7, D = 1.4, D l = 1.15, D r = 1.65, d = 0.25, β = 1, ν = 1 and the initial values s 0 = 16.7, x 0 = 0.01

Figure 8 .Figure 9 . 5 Figure 10 .

 89510 Figure 8. Extinction of species. Values of the parameters: s in = 20.1, µ 0 = 4, k s = 7, k i = 7, D = 1, D l = 0.75, D r = 1.25, d = 0.25, β = 1, ν = 1 and the initial values s 0 = 20.1, x 0 = 0.01

Acknowledgments

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).