QUICK迎风方案的速度自平流稳定性问题

计算科学 流体动力学
2021-11-27 11:31:50

我正在尝试使用 QUICK 对流方案来实现有限差分流体求解器(不可压缩、无粘性)。大多数情况下,我一直在关注诺里斯论文第 2 章 ( https://ses.library.usyd.edu.au/handle/2123/376 ) 的数学。

通过使用 Adams-Bashforth (AB2) 时间积分和 Leonard 的 ULTRA 限流器,使标量场平流稳定且不振荡。

但是,当将这种平流方案应用于速度场本身时,代码就会崩溃。特别是在域边界(使用开放/零压力边界条件)处会出现振荡。在一维中,振荡幅度似乎保持有限,但在二维中,发散模式出现在域角。

我尝试了多种边界单元处理方法,但均未成功。

是否有必要修改 QUICK/AB2 方案以实现速度自平流?

如果有帮助,我可以发布可视化。

0个回答
没有发现任何回复~