Pith. sign in

REVIEW 2 minor 34 references

Enriched higher-order LOD achieves optimal convergence rates for the multiscale wave equation

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-06-29 05:54 UTC pith:J2JFVO6F

load-bearing objection This extends enriched LOD corrections from parabolic to wave problems and claims to remove the order-2 saturation via a priori estimates plus numerics.

arxiv 2605.30118 v1 pith:J2JFVO6F submitted 2026-05-28 math.NA cs.NA

Enriched higher-order multiscale approaches with applications to wave propagation

Balaje Kalyanaraman , Felix Krumbiegel , Roland Maier , Siyang Wang This is my paper
classification math.NA cs.NA
keywords localized orthogonal decompositionmultiscale methodswave equationhigh-order convergenceenriched correctionsa priori error estimates
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper extends higher-order localized orthogonal decomposition methods from elliptic problems to linear hyperbolic ones, specifically the wave equation with strongly heterogeneous coefficients. Earlier versions of the method were limited to second-order spatial convergence regardless of polynomial degree, but the new enriched corrections remove this saturation. Under minimal assumptions on the coefficient and standard well-prepared data conditions, the authors prove a priori error estimates that recover the full high-order rates. The corrections decay exponentially and are localized to patches, allowing practical computation. Numerical tests using a fifth-order Rosenbrock-Wanner time integrator confirm the predicted spatial and temporal rates for varying polynomial degrees.

Core claim

The enriched higher-order localized orthogonal decomposition method for the wave equation produces a priori error estimates with optimal high-order convergence rates in space. The enriched corrections exhibit exponential decay, are computed on patches, and overcome the second-order saturation seen in prior constructions, provided the data satisfy standard well-preparedness conditions and the coefficient meets only minimal assumptions.

What carries the argument

The enriched corrections that are added to the higher-order LOD basis functions; they carry the additional information needed to restore optimal polynomial convergence rates while retaining exponential decay for localization.

Load-bearing premise

The initial data and right-hand side must satisfy the standard well-preparedness conditions that the error analysis invokes.

What would settle it

A concrete computation on a heterogeneous wave problem with data that deliberately violate well-preparedness, checking whether the observed spatial convergence rate remains capped at two or drops below the predicted high-order rate.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Optimal high-order spatial rates hold for the wave equation once the enriched corrections are included.
  • Localization to patches remains valid because the corrections still decay exponentially.
  • The combination with a fifth-order Rosenbrock-Wanner integrator preserves optimal temporal accuracy in the numerical tests.
  • The same enrichment strategy that worked for parabolic problems transfers directly to the hyperbolic setting.

Where Pith is reading between the lines

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

  • The patch-wise construction suggests the method can be parallelized efficiently for large three-dimensional domains.
  • Similar enrichment ideas could be tested on other time-dependent multiscale equations such as Maxwell or elasticity systems.
  • If the well-preparedness assumption can be relaxed or replaced by weaker conditions, the applicability of the method would widen considerably.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 2 minor

Summary. The manuscript provides an overview of higher-order localized orthogonal decomposition (LOD) methods for elliptic PDEs with heterogeneous coefficients, then generalizes the approach to linear hyperbolic multiscale problems (the wave equation). It proposes an enriched higher-order LOD construction that incorporates enriched corrections with exponential decay (computable on patches), derives a priori error estimates achieving optimal high-order spatial convergence rates under minimal coefficient assumptions plus standard well-preparedness conditions on the data, and verifies the claims numerically using a fifth-order Rosenbrock-Wanner time integrator, including localization-error studies for varying polynomial orders.

Significance. If the a priori estimates hold, the work is significant for multiscale wave propagation: it removes the second-order saturation barrier of prior LOD constructions for hyperbolic problems while retaining the localization and patch-computability advantages of the method. The explicit use of enriched corrections (building on the parabolic case) and the numerical verification with a high-order ROW integrator are strengths that support practical use in applications such as acoustics or seismics.

minor comments (2)
  1. [Abstract] Abstract: the phrase 'including recent advancements' for elliptic LOD is vague; adding one or two specific citations or a one-sentence characterization of the key prior results would improve readability without lengthening the abstract.
  2. The well-preparedness conditions on the data are invoked for the high-order rates but are described only as 'standard'; a brief reminder of the precise conditions (e.g., compatibility with the initial data or source term) in the statement of the main theorem would help readers assess applicability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. The referee's description of the manuscript is accurate. No major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The derivation builds a priori error estimates for the enriched higher-order LOD method on established LOD theory and standard well-preparedness conditions on the data, without any reduction of the central claims to self-citations, fitted inputs renamed as predictions, or self-definitional steps. The abstract and overview explicitly position the work as a generalization of prior LOD constructions with a new enrichment step for hyperbolic problems, and the error estimates are derived under explicit minimal assumptions rather than by construction from the inputs. No load-bearing uniqueness theorems or ansatzes are imported via self-citation in a way that collapses the result.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on domain assumptions about coefficients and data preparation that are standard in multiscale analysis but not independently verified here; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption minimal assumptions on the coefficient
    Invoked to derive the a priori error estimates for high-order convergence
  • domain assumption standard well-preparedness conditions on the data
    Required for the error estimates to achieve optimal rates

pith-pipeline@v0.9.1-grok · 5741 in / 1295 out tokens · 30951 ms · 2026-06-29T05:54:56.299719+00:00 · methodology

0 comments
read the original abstract

We consider the numerical solution of partial differential equations with coefficients that are strongly heterogeneous in space. We provide an overview of higher-order localized orthogonal decomposition (LOD) methods for the elliptic setting, including recent advancements, and then present a generalization of the strategy to linear hyperbolic multiscale problems. We address the limitations of earlier constructions for the wave equation, which only achieve second-order convergence in space, independent of the chosen polynomial degree. Building on the methodology of enriched corrections recently developed for parabolic multiscale problems, we motivate and propose an enriched higher-order LOD method for the wave equation. The enriched corrections exhibit exponential decay and can be computed on patches. Under minimal assumptions on the coefficient and standard well-preparedness conditions on the data, we derive a priori error estimates that achieve optimal high-order convergence rates, thereby overcoming the previously observed saturation of the convergence rate. With the fifth-order Rosenbrock-Wanner (ROW) time integrator, we conduct a series of numerical examples to verify our theoretical results. We provide examples showing the optimal spatial convergence of the method including the localization errors for different polynomial orders. We also present examples showing the optimal convergence rates of the time discretization.

Figures

Figures reproduced from arXiv: 2605.30118 by Balaje Kalyanaraman, Felix Krumbiegel, Roland Maier, Siyang Wang.

Figure 2.1
Figure 2.1. Figure 2.1: Illustration of a zero-order bubble function (left) and its sta￾bilized (extended) version (right) for p = 0, and H = 1 5 in one dimension. continuous and fulfill zero boundary conditions. More precisely, for any v ∈ L 2 (Ω) at an interior node z ∈ TH, we have with the definition of patches in (2.12) that (2.17) (IHv)(z) := X K∈N({z}) |K| |N({z})| Z K v dx, see [EG17] for further details on such construc… view at source ↗
Figure 3.1
Figure 3.1. Figure 3.1: Illustration of two localization strategies for iterative applications of elliptic operators: a classical localization strategy (left) and the more involved localization used in [KKMW25] that yields much smaller patches (right). can be combined. Note that λ can be decreased with an increasing distance (in terms of layers of elements) between G1 and K. This argument can be repeated also to correct the con… view at source ↗
Figure 4.1
Figure 4.1. Figure 4.1: Energy errors for Example 4.1 with different localization pa￾rameters for p = 1 (left) and p = 3 (right). The stabilization strategy proposed in [HLM26] is used in both cases. correction level j = ⌈p/2⌉ (see Theorem 3.1) for selected choices of the localization pa￾rameter ℓ. The results are shown in [PITH_FULL_IMAGE:figures/full_fig_p016_4_1.png] view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Energy errors for the enriched multiscale method applied to Example 4.2 with p = 1, j = 1 (left) and p = 3, j = 2 (right), different localization parameters, and different localizations strategies. performs significantly worse, partially by almost two orders of magnitude. For complete￾ness, [PITH_FULL_IMAGE:figures/full_fig_p017_4_2.png] view at source ↗
Figure 4.3
Figure 4.3. Figure 4.3: Optimal choice of the localization parameter ℓ for p = 2, 3, 4 and optimal j in Example 4.2 (left) and the energy errors for Exam￾ples 4.1 and 4.2 showing the temporal convergence of the Rodas5P scheme (right). with the theoretical results, confirming both the accuracy and robustness of the proposed approach. Acknowledgments F. Krumbiegel and R. Maier acknowledge funding from the Deutsche Forschungsgemei… view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

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

  1. [1]

    Abdulle and M

    A. Abdulle and M. Grote. Finite element heterogeneous multiscale method for the wave equation. Multiscale Model. Simul. , 9:766--792, 2011

    2011

  2. [2]

    Abdulle and P

    A. Abdulle and P. Henning. Localized orthogonal decomposition method for the wave equation with a continuum of scales. Math. Comp. , 86:549--587, 2017

    2017

  3. [3]

    Altmann, P

    R. Altmann, P. Henning, and D. Peterseim. Numerical homogenization beyond scale separation. Acta Numer. , 30:1--86, 2021

    2021

  4. [4]

    Badia and F

    S. Badia and F. Verdugo. Gridap: An extensible finite element toolbox in J ulia. J. Open Source Software , 5(52):2520, 2020

    2020

  5. [5]

    Z. Dong, M. Hauck, and R. Maier. An improved high-order method for elliptic multiscale problems. SIAM J. Numer. Anal. , 61(4):1918--1937, 2023

    1918

  6. [6]

    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

    2017

  7. [7]

    Engquist, H

    B. Engquist, H. Holst, and O. Runborg. Multi-scale methods for wave propagation in heterogeneous media. Commun. Math. Sci. , 9:33--56, 2011

    2011

  8. [8]

    Ekanathan, O

    S. Ekanathan, O. Smith, and C. Rackauckas. A fully adaptive radau method for the efficient solution of stiff ordinary differential equations at low tolerances. In 2025 IEEE High Performance Extreme Computing Conference (HPEC) , pages 1--9, 2025

    2025

  9. [9]

    L. C. Evans. Partial Differential Equations , volume 19 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI, second edition edition, 2010

    2010

  10. [10]

    om -formeln f\

    E. Fehlberg. Klassische Runge-Kutta-Nystr\"om -formeln f\"unfter und siebenter O rdnung mit S chrittweiten- K ontrolle. Computing , 4(2):93--106, 1969

    1969

  11. [11]

    C. W. Gear. Simultaneous numerical solution of differential-algebraic equations. IEEE Trans. Circuit Theory , 18(1):89--95, 1971

    1971

  12. [12]

    Geevers and R

    S. Geevers and R. Maier. Fast mass lumped multiscale wave propagation modelling. IMA J. Numer. Anal. , 43(1):44--72, 2023

    2023

  13. [13]

    Hauck, A

    M. Hauck, A. Lozinski, and R. Maier. A generalized framework for higher-order localized orthogonal decomposition methods. ESAIM Math. Model. Numer. Anal. , 60(1):445--471, 2026

    2026

  14. [14]

    Henning and D

    P. Henning and D. Peterseim. Oversampling for the multiscale finite element method. Multiscale Model. Simul. , 11(4):1149--1175, 2013

    2013

  15. [15]

    Hauck and D

    M. Hauck and D. Peterseim. Multi-resolution localized orthogonal decomposition for H elmholtz problems. Multiscale Model. Simul. , 20(2):657--684, 2022

    2022

  16. [16]

    Hairer and G

    E. Hairer and G. Wanner. Stiff differential equations solved by radau methods. J. Comput. Appl. Math. , 111(1-2):93–111, 1999

    1999

  17. [17]

    Jiang, Y

    L. Jiang, Y. Efendiev, and V. Ginting. Analysis of global multiscale finite element methods for wave equations with continuum spatial scales. Appl. Numer. Math. , 60:862--876, 2010

    2010

  18. [18]

    Optimal higher-order convergence rates for parabolic multiscale problems

    B. Kalyanaraman, F. Krumbiegel, R. Maier, and S. Wang. Optimal higher-order convergence rates for parabolic multiscale problems. ArXiv Preprint , 2510.09514, 2025

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  19. [19]

    Krumbiegel and R

    F. Krumbiegel and R. Maier. A higher order multiscale method for the wave equation. IMA J. Numer. Anal. , 45(4):2248--2273, 2025

    2025

  20. [20]

    Krumbiegel

    F. Krumbiegel. Higher-order numerical homogenization for wave propagation in strongly heterogeneous media . PhD thesis, Karlsruher Institut für Technologie (KIT), 2026

    2026

  21. [21]

    R. Maier. A high-order approach to elliptic multiscale problems with general unstructured coefficients. SIAM J. Numer. Anal. , 59(2):1067--1089, 2021

    2021

  22. [22]

    Montijano, Rández L., and Calvo M

    J.I. Montijano, Rández L., and Calvo M. Explicit runge–kutta–nyström methods for the numerical solution of second order linear inhomogeneous ivps. J. Comput. Appl. Math. , 438:115533, 2024

    2024

  23. [23]

    M lqvist and D

    A. M lqvist and D. Peterseim. Localization of elliptic multiscale problems. Math. Comp. , 83(290):2583--2603, 2014

    2014

  24. [24]

    Maier and D

    R. Maier and D. Peterseim. Explicit computational wave propagation in micro-heterogeneous media. BIT , 59(2):443--462, 2019

    2019

  25. [25]

    M lqvist and D

    A. M lqvist and D. Peterseim. Numerical homogenization by localized orthogonal decomposition , volume 5 of SIAM Spotlights . Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2020

    2020

  26. [26]

    Maier and B

    B. Maier and B. Verf\" u rth. Numerical upscaling for wave equations with time-dependent multiscale coefficients. Multiscale Model. Simul. , 20(4):1169--1190, 2022

    2022

  27. [27]

    Owhadi and C

    H. Owhadi and C. Scovel. Operator-adapted wavelets, fast solvers, and numerical homogenization , volume 35 of Cambridge Monographs on Applied and Computational Mathematics . Cambridge University Press, Cambridge, 2019

    2019

  28. [28]

    H. Owhadi. Multigrid with rough coefficients and multiresolution operator decomposition from hierarchical information games. SIAM Rev. , 59(1):99--149, 2017

    2017

  29. [29]

    Owhadi and L

    H. Owhadi and L. Zhang. Gamblets for opening the complexity-bottleneck of implicit schemes for hyperbolic and parabolic ODEs/PDEs with rough coefficients. J. Comput. Phys. , 347:99--128, 2017

    2017

  30. [30]

    Rackauckas and Q

    C. Rackauckas and Q. Nie. Differential E quations.jl--a performant and feature-rich ecosystem for solving differential equations in J ulia. J. Open Res. Software , 5(1), 2017

    2017

  31. [31]

    H. H. Rosenbrock. Some general implicit processes for the numerical solution of differential equations. Comput. J. , 5(4):329--330, 1963

    1963

  32. [32]

    L. F. Shampine and M. W. Reichelt. The MATLAB ODE suite. SIAM J. Sci. Comput. , 18(1):1--22, 1997

    1997

  33. [33]

    Steinebach

    G. Steinebach. Construction of rosenbrock–wanner method rodas5p and numerical benchmarks within the julia differential equations package. BIT , 63(2):27, 2023

    2023

  34. [34]

    Verdugo and S

    F. Verdugo and S. Badia. The software design of gridap: a finite element package based on the J ulia JIT compiler. Comput. Phys. Commun. , 276:108341, 2022

    2022