Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

SUMMARY This paper reports numerical convergence study for simulations of steady shock-induced combustion problems with high-resolution shock-capturing schemes. Five typical schemes are used: the Roe flux-based monotone upstream-centered scheme for conservation laws (MUSCL) and weighted essentially non-oscillatory (WENO) schemes, the Lax–Friedrichs splitting-based non-oscillatory no-free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension-by-dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second-order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth-order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady-state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.

[1]  Seokkwan Yoon,et al.  Numerical study of chemically reacting flows using a lower-upper symmetric successive overrelaxation scheme , 1989 .

[2]  B C Khoo,et al.  The cellular structure of a two-dimensional H2/O2/Ar detonation wave , 2004 .

[3]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[4]  Guang-Shan Jiang,et al.  A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics , 1999 .

[5]  Sin-Chung Chang,et al.  Regular Article: The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws , 1999 .

[6]  Wei Shyy,et al.  Combustion and Turbulence Effects in Hypersonic Projectile Flows , 1997 .

[7]  C. J. Schexnayder,et al.  Influence of Chemical Kinetics and Unmixedness on Burning in Supersonic Hydrogen Flames , 1980 .

[8]  Jianxian Qiu,et al.  Trigonometric WENO Schemes for Hyperbolic Conservation Laws and Highly Oscillatory Problems , 2010 .

[9]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[10]  Hanxin Zhang,et al.  NND Schemes and Their Applications to Numerical Simulation of Two- and Three-Dimensional Flows , 1991 .

[11]  Tariq D. Aslam,et al.  High resolution numerical simulation of ideal and non-ideal compressible reacting flows with embedded internal boundaries , 1997 .

[12]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[13]  Chi-Wang Shu,et al.  High order finite difference methods with subcell resolution for advection equations with stiff source terms , 2012, J. Comput. Phys..

[14]  A. Gross,et al.  High-Order-Accurate Numerical Method for Complex Flows , 2008 .

[15]  Gang Wang,et al.  An improved CE/SE scheme for numerical simulation of gaseous and two-phase detonations , 2010 .

[16]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[17]  Wang Chi-Shu,et al.  Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws , 1997 .

[18]  Paul E. Dimotakis,et al.  An adaptive Lagrangian method for computing 1D reacting and non-reacting flows , 1993 .

[19]  Tao Tang,et al.  Resolving the shock-induced combustion by an adaptive mesh redistribution method , 2007, J. Comput. Phys..

[20]  De-Fu Zhang,et al.  A Bricklaying Heuristic Algorithm for the Orthogonal Rectangular Packing Problem: A Bricklaying Heuristic Algorithm for the Orthogonal Rectangular Packing Problem , 2009 .

[21]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[22]  Adam P. Bruckner,et al.  Numerical simulation of hypervelocity projectiles in detonable gases , 1991 .

[23]  Robert J. Kee,et al.  CHEMKIN-III: A FORTRAN chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics , 1996 .

[24]  Cheatwood F. McNeil,et al.  User''s Manual for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) , 1996 .

[25]  D. Scott Stewart,et al.  The Dynamics of Detonation in Explosive Systems , 2007 .

[26]  H. Lehr,et al.  Experiments on Shock-Induced Combustion , 1972 .

[27]  Sergio Pirozzoli,et al.  Numerical Methods for High-Speed Flows , 2011 .

[28]  Ralf Deiterding,et al.  An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry , 2011, J. Comput. Phys..

[29]  Kozo Fujii,et al.  Flow features of shock-induced combustion around projectile traveling at hypervelocities , 1993 .

[30]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[31]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[32]  J. Mcvey,et al.  Mechanism of Instabilities of Exothermic Hypersonic Blunt-Body Flows , 1971 .

[33]  A. Jameson,et al.  Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations , 1988 .

[34]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[35]  J. Ryan,et al.  A Runge-Kutta discontinuous Galerkin approach to solve reactive flows: The hyperbolic operator , 2011, J. Comput. Phys..

[36]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[37]  C. Lewis,et al.  A new technique for low-to-high altitude predictions of ablative hypersonic flowfields , 1991 .

[38]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[39]  Richard Sanders,et al.  Multidimensional Dissipation for Upwind Schemes , 1998 .

[40]  Richard A. Thompson,et al.  A review of reaction rates and thermodynamic and transport properties for the 11-species air model for chemical and thermal nonequilibrium calculations to 30000 K , 1989 .

[41]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[42]  Francis A. Greene,et al.  Viscous equilibrium computations using program LAURA , 1991 .

[43]  Joseph E. Shepherd,et al.  Detonation Waves and Propulsion , 1994 .

[44]  Meng-Sing Liou,et al.  Ten Years in the Making: AUSM-Family , 2001 .