On the sonic point glitch

This paper presents theoretical and numerical analyses of the sonic point glitch based on some numerical schemes for the Burgers' equation and the Euler equations in fluid mechanics. The sonic glitch is formed in the sonic rarefaction fan. It has no any direct connection with the violation of the entropy condition or the size of numerical viscosity of a finite-difference scheme. Our results show that it is mainly coming from a disparity in wave speeds across the sonic point. If numerical viscosity depends on the characteristic direction, then the disparity may be formed between the numerical and physical wave speeds around the sonic point, and triggers the sonic wiggle in the numerical solution. We also find that the initial data reconstruction technique of van Leer can effectively eliminate the flaw around the sonic point for the Burgers' equation. Some other possible cures are also suggested.

[1]  Jean-Marc Moschetta,et al.  A Robust Low Diffusive Kinetic Scheme for the Navier-Stokes/Euler Equations , 1997 .

[2]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[3]  Philip L. Roe,et al.  Sonic Flux Formulae , 1992, SIAM J. Sci. Comput..

[4]  S. F. Davis,et al.  A simplified TVD finite difference sheme via artificial viscousity , 1987 .

[5]  J. Quirk.  A Contribution to the Great Riemann Solver Debate , 1994 .

[6]  Huazhong Tang,et al.  Pseudoparticle representation and positivity analysis of explicit and implicit Steger-Warming FVS schemes , 2001 .

[7]  D. Pullin,et al.  Direct simulation methods for compressible inviscid ideal-gas flow , 1980 .

[8]  Stanley Osher,et al.  Upwind schemes and boundary conditions with applications to Euler equations in general geometries , 1983 .

[9]  Bram van Leer,et al.  On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe , 1984 .

[10]  B. Leer,et al.  Flux-vector splitting for the Euler equations , 1997 .

[11]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[12]  R. H. Sanders,et al.  The possible relation of the 3-kiloparsec arm to explosions in the galactic nucleus , 1974 .

[13]  J. Steger,et al.  Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods , 1981 .

[14]  Kun Xu,et al.  Gas-kinetic schemes for unsteady compressible flow simulations , 1998 .

[15]  S. Osher,et al.  One-sided difference approximations for nonlinear conservation laws , 1981 .

[16]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[17]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[18]  J. Anderson,et al.  Modern Compressible Flow: With Historical Perspective , 1982 .

[19]  Charles Hirsch,et al.  Numerical computation of internal & external flows: fundamentals of numerical discretization , 1988 .

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

[21]  Jean-Marc Moschetta,et al.  A Cure for the Sonic Point Glitch , 2000 .

[22]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

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

[24]  Tao Tang,et al.  Adaptive Mesh Methods for One- and Two-Dimensional Hyperbolic Conservation Laws , 2003, SIAM J. Numer. Anal..

[25]  Wai How Hui,et al.  ON CONTACT OVERHEATING AND OTHER COMPUTATIONAL DIFFICULTIES OF SHOCK-CAPTURING METHODS , 2002 .

[26]  M. Liou,et al.  A New Flux Splitting Scheme , 1993 .