A phase-field model for large-density-ratio two-phase flows based on discrete unified gas-kinetic scheme

TL;DR

This paper introduces a phase-field model based on the discrete unified gas-kinetic scheme (DUGKS) for simulating large-density-ratio two-phase flows, demonstrating high accuracy and stability in complex interface dynamics.

Contribution

The paper develops a novel DUGKS-based phase-field model for multiphase flows with high density ratios, improving interface tracking and computational efficiency.

Findings

Achieved second-order convergence in interface translation tests.

Successfully simulated flows with density ratios up to 1000.

Predicted interface evolution accurately in Rayleigh-Taylor instability.

Abstract

In this paper, a phase-field based model under the framework of discrete unified gas-kinetic scheme (DUGKS) for incompressible multiphase fluid flows is proposed. Two kinetic models are constructed to solve the conservative Allen-Cahn (A-C) equation that accounts for the interface behavior and the incompressible hydrodynamic equations that govern the flow field, respectively. With a truncated equilibrium distribution function as well as a temporal derivative added to the source term, the macroscopic governing equations can be exactly recovered from the kinetic models through the Chapmann-Enskog analysis. Calculation of source terms involving high-order derivatives existed in the quasi-incompressible model is simplified. A series of benchmark cases including four interface-capturing tests and four binary flow tests are carried out. Results compared to that of lattice Boltzmann method (LBM) have been obtained. A convergence rate of second-order can be guaranteed in the test of interface diagonal translation. The capability of present method in interface tracking that undergoes a severe deformation has been verified. Stationary bubble and spinodal decomposition problems, both with a density ratio as high as 1000, are conducted and reliable solutions have been provided. The layered Poiseuille flow with a large viscosity ratio is simulated and numerical results agree well with the analytical solutions. Variation of positions of the bubble front and spike tip during the evolution of Rayleigh-Taylor instability (RTI) has been predicted precisely. However, the detailed depiction of complicated interface patterns appeared in the evolution process is failed, which is mainly caused by the relatively large numerical dissipation of DUGKS compared to that of LBM. A high-order DUGKS is needed to overcome this problem.