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 numerical investigation on the visco-plastic response of structures to different pulse loading shapes

The time varying intensity character of a load applied to a structure poses many difficulties in analysis. A remedy to this situation is to substitute a complex pulse shape by a rectangular equivalent one. It has been shown by others that this procedure works well for perfectly plastic elementary structures. This paper applies the concept of equivalent pulse to more complex structures. Special attention is given to the material behavior, which is allowed to be strain rate and strain hardening sensitive. Thanks to the explicit finite element solution, it is shown in this article that blast loads applied to complex structures made of real materials can be substituted by equivalent rectangular loads with both responses being practically the same. (c) 2007 Elsevier Ltd. All rights reserved.

Synthesis, characterization and evaluation of phosphorylated resins in the removal of Pb(2+) from aqueous solutions

Phosphinic-derivative poly(styrene-co-divinylbenzene)-based on PS-DVB copolymers with different porosity degrees have been prepared by aromatic electrophilic substitution reaction using PCl(3)/AlCl(3) followed by base-promoted hydrolysis. The phosphorylation reaction was analyzed by infra-red spectroscopy (FTIR), scanning electron microscopy (SEM), and thermogravimetry (TG/DTG). In addition, the phosphorous content of the phosphorylated copolymers was determined by spectrophotometry using the method based on sodium molybdate reactant so that the extension of that modification could be assessed. The performance of the phosphorylated resins in the extraction of Pb(2+) from aqueous solutions in a batch system was also evaluated. The Pb(2+) content was determined by atomic absorption spectrometry (AAS). These materials presented excellent extraction capacity under the contact time of 30 min and pH 6.

The electrochemical behaviour of Ti-13Nb-13Zr alloy in various solutions

The electrochemical behaviour of a near-beta Ti-13Nb-13Zr alloy for the application as implants was investigated in various solutions. The electrolytes used were 0.9 wt% NaCl solution, Hanks` solution and a culture medium known as minimum essential medium (MEM) composed of salts, vitamins and amino acids, all at 37 degrees C. The electrochemical behaviour was investigated by the following electrochemical techniques: open circuit potential measurements as a function of time, electrochemical impedance spectroscopy (EIS) and determination of polarisation curves. The obtained results showed that the Ti alloy was passive in all electrolytes. The EIS results were analysed using an equivalent electrical circuit representing a duplex structure oxide layer, composed of an inner barrier layer, mainly responsible for the alloy corrosion resistance, and an outer and porous layer that has been associated to osteointegration ability. The properties of both layers were dependent on the electrolyte used. The results suggested that the thickest porous layer is formed in the MEM solution whereas the impedance of the barrier layer formed in this solution was the lowest among the electrolytes used. The polarisation curves showed a current increase at potentials around 1300 mV versus saturated calomel electrode (SCE), and this increase was also dependent on the electrolyte used. The highest increase in current density was also associated to the MEM solution suggesting that this is the most aggressive electrolyte to the Ti alloy among the three tested solutions.

Failure assessments of corroded pipelines with axial defects using stress-based criteria: Numerical studies and verification analyses

Conventional procedures used to assess the integrity of corroded piping systems with axial defects generally employ simplified failure criteria based upon a plastic collapse failure mechanism incorporating the tensile properties of the pipe material. These methods establish acceptance criteria for defects based on limited experimental data for low strength structural steels which do not necessarily address specific requirements for the high grade steels currently used. For these cases, failure assessments may be overly conservative or provide significant scatter in their predictions, which lead to unnecessary repair or replacement of in-service pipelines. Motivated by these observations, this study examines the applicability of a stress-based criterion based upon plastic instability analysis to predict the failure pressure of corroded pipelines with axial defects. A central focus is to gain additional insight into effects of defect geometry and material properties on the attainment of a local limit load to support the development of stress-based burst strength criteria. The work provides an extensive body of results which lend further support to adopt failure criteria for corroded pipelines based upon ligament instability analyses. A verification study conducted on burst testing of large-diameter pipe specimens with different defect length shows the effectiveness of a stress-based criterion using local ligament instability in burst pressure predictions, even though the adopted burst criterion exhibits a potential dependence on defect geometry and possibly on material`s strain hardening capacity. Overall, the results presented here suggests that use of stress-based criteria based upon plastic instability analysis of the defect ligament is a valid engineering tool for integrity assessments of pipelines with axial corroded defects. (C) 2008 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,

Experimental and numerical investigation of dynamic heat transfer parameters in packed bed

Dynamic experiments in a nonadiabatic packed bed were carried out to evaluate the response to disturbances in wall temperature and inlet airflow rate and temperature. A two-dimensional, pseudo-homogeneous, axially dispersed plug-flow model was numerically solved and used to interpret the results. The model parameters were fitted in distinct stages: effective radial thermal conductivity (K (r)) and wall heat transfer coefficient (h (w)) were estimated from steady-state data and the characteristic packed bed time constant (tau) from transient data. A new correlation for the K (r) in packed beds of cylindrical particles was proposed. It was experimentally proved that temperature measurements using radially inserted thermocouples and a ring-shaped sensor were not distorted by heat conduction across the thermocouple or by the thermal inertia effect of the temperature sensors.

Residence time distribution in holding tubes using generalized convection model and numerical convolution for non-ideal tracer detection

A procedure is proposed for the determination of the residence time distribution (RTD) of curved tubes taking into account the non-ideal detection of the tracer. The procedure was applied to two holding tubes used for milk pasteurization in laboratory scale. Experimental data was obtained using an ionic tracer. The signal distortion caused by the detection system was considerable because of the short residence time. Four RTD models, namely axial dispersion, extended tanks in series, generalized convection and PER + CSTR association, were adjusted after convolution with the E-curve of the detection system. The generalized convection model provided the best fit because it could better represent the tail on the tracer concentration curve that is Caused by the laminar velocity profile and the recirculation regions. Adjusted model parameters were well cot-related with the now rate. (C) 2010 Elsevier Ltd. All rights reserved.

Batch Liquid-Liquid Extraction of Phenol from Aqueous Solutions

The aim of this work is the study of batch liquid-liquid extraction of phenol from aqueous solutions in a bench-scale well-mixed reactor. The influence of the ratio of phase volumes, temperature, and rotational speed on phenol removal (0.72-1.1% w/w) was investigated using methyl isobutyl ketone as an extracting solvent. For this purpose, the ratio of phase volumes were set at 0.1 and 0.2, the temperature at 10, 20, and 30 degrees C, and the rotational speed at 300, 400, and 500 rpm. A physical model based on the material balance of the phases as well as the equation of mass flux between the phases allowed the estimation of the overall coefficient of mass transfer coupled with the superficial area. Moreover, it proved to fit, satisfactorily well, the experimental data of residual phenol concentration in the organic phase versus time under all the conditions investigated.

WATER ACTIVITY OF AQUEOUS SOLUTIONS OF ETHYLENE OXIDE-PROPYLENE OXIDE BLOCK COPOLYMERS AND MALTODEXTRINS

The water activity of aqueous solutions of EO-PO block copolymers of six different molar masses and EO/PO ratios and of maltodextrins of three different molar masses was determined at 298.15 K. The results showed that these aqueous solutions present a negative deviation from Raoult`s law. The Flory-Huggins and UNIFAC excess Gibbs energy models were employed to model the experimental data. While a good agreement was obtained with the Flory-Huggins equation, discrepancies were observed when predicting the experimental behavior with the UNIFAC model. The water activities of ternary systems formed by a synthetic polymer, maltodextrin and water were also measured and used to test the predictive capability of both models.

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.

Threshold voltages of SOI MuGFETs

The multiple-gate field-effect transistor (MuGFET) is a device with a gate folded on different sides of the channel region. They are one of the most promising technological solutions to create high-performance ultra-scaled SOI CMOS. In this work, the behavior of the threshold voltage in double-gate, triple-gate and quadruple-gate SOI transistors with different channel doping concentrations is studied through three-dimensional numerical simulation. The results indicated that for double-gate transistors, one or two threshold voltages can be observed, depending on the channel doping concentration. However, in triple-gate and quadruple-gate it is possible to observe up to four threshold voltages due to the corner effect and the different doping concentration between the top and bottom of the Fin. (C) 2008 Elsevier Ltd. All rights reserved.

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.