Stabilizing the unstructured Volume-of-Fluid method for capillary flows in microstructures using artificial viscosity

Parasitic currents still pose a significant challenge for the investigation of two-phase flow in Lab-on-Chip (LoC) applications with Volume-of-Fluid (VoF) simulations. To counter the impact of such spurious velocity fields in the vicinity of the fluid interface, this work presents an implementation of an artificial interface viscosity model in OpenFOAM. The model is introduced as an additional damping term in the momentum conservation equation. It is implemented as an fvOption, allowing for its simple application to existing VoF solvers. Validation is performed with hydrodynamic and wetting cases, in which constant artificial viscosity values are prescribed to examine the sensitivity of the solution to the artificial damping. The artificial viscosity model shows promising results in reducing spurious currents for two considered geometrical VoF solvers, namely interIsoFoam and InterFlow. It is found that the influence of the artificial viscosity heavily depends on the fluid properties. Applying the model to simulations of an interface traversing through microcavities relevant in LoC applications, experimental results of the interface progression are predicted well, while spurious currents are effectively reduced by approximately one order of magnitude due to the artificial viscosity model. The code is publicly available on GitHub at https://github.com/boschresearch/sepMultiphaseFoam/tree/publications/ArtificialInterfaceViscosity.