扩展传热面综合评估
理论与物理
翅片阵列的性能评估
老师，我明白单个翅片的效率，但整个翅片阵列的性能该如何评估呢？
用总表面效率（Overall Surface Efficiency）$\eta_o$ 来评估。这是翅片表面和基座暴露表面合计的有效散热能力。
其中 $N$ 是翅片数量，$A_f$ 是单个翅片面积，$A_t = N A_f + A_b$ 是总表面积，$A_b$ 是基座暴露面积。
如果 $\eta_f = 1$，那么 $\eta_o = 1$ 呢。
是的。翅片效率越差，$\eta_o$ 就越低。实用的散热器中 $\eta_o$ 大约在 0.7 到 0.9 之间。
整体热阻
散热器整体的热阻是
$R_{\text{spread}}$ 是扩展热阻（当热源小于基座时产生），$R_{\text{base}}$ 是通过基座厚度的传导热阻。
扩展热阻是什么？
当CPU芯片（30mm见方）安装在散热器基座（60mm见方）中央时，热量在基座内横向扩散所产生的阻力。可以用Song-Lee-Au公式估算。
$a$ 是热源的等效半径，$\Psi$ 是基座厚度和毕渥数的函数。
翅片有效性
除了翅片效率，翅片有效性 $\varepsilon_f$ 也是一个重要指标。
表示与无翅片情况（仅靠翅根面积 $A_c$ 散热）相比的散热倍数。一般认为 $\varepsilon_f > 2$ 时添加翅片才有益处。
翅片效率和有效性是不同的东西呢。容易混淆。
效率是“翅片面积的有效百分比”，有效性是“添加翅片后散热能力提高了多少倍”。有时一个高而另一个低。
多种翅片的分类体系
扩展传热面可分为针状、矩形、三角形、环状四大形状。20世纪80年代Kays和伦敦在《紧凑型热交换器》中总结的NTU-ε法成为设计基础，至今仍在航空航天、化工厂领域被标准性地参考。
各项的物理意义
- 蓄热项 $\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$）变暖。锂离子电池充放电时的发热、刹车片的摩擦热在分析中也作为热源考虑。与从外部向“表面”施加热量的边界条件不同，热源项表示“内部”的能量生成。
假设条件与适用范围
数值解法与实现
各种翅片形状的效率
翅片的截面形状会影响效率吗？
影响很大。比较一下代表性翅片形状的效率。
| 翅片形状 | 效率公式 | 材料用量 | 制造性 |
|---|---|---|---|
| 矩形 | $\eta_f = \tanh(mL_c)/(mL_c)$ | 基准 | 容易 |
| 三角形 | $\eta_f = I_1(2mL)/(mL \cdot I_0(2mL))$ | −50% | 较困难 |
| 抛物线形 | 贝塞尔函数 | −67% | 困难 |
| 环状（圆盘） | 修正贝塞尔函数 | 视情况而定 | 容易 |
三角形翅片能节省一半材料呢。
根部厚、尖端薄的三角形翅片，其材料分配符合温度分布，效率更高。矩形翅片的尖端部分温度低，对散热贡献不大，因此浪费较多。
环状翅片的解析
安装在管道或圆筒外表面的环状翅片（环形翅片）需要用贝塞尔函数求解。
$r_1$ 是内半径（根部），$r_2$ 是外半径（尖端）。$I_0, I_1, K_0, K_1$ 是修正贝塞尔函数。
贝塞尔函数手算起来很麻烦呢。
使用效率图表（Incropera教材中的Fig. 3.20）更实用。根据 $r_2/r_1$ 和 $mL_c$ 两个参数从图表中读取。用Python的话，可以直接用scipy.special.iv/kv计算。
CHT分析与比较
解析解基于 $h$ 均匀的假设。实际的翅片阵列中，由于流动发展和涡流的产生，$h$ 变化很大。通过CFD进行共轭传热（CHT）分析自动计算局部 $h$，得到的散热量通常比解析解低10〜20%。
解析解是乐观的估算呢。
所以，先用解析解进行初步设计，再用CHT分析进行最终确认的两阶段方法更实用。
环状翅片的解析方法
环状翅片的效率由包含贝塞尔函数I₀, I₁, K₀, K₁的公式给出。当外半径/内半径比r₂/r₁为2.0时，效率比矩形翅片低约5〜8%。用MATLAB的besselj函数只需几行代码即可计算。
线性单元 vs 二次单元
在热传导分析中，线性单元通常足以获得足够的精度。对于温度梯度陡峭的区域（如热冲击），推荐使用二次单元。
热流密度的评估
根据单元内的温度梯度计算得出。有时需要像节点应力那样进行平滑处理。
对流-扩散问题
当佩克莱特数高（对流主导）时，需要迎风稳定化（如SUPG）。纯热传导问题则不需要。
非稳态分析的时间步长
时间步长应相对于热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$：热扩散率）足够小。对于急剧的温度变化，自动时间步长控制有效。
非线性收敛
由温度依赖性物性值引起的非线性通常较温和，皮卡迭代（直接替换法）通常就足够了。对于辐射的强非线性，推荐使用牛顿法。
稳态分析的判定
当所有节点的温度变化低于阈值（例如 $|\Delta T| / T_{max} < 10^{-5}$）时，判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快，但时间步长大时不稳定（会漏掉风暴）。隐式法是“也考虑未来状态的预测”——即使时间步长大也能稳定，但每个时间步都需要解方程，比较耗时。对于没有急剧温度变化的问题，使用隐式法配合较大的时间步长更高效。
实践指南
散热器选型流程
在实际工作中选择散热器时，应该按什么步骤进行？
标准流程是这样的。
1. 确定发热量：确认TDP（热设计功耗）
2. 确认允许温度：结温上限（例如：$T_j \leq 105$℃）
3. 计算允许热阻：$R_{\text{hs}} \leq (T_j - T_a)/Q - R_{jc} - R_{\text{TIM}}$
4. 选择散热器候选：从产品目录中选择 $R_{\text{hs}}$ 以下的型号
5. CFD验证：结合实际安装环境（风速、相邻部件）进行验证
$R_{jc}$ 是封装的热阻吧。
是的。$R_{jc}$（结到外壳）在IC制造商的数据手册中有记载。$R_{\text{TIM}}$（热界面材料）大约在0.1〜0.5 K/W 左右。
与实测数据的比较
なった
詳しく
報告