Next Article in Journal
On the Effectiveness of Scale-Averaged RANS and Scale-Resolved IDDES Turbulence Simulation Approaches in Predicting the Pressure Field over a NASCAR Racecar
Previous Article in Journal
Nanobubble Collapse Induced Erosion near Flexible and Rigid Boundaries: A Molecular Dynamics Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation and Numerical CFD Assessment of a Thermodynamic Breakup Model for Superheated Sprays with Injection Pressure up to 700 Bar

1
Industrial Engineering, Information and Economics Department, Università degli Studi dell’Aquila, Piazzale Ernesto Pontieri, Monteluco di Roio, 67100 L’Aquila, Italy
2
Consiglio Nazionale delle Ricerche, Istituto di Scienze e Tecnologie per l’Energia e la Mobilità Sostenibili (STEMS), Via G. Marconi 4, 80125 Napoli, Italy
*
Author to whom correspondence should be addressed.
Fluids 2023, 8(5), 155; https://doi.org/10.3390/fluids8050155
Submission received: 8 April 2023 / Revised: 26 April 2023 / Accepted: 10 May 2023 / Published: 14 May 2023
(This article belongs to the Topic Computational Fluid Dynamics (CFD) and Its Applications)

Abstract

:
Among the most relevant fields of research recently investigated for improving the performance of gasoline direct injection (GDI) engines, there are ultrahigh injection pressures and the flash-boiling phenomenon. Both perform relevant roles in improving the air/fuel mixing process, reducing tailpipe emissions and implementing new combustion methods. When a high-temperature fuel is released into an environment with a pressure lower than the fuel’s saturation pressure, flash boiling occurs. Due to complex two-phase flow dynamics and quick droplet vaporization, flash boiling can significantly modify spray formation. Specifically, if properly controlled, flash boiling produces important benefits for the fuel–air mixture formation, the combustion quality and, in general, for overall engine operation. Flash boiling was broadly investigated for classical injection pressure, but few works concern ultrahigh injection pressure. Here, the investigation of the spray produced by a multihole injector was performed using both experimental imaging techniques and CFD simulations aiming to highlight the combined impact of the injection pressure and the flash boiling occurrence on the spray morphology. The shadowgraph method was employed to observe the spray experimentally. The information gathered allows for assessing the performances of an Eulerian–Lagrangian algorithm purposely developed. Breakup and evaporation models, appropriate for flashing sprays, were implemented in a CFD (Computational Fluid Dynamics) code. The experimental results and the CFD simulations demonstrate a good agreement, demonstrating that through adoption of a flash-boiling breakup model, it is possible to reproduce non-evaporating and superheated sprays while changing few simulation parameters. Finally, the results also show the significance of injection pressure in preventing spray collapse.

1. Introduction

One of the most promising approaches to reduce the environmental impact of the light-duty transportation sector is generally acknowledged to be the employment of direct-injection spark-ignition propulsion systems gasoline-fueled (GDI), which enable engine downsizing and the use of improved combustion techniques [1,2]. However, further improvements are required to boost GDI engines performance and limit, among the other, the liquid film formation that triggers soot production [3]. The injection process plays a significant role in GDI engines and significantly impacts the chain of events that leads to mixture formation, combustion and pollution creation. Direct injection techniques into the cylinder require that the liquid fuel is broken up into tiny droplets and quickly evaporated in order to ensure an effective mixture formation [4,5,6].
It was shown that both high injection pressure as well as flash boiling improve spray characteristics [7]. According to previous scientific studies, fuel injection at high pressure levels (up to 1500 bar) enhances both fuel evaporation and in-cylinder turbulence giving, at the end, a more stable combustion [8]. The injection pressure affects the breakup dynamics and droplet dimensions, as shown by Yamaguchi in an in-depth investigation on spray evolution [9]. In another work, spray penetration and breakup time of high-pressure sprays were related to literature correlations for diesel injections [10]. On the other hand, liquid atomization by flash boiling mechanism is another potential method for achieving high-quality spray with a small and narrow drop-size distribution [3,11,12,13,14,15,16].
When the liquid fuel at high pressure and temperature (near its boiling point) is discharged through an orifice to a low ambient pressure (which is below its saturation pressure), flash boiling sprays are generated. The compressed liquid is expanded below its saturation pressure; vapor bubbles arise, and their continued expansion will produce severe spray atomization.
However, when flashing sprays are issued from multihole injectors, the single plumes deviate toward one another and combine into a single jet. This phenomenon is the so-called “spray collapse” and frequently takes place in flashing sprays. A collapsed spray causes spray-wall impingement (which increases the generation of exhaust soot) and it greatly reduces the combustion quality [7,17,18].
Flash boiling onset and injection pressure are intimately coupled. High injection pressure does more than just produce sprays with small drop sizes; it also prevents spray collapse and facilitates the vaporization of liquid fuel [19,20]. Therefore, it seems inevitable to develop a detailed understanding of the flash boiling injections, exploiting both experimental and numerical methods—in particular, when extremely high injection pressures are adopted. Eulerian–Lagrangian spray simulations are a smart and efficient tool for simulating the evolution of liquid sprays with limited calculation times that are in line with industry standards and that fit with the integration of such spray models in whole engine simulations. The breakup model is an important component of this approach. Conventional approaches, such as Kelvin–Helmholtz–Rayleigh–Taylor (KH–RT) and Reitz–Diwakar (RD), are based on the idea that liquid droplet breakage is caused by superficial aerodynamic instabilities and is primarily dependent on the relative velocity between liquid and gas phases. However, when a superheated liquid is injected, the vapor bubbles that form inside the liquid core also contribute to the breakup. Therefore, without adjusting certain constants, aerodynamic breakup techniques alone cannot replicate the various flash-boiling situations and, in particular, provide unrealistic droplet sizes [21,22]. An effervescent breakup model developed expressly to account for bubble expansion within droplets in flash-boiling sprays exhibited significant improvement for Lagrangian spray simulations [23] and so is strongly advised. It is a distinct primary breakup mechanism that interacts with the aerodynamic breakup and considers droplet shattering as a result of thermodynamic nonequilibrium conditions.
So, this paper focuses on the investigation of high-pressure liquid sprays in both non-evaporative and flash boiling conditions performed with experimental and numerical methods. The spray morphology was investigated under various thermodynamic conditions, and, in particular, the performance of the developed effervescent breakup model in simulating the behavior of high-pressure spray under flash boiling conditions was assessed. Common GDI injectors operate with a rail pressure of roughly 200 bar while current contribution deals with values up to 700 bar. The discussed CFD methodology was validated using the experimental images obtained using nonintrusive techniques. The outcomes of the research highlight that the adoption of the effervescent breakup model for flashing spray improves the result’s quality and avoid further tuning of the breakup parameters. The work is organized into three sections: Research Methodology, Discussion of the Result and Conclusions. The Research Methodology is structured into two subsections. The first one illustrates the experimental measurement performed; next, the CFD code and the effervescent breakup model are presented. Results are organized into three subparagraphs. Initially, the modeling approach adopted was assessed comparing the numerical results with the experiments. Next, details concerning the effervescent breakup model are discussed. Last, deepenings about the spray morphology and the role played by flash boiling onset are presented.

2. Research Methodology

A stepped-hole valve-covered orifice (SVCO) injector with 10 holes and 90 ° outer cone angle was considered in this study. Table 1 lists the key characteristics of the device.
The conditions investigated are summarized in Table 2.
Cases HP (high-pressure) and vHP (very high-pressure) examine nonevaporative gasoline sprays, while cases HP-fb and vHP-fb refer to intense evaporation conditions or, more easily, flash boiling. The fuel temperature is increased to 363 K and ambient pressure is lowered to 0.2 bar. Nitrogen inert gas was used to fill the combustion vessel and realize the desired back-pressures. CFD simulations were validated comparing spray penetration length and shape exploiting experimental images acquired using the shadowgraph optical technique. The flash-boiling breakup model developed by the authors for the CFD simulations was implemented in CONVERGE 3.0 as a User Defined Function (UDF).

2.1. Experimental Equipment and Setup

The spray evolution was investigated using Z-type shadowgraph imaging technique, based on the variations of the refractive index. The light source is a pulsed LED (Omicron LEDMOD V2—470 nm/450 mW). Figure 1 presents a scheme of the optical setup: two off-axis parabolic mirrors (f = 508 mm, d = 101.6 mm) are mounted in front of the optical windows of the chamber.
High-intensity flash was used as light source. An high-speed C-Mos camera (Photron FASTCAM SA4) equipped with a 90 mm optics lens was used to acquire spray images. The camera was synchronized with both the injection system and the LED. The imaging was performed at a shutter speed of 16.000 frames per second (fps) with an image dimension of 512 per 448 pixels. The exposure time of the high-speed camera was set to 3.26 μs, while a spatial resolution of 5.9 pixel/mm was obtained by the described setup. The spray images were processed by ad-hoc developed postprocessing procedure [20,24] allowed to identify the contours of the jet in the nitrogen gas and measure the main parameters characterizing the spray structure. The investigations were performed on an optically accessible constant-volume vessel where the thermodynamic conditions of a combustion chamber, at injection time, were reproduced. Full details of the experimental apparatus are reported in [25,26] and hereafter, only a brief description is provided. Fuel injection mainly includes a fuel supply system, which is used to supply fuel from a fuel storage container to the injector. The pressurized fuel was conveyed inside the common rail, where an electronic pressure regulator and a pressure transducer mounted on the rail allowed the setting and reading of the pressure value, respectively, set through a homemade electronic control unit (ECU). The injector was fixed at the top of the CVC in a specific holder, featuring a jacket for the temperature control of the nozzle nose, connected to a chiller for cooling liquid circulation. The fuel was heated by electrical resistance and controlled in temperature by a J-type thermocouple. Nozzle and fuel temperature management was performed with a governor via remote computer. A customized injector equipped with a thin thermocouple was used to calibrate liquid temperature for the jacket and the nozzle nose. Both injector and fuel temperature were kept constant at 293 K for cases HP and vHP while they were increased at 363 K for flashing conditions (Table 2). The injector and the acquisition system were activated in a synchronous mode with an external TTL signal that previously starts the ECU. The acquisitions were repeated five times to take into account for cycle-to-cycle variations.

2.2. Computational Method

The simulations were carried out with an Eulerian–Lagrangian approach. A modified PISO algorithm was adopted for the solution of the governing equations of the continuum-phase [27]. The idea behind the present contribution is to use a fixed parameter setup to simulate both nonevaporative and flash-boiling conditions. For this purpose, the developed effervescent breakup model is activated only in flash-boiling conditions to account for the thermodynamic breakup mechanism together with an improved vaporization model. Therefore, it is appropriate to provide a detailed description of the phase change, heat transfer and flash-boiling breakup models essential for simulating flash-boiling conditions. Among other effects, the nucleation of vapor bubbles within the liquid droplets triggers a thermodynamic breakup mechanism that is approximately ten times faster in comparison with the classical aerodynamic breakup [28,29]. Based on this consideration, an effervescent thermodynamic breakup model, named fbBreakup, was developed as an UDF within CONVERGE CFD code [30]. The secondary aerodynamic breakup was modelled with the Reitz–Diwakar model [31,32].
Figure 2 shows the effervescent breakup of the bubble–droplet system.
The liquid is injected in superheated conditions and so vapor bubbles nucleate within the liquid ligaments [31]. Their initial radius, R i , 0 , is described by Equation (1) and represents a dynamic equilibrium between the pressure forces and the liquid surface tension ( σ ):
R i , 0 = 2 σ P s a t ( T 0 ) P
where P is the ambient pressure computed in the Eulerian framework and P s a t is the liquid saturation pressure. T o is the liquid droplet temperature. If the vapor bubbles have a radius smaller than this value, they collapse under the force of the surface tension. R i , 0 is called the critical radius because any variation of the pressure difference “ P s a t ( T 0 ) P ” will cause the bubble to grow or collapse. According to the Rayleigh–Plesset equation and neglecting cooling effects from evaporation and thermal diffusion to the bubble surface (inertial-dominated growth), the bubble growth rate was calculated as follows.
V i = 3 ρ l p s a t p
This assumption can be accepted because for a bubble growing in a superheated droplet, the timescales from initiation to droplet breakup are in the order of tens of microseconds and so faster than vaporization and heat exchange phenomena [33] The droplet growth rate is related to the bubble growth rate as follows [29]:
V o = V i · Δ 2
The breakup of the bubble + droplet system is driven by instabilities that grow along the droplet and bubble surfaces represented in Figure 2a. The instability growth rate ω is the largest real root of the following normalized growth rate equation [29]:
Δ Δ 2 ψ o Δ Ω 2 + 1 + Δ 4 + ψ o W e o 1 / 2 Ω + 2 Δ 2 + 2 Δ 2 3 ψ i W e i M a i 2 Ω Ω + 3 W e i 1 / 2 Δ 2 = 0
where:
Ω = ρ R i 3 σ ω , W e o = ρ V o 2 R i σ , W e i = ρ V i 2 R i σ M a i = V i c Δ = R o R i , ψ o = ρ g o ρ , ψ i = ρ g i ρ
V i and V o are the bubble and droplet growth rate previously discussed, while ρ g i and ρ g o are the gas densities within and outside the bubble + droplet system. The following relationship describes the temporal evolution of this instability, denoted with η , which grows with the rate ω :
ln η η 0 = 0 t ω d t
When the instability amplitude becomes larger than a characteristic length of the spray (Figure 2b) [34], which can be chosen equal to the film thickness (i.e., the difference between droplet and bubble radius R o ( t ) R i ( t ) ), breakup takes place. It follows that the breakup criterion is equal to:
k b ( t ) = η 0 e 0 t ω d t R o ( t ) R i ( t )
where the initial disturbance η 0 is: η 0 = 0.05 · R o t 0 . Breakup occurs at t = t b when the instability amplitude equals the film thickness k b t b = k ¯ b = 1 . k ¯ b is a model constant that acknowledges the complexities that can be found in reality. In the present investigation, it is assumed equal to the unity. The flash boiling breakup is a catastrophic breakup because it leads to an explosion of the droplet, the expansions of the plumes and thus the plume-to-plume interaction. To reproduce this behavior, once the breakup occurs, a radial velocity component is added to the original velocity of the droplet as shown in Figure 3.
The magnitude of the radial velocity U r is obtained from mass and momentum conservation before and after breakup [29]:
U r = k v 3 R i 2 V i R o R i R o 3 R i 3
If follows that the new velocity components, after breakup, are equal to:
U x = U x , o l d + U x r U y = U y , o l d + U y r U z = U z , o l d + U z r
U r vector lies on a plane perpendicular to the initial velocity vector of the parcel and it has a random orientation. k v is a coefficient for scaling the radial velocity. Here, it is set equal to the unity. The size of droplets after breakup can be calculated with following relation obtained from the conservation of mass and energy [29].
S M R 32 1 = 1 2 R 0 2 + R i 2 R o 3 R i 3 + 3 2 R i 4 R i 1 R o 1 R o 3 R i 3 V i 2 U r 2 2 ρ 3 σ
where S M R 32 is the Sauter Mean Radius.
Mass conservation during the breakup process is obtained updating the number of particles is each parcel:
n n e w d n e w 3 = n o l d d o l d 3
where n o l d and n n e w are the number of particles, respectively, before and after breakup.
Each parcel can experience the flash-boiling breakup only one time. After that, it undergoes on aerodynamic breakup. Differently from the flash-boiling breakup, the latter can be applied multiple times. A specific input parameter allow to activate the fbBreakup model when flashing conditions were simulated. This combined thermodynamic–aerodynamic breakup model was validated in previous works of the authors, but in the present contribution, it will be extended to ultrahigh injection pressure conditions [23,28].
Further than modifying the breakup mechanism, flash boiling strongly influences the phase change process. Amsden’s energy balance was used to simulate heat transfer within droplets. The Ranz–Marshall correlation is applied to determine the rate of heat transfer per unit area. Full details and complete descriptions can be found in [30,35,36].
Price models was adopted for describing superheated liquid phase transition [37]. In such a model, the evaporation rate is given by two components: the subcooled and superheated terms. The subcooled term is computed as:
d M S c d t = 2 π R P Sh D i T f R f ln P P v P P S
where:
  • R: drop radius;
  • P : ambient pressure;
  • S h : nondimensional Sherwood number, function of Reynolds and Schmidt numbers:
    S h = 2 + 0.6 R e 0.5 S c 0.33
  • D i : binary diffusivity coefficient;
  • T f : vapor film temperature (average of the temperatures of the droplet and of the surrounding gas);
  • R f : vapor film specific gas constant;
  • P v : partial vapor pressure of the liquid species;
  • P s a t : saturation pressure.
The superheated term is:
d M s h d t = 4 π R 2 α T d T b H L
where:
  • H L : latent heat of vaporization of the liquid;
  • T b : boiling temperature of the fuel at the specific ambient pressure;
  • T d : droplet temperature;
  • α : heat transfer empirical coefficient [31].
It follows that the total evaporation rate is:
d M t d t = d M S C d t + d M s h d t
The other physical submodels, present in the Lagrangian simulation, are reported in Table 3.
Figure 4 shows the mass flow rate profiles measured by the injector manufacturer.
The parcel count for a grid of 0.125 mm is equal 200,000 parcels per nozzle accordingly with [30,38]. Parcels were initialized with a diameter equal to the nozzle one as reported in Table 1. The cone angle of the spray was measured from experimental images and set equal to 16 ° and 24 ° for, respectively, cases HP/vHP and cases HP-fb/vHPfb. Table 4 reports the parameters for the aerodynamic breakup. The Reitz–Diwakar’s constants C b a g , C b , C s t r i p and C s are set equal for both flash boiling and ambient conditions.
The injection environment is a cylinder with a height of 0.75 m and a diameter of 0.150 m, which matches the constant volume chamber used in the experiments. The domain was discretized using the internal CONVERGE mesh tool [30,39] with a 2 mm base grid, a fixed embedding refinement zone in the near nozzle area and an adaptive mesh refinement (AMR) depending on the velocity gradient for the spray development. Three different grid levels were tested: coarse, base and fine, characterized by different embedding and AMR levels. Table 5 reports all the details.
The injection environment is filled with nitrogen and RANS RNG-k- ϵ model used for turbulence modeling [21,32]. Numerical results were postprocessed using Paraview software.

3. Discussion of the Results

Experimental images were used to validate qualitatively and quantitatively the CFD simulations. The performances of the effervescent breakup model were assessed. Additional information not available from the shadowgraph images, such as droplet velocity and size, were also investigated to provide insights into spray structure.

3.1. Validation with Shadowgraph Acquisitions

Figure 5 reports, for case vHP, a comparison of the CFD spray shape (on the left) against the shadowgraph images (on the right) at different temporal instants after the start of injection (SOI).
The CFD model properly reproduce the V-shaped plumes experimentally observed and the evolution of the spray is correctly reproduced.
Figure 6 shows the same comparison but for the case vHP-fb.
The single plumes of the spray are clearly recognizable from the case vHP imsages (Figure 5), while, when the liquid is injected in superheated conditions (case vHP-fb-Figure 6), the plumes interact with each other, making it almost impossible to distinguish their profiles and appearing as a single large jet. This is an effect of the strong effervescent breakup, which promotes the interaction between the single plumes spreading liquid droplets in every direction and entailing greater compactness of the spray. It is interesting to observe that the spray collapse is almost completely avoided, despite being the superheat degree really high. This is a consequence of the high spray momentum owned by the droplet considering the 700 bar of injection pressure. As reported in Figure 7, the agreement between experiments and simulations is quite good also for the cases HP and HP-fb characterized by a lower injection pressure (equal to 400 bar). The spray shapes obtained from CFD simulations are compared against the shadowgraph images taken at 0.43   m s times after the start of injection (SOI).
A quantitative validation can performed comparing the spray penetration measured from the experiments and the one computed from the simulations.
Figure 8 and Figure 9 depicts a comparisons of experimental and numerical liquid penetration for all the simulated cases.
The penetration length was assumed as the maximum axial distance of any plume, measured from the injector tip, in the direction of the injector axis. The numerical penetration was computed as the axial distance from the injector where 99% of the liquid mass. The experiments were repeated five times to take into account for cycle-to-cycle variations. The uncertainty bars in the graphs report the variance from the average value. The simulations performed with the three different grid levels were used to assess the uncertainty of the numerical results. The criterion illustrated by Roache in [40], where the assessment procedures adopted by the AIAA, ITTC and IEEE are described and discussed in detail, was adopted for the purpose. The penetrations measured from the experimental images and those computed from the simulations are in quite good agreement, and the differences are almost always limited within the uncertainty bars. Only few differences can be found in the initial stages of the injections. These can be related to cavitation or flash-boiling phenomena that could take place within the injector. The jets of cases vHP and vHP-fb penetrate faster if compared to the ones of cases HP and HP-fb, respectively. The higher injection pressure gives an higher momentum to the liquid droplets. It is important to underline that, with this modeling approach, the same breakup constant can be used for simulating both non-evaporative and flash-boiling conditions.

3.2. Insights about the Flash Boiling Breakup Model

This section aims to assess the benefits of using the effervescent breakup model in the simulation of superheated sprays. The droplets’ size and the vaporization process were investigated. A further simulation, concerning case HP-fb, was performed adopting the Reitz-Diwakar model alone deactivating the fbBreakup model. The penetration lengths computed using Reitz-Diwakar and fbBreakup+Reitz-Diwakar breakup models are compared with the one experimentally measured in Figure 10.
The graph highlights that penetration length is overestimated using only the aerodynamic breakup model. This result has been already observed by the authors studying the Spray G with the same approach [23].
Figure 11 depicts the spray morphology at 0.4375 m s after the SOI; the droplets are represented, and they are colored accordingly with their diameter.
The effervescent breakup model predicts droplets approximately one order of magnitude smaller than what computed using only the Reitz–Diwakar model. Figure 12 shows the droplet size distributions obtained with fbBreakup+Reitz–Diwakar model and with only the Reitz–Diwakar model.
The parcels obtained from the fbBreakup+Reitz–Diwakar model have an average value equal to approximately 10 μ m while using the Reitz–Diwakar model parcel’s diameter ranges between 40 μ m – 80 μ m without a pseudo-average value. The different droplet size distribution influences also fuel vaporization. Figure 13 reports the vapor fuel mass in the computational domain as a function of the time elapsed from the start of the injection event.
The fbBreakup model promotes fuel vaporization indeed, after 750 μ s , the fuel evaporated is 20% greater.

3.3. Details about the Spray Morphology

After assessing the accuracy of CFD results with experiments, as well as the role of the effervescent breakup model, the simulation results were used to shed more light on other spray characteristics. Figure 14 is a snapshot taken at 0.4375   m s after the where spray’s droplets are represented and colored proportionately to their diameter for case vHP (top) and for case vHP-fb (bottom).
Flash boiling promotes both breakup and evaporation radically changing the morphology of the spray. The plumes of the flashing spray exhibit the interaction typical of a flash-boiling spray and are composed by smaller droplets. Figure 15 shows a volumetric rendering of the fuel mass fraction ( Y I C 8 H 18 ) evaporated after 0.4375   m s from the SOI for case vHP.
The amount of fuel vaporized is negligible because the injection takes place in non-evaporative conditions. More interesting is the case vHP-fb. Figure 16 reports the temporal evolution of the evaporated fuel for case vHP-fb.
The vaporization process regards only the core of the single plumes because the ambient temperature is lower (only 293 K) than the fuel one. The images indicate a not-severe plume-to-plume interaction but intense enough to entail great compactness of the spray, also in terms of vaporized fuel. Figure 17 shows the droplet’s Sauter mean Diameter ( S M D 32 ) as a function of the axial distance for case vHP and case vHP-fb.
The effervescent breakup strongly influences the morphology of the spray issued in the flash boiling condition. At a distance between 5 and 10 mm, the effervescent breakup reduces of approximately one order of magnitude the droplet’s dimension. The decrease of the SMD is much slighter in non-evaporative conditions. This is expected because the aerodynamic breakup is weaker if compared to the thermodynamic breakup mechanism.
Figure 18 shows the droplet’s Sauter mean Diameter ( S M D 32 ) vs the axial distance for the two injection pressures tested in flashing conditions.
The higher the injection pressure, the more downstream flash boiling breakup happens. This is a consequence of the higher injection pressure which gives to the jet an higher momentum leading to an extension of the flash-boiling breakup length. Furthermore, the higher injection pressure has two positive effects on the spray morphology. Figure 19 reports the droplet’s size distribution for the two injection pressures investigated.
An increase in such a parameter brings in having a more atomized spray, as demonstrated by the plot where the mean value of the distribution decrease as the injection pressure increases. The high injection pressure permits also to avoid an important undesired effect of flash-boiling: the spray collapse. Figure 20 is a representation of the spray where the interaction between air and fuel was highlighted over-imposing the velocity vectors on the vaporized fuel mass fraction. Case HP-fb is reported on the left and case HP-fb is reported on the right.
The velocity vectors demonstrate how the air flows from the bottom of the spray toward the injector’s nozzle avoiding the pressure reduction in the central region of the spray and promoting the mixture formation. This flow pattern inhibits also the spray collapse, which is one of the major causes of spray-wall impingement.
Further than improving spray breakup and vaporization, flash-boiling promotes also the turbulence generation and, consequently, air-fuel mixing [8]. Here, to evaluate the self-generated turbulence, vortexes were represented as Q-criterion isosurfaces and colored accordingly with total turbulent kinetic energy (TKE). Figure 21 reports such elaboration for cases vHP and vHP-fb.
The vortices spin around the axis of the single plume; their development is critical for air entrainment in the plume’s core to support the mixture formation. It is important to highlight also how moving from non-evaporative to flash-boiling conditions, the vortices dimension as well as TKE (turbulent kinetic energy) severally increase.

4. Conclusions

The analysis of the injection process in flash boiling conditions of a spray issued from an innovative multi-hole operating at of ultrahigh injection pressure was carried out. Experimental images of the spray acquired using shadowgraph technique were exploited to validate Eulerian–Lagrangian simulations of nonevaporative and flash-boiling conditions. The flash-boiling conditions considered were simulated using an effervescent breakup model, based on vapor bubbles growing within the liquid droplets, together with Reitz–Diwakar aerodynamic breakup model.
The main outcomes can be summarized as follows:
  • Subcooled and flash-boiling sprays can be reproduced with no changes the breakup constants by using an effervescent breakup model when superheated conditions are simulated.
  • Compared to an aerodynamic breakup model alone, the effervescent breakup model predicts droplets with smaller dimensions. The size distribution resemble a Gaussian distribution while, using Reitz–Diwakar model, the droplets’ diameters range from 40 μ m to 80 μ m .
  • The Reitz–Diwakar breakup model alone overestimates spray penetration when the constants used for nonevaporative conditions are used for flash-boiling conditions.
  • The effervescent breakup reduces the droplet diameter of about one order of magnitude (catastrophic breakup) at a distance which, depending on the injection pressure, ranges between 5 and 10 mm from the injector tip.
  • For the investigated cases, plume-to-plume interaction is significantly promoted by flash-boiling occurrence, but it is not strong enough to cause spray collapse.
  • Adopting an elevated values of injection pressures spray-collapse is avoided, and so air/fuel interaction is promoted as well as ambient gas recirculation.
  • Switching from nonevaporative to flash-boiling conditions both the turbulent kinetic energy and the dimensions of the vortical structures increase.
Future works will regard coupling of Eulerian near-nozzle and Eulerian–Lagrangian simulations to initialize the spray simulation properly and the investigation of even higher injection pressures.

Author Contributions

F.D.: Conceptualisation, methodology, software, data curation, validation, visualization, writing—original draft. A.D.V.: Writing—review and editing, supervision, project administration. A.M.: Investigation, resources, writing—review and editing. L.A.: Investigation, resources, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Acknowledgments

All the simulations were performed with the developed solver on the Galileo100 cluster of HPC CINECA facilities within the agreement between DIIIE—Università degli Studi dell’Aquila and CINECA.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
A M R Adaptive Mesh Refinement
C F D Computational Fluid Dynamic
D D M Discrete Droplet Model
E C N Engine Combustion Network
G D I Gasoline Direct Injection
H R M Homogeneous Relaxation Model
k b Breakup criterion
k v Constant for spray expansion
K H R T Kelvin–Helmholtz–Rayleigh–Taylor
M a i Bubble Mach number
nNumber of particles
P Ambient pressure
p a m b Nominal ambient pressure
P s a t Liquid saturation pressure
P I S O Pressure Implicit with splitting of operator
P I V Particle Image Velocimetry
R i , 0 Initial vapor bubble radius
R 0 Droplet radius
R i Vapor bubble radius
R p Ambient to saturation pressure ratio
R A N S Reynolds-Averaged Navier–Stokes
R D Reitz–Diwakar
S M D 32 Sauter Mean Diameter
S M R 32 Sauter Mean Radius
S O I Start Of Injection
S V C O Stepped Hole Valve Covered Orifice
T s a t Saturation temperature
T 0 Liquid temperature
t b Bubble breakup time
U r Magnitude of radial velocity
V i Bubble grow rate
V o Droplet grow rate
V O F Volume Of Fluid
W e i Bubble Weber number
W e o Droplet Weber number
η Disturbance
η 0 Initial disturbance
μ Dynamic viscosity
ω Instability grow rate
ψ o Air–Liquid density ratio
ψ i Bubble–Liquid density ratio
ρ g i Inner gas density
ρ g o Outer gas density
ρ l Liquid density
σ Liquid surface tension

References

  1. Payri, R.; Salvador, F.J.; Martí-Aldaraví, P.; Vaquerizo, D. ECN Spray G external spray visualization and spray collapse description through penetration and morphology analysis. Appl. Therm. Eng. 2017, 112, 304–316. [Google Scholar] [CrossRef]
  2. Di Ilio, G.; Krastev, V.K.; Falcucci, G. Evaluation of a Scale-Resolving Methodology for the Multidimensional Simulation of GDI Sprays. Energies 2019, 12, 2699. [Google Scholar] [CrossRef]
  3. Duronio, F.; De Vita, A.; Allocca, L.; Anatone, M. Gasoline direct injection engines—A review of latest technologies and trends. Part 1: Spray breakup process. Fuel 2020, 265, 116948. [Google Scholar] [CrossRef]
  4. Nocivelli, L.; Sforzo, B.A.; Tekawade, A.; Yan, J.; Powell, C.F.; Chang, W.; Lee, C.F.; Som, S. Analysis of the Spray Numerical Injection Modeling for Gasoline Applications; SAE Technical Papers; SAE International: Warrendale, PA, USA, 2020. [Google Scholar]
  5. Zembi, J.; Mariani, F.; Battistoni, M.; Irimescu, A.; Merola, S. Numerical Investigation of Water Injection Effects on Flame Wrinkling and Combustion Development in a GDI Spark Ignition Optical Engine; SAE WCX Digital Summit; SAE International: Warrendale, PA, USA, 2021. [Google Scholar]
  6. Oh, H.; Hwang, J.; Pickett, L.M.; Han, D. Machine-learning based prediction of injection rate and solenoid voltage characteristics in GDI injectors. Fuel 2022, 311, 122569. [Google Scholar] [CrossRef]
  7. Duronio, F.; De Vita, A.; Montanaro, A.; Villante, C. Gasoline direct injection engines—A review of latest technologies and trends. Part 2. Fuel 2020, 265, 116947. [Google Scholar] [CrossRef]
  8. Youso, T.; Fujikawa, T.; Yamakawa, M.; Kaminaga, T.; Yamaguchi, K.; Ratnak, S.; Kusaka, J. A Study on Combustion Characteristics of a High Compression Ratio SI Engine with High Pressure Gasoline Injection. In Proceedings of the 14th International Conference on Engines & Vehicles, Beijing, China, 28 July–1 August 2019; SAE International: Warrendale, PA, USA, 2019. [Google Scholar]
  9. Karrholm, F.P.; Helmantel, A.; Koopmans, L.; Dahlander, P.; Yamaguchi, A. Spray Characterization of Gasoline Direct Injection Sprays Under Fuel Injection Pressures up to 150 MPa with Different Nozzle Geometries. In Proceedings of the International Powertrains, Fuels & Lubricants Meeting, Kyoto, Japan, 26–29 August 2019; SAE International: Warrendale, PA, USA, 2019. [Google Scholar]
  10. Fatouraie, M.; Medina, M.; Wooldridge, M. High-Speed Imaging Studies of Gasoline Fuel Sprays at Fuel Injection Pressures from 300 to 1500 bar. In Proceedings of the WCX World Congress Experience, Detroit, MI, USA, 10–12 April 2018; SAE International: Warrendale, PA, USA, 2018. [Google Scholar]
  11. Guo, H.; Li, Y.; Lu, X.; Zhou, Z.; Xu, H.; Wang, Z. Radial expansion of flash boiling jet and its relationship with spray collapse in gasoline direct injection engine. Appl. Therm. Eng. 2019, 146, 515–525. [Google Scholar] [CrossRef]
  12. Xu, M.; Zhang, Y.; Zeng, W.; Zhang, G.; Zhang, M. Flash Boiling: Easy and Better Way to Generate Ideal Sprays than the High Injection Pressure. SAE Int. J. Fuels Lubr. 2013, 6, 137–148. [Google Scholar] [CrossRef]
  13. Sher, E.; Bar-Kohany, T.; Rashkovan, A. Flash-boiling atomization. Prog. Energy Combust. Sci. 2008, 34, 417–439. [Google Scholar] [CrossRef]
  14. Bar-Kohany, T.; Sher, E. Flash boiling atomization under negative pressure conditions. At. Sprays 2021, 31, 1–8. [Google Scholar] [CrossRef]
  15. Lefebvre, A.H.; McDonell, V.G. Atomization and Sprays; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  16. Brennen, E. Cavitation and Bubble Dynamics Christopher; California Institute of Technology Pasadena: Pasadena, CA, USA, 1995. [Google Scholar]
  17. Wang, B.; Mosbach, S.; Schmutzhard, S.; Shuai, S.; Huang, Y.; Kraft, M. Modelling soot formation from wall films in a gasoline direct injection engine using a detailed population balance model. Appl. Energy 2016, 163, 154–166. [Google Scholar] [CrossRef]
  18. Zeng, W.; Xu, M.; Zhang, G.; Zhang, Y.; Cleary, D.J. Atomization and vaporization for flash-boiling multi-hole sprays with alcohol fuels. Fuel 2012, 95, 287–297. [Google Scholar] [CrossRef]
  19. Chan, Q.N.; Bao, Y.; Kook, S. Effects of injection pressure on the structural transformation of flash-boiling sprays of gasoline and ethanol in a spark-ignition direct-injection (SIDI) engine. Fuel 2014, 130, 228–240. [Google Scholar] [CrossRef]
  20. Allocca, L.; Montanaro, A.; Meccariello, G. Effects of the Ambient Conditions on the Spray Structure and Evaporation of the ECN Spray G. In Proceedings of the WCX SAE World Congress Experience, Detroit, MI, USA, 9–11 April 2019; SAE International: Warrendale, PA, USA, 2019. [Google Scholar]
  21. Migliaccio, M.; Montanaro, A.; Paredi, D.; Lucchini, T.; Allocca, L.; D’Errico, G. CFD Modeling and Validation of the ECN Spray G Experiment under a Wide Range of Operating Conditions; SAE Technical Papers; SAE International: Warrendale, PA, USA, 2019. [Google Scholar]
  22. Duronio, F.; De Vita, A.; Allocca, L.; Montanaro, A.; Ranieri, S.; Villante, C. CFD Numerical Reconstruction of the Flash Boiling Gasoline Spray Morphology; SAE Technical Papers; SAE International: Warrendale, PA, USA, 2020. [Google Scholar]
  23. Duronio, F.; Mascio, A.D.; Villante, C.; Anatone, M.; Vita, A.D. ECN Spray G: Coupled Eulerian internal nozzle flow and Lagrangian spray simulation in flash boiling conditions. Int. J. Engine Res. 2023, 24, 14680874221090732. [Google Scholar] [CrossRef]
  24. Duronio, F.; Montanaro, A.; Ranieri, S.; Allocca, L.; De Vita, A. Under-Expanded Jets Characterization by Means of CFD Numerical Simulation Using an Open FOAM Density-Based Solver. In Proceedings of the 15th International Conference on Engines & Vehicles, Napoli, Italy, 12–16 September 2021; SAE International: Warrendale, PA, USA, 2021. [Google Scholar]
  25. Montanaro, A.; Allocca, L.; Ranieri, S.; Beatrice, C. Computational-Experimental Framework for Realizing a Novel Apparatus for Supercritical Water by Induction Heating. Heat Transf. Eng. 2022, 44, 785–802. [Google Scholar] [CrossRef]
  26. Duronio, F.; Montanaro, A.; Allocca, L.; Ranieri, S.; De Vita, A. Effects of Thermodynamic Conditions and Nozzle Geometry on the Methane Direct Injection Process in Internal Combustion Engines. In Proceedings of the WCX SAE World Congress Experience, Detroit, MI, USA, 5–7 April 2022; SAE International: Warrendale, PA, USA, 2022. [Google Scholar]
  27. Stiesch, G. Modeling Engine Spray and Combustion Processes; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  28. Duronio, F.; Ranieri, S.; Montanaro, A.; Allocca, L.; De Vita, A. ECN Spray G injector: Numerical modelling of flash-boiling breakup and spray collapse. Int. J. Multiph. Flow 2021, 145, 103817. [Google Scholar] [CrossRef]
  29. Zeng, Y. Modelling of Multicomponent Fuel Vaporization in Internal Combustion Engines. Ph.D. Thesis, University of Illinois at Urbana-Champaign, Champaign, IL, USA, 2000. [Google Scholar]
  30. Richards, K.; Senecal, P.; Pomraning, E. Converge 3.0; Convergent Science: Madison, WI, USA, 2021. [Google Scholar]
  31. Price, C.; Hamzehloo, A.; Aleiferis, P.; Richardson, D. Numerical modelling of fuel spray formation and collapse from multi-hole injectors under flash-boiling conditions. Fuel 2018, 22, 518–541. [Google Scholar] [CrossRef]
  32. Paredi, D.; Lucchini, T.; D’Errico, G.; Onorati, A.; Pickett, L.; Lacey, J. Validation of a comprehensive computational fluid dynamics methodology to predict the direct injection process of gasoline sprays using Spray G experimental data. Int. J. Engine Res. 2020. [Google Scholar] [CrossRef]
  33. Bilicki, Z.; Kestin, J.; Stuart, J.T. Physical aspects of the relaxation model in two-phase flow. Proc. R. Soc. Lond. Math. Phys. Sci. 1990, 428, 379–397. [Google Scholar]
  34. Reitz, R. Modeling atomization processes in high-pressure vaporizing sprays. At. Spray Technol. 1987, 3, 309–337. [Google Scholar]
  35. Faeth, G. Current status of droplet and liquid combustion. Prog. Energy Combust. Sci. 1977, 3, 191–224. [Google Scholar] [CrossRef]
  36. Amsden, A.A.; O’Rourke, P.J.; Butler, T.D. KIVA-II: A Computer Program for Chemically Reactive Flows with Sprays; Los Alamos National Lab.: Los Alamos, NM, USA, 1989. [Google Scholar]
  37. Price, C.; Hamzehloo, A.; Aleiferis, P.; Richardson, D. An approach to modeling flash-boiling fuel sprays for direct-injection spark-ignition engines. At. Sprays 2016, 26, 1197–1239. [Google Scholar] [CrossRef]
  38. Senecal, P.K.; Pomraning, E.; Richards, K.J.; Som, S. Grid-Convergent Spray Models for Internal Combustion Engine Computational Fluid Dynamics Simulations. J. Energy Resour. Technol. 2013, 136, 012204. [Google Scholar] [CrossRef]
  39. Di Angelo, L.; Duronio, F.; De Vita, A.; Di Mascio, A. Cartesian Mesh Generation with Local Refinement for Immersed Boundary Approaches. J. Mar. Sci. Eng. 2021, 9, 572. [Google Scholar] [CrossRef]
  40. Roache, P.J. Quantification of uncertainty in computational fluid dynamics. Annu. Rev. Fluid Mech. 1997, 29, 123–160. [Google Scholar] [CrossRef]
Figure 1. Scheme of the experimental apparatus highlighting the optical shadowgraph configuration.
Figure 1. Scheme of the experimental apparatus highlighting the optical shadowgraph configuration.
Fluids 08 00155 g001
Figure 2. Time evolution of bubble–droplet system. (a) The breakup of the bubble + droplet system is driven by instabilities that grow along the droplet and bubble surfaces. (b) When the instability amplitude becomes larger than a characteristic length of the spray.
Figure 2. Time evolution of bubble–droplet system. (a) The breakup of the bubble + droplet system is driven by instabilities that grow along the droplet and bubble surfaces. (b) When the instability amplitude becomes larger than a characteristic length of the spray.
Fluids 08 00155 g002
Figure 3. Droplet velocity after breakup updating.
Figure 3. Droplet velocity after breakup updating.
Fluids 08 00155 g003
Figure 4. Mass flow-rate profiles for the two injection pressures.
Figure 4. Mass flow-rate profiles for the two injection pressures.
Fluids 08 00155 g004
Figure 5. Comparison between numerical and experimental results for case vHP.
Figure 5. Comparison between numerical and experimental results for case vHP.
Fluids 08 00155 g005
Figure 6. Comparison between numerical and experimental results for case vHP-fb.
Figure 6. Comparison between numerical and experimental results for case vHP-fb.
Fluids 08 00155 g006
Figure 7. CFD results and experimental shadowgraph images for cases HP (a) and vHPfb (b). t = 0.43 m s after the SOI.
Figure 7. CFD results and experimental shadowgraph images for cases HP (a) and vHPfb (b). t = 0.43 m s after the SOI.
Fluids 08 00155 g007aFluids 08 00155 g007b
Figure 8. Numerical and experimental liquid penetration for case HP and vHP vs. time.
Figure 8. Numerical and experimental liquid penetration for case HP and vHP vs. time.
Fluids 08 00155 g008
Figure 9. Numerical and experimental liquid penetration for case vHP-fb and 2fb vs. time.
Figure 9. Numerical and experimental liquid penetration for case vHP-fb and 2fb vs. time.
Fluids 08 00155 g009
Figure 10. Comparison of the numerical penetration curves computed using Reitz–Diwakar and fbBreakup breakup models with experimental measurements for case HP-fb.
Figure 10. Comparison of the numerical penetration curves computed using Reitz–Diwakar and fbBreakup breakup models with experimental measurements for case HP-fb.
Fluids 08 00155 g010
Figure 11. Spray morphology represented coloring liquid droplets accordingly with their diameter. t = 0.4375 m s after the SOI. Case HP-fb simulated using fbBreakup+Reitz–Diwakar model is reported on top while simulation with only the Reitz–Diwakar model on bottom.
Figure 11. Spray morphology represented coloring liquid droplets accordingly with their diameter. t = 0.4375 m s after the SOI. Case HP-fb simulated using fbBreakup+Reitz–Diwakar model is reported on top while simulation with only the Reitz–Diwakar model on bottom.
Fluids 08 00155 g011
Figure 12. Predicted droplet size distributions for case HP-fb at 0.4 m s after SOI.
Figure 12. Predicted droplet size distributions for case HP-fb at 0.4 m s after SOI.
Fluids 08 00155 g012
Figure 13. Vaporized fuel mass as function of elapsed time from the beginning of the injection.
Figure 13. Vaporized fuel mass as function of elapsed time from the beginning of the injection.
Fluids 08 00155 g013
Figure 14. Spray morphology for the cases vHP-fb (top) and vHP (bottom) represented with droplets colored accordingly with their diameter at 0.4375 m s after the SOI.
Figure 14. Spray morphology for the cases vHP-fb (top) and vHP (bottom) represented with droplets colored accordingly with their diameter at 0.4375 m s after the SOI.
Fluids 08 00155 g014
Figure 15. 3D visualization of the concentration of the vaporized fuel.
Figure 15. 3D visualization of the concentration of the vaporized fuel.
Fluids 08 00155 g015
Figure 16. 3D visualization of the concentration of the vaporized fuel.
Figure 16. 3D visualization of the concentration of the vaporized fuel.
Fluids 08 00155 g016
Figure 17. Average SMD as a function of the axial distance from injector tip for cases vHP and case vHP-fb.
Figure 17. Average SMD as a function of the axial distance from injector tip for cases vHP and case vHP-fb.
Fluids 08 00155 g017
Figure 18. Average SMD as a function of the axial distance from injector tip for cases vHP and case vHP-fb.
Figure 18. Average SMD as a function of the axial distance from injector tip for cases vHP and case vHP-fb.
Fluids 08 00155 g018
Figure 19. Predicted droplet size distributions for case HP and vHP at 0.43 ms after SOI.
Figure 19. Predicted droplet size distributions for case HP and vHP at 0.43 ms after SOI.
Fluids 08 00155 g019
Figure 20. Velocity vectors and vapor mass fraction at t S O I = 0.43 ms for cases HP-fb (a) and vHPfb (b).
Figure 20. Velocity vectors and vapor mass fraction at t S O I = 0.43 ms for cases HP-fb (a) and vHPfb (b).
Fluids 08 00155 g020
Figure 21. Q-criterion isosurfaces colored by total TKE for case vHP (a) and vHP-fb (b).
Figure 21. Q-criterion isosurfaces colored by total TKE for case vHP (a) and vHP-fb (b).
Fluids 08 00155 g021
Table 1. Principal features of the GDI injector.
Table 1. Principal features of the GDI injector.
Spray shapeCircular
Nozzle typeVCO (Valve-covered orifice)
Nozzle shapeStepped-hole
Orifice length700 μ m
Orifice diameter100 μ m
L/D ratio7
Orifice drill angle 35.5 ° / 45 °
Full outer angle90 °
Flow rate10.55 cc/s @ 100 bar
Discharge Coeffient0.737 @ 400 bar/0.937 @ 700 bar
Table 2. Investigated operating conditions.
Table 2. Investigated operating conditions.
Case HPCase HP-fbCase vHPCase vHP-fb
Injection Pressure [bar]400400700700
Fuel Temp. [K]293363293363
Ambient Pressure [bar]10.210.2
Ambient Temp. [K]293293293293
Table 3. Lagrangian framework physical submodels.
Table 3. Lagrangian framework physical submodels.
Injection modelFlow Rate + Discharge
Effervescent BreakupfbBreakup
Aereodynamic BreakupReitz–Diwakar
Heat transferRanz–Marshall
VaporizationFrossling + Flash Boiling Vaporization
CollisionNo Time Counter (NTC)
CoalescencePost and Abraham
DispersionO’Rourke
Table 4. Breakup model parameters.
Table 4. Breakup model parameters.
k b 1
k v 1
C bag 6
C b 1.5
C strip 1
C s 10
Table 5. Description of the different grids adopted for CFD simulations.
Table 5. Description of the different grids adopted for CFD simulations.
AMR and Fixed Embedding LevelGrid Size
Coarse2500 μ m
Base3250 μ m
Fine4125 μ m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Duronio, F.; De Vita, A.; Montanaro, A.; Allocca, L. Experimental Investigation and Numerical CFD Assessment of a Thermodynamic Breakup Model for Superheated Sprays with Injection Pressure up to 700 Bar. Fluids 2023, 8, 155. https://doi.org/10.3390/fluids8050155

AMA Style

Duronio F, De Vita A, Montanaro A, Allocca L. Experimental Investigation and Numerical CFD Assessment of a Thermodynamic Breakup Model for Superheated Sprays with Injection Pressure up to 700 Bar. Fluids. 2023; 8(5):155. https://doi.org/10.3390/fluids8050155

Chicago/Turabian Style

Duronio, Francesco, Angelo De Vita, Alessandro Montanaro, and Luigi Allocca. 2023. "Experimental Investigation and Numerical CFD Assessment of a Thermodynamic Breakup Model for Superheated Sprays with Injection Pressure up to 700 Bar" Fluids 8, no. 5: 155. https://doi.org/10.3390/fluids8050155

APA Style

Duronio, F., De Vita, A., Montanaro, A., & Allocca, L. (2023). Experimental Investigation and Numerical CFD Assessment of a Thermodynamic Breakup Model for Superheated Sprays with Injection Pressure up to 700 Bar. Fluids, 8(5), 155. https://doi.org/10.3390/fluids8050155

Article Metrics

Back to TopTop