Search Results

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.

Three time integration schemes and four finite volume interpolation schemes for the convection term in momentum equation were tested under turbulent planar gas jet and Sandia non-reacting vaporizing Spray-H cases. The three time integration schemes are the first-order Euler implicit scheme, the second-order backward scheme, and the second-order Crank-Nicolson scheme. The four spatial interpolation schemes are cubic central, linear central, upwind, and vanLeer schemes. Velocity magnitude contour, centerline and radial mean velocity and Reynolds stress profiles for the planar turbulent gas jet case, and fuel vapor contour and liquid and vapor penetrations for the Diesel spray case predicted by the different numerical schemes were compared. The sensitivity of the numerical schemes to mesh resolution was also investigated. The non-viscosity based dynamic structure subgrid model was used. The numerical tool used in this study was OpenFOAM.

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.

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.

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.

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.

Using non-viscosity dynamic structure Large Eddy Simulations (LES) turbulence model, spray=induced turbulence is investigated on a number of different Computational Fluid Dynamics (CFD) grids of varying mesh sizes (from 0.5 to 2 mm mesh). Turbulent flow is induced inside a quiescent chamber by liquid fuel spray and then left to decay after end of injection by virtue of its molecular viscosity and turbulent dissipation. Coherent structures (CS) of this turbulent flow are constructed and visualized using λ2 definition. Using CS, analysis is performed on the turbulent flow around the liquid spray jet. These CS from LES are then compared against the results from RANS calculations as well. The visualization of CS helps to explain the mechanism of fuel-air mixing obtained from LES results and its difference with RANS calculations.

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

In-cylinder velocity measurements were acquired in a two-valve, single-cylinder research engine to study the bulk fluid motion and small-scale turbulence. Different port geometries (two), different port orientations (two) and both shrouded and non-shrouded intake valves were tested to vary the intake-generated flow. Tests were performed at engine speeds of 300, 600, 900 and 1200 RPM with an atmospheric intake pressure. Prior to testing on the engine, the different head configurations were tested on a steady flow bench. Particle image velocimetry data were taken on a single plane, parallel to the piston surface, in the engine using both a low magnification to characterize the large-scale flow phenomena, and a high magnification to characterize the turbulence field. The low-magnification results showed that the swirl center location was relatively insensitive to engine speed, but did change position throughout the cycle.

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.

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.

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.

For the application of advanced clean combustion technologies, such as diesel HCCI/LTC, a compressor with high efficiency over a broad operation range is required to supply a high amount of EGR with minimum pumping loss. A compressor with high pitch of vaneless diffuser would substantially improve the flow range of the compressor, but it is at the cost of compressor efficiency, especially at low mass flow area where most of the city driving cycles resides. In present study, an ultra low solidity compressor vane diffuser was numerically investigated. It is well known that the flow leaving the impeller is highly distorted, unsteady and turbulent, especially at relative low mass flow rate and near the shroud side of the compressor. A conventional vaned diffuser with high stagger angle could help to improve the performance of the compressor at low end. However, adding diffuser vane to a compressor typically restricts the flow range at high end.

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.

Prior HCCI optical engine experiments utilizing laser-induced fluorescence (LIF) measurements of stratified fuel-air mixtures have demonstrated the utility of probability density function (PDF) statistics for correlating mixture preparation with combustion. However, PDF statistics neglect all spatial details of in-cylinder fuel distribution. The current computational paper examines the effects of spatial fuel distribution on combustion using a novel combination of a 3-D CFD model with a 1-D linear-eddy model of turbulent mixing. In the simulations, the spatial coarseness of initial fuel distribution prior to the start of heat release is varied while keeping PDF statistics constant. Several cases are run, and as the initial mixture is made coarser, combustion phasing monotonically advances due to high local equivalence ratios that persist longer. The effect of turbulent mixing is more complex.

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.

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.

The reaction zone of a diesel fuel jet stabilizes at a location downstream of the fuel injector once the initial autoignition phase is over. This distance is referred to as flame lift-off length. Recent investigations have examined the effects of a wide range of parameters (injection pressure, orifice diameter, and ambient gas temperature, density and oxygen concentration) on lift-off length under quiescent diesel conditions. Many of the experimental trends in lift-off length were in agreement with scaling laws developed for turbulent, premixed flame propagation in gas-jet lifted flames at atmospheric conditions. However, several effects did not correlate with the gas-jet scaling laws, suggesting that other mechanisms could be important to lift-off stabilization at diesel conditions. This paper shows experimental evidence that ignition processes affect diesel lift-off stabilization.

Wall-wetting due to liquid fuel film motion and fuel droplet impingement on combustion chamber walls is a major source of unburned hydrocarbons (UBHC), and is a concern for oil dilution in PFI engines. An experimental study was carried out to investigate the effects of injection timing, a charge motion control device, and the matching of injector with port geometry, on the “footprints” of liquid fuel inside the combustion chamber during the PFI engine starting process. Using a gasoline-soluble dye and filter paper deployed on the cylinder liner and piston top land surfaces to capture the liquid fuel footprints, the effects of the mixture formation processes on the wetted footprints can be qualitatively and quantitatively examined by comparing the wetted footprint locations and their color intensities. Real-time filming of the development of wetted footprints using a high-speed camera can also show the time history of the fuel wetting process inside an optically accessible engine.