pith. sign in

arxiv: 2512.24121 · v2 · submitted 2025-12-30 · 🧮 math.NA · cs.NA· gr-qc

High order numerical discretizations of the Einstein-Euler equations in the Generalized Harmonic formulation

Pith reviewed 2026-05-16 19:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NAgr-qc
keywords Einstein-Euler equationsGeneralized Harmonic formulationwell-balancing propertyCWENO schemeADER discontinuous Galerkinnumerical relativityunstructured meshesblack hole accretion
0
0 comments X

The pith

Two high-order schemes with well-balancing exactly preserve stationary solutions for Einstein-Euler equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops two new high-order numerical schemes for the coupled Einstein-Euler equations in the Generalized Harmonic formulation. A finite difference central weighted essentially non-oscillatory scheme operates on Cartesian meshes, while an ADER discontinuous Galerkin scheme works on two-dimensional unstructured polygonal meshes. Both include a well-balancing property to keep known stationary equilibria exact at the discrete level. Validation covers vacuum tests, stable black hole evolutions including extreme Kerr, and matter interactions like accretion and neutron star perturbations. These techniques lay groundwork for complex three-dimensional simulations of gravitational systems with fluids.

Core claim

We propose two new alternative numerical schemes to solve the coupled Einstein-Euler equations in the Generalized Harmonic formulation. The first one is a finite difference central weighted essentially non-oscillatory scheme on a traditional Cartesian mesh, while the second one is an ADER discontinuous Galerkin scheme on 2D unstructured polygonal meshes, both equipped with a well-balancing property that preserves the equilibrium of a priori known stationary solutions exactly at the discrete level.

What carries the argument

The well-balancing property integrated into the FD CWENO and ADER DG discretizations, which ensures exact discrete preservation of stationary solutions for the Einstein-Euler system.

Load-bearing premise

The demonstrated well-balancing and stability in selected vacuum and matter tests will extend to the full nonlinear three-dimensional system on unstructured moving meshes without new instabilities.

What would settle it

Observing the development of instabilities or failure to preserve equilibrium in a three-dimensional simulation of a binary merger or similar dynamical system using the proposed schemes would indicate the assumption does not hold.

Figures

Figures reproduced from arXiv: 2512.24121 by Elena Gaburro, Michael Dumbser, Olindo Zanotti, Stefano Muzzolon.

Figure 1
Figure 1. Figure 1: An example of a rectangular 2D unstructured mesh covered by non-overlapping polygons [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the constraints for the robust stability test case with a random initial perturbation of amplitude [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Solution of the linearized wave test using the fourth order ADER-DG scheme on an unstructured polygonal mesh. Left [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Solution of the gauge wave test with amplitude [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Solution of the gauge wave test with high amplitude [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: 2D Schwarzschild black hole subject to an initial Gaussian perturbation in the variable [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: 2D Schwarzschild black hole in 3D Cartesian harmonic coordinates evolved over the equatorial plane [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: 2D Kerr black holes subject to an initial Gaussian perturbation in the variable [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: 3D Schwarzschild black hole in Cartesian Kerr-Schild coordinates using the seventh order CWENO-FD scheme. Left panel: [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: 3D Kerr black hole in Cartesian Kerr-Schild coordinates using the seventh order CWENO-FD scheme. Top panels: spin [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Spherical accretion onto Schwarzschild black hole in 2D spheroidal Kerr-Schild coordinates using the fifth order CWENO [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Time evolution of the normalized central density for the perturbed non-rotating neutron star solved with the fifth order [PITH_FULL_IMAGE:figures/full_fig_p021_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Time evolution of the constraints for the perturbed non-rotating neutron star. Comparison between the well-balanced (WB) [PITH_FULL_IMAGE:figures/full_fig_p022_13.png] view at source ↗
read the original abstract

We propose two new alternative numerical schemes to solve the coupled Einstein-Euler equations in the Generalized Harmonic formulation. The first one is a finite difference (FD) Central Weighted Essentially Non-Oscillatory (CWENO) scheme on a traditional Cartesian mesh, while the second one is an ADER (Arbitrary high order Derivatives) discontinuous Galerkin (DG) scheme on 2D unstructured polygonal meshes. The latter, in particular, represents a preliminary step in view of a full 3D numerical relativity calculation on moving meshes. Both schemes are equipped with a well-balancing (WB) property, which allows to preserve the equilibrium of a priori known stationary solutions exactly at the discrete level. We validate our numerical approaches by successfully reproducing standard vacuum test cases, such as the robust stability, the linearized wave, and the gauge wave tests, as well as achieving long-term stable evolutions of stationary black holes, including Kerr black holes with extreme spin. Concerning the coupling with matter, modeled by the relativistic Euler equations, we perform some special relativistic Riemann problems, a classical test of spherical accretion onto a Schwarzschild black hole, as well as an evolution of a perturbed non-rotating neutron star, demonstrating the capability of our schemes to operate also on the full Einstein-Euler system. Altogether, these results provide a solid foundation for addressing more complex and challenging simulations of astrophysical sources through DG schemes on unstructured 3D meshes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper proposes two new high-order numerical schemes for the coupled Einstein-Euler equations in the Generalized Harmonic formulation: a finite-difference Central Weighted Essentially Non-Oscillatory (CWENO) scheme on Cartesian meshes and an ADER discontinuous Galerkin (DG) scheme on 2D unstructured polygonal meshes. Both schemes incorporate a well-balancing property that exactly preserves a priori known stationary equilibria at the discrete level. Validation includes vacuum tests (robust stability, linearized waves, gauge waves), long-term evolutions of stationary black holes including extreme-spin Kerr, and matter-coupled tests (special-relativistic Riemann problems, spherical accretion onto Schwarzschild, and a perturbed non-rotating neutron star). The 2D DG scheme is presented as a preliminary step toward future 3D moving-mesh calculations.

Significance. If the reported results hold, the work supplies concrete, high-order, well-balanced discretizations that advance numerical relativity toward unstructured-mesh methods capable of handling the full nonlinear Einstein-Euler system. The successful long-term stability on extreme Kerr and selected matter tests, together with the explicit well-balancing construction, provides a reproducible foundation that can be directly compared against existing codes and extended to more complex astrophysical configurations.

minor comments (2)
  1. [Abstract] Abstract: the phrasing 'solid foundation for addressing more complex... simulations... through DG schemes on unstructured 3D meshes' should be qualified to reflect that all reported results are either vacuum or 2D; a single sentence clarifying the scope of the current validation would prevent overstatement.
  2. [Introduction] The manuscript would benefit from an explicit statement (perhaps in §1 or the conclusions) of which test problems are performed in 3D Cartesian FD-CWENO versus 2D polygonal ADER-DG, so that readers can immediately map the reported stability results to the two distinct discretizations.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The report correctly summarizes the contributions of the two high-order well-balanced schemes for the Einstein-Euler system in generalized harmonic gauge. Below we provide point-by-point responses to the major comments.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The manuscript constructs two original high-order discretizations (Cartesian FD-CWENO and 2D polygonal ADER-DG) for the Einstein-Euler system in generalized harmonic form and equips them with a well-balancing property by design. Validation proceeds via direct comparison to independent, externally known solutions (robust stability, linearized waves, gauge waves, stationary black holes including extreme Kerr, special-relativistic Riemann problems, spherical accretion, perturbed neutron stars). No equation or claim reduces by construction to a fitted parameter, self-referential definition, or self-citation chain; the central results are new numerical schemes whose correctness is assessed against benchmarks outside the paper's own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard numerical analysis assumptions for hyperbolic systems and the well-posedness of the generalized harmonic formulation; no new free parameters, ad-hoc entities, or non-standard axioms are introduced in the abstract.

axioms (2)
  • standard math The generalized harmonic formulation yields a strongly hyperbolic system suitable for high-order numerical discretization.
    Invoked implicitly when applying FD and DG schemes to the Einstein-Euler system.
  • domain assumption Well-balancing can be achieved by exact discrete preservation of known stationary solutions without affecting the scheme's order of accuracy.
    Central design choice stated in the abstract for both proposed methods.

pith-pipeline@v0.9.0 · 5569 in / 1577 out tokens · 52013 ms · 2026-05-16T19:25:56.685599+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

107 extracted references · 107 canonical work pages · 1 internal anchor

  1. [1]

    Alcubierre.Introduction to 3+1 numerical relativity, volume 140

    M. Alcubierre.Introduction to 3+1 numerical relativity, volume 140. Oxford University Press, 2008. 22

  2. [2]

    Alcubierre, G

    M. Alcubierre, G. Allen, C. Bona, D. Fiske, T. Goodale, F. S. Guzmán, I. Hawke, S. H. Haw- ley, S. Husa, M. Koppitz, C. Lechner, D. Pollney, D. Rideout, M. Salgado, E. Schnetter, E. Seidel, H. Shinkai, D. Shoemaker, B. Szilágyi, R. Takahashi, and J. Winicour. Towards standard testbeds for numerical relativity.Classical and Quantum Gravity, 21(2):589–613, ...

  3. [3]

    D. Alic, C. Bona-Casas, C. Bona, L. Rezzolla, and C. Palenzuela. Conformal and covariant for- mulation of the Z4 system with constraint-violation damping.Physical Review D, 85(6):064040, 2012

  4. [4]

    D. Alic, W. Kastaun, and L. Rezzolla. Constraint damping of the conformal and covariant formu- lation of the Z4 system in simulations of binary neutron stars.Physical Review D, 88(6):064049, 2013

  5. [5]

    K. Alvi. First-order symmetrizable hyperbolic formulations of Einstein’s equations including lapse and shift as dynamical fields.Classical and Quantum Gravity, 19(20):5153–5162, October 2002

  6. [6]

    Antón, O

    L. Antón, O. Zanotti, J. A. Miralles, J. M. Martí, J. M. Ibáñez, J. A. Font, and J. A. Pons. Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach.The Astrophysical Journal, 637(1):296, 2006

  7. [7]

    Audusse, F

    E. Audusse, F. Bouchut, M. O. Bristeau, R. Klein, and B. Perthame. A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows.SIAM Journal on Scientific Com- puting, 25(6):2050–2065, 2004

  8. [8]

    M. C. Babiuc, B. Szilágyi, and J. Winicour. Harmonic initial-boundary evolution in general relativity. Phys. Rev. D, 73:064017, Mar 2006

  9. [9]

    D. S. Balsara, D. Bhoriya, C.-W. Shu, and H. Kumar. Efficient Finite Difference WENO Scheme for Hyperbolic Systems with Non-Conservative Products.arXiv e-prints, page arXiv:2303.17672, March 2023

  10. [10]

    D. S. Balsara, D. Bhoriya, C.-W. Shu, and H. Kumar. Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Systems with Non-Conservative Products.arXiv e-prints, page arXiv:2403.01266, March 2024

  11. [11]

    D. S. Balsara, D. Bhoriya, O. Zanotti, and M. Dumbser. Well-balanced High-order Finite Difference Weighted Essentially Nonoscillatory Schemes for a First-order Z4 Formulation of the Einstein Field Equations.The Astrophysical Journal Supplement Series, 275(1):18, November 2024

  12. [12]

    D. S. Balsara and M. Dumbser. Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers.Journal of Computational Physics, pages 687–715, 2015

  13. [13]

    D. S. Balsara, M. Dumbser, and R. Abgrall. Multidimensional HLLC Riemann Solver for Unstruc- tured Meshes - With Application to Euler and MHD Flows.Journal of Computational Physics, 261:172–208, 2014

  14. [14]

    Banyuls, J

    F. Banyuls, J. A. Font, J. M. Ibáñez, J. M. Martí, and J. A. Miralles. Numerical 3+1 general relativistic hydrodynamics: A local characteristic approach.The Astrophysical Journal, 476(1):221, 1997

  15. [15]

    T. W. Baumgarte and S. L. Shapiro. Numerical integration of Einstein’s field equations.Physical Review D, 59(2):024007, 1998

  16. [16]

    T. W. Baumgarte and S. L. Shapiro.Numerical relativity: solving Einstein’s equations on the com- puter. Cambridge University Press, 2010. 23

  17. [17]

    J. P. Berberich, P. Chandrashekar, and C. Klingenberg. High order well-balanced finite volume meth- ods for multi-dimensional systems of hyperbolic balance laws.Computers and Fluids, 219:104858, 2021

  18. [18]

    Bermudez and M

    A. Bermudez and M. E. Vázquez-Cendón. Upwind methods for hyperbolic conservation laws with source terms.Computers&Fluids, 23(8):1049–1071, 1994

  19. [19]

    Birke, W

    C. Birke, W. Boscheri, and C. Klingenberg. A well-balanced semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach numbers.Journal of Scientific Computing, 98(2):34, 2024

  20. [20]

    C. Bona, T. Ledvinka, C. Palenzuela, and M. Zácek. Symmetry-breaking mechanism for the Z4 general-covariant evolution system.Phys. Rev. D, 69(6):064036, March 2004

  21. [21]

    Bona and C

    C. Bona and C. Palenzuela-Luque.Elements of Numerical Relativity. Springer-Verlag, Berlin, 2005

  22. [22]

    H. Bondi. On spherically symmetrical accretion.Monthly Notices of the Royal Astronomical Society, 112(2):195–204, 1952

  23. [23]

    Boscheri and M

    W. Boscheri and M. Dumbser. Arbitrary–Lagrangian–Eulerian One–Step WENO Finite V olume Schemes on Unstructured Triangular Meshes.Communications in Computational Physics, 14:1174– 1206, 2013

  24. [24]

    Boscheri and M

    W. Boscheri and M. Dumbser. A direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume scheme on unstructured tetrahedral meshes for conservative and non-conservative hyperbolic systems in 3D.Journal of Computational Physics, 275:484 – 523, 2014

  25. [25]

    Boscheri and M

    W. Boscheri and M. Dumbser. Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes.Journal of Computational Physics, 346:449 – 479, 2017

  26. [26]

    Boscheri, M

    W. Boscheri, M. Dumbser, and O. Zanotti. High Order Cell-Centered Lagrangian-Type Finite V ol- ume Schemes with Time-Accurate Local Time Stepping on Unstructured Triangular Meshes.Journal of Computational Physics, 291:120–150, 2015

  27. [27]

    Well balanced finite volume methods for nearly hydrostatic flows.Journal of Computational Physics, 196(2):539–565, 2004

    N Botta, R Klein, S Langenberg, and S Lützenkirchen. Well balanced finite volume methods for nearly hydrostatic flows.Journal of Computational Physics, 196(2):539–565, 2004

  28. [28]

    Boyle, L

    M. Boyle, L. Lindblom, H. P. Pfeiffer, M. A. Scheel, and L. E. Kidder. Testing the accuracy and stability of spectral methods in numerical relativity.Phys. Rev. D, 75:024006, Jan 2007

  29. [29]

    J. D. Brown, P. Diener, S. E. Field, J. S. Hesthaven, F. Herrmann, A. H. Mroué, O. Sarbach, E. Schnet- ter, M. Tiglio, and M. Wagman. Numerical simulations with a first-order BSSN formulation of Einstein’s field equations.Phys. Rev. D, 85(8):084004, April 2012

  30. [30]

    Bugner.Discontinuous Galerkin methods for general relativistic hydrodynamics

    M. Bugner.Discontinuous Galerkin methods for general relativistic hydrodynamics. PhD thesis, Friedrich-Schiller-Universität Jena, 2018

  31. [31]

    Busto, S

    S. Busto, S. Chiocchetti, M. Dumbser, E. Gaburro, and I. Peshkov. High order ADER schemes for continuum mechanics.Frontiers in Physiccs, 8:32, 2020

  32. [32]

    M. J. Castro, J. Gallardo, and C. Parés. High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems.Mathematics of computation, 75(255):1103–1134, 2006. 24

  33. [33]

    M. J. Castro, J. M. Gallardo, J. A. López, and C. Parés. Well-balanced high order extensions of Godunov’s method for semilinear balance laws.SIAM Journal on Numerical Analysis, 46(2):1012– 1039, 2008

  34. [34]

    M. J. Castro and C. Parés. Well-balanced high-order finite volume methods for systems of balance laws.Journal of Scientific Computing, 82(2):1–48, 2020

  35. [35]

    G. B. Cook and M. A. Scheel. Well-behaved harmonic time slices of a charged, rotating, boosted black hole.Physical Review D, 56(8):4775, 1997

  36. [36]

    Cordero-Carrión, J

    I. Cordero-Carrión, J. M. Ibáñez, E. Gourgoulhon, J. L. Jaramillo, and J. Novak. Mathematical issues in a fully constrained formulation of the Einstein equations.Phys. Rev. D, 77(8):084007, April 2008

  37. [37]

    T. G. Cowling. The non-radial oscillations of polytropic stars.Mon. Not. R. Astron. Soc., 101:367, 1941

  38. [38]

    Cravero, G

    I. Cravero, G. Puppo, M. Semplice, and G. Visconti. CWENO: uniformly accurate reconstructions for balance laws.Mathematics of Computation, 87(312):1689–1719, 2018

  39. [39]

    Cravero and M

    I. Cravero and M. Semplice. On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes.Journal of Scientific Computing, 67(3):1219–1246, 2016

  40. [40]

    De Donder

    T. De Donder. The mathematical theory of relativity. 1927

  41. [41]

    Del Zanna, O

    L. Del Zanna, O. Zanotti, N. Bucciantini, and P. Londrillo. ECHO: a Eulerian conservative high- order scheme for general relativistic magnetohydrodynamics and magnetodynamics.Astronomy& Astrophysics, 473(1):11–30, 2007

  42. [42]

    Deppe, F

    N. Deppe, F. Foucart, M. S. Bonilla, M. Boyle, N. J. Corso, M. D. Duez, M. Giesler, F. Hébert, L. E. Kidder, Y . Kim, P. Kumar, I. Legred, G. Lovelace, E. R. Most, J. Moxon, K. C. Nelli, H. P. Pfeiffer, M. A. Scheel, S. A. Teukolsky, W. Throwe, and N. L. Vu. Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method.Classic...

  43. [43]

    Deppe, W

    N. Deppe, W. Throwe, L. E. Kidder, N. L. Vu, K. C. Nelli, C. Armaza, M. S. Bonilla, F. Hébert, Y . Kim, P. Kumar, G. Lovelace, A. Macedo, J. Moxon, E. O’Shea, H. P. Pfeiffer, M. A. Scheel, S. A. Teukolsky, N. A. Wittek, I. Anantpurkar, C. Anderson, M. Boyle, A. Carpenter, A. Ceja, H. Chaudhary, N. Corso, L. N. Fayyazuddin, F. Foucart, N. Ghadiri, M. Giesl...

  44. [44]

    M. Dumbser. Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier–Stokes equations.Computers&Fluids, 39:60–76, 2010

  45. [45]

    Dumbser, D

    M. Dumbser, D. S. Balsara, E. F. Toro, and C.-D. Munz. A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes.Journal of Computational Physics, 227(18):8209–8253, 2008

  46. [46]

    Dumbser, W

    M. Dumbser, W. Boscheri, M. Semplice, and G. Russo. Central weighted ENO schemes for hy- perbolic conservation laws on fixed and moving unstructured meshes.SIAM Journal on Scientific Computing, 39(6):A2564–A2591, 2017. 25

  47. [47]

    Dumbser, M

    M. Dumbser, M. J. Castro, C. Parés, and E. F. Toro. ADER schemes on unstructured meshes for Non-Conservative hyperbolic systems: Applications to geophysical flows.Computers and Fluids, 38:1731––1748, 2009

  48. [48]

    Dumbser, C

    M. Dumbser, C. Enaux, and E. F. Toro. Finite volume schemes of very high order of accuracy for stiffhyperbolic balance laws.Journal of Computational Physics, 227:3971–4001, 2008

  49. [49]

    Dumbser, F

    M. Dumbser, F. Fambri, E. Gaburro, and A. Reinarz. On GLM curl cleaning for a first order re- duction of the CCZ4 formulation of the Einstein field equations.Journal of Computational Physics, 404:109088, 2020

  50. [50]

    Dumbser, F

    M. Dumbser, F. Guercilena, S. Köppel, L. Rezzolla, and O. Zanotti. Conformal and covariant Z4 formulation of the Einstein equations: Strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes.Physical Review D, 97(8):084053, 2018

  51. [51]

    Dumbser, A

    M. Dumbser, A. Hidalgo, M. J. Castro, C. Parés, and E. F. Toro. FORCE schemes on unstructured meshes II: Non-conservative hyperbolic systems.Computer Methods in Applied Mechanics and Engineering, 199(9):625–647, 2010

  52. [52]

    Dumbser and C.-D

    M. Dumbser and C.-D. Munz. Building blocks for arbitrary high order discontinuous Galerkin schemes.Journal of Scientific Computing, 27:215–230, 2006

  53. [53]

    Dumbser and O

    M. Dumbser and O. Zanotti. Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations.Journal of Computational Physics, 228(18):6991–7006, 2009

  54. [54]

    Dumbser, O

    M. Dumbser, O. Zanotti, E. Gaburro, and I. Peshkov. A well-balanced discontinuous Galerkin method for the first–order Z4 formulation of the Einstein–Euler system.J. Comput. Phys., 504:112875, 2024

  55. [55]

    Dumbser, O

    M. Dumbser, O. Zanotti, and I. Peshkov. 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.Phys. Rev. D, 110(8):084015, October 2024

  56. [56]

    Dumbser, O

    M. Dumbser, O. Zanotti, and G. Puppo. Monolithic first-order BSSNOK formulation of the Einstein- Euler equations and its solution with path-conservative finite difference central WENO schemes. Phys. Rev. D, 111(10):104072, May 2025

  57. [57]

    Fambri, M

    F. Fambri, M. Dumbser, S. Köppel, L. Rezzolla, and O. Zanotti. ADER discontinuous Galerkin schemes for general-relativistic ideal magnetohydrodynamics.Monthly Notices of the Royal Astro- nomical Society, 477(4):4543–4564, 2018

  58. [58]

    J. A. Font, T. Goodale, S. Iyer, M. Miller, L. Rezzolla, E. Seidel, N. Stergioulas, W.-M. Suen, and M. Tobias. Three-dimensional numerical general relativistic hydrodynamics. II. Long-term dynamics of single relativistic stars.Physical Review D, 65(8):084024, April 2002

  59. [59]

    E. Gaburro. A unified framework for the solution of hyperbolic PDE systems using high order di- rect Arbitrary–Lagrangian–Eulerian schemes on moving unstructured meshes with topology change. Archives of Computational Methods in Engineering, 28(3):1249–1321, 2021

  60. [60]

    E. Gaburro. High order well-balanced Arbitrary-Lagrangian-Eulerian ADER discontinuous Galerkin schemes on general polygonal moving meshes.Computers&Fluids, 301:106764, 2025

  61. [61]

    Gaburro, W

    E. Gaburro, W. Boscheri, S. Chiocchetti, C. Klingenberg, V . Springel, and M. Dumbser. High order direct Arbitrary-Lagrangian-Eulerian schemes on moving V oronoi meshes with topology changes. Journal of Computational Physics, 407:109167, 2020. 26

  62. [62]

    Gaburro, M

    E. Gaburro, M. J. Castro, and M. Dumbser. A Well Balanced Finite V olume Scheme for General Relativity.SIAM Journal on Scientific Computing, 43(6):B1226–B1251, January 2021

  63. [63]

    Gaburro and S

    E. Gaburro and S. Chiocchetti. High-order Arbitrary-Lagrangian-Eulerian schemes on crazy moving V oronoi meshes. InYoung Researchers Conference, pages 99–119. Springer, 2021

  64. [64]

    Gaburro, M

    E. Gaburro, M. Dumbser, and M. J. Castro. Direct Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming unstructured meshes.Computers&Fluids, 159:254–275, 2017

  65. [65]

    Gassner, M

    G. Gassner, M. Dumbser, F. Hindenlang, and C.-D. Munz. Explicit One–Step time discretizations for discontinuous Galerkin and finite volume schemes based on local predictors.Journal of Compu- tational Physics, 230:4232–4247, 2011

  66. [66]

    Ghosh and E

    D. Ghosh and E. M. Constantinescu. Well-Balanced, Conservative Finite Difference Algorithm for Atmospheric Flows.AIAA Journal, 54(4):1370–1385, 2016

  67. [67]

    L. Gosse. A well-balanced scheme using non-conservative products designed for hyperbolic systems of conservation laws with source terms.Mathematical Models and Methods in Applied Sciences, 11(02):339–365, 2001

  68. [68]

    Gourgoulhon.3+1 Formalism in General Relativity, volume 846

    E. Gourgoulhon.3+1 Formalism in General Relativity, volume 846. 2012

  69. [69]

    Hilditch, S

    D. Hilditch, S. Bernuzzi, M. Thierfelder, Z. Cao, W. Tichy, and B. Brügmann. Compact binary evolutions with the Z4c formulation.Phys. Rev. D, 88:084057, Oct 2013

  70. [70]

    G. S. Jiang and C.-W. Shu. Efficient implementation of weighted ENO schemes.Journal of Compu- tational Physics, 126:202–228, 1996

  71. [71]

    R. P. Kerr. Gravitational field of a spinning mass as an example of algebraically special metrics. Phys. Rev. Lett., 11:237–238, 1963

  72. [72]

    L. E. Kidder, S. E. Field, F. Foucart, E. Schnetter, S. A. Teukolsky, A. Bohn, N. Deppe, P. Diener, F. Hébert, J. Lippuner, J. Miller, C. D. Ott, M. A. Scheel, and T. Vincent. SpECTRE: A task- based discontinuous Galerkin code for relativistic astrophysics.Journal of Computational Physics, 335:84–114, April 2017

  73. [73]

    S. S. Komissarov. Electrodynamics of black hole magnetospheres.Monthly Notices of the Royal Astronomical Society, 350(2):427–448, 2004

  74. [74]

    R. J. LeVeque. Balancing source terms and flux gradients in high-resolution Godunov methods: the quasi-steady wave-propagation algorithm.Journal of Computational Physics, 146(1):346–365, 1998

  75. [75]

    D. Levy, G. Puppo, and G. Russo. Central WENO schemes for hyperbolic systems of conservation laws.Mathematical Modelling and Numerical Analysis, 33(3):547 – 571, 1999

  76. [76]

    D. Levy, G. Puppo, and G. Russo. Third order central WENO scheme for 2D conservation laws. Applied Numerical Mathematics, 33(1):415 – 421, 2000

  77. [77]

    D. Levy, G. Puppo, and G. Russo. Compact central WENO schemes for multidimensional conserva- tion laws.SIAM Journal on Scientific Computing, 22(2):656 – 672, 2001

  78. [78]

    Lindblom, K

    L. Lindblom, K. D. Matthews, O. Rinne, and M. A. Scheel. Gauge drivers for the generalized harmonic Einstein equations.Physical Review D—Particles, Fields, Gravitation, and Cosmology, 77(8):084001, 2008. 27

  79. [79]

    Lindblom, M

    L. Lindblom, M. A. Scheel, L. E. Kidder, R. Owen, and O. Rinne. A new generalized harmonic evolution system.Classical and Quantum Gravity, 23(16):S447–S462, August 2006

  80. [80]

    Lindblom and B

    L. Lindblom and B. Szilágyi. Improved gauge driver for the generalized harmonic Einstein system. Phys. Rev. D, 80:084019, Oct 2009

Showing first 80 references.