单一,广泛评论的文件,解决了假设小密度和温度变化的NS方程(Boussinesq)。我选择了模拟方法(通过MOLE)和显式前向时间方案。更详细和精确的时间离散方案对这种特殊情况的贡献很小。采用辛方法如Leapfrog计算瞬态项可以避免数值扩散。
对于尖锐的开尔文-亥姆霍兹波,使用dx == 0.0625,记住在二维平流问题中CFL条件对dt施加的影响。图片由m = 1600, n = 320, dt = 0.1得到
最耗时的部分(如预期)是压力场的计算,拉普拉斯矩阵是高度稀疏的,但不是正定的。
引用作为
约翰尼·科尔比诺·德尔加多(2021)。Navier-Stokes与模拟方法(//www.tianjin-qmedu.com/matlabcentral/fileexchange/87402-navier-stokes-with-mimetic-methods), MATLAB中央文件交换。检索.