圆柱坐标系热传导
理论与物理
圆柱坐标系的热传导方程
老师，用笛卡尔坐标分析管道或轴的热问题不行吗？
对于圆柱形状，使用笛卡尔坐标表示形状效率低，网格会变得复杂。使用圆柱坐标 $(r, \phi, z)$ 可以利用对称性降低维度。一维径向稳态热传导方程为
和笛卡尔坐标的傅里叶定律相比，多了 $1/r$ 这一项呢。
这是因为体积元与 $r$ 成正比，是反映面积变化的项。这也是临界绝热半径的物理原因。
解析解（无内部发热）
当 $\dot{q}_v = 0$ 时，解如下。
温度分布是对数函数，与平板的线性分布不同。热流量为
出现 $\ln$ 是因为面积会变化吧。
是的。因为截面积 $A = 2\pi r L$ 与 $r$ 成正比，对于恒定的热流量 $q$，热流密度 $q'' = q/A$ 与 $r$ 成反比。内表面的热流密度更大。
有内部发热的情况
存在均匀内部发热 $\dot{q}_v$ 时的解为
在中心 $r=0$ 处达到最高温度 $T_{\max} = T_s + \dot{q}_v R^2/(4k)$。
电线的焦耳发热就是这种情况吧。
没错。AWG18铜线（直径1.02mm）通过10A电流时发热约10 W/m。$\dot{q}_v \approx 1.2 \times 10^7$ W/m$^3$，可以用这个公式估算中心温升。
圆柱坐标拉普拉斯算子的由来
圆柱坐标系中的热传导方程表示为 ∂²T/∂r² + (1/r)∂T/∂r + (1/r²)∂²T/∂θ² + ∂²T/∂z² + q̇/λ = 0。这个拉普拉斯算子中的 (1/r)∂T/∂r 项源于坐标变换的雅可比行列式。拉梅在1833年完善了弹性论中的圆柱坐标拉普拉斯算子，20年后才被应用于傅里叶热方程，这是其历史背景。
各项的物理意义
- 蓄热项 $\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$）变暖。锂离子电池充放电时的发热、刹车片的摩擦热在分析中也作为热源考虑。与从外部对“表面”加热的边界条件不同，热源项表示“内部”的能量生成。
假设条件与适用范围
数值解法与实现
轴对称FEM模型
用FEM求解圆柱热传导时，需要建立3D模型吗？
如果周向对称，2D轴对称模型就足够了。计算成本可降至1/100以下。
| 模型类型 | 自由度 | 计算时间 | 精度 |
|---|---|---|---|
| 3D整体 | 约100万 | 分钟量级 | 基准 |
| 3D 1/4对称 | 约25万 | 秒量级 | 同等 |
| 2D轴对称 | 约1000 | 瞬时 | 同等 |
2D轴对称就足够却用3D计算，真是浪费啊。
利用对称性是计算效率的基本。Ansys Mechanical中使用PLANE55（KEYOPT(3)=1），Abaqus中使用DCAX4（轴对称4节点单元）。
有限差分法的离散化
圆柱坐标的中心差分需要注意。
在 $r = 0$（中心轴）处，$1/r$ 成为奇点，需应用洛必达法则，得到
。
中心的奇点需要特殊处理呢。
FEM中中心轴的处理同样重要。轴对称单元在 $r = 0$ 节点处会自动进行正确处理，但在3D模型中对轴线上使用楔形单元时需要注意。
多层圆柱的解法
多层圆柱（管道+隔热材料+包覆层等）需连接各层的解析解。
若包含层间接触热阻，则在分母中添加 $1/(2\pi r_i h_{c,i} L)$。
和电路串联电阻的加法一样呢。
正是热阻网络的思想。即使是5层以上的多层结构，也能用于计算，这是圆柱坐标的优势。
对数温度分布的求解步骤
无限长空心圆柱（内半径ri、外半径ro）的稳态温度分布为 T(r)=Ti + (Ti−To)·ln(r/ro)/ln(ri/ro)，呈对数型。其解法分为变量分离→ODE积分→由两个边界条件确定常数C三个步骤完成。日本大学入学考试中也出现过类似题目，国内的热工学教材必定将其作为章首例题。
线性单元 vs 二次单元
热传导分析中，线性单元通常也能获得足够精度。在温度梯度陡峭的区域（如热冲击等）推荐使用二次单元。
热流密度的评估
根据单元内的温度梯度计算得出。与节点应力类似，有时需要进行平滑处理。
对流-扩散问题
当佩克莱数较高（对流主导）时，需要迎风稳定化（如SUPG等）。纯热传导问题则不需要。
瞬态分析的时间步长
设定时间步长应远小于热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$：热扩散率）。对于急剧的温度变化，自动时间步长控制有效。
非线性收敛
由温度相关物性值引起的非线性通常较温和，皮卡德迭代（直接替代法）通常足够。对于辐射的强非线性，推荐牛顿法。
稳态分析的判定
当所有节点的温度变化低于阈值（如 $|\Delta T| / T_{max} < 10^{-5}$ 等）时判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快，但时间步长过大则不稳定（会错过风暴）。隐式法是“也考虑未来状态的预测”——即使时间步长较大也稳定，但每个步长都需要解方程，计算量大。对于没有急剧温度变化的问题，使用隐式法配合较大的时间步长更高效。
实践指南
管道的散热计算
请告诉我实际工作中应用圆柱热传导的典型场景。
最常见的是工厂管道的热损失计算。以蒸汽管道（外径114.3mm、SUS304、隔热材料50mm）为例。
| 层 | $r$ [mm] | $k$ [W/(m K)] | 热阻 [m K/W] |
|---|---|---|---|
| 管壁 | 48.6→57.15 | 16.3 | 0.00160 |
| 隔热材料 | 57.15→107.15 | 0.05 | 2.016 |
| 内表面对流 | — | — | 0.00033 |
| 外表面对流 | — | — | 0.0149 |
总热阻 $R_{\text{total}} = 2.033$ m K/W。隔热材料占总热阻的99%。
管壁的热阻几乎可以忽略不计呢。
是的。金属管的热阻不到隔热材料的1/1000，实际上几乎没有温降。所以管道材质是铜还是不锈钢，对散热量几乎没有影响。
对数平均半径
圆柱单位长度的热阻 $R = \ln(r_2/r_1)/(2\pi k)$ 可以表示为等效平板的 $R = t/(k \cdot A_{lm})$。$A_{lm}$ 是对数平均面积
使用对数平均就能沿用平板的公式了呢。
如果 $r_2/r_1 < 2$，用算术平均 $(r_1+r_2)/2$ 误差也在4%以下。简易计算时通常用算术平均就足够了。
网格划分注意事项
列举圆柱轴对称网格需要注意的几点。
- 径向使用偏置网格（内表面侧更细）
- 薄壁管径向至少3个单元
- 多层结构时，层间共享节点或使用绑定接触
使用偏置网格是因为内表面的热流密度更大吧。
相关主题
なった
詳しく
報告