Search Results

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.

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.

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

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.

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.

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.

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.

Shock waves have been recently observed in high-pressure diesel sprays. In this paper, three-dimensional numerical simulations of supersonic diesel spray injection have been performed to investigate the underlying dynamics of the induced shock waves and their interactions with the spray. A Volume-of-Fluid based method in the CFD software (CONVERGE) is used to model this multiphase phenomena. An adaptive Mesh Refinement (AMR) scheme is employed to capture the front of the spray and the shock waves with high fidelity. Simulation results are compared to the available experimental observations to validate the numerical procedure. Parametric studies with different injection and ambient conditions are conducted to examine the effect of these factors on the generation of shock waves and their dynamics.

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.

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

A state-of-the-art spray modeling methodology, recently applied to RANS simulations, is presented for LES calculations. Key features of the methodology, such as Adaptive Mesh Refinement (AMR), advanced liquid-gas momentum coupling, and improved distribution of the liquid phase, are described. The ability of this approach to use cell sizes much smaller than the nozzle diameter is demonstrated. Grid convergence of key parameters is verified for non-evaporating and evaporating spray cases using cell sizes down to 1/32 mm. It is shown that for global quantities such as spray penetration, comparing a single LES simulation to experimental data is reasonable, however for local quantities the average of many simulated injections is necessary. Grid settings are recommended that optimize the accuracy/runtime tradeoff for LES-based spray simulations.

A new Computational Fluid Dynamics (CFD) code has been developed in order to overcome the deficiencies of traditional grid generation and mesh motion methods. The new code uses a modified cut-cell Cartesian technique that eliminates the need for the computational grid to coincide with the geometry of interest. The code also includes state-of-the-art numerical techniques and sub-models for simulating the complex physical and chemical processes that occur in engines. Features such as shared and distributed memory parallelization, a multigrid pressure solver and user-specified grid embedding allow for efficient simulations while maintaining the grid resolution necessary for accurate engine modeling. In addition, a new Adaptive Grid Embedding (AGE) technique has been developed and implemented. Sub-models for turbulence, spray injection, spray breakup, liquid drop dynamics, ignition, combustion and emissions are also included in the code.

Recent measurements by Siebers et al. have shown that the flame of a high pressure Diesel spray stabilizes under quiescent conditions at a location downstream of the fuel injector. The effects of various ambient and injection parameters on the flame “lift-off” length have been investigated under typical Diesel conditions in a constant-volume combustion vessel. In the present study, the experiments of Siebers et al. have been modeled using a modified version of the KIVA-3V engine simulation code. Fuel injection and spray breakup are modeled using the KH-RT model that accounts for liquid surface instabilities due to the Kelvin-Helmholtz and Rayleigh-Taylor mechanisms. Combustion is simulated using Convergent Thinking's recently developed detailed transient chemistry solver (SAGE) that allows for any number of chemical species and reactions to be modeled.

In this study, a one-equation LES sub-grid model from Menon, et al. [5] is used in simulating the diesel combustion process. In addition, based on the one-equation methodology of Menon et al., a new one-equation LES scalar transport model is formulated. These models allow for the turbulent transfer coefficients for both momentum and scalar flux to be determined independent of each other. The turbulent viscosity, μt, is determined as a function of the sub-grid kinetic energy, which is in turn determined from the one-equation model. The formulation for the scalar transfer coefficient, μs, is similar to that of the turbulent viscosity, yet is made to be consistent with scalar transport. Results for the LES momentum transfer are compared to experimental data of a backward facing step. This model, in conjunction with the LES scalar flux model, is verified by comparing with experimental data for a non-reacting turbulent jet.

The present study extends the recently developed KIVA-GA computer code to incorporate a generalized piston bowl geometry parameterization and multi-mode optimization. The new code was used to optimize the combustion chamber geometry of a small-bore automotive Diesel engine. The Genetic Algorithm (GA) merit function, which was calculated with a modified version of the KIVA-3V engine simulation code, included NOx, unburned HC, soot, and fuel consumption. A novel parameterization was included in KIVA-GA that allows for a variable number of parameters to define the bowl shape. The in-house G-Smooth grid generation package was used to create the KIVA grids with a specified compression ratio and mesh resolution. The improved KIVA-GA methodology was used to optimize engine emissions and performance simultaneously for two operating conditions.