X img





作者投稿 专家审稿 编辑办公 编委办公 主编办公

冰川冻土 ›› 2021, Vol. 43 ›› Issue (4): 948-963.doi: 10.7522/j.issn.1000-0240.2021.0054

• 第二次青藏高原综合科学考察研究 • 上一篇    下一篇


胡锦华1,2(), 陆峥1,2, 仝金辉1,2, 李小雁1,2, 刘绍民1,2, 杨晓帆1,2()   

  1. 1.北京师范大学 地理科学学部 地表过程与资源生态国家重点实验室, 北京 100875
    2.北京师范大学 地理科学学部 自然资源学院, 北京 100875
  • 收稿日期:2021-02-08 修回日期:2021-06-20 出版日期:2021-08-31 发布日期:2021-09-09
  • 通讯作者: 杨晓帆 E-mail:hujinhua@mail.bnu.edu.cn;xfyang@bnu.edu.cn
  • 作者简介:胡锦华,博士研究生,主要从事寒区水文模拟研究. E-mail: hujinhua@mail.bnu.edu.cn
  • 基金资助:

Simulating thermo-hydrologic processes in cold region soil system: a computational fluid dynamics study

Jinhua HU1,2(), Zheng LU1,2, Jinhui TONG1,2, Xiaoyan LI1,2, Shaomin LIU1,2, Xiaofan YANG1,2()   

  1. 1.State Key Laboratory of Earth Surface Processes and Resource Ecology,Faculty of Geographical Science,Beijing Normal University,Beijing 100875,China
    2.School of Natural Resources,Faculty of Geographical Science,Beijing Normal University,Beijing 100875,China
  • Received:2021-02-08 Revised:2021-06-20 Online:2021-08-31 Published:2021-09-09
  • Contact: Xiaofan YANG E-mail:hujinhua@mail.bnu.edu.cn;xfyang@bnu.edu.cn



关键词: 寒区水文学, 地下水模型, 水热耦合过程, 计算流体力学, 基准测试算例


Cold regions, which supply freshwater to downstream communities, have been significantly disturbed by anthropogenic climate change that is pronounced to continue in future. As unique and critical elements of the cold regions, frozen soil plays a vital role in cold region hydrology. Although frozen soil is considered relatively impermeable to groundwater flow, the freezing and thawing of the frozen soils under seasonal and climate change may influence the hydrologic components such as surface water infiltration, groundwater recharge, and hydrogeologic connectivity that are important for water resources. Under the pronounced climate warming, increased soil temperature expedites the thawing of frozen soil in cold regions. Therefore, understanding the complex thermo-hydrologic processes in frozen soils is one of key issues in studying cold region hydrology. However, due to the harsh environment, it is difficult to conduct hydrometeorological and geophysical field observations yet collect gap-free, high-resolution datasets. With the recent development of computational methods and resources and based on the observational data in recent decades, it is essential to enable modeling and simulation of the thermo-hydrologic processes in frozen soils to study cold region hydrology. The thermo-hydrologic processes in frozen soils is fundamentally interpreted into multi-phase flow and heat transfer in porous media with phase change. The current challenges in simulating such processes, especially in cold regions include but not limited to the interactions among liquid water, vapor, ice and soil grains and the phase change between ice and liquid water, which introduce high non-linearities and uncertainties in the processes. However, traditional distributed or hydrogeological models were limited in terms of resolution, accessibility and capability of solving the coupled thermo-hydrologic processes with high fidelity. To account for these feedbacks as far as possible, a physically-based, massively-parallel, fine-resolution and fully-saturated cryo-hydrogeological model was in-house developed in OpenFOAM?http://www.openfoam.com), an open-source computational fluid dynamics (CFD) solver, which has been widely used in groundwater modeling and computational hydrology. The cryo-hydrogeological model, named darcyTHFOAM, was formulated by mass conservation, Darcy’s law, energy conservation and soil freezing function. To ensure the convenience and joy for users, an interface-based software was developed correspondingly using Python 3.5 (www.python.org) and designed with QT Designer (www.qt.io) in Linux environment. In general, rigorous validation and benchmarking of the numerical models are prerequisites and critical for their validation and applications. Hence, a series of benchmarking simulations under fully-saturated condition were performed to validate the current model, including the classical analytical solution of Stefan’s equation, two community-recognized benchmark cases (under thawing) using synthetic porous media samples in the Interfrost Project (wiki.lsce.ipsl.fr/interfrost) and laboratory freezing experiment with a 25 cm-long soil column. Simulated results include the evolution of temperature, liquid water content and ice content. Frozen depth, temperature profiles at fixed points or lines, the minimum of the temperature field were selected for quantitative comparisons. Simulated results were in excellent agreements with those obtained from the analytical solution, other cryo-hydrogeological models and experimental data, respectively, which demonstrated the reliability and accuracy of the current model. The study revealed that the proposed model is capable of simulating coupled thermo-hydrologic processes under fully-saturated condition in cold regions with high spatial resolutions and efficiency. Overall, the in-house developed cryo-hydrogeological model is expected to serve as a powerful tool for studying subsurface hydrology in cold regions. Further code developments involve coupling the surface processes, especially snow, and transport processes (e.g. contaminants) for studying groundwater-surface water interactions and hydro-biogeochemical processes in cold regions.

Key words: cold region hydrology, groundwater model, coupled thermo-hydrologic processes, computational fluid dynamics, benchmark case


  • P642.14