临界绝热厚
临界绝热厚的理论基础
临界绝热厚的概念
老师，临界绝热厚和临界绝热半径是同一回事吗？
本质上是同一现象。临界绝热厚 $t_{cr}$ 是圆筒外面包覆绝热材料时放热量最大的厚度，用 $t_{cr} = r_{cr} - r_i = k/h - r_i$ 表示。
平板没有临界绝热厚吧？
正确。平板由于绝热材厚度增加时外表面积不变，全热阻单调增加。只有圆筒和球等面积变化的几何形状才存在这一现象。
圆筒的全热阻
以绝热材厚度 $t$ 为变量的全热阻为
从 $dR/dt = 0$ 得到 $r_i + t = k/h$。因此 $t_{cr} = k/h - r_i$。当 $r_i > k/h$ 时，$t_{cr} < 0$，即不存在临界绝热厚问题。
球的情况
在球坐标中，全热阻为
临界半径为 $r_{cr} = 2k/h$，是圆筒的两倍。
球的临界半径比圆筒大是因为面积增加方式不同吗？
是的。球的表面积为 $4\pi r^2$，与半径的平方成正比；圆筒的侧面积 $2\pi r L$ 与半径的一次方成正比。球的面积增加效应更强，因此临界半径也更大。
设计判断
| 形状 | $r_{cr}$ | 实际影响 |
|---|---|---|
| 平板 | 无 | 增加绝热材总是有效的 |
| 圆筒 | $k/h$ | 小径配管、电线需注意 |
| 球 | $2k/h$ | 球形容器绝热设计需注意 |
增加绝热材反而增加热损失
平板绝热中增加厚度会降低热损失，但圆筒配管包覆绝热材料时，外表面积的增大会加剧对流，导致在"临界半径" rcrit = λ/h 超过之前热损失反而增加这一悖论。这一现象最初由1896年Rayleigh在圆柱体热辐射论文中数学证明。
临界绝热厚的数值计算方法
计算方法选择
临界绝热厚应该用什么方法计算比较恰当？
根据问题复杂度选择合适方法。
| 条件 | 手工计算 | 数值分析 |
|---|---|---|
| 简单圆筒、恒定 $k$、恒定 $h$ | 足够 | 无需 |
| 温度相关 $k(T)$ | 近似可行 | 推荐 |
| 包含辐射 | 线性化近似 | 推荐 |
| 非圆形截面 | 不可行 | 必需 |
| 多层绝热 | 复杂但可行 | 推荐 |
放热量的精确计算
单位长度的放热量作为绝热厚度的函数为
在Excel中制表可以轻松找到临界点。无需求导就能数值地找到最大值。
Excel就足够了吗？
基本情况用Excel的目标求解或规划求解就足够。但包含温度依赖性时，应该用Python的SciPy.optimize或FEM更保险。
FEM验证
Ansys Mechanical的验证步骤如下。
1. 建立2D轴对称模型（PLANE55, KEYOPT(3)=1）
2. 对圆环截面建模，内径 $r_i$
3. 使用APDL DO循环改变外径作为参数
4. 每个工况提取 Total Heat Flow
5. 与理论值 $r_{cr} = k/h$ 比较
薄绝热层的网格如何处理？
径向至少3个单元。即使绝热厚1mm，单元大小0.3mm也足够。单元长宽比以5以下为目标。
最优绝热厚的经济计算
最优绝热厚是通过"绝热材投资增量=热损失节能成本削减"经济优化问题求解。日本机械学会1980年代制定的配管绝热设计标准（JSME S010）规定了以LCC（生命周期成本）最小化为条件的厚度算出步骤，至今仍作为化工装置设计的标准被采用。
临界绝热厚的实务应用
实务应用场景
请告诉我临界绝热厚在实际中成为问题的情况。
我列举几个代表性案例。
电线的被覆设计
AWG24铜线（外径0.56mm）包覆PVC被覆（$k = 0.16$ W/(m K)），自然对流 $h = 10$ W/(m$^2$ K) 时，$r_{cr} = 16$ mm。被覆外径达32mm之前，厚度增加能改善放热。
增加被覆反而变冷，直觉违反啊。
电力电缆的允许电流计算（IEC 60287）考虑了被覆的放热改善效果。但实际被覆厚度由机械保护和绝缘决定。
小径冷媒管的保温
冷媒配管（外径6.35mm）包覆橡胶系绝热材（$k = 0.04$）时，$r_{cr} = 4$ mm。配管半径3.2mm与之接近，所以薄绝热效果不佳。
| 绝热厚 [mm] | 放热量比 | 表面温度 |
|---|---|---|
| 0 | 1.00 | 低 |
| 1 | 1.01 | 略微下降 |
| 5 | 0.92 | 开始上升 |
| 15 | 0.70 | 绝热效果充分 |
最少要5mm以上才有效啊。
实务中一般采用10mm以上绝热厚，所以通常不成问题。但"薄绝热就够"这种判断很危险。
经济绝热厚
实务中通过比较绝热材投资成本和省能成本来确定最优厚度。JIS A 9501"保温保冷工程施工标准"规定了经济厚的算出方法。
这与临界绝热厚是不同的优化问题啊。
临界绝热厚是物理限制，经济绝热厚是成本优化。需要同时理解两者才能做好设计判断。
工业装置的保温标准
石油化学装置中常见的蒸汽配管（150°C、100mm径）用玻璃棉保温材通常设计为50~75mm厚。JIS A 9504（2018年版）规定了使用温度对应的最小厚度，低于此值会导致表面温度过高，违反劳动安全卫生法"可接触表面≤60℃"的要求。
临界绝热厚的软件比较
商用工具验证
请告诉各工具的实现方法。
用基准问题在各工具中验证的步骤如下。
Ansys Mechanical
```
/PREP7
ET,1,PLANE55,,,1 ! 轴对称
MP,KXX,1,0.2 ! k=0.2 W/(mK)
CYL4,0,0,5,0,25,90 ! ri=5mm, ro=25mm
ESIZE,0.5
AMESH,ALL
/SOL
D,NODE(5,0,0),,100 ! 内表面T=100
SFL,LINE(外表面),CONV,10,25 ! h=10, Tinf=25
SOLVE
```
COMSOL
1. 2D轴对称 > 固体热传导
2. 几何：矩形（宽=绝热厚，高=L）
3. 内边界：温度 = $T_i$
4. 外边界：对流热Flux
5. 参数扫描：外半径5~50mm
Abaqus
用DCAX4单元建立轴对称模型。*HEAT TRANSFER, STEADY STATE步求解。Python脚本自动化参数化研究。
各工具都能轻松验证啊。
临界绝热厚的理论解明确，是工具导入初期的验证案例最佳选择。不吻合就知道有设置错误。
配管设计专用工具
| 工具 | 特点 |
|---|---|
| 3E Plus (NAIMA) | 免费。计算省能量、CO2削减量 |
| ISO 12241计算工具 | 欧洲标准准拼的绝热厚计算 |
| COMSOL Pipe Flow | 配管系统全体热损失1D模型评估 |
| Flownex | 系统分析。可处理整个配管网络 |
装置配管的绝热设计用专用工具更高效啊。
FEM用于局部验证，系统整体用1D工具评估。各得其所很重要。
Rockwool（Paroc）的绝热模拟器
Rockwool集团（丹麦）2018年公开了在线绝热计算工具"Rockwool Insulation Advisor"，只需输入配管径、温度、使用环境就能自动输出推荐绝热厚和LCC削减效果。采用ISO 12241规范的计算逻辑，广泛用于设计初期筛选。
临界绝热厚的先端研究
可变 $h$ 的临界绝热厚
自然对流时 $h$ 不是常数吧。
圆筒外表面自然对流采用Churchill-Chu相关式。
由于 $h$ 依赖外径 $D$，全热阻的导数变复杂。数值求解显示，临界半径比恒定 $h$ 假设更大。
恒定 $h$ 假设不是保守侧啊。
临界半径被低估，不安全。需要精度时推荐用 $h(D)$ 的计算。
绝热材的拓扑优化
近期研究课题是优化绝热材空间分布。非均匀厚，在热流大的地方集中配置。
均匀缠绕相比效率更高吗？
优先在弯曲部和支撑夹具（热桥）配置能更好降低总体热损失。COMSOL拓扑优化或自编伴随法都可实现。
极低温绝热
LNG（-162℃）或液态氢（-253℃）配管绝热是另一个层次的课题。真空多层绝热（MLI）有效导热率 $10^{-4}$ W/(m K) 量级，临界绝热厚可完全忽视。
极低温主要不是结露而是着霜问题吧。
绝热不良处着霜后绝热性能进一步恶化，陷入恶性循环。极低温设备中绝热接缝处理是最重要的施工品质指标。
极低温配管绝热设计
液态氮（−196℃）或液态氢（−253℃）配管为了最小化外部热侵入，采用真空多层绝热（MLI: Multi-Layer Insulation）。NASA在卡纳维拉尔角Artemis火箭液态氢管线采用的MLI绝热材有效导热率低于0.00003 W/m·K，仅为常规玻璃棉的十万分之一。
临界绝热厚的故障排除
常见问题与应对
请告诉与临界绝热厚相关的实务问题。
整理一些典型故障。
1. 绝热效果不理想（细径管）
现象：小径冷媒管套绝热管后结露不止。
原因：$r_i < r_{cr}$ 时绝热材太薄导致放热反而增加。表面温度下降到露点以下。
对策：确保 $r_o > 3 \times r_{cr}$。
2. 样本值与实效值的偏差
绝热材实际性能有时与样本不符吗？
经常出现。主要原因如下。
| 因素 | 样本条件 | 实际 |
|---|---|---|
| 含水 | 干燥 | 含水使 $k$ 增加2~5倍 |
| 经年劣化 | 新品 | 收缩、裂缝导致间隙 |
| 施工品质 | 理想 | 接缝、支撑部位有热桥 |
| 温度 | 常温 | 高温时 $k$ 增大 |
3. FEM结果与理论不符
最常见原因如下。
- 轴对称设置遗漏：Ansys PLANE55需KEYOPT(3)=1，Abaqus选择轴对称单元
- 单位系混乱：mm系需改换 $k$ = 0.0002 W/(mm K)
- 网格不足：薄绝热层仅1层单元
单位系转换错误不易发现啊。
通则是用理论解检算。$q = 2\pi k L(T_i - T_o)/\ln(r_o/r_i)$ 手算值与FEM吻合就没问题。
绝热材含水使λ剧增
玻璃棉绝热材吸水后导热率从干燥时0.04 W/m·K急剧增至0.5~0.6 W/m·K（接近水的值）。化工厂屋外配管保温材因裂缝吸收雨水的案例中，绝热效果降到设计值的1/10以下，结合结露腐蚀最终导致配管更换的大规模维修。
更详细
错误