pith. sign in

arxiv: 2510.05597 · v2 · submitted 2025-10-07 · 🧮 math.NA · cs.NA

Optimal L²-error estimates for the nonsymmetric Nitsche method in two dimensions

Gang Chen , Chaoran Liu , Yangwen Zhang This is my paper

Pith reviewed 2026-05-18 09:29 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Nitsche methodnonsymmetric formulationL2 error estimatesfinite element methodoptimal convergenceDirichlet boundary conditionsPoisson equationtwo dimensions
0
0 comments X

The pith

The stabilized nonsymmetric Nitsche method achieves optimal L2 error bounds of order k+1 for kth-order finite elements on convex polygonal domains in two dimensions.

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

The paper proves that the stabilized nonsymmetric Nitsche method for weakly enforcing Dirichlet conditions on Poisson's equation delivers the full optimal L2 convergence rate. Earlier analyses had predicted a possible half-order loss, yet this work shows the error satisfies a bound of C h^{k+1} times the W^{k+1,∞} norm of the solution. This matters for a sympathetic reader because it resolves a gap between theory and observed numerical behavior, supporting reliable use of the method in practice. The argument compares the Nitsche solution to an auxiliary conforming finite element solution with projected boundary data and relies on a two-layer boundary-strip lifting, an exact boundary identity on the one-dimensional boundary mesh, and localized residual estimates. The authors also supply a revised proof of the needed auxiliary W^{1,∞} estimate using L^∞ stability of the boundary L2 projection and a weak discrete maximum principle.

Core claim

For conforming kth-order finite elements on quasi-uniform triangulations of convex polygonal domains in two dimensions, the stabilized nonsymmetric Nitsche approximation satisfies ||u - u_h||_{L^2(Ω)} ≤ C h^{k+1} ||u||_{W^{k+1,∞}(Ω)}. The proof proceeds by comparing the Nitsche solution with an auxiliary conforming finite element solution that strongly imposes projected boundary data. It combines three ingredients: a two-layer boundary-strip lifting, an exact boundary identity on the one-dimensional boundary mesh, and localized residual estimates. The authors isolate the auxiliary W^{1,∞} estimate and give a revised proof based on the L^∞-stability of the boundary L2-projection together with

What carries the argument

The two-layer boundary-strip lifting together with the exact boundary identity on the one-dimensional boundary mesh, which together enable comparison to an auxiliary conforming solution.

Load-bearing premise

The domain must be a convex polygon and the mesh must be quasi-uniform in two dimensions for the boundary lifting and identity steps to apply.

What would settle it

Numerical computation of L2 errors on a sequence of successively refined quasi-uniform triangulations of a square with a smooth exact solution, checking whether the observed rate equals k+1, would settle the claim.

Figures

Figures reproduced from arXiv: 2510.05597 by Chaoran Liu, Gang Chen, Yangwen Zhang.

Figure 1
Figure 1. Figure 1: Four macro elements Lemma 12. If K is big enough, for each uh ∈ Vh, there exists wh ∈ H1 (Ω) ∩ P1(Th) such that ⟨uh, ∇wh · n⟩Γ ≥ C∥h − 1 2 uh∥ 2 0,Γ , (6.25) 21 [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The first level of grid for numerical experiment [PITH_FULL_IMAGE:figures/full_fig_p032_2.png] view at source ↗
read the original abstract

Nitsche's method is a standard device for weakly imposing Dirichlet boundary conditions, but for the stabilized nonsymmetric formulation the available $L^2$-error analysis for Poisson's equation still predicts a half-order loss, whereas numerical evidence indicates optimal convergence. We prove that, for conforming $k$th-order finite elements on quasi-uniform triangulations of convex polygonal domains in two dimensions, the stabilized nonsymmetric Nitsche approximation satisfies \[ \|{u-u_h}\|_{L^2(\Omega)} \le C h^{k+1}\|{u}\|_{W^{k+1,\infty}(\Omega)}. \] The proof compares the Nitsche solution with an auxiliary conforming finite element solution with strongly imposed projected boundary data and combines three ingredients: a two-layer boundary-strip lifting, an exact boundary identity on the one-dimensional boundary mesh, and localized residual estimates. In addition, we isolate the auxiliary $W^{1,\infty}$ estimate needed in the argument and provide a revised proof based on the $L^\infty$-stability of the boundary $L^2$-projection together with a weak discrete maximum principle for discrete harmonic functions. The analysis is intrinsically two-dimensional and clarifies why the stronger assumption $u\in W^{k+1,\infty}(\Omega)$ enters the proof.

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 manuscript proves that the stabilized nonsymmetric Nitsche method for the Poisson equation on convex polygonal domains in two dimensions achieves optimal L² convergence: for conforming kth-order finite elements on quasi-uniform triangulations, ||u - u_h||_{L²(Ω)} ≤ C h^{k+1} ||u||_{W^{k+1,∞}(Ω)}. The argument proceeds by comparing the Nitsche solution to an auxiliary conforming Galerkin solution with strongly imposed L²-projected boundary data, then controlling the difference via a two-layer boundary-strip lifting, an exact identity on the one-dimensional boundary mesh, and localized residual estimates. A revised proof of the auxiliary W^{1,∞} bound is supplied using L^∞-stability of the boundary projection together with a weak discrete maximum principle for discrete harmonic functions.

Significance. If the result holds, the paper resolves the longstanding discrepancy between the half-order loss predicted by existing analyses and the optimal rates observed numerically for the nonsymmetric Nitsche formulation. The explicit decomposition into three 2D-specific ingredients and the revised auxiliary estimate constitute a clear technical advance for the analysis of weakly imposed boundary conditions. The work is grounded in standard tools of finite-element theory and makes the dependence on convexity, quasi-uniformity, and the W^{k+1,∞} regularity assumption transparent.

minor comments (3)
  1. The abstract and introduction should explicitly reference the section containing the statement of the main theorem (presumably Theorem 1.1 or 2.1) so that readers can locate the precise hypotheses on the mesh and domain without searching.
  2. Notation for the boundary projection operator and the two-layer strip should be introduced once in a preliminary section and then used consistently; occasional re-definition in later sections risks confusion.
  3. A short remark clarifying why the exact boundary identity (the 1D mesh identity) does not hold in three dimensions would strengthen the discussion of the method’s scope without lengthening the paper.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive evaluation of our manuscript. The referee summary accurately captures the main result and the structure of the proof. We appreciate the recognition of the technical advance in resolving the discrepancy between analysis and numerics for the nonsymmetric Nitsche method.

Circularity Check

0 steps flagged

No significant circularity: derivation proceeds by direct comparison and standard estimates

full rationale

The paper proves the optimal L² error bound by comparing the nonsymmetric Nitsche solution to an auxiliary conforming finite-element solution with strongly imposed L²-projected boundary data. The argument relies on a two-layer boundary-strip lifting, an exact identity on the one-dimensional boundary mesh, and localized residual estimates, all constructed from standard finite-element techniques and Sobolev-space arguments without any fitted parameters or self-referential definitions. The auxiliary W^{1,∞} bound is obtained independently via L^∞-stability of the boundary projection together with a weak discrete maximum principle for discrete harmonic functions. No equation or step reduces the claimed h^{k+1} rate to the input data by construction, and the analysis is explicitly presented as self-contained within the stated two-dimensional hypotheses.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard Sobolev-space properties, finite-element approximation theory, and the validity of the three 2D-specific technical ingredients (boundary-strip lifting, exact boundary identity, localized residuals) whose proofs are not supplied in the abstract.

axioms (2)
  • standard math Standard approximation properties of conforming finite-element spaces and Sobolev embeddings hold on convex polygonal domains in 2D
    Invoked to control the auxiliary solution and boundary terms.
  • domain assumption The two-layer boundary-strip lifting and exact one-dimensional boundary identity exist and satisfy the required estimates
    These are the load-bearing 2D constructions named in the abstract.

pith-pipeline@v0.9.0 · 5766 in / 1397 out tokens · 59354 ms · 2026-05-18T09:29:30.935854+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

12 extracted references · 12 canonical work pages

  1. [1]

    Mfem: A modular finite element methods library.Computers & Mathematics with Applications, 81:42– 74, 2021

    Robert Anderson, Julian Andrej, Andrew Barker, Jamie Bramwell, Jean-Sylvain Camier, Jakub Cerveny, Veselin Dobrev, Yohann Dudouit, Aaron Fisher, Tzanio Kolev, et al. Mfem: A modular finite element methods library.Computers & Mathematics with Applications, 81:42– 74, 2021

    work page 2021

  2. [2]

    Arnold, Franco Brezzi, Bernardo Cockburn, and L

    Douglas N. Arnold, Franco Brezzi, Bernardo Cockburn, and L. Donatella Marini. Unified analysis of discontinuous galerkin methods for elliptic problems.SIAM Journal on Numerical Analysis, 39(5):1749–1779, 2002

    work page 2002

  3. [3]

    Brenner and L

    Susanne C. Brenner and L. Ridgway Scott.The mathematical theory of finite element methods, volume 15 ofTexts in Applied Mathematics. Springer, New York, third edition, 2008

    work page 2008

  4. [4]

    A penalty-free nonsymmetric Nitsche-type method for the weak imposition of boundary conditions.SIAM J

    Erik Burman. A penalty-free nonsymmetric Nitsche-type method for the weak imposition of boundary conditions.SIAM J. Numer. Anal., 50(4):1959–1981, 2012

    work page 1959

  5. [5]

    Springer-Verlag, Berlin, 1988

    Monique Dauge.Elliptic boundary value problems on corner domains, volume 1341 ofLecture Notes in Mathematics. Springer-Verlag, Berlin, 1988. Smoothness and asymptotics of solutions

    work page 1988

  6. [6]

    Springer, Cham, [2021]©2021

    Alexandre Ern and Jean-Luc Guermond.Finite elements I—Approximation and interpolation, volume 72 ofTexts in Applied Mathematics. Springer, Cham, [2021]©2021

    work page 2021

  7. [7]

    On weakly imposed boundary conditions for second order problems

    Jouni Freund and Rolf Stenberg. On weakly imposed boundary conditions for second order problems. 1995

    work page 1995

  8. [8]

    Nitsche’s method for general boundary conditions.Math

    Mika Juntunen and Rolf Stenberg. Nitsche’s method for general boundary conditions.Math. Comp., 78(267):1353–1374, 2009

    work page 2009

  9. [9]

    J. Nitsche. ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind.Abh. Math. Sem. Univ. Ham- burg, 36:9–15, 1971

    work page 1971

  10. [10]

    Society for Industrial and Applied Mathemat- ics (SIAM), Philadelphia, PA, 2008

    B´ eatrice Rivi` ere.Discontinuous Galerkin methods for solving elliptic and parabolic equations, volume 35 ofFrontiers in Applied Mathematics. Society for Industrial and Applied Mathemat- ics (SIAM), Philadelphia, PA, 2008. Theory and implementation

    work page 2008

  11. [11]

    Ridgway Scott and Shangyou Zhang

    L. Ridgway Scott and Shangyou Zhang. Finite element interpolation of nonsmooth functions satisfying boundary conditions.Math. Comp., 54(190):483–493, 1990

    work page 1990

  12. [12]

    Wloka.Partial differential equations

    J. Wloka.Partial differential equations. Cambridge University Press, Cambridge, 1987. Trans- lated from the German by C. B. Thomas and M. J. Thomas. 36

    work page 1987