地源热泵地埋管数值模拟研究报告.docVIP

天津联合利华地源热泵埋管 数值模拟研究报告 2014年1月 目 录 一、模拟方法 二、模拟参数说明 三、区域温度变化模拟研究 四、冷却塔容量选择分析 五、主要结论与建议 一、模拟方法 本次模拟所建立的模型是以圆柱源理论模型为依据的。该模型是本质上是基于线源模型，是线热源理论模型的一种延伸与发展，而线热源模型是圆柱热源模型的一种简化形式。该模型将U型埋管用当量直径等价为一根圆管，对于恒热流情况，可以给出其精确解。 具体解法如下： 1.1恒热流条件下钻孔壁温的确定 由圆柱源理论模型可以直接得到在恒热流条件下，土壤任意点处任意时刻时的温度值，其圆柱源解析解为： （1-1） 式中，为时刻半径为r处的土壤与远边界温度之差，℃。当r为孔洞半径时，即为孔洞壁面与土壤远边界的温度差；为埋管热流量，W/m，管吸热时为正，放热时为负； 为土壤导热系数，W/m℃；，分别为时刻半径r处的土壤温度及远边界的初始温度，℃；为理论积分解G函数，Fo=4at/d2。当r为钻孔半径时，对应于p=1。此时，G（Fo, 1）的表达式如下： （1-2） 把式（1-2）代入式（1-1）得到： （1-3） 式中，为钻孔壁温，℃；为孔洞壁面与土壤远边界的温差，℃。 1.2变热流条件下钻孔壁温的确定 在实际运行过程中，地下埋管从土壤中的取放热量是随建筑负荷的变化而动态改变的。这样随时间变化的热流可以看作是一系列连续作用于孔洞中的矩形阶跃负荷，如图1-1示。 图1-1 用矩形阶跃负荷近似连续热流 每一时刻的埋管阶跃负荷都会给孔洞壁面产生一个温度响应，则在某一时刻末对埋管孔洞的温度响应可应用叠加原理得到。根据叠加原理，将某一时刻以前各时间段埋管取放热量对孔洞周围土壤温度响应的作用都叠加到该时刻，如图1-2所示。 图1-2 负荷叠加说明图 对于具有四个不同埋管热流量时间段的情形，t4时刻末的土壤远边界与孔洞壁面间的温差应该由q1,q2,q3,q4共同决定，根据叠加原理有： （2-1） 可以看出，起始时刻的热流q1对整个运行段t4的温度响应均起作用，随后第二个时间段的热流q2只对随后的三个时间段（t4-t1）的温度响应有效，但此时必须抵消前一时间段q1热流对该段的影响，即该时间段作用的等效热流值应为（q2-q1）。依此类推。这样可以得到第tn时刻末土壤远边界与孔洞壁面的温差为： （2-2） 1.3 埋管内流体平均温度的确定 基于传热学理论，在已知孔洞壁面温度时只要能够计算出埋管内流体至孔洞壁面间的传热热阻，便可由下式确定出埋管内流体的平均温度，即 （3-1） 式中，为U型管内流体的平均温度，℃，定义为U 型管进出口流体温度的算术平均值；为孔洞壁面温度，℃。为孔洞热阻，m℃/W 图1－3 U型埋管孔洞结构 于是有： （3-2） 式中：，，分别为管内流体对流换热热阻、管壁导热热阻及灌浆热阻，其计算表达式分别对应于等式右边的三项。，分别为U型管的内、外直径，m；为U型管的当量直径，m；为管材的导热系数，W/m.℃为灌浆材料的导热系数，W /m oC; 为流体与管内壁面间的对流换热系数，W/m2 oC （4-1） 式中，为埋管出口流体温度，℃；H为钻孔深度，m；m为循环流体的流量，kg/s；为循环流体的质量比热容，kJ/kg.oC。 1.5 孔洞以外土壤温度的确定 根据以上的计算我们得到了变热流条件下孔洞壁面温度的变化情况。孔洞以外土壤的传热过程是非稳态的。其传热过程遵循半无限大圆柱形物体的非稳态导热规律。传热方程为： （5-1） 式中，为土壤密度，kg/m3；为土壤的比热容，kJ/kg.℃；为距离孔洞中心半径为r处，时刻的土壤温度与无穷远处土壤温度的差值；为土壤的导热系数，W/m.oC。由于在整个传热过程中，土壤的可以看成是定值，所以可以提到括号外面来。方程变为： （5-2） 式中，,为土壤的导温系数，m2/s。 1.6 程序模拟界面 1.6.1 参数输入界面： 图1－4 输入参数界面 1.6.2埋管出口流体温度计算界面 图1－5 出口流体温度输出界面 1.6.3土壤温度场计算界面