Discrete Dynamics in Nature and Society

Volume 2015, Article ID 306302, 7 pages

http://dx.doi.org/10.1155/2015/306302

## Periodic Oscillations in a Chemostat Model with Two Discrete Delays

College of Science, University of Shanghai for Science and Technology, Shanghai 200093, China

Received 22 October 2014; Accepted 30 January 2015

Academic Editor: Daniele Fournier-Prunaret

Copyright © 2015 Tiansi Zhang 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.

#### Abstract

Periodic oscillations of solutions of a chemostat-type model in which a species feeds on a limiting nutrient are considered. The model incorporates two discrete delays representing the lag in nutrient recycling and nutrient conversion. Through the study of characteristic equation associated with the linearized system, a unique positive equilibrium is found and proved to be locally asymptotically stable under some conditions. Meanwhile, a Hopf bifurcation causing periodic solutions is also obtained. Numerical simulations illustrate the theoretical results.

#### 1. Introduction

The chemostat is very important to imitate ecosystem and is used extensively in microbiology and population biology for its measurable parameters and maneuverable experiments; see [1, 2], for example. Many models have been constructed based on chemostat-type equations such as the growth of plankton of unicellular algae and plasmid-bearing and plasmid-free competition; see [3–5] and the references therein. Amongst them, regeneration of nutrient decomposed by bacterial played an indispensable role for the stability and development of the nature and thus was taken into account; see [6–8], but no one treated it with a time delay. In fact, a delay in nutrient recycling on ecosystem always exists and is influenced by the change of temperature [9]. In view of this, Bischi [10] embedded this time delay and discussed its stabilizing effect for a chemostat model with oscillations. After that, He, Ruan and Wang et al. studied the global behaviors of a chemostat model with delayed nutrient recycling and periodically pulsed input in [11] and [12], respectively. Freedman and Xu [13] turned to consider a competition model with delayed nutrient recycling which generalized the problem in [14] and gave some persistence and extinction conditions for the competing populations. For some other excellent work about the delayed nutrient recycling, one may refer to [15–20].

On the other hand, environment is a vital factor that affects the growth of microorganisms. When it changes, a lag phase of the growth would appear. In 1969, when studying the growth of* Isochrysis galbana*, Caperon [21] used a time lag in his model with alterable environment and his research well met the corresponding experiment results. Besides it, Ellermeyer [22] and Xia et al. [23] also investigated the delayed response in growth and got global asymptotic behaviors of a competition and transient oscillations, respectively.

From the above we see that time delays have biological significance since they can alter dynamics of solutions of a system and display more transient oscillations especially in numerical simulation. So in this paper, we confine us to a single species model with a limited nutrient at a constant rate together with two delays in nutrient recycling and nutrient conversion, which takes the following form: with initial value conditions where denotes the limiting nutrient concentration and denotes the concentration of the microorganism at time . , the Banach space of all continuous bounded functions, and all parameters are positive constants. is the input concentration of the limiting nutrient, is the washout rate of the chemostat, is the death rate of the microorganisms, is the fraction of nutrient conversion, and is the maximal nutrient uptake rate for the microorganisms. There are two delays of and due to nutrient conversion and nutrient recycling, respectively. is the nutrient uptake function and has the following properties: This is the first time to study such a problem. The method used here is established by Ruan and Wei [24]. Our analysis will develop from discussing roots of the associated characteristic equation of the linearized system, whose coefficients depend heavily on delays.

The paper is organized as follows. In the next section, we investigate the existence of a positive equilibrium. Then in Section 3, we study the asymptotic stability of the equilibrium by first linearizing the system and then analyzing the roots of the associated characteristic equation, which is a second degree exponential polynomial. This process is divided into three cases of which two delays are zero, only one delay is zero, and no delays are zero. Finally a Hopf bifurcation causing the appearance of periodic solutions is found to exist. The last two sections are to numerically illustrate our findings and to summarize the discussion in the paper, respectively.

#### 2. The Boundedness and Existence of Nonnegative Equilibria of Model 1

Firstly, by studying model 1 with 3, we have the following result about the boundedness of solutions.

Theorem 1. *All solutions of system 1 with 3 are bounded.*

*Proof. *Define a Lyapunov function
Obviously, and as . Differentiate along trajectories of 1; that is,
It is easy to see that the coefficient of is negative since . So outside the region bounded by the positive orthant and the line
is negative. This completes the proof.

Secondly, note that is always an equilibrium for system 1, and from the following equations system 1 has a positive equilibrium given by provided that where means the inverse function of . Moreover, since the function is increasing, the equilibrium is unique. So we can conclude the following.

Theorem 2. *System 1 has a unique trivial equilibrium , and has a unique positive equilibrium defined by 8 if condition 9 holds.*

#### 3. Asymptotic Stability of

In this section, we are going to discuss the asymptotic stability of the positive equilibrium of system 1. Throughout this section, we assume that condition 9 holds. To begin we linearize 1 at and then construct the associated characteristic equation.

Setting and and putting them into system 1, we have Denote Then the characteristic equation associated with 10 is given by or equivalently by the second degree exponential polynomial in ,

From now we restrict our attention to discuss the local asymptotic stability of by analyzing the sign of the real parts of roots of the characteristic equation 13. Note that the coefficients of 13 depend heavily on the time delays and ; there is not a direct way to go, so by use of the method proposed by Ruan and Wei in [24], we consider three cases: ; ; and .

##### 3.1. The Case

Assume in 13 that ; then we have According to the Routh-Hurwith criterion, the real parts of all eigenvalues of 14 are negative if and only if 9 is valid, which implies the asymptotic stability of . Hence, we have the following result.

Theorem 3. *Assume that condition 9 and hold. Then all roots of 13 have negative real parts, and the positive equilibrium of system 1 is locally asymptotically stable.*

##### 3.2. The Case ,

Setting in 13, since , the characteristic equation is now where It is not hard to see that and .

We know from Theorem 3 that the is locally asymptotically stable if set in 15, while when increases from zero, the equilibrium could lose its stability if pure imaginary roots of 13 appear. Therefore it is convenient to look for pure imaginary roots of 15.

Put , into 15 and separate the real and the imaginary parts of the equation; then we have Taking the sum of squares in 17 yields It has real solutions provided Notice that if is a solution of 17, so is . Accordingly 15 has two pure imaginary solutions with .

Now differentiate both sides of 15 with respect to , and then Therefore, Clearly the sign of is positive at or is negative at .

Furthermore, we see that there are two sets of satisfying 17: where , , and where , and To sum up, we have the following.

Lemma 4. *Let be an element of the sequences , . For , the characteristic equation 15 has a pair of simple conjugate pure imaginary roots , satisfying
*

*From the above analysis, the stability of the positive equilibrium of system 1 for is obvious.*

*Theorem 5. Assume that conditions 9 and 20 hold and . Let
Then the equilibrium of system 1 is locally asymptotically stable when and a Hopf bifurcation occurs at when if and only if .*

*3.3. The Case *

*The case can be regarded as a general one of that in the last section, so our discussion here is based on the result in Section 3.2. The asymptotic stability of can be derived from Theorem 5. To state it, we first introduce a lemma about the sign of the real parts of roots of the characteristic equation 13.*

*Lemma 6. If all roots of 15 have negative real parts for , then there exists a , such that all roots of 13 have negative real parts for .*

*Proof. *Suppose that 15 has no roots with nonnegative real parts for ; that is, all the roots of 13 have negative real parts for and . According to Theorem 2.1 in [21], on the left side of 13, the sum of the multiplicity of zeros may change if a zero locates on or passes through the imaginary axis in the open right half-plane when varies.

Recall that 13 with has no root with nonnegative real part; then there must exist a such that all roots of 13 have negative real parts for .

*Theorem 7. Assume that conditions 9 and 20 hold. Then for any , where is defined in Theorem 5, there exists a such that the positive equilibrium of system 1 is locally asymptotically stable when .*

*Proof. *From Theorem 5, it follows that all roots of 15 have negative real parts when . Then Lemma 6 indicates the validity of the conclusion.

*4. Numerical Simulations*

*In this section, we give some numerical tests to support our findings in the previous section. Besides this, we observe the existence of periodic solutions arising from the Hopf bifurcation.*

*In system 1, we assume that the nutrient uptake function is of Michaelis-Menten type
where is the half-saturation constant [25, 26]. We take the parameters , , , , , and . Figure 1 shows the existence and the stability of for the case .*