An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model
Pith reviewed 2026-05-08 05:26 UTC · model grok-4.3
The pith
A four-equation multiphase model with added subgrid spray and finite-rate phase change predicts evaporating spray behavior accurately.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a robust four-equation multiphase model, when augmented with a new phase-confined Eulerian Σ spray equation and an improved finite-rate phase change model thermodynamically bounded by Gibbs-free energy equilibration, satisfies interface equilibrium and phase immiscibility conditions while delivering large computational savings and accurate predictions of subgrid surface area and evaporation rates.
What carries the argument
The phase-confined Eulerian Σ spray model coupled to a Gibbs-free-energy-bounded finite-rate phase change model inside the four-equation multiphase framework.
If this is right
- The equilibrium assumptions produce substantial computational savings that permit larger domains or finer grids in LES of sprays.
- The phase-confined Σ model supplies subgrid interfacial area while preventing artificial mass or momentum transfer across resolved interfaces.
- The bounded phase-change model captures finite-rate evaporation and condensation across a range of spray regimes.
- Validation on the ECN Spray A case confirms that the combined closures reproduce measured spray penetration and vaporization under both non-evaporating and evaporating conditions.
Where Pith is reading between the lines
- If the subgrid equilibrium premise proves robust, analogous closures could be inserted into other four-equation or diffuse-interface models used for boiling or cavitation.
- The framework suggests that targeted subgrid models for surface area and thermodynamic driving force can compensate for the loss of explicit nonequilibrium variables.
- Direct comparison against resolved-interface DNS in a regime with strong local temperature or velocity slip would quantify the range of validity of the equilibrium assumption.
Load-bearing premise
The four-equation multiphase model assumes strict subgrid equilibrium of pressure, temperature, and velocity, so the new Σ and phase-change models must restore predictive power under these restrictive conditions.
What would settle it
If the predicted liquid penetration length, vapor mass fraction, or droplet size distribution deviates systematically from experimental data in an evaporating spray test case where subgrid velocity or temperature nonequilibrium is known to be strong, the claim that the equilibrium-based framework remains accurate would be falsified.
Figures
read the original abstract
In this work, an LES model with finite-rate phase change and subgrid spray based on a high-resolution numerical scheme for multiphase multi-component simulations which satisfies interface equilibrium and phase immiscibility conditions is proposed. The multiphase model is based on a robust implementation of the four-equation multiphase model which assumes a strict subgrid equilibrium of pressure, temperature, and velocity. Critically, the equilibrium assumptions of the four-equation model provide large computational savings compared to modeling the full non-equilibrium multiphase system. To obtain predictive capabilities with these restrictive equilibrium assumptions, a new phase-confined form of the Eulerian $\Sigma$ spray model is proposed to predict subgrid interfacial surface area while avoiding unphysical leakage across interfaces. Additionally, an improved finite rate phase change model which is thermodynamically bounded by the equilibration of the Gibbs-free energy is coupled with the $\Sigma$ equation to model complex phase change regimes. The full modeling framework is validated using the Engine Combustion Network (ECN) Spray A case in non-evaporating and evaporating conditions and shows excellent agreement with experimental measurements.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes an LES framework for multiphase multi-component flows with phase change, built on a four-equation equilibrium model that assumes strict subgrid equilibrium of pressure, temperature, and velocity for computational efficiency. To restore predictive capability under these restrictive assumptions, the authors introduce a phase-confined Eulerian Σ spray model that avoids unphysical interfacial leakage and a finite-rate phase-change model thermodynamically bounded by Gibbs-free energy equilibration. The full framework is validated on the ECN Spray A case under both non-evaporating and evaporating conditions, with the abstract claiming excellent quantitative agreement with experimental measurements.
Significance. If the reported agreement with the ECN Spray A benchmark is robust, the work demonstrates a practical route to accurate LES of sprays and phase change at reduced cost relative to fully non-equilibrium multiphase solvers. The targeted closures for subgrid surface area and bounded phase change directly address the limitations of the equilibrium assumptions while preserving thermodynamic consistency, which is a notable strength for applications in combustion and fluid dynamics.
minor comments (2)
- The abstract states 'excellent agreement' but does not specify the quantitative metrics (e.g., integrated error norms, pointwise comparisons, or visual metrics) used to support this claim; adding this detail would improve clarity.
- The manuscript would benefit from explicit discussion of mesh convergence and sensitivity to the equilibrium assumptions in the validation section, given that these assumptions are described as restrictive.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of the manuscript and the recommendation for minor revision. We are encouraged that the significance of combining the four-equation equilibrium model with the phase-confined Eulerian Σ spray closure and Gibbs-free-energy-bounded phase change is recognized as a practical route to accurate yet efficient LES of sprays.
Circularity Check
No significant circularity detected
full rationale
The paper develops a new modeling framework by extending the four-equation multiphase model with a phase-confined Eulerian Σ spray closure and a thermodynamically bounded finite-rate phase-change model. These additions are introduced to compensate for the restrictive subgrid equilibrium assumptions, and the framework is then validated directly against independent external experimental measurements from the public ECN Spray A dataset under both non-evaporating and evaporating conditions. No load-bearing derivation step, equation, or prediction is shown to reduce by construction to a parameter fit, self-definition, or self-citation chain within the paper itself; the reported agreement therefore constitutes external evidence rather than tautological reproduction of inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Strict subgrid equilibrium of pressure, temperature, and velocity
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
-
[4]
A. Chiapolino, P. Boivin, R. Saurel, A simple and fast phase transition relaxation solver for compressible multicomponent two-phase flows, Computers & Fluids 150 (2017) 31–45
work page 2017
-
[5]
M. Pelanti, Arbitrary-rate relaxation techniques for the numerical modeling of compressible two-phase flows with heat and mass transfer, International Journal of Multiphase Flow 153 (2022) 104097
work page 2022
-
[6]
A. D. Demou, N. Scapin, M. Pelanti, L. Brandt, A pressure-based diffuse interface method for low-mach multiphase flows with mass transfer, Journal of Computational Physics 448 (2022) 110730
work page 2022
-
[7]
P. Downar-Zapolski, Z. Bilicki, L. Bolle, J. Franco, The non-equilibrium relaxation model for one- dimensional flashing liquid flow, International journal of multiphase flow 22 (3) (1996) 473–483
work page 1996
-
[8]
J. W. Gaertner, A. Kronenburg, A. Rees, J. Sender, M. Oschwald, G. Lamanna, Numerical and exper- imental analysis of flashing cryogenic nitrogen, International Journal of Multiphase Flow 130 (2020) 103360
work page 2020
-
[9]
K. G. Lyras, S. Dembele, J. X. Wen, Numerical simulation of flashing jets atomisation using a unified approach, International Journal of Multiphase Flow 113 (2019) 45–58
work page 2019
-
[10]
B. Abramzon, W. A. Sirignano, Droplet vaporization model for spray combustion calculations, Inter- national journal of heat and mass transfer 32 (9) (1989) 1605–1618
work page 1989
-
[11]
H. Hertz, Ueber die verdunstung der flüssigkeiten, insbesondere des quecksilbers, im luftleeren raume, Annalen der Physik 253 (10) (1882) 177–193
-
[12]
M. Knudsen, Die maximale verdampfungsgeschwindigkeit des quecksilbers, Annalen der Physik 352 (13) (1915) 697–708
work page 1915
-
[13]
S. S. Sazhin, Advanced models of fuel droplet heating and evaporation, Progress in energy and com- bustion science 32 (2) (2006) 162–214
work page 2006
- [14]
- [15]
-
[16]
J. K. Dukowicz, A particle-fluid numerical model for liquid sprays, Journal of computational Physics 35 (2) (1980) 229–253
work page 1980
-
[17]
R. Reitz, et al., Modeling atomization processes in high-pressure vaporizing sprays, Atomisation and Spray technology 3 (4) (1987) 309–337
work page 1987
-
[18]
B. M. Devassy, C. Habchi, E. Daniel, Atomization modelling of liquid jets using a two-surface-density approach, Atomization and Sprays 25 (1)
-
[19]
Q. Xue, S. Som, P. K. Senecal, E. Pomraning, Large eddy simulation of fuel-spray under non-reacting ic engine conditions, Atomization and Sprays 23 (10)
-
[20]
Q. Xue, M. Battistoni, C. Powell, D. Longman, S. Quan, E. Pomraning, P. Senecal, D. Schmidt, S. Som, An eulerian cfd model and x-ray radiography for coupled nozzle flow and spray in internal combustion engines, International Journal of Multiphase Flow 70 (2015) 77–88. 27
work page 2015
-
[21]
Desantes, J.M.Garcia-Oliver, J.M
J.M. Desantes, J.M.Garcia-Oliver, J.M. Pastor, A.Pandal, Acomparison of dieselsprayscfdmodeling approaches: Ddm versusσ-y eulerian atomization model, Atomization and Sprays 26 (7)
-
[22]
G. M. Magnotti, C. L. Genzale, Detailed assessment of diesel spray atomization models using visible and x-ray extinction measurements, International Journal of Multiphase Flow 97 (2017) 33–45
work page 2017
-
[23]
M.Jia, H.Pan, Y.Bian, Z.Zhang, Y.Chang, H.Liu, Calibrationoftheconstantsinthekelvin-helmholtz rayleigh-taylor (kh-rt) breakup model for diesel spray under wide conditions based on advanced data analysis techniques, Atomization and Sprays 32 (6)
-
[24]
G. Tryggvason, R. Scardovelli, S. Zaleski, Direct numerical simulations of gas–liquid multiphase flows, Cambridge university press, 2011
work page 2011
-
[25]
M. Herrmann, On simulating primary atomization using the refined level set grid method, Atomization and Sprays 21 (4)
- [26]
- [27]
-
[28]
A. Burluka, R. Borghi, et al., Development of a eulerian model for the “atomization” of a liquid jet, Atomization and sprays 11 (6)
- [29]
- [30]
-
[31]
J. Anez, A. Ahmed, N. Hecht, B. Duret, J. Reveillon, F. Demoulin, Eulerian–lagrangian spray atom- ization model coupled with interface capturing method for diesel injectors, International Journal of Multiphase Flow 113 (2019) 325–342
work page 2019
-
[32]
J. Desantes, J. M. García-Oliver, J. Pastor, A. Pandal, E. Baldwin, D. Schmidt, Coupled/decoupled spray simulation comparison of the ecn spray a condition with theσ-y eulerian atomization model, International Journal of Multiphase Flow 80 (2016) 89–99
work page 2016
-
[33]
A. Pandal, J. M. Pastor, R. Payri, A. Kastengren, D. Duke, K. Matusik, J. S. Giraldo, C. Powell, D. Schmidt, Computational and experimental investigation of interfacial area in near-field diesel spray simulation, SAE International Journal of Fuels and Lubricants 10 (2) (2017) 423–431
work page 2017
-
[34]
J. Desantes, J. M. García-Oliver, J. Pastor, I. Olmeda, A. Pandal, B. Naud, Les eulerian diffuse-interface modeling of fuel dense sprays near-and far-field, International Journal of Multiphase Flow 127 (2020) 103272
work page 2020
-
[35]
G. Nykteri, M. Gavaises, Droplet aerobreakup under the shear-induced entrainment regime using a multiscale two-fluid approach, Physical Review Fluids 6 (8) (2021) 084304
work page 2021
-
[36]
H. Gaballa, C. Habchi, J.-C. de Hemptinne, Modeling and les of high-pressure liquid injection un- der evaporating and non-evaporating conditions by a real fluid model and surface density approach, International Journal of Multiphase Flow 160 (2023) 104372. 28
work page 2023
- [37]
-
[38]
C. A. Kennedy, A. Gruber, Reduced aliasing formulations of the convective terms within the navier– stokes equations for a compressible fluid, Journal of Computational Physics 227 (3) (2008) 1676–1700
work page 2008
-
[39]
Y. Kuya, S. Kawai, High-order accurate kinetic-energy and entropy preserving (keep) schemes on curvi- linear grids, Journal of Computational Physics 442 (2021) 110482.doi:10.1016/J.JCP.2021.110482
-
[40]
S. S. Jain, P. Moin, A kinetic energy–and entropy-preserving scheme for compressible two-phase flows, Journal of Computational Physics 464 (2022) 111307.doi:10.1016/J.JCP.2022.111307
- [41]
-
[42]
H. Collis, A robust high order framework for compressible multi-phase multi-component flows with interface regularization, phase change, and spray modeling, Stanford University, 2025
work page 2025
-
[43]
H. Collis, S. Mirjalili, A. Mani, Diffuse interface treatment in generalized curvilinear coordinates with grid-adapting interface thickness, Journal of Computational Physics 544 (2026) 114440.doi:10.1016/ j.jcp.2025.114440. URLhttps://www.sciencedirect.com/science/article/pii/S0021999125007223
-
[44]
S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM review 43 (1) (2001) 89–112
work page 2001
-
[45]
K. W. Thompson, Time-dependent boundary conditions for hyperbolic systems, ii, Journal of compu- tational physics 89 (2) (1990) 439–461
work page 1990
-
[46]
T. J. Poinsot, S. Lelef, Boundary conditions for direct simulations of compressible viscous flows, Journal of computational physics 101 (1) (1992) 104–129
work page 1992
-
[47]
N. Okong’o, J. Bellan, Consistent boundary conditions for multicomponent real gas mixtures based on characteristic waves, Journal of Computational Physics 176 (2) (2002) 330–344
work page 2002
- [48]
-
[49]
M. R. Baer, J. W. Nunziato, A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials, International journal of multiphase flow 12 (6) (1986) 861–889
work page 1986
-
[50]
Linstrom, Nist chemistry webbook, nist standard reference database 69, (No Title) (1997) 20899
P. Linstrom, Nist chemistry webbook, nist standard reference database 69, (No Title) (1997) 20899
work page 1997
-
[51]
X. Deng, P. Boivin, Diffuse interface modelling of reactive multi-phase flows applied to a sub-critical cryogenic jet, Applied Mathematical Modelling 84 (2020) 405–424
work page 2020
-
[52]
J. Chesnel, J. Reveillon, T. Menard, F.-X. Demoulin, Large eddy simulation of liquid jet atomization, Atomization and Sprays 21 (9)
-
[53]
S. Mirjalili, M. Khanwale, A. Mani, Consistent modeling of scalar transport in multiphase flows using conservative phase field methods, Proceedings of the Summer Program 2022, Stanford University Center for Turbulence Research
work page 2022
- [54]
-
[55]
A. Yoshizawa, K. Horiuti, A statistically-derived subgrid-scale kinetic energy model for the large-eddy simulation of turbulent flows, Journal of the Physical Society of Japan 54 (8) (1985) 2834–2839. 29
work page 1985
-
[56]
URLhttps://www.openfoam.com/documentation/guides/v2206/api/classFoam_1_1LESModels_ 1_1WALE.html
OpenCFD Ltd., OpenFOAM v2206 API Guide: Foam::LESModels::WALE Class, [Online; accessed 14-August-2025] (2022). URLhttps://www.openfoam.com/documentation/guides/v2206/api/classFoam_1_1LESModels_ 1_1WALE.html
work page 2025
-
[57]
J. W. Gärtner, A. Kronenburg, A novel elsa model for flash evaporation, International Journal of Multiphase Flow 174 (2024) 104784
work page 2024
-
[58]
Prosperetti, Vapor bubbles, Annual review of fluid mechanics 49 (1) (2017) 221–248
A. Prosperetti, Vapor bubbles, Annual review of fluid mechanics 49 (1) (2017) 221–248
work page 2017
- [59]
- [60]
-
[61]
Engine Combustion Network (ECN), ECN,https://ecn.sandia.gov, an international collaboration among experimental and computational researchers in engine combustion, maintained by the Applied Combustion Department, Sandia National Laboratories (2025)
work page 2025
- [62]
-
[63]
Engine Combustion Network (ECN), Spray a geometry,https://ecn.sandia.gov/ diesel-spray-combustion/target-condition/spray-a-nozzle-geometry(2022)
work page 2022
-
[64]
CMT, Virtual injection rate generator,https://www.cmt.upv.es/#/ecn/download/ InjectionRateGenerator/InjectionRateGenerator
- [65]
-
[66]
A. Kastengren, F. Tilocco, D. Duke, C. Powell, S. Moon, X. Zhang, et al., Time-resolved x-ray radiog- raphy of diesel injectors from the engine combustion network, ICLASS Paper 1369
-
[67]
Engine Combustion Network (ECN), Near-nozzle mixture derived from x-ray radiography,https:// ecn.sandia.gov/rad675/(2022)
work page 2022
-
[68]
I. Karathanassis, P. Koukouvinis, M. Gavaises, Comparative evaluation of phase-change mechanisms for the prediction of flashing flows, International Journal of Multiphase Flow 95 (2017) 257–270
work page 2017
-
[69]
J. Manin, M. Bardi, L. M. Pickett, Sp2-4 evaluation of the liquid length via diffused back-illumination imaging in vaporizing diesel sprays (sp: Spray and spray combustion, general session papers), in: The Proceedings of the International symposium on diagnostics and modeling of combustion in internal combustion engines 2012.8, The Japan Society of Mechan...
work page 2012
-
[70]
Engine Combustion Network (ECN), Measurements of the liquid length for nozzle 210675 using diffused back-illumination (dbi) at spray a and other conditions,https://ecn.sandia.gov/data/dbi675/ (2022)
work page 2022
- [71]
-
[72]
Engine Combustion Network (ECN), Diesel data search,https://ecn.sandia.gov/ecn-data-search (2022). 30
work page 2022
-
[73]
L. M. Pickett, C. L. Genzale, J. Manin, Uncertainty quantification for liquid penetration of evaporating sprays at diesel-like conditions, Atomization and Sprays 25 (5)
-
[74]
gov/diesel-spray-combustion/computational-method/modeling-standards/(2022)
EngineCombustionNetwork(ECN),Modelingstandardsandrecommendations,https://ecn.sandia. gov/diesel-spray-combustion/computational-method/modeling-standards/(2022)
work page 2022
-
[75]
E. F. Toro, Riemann solvers and numerical methods for fluid dynamics: a practical introduction, Springer Science & Business Media, 2013
work page 2013
-
[76]
L. Fu, X. Y. Hu, N. A. Adams, A new class of adaptive high-order targeted eno schemes for hyperbolic conservation laws, Journal of Computational Physics 374 (2018) 724–751. 31
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.