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.

The Reynolds stress turbulence model (RSM) has been developed to go beyond the Boussinesq hypothesis and to improve turbulence modeling of flows with significant mean streamline curvature and secondary flow. In this paper the RSM in commercial CFD software CONVERGE is tested for its performance and robustness when applying to complex flows. Several validation cases including flow over flat plate, vortex combustor, diesel engine spray and combustion were selected to test the RSM. The swirling flow in vortex combustor, non-reacting but vaporizing ECN Spray A (free jet) and Sandia small bore diesel engine case are used to demonstrate the benefits of the RSM over the widely used RNG k-epsilon model without model tuning. The vortex combustor case shows the RSM can provide good prediction for strong swirling flow. ECN spray A case was used to demonstrate that the RSM can accurately predict the liquid and vapor penetration lengths of a free jet under diesel engine conditions.

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 generalized re-normalization group (RNG) turbulence model based on the local "dimensionality" of the flow field is proposed. In this modeling approach the model coefficients C₁, C₂, and C₃ are all constructed as functions of flow strain rate. In order to further validate the proposed turbulence model, the generalized RNG closure model was applied to model the backward facing step flow (a classic test case for turbulence models). The results indicated that the modeling of C₂ in the generalized RNG closure model is reasonable, and furthermore, the predictions of the generalized RNG model were in better agreement with experimental data than the standard RNG turbulence model. As a second step, the performance of the generalized RNG closure was investigated for a complex engine flow.

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.

An integrated numerical model has been developed for diesel engine computations based on the KIVA-II code. The model incorporates a modified RNG k-ε, turbulence model, a ‘wave’ breakup spray model, the Shell ignition model, the laminar-and-turbulent characteristic-time combustion model, a crevice flow model, a spray/wall impingement model that includes rebounding and breaking-up drops, and other improved submodels in the KIVA code. The model was validated and applied to model successfully different types of diesel engines under various operating conditions. These engines include a Caterpillar engine with different injection pressures at different injection timings, a small Tacom engine at different loads, and a Cummins engine modified by Sandia for optical experiments. Good levels of agreement in cylinder pressures and heat release rate data were obtained using the same computer model for all engine cases.

A computational methodology has been developed for loss prediction in intake regions of internal combustion engines. The methodology consists of a hierarchy of four major tasks: (1) proper computational modeling of flow physics; (2) exact geometry and high quality and generation; (3) discretization schemes for low numerical viscosity; and (4) higher order turbulence modeling. Only when these four tasks are dealt with properly will a computational simulation yield consistently accurate results. This methodology, which is has been successfully tested and validated against benchmark quality data for a wide variety of complex 2-D and 3-D laminar and turbulent flow situations, is applied here to a loss prediction problem from industry. Total pressure losses in the intake region (inlet duct, manifold, plenum, ports, valves, and cylinder) of a Caterpillar diesel engine are predicted computationally and compared to experimental data.

Combustion simulations utilizing the characteristic time combustion model have been performed for four DI diesel engines ranging in size from heavy-duty to large-bore designs. It has been found that the pre-factor to the turbulent characteristic time acts as a scaling parameter between the engines. This phenomenon is explained in terms of the non-equilibrium behavior of the turbulent time and length scales, as is encountered in the rapidly distorting, spray-induced flows of DI diesel engines. In fact, the equilibrium assumption between turbulence production and dissipation, which forms the basis for the employed k-ε-type turbulence models, does not hold in these situations. For such flows, the real turbulent dissipation time scale is locally proportional to the turbulent characteristic time scale which is determined by a typical eddy turnover time.

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