Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes
arXiv:1707.09910 · doi:10.1103/PhysRevD.97.084053
Abstract
We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of standard gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.
27 pages, 8 figures; v2 matches version published in PRD. Discussion improved and proof of strong hyperbolicity is now given for standard gauges
References in corpus (12)
- Extraction of Gravitational Waves in Numerical Relativity
- Covariant formulations of BSSN and the standard gauge
- On the final fate of compact boson star mergers
- Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
- An explicit harmonic code for black-hole evolution using excision
- High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
- Implementation of standard testbeds for numerical relativity
- Mathematical Issues in a Fully-Constrained Formulation of Einstein Equations
- Discontinuous Galerkin method for the spherically reduced BSSN system with second-order operators
- Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry
- Towards a gauge-polyvalent Numerical Relativity code
- Hyperbolicity of High Order Systems of Evolution Equations
Cited by in corpus (34)
- Black holes, gravitational waves and fundamental physics: a roadmap
- Neutron Star Merger Remnants
- High order ADER schemes for continuum mechanics
- ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
- Dynamical ejecta and nucleosynthetic yields from eccentric binary neutron-star mergers
- On GLM curl cleaning for a first order reduction of the CCZ4 formulation of the Einstein field equations
- High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension
- Shocks and quark-meson scatterings at large density
- General-relativistic neutron star evolutions with the discontinuous Galerkin method
- The new discontinuous Galerkin methods based numerical relativity program Nmesh
- A new continuum model for general relativistic viscous heat-conducting media
- Simulating magnetized neutron stars with discontinuous Galerkin methods
- Summation by Parts and Truncation Error Matching on Hyperboloidal Slices
- On numerical methods for hyperbolic PDE with curl involutions
- A spectral method algorithm for numerical simulations of gravitational fields
- Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell WENO finite-volume limiting for simulation of non-stationary compressible multicomponent reactive flows
- A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
- Well balanced finite volume schemes for shallow water equations on manifolds
- Extreme Gravitational Interactions in the Problem of Three Black Holes in General Relativity
- Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting for multidimensional detonation waves simulation
- Well-balanced high order finite difference WENO schemes for a first-order Z4 formulation of the Einstein field equations
- High order well-balanced Arbitrary-Lagrangian-Eulerian ADER discontinuous Galerkin schemes on general polygonal moving meshes
- First-order hyperbolic formulation of the pure tetrad teleparallel gravity theory
- Arbitrary high order ADER-DG method with local DG predictor for solutions of initial value problems for systems of first-order ordinary differential equations
- High-order discontinuous Galerkin schemes with subcell finite volume limiter and adaptive mesh refinement for a monolithic first-order BSSNOK formulation of the Einstein-Euler equations
- Apples with Apples comparison of 3+1 conformal numerical relativity schemes
- ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE
- High order ADER-DG method with local DG predictor for solutions of differential-algebraic systems of equations
- Compact Runge-Kutta flux reconstruction methods for non-conservative hyperbolic equations
- Neutron star evolution by combining discontinuous Galerkin and finite volume methods
- Aliasing Instabilities in the Numerical Evolution of the Einstein Field Equations
- Improving the local solution of the DG predictor of the ADER-DG method for solving systems of ordinary differential equations and its applicability to systems of differential-algebraic equations
- Detrimental task execution patterns in mainstream OpenMP runtimes
- Efficient implementation of ADER discontinuous Galerkin schemes for a scalable hyperbolic PDE engine