Search Results

Traditional Lagrangian spray modeling approaches for internal combustion engines are highly grid-dependent due to insufficient resolution in the near nozzle region. This is primarily because of inherent restrictions of volume fraction with the Lagrangian assumption together with high computational costs associated with small grid sizes. A state-of-the-art grid-convergent spray modeling approach was recently developed and implemented by Senecal et al., (ASME-ICEF2012-92043) in the CONVERGE software. The key features of the methodology include Adaptive Mesh Refinement (AMR), advanced liquid-gas momentum coupling, and improved distribution of the liquid phase, which enables use of cell sizes smaller than the nozzle diameter. This modeling approach was rigorously validated against non-evaporating, evaporating, and reacting data from the literature.

Nowadays Spark Ignition (SI) engine developments focus on downsizing, in order to increase the engine load level and consequently its efficiency. As a side effect, knock occurrence is strongly increased. The current strategy to avoid knock is to reduce the spark advance which limits the potential of downsizing in terms of consumption reduction. Reducing the engine propensity to knock is therefore a first order subject for car manufacturers. Engineers need competitive tools to tackle such a complex phenomenon. In this paper the 3D RANS simulations ability to satisfactorily represent knock tendencies is demonstrated. ECFM (Extended Coherent Flame Model) has been recently implemented by IFPEN in CONVERGE and coupled with TKI (Tabulated Kinetics Ignition) to represent Auto-Ignition in SI engine. These models have been applied on a single cylinder engine configuration dedicated to abnormal combustion study.

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

Using natural gas in an internal combustion engine (ICE) is emerging as a promising way to improve thermal efficiency and reduce exhaust emissions. In the development of such engine platforms, computational fluid dynamics (CFD) plays a fundamental role in the optimization of geometries and operating parameters. One of the most relevant issues in the simulation of direct injection (DI) gaseous processes is the accurate prediction of the gas jet evolution. The simulation of the injection process for a gaseous fuel does not require complex modeling, nevertheless properly describing high-pressure gas jets remains a challenging task. At the exit of the nozzle, the injected gas is under-expanded, the flow becomes supersonic and shocks occur due to compressibility effects. These phenomena lead to challenging computational requirements resulting from high grid resolution and low computational time-steps.

A computational fluid dynamics (CFD) guided combustion system optimization was conducted for a heavy-duty compression-ignition engine with a gasoline-like fuel that has an anti-knock index (AKI) of 58. The primary goal was to design an optimized combustion system utilizing the high volatility and low sooting tendency of the fuel for improved fuel efficiency with minimal hardware modifications to the engine. The CFD model predictions were first validated against experimental results generated using the stock engine hardware. A comprehensive design of experiments (DoE) study was performed at different operating conditions on a world-leading supercomputer, MIRA at Argonne National Laboratory, to accelerate the development of an optimized fuel-efficiency focused design while maintaining the engine-out NOx and soot emissions levels of the baseline production engine.

An extensive numerical study of two-phase flow inside the nozzle holes and the issuing jets for a multi-hole direct injection gasoline injector is presented. The injector geometry is representative of the Spray G nozzle, an eight-hole counter-bored injector, from the Engine Combustion Network (ECN). Homogeneous Relaxation Model (HRM) coupled with the mixture multiphase approach in the Eulerian framework has been utilized to capture the phase change phenomena inside the nozzle holes. Our previous studies have demonstrated that this approach is capable of capturing the effect of injection transients and thermodynamic conditions in the combustion chamber, by predicting phenomenon such as flash boiling. However, these simulations were expensive, especially if there is significant interest in predicting the spray behavior as well.

The current trend of downsizing used in gasoline engines, while reducing fuel consumption and CO2 emissions, imposes severe thermal loads inside the combustion chamber. These critical thermodynamic conditions lead to the possible auto-ignition (AI) of fresh gases hot-spots around Top-Dead-Center (TDC). At this very moment where the surface to volume ratio is high, wall heat transfer influences the temperature field inside the combustion chamber. The use of a realistic wall temperature distribution becomes important in the case of a downsized engine where fresh gases hot spots found near high temperature walls can initiate auto-ignition. This paper presents a comprehensive numerical methodology for an accurately prediction of thermodynamic conditions inside the combustion chamber based on Conjugate Heat Transfer (CHT).

Reynolds-averaged Navier-Stokes (RANS) modeling is expected to deliver an ensemble-averaged result for the majority of turbulent flows. This could lead to the conclusion that multi-cycle internal combustion engine (ICE) simulations performed using RANS must exhibit a converging numerical solution after a certain number of consecutive cycles. However, for some engine configurations unsteady RANS simulations are not guaranteed to deliver an ensemble-averaged result. In this paper it is shown that, when using RANS modeling to simulate multiple engine cycles, the cycle-to-cycle variations (CCV) generated from different initial conditions at each cycle are not damped out even after a large number of cycles. A single-cylinder GDI research engine is simulated using RANS modeling and the numerical results for 20 consecutive engine cycles are evaluated for two specific operating conditions.

A 500-hour test cycle has been used to evaluate the durability of a prototype high pressure common rail injection system operating up to 1800 bar with E10 gasoline. Some aspects of the original diesel based hardware design were optimized in order to accommodate an opposed-piston, two-stroke engine application and also to mitigate the impacts of exposure to gasoline. Overall system performance was maintained throughout testing as fueling rate and rail pressure targets were continuously achieved and no physical damage was observed in the low-pressure components. Injectors showed no deviation in their flow characteristics after exposure to gasoline and high resolution imaging of the nozzle spray holes and pilot valve assemblies did not indicate the presence of cavitation damage. The high pressure pump did not exhibit any performance degradation during gasoline testing and teardown analysis after 500 hours showed no evidence of cavitation erosion.

This paper presents on-the-fly chemical mechanism reduction (termed as dynamic mechanism reduction) for speeding up the chemistry solution for practical internal combustion engine simulations. Small mechanisms are built at each time step which are valid under the local conditions of each cell. The Directed Relation Graph with Error Propagation (DRGEP) algorithm is used for generating local skeletal mechanisms. Dynamic mechanism reduction is combined with adaptive zoning (termed as multi-zone) to achieve good computational speed-up for engine simulations. The accuracy and efficiency of dynamic mechanism reduction is evaluated for a wide range of scenarios including (a) Diesel combustion, (b) Homogeneous Charge Compression Ignition (HCCI) combustion, and (c) Dual fuel combustion.

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 order to improve understanding of the primary atomization process for diesel-like sprays, a collaborative experimental and computational study was focused on the near-nozzle spray structure for the Engine Combustion Network (ECN) Spray D single-hole injector. These results were presented at the 5th Workshop of the ECN in Detroit, Michigan. Application of x-ray diagnostics to the Spray D standard cold condition enabled quantification of distributions of mass, phase interfacial area, and droplet size in the near-nozzle region from 0.1 to 14 mm from the nozzle exit. Using these data, several modeling frameworks, from Lagrangian-Eulerian to Eulerian-Eulerian and from Reynolds-Averaged Navier-Stokes (RANS) to Direct Numerical Simulation (DNS), were assessed in their ability to capture and explain experimentally observed spray details. Due to its computational efficiency, the Lagrangian-Eulerian approach was able to provide spray predictions across a broad range of conditions.

A set of reduced chemical mechanisms was developed for use in multi-dimensional engine simulations of premixed gasoline combustion. The detailed Primary Reference Fuel (PRF) mechanism (1034 species, 4236 reactions) from Lawrence Livermore National Laboratory (LLNL) was employed as the starting mechanism. The detailed mechanism, referred to here as LLNL-PRF, was reduced using a technique known as Parallel Direct Relation Graph with Error Propagation and Sensitivity Analysis. This technique allows for efficient mechanism reduction by parallelizing the ignition delay calculations used in the reduction process. The reduction was performed for a temperature range of 800 to 1500 K and equivalence ratios of 0.5 to 1.5. The pressure range of interest was 0.75 bar to 40 bar, as dictated by the wide range in spark timing cylinder pressures for the various cases. In order to keep the mechanisms relatively small, two reductions were performed.

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.

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

A numerical study has been carried out to assess the effects of needle movement and internal nozzle flow on spray formation for a multi-hole Gasoline Direct Injection system. The coupling of nozzle flow and spray formation is dynamic in nature and simulations with pragmatic choice of spatial and temporal resolutions are needed to analyze the sprays in a GDI system. The dynamic coupling of nozzle flow and spray formation will be performed using an Eulerian-Lagrangian Spray Atomization (ELSA) approach. In this approach, the liquid fuel will remain in the Eulerian framework while exiting the nozzle, while, depending on local instantaneous liquid concentration in a given cell and amount of liquid in the neighboring cells, part of the liquid mass will be transferred to the Lagrangian framework in the form of Lagrangian parcels.

Using gasoline and diesel simultaneously in a dual-fuel combustion system has shown effective benefits in terms of both brake thermal efficiency and exhaust emissions. In this study, the dual-fuel approach is applied to a light-duty spark ignition (SI) gasoline direct injection (GDI) engine. Three combustion modes are proposed based on the engine load, diesel micro-pilot (DMP) combustion at high load, SI combustion at low load, and diesel assisted spark-ignition (DASI) combustion in the transition zone. Major focus is put on the DMP mode, where the diesel fuel acts as an enhancer for ignition and combustion of the mixture of gasoline, air, and recirculated exhaust gas. Computational fluid dynamics (CFD) is used to simulate the dual-fuel combustion with the final goal of supporting the comprehensive optimization of the main engine parameters.

This work involves modeling internal and near-nozzle flows of a gasoline direct injection (GDI) nozzle. The Engine Combustion Network (ECN) Spray G condition has been considered for these simulations using the nominal geometry of the Spray G injector. First, best practices for numerical simulation of the two-phase flow evolution inside and the near-nozzle regions of the Spray G injector are presented for the peak needle lift. The mass flow rate prediction for peak needle lift was in reasonable agreement with experimental data available in the ECN database. Liquid plume targeting angle and liquid penetration estimates showed promising agreement with experimental observations. The capability to assess the influence of different thermodynamic conditions on the two-phase flow nature was established by predicting non-flashing and flashing phenomena.

Understanding the short-lived structure of the plasma that forms between the electrodes of a spark plug is crucial to the development of improved ignition models for SI engines. However, measuring the amount of energy deposited in the gas directly and non-intrusively is difficult, due to the short time scales and small length scales involved. The breakdown of the spark gap occurs at nanosecond time scales, followed by an arc phase lasting a few microseconds. Finally, a glow discharge phase occurs over several milliseconds. It is during the arc and glow discharge phases that most of the heat transfer from the plasma to the electrodes and combustion gases occurs. Light emission can be used to measure an average temperature, but micron spatial resolution is required to make localized measurements.