An electromechanical finite element model for piezoelectric energy harvester plates

Vibration-based energy harvesting has been investigated by several researchers over the last decade. The goal in this research field is to power small electronic components by converting the waste vibration energy available in their environment into electrical energy. Recent literature shows that piezoelectric transduction has received the most attention for vibration-to-electricity conversion. In practice, cantilevered beams and plates with piezoceramic layers are employed as piezoelectric energy harvesters. The existing piezoelectric energy harvester models are beam-type lumped parameter, approximate distributed parameter and analytical distributed parameter solutions. However, aspect ratios of piezoelectric energy harvesters in several cases are plate-like and predicting the power output to general (symmetric and asymmetric) excitations requires a plate-type formulation which has not been covered in the energy harvesting literature. In this paper. an electromechanically coupled finite element (FE) plate model is presented for predicting the electrical power output of piezoelectric energy harvester plates. Generalized Hamilton`s principle for electroelastic bodies is reviewed and the FE model is derived based on the Kirchhoff plate assumptions as typical piezoelectric energy harvesters are thin structures. Presence of conductive electrodes is taken into account in the FE model. The predictions of the FE model are verified against the analytical solution for a unimorph cantilever and then against the experimental and analytical results of a bimorph cantilever with a tip mass reported in the literature. Finally, an optimization problem is solved where the aluminum wing spar of an unmanned air vehicle (UAV) is modified to obtain a generator spar by embedding piezoceramics for the maximum electrical power without exceeding a prescribed mass addition limit. (C) 2009 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Vehicle Dynamic Safety in Measured Rough Pavement

Dynamic vehicle behavior is used to identify safe traffic speed limits. The proposed methodology is based on the vehicle vertical wheel contact force response excited by measured pavement irregularities on the frequency domain. A quarter-car model is used to identify vehicle dynamic behavior. The vertical elevation of an unpaved road surface has been measured. The roughness spectral density is quantified as ISO Level C. Calculations for the vehicle inertance function were derived by using the vertical contact force transfer function weighed by the pavement spectral density roughness function in the frequency domain. The statistical contact load variation is obtained from the vehicle inertance density function integration. The vehicle safety behavior concept is based on its handling ability properties. The ability to generate tangential forces on the wheel/road contact interface is the key to vehicle handling. This ability is related to tire/pavement contact forces. A contribution to establish a traffic safety speed limit is obtained from the likelihood of the loss of driveability. The results show that at speeds faster than 25 km/h the likelihood of tire contact loss is possible when traveling on the measured road type. DOI: 10.1061/(ASCE)TE.19435436.0000216. (C) 2011 American Society of Civil Engineers.

Metodologia para análise da confiabilidade estrutural de escavações em rocha.

A methodology for rock-excavation structural-reliability analysis that uses Distinct Element Method numerical models is presented. The methodology solves the problem of the conventional numerical models that supply only punctual results and use fixed input parameters, without considering its statistical errors. The analysis of rock-excavation stability must consider uncertainties from geological variability, from uncertainty in the choice of mechanical behaviour hypothesis, and from uncertainties in parameters adopted in numerical model construction. These uncertainties can be analyzed in simple deterministic models, but a new methodology was developed for numerical models with results of several natures. The methodology is based on Monte Carlo simulations and uses principles of Paraconsistent Logic. It will be presented in the analysis of a final slope of a large-dimensioned surface mine.

Mechanical Blocking Mechanism for the Columnar to Equiaxed Transition

Mechanical blocking of the columnar front during the columnar to equiaxed transition (CET) is studied by quantitatively comparing the CET positions obtained with one stochastic model and two deterministic models for the unidirectional solidification of an Al-7 (wt pct) Si alloy. One of the deterministic models is based on the solutal blocking of the columnar front, whereas the other model is based on the mechanical blocking. The solutal-blocking model and the mechanical-blocking model with the traditional blocking fraction of 0.49 give columnar zones larger than those predicted with the stochastic model. When a blocking fraction of 0.2 is adopted, however, the agreement is very good for a range of nucleation undercoolings and number density of equiaxed grains. Therefore, changing the mechanical-blocking fraction in deterministic models from 0.49 to 0.2 seems to model more accurately the mechanical-blocking process that can lead to the CET.

Photoinduced Degradation of the Herbicide Clomazone Model Reactions for Natural and Technical Systems

The photodegradation of the herbicide clomazone in the presence of S(2)O(8)(2-) or of humic substances of different origin was investigated. A value of (9.4 +/- 0.4) x 10(8) m(-1) s(-1) was measured for the bimolecular rate constant for the reaction of sulfate radicals with clomazone in flash-photolysis experiments. Steady state photolysis of peroxydisulfate, leading to the formation of the sulfate radicals, in the presence of clomazone was shown to be an efficient photodegradation method of the herbicide. This is a relevant result regarding the in situ chemical oxidation procedures involving peroxydisulfate as the oxidant. The main reaction products are 2-chlorobenzylalcohol and 2-chlorobenzaldehyde. The degradation kinetics of clomazone was also studied under steady state conditions induced by photolysis of Aldrich humic acid or a vermicompost extract (VCE). The results indicate that singlet oxygen is the main species responsible for clomazone degradation. The quantum yield of O(2)(a(1)Delta(g)) generation (lambda = 400 nm) for the VCE in D(2)O, Phi(Delta) = (1.3 +/- 0.1) x 10(-3), was determined by measuring the O(2)(a(1)Delta(g)) phosphorescence at 1270 nm. The value of the overall quenching constant of O(2)(a(1)Delta(g)) by clomazone was found to be (5.7 +/- 0.3) x 10(7) m(-1) s(-1) in D(2)O. The bimolecular rate constant for the reaction of clomazone with singlet oxygen was k(r) = (5.4 +/- 0.1) x 10(7) m(-1) s(-1), which means that the quenching process is mainly reactive.

Effects of mechanical properties, residual stress and indenter tip geometry on instrumented indentation data in thin films

In this work, an axisymmetric two-dimensional finite element model was developed to simulate instrumented indentation testing of thin ceramic films deposited onto hard steel substrates. The level of film residual stress (sigma(r)), the film elastic modulus (E) and the film work hardening exponent (n) were varied to analyze their effects on indentation data. These numerical results were used to analyze experimental data that were obtained with titanium nitride coated specimens, in which the substrate bias applied during deposition was modified to obtain films with different levels of sigma(r). Good qualitative correlation was obtained when numerical and experimental results were compared, as long as all film properties are considered in the analyses, and not only sigma(r). The numerical analyses were also used to further understand the effect of sigma(r) on the mechanical properties calculated based on instrumented indentation data. In this case, the hardness values obtained based on real or calculated contact areas are similar only when sink-in occurs, i.e. with high n or high ratio VIE, where Y is the yield strength of the film. In an additional analysis, four ratios (R/h(max)) between indenter tip radius and maximum penetration depth were simulated to analyze the combined effects of R and sigma(r) on the indentation load-displacement curves. In this case, or did not significantly affect the load curve exponent, which was affected only by the indenter tip radius. On the other hand, the proportional curvature coefficient was significantly affected by sigma(r) and n. (C) 2010 Elsevier B.V. All rights reserved.

MULTIVARIATE STATISTICAL METHOD IN OPTIMIZATION OF PROTEIN EXTRACTION FROM SOYBEAN FLOUR WITH DAIRY WHEY

The aim objective of this project was to evaluate the protein extraction of soybean flour in dairy whey, by the multivariate statistical method with 2(3) experiments. Influence of three variables were considered: temperature, pH and percentage of sodium chloride against the process specific variable ( percentage of protein extraction). It was observed that, during the protein extraction against time and temperature, the treatments at 80 degrees C for 2h presented great values of total protein (5.99%). The increasing for the percentage of protein extraction was major according to the heating time. Therefore, the maximum point from the function that represents the protein extraction was analysed by factorial experiment 2(3). By the results, it was noted that all the variables were important to extraction. After the statistical analyses, was observed that the parameters as pH, temperature, and percentage of sodium chloride, did not sufficient for the extraction process, since did not possible to obtain the inflection point from mathematical function, however, by the other hand, the mathematical model was significant, as well as, predictive.

The generalized log-gamma mixture model with covariates: local influence and residual analysis

In a sample of censored survival times, the presence of an immune proportion of individuals who are not subject to death, failure or relapse, may be indicated by a relatively high number of individuals with large censored survival times. In this paper the generalized log-gamma model is modified for the possibility that long-term survivors may be present in the data. The model attempts to separately estimate the effects of covariates on the surviving fraction, that is, the proportion of the population for which the event never occurs. The logistic function is used for the regression model of the surviving fraction. Inference for the model parameters is considered via maximum likelihood. Some influence methods, such as the local influence and total local influence of an individual are derived, analyzed and discussed. Finally, a data set from the medical area is analyzed under the log-gamma generalized mixture model. A residual analysis is performed in order to select an appropriate model.

A split-pot experiment with sorghum to test a root water uptake partitioning model

Correct modeling of root water uptake partitioning over depth is an important issue in hydrological and crop growth models. Recently a physically based model to describe root water uptake was developed at single root scale and upscaled to the root system scale considering a homogeneous distribution of roots per soil layer. Root water uptake partitioning is calculated over soil layers or compartments as a function of respective soil hydraulic conditions, specifically the soil matric flux potential, root characteristics and a root system efficiency factor to compensate for within-layer root system heterogeneities. The performance of this model was tested in an experiment performed in two-compartment split-pot lysimeters with sorghum plants. The compartments were submitted to different irrigation cycles resulting in contrasting water contents over time. The root system efficiency factor was determined to be about 0.05. Release of water from roots to soil was predicted and observed on several occasions during the experiment; however, model predictions suggested root water release to occur more often and at a higher rate than observed. This may be due to not considering internal root system resistances, thus overestimating the ease with which roots can act as conductors of water. Excluding these erroneous predictions from the dataset, statistical indices show model performance to be of good quality.

Fragmentation studies and electrospray ionization mass spectrometry of lapachol: protonated, deprotonated and cationized species

Electrospray ionization mass spectrometric analysis of lapachol (2-hydroxy-3-(3-methy1-2-butenyl)-1,4-naphthoquinone) was accomplished in order to elucidate the gas-phase dissociation reactions of this important biologically active natural product. The occurrence of protonated and cationized species in the positive mode and of deprotonated species in the negative mode was explored by means of collision-induced dissociation (CID) experiments. For the protonated molecule, the H(2)O and C(4)H(8) losses occur by two competitive channels. For the deprotonated molecule, the even-electron rule is not conserved, and the radicalar species are eliminated by formation of distonic anions. The fragmentation mechanism for each ion was suggested on the basis of computational thermochemistry. Atomic charges, relative energies, and frontier orbitals were employed aiming at a better understanding of the gas-phase reactivity of lapachol. Potential energy surfaces for fragmentation reactions were obtained by the B3LYP/6-31+G(d,p) model. Copyright (C) 2010 John Wiley & Sons, Ltd.

Quantum-chemical, NMR and X-ray diffraction studies on (+/-)-1-[3,4-(methylenedioxy)phenyl]-2-methylaminopropane

Time-averaged conformations of (+/-)-1-[3,4-(methylenedioxy)phenyl]-2-methylaminopropane hydrochloride (MDMA, ""ecstasy"") in D(2)O, and of its free base and trifluoroacetate in CDCl(3), were deduced from their (1)H NMR spectra and used to calculate their conformer distribution. Their rotational potential energy surface (PES) was calculated at the RHF/6-31G(d,p), 133LYP/6-31G(d,p), B3LYP/cc-pVDZ and AM1 levels. Solvent effects were evaluated using the polarizable continuum model. The NMR and theoretical studies showed that, in the free base, the N-methyl group and the ring are preferentially trans. This preference is stronger in the salts and corresponds to the X-ray structure of the hydrochloride. However, the energy barriers separating these forms are very low. The X-ray diffraction crystal structures of the anhydrous salt and its monohydrate differed mainly in the trans or cis relationship of the N-methyl group to the a-methyl, although these two forms interconvert freely in solution. (C) 2007 Elsevier Inc. All rights reserved.

Loss of Ergodicity in the Transition from Annealed to Quenched Disorder in a Finite Kinetic Ising Model

We consider a kinetic Ising model which represents a generic agent-based model for various types of socio-economic systems. We study the case of a finite (and not necessarily large) number of agents N as well as the asymptotic case when the number of agents tends to infinity. The main ingredient are individual decision thresholds which are either fixed over time (corresponding to quenched disorder in the Ising model, leading to nonlinear deterministic dynamics which are generically non-ergodic) or which may change randomly over time (corresponding to annealed disorder, leading to ergodic dynamics). We address the question how increasing the strength of annealed disorder relative to quenched disorder drives the system from non-ergodic behavior to ergodicity. Mathematically rigorous analysis provides an explicit and detailed picture for arbitrary realizations of the quenched initial thresholds, revealing an intriguing ""jumpy"" transition from non-ergodicity with many absorbing sets to ergodicity. For large N we find a critical strength of annealed randomness, above which the system becomes asymptotically ergodic. Our theoretical results suggests how to drive a system from an undesired socio-economic equilibrium (e. g. high level of corruption) to a desirable one (low level of corruption).

A Comparison Between Pulse Pressure Variation and Right End Diastolic Volume Index as Guides to Resuscitation in a Model of Hemorrhagic Shock in Pigs

Background: Different hemodynamic parameters including static indicators of cardiac preload as right ventricular end-diastolic volume index (RVEDVI) and dynamic parameters as pulse pressure variation (PPV) have been used in the decision-making process regarding volume expansion in critically ill patients. The objective of this study was to compare fluid resuscitation guided by either PPV or RVEDVI after experimentally induced hemorrhagic shock. Methods: Twenty-six anesthetized and mechanically ventilated pigs were allocated into control (group I), PPV (group II), or RVEDVI (group III) group. Hemorrhagic shock was induced by blood withdrawal to target mean arterial pressure of 40 mm Hg, maintained for 60 minutes. Parameters were measured at baseline, time of shock, 60 minutes after shock, immediately after resuscitation with hydroxyethyl starch 6% (130/0.4), 1 hour and 2 hours thereafter. The endpoint of fluid resuscitation was determined as the baseline values of PPV and RVEDVI. Statistical analysis of data was based on analysis of variance for repeated measures followed by the Bonferroni test (p < 0.05). Results: Volume and time to resuscitation were higher in group III than in group II (group III = 1,305 +/- 331 mL and group II = 965 +/- 245 mL, p < 0.05; and group III = 24.8 +/- 4.7 minutes and group II = 8.8 +/- 1.3 minutes, p < 0.05, respectively). All static and dynamic parameters and biomarkers of tissue oxygenation were affected by hemorrhagic shock and nearly all parameters were restored after resuscitation in both groups. Conclusion: In the proposed model of hemorrhagic shock, resuscitation to the established endpoints was achieved within a smaller amount of time and with less volume when guided by PPV than when guided by pulmonary artery catheter-derived RVEDVI.