A novel momentum-conserving, mass-momentum consistent method for interfacial flows involving large density contrasts

TL;DR

This paper introduces a new volume-of-fluid numerical method for simulating interfacial flows with large density differences, ensuring mass-momentum consistency and high accuracy in complex turbulent scenarios.

Contribution

It presents a novel, conservative, mass-momentum consistent VOF method with geometric mass advection and enhanced stability for large density contrast flows.

Findings

Accurately simulates interfacial flows with large density contrasts.

Demonstrates stability and robustness in turbulent flow simulations.

Achieves high precision in mass and momentum transport.

Abstract

We propose a novel method for the direct numerical simulation of interfacial flows involving large density contrasts, using a Volume-of-Fluid method. We employ the conservative formulation of the incompressible Navier-Stokes equations for immiscible fluids in order to ensure consistency between the discrete transport of mass and momentum in both fluids. This strategy is implemented on a uniform 3D Cartesian grid with a staggered configuration of primitive variables, wherein a geometrical reconstruction based mass advection is carried out on a grid twice as fine as that for the momentum. The implementation is in the spirit of Rudman (1998) [41], coupled with the extension of the direction-split time integration scheme of Weymouth & Yue (2010) [46] to that of conservative momentum transport. The resulting numerical method ensures discrete consistency between the mass and momentum propagation, while simultaneously enforcing conservative numerical transport to arbitrary levels of precision in 3D. We present several quantitative comparisons with benchmarks from existing literature in order to establish the accuracy of the method, and henceforth demonstrate its stability and robustness in the context of a complex turbulent interfacial flow configuration involving a falling raindrop in air.