{"id":"43529387-2437-44f4-84ef-360f5d246319","arxiv_id":"2608.05994","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors construct steady MHD disk flows whose inviscid interior is a rigidly rotating core, with the two rotation rates fixed by an explicit MHD-Wood law.","lead":"This paper proves a rigorous MHD counterpart of the Prandtl-Batchelor selection principle for steady flows in a disk, together with flux expulsion. A smart generalist should read it because it couples two classical fluid-dynamical mechanisms that were previously studied with the velocity prescribed, and it produces explicit formulas for the selected interior flow and magnetic field.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the non-Alfvénic restriction is an explicit scope condition, not a flaw in the stated theorem.","rationale":"The reader's weakest assumption identifies the non-Alfvénic condition, but this is a theorem hypothesis, not a hidden premise. The paper explicitly states the degeneration and uses it consistently: the same factor a² − b² appears as the determinant in the leading Prandtl coercivity, as the solvability factor in the Euler matching, and as the stability coefficient in the energy norm. I checked the algebraic structure of the Prandtl linearization, the derivation of the MHD–Wood identities, the Euler recursion, and the closing of the energy estimates. The powers of ε close correctly: the residual is ε^10, the linear estimate contributes ε^{-1}, and the nonlinear contraction gives error size ε^9. The refined estimate's delicate term involving Πθ is handled by substituting the magnetic error equation, and the residual bound (2.53) includes the needed control of ∂_θ R_g. I therefore recommend no change to the ACCEPT verdict. The only reservation is that the Alfvénic regime |α| = |β| remains untreated, but that is an explicitly stated limitation of the theorem rather than a flaw in the argument.","tokens_in":73993,"tokens_out":35923,"duration_ms":332283,"concrete_test":"As a verification step, independently re-derive the refined estimate in Lemma 3.9, in particular the bound for I1,8, where Πθ is eliminated through the second magnetic-potential error equation. Confirm that the resulting F2,θ term is controlled by the residual bound (2.53) with no loss of a power of ε, since this is the step that prevents derivative loss and is the most intricate point in the stability estimate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I reviewed the proof with attention to the places where the central claim could fail internally. The non-Alfvénic condition |α| ≠ |β| is explicitly assumed in Theorem 1.1 and is exactly the degeneracy of the three linearized mechanisms: det M_{α,β} = 1 − β²/α² in the Prandtl coercivity, the factor a² − b² in the Euler matching, and the coefficient a − b²/a in the stability norm. The proof therefore does not silently exclude a regime it purports to cover. Within the stated assumptions, the matched asymptotic construction closes: the MHD–Wood identities are derived from the von Mises system without circular solvability; the Prandtl fixed point is coercive for |α| ≠ |β|; the Euler recursion uses only ∆(rv_e^(k)) = 0, which follows from a² − b² ≠ 0; and the nonlinear error estimate closes with the three-scale norm E + P + L, with residual powers ε^10 mapping to ε^9 via the linear estimate. I did not find an unjustified identity or a missing boundary condition. The main limitation is the unaddressed Alfvénic regime, which is a scope restriction rather than an internal inconsistency.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the simultaneous vanishing-viscosity and vanishing-resistivity limit of steady incompressible MHD flows in the unit disk. For boundary velocity a small non-axisymmetric perturbation of rigid rotation with mean angular speed α, and tangential magnetic trace with mean β, the authors construct exact steady solutions of system (1.3) under the assumptions α≠0 and |α|≠|β|. The construction uses a matched asymptotic expansion consisting of a periodic MHD–Prandtl boundary layer, global Euler–MHD core profiles through order ε^11, residual estimates of order ε^10, a three-scale linear stability estimate, and a contraction argument for the nonlinear error equations. The paper proves convergence on compact interior subdisks to the Couette-type ideal core u_e=ar e_θ, B_e=br e_θ, where b=β and a²=α²+(η²/2π)∫₀^{2π}(f²−κ^{-1}hat_f²)dθ. An appendix presents a separate conditional rigidity argument for single-eddy ideal limits.","tokens_in":74170,"tokens_out":21756,"duration_ms":197116,"significance":"If correct, this is the first fully coupled steady MHD realization of Prandtl–Batchelor selection together with selective flux expulsion in a closed single-eddy geometry. The Lorentz force is retained at leading order, so the result goes beyond kinematic flux-expulsion models, and it identifies the boundary-circulation mode that survives expulsion. The proof has several checkable and genuinely structural components: the MHD–Wood identities are derived from the von Mises formulation without circular solvability; the non-Alfvénic coercivity condition det M_{α,β}=1−β²/α² is explicit and is exactly the degeneracy of the linearized mechanisms; and the stability norm uses the coupled quantity G whose L-scale is two powers of ε stronger than the naive diffusive control. The exclusion of the Alfvénic regime |α|=|β| is an explicitly stated hypothesis rather than a hidden failure, and it is consistently the source of the a²−b² factors in the Prandtl coercivity, the Euler matching, and the final stability estimate.","major_comments":[],"minor_comments":[{"comment":"The sign convention in B^ε=∇⊥Π^ε and the derivation of u^ε·∇Π^ε−ε²∆Π^ε=C₀ should be accompanied by an explicit verification of the vector identity, since the same identity is later used in Appendix B in the form B·∇Φ=ε²(j−2β).","section":"§2.1, Eq. (2.1)"},{"comment":"The notation w_{≠,ψ} and g_{≠,ψ} is used in the weighted estimates without definition; the authors should define w_{≠}=w−w₀ and w_{≠,ψ}=∂_ψ w_{≠}.","section":"§2.2.2, Step 5"},{"comment":"The cancellation b∫r(∆Π)_θ Φ = −b∫r∆Π Φ_θ is used silently; citing periodicity and the self-adjointness of the polar Laplacian would improve readability.","section":"§3.2.2, proof of Lemma 3.5"},{"comment":"The notation for the magnetic perturbation is inconsistently rendered: the abstract displays “fhat²” while the displayed formula uses \\(\\hat f^2\\). These should be unified.","section":"Abstract and after Eq. (1.2)"},{"comment":"In the Fourier representation of v_e^{(1)}, the zero mode is asserted but not shown to vanish; this follows from the zero mean of v_p^{(1)}(·,0) and should be stated explicitly.","section":"§2.2.3, Proposition 2.4"},{"comment":"The argument that regular levels are connected would benefit from a sentence explaining why a nested pair of level curves would force an additional critical point in the annulus between them.","section":"Appendix B, Lemma B.1"}],"recommendation":"accept","confidential_remarks":"The non-Alfvénic restriction is the main scope condition of the paper. I do not regard it as a flaw of the stated theorem, but a brief discussion of the degenerate Alfvénic case as an open problem would help readers. The conditional rigidity appendix is not used in the main proof and is clearly labeled as a separate result; this is acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is the bottom line: this is a genuinely new rigorous result, and the proof looks credible within its stated assumptions. The paper constructs steady MHD solutions in a disk whose interior limits are a rigid rotation and a uniform magnetic rotation, with the velocity rotation selected by a new coupled MHD–Wood law and the magnetic rotation fixed by the boundary circulation. That is the first fully coupled steady realization of Prandtl–Batchelor selection plus flux expulsion, and it is not the passive-field or small-field regime that most prior work stays in. The magnetic field back-reacts at leading order; the paper handles that rather than assuming it away.\n\nTechnically, the strongest part is the MHD–Prandtl layer treated as a genuinely two-field problem, with nonlocal normal components and a magnetic potential satisfying only a Neumann condition. The new quantity G = aΠ_θ - bΦ_θ - ... is a smart way to extract the strong scale L that replaces the missing magnetic boundary control; that is the hardest part of the analysis, and it is exactly where the non-Alfvénic condition enters. I also looked for circularity: Lemma 2.1 derives the Wood identities by exact integration, and the fixed point in Proposition 2.2 is constructed before the identities identify a and b. The circularity concern does not land.\n\nSoft spots, in proportion. The condition |α| ≠ |β| is real and load-bearing in three places (Prandtl coercivity, Euler matching, stability coefficient). The authors state it openly and explain it as the Alfvénic degeneracy; it is a scope restriction, not an inconsistency. But if you want the Alfvénic regime, you get nothing from this paper. Also, the proof is long and dense; I did not machine-check the energy estimates, and a referee should expect to grind through the residual and stability estimates. The boundary condition on B is a prescribed tangential trace (Neumann for the potential), justified via exterior matching; that is fine but narrower than a complete insulating-wall model.\n\nWho this is for: people working in mathematical MHD boundary layers, Prandtl–Batchelor theory, and vanishing-viscosity limits. It deserves a serious referee. I would send it out.","headline":"A genuinely new rigorous construction of coupled Prandtl–Batchelor selection and flux expulsion in steady MHD; the non-Alfvénic restriction is explicit and the proof holds up on reading.","tokens_in":74746,"tokens_out":2361,"would_cite":true,"duration_ms":24591,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76W05","76D10","35Q35","35B25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper constructs steady MHD disk flows whose interior limits are rigid rotations with constant vorticity and current, selected by a new MHD–Wood law that couples the two classical mechanisms.","keywords":["Prandtl–Batchelor principle","flux expulsion","MHD boundary layer","Wood formula","vanishing viscosity and resistivity","steady MHD in a disk","non-Alfvénic condition","magnetohydrodynamics"],"falsifier":"Solve the steady MHD system (1.3) computationally in the disk at small $\\varepsilon$ and $\\eta$ with smooth zero-mean $f$ and $\\hat f$, and compare the measured interior vorticity and current density with the constants $2a$ and $2b$ given by (1.2); a systematic mismatch beyond $O(\\varepsilon)$ would refute the selection law. Repeating the same computation at $\\alpha = \\beta$ would test whether the excluded Alfvénic regime still admits a rigid rotating core or genuinely degenerates.","tokens_in":2261,"feed_emoji":"🧲","tokens_out":8981,"duration_ms":131360,"temperature":0.7,"pith_summary":"This paper tries to prove that a magnetized eddy in a disk has a single, jointly selected inviscid core: as viscosity and resistivity vanish simultaneously, every small nonaxisymmetric perturbation of a rotating boundary pulls the interior into rigid rotation about the center, with both vorticity and out-of-plane current constant. The two constants are not free data. A new MHD–Wood law fixes the magnetic constant as the imposed mean boundary circulation and fixes the velocity constant by a quadratic balance of the kinetic and magnetic boundary oscillations. If true, this is the first steady, fully coupled instance in which Prandtl–Batchelor vorticity selection and magnetic flux expulsion are resolved together, with the field allowed to be as strong as the flow. It also explains what survives expulsion: the circulation-carrying uniform magnetic rotation survives, while oscillatory magnetic modes are confined to a thin layer, so expulsion is selective rather than absolute.","feed_headline":"One two-term law sets both cores in a magnetized disk eddy","feed_subtitle":"Boundary data fix the magnetic rotation and balance the velocity rotation; the rest of the field is expelled to a thin layer.","key_machinery":"The load-bearing object is a periodic steady MHD–Prandtl boundary layer solved in von Mises coordinates, with the square velocity $u_p^2$ and the tangential magnetic field $g_p$ as unknowns. Around the constant state $(\\alpha^2, \\beta)$ the principal tangential combinations are $\\partial_\\theta w - 2\\beta\\, \\partial_\\theta g$ and $\\partial_\\theta g - \\frac{\\beta}{2\\alpha^2}\\,\\partial_\\theta w$, whose coefficient matrix has determinant $1 - \\beta^2/\\alpha^2$; invertibility of this matrix is exactly the non-Alfvénic condition, and the same factor $a^2 - b^2$ appears later in Euler matching and in the stability energy. Two averaged identities across the layer, $\\frac{d^2}{d\\psi^2}\\int g\\, d\\theta = 0$ and $\\frac{d^2}{d\\psi^2}\\int (\\kappa u^2 - g^2)\\,d\\theta = 0$, produce the two scalars of the MHD–Wood law. The stability argument is carried by the coupled quantity $G = a\\Pi_\\theta - b\\Phi_\\theta$ plus zero-mode correction terms, which gives strong magnetic control (an $\\varepsilon^{-2}$ weighted norm) that the Neumann magnetic boundary condition does not supply; three energies $E, P, L$ with different $\\varepsilon$ scales close the estimate through the coercive factor $|a^2 - b^2|/|a|$.","core_discovery":"The central claim, stated as Theorem 1.1, is that for $\\alpha \\neq 0$, $|\\alpha| \\neq |\\beta|$ and sufficiently small $\\varepsilon, \\eta$, the steady MHD system with boundary traces $u\\cdot t = \\alpha+\\eta f(\\theta)$ and $B\\cdot t = \\beta+\\eta \\hat f(\\theta)$ has exact solutions that, on every compact interior subdisk, converge to $u_e = ar\\, e_\\theta$ and $B_e = br\\, e_\\theta$ with $b = \\beta$ and $a^2 = \\alpha^2 + \\frac{\\eta^2}{2\\pi}\\int_0^{2\\pi}\\!(f^2 - \\kappa^{-1}\\hat f^2)\\,d\\theta$, $a$ having the sign of $\\alpha$. Equivalently, the limiting vorticity is $2a$ and the limiting current density is $2b$. The velocity core is the output of a coupled kinetic–magnetic boundary-layer balance; the magnetic core is the output of the global circulation constraint, and the constant $2\\beta\\varepsilon^2$ in the magnetic-flux equation independently confirms $b = \\beta$. The construction is genuinely two-field: the Lorentz force enters the leading-order momentum balance, so the magnetic field actively changes the flow and need not be weak.","pith_inferences":["The formula suggests a practical test that the paper does not run: at moderate $\\varepsilon$ and $\\eta$ with measured boundary traces $f, \\hat f$, compare the interior vorticity and current against (1.2); agreement would confirm the selection law, disagreement would localize where the asymptotics break.","Because the MHD–Wood law is derived from two layer-averaged balances rather than from a specific trace shape, the same two-constant selection should extend to other simply connected eddies and to annuli, where a second boundary would supply a second magnetic datum; this is an inference, not a result of the paper.","The conditional rigidity argument of Appendix B suggests that, under a single-eddy and non-Alfvénic hypothesis, any locally $C^1$-convergent ideal limit must be the same rigid core without needing the explicit expansion; if combined with a uniqueness statement for the boundary layer, this could yield a selection theorem that does not construct the layer at all.","The Alfvénic endpoint $|\\alpha|=|\\beta|$, excluded here, may be a genuine bifurcation point rather than a removable technical assumption: at $a = \\pm b$ one Elsässer component vanishes, and the two-field interaction degenerates, so the selected-core structure may change or fail there; this is left open by the paper."],"forward_implications":["When the imposed magnetic circulation vanishes ($\\beta = 0$), every compact interior subdisk has magnetic field tending to zero: complete flux expulsion.","When $\\beta \\neq 0$, the uniform magnetic rotation $\\beta r e_\\theta$ survives in the core while all nonaxisymmetric magnetic modes are confined to an $O(\\varepsilon)$ boundary layer: selective, not absolute, expulsion.","The selected velocity rotation $a^2$ is raised by the kinetic boundary oscillation $\\int f^2$ and lowered by the magnetic boundary oscillation $\\kappa^{-1}\\int \\hat f^2$, so boundary magnetic wiggles act on the core through the diffusivity ratio $\\kappa$.","The limiting vorticity and current density are constants $2a$ and $2b$, with uniform convergence on compact subdisks and an $O(\\varepsilon)$ error in the full disk including the boundary layer.","The construction requires $|\\alpha| \\neq |\\beta|$; at the Alfvénic balance the coercivity factor $a^2 - b^2$ vanishes in the Prandtl estimate, the Euler matching, and the final stability energy, so the theorem leaves the Alfvénic regime open."],"supporting_citations":[{"why":"Provides the rigorous hydrodynamic Prandtl–Batchelor disk construction whose scalar Wood formula the present MHD–Wood law generalizes, and contributes the von Mises fixed-point strategy.","marker":"[11]"},{"why":"Gives the classical Wood formula for the selected rotation in hydrodynamic closed-streamline flow, the limiting case that (1.2) extends to two fields.","marker":"[55]"},{"why":"States the Prandtl–Batchelor principle that weak viscosity homogenizes vorticity in closed eddies, the velocity-side mechanism that the paper couples to the magnetic field.","marker":"[3]"},{"why":"Supplies the classical kinematic prediction of magnetic flux expulsion by an eddy that the zero-circulation case of the theorem realizes as a rigorous steady coupled limit.","marker":"[51]"},{"why":"Sets out the circular-eddy and shear configurations for flux-expulsion time scales whose uniform-field setup corresponds to the zero-circulation boundary data.","marker":"[36]"},{"why":"Provides the radial symmetry theorem used in the conditional rigidity argument of Appendix B to derive constant vorticity and current from a semilinear elliptic equation.","marker":"[14]"},{"why":"Introduces Elsässer variables, in which the degeneracy $a = \\pm b$ corresponds to one vanishing characteristic field and motivates the non-Alfvénic hypothesis.","marker":"[8]"},{"why":"Gives the classical flux-expulsion picture in rotating cylinders and spheres with a quasi-uniform field, another zero-circulation counterpart of the selective result.","marker":"[38]"}],"fun_headline_variants":["Magnetized disk eddy: single law nails both cores","Active magnetic field selects twin cores in disk MHD","Flux expulsion meets Batchelor: one law for two cores","One balance sets velocity and magnetic cores in eddy","Disk MHD: coupled law fixes vorticity and current"],"cache_read_input_tokens":76928,"weakest_assumption_plain":"Everything depends on the boundary means being unequal: if $|\\alpha| = |\\beta|$, the key coercivity factor $a^2 - b^2$ vanishes and the construction stops, so the theorem says nothing about the Alfvénic case.","fun_headline_variants_meta":{"raw":{"variants":["Magnetized disk eddy: single law nails both cores","Active magnetic field selects twin cores in disk MHD","Flux expulsion meets Batchelor: one law for two cores","One balance sets velocity and magnetic cores in eddy","Disk MHD: coupled law fixes vorticity and current"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000263,"raw_usage":{"total_tokens":1678,"prompt_tokens":1101,"completion_tokens":577,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":717,"completion_tokens_details":{"reasoning_tokens":495}},"tokens_in":717,"tokens_out":577,"duration_ms":6625,"temperature":1.0,"reasoning_tokens":495,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T19:36:51.571275+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the steady MHD system (1.3) computationally in the disk at small $\\varepsilon$ and $\\eta$ with smooth zero-mean $f$ and $\\hat f$, and compare the measured interior vorticity and current density with the constants $2a$ and $2b$ given by (1.2); a systematic mismatch beyond $O(\\varepsilon)$ would refute the selection law. Repeating the same computation at $\\alpha = \\beta$ would test whether the excluded Alfvénic regime still admits a rigid rotating core or genuinely degenerates.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the rigorous hydrodynamic Prandtl–Batchelor disk construction whose scalar Wood formula the present MHD–Wood law generalizes, and contributes the von Mises fixed-point strategy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the classical Wood formula for the selected rotation in hydrodynamic closed-streamline flow, the limiting case that (1.2) extends to two fields."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"States the Prandtl–Batchelor principle that weak viscosity homogenizes vorticity in closed eddies, the velocity-side mechanism that the paper couples to the magnetic field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical kinematic prediction of magnetic flux expulsion by an eddy that the zero-circulation case of the theorem realizes as a rigorous steady coupled limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets out the circular-eddy and shear configurations for flux-expulsion time scales whose uniform-field setup corresponds to the zero-circulation boundary data."},{"cited_title":"Gidas, W","cited_arxiv_id":null,"evidence_quote":"Provides the radial symmetry theorem used in the conditional rigidity argument of Appendix B to derive constant vorticity and current from a semilinear elliptic equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces Elsässer variables, in which the degeneracy $a = \\pm b$ corresponds to one vanishing characteristic field and motivates the non-Alfvénic hypothesis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the classical flux-expulsion picture in rotating cylinders and spheres with a quasi-uniform field, another zero-circulation counterpart of the selective result."}],"review_version":1}