Modelling soot formation in a DISI engine

In this work, the formation of soot in a Direct Injection Spark Ignition (DISI) engine is simulated using the Stochastic Reactor Model (SRM) engine code. Volume change, convective heat transfer, turbulent mixing, direct injection and flame propagation are accounted for. In order to simulate flame propagation, the cylinder is divided into an unburned, entrained and burned zone, with the rate of entrainment being governed by empirical equations but combustion modelled with chemical kinetics. The model contains a detailed chemical mechanism as well as a highly detailed soot formation model, however computation times are relatively short. The soot model provides information on the morphology and chemical composition of soot aggregates along with bulk quantities, including soot mass, number density, volume fraction and surface area. The model is first calibrated by simulating experimental data from a Gasoline Direct Injection (GDI) Spark Ignition (SI) engine. The model is then used to simulate experimental data from the literature, where the numbers, sizes and derived mass particulate emissions from a 1.83 L, 4-cylinder, 4 valve production DISI engine were examined. Experimental results from different injection and spark timings are compared with the model and the qualitative trends in aggregate size distribution and emissions match the exhaust gas measurements well. © 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Time resolved measurement of cold start hc concentration using the fast FID

Understanding mixture formation phenomena during the first few cycles of an engine cold start is extremely important for achieving the minimum engine-out emission levels at the time when the catalytic converter is not yet operational. Of special importance is the structure of the charge (film, droplets and vapour) which enters the cylinder during this time interval as well as its concentration profile. However, direct experimental studies of the fuel behaviour in the inlet port have so far been less than fully successful due to the brevity of the process and lack of a suitable experimental technique. We present measurements of the hydrocarbon (HC) concentration in the manifold and port of a production SI engine using the Fast Response Flame Ionisation Detector (FRFID). It has been widely reported in the past few years how the FRFID can be used to study the exhaust and in-cylinder HC concentrations with a time resolution of a few degrees of crank angle, and the device has contributed significantly to the understanding of unburned HC emissions. Using the FRFID in the inlet manifold is difficult because of the presence of liquid droplets, and the low and fluctuating pressure levels, which leads to significant changes in the response time of the instrument. However, using recently developed procedures to correct for the errors caused by these effects, the concentration at the sampling point can be reconstructed to align the FRFID signal with actual events in the engine. © 1996 Society of Automotive Engineers, Inc.

Investigation of the dilution process for measurement of particulate matter from spark-ignition engines

Measurements of particulate matter (PM) from spark ignition (SI) engine exhaust using dilution tunnels will become more prevalent as emission standards are tightened. Hence, a study of the dilution process was undertaken in order to understand how various dilution related parameters affect the accuracy with which PM sizes and concentrations can be determined. A SI and a compression ignition (CI) engine were separately used to examine parameters of the dilution process; the present work discusses the results in the context of SI exhaust dilution. A Scanning Mobility Particle Sizer (SMPS) was used to measure the size distribution, number density, and volume fraction of PM. Temperature measurements in the exhaust pipe and dilution tunnel reveal the degree of mixing between exhaust and dilution air, the effect of flowrate on heat transfer from undiluted and diluted exhaust to the environment, and the minimum permissible dilution ratio for a maximum sample temperature of 52°C. Measurements of PM concentrations as a function of dilution ratio show the competing effects of temperature and particle/vapor concentrations on particle growth dynamics, which result in a range of dilution ratios-from 13 to 18-where the effect of dilution ratio, independent of flowrate, is kept to a minimum. This range of dilution ratios is therefore optimal in order to achieve repeatable PM concentration measurements. Particle dynamics during transit through the tunnel operating at the optimal dilution ratio was found statistically insignificant compared to data scatter. Such small differences in number concentration may be qualitatively representative of particle losses for SI exhaust, but small increases in PM volume fraction during transit through the tunnel may significantly underestimate accretion of mass due to unburned hydrocarbons (HCs) emitted by SI engines. The fraction of SI-derived PM mass due to adsorbed/absorbed vapor, estimated from these data, is consistent with previous chemical analyses of PM. © 1998 Society of Automotive Engineers, Inc.

Development of a time and space resolved sampling probe diagnostic for engine exhaust hydrocarbons

In order to understand how unburned hydrocarbons emerge from SI engines and, in particular, how non-fuel hydrocarbons are formed and oxidized, a new gas sampling technique has been developed. A sampling unit, based on a combination of techniques used in the Fast Flame Ionization Detector (FFID) and wall-mounted sampling valves, was designed and built to capture a sample of exhaust gas during a specific period of the exhaust process and from a specific location within the exhaust port. The sampling unit consists of a transfer tube with one end in the exhaust port and the other connected to a three-way valve that leads, on one side, to a FFID and, on the other, to a vacuum chamber with a high-speed solenoid valve. Exhaust gas, drawn by the pressure drop into the vacuum chamber, impinges on the face of the solenoid valve and flows radially outward. Once per cycle during a specified crank angle interval, the solenoid valve opens and traps exhaust gas in a storage unit, from which gas chromatography (GC) measurements are made. The port end of the transfer tube can be moved to different locations longitudinally or radially, thus allowing spatial resolution and capturing any concentration differences between port walls and the center of the flow stream. Further, the solenoid valve's opening and closing times can be adjusted to allow sampling over a window as small as 0.6 ms during any portion of the cycle, allowing resolution of a crank angle interval as small as 15°CA. Cycle averaged total HC concentration measured by the FFID and that measured by the sampling unit are in good agreement, while the sampling unit goes one step further than the FFID by providing species concentrations. Comparison with previous measurements using wall-mounted sampling valves suggests that this sampling unit is fully capable of providing species concentration information as a function of air/fuel ratio, load, and engine speed at specific crank angles. © Copyright 1996 Society of Automotive Engineers, Inc.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.