Interface tracking characteristics of color-gradient lattice Boltzmann model for immiscible fluids

Generalized equilibria for color-gradient lattice Boltzmann model based on higher-order Hermite polynomials: A simplified implementation with central moments

We propose generalized equilibria of a three-dimensional color-gradient lattice Boltzmann model for two-component two-phase flows using higher-order Hermite polynomials. Although the resulting equilibrium distribution function, which includes a sixth-order term on the velocity, is computationally cumbersome, its equilibrium central moments (CMs) are velocity-independent and have a simplified form. Numerical experiments show that our approach, as in Wen et al. [Phys. Rev. E 100, 023301 (2019)2470-004510.1103/PhysRevE.100.023301] who consider terms up to third order, improves the Galilean invariance compared to that of the conventional approach. Dynamic problems can be solved with high accuracy at a density ratio of 10; however, the accuracy is still limited to a density ratio of 1000. For lower density ratios, the generalized equilibria benefit from the CM-based multiple-relaxation-time model, especially at very high Reynolds numbers, significantly improving the numerical stability.

Color-gradient lattice Boltzmann model for immiscible fluids with density contrast

We present a color-gradient-based lattice Boltzmann model for immiscible fluids with a large density contrast. The model employs the velocity-based equilibrium distribution function, initially proposed for the phase-field-based model by Zu and He [Phys. Rev. E 87, 043301 (2013)1539-375510.1103/PhysRevE.87.043301], with a modification necessary to satisfy the kinematic condition at the interface. Different from the existing color-gradient models, the present model allows to specify interface mobility that is independent of the fluid density ratio. Further, we provide a unified framework, which uses the recursive representation of the lattice Boltzmann equation, to derive the governing equations of the system. The emergent color dynamics thus obtained, through an analysis of the segregation operator, is shown to obey the locally conservative Allen-Cahn equation. We use a series of benchmarks, which include a stationary drop, a layered Poiseuille flow, translation of a drop under a forced velocity field, the Rayleigh-Taylor instability, and the capillary intrusion test to demonstrate the model's ability in dealing with complex flow problems.

Chromodynamic multirelaxation-time lattice Boltzmann scheme for fluids with density difference

We develop, after Dellar [Phys. Rev. E. 65, 036309 (2002)10.1103/PhysRevE.65.036309; J. Comput. Phys. 190, 351 (2003)10.1016/S0021-9991(03)00279-1], a multiple-relaxation-time (MRT), chromodynamic, multicomponent lattice Boltzmann equation (MCLBE) scheme for simulation of isothermal, immiscible fluid flow with a density contrast. It is based on Lishchuk's method [Brackbill, Kothe, and Zemach, J. Comp. Phys. 100, 335 (1992)10.1016/0021-9991(92)90240-Y; Lishchuk, Care, and Halliday, Phys. Rev. E. 67, 036701, (2003)10.1103/PhysRevE.76.036701] and the segregation of d'Ortona et al. [Phys. Rev. E. 51, 3718, (1995)10.1103/PhysRevE.51.3718]. We focus on fundamental model verifiability but do relate some of our data to that from previous approaches, due to Ba et al. [Phys. Rev. E 94, 023310 (2016)10.1103/PhysRevE.94.023310] and earlier Liu et al. [Phys. Rev. E 85, 046309 (2012)10.1103/PhysRevE.85.046309], who pioneered large density difference chromodynamic MCLBE and showed the practical benefits of an MRT collision model. Specifically, we test the extent to which chromodynamic MCLBE MRT schemes comply with the kinematic condition of mutual impenetrability and the continuous traction condition by developing analytical benchmarking flows. We conclude that our data, taken with those of Ba et al., verify the utility of MRT chromodynamic MCLBE.