边界层网格、首层高度与 y+ 设计
说明壁面摩擦速度、无量纲壁面距离 y+、首层高度、棱柱层数量和增长率之间的关系,并区分壁函数与低雷诺近壁解析策略。
边界层网格y+首层高度壁函数棱柱层工程
边界层网格、首层高度与 y+ 设计
近壁速度和温度梯度决定摩擦阻力与换热。首层高度必须与湍流模型的近壁处理相匹配,不能在计算后才把 y+ 当作装饰性云图。
1. y+ 定义
y+ = ρ uτ y / μ
uτ = sqrt(τw / ρ)
其中 y 是壁面到首个计算点的法向距离。设计阶段可用摩擦因子关联式估算壁面剪切和 uτ,计算后再用实际壁面剪切校核。
2. 两类近壁策略
低雷诺解析通常要求首点位于黏性底层,常以 y+ 约 1 为目标;传统壁函数把首点放在对数区,具体建议范围取决于软件和壁函数实现。应避免大量单元落在既不适合壁函数、又未解析黏性底层的过渡区。
3. 棱柱层设计
除了首层高度,还要保证总厚度覆盖主要边界层、层数足以表达梯度、增长率平滑,并在尖角和间隙处避免层塌陷。热边界层与速度边界层厚度可能不同,需结合 Pr 数判断。
4. 计算后检查
统计目标壁面的 y+ 最小值、最大值、面积分布和异常位置。若模型采用壁函数但大部分面积 y+ 很低,应重新选择近壁处理或网格;若局部高 y+ 位于分离、冲击或高热流区,应优先修复。
5. 常见错误
只给一个平均 y+、用入口速度估算所有壁面、棱柱层总厚度过小,以及在几何尖角处保留极差单元,都会降低结果可信度。
6. 参考资料
- Schlichting & Gersten, Boundary-Layer Theory.
- Wilcox, Turbulence Modeling for CFD.