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.

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.

Partially premixed combustion (PPC) is an operating mode that lies between the conventional compression ignition (CI) mode and homogeneous charge compression ignition (HCCI) mode. The combustion in this mixed mode is complex as it is neither diffusion-controlled (CI mode) nor governed solely by chemical kinetics (HCCI mode). In this study, CFD simulations were performed to evaluate flame index, which distinguishes between zones having a premixed flame and non-premixed flame. Experiments performed in the optical engine supplied data to validate the model. In order to realize PPC, the start of injection (SOI) was fixed at −40 CAD (aTDC) so that a required ignition delay is created to premix air/fuel mixture. The reference operating point was selected to be with 3 bar IMEP and 1200 rpm. Naphtha with a RON of 77 and its corresponding PRF surrogate were tested. The simulations captured the general trends observed in the experiments well.

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.

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.

Investigating combustion characteristics of oxygenated gasoline and gasoline blended ethanol is a subject of recent interest. The non-linearity in the interaction of fuel components in the oxygenated gasoline can be studied by developing chemical kinetics of relevant surrogate of fewer components. This work proposes a new reduced four-component (isooctane, heptane, toluene, and ethanol) oxygenated gasoline surrogate mechanism consisting of 67 species and 325 reactions, applicable for dynamic CFD applications in engine combustion and sprays. The model introduces the addition of eight C1-C3 species into the previous model (Li et al; 2019) followed by extensive tuning of reaction rate constants of C7 - C8 chemistry. The current mechanism delivers excellent prediction capabilities in comprehensive combustion applications with an improved performance in lean conditions.

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.

Stochastic pre-ignition remains one of the major barriers limiting further engine downsizing and down-speeding; two widely used strategies for improving the efficiency of spark-ignited engines. One of the most cited mechanisms thought to be responsible for pre-ignition is the ignition of a rogue droplet composed of lubricant oil and fuel. This originates during mixture formation from interactions between the fuel spray and oil on the cylinder liner. In the present study, this hypothesis is further examined using a single cylinder supercharged engine which employs a range of air-fuel mixture formation strategies. These strategies include port-fuel injection (PFI) along with side and central direct injection (DI) of an E5 gasoline (RON 97.5) using single and multiple injection events. Computational fluid dynamic (CFD) calculations are then used to explain the observed trends.

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.

The experimental in-cylinder combustion process was compared with the numerical simualtion for naphtha fuel under conventional compression ignition (CI) and partially premixed combustion (PPC) conditions. The start of injection timing (SOI) with the single injection strategy was changed from late of −10 CAD aTDC to early of −40 CAD aTDC. The three-dimensional full cycle engine combustion simulation was performed coupling with gas phase chemical kinetics by the CFD code CONVERGE™. The flame index was used for evaluating the combustion evolution of premixed flame and diffusion flame. The results show that the flame index could be used as an indicator for in-cylinder homogeneity evaluation. Hydroperoxyl shows a similar distribution with the premixed combustion. Formaldehyde could be used as an indicator for low temperature combustion.

In modern compression ignition engines, the dense liquid fuel is directly injected into high pressure and temperature atmosphere, so the spray transitions from subcritical to supercritical conditions. To gain better control of the spray-combustion heat release process, it is important to have a physically accurate description of the spray development process. This work explored the effect of real-fluid thermodynamics in the computational prediction of multiphase flow for two non-ideal situations: the cryogenic nitrogen and non-cryogenic n-dodecane and ammonia sprays. Three real-fluid equations of state (EoS) such as the Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), and Redlich-Kwong-Peng-Robinson (RKPR) coupled with the real-fluid Chung transport model were implemented in OpenFoam to predict the real-fluid thermodynamic properties. Validations against the CoolProp database were conducted.

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 study showcases the strength of machine learning (ML) models in imitating the operation of an advanced engine concept - the gasoline compression ignition (GCI) - at low loads. The GCI engine is prone to exceeding the limits of criteria emissions at such loads, especially at the cold start when the catalyst is not activated. One proposition to accelerate catalyst light-off is using spark-ignition. This, however, adds an extra level of complexity in identifying an optimum operation point. The ML models can be a useful tool in guiding the engine calibration process. In this study, the ML models are trained on GCI engine experiments, covering different intake conditions, injection strategies, and spark settings. The models can predict seven engine performance parameters: fuel consumption, four engine-out emissions, exhaust temperature, and coefficient of variation (COV) in indicated mean effective pressure (IMEP).

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.

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.