pith. sign in

arxiv: 2504.18322 · v2 · submitted 2025-04-25 · 🧮 math.NA · cs.NA

Stable localized orthogonal decomposition in Raviart-Thomas spaces

Patrick Henning , Hao Li , Timo Sprekeler This is my paper

Pith reviewed 2026-05-22 17:52 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords localized orthogonal decompositionRaviart-Thomas spacesmixed finite elementsmultiscale methodnumerical homogenizationelliptic equationsheterogeneous coefficientsfinite element method
0
0 comments X

The pith

Stable localized orthogonal decomposition in Raviart-Thomas spaces yields multiscale approximations for mixed elliptic problems free of pollution terms.

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

The paper introduces a multiscale method for the mixed formulation of second-order linear elliptic equations with heterogeneous coefficients and homogeneous Neumann boundary conditions. It employs a stable localized orthogonal decomposition directly in Raviart-Thomas finite element spaces, constructing low-dimensional coarse approximation spaces by solving local patch problems on a fine mesh. This incorporates fine-scale information from the coefficients without requiring periodicity or scale separation and applies in two and three dimensions. The approach includes a rigorous error analysis and specifically avoids pollution terms that appeared in earlier LOD constructions for mixed problems. Numerical experiments illustrate the method's performance and convergence behavior.

Core claim

By realizing a stable localized orthogonal decomposition in Raviart-Thomas spaces, the method generates coarse approximation spaces that capture fine-scale coefficient variations through local computations, enabling accurate solutions to the mixed elliptic problem while eliminating pollution effects seen in previous constructions.

What carries the argument

Stable localized orthogonal decomposition (LOD) in Raviart-Thomas finite element spaces, which solves local patch problems to transfer fine-scale coefficient information stably into the coarse spaces.

If this is right

  • Coarse-scale approximations achieve optimal convergence rates independent of fine-scale oscillations.
  • The method applies to general heterogeneous coefficients in two and three dimensions without scale separation assumptions.
  • Rigorous a priori error estimates bound the approximation error for the mixed problem.
  • Numerical experiments confirm the theoretical predictions and practical efficiency of the scheme.
  • The absence of pollution terms improves reliability compared to prior LOD variants for mixed formulations.

Where Pith is reading between the lines

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

  • This approach could enable more efficient modeling of fluid flow through heterogeneous porous media using mixed finite element methods.
  • The stability property might facilitate extensions to adaptive mesh refinement or time-dependent problems by reusing the local patch constructions.
  • Similar LOD techniques could be tested on other mixed finite element spaces or for problems with different boundary conditions.
  • The pollution-free property suggests potential gains in efficiency when coupling to iterative solvers at the coarse scale.

Load-bearing premise

The local patch problems solved on the fine mesh are assumed to fully capture and transfer the fine-scale coefficient information into the coarse Raviart-Thomas spaces without introducing instability or requiring additional stabilization.

What would settle it

Numerical experiments on a highly oscillatory coefficient problem that reveal persistent pollution terms or failure of the expected convergence rates in the mixed formulation would disprove the central stability claim.

Figures

Figures reproduced from arXiv: 2504.18322 by Hao Li, Patrick Henning, Timo Sprekeler.

Figure 1
Figure 1. Figure 1: Plots of the coefficient and the magnitude of the reference flux for Experiment 1. [PITH_FULL_IMAGE:figures/full_fig_p023_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Convergence Test. Relative errors for the LOD approximations of Experiment 1. [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Plots of the coefficient and the magnitude of the reference flux for Experiment 2. [PITH_FULL_IMAGE:figures/full_fig_p024_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: SPE10-85 test. Figure 4d shows the magnitude of the reference flux solution. Figures 4a to 4c display the magnitudes of the multiscale flux solutions for m = 2, 3, 4 with ℓ = m+1, respectively. References [1] J. E. Aarnes. On the use of a mixed multiscale finite element method for greater flexibility and increased speed or improved accuracy in reservoir simulation. Multiscale Model. Simul., 2(3):421– 439, … view at source ↗
read the original abstract

This work proposes a computational multiscale method for the mixed formulation of a second-order linear elliptic equation subject to a homogeneous Neumann boundary condition, based on a stable localized orthogonal decomposition (LOD) in Raviart-Thomas finite element spaces. In the spirit of numerical homogenization, the construction provides low-dimensional coarse approximation spaces that incorporate fine-scale information from the heterogeneous coefficients by solving local patch problems on a fine mesh. The resulting numerical scheme is accompanied by a rigorous error analysis, and it is applicable beyond periodicity and scale-separation in spatial dimensions two and three. In particular, this novel realization circumvents the presence of pollution terms observed in a previous LOD construction for elliptic problems in mixed formulation. Finally, various numerical experiments are provided that demonstrate the performance of the method.

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 / 2 minor

Summary. This paper claims to develop a stable localized orthogonal decomposition (LOD) in Raviart-Thomas finite element spaces for the mixed formulation of second-order linear elliptic equations with homogeneous Neumann boundary conditions. By solving local patch problems on a fine mesh, it constructs coarse spaces that incorporate fine-scale information from heterogeneous coefficients. The method is supported by a rigorous a priori error analysis and numerical experiments, and it is asserted to avoid the pollution terms present in a previous mixed LOD construction, while being applicable in 2D and 3D without periodicity or scale separation.

Significance. If the result holds, this work offers a significant improvement in multiscale modeling for mixed elliptic problems by providing a pollution-free LOD approach in Raviart-Thomas spaces. This could enhance the accuracy of coarse-scale simulations for problems with highly varying coefficients, such as in subsurface flow modeling. The combination of theoretical analysis and numerical validation adds to its potential impact in the numerical analysis community.

major comments (1)
  1. [§4] §4 (Error Analysis): The central claim of stability without pollution terms rests on the a priori estimates. The analysis should explicitly derive how the Raviart-Thomas degrees of freedom in the local patch problems bound the consistency error independently of the coefficient contrast, as this underpins the absence of pollution relative to prior mixed LOD work.
minor comments (2)
  1. The introduction would benefit from a brief table or paragraph contrasting the new construction with the referenced previous mixed LOD method to make the novelty clearer.
  2. [Numerical Experiments] In the numerical experiments section, include observed convergence rates alongside the error tables for the different test cases.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the positive recommendation for minor revision. The single major comment concerns the clarity of the error analysis in Section 4, and we address it directly below. We will revise the manuscript to incorporate the requested explicit derivation.

read point-by-point responses

  1. Referee: [§4] §4 (Error Analysis): The central claim of stability without pollution terms rests on the a priori estimates. The analysis should explicitly derive how the Raviart-Thomas degrees of freedom in the local patch problems bound the consistency error independently of the coefficient contrast, as this underpins the absence of pollution relative to prior mixed LOD work.

    Authors: We appreciate the referee's suggestion to strengthen the presentation of the error analysis. In the current manuscript, Theorems 4.1 and 4.2 establish the a priori estimates for the stable LOD method in Raviart-Thomas spaces by constructing the coarse spaces through local patch problems on the fine mesh. The key property is that the Raviart-Thomas degrees of freedom in these local problems allow the consistency error to be controlled via the orthogonal decomposition without introducing contrast-dependent pollution terms, unlike earlier mixed LOD constructions. To make this explicit as requested, we will add a dedicated remark (or short subsection) in the revised Section 4 that walks through the bounding argument step by step, showing how the local RT degrees of freedom yield contrast-independent estimates on the consistency error. This addition will clarify the distinction from prior work and improve accessibility without altering the overall proof structure. We agree that this clarification is beneficial. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper constructs stable LOD coarse spaces in Raviart-Thomas elements by solving local patch problems on a fine mesh to transfer heterogeneous coefficient information, then performs a rigorous a priori error analysis for the mixed elliptic problem with homogeneous Neumann conditions. This approach is presented as a direct extension of standard finite-element localization techniques that avoids pollution terms seen in earlier mixed LOD formulations. No load-bearing step reduces by construction to fitted parameters, self-definitions, or unverified self-citation chains; the central claims rest on explicit local solves and independent error estimates that are externally falsifiable. The derivation is therefore self-contained against standard FEM benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method relies on standard assumptions from finite element theory and numerical homogenization; no free parameters, ad-hoc axioms, or new invented entities are indicated in the abstract.

axioms (1)
  • standard math Standard properties of Raviart-Thomas finite element spaces and localized orthogonal decomposition hold for the mixed elliptic problem.
    Invoked implicitly when constructing the coarse spaces and claiming stability.

pith-pipeline@v0.9.0 · 5651 in / 1303 out tokens · 38062 ms · 2026-05-22T17:52:28.672946+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

53 extracted references · 53 canonical work pages

  1. [1]

    J. E. Aarnes. On the use of a mixed multiscale finite element method for greater flexibility and increased speed or improved accuracy in reservoir simulation. Multiscale Model. Simul., 2(3):421– 439, 2004

    work page 2004

  2. [2]

    Abdulle, W

    A. Abdulle, W. E, B. Engquist, and E. Vanden-Eijnden. The heterogeneous multiscale method. Acta Numer., 21:1–87, 2012

    work page 2012

  3. [3]

    Alber, C

    C. Alber, C. Ma, and R. Scheichl. A mixed multiscale spectral generalized finite element method. Numer. Math., 157(1):1–40, 2025

    work page 2025

  4. [4]

    Altmann, P

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

    work page 2021

  5. [5]

    Arbogast

    T. Arbogast. Analysis of a two-scale, locally conservative subgrid upscaling for elliptic problems. SIAM J. Numer. Anal. , 42(2):576–598, 2004

    work page 2004

  6. [6]

    Arbogast

    T. Arbogast. Homogenization-based mixed multiscale finite elements for problems with anisotropy. Multiscale Model. Simul., 9(2):624–653, 2011

    work page 2011

  7. [7]

    Arbogast and K

    T. Arbogast and K. J. Boyd. Subgrid upscaling and mixed multiscale finite elements. SIAM J. Numer. Anal., 44(3):1150–1171, 2006

    work page 2006

  8. [8]

    D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications. Acta Numer., 15:1–155, 2006

    work page 2006

  9. [9]

    D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus: from Hodge theory to numerical stability. Bull. Amer. Math. Soc. (N.S.) , 47(2):281–354, 2010

    work page 2010

  10. [10]

    Babuska and R

    I. Babuska and R. Lipton. Optimal local approximation spaces for generalized finite element methods with application to multiscale problems. Multiscale Model. Simul., 9(1):373–406, 2011

    work page 2011

  11. [11]

    Babuˇ ska, R

    I. Babuˇ ska, R. Lipton, P. Sinz, and M. Stuebner. Multiscale-spectral GFEM and optimal over- sampling. Comput. Methods Appl. Mech. Engrg. , 364:112960, 28, 2020

    work page 2020

  12. [12]

    Boffi, F

    D. Boffi, F. Brezzi, and M. Fortin. Mixed finite element methods and applications , volume 44 of Springer Series in Computational Mathematics . Springer, Heidelberg, 2013

    work page 2013

  13. [13]

    Chen and T

    Z. Chen and T. Y. Hou. A mixed multiscale finite element method for elliptic problems with oscillating coefficients. Math. Comp., 72(242):541–576, 2003

    work page 2003

  14. [14]

    S. W. Cheung, E. Chung, Y. Efendiev, W. T. Leung, and S.-M. Pun. Iterative oversampling technique for constraint energy minimizing generalized multiscale finite element method in the mixed formulation. Appl. Math. Comput. , 415:Paper No. 126622, 17, 2022

    work page 2022

  15. [15]

    M. A. Christie and M. J. Blunt. Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques. SPE Reservoir Evaluation & Engineering , 4(04):308–317, 2001

    work page 2001

  16. [16]

    Chung, Y

    E. Chung, Y. Efendiev, and T. Y. Hou. Multiscale model reduction—multiscale finite element methods and their generalizations, volume 212 of Applied Mathematical Sciences. Springer, Cham, 2023

    work page 2023

  17. [17]

    E. T. Chung, Y. Efendiev, and C. S. Lee. Mixed generalized multiscale finite element methods and applications. Multiscale Model. Simul., 13(1):338–366, 2015

    work page 2015

  18. [18]

    D¨ oding, P

    C. D¨ oding, P. Henning, and J. W¨ arneg˚ ard. A two level approach for simulating Bose-Einstein con- densates by localized orthogonal decomposition. ESAIM Math. Model. Numer. Anal., 58(6):2317– 2349, 2024

    work page 2024

  19. [19]

    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

    work page 1918

  20. [20]

    W. E and B. Engquist. The heterogeneous multiscale methods. Commun. Math. Sci., 1(1):87–132, 2003

    work page 2003

  21. [21]

    Efendiev, J

    Y. Efendiev, J. Galvis, and T. Y. Hou. Generalized multiscale finite element methods (GMsFEM). J. Comput. Phys. , 251:116–135, 2013

    work page 2013

  22. [22]

    Efendiev and T

    Y. Efendiev and T. Y. Hou. Multiscale finite element methods, volume 4 of Surveys and Tutorials in the Applied Mathematical Sciences . Springer, New York, 2009. Theory and applications

    work page 2009

  23. [23]

    A. Ern, T. Gudi, I. Smears, and M. Vohral´ ık. Equivalence of local- and global-best approximations, a simple stable local commuting projector, and optimal hp approximation estimates in H(div). IMA J. Numer. Anal. , 42(2):1023–1049, 2022

    work page 2022

  24. [24]

    Ern and J.-L

    A. Ern and J.-L. Guermond. Finite elements II—Galerkin approximation, elliptic and mixed PDEs, volume 73 of Texts in Applied Mathematics . Springer, Cham, 2021. 25

    work page 2021

  25. [25]

    Freese, D

    P. Freese, D. Gallistl, D. Peterseim, and T. Sprekeler. Computational multiscale methods for nondivergence-form elliptic partial differential equations. Comput. Methods Appl. Math. , 24(3):649–672, 2024

    work page 2024

  26. [26]

    Gallistl, P

    D. Gallistl, P. Henning, and B. Verf¨ urth. Numerical homogenization of H(curl)-problems. SIAM J. Numer. Anal. , 56(3):1570–1596, 2018

    work page 2018

  27. [27]

    Gallistl, T

    D. Gallistl, T. Sprekeler, and E. S¨ uli. Mixed finite element approximation of periodic Hamilton- Jacobi-Bellman problems with application to numerical homogenization.Multiscale Model. Simul., 19(2):1041–1065, 2021

    work page 2021

  28. [28]

    Girault and P.-A

    V. Girault and P.-A. Raviart. Finite element methods for Navier-Stokes equations , volume 5 of Springer Series in Computational Mathematics . Springer-Verlag, Berlin, 1986. Theory and algorithms

    work page 1986

  29. [29]

    Hauck and A

    M. Hauck and A. Lozinski. A localized orthogonal decomposition method for heterogeneous stokes problems. arXiv preprint arXiv:2410.14514 , 2024

    work page arXiv 2024

  30. [30]

    Hauck and D

    M. Hauck and D. Peterseim. Super-localization of elliptic multiscale problems. Math. Comp. , 92(341):981–1003, 2023

    work page 2023

  31. [31]

    Hellman, P

    F. Hellman, P. Henning, and A. M˚ alqvist. Multiscale mixed finite elements. Discrete Contin. Dyn. Syst. Ser. S , 9(5):1269–1298, 2016

    work page 2016

  32. [32]

    R. Helmig. Multiphase flow and transport processes in the subsurface: A contribution to the modeling of hydrosystems. Springer Berlin Heidelberg, 1997

    work page 1997

  33. [33]

    Henning and A

    P. Henning and A. Persson. Computational homogenization of time-harmonic Maxwell’s equations. SIAM J. Sci. Comput. , 42(3):B581–B607, 2020

    work page 2020

  34. [34]

    Henning and D

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

    work page 2013

  35. [35]

    Hiptmair and J

    R. Hiptmair and J. Xu. Nodal auxiliary space preconditioning in H(curl) and H(div) spaces. SIAM J. Numer. Anal. , 45(6):2483–2509, 2007

    work page 2007

  36. [36]

    T. Y. Hou and X.-H. Wu. A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys. , 134(1):169–189, 1997

    work page 1997

  37. [37]

    T. J. R. Hughes. Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods.Comput. Methods Appl. Mech. Engrg., 127(1-4):387–401, 1995

    work page 1995

  38. [38]

    T. J. R. Hughes, G. R. Feij´ oo, L. Mazzei, and J.-B. Quincy. The variational multiscale method— a paradigm for computational mechanics. Comput. Methods Appl. Mech. Engrg. , 166(1-2):3–24, 1998

    work page 1998

  39. [39]

    Kornhuber, D

    R. Kornhuber, D. Peterseim, and H. Yserentant. An analysis of a class of variational multiscale methods based on subspace decomposition. Math. Comp., 87(314):2765–2774, 2018

    work page 2018

  40. [40]

    M. G. Larson and A. M˚ alqvist. A mixed adaptive variational multiscale method with applications in oil reservoir simulation. Math. Models Methods Appl. Sci. , 19(7):1017–1042, 2009

    work page 2009

  41. [41]

    Ma and R

    C. Ma and R. Scheichl. Error estimates for discrete generalized FEMs with locally optimal spectral approximations. Math. Comp., 91(338):2539–2569, 2022

    work page 2022

  42. [42]

    C. Ma, R. Scheichl, and T. Dodwell. Novel design and analysis of generalized finite element methods based on locally optimal spectral approximations. SIAM J. Numer. Anal. , 60(1):244– 273, 2022

    work page 2022

  43. [43]

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

    work page 2021

  44. [44]

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

    work page 2021

  45. [45]

    M˚ alqvist

    A. M˚ alqvist. Multiscale methods for elliptic problems.Multiscale Model. Simul., 9(3):1064–1086, 2011

    work page 2011

  46. [46]

    M˚ alqvist and D

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

    work page 2014

  47. [47]

    M˚ alqvist and D

    A. M˚ alqvist and D. Peterseim.Numerical homogenization by localized orthogonal decomposition , volume 5 of SIAM Spotlights. Society for Industrial and Applied Mathematics (SIAM), Philadel- phia, PA, 2021

    work page 2021

  48. [48]

    M˚ alqvist and B

    A. M˚ alqvist and B. Verf¨ urth. An offline-online strategy for multiscale problems with random defects. ESAIM Math. Model. Numer. Anal. , 56(1):237–260, 2022

    work page 2022

  49. [49]

    P. A. Raviart and J. M. Thomas. A mixed finite element method for 2-nd order elliptic problems. In I. Galligani and E. Magenes, editors, Mathematical Aspects of Finite Element Methods , pages 292–315, Berlin, Heidelberg, 1977. Springer Berlin Heidelberg

    work page 1977

  50. [50]

    Sch¨ oberl

    J. Sch¨ oberl. A posteriori error estimates for Maxwell equations. Math. Comp., 77(262):633–649, 2008

    work page 2008

  51. [51]

    Verf¨ urth

    B. Verf¨ urth. Numerical homogenization for nonlinear strongly monotone problems. IMA J. Numer. Anal., 42(2):1313–1338, 2022

    work page 2022

  52. [52]

    Y. Wang, E. Chung, and L. Zhao. Constraint energy minimization generalized multiscale fi- nite element method in mixed formulation for parabolic equations. Math. Comput. Simulation , 188:455–475, 2021

    work page 2021

  53. [53]

    Wendland

    H. Wendland. Divergence-free kernel methods for approximating the Stokes problem. SIAM J. Numer. Anal., 47(4):3158–3179, 2009. A Inf-sup stability in classical Raviart–Thomas spaces In the following, we give a short proof of the classical inf-sup stability in (16). LetqH ∈ Qk H ∩L2 0(Ω)\{0} and let φqH ∈ H1(Ω) ∩ L2 0(Ω) be the unique solution to the Neuma...

    work page 2009