基于构形定律的散热器翅片设计：用于风冷电子设备的CFD方法

《International Communications in Heat and Mass Transfer》：Heat sink fins design based on the Constructal Law: a CFD approach for air-cooled electronic devices

　　

热管理对现代电力电子设备（如中等功率光伏（PV）逆变器）的可靠性和使用寿命构成关键挑战。本研究提出一种自组织生成式设计方法，用于嵌入真实光伏逆变器外壳内的风冷散热器，该方法由构形定律（Constructal Law）确定性驱动。数值模拟采用OpenFOAM框架中的chtMultiRegionFoam求解器进行共轭传热和湍流计算，并以标准κ??湍流模型封闭。研究人员未对预定义几何形状进行调优，而是通过一种专用构造算法，基于局部热力学机会函数逐步扩展换热表面。在四种不同翅片厚度约束下，研究人员系统考察了流体-固体界面的演化过程，涵盖一系列求解路径。结果表明，所有构型的全局热阻均呈现显著持续下降。该趋势建立了一种性能权衡关系：润湿周长的巨大增益完全抵消了因通道水力直径衰减导致的局部对流恶化。在分析的路径中，较薄翅片厚度（1.5 mm和2.5 mm）的构型产生了更优的物理设计，与传统商业基准相比，峰值工作温度显著降低，且总材料体积也大幅减少。本研究证明，允许几何结构为追求更大流通能力而自由生长，为最小化系统热阻提供了一条自然路径，直接体现了构形定律——通过调整流动架构以最大化性能，同时确保材料使用的合理性与减量化。

以下为论文解读文章：

**基于构形定律的风冷散热器翅片设计：面向光伏逆变器的CFD生成式方法**

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

现代电力电子设备持续小型化和功率密度不断提升，对热管理提出了严峻挑战。维持器件工作温度在安全范围内是保证性能与可靠性的关键。被动冷却因其简单和稳健而具有吸引力，但传统被动散热器的效率受限于几何刚性，无法适应运行中不断变化的热流条件。在电子和电力系统中，最普遍且经济有效的热控制方案是混合空气对流（自然对流+强制对流），并通过金属散热器增加有效换热面积。散热器性能与其几何结构密切相关：翅片高度、间距和厚度直接影响对流系数和整体压降。然而，现有优化大多依赖预定义形状，不能反映实际流动系统中观察到的自然几何演化。构形定律（Constructal Law）为热系统流动架构的自然演化提供了理论基础：有限尺寸流动系统若要持续存在，必须自由演化以提供更容易、更大的流动通道。尽管已有研究将构形理论应用于简化或理想化系统，但将其推广至真实风冷电子设备仍存在明显空白。主要研究缺口包括：第一，大多数基于构形理论的分析局限于简化域（如层流通道或均匀边界条件），很少扩展至真实电子外壳中复杂的非均匀流场；第二，基于CFD的优化方法虽能改进预定义散热器几何，但缺乏反映构形定律所预测的自然几何演化的自组织生长机制。为填补这些缺口，本研究提出一种基于构形定律的生长算法，使翅片能够根据局部速度场与温度场乘积所定义的生长函数在电子散热器中逐步涌现，将构形原理与精细CFD建模相结合，使散热器形态自然演化，在真实强制风冷条件下增强对流散热。

**2. 研究内容与方法**

研究对象为一台中等功率光伏逆变器，计算域分为下部流体区域（含风扇和散热器，负责强制对流）和上部流体区域（包围电子板和无源元件，主要受自然对流影响），两区域由铝基板隔开。下部区域包含两个轴流风扇和两个散热器（HS），IGBT模块作为主要热源位于固体结构下部，与散热器紧密热接触。散热器翅片并非预定义，而是由构形生长算法逐步生成。

研究人员采用OpenFOAM v10软件，利用snappyHexMesh生成网格，使用chtMultiRegionFoam求解器求解共轭传热问题。流体域假设为不可压缩湍流，采用RANS方程结合Boussinesq近似处理浮力，标准κ??湍流模型封闭。固体域求解导热方程，流体-固体界面通过共轭传热（CHT）条件耦合。边界条件设置中，风扇1采用恒定质量流量，风扇2和开放口采用总压力条件，使流动能响应翅片几何演化引起的压降变化。热载荷包括IGBT 2（70 W）、IGBT 3和4（各35 W）、电感线圈（2.62 W）和电感器（18 W）。

核心方法是一种进化式构造算法。该算法从初始几何（商用散热器基座）出发，每步进行CFD模拟，然后计算构造函数F_obj = U_local × T_local（局部速度与局部温度的乘积），该函数基于牛顿冷却定律的物理启发，指示局部换热改善的机会。算法映射F_obj场，选择最大值作为候选点，并通过决策树判断：若候选点位于已有翅片正上方则垂直延伸；否则检查是否有空间和阈值（F_obj > 400 K·m/s）创建新翅片，并验证最小间距约束（d_bf > 1 mm）。循环重复直到无新点满足条件。研究考察了四种翅片厚度（1.5 mm、2.5 mm、4 mm、6 mm），每种厚度产生唯一的确定性生长路径。

**3. 主要结果**

**3.1 形态演化与生长动力学**
从仅有铝基板的初始状态出发，算法首先在最强的风扇射流路径中创建翅片。生长过程中，翅片高度和位置非均匀演化，适应逆变器内部的非均匀流场。随着翅片增加，基板温度从约157 °C降至40 °C左右，最终形态呈现冷蓝色调，验证了构造函数逻辑的有效性。

**3.2 内部温度分布与参数比较**
较薄的翅片（1.5 mm）允许更高密度的换热表面而不违反间距约束，冷空气能更深入穿透翅片间，温度分布更均匀。6 mm和4 mm几何在结构间出现热空气积聚，局部热饱和。1.5 mm构型中，热点转移到电感器而非IGBT，表明优化成功降低了半导体峰值温度。与商业基准相比，1.5 mm设计在约第10次迭代实现性能反超，最终最大温度降低17 °C。

**3.3 性能演化与无量纲数**
所有构型的全局热阻随构造步骤持续下降，但较薄翅片（1.5 mm和2.5 mm）的下降曲线更陡且持续更久。Reynolds数（Re）和Nusselt数（Nu）随分支增加而下降，因为翅片增加导致流动截面减小、水力直径（D_h）降低。然而，总传热量从约54 W增至69 W，表明换热面积增加带来的收益完全补偿了局部对流系数下降。对于2.5 mm厚度，总热流率持续增长而热流密度从2500 W/m2降至约400 W/m2，表明面积扩大降低了单位面积热负荷，有利于减少热梯度。

**3.4 热效率与材料用量的权衡**
1.5 mm构型用最少的材料（约800 cm3）实现了最佳热性能（最大温度40.3 °C），而6 mm构型用了约1000 cm3材料，性能却最差（57.9 °C）。效率指标（EI）定义为传热量与最大温度和材料体积乘积的加权平均，1.5 mm构型在步骤21达到最高值约4.55×103 W/(m3·°C)，远超商业基准（2.36×103）。热阻-体积乘积（R_thV）随演化下降，但其局部最小值不一定对应足够低的峰值温度，不能单独作为停止准则。

**3.5 流动动力学与速度场重构**
最终构型在薄翅片间形成多个平行的窄通道，速度场更均匀且峰值速度约6 m/s，而初始构型存在大面积滞止区。流线显示最终构型的气流高效穿透翅片，商业散热器则表现出更混沌的流线和回流区。体积流量分析表明，算法不仅扩大了传热面积，还最小化了全局流动阻力：初始时风扇2出现负流量（回流），随着翅片演化，流量逐渐转为正值并稳定。1.5 mm设计稳定在约6 L/s，而6 mm模型在步骤10达到峰值6.3 L/s后骤降至5.6 L/s，表明过厚翅片导致气动壅塞，削弱冷却效果。

**4. 讨论与结论**

讨论部分指出，最终散热器几何取决于边界条件、热源分布、流场和几何约束，但构造程序具有通用性：材料根据每步计算的温度和速度场逐步分布。非均匀翅片分布是耦合求解过程的自然结果。局部传热参数（如Nu）下降并不妨碍整体热阻降低，因为面积增加补偿了局部性能损失。翅片厚度影响几何自由度：较薄元素允许更精细的材料分配，但并不意味着在所有应用中必然更优。材料数量的增加并不自动导致温度降低，材料投放位置对性能的影响更为关键。因此，可推广的贡献不在于具体翅片尺寸或布置，而在于确定性构造程序。

研究结论：本研究提出了一种基于构形定律的新型生长算法，用于光伏逆变器风冷散热器的生成式设计。通过将进化式翅片构造算法与高保真CFD模拟（OpenFOAM中的chtMultiRegionFoam求解器）耦合，该方法使散热器形态从初始原始几何确定性演化，逐步促进热流通道，符合构形定律的基本原理。结果表明，算法成功生成了适应逆变器外壳复杂非均匀流场的有效翅片架构。在评估的翅片厚度中，1.5 mm构型表现出最佳整体热性能，最高温度达40.3 °C，比商业散热器降低约17 °C（对关键IGBT模块），同时温度分布更均匀、有效换热面积显著增加且未造成过度流动阻塞。未来工作应扩展至分支（树状）架构，并评估实际光伏系统典型的瞬态和变负荷条件。该研究证明，将构形定律与CFD驱动的几何生长相结合，是设计高性能风冷散热器的有效方法，超越传统参数优化，为电力电子领域更高效、紧凑和可持续的冷却方案开辟了新途径。