非饱和土壤高温热湿迁移的数值模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号： 083029 非饱和土壤高温热湿迁移的数值模拟 王华军，齐承英 （河北工业大学能源与环境工程学院，天津 300401） (Tel: 022 Email: huajunwang@126.com) 摘要：对非饱和土壤高温储热热湿迁移过程进行了模拟计算。结果表明，土壤热湿迁移中的湿度场明 显滞后于温度场，且与土壤特性有很大关系。对于水力传导性较差的土壤，在高温储热初期，靠近热 源的地方容易产生湿份聚集，使得湿度曲线出现一个短期峰值。与低温储热相比，土壤高温储热时湿 度迁移对温度场的影响较大，温度场呈整体降低趋势，变化幅度顺序为：砂土壤土粘土。 关键词：土壤，高温热储，热湿迁移，数值模拟 0 前言 近年来，土壤储热技术得到了较为广泛的关注和应用，并成为一种提高系统能源利 用率的有效方式[1] 。按照热源温度差异，土壤储热可分为低温（低于40 ℃）和高温（40 －90℃）两种类型。对于低温土壤储热（如地源热泵），通常将地下换热器与周围土壤 之间的热交换视为纯导热过程。实践表明，基于经典导热理论的线/柱模型可以满足工程 应用的基本要求。尽管如此，土壤储热与热湿迁移问题仍倍受关注和争议[2-3] 。 与低温储热相比，高温土壤储热过程中的温度场和湿度场通常显得更为关联。一方 面，土壤的热物性（如导热系数、比热等）随温度变化较大；另一方面，由于温度梯度 较大，使得土壤中水分、空气甚至水蒸气对流扩散运动趋于强烈。此时的土壤储热模式 实质上是一个较为复杂的热湿耦合传热传质过程。Reuss[4]等研究表明，地下换热器附近 土壤热湿迁移现象可能导致局部干涸和裂隙，从而影响传热速率。因此，就土壤储热系 统优化设计而言，研究土壤高温热湿迁移规律具有重要意义。 1 计算模型 1.1 控制方程 多孔土壤中的一维水分运动满足流动方程 ∂θ(h) ∂ ⎡ ∂h ∂T ⎤ K K S h （1） ∂τ ∂x ⎢⎣ ( ∂x =+cosα) + LT ∂x ⎥⎦− ( ) 其中， h 为水压头，K 和KLT 土壤固液相的水力传导率。α 为流动方向与重力方 向夹角（0º 为垂直流动，90º 为水平流动），S 为源项。 [5] 土壤的水力传导率K (h) 可采用 Brooks-Corey 模型 ，即 2/ n+3 ⎧⎪αh -n h -1/α K K S Se ⎨ （2 ） s e 1 h ≥-1/α ⎪⎩ n α 其中，K 为饱和水力传导率， 为粒径分布指数， 为持水系数。 s 多孔土壤热湿迁移中传导-对流-扩散过程的能量守恒过程满足 ∂T ∂ ⎡ ∂T ⎤ ∂T c (θ) λθ( ) =−c q （3 ） p ∂ ∂τ x ⎢⎣ ∂x ⎥⎦ w ∂x 土壤的容积比热c (θ) 可以通过其固相、有机质和液相的体积比得到，即