Application Features of a VOF Method for Simulating Boiling and Condensation Processes

This article presents the results of a study on the possibility of using a single-speed multiphase model with free surface allowance for simulating boiling and condensation processes. The simulation is based on the VOF method, which allows the position of the interphase boundary to be tracked. To increase the stability of the iterative procedure for numerically solving volume fraction transfer equations using a finite volume discretization method on arbitrary unstructured grids, the basic VOF method is been modified by writing these equations in a semi-divergent form. The models of Tanasawa, Lee, and Rohsenow are considered models of interphase mass transfer, in which the evaporated or condensed mass linearly depends on the difference between the local temperature and the saturation temperature with accuracy in empirical parameters. This paper calibrates these empirical parameters for each mass transfer model. The results of our study of the influence of the values of the empirical parameters of models on the intensity of boiling and evaporation, as well as on the dynamics of the interphase boundary, are presented. This research is based on Stefan’s problem of the movement of the interphase boundary due to the evaporation of a liquid and the problem of condensation of vapor bubbles water columns. As a result of a series of numerical experiments, it is shown that the average error in the position of the interfacial boundary for the Tanasawa and Lee models does not exceed 3–6%. For the Rohsenow model, the result is somewhat worse, since the interfacial boundary moves faster than it should move according to calculations based on analytical formulas. To investigate the possibility of condensation modeling, the results of a numerical solution of the problem of an emerging condensing vapor bubble are considered. A numerical assessment of its position in space and the shape and dynamics of changes in its diameter over time is carried out using the VOF method, taking into account the free surface. It is shown herein that the Tanasawa model has the highest accuracy for modeling the condensation process using a VOF method taking into account the free surface, while the Rohsenow model is most unstable and prone to deformation of the bubble shape. At the same time, the dynamics of bubble ascent are modeled by all three models. The results obtained confirm the fundamental possibility of using a VOF method to simulate the processes of boiling and condensation and taking into account the dynamics of the free surface. At the same time, the problem of the studied models of phase transitions is revealed, which consists of the need for individual selection of optimal values of empirical parameters for each specific task.