MPC for tracking zone regions

This paper deals with the problem of tracking target sets using a model predictive control (MPC) law. Some MPC applications require a control strategy in which some system outputs are controlled within specified ranges or zones (zone control), while some other variables - possibly including input variables - are steered to fixed target or set-point. In real applications, this problem is often overcome by including and excluding an appropriate penalization for the output errors in the control cost function. In this way, throughout the continuous operation of the process, the control system keeps switching from one controller to another, and even if a stabilizing control law is developed for each of the control configurations, switching among stable controllers not necessarily produces a stable closed loop system. From a theoretical point of view, the control objective of this kind of problem can be seen as a target set (in the output space) instead of a target point, since inside the zones there are no preferences between one point or another. In this work, a stable MPC formulation for constrained linear systems, with several practical properties is developed for this scenario. The concept of distance from a point to a set is exploited to propose an additional cost term, which ensures both, recursive feasibility and local optimality. The performance of the proposed strategy is illustrated by simulation of an ill-conditioned distillation column. (C) 2010 Elsevier Ltd. All rights reserved.

The Policy Iteration Algorithm for Average Continuous Control of Piecewise Deterministic Markov Processes

The main goal of this paper is to apply the so-called policy iteration algorithm (PIA) for the long run average continuous control problem of piecewise deterministic Markov processes (PDMP`s) taking values in a general Borel space and with compact action space depending on the state variable. In order to do that we first derive some important properties for a pseudo-Poisson equation associated to the problem. In the sequence it is shown that the convergence of the PIA to a solution satisfying the optimality equation holds under some classical hypotheses and that this optimal solution yields to an optimal control strategy for the average control problem for the continuous-time PDMP in a feedback form.

THE VANISHING DISCOUNT APPROACH FOR THE AVERAGE CONTINUOUS CONTROL OF PIECEWISE DETERMINISTIC MARKOV PROCESSES

This work is concerned with the existence of an optimal control strategy for the long-run average continuous control problem of piecewise-deterministic Markov processes (PDMPs). In Costa and Dufour (2008), sufficient conditions were derived to ensure the existence of an optimal control by using the vanishing discount approach. These conditions were mainly expressed in terms of the relative difference of the alpha-discount value functions. The main goal of this paper is to derive tractable conditions directly related to the primitive data of the PDMP to ensure the existence of an optimal control. The present work can be seen as a continuation of the results derived in Costa and Dufour (2008). Our main assumptions are written in terms of some integro-differential inequalities related to the so-called expected growth condition, and geometric convergence of the post-jump location kernel associated to the PDMP. An example based on the capacity expansion problem is presented, illustrating the possible applications of the results developed in the paper.

Singular Perturbation for the Discounted Continuous Control of Piecewise Deterministic Markov Processes

This paper deals with the expected discounted continuous control of piecewise deterministic Markov processes (PDMP`s) using a singular perturbation approach for dealing with rapidly oscillating parameters. The state space of the PDMP is written as the product of a finite set and a subset of the Euclidean space a""e (n) . The discrete part of the state, called the regime, characterizes the mode of operation of the physical system under consideration, and is supposed to have a fast (associated to a small parameter epsilon > 0) and a slow behavior. By using a similar approach as developed in Yin and Zhang (Continuous-Time Markov Chains and Applications: A Singular Perturbation Approach, Applications of Mathematics, vol. 37, Springer, New York, 1998, Chaps. 1 and 3) the idea in this paper is to reduce the number of regimes by considering an averaged model in which the regimes within the same class are aggregated through the quasi-stationary distribution so that the different states in this class are replaced by a single one. The main goal is to show that the value function of the control problem for the system driven by the perturbed Markov chain converges to the value function of this limit control problem as epsilon goes to zero. This convergence is obtained by, roughly speaking, showing that the infimum and supremum limits of the value functions satisfy two optimality inequalities as epsilon goes to zero. This enables us to show the result by invoking a uniqueness argument, without needing any kind of Lipschitz continuity condition.

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.

Hardware/firmware implementation of a soft sensor using an improved version of a fuzzy identification algorithm

This paper presents the design and implementation of an embedded soft sensor, i. e., a generic and autonomous hardware module, which can be applied to many complex plants, wherein a certain variable cannot be directly measured. It is implemented based on a fuzzy identification algorithm called ""Limited Rules"", employed to model continuous nonlinear processes. The fuzzy model has a Takagi-Sugeno-Kang structure and the premise parameters are defined based on the Fuzzy C-Means (FCM) clustering algorithm. The firmware contains the soft sensor and it runs online, estimating the target variable from other available variables. Tests have been performed using a simulated pH neutralization plant. The results of the embedded soft sensor have been considered satisfactory. A complete embedded inferential control system is also presented, including a soft sensor and a PID controller. (c) 2007, ISA. Published by Elsevier Ltd. All rights reserved.

Mechanical properties, drying and autogenous shrinkage of blast furnace slag activated with hydrated lime and gypsum

This article reports the characteristics of blast furnace slag (BFS) pastes activated with hydrated lime (5%) and hydrated lime (2%) plus gypsum (6%) in relation to compressive strength, shrinkage (autogenous and drying) and microstructure (porosity, hydrated products). The paste mixtures were characterized using powder X-ray diffraction (XRD), mercury intrusion porosimetry (MIP) and thermogravimetric analysis (TG/DTG). BSF activated with lime and gypsum (LG) results in larger amounts of ettringite when compared with BFS activated with lime (L). Although the porosities of the L and LG mixtures were about the same, there was a greater pore refinement for the BFS activated with lime, with an increase in mesopores volume with age. The presence of ettringite and the higher volumes of macropores cause the compressive strength of BSF activated with hydrated lime plus gypsum to be smaller than that of BFS activated with lime. For both chemical activators, compressive strength developed slowly at early ages. Autogenous and drying shrinkage were greater for the BFS activated with lime, believed to result from the more refined porous structure in comparison with the mixture activated with gypsum plus lime. (c) 2010 Elsevier Ltd. All rights reserved.

Occurrence of AFM(1) in urine samples of a Brazilian population and association with food consumption

This survey evaluated the presence of AFM(1) in human urine samples from a specific Brazilian population, as well as corn, peanut, and milk consumption measured by two types of food inquiry. Urine samples from donors who live in the city of Piracicaba, State of Sao Paulo, Brazil were analyzed to detect the presence of aflatoxin M(1) (AFM(1)). an aflatoxin B(1) metabolite, which may be used as aflatoxin B(1) exposure biomarker. The AFM(1) analysis was performed using immunoaffinity clean-up and detection by high-performance-liquid chromatography with fluorescence detector. A total of 69 samples were analyzed and 45 of them (65%) presented contaminations >= 1.8 pg ml(-1), which was the limit of quantification (LOQ). Seventy eight percent (n = 54) of the samples presented detectable concentrations of AFM(1) (>0.6 pg ml(-1)). The AFM(1) concentration among samples above LOQ ranged from 1.8 to 39.9 pg ml(-1). There were differences in food consumption profile among donors, although no association was found between food consumption and AFM(1) concentration in urine. The high frequency of positive samples suggests exposure of the populations studied to aflatoxins. (C) 2009 Elsevier Ltd. All rights reserved.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

Inorganic elements in the fat bodies of Diatraea saccharalis (Lepidoptera: Crambidae) larvae parasitized by Cotesia flavipes (Hymenoptera: Braconidae)

Koinobiont parasitoids use several strategies to regulate the host`s physiological processes during parasitism. Although many aspects of host-parasitoid interactions have been explored, studies that attempted to assess the effects of parasitism on the availability of inorganic elements in the host are virtually nonexistent. Therefore, we aimed to evaluate the effects of parasitism on the concentrations of inorganic elements in the fat bodies of larvae of Diatraea saccharalis (Lepidoptera: Crambidae) during the development of the parasitoid Cotesia flavipes (Hymenoptera: Braconidae), by using total reflection X-ray fluorescence (TXRF). TXRF analysis allowed comparisons of the changes in the availability of the elements P. S. K, Ca, Cr, Fe, Ni, Cu, and Zn in the fat body tissues of D. saccharalis larvae parasitized by C. flavipes. Overall, the concentration of inorganic elements was higher early in parasitoid development (1 and 3 days after parasitism) compared to non-parasitized larvae, but much lower towards the end of parasitoid development (7 and 9 days after parasitism). Ca, K, and S were reduced after the fifth day of parasitism, which affected the total abundance of inorganic elements observed in the fat bodies of the parasitized hosts. The regulatory mechanisms or pathological effects related to the observed variation of the host inorganic elements induced by the parasitoid remain unknown, but there might be a strategy to make these elements available to the parasitoid larvae at the end of their development, when higher metabolic activity of the host fat body is required to sustain parasitoid growth. The observed variation of the host`s inorganic elements could also be related to the known effects of parasitism on the host`s immune response. (C) 2010 Elsevier Inc. All rights reserved.

Predation potential and biology of Protogamasellopsis posnaniensis Wisniewski & Hirschmann (Acari: Rhodacaridae)

Rhodacaridae are cosmopolitan mites mentioned as predators, although nothing is known about their potential as biological control agents. One of the objectives of the work reported in this paper was to evaluate the potential of Protogamasellopsis posnaniensis (Acari: Rhodacaridae) as predator of representative species of insects of the families Sciaridae (Bradysia matogrossensis (Lane)) and Thripidae (Frankliniella occidentalis (Pergande)), of mites of the family Acaridae (Tyrophagus putrescentiae (Schrank) and Rhizoglyphus echinopus (Fumouze & Robin) and of nematodes of the family Rhabditidae (Protorhabditis sp.). Another objective was to determine the biological cycle of P. posnaniensis when fed the prey on which it performed best in the preceding predation test. The study was conducted in a laboratory where the experimental units were maintained at 25 +/- 1 degrees C, 97 +/- 3% RH and in the dark. Although the predator was able to kill all prey species considered in this study, the most favorable prey were T. putrescentiae, F. occidentalis and Protorhabditis sp. Survivorship of the predator in predation tests was always 98% or higher. Life table biological parameters when the predator was fed T. putrescentiae were: R(o) = 109.29; T = 19.06 days; lambda = 1.28 e r(m) = 0.32 female/female/day. Despite preying upon larvae of B. matogrossensis, eggs of the former can also be killed by the latter. The results indicated that A posnaniensis is a promising biological control agent, deserving additional studies on its possible use for the control of soil pests. (C) 2008 Elsevier Inc. All rights reserved.

Validation of a root water uptake model to estimate transpiration constraints

Experimental results obtained from a greenhouse trial with common bean (Phaseolus vulgaris L) plants performed to test model hypotheses regarding the onset of limiting hydraulic conditions and the shape of the transpiration reduction curve in the falling rate phase are presented. According to these hypotheses based on simulations with an upscaled single-root model, the matric flux potential at the onset of limiting hydraulic conditions is as a function of root length density and potential transpiration rate, while the relative transpiration in the falling rate phase equals the relative matric flux potential. Transpiration of bean plants in water stressed pots with four different soils was determined daily by weighing and compared to values obtained from non-stressed pots. This procedure allowed determining the onset of the falling rate phase and corresponding soil hydraulic conditions. At the onset of the falling rate phase, the value of matric flux potential M(I) showed to differ in order of magnitude from the model predicted value for three out of four soils. This difference between model and experiment can be explained by the heterogeneity of the root distribution which is not considered by the model. An empirical factor to deal with this heterogeneity should be included in the model to improve predictions. Comparing the predictions of relative transpiration in the falling rate phase using a linear shape with water content, pressure head or matric flux potential, the matric flux potential based reduction function, in agreement with the hypothesis, showed the best performance, while the pressure head based equation resulted in the highest deviations between observed and predicted values of relative transpiration rates. (C) 2010 Elsevier B.V. All rights reserved.

Evaluation of citrus somatic hybrids for tolerance to Phytophthora nicotianae and citrus tristeza virus

Somatic hybridization is a biotechnology tool that can be used in citrus breeding programs to produce somatic hybrids with the complete genetic combination of both parents. The goal of this work was to test the reaction of citrus somatic hybrids that may be useful as rootstocks to trunk and root infections caused by Phytophthora nicotianae van Breda de Haan (P parasitica Dastur) and to citrus tristeza virus (CTV). The somatic hybrids evaluated were `Caipira` sweet orange (Citrus sinensis L. Osbeck) + `Rangpur` lime (C. limonia Osbeck), `Caipira` sweet orange + `Cleopatra` mandarin (C. reshni hort. ex Tanaka), `Caipira` sweet orange + `Volkamer` lemon (C. volkameriana V Ten. & Pasq.), `Caipira` sweet orange + rough lemon (C. jambhiri Lush.), `Cleopatra` mandarin + `Volkamer` lemon, `Cleopatra` mandarin + sour orange (C. aurantium L.), `Rangpur` lime + `Sunki` mandarin (C. sunki (Hayata) hort. ex Tanaka), `Ruby Blood` sweet orange (C. sinensis L. Osbeck) + `Volkamer` lemon, `Rohde Red` sweet orange (C. sinensis L. Osbeck) + `Volkamer` lemon, and `Valencia` sweet orange + Fortunella obovata hort. ex Tanaka. For P. nicotianae trunk and root infection assays, plants of the somatic hybrids, obtained from 9-month semi-hardwood cuttings, were evaluated and compared with diploid citrus rootstock cultivars after mycelia inoculation in the trunk or spore infestation in the substrate, respectively. `Cleopatra` mandarin + sour orange, `Rangpur` lime + `Sunki` mandarin, `Cleopatra` mandarin + `Volkamer` lemon, `Ruby Blood` sweet orange + `Volkamer` lemon, `Rohde Red` sweet orange + `Volkamer` lemon, and `Caipira` sweet orange + `Volkamer` lemon had less trunk rot occurrence, whereas the somatic hybrids `Cleopatra` mandarin + `Volkamer` lemon, `Cleopatra` mandarin + sour orange, `Caipira` sweet orange + `Volkamer` lemon, and `Caipira` sweet orange + `Rangpur` lime were tolerant to root rot. For CTV assays, plants of the somatic hybrids along with tolerant and intolerant rootstocks were budded with a mild strain CTV-infected or healthy `Valencia` sweet orange budwood. Differences in average scion shoot length indicated that the hybrids `Cleopatra` mandarin + sour orange and `Valencia` sweet orange + Fortunella obovata were intolerant to CTV (c) 2007 Elsevier B.V. All rights reserved.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.