第 1 章 · CFD 的哲学：为什么用计算机算流体

这一节回答一个最根本的问题：流体方程几百年前就写出来了，为什么还要用计算机去"算"它？答案是：纳维-斯托克斯方程是非线性偏微分方程，除了极少数特例，解析解根本不存在。飞机整机绕流、燃烧室内的湍流燃烧，这些问题的解只能"数值逼近"。CFD 的哲学就是：把连续的流动空间切成有限个点，把连续的导数变成离散的加减乘除，让计算机一步步把解"走"出来。

计算机算流体，到底在算什么

流体力学的三条守恒律——质量守恒、动量守恒、能量守恒——写成数学形式就是一组偏微分方程。方程里的速度 $u$、压力 $p$、密度 $\rho$ 都是空间和时间的连续函数。计算机存不下"连续函数"，只能存有限个网格点上的数值。所以 CFD 的第一步永远是离散化：把"在每一点都成立"改成"在网格点上近似成立"。理解了这一点，全书的所有技术都是在回答两个问题：导数怎么近似？近似会带来什么误差？

解析解与数值解的分界线

• 能求解析解的：驻点流、库埃特流、泊肃叶流等极少数简单外形——它们靠线性化和特殊几何才有解。
• 只能数值解的：任意外形的绕流、有激波的跨声速流、湍流。工程上 99% 的问题在这里。
• 数值解不是"近似公式的解"，而是一个可以控制精度的逼近过程：网格加密一倍，误差应该按理论阶数下降。这个"收敛性"的观念比任何一个具体算法都重要。

CFD 的三步工作流

1. 前处理：画几何、切网格。网格质量直接决定结果可信度。 2. 求解：选方程模型（欧拉 / N-S / 雷诺平均）、选离散格式、迭代到收敛。 3. 后处理：看云图、流线、力的积分，判断物理上是否合理。

CFD 结果永远要经过验证与确认：验证（verification）检查代码是否正确解了方程，确认（validation）检查方程是否正确描述了物理。两者都不能省。

老师的提醒

• 不要一上来就调软件画云图。先手推一遍最简单的问题（一维巴格列夫斯基问题），才知道软件在替你做什么。
• "收敛"有两个含义：迭代残差下降（数值收敛）和结果不再随网格变化（网格收敛）。只看前者会自欺欺人。

自检：能答上来才算过关

• 为什么非线性是解析解的拦路虎？（叠加原理失效，无法用基本解构造复杂解）
• 离散化的本质代价是什么？（把无限维问题变成有限维，引入截断误差）
• 验证与确认的区别是什么？

更深一层：连续介质假设是整座大厦的第一块砖

CFD 的一切合法性都始于一个看似理所当然的断言：流体可以被当作连续介质描述。空气由分子组成，标准状态下每立方厘米约有 $2.7\times10^{19}$ 个分子，分子间平均自由程约 $68\ \mathrm{nm}$。只要我们关心的流动尺度远大于自由程（比如飞机机翼的米级、管道的厘米级），微观离散性就完全被统计平均抹平——密度、速度、温度这些宏观量才能定义为"体积平均后的光滑场"。无量纲判据是克努森数 $Kn = \lambda/L$（自由程与特征尺度之比）：$Kn < 0.01$ 时连续介质假设成立，纳维-斯托克斯方程可靠；$Kn \approx 0.1\sim1$（高空稀薄大气、微机电系统、真空设备）滑移效应显现，需要速度滑移与温度跳跃边界条件修正；$Kn > 10$ 则分子自由流，连续方程彻底失效，必须用玻尔兹曼方程或直接模拟蒙特卡洛（DSMC）。

把这条判据记在手边，你就拥有了判断"CFD 到底能不能描述眼前这个问题"的第一把尺子。工程上最常见的翻车不是方程解错，而是在连续介质假设失效的场合硬用 CFD：高真空的分子泵、纳米间隙的硬盘磁头、100 km 高空的再入初段——这些领域的"仿真结果"若不看 $Kn$，看起来再漂亮也只是数字艺术。

推导细读：从守恒到偏微分方程只要三行

很多人第一次接触控制方程时被积分符号吓退，其实从守恒律到 PDE 的过渡极其朴素。以质量守恒为例：取固定控制体 $\Omega$，"体内质量的增长 = 净流入"，

右侧用高斯定理转为体积分，移项：

由 $\Omega$ 的任意性，被积函数必须逐点为零——偏微分方程就此诞生。"任意性"这一步是整个推导的灵魂：正是它把一个积分陈述升级成了逐点陈述，也正因为如此，弱解（第 7 章激波）只是回到积分形式、放弃逐点要求而已。理解了这条退路，激波捕捉格式的合法性便不再神秘。

在真实世界里：CFD 的成本结构

一个真实项目的算力账单值得初学者早看：波音级整机外流场 2 亿网格、雷诺平均方程、每个物理时间步内迭代 50 次收敛伪时间、模拟 5 秒飞行——单 CPU 核心要算上百年；上千核并行集群也要数天。这笔账决定了工程实践的三条铁律：网格无关性研究只对关键部位局部加密而非全域加密；湍流用 RANS 模型而不是 LES 除非预算充裕；几何能简化就简化（去掉铆钉、小凸台）。CFD 工程师一半的功力在"少算"，另一半才在"算对"。

常见误区与辨析

• "CFD 是精确科学"：它是近似方法的系统化，误差来源至少有四层——模型误差（方程本身，如湍流模型）、离散误差（网格与格式）、迭代误差（收敛不充分）、几何误差（CAD 简化）。报告结果时应能说出每一层的量级。
• "网格越密结果越好"：网格密度与边界条件质量、湍流模型适配性共同决定精度；在错误的边界条件上加密网格只是更快地得到错误答案。
• "验证过商业软件就不用验证了"：软件验证 ≠ 你的算例验证。求解器对，你的网格、物理模型、边界条件仍可能全错。

与其他节的联系

本节的"三步工作流"（前处理—求解—后处理）与第 5 节网格部分、第 8 节边界条件部分一一呼应；"验证与确认"的思想将在第 7 节（与精确黎曼解对比）和第 9 节末（验证清单）反复回响。学完整本书回头看，第 1 章勾勒的地图上每一条路线都走过了——那时你对"CFD 的哲学"会有第二遍的理解，而第二遍才是真的。

专题深入：CFD 的三代求解器架构史

CFD 求解器的架构演进本身就是一部"硬件约束下的算法史"。第一代（1970s）：结构化网格 + 有限差分 + 显式推进——内存按 KB 计，几何只能是简单外形，代表是早期机翼跨声速程序。第二代（1980s–90s）：非结构网格 + 有限体积 + 隐式/多重网格——几何自由 + 收敛加速，Fluent/CFX 在这一代定型。第三代（2000s–）：任意多面体网格 + 自动网格生成 + 大规模并行 + GPU——自动化与规模是主题。每一代的"杀手级进步"都来自一个具体瓶颈的突破：显式的 CFL 束缚 → 隐式；结构化几何限制 → 非结构；CPU 串行极限 → 并行/GPU。读懂求解器架构史，就读懂了"为什么今天的 CFD 长这样"——每个现代特性都对应一个被解决的历史痛点。

数学细节展开：偏微分方程的三种类型与它们的"信息传播方式"

双曲、抛物、椭圆三分法的物理内核是信息传播方式。双曲（波动方程）：信息沿特征线以有限速度传播——下游不能影响上游、间断可以存活、边界条件按特征方向指定。抛物（热传导）：信息瞬时传遍全域（理论上无穷速度）但强度随距离指数衰减——解瞬间光滑化、时间上有"方向性"（正时间演化稳定、反时间病态）。椭圆（拉普拉斯）：无时间之箭、全局耦合——任何局部扰动瞬间影响全域、边界条件必须完整包围求解域。三分法决定了离散方法的性格：双曲用显式推进 + 迎风、抛物用大步长隐式、椭圆用全局迭代（多重网格是VIP待遇）。CFD 里这三种类型混合出现（N-S 方程 = 对流双曲 + 扩散抛物 + 压力椭圆），混合型的数值处理（分步/算子分裂）是高级章节的核心难点。

工程案例细节：一次机翼 CFD 分析的错误排查实录

一次典型的机翼 CFD 排障实录，展示"四本账"的实战用法：现象：升力系数比风洞数据高 8%。排查顺序：先查离散误差（网格加密 → 升力变化 1.2% —— 不是主因）→ 查迭代收敛（残差降了 4 个量级 ✓）→ 查湍流模型（SA → SST 后升力降 3% —— 部分原因）→ 查边界条件（远场距离 20 倍弦长改成 50 倍，升力降 2% —— 剩余原因）→ 查几何（忽略了翼尖小翼的几何细节 —— 补上后基本吻合）。教训：误差从不来自单一账本——这份排障清单值得贴在工位上。

常见误区与辨析（补充）

• "CFD 算的是'真实'流场"：CFD 算的是离散化后的模型方程的解——模型误差（湍流模型、化学模型）永远存在；把 CFD 当"虚拟风洞"是对的，把它当"真相"是危险的。
• "后处理云图越漂亮越可信"：云图的颜色映射可以掩盖一切；关键量（力、力矩、流量）的数值与网格无关性才是可信度的证据。