岩土力学 ›› 2018, Vol. 39 ›› Issue (12): 4527-4536.doi: 10.16285/j.rsm.2017.0827

• 基础理论与实验研究 • 上一篇    下一篇

基于SPH方法的非饱和多孔介质相变耦合研究

饶登宇1,白 冰1,陈佩佩2   

  1. 1. 北京交通大学 土木建筑工程学院,北京 100044;2. 北京建筑大学 理学院,北京 102616
  • 收稿日期:2017-04-27 出版日期:2018-12-11 发布日期:2018-12-31
  • 通讯作者: 白冰,男,1966 年生,博士,教授,博士生导师,主要从事岩土介质热力学特性方面的研究工作。E-mail: bbai@bjtu.edu.cn E-mail:raodengyu@sina.com
  • 作者简介:饶登宇,男,1995年生,博士研究生,主要从事岩土工程数值传热学等方面的研究工作。
  • 基金资助:
    国家自然科学基金资助项目(No.51478034,No.51678043);北京建筑大学科学研究基金资助项目(No.KYJJ2017018)

Simulation of hydro-thermal coupling with phase-change in unsaturated porous media by SPH method

RAO Deng-yu1, BAI Bing1, CHEN Pei-pei2   

  1. 1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China; 2. School of Science, Beijing University of Civil Engineering and Architecture, Beijing 102616, China
  • Received:2017-04-27 Online:2018-12-11 Published:2018-12-31
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (51478034, 51678043) and the Beijing University of Civil Engineering and Architecture Research Foundation (KJJ2017018).

摘要: 在考虑相变的热能平衡方程和非饱和水分迁移质量控制方程的基础上,建立温度场-水分场的耦合模型,并采用一种无网格粒子算法(SPH)进行数值求解。其中,耦合方程中考虑了水流传热以及温度势对水流的直接驱动,在不考虑相变的情况下,该耦合模型可退化为常温下的水-热耦合模型,故可用于模拟冻融循环的相关问题。从求解热能平衡方程中的含冰量出发,实现解耦并对半无限单向冻结条件下介质内非稳态温度场和体积含水率分布场进行模拟,将耦合作用下的温度场与不耦合的解析解进行对比,反映出水分迁移对温度场存在较大影响。最后,求解了路基边坡在季节性周期温度边界下,温度场、水分场分布的演变规律,并评估了边坡阴阳面受热不均对水热两场分布的影响。计算结果基本能反映土冻结相变的实际物理过程,光滑粒子算法可以用于尝试解决冻土领域的其他相关问题。

关键词: 非饱和土, 相变, 水-热耦合, 光滑粒子法, 路基

Abstract: Moisture-temperature coupling control equations are introduced into the heat and mass balance equation for the consideration of phase change. Smoothed Particle Hydrodynamics method is used to solve the equations, where both moisture transfer driven by temperature and heat transfer by water are considered. While the phase transition process is ignored, the model can directly degrade to a generic coupled model of moisture and heat at normal temperature. Therefore it can be used to study freezing and thawing cycle of frozen soils (e. g., seasonal or daily variation in ground temperature). To validate the model, the solution of a semi-infinite unsaturated porous medium subject to a frozen wall is studied. The comparison result shows that the process of moisture migration impacts on heat transfer, with the degree intensifying as time. In addition, the model is used in an unsaturated subgrade with periodical temperature boundary. The temperature and moisture distribution of subgrade are different when considering different solar radiation in sunny and shady slope. The evolution of water content field and thermal field is given at the end. Since the algorithm can basically reflect the real physical process of phase change, it’s credible to study other problems in the field of frozen soil by SPH method.

Key words: unsaturated porous media, phase change, hydro-thermal coupling, SPH method, subgrade

中图分类号: 

  • TU 43
[1] 杨长卫, 童心豪, 王栋, 谭信荣, 郭雪岩, 曹礼聪, . 地震作用下有砟轨道路基动力响应 规律振动台试验[J]. 岩土力学, 2020, 41(7): 2215-2223.
[2] 杨志浩, 岳祖润, 冯怀平, . 非饱和粉土路基内水分迁移规律试验研究[J]. 岩土力学, 2020, 41(7): 2241-2251.
[3] 陈昊, 胡小荣. 非饱和土三剪强度准则及验证[J]. 岩土力学, 2020, 41(7): 2380-2388.
[4] 文伟, 赖远明, 尤哲敏, 李积锋, . 基于Pitzer离子模型的盐渍非饱和土孔隙 相对湿度计算[J]. 岩土力学, 2020, 41(6): 1944-1952.
[5] 陶帅, 董毅, 韦昌富, . 环境湿度可控的土体小应变刚度试验系统[J]. 岩土力学, 2020, 41(6): 2132-2142.
[6] 柳鸿博, 周凤玺, 岳国栋, 郝磊超. 非饱和土中热弹性波的传播特性分析[J]. 岩土力学, 2020, 41(5): 1613-1624.
[7] 孙银磊, 汤连生, 刘洁, . 非饱和土微观结构与粒间吸力的研究进展[J]. 岩土力学, 2020, 41(4): 1095-1122.
[8] 李潇旋, 李涛, 李舰, 张涛. 循环荷载下非饱和结构性黏土的弹塑性双面模型[J]. 岩土力学, 2020, 41(4): 1153-1160.
[9] 李华, 李同录, 江睿君, 范江文. 基于滤纸法的非饱和渗透性曲线测试[J]. 岩土力学, 2020, 41(3): 895-904.
[10] 陈仁朋, 王朋飞, 刘鹏, 程威, 康馨, 杨微, . 路基煤矸石填料土-水特征曲线试验研究[J]. 岩土力学, 2020, 41(2): 372-378.
[11] 程涛, 晏克勤, 胡仁杰, 郑俊杰, 张欢, 陈合龙, 江志杰, 刘强, . 非饱和土拟二维平面应变固结问题的解析计算方法[J]. 岩土力学, 2020, 41(2): 453-460.
[12] 邓子千, 陈嘉帅, 王建伟, 刘小文, . 基于SFG模型的统一屈服面本构模型与试验研究[J]. 岩土力学, 2020, 41(2): 527-534.
[13] 李潇旋, 李涛, 彭丽云, . 控制吸力循环荷载下非饱和黏性土 的弹塑性双面模型[J]. 岩土力学, 2020, 41(2): 552-560.
[14] 程昊, 唐辉明, 吴琼, 雷国平. 一种考虑水力滞回效应的非饱和土弹塑性扩展 剑桥本构模型显式算法有限元实现[J]. 岩土力学, 2020, 41(2): 676-686.
[15] 艾希, 冷伍明, 徐方, 张期树, 翟斌, . 新型预应力路基水平附加应力计算的图表法[J]. 岩土力学, 2020, 41(1): 253-266.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!