套管式地下换热器传热模型及其传热分析.docVIP

离散化数值模型 Eskilson方法 Eskilson模型基于有限差数值离散方法求解。考虑地表面温度及钻孔之间的相互影响，考虑热泵机组间歇运行工况，考虑土壤土质的影响。由于数值求解需要大量时间，此模型采用无因子反应系数方法，即用g函数法来确定钻孔周围的温度分布。每个g函数代表一个钻孔域的结构。对于多个钻孔情况，采用单个钻孔温度相应叠加的方法得到。 该模型可以求解支管内流体温度和热流情况。尤其对于判断系统长期的运行的稳定性，保证冷热负荷的平衡有很好的效果。但因为g函数没有考虑钻孔局部的几何形状，所以当时间步长较小的时候，此模型不能提供准确的温度响应。 V.C.Mei方法（有限差分法） 1986年，V.C.Mei提出了三维瞬态远边界传热模型，该理论是建立在能量平衡基础上，由系统能量平衡结合热传导方程构成。该模型的假设为：⑴岩土是均匀的；⑵埋管内同一截面流体温度、速度相同；⑶大地土壤热物性参数不变；⑷不考虑热湿迁移影响⑸忽略埋管与岩土的接触热阻。在这些假设的基础上，可以对各截面的径向传热建立方程，通过截面推移得到三维温度场。 热泵运行时，流体从内管进入，外腔流出，与内外管壁发生对流换热，如下图所示： 套管式换热器模型 进入内管的流体与内管内壁的传热： （2-1） 内管管壁的导热： （2-2） 环腔内流体与套管管壁的传热; （2-3） 外管管壁的导热： （2-4） 大地中的导热： （2-5） 初始条件（t=0）：(Z)=(Z)=(Z)=(Z)=(Z), 进口处流体表面的初始条件（Z=0）：(t,0)=f(t),是已知时间的函数； 换热器管底末端的边界条件（Z=L）： (t,L)=(t,L) （2-5） 内管中流体与内管管壁处的边界条件： （2-6） 环腔内流体与内管外管壁处的边界条件： （2-7） 环腔内流体与套管管壁处的边界条件： （2-8） 外管与大地交界处的边界条件： （2-9） 假设的远边界处边界条件：，（） 式中： ------内管和环腔内的流体速度，m/s； Z------沿管长方向坐标，m； ------内管和环腔内的流体温度，℃； ～------径向坐标，m；r------径向坐标，m； ，，------内管管壁、外管管壁和大地的导热系数，w/（m·℃）; ------流体的密度，Kg/； ------流体的比热，J/（Kg·℃）； ，，------内管管壁、外管管壁和埋管附近大地的温度，℃； ，，------内管管壁、外管管壁、土壤的导温系数，； ，，------内管中流体与内管内壁、环腔中流体与内管外壁及套管内壁的对流换热系数，； t------从运行开始计算的时间，s。 3.Rottmayer方法（有限差分法） Rottmayer方法建立了准三维的传热模型，考虑流体温度沿埋管深度方向的变化，考虑间歇运行工况，考虑灌浆材料的影响。在水平方向在极坐标下划分扇形网格，不考虑埋管深度方向导热的影响，利用能量守恒建立二维传热模型，并用显式差分格式离散求解。U型管处用扇形近似，使扇形的周长与U型管周长相同，并提出采用几何因子来修正由扇形区域代替圆管产生的误差。 4.Yavuzturk方法（有限容积法） Yavuzturk等建立的方法书瞬态二维隐式有限容积模型，时间采用向后差分，空间采用中心差分的离散方法。与Rottmayer方法相同，Yavuzturk方法也采用扇形面积代替圆管，水平方向在极坐标下划分扇形网格。对于U型埋管，温度场以两支管中心连线对称，因此水平方向只对半圆形区域划分网格，研究土壤传热。为了地源热泵系统能耗分析和混合式地源热泵的设计，作者利用此模型计算无因此温度影响因子。