Abstract:We present the detailed analysis of the diffusive transport of spatially inhomogeneous fluid mixtures and the interplay between structural and dynamical properties varying on the atomic scale. The present treatment is based on different areas of liquid state theory, namely kinetic and density functional theory and their implementation as an effective numerical method via the Lattice Boltzmann approach. By combining the first two methods it is possible to obtain a closed set of kinetic equations for the singlet phase space distribution functions of each species. The interactions among particles are considered within a self-consistent approximation and the resulting effective molecular fields are analyzed. We focus on multispecies diffusion in systems with short-range hard-core repulsion between particles of unequal sizes and weak attractive long-range interactions. As a result, the attractive part of the potential does not contribute explicitly to viscosity but to diffusivity and the thermodynamic properties. Finally, we obtain a practical scheme to solve the kinetic equations by employing a discretization procedure derived from the Lattice Boltzmann approach. Within this framework, we present numerical data concerning the mutual diffusion properties both in the case of a quiescent bulk fluid and shear flow inducing Taylor dispersion.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the understanding and modeling of diffusion behavior in multi - component nanosystems. Specifically, the authors focus on the diffusion - transport characteristics of spatially inhomogeneous fluid mixtures and the interaction of their structural and kinetic properties at the atomic scale. The main challenge of the research lies in how to describe the diffusion behavior of complex fluid systems in combination with microscopic - structure details while maintaining the feasibility of numerical calculations.
### Research Background
Modern scientific, medical, and technological applications require a better understanding of the molecular mechanisms of liquids near solid substrates or at interfaces. The structures and transport characteristics of highly - confined fluids or fluids near free surfaces are different from their bulk behaviors because these systems have a large surface - to - volume ratio. Many phenomena occurring at the molecular scale, such as diffusion, mixing, shear - thinning, and lane formation, involve the interaction of microstructure and transport characteristics, which requires the study of long - time flow behaviors. Since the molecular - dynamics method is computationally very expensive, alternative methods need to be explored.
### Research Methods
The authors adopted several theoretical frameworks to solve this problem:
1. **Kinetic theory**: Extend the Boltzmann equation to the dense - fluid region to describe the thermodynamic and hydrodynamic behaviors of simple fluids.
2. **Density - functional theory (DDFT)**: Assume that the evolution of the system is driven by "thermodynamic forces", that is, the functional derivative of the free - energy functional with respect to the local density.
3. **Lattice - Boltzmann method (LBM)**: Achieve efficient numerical solutions by discretizing position and velocity degrees of freedom.
### Main Contributions
- Proposed a detailed analysis of multi - component fluid diffusion, especially the roles of short - range hard - core repulsion and weak long - range attraction.
- Derived a closed - form kinetic equation applicable to multi - component fluids and achieved numerical solutions through LBM.
- Verified theoretical predictions through numerical simulations, especially the mutual - diffusion characteristics in the case of stationary bulk fluids and shear flows inducing Taylor dispersion.
### Mathematical Formulas
Several key formulas are used in the paper to describe the kinetic behavior of the system, for example:
- Transport equation for the single - particle phase - space distribution function:
\[
\frac{\partial f_{\alpha}(\mathbf{r}, \mathbf{v}, t)}{\partial t}+\mathbf{v} \cdot \nabla f_{\alpha}(\mathbf{r}, \mathbf{v}, t)+\frac{\mathbf{F}_{\alpha}(\mathbf{r})}{m_{\alpha}} \cdot \frac{\partial f_{\alpha}(\mathbf{r}, \mathbf{v}, t)}{\partial \mathbf{v}}=-\omega\left[f_{\alpha}(\mathbf{r}, \mathbf{v}, t)-\psi_{\alpha}^{\perp}(\mathbf{r}, \mathbf{v}, t)\right]+\frac{\Phi_{\alpha}(\mathbf{r}, t)}{k_{B} T} \cdot(\mathbf{v}-\mathbf{u}(\mathbf{r}, t)) \psi_{\alpha}(\mathbf{r}, \mathbf{v}, t)
\]
- Conservation equation for local concentration:
\[
\frac{\partial c(\mathbf{r}, t)}{\partial t}+\mathbf{u}(\mathbf{r}, t) \cdot \nabla c(\mathbf{r}, t)+\frac{1}{\rho} \nabla \cdot[\rho(\mathbf{r}, t) c(\mathbf{r}, t)(1 - c(\mathbf{r}, t)) \mathbf{w}(\mathbf{r}, t)] = 0
\]
- Expression for the diffusion coefficient \(D_{AB}\):
\[
D_{AB}=\frac{3}{8n}(k_{B}