Stochastic differential equations (SDEs) have been intensively used to model the natural phenomena in the last decades and these equations play a prominent applied role in various fields  . Most SDEs do not have explicit solution but these equations can be solved numerically to approximate their solutions . We use the SRK method to approximate the solution numerically (  , p. 69)    . Whereas, Rüemelin in  proposed the S-stage stochastic Runge-Kutta explicit method. Moreover, Xiao in  introduced the High strong order stochastic Runge-Kutta methods for stochastic differential equations in case of Stratonovich since our model is in Stratonovich sense.
This research deals with the stability of stochastic logistic model with Ornstein-Uhlenbeck process. This process is a batter model of Brownian motion (  , p. 86). The Ornstein-Uhlenbeck process is used to model stock price, currency exchange rates and velocity.
The classical Brownian motion was introduced by Scottish Botanist Robert Brown in (1827). He described this motion based on random movements of pollen grains in liquid or gas . This theory was further supported by Norbert Wiener (1923) who explained full mathematical theory of Brownian motion which existed as a rigorously defined mathematical objective to recognize his contribution . Brownian motion is a simple continuous stochastic process which is extensively used to model phenomena in various fields such as industry, dynamic process, physics, finance and fermentation process   . This Brownian motion is the cause of instability (turbulence) in fermentation process .
Fermentation process converts sugar into alcohol with the help of yeast . Madihah, in  used the logistic model to elaborate the cell growth of Clostridium acetobutylicum p. 262 in fermentation process proposed by Verhulst (1838). This model is inadequate to describe the cell growth of microorganism in fermentation process because fermentation process contains random fluctuations. Bahar and Mao in  introduced new stochastic logistic model for population dynamics. The stochastic logistic model was used in  to elucidate the cell growth of Clostridium acetobutylicum p. 262 in fermentation process. This cell growth of microorganism in fermentation process is affected by environmental noise, a batter model for modelling environmental noise is Ornstein-Uhlenbeck process. However, no specific research was conducted on stability stochastic logistic model with Ornstein-Uhlenbeck process for cell growth of microorganism in fermentation process. Nevertheless, there are some previous research works which have been done on stability of logistic equations with white noise, Lévy Jumps and other different aspects  - . None of these researches were investigated the stability of stochastic logistic model with Ornstein-Uhlenbeck process for cell growth of microorganism in fermentation process. This research establishes the sufficient conditions for stochastic logistic equation with Ornstein-Uhlenbeck process for positive equilibrium point and zero solution by using Lyapunov function. Since, the Ornstein-Uhlenbeck process is a better model of Brwonian motion (  , p. 86). In addition, our current study will indicate that the Ornstein-Uhlenbeck process is unfavorable for stability of cell growth in fermentation process. Moreover, we apply the SRK4 method to evaluate the numerical solution.
This paper is organized in five main sections; Introduction, Preliminaries and Models Description, Main Results, Numerical Simulation and Conclusion.
2. Preliminaries and Models Description
Throughout this paper; the notation is the expectation of , T shows the terminal time, t is time, illustrates the initial time, corresponds to the highest cell size, is initial cell size, denotes the maximum specific growth rate, illustrates a carrying capacity, indicates the random fluctuation, shows Brownian motion and b is a constant.
The simplest mathematical model to illustrate the exponential phase for cell growth in fermentation process is:
The solution of Equation (1) is:
where t is time, corresponds to the initial cell size and denotes the maximum specific growth rate. If , Equation (1) is strongly ascending, and if, , strongly descending. Hence, Equation (1) is not adequate to model stationary phase in fermentation. Therefore, Pierre Francois Verhulst (1838) introduced new model containing stationary phase in fermentation process. Thus, the exponential growth model (1) is augmented by the inclusion of multiplicative
factor of . Hence, the logistic ordinary differential equation is:
where is a carrying capacity of a microbial species and Equation (3) can be solved analytically by determining the solution is:
where is the highest cell size. Madihah, in  used Equation (3) to elaborate cell growth of Clostridium acetobutylicum p. 262 in batch fermentation process. Murray in  and May in  were proved that Equation (3) is stable. This model is inadequate to describe the cell growth of Clostridium acetobutylicum p. 262 in fermentation process because fermentation process suffers from random fluctuation. Bahar and Mao in  introduced new stochastic logistic equation to model uncontrolled fluctuation in stationary phase at fermentation process by using new perturbation very randomly which is:
where is diffusion coefficient. Model (3) with perturbation (5) is:
where is m-dimensional Weiner process and is corresponding to the random fluctuation. Bazli, in  used Equation (6) to evaluate the cell growth of Clostridium acetobutylicum p. 262 in fermentation process.
Remarks 2.1: Liu et al., in  studied the stability of stochastic logistic with white noise which cannot model random fluctuations in stationary phase. Therefore, this research introduces new stochastic logistic model to cope with this problem and established the sufficient condition for positive equilibrium point. The new model is Equation (6) with Ornstein-Uhlenbeck process. The Ornstein-Uhlenbeck process is:
where stands for positive Brownian motion at time t, shows the coefficient friction, indicates the diffusion coefficient and is white noise.
into Equation (7) hence, the Ornstein-Uhlenbeck process becomes;
Substituting Equation (9) into Equation (6) yields
Equation (10) is stochastic logistic model with Ornstein-Uhlenbeck process.
3. Main Results
First, it is needed to prove that Equation (10) has a unique positive solution and later on we will pay attention on stability.
Theorem 3.1: For all and for any Equation (10) has a unique positive solution.
Proof: For any given initial value , the coefficients of Equation (10) are locally Lipschitz continuous. Hence, there is a unique locally solution Where shows the explosion time Arnold in  and Friedman . To show that the is a unique positive solution. So, it is
really necessary to determine that . Let and
is sufficiently large for each component of , the stopping time for every integer is:
is sufficiently large means that the is increasing as . Hence, , whereas .
If we can show that a.s., then a.s. and a.s. for all . In other words, to complete the proof all we need to show is that a.s. For if this statement is false, then there is a pair of constants and such that
Therefore, there exist an integer
Define for all and the nonnegativity of this function can be seen from
Itô (  , p. 95)    (  , p. 54) formula shows that
To simplify A and B separately yields
An application of these facts and (  , p. 87) and somewhat lengthy calculation we have
By substituting the value of A and B into Equation (12) yields
If then , it follows the boundedness . Therefore, there is nonnegative number which independent of x and t.
If , obviously there is a positive number which is independent of x and t, and it follow bounded , such that
On the other hand, we show that there is a nonnegative number Q which independent of x and t such that
By substituting the inequality (19) into Equation (13) yields
where is Brownian motion, and taking integral and expectation
By taking into account, the inequality (11) then we have
For every , or . Therefore, is no less than either
Based on (20) we have
where illustrates the function . Letting n tends to infinity leads to dissidence
Therefore, we need to have . Eventually, our claim is proved truly.
In the next section we prove that Equation (10) with Ornstein-Uhlenbeck process is globally asymptotically stable in positive solution which illustrates the highly randomness in stationary phase.
Theorem 3.2: Equation (10) is globally asymptotically stable in positive equilibrium points , under following assumptions a.s.
where, and are positive constant then the solution satisfies the .
Proof: An application Itô’s formula (  , p. 95)    for Equation (10) and somewhat lengthy calculation leads to
For the convenience we set the denominator of Equation (22) to and according to the (  , p. 69), we have:
Taking limit from both sides of above equation yields:
is a Martingale with quadratic variation
Based on (  , p. 65) the
where then Equation (23) becomes
By substituting Equation (29) into Equation (22) we have:
Ultimately, prove is completed.
Remark 3.1: Liu et al., in  studied the stability of stochastic logistic model with white noise which illustrates highly randomness in lag phase and the stationary phase is steady means there is no fluctuations. In lag phase microbes are adapted with each other and as a result no growth takes place. But, the model, conditions and approach which are presented in this study quite different and easy respectively. This research considers Equation (6) with Ornstein-Uhlenbeck process which highly presents the randomness in stationary phase. In stationary phase microbes fight with each other to survive (uncontrolled fluctuations reached maximum level). In addition, concentration of toxin reaches a value which is unable to maintain the maximum cell growth of microorganism rate in fermentation process. These agents become the cause of instability at stationary phase in fermentation process. In conclusion, the stochastic logistic model which presented in Lui et al.,  is not suitable model for cell growth of microorganism in fermentation process. Since, Lui et al., in  narrates in zero solution (almost lag phase in fermentation process) the population (cell growth) is extinctive and in positive equilibrium point (almost stationary phase) the population will be permanent. It is contrary to the features phases in fermentation process. To support our theory, we use the SRK4 method for numerical simulations which presented in below section.
4. Numerical Simulation
In this section we consider a strong and accurate numerical method (SRK4) to elaborate the analytical results (  , p. 69)    . Worth mentioning, we cannot use the stochastic Runge-Kutta to approximate the numerical solution of Equation (10), due to it is in Itô sense . Thus, it is convertible to Stratonovich sense by using the below formula
where is used to denote the Stratonovich form of SDE (i.e. . Equation (10) and Equation (31) present some solution under different approach. We use the SRK4 method for numerical approximation. This method was introduced by Rümellin . SRK was developed based on the increment of Wiener
process, . Furthermore, Xiao in  introduced the High strong order
stochastic Runge-Kutta methods for stochastic differential equations in case of Stratonovich with scalar noise. Since our model is in Stratonovich sense. Figure 1 and Figure 2 show the stability of Equation (3) and Equation (10) in positive equilibrium for cell growth of microorganism in fermentation process respectively.
Figure 1. Shows the stability of Equation (10). For deterministic part we use , , , , , and and for stochastic part we chose same value just difference are , and respectively.
Figure 2. Illustrates the stability of Equation (10). For stochastic part we chose the , , , , , and . For deterministic part we use , , , , , .
The blue line indicates the simple path of Equation (10) and the red line shows the simple path of Equation (3) respectively. In Figure 1 and Figure 2, the difference are between the value of . The blue line indicates the simple path of Equation (10) and the red line shows the simple path of Equation (3) respectively. Figure 1 and Figure 2 are plotted under H1.
In both Figure 1 and Figure 2 lag phase is the luck of random fluctuations due to microbes are adopt with each other. Therefore, no growth takes place in this phase. Contrariwise, in stationary phase microbes fight with each other for food and space and the behavior of cell proliferation of microorganism is different in phases at fermentation process. As time evolves, the system illustrates the intrinsic variability of the competing within species and deviations from exponential growth arise. It happens as a result of the nutrient level and toxin concentration achieves a value which can no longer support the maximum growth rate. So, the stochastic fluctuations mainly affect the logistic growth of microorganism and for more see Remarks 3.1.
Remark 4.1: The red line shows the simple path of Equation (3) and the blue line indicates the simple path of Equation (10) respectively. If , there is no phases for cell growth which means this random fluctuations will destroy the phases in fermentation process on the other hands the process will tend to infinity.
Remark 4.2: Figures 1-4 show the stability of Equation (3) and Equation (10) in different equilibrium points. Figure 3 and Figure 4 are plotted under H1 option 2. The thick red line illustrates the stability of Equation (10) without Ornstein-Uhlenbeck process. It is thick in stationary phase due to the growth rate , is increased, compare the value of with three others Figures.
This research is conducted on stability of stochastic logistic model with Ornstein-Uhlenbeck process for cell growth of microorganism in fermentation.
Figure 3. Shows the stability of Equation (10). For deterministic part we use , , , , , and and for stochastic part we chose same value just difference are , and respectively.
Figure 4. Illustrates the stability of Equation (10). For stochastic part we chose the , , , , , , and . For deterministic part we use , , , , , .
This research proved that Equation (10) has a unique positive solution (see Theorem 3.1). Moreover, we proved that Equation (10) was stochastically stable in zero solution and positive equilibrium point (see Theorem 3.2 and Appendix A Theorem 1) and viewed in Figure 1 and Figure 2 and Appendix A. For numerical simulation we used SRK4 method to show the reality of this research. If the , the result will be infinity which means that there are no phases for cell growth. On the other hand, it will not preserve the stochastically stable. If , what will happen? Unfortunately, there are some impediments here which need some further investigation in future perspective.
This work is jointly supported by the National Natural Science Foundation of China under Grant Nos. 61573291, the Fundamental Research Funds for Central Universities XDJK2016B036.
Theorem 1: Equation (10) under the following assumption is globally asymptotically stable in zero solution a.s (almost surely).
(H2): For , then the .
Proof: An application of Itô’s formula for Equation (10) (  , p. 95) leads to
Equation (32) has two parts A and B.
To simplify the B we have
Using these facts and (  , p. 87) and somewhat lengthy calculation we have
By substituting the value of A and B into Equation (32) yields:
where is Brownian motion and taking expectation from both sides of above equation. We know the expectation of Brownian motion is zero (  , Section 3.2.3, p. 43), we get
For and based on (  , p. 66 and p. 69) the and . Therefore
Under H2 yields
Therefore, Equation (38) becomes
Ultimately, our claim is proved.
Figure 5. Shows the stability of Equation (10) with different value of parameters.
Figure 6. Illustrates the stability of Equation (10) with different value of parameters.
Note: The thick red line shows the stability of Equation (10) without Ornstein-Uhlenbeck process. It is thick due to the growth rate , is increased, compare the value of with Figure 5.
 Ayoubi, T., Rosli, N., Bahar, A. and Salleh, M.M. (2015) Time Delay and Noise Explaining the Behaviour of the Cell Growth in Fermentation Process. AIP Conference Proceedings, 1643, 547-554.
 Burrage, K. and Burrage, P.M. (1996) High Strong Order Explicit Runge-Kutta Methods for Stochastic Ordinary Differential Equations. Applied Numerical Mathematics, 22, 81-101.
 Ariffin, N.N.A., Rosli, N. and Kasim, A.M. (2017) Stability Analysis of 4-Stage Stochastic Runge-Kutta Method (SRK4) and Specific Stochastic Runge-Kutta Method (SRKS1.5) for Stochastic Differential Equations. In: Kor, L.K., Ahmad, A.R., Idrus, Z. and Mansor, K., Eds., Proceedings of the 3rd International Conference on Computing, Mathematics and Statistics, Springer, Singapore, 187-194.
 Xiao, A. and Tang, X. (2016) High Strong Order Stochastic Runge-Kutta Methods for Stratonovich Stochastic Differential Equations with Scalar Noise. Numerical Algorithms, 72, 259-296.
 Mishura, I.S. and Mishura, Y. (2008) Stochastic Calculus for Fractional Brownian Motion and Related Processes. Springer Science and Business Media, Berlin, Heidelberg.
 Madihah, M.S., Ariff, A.B., Khalil, M.S., Suraini, A.A. and Karim, M.I.A. (2001) Anaerobic Fermentation of Gelatinized Sago Starch-Derived Sugars to Acetone—1-Butanol—Ethanol Solvent by Clostridium acetobutylicum. Folia Microbiologica, 46, 197-204.
 Madihah, M.S. (2002) Direct Fermentation of Gelatinised Sago Starch to Solvent (Actone, Butanol and Ethanol) by Clostridium acetobutylicum. Ph.D. Thesis, Universiti Putra Malaysia, Serdang, Selangor, Malaysia.
 Arugaslan, D. and Güzel, L. (2015) Stability of the Logistic Population Model with Generalized Piecewise Constant Delays. Advances in Difference Equations, 1, 173.
 Jiang, D., Shi, N. and Li, X. (2008) Global Stability and Stochastic Permanence of a Non-Autonomous Logistic Equation with Random Perturbation. Journal of Mathematical Analysis and Applications, 340, 588-597.
 Liu, Q. and Chen, Q. (2016) Analysis of a General Stochastic Non-Autonomous Logistic Model with Delays and Lévy Jumps. Journal of Mathematical Analysis and Applications, 433, 95-120.
 Liu, Q., Jiang, D., Hayat, T. and Alsaedi, A. (2018) Long-Time Behavior of a Stochastic Logistic Equation with Distributed Delay and Nonlinear Perturbation. Physic A: Statistical Mechanics and Its Applications, 508, 289-304.
 Mateos-Núñez, D. and Cortés, J. (2013) Stability of Stochastic Differential Equations with Additive Persistent Noise. 2013 American Control Conference, Washington DC, 17-19 June 2013, 5427-5432.
 Rosli, N., Bahar, A., Su Hoe, Y. and Abd Rahman, H. (2010) Stochastic Model of Gelatinised Sago Starch to Solvent Production by C. Acetobutylicum P262. In: Proceedings of the Regional Conference on Statistical Sciences 2010, Shah Alam, Selangor, 9-20.