我正在尝试使用 QUICK 对流方案来实现有限差分流体求解器(不可压缩、无粘性)。大多数情况下,我一直在关注诺里斯论文第 2 章 ( https://ses.library.usyd.edu.au/handle/2123/376 ) 的数学。
通过使用 Adams-Bashforth (AB2) 时间积分和 Leonard 的 ULTRA 限流器,使标量场平流稳定且不振荡。
但是,当将这种平流方案应用于速度场本身时,代码就会崩溃。特别是在域边界(使用开放/零压力边界条件)处会出现振荡。在一维中,振荡幅度似乎保持有限,但在二维中,发散模式出现在域角。
我尝试了多种边界单元处理方法,但均未成功。
是否有必要修改 QUICK/AB2 方案以实现速度自平流?
如果有帮助,我可以发布可视化。