Application of a novel optimal control algorithm to a benchmark fed-batch fermentation process

A novel algorithm for solving nonlinear discrete time optimal control problems with model-reality differences is presented. The technique uses Dynamic Integrated System Optimisation and Parameter Estimation (DISOPE) which has been designed to achieve the correct optimal solution in spite of deficiencies in the mathematical model employed in the optimisation procedure. A method based on Broyden's ideas is used for approximating some derivative trajectories required. Ways for handling con straints on both manipulated and state variables are described. Further, a method for coping with batch-to- batch dynamic variations in the process, which are common in practice, is introduced. It is shown that the iterative procedure associated with the algorithm naturally suits applications to batch processes. The algorithm is success fully applied to a benchmark problem consisting of the input profile optimisation of a fed-batch fermentation process.

Novel developments in process optimisation using predictive control

In industrial practice, constrained steady state optimisation and predictive control are separate, albeit closely related functions within the control hierarchy. This paper presents a method which integrates predictive control with on-line optimisation with economic objectives. A receding horizon optimal control problem is formulated using linear state space models. This optimal control problem is very similar to the one presented in many predictive control formulations, but the main difference is that it includes in its formulation a general steady state objective depending on the magnitudes of manipulated and measured output variables. This steady state objective may include the standard quadratic regulatory objective, together with economic objectives which are often linear. Assuming that the system settles to a steady state operating point under receding horizon control, conditions are given for the satisfaction of the necessary optimality conditions of the steady-state optimisation problem. The method is based on adaptive linear state space models, which are obtained by using on-line identification techniques. The use of model adaptation is justified from a theoretical standpoint and its beneficial effects are shown in simulations. The method is tested with simulations of an industrial distillation column and a system of chemical reactors.

Dynamic integrated system optimization and parameter estimation for discrete time optimal control of nonlinear systems

An algorithm for solving nonlinear discrete time optimal control problems with model-reality differences is presented. The technique uses Dynamic Integrated System Optimization and Parameter Estimation (DISOPE), which achieves the correct optimal solution in spite of deficiencies in the mathematical model employed in the optimization procedure. A version of the algorithm with a linear-quadratic model-based problem, implemented in the C+ + programming language, is developed and applied to illustrative simulation examples. An analysis of the optimality and convergence properties of the algorithm is also presented.

Genetic algorithms for optimal control of beer fermentation

This paper uses genetic algorithms to optimise the mathematical model of a beer fermentation process that operates in batch mode. The optimisation is based in adjusting the temperature profile of the mixture during a fixed period of time in order to reach the required ethanol levels but considering certain operational and quality restrictions.

Integrating predictive control and economic optimisation

In most commercially available predictive control packages, there is a separation between economic optimisation and predictive control, although both algorithms may be part of the same software system. This method is compared in this article with two alternative approaches where the economic objectives are directly included in the predictive control algorithm. Simulations are carried out using the Tennessee Eastman process model.

Optimal control of nonlinear systems represented by differential algebraic equations

An iterative procedure is described for solving nonlinear optimal control problems subject to differential algebraic equations. The procedure iterates on an integrated modified simplified model based problem with parameter updating in such a manner that the correct solution of the original nonlinear problem is achieved.

Development and application of a novel algorithm for steady-state optimising control using dynamic information

A novel optimising controller is designed that leads a slow process from a sub-optimal operational condition to the steady-state optimum in a continuous way based on dynamic information. Using standard results from optimisation theory and discrete optimal control, the solution of a steady-state optimisation problem is achieved by solving a receding-horizon optimal control problem which uses derivative and state information from the plant via a shadow model and a state-space identifier. The paper analyzes the steady-state optimality of the procedure, develops algorithms with and without control rate constraints and applies the procedure to a high fidelity simulation study of a distillation column optimisation.

Hierarchical control and optimisation of the Tennessee Eastman process

Based on integrated system optimisation and parameter estimation a method is described for on-line steady state optimisation which compensates for model-plant mismatch and solves a non-linear optimisation problem by iterating on a linear - quadratic representation. The method requires real process derivatives which are estimated using a dynamic identification technique. The utility of the method is demonstrated using a simulation of the Tennessee Eastman benchmark chemical process.

Optimal control of nonlinear differential algebraic equation systems

A novel iterative procedure is described for solving nonlinear optimal control problems subject to differential algebraic equations. The procedure iterates on an integrated modified linear quadratic model based problem with parameter updating in such a manner that the correct solution of the original non-linear problem is achieved. The resulting algorithm has a particular advantage in that the solution is achieved without the need to solve the differential algebraic equations . Convergence aspects are discussed and a simulation example is described which illustrates the performance of the technique. 1. Introduction When modelling industrial processes often the resulting equations consist of coupled differential and algebraic equations (DAEs). In many situations these equations are nonlinear and cannot readily be directly reduced to ordinary differential equations.

A comparison of extreme European daily precipitation simulated by a global and a regional climate model for present and future climates

The intensity and distribution of daily precipitation is predicted to change under scenarios of increased greenhouse gases (GHGs). In this paper, we analyse the ability of HadCM2, a general circulation model (GCM), and a high-resolution regional climate model (RCM), both developed at the Met Office's Hadley Centre, to simulate extreme daily precipitation by reference to observations. A detailed analysis of daily precipitation is made at two UK grid boxes, where probabilities of reaching daily thresholds in the GCM and RCM are compared with observations. We find that the RCM generally overpredicts probabilities of extreme daily precipitation but that, when the GCM and RCM simulated values are scaled to have the same mean as the observations, the RCM captures the upper-tail distribution more realistically. To compare regional changes in daily precipitation in the GHG-forced period 2080-2100 in the GCM and the RCM, we develop two methods. The first considers the fractional changes in probability of local daily precipitation reaching or exceeding a fixed 15 mm threshold in the anomaly climate compared with the control. The second method uses the upper one-percentile of the control at each point as the threshold. Agreement between the models is better in both seasons with the latter method, which we suggest may be more useful when considering larger scale spatial changes. On average, the probability of precipitation exceeding the 1% threshold increases by a factor of 2.5 (GCM and RCM) in winter and by I .7 (GCM) or 1.3 (RCM) in summer.

Evaluating the potential for statistical decadal predictions of sea surface temperatures with a perfect model approach

We explore the potential for making statistical decadal predictions of sea surface temperatures (SSTs) in a perfect model analysis, with a focus on the Atlantic basin. Various statistical methods (Lagged correlations, Linear Inverse Modelling and Constructed Analogue) are found to have significant skill in predicting the internal variability of Atlantic SSTs for up to a decade ahead in control integrations of two different global climate models (GCMs), namely HadCM3 and HadGEM1. Statistical methods which consider non-local information tend to perform best, but which is the most successful statistical method depends on the region considered, GCM data used and prediction lead time. However, the Constructed Analogue method tends to have the highest skill at longer lead times. Importantly, the regions of greatest prediction skill can be very different to regions identified as potentially predictable from variance explained arguments. This finding suggests that significant local decadal variability is not necessarily a prerequisite for skillful decadal predictions, and that the statistical methods are capturing some of the dynamics of low-frequency SST evolution. In particular, using data from HadGEM1, significant skill at lead times of 6–10 years is found in the tropical North Atlantic, a region with relatively little decadal variability compared to interannual variability. This skill appears to come from reconstructing the SSTs in the far north Atlantic, suggesting that the more northern latitudes are optimal for SST observations to improve predictions. We additionally explore whether adding sub-surface temperature data improves these decadal statistical predictions, and find that, again, it depends on the region, prediction lead time and GCM data used. Overall, we argue that the estimated prediction skill motivates the further development of statistical decadal predictions of SSTs as a benchmark for current and future GCM-based decadal climate predictions.

Correlations of control variables in variational data assimilation

Variational data assimilation systems for numerical weather prediction rely on a transformation of model variables to a set of control variables that are assumed to be uncorrelated. Most implementations of this transformation are based on the assumption that the balanced part of the flow can be represented by the vorticity. However, this assumption is likely to break down in dynamical regimes characterized by low Burger number. It has recently been proposed that a variable transformation based on potential vorticity should lead to control variables that are uncorrelated over a wider range of regimes. In this paper we test the assumption that a transform based on vorticity and one based on potential vorticity produce an uncorrelated set of control variables. Using a shallow-water model we calculate the correlations between the transformed variables in the different methods. We show that the control variables resulting from a vorticity-based transformation may retain large correlations in some dynamical regimes, whereas a potential vorticity based transformation successfully produces a set of uncorrelated control variables. Calculations of spatial correlations show that the benefit of the potential vorticity transformation is linked to its ability to capture more accurately the balanced component of the flow.

Northern Hemisphere extratropical cyclones in a warming climate in the HiGEM high-resolution climate model

Changes to the Northern Hemisphere winter (December, January and February) extratropical storm tracks and cyclones in a warming climate are investigated. Two idealised climate change experiments with HiGEM1.1, a doubled CO2 and a quadrupled CO2 experiment, are compared against a present day control run. An objective feature tracking method is used and a focus given to regional changes. The climatology of extratropical storm tracks from the control run is shown to be in good agreement with ERA-40, while the frequency distribution of cyclone intensity also compares well. In both simulations the mean climate changes are generally consistent with the simulations of the IPCC AR4 models, with a strongly enhanced surface warming at the winter pole and the reduced lower tropospheric warming over the North Atlantic Ocean associated with the slowdown of the Meridional Overturning Circulation. The circulation changes in the North Atlantic are different between the two idealised simulations with different CO2 forcings. In the North Atlantic the storm tracks are influenced by the slowdown of the MOC, the enhanced surface polar warming, and the enhanced upper tropical troposphere warming, giving a north eastward shift of the storm tracks in the 2XCO2 experiment, but no shift in the 4XCO2 experiment. Over the Pacific, in the 2XCO2 experiment, changes in the mean climate are associated with local temperature changes, while in the 4XCO2 experiment the changes in the Pacific are impacted by the weakened tropical circulation. The storm track changes are consistent with the shifts in the zonal wind. Total cyclone numbers are found to decrease over the Northern Hemisphere with increasing CO2 forcing. Changes in cyclone intensity are found using 850hPa vorticity, mean sea level pressure, and 850hPa winds. The intensity of the Northern Hemisphere cyclones is found to decrease relative to the control.

Limitations of the MRL mouse as a model for cardiac regeneration

Objective Myocardial repair following injury in mammals is restricted such that damaged areas are replaced by scar tissue, impairing cardiac function. MRL mice exhibit exceptional regenerative healing in an ear punch wound model. Some myocardial repair with restoration of heart function has also been reported following cryoinjury. Increased cardiomyocyte proliferation and a foetal liver stem cell population were implicated. We investigated molecular mechanisms facilitating myocardial repair in MRL mice to identify potential therapeutic targets in non-regenerative species. Methods Expressions of specific cell-cycle regulators that might account for regeneration (CDKs 1, 2, 4 and 6; cyclins A, E, D1 and B1; p21, p27 and E2F5) were compared by immunoblotting in MRL and control C57BL/6 ventricles during development. Flow cytometry was used to investigate stem cell populations in livers from foetal mice, and infarct sizes were compared in coronary artery-ligated and sham-treated MRL and C57BL/6 adult mice. Key findings No differences in the expressions of cell cycle regulators were observed between the two strains. Expressions of CD34+Sca1+ckit-, CD34+Sca1+ckit+ and CD34+Sca1-ckit+ increased in livers from C57BL/6 vs MRL mice. No differences were observed in infarct sizes, levels of fibrosis, Ki67 staining or cardiac function between MRL and C57BL/6 mice. Conclusions No intrinsic differences were observed in cell cycle control molecules or stem cell populations between MRL and control C57BL mouse hearts. Pathophysiologically relevant ischaemic injury is not repaired more efficiently in MRL myocardium, questioning the use of the MRL mouse as a reliable model for cardiac regeneration in response to pathophysiologically relevant forms of injury.

Decentralized PID self-tuning control of industrial robots

A three degrees of freedom industrial robot is controlled by applying PID self-tuning (PID/ST) controllers. This control is considered as a corrective term to a nominal value, centrally computed from an inaccurate and/ or simplified dynamic model. An identification scheme on an assumed linear plant describing the deviation from the desired trajectory is employed in order to tune the controller coefficients and thus accomplish a behaviour prescribed through a desired pole placement. A salient feature of our approach is the decentralized nature of the controllers producing the corrective term for each joint. This opens the way to practical implementation, as recent computing requirement calculations for similar set-ups have shown in the literature. Numerical results are presented.