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.

In this work, a turbulent combustion model is developed for large eddy simulation (LES) using a novel flamelet tabulation technique based on the framework of the multi-flamelet representative interactive flamelet (RIF) model. The overall aim is to develop a detailed model with elaborate chemistry mechanisms, LES turbulence models and highly resolved grids leveraging the computational cost advantage of a tabulated model. A novel technique of implementing unsteady flamelet libraries by using the residence time instead of the progress variables is proposed. In this study, LES of n-dodecane spray flame is performed using the tabulated turbulent combustion model along with a dynamic structure subgrid model. A high-resolution mesh is employed with a cell size of 62.5 microns in the entire spray and combustion regions. This model is then validated against igniting n-dodecane sprays under diesel engine conditions.

This paper focuses on detailed numerical simulations of direct injection diesel and gasoline sprays from production grade, multi-hole injectors. In a dual-fuel engine the direct injection of both the fuels can facilitate appropriate mixture preparation prior to ignition and combustion. Diesel and gasoline sprays were simulated using high-fidelity Large Eddy Simulations (LES) with the dynamic structure sub-grid scale model. Numerical predictions of liquid penetration, fuel density distribution as well as transverse integrated mass (TIM) at different axial locations versus time were compared against x-ray radiography data obtained from Argonne National Laboratory. A necessary, but often overlooked, criterion of grid-convergence is ensured by using Adaptive Mesh Refinement (AMR) for both diesel and gasoline. Nine different realizations were performed and the effects of random seeds on spray behavior were investigated.

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.

The primary strength of large eddy simulation (LES) is in directly resolving the instantaneous large-scale flow features which can then be used to study critical flame properties such as ignition, extinction, flame propagation and lift-off. However, validation of the LES results with experimental or direct numerical simulation (DNS) datasets requires the determination of statistically-averaged quantities. This is typically done by performing multiple realizations of LES and performing a statistical averaging among this sample. In this study, LES of n-dodecane spray flame is performed using a well-mixed turbulent combustion model along with a dynamic structure subgrid model. A high-resolution mesh is employed with a cell size of 62.5 microns in the entire spray and combustion regions. The computational cost of each calculation was in the order of 3 weeks on 200 processors with a peak cell count of about 22 million at 1 ms.

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.

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.

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.

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.

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.

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.