Search Results

A computational methodology has been developed for loss prediction in intake regions of internal combustion engines. The methodology consists of a hierarchy of four major tasks: (1) proper computational modeling of flow physics; (2) exact geometry and high quality and generation; (3) discretization schemes for low numerical viscosity; and (4) higher order turbulence modeling. Only when these four tasks are dealt with properly will a computational simulation yield consistently accurate results. This methodology, which is has been successfully tested and validated against benchmark quality data for a wide variety of complex 2-D and 3-D laminar and turbulent flow situations, is applied here to a loss prediction problem from industry. Total pressure losses in the intake region (inlet duct, manifold, plenum, ports, valves, and cylinder) of a Caterpillar diesel engine are predicted computationally and compared to experimental data.

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.

This work presents turbulent premixed combustion modeling in spark ignition engines using G-equation based turbulent combustion model. In present study, a turbulent flame speed expression proposed and validated in recent years by two co-authors of this paper is applied to the combustion simulation of spark ignition engines. This turbulent flame speed expression has no adjustable parameters and its constants are closely tied to the physics of scalar mixing at small scales. Based on this flame speed expression, a minor modification is introduced in this paper considering the fact that the turbulent flame speed changes to laminar flame speed if there is no turbulence. This modified turbulent flame speed expression is implemented into Ford in-house CFD code MESIM (multi-dimensional engine simulation), and is validated extensively.

A KIVA-based CFD tool has been utilized to simulate the effect of a Common-Rail injection system applied to a large, uniflow-scavenged, two-stroke diesel engine. In particular, predictions for variations of injection pressure and injection duration have been validated with experimental data. The computational models have been evaluated according to their predictive capabilities of the combustion behavior reflected by the pressure and heat release rate history and the effects on nitric oxide formation and wall temperature trends. In general, the predicted trends are in good agreement with the experimental observations, thus demonstrating the potential of CFD as a design tool for the development of large diesel engines equipped with Common-Rail injection. Existing deficiencies are identified and can be explained in terms of model limitations, specifically with respect to the description of turbulence and combustion chemistry.

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.

Well-mixed lean SI engine operation can provide improvements of the fuel economy relative to that of traditional well-mixed stoichiometric SI operation. This work examines the use of two methods for improving the stability of lean operation, namely multi-pulse transient plasma ignition and intake air preheating. These two methods are compared to standard SI operation using a conventional high-energy inductive ignition system without intake air preheating. E85 is the fuel chosen for this study. The multi-pulse transient plasma ignition system utilizes custom electronics to generate 10 kHz bursts of 10 ultra-short (12ns), high-amplitude pulses (200 A). These pulses were applied to a custom spark plug with a semi-open ignition cavity. High-speed imaging reveals that ignition in this cavity generates a turbulent jet-like early flame spread that speeds up the transition from ignition to the main combustion event.

A recently developed spark plug equipped with fiber-optic flame-arrival detectors has been used to measure the motion and rate of growth of the early flame kernel. The cylinder pressure and gas velocity in the spark gap were measured simultaneously with the flame kernel measurements, permitting the data to be analyzed on a cycle-by-cycle basis to identify cause-and-effect correlations between the measured parameters. The data were obtained in a homogeneous-charge research engine that could be modified to produce three very different flow fields: (1) high swirl with high turbulence intensity, (2) tumble vortex with moderate turbulence intensity, and (3) negligible bulk motion with low turbulence intensity. The results presented show a moderate correlation between the combustion duration and the rate of growth of the flame kernel, but virtually no correlation with either the magnitude or direction of movement of the flame kernel away from the spark gap.

Recently, high-speed optical imaging data for a single operating point of a spray-guided gasoline engine has, along with the flamelet model and the G-equation theory, enabled the development of the new spark-ignition model SparkCIMM. Within its framework, detailed chemistry flamelet models capture the experimental feature of multiple localized ignition events along the excessively stretched and restriking spark channel, as well as the observations of non-spherical highly corrugated early turbulent flame fronts. The developed flamelet models account for the substantial turbulent fluctuations in equivalence ratio and enthalpy present under spray-guided conditions. A non-unity Lewis number formulation captures the deficient species diffusion into the highly curved flame reaction zone.

A one-dimensional model for crevice HC post-flame oxidation is used to calculate and understand the effect of operating parameters and fuel type (propane and isooctane) on the extent of crevice hydrocarbon and the product distribution in the post flame environment. The calculations show that the main parameters controlling oxidation are: bulk burned gas temperatures, wall temperatures, turbulent diffusivity, and fuel oxidation rates. Calculated extents of oxidation agree well with experimental values, and the sensitivities to operating conditions (wall temperatures, equivalence ratio, fuel type) are reasonably well captured. Whereas the bulk gas temperatures largely determine the extent of oxidation, the hydrocarbon product distribution is not very much affected by the burned gas temperatures, but mostly by diffusion rates. Uncertainties in both turbulent diffusion rates as well as in mechanisms are an important factor limiting the predictive capabilities of the model.

Formaldehyde (HCHO) near the lift-off location in a reacting diesel jet was visualized using planar laser-induced fluorescence (PLIF). Simultaneous imaging of OH chemiluminescence identified the high-temperature combustion region (lift-off). Experiments were performed in a constant-volume combustion vessel at ambient gas conditions (temperature and oxygen concentration) that generate no-soot, low-soot and moderate-soot diesel jets during mixing-controlled combustion. For no-soot conditions, results show that HCHO is formed upstream of the lift-off location and is consumed downstream of the lift-off length in fuel-rich premixed reaction zones at the jet center. Despite the fuel-rich combustion, and downstream regions that are surrounded by a high-temperature diffusion flame, there is no detectable PAH formation in the no-soot condition.

High resolution planar laser-induced fluorescence (PLIF) measurements were performed in an optically accessible internal combustion (IC) engine to investigate the behavior of scalar dissipation and the fine-scale structures of the turbulent scalar field. The fluorescent tracer fluorobenzene was doped into one of the two intake streams and nitrogen was used as the carrier gas to permit high signal-to-noise ratio fluorescence measurements without oxygen quenching effects. The resulting two-dimensional images allowed for an analysis of the structural detail of the scalar and scalar dissipation fields defined by the mixing of the two adjacent intake streams. High levels of scalar dissipation were found to be located within convoluted, sheet-like structures in accordance with previous studies. The fluorescence data, which were acquired during the intake stroke, were also used to examine the scalar energy and dissipation spectra.

Diesel combustion and emissions formation is largely spray and mixing controlled and hence understanding spray parameters, specifically vaporization, is key to determine the impact of fuel injector operation and nozzle design on combustion and emissions. In this study, an eight-hole common rail piezoelectric injector was tested in an optically accessible constant volume combustion vessel at charge gas conditions typical of full load boosted engine operation. Liquid penetration of the eight sprays was determined via processing of images acquired from Mie back scattering under vaporizing conditions by injecting into a charge gas at elevated temperature with 0% oxygen. Conditions investigated included a charge temperature sweep of 800 to 1300 K and injection pressure sweep of 1034 to 2000 bar at a constant charge density of 34.8 kg/m₃.

Combustion cycle-to-cycle variation (CCV) of Spark-Ignition (SI) engines can be influenced by the cyclic variations in charge motion, trapped mass and mixture composition inside the cylinder. A high CCV leads to misfire or knock, limiting the engine’s operating regime. To understand the mechanism of the effect of flow field and mixture compositions on CCV, the present numerical work was performed in a single cylinder Direct Injection Spark-Ignition (DISI) engine. A large eddy simulation (LES) approach coupled with the G-equation combustion model was developed to capture the CCV by accurately resolving the turbulent flow field spatially and temporally. Further, the ignition process was modeled by sourcing energy during the breakdown and arc phases with a line-shape ignition model which could move with the local flow. Detailed chemistry was solved both inside and outside the flame front. A compact 48-species 152-reactions primary reference fuel (PRF) reduced mechanism was used.

A novel algebraic similarity model for subgrid scalar dissipation rate has been developed as part of the Large Eddy Simulation (LES) package KIVA3V-LES for diesel engine study. The model is proposed from an a priori study using Direct Numerical Simulation (DNS) of forced isotropic turbulence. In the a posteriori test, fully resolved turbulent passive scalar field measurements are used to validate the model in actual engine flows. For reason of the length limit by SAE and the specific interest in engine applications, only a prior test and a posteriori test in engine flows are included in this paper. A posteriori tests in isotropic cube flow, turbulent round jet and flame cases will be presented in separate papers. An engine LES simulation of multi consecutive cycles was performed in this study.

Recent experiments of diesel injection processes have demonstrated that mixing accelerates after the end of injection (EOI). This finding has significant implications for low-temperature combustion (LTC) diesel engines. Previous simulations using a one-dimensional model of a single-pulsed air jet, which is analogous in many aspects to diesel jets, suggest that the rapid mixing after EOI in diesel jets is due to a temporary increase in the entrainment rate as the jet decelerates. In the present study, we performed a high-fidelity large eddy simulation (LES) of an unsteady air jet identical to that used for the one-dimensional model. The LES calculation agrees well with available experimental data and provides both spatially and temporally resolved details of the three-dimensional transient mixing field. Results show that entrainment increases during deceleration.

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.

Objective of this work is the incorporation of the flame stretch effects in an Eulerian-Lagrangian model for premixed SI combustion in order to describe ignition and flame propagation under highly inhomogeneous flow conditions. To this end, effects of energy transfer from electrical circuit and turbulent flame propagation were fully decoupled. The first ones are taken into account by Lagrangian particles whose main purpose is to generate an initial burned field in the computational domain. Turbulent flame development is instead considered only in the Eulerian gas phase for a better description of the local flow effects. To improve the model predictive capabilities, flame stretch effects were introduced in the turbulent combustion model by using formulations coming from the asymptotic theory and recently verified by means of DNS studies. Experiments carried out at Michigan Tech University in a pressurized, constant-volume vessel were used to validate the proposed approach.

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.

This paper presents a study of the turbulence field in an optical diesel engine operated under motored conditions using both large eddy simulation (LES) and Particle Image Velocimetry (PIV). The study was performed in a laboratory optical diesel engine based on a recent production engine from VOLVO Car. PIV is used to study the flow field in the cylinder, particularly inside the piston bowl that is also optical accessible. LES is used to investigate in detail the structure of the turbulence, the vortex cores, and the temperature field in the entire engine, all within a single engine cycle. The LES results are compared with the PIV measurements in a 40 × 28 mm domain ranging from the nozzle tip to the cylinder wall. The LES grid consists of 1283 cells. The grid dynamically adjusts itself as the piston moves in the cylinder so that the engine cylinder, including the piston bowl, is described by the grid.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.