A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
arXiv:1907.01572 · doi:10.1103/PhysRevD.100.084052
Abstract
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.
arXiv comment20 pages, 21 figures
References in corpus (34)
- A simple construction of initial data for multiple black holes
- A single-domain spectral method for black hole puncture data
- Initial Data for Numerical Relativity
- A multidomain spectral method for solving elliptic equations
- Binary-black-hole initial data with nearly-extremal spins
- Excision boundary conditions for black hole initial data
- Three-dimensional initial data for the collision of two black holes II: Quasi-circular orbits for equal-mass black holes
- Circular orbits and spin in black-hole initial data
- SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
- Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement
- Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
- A performance comparison of continuous and discontinuous Galerkin methods with fast multigrid solvers
- ADER discontinuous Galerkin schemes for general-relativistic ideal magnetohydrodynamics
- Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
- A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
- Discontinuous Galerkin method for computing gravitational waveforms from extreme mass ratio binaries
- Solving 3D relativistic hydrodynamical problems with WENO discontinuous Galerkin methods
- The Einstein constraints: uniqueness and non-uniqueness in the conformal thin sandwich approach
- On thin plate spline interpolation
- Finite Element, Discontinuous Galerkin, and Finite Difference Evolution Schemes in Spacetime
- Discontinuous Galerkin method for the spherically reduced BSSN system with second-order operators
- Numerical simulations with a first order BSSN formulation of Einstein's field equations
- Formulation of discontinuous Galerkin methods for relativistic astrophysics
- Multi-Domain Spectral Method for Initial Data of Arbitrary Binaries in General Relativity
- Approximate initial data for binary black holes
- General-relativistic neutron star evolutions with the discontinuous Galerkin method
- An Operator-Based Local Discontinuous Galerkin Method Compatible With the BSSN Formulation of the Einstein Equations
- Initial-data contribution to the error budget of gravitational waves from neutron-star binaries
- Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates
- Implicit-explicit (IMEX) evolution of single black holes
- Robust Approaches to Handling Complex Geometries with Galerkin Difference Methods
- IMEX evolution of scalar fields on curved backgrounds
- An Energy Stable Approach for Discretizing Hyperbolic Equations with Nonconforming Discontinuous Galerkin Methods
- Initial data for high-compactness black hole-neutron star binaries
Cited by in corpus (5)
- Constructing Binary Neutron Star Initial Data with High Spins, High Compactness, and High Mass-Ratios
- A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
- thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state
- Unified discontinuous Galerkin scheme for a large class of elliptic equations
- Discontinuous Galerkin scheme for elliptic equations on extremely stretched grids