An auto-adaptive approximate Riemann solver for non-linear Euler equations

We present a new HLL-type approximate Riemann solver that aims at capturing any isolated discontinuity without necessitating extensive characteristic analysis of governing partial differential equations. This property is especially attractive for complex hyperbolic systems with more than two equations. Following Linde’s approach [6], we introduce a generic middle wave into the classical two-state HLL solver. The property of this third wave is typified by the way of a “strength indicator” that is derived from polynomial considerations. The polynomial that constitutes the basis of the procedure is made non-oscillatory by an adapted fourth-order WENO algorithm (CWENO4). This algorithm makes it possible to derive an expression for the strength indicator. According to the size of this latter parameter, the resulting solver (HLL-RH), either computes the multi-dimensional Rankine-Hugoniot equations if an isolated discontinuity appears in the Riemann fan, or asymptotically tends towards the two-state HLL solver if the solution is locally smooth. The asymptotic version of the HLL-RH solver is demonstrated to be positively conservative and entropy satisfying in its first-order multi-dimensional form provided that a relevant and not too restrictive CFL condition is considered; specific limitations of the conservative increments of the numerical solution and a suited entropy condition enable to maintain these properties in its high-order version.With a monotonicity-preserving algorithm for the time integration, the numerical method so generated, is third order in time and fourth-order accurate in space for the smooth part of the solution; moreover, the scheme is stable and accurate when capturing a shock wave, whatever the complexity of the underlying differential system.Extensive numerical tests for the one- and two-dimensional Euler equation of gas dynamics and comparisons with classical Godunov-type methods help to point out the potentialities and insufficiencies of the method.

► The velocity of an isolated discontinuity is accurately computed for problems of aerodynamics.
► The numerical procedure enables to compute complex flows while maintaining a high degree of accuracy.
► A “strength indicator” characterizes this procedure.