# A fractional step lattice Boltzmann model for two phase flows with large density differences

## TL;DR

This paper introduces a fractional step lattice Boltzmann model for accurately simulating two-phase flows with large density differences, improving interface tracking and numerical stability.

## Contribution

The proposed method combines fractional step schemes with correction terms and advanced collision models to enhance accuracy and stability in high density ratio two-phase flow simulations.

## Key findings

- Successfully simulates flows with density ratios up to 1000.
- Accurately predicts interface deformation and flow dynamics.
- Validated against analytical solutions and literature data.

## Abstract

In this paper, a fractional step lattice Boltzmann method is proposed to model two-phase flows with large density differences by solving Cahn-Hilliard phase-field equation and the incompressible Navier-Stokes equations.In order to maintain a hyperbolic tangent property of the interface profile and conserve the volume, an interfacial profile correction term and a flux correction term are added into the original Cahn-Hilliard equation respectively. By using a fractional step scheme, the modified Cahn-Hilliard equation is split into two sub-equations. One is solved in the framework of lattice Boltzmann equation method. The other is solved by the finite difference method. Compared with the previous lattice Boltzmann methods, the proposed method is able to maintain the order parameter within a physically meaningful range, which is conductive to track the interface accurately. In addition, the multi-relaxation-time collision model and a high-order compact selective filter operation are employed to enhance the numerical stability. The proposed method can simulate two-phase fluid flows with the density ratio up to $1000$. In order to validate the accuracy and capability of the method, several benchmark problems, including single vortex deform of a circle, translation of a drop, Laplace-Young law, capillary wave and rising bubble with large density ratios, are presented. The results are in good agreement with the analytical solutions and the data in the literature for the investigated benchmarks.