Global Optimization of Continuous Fermentation Involving Zymomonas mobilis

Global steady-state optimization of Zymomonas mobilis fermentation process problems is performed to demonstrate the existence of multiple optimum solutions necessitating the use of a global optimization strategy. It is shown that the steady-state equations for the Zymomonas mobilis fermentation process can be reduced to a single equation which when used as a constraint results in yielding only one optimum solution, the global one.


Introduction
Production of ethanol from a variety of biomass feedstocks, solid and liquid waste resources is a national priority for a variety of reasons but most importantly for energy independence and sustainability.The biggest hurdle in the commercialization of any of the processes for making ethanol is the cost of production.Hence, it is important to develop strategies to produce clean ethanol with the least amount of expense.Ethanol can be produced by the fermentation of Zymomonas mobilis and it is important to ensure that that the process of fermentation operates in an optimal fashion.
A lot of work has been done in dealing with the nonlinearity of the Zymomonas Mibilis fermentation process [1][2][3][4][5][6].The most comprehensive work that demonstrates the nonlinear behavior of the Zymomonas mobilis fermentation process is done by [7][8][9].These articles show that the model presented by [7] is excellent not only in predicting the output concentrations but also the nonlinear behavior that results in multiple steady-states and oscillatory behavior.
While almost all of the computational work regarding the fermentation of Zymomonas mobilis involves demonstration of bifurcations and multiple steady-states, there has been little work concerning global optimization of these processes.
This paper deals with 1) a computational demonstration of the existence of multiple optimum solutions for the fermentation of Zymomonas mobilis requiring the use of a global optimum strategy to obtain the most optimum solution, and 2) using the process of elimination to reduce the steady-state equations for the Zymomonas mo-bilis fermentation process to a single equation and 2) demonstrating that if this equation is used instead of the original set of equations only one minimum (the global one) is obtained and not the unnecessary misleading local minima.This paper is organized as follows.First the model equations describing the steady-state fermentation process involving Zymomonas mobilis is described.In the next section a process of elimination is used to reduce these equations to a single equation.The optimization of the Zymomonas mobilis fermentation process using both these formulations are performed and compared.It is demonstrated that the original model optimization yields both local and global optimum solution, while the reduced model optimization yields only the global optimum.

Zymomobilis mobilis Fermentation
The steady-state equations governing the Zymomonas mobils fermentation process are given by [7] as 1 0 These equations can be succinctly represented as The variables are C X , C e , C S & C P , which are the concentrations of the biomass, key component, substrate and product.Table 1 has the values of the parameters and constants.The product concentration C P must be maximized which implies that the optimization minimization objective function is M P f C   min subject to 0 and the equality constraints are given by Equations ( 1)-( 4).Hence the optimization problem can be expressed as Garhyan and Elnshaie [7] studied the fermentation for C S0 values of 150.3 and 200 and performed a detailed bifurcation analysis for these two values of substrate concentrations.In this paper the MATLAB global optimization routines Global-Search and Multistart routines were used to determine all the optimum solutions satisfying Equation ( 6) for the C S0 values from 150 to 200.More one minima were found for each of these C S0 values.This demonstrates that to obtain the optimum solution for the Zymomonas mobilis fermentation process, it is necessary to use a global strategy if the equations are used in the form given by ( 1)-( 4).Table 2 gives the optimum values of C P for the various values of the substrate concentration C S0 .The global optimum solution is indicated as "g".
The entire range of C S0 values for which the Zymomonas mobilis fermentation process has been tested theoretically and experimentally [7,8] has been considered and two optimum solutions were found for the entire range when Equations ( 1)-(4) were used as constraints.

Reformulation of the Zymomonas mobilis
Fermentation Steady-State Equations Similar substitution for Equating ( 9) and (10) yields Simplification yields The optimization problem can be written as reduc min subje (15)

Global Optimization of the Zymomonas Fermentation Process
The global optimization of the Zymomonas mobilis fermentation process is performed using the global search and multistart routines in the global optimization MAT-LAB toolbox.Both these solvers use gradient based methods and start a local solver (fmincon) from multiple starting points to obtain local and global minima.The global search optimization uses a scatter search algorithm to generate multiple starting points while the multistart algorithm uses uniformly distributed starting points with predefined bounds.Both these solvers are used to ensure the correctness of the local and global optimum solutions.
Table 2 provides the optimum values of the product concentration that were obtained for each of the input substrate concentrations.The whole range of input concentrations for which this model is theoretically studied and experimentally verified [7,8] is investigated for the existence of multiple optimum values of C P .It is shown that the original formulation given by Equation ( 6) yields 2 optimum values of C P .The global value is labeled (g).In all the cases when the formulation given by Equation ( 15) is used only the global value is obtained.Even when the local value is used as a starting point the MATLAB routine fmincon converges to the global value when the reduced formulation is used.This concludes that the reduced formulation (given by Equation ( 15)) yields only one optimum solution (the global one).This reduced formulation does not admit the unnecessary and misleading global solution that the original formulation admits.

Conclusions
The steady-state optimization of continuous fermentation processes involving both Zymomonas mobilis is shown to have multiple solutions satisfying the mathematical requirements for optimization.This makes it necessary to use global optimization strategies to obtain the most optimum solution for the optimization of the fermentation process.The steady-state equations for the Zymomonas mobilis fermentation process can be reduced to a single equation which when used as a constraint, yields only the required global optimum solution and if this new formulation is used, even a local optimization strategy will produce the most optimum solution.