Theory and internal structure of ADER-DG method for ordinary differential equations
Pith reviewed 2026-05-18 22:43 UTC · model grok-4.3
The pith
The ADER-DG method for ODE systems is A-, AN-, L-, B-, BN-, and algebraically stable.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The ADER-DG method is A- and AN-stable, L-stable, B- and BN-stable, and algebraically stable. Several other relations useful for an application and implementation of the ADER-DG method are proved. Applications of the ADER-DG method demonstrated compliance with the expected theoretical results.
What carries the argument
The ADER-DG time-stepping scheme whose stability is established through direct analysis of its amplification matrix and Butcher tableau properties for general ODE right-hand sides.
If this is right
- The method can be used on stiff ODEs without artificial step-size limits from stability concerns.
- Algebraic stability guarantees bounded long-term behavior for nonlinear problems.
- Convergence follows from the established stability properties for sufficiently smooth solutions.
- Implementation can rely on the derived relations for efficient coefficient computation.
- Applications match theory, confirming the method works as predicted in practice.
Where Pith is reading between the lines
- Similar stability might extend to the method when applied to semi-discretized PDEs.
- Direct comparison with high-order Runge-Kutta schemes could quantify accuracy gains per computational cost.
- The algebraic stability property may imply preservation of quadratic invariants in certain conservative systems.
- Testing on specific stiff models such as chemical reaction networks would provide practical validation.
Load-bearing premise
The stability proofs hold for arbitrary ODE systems without hidden restrictions on the right-hand side that would invalidate the algebraic stability arguments.
What would settle it
Numerical integration of a stiff linear test equation with large time steps that shows the method violating L-stability or algebraic stability bounds.
read the original abstract
Investigation of the approximation properties, convergence, and stability of the ADER-DG method for solving an ODE system is carried out. The ADER-DG method is $A$- and $AN$-stable, $L$-stable, $B$- and $BN$-stable, and algebraically stable. Several other relations useful for an application and implementation of the ADER-DG method are proved. Applications of the ADER-DG method demonstrated compliance with the expected theoretical results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the approximation properties, convergence, and stability of the ADER-DG method for systems of ordinary differential equations. It claims that the method is A- and AN-stable, L-stable, B- and BN-stable, and algebraically stable, proves several auxiliary relations for implementation, and reports that numerical applications comply with the theoretical predictions.
Significance. If the stability claims hold for general nonlinear ODEs without unstated restrictions, the work would be a useful contribution to numerical ODE methods by providing a DG-based scheme with algebraic stability, a property that guarantees favorable behavior for stiff problems and is usually associated with certain implicit Runge-Kutta methods. The internal-structure analysis could also facilitate implementation and further extensions.
major comments (2)
- [§4, Theorem 4.3] §4 (Stability Properties), Theorem 4.3 on algebraic stability: the derivation of the energy estimate leading to algebraic stability for arbitrary f does not state or verify the required bound on the one-sided Lipschitz constant of the right-hand side or the restriction on the time-step size h needed to preserve monotonicity; without this, the extension from linear to general nonlinear systems is not justified and directly affects the central stability claims.
- [§3.2] §3.2 (B- and BN-stability): the proof invokes an identity that appears to hold only when the local truncation error satisfies a specific dissipativity condition; no explicit assumption or counter-example check is provided for general f, undermining the claim that B/BN-stability holds unconditionally.
minor comments (2)
- [§2] The notation for the DG polynomial space and the choice of basis in §2 is introduced without a clear statement of the reference element or quadrature rule, which affects readability of later stability estimates.
- [Figure 1] Figure 1 caption does not specify the test problem or the exact parameter values used, making it hard to reproduce the compliance demonstration.
Simulated Author's Rebuttal
We thank the referee for the thorough review and insightful comments on the stability analysis. We address each major comment below and have revised the manuscript to incorporate clarifications on the assumptions required for the nonlinear case.
read point-by-point responses
-
Referee: [§4, Theorem 4.3] §4 (Stability Properties), Theorem 4.3 on algebraic stability: the derivation of the energy estimate leading to algebraic stability for arbitrary f does not state or verify the required bound on the one-sided Lipschitz constant of the right-hand side or the restriction on the time-step size h needed to preserve monotonicity; without this, the extension from linear to general nonlinear systems is not justified and directly affects the central stability claims.
Authors: We agree that the proof of algebraic stability in Theorem 4.3 requires explicit mention of the one-sided Lipschitz condition on f (with constant μ ≤ 0) and a corresponding restriction on the time step h to guarantee the energy estimate. These are standard assumptions for extending algebraic stability results to nonlinear problems. We have added a clear statement of these conditions at the start of Section 4 and inserted a short verification step in the proof of Theorem 4.3 to confirm monotonicity preservation under the stated bound on h. The core derivation remains unchanged, but the scope of the result is now precisely delineated. revision: yes
-
Referee: [§3.2] §3.2 (B- and BN-stability): the proof invokes an identity that appears to hold only when the local truncation error satisfies a specific dissipativity condition; no explicit assumption or counter-example check is provided for general f, undermining the claim that B/BN-stability holds unconditionally.
Authors: The identity invoked in the B- and BN-stability proof of §3.2 is valid when the local truncation error inherits the dissipativity property from the underlying ADER-DG discretization, which holds unconditionally for linear f and under the one-sided Lipschitz condition for nonlinear f. We have revised §3.2 to state this assumption explicitly at the outset of the subsection and added a remark clarifying that unconditional B/BN-stability is claimed only for the linear case, while the nonlinear extension requires the dissipativity condition on the truncation error. No counter-example is needed once the assumption is stated, as the proof then follows directly. revision: yes
Circularity Check
No significant circularity in ADER-DG stability derivation
full rationale
The paper derives stability properties (A/AN/L/B/BN/algebraic) for the ADER-DG scheme through direct mathematical analysis of its internal structure and approximation properties for general ODE systems. No steps reduce by construction to fitted parameters, self-definitions, or load-bearing self-citations; the claims rest on standard energy estimates and Butcher tableau conditions that are independently verifiable. The provided abstract and context show a self-contained proof chain without renaming known results or smuggling ansatzes via prior author work. This is the expected honest non-finding for a paper whose central results have external mathematical content.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The IRK-GL-ADER-DG method is algebraically stable. ... Matrix M is non-negative definite ... μ is positive definite
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.
Forward citations
Cited by 1 Pith paper
-
Improving the local solution of the DG predictor of the ADER-DG method for solving systems of ordinary differential equations and its applicability to systems of differential-algebraic equations
An improved local solution for the ADER-DG DG predictor converges one order higher than the original, is continuous at nodes, and applies to DAEs with rigorous proofs.
Reference graph
Works this paper leans on
-
[1]
I. Babuˇska and T. Strouboulis , The Finite Element Method and Its Reliability, Numerical Mathematics and Scientific Computation , Clarendon Press, Oxford, 2001
work page 2001
-
[2]
M. Baccouch, Analysis of a posteriori error estimates of the discontinuous Galerkin method for nonlinear ordinary differential equations, Appl. Numer. Math., 106 (2016), pp. 129–153
work page 2016
-
[3]
M. Baccouch , A posteriori error estimates and adaptivity for the discontinuous Galerkin solutions of nonlinear second-order initial-value problems, Appl. Numer. Math., 121 (2017), pp. 18–37
work page 2017
-
[4]
M. Baccouch, Superconvergence of the discontinuous Galerkin method for nonlinear second- order initial-value problems for ordinary differential equations , Appl. Numer. Math., 115 (2017), pp. 160–179. THEORY OF ADER-DG METHOD FOR ODES 19
work page 2017
-
[5]
M. Baccouch , Analysis of optimal superconvergence of the local discontinuous Galerkin method for nonlinear fourth-order boundary value problems , Numerical Algorithms, 86 (2021), pp. 1615–1650
work page 2021
-
[6]
M. Baccouch, The discontinuous Galerkin method for general nonlinear third-order ordinary differential equations, Appl. Numer. Math., 162 (2021), pp. 331–350
work page 2021
-
[7]
M. Baccouch, Superconvergence of an ultra-weak discontinuous Galerkin method for nonlin- ear second-order initial-value problems, International Journal of Computational Methods, 20(2) (2023), p. 2250042
work page 2023
- [8]
- [9]
-
[10]
Butcher, Numerical Methods for Ordinary Differential Equations , Wiley, UK, 2016
J. Butcher, Numerical Methods for Ordinary Differential Equations , Wiley, UK, 2016
work page 2016
-
[11]
B. Cockburn, S. Hou, , and C.-W. Shu , TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , Math. Comp., 54 (1990), pp. 545–581
work page 1990
-
[12]
B. Cockburn, S.-Y. Lin, and C.-W. Shu , TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One-dimensional systems , J. Comput. Phys., 84 (1989), pp. 90–113
work page 1989
-
[13]
B. Cockburn and C.-W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework , Math. Comp., 52 (1989), pp. 411–435
work page 1989
-
[14]
B. Cockburn and C.-W. Shu , The Runge-Kutta local projection P 1-discontinuous-Galerkin finite element method for scalar conservation laws, ESAIM: M2AN, 25 (1991), pp. 337–361
work page 1991
-
[15]
B. Cockburn and C.-W. Shu , TVB Runge-Kutta local projection discontinuous Galerkin Method for Conservation Laws V: Multidimensional Systems , J. Comput. Phys., 141 (1998), pp. 199–224
work page 1998
-
[16]
K. Dekker and J. Verwer, Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations, North-Holland, Amsterdam, 1984
work page 1984
-
[17]
M. Delfour and F. Dubeau, Discontinuous polynomial approximations in the theory of one- step, hybrid and multistep methods for nonlinear ordinary differential equations , Math. Comp., 47 (1986), p. 169
work page 1986
-
[18]
M. Delfour, W. Hager, and F. Trochu , Discontinuous Galerkin methods for ordinary dif- ferential equations, Math. Comp., 36 (1981), pp. 455–473
work page 1981
-
[19]
M. Dumbser, Arbitrary high order PN PM schemes on unstructured meshes for the compressible Navier-Stokes equations, Computers & Fluids, 39 (2010), pp. 60–76
work page 2010
-
[20]
M. Dumbser, C. Enaux, and E. Toro , Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , J. Comput. Phys., 227 (2008), p. 3971
work page 2008
-
[21]
M. Dumbser and O. Zanotti , Very high order PN PM schemes on unstructured meshes for the resistive relativistic mhd equations , J. Comput. Phys., 228 (2009), p. 6991
work page 2009
-
[22]
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. Comp. Phys., 504 (2024), p. 112875
work page 2024
-
[23]
E. Fernandez, M. Diaz, M. Dumbser, and T. de Luna, An arbitrary high order well-balanced ADER-DG numerical scheme for the multilayer shallow-water model with variable density , J. Sci. Comput., 90 (2022), p. 52
work page 2022
-
[24]
E. Gaburro and M. Dumbser , A posteriori subcell finite volume limiter for general PN PM schemes: Applications from gasdynamics to relativistic magnetohydrodynamics , J. Sci. Comput., 86 (2021), p. 37
work page 2021
-
[25]
E. Gaburro, P. Offner, M. Ricchiuto, and D. Torlo , High order entropy preserving ADER-DG schemes, Applied Mathematics and Computation, 440 (2023), p. 127644
work page 2023
- [26]
-
[27]
E. Hairer and G. Wanner, Solving Ordinary Differential Equations II: Stiff and Differential- Algebraic Problems, Springer-Verlag, Berlin, Heidelberg, 1996
work page 1996
-
[28]
M. Han Veiga, P. Offner, and D. Torlo , DeC and ADER: Similarities, differences and a unified framework, Journal of Scientific Computing, 87 (2021), p. 2
work page 2021
-
[29]
A. Hidalgo and M. Dumbser , ADER schemes for nonlinear systems of stiff advection- diffusion-reaction equations, J. Sci. Comput., 48 (2011), p. 173
work page 2011
-
[30]
A. Jorba and M. Zou , A software package for the numerical integration of ODEs by means of high-order Taylor methods , Experiment. Math., 14 (2005), pp. 99–117. 20 I. S. POPOV
work page 2005
-
[31]
S.-C. Klein, Stabilizing discontinuous Galerkin methods using Dafermos’ entropy rate crite- rion: I — one-dimensional conservation laws , J. Sci. Comput., 95 (2023), p. 55
work page 2023
-
[32]
L. Micalizzi, D. Torlo, and W. Boscheri , Efficient iterative arbitrary high-order methods: an adaptive bridge between low and high order , Commun. Appl. Math. Comput., 7 (2023), p. 40
work page 2023
-
[33]
I. Popov, Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell WENO finite-volume limiting for simulation of non-stationary com- pressible multicomponent reactive flows, J. Sci. Comput., 95 (2023), p. 44
work page 2023
-
[34]
I. Popov , Arbitrary high order ADER-DG method with local DG predictor for solutions of initial value problems for systems of first-order ordinary differential equations , J. Sci. Comput., 100 (2024), p. 22
work page 2024
-
[35]
I. Popov, Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting for multidimensional detonation waves simulation, Computers & Fluids, 284 (2024), p. 106425
work page 2024
-
[36]
I. Popov, High order ADER-DG method with local DG predictor for solutions of differential- algebraic systems of equations , J. Sci. Comput., 102 (2025), p. 48
work page 2025
-
[37]
W. Reed and T. Hill , Triangular mesh methods for the neutron transport equation , Tech. Report LA-UR-73-479, Los Alamos Scientific Laboratory, 1973
work page 1973
-
[38]
V. Titarev and E. Toro , ADER: arbitrary high order Godunov approach , J. Sci. Comput., 17 (2002), p. 609
work page 2002
-
[39]
V. Titarev and E. Toro, ADER schemes for three-dimensional nonlinear hyperbolic systems, J. Comput. Phys., 204 (2005), p. 715
work page 2005
-
[40]
E. Toro, V. Titarev, M. Dumbser, A. Iske, C. Goetz, C. Castro, G. Montecinos, and D. R., The ADER approach for approximating hyperbolic equations to very high accuracy , in Hyperbolic Problems: Theory, Numerics, Applications. Volume I. HYP 2022. SEMA SIMAI Springer Series, vol 34., C. Pares, M. Castro, T. de Luna, and M. Mu˜ noz-Ruiz, eds., Springer, Cham, 2024
work page 2022
-
[41]
L. Wahlbin, Superconvergence in Galerkin Finite Element Methods , Springer-Verlag, Berlin Heidelberg, 1995
work page 1995
-
[42]
S. Wolf, M. Galis, C. Uphoff, A.-A. Gabriel, M. P., D. Gregor, and M. Bader , An efficient ADER-DG local time stepping scheme for 3D HPC simulation of seismic waves in poroelastic media, J. Comp. Phys., 455 (2022), p. 110886
work page 2022
-
[43]
O. Zanotti, F. Fambri, M. Dumbser, and A. Hidalgo , Space-time adaptive ADER discon- tinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting , Computers & Fluids, 118 (2015), p. 204
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.