Abstract
The Generalized Riemann Problem (GRP) for a nonlinear hyperbolic system of m balance laws (or alternatively “quasi-conservative” laws) in one space dimension is now well-known and can be formulated as follows: Given initial-data which are analytic on two sides of a discontinuity, determine the time evolution of the solution at the discontinuity. In particular, the GRP numerical scheme (second-order high resolution) is based on an analytical evaluation of the first time derivative. It turns out that this derivative depends only on the first-order spatial derivatives, hence the initial data can be taken as piecewise linear. The analytical solution is readily obtained for a single equation (m = 1) and, more generally, if the system is endowed with a complete (coordinate) set of Riemann invariants. In this case it can be “diagonalized” and reduced to the scalar case. However, most systems with m > 2 do not admit such a set of Riemann invariants. This paper introduces a generalization of this concept: weakly coupled systems (WCS). Such systems have only “partial set” of Riemann invariants, but these sets are weakly coupled in a way which enables a “diagonalized” treatment of the GRP. An important example of a WCS is the Euler system of compressible, nonisentropic fluid flow (m = 3). The solution of the GRP discussed here is based on a careful analysis of rarefaction waves. A “propagation of singularities” argument is applied to appropriate Riemann invariants across the rarefaction fan. It serves to “rotate” initial spatial slopes into “time derivative”. In particular, the case of a “sonic point” is incorporated easily into the general treatment. A GRP scheme based on this solution is derived, and several numerical examples are presented. Special attention is given to the “acoustic approximation” of the analytical solution. It can be viewed as a proper linearization (different from the approach of Roe) of the nonlinear system. The resulting numerical scheme is the simplest (second-order, high-resolution) generalization of the Godunov scheme.
Similar content being viewed by others
References
Ben-Artzi M. and Falcovitz J. (1984). A second-order Godunov-type scheme for compressible fluid dynamics. J. Comput. Phys. 55(1): 1–32
Ben-Artzi M. and Falcovitz J. (1986). An upwind second-order scheme for compressible duct flows. SIAM J. Sci. Statist. Comput. 7(3): 744–768
Ben-Artzi M. (1989). The generalized Riemann problem for reactive flows. J. Comput. Phys. 81(1): 70–101
Ben-Artzi M. and Falcovitz J. (2003). Generalized Riemann Problems in Computational Gas Dynamics. Cambridge University Press, Cambridge
Ben-Artzi M., Li J. and Warnecke G. (2006). A direct Eulerian GRP scheme for compressible fluid flows. J. Comp. Phys. 218: 19–43
Bouchut F., Le Sommer J. and Zeitlin V. (2004). Frontal geostrophic adjustment and nonlinear wave phenomena in one dimensional rotating shallow water. Part 2: High resolution numerical simulations. J. Fluid Mech. 514: 35–63
Bourgeade A., LeFloch Ph. and Raviart P.-A. (1989). An asymptotic expansion for the solution of the generalized Riemann problem. II. Application to the equations of gas dynamics. Ann. Inst. H. Poincar Anal. Non Linaire 6(6): 437–480
Dafermos C.M. (2000). Hyperbolic Conservation Laws in Continuum Physics. Springer, Heidelberg
Godunov S.K. (1959). A finite difference method for the numerical computation and disontinuous solutions of the equations of fluid dynamics. Mat. Sb. 47: 271–295
Godlewski, E., Raviart, P.-A.: Numerical approximation of hyperbolic systems of conservation laws. Appl. Math. Sci. 118, Springer, Heidelberg (1996)
Jin S. (2001). A steady-state capturing method for hyperbolic systems with geometrical source terms. M2 AN Math. Model. Numer. Anal. 35(4): 631–645
Leer B. (1979). Towards the ultimate conservative difference scheme, V. A second-order sequel to Godunov’s method. J. Comput. Phys. 32(1): 101–136
Li J. and Chen G. (2006). The generalized Riemann problem method for the shallow water equations with bottom topography. Int. J. Numer. Methods Eng. 65(6): 834–862
Li, T.T.: Global Classical Solutions for Quasilinear Hyperbolic Systems. Research in Applied Mathematics. Wiley, Chichester/Masson, Paris (1994)
LeFloch Ph. and Raviart P.-A. (1988). An asymptotic expansion for the solution of the generalized Riemann problem. I. General theory. Ann. Inst. H. Poincar Anal. Non Linaire 5(2): 179–207
Smoller J. (1994) Shock Waves and Reaction-Diffusion Equations, 2nd edn, vol. 258. Springer, New York
Sweby P.K. (1984). High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Anal. 21(5): 995–1011
Toro E. (1997). Riemann solvers and numerical methods for fluid dynamics. A practical introduction. Springer, Berlin
Toro E. (2006). Derivative Riemann solvers for systems of conservation laws and ADER methods. J. Comp. Phys. 212: 150–165
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Ben-Artzi, M., Li, J. Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem. Numer. Math. 106, 369–425 (2007). https://doi.org/10.1007/s00211-007-0069-y
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-007-0069-y