Discontinuous Galerkin scheme for elliptic equations on extremely stretched grids
arXiv:2405.06120 · doi:10.1103/PhysRevD.110.084062
Abstract
Discontinuous Galerkin (DG) methods for solving elliptic equations are gaining popularity in the computational physics community for their high-order spectral convergence and their potential for parallelization on computing clusters. However, problems in numerical relativity with extremely stretched grids, such as initial data problems for binary black holes that impose boundary conditions at large distances from the black holes, have proven challenging for DG methods. To alleviate this problem we have developed a primal DG scheme that is generically applicable to a large class of elliptic equations, including problems on curved and extremely stretched grids. The DG scheme accommodates two widely used initial data formulations in numerical relativity, namely the puncture formulation and the extended conformal thin-sandwich (XCTS) formulation. We find that our DG scheme is able to stretch the grid by a factor of and hence allows to impose boundary conditions at large distances. The scheme converges exponentially with resolution both for the smooth XCTS problem and for the nonsmooth puncture problem. With this method we are able to generate high-quality initial data for binary black hole problems using a parallelizable DG scheme. The code is publicly available in the open-source SpECTRE numerical relativity code.
arXiv comment14 pages, 10 figures. Results are reproducible with the ancillary input files
References in corpus (33)
- A simple construction of initial data for multiple black holes
- The SXS Collaboration catalog of binary black hole simulations
- A single-domain spectral method for black hole puncture data
- Conformal ``thin sandwich'' data for the initial-value problem of general relativity
- A multidomain spectral method for solving elliptic equations
- Initial Data for Numerical Relativity
- Extrinsic Curvature and the Einstein Constraints
- Binary-black-hole initial data with nearly-extremal spins
- Excision boundary conditions for black hole initial data
- Kadath: a spectral solver for theoretical physics
- Extra-Large Remnant Recoil Velocities and Spins from Near-Extremal-Bowen-York-Spin Black-Hole Binaries
- A new public code for initial data of unequal-mass, spinning compact-object binaries
- New conformally flat initial data for spinning black holes
- Key Elements of Robustness in Binary Black Hole Evolutions using Spectral Methods
- ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
- A new numerical method to construct binary neutron star initial data
- Massively Parallel Simulations of Binary Black Hole Intermediate-Mass-Ratio Inspirals
- GRAthena++: puncture evolutions on vertex-centered oct-tree AMR
- Improvements to the construction of binary black hole initial data
- Assessing the Readiness of Numerical Relativity for LISA and 3G Detectors
- Formulation of discontinuous Galerkin methods for relativistic astrophysics
- Black hole initial data on hyperboloidal slices
- Puncture Initial Data for Black-Hole Binaries with High Spins and High Boosts
- Comparison of binary black hole initial data sets
- The new discontinuous Galerkin methods based numerical relativity program Nmesh
- Initial data for perturbed Kerr black holes on hyperboloidal slices
- A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
- Spherically symmetric black hole spacetimes on hyperboloidal slices
- Unified discontinuous Galerkin scheme for a large class of elliptic equations
- A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
- Efficient simulations of high-spin black holes with a new gauge
- Extending superposed harmonic initial data to higher spin
- High-accuracy numerical models of Brownian thermal noise in thin mirror coatings