在计算流体动力学（CFD）领域，描述流体流动、混合及其与环境相互作用的方式不止一种。数十年来，**有限体积法（FVM）**一直是行业标准，构成了大多数商业和学术CFD求解器的核心基础。然而，随着应用场景日趋复杂——涉及剧烈晃动、自由液面或运动边界等问题——基于粒子的CFD方法已作为一种极具竞争力的替代方案崭露头角。
两种方法求解的物理方程相同，即纳维–斯托克斯方程。然而，它们表征流体、计算相互作用以及处理几何边界的方式存在本质差异。理解这些差异，有助于工程师针对所需捕捉的物理现象选择合适的数值工具。
从网格到粒子——视角的转变
FVM从欧拉视角描述流体运动。它将仿真域划分为控制体网格，并计算质量、动量和能量在各控制体之间的通量。与之相对，基于粒子的CFD采用拉格朗日方法，直接追踪离散流体微团的运动及其相互作用。
在FVM中，流体在固定或运动网格中流动；而在粒子CFD中，网格不复存在，由动态交互的粒子云来表征流体。这一差异看似细微，却对数值稳定性、灵活性和可扩展性产生了深远影响。
有限体积法
FVM是CFD中的基础性技术，因为它严格遵循质量、动量和能量守恒原理，确保进入某一控制体的物质量等于离开该控制体的物质量。通过对每个网格单元积分守恒方程并利用散度定理，FVM将微分方程转化为求解器可高效处理的代数形式。即便面对复杂的非结构化网格，该框架也能保证全局守恒性。
FVM的优势已得到广泛认可：
- 数值精度高，通量守恒严格。
- 拥有经过数十年验证的丰富湍流与燃烧模型。
- 对壁面、入口和出口等边界条件的处理稳健可靠。
- 在并行计算架构上具有良好的可扩展性。
这些特点使FVM尤其适用于计算域几何形状固定、流动保持连续的场景，例如空气动力学、透平机械或电子散热等领域。然而，FVM对网格质量的依赖也带来了挑战。当几何形状发生变形或界面破碎——如飞溅、破碎或自由液面流动的情况——维持网格质量的难度将大幅增加。
基于粒子的方法
粒子CFD则完全摒弃了网格。每个粒子代表一小团流体，携带各自的速度、压力和热力学属性，粒子之间通过核函数定义局部邻域并进行相互作用。
相邻粒子之间的影响程度可通过如下简单权重函数计算：
拉格朗日粒子CFD的常见变体包括光滑粒子流体动力学法（SPH）和移动粒子半隐式法（MPS）。 由于该方法完全基于拉格朗日框架，粒子的运动自然地随流体流动。无需对网格进行变形，无需重新划分网格，也不存在拓扑约束。
这使得粒子CFD在以下场景中尤为强大：
- 自由液面与晃动流动（如液态燃料箱、波浪–结构物相互作用）。
- 固体边界运动或变形的流固耦合（FSI）问题。
- 多相流与颗粒系统，如悬浮液或粉体流动。
- 相变过程，如熔化、凝固或喷雾动力学。
无网格的特性使前处理大为简化，几何建模更加快捷，尤其适用于复杂或动态变化的计算域。
选择合适的方法
粒子CFD与有限体积CFD之间的选择，更多取决于问题本身的物理特性，而非个人偏好：
- 对于稳态、高雷诺数或热流问题，FVM具有无可比拟的优势。
- 对于高度瞬态、多相流或变形域问题，粒子CFD能有效捕捉传统网格方法难以处理的复杂现象。
经验丰富的仿真工程师越来越倾向于将这两种方法视为互补关系——如同同一数值"乐团"中的不同乐器，各自擅长流体运动的不同方面。
突破CFD的传统边界
目前，粒子CFD和FVM都在积极适配新型计算架构，使两种方法均能实现更大规模、更精细的仿真。FVM领域在运动自适应网格方面的进展，以及拉格朗日CFD中可变粒子尺寸技术的发展，正在进一步拓展仿真的可能性边界，并在一定程度上预示着两种方法在应用领域上的相互融合趋势。
归根结底，有限体积法和粒子法都旨在帮助工程师"看见不可见之物"——准确预测塑造我们所处环境的流体行为。