Numerical simulation of saturated film boiling on a horizontal surface

The past efforts in applying linear Taylor instability theory to the prediction of heat transfer during film boiling on a horizontal surface have suffered from the fact that empirical correlations must be used to define the shape of vapor-liquid interfaces and to determine the transport of mass and heat across these interfaces. The objective of this study is to clarify the physics of film boiling and to predict heat transfer coefficients through complete numerical simulation of the evolving interface between superposed layers of immiscible fluids. A coordinate transformation technique supplemented by a numerical grid generation method and a second-order projection method are combined to solve for the flow and temperature fields associated with an evolving interface. From the numerical simulation, the film thickness and, in turn, the heat transfer coefficient are found to vary both spatially and temporally. Increased wall superheat not only thickens the vapor film in the valley but also enlarges the vapor bulge. The effect of increased system pressure is to slow down the growth of the interface.