Pith. sign in

REVIEW 2 major objections 3 minor 28 references

Non-uniqueness of weak solutions to 2D hypoviscous Navier-Stokes equations

T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read This paper proves that any smooth zero-mean divergence-free flow on the 2D torus is $L^1$-approximable by a weak solution of the hypoviscous Navier–Stokes equations, for every $\theta\in[0,1)$.

desk verdict Solid convex integration paper: proves an h-principle and compact temporal support for 2D hypoviscous Navier-Stokes below theta=1, even though plain non-uniqueness was already available from the authors' Boussinesq paper. read the letter →

arxiv 1908.06005 v2 pith:JR5PO67M submitted 2019-08-16 math.AP

Tianwen Luo , Peng Qu This is my paper · ORCID
classification math.AP MSC 35Q3035D3035R1176D0376D05
keywords non-uniquenessweaksolutionstwo-dimensionalhypoviscousNavier-StokesequationsconvexintegrationintermittentstationaryflowfractionalLaplacianh-principlecompacttemporalsupport
topics Navier-Stokes Regularity
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 establishes that the two-dimensional hypoviscous incompressible Navier–Stokes equations, with fractional dissipation exponent $\theta\in[0,1)$, have a highly flexible set of weak solutions. Its main theorem is an h-principle: any smooth divergence-free velocity field with zero spatial mean can be approximated as closely as desired, in the $L^\infty_tL^1_x$ norm, by a genuine weak solution that is continuous in time and square-integrable in space, with the solution's temporal support kept inside a small neighbourhood of the target field's temporal support. Choosing the target field to be supported on a short time interval then yields nontrivial weak solutions with compact temporal support, so weak solutions to the Cauchy problem are generally not unique. The proof works by a convex-integration iteration that adds high-frequency, spatially concentrated two-dimensional stationary waves to cancel the Reynolds stress step by step.

What carries the argument

The load-bearing object is the two-dimensional intermittent stationary flow $W_\xi(t,x)=\eta_\xi(t,x)b_{\xi,\lambda}(x)$, where $b_{\xi,\lambda}(x)=i\xi^\perp e^{i\lambda\xi\cdot x}$ is a stationary Euler flow and $\eta_\xi$ is a directed, rescaled Dirichlet kernel with a temporal shift that carries the concentration along characteristics. A geometric lemma decomposes every symmetric trace-free $2\times2$ matrix $\mathring R$ as a sum of squares $\sum_{\xi}(\gamma_\xi(\mathring R))^2(\xi\,\hat\otimes\,\xi)$, so the coefficients $a_\xi$ chosen from this decomposition make the principal perturbation cancel the current Reynolds stress. The argument hinges on the two-dimensional $L^p$ bound $\|D_r\|_{L^p}\lesssim r^{1-2/p}$ for the Dirichlet kernel, which is dimensionally different from the three-dimensional case and forces the parameter scaling $r=\lambda^{1-6\alpha}$, $\mu=\lambda^{1-4\alpha}$, $\sigma=\lambda^{-(1-2\alpha)}$; with these scales the Reynolds-stress estimates close and the iteration converges.

What would settle it

A reader could settle the mechanism by computing or rigorously bounding the exact $L^p$ norm of the two-dimensional Dirichlet kernel (4.9) and checking whether the bound $\|D_r\|_{L^p}\lesssim r^{1-2/p}$ really holds for all $1<p\le\infty$; if the norm grew like $r^{1-2/p+\delta}$ for any $\delta>0$, the parameter constraints (7.20) could not close. More directly, if for some $\theta\in[0,1)$ one produced a smooth zero-mean field $u$ that no weak solution $v$ can approximate in the sense of (1.3)–(1.4), Theorem 1.1 would be false.

Watch

Extended reading notes

Core claim

The central claim is Theorem 1.1: for every $\theta\in[0,1)$, every $T>0$, and every smooth zero-mean divergence-free field $u$ on $[0,T]\times\mathbb{T}^2$, and for every $\varepsilon_*>0$, there exists a weak solution $v\in C^0_tL^2_x$ to (1.1) with zero spatial mean such that $\|v-u\|_{L^\infty_tL^1_x}\le\varepsilon_*$ and the temporal support of $v$ is contained in the $\varepsilon_*$-neighbourhood of the temporal support of $u$. In particular, by taking $u$ with compact temporal support, the system admits nontrivial compactly supported weak solutions, and therefore the weak solutions of the Cauchy problem for (1.1) are not unique. The proof constructs $v$ as the strong limit of an iteration in $C^0_tH^{\beta'}$ for small $\beta'>0$, with each step adding an intermittent two-dimensional stationary wave that eliminates the current Reynolds stress while keeping the perturbation small in $L^2$.

Load-bearing premise

The iteration shrinks the Reynolds stress only if the two-dimensional building blocks concentrate exactly as fast as the Dirichlet-kernel bound $\|D_r\|_{L^p}\lesssim r^{1-2/p}$; if the true growth were worse by any positive power of $r$, the parameter choices in (7.20) would fail and the correction step would stop shrinking.

Editorial extensions

If this is right

  • For every $\theta\in[0,1)$, the Cauchy problem for (1.1) admits weak solutions in $C^0_tL^2_x$ that are not unique, because nontrivial compactly supported weak solutions exist.
  • The h-principle holds: the set of weak solutions is dense in the $L^\infty_tL^1_x$ topology in the space of smooth zero-mean divergence-free fields, with the temporal support of the approximation controlled by the chosen accuracy.
  • The constructed weak solutions are strong limits in $C^0_tH^{\beta'}_x$ for small $\beta'>0$, so they are genuine continuous-in-time $L^2$ functions, not merely formal distributional limits.
  • The result covers the whole subcritical range $\theta\in[0,1)$, including the damping case $\theta=0$; only the classical case $\theta=1$ is excluded, where weak solutions are known to be unique.

Reading between the lines

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

  • Editorial inference: the support-containment estimate suggests a time-patching procedure not stated in the paper: one could prescribe any smooth zero-mean flow on a short time interval, allow wild behaviour on a later interval, and then return to a smooth flow, because the proof controls exactly the temporal support of the approximation.
  • Editorial inference: the endpoint $\theta=1$ is not merely a technical cutoff, since for $\theta=1$ the two-dimensional weak solutions are unique; the construction therefore works uniformly below the classical threshold, and it remains open whether some smaller threshold inside $[0,1)$ separates flexibility from rigidity.
  • Editorial inference: because the proof's quantitative mechanism is the two-dimensional Dirichlet-kernel bound, one could test whether other concentration profiles with better $L^p$ growth would allow larger ranges of parameters or slightly smoother weak solutions within the same convex-integration framework.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 3 minor

Summary. The manuscript adapts the intermittent convex-integration scheme of Buckmaster and Vicol to the two-dimensional setting. For any fractional-viscosity exponent θ in [0,1), it claims an h-principle: every smooth divergence-free zero-mean field u on [0,T]×T^2 can be approximated in L∞_t L^1_x by a C^0_t L^2_x weak solution v to the 2D hypoviscous Navier–Stokes equations (1.1), with the temporal support of v contained in a prescribed neighbourhood of supp_t u. The core of the proof is an iteration lemma (Lemma 2.1) built on a 2D intermittent stationary flow obtained by modulating the stationary flows of Choffrut–De Lellis–Szekelyhidi with a rescaled Dirichlet kernel. A geometric lemma (Appendix A) and a sequence of a priori estimates for the perturbation and the Reynolds stress are supplied.

Significance. If the result were correct, it would be a substantial extension of high-dimensional convex-integration machinery to 2D and would yield nontrivial compactly-supported weak solutions and non-uniqueness for the full hypoviscous range θ∈[0,1). The paper is clearly written, the geometric lemma is explicit, and the iteration is organised so that the quantitative closure can be checked. Several auxiliary lemmas are imported from previous works, which is acceptable if the cited statements are standard. However, the central Reynolds-stress estimate contains a term with a positive power of the frequency parameter, so the iteration does not close as written; this affects the main theorem and corollary.

major comments (2)
  1. [§7, Eq. (7.18)] The second term on the right-hand side of (7.18), namely ℓ^{-4}λ_{q+1}^{θ*} r^{1-2/p}, is not controlled by the parameter choice (7.20). With ℓ=λ_q^{-20}=λ_{q+1}^{-20/B} and r^{1-2/p}=λ_{q+1}^{α/2}, this term equals λ_{q+1}^{80/B+θ*+α/2}. Its exponent is positive for every θ∈[0,1): for θ≤1/2 one has θ*=0 and the exponent is 80/B+α/2>0, while for θ>1/2 the exponent is even larger. The desired estimate (2.15) requires this term to be bounded by Aε_{q+2}=Aλ_{q+1}^{-2βB}, which is impossible for large q because the left-hand side grows while the right-hand side tends to zero. The origin of the difficulty is the bound in (7.11), where R((−Δ)^θ w_{q+1}) is estimated through ‖w_{q+1}‖_{L^p}^{1−θ*}‖∇w_{q+1}‖_{L^p}^{θ*} together with the lossy L^p bounds of Proposition 6.3; the λ^{-1} gained from the anti-divergence and the ε_{q+1}^{1/2} smallness of the coefficients are not exploited. This is a load-bearing gap in the proof of Lemma 2.1 and hence in Theorem 1.1.
  2. [§2 and §7, parameter ranges] The parameter choices contain two inconsistencies that are part of the same closure problem. First, the text after (7.20) asserts that p=(2−12α)/(2−13α) lies in (1,2) for α satisfying (2.3); this is false when α≥2/25, and (2.3) permits α as large as 1/8. Second, controlling the fractional-viscosity term in the natural sharpened form requires, for θ>1/2, the lower bound β(2B−1) ≥ 2θ−1, i.e. β≳θ*/(2B), whereas (2.7) imposes β<1/(100B^2). The two constraints are incompatible as stated. The parameter ranges and the exponent bookkeeping in (7.18)–(7.20) must be reworked before the iteration can close.
minor comments (3)
  1. [Proof of Theorem 1.1] The step from C^0_t H^{β'} convergence and L∞_t L^2 convergence to v∈C^0_t L^2 should be made explicit: because the approximating sequence is uniformly L^2-continuous and converges uniformly in L^2, the limit is also L^2-continuous. This is true but is not spelled out.
  2. [§2, Eq. (2.10)] Equation (2.10) writes ε_{q+2}^{-1}=λ_q^{2βB^2}; it would be clearer to state explicitly that λ_{q+2}=λ_q^{B^2} follows from (2.4).
  3. [Throughout] The text contains several typesetting artifacts, for example '/upslope' before integrals and 'Na vier-Stokes' in the title block; these should be cleaned in the final version.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the convex-integration construction is self-contained and the cited technical lemmas are external, not load-bearing self-citations.

full rationale

The paper's central claim is established by a convex-integration iteration whose quantitative closure is checked in the paper itself. Lemma 2.1 starts from a Reynolds stress bounded by Aε_{q+1} and produces a new stress bounded by Aε_{q+2}; the parameter choices in (7.20) and the subsequent exponent bookkeeping close the iteration without fitting any parameter to the target solution u or encoding the conclusion. The geometric lemma (Lemma 4.1) is proved in Appendix A, and the identity (4.21) follows algebraically from it and the normalization (4.13). The analytic tools imported from the literature — the L^p product estimate (Lemma 6.2) and the λ^{-1} frequency-localization estimate (Lemma 7.4) — are cited from external works [5, 23], not from the authors' own prior papers, and they are standard. The authors' self-citations to [20] and [21] are contextual: [20] is mentioned only to explain what is new in this note, and [21] is used only for the idea of a temporal cut-off function, which is then explicitly constructed and estimated in Section 5. No claim in the paper reduces to a fitted value, a renamed known result, or an unverified self-citation chain. Therefore no circular step is present.

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

All parameters in the iteration are chosen to satisfy inequalities, not fitted to any data. The background results cited from [5,7,13,23] are standard in convex integration. No new physical entities are introduced.

free parameters (4)
  • alpha = any rational in (0, (1-theta*)/8]
    Controls the frequency and amplitude scaling; fixes the L^p exponent p = (2-12alpha)/(2-13alpha) in the parameter choice (7.20).
  • B = any integer > 320/alpha
    Sets the growth rate lambda_q = A^{B^q} so that consecutive frequencies are far apart and error decays fast.
  • beta = any real in (0, 1/(100 B^2))
    Sets the amplitude epsilon_q = lambda_q^{-2beta}; needs to be small so the series converges in C^0_t H^{beta'} and the limit is a weak solution.
  • A = any sufficiently large multiple of 5 with A^alpha in 5N
    Absorbs absolute constants in the estimates and ensures lambda_q in 5N, which is used for the frequency localization of the building blocks.
assumptions (4)
  • standard math L^p product estimate (Lemma 6.2) and frequency localization estimates (4.16)-(4.18) from [5,23] hold.
    Central to the decoupling of high-frequency oscillations; the paper cites rather than proves them.
  • standard math The anti-divergence operator R satisfies the estimates in Lemma 7.3 and the lambda^{-1} gain in Lemma 7.4.
    These are standard Calderon-Zygmund and Schauder estimates, cited from [7,13,5].
  • standard math Fractional Laplacian interpolation for R(-Delta)^theta: ||R(-Delta)^theta w||_{L^p} <= C ||w||_{L^p}^{1-theta*} ||grad w||_{L^p}^{theta*} with theta* = max(2theta-1,0).
    Used in Section 7 before (7.11) to estimate the viscous term; follows from Fourier multiplier order and complex interpolation.
  • standard math Sobolev embedding H^{beta'} subset L^2 and strong convergence in C^0_t H^{beta'} imply the limit is in C^0_t L^2_x.
    Used in the proof of the main theorem to identify the limit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Non-uniqueness of weak solutions to 2D hypoviscous Navier-Stokes equations." pith.science (2026). https://pith.science/paper/JR5PO67M

@misc{pith2026190806005,
  author       = {Pith},
  title        = {Pith review of: Non-uniqueness of weak solutions to 2D hypoviscous Navier-Stokes equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JR5PO67M}},
  note         = {Machine review of arXiv:1908.06005}
}
abstract

Through an adaption of the convex integration scheme in the two dimensional case, the non-uniqueness of $C^0_t L^2_x$ weak solutions is presented for the two-dimensional hypoviscous incompressible Navier-Stokes equations.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 24 canonical work pages

  1. [20]

    T. Luo, T. Tao & L. Zhang, Finite energy weak solutions of 2d Boussinesq Equations wit h diffusive temperature, arXiv:1901.09179

    work page Pith review arXiv 1901

  2. [5]

    Buckmaster & V

    T. Buckmaster & V. Vicol, Nonuniqueness of weak solutions to the Navier–Stokes equat ion, Ann. of Math. 189:1 (2019), 101–144

    work page 2019

  3. [1]

    Buckmaster, M

    T. Buckmaster, M. Colombo & V. Vicol, Wild solutions of the Navier–Stokes equations whose singular sets in time have Hausdorff dimension strictl y less than 1 , arXiv:1809.00600

    arXiv

  4. [2]

    Buckmaster, C

    T. Buckmaster, C. De Lellis, P. Isett & L. Sz´ ekelyhidi, J r., Anomalous dissipation for 1/5-H¨ older Euler flows, Ann. of Math. 182:1 (2015), 127–172

    work page 2015

  5. [3]

    Buckmaster, C

    T. Buckmaster, C. De Lellis, L. Sz´ ekelyhidi & V. Vicol, Onsager’s conjecture for admissible weak solutions, Comm. Pure Appl. Math. 72:2 (2019), 229–274

    work page 2019

  6. [4]

    Buckmaster, C

    T. Buckmaster, C. De Lellis & L. Sz´ ekelyhidi, Dissipative Euler flows with Onsager-critical spatial regularity, Comm. Pure Appl. Math. 69:9 (2016), 1613–1670

    work page 2016

  7. [6]

    Cheskidov & X

    A. Cheskidov & X. Luo, Stationary and discontinuous weak solutions of the Navier– Stokes equations, arXiv:1901.07485

    arXiv 1901

  8. [7]

    Choffrut, C

    A. Choffrut, C. De Lellis & L. Sz´ ekelyhidi Jr., Dissipative continuous Euler flows in two and three dimensions , arXiv:1205.1226

    arXiv

Show all 28 references
  1. [8]

    Colombo, C

    M. Colombo, C. De Lellis & L. De Rosa, Ill-Posedness of Leray solutions for the hypodis- sipative Navier–Stokes equations , Comm. Math. Phys. 362:2 (2018), 659–688

    2018

  2. [9]

    Colombo, C

    M. Colombo, C. De Lellis & A. Massaccesi, The generalized Caffarelli–Kohn–Nirenberg theorem for the hyperdissipative Navier–Stokes system , arXiv:1712.07015

    arXiv

  3. [10]

    Constantin, W

    P. Constantin, W. E & E. S. Titi; Onsager’s conjecture on the energy conservation for solutions of Euler’s equation , Comm. Math. Phys. 165:1 (1994), 207–209

    1994

  4. [11]

    Caffarelli, R

    L. Caffarelli, R. Kohn & L. Nirenberg, Partial regularity of suitable weak solutions to the Navier–Stokes equations , Comm. Pure Appl. Math. 35 (1982), 771–831

    1982

  5. [12]

    De Lellis & L

    C. De Lellis & L. Sz´ ekelyhidi, Jr., The Euler equations as a differential inclusion . Ann. of Math. 170:3 (2009), 1417–1436

    2009

  6. [13]

    De Lellis & L

    C. De Lellis & L. Sz´ ekelyhidi Jr., Dissipative continu ous Euler flows, Invent. Math. 193(2), 2013, 377–407

    2013

  7. [14]

    Isett, H¨ older continuous Euler flows with compact support in time , Doctoral thesis, Princeton University, 2013

    P. Isett, H¨ older continuous Euler flows with compact support in time , Doctoral thesis, Princeton University, 2013

    2013

  8. [15]

    Isett, A Proof of Onsager’s Conjecture , Ann

    P. Isett, A Proof of Onsager’s Conjecture , Ann. of Math. 188:3 (2018), 1–93

    2018

  9. [16]

    Q. Jiu, Y. Wang, On possible time singular points and eventual regularity of weak solutions to the fractional Navier-Stokes equations , Dyn. Partial Differ. Equ. 11:4 (2014), 321–343

    2014

  10. [17]

    N. H. Katz, N. A. Pavlovi´ c, A cheap Caffarelli–Kohn–Nirenberg inequality for the Navier – Stokes equation with hyper-dissipation , Geom. Funct. Anal. 12:2 (2002), 355–379

    2002

  11. [18]

    J. L. Lions, Quelques r´ esultats d’existence dans des ´ equations aux d´eriv´ ees partielles non lin´ eaires. Bull. Soc. Math. France 87 (1959), 245–273

    1959

  12. [19]

    J. L. Lions, Quelques M´ ethodes de Resolution des Probl´ emes aux Limites Non lin´ eaires, Vol 1. Dunod, Paris, 1969

    1969

  13. [21]

    T. Luo & E. S. Titi; Non-uniqueness of weak solutions to hyperviscous Navier–S tokes equa- tions – on sharpness of J.-L. Lions exponent , arXiv:1808.07595

    arXiv

  14. [22]

    Xiaoyutao Luo, Stationary solution and nonuniquenes of weak solution for t he Navier- Stokes euation on high dimensions , Arch. Ration. Mech. Anal., 233:2 (2019), 701–747

    2019

  15. [23]

    Modena & L

    S. Modena & L. Sz´ ekelyhidi Jr.; Non-uniqueness for the transport equation with Sobolev vector fields , Ann. PDE 4:2 (2018), Art. 18, 38 pp

    2018

  16. [24]

    Olson, E

    E. Olson, E. S. Titi, Viscosity versus vorticity stretching: Global well-posed ness for a family of NavierStokes-alpha-like models , Nonlinear Anal. 66:11 (2007), 2427–2458

    2007

  17. [25]

    Temam, Navier–Stokes equations

    R. Temam, Navier–Stokes equations. Theory and numerical analysis , North Holland, Am- sterdam, 1977

    1977

  18. [26]

    Wu, Generalized MHD equations , J

    J. Wu, Generalized MHD equations , J. Differential Equations 195 (2003), 284–312

    2003

  19. [27]

    Tao, Global regularity for a logarithmically supercritical hyp erdissipative Navier-Stokes equation, Anal

    T. Tao, Global regularity for a logarithmically supercritical hyp erdissipative Navier-Stokes equation, Anal. PDE, 3 (2009), 361-366. 17

    2009

  20. [28]

    Tang & Y

    L. Tang & Y. Yu. Partial regularity of suitable weak solutions to the fracti onal Navier-Stokes equations. Comm. Math. Phys., 334(3):1455–1482, 2015. Yau Mathematical Sciences Center, Tsinghua University, Ch ina. E-mail address, T. Luo: twluo@mail.tsinghua.edu.cn School of Mat...

    2015

Pith tools

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