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.

Urea water solution-based Selective Catalytic Reduction (SCR) of NOx emissions from vehicular diesel engines is now widely used world-wide to meet strict health and environmental protection regulations. While urea-based SCR is proven effective, urea-derived deposits often form near injectors, on mixers and pipes, and on the SCR catalyst face. Further understanding of these deposit-formation processes is needed to design aftertreatment system hardware and control systems capable of avoiding severe urea-derived deposits. Computational Fluid Dynamics (CFD) is widely used in SCR aftertreatment design. Film formation, movement, solid wall cooling and deposit initiation/growth time-scales are in the range of minutes to hours, but traditional CFD simulations take too long to reach these time-scales. Here, we propose and demonstrate the frozen flow approach for pulsed sprays and conjugate heat transfer to reduce computation time while maintaining accuracy of key physics.

Detailed chemical kinetics, although preferred due to increased accuracy, can significantly slow down CFD combustion simulations. Chemistry solutions are typically the most computationally costly step in engine simulations. The calculation time can be significantly accelerated using a multi-zone combustion model. The multi-zone model is integrated into the CONVERGE CFD code. At each time-step, the CFD cells are grouped into zones based on the cell temperature and equivalence ratio. The chemistry solver is invoked only on each zone. The zonal temperature and mass fractions are remapped onto the CFD cells, such that the temperature and composition non-uniformities are preserved. Two remapping techniques published in the literature are compared for their relative performance. The accuracy and speed-up of the multi-zone model is improved by using variable bin sizes at different temperature and equivalence ratios.

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.

Ammonia is a promising carbon-free alternative fuel for use in combustion systems. The main associated challenges are its relatively low reactivity and high NOx emissions compared to conventional fuels. Therefore, the combustion behaviour of ammonia and ammonia blends still needs to be better understood over a wide range of conditions. To this end, a comprehensive chemical kinetic mechanism C3MechV3.4, which is an update of C3MechV3.3, has been developed for improved predictions of the combustion of ammonia and ammonia blends. C3MechV3.4 has been validated using a wide range of experimental results for pure ammonia and ammonia/hydrogen, ammonia/methanol and ammonia/n-heptane blends. These validations target different data sets including ignition delay times, species profiles measured as a function of time, and/or temperature and laminar flame speeds over a wide range of conditions.

Interest in operational cost reduction is driving engine manufacturers to consider low-cost fuel substitution in heavy-duty diesel engines. These dual-fuel (DF) engines could be operated either in diesel-only mode or operated with premixed natural gas (NG) ignited by a pilot flame of compression-ignited direct-injected diesel fuel. Under certain conditions, dual-fuel operation can result in increased cycle-to-cycle variability (CCV) during combustion. CFD can greatly help in understanding and identifying critical parameters influencing CCV. Innovative modelling techniques and large computing resources are needed to investigate the factors affecting CCV in dual-fuel engines. This paper discusses the use of the High Performance Computing resource Titan, at Oak Ridge National Laboratory, to investigate CCV of a dual-fuel engine.

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.

The Computational Chemistry Consortium (C3) is dedicated to leading the advancement of combustion and emissions modeling. The C3 cluster combines the expertise of different groups involved in combustion research aiming to refine existing chemistry models and to develop more efficient tools for the generation of surrogate and multi-fuel mechanisms, and suitable mechanisms for CFD applications. In addition to the development of more accurate kinetic models for different components of interest in real fuel surrogates and for pollutants formation (NOx, PAH, soot), the core activity of C3 is to develop a tool capable of merging high-fidelity kinetics from different partners, resulting in a high-fidelity model for a specific application. A core mechanism forms the basis of a gasoline surrogate model containing larger components including n-heptane, iso-octane, n-dodecane, toluene and other larger hydrocarbons.

The internal combustion engine (ICE) has long dominated the heavy-duty sector by using liquid fossil fuels such as diesel but global commitments by countries and OEMs to reduce lifecycle carbon dioxide (CO2) emissions has garnered interest in alternative fuels like hydrogen. Hydrogen is a unique gaseous fuel that contains zero carbon atoms and has desired thermodynamic properties of high energy density per unit mass and high flame speeds. However, there are challenges related to its adoption to the heavy-duty sector as a drop-in fuel replacement for compression ignition (CI) diesel combustion given its high autoignition resistance. To overcome this fundamental barrier, engine manufacturers are exploring dual fuel combustion engines by substituting a fraction of the diesel fuel with hydrogen which enables fuel flexibility when there is no infrastructure and retrofittability to existing platforms.

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

The standard capability of engine experimental studies is that ensemble averaged quantities like in-cylinder pressure from multiple cycles and emissions are reported and the cycle to cycle variation (CCV) of indicated mean effective pressure (IMEP) is captured from many consecutive combustion cycles for each test condition. However, obtaining 3D spatial distribution of all the relevant quantities such as fuel-air mixing, temperature, turbulence levels and emissions from such experiments is a challenging task. Computational Fluid Dynamics (CFD) simulations of engine flow and combustion can be used effectively to visualize such 3D spatial distributions. A dual fuel engine is considered in the current study, with manifold injected natural gas (NG) and direct injected diesel pilot for ignition.

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.

The internal structure of Diesel fuel injectors is known to have a significant impact on the nozzle flow and the resulting spray emerging from each hole. In this paper the three-dimensional transient flow structures inside a Diesel injector is studied under nominal (in-axis) and realistic (including off-axis lateral motion) operating conditions of the needle. Numerical simulations are performed in the commercial CFD code CONVERGE, using a two-phase flow representation based on a mixture model with Volume of Fluid (VOF) method. Moving boundaries are easily handled in the code, which uses a cut-cell Cartesian method for grid generation at run time. First, a grid sensitivity study has been performed and mesh requirements are discussed. Then the results of moving needle calculations are discussed. Realistic radial perturbations (wobbles) of the needle motion have been applied to analyze their impact on the nozzle flow characteristics.

Fluid-Structure Interaction (FSI) simulation approach can be used to simulate a turbocharger. However, this predictive 3D simulation encounters the challenge of a long computational time. The impeller speed can be above 100,000 rpm, and generally a CFD solver limits the maximum movement of the impeller surface per time step. The maximum movement must be a fraction (~0.3) of the cell length, thus the time step will be very small. A Multiple Reference Frame (MRF) approach can reduce computational time by eliminating the need to regenerate the mesh at each time-step to accommodate the moving geometry. A static local reference zone encompassing the impeller is created and the impact of the impeller movement is modeled via a momentum source. However, the MRF approach is not a predictive simulation because the impeller speed must be given by the User. A new simulation approach was introduced that coupled the FSI and MRF approach.

Gas turbines play an extremely important role in fulfilling a variety of power needs and are mainly used for power generation and propulsion applications. The performance and efficiency of gas turbine engines are to a large extent dependent on turbine rotor inlet temperatures: typically, the hotter the better. In gas turbines, the combustion temperature and the fuel efficiency are limited by the heat transfer properties of the turbine blades. However, in pushing the limits of hot gas temperatures while preventing the melting of blade components in high-pressure turbines, the use of effective cooling technologies is critical. Increasing the turbine inlet temperature also increases heat transferred to the turbine blade, and it is possible that the operating temperature could reach far above permissible metal temperature. In such cases, insufficient cooling of turbine blades results in excessive thermal stress on the blades causing premature blade failure.