📚 课程目录第4章计算域与网格划分策略

4. 计算域与网格划分策略

微通道冷却仿真的结果可信度，很大程度上在网格阶段就已经被决定。微通道的典型几何特征是：通道水力直径常在 100～1000 μm 量级，长径比可达 50～200，固体与流体域紧密耦合，壁面附近存在速度与温度梯度极大的边界层。这意味着网格策略必须同时满足三个要求：通道内部轴向流动能被低耗散地解析、近壁面法向梯度能被足够密的层状网格捕捉、固体导热与流固界面热流能被连续传递。本节围绕结构化网格、非结构化网格、边界层网格与多面体网格四类方案，从几何适配、网格质量、计算代价与结果可信度四个维度展开对比。

四类网格的本质差异

结构化网格指每个内部节点在拓扑上具有相同数量的相邻单元，通常表现为四边形（2D）或六面体（3D）的规则排列。在微通道中，若通道为直矩形或直圆形，可沿轴向、展向、法向建立 i-j-k 索引，生成高质量六面体网格。其优势是数值耗散低、轴向流动方向与网格线对齐、可用较少的单元获得较高的精度，且便于做网格无关性验证。局限在于对复杂几何（弯折、歧管、针肋、渐变截面）适配性差，往往需要分块拼接，块间界面处理不当会引入插值误差。

非结构化网格以四面体、棱柱或混合单元为主，节点连接关系不固定，能自动填充任意复杂几何。微通道歧管、进出口集流腔、针肋阵列等区域用四面体网格可快速生成。代价是相同单元数下精度通常低于六面体，且四面体在近壁面若不做棱柱层，会因单元形状与流动方向不匹配而严重高估壁面剪切与换热。

边界层网格不是独立网格类型，而是针对近壁区的加密策略：在壁面法向堆叠若干层棱柱或六面体单元，第一层高度由目标 y+ 反推，层间以 1.1～1.3 的增长率向外扩展。微通道中边界层厚度与通道半径同量级，若第一层过高，会把整个通道“吃掉”，因此边界层网格必须与通道核心区网格协同设计。

多面体网格由多面体单元（常见为十二面体、十四面体）构成，每个单元有更多相邻面，梯度重构更充分。工程经验表明，在相同精度下多面体网格单元数约为四面体的 1/3～1/2，收敛更快，对复杂几何适应性强。但在微通道这类强方向性流动中，多面体网格的随机取向会削弱轴向解析能力，通常需要与边界层棱柱层配合使用。

微通道中的关键参数：y+ 与第一层高度

边界层网格的核心控制量是壁面第一层网格的无量纲高度 y+。定义 y+ = ρ·u_τ·y/μ，其中 u_τ = √(τ_w/ρ) 为摩擦速度，y 为第一层单元中心到壁面的距离。若采用低雷诺数模型或增强壁面处理，通常要求 y+ ≈ 1；若采用壁面函数，则要求 30 < y+ < 300。微通道中流速较高、通道尺寸小，u_τ 往往较大，因此第一层高度常在 0.1～1 μm 量级。

以水为例，取通道水力直径 D_h = 500 μm，入口平均流速 2 m/s，20 ℃ 水物性 ρ = 998 kg/m³、μ = 1.002×10⁻³ Pa·s，则 Re ≈ 998×2×500×10⁻⁶/1.002×10⁻³ ≈ 996，属层流。层流下壁面剪切可由解析式估算，若取 τ_w ≈ 8 μU/D_h ≈ 8×1.002×10⁻³×2/500×10⁻⁶ ≈ 32 Pa，则 u_τ = √(32/998) ≈ 0.179 m/s，y+ = 1 对应 y ≈ μ/(ρ·u_τ) ≈ 1.002×10⁻³/(998×0.179) ≈ 5.6×10⁻⁶ m，即约 5.6 μm。若目标 y+ = 1，第一层高度应取约 2.8 μm（单元中心距壁面 2.8 μm，单元高度约 5.6 μm）。这一量级说明：微通道边界层网格必须精细到微米级，否则近壁梯度会被严重抹平。

层流微通道中，若采用解析壁面剪切估算 u_τ 有困难，可先用较粗网格试算一轮，读取壁面剪切应力分布，再反推 u_τ 并重设第一层高度。这是工程上最稳妥的迭代做法。

四类网格在微通道不同区域的组合策略

实际微通道冷却系统很少只用一种网格。合理的做法是按区域分工：

直通道段：优先采用结构化六面体，沿轴向拉伸，展向与法向按边界层要求加密。轴向单元尺寸可放宽至 20～50 μm，法向第一层控制在 2～5 μm，增长率 1.15～1.2，层数 10～15。
进出口歧管与集流腔：几何复杂、流动三维性强，采用非结构化四面体或六面体核心+四面体过渡，单元尺寸 30～80 μm，局部在拐角与射流冲击区加密。
针肋/扰流柱阵列：柱面附近必须布置边界层棱柱层，柱后尾流区可用多面体或四面体，避免因网格取向与分离流方向不一致导致尾流解析失真。
固体导热域：若只关心温度场，可用较粗的六面体或四面体，但流固界面两侧网格尺寸比不宜超过 3:1，否则界面热流插值误差会明显增大。

网格无关性与质量判据

微通道仿真必须做网格无关性验证，建议至少三套网格：粗、中、细，单元数按约 1.5～2 倍递增。观察量优先选压降 Δp、壁面平均温度 T_w,avg、热阻 R_th 和最大壁温 T_w,max。当细网格相对中网格的上述指标变化小于 2%～3%，且中网格相对粗网格变化明显更大时，可认为中网格已进入无关区。

质量判据方面：六面体网格正交质量（Orthogonal Quality）应大于 0.2，最好大于 0.3；四面体网格偏斜度（Skewness）应小于 0.85；边界层棱柱层要检查是否存在负体积或高度突变。微通道中常见的问题是边界层与核心区过渡处增长率过大，导致相邻层高度比超过 1.5，此时应降低增长率或增加层数。

不要为了追求 y+ = 1 而把第一层高度设得过小。若第一层高度远小于分子尺度或求解器双精度极限，会导致矩阵条件数恶化、收敛困难。微通道中第一层高度一般不低于 0.5 μm，若计算得到的 y+ 仍偏大，应优先检查流速与物性输入是否合理。

选型建议汇总

综合来看：直矩形/圆形微通道、以压降与平均换热为目标时，结构化六面体+边界层棱柱层是精度与效率的最佳组合；含歧管、针肋、弯折的复杂微通道，采用多面体+边界层棱柱层或六面体核心+四面体过渡，可在可接受的计算量下获得稳定结果；纯四面体网格仅在几何极度复杂且对近壁精度要求不高时使用，且必须配合棱柱层。无论选择哪种方案，边界层网格在微通道中都不是可选项，而是决定换热与压降预测是否可信的关键环节。