Evaluation of a lattice Boltzmann method in a complex nanoflow

Multiscale modeling of gas flow behaviors in nanoporous shale matrix considering multiple transport mechanisms

This study proposes a multiscale model combining molecular simulation and the lattice Boltzmann method (LBM) to explore gas flow behaviors with multiple transport mechanisms in nanoporous media of shale matrix. The gas adsorption characteristics in shale nanopores are first investigated by molecular simulations, which are then integrated and upscaled into the LBM model through a local adsorption density parameter. In order to adapt to high Knudsen number and nanoporous shale matrix, a multiple-relaxation-time pore-scale LBM model with a regularization procedure is developed. The combination of bounce-back and full diffusive boundary condition is adopted to take account of gas slippage and surface diffusion induced by gas adsorption. Molecular simulation results at the atomic scale show that gas adsorption behaviors are greatly affected by the pressure and pore size of the shale organic nanopore. At the pore scale, the gas transport behaviors with multiple transport mechanisms in nanoporous shale matrix are explored by the developed multiscale model. Simulation results indicate that pressure exhibits more significant influences on the transport behaviors of shale gas than temperature does. Compared with porosity, the average pore size of nanoporous shale matrix plays a more significant role in determining the apparent permeability of gas transport. The roles of the gas adsorption layer and surface diffusion in shale gas transport are discussed. It is observed that under low pressure, the gas adsorption layer has a positive influence on gas transport in shale matrix due to the strong surface diffusion effect. The nanoporous structure with the anisotropy characteristic parallel to the flow direction can enhance gas transport in shale matrix. The obtained results may provide underlying and comprehensive understanding of gas flow behaviors considering multiple transport mechanisms in shale matrix. Also, the proposed multiscale model can be considered as a powerful tool to invesigate the multiscale and multiphysical flow behaviors in porous media.

Corner-transport-upwind lattice Boltzmann model for bubble cavitation

Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann model that describes a two-dimensional (2D) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner-transport-upwind (CTU) numerical scheme on large square lattices (up to 6144×6144 nodes). The numerical viscosity and the regularization of the model are discussed for first- and second-order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows us to recover the solution of the 2D Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation, and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient D and the capillary number Ca is found at small Ca but with a different factor than in equilibrium liquids. A nonlinear regime is observed for Ca≳0.2.

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.

Comparative study of the discrete velocity and lattice Boltzmann methods for rarefied gas flows through irregular channels

Rooted from the gas kinetics, the lattice Boltzmann method (LBM) is a powerful tool in modeling hydrodynamics. In the past decade, it has been extended to simulate rarefied gas flows beyond the Navier-Stokes level, either by using the high-order Gauss-Hermite quadrature, or by introducing the relaxation time that is a function of the gas-wall distance. While the former method, with a limited number of discrete velocities (e.g., D2Q36), is accurate up to the early transition flow regime, the latter method (especially the multiple relaxation time (MRT) LBM), with the same discrete velocities as those used in simulating hydrodynamics (i.e., D2Q9), is accurate up to the free-molecular flow regime in the planar Poiseuille flow. This is quite astonishing in the sense that less discrete velocities are more accurate. In this paper, by solving the Bhatnagar-Gross-Krook kinetic equation accurately via the discrete velocity method, we find that the high-order Gauss-Hermite quadrature cannot describe the large variation in the velocity distribution function when the rarefaction effect is strong, but the MRT-LBM can capture the flow velocity well because it is equivalent to solving the Navier-Stokes equations with an effective shear viscosity. Since the MRT-LBM has only been validated in simple channel flows, and for complex geometries it is difficult to find the effective viscosity, it is necessary to assess its performance for the simulation of rarefied gas flows. Our numerical simulations based on the accurate discrete velocity method suggest that the accuracy of the MRT-LBM is reduced significantly in the simulation of rarefied gas flows through the rough surface and porous media. Our simulation results could serve as benchmarking cases for future development of the LBM for modeling and simulation of rarefied gas flows in complex geometries.

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.

Study of Gas Flow Characteristics in Tight Porous Media with a Microscale Lattice Boltzmann Model

To investigate the gas flow characteristics in tight porous media, a microscale lattice Boltzmann (LB) model with the regularization procedure is firstly adopted to simulate gas flow in three-dimensional (3D) digital rocks. A shale digital rock and a sandstone digital rock are reconstructed to study the effects of pressure, temperature and pore size on microscale gas flow. The simulation results show that because of the microscale effect in tight porous media, the apparent permeability is always higher than the intrinsic permeability, and with the decrease of pressure or pore size, or with the increase of temperature, the difference between apparent permeability and intrinsic permeability increases. In addition, the Knudsen numbers under different conditions are calculated and the results show that gas flow characteristics in the digital rocks under different Knudsen numbers are quite different. With the increase of Knudsen number, gas flow in the digital rocks becomes more uniform and the effect of heterogeneity of the porous media on gas flow decreases. Finally, two commonly used apparent permeability calculation models are evaluated by the simulation results and the Klinkenberg model shows better accuracy. In addition, a better proportionality factor in Klinkenberg model is proposed according to the simulation results.

Thermal lattice Boltzmann method for complex microflows

A methodology to simulate thermal fields in complex microflow geometries is proposed. For the flow fields, the regularized multiple-relaxation-time lattice Boltzmann method (LBM) is applied coupled with the diffusive-bounce-back boundary condition for wall boundaries. For the thermal fields, the regularized lattice Bhatnagar-Gross-Krook model is applied. For the thermal wall boundary condition, a newly developed boundary condition, which is a mixture of the diffuse scattering and constant temperature conditions, is applied. The proposed set of schemes is validated by reference data in the Fourier flows and square cylinder flows confined in a microchannel. The obtained results confirm that it is essential to apply the regularization to the thermal LBM for avoiding kinked temperature profiles in complex thermal flows. The proposed wall boundary condition is successful to obtain thermal jumps at the walls with good accuracy.

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.