On GLM curl cleaning for a first order reduction of the CCZ4 formulation of the Einstein field equations

Abstract

In this paper we propose an extension of the generalized Lagrangian multiplier method (GLM) of Munz et al. (JCP 2000, JCP 2002), which was originally conceived for the numerical solution of the Maxwell and MHD equations with divergence-type involutions, to the case of hyperbolic PDE systems with curl-type involutions. The key idea here is to solve an augmented PDE system, in which curl errors propagate away via a Maxwell-type evolution system. The new approach is first presented on a simple model problem, in order to explain the basic ideas. Subsequently, we apply it to a strongly hyperbolic first order reduction of the CCZ4 formulation (FO-CCZ4) of the Einstein field equations of general relativity, which is endowed with 11 curl constraints. Several numerical examples, including the long-time evolution of a stable neutron star in anti-Cowling approximation, are presented in order to show the obtained improvements with respect to the standard formulation without special treatment of the curl involution constraints. The main advantages of the proposed GLM approach are its complete independence of the underlying numerical scheme and grid topology and its easy implementation into existing computer codes. However, this flexibility comes at the price of needing to add for each curl involution one additional 3 vector plus another scalar in the augmented system for homogeneous curl constraints, and even two additional scalars for non-homogeneous curl involutions. For the FO-CCZ4 system with 11 homogeneous curl involutions, this means that additional 44 evolution quantities need to be added.

References in corpus (33)

• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• On the Numerical Integration of Einstein's Field Equations
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• New Formalism for Numerical Relativity
• A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws
• Constraint damping in the Z4 formulation and harmonic gauge
• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Conformal and covariant formulation of the Z4 system with constraint-violation damping
• Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
• Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
• An efficient shock-capturing central-type scheme for multidimensional relativistic flows. I. Hydrodynamics
• ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
• A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
• General-covariant evolution formalism for Numerical Relativity
• High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids
• Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting
• High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
• Toward standard testbeds for numerical relativity
• Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement
• A symmetry-breaking mechanism for the Z4 general-covariant evolution system
• Continuum Mechanics and Thermodynamics in the Hamilton and the Godunov-type Formulations
• High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
• Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes
• Hyperbolicity of second-order in space systems of evolution equations
• A 3+1 covariant suite of Numerical Relativity Evolution Systems
• Numerical simulations with a first order BSSN formulation of Einstein's field equations
• Dynamical shift conditions for the Z4 and BSSN hyperbolic formalisms
• Continuum mechanics with torsion
• Towards a gauge-polyvalent Numerical Relativity code
• von Neumann Stability Analysis of Globally Constraint-Preserving DGTD and PNPM Schemes for the Maxwell Equations using Multidimensional Riemann Solvers
• Globally constraint-preserving FR/DG scheme for Maxwell's equations at all orders
• Efficient, Divergence-Free, High Order MHD on 3D Spherical Meshes with Optimal Geodesic Meshing
• Fifth order finite volume WENO in general orthogonally-curvilinear coordinates

Cited by in corpus (17)

• High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
• The Science of the Einstein Telescope
• High order ADER schemes for continuum mechanics
• A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
• High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension
• High order entropy preserving ADER scheme
• On numerical methods for hyperbolic PDE with curl involutions
• A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
• Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes
• High-order ADER Discontinuous Galerkin schemes for a symmetric hyperbolic model of compressible barotropic two-fluid flows
• Well-balanced high order finite difference WENO schemes for a first-order Z4 formulation of the Einstein field equations
• An exactly curl-free finite-volume scheme for a hyperbolic compressible barotropic two-phase model
• 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
• 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
• ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE
• Von Neumann Stability Analysis of DG-like and PNPM-like Schemes for PDEs that have Globally Curl-Preserving Evolution of Vector Fields