A Note on Singularity Formation for a Nonlocal Transport Equation

Maria Radosz; Vu Hoang

REVIEW 2 major objections 3 minor 14 references

A Note on Singularity Formation for a Nonlocal Transport Equation

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 Odd unbounded solutions of the one-dimensional inviscid α-patch transport equation develop a singularity in finite time, ending in at least an odd cusp with profile sign(x)|x|^p.

desk verdict New finite-time singularity for the inviscid α-patch model, proved for a specific unbounded-data class; the abstract oversimplifies, but the core comparison argument is sound. read the letter →

arxiv 1908.09424 v2 pith:L44IW2XX submitted 2019-08-26 math.AP

Vu Hoang , Maria Radosz This is my paper · ORCID

classification math.AP MSC 35Q3535B44

keywords alpha-patchmodelnonlocaltransportequationfinite-timesingularityoddcuspbarriermethodactivescalarfractionalLaplacianunboundedinitialdata

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

This paper studies the one-dimensional inviscid α-patch model, a nonlocal transport equation in which the velocity is obtained from the transported quantity by a fractional integration law. The authors prove that for odd initial data that grow at least like $x^p$ at infinity, with $p=\gamma/2$ and $\gamma=1-\alpha$, the smooth solution cannot persist forever: its maximal lifetime is bounded by an explicit time $T(a_0)$ determined by a simple ODE. At that time the solution has developed at least an odd cusp, meaning $\omega(t,x)\sim \mathrm{sign}(x)|x|^p$ as $t\to T(a_0)$. This matters because one-dimensional active scalar models of this type are used as testing grounds for singularity and cusp formation in more realistic fluid equations such as surface quasi-geostrophic flow.

What carries the argument

The load-bearing object is the time-dependent barrier $\varphi(t,x)=a(t)^p\left((x/a(t))+1\right)^p-a(t)^p$, moving under the ODE $\dot a=-c_0 a^{1-p}$. The barrier is a strict subsolution: $\varphi_t+u[\varphi]\varphi_x<0$ for $x>0$, which follows from a ratio inequality $U(z)f'(z)/(-p f(z)+z f'(z))\ge c>0$ together with the chosen ODE. Because the equation is pure transport, comparison along particle trajectories lets the authors propagate $\omega>\varphi$ from the initial data all the way to the singular time. Since $\varphi(T(a_0),x)=x^p$ while $\omega$ is odd and $\omega(T(a_0),0)=0$, the inequality forces an infinite slope at the origin: at least an odd cusp.

What would settle it

Take any smooth odd solution satisfying all hypotheses of Theorem 2 and track the origin slope: compute $\limsup_{t\to T(a_0)}|\omega_x(t,0)|$. The theorem predicts this limit is infinite (at least an odd cusp); if for one such solution the slope remains finite through $T(a_0)$, the cusp claim is false. A complementary check: for compactly supported smooth initial data the inequality $\omega_0(x)>(1+\epsilon)\varphi(0,x)$ fails at large $x$, so testing localized data isolates whether the unbounded-growth assumption is essential.

Watch

Extended reading notes

Core claim

Theorem 2 states: set $p=\gamma/2$; there is a constant $c_0>0$ such that if $a(t)$ solves $\dot a=-c_0 a^{1-p}$ with $a(0)=a_0<1$, and if a smooth odd solution $\omega$ of $\omega_t+u[\omega]\omega_x=0$ starts above the barrier, $\omega(0,x)>(1+\epsilon)\varphi(0,x)$ for $x>0$ with $\epsilon$ satisfying (9) and finite weighted norms, then the solution's maximal lifetime $\bar T$ is at most $T(a_0)$. For all $t<\min\{\bar T,T(a_0)\}$, the comparison $\omega(t,x)>\varphi(t,x)$ holds; at $t=T(a_0)$ the barrier collapses to $x^p$, so if the solution has not broken down earlier it must have at least an odd cusp, $\omega\sim \mathrm{sign}(x)|x|^p$, or a stronger singularity such as a shock.

Load-bearing premise

The load-bearing premise is that the initial vorticity is unbounded and grows at least like $x^p$ at infinity—specifically $\omega_0(x)>(1+\epsilon)((x+a_0)^p-a_0^p)$ for all $x>0$, with the weighted norms finite—because every later comparison $\omega>\varphi$, and hence the forced infinite slope at the origin, traces back to that growth through the choice of $\epsilon$.

Editorial extensions

If this is right

The inviscid α-patch model has finite-time singularity formation for the whole class of odd, unbounded initial data described in Theorem 2; smooth solutions cannot be continued past $T(a_0)$.
The singular profile is pinned down from below: the solution must be at least an odd cusp with exponent $p=\gamma/2$, tied to the fractional order of the nonlocal velocity law, while stronger singularities such as shocks are not excluded.
The singularity time has an explicit upper bound set by the ODE $\dot a=-c_0 a^{1-p}$ and the initial scale $a_0$, giving quantitative control independent of the detailed shape of the initial data.
The comparison with $\varphi$ holds throughout the maximal smooth interval, connecting loss of smoothness directly to the collapse of the barrier rather than to estimates on $\omega$ alone.

Reading between the lines

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

Editorial: The theorem deliberately restricts to initial data that grow like $x^p$ at infinity; a natural testable extension is to determine whether data decaying faster than $x^p$ (for instance Schwartz-class data) can still form cusps by another mechanism, or whether the unbounded growth is essential to this singularity route.
Editorial: The barrier construction isolates a calculable property of the nonlocal kernel—positivity of $U$ and the ratio bound in Proposition 6—so the same proof scheme could transfer to other one-dimensional nonlocal transport models if a self-similar subsolution with the same inequality can be built.
Editorial: Since the proof produces the constant $c_0$ through compactness, numerically evaluating the ratio $U(z)f'(z)/(-p f(z)+z f'(z))$ could give sharper values of $c_0$ and hence quantitative lower bounds on the blowup time for given $a_0$.

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.

Desk Editor's Note

read the letter

Short version: this is a real result, not a hollow note. The authors construct a shrinking barrier φ(t,x)=a(t)^p((x/a(t)+1)^p -1) and show that any odd solution of the inviscid α-patch model that starts above it must stay above it until the barrier collapses. The collapse time is finite and the solution must develop at least a cusp of order x^p, p=γ/2. The barrier method is adapted from their ARMA 2017 paper, but the application to this equation with kernel |x|^{-γ} is new, and the exponent p=γ/2 is forced by scaling in Proposition 6, not fitted. The computations in Propositions 4 and 6 check out, and the comparison principle through first touching is coherent. The acknowledgements show the authors found and fixed gaps from the first version; that is honest work.

The soft spots are real but narrower than the abstract suggests. Theorem 2 requires initial data to be odd, smooth in the weighted norm, and strictly above φ(0,x) ~ x^p at infinity. That is a genuine unbounded-data class. For localized or Schwartz data, the theorem says nothing. The abstract's 'solutions of this model form singularities' is fair only if you read 'solutions' as 'solutions in this class.' I'd want the abstract to state the data condition. Second, the conclusion is a one-sided bound: at the singular time ω(T,x) ≥ x^p, not ω ~ sign(x)|x|^p as a two-sided asymptotic. The authors say 'at least a cusp (or potentially stronger singularity)' in the theorem, so they are precise in the statement; the abstract's 'characterization' overpromises. Both issues are presentation, not mathematics. The proof within the stated hypotheses is solid, and I did not find a load-bearing gap. The compactness step in Theorem 1 is sketched, but for a short note in this area that is acceptable.

Who is this for? Anyone working on blowup for active scalars, especially one-dimensional models and the α-patch family. It's a niche but useful increment, not a resolution of SQG or 2D Euler. For a referee, I'd send it out. The main theorem is clearly stated, the proof is mostly rigorous, and the remaining issues are wording and a few details a referee can ask to expand. With an abstract revision, it deserves publication in a good PDE journal.

Referee Report

2 major / 3 minor

Summary. The paper considers the one-dimensional inviscid α-patch transport equation ω_t + u[ω]ω_x = 0 with the nonlocal Biot–Savart law u[ω](x) = -∫_R |y-x|^{-(1-α)} ω(y) dy, and studies finite-time singularity formation for odd solutions. The main result, Theorem 2, constructs a self-similar barrier φ(t,x) = a(t)^p f(x/a(t)) with f(z) = (z+1)^p - 1 and p = γ/2, γ = 1-α. If a solves ȧ = -c_0 a^{1-p} with a(0) = a_0 < 1, and if the initial datum is odd, satisfies the weighted growth conditions in (5), and lies strictly above (1+ε)φ(0,x) for all x>0, then the solution remains above φ(t,x) up to time T(a_0); hence its maximal lifespan is at most T(a_0), so a singularity forms in finite time, and, if no earlier breakdown occurs, the solution is bounded below by an x^p cusp at the singular time. The proof uses weighted estimates for u[ω] and its derivatives, a regularized local-existence argument, and a comparison principle for the barrier.

Significance. The barrier argument is self-contained and has a genuine parameter-free feature: the exponent p = γ/2 is forced by the scaling in Proposition 6, and the constant c_0 is chosen below a universal constant rather than fitted to a particular solution. The comparison mechanism through Lemmas 1 and 2 and Proposition 5 is coherent, and the local-existence framework is standard. If the claims are restricted to the actual theorem hypotheses, this is a useful addition to the literature on singularity formation for nonlocal transport models, showing that a class of unbounded odd data can produce a cusp-type lower bound at breakdown. However, the paper's abstract and introduction claim more than the theorem proves: the theorem covers only initial data with at least x^p growth at infinity, and the profile statement is a one-sided lower bound, not the asymptotic cusp characterization advertised in the abstract.

major comments (2)

[Abstract and Section 2, Theorem 2] The abstract states without qualification that 'solutions of this model form singularities in finite time', but Theorem 2 applies only to odd initial data satisfying ω(0,x) > (1+ε)φ(0,x) for all x>0, where φ(0,x) = (x+a_0)^p - a_0^p. This hypothesis forces ω(0,x) ≥ const·x^p as x→∞, so the theorem says nothing about localized, Schwartz, or compactly supported data. The unbounded-growth condition is used essentially in Lemma 1 for x≥1: the proof uses particle origins X_0 ≥ x and requires (1+ε)φ(0,X_0) > (X_0+a_0)^p for all large X_0, an inequality that fails if ω_0 decays at infinity. The abstract and the introductory sentence 'We show that solutions of this model form singularities in finite time' must be qualified to this unbounded-data class; otherwise the stated generality is not established.
[Section 3.2.2, end of proof of Theorem 2] The paper claims in the abstract and in Section 1 a 'characterization of the solution profile at the singular time' and defines an odd cusp via the asymptotic equivalence ω ∼ sign(x)|x|^p as t→T_s. Theorem 2, however, proves only the lower bound ω(T(a_0),x) ≥ φ(T(a_0),x) = x^p for x>0, and the text itself says the solution forms 'at least a cusp (or a potentially stronger singularity)'. No matching upper bound is proved, so the exact asymptotic cusp profile advertised in the abstract is not obtained. The characterization language should be replaced by a statement that the singular profile dominates an x^p cusp, or a proof of the matching upper bound should be supplied.

minor comments (3)

[Section 3.2.1, Lemma 1] The derivation of the condition on ε from (9) is compressed: the proof states that (9) guarantees ε/(1+ε) > a_0^p/(1+a_0)^p, but the displayed inequality is not immediately equivalent. A one-line derivation would improve readability.
[Section 3.1] There is a typo in the sentence 'We first construct global solutions of of an approximate problem'; the duplicated 'of' should be removed.
[Section 3.1, regularized kernel] The definition k_ε(z) = η_ε^{-γ}(|z|) is typographically ambiguous; writing η_ε(|z|)^{-γ} or explicitly stating the exponent would avoid confusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the comparison proof is self-contained and the only self-citation is a non-load-bearing remark.

full rationale

The paper's derivation chain is self-contained and does not reduce any claimed prediction to fitted inputs or definitional identities. The barrier phi(t,x)=a(t)^p((x/a(t))+1)^p-a(t)^p is an explicit constructed function, and the comparison theorem is proven from the equation itself. Proposition 4 computes u[phi] by a direct change of variables, giving the exact scaling u[phi]=a^{1-gamma+p}U(x/a). Proposition 6 establishes a positive lower bound for the ratio U(z)f'(z)/(-p f(z)+z f'(z)) using the condition p=gamma/2; the constant c0 in Theorem 2 is any positive number below that universal bound, so it is not fitted to the solution omega or to any numerical blowup data. The assumption omega(0,x)>(1+epsilon)phi(0,x) is an external initial-data hypothesis, not a restatement of the conclusion. Lemma 1 uses that hypothesis and the sign of u to propagate the comparison from x>=1 and near x=0; Lemma 2 and Proposition 5 convert the strict subsolution condition (26), verified from Proposition 6, into a contradiction at any first contact time. The proof of finite-time singularity then follows structurally: if the smooth solution survived past T(a0), the comparison would force an infinite slope at x=0, contradicting smoothness. No step in this chain is equivalent to its input by construction. The only self-citation is [7], cited in the introduction with the remark 'Another result on cusp formation can be found in [7].' That citation is purely contextual and is not used as evidence for any estimate or theorem in this paper; all load-bearing inequalities are proved in Sections 3.1 and 3.2. The acknowledged limitations, such as the unbounded growth class for the initial data and the one-sided nature of the cusp lower bound, concern the scope of the theorem's applicability and its advertised strength, not circularity. They do not constitute a fitted-input-called-prediction or a self-citation chain. Accordingly, no circular step is present.

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

The paper's central claim rests on the model definition (a domain assumption), the odd-symmetry class of solutions, and standard weighted-norm and compactness arguments. There are no fitted free parameters and no invented physical entities; the barrier φ is a mathematical construction. The only potentially restrictive axiom is the unbounded-growth initial data class, which is explicit in Theorem 2 but omitted from the abstract.

assumptions (4)

domain assumption The 1D α-patch model with velocity u = -∫ |y-x|^{-γ}ω(y)dy, γ=1-α∈(0,1), and normalized constant is the object of study.
Defined in Section 2 before Eq. (4); the paper assumes this Biot-Savart law rather than deriving it. The physical relevance to 2D Euler/SQG is motivational.
domain assumption Solutions are odd and the velocity kernel K(x,y)=|y-x|^{-γ}-|x+y|^{-γ} is nonnegative for x,y>0.
Oddness is imposed on initial data after (5) and shown to propagate; nonnegativity of K is used for monotonicity of u[ω] in Lemmas 2 and Proposition 5.
standard math Weighted norm estimates (10) for u[ω] and its derivatives hold with universal constant C.
Proposition 1 states these; they are the technical backbone of the a priori bounds and are proved by direct integral estimates using 0<q<γ.
standard math Arzelà-Ascoli compactness in the weighted spaces recovers a C1 solution in the ε→0 limit of regularized problems.
Used at the end of Section 3.1 to complete the local existence proof; the details of the weighted-space compactness are not written out.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A Note on Singularity Formation for a Nonlocal Transport Equation." pith.science (2026). https://pith.science/paper/L44IW2XX

@misc{pith2026190809424,
  author       = {Pith},
  title        = {Pith review of: A Note on Singularity Formation for a Nonlocal Transport Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/L44IW2XX}},
  note         = {Machine review of arXiv:1908.09424}
}

abstract

The $\alpha$-patch model is used to study aspects of fluid equations. We show that solutions of this model form singularities in finite time and give a characterization of the solution profile at the singular time.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 12 canonical work pages

[1]

Choi, T.Y

K. Choi, T.Y. Hou, A. Kiselev, G. Luo, V. ’verák and Y. Yao, On the Finite-Time Blowup of a One-Dimensional Model for the Three-Dimensional Axisymmetric Euler Equations. CommunicationsonPureandAppliedMathematics,vol.70,no.11,Nov.2017,pp.2218–43. Scopus, doi:10.1002/cpa.21697

doi:10.1002/cpa.21697 2017
[2]

K.Choi,A.KiselevandY.Yao, Finitetimeblowupfora1dmodelof2dBoussinesqsystem. Comm. Math. Phys., 334(3):1667–1679, 2015

work page 2015
[3]

Constantin, P.D

P. Constantin, P.D. Lax and A. Majda, A simple one-dimensional model for the three- dimensional vorticity equation,Comm. Pure Appl. Math., 38:715–724, 1985

work page 1985
[4]

Constantin, A

P. Constantin, A. Majda and E. Tabak Formation of strong fronts in the 2-d quasigeostrophic thermal active scalar,Nonlinearity, 7(6):1495–1533, 1994

work page 1994
[5]

of Math.(2), 162(3):1377–1389, 2005

A.Córdoba,D.CórdobaandM.A.Fontelos,Formationofsingularitiesforatransportequation with nonlocal velocity,Ann. of Math.(2), 162(3):1377–1389, 2005

work page 2005
[6]

On a one-dimensionalα-patch model with nonlocal drift and fractional dissipation.Trans

Hongjie Dong and Dong Li. On a one-dimensionalα-patch model with nonlocal drift and fractional dissipation.Trans. Amer. Math. Soc., 366(4):2041–2061, 2014

work page 2014
[7]

Hoang, M

V. Hoang, M. Radosz, Cusp formation for a nonlocal evolution equation, Archive for Rational Mechanics and Analysis, June 2017, Volume 224, Issue 3, pp 1021–1036, https://doi.org/10.1007/s00205-017-1094-3

work page doi:10.1007/s00205-017-1094-3 2017
[8]

Hoang, B

V. Hoang, B. Orcan-Ekmeckci, M. Radosz and H. Yang, Blowup with vorticity control for a 2D model of the Boussinesq equations. Journal of Diﬀerential Equations, Vol. 264 (12), p.7328–7356 (2018)

work page 2018

Show all 14 references

[9]

T. Y. Hou and, G. Luo, On the ﬁnite-time blowup of a 1D model for the 3D incompressible Euler equations,arXiv:1311.2613

arXiv
[10]

Hou and G

T.Y. Hou and G. Luo, Potentially singular solutions of the 3D axisymmetric Euler equations. PNAS, vol. 111 no. 36, 12968–12973,DOI 10.1073/pnas.1405238111

doi
[11]

Hou and C

T.Y. Hou and C. Li, Dynamic stability of the three-dimensional axisymmetric Navier-Stokes equations with swirl, Commun. Pure Appl. Math, vol. LXI, 661-697 (2008)

2008
[12]

Kiselev and V

A. Kiselev and V. ’verák. Small scale creation for solutions of the incompressible two dimensional Euler equation.Ann. of Math.(2), 180(3):1205–1220, 2014

2014
[13]

Madja, A

A. Madja, A. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press, 2002

2002
[14]

Pedlosky, Geophysical ﬂuid dynamics, Springer Verlag, New York (1979)

J. Pedlosky, Geophysical ﬂuid dynamics, Springer Verlag, New York (1979)

1979

Pith tools

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

[1] [1]

Choi, T.Y

K. Choi, T.Y. Hou, A. Kiselev, G. Luo, V. ’verák and Y. Yao, On the Finite-Time Blowup of a One-Dimensional Model for the Three-Dimensional Axisymmetric Euler Equations. CommunicationsonPureandAppliedMathematics,vol.70,no.11,Nov.2017,pp.2218–43. Scopus, doi:10.1002/cpa.21697

doi:10.1002/cpa.21697 2017

[2] [2]

K.Choi,A.KiselevandY.Yao, Finitetimeblowupfora1dmodelof2dBoussinesqsystem. Comm. Math. Phys., 334(3):1667–1679, 2015

work page 2015

[3] [3]

Constantin, P.D

P. Constantin, P.D. Lax and A. Majda, A simple one-dimensional model for the three- dimensional vorticity equation,Comm. Pure Appl. Math., 38:715–724, 1985

work page 1985

[4] [4]

Constantin, A

P. Constantin, A. Majda and E. Tabak Formation of strong fronts in the 2-d quasigeostrophic thermal active scalar,Nonlinearity, 7(6):1495–1533, 1994

work page 1994

[5] [5]

of Math.(2), 162(3):1377–1389, 2005

A.Córdoba,D.CórdobaandM.A.Fontelos,Formationofsingularitiesforatransportequation with nonlocal velocity,Ann. of Math.(2), 162(3):1377–1389, 2005

work page 2005

[6] [6]

On a one-dimensionalα-patch model with nonlocal drift and fractional dissipation.Trans

Hongjie Dong and Dong Li. On a one-dimensionalα-patch model with nonlocal drift and fractional dissipation.Trans. Amer. Math. Soc., 366(4):2041–2061, 2014

work page 2014

[7] [7]

Hoang, M

V. Hoang, M. Radosz, Cusp formation for a nonlocal evolution equation, Archive for Rational Mechanics and Analysis, June 2017, Volume 224, Issue 3, pp 1021–1036, https://doi.org/10.1007/s00205-017-1094-3

work page doi:10.1007/s00205-017-1094-3 2017

[8] [8]

Hoang, B

V. Hoang, B. Orcan-Ekmeckci, M. Radosz and H. Yang, Blowup with vorticity control for a 2D model of the Boussinesq equations. Journal of Diﬀerential Equations, Vol. 264 (12), p.7328–7356 (2018)

work page 2018