深入拆解OpenFOAM的pitzDaily案例：从网格划分到湍流模型实战解析

当你第一次成功运行OpenFOAM的pitzDaily案例时，那种成就感可能让你兴奋不已。但很快你会发现，仅仅会复制粘贴运行案例远远不够——真正的价值在于理解每个参数背后的工程逻辑。本文将带你深入这个经典案例的底层设计，让你从"会运行"进阶到"会设计"。

1. pitzDaily案例的工程背景与几何简化

pitzDaily是一个经典的二维后向台阶流动案例，得名于研究者Daily在Pitz实验室的工作。这个看似简单的几何形状，却包含了分离流、再附着流、剪切层发展等丰富的流体力学现象。

为什么选择这个几何？ 在工程实践中，后向台阶流动常见于管道扩张、燃烧室入口等场景。pitzDaily案例的几何尺寸（台阶高度H=12.7mm，入口高度2H，出口高度3H）经过精心设计，使得雷诺数在合理范围内（约10^4量级），既能展现湍流特性又不会过于复杂。

几何定义在blockMeshDict中体现为6个关键顶点：

cpp复制vertices
(
    (0 0 0)       // 0
    (0 0.0127 0)  // 1
    (0.0508 0 0)  // 2
    (0.0508 0.0127 0) // 3
    (0.3048 0 0)  // 4
    (0.3048 0.0381 0) // 5
);

这种参数化设计让案例既保留了物理本质，又避免了不必要的计算开销。实际工程中，我们常需要做类似的简化：

• 将三维问题降维处理（如对称性利用）
• 合理截断计算域（如出口位置选择）
• 特征尺寸标准化（以台阶高度为基准）

提示：修改几何时，保持特征尺寸的比例关系比绝对尺寸更重要。例如改变台阶高度时，应同步调整其他尺寸。

2. 网格生成策略与blockMeshDict解析

pitzDaily的网格设计体现了对流动特性的预判。打开constant/polyMesh/blockMeshDict，你会发现它采用了分块结构化网格：

cpp复制blocks
(
    hex (0 1 3 2 0 1 3 2) (20 10 1) simpleGrading (1 1 1) // 入口段
    hex (2