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)] and Ba et al. [Phys. Rev. E 94, 023310 (2016)], who advanced earlier work by Liu et al. [Phys. Rev. E 85, 046309 (2012)] 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)] and the free-energy scheme [see, e.g., Phys. Rev. E. 76, 045702(R) (2007), and references therein] which, in contradistinction to the present scheme and perhaps wisely, postulate appropriate kinematics a priori.

No abstract

We present a challenging validation of phase field multicomponent lattice Boltzmann equation (MCLBE) simulation against the Re=0 Stokes flow regime Taylor-Einstein theory of dilute suspension viscosity. By applying a number of recent advances in the understanding and the elimination of the interfacial microcurrent artefact, extending to a three-dimensional class of stability-enhancing multiple relaxation time collision models (which require no explicit collision matrix, note) and developing new interfacial interpolation schemes, we are able to obtain data that show that MCLBE may be applied in new flow regimes. Our data represent one of the most stringent tests yet attempted on LBE-one which received wisdom would preclude on grounds of overwhelming artefact flow.

An extended Benzi-Dellar lattice Boltzmann equation scheme [R. Benzi, S. Succi, and M. Vergassola, Europhys. Lett. 13, 727 (1990)EULEEJ0295-507510.1209/0295-5075/13/8/010; R. Benzi, S. Succi, and M. Vergassola, Phys. Rep. 222, 145 (1992)PRPLCM0370-157310.1016/0370-1573(92)90090-M; P. J. Dellar, Phys. Rev. E 65, 036309 (2002)1063-651X10.1103/PhysRevE.65.036309] is developed and applied to the problem of confirming, at low Re and drop fluid concentration, c, the variation of effective shear viscosity, η_{eff}=η_{1}[1+f(η_{1},η_{2})c], with respect to c for a sheared, two-dimensional, initially crystalline emulsion [here η_{1} (η_{2}) is the fluid (drop fluid) shear viscosity]. Data obtained with our enhanced multicomponent lattice Boltzmann method, using average shear stress and hydrodynamic dissipation, agree well once appropriate corrections to Landau's volume average shear stress [L. Landau and E. M. Lifshitz, Fluid Mechanics, 6th ed. (Pergamon, London, 1966)] are applied. Simulation results also confirm the expected form for f(η_{i},η_{2}), and they provide a reasonable estimate of its parameters. Most significantly, perhaps, the generality of our data supports the validity of Taylor's disputed simplification [G. I. Taylor, Proc. R. Soc. London, Ser. A 138, 133 (1932)1364-502110.1098/rspa.1932.0175] to reduce the effect of one hydrodynamic boundary condition (on the continuity of the normal contraction of stress) to an assumption that interfacial tension is sufficiently strong to maintain a spherical drop shape.

Emulsions are found throughout various industries including oil extraction, biological materials, and food products such as milk, condiments, and spreads. The study of their rheology is therefore important due to its impact on manufacturing efficiency and end product desirability. A key rheological measure is the emulsion viscosity, the fluid's resistance to flow, which affects the power required in production as well as the taste and texture. An emulsion's viscosity displays complex behaviour due to the droplet interfaces and interactions. Similarly, the sheared self-diffusion coefficient measures the amount of movement the droplets exhibit, due to the interactions between droplets. The presented mesoscopic lattice-Boltzmann models allow for these macroscopic properties to emerge from the simulations due to the explicit modelling of the droplets. A continuous surface force is applied to the lattice fluids to model droplet interfaces. The model is implemented in such a way as to allow the simulation of hundreds of droplets with limited computing power. The model is initially applied to a pipe flow, with the development of a pressure boundary condition. Boundary effects from the solid walls require their removal, using Lees-Edwards boundary conditions to represent bulk flow in a sheared system. The boundary conditions are extended to the multi-component flow, which allowed simulations to provide results for various emulsion systems with varying droplet concentrations, surface tensions, viscosity ratios, and shear rates. Trends and results from experimental and theoretical literature are recovered and constitutive models of emulsion viscosity have been evaluated. The agreement of these two dimensional lattice-Boltzmann models with three dimensional experimental results shows the usefulness of the method. The structure of the droplets and clustering behaviour they exhibit are examined and compared to solid particle suspension literature. Finally, the model is used in exploratory simulations to examine the effect of droplet bidispersity on the macroscopic properties; the witnessed effect agrees well with solid suspension literature. This mesoscopic model will allow for phenomenon on this scale to be more easily studied and may provide more accurate information for multi-scale analysis.