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.

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.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

Absolute Atomic Oxygen Density Distributions in the Effluent of a Micro Scale Atmospheric Pressure Plasma Jet

The coplanar microscale atmospheric pressure plasma jet (µ-APPJ) is a capacitively coupled radio frequency discharge (13.56 MHz, ~15W rf power) designed for optimized optical diagnostic access. It is operated in a homogeneous glow mode with a noble gas flow (1.4 slm He) containing a small admixture of molecular oxygen (~0.5%). Ground state atomic oxygen densities in the effluent up to 2 × 1014 cm-3 are measured by two-photon absorption laser-induced fluorescence spectroscopy (TALIF) providing space resolved density maps. The quantitative calibration of the TALIF setup is performed by comparative measurements with xenon. A maximum of the atomic oxygen density is observed for 0.6% molecular oxygen admixture. Furthermore, an increase in the rf power up to about 15W (depending on gas flow and mixture) leads to an increase in the effluent’s atomic oxygen density, then reaching a constant level for higher powers.

A survey of molecular line emission towards ultracompact HII regions

We have used the JCMT to survey molecular line emission towards 14 ultracompact HII regions (G5.89, G9.62, G10.30, G10.47, G12.21, G13.87, G29.96, G31.41, G34.26, G43.89, G45.12, G45.45, G45.47, and G75.78). For each source, we observed up to ten 1 GHz bands between 200 and 350 GHz, covering lines of more than 30 species including multiple transitions of CO isotopes, CH3OH, CH3CCH, CH3CN and HCOOCH3, and sulphuretted molecules. The number of transitions detected varied by a factor of 20 between sources; which were chosen following observations of high-excitation ammonia (Cesaroni et al. 1994a) and methyl cyanide (Olmi et al. 1993). In half our sample (the line-poor sources), only (CO)-O-17: (CO)-O-18, SO, (CS)-S-34 and CH3OH were detected. In the line-rich sources, we detected over 150 lines, including high excitation lines of CH3CN, HCOOCH3; C2H5CN, CH3OH, and CH3CCH. We have calculated the physical conditions of the molecular gas. To reproduce the emission from the line-rich sources requires both a hot, dense compact core and an ambient cloud consisting of less dense, cooler gas. The hot cores, which are less than 0.1 pc in size; reach densities of at least 10(8) cm(-3) and temperatures of more than 80 K. The line-poor sources can be modelled without a hot core by a 20-30 K, 10(5) cm(-3) cloud. We find no correlation between the size of the HII region and the current physical conditions in the molecular environment. A comparison with chemical models (Millar et al. 1997) confirms that grain surface chemistry is important in hot cores.

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Diagnostics on an atmospheric pressure plasma jet

The atmospheric pressure plasma jet (APPJ) is a homogeneous non-equilibrium discharge at ambient pressure. It operates with a noble base gas and a percentage-volume admixture of a molecular gas. Applications of the discharge are mainly based on reactive species in the effluent. The effluent region of a discharge operated in helium with an oxygen admixture has been investigated. The optical emission from atomic oxygen decreases with distance from the discharge but can still be observed several centimetres in the effluent. Ground state atomic oxygen, measured using absolutely calibrated two-photon laser induced fluorescence spectroscopy, shows a similar behaviour. Detailed understanding of energy transport mechanisms requires investigations of the discharge volume and the effluent region. An atmospheric pressure plasma jet has been designed providing excellent diagnostics access and a simple geometry ideally suited for modelling and simulation. Laser spectroscopy and optical emission spectroscopy can be applied in the discharge volume and the effluent region.

Absolute atomic oxygen density profiles in the discharge core of a microscale atmospheric pressure plasma jet

The micro atmospheric pressure plasma jet is an rf driven (13.56 MHz, ~20 W) capacitively coupled discharge producing a homogeneous plasma at ambient pressure when fed with a gas flow of helium (1.4 slm) containing small admixtures of oxygen (~0.5%). The design provides excellent optical access to the plasma core. Ground state atomic oxygen densities up to 3x1016 cm-3 are measured spatially resolved in the discharge core by absolutely calibrated two-photon absorption laser-induced fluorescence spectroscopy. The atomic oxygen density builds up over the first 8 mm of the discharge channel before saturating at a maximum level. The absolute value increases linearly with applied power.

Sulphur chemistry and evolution in hot cores

We compare the results of our JCMT spectral line survey of molecular gas towards ultracompact HII regions with the predictions of models of sulphur chemistry in hot cores. We investigate the range of evolutionary models that are consistent with the observed physical conditions and chemical abundances, and see to what extent it is possible to constrain core ages by comparing abundances with the predictions of chemical models. The observed abundance ratios vary little from source to source, suggesting that all the sources are at a similar evolutionary stage. The models are capable of predicting the observed abundances of H2S, SO, SO2, and CS. The models fail to predict the amount of OCS observed, suggesting that an alternative formation route is required. An initial H2S abundance from grain mantle evaporation of similar to 10(-7) is preferred.

Polycrystalline silicon film growth in a SiF4/SiH4/H2 plasma

The growth of polycrystalline silicon (polysilicon) films from SiF4/SiH4/H2 gas mixtures is reported. The polysilicon films have been deposited in a multi process reactor by a PECVD process. The effect of r.f. power, chamber temperature and gas flow ratios on grain size and deposition rate have been determined. The fluorine concentration and the grain sizes of the films have been determined by SIMS and atomic force microscopy (AFM), respectively. Grain sizes in excess of 900 A are reported for layers deposited at 300°C. © 1999 Published by Elsevier Science S.A. All rights reserved.

Parameter studies of fibre laser micro-welding of AISI 316L using Taguchi method

This paper discusses the application of the Taguchi experimental design approach in optimizing the key process parameters for micro-welding of thin AISI 316L foil using the 100W CW fibre laser. A L16 Taguchi experiment was conducted to systematically understand how the power, scanning velocity, focus position, gas flow rate and type of shielding gas affect the bead dimensions. The welds produced in the L16 Taguchi experiment was mainly of austenite cellular-dendrite structure with an average grain size of 5µm. An exact penetration weld with the largest penetration to fusion width ratio was obtained. Among those process parameters, the interaction between power and scanning velocity presented the strongest effect to the penetration to fusion width ratio and the power was found to be the predominantly important factor that drives the interaction with other factors to appreciably affect the bead dimensions.

Understanding the synthesis of directly spinnable carbon nanotube forests

Carbon nanotube forests that can be spun directly from the growth substrate into pure, highly aligned webs, ribbons or yarn promise novel applications that capture the strength and other characteristics of this material. The precise conditions for high spinnability over a maximum proportion of the reactor space are extremely sensitive. The roles of catalyst, substrate, temperature, gas flow rates, reaction time with acetylene etc. were studied to identify and understand the key parameters and develop a robust, scalable process. Using a 44 mm (id) reactor, the optimum values for these variables were determined as comprising a 2.3 nm thick iron catalyst layer on a silicon substrate with 50 nm of thermal oxide; 670 °C running temperature; 650 sccm helium and 34 sccm acetylene for 20 min. The effects of deviating from these optima were explored and the role of amorphous carbon deposition clarified. Crown Copyright © 2009.

Multiphysics Simulation of Thermal Plasma Produced by Lightning Strike on Aircraft Structures

A novel numerical technique is proposed to model thermal plasma of microseconds/milliseconds time-scale effect. Modelling thermal plasma due to lightning strike will allow the estimation of electric current density, plasma pressure, and heat flux at the surface of the aircraft structure. These input data can then be used for better estimation of the mechanical/thermal induced damage on the aircraft structures for better protection systems design. Thermal plasma generated during laser cutting, electric (laser) welding and other plasma processing techniques have been the focus of many researchers. Thermal plasma is a gaseous state that consists from a mixture of electrons, ions, and natural particles. Thermal plasma can be assumed to be in local thermodynamic equilibrium, which means the electrons and the heavy species have equal temperature. Different numerical techniques have been developed using a coupled Navier Stokes – Heat transfer – Electromagnetic equations based on the assumption that the thermal plasma is a single laminar gas flow. These previous efforts focused on generating thermal plasma of time-scale in the range of seconds. Lighting strike on aircraft structures generates thermal plasma of time-scale of milliseconds/microseconds, which makes the previous physics used not applicable. The difficulty comes from the Navier-Stokes equations as the fluid is simulated under shock load, this introducing significant changes in the density and temperature of the fluid.