  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., 400: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.
