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.

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

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.

Two-dimensional computations of premixed-charge engine combustion were made using the KIVA-II code. The purpose of the study was to assess the influence of heat transfer and turbulence model boundary conditions on engine combustion predictions. Combustion was modeled using a laminar- and turbulent-characteristic-time model. Flow through the piston-cylinder-ring crevice was accounted for using a phenomenological crevice-flow model. The predictions were compared to existing cylinder pressure and wall heat transfer experimental data under motoring and fired conditions, at two engine speeds. Two different wall heat transfer model formulations were considered. The first is the standard wall function method. The second is based on solutions to the one-dimensional unsteady energy equation, formulated such that the standard wall function method is recovered in the quasi-steady limit. Turbulence was modeled using the standard k-ε turbulence model equations.

Fuel spray atomization and breakup processes within a direct-injection spark-ignition (DISI) engine and outside the engine were modeled using a modified KIVA-3V code with improved spray models. The structures of the predicted sprays were qualitatively compared with planar images. The considered sprays were created by a prototype pressure-swirl injector and the planar images were obtained by laser sheet imaging in an optical DISI engine. In the out-of-engine case, the spray was injected into atmospheric air, and was modeled in a two dimensional bomb. In the engine case, the injection started from 270° ATDC, and full 3-D computations in the same engine were performed. In both cases, two liquid injection pressure conditions were applied, that is, 3.40 MPa and 6.12 MPa. The model gives good prediction of the tip penetration, and external spray shape, but the internal structure prediction has relatively lower accuracy, especially near the spray axis.

Multidimensional computations were made of spark-ignited premixed-charge combustion in a pancake-combustion-chamber engine with a centrally located spark plug and in two pent-roof-chamber engines, one with a central spark plug and the other with dual lateral spark plugs. A global combustion submodel was used that accounts for laminar kinetics and turbulent mixing effects. The predictions were compared with available measurements in the pancake-chamber engine over a range of loads, speeds, and equivalence ratios. In all cases the computed and measured cylinder pressures agreed well in trends and magnitudes (within 8%) for the entire duration of combustion. Fair agreements were also obtained between predicted and measured values of wall heat flux and emission index of nitric oxide. In the pent-roof-chamber engines the predicted maximum cylinder pressures also agreed well with measurements (within 12%) in cases with MBT (Minimum spark advance for Best Torque) or advanced spark timing.

A new ignition and combustion model has been developed and tested for use in premixed spark-ignition engines. The ignition model is referred to as the Discrete Particle Ignition Kernel (DPIK) model, and it uses Lagrangian markers to track the flame-front growth. The model includes the effects of electrode heat transfer on the early flame kernel growth process, and it is used in conjunction with a characteristic-time-scale combustion model once the ignition kernel has grown to a size where the effects of turbulence on the flame must be considered. A new term which accounts for the effect of air-fuel ratio, was added to the combustion model for modeling combustion in very lean and very rich mixtures. The flame kernel size predicted by the DPIK model was compared with measurements of Maly and Vogel. Furthermore, predictions of the electrode heat transfer were compared with data of Kravchik and Heywood. In both comparisons the model predictions were in good agreement with the experiments.

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.

A study was performed to correlate the Sauter Mean Diameter (SMD), NOx and particulate emissions of a direct injection diesel engine with various injection pressures and different nozzle geometry. The spray experiments and engine emission tests were conducted in parallel using the same fuel injection system and same operating conditions. With high speed photography and digital image analysis, a light extinction technique was used to obtain the spray characteristics which included spray tip penetration length, spray angle, and overall average SMD for the entire spray. The NOx and particulate emissions were acquired by running the tests on a fully instrumented Caterpillar 3406 heavy duty engine. Experimental results showed that for higher injection pressures, a smaller SMD was observed, i.e. a finer spray was obtained. For this case, a higher NOx and lower particulate resulted.

Multidimensional modeling is used to study air-fuel mixing in a direct-injection spark-ignition engine. Emphasis is placed on the effects of the start of fuel injection on gas/spray interactions, wall wetting, fuel vaporization rate and air-fuel ratio distributions in this paper. It was found that the in-cylinder gas/spray interactions vary with fuel injection timing which directly impacts spray characteristics such as tip penetration and spray/wall impingement and air-fuel mixing. It was also found that, compared with a non-spray case, the mixture temperature at the end of the compression stroke decreases substantially in spray cases due to in-cylinder fuel vaporization. The computed trapped-mass and total heat-gain from the cylinder walls during the induction and compression processes were also shown to be increased in spray cases.

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.

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.

Practical diesel fuel sprays under high-pressure conditions were investigated by using multidimensional modeling combined with continuous thermodynamics and high-pressure multicomponent fuel vaporization models. Transport equations, which are general for the moments of the distributions and independent of the distribution function, are derived for the continuous system consisting of the both gas and liquid phases. A general treatment of the vapor-liquid equilibrium (VLE) is conducted, and the Peng-Robinson Equation of State (EOS) is used to find the surface equilibrium composition. Relations for the properties of the continuous species are formulated. The KH-RT model is used for spray breakup prediction. The fuel droplets are assumed to be well mixed with uniform temperature and composition within each droplet. The turbulent flow field is calculated using the RNG k -ε turbulence model.

A three-dimensional computer code (KIVA) is being modified to include state-of-the-art submodels for diesel engine flow and combustion: spray atomization, drop breakup/coalescence, multi-component fuel vaporization, spray/wall interaction, ignition and combustion, wall heat transfer, unburned HC and NOx formation, soot and radiation and the intake flow process. Improved and/or new submodels which have been completed are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NOx, and spray/wall impingement with rebounding and sliding drops.

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.

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.

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.

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.