A Modified Roe Scheme and Stability Analysis
-
摘要: 低耗散的激波捕捉方法,包括流行的Roe格式,在计算多维强激波问题时会遭遇激波不稳定现象的困扰,这会严重影响格式对于高超声速流动问题的精确模拟.对Roe格式进行小扰动分析,结果表明:激波面纵向所有物理量的扰动均会衰减,而横向的密度扰动和剪切速度扰动不会衰减.在横向数值通量上增加与熵波和剪切波相对应的黏性来抑制Roe格式不稳定现象的发生.为了防止不合适的黏性影响格式对于接触间断和剪切层的分辨率,定义两个开关函数,使得黏性仅仅添加在激波层亚声速区的横向数值通量上.数值测试的结果表明:改进的Roe格式不仅保留了原始Roe格式高分辨率的优点,而且具有更好的鲁棒性,消除了激波不稳定现象.Abstract: Low dissipation shock-capturing methods, including the popular Roe scheme, will encounter the shock instability phenomenon in the computation of the multidimensional strong shock wave problems. This will seriously affect the schemes’ accurate simulation of the hypersonic flow problems. The small perturbation analysis of the Roe scheme was carried out. The results show that, all perturbations in the longitudinal direction of the shock front are damped, but the perturbations of density and shear velocity in the transverse direction are undamped. The viscosities corresponding to the entropy wave and shear wave were added to the flux transverse to the shock front to suppress instability of the Roe scheme. To prevent the improper viscosity from influencing the resolution of contact discontinuity and shear layers, 2 switching functions were defined, so that the viscosity was only added to the transverse flux in the subsonic region of the shock layer. The numerical tests show that, the modified Roe scheme not only retains the merit of high resolution of the original Roe scheme, but also has better robustness and eliminates the shock wave instability.
-
Key words:
- hypersonic flow /
- Roe scheme /
- small perturbation analysis /
- numerical viscosity /
- shock instability
-
[1] ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics,1981,43(2): 357-372. [2] TORO E F, SPRUCE M, SPEARES W. Restoration of the contact surface in the HLL Riemann solver[J]. Shock Waves,1994,4(1): 25-34. [3] QUIRK J J. A contribution to the great Riemann solver debate[J]. International Journal for Numerical Methods in Fluids,1994,18(6): 555-574. [4] GRESSIER J, MOSCHETTA J M. Robustness versus accuracy in shock-wave computations[J]. International Journal for Numerical Methods in Fluids,2000,33(3): 313-332. [5] SHEN Z J, YAN W, YUAN G W. A robust HLLC-type Riemann solver for strong shock[J]. Journal of Computational Physics,2016,309: 185-206. [6] XIE W J, LI H, TIAN Z Y, et al. A low diffusion flux splitting method for inviscid compressible flows[J]. Computers & Fluids,2015,112(2): 83-93. [7] LIOU M S. Mass flux scheme and connection to shock instability[J]. Journal of Computational Physics,2000,160(2): 623-648. [8] LIOU M S. A sequel to AUSM: AUSM+[J]. Journal of Computational Physics,1996,129(2): 364-382. [9] XU K, LI Z W. Dissipative mechanism in Godunov-type schemes[J]. International Journal of Numerical Methods in Fluids,2001,37(1): 1-22. [10] DUMBSER M, MOSCHETTA J M, GRESSIER J. A matrix stability analysis of the carbuncle phenomenon[J]. Journal of Computational Physics,2004,197(2): 647-670. [11] CHEN S S, YAN C, LIN B X, et al. Affordable shock-stable item for Godunov-type schemes against carbuncle phenomenon[J]. Journal of Computational Physics,2018,373: 662-672. [12] XIE W J, LI W, LI H, et al. On numerical instabilities of Godunov-type schemes for strong shocks[J]. Journal of Computational Physics,2017,350: 607-637. [13] SIMON S, MANDAL J C. A cure for numerical shock instability in HLLC Riemann solver using antidiffusion control[J]. Computers & Fluids,2018,174: 144-166. [14] WU H, SHEN L J, SHEN Z J. A hybrid numerical method to cure numerical shock instability[J]. Communications in Computational Physics,2010,8: 1264-1271. [15] REN Y X. A robust shock-capturing scheme based on rotated Riemann solvers[J]. Computers & Fluids,2003,32(10): 1379-403. [16] NISHIKAWA H, KITAMURA K. Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers[J]. Journal of Computational Physics,2008,227(4): 2560-2581. [17] PEERY K M, IMLAY S T. Blunt-body flow simulations[C]//24th Joint Propulsion Conference . Boston, MA, 1988. [18] FLEISCHMANN N, ADAMI S, HU X Y, et al. A low dissipation method to cure the grid-aligned shock instability[J]. Journal of Computational Physics,2020,401: 109004. [19] KIM S, KIM C, RHO O H, et al. Cures for the shock instability: development of a shock-stable Roe scheme[J]. Journal of Computational Physics,2003,185(2): 342-374. [20] RODIONOV A. Artificial viscosity in Godunov-type schemes to cure the carbuncle phenomenon[J]. Journal of Computational Physics,2017,345: 308-329. [21] RODIONOV A. Artificial viscosity to cure the shock instability in high-order Godunov-type schemes[J]. 〖JP〗 Computers & Fluids,2019,190: 77-97. [22] RODIONOV A. Artificial viscosity to cure the carbuncle phenomenon: the three-dimensional case[J]. Journal of Computational Physics,2018,361: 50-55. [23] WOODWARD P, COLELLA P. The numerical simulation of two-dimensional fluid flow with strong shocks[J]. Journal of Computational Physics,1984,54(1): 115-173. [24] HUANG K B, WU H, YU H, et al. Cures for numerical shock instability in HLLC solver[J]. International Journal of Numerical Methods in Fluids,2011,65(9): 1026-1038.
点击查看大图
计量
- 文章访问数: 4241
- HTML全文浏览量: 265
- PDF下载量: 339
- 被引次数: 0