本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第十二篇（完结篇），覆盖教材第13-14章「湍流模型与现代数值方法」的核心内容。从工程CFD的奠基到高保真仿真前沿，为这个系列画上句号。

第十二篇：湍流模型与现代数值方法——从工程 CFD 到高保真仿真

一、CFD 的百年演变

1.1 概念引入

如果你在 1950 年问一位空气动力学家如何计算一个翼型的升力和阻力，答案会是：去风洞。如果你在 2026 年问同样的问题，答案会是：先用 CFD 算一遍，再用风洞验证几个关键工况。

计算流体力学（Computational Fluid Dynamics, CFD） 已经从学术工具变成了工程标准。

1.2 历史节点

年份 事件 影响
1910 Richardson 提出”数值天气预报”概念 CFD 最早的思想
1930s 手算松弛法处理势流 电脑出现前的数值方法
1946 第一台电子计算机 ENIAC 诞生 开始有了”计算”的硬件
1960s NASA Ames 的团队开发隐式格式 可压缩流 NS 方程首次数值求解
1970s Jameson 的有限体积法 + Runge-Kutta 工业 CFD 的基石
1980s Spalart-Allmaras 湍流模型 工程用单方程模型
1990s-今 LES/DES/DNS 高保真模拟

二、湍流模型——工程的”权宜之计”

2.1 为什么需要湍流模型

直接求解湍流的最小尺度需要网格量 N∝Re9/4N \propto Re^{9/4}。对于 Re=107Re=10^7 的机翼，需要约 101510^{15} 个网格点——这远远超出了目前任何一台计算机的能力。

因此，工程上采用湍流模型（Turbulence Model）来近似湍流效应——用 Reynolds 平均 Navier-Stokes 方程（RANS），把湍流的影响包在”涡黏性”中。

2.2 湍流模型的层次

1
2
3
4
5
6
7
DNS (直接数值模拟)                    精度最高，成本最高
         ↑
LES (大涡模拟)
         ↑
DES (分离涡模拟)
         ↑
RANS + 湍流模型                       精度最低，成本最低

2.3 主要工程湍流模型

模型 方程数 特点 典型应用
混合长度模型 0 简单，但无法处理分离 平板、简单管道
Spalart-Allmaras (SA) 1 适合外气动，鲁棒性好 机翼、飞机外流
k−εk-\varepsilon 2 高雷诺数标准模型 内流、燃烧
k−ωk-\omega SST 2 结合 k−ωk-\omega 近壁 + k−εk-\varepsilon 远场 分离流、机翼（推荐）
SST + 转捩模型 4 可预测层流→湍流转捩 低速翼型、层流翼型
v2−fv^2-f 4 壁面处理精确 复杂几何

2.4 工程中最常用的推荐

对于大多数空气动力学外流问题（机翼、机身、无人机），推荐使用 k−ωk-\omega SST 模型（Menter, 1994）：

  • 近壁区：k−ωk-\omega 行为（对逆压梯度敏感，适合分离预测）
  • 远场：k−εk-\varepsilon 行为（对自由来流参数不敏感）
  • 通过混合函数 F~1\tilde{F}_1 平滑过渡

三、数值方法——从偏微分方程到代数方程

3.1 离散化方法

把偏微分方程转化为计算机能求解代数方程的方法有三种：

方法 原理 优/缺点
有限差分法 (FDM) 用差商代替偏导数 简单，但复杂几何难
有限体积法 (FVM) 对控制体积分 守恒性好，适合复杂几何
有限元法 (FEM) 用基函数逼近 适合结构分析，空气动力学中较少用

现代 CFD 空气动力学求解器几乎都基于有限体积法（FVM）。

3.2 网格类型

网格类型 特征 典型软件
结构化网格 六面体，坐标线对齐 ICEM, Pointwise
非结构化网格 四面体/多面体 ANSYS Meshing, Fluent
笛卡尔网格 各向同性，无需复杂几何适应 OpenFOAM, STAR-CCM+
重叠网格 多块网格叠加 OVERFLOW, FUN3D
自适应网格 按误差指示自动加密 Fluent, OpenFOAM

3.3 时间推进

隐式格式（Implicit）：允许大时间步，稳定性好，适合稳态问题
显式格式（Explicit）：时间步受限（CFL 条件），适合瞬态/高速问题

对于定常流动的空气动力学计算，几乎都使用隐式格式。

四、从 CFD 到真实飞行——验证与确认

4.1 误差源

CFD 结果与实验之间的差距来自：

误差源 性质 如何减小
离散误差 网格不够密 网格收敛性研究
模型误差 湍流模型近似 更精细的模型/LES
数值误差 格式精度 高阶格式（2→4阶）
边界条件 入口/出口不准确 更接近实际的风洞条件

4.2 网格收敛性

Anderson 强调：没有网格收敛性研究的 CFD 结果没有意义。

标准做法：

  1. 粗、中、细三套网格
  2. 目标量（CL,CDC_L, C_D）随网格加密趋于稳定
  3. 用 Richardson 外推估计离散误差

4.3 验证 vs 确认

  • 验证（Verification）：你的代码是否正确解出了方程？
  • 确认（Validation）：你的方程是否正确描述了现实？

两者缺一不可。

五、工业级 CFD 的工作流程

5.1 一个典型的外气动计算流程

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
步骤1：几何准备
    └── CAD 模型 → 清理 → 简化

步骤2：网格生成
    └── 表面网格 → 体网格 → 边界层加密
    └── 检查: y+, 偏度, 正交性

步骤3：求解器设置
    ├── 边界条件: 入口(速度/总压), 出口(静压/远场)
    ├── 湍流模型: k-ω SST (推荐)
    ├── 离散格式: 2阶迎风 (默认)
    ├── 收敛判据: 残差 < 1e-6
    └── 监控: CL, CD 收敛

步骤4：后处理
    ├── 压力分布 Cp
    ├── 摩擦阻力 Cf
    ├── 流线/涡结构
    └── 升力/阻力/力矩

步骤5：验证
    └── 网格收敛性
    └── 与实验/基准解对比

5.2 常用工业软件

软件 类型 特点
ANSYS Fluent 商用 功能最全面
STAR-CCM+ 商用 多物理场耦合
OpenFOAM 开源 高度自定义
SU2 开源 气动优化设计
OVERFLOW NASA 重叠网格强项
FUN3D NASA 非结构网格强项

六、高保真仿真前沿——LES/DNS/DES

6.1 大涡模拟（LES）

对Re=106Re=10^6 的机翼，LES 需要约 10910^9 网格点。这比 RANS 的 106−10710^6-10^7 高了 2-3 个数量级，但比 DNS 的 101210^{12} 现实得多。

LES 的原理：

  • 大尺度涡：直接求解（受几何边界影响大，需要解析）
  • 小尺度涡：次网格模型（”通用”的，可以模化）

6.2 分离涡模拟（DES）

DES 是一种 RANS/LES 混合方法：

  • 近壁面（边界层内）：RANS（k−ωk-\omega SST）
  • 远离壁面（分离区、尾迹）：LES

这种方法对分离流特别有效，计算成本远低于完整 LES。

6.3 直接数值模拟（DNS）

DNS 求解所有尺度的湍流——从最大的流动尺度到最小的 Kolmogorov 尺度。

DNS 的网格量：N∝Re9/4N \propto Re^{9/4}

场景 ReRe 所需网格 目前的可行性
平板管道流 10410^4 10810^8 ✅ 可行
低雷诺数翼型 10510^5 101010^{10} ⚠️ 大型集群
全机（巡航） 10710^7 101510^{15} ❌ 遥不可及

DNS 在可预见的未来仍将局限于低雷诺数的基础研究。

七、完整概念地图

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
空气动力学基础系列完整知识体系

第一篇: 流体力学基础
└── 连续介质、伯努利方程、马赫数
    ↓
第二篇: 势流理论(上)
└── 拉普拉斯方程、四种基本解
    ↓
第三篇: 势流理论(下)
└── 库塔-儒科夫斯基定理、保角变换
    ↓
第四篇: 翼型空气动力学
└── 薄翼理论、NACA 翼型几何
    ↓
第五篇: 翼型特性与数据
└── 真实风洞数据、选型指南
    ↓
第六篇: 有限翼理论
└── 升力线、诱导阻力、展弦比
    ↓
第七篇: 有限翼实验与修正
└── 后掠翼、三角翼、地面效应
    ↓
第八篇: 可压缩流基础
└── 等熵流、正激波、拉瓦尔喷管
    ↓
第九篇: 超声速线化理论
└── 斜激波 θ-β-Ma、膨胀波
    ↓
第十篇: 超声速翼型设计
└── 菱形/超临界翼型、波阻
    ↓
第十一篇: 边界层理论
└── 层流/湍流、转捩、分离
    ↓
第十二篇: CFD 与湍流模型
└── RANS/LES/DNS、有限体积法

八、核心公式速查卡

公式/方法 含义 层级
RANS: ui′uj′‾=νt(∂ui‾∂xj+∂uj‾∂xi)−23kδij\overline{u_i'u_j'} = \nu_t \left(\frac{\partial \overline{u_i}}{\partial x_j} + \frac{\partial \overline{u_j}}{\partial x_i}\right) - \frac{2}{3}k\delta_{ij} 雷诺应力的Boussinesq假设 工程 CFD
k−ωk-\omega SST: DkDt=Pk−β∗kω+∇⋅[(ν+σkνt)∇k]\frac{Dk}{Dt} = P_k - \beta^* k\omega + \nabla\cdot[(\nu+\sigma_k\nu_t)\nabla k] 湍动能方程 工程 CFD
LES: τijsgs=−2νsgsS‾ij\tau_{ij}^{\text{sgs}} = -2\nu_{\text{sgs}}\overline{S}_{ij} 次网格应力模型 高保真
DNS: Ngrid∝Re9/4N_{\text{grid}} \propto Re^{9/4} DNS 网格需求估计 前沿
CFL: Δt≤Δx∣V∣+a\Delta t \le \frac{\Delta x}{|\mathbf{V}| + a} 显式时间步限制 数值方法

九、给自学者的最终建议

12 篇笔记到此结束。如果要从头到尾自学 Anderson《Fundamentals of Aerodynamics》，建议：

  1. 第一篇到第六篇：这是核心。势流 → 薄翼 → 有限翼。理解和掌握 80% 的公式推导。
  2. 第八篇到第十篇：可压缩流。至少理解等熵关系和激波/膨胀波的基本概念。
  3. 第十一篇：边界层。理解转捩和分离的物理机制。
  4. 第七与第十二篇：有兴趣再深入。前面不掌握，直接读涡面法和湍流模型意义不大。

手算练习是必须的——用纸笔推一次 Cl=2παC_l = 2\pi\alpha，用计算器算一次激波关系式，比读十遍理论更有意义。

参考文献

  1. Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapters 13-14). McGraw-Hill.
  2. Menter, F. R. (1994). “Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications”. AIAA Journal, 32(8), 1598-1605.
  3. Spalart, P. R., & Allmaras, S. R. (1992). “A One-Equation Turbulence Model for Aerodynamic Flows”. AIAA Paper 92-0439.
  4. Jameson, A., Schmidt, W., & Turkel, E. (1981). “Numerical Solutions of the Euler Equations by Finite Volume Methods Using Runge-Kutta Time-Stepping Schemes”. AIAA Paper 81-1259.
  5. Pope, S. B. (2000). Turbulent Flows. Cambridge University Press.
  6. Ferziger, J. H., & Perić, M. (2002). Computational Methods for Fluid Dynamics (3rd ed.). Springer.
  7. Wilcox, D. C. (2006). Turbulence Modeling for CFD (3rd ed.). DCW Industries.

📚 Anderson《空气动力学基础》读书笔记系列完

12 篇笔记覆盖了从低速到超声速、从势流到黏流、从翼型到有限翼、从理论到 CFD 的完整知识体系。感谢阅读。

下一篇系列预告：可能的方向包括——

  • 《飞行器稳定性与控制》
  • 《螺旋桨与旋翼空气动力学》
  • 《推进系统基础》

感谢你的阅读，如果哪一篇有疑问，随时留言讨论。