不久前，我在 Rust 中使用 BGK 近似和 D2Q9 晶格上的 TRT（双弛豫时间）方法编写了一个 Lattice Boltzmann Solver。在这两种情况下，即使在简单的示例中，我也会遇到主要的稳定性问题。这是一个模拟带有圆柱形障碍物的从左到右流动的输出视频（使用反弹方法）：链接。该视频显示了使用亮度作为大小和颜色作为方向的动量密度。我不记得它是使用 BGK 还是 TRT 生成的，但我在使用两个碰撞运算符时都遇到了类似的故障。我正在使用带有[[1, 1, 1], [1, 0, 1], [1, 1, 1]]模板的卷积来实现流式处理步骤。我尝试按照格子玻尔兹曼方法中描述的步骤：原理与实践在 BGK 和 TRT 模拟中获得稳定性，但没有任何运气。我在做一些明显错误的事情吗？