Search Results

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.

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

The design and testing of the valve train for a new two-stroke diesel engine concept [1,2] is presented. The gas exchange of this process requires extremely fast-acting inlet valves, which constituted a very demanding designing task. A simulation model of the prototype valve train was constructed with commercially available software. The simulation program served as the main tool for optimizing the dynamic behavior of the valve train. The prototype valve train was built according to the simulations and valve acceleration measurements were performed in order to validate the simulation results. The simulations and measurements are presented in detail in this paper.

A novel two-stroke engine concept is introduced. The cylinder scavenging takes place during the upward motion of the piston. The gas exchange valves are similar to typical four-stroke valves, but the intake valves are smaller and lighter. The scavenging air pressure is remarkably higher than in present-day engines. The high scavenging air pressure is produced by an external compressor. The two-stroke operation is achieved without the drawbacks of port scavenged engines. Moreover, the combustion circumstances, charge pressure and temperature and internal exhaust gas re-circulation (EGR) can be controlled by using valve timings. There is good potential for a substantial reduction in NOx emissions through the use of adjustable compression pressure and temperature and by using the adjustable amount of exhaust gas re-circulation.

Three-dimensional diesel engine combustion simulations with single-step chemistry have been compared with two-step and three-step chemistry by means of the Laminar and Turbulent Characteristic Time Combustion model using the Star-CD program. The second reaction describes the oxidation of CO and the third reaction describes the combustion of H2. The comparisons have been performed for two heavy-duty diesel engines. The two-step chemistry was investigated for a purely kinetically controlled, for a mixing limited and for a combination of kinetically and mixing limited oxidation. For the latter case, two different descriptions of the laminar reaction rates were also tested. The best agreement with the experimental cylinder pressure has been achieved with the three-step mechanism but the differences with respect to the two-step and single-step reactions were small.

A modified version of the Laminar and Turbulent Characteristic Time combustion model and the Hiroyasu-Magnussen soot model have been implemented in the flow solver Star-CD. Combustion simulations of three DI diesel engines, utilizing the standard k-ε turbulence model and a modified version of the RNG k-ε turbulence model, have been performed and evaluated with respect to combustion performance and emissions. Adjustments of the turbulent characteristic combustion time coefficient, which were necessary to match the experimental cylinder peak pressures of the different engines, have been justified in terms of non-equilibrium turbulence considerations. The results confirm the existence of a correlation between the integral length scale and the turbulent time scale. This correlation can be used to predict the combustion time scale in different engines.

Non-evaporating diesel sprays have been simulated utilizing the ETAB and the WAVE atomization and breakup models and have been compared with experimental data. The experimental penetrations and widths were determined from back-lit spray images and the droplet sizes have been measured by means of a Malvern particle sizer. The model evaluation criteria include the spray penetration, the spray width and the local droplet size. The comparisons have been performed for variations of the injection pressure, the gas density and the fuel viscosity. The fuel nozzle exit velocities used in the simulations have been computed with a special code that considers the effect of in-nozzle cavitation. The simulations showed good overall agreement with experimental data. However, the capabilities of the models to predict the droplet size for different fuels could be improved.

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.

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.

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

In this paper the KH-RT and the CAB droplet breakup models are analyzed. The focus is on near nozzle spray simulation data that will be qualitatively compared with results obtained from x-ray experiments. Furthermore, the suitability of the x-ray method for spray studies is assessed and its importance for droplet breakup modeling is discussed. The simulations have been carried out with the Kiva3VRel2 CFD-code into which the KH-RT- and the CAB- droplet breakup models have been implemented. Since the x-ray method gives an integrated line-of-sight mass distribution of the spray, a suitable comparison of the experimental distributions and the simulated ones is made. Additionally, modeling aspects are discussed and the functioning of the models demonstrated by illustrating how the parcel Weber numbers and radii vary spatially. The transient nature of the phenomenon is highlighted and the influence of the breakup model parameters is discussed.

Reduced-order models typically assume that the flow through the injector orifice is quasi-steady. The current study investigates to what extent this assumption is true and what factors may induce large-scale variations. Experimental data were collected from a single-hole metal injector with a smoothly converging hole and from a transparent facsimile. Gas, likely indicating cavitation, was observed in the nozzles. Surface roughness was a potential cause for the cavitation. Computations were employed using two engineering-level Computational Fluid Dynamics (CFD) codes that considered the possibility of cavitation. Neither computational model included these small surface features, and so did not predict internal cavitation. At steady state, it was found that initial conditions were of little consequence, even if they included bubbles within the sac. They however did modify the initial rate of injection by a few microseconds.

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