An Aubin-Nitsche Lemma for a positivity-preserving finite element method

Gabriel R. Barrenechea; Theophile Chaumont-Frelet

arxiv: 2606.21560 · v1 · pith:RZOVSP5Rnew · submitted 2026-06-19 · 🧮 math.NA · cs.NA

An Aubin-Nitsche Lemma for a positivity-preserving finite element method

Gabriel R. Barrenechea , Theophile Chaumont-Frelet This is my paper

Pith reviewed 2026-06-26 13:37 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Aubin-Nitsche lemmapositivity-preserving finite elementL2 error estimatelinearised adjoint problemelliptic problemsnonlinear discretization

0 comments

The pith

A linearised adjoint problem yields an Aubin-Nitsche lemma for positivity-preserving finite element discretizations of elliptic problems.

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

The paper proves an Aubin-Nitsche lemma that delivers optimal-order L2 error estimates for a nonlinear positivity-preserving finite element method on elliptic problems. Standard duality arguments fail because the discretization is nonlinear, so the authors first construct a linearised adjoint problem. They link this adjoint to the original method by choosing appropriate weights, then invoke a regularity result on the adjoint to close the argument and obtain the L2 bound.

Core claim

The central claim is that the linearised adjoint problem, connected to the positivity-preserving discretization through a suitable selection of weights and supported by a regularity result, permits the proof of an optimal-order error estimate in the L2-norm of the error.

What carries the argument

The linearised adjoint problem linked to the nonlinear method by weight selection, which enables the duality argument for the L2 estimate.

If this is right

Optimal L2 error estimates hold for the positivity-preserving discretization of elliptic problems.
The Aubin-Nitsche technique extends to this class of nonlinear finite element methods.
Error control in the L2 norm achieves the same order as in the energy norm for these schemes.

Where Pith is reading between the lines

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

The weight-selection technique for linking the adjoint may extend to other nonlinear positivity-preserving methods.
The result could support improved a posteriori error indicators that use the L2 norm for problems with positivity constraints.
Testing the lemma on obstacle-type elliptic problems would check whether the regularity assumption holds in practice.

Load-bearing premise

A regularity result for the linearised adjoint problem exists that is strong enough to close the duality argument after the weights are chosen.

What would settle it

A concrete numerical example on an elliptic problem where the computed L2 convergence rate of the positivity-preserving method falls below the optimal order predicted by the lemma.

read the original abstract

In this work we prove an Aubin-Nitsche Lemma for a positivity-preserving discretisation of an elliptic problem. Due to the nonlinearity of the discretisation, the result requires as a first step the proposal of a linearised adjoint problem that can be linked to the method of choice by appropriately selecting weights. This linearised adjoint problem, together with a regularity result, allow the proof of an optimal-order error estimate in the $L^2$-norm of the error.

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.

Desk Editor's Note private letter to a colleague

The paper gives a workable construction for an Aubin-Nitsche lemma on a nonlinear positivity-preserving scheme via a weighted linearised adjoint, but the regularity step for that adjoint is the part that still needs explicit checking.

read the letter

The main takeaway is that they have built a linearised adjoint problem whose weights are picked so the discrete relation matches the original nonlinear positivity-preserving method, then used that to run a duality argument for an optimal L2 error bound.

The weight-selection step is the actual new piece. Most Aubin-Nitsche proofs work directly with the linear adjoint; here the scheme is nonlinear, so the linearisation plus the specific weights let them recover the necessary identity without leaving the positivity framework. That is a concrete technical fix for a class of methods that otherwise lack L2 estimates.

The argument is presented as a direct derivation from the linearised problem and a regularity result, with no obvious circularity or reduction to fitted quantities. The abstract is clear about the strategy.

The soft spot is exactly the regularity claim. The proof needs an H2 bound on the linearised adjoint to close the duality argument. Standard elliptic regularity theorems require sufficient smoothness on the coefficients and usually a convex domain. The weights are solution-dependent, so the resulting operator has variable coefficients whose smoothness is not automatic. The abstract simply says the linearised problem together with a regularity result works; it does not show whether the weights preserve the hypotheses of the cited theorem or whether a new regularity proof is supplied. If the weights drop the coefficients below the required class, the H2 estimate fails and the L2 rate does not follow. That is the only load-bearing assumption that is not yet visible from the given material.

This is narrow but useful work for people who already use positivity-preserving finite elements for elliptic problems and need L2 control. It is the sort of targeted technical result that deserves a serious referee so the regularity application can be verified in the full text. I would send it to review.

Referee Report

2 major / 0 minor

Summary. The manuscript claims to prove an Aubin-Nitsche lemma for a positivity-preserving finite element discretization of an elliptic problem. Because the discretization is nonlinear, the argument first constructs a linearised adjoint problem linked to the scheme by a suitable choice of weights; this linearised problem, together with a regularity result, is then used to obtain an optimal-order L2 error estimate via a duality argument.

Significance. If the central argument is complete and the regularity result applies to the weighted linearised operator, the result would extend duality techniques to a class of nonlinear positivity-preserving schemes, which is useful for L2-error analysis in applications where standard linear theory does not directly apply.

major comments (2)

[Abstract / proof outline] The abstract states that the linearised adjoint 'together with a regularity result' closes the proof, but supplies neither the explicit construction of the weights nor a verification that the resulting operator satisfies the hypotheses (C^1 coefficients, uniform ellipticity, convex domain) of the invoked regularity theorem. This verification is load-bearing for the H^2 estimate required by the Aubin-Nitsche step.
[Abstract / proof outline] No derivation steps are given showing that the chosen weights produce the required discrete adjoint identity that links the linearised problem back to the nonlinear positivity-preserving method. Without this identity, the duality argument does not connect to the discretization error.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. The points raised concern the level of detail in the abstract and proof outline; we address them point by point below and will revise the manuscript to improve clarity.

read point-by-point responses

Referee: [Abstract / proof outline] The abstract states that the linearised adjoint 'together with a regularity result' closes the proof, but supplies neither the explicit construction of the weights nor a verification that the resulting operator satisfies the hypotheses (C^1 coefficients, uniform ellipticity, convex domain) of the invoked regularity theorem. This verification is load-bearing for the H^2 estimate required by the Aubin-Nitsche step.

Authors: The weights are constructed explicitly in Section 3.2 by setting them equal to the reciprocal of the positive part of the discrete solution; this choice is used to obtain the linearised operator. The verification that the resulting coefficients are C^1 and the operator remains uniformly elliptic (under the standing assumption of a convex domain) appears in the proof of Lemma 4.3. We agree that these elements are not referenced in the abstract and will revise the abstract to include a short statement on the weight construction and the applicability of the regularity result. revision: yes
Referee: [Abstract / proof outline] No derivation steps are given showing that the chosen weights produce the required discrete adjoint identity that links the linearised problem back to the nonlinear positivity-preserving method. Without this identity, the duality argument does not connect to the discretization error.

Authors: The discrete adjoint identity is obtained in the proof of Theorem 3.1 by substituting the chosen weights into the bilinear form of the linearised problem and verifying that the resulting expression equals the difference between the nonlinear residual and its linearisation. This identity is what allows the duality argument to bound the L2 error by the consistency term. We acknowledge that the intermediate algebraic steps could be expanded for readability and will add them in the revised version. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation uses external regularity result on constructed adjoint

full rationale

The paper proposes a linearised adjoint problem linked to the positivity-preserving scheme by weight selection, then invokes a regularity result to close the Aubin-Nitsche duality argument for the L2 estimate. No quoted step shows the L2 result reducing to a fitted quantity, self-definition, or load-bearing self-citation chain; the regularity is treated as an external input whose applicability is a separate correctness question rather than a circularity reduction. The central claim therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on (i) the existence of a suitable linearised adjoint whose weights recover the discrete method and (ii) an unstated regularity result for that adjoint. No free parameters or invented entities are mentioned.

axioms (1)

domain assumption A regularity result exists for the linearised adjoint problem that is sufficient to obtain the optimal L2 rate.
Invoked explicitly in the abstract as the second ingredient after the linearised adjoint is constructed.

pith-pipeline@v0.9.1-grok · 5614 in / 1181 out tokens · 15722 ms · 2026-06-26T13:37:16.111331+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

25 extracted references · 2 canonical work pages

[1]

Goal-adaptive discretization of fluid-structure interaction , author =
[2]

Proximal

Keith, Brendan and Surowiec, Thomas M , date-added =. Proximal. Foundations of Computational Mathematics , pages =
[3]

and Pryer, Tristan , date-added =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , date-added =. A nodally bound-preserving finite element method for reaction--convection--diffusion equations , volume =. Mathematical Models and Methods in Applied Sciences , number =
[4]

The discrete maximum principle for linear simplicial finite element approximations of a reaction--diffusion problem , volume =

Brandts, Jan H and Korotov, Sergey and K. The discrete maximum principle for linear simplicial finite element approximations of a reaction--diffusion problem , volume =. Linear Algebra and its Applications , number =
[5]

and Pryer, Tristan , doi =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , doi =. A nodally bound-preserving finite element method for reaction--convection--diffusion equations , url =. Math. Models Methods Appl. Sci. , mrclass =. 2024 , bdsk-url-1 =

2024
[6]

and Pryer, Tristan , journal =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , journal =. A nodally bound-preserving finite element method for time-dependent convection--diffusion equations , year =
[7]

and Burman, Erik and Karakatsani, Fotini , fjournal =

Barrenechea, Gabriel R. and Burman, Erik and Karakatsani, Fotini , fjournal =. Edge-based nonlinear diffusion for finite element approximations of convection-diffusion equations and its relation to algebraic flux-correction schemes , volume =. Numer. Math. , number =
[8]

and Georgoulis, Emmanuil H

Barrenechea, Gabriel R. and Georgoulis, Emmanuil H. and Pryer, Tristan and Veeser, Andreas , doi =. A nodally bound-preserving finite element method , url =. IMA J. Numer. Anal. , mrclass =. 2024 , bdsk-url-1 =

2024
[9]

and John, Volker and Knobloch, Petr , doi =

Barrenechea, Gabriel R. and John, Volker and Knobloch, Petr , doi =. Finite element methods respecting the discrete maximum principle for convection-diffusion equations , url =. SIAM Rev. , mrclass =. 2024 , bdsk-url-1 =

2024
[10]

Brenner, S. C. and Scott, L. R. , edition =. The Mathematical Theory of Finite Element Methods , volume =
[11]

A nodal based a posteriori error estimator preserving the positivity for a singularly perturbed reaction-diffusion problem , volume =

Burman, Erik and Guzm. A nodal based a posteriori error estimator preserving the positivity for a singularly perturbed reaction-diffusion problem , volume =. 2016 , bdsk-url-1 =. doi:10.1007/s10915-016-0168-y , journal =

work page doi:10.1007/s10915-016-0168-y 2016
[12]

Stabilized

Burman, Erik and Ern, Alexandre , coden =. Stabilized. Math. Comp. , mrclass =. 2005 , bdsk-url-1 =. doi:10.1090/S0025-5718-05-01761-8 , fjournal =

work page doi:10.1090/s0025-5718-05-01761-8 2005
[13]

, publisher =

Renardy, Michael and Rogers, Robert C. , publisher =. An introduction to partial differential equations , volume =
[14]

and Raviart, Pierre-Arnaud , fjournal =

Ciarlet, Philippe G. and Raviart, Pierre-Arnaud , fjournal =. Maximum principle and uniform convergence for the finite element method , volume =. Comput. Methods Appl. Mech. Engrg. , mrclass =
[15]

Finite Elements I , year =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite Elements I , year =
[16]

Finite Elements II , year =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite Elements II , year =
[17]

Finite elements III: first-order and time-dependent PDEs , volume =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite elements III: first-order and time-dependent PDEs , volume =
[18]

, edition =

Evans, Lawrence C. , edition =. Partial differential equations , volume =
[19]

A note on why enforcing discrete maximum principles by a simple a posteriori cutoff is a good idea , volume =

Kreuzer, Christian , doi =. A note on why enforcing discrete maximum principles by a simple a posteriori cutoff is a good idea , volume =. Numer. Methods Partial Differential Equations , mrclass =. 2014 , bdsk-url-1 =

2014
[20]

Functional analysis,

Brezis, Haim , isbn =. Functional analysis,
[21]

An introduction to variational inequalities and their applications , year =

Kinderlehrer, David and Stampacchia, Guido , publisher =. An introduction to variational inequalities and their applications , year =
[22]

A monotone finite element scheme for convection-diffusion equations , url =

Xu, Jinchao and Zikatanov, Ludmil , doi =. A monotone finite element scheme for convection-diffusion equations , url =. Math. Comp. , mrclass =. 1999 , bdsk-url-1 =

1999
[23]

and John, Volker and Knobloch, Petr , publisher =

Barrenechea, Gabriel R. and John, Volker and Knobloch, Petr , publisher =. Monotone Discretizations for Elliptic Second Order Partial Differential Equations , year =
[24]

2023 , doi =

Kuzmin, Dmitri and Hajduk, Hennes , title =. 2023 , doi =

2023
[25]

and Wahlbin, Lars B

Johnson, Claes and Schatz, Alfred H. and Wahlbin, Lars B. , doi =. Crosswind smear and pointwise errors in streamline diffusion finite element methods , url =. Math. Comp. , OPTnumber =

[1] [1]

Goal-adaptive discretization of fluid-structure interaction , author =

[2] [2]

Proximal

Keith, Brendan and Surowiec, Thomas M , date-added =. Proximal. Foundations of Computational Mathematics , pages =

[3] [3]

and Pryer, Tristan , date-added =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , date-added =. A nodally bound-preserving finite element method for reaction--convection--diffusion equations , volume =. Mathematical Models and Methods in Applied Sciences , number =

[4] [4]

The discrete maximum principle for linear simplicial finite element approximations of a reaction--diffusion problem , volume =

Brandts, Jan H and Korotov, Sergey and K. The discrete maximum principle for linear simplicial finite element approximations of a reaction--diffusion problem , volume =. Linear Algebra and its Applications , number =

[5] [5]

and Pryer, Tristan , doi =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , doi =. A nodally bound-preserving finite element method for reaction--convection--diffusion equations , url =. Math. Models Methods Appl. Sci. , mrclass =. 2024 , bdsk-url-1 =

2024

[6] [6]

and Pryer, Tristan , journal =

Amiri, Abdolreza and Barrenechea, Gabriel R. and Pryer, Tristan , journal =. A nodally bound-preserving finite element method for time-dependent convection--diffusion equations , year =

[7] [7]

and Burman, Erik and Karakatsani, Fotini , fjournal =

Barrenechea, Gabriel R. and Burman, Erik and Karakatsani, Fotini , fjournal =. Edge-based nonlinear diffusion for finite element approximations of convection-diffusion equations and its relation to algebraic flux-correction schemes , volume =. Numer. Math. , number =

[8] [8]

and Georgoulis, Emmanuil H

Barrenechea, Gabriel R. and Georgoulis, Emmanuil H. and Pryer, Tristan and Veeser, Andreas , doi =. A nodally bound-preserving finite element method , url =. IMA J. Numer. Anal. , mrclass =. 2024 , bdsk-url-1 =

2024

[9] [9]

and John, Volker and Knobloch, Petr , doi =

Barrenechea, Gabriel R. and John, Volker and Knobloch, Petr , doi =. Finite element methods respecting the discrete maximum principle for convection-diffusion equations , url =. SIAM Rev. , mrclass =. 2024 , bdsk-url-1 =

2024

[10] [10]

Brenner, S. C. and Scott, L. R. , edition =. The Mathematical Theory of Finite Element Methods , volume =

[11] [11]

A nodal based a posteriori error estimator preserving the positivity for a singularly perturbed reaction-diffusion problem , volume =

Burman, Erik and Guzm. A nodal based a posteriori error estimator preserving the positivity for a singularly perturbed reaction-diffusion problem , volume =. 2016 , bdsk-url-1 =. doi:10.1007/s10915-016-0168-y , journal =

work page doi:10.1007/s10915-016-0168-y 2016

[12] [12]

Stabilized

Burman, Erik and Ern, Alexandre , coden =. Stabilized. Math. Comp. , mrclass =. 2005 , bdsk-url-1 =. doi:10.1090/S0025-5718-05-01761-8 , fjournal =

work page doi:10.1090/s0025-5718-05-01761-8 2005

[13] [13]

, publisher =

Renardy, Michael and Rogers, Robert C. , publisher =. An introduction to partial differential equations , volume =

[14] [14]

and Raviart, Pierre-Arnaud , fjournal =

Ciarlet, Philippe G. and Raviart, Pierre-Arnaud , fjournal =. Maximum principle and uniform convergence for the finite element method , volume =. Comput. Methods Appl. Mech. Engrg. , mrclass =

[15] [15]

Finite Elements I , year =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite Elements I , year =

[16] [16]

Finite Elements II , year =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite Elements II , year =

[17] [17]

Finite elements III: first-order and time-dependent PDEs , volume =

Ern, Alexandre and Guermond, Jean-Luc , publisher =. Finite elements III: first-order and time-dependent PDEs , volume =

[18] [18]

, edition =

Evans, Lawrence C. , edition =. Partial differential equations , volume =

[19] [19]

A note on why enforcing discrete maximum principles by a simple a posteriori cutoff is a good idea , volume =

Kreuzer, Christian , doi =. A note on why enforcing discrete maximum principles by a simple a posteriori cutoff is a good idea , volume =. Numer. Methods Partial Differential Equations , mrclass =. 2014 , bdsk-url-1 =

2014

[20] [20]

Functional analysis,

Brezis, Haim , isbn =. Functional analysis,

[21] [21]

An introduction to variational inequalities and their applications , year =

Kinderlehrer, David and Stampacchia, Guido , publisher =. An introduction to variational inequalities and their applications , year =

[22] [22]

A monotone finite element scheme for convection-diffusion equations , url =

Xu, Jinchao and Zikatanov, Ludmil , doi =. A monotone finite element scheme for convection-diffusion equations , url =. Math. Comp. , mrclass =. 1999 , bdsk-url-1 =

1999

[23] [23]

and John, Volker and Knobloch, Petr , publisher =

Barrenechea, Gabriel R. and John, Volker and Knobloch, Petr , publisher =. Monotone Discretizations for Elliptic Second Order Partial Differential Equations , year =

[24] [24]

2023 , doi =

Kuzmin, Dmitri and Hajduk, Hennes , title =. 2023 , doi =

2023

[25] [25]

and Wahlbin, Lars B

Johnson, Claes and Schatz, Alfred H. and Wahlbin, Lars B. , doi =. Crosswind smear and pointwise errors in streamline diffusion finite element methods , url =. Math. Comp. , OPTnumber =