格子玻尔兹曼方法

格子玻尔兹曼方法（lattice Boltzmann method, LBM）不直接离散 Navier-Stokes 方程，而是离散更底层的 Boltzmann 方程：在每个格点上维护一组离散速度的粒子分布函数，"碰撞"在格点内局部发生，"流送"把分布函数显式搬到邻格。它由格子气自动机的统计平均演化而来（McNamara-Zanetti 1988McNamara & Zanetti 1988 (PRL)，见 lattice-gas-automata），经典综述见 Benzi-Succi-VergassolaBenzi et al. 1992 (Phys. Rep.)。

图:AA-pattern 就地流送格式——同一数据结构在奇偶时间步交替充当"发射"与"接收"数组,免去双缓冲(来源:QC068)

从连续 Boltzmann 到格子

出发点是连续 Boltzmann 方程与 BGK 碰撞项：

$$\frac{\partial f}{\partial t}+\mathbf{\xi }\cdot \nabla f=-\frac{1}{{\tau }_c}\left(f-f^{\mathrm{eq}}\right),$$

其中 $f(\mathbf{x},\mathbf{\xi },t)$ 是单粒子分布函数，$f^{\mathrm{eq}}$ 是 Maxwell-Boltzmann 平衡态，${\tau }_c$ 是碰撞松弛时间。做三件事：（1）把速度空间截断为有限离散速度集 $\{{\mathbf{c}}_i\}$（如 D2Q9，见 d2q9-lattice）；（2）令速度模长恰好使 $\mathbf{x}+{\mathbf{c}}_i\Delta t$ 落在邻格上；（3）时间用显式欧拉。得到格子 Boltzmann 方程（LBE）：

$$f_i(\mathbf{x}+{\mathbf{c}}_i\Delta t,t+\Delta t)=f_i(\mathbf{x},t)-\frac{\Delta t}{{\tau }_c}\left[f_i(\mathbf{x},t)-f_i^{\mathrm{eq}}(\mathbf{x},t)\right].$$

宏观量按矩恢复：$\rho ={\sum }_i{f}_i$、$\rho \mathbf{u}={\sum }_i{\mathbf{c}}_i{f}_i$。通过 Chapman-Enskog 多尺度展开（见 chapman-enskog）可以证明，宏观极限正是 Navier-Stokes 方程，粘度为 $\nu =c_s^2({\tau }_c-\frac{1}{2})\Delta t$。

介观图像：为什么这是”另一个世界”

• 变量是分布函数而非场量：未知量个数是”格点数 × 速度数”，比 NS 多一个因子（D2Q9 是 9 倍），但每个分量都活在规则格子上、无几何复杂度；
• 碰撞完全局部：$f_i^{\mathrm{coll}}=f_i-\frac{1}{\tau }(f_i-f_i^{\mathrm{eq}})$ 只依赖本格点的矩，两格点之间没有任何通信；
• 流送是完全显式的数据搬移：右端就是”把数组沿固定方向平移一格”（上图 AA-pattern 即其就地实现），不含任何浮点通信。

后两条使 LBM 在 GPU 与 Exascale 机器上接近峰值性能，Succi 组的软物质综述系统讨论了混合精度（population shifting）等实现技术Tiribocchi et al. 2025 (Phys. Rep.)。

与有限差分的结构差异

把 LBE 看成”对流项的有限差分”会错过要点。对 NS 格式（见 finite-difference-method），非线性在物理空间的通量 $\mathbf{u}\cdot \nabla \mathbf{u}$ 里，需要邻点通信与稳定化（迎风、限制器）。对 LBM：

• 对流（流送）是线性且精确的——分布函数沿特征线精确移位，无数值扩散；
• 全部非线性集中在局部平衡态 $f_i^{\mathrm{eq}}\propto \rho [1+({\mathbf{c}}_i\cdot \mathbf{u})/c_s^2+({\mathbf{c}}_i\cdot \mathbf{u}{)}^2/2c_s^4-u^2/2c_s^2]$ 里，其强度由马赫数 $Ma=u/c_s$ 控制；
• 因此 LBM 的非线性度由 $Ma$ 而非 $Re$ 决定：固定 $Re$ 降 $Ma$（小速度或重标定），非线性可以任意减弱而物理不变。

这条性质是 Carleman 线性化量子 LBM 的理论根基——量子 Carleman 嵌入的维度代价随非线性强度增长，而 LBM 允许在不动物理的前提下把非线性压到机器可以处理的水平Sanavio & Succi 2024 (IntechOpen)，详见 carleman-qlbm 与 carleman-convergence-analysis。

局限与代价

LBM 的标准格子是均匀笛卡儿的，贴体复杂几何需要 immersed boundary、cut-cell 或格子细化弥补；标准 LBM 只到 $O(M{a}^2)$ 精度，不适合高亚声速以上；低粘度（$\tau \to 1/2^{+}$）时 BGK 碰撞数值失稳，需要 MRT 碰撞（见 bgk-mrt-collision）。量子化后的图景见 quantum-lbm-overview：流送是格点平移的酉算子、碰撞是非酉非线性映射，两者的组合构成结构性矛盾。

延伸阅读

• lattice-gas-automata
• bgk-mrt-collision
• d2q9-lattice
• chapman-enskog
• quantum-lbm-overview
• carleman-qlbm