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.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

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

A computational fluid dynamics (CFD) guided combustion system optimization was conducted for a heavy-duty diesel engine running with a gasoline fuel that has a research octane number (RON) of 80. The goal was to optimize the gasoline compression ignition (GCI) combustion recipe (piston bowl geometry, injector spray pattern, in-cylinder swirl motion, and thermal boundary conditions) for improved fuel efficiency while maintaining engine-out NOx within a 1-1.5 g/kW-hr window. The numerical model was developed using the multi-dimensional CFD software CONVERGE. A two-stage design of experiments (DoE) approach was employed with the first stage focusing on the piston bowl shape optimization and the second addressing refinement of the combustion recipe. For optimizing the piston bowl geometry, a software tool, CAESES, was utilized to automatically perturb key bowl design parameters. This led to the generation of 256 combustion chamber designs evaluated at several engine operating conditions.

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.

In this study, the combustion system of a light-duty compression ignition engine running on a market gasoline fuel with Research Octane Number (RON) of 91 was optimized using computational fluid dynamics (CFD) and Machine Learning (ML). This work was focused on optimizing the piston bowl geometry at two compression ratios (CR) (17 and 18:1) and this exercise was carried out at full-load conditions (20 bar indicated mean effective pressure, IMEP). First, a limited manual piston design optimization was performed for CR 17:1, where a couple of pistons were designed and tested. Thereafter, a CFD design of experiments (DoE) optimization was performed where CAESES, a commercial software tool, was used to automatically perturb key bowl design parameters and CONVERGE software was utilized to perform the CFD simulations. At each compression ratio, 128 piston bowl designs were evaluated.

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.

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.

The accurate modeling of fuel spray behavior under diesel engine conditions requires well-characterized boundary conditions. Among those conditions, the spray cone angle is important due to its impact on the spray mixing process, flame lift-off locations and subsequent soot formation. The spray cone angle is a highly dynamic variable, but existing correlations have been developed mainly for diesel fuels at quasi-steady state and relatively low injection pressures. The objective of this study was to develop spray cone angle correlations for both diesel and a light-end gasoline fuel over a wide range of diesel-engine operating conditions that are capable of capturing both the transient and quasi-steady state processes. Two important macroscopic characteristics of solid cone sprays, the spray cone angle and spray penetration, were measured using a single-hole heavy-duty injector using two fuels at diesel engine conditions in an optical constant volume vessel.

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.

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.

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.

In internal combustion engine computational fluid dynamics (CFD) simulations, uncertainties arise from various sources, such as estimates of model parameters, experimental boundary conditions, estimates of chemical kinetic rates, etc. These uncertainties propagate through the model and may result in discrepancies compared to experimental measurements. The relative importance of the various sources of uncertainty can be quantified by performing a sensitivity analysis. In this work, global sensitivity analysis (GSA) was applied to engine CFD simulations of a low-temperature combustion concept called gasoline compression ignition, to understand the influence of experimental measurement uncertainties from various sources on specific targets of interest-spray penetration, ignition timing, combustion phasing, combustion duration, and emissions. The sensitivity of these targets was evaluated with respect to imposed uncertainties in experimental boundary conditions and fuel properties.

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.

Reynolds-averaged Navier-Stokes (RANS) turbulence model has been used extensively for diesel engine simulations due to its computational efficiency and is expected to remain the workhorse computational fluid dynamics (CFD) tool for industry in the near future. Alternatively, large eddy simulations (LES) can potentially deal with complex flows and cover a large disparity of turbulence length scales, which makes this technique more and more attractive in the engine community. An n-dodecane spray flame (Spray A from Engine Combustion Network) was simulated using a dynamic structure LES model to understand the transient behavior of this turbulent flame. The liquid spray was treated with a traditional Lagrangian method and the gas-phase reaction was closed using a delta probability density function (PDF) combustion model. A 103-species skeletal mechanism was used for n-dodecane chemical kinetic model.

The combustion system of a heavy-duty diesel engine operated in a gasoline compression ignition mode was optimized using a CFD-based response surface methodology and a machine learning genetic algorithm. One common dataset obtained from a CFD design of experiment campaign was used to construct response surfaces and train machine learning models. 128 designs were included in the campaign and were evaluated across three engine load conditions using the CONVERGE CFD solver. The design variables included piston bowl geometry, injector specifications, and swirl ratio, and the objective variables were fuel consumption, criteria emissions, and mechanical design constraints. In this study, the two approaches were extensively investigated and applied to a common dataset. The response surface-based approach utilized a combination of three modeling techniques to construct response surfaces to enhance the performance of predictions.