Investigations on High-Speed Flash Boiling Atomization of Fuel Based on Numerical Simulations

Wei Zhong,Zhenfang Xin,Lihua Wang,⋆and Haiping Liu

1School of Aerospace Engineering and Applied Mechanics,Tongji University,Shanghai,200092,China

2Beijing Mechanical Equipment Research Institute,Beijing,100854,China

ABSTRACT Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises from the involvement of many complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study, an HFBA model for gasoline direct injection (GDI) is established.This model incorporates primary and secondary atomization, as well as vaporization and boiling models,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.

KEYWORDS High-speed flash boiling atomization; numerical simulations; Eulerian description; Lagrangian description;gasoline direct injection

1 Introduction

Over the last few decades,significant efforts have been made to increase the efficiency of internal combustion engines and simultaneously reduce tailpipe emissions.The spray atomization process in internal combustion engines directly affects the engine’s thermal efficiency and emissions [1].In gasoline direct injection (GDI) engines, rapidly forming uniform and finely atomized particles is essential,which demands a spray with proper penetration length,dispersed droplet distribution,and small droplets[2].Although increasing the injection pressure can meet these requirements[3],it leads to a significant shortcoming known as wall-wetting caused by its longer penetration [4,5].On the other hand,flash boiling atomization(FBA)is a promising approach for improving the atomization process.Flashing of hot liquid resembles a boiling process characterized by thermal non-equilibrium(temperature difference across the interface).When the liquid’s temperature is higher than its boiling point, or the local pressure is lower than its saturation pressure,the vapor bubble explosion formed in the superheated liquid can promote the atomization process, thus significantly improving the atomization performance.As the viscosity, surface tension, and latent heat of vaporization of the fuel droplets all decrease with the increase in temperature [6–9], flash boiling causes a rapid phase transition from liquid to vapor.This phase transition leads to volume expansion and changes in spray angle,spray width,droplet size,and penetration,enabling flash boiling to produce finer droplets and even particle distribution.

Many researchers have studied FBA from theoretical,experimental,and numerical aspects.The theoretical studies have investigated the flash boiling process’s theories,including flashing inception,nucleation,bubble growth[10,11],evaporation wave[12,13],and spray collapse[14,15].These works explored the physics and modeling of flashing flows, investigated evaporation waves in liquid jet flashing regimes, and examined spray collapse mechanisms under varied conditions.Despite these efforts, understanding physical processes, such as bubble formation in superheated liquid, remains insufficient.In experimental examinations,numerous studies have focused on the effects of superheating degree and injector configuration on flash boiling spray characteristics,visualizing the atomization process to study motion and features.The most commonly used experimental methods include phase Doppler anemometry (PDA) [16], phase Doppler particle analyzer (PDPA) [17], schlieren[18], Mie scattering [19], laser-induced exciplex fluorescence (LIEF) [20,21].These investigations shed light on flash-boiling spray collapse mechanisms, plume-to-plume interactions, and distinctive characteristics of flash-boiling sprays compared to traditional liquid sprays.The findings contribute valuable insights for optimizing fuel atomization and combustion processes in direct-injection engines.However, experiments were hindered by the high cost of equipment and materials.Limited samples made predicting flash boiling phenomena under unknown conditions challenging,and preparing for experiments required significant time investment.

Since the theory of FBA is quite complicated,covering many physical processes,and experiments under high temperature and high pressure are hard to control,more investigations have been achieved through numerical simulations [22,23].These works collectively explored various aspects of flash boiling phenomena in fuel sprays,focusing on gasoline direct injection(GDI)injectors.They covered theoretical, experimental, and numerical studies, analyzing factors such as superheating degree and injector configuration [24,10].The collapse mechanisms of flash boiling sprays were investigated through multiple perspectives, including the velocity field, temperature field, vapor concentration field, and droplet diameter.Additionally, the review highlighted the importance of understanding flash boiling spray behaviors for applications like injector design and injection strategies.Another study addressed the modeling challenges of flash boiling effects in sprays,emphasizing the need for an economical and computationally efficient methodology.In these studies, the Eulerian-Eulerian and Eulerian-Lagrangian methods are typical for simulating fuel spray.In Eulerian-Eulerian methods[25,26],the gas and liquid phases are calculated based on the Eulerian description,where the internal flow and near-field spray can be obtained directly.The studied problem domain includes the nozzle and spray field to predict the flash boiling process considering nucleation,bubble growth,and evaporation[27,28].Numerous interface tracking methods can be utilized to capture the spray,such as the volume of fluid (VOF) [29,30], piecewise-linear interface calculation (PLIC) method [31], level-set method,and VOF level-set (VOF-LS) [32] methods, and others.However, these methods require very dense meshes to clearly trace the gas-liquid interfaces, making them computationally expensive.In recent years,a method based on the Euler framework has also experienced rapid development,providing the capability to combine the Euler and Lagrangian methods using a transition mechanism[33–35].The initial jet and primary breakup are predicted using the Euler method on a sufficiently fine mesh; if a liquid block satisfies the transformation conditions,it will be converted into Lagrangian particles,simplifying calculations.However,this method still involves significant computational complexity and is prone to divergence.Moreover, most liquid is atomized into droplets at the outlet in the flash atomization process, with almost no continuous liquid phase present.Therefore, this method is not suitable for the current simulation work.

In order to save computational costs,more simulations of flash boiling spray resort to Eulerian-Lagrangian methods.Models based on Eulerian-Lagrangian descriptions include primary atomization models[36,37],secondary atomization models[35],evaporation[38,39]and boil models[40,41],and heat and mass transfer models [42,43], exploring phenomena like flash-boiling, droplet dynamics,spray structure,and the intricate interplay between liquid and gas phases[44,45].The methodologies included numerical simulations,machine-learning algorithms,parametric studies,and novel models.An et al.[46], Price et al.[40], and Zhou et al.[39] investigated the influence of the evaporation model on numerical results and conducted qualitative studies on the results,providing a reference for subsequent simulation work from an academic perspective.Liao and Lucas[10]reviewed the models for flashing inception,nucleation,bubble growth,vapor generation,and interfacial area density.They noted that the models for representing the formation,growth,and breakage of bubbles in superheated liquid are still far from the requirements of engineering applications.For instance,many models are only applicable under ideal conditions or are more or less empirical, resulting in low accuracy of numerical simulations.For instance,the Homogeneous Relaxation Model(HRM)treats the transition of the thermodynamic system from non-equilibrium to equilibrium as a relaxation process.The two states are bridged by an empirical coefficient,i.e.,the relaxation time[10].

When the outlet speed of the nozzle exceeds 400 m/s,the high-speed FBA requires a very dense mesh to trace the flow field information of the gas and droplets and the gas-liquid interfaces.Millions or even billions of grids should be introduced,and the data post-processing becomes quite burdensome[25,47].Due to the violent liquid flashing,a high-speed gas-liquid mixture is formed near the exit of the nozzle.Moreover, the physical mechanism of high-speed flash boiling is very complicated.The existing models often suffer from low simulation accuracy for high-speed FBA.Generally,the errors may exceed 15%[37].Besides,ambient pressure and temperature are believed to play important roles in atomizing superheated liquids.Nevertheless,most studies only consider one of these factors[10].

This study establishes a high-speed flash boiling atomization(HFBA)model for fuel atomization,catering to scenarios where the outlet speed exceeds 400 m/s,reaching up to approximately 700 m/s.This model encompasses the processes of primary atomization, secondary atomization, and the interaction between gas and liquid during atomization.The Eulerian method is employed for modeling the gas, while the Lagrangian method is utilized for simulating the liquid.These methods enable effective capture of the droplets and gas-liquid interfaces.Implementing an adaptive mesh enhances the resolution of the gas-liquid interfaces, improving the model’s accuracy.The model operates with fewer than one hundred thousand grids, maintaining a solution error below 10% compared to experimental data [48].Subsequently, this validated model is applied to examine the effects of fuel injection temperature and ambient pressure on various aspects,such as the velocity field,temperature field, vapor mass fraction distribution, particle size distribution, and spray penetration length.This comprehensive application delineates the entire process and intricacies of the HFBA.

2 Theoretical Model

Considering the computational costs,the Eulerian-Lagrangian method is adopted for fuel atomization simulations.Compared to the Eulerian-Eulerian method, this method simplifies capturing the flow field, particle size, and other information post-atomization and facilitates the design and optimization of the numerical model[36,37].The spray mixing of liquid fuel in a combustion chamber is a transient, complex, multiphase flow process.This process involves the formation of droplets,fragmentation, collision, aggregation, boiling, evaporation, vaporization, and the mutual coupling with air,integrating and superimposing many processes[41].

2.1 Eulerian Description for the Gas

The flash boiling spray of fuel is a complex gas-liquid two-phase flow process, encompassing turbulent flow,gas-liquid heat transfer,and mass transfer processes.These physical processes adhere to basic conservation laws,including mass and momentum conservation.For the heat transfer problem,the law of energy conservation and the equation of state of the ideal gas are also applicable[49,50].

(1)Mass conservation equation

In this study,the three-dimensional problem is simplified to be a 2D axisymmetric problem,wherexis the axial coordinate,yis the radial coordinate, as shown in Fig.1.For this 2D axisymmetric geometry problem,the continuity equation is expressed as

whereu=is the velocity vector of the gas in this 2D problem;tis the time,ρis density of the gas;Sgis the source term for the gas-liquid interactions.

Figure 1:Coordinate conversion

(2)Momentum conservation equations

For the 2D axisymmetric problem,the momentum conservation equation is given by

wherepis the pressure of the gas;Fis the body force vector representing the external volume forces and the other model-dependent source terms such as porous media and user-defined sources.μeff=μ+μt,whereμis dynamic viscosity of the gas, andμt=ρCμk2/εis the turbulence viscosity of the gas in whichCμ= 0.09 is a constant,andkandεcan be solved from the turbulent equation which will be introduced in the following paper.

(3)Energy conservation equation

The energy conservation equation can be expressed as

whereiis theith species,Eis the fluid energy;Jiis the diffusion flux of theith species;biis the enthalpy;Scis the absorption and exotherm from the chemical reaction and other specified volumetric heat sources.

(4)Vapor-liquid equilibrium equation

The vapor-liquid equilibrium equation is described as follows:

whereKiis the volume concentration of theith species,which includes vapor,nitrogen,and n-pentane in this paper;Riis the net production rate due to the chemical reaction;Siis the additional production rate due to the discrete phases(droplets)and the specified source terms.

(5)Equation of state of the ideal gas

For the ideal gas,the equation of state can be expressed by

whereξ0= 8.314 is a constant;Tis the temperature of mixed gas;ldenotes the nitrogen and fuel vapor;Wlis the molecular weight.

2.2 Lagrangian Description for the Droplets

This study simulates the FBA and evaporation processes of fuel using the discrete phase model(DPM)in the Lagrangian coordinate system.The DPM model captures various physical processes in the FBA process,including liquid spraying into the combustion chamber,droplet formation,and the exchange of momentum,mass,energy,and species transfer between droplets and gas[36,38,41].

2.2.1 Motion of Droplets

The droplet is considered a discrete particle,and the particle trajectory can be predicted by solving the equation of motion under the Lagrangian reference system.The equation of motion is given as follows[51]:

whereuis the velocity of the gas;mpis the mass of the droplet particle;upis the velocity of the particle;ρpis the density of the particle, which changes with the ambient temperature;Fpis the additional force on the particle.In fact, additional forces include the “virtual mass” force, which is the force required to accelerate the fluid surrounding the particle.This force can be written as:whereCVMis the virtual mass factor with a default value of 0.5,an additional force arises due to the pressure gradient in the fluid:F=mpρ/ρpu∇ump(u-up)/τris the drag force.Hereτris the relaxation time of the particle which can be calculated as follows:

whereμpis the dynamic viscosity of the droplet;Dpis the diameter of the particle;Cdis the drag coefficient;Re is the Reynolds number of gas denoted as follows:

2.2.2 Heat and Mass Transfer Model

(1)Droplet heating and cooling

When the temperature of the droplet is below the defined evaporation temperatureTe, the law of heating or cooling should be satisfied as follows[52].The evaporation temperature is set to 223 K[53].WhenTp

wherecpis the heat capacity of the particle;Apis the surface area of the particle;ζis the coefficient of the convective heat transfer;εpis the emissivity of the particle;σis the Stefan-Boltzmann constant;TRis the radiation temperature.In Eq.(7), the coefficient of the convective heat transferζcan be calculated by the following equation[53]:

where Nu is Nusselt number;qtis the thermal conductivity of the gas; Pr is Prandtl number of gas;Re is the Reynolds number based on the particle diameter and the relative velocity which has been defined in Eq.(6).When a droplet passes through the Eulerian mesh, it causes a change in the heat content of the surrounding gas.This heat exchange acts as a source termScin the gas energy equation,as presented in Eq.(1).

(2)Evaporation of droplets

The evaporation law can be utilized to evaluate the evaporation of a droplet in the discrete phase when the temperature of the droplet satisfies the following conditionsTe≤Tp≤Tb.HereTbis the boiling temperature.The flux of droplet vapor into the gas phase depends on the difference of the vapor concentration at the droplet surface and the bulk gas,which is expressed as

whereNis the molar flux of gas;ςis the coefficient of mass transfer;Cis the vapor concentration in the gas phase;Csis the vapor concentration at the surface of droplet.The mass transfer coefficientςin Eq.(9)can be calculated from the Sherwood number as below[53]:

where Sh is the Sherwood number;βis the diffusion coefficient of gas; Sc is the Schmidt number,which can be expressed as Sc=μ/ρβ.

During the evaporation of the droplet,the reduction of the mass of the droplet is according to the following equation:

whereΔtis the time increment for each time step;Mis the molecular weight of gas.

WhenTe≤Tp≤Tb, finally, the droplet temperature is updated based on a heat balance that relates the sensible heat change in the droplet to the convective and latent heat transfer between the droplet and the continuous phase[52,53]

where dmp/dtis the evaporation rate of the droplet;ϑprepresents the blackness of the droplet;hrepresents the latent heat of the droplet.

(3)Boiling of droplets

When the boiling conditionTp≥Tbis satisfied, the boiling evaporation rate equation can be expressed as follows[54]:

wherecis the heat capacity of gas.By solving Eq.(13),Dpin Eq.(6)can be obtained.

2.3 Coupling between Gas and Droplets

In the interactions between gas and droplets, exchanges of mass and momentum occur, as illustrated in Fig.2.Regarding mass exchange,evaporation constantly occurs on droplet surfaces,with vapor entering the gas field as a source term for the continuous gas phase.In terms of momentum exchange, the droplet’s motion trajectory is influenced by the continuous gas field, and the droplet motion also impacts the continuous gas field flow.The equations for mass and momentum exchanges are detailed as follows.

Figure 2:The mass and momentum exchange between the dispersed droplet and the continuous gas

(1)Momentum exchange

The transfer of momentum from the gas to the droplets can be calculated by examining the change in momentum of the droplets when they pass through each control volume.The change in momentum can be described as follows[55]:

(2)Mass exchange

The mass transfer from the droplets to the gas can be evaluated by computing the change in mass of the droplets when they pass through each control volume.The mass changeMecan be expressed as

wheremp0is the initial mass of the particle;dmp0/dtis the initial mass flow rate of the particle injection.The process of mass exchange between two phases can be modeled by including it as a source term in the continuous phase continuity equation.

2.4 Turbulence Model

Turbulence is a nonlinear and complex flow that induces turbulent fluctuations, resulting in mass and momentum exchange and concentration changes between fluid media.Common numerical methods for modeling turbulence include the direct simulation method(DNS),large eddy simulation(LES),and Reynolds-averaged Navier-Stokes(RANS)method.DNS directly solves the Navier-Stokes equation without a turbulence model, encompassing the flow field from the turbulent dissipation scale (Kolmogorov scale) to the large eddy scale, requiring a very small time step and refined grid scale [56,57].LES solves the equation directly for large eddies, using models for small-scale eddies[58].Although DNS and LES can achieve accurate simulations, their computational processes are complex and costly.Hence, they are typically used for simple and ideal flow fields.The unsteady Reynolds-averaged Navier-Stokes method(URANS)decomposes flow parameters into time-averaged and pulsation quantities, then reintegrates these values into the Navier-Stokes equation solution,offering cost-effectiveness and good accuracy[59,60].This study selects the k-epsilon turbulence model in URANS to solve two transport equations for turbulence determination[61].This model effectively simulates the diffusion process in jet phenomena and handles locations with large flow field parameter gradients[60].The transport equations of the k-epsilon model are presented as follows:

whereC1= max[0.43,η/η+5],η=Sk/ε,S=Sij=(∂uj/∂xi+∂ui/∂xj)/2 defines the deformation rate tensor;kis the turbulence kinetic energy of the gas,εis the turbulence dissipation rate of the gas;νis the kinematic viscosity of the gas;kgis the turbulent kinetic energy due to the mean velocity gradient,kbis the turbulent kinetic energy due to the effect of buoyancy;KMis the effect of the amount of expansion pulsation on the turbulent kinetic energy in the compressible turbulence;C1ε= 1.44,C2ε= 1.9,andC3εare the constants;αk=αε≈1.393 are the turbulent Prandtl numbers,respectively;SkandSεare the source terms[26].

2.5 Initial Conditions Based on the Primary Atomization Model

In the HFBA process,the fuel undergoes intense heat and momentum exchange with surrounding air during the primary atomization period.As a critical component of the atomization process, the accuracy of the entire simulation depends on the model selection for this period[36,62,43].This work chooses the effervescent atomizer model for primary atomization,accurately presenting information on velocity, temperature, droplet size distribution, and the flow field around the nozzle post-initial atomization.

Effervescent atomization is a process where a jet of superheated liquid is injected, and upon leaving the nozzle, the volatile liquid rapidly boils and undergoes phase change.This sudden phase change effectively breaks the continuous liquid stream into small droplets with a wide atomization angle,forming a droplet group.In the model,the initial droplet velocity is calculated using the law of conservation of mass.It is assumed that the nozzle exit’s cross-sectional area isCcttimes the area of the nozzle area,where times the area of the nozzle area whereCctis a constant 0.611,as shown in the following equation:

whereu0=u(t=0);Ais the cross-sectional area of the nozzle;meffis the effective mass flow of the nozzle.

The maximum droplet diameter is set to be the effective diameter of the nozzle

whereDnis the nozzle outlet diameter;is the median diameter which depends on the angle between the orbit of the droplet and the direction of the jet and expressed by

whereθis the spray cone angle;ψis calculated from the gas mass fractionfand the dispersion coefficientCeff.

and

2.6 Secondary Atomization

Post-primary atomization,large droplets may further fragment in subsequent motions.The WAVE model, which is applicable to high Weber number problems, is introduced in this paper to simulate the HFBA process of fuel [63–65].The WAVE model posits that droplet breakup time and the size distribution of broken droplets relate to the fast-growing Kelvin-Helmholtz instability wave.The dispersion equation is presented as follows[26]:

whereωis wave growth rate;νis the kinematic viscosity of particle;κis the wavenumber defined asκ= 2π/λ;λis the wave length;Lis another wavenumber defined asL= √Inis the nth order modified Bessel function of the first kind;rjis the radius of liquid jet;σsis the surface tension coefficient of liquid;ι=√

The most unstable growth rateΩand wavelengthΛ/rjcan be obtained as

The radius of the newly-formed droplets is proportional to the wavelength of the fastest-growing unstable surface wave on the parent droplet

whereB0is a model constant set to be 0.61.

3 Numerical Simulations

In this section, several examples of numerical simulations for the HFBA process are presented.The experimental data from Li et al.’s study[48]serves as the reference solution.In these experiments,a single-component fuel,n-pentane,is used for the HFBA tests.The fuel line is heated around the fuel rail using a high-pressure nozzle from a gasoline direct injection (GDI) engine, achieving fuel flash boiling in the test chamber by heating the fuel jet.The experimental conditions are listed in Table 1.

Table 1: Experimental conditions for flash boiling injections

3.1 Mesh Independent Analysis and Model Validation

In the numerical simulations, a 2D axisymmetric geometry problem is studied, as depicted in Fig.3.The problem domain extends 70 mm in length and 50 mm in width.The standard wall function [50,64] is applied to the boundary, which is positioned far from the nozzle and minimally affects the flash spray process.The detailed numerical settings for the simulations are listed in Table 2.Superheated fuel is injected from the fuel inlet into a test chamber filled with nitrogen at room temperature.Under varying ambient pressures and injection temperatures, the fuel undergoes changes like boiling,fragmentation,droplet collision,and evaporation,causing it to break into smaller droplets,ejected from the spray field outlet.

Table 2: Numerical setting

Figure 3:Coordinate conversion

For the finite volume method,the accuracy of simulation results depends on mesh distribution.Near the nozzle, where a large number of liquids boil and vaporize instantaneously, forming an extremely high-velocity jet stream with drastic changes in temperature, momentum, and species concentration,mesh refinement is crucial.As shown in Fig.4,locally encrypted mesh with a growth rate of about 1.05 is employed for this example,covering an area approximately 15 mm long and 10 mm wide.The solutions converge under this mesh refinement, as demonstrated in Fig.5, where the injection temperature is 333 K and the ambient pressure is 20 kPa.When the number of elements exceeds 70200, the solutions become convergent, and beyond this point, they are almost meshindependent.Therefore,the problem domain is discretized in the simulations with 91,200 elements.

Figure 4:Mesh distribution

Figure 5:Solutions of pressure and temperature at coordinate(x=50 mm,y=0 mm)with different quantities of mesh:(a)Pressure;(b)Temperature

In order to validate the proposed atomization model, experimental results of n-pentane fuel atomization from the literature[48]are introduced as the reference solution,and the entire simulation process is illustrated in Fig.6.The numerical solutions are compared with the experimental results in Fig.7,where the error is within±10%.The SMD is the average equivalent diameter,as defined in reference[48].

wherenis the number of droplets.The solution accuracy can satisfy the engineering problems’requirement,proving that the proposed HFBA model is reliable.

Figure 6:Simulation flow chart

Figure 7:SMD in different ambient pressures

Fig.8 displays the HFBA of fuel under an injection temperature of 379 K and ambient pressure of 60 kPa,the highest injection temperature in all the simulations presented in this work.The velocity contour reveals a microburst observable in the experiments[48],achieving a maximum exit velocity of 661.7 m/s.

Figure 8:The velocity contour downstream of the nozzle under T=379 K and P=60 kPa

3.2 Analysis of the Influence Factors on Atomization Performance

For the FBA,the ambient-to-saturation pressure and the degree of superheating can be utilized as characteristic parameters to describe the phase transformation of the flash boiling spray[6,20].

3.2.1 Different Ambient Pressures

(1)Solutions of gas under different ambient pressures

As depicted in Fig.9,when the injection temperature is set at 333 K,the velocity and temperature of the gas at positions along the x-axis gradually decrease as the ambient pressure decreases.A noticeable change is observed near the nozzle due to the intensified flash evaporation caused by the decrease in ambient pressure.Near the nozzle, more diffusion of liquid fuel occurs in the radial direction, and the droplet momentum has a larger component in the radial direction compared with traditional atomization.These dynamics result in reductions and fluctuations of velocity in the axial direction in a small region.Higher ambient pressures introduce higher velocities,although the differences are not quite pronounced.Moreover,the droplets become more dispersed due to the drop in velocity, leading to a temperature drop at the positions along the x-axis.The lower the ambient pressure,the greater the temperature change.

Figure 9: Velocity and temperature of the positions on the x-axis under different ambient pressures:(a)Velocity;(b)Temperature

As shown in Fig.10, under conditions of an injection temperature of 333 K and an ambient pressure of 10 kPa, when the liquid fuel is sprayed from the nozzle into the flash chamber, the gas velocity and the fuel concentration reach their maximum values during a brief period.This occurs because the superheated liquids in the nozzle contain a large amount of internal energy,which causes the violent boiling of a large amount of liquid fuel within a specific range.The boiling generates an“explosion”that leads to a sharp increase in the velocity of the surrounding gas and enhances the fuel vapor concentration near the x-axis.

Figure 10:Velocity and mass fraction of n-pentane vapor of gas for different cross sections:(a)Velocity;(b)Mass fraction of n-pentane vapor

As the spray continues to move,the maximum velocity and vapor concentration along the x-axis direction gradually decrease while the mass fraction increases in the y-axis direction(radial direction).This phenomenon occurs since work is required to overcome air resistance during the spray motion,resulting in energy consumption.As a result,the velocity decays in both the axial and radial directions,accompanied by the diffusion of vapor in the radial direction.

(2)Solutions of droplets under different ambient pressures

In order to investigate the influence of ambient pressure on the size of atomized droplets,the SMD of droplets under an injection temperature of 333 K and ambient pressures of 5,10,20,30,60,and 100 kPa are compared,as shown in Fig.11.The results demonstrate that the SMD of droplets tends to increase with the growth of ambient pressure.

The size distribution of the droplets under ambient pressures of 5,10,20,30,60,and 100 kPa is presented in Fig.12.The results indicate that lower ambient pressures lead to smaller droplet sizes and more concentrated particle size distributions.This is because instantaneous boiling occurs when the nozzle sprays the superheated jet liquid.Under lower ambient pressures,the saturation temperature decreases,and the liquid superheat increases,driving the flash vaporization process.In addition,the lower resistance to gas expansion in the atmospheric environment at lower pressures promotes bubble formation through boiling.The expansion force generated by these bubbles on the liquid fuel causes strong disturbance forces,further enhancing the jet’s expansion,leading to easier droplet breakup and smaller atomized SMD particle sizes.

Figure 11:Variation of atomization particle size with ambient pressures

Figure 12: (Continued)

Fig.13 illustrates the characteristic particle size parameters of n-pentane under different ambient pressures.D10 represents the diameter below which 10%of the liquid volume is atomized into smaller droplets,and the same applies to D30 and D50.A comparison of the results reveals that the particle size parameters, including D10, D30, and D50, are significantly smaller at low ambient pressures compared to high ambient pressures.This observation indicates that decreasing the ambient pressure is an effective way to reduce the size of atomized particles.

Figure 13:Characteristic parameter values of particle size under different ambient pressures

Furthermore, as displayed in Fig.14, under an ambient pressure of 100 kPa, the percentage of droplets with a diameter of less than 5 μm is 77.1%of the total droplets formed by atomization.Under an ambient pressure of 60 kPa,this proportion increases to 85.7%and 92.2%at an ambient pressure of 10 kPa.The number of small droplets generated by the HFBA of fuel increases with decreasing ambient pressure, resulting in a higher proportion of smaller droplets.Conversely, increasing the ambient pressure leads to a higher proportion of large droplets and raises the average particle size of atomization.These results agree with the statistical results for SMD presented in Fig.11.Therefore,reducing the ambient pressure is beneficial for liquid fragmentation atomization and can improve atomization efficiency.

Figure 14:Cumulative distribution of droplet particle size under different ambient pressures

Fig.15 exhibits the spray penetration depth variation in HFBA with time under different ambient pressures.As time elapses,the variation trends remain the same,but the growth rates decelerate in all three cases.As the ambient pressure decreases,the atomization becomes more violent,resulting in a“microburst effect”that causes the fuel to break into smaller droplets.These smaller droplets have less mass and volume,which disrupts the original trend of the airflow movement and reduces the axial penetration velocity.Therefore,the penetration depth of atomization spray under low ambient pressure is much smaller than that under high ambient pressure.In addition,as the ambient pressure decreases,the surface area of the total droplets increases,and overall air resistance also increases,which reduces droplet momentum and velocity.Consequently,the penetration depth is naturally reduced.

Figure 15: Penetration distance of droplets on the position of the x-axis with time under different ambient pressures

3.2.2 Different Injection Temperatures

(1)Solutions of gas under different injection temperatures

As presented in Fig.16, when the ambient pressure is set to 60 kPa, the fuel injection velocity(velocity of the gas)increases significantly with the rise of injection temperature.At the same ambient pressure,a higher temperature results in a more pronounced flash effect,and the gas expansion leads to a higher exit velocity of the fuel.However,due to the momentum exchange with the surrounding air,this velocity gradually decreases during subsequent motion and ultimately converges to a similar velocity at x=50 mm under different injection temperatures.Regarding the temperature distribution,liquid fuel atomization is more complete at higher injection temperatures,and the dispersed droplets accelerate the heat exchange efficiency between the droplets and surrounding air.As a result, under higher injection temperatures,a more pronounced temperature drop of the atomized droplets along the x-axis can be observed.As shown in Fig.17,under conditions of an injection temperature of 379 K and an ambient pressure of 60 kPa,the gas velocity and the fuel concentration reach their maximum values during a brief period, where the trend is the same as in the previous section under different ambient pressures.

Figure 16:Velocity and temperature of the positions on the x-axis under different injection temperatures:(a)Velocity;(b)Temperature

Figure 17:Velocity and mass fraction of n-pentane vapor of gas for different cross sections:(a)Velocity;(b)Mass fraction of n-pentane vapor

(2)Solutions of droplets under different injection temperatures

In order to investigate the differences in atomized droplet size at 50 mm from the nozzle under different injection temperatures, the droplet SMD at an ambient pressure of 60 kPa and injection temperatures of 293,313,333,343,353,and 397 K are compared,as presented in Fig.18.The results agree with the available experimental data[48],which once again validates the high simulation accuracy of the proposed HFBA model.

Figure 18:Variation of atomization particle size with injection temperature

As shown in Fig.19,the SMD values of droplets decrease gradually with the increase in injection temperature.This indicates that the injection temperature can affect the atomization fragmentation.This is because a higher liquid injection temperature can result in larger evaporation under a certain flash ambient pressure.In addition, increased injection temperature promotes bubble formation and growth, which helps break up the liquid, improve liquid atomization quality, reduce the spray droplet size,and increase the flash evaporation degree.Furthermore,when the injection temperature is increased,the latent heat of the vaporization of water decreases and further facilitates the evaporation process.Therefore, increasing the injection temperature can cause the liquid fuel to break up and strengthen the atomization process.Fig.20 displays the particle size characteristic parameters of npentane fuel at varying injection temperatures.The results reveal that at higher injection temperatures,the particle size parameters, including D10, D30, D50, and SMD, are notably smaller than those achieved at lower injection temperatures.This indicates that increasing the injection temperature can effectively enhance liquid fuel atomization.

As displayed in Fig.21, at an injection temperature of 293 K, the droplets, with a diameter of less than 8 μm formed by atomization, account for approximately 80.4% of the total number of droplets.At an injection temperature of 333 K, about 95.9% of the droplets are less than 8 μm in diameter.As the injection temperature increases to 379 K,the proportion of droplets less than 8 μm in diameter increases to approximately 98.7%.The cumulative curves for the higher temperature working conditions of 333 K and 379 K show that more than 90%of droplets are less than 6 μm in diameter,which is significantly larger than the diameter of droplets under the working condition at 293 K.These results suggest that increasing the injection temperature leads to more complete atomization,and the number of small diameter droplets continues to increase,accounting for a larger proportion.This is beneficial for the atomization process.Fig.22 shows the atomization spray penetration depth changes over time under various injection temperatures.Although higher temperatures result in more intense flash vaporization, more pronounced radial expansion, and weaker axial expansion, the ambient temperature does not significantly affect the droplet penetration distance.

Figure 19:Distribution of particle size ratio under different injection temperatures

Figure 20:Characteristic parameter values of particle size under different injection temperatures

Figure 21:Cumulative distribution of droplet particle size under different injection temperatures

Figure 22:Penetration distance of droplets on the position of x-axis with time under different injection temperatures

4 Conclusion

This study proposes an HFBA model to simulate the atomization process under high outlet speed.Compared to the experimental results,the maximum SMD error from the proposed model is less than 8%when the numerical model’s boundary conditions align with the experiments’operating conditions[48].This error rate is significantly lower than the 15–20% error observed in traditional models.Utilizing fewer than one hundred thousand grids can achieve mesh-independent results, markedly more efficient than conventional simulations that require over a million grids.

The results demonstrate that HFBA can generate extremely fine and uniformly distributed spray.During flash evaporation, vapor generation and bubble growth stem from mass transfer induced by interphase heat transfer, with the vapor generation rate being regulated by thermal imbalance.Increasing the initial injection temperature or decreasing the ambient pressure enhances the liquid’s superheat, intensifies the thermal imbalance, and amplifies energy instability within the liquid, thus accelerating vapor generation.Additionally, higher superheat fosters bubble nucleation density and frequency, precipitates the onset of flash vaporization, and reduces droplet size, thereby driving the flash vapor.Incremental increases in injection temperature or reductions in ambient pressure lead to a gradual increase in the proportion of smaller droplets and a corresponding decrease in D10,D30,D50, and SMD, culminating in finer and more uniform atomized droplets.Furthermore, lowering the ambient pressure diminishes the axial velocity of atomized droplets,enhancing their combustion efficiency in the chamber.This HFBA model shows significant potential for extensive application in simulating high-velocity fuel atomization challenges.

Acknowledgement:The authors express their sincere gratitude to Professor Xuesong Li of the School of Mechanical Engineering, Shanghai Jiao Tong University, for providing invaluable insights and details crucial to the success of the experimental work in this paper.Special thanks also extend to the dedicated staff at Beijing Mechanical Equipment Research Institute for their valuable input and guidance, which contributed to the enhancement of this research.Their collaborative efforts have greatly enriched the quality and depth of their work.

Funding Statement:This work is supported by the National Natural Science Foundation of China(Project Nos.12272270,11972261).

Author Contributions:The authors confirm their contribution to the paper as follows: Study conception and design: Wei Zhong, Lihua Wang; data collection: Wei Zhong, Zhenfang Xin, Lihua Wang,Haiping Liu;analysis and interpretation of results:Wei Zhong,Lihua Wang;draft manuscript preparation:Wei Zhong,Lihua Wang.All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials:The material data utilized in this study have been sourced from Google,which provides a comprehensive platform for various datasets,and the experimental data[48]presented in this article have been cited in the main text.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.