2025-11-14T05:04:11.127420

A combined thermal-resistance-capacity and finite-element model for very fast and accurate short- and medium-term simulations of single U-tube borehole heat exchangers

Zanchini, Zanchini, Naldi
An accurate design of a ground-coupled heat pump system requires the knowledge of the outlet fluid temperature from the borehole heat exchangers (BHEs), both in the short and long term. This paper fucuses on the short and medium term. In this time range, either 3D finite-element simulations or Thermal Resistance Capacity Models (TRCMs) can be applied. The former can yield very accurate results but require long computation times. The latter are much faster but cannot be fully precise, because they require simplifying assumptions. In this paper, we present a new method for the short-term and medium-term simulation of single U-tube BHEs, which combines the speed of TRCMs and the accuracy of finite-elements simulations. The method uses a TRCM to estimate the thermal response of the BHE, then corrects the results by interpolation with a dataset, which has been produced by running 54 finite-element simulations in various configurations. The model is implemented in a C++ program, available at the open-source online data repository of the University of Bologna. The program provides, within two seconds, the time evolution of the inlet, outlet and mean fluid temperature, of the mean BHE surface temperature, of the 3D and the effective borehole thermal resistance. It can be easily connected to long-term simulation tools to obtain the full-time-scale thermal response of a bore field.
academic

A combined thermal-resistance-capacity and finite-element model for very fast and accurate short- and medium-term simulations of single U-tube borehole heat exchangers

基本信息

  • 论文ID: 2510.14421
  • 标题: A combined thermal-resistance-capacity and finite-element model for very fast and accurate short- and medium-term simulations of single U-tube borehole heat exchangers
  • 作者: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (University of Bologna)
  • 分类: physics.app-ph
  • 发表时间: 2024年
  • 论文链接: https://arxiv.org/abs/2510.14421

摘要

地源热泵系统的精确设计需要了解钻孔换热器(BHE)在短期和长期的出口流体温度。本文专注于短期和中期分析。在此时间范围内，可采用3D有限元仿真或热阻容量模型(TRCM)。前者精度高但计算耗时长，后者计算快速但由于简化假设而精度有限。本文提出了一种结合TRCM速度和有限元仿真精度的新方法，用于单U型钻孔换热器的短期和中期仿真。该方法使用TRCM估算BHE的热响应，然后通过54组不同配置有限元仿真生成的数据集进行插值修正。模型已实现为C++程序并在博洛尼亚大学开源数据库中提供，能在两秒内提供进出口和平均流体温度、平均BHE表面温度、3D和有效钻孔热阻的时间演化。

研究背景与动机

问题描述

地源热泵系统作为高效的建筑供暖制冷技术，其设计需要准确预测钻孔换热器的热响应特性。关键挑战在于：

  1. 计算精度与效率的矛盾：3D有限元仿真精度高但耗时数小时；TRCM快速但精度有限
  2. 短期和中期仿真需求：现有方法在10⁻²·⁵小时到10³小时时间范围内存在精度或效率问题
  3. 实际工程应用需求：需要既快速又准确的仿真工具用于系统设计和优化

现有方法局限性

  1. 解析解方法：将流体替换为发热固体或表面，无法考虑上下行流体间的热交换效应
  2. 3D数值仿真：精度高但单次仿真需数小时，不适合工程快速设计
  3. 现有TRCM模型：采用水平均匀表面温度(HUST)假设，导致中期热响应低估

研究动机

开发一种兼具TRCM速度和有限元精度的混合仿真方法，满足工程实际需求。

核心贡献

  1. 发现并量化了HUST假设的影响：通过对比有无高导热层的有限元仿真，证明HUST假设会低估T_fm、T_in、T_out、R_b3D和R_beff
  2. 开发了改进的TRCM模型：采用分层热路网络，包含流体、管道、灌浆和地层节点，提高了仿真精度
  3. 建立了修正系数数据库：通过54组不同配置的有限元仿真，构建了时间相关的修正系数数据集
  4. 实现了高效的混合仿真方法：结合TRCM基础计算和插值修正，在2秒内达到有限元仿真精度
  5. 提供了开源仿真工具：C++程序可用于快速准确的BHE短中期仿真

方法详解

任务定义

输入：BHE几何参数(长度L、半径r_b、管间距s等)、热物性参数、流量、热负荷 输出：时间演化的进出口温度(T_in, T_out)、平均流体温度(T_fm)、BHE表面温度(T_b)、3D热阻(R_b3D)、有效热阻(R_beff) 时间范围：10⁻²·⁵小时至10³小时

模型架构

1. 有限元基准模型

  • 使用COMSOL Multiphysics的传热固体界面模拟固体域热传导
  • 采用管流模块模拟流体能量平衡
  • 网格独立性验证确保精度
  • 垂直坐标缩放(因子20)提高网格质量

2. 改进TRCM模型

分层结构：BHE和地层分为m个水平层，每层高度l = L/m

热路网络：每层包含以下节点

  • 流体节点：f1(下行)、f2(上行)
  • 管道节点：p1、p2
  • 灌浆节点：gt1、gt2
  • 钻孔节点：b
  • 地层节点：g1, g2, ..., gn

热阻计算：

R₁^Δ = 2R_b                    (流体到钻孔壁)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (流体间热阻)
R_conv = 1/(2πr_pi h)           (对流热阻)

能量平衡方程： 对每个节点建立瞬态能量平衡方程，形成线性方程组求解。

3. 修正系数方法

修正系数定义：

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

修正过程：

  1. TRCM计算基础热响应
  2. 根据参数插值获得修正系数
  3. 修正R_b3D并重新计算温度场

技术创新点

  1. HUST影响的定量分析：通过引入高导热层强制等温边界，定量分析了HUST假设对仿真精度的影响
  2. 分层热路网络设计：相比传统TRCM，增加了管道和灌浆节点，更精确描述BHE内部传热
  3. 后验修正策略：利用TRCM准确预测T_b的特性，通过修正R_b3D间接修正其他参数
  4. 多参数插值技术：基于钻孔半径、管间距、灌浆导热系数、管外径四个关键参数进行抛物线插值

实验设置

仿真配置

几何参数范围：

  • 钻孔半径r_b：64-76 mm
  • 管间距s：27-47 mm
  • 管外径r_pe：16, 20 mm
  • 灌浆导热系数k_gt：1.0-2.2 W/(m·K)

固定参数：

  • BHE长度：100 m
  • 埋深：1.8 m
  • 地层导热系数：1.8 W/(m·K)
  • 热负荷：50 W/m

数据集构建

  • 54组不同配置的BHE参数组合
  • 每组进行COMSOL和TRCM对比仿真
  • 时间点：56个对数等间距点(10⁻²·⁵至10³小时)
  • 生成修正系数数据库用于插值

验证案例

选择3个典型BHE配置验证：

  • BHE 1：标准配置，用于修正系数计算
  • BHE 2：参数需插值的配置
  • BHE 3：需外推的小半径配置

实验结果

主要结果

HUST假设影响分析

通过对比有无高导热层的仿真结果：

  • T_fm、T_in、T_out被低估，T_b几乎无影响
  • R_b3D和R_beff显著低估
  • 低估程度随时间增加，在中期(>4小时)尤为明显

修正效果验证

BHE 1结果：

  • T_fm均方根误差：0.037°C
  • T_in、T_out均方根误差：<0.038°C
  • R_beff最终时刻误差：0.066%

BHE 2结果(需插值)：

  • T_fm均方根误差：0.051°C
  • R_b3D最终时刻误差：0.13%
  • R_beff最终时刻误差：0.23%

BHE 3结果(需外推)：

  • T_fm均方根误差：0.052°C
  • R_b3D最终时刻误差：0.20%
  • R_beff最终时刻误差：0.40%

计算效率

  • 修正TRCM：2秒内完成
  • COMSOL仿真：数小时
  • 精度提升：达到有限元仿真精度

对比分析

与Pasquier-Marcotte人工神经网络对比：

  • 两种方法都存在中期低估问题
  • 本文方法通过修正系数有效解决了此问题

相关工作

解析解方法

  • Lamarche & Beauchamp：短期解析解，但无法考虑流体间热交换
  • Beier：考虑流体能量平衡的解析模型，需拉普拉斯逆变换

数值仿真方法

  • Li & Zheng：3D固体域+1D流体能量平衡
  • Marcotte & Pasquier：3D固体和流体域完全耦合仿真

TRCM方法

  • Bauer等：标准三角热阻网络，采用HUST假设
  • Pasquier & Marcotte：改进TRCM并发展为准3D模型
  • Ruiz-Calvo等：钻孔-地层模型，避免等温假设但精度有限

机器学习方法

  • Pasquier等：基于15000组TRCM仿真训练的人工神经网络

结论与讨论

主要结论

  1. HUST假设是TRCM精度限制的根本原因，会导致中期热响应系统性低估
  2. 后验修正策略有效：利用TRCM准确预测T_b的优势，通过修正R_b3D实现整体精度提升
  3. 混合方法兼具速度和精度：2秒内达到有限元仿真精度，满足工程应用需求
  4. 插值外推效果良好：即使在参数范围外仍保持较高精度

局限性

  1. 仅适用于单U型BHE：双U型或同轴型需要重新建立修正数据库
  2. 参数范围限制：超出插值范围的参数可能影响精度
  3. 修正系数的物理解释：缺乏修正系数与物理机制的深层关联分析
  4. 长期精度未验证：超过1000小时的长期仿真精度有待验证

未来方向

  1. 扩展到其他类型BHE(双U型、同轴型)
  2. 考虑地下水渗流等复杂传热机制
  3. 与长期g函数方法耦合，实现全时间尺度仿真
  4. 开发自适应修正策略，减少对预构建数据库的依赖

深度评价

优点

  1. 问题识别准确：深入分析了HUST假设的物理机制和影响程度
  2. 方法设计巧妙：充分利用TRCM优势，通过后验修正弥补不足
  3. 实验设计完备：54组配置覆盖主要参数空间，验证充分
  4. 工程价值显著：开源C++程序具有实际应用价值
  5. 精度验证严格：多个验证案例证明方法有效性

不足

  1. 理论分析不足：对修正系数的物理含义缺乏深入解释
  2. 参数敏感性分析有限：未充分分析各参数对修正效果的影响
  3. 适用范围限制：仅针对单U型BHE，通用性有限
  4. 长期验证缺失：超长期仿真精度有待验证

影响力

  1. 学术贡献：为BHE仿真领域提供了新的混合建模思路
  2. 工程价值：显著提升了短中期BHE仿真的效率和精度
  3. 开源贡献：提供的C++程序促进了方法的推广应用
  4. 方法论启发：后验修正策略可推广到其他工程仿真问题

适用场景

  1. 地源热泵系统设计：快速准确评估BHE短中期热性能
  2. 热响应测试分析：辅助现场测试数据解释和参数识别
  3. 系统优化研究：为BHE参数优化提供高效仿真工具
  4. 教学科研：作为BHE传热分析的标准工具

参考文献

论文引用了45篇相关文献，涵盖了BHE仿真的主要方法和最新进展，为本研究提供了坚实的理论基础。关键文献包括：

  • Claesson & Javed (2011)：g函数分析方法
  • Bauer et al. (2011)：TRCM标准方法
  • Pasquier & Marcotte (2014, 2020)：改进TRCM和机器学习方法
  • Zanchini (2024)：BHE仿真方法对比分析

总体评价：这是一篇高质量的应用物理论文，针对地源热泵BHE仿真的实际工程需求，提出了创新的混合建模方法。论文理论分析深入，实验验证充分，具有显著的学术价值和工程应用前景。开源程序的提供进一步增强了其影响力和实用价值。