Navier-Stokes with Mimetic Methods:2D Lock exchange test case us...

单一的、广泛评论的文件，假设密度和温度变化很小 (Boussinesq)，它可以求解 NS 方程。 我选择了模拟方法（通过 MOLE）和明确的前向时间方案。 更详细和精确的时间离散化方案对这种特定场景的贡献很小。 通过使用 Leapfrog 等辛方法计算瞬态项，可以避免数值扩散。 对于锐利的Kelvin-Helmholtz滚滚，使用dx == 0.0625，请记住对于二维对流问题，在dt上施加CFL条件。 图片是使用 m = 1600、n = 320 和 dt = 0.1 获得的 最耗时的部分（如预期的那样）是压力场的计算，拉普拉斯矩阵高度稀疏但不是正定的。