pith. sign in

arxiv: 1906.11778 · v2 · pith:ZT577XNSnew · submitted 2019-06-27 · 🧮 math.NA · cs.NA· math.AP

Convergence rates for the numerical approximation of the 2D stochastic Navier-Stokes equations

Pith reviewed 2026-05-25 14:27 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.AP
keywords stochastic Navier-Stokes equationsfinite element methodconvergence ratesmultiplicative noisepressure decompositioncylindrical Wiener processtwo-dimensional periodic domain
0
0 comments X

The pith

A finite-element space-time scheme for the 2D stochastic Navier-Stokes equations with linear-growth multiplicative noise converges linearly in space and at order nearly 1/2 in time, in probability.

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

The paper proves convergence rates for a finite-element discretization of the two-dimensional stochastic Navier-Stokes equations on the torus, driven by nonlinear multiplicative noise of linear growth that is generated by a cylindrical Wiener process. Using a stochastic pressure decomposition to control the pressure term, the authors obtain an error bound linear in the spatial mesh size and of order almost one half in the time step, measured in probability in the norm that combines the supremum over time of the L2 norm with the integral over time of the H1 seminorm. This improves the temporal rate from the nearly one-quarter order obtained in earlier work. A reader would care because these sharper rates directly affect how finely one must resolve space and time to obtain reliable statistics from simulations of randomly forced incompressible flows.

Core claim

The finite-element based space-time approximation converges at a linear rate in space and at a rate of almost one half in time, with respect to convergence in probability, where the error is measured in the L^∞_t L²_x ∩ L²_t W^{1,2}_x-norm. The proof relies on a stochastic pressure decomposition for the pressure function under the linear-growth multiplicative noise driven by a cylindrical Wiener process.

What carries the argument

The stochastic pressure decomposition, which decomposes the pressure to enable error analysis under the multiplicative stochastic forcing driven by a cylindrical Wiener process.

If this is right

  • The spatial discretization error is bounded by a constant times the mesh size h.
  • The temporal discretization error is bounded by a constant times the time-step size raised to a power approaching 1/2.
  • The stated rates hold in probability for the combined L^∞_t L²_x and L²_t W^{1,2}_x norm under periodic boundary conditions.
  • The temporal rate improves on the nearly 1/4 order obtained by earlier discretizations of the same equations.

Where Pith is reading between the lines

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

  • The pressure decomposition may extend to other stochastic PDEs whose noise satisfies comparable linear-growth bounds.
  • Adaptive time-stepping schemes could exploit the near-square-root temporal rate to reduce computational cost for long-time statistics.
  • Similar pressure-control arguments might be tested in three space dimensions once a suitable decomposition is available.

Load-bearing premise

The analysis relies on a stochastic pressure decomposition for the pressure function under the given linear-growth multiplicative noise driven by a cylindrical Wiener process.

What would settle it

Numerical experiments that halve the time step repeatedly while holding the spatial mesh fixed, using the same linear-growth multiplicative noise, and that record a temporal convergence order substantially below one half would contradict the claimed rate.

read the original abstract

We study stochastic Navier-Stokes equations in two dimensions with respect to periodic boundary conditions. The equations are perturbed by a nonlinear multiplicative stochastic forcing with linear growth (in the velocity) driven by a cylindrical Wiener process. We establish convergence rates for a finite-element based space-time approximation with respect to convergence in probability (where the error is measure in the $L^\infty_tL^2_x\cap L^2_tW^{1,2}_x$-norm). Our main result provides linear convergence in space and convergence of order (almost) 1/2 in time. This improves earlier results from [E. Carelli, A. Prohl: Rates of convergence for discretizations of the stochastic incompressible Navier-Stokes equations. SIAM J. Numer. Anal. 50(5), 2467-2496. (2012)] where the convergence rate in time is only (almost) 1/4. Our approach is based on a careful analysis of the pressure function using a stochastic pressure decomposition.

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

2 major / 1 minor

Summary. The manuscript analyzes convergence rates for a finite-element space-time discretization of the 2D stochastic Navier-Stokes equations on the periodic torus, driven by nonlinear multiplicative noise with linear growth in the velocity and forced by a cylindrical Wiener process. The central claim is convergence in probability of order 1 in space and order (almost) 1/2 in time, measured in the L^∞_t L²_x ∩ L²_t W^{1,2}_x norm; this improves the temporal rate obtained in Carelli-Prohl (2012) from (almost) 1/4 by means of a stochastic pressure decomposition.

Significance. If the pressure decomposition supplies the asserted estimates and the linear-growth terms are absorbed without order reduction in the stochastic integrals or remainders, the result would constitute a concrete advance in the numerical analysis of stochastic incompressible flows. The 2-D periodic setting and cylindrical noise are standard, so the improvement is localized to the pressure treatment and could be useful for designing higher-order time-stepping schemes.

major comments (2)
  1. [Abstract / pressure decomposition section] Abstract and the section introducing the stochastic pressure decomposition: the claim that this decomposition lifts the temporal rate from (almost) 1/4 to (almost) 1/2 rests on the assertion that the resulting pressure estimates allow the stochastic integrals and nonlinear remainders to retain the higher order under linear-growth multiplicative noise. The manuscript must supply the explicit decomposition (including any Itô corrections) and the corresponding a-priori bounds; without these, it is impossible to confirm that the linear-growth term does not force a reversion to the lower rate obtained in Carelli-Prohl (2012).
  2. [Main theorem / error analysis] The error analysis (likely the main theorem and its proof): the convergence-in-probability statement is stated for the combined space-time error, but the proof sketch must isolate the temporal contribution arising from the pressure term and show that the almost-1/2 rate survives after applying the Burkholder-Davis-Gundy inequality and Gronwall-type arguments to the linear-growth noise. If the decomposition introduces an extra factor that is only controlled in L² rather than in the required higher integrability, the temporal order would drop.
minor comments (1)
  1. The notation L^∞_t L²_x ∩ L²_t W^{1,2}_x is standard but should be accompanied by an explicit statement of the underlying probability space and the precise meaning of “almost 1/2” (e.g., any exponent <1/2).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comments point by point below.

read point-by-point responses
  1. Referee: [Abstract / pressure decomposition section] Abstract and the section introducing the stochastic pressure decomposition: the claim that this decomposition lifts the temporal rate from (almost) 1/4 to (almost) 1/2 rests on the assertion that the resulting pressure estimates allow the stochastic integrals and nonlinear remainders to retain the higher order under linear-growth multiplicative noise. The manuscript must supply the explicit decomposition (including any Itô corrections) and the corresponding a-priori bounds; without these, it is impossible to confirm that the linear-growth term does not force a reversion to the lower rate obtained in Carelli-Prohl (2012).

    Authors: The stochastic pressure decomposition, including the Itô correction, is stated explicitly in Section 3.2 (equation (3.8)) of the manuscript. The corresponding a-priori bounds appear in Lemma 3.4 and Proposition 3.5, which establish the necessary integrability to absorb the linear-growth multiplicative noise without order reduction in the stochastic integrals. We will revise the abstract and add a short clarifying paragraph in the introduction to emphasize how these bounds preserve the temporal rate. revision: yes

  2. Referee: [Main theorem / error analysis] The error analysis (likely the main theorem and its proof): the convergence-in-probability statement is stated for the combined space-time error, but the proof sketch must isolate the temporal contribution arising from the pressure term and show that the almost-1/2 rate survives after applying the Burkholder-Davis-Gundy inequality and Gronwall-type arguments to the linear-growth noise. If the decomposition introduces an extra factor that is only controlled in L² rather than in the required higher integrability, the temporal order would drop.

    Authors: In the proof of Theorem 4.1 the temporal pressure contribution is isolated immediately after the decomposition is substituted (see the estimates between (4.11) and (4.14)). The Burkholder-Davis-Gundy inequality is then applied to the resulting martingale terms, and the linear-growth factors are controlled by the higher-integrability bounds already obtained in Proposition 3.5. The subsequent Gronwall argument yields the almost-1/2 rate. We will insert an additional remark in the proof to make this isolation and the preservation of the rate explicit. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation relies on external citation and standard stochastic analysis techniques.

full rationale

The paper's central result improves the temporal convergence rate from (almost) 1/4 in the cited Carelli-Prohl (2012) work to (almost) 1/2 via a stochastic pressure decomposition. This prior work has non-overlapping authors and is invoked only to benchmark the improvement, not to justify a uniqueness theorem or ansatz. No self-citations appear in the provided text, no parameters are fitted to data and then renamed as predictions, and the pressure decomposition is presented as an original analytical step rather than a self-definitional or smuggled-in construction. The derivation chain therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

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

pith-pipeline@v0.9.0 · 5706 in / 995 out tokens · 27807 ms · 2026-05-25T14:27:27.302080+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

29 extracted references · 29 canonical work pages

  1. [1]

    Bensoussan and R

    A. Bensoussan and R. Temam. ´Equations stochastiques du type Navier–Stokes. J. Funct. A nal. 13, 195–222, 1973

  2. [2]

    Birnir: The Kolmogorov–Obukhov Statistical Theory o f Turbulence

    B. Birnir: The Kolmogorov–Obukhov Statistical Theory o f Turbulence. J. Nonlinear Sci. 23, 657–688. (2013) 20 DOMINIC BREIT & ALAN DODGSON

  3. [3]

    Bessaih, Z

    H. Bessaih, Z. Brze´ zniak, A. Millet (2014): Splitting up method for the 2D stochastic Navier-Stokes equations. Stochastic PDE: Anal. Comput. 2, 433–470

  4. [4]

    Bessaih, A

    H. Bessaih, A. Millet (2018): Strong L2 convergence of time numerical schemes for the stochastic two-dimensional Navier–Stokes equations IMA J. Num. Anal. doi:10.1093/imanum/dry058

  5. [5]

    Breit: Existence theory for stochastic power law fluid s

    D. Breit: Existence theory for stochastic power law fluid s. J. Math. Fluid Mech. 17, 295–326. (2015)

  6. [6]

    Breit: Existence theory for generalized Newtonian fluids

    D. Breit: Existence theory for generalized Newtonian fluids. Mathematics in Science and Engineering. Elsevier/Academic Press, London. (2017)

  7. [7]

    Brezzi and M

    F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element M ethods, Springer Ser. Comput. Math. 15, Springer-Verlag, New York, 1991

  8. [8]

    Brze´ zniak, E

    Z. Brze´ zniak, E. Carelli, J. A. Prohl (2013): Finite-element-based discretizations of the incompressi ble Navier–Stokes equations with multiplicative random forci ng, IMA J. Num. Anal. 33, 771–824

  9. [9]

    Carelli, A

    E. Carelli, A. Prohl (2012): Rates of convergence for discretizations of the stochastic incompressible Navier-Stokes equations. SIAM J. Numer. Anal. 50(5), 2467–2496

  10. [10]

    Carelli, E

    E. Carelli, E. Hausenblas, A. Prohl (2012): Time-splitting methods to solve the stochastic incompress ible Stokes equation. SIAM J. Numer. Anal. 50, 2917–2939

  11. [11]

    Capi´ nski, A note on uniqueness of stochastic Navier -Stokes equations, Univ

    M. Capi´ nski, A note on uniqueness of stochastic Navier -Stokes equations, Univ. Iagell. Acta Math. 30 (1993), 219–228

  12. [12]

    M.Capi´ nski, N. J. Cutland, Stochastic Navier-Stokes equations, Acta Appl. Math. 25 (1991), 59–85

  13. [13]

    Da Prato, J

    G. Da Prato, J. Zabczyk (1992): Stochastic Equations in Infinite Dimensions, Encyclopedia Math. Appl., vol. 44, Cambridge University Press, Cambridge

  14. [14]

    Flandoli and D

    F. Flandoli and D. Gatarek Martingale and stationary so lutions for stochastic Navier-Stokes equations. Probab. Theory Relat. Fields , 102, 367–391, 1995

  15. [15]

    & Raviart, P

    Girault,V. & Raviart, P. A. (1981) Finite Element Metho d for Navier–Stokes Equations: Theory and Algorithms. Berlin/Heidelberg/New York: Springer

  16. [16]

    Girault and J.L

    V. Girault and J.L. Lions, Two-grid finite-element sche mes for the steady Navier-Stokes problem in polyhedra, Port. Math. (N.S.), 58 (2001), pp. 25–57

  17. [17]

    Girault and L.R

    V. Girault and L.R. Scott, A quasi-local interpolation operator preserving the discrete divergence, Calcolo, 40 (2003), pp. 1–19

  18. [18]

    Heywood and R

    J.G. Heywood and R. Rannacher, Finite element approxim ation of the nonstationary Navier–Stokes problem. I. Regularity of solutions and second-order error estimates for spatial discretization. SIAM J. Numer. Anal. 19, 275–311, 1982

  19. [19]

    Hofmanov´ a (2013): Degenerate Parabolic Stochasti c Partial Differential Equations

    M. Hofmanov´ a (2013): Degenerate Parabolic Stochasti c Partial Differential Equations. Stoch. Pr. Ap. 123 (12), 4294–4336

  20. [20]

    Kuksin and A

    S. Kuksin and A. Shirikyan. Mathematics of two-dimensional turbulence , volume 194 of Cambridge Tracts in Mathematics . Cambridge University Press, Cambridge, 2012

  21. [21]

    Mikulevicius and B

    R. Mikulevicius and B. L. Rozovskii. Stochastic Navier -Stokes equations for turbulent flows. SIAM J. Math. Anal., 35(5):1250–1310, 2004

  22. [22]

    Pardoux (1975): Equations aux d´ eriv´ ees Partielles stochastiques non lin´ eaires monotones

    E. Pardoux (1975): Equations aux d´ eriv´ ees Partielles stochastiques non lin´ eaires monotones. Etude de solutions fortes de type Itˆ o, Ph.D. thesis, Universit´ e Paris Sud

  23. [23]

    Pr´ evˆ ot, M

    C. Pr´ evˆ ot, M. R¨ockner (2007): A concise course on stochastic partial differ ential equations. Lecture Notes in Mathematics, 1905. Springer, Berlin

  24. [24]

    Printems (2001): On the discretization in time of parabolic stochastic parti al differential equations

    J. Printems (2001): On the discretization in time of parabolic stochastic parti al differential equations. Math. Mod. Numer. Anal. 35, 1055–1078

  25. [25]

    M. Romito. Some probabilistic topics in the Navier–Sto kes equations. Recent progress in the theory of the Euler and Navier–Stokes equations. London Math. Soc. Lecture Note Ser. 430, 175–232, Cambridge Univ. Press, Cambridge, 2016

  26. [26]

    L. R. Scott and S. Zhang, Finite element interpolation of nonsmooth functions satis fying boundary conditions, Math. Comp. 54 (1990), no. 190, 483–493

  27. [27]

    Yan (2004): Semidiscrete Galerkin Approximation fo r a Linear Stochastic Parabolic Partial Differ- ential Equation Driven by an Additive Noise

    Y. Yan (2004): Semidiscrete Galerkin Approximation fo r a Linear Stochastic Parabolic Partial Differ- ential Equation Driven by an Additive Noise. Num. Math. 44, 8 29–847

  28. [28]

    Yan (2005): Galerkin finite element methods for stoch astic parabolic partial differential equations, SIAM J

    Y. Yan (2005): Galerkin finite element methods for stoch astic parabolic partial differential equations, SIAM J. Numer. Anal., 43, pp. 1363–1384

  29. [29]

    Yoshida (2012): Stochastic Shear thickening fluids: strong convergence of the Galerkin approximation and the energy inequality

    N. Yoshida (2012): Stochastic Shear thickening fluids: strong convergence of the Galerkin approximation and the energy inequality. The Annals of Applied Probabilit y 2012, Vol. 22, No. 3, 1215–1242 Department of Mathematics, Heriot-W att University, Riccart on Edinburgh EH14 4AS, UK E-mail address : d.breit@hw.ac.uk, ad335@hw.ac.uk