3. 换热器网格划分策略

换热器仿真的网格策略，本质上是在“几何复杂度”与“流动方向性”之间做取舍。管壳式换热器有折流板、进出口接管、管束与壳体间隙，几何高度不规则；板翅式换热器则包含大量薄翅片、隔板和封条，通道狭长且周期性重复。网格划分不是把几何填满就结束，而是要让数值格式在关键区域有足够的解析能力，同时把网格量控制在可计算范围内。

结构化与非结构化网格的选择依据

结构化网格的节点在计算域内按(i, j, k)索引排列，每个内部节点有固定数量的邻居，网格线沿流动方向自然对齐。它的优势是数值扩散小、离散精度高、内存占用低，特别适合管程直管段、板式换热器单通道、套管环隙这类几何规则、流动方向明确的区域。缺点是贴体能力弱，遇到折流板缺口、接管相贯线时难以保持正交性。

非结构化网格以四面体、棱柱或六面体核心加边界层棱柱的形式填充任意几何，贴体性好，适合管壳式换热器整体、带折流板的壳程、复杂进出口段。代价是同等网格量下数值扩散更大，四面体网格在边界层内若不做棱柱层，壁面法向分辨率会严重不足。工程上常见的折中方案是：壳程主体用多面体或六面体核心网格，管束表面和壁面附近用棱柱层，接管和折流板附近局部用四面体过渡。

选择时可参考一个简单判据：若几何中超过70%的壁面可以沿同一方向拉伸成扫掠体，优先考虑结构化或扫掠网格；若折流板、接管、支撑板造成大量非规则相交，则采用非结构化为主、局部结构化的混合策略。对于板翅式换热器，由于翅片通道高度重复，可先对单个周期通道生成结构化网格，再通过周期性复制或扫掠扩展到整个芯体，能显著降低网格量。

边界层网格设置

换热器传热系数对壁面法向速度梯度和温度梯度极其敏感，边界层网格直接决定Nu数预测精度。设置时先确定壁面处理方式：若采用低雷诺数模型（如k-ω SST的增强壁面处理），要求第一层网格中心到壁面的无量纲距离y⁺≈1；若采用标准壁面函数，y⁺应落在30～300之间，且第一层节点应位于对数律层内。

第一层高度可用平板湍流边界层估算：Δy₁ = y⁺·ν / u_τ，其中u_τ = U∞·√(Cf/2)，Cf ≈ 0.058·Re_x^(-0.2)。以壳程水侧为例，U∞=1.2 m/s，特征长度0.02 m，ν=1.0×10⁻⁶ m²/s，Re≈2.4×10⁴，Cf≈0.0058，u_τ≈0.065 m/s。若目标y⁺=1，则Δy₁≈1.5×10⁻⁵ m；若目标y⁺=50，则Δy₁≈7.7×10⁻⁴ m。增长率通常取1.1～1.2，层数15～25层，总厚度应覆盖边界层厚度的1.2倍以上，避免边界层被截断。

管束外壁、折流板两侧、壳体壁面都需要边界层。对于管束，若管间距很小，相邻管的边界层可能重叠，此时应适当减小增长率并增加层数，保证管间速度剖面被充分解析。板翅式通道高度通常只有1～3 mm，边界层可能占满整个通道，此时应保证通道高度方向至少15～20个节点，且第一层高度按y⁺目标反算后检查是否与通道高度匹配。

管束与翅片区域网格加密

管束区域加密的核心是控制管壁周向节点数和管间径向节点数。周向节点数建议不少于40～60个，保证圆柱绕流分离点附近有足够分辨率；管间径向至少8～12个节点，用于捕捉尾流和涡脱落。若管束排列为三角形或旋转正方形，加密时应沿管间最短距离方向布置节点，避免出现长宽比过大的单元。

翅片区域加密要区分翅片厚度方向和通道高度方向。翅片厚度通常为0.1～0.3 mm，若直接按厚度划分，网格量会急剧膨胀。工程上常将翅片简化为零厚度壁面或薄壁热阻模型，只在翅片表面布置边界层，通道高度方向按边界层要求加密。若必须解析翅片厚度，则厚度方向至少3～5层，且与通道方向网格过渡要平滑，避免相邻单元体积突变超过2倍。

局部加密的过渡区应设置增长率控制，一般不超过1.3。加密区与粗网格区之间建议保留2～3层过渡单元，防止数值扩散突然增大。对于周期性翅片通道，可在单个周期内完成加密后通过周期性边界复制，既保证加密一致性，又减少人工操作。

网格质量评价指标与检查方法

Fluent中常用的网格质量指标包括正交质量（Orthogonal Quality）、偏斜度（Skewness）、长宽比（Aspect Ratio）和最小/最大体积。正交质量取值范围0～1，越接近1越好，一般要求最小值大于0.15，关键传热面附近建议大于0.3。偏斜度0～1，越接近0越好，一般要求最大值小于0.85，边界层内可放宽至0.9。长宽比在边界层内允许较大，但核心区建议小于100，否则离散误差会明显增大。

检查方法分三步：第一，在Fluent Meshing或ANSYS Meshing中查看质量直方图，确认没有负体积和零体积单元；第二，用“Mesh Check”检查网格连通性和边界条件映射，确认没有孤立单元和未映射面；第三，在关键截面（如管间最小截面、翅片通道入口）做网格无关性验证，逐步加密至传热系数和压降变化小于2%～3%，即可认为网格收敛。

实际工程中，网格策略往往需要迭代调整。建议先按上述原则生成初始网格，运行100～200步获得初步流场，再根据y⁺云图、温度梯度和速度梯度分布，对梯度大的区域做自适应加密或手动局部加密。每次加密后重新检查质量指标，确保没有因加密引入新的畸形单元。最终网格应在质量、精度和计算成本之间取得平衡，而不是一味追求最密网格。