基于CFD的离子风用于高功率电子器件高效热管理研究

《Electronics》：CFD-Based Study of Ionic Wind for Efficient Thermal Management of High-Power Electronics

　　

高功率电子元件的高效热管理已成为一项关键的工程挑战，随着功率密度增长和器件几何尺寸缩小。基于被动散热器和机械风扇的传统解决方案在噪声、可靠性和紧凑性至关重要的应用中日益不足。本论文提出了一项全面的计算流体动力学（CFD）研究，探讨离子风——一种电液动力学（EHD）现象，其中非对称电极之间的电晕放电产生无移动部件的定向气流——作为高功率电子器件冷却的高能效替代方案。研究人员在COMSOL Multiphysics中开发了一个完全耦合的二维多物理场模型，集成了静电场、离子输运（Nernst-Planck方程）、Navier-Stokes流体动力学和对流换热。二维公式虽然计算效率高且与先前的EHD建模研究一致，但忽略了真实三维针形配置固有的横向射流扩散，因此预计会高估峰值冲击速度；研究人员在解释与实验温度的定量比较时考虑了这一局限性。研究聚焦于应用于加热铝板的针-收集器配置，该铝板代表高功率电子元件，如发光二极管（LED）、功率晶体管或微处理器芯片。模拟结果通过与高功率板上芯片（COB）LED系统的纹影光学实验数据进行了间接验证。离子风将最高表面温度降低了8.1 K，并显著削弱了中心热点，使热量横向重新分布。系统的参数研究表明，施加电压和针尖距热源的高度是主要设计参数，而能量平衡显示，EHD射流直接排除了约1.8%的生成热量——主要作为表面对流增强器，而非整体热提取器。这些发现为任何由EHD系统冷却的功率电子元件提供了定量设计指南。

**基于CFD的离子风冷却高功率电子器件研究解读**

**1. 研究背景与问题**

高功率电子器件（如LED模块、功率晶体管、微处理器芯片）的热管理已成为现代电子学中的关键约束。随着功率密度持续增长和器件几何尺寸缩小，传统被动散热器（基于传导和自然对流）和机械风扇主动冷却方案在噪声、可靠性、紧凑性和功耗方面日益不足。尤其在高功率照明应用中，LED可将其输入功率的50%至85%转化为热量，需高效散热以维持发光效率、色稳定性和寿命。国际能源署估计仅照明就占全球电力消耗约15%，而热管理及相关冷却基础设施占高密度功率电子系统总成本的30-40%。半导体结温升高会降低效率、加速老化并增加早期失效风险。

离子风，即电液动力学（EHD）流动，提供了一种根本不同的路径：通过施加高压直流电于尖锐发射电极（针、线或尖端）与接地收集器（板、环或圆柱）之间，在发射极附近产生强电场电离空气分子，离子被加速向收集器，与中性分子碰撞传递动量产生定向气流。该过程无移动部件、无噪声，且具有最高效率。然而，现有研究存在以下空白：缺乏统一的多物理场框架系统研究EHD详细流场结构与热性能的关系；关键设计参数（电压、电极间距、针高、功率密度）的相对贡献不明确；冷却机制中直接热移除与对流阻降低的物理分解未量化；缺乏针对针-收集器配置用于平板功率元件的纹影光学验证。本研究旨在通过完全耦合的COMSOL多物理场模型、完整参数化研究、定量能量平衡及与高功率LED系统的纹影实验对比，填补这些空白。

**2. 主要关键技术与方法**

研究人员采用以下关键技术方法：基于COMSOL Multiphysics开发完全耦合的二维多物理场模型，耦合静电场泊松方程、离子输运Nernst-Planck方程、不可压缩Navier-Stokes方程组（含EHD库仑体力项ρ_eE）和能量方程（对流与传导）。使用PARDISO直接求解器，相对收敛容差10^-5，实现全双向耦合。采用非结构化三角网格，并在针尖、收集器表面和加热板表面局部加密，通过网格独立性验证确保差异小于0.3 K。离子注入采用Peek经验定律。实验参考来自Bridgelux BRXA-W3500 COB LED模块（11 W耗散），使用纹影光学可视化流场，热电偶记录壳温，并评估热阻。通过参数化扫描研究施加电压（8-16 kV）、电极间距（4-6 mm）、针尖高度（4-6 mm）和体积热源功率密度的影响，并进行能量平衡分析。

**3. 研究结果**

**3.1 气流结构：自然对流与离子风对比**
模拟显示，离子风激活后，EHD体力从针尖向下加速空气，产生定向射流，峰值速度1.5-2 m/s，冲击加热板后形成径向壁面射流。该结构匹配纹影图像中观察到的向下射流和横向扩散特征，与现有文献报道一致。无离子风时，空气近乎静止，仅依赖传导。

**3.2 温度分布与热点抑制**
无离子风时，板中心出现强烈热点，垂直热羽流；有离子风时，冲击射流偏转热羽流，热点被削弱，温度分布更均匀，梯度降低。最大表面温度降低8.1 K（从459.4 K降至451.3 K）。根据阿伦尼乌斯模型，10 K降幅可使器件寿命约延长一倍，因此该温度降低对可靠性有显著意义。

**3.3 能量平衡**
通过控制面通量积分计算，离子风射流直接携带的热量仅占板内总热生成的约1.8%。其余98.2%的热量通过增强表面对流（降低热阻）而非直接体热移除得到改善。这表明EHD冷却本质上是边界层管理技术，核心机制是减薄热边界层，而非直接带走大量热量。

**3.4 参数化研究：施加电压**
电压从8 kV升至16 kV，T_max从454 K降至449 K，但12 kV以上增益饱和，归因于空间电荷限制效应。存在最佳工作窗口（8-12 kV）。

**3.5 参数化研究：电极间距与针尖高度**
电极间距4 mm时T_max=450.5 K，5 mm时451.3 K，6 mm时452 K，存在最优间距以平衡射流集中性与覆盖范围。针尖高度4 mm时T_max=450 K，6 mm时452 K，较低针尖提供更聚焦的冲击。

**3.6 参数化研究：耗散功率**
随功率增加，T_max上升，但离子风带来的相对降温（约8 K）在所有测试功率水平下保持恒定，说明其冷却效益与功率水平自然成比例。

**3.7 对流换热增强定量评估**
使用层流冲击射流关联式（Nu_EHD=0.5 Re^0.5 Pr^0.37）估算，EHD贡献的对流换热系数h_EHD从8 kV时的15.3 W/(m²·K)增至16 kV时的18.2 W/(m²·K)，总换热系数h_total增强12.7%-15.2%。对应的等效温度降幅为18.7-21.9 K（基于平均系数），高于CFD热点降幅（8.1 K），反映了空间平均与局部热点的差异。

**4. 讨论与结论**

讨论部分指出，离子风并非通过热力学循环带走热量，而是作为主动边界层管理工具，用结构化快速流动替代缓慢的自然对流层。设计时应关注热阻而非总功耗。电极间距和针尖高度对冷却空间定位的影响超过电压饱和后的增量。研究结论适用于任何产生集中热量的平面源、目前依赖自然对流、需静音或密封运行的电子器件（如功率晶体管、IGBT模块、GaN-on-SiC放大器等）。主要局限性是二维模型忽略了三维射流横向扩散，导致高估中心速度并低估侧向扩散，因此绝对温度降幅需通过三维模拟或粒子图像测速（PIV）测量验证。研究结论如下：

（1）离子风重构近表面气流，以定向冲击射流替代自然对流，减薄热边界层并抑制热点，模拟流场（向下射流及径向壁面流，速度1.5-2 m/s）与纹影观测定性一致。

（2）主要机制为对流阻降低，而非体热移除；仅1.8%生成热量直接被射流带走，其余通过增强表面换热系数实现。

（3）施加电压和针尖高度是主要设计参数；8-12 kV电压、4-5 mm针高提供最佳平衡。

（4）冷却效益与耗散功率成比例，使离子风成为自然可扩展的解决方案。

（5）实验参考证实了卓越能效：稳定温升降低56.4%，热阻从5.2 °C/W降至2.73 °C/W，EHD功耗低于1 mW。这些结果确立了离子风作为中等功率电子器件（约10 W，~2-5 W/cm²热流密度）技术成熟且能效高的冷却策略。未来工作将开发全三维耦合EHD-热模型，进行PIV速度测量，并探索与优化散热器几何的集成。