A volume of fluid method for three dimensional direct numerical simulations of immiscible droplet collisions

This paper presents an advanced Volume of Fluid (VOF) method that enables performant three dimensional Direct Numerical Simulations (DNS) of the interaction of two immiscible fluids in a gaseous environment with large topology changes, e.g., binary droplet collisions. One of the challenges associated with the introduction of a third immiscible phase into the VOF method is the reconstruction of the phase boundaries near the triple line in arbitrary arrangements. For this purpose, an efficient method based on a Piecewise Linear Interface Calculation (PLIC) is shown. Moreover, the surface force modeling with the robust Continuous Surface Stress (CSS) model was enhanced to treat such three-phase situations with large topology changes and thin films. A consistent scaling of the fluid properties at the interfaces ensures energy conservation.The implementation of these methods in the multi-phase flow solver Free Surface 3D (FS3D) allowed a successful validation. A qualitative comparison of the morphology in binary collisions of immiscible droplets as well as a quantitative comparison regarding the threshold velocities that distinguish different collision regimes shows excellent agreement with experimental results.These simulations enable the evaluation of experimentally inaccessible data like the contributions of kinetic, surface and dissipative energy of both immiscible liquids during the collision process. Furthermore, the comparison against binary collisions of the same liquids highlights similarities and differences between immiscible and equal droplet collisions. Both can support the modeling of the immiscible liquid interaction in the future.

Automated summary

This paper presents an advanced Volume of Fluid (VOF) method for performing three-dimensional Direct Numerical Simulations (DNS) of the interaction of two immiscible fluids in a gaseous environment with large topology changes. The method includes efficient reconstruction of phase boundaries near the triple line using a Piecewise Linear Interface Calculation (PLIC) method and enhanced surface force modeling with the Continuous Surface Stress (CSS) model. Implementation of these methods in the multi-phase flow solver Free Surface 3D (FS3D) yielded successful validation. The simulations provide valuable insights into the collision process and can support future modeling of immiscible liquid interaction.