Thermal Bubble Algorithms

• Thermal bubble algorithms are computational frameworks that model bubble nucleation, growth, and interactions driven by temperature gradients and latent heat.
• These methods employ sharp-interface, diffuse-interface, and lattice kinetic schemes to couple fluid dynamics, energy conservation, and interfacial forces.
• Applications span boiling heat transfer, cosmological phase transitions, microgravity bubble dynamics, and turbulent nucleation, validated by rigorous benchmarks.

Thermal bubble algorithms are a class of computational methods and theoretical frameworks for modeling, simulating, and predicting the nucleation, growth, migration, and interaction of bubbles in systems driven by thermal effects. These algorithms underlie critical research areas in multiphase fluid dynamics, phase transitions in condensed matter and cosmology, boiling heat transfer, and interface-driven flows. They encompass sharp- and diffuse-interface formalisms, volume- and boundary-integral methods, fluctuating hydrodynamics, effective field theory, and lattice kinetic schemes, attaining high-fidelity treatment of interface evolution, latent heat, Marangoni stresses, and thermally activated nucleation.

1. Fundamental Principles and Classes of Thermal Bubble Algorithms

Thermal bubble algorithms are unified by their treatment of bubbles as interfacial structures in a medium where temperature gradients, phase change, and thermodynamic fluctuations play central roles. The dominant physical phenomena include:

• Thermodynamic instability, nucleation, and growth: Bubbles nucleate when fluctuations surmount a free-energy barrier, with subsequent expansion governed by the balance of capillarity, vapor recoil, and latent heat transport (Gallo et al., 2018, Gould et al., 2021).
• Interfacial transport and forces: Surface-tension gradients (Marangoni stresses) drive migration; buoyancy contributes under gravity (Tripathi et al., 2018).
• Latent-heat coupling: Energy conservation at the interface includes latent heat absorption or release, which drives the Stefan condition and dynamically reconfigures temperature and flow fields (Malan et al., 2020).
• Stochastic and deterministic evolution: Both deterministic partial differential equations (PDEs) and stochastic fluctuating hydrodynamics are used, depending on whether thermal noise and rare nucleation events are critical (Gallo et al., 2018).

Algorithmic frameworks are delineated as follows:

Methodology	Physical Regime	Interface Treatment
Volume-of-Fluid (VoF)	Two-phase thermal convection, phase change	Sharp, reconstructed
Boundary Element Method (BEM)	Thermal diffusion–dominated boiling	Sharp, moving contour
Diffuse Interface/Fluctuating Hydro	Thermally activated nucleation	Diffuse, stochastic
Effective Field Theory	Quantum/classical phase transitions	Phase field, bounce soln
Lattice Boltzmann	Boiling, multiphase convection	Emergent, kinetic grid
Multigrid All-Mach Solvers	Compressible phase-change, sonoluminescence	VoF, fully coupled

This diversity reflects the range of application scales, from molecular and statistical physics to continuum thermal-fluid mechanics and cosmological or quantum field environments.

2. Mathematical Formulations and Governing Equations

Thermal bubble algorithms are typically based on the coupled evolution of mass, momentum, energy (temperature), and interface geometry. Key mathematical structures include:

Navier–Stokes and Transport Equations with Interfacial Forces

$$\rho(\partial_t\mathbf{u} + \mathbf{u}\cdot\nabla\mathbf{u}) = -\nabla p + \nabla\cdot[\mu(\nabla\mathbf{u} + \nabla\mathbf{u}^T)] + \rho\mathbf{g} + \mathbf{F}_\sigma$$

$$\rho c_p \left( \frac{\partial T}{\partial t} + \mathbf{u}\cdot\nabla T \right) = \nabla\cdot(k\nabla T) + S_{\text{pc}}$$

where $\mathbf{F}_\sigma$ encompasses curvature and Marangoni forces (Tripathi et al., 2018), and $S_{\text{pc}}$ is the latent-heat source/sink (Malan et al., 2020).

Sharp-interface bubble evolution (Stefan condition)

$$V_n = \frac{[k_\ell (\partial_n T)_\ell - k_g (\partial_n T)_g]}{\rho_\ell h_{fg}}$$

(Malan et al., 2020)

Thermal nucleation rates—free-energy barriers and stochastic PDEs

The nucleation rate $J$ for thermal bubbles is governed by the Arrhenius–like escape over a free-energy barrier $\Delta\Omega$, often modeled by Langevin/field-theoretical or classical density functional equations (Gallo et al., 2018, Gould et al., 2021):

$J \propto \exp(-\Delta\Omega / k_BT)$

Phase-field/bounce formalisms for first-order transitions:

The nucleation rate per unit volume is given by

$$\Gamma = \frac{\kappa}{2\pi} \Sigma, \quad \Sigma = V\left( \frac{S_3}{2\pi T} \right)^{3/2} \left| \frac{\det[-\nabla^2 + V^{\prime\prime}_{3,\text{meta}}]}{\det'[-\nabla^2 + V^{\prime\prime}_{3,\text{bounce}}]} \right|^{1/2} e^{-S_3/T}$$

where $S_3$ is the three-dimensional Euclidean bounce action, and the determinants account for quantum/thermal fluctuations around the critical bubble (Gould et al., 2021).

3. Numerical Methodologies and Discretization Techniques

Algorithmic innovations specifically tailored for thermal bubble problems include:

Volume-of-Fluid schemes with conserved thermal energy and sharp interface:

• Geometric PLIC interface reconstruction.
• Two-step VOF advection to handle discontinuities in interface velocities due to phase change—liquid extension via local Poisson solve, followed by explicit interface shift (Malan et al., 2020).
• Conservative advection of $\rho c_p$ ensures consistent energy fluxes; consistent with changes in phase fraction.

Height-function/column-averaged surface tension gradients:

• Marangoni forces sharply localized by computing column-averaged $\sigma_c$ and high-order finite differences in the direction tangential to the interface, avoiding spurious smoothing (Tripathi et al., 2018).

Boundary Element Method for heat-diffusion-limited growth:

• Time-dependent single-layer potential formulation for the heat equation in multiply connected, moving boundaries (Nikolayev et al., 2016).
• Analytical and quadrature-based treatment of singular kernels at the triple point (liquid–solid–vapor contact).
• Iterative mesh adaptation and contour reconstruction based on vapor recoil and capillarity.

Fluctuating hydrodynamics for stochastic nucleation:

• Second-order stochastic Runge–Kutta integration, spatial discretization with central finite-volume/difference (Gallo et al., 2018).
• Thermal noise terms constructed to satisfy fluctuation-dissipation balance.
• Automated cluster analysis for bubble identification and nucleation-rate extraction.

Multigrid solvers for implicit coupling of thermal, pressure, and interface dynamics:

• Fully coupled Poisson-Helmholtz solves for pressure and temperature using V-cycle multigrid with $O(N)$ complexity (Saade et al., 2022).

Lattice Boltzmann kinetic methods with phase change:

• D3Q19 or similar stencils for advection and collision; non-ideal equation of state for phase separation and Clausius–Clapeyron relation (Biferale et al., 2011).
• Latent heat captured by source terms in the energy population, with correct thermodynamic response to boiling and condensation.

4. Validation, Precision, and Benchmarking

Thermal bubble algorithms have undergone rigorous validation:

• Terminal velocity and migration: Accurate recovery of theoretical terminal velocity for thermocapillary migration under imposed temperature gradients, with errors below 1.6% and better performance than previous VoF/front-tracking results (Tripathi et al., 2018).
• Bubble growth in superheated liquid: For 3D PLIC-VOF, errors in final radius scale down to $<1\%$ on fine grids, and convergence is monotonic (Malan et al., 2020).
• Epstein–Plesset and Rayleigh collapse: Multigrid VoF solver quantitatively matches analytical solutions for both isothermal and adiabatic transition regimes (Saade et al., 2022).
• Boiling convection: Lattice Boltzmann framework quantitatively recovers Clausius–Clapeyron and exhibits $20\%$ Nusselt number enhancement and strongly non-Gaussian statistics in the presence of bubbles (Biferale et al., 2011).
• Stochastic nucleation: Fluctuating hydrodynamics approach yields nucleation rates $J$ and barrier heights $\Delta\Omega$ in agreement with classical nucleation theory and molecular dynamics, within an order of magnitude (Gallo et al., 2018).
• Cosmological bubble wall velocities: Local thermal equilibrium (LTE) approaches with uniform-sound-speed ansatz provide model-independent estimates of wall speed with $<$10\% error for moderate parameter regimes, validated using fit formulas and code snippets (Ai et al., 2023).

5. Applications and Physical Phenomena Captured

Thermal bubble algorithms have enabled predictive study in domains including:

• Boiling and heat transfer: Modeling heat-flux-driven bubble growth, dry spot formation, and boiling crisis (CHF) in pool boiling scenarios (Nikolayev et al., 2016).
• Thermocapillary migration: Simulations of Marangoni-driven bubble/droplet motion in both microgravity and gravitational environments, including high density and viscosity ratio regimes (Tripathi et al., 2018).
• Phase transition dynamics in the early universe: Calculation of bubble wall velocities, nucleation kinetics, and absorption of latent heat during first-order phase transitions (e.g., baryogenesis and gravitational wave generation) (Ai et al., 2023, Gould et al., 2021).
• Sonoluminescence and collapse: Accurate reproduction of SBSL and Rayleigh collapse waveforms, quantification of heat flux and damping dependent on collapse intensity and proximity to walls (Saade et al., 2022).
• Turbulent boiling: Analysis of turbulence statistics, intermittency, and enhancement of heat transport due to nucleate boiling (Biferale et al., 2011).
• Thermally activated nucleation in statistical fluids: Calculation of rates, lifetimes, and size spectra in metastable, overheated liquids (Gallo et al., 2018).

6. Algorithmic Innovations and Performance Considerations

Salient algorithmic advances include:

• Non-smeared Marangoni forcing: Elimination of artificial smoothing in computation of $\nabla_s\sigma$ yields formally second-order sharp localization and high accuracy even at strong surface tension gradients (Tripathi et al., 2018).
• Discontinuous-velocity two-step VOF advection: Extension-shift scheme handles divergence in velocity field near phase boundaries during evaporation or condensation (Malan et al., 2020).
• Adaptive mesh refinement and graded meshing: Concentrates resolution near moving interfaces and critical points (e.g., contact lines), essential for resolving thin boundary layers (Tripathi et al., 2018, Nikolayev et al., 2016).
• Coupled multigrid for strongly stiff systems: Enables robust and efficient solution of the coupled pressure–temperature field, bypassing explicit diffusion timestep restrictions (Saade et al., 2022).
• High-parallel scaling: For instance, in geometric VOF-PLIC routines, near-ideal strong scaling to $\sim4000$ cores for large 3D problems (Malan et al., 2020).
• Open-source implementations: Many methods are implemented in codes such as Basilisk, PARIS, and actively compared against alternative solvers and analytical benchmarks.

7. Limitations, Extensions, and Outlook

• Physical modeling limits: Most current sharp-interface methods assume measurable surface tension and neglect interface diffuseness at molecular scales; diffusive and stochastic approaches can bridge to mesoscale and nanoscopic regimes (Gallo et al., 2018).
• Treatment of turbulence and fluctuations: While fluctuating hydrodynamics and LB approaches naturally accommodate noise and intermittent events, most deterministic CFD frameworks neglect rare nucleation unless explicitly coupled with stochastic modeling (Gallo et al., 2018, Biferale et al., 2011).
• Quantum and cosmological regimes: Effective field theory and field-theoretic bounce algorithms provide first-principles predictions in contexts where traditional hydrodynamics fails, e.g., electroweak baryogenesis or gravitational wave backgrounds (Gould et al., 2021, Ai et al., 2023).
• Algorithmic improvements: Recommendations include higher-order discretizations, more robust coupling of interface geometry and thermodynamic fields, and systematic inclusion of fluctuation-driven nucleation and interface instabilities (Gallo et al., 2018, Gould et al., 2021, Saade et al., 2022).
• Model independence and universality: The transition to fit-based, input-invariant approaches (e.g., for wall velocities in phase transitions) promotes rapid assessment and cross-model comparison (Ai et al., 2023).

Thermal bubble algorithms are essential computational and theoretical tools for quantitatively accurate, first-principles modeling of phase-change phenomena driven by thermal mechanisms, and their continued development, validation, and extension underpin advances across fluid dynamics, statistical mechanics, material science, and cosmology.