Simulation of the Unsteady Gas Flow through a Three-way Automotive Catalyst

This paper describes a model of a 1.8-litre four-cylinder four-stroke gasoline engine fitted with a close-coupled three-way catalyst (TWC). Designed to meet EURO 3 emissions standards, the engine includes some advanced emission control features in addition to the TWC, namely: variable valve timing (VVT), swirl control plates, and exhaust gas recirculation (EGR). Gas flow is treated as one-dimensional (1D) and unsteady in the engine ducting and in the catalyst. Reflection and transmission of pressure waves at the boundaries of the catalyst monolith are modelled. In-cylinder combustion is represented by a two-zone burn model with dissociation and reaction kinetics. A single Wiebe analysis of measured in-cylinder pressure data is used to determine the mass fraction burned as a function of crank angle (CA) at each engine speed. Measured data from steady-state dynamometer tests are presented for operation at wide open throttle (WOT) over a range of engine speeds. These results include CA-resolved traces of pressure at various locations throughout the engine together with cycle-averaged traces of gas composition entering the catalyst as indicated by a fast-response emissions analyser. Simulated engine performance and pressure wave action throughout the engine are well validated by the measured data.

Modelling of Unsteady Gas-Dynamic Flow in a Pipe and at its Exit Using CFD

Traditionally the simulation of the thermodynamic aspects of the internal combustion engine has been undertaken using one-dimensional gas-dynamic models to represent the intake and exhaust systems. CFD analysis of engines has been restricted to modelling of in-cylinder flow structures. With the increasing accessibility of CFD software it is now worth considering its use for complete gas-dynamic engine simulation. This paper appraises the accuracy of various CFD models in comparison to a 1D gas-dynamic simulation. All of the models are compared to experimental data acquired on an apparatus that generates a single gas-dynamic pressure wave. The progress of the wave along a constant area pipe and its subsequent reflection from the open pipe end are recorded with a number of high speed pressure transducers. It was found that there was little to choose between the accuracy of the 1D model and the best CFD model. The CFD model did not require experimentally derived loss coefficients to accurately represent the open pipe end; however, it took several hundred times longer to complete its analysis. The best congruency between the CFD models and the experimental data was achieved using the RNG k-e turbulence model. The open end of the pipe was most effectively represented by surrounding it with a relatively small volume of cells connected to the rest of the environment using a pressure boundary.

An Experimental and Predictive Evaluation of Unsteady Gas Flow through Automotive Catalyst Elements

The incorporation of one-dimensional simulation codes within engine modelling applications has proved to be a useful tool in evaluating unsteady gas flow through elements in the exhaust system. This paper reports on an experimental and theoretical investigation into the behaviour of unsteady gas flow through catalyst substrate elements. A one-dimensional (1-D) catalyst model has been incorporated into a 1-D simulation code to predict this behaviour.

Experimental data was acquired using a ‘single pulse’ test rig. Substrate samples were tested under ambient conditions in order to investigate a range of regimes experienced by the catalyst during operation. This allowed reflection and transmission characteristics to be quantified in relation to both geometric and physical properties of substrate elements. Correlation between measured and predicted results is demonstrably good and the model provides an effective analysis tool for evaluating unsteady gas flow through different catalytic converter designs.

Simple models of complex aggregation: Vesicle formation by soft repulsive spheres with dipolelike interactions

Structural and thermodynamic properties of spherical particles carrying classical spins are investigated by Monte Carlo simulations. The potential energy is the sum of short range, purely repulsive pair contributions, and spin-spin interactions. These last are of the dipole-dipole form, with however, a crucial change of sign. At low density and high temperature the system is a homogeneous fluid of weakly interacting particles and short range spin correlations. With decreasing temperature particles condense into an equilibrium population of free floating vesicles. The comparison with the electrostatic case, giving rise to predominantly one-dimensional aggregates under similar conditions, is discussed. In both cases condensation is a continuous transformation, provided the isotropic part of the interatomic potential is purely repulsive. At low temperature the model allows us to investigate thermal and mechanical properties of membranes. At intermediate temperatures it provides a simple model to investigate equilibrium polymerization in a system giving rise to predominantly two-dimensional aggregates.

Power dissipation in nanoscale conductors

A previous tight-binding model of power dissipation in a nanoscale conductor under an applied bias is extended to take account of the local atomic topology and the local electronic structure. The method is used to calculate the power dissipated at every atom in model nanoconductor geometries: a nanoscale constriction, a one-dimensional atomic chain between two electrodes with a resonant double barrier, and an irregular nanowire with sharp corners. The local power is compared with the local current density and the local density of states. A simple relation is found between the local power and the current density in quasiballistic geometries. A large enhancement in the power at special atoms is found in cases of resonant and anti-resonant transmission. Such systems may be expected to be particularly unstable against current-induced modifications.

Space and phase resolved plasma parameters in an industrial dual-frequency capacitively coupled radio-frequency discharge

The dynamics of high energetic electrons (>= 11.7 eV) in a modified industrial confined dual-frequency capacitively coupled RF discharge (Exelan, Lam Research Inc.), operated at 1.937 MHz and 27.118 MHz, is investigated by means of phase resolved optical emission spectroscopy. Operating in a He-O-2. plasma with small rare gas admixtures the emission is measured, with one-dimensional spatial resolution along the discharge axis. Both the low and high frequency RF cycle are resolved. The diagnostic is based on time dependent measurements of the population densities of specifically chosen excited rare gas states. A time dependent model, based on rate equations, describes the dynamics of the population densities of these levels. Based on this model and the comparison of the excitation of various rare gas states, with different excitation thresholds, time and space resolved electron temperature, propagation velocity and qualitative electron density as well as electron energy distribution functions are determined. This information leads to a better understanding of the dual-frequency sheath dynamics and shows, that separate control of ion energy and electron density is limited.

A 3-D well-differentiated model of pediatric bronchial epithelium demonstrates unstimulated morphological differences between asthmatic and nonasthmatic cells

There is a need for reproducible and effective models of pediatric bronchial epithelium to study disease states such as asthma. We aimed to develop, characterize, and differentiate an effective, an efficient, and a reliable three-dimensional model of pediatric bronchial epithelium to test the hypothesis that children with asthma differ in their epithelial morphologic phenotype when compared with nonasthmatic children. Primary cell cultures from both asthmatic and nonasthmatic children were grown and differentiated at the air-liquid interface for 28 d. Tight junction formation, MUC5AC secretion, IL-8, IL-6, prostaglandin E2 production, and the percentage of goblet and ciliated cells in culture were assessed. Well-differentiated, multilayered, columnar epithelium containing both ciliated and goblet cells from asthmatic and nonasthmatic subjects were generated. All cultures demonstrated tight junction formation at the apical surface and exhibited mucus production and secretion. Asthmatic and nonasthmatic cultures secreted similar quantities of IL-8, IL-6, and prostaglandin E2. Cultures developed from asthmatic children contained considerably more goblet cells and fewer ciliated cells compared with those from nonasthmatic children. A well-differentiated model of pediatric epithelium has been developed that will be useful for more in vivo like study of the mechanisms at play during asthma.

Laser pulse compression and amplification via Raman backscattering in plasma

A simple theoretical model is proposed for the interaction between two counter-propagating laser pulses (a pump and a seed pulse) in unmagnetized plasma. Pulse compression and amplification are observed via numerical simulation. A one dimensional fluid model for stimulated Raman backscattering is proposed to investigate the pulse compression and pulse amplification mechanisms. To accomplish this, energy is transferred from the long pump pulse to a seed pulse, with a Langmuir plasma wave mediating the transfer. The study focuses on the intensity profile of the pump laser pulse. A Gaussian and a ring intensity profile are, separately, considered for the pump laser pulse.

A two compartment in vitro release model for the testing of HIV microbicide formulations

BACKGROUND: In vitro release testing of vaginal formulations is usually performed in a one-compartment model (OCM) where the release medium, usually comprising pH-adjusted water, an aqueous surfactant solution or a solvent-water solution, provides sink conditions throughout the release experiment. Although this model is useful in evaluating the effect of formulation parameters upon release, it rarely reflects in vivo conditions. Here we report use of a two-compartment diffusion cell model (TCM, comprising a small volume donor, a large volume receptor, and separated by a model epithelial membrane) to more closely mimic in vivo vaginal release and tissue absorption following administration of a UC781 vaginal ring.

METHODS: Macaque-sized matrix silicone elastomer vaginal rings containing 100mg UC781 were prepared by injection molding, and in vitro release testing performed using both OCM (20mL simulated vaginal fluid, SVF) and TCM (5mL SVF in donor cell and variable quantities of Tween 80; silicone elastomer membrane; 100mL 3:2 ethanol/water in receptor cell). In the TCM, drug levels were measured by HPLC in both donor and receptor cells, representing fluid and tissue levels respectively. Rings containing 100mg UC781 and 10% w/w Tween 80 were also manufactured and tested.

RESULTS: The amount of UC781 released from rings was significantly influenced by the choice of release model. Greatest release (56mg/14 days) was measured in the ethanol/water OCM, compared with no measurable release into SVF only. Increasing the concentration of Tween 80 in the SVF medium (1, 3 and 5% w/w) led to increased UC781 release (11, 16 and 18mg, respectively), demonstrating that vaginal fluid solubility of UC781 may be rate-determining in vivo. In the TCM, UC781 accumulates in the receptor cell medium over time, despite not being measured in the donor medium containing the ring device. Incorporation of Tween 80 directly into the ring provided enhanced release in both donor and receptor cells.

CONCLUSIONS: Release of UC781 was influenced by the choice of release medium and the inclusion of Tween 80 in the ring. Although use of SVF-only in the OCM indicated no measurable UC781 release from rings, data from the TCM confirms that UC781 is not only released but is also capable of penetrating across the model epithelial membrane. The TCM may therefore provide a more representative in vitro release model for mimicking in vivo absorption.

A numerical investigation of the flow structures and losses for turbulent flow in 90 degrees elbow bends

Computational fluid dynamic modelling was carried out on a series of pipe bends having R/r values of 1.3, 5, and 20, with the purpose of determining the accuracy of numerical models in predicting pressure loss data from which to inform one-dimensional loss models. Four separate turbulence models were studied: the standard k-epsilon model, realizable k-epsilon model, k-omega model, and a Reynolds stress model (RSM). The results are presented for each bend in the form of upstream and downstream pressure profiles, pressure distributions along the inner and outer walls, detailed pressure and velocity fields as well as overall loss values. In each case, measured data were presented to evaluate the predictive ability of each model. The RSM was found to perform the best, producing accurate pressure loss data for bends with R/r values of 5 and 20. For the tightest bend with an R/r value of 1.3, however, predictions were significantly worse due to the presence of flow separation, stronger pressure gradients, and high streamline curvature.

Food-web structure in low- and high-dimensional trophic niche spaces

A question central to modelling and, ultimately, managing food webs concerns the dimensionality of trophic niche space, that is, the number of independent traits relevant for determining consumer-resource links. Food-web topologies can often be interpreted by assuming resource traits to be specified by points along a line and each consumer's diet to be given by resources contained in an interval on this line. This phenomenon, called intervality, has been known for 30 years and is widely acknowledged to indicate that trophic niche space is close to one-dimensional. We show that the degrees of intervality observed in nature can be reproduced in arbitrary-dimensional trophic niche spaces, provided that the processes of evolutionary diversification and adaptation are taken into account. Contrary to expectations, intervality is least pronounced at intermediate dimensions and steadily improves towards lower- and higher-dimensional trophic niche spaces.

Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

Collective decoherence of cold atoms coupled to a Bose-Einstein condensate

We examine the time evolution of cold atoms (impurities) interacting with an environment consisting of a degenerate bosonic quantum gas. The impurity atoms differ from the environment atoms, being of a different species. This allows one to superimpose two independent trapping potentials, each being effective only on one atomic kind, while transparent to the other. When the environment is homogeneous and the impurities are confined in a potential consisting of a set of double wells, the system can be described in terms of an effective spin-boson model, where the occupation of the left or right well of each site represents the two (pseudo)-spin states. The irreversible dynamics of such system is here studied exactly, i.e. not in terms of a Markovian master equation. The dynamics of one and two impurities is remarkably different in respect of the standard decoherence of the spin-boson system. In particular, we show: (i) the appearance of coherence oscillations, (ii) the presence of super and subdecoherent states that differ from the standard ones of the spin-boson model, and (iii) the persistence of coherence in the system at long times. We show that this behaviour is due to the fact that the pseudospins have an internal spatial structure. We argue that collective decoherence also prompts information about the correlation length of the environment. In a one-dimensional (1D) configuration, one can change even more strongly the qualitative behaviour of the dephasing just by tuning the interaction of the bath.

Ground state of low-dimensional dipolar gases: Linear and zigzag chains

We study the ground-state phase diagram of ultracold dipolar gases in highly anisotropic traps. Starting from a one-dimensional geometry, by ramping down the transverse confinement along one direction, the gas reaches various planar distributions of dipoles. At large linear densities, when the dipolar gas exhibits a crystal-like phase, critical values of the transverse frequency exist below which the configuration exhibits transverse patterns. These critical values are found by means of a classical theory, and are in full agreement with classical Monte Carlo simulations. The study of the quantum system is performed numerically with Monte Carlo techniques and shows that the quantum fluctuations smoothen the transition and make it completely disappear in a gas phase. These predictions could be experimentally tested and would allow one to reveal the effect of zero-point motion on self-organized mesoscopic structures of matter waves, such as the transverse pattern of the zigzag chain.

In Silico Prediction of the Mechanobiological Response of Arterial Tissue: Application to Angioplasty and Stenting

One way to restore physiological blood flow to occluded arteries involves the deformation of plaque using an intravascular balloon and preventing elastic recoil using a stent. Angioplasty and stent implantation cause unphysiological loading of the arterial tissue, which may lead to tissue in-growth and reblockage; termed “restenosis.” In this paper, a computational methodology for predicting the time-course of restenosis is presented. Stress-induced damage, computed using a remaining life approach, stimulates inflammation (production of matrix degrading factors and growth stimuli). This, in turn, induces a change in smooth muscle cell phenotype from contractile (as exists in the quiescent tissue) to synthetic (as exists in the growing tissue). In this paper, smooth muscle cell activity (migration, proliferation, and differentiation) is simulated in a lattice using a stochastic approach to model individual cell activity. The inflammation equations are examined under simplified loading cases. The mechanobiological parameters of the model were estimated by calibrating the model response to the results of a balloon angioplasty study in humans. The simulation method was then used to simulate restenosis in a two dimensional model of a stented artery. Cell activity predictions were similar to those observed during neointimal hyperplasia, culminating in the growth of restenosis. Similar to experiment, the amount of neointima produced increased with the degree of expansion of the stent, and this relationship was found to be highly dependant on the prescribed inflammatory response. It was found that the duration of inflammation affected the amount of restenosis produced, and that this effect was most pronounced with large stent expansions. In conclusion, the paper shows that the arterial tissue response to mechanical stimulation can be predicted using a stochastic cell modeling approach, and that the simulation captures features of restenosis development observed with real stents. The modeling approach is proposed for application in three dimensional models of cardiovascular stenting procedures.