Related Experiment Video
Updated: Jan 3, 2026

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
Kinematics of chromodynamic multicomponent lattice Boltzmann simulation with a large density contrast
K Burgin1, J Spendlove1, X Xu1,2
1Materials & Engineering Research Institute, Sheffield Hallam University, Howard Street, S1 1WB, United Kingdom.
Abstract:
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.
More Related Videos
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
11:03An Analog Macroscopic Technique for Studying Molecular Hydrodynamic Processes in Dense Gases and Liquids
Published on: December 4, 2017
Related Concept Videos
Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
Trends in Lattice Energy: Ion Size and Charge
Equilibrium Conditions for a Particle
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
First Law: Particles in Two-dimensional Equilibrium
Newton's first law tells us about...
Distribution of Molecular Speeds
First Law: Particles in One-dimensional Equilibrium