paper

Unified discontinuous Galerkin scheme for a large class of elliptic equations

arXiv:2108.05826 · doi:10.1103/PhysRevD.105.024034

Abstract

We present a discontinuous Galerkin internal-penalty scheme that is applicable to a large class of linear and nonlinear elliptic partial differential equations. The unified scheme can accommodate all second-order elliptic equations that can be formulated in first-order flux form, encompassing problems in linear elasticity, general relativity, and hydrodynamics, including problems formulated on a curved manifold. It allows for a wide range of linear and nonlinear boundary conditions, and accommodates curved and nonconforming meshes. Our generalized internal-penalty numerical flux and our Schur-complement strategy of eliminating auxiliary degrees of freedom make the scheme compact without requiring equation-specific modifications. We demonstrate the accuracy of the scheme for a suite of numerical test problems. The scheme is implemented in the open-source SpECTRE numerical relativity code.

22 pages, 13 figures, 2 tables, published version. Results are reproducible with the ancillary input files

References in corpus (6)

Cited by in corpus (6)