A robust incompressible Navier-Stokes solver for high density ratio multiphase flows

TL;DR

This paper introduces a stable, adaptive Navier-Stokes solver for high density ratio multiphase flows that maintains accuracy and stability through consistent mass-momentum treatment and advanced numerical techniques.

Contribution

The paper develops a robust, second-order accurate incompressible Navier-Stokes solver with adaptive refinement and consistent mass-momentum transport for high density ratio multiphase flows.

Findings

Stable resolution of material contrast ratios up to six orders of magnitude

Demonstrated second-order accuracy and scalability of the flow solver

Highlighting the importance of consistent mass and momentum transport for accuracy

Abstract

This paper presents a robust, adaptive numerical scheme for simulating high density ratio and high shear multiphase flows on locally refined Cartesian grids that adapt to the evolving interfaces and track regions of high vorticity. The algorithm combines the interface capturing level set method with a variable-coefficient incompressible Navier-Stokes solver that is demonstrated to stably resolve material contrast ratios of up to six orders of magnitude. The discretization approach ensures second-order pointwise accuracy for both velocity and pressure with several physical boundary treatments, including velocity and traction boundary conditions. The paper includes several test cases that demonstrate the order of accuracy and algorithmic scalability of the flow solver. To ensure the stability of the numerical scheme in the presence of high density and viscosity ratios, we employ a consistent treatment of mass and momentum transport in the conservative form of discrete equations. This consistency is achieved by solving an additional mass balance equation, which we approximate via a strong stability preserving Runga-Kutta time integrator and by employing the same mass flux (obtained from the mass equation) in the discrete momentum equation. The scheme uses higher-order total variation diminishing (TVD) and convection-boundedness criterion (CBC) satisfying limiter to avoid numerical fluctuations in the transported density field. The high-order bounded convective transport is done on a dimension-by-dimension basis, which makes the scheme simple to implement. We also demonstrate through several test cases that the lack of consistent mass and momentum transport in non-conservative formulations, which are commonly used in practice, or the use of non-CBC satisfying limiters can yield very large numerical error and very poor accuracy for convection-dominant high density ratio flows.