REVIEW 3 major objections 3 minor 40 references
Finite-time singularity via multi-layer degenerate pendula for the 2D Boussinesq equation with uniform $C^{1,\sqrt{\frac{4}{3}}-1-\epsilon}\cap L^2$ force
T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read For every α below √(4/3)−1, smooth classical solutions of the forced 2D Boussinesq system blow up at t = 1 — the density-gradient integral diverges — while both forces keep uniform Hölder regularity up to the singular time.
desk verdict First whole-plane forced 2D Boussinesq blow-up with uniformly C^{1,alpha} force for all alpha below sqrt(4/3)-1; the construction is explicit and honest, but the load-bearing parameter feasibility in Section 6 needs referee scrutiny. 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
The load-bearing mechanism is a hierarchy of fluid layers, each prescribed by a stream function $\widetilde{\psi}^{(n)}_n(t,x) = B_n(t)\,\varphi(\lambda_n x_1)\varphi(\lambda_n x_2)\sin(x_1)\sin(x_2)$ written in coordinates that move and deform with the layer, so that the nonlinear transport is absorbed by a change of variables and the leading dynamics is read off from ODEs. The central identity is the degenerate-pendulum profile $\sin F(t) = 1/\cosh(t_{\max}-t)$: after rescaling time, the relative displacement of consecutive layer centers solves $\dot F = \sin F$ with the distinguished energy $E=1$, so the sine of the displacement is a single hyperbolic pulse that peaks exactly at the midpoint of the layer's life. This pulse synchronizes the growth of the vorticity amplitude — the growth of a layer is the integral of the density pulse, which is why the vorticity increase survives after the density is switched off (hysteresis) — and it also fixes the shape parameter $k_n(t)$ through $\int \cos F\,dt = \ln(\sin F(t)/\sin F(0))$. The exponents $\alpha^* = \sqrt{4/3}-1$ come out of the balance among the superexponential layer scalings $a_n = C^{(1-k_n)(1/(1-\gamma))^n}$, $b_n = C^{(1+k_n)(1/(1-\gamma))^n}$, the superexponentially small cutoffs $\lambda_n = C^{-\Lambda(1/(1-\gamma))^n}$ and the superexponentially short time steps $1-t_n \sim C^{-\delta(1/(1-\gamma))^n}$, chosen so that every error term stays subordinate to the leading term up to $t=1$.
What would settle it
Evaluate the parameter constraints of Proposition 10 and Section 6 numerically for a fine grid of $\alpha$ approaching $\sqrt{4/3}-1$: if the admissible region (say, the maximum allowed $k_{\max}$ as a function of $\alpha$) closes strictly below $\alpha^*$, the stated range of Hölder exponents collapses. Independently, simulate the two-layer toy model and compare the real displacement $\Xi^{(n)}(t)$ with the explicit pendulum profile $\sin(a_{n-1}(1)\Xi^{(n)}_0(t)) = 1/\cosh(\hat t_{\max}-\hat t)$: the error should shrink as $n$ grows, and any deviation that grows with $n$ would contradict Proposition 10.
Extended reading notes
Core claim
The paper's central claim is Theorem 1: for every $\alpha \in (0, \alpha^*)$ with $\alpha^* = \sqrt{4/3}-1$, there exist classical solutions $(u,\rho)$ of the forced Boussinesq system on $[0,1)\times\mathbb{R}^2$ such that $u$ and $\rho$ are smooth and compactly supported on every slab $[0,1-\varepsilon]$, the vorticity force satisfies $f_\omega \in C^0_t C^{\alpha}_{x,c}$ and the density force $f_\rho \in C^0_t C^{1,\alpha}_{x,c}$ uniformly up to $t=1$, and $\int_0^T \|\nabla\rho(t,\cdot)\|_{L^\infty}\,dt \to \infty$ as $T\to 1$. Since the divergence of that integral is the classical blow-up criterion for Boussinesq solutions, the singularity is genuine, and it occurs while both forces preserve their Hölder regularity at the blow-up time itself. The singular profile is odd-odd in the spatial variables, the forces are compactly supported, and the vorticity force is odd in $x_2$ so that the corresponding velocity force lies in $C^{1,\alpha} \cap L^2$.
Load-bearing premise
The load-bearing premise is that a non-empty parameter regime exists in which the real layer-center dynamics stays close to the degenerate-pendulum toy model on every interval where the error estimates are made, while all the intertwined force-regularity inequalities hold at once — a feasibility check deferred to the final parameter optimization in Section 6 rather than exhibited in full.
Editorial extensions
If this is right
- The integral of $\|\nabla \rho\|_{L^\infty}$ — the established blow-up criterion for the Boussinesq system — is shown to diverge for classical solutions whose external forces stay uniformly Hölder ($C^\alpha$ for the vorticity force, $C^{1,\alpha}$ for the density force) up to the singular time, so the singularity is driven by the solution itself, not by the regularity of the forcing.
- The blow-up happens in the well-posedness regime of the Boussinesq equations, so it is not an artifact of ill-posed initial data; the solution is smooth on every time slab $[0,1-\varepsilon]$ and only loses regularity as $t\to 1$.
- At the blow-up time the solution loses regularity by an amount $r_{\mathrm{loss}}$ that can be made arbitrarily small by choosing parameters, but with an upper bound that tends to $\alpha^*/2$ as $\alpha\to 0$ and to $0$ as $\alpha\to\alpha^*$.
- The blow-up rate is not well defined: there are time sequences accumulating at $t=1$ along which $\|\omega\|_{L^\infty}$ and $\|\partial\rho/\partial x_2\|_{L^\infty}$ grow like different powers of $1/(1-t)$, including instants where $\|\partial\rho/\partial x_2\|_{L^\infty} = 0$.
- The same layer construction, with the same degenerate pendulums, yields a finite-time blow-up for the axisymmetric 3D Euler equations, as the paper announces will appear in a forthcoming companion.
Reading between the lines
- The threshold $\sqrt{4/3}-1$ is derived from the parameter balance rather than from the equation's scaling; it would be informative to test numerically whether the admissible parameter region (for instance, the maximum allowed $k_{\max}$ as a function of $\alpha$) actually closes up exactly at $\alpha^*$, or whether a refined optimization could push the exponent higher.
- The hysteresis principle suggests a transferable design rule: any forcing protocol that delivers the same localized pendulum pulse to each successive layer with well-separated scales should produce the same finite-time blow-up, so the specific density profile (the choice of the activation function $h^{(n)}$) is likely not essential.
- Because the forces are compactly supported and the solution is classical away from the singular point, the construction offers a concrete target for numerical verification: the predicted layer displacement profile $\sin(a_{n-1}(1)\Xi^{(n)}_0(t)) = 1/\cosh(\hat t_{\max} - \hat t)$ is explicit enough to check in a two-layer simulation before testing the full cascade.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs compactly supported classical solutions of the forced inviscid 2D Boussinesq system on [0,1)×R^2 that develop a finite-time singularity at t=1, while the external forces retain Hölder regularity up to the singular time: the vorticity force belongs to C^0_t C^α_x and the density force to C^0_t C^{1,α}_x for every α<√(4/3)-1, with both forces compactly supported. The construction uses an infinite sequence of vorticity/density layers of decreasing spatial scale, introduced at times t_n accumulating at 1. The layer centers are designed to follow approximately the dynamics of a chain of degenerate inverted half-pendula, and the density is turned on and off in each layer so that the vorticity accumulates a hysteresis-like growth. The forces are defined as residuals of the layer equations, so their regularity is derived rather than assumed. The proof is organized as follows: Section 2 defines the stream function, density, and parameter choices; Section 3 proves that the actual layer separations Ξ^(n)(t) and the shape parameters k_n(t) stay close to the toy model; Sections 4 and 5 bound the residual forces; Section 6 performs the final parameter optimization and closes the proof of Theorem 1.
Significance. If the proof is completed, this would be a notable result: a pen-and-paper construction of a forced Boussinesq blow-up in the whole plane with forces that remain Hölder regular at the singular time, have compact support, and achieve the non-small Hölder range α<√(4/3)-1. The mechanism, based on an infinite chain of degenerate pendula and flickering density, is original and differs from the self-similar constructions of Elgindi–Pasqualotto and from the earlier forced constructions for 3D Euler and IPM. A strength of the paper is that the solution and the forces are explicitly constructed, with the forces obtained as residuals and no computer-assisted verification. The extension to axisymmetric 3D Euler announced in Remark 6, if realized, would increase the impact of the method. However, the current text leaves two load-bearing pieces incomplete: the proof of the uniform closeness estimate in Proposition 10 and the explicit parameter feasibility in Section 6. For this reason my assessment is conditional: the strategy is plausible and no internal contradiction is visible, but the main theorem is not yet fully established in the manuscript under review.
major comments (3)
- [§3.2, Proposition 10 and Eq. (108)] Proposition 10 is the load-bearing bridge between the actual layer separations Ξ^(n)(t) and the toy model Ξ_0^(n)(t). The proof of point 3 defines t*_n using inequality (108), 1/(π−a_{n-1}(1)Ξ^(n)(t)) ≤ K_n, and then asserts a uniform error 14C^{−ββ′min{β″,1−β″}δγ(1/(1−γ))^{n−1}} on [t*_n,1]. To close the argument one must show, with K_n = C^{ββ′β″δγ(1/(1−γ))^{n−1}} as proposed in the text, that the interval [t*_n,1] is nonempty and that this uniform error is small compared with the distance from a_{n-1}(1)Ξ_0^(n)(1) to π, which is C^{−k_max(1/(1−γ))^n} by Choices 12–14 and Remark 20. The manuscript does not verify these inequalities; it defers the admissible parameter set to Section 6. This is not a cosmetic omission, because the force bounds in Sections 4–5 are computed for the toy configurations and the entry conditions for layer n+1 are taken from Choice 13/Remark 20 on the basis of this closeness. I do not see a contradiction in the claimed estimates, but the parameter regime that makes them simultaneously true has not been exhibited.
- [§6, final parameter optimization] The theorem's quantitative range α∈(0,√(4/3)−1) depends on the existence of a nonempty admissible parameter set for C, γ, k_max, δ, Λ, Y, ε, ζ. Section 6 is presented as the place where the remaining parameters are optimized and the proof is closed, but the text provided stops short of displaying an explicit tuple of parameters or a chain of inequalities proving that all constraints coming from Proposition 9, Proposition 10, and the force bounds of Sections 4–5 are simultaneously satisfiable for every α below the stated threshold. In particular, the derivation of the critical exponent √(4/3)−1 from the parameter constraints is not shown. Since this exponent is the paper's headline improvement over previous forced constructions, the missing verification is load-bearing for Theorem 1.
- [§1.3, Remark 1] Remark 1 asserts that the local well-posedness argument of Chae–Kim–Nam extends to the forced case, but no proof is provided. If the theorem is meant to place the singularity 'in the well-posedness regime', the extension should be either proved or verified as a separate lemma; otherwise the claim in Remark 1 and the appeal to the blow-up criterion of [6] are not supported. The explicit construction itself may suffice to justify the existence of classical solutions on [0,1−ε], but the relationship between the constructed solutions and the asserted well-posedness theory should be made precise.
minor comments (3)
- [Choice 11, p. 55] The line 'we choose kn(1) = kn(1) = 0' appears to contain a typo; it should presumably read 'kn(1)=0'.
- [Abstract and Theorem 1] The abstract writes the Hölder exponent as C^{1,√(4/3)−1−ε}, while Theorem 1 states α∈(0,√(4/3)−1). It would be helpful to state explicitly that the ε in the abstract corresponds to the gap between α and the critical value.
- [§2.6, around Eq. (49)] The informal notation '∼' is used in several places to indicate order-of-magnitude equivalence without a precise definition. A short explanation of the convention would improve readability, especially for readers trying to follow the heuristic derivation of the toy model.
Circularity Check
No significant circularity: the construction is self-contained and the forces are defined as residuals, making their regularity an output of the estimates rather than an input.
full rationale
The paper explicitly prescribes the stream function and density (Choice 3 and Choice 8) and then defines the forces as the residuals of the layer equations (Proposition 1, equations (43)-(44)). The stated force regularity is therefore an output of the bounds in Sections 4-5, not a fitted or assumed input. The toy models for Ξ^(n)(t) and k_n(t) are derived in-paper from the ODEs in Choice 4/Proposition 6, and Proposition 10 proves the closeness of the real dynamics to these toy models by direct estimates that rely on Proposition 9, which is itself an independent induction. The citations to the authors' prior works [15], [17], [18] appear only as contextual 'siblings' in the introduction and comparison section; no load-bearing theorem, uniqueness statement, or ansatz is imported from them. The main caveat is that some parameter feasibility and the final optimization are deferred to Section 6, which is a completeness/correctness concern rather than a circularity: the paper does not assume the target force regularity or the blow-up rate in order to prove Proposition 10 or the force bounds. No equation is equivalent to its own input by construction, and no fitted parameter is renamed as a prediction.
Assumptions & free parameters
free parameters (8)
- C =
C > 2, taken sufficiently large (C >= Upsilon(...))
- gamma =
in (1/2,1), tends to 1 as alpha tends to alpha*
- k_max =
in [1/100,1), chosen close to 1
- delta =
> 0 small
- Lambda =
> 0, to be determined
- Y =
> 0, chosen so that t_1 = 0
- zeta =
in (0,1/4)
- epsilon =
in (0,1/4)
assumptions (4)
- standard math Hölder norm composition and product estimates (equations (22)-(25))
- standard math Existence and uniqueness for the pendulum ODE dF/dt = sin(F) with F(0) in [0,pi]
- domain assumption Blow-up criterion from Chae-Kim-Nam [6]: divergence of the integral of ||grad rho||_Linf implies singularity
- ad hoc to paper Extension of local well-posedness to the forced case (Remark 1)
Cite this review
Pith. "Pith review of Finite-time singularity via multi-layer degenerate pendula for the 2D Boussinesq equation with uniform $C^{1,\sqrt{\frac{4}{3}}-1-\epsilon}\cap L^2$ force." pith.science (2026). https://pith.science/paper/LXUHNVLP
@misc{pith2026250520988,
author = {Pith},
title = {Pith review of: Finite-time singularity via multi-layer degenerate pendula for the 2D Boussinesq equation with uniform $C^1,\sqrt\frac43-1-\epsilon\cap L^2$ force},
year = {2026},
howpublished = {\url{https://pith.science/paper/LXUHNVLP}},
note = {Machine review of arXiv:2505.20988}
}
abstract
We establish the existence of compactly supported solutions of the inviscid incompressible 2D Boussinesq equation with $C^{1,\sqrt{\frac{4}{3}}-1-\varepsilon}\cap L^{2}$ force that develop a singularity in finite time. Importantly, the force preserves this regularity at the blow-up time. Moreover, the forces in the vorticity and density equations have compact support. The mechanism behind the blow-up is an accumulated hysteresis effect on the vorticity caused by an infinite chain of "degenerate" pendula and flickering density.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[15]
D. C´ ordoba and L. Mart ´ ınez-Zoroa. Blow-up for the incompressible 3D-Euler equations with uniform C 1, 1 2 −ϵ ∩ L2 force. arXiv preprint arXiv:2309.08495, (2023)
arXiv 2023
-
[17]
D. C´ ordoba, L. Mart ´ ınez-Zoroa and F. Zheng. Finite time blow-up for the hypodissipative Navier Stokes equations with a force in L1 t C 1,ϵ x ∩ L∞ t L2 x. arXiv preprint arXiv:2407.06776, (2024)
arXiv 2024
-
[18]
D. C´ ordoba and L. Mart ´ ınez-Zoroa. Finite time singularities of smooth solutions for the 2D incompressible porous media (IPM) equation with a smooth source. arXiv preprint arXiv:2410.22920, (2024)
arXiv 2024
- [6]
-
[1]
J. Bedrossian, R. Bianchini, M.C. Zelati and M. Dolce. Nonlinear inviscid damping and shear-buoyancy instability in the two-dimensional Boussinesq equations. Comm. Pure. Appl. Math., 76:3585-3768, (2023)
work page 2023
-
[2]
R. Bianchini, L.E. Hientzsch, F. Iandoli. Strong ill-posedness in L∞ of the 2d Boussinesq equations in vorticity form and application to the 3d axisymmetric Euler Equations. SIAM Journal on Mathematical Analysis, 56(5), (2024)
work page 2024
-
[3]
J. Boussinesq. Th´ eorie Analytique de la Chaleur, volume Vol. II. Gauthier-Villars, Paris, (1903)
work page 1903
-
[4]
D. Chae and P. Constantin: Remarks on Type I Blow-Up for the 3D Euler Equations and the 2D Boussinesq Equations. Journal of Nonlinear Science, 31(77), (2021)
work page 2021
Show all 40 references
-
[5]
Chae and H.-S
D. Chae and H.-S. Nam. Local existence and blow-up criterion for the Boussinesq equations. Proc. Roy. Soc. Edinburgh Sect. A, 127(5):935–946, (1997)
1997
-
[7]
Chae and J
D. Chae and J. Wolf. On the regularity of solutions to the 2D Boussinesq equations satisfying type I conditions. Journal of Nonlinear Sciences. 29:643-654 (2019)
2019
-
[8]
J. Chen. Remarks on the smoothness of the C 1,α asymptotically self-similar singularity in the 3D Euler and 2D Boussinesq equations. Nonlinearity, 37(6), (2024)
2024
-
[9]
Chen and T
J. Chen and T. Hou. Finite time blowup of 2D Boussinesq and 3D Euler equations with C 1,α velocity and boundary. Communications in Mathematical Physics, 383(3):1559–1667, (2021)
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. arXiv:2210.07191, (2022)
2022 arXiv
-
[11]
Chen and T
J. Chen and T. Hou. Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data II: Rigorous numerics. Multiscale Modeling and Simulation, 23(1), (2025)
2025
-
[12]
K. Choi, T. Hou, A. Kiselev, G. Luo, V. Sverak, and Y. Yao. On the finite-time blowup of a one-dimensional model for the three-dimensional axisymmetric euler equations. Communications on Pure and Applied Math- ematics, 70(11):2218–2243, (2017)
2017
-
[13]
Constantin
P. Constantin. Analysis of hydrodynamic models. CBMS-NSF Regional Conf. Ser. in Appl. Math., 90 Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, (2017)
2017
-
[14]
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 axisymmetric incompressible 3D Euler equation with uniform C 1,α ∩ L2 force. To appear
-
[16]
C´ ordoba, L
D. C´ ordoba, L. Mart ´ ınez-Zoroa and F. Zheng. Finite time singularities to the 3D incompressible Euler equations for solutions in C ∞(R3 \ 0) ∩ C 1,α ∩ L2 . Annals of PDE, 11(19), (2025)
2025
-
[19]
R. Danchin. Remarks on the lifespan of the solutions to some models of incompressible fluid mechanics. Proc. Amer. Math. Soc., 141(6):1979–1993, (2013). 127
2013
-
[20]
E and C.-W
W. E and C.-W. Shu. Small-scale structures in Boussinesq convection. Phys. Fluids, 6(1):49–58, (1994)
1994
-
[21]
T.M. Elgindi. Finite-time singularity formation for C 1,α solutions to the incompressible Euler equations on R3. Ann. Math. (2) 194.3 (2021), pp. 647–727
2021
-
[22]
Elgindi, T.-E
T.M. Elgindi, T.-E. Ghoul and N. Masmoudi. On the Stability of Self-similar Blow-up for C 1,α Solutions to the Incompressible Euler Equations on R3. Camb. J. Math. 9 (2021), no. 4, 1035–1075
2021
-
[23]
Elgindi and I.-J
T.M. Elgindi and I.-J. Jeong. Finite-time singularity formation for strong solutions to the Boussinesq system. Annals of PDE 6,5 (2020)
2020
-
[24]
Elgindi and I.-J
T.M. Elgindi and I.-J. Jeong. Finite-time singularity formation for strong solutions to the axisymmetric 3D Euler equations. Annals of PDE 5.2 (2019): 1-51
2019
-
[25]
Elgindi and F
T.M. Elgindi and F. Pasqualotto. From Instability to Singularity Formation in Incompressible Fluids. arXiv:2310.19780, (2023)
2023 arXiv
-
[26]
Elgindi and F
T.M. Elgindi and F. Pasqualotto. Invertibility of a linearized Boussinesq flow: a symbolic approach. ArXiv preprint, (2023)
2023
-
[27]
Elgindi and K
T.M. Elgindi and K. Widmayer. Sharp decay estimates for an anisotropic linear semigroup and applications to the surface quasi-gestrophic and inviscid Boussinesq systems. SIAM J. Math. Anal. 47(6):4672-4684, (2015)
2015
-
[28]
Huang, X
D. Huang, X. Qin, X. Wang, D. Wei. On the exact self-similar blowup of the Hou-Luo model with smooth profiles. Available at SSRN: https://ssrn.com/abstract=4602516 or http://dx.doi.org/10. 2139/ssrn.4602516. (2023)
2023
-
[29]
Kiselev, J
A. Kiselev, J. Park, and Y. Yao. Small scale formation for the 2d Boussinesq equation. Analysis and PDE, 18(1), (2025)
2025
-
[30]
Jurja and K
C. Jurja and K. Widmayer. Long-time stability of a stably stratified rest state in the inviscid 2D Boussinesq equation. arXiv:2408.15154v1, (2024)
2024 arXiv
-
[31]
Luo and T
G. Luo and T. Hou. Toward the finite-time blowup of the 3D axisymmetric Euler equations: a numerical investigation. Multiscale Model. Simul., 12(4):1722–1776, (2014)
2014
-
[32]
A. Majda. Introduction to PDEs and waves for the atmosphere and ocean, volume 9 of Courant Lec- ture Notes in Mathematics. New York University, Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI, (2003)
2003
-
[33]
Majda and A.L
A.J. Majda and A.L. Bertozzi. Vorticity and Incompressible Flow. Cambridge Texts in Applied Mathemat- ics. Cambridge University Press, (2001)
2001
-
[34]
Masmoudi, B
N. Masmoudi, B. Said-Houari and W. Zhao. Stability of the Couette flow for a 2D Boussinesq system without thermal diffusivity. ARch. Ration. Mech. Anal., 245:645-752, (2022)
2022
-
[35]
Oberbeck
A. Oberbeck. ¨Uber die W¨ armeleitung der Fl¨ ussigkeiten bei der Ber¨ ucksichtigung der Str¨ omungen infolge von Temperaturdifferenzen. Annual Review of Physical Chemistry, 7((6)):271-292, (1879)
-
[36]
Pumir and E.D
A. Pumir and E.D. Siggia. Development of singular solutions to the axisymmetric Euler equations. Phys. Fluids A, 4(7):1472–1491, (1992)
1992
-
[37]
Rieutord
M. Rieutord. Fluid dynamics. An introduction. Grad. Texts Phys., Cham: Springer, (2015)
2015
-
[38]
Sarria and J
A. Sarria and J. Wu. Blowup in stagnation-point form solutions of the inviscid 2d Boussinesq equations. J. Differential Equations, 259(8):3559-3576, (2015)
2015
-
[39]
L. Tao, J. Wu, K. Zhao and X. Zheng. Stability near hydrostatic equilibrium to the 2D Boussinesq equations without thermal diffusion. Arch. Ration. Mech. Anal. 237:585-630, (2020)
2020
-
[40]
Y. Wang, C. Lai, J. G´ omez-Serrano and T. Buckmaster. Asymptotic self-similar blow up profile for 3-D Euler via physicsinformed neural networks, Phys. Rev. Lett. 130, (2023). 128
2023
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.