A Stability Analysis of Hybrid Schemes to Cure Shock Instability
Zhijun Shen 1*, Wei Yan 2, Guangwei Yuan 21 National Key Laboratory of Science and Technology on Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, China; Center for Applied Physics and Technology, HEDPS, Peking University, Beijing 100871, China.
2 National Key Laboratory of Science and Technology on Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, China.
Received 21 May 2013; Accepted (in revised version) 9 October 2013
Available online 14 February 2014
The carbuncle phenomenon has been regarded as a spurious solution produced by most of contact-preserving methods. The hybrid method of combining high resolution flux with more dissipative solver is an attractive attempt to cure this kind of non-physical phenomenon. In this paper, a matrix-based stability analysis for 2-D Euler equations is performed to explore the cause of instability of numerical schemes. By combining the Roe with HLL flux in different directions and different flux components, we give an interesting explanation to the linear numerical instability. Based on such analysis, some hybrid schemes are compared to illustrate different mechanisms in controlling shock instability. Numerical experiments are presented to verify our analysis results. The conclusion is that the scheme of restricting directly instability source is more stable than other hybrid schemes.AMS subject classifications: 35L65, 65M08, 76M12, 76L05
Notice: Undefined variable: pac in /var/www/html/readabs.php on line 165
Key words: Godunov methods, numerical shock instability, hybrid scheme.
Email: firstname.lastname@example.org (Z. Shen), email@example.com (W. Yan), firstname.lastname@example.org (G. Yuan)