Search Results

A three-dimensional homogeneous equilibrium model (HEM) has been developed and implemented into an engine computational fluid dynamics (CFD) code KIVA-3V. The model was applied to simulate cavitating flow within injector nozzle passages. The effects of nozzle passage geometry and injection conditions on the development of cavitation zones and the nozzle discharge coefficient were investigated. Specifically, the effects of nozzle length (L/D ratio), nozzle inlet radius (R/D ratio) and K or KS factor (nozzle passage convergence) were simulated, and the effects of injection and chamber pressures, and time-varying injection pressure were also investigated. These effects are well captured by the nozzle flow model, and the predicted trends are consistent with those from experimental observations and theoretical analyses.

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.

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.

An emissions and performance study was performed to show the effects of injection pressure, nozzle hole inlet condition (sharp and rounded edge) and nozzle included spray angle on particulate, NOx, and BSFC. The tests were conducted on a fully instrumented single-cylinder version of the Caterpillar 3406 heavy duty engine at 75% and 25% load at 1600 RPM. The fuel system consisted of an electronically controlled, hydraulically actuated, unit injector capable of injection pressures up to 160 MPa. Particulate versus NOx trade-off curves were generated for each case by varying the injection timing. The 75% load results showed the expected decrease in particulate and flattening of the trade-off curve with increased injection pressure. However, in going from 90 to 160 MPa, the timing had to be retarded to maintain the same NOx level, and this resulted in a 1 to 2% increase in BSFC. The rounded edged nozzles were found to have an increased discharge coefficient.

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.

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 goal of this research is to investigate the physical parameters of stoichiometric operation of a diesel engine under a light load operating condition (6∼7 bar IMEP). This paper focuses on improving the fuel efficiency of stoichiometric operation, for which a fuel consumption penalty relative to standard diesel combustion was found to be 7% from a previous study. The objective is to keep NOx and soot emissions at reasonable levels such that a 3-way catalyst and DPF can be used in an aftertreatment combination to meet 2010 emissions regulation. The effects of intake conditions and the use of group-hole injector nozzles (GHN) on fuel consumption of stoichiometric diesel operation were investigated. Throttled intake conditions exhibited about a 30% fuel penalty compared to the best fuel economy case of high boost/EGR intake conditions. The higher CO emissions of throttled intake cases lead to the poor fuel economy.

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.

Fuel Injection System (FIS) research on injection pressure, timing control, and rate shaping, and studies on the modeling of injector nozzle flows and their effect on fuel spray characteristics are usually conducted separately. Only recently has the fuel injection and spraying process been studied as a complete system, i.e., including both the high-pressure fuel delivery and its effect on the nozzle flow characteristics, including nozzle cavitation. A methodology for coupling the fuel injection system and its effect on spray characteristics is presented here. The method is applied to an example case of a conventional pump-line-nozzle system. Mathematical models for characterizing the flows from the pump to the nozzle are formulated and solved using the Method of Characteristics and finite difference techniques. The nozzle internal flow is modeled using zero-dimensional flow models, in which the nozzle cavitation and its effect on the nozzle exit flow are accounted for.

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.

Managing the injection rate profile is a powerful tool to control engine performance and emission levels. In particular, Common Rail (C.R.) injection systems allow an almost completely flexible fuel injection event in DI-diesel engines by permitting a free mapping of the start of injection, injection pressure, rate of injection and, in the near future, multiple injections. This research deals with the development of a network-based numerical tool for understanding operating condition limits of the Common Rail injector. The models simulate the electro-fluid-mechanical behavior of the injector accounting for cavitation in the nozzle holes. Validation against experiments has been performed. The model has been used to provide insight into the operating conditions of the injector and in order to highlight the application to injection system design.

In the Atomization regime, liquid jets breakup either within the nozzle or immediately upon entering the chamber gas and drops much smaller than the jet diameter are formed. The mechanism of Atomization, which is presently unknown, was investigated by the simultaneous use of two photographic techniques. The initial transient was observed with a 106 frames/s camera and the steady state by a technique similar to spark photography. The experiment range was: liquid pressure 500 to 2500 psia; five mixtures of water and glycerol to vary the liquid viscosity; air, nitrogen, helium, and xenon at up to 600 psia as chamber gases to separate gas pressure from gas density effects; and 14 nozzle designs. Not changed were the temperature (room value), the nozzle diameter (340 μ), and the surface tension (70 dyne/cm).

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.