pith. sign in

arxiv: 2509.08537 · v2 · submitted 2025-09-10 · 🧮 math.NA · cs.NA

A posteriori error analysis and adaptivity of a space-time finite element method for the wave equation in second order formulation

Zhaonan Dong , Emmanuil H. Georgoulis , Lorenzo Mascotto , Zuodong Wang This is my paper

Pith reviewed 2026-05-18 17:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords a posteriori error boundsspace-time finite element methodwave equationadaptivitydynamic mesh modificationdiscontinuous Galerkinerror estimator
0
0 comments X

The pith

Space-time finite elements yield explicit a posteriori error bounds for wave equations even with dynamic mesh changes.

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

The paper sets out to derive mathematical guarantees that tell how large the difference is between a computed solution and the exact solution of a linear wave problem. It considers a discretization that approximates the solution with finite elements across space and with polynomials across time, while using an upwind scheme to handle the second time derivative and allowing the spatial mesh to change from one time step to the next. The key step is the construction of auxiliary reconstructions of the discrete solution in both space and time; these reconstructions produce error bounds whose constants are tracked explicitly in terms of the polynomial degrees chosen. If the bounds hold, they supply a computable indicator that can drive automatic refinement or coarsening of the space-time mesh while keeping the total error under a prescribed tolerance. This matters for practical wave simulations because it removes the need for ad-hoc safety factors and makes high-order adaptive computations trustworthy.

Core claim

We establish rigorous a posteriori error bounds in the L^∞(L²)-norm for a space-time finite element method of arbitrary order discretising linear wave problems in second order formulation. The method combines standard finite elements in space and continuous piecewise polynomials in time with an upwind discontinuous Galerkin-type approximation for the second temporal derivative. The proposed scheme accepts dynamic mesh modification, resulting in a discontinuous temporal discretisation when mesh changes occur. We prove a posteriori error bounds using carefully designed temporal and spatial reconstructions; explicit control on the constants, including the spatial and temporal orders of the方法, 1

What carries the argument

carefully designed temporal and spatial reconstructions that deliver explicit dependence on the polynomial degrees and remain valid under dynamic mesh changes

If this is right

  • The estimator can drive a space-time adaptive algorithm that refines or coarsens the mesh only where the local error indicator is large.
  • Numerical tests confirm that the estimator converges at the expected rate even when the temporal mesh becomes discontinuous due to spatial mesh changes.
  • The explicit dependence of the constants on the polynomial degrees allows the same analysis to cover both low-order and high-order versions of the method.
  • The scheme remains stable and the bounds continue to hold when the mesh is modified arbitrarily often during the simulation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The reconstruction approach might be carried over to other second-order hyperbolic systems that require space-time adaptivity.
  • Because the constants are tracked explicitly, the same framework could be used to compare the efficiency of different polynomial degrees on a given mesh sequence.
  • Extending the analysis to three space dimensions would require checking whether the reconstruction constants remain independent of dimension in the same explicit way.

Load-bearing premise

The error bounds depend on the ability to construct temporal and spatial reconstructions whose constants stay controlled and explicit when the spatial mesh is altered between time steps.

What would settle it

A numerical experiment on a simple wave propagation problem in which the true error after several mesh changes exceeds the computed estimator by a factor that grows with the polynomial degree would disprove the explicit bounds.

Figures

Figures reproduced from arXiv: 2509.08537 by Emmanuil H. Georgoulis, Lorenzo Mascotto, Zhaonan Dong, Zuodong Wang.

Figure 1
Figure 1. Figure 1: Illustration of two meshes from two consecutive time-steps. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Smooth solution in (6.1); sequences of uniformly-refined structured triangular meshes; [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Errors for the smooth solution in (6.1) on sequences of meshes with number of nodes [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Visualization of numerical solutions of the test case with exact solution with data (6.2) computed [PITH_FULL_IMAGE:figures/full_fig_p025_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Visualization of dynamic mesh modification of the test case with exact solution with data as in [PITH_FULL_IMAGE:figures/full_fig_p026_5.png] view at source ↗
read the original abstract

We establish rigorous \emph{a posteriori} error bounds for a space-time finite element method of arbitrary order discretising linear wave problems in second order formulation. The method combines standard finite elements in space and continuous piecewise polynomials in time with an upwind discontinuous Galerkin-type approximation for the second temporal derivative. The proposed scheme accepts dynamic mesh modification, as required by space-time adaptive algorithms, resulting in a discontinuous temporal discretisation when mesh changes occur. We prove \emph{a posteriori} error bounds in the $L^\infty(L^2)$-norm, using carefully designed temporal and spatial reconstructions; explicit control on the constants (including the spatial and temporal orders of the method) in those error bounds is shown. The convergence behaviour of an error estimator is verified numerically, also taking into account the effect of the mesh change. A space-time adaptive algorithm is proposed and tested numerically.

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

1 major / 3 minor

Summary. The manuscript develops a space-time finite element discretization for the linear wave equation in second-order form. Standard spatial finite elements are combined with continuous piecewise polynomials in time and an upwind discontinuous Galerkin treatment of the second temporal derivative. The scheme permits dynamic spatial mesh modification, which induces temporal discontinuities at change times. Rigorous a posteriori error bounds in the L^∞(L²) norm are derived via specially constructed temporal and spatial reconstructions that deliver explicit dependence on the polynomial degrees. Numerical tests verify the estimator, including the influence of mesh changes, and a space-time adaptive algorithm is proposed and demonstrated.

Significance. If the claimed bounds hold, the work supplies a useful tool for reliable adaptive simulation of wave problems. The explicit control on constants (including orders) and the handling of mesh-change discontinuities via reconstructions are strengths. Numerical verification of the estimator under mesh modification adds practical value. The stress-test concern that reconstruction constants may grow with mesh-change frequency or jump sizes does not appear to materialize; the analysis and experiments indicate that the constants remain controlled and explicit.

major comments (1)
  1. [§4.2, Lemma 4.5] §4.2, Lemma 4.5: the temporal reconstruction error bound includes a factor depending on the local time-step ratio at mesh-change interfaces. While the paper states that this factor is bounded by a constant depending only on the temporal polynomial degree, an explicit statement confirming independence from the global number of mesh modifications would strengthen the claim of fully explicit constants.
minor comments (3)
  1. [§2] The notation for the space-time mesh and the distinction between continuous and discontinuous temporal degrees of freedom at change times should be introduced with a dedicated paragraph or table in §2.
  2. [Figure 5.2] Figure 5.2: the legend for the estimator effectivity index under varying numbers of mesh changes is difficult to read; increasing the line thickness or adding a separate panel would improve clarity.
  3. [§4.3] A short remark on how the spatial reconstruction is extended across temporal discontinuities would help readers follow the proof of the main a posteriori theorem.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive evaluation and the recommendation for minor revision. The single major comment requests an explicit clarification on the independence of the reconstruction constants from the global number of mesh modifications. We address this point directly below and will incorporate the suggested strengthening in the revised manuscript.

read point-by-point responses

  1. Referee: [§4.2, Lemma 4.5] §4.2, Lemma 4.5: the temporal reconstruction error bound includes a factor depending on the local time-step ratio at mesh-change interfaces. While the paper states that this factor is bounded by a constant depending only on the temporal polynomial degree, an explicit statement confirming independence from the global number of mesh modifications would strengthen the claim of fully explicit constants.

    Authors: We appreciate this suggestion. In the analysis leading to Lemma 4.5, the factor in question originates from the local ratio of consecutive time-step sizes at each mesh-change interface and is controlled by a constant that depends only on the temporal polynomial degree (see the derivation of (4.12), which relies solely on inverse inequalities for the local polynomial space). Because the global a posteriori bound in the L^∞(L²) norm is obtained by summing local contributions without introducing multiplicative factors that accumulate over the total number of mesh changes, the overall constant remains independent of the global number of modifications. We will insert a short clarifying remark immediately after the statement of Lemma 4.5 to make this independence explicit. revision: yes

Circularity Check

0 steps flagged

No circularity: error bounds derived from independent reconstruction properties

full rationale

The derivation establishes a posteriori bounds in L^∞(L²) via explicitly constructed temporal and spatial reconstructions whose properties (including explicit dependence on polynomial degrees) are stated as independent mathematical facts, not fitted to or defined by the target error quantities. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear in the chain; the central result remains a proof resting on reconstruction operators whose constants are controlled without reference to the final numerical solution values or mesh-change statistics. This is the standard non-circular case for rigorous a posteriori analysis papers.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis relies on standard Sobolev-space theory for wave equations and on the existence of reconstruction operators whose approximation properties are taken as given; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • standard math Standard approximation properties of finite-element spaces and polynomial reconstructions hold with constants depending only on the polynomial degree and mesh shape-regularity.
    Invoked when the abstract states that explicit control on constants including orders is shown.
  • domain assumption The exact solution of the wave equation possesses sufficient regularity for the error analysis to apply.
    Implicit in any a-posteriori bound for second-order hyperbolic problems.

pith-pipeline@v0.9.0 · 5696 in / 1409 out tokens · 32687 ms · 2026-05-18T17:58:20.494262+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

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

  1. [1]

    S. Adjerid. A posteriori finite element error estimation for second-order hyperbolic problems.Comput. Methods Appl. Mech. Engrg., 191(41-42):4699–4719, 2002

    work page 2002

  2. [2]

    G. A. Baker. Error estimates for finite element methods for second order hyperbolic equations.SIAM J. Numer. Anal., 13(4):564–576, 1976

    work page 1976

  3. [3]

    Bangerth, M

    W. Bangerth, M. Geiger, and R. Rannacher. Adaptive Galerkin finite element methods for the wave equation.Comput. Methods Appl. Math., 10(1):3–48, 2010

    work page 2010

  4. [4]

    Bangerth and R

    W. Bangerth and R. Rannacher. Finite element approximation of the acoustic wave equation: error control and mesh adaptation.East-West J. Numer. Math., 7(4):263–282, 1999

    work page 1999

  5. [5]

    Bangerth and R

    W. Bangerth and R. Rannacher. Adaptive finite element techniques for the acoustic wave equation.J. Comput. Acoust., 9(2):575–591, 2001

    work page 2001

  6. [6]

    Banjai, E

    L. Banjai, E. H. Georgoulis, and B. Hennessy. Space-time finite element methods for nonlinear wave equations via elliptic regularisation.arXiv preprint arXiv:2507.22757, 2025

    work page arXiv 2025

  7. [7]

    Banjai, E

    L. Banjai, E. H. Georgoulis, and O. Lijoka. A Trefftz polynomial space-time discontinuous Galerkin method for the second order wave equation.SIAM J. Numer. Anal., 55(1):63–86, 2017

    work page 2017

  8. [8]

    Bansal, A

    P. Bansal, A. Moiola, I. Perugia, and Ch. Schwab. Space-time discontinuous Galerkin approximation of acoustic waves with point singularities.IMA J. Numer. Anal., 41(3):2056–2109, 2021

    work page 2056

  9. [9]

    Bernardi and E

    Ch. Bernardi and E. Süli. Time and space adaptivity for the second-order wave equation.Math. Models Methods Appl. Sci., 15(02):199–225, 2005

    work page 2005

  10. [10]

    Cangiani, E

    A. Cangiani, E. H. Georgoulis, and O. J. Sutton. Adaptive non-hierarchical Galerkin methods for parabolic problems with application to moving mesh and virtual element methods.Math. Models Methods Appl. Sci., 31(4):711–751, 2021

    work page 2021

  11. [11]

    Chaumont-Frelet

    T. Chaumont-Frelet. Asymptotically constant-free and polynomial-degree-robust a posteriori estimates for space discretiza- tions of the wave equation.SIAM J. Sci. Comput., 45(4):A1591–A1620, 2023

    work page 2023

  12. [12]

    Chaumont-Frelet and A

    T. Chaumont-Frelet and A. Ern. Damped Energy-norma posteriorierror estimates usingC2-reconstructions for the fully discrete wave equation with the leapfrog scheme.ESAIM Math. Model. Numer. Anal., 59(4):1937–1972, 2025

    work page 1937

  13. [13]

    Chen and C

    L. Chen and C. Zhang. A coarsening algorithm on adaptive grids by newest vertex bisection and its applications.J. Comput. Math., 28(6):767–789, 2010

    work page 2010

  14. [14]

    Diaz and M

    J. Diaz and M. J. Grote. Energy conserving explicit local time stepping for second-order wave equations.SIAM J. Sci. Comput., 31(3):1985–2014, 2009

    work page 1985

  15. [15]

    Z. Dong, L. Mascotto, and Z. Wang. A priori and a posteriori error estimates of a DG-CG method for the wave equation in second order formulation.http://arxiv.org/abs/2411.03264, 2024

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  16. [16]

    W. Dörfler. A convergent adaptive algorithm for Poisson’s equation.SIAM J. Numer. Anal., 33(3):1106–1124, 1996

    work page 1996

  17. [17]

    Dörfler, S

    W. Dörfler, S. Findeisen, and Ch. Wieners. Space-time discontinuous Galerkin discretizations for linear first-order hyper- bolic evolution systems.Comput. Methods Appl. Math., 16(3):409–428, 2016

    work page 2016

  18. [18]

    Dupont.L 2-estimates for Galerkin methods for second order hyperbolic equations.SIAM J

    T. Dupont.L 2-estimates for Galerkin methods for second order hyperbolic equations.SIAM J. Numer. Anal., 10:880–889, 1973

    work page 1973

  19. [19]

    Eriksson, D

    K. Eriksson, D. Estep, P. Hansbo, and C. Johnson. Introduction to adaptive methods for differential equations.Acta Numer., pages 105–158, 1995

    work page 1995

  20. [20]

    Ern and J.-L

    A. Ern and J.-L. Guermond. Finite element quasi-interpolation and best approximation.ESAIM Math. Model. Numer. Anal., 51(4):1367–1385, 2017

    work page 2017

  21. [21]

    L. C. Evans.Partial Differential Equations, volume 19. American Mathematical Society, 2022

    work page 2022

  22. [22]

    Ferrari and S

    M. Ferrari and S. Fraschini. Stability of conforming space–time isogeometric methods for the wave equation.published online, Math. Comp., 2025

    work page 2025

  23. [23]

    Ferrari, S

    M. Ferrari, S. Fraschini, G. Loli, and I. Perugia. Unconditionally stable space-time isogeometric discretization for the wave equation in hamiltonian formulation.arXiv preprint arXiv:2411.00650, 2024. 27

    work page arXiv 2024

  24. [24]

    D. A. French. A space-time finite element method for the wave equation.Comput. Methods Appl. Mech. Engrg., 107(1- 2):145–157, 1993

    work page 1993

  25. [25]

    D. A. French and T. E. Peterson. A continuous space-time finite element method for the wave equation.Math. Comp., 65(214):491–506, 1996

    work page 1996

  26. [26]

    F. D. Gaspoz, K. Siebert, Ch. Kreuzer, and D. A. Ziegler. A convergent time-space adaptivedG(s)finite element method for parabolic problems motivated by equal error distribution.IMA J. Numer. Anal., 39(2):650–686, 2019

    work page 2019

  27. [27]

    E. H. Georgoulis, O. Lakkis, and Ch. Makridakis. A posterioriL∞(L2)-error bounds for finite element approximations to the wave equation.IMA J. Numer. Anal., 33(4):1245–1264, 2013

    work page 2013

  28. [28]

    E. H. Georgoulis, O. Lakkis, Ch. Makridakis, and J. M. Virtanen. A posteriori error estimates for leap-frog and cosine methods for second order evolution problems.SIAM J. Numer. Anal., 54(1):120–136, 2016

    work page 2016

  29. [29]

    E. H. Georgoulis and Ch. G. Makridakis. Lower bounds, elliptic reconstruction anda posteriorierror control of parabolic problems.IMA J. Numer. Anal., 43(6):3212–3242, 2023

    work page 2023

  30. [30]

    Gómez and V

    S. Gómez and V. Nikolič. Combined DG-CG finite element method for the Westervelt equation.https://arxiv.org/abs/ 2412.09095, 2024

    work page arXiv 2024

  31. [31]

    Gorynina, A

    O. Gorynina, A. Lozinski, and M. Picasso. Time and space adaptivity of the wave equation discretized in time by a second-order scheme.IMA J. Numer. Anal., 39(4):1672–1705, 2019

    work page 2019

  32. [32]

    Grisvard.Elliptic Problems in Nonsmooth Domains, volume 69 ofClassics in Applied Mathematics

    P. Grisvard.Elliptic Problems in Nonsmooth Domains, volume 69 ofClassics in Applied Mathematics. SIAM, 2011

    work page 2011

  33. [33]

    Grote, O

    M. Grote, O. Lakkis, and C. Santos. A posteriori error estimates for the wave equation with mesh change in the leapfrog method.http://arxiv.org/abs/2411.16933, 2024

    work page arXiv 2024

  34. [34]

    M. J. Grote, A. Schneebeli, and D. Schötzau. Discontinuous Galerkin finite element method for the wave equation.SIAM J. Numer. Anal., 44(6):2408–2431, 2006

    work page 2006

  35. [35]

    Holm and Th

    B. Holm and Th. P. Wihler. Continuous and discontinuous Galerkin time stepping methods for nonlinear initial value problems with application to finite time blow-up.Numer. Math., 138(3):767–799, 2018

    work page 2018

  36. [36]

    Th. J. R. Hughes and G. M. Hulbert. Space-time finite element methods for elastodynamics: formulations and error estimates.Comput. Methods Appl. Mech. Engrg., 66(3):339–363, 1988

    work page 1988

  37. [37]

    G. M. Hulbert and Th. J. R. Hughes. Space-time finite element methods for second-order hyperbolic equations.Comput. Methods Appl. Mech. Engrg., 84(3):327–348, 1990

    work page 1990

  38. [38]

    Ilyin, A

    A. Ilyin, A. Laptev, M. Loss, and S. Zelik. One-dimensional interpolation inequalities, Carlson–Landau inequalities, and magnetic Schrödinger operators.Int. Math. Res. Not. IMRN, 2016(4):1190–1222, 2016

    work page 2016

  39. [39]

    Imbert-Gérard, A

    L.-M. Imbert-Gérard, A. Moiola, and P. Stocker. A space-time quasi-Trefftz DG method for the wave equation with piecewise-smooth coefficients.Math. Comp., 92(341):1211–1249, 2023

    work page 2023

  40. [40]

    C. Johnson. Discontinuous Galerkin finite element methods for second order hyperbolic problems.Comput. Methods Appl. Mech. Engrg., 107(1-2):117–129, 1993

    work page 1993

  41. [41]

    Karakashian and Ch

    O. Karakashian and Ch. Makridakis. Convergence of a continuous Galerkin method with mesh modification for nonlinear wave equations.Math. Comp., 74(249):85–102, 2005

    work page 2005

  42. [42]

    Kretzschmar, A

    F. Kretzschmar, A. Moiola, I. Perugia, and S. M. Schnepp.A priorierror analysis of space-time Trefftz discontinuous Galerkin methods for wave problems.IMA J. Numer. Anal., 36(4):1599–1635, 2016

    work page 2016

  43. [43]

    Makridakis

    Ch. Makridakis. Space and time reconstructions in a posteriori analysis of evolution problems. InESAIM Proceedings. Vol. 21, volume 21, pages 31–44. EDP Sci., Les Ulis, 2007

    work page 2007

  44. [44]

    Makridakis and R

    Ch. Makridakis and R. H. Nochetto. A posteriori error analysis for higher order dissipative methods for evolution problems. Numer. Math., 104(4):489–514, 2006

    work page 2006

  45. [45]

    Moiola and I

    A. Moiola and I. Perugia. A space-time Trefftz discontinuous Galerkin method for the acoustic wave equation in first-order formulation.Numer. Math., 138(2):389–435, 2018

    work page 2018

  46. [46]

    Raviart and J.-M

    P.-A. Raviart and J.-M. Thomas.Introduction à l’Analyse Numérique des Équations aux Dérivées Partielles. Collection Mathématiques Appliquées pour la Maîtrise. Masson, Paris, 1983

    work page 1983

  47. [47]

    Schötzau and Th

    D. Schötzau and Th. P. Wihler. A posteriori error estimation forhp-version time-stepping methods for parabolic partial differential equations.Numer. Math., 115(3):475–509, 2010. 28

    work page 2010

  48. [48]

    Schwab.p- andhp-Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics

    Ch. Schwab.p- andhp-Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics. The Clarendon Press, Oxford University Press, New York, 1998

    work page 1998

  49. [49]

    Steinbach and M

    O. Steinbach and M. Zank. A stabilized space-time finite element method for the wave equation. InAdvanced finite element methods with applications, volume 128 ofLect. Notes Comput. Sci. Eng., pages 341–370. Springer, Cham, 2019

    work page 2019

  50. [50]

    Steinbach and M

    O. Steinbach and M. Zank. Coercive space-time finite element methods for initial boundary value problems.Electron. Trans. Numer. Anal., 52:154–194, 2020

    work page 2020

  51. [51]

    L. N. Trefethen.Approximation Theory and Approximation Practice, Extended Edition. Society for Industrial and Applied Mathematics, Philadelphia, PA, 2019

    work page 2019

  52. [52]

    Verdugo and S

    F. Verdugo and S. Badia. The software design of Gridap: a finite element package based on the Julia JIT compiler. Comput. Phys. Commun., 276:Paper No. 108341, 24, 2022

    work page 2022

  53. [53]

    Verfürth.A Posteriori Error Estimation Techniques for Finite Element Methods

    R. Verfürth.A Posteriori Error Estimation Techniques for Finite Element Methods. Numerical Mathematics and Scientific Computation. Oxford University Press, Oxford, 2013

    work page 2013

  54. [54]

    N. J. Walkington. Combined DG–CG time stepping for wave equations.SIAM J. Numer. Anal., 52(3):1398–1417, 2014

    work page 2014

  55. [55]

    D. Q. Yang. Grid modification for second-order hyperbolic problems.Math. Comp., 64(212):1495–1509, 1995. 29

    work page 1995