Developing a hybrid flux function suitable for hypersonic flow simulation with high-order methods

In this paper, we develop a new hybrid Euler flux function based on Roe's flux difference scheme, which is free from shock instability and still preserves the accuracy and efficiency of Roe's flux scheme. For computational cost, only 5% extra CPU time is required compared with Roe's FDS. In hypersonic flow simulation with high-order methods, the hybrid flux function would automatically switch to the Rusanov flux function near shock waves to improve the robustness, and in smooth regions, Roe's FDS would be recovered so that the advantages of high-order methods can be maintained. Multidimensional dissipation is introduced to eliminate the adverse effects caused by flux function switching and further enhance the robustness of shock-capturing, especially when the shock waves are not aligned with grids. A series of tests shows that this new hybrid flux function with a high-order weighted compact nonlinear scheme is not only robust for shock-capturing but also accurate for hypersonic heat transfer prediction. Copyright (C) 2015 John Wiley & Sons, Ltd.