Boiling Flow Algorithm Overview

• Boiling flow algorithms are computational methods that simulate multiphase liquid-vapor transitions by coupling nucleation, bubble dynamics, and convective heat transfer.
• They integrate techniques such as the lattice Boltzmann method, finite-element VOF/level-set approaches, and machine-learned surrogates to balance fidelity and computational cost.
• Recent advancements enable robust simulation of nucleate, transition, and film boiling regimes while ensuring thermodynamic consistency and accurate interfacial dynamics.

Boiling flow algorithms are computational methodologies for simulating the dynamic, multiphase, and thermally driven process of liquid-vapor phase change in the presence of convective flow. They integrate the fundamental physics of nucleation, bubble growth, departure, phase interfaces, latent heat transfer, and their coupling to fluid and thermal fields. The field encompasses kinetic-theory-based mesoscopic methods such as the lattice Boltzmann method (LBM), sharp and diffuse-interface finite-volume and finite-element models, one-fluid and two-fluid Navier–Stokes approaches, fully machine-learned surrogates, and spectral stochastic models—each providing different balances of fidelity, computational cost, and generalizability.

1. Mathematical Formulations and Physical Principles

Boiling flow algorithms numerically represent mass, momentum, and energy transport equations augmented by models of phase change and interface dynamics. Key formulations can be classified as follows:

• Mesoscopic Kinetic Theory Models: The pseudopotential lattice Boltzmann method (LBM), with D2Q9 or D3Q19 lattice stencils, is frequently employed. The collision-streaming equation

$$f_i(\mathbf{x}+\mathbf{c}_i\Delta t, t+\Delta t) = f_i(\mathbf{x},t) - \frac{1}{\tau}\left[ f_i(\mathbf{x},t) - f_i^{\mathrm{eq}}(\mathbf{x},t) \right] + \Delta t F_i(\mathbf{x},t)$$

is coupled to a fluid-fluid pseudopotential force, e.g., of Shan–Chen type, which induces spontaneous liquid-vapor phase separation and interface tension without explicit interface tracking (Li et al., 2020, Li et al., 2015, Mukherjee et al., 2020, Biferale et al., 2011, Saito et al., 2019).

• Energy Equation and Latent Heat: Energy transport is consistently modeled by including the latent heat effect as a source or flux term, as in

$$\rho c_v ( \partial_t T + \mathbf{u} \cdot \nabla T ) = \nabla \cdot (\lambda \nabla T) - T \left( \frac{\partial p_{EOS}}{\partial T} \right)_\rho \nabla \cdot \mathbf{u}$$

This framework automatically captures entropy generation and phase change at spontaneously evolving interfaces (Li et al., 2015, Li et al., 2020, Saito et al., 2019).

• Multiphase Interface Models: Approaches include:
  • Diffuse-interface: Smoothed or signed-distance level set or volume fraction fields characterize the interface, allowing for surface tension forces (CSF model), phase change localization, and mass transfer sources (Iskhakova et al., 2022, Poblador-Ibanez et al., 2024, Dhruv, 2023).
  • Sharp-interface (VOF/Level-Set): Interfacial jump conditions are enforced directly, often with specialized numerical stencils for mass, momentum, and energy exchange (Poblador-Ibanez et al., 2024).
• One-fluid vs. Two-fluid Formulations: The one-fluid model blends properties across phases and resolves interfaces via additional source terms or steep property transitions; the two-fluid approach treats liquid and vapor as interpenetrating continua with distinct velocity, temperature, and volume fraction fields, requiring closure laws for interfacial transfer of mass, momentum, and energy (e.g., RPI partitioning model) (Fu et al., 2017, Končar et al., 22 Sep 2025).

A prevailing theme is the necessity to recover the correct macroscopic jump conditions for velocity, pressure, and temperature at the interface, either implicitly via kinetic models or explicitly via local source terms and corrective forces (Poblador-Ibanez et al., 2024).

2. Algorithmic Strategies and Numerical Methods

Boiling flow algorithms are implemented in both explicit and implicit time integration schemes and across a range of discretizations:

• Central-Moment Lattice Boltzmann: Advanced schemes operate in central-moment (CM) space for robust high-Reynolds number stability and precise thermodynamic consistency. For example, collision and forcing are projected onto a non-orthogonal CM basis, with relaxation parameters tuned for conserved, shear, and higher moments (Saito et al., 2019).
• Pseudopotential-Driven Phase Change: In kinetic-theory models, phase change, nucleation, and bubble dynamics arise naturally from the non-ideal equation of state and the pseudopotential coupling without explicit source terms or nucleation-site "seeding" (Li et al., 2020, Mukherjee et al., 2020, Biferale et al., 2011, Saito et al., 2019, Li et al., 2015).
• Finite-Element and Volume-of-Fluid (VOF): High-fidelity boiling simulations with complex geometries utilize stabilized finite-element methods coupled with level-set or VOF for interface capturing and interface-conformal application of jump conditions. Operator-split and monolithic solvers are used for Navier–Stokes and energy equations, with phase change and surface tension entered as immersed or interface-localized sources (Iskhakova et al., 2022, Poblador-Ibanez et al., 2024).
• Heat Flux Partitioning Models: In two-fluid and industrial CFD codes (OpenFOAM, ANSYS CFX), the boiling wall heat flux is partitioned into convective, quenching, and evaporation terms. Closure relations for nucleation site density, bubble departure diameter, and frequency are required (e.g., Lemmert–Chawla, Tolubinsky–Kostanchuk, Cole models) (Fu et al., 2017, Končar et al., 22 Sep 2025).
• Hybrid and Machine-Learned Surrogates: Recent advances include data-driven surrogates such as transformer-based PDE predictors (Bubbleformer) that learn spatiotemporal evolution of interface, temperature, and velocity from high-fidelity simulation data, and latent diffusion models for mapping phase indicator fields to temperature distributions at a fraction of the cost of direct simulation (Hassan et al., 28 Jul 2025, Na et al., 27 Jan 2025).

3. Physical Regimes and Validation Metrics

Boiling flow algorithms are systematically validated against canonical regimes and metrics:

• Boiling Regimes: Accurate reproduction of nucleate, transition, and film boiling is a standard for credibility. Algorithms must recover the full boiling (Nukiyama) curve: critical heat flux (CHF), minimum heat flux (MHF), boiling crisis/Leidenfrost point, and regime transitions, including spontaneous nucleation, bubble ebullition, coalescence, and film insulation (Saito et al., 2019, Li et al., 2015, Li et al., 2020, Mukherjee et al., 2020).
• Thermodynamic Consistency: Recovery of the Maxwell coexistence curve, capturing the correct latent heat, and physical surface tension is assessed via comparison with thermodynamic predictions, such as Clausius–Clapeyron relations and binodal lever rules (Biferale et al., 2011, Saito et al., 2019).
• Hydrodynamic and Heat Transfer Benchmarks: Quantitative accuracy in bubble departure diameter vs. gravity (Fritz law), Nusselt vs. Jacob number scaling, mean and fluctuating heat-transfer coefficients, and pressure drop evolution are central benchmarks (Saito et al., 2019, Mukherjee et al., 2020, Iskhakova et al., 2022, Fu et al., 2017, Končar et al., 22 Sep 2025).
• Interfacial Dynamics: Models are probed for their ability to predict nucleation site sensitivity to contact angle/wettability, density depletion rates near solid boundaries, and preferential nucleation at concave corners (Li et al., 2020, Li et al., 2015).
• Statistical and Spectral Diagnostics: Machine-learned and surrogate models are validated by energy spectrum comparisons (e.g., $E_T(k) \sim k^{-5/3}$ scaling in inertial range), mass conservation (relative vapor volume error), and interface regularity (Eikonal loss) (Na et al., 27 Jan 2025, Hassan et al., 28 Jul 2025).

4. Boundary and Initial Conditions, Wettability, and Surface Effects

Application of boundary and initial conditions is crucial:

• Thermal BCs: Fixed (Dirichlet) or flux (Neumann) temperature boundary conditions at heaters, and isothermal or adiabatic walls, are implemented depending on the scenario (Li et al., 2015, Li et al., 2020, Saito et al., 2019).
• Contact Angle Modeling: Wettability enters via wall-pseudopotential (LBM), fluid-solid adhesion forces, or subgrid contact angle models (level-set). Implementation is typically through an adjusted wall density or surface force parameter to enforce Young's law (Li et al., 2020, Li et al., 2015, Iskhakova et al., 2022).
• Structured Surfaces: In LBM, the density depletion at structured/cavitated surfaces is enhanced at concave corners by geometric amplification of the pseudopotential force, explaining observed preferential nucleation phenomena (Li et al., 2020).
• Inflow/Outflow and Numerical Stability: Convective, periodic, or sponge (vortex-damping) outflow conditions are adopted to ensure physically realistic boundary behavior and suppress upstream-propagating disturbances due to outgoing bubbles or phase change (Dhruv, 2023, Saito et al., 2019).

5. Recent Developments, Limitations, and Future Directions

Boiling flow algorithm research addresses both physical fidelity and numerical robustness:

• Advanced Stability and Consistency: The central-moment LBM and thermodynamic-consistency corrections improve algorithmic stability at high Reynolds numbers and ensure correct interface dynamics, crucial for simulating forced-convection boiling (Saito et al., 2019).
• Interface Momentum Correction: Recent work formalizes and corrects the momentum imbalance arising in non-conservative one-fluid VOF formulations with phase change. Corrective interface-localized forces and predictor-projection decompositions are introduced to precisely recover interfacial momentum jumps due to evaporation—a requirement for accurate film-boiling instability simulation at high viscosity (Poblador-Ibanez et al., 2024).
• Two-Phase/Two-Fluid Closure Model Limitations: In RPI heat-flux partitioning models, deficiencies in the closure for bubble departure frequency, bubble size, and wall area fractions at extreme velocities and heat fluxes lead to underprediction of heat transfer and instability. Remedial adjustments, such as boosting bubble frequency and capping area fractions, are practical but highlight a need for improved mechanistic models (Končar et al., 22 Sep 2025).
• Data-Driven and Surrogate Modeling: ML surrogates have demonstrated the ability to match physical heat flux, interface dynamics, and evolution metrics at orders-of-magnitude lower computational cost and without requiring explicit specification of nucleation events or interface-seeding, provided sufficient diversity and quality of simulation data (Hassan et al., 28 Jul 2025, Na et al., 27 Jan 2025).
• Fidelity Trade-Offs: Simplified, computationally efficient boiling flow algorithms (e.g., for aero-optic turbulence) trade strict physical accuracy for statistical or spectral concordance, with tunable parameters to match observed temporal or spatial statistics (structure functions, power spectra) (Utley et al., 17 Jan 2026).

6. Summary Table: Key Algorithmic Paradigms in Boiling Flow Simulation

Algorithm Class	Key Features	References
Pseudopotential LBM (BGK/MRT/CM)	Mesoscopic, spontaneous phase change, thermal LBE	(Saito et al., 2019, Li et al., 2020, Li et al., 2015, Mukherjee et al., 2020)
Finite-Element Level-Set	High-res; explicit interface, bubble tracking, surface tension, contact angle	(Iskhakova et al., 2022)
One/Two-Fluid Navier-Stokes	Phase fraction, wall heat partitioning, CFD closures	(Fu et al., 2017, Končar et al., 22 Sep 2025, Poblador-Ibanez et al., 2024)
VOF with Corrected Momentum	Exact interface momentum/pressure jump, phase-wise stencils	(Poblador-Ibanez et al., 2024)
ML Surrogate (Transformer, Diffusion)	Forecast interface, heat transfer, flow states	(Hassan et al., 28 Jul 2025, Na et al., 27 Jan 2025)
Spectral "Boiling Flow" (aero-optics)	AR(1) in Fourier domain; match structure function/PSD	(Utley et al., 17 Jan 2026)

Each algorithm exhibits trade-offs between physical interpretability, computational cost, scalability, and generalizability. State-of-the-art methods now routinely achieve physically consistent boiling curves, interface instabilities, nucleation, and mass/energy conservation across both DNS-style and surrogate-model domains. Remaining challenges include subgrid modeling of turbulence and interface phenomena, robust energetic closure at high heat flux and complex geometries, and reliable extrapolation to diverse fluids and operating conditions.