Discrete effects on boundary conditions for the lattice Boltzmann equation in simulating microscale gas flows

Novel Schemes of No-Slip Boundary Conditions for the Discrete Unified Gas Kinetic Scheme Based on the Moment Constraints

The boundary conditions are crucial for numerical methods. This study aims to contribute to this growing area of research by exploring boundary conditions for the discrete unified gas kinetic scheme (DUGKS). The importance and originality of this study are that it assesses and validates the novel schemes of the bounce back (BB), non-equilibrium bounce back (NEBB), and Moment-based boundary conditions for the DUGKS, which translate boundary conditions into constraints on the transformed distribution functions at a half time step based on the moment constraints. A theoretical assessment shows that both present NEBB and Moment-based schemes for the DUGKS can implement a no-slip condition at the wall boundary without slip error. The present schemes are validated by numerical simulations of Couette flow, Poiseuille flow, Lid-driven cavity flow, dipole-wall collision, and Rayleigh-Taylor instability. The present schemes of second-order accuracy are more accurate than the original schemes. Both present NEBB and Moment-based schemes are more accurate than the present BB scheme in most cases and have higher computational efficiency than the present BB scheme in the simulation of Couette flow at high Re. The present Moment-based scheme is more accurate than the present BB, NEBB schemes, and reference schemes in the simulation of Poiseuille flow and dipole-wall collision, compared to the analytical solution and reference data. Good agreement with reference data in the numerical simulation of Rayleigh-Taylor instability shows that they are also of use to the multiphase flow. The present Moment-based scheme is more competitive in boundary conditions for the DUGKS.

Three-Dimensional Simulations of Anisotropic Slip Microflows Using the Discrete Unified Gas Kinetic Scheme

The specific objective of the present work study is to propose an anisotropic slip boundary condition for three-dimensional (3D) simulations with adjustable streamwise and spanwise slip length by the discrete unified gas kinetic scheme (DUGKS). The present boundary condition is proposed based on the assumption of nonlinear velocity profiles near the wall instead of linear velocity profiles in a unidirectional steady flow. Moreover, a 3D corner boundary condition is introduced to the DUGKS to reduce the singularities. Numerical tests validate the effectiveness of the present method, which is more accurate than the bounce-back and specular reflection slip boundary condition in the lattice Boltzmann method. It is of significance to study the lid-driven cavity flow due to its applications and its capability in exhibiting important phenomena. Then, the present work explores, for the first time, the effects of anisotropic slip on the two-sided orthogonal oscillating micro-lid-driven cavity flow by adopting the present method. This work will generate fresh insight into the effects of anisotropic slip on the 3D flow in a two-sided orthogonal oscillating micro-lid-driven cavity. Some findings are obtained: The oscillating velocity of the wall has a weaker influence on the normal velocity component than on the tangential velocity component. In most cases, large slip length has a more significant influence on velocity profiles than small slip length. Compared with pure slip in both top and bottom walls, anisotropic slip on the top wall has a greater influence on flow, increasing the 3D mixing of flow. In short, the influence of slip on the flow field depends not only on slip length but also on the relative direction of the wall motion and the slip velocity. The findings can help in better understanding the anisotropic slip effect on the unsteady microflow and the design of microdevices.

Simulation of Water Flow in a Nanochannel with a Sudden Contraction or Expansion

Water flow in a nanoscale channel is known to be affected by strong water-wall interactions; as a result, the flow considerably deviates from the conventional continuum flow. Nanochannel with a sudden contraction or expansion is the most fundamental morphological nanostructure in many nanoporous systems such as shale matrix, mudrock, membrane, etc. However, the nanoconfinement effects of water flow in nanoporous systems and its effect on macroscopic flow behavior are still evolving research topics. In this work, our recently developed pore-scale lattice Boltzmann method (LBM) considering the nanoscale effects is extended to directly simulate water flow in nanoporous systems. The results show that the flow rate is dramatically decreased in hydrophobic nanopores because of additional flow resistances at the contraction and expansion junctions. This indicates that the bundle of capillary models or the permeability averaging method overestimates the water flow rate in nanoporous media if the contraction/expansion effects between different nanopores are ignored. This work highlights the importance of wettability of the nanochannel in the determination of dynamic water flow behaviors in the contraction/expansion nanosystem. Other important aspects, including velocity distribution, flow patterns, and vortex characteristics as well as pressure variation along the flow direction, are for the first time revealed and quantified. Large differences can be found comparing gas or larger-scale water flows through the same system. A new type of pressure variation curve along the axis of flow direction is found in the hydrophobic nanochannel with a sudden contraction/expansion. This work provides the fundamental understanding of water transport through the nanoscale system with contraction and expansion effects, giving implications to a wide range of industry applications.

Onsager-regularized lattice Boltzmann method: A nonequilibrium thermodynamics-based regularized lattice Boltzmann method

The regularized class of lattice Boltzmann methods (LBMs) leverage the potency of the standard lattice-Bhatnagar-Gross-Krook method by filtering out spurious nonhydrodynamic moments from the moment space; this is achieved through evaluating regularized populations via a multiscale or a Hermite polynomial expansion approach. In this paper, we propose an alternative approach for evaluating the lattice populations. This approach is based on a kinetic theory that is consistent with nonequilibrium thermodynamics and obeys the Onsager-reciprocity principle. The proposed method is verified and validated for a number of canonical problems such as the athermal shock tube, the double periodic shear layer, the lid driven cavity, flow past square cylinder, and Poiseuille flow at nonvanishing Knudsen numbers. Additionally, the proposed method is compared to existing regularized LBM schemes and is shown to yield significant improvement in the stability and accuracy of the simulations.

Influence of numerical resolution on the dynamics of finite-size particles with the lattice Boltzmann method

We investigate and compare the accuracy and efficiency of different numerical approaches to model the dynamics of finite-size particles using the lattice Boltzmann method (LBM). This includes the standard bounce-back (BB) and the equilibrium interpolation (EI) schemes. To accurately compare the different implementations, we first introduce a boundary condition to approximate the flow properties of an unbounded fluid in a finite simulation domain, taking into account the perturbation induced by a moving particle. We show that this boundary treatment is efficient in suppressing detrimental effects on the dynamics of spherical and ellipsoidal particles arising from the finite size of the simulation domain. We then investigate the performances of the BB and EI schemes in modeling the dynamics of a spherical particle settling under Stokes conditions, which can now be reproduced with great accuracy thanks to the treatment of the exterior boundary. We find that the EI scheme outperforms the BB scheme in providing a better accuracy scaling with respect to the resolution of the settling particle, while suppressing finite-size effects due to the particle discretization on the lattice grid. Additionally, in order to further increase the capability of the algorithm in modeling particles of sizes comparable to the lattice spacing, we propose an improvement to the EI scheme, the complete equilibrium interpolation (CEI). This approach allows us to accurately capture the boundaries of the particle also when located between two fluid nodes. We evaluate the CEI performance in solving the dynamics of an under-resolved particle under analogous Stokes conditions and also for the case of a rotating ellipsoid in a shear flow. Finally, we show that EI and CEI are able to recover the correct flow solutions also at small, but finite, Reynolds number. Adopting the CEI scheme it is not only possible to detect particles with zero lattice occupation, but also to increase up to one order of magnitude the accuracy of the dynamics of particles with a size comparable to the lattice spacing with respect to the BB and the EI schemes.

Mesoscopic method to study water flow in nanochannels with different wettability

Molecular dynamics (MD) simulations is currently the most popular and credible tool to model water flow in nanoscale where the conventional continuum equations break down due to the dominance of fluid-surface interactions. However, current MD simulations are computationally challenging for the water flow in complex tube geometries or a network of nanopores, e.g., membrane, shale matrix, and aquaporins. We present a novel mesoscopic lattice Boltzmann method (LBM) for capturing fluctuated density distribution and a nonparabolic velocity profile of water flow through nanochannels. We incorporated molecular interactions between water and the solid inner wall into LBM formulations. Details of the molecular interactions were translated into true and apparent slippage, which were both correlated to the surface wettability, e.g., contact angle. Our proposed LBM was tested against 47 published cases of water flow through infinite-length nanochannels made of different materials and dimensions-flow rates as high as seven orders of magnitude when compared with predictions of the classical no-slip Hagen-Poiseuille (HP) flow. Using the developed LBM model, we also studied water flow through finite-length nanochannels with tube entrance and exit effects. Results were found to be in good agreement with 44 published finite-length cases in the literature. The proposed LBM model is nearly as accurate as MD simulations for a nanochannel, while being computationally efficient enough to allow implications for much larger and more complex geometrical nanostructures.

Analysis and assessment of the no-slip and slip boundary conditions for the discrete unified gas kinetic scheme

The discrete unified gas kinetic scheme (DUGKS) with a force term is a finite volume solver for the Boltzmann equation. Unlike the standard lattice Boltzmann method (LBM), DUGKS can be applied on nonuniform grids. For both the LBM and DUGKS, the boundary conditions need to be processed through the density distribution function. So researchers introduced the boundary conditions from the LBM frame into the DUGKS. However, the accuracy of these boundary conditions in the DUGKS has not been studied thoroughly. Through strict theoretical deduction, we find that the bounce-back (BB) scheme leads to a different dependence of the numerical error term in the DUGKS as compared to the LBM. The error term is influenced by the relaxation time and the body force. And it can be reduced by lowering the kinetic viscosity. Unlike the BB scheme, the nonequilibrium bounce-back scheme has the ability to implement real no-slip boundary condition. Furthermore, two slip boundary conditions incorporated with Navier's slip model are introduced from the LBM framework into the DUGKS. The tangential momentum change-based (TMAC) scheme can be used directly in the DUGKS because it generates no numerical error term in the DUGKS. For the combination of the bounce-back and specular reflection schemes (BSR), the relation between the slip length and the combination parameter should be modified in accordance with the numerical error term. Analysis shows that the TMAC scheme can simulate a wider range of slip length than the BSR scheme. Numerical simulations of the Couette flow and the Poiseuille flow confirm our theoretical analysis.

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime. II. Application to curved boundaries

Gaseous flows inside microfluidic devices often fall in the slip-flow regime. According to this theoretical description, the Navier-Stokes model remains applicable in bulk, while at solid walls a slip velocity boundary model shall be considered. Physically, it is well established that, to properly account for the wall curvature, the wall slip velocity must be determined by the shear stress, rather than the normal component of the velocity derivative alone, as commonly applied to planar surfaces. It follows that the numerical transcription of this type of boundary condition is generally a challenging task for standard computational fluid dynamics (CFD) techniques. This paper aims to show that the realization of the slip velocity condition on arbitrarily shaped boundaries can be accomplished in a natural way with the lattice Boltzmann method (LBM). To substantiate this conclusion, this work undertakes the following three studies. First, we examine the conditions under which the generic reflection-type boundary rules used by LBM become consistent models for the slip velocity boundary condition. This effort makes use of the second-order Chapman-Enskog expansion method, where we address both planar and curved boundaries. The analysis also clarifies the capabilities and limitations behind the considered reflection-type slip schemes. Second, we revisit the family of parabolic accurate LBM slip boundary schemes, originally formulated in [Phys. Rev. E 96, 013311 (2017)2470-004510.1103/PhysRevE.96.013311] on the basis of the multireflection framework, and discuss their characteristics when operating on curved boundaries as well as the limitations of other less accurate LBM slip boundary formulations, such as the linearly accurate slip schemes and the widely popular "kinetic-based" boundary schemes. In addition, we also discuss the numerical stability of the parabolic slip schemes previously developed, providing an heuristic strategy to improve their stable range of operation. Third, we evaluate the performance of the several slip boundary schemes debated in this paper. The numerical tests correspond to two classical 2D benchmark flow problems of slip over non-planar solid surfaces, namely: (i) the velocity profile of the cylindrical Couette flow, and (ii) the permeability of a slow rarefied gas over a periodic array of circular cylindrical obstacles. The obtained numerical results confirm the competitiveness of the LBM when equipped with slip boundary schemes of parabolic accuracy as CFD tool to simulate slippage phenomena over arbitrarily non-planar surfaces. Indeed, although operating on a simple uniform mesh discretization, the LBM yields a similar, or even superior, level of accuracy compared to state-of-the-art FEM simulations conducted on hardworking body-fitted meshes. This conclusion establishes the LBM as a very appealing CFD technique for simulating microfluidic flows in the slip-flow regime, a result that deserves further exploration in future studies.

Lattice Boltzmann Simulation of the Hydrodynamic Entrance Region of Rectangular Microchannels in the Slip Regime

Developing a three-dimensional laminar flow in the entrance region of rectangular microchannels has been investigated in this paper. When the hydrodynamic development length is the same magnitude as the microchannel length, entrance effects have to be taken into account, especially in relatively short ducts. Simultaneously, there are a variety of non-continuum or rarefaction effects, such as velocity slip and temperature jump. The available data in the literature appearing on this issue is quite limited, the available study is the semi-theoretical approximate model to predict pressure drop of developing slip flow in rectangular microchannels with different aspect ratios. In this paper, we apply the lattice Boltzmann equation method (LBE) to investigate the developing slip flow through a rectangular microchannel. The effects of the Reynolds number (1 < Re < 1000), channel aspect ratio (0 < ε < 1), and Knudsen number (0.001 < Kn < 0.1) on the dimensionless hydrodynamic entrance length, and the apparent friction factor, and Reynolds number product, are examined in detail. The numerical solution of LBM can recover excellent agreement with the available data in the literature, which proves its accuracy in capturing fundamental fluid characteristics in the slip-flow regime.

Impact of the kinetic boundary condition on porous media flow in the lattice Boltzmann formulation

To emphasize the importance of the kinetic boundary condition for micro- to nanoscale flow, we present an ad hoc kinetic boundary condition suitable for torturous geological porous media. We found that the kinetic boundary condition is one of the essential features which should be supplemented to the standard lattice Boltzmann scheme in order to obtain accurate continuum observables. The claim is validated using a channel flow setup by showing the agreement of mass flux with analytical value. Further, using a homogeneous porous structure, the importance of the kinetic boundary condition is shown by comparing the permeability correction factor with the analytical value. Finally, the proposed alternate to the kinetic boundary condition is validated by showing its capability to capture the basic feature of the kinetic boundary condition.

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime: Application to plane boundaries

The first nonequilibrium effect experienced by gaseous flows in contact with solid surfaces is the slip-flow regime. While the classical hydrodynamic description holds valid in bulk, at boundaries the fluid-wall interactions must consider slip. In comparison to the standard no-slip Dirichlet condition, the case of slip formulates as a Robin-type condition for the fluid tangential velocity. This makes its numerical modeling a challenging task, particularly in complex geometries. In this work, this issue is handled with the lattice Boltzmann method (LBM), motivated by the similarities between the closure relations of the reflection-type boundary schemes equipping the LBM equation and the slip velocity condition established by slip-flow theory. Based on this analogy, we derive, as central result, the structure of the LBM boundary closure relation that is consistent with the second-order slip velocity condition, applicable to planar walls. Subsequently, three tasks are performed. First, we clarify the limitations of existing slip velocity LBM schemes, based on discrete analogs of kinetic theory fluid-wall interaction models. Second, we present improved slip velocity LBM boundary schemes, constructed directly at discrete level, by extending the multireflection framework to the slip-flow regime. Here, two classes of slip velocity LBM boundary schemes are considered: (i) linear slip schemes, which are local but retain some calibration requirements and/or operation limitations, (ii) parabolic slip schemes, which use a two-point implementation but guarantee the consistent prescription of the intended slip velocity condition, at arbitrary plane wall discretizations, further dispensing any numerical calibration procedure. Third and final, we verify the improvements of our proposed slip velocity LBM boundary schemes against existing ones. The numerical tests evaluate the ability of the slip schemes to exactly accommodate the steady Poiseuille channel flow solution, over distinct wall slippage conditions, namely, no-slip, first-order slip, and second-order slip. The modeling of channel walls is discussed at both lattice-aligned and non-mesh-aligned configurations: the first case illustrates the numerical slip due to the incorrect modeling of slippage coefficients, whereas the second case adds the effect of spurious boundary layers created by the deficient accommodation of bulk solution. Finally, the slip-flow solutions predicted by LBM schemes are further evaluated for the Knudsen's paradox problem. As conclusion, this work establishes the parabolic accuracy of slip velocity schemes as the necessary condition for the consistent LBM modeling of the slip-flow regime.

Discrete effect on the halfway bounce-back boundary condition of multiple-relaxation-time lattice Boltzmann model for convection-diffusion equations

In this paper, we will focus on the multiple-relaxation-time (MRT) lattice Boltzmann model for two-dimensional convection-diffusion equations (CDEs), and analyze the discrete effect on the halfway bounce-back (HBB) boundary condition (or sometimes called bounce-back boundary condition) of the MRT model where three different discrete velocity models are considered. We first present a theoretical analysis on the discrete effect of the HBB boundary condition for the simple problems with a parabolic distribution in the x or y direction, and a numerical slip proportional to the second-order of lattice spacing is observed at the boundary, which means that the MRT model has a second-order convergence rate in space. The theoretical analysis also shows that the numerical slip can be eliminated in the MRT model through tuning the free relaxation parameter corresponding to the second-order moment, while it cannot be removed in the single-relaxation-time model or the Bhatnagar-Gross-Krook model unless the relaxation parameter related to the diffusion coefficient is set to be a special value. We then perform some simulations to confirm our theoretical results, and find that the numerical results are consistent with our theoretical analysis. Finally, we would also like to point out the present analysis can be extended to other boundary conditions of lattice Boltzmann models for CDEs.

Gas flow through rough microchannels in the transition flow regime

A multiple-relaxation-time lattice Boltzmann model of Couette flow is developed to investigate the rarified gas flow through microchannels with roughness characterized by fractal geometry, especially to elucidate the coupled effects of roughness and rarefaction on microscale gas flow in the transition flow regime. The results indicate that the surface roughness effect on gas flow behavior becomes more significant in rarefied gas flow with the increase of Knudsen number. We find the gas flow behavior in the transition flow regime is more sensitive to roughness height than that in the slip flow regime. In particular, the influence of fractal dimension on rarefied gas flow behavior is less significant than roughness height.

Numerical simulation of gas flow and heat transfer in a rough microchannel using the lattice Boltzmann method

In microfluidics, two important factors responsible for the differences between the characteristics of the flow and heat transfer in microchannels and conventional channels are rarefaction and surface roughness which are studied in the present work. An incompressible gas flow in a microchannel is simulated two dimensionally using the lattice Boltzmann method. The flow is in the slip regime and surface roughness is modeled by both regular and Gaussian random distribution of rectangular modules. The effects of relative surface roughness height and Knudsen number on gaseous flow and heat transfer are studied. It was shown that as the relative roughness height increases, the Poiseuille number increases and the Nusselt number has a decreasing or increasing trend, depending on the degree of rarefaction. A comparison between the flow and heat transfer characteristics in regular and random distribution of surface roughness demonstrates that in regular roughness, circular flows are more pronounced; Poiseuille number is higher and Nusselt number is lower than that of its equivalent random roughness.

Boundary condition for lattice Boltzmann modeling of microscale gas flows with curved walls in the slip regime

The lattice Boltzmann method (LBM) has been widely used to simulate microgaseous flows in recent years. However, it is still a challenging task for LBM to model that kind of microscale flow involving complex geometries, owing to the use of uniform Cartesian lattices in space. In this work, a boundary condition for microflows in the slip regime is developed for LBM in which the shape of a solid wall is well considered. The proposed treatment is a combination of the Maxwellian diffuse reflection scheme and the simple bounce-back method. A portion of each part is determined by the relative position between the boundary node and curved walls, which is the key point that distinguishes this method from some previous schemes where the smooth curved surface was assumed to be zigzag lines. The present curved boundary condition is implemented with the multiple-relaxation-times model and verified for several established cases, including the plane microchannel flow (aligned and inclined), microcylindrical Couette flow, and the flow over an inclined microscale airfoil. The numerical results agree well with those predicted by the direct simulation Monte Carlo method.

Lattice Boltzmann simulation of shale gas transport in organic nano-pores

Permeability is a key parameter for investigating the flow ability of sedimentary rocks. The conventional model for calculating permeability is derived from Darcy's law, which is valid only for continuum flow in porous rocks. We discussed the feasibility of simulating methane transport characteristics in the organic nano-pores of shale through the Lattice Boltzmann method (LBM). As a first attempt, the effects of high Knudsen number and the associated slip flow are considered, whereas the effect of adsorption in the capillary tube is left for future work. Simulation results show that at small Knudsen number, LBM results agree well with Poiseuille's law, and flow rate (flow capacity) is proportional to the square of the pore scale. At higher Knudsen numbers, the relaxation time needs to be corrected. In addition, velocity increases as the slip effect causes non negligible velocities on the pore wall, thereby enhancing the flow rate inside the pore, i.e., the permeability. Therefore, the LBM simulation of gas flow characteristics in organic nano-pores provides an effective way of evaluating the permeability of gas-bearing shale.

Filter-matrix lattice Boltzmann model for microchannel gas flows

The lattice Boltzmann method has been shown to be successful for microscale gas flows, and it has attracted significant research interest. In this paper, the recently proposed filter-matrix lattice Boltzmann (FMLB) model is first applied to study the microchannel gas flows, in which a Bosanquet-type effective viscosity is used to capture the flow behaviors in the transition regime. A kinetic boundary condition, the combined bounce-back and specular-reflection scheme with the second-order slip scheme, is also designed for the FMLB model. By analyzing a unidirectional flow, the slip velocity and the discrete effects related to the boundary condition are derived within the FMLB model, and a revised scheme is presented to overcome such effects, which have also been validated through numerical simulations. To gain an accurate simulation in a wide range of Knudsen numbers, covering the slip and the entire transition flow regimes, a set of slip coefficients with an introduced fitting function is adopted in the revised second-order slip boundary condition. The periodic and pressure-driven microchannel flows have been investigated by the present model in this study. The numerical results, including the velocity profile and the mass flow rate, as well as the nonlinear pressure distribution along the channel, agree fairly well with the solutions of the linearized Boltzmann equation, the direct simulation Monte Carlo results, the experimental data, and the previous results of the multiple effective relaxation lattice Boltzmann model. Also, the present results of the velocity profile and the mass flow rate show that the present model with the fitting function can yield improved predictions for the microchannel gas flow with higher Knudsen numbers in the transition flow regime.

Role of rough surface topography on gas slip flow in microchannels

We conduct a lattice Boltzmann simulation of gas slip flow in microchannels incorporating rough surface effects as characterized by fractal geometry with a focus on gas-solid interaction. The gas slip flow in rough microchannels, which is characterized by Poiseuille number and mass flow rate, is evaluated and compared with smooth microchannels. The effects of roughness height, surface fractal dimension, and Knudsen number on slip behavior of gas flow in microchannels are all investigated and discussed. The results indicate that the presence of surface roughness reduces boundary slip for gas flow in microchannels with respect to a smooth surface. The gas flows at the valleys of rough walls are no-slip while velocity slips are observed over the top of rough walls. We find that the gas flow behavior in rough microchannels is insensitive to the surface topography irregularity (unlike the liquid flow in rough microchannels) but is influenced by the statistical height of rough surface and rarefaction effects. In particular, decrease in roughness height or increase in Knudsen number can lead to large wall slip for gas flow in microchannels.

Analytical modeling for heat transfer in sheared flows of nanofluids

We developed a model for the enhancement of the heat flux by spherical and elongated nanoparticles in sheared laminar flows of nanofluids. Besides the heat flux carried by the nanoparticles, the model accounts for the contribution of their rotation to the heat flux inside and outside the particles. The rotation of the nanoparticles has a twofold effect: it induces a fluid advection around the particle and it strongly influences the statistical distribution of particle orientations. These dynamical effects, which were not included in existing thermal models, are responsible for changing the thermal properties of flowing fluids as compared to quiescent fluids. The proposed model is strongly supported by extensive numerical simulations, demonstrating a potential increase of the heat flux far beyond the Maxwell-Garnett limit for the spherical nanoparticles. The road ahead, which should lead toward robust predictive models of heat flux enhancement, is discussed.

Microscale boundary conditions of the lattice Boltzmann equation method for simulating microtube flows

The lattice Boltzmann equation (LBE) method has been shown to be a promising tool for microscale gas flows. However, few works focus on the microtube flows, and there still are some fundamental problems for the LBE to such flows. In this paper, a recently proposed axisymmetric LBE with three kinetic boundary conditions, i.e., the combination of bounceback and specular reflection scheme, the combination of the Maxwell and specular-reflection scheme, and the combination of the Maxwell and bounceback scheme, have been investigated in detail. By analyzing the micro-Hagen-Poiseuille flow, we observed the discrete boundary condition effect and provided a revised boundary scheme to overcome such effect near the slip flow regime. Some numerical tests for the micro-Hagen-Poiseuille have been carried out to validate the analysis, and the numerical results of the revised boundary schemes agree well with the analytic solutions which confirmed our theoretical analysis. In addition, we also applied the revised combination of the Maxwell and bounceback scheme to microtube flow with sudden expansion and contraction, the numerical results of the pressure distribution and normalized slip velocity agree well with the theoretical ones.

General bounce-back scheme for concentration boundary condition in the lattice-Boltzmann method

In this paper, a general bounce-back scheme is proposed to implement concentration or thermal boundary conditions of convection-diffusion equation with the lattice Boltzmann method (LBM). Using this scheme, the general concentration boundary conditions, i.e., b1(∂Cw/∂n) + b2Cw = b3, can be easily implemented at boundaries with complex geometry structure like that in porous media. The numerical results obtained using the present scheme are in excellent agreement with the analytical solutions of flows with both stationary and moving interfaces. Furthermore, to better understand the halfway bounce-back scheme, an analytical study of the concentration jump is presented. The studies of theoretical analysis and numerical experiments demonstrate that the proposed scheme has second-order accuracy.

Higher order slip according to the linearized Boltzmann equation with general boundary conditions

In the present paper, we provide an analytical expression for the first- and second-order velocity slip coefficients by means of a variational technique that applies to the integrodifferential form of the Boltzmann equation based on the true linearized collision operator and the Cercignani-Lampis scattering kernel of the gas-surface interaction. The polynomial form of the Knudsen number obtained for the Poiseuille mass flow rate and the values of the velocity slip coefficients are analysed in the frame of potential applications of the lattice Boltzmann methods in simulations of microscale flows.

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results give assistance in the optimization of microchannel flows for large throughput.

Robust thermal boundary conditions applicable to a wall along which temperature varies in lattice-gas cellular automata

We show that the heat exchange between fluid particles and boundary walls can be achieved by controlling the velocity change rate following the particles' collision with a wall in discrete kinetic theory, such as the lattice-gas cellular automata and the lattice Boltzmann method. We derive a relation between the velocity change rate and temperature so that we can control the velocity change rate according to a given temperature boundary condition. This relation enables us to deal with the thermal boundary whose temperature varies along a wall in contrast to the previous works of the lattice-gas cellular automata. In addition, we present simulation results to compare our method to the existing and give an example in a microchannel with a high temperature gradient boundary condition by the lattice-gas cellular automata.

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.

Slip velocity and Knudsen layer in the lattice Boltzmann method for microscale flows

We present mesoscopic fluid-wall interaction models for lattice Boltzmann (LB) model simulations of microscale flows. The exact solution of the slip velocity for the LB equation with the Bhatnagar-Gross-Krook collision operator is obtained for Poiseuille flow at finite Knudsen numbers. With a consistent definition of the Knudsen number, the slip coefficients of the LB equation with the standard D2Q9 scheme are found to be slightly larger than those of the Boltzmann equation with the same boundary condition, which makes the standard LB method remain quantitatively accurate only for small Knudsen numbers. By modifying the nonequilibrium energy flux or introducing the effective relaxation time, the LB method is analytically shown to reproduce the slip phenomena up to second order in the Knudsen number. For the standard LB method, the Knudsen layer is captured only with modification of the relaxation dynamics such as in the effective relaxation time model.