临界绝热半径
临界绝热半径的理论基础
临界绝热半径是什么
老师，我原以为缠绕绝热材料一定会减少散热，但有时反而会增加散热吗？
确实如此。在圆筒或球的外表面缠绕绝热材料时，虽然导热热阻会增加，但外表面积也会增大。对流热阻与面积成反比，所以在某个半径范围内，对流阻力的减少超过导热阻力的增加，反而导致总体散热量增加。
这与直觉相悖呢。
散热量最大的半径称为临界绝热半径 $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。许多工业管道的外半径远超这个值，所以在实际问题中很少出现，但在电线被覆和医疗管道设计中变得重要。
临界绝热半径的数值计算方法
解析解
临界绝热半径是可以得到漂亮解析解的问题呢。
对。作为一维圆筒坐标稳态热传导可以得到严格解。温度分布为
散热量 $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℃时 $k$ 约为常温的两倍。高温管道的绝热设计中忽视温度依存性会导致设计过程不足。
通过绘制热损失曲线进行视觉确认
要体验临界半径的存在，可以绘制以绝热外径为横轴、热损失Q为纵轴的图表。Q(r)=2πLΔT/[ln(r/ri)/λ + 1/(h·r)]的导数令其为零可求出临界点。1970年代Incropera & DeWitt的传热学教科书（现已更新到第7版）将此图作为标准例题，在全球热工程教育中广泛传播。
临界绝热半径的实务应用
设计应用
临界绝热半径在实务中怎样应用呢？
主要有两种场景。
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}$ 处是否有峰值，前后是否单调变化
- 能量守恒：内表面的入热是否与外表面的散热相等
FEM结果如果不是光滑曲线怎么办？
特别是薄绝热层只有一层要素时精度会下降。径向至少要3个要素以上。
医疗导管的绝热设计
通过血管的加温导管（外径2〜4 mm）是临界半径影响的典型例。医用硅橡胶（λ≈0.2 W/m·K）的被覆如果过厚，体液的对流传热面积增大会导致热损失增加。美国FDA提交的510(k)申请文件中已要求附加Fourier热传导分析作为绝热设计的根据。
临界绝热半径的软件比较
商用工具的参数化分析
用商用工具确认临界绝热半径应该怎样做？
以绝热材料厚度为参数进行扫描分析最直接。
| 工具 | 参数化方法 | 所需时间目安 |
|---|---|---|
| Ansys Mechanical | APDL宏中的DO循环 | 10个工况×数秒 |
| Ansys Workbench | Design Explorer中的参数扫描 | 自动更新网格等 |
| COMSOL | 参数化扫描研究 | GUI中定义绝热厚度变量 |
| Abaqus | Python脚本 + DCAX4要素 | 参数研究自动化 |
Ansys Workbench的参数化分析很方便呢。
Workbench中CAD参数（SpaceClaim中的绝热材料厚度尺寸）可直接作为设计变量。从网格更新到求解全部自动化，结果可用响应曲面可视化。
APDL实现例
最小化的APDL参数化分析代码如下。
```
*DO,t_ins,1,50,1 ! 绝热材厚1〜50mm
/PREP7
ET,1,PLANE55,,,1 ! 轴对称
MP,KXX,1,0.2 ! 绝热材
CYL4,0,0,5,0,5+t_ins,90
ESIZE,0.5
AMESH,ALL
/SOL
D,NODE(5,0,0),,100 ! 内面100℃
SFL,LINE(外面),CONV,10,25
SOLVE
*GET,Q_total,FSUM,,HEAT
*ENDDO
```
以1mm间隔运行50个工况呢。
这程度的问题总共只需几分钟。绘制结果图可一目了然地看出 $r_{cr}$ 的位置。
工具选择指南
| 场景 | 推荐 |
|---|---|
| 简单圆筒的概算 | 手计算（$r_{cr}=k/h$）足够 |
| 考虑温度依存性、辐射 | Ansys Mechanical, COMSOL |
| 3D管道绝热优化 | Workbench + Design Explorer |
| 配管系统整体热损失评估 | Flownex, COMSOL Pipe Flow |
能用手算解决的问题就不必非得用FEM不可。
对。工程判断力就是为问题选择合适工具的能力。
用Wolfram Mathematica即时可视化解析解
Wolfram Mathematica 12（2019年）的Manipulate函数可将rcrit = λ/h的参数用滑块动态改变，同时实时绘制热损失曲线。许多大学的传热工程讲义中现在都发放这样的笔记本。MIT的OCW也从2021年起公开了相关的Mathematica笔记本。
临界绝热半径的先端研究
包含辐射的临界绝热半径
高温环境下辐射的影响也会出现吧？
是的。包含外表面辐射时全热阻的分析变复杂。辐射的线性化传热系数为
以 $h_{\text{eff}} = h_{\text{conv}} + h_r$ 作为全外表面传热系数求临界半径时
因为 $h_r$ 随温度变化，需要反复计算，但相比仅对流情况，$r_{cr}$ 会变小。
加上辐射的话临界半径会变小呢。
根据Stefan-Boltzmann定律 $q = \varepsilon \sigma A T^4$，辐射与面积成正比。表面积增加导致的散热增加效果更大。300℃以上的高温管道中，辐射的贡献与对流相当，不可忽视。
多层绝热的优化
多种绝热材料堆叠时，材料的配置顺序会影响性能。
一般在高温侧放耐热材料（陶瓷纤维），低温侧放低热导率材料（气凝胶）。
多层结构的优化用手算会很复杂吧。
用Excel的Solver或Python的SciPy.optimize就能对付。3层以上的话用COMSOL的Optimization模块自动优化也很高效。
真空绝热的发展
真空绝热板（VIP）的 $k = 0.004$ W/(m K) 极其之小，$r_{cr} = 0.4$ mm。几乎所有形状都没问题。在LNG运输坦克和航天器多层绝热（MLI）中有效热导率可达 $10^{-4}$ W/(m K) 量级。
航天器的绝热设计真是另一个世界呢。
真空环境中没有对流，辐射遮蔽是主要绝热机制。地面的临界绝热半径概念不适用，但遮蔽片数优化的问题相类似。
球壳绝热没有临界半径
与圆筒绝热不同，球壳绝热中增大外半径时热损失单调递减（无极值）。这是因为球的表面积随半径平方增长，导热阻力也按平方项增长，两者保持平衡，不会产生极值点。Cengel & Ghajar的传热学教科书（2020年第6版）用比较图表清晰地说明了球壳与圆筒的这种区别。
临界绝热半径的故障处理
常见问题与对策
临界绝热半径方面在实务中会遇到什么困难？
整理一下典型的问题。
1. 增加绝热材料后温度反而上升
原因：原管径在 $r_{cr}$ 以下，增加绝热材料反而增加了散热。
对策：计算 $r_{cr} = k/h$ 与现有管径比较。若 $r_i < r_{cr}$，选择热导率更低的绝热材料，或在外表面加风挡以降低 $h$。
2. 绝热材料结露
低温管道的绝热中结露会成为问题吧。
绝热材料内部温度降到露点以下就会结露。纤维型绝热材料吸水后 $k$ 急剧增加，结露进一步加剧，形成恶性循环。
对策：在绝热材料外表面（温度较高一侧）加防潮层（PE/AL膜）。低温配管绝对必须。
3. 分析结果与理论值不符
FEM结果与理论值的不一致检查清单。
| 检查项目 | 常见原因 |
|---|---|
| 轴对称设置 | 2D要素未设置轴对称 |
| 对流面 | 内表面错误设置对流条件 |
| 单位制 | mm单位下 $k$ 的值未换算 |
| 网格 | 薄绝热层只有一层要素 |
| 端面条件 | 端面是否绝热（默认） |
轴对称设置忘记很常见呢。
2D平面默认计算的话面积变化的效果无法包含，无法重现临界绝热半径现象。Ansys中PLANE55须设KEYOPT(3)=1，Abaqus中DCAX4要素要明确选择轴对称。
细电线过热故障要小心
直径3mm以下的电线缠绕聚氯乙烯绝缘（λ≈0.17 W/m·K）时，增加被覆厚度可能导致发热体（电线）温度升高。UL规范中规定了不同电流容量下的最优被覆厚度，但选择规范外的组合时会出现过热故障。UL White Paper中记录了此类实际故障案例。
更详细
错误