Formulation of discontinuous Galerkin methods for relativistic astrophysics

JOURNAL OF COMPUTATIONAL PHYSICS

卷:312

Article

Teukolsky, Saul A.^1,2 

[1] Cornell Univ, Dept Phys & Astron, Space Sci Bldg, Ithaca, NY 14853 USA

[2] CALTECH, TAPIR, Walter Burke Inst Theoret Phys, MC 350-17, Pasadena, CA 91125 USA

关键词: Discontinuous Galerkin;    Hydrodynamics;    Magnetohydrodynamics;    Einstein's equations;    Moving mesh;    Arbitrary Lagrangian-Eulerian (ALE);    Metric identities;    Geometric conservation law;   

DOI  :  10.1016/j.jcp.2016.02.031

来源: Elsevier

【 摘 要 】

The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to arbitrary time-dependent curved spacetimes. We show how to formulate the algorithm in such spacetimes for applications in relativistic astrophysics. We also show how to formulate the algorithm for equations in non-conservative form, such as Einstein's field equations themselves. We find two computationally distinct formulations in both cases, one of which has seldom been used before for flat space in curvilinear coordinates but which may be more efficient. We also give a new derivation of the ALE algorithm (Arbitrary Lagrangian-Eulerian) using 4-vector methods that is much simpler than the usual derivation and explains why the method preserves the conservation form of the equations. The various formulations are explored with some simple numerical experiments that also investigate the effect of the metric identities on the results. The results of this paper may also be of interest to practitioners of DG working with curvilinear elements in flat space. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free