Fluid Simulation for Dummies：面向初学者的流体模拟深度解读

背景

这篇文章源自 Hacker News 社区，作者回顾了自己于 2005 年春夏季撰写的硕士论文，主题涉及高性能实时 3D 流体模拟与体渲染。作者在写作过程中发现，尽管现有的参考材料质量很高，但它们过于侧重物理学和数学推导，对于缺乏深厚物理和微积分背景的学习者来说，理解门槛极高。

为了弥补这一空白，作者希望撰写一篇“如果一年前有人能给我看就好了”的指南。本文旨在以编程实现为核心，通过 Jos Stam 的经典论文《Real-Time Fluid Dynamics for Games》（面向游戏的实时流体动力学）中的基础思想，分步骤展示如何从零开始构建一个简单的 3D 流体模拟系统。

核心内容

1. 基础概念：从纳维-斯托克斯方程到离散网格

流体模拟的数学基础是纳维-斯托克斯方程（Navier-Stokes equations）。这是一组描述流体运动的基本方程，虽然有趣且基础，但如果没有良好的物理和微分方程背景，极难理解。

作者建议初学者暂时忽略复杂的方程推导，转而采用**离散网格（Grid-based）**的视角：

• 网格单元：将流体想象为由无数个微小的盒子（或称单元格）组成。
• 属性：每个盒子拥有速度（velocity）和密度（density）等属性。
• 交互：盒子与其邻居相互作用，影响彼此的速度和密度。

在现实世界中，流体包含极多微观粒子，每秒进行无数次交互。但在计算机模拟中，我们必须妥协：将流体划分为合理数量的盒子，并每秒执行有限次数的交互计算。

2. 简化假设：不可压缩流体

为了降低难度，本文仅关注不可压缩流体（Incompressible Fluid）。

• 可压缩流体：如空气，受压后体积缩小。
• 不可压缩流体：如水，受压后体积不变，而是产生反作用力。

不可压缩流体的密度和压力始终保持恒定，这使得模拟变得简单。

3. 可视化技巧：染料密度 vs. 流体密度

这是作者强调的一个关键认知误区。纯水是均匀且透明的，移动时难以肉眼观察。因此，模拟中引入“染料”（dye）来可视化流动。

• 密度（Density）：在代码和模拟中，所谓的“密度”实际上指的是染料的浓度，而非流体本身的物质密度。
• 均匀性：流体本身密度处处相等，但染料分布不均，这种差异让我们看到了流体的运动。

4. 数据结构实现

由于 C 语言对多维数组支持不佳，作者使用一维数组来模拟三维网格。

• 索引宏定义：

  #define IX(x, y, z) ((x) + (y) * N + (z) * N * N)
  

  该宏将三维坐标 $(x, y, z)$ 映射到一维数组的索引。

• FluidCube 结构体： 包含以下核心成员：

  • size: 网格尺寸 $N$。
  • dt: 时间步长。
  • diff: 扩散系数。
  • visc: 粘度系数。
  • s, density: 用于存储染料密度的临时数组和当前密度数组。
  • Vx, Vy, Vz: 当前速度场的三个分量。
  • Vx0, Vy0, Vz0: 上一帧或中间状态的速度分量（用于半拉格朗日法中的回溯）。

  提供了创建 (FluidCubeCreate)、释放 (FluidCubeFree)、添加密度 (FluidCubeAddDensity) 和添加速度 (FluidCubeAddVelocity) 的基础函数。

5. 模拟核心步骤

一个完整的模拟步长（Simulation Step）主要包含三个核心操作和两个辅助子程序。

核心操作

1. 扩散 (Diffuse)：
   • 现象：如同酱油滴入水中，即使水静止，酱油也会逐渐散开。
   • 作用：不仅用于扩散染料，也用于扩散流体的速度场（模拟粘度）。
2. 投影 (Project)：
   • 目的：确保流体保持不可压缩性。
   • 原理：每个网格单元内的流体总量必须守恒（流入量等于流出量）。之前的操作可能导致某些单元格出现净流入或流出，project 操作通过求解压力场，修正速度场，使其满足散度为零的条件，达到平衡状态。
3. 平流/对流 (Advect)：
   • 现象：流体中的粒子随速度场移动。
   • 作用：将染料和速度场本身沿着速度方向进行迁移。

辅助子程序

1. 设置边界 (set_bnd)：
   • 功能：防止流体从模拟容器边缘泄漏。
   • 实现：将最外层单元格视为墙壁。通过镜像反射紧邻墙壁内侧的速度值，使得墙壁处的速度恰好抵消内侧速度，从而形成刚性边界。
2. 线性求解器 (lin_solve)：
   • 功能：用于求解扩散和投影步骤中产生的线性微分方程组。作者坦言其具体数学细节复杂，但它是实现物理真实感的关键数值方法。

模拟主循环 (FluidCubeStep)

完整的模拟步骤按以下顺序执行：

1. 扩散速度：对 $Vx, Vy, Vz$ 三个分量分别进行扩散。
2. 投影速度：修正速度场，确保不可压缩性。
3. 平流速度：根据当前速度移动速度场本身（注意：这里使用的是半拉格朗日法，需要回溯位置）。
4. 再次投影：由于平流操作可能破坏不可压缩性，需再次修正速度场。
5. 扩散染料：对染料密度进行扩散。
6. 平流染料：根据速度场移动染料密度。

关键要点

• 物理简化：通过忽略纳维-斯托克斯方程的复杂推导，采用离散网格和不可压缩假设，将复杂的流体力学转化为可编程的数值计算问题。
• 密度即染料：明确区分流体物质密度（恒定）和染料浓度（变化），后者是可视化的关键。
• 三大算子：流体模拟的核心逻辑由扩散（解决粘性和混合）、投影（解决不可压缩性和质量守恒）、平流（解决物质随流运动）构成。
• 边界处理：通过镜像反射速度值来模拟刚性墙壁，防止数值泄漏。
• 数值方法：使用一维数组模拟三维网格以优化内存访问和兼容性；利用线性求解器处理隐式扩散和压力投影步骤。
• 半拉格朗日法：在平流步骤中，通过回溯粒子位置来确定其携带的属性，这种方法通常比欧拉法更稳定且允许较大的时间步长。

意义与影响

这篇文章及其背后的技术思想对计算机图形学、游戏开发以及科学可视化领域产生了深远影响。

1. 降低入门门槛：Jos Stam 的论文及其衍生教程（如本文） democratized（民主化）了流体模拟技术。它证明了无需深厚的数学背景，程序员即可通过数值方法实现逼真的流体效果，极大地推动了实时流体模拟在游戏中的普及。
2. 实时模拟的基石：通过简化物理模型并采用高效的数值算法（如线性求解器和半拉格朗日平流），使得在普通硬件上进行实时 3D 流体模拟成为可能。这为后来的许多游戏引擎（如 Unity 和 Unreal Engine 中的流体插件）和视觉效果软件奠定了理论基础。
3. 编程与理论的桥梁：文章强调了代码实现的重要性，展示了如何将抽象的物理方程转化为具体的数组操作和循环逻辑，为后续研究者提供了从理论到实践的清晰路径。
4. 可扩展性：虽然本文仅展示了基础的 3D 不可压缩流体，但其