A geometric VOF method for interface resolved phase change and conservative thermal energy advection

JOURNAL OF COMPUTATIONAL PHYSICS

卷:426

Article

Malan, L. C.^1,2,3  Malan, A. G.³  Zaleski, S.^1,2  Rousseau, P. G.⁴ 

[1] Sorbonne Univ, Paris, France

[2] CNRS, Inst Jean Rond Alembert, UMR 7190, Paris, France

[3] Univ Cape Town, Ind CFD Res Grp, Dept Mech Engn, Private Bag X3, ZA-7701 Rondebosch, South Africa

[4] Univ Cape Town, Dept Mech Engn, CATPRoM Res Grp, Private Bag X3, ZA-7701 Rondebosch, South Africa

关键词: VOF;    Incompressible;    Phase change;    Thermal energy conservation;    Finite volume;   

DOI  :  10.1016/j.jcp.2020.109920

来源: Elsevier

【 摘 要 】

We present a novel numerical method to solve the incompressible Navier-Stokes equations for two-phase flows with phase change. Separate phases are tracked using a geometric Volume-Of-Fluid (VOF) method with piecewise linear interface construction (PLIC). Thermal energy advection is treated in conservative form and the geometric calculation of VOF fluxes at computational cell boundaries is used consistently to calculate the fluxes of heat capacity. The phase boundary is treated as sharp (infinitely thin), which leads to a discontinuity in the velocity field across the interface in the presence of phase change. The numerical difficulty of this jump is accommodated with the introduction of a novel two-step VOF advection scheme. The method has been implemented in the open source code PARIS and is validated using well-known test cases. These include an evaporating circular droplet in microgravity (2D), the Stefan problem and a 3D bubble in superheated liquid. The accuracy shown in the results were encouraging. The 2D evaporating droplet showed excellent prediction of the droplet volume evolution as well as preservation of its circular shape. A relative error of less than 1% was achieved for the Stefan problem case, using water properties at atmospheric conditions. Two cases of a bubble in superheated liquid was performed, at respective Jacob numbers of 0.5 and 2.15. For the final radius of the bubbles in both cases, a relative error of less than 7% was obtained on the coarsest grid, with less than 1% on the finest. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free