北京
卡门涡街——钝体后方交替脱落的那串漩涡，你在桥墩、机翼、烟囱背后都见过。它是流体力学里辨识度最高的现象之一，也是检验数值方法的经典试金石。
通常你要老老实实求解纳维-斯托克斯（NS）方程，但这就进入了一个工程噩梦：复杂几何需要贴体网格，每次改外形都得重新剖分，除了极少数玩具形状，解析解根本不存在。有限体积法、有限差分法、有限元法虽然管用，但处理弯曲边界时，要么付出高昂的网格生成代价，要么在壁面附近拼凑特殊的模板修正。几何每变一次，网格就得重来——这是物理问题上面叠了一层工程问题。
现在有一种方法几乎绕开了这一切。它不把 NS 方程当作假设，而是从更底层的统计力学出发，把流体看成大量粒子碰撞的统计分布。你看到的涡量、压力梯度、粘性阻尼，全都作为自然结果涌现出来，而非事先写进方程里。
这个方法叫格子玻尔兹曼方法（LBM），它位于分子动力学和 NS 求解器之间，操作在介观尺度。它的核心对象是单粒子分布函数 f(x, v, t)：不直接问“这一点速度场是多少”，而是问“向各个方向运动的粒子各有多少”。流体行为由此在一个离散的速度格子上演化，每一步只是一个简单的更新规则，却能在恰当的极限下精确恢复 NS 方程——不是假设，而是推导出来的结果。
而这带来的实际好处相当直接。第一，告别复杂网格。因为 LBM 天然在规则的格点上运作，复杂的曲面边界可以用简单的反弹规则处理，无需为每一处几何变动重新剖分网格。第二，极为适合并行。每个格点的碰撞和迁移操作只依赖局部邻居，计算几乎可以无脑拆分到成千上万个核心上，这也是作者后来把它搬到 MareNostrum 5 超级计算机上的原因。
第三，代码量少得出奇。一整套模拟流程——包括边界条件、碰撞算子、宏观量提取——仅用大约 200 行 C++ 就能完成，同时还有一份更简化的 Python/NumPy 版本放在 GitHub 上。你不用在求解器框架上堆砖，就能亲眼看到卡门涡街从数值噪声中自发成形。
当然 LBM 也有自己的适用边界：它在低马赫数、亚声速区间表现最好，在极低粘性下会出现数值不稳定。但作为一种理论优雅、工程轻量的方法，它正在让流体模拟从“解方程”转向“让物理自己演化”。
完整的 C++ 源代码和 Python 版本均已开源，你可以直接跑起来，在不碰一下 NS 方程的情况下，生成属于自己的涡街。
特别声明：以上内容(如有图片或视频亦包括在内)为自媒体平台“网易号”用户上传并发布，本平台仅提供信息存储服务。
Notice: The content above (including the pictures and videos if any) is uploaded and posted by a user of NetEase Hao, which is a social media platform and only provides information storage services.
