Multiphase Flows of N Immiscible Incompressible Fluids: A Reduction-Consistent and Thermodynamically-Consistent Formulation and Associated Algorithm

TL;DR

This paper introduces a new formulation and algorithm for simulating multiple immiscible incompressible fluids that is both reduction-consistent and thermodynamically consistent, ensuring accurate modeling across different fluid combinations.

Contribution

The paper develops a systematic method for constructing reduction-consistent, thermodynamically consistent N-phase flow models, applicable to fluids with varying properties.

Findings

Accurately simulates flows with multiple fluids and large property ratios.

Reduces correctly to fewer-phase models when some fluids are absent.

Produces results consistent with physical theories and solutions.

Abstract

We present a reduction-consistent and thermodynamically consistent formulation and an associated numerical algorithm for simulating the dynamics of an isothermal mixture consisting of $N$ ($N\geqslant 2$) immiscible incompressible fluids with different physical properties (densities, viscosities, and pair-wise surface tensions). By reduction consistency we refer to the property that if only a set of $M$ ($1\leqslant M\leqslant N-1$) fluids are present in the system then the N-phase governing equations and boundary conditions will exactly reduce to those for the corresponding $M$-phase system. By theromdynamic consistency we refer to the property that the formulation honors the thermodynamic principles. Our N-phase formulation is developed based on a more general method that allows for the systematic construction of reduction-consistent formulations, and the method suggests the existence of many possible forms of reduction-consistent and thermodynamically consistent N-phase formulations. Extensive numerical experiments have been presented for flow problems involving multiple fluid components and large density ratios and large viscosity ratios, and the simulation results are compared with the physical theories or the available physical solutions. The comparisons demonstrate that our method produces physically accurate results for this class of problems.