pith. sign in

arxiv: 2607.02362 · v1 · pith:4KGGE3TXnew · submitted 2026-07-02 · 🧮 math.NA · cs.NA

Sobolev stability of the L²-projection on hybrid meshes

Lars Diening , Viktoria Lingert , Tabea Tscherpel This is my paper

Pith reviewed 2026-07-03 07:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords L2-projectionSobolev stabilityhybrid meshesadaptive refinementfinite elementsconvex quadrilateralsLagrange elementsstability estimates
0
0 comments X

The pith

The L²-projection onto mapped Lagrange finite elements remains stable in Lp and W¹,p norms on hybrid meshes of triangles and convex quadrilaterals from adaptive refinement.

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 prove that the L²-projection onto the finite element space stays bounded when measured in Lp and first-order Sobolev Lp norms on meshes that combine triangles with convex quadrilaterals produced by local adaptive refinement. A reader would care because this bound is needed to carry standard error estimates over to these meshes without introducing factors that grow with element distortion. The result holds in particular for W¹,²-stability at every polynomial degree K at least two for two common quadrilateral refinement patterns, by extending a geometric argument that controls the constants on general convex quadrilaterals rather than only on parallelograms. If correct, the projection can be inserted into a priori and a posteriori analyses on such hybrid meshes for all degrees starting from two.

Core claim

The central claim is that the L²-projection onto mapped Lagrange finite elements on hybrid meshes is Lp- and W¹,p-stable in general, and in particular W¹,²-stable for all tensor-product degrees K greater than or equal to 2 on the two specified quadrilateral refinement patterns. This follows from extending a geometric control technique to arbitrary convex quadrilaterals that arise during adaptive refinement.

What carries the argument

Extension of a geometric technique that controls the constants appearing in the stability estimate for the L²-projection when the elements are general convex quadrilaterals rather than parallelograms.

If this is right

  • Standard a priori and a posteriori error estimates for finite element methods carry over directly to adaptively refined hybrid meshes without extra dependence on mesh distortion for degrees K at least 2.
  • The projection remains usable inside analyses of nonlinear problems that rely on W¹,p bounds.
  • Both Lp and W¹,p stability are obtained, so the result applies to a wider class of estimates than the Sobolev case alone.
  • The same meshes can be used for higher-order approximations without loss of the stability property needed for convergence proofs.

Where Pith is reading between the lines

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

  • If the geometric technique extends further, similar stability statements could be checked on three-dimensional hybrid meshes.
  • The same control might be reusable for other canonical projections that appear in mixed finite element methods on the same meshes.
  • Numerical tests that track the stability constant on a sequence of convex quadrilaterals with increasing aspect ratio would directly test the claimed independence from shape.

Load-bearing premise

The geometric technique that works for parallelograms continues to control the constants when the quadrilaterals become arbitrary convex shapes produced by adaptive refinement.

What would settle it

A concrete family of increasingly distorted convex quadrilaterals together with a fixed polynomial degree K=2 on which the ratio of the W¹,² norm of the projected function to the original function grows without bound.

Figures

Figures reproduced from arXiv: 2607.02362 by Lars Diening, Tabea Tscherpel, Viktoria Lingert.

Figure 1
Figure 1. Figure 1: Transformation of a triangle (left) and a convex quadri￾lateral (right). The shape-regularity constant of an element T ∈ T is defined as χ(T) := max xb∈Tb ∥∇BT (xb)∥ max x∈T ∥∇B −1 T (x)∥, (2.1) where ∥·∥ is the spectral norm which is induced by the Euclidean vector norm. Note that χ(T) ≥ 1, and that χ(T) increases with the deformation of T. The mapping BT is not unique, since the vertices can be rotated. … view at source ↗
Figure 2
Figure 2. Figure 2: Local refinement of [0, 1]2 (first line) and of a general convex quadrilateral (second line) as used in Q-RG and Q-RB re￾finement, see [AFS22] and [BSW83; Kob96] straight lines in T only if they are axis-parallel. Hence, mapping the new edges in [0, 1]2 does not lead to quadrilaterals with straight edges after mapping them. The global adaptive refinement routine uses red refinements, and temporary green re… view at source ↗
Figure 3
Figure 3. Figure 3: Shape-regularity constants of selected elements occur￾ring in the refinement of [0, 1]2 , see [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Convex set C as defined in (3.19) and truncated set Cτ as in (3.38). (b) is a general convex quadrilateral, i.e., (c1, c2) ⊤ ∈ C. Proposition 3.4. If T ∈ T is a parallelogram, then the smallest constant K1 in (3.9) is given by K1 = [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Decay rate q plotted over polynomial degree K for different types of meshes, see Remark 3.12. Note that Q-RG gen￾erates SP-hybrid and Q-RB generates general hybrid meshes for P-quadrilateral initial mesh, and that both generate general hy￾brid meshes for quadrilateral initial mesh, see [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
read the original abstract

We establish $L^p$- and $W^{1,p}$-stability of the $L^2$-projection onto mapped Lagrange finite elements on hybrid meshes consisting of triangles and convex quadrilaterals arising from adaptive mesh refinement. If $K$ is the (tensor product) degree of polynomials of the discretisation, then we show, in particular, $W^{1,2}$-stability for all $K\geq 2$ for the Q-RG and Q-RB refinements. This extends results by Ali, Funken, and Schmidt (2022) which hold for the range $2 \leq K \leq 9$ for initial meshes consisting of parallelograms. Our proof relies on an extension of the technique by Diening, Storn and Tscherpel (2021) to general convex quadrilaterals.

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. The manuscript establishes L^p- and W^{1,p}-stability of the L^2-projection onto mapped Lagrange finite elements on hybrid meshes consisting of triangles and convex quadrilaterals arising from adaptive mesh refinement. In particular, it proves W^{1,2}-stability for all K ≥ 2 for the Q-RG and Q-RB refinements by extending the technique of Diening, Storn and Tscherpel (2021) to general convex quadrilaterals, thereby improving on the range 2 ≤ K ≤ 9 obtained by Ali, Funken and Schmidt (2022) for parallelograms.

Significance. If the extension is carried out with uniform control of all constants, the result would provide useful stability guarantees for higher-order elements on adaptively refined hybrid meshes, which are common in practical FEM computations.

major comments (1)
  1. [Proof of the extension to general convex quadrilaterals (the section containing the adaptation of the 2021 technique)] The central claim of W^{1,2}-stability for K ≥ 2 on general convex quadrilaterals rests on the extension of the Diening-Storn-Tscherpel (2021) argument. On such elements the Jacobian of the bilinear map is non-constant, so the pullbacks of tensor-product polynomials are no longer polynomials and the norm equivalences acquire additional factors. The manuscript must demonstrate explicitly (with mesh-independent bounds) that these factors remain controlled uniformly for the Q-RG and Q-RB families; otherwise the stability constant may depend on the local geometry and the result fails to hold as stated.
minor comments (2)
  1. [Introduction] Clarify the precise definitions of the Q-RG and Q-RB refinement patterns at the beginning of the paper so that the reader can follow the geometric assumptions without searching later sections.
  2. [Main theorems] Ensure that all constants appearing in the stability estimates are stated to be independent of the mesh size and the local aspect ratios; this is implied but not always written explicitly in the statements.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive feedback. We address the single major comment below.

read point-by-point responses

  1. Referee: The central claim of W^{1,2}-stability for K ≥ 2 on general convex quadrilaterals rests on the extension of the Diening-Storn-Tscherpel (2021) argument. On such elements the Jacobian of the bilinear map is non-constant, so the pullbacks of tensor-product polynomials are no longer polynomials and the norm equivalences acquire additional factors. The manuscript must demonstrate explicitly (with mesh-independent bounds) that these factors remain controlled uniformly for the Q-RG and Q-RB families; otherwise the stability constant may depend on the local geometry and the result fails to hold as stated.

    Authors: We agree that an explicit verification of the uniform control on the additional factors is required for the claim to be fully rigorous. In extending the Diening–Storn–Tscherpel argument we track the Jacobian contributions arising from the bilinear mapping on convex quadrilaterals and invoke the geometric properties of the Q-RG and Q-RB families (convexity preservation together with angle and aspect-ratio bounds inherited from the initial mesh) to obtain K-dependent but mesh-independent constants. Nevertheless, these steps are currently distributed throughout the proof rather than isolated. We will therefore add a dedicated lemma in the revised manuscript that isolates the Jacobian-related equivalence constants and proves their uniformity for the two refinement families. revision: yes

Circularity Check

0 steps flagged

No circularity; proof provides the claimed extension

full rationale

The manuscript is a mathematical proof paper establishing Sobolev stability by extending a prior technique to convex quadrilaterals with bilinear maps. The abstract explicitly states that the proof consists of this extension, so the derivation chain is self-contained within the present work rather than reducing to an unverified self-citation. No self-definitional loops, fitted parameters renamed as predictions, or ansatzes smuggled via citation appear. The 2021 citation supplies the base technique being extended, which is standard and does not create load-bearing circularity when the extension is proved here. This is the normal, non-circular outcome for a direct proof paper.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Based on abstract only; the result rests on extending a 2021 technique to general convex quads and standard FEM assumptions. No free parameters or invented entities are indicated.

axioms (2)
  • domain assumption The technique by Diening, Storn and Tscherpel (2021) extends to general convex quadrilaterals
    Key extension claimed to achieve the result for hybrid meshes.
  • standard math Standard properties of mapped Lagrange finite elements and Sobolev spaces on convex domains
    Background assumptions for L^p and W^{1,p} stability proofs.

pith-pipeline@v0.9.1-grok · 5680 in / 1448 out tokens · 35107 ms · 2026-07-03T07:25:50.200887+00:00 · methodology

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

36 extracted references · 31 canonical work pages · 1 internal anchor

  1. [1]

    10.1051/m2an/197509R200771

    Approximation by finite element functions using local regularization , DOI= "10.1051/m2an/197509R200771", url= "https://doi.org/10.1051/m2an/197509R200771", journal =

    work page doi:10.1051/m2an/197509r200771

  2. [2]

    Ciarlet.The Finite Element Method for Elliptic Problems

    Ciarlet, P. , publisher =. 2002 , isbn =. doi:10.1137/1.9780898719208 , owner =

    work page doi:10.1137/1.9780898719208 2002

  3. [3]

    Interpolatory Subdivision on Open Quadrilateral Nets with Arbitrary Topology , year =

    Kobbelt, Leif , journal =. Interpolatory Subdivision on Open Quadrilateral Nets with Arbitrary Topology , year =. doi:10.1111/1467-8659.1530409 , fjournal =

    work page doi:10.1111/1467-8659.1530409

  4. [4]

    , journal =

    Carstensen, C. , journal =. Merging the. 2002 , issn =. doi:10.1090/S0025-5718-01-01316-3 , fjournal =

    work page doi:10.1090/s0025-5718-01-01316-3 2002

  5. [5]

    and Thom\'ee, V

    Crouzeix, M. and Thom\'ee, V. , journal =. The Stability in. 1987 , number =. doi:10.2307/2007825 , fjournal =

    work page doi:10.2307/2007825 1987

  6. [6]

    Generalized convexity and optimization , year =

    Cambini, Alberto and Martein, Laura , publisher =. Generalized convexity and optimization , year =

  7. [7]

    and Pasciak, Joseph E

    Bramble, James H. and Pasciak, Joseph E. and Steinbach, Olaf , journal =. On the stability of the. 2002 , issn =. doi:10.1090/S0025-5718-01-01314-X , fjournal =

    work page doi:10.1090/s0025-5718-01-01314-x 2002

  8. [8]

    and Sherman, Andrew H

    Bank, Randolph E. and Sherman, Andrew H. and Weiser, Alan , booktitle =. Refinement algorithms and data structures for regular local mesh refinement , year =

  9. [9]

    An adaptive mesh-refining algorithm allowing for an

    Carstensen, Carsten , journal =. An adaptive mesh-refining algorithm allowing for an. 2004 , issn =. doi:10.1007/s00365-003-0550-5 , fjournal =

    work page doi:10.1007/s00365-003-0550-5 2004

  10. [10]

    doi:10.1090/mcom/4102 , fjournal =

    Diening, Lars and Storn, Johannes and Tscherpel, Tabea , title =. doi:10.1090/mcom/4102 , fjournal =

    work page doi:10.1090/mcom/4102

  11. [11]

    and Boffi, Daniele and Falk, Richard S

    Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S. , journal =. Approximation by quadrilateral finite elements , year =. doi:10.1090/S0025-5718-02-01439-4 , fjournal =

    work page doi:10.1090/s0025-5718-02-01439-4

  12. [12]

    On the finite element method on quadrilateral meshes , year =

    Boffi, Daniele , journal =. On the finite element method on quadrilateral meshes , year =. doi:10.1016/j.apnum.2006.03.027 , fjournal =

    work page doi:10.1016/j.apnum.2006.03.027 2006

  13. [13]

    and Boffi, Daniele and Falk, Richard S

    Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S. and Gastaldi, Lucia , journal =. Finite element approximation on quadrilateral meshes , year =. doi:10.1002/cnm.450 , fjournal =

    work page doi:10.1002/cnm.450

  14. [14]

    Breit, Dominic and Diening, Lars and Storn, Johannes and Wichmann, J\". IMA J. Numer. Anal. , title =. 2021 , issn =. doi:10.1093/imanum/draa081 , fjournal =

    work page doi:10.1093/imanum/draa081 2021

  15. [15]

    Bramble, J. H. and Xu, J. , journal =. Some estimates for a weighted. 1991 , issn =. doi:10.2307/2008391 , fjournal =

    work page doi:10.2307/2008391 1991

  16. [16]

    , journal =

    Steinbach, O. , journal =. On the stability of the. 2001 , issn =. doi:10.1007/PL00005449 , fjournal =

    work page doi:10.1007/pl00005449 2001

  17. [17]

    , journal =

    Steinbach, O. , journal =. On a generalized. 2002 , issn =. doi:10.1007/s002110100329 , fjournal =

    work page doi:10.1007/s002110100329 2002

  18. [18]

    and Heine, Claus-Justus and Siebert, Kunibert G

    Gaspoz, Fernando D. and Heine, Claus-Justus and Siebert, Kunibert G. , journal =. Optimal grading of the newest vertex bisection and. 2016 , issn =. doi:10.1093/imanum/drv044 , fjournal =

    work page doi:10.1093/imanum/drv044 2016

  19. [19]

    Gaspoz, F. D. and Heine, C.-J. and Siebert, K. G. , note =. An Alternative Proof of the H^1 -Stability of the L_2 -Projection on Graded Meshes , year =. https://www.f08.uni-stuttgart.de/.content/media/downloads/Mathematik/mathematische_berichte/2019/2019-001.pdf , owner =

    2019

  20. [20]

    Maubach, J. M. , journal =. Local bisection refinement for. 1995 , issn =. doi:10.1137/0916014 , fjournal =

    work page doi:10.1137/0916014 1995

  21. [21]

    Anisotropic finite elements: local estimates and applications , year =

    Apel, Thomas , publisher =. Anisotropic finite elements: local estimates and applications , year =

  22. [22]

    Adaptive finite element methods on quadrilateral meshes without hanging nodes , year =

    Zhao, Xuying and Mao, Shipeng and Shi, Zhongci , journal =. Adaptive finite element methods on quadrilateral meshes without hanging nodes , year =. doi:10.1137/090772022 , fjournal =

    work page doi:10.1137/090772022

  23. [23]

    and Girault, V

    Bernardi, C. and Girault, V. , journal =. A local regularization operator for triangular and quadrilateral finite elements , year =. doi:10.1137/S0036142995293766 , fjournal =

    work page doi:10.1137/s0036142995293766

  24. [24]

    and Melenk, J

    Apel, T. and Melenk, J. M. , pages =. Interpolation and Quasi-Interpolation in h - and hp - Version Finite Element Spaces , year =. Encyclopedia of Computational Mechanics Second Edition , doi =

  25. [25]

    Exponential convergence of

    Faustmann, Markus and Marcati, Carlo and Melenk, Jens Markus and Schwab, Christoph , journal =. Exponential convergence of. 2023 , issn =. doi:10.1137/22M152493X , fjournal =

    work page doi:10.1137/22m152493x 2023

  26. [26]

    Theory of two-level Schwarz preconditioners with piecewise-polynomial coarse spaces for the high-frequency Helmholtz equation

    Ivan G. Graham and Euan A. Spence , title =. 2025 , archiveprefix =. doi:10.48550/arXiv.2501.15976 , url =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2501.15976 2025

  27. [27]

    and Guermond, J.-L

    Ern, A. and Guermond, J.-L. , publisher =. Finite elements. 2021 , isbn =. doi:10.1007/978-3-030-56341-7 , pages =

    work page doi:10.1007/978-3-030-56341-7 2021

  28. [28]

    Ali, Mazen and Funken, Stefan A and Schmidt, Anja , journal =. H^. 2022 , issn =. doi:10.1093/imanum/drab048 , publisher =

    work page doi:10.1093/imanum/drab048 2022

  29. [29]

    Estimates for the

    Boman, Mats , journal =. Estimates for the. 2006 , issn =. doi:10.1007/s10543-006-0062-3 , fjournal =

    work page doi:10.1007/s10543-006-0062-3 2006

  30. [30]

    and Yserentant, Harry , journal =

    Bank, Randolph E. and Yserentant, Harry , journal =. On the. 2014 , issn =. doi:10.1007/s00211-013-0562-4 , fjournal =

    work page doi:10.1007/s00211-013-0562-4 2014

  31. [31]

    Diening, Lars and Storn, Johannes and Tscherpel, Tabea , journal =. On the. 2021 , issn =. doi:10.1137/20M1358013 , fjournal =

    work page doi:10.1137/20m1358013 2021

  32. [32]

    Adaptive Finite Element Methods for Parabolic Problems

    Eriksson, Kenneth and Johnson, Claes , journal =. Adaptive Finite Element Methods for Parabolic Problems. 1995 , issn =. doi:10.1137/0732033 , fjournal =

    work page doi:10.1137/0732033 1995

  33. [33]

    Ridgway and Zhang, Shangyou , journal =

    Scott, L. Ridgway and Zhang, Shangyou , journal =. Finite element interpolation of nonsmooth functions satisfying boundary conditions , year =. doi:10.2307/2008497 , fjournal =

    work page doi:10.2307/2008497

  34. [34]

    Traxler, C. T. , journal =. An algorithm for adaptive mesh refinement in. 1997 , issn =. doi:10.1007/BF02684475 , fjournal =

    work page doi:10.1007/bf02684475 1997

  35. [35]

    Tantardini, Francesca and Veeser, Andreas , journal =. The. 2016 , issn =. doi:10.1137/140996811 , fjournal =

    work page doi:10.1137/140996811 2016

  36. [36]

    , journal =

    Knupp, Patrick M. , journal =. Algebraic mesh quality metrics , year =. doi:10.1137/S1064827500371499 , fjournal =

    work page doi:10.1137/s1064827500371499