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.

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.

An experimental study was conducted on an air cooled high-speed, direct-injection diesel generator that investigated the use of gasoline in a dual fuel PCCI strategy. The single-speed generator used in the study has an effective compression ratio of 17 and runs at 3600 rev/min. Varying amounts of gasoline were introduced into the combustion chamber through a port injection system. The generator uses an all-mechanical diesel fuel injection system that has a fixed injection timing. The experiments explored variable intake temperatures and fuel split quantities to investigate different combustion phasing regimes. Results from the study showed low combustion efficiency at low load. Low load operation was also characterized by high levels of HC and CO (in excess of 20 g/kwh and 50 g/kwh respectively). Operation at 75% load was more efficient, and displayed three different combustion regimes that are possible with PIG (port injected gasoline) dual fuel PCCI.

The possibility to operate current diesel engines in dual-fuel mode with the addition of an alternative fuel is fundamental to accelerate the energy transition to achieve carbon neutrality. The simulation of the dual- fuel combustion process with 0D/1D combustion models is fundamental for the performance prediction, but still particularly challenging, due to chemical interactions of the mixture. The authors defined a novel data-driven workflow for the development of combustion reaction mechanisms and used it to generate a dual-fuel mechanism for Ammonia and Diesel Primary Reference Fuels (DPRF) suitable for efficient combustion simulations in heavy duty engines, with variable cetane number Diesel fuels. A baseline reaction mechanism was created by merging the detailed ammonia mechanism by Glarborg et al. with reaction pathways for n- hexadecane and 2,2,4,4,6,8,8-heptamethylnonane from a well-established multi-component fuel mechanism.

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.

Contamination from water and particulate matter in diesel or aviation fuel can cause complications even under normal operating conditions. Fuel contamination becomes a major problem in air transportation, where engine flame-out due to injection system blockage or malfunction might have catastrophic consequences. The presence of contaminants in fuel has been linked to several incidents including commercial and military aircraft over the past decade. Unless they are detected and separated from the fuel, water or solid particles have several ways to reach the engine and cause troubles. Aviation fuel is commonly stored in large tanks and transferred frequently before it reaches the aircraft, increasing the risk for contamination. At the same time, gas bubbles may be present in the system. While harmless to the aircraft, gas bubbles have been the reason why contamination monitoring systems would fail, as the system would be triggered by gas bubbles, instead of detrimental contaminants.

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.

The present study investigated HCCI combustion in a heavy-duty diesel engine both experimentally and numerically. The engine was equipped with a hollow-cone pressure-swirl injector using gasoline direct injection. Characteristics of HCCI combustion were obtained by very early injection with a heated intake charge. Experimental results showed an increase in NOx emission and a decrease in UHC as the injection timing was retarded. It was also found that optimization can be achieved by controlling the intake temperature together with the start-of-injection timing. The experiments were modeled by using an engine CFD code with detailed chemistry. The CHEMKIN code was implemented into KIVA-3V such that the chemistry and flow solutions were coupled. The model predicted ignition timing, cylinder pressure, and heat release rates reasonably well. The NOx emissions were found to increase as the injection timing was retarded, in agreement with experimental results.

An image acquisition-and-processing camera system was developed for in-cylinder diagnostics of a single-cylinder heavy duty diesel engine. The engine was equipped with an electronically-controlled common-rail fuel injection system that allowed both single and split (multiple) injections to be studied. The imaging system uses an endoscope to acquire luminous flame images from the combustion chamber and ensures minimum modification to the engine geometry. The system also includes an optical linkage, an image intensifier, a CID camera, a frame grabber, control circuitry and a computer. Experiments include both single and split injection cases at 90 MPa and 45 MPa injection pressures at 3/4 load and 1600 rev/min with simulated turbocharging. For the single injection at high injection pressure (90 MPa) the results show that the first luminous emissions from the ignition zone occur very close to the injector exit followed by rapid luminous flame spreading.

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 continuous multi-component fuel evaporation model has been integrated with an improved G-equation combustion and detailed chemical kinetics model. The integrated code has been successfully used to simulate a gasoline direct injection engine. In the multi-component fuel model, the theory of continuous thermodynamics is used to model the properties and composition of multi-component fuels such as gasoline. In the improved G-equation combustion model a flamelet approach based on the G-equation is used that considers multi-component fuel effects. To precisely calculate the local and instantaneous residual which has a great effect on the laminar flame speed, a “transport equation residual” model is used. A Damkohler number criterion is used to determine the combustion mode in flame containing cells.

A numerical study was performed to compare the formation of nitric oxide (NO) and nitrogen dioxide (NO₂), collectively termed NOx, resulting from biodiesel and diesel combustion in an internal combustion engine. It has been shown that biodiesel tends to increase NOx compared to diesel, and to-date, there is no widely accepted explanation. Many factors can lead to increased NOx formation and it was of interest to determine if fuel chemistry plays a significant role. Therefore, in order to isolate the fuel chemistry from mixing processes typical in a compression ignition engine, sprays were not considered in the present investigation. The current study compares the NOx formation of surrogates for biodiesel (as represented by methyl butanoate and n-heptane) and diesel (n-heptane) under completely homogeneous conditions. Combustion of each fuel was simulated using the Senkin code for both an adiabatic, constant volume reactor, and an adiabatic, single-zone HCCI engine 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.

An improved discrete particle ignition kernel (DPIK) model and the G-equation combustion model have been developed and implemented in KIVA-3V. In the ignition model, the spark ignition kernel growth is tracked by Lagrangian markers and the spark discharge energy and flow turbulence effects on the ignition kernel growth are considered. The predicted ignition kernel size was compared with the available measurements and good agreement was obtained. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the level set method (G-equation) is used to track the mean turbulent flame propagation. It is shown that, by ignoring the detailed turbulent flame brush structure, fine numerical resolution is not needed, thus making the models suitable for use in multidimensional modeling of SI engine combustion.

A level set G-equation model has been developed to model the combustion process in spark ignition engines. The spark ignition process was modeled using an improved version of the Discrete Particle Ignition Kernel (DPIK) model. The two models were implemented into the KIVA-3V code to simulate SI engine combustion under both premixed and direct injection conditions. In the ignition model, the ignition kernel growth is tracked by Lagrangian markers, and spark discharge energy and flow turbulence effects on the kernel growth are considered. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the G-equation model is used to track the mean turbulent flame evolution. When combined with a characteristic time scale combustion model, the models were also used to simulate stratified combustion in DISI engines, where the triple flame structure must be considered.

Detailed chemical kinetics was implemented in the KIVA-3V multidimensional CFD code to study the combustion process in Homogeneous Charge Compression Ignition (HCCI) engines. The CHEMKIN code was implemented such that the chemistry and flow solutions were coupled. Detailed reaction mechanisms were used to simulate the fuel chemistry of ignition and combustion. Effects of turbulent mixing on the reaction rates were also considered. The model was validated using the experimental data from two modified heavy-duty diesel engines, including a Volvo engine and a Caterpillar engine operated at the HCCI mode. The results show that good levels of agreement were obtained using the present KIVA/CHEMKIN model for a wide range of engine conditions, including various fuels, injection systems, engine speeds, and EGR levels. Ignition timings were predicted well without the need to adjust any kinetic constants.

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.