Bending of stochastic Kirchhoff plates on Winkler foundations via the Galerkin method and the Askey-Wiener scheme

In this paper, the method of Galerkin and the Askey-Wiener scheme are used to obtain approximate solutions to the stochastic displacement response of Kirchhoff plates with uncertain parameters. Theoretical and numerical results are presented. The Lax-Milgram lemma is used to express the conditions for existence and uniqueness of the solution. Uncertainties in plate and foundation stiffness are modeled by respecting these conditions, hence using Legendre polynomials indexed in uniform random variables. The space of approximate solutions is built using results of density between the space of continuous functions and Sobolev spaces. Approximate Galerkin solutions are compared with results of Monte Carlo simulation, in terms of first and second order moments and in terms of histograms of the displacement response. Numerical results for two example problems show very fast convergence to the exact solution, at excellent accuracies. The Askey-Wiener Galerkin scheme developed herein is able to reproduce the histogram of the displacement response. The scheme is shown to be a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2009 Elsevier Ltd. All rights reserved.

Chaos-Galerkin solution of stochastic Timoshenko bending problems

This paper presents an accurate and efficient solution for the random transverse and angular displacement fields of uncertain Timoshenko beams. Approximate, numerical solutions are obtained using the Galerkin method and chaos polynomials. The Chaos-Galerkin scheme is constructed by respecting the theoretical conditions for existence and uniqueness of the solution. Numerical results show fast convergence to the exact solution, at excellent accuracies. The developed Chaos-Galerkin scheme accurately approximates the complete cumulative distribution function of the displacement responses. The Chaos-Galerkin scheme developed herein is a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2011 Elsevier Ltd. All rights reserved.

Timoshenko versus Euler beam theory: Pitfalls of a deterministic approach

The selection criteria for Euler-Bernoulli or Timoshenko beam theories are generally given by means of some deterministic rule involving beam dimensions. The Euler-Bernoulli beam theory is used to model the behavior of flexure-dominated (or ""long"") beams. The Timoshenko theory applies for shear-dominated (or ""short"") beams. In the mid-length range, both theories should be equivalent, and some agreement between them would be expected. Indeed, it is shown in the paper that, for some mid-length beams, the deterministic displacement responses for the two theories agrees very well. However, the article points out that the behavior of the two beam models is radically different in terms of uncertainty propagation. In the paper, some beam parameters are modeled as parameterized stochastic processes. The two formulations are implemented and solved via a Monte Carlo-Galerkin scheme. It is shown that, for uncertain elasticity modulus, propagation of uncertainty to the displacement response is much larger for Timoshenko beams than for Euler-Bernoulli beams. On the other hand, propagation of the uncertainty for random beam height is much larger for Euler beam displacements. Hence, any reliability or risk analysis becomes completely dependent on the beam theory employed. The authors believe this is not widely acknowledged by the structural safety or stochastic mechanics communities. (C) 2010 Elsevier Ltd. All rights reserved.

Galerkin Solution of Stochastic Beam Bending on Winkler Foundations

In this paper, the Askey-Wiener scheme and the Galerkin method are used to obtain approximate solutions to stochastic beam bending on Winkler foundation. The study addresses Euler-Bernoulli beams with uncertainty in the bending stiffness modulus and in the stiffness of the foundation. Uncertainties are represented by parameterized stochastic processes. The random behavior of beam response is modeled using the Askey-Wiener scheme. One contribution of the paper is a sketch of proof of existence and uniqueness of the solution to problems involving fourth order operators applied to random fields. From the approximate Galerkin solution, expected value and variance of beam displacement responses are derived, and compared with corresponding estimates obtained via Monte Carlo simulation. Results show very fast convergence and excellent accuracies in comparison to Monte Carlo simulation. The Askey-Wiener Galerkin scheme presented herein is shown to be a theoretically solid and numerically efficient method for the solution of stochastic problems in engineering.

FOREIGN INVESTMENT AND THE BRAZILIAN REAL ESTATE MARKET

This article discusses the main aspects of the Brazilian real estate market in order to illustrate if it would be attractive for a typical American real estate investor to buy office-building portfolios in Brazil. The article emphasizes: [i] - the regulatory frontiers, comparing investment securitization, using a typical American REIT structure, with the Brazilian solution, using the Fundo de Investimento Imobiliario - FII; [ii] - the investment quality attributes in the Brazilian market, using an office building prototype, and [iii] - the comparison of [risk vs. yield] generated by an investment in the Brazilian market, using a FII, benchmarked against an existing REIT (OFFICE SUB-SECTOR) in the USA market. We conclude that investing dollars exchanged for Reais [the Brazilian currency] in a FII with a triple A office-building portfolio in the Sao Paulo marketplace will yield an annual income and a premium return above an American REIT investment. The highly aggressive scenario, along with the strong persistent exchange rate detachment to the IGP-M variations, plus instabilities affecting the generation of income, and even if we adopt a 300-point margin for the Brazil-Risk level, demonstrates that an investment opportunity in the Brazilian market, in the segment we have analyzed, outperforms an equivalent investment in the American market.

Reconfiguration of distribution networks to minimize loss and disruption costs using genetic algorithms

In this paper a computational implementation of an evolutionary algorithm (EA) is shown in order to tackle the problem of reconfiguring radial distribution systems. The developed module considers power quality indices such as long duration interruptions and customer process disruptions due to voltage sags, by using the Monte Carlo simulation method. Power quality costs are modeled into the mathematical problem formulation, which are added to the cost of network losses. As for the EA codification proposed, a decimal representation is used. The EA operators, namely selection, recombination and mutation, which are considered for the reconfiguration algorithm, are herein analyzed. A number of selection procedures are analyzed, namely tournament, elitism and a mixed technique using both elitism and tournament. The recombination operator was developed by considering a chromosome structure representation that maps the network branches and system radiality, and another structure that takes into account the network topology and feasibility of network operation to exchange genetic material. The topologies regarding the initial population are randomly produced so as radial configurations are produced through the Prim and Kruskal algorithms that rapidly build minimum spanning trees. (C) 2009 Elsevier B.V. All rights reserved.

Fuzzy risk index for power transformer failures due to external short-circuits

A novel methodology to assess the risk of power transformer failures caused by external faults, such as short-circuit, taking the paper insulation condition into account, is presented. The risk index is obtained by contrasting the insulation paper condition with the probability that the transformer withstands the short-circuit current flowing along the winding during an external fault. In order to assess the risk, this probability and the value of the degree of polymerization of the insulating paper are regarded as inputs of a type-2 fuzzy logic system (T2-FLS), which computes the fuzzy risk level. A Monte Carlo simulation has been used to find the survival function of the currents flowing through the transformer winding during a single-phase or a three-phase short-circuit. The Roy Billinton Test System and a real power system have been used to test the results. (C) 2008 Elsevier B.V. All rights reserved.

Evaluating the efficiency of fractional integration parameter estimators

This article deals with the efficiency of fractional integration parameter estimators. This study was based on Monte Carlo experiments involving simulated stochastic processes with integration orders in the range]-1,1[. The evaluated estimation methods were classified into two groups: heuristics and semiparametric/maximum likelihood (ML). The study revealed that the comparative efficiency of the estimators, measured by the lesser mean squared error, depends on the stationary/non-stationary and persistency/anti-persistency conditions of the series. The ML estimator was shown to be superior for stationary persistent processes; the wavelet spectrum-based estimators were better for non-stationary mean reversible and invertible anti-persistent processes; the weighted periodogram-based estimator was shown to be superior for non-invertible anti-persistent processes.

Molecular Characterization, Distribution, and Dynamics of Hepatitis C Virus Genotypes in Blood Donors in Colombia

Hepatitis C virus (HCV) is a frequent cause of acute and chronic hepatitis and a leading cause for cirrhosis of the liver and hepatocellular carcinoma. HCV is classified in six major genotypes and more than 70 subtypes. In Colombian blood banks, serum samples were tested for anti-HCV antibodies using a third-generation ELISA. The aim of this study was to characterize the viral sequences in plasma of 184 volunteer blood donors who attended the ""Banco Nacional de Sangre de la Cruz Roja Colombiana,`` Bogota, Colombia. Three different HCV genomic regions were amplified by nested PCR. The first of these was a segment of 180 bp of the 5`UTR region to confirm the previous diagnosis by ELISA. From those that were positive to the 5`UTR region, two further segments were amplified for genotyping and subtyping by phylogenetic analysis: a segment of 380 bp from the NS5B region; and a segment of 391 bp from the E1 region. The distribution of HCV subtypes was: 1b (82.8%), 1a (5.7%), 2a (5.7%), 2b (2.8%), and 3a (2.8%). By applying Bayesian Markov chain Monte Carlo simulation, it was estimated that HCV-1b was introduced into Bogota around 1950. Also, this subtype spread at an exponential rate between about 1970 to about 1990, after which transmission of HCV was reduced by anti-HCV testing of this population. Among Colombian blood donors, HCV genotype 1b is the most frequent genotype, especially in large urban conglomerates such as Bogota, as is the case in other South American countries. J. Med. Virol. 82: 1889-1898, 2010. (C) 2010 Wiley-Liss, Inc.

A portable Compton spectrometer for clinical X-ray beams in the energy range 20-150 keV

Primary beam spectra were obtained for an X-ray industrial equipment (40-150 kV), and for a clinical mammography apparatus (25-35 kV) from beams scattered at angles close to 90 degrees, measured with a CdTe Compton spectrometer. Actual scattering angles were determined from the Compton energy shift of characteristic X-rays or spectra end-point energy. Evaluated contribution of coherent scattering amounts to more than 15% of fluence in mammographic beams. This technique can be used in clinical environments. (C) 2010 Elsevier Ltd. All rights reserved.

Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Electronic spectroscopy of biomolecules in solution: fluorescein dianion in water

A combined and sequential Monte Carlo-quantum mechanics methodology is used to describe the electronic absorption spectrum of the fluorescein dianion in water. Different sets of 100 statistically relevant configurations composed of the solute and several solvent molecules are sampled from the Monte Carlo simulation for a posteriori quantum mechanical calculations of the spectra. In the largest case the configurations are composed of fluorescein and 90 explicit water molecules embedded in the electrostatic field of all remaining water molecules within a distance of 11.3 angstrom. These configurations include 305 atoms and 842 valence electrons, justifying the use of a semi-empirical approach. The electronic spectrum is then calculated using the INDO/CIS method. The solvatochromic shift of fluorescein in water, compared with in isolation, is calculated using the discrete and explicit solvent models. The use of electrostatically embedded explicit water molecules, in INDO/CIS calculations, gives a good description of the spectral shift of the fluorescein dianion in aqueous environment. The results are verified to converge both statistically and with respect to the number of explicit solvent molecules used.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.