Abstract

The numerical simulation of multiphase flows presents several challenges, namely the transport of different phases within de domain and the inclusion of capillary effects. Here, these are approached by enforcing a discrete physics-compatible solution. Extending our previous work on the discretization of surface tension [N. Valle, F. X. Trias, and J. Castro, 'An energy-preserving level set method for multiphase flows,' J. Comput. Phys., vol. 400, p. 108991, 2020] with a consistent mass and momentum transfer a fully energy-preserving multiphase flow method is presented. This numerical technique is showcased within the simulation of a falling film under several working conditions related to the normal operation of LiBr absorption chillers.