_{1}

^{*}

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.

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-6]. The most comprehensive work that demonstrates the nonlinear behavior of the Zymomonas mobilis fermentation process is done by [7-9]. These articles show that the model presented by [

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 mobilis 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.

The steady-state equations governing the Zymomonas mobils fermentation process are given by [

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. _{P} must be maximized which implies that the optimization minimization objective function is and the equality constraints are given by Equations (1)-(4). Hence the optimization problem can be expressed as

Garhyan and Elnshaie [_{S}_{0} 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_{S}_{0} values from 150 to 200. More than one minima were found for each of these C_{S}_{0} 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). _{P} for the various values of the substrate concentration C_{S}_{0}. The global optimum solution is indicated as “g”.

The entire range of C_{S}_{0} 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.

Let

This implies that

and

Substituting for from Equation (2) into Equation (3) we take

_{P} Values for Zymomonas mobilis fermentation which are obtained by using the original model in [

and substituting for C_{S} from Equation (6) yields

Similar substitution for from Equation (2) into Equation (4) yields

Equating (9) and (10) yields

Simplification yields

Defining we get

The optimization problem can be written as

The global optimization of the Zymomonas mobilis fermentation process is performed using the global search and multistart routines in the global optimization MATLAB 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. _{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.

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.

C Concentration D Dilution Rate V Volume E ethanol X cell mass concentration in K_{1}, K_{2}, K_{3 } empirical constants P Permeability Y yield factor

e key component of biomas 0 input P product S substrate X biomass