Modified phase-field-based lattice Boltzmann model for incompressible multiphase flows

Based on the phase-field theory, a multiple-relaxation-time (MRT) lattice Boltzmann model is proposed for the immiscible multiphase fluids. In this model, the local Allen-Chan equation is chosen as the target equation to capture the phase interface. Unlike previous MRT schemes, an off-diagonal relaxation matrix is adopted in the present model so that the target phase-field equation can be recovered exactly without any artificial terms. To check the necessity of removing those artificial terms, comparative studies were carried out among different MRT schemes with or without correction. Results show that the artificial terms can be neglected at low March number but will cause unphysical diffusion or interface undulation instability for the relatively large March number cases. The present modified model shows superiority in reducing numerical errors by adjusting the free parameters. As the interface transport coupled to the fluid flow, a pressure-evolution lattice Boltzmann equation is adopted for hydrodynamic properties. Several benchmark cases for multiphase flow were conducted to test the validity of the present model, including the static drop test, Rayleigh-Taylor instability, and single rising bubble test. For the rising bubble simulation at high density ratios, bubble dynamics obtained by the present modified MRT lattice Boltzmann model agree well with those obtained by the FEM-based level set and FEM-based phase-field models.