Search Results

This paper implements a coupled approach to integrate the internal nozzle flow and the ensuing fuel spray using a Volume-of-Fluid (VOF) method in the CONVERGE CFD software. A VOF method was used to model the internal nozzle two-phase flow with a cavitation description closed by the homogeneous relaxation model of Bilicki and Kestin [1]. An Eulerian single velocity field approach by Vallet et al. [2] was implemented for near-nozzle spray modeling. This Eulerian approach considers the liquid and gas phases as a complex mixture with a highly variable density to describe near nozzle dense sprays. The mean density is obtained from the Favreaveraged liquid mass fraction. The liquid mass fraction is transported with a model for the turbulent liquid diffusion flux into the gas.

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.

This work involves modeling internal and near-nozzle flows of a gasoline direct injection (GDI) nozzle. The Engine Combustion Network (ECN) Spray G condition has been considered for these simulations using the nominal geometry of the Spray G injector. First, best practices for numerical simulation of the two-phase flow evolution inside and the near-nozzle regions of the Spray G injector are presented for the peak needle lift. The mass flow rate prediction for peak needle lift was in reasonable agreement with experimental data available in the ECN database. Liquid plume targeting angle and liquid penetration estimates showed promising agreement with experimental observations. The capability to assess the influence of different thermodynamic conditions on the two-phase flow nature was established by predicting non-flashing and flashing phenomena.

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

Consideration of the heat transfer effects in low-heat-rejection engines has prompted further study into engine heat transfer phenomena. In a previous study, an approximate solution of the one-dimensional energy equation was acquired for transient, compressible, low-Mach number, turbulent boundary layers typical of those found in engines. The current study shows that an approximate solution of the one-dimensional energy equation with arbitrarily-distributed heat release can also be obtained. Using this model, the effects of high temperature walls, combustion, and autoignition on heat transfer can be studied. In the case of high temperature walls, the model predicts the expected behavior unless the quench distance gets very small. For combustion, the reaction must occur close to the wall for a direct effect on the heat transfer to be observed. With autoignition, instantaneous values of heat flux reach levels as high as 6 MW/m2, and oscillate in phase with the pressure wave.

Engine design is a time consuming process in which many costly experimental tests are usually conducted. With increasing prediction ability of engine simulation tools, engine design aided by CFD software is being given more attention by both industry and academia. It is also of much interest to be able to use design information gained from an existing engine design of one size in the design of engines of other sizes to reduce design time and costs. Therefore it is important to study size-scaling relationships for engines over wide range of operating conditions. This paper presents CFD studies conducted for two production diesel engines - a light-duty GM-Fiat engine (0.5L displacement) and a heavy-duty Caterpillar engine (2.5L displacement). Previously developed scaling arguments, including an equal spray penetration scaling model and an extended, equal flame lift-off length scaling model were employed to explore the parametric scaling connections between the two engines.

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.

Multi-Dimensional modeling was carried out for a Mercury Marine two-stroke DISI engine. Recently developed spray, ignition, and combustion models were applied to medium load cases with an air-fuel ratio of 30:1. Three injection timings, 271, 291 and 306 ATDC were selected to investigate the effects of the injection timing on mixture formation, ignition and combustion. The results indicate that at this particular load condition, earlier injection timing allows more fuel to evaporate. However, because the fuel penetrates further toward the piston, a leaner mixture is created near the spark plug; thus, a slower ignition process with a weaker ignition kernel was found for the SOI 271 ATDC case. The measured and computed combustion results such as average in-cylinder pressure and NOx are in good agreements. The later injection case produces lower NOx emission and higher CO emission; this is due to poor mixing and is in agreement with experimental measurements.

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

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.

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.

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.

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.

The 3-Zones Extended Coherent Flame Model (ECFM3Z) and the Tabulated Kinetics for Ignition (TKI) auto-ignition model are widely used for RANS simulations of reactive flows in Diesel engines. ECFM3Z accounts for the turbulent mixing between one zone that contains compressed air and EGR and another zone that contains evaporated fuel. These zones mix to form a reactive zone where combustion occurs. In this mixing zone TKI is applied to predict the auto-ignition event, including the ignition delay time and the heat release rate. Because it is tabulated, TKI can model complex fuels over a wide range of engine thermodynamic conditions. However, the ECFM3Z/TKI combustion modeling approach requires an efficient predictive spray injection calculation. In a Diesel direct injection engine, the turbulent mixing and spray atomization are mainly driven by the liquid/gas coupling phenomenon that occurs at moving liquid/gas interfaces.

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.

Relationships between radiant emissions, as measured by an in-cylinder optical sensor, and spark-ignition engine combustion parameters are presented for possible use in engine combustion diagnostics and future engine control strategies. A monochromatic gas radiation model, developed in a previous study, was used to derive a series of relationships between the measured radiant emission characteristics and several spark-ignition engine combustion parameters, such as the amplitude and phasing of the peak heat-release rate, combustion duration, IMEP, NOx emission, pressure, trapped mass and exhaust-gas temperature. In addition, many engine parameters of interest can be estimated indirectly from the radiation signal using empirical models. Correlations of air-fuel ratio and exhaust emissions are presented which contain a combination of radiant emission parameters and known base-engine operating parameters, such as intake manifold pressure, etc.

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.

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.

Current demands for high fuel efficiency and low emissions in automotive powerplants have drawn attention to the two-stroke engine configuration. The present study measured trapping and scavenging efficiencies of a firing two-stroke spark-ignition engine by in-cylinder gas composition analysis. Intermediate results of the procedure included the trapped air-fuel ratio and residual exhaust gas fraction. Samples, acquired with a fast-acting electromagnetic valve installed in the cylinder head, were taken of the unburned mixture without fuel injection and of the burned gases prior to exhaust port opening, at engine speeds of 1000 to 3000 rpm and at 10 to 100% of full load. A semi-empirical, zero-dimensional scavenging model was developed based on modification of the non-isothermal, perfect-mixing model. Comparison to the experimental data shows good agreement.