Search

Numerical predictions of the fuel heating and cavitation erosion location indicators occurring during the opening and closing periods of the needle valve inside a five-hole common rail Diesel fuel injector are presented. These have been obtained using an explicit density-based solver of the compressible Navier-Stokes (NS) and energy conservation equations; the flow solver is combined with two thermodynamic closure models for the liquid, vapour and vapour-liquid equilibrium (VLE) property variation as function of pressure and temperature. The first is based on tabulated data for a 4-component Diesel fuel surrogate, derived from the Perturbed-Chain, Statistical Associating Fluid Theory (PC-SAFT) Equation of State (EoS), allowing for the variation of the physical and transport properties of the fuel with the local pressure and temperature to be quantified. The second thermodynamic closure is based on the widely used barotropic Equation of State (EoS) approximation between density and pressure only and neglects viscous heating. The Wall Adapting Local Eddy viscosity (WALE) LES model was used to resolve sub-grid scale turbulence while a cell-based mesh deformation Arbitrary Lagrangian–Eulerian (ALE) formulation is used for modelling the injector's needle valve movement. Model predictions are found in close agreement against 0-D estimates of the temporal variation of the fuel temperature difference between the feed and hole exit during the injection period. Two mechanisms affecting the temperature distribution within the fuel injector have been revealed and quantified. The first is ought to wall friction-induced heating, which may result to local liquid temperature increase up to fuel's boiling point while superheated vapour is formed. At the same time, liquid expansion due to the depressurisation of the injected fuel results to liquid cooling relative to the fuel's feed temperature; this is occurring at the central part of the injection orifice. The spatial and temporal temperature and pressure gradients induce significant variations in the fuel density and viscosity, which in turn, affect the formed coherent vortical flow structures. It is found, in particular, that these affect the locations of cavitation formation and collapse, that may lead to erosion of the surfaces of the needle valve, sac volume and injection holes. Model predictions are compared against corresponding X-ray surface erosion images obtained from injector durability tests, showing good agreement.

An investigation of the fuel heating, vapor formation, and cavitation erosion location patterns inside a five-hole common rail diesel fuel injector, occurring during the early opening period of the needle valve (from 2 μm to 80 μm), discharging at pressures of up to 450 MPa, is presented. Numerical simulations were performed using the explicit density-based solver of the compressible Navier–Stokes (NS) and energy conservation equations. The flow solver was combined with tabulated property data for a four-component diesel fuel surrogate, derived from the perturbed chain statistical associating fluid theory (PC-SAFT) equation of state (EoS), which allowed for a significant amount of the fuel’s physical and transport properties to be quantified. The Wall Adapting Local Eddy viscosity (WALE) Large Eddy Simulation (LES) model was used to resolve sub-grid scale turbulence, while a cell-based mesh deformation arbitrary Lagrangian–Eulerian (ALE) formulation was used for modelling the injector’s needle valve movement. Friction-induced heating was found to increase significantly when decreasing the pressure. At the same time, the Joule–Thomson cooling effect was calculated for up to 25 degrees K for the local fuel temperature drop relative to the fuel’s feed temperature. The extreme injection pressures induced fuel jet velocities in the order of 1100 m/s, affecting the formation of coherent vortical flow structures into the nozzle’s sac volume.

In the present numerical study, implicit large eddy simulations (LES) of non-reacting multi-components mixing processes of cryogenic nitrogen and n-dodecane fuel injections under transcritical and supercritical conditions are carried out, using a modified reacting flow solver, originally available in the open source software OpenFOAM®. To this end, the Peng-Robinson (PR) cubic equation of state (EOS) is considered and the solver is modified to account for the real-fluid thermodynamics. At high pressure conditions, the variable transport properties such as dynamic viscosity and thermal conductivity are accurately computed using the Chung transport model. To deal with the multicomponent species mixing, molar averaged homogeneous classical mixing rules are used. For the velocity-pressure coupling, a PIMPLE based compressible algorithm is employed. For both cryogenic and non-cryogenic fuel injections, qualitative and quantitative analyses are performed, and the results show significant effects of the chamber pressure on the mixing processes and entrainment rates. The capability of the proposed numerical model to handle multicomponent species mixing with real-fluid thermophysical properties is demonstrated, in both supercritical and transcritical regimes.

This paper introduces a new model for real gas thermodynamics, with improved accuracy, performance, and robustness compared to state-of-the-art models. It is motivated by the physical insight that in non-premixed flames, as encountered in high pressure liquid propellant rocket engines, mixing takes place chiefly in the hot reaction zone among ideal gases. We developed a new model taking advantage of this: When real fluid behavior only occurs in the cryogenic oxygen stream, this is the only place where a real gas equation of state (EOS) is required. All other species and the thermodynamic mixing can be treated as ideal. Real fluid properties of oxygen are stored in a library; the evaluation of the EOS is moved to a preprocessing step. Thus decoupling the EOS from the runtime performance, the method allows the application of accurate high quality EOS or tabulated data without runtime penalty. It provides fast and robust iteration even near the critical point and in the multiphase coexistence region. The model has been validated and successfully applied to the computation of 0D phase change with heat addition, and a supercritical reactive coaxial LOX/GH 2 single injector.