竖直地埋管动态负荷下换热特性三维数值分析.docxVIP

竖直地埋管动态负荷下换热特性三维数值分析 0 基于三维数值模拟的土壤温度分析 垂直地埋式换热器的传热模型主要包括无限线性光源模型、无限长圆柱形水源模型和有限长度喷泉模型。这3种模型描述的都是地埋管换热器在恒定负荷下的温度响应。 地埋管换热器的传热过程是非稳态过程。由于计算时刻的温度响应是过去时刻负荷综合作用的结果, 要准确评价地埋管的性能必须进行逐时计算。三维数值模拟能考虑更多和更真实的边界条件, 物理过程的假设也较少, 其预测结果能更好地反映地埋管的真实换热性能。然而, 全尺寸三维数值模拟方法会受到计算机性能的限制, 将耗费大量的计算资源和时间, 在工程中很难实现。相比数值计算, 解析计算通用性强, 建模简单, 对计算资源的要求较低。在地埋管换热器的3种解析传热模型中, 无限长线热源模型最为简单, 使用也最为广泛, 若将该模型用于变负荷下的地埋管换热器性能分析, 可节约大量的计算时间。 对于动态负荷问题, 工程中大都采用Eskilson提出的脉冲负荷叠加原理, 即将连续动态负荷近似为由不同起始时刻的恒定负荷叠加构成。本文基于无限长线热源模型, 应用负荷叠加原理, 建立了竖直地埋管动态负荷换热二维解析计算模型, 并分别应用该模型和已建立的全尺寸三维地埋管换热数值计算模型对地埋管换热性能进行了全年计算, 并对两者计算结果进行了对比。 1 在分析变量负载的性能时，计算模型 应用无限长线热源模型和脉冲负荷叠加原理, 得到二维单U形地埋管变负荷运行换热特性解析计算模型。 1.1 巷道尾数tr,的岩土导热系数 将竖直U形地埋管看成是置于温度均匀的土壤中, 且为恒热流的无限长线热源, 其周边岩土区域无限大, 从零时刻启动, 任意时刻的岩土温度分布为 式中t (r, τ) 为与线热源中心距离为r处的岩土温度, ℃;τ为时间, s;t0为岩土初始温度, ℃;q为线热源散热强度, W/m;k为岩土导热系数, W/ (m·℃) ;, 其中a为岩土热扩散率, m2/s;β为积分变量。 令, 对式 (1) 进行积分变换, 得 由于实际工程问题通常为动态负荷问题, 然而式 (1) , (2) 中的q为常量, 因此提出由此推导变负荷运行换热特性解析计算模型。 1.2 动态负荷分析的计算模型 1.2.1 不同时段开始的连续负荷 近似认为随时间连续变化的动态负荷在较小的时间段内是稳定的, 每个时间段内的稳定负荷可用这个时间段的动态负荷的平均值表示, 这样, 可将动态负荷分解为阶梯状的离散负荷, 如图1a所示, 将一个随时间连续变化的负荷分成5个连续时间段的离散负荷, 这5个连续时间段用 (0, τ1) , (τ1, τ2) , …, (τ4, τ5) 表示, 每个时间段内的负荷值为实际负荷的平均值, 用q1, q2, …, q5表示。 在此基础上, 设置不同时刻开始的连续负荷Q1至Q5, Q1=q1, 从0时刻开始到τ5;Q2=q2-q1, 从τ1时刻开始到τ5;…;Q5=q5-q4, 从τ4时刻开始到τ5。对这些连续负荷进行叠加即得到所对应的离散负荷q1, q2, …, q5, 见图1b。 1.2.2 岩土温升lb 基于无限长线热源计算模型, 即式 (2) , 可以得到从τi-1时刻开始的第i个负荷Qi持续到τn时刻所引起的岩土温升: 式中τn为终点时刻;rb为地埋管钻孔半径, m。 应用叠加原理, 将不同时刻开始的连续负荷所引起的岩土温升叠加起来, 即得到单根竖直地埋管在动态负荷下周围岩土温度场的计算式, 为 1.2.3 流体温度tf 由于钻孔内的热容量较钻孔外小得多, 在同一时刻, 假设地埋管的换热量等于钻孔内壁的换热量。考虑钻孔内热阻, 即可得到U形管内流体的平均温度: 式中tf(τn) 为τn时刻U形管内流体的平均温度, ℃;t (rb, τn) 为τn时刻地埋管钻孔壁温度, ℃;q (τn) 为τn时刻地埋管延米换热量, W/m;Rb为钻孔内热阻, m2·℃/W, 根据文献计算。 由式 (4) , (5) 可得τn时刻地埋管内流体的平均温度为 从而, 根据热平衡原理, 并假设单U形地埋管的上、下竖管换热量相同, 可以得到τn时刻地埋管进、出水温度分别为 式 (7) , (8) 中tf, o为U形管出口流体温度, ℃;tf, i为U形管进口流体温度, ℃;c为循环流体的比热容, J/ (kg·K) ;m为循环流体的质量流量, kg/s;H为埋管深度, m。 2 数值模型验证 分别应用以上得到的解析计算模型和文献建立的地埋管全尺寸三维数值计算模型 (通过了实验验证) , 对单U形竖直地埋管进行全年动态负荷计算, 对地埋管进、出口流体温度和周围岩土温度进行预测, 并对两种方法计算的结果进行对比。 2.1 单根单u形地埋管逐时负荷 以西安地区某一办公建