计算流体力学(computational fluid dynamics, CFD)把连续的流体力学方程变成计算机可解的大规模代数方程组,是航空航天、气象、能源与生物医学工程的核心设计工具。理解它的流水线与成本结构,是判断"量子计算机该从哪一步切入"的前提——量子 CFD 的几乎全部提案,都是瞄准这条流水线中的某一站。

流水线:从物理到代数
经典 CFD 的一次求解可以拆成五站:
- 物理建模:从 Navier-Stokes 方程(见 navier-stokes)出发,决定是否忽略压缩性、热效应、湍流脉动(DNS/LES/RANS 三级取舍见 turbulence-simulation),或改用介观模型(见 lattice-boltzmann)。
- 几何与网格:把连续的计算域剖分成有限个节点或单元。网格数 直接决定未知量个数,也决定量子算法需要编码的态空间维度 。
- 离散化:把偏微分方程变成代数方程组 。四条主流路线——finite-difference-method、finite-volume-method、finite-element-method、spectral-methods,加上格子玻尔兹曼这条”不走宏观方程”的介观路线,各有自己的代数结构。
- 代数求解:稳态问题得到一个大稀疏线性系统;非稳态问题每个时间步都要解一次(压力泊松方程是其中的大头,见 projection-method 与 multigrid)。线性求解通常占不可压求解器的主要计算时间。
- 后处理、验证与确认(verification & validation, V&V):验证问”方程解对了吗”(网格收敛性、与解析基准对比,见 classical-benchmarks),确认问”方程本身对吗”(与实验对比)。
离散化是第一个分岔点
对量子 CFD 而言,离散化选择比经典计算更重要,因为它同时决定了矩阵 的结构与编码成本:
- 有限差分给出规则的稀疏模式(如三维泊松的 7 点模板),天然适配稀疏哈密顿量模拟与块编码;
- 有限体积是工业 CFD 的主力,本源量子的混合框架即用 FVM 离散加变分量子线性求解器求解Ye et al. 2024 (Phys. Fluids);
- 有限元带来灵活几何但矩阵结构不规则,量子有限元需专门处理刚度矩阵装配Montanaro & Pallister 2016 (PRA);
- 谱方法的基函数与量子傅里叶变换天然同构(见 quantum-spectral-method);
- LBM 则把非线性从对流项搬到局部碰撞项,成为量子线性化路线的最爱(见 carleman-qlbm)。
成本结构与 Exascale 现状
CFD 的计算成本由湍流主导。直接数值模拟(DNS)的总代价按 标度:三维网格数约 (Kolmogorov 耗散尺度 ,见 kolmogorov-cascade),时间步数约 ,二者相乘即 Sanavio & Succi 2024 (IntechOpen)。据此估计,全尺寸飞机设计()所需的自由度约为 ——存储这些自由度只需约 80 个逻辑量子比特(),而经典机器即使进入每秒 次浮点运算的 Exascale 时代(2022 年前后 Frontier 等系统达成),全尺寸湍流模拟仍在可望而不可即之列;工业界的量子路线图已经以 Frontier 为对标基准评估量子加速空间Syamlal et al. 2024 (AIAA)。
需要立刻泼冷水的是:80 量子比特只是存储账。读取、态制备、时间推进的重复开销会把实际资源推高若干数量级(见 io-problem 与 quantum-advantage-debate)。
量子切入点的谱系
按切入站点不同,量子 CFD 提案自然分族:切进”代数求解”站的是 HHL 系线性系统算法与变分求解器(见 hhl、vqls);把整条演化看成动力学去模拟的是 Carleman 线性化与 Schrödinger 化(见 carleman-linearization、schrodingerization);绕开宏观方程直接量化介观分布函数的是量子 LBM 家族(见 quantum-lbm-overview)。这条谱系的完整地图见 qcfd-overview 与 qcfd-taxonomy。