Lattice Boltzmann model for the simulation of multicomponent mixtures

Lattice Boltzmann modeling and simulation of velocity and concentration slip effects on the catalytic reaction rate of strongly nonequimolar reactions in microflows

A lattice Boltzmann (LB) interfacial gas-solid three-dimensional model is developed for isothermal multicomponent flows with strongly nonequimolar catalytic reactions, further accounting for the presence of velocity slips and concentration jumps. The model includes diffusion coefficients of all reactive species in the calculation of the catalytic reaction rates as well as an updated velocity at the reactive boundary node. Lattice Boltzmann simulations are performed in a catalytic channel-flow geometry under a wide range of Knudsen (Kn) and surface Damköhler (Da_{s}) numbers. Comparisons with simulations from a computational fluid dynamics (CFD) Navier-Stokes solver show good agreement in the continuum regime (Kn<0.01) in terms of flow velocity and reactive species distributions, while comparisons with direct simulation Monte Carlo results from the literature attest to the model's applicability in capturing the correct slip velocity at Kn as high as 0.1, even with a significantly reduced number of grid points (N=10) in the cross-flow direction. Theoretical and numerical results demonstrate that the term Da_{s}×Kn×A_{2} (where A_{2} is a function of the mass accommodation coefficient) determines the significance of the concentration jump on the catalytic reaction rate. The developed model is applicable for many catalytic microflow systems with complex geometries (such as reactors with porous networks) and large velocity/concentration slips (such as catalytic microthrusters for space applications).

Two-fluid kinetic theory for dilute polymer solutions

We provide a Boltzmann-type kinetic description for dilute polymer solutions based on two-fluid theory. This Boltzmann-type description uses a quasiequilibrium based relaxation mechanism to model collisions between a polymer dumbbell and a solvent molecule. The model reproduces the desired macroscopic equations for the polymer-solvent mixture. The proposed kinetic scheme leads to a numerical algorithm which is along the lines of the lattice Boltzmann method. Finally, the algorithm is applied to describe the evolution of a perturbed Kolmogorov flow profile, whereby we recover the major elastic effect exhibited by a polymer solution, specifically, the suppression of the original inertial instability.

Lattice Boltzmann model with generalized wall boundary conditions for arbitrary catalytic reactivity

A lattice Boltzmann model for multispecies flows with catalytic reactions is developed, which is valid from very low to very high surface Damköhler numbers (Da_{s}). The previously proposed model for catalytic reactions [S. Arcidiacono, J. Mantzaras, and I. V. Karlin, Phys. Rev. E 78, 046711 (2008)PLEEE81539-375510.1103/PhysRevE.78.046711], which is applicable for low-to-moderate Da_{s} and encompasses part of the mixed kinetics and transport-controlled regime, is revisited and extended for the simulation of arbitrary kinetics-to-transport rate ratios, including strongly transport-controlled conditions (Da_{s}→∞). The catalytic boundary condition is modified by bringing nonlocal information on the wall reactive nodes, allowing accurate evaluation of chemical rates even when the concentration of the deficient reactant at the wall becomes vanishingly small. The developed model is validated against a finite volume Navier-Stokes CFD (Computational Fluid Dynamics) solver for the total oxidation of methane in an isothermal channel-flow configuration. CFD simulations and lattice Boltzmann simulations with the old and new catalytic reaction models are compared against each other. The new model demonstrates a second order accuracy in space and time and provides accurate results at very high Da_{s} (∼10^{9}) where the old model fails. Moreover, to achieve the same accuracy at moderate-to-high Da_{s} of O(1), the new model requires ∼2^{d}×N coarser grid than the original model, where d is the spatial dimension and N the number of species.

Multiple-relaxation-time discrete Boltzmann modeling of multicomponent mixture with nonequilibrium effects

A multiple-relaxation-time discrete Boltzmann model (DBM) is proposed for multicomponent mixtures, where compressible, hydrodynamic, and thermodynamic nonequilibrium effects are taken into account. It allows the specific heat ratio and the Prandtl number to be adjustable, and is suitable for both low and high speed fluid flows. From the physical side, besides being consistent with the multicomponent Navier-Stokes equations, Fick's law, and Stefan-Maxwell diffusion equation in the hydrodynamic limit, the DBM provides more kinetic information about the nonequilibrium effects. The physical capability of DBM to describe the nonequilibrium flows, beyond the Navier-Stokes representation, enables the study of the entropy production mechanism in complex flows, especially in multicomponent mixtures. Moreover, the current kinetic model is employed to investigate nonequilibrium behaviors of the compressible Kelvin-Helmholtz instability (KHI). The entropy of mixing, the mixing area, the mixing width, the kinetic and internal energies, and the maximum and minimum temperatures are investigated during the dynamic KHI process. It is found that the mixing degree and fluid flow are similar in the KHI process for cases with various thermal conductivity and initial temperature configurations, while the maximum and minimum temperatures show different trends in cases with or without initial temperature gradients. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI morphological structure.

Multiple-relaxation-time lattice Boltzmann method for the Navier-Stokes and nonlinear convection-diffusion equations: Modeling, analysis, and elements

In this paper, we first present a unified framework of multiple-relaxation-time lattice Boltzmann (MRT-LB) method for the Navier-Stokes and nonlinear convection-diffusion equations where a block-lower-triangular-relaxation matrix and an auxiliary source distribution function are introduced. We then conduct a comparison of the four popular analysis methods (Chapman-Enskog analysis, Maxwell iteration, direct Taylor expansion, and recurrence equations approaches) that have been used to obtain the macroscopic Navier-Stokes and nonlinear convection-diffusion equations from the MRT-LB method and show that from mathematical point of view, these four analysis methods can give the same equations at the second-order of expansion parameters. Finally, we give some elements that are needed in the implementation of the MRT-LB method and also find that some available LB models can be obtained from this MRT-LB method.

Lattice Boltzmann method for miscible gases: A forcing-term approach

A lattice Boltzmann method for miscible gases is presented. In this model, the standard lattice Boltzmann method is employed for each species composing the mixture. Diffusion interaction among species is taken into account by means of a force derived from kinetic theory of gases. Transport coefficients expressions are recovered from the kinetic theory. Species with dissimilar molar masses are simulated by also introducing a force. Finally, mixing dynamics is recovered as shown in different applications: an equimolar counterdiffusion case, Loschmidt's tube experiment, and an opposed jets flow simulation. Since collision is not altered, the present method can easily be introduced in any other lattice Boltzmann algorithms.

Maxwell-Stefan-theory-based lattice Boltzmann model for diffusion in multicomponent mixtures

The phenomena of diffusion in multicomponent (more than two components) mixtures are universal in both science and engineering, and from the mathematical point of view, they are usually described by the Maxwell-Stefan (MS)-theory-based diffusion equations where the molar average velocity is assumed to be zero. In this paper, we propose a multiple-relaxation-time lattice Boltzmann (LB) model for the mass diffusion in multicomponent mixtures and also perform a Chapman-Enskog analysis to show that the MS continuum equations can be correctly recovered from the developed LB model. In addition, considering the fact that the MS-theory-based diffusion equations are just a diffusion type of partial differential equations, we can also adopt much simpler lattice structures to reduce the computational cost of present LB model. We then conduct some simulations to test this model and find that the results are in good agreement with the previous work. Besides, the reverse diffusion, osmotic diffusion, and diffusion barrier phenomena are also captured. Finally, compared to the kinetic-theory-based LB models for multicomponent gas diffusion, the present model does not include any complicated interpolations, and its collision process can still be implemented locally. Therefore, the advantages of single-component LB method can also be preserved in present LB model.

In silico modeling of magnetic resonance flow imaging in complex vascular networks

The paper presents a computational model of magnetic resonance (MR) flow imaging. The model consists of three components. The first component is used to generate complex vascular structures, while the second one provides blood flow characteristics in the generated vascular structures by the lattice Boltzmann method. The third component makes use of the generated vascular structures and flow characteristics to simulate MR flow imaging. To meet computational demands, parallel algorithms are applied in all the components. The proposed approach is verified in three stages. In the first stage, experimental validation is performed by an in vitro phantom. Then, the simulation possibilities of the model are shown. Flow and MR flow imaging in complex vascular structures are presented and evaluated. Finally, the computational performance is tested. Results show that the model is able to reproduce flow behavior in large vascular networks in a relatively short time. Moreover, simulated MR flow images are in accordance with the theoretical considerations and experimental images. The proposed approach is the first such an integrative solution in literature. Moreover, compared to previous works on flow and MR flow imaging, this approach distinguishes itself by its computational efficiency. Such a connection of anatomy, physiology and image formation in a single computer tool could provide an in silico solution to improving our understanding of the processes involved, either considered together or separately.

Thermal multicomponent lattice Boltzmann model for catalytic reactive flows

Catalytic reactions are of great interest in many applications related to power generation, fuel reforming and pollutant abatement, as well as in various biochemical processes. A recently proposed lattice Boltzmann model for thermal binary-mixture gas flows [J. Kang, N. I. Prasianakis, and J. Mantzaras, Phys. Rev. E. 87, 053304 (2013)] is revisited and extended for the simulation of multispecies flows with catalytic reactions. The resulting model can handle flows with large temperature and concentration gradients. The developed model is presented in detail and validated against a finite volume Navier-Stokes solver in the case of channel-flow methane catalytic combustion. The surface chemistry is treated with a one-step global reaction for the catalytic total oxidation of methane on platinum. In order to take into account thermal effects, the catalytic boundary condition of S. Arcidiacono, J. Mantzaras, and I. V. Karlin [Phys. Rev. E 78, 046711 (2008)] is adapted to account for temperature variations. Speed of sound simulations further demonstrate the physical integrity and unique features of the model.

Lattice Boltzmann scheme for electrolytes by an extended Maxwell-Stefan approach

This paper presents an extended multicomponent lattice Boltzmann model for the simulation of electrolytes. It is derived by means of a finite discrete velocity model and its discretization. The model recovers momentum and mass transport according to the incompressible Navier-Stokes equation and Maxwell-Stefan formulation, respectively. It includes external driving forces (e.g., electric field) on diffusive and viscous scales, concentration-dependent Maxwell-Stefan diffusivities, and thermodynamic factors. The latter take into account nonideal diffusion behavior, which is essential as electrolytes involve charged species and therefore nonideal long and short-range interactions among the molecules of the species. Furthermore, we couple our scheme to a finite element method to include electrostatic interactions on the macroscopic level. Numerical experiments show the validity of the presented model.

Lattice-Boltzmann method for the simulation of multiphase mass transfer and reaction of dilute species

Despite the popularity of the lattice-Boltzmann method (LBM) in simulating multiphase flows, a general approach for modeling dilute species in multiphase systems is still missing. In this report we propose to modify the collision operator of the solute by introducing a modified redistribution scheme. This operator is based on local fluid variables and keeps the parallelism inherent to LBM. After deriving macroscopic transport equations, an analytical equation of state of the solute is exhibited and the method is proven constituting a unified framework to simulate arbitrary solute distribution between phases, including single-phase soluble compounds, amphiphilic species with a partition coefficient, and surface-adsorbed compounds.

Interaction pressure tensor for a class of multicomponent lattice Boltzmann models

We present a theory to obtain the pressure tensor for a class of nonideal multicomponent lattice Boltzmann models, thus extending the theory presented by X. Shan [Phys. Rev. E 77, 066702 (2008)] for single-component fluids. We obtain the correct form of the pressure tensor directly on the lattice and the resulting equilibrium properties are shown to agree very well with those measured from numerical simulations. Results are compared with those of alternative theories.

Lattice Boltzmann model for thermal binary-mixture gas flows

A lattice Boltzmann model for thermal gas mixtures is derived. The kinetic model is designed in a way that combines properties of two previous literature models, namely, (a) a single-component thermal model and (b) a multicomponent isothermal model. A comprehensive platform for the study of various practical systems involving multicomponent mixture flows with large temperature differences is constructed. The governing thermohydrodynamic equations include the mass, momentum, energy conservation equations, and the multicomponent diffusion equation. The present model is able to simulate mixtures with adjustable Prandtl and Schmidt numbers. Validation in several flow configurations with temperature and species concentration ratios up to nine is presented.

High-order thermal lattice Boltzmann models derived by means of Gauss quadrature in the spherical coordinate system

We use the spherical coordinate system in the momentum space and an appropriate discretization procedure to derive a hierarchy of lattice Boltzmann (LB) models with variable temperature. The separation of the integrals in the momentum space into angular and radial parts allows us to compute the moments of the equilibrium distribution function by means of Gauss-Legendre and Gauss-Laguerre quadratures, as well as to find the elements of the discrete momentum set for each LB model in the hierarchy. The capability of the high-order models in this hierarchy to capture specific effects in microfluidics is investigated through a computer simulation of Couette flow by using the Shakhov collision term to get the right value of the Prandtl number.

Lattice Boltzmann study of pattern formation in reaction-diffusion systems

Pattern formation in reaction-diffusion systems is of great importance in surface micropatterning [Grzybowski et al., Soft Matter 1, 114 (2005)], self-organization of cellular micro-organisms [Schulz et al., Annu. Rev. Microbiol. 55, 105 (2001)], and in developmental biology [Barkai et al., FEBS Journal 276, 1196 (2009)]. In this work, we apply the lattice Boltzmann method to study pattern formation in reaction-diffusion systems. As a first methodological step, we consider the case of a single species undergoing transformation reaction and diffusion. In this case, we perform a third-order Chapman-Enskog multiscale expansion and study the dependence of the lattice Boltzmann truncation error on the diffusion coefficient and the reaction rate. These findings are in good agreement with numerical simulations. Furthermore, taking the Gray-Scott model as a prominent example, we provide evidence for the maturity of the lattice Boltzmann method in studying pattern formation in nonlinear reaction-diffusion systems. For this purpose, we perform linear stability analysis of the Gray-Scott model and determine the relevant parameter range for pattern formation. Lattice Boltzmann simulations allow us not only to test the validity of the linear stability phase diagram including Turing and Hopf instabilities, but also permit going beyond the linear stability regime, where large perturbations give rise to interesting dynamical behavior such as the so-called self-replicating spots. We also show that the length scale of the patterns may be tuned by rescaling all relevant diffusion coefficients in the system with the same factor while leaving all the reaction constants unchanged.

Multicomponent lattice Boltzmann model from continuum kinetic theory

We derive from the continuum kinetic theory a multicomponent lattice Boltzmann model with intermolecular interaction. The resulting model is found to be consistent with the model previously derived from a lattice-gas cellular automaton [X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993)] but applies in a much broader domain. A number of important insights are gained from the kinetic theory perspective. First, it is shown that even in the isothermal case, the energy equipartition principle dictates the form of the equilibrium distribution function. Second, thermal diffusion is shown to exist and the corresponding diffusivities are given in terms of macroscopic parameters. Third, the ordinary diffusion is shown to satisfy the Maxwell-Stefan equation at the ideal-gas limit.

Finite-difference-based multiple-relaxation-times lattice Boltzmann model for binary mixtures

In this paper, we propose a finite-difference-based lattice Boltzmann equation (LBE) model with multiple-relaxation times (MRT), in which the distribution functions of individual species evolve on a same regular lattice without any interpolations. Furthermore, the use of the MRT enables the model more flexible so that it can be applied to mixtures of species with different viscosities and adjustable Schmidt number. Some numerical tests are conducted to validate the model, the numerical results are found to agree well with analytical solutions/or other numerical results, and good numerical stability of the proposed LBE model is also observed.

Three ways to lattice Boltzmann: a unified time-marching picture

It is shown that the lattice Boltzmann equation (LBE) corresponds to an explicit Verlet time-marching scheme for a continuum generalized Boltzmann equation with a memory delay equal to a half time step. This proves second-order accuracy of LBE with respect to this generalized equation, with no need of resorting to any implicit time-marching procedure (Crank-Nicholson) and associated nonlinear variable transformations. It is also shown, and numerically demonstrated, that this equivalence is not only formal, but it also translates into a complete equivalence of the corresponding computational schemes with respect to the hydrodynamic equations. Second-order accuracy with respect to the continuum kinetic equation is also numerically demonstrated for the case of the Taylor-Green vortex. It is pointed out that the equivalence is however broken for the case in which mass and/or momentum are not conserved, such as for chemically reactive flows and mixtures. For such flows, the time-centered implicit formulation may indeed offer a better numerical accuracy.

Lattice Boltzmann scheme for mixture modeling: analysis of the continuum diffusion regimes recovering Maxwell-Stefan model and incompressible Navier-Stokes equations

A finite difference lattice Boltzmann scheme for homogeneous mixture modeling, which recovers Maxwell-Stefan diffusion model in the continuum limit, without the restriction of the mixture-averaged diffusion approximation, was recently proposed [P. Asinari, Phys. Rev. E 77, 056706 (2008)]. The theoretical basis is the Bhatnagar-Gross-Krook-type kinetic model for gas mixtures [P. Andries, K. Aoki, and B. Perthame, J. Stat. Phys. 106, 993 (2002)]. In the present paper, the recovered macroscopic equations in the continuum limit are systematically investigated by varying the ratio between the characteristic diffusion speed and the characteristic barycentric speed. It comes out that the diffusion speed must be at least one order of magnitude (in terms of Knudsen number) smaller than the barycentric speed, in order to recover the Navier-Stokes equations for mixtures in the incompressible limit. Some further numerical tests are also reported. In particular, (1) the solvent and dilute test cases are considered, because they are limiting cases in which the Maxwell-Stefan model reduces automatically to Fickian cases. Moreover, (2) some tests based on the Stefan diffusion tube are reported for proving the complete capabilities of the proposed scheme in solving Maxwell-Stefan diffusion problems. The proposed scheme agrees well with the expected theoretical results.

Lattices for the lattice Boltzmann method

A recently introduced theory of higher-order lattice Boltzmann models [Chikatamarla and Karlin, Phys. Rev. Lett. 97, 190601 (2006)] is elaborated in detail. A general theory of the construction of lattice Boltzmann models as an approximation to the Boltzmann equation is presented. New lattices are found in all three dimensions and are classified according to their accuracy (degree of approximation of the Boltzmann equation). The numerical stability of these lattices is argued based on the entropy principle. The efficiency and accuracy of many new lattices are demonstrated via simulations in all three dimensions.

Lattice Boltzmann equation for microscale gas flows of binary mixtures

Modeling and simulating gas flows in and around microdevices are a challenging task in both science and engineering. In practical applications, a gas is usually a mixture made of different components. In this paper we propose a lattice Boltzmann equation (LBE) model for microscale flows of a binary mixture based on a recently developed LBE model for continuum mixtures [P. Asinari and L.-S. Luo, J. Comput. Phys. 227, 3878 (2008)]. A consistent boundary condition for gas-solid interactions is proposed and analyzed. The LBE is validated and compared with theoretical results or other reported data. The results show that the model can serve as a potential method for flows of binary mixture in the microscale.

Lattice Boltzmann modeling of multicomponent diffusion in narrow channels

We investigate lattice Boltzmann (LB) modeling of multicomponent diffusion for finite Knudsen numbers. Analytic solutions for binary diffusion in narrow channels, where both molecular and Knudsen diffusion are of importance, are obtained for the standard and higher-order LB methods and validated against the results from the direct simulation Monte Carlo (DSMC) method. The LB methods are shown to reproduce the diffusion slip phenomena. In the DSMC method, while fluid particles are diffusely reflected on a wall, significant component slip and a kinetic boundary layer are observed. It is shown that a higher-order LB method accurately captures the characteristics observed in the DSMC method.

Lattice Boltzmann simulation of catalytic reactions

A lattice Boltzmann model is developed to simulate finite-rate catalytic surface chemistry. Diffusive wall boundary conditions are established to account for catalytic reactions in multicomponent mixtures. Implementation of wall boundary conditions with chemical reactions is based on a general second-order accurate interpolation scheme. Results of lattice Boltzmann simulations for a four-component mixture with a global catalytic methane oxidation reaction in a straight channel are in excellent agreement with a finite volume Navier-Stokes solver in terms of both the flow field and species concentrations.

Capillary filling with pseudo-potential binary Lattice-Boltzmann model

We present a systematic study of capillary filling for a binary fluid by using a mesoscopic lattice Boltzmann model for immiscible fluids describing a diffusive interface moving at a given contact angle with respect to the walls. The phenomenological way to impose a given contact angle is analysed. Particular attention is given to the case of complete wetting, that is contact angle equal to zero. Numerical results yield quantitative agreement with the theoretical Washburn's law, provided that the correct ratio of the dynamic viscosities between the two fluids is used. Finally, the presence of precursor films is experienced and it is shown that these films advance in time with a square-root law but with a different prefactor with respect to the bulk interface.

Multiple-relaxation-time lattice Boltzmann scheme for homogeneous mixture flows with external force

A lattice Boltzmann scheme is developed for homogeneous mixture modeling, based on the multiple-relaxation-time (MRT) formulation, which fully recovers the Maxwell-Stefan diffusion model in the continuum limit with (a) external force and (b) tunable Schmidt number. The theoretical basis of the proposed MRT formulation is a recently proposed Bhatnagar-Gross-Krook-type kinetic model for gas mixtures [Andries, J. Stat. Phys. 106, 993 (2002)] and it substantially extends the applicability of a scheme already proposed by the same author, which used only one relaxation parameter. The recovered equations at the macroscopic level are derived by an innovative expansion technique, based on the Grad moment system. Some numerical simulations are reported for the solvent test case with external force, aiming to find the numerical ranges for the transport coefficients that ensure acceptable accuracies. The numerical results reduce the theoretical expectations, which are based on a strong separation among the characteristic scales.