Abstract
This study firstly proposes a simple recursive method for deriving the macroscale equations from lattice Boltzmann equations. Similar to the Maxwell iteration based on the convective scaling, this method is used to expand the lattice Boltzmann (LB) equations with the time step \(\delta _{t}\). It is characterised by the incorporation of a nonequilibrium distribution function not appearing in the Maxwell iteration to considerably reduce the mathematical manipulations required. Next, we define the kinetic equations of a multicomponent (i.e. N-component) system based on a model using the Maxwell velocity distribution law for the equilibrium distribution function appearing in the cross-collision terms. Then, using this simple recursive method, we derive the generalized Stefan–Maxwell equation, which is the macroscale governing equation of a multicomponent system while ensuring the mass conservation. In short, our objective is to firstly define the kinetic equations of a multi-component system having a clear physical interpretation and then formulate the LB equations of any N-component system deductively.
Similar content being viewed by others
References
Ryan, E.M., Mukherjee, P.P.: Mesoscale modeling in electrochemical devices-a critical perspective. Prog. Energy Combust. Sci. 71, 118 (2019)
Jiang, Z.Y., Qu, Z.G., Zhou, L.: Lattice Boltzmann simulation of ion and electron transport during the discharge process in a randomly reconstructed porous electrode of a lithium-ion battery. Int. J. Heat Mass Transf. 123, 500 (2018)
Ishikawa, H., Sugawara, Y., Inoue, G., Kawase, M.: Effects of Pt and ionomer ratios on the structure of catalyst layer: a theoretical model for polymer electrolyte fuel cells. J. Power Sources 374, 196 (2018)
He, A., Gong, J., Shikazono, N.: Three dimensional electrochemical simulation of solid oxide fuel cell cathode based on microstructure reconstructed by marching cubes method. J. Power Sources 385, 91 (2018)
Chen, S., Doolen, G.D.: lattice Boltzmann method for fluid flows. Annu. Rev. Fluid Mech. 30, 329 (1998)
Fares, E.: Unsteady flow simulation of the Ahmed reference body using a lattice Boltzmann approach. Comput. Fluids 35, 940 (2006)
Luo, L.S., Girimaji, S.S.: Theory of the lattice Boltzmann method: two-fluid model for binary mixtures. Phys. Rev. 67, 036302 (2003)
Gross, E.P., Krook, M.: Model for collision processes in gases: small-amplitude oscillations of charged two-component systems. Phys. Rev. E 102, 593 (1956)
Krishna, R., Wesselingh, J.A.: The Maxwell–Stefan approach to mass transfer. Chem. Eng. Sci. 52, 861 (1997)
Tong, Z.X., He, Y.L., Chen, L., Xie, T.: A multi-componet lattice Boltzmann method in consistent with Stefan–Maxwell equations. Comput. Fluids 105, 155 (2014)
Shan, X., Doolen, G.: Multicomponent lattice-Boltzmann model with interparticle interaction. J. Stat. Phys. 81, 379 (1995)
Chai, Z., Zhao, T.S.: A pseudopotential-based multiple-relaxation-time lattice Boltzmann model for multicomponent/multiphase flows. Acta. Mech. Sin. 28, 983 (2012)
Swift, M.R., Orlandini, E., Osborn, W.R., Yeomans, J.M.: Lattice Boltzmann simulations of liquid–gas and binary fluid systems. Phys. Rev. E 54, 5041 (1996)
Ridl, K.S., Wagner, A.J.: Lattice Boltzmann simulation of mixtures with multicomponent van der Waals equation of state. Phys. Rev. E 98, 043305 (2018)
Hosseini, S.A., Darabiha, N., Thévenin, D.: Mass-conserving advection-diffusion Lattice Boltzmann model for multi-species reacting flows. Physica A 499, 40 (2018)
Chai, Z., Guo, X., Wang, L., Shi, B.: Maxwell-Stefan-theory-based lattice Boltzmann model for diffusion in multicomponent mixtures. Phys. Rev. E 99, 023312 (2019)
He, X., Luo, L.S.: Theory of the lattice Boltzmann method: from the Boltzmann equation to the lattice Boltzmann equation. Phys. Rev. E 56, 6811 (1997)
Ansumali, S., Karlin, I.V., Arcidiacono, S., Abbas, A., Prasianakis, N.I.: Hydrodynamics beyond Navier–Stokes: exact Solution to the Lattice Boltzmann Hierarchy. Phys. Rev. Lett. 98, 124502 (2007)
Chapman, S., Cowling, T.G.: Theory of Non-uniform Gases, 3rd edn. Cambridge University Press, Cambridge (1970)
Zhao, W., Yong, W.A.: Maxwell iteration for the lattice Boltzmann method with diffusive scaling. Phys. Rev. E 95, 033311 (2017)
Vienne, L., Marié, S., Grasso, F.: Lattice Boltzmann method for miscible gases: a forcing-term approach. Phys. Rev. 100, 023309 (2019)
Krishna, R.: Uphill diffusion in multicomponent mixtures. Chem. Soc. Rev. 44, 2812 (2015)
Hou, S., Zou, Q., Chen, S., Cogley, A.C.: Simulation of Cavity Flow by the Lattice Boltzmann Method. J. Comput. Phys. 118, 329 (1995)
Garzó, V., Santos, A., Brey, J.J.: A kinetic model for a multicomponent gas. Phys. Fluids A 1, 380 (1989)
Sirovich, L.: Kinetic modeling of gas mixtures. Phys. Fluids 5, 908 (1962)
Jianfei, X., Matthew, K.B., Livio, G., Oliver, L.H., Duncan, A., Jason, M.R.: Effective mean free path and viscosity of confined gases. Phys. Fluids 31, 072002 (2019)
Joshi, A.S., Peracchio, A.A., Grew, K.N., Chiu, W.K.S.: Lattice Boltzmann method for continuum, multi-component mass diffusion in complex 2D geometries. J. Phys. D 40, 2961 (2007)
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Shin-ichi Sasa.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A
Appendix A proves that coefficients of the same order must be equal on both sides for (6) to hold true for any \(\delta _{t}\). The zero-th order terms of (6) are gathered on the left-hand side and all other terms are gathered on the right-hand side. The absolute value is taken on both sides, and the right-hand side is evaluated with an inequality sign.
Here, when \(\delta _{t} \rightarrow 0\), the right-hand side converges to 0. Therefore, the necessary conditions for (6) to hold true for any \(\delta _{t}\) can be obtained as:
In addition, (81) is substituted in (6), both sides are divided by \(\delta _{t}\), and an evaluation similar to (80) is performed to obtain (82).
The proof is completed by repeating this operation in the same manner.
Appendix B
Appendix B shows the detailed derivation of the mass conservation law, the Euler equation and the Navier–Stokes equation from (21) to (29).
First, the mass conservation law is derived. Concerning (8), the sum is taken with i as below.
The lef-hand side of (83) is transformed as:
(84) is substituted into (83), the mass conservation law (85) is derived.
Next, the Euler equation is derived. Concerning (8), both sides of (8) are multiplied by the discrete velocity and the sum is taken with i as below.
The left-hand side of (86) is transformed as:
(87) is substituted into (86), and using (85) the Euler equation (88) is derived.
Finally, the Navier–Stokes equation is derived. As the derivation process up to (21) have been shown, we begin the calculation from that equation below.
Here, the Euler equation (88) is substituted into the above, the right-hand side is further transformed as:
To continue the calculations, it is necessary to use the equation of state (23) and the equilibrium distribution function, so we take the equilibrium distribution function to be the conventional expression [26] :
where, \(\sigma \) represents the type of the lattice velocity vector.
(23) and (90) are substituted into the right-hand side of (89) and the time derivative is rewritten as the space derivative using the mass conservation law (85), the right-hand side of the equation is transformed as:
Here, the last term tensor in (91), is calculated about the D2Q9, for example, as follows:
By setting \(\frac{2{{c}^{4}}}{3}{{A}_{1}}{{B}_{1}}+\frac{4{{c}^{4}}}{3}{{A}_{2}}{{B}_{2}} =a \rho \), (92) is calculated as below:
Substituting this calculated tensor into (91) yields,
Then, by setting \(a=C_{s}^{2}\), we confirm that the the Navier–Stokes equation (28) can be derived.
By the way, about the D2Q9 equilibrium distribution functions (25), we can confirm \(a=C_{s}^{2}\) directly. From the equilibrium distribution functions (25), \({{A}_{1}}{{B}_{1}}\), \({{A}_{2}}{{B}_{2}}\) are defined by (95), (96).
From these relations, \(a=C_{s}^{2}\) is confirmed as below.
As a matter of fact, (97) is a necessary condition in order to construct the equilibrium distribution functions [26].
Rights and permissions
About this article
Cite this article
Yamamoto, K., Seta, T. Derivation of Multicomponent Lattice Boltzmann Equations by Introducing a Nonequilibrium Distribution Function into the Maxwell Iteration Based on the Convective Scaling. J Stat Phys 182, 4 (2021). https://doi.org/10.1007/s10955-020-02686-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10955-020-02686-x