求做空气动力学的大佬解答两个问题
-
最近要做一些跨临界的问题,但是在跨临界中,如下图,由于伪沸腾的影响,流体的热力学性质会在伪沸腾区呈现剧烈的变化,形成类似于激波的区域。但是在亚临界区域,压力对密度的变化极为敏感,所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如:求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐,且由于通过压力反推能量会引起能量不守恒。主要有两个问题,非常感谢您的解答:
- 求解压力传播方程的那个算法很类似simple算法的pressure based方法,为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法,更倾向于TVD RK这种显示迭代?
- 显式格式很适合GPU加速计算,想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算,截断误差是否会对计算结果产生比较大的影响?