An operator-based local discontinuous Galerkin method compatible with the BSSN formulation of the Einstein equations. (8th December 2016)

• Record Type: Journal Article
• Authors: Miller, Jonah M; Schnetter, Erik
• Is Part Of: Classical and quantum gravity. Volume 34:Number 1(2017:Jan.)
• Publication Date: 2016-12-08
• Subjects: numerical relativity -- einstein toolkit -- discontinuous Galerkin finite elements -- BSSN; Quantum gravity -- Periodicals; Gravitation -- Periodicals; Relativity (Physics) -- Periodicals; Space and time -- Periodicals; Periodicals; 521.1
• DOI: 10.1088/1361-6382/34/1/015003
• Languages: English
• ISSNs: 0264-9381

Abstract

Abstract: Discontinuous Galerkin finite element (DGFE) methods offer a mathematically beautiful, computationally efficient, and efficiently parallelizable way to solve partial differential equations (PDEs). These properties make them highly desirable for numerical calculations in relativistic astrophysics and many other fields. The BSSN formulation of the Einstein equations has repeatedly demonstrated its robustness. The formulation is not only stable but allows for puncture-type evolutions of black hole systems. To-date no one has been able to solve the full (3 + 1)-dimensional BSSN equations using DGFE methods. This is partly because DGFE discretization often occurs at the level of the equations, not the derivative operator, and partly because DGFE methods are traditionally formulated for manifestly flux-conservative systems. By discretizing the derivative operator, we generalize a particular flavor of DGFE methods, Local DG methods, to solve arbitrary second-order hyperbolic equations. Because we discretize at the level of the derivative operator, our method can be interpreted as either a DGFE method or as a finite differences stencil with non-constant coefficients.