A Convergent Hybridizable Discontinuous Galerkin Method for Einstein--Scalar Equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Dwivedi, Mukul, Rupp, Andreas
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911569438310400
author Dwivedi, Mukul
Rupp, Andreas
author_facet Dwivedi, Mukul
Rupp, Andreas
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2604_04613
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle A Convergent Hybridizable Discontinuous Galerkin Method for Einstein--Scalar Equations
Dwivedi, Mukul
Rupp, Andreas
Numerical Analysis
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.
title A Convergent Hybridizable Discontinuous Galerkin Method for Einstein--Scalar Equations
topic Numerical Analysis
url https://arxiv.org/abs/2604.04613