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.

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.

The effects of the temporal and spatial distributions of ignition timings of combustion zones on combustion noise in a Direct Injection Compression Ignition (DICI) engine were studied using experimental tests and numerical simulations. The experiments were performed with different fuel injection strategies on a heavy-duty diesel engine. Cylinder pressure was measured with the sampling intervals of 0.1°CA in order to resolve noise components. The simulations were performed using the KIVA-3V code with detailed chemistry to analyze the in-cylinder ignition and combustion processes. The experimental results show that optimal sequential ignition and spatial distribution of combustion zones can be realized by adopting a two-stage injection strategy in which the proportion of the pilot injection fuel and the timings of the injections can be used to control the combustion process, thus resulting in simultaneously higher thermal efficiency and lower noise emissions.

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.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline (termed GDICI for Gasoline Direct-Injection Compression Ignition) in the low temperature combustion (LTC) regime is presented. As an aid to plan engine experiments at full load (16 bar IMEP, 2500 rev/min), exploration of operating conditions was first performed numerically employing a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion. Operation ranges of a light-duty diesel engine operating with GDICI combustion with constraints of combustion efficiency, noise level (pressure rise rate) and emissions were identified as functions of injection timings, exhaust gas recirculation rate and the fuel split ratio of double-pulse injections.

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.

This paper concerns a study of an innovative concept to control HCCI combustion in diesel-fueled engines. The concept consists in forming a pre-compressed homogeneous charge outside the cylinder and in gradually admitting it into the cylinder during the combustion process. This new combustion concept has been called Homogeneous Charge Progressive Combustion (HCPC). CFD analysis was conducted to understand the feasibility of the HCPC concept and to identify the parameters that control and influence this novel HCCI combustion. A CFD code with detailed kinetic chemistry (AVL FIRE) was used in the study. Results in terms of pressure, heat release rate, temperature, and emissions production are presented that demonstrate the validity of the HCPC combustion concept.

In the present paper optimization tools are used to recommend low-emission engine combustion chamber designs, spray targeting and swirl ratio levels for a heavy-duty diesel engine operated at high-load. The study identifies aspects of the combustion and pollution formation that are affected by mixing processes, and offers guidance for better matching of the piston geometry with the spray plume geometry for enhanced mixing. By coupling a GA (genetic algorithm) with the KIVA-CFD code, and also by utilizing an automated grid generation technique, multi-objective optimizations with goals of low emissions and fuel economy were achieved. Three different multi-objective genetic algorithms including a Micro-Genetic Algorithm (μGA), a Nondominated Sorting Genetic Algorithm II (NSGA II) and an Adaptive Range Multi-Objective Genetic Algorithm (ARMOGA) were compared for conducting the optimization under the same conditions.

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.

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].

High-speed video imaging in a swirl-supported (Rs = 1.7), direct-injection heavy-duty diesel engine operated with moderate-to-high EGR rates reveals a distinct correlation between the spatial distribution of luminous soot and mean flow vorticity in the horizontal plane. The temporal behavior of the experimental images, as well as the results of multi-dimensional numerical simulations, show that this soot-vorticity correlation is caused by the presence of a greater amount of soot on the windward side of the jet. The simulations indicate that while flow swirl can influence pre-ignition mixing processes as well as post-combustion soot oxidation processes, interactions between the swirl and the heat release can also influence mixing processes. Without swirl, combustion-generated gas flows influence mixing on both sides of the jet equally. In the presence of swirl, the heat release occurs on the leeward side of the fuel sprays.

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.

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

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.

In-cylinder fluid velocity is measured in an optically accessible, fired HSDI engine at idle. The velocity field is also calculated, including the full induction stroke, using multi-dimensional fluid dynamics and combustion simulation models. A detailed comparison between the measured and calculated velocities is performed to validate the computed results and to gain a physical understanding of the flow evolution. Motored measurements are also presented, to clarify the effects of the fuel injection process and combustion on the velocity field evolution. The calculated mean in-cylinder angular momentum (swirl ratio) and mean flow structures prior to injection agree well with the measurements. Modification of the mean flow by fuel injection and combustion is also well captured.

The effects of numerical methodology in defining the initial conditions and simulating the compression stroke in D.I. diesel engine CFD computations are studied. Lumped and pointwise approaches were adopted in assigning the initial conditions at IVC. The lumped approach was coupled with a two-dimensional calculation of the compression stroke. The pointwise methodology was based on the results of an unsteady calculation of the intake stroke performed by using the STAR-CD code in the realistic engine and port geometry. Full engine and 60 deg. sector meshes were used in the compression stroke calculations in order to check the accuracy of the commonly applied axi-symmetric fluid dynamics assumption. Analysis of the evolution of the main fluid dynamics parameters revealed that local conditions at the time of injection strongly depend on the numerical procedure adopted.

To model sprays from pressure-swirl atomizers, the connection between the injector and the downstream spray must be considered. A new model for pressure-swirl atomizers is presented which assumes little knowledge of the internal details of the injector, but instead uses available observations of external spray characteristics. First, a correlation for the exit velocity at the injector exit is used to define the liquid film thickness. Next, the film must be modeled as it becomes a thin, liquid sheet and breaks up, forming ligaments and droplets. A linearized instability analysis of the breakup of a viscous, liquid sheet is used as part of the spray boundary condition. The spray angle is estimated from spray photographs and patternator data. A mass averaged spray angle is calculated from the patternator data and used in some of the calculations.

In multidimensional modeling, fuels have been represented predominantly by single components, such as octane for gasoline. Several bicomponent studies have been performed, but these are still limited in their ability to represent real fuels, which are blends of as many as 300 components. This study outlines a method by which the fuel composition is represented by a distribution function of the fuel molecular weight. This allows a much wider range of compositions to be modeled, and only requires including two additional “species” besides the fuel, namely the mean and second moment of the distribution. This approach has been previously presented but is applied here to multidimensional calculations. Results are presented for single component droplet vaporization for comparison with single component fuel predictions, as well as results for a multicomponent gasoline and a diesel droplet.