1. COMSOL 中一维稳态热传导解析验证
问题定义与解析基准
考虑一根长度 L = 0.2 m 的均质导热杆，导热系数 k = 45 W/(m·K)，无内热源。左端 x = 0 为第一类边界（定温 T₀ = 300 K），右端 x = L 为第二类边界（定热流 q = 5000 W/m²，沿 +x 方向流出）。一维稳态导热控制方程为：
d/dx( k · dT/dx ) = 0 → T(x) = T₀ + (q/k)·x代入参数可得解析解：T(x) = 300 + (5000/45)·x，右端温度 T(L) = 300 + 5000×0.2/45 ≈ 322.22 K。该线性温度分布是后续数值解对比的唯一基准。
COMSOL 建模关键设置
在 COMSOL 中新建「传热传质 → 固体传热（ht）」稳态研究，空间维度选一维。几何用「区间」从 0 到 0.2 m。材料参数直接输入 k = 45，密度与热容在稳态下不参与计算，可保留默认值。边界条件设置如下：
- 左端点（x=0）：选「温度」边界，
T₀ = 300 K。 - 右端点（x=L）：选「热通量」边界，
q₀ = 5000 W/m²，方向为外法向（即 +x 流出）。 - 初始值仅作迭代起点，稳态结果不依赖它，可设
T_init = 300 K。
网格收敛性验证方法
一维线性温度场的数值解在理论上对任意网格都应精确，但 COMSOL 默认的二次单元在单元数极少时仍可能因边界热流积分精度产生微小偏差。验证时依次设置单元数为 2、4、8、16、32，记录右端温度 T(L)，并与解析值 322.22 K 比较。判断标准：相邻网格规模间 T(L) 变化小于 0.01 K，且与解析解偏差小于 0.05 K，即可认为网格收敛。
数值解与解析解对比
在 COMSOL 中计算后，用「点图」提取各节点温度，或用「一维绘图组」绘制温度曲线。下表为不同单元数下右端温度 T(L) 的示例结果（教学示例数据，非实测）：
| 单元数 | 节点数 | T(L) 数值解 (K) | 与解析解偏差 (K) | 相邻变化 (K) |
|---|---|---|---|---|
| 2 | 3 | 322.18 | 0.04 | — |
| 4 | 5 | 322.21 | 0.01 | 0.03 |
| 8 | 9 | 322.22 | 0.00 | 0.01 |
| 16 | 17 | 322.22 | 0.00 | 0.00 |
| 32 | 33 | 322.22 | 0.00 | 0.00 |
从 8 单元起，数值解与解析解偏差已小于 0.01 K，相邻网格变化为 0，说明网格已收敛。若继续加密，结果不再变化，可停止。
边界条件正确性检查
除温度对比外，还应检查热流是否守恒。在 COMSOL 中计算左端热流 q_left = -k·dT/dx，理论上应等于右端流出热流 5000 W/m²。若左端热流与右端不相等，说明边界条件设置或材料参数有误。常见错误包括：热通量方向设为流入、导热系数单位错误、左端温度未固定等。
可复现操作清单
- 新建一维稳态「固体传热」模型，几何区间 [0, 0.2] m。
- 材料导热系数设为 45 W/(m·K)，其余参数保持默认。
- 左端点设「温度」300 K；右端点设「热通量」5000 W/m²（外法向流出）。
- 网格序列设为 2、4、8、16、32 单元，逐次计算。
- 每次计算后提取右端温度，填入对比表，观察偏差与相邻变化。
- 确认收敛后，检查左端热流是否等于 5000 W/m²，完成验证。