Numerical simulation of boiling on unstructured grids

• Published in: Journal of computational physics 2021-05, Vol.432, p.110161, Article 110161
• Main Authors: Sahut, Guillaume, Ghigliotti, Giovanni, Balarac, Guillaume, Bernard, Manuel, Moureau, Vincent, Marty, Philippe
• Format: Article
• Language: English
• Subjects: Boiling; Boundary conditions; Bubble growth; Cartesian coordinates; Computational fluid dynamics; Computational physics; Differential equations; Dirichlet problem; Engineering Sciences; Fluid mechanics; Fluids mechanics; Grid refinement (mathematics); Heat flux; Heat transmission; Latent heat; Mass transfer; Mass transfer rate; Mechanics; Multiphase flow; Numerical methods; Operators (mathematics); Phase change; Physics; Simulation; Two-phase flows; Unstructured grids; Unstructured grids (mathematics)
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2021.110161

A numerical method for the simulation of multiphase flows with phase change on unstructured grids is presented. Based on the work of Tanguy et al. [31] using a Level Set/Ghost Fluid Method coupling for two-dimensional axisymmetric cartesian grids, we extend the method to two- and three-dimensional unstructured grids with the aim of taking a step towards realistic boiling simulations in industrial context. The mass transfer rate at the interface accounting for phase change is computed as a function of the liquid and vapor heat fluxes at the interface and of the latent heat of the fluid by means of a new framework improving the accuracy of differential operators. The mass transfer rate is then used in the projection method to solve Navier-Stokes equations and in the advection equation of the Level Set function to account for the interface movement due to phase change. An immersed Dirichlet boundary condition is imposed at the interface to ensure that boiling always occurs at saturation temperature. We demonstrate the accuracy of our method first in the case of a static growing bubble with a fixed constant mass transfer rate, and second with a mass transfer rate computed from the thermal fluxes at the interface. In both cases, the bubble radius at final time converges with grid refinement towards the theoretical value. •Boiling is simulated as a coupling between the Navier-Stokes and heat equations.•Boiling is driven by the difference of the heat fluxes at the two-phase interface.•The mass transfer rate is computed precisely at the subgrid interface.•The interface is captured by a Level Set function taking phase change into account.•The velocity field is discontinuous at the interface due to phase change.