土壤冻结对地源热泵地埋管换热器换热特性的影响.docxVIP

土壤冻结对地源热泵地埋管换热器换热特性的影响 0 地埋管周围土壤的冻结 与传统的空调系统相比，内部空调系统（pgi）获得了更好的推广和应用，因为它具有良好的节能性、稳定性和环境性。地源热泵研究的核心课题之一就是要解决地埋管换热器与周围土壤间的换热问题。在冬季供热工况下,特别是在北方寒冷地区,当埋管从地下连续取热时,地源热泵蒸发器出口(埋管进口)流体温度会低于0℃,从而导致埋管周围含湿土壤有可能冻结,这在一定程度上直接影响了埋管的换热特性。由于地埋管与周围土壤间传热及相变传热机理的复杂性,关于地埋管周围土壤的冻结及其对埋管换热特性影响的研究很少有报道。 土壤的冻结过程是一个复杂的含湿多孔介质体的非稳态相变传热过程,由于土壤水是含有各种杂质的非纯净物质,从而使得土壤水的相变过程不是严格地在某一温度下发生,而是在一个小的温度范围内完成。相变问题又称为Stefan问题,其特点是计算区域内存在一个位置随时间移动的两相界面,在该界面上放出或吸收潜热,因此亦称为“移动边界问题”。此类问题除极少数简单情况能进行分析与近似法求解外,主要依靠数值解法。迄今为止,研究者已提出了多种求解相变问题的数值计算方法,如变时间步长法、有限差分-有限元法、焓法和显热容法等。其中以显热容法较为简单实用、易于进行多维推广,并特别适用于求解发生在某一温度范围的相变传热问题。文献比较了显热容法与精确解和有限元方法,结果表明显热容法的计算精度令人满意,在处理多维问题时优势明显。文献对相变潜热的不同分布情况进行了研究,分析表明显热容法是处理较复杂相变问题的有效方法。文献对显热容法和焓法进行了对比研究,表明显热容法更为简便。 基于显热容法思想,本文重点对地源热泵地埋管换热器周围土壤的冻结特性及其对埋管换热的影响进行研究,探讨土壤含水率、土壤原始温度、土壤热扩散率及Stefan数对地埋管周围土壤温度分布及土壤冻结特性的影响。 1 计算模型 1.1 土壤冻结网络的热物性参数 由于地埋管周围土壤的冻结是一个极其复杂的相变传热过程,为了简化问题以使数学研究成为可能,需作以下近似简化处理:1)土壤为各向同性、均质、孔隙分布均匀的刚性含湿多孔介质体,且孔隙中充满水分;2)土壤初始温度均匀一致,且不考虑深度方向的传热;3)水结冰前后密度保持不变,即忽略水在相变时的容积变化;4)忽略土壤热湿迁移耦合作用及自然对流效应的影响;5)冻结的和未冻结的土壤热物性参数均为常数;6)土壤的冻结相变过程发生在一个小的温度范围内,且认为在土壤相变中存在三个区域:冻结区、未冻结区及介于两区之间的模糊两相区;7)假设钻孔中回填材料的热物性参数与周围土壤保持一致,因为钻孔的直径与所讨论的传热区域相比是个小量;8)认为地埋管是一条线热汇/源,实际工程中,这种假设适用于数小时至数月的运行时间范围。 根据以上简化假设,可将地埋管与周围土壤的换热过程描述为:一条线热汇(地埋管)埋设于初始温度均匀为t0的无限大空间(土壤)中,土壤的含水率为ε。自τ=0时刻起,该线热汇开始从空间吸热,其单位长度热流强度为q。随着吸热过程的进行,埋管周围土壤的温度逐渐降低,并沿径向方向依次呈现出三个相区:冻结区、模糊两相区及未冻结区。如图1所示,图中下标s表示冻结,1表示未冻结,m表示模糊两相区;tm为冻结中心温度,△t为土壤冻结相变温度区间的一半;S1(τ),S2(τ)及Sm、(τ)分别为土壤完全冻结、开始冻结时及土壤处于相变温度区间中心温度时的锋面半径;r为径向半径。 1.2 显热容法的基本原理 由以上物理模型可得如下各相区的控制方程。 冻结区 未冻结区 模糊两相区 相变界面 式(1)-(6)中Cs,C1,Cm分别为冻结区、未冻结区及模糊两相区的土壤容积比热容,J/(m3·℃);λs,λ1,λm分别为冻结区、未冻结区及模糊两相区的土壤导热系数,W/(m·℃);f5为冻结率;L为单位容积土壤的相变潜热,J/m3。 分析可以看出,该导热问题较之一般的变系数导热问题要复杂得多,主要在于求解区域中存在着随时间移动的相变界面。为了处理这个问题,可以采用把它视为“单相”的非线性导热问题的处理方法——显热容法,其思想是:由于土壤的冻结相变不是严格地在某一特定的温度下发生,而是在一个小的温度范围内发生,因此可以把土壤中水的相变潜热看作是在这个小温度范围内有一个很大的显热容,显热容的大小由相变潜热和相变温度范围来决定。这样就可以把上述分区描述的控制方程及界面守恒条件转化成在整个区域上适用的一个单相非线性瞬态导热问题,不需要跟踪两相界面,待求出温度场后,再根据温度确定相变界面的位置。 根据显热容法,认为模糊两相区内导热系数随温度呈线性分布,热容量按平均法取值,且相变潜热不随温度变化(文献研究表明用相变温度范围内的平均潜热代替实际相变潜热且视为常数是可行的