欧拉方程数值求解的高精度通量分裂方法

针对欧拉方程,提出一种将总能对流迎风和分压(E-CUSP)格式与加权本质无振荡(WENO)格式相耦合的新格式。在空间方向上,通过对低耗散E-CUSP格式的通量,采用高精度WENO格式进行重构;在时间方向上,使用4阶总变差递减(TVD)的Runge-Kutta方法进行推进,由此得到求解欧拉方程的高精度通量分裂方法。考虑E-CUSP格式与WENO重构进行耦合得到新格式,使其空间精度进一步提高。通过对激波管问题进行数值模拟发现,新的格式相对于E-CUSP格式对激波和接触间断捕捉的效果更加精准。数值结果表明:耦合得到的新格式具有更高的准确性和稳健性。

一种求解Euler方程的高阶格式

将总能对流迎风和分压(E-CUSP)格式与WENO-η-格式耦合给出一种新格式(记为E-CUSP-WENO-η),解决了一维激波管问题.数值模拟结果表明:新格式对接触间断和激波的捕捉更精确,其中带有高阶最优全局光滑因子的新格式(E-CUSP-WENO-η(τ7^(opt)))数值耗散最小;耦合的新格式可以更陡峭地捕捉到激波,且计算结果精度高,稳定性好.

一种求解欧拉方程的高精度数值格式

为了优化加权本质无振荡(weighted essentially non-oscillatory,WENO)格式在流体力学中的计算性能,在WENO-JS格式和WENO-M格式加权方法研究分析的基础上,基于广义不可分差分算子推导出新局部光滑因子,构造出WENO-NS格式。将总能对流迎风和分压(E-CUSP)格式与七阶WENO-NS格式耦合得到新格式。在空间层上,采用七阶WENO-NS格式重构低耗散E-CUSP格式中的通量;在时间层上,采用四阶总变差递减的龙格-库塔(total variation decreasing Runge-Kutta,TVD-Runge-Kutta)方法推进,对一维激波管问题进行了数值模拟。观察可得,耦合的E-CUSP-WENO-NS7对接触间断和激波的捕捉更为精确。结果分析表明,耦合的新格式可以更陡峭地捕捉到激波,计算结果精度高,稳定性能好。