Optimal Management during the Microorganism Culture Based on the Continuous Purifying Effort

This paper deals with the problem of selective harvesting in a chemostat model. Here, we have taken the purifying effort as a dynamic variable and tax as a control instrument. The existence of the possible steady states along with their globally stable equilibrium is discussed. The optimal tax policy is also discussed with the help of Pontryagin’s maximum principle. Finally, numerical examples are taken to illustrate some of the key results.


Introduction
The chemostat is an important laboratory apparatus used to culture microorganisms 1-4 .Species grow in continuously stirred homogenous fermenters which are fed continuously by a nutrient and the cells are drawn off continuously.Such models have applications in ecology to model biological behavior of a simple lake and in biotechnology to model bioreaction in commercial bioreactors.Predictions based on parameters in the model that can be measured have been tested experimentally and outcomes have shown to agree rather well with the theory.Therefore, it has been extensively used in agriculture and many industrial applications i.e., pharmaceuticals, nutraceuticals, hydrogen production, and waste treatment .
It is well known that extracting microorganisms has a strong impact on the dynamic behavior of chemostat, and microorganism resources in the chemostat are usually harvested with the purpose of achieving the economic interest.Achieving a high productivity in a bioreactor plays a crucial role in determining the economics of bulk biochemical products such as ethanol.For the purpose of continuously culturing the microorganism and reaching the maximum profit, it is necessary to establish a constructive management of commercial extraction of the microorganism resources.The techniques and issues associated with the bioeconomic exploitation have been discussed in detail by Clark 5 .Taxation and market price usually are considered as possible factors affecting the producers' profit.However, because of the economic flexibility of the Taxation 5 , economists are particularly attracted to taxation because a competitive system can be better maintained under taxation rather than other regulatory methods.
Recently, there has been a considerable interest in the modeling of harvesting of biological resources 6-9 .The harvest effort is considered to be a dynamic variable; several kinds of harvesting policies are utilized to study the dynamical behavior of the model system.Furthermore, the optimal harvesting dynamical behavior of a harvested prey-predator model policies with taxation is also discussed.However, to the authors' best knowledge, there is no attempt focusing on discussing the optimal tax policy in the chemostat.In this paper, an extracted chemostat model is established in Section 2. The stability analysis of the equilibria is performed in Section 3. Furthermore, an optimal extracting policy for microorganism is also discussed in Section 4.

Model Formulation and Stability of the Equilibria
The general model of continuously culturing microorganism in a chemostat is given by the following form of differential equations 10 : where S t denotes the concentration of the substrate and x t denotes the concentration of the microorganism in the chemostat at time t.D is the dilution rate of the chemostat.S 0 is the concentration of the input substrate.The constant δ is the yield constant.μ is called the maximal specific growth rate of the microorganisms.K s is the self-saturation constant.In practice, when the microorganism has been cultured, it should be purified from the vessel.Hence, we suppose E t represents the harvest or extraction effort in order to purify the microorganism.To conserve the resource, the regulatory agency imposes a tax τ > 0 per unit biomass τ < 0 denotes the subsidies given to the harvest or extraction effort .Based on the above aspects, the model can be governed by the following differential equations: where p is the price of the unit harvest effort and c is the cost of the unit effort.Other parameters are the same as system 2.1 .By simple computation, we obtain that system 2.2 has a trivial equilibrium P 0 S 0 , 0, 0 and Owing to ∂S * /∂τ < 0, there exists a maximum value τ max such that 0 ≤ τ ≤ τ max .It provides the range of tax for the existence of interior equilibrium, which is of inspiration for people to regulate the extraction effort by means of economic instrument.Furthermore, it is easy to show the positivity and boundedness of solutions of the model system 2.2 .
Remark 2.1.It is biologically meaningful to interpret the positivity and boundedness of solutions of the model system 2.2 .Since the component S t , x t of the solution of system 2.2 represents the relation between the microorganism and substrate, the positivity of solutions reflects the survival of microorganism in the cultured vessel and the boundedness of solutions reveals a natural restriction to growth as a consequence of limited resources.Furthermore, due to the limitation of the microorganism, the extraction effort cannot increase without any restriction.
Next, we begin to analyze the stability of the equilibria.The characteristic equation of equilibrium P 0 is The stability of the equilibrium P 0 S 0 , 0, 0 is determined by λ μS 0 / K S S 0 −D.Therefore, we have the following theorem.Theorem 2.2.If μS 0 / K S S 0 < D holds, then the microorganism-free equilibrium The characteristic equation of equilibrium P 1 is The stability of the equilibrium P 1 is determined by and λ k p −τ x * −c .According to the relation between roots and coefficients, all the eigenvalues of system 2.5 has negative real parts.If Proof.The local stability has been proved above.Next, we prove the attractivity.Construct a Lyapunov function where the nonnegative constants c i i 1, 2, 3 will be determined later.
We compute the derivative V S, x, E along the system 2.2 Hence, we obtain the equilibrium P 1 which is globally asymptotically stable for p − τ δS 0 < p − τ δK S D/ μ − D c.
The proof is completed.
Next, we consider the stability of the positive equilibrium.Proof.The characteristic equation of equilibrium P 2 is where According to the Routh-Hurwitz criterion 11 , all the roots of the characteristic equation have negative real parts for the above cubic equation and the following criteria a 1 > 0, a 1 a 2 − a 3 > 0 should be satisfied.Hence, the equilibrium P 2 S * , x * , E * is locally asymptotically stable.
Define a function where the nonnegative constants c i i 1, 2, 3 are positive.We compute the derivative V S, x, E along the system 2.2

Optimal Extraction Policy during the Bioprocess
Fermentation technology is a response for producing the majority of bioproducts.Since some substrates of bioprocess are expensive, it is important to optimize the process to maximize the desired products and profits.The objective of the regulatory agency is to maximize the total discounted net revenues that the factory derives from the microorganism fermentation.Symbolically, this objective amounts to maximizing the present value J of a continuous timestream of revenues given by where δ denotes the instantaneous annual rate of discount.c is the extraction cost per unit effort.p is the price per unit biomass of x.To solve this optimization problem, we utilize Pontryagin's maximal principle 12 .
We treat τ as the control variable and wish to determine a suitable effort which maximizes J subject to the system 2.2 and the control constraints 3.2 τ min and τ max represent a feasible upper and lower limit of tax for the harvest effort, respectively.Specially τ min < 0 implies that subsidies have the effect of increasing the rate of expansion of the extraction.The Hamiltonian function is given by where λ i t i 1, 2, 3 are additional unknown functions called the adjoint variables.The Hamiltonian 3.3 must be maximized for τ ∈ τ min , τ max .Assuming that the control constraints are not binding i.e., the optimal solution does not occur at τ τ min or τ τ max .we have a singular control given by ∂H/∂τ 0. Hence, we can obtain λ 3 0.
The adjoint equations are We obtain from 3.6 that In order to obtain an optimal equilibrium solution, by considering the interior equilibrium P 2 , 3.4 can be rewritten as 3.8 we can obtain A e −δt , 3.9 where A D K S μx * / δ K S S * 2 .Substituting 3.7 and 3.9 into 3.5 , we have which provides an equation to the singular path and gives the optimal equilibrium levels of microorganism S * S * δ , x * x * δ .Then the optimal equilibrium levels of the harvest effort and tax can be obtained as follows:

3.11
Next, with the help of MATLAB, a simulation work with a hypothetical set of parameters is performed to understand the theoretical results which have been established.Let the parameters be S 0 5, K S 0.2, μ 1, p 3, D 0.6, δ 0.6, and c 2. For the system 2.2 , the range of the taxation can be obtained τ ∈ 0, 2.267399 in view of the positive equilibrium P 2 S * , x * , E * .According to the parameters given above, 3.10 can be numerically computed as follows:

Discussion
In this paper, a bioeconomic model is established to investigate the effects of the harvest effort on the dynamic behavior of the chemostat.In Theorem 2.2, we obtain the microorganism-free equilibrium P 0 S 0 , 0, 0 is stable if μS 0 / K S S 0 < D holds, which is simulated in Figure 2. Theorem 2.3 shows that the equilibrium P 1 is globally asymptotically stable, if p − τ δS 0 < p − τ δK S D / μ − D c holds see Figure 3 .The existence and global stability of the positive equilibrium is proven in Theorem 2.4.
Nowadays, the biological resources in the chemostat model are mostly harvested with the aim of achieving economic interest and the taxation is used as an economic control instrument to protect the resources from overexploitation, which motivates the introduction of the harvest effort and tax into the proposed model.The application of the control theory enabled us to show the existence of a unique optimal equilibrium point which is stable.These results can be used as a microorganism culture such as ethanol fermentation and lactic acid fermentation to obtain a more economic profit.

Figure 1 :
Figure 1: The positive equilibrium of the optimal purification.a Time series of the substrate concentration.b Time series of the microorganism concentration.c Time series of the purification effort.d Phase space trajectories corresponding to the optimal τ 0.2443807704.

Figure 2 :
Figure 2: The dynamic behavior of the microorganism-free equilibrium with the parameters D 0.4, S 0 0.3, δ 1, K s 0.9, c 1, k 0.1, and p 0.3.a Time series of the substrate concentration.b Time series of the microorganism extinction.c Time series of the effort loss.d Phase space trajectories of the microorganism-free equilibrium.

Figure 3 :
Figure 3: The dynamic behavior of the effort-free equilibrium with the parameters D 0.4, S 0 2, δ 1, K s 0.5, c 3, k 0.3, and p 3. a Time series of the substrate concentration.b Time series of the microorganism concentration.c Time series of the effort loss.d Phase space trajectories of the effort-free equilibrium.