土壤冻结对地源热泵系统中地热换热器热交换过程的影响.docxVIP

土壤冻结对地源热泵系统中地热换热器热交换过程的影响 地源热泵系统与其他空气调节系统相比具有许多优点。由于地层深处(约7m以下)温度常年维持不变,远远高于冬季的室外温度,而又大大低于夏季室外温度,因此地源热泵可以克服空气源热泵的技术障碍,且效率大大提高。地源热泵具有较高的“能效比(EER)”,供热时产生的热量是消耗电能的3倍以上。另外它还具有低噪音、占地面积少、不排放污染物、不抽取地下水、运行及维修费用低廉、设备使用寿命长等许多优点。 在供热工况下,当地源热泵系统蒸发器出口,即地热换热器入口的温度低于0℃时,换热器周围含湿土壤有可能冻结。由于地热换热器与周围土壤的换热计算及相变传热机理研究均十分复杂,因此关于土壤冻结对地热换热器与其周围换热过程的研究鲜有报道。 1 模型 1.1 城市地下埋管的热物性 为了简化问题使数学研究成为可能,作了以下一些假设: (1) 冻结的和未冻结的土壤热物性参数分别为常数。 (2) 土壤的孔隙分布均匀,孔隙中充满水分。 (3) 土壤中的传热方式为导热。 (4) 水结冰前后密度保持不变, 即忽略水在相变时的体积变化。 (5) 假设钻孔中回填材料的热物性参数与周围土壤保持一致,因为钻孔的直径与所讨论的传热区域相比是个小量。 (6) 认为地下埋管是一条线热汇/源,实际工程中,这种假设适用于数小时至数月的运行时间范围。 根据上述假设,可以将地下埋管与周围土壤换热过程描述为图1:图中下标s表示冻结,l表示未冻结。一条线热汇(地下埋管)埋设在初始温度为Ti的无限大的空间(大地)中,其中下标i表示初始时刻,自t=0时刻起,该线热汇开始从空间吸热,其单位长度热流强度为Q。空间中充满孔隙率为ε的土壤,孔隙中充满水。干土壤的物性参数分别为:密度ρg,导热系数kg,比热Cg;水与冰的物性参数分别为:kw,ki,ρ,Cw,Ci,冰点为Tm,潜热为L。取t时刻冻结锋面半径为S(t)。其中下标g表示土壤,w表示水,i表示冰,m表示冰点。 1.2 指数积分函数 冻结区域的当量导热系数、比热和密度为: ks=kg(1-ε)+εki(1)Cs=Cg(1-ε)+εCi(2)ρs=ρg(1-ε)+ερi(3)ks=kg(1?ε)+εki(1)Cs=Cg(1?ε)+εCi(2)ρs=ρg(1?ε)+ερi(3) 未冻结区域的当量导热系数、比热和密度为: kl=kg(1-ε)+εkw(4)Cl=Cg(1-ε)+εCw(5)ρl=ρg(1-ε)+ερw(6)kl=kg(1?ε)+εkw(4)Cl=Cg(1?ε)+εCw(5)ρl=ρg(1?ε)+ερw(6) 为了减少参数数量,本文采用无量纲形式。取:υ=(αl/αs)1/2,β=kl/ks,Stel=(Ti-Tm)Cl/(εL),其中r0是选取的特征长度,一般可选钻孔半径。另外引入无量纲变量:ξ=r/r0,Fo=αlt/r2020,Σ(Fo)=S(t)/r0,Um=(Tm-Ti)4πkl/Q,以及 Us(ξ?Fo)=[Τs(r?t)-Τi]4πklQ,Ul(ξ?Fo)=[Τl(r?t)-Τi]4πklQ(7)Us(ξ?Fo)=[Ts(r?t)?Ti]4πklQ,Ul(ξ?Fo)=[Tl(r?t)?Ti]4πklQ(7) 这样该问题的数学表述为: 控制方程:1ξ??ξ[ξ?Us?ξ]=υ2?Us?Fo1ξ??ξ[ξ?Us?ξ]=υ2?Us?Fo(冻结区) (8) 1ξ??ξ[ξ?Ul?ξ]=?Ul?Fo(未冻结区) (9) 边界条件∶ξ2β?Us?ξ=1(ξ→0?Fo＞0)(10)Ul(ξ?Fo)=0(ξ→∞?Fo＞0)(11) 初始条件:Ul(ξ,0)=0 (ξ0) (12) 冻结锋面条件:Us(Σ(Fo),Fo)=Ul(Σ(Fo),Fo)=Um(13) Σ′(Fo)=StelUm[?Ul?ξ-?Us?ξ1β](14) 解上述方程可得: 无量纲冻结锋面半径∶Σ(Fo)=2λ√Fo(15) 冻结区域与未冻结区域的无量纲温度分布: Us(ξ?Fo)=Um+β[Ei(-υ2ξ24Fo)-Ei(-λ2υ2)](16)Ul(ξ?Fo)=UmEi(-ξ2/(4Fo))Ei(-λ2)(17) 式中:Ei(x)=∫∞xe-SSdS,称为指数积分函数;λ是超越方程e-λ2Ei(-λ2)-e-λ2υ2Um=λ2Stel的解。 若不考虑土壤中水分的冻结,则土壤中温度分布为: U(ξ?Fo)=Ei(-ξ24Fo)(18) 2 结果与讨论 2.1 土壤冻结对热泵设计的影响 为了更清楚的说明问题,我们采用有量纲数据进行分析。取kw=0.5644W/(m·K),ki=2.24W/(m·K),kg=2.0W/(m·K),Cw=4.1868kJ/(kg·K),Ci=2.0927kJ/(kg·