基于COMSOL Multiphysics的不同水头压力下水热耦合分析
Hydrothermal Coupling Analysis under Different Head Pressures Based on COMSOL Multiphysics
摘要: 冻土的水热耦合问题一直是冻胀融沉的主要原因之一,涉及到土壤中水分的迁移、热量的传递以及相变过程。在多年冻土和季节冻土区,由于环境因素的影响,冻土路基容易发生不均匀沉降、冻胀融沉等灾害,这些灾害都与水分迁移、相变以及温度变化息息相关。为了解决冻土水热耦合问题,本文通过COMSOL Multiphysic软件建立水热耦合的二维模型,研究不同水头压力作用下同一时间段温度和压力的影响、对流出边界总热通量的影响、对总液态水体积的影响,以及对最低温度的影响进行分析,模拟了3种不同水头压力的工况。结果表明:同一时间段的温度和压力随着水头压力的增加而变化;流出边界离开系统的总热通量也随着水头压力的增大而增大;不同工况下水头梯度越大,冰块融化速度越快,达到总液态水体积最大值的时间越短。
Abstract: The hydrothermal coupling problem of frozen soil has always been one of the main reasons for frost heave and thaw settlement, involving the migration of water in the soil, heat transfer, and phase change process. In permafrost and seasonal permafrost areas, due to environmental factors, permafrost roadbeds are prone to disasters such as uneven settlement, frost heave and thaw settlement, which are closely related to water migration, phase change, and temperature changes. In order to solve the problem of hydrothermal coupling in permafrost, this paper establishes a two-dimensional model of hydrothermal coupling using COMSOL Multiphysic software. The influence of temperature and pressure at the same time period under different head pressures, the influence on the total heat flux at the outflow boundary, the influence on the total liquid water volume, and the influence on the minimum temperature are studied. Three different head pressure working conditions are simulated. The results show that the temperature and pressure during the same time period change with the increase of head pressure; The total heat flux leaving the system at the outflow boundary also increases with the increase of head pressure; The larger the head gradient under different working conditions, the faster the melting rate of ice, and the shorter the time to reach the maximum total liquid water volume.
文章引用:刘海圳, 苗志琪. 基于COMSOL Multiphysics的不同水头压力下水热耦合分析[J]. 土木工程, 2024, 13(8): 1319-1331. https://doi.org/10.12677/hjce.2024.138144

参考文献

[1] 路建国, 张明义, 张熙胤, 等. 冻土水热力耦合研究现状及进展[J]. 冰川冻土, 2017, 39(1): 102-111.
[2] 张明礼, 郭宗云, 韩晓斌, 等. 基于COMSOL Multiphysics数学模块的冻土水热耦合分析[J]. 科学技术与工程, 2018, 18(33): 7-12.
[3] 尚松浩, 雷志栋, 杨诗秀, 等. 冻结条件下土壤水热耦合迁移数值模拟的改进[J]. 清华大学学报(自然科学版), 1997, 37(8): 32-34.
[4] 白青波, 李旭, 田亚护, 等. 冻土水热耦合方程及数值模拟研究[J]. 岩土工程学报, 2015, 37(S2): 131-136.
[5] 刘思源. 多年冻土区不同结构路基温度、水分及变形特性研究[D]: [硕士学位论文]. 合肥: 安徽建筑大学, 2023.
[6] 韩昀希. 多场耦合作用下富水砂层人工冻结状态研究[D]: [硕士学位论文]. 石家庄: 石家庄铁道大学, 2022.
[7] 王蕴嘉, 郭惠芹, 叶阳升, 等. 多年冻土区路基填筑短期内冻土温度场数值分析[J]. 中国铁道科学, 2021, 42(4): 9-18.
[8] 曾照军, 唐朝生, 程青, 等. 非饱和土水热耦合模型中水分相变/迁移因素影响分析[J]. 岩土工程学报, 2022, 42(z1): 40-45.
[9] 申林方, 王志良, 李腾风, 等. 非饱和土体水平冻结过程中的水热耦合作用机制研究[J]. 铁道学报, 45(10): 135-143.
[10] 吴东军. 人工冻土水热迁移试验及水热耦合数值模拟研究[D]: [硕士学位论文]. 淮南: 安徽理工大学, 2019.
[11] 王青林, 陈磊, 明锋, 等. 水热耦合模拟中的冻土渗透系数模型选择[J]. 冰川冻土, 2023, 45(5): 1536-1543.
[12] 蒋瑶钰, 王彬, 陈祖明, 等. 基于COMSOL模型的黑土农耕地剖面冻结水热耦合分析[J]. 水土保持学报, 2023, 37(6): 187-193.
[13] 郭常凯, 王磊, 章青, 等. 基于COMSOL平台的水工混凝土冻融损伤模型试验验证与应用[J]. 混凝土, 2023(2): 40-48.
[14] 张东, 刘晓丽, 王恩志, 等. 非均匀多孔介质等效渗透率的普适表达式[J]. 水文地质工程地质, 2020, 47(4): 35-42.