pith. sign in

arxiv: 2604.23453 · v1 · submitted 2026-04-25 · 🧮 math.NA · cs.NA

A robust a posteriori error estimator for the Oseen problem

Muhammad Afzal , Naveed Ahmed , Volker John This is my paper

Pith reviewed 2026-05-08 07:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords a posteriori error estimatorOseen problemconvection-dominated regimestabilized finite elementsSUPGrobustnessNavier-Stokes equations
0
0 comments X

The pith

A residual-based a posteriori error estimator remains robust for the convection-dominated Oseen problem discretized by stabilized finite elements.

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

The paper proposes a residual-based a posteriori error estimator for the incompressible Oseen problem in the convection-dominated regime using a SUPG/PSPG/grad-div stabilized finite element discretization. The estimator bounds the global error in the same norm employed in the a priori analysis of the method. Robustness of the estimator with respect to the convection parameter is proved under several hypotheses on the error and interpolation errors, supported by numerical studies. The work also discusses extending the estimator to the steady-state Navier-Stokes equations.

Core claim

The central claim is that the proposed residual-based a posteriori error estimator controls the discretization error for the Oseen problem in the convection-dominated regime. The estimator is constructed for the SUPG/PSPG/grad-div stabilized finite element method and is shown to be robust in the norm of the a priori error analysis, with the proof relying on hypotheses concerning the error and interpolation errors. Numerical experiments confirm the robustness, and the estimator is extended to the steady Navier-Stokes equations.

What carries the argument

The residual-based a posteriori error estimator built from local residuals and interelement jumps, designed to estimate the error in the norm of the a priori analysis and proven robust under hypotheses on approximation errors.

Load-bearing premise

The proof of robustness in the convection-dominated regime rests on several hypotheses concerning the error and interpolation errors.

What would settle it

A sequence of numerical computations with successively larger convection parameters where the ratio of the estimator to the true error becomes unbounded would falsify the robustness claim.

Figures

Figures reproduced from arXiv: 2604.23453 by Muhammad Afzal, Naveed Ahmed, Volker John.

Figure 1
Figure 1. Figure 1: Example 4.1. Errors }pu ´ uh, p ´ phq}spg for different pairs of finite element spaces and different values of the viscosity coefficient. where the individual contributions are defined in Theorem 5 and always the constant C “ 1 was used. Two examples in two dimensions are considered. The first one has a smooth solution without layers and the solution of the second example exhibits boundary layers. Both exa… view at source ↗
Figure 2
Figure 2. Figure 2: Example 4.1. The a posteriori error estimator view at source ↗
Figure 3
Figure 3. Figure 3: Example 4.1. Effectivity indices for different pairs of finite element spaces and different values of the viscosity view at source ↗
Figure 4
Figure 4. Figure 4: Example 4.1. Contributions of the individual parts of view at source ↗
Figure 5
Figure 5. Figure 5: Example 4.1. Comparison of 2 }pu ´ uh, p ´ phq}2 spg (blue line) and of 2}u ´ I V h u} 2 L 2 pΩq (gray line) with terms that appear in the hypotheses from Hypothesis 4 as well as with the sum of the two last terms from the error bound (29), simulations with ν “ 10´6 . (18), since in all cases continuous pressure approximations were used, so that the term rF puh, phq, and with that ηF , contains the small f… view at source ↗
Figure 6
Figure 6. Figure 6: Example 4.2. Adaptively refined meshes for the view at source ↗
Figure 7
Figure 7. Figure 7: Example 4.2. Errors }pu ´ uh, p ´ phq}spg for different pairs of finite element spaces and different values of the viscosity coefficient view at source ↗
Figure 8
Figure 8. Figure 8: Example 4.2. The a posteriori error estimator view at source ↗
Figure 9
Figure 9. Figure 9: Example 4.2. Effectivity indices for different pairs of finite element spaces and different values of the viscosity view at source ↗
Figure 10
Figure 10. Figure 10: Example 4.2. Contributions of the individual parts of view at source ↗
Figure 11
Figure 11. Figure 11: Navier–Stokes problem. Errors }pu ´ uh, p ´ phq}spg,nse for different pairs of finite element spaces and different values of the viscosity coefficient. think that an appropriate norm is }pv, qq}spg,nse “ ˜ ν }∇v} 2 L 2 pΩq ` ν }q} 2 L 2 pΩq ` ÿ KPTh µK }∇ ¨ v} 2 L 2 pKq ` ÿ KPTh δK }∇q} 2 L 2 pKq ¸1{2 . The numerical studies considered the Navier–Stokes problem with the right-hand side and the Dirichlet b… view at source ↗
Figure 12
Figure 12. Figure 12: Navier–Stokes problem. The a posteriori error estimator view at source ↗
Figure 13
Figure 13. Figure 13: Navier–Stokes problem. Effectivity indices for different pairs of finite element spaces and different values of view at source ↗
read the original abstract

A residual-based a posteriori error estimator is proposed for the incompressible Oseen problem in the convection-dominated regime. The SUPG/PSPG/grad-div stabilized finite element method is used as discretization. The error estimator estimates the global error in a norm that is used in the a priori error analysis of the method. Based on several hypotheses concerning the error and interpolation errors, the robustness of the estimator in the convection-dominated regime is proved. Numerical studies support the analytic results. Finally, the extension of the a posteriori error estimator to the steady-state Navier--Stokes equations is discussed.

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 paper proposes a residual-based a posteriori error estimator for the incompressible Oseen problem discretized using the SUPG/PSPG/grad-div stabilized finite element method. The estimator is constructed to bound the global error in the norm employed in the corresponding a priori analysis. Robustness of the estimator in the convection-dominated regime is proved under several hypotheses on the error and interpolation errors. Numerical studies are cited in support of the analysis, and an extension of the estimator to the steady-state Navier-Stokes equations is discussed.

Significance. If the hypotheses hold and the estimator is robust without further restrictions, the work would provide a practical tool for adaptive mesh refinement in convection-dominated incompressible flow simulations, where standard residual estimators often lose robustness. The explicit alignment between the estimator and the a priori error norm is a methodological strength, as is the inclusion of numerical validation and the brief discussion of the Navier-Stokes extension. These elements would enhance the utility of stabilized methods in applications.

major comments (1)
  1. [Robustness proof] Abstract and robustness proof: The central claim that the estimator is robust in the convection-dominated regime is explicitly conditioned on 'several hypotheses concerning the error and interpolation errors' invoked to control the convection term and relate the estimator to the a priori norm. These hypotheses are not stated explicitly, nor is their validity verified or relaxed for the relevant mesh and parameter regimes. This conditionality is load-bearing for the main result and must be addressed by either removing the hypotheses, proving them, or providing concrete verification.
minor comments (2)
  1. [Abstract] The abstract refers to numerical studies supporting the results but provides no quantitative details (e.g., error tables, specific convection parameters, or mesh sizes). Adding a concise summary or pointer to the relevant tables/figures would improve the abstract's informativeness.
  2. [Notation and preliminaries] Ensure consistent notation for all norms, constants, and stabilization parameters upon first use; a short notation table or list would aid readers unfamiliar with the specific stabilized formulation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We appreciate the recognition of the estimator's potential utility for adaptive mesh refinement in convection-dominated incompressible flow simulations. We address the major comment below.

read point-by-point responses

  1. Referee: Abstract and robustness proof: The central claim that the estimator is robust in the convection-dominated regime is explicitly conditioned on 'several hypotheses concerning the error and interpolation errors' invoked to control the convection term and relate the estimator to the a priori norm. These hypotheses are not stated explicitly, nor is their validity verified or relaxed for the relevant mesh and parameter regimes. This conditionality is load-bearing for the main result and must be addressed by either removing the hypotheses, proving them, or providing concrete verification.

    Authors: We agree that the hypotheses are central to the robustness result and that their presentation requires greater transparency. In the revised manuscript we will explicitly list all hypotheses both in the abstract and in a dedicated paragraph immediately preceding the robustness theorem. We will also add a new subsection that discusses their validity, including additional numerical experiments that verify the hypotheses hold in the convection-dominated regime for a range of mesh sizes and parameter values. We have chosen not to remove the hypotheses, as they are essential to the current proof technique, nor to claim a full proof of their validity in complete generality, which lies beyond the scope of this work. The added numerical verification will make the conditionality concrete and the result more usable. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation remains self-contained

full rationale

The paper proposes a residual-based a posteriori error estimator for the SUPG/PSPG/grad-div stabilized discretization of the Oseen problem and proves its robustness in the convection-dominated regime under several hypotheses on the error and interpolation errors. These hypotheses are invoked as assumptions to control terms in the proof and relate the estimator to the a priori error norm, but they do not reduce the central claim to a self-referential definition or a fitted input renamed as a prediction. No load-bearing self-citations, uniqueness theorems imported from the authors' prior work, or ansatzes smuggled via citation are visible in the abstract or context. The estimator is constructed directly from residuals and tied to an existing norm, which is a standard non-circular construction. The derivation chain is therefore independent and does not collapse to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The robustness proof rests on unstated hypotheses about the error and interpolation errors; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption Several hypotheses concerning the error and interpolation errors
    Invoked explicitly to establish robustness of the estimator in the convection-dominated regime.

pith-pipeline@v0.9.0 · 5386 in / 1156 out tokens · 69955 ms · 2026-05-08T07:20:13.562132+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

32 extracted references · 13 canonical work pages

  1. [1]

    On the grad-div stabilization for the steady Oseen and Navier-Stokes equations.Calcolo, 54(1): 471–501, 2017

    Naveed Ahmed. On the grad-div stabilization for the steady Oseen and Navier-Stokes equations.Calcolo, 54(1): 471–501, 2017. ISSN 0008-0624,1126-5434. doi:10.1007/s10092-016-0194-z. URL https://doi.org/10. 1007/s10092-016-0194-z

    work page doi:10.1007/s10092-016-0194-z 2017

  2. [2]

    Tinsley Oden

    Mark Ainsworth and J. Tinsley Oden. A posteriori error estimators for the Stokes and Oseen equations.SIAM J. Numer. Anal., 34(1):228–245, 1997. ISSN 0036-1429. doi:10.1137/S0036142994264092. URL https: //doi.org/10.1137/S0036142994264092

    work page doi:10.1137/s0036142994264092 1997

  3. [3]

    Tinsley Oden.A posteriori error estimation in finite element analysis

    Mark Ainsworth and J. Tinsley Oden.A posteriori error estimation in finite element analysis. Pure and Applied Mathematics (New York). Wiley-Interscience [John Wiley & Sons], New York, 2000. ISBN 0-471-29411-X. 20 A robust a posteriori error estimator for the Oseen problem

    2000

  4. [4]

    Error estimation for low-order adaptive finite element approximations for fluid flow problems.IMA J

    Alejandro Allendes, Francisco Durán, and Richard Rankin. Error estimation for low-order adaptive finite element approximations for fluid flow problems.IMA J. Numer. Anal., 36(4):1715–1747, 2016. ISSN 0272-4979. doi:10.1093/imanum/drv031. URLhttp://dx.doi.org/10.1093/imanum/drv031

    work page doi:10.1093/imanum/drv031 2016

  5. [5]

    Barrios, J

    Tomás P. Barrios, J. Manuel Cascón, and Marí a González. Augmented mixed finite element method for the Oseen problem: a priori and a posteriori error analyses.Comput. Methods Appl. Mech. Engrg., 313:216–238, 2017. ISSN 0045-7825. doi:10.1016/j.cma.2016.09.012. URLhttp://dx.doi.org/10.1016/j.cma.2016.09.012

    work page doi:10.1016/j.cma.2016.09.012 2017

  6. [6]

    Robustness in a posteriori error analysis for FEM flow models.Numer

    Stefano Berrone. Robustness in a posteriori error analysis for FEM flow models.Numer. Math., 91(3):389–422,

  7. [7]

    doi:10.1007/s002110100370

    ISSN 0029-599X. doi:10.1007/s002110100370. URL http://0-dx.doi.org.library.ucc.ie/10. 1007/s002110100370

    work page doi:10.1007/s002110100370

  8. [8]

    Braack and E

    M. Braack and E. Burman. Local projection stabilization for the Oseen problem and its interpretation as a variational multiscale method.SIAM J. Numer. Anal., 43(6):2544–2566 (electronic), 2006. ISSN 0036-1429

    2006

  9. [9]

    Robust error estimates for stabilized finite element approximations of the two dimensional Navier- Stokes’ equations at high Reynolds number.Comput

    Erik Burman. Robust error estimates for stabilized finite element approximations of the two dimensional Navier- Stokes’ equations at high Reynolds number.Comput. Methods Appl. Mech. Engrg., 288:2–23, 2015. ISSN 0045-7825. doi:10.1016/j.cma.2014.11.006. URL http://0-dx.doi.org.library.ucc.ie/10.1016/j. cma.2014.11.006

    work page doi:10.1016/j.cma.2014.11.006 2015

  10. [10]

    Fernández, and Peter Hansbo

    Erik Burman, Miguel A. Fernández, and Peter Hansbo. Continuous interior penalty finite element method for Oseen’s equations.SIAM J. Numer. Anal., 44(3):1248–1274, 2006. ISSN 0036-1429

    2006

  11. [11]

    Ciarlet.The finite element method for elliptic problems

    Philippe G. Ciarlet.The finite element method for elliptic problems. North-Holland Publishing Co., Amsterdam,

  12. [12]

    Studies in Mathematics and its Applications, V ol

    ISBN 0-444-85028-7. Studies in Mathematics and its Applications, V ol. 4

  13. [13]

    Analysis of the grad-div stabilization for the time-dependent Navier-Stokes equations with inf-sup stable finite elements.Adv

    Javier de Frutos, Bosco Garcí a Archilla, V olker John, and Julia Novo. Analysis of the grad-div stabilization for the time-dependent Navier-Stokes equations with inf-sup stable finite elements.Adv. Comput. Math., 44 (1):195–225, 2018. ISSN 1019-7168. doi:10.1007/s10444-017-9540-1. URL https://doi.org/10.1007/ s10444-017-9540-1

    work page doi:10.1007/s10444-017-9540-1 2018

  14. [14]

    Franca and Sérgio L

    Leopoldo P. Franca and Sérgio L. Frey. Stabilized finite element methods. II. The incompressible Navier-Stokes equations.Comput. Methods Appl. Mech. Engrg., 99(2-3):209–233, 1992. ISSN 0045-7825

    1992

  15. [15]

    Franca and Thomas J

    Leopoldo P. Franca and Thomas J. R. Hughes. Two classes of mixed finite element methods.Comput. Methods Appl. Mech. Engrg., 69(1):89–129, 1988. ISSN 0045-7825

    1988

  16. [16]

    On the convergence order of the finite element error in the kinetic energy for high Reynolds number incompressible flows.Comput

    Bosco García-Archilla, V olker John, and Julia Novo. On the convergence order of the finite element error in the kinetic energy for high Reynolds number incompressible flows.Comput. Methods Appl. Mech. Engrg., 385:Paper No. 114032, 54, 2021. ISSN 0045-7825. doi:10.1016/j.cma.2021.114032. URLhttps://doi.org/10.1016/ j.cma.2021.114032

    work page doi:10.1016/j.cma.2021.114032 2021

  17. [17]

    A velocity-pressure streamline diffusion finite element method for the incompressible Navier-Stokes equations.Comput

    Peter Hansbo and Anders Szepessy. A velocity-pressure streamline diffusion finite element method for the incompressible Navier-Stokes equations.Comput. Methods Appl. Mech. Engrg., 84(2):175–192, 1990. ISSN 0045-7825

    1990

  18. [18]

    A numerical study of a posteriori error estimators for convection-diffusion equations.Comput

    V olker John. A numerical study of a posteriori error estimators for convection-diffusion equations.Comput. Methods Appl. Mech. Engrg., 190(5-7):757–781, 2000. ISSN 0045-7825

    2000

  19. [19]

    Melin, O

    V olker John.Finite element methods for incompressible flow problems, volume 51 ofSpringer Series in Computa- tional Mathematics. Springer, Cham, 2016. ISBN 978-3-319-45749-9; 978-3-319-45750-5. doi:10.1007/978-3- 319-45750-5. URLhttp://dx.doi.org/10.1007/978-3-319-45750-5

    work page doi:10.1007/978-3- 2016

  20. [20]

    MooNMD—a program package based on mapped finite element methods

    V olker John and Gunar Matthies. MooNMD—a program package based on mapped finite element methods. Comput. Vis. Sci., 6(2-3):163–169, 2004. ISSN 1432-9360

    2004

  21. [21]

    Error analysis of the SUPG finite element discretization of evolutionary convection- diffusion-reaction equations.SIAM J

    V olker John and Julia Novo. Error analysis of the SUPG finite element discretization of evolutionary convection- diffusion-reaction equations.SIAM J. Numer. Anal., 49(3):1149–1176, 2011. ISSN 0036-1429

    2011

  22. [22]

    A robust SUPG norm a posteriori error estimator for stationary convection- diffusion equations.Comput

    V olker John and Julia Novo. A robust SUPG norm a posteriori error estimator for stationary convection- diffusion equations.Comput. Methods Appl. Mech. Engrg., 255:289–305, 2013. ISSN 0045-7825. doi:10.1016/j.cma.2012.11.019. URL http://0-dx.doi.org.library.ucc.ie/10.1016/j.cma.2012.11. 019. 21 A robust a posteriori error estimator for the Oseen problem

    work page doi:10.1016/j.cma.2012.11.019 2013

  23. [23]

    Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?Comput

    V olker John, Petr Knobloch, and Julia Novo. Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?Comput. Vis. Sci., 19(5-6):47–63, 2018. ISSN 1432-9360. doi:10.1007/s00791-018-0290-5. URLhttps://doi.org/10.1007/s00791-018-0290-5

    work page doi:10.1007/s00791-018-0290-5 2018

  24. [24]

    A robust a posteriori error estimator for divergence-conforming discontinuous Galerkin methods for the Oseen equation.SIAM J

    Arbaz Khan and Guido Kanschat. A robust a posteriori error estimator for divergence-conforming discontinuous Galerkin methods for the Oseen equation.SIAM J. Numer. Anal., 58(1):492–518, 2020. ISSN 0036-1429,1095-

    2020

  25. [25]

    URLhttps://doi.org/10.1137/18M1169072

    doi:10.1137/18M1169072. URLhttps://doi.org/10.1137/18M1169072

    work page doi:10.1137/18m1169072

  26. [26]

    Springer-Verlag, Berlin, second edition,

    Hans-Görg Roos, Martin Stynes, and Lutz Tobiska.Robust numerical methods for singularly perturbed differential equations, volume 24 ofSpringer Series in Computational Mathematics. Springer-Verlag, Berlin, second edition,

  27. [27]

    ISBN 978-3-540-34466-7

  28. [28]

    A modified streamline diffusion method for solving the stationary Navier-Stokes equation.Numer

    Lutz Tobiska and Gert Lube. A modified streamline diffusion method for solving the stationary Navier-Stokes equation.Numer. Math., 59(1):13–29, 1991. ISSN 0029-599X

    1991

  29. [29]

    Analysis of a streamline diffusion finite element method for the Stokes and Navier-Stokes equations.SIAM J

    Lutz Tobiska and Rüdiger Verfürth. Analysis of a streamline diffusion finite element method for the Stokes and Navier-Stokes equations.SIAM J. Numer. Anal., 33(1):107–127, 1996. ISSN 0036-1429

    1996

  30. [30]

    Verfürth

    R. Verfürth. A posteriori error estimation and adaptive mesh-refinement techniques.J. Comput. Appl. Math., 50 (1-3):67–83, 1994. ISSN 0377-0427

    1994

  31. [31]

    Verfürth

    R. Verfürth. A posteriori error estimators for convection-diffusion equations.Numer. Math., 80(4):641– 663, 1998. ISSN 0029-599X,0945-3245. doi:10.1007/s002110050381. URL https://doi.org/10.1007/ s002110050381

    work page doi:10.1007/s002110050381 1998

  32. [32]

    Numerical Mathematics and Scientific Computation

    Rüdiger Verfürth.A posteriori error estimation techniques for finite element methods. Numerical Mathematics and Scientific Computation. Oxford University Press, Oxford, 2013. ISBN 978-0-19-967942-3. 22

    2013