Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks

In this paper, an entropy-consistent flux is developed, continuing from the work of the previous paper. To achieve entropy consistency, a second and third-order differential terms are added to the entropy-conservative flux. This new flux function is tested on several one dimensional problems and compared with the original Roe flux. The new flux function exactly preserves the stationary contact discontinuity and does not capture the unphysical rarefaction shock. For steady shock problems, the new flux predicts a slightly more diffused profile whereas for unsteady cases, the captured shock is very similar to those produced by the Roe- flux. The shock stability is also studied in one dimension. Unlike the original Roe flux, the new flux is completely stable which will provide as a candidate to combat multidimensional shock instability, particularly the carbuncle phenomenon.

[1]  Domenic D'Ambrosio,et al.  Numerical Instablilities in Upwind Methods: Analysis and Cures for the “Carbuncle” Phenomenon , 2001 .

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

[3]  Keiichi Kitamura,et al.  An Evaluation of Euler Fluxes for Hypersonic Flow Computations , 2007 .

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

[5]  R. LeVeque.  Approximate Riemann Solvers , 1992 .

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

[7]  Keiichi Kitamura,et al.  Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers , 2008, J. Comput. Phys..

[8]  Eitan Tadmor,et al.  COMPENSATED COMPACTNESS FOR 2D CONSERVATION LAWS , 2005 .

[9]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[10]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier—Stokes equations and the second law of thermodynamics , 1986 .

[11]  Farzad Ismail,et al.  Toward a reliable prediction of shocks in hypersonic flow: Resolving carbuncles with entropy and vorticity control , 2006 .

[12]  Eitan Tadmor,et al.  ENTROPY STABLE APPROXIMATIONS OF NAVIER-STOKES EQUATIONS WITH NO ARTIFICIAL NUMERICAL VISCOSITY , 2006 .

[13]  Michael Dumbser,et al.  A matrix stability analysis of the carbuncle phenomenon , 2004 .

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

[15]  P. Lax.  Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves , 1987 .

[16]  Keiichi Kitamura,et al.  Artificial Surface Tension to Stabilize Captured Shockwaves , 2008 .

[17]  S. Osher.  Riemann Solvers, the Entropy Condition, and Difference , 1984 .