×

A cure for numerical shock instability in HLLC Riemann solver using antidiffusion control. (English) Zbl 1410.76334

Summary: Various forms of numerical shock instabilities are known to plague many contact and shear preserving approximate Riemann solvers, including the popular Harten-Lax-van Leer with contact (HLLC) scheme, during high speed flow simulations governed by the Euler system of equations. In this paper, we propose a simple and inexpensive novel strategy to prevent the HLLC scheme from developing such spurious solutions without compromising on its linear wave resolution ability. The cure is primarily based on a reinterpretation of the HLLC scheme as a combination of its well-known diffusive counterpart, the HLL scheme, and an antidiffusive term responsible for its accuracy on linear wavefields. In our study, a linear analysis of this alternate form indicates that shock instability in the HLLC scheme could be triggered due to the unwanted activation of the antidiffusive term appearing in its mass and interface-normal momentum flux component discretizations on interfaces that are not aligned with the shock front. This inadvertent activation results in weakening of the favourable dissipation provided by its inherent HLL scheme and causes unphysical mass flux variations along the shock front. To mitigate this, we propose a modified HLLC scheme that employs a simple differentiable pressure-ratio based multidimensional shock sensor to achieve smooth control of these critical antidiffusive terms near shocks. Using a linear perturbation analysis and a matrix based stability analysis, we establish that the resulting scheme, called HLLC-ADC (anti-diffusion control), is shock stable over a wide range of freestream Mach numbers. Results from standard numerical test cases demonstrate that the HLLC-ADC scheme is indeed free from the most common manifestations of shock instability including the carbuncle phenomenon without significant loss of accuracy on shear dominated viscous flows.

MSC:

76M25 Other numerical methods (fluid mechanics) (MSC2010)
65M25 Numerical aspects of the method of characteristics for initial value and initial-boundary value problems involving PDEs
76L05 Shock waves and blast waves in fluid mechanics

Software:

HLLE

References:

[1] Roe, P., Approximate Riemann solvers, parameter vectors, and difference schemes, J Comput Phys, 43, 357-372, (1981) · Zbl 0474.65066
[2] Harten, A., High resolution schemes for hyperbolic conservation laws, J Comput Phys, 49, 3, 357-393, (1983) · Zbl 0565.65050
[3] Toro, E. F.; Spruce, M.; Speares, W., Restoration of the contact surface in the HLL-Riemann solver, Shock Waves, 4, 25-34, (1994) · Zbl 0811.76053
[4] Einfeldt, B.; Munz, C. D.; Roe, P. L.; Sjogreen, B., On gudonov type methods near low densities, J Comput Phys, 92, 273-295, (1991) · Zbl 0709.76102
[5] Liou, M. S.; Steffen, C. J., A new flux splitting scheme, J Comput Phys, 107, 23-39, (1993) · Zbl 0779.76056
[6] Park, S. H.; Kwon, J. H., On the dissipation mechanism of Godunov-type schemes, 188, (2003) · Zbl 1022.76037
[7] Dumbser, M.; Balsara, D. S., A new efficient formulation of the hllem Riemann solver for general conservative and non-conservative hyperbolic systems, J Comput Phys, 304, 275-319, (2016) · Zbl 1349.76603
[8] Xie, W.; Li, W.; Li, H.; Tian, Z.; Pan, S., On numerical instabilities of Godunov-type schemes for strong shocks, J Comput Phys, 350, 607-637, (2017) · Zbl 1380.76073
[9] Einfeldt, B., On Godunov-type methods for gas dynamics, SIAM J Numer Anal, 25, 2, 294-318, (1988) · Zbl 0642.76088
[10] Shen, Z.; Yan, W.; Yuan, G., A robust hllc-type Riemann solver for strong shock, J Comput Phys, 309, 185-206, (2016) · Zbl 1351.76043
[11] Mandal, J. C.; Rao, S. P., High resolution finite volume computations on unstructured grids using solution dependent weighted least squares gradients, Comput Fluids, 44, 1, 23-31, (2011) · Zbl 1271.76189
[12] van Leer, B.; Thomas, J. L.; Roe, P. L.; Newsome, R. W., A comparison of numerical flux formulas for the Euler and Navier-Stokes equations, 36-41, (1987), Am Instit Aeronaut Astronaut Inc, AIAA
[13] Henderson, S.; Menart, J., Grid study on blunt bodies with the carbuncle phenomenon, 39th AIAA Thermophysics conference, Miami, Florida, (2007)
[14] Dumbser, M.; Moschetta, J. M.; Gressier, J., A matrix stability analysis of the carbuncle phenomenon, J Comput Phys, 197, 647-670, (2004) · Zbl 1079.76607
[15] Chauvat, Y.; Moschetta, J. M.; Gressier, J., Shock wave numerical structure and the carbuncle phenomenon, Int J Numer Methods Fluids, 47, 903-909, (2005) · Zbl 1134.76372
[16] Gressiera, J.; Moschetta, J.-M., Robustness versus accuracy in shock-wave computations, Int J Numer Methods Fluids, 33, 313-332, (2000) · Zbl 0980.76072
[17] Peery, K. M.; Imlay, S. T., Blunt-body flow simulations, AIAA Paper 88-2904, (1988)
[18] Quirk, J. J., A contribution to the great Riemann solver debate, Int J Numer Methods Fluids, 18, 555-574, (1994) · Zbl 0794.76061
[19] Sanders, R.; Morano, E.; Druguet, M. C., Multidimensional dissipation for upwind schemes:stability and applications to gas dynamics, J Comput Phys, 145, 511-537, (1998) · Zbl 0924.76076
[20] Shen, Z.; Yan, W.; Yuan, G., A stability analysis of hybrid schemes to cure shock instability, Commun Comput Phys, 15, 1320-1342, (2014) · Zbl 1373.76144
[21] Pandolfi, M.; DAmbrosio, D., Numerical instabilities in upwind methods: analysis and cures for the carbuncle phenomenon, J Comput Phys, 166, 271-301, (2001) · Zbl 0990.76051
[22] Ismail, F., Toward a reliable prediction of shocks in hypersonic flow: Resolving carbuncles with entropy and vorticity control, Ph.D. thesis, (2006), University of Michigan
[23] Lin, H.-C., Dissipation additions to flux-difference splitting, J Comput Phys, 117, 1, 20-27, (1995) · Zbl 0836.76061
[24] Ren, Y.-X., A robust shock-capturing scheme based on rotated Riemann solvers, Comput Fluids, 32, 1379-1403, (2003) · Zbl 1034.76035
[25] Nishikawa, H.; Kitamura, K., Very simple carbuncle free boundary layer resolving rotated hybrid Riemann solvers, J Comput Phys, 227, 2560-2581, (2008) · Zbl 1388.76185
[26] Zhang, F.; Liu, J.; Chen, B.; Zhong, W., Evaluation of rotated upwind schemes for contact discontinuity and strong shock, Comput Fluids, 134, 11-22, (2016) · Zbl 1390.76541
[27] Huang, K.; Wu, H.; Yu, H.; Yan, D., Cures for numerical shock instability in hllc solver, Int J Numer Methods Fluids, 65, 1026-1038, (2011) · Zbl 1428.76120
[28] Wu, H.; Shen, L.; Shen, Z., A hybrid numerical method to cure numerical shock instability, Commun Comput Phys, 8, 1264-1271, (2010) · Zbl 1364.76163
[29] Wang, D.; Deng, X.; Wang, G.; Dong, Y., Developing a hybrid flux function suitable for hypersonic flow simulation with high-order methods, Int J Numer Methods Fluids, 81, 5, 309-327, (2016)
[30] Kim, S.; Lee, B.; Lee, H.; Jeung, I.-S., Robust hllc Riemann solver with weighted average flux scheme for strong shock, J Comput Phys, 228, 7634-7642, (2009) · Zbl 1391.76556
[31] Kim, S. D.; Lee, B. J.; Lee, H. J.; Jeung, I.-S.; Choi, J.-Y., Realization of contact resolving approximate Riemann solvers for strong shock and expansion flows, Int J Numer Methods Fluids, 62, 1107-1133, (2010) · Zbl 1423.76340
[32] Zhang, F.; Liu, J.; Chen, B.; Zhong, W., A robust low-dissipation ausm-family scheme for numerical shock stability on unstructured grids, Int J Numer Methods Fluids, 84, 135-151, (2017)
[33] Kim, K.; Lee, J.; Rho, O., An improvement of ausm schemes by introducing the pressure-based weight functions, Comput Fluids, 27, 311-346, (1998) · Zbl 0964.76064
[34] Wada, Y.; Liou, M.-S., An accurate and robust flux splitting scheme for shock and contact discontinuities, SIAM J Sci Comput, 18, 633-657, (1997) · Zbl 0879.76064
[35] Zoltak, J.; Drikakis, D., Hybrid upwind methods for the simulation of unsteady shock-wave diffraction over a cylinder, Comput Method Appl M, 162, 97, 165-185, (1998) · Zbl 0946.76068
[36] Liou, M.-S., Mass flux schemes and connection to shock instability, J Comput Phys, 160, 623-648, (2000) · Zbl 0967.76062
[37] Harten, A.; Lax, P. D.; van Leer, B., On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev, 25, 35-61, (1983) · Zbl 0565.65051
[38] Toro, E.; Cendon, V., Flux splitting schemes for the Euler equations, Comput Fluids, 70, 1-12, (2012) · Zbl 1365.76243
[39] Batten, P.; Clarke, N.; Lambert, C.; Causon, D. M., On the choice of wavespeeds for the hllc Riemann solver, SIAM J Sci Comput, 18, 6, 1553-1570, (1997) · Zbl 0992.65088
[40] Phongthanapanich, S., A modified multidimensional dissipation technique for ausm+ on triangular grids, Int J Comput Fluid Dyn, 29, 1, 1-11, (2015) · Zbl 07514811
[41] Blazek, J., Computational fluid dynamics: principles and applications (second edition), (2005), Elsevier Science
[42] Barth, T.; Jesperson, D., The design and application of upwind schemes on unstructured meshes, (1989), Am Instit Aeronaut Astronaut Inc, AIAA
[43] Gottlieb, S.; Shu, C.-W.; Tadmor, E., Strong stability-preserving high-order time discretization, SIAM Rev, 43, 1, 89-112, (2001) · Zbl 0967.65098
[44] Ohwada, T.; Adachi, R.; Xu, K.; Luo, J., On the remedy against shock anomalies in kinetic schemes, J Comput Phys, 255, 106-129, (2013) · Zbl 1349.76722
[45] Woodward, P.; Colella, P., The numerical simulation of two-dimensional fluid flow with strong shocks, J Comput Phys, 54, 1, 115-173, (1984) · Zbl 0573.76057
[46] Coirier, W., An adaptively refined cartesian, cell based scheme for the euler and navier-stokes equations, Ph.D. thesis, (1994), University of Michigan
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.