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