A novel approach for the solution of Stefan problem within the framework of the multi fluid model supplemented with Volume of Fluid (VOF) method, i.e. two-fluid VOF, is presented in this paper. The governing equation set is comprised of mass, momentum and energy conservation equations, written on a per phase basis and supplemented with closure models via the source terms. In our method, the heat and mass transfer is calculated from the heat transfer coefficient, which has a fictitious function and depends on the local cell size and the thermal conductivity, and the implementation is straightforward because of the usage of the local value instead of a global parameter. The interface sharpness is ensured by the application of the geometrical reconstruction scheme implemented in VOF. The model is verified for three types of computational meshes including triangular cells, and good agreement was obtained for the interface position and the temperature field. Although the developed method was validated only for Stefan problem, the application of the method to engineering problems is considered to be straightforward since it is implemented to a commercial CFD code only using a local value; especially in the field of naval hydrodynamics wherein the reduction of ship resistance using boiling flow can be computed efficiently since the method handles phase change processes using low resolution meshes.

The DNS of film boiling requires strong computational resources that are difficult to obtain for daily CFD use by expert practitioners of industrial R&D. On the other hand, film boiling experiments are associated with the usage of expensive and highly sophisticated apparatus, and research to this end is relatively difficult due to high heat flow rates that are present in the process itself. When combined with transient heat conduction in a solid, the problem becomes significantly difficult. Therefore, a novel method in computation of conjugate heat transfer during film boiling in a quiescent liquid is proposed in this paper. The method relies on the solution of mass, momentum and energy conservation equations in a two-fluid framework, supplemented with the appropriate closures. Furthermore, turbulent flow was determined as an important parameter in obtaining an accurate solution to temperature field evolution in a solid specimen, via the proper modeling of the turbulent kinetic energy (TKE) value, that was imposed as a constant value, i.e., the frozen turbulence approach. It was found, in addition, that the appropriate TKE value can be obtained by use of Kelvin–Helmholtz instability theory in conjunction with boundary layer theory. The obtained results show excellent agreement with the experimental data within the first 15 s of the experiment, i.e., the first ca. 10% of the total duration of the film boiling mode of heat transfer. Furthermore, the heat transfer coefficient matched the error bands prescribed by the authors of this paper, which presented the correlations, whilst the averaged values are far beyond this band, i.e., are slightly more than 30% higher. Further inspection revealed a measure of similarity between the computational result of the volume fraction field distribution and the experiment, thus confirming the capability of the method to obtain realistic interface evolution in time. The method shows full capability for further pursuing industrial-scale film boiling problems that involve turbulent flow and the conjugate heat transfer approach.

The DNS of film boiling requires strong computational resources that are difficult to obtain for a daily CFD use by the practitioners of the industrial R&D experts. On the other hand, the film boiling experiments are associated with the usage of the expensive and highly sophisticated apparatus, and the research to this end is found to be relatively difficult due to high heat flow rates that are present in the process itself. When combined with a transient heat conduction in a solid, the problem becomes significantly difficult. Therefore, a novel method in computation of conjugate heat transfer during film boiling in a quiescent liquid has been proposed in this paper. The method relies on the solution of mass, momentum and energy conservation equations in a two-fluid framework, supplemented with the appropriate closures. Furthermore, the turbulent flow has been found as an important parameter in obtaining the accurate solution of the temperature field evolution in a solid specimen, via the proper modeling of turbulent kinetic energy (TKE) value, that has been imposed as a constant value, i.e., the frozen turbulence approach. It was found, in addition, that the appropriate TKE value can be obtained by use of Kelvin-Helmholtz instability theory in conjunction with the boundary layer theory. The obtained results show excellent agreement with the experimental data within the first 15 s of the experiment, i.e., the first ca. 10 % of the total duration of the film boiling mode of heat transfer. Furthermore, the heat transfer coefficient has matched the error bands prescribed by the authors of the paper that has presented the correlations, whilst the averaged values are far beyond this band, i.e., are slightly more than 30 % higher. The further inspection has revealed a measure of similarity between the computational result of the volume fraction field distribution and the experiment, thus confirming the capability of the method to obtain realistic interface evolution in time. The method has shown full capability for further pursuing the industrial-scale film boiling problems that involve turbulent flow and the conjugate heat transfer approach.