REVIEW 3 major objections 4 minor 1 cited by
Norm Inflation for Inviscid and Fully Dissipative Boussinesq Systems in Supercritical Spaces
T0 review · 3 major / 4 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read The 2D Boussinesq system is strongly ill-posed in nearly every supercritical Besov space: arbitrarily small smooth data can grow the density's norm past 1/ε in time <ε, while velocity stays bounded.
desk verdict Solid, likely correct norm-inflation construction for 2D Boussinesq in supercritical Besov spaces, with sharp thresholds; the dissipative stability step needs more detail before publication. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Central object: the approximate solution is a stationary radial vorticity $\bar w = \lambda^{2/p+1-\beta}(\log\log\lambda)^{-|\beta|/(2+|\beta|)} f(\lambda r)$, whose Biot–Savart velocity $\bar u$ is time-independent, purely angular, with angular frequency $\Omega_\lambda(r) = \lambda^{2/p-\beta}(\log\log\lambda)^{-|\beta|/(2+|\beta|)} u_\theta[f](\lambda r)/r$. The density is the initial profile advected by this flow; the profiles $f,g$ have disjoint supports and nondegenerate shear $|h'|\geq c_0$ on supp $g$, so at $t_* = \lambda^{-2/p-1+\beta}\log\log\lambda$ the shear creates or unwinds rapid radial oscillations and yields the Besov lower bound. A stability estimate transfers the growth to the exact solution; in the dissipative case the Laplacians are perturbative errors, shifting the range by two derivativ
What would settle it
Run a high-resolution numerical simulation of the inviscid Boussinesq system with the paper's compactly supported initial data and measure $\|\rho(t_*)\|_{\dot B^{\beta}_{p,\infty}}$ at $t_* = \lambda^{\beta - 2/p - 1} \log\log\lambda$: if the density norm fails to exceed $1/\epsilon$ for arbitrarily small $\epsilon$ while the velocity stays bounded, the stability transfer is wrong. A targeted analytic check is whether the bootstrap bound (3.31) closes for $\beta - 2/p$ just below $1$; if the Gronwall factor no longer tends to $0$, the transfer step is the point of failure.
Extended reading notes
Core claim
Inviscid case ($\mu=\nu=0$): norm inflation in $\dot B^{\beta}_{p,q}\times \dot B^{\beta}_{p,r}$ for $\beta\neq 0$, $1<p\leq\infty$, $-2<\beta-2/p<1$. Fully dissipative case ($\mu=\nu=1$): same conclusion for $-2<\beta-2/p<-1$. In both, for every $\epsilon>0$ there is a $C_c^\infty$ datum $(u_0,\rho_0)$ with $\|u_0\|_{\dot B^{\beta}_{p,1}}+\|\rho_0\|_{\dot B^{\beta}_{p,1}}<\epsilon$ and a time $0<t_*<\epsilon$ with $\|\rho(t_*)\|_{\dot B^{\beta}_{p,\infty}}>1/\epsilon$. Embedding $\dot B^{\beta}_{p,1}\hookrightarrow \dot B^{\beta}_{p,q}$ and $\|\cdot\|_{\dot B^{\beta}_{p,r}}\geq \|\cdot\|_{\dot B^{\beta}_{p,\infty}}$ give the result for every $q,r$. The velocity stays bounded in $\dot B^{\beta}_{p,\infty}$; only the density inflates.
Load-bearing premise
The load-bearing premise is that the exact Boussinesq solution with the paper's initial data stays within $o(1)$ of the explicit shear-transport approximate solution in the inflation norm up to the inflation time $t_*$; the Gronwall/bootstrap estimates that enforce this close only under the stated regime conditions, and the transfer also assumes smooth solutions exist and extend through $[0,t_*]$.
Editorial extensions
If this is right
- The solution map of the Boussinesq system is discontinuous at the origin in every Besov space covered: strong ill-posedness holds throughout the supercritical locally integrable range.
- For the inviscid system the upper threshold β−2/p<1 is sharp against the known local well-posedness at one derivative above, and for the dissipative system β−2/p<−1 matches the velocity-critical parabolic scaling.
- Norm inflation is driven purely by the density transport: even when the density's own scaling is subcritical (dissipative case), a velocity in the supercritical regime amplifies ρ above 1/ε while u stays small.
- The same construction yields norm inflation in H^s for the fractionally dissipative system in the range −1<s<2−max{α,γ}, s≠0.
Reading between the lines
- The paper's mechanism suggests a general criterion: any equation in which a passive scalar is advected by a flow with supercritical Besov regularity can exhibit density norm inflation, even when the scalar's own scaling looks subcritical; testing this on other stratified-flow models (anelastic, primitive) would show whether the phenomenon is generic.
- A quantitative next step is to determine whether the rate (log log λ)^{|β|/(2+|β|)} is optimal for these spaces; if a faster growth rate can be constructed, the current bound is not sharp.
- Because t*→0 as λ→∞, the construction implies instantaneous loss of regularity at t=0 in the supercritical regime; this could be made explicit by extracting a sequence of solutions with blow-up times tending to zero.
- The condition β−2/p>−2 confines the result to locally integrable data; the complementary distributional range is untouched by this construction and may require a genuinely different mechanism.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proves norm inflation (strong ill-posedness) for the 2D Boussinesq system in supercritical homogeneous Besov spaces. For the inviscid system, Theorem 1.1 asserts that for β≠0, 1<p≤∞, 1≤q,r≤∞, and −2<β−2/p<1, for every ε>0 there are C_c^∞ initial data with sum of the B^β_{p,1} norms <ε and a time t*<ε such that the density's B^β_{p,∞} norm exceeds 1/ε, while the velocity remains bounded. For the fully dissipative system, Theorem 1.2 asserts the same conclusion in the range −2<β−2/p<−1. The construction is explicit: a stationary radial vorticity creates an angular shear, and the density is transported by the associated velocity. For β>0 the shear creates rapid radial oscillations by time t*; for β<0 the initial density is premixed and then unwound. A perturbation (stability) argument transfers the growth from the approximate solution to the exact solution. The main λ-exponent bookkeeping appears consistent, with the admissible ranges exactly where the Gronwall factors and the final error terms decay.
Significance. If correct, the paper establishes a strong ill-posedness result in almost all supercritical locally integrable Besov spaces for the inviscid Boussinesq system and in a substantial supercritical range for the fully dissipative system, including negative regularity. This is a significant advance over the existing critical-space ill-posedness results and is directly in line with the recent Euler/Navier–Stokes norm-inflation program of Luo and others. The construction is explicit and largely parameter-free: the approximate solution is given in closed form, the smallness of the premixed negative-β data is proved by duality rather than assumed, and the lower bound at t* is computed directly. These are genuine strengths. The main weaknesses are two load-bearing technical points that are not fully written out: the endpoint case of the Besov duality statement used for p=∞, and the compressed stability proof in the fully dissipative case. Both are fixable, but they need to be addressed before the claims are fully verified.
major comments (3)
- [§2.1, Proposition 2.3] The stated Besov duality “for 1≤p,r≤∞” is not standard at p=∞ or r=∞. In the usual form, ∙B^s_{p,r}∙ is the dual of ∙B^{-s}_{p',r'}∙ only in the reflexive range; at p=∞ or r=∞ one must work with the appropriate “small Besov” completion. This proposition is used one-directionally at endpoints in Lemma 3.1 (q=∞) and Lemma 3.5 (p=∞), and those endpoint cases are part of Theorems 1.1 and 1.2. Please state the precise predual characterization, give a proof or a precise reference, and confirm that the estimates (3.10), (3.23)–(3.24) indeed yield the claimed Besov bounds when p=∞ or q=∞.
- [§4.1, Corollary 4.1] The stability transfer in the fully dissipative case is the most compressed step of the paper. The proof asserts that “using the standard commutator estimate gives inductively” the bound (4.15), but it does not display the commutator expansions, the pressure contributions, or the precise role of the bootstrap assumption (4.13). Since (4.12) is exactly what forces the admissible range β−2/p<−1, this is load-bearing. A complete induction with all product, commutator, pressure, and diffusion error terms should be supplied, at least for the first nontrivial derivative order.
- [§3.2, Lemma 3.6] In the inviscid stability lemma, the high-order L^2 induction and its closure are also only sketched. In particular, (3.35) is asserted rather than derived, and the final closure of the L^∞ bootstrap (3.31) relies on the Gagliardo–Nirenberg interpolation in Step 4 without displaying the interpolation parameters for all needed ranges of k and q. Since this lemma is the mechanism that transfers the approximate norm inflation to the exact solution, the missing details should be written out. The argument is plausible, but as written it is not fully verifiable.
minor comments (4)
- [§2.2, Lemma 2.6 / Appendix A] The construction of f and g is correct, but the verification that f^{(-2)} is compactly supported uses the convention f^{(-1)}(r)=∫_{-∞}^r f(s)ds; since f vanishes for r<1/2, this is fine but should be stated explicitly once for clarity.
- [§3.2, Lemma 3.6] The notation “Eλ≤exp(CL^2)=λ^{o(1)}” is used without tracking the constants C. It would help to state that C may depend on k and q but not on λ, and that all λ-powers are uniform.
- [§4.1, Corollary 4.1] At the end of the proof, the phrase “standard commutator estimates” should be accompanied by a precise reference to a textbook inequality (e.g., Kato–Ponce or the commutator estimates in [1]) for the Besov/Leibniz products used in both the inviscid and dissipative stability proofs.
- [§1, Remark 1.3] The claim that the argument also applies to the fractionally dissipative system and gives ill-posedness in H^s(R^2) for −1<s<2−max{α,γ}, s≠0, is not proved or even sketched. This is a minor issue if the remark is intended only as an outlook, but it should be labeled as such.
Circularity Check
No circularity: the construction is explicit and the stability transfer is independent of the inflation claim.
full rationale
The derivation chain is self-contained. The approximate solution is explicitly defined in (3.1)–(3.5); the radial profiles f and g are constructed in Lemma 2.6 from moment conditions and a shear nondegeneracy condition, with the proof in Appendix A. The smallness of the initial data is proved by direct duality and interpolation in Lemma 3.5, not assumed. The lower bound at t* is obtained by differentiating the explicit oscillatory profile and by interpolation, again with explicit estimates. The transfer to exact solutions is performed through the error systems (3.26) and (4.6), with zero initial error, and the estimates in Lemma 3.6 and Corollary 4.1 close by Gronwall arguments; no fitted parameter is renamed as a prediction. The only compressed point is Corollary 4.1's phrase 'standard commutator estimate' for the inductive higher-derivative bounds, but this is an omitted technical verification, not a circular step: the asserted bound (4.15) is not used as an input to itself nor derived from the norm-inflation conclusion. There is no load-bearing self-citation and no imported uniqueness theorem; the cited shear and cascade mechanisms are motivational, not the source of the theorem's content. Accordingly, the derivation does not reduce, by the paper's own equations, to the result it claims to prove.
Assumptions & free parameters
free parameters (2)
- λ (large frequency/size parameter) =
λ → ∞; chosen so that ‖u₀,ρ₀‖_{\dot B^β_{p,1}} < ε and ‖ρ(t*)‖_{\dot B^β_{p,∞}} > 1/ε
- radial profiles f, g =
explicit C_c^∞ bumps chosen in Lemma 2.6 and constructed in Appendix A
assumptions (4)
- standard math Besov-space toolbox: embeddings, interpolation, duality (Prop. 2.3), Bernstein inequalities, Fourier-multiplier boundedness (Lemmas 2.1–2.5)
- domain assumption Classical local well-posedness and continuation for smooth inviscid and dissipative Boussinesq solutions
- standard math Logarithmic-potential mean-value identity ∫_{−π}^{π} log|r−ρe^{iα}|dα = 2π log max{r,ρ}
- standard math Calderón–Zygmund boundedness of the Biot–Savart operator and zero-order Fourier multipliers on L^p (1<p<∞) and on homogeneous Besov spaces
Cite this review
Pith. "Pith review of Norm Inflation for Inviscid and Fully Dissipative Boussinesq Systems in Supercritical Spaces." pith.science (2026). https://pith.science/paper/WZQLFJ4J
@misc{pith2026260713694,
author = {Pith},
title = {Pith review of: Norm Inflation for Inviscid and Fully Dissipative Boussinesq Systems in Supercritical Spaces},
year = {2026},
howpublished = {\url{https://pith.science/paper/WZQLFJ4J}},
note = {Machine review of arXiv:2607.13694}
}
abstract
We prove norm inflation, in the sense of strong ill-posedness, for the two-dimensional Boussinesq system in supercritical Besov spaces. For the inviscid system, norm inflation holds in \(\dot B^\beta_{p,q}(\mathbb R^2)\times \dot B^\beta_{p,r}(\mathbb R^2)\) for \(\beta\neq0\), \(1<p\leq\infty\), \(1\leq q,r\leq\infty\), and \(-2<\beta-\frac{2}{p}<1\). For the fully dissipative system, the same conclusion holds in the range \(-2<\beta-\frac{2}{p}<-1\). In both cases, the results cover almost all supercritical Besov spaces satisfying the local integrability condition. Norm inflation occurs in the density component \(\rho\), while the velocity component \(u\) remains bounded. In the fully dissipative case, the inflation space is supercritical for \(u\), but subcritical for \(\rho\) with respect to its own scaling. This is not a contradiction: the density is transported by a velocity field in a supercritical regime, and this transport mechanism is precisely what produces norm inflation in \(\rho\).
Forward citations
Cited by 1 Pith paper
-
Strong ill-posedness for the MHD system in the supercritical regime: inviscid and viscous
Arbitrarily small smooth data for 3D incompressible MHD produce norm inflation in every supercritical Sobolev space (e.g., H^s, 0<s<5/2, for ideal MHD), with the velocity staying bounded in the ideal case.
Reference graph
Works this paper leans on
-
[1]
Bahouri, J
H. Bahouri, J. Chemin and R. Danchin,Fourier analysis and nonlinear partial dif- ferential equations, Grundlehren der mathematischen Wissenschaften, 343, Springer, Heidelberg, 2011. 33
2011
-
[2]
Bianchini, L
R. Bianchini, L. E. Hientzsch and F. Iandoli,Strong ill-posedness inL ∞ of the 2D Boussinesq equations in vorticity form and application to the 3D axisymmetric Euler equations, SIAM J. Math. Anal. 56 (2024), 5915–5968
2024
-
[3]
Bourgain and D
J. Bourgain and D. Li,Strong ill-posedness of the incompressible Euler equation in borderline Sobolev spaces, Invent. Math. 201 (2015), 97–157
2015
-
[4]
Bourgain and D
J. Bourgain and D. Li,Strong ill-posedness of the 3D incompressible Euler equation in borderline spaces, Int. Math. Res. Not. (2021), 12155–12264
2021
-
[5]
Bourgain and N
J. Bourgain and N. Pavlovi´ c,Ill-posedness of the Navier-Stokes equations in a critical space in 3D, J. Funct. Anal. 255 (2008), 2233–2247
2008
-
[6]
Chae,Local existence and blowup criterion for the Euler equations in the Besov spaces, Asymptotic Anal
D. Chae,Local existence and blowup criterion for the Euler equations in the Besov spaces, Asymptotic Anal. 38 (2004), 339–358
2004
-
[7]
Chae,Global regularity for the 2D Boussinesq equations with partial viscosity terms, Adv
D. Chae,Global regularity for the 2D Boussinesq equations with partial viscosity terms, Adv. Math. 203 (2006), 497–513
2006
-
[8]
Chae and H
D. Chae and H. Nam,Local existence and blow-up criterion for the Boussinesq equa- tions, Proc. R. Soc. Edinb., Sect. A, Math. 127 (1997), 935–946
1997
Show all 33 references
-
[9]
Chen and T
J. Chen and T. Hou,Finite time blowup of 2D Boussinesq and 3D Euler equations with C1,α velocity and boundary, Commun. Math. Phys. 383 (2021), 1559–1667
2021
-
[10]
Chen and T
J. Chen and T. Hou,Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data I: Analysis, Preprint, arXiv:2210.07191, 2022
2022 arXiv
-
[11]
C´ ordoba and L
D. C´ ordoba and L. Mart ´ ınez-Zoroa,Non existence and strong ill-posedness inCk and Sobolev spaces for SQG, Adv. Math. 407 (2022), Article ID 108570, 74 p
2022
-
[12]
C´ ordoba and L
D. C´ ordoba and L. Mart ´ ınez-Zoroa,Non-existence and strong ill-posedness inCk,β for the generalized surface quasi-geostrophic equation, Commun. Math. Phys. 405 (2024), Paper No. 170, 53 p
2024
-
[13]
C´ ordoba and L
D. C´ ordoba and L. Mart ´ ınez-Zoroa,Global unique solutions with instantaneous loss of regularity for SQG with fractional diffusion, Ann. PDE 10 (2024), No. 2, Paper No. 21, 52 p
2024
-
[14]
C´ ordoba, L
D. C´ ordoba, L. Mart ´ ınez-Zoroa and W. S. O˙ za´ nski,Instantaneous gap loss of Sobolev regularity for the 2D incompressible Euler equations, Duke Math. J. 173 (2024), 1931– 1971
2024
-
[15]
C´ ordoba, A
D. C´ ordoba, A. La ´ ın-Sanclemente and L. Mart ´ ınez-Zoroa,Finite-time singularity via multi-layer degenerate pendula for the 2D Boussinesq equation with uniform C1, q 4 3 −1−ε ∩L 2 force, Adv. Math. 480 (2025), Article ID 110480, 205 p
2025
-
[16]
Elgindi and I
T. Elgindi and I. Jeong,Ill-posedness for the incompressible Euler equations in critical Sobolev spaces, Ann. PDE 3 (2017), Paper No. 7, 19 p
2017
-
[17]
Elgindi and I
T. Elgindi and I. Jeong,Finite-time singularity formation for strong solutions to the Boussinesq system, Ann. PDE 6 (2020), Paper No. 5, 50 p. 34
2020
-
[18]
Elgindi and N
T. Elgindi and N. Masmoudi,L ∞ ill-posedness for a class of equations arising in hydrodynamics, Arch. Ration. Mech. Anal. 235 (2020), 1979–2025
2020
-
[19]
Guo and X
D. Guo and X. Luo,Norm inflation for the critical SQG equation, Int. Math. Res. Not. IMRN (2026), no. 11, rnag108
2026
-
[20]
Q. Jiu, C. Miao, J. Wu and Z. Zhang,The 2D incompressible Boussinesq equations with general critical dissipation, SIAM J. Math. Anal. 46 (2014), 3426–3454
2014
-
[21]
Kato and G
T. Kato and G. Ponce,Commutator estimates and the Euler and Navier-Stokes equa- tions, Comm. Pure Appl. Math. 41 (1988), 891–907
1988
-
[22]
Kim and I
J. Kim and I. Jeong,A simple ill-posedness proof for incompressible Euler equations in critical Sobolev spaces, J. Funct. Anal. 283 (2022), Article ID 109673, 34 p
2022
-
[23]
Li and W
Z. Li and W. Wang,Norm inflation for the Boussinesq system, Discrete Contin. Dyn. Syst., Ser. B 26 (2021), 5449–5463
2021
-
[24]
Luo,Ill-posedness of incompressible fluids in supercritical Sobolev spaces, arXiv preprint, arXiv:2404.07813, 2024
X. Luo,Ill-posedness of incompressible fluids in supercritical Sobolev spaces, arXiv preprint, arXiv:2404.07813, 2024
2024 arXiv
-
[25]
Luo,Sharp norm inflation for 3D Navier-Stokes equations in supercritical spaces, arXiv preprint, arXiv:2504.08288, 2025
X. Luo,Sharp norm inflation for 3D Navier-Stokes equations in supercritical spaces, arXiv preprint, arXiv:2504.08288, 2025
2025 arXiv
-
[26]
A. J. Majda,Introduction to PDEs and Waves for the Atmosphere and Ocean, Courant Lect. Notes Math. 9, AMS, New York, 2003
2003
-
[27]
A. J. Majda and A. L. Bertozzi,Vorticity and Incompressible Flow, Cambridge Uni- versity Press, Cambridge, 2002
2002
-
[28]
Misio lek and T
G. Misio lek and T. Yoneda,Local ill-posedness of the incompressible Euler equations inC 1 andB 1 ∞,1, Math. Ann. 364 (2016), 243–268
2016
-
[29]
Stefanov and J
A. Stefanov and J. Wu,A global regularity result for the 2D Boussinesq equations with critical dissipation, J. Anal. Math. 137 (2019), 269–290
2019
-
[30]
Stefanov, J
A. Stefanov, J. Wu, X. Xu and Z. Ye,Global regularity results of the 2D fractional Boussinesq equations, Math. Ann. 391 (2025), 5965–6012
2025
-
[31]
G. K. Vallis,Atmospheric and oceanic fluid dynamics. Fundamentals and large-scale circulation, 2nd edition. Cambridge: Cambridge University Press, 2017
2017
-
[32]
Wang,Ill-posedness for the Navier-Stokes equations in critical Besov spaces ˙B−1 ∞,q, Adv
B. Wang,Ill-posedness for the Navier-Stokes equations in critical Besov spaces ˙B−1 ∞,q, Adv. Math. 268 (2015), 350–372
2015
-
[33]
Yoneda,Ill-posedness of the 3D-Navier-Stokes equations in a generalized Besov space nearBM O−1, J
T. Yoneda,Ill-posedness of the 3D-Navier-Stokes equations in a generalized Besov space nearBM O−1, J. Funct. Anal. 258 (2010), 3376–3387. 35
2010
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.