pith. sign in

arxiv: 2508.06235 · v4 · pith:UPZ6TWBBnew · submitted 2025-08-08 · 🧮 math.NA · cs.NA

Fully discrete error analysis of finite element discretizations of time-dependent Stokes equations in a stream-function formulation

Dmitriy Leykekhman , Boris Vexler , Jakob Wagner This is my paper

Pith reviewed 2026-05-21 23:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords finite element methodstream-function formulationtime-dependent Stokes equationserror estimatesdiscontinuous GalerkinGalerkin orthogonalitybest approximation
0
0 comments X

The pith

Error estimates for fully discrete solutions of time-dependent Stokes equations in stream-function formulation achieve best approximation rates.

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

The paper establishes best approximation type error estimates for the fully discrete Galerkin solutions of the time-dependent Stokes problem using the stream-function formulation. Time discretization employs the discontinuous Galerkin method of arbitrary degree, while space discretization is handled through a general framework that covers methods satisfying specific Galerkin orthogonality conditions. This approach requires only the natural regularity given by the domain and data, without additional assumptions. The results apply to conformal C1 and C0 interior penalty methods and can be used for optimal control problems.

Core claim

We establish best approximation type error estimates for the fully discrete Galerkin solutions of the time-dependent Stokes problem using the stream-function formulation. For the time discretization we use the discontinuous Galerkin method of arbitrary degree, whereas we present the space discretization in a general framework. This makes our result applicable for a wide variety of space discretization methods, provided some Galerkin orthogonality conditions are satisfied. As an example, conformal C^1 and C^0 interior penalty methods are covered by our analysis. The results do not require any additional regularity assumptions beyond the natural regularity given by the domain and data and can

What carries the argument

General framework for space discretization requiring Galerkin orthogonality conditions, paired with discontinuous Galerkin time discretization of arbitrary degree, to obtain best approximation error estimates.

Load-bearing premise

The space discretization must satisfy the Galerkin orthogonality conditions required by the general framework.

What would settle it

Numerical experiments on a covered method such as the C0 interior penalty scheme where observed convergence rates fall short of the predicted best approximation rates would falsify the error estimates.

read the original abstract

In this paper we establish best approximation type error estimates for the fully discrete Galerkin solutions of the time-dependent Stokes problem using the stream-function formulation. For the time discretization we use the discontinuous Galerkin method of arbitrary degree, whereas we present the space discretization in a general framework. This makes our result applicable for a wide variety of space discretization methods, provided some Galerkin orthogonality conditions are satisfied. As an example, conformal $C^1$ and $C^0$ interior penalty methods are covered by our analysis. The results do not require any additional regularity assumptions beyond the natural regularity given by the domain and data and can be used for optimal control problems.

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 establishes best-approximation error estimates for the fully discrete Galerkin solutions of the time-dependent Stokes problem in stream-function formulation. Time discretization employs discontinuous Galerkin methods of arbitrary degree, while spatial discretization is treated in a general framework that requires only certain Galerkin orthogonality conditions. The analysis covers conforming C¹ elements and C⁰ interior-penalty methods and holds under the natural regularity induced by the data and domain, with potential application to optimal control problems.

Significance. If the central estimates hold, the work supplies a flexible, unified error analysis for multiple spatial discretizations of the stream-function form of the time-dependent Stokes equations. The reduction to Galerkin orthogonality plus standard approximation theory, together with explicit verification for both conforming and nonconforming methods, is a clear strength; the lack of extra regularity assumptions further enhances practical utility, especially for optimal-control contexts.

major comments (1)
  1. [§3] §3 (general framework): the fully discrete error bound is derived by combining the spatial orthogonality with DG time-stepping; the manuscript should explicitly verify that the time-discretization consistency terms do not introduce additional non-orthogonal contributions when the spatial operator is only weakly consistent (as in the C⁰-IP case).
minor comments (2)
  1. [Introduction] The introduction would benefit from a short paragraph recalling the precise natural regularity class (e.g., the precise Sobolev indices for velocity and pressure) that is assumed throughout the estimates.
  2. [§4.2] In the verification for the C⁰ interior-penalty method, the dependence of the hidden constants on the penalty parameter should be stated explicitly so that the reader can see how the error bound scales with this parameter.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the thorough review and the positive assessment of our manuscript. We are pleased that the work is viewed as providing a flexible unified error analysis. We address the major comment below and will incorporate the suggested clarification in the revised version.

read point-by-point responses

  1. Referee: [§3] §3 (general framework): the fully discrete error bound is derived by combining the spatial orthogonality with DG time-stepping; the manuscript should explicitly verify that the time-discretization consistency terms do not introduce additional non-orthogonal contributions when the spatial operator is only weakly consistent (as in the C⁰-IP case).

    Authors: We appreciate this observation. In the general framework of Section 3, the error analysis proceeds by testing the error equation with the discrete solution and exploiting the Galerkin orthogonality in space at each time step. The discontinuous Galerkin time discretization introduces consistency terms that are estimated using the properties of the DG method, which do not rely on strong orthogonality but rather on the weak form. For methods with weak consistency such as the C⁰ interior penalty method, the spatial consistency error is already bounded by the best-approximation term in our estimates. The time-discretization terms remain orthogonal in the sense that they do not couple with the spatial weak consistency in a way that introduces new non-orthogonal contributions, as the spatial operator is evaluated consistently within the time-stepping scheme. To address the referee's request for explicit verification, we will add a dedicated remark in §3 clarifying this point with reference to the C⁰-IP example, ensuring the argument is transparent for weakly consistent spatial discretizations. revision: yes

Circularity Check

0 steps flagged

Error estimates derive from standard approximation theory and verified Galerkin orthogonality

full rationale

The manuscript establishes best-approximation fully discrete error bounds by reducing the analysis to a general framework that requires only Galerkin orthogonality of the spatial method plus standard approximation properties of the finite element spaces. It then applies discontinuous Galerkin time discretization of arbitrary degree. The paper explicitly checks that both conforming C¹ elements and C⁰ interior-penalty methods satisfy the needed orthogonality identities for the fourth-order stream-function formulation, using only the natural regularity of the data and domain. No step equates a derived quantity to a fitted parameter by construction, renames a known result, or relies on a load-bearing self-citation whose validity is assumed rather than independently verified. The derivation is therefore self-contained against external benchmarks of approximation theory and does not reduce to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard Sobolev-space theory for the Stokes problem and the assumption that chosen space methods obey Galerkin orthogonality; no free parameters or new entities are introduced.

axioms (2)
  • domain assumption The domain and data supply only the natural regularity of the time-dependent Stokes problem.
    Explicitly invoked to avoid additional smoothness requirements.
  • domain assumption Space discretizations satisfy Galerkin orthogonality conditions.
    Required for the general framework to cover C1 and C0 interior penalty methods.

pith-pipeline@v0.9.0 · 5645 in / 1266 out tokens · 37069 ms · 2026-05-21T23:58:07.722113+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

33 extracted references · 33 canonical work pages

  1. [1]

    A DAK, D

    D. A DAK, D. M ORA , S. N ATARAJAN , AND A. S ILGADO , A virtual element discretization for the time dependent Navier-Stokes equations in stream-function formulation, ESAIM Math. Model. Numer. Anal., 55 (2021), pp. 2535–2566

    work page 2021

  2. [2]

    A HMED , G

    N. A HMED , G. R. B ARRENECHEA , E. B URMAN , J. G UZM ´AN, A. L INKE , AND C. M ERDON , A pressure-robust discretization of Oseen’s equation using stabilization in the vorticity equation, SIAM J. Numer. Anal., 59 (2021), pp. 2746–2774

    work page 2021

  3. [3]

    A HMED , A

    N. A HMED , A. L INKE , AND C. M ERDON , Towards pressure-robust mixed methods for the incompressible Navier-Stokes equations, Comput. Methods Appl. Math., 18 (2018), pp. 353–372

    work page 2018

  4. [4]

    A INSWORTH AND C

    M. A INSWORTH AND C. PARKER , Mass conserving mixed ℎ𝑝-FEM approximations to Stokes flow. Part I: Uniform stability, SIAM J. Numer. Anal., 59 (2021), pp. 1218–1244

    work page 2021

  5. [5]

    H. W. A LT, Linear Functional Analysis, Universitext, Springer-Verlag London, Ltd., London, 2016. An application-oriented introduc- tion, Translated from the German edition by Robert N¨urnberg

    work page 2016

  6. [6]

    P. R. S. A NTUNES , On the buckling eigenvalue problem, J. Phys. A, 44 (2011), pp. 215205, 13

    work page 2011

  7. [7]

    B ACUTA , J

    C. B ACUTA , J. H. B RAMBLE , AND J. E. P ASCIAK , Shift theorems for the biharmonic dirichlet problem, in Recent Progress in Com- putational and Applied PDES, Springer US, pp. 1–26

    work page

  8. [8]

    F. B. B ELGACEM AND S. C. B RENNER , Some nonstandard finite element estimates with applications to 3D Poisson and Signorini problems, Electronic Transactions on Numerical Analysis, 12 (2001), pp. 134–148

    work page 2001

  9. [9]

    B LUM AND R

    H. B LUM AND R. R ANNACHER , On the boundary value problem of the biharmonic operator on domains with angular corners , 2, pp. 556–581. eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/mma.1670020416

    work page doi:10.1002/mma.1670020416

  10. [10]

    S. C. B RENNER AND L.-Y. S UNG, 𝐶0 interior penalty methods for fourth order elliptic boundary value problems on polygonal domains, J. Sci. Comput., 22/23 (2005), pp. 83–118

    work page 2005

  11. [11]

    C IARLET , Basic error estimates for elliptic problems, in Finite Element Methods (Part 1), vol

    P. C IARLET , Basic error estimates for elliptic problems, in Finite Element Methods (Part 1), vol. 2 of Handbook of Numerical Analysis, Elsevier, 1991, pp. 17–351

    work page 1991

  12. [12]

    P. G. C IARLET AND J. L. L IONS , eds., Handbook of numerical analysis. Vol. IX, vol. IX of Handbook of Numerical Analysis, North- Holland, Amsterdam, 2003. Numerical methods for fluids. Part 3

    work page 2003

  13. [13]

    L. C. E VANS, Partial Differential Equations, vol. 19 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, second ed., 2010

    work page 2010

  14. [14]

    R. S. F ALK AND M. N EILAN , Stokes complexes and the construction of stable finite elements with pointwise mass conservation, SIAM J. Numer. Anal., 51 (2013), pp. 1308–1326. 22 DMITRIY LEYKEKHMAN AND BORIS VEXLER and JAKOB W AGNER

    work page 2013

  15. [15]

    G IRAULT AND P

    V. G IRAULT AND P. A. R AVIART, Finite element methods for Navier-Stokes equations , vol. 5 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 1986

    work page 1986

  16. [16]

    G UERMOND AND L

    J. G UERMOND AND L. Q UARTAPELLE , Equivalence of u– p and 𝜁 − 𝜓 formulations of the time-dependent Navier–Stokes equations, International Journal for Numerical Methods in Fluids, 18 (1994), pp. 471–487

    work page 1994

  17. [17]

    G UERMOND AND L

    J.-L. G UERMOND AND L. Q UARTAPELLE , Equivalence of u-𝑝 and 𝜁-𝜓 formulations of the time-dependent Navier-Stokes equations, Internat. J. Numer. Methods Fluids, 18 (1994), pp. 471–487

    work page 1994

  18. [18]

    M. D. G UNZBURGER , Finite element methods for viscous incompressible flows: a guide to theory, practice, and algorithms, Elsevier, 2012

    work page 2012

  19. [19]

    G UZM ´AN AND M

    J. G UZM ´AN AND M. N EILAN , Conforming and divergence-free Stokes elements in three dimensions, IMA J. Numer. Anal., 34 (2014), pp. 1489–1508

    work page 2014

  20. [20]

    Comp., 83 (2014), pp

    , Conforming and divergence-free Stokes elements on general triangular meshes, Math. Comp., 83 (2014), pp. 15–36

    work page 2014

  21. [21]

    J OHN , A

    V. J OHN , A. L INKE , C. M ERDON , M. N EILAN , AND L. G. R EBHOLZ , On the divergence constraint in mixed finite element methods for incompressible flows, SIAM Rev., 59 (2017), pp. 492–544

    work page 2017

  22. [22]

    V. A. K ONDRAT ’EV, Boundary value problems for elliptic equations in domains with conical or angular points , Trudy Moskov. Mat. Obˇsˇc., 16 (1967), pp. 209–292

    work page 1967

  23. [23]

    L EYKEKHMAN AND B

    D. L EYKEKHMAN AND B. V EXLER , Pointwise best approximation results for Galerkin finite element solutions of parabolic problems, SIAM J. Numer. Anal., 54 (2016), pp. 1365–1384

    work page 2016

  24. [24]

    Math., 135 (2017), pp

    , Discrete maximal parabolic regularity for Galerkin finite element methods, Numer. Math., 135 (2017), pp. 923–952

    work page 2017

  25. [25]

    , Discrete maximal parabolic regularity for Galerkin finite element methods, Numerische Mathematik, 135 (2017), pp. 923–952

    work page 2017

  26. [26]

    L INKE , C

    A. L INKE , C. M ERDON , AND M. N EILAN , Pressure-robustness in quasi-optimal a priori estimates for the Stokes problem, Electron. Trans. Numer. Anal., 52 (2020), pp. 281–294

    work page 2020

  27. [27]

    M EIDNER , Adaptive Space-Time Finite Element Methods for Optimization Problems Governed by Nonlinear Parabolic Systems , PhD thesis, Ruprecht-Karls-Universit¨at, Heidelberg, 2008

    D. M EIDNER , Adaptive Space-Time Finite Element Methods for Optimization Problems Governed by Nonlinear Parabolic Systems , PhD thesis, Ruprecht-Karls-Universit¨at, Heidelberg, 2008

    work page 2008

  28. [28]

    M EIDNER AND B

    D. M EIDNER AND B. V EXLER , A priori error estimates for space-time finite element discretization of parabolic optimal control problems. I. Problems without control constraints, SIAM J. Control Optim., 47 (2008), pp. 1150–1177

    work page 2008

  29. [29]

    L. R. S COTT AND M. V OGELIUS , Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials, RAIRO Mod´el. Math. Anal. Num´er., 19 (1985), pp. 111–143

    work page 1985

  30. [30]

    T EMAM , Navier-Stokes equations: theory and numerical analysis, AMS Chelsea Pub

    R. T EMAM , Navier-Stokes equations: theory and numerical analysis, AMS Chelsea Pub

    work page

  31. [31]

    T EMAM , Navier-Stokes equations

    R. T EMAM , Navier-Stokes equations. Theory and numerical analysis , North-Holland Publishing Co., Amsterdam-New York-Oxford,

    work page

  32. [32]

    Studies in Mathematics and its Applications, V ol. 2

    work page

  33. [33]

    W IENERS , A numerical existence proof of nodal lines for the first eigenfunction of the plate equation, Arch

    C. W IENERS , A numerical existence proof of nodal lines for the first eigenfunction of the plate equation, Arch. Math. (Basel), 66 (1996), pp. 420–427

    work page 1996