手机版
           

Fluent传热仿真:共轭传热建模与热管理设计实战

发布时间:2026-06-24   来源:科研学术网    
字号:

Fluent传热仿真中,项目组在IGBT功率模块散热设计上经历了一次方案反转:初始设计的强制风冷方案仿真给出芯片结温约110°C——低于150°C的安全限值。但客户提出要评估自然对流方案以降低风扇噪音和故障率。自然对流的仿真给出结温约165°C——超标。项目组通过优化散热器翅片参数(翅高增加、翅间距减小),在自然对流下将结温降到约140°C。两种方案的仿真精度通过红外热像仪实测验证,偏差在数°C以内。

共轭传热(CHT)的建模要点

功率模块散热涉及固体(芯片、基板、散热器)和流体(空气)两个域的耦合传热——这就是共轭传热(Conjugate Heat Transfer, CHT)。Fluent中CHT的关键设置:

1. 计域定义

  • 固体域:芯片(Si)、焊料层(SnAgCu)、铜基板、铝散热器
  • 流体域:散热器周围空气
  • 耦合壁面:固体-流体界面(Fluent中用”coupled”壁面热边界条件)

2. 材料属性: 各材料的热导率对结果影响最大。项目组使用的参数:

材料 热导率(W/m·K) 比热(J/kg·K) 密度(kg/m³)
Si (芯片) 148 700 2329
SnAgCu (焊料) 55 230 7400
Cu (基板) 398 385 8978
Al (散热器) 201 871 2719
空气 温度依赖 1006 Boussinesq近似

注意空气的热导率是温度的函数——在自然对流中温差大(ΔT>50K),如果用常数热导率会引入误差。Fluent中设置air为”ideal gas”或使用温度相关的多项式热导率。

3. 热源定义: IGBT芯片的功耗为150 W(稳态),分布在8.5×8.5 mm的芯片面积上。热流密度q”=150/(8.5e-3)²=2.08×10⁶ W/m²。在Fluent中设为固体域内的体积热源。

自然对流 vs 强制对流:方案对比

强制对流方案

  • 风扇:40mm轴流风扇,风量25 CFM
  • 入口边界:速度入口(3.5 m/s,298K)
  • 出口边界:压力出口(0 Pa表压)
  • 湍流模型:k-ω SST(Re≈8000,过渡区)
  • 辐射模型:S2S(Surface-to-Surface),表面发射率ε=0.9(阳极氧化铝)

仿真结果:芯片结温约110°C,散热器基板温度约80°C,出口空气温升约45K。散热器的总热阻约0.58 K/W。

自然对流方案

  • 无风扇,靠浮力驱动
  • 入口/出口:压力边界(0 Pa表压)
  • 流体密度:Boussinesq近似(ρ = ρ₀[1-β(T-T₀)]),β=3.3×10⁻³ K⁻¹
  • 重力:-9.81 m/s²(z方向)
  • 湍流模型:层流(Gr≈10⁶,自然对流处于层流-过渡区)
  • 辐射模型:S2S

仿真结果(优化前):芯片结温约165°C,散热器基板温度约130°C。超标。

自然对流的传热系数远低于强制对流(自然对流h≈5-15 W/m²K,强制对流h≈30-80 W/m²K)。要弥补这个差距,必须增大散热面积。

散热器翅片优化

项目组用参数化扫描优化翅片几何:

翅高(mm) 翅间距(mm) 翅片数 结温(°C, 相对趋势)
30 8 12 最高
30 6 16 降低
35 6 16 继续降低
40 5 19 继续降低
45 5 19 最优
50 5 19 边际改善
45 4 24 反弹升高

从30mm/8mm到45mm/5mm,结温持续下降。翅片进一步增加(间距从5mm降到4mm),结温反而升高——翅间距太小时,翅间的空气流动受阻(粘性效应主导),有效传热面积虽然增大但传热效率下降。

最优参数:翅高45mm、翅间距5mm。结温降至约140°C,低于150°C限值约10°C——安全裕量足够。

翅高从45mm增到50mm改善有限——边际收益递减。翅高再增加会增加体积和重量,不划算。

热接触阻的模拟

IGBT模块中芯片-焊料-基板之间存在界面热阻,简称热接触阻(Thermal Contact Resistance, TCR)。如果不模拟TCR,仿真会低估结温10-20°C。

项目组在Fluent中用薄壁热阻层(Thin Wall Thermal Resistance)模拟TCR:

R_tc = t_eff / k_eff

各界面热阻值:

  • Si-焊料:0.5×10⁻⁴ m²K/W
  • 焊料-Cu基板:0.8×10⁻⁴ m²K/W
  • Cu基板-导热硅脂:2.0×10⁻⁴ m²K/W
  • 导热硅脂-散热器:2.0×10⁻⁴ m²K/W

导热硅脂的接触阻最大(4.0×10⁻⁴ m²K/W),是因为硅脂层厚度约0.1mm,热导率仅1.5 W/mK。这个接触阻贡献了约0.27 K/W的热阻——占总热阻的35%。客户后来换用了热导率3.5 W/mK的相变材料,接触阻降到1.7×10⁻⁴ m²K/W,结温再降4°C。

瞬态传热分析

稳态仿真给出的是最坏情况结温(热平衡态),但实际工况中功率是脉冲的。项目组做了瞬态仿真分析热响应时间:

  • 时间步长:0.1 s(自然对流),0.01 s(强制对流)
  • 总模拟时间:600 s
  • 功率曲线:0-60s为150W,60-120s为0W(待机),周期重复

瞬态结果显示:强制对流方案的热时间常数约120s(结温从25°C升到稳态112°C的63%用时120s),自然对流方案约300s。这意味着在间歇工作模式下(占空比50%),强制对流的峰值结温约95°C(低于稳态值17°C),自然对流约128°C(低于稳态值14°C)——间歇工况比连续工况有利。

辐射传热的作用

在自然对流方案中,辐射传热占总散热量的15-20%——不可忽略。项目组对比了开启/关闭辐射模型的结果:

  • 不开辐射:结温偏高约10°C
  • 开辐射(ε=0.9):结温正常

差距约10°C——辐射对自然对流散热方案至关重要。但如果散热器表面是抛光铝(ε=0.05),辐射贡献几乎为零。项目组建议客户保持阳极氧化处理(ε=0.9)以利用辐射散热。

反思:仿真的验证与可信度

CFD传热仿真的可信度依赖实验验证。项目组用红外热像仪(FLIR T660)实测了散热器表面温度:

测点 仿真温度(°C) 实测温度(°C) 偏差
散热器基板中心 ~78 ~80 ~2
翅片根部 ~76 ~78 ~2
翅片顶端 ~62 ~64 ~2
基板边缘 ~68 ~67 ~1

最大偏差约2°C——仿真精度在工程可接受范围内。偏差来源包括:热接触阻值的不确定性、环境温度波动、红外发射率标定误差。

Fluent传热仿真的价值不在于”绝对精确”,而在于”相对正确”——方案对比和参数优化的方向是可靠的。更多热管理设计的仿真经验,可以参考有限元仿真栏目,或返回科研学术网首页。

图说天下

×
abaqus仿真
ansys仿真
comsol仿真
fluent仿真
力学仿真
多相流仿真
流体/流动仿真