DeepAI AI Chat
Log In Sign Up

A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity

by   Trevor Vincent, et al.

A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other discontinuities. This paper investigates discontinuous Galerkin methods for the solution of elliptic problems in numerical relativity. We present a novel hp-adaptive numerical scheme for curvilinear and non-conforming meshes. It uses a multigrid preconditioner with a Chebyshev or Schwarz smoother to create a very scalable discontinuous Galerkin code on generic domains. The code employs compactification to move the outer boundary near spatial infinity. We explore the properties of the code on some test problems, including one mimicking Neutron stars with phase transitions. We also apply it to construct initial data for two or three black holes.


page 12

page 14

page 15

page 17

page 18


A discontinuous Galerkin scheme for a large class of elliptic equations

We present a discontinuous Galerkin internal-penalty scheme that is appl...

A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework

In this paper, we present a discontinuous Galerkin solver based on previ...

Robust Multigrid Methods for Discontinuous Galerkin Discretizations of an Elliptic Optimal Control Problem

We consider discontinuous Galerkin methods for an elliptic distributed o...

A DPG method for shallow shells

We develop and analyze a discontinuous Petrov–Galerkin method with optim...

Yet Another Tensor Toolbox for discontinuous Galerkin methods and other applications

The numerical solution of partial differential equations is at the heart...