pith. sign in

arxiv: 2604.21887 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

Guaranteed inf-sup bounds and existence verification for semilinear elliptic problems via nonconforming finite elements

Pith reviewed 2026-05-09 20:38 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords inf-sup constantnonconforming finite elementsa posteriori verificationsemilinear elliptic problemsNavier-Stokes equationsNewton-Kantorovich argumentexistence verification
0
0 comments X

The pith

Nonconforming finite elements provide a guaranteed lower bound on the continuous inf-sup constant, enabling a posteriori verification of unique solutions for semilinear elliptic problems.

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

The paper develops a method to verify the existence of a unique regular solution near a computed numerical approximation using only quantities that can be computed from the discretisation itself. It relies on a Newton-Kantorovich-type argument that works for non-selfadjoint problems and extends previous theory to nonconforming finite element methods. The central step is deriving a guaranteed lower bound for the continuous inf-sup constant from a quasi-optimal nonconforming approximation, which also supports a new a priori error estimator. This allows all verification data to come from post-processing a single discrete solution, with proved convergence rates. The approach is tested on a fourth-order form of the stationary two-dimensional Navier-Stokes equations.

Core claim

A Newton--Kantorovich-type argument enables the a posteriori existence verification of a unique regular root near a computed approximation, purely from computable data. This framework allows for non-selfadjoint problems and extends the existing verification theory to nonconforming discretisations. A key ingredient is a guaranteed lower bound on the continuous inf-sup constant from a quasi-optimal nonconforming discretisation that enables a novel a priori error estimator. All quantities are obtained by post-processing a single discretisation; convergence rates are proved. The theory is applied to a fourth-order formulation of the stationary two-dimensional Navier--Stokes equations.

What carries the argument

Quasi-optimal nonconforming discretisation that yields a guaranteed lower bound on the continuous inf-sup constant and supports a novel a priori error estimator for the verification procedure.

Load-bearing premise

The Newton-Kantorovich-type argument must apply and the nonconforming discretisation must be quasi-optimal to produce the guaranteed inf-sup lower bound.

What would settle it

A concrete case in which the computed inf-sup lower bound is positive yet no unique regular solution exists near the approximation, or in which the nonconforming scheme fails to deliver a positive bound on a known well-posed problem.

Figures

Figures reproduced from arXiv: 2604.21887 by Benedikt Gr\"a{\ss}le.

Figure 1
Figure 1. Figure 1: Initial triangulation T0 for the unit square and L-shaped domain. 5.2 Academic example on the unit square This benchmark considers the source Fλ ≡ fλ ∈ L 2 (Ω) matching the smooth solution u(x, y) ≡ uλ(x, y) = λ x2 (1 − x) 2 y 2 (1 − y) 2 ∈ H 2 0 (Ω) for some parameter λ ∈ {1, 10, 100} on the unit square Ω = (0, 1)2 with initial triangula￾tion displayed in figure 1. As shown in figure 2, all three refineme… view at source ↗
Figure 2
Figure 2. Figure 2: Convergence history plots of the error |||u − unc|||pw, the error estimator η, and the guaranteed error bound ϱ− from theorem 2.7.a under the three refinement strategies λ = 1 (left) and λ = 100 (right) in subsection 5.2. uniform adaptive+reduction of hmax adaptive 100 101 102 103 104 105 106 107 10−4 10−2 100 102 1 1 1 0.5 1 0.3 ndof Cβ,3 κ(T , 1) µres µb |||(1−J)unc|||pw 100 101 102 103 104 105 106 107 0… view at source ↗
Figure 3
Figure 3. Figure 3: Convergence history plot of several quantities of table 1 for [PITH_FULL_IMAGE:figures/full_fig_p027_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Convergence history of the (lower bounds on the) inf-sup constants under the [PITH_FULL_IMAGE:figures/full_fig_p028_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Convergence history plot of different error measures (top left), the (lower [PITH_FULL_IMAGE:figures/full_fig_p029_5.png] view at source ↗
read the original abstract

A Newton--Kantorovich-type argument enables the a posteriori existence verification of a unique regular root near a computed approximation, purely from computable data. This framework allows for non-selfadjoint problems and extends the existing verification theory to nonconforming discretisations. A key ingredient is a guaranteed lower bound on the continuous inf-sup constant from a quasi-optimal nonconforming discretisation that enables a novel a priori error estimator. All quantities are obtained by post-processing a single discretisation; convergence rates are proved. The theory is applied to a fourth-order formulation of the stationary two-dimensional Navier--Stokes equations and illustrated by numerical experiments.

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

0 major / 3 minor

Summary. The paper develops a Newton-Kantorovich-type a posteriori framework for verifying existence and uniqueness of a unique regular root for semilinear elliptic problems near a computed approximation, using only computable quantities. A central ingredient is a guaranteed lower bound on the continuous inf-sup constant extracted from a quasi-optimal nonconforming finite element discretization; this bound also produces a novel a priori error estimator. The theory extends prior verification results to non-selfadjoint operators and nonconforming spaces, proves convergence rates, and is applied to a fourth-order formulation of the stationary two-dimensional Navier-Stokes equations with supporting numerical experiments.

Significance. If the derivations hold, the work supplies a practical, computable route to rigorous existence verification for nonlinear problems where conforming methods are inconvenient or unavailable. The explicit construction of the inf-sup lower bound from a single discretization, the mesh-regularity hypotheses that guarantee positivity independent of mesh size, and the proved convergence rates under those hypotheses are concrete strengths that enhance reliability in applications such as fluid dynamics.

minor comments (3)
  1. [§2.2] §2.2, definition of the nonconforming space: the quasi-optimality statement would be easier to verify if the precise mesh-regularity assumption (e.g., shape-regularity constant) were stated explicitly rather than referenced to an external lemma.
  2. [Table 1] Table 1 (Navier-Stokes example): the reported inf-sup lower bounds decrease slightly with refinement; a short remark explaining whether this is consistent with the proved mesh-independent positivity would prevent reader confusion.
  3. [Abstract / Introduction] The abstract claims 'convergence rates are proved' but does not indicate the norms or the order; adding one sentence in the introduction summarizing the rates would improve accessibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, the accurate summary of its contributions, and the recommendation for minor revision. The referee's description correctly identifies the key elements: the Newton-Kantorovich a posteriori verification framework, the guaranteed inf-sup lower bound from a nonconforming discretization, the extension to non-selfadjoint operators, the convergence rates, and the application to the stationary Navier-Stokes equations.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper obtains a guaranteed lower bound on the continuous inf-sup constant via post-processing of a single quasi-optimal nonconforming discretization, then applies a Newton-Kantorovich argument for a posteriori existence verification of a unique regular root. Both steps rest on explicitly stated mesh-regularity hypotheses, standard functional-analytic estimates, and proved a priori error bounds that are independent of the target existence result. No equation or claim reduces by construction to a fitted parameter, self-definition, or self-citation chain; the cited verification theory is extended rather than presupposed, and the Navier-Stokes example satisfies the stated regularity without internal inconsistency.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard functional analysis results for inf-sup conditions and the Newton-Kantorovich theorem in Banach spaces, plus domain assumptions on the semilinear elliptic operator and the nonconforming finite element spaces.

axioms (2)
  • domain assumption The continuous problem satisfies the conditions for the Newton-Kantorovich theorem to apply near the approximation.
    Invoked to guarantee existence and uniqueness of the regular root from computable data.
  • domain assumption The nonconforming discretization is quasi-optimal for the inf-sup constant.
    Central to obtaining the guaranteed lower bound and the a priori estimator.

pith-pipeline@v0.9.0 · 5401 in / 1431 out tokens · 35939 ms · 2026-05-09T20:38:33.007629+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

52 extracted references · 52 canonical work pages

  1. [1]

    T. Aubin. Probl \`e mes isop \'e rim \'e triques et espaces de Sobolev . J. Differential Geom. , 11(4), January 1976

  2. [2]

    Bertrand, C

    F. Bertrand, C. Carstensen, B. Gr \"a le, and N. T. Tran. Stabilization-free HHO a posteriori error control. Numer. Math. , 154(3-4):369--408, August 2023

  3. [3]

    Bespalov, A

    A. Bespalov, A. Haberl, and D. Praetorius. Adaptive FEM with coarse initial mesh guarantees optimal convergence rates for compactly perturbed elliptic problems. Comput. Methods Appl. Mech. Eng. , 317:318--340, April 2017

  4. [4]

    S. C. Brenner, M. Neilan, A. Reiser, and L.-Y. Sung. A C ^0 interior penalty method for a von K \'a rm \'a n plate. Numer. Math. , 135(3):803--832, March 2017

  5. [5]

    D. Braess. Finite Elements : Theory , Fast Solvers , and Applications in Elasticity Theory . Cambridge University Press, Cambridge, 3 edition, 2007

  6. [6]

    S. C. Brenner and L. R. Scott. The Mathematical Theory of Finite Element Methods , volume 15 of Texts in Applied Mathematics . Springer New York, New York, NY, 2008

  7. [7]

    Carstensen and D

    C. Carstensen and D. Gallistl. Guaranteed lower eigenvalue bounds for the biharmonic equation. Numer. Math. , 126(1):33--51, January 2014

  8. [8]

    Carstensen and B

    C. Carstensen and B. Gr \"a le. Rate-optimal higher-order adaptive conforming FEM for biharmonic eigenvalue problems on polygonal domains. Comput. Methods Appl. Mech. Eng. , 425:116931, May 2024

  9. [9]

    Carstensen and B

    C. Carstensen and B. Gr \"a le. Adaptive Morley FEM for 2D stationary Navier-Stokes . Math. Comp. , 95(358):613--645, March 2025

  10. [10]

    Carstensen and B

    C. Carstensen and B. Gr \"a le. Optimal convergence rates of adaptive WOPSIP from superclose nonconforming FEM . Math. Comp. , March 2025

  11. [11]

    P. G. Ciarlet, L. Gratie, and S. Kesavan. On the generalized von Karman equations and their approximation. Math. Models Methods Appl. Sci. , 17(04):617--633, April 2007

  12. [12]

    Carstensen, B

    C. Carstensen, B. Gr \"a le, and N. Nataraj. A posteriori error control for fourth-order semilinear problems with quadratic nonlinearity. SIAM J. Numer. Anal. , 62(2):919--945, April 2024

  13. [13]

    Carstensen, B

    C. Carstensen, B. Gr \"a le, and N. Nataraj. Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation. J. Numer. Math. , 32(1):77--109, 2024

  14. [14]

    P. G. Ciarlet. The Finite Element Method for Elliptic Problems . Classics in Applied Mathematics . Society for Industrial and Applied Mathematics , January 2002

  15. [15]

    P. G. Ciarlet and C. Mardare. On the Newton-Kantorovich theorem. Anal. Appl. , 10(03):249--269, July 2012

  16. [16]

    Carstensen, G

    C. Carstensen, G. Mallik, and N. Nataraj. A priori and a posteriori error control of discontinuous Galerkin finite element methods for the von K \'a rm \'a n equations. IMA J. Numer. Anal. , 39:167--200, 2019

  17. [17]

    Carstensen, G

    C. Carstensen, G. Mallik, and N. Nataraj. Nonconforming finite element discretization for semilinear problems with trilinear nonlinearity. IMA J. Numer. Anal. , 41(1):164--205, January 2021

  18. [18]

    Cayco and R

    M. Cayco and R. Nicolaides. Analysis of nonconforming stream function and pressure finite element spaces for the Navier-Stokes equations. Comput. Math. Appl. , 18(8):745--760, 1989

  19. [19]

    Carstensen and N

    C. Carstensen and N. Nataraj. A priori and a posteriori error analysis of the Crouzeix -- Raviart and Morley FEM with original and modified right-hand sides. Comput. Methods Appl. Math. , 21(2):289--315, 2021

  20. [20]

    Carstensen and N

    C. Carstensen and N. Nataraj. Lowest-order equivalent nonstandard finite element methods for biharmonic plates. ESAIM: M2AN , 56(1):41--78, January 2022

  21. [21]

    Carstensen, N

    C. Carstensen, N. Nataraj, G. C. Remesan, and D. Shylaja. Unified a priori analysis of four second-order FEM for fourth-order quadratic semilinear problems. Numer. Math. , 154(3-4):323--368, August 2023

  22. [22]

    Foias and G

    C. Foias and G. Prodi. Sur le comportement global des solutions non-stationnaires des \'e quations de Navier-Stokes en dimension 2 . Rend. Semin. Mat. Univ. Padova , 39:1--34, 1967

  23. [23]

    Foias and R

    C. Foias and R. Temam. Structure of the set of stationary solutions of the Navier - Stokes equations. Comm. Pure Appl. Math. , 30(2):149--164, March 1977

  24. [24]

    Gallistl

    D. Gallistl. Morley finite element method for the eigenvalues of the biharmonic operator. IMA J. Numer. Anal. , 35(4):1779--1811, October 2015

  25. [25]

    K. Geba. The Leray-Schauder degree and framed bordism. Th \'e orie Points Fixes Ses Appl. \`A L'analyse , 1975

  26. [26]

    Girault and P.-A

    V. Girault and P.-A. Raviart. Finite Element Methods for Navier-Stokes Equations : Theory and Algorithms . Number 5 in Springer series in computational mathematics. Springer, Berlin, softcover reprint of the hardcover 1. ed. 1986 edition, 1986

  27. [27]

    Gr \"a le

    B. Gr \"a le. AFEM .jl (v0.1). Julia library for adaptive finite element methods, October 2024

  28. [28]

    Grisvard

    P. Grisvard. Singularities in Boundary Value Problems . Number 22 in Research notes in applied mathematics. Masson, Paris, 1992

  29. [29]

    Haegeman

    J. Haegeman. KrylovKit (v0.7.1). https://doi.org/10.5281/zenodo.10884302, 2024

  30. [30]

    Jamet and P

    P. Jamet and P. A. Raviart. Numerical solution of the stationary Navier-Stokes equations by finite element methods. In Computing Methods in Applied Sciences and Engineering Part 1 , volume 10, pages 193--223. Springer Berlin Heidelberg, Berlin, Heidelberg, 1974

  31. [31]

    G. H. Knightly. An existence theorem for the von K \'a rm \'a n equations. Arch. Ration. Mech. Anal. , 27(3):233--242, 1967

  32. [32]

    X. Liu, M. T. Nakao, and S. Oishi. Computer-assisted proof for the stationary solution existence of the Navier -- Stokes equation over 3D domains. Commun. Nonlinear Sci. Numer. Simul. , 108:106223, May 2022

  33. [33]

    Liu and S

    X. Liu and S. Oishi. Verified eigenvalue evaluation for the Laplacian over polygonal domains of arbitrary shape. SIAM J. Numer. Anal. , 51(3):1634--1654, January 2013

  34. [34]

    Liao, Y.-C

    S.-K. Liao, Y.-C. Shu, and X. Liu. Optimal estimation for the Fujino -- Morley interpolation error constants. Japan J. Indust. Appl. Math. , 36(2):521--542, July 2019

  35. [35]

    R. Ma. The constants for the L ^2 error estimates of the interpolation of Morley-Wang-Xu elements . Research seminar talk at Humboldt-Universit \"a t zu Berlin, 2021

  36. [36]

    Mallik and N

    G. Mallik and N. Nataraj. A nonconforming finite element approximation for the von Karman equations. ESAIM: M2AN , 50(2):433--454, March 2016

  37. [37]

    M. T. Nakao. A numerical approach to the proof of existence of solutions for elliptic problems. Japan J. Appl. Math. , 5(2):313--332, June 1988

  38. [38]

    M. T. Nakao, K. Hashimoto, and Y. Watanabe. A numerical method to verify the invertibility of linear elliptic operators with applications to nonlinear problems. Computing , 75(1):1--14, July 2005

  39. [39]

    M. T. Nakao, M. Plum, and Y. Watanabe. Numerical Verification Methods and Computer-Assisted Proofs for Partial Differential Equations , volume 53 of Springer Series in Computational Mathematics . Springer Singapore, Singapore, 2019

  40. [40]

    M. Plum. Computer-assisted proofs for semilinear elliptic boundary value problems. Japan J. Indust. Appl. Math. , 26(2-3):419--442, October 2009

  41. [41]

    a t zu Berlin, Mathematisch-Naturwissenschaftliche Fakult \

    S. L. Puttkammer. Direct Guaranteed Lower Eigenvalue Bounds with Quasi-Optimal Adaptive Mesh-Refinement . PhD thesis, Humboldt-Universit \"a t zu Berlin, Mathematisch-Naturwissenschaftliche Fakult \"a t, 2024

  42. [42]

    Quinn and A

    F. Quinn and A. Sard. Hausdorff conullity of critical images of Fredholm maps. Amer. J. Math. , 94(4):1101, October 1972

  43. [43]

    S. Smale. An infinite dimensional version of Sard 's theorem. Amer. J. Math. , 87(4):861--866, 1965

  44. [44]

    Saut and R

    J. Saut and R. Temam. Generic properties of nonlinear boundary value problems. Commun. Partial Differ. Equ. , 4(3):293--319, January 1979

  45. [45]

    D. B. Szyld. The many proofs of an identity on the norm of oblique projections. Numer. Algorithms , 42(3-4):309--323, November 2006

  46. [46]

    G. Talenti. Best constant in Sobolev inequality. Ann. Mat. Pura Appl. , 110(1):353--372, December 1976

  47. [47]

    R. Temam. Navier-- Stokes Equations and Nonlinear Functional Analysis : Second Edition . Society for Industrial and Applied Mathematics , January 1995

  48. [48]

    Takayasu, X

    A. Takayasu, X. Liu, and S. Oishi. Verified computations to semilinear elliptic boundary value problems on arbitrary polygonal domains. NOLTA , 4(1):34--61, 2013

  49. [49]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. II---overconsistency and classical nonconforming elements. SIAM J. Numer. Anal. , 57(1):266--292, January 2019

  50. [50]

    Watanabe, M

    Y. Watanabe, M. T. Nakao, and K. Nagatou. On the compactness of a nonlinear operator related to stream function-vorticity formulation for the Navier-Stokes equations. JSIAM Letters , 9(0):77--80, 2017

  51. [51]

    E. Zeidler. Nonlinear Functional Analysis and Its Applications . 1: Fixed-point Theorems . Springer, New York Berlin Heidelberg, corr. 2. print edition, 1992

  52. [52]

    E. Zeidler. Nonlinear Functional Analysis and Its Applications . 4: Applications to Mathematical Physics . Springer, New York Berlin Heidelberg, corr. 2. print edition, 1997