Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Toward Optimal Design of Piezoelectric Transducers Based on Multifunctional and Smoothly Graded Hybrid Material Systems

This work explores the design of piezoelectric transducers based on functional material gradation, here named functionally graded piezoelectric transducer (FGPT). Depending on the applications, FGPTs must achieve several goals, which are essentially related to the transducer resonance frequency, vibration modes, and excitation strength at specific resonance frequencies. Several approaches can be used to achieve these goals; however, this work focuses on finding the optimal material gradation of FGPTs by means of topology optimization. Three objective functions are proposed: (i) to obtain the FGPT optimal material gradation for maximizing specified resonance frequencies; (ii) to design piezoelectric resonators, thus, the optimal material gradation is found for achieving desirable eigenvalues and eigenmodes; and (iii) to find the optimal material distribution of FGPTs, which maximizes specified excitation strength. To track the desirable vibration mode, a mode-tracking method utilizing the `modal assurance criterion` is applied. The continuous change of piezoelectric, dielectric, and elastic properties is achieved by using the graded finite element concept. The optimization algorithm is constructed based on sequential linear programming, and the concept of continuum approximation of material distribution. To illustrate the method, 2D FGPTs are designed for each objective function. In addition, the FGPT performance is compared with the non-FGPT one.

Optimal design of periodic functionally graded composites with prescribed properties

The computational design of a composite where the properties of its constituents change gradually within a unit cell can be successfully achieved by means of a material design method that combines topology optimization with homogenization. This is an iterative numerical method, which leads to changes in the composite material unit cell until desired properties (or performance) are obtained. Such method has been applied to several types of materials in the last few years. In this work, the objective is to extend the material design method to obtain functionally graded material architectures, i.e. materials that are graded at the local level (e.g. microstructural level). Consistent with this goal, a continuum distribution of the design variable inside the finite element domain is considered to represent a fully continuous material variation during the design process. Thus the topology optimization naturally leads to a smoothly graded material system. To illustrate the theoretical and numerical approaches, numerical examples are provided. The homogenization method is verified by considering one-dimensional material gradation profiles for which analytical solutions for the effective elastic properties are available. The verification of the homogenization method is extended to two dimensions considering a trigonometric material gradation, and a material variation with discontinuous derivatives. These are also used as benchmark examples to verify the optimization method for functionally graded material cell design. Finally the influence of material gradation on extreme materials is investigated, which includes materials with near-zero shear modulus, and materials with negative Poisson`s ratio.

A ductile damage criterion at various stress triaxialities

The paper discusses the effect of stress triaxiality on the onset and evolution of damage in ductile metals. A series of tests including shear tests and experiments oil smooth and pre-notched tension specimens wits carried Out for it wide range of stress triaxialities. The underlying continuum damage model is based oil kinematic definition of damage tensors. The modular structure of the approach is accomplished by the decomposition of strain rates into elastic, plastic and damage parts. Free energy functions with respect to fictitious undamaged configurations as well as damaged ones are introduced separately leading to elastic material laws which are affected by increasing damage. In addition, a macroscopic yield condition and a flow rule are used to adequately describe the plastic behavior. Numerical simulations of the experiments are performed and good correlation of tests and numerical results is achieved. Based oil experimental and numerical data the damage criterion formulated in stress space is quantified. Different branches of this function are taken into account corresponding to different damage modes depending oil stress triaxiality and Lode parameter. In addition, identification of material parameters is discussed ill detail. (C) 2007 Elsevier Ltd. All rights reserved.

Application of the essential work of fracture method in ranking the performance in service of high-density polyethylene resins employed in pressure pipes

High-density polyethylene resins have increasingly been used in the production of pipes for water- and gas-pressurized distribution systems and are expected to remain in service for several years, but they eventually fail prematurely by creep fracture. Usual standard methods used to rank resins in terms of their resistance to fracture are expensive and non-practical for quality control purposes, justifying the search for alternative methods. Essential work of fracture (EWF) method provides a relatively simple procedure to characterize the fracture behavior of ductile polymers, such as polyethylene resins. In the present work, six resins were analyzed using the EWF methodology. The results show that the plastic work dissipation factor, beta w(p), is the most reliable parameter to evaluate the performance. Attention must be given to specimen preparation that might result in excessive dispersion in the results, especially for the essential work of fracture w(e).

Effect of alpha prime due to 475 degrees C aging on fracture behavior and corrosion resistance of DIN 1.4575 and MA 956 high performance ferritic stainless steels

The 475 degrees C embrittlement in stainless steels is a well-known phenomenon associated to alpha prime (alpha`) formed by precipitation or spinodal decomposition. Many doubts still remain on the mechanism of alpha` formation and its consequence on deformation and fracture mechanisms and corrosion resistance. In this investigation, the fracture behavior and corrosion resistance of two high performance ferritic stainless steels were investigated: a superferritic DIN 1.4575 and MA 956 superalloy were evaluated. Samples of both stainless steels (SS) were aged at 475 degrees C for periods varying from 1 to 1,080 h. Their fracture surfaces were observed using scanning electron microscopy (SEM) and the cleavage planes were determined by electron backscattering diffraction (EBSD). Some samples were tested for corrosion resistance using electrochemical impedance spectroscopy (EIS) and potentiodynamic polarization. Brittle and ductile fractures were observed in both ferritic stainless steels after aging at 475 degrees C. For aging periods longer than 500 h, the ductile fracture regions completely disappeared. The cleavage plane in the DIN 1.4575 samples aged at 475 degrees C for 1,080 h was mainly {110}, however the {102}, {314}, and {131} families of planes were also detected. The pitting corrosion resistance decreased with aging at 475 degrees C. The effect of alpha prime on the corrosion resistance was more significant in the DIN 1.4575 SS comparatively to the Incoloy MA 956.

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.

Stress relaxation behavior of PMMA/PS polymer blends

In this work, the stress relaxation behavior of PMMA/PS blends, with or without random copolymer addition, submitted to step shear strain experiments in the linear and nonlinear regime was studied. The effect of blend composition (ranging from 10 to 30 wt.% of dispersed phase), viscosity ratio (ranging from 0.1 to 7.5), and random copolymer addition (for concentrations up to 8 wt.% with respect to the dispersed phase) was evaluated and correlated to the evolution of the morphology of the blends. All blends presented three relaxation stages: a first fast relaxation which was attributed to the relaxation of the pure phases, a second one which was characterized by the presence of a plateau, and a third fast one. The relaxation was shown to be faster for less extended and smaller droplets and to be influenced by coalescence for blends with a dispersed phase concentration larger than 20 wt.%. The relaxation of the blend was strongly influenced by the matrix viscosity. The addition of random copolymer resulted in a slower relaxation of the droplets.

Experimental investigation of Vortex-Induced Vibration on rigid, smooth and inclined cylinders

This paper presents new experimental results of Vortex-Induced Vibration (VIV) on inclined cylinders. Models are mounted on a low damping air-bearing elastic base with one degree-of-freedom, constrained to oscillate only in the transverse direction to a free stream. The Reynolds number varied in the range 2000 less than or similar to Re less than or similar to 8000. New measurements on the dynamic response oscillations of inclined cylinders, due to VIV, are compared with previous experiments of a vertical cylinder. Models with circular and elliptical cross sections have been tested. The purpose of this work is to check the validity of the normal velocity correction of VIV studies of inclined structures. The results show that the reduced velocity range, in which the upper and lower branches of VIV occurs, is similar to the vertical cylinder case if the proper projected velocity is considered. Tests have been conducted to support this observation with inclinations up to 45 degrees. We have also observed that the amplitudes of oscillation of the inclined circular cylinder are comparable, but slightly lower than, to the amplitudes observed in the vertical cylinder experiments. Measured forces and added mass also show similar behaviour. However, for cases with an elliptical cylinder, the amplitudes of oscillation are considerably lower than those observed for a circular cylinder. This difference is explained by the higher added mass of the elliptical cylinder. (C) 2009 Elsevier Ltd. All rights reserved.

Numerical modeling of ductile crack extension in high pressure pipelines with longitudinal flaws

This study examines the applicability of a micromechanics approach based upon the computational cell methodology incorporating the Gurson-Tvergaard (GT) model and the CTOA criterion to describe ductile crack extension of longitudinal crack-like defects in high pressure pipeline steels. A central focus is to gain additional insight into the effectiveness and limitations of both approaches to describe crack growth response and to predict the burst pressure for the tested cracked pipes. A verification study conducted on burst testing of large-diameter, precracked pipe specimens with varying crack depth to thickness ratio (a/t) shows the potential predictive capability of the cell approach even though both the CT model and the CTOA criterion appear to depend on defect geometry. Overall, the results presented here lend additional support for further developments in the cell methodology as a valid engineering tool for integrity assessments of pipelines with axial defects. (C) 2011 Elsevier Ltd. All rights reserved,

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.

Galerkin Method and Weighting Functions Applied to Nonlinear H(infinity) Control with Output Feedback

This paper considers two aspects of the nonlinear H(infinity) control problem: the use of weighting functions for performance and robustness improvement, as in the linear case, and the development of a successive Galerkin approximation method for the solution of the Hamilton-Jacobi-Isaacs equation that arises in the output-feedback case. Design of nonlinear H(infinity) controllers obtained by the well-established Taylor approximation and by the proposed Galerkin approximation method applied to a magnetic levitation system are presented for comparison purposes.

Using bifurcations in the determination of lock-in ranges for third-order phase-locked loops

Transmission and switching in digital telecommunication networks require distribution of precise time signals among the nodes. Commercial systems usually adopt a master-slave (MS) clock distribution strategy building slave nodes with phase-locked loop (PLL) circuits. PLLs are responsible for synchronizing their local oscillations with signals from master nodes, providing reliable clocks in all nodes. The dynamics of a PLL is described by an ordinary nonlinear differential equation, with order one plus the order of its internal linear low-pass filter. Second-order loops are commonly used because their synchronous state is asymptotically stable and the lock-in range and design parameters are expressed by a linear equivalent system [Gardner FM. Phaselock techniques. New York: John Wiley & Sons: 1979]. In spite of being simple and robust, second-order PLLs frequently present double-frequency terms in PD output and it is very difficult to adapt a first-order filter in order to cut off these components [Piqueira JRC, Monteiro LHA. Considering second-harmonic terms in the operation of the phase detector for second order phase-locked loop. IEEE Trans Circuits Syst [2003;50(6):805-9; Piqueira JRC, Monteiro LHA. All-pole phase-locked loops: calculating lock-in range by using Evan`s root-locus. Int J Control 2006;79(7):822-9]. Consequently, higher-order filters are used, resulting in nonlinear loops with order greater than 2. Such systems, due to high order and nonlinear terms, depending on parameters combinations, can present some undesirable behaviors, resulting from bifurcations, as error oscillation and chaos, decreasing synchronization ranges. In this work, we consider a second-order Sallen-Key loop filter [van Valkenburg ME. Analog filter design. New York: Holt, Rinehart & Winston; 1982] implying a third order PLL The resulting lock-in range of the third-order PLL is determined by two bifurcation conditions: a saddle-node and a Hopf. (C) 2008 Elsevier B.V. All rights reserved.

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.