地埋管换热器能效系数影响因素的三维数值模拟研究.docxVIP

地埋管换热器能效系数影响因素的三维数值模拟研究 0 地埋管传热模型 地埋地源热泵系统利用放置在岩石体上的垂直地埋地幔换热器，实现清热装置和岩石体之间的换热，并利用低附加值地下能耗。如果地埋管换热器设计不当,将会很大程度上损害地埋管地源热泵理论上的潜在优势。Cane和Forgas曾经就北美地埋管标准设计方法,即国际地源热泵协会(IGSHPA)推荐的方法进行了评估,基于该简化方法的地埋管换热器设计管长都超过实际所需的10%～33%,这种比例足以使地埋管地源热泵技术的应用在短期内丧失经济性优势。因而,要提高地埋管地源热泵系统的经济性和运行的可靠性,就需要建立较精确的地埋管传热模型,并分析地埋管换热器在不同条件下的换热过程。 U形地埋管换热器与周围岩土体的传热过程实际上是一个多层介质的导热过程,多为无渗流的传热过程,会受到流体流动特性、岩土体物性、埋管几何结构及地埋管换热负荷变化等诸多因素影响的复杂过程。地埋管地源热泵系统运行过程中,在进口温度一定的条件下,地埋管循环流体出口温度不仅反映地埋管与岩土体的换热效果,同时极大地影响着热泵机组的运行效率,是地埋管换热过程的一个重要能效特性参数。本文将基于能效系数概念,利用文献无渗流条件下的地埋管传热模型模拟分析地埋管传热特性及其在各种影响因素下的变化规律,为地埋管换热器的结构优化及运行模式合理化奠定理论基础,从而给地埋管地源热泵系统的推广和应用提供技术参考。 1 计算嵌入热交换器换热器的变换值的条件 1.1 埋地管道的结构和岩石的物理条件 为了研究地埋管换热器在无渗流条件下的换热特性,本文针对地埋管地源热泵夏季空调工况进行计算分析,采用的基本模拟条件如表1所示。 1.2 模拟无穷远边界值 埋设于岩土体中的单U形地埋管实际上是位于无限大介质中,无边界条件限制。要在数值模拟中完全实现工程应用的地埋管传热物理模型相当困难,因此必须选择适当的边界值以保证在模拟计算中边界不受地埋管排热或放热的影响。 无渗流条件下的地埋管传热是向四周呈圆形辐射状传热,以半径rs表示模拟无穷远边界值,一般取值范围为几百倍到几千倍的地埋管外径,考虑对称性,计算域可取圆的一半,如图1所示。rs越大,边界发生热扰动的可能性越小,计算的精确度就越高,但这对计算机的性能要求较高,甚至无法计算,而rs取值较小的话就会导致计算结果偏离实际,因而有必要分析地埋管传热过程对数值模拟边界的要求。 根据文献对量纲一边界热扰动因子影响因素的分析可知,地埋管传热过程模拟的无穷远边界取值主要依赖于岩土体的热物性特别是变化幅度较大的导热系数以及模拟时间。为保证数值计算结果的正确可靠,需要根据模拟时间来选取适当的计算边界值。在大多数情形下,夏季或冬季空调系统运行时间一般为2 000～5 000 h,rs取1 000ro～1 500ro可以满足模拟要求。 2 地埋管换热器能效系数 地埋管换热器作为热泵机组和岩土体之间的传热通道,其换热性能将不可避免地受到热泵机组的取热/排热负荷、岩土体状况等诸多因素的影响。地埋管换热器在夏季放热状态下,较高温度的冷凝器出水进入到地埋管换热器的下降管,与土壤换热后温度下降;地埋管换热器在冬季吸热状态下,较低温度的蒸发器出水进入到地埋管换热器的下降管,与土壤换热后温度上升。从理论上来说,地埋管出水温度可最大限度地接近土壤初温,此时的地埋管换热量最大,但实际由于岩土体的蓄热和热堆积效应并不可能实现。为量化地埋管与岩土体换热的充分程度,借鉴传统换热器效能的概念,定义U形地埋管换热器能效系数E为地埋管换热器实际换热量Q与最大理论换热量Q′的比值,其计算式为 E=QQ′=Gc(tin?tout)Gc(tin?t0)=tin?touttin?t0(1) 式中G为埋管内流体流量,kg/s;c为埋管内流体比热容,J/(kg·℃);tin为地埋管换热器进水温度,℃;tout为地埋管换热器出水温度,℃;t0为岩土体初始温度,℃。 地埋管换热器能效系数是一个量纲一的瞬时变化量,其取值范围为0～1。在给定的地埋管进口温度和岩土体初始温度条件下,能效系数与地埋管的出口温度密切相关,表征了U形地埋管与周围岩土体连续换热后管中流体出口温度能够达到的最低(夏季工况)或最高(冬季工况)的能力,以及流体进出口温差达到的最大值,也即是进入到热泵机组冷凝器(夏季工况)或蒸发器(冬季工况)中流体的温度状态,影响着机组的运行效率,与岩土体的热物性及地埋管换热器的设计参数相关。 地埋管单位井深换热量可改写为 Q=Gc(tin?tout)H=EGc(tin?t0)H(2) 如果流量G保持不变,地埋管换热器能效系数在一定程度上反映了地埋管的换热效果。本文将有针对性地分析地埋管换热器能效系数的影响因素及其变化规律,以利于增强地埋管地源热泵系统的换热性能。 2.1 岩石物理