扩大传热面(翅片)
理论与物理
扩展传热面的基础
老师,散热片的鳍片为什么有效?
根据对流换热的基本公式 $q = hA(T_s - T_\infty)$,当难以增大 $h$ 时,增加面积 $A$ 就是扩展传热面(Extended Surface)的策略。仅靠基面散热不足时,通过鳍片将表面积扩大10~100倍来实现散热。
只要单纯增加面积就行了吧。
但是,由于越靠近鳍片尖端温度越低,所以并非整个表面都以基座温度散热。需要用鳍片效率 $\eta_f$ 来评估其有效性。
鳍片的控制方程
根据能量守恒,均匀截面的直鳍片温度分布为
其中 $\theta = T(x) - T_\infty$,$m = \sqrt{hP/(kA_c)}$。$P$ 是鳍片周长,$A_c$ 是鳍片横截面积。
通解为 $\theta(x) = C_1 e^{mx} + C_2 e^{-mx}$,常数由边界条件决定。
$m$ 是鳍片参数呢。$m$ 越大温度下降越快。
是的。$m$ 表示“鳍片的细长程度”。鳍片越薄越长($A_c$ 小,$P$ 大),$m$ 越大,尖端温度越低。
边界条件与解
| 尖端条件 | 温度分布 | 散热量 |
|---|---|---|
| 绝热尖端 | $\theta = \theta_b \frac{\cosh m(L-x)}{\cosh mL}$ | $q = \sqrt{hPkA_c}\,\theta_b \tanh mL$ |
| 恒定温度尖端 | 双曲函数的线性组合 | 根据情况而定 |
| 对流尖端 | 用修正长度 $L_c = L + A_c/P$ 近似 | 将绝热尖端的公式应用于 $L_c$ |
| 无限长鳍片 | $\theta = \theta_b e^{-mx}$ | $q = \sqrt{hPkA_c}\,\theta_b$ |
实际工作中多用绝热尖端近似吗?
尖端散热只占总量的百分之几,因此使用修正长度 $L_c$ 的绝热尖端近似足以获得足够的精度。
鳍片诞生的历史
扩展传热面(鳍片)的概念由阿尔弗雷德·哈珀于1922年系统化。如今,它也被应用于英特尔的CPU散热器,通过将表面积最大扩展10倍,成为大幅提升散热量的核心技术基础。
各项的物理意义
- 蓄热项 $\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$ 是鳍片整个表面处于基座温度时的最大散热量与实际散热量的比值。
$$\eta_f = \frac{q_{\text{actual}}}{q_{\text{max}}} = \frac{q_{\text{actual}}}{hA_f \theta_b}$$
对于绝热尖端的直鳍片,$\eta_f = \tanh(mL)/(mL)$。当 $mL < 1$ 时,鳍片效率在90%以上;$mL > 3$ 时急剧下降。
$mL$ 是鳍片设计的关键参数呢。
一般认为 $mL \approx 1$ 是成本效益的最优点。超过这个长度,材料增加带来的散热提升就不明显了。
FEM中的鳍片分析
鳍片的FEM分析因其薄壁结构,有时可以用壳单元或梁单元建模,但在热分析中通常使用实体单元。
方法 优点 缺点
3D实体 精度最高 网格数量多
2D截面 高效分析重复结构 忽略3D效应
1D解析解 快速,参数化容易 忽略2D/3D热流
如果散热片有100片鳍片,需要全部建模吗?
利用对称性,只需一片鳍片的模型即可。根据鳍片间距设置对称面,并在对称面上施加绝热条件。FloTHERM或Icepak等软件具有将鳍片阵列作为参数化模型自动生成的功能。
整体鳍片效率
散热片整体的性能用结合了鳍片效率和基座面的整体效率 $\eta_o$ 来评估。
$$\eta_o = 1 - \frac{A_f}{A_t}(1 - \eta_f)$$
$A_f$ 是鳍片表面的总面积,$A_t$ 是总表面积(鳍片+基座裸露部分)。
基座裸露部分按鳍片效率100%处理呢。
没错。使用整体效率,总散热量可以简洁地表示为 $q = \eta_o h A_t \theta_b$。
Coffee Break 闲谈
鳍片效率的计算步骤
鳍片效率η由双曲函数tanh(mL)/(mL)表示。m是鳍片参数,等于√(hP/kA)。铝制针状鳍片(k=237 W/m·K)可实现95%以上的效率,铜(k=401 W/m·K)也能实现类似的高效率。
线性单元 vs 二次单元
在热传导分析中,线性单元通常也能获得足够的精度。对于温度梯度剧烈的区域(如热冲击),推荐使用二次单元。
热流密度的评估
根据单元内的温度梯度计算得出。与节点应力类似,有时需要进行平滑处理。
对流-扩散问题
当佩克莱特数较高(对流主导)时,需要迎风稳定化(如SUPG)。纯热传导问题则不需要。
瞬态分析的时间步长
时间步长应远小于热扩散的特征时间 $\tau = L^2 / \alpha$($\alpha$:热扩散率)。对于急剧的温度变化,自动时间步长控制有效。
非线性收敛
由温度相关物性值引起的非线性通常较为温和,皮卡德迭代(直接替代法)往往就足够了。对于辐射的强非线性,推荐使用牛顿法。
稳态分析的判定
当所有节点的温度变化低于阈值(例如 $|\Delta T| / T_{max} < 10^{-5}$)时,判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快,但时间步长过大则不稳定(会错过风暴)。隐式法是“也考虑未来状态的预测”——即使时间步长较大也能稳定,但每个时间步都需要解方程,计算量大。对于没有急剧温度变化的问题,使用隐式法并采用较大的时间步长更高效。
鳍片效率具体表示什么?
鳍片效率 $\eta_f$ 是鳍片整个表面处于基座温度时的最大散热量与实际散热量的比值。
对于绝热尖端的直鳍片,$\eta_f = \tanh(mL)/(mL)$。当 $mL < 1$ 时,鳍片效率在90%以上;$mL > 3$ 时急剧下降。
$mL$ 是鳍片设计的关键参数呢。
一般认为 $mL \approx 1$ 是成本效益的最优点。超过这个长度,材料增加带来的散热提升就不明显了。
鳍片的FEM分析因其薄壁结构,有时可以用壳单元或梁单元建模,但在热分析中通常使用实体单元。
| 方法 | 优点 | 缺点 |
|---|---|---|
| 3D实体 | 精度最高 | 网格数量多 |
| 2D截面 | 高效分析重复结构 | 忽略3D效应 |
| 1D解析解 | 快速,参数化容易 | 忽略2D/3D热流 |
如果散热片有100片鳍片,需要全部建模吗?
利用对称性,只需一片鳍片的模型即可。根据鳍片间距设置对称面,并在对称面上施加绝热条件。FloTHERM或Icepak等软件具有将鳍片阵列作为参数化模型自动生成的功能。
散热片整体的性能用结合了鳍片效率和基座面的整体效率 $\eta_o$ 来评估。
$A_f$ 是鳍片表面的总面积,$A_t$ 是总表面积(鳍片+基座裸露部分)。
基座裸露部分按鳍片效率100%处理呢。
没错。使用整体效率,总散热量可以简洁地表示为 $q = \eta_o h A_t \theta_b$。
鳍片效率的计算步骤
鳍片效率η由双曲函数tanh(mL)/(mL)表示。m是鳍片参数,等于√(hP/kA)。铝制针状鳍片(k=237 W/m·K)可实现95%以上的效率,铜(k=401 W/m·K)也能实现类似的高效率。
线性单元 vs 二次单元
在热传导分析中,线性单元通常也能获得足够的精度。对于温度梯度剧烈的区域(如热冲击),推荐使用二次单元。
热流密度的评估
根据单元内的温度梯度计算得出。与节点应力类似,有时需要进行平滑处理。
对流-扩散问题
当佩克莱特数较高(对流主导)时,需要迎风稳定化(如SUPG)。纯热传导问题则不需要。
瞬态分析的时间步长
时间步长应远小于热扩散的特征时间 $\tau = L^2 / \alpha$($\alpha$:热扩散率)。对于急剧的温度变化,自动时间步长控制有效。
非线性收敛
由温度相关物性值引起的非线性通常较为温和,皮卡德迭代(直接替代法)往往就足够了。对于辐射的强非线性,推荐使用牛顿法。
稳态分析的判定
当所有节点的温度变化低于阈值(例如 $|\Delta T| / T_{max} < 10^{-5}$)时,判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快,但时间步长过大则不稳定(会错过风暴)。隐式法是“也考虑未来状态的预测”——即使时间步长较大也能稳定,但每个时间步都需要解方程,计算量大。对于没有急剧温度变化的问题,使用隐式法并采用较大的时间步长更高效。
实践指南
散热片的设计指南
设计散热片的鳍片时,应该以什么为基准?
我们来整理一下设计参数及其影响。
| 参数 | 增加时 | 权衡 |
|---|---|---|
| 鳍片高度 | 面积增加,散热提升 | $mL$ 增大,鳍片效率下降 |
| 鳍片数量 | 面积增加,散热提升 | 流道变窄,压力损失增大 |
| 鳍片厚度 | $A_c$ 增加,鳍片效率提升 | 重量增加,流通面积减少 |
| 鳍片间距 | — | 存在最优值(自然对流下为6~12mm) |
鳍片间距有最优值呢。
自然对流下,间距过小则空气无法流动,过大则面积不足。Bar-Cohen和Rohsenow的最优间距关联式为
$L$ 是鳍片高度,$\text{Ra}_L$ 是瑞利数。
材料选择
| 材料 | $k$ [W/(m K)] | 密度 [kg/m$^3$] | 用途 |
|---|---|---|---|
| 铜 C1100 | 398 | 8,960 | 高性能散热片 |
| 铝 A6063 | 200 | 2,700 | 通用散热片 |
| 铝 A1050 | 230 | 2,710 | 压铸鳍片 |
| 石墨 | 150~400(面内) | 2,200 | 薄型散热片 |
铜和铝的 $k$ 值相差约2倍,但密度相差3倍以上呢。
就单位重量的散热性能而言,铝更优。航空和车载领域以铝为主流,在数据中心服务器冷却等重量限制较宽松的场合会使用铜。
制造方法与形状限制
制造方法会限制鳍片的形状。
| 制 |
|---|
相关主题
なった
詳しく
報告