Robust integration of real time optimization with linear model predictive control

Here, we study the stable integration of real time optimization (RTO) with model predictive control (MPC) in a three layer structure. The intermediate layer is a quadratic programming whose objective is to compute reachable targets to the MPC layer that lie at the minimum distance to the optimum set points that are produced by the RTO layer. The lower layer is an infinite horizon MPC with guaranteed stability with additional constraints that force the feasibility and convergence of the target calculation layer. It is also considered the case in which there is polytopic uncertainty in the steady state model considered in the target calculation. The dynamic part of the MPC model is also considered unknown but it is assumed to be represented by one of the models of a discrete set of models. The efficiency of the methods presented here is illustrated with the simulation of a low order system. (C) 2010 Elsevier Ltd. All rights reserved.

The importance of the industrialization of Brazilian shale when faced with the world energy scenario

This article discusses the importance of the industrialization of Brazilian shale based on factors such as: security of the national energy system security, global oil geopoliticsl, resources available, production costs, oil prices, environmental impacts and the national oil reserves. The study shows that the industrialization of shale always arises when issues such as peak oil or its geopolitics appear as factors that raise the price of oil to unrealistic levels. The article concludes that in the Brazilian case, shale oil may be classified as a strategic resource, economically viable, currently in development by the success of the retorting technology for extraction of shale oil and the price of crude oil. The article presents the conclusion that shale may be the driving factor for the formation of a technology park in Sao Mateus do Sul, due to the city`s economic dependence on Petrosix.

Performance of a solar energy powered falling film evaporator with film promoter

A solar energy powered failing film evaporator with film promoter was developed for concentrating diluted solutions (industrial effluents). The procedure proposed here does not emit CO(2), making it a viable alternative to the method of concentrating solutions that uses vapor as a heat source and releases CO(2) from burning fuel oil in a furnace, in direct opposition to the carbon reduction agreement established by the Kyoto protocol. This novel device consists of the following components: a flat plate solar collector with adjustable inclination, a film promoter (adhering to the collector), a liquid distributor, a concentrate collector. and accessories. The evaporation rate of the device was found to be affected both by the inclination of the collector and by the feed flow. The meteorological variables cannot be controlled, but were monitored constantly to ascertain the behavior of the equipment in response to the variations occurring throughout the day. Higher efficiencies were attained when the inclination of the collector was adjusted monthly, showing up to 36.4% higher values than when the collector remained in a fixed position. (c) 2008 Elsevier Ltd. All rights reserved.

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.

Optimization of the Topology of Electric Energy Distribution Networks by Using Algorithm Inspired on Ant Behaviors

An algorithm inspired on ant behavior is developed in order to find out the topology of an electric energy distribution network with minimum power loss. The algorithm performance is investigated in hypothetical and actual circuits. When applied in an actual distribution system of a region of the State of Sao Paulo (Brazil), the solution found by the algorithm presents loss lower than the topology built by the concessionary company.

On the Filtering Problem for Continuous-Time Markov Jump Linear Systems with no Observation of the Markov Chain

We consider in this paper the optimal stationary dynamic linear filtering problem for continuous-time linear systems subject to Markovian jumps in the parameters (LSMJP) and additive noise (Wiener process). It is assumed that only an output of the system is available and therefore the values of the jump parameter are not accessible. It is a well known fact that in this setting the optimal nonlinear filter is infinite dimensional, which makes the linear filtering a natural numerically, treatable choice. The goal is to design a dynamic linear filter such that the closed loop system is mean square stable and minimizes the stationary expected value of the mean square estimation error. It is shown that an explicit analytical solution to this optimal filtering problem is obtained from the stationary solution associated to a certain Riccati equation. It is also shown that the problem can be formulated using a linear matrix inequalities (LMI) approach, which can be extended to consider convex polytopic uncertainties on the parameters of the possible modes of operation of the system and on the transition rate matrix of the Markov process. As far as the authors are aware of this is the first time that this stationary filtering problem (exact and robust versions) for LSMJP with no knowledge of the Markov jump parameters is considered in the literature. Finally, we illustrate the results with an example.

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.

Comments on energy confinement in non-periodic structures

In this study, we investigate the possibility of mode localization occurrence in a non-periodic Pfluger`s column model of a rocket with an intermediate concentrated mass at its middle point. We discuss the effects of varying the intermediate mass magnitude and its position and the resulting energy confinement for two cases. Free vibration analysis and the severity of mode localization are appraised, without decoupling the system, by considering as a solution basis the fundamental free response or dynamical solution. This allows for the reduction of the dimension of the algebraic modal equation that arises from satisfying the boundary and continuity conditions. By using the same methodology, we also consider the case of a cantilevered Pluger`s column with rotational stiffness at the middle support instead of an intermediate concentrated mass. (c) 2008 Elsevier Ltd. All rights reserved.

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.

Linear minimum mean square filter for discrete-time linear systems with Markov jumps and multiplicative noises

In this paper we obtain the linear minimum mean square estimator (LMMSE) for discrete-time linear systems subject to state and measurement multiplicative noises and Markov jumps on the parameters. It is assumed that the Markov chain is not available. By using geometric arguments we obtain a Kalman type filter conveniently implementable in a recurrence form. The stationary case is also studied and a proof for the convergence of the error covariance matrix of the LMMSE to a stationary value under the assumption of mean square stability of the system and ergodicity of the associated Markov chain is obtained. It is shown that there exists a unique positive semi-definite solution for the stationary Riccati-like filter equation and, moreover, this solution is the limit of the error covariance matrix of the LMMSE. The advantage of this scheme is that it is very easy to implement and all calculations can be performed offline. (c) 2011 Elsevier Ltd. All rights reserved.

Double generalized linear model for tissue culture proportion data: a Bayesian perspective

Joint generalized linear models and double generalized linear models (DGLMs) were designed to model outcomes for which the variability can be explained using factors and/or covariates. When such factors operate, the usual normal regression models, which inherently exhibit constant variance, will under-represent variation in the data and hence may lead to erroneous inferences. For count and proportion data, such noise factors can generate a so-called overdispersion effect, and the use of binomial and Poisson models underestimates the variability and, consequently, incorrectly indicate significant effects. In this manuscript, we propose a DGLM from a Bayesian perspective, focusing on the case of proportion data, where the overdispersion can be modeled using a random effect that depends on some noise factors. The posterior joint density function was sampled using Monte Carlo Markov Chain algorithms, allowing inferences over the model parameters. An application to a data set on apple tissue culture is presented, for which it is shown that the Bayesian approach is quite feasible, even when limited prior information is available, thereby generating valuable insight for the researcher about its experimental results.

Alterations in energy properties of eucalyptus wood and bark subjected to torrefaction: The potential of mass loss as a synthetic indicator

Torrefaction is a mild pyrolysis process (usually up to 300 degrees C) that changes the chemical and physical properties of biomass. This process is a possible pre-treatment prior to further processes (transport, grinding, combustion, gasification, etc) to generate energy or biofuels. In this study, three eucalyptus wood species and bark were subjected to different torrefaction conditions to determine the alterations in their structural and energy properties. The most severe treatment (280 degrees C, 5 h) causes mass losses of more than 35%, with severe damage to anatomical structure, and an increase of about 27% in the specific energy content. Bark is more sensitive to heat than wood. Energy yields are always higher than mass yields, thereby demonstrating the benefits of torrefaction in concentrating biomass energy. The overall mass loss is proposed as a relevant parameter to synthesize the effect of torrefaction conditions (temperature and duration). Accordingly, all results are summarised by analytical expressions able to predict the energy properties as a function of the overall mass loss. These expressions are intended to be used in any optimization procedure, from production in the field to the final use. (c) 2010 Elsevier Ltd. All rights reserved.

ECONOMICAL OCCUPATION OF AGRICULTURAL AREA OF PIRACICABA (SAO PAULO STATE, BRAZIL), INCLUDING RISK THROUGH LINEAR PROGRAMMING

The economic occupation of an area of 500 ha for Piracicaba was studied with the irrigated cultures of maize, tomato, sugarcane and beans, having used models of deterministic linear programming and linear programming including risk for the Target-Motad model, where two situations had been analyzed. In the deterministic model the area was the restrictive factor and the water was not restrictive for none of the tested situations. For the first situation the gotten maximum income was of R$ 1,883,372.87 and for the second situation it was of R$ 1,821,772.40. In the model including risk a producer that accepts risk can in the first situation get the maximum income of R$ 1,883,372. 87 with a minimum risk of R$ 350 year(-1), and in the second situation R$ 1,821,772.40 with a minimum risk of R$ 40 year(-1). Already a producer averse to the risk can get in the first situation a maximum income of R$ 1,775,974.81 with null risk and for the second situation R$ 1.707.706, 26 with null risk, both without water restriction. These results stand out the importance of the inclusion of the risk in supplying alternative occupations to the producer, allowing to a producer taking of decision considered the risk aversion and the pretension of income.

Growth performance and body composition of pacu Piaractus mesopotamicus (Holmberg 1887) in response to dietary protein and energy levels

Improper dietary protein and energy levels and their ratio will lead to increased fish production cost. This work evaluated effects of dietary protein : energy ratio on growth and body composition of pacu, Piaractus mesopotamicus. Fingerling pacu (15.5 +/- 0.4 g) were fed twice a day for 10 weeks until apparent satiation with diets containing 220, 260, 300, 340 or 380 g kg-1 crude protein (CP) and 10.9, 11.7, 12.6, 13.4 or 14.2 MJ kg-1 digestible energy (DE) in a totally randomized experimental design, 5 x 5 factorial scheme (n = 3). Weight gain, specific growth rate increased and feed conversion ratio (FCR) decreased significantly (P < 0.05) when CP increased from 220 to 271, 268 and 281 g kg-1 respectively. Pacu was able to adjust feed consumption in a wide range of dietary DE concentration. Fish fed 260 CP diets showed best (P < 0.05) protein efficiency ratio and FCR with 11.7-12.6 MJ kg-1; but for the 380 CP-diets group, significant differences were observed only at 14.2 MJ kg-1 dietary energy level, suggesting that pacu favours protein as energy source. DE was the chief influence on whole body chemical composition. Minimum dietary protein requirement of pacu is 270 g kg-1, with an optimum CP : DE of 22.2 g MJ-1.

Growth and haematology of pacu, Piaractus mesopotamicus, fed diets with varying protein to energy ratio

Haematopoiesis and blood cells` functions can be influenced by dietary concentration of nutrients. This paper studied the effects of dietary protein:energy ratio on the growth and haematology of pacu, Piaractus mesopotamicus. Fingerling pacu (15.5 +/- 0.4 g) were fed twice a day for 10 weeks until apparent saciety with diets containing 220, 260, 300, 340 or 380 g kg(-1) crude protein (CP) and 10.88, 11.72, 12.55, 13.39, 14.22 MJ kg(-1) digestible energy (DE) in a totally randomized experimental design, 5 x 5 factorial scheme (n=3). Weight gain and specific growth rate were affected (P < 0.05) by protein level only. Protein efficiency ratio decreased (P < 0.05) with increasing dietary protein at all levels of dietary energy. Daily feed intake decreased (P < 0.05) with increasing dietary energy. Mean corpuscular haemoglobin concentration was affected (P < 0.05) by DE and interaction between dietary CP and DE. Total plasma protein increased (P < 0.05) with dietary protein and energy levels. Plasma glucose decreased (P < 0.05) with increasing dietary protein. The CP requirement and optimum protein:energy ratio for weight gain of pacu fingerlings, determined using broken-line model, were 271 g kg(-1) and 22.18 g CP MJ(-1) DE respectively. All dietary CP and DE levels studied did not pose damages to fish health.