In this work, we present a dynamical theory of boiling based on fluctuating hydrodynamics and the diffuse interface approach. The model is able to describe boiling from the stochastic nucleation up to the macroscopic bubble dynamics. It covers, with a modest computational cost, the mesoscale area from nano to micrometers, where most of the controversial observations related to the phenomenon originate. In particular, the role of wettability in the macroscopic observables of boiling is elucidated. In addition, by comparing the ideal case of boiling on ultra-smooth surfaces with a chemically heterogeneous wall, our results will definitively shed light on the puzzling low onset temperatures measured in experiments. Sporadic nanometric spots of hydrophobic wettability will be shown to be enough to trigger the nucleation at low superheat, significantly reducing the temperature of boiling onset, in line with experimental results. The proposed mesoscale approach constitutes the missing link between macroscopic approaches and molecular dynamics simulations and will open a breakthrough pathway toward accurate understanding and prediction.

The present investigation, identifies the exact quantitative effects of fundamental parameters, on the detachment characteristics of isolated bubbles, emanating quasi-statically from submerged orifices into isothermal liquid pools. For this purpose, a Volume of Fluid (VOF) based interface capturing approach is further improved, for the conduction of axisymmetric and 3D numerical experiments on adiabatic bubble growth dynamics. The predictions of the model, are quantitatively validated against literature available experimental data, showing excellent agreement. Two series of numerical experiments are performed, quantitatively exploring the parametric effects of the liquid phase properties in five different gravity levels, and the effect of the gravity vector direction inclination angle, respectively. It is found that the bubble detachment characteristics, are more sensitive in the variation of the surface tension, liquid phase density and gravity, while the effect of liquid phase dynamic viscosity is generally minimal. From dimensionless analysis, two correlations are derived, which for the examined range of Eötvos numbers, are able to predict the equivalent bubble detachment diameter and the bubble detachment time, respectively. It is also found that the bubble detachment characteristics, reduce significantly as the gravity vector direction gradually deviates from being parallel to the bubble injection orifice, following a non-linear decrease