The analytical solution of two interesting hyperbolic problems as a test case for a finite volume method with a new grid refinement technique

• Published in: Journal of computational and applied mathematics 2008-05, Vol.214 (2), p.509-532
• Main Authors: Heineken, Wolfram, Kunik, Matthias
• Format: Article
• Language: English
• Subjects: Exact sciences and technology; Finite volume method; Global analysis, analysis on manifolds; Grid adaption; Lorentz transformation; Mathematical analysis; Mathematics; Measure and integration; Min-mod limiter; Numerical analysis; Numerical analysis. Scientific computation; Partial differential equations; Riemann problem; Scalar conservation law; Sciences and techniques of general use; Shock waves; Topology. Manifolds and cell complexes. Global analysis and analysis on manifolds; Utra-relativistic Euler equations
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/j.cam.2007.03.008

A finite volume method with grid adaption is applied to two hyperbolic problems: the ultra-relativistic Euler equations, and a scalar conservation law. Both problems are considered in two space dimensions and share the common feature of moving shock waves. In contrast to the classical Euler equations, the derivation of appropriate initial conditions for the ultra-relativistic Euler equations is a non-trivial problem that is solved using one-dimensional shock conditions and the Lorentz invariance of the system. The discretization of both problems is based on a finite volume method of second order in both space and time on a triangular grid. We introduce a variant of the min-mod limiter that avoids unphysical states for the Euler system. The grid is adapted during the integration process. The frequency of grid adaption is controlled automatically in order to guarantee a fine resolution of the moving shock fronts. We introduce the concept of “width refinement” which enlarges the width of strongly refined regions around the shock fronts; the optimal width is found by a numerical study. As a result we are able to improve efficiency by decreasing the number of adaption steps. The performance of the finite volume scheme is compared with several lower order methods.