格子玻尔兹曼方法及其在大气湍流研究中的应用

文章的目的是对格子玻尔兹曼方法进行系统的介绍，格子玻尔兹曼方法(Lattice Boltzmann Method)的出现直接来源于20世纪60年代的元胞自动机(Cellular Automata)思想，而这一方法用于解决流动现象时，又可以追溯到19世纪的分子运动论，求解的是Boltzmann提出的玻尔兹曼输运方程，因此将这一方法称为格子玻尔兹曼方法，之前也被称为格子气自动机(Lattice Gas Automaton)。该方法多用于研究复杂现象，如材料晶体凝聚时的生长过程、城市土地利用的演化等方面。在20世纪70年代由Hardy、Pomeau和Pazzis建立了第一个用于研究流体运动的格子气自动机，此后，这一方法被广泛用来模拟各种流动问题，诸如二相流、孔隙介质中的渗流等，并根据这一方法开发了相应的商业软件PowerFlow。同时，格子玻尔兹曼方法由于其在微观水平描述运动的特点，成为研究湍流的一个很好的数值计算工具，特别是用其进行直接数值模拟(DNS)计算，成为继传统的差分法、有限体积法和谱方法之后的又一有力的手段。而作为大气运动的一个主要现象的大气湍流，比普通湍流更加复杂，在这里着重介绍了大气湍流的特点和应用格子玻尔兹曼方法模拟湍流的发展过程。     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体（分辨率3.7 μm）为研究对象，在空间节点数多达64 000 000的情况下，基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数，并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现，结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD（即10^-13m²）量级，远低于规则堆叠的圆球结构（>20 000 mD）；土壤的迂曲度为1.40~1.60，明显高于规则堆叠的圆球结构。研究结果表明，渗透率大的样品具有较小的迂曲度，这与结构的空间异质性有较强的关系；土壤的渗透率和迂曲度呈现各向异性；在水力梯度一定的前提下，渗透率较大的样品，纵向弥散系数也较大；同时，结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟，可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

基岩裂隙水运移模型评价

结合近几年的研究成果，按照建立模型时的能量方程不同，将基岩裂隙水运移模型划分为以Darcy定律为基础的模型和以Chezy公式为基础的模型；按照介质概化条件的差异，进一步划分为孔隙介质、裂隙介质和双重介质数学模型，并评价了各种模型的原理、适用条件及其优缺点。     

可压缩流体饱和孔隙介质中孔隙压力波传播数值分析

首次将高精度时空守恒元/解元方法推广到可压缩流体饱和孔隙介质中孔隙压力波传播的数值计算中。将孔隙度梯度从源（汇）项中分离,直接引入流通量,改进了理论模型。通过对孔隙介质激波问题的数值模拟,验证了方法的精度和有效性。在此基础上,提出了孔隙介质中二维黎曼问题,并揭示了孔隙压力波存在接触间断、激波、膨胀波、压缩波等复杂的结构特征。该成果对二氧化碳地质封存、二氧化碳提高石油采收率、页岩气压裂开采以及地震破裂过程的研究具有重要的理论与应用意义。     

基于孔隙网络模型的非水溶相液体运移实验研究进展

进行多孔介质中非水溶相液体（Non Aqueous Phase Liquids，NAPLs）运移的微观机理研究，微观孔隙网络模型实验是目前应用比较广泛且行之有效的方法。通过网络模型实验，获得对NAPLs在多孔介质中运移更深入的认识。从多孔介质孔隙结构测量、孔隙网络模型制作、NAPLs运移网络模型实验和数值模拟4个方面评述了该方向的研究进展，结果显示测量孔隙结构方法、图像刻蚀技术、可视化测量实验数据方法等有力地促进了本实验研究的发展。分析了孔隙网络模型实验存在的问题以及未来的发展趋势，对开展孔隙网络模型实验研究有一定的启发作用。     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

裂缝诱导双相HTI介质模型及其弹性波传播方程

将Biot双相介质理论与Gurevich裂缝各向异性理论相结合,建立了能够同时考虑实际裂缝性储层孔隙性和各向异性的裂缝诱导双相HTI介质模型。从本构方程、动力学方程和动力学达西定律出发,推导出了裂缝诱导双相HTI介质中弹性波传播的一阶速度-应力方程,并针对方程的刚性问题,给出了利用显式二阶时间积分法数值求解该方程时所需要满足的稳定性条件。该方程能够定量地给出双相HTI介质的波场特征与裂缝参数、背景孔隙介质参数之间的关系,描述弹性波在这种介质中的传播机理。     

非饱和土中热弹性波的传播特性分析

基于非饱和多孔介质的研究成果,考虑热效应和孔隙流体迂曲度的影响,研究了非饱和土中热弹性波的传播特性。利用非饱和土中耦合热的固-液-气三相介质的质量平衡方程、渗流连续方程、动量平衡方程和广义非Fourier热传导定律,建立了问题的热弹性波动方程。通过引入势函数,经过理论推导给出了非饱和土中热弹性波的弥散特征方程。结合数值算例,分析了几类热弹性波的波速和衰减系数随迂曲度、热膨胀系数和介质温度等热物理参数的变化规律。结果表明：孔隙水迂曲度的增大将引起P1波、P3波和S波的波速增大,而孔隙气体迂曲度的增大仅使得P2波的波速增大;热膨胀系数的增大将造成P1波波速的增大和热（T）波波速的减小;介质温度的升高将引起各类热弹性波波速的增大;频率、热膨胀系数和介质温度的变化对各类热弹性波的衰减系数均有较大影响,不可忽视。     

考虑固结效应的溶质传输研究

传统的孔隙介质水动力学采用对流-扩散方程,研究溶质在流体中的迁移。在这个过程中,孔隙介质被认为是不变形的,因而是一个稳态问题。针对二维情况下孔隙介质变形对溶质传输的影响,给出了考虑孔隙介质固结效应的溶质传输方程,并且探讨了该类问题的求解方法。     

非饱和土中热弹性波的传播特性分析

基于非饱和多孔介质的研究成果,考虑热效应和孔隙流体迂曲度的影响,研究了非饱和土中热弹性波的传播特性。利用非饱和土中耦合热的固-液-气三相介质的质量平衡方程、渗流连续方程、动量平衡方程和广义非Fourier热传导定律,建立了问题的热弹性波动方程。通过引入势函数,经过理论推导给出了非饱和土中热弹性波的弥散特征方程。结合数值算例分析了几类热弹性波的波速和衰减系数随迂曲度、热膨胀系数和介质温度等热物理参数的变化规律。结果表明：孔隙水迂曲度的增大将引起P_(1)波、P_(3)波和S波的波速增大,而孔隙气体迂曲度的增大仅使得P_(2)波的波速增大;热膨胀系数的增大将造成P_(1)波波速的增大和热（T）波波速的减小;介质温度的升高将引起各类热弹性波波速的增大;频率、热膨胀系数和介质温度的变化对各类热弹性波的衰减系数均有较大影响,不可忽视。     

A double medium model for diffusion in fluid-bearing rock

The concept of a double porosity medium to model fluid flow in fractured rock has been applied to model diffusion in rock containing a small amount of a continuous fluid phase that surrounds small volume elements of the solid matrix. The model quantifies the relative role of diffusion in the fluid and solid phases of the rock. The fluid is the fast diffusion path, but the solid contains the volumetrically significant amount of the diffusing species. The double medium model consists of two coupled differential equations. One equation is the diffusion equation for the fluid concentration; it contains a source term for change in the average concentration of the diffusing species in the solid matrix. The second equation represents the assumption that the change in average concentration in a solid element is proportional to the difference between the average concentration in the solid and the concentration in the fluid times the solid-fluid partition coefficient. The double medium model is shown to apply to laboratory data on iron diffusion in fluid-bearing dunite and to measured oxygen isotope ratios at marble-metagranite contacts. In both examples, concentration profiles are calculated for diffusion taking place at constant temperature, where a boundary value changes suddenly and is subsequently held constant. Knowledge of solid diffusivities can set a lower bound to the length of time over which diffusion occurs, but only the product of effective fluid diffusivity and time is constrained for times longer than the characteristic solid diffusion time. The double medium results approach a local, grain-scale equilibrium model for times that are large relative to the time constant for solid diffusion.     

Hydro-micromechanical modeling of wave propagation in saturated granular crystals

Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.     

格子Boltzmann方法地震波场模拟

格子Ｂｏｌｔｚｍａｎｎ方法是细胞自动机在某些学科中的具体化和应用。它根据微观运动过程的某些基本特征建立简化的、时间和空间完全离散的动力学模型，这种模型的平行行为符合宏观的微分方程。     

基于格子博尔兹曼方法表征体元尺度土体细观渗流场的数值模拟

为了研究土体的细观渗流特性,假设土体是完全饱和且在渗流过程中水分的流动始终处于层流状态。考虑宏观统计参数（孔隙率、渗透率及有效黏滞系数等）的影响,基于表征体元（REV）尺度的格子博尔兹曼（Boltzmann）方法,建立了压力作用下土体细观渗流的数值模型。采用D2Q9模型考虑水分流动的离散速度分布,宏观边界条件为左右侧面为不透水边界 ,上下边界设置不同的密度来控制压力边界,在微观边界条件上采用非平衡态外推格式。编制相应的计算程序,将计算区域内的多孔介质材料设置成流体（孔隙率 1.0,渗透率 ）,验证了经典的Poiseuille流。此外,结合算例分别讨论了土体在压力作用下孔隙率、渗透率及渗透压力等影响因素与渗流速度的相互关系,研究表明该数值方法与Darcy定律得到的计算结果较为吻合。因此,基于REV尺度的格子Boltzmann方法可以有效地模拟土体的渗流机制,为进一步研究土体渗流特性提供了一种新的研究手段。     

求解偏微分方程的格点法

格点法是在计算流体力学中首先发展起来的数值模拟新方法，共根本思想是对问题重新建模，建立直接模拟流体运动的离散格点模型。本文以一阶拟线性双典型方程及Ｋｄｖ方程为例，推广这种求解一般偏微分方程的。它是运用多尺度分析方法，构造出格点模型的演化方程的局部平衡分布函数。数值试验表明，该方法程序实现简单，求解速度快，数值结果令人满意。     

Free-surface Simulations of Newtonian and Non-Newtonian Fluids with the Lattice Boltzmann Method

This paper describes the application of a three-dimensional lattice Boltzmann method (LBM) to Newtonian and non-Newtonian (Bingham fluid in this work) flows with free surfaces. A mass tracking algorithm was incorporated to capture the free surface, whereas Papanastasiou’s modified model was used for Bingham fluids. The lattice Boltzmann method was first validated using two benchmarks: Newtonian flow through a square cross-section tube and Bingham flow through a circular cross-section tube. Afterward, the dam-break problem for the Newtonian fluid and the slump test for Bingham fluid were simulated to validate the free-surface-capturing algorithm. The numerical results were in good agreement with analytical results, as well as other simulations, thereby proving the validity and correctness of the current method. The proposed method is a promising substitute for time-consuming and costly physical experiments to solve problems encountered in geotechnical and geological engineering, such as the surge and debris flow induced by a landslide or earthquake.     

A coupled 3‐dimensional bonded discrete element and lattice Boltzmann method for fluid‐solid coupling in cohesive geomaterials

This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams.     

Limitations of Lattice Boltzmann Modeling of Micro‐Flows in Complex Nanopores

The multiscale transport mechanism of methane in unconventional reservoirs is dominated by slip and transition flows resulting from the ultra-low permeability of micro/nano-scale pores, which requires consideration of the microscale and rarefaction effects. Traditional continuum-based computational fluid dynamics (CFD) becomes problematic when modeling micro-gaseous flow in these multiscale pore networks because of its disadvantages in the treatment of cases with a complicated boundary. As an alternative, the lattice Boltzmann method (LBM), a special discrete form of the Boltzmann equation, has been widely applied to model the multi-scale and multi-mechanism flows in unconventional reservoirs, considering its mesoscopic nature and advantages in simulating gas flows in complex porous media. Consequently, numerous LBM models and slip boundary schemes have been proposed and reported in the literature. This study investigates the predominately reported LBM models and kinetic boundary schemes. The results of these LBM models systematically compare to existing experimental results, analytical solutions of Navier-Stokes, solutions of the Boltzmann equation, direct simulation of Monte Carlo (DSMC) and information-preservation DSMC (IP_DSMC) results, as well as the numerical results of the linearized Boltzmann equation by the discrete velocity method (DVM). The results point out the challenges and limitations of existing multiple-relaxation-times LBM models in predicting micro-gaseous flow in unconventional reservoirs.     

On the implicit immersed boundary method in coupled discrete element and lattice Boltzmann method

The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity.     

The intrinsic permeability of microcracks in porous solids: Analytical models and Lattice Boltzmann simulations

The paper presents a synthesis of analytical modeling and computational simulations of the intrinsic permeability of microcracks, embedded in porous materials taking into account the interaction of the fluid flow in the microcrack with the surrounding porous material. In the first part of the paper, using the DARCY , STOKES , BRINKMAN , and the BEAVERS–JOSEPH approximations, we derive the intrinsic permeability of a plain non‐rough microcrack in terms of the microcrack geometry and the permeability of the porous material surrounding the microcrack. In the second part of the paper, the intrinsic permeability of a microcrack is determined by means of computational simulations using the framework of the lattice Boltzmann method with partial bounceback conditions. The comparison of predictions from the analytical model and the numerical simulations show an excellent agreement. Copyright © 2017 John Wiley & Sons, Ltd.