Stability of Poiseuille Flow of Navier-Stokes Equations on $\mathbb{R}^2$

Zhile Li

Stability of Poiseuille Flow of Navier-Stokes Equations on $\mathbb{R}^2$

T0 review · 0 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper proves nonlinear stability of the Poiseuille shear flow on $\mathbb{R}^2$: initial vorticity perturbations of size at most $\delta\nu^{7/3}$ stay within twice their size, and linearized high-frequency perturbations decay at rate…

desk verdict A credible extension of the Couette-flow energy method to Poiseuille flow on R^2; the reviewer's main objection does not survive contact with the weighted identity in Section 2. read the letter →

arxiv 2411.19716 v4 pith:JNKC3IWG submitted 2024-11-29 math.AP

Zhile Li This is my paper · ORCID

classification math.AP MSC 35Q3076E0535B35

keywords Navier–StokesequationsPoiseuilleflowenhanceddissipationnonlinearstabilityvorticityanisotropicSobolevspaceunboundeddomain

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish that the Poiseuille shear flow $U=(y^2,0)$ on the whole plane $\mathbb{R}^2$ is quantitatively stable: if the initial vorticity perturbation has size at most $\delta\nu^{7/3}$ in a weighted anisotropic Sobolev space, then the perturbation energy stays bounded by twice its initial value for all time. It also proves enhanced dissipation for the linearized flow, with $x$-frequencies $|k|\ge\nu^{-1/3}$ decaying at rate $\nu^{1/2}|k|^{1/2}$, which is faster than the heat equation. The contribution is a parameter-free $\nu$-scaling threshold for nonlinear stability on an unbounded domain, extending results that were previously available mostly on periodic or bounded domains. The argument works by building a frequency-dependent energy that includes stream-function terms and by absorbing all nonlinear error terms into the dissipation.

What carries the argument

The machinery is a frequency-dependent energy-dissipation pair $(E_k,D_k)$ whose coefficients $\alpha_k,\beta_k,\gamma_k$ switch at the cutoff $|k|=\nu^{-1/3}$. The linear evolution satisfies $E_k\approx \|\omega_k\|_2^2+\alpha_k\|\nabla_k\omega_k\|_2^2+\gamma_k\|y\omega_k\|_2^2+\gamma_k\|\partial_y\psi_k\|_2^2$ and $\frac{d}{dt}E_k\le -4cD_k-4c\lambda_kE_k$. The nonlinear argument sums these with weights $\langle k\rangle^{2m}\langle c\lambda_k t\rangle^{2J}$ and a correction factor $M_k(t)$ whose ODE absorbs derivatives of the time multiplier; the low-frequency control of the stream function is what makes the nonlinear terms finite.

What would settle it

A direct Fourier check settles the closure: at y-frequency $\eta = k$, $|k|\,\lVert \psi_k\rVert_2\,\lVert \partial_y \psi_k\rVert_2 = (4|k|^2)^{-1}\lVert \omega_k\rVert_2^2$, so the ratio blows up as $k \to 0$. Because this inequality is what absorbs the nonlinear term $T_6$ in Section 3.5, observing the ratio diverge on low-frequency modes would refute the closing argument and hence the proof of Theorem 1.1.

Watch

Extended reading notes

Core claim

The central claim is Theorem 1.1: for $J\ge1$ and $m>3/4$, if the initial weighted norm $\epsilon$ defined from $\omega_{\mathrm{in}}$ is at most $\delta\nu^{7/3}$, then the solution satisfies $E(t)\le 2E(0)$ for all $t$, with the same bound after applying the multiplier $\langle c\lambda^{\mathrm{pl}}_\nu(\partial_x)t\rangle^J$. For the linearized equation the energy $E_k$ obeys $\frac{d}{dt}E_k\le -4cD_k-4c\lambda_k E_k$, giving the decay rate $\lambda_k=\nu^{1/2}|k|^{1/2}$ for $|k|\ge\nu^{-1/3}$. The nonlinear proof sums these energies and establishes the bootstrap inequality $E(t)\le 2E(0)-4cD(t)+C\nu^{-7/6}D(t)\sup_{s\in[0,t]}E(s)^{1/2}$, which closes when $E(0)\le c^2C^{-2}\nu^{7/3}$.

Load-bearing premise

The load-bearing premise is a low-frequency bound on the stream function: the proof needs $|k|\,\lVert\psi_k\rVert_\infty^2 \lesssim \lVert\omega_k\rVert_2^2$ to absorb the term $T_6$ in Section 3.5, and this is the point at which Fourier modes with $k\to0$ can make the argument fail.

Editorial extensions

If this is right

Above the cutoff $|k|\ge\nu^{-1/3}$, linearized vorticity decays like $e^{-c\nu^{1/2}|k|^{1/2}t}$, so high x-frequency structures relax on the time scale $\nu^{-1/2}|k|^{-1/2}$, faster than the heat equation's $\nu^{-1}$ scale.
The nonlinear theorem yields a quantitative stability threshold with exponent $\gamma=7/3$ in the stability definition.
The bootstrap inequality gives control of the dissipation integral $D(t)$ by the initial energy, not just pointwise energy bounds.
The result holds on the unbounded domain $\mathbb{R}^2$ with no periodicity in $x$, using an anisotropic Sobolev norm adapted to the shear flow.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

An editorial inference: if the Section 3.5 stream-function estimate needs repair, the natural fix is an additional low-frequency weight on $\psi$, which would likely change the power of $\nu$ in the threshold.
The multiplier construction suggests a template for other power-law shear profiles $y^p$ on $\mathbb{R}^2$: choose the weights from $\lambda\sim\nu^{1/2}|k|^{1/2}$ and isolate the stream-function coupling in a separate supremum-in-$k$ term.
One testable extension is to run the closure inequality on a single Fourier mode with small $k$; the ratio $(4|k|^2)^{-1}$ would show immediately whether the absorption argument can hold without modification.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Quick take: this is a real extension of the Arbon–Bedrossian energy method to Poiseuille flow on the whole plane, and the specific objection raised in the review—the alleged failure of the stream-function bound in §3.5—does not land. The paper deserves a serious referee, not a desk reject.

What is actually new: linear enhanced dissipation at rate ν^{1/2}|k|^{1/2} for |k| ≥ ν^{-1/3}, and a nonlinear stability threshold of ν^{7/3} in an anisotropic norm. The linear energy identities in Section 2 check out, and the piecewise definition of λ_k is handled consistently. The extra ∂_xψ term in the vorticity equation is the genuine new difficulty, and the paper addresses it by putting γ_k||∂_yψ_k||^2 and γ_k||yω_k||^2 into the energy. That structure is what makes the nonlinear closure plausible.

On the reader's strongest claim: the counterexample for §3.5 omits the y-weighted term. The identity from the proof of Theorem 2.1 gives ||kψ_k||_2^2 ≤ 2||∂_yψ_k||_2^2 + ||yω_k||_2^2, and therefore |k| ||ψ_k||_2 ||∂_yψ_k||_2 is controlled by ||∂_yψ_k||_2^2 + ||yω_k||_2^2. In the low-frequency region γ_k = ν^{-2/3}, and the energy stores exactly those two weighted pieces. So the disputed T6,HL,L,HL estimate closes; an unweighted single-mode calculation is not the relevant quantity.

Soft spots are minor. In §3.4 the I3 factor is written as ν^{-1/6}; a direct low-frequency estimate gives ν^{-1/3}, but that changes the final bound from ν^{-1/2} to ν^{-2/3}, still well below the dominant ν^{-7/6}. The paper also does not discuss local well-posedness or the k = 0 mode in detail; for this style of energy argument that is an expositional gap, not a fatal flaw. The constants in Theorem 2.1 are stated with inequalities that ensure positivity but the proof does not explicitly show D_k ≲ I; choosing cα small and cγ large makes it true, so this is also repairable.

The citation pattern is fine: the method is openly borrowed from [18], and the new flow is the contribution. This paper is for researchers working on quantitative shear-flow stability, and it gives them a solid new example with the same threshold machinery. I would send it to peer review and ask the author to tighten the few loose estimates and add a word on well-posedness.

Referee Report

0 major / 5 minor

Summary. The paper studies the 2D Navier-Stokes equations linearized and nonlinearized around the Poiseuille flow (y^2,0) on the whole plane. For the linearized vorticity equation, the author proves an energy inequality in which each x-frequency k has a dissipation rate λ_k = ν^{1/2}|k|^{1/2} for |k| ≥ ν^{-1/3} and λ_k = ν|k|^2 for small |k|, yielding enhanced dissipation at high frequencies. The energy E_k includes the weighted quantities ||yω_k||_2 and ||∇ψ_k||_2, which are needed because the linearized equation contains the extra term 2ikψ_k. For the nonlinear problem, the paper defines a time-weighted energy E(t) with a bootstrap multiplier ⟨cλ(∂x)t⟩^J and proves that if an anisotropic Sobolev norm of the initial vorticity perturbation is at most δν^{7/3}, then E(t) ≤ 2E(0) for all t ≥ 0. The proof is a continuation argument in which all nonlinear terms are bounded by C ν^{-7/6} D(t) sup_{s≤t} E(s)^{1/2}.

Significance. If the bounds are correct, this is a solid contribution: it appears to give the first quantitative nonlinear stability threshold for Poiseuille flow in the unbounded domain R^2, extending the Arbon-Bedrossian approach for Couette flow to a setting with a non-constant background shear and the additional stream-function coupling 2∂_xψ. The linear energy identities are explicit, the decay rate is derived rather than fitted, and the nonlinear estimates are fully itemized; there are no free parameters and no circular construction. The claimed threshold ν^{7/3} is very small, but the paper's goal is qualitative stability rather than optimality. The stress-test objection to §3.5 does not land: the unweighted counterexample with a single y-Fourier mode is bypassed by the weighted identity in the proof of Theorem 2.1, so the low-frequency nonlinear closure is valid once that identity is invoked explicitly.

minor comments (5)

[§3.5] The bound on T_{6,HL,L,HL} uses the step |k| ||ψ_k||_∞^2 ≲ |k| ||ψ_k||_2 ||∂_yψ_k||_2 ≲ integrated against ν^{-2/3}γ_k. As written, the second inequality is not immediate: for a single y-Fourier mode with η = k, one has |k| ||ψ_k||_2 ||∂_yψ_k||_2 = (4k^2)^{-1} ||ω_k||_2^2, which is unbounded as k → 0. The intended bound follows from the identity ||kψ_k||_2^2 - ||∂_yψ_k||_2^2 = 2 Re⟨∂_yψ_k, yω_k⟩ proved in Theorem 2.1, which gives |k| ||ψ_k||_2 ||∂_yψ_k||_2 ≤ C(||∂_yψ_k||_2^2 + ||yω_k||_2^2). Please cite this identity explicitly at the point of use in §3.5.
[§3.1 / Lemma 3.4] In the proof of Lemma 3.4, the replacement ||∂_y∇_kω_k||_2 ≤ ||Δ_kω_k||_2 is used without comment. This is valid in y-Fourier variables, but it is a non-obvious step; adding one sentence with the Fourier multiplier comparison would improve readability.
[§3.4] In the bound on I3 in the estimate of T_{5,HL,H''}, the inequality ∫⟨k⟩^{-2m}(γ_kν)^{-1}dk ≲ ν^{-1/3} requires splitting the integral into |k| ≥ ν^{-1/3} and |k| < ν^{-1/3}; the two regimes give ν^{-1/6} and ν^{-1/3} respectively, so the low-frequency part dominates. The text currently says this follows from γ_k ≥ ν^{-2/3} alone, which is slightly terse; please spell out the split.
[Theorem 1.1] The statement refers to 'the corresponding solution' of the nonlinear equation but does not cite or state a local well-posedness theorem for (1.3) on the relevant anisotropic spaces. Adding a reference or a short existence/uniqueness remark would make the theorem self-contained.
[Throughout] There are several typographical errors, e.g. 'eatimates' in §3.1, 'obatain' in §3.5, and 'separat' in §3.5. These should be corrected in a final pass.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the derivation is self-contained; Theorem 2.1's decay estimate is proven from the linearized PDE, and the nonlinear bootstrap closes the same energy with no fitted parameters or self-citation chain.

full rationale

The paper contains no fitted parameters, no empirical data, and no load-bearing self-citation: the only methodological citation, [18] (Arbon and Bedrossian), is external and used as inspiration, while all weighted identities and inequalities are proved in Section 2. The linear decay rate λ_k is not assumed: Theorem 2.1 derives dE_k/dt ≤ -4cD_k - 4cλ_kE_k from Lemma 2.3's exact time-derivative identities and from explicit verification that λ_kE_k ≲ I in both frequency regimes. The nonlinear stability statement in Theorem 1.1 has the same norms on both sides by design, as a standard a priori bootstrap, and Theorem 3.7 closes it via the contradiction argument from (3.5), with the threshold δν^{7/3} arising from the ν-exponents in the estimates of T1–T8, not from a fit. The possible failure of the §3.5 estimate flagged by the reader is a mathematical correctness concern, not a circularity: no equation is defined in terms of the conclusion, and no error estimate is renamed as a prediction. Therefore no circular step is present.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

No new physical entities and no fitted parameters are introduced. The central claim rests on standard analytic tools plus the specific anisotropic smallness class; the fragile input is the stream-function estimate near k=0.

assumptions (3)

standard math 2D Navier-Stokes global well-posedness for the perturbed vorticity equation on R^2
The proof freely uses time-continuous energy evolution and global solutions; this is standard for 2D Navier-Stokes.
domain assumption Elliptic estimates for Δ_k^{-1} on R_y in the k→0 limit, including the mean-zero condition ∫ω_0 = 0 required for ∂_yψ_0 ∈ L^2
The theorem's norms include ∂_yΔ^{-1}ω and ∇ψ; finiteness implicitly restricts the k=0 mode. This is not proven from the PDE but assumed by the choice of the norm.
standard math Gagliardo-Nirenberg-Sobolev inequalities in one dimension applied mode-by-mode in y
Used in Lemmas 3.4, 3.5 and 3.6; these are standard, but the specific application to ψ_k near k=0 in Section 3.5 is not valid without an additional argument.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Stability of Poiseuille Flow of Navier-Stokes Equations on $\mathbb{R}^2$." pith.science (2026). https://pith.science/paper/JNKC3IWG

@misc{pith2026241119716,
  author       = {Pith},
  title        = {Pith review of: Stability of Poiseuille Flow of Navier-Stokes Equations on $\mathbbR^2$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JNKC3IWG}},
  note         = {Machine review of arXiv:2411.19716}
}

abstract

We consider solutions to the Navier-Stokes equations on $\mathbb{R}^2$ close to the Poiseuille flow with viscosity $0< \nu < 1$. For the linearized problem, we prove that when the $x$-frequency satisfy $|k| \ge \nu^{-\frac{1}{3}}$, the perturbation decays on a time-scale proportional to $\nu^{-\frac{1}{2}}|k|^{-\frac{1}{2}}$. Since it decays faster than the heat equation, this phenomenon is referred to as enhanced dissipation. Then we concern the non-linear equations. We show that if the initial perturbation $\omega_{in}$ is at most of size $\nu^\frac{7}{3}$ in an anisotropic Sobolev space, then the size of the perturbation remains no more than twice the size of its initial value.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 15 canonical work pages

[18]

Quantitative hydrodynamic stability for couette flow on unbounded domains with navier boundary conditions,

R. Arbon and J. Bedrossian, “Quantitative hydrodynamic stability for couette flow on unbounded domains with navier boundary conditions,” arXiv preprint arXiv:2404.02412 , 2024. 18

arXiv 2024
[1]

Stability of the couette flow at high reynolds numbers in two dimensions and three dimensions,

J. Bedrossian, P. Germain, and N. Masmoudi, “Stability of the couette flow at high reynolds numbers in two dimensions and three dimensions,” Bulletin of the American Mathematical Society , vol. 56, no. 3, pp. 373–414, 2019

work page 2019
[2]

Rapid relaxation of an axisymmetric vortex,

A. J. Bernoff and J. F. Lingevitch, “Rapid relaxation of an axisymmetric vortex,” Physics of Fluids, vol. 6, no. 11, pp. 3717–3723, 1994

work page 1994
[3]

Transient anomalous diffusion in poiseuille flow,

M. Latini and A. J. Bernoff, “Transient anomalous diffusion in poiseuille flow,”Journal of Fluid Mechanics, vol. 441, pp. 399–411, 2001

work page 2001
[4]

Metastability and rapid convergence to quasi-stationary bar states for the two-dimensional navier–stokes equations,

M. Beck and C. E. Wayne, “Metastability and rapid convergence to quasi-stationary bar states for the two-dimensional navier–stokes equations,” Proceedings of the Royal Society of Edinburgh Section A: Math- ematics, vol. 143, no. 5, pp. 905–927, 2013

work page 2013
[5]

Enhanced dissipation, hypoellipticity, and anomalous small noise inviscid limits in shear flows,

J. Bedrossian and M. Coti Zelati, “Enhanced dissipation, hypoellipticity, and anomalous small noise inviscid limits in shear flows,” Archive for Rational Mechanics and Analysis , vol. 224, no. 3, pp. 1161– 1204, 2017

work page 2017
[6]

Invariant measures for passive scalars in the small noise inviscid limit,

J. Bedrossian, M. Coti Zelati, and N. Glatt-Holtz, “Invariant measures for passive scalars in the small noise inviscid limit,” Communications in Mathematical Physics , vol. 348, pp. 101–127, 2016

work page 2016
[7]

Enhanced dissipation and inviscid damping in the invis- cid limit of the navier–stokes equations near the two dimensional couette flow,

J. Bedrossian, N. Masmoudi, and V. Vicol, “Enhanced dissipation and inviscid damping in the invis- cid limit of the navier–stokes equations near the two dimensional couette flow,” Archive for Rational Mechanics and Analysis , vol. 219, pp. 1087–1159, 2016

work page 2016

Show all 18 references

[8]

Transition threshold for the 3d couette flow in sobolev space,

D. Wei and Z. Zhang, “Transition threshold for the 3d couette flow in sobolev space,” Communications on Pure and Applied Mathematics , vol. 74, no. 11, pp. 2398–2479, 2021

2021
[9]

Pseudospectral and spectral bounds for the oseen vortices operator,

T. Li, D. Wei, and Z. Zhang, “Pseudospectral and spectral bounds for the oseen vortices operator,” arXiv preprint arXiv:1701.06269, 2017. 17

2017 arXiv
[10]

Enhanced dissipation and taylor dispersion in higher-dimensional parallel shear flows,

M. Coti Zelati and T. Gallay, “Enhanced dissipation and taylor dispersion in higher-dimensional parallel shear flows,” Journal of the London Mathematical Society , vol. 108, no. 4, pp. 1358–1392, 2023

2023
[11]

Stability of fluid motion: rectilinear motion of viscous fluid between two parallel plates,

L. Kelvin, “Stability of fluid motion: rectilinear motion of viscous fluid between two parallel plates,” Phil. Mag, vol. 24, no. 5, pp. 188–196, 1887
[12]

Enhanced dissipation in the navier–stokes equations near the poiseuille flow,

M. Coti Zelati, T. M. Elgindi, and K. Widmayer, “Enhanced dissipation in the navier–stokes equations near the poiseuille flow,” Communications in Mathematical Physics , vol. 378, no. 2, pp. 987–1010, 2020

2020
[13]

The sobolev stability threshold for 2d shear flows near couette,

J. Bedrossian, V. Vicol, and F. Wang, “The sobolev stability threshold for 2d shear flows near couette,” Journal of Nonlinear Science , vol. 28, pp. 2051–2075, 2018

2018
[14]

A dynamical approach to the study of instability near couette flow,

H. Li, N. Masmoudi, and W. Zhao, “A dynamical approach to the study of instability near couette flow,” Communications on Pure and Applied Mathematics , vol. 77, no. 6, pp. 2863–2946, 2024

2024
[15]

Stability threshold of two-dimensional couette flow in sobolev spaces,

N. Masmoudi and W. Zhao, “Stability threshold of two-dimensional couette flow in sobolev spaces,” in Annales de l’Institut Henri Poincar´ e C, Analyse non lin´ eaire, vol. 39, no. 2, 2022

2022
[16]

Stability threshold of nearly-couette shear flows with navier boundary conditions in 2d,

J. Bedrossian, S. He, S. Iyer, and F. Wang, “Stability threshold of nearly-couette shear flows with navier boundary conditions in 2d,” arXiv preprint arXiv:2311.00141 , 2023

2023 arXiv
[17]

Transition threshold for the 2-d couette flow in a finite channel,

Q. Chen, T. Li, D. Wei, and Z. Zhang, “Transition threshold for the 2-d couette flow in a finite channel,” Archive for rational mechanics and analysis , vol. 238, no. 1, pp. 125–183, 2020

2020

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [18]

Quantitative hydrodynamic stability for couette flow on unbounded domains with navier boundary conditions,

R. Arbon and J. Bedrossian, “Quantitative hydrodynamic stability for couette flow on unbounded domains with navier boundary conditions,” arXiv preprint arXiv:2404.02412 , 2024. 18

arXiv 2024

[2] [1]

Stability of the couette flow at high reynolds numbers in two dimensions and three dimensions,

J. Bedrossian, P. Germain, and N. Masmoudi, “Stability of the couette flow at high reynolds numbers in two dimensions and three dimensions,” Bulletin of the American Mathematical Society , vol. 56, no. 3, pp. 373–414, 2019

work page 2019

[3] [2]

Rapid relaxation of an axisymmetric vortex,

A. J. Bernoff and J. F. Lingevitch, “Rapid relaxation of an axisymmetric vortex,” Physics of Fluids, vol. 6, no. 11, pp. 3717–3723, 1994

work page 1994

[4] [3]

Transient anomalous diffusion in poiseuille flow,

M. Latini and A. J. Bernoff, “Transient anomalous diffusion in poiseuille flow,”Journal of Fluid Mechanics, vol. 441, pp. 399–411, 2001

work page 2001

[5] [4]

Metastability and rapid convergence to quasi-stationary bar states for the two-dimensional navier–stokes equations,

M. Beck and C. E. Wayne, “Metastability and rapid convergence to quasi-stationary bar states for the two-dimensional navier–stokes equations,” Proceedings of the Royal Society of Edinburgh Section A: Math- ematics, vol. 143, no. 5, pp. 905–927, 2013

work page 2013

[6] [5]

Enhanced dissipation, hypoellipticity, and anomalous small noise inviscid limits in shear flows,

J. Bedrossian and M. Coti Zelati, “Enhanced dissipation, hypoellipticity, and anomalous small noise inviscid limits in shear flows,” Archive for Rational Mechanics and Analysis , vol. 224, no. 3, pp. 1161– 1204, 2017

work page 2017

[7] [6]

Invariant measures for passive scalars in the small noise inviscid limit,

J. Bedrossian, M. Coti Zelati, and N. Glatt-Holtz, “Invariant measures for passive scalars in the small noise inviscid limit,” Communications in Mathematical Physics , vol. 348, pp. 101–127, 2016

work page 2016

[8] [7]

Enhanced dissipation and inviscid damping in the invis- cid limit of the navier–stokes equations near the two dimensional couette flow,

J. Bedrossian, N. Masmoudi, and V. Vicol, “Enhanced dissipation and inviscid damping in the invis- cid limit of the navier–stokes equations near the two dimensional couette flow,” Archive for Rational Mechanics and Analysis , vol. 219, pp. 1087–1159, 2016

work page 2016