Journal of Applied Mathematics

Volume 2013 (2013), Article ID 354095, 10 pages

http://dx.doi.org/10.1155/2013/354095

## Impulsive State Feedback Control of Cheese Whey Fermentation for Single-Cell Protein Production

^{1}Science College, Jimei University, Xiamen, Fujian 361021, China^{2}Academy of Mathematics and System Sciences, Chinese Academy of Sciences, Beijing 100080, China

Received 22 March 2013; Accepted 18 June 2013

Academic Editor: Yongkun Li

Copyright © 2013 Chunjin Wei 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

The work is the analysis of a mathematical model of cheese whey fermentation for single-cell protein production with impulsive state feedback control. Through the analysis, the sufficient conditions of existence and stability of positive order-1 periodic solution are obtained. It is shown that the system either tends to a stable state or has a periodic solution, which depends on the feedback state, the control parameter of the dilution rate, and the initial concentrate of microorganism and substrate. For some special cases, it is also shown that the system may exist in order-2 periodic solution. Furthermore, our findings are confirmed by means of numerical simulations.

#### 1. Introduction

Single-cell protein (SCP) is the protein extracted from cultivated microbial biomass. It can be used for protein supplementation of a staple diet by replacing costly conventional sources like soymeal and fishmeal to alleviate the problem of protein scarcity. Increasing demand for protein sources of high nutrition value have stimulated the application of single-cell proteins, yeast, bacteria, and algae in human foods or animal feeds [1]. The interest in agricultural and industrial wastes as substrates for the production of single-cell proteins has increased recently (see [2–5] and the references therein). For example, Kurbanoglu and Algur [6] studied single-cell protein production from ram horm hydrolysate by bacteria, and El-Saadany et al. [7] studied the production of single cell protein from agricultural wastes by fungi.

The production of microbial biomass is done mainly by a continuous fermentation system [8]. According to different reactions and differential control technologies, many dynamic models concerning the culture of microorganisms in the continuous bioreactor have been established [9–13]. However, there are a lot of factors affecting the growth and reproduction of the microorganisms in the process of bioreacts [14]. For example, for some aerobic microorganisms, the dissolved oxygen (DO) is a key factor to the growth of the aerobic microorganisms. In order to maintain the dissolved oxygen concentration in an appropriate range, it is necessary to regulate the microorganisms’ concentration to not exceed a set level [15].

As shown in Figure 1, if the microorganisms’ concentration reaches the preset value, the valve is opened by the photoelectric control devices, and a certain amount of regulator solution is input. Then, the microorganisms’ concentration is decreased and less than the preset value. At this time, the valve is closed by the control system until the next moment at which the microorganisms’ concentration reaches the preset value once again. Therefore, the microorganism concentration can be decreased through the control methods such as replenishing water or other regulator solution to reduce the concentration of the microorganism.

In recent years, in order to obtain the optimal conditions for microorganism growth and decrease the inhibition of the microorganism concentration or other negative effects, many authors introduced the state-dependent impulsion into the microorganism culture process and analyzed the system's dynamic behavior [16–18]. However, so far, no papers have discussed the mathematical model of cheese whey fermentation for single-cell protein production with impulsive feedback control.

This paper aims to propose a mathematical model of cheese whey fermentation for single-cell protein production with feedback control. We will discuss the existence and stability of periodic solution of the system with impulsive feedback control by using the existence criteria [19] and the stability theorem [20]. The paper is organized as follows. The mathematical model and some preliminary results are introduced in Section 2. In Section 3, the qualitative analysis of the system without impulsive effect is given. In Section 4, the existence and stability of order-1 periodic solution of the system with impulsive state feedback control are investigated. Numerical simulations and some discussion are provided in Section 5.

#### 2. Model Formulation

Ghaly et al. [8] proposed the following mathematical model of cheese whey fermentation for the production of single-cell protein: where denotes the concentration of the microorganism (g/L) and denotes the concentration of the limiting substrate (g/L) in the bioreactor at time . denotes the influent substrate concentration (g/L). is the dilution rate; is the maximum special growth rate; is the saturation constant, equal to the substrate concentration at one-half the maximum special growth rate (g/L); is the endogenous microbial decay coefficient (); is the yield coefficient; is the maintenance coefficient; and denotes the substrate uptake rate for the microorganism maintenance.

In some cases, if the concentration of microorganism is lower than a critical value (a preset value given by experiments), it is not necessary to take the control measures. But if the concentration reaches the critical value, some negative effects such as airborne inhibition and production inhibition and so on may happen, and the control methods such as liquefying medium and replenishing water should be taken to decrease the concentration of the microorganism. Once the concentration is decreased, it also should take some time to reach the critical value. Therefore, the system (1) can be modified as follows by introducing the impulsive state feedback control: where , , , is the fraction of the concentration of microorganism and substrate that decreases due to the feedback control when the microorganism concentrate reaches the critical value .

In the following, we mainly discuss the existence of periodic solution of (2) by the existence criteria of the general impulsive autonomous system. Before starting our main results, we give some topological definitions about impulsive differential equations.

*Definition 1 (Laksmikantham et al. [21]). *A triple is said to be a semidynamical system if is metric space, is the set of all nonnegative reals, and is a continuous function such that(i) for all ;(ii) for all and .

We denote sometimes a semidynamical system by .

For any , the function defined by is continuous, and we call the trajectory of . The set is called the positive orbit of . For any subset of , we let where is the attainable set of at . Finally, we set .

*Definition 2 (Laksmikantham et al. [21]). *An impulsive semidynamical system consists of a semidynamical system together with a nonempty closed subset of and a continuous function such that the following properties hold:(i)no point is a limit point of ,(ii) is a closed subset of .

We write and for any , . We call the impulsive set and the impulsive function.

We define a function as follows: here, we call the time without impulse of , that is is the first time when hits .

*Definition 3 (Laksmikantham et al. [21]). * Let be an impulsive semidynamical system, and let and . The trajectory of is a function defined on subset of ( may be ) to inductively as follows:
where is the sequence of impulse points of , which satisfied . is the sequence of time of impulses relative to , .

*Definition 4 (Laksmikantham et al. [21]). * A trajectory is said to be periodic of period and order if there exist positive integers and such that is the smallest integer for which and .

Theorem 5 (Brouwer's fixed-point theorem (Griffel [22])). * Every continuous mapping of a closed bounded convex set in into itself has a fixed point. *

Zeng et al. [19] have proved the existence criteria of periodic solution of an impulsive autonomous system by means of Browers’ fixed-point theorem. For convenience of reading, we repeat the main results as follows.

Consider the following general autonomous impulsive differential equations: Here, , , , , and are all functions mapping into , and is the set of impulse, and we assume(H1), are all continuous with respect to , in ;(H2) is a line, and and are linear functions of and . For each point, , and we define that :

Obviously, is also a line of or a subset of a line, and we assume that . From Definition 2, we know that (8) is an impulsive semidynamical system. The following theorem gives the conditions on which (8) has an order-1 periodic solution defined by Definition 4.

Theorem 6 (see [19]). * If system (8) satisfies assumptions (H1) and (H2), there exists a bounded closed simply connected region which has the following properties:*(i)*there is no singularity in it, and the boundary of satisfies ;*(ii)* cannot be tangent with trajectories of (8) except at end points and ;*(iii)*trajectories with initial point in will enter into interior of ; then, there must exist an order-1 periodic solution of system (8) in region . *

#### 3. Qualitative Analysis of System (1)

Before discussing the periodic solution of system (2), we should consider the qualitative characteristic of system (2) without impulsive effect, that is, consider the qualitative characteristic of system (1).

Lemma 7. *The system (1) is uniformly bounded. *

*Proof. *From the first equation of system (1), we have
The isocline , that is, , which intersects the -axis at points and , and , . Furthermore, , when from the right; therefore, there exists and such that . Let the straight line , and intersect at the points , , respectively. Define a function , where . The function intersects the line and at the points and , respectively. Then, we have , , . Hence, the system (1) is uniformly bounded. The proof is completed. We denote that the region is formed by four points , , , and .

The equilibrium points of system (1) satisfies the following equation. It can be seen that system (11) has a boundary equilibrium and a positive equilibrium if , where , .

In the following, we will analyze the stability of equilibrium and of system (1). If , the positive equilibrium point does not exist, and the boundary equilibrium is stable. That is, we have the following result.

Theorem 8. *The equilibrium is locally asymptotically stable if , where . Otherwise, it is unstable. *

*Proof. * The Jacobian matrix of system (1) at takes the form of

It is easy to obtain that the eigenvalues for the matrix are , . For , then , that is, ; hence, is locally asymptotically stable. Otherwise, is unstable. This completes the proof.

Theorem 9. *The equilibrium point is globally asymptotically stable if , where . *

*Proof. * Consider the Lyapunov function ; we can obtain the result immediately.

Theorem 10. * If , then the positive equilibrium point is locally asymptotically stable. *

*Proof. * The Jacobian matrix of system (1) at takes the form of
The eigenvalue problem for the provides the following characteristic equation:
where the coefficients , are
Note that ; , then, we have that is locally asymptotically stable.

In the following, we will discuss the global stability of system (1). Firstly, we will give the following lemma.

Lemma 11. *Suppose that is a periodic orbit with of system (1), and is the set which consists of all the points in Phase plane .**Denote
**
where , ; then, we can obtain .*

*Proof. *
for is a period function with , so ; hence,
It is evident that ; the proof is completed.

Theorem 12. *If , the positive equilibrium of system (1) is globally asymptotically stable. *

*Proof. *From Theorem 10, we know that is locally stable. According to Lemma 11, we can obtain if there exists periodic solution around ; then, it is stable for any periodic solution. This is impossible. According to the Poincare-Bendixson theorem, limit set of all orbits must be equilibrium point . This implies that is globally asymptotically stable in . This completes the proof.

#### 4. Existence and Stability of Periodic Solution of System (2)

##### 4.1. Existence of Order-1 Periodic Solution

From discussions of the qualitative characteristic of system (2) without the impulsive effects, we can see that if , is globally stable. For the initial points which satisfy and , if , then all the solutions of system (2) will tend to the positive equilibrium after at most finite times impulsive effects. Therefore, we mainly focus our attention on the cases , and .

According to the existence criteria (Theorem 6), our ideas to prove the existence of periodic solution is to construct a closed region such that all the solutions of (2) enter the closed region and retain there. The ideas will be illustrated as follows using Figures 2 and 3.

From Figure 2, we can see that the line interacts with the line at the point , where . The impulsive set , . The impulsive functions and map the impulsive set into , , where , , , and . From the third equation of (2), we know that for , and furthermore, . Since the straight lines and satisfy and , respectively, and both pass through the point , then it is clear that there are two cases of the point , that is, and .

Theorem 13. * Suppose that , , , then system (2) has an order-1 periodic solution; furthermore, this order-1 periodic solution is unique. *

*Proof. * Suppose that , then the trajectories of (2) starting from the region must interact with the segment . Next, we construct the closed region . Let be the intersection point of the line , and the trajectories of (2) starting from the point , let be the intersection point of the line , and , let be the intersection point of the line and the segment . From the qualitative characteristic of (1), it is easily obtained that , , , ; thus, the closed region consists of , , , , and (see Figure 2), where the arc is the part of trajectory passing through and . In addition, there is no singularity in the region and all the trajectories satisfying the conditions of the theorem enter the closed region and retain there; therefore, from Theorem 6, we know that system (2) has an order-1 periodic solution.

Suppose that . Similar to the discussions of the case , let be the intersection point of the line and the trajectories of (2) starting from the point ; we can take the arc as a part of the boundaries of the closed region . Then, we obtain the closed region which consists of , , , , , and (see Figure 3); therefore, from Theorem 6, we know that system (2) has an order-1 periodic solution. To sum up, the system (2) has an order-1 periodic solution under the condition of Theorem 13. This completes the proof.

*Remark 14. *If system (2) has an order-1 periodic solution, then the order-1 periodic solution is unique.

We can prove the uniqueness of order-1 periodic solution in system (2) by using the method of successor functions which was introduced in [23], here omitted.

##### 4.2. Stability of Order-1 Periodic Solution

In the following, we will analyze the stability of order-1 periodic solution in system (2). Firstly, we give one lemma to discuss the stability of this periodic solution of system (2).

Lemma 15. * The -periodic solution of the following system
**
is orbitally asymptotically stable if the Floquet multiplier satisfies the condition , where
**
with
**
and , , , , , , and are calculated at the point , , . is a sufficiently smooth function with , and is the time of the th jump. *

The proof of this lemma referres to Simeonov and Bainov [20].

In the following, we suppose that this periodic solution of system (2) with period passes through the points and (see Figures 2 and 3). As the expression and the period of this solution are unknown, we discuss the stability of this periodic solution by Lemma 15. In our case,

Then,

Set ; then;

Because is a periodic solution of system (2), we know that , by the proof of Lemma 11, that is, . Obviously, if . Therefore, we have the following theorem.

Theorem 16. * If , , and , then the order-1 periodic solution of system (2) is orbitally asymptotically stable. *

##### 4.3. Order-2 Periodic Solution

From Theorem 13, we know that system (2) has an order-1 periodic solution. In this subsection, we will discuss the existence of order-2 periodic solution.

Suppose that is a periodic solution of system (2), then and . It is easy to obtain that for by the third equation of system (2). Let be the arbitrary solution of system (2) starting from the point , and let the first interaction point of trajectory and the set be and let the corresponding interaction points of every impulse be , respectively. Therefore, under the effect of impulsive function , the corresponding points after impulse are , respectively.

By the qualitative analysis of system (2), if there exists the th impulse effect such that for all , then the sequence will be monotone, that is, or . In this case, system (2) has no order periodic solution. Therefore, if system (2) has order periodic solution, there must exist a such that . Besides, when , we can know that for . So, in the following, we mainly discuss the existence of order-2 periodic solution for .

If , then system (2) has an order-1 periodic solution.

If , then or . Suppose that , for and the qualitative properties of system (2), we know that and or . Therefore, there exists an order-2 periodic solution if holds.

If holds, we can obtain(i). While if holds, we have (ii). The same discussion for the case , that is, we have or . If holds, we have(iii). While if holds, we have(iv).

If the case of (ii) and (iv) holds, system (2) has an order-1 periodic solution. It follows from (i) (or (iii)) that , , which implies that system (2) exists with an order-2 periodic solution. Similar to the previous discussion and the proof of Proposition 3.2 in [24], we can obtain that there is no order periodic solution in system (2) for .

*Remark 17. *Since the exact expression of the solution to system (2) is unknown, it is difficult to give the sufficient conditions to guarantee that the system (2) has an order-2 periodic solution or not.

#### 5. Numerical Analysis and Discussion

In this paper, a mathematical model of cheese whey fermentation for single-cell protein production with impulsive state feedback control is proposed. We investigated the qualitative characteristic of the system without impulsive effect and show that the system is globally asymptotical stability. It is shown that the system with impulsive state feedback control has an order-1 periodic solution, and sufficient conditions for existence and stability of order-1 periodic solution are also obtained. The case, in which it is possible that there is an order-2 periodic solution, is also observed; furthermore, it is also pointed out that the system has no order periodic solution for . These results show that the cheese whey fermentation system for single-cell protein production with impulsive state feedback control either tends to a stable state or has a periodic solution.

To verify the theoretical results obtained in this paper, we give the numerical simulations of system (2). Let , , , , and ; then, we have , . If , , , and , then the system (2) tends to the positive equilibrium point, which can be seen in Figure 4. Figure 4 shows that no impulse occurs when . Figure 5 gives the time series and phase portrait when , , , and and show the trajectory tends to be periodic. Figure 6 is the time series and phase portrait when , , , and . Figures 5 and 6 show different positions of the periodic solution under different parameter value and different initial values. Furthermore, the phase portrait of Figure 5 indicates that the microorganism concentration always keeps increasing, but Figure 6 indicates that the microorganism concentration firstly decreases and then begins to increase. Therefore, a different control parameter can result in different bioprocesses and different efficiencies of microorganism culture. Which will give a conclusion theoretically to the researchers in the field of single-cell protein production. Researchers should give suitable control parameter and appropriate initial concentration in order to obtain a steady and optimal production.

In conclusion, we give a mathematical conclusion to the researchers in the field as follows: researchers of cheese whey fermentation for single-cell protein production should give suitable feedback state (the value of ), control parameters (), and appropriate initial concentration of microorganism and substrate in order to obtain a steady and optimal production.

#### Acknowledgments

This paper supported by National Natural Science Foundation of China (11171284, 31272653), Fujian Provincial Natural Science Foundation of China (2012J01012), the Fujian Provincial Education Foundation (JA12198), and the Scientific Research Foundation of Jimei University of China (ZC2011003).

#### References

- G. Emtiazi, Z. Etemadifar, and M. Tavassoli, “A novel nitrogen-fixing cellulytic bacterium associated with root of corn is a candidate for production of single cell protein,”
*Biomass and Bioenergy*, vol. 25, no. 4, pp. 423–426, 2003. View at Publisher · View at Google Scholar · View at Scopus - R. B. Vasey and K. A. Powell, “Single-cell protein,”
*Biotechnology and Genetic Engineering Reviews*, vol. 2, pp. 285–311, 1984. View at Publisher · View at Google Scholar · View at Scopus - R. D. MacElroy and D. Wang, “Waste recycling issues in bioregenerative life support,”
*Advances in Space Research*, vol. 9, no. 8, pp. 75–84, 1989. View at Google Scholar · View at Scopus - J. R. Mielenz, “Ethanol production from biomass: technology and commercialization status,”
*Current Opinion in Microbiology*, vol. 4, no. 3, pp. 324–329, 2001. View at Publisher · View at Google Scholar · View at Scopus - M. I. Rajoka, S. H. Khan, M. A. Jabbar, M. S. Awan, and A. S. Hashmi, “Kinetics of batch single cell protein production from rice polishings with
*Candida utilis*in continuously aerated tank reactors,”*Bioresource Technology*, vol. 97, no. 15, pp. 1934–1941, 2006. View at Publisher · View at Google Scholar · View at Scopus - E. B. Kurbanoglu and O. F. Algur, “Single-cell protein production from ram horn hydrolysate by bacteria,”
*Bioresource Technology*, vol. 85, no. 2, pp. 125–129, 2002. View at Publisher · View at Google Scholar · View at Scopus - R. El-Saadany, H. Khalaf, H. El-Manawaty, and F. Salom, “The production of single cell protein from agricultural wastes by fungi,”
*Acta Aliment Acad Sci Hung*, vol. 17, no. 4, pp. 376–377, 1988. View at Google Scholar - A. E. Ghaly, M. Kamal, and L. R. Correia, “Kinetic modelling of continuous submerged fermentation of cheese whey for single cell protein production,”
*Bioresource Technology*, vol. 96, no. 10, pp. 1143–1152, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. E. Ghaly and A. A. El-Taweel, “Kinetic modelling of continuous production of ethanol from cheese whey,”
*Biomass and Bioenergy*, vol. 12, no. 6, pp. 461–472, 1997. View at Publisher · View at Google Scholar · View at Scopus - K. Schügerl and K. H. Bellgardt, Eds.,
*Bioreaction Engineering: Modeling and Control*, Springer, Berlin, Germany, 2000. - S. Sun and L. Chen, “Dynamic behaviors of Monod type chemostat model with impulsive perturbation on the nutrient concentration,”
*Journal of Mathematical Chemistry*, vol. 42, no. 4, pp. 837–847, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. de Leenheer and H. Smith, “Feedback control for chemostat models,”
*Journal of Mathematical Biology*, vol. 46, no. 1, pp. 48–70, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Fu, W. Ma, and S. Ruan, “Qualitative analysis of a chemostat model with inhibitory exponential substrate uptake,”
*Chaos, Solitons & Fractals*, vol. 23, no. 3, pp. 873–886, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Sun, A. Kasperski, Y. Tian, and L. Chen, “Modelling and optimization of a continuous stirred tank reactor with feedback control and pulse feeding,”
*Chemical Engineering and Processing*, vol. 50, no. 7, pp. 675–686, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. Guo and L. Chen, “Periodic solution of a chemostat model with Monod growth rate and impulsive state feedback control,”
*Journal of Theoretical Biology*, vol. 260, no. 4, pp. 502–509, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - Z. Zhao, L. Yang, and L. Chen, “Impulsive state feedback control of the microorganism culture in a turbidostat,”
*Journal of Mathematical Chemistry*, vol. 47, no. 4, pp. 1224–1239, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Tian, A. Kasperski, K. Sun, and L. Chen, “Theoretical approach to modelling and analysis of the bioprocess with product inhibition and impulse effect,”
*BioSystems*, vol. 104, no. 2-3, pp. 77–86, 2011. View at Publisher · View at Google Scholar · View at Scopus - Z. Li, L. Chen, and Z. Liu, “Periodic solution of a chemostat model with variable yield and impulsive state feedback control,”
*Applied Mathematical Modelling*, vol. 36, no. 3, pp. 1255–1266, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Zeng, L. Chen, and L. Sun, “Existence of periodic solution of order one of planar impulsive autonomous system,”
*Journal of Computational and Applied Mathematics*, vol. 186, no. 2, pp. 466–481, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. E. Simeonov and D. D. Bainov, “Orbital stability of periodic solutions autonomous systems with impulse effect,”
*International Journal of Systems Science*, vol. 19, no. 12, pp. 2562–2585, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - V. Laksmikantham, D. D. Bainov, and P. S. Simeonov,
*Theory of Impulsive Differential Equations*, World Scientific, Singapore, 1989. - D. H. Griffel,
*Applied Functional Analysis*, Ellis Horwood, Chichester, UK, 1981. View at MathSciNet - B. Liu, Y. Tian, and B. Kang, “Dynamics on a Holling II predator-prey model with state-dependent impulsive control,”
*International Journal of Biomathematics*, vol. 5, no. 3, article 1260006, 18 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - G. Jiang, Q. Lu, and L. Qian, “Complex dynamics of a Holling type II prey-predator system with state feedback control,”
*Chaos, Solitons & Fractals*, vol. 31, no. 2, pp. 448–461, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet