Affiliation:
1. NLCFD, School of Aeronautic Science and Engineering Beihang University Beijing People's Republic of China
2. School of Aeronautics Northwestern Polytechnical University Xi'an People's Republic of China
3. China Aerodynamic Research and Development Center Mianyang People's Republic of China
Abstract
AbstractFollowing the solution formula method given in Dong et al. (High order discontinuities decomposition entropy condition schemes for Euler equations. CFD J. 2002;10(4): 448–457), this article studies a type of one‐step fully‐discrete scheme, and constructs a third‐order scheme which is written into a compact form via a new limiter. The highlights of this study and advantages of new third‐order scheme are as follows: ① We proposed a very simple new methodology of constructing one‐step, consistent high‐order and non‐oscillation schemes that do not rely on Runge–Kutta method; ② We systematically studied new scheme's theoretical problems about entropy conditions, error analysis, and non‐oscillation conditions; ③ The new scheme achieves exact solution in linear cases and performing better in nonlinear cases when CFL → 1; ④ The new scheme is third order but high resolution with excellent shock‐capturing capacity which is comparable to fifth order WENO scheme; ⑤ CPU time of new scheme is only a quarter of WENO5 + RK3 under same computing condition; ⑥ For engineering applications, the new scheme is extended to multi‐dimensional Euler equations under curvilinear coordinates. Numerical experiments contain 1D scalar equation, 1D,2D,3D Euler equations. Accuracy tests are carried out using 1D linear scalar equation, 1D Burgers equation and 2D Euler equations and two sonic point tests are carried out to show the effect of entropy condition linearization. All tests are compared with results of WENO5 and finally indicate EC3 is cheaper in computational expense.