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.

The importance of radiative heat transfer on the combustion and soot formation characteristics under nominal ECN Spray A conditions has been studied numerically. The liquid n-dodecane fuel is injected with 1500 bar fuel pressure into the constant volume chamber at different ambient conditions. Radiation from both gas-phase as well as soot particles has been included and assumed as gray. Three different solvers for the radiative transfer equation have been employed: the discrete ordinate method, the spherical-harmonics method and the optically thin assumption. The radiation models have been coupled with the transported probability density function method for turbulent reactive flows and soot, where unresolved turbulent fluctuations in temperature and composition are included and therefore capturing turbulence-chemistry-soot-radiation interactions. Results show that the gas-phase (mostly CO2 ad H2O species) has a higher contribution to the net radiation heat transfer compared to soot.

Numerical simulations of soot formation were performed for n-dodecane spray using the transported probability density function (TPDF) method. Liquid n-dodecane was injected with 1500 bar fuel pressure into a constant-volume vessel with an ambient temperature, oxygen volume fraction and density of 900 K, 15% and 22.8 kg/m3, respectively. The interaction by exchange with the mean (IEM) model was employed to close the micro-mixing term. The unsteady Reynolds-averaged Navier-Stokes (RANS) equations coupled with the realizable k-ε turbulence model were used to provide turbulence information to the TPDF solver. A 53-species reduced n-dodecane chemical mechanism was employed to evaluate the reaction rates. Soot formation was modelled with an acetylene-based two-equation model which accounts for simultaneous soot particle inception, surface growth, coagulation and oxidation by O2 and OH.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.

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.

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.

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-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 design is a time consuming process in which many costly experimental tests are usually conducted. With increasing prediction ability of engine simulation tools, engine design aided by CFD software is being given more attention by both industry and academia. It is also of much interest to be able to use design information gained from an existing engine design of one size in the design of engines of other sizes to reduce design time and costs. Therefore it is important to study size-scaling relationships for engines over wide range of operating conditions. This paper presents CFD studies conducted for two production diesel engines - a light-duty GM-Fiat engine (0.5L displacement) and a heavy-duty Caterpillar engine (2.5L displacement). Previously developed scaling arguments, including an equal spray penetration scaling model and an extended, equal flame lift-off length scaling model were employed to explore the parametric scaling connections between the two engines.

The influence of the constant C3, which multiplies the mean flow divergence term in the model equation for the turbulent kinetic energy dissipation, is examined in a motored diesel engine for three different swirl ratios and three different spatial locations. Predicted temporal histories of turbulence energy and its dissipation are compared with experimentally-derived estimates. A “best-fit” value of C3 = 1.75, with an approximate uncertainty of ±0.3 is found to minimize the error between the model predictions and the experiments. Using this best-fit value, model length scale behavior corresponds well with that of measured velocity-correlation integral scales during compression. During expansion, the model scale grows too rapidly. Restriction of the model assessment to the expansion stroke suggests that C3 = 0.9 is more appropriate during this period.

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.

The present study uses a numerical model to investigate the effects of flow turbulence on premixed iso-octane HCCI engine combustion. Different levels of in-cylinder turbulence are generated by using different piston geometries, namely a disc-shape versus a square-shape bowl. The numerical model is based on the KIVA code which is modified to use CHEMKIN as the chemistry solver. A detailed reaction mechanism is used to simulate the fuel chemistry. It is found that turbulence has significant effects on HCCI combustion. In the current engine setup, the main effect of turbulence is to affect the wall heat transfer, and hence to change the mixture temperature which, in turn, influences the ignition timing and combustion duration. The model also predicts that the combustion duration in the square bowl case is longer than that in the disc piston case which agrees with the measurements.

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.

Cycle-resolved analysis of velocity data obtained in the re-entrant bowl of a fired high-;speed, direct-injection diesel engine, demonstrates an unambiguous, approximately 100% increase in late-cycle turbulence levels over the levels measured during motored operation. Model predictions of the flow field, obtained employing RNG k-ε turbulence modeling in KIVA-3V, do not capture this increased turbulence. A combined experimental and computational approach is taken to identify the source of this turbulence. The results indicate that the dominant source of the increased turbulence is associated with the formation of an unstable distribution of mean angular momentum, characterized by a negative radial gradient. The importance of this source of flow turbulence has not previously been recognized for engine flows. The enhanced late-cycle turbulence is found to be very sensitive to the flow swirl level.

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.

A correction for the turbulence dissipation, based on non-equilibrium turbulence considerations from rapid distortion theory, has been derived and implemented in combination with the RNG k - ε model in a KIVA-based code. This model correction has been tested and compared with the standard RNG k - ε model for the compression and the combustion phase of two heavy duty DI diesel engines. The turbulence behavior in the compression phase shows clear improvements over the standard RNG k - ε model computations. In particular, the macro length scale is consistent with the corresponding time scale and with the turbulent kinetic energy over the entire compression phase. The combustion computations have been performed with the characteristic time combustion model. With this dissipation correction no additional adjustments of the turbulent characteristic time model constant were necessary in order to match experimental cylinder pressures and heat release rates of the two engines.

A study was performed that takes into account both turbulence and chemical kinetic effects in the numerical simulation of diesel engine combustion in order to better understand the importance of their respective roles at changing operating conditions. An approach was developed which combines the simplicity and low computational and storage requests of the laminar-and-turbulent characteristic-time model with a detailed combustion chemistry model based on well-known simplified mechanisms. Assuming appropriate simplifications such as steady state or equilibrium for most of the radicals and intermediate species, the kinetics of hydrocarbons can be described by means of three overall steps. This approach was integrated in the KIVA-II code. The concept was validated and applied to a single-cylinder, heavy-duty engine. The simulation covers a wide range of operating conditions.