Search Results

The basic idea behind large-eddy simulation (LES) is to accurately resolve the large energy-containing scales and to use subgrid-scale (SGS) models for the smaller scales. The accuracy of LES can be significantly impacted by the numerical discretization schemes and the choice of the SGS model. This work investigates the accuracy of low-order LES codes in the simulation of a turbulent round jet which is representative of fuel jets in engines. The turbulent jet studied is isothermal with a Reynolds number of 6800. It is simulated using Converge, which is second-order accurate in space and first-order in time, and FLEDS, developed at Purdue University, which is sixth-order accurate in space and fourth-order in time. The high-order code requires the resolution of acoustic time-scales and hence is approximately 10 times more expensive than the low-order code.

Accurate modeling of the transient structure of reacting diesel jets is important as transient features like autoignition, flame propagation, and flame stabilization have been shown to correlate with combustion efficiency and pollutant formation. In this work, results from Reynolds-averaged Navier-Stokes (RANS) simulations of flame lift-off in diesel jets are examined to provide insight into the lift-off physics. The large eddy simulation (LES) technique is also used to computationally model a lifted jet flame at conditions representative of those encountered in diesel engines. An unsteady flamelet progress variable (UFPV) model is used as the turbulent combustion model in both RANS simulations and LES. In the model, a look-up table of reaction source terms is generated as a function of mixture fraction Z, stoichiometric scalar dissipation rate Xst, and progress variable Cst by solving the unsteady flamelet equations.

In Diesel engines, the vapor phase of the fuel jet is known to impinge on the walls. This impingement is likely to have an effect on mixing characteristics, the structure of the diffusion flame and on pollutant formation and oxidation. These effects have not been studied in detail in the literature. In this work, the structure of a laminar wall jet that is generated from the impingement of a free laminar jet on a wall is discussed. We study the laminar jet with the belief that the local structure of the reaction zone in the turbulent reacting jet is that of a laminar flame. Results from non-reacting and reacting jets will be presented. In the case of the non-reacting jets, the focus of the inquiry is on assessing the accuracy of the computed results by comparing them with analytical results. Velocity profiles in the wall jet, growth rates of the half-width of the jet and penetration rates are presented.

In this paper, computations of pulsating flows in a duct with multiple inlets using the lattice Boltzmann method (LBM) are reported. As future emissions standards present a significant challenge for Diesel engine manufacturers, several options are being investigated to identify strategies to meet such regulations. Exhaust gas aftertreatment is one of the most important among them. As the performance of the various aftertreatment devices is sensitive to the flow conditions in the exhaust, a greater understanding of the flows under pulsating conditions in the presence of multiple cylinders is needed. The Lattice Boltzmann Method (LBM) is a relatively new and promising computational approach for applications to fluid dynamics problems. Two advantages of the method relative to traditional methods are ease of implementation and ease of parallelization and performance on parallel computers.

The spray submodel is an important component in multidimensional models for Diesel engines. The satisfactory representation of the spray is dependent on adequate representation of turbulence in the jet which, in part, determines its spread and penetration. In this work, the RNG k-ϵ model is evaluated relative to the standard k-ϵ model for computing turbulent jets. Computations are made for both gas jets and sprays. The gas jet is computed with an adequately high degree of numerical spatial resolution of the order of the orifice diameter. In the case of the spray, achieving such a high resolution would be challenging. Since the spray has similarities to the gas jet, and the gas jet may be computed with such high resolution and adequate accuracy, firm conclusions may be drawn for it and they may be applicable to sprays. It is concluded that the RNG k - ϵ model, in general, results in predictions of greater mixing in the jets relative to the standard model.

It is generally agreed that adequate resolution is required to reproduce the structure of spray and gas jets in numerical computations. It has not been clarified what this resolution should be although it would appear reasonable to assume that it should be such that the physical scales of the problem are resolved. In the case of a jet, this implies that near the orifice, the jet diameter has to be resolved since this is the appropriate length scale. It is shown in this work that if such a resolution is not used in computing transient jets, the structure of the jet is not reproduced with adequate accuracy. In fact, unexpected, erroneous and misleading dependence on ambient turbulence length and time scales will be predicted when the initial ambient turbulence diffusivity is small relative to the jet diffusivity. When the ambient turbulence diffusivity is of the same order as the jet diffusivity or greater, entrainment rates are significantly underpredicted.

The first comparisons of measured and computed mean velocities, and of measured fluctuation intensities and computed turbulence intensities in a motored rotary engine are presented. The computations were performed with a recently developed three–dimensional model. The measurements were made at the Sandia National Laboratories by Dimpelfeld and Witze at several locations along the rotor housing and at two engine speeds. The measured and computed mean velocities agree to within 15% whereas the computed turbulence intensities correctly are lower than the measured fluctuation intensities by about 30% to 50%, as anticipated. Under motored conditions, the turbulence intensity tends to be rather homogeneous and of similar magnitude somewhat before and after top dead center but significantly inhomogeneous and of greater magnitude around top dead center. The comparisons suggest that predictions of mean gas velocity and of turbulence intensity can be made with the available model.

Propagation of turbulent flames in spark-ignition engines is considered from the viewpoint of the different possible regimes of premixed turbulent combustion. Nondimensional parameters defining known combustion regimes are reviewed, and numerical values of these parameters are estimated for both research and production engines. The reaction-sheet regime is inferred to apply at least for some operating conditions, and therefore literature on turbulent flame propagation in the reaction-sheet regime is reviewed. Implications of these results on interpretations of existing experimental observations of combustion in engine cylinders and on modeling of turbulent flame propagation in engines are discussed.