文件名称:Navier-Stokes with Mimetic Methods:2D Lock exchange test case using mimetic methods-matlab开发
文件大小:4KB
文件格式:ZIP
更新时间:2024-06-20 11:19:24
matlab
单一的、广泛评论的文件,假设密度和温度变化很小 (Boussinesq),它可以求解 NS 方程。 我选择了模拟方法(通过 MOLE)和明确的前向时间方案。 更详细和精确的时间离散化方案对这种特定场景的贡献很小。 通过使用 Leapfrog 等辛方法计算瞬态项,可以避免数值扩散。 对于锐利的Kelvin-Helmholtz滚滚,使用dx == 0.0625,请记住对于二维对流问题,在dt上施加CFL条件。 图片是使用 m = 1600、n = 320 和 dt = 0.1 获得的 最耗时的部分(如预期的那样)是压力场的计算,拉普拉斯矩阵高度稀疏但不是正定的。
【文件预览】:
lock_exchange.m.zip