五阶FD-WENO格式与二阶Godunov格式MUSCL的数值测试与比较硕士论文.pdfVIP

℄ FD- WENO  70-80  Godunov ICF  2 Godunov MUSCL Euler ℄ FD-WENO ℄ WENO5 MUSCL ℄ Burgers Riemann ℄ Richtmyer-Meshkov(RM Rayleigh-Taylor(RT) ℄ ℄ Sod Riemann Lax Riemann ℄ WENO5 MUSCL Peak Riemann Woodward-Colella WENO5 MUSCL 81 × 81 WENO5 t=10 161 × 161 MUSCL t=10 RM RT WENO5 MUSCL ℄ Riemann MUSCL Æ ICF FD-WENO WENO Godunov I ABSTRACT A numerical study is undertaken comparing the fifth-order finite difference weighted essentially non-oscillatory scheme (FD-WENO5) to the second-order Godunov scheme MUSCL (monotonic upstream-centered scheme for conserva- tion laws). For quantitative comparison purpose, these methods are tested on a series of problems whose true solutions are known or can be computed with high accuracy, such as linear advection problem, Sod Riemann, ”peak” Riemann problem, Woodward and Colella interacting shock wave problem and the two- dimensional periodic vortex problem. For qualitative comparison purpose, these methods are tested on Rayleigh-Taylor instability and Richtmyer-Meshkov insta- bility problems. The numerical results show that: for the Sod and Lax Riemann problems, high accuracy and high resolution results can be easily achieved by both schemes, but MUSCL is performed much faster than WENO5. For the ”Peak” Riemann problem and the Woodward and Colella interacting shock wave problem, when both methods are performed on same space mesh, we find that WENO5 is better than MUSCL around the shock and the rarefaction, in more de- tail, WENO5 has not only higher resolution but also smaller errors than MUSCL. For the two-dimensional periodic vortex problem, W