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.

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.

A new technique which is based on optoacoustic phenomena has been developed for measuring in-cylinder gas temperature and turbulent diffusivity. In the experiments, a high energy Nd:YAG pulsed laser beam was focused to cause local ionization of air at a point in the combustion chamber. This initiates a shock wave and creates a hot spot. The local temperature and turbulent diffusivity are determined by monitoring the shock propagation and the hot spot growth, respectively, with a schlieren photography system. In order to assess the validity and accuracy of the measurements, the technique was also applied to a turbulent jet. The temperature measurements were found to be accurate to within 3%. Results from the turbulent jet measurements also showed that the growth rate of the hot spot diameter can be used to estimate the turbulent diffusivity. In-cylinder gas temperature measurements were made in a motored single cylinder Caterpillar diesel engine, modified for optical access.

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.

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.

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.

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.

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.

The interaction of fuel sprays and in-cylinder flow in direct-injection engines is expected to alter kinetic energy and integral length scales at least during some portions of the engine cycle. High-speed particle image velocimetry was implemented in an optical four-valve, pent-roof spark-ignition direct-injection single-cylinder engine to quantify this effect. Non-firing motored engine tests were performed at 1300 RPM with and without fuel injection. Two fuel injection timings were investigated: injection in early intake stroke represents quasi-homogenous engine condition; and injection in mid compression stroke mimics the stratified combustion strategy. Two-dimensional crank angle resolved velocity fields were measured to examine the kinetic energy and integral length scale through critical portions of the engine cycle. Reynolds decomposition was applied on the obtained engine flow fields to extract the fluctuations as an indicator for the turbulent flow.

The turbulent in-cylinder air flow and the unsteady high-pressure fuel injection lead to a highly transient air fuel mixing process in spark-ignition direct-injection (SIDI) engines, which is the leading cause for combustion cycle-to-cycle variation (CCV) and requires further investigation. In this study, crank-angle resolution particle image velocimetry (PIV) was employed to simultaneously measure the air flow and fuel spray structure at 1300 rpm in an optically accessible single-cylinder SIDI engine. The measurement was conducted at the center tumble plane of the four-valve pent-roof engine, bisecting the spark plug and fuel injector. 84 consecutive cycles were recorded for three engine conditions, i.e. (1) none-fueled motored condition, (2) homogeneous-charge mode with start of injection (SOI) during intake (50 crank-angle degree (CAD) after top dead center exhaust, aTDCexh), and (3) stratified-charge mode with SOI during mid compression (270 aTDCexh).

Intake flow simulations were carried out for a prototype DISI engine using the standard k-ε model and the RNG k-ε model. The results were compared with PTV (transient water analog) measurements. The study was focused on low load operations with engine speed at 400 rev/min. Two cases were studied, a single intake case in which one intake port was blocked and a dual intake port case. In the computations, the results show that the standard k-ε model tends to produce higher turbulence levels when turbulence is generated and decays faster when turbulence dissipates. Different turbulence models predict almost the same flow structures. However, the effects of the turbulence model on the predicted tumble and swirl ratios are significant. The TKE distributions at BDC predicted by the two models are also different. The standard k-ε model seems to be more diffusive. Good agreements with PTV data were obtained in the single valve case with the RNG k-ε model.

The influence of the bulk in-cylinder flow on the spray evolution, evaporation, fuel-air mixing and subsequent flame propagation has been studied in an optical SIDI engine. Quantitative LIF imaging of equivalence ratios was used to characterize the mixture formation and the influence of the local equivalence ratio at the time of spark on the flame propagation. Two extreme bulk flow conditions - high and low swirl - were investigated and pronounced differences in mixture homogeneity and flame propagation were found and characterized.

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.

A primary goal of large eddy simulation, LES, is to capture in-cylinder cycle-to-cycle variability, CCV. This is a first step to assess the efficacy of 35 consecutive computed motored cycles to capture the kinetic energy in the TCC-III engine. This includes both the intra-cycle production and dissipation as well as the kinetic energy CCV. The approach is to sample and compare the simulated three-dimensional velocity equivalently to the available two-component two-dimensional PIV velocity measurements. The volume-averaged scale-resolved kinetic energy from the LES is sampled in three slabs, which are volumes equal to the two axial and one azimuthal PIV fields-of-view and laser sheet thickness. Prior to the comparison, the effects of sampling a cutting plane versus a slab and slabs of different thicknesses are assessed. The effects of sampling only two components and three discrete planar regions is assessed.

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.

Multidimensional computations were carried out to simulate the in-cylinder fuel/air mixing process of a direct-injection spark-ignition engine using a modified version of the KIVA-3 code. A hollow cone spray was modeled using a Lagrangian stochastic approach with an empirical initial atomization treatment which is based on experimental data. Improved Spalding-type evaporation and drag models were used to calculate drop vaporization and drop dynamic drag. Spray/wall impingement hydrodynamics was accounted for by using a phenomenological model. Intake flows were computed using a simple approach in which a prescribed velocity profile is specified at the two intake valve openings. This allowed three intake flow patterns, namely, swirl, tumble and non-tumble, to be considered. It was shown that fuel vaporization was completed at the end of compression stroke with early injection timing under the chosen engine operating conditions.

A numerical study of air-fuel mixing in a direct-injection spark-ignition engine was carried out. In this paper, the numerical models are described and grid generation methods to represent a realistic port-valve-chamber geometry is discussed. To model a vaporizing hollow-cone spray resulting from an automotive pressure-swirl injector, a newly developed sheet spray atomization model was used to compute the processes of disintegration of the liquid sheet and breakup of the subsequent drops. Computations were performed of a particular 4-valve pent-roof engine configuration in which the intake process and an early fuel injection scheme were considered. After an analysis of the intake-generated flow structures in this engine configuration, the spray behavior and the spatial and temporal evolution of fuel liquid and vapor phases are characterized.

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.