Search Results

A zero-dimensional code is developed to simulate turbulent spray combustion and NOx and soot emission in direct injection diesel engines. The code consists of two major parts; mixing calculation for the probability density function (PDF) based on the multi-zone model by Hiroyasu et al., (1983) and the flame structure by the conditional moment closure (CMC) model (Klimenko & Bilger, 1999). The skeletal mechanism of n-heptane is employed with the elementary reaction steps for heat release and the NOx chemistry in GRI 3.0. The spray model accounts for evaporation and mixing based on momentum balance of the spray zones, while the CMC model incorporates the conditional flame structures with one fuel group or flame structure for each injection. The spatially integrated density-weighted PDF, F(η), is defined to represent inhomogeneous mixture distribution in the cylinder. The one-equation soot model is employed for prediction of the soot emission.

A spark ignition engine is simulated for its flow field and combustion using a modified version of KIVA-II. The code handles an arbitrary 3-D computational domain and boundary conditions in the Cartesian coordinate. The valve is fixed at a half open position and the pressure boundary condition is specified at the port upstream. The grid is generated by a 2-D CAD package and later combined to a 3-D. Post-processing is performed by a finite element postprocessor package. Results for a typical spark ignition engine show reasonable trends for the velocity field and flame propagation throughout the intake, compression, combustion and expansion stroke.

Control of turbulence during the compression stroke is suggested by both theoretical calculations and experimental results obtained with an LDV measurement in a motored engine. The authors have found experimentally that when an axial distribution of swirl intensity exists, a large-scale annular vortex is formed inside the cylinder during the compression stroke and this vortex generates and transports turbulence energy. A numerical calculation is adopted to elucidate this phenomenon. Then, an axial stratification of swirl intensity is found to generate a large-scale annular vortex during the compression stroke by an interaction between the piston motion and the axial pressure gradient. The initial swirl profile is parametrically varied to assess its effect on the turbulence parameters. Among calculated results, turbulence energy is enhanced strongest when the swirl intensity is highest at the piston top surface and lowest at the bottom surface of the cylinder head.

According to previous papers on the combustion process in LHR diesel engines the combustion seems to deteriorate in LHR diesel engines. However it has been unclear whether this was caused by the high temperature gas or high temperature combustion chamber walls. This study was intended to investigate the effect of gas temperature on the rate of heat release through the heat release analysis and other measurements using a rapid compression-expansion machine. Experiments conducted at high gas temperatures which was achieved by the employment of oxygen-argon-helium mixture made it clear that the combustion at a high gas temperature condition deteriorated actually and this was probably due to the poorer mixing rate because of the increase in gas viscosity at a high gas temperature condition.

The applicability of several popular diesel particulate matter (PM) measurement techniques to low temperature combustion is examined. The instruments' performance in measuring low levels of PM from advanced diesel combustion is evaluated. Preliminary emissions optimization of a high-speed light-duty diesel engine was performed for two conventional and two advanced low temperature combustion engine cases. A low PM (<0.2 g/kg_fuel) and NOx (<0.07 g/kg_fuel) advanced low temperature combustion (LTC) condition with high levels of exhaust gas recirculation (EGR) and early injection timing was chosen as a baseline. The three other cases were selected by varying engine load, injection timing, injection pressure, and EGR mass fraction. All engine conditions were run with ultra-low sulfur diesel fuel. An extensive characterization of PM from these engine operating conditions is presented.

In this study a new model is proposed for turbulent premixed combustion in a spark-ignition engine. An independent transport equation is solved for the mean reaction progress variable in a propagation form in KIVA-3V. An expression for turbulent burning velocity was previously given as a product of turbulent diffusivity in unburned gas, laminar flame speed and maximum flame surface density. The model has similarity with the G equation approach, but originates from zone conditionally averaged formulation for unburned gas. A spark kernel grows initially as a laminar flame and becomes a fully developed turbulent flame brush according to a transition criterion in terms of the kernel size and the integral length scale. Simulation of a homogeneous charge pancake chamber engine showed good agreement with measured flame propagation and pressure trace. The model was also applied against experimental data of Hyundai θ-2.0L SI engine.

This study presents an extended Conditional Moment Closure (CMC) model [1] for turbulent spray combustion of a diesel engine. A spatially integrated CMC equation involves multiple flame structures, evaporation source terms for mixture fraction variance and scaled conditional scalar dissipation rates (CSDR) and probability density functions (PDF). An independent transport equation is solved for each flame group with equal mass of sequentially injected and evaporated fuel group. The implementation strategy in KIVA is to determine mean scalar variables instead of mean reaction rates in terms of conditional flame structures and local PDF's. Simulation is performed to test the suggested CMC model for a heavy duty diesel engine with early injection timings. NOx chemistry [2] is combined with skeletal chemistry of n-heptane [3] while soot is modeled by both one equation model [4] and two equation model [5].

Combustion characteristics of the emulsified diesel fuels are investigated in a Rapid Compression and Expansion Machine(RCEM). Among the test cases the 40 W/O fuel injected at BTDC 20° has shown the best performance with respect to the efficiency and NOx and soot emissions. The pressure trace of the 40 W/O fuel is characterized by a longer ignition delay and a lower rate of pressure rise in premixed combustion. High speed photographs show reduced flame luminosity and lower flame temperature with the increasing W/O ratio. Micro-explosions of emulsified fuel droplets which affect the local shape and brightness of the flame are identified in magnified flame images.

The accuracy of the injection rate meter based on W. Zeuch's method in the measurement of multiple injection rate and amount was calibrated using a small cam driven piston that is driven by an electric motor. For the pre- or early-injection, a sensor with a high sensitivity can be applied to measure the small pressure increase due to the small injection amount. In case of the multiple injection that has the post and/or late injection, a pressure sensor with a low sensitivity must cover not only the large pressure increase due to the main injection but also the small pressure increase due to the post and/or late injection because the output of the high sensitivity sensor is saturated after the main injection. So the linearity of the low sensitivity pressure sensor was calibrated with the cam driven piston prior to the experiment with the actual injection system.

Analysis of flow field and charge distribution in a gasoline direct-injection (GDI) engine is performed by a modified version of the KIVA code. A particle-based spray model is proposed to simulate a swirl-type hollow-cone spray in a GDI engine. Spray droplets are assumed to be fully atomized and introduced at the sheet breakup locations as determined by experimental correlations and energy conservation. The effects of the fuel injection parameters such as spray cone angle and ambient pressure are examined for different injectors and injection conditions. Results show reasonable agreement with the measurements for penetration, dispersion, global shape, droplet velocity and size distribution by Phase Doppler Particle Anemometry(PDPA) in a constant-volume chamber. The test engine is a 4-stroke 4-valve optically accessible single-cylinder engine with a pent-roof head and tumble ports.

Submodels are developed for injection, evaporation and wall impingement of a liquid LPG spray. The injection model determines the quality of fuel as two-phase choke flow at the nozzle exit. Wind tunnel experiments show the spray penetration more sensitive to ambient flow velocity than to injection pressure. Most evaporation occurs during choking, while heat transfer from surrounding air has a negligible effect on downstream droplet sizes. Three dimensional simulation shows that the bathtub cavity is better than the dog-dish cavity for stable flame propagation in lean-burn conditions. The injection timing during the IVC period has a negligible effect, while injection during an intake stroke enhances fuel/air mixing to result in more homogeneous cylinder charge.

The relation between the characteristics of a non-evaporating spray and those of a corresponding frame achieved in a rapid compression machine was investigated experimentally. The fuel injection pressure was changed in a range of 55 to 260 MPa and the other injection parameters such as orifice diameter and injection duration were changed systematically. The characteristics of the non-evaporating spray such as the Sauter mean diameter and the mean excess air ratio of the spray were measured by an image analysis technique. The time required for a pressure rise due to combustion was taken as an index to characterize the flame. It was concluded that the mean excess air ratio of a spray is the major factor which controls the burning rate and that the high injection pressure is effective in shortening the combustion duration and reducing soot formation.

The two-dimensional distribution of a soot cloud in an unsteady spray flame in a rapid compression machine(RCM) was visualized using the laser sheet scattering technique. A 40 mm x 50 mm cross section on the flame axis was illuminated by a thin laser sheet from a single pulsed Nd:YAG laser(wavelength 532 nm). Scattered light from soot particles was taken by a CCD camera via a high speed gated image intensifier. The temporal variation of the scattered light images were presented with the injection pressure as a parameter. The results showed that scattered light was intense near the periphery of the flame tip and that the scattered light becomes weaker significantly and disappears fast after the end of injection as injection pressure is increased. This technique was also applied to the visualization of the two-dimensional distribution of liquid droplets in the non-evaporating spray to correlate it with the soot concentration distribution.

Multidimensional cycle simulation is performed for a spark ignition engine for the entire four stroke cycle of intake, compression, expansion and exhaust stroke. Combustion is divided into ignition period and flame propagation period. Simple Arrhenius form is used in ignition and Bray-Moss model is used for turbulence controlled flame propagation. Heat flux is calculated by the wall function method in k-ε turbulence model. Grid is generated by a two dimensional CAD package. Results show proper trends for all flow variables, burning rate and heat flux.

A low temperature combustion (LTC) diesel engine has been under investigation for reduction of NOx and soot with acceptable compromise in the efficiency through modification of the combustion process. In this paper computational simulation is performed as a preliminary step for development of an LTC diesel engine for off-highway construction vehicles. Validation is performed for major physical models against measurements in LTC conditions. The conditional moment closure (CMC) is employed to address coupling between chemistry and turbulence in KIVA-CMC. The Kelvin-Helmholtz/Rayleigh-Taylor (KH-RT) model is employed for spray breakup and a skeletal n-heptane mechanism for both low and high temperature chemistry. Parametric evaluation is performed for design and operating conditions including EGR rate and injection timing. Results are obtained for efficiency, IMEP, CO, NOx and PM emissions at intake boost pressures of 1, 2 and 3 bar.

Computation is performed to predict number density, volume fraction and size distribution of soot particles in typical operating conditions of a diesel engine. KIVA has been integrated with the CMC routine to consider turbulence/chemistry coupling and gas phase kinetics for heat release and soot precursors. The compositions of soot precursors are estimated by tracking Lagrangian particles to consider spatial inhomogeneity and differential diffusion in KIVA. The soot simulator SWEEP is employed as a postprocessing step to calculate conditional and integral quantities of soot particles. There are larger particles produced at a higher load or a lower rpm, but no consistent trend for injection timing in the conditional size distribution at the mixture fraction of 0.12. The integral results are obtained for number density, total mass and size distribution by summing up the histories of all tracked particles in the cylinder.

The conditional moment closure (CMC) model has successfully been employed for simulation of turbulent nonpremixed combustion in diverse engineering problems. The solution procedure was developed for Lagrangian CMC with a single or multiple fuel or flame groups for a diesel engine in our previous work. Multiple flame groups may consider different histories and residence times through the sequence of evaporation, ignition and mixing controlled combustion. The Lagrangian CMC model has an advantage of less computational burden than Eulerian CMC with fixed spatially resolved grids. However, it cannot allow interaction among different flame groups through premixed combustion, while such interaction occurs through spatial transport terms in Eulerian CMC. In this work an improved version of the Lagrangian CMC model is proposed to consider interaction by premixed flame propagation along a constant mixture fraction contour between neighboring flame groups.

A thermodynamic two-zone model which assumes a stoichiornetric burned gas region and unburned air region is presented in an attempt to calculate more precise rate of heat release of diesel combustion. A comparison is made of the rate of heat release obtained by the two-zone model with that obtained by the conventional single-zone model. It shows around 10 % increase in the rate of heat release with the two-zone model. The effect of state equation of gas is also examined with the single-zone model and the use of a real gas law in stead of the perfect gas law is found to yield minor difference in the rate of heat release at a high boost operating condition.

A modified version of KIVA-II code is applied to scavenging flow simulation of a four-poppet-valved two-stroke engine. The standard k-ε turbulence model is used with no slip on the wall. Grid is generated through direct interface with three dimensional CAD data using a commercial CAE package. Valve motion is taken into account by identifying the cells occupied by the valve structure and putting them as solid obstacle cells. Results show reasonable trends for variation of the velocity field and fresh air mass fraction distribution with crank angle. Parametric study shows that the intake port orientation does not have much effect on the cylinder flow and scavenging efficiency due to strong flow diversion by the valves.

To understand further the mixing process between the injected fuel and air in the combustion chamber of a diesel engine, the turbulent mixing process in a one-phase, two-dimensional transient jet was theoretically studied using the discrete vortex simulation. First, the simulation model was evaluated by comparisons between calculated and experimental data on two-dimensional turbulent jets. Second, the trajectories of the injected fluid elements marked with different colors were graphically demonstrated. Also the process of entrainment of the surrounding fluid into the jet was visually presented using colored tracers.