๐”– Bobbio Scriptorium
โœฆ   LIBER   โœฆ

Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws

โœ Scribed by Michael Dumbser; Cedric Enaux; Eleuterio F. Toro


Publisher
Elsevier Science
Year
2008
Tongue
English
Weight
525 KB
Volume
227
Category
Article
ISSN
0021-9991

No coin nor oath required. For personal study only.

โœฆ Synopsis


In this article, we propose a new class of finite volume schemes of arbitrary accuracy in space and time for systems of hyperbolic balance laws with stiff source terms. The new class of schemes is based on a three stage procedure. First a highorder WENO reconstruction procedure is applied to the cell averages at the current time level. Second, the temporal evolution of the reconstruction polynomials is computed locally inside each cell using the governing equations. In the original ENO scheme of Harten et al. and in the ADER schemes of Titarev and Toro, this time evolution is achieved via a Taylor series expansion where the time derivatives are computed by repeated differentiation of the governing PDE with respect to space and time, i.e. by applying the so-called Cauchy-Kovalewski procedure. However, this approach is not able to handle stiff source terms. Therefore, we present a new strategy that only replaces the Cauchy-Kovalewski procedure compared to the previously mentioned schemes. For the time-evolution part of the algorithm, we introduce a local space-time discontinuous Galerkin (DG) finite element scheme that is able to handle also stiff source terms. This step is the only part of the algorithm which is locally implicit. The third and last step of the proposed ADER finite volume schemes consists of the standard explicit space-time integration over each control volume, using the local space-time DG solutions at the Gaussian integration points for the intercell fluxes and for the space-time integral over the source term. We will show numerical convergence studies for nonlinear systems in one space dimension with both non-stiff and with very stiff source terms up to sixth order of accuracy in space and time. The application of the new method to a large set of different test cases is shown, in particular the stiff scalar model problem of LeVeque and Yee [R.J. LeVeque, H.C. Yee, A study of numerical methods for hyperbolic conservation laws with stiff source terms, Journal of Computational Physics 86 (1) (1990) 187-210], the relaxation system of Jin and Xin [S. Jin, Z. Xin, The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications on Pure and Applied Mathematics 48 (1995) 235-277] and the full compressible Euler equations with stiff friction source terms.


๐Ÿ“œ SIMILAR VOLUMES


A Simple Riemann Solver and High-Order G
โœ Wenlong Dai; Paul R. Woodward ๐Ÿ“‚ Article ๐Ÿ“… 1995 ๐Ÿ› Elsevier Science ๐ŸŒ English โš– 517 KB

A simple approximate Riemann solver for hyperbolic systems of conservation laws is developed for its use in Godunov schemes. The solver is based on characteristic formulations and is illustrated through Euler and ideal magnetohydrodynamical (MHD) equations. The procedure of a high-order Godunov sche

Strict Stability of High-Order Compact I
โœ Saul S. Abarbanel; Alina E. Chertock ๐Ÿ“‚ Article ๐Ÿ“… 2000 ๐Ÿ› Elsevier Science ๐ŸŒ English โš– 248 KB

Temporal, or "strict," stability of approximation to PDEs is much more difficult to achieve than the "classical" Lax stability. In this paper, we present a class of finitedifference schemes for hyperbolic initial boundary value problems in one and two space dimensions that possess the property of st