pith. sign in

arxiv: 2605.04455 · v1 · submitted 2026-05-06 · 🧮 math.NA · cs.NA

Long-time L²&H¹-stability of the Family of DLN Methods for the Two-dimensional Incompressible Navier-Stokes Equations

Isabel Barrio Sanchez , Wenlong Pei , Catalin Trenchea This is my paper

Pith reviewed 2026-05-08 16:51 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords DLN methodsNavier-Stokes equationslong-time stabilityG-stabilityuniform time gridsGrönwall inequalityincompressible flowsnumerical methods for PDEs
0
0 comments X

The pith

DLN methods for the two-dimensional incompressible Navier-Stokes equations possess uniform-in-time L² and H¹ stability under uniform time grids and mild step-size constraints.

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

The paper examines the long-time behavior of the one-leg DLN family of methods, parameterized by θ, when applied to the 2D incompressible Navier-Stokes equations. It derives a new version of the G-stability identity that holds specifically on uniform time grids satisfying a mild restriction on the step size. This identity is then combined with the discrete uniform Grönwall inequality to establish that the numerical solutions remain bounded in both the L² and H¹ norms, with the bounds independent of the length of the time interval and of the initial data. The resulting stability statements reproduce the long-time behavior already known for the continuous problem in two dimensions.

Core claim

Under uniform time grids and mild time-step constraints, the family of one-leg DLN methods satisfies a new G-stability identity. This identity, together with the discrete uniform Grönwall inequality lemma, yields L² and H¹ bounds on the numerical solutions that remain independent of the time-interval length and of the initial conditions, in agreement with the known long-time stability theory for the continuous two-dimensional incompressible Navier-Stokes equations.

What carries the argument

The new G-stability identity for DLN methods on uniform time grids, which supplies the key auxiliary relation needed to close the discrete uniform Grönwall argument.

Load-bearing premise

The new G-stability identity holds only when the time grid is uniform and the time steps satisfy a mild size restriction; without these conditions the identity and the subsequent Grönwall argument may fail to produce time-uniform bounds.

What would settle it

A concrete numerical solution computed on a uniform time grid obeying the mild step-size restriction whose L² or H¹ norm grows without bound as the number of steps increases would falsify the claimed uniform stability.

read the original abstract

In this report, we study the long-time stability of the family of one-leg DLN methods for the two-dimensional incompressible Navier-Stokes equations. The family of DLN methods (with one parameter $\theta$), non-linear energy stable ($G$-stable) and second-order accurate under arbitrary time grids, has been widely applied to the simulations of various fluid models with success. We derive a new version of the $G$-stability identity for the family of DLN methods under uniform time grids and mild time constraints. Then we utilize this crucial auxiliary tool and the discrete uniform Gr\"onwall inequality lemma to prove the uniform-in-time stability of the numerical solutions. Essentially, the bounds are independent of the time interval and the initial conditions, consistent with the theories of the continuous case.

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 investigates the long-time L² and H¹ stability of the family of one-leg DLN methods (parameterized by θ) for the 2D incompressible Navier-Stokes equations. It derives a new G-stability identity valid under uniform time grids and mild time-step constraints, then applies the discrete uniform Grönwall inequality to obtain uniform-in-time bounds on the numerical solutions that depend only on viscosity, forcing, and the fixed time-step size, independent of the time interval length and initial data, consistent with continuous 2D NSE theory.

Significance. If the central claims hold, the work supplies rigorous long-time analysis for a family of second-order, G-stable methods already used in fluid simulations. The tailored G-stability identity on uniform grids enables the uniform Gronwall argument, yielding an absorbing ball whose radius matches the continuous case; this is a concrete strength for numerical analysis of NSE and supports reliable long-time computations.

minor comments (2)
  1. [Abstract] Abstract: the notation 'L^2&$H^1$-stability' is awkward; replace with 'L² and H¹ stability' for readability.
  2. [Introduction / Theorem 3.1] The precise form of the 'mild time constraints' (e.g., the explicit bound on Δt in terms of ν or data) should be stated once in the introduction and repeated in the statement of the main theorem.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive assessment of our work on the long-time L² and H¹ stability of the DLN family for the 2D incompressible Navier-Stokes equations. The referee's summary correctly identifies the key technical contributions: the new G-stability identity under uniform time grids and the subsequent application of the discrete uniform Grönwall inequality to obtain bounds independent of the time interval length.

Circularity Check

0 steps flagged

No significant objection identified

full rationale

The paper derives a new G-stability identity under uniform time grids and mild constraints, then feeds it into the standard discrete uniform Grönwall lemma to obtain time-uniform L2 and H1 bounds. This chain is self-contained: the identity follows from the DLN one-leg scheme definition on uniform steps, the Grönwall tool is external and parameter-free, and the resulting absorbing ball depends only on viscosity, forcing, and fixed Δt. Prior citations to DLN methods exist but are not load-bearing for the stability result; no equation reduces by construction to a fitted input or self-citation encoding the target bound.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the existence of a new G-stability identity valid only for uniform grids plus the applicability of the discrete uniform Gronwall lemma to the resulting energy inequality. No free parameters or invented entities are introduced.

axioms (2)
  • domain assumption The family of DLN methods satisfies a G-stability identity under uniform time steps and mild step-size restrictions.
    Invoked in the derivation step that precedes application of Gronwall.
  • standard math The discrete uniform Gronwall inequality applies directly to the energy estimate obtained from the new identity.
    Standard lemma used to absorb the time-dependent terms and obtain time-uniform bounds.

pith-pipeline@v0.9.0 · 5450 in / 1353 out tokens · 19698 ms · 2026-05-08T16:51:55.583076+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

44 extracted references · 44 canonical work pages

  1. [1]

    Akbas, S

    M. Akbas, S. Kaya, and L. G. Rebholz. On the stability at all times of linearly extrapolated BDF2 timestepping for multiphysics incompressible flow problems.Numer. Methods Partial Differential Equations, 33(4):999–1017, 2017

    work page 2017

  2. [2]

    Akbas Belenli, L

    M. Akbas Belenli, L. G. Rebholz, and F. Tone. A note on the importance of mass conservation in long-time stability of Navier-Stokes simulations using finite elements.Appl. Math. Lett., 45:98–102, 2015

    work page 2015

  3. [3]

    Breckling and S

    S. Breckling and S. Shields. The long-timeL 2 andH 1 stability of linearly extrapolated second- order time-stepping schemes for the 2D incompressible Navier-Stokes equations. Appl. Math. Comput., 342:263–279, 2019

    work page 2019

  4. [4]

    H. Chen, S. Sun, and T. Zhang. Energy stability analysis of some fully discrete numerical schemes for incompressible Navier-Stokes equations on staggered grids. J. Sci. Comput., 75(1):427–456, 2018

    work page 2018

  5. [5]

    Y. Chen, D. Luo, W. Pei, and Y. Xing. Efficient variable time-stepping adaptive DLN algo- rithms for the Allen-Cahn equation. J. Sci. Comput., 104(3):Paper No. 67, 28, 2025

    work page 2025

  6. [6]

    Cheng and C

    K. Cheng and C. Wang. Long time stability of high order multistep numerical schemes for two- dimensional incompressible Navier-Stokes equations. SIAM J. Numer. Anal., 54(5):3123– 3144, 2016

    work page 2016

  7. [7]

    G. G. Dahlquist, W. Liniger, and O. Nevanlinna. Stability of two-step methods for variable integration steps. SIAM J. Numer. Anal., 20(5):1071–1085, 1983

    work page 1983

  8. [8]

    Foias, O

    C. Foias, O. Manley, R. Rosa, and R. Temam. Navier-Stokes equations and turbulence, vol- ume 83 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 2001

    work page 2001

  9. [9]

    Formaggia, A

    L. Formaggia, A. Quarteroni, and A. Veneziani, editors. Cardiovascular mathematics, volume 1 of MS&A. Modeling, Simulation and Applications. Springer-Verlag, Milan, 2009. Modeling and simulation of the circulatory system

    work page 2009

  10. [10]

    G. P. Galdi. An introduction to the mathematical theory of the Navier-Stokes equations. Steady-state problems. Springer Monographs in Mathematics. Springer, New York, second edition, 2011

    work page 2011

  11. [11]

    Gottlieb, F

    S. Gottlieb, F. Tone, C. Wang, X. Wang, and D. Wirosoetisno. Long time stability of a classical efficient scheme for two-dimensional Navier-Stokes equations. SIAM J. Numer. Anal., 50(1):126–150, 2012

    work page 2012

  12. [12]

    Grisvard

    P. Grisvard. Elliptic problems in nonsmooth domains, volume 24 of Monographs and Studies in Mathematics. Pitman (Advanced Publishing Program), Boston, MA, 1985

    work page 1985

  13. [13]

    M. D. Gunzburger. Perspectives in flow control and optimization, volume 5 of Advances in Design and Control. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2003

    work page 2003

  14. [14]

    Heister, M

    T. Heister, M. A. Olshanskii, and L. G. Rebholz. Unconditional long-time stability of a velocity- vorticity method for the 2D Navier-Stokes equations. Numer. Math., 135(1):143–167, 2017

    work page 2017

  15. [15]

    in the large

    O. A. Ladyˇ zenskaja. Solution “in the large” of the boundary value problem for the Navier-Stokes equations in two space variables. Soviet Physics. Dokl., 3:1128–1131, 1958. Translation of Dokl. Akad. Nauk SSSR 123 (1958), 427–429

    work page 1958

  16. [16]

    in the large

    O. A. Ladyˇ zenskaja. Solution “in the large” of the nonstationary boundary value problem for the Navier-Stokes system with two space variables. Comm. Pure Appl. Math., 12:427–433, 1959

    work page 1959

  17. [17]

    O. A. Ladyˇ zenskaja.The mathematical theory of viscous incompressible flow. Mathematics and its Applications, Vol. 2. Gordon and Breach Science Publishers, New York-London-Paris, second English edition, 1969

    work page 1969

  18. [18]

    Layton, W

    W. Layton, W. Pei, Y. Qin, and C. Trenchea. Analysis of the variable step method of Dahlquist, Liniger and Nevanlinna for fluid flow. Numer. Methods Partial Differential Equations, 38(6):1713–1737, 2022

    work page 2022

  19. [19]

    Layton, W

    W. Layton, W. Pei, and C. Trenchea. Refactorization of a variable step, unconditionally stable method of Dahlquist, Liniger and Nevanlinna. Appl. Math. Lett., 125:Paper No. 107789, 7, 2022

    work page 2022

  20. [20]

    Layton, W

    W. Layton, W. Pei, and C. Trenchea. Time step adaptivity in the method of Dahlquist, Liniger and Nevanlinna. Advances in Computational Science and Engineering, 1(3):320–350, 2023

    work page 2023

  21. [21]

    J. Leray. ´Etude de diverses ´ equations int´ egrales non lin´ eaires et de quelques probl` emes que pose l’hydrodynamique. J. Math. Pures Appl. (9), 12:1–82, 1933

    work page 1933

  22. [22]

    W. Pei. The semi-implicit DLN algorithm for the Navier-Stokes equations. Numerical Algorithms, 97(4):1673–1713, 2024

    work page 2024

  23. [23]

    W. Pei. The variable time-stepping DLN-ensemble algorithms for incompressible Navier-Stokes equations. Numerical Algorithms, 2025

    work page 2025

  24. [24]

    Y. Qin, L. Chen, Y. Wang, Y. Li, and J. Li. An adaptive time-stepping DLN decoupled 23 algorithm for the coupled Stokes-Darcy model. Appl. Numer. Math., 188:106–128, 2023

    work page 2023

  25. [25]

    Y. Qin, Y. Hou, W. Pei, and J. Li. A variable time-stepping algorithm for the unsteady Stokes/Darcy model. J. Comput. Appl. Math., 394:Paper No. 113521, 14, 2021

    work page 2021

  26. [26]

    Rebholz and F

    L. Rebholz and F. Tone. Long-timeH 1-stability of BDF2 time stepping for 2D Navier-Stokes equations. Appl. Math. Lett., 141:Paper No. 108624, 8, 2023

    work page 2023

  27. [27]

    I. B. Sanchez, C. Trenchea, and W. Pei. Long-time H1-stability of the Cauchy one-leg Theta- method for the Navier-Stokes equations. arXiv preprint arXiv:2603.27861, 2026

    work page arXiv 2026

  28. [28]

    J. Shen. Long time stability and convergence for fully discrete nonlinear Galerkin methods. Appl. Anal., 38(4):201–229, 1990

    work page 1990

  29. [29]

    Siddiqua and W

    F. Siddiqua and W. Pei. Variable time step method of Dahlquist, Liniger and Nevanlinna (DLN) for a corrected Smagorinsky model. International Journal of Numerical Analysis and Modeling, 21(6):879–909, 2024

    work page 2024

  30. [30]

    R. Temam. Navier-Stokes equations. Theory and numerical analysis. Studies in Mathematics and its Applications. Vol. 2. North-Holland Publishing Co., Amsterdam-New York-Oxford, 1977

    work page 1977

  31. [31]

    R. Temam. Infinite-dimensional dynamical systems in mechanics and physics, volume 68 of Applied Mathematical Sciences. Springer-Verlag, New York, second edition, 1997

    work page 1997

  32. [32]

    F. Tone. On the long-time stability of the Crank-Nicolson scheme for the 2D Navier-Stokes equations. Numer. Methods Partial Differential Equations, 23(5):1235–1248, 2007

    work page 2007

  33. [33]

    Tone and D

    F. Tone and D. Wirosoetisno. On the long-time stability of the implicit Euler scheme for the two-dimensional Navier-Stokes equations. SIAM J. Numer. Anal., 44(1):29–40, 2006

    work page 2006

  34. [34]

    G. K. Vallis. Atmospheric and oceanic fluid dynamics. Cambridge University Press, Cambridge,

    work page

  35. [35]

    Fundamentals and large-scale circulation

    work page

  36. [36]

    H. Wang, X. Wang, and M. Zhang. Unconditionally long-time stable variable-step second- order exponential time-differencing schemes for the incompressible NSE. arXiv preprint arXiv:2602.10268, 2026

    work page arXiv 2026

  37. [37]

    X. Wang. Approximation of stationary statistical properties of dissipative dynamical systems: time discretization. Math. Comp., 79(269):259–280, 2010

    work page 2010

  38. [38]

    X. Wang. An efficient second order in time scheme for approximating long time statistical properties of the two dimensional Navier-Stokes equations. Numer. Math., 121(4):753– 779, 2012

    work page 2012

  39. [39]

    T. T. Warner. Numerical weather and climate prediction. Cambridge University Press, Cam- bridge, 2011

    work page 2011

  40. [40]

    W. M. Washington and C. L. Parkinson. Introduction to three-dimensional climate modeling. G—Reference, Information and Interdisciplinary Subjects Series. University Science Books, 2005

    work page 2005

  41. [41]

    Xie and K

    C. Xie and K. Wang. Uniform finite element error estimates with power-type asymptotic constants for unsteady Navier-Stokes equations. Entropy, 24(7):Paper No. 948, 21, 2022

    work page 2022

  42. [42]

    Zhang and J

    T. Zhang and J. Yuan. Unconditional stability and optimal error estimates of Euler implicit/explicit-SAV scheme for the Navier-Stokes equations.J. Sci. Comput., 90(1):Paper No. 1, 20, 2022

    work page 2022

  43. [43]

    Zheng, Q

    N. Zheng, Q. Guan, W. Pei, and W. Zhao. A divergence constrained mixed finite element method for thermally driven active fluid model. Applied Numerical Mathematics, 227:81– 107, 2026

    work page 2026

  44. [44]

    Zheng, X

    N. Zheng, X. Guo, W. Pei, and W. Zhao. Mixed finite element methods for an incompressible fourth-order active fluid model. Journal of Scientific Computing, 106(3):65, 2026. 24

    work page 2026