Axisymmetric self-similar solutions to the MHD equations without magnetic diffusion

Shaoheng Zhang

REVIEW 2 major objections 6 minor 11 references

Axisymmetric self-similar solutions to the MHD equations without magnetic diffusion

T0 review · 2 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read Smooth axisymmetric self-similar solutions of the MHD equations without magnetic diffusion, with a purely azimuthal magnetic field, are exactly the Landau solutions with zero magnetic field; in the half-space the only such solution is zero.

desk verdict Plausible extension of the Landau classification to ideal MHD with a swirl-only magnetic field; the proof has a repairable gap in the Riccati step that must be fixed before acceptance. read the letter →

arxiv 2506.20131 v1 pith:PIRDG2OB submitted 2025-06-25 math.AP

Shaoheng Zhang This is my paper · ORCID

classification math.AP MSC 35Q3576W05

keywords magnetohydrodynamicsself-similarsolutionLandauno-slipboundaryconditionNavier-slipaxisymmetricmagneticdiffusionpointsingularity

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 studies steady, incompressible magnetohydrodynamic flows with no magnetic diffusion, looking for solutions that are both axisymmetric and scale-invariant (self-similar), with a magnetic field pointing only in the azimuthal direction. Its first theorem states that in $\mathbb{R}^3\setminus\{0\}$, the only such smooth solutions are the classical Landau jets of the Navier-Stokes equations paired with zero magnetic field; the magnetic field cannot survive the self-similar ansatz. The second theorem states that in the half-space, with either no-slip or Navier-slip boundary conditions, the only such solution is $(\mathbf{u},\mathbf{B})=(0,0)$. A sympathetic reader would care because this gives a complete rigidity statement: within this symmetry class, the diffusion-free MHD system has no new point-singular steady solutions beyond those already known for Navier-Stokes, and the magnetized swirl component is always killed by the equations and axis regularity.

What carries the argument

The working object is the angular-profile ODE system obtained by substituting the self-similar axisymmetric ansatz into (1), reducing the three-dimensional PDE to a one-dimensional boundary-value problem on $\varphi\in(0,\pi)$ (or $(0,\pi/2)$). The two identities that carry the rigidity are $H'=gH$ for $H=(h\sin\varphi)'/\sin\varphi$, which forces $h\equiv0$, and $(Bg^2\sin\varphi)'=0$, which, together with $g(0)=B(0)=0$, forces $Bg^2=0$ and then kills $B$ on each maximal interval where it is nonzero via the pressure relations $B^2-2P=0$ and $P'+B^2\cot\varphi=0$. Once those components vanish, the remaining ODE for the angular profiles is the Landau equation $(1-t^2)L''+2L+LL'=0$ with $L=K(\varphi)$, and its explicit solution $L(t)=2(1-t^2)/(t-a)$, $|a|>1$, yields the Landau jet, the classical explicit $(-1)$-homogeneous Navier-Stokes solution.

What would settle it

Solve the angular ODE system (6)-(10) on $(0,\pi)$ with boundary conditions (11) and a nonzero $B$ profile supported in a proper subinterval $(\alpha,\beta)$: the theorem predicts that the endpoint relation $P(\alpha)=B(\alpha)^2/2$ cannot coexist with $P=C/\sin^2\varphi$, so a numerical or symbolic solution with $B$ not identically zero would refute Theorem 1.1. In the half-space, the same check on $(0,\pi/2)$ with the boundary data (21)-(23) would test Theorem 1.2.

Watch

Extended reading notes

Core claim

The central discovery is Theorem 1.1: every smooth axisymmetric self-similar solution $(\mathbf{u},\mathbf{B})$ of (1) in $\mathbb{R}^3\setminus\{0\}$ with $\mathbf{B}=B(\varphi)/\rho\,\mathbf{e}_\theta$ has $\mathbf{u}$ a Landau solution and $\mathbf{B}=0$. Theorem 1.2 is the companion statement in $\mathbb{R}^3_+$ under no-slip or Navier-slip conditions, where the conclusion is $\mathbf{u}=\mathbf{B}=0$. The proof writes the unknown as $\mathbf{u}=f(\varphi)/\rho\,\mathbf{e}_\rho+g(\varphi)/\rho\,\mathbf{e}_\varphi+h(\varphi)/\rho\,\mathbf{e}_\theta$, $\mathbf{B}=B(\varphi)/\rho\,\mathbf{e}_\theta$, $p=P(\varphi)/\rho^2$, reduces (1) to an ODE system in $\varphi$, then shows $h\equiv 0$ by combining $H'=gH$ with the identity $\int_0^\pi H\sin\varphi\,d\varphi=0$; then shows $B\equiv0$ by a maximal-interval contradiction using $g=0$, $B^2-2P=0$, and $P'+B^2\cot\varphi=0$ with endpoint values $P(\alpha)=B(\alpha)^2/2$. With $B=h=0$, the remaining equations are exactly the Landau ODE, and the explicit solution gives the Landau jet. In the half-space, after the same steps the remaining Navier-Stokes system is known to have only the zero solution under the stated boundary conditions.

Load-bearing premise

The load-bearing premise is that the angular profiles are smooth up to the symmetry axis, so that $B(0)=B(\pi)=0$ and $P(\alpha)=B(\alpha)^2/2$ hold at the endpoints where $B$ vanishes; if the axis had unaccounted singular behaviour, or (for the half-space) if the known vanishing of axisymmetric self-similar Navier-Stokes solutions failed, the conclusion would not follow.

Editorial extensions

If this is right

Within the swirl-only magnetic ansatz, the diffusion-free MHD system has the same one-point singular steady solutions as Navier-Stokes: every such solution is a Landau jet with zero magnetic field.
In the half-space, no nontrivial axisymmetric self-similar steady solution exists under no-slip or Navier-slip conditions, so self-similar jets or blow-up profiles are excluded in that geometry.
The coupling term $(\mathbf{B}\cdot\nabla)\mathbf{B}$ cannot sustain a nonzero azimuthal magnetic field in a self-similar axisymmetric setting; the magnetic field is killed by the angular ODE structure even though the system is nonlinear.
The characterization is sharp: any Landau jet paired with $\mathbf{B}=0$ is already a solution of the MHD system, so the theorem identifies the full solution set, not just a necessary condition.

Reading between the lines

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

A natural test of the mechanism is to relax the smoothness-at-the-axis assumption: if profiles with $B(0)\neq 0$ or endpoint values other than $P(\alpha)=B(\alpha)^2/2$ can be constructed in a weaker regularity class, the rigidity would likely fail; the paper does not address that class.
The same angular ODE structure may extend to other symmetry reductions of the MHD equations, such as magnetic fields with poloidal components or flows in exterior domains, but those extensions are not established here.
One could probe Theorem 1.1 numerically by discretizing the angular ODE system and searching for nonzero $B$ profiles; the theorem predicts every such search converges to $B=0$, so a converged nonzero profile would pinpoint where the endpoint argument breaks.

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 plausible and useful extension of the Landau classification to the ideal MHD system with a swirl-only magnetic field, but the proof as written has a real gap in the Riccati reduction. I think it is fixable and the paper deserves a serious referee.

The genuinely new content is Theorem 1.1: under B = B(φ)/ρ e_θ, the only axisymmetric self-similar solutions in R^3\{0} are the Landau jets with no magnetic field. Theorem 1.2 extends this to the half-space. Both remove the smallness condition that appeared in the earlier diffusive result [ZWW25] and drop magnetic diffusion from the system. The ODE reduction is standard (following Tsai's book), and the appendix derivation of the profile equations via the Riemannian connection is careful. The arguments that h ≡ 0 (via H' = gH and an integral zero) and that Bg^2 = 0 are correct.

Where I have a problem: the step after Eq. (15). The ODE (1-t^2)L'' + 2L + LL' = 0 with L = (1-t^2)w does not immediately give w' + w^2/2 = 0. Direct computation shows q := w' + w^2/2 satisfies (1-t^2)q' - 4t q = 0, so q = C/(1-t^2)^2. To conclude q = 0 you need an endpoint argument: axis regularity gives g(φ)=O(sin φ), hence L(t) = O(1-t^2) near t=±1, so w = L/(1-t^2) is bounded; then C must vanish. The paper skips this, and the skipped argument is load-bearing for the uniqueness claim. I believe it is repairable, but it needs to be written.

Two smaller notes. In the maximal-interval argument for B, the identity P = C_1/sin^2φ and the conclusion C_1 = 0 from P(α) = 0 are fine when α ∈ (0,π), but in the half-space proof one endpoint can be α = 0, a coordinate singularity; this needs a short limit argument rather than an assertion. And the use of [KMT18, Theorem 5.1] at the end of Theorem 1.2 is legitimate, but the reduction to the Navier-Stokes case should state explicitly that the cited theorem applies after B and h vanish.

None of this shakes my view that the result is very likely true. The paper is worth refereeing; I would ask for the missing integration step and the endpoint regularity before accepting.

Referee Report

2 major / 6 minor

Summary. The paper studies (-1)-homogeneous axisymmetric self-similar solutions of the stationary incompressible MHD equations without magnetic diffusion, assuming the magnetic field has only the swirl component B = B(φ)/ρ e_θ. Theorem 1.1 asserts that every smooth such solution in R^3\{0} has u equal to a Landau solution and B ≡ 0. Theorem 1.2 asserts the analogous triviality in the half-space R^3_+ under either no-slip or Navier-slip boundary conditions: u ≡ 0 and B ≡ 0. The proof strategy is to reduce the PDE system to a system of ODEs for the angular profiles f, g, h, B, P (equations (6)-(10)), first proving h ≡ 0 from the identity for H = (h sin φ)'/sin φ, then proving B ≡ 0 from the conservation law (B g^2 sin φ)' = 0 together with an interval argument, and finally deriving the Landau profile from a scalar ODE for L(t). The manuscript also contains an alternative derivation of the ODE system via the Riemannian connection in Appendix A.

Significance. If the classification is correct, the paper gives a clean rigidity result: the only one-point singular, axisymmetric, self-similar steady states of the diffusion-free MHD system with swirl-only magnetic field are the classical Landau jets with zero magnetic field, and the half-space analogue is trivial. A particular strength is that the ODE reduction is explicit and parameter-free, and the proof does not use fitted parameters or the author's own prior results. However, one load-bearing step in the derivation of the Landau profile is not justified as written, and the endpoint treatment of the interval argument for B ≡ 0 has a gap. These issues are localized and repairable, but they currently prevent the central uniqueness claim from being fully established.

major comments (2)

[Proof of Theorem 1.1, after Eq. (15)] The Riccati reduction is not derived correctly. Substituting L(t) = (1-t^2)w(t) into (1-t^2)L'' + 2L + LL' = 0 does not directly give w' + w^2/2 = 0; direct algebra gives (1-t^2)q' - 4t q = 0 for q := w' + w^2/2, whose general solution is q = C/(1-t^2)^2. To conclude C = 0 one needs an endpoint argument: the boundary conditions g(0) = g(π) = 0 and smoothness of g imply L(t) = O(1-t) as t→1 and L(t) = O(1+t) as t→-1, so w = L/(1-t^2) is bounded at the endpoints. A nonzero C would make q nonintegrable at t = ±1, forcing w to be unbounded, a contradiction; but this argument is absent from the manuscript. Without it, the possibility of additional non-Landau solutions with the same boundary data is not excluded, so this is a load-bearing gap for Theorem 1.1.
[Proof of Theorem 1.1, endpoint case in the B ≡ 0 argument] In the maximal-interval argument proving B ≡ 0, the interval (α,β) may have α = 0 or β = π, since B(0) = B(π) = 0 but B may be nonzero arbitrarily close to the poles. For such an endpoint, the identity P(α) = B(α)^2/2 = 0 is not available: P(0) need not vanish for a smooth homogeneous solution, and indeed the Landau pressure has P(0) ≠ 0. Thus the conclusion C_1 = 0 from P(α) = 0 is unjustified when the interval touches φ = 0 or φ = π. The gap is repairable by noting that the local formulas P = C_1/sin^2 φ (equivalently B = C/sin φ) are incompatible with smoothness of B and P at the axis unless C_1 = 0, but the manuscript does not supply this argument. The same endpoint issue occurs in the proof of Theorem 1.2 when the maximal interval has α = 0.

minor comments (6)

[Proof of Theorem 1.1, after Eq. (15)] The boundary condition L(-1) = L(1) = 0 is stated without explanation; it follows from K(0) = K(π) = 0 via K(φ) = g(φ) sin φ and should be derived explicitly.
[Equations (12) and (19)] The notation Bg^2 = 0 is potentially ambiguous; it should be written as B(φ) g(φ)^2 = 0 or accompanied by a sentence clarifying that it is a product and not composition.
[Proof of Theorem 1.2, final step] The application of [KMT18, Theorem 5.1] is the sole external input for the conclusion u = 0 in the half-space case; the hypotheses of that theorem should be stated explicitly so that the reader can verify that the reduced Navier-Stokes solution obtained here satisfies them.
[Introduction] There is a typo in the sentence about the result of Sverak: 'axissymmetry' should be 'axisymmetry'.
[Section 2.2] In the formula for β, the expression '2a2 log' should be typeset as 2a^2 log to match the standard Landau formula; please check the transcription.
[References] The reference [ZWW25] is an arXiv preprint; if a journal version exists, it should be cited instead of or in addition to the preprint.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; the derivation is self-contained and relies only on standard external results.

full rationale

The proof of Theorem 1.1 reduces the MHD system to a set of ODEs for the profile functions f, g, h, B, P, then solves those ODEs explicitly. The assumption B = B(φ)/ρ eθ is a symmetry ansatz, not the conclusion; the paper then proves B ≡ 0 and recovers the Landau solution. No fitted parameters are introduced, no target result is assumed, and no author's own prior theorem is used as a load-bearing input. The only cited external results are the standard Landau solution formulas from [Lan44, Tsa18, CK04, TX98] and the half-space Navier-Stokes classification [KMT18, Theorem 5.1], both of which are independent benchmark results and are cited for classification, not for the derivation itself. The author's prior work [ZWW25] appears only in the introduction as background context and is not used in either proof. The skeptical note about the Riccati step — that q = w' + w²/2 satisfies (1−t²)q' − 4tq = 0, so q = 0 needs an endpoint or regularity argument — is a mathematical rigor concern, not a circularity concern; a gap in an explicit computation does not make the claim equivalent to its input by definition. Therefore the circularity score is 0.

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

No parameters are fitted to data; the Landau parameter a simply labels which Landau solution appears. No new physical entities are introduced. The analysis relies on standard ODE uniqueness, standard Riemannian geometry identities, the self-similar ansatz, and one external classification theorem.

assumptions (5)

standard math Standard ODE uniqueness for first-order equations (H' = gH), used to conclude H identically zero from a single zero.
Invoked twice in Section 3 to force h ≡ 0.
standard math Riemannian connection formulas in Appendix A, Lemma A.1 and Lemma A.2, for Laplacian, divergence, and convective terms in spherical coordinates.
Alternative derivation of ODE system (6)-(10); the formulas are standard differential geometry but are reproduced in the appendix.
domain assumption The self-similar axisymmetric ansatz u = f(φ)/ρ eρ + g(φ)/ρ eφ + h(φ)/ρ eθ, B = B(φ)/ρ eθ, p = P(φ)/ρ^2.
Reduces the PDE system to ODEs; follows from (-1)-homogeneity and axisymmetry, but it is a structural premise of the whole proof.
domain assumption Regularity of B(φ) and P(φ) on the closed interval [0,π] (or [0,π/2]), including vanishing boundary values at the axis and endpoint continuity used in the maximal-interval contradiction.
Needed for B(0)=B(π)=0, H integral identity, and P(α)=0; not explicitly discussed as a regularity assumption.
domain assumption The external classification Theorem 5.1 of Kang, Miura and Tsai [KMT18] for axisymmetric self-similar Navier-Stokes solutions in the half-space.
Used at the end of Theorem 1.2 to conclude u=0 after reducing the MHD problem to the Navier-Stokes problem; the paper does not reprove it.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Axisymmetric self-similar solutions to the MHD equations without magnetic diffusion." pith.science (2026). https://pith.science/paper/PIRDG2OB

@misc{pith2026250620131,
  author       = {Pith},
  title        = {Pith review of: Axisymmetric self-similar solutions to the MHD equations without magnetic diffusion},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PIRDG2OB}},
  note         = {Machine review of arXiv:2506.20131}
}

abstract

We study the axisymmetric self-similar solutions $(\mathbf{u},\mathbf{B})$ to the stationary MHD equations without magnetic diffusion, where $\mathbf{B}$ has only the swirl component. Our first result states that in $\mathbb{R}^3\setminus\{0\}$, $\mathbf{u}$ is a Landau solution and $\mathbf{B}=0$. Our second result proves the triviality of axisymmetric self-similar solutions in the half-space $\mathbb{R}^3_+$ with the no-slip boundary condition or the Navier slip boundary condition.

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

11 extracted references · 10 canonical work pages

[1]

G. K. Batchelor. An introduction to fluid dynamics . Cambridge Mathematical Library. Cambridge University Press, Cambridge, paperback edition, 1999

work page 1999
[2]

Smooth or singular solutions to the N avier- S tokes system? J

Marco Cannone and Grzegorz Karch. Smooth or singular solutions to the N avier- S tokes system? J. Differential Equations , 197(2):247--274, 2004

work page 2004
[3]

Hamilton's R icci flow , volume 77 of Graduate Studies in Mathematics

Bennett Chow, Peng Lu, and Lei Ni. Hamilton's R icci flow , volume 77 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI; Science Press Beijing, New York, 2006

work page 2006
[4]

Riemannian geometry

Manfredo Perdig\ ao do Carmo. Riemannian geometry . Mathematics: Theory & Applications. Birkh\"auser Boston, Inc., Boston, MA, portuguese edition, 1992

1992
[5]

Green tensor of the S tokes system and asymptotics of stationary N avier- S tokes flows in the half space

Kyungkeun Kang, Hideyuki Miura, and Tai-Peng Tsai. Green tensor of the S tokes system and asymptotics of stationary N avier- S tokes flows in the half space. Adv. Math. , 323:326--366, 2018

work page 2018
[6]

L. Landau. A new exact solution of N avier- S tokes equations. C. R. (Doklady) Acad. Sci. URSS (N.S.) , 43:286--288, 1944

work page 1944
[7]

H. B. Squire. The round laminar jet. Quart. J. Mech. Appl. Math. , 4:321--329, 1951

work page 1951
[8]

Lectures on N avier- S tokes equations , volume 192 of Graduate Studies in Mathematics

Tai-Peng Tsai. Lectures on N avier- S tokes equations , volume 192 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI, 2018

work page 2018

Show all 11 references

[9]

One-point singular solutions to the N avier- S tokes equations

Gang Tian and Zhouping Xin. One-point singular solutions to the N avier- S tokes equations. Topol. Methods Nonlinear Anal. , 11(1):135--145, 1998

1998
[10]

S ver\'ak

V. S ver\'ak. On L andau's solutions of the N avier- S tokes equations. volume 179, pages 208--228. 2011. Problems in mathematical analysis. No. 61

2011
[11]

Point singularities of solutions to the stationary incompressible MHD equations

Shaoheng Zhang, Kui Wang, and Yun Wang. Point singularities of solutions to the stationary incompressible MHD equations. arXiv:2504.09607 , 2025

2025 arXiv

Pith tools

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

[1] [1]

G. K. Batchelor. An introduction to fluid dynamics . Cambridge Mathematical Library. Cambridge University Press, Cambridge, paperback edition, 1999

work page 1999

[2] [2]

Smooth or singular solutions to the N avier- S tokes system? J

Marco Cannone and Grzegorz Karch. Smooth or singular solutions to the N avier- S tokes system? J. Differential Equations , 197(2):247--274, 2004

work page 2004

[3] [3]

Hamilton's R icci flow , volume 77 of Graduate Studies in Mathematics

Bennett Chow, Peng Lu, and Lei Ni. Hamilton's R icci flow , volume 77 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI; Science Press Beijing, New York, 2006

work page 2006

[4] [4]

Riemannian geometry

Manfredo Perdig\ ao do Carmo. Riemannian geometry . Mathematics: Theory & Applications. Birkh\"auser Boston, Inc., Boston, MA, portuguese edition, 1992

1992

[5] [5]

Green tensor of the S tokes system and asymptotics of stationary N avier- S tokes flows in the half space

Kyungkeun Kang, Hideyuki Miura, and Tai-Peng Tsai. Green tensor of the S tokes system and asymptotics of stationary N avier- S tokes flows in the half space. Adv. Math. , 323:326--366, 2018

work page 2018

[6] [6]

L. Landau. A new exact solution of N avier- S tokes equations. C. R. (Doklady) Acad. Sci. URSS (N.S.) , 43:286--288, 1944

work page 1944

[7] [7]

H. B. Squire. The round laminar jet. Quart. J. Mech. Appl. Math. , 4:321--329, 1951

work page 1951

[8] [8]

Lectures on N avier- S tokes equations , volume 192 of Graduate Studies in Mathematics

Tai-Peng Tsai. Lectures on N avier- S tokes equations , volume 192 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI, 2018

work page 2018