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.

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.

Light scattering in polycrystalline alumina with bi-dimensionally large surface grains

Alumina ceramics with high in-line transmittance at 0.5-1.0 mm-thickness were prepared with different doping additives by sintering at 1850 degrees C in vacuum for 1-8 h. Depending on the additive contents and sintering variables bi-dimensionally large surface grains, caused by surface evaporation of MgO, had grown parallel to the surface with similar to 100 mu m thickness and lateral sizes up to the millimeter range. The abnormal grain-growth process also resulted in the formation of pores entrapped inside the large surface grains within a narrow zone at 10-20 mu m distance from the surface. The fraction of these pores is thickness-invariant. Scattering factors associated to the pores entrapped inside the bi-dimensionally large surface grains, second-phase particles, grain-boundaries, and microstructural surface defects are derived from the results of in-line transmission (at 600 nm) and are used together with microstructural characteristics to explain the light transmittance in these materials. (C) 2008 Elsevier Ltd. All rights reserved.

Rotation-discriminating template matching based on Fourier coefficients of radial projections with robustness to scaling and partial occlusion

We consider brightness/contrast-invariant and rotation-discriminating template matching that searches an image to analyze A for a query image Q We propose to use the complex coefficients of the discrete Fourier transform of the radial projections to compute new rotation-invariant local features. These coefficients can be efficiently obtained via FFT. We classify templates in ""stable"" and ""unstable"" ones and argue that any local feature-based template matching may fail to find unstable templates. We extract several stable sub-templates of Q and find them in A by comparing the features. The matchings of the sub-templates are combined using the Hough transform. As the features of A are computed only once, the algorithm can find quickly many different sub-templates in A, and it is Suitable for finding many query images in A, multi-scale searching and partial occlusion-robust template matching. (C) 2009 Elsevier Ltd. All rights reserved.

A Separation Principle for the Continuous-Time LQ-Problem With Markovian Jump Parameters

In this paper, we devise a separation principle for the finite horizon quadratic optimal control problem of continuous-time Markovian jump linear systems driven by a Wiener process and with partial observations. We assume that the output variable and the jump parameters are available to the controller. It is desired to design a dynamic Markovian jump controller such that the closed loop system minimizes the quadratic functional cost of the system over a finite horizon period of time. As in the case with no jumps, we show that an optimal controller can be obtained from two coupled Riccati differential equations, one associated to the optimal control problem when the state variable is available, and the other one associated to the optimal filtering problem. This is a separation principle for the finite horizon quadratic optimal control problem for continuous-time Markovian jump linear systems. For the case in which the matrices are all time-invariant we analyze the asymptotic behavior of the solution of the derived interconnected Riccati differential equations to the solution of the associated set of coupled algebraic Riccati equations as well as the mean square stabilizing property of this limiting solution. When there is only one mode of operation our results coincide with the traditional ones for the LQG control of continuous-time linear systems.

Flatness-based control of a single qubit gate

This work considers the open-loop control problem of steering a two-level quantum system from any initial to any final condition. The model of this system evolves on the state space X = SU(2), having two inputs that correspond to the complex amplitude of a resonant laser field. A symmetry preserving flat output is constructed using a fully geometric construction and quaternion computations. Simulation results of this flatness-based open-loop control are provided.

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes. (C) 2008 Elsevier Ltd. All rights reserved.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

One-dimensional point interaction with three complex parameters

For a pair of non-Hermitian Hamiltonian H and its Hermitian adjoint H(dagger), there are situations in which their eigenfunctions form a biorthogonal system. We illustrate such a situation by means of a one-particle system with a one-dimensional point interaction in the form of the Fermi pseudo-potential. The interaction consists of three terms with three strength parameters g(i) (i = 1, 2 and 3), which are all complex. This complex point interaction is neither Hermitian nor PT-invariant in general. The S-matrix for the transmission reflection problem constructed with H (or with H(dagger)) in the usual manner is not unitary, but it conforms to the pseudo-unitarity that we define. The pseudounitarity is closely related to the biorthogonality of the eigenfunctions. The eigenvalue spectrum of H with the complex interaction is generally complex but there are cases where the spectrum is real. In such a case H and H(dagger) form a pseudo-Hermitian pair.

Hypomania: a transcultural perspective

This study examined the transcultural robustness of a screening instrument for hypomania, the Hypomania Checklist-32. first revised version(HCL-32 R1) It was carried out in 2606 patients from twelve countries in five geographic regions (Northern, Southern and Eastern Europe, South America and East Asia) In addition, GAMIAN Europe contributed data from its members Exploratory and confirmatory factor analyses were used to examine the transregional stability of the measurement properties of the HCL-32 R 1, including the influence of sex and age as covariates A cross cultures, a two-factor structure was confirmed the first factor (F 1) reflected the more positive aspects of hypomania (being more active, elated, self-confident, and cognitively, enhanced), the second factor (F2) reflected the more negative aspects (being irritable, impulsive, careless. more substance use). The measurement properties of the HCL-32 R1 were largely invariant across cultures Only few items showed transcultural differences in their relation to hypomania as measured by the test F2 was higher among men and in more severe manic syndromes, F1 was highest in North and East Europe and lowest in South America Vie scores decreased slightly with age The frequency of the 32 items showed remarkable similarities across geographic areas, with two exceptions. South Europeans had lower symptom frequencies in general and East Europeans higher rates of substance use These findings support the international applicability of the HCL-32 R1 as a screening instrument for hypomania

Lower cytokine secretion ex vivo by natural killer T cells in HIV-infected individuals is associated with higher CD161 expression

Objective: Natural killer T (NKT) cells are efficiently targeted by HIV and severely reduced in numbers in the circulation of infected individuals. The functional capacity of the remaining NKT cells in HIV-infected individuals is poorly characterized. This study measured NKT cell cytokine production directly ex vivo and compared these responses with both the disease status and NKT subset distribution of individual patients. Methods: NKT cell frequencies, subsets, and ex-vivo effector functions were measured in the peripheral blood mononuclear cells of HIV-infected patients and healthy controls by flow cytometry. We measured cytokines from NKT cells after stimulation with either a-galactosyl ceramide-loaded CD1d dimers (DimerX-alpha GalCer) or phorbol myristate acetate and ionomycin. Results: The frequencies of NKT cells secreting interferon-gamma and tumor necrosis factor-alpha were significantly lower in HIV-infected patients than healthy controls after DimerX-alpha GalCer treatment, but responses were similar after treatment with phorbol myristate acetate and ionomycin. The magnitude of the interferon-gamma response to DimerX-alpha GalCer correlated inversely with the number of years of infection. Both interferon-gamma and tumor necrosis factor-alpha production in response to DimerX-alpha GalCer correlated inversely with CD161 expression. Conclusion: The ex-vivo Th1 responses of circulating NKT cells to CD1d-glycolipid complexes are impaired in HIV-infected patients. NKT cell functions may be progressively lost over time in HIV infection, and CD161 is implicated in the regulation of NKT cell responsiveness. (C) 2009 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins

Study of the discovery potential of light custodians in a clean decay channel at the LHC

Models of warped extra dimensions with custodial symmetry usually predict the existence of a light Kaluza-Klein fermion arising as a partner of the right-handed top quark, sometimes called light custodians which we will denote (b) over tilde (R). The production of these particles at the LHC can give rise to multi-W events which could be observed in same-sign dilepton channels, but its mass reconstruction is challenging. In this paper we study the possibility of finding a signal for the pair production of this new particle at the LHC focusing on a rarer, but cleaner decay mode of a light custodian into a Z boson and a b-quark. In this mode it would be possible to reconstruct the light custodian mass. In addition to the dominant standard model QCD production processes, we include the contribution of a Kaluza-Klein gluon first mode. We find that (b) over tilde (R) stands out from the background as a peak in the bZ invariant mass. However, when taking into account only the electronic and muonic decay modes of the Z boson and b-tagging efficiencies, the LHC will have access only to the very light range of masses, m((b) over tilde) = O(500) GeV.