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.     

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.     

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.     

Optimal control for nonlinear dynamical system of 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. So a proper feeding rate is required during the process. Taking the concentration of 1, 3-PD at the terminal time as the performance index and the feeding rate of glycerol as the control function, we propose an optimal control model subject to a nonlinear dynamical system and constraints of continuous state and non-stationary control. A computational approach is constructed to seek the solution of the above model in two aspects. On the one hand we transcribe the optimal control model into an unconstrained one based on the penalty functions and an extension of the state space; on the other hand, by approximating the control function with simple functions, we transform the unconstrained optimal control problem into a sequence of nonlinear programming problems, which can be solved using gradient-based optimization techniques. The convergence analysis of this approximation is also investigated. Numerical results show that, by employing the optimal control policy, the concentration of 1, 3-PD at the terminal time can be increased considerably.     

Parameter Tuning of Multi-Proportional-Integral-Derivative Controllers Based on Optimal Switching Algorithms

Under the framework of switched systems, this paper considers a multi-proportional-integral-derivative controller parameter tuning problem with terminal equality constraints and continuous-time inequality constraints. The switching time and controller parameters are decision variables to be chosen optimally. Firstly, we transform the optimal control problem into an equivalent problem with fixed switching instants by introducing an auxiliary function and a time-scaling transformation. Because of the complexity of constraints, it is difficult to solve the problem by conventional optimization techniques. To overcome this difficulty, a novel exact penalty function is introduced for these constraints. Furthermore, the penalty function is appended to the cost functional to form an augmented cost functional, giving rise to an approximate nonlinear parameter optimization problem that can be solved using any gradient-based method. Convergence results indicate that any local optimal solution of the approximate problem is also a local optimal solution of the original problem as long as the penalty parameter is sufficiently large. Finally, an example is provided to illustrate the effectiveness of the developed algorithm.     

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.     

Optimal Control for Switched Systems with Pre-defined Order and Switch-Dependent Dynamics

The optimal control problem for switched systems, with a predefined order of switches, is considered. The differential equations may depend explicitly on previous switching instants, and the latter may be state dependent. The solution is based on the Calculus-of-Variations, which leads to a single two-point boundary-value problem. A new condition for the Hamiltonian jump at the switching instants is obtained. Simple numerical examples demonstrate the results.     

Optimal control of switched systems with time delay

In this note, we develop a computational method for solving an optimal control problem which is governed by a switched dynamical system with time delay. Our approach is to parameterize the switching instants as a new parameter vector to be optimized. Then, we derive the required gradient of the cost function which is obtained via solving a number of delay differential equations forward in time. On this basis, the optimal control problem can be solved as a mathematical programming problem.     

Stability and stabilization of impulsive switched system with inappropriate impulsive switching signals under asynchronous switching

The main objective of this paper is to study the stability and stabilization problems for a class of impulsive switched systems with inappropriate impulsive switching signals under asynchronous switching. Here, “inappropriate” means that the impulse jump moment may be inconsistent with the asynchronous switching moment or the system switching moment. And “asynchronous” implies that the switching of controller modes lags behind that of system modes. The hybrid case of stable or unstable subsystems combining with stable and unstable impulses is explored. A novel Lyapunov-like function is constructed, which is discontinuous at some special instants, including the switching instants, the instants when the system modes and filter modes are matched, and the impulse jump instants. Based on the novel multiple Lyapunov-like function, the sufficient conditions for the closed loop system to be globally uniformly exponentially stable (GUES) are obtained with admissible edge-dependent switching signals. Furthermore, by excogitating the state-feedback switching controller, the gain matrix of the controller can be obtained by solving the linear matrix inequalities. Finally, two numerical examples and simulation results are given to prove the effectiveness of our main results.     

Switching Time and Parameter Optimization in Nonlinear Switched Systems with Multiple Time-Delays

In this paper, we consider a dynamic optimization problem involving a general switched system that evolves by switching between several subsystems of nonlinear delay-differential equations. The optimization variables in this system consist of: (1) the times at which the subsystem switches occur; and (2) a set of system parameters that influence the subsystem dynamics. We first establish the existence of the partial derivatives of the system state with respect to both the switching times and the system parameters. Then, on the basis of this result, we show that the gradient of the cost function can be computed by solving the state system forward in time followed by a costate system backward in time. This gradient computation procedure can be combined with any gradient-based optimization method to determine the optimal switching times and parameters. We propose an effective optimization algorithm based on this idea. Finally, we consider three numerical examples, one involving the 1,3-propanediol fed-batch production process, to illustrate the effectiveness and applicability of the proposed algorithm.     

Nonlinear impulsive system of microbial production in fed-batch culture and its optimal control

In this study the optimal control of fed-batch glycerol fermentation is investigated based on an impulsive dynamical system. Considering the sudden increase of the glycerol and alkali in fed-batch culture of biodissimilation of glycerol to 1,3-propanediol, this paper proposes a nonlinear impulsive system of fed-batch culture. The existence, uniqueness and regularity properties of piecewise solution for the system are proved. In view of the controllability of volumes of glycerol added to the reactor instantaneously, the paper constructs an optimal control model based on the nonlinear impulsive system and the existence of the optimal control is obtained. The control variables here are the moments and the sizes of jumps in the states at the discrete instants and the objective is to maximize the productivity of 1,3-propanediol over one cycle.     

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.     

Complex metabolic network of glycerol fermentation by Klebsiella pneumoniae and its system identification via biological robustness

The bioconversion of glycerol to 1,3-propanediol (1,3-PD) by Klebsiella pneumoniae (K. pneumoniae) can be characterized by an intricate metabolic network of interactions among biochemical fluxes, metabolic compounds, key enzymes and genetic regulation. Since there are some uncertain factors in the fermentation, especially the transport mechanisms of substances across cell membrane, the metabolic network contains multiple possible metabolic systems. In this paper, we establish a complex metabolic network and the corresponding nonlinear hybrid dynamical system aiming to determine the most possible metabolic system. The existence, uniqueness and continuity of solutions are discussed. We quantitatively describe biological robustness and present a system identification model on the basis of robustness performance. The identification problem is decomposed into two subproblems and a procedure is constructed to solve them. Numerical results show that it is most possible that both glycerol and 1,3-PD pass the cell membrane by active transport coupling with passive diffusion under substrate-sufficient conditions, whereas, under substrate-limited conditions, glycerol passes cell membrane by active transport coupling with passive diffusion and 1,3-PD by active transport.