冻结条件下土壤水热耦合迁移数值模拟的改进3.PDF

( ) 清 华 大 学 学 报 自 然 科 学 版 1830　 　1997 年 第 37 卷 Journal of T singhua U niversity (Sci T ech) 第 8 期 第 62～ 64 页　 冻结条件下土壤水热耦合迁移数值模拟的改进 尚松浩, 　雷志栋, 　杨诗秀 清华大学 水利水电工程系, 北京 100084 文　摘　根据冻结条件下土壤水分运动和热量迁移的基本 分、温度的变化自动调节时间步长, 加速了迭代计算 方程, 推导出了冻土水热耦合方程, 并构造了采用变化的空 速度。同时分析了阻抗系数对冻土水分、温度变化的 间步长和自动调节时间步长的全隐式有限差分计算格式, 使 影响。 模型求解中的迭代计算过程大为简化, 计算速度有很大提 高。用本模型对张掖壤土的室内冻结过程进行了模拟, 土壤 1　数学模型 水分、温度的计算值和试验值吻合较好。计算结果表明, 该方 考虑均质各向同性的垂直土柱的一维冻结问 法不仅能模拟冻结过程中土壤水分、温度的变化, 且在计算 题。在忽略水汽迁移、热量对流作用的情况下, 非饱 中土壤含水量较低时不需要在冻结区引入冰的阻抗系数的 和土壤水分、热量运动方程分别为[ 5 ] 概念, 较好地改进了已有的数值模拟, 为模拟各种不同条件 u i i = D ( ) u - K ( ) - ( 1) 下的土壤冻结过程创造了条件。 t z u Z u t w 关键词　土壤冻结; 水热迁移; 水热耦合方程; 有限差分法 T T i Cv = + L i (2) 分类号　 641. 2 t Z Z t P 式中 u , i 分别为土壤未冻水、冰的体积含量, T 为 ( 土壤冻结是一种常见的自然现象, 在我国北方 土壤温度, t, Z 分别为时间、空间坐标 垂直向下为 ) ( ( ) 正 , u , u 分别为非饱和土壤水分扩散率、 地区永久性冻土和季节性冻土区的面积占国土面积