竖直地埋管换热器传热温度场数值模拟与层换热理论.docxVIP

竖直地埋管换热器传热温度场数值模拟与层换热理论 1 家庭管道三维热态模型 1.1 数值模型的建立 目前，我国对地源热泵垂直地埋地沟换热器的三维数值模拟较少。基于二维热场数值模拟，忽视了以下换热问题。 1)地埋管换热器进、出水管在水平方向的相互影响; 2)地埋管换热器进水管和出水管在竖直方向和岩土换热的影响; 3)地埋管换热器之间的相互影响。 地埋管换热器如果按照此换热理论进行分析,得到的结论会与实际换热情况偏差很大,也不能正确指导地埋管换热器的换热计算和分析。因此,对于地埋管换热器的换热分析应该建立三维传热模型。 本文利用以离散化数值计算为基础的数值解模型建立了U形竖直地埋管换热器进、出水管内的流体传热和流动模型,以及岩土的传热模型。 岩土传热模型见式(1)。 式中ρe,ce,Te,λe,qv分别为岩土的密度、比热容、热力学温度、导热系数、单位容积的热产生率;τ为时间。 U形管内流体流动是湍流流动,本文采用K-ε双方程湍流模型。U形管内的湍流流动的控制方程包括连续性方程、动量方程、能量方程、湍流动能方程和湍流动能耗散率方程,所有这些方程满足如下通用形式的方程: 式中φ为通用物理量;ρ,v,J,S分别表示U形管内流动介质的密度、速度、扩散通量和源项。 计算边界条件与初始条件如下。 (1)地下岩土的初始温度均匀,且近似为半无限大的传热介质; (2)地下岩土的热物性是均匀的,且不随岩土温度的变化而变化,即具有常物性; (3)地表为第三类边界条件,设定空气温度是当地的月平均温度。 2)U形管内流体流动的边界条件 对于能量方程,U形管外壁面温度与岩土传热耦合求解。假定在壁面上水流不可渗透,对于湍流动能K方程,采用在壁面处扩散通量为零的边界条件: 式中n表示法向方向。 进口边界:包含温度、速度、湍流动能和动能耗散率,进口流体的温度由实际地源热泵机组运行中冷凝器或蒸发器的出水温度决定。对于速度,按U形管的流体流量计算给出。 出口边界:根据计算流体力学和数值传热学的方法,假定在出口平面上,流动充分发展,不存在边界下游对上游的影响,即在开口断面上网络节点的参数值对开口边界内邻近节点上的参数值无影响 在此基础上进行耦合求解,可以得到各种工况下的管群或单管地埋管换热器中水温动态分布及U形管换热器和周围岩土的动态温度分布。 1.2 初支集体的水温分布 为对比分析不同计算模型的计算差异,在数值计算分析基础上得到相同条件下的计算结果。图1显示的是利用二维计算模型得到的线热源状况下的3孔温度场分布,图2显示的是三维数值计算基础上的单孔温度场分布。 从图1,2可以看出,利用二维计算模型得到的温度场分布是以换热孔为中心的圆柱状分布,而利用三维计算模型得到的温度场分布是不规则的,根本原因在于进水管内流体的温度分布和出水管内流体的温度分布是不一致的,这就决定了两管的换热条件不一致,与换热孔及其周围岩土的换热条件也不一致。因此,进水侧和出水侧的温度分布差异较大。 由三维计算模型得到的夏、冬季不同工况下地埋管换热器(流量1.2 m3/h,埋深50 m)内的水温分布结果见图3,4。图3中工况1,2,3的进水温度分别为32,36,38℃。图4中工况1,2,3的进水温度分别为8,10,12℃。岩土初始温度均设定为18℃。 从图3,4可以看出,竖直地埋管换热器中水温分布呈现一定的规律性。在地下换热系统运行到一定时间后,在冬季和夏季工况下均呈现出进水管在一定深度范围内水温变化较小(夏季工况下深度在0~5 m范围,冬季状况下深度在0~20 m范围),水温变化小即表示换热器的进水管在该区域换热能力弱。进水管内流体跨越上述范围后,水温变化很明显,即达到一定深度后换热器内流体换热能力增强(夏季状况下深度在5~10 m范围,冬季状况下深度在20~25 m范围)。但是,当进水管流体继续下降达到某深度后,水温接近岩土温度,换热器内流体与周围岩土不发生工程意义上的换热。以上换热机理可以用层换热理论进行分析。 2 层交换理论 2.1 饱和传热层扩散 用L表示地埋管换热器的埋设深度,则L=L饱和+L换热+L未,其中L饱和为饱和换热层深度,L换热为换热层深度,L未为未换热层深度,它们是动态变化的。这3个层是统一划分的,包括了进水管和出水管。饱和换热层形成的原因是进水管带入热量的叠加导致该区域的岩土温度不断升高,达到一定程度时,岩土温度和进水温度接近,无工程意义上的换热发生;而换热层是进水管穿越饱和换热层后,与该区域岩土温度形成了温差,大量换热在该区域发生;当换热达到一定程度后,进水管中最终的水温与岩土初始温度接近,岩土已无工程意义上的换热能力,该层为未换热层。如果换热量大,换热时间长,L饱和大,即饱和层深度逐渐向换热器下部推移。若运行参数不变,L换热不变,但换热层位置要随饱和换热层的发展而向下移动