Hypermesh体网格生成时如何控制边界层质量？

Question

在使用HyperMesh生成体网格时，边界层质量控制是影响CFD仿真精度的关键环节。常见的技术问题是：如何在复杂几何表面（如曲率变化大或狭小间隙区域）合理设置边界层参数，以避免出现扭曲、反转或高纵横比的棱柱层？用户常遇到首层厚度与增长比不匹配、棱柱层数不足或穿透核心流体区等问题，导致Y+值不达标或网格畸变。需结合流动特征、壁面分辨率要求及几何特点，通过调整“First Layer Height”、“Growth Rate”和“Number of Layers”等参数，并借助检查工具如Jacobian、skewness进行质量评估，确保边界层网格平滑过渡且满足求解器输入要求。

Answer 1

1. 边界层网格生成的基本概念与关键参数

在使用HyperMesh进行CFD体网格划分时，边界层（Prism Layer）的生成是确保壁面附近流动解析精度的核心环节。边界层网格通常由一系列沿壁面法向生长的棱柱单元构成，其质量直接影响Y+值的准确性及湍流模型的适用性。

• First Layer Height：首层高度决定了最靠近壁面的网格节点位置，需根据目标Y+值反推计算。
• Growth Rate：增长比控制每一层相对于前一层的厚度扩展速度，典型值为1.1~1.3。
• Number of Layers：层数影响边界层总厚度覆盖范围，应足以捕捉速度梯度变化。

这些参数的选择必须结合雷诺数、流动类型（如分离流、附着流）以及几何特征复杂度综合判断。

2. 常见技术问题分析

在实际操作中，用户常因参数设置不当导致以下典型问题：

1. 首层过厚或过薄，造成Y+偏离理想区间（如k-ε模型推荐30~300）；
2. 增长比过大引发棱柱层扭曲甚至反转；
3. 层数不足无法充分解析边界层发展；
4. 在狭小间隙区域发生棱柱穿透至对面壁面；
5. 高曲率表面出现网格折叠或Jacobian劣化；
6. 多部件交界处棱柱冲突或重叠；
7. 自动尺寸函数未考虑局部几何细节，导致分辨率不均；
8. 未能同步调整全局网格尺寸与边界层匹配；
9. 忽略近壁面曲率对法向方向的影响；
10. 缺乏后处理验证手段，难以定位质量问题根源。

3. 参数设定方法与工程实践建议

参数	推荐范围	影响因素	调整策略
First Layer Height	5×10⁻⁶ ~ 5×10⁻⁴ m	Re, μ, U∞, Y+	通过Y+估算公式反算
Growth Rate	1.1 ~ 1.3	曲率、压力梯度	高曲率区取低值
Number of Layers	5 ~ 20	边界层厚度δ	确保覆盖90%以上速度剖面
Total Prism Thickness	< δ	Re, 物体特征长度	避免侵入主流区

4. 质量评估指标与检查流程

完成边界层生成后，必须借助HyperMesh内置的质量检查工具进行验证：

// 示例：在Tcl脚本中批量输出棱柱层质量统计
*createmark elems 1 "boundary_layer_prisms"
set jacobian_min [hm_getentityvalue elems 1 Jacobian 6]
set skewness_max [hm_getentityvalue elems 1 Skewness 6]
puts "Min Jacobian: $jacobian_min | Max Skewness: $skewness_max"

重点关注以下指标：

• Jacobian > 0.6 表示单元变形可控；
• Skewness < 0.9 避免严重偏斜；
• Aspect Ratio 单层内不宜超过20；
• Check for inverted elements (negative volume)；
• Visual inspection of layer continuity and smooth transition.

5. 复杂几何下的高级处理策略

graph TD
    A[导入CAD几何] --> B{是否存在狭小间隙?}
    B -- 是 --> C[启用Small Feature Capture]
    B -- 否 --> D[标准面网格划分]
    C --> E[设置最小边长阈值]
    D --> F[应用Inflation Parameters]
    F --> G{曲率变化大?}
    G -- 是 --> H[降低Growth Rate至1.15]
    G -- 否 --> I[使用默认1.2]
    H --> J[执行Boundary Layer Meshing]
    I --> J
    J --> K[运行Quality Check]
    K --> L{满足标准?}
    L -- 否 --> M[调整First Layer或层数]
    L -- 是 --> N[导出至求解器]

针对高曲率区域，可采用局部坐标系定义法向方向，并利用“Curvature-based”尺寸函数增强分辨率。对于存在多个相邻边界层的情况，建议启用“Collision Detection”功能防止穿透。