pith. sign in

arxiv: 2408.15731 · v2 · submitted 2024-08-28 · 🧮 math.NA · cs.NA

Finite element discretization of the steady, generalized Navier-Stokes equations for small shear stress exponents

Alex Kaltenbach , Julius Je{\ss}berger This is my paper

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

classification 🧮 math.NA cs.NA
keywords finite element methodgeneralized Navier-Stokes equationsdivergence reconstruction operatorsa priori error estimatesshear stress exponentincompressible flownumerical analysis
0
0 comments X

The pith

Finite element discretization for generalized Navier-Stokes equations works for shear stress exponents p > 2d/(d+2).

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

The paper proposes a finite element discretization for the steady incompressible generalized Navier-Stokes equations that handles fully inhomogeneous data. Divergence reconstruction operators are used to obtain a formulation valid for all shear stress exponents above 2d/(d+2). A priori error estimates are derived for the velocity vector field and the kinematic pressure. Numerical experiments confirm quasi-optimality of the velocity error estimate in all cases and of the pressure estimate when p is at most 2.

Core claim

By the method of divergence reconstruction operators, a finite element formulation is obtained for the steady, incompressible, fully inhomogeneous generalized Navier-Stokes equations that remains valid for all shear stress exponents p > 2d/(d+2). The Dirichlet boundary condition is imposed strongly using any discretization of the boundary data that converges at a sufficient rate. A priori error estimates for the velocity vector field and kinematic pressure are derived, and these estimates are quasi-optimal for the velocity in all admissible cases.

What carries the argument

Divergence reconstruction operators that restore consistency and supply the approximation properties needed for well-posedness and error analysis at small p.

If this is right

  • The discretization applies to the entire range of shear stress exponents p > 2d/(d+2).
  • A priori error estimates hold for both the velocity vector field and the kinematic pressure.
  • The velocity error estimate is quasi-optimal for every admissible p.
  • The pressure error estimate is quasi-optimal whenever p ≤ 2.
  • Strong imposition of the Dirichlet boundary condition is compatible with the method.

Where Pith is reading between the lines

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

  • The reconstruction technique may simplify code reuse by allowing standard finite element spaces to be employed without additional stabilization terms.
  • Error estimates of this type could be used to design a posteriori indicators for adaptive mesh refinement in non-Newtonian flow simulations.
  • The same operator construction might apply to related nonlinear elliptic systems that suffer from low integrability at small exponents.

Load-bearing premise

Divergence reconstruction operators exist with the approximation properties needed to make the discrete formulation well-posed and to derive the error estimates for p > 2d/(d+2).

What would settle it

Numerical computation on a sequence of successively refined meshes showing that the discrete velocity or pressure error fails to converge at the predicted rate for some fixed p just above 2d/(d+2) would falsify the a priori estimates.

read the original abstract

A finite element (FE) discretization for the steady, incompressible, fully inhomogeneous, generalized Navier-Stokes equations is proposed. By the method of divergence reconstruction operators, the formulation is valid for all shear stress exponents $p > \tfrac{2d}{d+2}$. The Dirichlet boundary condition is imposed strongly, using any discretization of the boundary data which converges at a sufficient rate. $\textit{A priori}$ error estimates for the velocity vector field and kinematic pressure are derived and numerical experiments are conducted. These confirm the quasi-optimality of the $\textit{a priori}$ error estimate for the velocity vector field. The $\textit{a priori}$ error estimates for the kinematic pressure are quasi-optimal if $p \leq 2$.

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 / 2 minor

Summary. The paper proposes a finite element discretization for the steady, incompressible, fully inhomogeneous generalized Navier-Stokes equations. By the method of divergence reconstruction operators, the formulation is valid for all shear stress exponents p > 2d/(d+2). The Dirichlet boundary condition is imposed strongly using any discretization of the boundary data which converges at a sufficient rate. A priori error estimates for the velocity vector field and kinematic pressure are derived, and numerical experiments confirm the quasi-optimality of the a priori error estimate for the velocity vector field (with pressure estimates quasi-optimal if p ≤ 2).

Significance. If the claimed construction of divergence reconstruction operators with the required approximation properties holds, the work extends finite element methods for generalized Navier-Stokes to a wider and practically relevant range of small p, which is important for certain non-Newtonian fluid models. The derivation of a priori estimates together with numerical confirmation of quasi-optimality for the velocity constitutes a concrete strength; the explicit handling of the fully inhomogeneous case and strong boundary imposition are additional positive features.

minor comments (2)
  1. [Abstract] Abstract: the statement that experiments 'confirm the quasi-optimality' would benefit from a brief indication of the specific p-values and mesh families tested, to allow immediate assessment of the covered regime.
  2. The paper could add a short remark on how the divergence reconstruction operators are constructed in practice (e.g., reference to a specific section or algorithm) to aid readers interested in implementation.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of the manuscript, the recognition of its contributions to extending finite element methods for the generalized Navier-Stokes equations to the range p > 2d/(d+2), and the recommendation of minor revision. No major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper proposes a new finite-element discretization for the generalized Navier-Stokes equations that relies on constructing divergence reconstruction operators to achieve well-posedness and a priori error estimates for p > 2d/(d+2). The central steps are the definition of the discrete formulation, the proof that the operators deliver the required approximation properties, and the subsequent derivation of the error bounds. These steps are presented as original constructions and proofs rather than reductions to fitted parameters, self-citations, or renamed prior results. No load-bearing premise collapses to an input by definition or to a self-citation chain; the work is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no specific free parameters, axioms, or invented entities can be identified from the provided text.

pith-pipeline@v0.9.0 · 5651 in / 1025 out tokens · 23646 ms · 2026-05-23T22:23:02.792111+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

35 extracted references · 35 canonical work pages

  1. [1]

    P. R. Amestoy, I. S. Duff, J. Koster, and J.-Y. L’Excellent , A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM Journal on Matrix Analysis and Applications, 23 (2001), pp. 15–41

    work page 2001

  2. [2]

    D. N. Arnold, F. Brezzi, and M. Fortin , A stable finite element for the Stokes equations, Calcolo, 21 (1984), pp. 337–344 (1985), https://doi.org/10.1007/BF02576171. FE DISCRETIZATION OF STEADY, GENERALIZED NAVIER–STOKES EQUATIONS21

    work page doi:10.1007/bf02576171 1984

  3. [3]

    Bartels, Numerical methods for nonlinear partial differential equations, vol

    S. Bartels, Numerical methods for nonlinear partial differential equations, vol. 47 of Springer Ser. Comput. Math., Cham: Springer, 2015, https://doi.org/10.1007/978-3-319-13797-1_1

    work page doi:10.1007/978-3-319-13797-1_1 2015

  4. [4]

    Beirão da Veiga, P

    H. Beirão da Veiga, P. Kaplický, and M. Růžička , Boundary regularity of shear–thickening flows, J. Math. Fluid Mech., 13 (2011), pp. 387–404

    work page 2011

  5. [5]

    Belenki, L

    L. Belenki, L. Berselli, L. Diening, and M. Růžička , On the finite element approximation of p -stokes systems, SIAM Journal on Numerical Analysis, 50 (2012), pp. 373–397, https: //doi.org/10.2307/41582741

    work page doi:10.2307/41582741 2012

  6. [6]

    Bernardi and G

    C. Bernardi and G. Raugel , Analysis of some finite elements for the Stokes problem, Math. Comp., 44 (1985), pp. 71–79, https://doi.org/10.2307/2007793

    work page doi:10.2307/2007793 1985

  7. [7]

    Berselli, L

    L. Berselli, L. Diening, and M. Ruzicka , Existence of strong solutions for incompressible fluids with shear dependent viscosities, Journal of Mathematical Fluid Mechanics, 12 (2010), pp. 101–132, https://doi.org/10.1007/s00021-008-0277-y

    work page doi:10.1007/s00021-008-0277-y 2010

  8. [8]

    R. Bird, R. Armstrong, and O. Hassager , Dynamics of polymeric liquids, 2nd Edition, vol. Volume 1, Fluid Mechanics, Wiley, 2nd ed., 1987

    work page 1987

  9. [9]

    Analyse Numérique, 9 (1975), pp

    Clément, Ph., Approximation by finite element functions using local regularization, R.A.I.R.O. Analyse Numérique, 9 (1975), pp. 77–84, https://doi.org/10.1051/m2an/197509R200771

    work page doi:10.1051/m2an/197509r200771 1975

  10. [10]

    Crouzeix and P.-A

    M. Crouzeix and P.-A. Raviart , Conforming and nonconforming finite element methods for solving the stationary Stokes equations. I, Rev. Française Automat. Informat. Recherche Opérationnelle Sér. Rouge, 7 (1973), pp. 33–75

    work page 1973

  11. [11]

    Diening, C

    L. Diening, C. Kreuzer, and E. Süli , Finite element approximation of steady flows of incompressible fluids with implicit power-law-like rheology, SIAM Journal on Numerical Analysis, 51 (2012), https://doi.org/10.1137/120873133

    work page doi:10.1137/120873133 2012

  12. [12]

    Diening, J

    L. Diening, J. Málek, and M. Steinhauer , On Lipschitz truncations of Sobolev functions (with variable exponent) and their selected applications, ESAIM: Control, Opt. Calc. Var., 14 (2008), pp. 211–232

    work page 2008

  13. [13]

    Diening and M

    L. Diening and M. Růžička , Interpolation operators in Orlicz-Sobolev spaces, Numer. Math., 107 (2007), pp. 107–129, https://doi.org/10.1007/s00211-007-0079-9

    work page doi:10.1007/s00211-007-0079-9 2007

  14. [14]

    Diening and M

    L. Diening and M. Růžička , Non-newtonian fluids and function spaces, in Nonlinear Analysis, Function Spaces and Applications, Institute of Mathematics of the Academy of Sciences of the Czech Republic, 2007, pp. 95–143, http://eudml.org/doc/221529

    work page 2007

  15. [15]

    Diening, J

    L. Diening, J. Storn, and T. Tscherpel , On the sobolev andLp-stability of theL2-projection, SIAM Journal on Numerical Analysis, 59 (2021), pp. 2571–2607, https://doi.org/10.1137/ 20M1358013

    work page 2021

  16. [16]

    D. T. W. L. Douglas, J. jr. , The stability in lq of the l2-projection into finite element function spaces., Numerische Mathematik, 23 (1974/75), pp. 193–198, http://eudml.org/doc/132299

    work page 1974

  17. [17]

    Dziuk, Theorie und Numerik partieller Differentialgleichungen, Walter de Gruyter GmbH & Co

    G. Dziuk, Theorie und Numerik partieller Differentialgleichungen, Walter de Gruyter GmbH & Co. KG, Berlin, 2010, https://doi.org/10.1515/9783110214819

    work page doi:10.1515/9783110214819 2010

  18. [18]

    D. A. H. et al., Firedrake User Manual, Imperial College London and University of Oxford and Baylor University and University of Washington, first edition ed., 5 2023, https: //doi.org/10.25561/104839

    work page doi:10.25561/104839 2023

  19. [19]

    Ern and J.-L

    A. Ern and J.-L. Guermond , Theory and practice of finite elements, vol. 159 of Ap- plied Mathematical Sciences, Springer-Verlag, New York, 2004, https://doi.org/10.1007/ 978-1-4757-4355-5

    work page 2004

  20. [20]

    Ern and J

    A. Ern and J. L. Guermond , Finite Elements I: Approximation and Interpolation, no. 1 in Texts in Applied Mathematics, Springer International Publishing, 2021, https://doi.org/10. 1007/978-3-030-56341-7

    work page 2021

  21. [21]

    Ern and J.-L

    A. Ern and J.-L. Guermond , Finite elements II—Galerkin approximation, elliptic and mixed PDEs, vol. 73 of Texts in Applied Mathematics, Springer, Cham, [2021]©2021, https://doi.org/10.1007/978-3-030-56923-5

    work page doi:10.1007/978-3-030-56923-5 2021

  22. [22]

    F arrell, P

    P. F arrell, P. A. Gazca Orozco, and E. Süli , Finite element approximation and precondi- tioning for anisothermal flow of implicitly-constituted non-Newtonian fluids, Math. Comp., 91 (2022), pp. 659–697, https://doi.org/10.1090/mcom/3703

    work page doi:10.1090/mcom/3703 2022

  23. [23]

    Giaquinta and G

    M. Giaquinta and G. Modica , Remarks on the regularity of the minimizers of certain degenerate functionals, Manuscripta Math., 57 (1986), pp. 55–99, https://doi.org/10.1007/ BF01172492

    work page 1986

  24. [24]

    Girault and P.-A

    V. Girault and P.-A. Raviart , Finite element methods for Navier-Stokes equations, vol. 5 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 1986, https: //doi.org/10.1007/978-3-642-61623-5. Theory and algorithms

    work page doi:10.1007/978-3-642-61623-5 1986

  25. [25]

    Giusti, Direct methods in the calculus of variations, World Scientific Publishing Co., Inc., River Edge, NJ, 2003, https://doi.org/10.1142/9789812795557

    E. Giusti, Direct methods in the calculus of variations, World Scientific Publishing Co., Inc., River Edge, NJ, 2003, https://doi.org/10.1142/9789812795557

    work page doi:10.1142/9789812795557 2003

  26. [26]

    Jeßberger and A

    J. Jeßberger and A. Kaltenbach , Finite element discretization of the steady, generalized Navier-Stokes equations with inhomogeneous Dirichlet boundary conditions, SIAM J. Numer. 22 J. JEßBERGER AND A. KALTENBACH Anal., 62 (2024), pp. 1660–1686, https://doi.org/10.1137/23M1607398

    work page doi:10.1137/23m1607398 2024

  27. [27]

    Jeßberger and M

    J. Jeßberger and M. Růžička, Existence of weak solutions for inhomogeneous generalized navier–stokes equations, Nonlinear Analysis, 212 (2021), p. 112538, https://doi.org/https: //doi.org/10.1016/j.na.2021.112538

    work page doi:10.1016/j.na.2021.112538 2021

  28. [28]

    V. John, A. Linke, C. Merdon, M. Neilan, and L. G. Rebholz , On the divergence constraint in mixed finite element methods for incompressible flows, SIAM Rev., 59 (2017), pp. 492–544, https://doi.org/10.1137/15M1047696

    work page doi:10.1137/15m1047696 2017

  29. [29]

    Kaltenbach and M

    A. Kaltenbach and M. Ruzicka , A local discontinuous galerkin approximation for the p- navier–stokes system, part ii: Convergence rates for the velocity, SIAMJournalonNumerical Analysis, 61 (2023), pp. 1641–1663, https://doi.org/10.1137/22M1514751

    work page doi:10.1137/22m1514751 2023

  30. [30]

    Kaplický, J

    P. Kaplický, J. Málek, and J. Stará , C1,α-regularity of weak solutions to a class of nonlinear fluids in two dimensions - stationary Dirichlet problem, Zap. Nauchn. Sem. Pt. Odel. Mat. Inst., 259 (1999), pp. 89–121

    work page 1999

  31. [31]

    P. L. Lederer, A. Linke, C. Merdon, and J. Schöberl , Divergence-free reconstruction operators for pressure-robust Stokes discretizations with continuous pressure finite elements, SIAM J. Numer. Anal., 55 (2017), pp. 1291–1314, https://doi.org/10.1137/16M1089964

    work page doi:10.1137/16m1089964 2017

  32. [32]

    L. R. Scott and S. Zhang , Finite element interpolation of nonsmooth functions satisfying boundary conditions, Mathematics of Computation, 54 (1990), pp. 483–493, http://www. jstor.org/stable/2008497 (accessed 2024-01-23)

    work page arXiv 1990

  33. [33]

    Taylor and P

    C. Taylor and P. Hood , A numerical solution of the Navier-Stokes equations using the finite element technique, Internat. J. Comput. & Fluids, 1 (1973), pp. 73–100, https: //doi.org/10.1016/0045-7930(73)90027-3

    work page doi:10.1016/0045-7930(73)90027-3 1973

  34. [34]

    Temam, Navier-Stokes equations, vol

    R. Temam, Navier-Stokes equations, vol. 2 of Studies in Mathematics and its Applications, North-Holland Publishing Co., Amsterdam, third ed., 1984. Theory and numerical analysis, With an appendix by F. Thomasset

    work page 1984

  35. [35]

    T. Tscherpel, Finite element approximation for the unsteady flow of implicitly constituted incompressible fluids, PhD thesis, University of Oxford, 2018, https://ora.ox.ac.uk/objects/ uuid:01b4901c-9705-4087-80c1-4d656d160aed

    work page 2018