用极坐标中的冲击波求解可压缩无粘性欧拉方程

计算科学 双曲-pde
2021-12-17 14:21:17

在过去的几周里,我一直在尝试调整 Lax-Wendroff 或一些类似的极坐标方案。这个过程很复杂,因为我无法找到所用方案的逐步推导,所以我不得不猜测笛卡尔坐标的过程,并对极坐标重复同样的过程。

我设法为 MacCormack 的方案做到了,即使在稍微平稳、无冲击波的测试中它也爆炸了。

由于我的系统的几何形状以及控制它的其他方程(我试图将流体动力学添加到二维等离子体放电模型中),我试图坚持极坐标。然而,现在我很确定我无法为极坐标中的欧拉方程找到现成的方案。自己推导出一个似乎也不可信。

因此我的问题出现了。

对于极坐标中具有急剧梯度的无粘性可压缩欧拉方程,是否有现成的稳健求解器?

如果没有,我如何尝试推导出一个?我应该尝试推导出哪个方案?

尝试坚持极坐标是否可行?或者我应该在笛卡尔坐标中求解欧拉方程,然后为模型的其余部分插值?如果是这样,我应该在边界上使用带有裁剪单元格的统一网格,还是应该使用一些有限体积的方法?

如果我应该远离极坐标,我应该使用哪种方案?有很多方案,有不同类型的纠错机制等 - 但是有经验的人可以提出一些建议吗?

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