›› 2015, Vol. 36 ›› Issue (S2): 689-694.doi: 10.16285/j.rsm.2015.S2.098

• 数值分析 • 上一篇    下一篇

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

申林方1,王志良1,李邵军2   

  1. 1. 昆明理工大学 建筑工程学院,云南 昆明 650500; 2. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071
  • 收稿日期:2014-11-18 出版日期:2015-08-31 发布日期:2018-06-14
  • 通讯作者: 王志良,男,1982年生,博士,副教授,主要从事地下建筑渗流方面的研究工作。E-mail: wangzhiliangtj@126.com E-mail:shenlinfang@kmust.edu.cn
  • 作者简介:申林方,女,1982年生,博士后,主要从事岩土工程渗流方面的研究工作。
  • 基金资助:
    国家自然科学基金(No.51408284,No.51179187);昆明理工大学人才培养基金(No. KKSY201306023)

Numerical simulation for mesoscopic seepage field of soil based on lattice Boltzmann method at REV scale

SHEN Lin-fang1, WANG Zhi-liang1, LI Shao-jun2   

  1. 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:2014-11-18 Online:2015-08-31 Published:2018-06-14

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

关键词: 格子博尔兹曼(Boltzmann)方法, 表征体元(REV)尺度, 多孔介质, 渗流场, 数值模拟

Abstract: In order to study seepage mechanism of soil, there are some basic assumptions, the soil is saturated, and flow is laminar in the process of seepage. Based on lattice Boltzmann method at the representative elementary volume(REV) scale, the mesoscopic seepage numerical model of soil under pressure is established, considering the effects of macroscopic statistical parameters, i.e. porosity, permeability and effective viscosity coefficient. The D2Q9 model is applied for the discrete velocity direction. In the macroscopic boundary condition, it is impermeable in the left and right sides; and it is controlled by pressure in the upper and lower boundaries, which is decided by setting different density. And the non- equilibrium extrapolation scheme is used in the microscopic boundary. The porous media is set to be fluid in the study region by letting 1.0 and , the corresponding program is compiled to verify the Poiseuille flow. Combined with an example, the influence of porosity, permeability and seepage pressure on the seepage velocity driven by pressure is discussed for soil, and the results of this paper are in good agreement with that calculated by Darcy’s law. Therefore, lattice Boltzmann method at the REV scale could effectively simulate the seepage characteristics of soil, and it provides a new research approach for further study of seepage mechanism.

Key words: lattice Boltzmann method, representative elementary volume(REV) scale, porous media, seepage field, numerical simulation

中图分类号: 

  • O 357.3
[1] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[2] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[3] 周凤玺, 柳鸿博, 蔡袁强, . 饱和多孔热弹性介质中Rayleigh波 传播特性分析[J]. 岩土力学, 2020, 41(1): 315-324.
[4] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[5] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[6] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[7] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[8] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[9] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
[10] 吴凤元, 樊赟赟, 陈剑平, 李军, . 基于不同侵蚀模型的高速崩滑碎屑 流动力过程模拟分析[J]. 岩土力学, 2019, 40(8): 3236-3246.
[11] 孙峰, 薛世峰, 逄铭玉, 唐梅荣, 张翔, 李川, . 基于连续损伤的水平井射孔-近井筒三维破裂模拟[J]. 岩土力学, 2019, 40(8): 3255-3261.
[12] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[13] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[14] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[15] 张 聪, 梁经纬, 阳军生, 曹 磊, 谢亦朋, 张贵金, . 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 黄建华,宋二祥. 大型锚碇基础围护工程冻结帷幕力学性态研究[J]. , 2009, 30(11): 3372 -3378 .
[2] 王观石,李长洪,陈保君,李世海. 应力波在非线性结构面介质中的传播规律[J]. , 2009, 30(12): 3747 -3752 .
[3] 王朝阳,许 强,倪万魁. 原状黄土CT试验中应力-应变关系的研究[J]. , 2010, 31(2): 387 -391 .
[4] 邓 琴,郭明伟,李春光,葛修润. 基于边界元法的边坡矢量和稳定分析[J]. , 2010, 31(6): 1971 -1976 .
[5] 万少石,年廷凯,蒋景彩,栾茂田. 边坡稳定强度折减有限元分析中的若干问题讨论[J]. , 2010, 31(7): 2283 -2288 .
[6] 闫 铁,李 玮,毕雪亮. 基于分形方法的多孔介质有效应力模型研究[J]. , 2010, 31(8): 2625 -2629 .
[7] 刘 嘉,王 栋. 正常固结黏土中平板锚基础的吸力和抗拉力[J]. , 2009, 30(3): 735 -740 .
[8] 赵尚毅,郑颖人,李安洪,邱文平,唐晓松,徐 俊. 多排埋入式抗滑桩在武隆县政府滑坡中的应用[J]. , 2009, 30(S1): 160 -164 .
[9] 刘振平,贺怀建,朱发华. 基于钻孔数据的三维可视化快速建模技术的研究[J]. , 2009, 30(S1): 260 -266 .
[10] 魏厚振,颜荣涛,韦昌富,吴二林,陈 盼,田慧会. 含天然气水合物沉积物相平衡问题研究综述[J]. , 2011, 32(8): 2287 -2294 .