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.

The importance of radiative heat transfer on the combustion and soot formation characteristics under nominal ECN Spray A conditions has been studied numerically. The liquid n-dodecane fuel is injected with 1500 bar fuel pressure into the constant volume chamber at different ambient conditions. Radiation from both gas-phase as well as soot particles has been included and assumed as gray. Three different solvers for the radiative transfer equation have been employed: the discrete ordinate method, the spherical-harmonics method and the optically thin assumption. The radiation models have been coupled with the transported probability density function method for turbulent reactive flows and soot, where unresolved turbulent fluctuations in temperature and composition are included and therefore capturing turbulence-chemistry-soot-radiation interactions. Results show that the gas-phase (mostly CO2 ad H2O species) has a higher contribution to the net radiation heat transfer compared to soot.

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.

Fuel injection characteristics, in particular the atomization and penetration of the fuel droplets in the region close to the nozzle orifice, are known to affect emission and particulate formation in Diesel engines. It is also well established that the primary fuel atomization process is induced by aerodynamics in the near nozzle region as well as cavitation and turbulence from the injector nozzle. Typical breakup models in the literature however, do not consider the effects of cavitation and turbulence from nozzle injector. In this paper, a comprehensive primary breakup model incorporating the inner nozzle flow effects such as cavitation and turbulence along with aerodynamically induced breakup is developed and incorporated in the CONVERGE CFD code. This new primary breakup model is tested in a constant volume spray chamber against various spray data available in the literature.

By taking advantage of high-intensity and high-brilliance x-ray beams available at the Advanced Photon Source (APS), ultrafast (150 ps) propagation-based phase-enhanced imaging was developed to visualize high-pressure high-speed diesel sprays in the optically dense near-nozzle region. The sub-ns temporal and μm spatial resolution allows us to capture the morphology of the high-speed fuel sprays traveling at 500 m/s with a negligible motion blur. Both quality and quantitative information about the spray feature can be readily obtained. In the experiment, two types of single-hole nozzles have been used, one with a hydroground orifice inlet and the other with a sharp one. Within 3 mm from the nozzle, the sprays from these nozzles behave differently, ranging from laminar flow with surface instability waves to turbulent flow. The sprays are correlated with the nozzle internal geometry, which provides practical information for both nozzle design and supporting numerical simulation models.

Three-dimensional diesel engine combustion simulations with single-step chemistry have been compared with two-step and three-step chemistry by means of the Laminar and Turbulent Characteristic Time Combustion model using the Star-CD program. The second reaction describes the oxidation of CO and the third reaction describes the combustion of H2. The comparisons have been performed for two heavy-duty diesel engines. The two-step chemistry was investigated for a purely kinetically controlled, for a mixing limited and for a combination of kinetically and mixing limited oxidation. For the latter case, two different descriptions of the laminar reaction rates were also tested. The best agreement with the experimental cylinder pressure has been achieved with the three-step mechanism but the differences with respect to the two-step and single-step reactions were small.

A modified version of the Laminar and Turbulent Characteristic Time combustion model and the Hiroyasu-Magnussen soot model have been implemented in the flow solver Star-CD. Combustion simulations of three DI diesel engines, utilizing the standard k-ε turbulence model and a modified version of the RNG k-ε turbulence model, have been performed and evaluated with respect to combustion performance and emissions. Adjustments of the turbulent characteristic combustion time coefficient, which were necessary to match the experimental cylinder peak pressures of the different engines, have been justified in terms of non-equilibrium turbulence considerations. The results confirm the existence of a correlation between the integral length scale and the turbulent time scale. This correlation can be used to predict the combustion time scale in different engines.

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.

The fuel spray behavior in the near nozzle region of a gasoline injector is challenging to predict due to existing pressure gradients and turbulences of the internal flow and in-nozzle cavitation. Therefore, statistical parameters for spray characterization through experiments must be considered. The characterization of spray velocity fields in the near-nozzle region is of particular importance as the velocity information is crucial in understanding the hydrodynamic processes which take place further downstream during fuel atomization and mixture formation. This knowledge is needed in order to optimize injector nozzles for future requirements. In this study, the results of three experimental approaches for determination of spray velocity in the near-nozzle region are presented. Two different injector nozzle types were measured through high-speed shadowgraph imaging, Laser Doppler Anemometry (LDA) and X-ray imaging.

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.

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.

Multi-hole DI injectors are being adopted in the advanced downsized DISI ICE powertrain in the automotive industry worldwide because of their robustness and cost-performance. Although their injector design and spray resembles those of DI diesel injectors, there are many basic but distinct differences due to different injection pressure and fuel properties, the sac design, lower L/D aspect ratios in the nozzle hole, closer spray-to-spray angle and hense interactions. This paper used Phase-Contrast X ray techniques to visualize the spray near a 3-hole DI gasoline research model injector exit and compared to the visible light visualization and the internal flow predictions using with multi-dimensional multi-phase CFD simulations. The results show that strong interactions of the vortex strings, cavitation, and turbulence in and near the nozzles make the multi-phase turbulent flow very complicated and dominate the near nozzle breakup mechanisms quite unlike those of diesel injections.

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

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.

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.

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.

Combustion is governed by only two phenomena: chemical reactions and mixing (i.e., transport of energy, species, and momentum). A Reynolds Averaged Navier-Stokes (RANS) turbulence model is commonly employed to account for the enhanced mixing due to the presence of turbulence in fluid flow. A RANS turbulence model enhances mixing by introducing a turbulent viscosity. The addition of a turbulent viscosity not only enhances mixing but it also eliminates smaller scales in the CFD simulation. Even though the turbulent viscosity eliminates smaller scales, it is common for RANS engine combustion simulations to be under-resolved. The lack of sufficient mesh resolution to resolve the remaining scales in a RANS combustion simulation may result in a significant sub-grid term that needs to be modeled. In the context of combustion simulation, it is shown that frequently this sub-grid term is significantly more important than Turbulent Chemistry Interaction terms (TCI).

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.