我正在模拟雷诺数为 500 的气缸上的不可压缩流动。我正在使用压力校正方法求解纳维斯托克斯方程。我的解决方案在一段时间后（大约 5 秒）变得不稳定。

我尝试改进我的网格，步长（0.05）（确保我的 CFL < 1，即使我使用的是隐式方法）

我的边界条件、网格和不稳定结果显示在附图中。域大约是圆柱体直径的 25 倍。

我尝试过模拟这个问题 O 网格（几乎立即变得不稳定）。

以下链接包含边界条件和结果的图片。

如果有人能分享他们对这个问题的想法/经验，我将不胜感激。非常感谢。

编辑：

为打字错误道歉：

我正在使用以下边界条件：Neumann 边界

在狄利克雷边界上

编辑：

我在狄利克雷边界周围的节点上应用了速度边界条件。此外，右上角和右下角节点是速度为 1 的狄利克雷边界。

之后，我更深入地研究了模拟结果，我注意到不稳定开始在流入/流出交界处蔓延。