简介：这份PDF面向刚接触开源CFD工具的研究生、工程师与科研人员，系统梳理OpenFOAM的来龙去脉与核心能力，帮助读者在动手建模前建立整体认知。内容从帝国理工学院Gosman小组与Charlie Hill博士论文的源头讲起，延伸到Henry Weller、Hrvoje Jasak的数值工作及2004年公共域版本的演进，并逐项展开有限体积积分离散、非结构化多面体网格、以类封装离散与数值算法、大规模并行计算等特性，还涉及自动网格运动、拓扑变化、拉格朗日粒子追踪与实时选择表定制等高级功能，配合分层实现与求解器代码片段说明设计理念。资源包仅含1个PDF文件，大小约162KB，篇幅精炼、便于快速通读与随查随用。已有557人学习，适合作为风工程OpenFOAM基础培训的入门讲义，也可为后续求解器选型与二次开发提供背景参照。
1. 风工程入门第一课：OpenFOAM 到底是个什么东西
如果你刚听说“计算流体力学”这个词，脑子里冒出来的画面大概是满屏公式和看不懂的曲线，然后本能地想关掉页面。别急，我第一次接触 OpenFOAM 的时候也一样，看到“面向对象的计算力学数值模拟库”这串字，第一反应是——这到底是软件还是编程语言？其实你可以把它理解成一个“流体问题的乐高工具箱”：它不直接给你一辆拼好的车，而是给你轮子、底盘、马达，让你按自己的需求搭出想要的东西。OpenFOAM 就是这样一个开源的计算流体力学（CFD）工具，用 C++ 写成，专门用来模拟空气怎么流过建筑、水怎么绕过桥墩这类问题。风工程入门之所以绕不开它，是因为它免费、可改、能并行算大网格，而且全球做风工程和湍流研究的课题组几乎都在用。这份培训讲稿就是带你从“它是什么”开始，一步步走到“我能跑一个算例”的入门材料，零基础也能跟。
2. 为什么风工程新手值得学 OpenFOAM：先搞懂它能干什么
2.1 风工程里那些“看不见的麻烦”，它正好能接住
风工程要解决的问题，说白了就是“风碰到建筑、桥梁、地形之后会怎样”。现实里你不可能每次都去建个风洞吹一遍，成本高、周期长，而且有些场景根本没法做缩尺实验。OpenFOAM 的价值就在于，它把流体运动的基本方程——连续性方程和动量方程——用有限体积法离散到网格上，然后一步步迭代求解，让你在电脑里“吹一遍风”。这份讲稿里反复强调的一个能力是：它用积分形式对偏微分方程做有限体积离散，所以最多处理二阶方程，而传递方程刚好落在这个范围里。翻译成人话就是：风工程里最常见的流动和传热问题，它都能算，而且算得有理有据。你不需要先成为数学系毕业生，只要知道“它能把风场算出来”这件事就够了。
更实际的一点是，风工程里经常遇到复杂几何——不是简单的圆管或平板，而是带挑檐的建筑、带附属构件的桥梁断面。OpenFOAM 支持非结构化网格，多面体单元也能用，这意味着你可以把网格贴着复杂表面铺上去，而不是硬把几何简化成方块。讲稿里也提醒了一句：精度和网格相关性较大。这句话是血泪经验，后面避坑章节会展开说。对新手来说，先记住“它能处理复杂形状，但网格质量你得盯着”就够了。
2.2 它和商业软件最大的差别：你能改，而且改得动
很多人问，既然有成熟的商业 CFD 软件，为什么还要学 OpenFOAM？讲稿里给了一个很实在的答案：通用 CFD 软件不可能预先知道并满足所有求解器需求，所以写不出完全通用的求解器。OpenFOAM 换了个思路——它不给你一个万能黑匣子，而是给你一堆数值库和工具，让你通过代码重用来搭自己的求解器。讲稿里展示了一段方程表达，长这样：
solve ( fvm::ddt(k) + fvm::div(phi, k) - fvm::laplacian(nu() + nut, k) == nut*magSqr(symm(fvc::grad(U))) - fvm::Sp(epsilon/k, k) );这段代码是湍动能 k 的输运方程。你不用现在就看懂每一行，只需要感受一件事：方程在代码里几乎是一行一行写出来的，ddt 是时间项，div 是对流项，laplacian 是扩散项。这意味着当你想改一个湍流模型、换一种边界条件、加一个源项，你改的是方程本身，而不是在图形界面里翻半天找不到入口。讲稿里把这叫做“分层实现”：张量代数、场代数、离散方法、网格处理、线性方程组求解、物理模型，各层可以独立开发又相互依赖。对风工程新手来说，这个设计最大的好处是——你不需要一次搞懂全部，可以先用现成的求解器跑起来，再一层一层往下挖。
2.3 第一个例子怎么迈步：从“看得懂目录”开始
零基础最容易卡住的地方不是方程，而是“我装完了，然后呢”。讲稿里没有给安装步骤，但按新手最容易上手的路径，我建议你先做一件事：把 OpenFOAM 的案例目录翻一遍。安装完成后，在终端里输入：
# 进入 OpenFOAM 的用户案例目录（具体路径随安装方式不同略有差异） cd $FOAM_RUN # 把系统自带的经典算例复制到自己的运行目录 cp -r $FOAM_TUTORIALS/incompressible/simpleFoam/pitzDaily . cd pitzDaily这三行命令做了一件事：把一个叫 pitzDaily 的二维管道后台阶算例复制到你自己的目录里。为什么选它？因为它网格小、物理简单、跑得快，是很多人跑通的第一个算例。复制完之后，你会在目录里看到三个关键文件夹：0放初始条件和边界条件，constant放网格和物性参数，system放求解控制。先别急着运行，打开system/controlDict看一眼，里面写着时间步长、总时长、输出间隔。再打开0/U看一眼，里面写着入口速度、出口压力、壁面条件。这一步不需要你改任何东西，目的是建立“文件对应物理量”的直觉。讲稿里提到 OpenFOAM 用字典文件处理输入输出，关键字和值成对出现，用分号隔开，次序无关。你看到的每一个条目，都是在告诉求解器“这个边界上风怎么吹”。
跑起来只需要一行：
# 用 simpleFoam 求解器跑当前目录下的算例 simpleFoam如果终端开始刷迭代信息，说明你迈出了第一步。跑完之后用paraFoam打开看速度云图，你会看到台阶后面有一个回流区。这个回流区就是风工程里最常打交道的现象之一——钝体尾流。第一次跑通的意义不在于结果多精确，而在于你亲手走完了“网格→边界条件→求解→后处理”这条链路。
3. 零基础怎么读懂 OpenFOAM 的结构：从字典文件到求解器
3.1 先认识三个文件夹，别一上来就啃 C++
讲稿里花了不少篇幅讲元类、容器类、字典，这些是给要改代码的人准备的。对刚入门的风工程新手，我建议先把注意力放在三个文件夹上，而不是一头扎进 C++ 源码。0文件夹里放的是场文件，比如速度 U、压力 p、湍动能 k，每个文件开头写着边界类型：入口用 fixedValue 给定速度，出口用 zeroGradient 让流体自由流出，壁面用 noSlip 让速度归零。constant文件夹里放polyMesh（网格）和transportProperties（物性），风工程里空气的黏度和密度就在这里改。system文件夹里放controlDict（时间控制）、fvSchemes（离散格式）、fvSolution（线性求解器设置）。讲稿里说字典允许嵌入表，表的表也能嵌，所以你会在fvSolution里看到多层嵌套，别慌，一层一层展开看，每一层都在回答一个具体问题：用哪个求解器、容差多少、最多迭代几次。
3.2 求解器不是越高级越好，先选对“档位”
OpenFOAM 自带一堆求解器，名字看起来像乱码：simpleFoam、pisoFoam、pimpleFoam、buoyantSimpleFoam。讲稿里提到它提供面向问题的求解器框架，而不是一个万能求解器。对风工程新手，你只需要先记住两个：定常不可压用 simpleFoam，非定常不可压用 pisoFoam 或 pimpleFoam。建筑外流场如果只关心平均风压，simpleFoam 够用；如果要看涡脱落，就得用非定常求解器。选错求解器的典型症状是：残差一直不降，或者跑着跑着发散了。这时候先别怀疑网格，回头看一眼system/controlDict里的application字段写的是不是你以为的那个求解器。讲稿里强调“用户修改要预先定义”，意思是很多行为在字典里已经定好了，你改对地方才有效。
3.3 并行计算不是玄学，先在小算例上试一次
讲稿里提到 OpenFOAM 支持大规模并行化计算，底层做了并行化实施。对新手来说，并行听起来很高级，其实操作门槛不高。以 pitzDaily 为例，你先用decomposePar把网格切成几块：
# 在 system/decomposeParDict 里设置切分块数，比如 4 块 # 然后执行网格分解 decomposePar # 用 4 个进程并行运行求解器 mpirun -np 4 simpleFoam -parallel # 跑完后把结果合并回一个网格 reconstructPar这三步对应“切分→并行算→合并”。讲稿里提到 Halo-free 并行化，你不用深究它怎么实现，只需要知道并行之后每个进程算一部分网格，进程之间交换边界信息。新手最容易翻车的地方是decomposeParDict里的块数和mpirun -np的进程数对不上，或者切分方向设得不合理导致某个进程负载过重。第一次试并行，建议就用 4 块、4 进程，跑通再调。
4. 新手最容易踩的五个坑：现象、原因和解决
4.1 残差曲线一路往上飙，最后发散
现象：终端里初始残差越跑越大，几十步之后报浮点异常退出。原因通常不是软件坏了，而是边界条件或初始条件给得不合理。比如入口给了很大的速度，出口却设成了固定压力且方向反了，或者初始场里某个场是零导致除零。解决：先把时间步长调小，把fvSolution里的松弛因子调低（比如 U 从 0.7 降到 0.3），再看0文件夹里每个场的边界类型是否和物理场景匹配。讲稿里说“用户修改要预先定义”，边界条件就是最需要预先想清楚的部分。
4.2 网格检查报负体积，求解器直接拒绝启动
现象：运行checkMesh时提示负体积或高度扭曲单元。原因是非结构化网格在复杂几何附近容易生成质量差的单元，讲稿里也提醒精度与网格相关性较大。解决：回到网格生成环节，检查是否有单元翻转，适当加密或调整网格尺寸过渡。如果用的是别人给的网格，先用checkMesh过一遍再跑，别跳过这一步。
4.3 跑完了但云图一片均匀，看不出回流或尾流
现象：后处理打开发现速度场几乎没变化。原因多半是迭代步数不够，或者收敛判据设得太松，求解器还没算到有结构的流场就停了。解决：把controlDict里的endTime加大，把fvSolution里的残差目标调严一点，再跑一遍。另一种可能是后处理时选错了时间目录，看的是第 0 步的初始场。
4.4 并行跑完结果和单核不一样
现象：并行算出来的场和单核结果对不上。原因通常是并行分解后某些边界条件没有正确传递，或者reconstructPar合并时漏了时间步。解决：先用小算例对比单核和并行的结果，确认差异在可接受范围内再上大算例。讲稿里提到并行化是底层实施的，但用户侧的字典设置仍然要自己负责。
4.5 改了字典文件但求解器好像没反应
现象：改了fvSchemes里的离散格式，结果完全没变。原因可能是改错了文件路径，或者字典里有关键字拼写错误导致 OpenFOAM 静默忽略了那一项。解决：养成改完字典先跑一遍、看终端有没有警告的习惯。OpenFOAM 对未知关键字有时不报错，只是不生效，这是新手最容易郁闷的地方。
5. 下一步学什么：从跑通一个算例到能改一个模型
跑通 pitzDaily 只是起点。接下来我建议你按这个顺序往下走，每一步都别贪多。第一步，把 pitzDaily 的网格换成你自己关心的一个简单几何，比如一个二维方柱，用blockMesh画结构化网格，体会“几何→网格→边界条件”的对应关系。第二步，把 simpleFoam 换成 pisoFoam，把定常改成非定常，观察方柱后面的涡脱落，这一步会让你真正理解时间步长和库朗数的关系。第三步，打开fvSchemes把对流项从一阶迎风换成二阶迎风，对比结果差异，你会直观感受到离散格式对数值耗散的影响。第四步，如果你对湍流模型感兴趣，去翻讲稿里提到的 RANS 和 LES 部分，先改一个模型系数跑一遍，再试着加一个源项。
下面这张表可以作为你前两周的练习路线参考：
| 阶段 | 目标 | 关键动作 | 判断成功的标准 |
|---|---|---|---|
| 第 1 周 | 跑通自带算例 | 复制 pitzDaily，运行 simpleFoam，paraFoam 看云图 | 能看到台阶后回流区 |
| 第 2 周 | 换几何 | 用 blockMesh 画方柱网格，设边界条件 | checkMesh 通过，残差下降 |
| 第 3 周 | 改格式 | 对比一阶和二阶迎风结果 | 能说出两者尾流长度差异 |
| 第 4 周 | 试并行 | decomposePar 切 4 块，mpirun 跑通 | 并行结果与单核趋势一致 |
讲稿里有一句话我印象很深：如果重用很难，继续分解代码，并重写代码。这句话放在学习路径上也成立——如果某个算例你跑不通，就把它拆成更小的部分，先跑一个更简单的。我第一次学的时候，在一个复杂建筑外流场上卡了两周，后来退回二维方柱重新跑了一遍，才发现问题出在出口边界条件设反了。从那以后我每次新建算例，都强制先跑一个最简版本，确认边界条件和网格没问题，再往上加复杂度。希望这份入门路径能帮你少走一点弯路，先把第一个算例跑起来，剩下的边跑边学。
本文还有配套的精品资源，点击获取
