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.

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.

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.

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.

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.

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.

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

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.

Combustion simulations utilizing the characteristic time combustion model have been performed for four DI diesel engines ranging in size from heavy-duty to large-bore designs. It has been found that the pre-factor to the turbulent characteristic time acts as a scaling parameter between the engines. This phenomenon is explained in terms of the non-equilibrium behavior of the turbulent time and length scales, as is encountered in the rapidly distorting, spray-induced flows of DI diesel engines. In fact, the equilibrium assumption between turbulence production and dissipation, which forms the basis for the employed k-ε-type turbulence models, does not hold in these situations. For such flows, the real turbulent dissipation time scale is locally proportional to the turbulent characteristic time scale which is determined by a typical eddy turnover time.

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.

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.

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

An integrated numerical model has been developed for diesel engine computations based on the KIVA-II code. The model incorporates a modified RNG k-ε, turbulence model, a ‘wave’ breakup spray model, the Shell ignition model, the laminar-and-turbulent characteristic-time combustion model, a crevice flow model, a spray/wall impingement model that includes rebounding and breaking-up drops, and other improved submodels in the KIVA code. The model was validated and applied to model successfully different types of diesel engines under various operating conditions. These engines include a Caterpillar engine with different injection pressures at different injection timings, a small Tacom engine at different loads, and a Cummins engine modified by Sandia for optical experiments. Good levels of agreement in cylinder pressures and heat release rate data were obtained using the same computer model for all engine cases.

Particle Image Velocimetry (PIV) was used to make gas velocity and turbulence measurements in a motored diesel engine. The experiments were conducted using a single-cylinder version of the Caterpillar 3406 production engine. One of the exhaust valves and the fuel injector port were used to provide optical access to the combustion chamber so that modifications to the engine geometry were minimal, and the results are representative of the actual engine. Measurements of gas velocity were made in a plane in the piston bowl using TiO2 seed particles. The light sheet necessary for PIV was formed by passing the beam from a Nd:YAG laser through the injector port and reflecting the beam off a conical mirror at the center of the piston. PIV data was difficult to obtain due to significant out-of-plane velocities. However, data was acquired at 25° and 15° before top dead center of compression at 750 rev/min.

Combustion and flow were calculated in a spark-ignited two-stroke crankcase-scavenged engine using a laminar and turbulent characteristic-time combustion submodel in the three-dimensional KIVA code. Both premixed-charge and fuel-injected cases were examined. A multi-cylinder engine simulation program was used to specify initial and boundary conditions for the computation of the scavenging process. A sensitivity study was conducted using the premixed-charge engine data. The influence of different port boundary conditions on the scavenging process was examined. At high delivery ratios, the results were insensitive to variations in the scavenging flow or residual fraction details. In this case, good agreement was obtained with the experimental data using an existing combustion submodel, previously validated in a four-stroke engine study.