Search Results

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

Due to the new challenge of meeting number-based regulations for particulate matter (PM), a numerical and experimental study has been conducted to better understand particulate formation in engines fuelled with compressed natural gas. The study has been conducted on a Heavy-Duty, Euro VI, 4-cylinder, spark ignited engine, with multipoint sequential phased injection and stoichiometric combustion. For the experimental measurements two different instruments were used: a condensation particle counter (CPC) and a fast-response particle size spectrometer (DMS) the latter able also to provide a particle size distribution of the measured particles in the range from 5 to 1000 nm. Experimental measurements in both stationary and transient conditions were carried out. The data using the World Harmonized Transient Cycle (WHTC) were useful to detect which operating conditions lead to high numbers of particles. Then a further transient test was used for a more detailed and deeper analysis.

Reactivity Controlled Compression Ignition (RCCI) is an engine combustion strategy that utilizes in-cylinder fuel blending to produce low NOx and PM emissions while maintaining high thermal efficiency. The current study investigates RCCI and conventional diesel combustion (CDC) operation in a light-duty multi-cylinder engine over transient operating conditions using a high-bandwidth, transient capable engine test cell. Transient RCCI and CDC combustion and emissions results are compared over an up-speed change from 1,000 to 2,000 rev/min. and a down-speed change from 2,000 to 1,000 rev/min. at a constant 2.0 bar BMEP load. The engine experiments consisted of in-cylinder fuel blending with port fuel-injection (PFI) of gasoline and early-cycle, direct-injection (DI) of ultra-low sulfur diesel (ULSD) for the RCCI tests and the same ULSD for the CDC tests.

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

3-D Computational Fluid Dynamics (CFD) simulations have been performed to study particulate formation in a Spark-Ignition (SI) engine under premixed conditions. A semi-detailed soot model and a chemical kinetic model, including poly-aromatic hydrocarbon (PAH) formation, were coupled with a spark ignition model and the G equation flame propagation model for SI engine simulations and for predictions of soot mass and particulate number density. The simulation results for in-cylinder pressure and particle size distribution (PSDs) are compared to available experimental studies of equivalence ratio effects during premixed operation. Good predictions are observed with regard to cylinder pressure, combustion phasing and engine load. Qualitative agreements of in-cylinder particle distributions were also obtained and the results are helpful to understand particulate formation processes.

Reactivity Controlled Compression Ignition (RCCI) is an engine combustion strategy that utilizes in-cylinder fuel blending to produce low NOx and PM emissions, while maintaining high thermal efficiency. Previous RCCI steady-state performance studies provided a fundamental understanding of the RCCI combustion process in steady-state, single-cylinder and multi-cylinder engine tests. The current study investigates RCCI and conventional diesel combustion (CDC) operation in a light-duty multi-cylinder engine over transient operating conditions. In this study, a high-bandwidth, transient-capable engine test cell was used and multi-cylinder engine RCCI combustion is compared to CDC over a step load change from 1 to 4 bar BMEP at 1,500 rev/min. The engine experiments consisted of in-cylinder fuel blending using port fuel-injection (PFI) of gasoline and early-cycle, direct-injection (DI) of ultra-low sulfur diesel (ULSD) for the RCCI tests and used the same ULSD for the CDC tests.

Advanced combustion systems that simultaneously address PM and NOx while retaining the high efficiency of modern diesel engines, are being developed around the globe. One of the most difficult problems in the area of advanced combustion technology development is the control of combustion initiation and retaining power density. During the past several years, significant progress has been accomplished in reducing emissions of NOx and PM through strategies such as LTC/HCCI/PCCI/PPCI and other advanced combustion processes; however control of ignition and improving power density has suffered to some degree - advanced combustion engines tend to be limited to the 10 bar BMEP range and under. Experimental investigations have been carried out on a light-duty DI multi-cylinder diesel automotive engine. The engine is operated in low temperature combustion (LTC) mode using 93 RON (Research Octane Number) and 74 RON fuel.

A generalized re-normalization group (RNG) turbulence model based on the local "dimensionality" of the flow field is proposed. In this modeling approach the model coefficients C₁, C₂, and C₃ are all constructed as functions of flow strain rate. In order to further validate the proposed turbulence model, the generalized RNG closure model was applied to model the backward facing step flow (a classic test case for turbulence models). The results indicated that the modeling of C₂ in the generalized RNG closure model is reasonable, and furthermore, the predictions of the generalized RNG model were in better agreement with experimental data than the standard RNG turbulence model. As a second step, the performance of the generalized RNG closure was investigated for a complex engine flow.

The present study attempted to model experimental results obtained on an optical engine at the Sandia National Laboratory. Measurements of in-cylinder unburned hydrocarbon (UHC) distributions were provided using advanced optical diagnostics on a near production type piston. Previous multidimensional modeling provided accurate pressure profiles and heat release rate (HRR) predictions. However, the experimental UHC distribution was not matched, and the model predicted UHC extending from the bowl into the squish region in the expansion stroke. To explore the causes of this discrepancy a parametric study was performed using a variety of initial conditions, boundary conditions and model constants to explore their effects on the UHC distribution. Of the initial conditions, the swirl ratio was found to have the biggest impact on the UHC distribution.

RNG k-ε closure turbulence dissipation equations are evaluated employing the CFD code KIVA-3V Release 2. The numerical evaluations start by considering simple jet flows, including incompressible air jets and compressible helium jets. The results show that the RNG closure turbulence model predicts lower jet tip penetration than the "standard" k-ε model, as well as being lower than experimental data. The reason is found to be that the turbulence kinetic energy is dissipated too slowly in the downstream region near the jet nozzle exit. In this case, the over-predicted R term in RNG model becomes a sink of dissipation in the ε-equation. As a second step, the RNG turbulence closure dissipation models are further tested in complex engine flows to compare against the measured evolution of turbulence kinetic energy, and an estimate of its dissipation rate, during both the compression and expansion processes.

Sources of unburned hydrocarbon (UHC) emissions are examined for a highly dilute (10% oxygen concentration), moderately boosted (1.5 bar), low load (3.0 bar IMEP) operating condition in a single-cylinder, light-duty, optically accessible diesel engine undergoing partially-premixed low-temperature combustion (LTC). The evolution of the in-cylinder spatial distribution of UHC is observed throughout the combustion event through measurement of liquid fuel distributions via elastic light scattering, vapor and liquid fuel distributions via laser-induced fluorescence, and velocity fields via particle image velocimetry (PIV). The measurements are complemented by and contrasted with the predictions of multi-dimensional simulations employing a realistic, though reduced, chemical mechanism to describe the combustion process.

Improvements to combustion models for modeling spark ignition engines using the G-equation flame propagation model and detailed chemical kinetics have been performed. The improvements include revision of a PRF chemistry mechanism, precise calculation of “primary heat release” based on the sub-grid scale unburned/burnt volumes of flame-containing cells, modeling flame front quenching in highly stratified mixtures, introduction of a Damkohler model for assessing the combustion regime of flame-containing cells, and a better method of modeling the effects of the local residual value on the burning velocity. The validation of the revised PRF mechanism shows that the calculated ignition delay matches shock tube data very well. The improvements to the “primary heat release” model based on the cell unburned/burnt volumes more precisely consider the chemical kinetics heat release in unburned regions, and thus are thought to be physically reasonable.

In-cylinder flow and combustion processes simulated with the standard k-ε turbulence model and with an alternative model-employing a non-linear, quadratic equation for the turbulent stresses-are contrasted for both motored and fired engine operation at two loads. For motored operation, the differences observed in the predictions of mean flow development are small and do not emerge until expansion. Larger differences are found in the spatial distribution and magnitude of turbulent kinetic energy. The non-linear model generally predicts lower energy levels and larger turbulent time scales. With fuel injection and combustion, significant differences in flow structure and in the spatial distribution of soot are predicted by the two models. The models also predict considerably different combustion efficiencies and NOx emissions.

In this paper, knock in a Ford single cylinder direct-injection spark-ignition (DISI) engine was modeled and investigated using the KIVA-3V code with a G-equation combustion model coupled with detailed chemical kinetics. The deflagrative turbulent flame propagation was described by the G-equation combustion model. A 22-species, 42-reaction iso-octane (iC8H18) mechanism was adopted to model the auto-ignition process of the gasoline/air/residual-gas mixture ahead of the flame front. The iso-octane mechanism was originally validated by ignition delay tests in a rapid compression machine. In this study, the mechanism was tested by comparing the simulated ignition delay time in a constant volume mesh with the values measured in a shock tube under different initial temperature, pressure and equivalence ratio conditions, and acceptable agreements were obtained.

Engine development is both time consuming and economically straining. Therefore, efforts are being made to optimize the research and development process for new engine technologies. The ability to apply information gained by studying an engine of one size/application to an engine of a completely different size/application would offer savings in both time and money in engine development. In this work, a computational study of diesel engine size-scaling relationships was performed to explore engine scaling parameters and the fundamental engine operating components that should be included in valid scaling arguments. Two scaling arguments were derived and tested: a simple, equal spray penetration scaling model and an extended, equal lift-off length scaling model. The simple scaling model is based on an equation for the conservation of mass and an equation for spray tip penetration developed by Hiroyasu et al. [1].

A level set method (G-equation)-based combustion model incorporating detailed chemical kinetics has been developed and implemented in KIVA-3V for Spark-Ignition (SI) engine simulations for better predictions of fuel oxidation and pollutant formation. Detailed fuel oxidation mechanisms coupled with a reduced NOX mechanism are used to describe the chemical processes. The flame front in the spark kernel stage is tracked using the Discrete Particle Ignition Kernel (DPIK) model. In the G-equation model, it is assumed that after the flame front has passed, the mixture within the mean flame brush tends to local equilibrium. The subgrid-scale burnt/unburnt volumes of the flame containing cells are tracked for the primary heat release calculation. A progress variable concept is introduced into the turbulent flame speed correlation to account for the laminar to turbulent evolution of the spark kernel flame.

The effect of injection rate shape on spray evolution and emission characteristics is investigated and a methodology for active control of fuel injection is proposed. Extensive validation of advanced vaporization and primary jet breakup models was performed with experimental data before studying the effects of systematic changes of injection rate shape. Excellent agreement with the experiments was obtained for liquid and vapor penetration lengths, over a broad range of gas densities and temperatures. Also the predicted flame lift-off lengths of reacting diesel fuel sprays were in good agreement with the experiments. After the validation of the models, well-defined rate shapes were used to study the effect of injection rate shape on liquid and vapor penetration, flame lift-off lengths and emission characteristics.

A level set G-equation model has been developed to model the combustion process in spark ignition engines. The spark ignition process was modeled using an improved version of the Discrete Particle Ignition Kernel (DPIK) model. The two models were implemented into the KIVA-3V code to simulate SI engine combustion under both premixed and direct injection conditions. In the ignition model, the ignition kernel growth is tracked by Lagrangian markers, and spark discharge energy and flow turbulence effects on the kernel growth are considered. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the G-equation model is used to track the mean turbulent flame evolution. When combined with a characteristic time scale combustion model, the models were also used to simulate stratified combustion in DISI engines, where the triple flame structure must be considered.

An improved discrete particle ignition kernel (DPIK) model and the G-equation combustion model have been developed and implemented in KIVA-3V. In the ignition model, the spark ignition kernel growth is tracked by Lagrangian markers and the spark discharge energy and flow turbulence effects on the ignition kernel growth are considered. The predicted ignition kernel size was compared with the available measurements and good agreement was obtained. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the level set method (G-equation) is used to track the mean turbulent flame propagation. It is shown that, by ignoring the detailed turbulent flame brush structure, fine numerical resolution is not needed, thus making the models suitable for use in multidimensional modeling of SI engine combustion.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.