4. 网格划分策略与质量评估

风冷翅片与液冷微通道虽然都属于“固体壁面与流体换热”的问题,但几何尺度差异极大:翅片间距通常在 2~10 mm,而液冷微通道水力直径可能只有 0.3~1 mm,通道高度甚至低至 0.2 mm。这种尺度差异直接决定了网格策略必须分开设计,不能套用同一套尺寸参数。本节围绕几何分解、边界层处理、独立性验证和质量指标四个环节展开,并给出可直接参考的参数设置与判断依据。

一、几何分解与网格类型选择

风冷翅片散热器通常由基板、若干平行翅片和外部空气域组成。空气域外边界应距离翅片顶部和侧面至少 5 倍翅片高度,避免出口边界影响回流区。对于这种“开放外流+规则翅片”的结构,推荐采用多面体网格或六面体主导网格:翅片间流道规则,可先对单个流道做扫掠(Sweep)生成六面体,再通过周期或复制扩展到全部流道。多面体网格在相同单元数下对复杂外形的适应性更好,且收敛性通常优于四面体,适合翅片顶部圆角、基板倒角等细节。

液冷微通道则完全不同。微通道内部流动接近充分发展层流,壁面法向速度梯度极大,必须使用结构化六面体网格。若通道为矩形直通道,可沿流动方向扫掠;若带有扰流柱、凹槽或分形结构,则应在通道核心区保留六面体,仅在局部复杂区域使用多面体过渡。需要特别注意的是,微通道的宽高比可能达到 5:1 甚至 10:1,此时应控制单元长宽比不超过 20,否则离散误差会显著增大。

二、边界层处理:第一层高度与增长率

边界层网格是两类散热器网格划分中最容易出问题的环节。核心参数是壁面第一层高度 Δy,它由目标 y+ 决定。对于风冷翅片,若采用 k-ε 或 k-ω SST 湍流模型并配合壁面函数,目标 y+ 通常取 30~100;若希望解析到黏性底层,则需 y+≈1。对于液冷微通道,流动多为层流,不存在壁面函数适用性问题,但仍需保证壁面法向至少有 8~15 层网格,以捕捉抛物线速度剖面。

第一层高度可用平板边界层估算式反推:Δy = y+ · ν / uτ,其中 uτ = √(τw/ρ)。工程上更常用的是基于雷诺数的经验估算。以水为例,通道水力直径 0.5 mm、入口流速 1 m/s、20 ℃ 水 ν=1.0×10⁻⁶ m²/s,则 Re≈500,属于层流。若希望壁面第一层 y+≈1,可估算 Δy 约为 1~2 μm。这个尺度意味着微通道网格总量会迅速上升,因此必须配合增长率控制:边界层增长率取 1.1~1.2,总层数 10~15 层,既保证分辨率又避免单元数爆炸。

风冷翅片的边界层相对宽松。以翅片高度 30 mm、风速 3 m/s 为例,若采用 SST k-ω 模型并解析边界层,第一层高度可设为 0.05~0.1 mm,增长率 1.2,层数 8~12 层。若改用壁面函数,第一层高度可放宽到 0.3~0.5 mm,但必须确认翅片表面 y+ 落在 30~100 区间,否则壁面剪切和换热系数都会失真。

边界层网格生成后,务必在 Fluent 中查看壁面 y+ 云图,而不是只看网格软件里的统计值。风冷翅片重点检查翅片前缘和尾缘,液冷微通道重点检查通道入口段和拐角处,这些位置最容易出现 y+ 超限。

三、网格独立性验证

网格独立性验证不是“把网格加密一遍看结果变不变”,而是要有明确的加密序列和判据。推荐至少做三套网格:粗、中、细,单元数比例约为 1:2:4。对风冷翅片,监测量选散热器总热阻 Rth 和压降 Δp;对液冷微通道,监测量选平均努塞尔数 Nu 和通道压降 Δp。判据是:相邻两套网格的监测量相对偏差小于 2%,且细网格相对中网格的变化小于中网格相对粗网格的变化,即可认为中网格已进入独立区间。

实际操作中,可用 Richardson 外推估算离散误差。设细网格结果为 f3、中网格为 f2、粗网格为 f1,网格细化比 r=2,则表观收敛阶 p≈ln[(f1-f2)/(f2-f3)]/ln r。若 p 接近模型理论阶数,说明收敛正常;若 p 异常偏小,往往意味着边界层分辨率不足或局部网格畸变,需要回到边界层参数重新调整,而不是继续全局加密。

风冷翅片与液冷微通道网格策略对比与质量评估链 风冷翅片散热器 几何特征:翅片间距 2~10 mm,翅片高 20~50 mm 空气域外边界 ≥ 5 倍翅片高度 推荐:多面体 / 六面体主导 + 扫掠 边界层参数(示例) 湍流模型:SST k-ω 第一层高度 Δy:0.05~0.1 mm(解析边界层) 增长率 1.2,层数 8~12;壁面函数时 y+ 30~100 独立性验证监测量 总热阻 Rth、压降 Δp 粗:中:细 ≈ 1:2:4,相邻偏差 < 2% 重点检查:翅片前缘、尾缘 y+ 云图 质量指标 正交质量 > 0.15,偏斜度 < 0.85 长宽比 < 20,单元数约 300 万~800 万 液冷微通道冷板 几何特征:水力直径 0.3~1 mm,高 0.2~0.5 mm 宽高比可达 5:1~10:1 推荐:结构化六面体 + 局部多面体过渡 边界层参数(示例) 流动状态:层流,Re ≈ 500(D=0.5 mm,v=1 m/s) 第一层高度 Δy:1~2 μm(目标 y+≈1) 增长率 1.1~1.2,层数 10~15 独立性验证监测量 平均努塞尔数 Nu、通道压降 Δp Richardson 外推估算收敛阶 p 重点检查:入口段、拐角处 y+ 云图 质量指标 正交质量 > 0.2,偏斜度 < 0.8 长宽比 < 20,单元数约 500 万~1500 万 尺度差异

四、质量指标与常见问题处理

Fluent 中常用的网格质量指标包括正交质量(Orthogonal Quality)、偏斜度(Skewness)和长宽比(Aspect Ratio)。对风冷翅片,正交质量建议大于 0.15,偏斜度小于 0.85;对液冷微通道,由于结构化程度高,正交质量应大于 0.2,偏斜度小于 0.8。若出现大量低质量单元,优先检查翅片顶部圆角、微通道拐角和进出口连接段,这些位置往往因几何倒角过小或扫掠方向冲突而产生畸变。

另一个容易被忽视的问题是边界层与核心区网格的过渡。若增长率过大(如 1.5 以上),边界层最后一层与核心区第一层尺寸突变,会在交界处产生虚假扩散,表现为壁面换热系数偏低或压降偏高。处理方法是把增长率降到 1.2 以下,或在边界层与核心区之间插入一层过渡网格。对于微通道,若单元数受限,可适当减少边界层层数但保持第一层高度不变,因为第一层高度对 y+ 的影响远大于层数。

不要用“全局加密”代替“局部修正”。若 y+ 超限,应调整第一层高度或边界层层数;若某处偏斜度超标,应局部重构该区域网格。全局加密不仅计算成本高,还可能因为长宽比恶化而降低精度。

五、可参考的网格设置检查表

检查项 风冷翅片 液冷微通道
网格类型 多面体 / 六面体主导 结构化六面体为主
第一层高度 0.05~0.1 mm(解析边界层) 1~2 μm(目标 y+≈1)
边界层增长率 / 层数 1.2 / 8~12 1.1~1.2 / 10~15
独立性监测量 Rth、Δp Nu、Δp
正交质量下限 > 0.15 > 0.2
偏斜度上限 < 0.85 < 0.8
长宽比上限 < 20 < 20
重点检查位置 翅片前缘、尾缘 通道入口段、拐角

实际项目中,建议先按上表生成一套中等密度网格,完成初步求解并查看 y+ 云图和残差曲线,再根据监测量偏差决定是否加密。若中网格与细网格的 Rth 或 Nu 偏差已小于 2%,且压降偏差小于 3%,即可认为网格满足独立性要求,不必继续加密。这样既控制了计算成本,也保证了风冷与液冷两类散热器在后续性能优化中的结果可信度。