A Convergent Hybridizable Discontinuous Galerkin Method for Einstein--Scalar Equations
This work provides a novel numerical method for simulating gravitational collapse in astrophysics, though it is incremental as it adapts an existing technique to a specific system.
The authors tackled the numerical solution of the spherically symmetric Einstein--scalar equations by developing a hybridizable discontinuous Galerkin method, achieving an optimal order L^2 error bound for the main evolution variable and verifying spatial convergence rates in numerical experiments.
We propose and analyze a hybridized discontinuous Galerkin (HDG) method for the spherically symmetric Einstein--scalar system in Bondi gauge. After rewriting the model as a local first-order PDE--ODE system by introducing suitable scaled variables, we construct a semidiscrete scheme in which the element unknowns are computed locally and the coupling is carried by traces on the mesh skeleton. In the present radial setting, these traces can be eliminated recursively, so that only the main evolution variable is advanced in time, while the metric variables are recovered from discrete constraint relations. We prove local semidiscrete well-posedness, derive a global \(L^2\)--stability estimate, establish an optimal order \(L^2\) error bound for the main evolution variable for polynomial degree \(k\ge 1\), and obtain reconstruction error estimates for the metric variables and the associated mass functional. Numerical experiments verify the predicted spatial convergence rate and illustrate qualitative features of the Einstein--scalar dynamics, including large-data collapse profiles and smooth-pulse evolution.