用Lattice Boltzmann Method实现卡门涡街

不用去死磕复杂的纳维-斯托克斯方程（Navier-Stokes），竟然能模拟出完美的卡门涡街，这在接触LBM（格子玻尔兹曼法）之前我是不敢想的。这种从微观粒子碰撞分布切入，而非直接解宏观偏微分方程的思路，对开发者来说极其友好。

我在C++实现过程中最核心的逻辑就是处理分布函数的碰撞和迁移。简单来说，就是定义一个离散的速度集，让粒子在格点之间跳跃。

实操中的核心伪代码逻辑大概是这样：

// 简化版的LBM碰撞步骤
for (int i = 0; i < Q; ++i) {
    double feq = calculate_equilibrium(rho, ux, uy, i);
    f[x][y][i] = f[x][y][i] - (f[x][y][i] - feq) / tau;
}
// 迁移步骤：将f[x][y][i] 移动到相邻格点 f[x+cx[i]][y+cy[i]][i]

在部署到超算运行的时候，最容易踩坑的地方在于边界条件的处理。如果Bounce-back（反弹）边界没写好，流体在经过障碍物时会出现奇怪的数值震荡，导致涡流形状崩掉。我排查了很久才发现是索引偏移量算错了，导致粒子在碰撞时没有正确地在格点间对称反弹。

这个方案的优势在于并行化极其简单，每个格点的计算只依赖于周围几个邻居，非常适合用OpenMP或者CUDA去加速。对于想尝试流体模拟但不想被数学公式劝退的开发者，这种基于粒子的实战路径比传统CFD要直观得多。