NAFEMS T4：对流边界条件热传导
基准问题的规格与理论背景
NAFEMS基准问题是什么
NAFEMS（发源于英国的分析工程师协会）发布的标准基准问题集提供带参考解的小规模验证问题，是求解器验证的世界通用语言。T4是其传热篇的代表问题，设计成用1道题就能确认“对流边界条件（第3类、即Robin边界）是否被正确实现与设置”。在引入新求解器、版本升级以及传热分析教学时，它都是应当最先求解的经典题目。
T4的问题设定
稳态2D热传导，几何为宽0.6 m×高1.0 m的矩形板。
| 项目 | 数值 |
|---|---|
| 导热系数 | \( k = 52 \) W/(m·K) |
| 下边（y=0） | 给定温度 \( T = 100 \)°C |
| 左边（x=0） | 绝热（对称） |
| 右边与上边 | 对流：\( h = 750 \) W/(m²·K)，环境温度 \( T_\infty = 0 \)°C |
| 评价点E（x=0.6, y=0.2） | 参考解 \( T_E = 18.3 \)°C |
控制方程为拉普拉斯方程 \( \nabla^2 T = 0 \)，对流边上的边界条件为 \( -k\,\partial T/\partial n = h\,(T - T_\infty) \)。
为什么说它是“好的验证问题”
不就是一块矩形板的热传导吗？这么简单的问题究竟能验证出什么？
它比看上去刁钻得多。先看毕渥数——\( Bi = hL/k = 750\times0.6/52 \approx 8.7 \)。对流强于导热，温度场会在边界附近出现陡峭的梯度。而且评价点E靠近给定温度边（100°C）与对流边相交的右下角，正是梯度最陡的位置，网格一粗就轻易偏出几摄氏度。再加上对流边界是以“表面单元上的积分”实现的，h给错、面的朝向弄反、环境温度设错，都会立刻反映到数值上。Dirichlet、绝热、Robin这3种边界，用1道题全部走一遍——作为设置错误的检测器，它设计得实在漂亮。
实施步骤与收敛确认
标准的实施步骤
- 单元选择——推荐2阶四边形单元（8节点）。1阶单元也能求解，但收敛较慢
- 网格序列——例如6×10、12×20、24×40这3个水平（细分比2）。选择让评价点E落在节点上的划分，插值误差的争论就自然消失
- 边界设置——下边给100°C，右边与上边给对流（h=750、T∞=0°C），左边什么都不设（自然成为绝热）
- 评价——把各水平的 \( T_E \) 列成表，确认其向参考值18.3°C收敛
合格判定的标准
用2阶单元的细网格时 \( T_E \) 收敛到18.2～18.3°C。实务上的合格标准是单调收敛到参考值±0.5%（±0.1°C）以内。只解1个水平就以“算出18.6°C，大致吻合”收尾，作为验证是不充分的：只有把网格收敛的趋势也展示出来，才算真正解过T4。若再做到收敛阶的评估（GCI），它就能直接成为公司内部的验证模板。
要有意识地看待对流边界的离散化
在有限元法中，对流边界条件是以边界单元上的积分来实现的：\( \int_\Gamma h\,N_i N_j\,dS \) 进入刚度矩阵，\( \int_\Gamma h\,T_\infty N_i\,dS \) 进入载荷向量。也就是说h是“附着在面上”的量，误加到节点或体上、选错面、把单元的正反面弄反，都会原样出现在温度场里。若T4能重现参考值，就证明这条边界实现链路整体是正确的。
实务中的活用方法
作为验收试验与回归测试的用法
T4的实务价值不在于“解一次就结束”，而在于反复使用的基准器：①新求解器、新版本的验收试验（升级后边界处理有没有改变）；②分析环境迁移时的等效性确认；③新人培训的第1个课题（用身体记住3种边界的设置方法）；④输入模板的合理性确认。把输入文件与期望值（18.3°C±0.1）脚本化，此后每次环境更新只要几分钟就能完成回归确认。
单位制与温度基准的陷阱
T4以摄氏度定义，但由于是不含辐射的线性问题，改用绝对温度也会得到同样的答案（取T∞=0°C=273.15 K与100°C=373.15 K，则 \( T_E = 291.45 \) K = 18.3°C）。不过一旦与求解器的默认温度单位混用就会立刻出事——“给T∞填了0，但单位其实是K（＝−273°C）”是经典错误，结果大幅偏向负值时首先应当怀疑它。在非SI单位制（mm制）下，h的单位换算（W/m²K → mW/mm²K）同样需要注意。
扩展的方法——以T4为起点的分级验证
T4通过之后，可以用同一个模型把验证范围逐级扩大：①加入温度相关物性，确认非线性热传导的迭代收敛；②改为瞬态，验证时间积分（从初始温度0°C开始的过渡过程）；③让对流h随温度变化，验证边界的非线性；④把固体分成2种材料，验证界面连续性。“从有参考解的问题出发只走1步扩展”，是在复杂模型中把程序缺陷与物理现象区分开来的最短路径。
各工具的设置方法
对流边界条件的设置名对照
| 工具 | 对流边界的设置 | 备注 |
|---|---|---|
| Abaqus | *FILM（在表面上给出膜系数h与Sink温度） | 施加在面（surface）定义上。单位制在整个输入中按无量纲方式统一管理 |
| Ansys Mechanical | Convection载荷（Film Coefficient＋Ambient Temperature） | MAPDL中为SF,,CONV。注意确认边／面的选择 |
| Nastran（热） | CONV＋PCONV（对流边界），温度给定用相当于SPC的卡片 | 注意热分析的边界卡片体系与结构分析是分开的 |
| COMSOL | 热通量边界中的“对流热通量”（h与T_ext） | 2D的面外厚度设置不影响结果（线性稳态），但建议显式指定 |
| OpenFOAM | externalWallHeatFluxTemperature（mode: coefficient） | 固体导热用laplacianFoam系列＋边界类型来搭建 |
| 自研代码与教学用FEM | 自行实现边界积分项 | T4正是检验该实现的最佳验证问题 |
建模细节中必须对齐的事项
在工具之间比对结果时，要统一以下几点：①评价点数值的取法（节点值还是单元内插值）；②2阶单元中间节点的处理；③2D列式的类型（有没有在纯2D、单位厚度平面、轴对称之间选错）。像T4这样简单的问题若在工具间仍残留0.1°C以上的差异，原因几乎全在这些后处理与列式的差别，而不在物理。
基准问题文化的现状
标准基准问题的体系与T4的位置
NAFEMS基准问题集分为传热（T系列）、线弹性（LE系列）、自由振动（FV系列）等体系，全部按照“有参考解、并针对特定功能做定点验证”的思想设计。T4针对对流边界，LE5针对约束扭转的翘曲，LE11针对热应力，1题1功能的干脆取舍正是其特征。用ASME V&V规范的说法，它们属于代码验证（code verification）的一部分，位于依据实验比对的确认（validation）之前。
走向自动化的验证套件
在现代的求解器开发与运维中，把NAFEMS型基准问题接入CI（持续集成）、每次代码变更都自动执行，已经成为标准做法。用户一侧也是如此：越来越多的组织会在分析环境（求解器、操作系统、库）更新时自动运行内部基准问题集，并报告与参考值的差异。T4这类轻量问题最适合这个用途，它构成了用机械化手段保证“环境变了答案也不变”的安全网。
参考解的高精度化
经典基准问题的参考值是用当年的计算资源做出来的，如今的高精度计算（hp-FEM、级数解、外推）已经把更多位数确定下来。T4的18.3°C在实用上完全够用，但若想做更严格的收敛阶讨论，实务上的做法是同时准备一个“工作参考值”——把自家超密网格解做Richardson外推得到的数值。
故障排查
“算不出18.3°C”时的诊断表
| 症状 | 可能原因 | 对策 |
|---|---|---|
| 高或低几摄氏度（网格加密后逐渐接近） | 单纯的离散化误差 | 用2阶单元＋细分确认收敛。属于正常行为 |
| 数值大幅偏负、非物理 | 温度单位混用（T∞=0的单位被当作K） | 检查求解器的温度单位设置 |
| 整体温度偏高 | 对流边漏设（仍为绝热）、h的数量级写错 | 把边界的分配可视化后确认。核对h=750的单位 |
| 温度场左右不对称 | 绝热边与对流边弄反 | 重新确认边的对应关系（左＝绝热，右与上＝对流） |
| 评价点数值在工具间略有差异 | 取的是节点值还是插值，以及平均方式不同 | 让评价点与节点重合。统一提取方法 |
| 收敛出现振荡 | 自适应网格的非系统性细分 | 改为结构化网格的系统性细分（6×10→12×20→24×40） |
基准问题的真正用法
顺利算出18.3°C了！这样是不是就可以说“我们的传热分析已通过验证”？
能说的只是“稳态线性热传导与对流边界的实现和设置是正确的”——即便如此也已经是一大步了。不过实际项目里还会出现辐射、温度相关物性、接触热阻、瞬态等因素，因此正确的做法是按功能逐项累积对应的基准问题。我推荐把“要用到的功能清单”与“已验证的基准问题清单”做成对照表。表中的空白格，就是你的分析环境里尚未验证的功能——也就是风险之所在。T4正是填上这张表第1行的那道题。
帮助
更多
错误