Benchmarking of three-dimensional multicomponent lattice Boltzmann equation

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.

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

We study the interface tracking characteristics of a color-gradient-based lattice Boltzmann model for immiscible flows. Investigation of the local density change in one of the fluid phases, via a Taylor series expansion of the recursive lattice Boltzmann equation, leads to the evolution equation of the order parameter that differentiates the fluids. It turns out that this interface evolution follows a conservative Allen-Cahn equation with a mobility which is independent of the fluid viscosities and surface tension. The mobility of the interface, which solely depends upon lattice speed of sound, can have a crucial effect on the physical dynamics of the interface. Further, we find that, when the equivalent lattice weights inside the segregation operator are modified, the resulting differential operators have a discretization error that is anisotropic to the leading order. As a consequence, the discretization errors in the segregation operator, which ensures a finite interface width, can act as a source of the spurious currents. These findings are supported with the help of numerical simulations.

Kinematics of chromodynamic multicomponent lattice Boltzmann simulation with a large density contrast

The utility of an enhanced chromodynamic, color gradient or phase-field multicomponent lattice Boltzmann (MCLB) equation for immiscible fluids with a density difference was demonstrated by Wen et al. [Phys. Rev. E 100, 023301 (2019)2470-004510.1103/PhysRevE.100.023301] and Ba et al. [Phys. Rev. E 94, 023310 (2016)2470-004510.1103/PhysRevE.94.023310], who advanced earlier work by Liu et al. [Phys. Rev. E 85, 046309 (2012)PLEEE81539-375510.1103/PhysRevE.85.046309] by removing certain error terms in the momentum equations. But while these models' collision scheme has been carefully enhanced by degrees, there is, currently, no quantitative consideration in the macroscopic dynamics of the segregation scheme which is common to all. Here, by analysis of the kinetic-scale segregation rule (previously neglected when considering the continuum behavior) we derive, bound, and test the emergent kinematics of the continuum fluids' interface for this class of MCLB, concurrently demonstrating the circular relationship with-and competition between-the models' dynamics and kinematics. The analytical and numerical results we present in Sec. V confirm that, at the kinetic scale, for a range of density contrast, color is a material invariant. That is, within numerical error, the emergent interface structure is isotropic (i.e., without orientation dependence) and Galilean-invariant (i.e., without dependence on direction of motion). Numerical data further suggest that reported restrictions on the achievable density contrast in rapid flow, using chromodynamic MCLB, originate in the effect on the model's kinematics of the terms deriving from our term F_{1i} in the evolution equation, which correct its dynamics for large density differences. Taken with Ba's applications and validations, this result significantly enhances the theoretical foundation of this MCLB variant, bringing it somewhat belatedly further into line with the schemes of Inamuro et al. [J. Comput. Phys. 198, 628 (2004)JCTPAH0021-999110.1016/j.jcp.2004.01.019] and the free-energy scheme [see, e.g., Phys. Rev. E. 76, 045702(R) (2007)10.1103/PhysRevE.76.045702, and references therein] which, in contradistinction to the present scheme and perhaps wisely, postulate appropriate kinematics a priori.