Physical Formulation and Numerical Algorithm for Simulating N Immiscible Incompressible Fluids Involving General Order Parameters

Abstract

We present a physical formulation, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of $N$ ($N⩾2$) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The introduction of general order parameters leads to a more strongly coupled system of phase field equations, in contrast to that with certain special choice of the order parameters. However, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. From the simulation perspective, the increased complexity in the form of the phase field equations with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the ($N-1$) strongly-coupled phase field equations for general order parameters into $2(N-1)$ Helmholtz-type equations that are completely de-coupled from one another, leading to a computational complexity essentially the same as that of the simpler phase field equations associated with special choice of order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir-de Gennes theory of floating liquid lenses we show that the method produces physically accurate results for multiple fluid phases.

38 pages, 9 figures, 6 tables

References in corpus (1)

• Surface tension of multi-phase flow with multiple junctions governed by the variational principle

Cited by in corpus (13)

• Multiphase Flows of N Immiscible Incompressible Fluids: A Reduction-Consistent and Thermodynamically-Consistent Formulation and Associated Algorithm
• Wall-bounded multiphase flows of N immiscible incompressible fluids: consistency and contact-angle boundary condition
• An Unconditionally Energy-Stable Scheme Based on an Implicit Auxiliary Energy Variable for Incompressible Two-Phase Flows with Different Densities Involving Only Precomputable Coefficient Matrices
• A consistent and conservative Phase-Field method for multiphase incompressible flows
• Multiphase Allen-Cahn and Cahn-Hilliard Models and Their Discretizations with the Effect of Pairwise Surface Tensions
• A consistent and conservative model and its scheme for $N$-phase-$M$-component incompressible flows
• A Method for Computing Inverse Parametric PDE Problems with Random-Weight Neural Networks
• A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models
• Phase-field-based lattice Boltzmann model for immiscible incompressible N-phase flows
• Multiphase Flows of N Immiscible Incompressible Fluids: An Outflow/Open Boundary Condition and Algorithm
• Phase Field Modeling and Numerical Algorithm for Two-Phase Dielectric Fluid Flows
• A Rotational Pressure-Correction Scheme for Incompressible Two-Phase Flows with Open Boundaries
• Fitted front tracking methods for two-phase incompressible Navier--Stokes flow: Eulerian and ALE finite element discretizations