在计算流体力学（CFD）和复杂几何模拟中，处理运动部件或几何形状变化较大的场景时，传统单一网格方法常面临网格扭曲、重构耗时等挑战。重叠网格技术通过将计算域划分为多个独立网格区域，使各区域能够相对运动或局部加密，为这类问题提供了高效解决方案。本文将系统介绍重叠网格的核心原理、实施步骤及典型应用场景，通过完整案例演示其配置流程，并针对常见实施难点提供排查思路。

1. 重叠网格技术概述

1.1 基本概念与定义

重叠网格（Overset Grid），又称嵌套网格或复合网格，是一种将多个独立网格区域重叠组合的离散化方法。其核心思想是将复杂计算域分解为多个简单子域，每个子域生成独立网格，通过插值交换边界信息实现全场耦合求解。

与传统单一网格相比，重叠网格具备三大优势：

• 几何适应性：各子网格可针对局部几何特征优化（如边界层网格、运动部件网格）
• 运动处理：支持网格间相对运动，避免动态网格重构
• 并行效率：子网格可独立划分计算任务，提升并行计算规模

1.2 典型应用场景

重叠网格技术特别适用于以下三类问题：

1. 运动边界问题：如直升机旋翼、风力机叶片、船舶螺旋桨等旋转机械
2. 多尺度几何：如飞机襟翼缝翼、汽车后视镜等局部精细结构
3. 自适应加密：在激波、剪切层等物理现象显著区域局部加密网格

2. 重叠网格实施基础

2.1 核心组件与术语

实施重叠网格需理解以下关键组件：

• 背景网格：覆盖主要计算域的基础网格，通常较稀疏
• 组件网格：围绕运动部件或局部区域的精细网格
• 洞切割：标识被组件网格覆盖的背景网格区域，这些区域不参与计算
• 插值边界：组件网格与背景网格的数据交换界面

2.2 数据交换原理

重叠网格的核心是组件网格与背景网格间的数据交换，通过以下步骤实现：

1. 洞识别：自动识别被组件网格覆盖的背景网格单元
2. 插值点定位：为每个组件网格边界点定位其在背景网格中的宿主单元
3. 插值计算：基于宿主单元节点值通过权重插值获得边界值
4. 数据传递：将插值结果作为组件网格的边界条件

3. 重叠网格生成流程

3.1 几何准备与分解

首先需将计算几何分解为适合独立网格生成的子区域：

# 示例：几何分解策略 geometry_decomposition = { "background_region": "远场流域", "component_regions": [ {"name": "rotor", "type": "旋转部件", "refinement_level": 2}, {"name": "wake", "type": "尾迹区域", "refinement_level": 1} ] }

分解原则：各子区域间应有足够重叠区域（通常3-5层网格）确保插值稳定性。

3.2 独立网格生成

为每个子区域生成质量优化的独立网格：

# 背景网格生成命令示例 blockMesh -case background_region snappyHexMesh -case background_region -overwrite # 组件网格生成命令示例 surfaceFeatureExtract -case rotor_region snappyHexMesh -case rotor_region -overwrite

关键参数：重叠区域网格尺寸应渐变过渡，避免尺寸跳跃过大。

3.3 重叠关系建立

通过网格组装工具建立重叠关系：

# 重叠关系配置示例 overset_config = { "background_mesh": "background/constant/polyMesh", "component_meshes": [ "rotor/constant/polyMesh", "wake/constant/polyMesh" ], "hole_cut_method": "distance", "interpolation_tolerance": 1e-6 }

4. OpenFOAM中的重叠网格实施

4.1 环境准备与求解器选择

OpenFOAM提供了成熟的重叠网格支持，推荐环境：

• OpenFOAM v1912及以上版本
• 支持MPI的并行环境
• 适用于重叠网格的求解器：overPimpleDyMFoam（可压缩流）、overInterDyMFoam（不可压缩流）

4.2 基础目录结构

典型的重叠网格案例目录结构：

overset_case/ ├── constant/ │ ├── dynamicMeshDict # 动网格配置 │ └── oversetMeshDict # 重叠网格配置 ├── system/ │ ├── controlDict # 求解控制 │ ├── fvSchemes # 离散格式 │ └── fvSolution # 求解算法 └── 0/ # 初始条件

4.3 关键配置文件详解

dynamicMeshDict 配置：

dynamicFvMesh oversetFvMesh; oversetFvMeshCoeffs { // 背景网格定义 backgroundMesh { meshName background; cellZone backgroundZone; } // 组件网格定义 componentMeshes ( { meshName rotor; cellZone rotorZone; motion rotating; } ); }

oversetMeshDict 配置：

oversetInterpolation { method inverseDistance; tolerance 1e-6; nIterations 100; } holeCutting { method cellVolumeWeight; exposedPatchName overset; }

4.4 边界条件设置

重叠网格区域的边界条件需要特殊处理：

// 0/p 文件示例 boundaryField { overset { type overset; value uniform 0; } // 其他边界条件 inlet { type fixedValue; value uniform 1e5; } }

5. 旋转机械案例实战

5.1 案例背景与几何设置

以风力机叶片旋转为例，建立包含以下区域的网格系统：

• 背景网格：20R×20R×10R的矩形流域（R为叶片半径）
• 旋转网格：围绕叶片的圆柱形区域，半径1.2R，高度0.5R
• 尾迹网格：下游延伸区域，捕捉尾迹发展

5.2 网格生成与组装

# 步骤1：生成背景网格 blockMesh -case background # 步骤2：生成旋转区域网格 surfaceFeatureExtract -case rotor snappyHexMesh -case rotor -overwrite # 步骤3：组装重叠网格 createOversetMesh -case overset_case

5.3 运动定义与求解控制

旋转运动通过dynamicMeshDict定义：

rotorZone { type rotating; origin (0 0 0); axis (0 0 1); omega 10; // 转速10 rad/s }

求解器设置需特别注意重叠网格兼容性：

// system/fvSolution 配置 solvers { "(U|p|k|omega)" { solver GAMG; tolerance 1e-6; relTol 0.1; } } oversetInterpolation { correctPhi yes; massCorrection yes; }

5.4 计算结果验证

成功实施后应检查以下指标：

• 插值边界质量：插值权重应平滑分布
• 质量守恒：全局质量误差应小于1e-10
• 运动同步：各网格区域运动协调无冲突

6. 常见问题与解决方案

6.1 插值不稳定问题

现象：计算发散或出现非物理振荡原因：

• 重叠区域网格尺寸跳跃过大
• 插值点定位精度不足
• 时间步长设置不合理

解决方案：

// 优化插值参数 oversetInterpolation { method inverseDistance; tolerance 1e-8; // 提高定位精度 nIterations 200; // 增加迭代次数 acceptors all; // 接受所有插值点 } // 调整时间步长控制 controlDict { maxCo 0.5; // 降低库朗数 maxDeltaT 1e-4; // 限制最大时间步 }

6.2 质量守恒误差

现象：计算域内质量不守恒原因：

• 洞边界识别错误
• 插值边界通量校正不足

解决方案：

// 启用质量校正 oversetFvMeshCoeffs { massCorrection yes; correctPhi yes; } // 检查洞切割质量 holeCutting { method cellVolumeWeight; minVolumeRatio 0.01; // 提高切割精度 }

6.3 并行计算负载不均

现象：并行效率低下，部分进程空闲原因：网格分区未考虑重叠网格特性

解决方案：

# 使用加权分区 decomposePar -weighted -overwrite

7. 性能优化最佳实践

7.1 网格质量优化

• 尺寸过渡：重叠区域网格尺寸比控制在1.5以内
• 正交性：保持插值边界网格正交性>25°
• 层数保证：重叠区域至少保证3层网格单元

7.2 计算参数调优

// 优化求解策略 fvSolution { relaxationFactors { equations { U 0.7; p 0.3; } } // 使用GAMG预处理提高收敛性 p { solver GAMG; preconditioner DIC; tolerance 1e-6; relTol 0.1; } }

7.3 运动定义优化

对于复杂运动，建议使用函数定义：

rotorMotion { type coded; name prescribedMotion; codeInclude { #include "fvCFD.H" }; codeOptions { -I$(LIB_SRC)/finiteVolume/lnInclude }; codeLibs { -lfiniteVolume }; code { // 自定义运动函数 vector displacement = vector(0, 0, sin(constant::mathematical::pi*time()/10)); return displacement; }; }

8. 工程应用注意事项

8.1 网格独立性验证

实施重叠网格后必须进行网格独立性验证：

1. 逐步加密各区域网格，观察关键参数变化
2. 确保结果变化小于5%时认为网格独立
3. 特别注意重叠区域的网格敏感性

8.2 计算资源规划

重叠网格计算资源需求估算：

• 内存需求：各网格区域内存之和的1.2-1.5倍（考虑重叠开销）
• 计算核心：建议每个网格区域分配独立计算节点
• 存储需求：结果文件大小为传统网格的1.3-2.0倍

8.3 结果后处理技巧

重叠网格结果后处理需要特殊考虑：

# ParaView中合并显示多个网格区域 reader1 = OpenFOAMReader(FileName='background.foam') reader2 = OpenFOAMReader(FileName='rotor.foam') # 使用Append Datasets合并显示 append = AppendDatasets(Input=[reader1, reader2])

重叠网格技术为复杂流动模拟提供了强大工具，但成功实施需要系统的方法和细致的参数调试。建议从简单几何开始逐步积累经验，重点关注网格质量、插值稳定性和运动协调性三个关键环节。