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 区间,否则壁面剪切和换热系数都会失真。
三、网格独立性验证
网格独立性验证不是“把网格加密一遍看结果变不变”,而是要有明确的加密序列和判据。推荐至少做三套网格:粗、中、细,单元数比例约为 1:2:4。对风冷翅片,监测量选散热器总热阻 Rth 和压降 Δp;对液冷微通道,监测量选平均努塞尔数 Nu 和通道压降 Δp。判据是:相邻两套网格的监测量相对偏差小于 2%,且细网格相对中网格的变化小于中网格相对粗网格的变化,即可认为中网格已进入独立区间。
实际操作中,可用 Richardson 外推估算离散误差。设细网格结果为 f3、中网格为 f2、粗网格为 f1,网格细化比 r=2,则表观收敛阶 p≈ln[(f1-f2)/(f2-f3)]/ln r。若 p 接近模型理论阶数,说明收敛正常;若 p 异常偏小,往往意味着边界层分辨率不足或局部网格畸变,需要回到边界层参数重新调整,而不是继续全局加密。
四、质量指标与常见问题处理
Fluent 中常用的网格质量指标包括正交质量(Orthogonal Quality)、偏斜度(Skewness)和长宽比(Aspect Ratio)。对风冷翅片,正交质量建议大于 0.15,偏斜度小于 0.85;对液冷微通道,由于结构化程度高,正交质量应大于 0.2,偏斜度小于 0.8。若出现大量低质量单元,优先检查翅片顶部圆角、微通道拐角和进出口连接段,这些位置往往因几何倒角过小或扫掠方向冲突而产生畸变。
另一个容易被忽视的问题是边界层与核心区网格的过渡。若增长率过大(如 1.5 以上),边界层最后一层与核心区第一层尺寸突变,会在交界处产生虚假扩散,表现为壁面换热系数偏低或压降偏高。处理方法是把增长率降到 1.2 以下,或在边界层与核心区之间插入一层过渡网格。对于微通道,若单元数受限,可适当减少边界层层数但保持第一层高度不变,因为第一层高度对 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%,即可认为网格满足独立性要求,不必继续加密。这样既控制了计算成本,也保证了风冷与液冷两类散热器在后续性能优化中的结果可信度。