二维黎曼问题数值振荡

[Xia YC]

我在进行二维黎曼问题模拟时，随着时间步推进，温度云图如图所示，可以看到在间断位置出现了很多的振荡，这一般是什么原因呢，我的边界就是采用的外推边界，这应该不会导致反射吧

[李东岳]

云图看起来不太好看。最好用线图。
1D有震荡么

[Xia YC]

@李东岳 谢谢老师建议！然后我算了一下一维的，发现在间断位置也是存在振荡的。。。

[李东岳]

娜个求解器，还是你自己写的代码

[Xia YC]

是我自己写的代码，不过用的是DUGKS方法+Shakhov模型

[李东岳]

DUGKS我不懂了。
原生GKS我植入过，那个宏观变量太麻烦了我都要疯了。我一直琢磨大家植入GKS的时候，宏观通量都是这么复杂么

[Xia YC]

原生的GKS确实感觉挺麻烦，我弄明白里面的那些通量花了不少时间，不过DUGKS就简化很多了，宏观通量可以都直接从分布方程来。我以前一直用LBM的，感觉DUGKS共通的比较多也比较好理解没那么麻烦，就在用这个了，之前不可压问题复现都没问题，到可压这块就出问题了

[李东岳]

原来如此。我只看了GKS，给我搞疯了，后来就不看了。GKS那面跟矩方法也类似。都是从玻尔兹曼方程出来的，应该确实有共同点。

你要不要面插值的时候试试upwind，不过你应该已经试过了。

[Xia YC]

嗯嗯，还是谢谢老师了，我再比较比较文献和我程序实现里面的不同看看