简介OpenFOAM 是开源的计算流体力学工具自带算例是理解各类求解器最直观的入口。这份 PDF 面向刚接触 OpenFOAM 的工程师、科研人员和学生从中挑选了摩托车绕流、法兰冷热交接面热扩散、圆柱绕流、非牛顿流体流动、管道流、化学反应与甲烷燃烧、发动机内燃、顶盖驱动流、搅拌器、管内流、混合器、弯管流动等十余类典型算例逐一对应 simpleFoam、laplacianFoam、potentialFoam、nonNewtonianIcoFoam、chemFoam、engineFoam、fireFoam、icoFoam、pimpleFoam 等常用求解器并结合物理场景介绍求解思路帮助读者快速建立“问题—算例—求解器”的映射关系。资源为 1 个 PDF 文件大小约 1.73MB形式轻量、结构清晰除算例解析外还附有标准求解器分类可作为 OpenFOAM 入门阶段的速查手册。已有 427 人学习下载适合希望系统梳理自带算例、快速定位所需求解器并用于后续模拟的 CFD 学习者。1. 从 tutorials 目录读 OpenFOAM 算例的两种方式很多初学者把自带算例当成能跑就行的演示包跑完不改参数等到自己搭算例时才发现处处报错。实际上 tutorials 目录是 OpenFOAM 最值得反复读的配置档案每个子文件夹都对应一条完整的求解链路从网格生成、离散格式到收敛判据全部可以一键复现。读这类算例的姿势有两种按求解器横向读对比同一几何在不同算法下的行为差异或者按物理问题纵向读找出一个算例族里求解器、边界条件和网格要求是怎么配套的。这篇内容按不可压缩流、燃烧传热、初始场生成、算例验收四条线展开适合刚装完 OpenFOAM 想系统过一遍求解器的人也适合要迁移到新求解器但不想从空白 case 起步的熟手。跑任何算例之前先确认环境变量指向正确执行echo $FOAM_TUTORIALS能看到自带算例根目录即可后面所有路径都基于这个目录。下面先从不可压缩流这个最大分支开始拆。2. 不可压缩流求解器从 motorBike 到 cavity 的算法对照2.1 motorBike 的完整启动链路与 SIMPLE 的稳态含义motorBike 算例位于 tutorials/incompressible/simpleFoam/motorBike演示的是外部绕流的标准启动流程先建立矩形背景网格再通过 snappyHexMesh 把摩托车表面贴合到网格上最后交给 simpleFoam 迭代。手动执行的命令序列如下cd $FOAM_TUTORIALS/incompressible/simpleFoam/motorBike blockMesh surfaceFeatureExtract -inputDict system/surfaceFeatureExtractDict snappyHexMesh -overwrite simpleFoam这段链路里blockMesh 负责生成带加密区的背景网格surfaceFeatureExtract 从 STL 表面提取特征线供 snappyHexMesh 做表面贴合-overwrite表示直接覆盖 constant/polyMesh避免生成多层时间目录最后的 simpleFoam 做稳态迭代。snappyHexMesh 的 log 里如果出现大量负体积警告不要直接进入 simpleFoam先回 blockMeshDict 加宽加密区否则速度场里会留下不光滑的云图这类网格问题调求解器参数救不回来。simpleFoam 使用 SIMPLE 算法稳态求解不可压缩湍流核心是压力修正配合亚松弛迭代没有真实时间步。fvSolution 里的 relaxationFactors 控制收敛路径U 通常设 0.7p 设 0.3k 和 epsilon 设 0.5 左右。残差呈锯齿状振荡时先降 U 的松弛因子而不是加密网格。网格密度在 system/blockMeshDict 里调整每个 block 的子网格数就藏在 blocks 列表的最后三个数里blocks ( hex (0 1 2 3 4 5 6 7) (120 60 40) simpleGrading (1 1 1) );括号里前 8 个点是角点编号(120 60 40)是 x/y/z 三个方向的网格数simpleGrading 后面是各方向膨胀比。改成(180 90 60)就能得到更密的网格但注意外流算例远离壁面的区域加密收益很低优先加密车体附近的 block 更划算。2.2 用 cavityHighRe 对比 icoFoam 与 pimpleFoam 的时间步策略顶盖驱动流算例 cavityHighRe 是 icoFoam 的标配瞬态求解不可压缩牛顿流体层流。运行方式比 motorBike 简单得多cd $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavityHighRe blockMesh icoFoamicoFoam 基于 PISO 算法每个时间步内做多次压力修正时间步长受库朗数限制Co 必须小于 1。算例 controlDict 里 deltaT 取 0.005是根据 Co U*deltaT/deltaX 反推出来的改网格后必须重新估算这个值否则很容易出现压力震荡。与之对照的 channel395 算例使用 pimpleFoamPIMPLE 把 PISO 的内循环和 SIMPLE 的外循环合并允许 Co 远大于 1适合 channel395 这类需要长时间统计平均的高雷诺数槽道湍流。关键参数在 fvSolution 里PIMPLE { momentumPredictor no; nOuterCorrectors 3; nCorrectors 2; nNonOrthogonalCorrectors 1; }momentumPredictor 为 no 表示第一个外循环不求解动量预演nOuterCorrectors 是每个时间步的外循环次数值越大越稳但越贵nCorrectors 是每个外循环内的 PISO 压力修正次数nNonOrthogonalCorrectors 处理非正交网格修正通常 1 就够。大时间步算例发散时先降 deltaT再调大 nOuterCorrectors不要一上来就换湍流模型。求解器算法算例路径典型场景时间步限制simpleFoamSIMPLEmotorBike外流稳态无真实时间步icoFoamPISOcavityHighRe顶盖驱动流Co 1pimpleFoamPIMPLEchannel395高Re槽道流Co 可大于1nonNewtonianIcoFoamPISOoffsetCylinder非牛顿流体层流Co 1SRFPimpleFoamPIMPLE SRFrotor2D旋转参考系由旋转速度约束SRFSimpleFoamSIMPLE SRFmixer旋转稳态混合无真实时间步2.3 offsetCylinder 与 rotor2D非牛顿和旋转参考系的两个典型坑nonNewtonianIcoFoam 的算例 offsetCylinder 是偏心圆柱间的聚合物流体流动transportProperties 里要指定非牛顿本构模型常见的是 Bird-CarreautransportModel BirdCarreau; nu0 1e-04; nuInf 1e-06; k 0.01; n 0.4;nu0 是零剪切黏度nuInf 是极限黏度k 是时间常数n 是幂律指数。n 小于 1 就是剪切变稀计算收敛后不要只看残差还要看壁面附近的 nu 分布是否光滑突变点通常意味着边界层第一层网格高度没匹配好。rotor2D 属于 SRFPimpleFoam用旋转参考系把旋转机械问题转成壁面旋转坐标省去动网格。常见错误是把进出口边界写成 absolute 速度而旋转区内部是 relative叠加后速度场和压力场错位。收敛后在 ParaView 里看转子进出口的轴向速度剖面如果出现锯齿优先检查 U 的边界类型是 rotatingWallVelocity 还是 fixedValue。3. 传热与燃烧算例族flange、chemFoam 与 engineFoam 的边界3.1 flange 算例与拉普拉斯热扩散方程的参数迁移laplacianFoam 是 OpenFOAM 里代码量最小的求解器之一核心方程就是瞬态热扩散solve(fvm::ddt(T) - fvm::laplacian(DT, T));flange 算例模拟法兰冷热交接面模型左侧高温、右侧低温DT 是热扩散系数T 是温度场。这个算例的价值在于它足够短适合第一次改源码的人把 DT 从 uniform 改成随温度变化的多项式只需要在 createFields.H 里替换 DT 的初始化方式重编译后跑同一个算例对比温度场差异整个流程十分钟内能完成。运行命令是cd $FOAM_TUTORIALS/basic/laplacianFoam/flange blockMesh laplacianFoam稳态后计算温度梯度postProcess -func grad(T)grad(T) 会输出温度梯度场梯度最大的区域就是热流集中的位置。观察这个场可以帮你在自定义固体传热算例时判断网格是否需要局部加密而不是盲目全局加密。与之相关的还有 rhoPimpleDyMFoam 的 annularThermalMixer 算例那是 HVAC 场景下的动网格传热比静态 flange 多了一层网格运动与温度场的耦合跑之前先确认动网格边界在 dynamicMeshDict 里的设置和初始时间步匹配。3.2 chemFoam 与 GRI-Mech 3.0 的零维反应器对比chemFoam 是单单元化学求解器算例路径 tutorials/basic/chemFoam/gri用的机理是 GRI-Mech 3.0包含 53 个组分、325 个反应模拟甲烷燃烧。所谓零维反应器是计算域只有一个单元没有流动只解化学反应源项引起的组分和温度变化。这个设计让它非常适合与 CHEMKIN 这类专用化学动力学工具做交叉验证。运行方式cd $FOAM_TUTORIALS/basic/chemFoam/gri chemFoamchemistryProperties 里可以选择 ODE 积分器和容差对点火延迟影响最大的是绝对容差和相对容差默认值可能偏松导致点火延迟偏大。把容差收紧一个数量级看温度时间曲线是否明显变化就能判断当前结果是数值误差还是机理本身的属性。GRI 文件里反应速率常数常见单位是 cm3/mol/s和 OpenFOAM 内部单位制不一致时需要在导入机理时做单位换算这是 chemFoam 算例最容易出错的地方。求解器算例路径物理场景关键特征chemFoambasic/chemFoam/gri零维点火53组分325反应engineFoamcombustion/engineFoam/kivaTest内燃机燃烧活塞动网格fireFoamcombustion/fireFoam/les/smallPoolFire2D油池火大涡模拟XiFoamcombustion/XiFoam/ras/moriyoshiHomogeneous预混火焰湍流燃烧速度模型3.3 燃烧求解器的场景边界与可压缩算例的注意事项engineFoam 对应内燃机燃烧kivaTest 算例带活塞运动属于拓扑变化动网格fireFoam 对应开放空间扩散火焰smallPoolFire2D 是二维小油池火适合火灾模拟XiFoam 处理可压缩预混/部分预混湍流燃烧moriyoshiHomogeneous 是均匀预混气点火问题。区分这三个求解器的方法很直接有活塞运动选 engineFoam开放空间扩散火选 fireFoam预混火焰传播选 XiFoam。容易踩的坑是把 smallPoolFire2D 当成二维轴对称算例用它实际是笛卡尔网格需要轴对称几何时必须自己把 blockMeshDict 改成 wedge 结构。可压缩分支里rhoSimplecFoam/squareBend 是跨音速弯管sonicLiquidFoam/decompressionTankFine 是液体减压柜。这两个算例涉及激波和压力波传播跑之前要确认 thermo 字典里的状态方程与算例原始设置一致比如 perfectGas 还是 Peng-Robinson状态方程选错压力波传播速度会差一个量级表现出的振荡频率完全不对。4. potentialFoam 与 boundaryFoam初始场和入口边界条件的生成链路4.1 potentialFoam 为什么只需要解压力方程potentialFoam 是势流求解器读入标量 p、向量 U 和表面场 phi求解压力方程后根据表面场重建速度场。cylinder 算例路径是 tutorials/basic/potentialFoam/cylinder典型的圆柱绕流。势流假定流动无旋、无黏所以不需要解完整的 NS 方程只需要满足连续性计算代价极低。用势流场作为 NS 求解的初始场能避免从零场启动时压力-速度耦合带来的发散。运行方式和解普通算例一样先 blockMesh 再 potentialFoamcd $FOAM_TUTORIALS/basic/potentialFoam/cylinder blockMesh potentialFoam foamToVTK -latestTimefoamToVTK 导出后在 ParaView 里检查 U 场前驻点速度应为零后驻点应出现对称的绕流形态。如果势流场后缘出现过大的负压尖峰说明圆柱附近的网格分辨率不够直接在 blockMeshDict 里加密过渡区。边界条件的搭配是关键势流计算要求边界上 p 和 U 的类型配对正确典型的配对方式如下边界压力 p速度 U说明inletzeroGradientfixedValue给定来流速度outletfixedValuezeroGradient给定出口压力wallzeroGradientnoSlip无滑移壁面frontAndBackemptyempty二维算例必须保留自定义几何时照抄这种配对否则压力方程会出现奇异残差直接卡在初始量级不下降。4.2 boundaryFoam 生成 inlet 剖面再映射回主算例boundaryFoam 主要用于一维稳态不可压缩湍流模型求解算例路径是 tutorials/incompressible/boundaryFoam/boundaryLaunderSharma。它的特殊价值在于生成一个充分发展的 inlet 边界条件供后续算例使用。对于 LES 或 高质量 RANS 算例入口条件直接用 uniform 速度会导致转换延迟甚至层流化而 boundaryFoam 可以算出包含湍流边界层信息的入口剖面。cd $FOAM_TUTORIALS/incompressible/boundaryFoam/boundaryLaunderSharma blockMesh boundaryFoam foamToVTK -latestTime导出后重点看 U 沿壁面法向的剖面验证是否接近充分发展的湍流边界层对数律。LaunderSharma 是低雷诺数 k-epsilon 模型壁面第一层网格 y 必须小于 1否则生成的剖面是错的。确认剖面合理后用 mapFields 把它映射到目标算例的 inlet 边界这是把边界层信息迁移过去的常见做法。映射完成后还要检查一次入口附近的 U 分布是否连续因为两个算例的网格拓扑不同映射过程会产生插值误差。4.3 pitzDaily 把标量输运和伴随优化串起来pitzDaily 是后向台阶管道流同一个几何同时被 scalarTransportFoam 和 adjointShapeOptimizationFoam 使用。scalarTransportFoam 求解标量输运方程控制方程是solve(fvm::ddt(T) fvm::div(phi, T) - fvm::laplacian(DT, T));这个算例里 T 是被动标量只随给定速度场做对流和扩散不反馈影响流场。要得到有意义的标量输运结果必须先有收敛的速度场常见做法是用 simpleFoam 先算出稳态流场再切到 scalarTransportFoam 继续迭代。如果直接拿 uniform 速度场跑输运结果毫无参考价值。adjointShapeOptimizationFoam 走的是伴随法形状优化路线稳态求解不可压缩非牛顿流体在变形管道中的受阻湍流流动输出目标函数对形状的敏感度场。敏感度场正值表示该区域收缩可降低目标函数负值表示扩张。它假定湍流冻结只在形状小幅调整时有效大变形要配合动网格重新求解正问题这是使用这个求解器前必须明确的边界。5. 自带算例当验收基准checkMesh、库朗数与残差的联动检查5.1 换机器后先用自带算例做 smoke test跨机器或换版本之后不要急着跑业务算例先用自带算例验证环境。我用这一组命令在十分钟内判断 OpenFOAM 安装是否正常cd $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavityHighRe ./Allrun跑完之后检查两个点log 文件里残差是否单调下降到 1e-6 量级最后一个时间步是否完整写出。两者都满足说明 MPI 库、编译链路和基础求解器没问题。再跑一个 chemFoam/gri如果点火延迟曲线和已知结果一致说明化学反应相关的库也正常。业务算例出问题时先用这两个基准排除环境因素能省掉大量排查时间。5.2 改网格后必须检查的三个量修改任何算例的 blockMeshDict 之后都要运行 checkMesh 并看三个关键量checkMesh | grep -E non-orthogonality|skewness|minimum volume第一条看最大非正交角低于 70 度算健康高于 85 度基本直接发散第二条看最大歪斜率经验值应小于 4第三条看最小体积绝对不能为负。按这个顺序排查先看最小体积是否为负再看非正交角最后看歪斜率。三个量都过关才值得启动求解器。对于 motorBike 这类经过 snappyHexMesh 的算例网格统计里还会出现 cells 数量和 face 数量的比例异常也能从 checkMesh 输出里看出来。很多启动即发散的案例根因不是松弛因子没调好而是 snappyHexMesh 阶段生成的网格就没达到求解器的最低要求这个检查顺序在 OpenFOAM 里几乎是诊断的第一原则。本文还有配套的精品资源点击获取