地暖系统模拟 — 埋设管道热传导分析与温度分布预测

简单来说，我们要用FEM求解埋设管道的热传导，预测地板表面的温度分布。别看"只是铺管道"，实际上管道的布置模式会完全改变舒适性。

比如管道间距200mm和300mm，地板表面的温度不均匀度（最高温和最低温的差值 $\Delta T$）完全不同。200mm间距时 $\Delta T \approx 1\text{--}2°\text{C}$，而300mm间距时 $\Delta T \approx 3\text{--}5°\text{C}$。脚底能明显感觉到"这里冷，那里热"。

日本建筑规范推荐地板表面温度29°C以下。WHO指南也说居室地板应该在19~29°C范围内。太高容易有低温烫伤风险，太低没有暖房效果。为了命中这个"舒服的范围"，需要对管道间距、埋设深度、保温层厚度进行综合最优化，这就是CAE的用武之地。

好问题。地暖的热传递分为三个机制：

完全一样。空调主要靠对流暖空气，所以容易出现"头热脚冷"。地暖反过来实现"头凉脚暖"。Fanger舒适度研究（1972年）表明，脚部24°C、头部19°C是人体最舒适的状态。所以CAE里辐射的模型化特别关键。

基础是三维非定常热传导方程。地板楼板（混凝土+砂浆+面层）内的温度分布这样求解：

管道内的热水通常用一维流动方程处理。如果把管道的弯曲和分支全部用3D CFD求解，计算代价太大，所以实务中通常采用"1D管道+3D固体传导"耦合模型：

这里 $T_w$ 是水温，$v$ 是流速，$s$ 是沿管道长度的坐标，$A_p$ 是管道截面积，$d_i$ 是管道内径，$q_p$ 是管道壁面单位长度的放热量。水流过管道时逐渐冷却，这样可以追踪整个过程。

正是。热水式地暖的入口温度通常是40~60°C，到出口时会下降5~10°C。这成为管道配置设计的重要参数。比如螺旋形布局，入口和出口管道交替分布，这样能很好地平均化温度不均。

可以把热从热水到室内空气的总热阻理解为"串联电路"。就像电的欧姆定律一样：

命中了。实际上面层材料的选择是地暖设计中最容易被忽视但最关键的点。实木地板（$k \approx 0.12$ W/(m·K)）和瓷砖（$k \approx 1.0$ W/(m·K)）的热导率相差约8倍。瓷砖能更有效地把管道的热传到表面，所以即使热水温度相同，瓷砖地板的表面温度也会更高。欧洲浴室常见瓷砖配地暖的组合就是这个原因。

ISO 7730定义的PMV（预测平均意见值）指标是标准。从温度、湿度、气流速度、辐射温度、穿着量、代谢率这6个要素，量化为"过热（+3）～过冷（-3）"的感觉尺度。把地暖模拟结果（地板表面温度分布、室内气流场）输入PMV计算，就能空间映射出哪里舒适、哪里不舒适。

特别是地暖中，地板表面温度的不均一性（$\Delta T_{floor}$）受到ASHRAE 55的限制。坐着的人脚底感受到3°C以上的温差就会不舒服。通过CAE事先检查这一点，就能在施工前评估管道布置的好坏。

用有限元法（FEM）把空间离散化。地板楼板分成许多小单元，在每个单元内近似求解温度。通过弱形式（变分形式）和Galerkin法，最终得到这样的联立方程组：

这里 $[C]$ 是热容量矩阵，$[K]$ 是热传导矩阵，$\{F\}$ 是外部热输入向量。和结构分析的 $[K]\{u\} = \{F\}$ 形式很像吧？热分析里温度是未知数，只是用热传导矩阵代替刚度矩阵而已。

地暖分析的标准是8节点六面体单元（线性）。地板楼板是板状，六面体单元最合适。管道周围细分，远离管道的地方粗分——非均一网格是基本做法。

典型流程是这样的：

地暖网格的三大法则：

20叠米（约33m²）的客厅一间，2D对称截面模型的话数千单元足够。完整3D加管道路径的话，10万~50万单元的范围。网格收敛性一定要检查——粗网格→中→细，结果在2%以内收敛就OK。

地暖要设的边界条件整理一下：

管道内是强制对流。用Dittus-Boelter关联式算：$\text{Nu} = 0.023 \, \text{Re}^{0.8} \, \text{Pr}^{0.4}$。实际上管道内壁热阻占总阻的比例很小，精度影响不大。最敏感的（对结果影响最大）是面层热导率和管道埋设深度。这些错了温度分布完全变样。

参数化分析探索最优解。设计变量和目标函数这样设：

完全对。典型结果汇总：

保温层方面，30mm的XPS（挤塑聚苯乙烯泡沫，$k \approx 0.035$ W/(m·K)）能切断80%以上的向下热损失。50mm能达到90%以上，但成本和地板厚增加。一般住宅30~40mm最优。现在的高断热住宅如果基础保温做得好，20mm也够。

现场常见的三个失败：

正是。面热源模型可以用来估算"总散热量"，但没法评估"管道间的温度波纹"。目的决定模型细度——这是分析工程师最重要的原则。

主要选项按地暖功能对比：

各有所长。COMSOL的1D-3D耦合管道模型最顺畅，地板楼板内温度分布分析最理想。Fluent在完全耦合的室内空气CFD+固体传导（共轭热传递）上最强。学术论文中地暖的CFD分析用Fluent最多。TRNSYS特化于年度能耗评估，虽然没有FEM式的温度分布，但对热泵+管道+室温控制的系统整体优化不可或缺。

行得通。选项有：

按目的决策：

最热的是数字孪生。在建筑里装几十个IoT温度传感器，把实时计测数据和CAE模型同步。按室外温度、日射、在室人数变化，在线优化热水温度——"预测制御"变成现实。

对。大型办公楼商业设施，只把暖房启动时间提前30分钟，能源消费就能削减10~15%。混凝土楼板的蓄热效应（热的惯性）用得更巧妙，晚上便宜电价时蓄热，白天释热，成本优化的可能很大。

完全可以。用FEM生成20~50个学习数据，高斯过程回归（Kriging）或神经网络做代理模型。这样数千个参数组合秒级评估。管道间距、深度、保温厚、热水温的四变量优化现实时间内完成。

最近PINN（物理约束神经网络）也被用于地暖控制优化。物理法则（热传导方程）编入神经网络的损失函数，学习数据少也能给出物理上可信的预测。

PCM是石蜡等材料，融点（比如23°C）附近吸收大量潜热。地板材料里混进PCM微胶囊，能蓄热容量增几倍。模拟用等效比热法（焓法）标准做法：

实现上把delta函数用高斯分布近似，融点 $T_m$ 前后2~3°C范围内比热增大。COMSOL有内置的相变材料模型，设定很简单。PCM入地暖，关机后还能维持23°C的地板温2~4小时——深夜便宜电蓄热、白天用热的"峰值转移"有很大潜力。

先检查这些：

常见原因：