Stable MPC for tracking with maximal domain of attraction

In this work, a stable MPC that maximizes the domain of attraction of the closed-loop system is proposed. The proposed approach is suitable to real applications in the sense that it accounts for the case of output tracking, it is offset free if the output target is reachable and minimizes the offset if some of the constraints are active at steady state. The new approach is based on the definition of a Minkowski functional related to the input and terminal constraints of the stable infinite horizon MPC. It is also shown that the domain of attraction is defined by the system model and the constraints, and it does not depend on the controller tuning parameters. The proposed controller is illustrated with small order examples of the control literature. (C) 2011 Elsevier Ltd. All rights reserved.

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.

Transient and Steady-State Analysis of the Affine Combination of Two Adaptive Filters

In this paper, we propose an approach to the transient and steady-state analysis of the affine combination of one fast and one slow adaptive filters. The theoretical models are based on expressions for the excess mean-square error (EMSE) and cross-EMSE of the component filters, which allows their application to different combinations of algorithms, such as least mean-squares (LMS), normalized LMS (NLMS), and constant modulus algorithm (CMA), considering white or colored inputs and stationary or nonstationary environments. Since the desired universal behavior of the combination depends on the correct estimation of the mixing parameter at every instant, its adaptation is also taken into account in the transient analysis. Furthermore, we propose normalized algorithms for the adaptation of the mixing parameter that exhibit good performance. Good agreement between analysis and simulation results is always observed.

Improving the tracking capability of adaptive filters via convex combination

As is well known, Hessian-based adaptive filters (such as the recursive-least squares algorithm (RLS) for supervised adaptive filtering, or the Shalvi-Weinstein algorithm (SWA) for blind equalization) converge much faster than gradient-based algorithms [such as the least-mean-squares algorithm (LMS) or the constant-modulus algorithm (CMA)]. However, when the problem is tracking a time-variant filter, the issue is not so clear-cut: there are environments for which each family presents better performance. Given this, we propose the use of a convex combination of algorithms of different families to obtain an algorithm with superior tracking capability. We show the potential of this combination and provide a unified theoretical model for the steady-state excess mean-square error for convex combinations of gradient- and Hessian-based algorithms, assuming a random-walk model for the parameter variations. The proposed model is valid for algorithms of the same or different families, and for supervised (LMS and RLS) or blind (CMA and SWA) algorithms.

Jointly multi-user detection and channel estimation with genetic algorithm

This work aims at proposing the use of the evolutionary computation methodology in order to jointly solve the multiuser channel estimation (MuChE) and detection problems at its maximum-likelihood, both related to the direct sequence code division multiple access (DS/CDMA). The effectiveness of the proposed heuristic approach is proven by comparing performance and complexity merit figures with that obtained by traditional methods found in literature. Simulation results considering genetic algorithm (GA) applied to multipath, DS/CDMA and MuChE and multi-user detection (MuD) show that the proposed genetic algorithm multi-user channel estimation (GAMuChE) yields a normalized mean square error estimation (nMSE) inferior to 11%, under slowly varying multipath fading channels, large range of Doppler frequencies and medium system load, it exhibits lower complexity when compared to both maximum likelihood multi-user channel estimation (MLMuChE) and gradient descent method (GrdDsc). A near-optimum multi-user detector (MuD) based on the genetic algorithm (GAMuD), also proposed in this work, provides a significant reduction in the computational complexity when compared to the optimum multi-user detector (OMuD). In addition, the complexity of the GAMuChE and GAMuD algorithms were (jointly) analyzed in terms of number of operations necessary to reach the convergence, and compared to other jointly MuChE and MuD strategies. The joint GAMuChE-GAMuD scheme can be regarded as a promising alternative for implementing third-generation (3G) and fourth-generation (4G) wireless systems in the near future. Copyright (C) 2010 John Wiley & Sons, Ltd.

Fear of heights: cognitive performance and postural control

Fear of heights, or acrophobia, is one of the most frequent subtypes of specific phobia frequently associated to depression and other anxiety disorders. Previous evidence suggests a correlation between acrophobia and abnormalities in balance control, particularly involving the use of visual information to keep postural stability. This study investigates the hypotheses that (1) abnormalities in balance control are more frequent in individuals with acrophobia even when not exposed to heights, that (2) acrophobic symptoms are associated to abnormalities in visual perception of movement; and that (3) individuals with acrophobia are more sensitive to balance-cognition interactions. Thirty-one individuals with specific phobia of heights and thirty one non-phobic controls were compared using dynamic posturography and a manual tracking task. Acrophobics had poorer performance in both tasks, especially when carried out simultaneously. Previously described interference between posture control and cognitive activity seems to play a major role in these individuals. The presence of physiologic abnormalities is compatible with the hypothesis of a non-associative acquisition of fear of heights, i.e., not associated to previous traumatic events or other learning experiences. Clinically, this preliminary study corroborates the hypothesis that vestibular physical therapy can be particularly useful in treating individuals with fear of heights.

Additive Effects of Obstructive Sleep Apnea and Hypertension on Early Markers of Carotid Atherosclerosis

Obstructive sleep apnea (OSA) has emerged as an independent risk factor for atherosclerosis. However, OSA is frequently associated with several risk factors for atherosclerosis, including hypertension (HTN). The impact of OSA and HTN alone compared with the association of both conditions on carotid atherosclerosis is not understood. We studied 94 middle-aged participants free of smoking and diabetes mellitus who were divided into 4 groups: controls (n = 22), OSA (n = 25), HTN (n = 20), and OSA + HTN (n = 27). All of the participants underwent polysomnography and carotid measurements of intima-media thickness, diameter, and distensibility with an echo-tracking device. Compared with controls, intima-media thickness and carotid diameter were similarly higher in OSA (713 +/- 117 and 7117 +/- 805 mu m), and HTN groups (713 +/- 182 and 7191 +/- 818 mu m), with a further significant increase in OSA + HTN patients (837 +/- 181 and 7927 +/- 821 mu m, respectively; P < 0.01). Carotid distensibility was significantly lower in HTN (P < 0.05) and OSA + HTN subjects (P < 0.001) compared with controls. In the OSA + HTN group, carotid distensibility was significantly lower than in the OSA group and controls (P < 0.05 for each comparison). Multivariate analysis showed that intima-media thickness was positively related to systolic blood pressure and apnea-hypopnea index. Apnea-hypopnea index was the only factor related to carotid diameter. Age and systolic blood pressure were independently related to carotid distensibility. In conclusion, the association of OSA and HTN has additive effects on markers of carotid atherosclerosis. Because early markers of carotid atherosclerosis predict future cardiovascular events, including not only stroke but also myocardial infarction, these findings may help to explain the increased risk of cardiovascular disease in patients with OSA. (Hypertension. 2009; 53: 64-69.)

In vivo electrochemical characterization and inflammatory response of multiwalled carbon nanotube-based electrodes in rat hippocampus

Stimulating neural electrodes are required to deliver charge to an environment that presents itself as hostile. The electrodes need to maintain their electrical characteristics (charge and impedance) in vivo for a proper functioning of neural prostheses. Here we design implantable multi-walled carbon nanotubes coating for stainless steel substrate electrodes, targeted at wide frequency stimulation of deep brain structures. In well-controlled, low-frequency stimulation acute experiments, we show that multi-walled carbon nanotube electrodes maintain their charge storage capacity (CSC) and impedance in vivo. The difference in average CSCs (n = 4) between the in vivo (1.111 mC cm(-2)) and in vitro (1.008 mC cm(-2)) model was statistically insignificant (p > 0.05 or P-value = 0.715, two tailed). We also report on the transcription levels of the pro-inflammatory cytokine IL-1 beta and TLR2 receptor as an immediate response to low-frequency stimulation using RT-PCR. We show here that the IL-1 beta is part of the inflammatory response to low-frequency stimulation, but TLR2 is not significantly increased in stimulated tissue when compared to controls. The early stages of neuroinflammation due to mechanical and electrical trauma induced by implants can be better understood by detection of pro-inflammatory molecules rather than by histological studies. Tracking of such quantitative response profits from better analysis methods over several temporal and spatial scales. Our results concerning the evaluation of such inflammatory molecules revealed that transcripts for the cytokine IL-1 beta are upregulated in response to low-frequency stimulation, whereas no modulation was observed for TLR2. This result indicates that the early response of the brain to mechanical trauma and low-frequency stimulation activates the IL-1 beta signaling cascade but not that of TLR2.

Summertime moisture transport over Southeastern South America and extratropical cyclones behavior during inter-El Nio events

The impact of the inter-El Nio (EN) variability on the moisture availability over Southeastern South America (SESA) is investigated. Also, an automatic tracking scheme was used to analyze the extratropical cyclones properties (system density - SD and central pressure - CP) in this region. During the austral summer period from 1977-2000, the differences for the upper-level wave train anomaly composites seem to determine the rainfall composite differences. In fact, the positive rainfall anomalies over most of the SESA domain during the strong EN events are explained by an upper-level cyclonic center over the tropics and an anticyclonic center over the eastern subtropical area. This pattern seems to contribute to upward vertical motion at 500 hPa and reinforcement of the meridional moisture transport from the equatorial Atlantic Ocean and western Amazon basin to the SESA region. These features may contribute to the positive SD and negative CP anomalies explaining part of the positive rainfall anomalies found there. On the other hand, negative rainfall anomalies are located in the northern part of SESA for the weak EN years when compared to those for the strong events. Also, positive anomalies are found in the southern part, albeit less intense. It was associated with the weakening of the meridional moisture transport from the tropics to the SESA that seems have to contributed with smaller SD and CP anomalies over the most part of subtropics, when compared to the strong EN years.

Front tracking with moving-least-squares surfaces

The representation of interfaces by means of the algebraic moving-least-squares (AMLS) technique is addressed. This technique, in which the interface is represented by an unconnected set of points, is interesting for evolving fluid interfaces since there is]to surface connectivity. The position of the surface points can thus be updated without concerns about the quality of any surface triangulation. We introduce a novel AMLS technique especially designed for evolving-interfaces applications that we denote RAMLS (for Robust AMLS). The main advantages with respect to previous AMLS techniques are: increased robustness, computational efficiency, and being free of user-tuned parameters. Further, we propose a new front-tracking method based on the Lagrangian advection of the unconnected point set that defines the RAMLS surface. We assume that a background Eulerian grid is defined with some grid spacing h. The advection of the point set makes the surface evolve in time. The point cloud can be regenerated at any time (in particular, we regenerate it each time step) by intersecting the gridlines with the evolved surface, which guarantees that the density of points on the surface is always well balanced. The intersection algorithm is essentially a ray-tracing algorithm, well-studied in computer graphics, in which a line (ray) is traced so as to detect all intersections with a surface. Also, the tracing of each gridline is independent and can thus be performed in parallel. Several tests are reported assessing first the accuracy of the proposed RAMLS technique, and then of the front-tracking method based on it. Comparison with previous Eulerian, Lagrangian and hybrid techniques encourage further development of the proposed method for fluid mechanics applications. (C) 2008 Elsevier Inc. All rights reserved.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

The magnetized steel and scintillator calorimeters of the MINOS experiment

The Main Injector Neutrino Oscillation Search (MINOS) experiment uses an accelerator-produced neutrino beam to perform precision measurements of the neutrino oscillation parameters in the ""atmospheric neutrino"" sector associated with muon neutrino disappearance. This long-baseline experiment measures neutrino interactions in Fermilab`s NuMI neutrino beam with a near detector at Fermilab and again 735 km downstream with a far detector in the Soudan Underground Laboratory in northern Minnesota. The two detectors are magnetized steel-scintillator tracking calorimeters. They are designed to be as similar as possible in order to ensure that differences in detector response have minimal impact on the comparisons of event rates, energy spectra and topologies that are essential to MINOS measurements of oscillation parameters. The design, construction, calibration and performance of the far and near detectors are described in this paper. (C) 2008 Elsevier B.V. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.