Optimal switching control of a fed-batch fermentation process

Considering the hybrid nature in fed-batch culture of glycerol biconversion to 1,3-propanediol (1,3-PD) by Klebsiella pneumoniae, we propose a state-based switching dynamical system to describe the fermentation process. To maximize the concentration of 1,3-PD at the terminal time, an optimal switching control model subject to our proposed switching system and constraints of continuous state inequality and control function is presented. Because the number of the switchings is not known a priori, we reformulate the above optimal control problem as a two-level optimization problem. An optimization algorithm is developed to seek the optimal solution on the basis of a heuristic approach and control parametrization technique. Numerical results show that, by employing the obtained optimal control strategy, 1,3-PD concentration at the terminal time can be increased considerably.     

Discrete optimal control model and bound error for microbial continuous fermentation

A discrete optimal control model is given according to the discrete nonlinear dynamic system of continuous fermentations of glycerol to 1,3-propanediol (1,3-PD). The property of some major functions and their bounds on approximation errors are studied in this paper. Then, the conclusion that the optimality function of discrete optimal control model is the consistent approximation to one of the continuous optimal control model is proved. The results presented in this work can be used as guidelines for the optimal algorithm and its convergence.     

Optimal switching control for microbial fed-batch culture

In fed-batch culture of glycerol bio-dissimilation to 1,3-propanediol (1,3-PD), the aim of adding glycerol is to obtain as much 1,3-PD as possible. Hence, a proper feed strategy is required during the process. In this paper, we present an optimal switching control model based on our proposed controlled switching system. Some properties of the controlled switching system are obtained. Subsequently, we prove the existence of optimal control. In order to deduce the optimality conditions, we transcribe the optimal switching control model into an equivalent one with fixed switching instants and parameters. Finally, the optimality conditions of the equivalent problem are investigated by calculus of variations.     

Parameter identification and optimization algorithm in microbial continuous culture

In this study, a novelty mathematical model is established to formulate the continuous culture of glycerol to 1,3-Propanediol (1,3-PD) by Klebsiella pneumoniae, in which the inhibition of 3-hydroxypropionaldehyde (3-HPA) to cells growth and activity of some enzymes (such as glycerol dehydratase (GDHt) and 1,3-PD oxidoreductase (PDOR)), and the passive diffusion and active transport of glycerol and 1,3-PD across cell membrane are all taken into consideration. Taking the mean relative error between the experimental data and calculated values as the performance index, a parameter identification model involving multiple nonlinear dynamic systems is presented. The identifiability of the parameter identification model is also proved. Finally, an improved particle swarm optimization (PSO) algorithm is constructed to find the optimal parameters for the systems under substrate limitation and excess conditions, respectively. Numerical results not only show that the established model can be used to describe the continuous fermentation reasonably, but also the improved PSO algorithm is valid.     

An Exact Penalty Function Method for Continuous Inequality Constrained Optimal Control Problem

In this paper, we consider a class of optimal control problems subject to equality terminal state constraints and continuous state and control inequality constraints. By using the control parametrization technique and a time scaling transformation, the constrained optimal control problem is approximated by a sequence of optimal parameter selection problems with equality terminal state constraints and continuous state inequality constraints. Each of these constrained optimal parameter selection problems can be regarded as an optimization problem subject to equality constraints and continuous inequality constraints. On this basis, an exact penalty function method is used to devise a computational method to solve these optimization problems with equality constraints and continuous inequality constraints. The main idea is to augment the exact penalty function constructed from the equality constraints and continuous inequality constraints to the objective function, forming a new one. This gives rise to a sequence of unconstrained optimization problems. It is shown that, for sufficiently large penalty parameter value, any local minimizer of the unconstrained optimization problem is a local minimizer of the optimization problem with equality constraints and continuous inequality constraints. The convergent properties of the optimal parameter selection problems with equality constraints and continuous inequality constraints to the original optimal control problem are also discussed. For illustration, three examples are solved showing the effectiveness and applicability of the approach proposed.     

Pathway identification using parallel optimization for a complex metabolic system in microbial continuous culture

The bio-dissimilation of glycerol to 1,3-propanediol (1,3-PD) by Klebsiella pneumoniae (K. pneumoniae) is a complex bioprocess due to the multiple inhibitions of substrate and products onto the cell growth. In consideration of the fact that both the inhibition mechanisms of 3-hydroxypropionaldehyde (3-HPA) onto the cell growth and the transport systems of glycerol and 1,3-PD across the cell membrane are still unclear, we consider 72 possible metabolic pathways, and establish a novel mathematical model which is represented by an eight-dimensional nonlinear dynamical system. The existence, uniqueness, continuous dependence of solutions to the system and the compactness of the solution set are explored. On the basis of biological robustness, we give a quantitative definition of robustness index of the intracellular substances. Taking the robustness index of the intracellular substances together with the relative error between the experimental data and the computational values of the extracellular substances as a performance index, a parameter identification model is proposed for the nonlinear dynamical system, in which 43848 continuous variables and 1152 discrete variables are involved. A parallel particle swarm optimization — pathways identification algorithm (PPSO-PIA) is constructed to find the optimal pathway and parameters under various experiments conditions. Numerical results show that the optimal pathway and the corresponding dynamical system can describe the continuous fermentation reasonably.     

A computational method for free terminal time optimal control problem governed by nonlinear time delayed systems

This paper considers a free terminal time optimal control problem governed by nonlinear time delayed system, where both the terminal time and the control are required to be determined such that a cost function is minimized subject to continuous inequality state constraints. To solve this free terminal time optimal control problem, the control parameterization technique is applied to approximate the control function as a piecewise constant control function, where both the heights and the switching times are regarded as decision variables. In this way, the free terminal time optimal control problem is approximated as a sequence of optimal parameter selection problems governed by nonlinear time delayed systems, each of which can be viewed as a nonlinear optimization problem. Then, a fully informed particle swarm optimization method is adopted to solve the approximate problem. Finally, two free terminal time optimal control problems, including an optimal fishery control problem, are solved by using the proposed method so as to demonstrate its applicability.     

Optimal investment and premium control in a nonlinear diffusion model

This paper considers the optimal investment and premium control problem in a diffusion approximation to a non-homogeneous compound Poisson process. In the nonlinear diffusion model, it is assumed that there is an unspecified monotone function describing the relationship between the safety loading of premium and the time-varying claim arrival rate. Hence, in addition to the investment control, the premium rate can be served as a control variable in the optimization problem. Specifically, the problem is investigated in two cases: (i) maximizing the expected utility of terminal wealth, and (ii) minimizing the probability of ruin respectively. In both cases, some properties of the value functions are derived, and closed-form expressions for the optimal policies and the value functions are obtained. The results show that the optimal investment policy and the optimal premium control policy are dependent on each other. Most interestingly, as an example, we show that the nonlinear diffusion model reduces to a diffusion model with a quadratic drift coefficient when the function associated with the premium rate and the claim arrival rate takes a special form. This example shows that the model of study represents a class of nonlinear stochastic control risk model.     

Bilevel parameters identification for the multi-stage nonlinear impulsive system in microorganisms fed-batch cultures

Considering the abrupt jump of the substrate and different characters of the bacterial in the lag phase, the exponential phase and the stationary phase this paper proposes the multi-stage nonlinear impulsive system for the fed-batch fermentation from glycerol to 1,3-propanediol (1,3-PD) and establishes the bilevel identification system for its sensitive parameters. The properties of the solutions for the nonlinear multi-stage dynamical system are investigated and identifiability of the parameters is proved. Finally an optimal algorithm is constructed to obtain the optimal solution of the identification model and the numerical example is then discussed to illustrate the algorithm.     

Vector measure as controls for explicit nonlinear impulsive system of fed-batch culture

In this paper, we consider an optimal control problem of microbial fermentation process in which glycerol is converted to 1,3-propanediol by Klebsiella pneumoniae in fed-batch culture. During the period of reaction, the variation of pH value is monitored to determine glycerol replenishment quantity, guaranteeing that microorganism can always keep growing fast under enough nutrition. Every time pH value is lower than seven, the quantity of glycerol added is such that pH value returns seven again. Glycerol is poured into reactor at discrete time instant and the quantity is controllable. The problem is to determine for each discrete time instant the glycerol quantity to add and maximize the final concentration of 1,3-propanediol. We present a controlled explicit nonlinear impulsive dynamical system of fed-batch culture with state independent vector measures as controls and study the existence, uniqueness, boundedness, continuous dependence and Gâteaux differentiability of its solution with respect to controls. We then propose a multiple objective programming model and demonstrate the regularity of cost functionals and weak compactness of admissible control set. Finally we discuss the existence of optimal control and implement a hybrid particle swarm optimization algorithm to solve the model optimally. Computational results are presented on a numerical example.     

微生物批式流加发酵多阶段最优控制的最优性条件

在微生物批式流加发酵生产1,3一丙二醇(1,3-PD)过程中,关键是如何控制甘油和碱的流加速度.本文将流加速度看成一个随时间变化的控制函数,提出一个带控制的多阶段动力系统描述批式发酵过程,并证明了系统的一些性质.以终端时刻1,3-PD的生产强度最大为性能指标,以上述动力系统和连续状态不等式为约束条件建立了最优控制模型,最后利用不可微优化理论得到了最优控制问题的最优性条件,并证明了最优性条件和最优性函数零点的等价性.     

Global optimization of polynomial-expressed nonlinear optimal control problems with semidefinite programming relaxation

In this paper, we propose a new deterministic global optimization method for solving nonlinear optimal control problems in which the constraint conditions of differential equations and the performance index are expressed as polynomials of the state and control functions. The nonlinear optimal control problem is transformed into a relaxed optimal control problem with linear constraint conditions of differential equations, a linear performance index, and a matrix inequality condition with semidefinite programming relaxation. In the process of introducing the relaxed optimal control problem, we discuss the duality theory of optimal control problems, polynomial expression of the approximated value function, and sum-of-squares representation of a non-negative polynomial. By solving the relaxed optimal control problem, we can obtain the approximated global optimal solutions of the control and state functions based on the degree of relaxation. Finally, the proposed global optimization method is explained, and its efficacy is proved using an example of its application.     

Algorithms for Unconstrained Optimization Problems via Control Theory

Existing algorithms for solving unconstrained optimization problems are generally only optimal in the short term. It is desirable to have algorithms which are long-term optimal. To achieve this, the problem of computing the minimum point of an unconstrained function is formulated as a sequence of optimal control problems. Some qualitative results are obtained from the optimal control analysis. These qualitative results are then used to construct a theoretical iterative method and a new continuous-time method for computing the minimum point of a nonlinear unconstrained function. New iterative algorithms which approximate the theoretical iterative method and the proposed continuous-time method are then established. For convergence analysis, it is useful to note that the numerical solution of an unconstrained optimization problem is none other than an inverse Lyapunov function problem. Convergence conditions for the proposed continuous-time method and iterative algorithms are established by using the Lyapunov function theorem.