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 modeling of fuel jet atomization is key in the characterization of Internal Combustion (IC) engines, and 3D Computational Fluid Dynamics (CFD) is a recognized tool to provide insights for design and control purposes. Multi-hole injectors with counter-bored nozzle are the standard for Gasoline Direct Injection (GDI) applications and the Spray-G injector from the Engine Combustion Network (ECN) is considered the reference for numerical studies, thanks to the availability of extensive experimental data. In this work, the behavior of the Spray-G injector is simulated in a constant volume chamber, ranging from sub-cooled (nominal G) to flashing conditions (G2), validating the models on Diffused Back Illumination and Phase Doppler Anemometry data collected in vaporizing inert conditions.

This paper presents a computational fluid dynamics (CFD) study of the Engine Combustion Network (ECN) Spray G under non-vaporizing condition, focusing on the impacts of fuel properties as well as realistic geometry on the atomization process. The large-eddy-simulation method, coupled with the volume-of-fluid method, is used to model the high-speed turbulent two-phase flow. A moving-needle boundary condition is applied to capture the internal flow boundary condition accurately. The injector geometry was measured with micron-level resolution using x-ray tomographic imaging at the Advanced Photon Source at Argonne National Laboratory, providing detailed machining tolerance and defects from manufacturing and a realistic rough surface. A 2.5-μm fine mesh is used to sufficiently resolve the details of liquid-gas interface and the breakup process.

Knock is a major bottleneck to achieving higher thermal efficiency in spark ignition (SI) engines. The overall tendency to knock is highly dependent on fuel anti-knock quality as well as engine operating conditions. It is, therefore, critical to gain a better understanding of fuel-engine interactions in order to develop robust knock mitigation strategies. In the present work, a numerical model based on three-dimensional (3-D) computational fluid dynamics (CFD) was developed to capture knock in a Cooperative Fuel Research (CFR) engine. For combustion modeling, a hybrid approach incorporating the G-equation model to track turbulent flame propagation, and a homogeneous reactor multi-zone model to predict end-gas auto-ignition ahead of the flame front and post-flame oxidation in the burned zone, was employed.

Modeling plume interaction and collapse for direct-injection gasoline sprays is important because of its impact on fuel-air mixing and engine performance. Nevertheless, the aerodynamic interaction between plumes and the complicated two-phase coupling of the evaporating spray has shown to be notoriously difficult to predict. With the availability of high-speed (100 kHz) Particle Image Velocimetry (PIV) experimental data, we compare velocity field predictions between plumes to observe the full temporal evolution leading up to plume merging and complete spray collapse. The target “Spray G” operating conditions of the Engine Combustion Network (ECN) is the focus of the work, including parametric variations in ambient gas temperature. We apply both LES and RANS spray models in different CFD platforms, outlining features of the spray that are most critical to model in order to predict the correct aerodynamics and fuel-air mixing.

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.

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.

It is well known that in-nozzle flow behavior can significantly influence the near-nozzle spray formation and mixing that in turn affect engine performance and emissions. This in-nozzle flow behavior can, in turn, be significantly influenced by fuel properties. The goal of this study is to characterize the behavior of two different fuels, namely, a straight-run naphtha that has an anti-knock index of 58 (denoted as “Full-Range Naphtha”) and n-dodecane, in a simulated multi-hole common-rail diesel fuel injector. Simulations were carried out using a fully compressible multi-phase flow representation based on the mixture model assumption with the Volume of Fluid method. Our previous studies have shown that the characteristics of internal and near-nozzle flow are strongly related to needle motion in both the along- and off-axis directions.

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.

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. Hence, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region.

In this last decade, non-destructive X-ray measurement techniques have provided unique insights into the internal surface and flow characteristics of automotive injectors. This has in turn contributed to enhancing the accuracy of Computational Fluid Dynamics (CFD) models of these critical injection system components. By employing realistic injector geometries in CFD simulations, designers and modelers have identified ways to modify the injectors’ design to improve their performance. In recent work, the authors investigated the occurrence of cavitation in a heavy-duty multi-hole diesel injector operating with a high-volatility gasoline-like fuel for gasoline compression ignition applications. They proposed a comprehensive numerical study in which the original diesel injector design would be modified with the goal of suppressing the in-nozzle cavitation that occurs when gasoline fuels are used.

For several years there has been a great deal of effort made in researching ways to run a compression ignition engine with simultaneously high efficiency and low emissions. Recently much of this focus has been dedicated to using gasoline-like fuels that are more volatile and less reactive than conventional diesel fuel to allow the combustion to be more premixed. One of the key challenges to using fuels with such properties in a compression ignition engine is stable engine operation at low loads. This paper provides an analysis of how stable gasoline compression ignition (GCI) engine operation was achieved down to idle speed and load on a multi-cylinder compression ignition engine using only 87 anti-knock index (AKI) gasoline. The variables explored to extend stable engine operation to idle included: uncooled exhaust gas recirculation (EGR), injection timing, injection pressure, and injector nozzle geometry.

The occurrence of cavitation inside injectors is generally undesirable since it can cause material erosion and result in deviations from the expected injector performance. Previous numerical work employing an injector geometry measured with x-ray diagnostics and operating with a high-volatility straight-run gasoline (SRG) has shown that: (1) most of the cavitation is generally observed at low needle lifts, (2) needle motion is responsible for asymmetric structures in the internal flow as well as large pressure and velocity gradients that trigger phase transition at the orifice inlets, and (3) cavitation affects the injector discharge coefficient and distribution of injected fuel. To explore the potential for material damage within the injector orifices due to cavitation cloud collapse, the cavitation-induced erosion risk assessment (CIERA) tool has been applied for the first time to the realistic geometry of a heavy-duty injector using the CONVERGE software.

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.