Predictive control of transport-reaction processes

This work focuses on the development of computationally efficient predictive control algorithms for nonlinear parabolic and hyperbolic PDEs with state and control constraints arising in the context of transport-reaction processes. We first consider a diffusion-reaction process described by a nonlinear parabolic PDE and address the problem of stabilization of an unstable steady-state subject to input and state constraints. Galerkin’s method is used to derive finite-dimensional systems that capture the dominant dynamics of the parabolic PDE, which are subsequently used for controller design. Various model predictive control (MPC) formulations are constructed on the basis of the finite dimensional approximations and are demonstrated, through simulation, to achieve the control objectives. We then consider a convection-reaction process example described by a set of hyperbolic PDEs and address the problem of stabilization of the desired steady-state subject to input and state constraints, in the presence of disturbances. An easily implementable predictive controller based on a finite dimensional approximation of the PDE obtained by the finite difference method is derived and demonstrated, via simulation, to achieve the control objective.     

Predictive control of transport-reaction processes

This work focuses on the development of computationally efficient predictive control algorithms for nonlinear parabolic and hyperbolic PDEs with state and control constraints arising in the context of transport-reaction processes. We first consider a diffusion-reaction process described by a nonlinear parabolic PDE and address the problem of stabilization of an unstable steady-state subject to input and state constraints. Galerkin’s method is used to derive finite-dimensional systems that capture the dominant dynamics of the parabolic PDE, which are subsequently used for controller design. Various model predictive control (MPC) formulations are constructed on the basis of the finite dimensional approximations and are demonstrated, through simulation, to achieve the control objectives. We then consider a convection-reaction process example described by a set of hyperbolic PDEs and address the problem of stabilization of the desired steady-state subject to input and state constraints, in the presence of disturbances. An easily implementable predictive controller based on a finite dimensional approximation of the PDE obtained by the finite difference method is derived and demonstrated, via simulation, to achieve the control objective.     

Optimal boundary control of a diffusion–convection-reaction PDE model with time-dependent spatial domain: Czochralski crystal growth process

In this paper the optimal boundary control problem for diffusion–convection-reaction processes modeled by partial differential equations (PDEs) defined on time-dependent spatial domains is considered. The model of the transport system with time-varying domain arises in the context of high energy consuming Czochralski crystal growth process in which the crystal temperature regulation must successfully account for the change in the crystal spatial domain due to the crystal growth process realized by the pulling crystal out of melt. Starting from the first principles of continuum mechanics and transport theorem the time-varying parabolic PDE describing temperature evolution is derived and represented as a nonautonomous parabolic evolution system on an appropriately defined function space which is exactly transformed in the infinite-dimensional boundary control problem for which a boundary linear quadratic regulator is proposed. Properties of the solution of the time-varying parabolic PDEs given by the two-parameter evolutionary system are utilized in the synthesis of the optimal boundary regulator, and the control law is applied to the model given by a two-dimensional partial differential equation in the cylindrical coordinates representing the Czochralski crystal growth process with one-dimensional growth direction. Finally, numerical results demonstrate optimal stabilization of the two-dimensional temperature distribution in the crystal.     

Fault-tolerant control of fluid dynamic systems via coordinated feedback and switching

This work addresses the problem of designing a fault-tolerant control system for fluid dynamic systems modeled by highly-dissipative partial differential equations (PDEs) with constrained control actuators. The proposed approach is predicated upon the idea of coordinating feedback controller synthesis and switching between multiple, spatially-distributed control actuator configurations. Using appropriate finite-dimensional approximations of the PDE system, a stabilizing feedback controller is designed for a given actuator configuration, and an explicit characterization of the constrained stability region is obtained. Switching laws are then derived, on the basis of these stability regions, to orchestrate the switching between the control actuator configurations, in a way that guarantees constraint satisfaction and preserves closed-loop stability of the infinite-dimensional system in the event of actuator failures. The results are demonstrated through an application of the proposed methodology to the suppression of wave formation in falling liquid films via the stabilization of the zero solution of the one-dimensional Kuramoto–Sivashinsky equation (KSE), with periodic boundary conditions, subject to actuator constraints and failures.     

Model predictive control of Kuramoto–Sivashinsky equation with state and input constraints

This work focuses on the model predictive control design methodology that successfully accounts for the state and input constraints applied in the context of highly dissipative Kuramoto–Sivashinsky (KS) partial differential equation (PDE) describing stability of a thin film thickness in the two-phase annular flow in vertical pipes. The evolution of a linear dissipative KSE PDE state is given by an abstract evolution equation in an appropriate functional space. The proposed constrained predictive control law utilizes a low order modal representation in the optimization functional, while higher modes are included only in the PDE state constraints. Simulation results demonstrate a successful application of the proposed predictive control technique that achieves optimal stabilization of a spatially-uniform unstable steady state of Kuramoto–Sivashinsky equation in the presence of input and state constraints.     

基于PDE降阶模型的最优控制

针对带有约束条件的偏微分方程（PDE）模型最优控制的实时性要求和巨大的内存开销问题,提出了基于降阶模型的输入/状态约束的最优实时控制方法。采用特征正交分解和奇异值分解以及Galerkin投影方法导出了动态PDE具有高精度离散形式的状态空间低阶模型,提出了一定输入/状态约束条件下的基于二次规划单步滚动最优控制并与基于线性二次调节器的极值验证最优控制策略相互验证。通过对流-扩散-反应过程的控制仿真结果证明了所提方法的高效性和准确性。     

Nonlinear predictive control using multi-rate sampling

A nonlinear predictive control (NLPC) strategy based on a nonlinear, lumped parameter model of the process is developed in this paper. A constrained optimization approach is used to estimate unmeasured state variables and load disturbances. Additional model/process mismatch is handled by using an additive output term which is equivalent to the Internal Model Control approach. Similar to linear predictive control methods, an optimal sequence of future control moves is determined in order to minimize an objective function based on a desired output trajectory, subject to manipulated variable constraints (absolute and velocity). Deadtime is explicitly included in the model formulation, giving NLPC the same deadtime compensation feature of linear model-predictive techniques. The multi-rate sampling nature of most chemical processes is also used to improve estimates of process disturbances. Infrequent composition measurements in conjunction with frequent temperature measurements are used to improve the “inferential” control of the composition in a continuous flow stirred tank reactor (CSTR).     

Dynamic output feedback covariance control of stochastic dissipative partial differential equations

In this work, we develop a method for dynamic output feedback covariance control of the state covariance of linear dissipative stochastic partial differential equations (PDEs) using spatially distributed control actuation and sensing with noise. Such stochastic PDEs arise naturally in the modeling of surface height profile evolution in thin film growth and sputtering processes. We begin with the formulation of the stochastic PDE into a system of infinite stochastic ordinary differential equations (ODEs) by using modal decomposition. A finite-dimensional approximation is then obtained to capture the dominant mode contribution to the surface roughness profile (i.e., the covariance of the surface height profile). Subsequently, a state feedback controller and a Kalman-Bucy filter are designed on the basis of the finite-dimensional approximation. The dynamic output feedback covariance controller is subsequently obtained by combining the state feedback controller and the state estimator. The steady-state expected surface covariance under the dynamic output feedback controller is then estimated on the basis of the closed-loop finite-dimensional system. An analysis is performed to obtain a theoretical estimate of the expected surface covariance of the closed-loop infinite-dimensional system. Applications of the linear dynamic output feedback controller to both the linearized and the nonlinear stochastic Kuramoto-Sivashinsky equations (KSEs) are presented. Finally, nonlinear state feedback controller and nonlinear output feedback controller designs are also presented and applied to the nonlinear stochastic KSE.     

Model predictive control of a second-order hyperbolic transport-reaction process

The model predictive controller (MPC) design is developed for a tubular chemical reactor, considering a second-order hyperbolic partial differential equation as the model of the transport-reaction process with boundary actuation. Without loss of generality, closed–closed boundary conditions and relaxed total flux are assumed. At the same time, the model is discretized in time by the Cayley–Tustin method, and, under the assumption that only the reactor's output is measurable, the observer design for the state reconstruction is addressed and integrated with the MPC design. The Luenberger observer gain is obtained by solving the operator Ricatti equation in the discrete-time setting, while the MPC accounts for constrained and optimal control. The simulations show that the output-based MPC design stabilizes the system under the input and output constraints satisfaction. In addition, to address the models' disparities, the results for both parabolic and hyperbolic equations are presented and discussed.     

Optimal mechano-electric stabilization of cardiac alternans

Alternation of normal action-potential morphology in the myocardium is a condition with a beat-to-beat oscillation in the length of the electric wave which is linked through electromechanical coupling to the cardiac muscle contraction, and is believed to be the first manifestation of the onset of life threatening ventricular arrhythmias and sudden cardiac death. In this work, the effects of electrical and mechanical stimuli are utilized in alternans annihilation problem. Electrical stimuli that alter the action-potential morphology are represented by a pacer located at the domain's boundary, while mechanical stimuli are distributed within the spatial domain and affect the action potential by altering intracellular calcium kinetics. Alternation of action potential is described by the small amplitude of alternans parabolic partial differential equation (PDE). Spatially uniform unstable steady state of the alternans amplitude PDE is stabilized by optimal control methods through boundary and spatially distributed actuation. Mixed boundary and spatially distributed actuation is manipulated by a linear quadratic regulator (LQR) in the full-state-feedback control structure and in a compensator design with a finite-dimensional Luenberger-type observer, and it achieves exponential stabilization in a finite size tissue cable length. The proposed control problem formulation and the performance and robustness of the closed-loop system under the proposed linear controller are evaluated through simulations.     

Boundary control synthesis for a lithium‐ion battery thermal regulation problem

The thermal regulation problem for a lithium ion (Li‐ion) battery with boundary control actuation is considered. The model of the transient temperature dynamics of the battery is given by a nonhomogeneous parabolic partial differential equation (PDE) on a two‐dimensional spatial domain which accounts for the time‐varying heat generation during the battery discharge cycle. The spatial domain is given as a disk with radial and angular coordinates which captures the nonradially symmetric heat‐transfer phenomena due to the application of the control input along a portion of the spatial domain boundary. The Li‐ion battery model is formulated within an appropriately defined infinite‐dimensional function space setting which is suitable for spectral controller synthesis. The key challenges in the output feedback model‐based controller design addressed in this work are: the dependence of the state on time‐varying system parameters, the restriction of the input along a portion of the battery domain boundary, the observer‐based optimal boundary control design where the separation principle is utilized to demonstrate the stability of the closed loop system, and the realization of the outback feedback control problem based on state measurement and interpolation of the temperature field. Numerical results for simulation case studies are presented. © 2013 American Institute of Chemical Engineers AIChE J, 59: 3782–3796, 2013     

Predictive control of particle size distribution in particulate processes

In this work, we focus on the development and application of predictive-based strategies for control of particle size distribution (PSD) in continuous and batch particulate processes described by population balance models (PBMs). The control algorithms are designed on the basis of reduced-order models, utilize measurements of principle moments of the PSD, and are tailored to address different control objectives for the continuous and batch processes. For continuous particulate processes, we develop a hybrid predictive control strategy to stabilize a continuous crystallizer at an open-loop unstable steady-state. The hybrid predictive control strategy employs logic-based switching between model predictive control (MPC) and a fall-back bounded controller with a well-defined stability region. The strategy is shown to provide a safety net for the implementation of MPC algorithms with guaranteed stability closed-loop region. For batch particulate processes, the control objective is to achieve a final PSD with desired characteristics subject to both manipulated input and product quality constraints. An optimization-based predictive control strategy that incorporates these constraints explicitly in the controller design is formulated and applied to a seeded batch crystallizer. The strategy is shown to be able to reduce the total volume of the fines by 13.4% compared to a linear cooling strategy, and is shown to be robust with respect to modeling errors.     

基于LPV模型的燃料电池空气进气系统控制

质子交换膜燃料电池是一种通过氢气和氧气的电化学反应将化学能直接转化为电能的装置。提出一种改进的四阶燃料电池进气系统模型,分析了系统的约束性。针对系统模型所具有的非线性特性,提出建立线性变参数（LPV）模型用于对系统的控制。针对状态变量不可测的问题引入卡尔曼滤波器,同时通过可观性分析得出系统所需测量的最佳变量。在符合约束条件下设计基于线性变参数模型的状态空间模型预测控制器,控制空压机的工作电压保证氢气燃料的充分反应。仿真结果表明,基于LPV模型的模型预测控制器能够对空气进气系统进行有效的控制,且满足空压机喘振和阻塞边界等约束条件,与单模型预测控制相比具有更好的控制效果。     

Linear model predictive control for transport‐reaction processes

The article deals with systematic development of linear model predictive control algorithms for linear transport‐reaction models emerging from chemical engineering practice. The finite‐horizon constrained optimal control problems are addressed for the systems varying from the convection dominated models described by hyperbolic partial differential equations (PDEs) to the diffusion models described by parabolic PDEs. The novelty of the design procedure lies in the fact that spatial discretization and/or any other type of spatial approximation of the process model plant is not considered and the system is completely captured with the proposed Cayley‐Tustin transformation, which maps a plant model from a continuous to a discrete state space setting. The issues of optimality and constrained stabilization are addressed within the controller design setting leading to the finite constrained quadratic regulator problem, which is easily realized and is no more computationally intensive than the existing algorithms. The methodology is demonstrated for examples of hyperbolic/parabolic PDEs. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2644–2659, 2017     

Robust model predictive control of nonlinear process systems: Handling rate constraints

This work considers the problem of stabilization of control affine nonlinear process systems subject to constraints on the rate of change and magnitude of control inputs in the presence of uncertainty. We first handle rate constraints within a soft constraints framework. A new robust predictive controller formulation that minimizes rate constraint violation while guaranteeing stabilization and input constraint satisfaction from an explicitly characterized stability region is designed. We then derive conditions that allow for guaranteed satisfaction of hard rate constraints. Subsequently, a predictive controller is designed that ensures rate constraints satisfaction when the required conditions are satisfied, relaxing them otherwise to preserve feasibility and robust stability. The implementation of the proposed predictive controllers is illustrated via a chemical reactor example.     

LMI-based robust model predictive control for a continuous MMA polymerization reactor

A linear matrix inequality (LMI)-based robust model predictive control (MPC) is applied to a continuous stirred-tank reactor for the polymerization of methyl methacrylate (MMA). The polytopic model is constructed to predict the responses to various control input sequences by using Jacobians of uncertain nonlinear model at several operating points and the controller design is characterized as the problem of minimizing an upper bound on the ‘worst-case’ infinite horizon objective function subject to constraints on the control input and plant output. Simulation studies under different conditions are conducted to validate the feasibility of the optimization problem and evaluate the applicability of such a control scheme. Simulation results show that, despite the model uncertainty, the LMI-based robust model predictive controller performs quite satisfactorily for the property control of the continuous polymerization reactor and guarantees the robust stability.     

改进的模糊预测控制算法及其在CSTR过程中的应用

A finite horizon predictive control algorithm,which applies a saturated feedback control law as its local control law,is presented for nonlinear systems with time-delay subject to input constraints.In the algorithm,N free control moves,a saturated local control law and the terminal weighting matrices are solved by a minimization problem based on linear matrix inequality(LMI) constraints online.Compared with the algorithm with a nonsaturated local law,the presented algorithm improves the performances of the closed-loop systems such as feasibility and optimality.This model predictive control(MPC) algorithm is applied to an industrial continuous stirred tank reactor(CSTR) with explicit input constraint.The simulation results demonstrate that the presented algorithm is effective.     

Backstepping control of chemical tubular reactors

In this paper, a globally stabilizing boundary feedback control law for an arbitrarily fine discretization of a nonlinear PDE model of a chemical tubular reactor is presented. A model that assumes no radial velocity and concentration gradients in the reactor, the temperature gradient described by use of a proper value of the effective radial conductivity, a homogeneous reaction, the properties of the reaction mixture characterized by average values, the mechanism of axial mixing described by a single parameter model, and the kinetics of the first order is considered. Depending on the values of the nondimensional Peclet numbers, Damköhler number, the dimensionless adiabatic temperature rise, and the dimensionless activation energy, the coupled PDE equations for the temperature and concentration can have multiple equilibria that can be either stable or unstable. The objective is to stabilize an unstable steady state of the system using boundary control of temperature and concentration on the inlet side of the reactor. We discretize the original nonlinear PDE model in space using finite difference approximation and get a high order system of coupled nonlinear ODEs. Then, using backstepping design for parabolic PDEs we transform the original coupled system into two uncoupled target systems that are asymptotically stable in l²-norm with appropriate homogeneous boundary conditions. In the real system, the designed control laws would be implemented through small variations of the prescribed inlet temperature and prescribed inlet concentration. The control design is accompanied by a simulation study that shows the feedback control law designed with sensing only on a very coarse grid (using just a few measurements of the temperature and concentration fields) can successfully stabilize the actual system for a variety of different simulation settings (on a fine grid).     

Soft-constrained model predictive control based on data-driven distributionally robust optimization

This article proposes a novel distributionally robust optimization (DRO)-based soft-constrained model predictive control (MPC) framework to explicitly hedge against unknown external input terms in a linear state-space system. Without a priori knowledge of the exact uncertainty distribution, this framework works with a lifted ambiguity set constructed using machine learning to incorporate the first-order moment information. By adopting a linear performance measure and considering input and state constraints robustly with respect to a lifted support set, the DRO-based MPC is reformulated as a robust optimization problem. The constraints are softened to ensure recursive feasibility. Theoretical results on optimality, feasibility, and stability are further discussed. Performance and computational efficiency of the proposed method are illustrated through motion control and building energy control systems, showing 18.3% less cost and 78.8% less constraint violations, respectively, while requiring one third of the CPU time compared to multi-stage scenario based stochastic MPC.     

Transition from Batch to Continuous Operation in Bio‐Reactors: A Model Predictive Control Approach and Application

This work considers the problem of determining the transition of ethanol‐producing bio‐reactors from batch to continuous operation and subsequent control subject to constraints and performance considerations. To this end, a Lyapunov‐based non‐linear model predictive controller is utilized that stabilizes the bio‐reactor under continuous mode of operation. The key idea in the predictive controller is the formulation of appropriate stability constraints that allow an explicit characterization of the set of initial conditions from where feasibility of the optimization problem and hence closed‐loop stability is guaranteed. Additional constraints are incorporated in the predictive control design to expand on the set of initial conditions that can be stabilized by control designs that only require the value of the Lyapunov function to decay. Then, the explicit characterization of the set of stabilizable initial conditions is used in determining the appropriate time for which the reactor must be run in batch mode. Specifically, the predictive control approach is utilized in determining the appropriate batch length that achieves stabilizable values of the state variables at the end of the batch. Application of the proposed method to the ethanol production process using Zymomonas mobilis as the ethanol producing micro‐organism demonstrates the effectiveness of the proposed model predictive control strategy in stabilizing the bio‐reactor.