Abstract:In this work a mass-conserved diffuse interface method is proposed for simulating incompressible flows of binary fluids with large density ratio. In the method, a mass correction term is introduced into the Cahn–Hilliard equation to compensate the mass losses or offset the mass increases caused by the numerical and modeling diffusion. Since the mass losses or increases are through the phase interfaces and at each time step, their values are very small, to keep mass conservation, mass sources or sinks are introduced and uniformly distributed in the volume of diffuse layer. With the uniform distribution, the mass correction term representing mass sources or sinks is derived analytically by applying mass conservation principle. By including the mass correction, the modified Cahn–Hilliard equation is solved by the fifth-order upwind scheme to capture the phase field of the bindery fluids. The flow field is simulated by the newly-developed multiphase lattice Boltzmann flux solver [20]. The proposed approach is validated by simulating the Laplace law, the merging of two bubbles, Rayleigh–Taylor instability and bubble rising under gravity with density ratio of 1000 and viscosity ratio of 100. Numerical results of interface shapes and flow properties agree well with both analytical solutions and benchmark data in the literature. Numerical results also show that the mass is well-conserved in all cases considered. In addition, it is demonstrated that the mass correction term at each time step is in the order of 10−4∼10−5, which is a small number compared with the magnitude of order parameter.

A low diffusion flux splitting method for inviscid compressible flows

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

Numerical Simulation of Incompressible Interfacial Flows by a Level Set Re-Distancing Method with Improved Mass Conservation

On the Accuracy and Robustness of A New Flux Splitting Method

Solution to Euler Equations by High‐resolution Upwind Compact Scheme Based on Flux Splitting

A New Splitting Method for Time-dependent Convection-dominated Diffusion Problems

An Effective Low Dissipation Method for Compressible Flows

New Splitting Methods for Convection-Dominated Diffusion Problems and Navier-Stokes Equations

A low-diffusion robust flux splitting scheme towards wide-ranging Mach number flows

An Interface-Compressed Diffuse Interface Method and Its Application for Multiphase Flows

A hybrid flux splitting method for compressible flow

A Three-Dimensional Explicit Sphere Function-Based Gas-Kinetic Flux Solver for Simulation of Inviscid Compressible Flows.

A Dimension Splitting Method for 3-D Incompressible Thermal Flow

Flux limiter lattice Boltzmann for compressible flows

A mass-conserved diffuse interface method and its application for incompressible multiphase flows with large density ratio

A Simplified Artificial Compressibility Flux for the Discontinuous Galerkin Solution of the Incompressible Navier-Stokes Equations

A Dissipation‐free Numerical Method to Solve One‐dimensional Hyperbolic Flow Equations

Meshfree Lattice Boltzmann Flux Solver for Compressible Inviscid Flows

An implicit kinetic inviscid flux for predicting continuum flows in all speed regimes

A Finite Volume Scheme For Weakly Compressible Equations, Part2: Application In Incompressible Flows

A HLLC-type Finite Volume Method for Incompressible Two-Phase Flows