基于有限体积法的冻土水热耦合程序开发及验证

doi:10.16285/j.rsm.2019.1263

岩土力学 ›› 2020, Vol. 41 ›› Issue (5): 1781-1789.doi: 10.16285/j.rsm.2019.1263

基于有限体积法的冻土水热耦合程序开发及验证

胡田飞^{1, 2}，王天亮^{1, 2}，常键³，刘建勇^{1, 2}，卢玉婷^{1, 2}

1. 石家庄铁道大学省部共建交通工程结构力学行为与系统安全国家重点实验室，河北石家庄 050043； 2. 石家庄铁道大学土木工程学院，河北石家庄 050043；3. 北京交通大学土木建筑工程学院，北京 100044

收稿日期:2019-07-20 修回日期:2019-09-11 出版日期:2020-05-11 发布日期:2020-07-08
通讯作者: 王天亮，男，1981年生，博士，副教授，硕士生导师，主要从事路基工程方面的教学和科研工作。E-mail: wangtl@stdu.edu.cn E-mail:hutianfei@stdu.edu.cn
作者简介:胡田飞，男，1988年生，博士，讲师，硕士生导师，主要从事冻土力学与寒区路基工程方面的研究工作。
基金资助:
国家自然科学基金资助项目（No. 41731281），河北省高等学校科学技术研究项目（No. QN2020180）。

Code development and verification for coupled process of water migration and heat transfer of frozen soil based on finite volume method

HU Tian-fei^{1, 2}, WANG Tian-liang^{1, 2}, CHANG Jian³, LIU Jian-yong^{1, 2}, LU Yu-ting^{1, 2}

1. State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures, Shijiazhuang Tiedao University, Shijiazhuang, Hebei 050043, China; 2. School of Civil Engineering, Shijiazhuang Tiedao University, Shijiazhuang, Hebei 050043, China; 3. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China

Received:2019-07-20 Revised:2019-09-11 Online:2020-05-11 Published:2020-07-08
Supported by:
This work was supported by the National Natural Science Foundation of China (41731281) and the Science and Technology Research Project of Hebei Education Department, China (QN2020180).

摘要/Abstract

摘要： 土体冻结和融化时的水分迁移、相变与传热是一个相互影响的耦合过程。采用基于有限体积法的开源软件OpenFOAM，编制描述土体冻融过程的水热耦合计算程序。首先，基于土体水分和热量迁移基本方程、水分相变与温度的平衡方程，同时考虑相变对水分特征参数和热特性参数的影响以及相变潜热对传热过程的影响，建立冻土水热耦合数学模型。然后，采用基于多面体网格的有限体积方法对水热耦合控制方程进行空间离散，采用全隐式向后差分方法对方程进行时间离散，由此编制冻土水热耦合计算程序。该程序具有良好的几何适应性、质量和能量守恒性，具备面向复杂问题的并行计算功能。最后，采用该程序对两组不同温度边界条件的室内土体冻结试验进行数值模拟，并与试验结果进行对比，结果表明该程序可以较为准确地模拟土体冻结过程中温度场和水分场的演化特征。

关键词: 冻土, 水热耦合, OpenFOAM, 有限体积法, 温度, 未冻水含量, 含冰量

Abstract: During soil freezing and thawing, water migration, phase transition and heat transfer are a multi-field coupled process in which these factors affect each other. In this study, a mathematical model for the coupling process of water migration and heat transfer was constructed by solving a simultaneous equation group, which contains the basic equations of water migration and heat transfer, and the equilibrium equation of water phase change and temperature. In the mathematical model, the effects of water phase change on hydraulic and thermal parameters were considered, so was the effect of latent heat in water phase change. Also, a code for the coupled heat transfer and water migration simulation was developed based on the finite volume method for the spatial discretization and fully implicit backward difference scheme for the time discretization. The code supporting unstructured mesh can ensure the conservation of mass and energy. It also has a parallel function, which can be used for the calculation of complex problems. Furthermore, an indoor soil freezing tests with two different temperature boundary conditions were conducted, and the corresponding numerical calculations were subsequently conducted. The experimental results generally matched the calculated ones, so the developed numerical simulation program can well simulate the evolution characteristics of the temperature field and water field during the soil freezing process. Hence, the established OpenFOAM program is an effective method to simulate the coupling process of water migration and heat transfer for the frozen soil.

Key words: freezing soil, coupled process of water migration and heat transfer, OpenFOAM, finite volume method, temperature, unfrozen moisture content, ice content

中图分类号:

TU 445

胡田飞, 王天亮, 常键, 刘建勇, 卢玉婷, . 基于有限体积法的冻土水热耦合程序开发及验证[J]. 岩土力学, 2020, 41(5): 1781-1789.

HU Tian-fei, WANG Tian-liang, CHANG Jian, LIU Jian-yong, LU Yu-ting, . Code development and verification for coupled process of water migration and heat transfer of frozen soil based on finite volume method[J]. Rock and Soil Mechanics, 2020, 41(5): 1781-1789.

参考文献

相关文章 15

[1]	任连伟, 任军洋, 孔纲强, 刘汉龙, . 冷热循环下PHC能量桩热力响应和承载性能现场试验[J]. 岩土力学, 2021, 42(2): 529-536.
[2]	季伟伟, 孔纲强, 刘汉龙, 杨庆, . 软塑黄土地区隧道仰拱热力响应特性现场试验[J]. 岩土力学, 2021, 42(2): 558-564.
[3]	张泽, 马巍, ROMAN Lidia, MELNIKOV Andrey, 杨希, 李宏璧. 基于冻融次数−物理时间比拟理论的冻土长期强度预测方法[J]. 岩土力学, 2021, 42(1): 86-92.
[4]	杨高升, 白冰, 姚晓亮, 陈佩佩, . 非饱和冻土水汽迁移与相变过程的光滑粒子法模拟[J]. 岩土力学, 2021, 42(1): 291-300.
[5]	郤保平, 吴阳春, 王帅, 熊贵明, 赵阳升, . 热冲击作用下花岗岩力学特性及其随冷却温度演变规律试验研究[J]. 岩土力学, 2020, 41(S1): 83-94.
[6]	张科, 李娜, 陈宇龙, 刘文连, . 裂隙砂岩变形破裂过程中应变场及红外辐射温度场演化特征研究[J]. 岩土力学, 2020, 41(S1): 95-105.
[7]	徐衍, 周晓敏, 和晓楠, 吴涛, 张建岭, 李森. 矿山竖井井壁与围岩热−固耦合作用分析[J]. 岩土力学, 2020, 41(S1): 217-226.
[8]	周祥运, 孙德安, 罗汀. 核废料处置库近场温度半解析研究[J]. 岩土力学, 2020, 41(S1): 246-254.
[9]	张继文, 穆青翼, 廖红建, 刘芬良, . 考虑土体孔隙比和比表面积影响的未冻结体积含水率曲线模型[J]. 岩土力学, 2020, 41(9): 2913-2921.
[10]	骆赵刚, 汪时机, 杨振北, . 膨胀土湿干胀缩裂隙演化及其定量分析[J]. 岩土力学, 2020, 41(7): 2313-2323.
[11]	张明礼, 温智, 董建华, 王得楷, 岳国栋, 王斌, 高樯. 考虑降雨作用的多年冻土区不同地表土质活动层水热过程差异分析[J]. 岩土力学, 2020, 41(5): 1549-1559.
[12]	孟祥传, 周家作, 韦昌富, 张坤, 沈正艳, 杨周洁, . 盐分对土的冻结温度及未冻水含量的影响研究[J]. 岩土力学, 2020, 41(3): 952-960.
[13]	杨高升, 白冰, 姚晓亮, . 高含冰量冻土路基融化固结规律研究[J]. 岩土力学, 2020, 41(3): 1010-1018.
[14]	陈卫忠, 李翻翻, 雷江, 于洪丹, 马永尚, . 热−水−力耦合条件下黏土岩蠕变特性研究[J]. 岩土力学, 2020, 41(2): 379-388.
[15]	周智超, 王淼, 孟上九, 孙义强, . 低温影响下FBG永久变形计算方法研究[J]. 岩土力学, 2020, 41(12): 4005-4014.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

基于有限体积法的冻土水热耦合程序开发及验证

Code development and verification for coupled process of water migration and heat transfer of frozen soil based on finite volume method

PDF (Chinese)

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 10