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
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.
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
- 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.
Referee Report
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)
- [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).
- [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)
- 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
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
-
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
-
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
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
Reference graph
Works this paper leans on
-
[1]
A. Bensoussan and R. Temam. ´Equations stochastiques du type Navier–Stokes. J. Funct. A nal. 13, 195–222, 1973
work page 1973
-
[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
work page 2013
-
[3]
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
work page 2014
-
[4]
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]
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)
work page 2015
-
[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)
work page 2017
-
[7]
F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element M ethods, Springer Ser. Comput. Math. 15, Springer-Verlag, New York, 1991
work page 1991
-
[8]
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
work page 2013
-
[9]
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
work page 2012
-
[10]
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
work page 2012
-
[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
work page 1993
-
[12]
M.Capi´ nski, N. J. Cutland, Stochastic Navier-Stokes equations, Acta Appl. Math. 25 (1991), 59–85
work page 1991
-
[13]
G. Da Prato, J. Zabczyk (1992): Stochastic Equations in Infinite Dimensions, Encyclopedia Math. Appl., vol. 44, Cambridge University Press, Cambridge
work page 1992
-
[14]
F. Flandoli and D. Gatarek Martingale and stationary so lutions for stochastic Navier-Stokes equations. Probab. Theory Relat. Fields , 102, 367–391, 1995
work page 1995
-
[15]
Girault,V. & Raviart, P. A. (1981) Finite Element Metho d for Navier–Stokes Equations: Theory and Algorithms. Berlin/Heidelberg/New York: Springer
work page 1981
-
[16]
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
work page 2001
-
[17]
V. Girault and L.R. Scott, A quasi-local interpolation operator preserving the discrete divergence, Calcolo, 40 (2003), pp. 1–19
work page 2003
-
[18]
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
work page 1982
-
[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
work page 2013
-
[20]
S. Kuksin and A. Shirikyan. Mathematics of two-dimensional turbulence , volume 194 of Cambridge Tracts in Mathematics . Cambridge University Press, Cambridge, 2012
work page 2012
-
[21]
R. Mikulevicius and B. L. Rozovskii. Stochastic Navier -Stokes equations for turbulent flows. SIAM J. Math. Anal., 35(5):1250–1310, 2004
work page 2004
-
[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
work page 1975
-
[23]
C. Pr´ evˆ ot, M. R¨ockner (2007): A concise course on stochastic partial differ ential equations. Lecture Notes in Mathematics, 1905. Springer, Berlin
work page 2007
-
[24]
J. Printems (2001): On the discretization in time of parabolic stochastic parti al differential equations. Math. Mod. Numer. Anal. 35, 1055–1078
work page 2001
-
[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
work page 2016
-
[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
work page 1990
-
[27]
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
work page 2004
-
[28]
Y. Yan (2005): Galerkin finite element methods for stoch astic parabolic partial differential equations, SIAM J. Numer. Anal., 43, pp. 1363–1384
work page 2005
-
[29]
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
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.