Prediction of bubble departure in forced convection boiling: A mechanistic model

In the context of computational fluid dynamic simulations of boiling flows using time-averaged Eulerian multi-phase approaches, the many sub-models required to describe such a complex phenomena are of particular importance. Of interest here, wall boiling requires calculation of the contribution of evaporation to global heat transfer, which in turn relies on determination of the active nucleation site density, bubble departure diameter and frequency of bubble departure. In this paper, an improved mechanistic model for the bubble departure diameter during flow boiling is developed. The model is based on the balance of forces acting on a bubble at a single nucleation site, with a new equation governing bubble growth proposed. The formulation accounts for evaporation of the micro-layer under the bubble, heat transfer from superheated liquid around the bubble surface, and condensation on the bubble cap due to the presence of sub-cooled liquid. Validation of the growth equation is provided through comparison against experiments in both pool boiling and flow boiling conditions. Introduction of condensation on the bubble cap allows reproduction of the growth of the bubble for different sub-cooling temperatures of the surrounding liquid. In addition, a sensitivity study guarantees dependency of the bubble departure diameter on relevant physical quantities such as mass flow rate, heat flux, liquid sub-cooling and pressure, with any inclination of the channel walls correctly accounted for. Predictions of bubble departure diameter and bubble lift-off are validated against three different databases on sub-cooled flow boiling with water and an additional database on saturated boiling with refrigerant R113. The whole data set guarantees validation is performed over a range of parameters and operating conditions as broad as possible. Satisfactory predictive accuracy is obtained in all conditions. The present formulation provides an appropriate starting point for prediction of the behaviour of vapour bubbles under more general conditions which include lift-off after sliding, the frequency of bubble departure, bubble merging and bubble shrinking and collapse due to condensation.