A low-order model investigation of the analysis of gravity waves in the ensemble Kalman filter.

The behavior of the ensemble Kalman filter (EnKF) is examined in the context of a model that exhibits a nonlinear chaotic (slow) vortical mode coupled to a linear (fast) gravity wave of a given amplitude and frequency. It is shown that accurate recovery of both modes is enhanced when covariances between fast and slow normal-mode variables (which reflect the slaving relations inherent in balanced dynamics) are modeled correctly. More ensemble members are needed to recover the fast, linear gravity wave than the slow, vortical motion. Although the EnKF tends to diverge in the analysis of the gravity wave, the filter divergence is stable and does not lead to a great loss of accuracy. Consequently, provided the ensemble is large enough and observations are made that reflect both time scales, the EnKF is able to recover both time scales more accurately than optimal interpolation (OI), which uses a static error covariance matrix. For OI it is also found to be problematic to observe the state at a frequency that is a subharmonic of the gravity wave frequency, a problem that is in part overcome by the EnKF.However, error in themodeled gravity wave parameters can be detrimental to the performance of the EnKF and remove its implied advantages, suggesting that a modified algorithm or a method for accounting for model error is needed.

Impact of the modal aerosol scheme GLOMAP-mode on aerosol forcing in the Hadley Centre Global Environmental Model

The Hadley Centre Global Environmental Model (HadGEM) includes two aerosol schemes: the Coupled Large-scale Aerosol Simulator for Studies in Climate (CLASSIC), and the new Global Model of Aerosol Processes (GLOMAP-mode). GLOMAP-mode is a modal aerosol microphysics scheme that simulates not only aerosol mass but also aerosol number, represents internally-mixed particles, and includes aerosol microphysical processes such as nucleation. In this study, both schemes provide hindcast simulations of natural and anthropogenic aerosol species for the period 2000–2006. HadGEM simulations of the aerosol optical depth using GLOMAP-mode compare better than CLASSIC against a data-assimilated aerosol re-analysis and aerosol ground-based observations. Because of differences in wet deposition rates, GLOMAP-mode sulphate aerosol residence time is two days longer than CLASSIC sulphate aerosols, whereas black carbon residence time is much shorter. As a result, CLASSIC underestimates aerosol optical depths in continental regions of the Northern Hemisphere and likely overestimates absorption in remote regions. Aerosol direct and first indirect radiative forcings are computed from simulations of aerosols with emissions for the year 1850 and 2000. In 1850, GLOMAP-mode predicts lower aerosol optical depths and higher cloud droplet number concentrations than CLASSIC. Consequently, simulated clouds are much less susceptible to natural and anthropogenic aerosol changes when the microphysical scheme is used. In particular, the response of cloud condensation nuclei to an increase in dimethyl sulphide emissions becomes a factor of four smaller. The combined effect of different 1850 baselines, residence times, and abilities to affect cloud droplet number, leads to substantial differences in the aerosol forcings simulated by the two schemes. GLOMAP-mode finds a presentday direct aerosol forcing of −0.49Wm−2 on a global average, 72% stronger than the corresponding forcing from CLASSIC. This difference is compensated by changes in first indirect aerosol forcing: the forcing of −1.17Wm−2 obtained with GLOMAP-mode is 20% weaker than with CLASSIC. Results suggest that mass-based schemes such as CLASSIC lack the necessary sophistication to provide realistic input to aerosol-cloud interaction schemes. Furthermore, the importance of the 1850 baseline highlights how model skill in predicting present-day aerosol does not guarantee reliable forcing estimates. Those findings suggest that the more complex representation of aerosol processes in microphysical schemes improves the fidelity of simulated aerosol forcings.

Understanding advances in the simulation of the Madden-Julian Oscillation in a numerical weather prediction model.

The Madden-Julian Oscillation (MJO) is the dominant mode of intraseasonal variability in the Trop- ics. It can be characterised as a planetary-scale coupling between the atmospheric circulation and organised deep convection that propagates east through the equatorial Indo-Pacific region. The MJO interacts with weather and climate systems on a near-global scale and is a crucial source of predictability for weather forecasts on medium to seasonal timescales. Despite its global signifi- cance, accurately representing the MJO in numerical weather prediction (NWP) and climate models remains a challenge. This thesis focuses on the representation of the MJO in the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasting (ECMWF), a state-of-the-art NWP model. Recent modifications to the model physics in Cycle 32r3 (Cy32r3) of the IFS led to ad- vances in the simulation of the MJO; for the first time the observed amplitude of the MJO was maintained throughout the integration period. A set of hindcast experiments, which differ only in their formulation of convection, have been performed between May 2008 and April 2009 to asses the sensitivity of MJO simulation in the IFS to the Cy32r3 convective parameterization. Unique to this thesis is the attribution of the advances in MJO simulation in Cy32r3 to the mod- ified convective parameterization, specifically, the relative-humidity-dependent formulation for or- ganised deep entrainment. Increasing the sensitivity of the deep convection scheme to environmen- tal moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid-troposphere. Due to the modified precipitation-moisture relationship more moisture is able to build up which effectively preconditions the tropical atmosphere for the transition to deep convection. Results from this thesis suggest that a tropospheric moisture control on convection is key to simulating the interaction between the physics and large-scale circulation associated with the MJO.

Mode of cembranoid action on embryonic muscle acetylcholine receptor

The mechanism of eupalmerin acetate (EUAC) actions on the embryonic muscle nicotinic acetylcholine receptor (nAChR) in BC3H-1 cells was studied by using whole-cell and single-channel patch-clamp current measurements. With whole-cell currents, EUAC did not act as an agonist on this receptor. Coapplication of 30 mu M EUAC with 50 mu M, 100 N, or 500 mu M carbamoylcholine (CCh) reversibly inhibited the current amplitude, whereas, with 20 mu M CCh, current was increased above control values in the presence of EUAC. EUAC concentration curves (0.01-40 N) obtained with 100 mu M and 500 mu M CCh displayed slope coefficients, n(H), significantly smaller than one, suggesting that EUAC bound to several sites with widely differing affinities on the receptor molecule. The apparent rate of receptor desensitization in the presence of EUAC and CCh was either slower than or equal to that obtained with CCh alone. The major finding from single-channel studies was that EUAC did not affect single-channel conductance or the ability of CCh to interact with the receptor. Instead, EUAC acted by increasing the channel closing rate constant. The results are not consistent with the competitive model for EUAC inhibition, with the sequential open-channel block model, or with inhibition by increased desensitization. The data are best accounted for by a model in which EUAC acts by closed-channel block at low concentrations, by positive modulation at intermediate concentrations, and by negative allosteric modulation of the open channel at high concentrations. (c) 2007 Wiley-Liss, Inc.

Combination of Tevatron Searches for the Standard Model Higgs Boson in the W+W- Decay Mode

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

On non-ideal simple portal frame structural model: Experimental results under a non-ideal excitation

We present measurements of the non-linear oscillations of a portal frame foundation for a non-ideal motor. We consider a three-time redundant structure with two columns, clamped in their bases and a horizontal beam. An electrical unbalanced motor is mounted at mid span of the beam. Two non-linear phenomena are studied: a) mode saturation and energy transfer between modes; b) interaction between high amplitude motions of the structure and the rotation regime of a real limited power motor. The dynamic characteristics of the structure were chosen to have one-to-two internal resonance between the anti-symmetrical mode (sway motions) and the first symmetrical mode natural frequencies. As the excitation frequency reaches near resonance conditions with the 2nd natural frequency, the amplitude of this mode grows up to a certain level and then it saturates. The surplus energy pumped into the system is transferred to the sway mode, which experiences a sudden increase in its amplitude. Energy is transformed from low amplitude high frequency motion into high amplitude low frequency motion. Such a transformation is potentially dangerous.We consider the fact that real motors, such as the one used in this study, have limited power output. In this case, this energy source is said to be non-ideal, in contrast to the ideal source whose amplitude and frequency are independent of the motion of the structure. Our experimental research detected the Sommerfeld Effect: as the motor accelerates to reach near resonant conditions, a considerable part of its output energy is consumed to generate large amplitude motions of the structure and not to increase its own angular speed. For certain parameters of the system, the motor can get stuck at resonance not having enough power to reach higher rotation regimes. If some more power is available, jump phenomena may occur from near resonance to considerably higher motor speed regimes, no stable motions being possible between these two.

Pion-distribution amplitude within the instanton model

The leading-twist pion-distribution amplitude is obtained at a low normalization scale of order ρc (inverse average size of an instanton). Pion dynamics, consistent with gauge invariance and low-energy theorems, is considered within the instanton vacuum model. The results are QCD-evolved to higher momentum-transfer values and are in agreement with recent data from CLEO on the pion transition form factor. It is also shown that some previous calculations violate the axial Ward-Takahashi identity. © 2001 MAIK Nauka/Interperiodica.

Characterization, inclusion mode, phase-solubility and in vitro release studies of inclusion binary complexes with cyclodextrins and meglumine using sulfamerazine as model drug

In order to investigate the effect on the aqueous solubility and release rate of sulfamerazine (SMR) as model drug, inclusion complexes with beta-cyclodextrin (beta CD), methyl-beta-cyclodextrin (M beta CD) and hydroxypropyl-beta-cyclodextrin (HP beta CD) and a binary system with meglumine (MEG) were developed. The formation of 1: 1 inclusion complexes of SMR with the CDs and a SMR: MEG binary system in solution and in solid state was revealed by phase solubility studies (PSS), nuclear magnetic resonance (NMR), Fourier-transform infrared spectroscopy (FT-IR), thermal analysis and X-Ray diffractometry (XRD) studies. The CDs solubilization of SMR could be improved by ionization of the drug molecule through pH adjustments. The higher apparent stability constants of SMR:CDs complexes were obtained in pH 2.00, demonstrating that CDs present more affinity for the unionized drug. The best approach for SMR solubility enhancement results from the combination of MEG and pH adjustment, with a 34-fold increment and a S-max of 54.8 mg/ml. The permeability of the drug was reduced due to the presence of beta CD, M beta CD, HP beta CD and MEG when used as solubilizers. The study then suggests interesting applications of CD or MEG complexes for modulating the release rate of SMR through semipermeable membranes.

Model to describe the mode I fracture of steel fiber reinforced ultra-high performance concrete

Ultra-high performance fiber reinforced concrete (UHPFRC) has arisen from the implementation of a variety of concrete engineering and materials science concepts developed over the last century. This material offers superior strength, serviceability, and durability over its conventional counterparts. One of the most important differences for UHPFRC over other concrete materials is its ability to resist fracture through the use of randomly dispersed discontinuous fibers and improvements to the fiber-matrix bond. Of particular interest is the materials ability to achieve higher loads after first crack, as well as its high fracture toughness. In this research, a study of the fracture behavior of UHPFRC with steel fibers was conducted to look at the effect of several parameters related to the fracture behavior and to develop a fracture model based on a non-linear curve fit of the data. To determine this, a series of three-point bending tests were performed on various single edge notched prisms (SENPs). Compression tests were also performed for quality assurance. Testing was conducted on specimens of different cross-sections, span/depth (S/D) ratios, curing regimes, ages, and fiber contents. By comparing the results from prisms of different sizes this study examines the weakening mechanism due to the size effect. Furthermore, by employing the concept of fracture energy it was possible to obtain a comparison of the fracture toughness and ductility. The model was determined based on a fit to P-w fracture curves, which was cross referenced for comparability to the results. Once obtained the model was then compared to the models proposed by the AFGC in the 2003 and to the ACI 544 model for conventional fiber reinforced concretes.

Interfaces, strings, and a soft mode in the square lattice quantum dimer model

The quantum dimer model on the square lattice is a U(1) gauge theory that addresses aspects of the physics of high-Tc superconductors. Using a quantum Monte Carlo method, we show that the theory exists in a confining columnar valence bond solid phase. The interfaces separating distinct columnar phases display plaquette order, which, however, is not realized as a bulk phase. Static “electric” charges are confined by flux tubes that consist of multiple strands, each carrying a fractionalized flux ¼. A soft pseudo-Goldstone mode (which becomes exactly massless at the Rokhsar-Kivelson point) extends deep into the columnar phase, with potential implications for high-Tc physics.

Projector Model for Efficient List-Mode Reconstruction in PET Scanners with Parallel Planar Detectors

We have developed a new projector model specifically tailored for fast list-mode tomographic reconstructions in Positron emission tomography (PET) scanners with parallel planar detectors. The model provides an accurate estimation of the probability distribution of coincidence events defined by pairs of scintillating crystals. This distribution is parameterized with 2D elliptical Gaussian functions defined in planes perpendicular to the main axis of the tube of response (TOR). The parameters of these Gaussian functions have been obtained by fitting Monte Carlo simulations that include positron range, acolinearity of gamma rays, as well as detector attenuation and scatter effects. The proposed model has been applied efficiently to list-mode reconstruction algorithms. Evaluation with Monte Carlo simulations over a rotating high resolution PET scanner indicates that this model allows to obtain better recovery to noise ratio in OSEM (ordered-subsets, expectation-maximization) reconstruction, if compared to list-mode reconstruction with symmetric circular Gaussian TOR model, and histogram-based OSEM with precalculated system matrix using Monte Carlo simulated models and symmetries.

Massively parallelizable list-mode reconstruction using a Monte Carlo-based elliptical Gaussian model

Purpose: A fully three-dimensional (3D) massively parallelizable list-mode ordered-subsets expectation-maximization (LM-OSEM) reconstruction algorithm has been developed for high-resolution PET cameras. System response probabilities are calculated online from a set of parameters derived from Monte Carlo simulations. The shape of a system response for a given line of response (LOR) has been shown to be asymmetrical around the LOR. This work has been focused on the development of efficient region-search techniques to sample the system response probabilities, which are suitable for asymmetric kernel models, including elliptical Gaussian models that allow for high accuracy and high parallelization efficiency. The novel region-search scheme using variable kernel models is applied in the proposed PET reconstruction algorithm. Methods: A novel region-search technique has been used to sample the probability density function in correspondence with a small dynamic subset of the field of view that constitutes the region of response (ROR). The ROR is identified around the LOR by searching for any voxel within a dynamically calculated contour. The contour condition is currently defined as a fixed threshold over the posterior probability, and arbitrary kernel models can be applied using a numerical approach. The processing of the LORs is distributed in batches among the available computing devices, then, individual LORs are processed within different processing units. In this way, both multicore and multiple many-core processing units can be efficiently exploited. Tests have been conducted with probability models that take into account the noncolinearity, positron range, and crystal penetration effects, that produced tubes of response with varying elliptical sections whose axes were a function of the crystal's thickness and angle of incidence of the given LOR. The algorithm treats the probability model as a 3D scalar field defined within a reference system aligned with the ideal LOR. Results: This new technique provides superior image quality in terms of signal-to-noise ratio as compared with the histogram-mode method based on precomputed system matrices available for a commercial small animal scanner. Reconstruction times can be kept low with the use of multicore, many-core architectures, including multiple graphic processing units. Conclusions: A highly parallelizable LM reconstruction method has been proposed based on Monte Carlo simulations and new parallelization techniques aimed at improving the reconstruction speed and the image signal-to-noise of a given OSEM algorithm. The method has been validated using simulated and real phantoms. A special advantage of the new method is the possibility of defining dynamically the cut-off threshold over the calculated probabilities thus allowing for a direct control on the trade-off between speed and quality during the reconstruction.

Spatial linear global instability analysis of the HIFiRE-5 elliptic cone model flow

The linear instability of the three-dimensional boundary-layer over the HIFiRE-5 flight test geometry, i.e. a rounded-tip 2:1 elliptic cone, at Mach 7, has been analyzed through spatial BiGlobal analysis, in a effort to understand transition and accurately predict local heat loads on next-generation ight vehicles. The results at an intermediate axial section of the cone, Re x = 8x10 5, show three different families of spatially amplied linear global modes, the attachment-line and cross- ow modes known from earlier analyses, and a new global mode, peaking in the vicinity of the minor axis of the cone, termed \center-line mode". We discover that a sequence of symmetric and anti-symmetric centerline modes exist and, for the basic ow at hand, are maximally amplied around F* = 130kHz. The wavenumbers and spatial distribution of amplitude functions of the centerline modes are documented