数值分析

基于格子Boltzmann方法的岩体单裂隙面渗流特性研究

展开
  • 1. 昆明理工大学 建筑工程学院,云南 昆明 650500;2. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071
王志良,男,1982年生,博士,副教授,主要从事岩土工程渗流方面的研究工作

收稿日期: 2015-05-18

  网络出版日期: 2018-06-05

基金资助

国家自然科学基金资助(No. 51408284,No. 51508253,No. U1502231);云南省应用基础研究计划项目(No. 2016FB077)。

Seepage characteristics of a single fracture based on lattice Boltzmann method

Expand
  • 1. Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming, Yunnan 650500, China; 2. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China

Received date: 2015-05-18

  Online published: 2018-06-05

Supported by

This work was supported by the National Natural Science Foundation of China (51408284, 51508253, U1502231) and the Yunnan Applied Basic Research Project (2016FB077).

摘要

为了研究层流状态下粗糙单裂隙面的渗流特性,基于格子Boltzmann方法,建立了压力作用下单裂隙面渗流的数值模型。采用D2Q9模型模拟离散速度分布方向,在宏观尺度上,上、下边界设置为不透水边界(ux =uy =0),左、右边界为压力边界(左侧压力pin大于右侧压力pout);在微观尺度上,左、右边界及光滑裂隙表面采用非平衡态外推格式,粗糙裂隙表面采用标准反弹格式。编制了相应的计算程序,验证了光滑平板裂隙流的立方定律。采用分段随机裂隙长及裂隙宽的方法生成了粗糙裂隙面,并基于格子Boltzmann方法研究了不同粗糙裂隙面方案的渗流特性。研究结果表明:对于粗糙裂隙面,渗流特征极大程度取决于裂隙表面形貌,随着相对粗糙度的增加,裂隙渗流规律明显偏离立方定律。为此,考虑裂隙面相对粗糙度的影响,根据不同裂隙粗糙面方案的数值计算结果,提出了层流状态下粗糙裂隙面渗流的立方定律修正公式,为系统研究复杂粗糙裂隙的水力特性奠定了一定的基础。

本文引用格式

王志良,申林方,李邵军,徐则民, . 基于格子Boltzmann方法的岩体单裂隙面渗流特性研究[J]. 岩土力学, 2017 , 38(4) : 1203 -1210 . DOI: 10.16285/j.rsm.2017.04.035

Abstract

To investigate seepage characteristics of a single rough fracture, based on lattice Boltzmann method, a numerical model is established to analyze the water flow in laminar under pressure. The D2Q9 model is applied to simulate the discrete velocity direction. On the macroscopic-scale, the upper and lower boundaries (ux =uy =0) are assumed to be impermeable, and the left and right boundaries are controlled by pressure (i.e., the pressure on the left side is larger than that on the right side). On the microscopic-scale, the non-equilibrium extrapolation scheme is set on the pressure boundary and smooth fracture surface boundary, and the standard bounce-back scheme is set on the rough fracture surface boundary. The corresponding program is compiled to verify the classical cubic law of smooth plate fracture flow. The rough fracture surface is generated by setting piecewise random length and random width, and seepage characteristics of different rough fracture surfaces are discussed in detail. The results show that fracture surface roughness greatly affects seepage characteristics. With the increase of relative roughness, the deviation from cubic law becomes more obvious. Therefore, by considering the effect of relative fracture surface roughness, the modified cubic law is proposed according to numerical results of different schemes of rough fracture surface. This study lays the foundation for further research on the complex hydraulic characteristics of rough fracture.
文章导航

/