Stable Model Predictive Control for Integrating Systems with Optimizing Targets

This paper concern the development of a stable model predictive controller (MPC) to be integrated with real time optimization (RTO) in the control structure of a process system with stable and integrating outputs. The real time process optimizer produces Optimal targets for the system inputs and for Outputs that Should be dynamically implemented by the MPC controller. This paper is based oil a previous work (Comput. Chem. Eng. 2005, 29, 1089) where a nominally stable MPC was proposed for systems with the conventional control approach where only the outputs have set points. This work is also based oil the work of Gonzalez et at. (J. Process Control 2009, 19, 110) where the zone control of stable systems is studied. The new control for is obtained by defining ail extended control objective that includes input targets and zone controller the outputs. Additional decision variables are also defined to increase the set of feasible solutions to the control problem. The hard constraints resulting from the cancellation of the integrating modes Lit the end of the control horizon are softened,, and the resulting control problem is made feasible to a large class of unknown disturbances and changes of the optimizing targets. The methods are illustrated with the simulated application of the proposed,approaches to a distillation column of the oil refining industry.

A systemic approach for optimal cooling tower operation

The thermal performance of a cooling tower and its cooling water system is critical for industrial plants, and small deviations from the design conditions may cause severe instability in the operation and economics of the process. External disturbances such as variation in the thermal demand of the process or oscillations in atmospheric conditions may be suppressed in multiple ways. Nevertheless, such alternatives are hardly ever implemented in the industrial operation due to the poor coordination between the utility and process sectors. The complexity of the operation increases because of the strong interaction among the process variables. In the present work, an integrated model for the minimization of the operating costs of a cooling water system is developed. The system is composed of a cooling tower as well as a network of heat exchangers. After the model is verified, several cases are studied with the objective of determining the optimal operation. It is observed that the most important operational resources to mitigate disturbances in the thermal demand of the process are, in this order: the increase in recycle water flow rate, the increase in air flow rate and finally the forced removal of a portion of the water flow rate that enters the cooling tower with the corresponding make-up flow rate. (C) 2009 Elsevier Ltd. All rights reserved.

Optimal synthesis of anaerobic digester networks

The objective of this paper is to develop a mathematical model for the synthesis of anaerobic digester networks based on the optimization of a superstructure that relies on a non-linear programming formulation. The proposed model contains the kinetic and hydraulic equations developed by Pontes and Pinto [Chemical Engineering journal 122 (2006) 65-80] for two types of digesters, namely UASB (Upflow Anaerobic Sludge Blanket) and EGSB (Expanded Granular Sludge Bed) reactors. The objective function minimizes the overall sum of the reactor volumes. The optimization results show that a recycle stream is only effective in case of a reactor with short-circuit, such as the UASB reactor. Sensitivity analysis was performed in the one and two-digester network superstructures, for the following parameters: UASB reactor short-circuit fraction and the EGSB reactor maximum organic load, and the corresponding results vary considerably in terms of digester volumes. Scenarios for three and four-digester network superstructures were optimized and compared with the results from fewer digesters. (C) 2009 Elsevier B.V. All rights reserved.

Thermodynamics of aqueous solutions of polyelectrolytes: Experimental results for the activity of water in aqueous solutions of (a single synthetic polyelectrolyte and sodium chloride)

Experimental results for the activity of water in aqueous solutions of 10 single, synthetic polyelectrolytes (polysodium acrylate, polysodium methacrylate, polyammonium acrylate, polysodium ethylene sulfonate, and polysodium styrene sulfonate) and sodium chloride at 298.2 K are presented. The experimental work was performed by applying the isopiestic method with sodium chloride as a reference substance. As expected, the activity of water decreases when the concentration of a polyelectrolyte and/or sodium chloride increases. At constant concentration of a polyelectrolyte and sodium chloride, the activity of water depends on the monomer unit and the molecular mass of the polyelectrolyte. The new data are to be used in future work to develop and test models for the Gibbs excess energy of aqueous solutions of polyelectrolytes.

An extension of the Pitzer equation for the excess Gibbs energy of aqueous electrolyte systems to aqueous polyelectrolyte solutions

Pitzer`s equation for the excess Gibbs energy of aqueous solutions of low-molecular electrolytes is extended to aqueous solutions of polyelectrolytes. The model retains the original form of Pitzer`s model (combining a long-range term, based on the Debye-Huckel equation, with a short-range term similar to the virial equation where the second osmotic virial coefficient depends on the ionic strength). The extension consists of two parts: at first, it is assumed that a constant fraction of the monomer units of the polyelectrolyte is dissociated, i.e., that fraction does not depend on the concentration of the polyelectrolyte, and at second, a modified expression for the ionic strength (wherein each charged monomer group is taken into account individually) is introduced. This modification is to account for the presence of charged polyelectrolyte chains, which cannot be regarded as punctual charges. The resulting equation was used to correlate osmotic coefficient data of aqueous solutions of a single polyelectrolyte as well as of binary mixtures of a single polyelectrolyte and a salt with low-molecular weight. It was additionally applied to correlate liquid-liquid equilibrium data of some aqueous two-phase systems that might form when a polyelectrolyte and another hydrophilic but neutral polymer are simultaneously dissolved in water. A good agreement between the experimental data and the correlation result is observed for all investigated systems. (c) 2008 Elsevier B.V. All rights reserved.

Thermodynamics of aqueous solutions of polyelectrolytes: Experimental results for the activity of water in aqueous solutions of some single synthetic polyelectrolytes

Experimental results for the activity of water in aqueous solutions of 10 single polyelectrolytes (two polysodium acrylates, two polysodium methacrylates, three polyammonium acrylates, two polysodium ethylene sulfonates, and one polysodium styrene sulfonate) at (298.2 and 323.2) K are reported. The isopiestic method was employed in these experiments with aqueous solutions of sodium chloride as references. The polyelectrolytes were characterized by three averaged molecular masses determined by gel permeation chromatography. Furthermore, the density and the refractive index increments of the aqueous polyelectrolyte solutions are reported. Although a similar pattern for the activity of water was observed for all systems (i.e., the osmotic coefficient increases with rising polyelectrolyte concentration), the experimental results show that this property depends on the monomer type as well as on the size of the polymer chain. The temperature (varied from (298.2 to 323.2) K) has only a small influence on the activity of water.

Discrete-time mean variance optimal control of linear systems with Markovian jumps and multiplicative noise

In this article, we consider the stochastic optimal control problem of discrete-time linear systems subject to Markov jumps and multiplicative noise under three kinds of performance criterions related to the final value of the expectation and variance of the output. In the first problem it is desired to minimise the final variance of the output subject to a restriction on its final expectation, in the second one it is desired to maximise the final expectation of the output subject to a restriction on its final variance, and in the third one it is considered a performance criterion composed by a linear combination of the final variance and expectation of the output of the system. We present explicit sufficient conditions for the existence of an optimal control strategy for these problems, generalising previous results in the literature. We conclude this article presenting a numerical example of an asset liabilities management model for pension funds with regime switching.

Generalized Coupled Algebraic Riccati Equations for Discrete-time Markov Jump with Multiplicative Noise Systems

In this paper we consider the existence of the maximal and mean square stabilizing solutions for a set of generalized coupled algebraic Riccati equations (GCARE for short) associated to the infinite-horizon stochastic optimal control problem of discrete-time Markov jump with multiplicative noise linear systems. The weighting matrices of the state and control for the quadratic part are allowed to be indefinite. We present a sufficient condition, based only on some positive semi-definite and kernel restrictions on some matrices, under which there exists the maximal solution and a necessary and sufficient condition under which there exists the mean square stabilizing solution fir the GCARE. We also present a solution for the discounted and long run average cost problems when the performance criterion is assumed be composed by a linear combination of an indefinite quadratic part and a linear part in the state and control variables. The paper is concluded with a numerical example for pension fund with regime switching.

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes. (C) 2008 Elsevier Ltd. All rights reserved.

Extension of a Recent Method for Obtaining Exact Solutions of the Bruce and Klute Equation

A method based on a specific power-law relationship between the hydraulic head and the Boltzmann variable was recently presented. We generalized this relationship to a range of powers and extended the solution to include the saturated zone. As a result, the new solution satisfies the Bruce and Klute equation exactly.

Optimal irrigation depth for bean crops considering water constraints and farmer`s risk aversion level

Nowadays, the rising competition for the use of water and environmental resources with consequent restrictions for farmers should change the paradigm in terms of irrigation concepts, or rather, in order to attain economical efficiency other than to supply water requirement for the crop. Therefore, taking into account the social and economical role of bean activity in Brazil, as well as the risk inherent to crop due to its high sensibility to both deficit and excessive water, the optimization methods regarding to irrigation management have become more interesting and essential. This study intends to present a way to determine the optimal water supply, considering different combinations between desired bean yield and level of risk, bringing as a result a graph with the former associated with the latter, depending on different water depths.

Spatial Modeling of a Soil Fertility Index using Visible-Near-Infrared Spectra and Terrain Attributes

Our objective was to develop a methodology to predict soil fertility using visible near-infrared (vis-NIR) diffuse reflectance spectra and terrain attributes derived from a digital elevation model (DEM). Specifically, our aims were to: (i) assemble a minimum data set to develop a soil fertility index for sugarcane (Sarcharum officinarum L.) (SFI-SC) for biofuel production in tropical soils; (ii) construct a model to predict the SFI-SC using soil vis-NIR spectra and terrain attributes; and (iii) produce a soil fertility map for our study area and assess it by comparing it with a green vegetation index (GVI). The study area was 185 ha located in sao Paulo State, Brazil. In total, 184 soil samples were collected and analyzed for a range of soil chemical and physical properties. Their vis-NIR spectra were collected from 400 to 2500 nm. The Shuttle Radar Topographic Mission 3-arcsec (90-m resolution) DEM of the area was used to derive 17 terrain attributes. A minimum data set of soil properties was selected to develop the SFI-SC. The SFI-SC consisted of three classes: Class 1, the highly fertile soils; Class 2, the fertile soils; and Class 3, the least fertile soils. It was derived heuristically with conditionals and using expert knowledge. The index was modeled with the spectra and terrain data using cross-validated decision trees. The cross-validation of the model correctly predicted Class 1 in 75% of cases, Class 2 in 61%, and Class 3 in 65%. A fertility map was derived for the study area and compared with a map of the GVI. Our approach offers a methodology that incorporates expert knowledge to derive the SFI-SC and uses a versatile spectro-spatial methodology that may be implemented for rapid and accurate determination of soil fertility and better exploration of areas suitable for production.

Determination of Atropine Enantiomers in Ophthalmic Solutions by Liquid Chromatography Using a Chiral AGP (R) Column

Many therapeutic agents are commercialized under their racemic form. The enantiomers can show differences in the pharmacokinetic and pharmacodynamic profile. The use of a pure enantiomer in pharmaceutical formulations may result in a better therapeutic index and fewer adverse effects. Atropine, an alkaloid of Atropa belladonna, is a racemic mixture of l-hyoscyamine and d-hyoscyamine. It is widely used to dilate the pupil. To quantify these enantiomers in ophthalmic solutions, an HPLC method was developed and validated using a Chiral AGP (R) column at 20 degrees C. The mobile phase consisted of a buffered phosphate solution (containing 10 mM 1-octanesulfonic acid sodium salt and 7.5 mM triethylamine, adjusted to pH 7.0 with orthophosphoric acid) and acetonitrile (99 + 1, v/v). The flow rate was 0.6 mL/min, with UV detection at 205 nm. In the concentration range of 14.0-26.0 mu g/mL, the method was found to be linear (r > 0.9999), accurate (with recovery of 100.1-100.5%), and precise (RSD system: <= 0.6%; RSD intraday: <= 1.1%; RSD interday: <= 0.9%). The method was specific, and the standard and sample solutions were stable for up to 72 h. The factorial design assures robustness with a variation of +/-10% in the mobile phase components and 2 degrees C of column temperature. The complete validation, including stress testing and factorial design, was studied and is presented in this research.

Simultaneous Determination of Choline Citrate and Acetylmethionine in Injectable Solutions by High Performance Liquid Chromatography

Choline citrate (CC) and acetylmethionine (AM) are lipotropic drugs used in several pharmaceutical formulations. The objective of this research was to develop and validate a high performance liquid chromatographic (HPLC) method for simultaneous determination of CC and AM in injectable solutions, aiming its application in routine analysis for quality control of these pharmaceutical formulations. The method was validated using a Shim-Pack (R) C18 (250 x 4.6 mm, 5 mu m) column. The mobile phase was constituted of 25 mM potassium phosphate buffer solution, pH 5.7, adjusted with 10 % orthophosphoric acid, acetonitrile and methanol (88:10:2, v/v/v). The flow rate was 1.1 mL.min(-1) and the UV detection was made at 210 nm. The analyses were made at room temperature (25 +/- 1 degrees C). The method is precise, selective, accurate and robust, and was successfully applied for simultaneous quantitative determination of CC and AM in injectables.

Influence of Pluronic (R) F68 on Ceftazidime Biological Activity in Parenteral Solutions

beta-Lactam antimicrobials are known to have a low concentration/therapeutic response. However, extending the period in which beta-lactam are free in the plasma does directly influence therapeutic outcomes. The objective of this study was to evaluate the influence of Pluronic (R) F68 on the antimicrobial activity of ceftazidime when admixed with aminophylline in parenteral solutions by the evaluation of its minimal inhibitory concentration (MIC) within 24 h. Ceftazidime, aminophylline, and Pluronics (R) F68 were evaluated using the MIC method against Escherichia coli and Pseudomonas aeruginosa, with these compounds individually and associated in the same parenteral solutions. When Pluronics (R) F68 was admixtured with ceftazidime alone or with ceftazidime and aminophylline, it was possible to observe lower MIC values not only at 24 h but also at 0 h for both microorganisms. This indicates that Pluronics (R) F68 may be able to enhance ceftazidime antimicrobial activity in the presence or absence of aminophylline. This fact suggests that Pluronics (R) F68 can be applied to allow the administration of ceftazidime under continuous infusion in parenteral solutions, beneficiating hospital pharmacotherapy. It may also be possible to reduce ceftazidime doses in formulations achieving the same therapeutic results. (C) 2010 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 100:715-720, 2011