Pith. sign in

REVIEW 3 major objections 6 minor 33 references

Non-uniqueness of Oberbeck-Boussinesq System with Gravitational Field

T0 review · 3 major / 6 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 constructs force and heat-source terms for which the Oberbeck-Boussinesq system with Newtonian gravity has infinitely many weak solutions from zero data.

desk verdict Genuine extension of Navier-Stokes non-uniqueness to Boussinesq with gravitational forcing, but the key semigroup estimate has sign and scaling errors; worth a referee, not acceptable as-is. read the letter →

arxiv 2411.17929 v1 pith:GVBWVQYN submitted 2024-11-26 math.AP

Xiaoran Liu This is my paper · ORCID
classification math.AP MSC 35L7135B4037K40
keywords Oberbeck-Boussinesqsystemnon-uniquenessweaksolutionsNavier-Stokesunstableprofileself-similarcoordinatescontractionmappinggravitationalfieldL1_tL2_xforcing
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 proves a non-uniqueness theorem for the Oberbeck-Boussinesq system with a Newtonian gravitational field on $\mathbb{R}^3$: there is a short time interval and external forces in $L^1_t L^2_x$ for which the system has infinitely many distributional solutions starting from rest. If correct, this shows that coupling the forced Navier-Stokes equations to a temperature equation and a gravitational term does not restore uniqueness in the weak class. The proof transplants an unstable mode of the linearized Navier-Stokes operator into self-similar coordinates, adds a non-steady temperature profile, and closes with a contraction mapping. A reader should care because the result makes weak solutions of thermal convection genuinely ambiguous, not merely formally underdetermined.

What carries the argument

The load-bearing object is the linearized temperature semigroup $e^{\tau L}$ with $L=\Delta+\tfrac12(1+\xi\cdot\nabla)+\bar U\cdot\nabla$ in self-similar variables. Lemma 3.1 claims it satisfies the heat-type smoothing estimate $\|e^{\tau L}\Theta_0\|_{H^k}\lesssim \max\{\tau^{-(k-m)/2},1\}\|\Theta_0\|_{H^m}$ for $k\ge m\ge0$, which is what lets the convection terms be absorbed in the contraction argument. The second ingredient is the unstable mode of the Navier-Stokes linearized operator $L_{ss}$: a smooth divergence-free profile with a positive eigenvalue and a parabolic semigroup bound, which seeds the velocity perturbation and forces the exponential weights in the solution spaces.

What would settle it

A direct check is to take a simple Gaussian profile $\bar U$ and initial temperature $\Theta_0$, write the linearized temperature equation in the paper's self-similar coordinates, and test numerically whether $\|e^{\tau L}\Theta_0\|_{H^k}$ grows at most like $\max\{\tau^{-(k-m)/2},1\}\|\Theta_0\|_{H^m}$ and whether the apparent drift term has the sign the paper assigns it. If the estimate fails for one such profile, Lemma 3.1 and the contraction step fail with it.

Watch

Extended reading notes

Core claim

Theorem 1.1 asserts that there exists $t_0>0$ and forces $f,h\in L^1_t L^2_x((0,t_0)\times\mathbb{R}^3)$ such that (1.1) has infinitely many distributional solutions $(u,\theta)\in L^2_t L^2_x((0,t_0)\times\mathbb{R}^3)$ with initial data $(u_0,\theta_0)=(0,0)$. The construction linearizes the self-similar form of the system around a compactly supported unstable profile $\bar U$ of the Navier-Stokes linearized operator, chooses the linear velocity perturbation as the real part of the unstable mode $e^{\lambda\tau}\rho$, sets a temperature profile $\bar\Theta(\xi,\tau)=e^{b\tau}\Theta(\xi)$, and defines the force and heat source by substitution. The remaining perturbation is found by a contraction map, and multiplying the unstable mode by different constants yields infinitely many distinct solutions.

Load-bearing premise

The proof rests on a single estimate: that the temperature equation keeps its heat-like smoothing property when the unstable background flow is added. The demonstration of that estimate in the paper is internally inconsistent — the drift term has conflicting signs, the claimed scaling of the background flow's derivatives does not match the coordinate change, and the parameter choices for the contraction step point the wrong way — so the theorem stands or falls on whether that estimate can be fixed.

Editorial extensions

If this is right

  • If the theorem is correct, the forced Oberbeck-Boussinesq system on $\mathbb{R}^3$ is non-unique in the weak class for zero initial data: the same force and heat source admit infinitely many distributional solutions.
  • The temperature field is not determined by the data in that class, so weak formulations of thermal convection need extra admissibility criteria to be predictive.
  • By the scaling property of the system, the construction gives examples on arbitrarily short time intervals, not just one fixed $t_0$.
  • The proof identifies a reusable mechanism: an unstable Navier-Stokes mode drives the linear velocity, and a heat-like semigroup estimate for the temperature equation with drift closes the contraction.
  • The result extends the known non-uniqueness of forced Navier-Stokes to a coupled system with gravity and thermal transport.

Reading between the lines

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

  • A natural next test is whether the same construction works for a bounded gravity vector such as $(0,0,1)$; the proof uses the Hardy inequality for $1/|\xi|$, and a bounded $\nabla G$ would need a different estimate for the temperature-to-velocity coupling.
  • The theorem does not address the unforced case $f=h=0$; here the forces are chosen to cancel the profile, so an unforced non-uniqueness, if possible, likely needs a new mechanism.
  • If the key temperature estimate is repaired, it could serve as a general tool for active-scalar systems with drift, making the whole strategy a template for transporting Navier-Stokes non-uniqueness into other coupled fluid models.
  • Because the solutions are built from an exponentially growing unstable mode, they probably violate the energy inequality; a testable conjecture is that energy-admissible weak solutions of this system remain unique even though the broader weak class is not.
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.

Referee Report

3 major / 6 minor

Summary. The paper constructs non-unique weak solutions to the three-dimensional Oberbeck-Boussinesq system with Newtonian gravitational field and zero initial data, with external force and heat source in L^1_t L^2_x. The construction follows the unstable-profile strategy of Albritton, Brué, and Colombo for forced Navier-Stokes: in self-similar variables the solution is split into an unstable profile, an unstable linear mode, and a perturbation, the force and heat source are solved from the profile, and the perturbation is obtained by a contraction mapping. The main novelty is to absorb the drift term \bar U·∇ into the heat semigroup for the temperature equation, leading to Lemma 3.1, which is then used to close the fixed-point argument for the temperature perturbation.

Significance. If Theorem 1.1 can be established, it would be a meaningful extension of the known non-uniqueness phenomenon for forced Navier-Stokes equations to the coupled Oberbeck-Boussinesq system with a Newtonian gravitational potential. The strategy of incorporating the drift term into the semigroup is natural, and the paper correctly identifies the temperature equation as the part where the standard heat-semigroup control is insufficient. The manuscript is not circular: it imports the unstable profile and semigroup estimates from [3] and constructs the forces explicitly. However, the key semigroup estimate for the operator actually used in the contraction is not established as written, so the current proof does not certify Theorem 1.1.

major comments (3)
  1. [§3.1, Eq. (3.3), and Lemma 3.1] The sign of the drift term in the operator L is inconsistent. The equation solved by Θ_p is ∂_τ Θ_p − ∆Θ_p − 1/2(1+ξ·∇)Θ_p + \bar U·∇Θ_p + (remaining nonlinear terms) = 0, so the relevant generator is ∆ + 1/2(1+ξ·∇) − \bar U·∇. The paper defines L = ∆ + 1/2(1+ξ·∇) + \bar U·∇, but the proof of Lemma 3.1 verifies the estimate for the solution of ∂_τ Θ − ∆Θ − 1/2(1+ξ·∇)Θ = −\bar U·∇Θ, which is the operator with the opposite drift sign. Thus the semigroup estimate proved in §3.2 does not apply to the operator named L in (3.3), and the Duhamel representation of Θ_p is not justified by the stated lemma.
  2. [§3.2, proof of Lemma 3.1] The scaling of the derivatives of \bar u is incorrect. For \bar u(x,t) = t^{-1/2}\bar U(x/√t), one has ‖D^k \bar u‖_{L∞} ≍ t^{-(k+1)/2}, not t^{(k+1)/2} as written. Consequently the displayed Leibniz-rule estimate involving the factor t^{-1} does not follow from the preceding line, and the estimate leading to (3.5) is not established. This is a separate defect from the sign issue above, and it concerns a load-bearing step in the proof of Lemma 3.1.
  3. [§3.4, contraction mapping] The statement that Φ maps B to itself and is a contraction when τ0 is sufficiently small is not supported. Under the parameter conditions (3.16), the exponents 2a−β, a+b−γ, −β+min{a+β,2β,γ}, and −γ+min{a+γ,β+b,β+γ} are positive, so the displayed factors e^{(2a−β)τ0}, e^{(a+b−γ)τ0}, and the corresponding M-factors are small only when τ0 is sufficiently negative with large magnitude, not when τ0 is close to zero. Since t0 = e^{τ0}, this means t0 should be chosen very small, but the argument must be rewritten to require a sufficiently negative τ0 and to verify that all estimates are uniform on (−∞, τ0].
minor comments (6)
  1. [§3.4, definition of Φ1] The displayed definition of Φ1 appears with an integral from τ to ∞, whereas the Duhamel formula (3.2) and the estimates (3.6)–(3.10) use integration from −∞ to τ. As printed, the map does not match the intended fixed-point equation; this should be corrected.
  2. [§3.4, Eq. (3.15)] The displayed estimate for the term Up·∇Θp is missing the factor ‖Up‖_X ‖Θp‖_Y; as written the line contains only the exponential factor and so cannot imply the stated contraction bound.
  3. [Theorem 1.1 and §3] Theorem 1.1 uses g for the heat source, while system (1.1) and the rest of the paper use h; this should be harmonized.
  4. [§3.2, proof of Lemma 3.1] The definition of φ_k(t) omits the volume element dx, and several integration-by-parts steps, including the term (−1)^k ∫ ∆^k θ · ∂_t θ dx, are not explained; these should be written out or justified.
  5. [§3.3, choice of \bar Θ] The paper should specify that the compactly supported profile Θ(ξ) in \bar Θ(ξ,τ)=e^{bτ}Θ(ξ) vanishes near ξ=0, or otherwise justify that the term \bar Θ ∇(1/|ξ|) belongs to L^2_x pointwise in time, since the Hardy-inequality step in §3.3 is only applied to Θ_p.
  6. [Throughout] There are several typographical errors and misspellings, for example 'sysem' in the abstract, 'Boussiness' in the introduction, 'Nest' on page 4, and 'τ′ is sufficiently small' in the last line of §3.4 instead of τ0; these should be corrected.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found; the derivation is self-contained given external Navier-Stokes instability input, and the internal proof gaps are correctness issues rather than circularity.

full rationale

The paper's central claim is not circular. The non-uniqueness result imports two external ingredients from Albritton, Brué and Colombo [3]: the unstable eigenfunction for the linearized Navier-Stokes operator (Lemma 2.1) and the associated semigroup regularization estimate (Lemma 2.3). These are cited theorems from other authors, not from this paper, and they do not assume the Oberbeck-Boussinesq conclusion. The forcing terms f and h are not fitted to the target multiplicity; they are defined explicitly in Section 3.1 by substituting the profiles Ubar and Thetabar into (1.1), and the claimed integrability is then checked directly. The perturbation equations (3.2)-(3.3) are solved by a contraction mapping in the spaces X and Y, and the non-uniqueness comes from the unstable mode Ul whose amplitude can be varied; no parameter is renamed as a prediction. The internal sign mismatch in Lemma 3.1 (the operator L is defined with a plus drift while the proof treats a minus drift), the scaling claim ||D^k ubar||_{L-infinity} ~ t^{(k+1)/2}, and the 'tau0 sufficiently small' versus 'sufficiently negative' discrepancies are genuine proof gaps, but they do not make any load-bearing step equivalent to its own input by construction. No self-citation is load-bearing, no fitted input is called a prediction, and no known result is merely renamed. Therefore no circularity step is identified.

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

No empirical fitting is involved. The numbers a, beta, gamma, b, tau0, M are proof parameters, not data-fitted constants; the only external input is the unstable profile from [3]. The paper introduces no new physical entities.

free parameters (3)
  • unstable eigenvalue a = a > 1
    Growth rate of the unstable Navier-Stokes mode taken from Lemma 2.1 and Remark 2.2. The proof requires a > 1 and a > delta; it is not derived in this paper.
  • growth exponents beta, gamma, b = satisfy the inequalities (3.11) and (3.16)
    These exponents define the weighted X and Y norms and are chosen by hand to make the integral estimates converge and the map contractive.
  • time horizon tau0 (with t0 = exp(tau0)) = negative with large magnitude
    The contraction factors depend on e^{c*tau0}; the proof text says sufficiently small but the argument requires -tau0 large, i.e., t0 very small.
assumptions (3)
  • domain assumption Unstable profile and linearized semigroup estimates from Albritton-Brue-Colombo (Lemmas 2.1 and 2.3) are valid.
    The entire construction is built on these imported results; they are not reproved.
  • ad hoc to paper The semigroup e^{tau L} for the temperature linearization satisfies Lemma 3.1.
    This is the new technical estimate; its proof in Section 3.2 contains sign and scaling inconsistencies, so as stated it is an unproven assumption.
  • standard math Hardy inequality bounds the L2 norm of (1/|xi|) grad Theta by a Sobolev norm of Theta.
    Used in Section 3.3 to control the gravity coupling term in the velocity equation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Non-uniqueness of Oberbeck-Boussinesq System with Gravitational Field." pith.science (2026). https://pith.science/paper/GVBWVQYN

@misc{pith2026241117929,
  author       = {Pith},
  title        = {Pith review of: Non-uniqueness of Oberbeck-Boussinesq System with Gravitational Field},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GVBWVQYN}},
  note         = {Machine review of arXiv:2411.17929}
}
read the original abstract

We consider the Oberbeck-Boussinesq system with gravitational force on the whole space. We prove the non-uniqueness of the system applying the unstable profile of Navier-Stokes equation.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 30 canonical work pages

  1. [3]

    Albritton, E

    D. Albritton, E. Bru´ e, and M. Colombo. Non-uniqueness o f Leray solutions of the forced Navier-Stokes equations. Annals of Mathematics , 196(1):415-455, 2022

    work page 2022

  2. [1]

    Abbatiello, M

    A. Abbatiello, M. Bul´ aˇ cek, and D. Lear. On the existenc e of solutions to generalized Navier- Stokes-Fourier system with dissipative heating. Meccanica (2024)

    work page 2024

  3. [2]

    Abbatiello and E

    A. Abbatiello and E. Feireisl. The Oberbeck-Boussinesq system with non-local boundary conditions. Quart. Appl. Math. 81(2),297–306(2023)

    work page 2023

  4. [4]

    Bella, E

    P. Bella, E. Feireisl and F. Oschmann. Rigorous derivati on of the Oberbeck-Boussinesq ap- proximation revealing unexpected term. Comm. Math. Phys. 403(3):1245–1273, 2023

    work page 2023

  5. [5]

    Boltzman

    B. Boltzman. Finite Amplitude Free Convection as an Init ial Value Problem-I. J.Atmos.Sci. 19,329(1962)

    work page 1962

  6. [6]

    Boussinesq

    J. Boussinesq. The´ eorie analytique de la chaleur (Gauthier-Villars, 1903)

    work page 1903

  7. [7]

    Large self-similar solutions to Oberbeck-Boussinesq system with Newtonian gravitational field

    L. Brandolese and G. Karch. Large self-similar solution s to Oberbeck-Boussinesq system with Newtonian gravitational field. arXiv : 2311.01093v1

    work page Pith review arXiv

  8. [8]

    Brandolese and M.E

    L. Brandolese and M.E. Schonbek. Large time decay and gro wth for solutions of a viscous Boussinesq system. Trans. Amer. Math. Soc. 364, 10 (2012), 5057–5090

    work page 2012

Show all 33 references
  1. [9]

    Bru´ e, M

    E. Bru´ e, M. Colombo and A. Kumar. Flexibility of Two-Dim ensional Euler Flows with Inte- grable Vorticity. arXiv : 2408.07934

    arXiv

  2. [10]

    Buckmaster and V

    T. Buckmaster and V. Vicol. Nonuniqueness of weak solut ions to the Navier-Stokes equation. Ann. of Math. (2) 189(1):101–144, 2019

    2019

  3. [11]

    Climent-Ezquerra, E

    B. Climent-Ezquerra, E. Ortega-Torres, M. Rojas-Meda r. On the regularity of weak solution for the Oberbeck-Boussinesq equations. arXiv : 2402.15207v1

    arXiv

  4. [12]

    Asymptot

    Raph¨ ael Danchin and Lingbing He, The Oberbeck-Boussi nesq approximation in critical spaces. Asymptot. Anal. 84(2013), no. 1-2, 61-102. MR 3134744

    2013

  5. [13]

    Danchin and M

    R. Danchin and M. Paicu. Existence and uniqueness resul ts for the Boussinesq system with data in Lorentz spaces. Physica D 237 (10-12) (2008), pp. 1444-1460

    2008

  6. [14]

    Diaz and G

    J.I. Diaz and G. Galiano. Existence and uniqueness of so lutions of the Boussinesq system with nonlinear thermal diffusion. Topological Methods Nonlinear Anal. 11(1): 59-82 (1998)

    1998

  7. [15]

    Feireisl and A

    E. Feireisl and A. Novotn´ y. Singular limits in thermod ynamics of viscous fluids. Advancesin Mathematical Fluid Mechanics. Birkh¨ auser/Springer, Cham, 2017. Second edition

    2017

  8. [16]

    Gao and H

    H. Gao and H. Liu. Global well-posedness and long time de cay of the 3D Boussinesq equations. J. Differ. Equ. 263 (2017), pp. 8649-8665

    2017

  9. [17]

    Grandi, A

    D. Grandi, A. Passerini. On the oberbeck-Boussinesq ap proximation for gases. International Journal of Nonlinear Mechanics 134: 103738, (2021)

    2021

  10. [18]

    Jia and V

    H. Jia and V. ˇSve´ rak. Local-in-space estimates near initial time for we ak solutions of the Navier-Stokes equations and forward self-similar solutio ns. Invent. Math. 196, 1(2014), 233- 265

    2014

  11. [19]

    Jia and V

    H. Jia and V. ˇSve´ rak. Are the incompressible 3d Navier-Stokes equation s locally ill-posed in the natural energy space? J. Funct. Anal. 268(12):3734-3766, 2015

    2015

  12. [20]

    Karch and N

    G. Karch and N. Prioux. Self-similarity in viscous Bous sinesq equations. Proc. Amer. Math. Soc. 136, 3 (2008), 879-888

    2008

  13. [21]

    Komo, Optimal initial value conditions for the exist ence of strong solutions of the Boussi- nesq equations

    C. Komo, Optimal initial value conditions for the exist ence of strong solutions of the Boussi- nesq equations. Ann Univ Ferrara 60 (2014), 377-396

    2014

  14. [22]

    B. Lai, C. Miao, and X. Zheng. Forward self-similar solu tions of the fractional Navier Stokes equations. Adv. Math. 352:981–1043, 2019

    2019

  15. [23]

    C.-C. Lai. Forward discretely self-similar solutions of the MHD equations and the viscoelastic Navier-Stokes equations with damping. J. Math. Fluid Mech. 21(3):Paper No. 38, 28, 2019

    2019

  16. [24]

    E.N. Lorenz. Deterministic Nonperiodic Flow. J.Atmos.Sci. 20,130(1963)

    1963

  17. [25]

    Oberbeck

    A. Oberbeck. Ueber die W¨ armleitung der Fl¨ ussigkeite n bei Ber¨ ucksichtigung der Str¨ umungeninfolge von Temperaturdifferenzen.Ann. Phys. Chem. 243(6), 271–292 (1879)

  18. [26]

    Rajagopal, G

    K.R. Rajagopal, G. Saccomandi, L. Vergori. On the appro ximation of isochoric motions of fluids under different flow conditions. Proceedings of the Royal Society A 471: 20150159, (2015)

    2015

  19. [27]

    V. Scheffer. An inviscid flow with compact support in spac e-time.J. Geom. Anal. 3 (1993), 343–401, 1993

    1993

  20. [28]

    Shnirelman

    A. Shnirelman. On the nonuniqueness of weak solution of the Euler equation. Comm. Pure Appl. Math. 50, 1997

    1997

  21. [29]

    Shnirelman

    A. Shnirelman. W eak solutions with decreasing energy o f incompressible Euler equa- tions.Comm. Math. Phys. 210, 541–603, 2000

    2000

  22. [30]

    M. Vishik. Instability and non-uniqueness in the Cauch y problem for the Euler equations of an ideal incompressible fluid. Part I. arXiv : 1805.09426. 10 XIAORAN LIU

    arXiv

  23. [31]

    M. Vishik. Instability and non-uniqueness in the Cauch y problem for the Euler equations of an ideal incompressible fluid. Part II. arXiv : 1805.09440

    arXiv

  24. [32]

    V. I. Yudovich. Non-stationary flows of an ideal incompr essible fluid. Z. Vychisl. Mat. i Mat. Fiz. 3, 1032–1066 (Russian), 1963

    1963

  25. [33]

    Zeytounian

    R.K. Zeytounian. Joseph Boussinesq and his approximat ion: a contemporary view. Compt. Rendus Mec. 331 (2003) 575-586. School of Mathematical Sciences, University of Science and T echnology of China, Hefei 230026, Anhui, China Email address : liuxiaoran@mail.ustc.edu.cn

    2003

Pith tools

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