臨界断熱半径
理论与物理
什么是临界隔热半径
老师，我以为只要包上隔热材料就一定能减少散热，但听说有时反而会增加，这是真的吗？
是真的。在圆柱或球的外面包覆隔热材料时，传导热阻会增加，但外表面积也会增大。由于对流热阻与面积成反比，在达到某个半径之前，对流热阻的减少会超过传导热阻的增加，导致整体散热量反而增加。
这有悖直觉呢。
散热量达到最大时的半径称为临界隔热半径 $r_{cr}$。对于圆柱，总热阻为
将其对 $r$ 求导并令其为零，可得
对于球体，则为 $r_{cr} = 2k/h$。
$k$ 越小、$h$ 越大，临界半径就越小呢。
是的。对于自然对流（$h \approx 5$ W/(m$^2$ K)）和环氧树脂隔热材料（$k \approx 0.2$ W/(m K)），$r_{cr} = 0.04$ m = 40 mm。如果原始管道外径小于40mm，那么包覆隔热材料后，会存在一个效果适得其反的区域。
物理解释
将总热阻分解开来考虑会更容易理解。
| 半径 $r$ | 传导热阻 | 对流热阻 | 总热阻 | 散热量 |
|---|---|---|---|---|
| $r < r_{cr}$ | 增加（小） | 减少（大） | 减少 | 增大 |
| $r = r_{cr}$ | — | — | 最小 | 最大 |
| $r > r_{cr}$ | 增加（大） | 减少（小） | 增加 | 减少 |
原来是对流热阻的减少能否胜过传导热阻的增加的较量啊。
正是如此。在实际应用中，电线包覆和管道隔热会出现这个问题。设计时需确认 $r_i > r_{cr}$，然后再决定隔热材料厚度。
rcrit = λ/h 的直观含义
圆柱隔热的临界半径 rcrit = λins/ho 是“隔热材料热传导阻力的增加速度”与“外表面热对流阻力增加速度”达到平衡的点。典型的室外条件（ho = 10 W/m²·K）和玻璃棉（λ=0.04 W/m·K）下，rcrit = 4 mm。许多工业管道的外半径远大于此，因此在实践中不会出现此问题，但在电线包覆或医用导管设计中则很重要。
各项的物理意义
- 蓄热项 $\rho c_p \partial T/\partial t$：单位体积的热能储存率。【日常例子】铁锅不易加热也不易冷却，而铝锅则易加热易冷却——这是密度 $\rho$ 和比热 $c_p$ 的乘积（热容）不同所致。热容大的物体温度变化缓慢。水的比热非常大（4,186 J/(kg·K)），因此沿海地区气温比内陆更稳定。在非稳态分析中，此项决定了温度随时间的变化速率。
- 热传导项 $\nabla \cdot (k \nabla T)$：基于傅里叶定律的热传导。与温度梯度成比例的热流。【日常例子】将金属勺放入热锅中，勺柄也会变热——因为金属的热导率 $k$ 高，热量能迅速从高温侧传到低温侧。木勺不会变热是因为其 $k$ 值小。隔热材料（如玻璃棉）的 $k$ 值极低，即使存在温度梯度，热量也难以传递。这是将“有温差的地方就有热流”这一自然趋势公式化的结果。
- 对流项 $\rho c_p \mathbf{u} \cdot \nabla T$：伴随流体运动的热量输送。【日常例子】吹风扇感到凉爽，是因为风（流体流动）带走了体表附近的热空气，并供应了新鲜的冷空气——这就是强制对流。暖气使房间天花板附近变暖，是因为受热空气因浮力上升的自然对流。PC的CPU散热风扇也是通过强制对流来散热。对流是比热传导效率高得多的热量输送方式。
- 热源项 $Q$：内部发热（焦耳热、化学反应热、辐射吸收等）。单位：W/m³。【日常例子】微波炉通过食品内部的微波吸收（体积发热）来加热。电热毯的加热线通过焦耳发热（$Q = I^2 R / V$）来产生热量。锂离子电池充放电时的发热、刹车片的摩擦热在分析中也作为热源考虑。与外部向“表面”提供热量的边界条件不同，热源项表示“内部”的能量生成。
假设条件与适用范围
数值解法与实现
解析解法
临界隔热半径是一个能得到漂亮解析解的问题呢。
是的。作为一维圆柱坐标下的稳态热传导问题，可以得到精确解。温度分布为
散热量 $q$ 为
将 $r_o$ 作为变量，令 $dq/dr_o = 0$ 即可导出临界半径。
通过数值分析进行验证
用FEM也能得到相同结果吗？
当然可以。实际上，与理论解的比较常被用作求解器的基准验证。在Ansys Mechanical中创建SOLID70的圆柱网格，在内表面施加温度约束，在外表面设置对流条件即可。
| 参数 | 值 |
|---|---|
| 内径 $r_i$ | 5 mm |
| 隔热材料 $k$ | 0.2 W/(m K) |
| 外部 $h$ | 10 W/(m$^2$ K) |
| $r_{cr}$ 理论值 | 20 mm |
将隔热层厚度从5mm变化到50mm进行参数化分析，可以确认散热量在 $r_o = 20$ mm 处达到最大。与理论值的误差在0.1%以下。
这真是个绝佳的验证问题。
用APDL宏编写DO循环，就能自动运行所有工况。在Abaqus中同样可以用Python脚本，以 $r_o$ 为参数进行研究。
包含温度依赖性的情况
当隔热材料的 $k$ 依赖于温度时，$r_{cr} = k(T)/h$ 中的 $k$ 应在哪个温度下评估就成了问题。需要使用Newton-Raphson迭代法求自洽解。
现实中的隔热材料，温度升高时 $k$ 通常也会变大，对吧？
玻璃棉在200℃时约为常温的2倍。高温管道的隔热设计如果忽略温度依赖性，可能会导致设计不足。
通过绘制热损失曲线进行视觉确认
要直观感受临界半径的存在，可以以隔热层外径为横轴，热损失Q为纵轴绘图。对 Q(r)=2πLΔT/[ln(r/ri)/λ + 1/(h·r)] 求导，令 dQ/dr=0 即可求得临界点。1970年代的Inciropera & DeWitt热传导教科书（现已修订至第7版）将此图作为标准例题，推广至全球的热工学教育中。
线性单元 vs 二次单元
在热传导分析中，线性单元通常足以获得足够的精度。在温度梯度陡峭的区域（如热冲击等），推荐使用二次单元。
热流密度评估
根据单元内的温度梯度计算得出。有时需要像节点应力那样进行平滑处理。
对流-扩散问题
当佩克莱特数较高（对流主导）时，需要迎风稳定化（如SUPG等）。纯热传导问题则不需要。
非稳态分析的时间步长
时间步长应相对于热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$：热扩散率）足够小。对于急剧的温度变化，自动时间步长控制是有效的。
非线性收敛
由温度依赖物性值引起的非线性通常比较温和，皮卡德迭代法（直接替换法）通常就足够了。对于辐射的强非线性，推荐使用牛顿法。
稳态分析判定
当所有节点的温度变化低于阈值（例如 $|\Delta T| / T_{max} < 10^{-5}$）时，判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算速度快，但时间步长过大时会不稳定（漏掉风暴）。隐式法是“也考虑未来状态的预测”——即使时间步长较大也能保持稳定，但每个时间步都需要解方程，比较耗时。对于没有急剧温度变化的问题，使用隐式法并采用较大的时间步长更高效。
实践指南
设计应用
临界隔热半径在实际工作中如何应用呢？
主要有两种场景。
1. 隔热设计：决定管道或风管的隔热层厚度时，先确认 $r_i > r_{cr}$，再进行厚度优化
2. 散热设计：决定电线包覆厚度时，如果 $r_i < r_{cr}$，则增加包覆厚度反而能改善散热
反向应用于散热设计，这很有趣。
电线包覆是典型例子。对于AWG24铜线（外径0.56mm）包覆PVC（$k = 0.16$ W/(m K)）的情况，自然对流 $h = 10$ W/(m$^2$ K) 下，$r_{cr} = 16$ mm。在包覆外径达到32mm之前，增加包覆厚度反而能改善散热。
典型隔热材料的物性
| 隔热材料 | $k$ [W/(m K)] | 使用温度范围 | $r_{cr}$（$h$=10） |
|---|---|---|---|
| 玻璃棉 | 0.04 | 〜450℃ | 4 mm |
| 岩棉 | 0.04 | 〜700℃ | 4 mm |
| 聚氨酯泡沫 | 0.02 | 〜100℃ | 2 mm |
| 二氧化硅气凝胶 | 0.015 | 〜650℃ | 1.5 mm |
| 陶瓷纤维 | 0.08 | 〜1200℃ | 8 mm |
性能越好的隔热材料 $r_{cr}$ 越小，所以即使是细管道也能放心使用呢。
气凝胶的 $r_{cr} = 1.5$ mm，所以对于外径3mm以上的管道，几乎不成问题。虽然成本高，但在航天、LNG领域有实际应用。
结果验证
临界隔热半径分析结果的验证要点如下。
- 与理论值比较：是否与 $r_{cr} = k/h$（圆柱）或 $2k/h$（球体）一致
- 散热量曲线：在 $r = r_{cr}$ 处是否出现峰值，前后是否单调变化
- 能量平衡：内表面的入热与外表面的散热是否一致
なった
詳しく
報告