Stability of vertically charged steady magnetic field in 3D incompressible magneto-micropolar fluids without magnetic and angular viscosity in a strip domain
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-08 01:20 UTCglm-5.2pith:2JRGAP32record.jsonopen to challenge →
The pith
Vertically charged magnetic field stabilizes 3D micropolar fluid
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors prove that for the 3D incompressible magneto-micropolar system with zero magnetic and angular viscosity in a strip domain, a vertically oriented constant magnetic field is a stable equilibrium. Specifically, if the initial perturbation is sufficiently small in a high-order Sobolev norm (with the regularity index N ≥ 4), then there exists a unique global classical solution whose energy remains bounded by the initial data for all time, and the perturbation decays to zero at a rate of (1+t)^{4−2N}. This is the first global well-posedness and stability result for this system under these physical conditions.
What carries the argument
The Lagrangian coordinate reformulation reduces the magnetic field to b = e₃ + ∂₃η, which produces the wave-structure term ∂²₃η in the momentum equation. The Helmholtz projection P eliminates pressure from the velocity equation, yielding a projected coupled system for Q = P(∂²₃η) and V = P(∇×w). A novel anisotropic estimate (Lemma 3.2) controls the horizontal gradient of the potential ϕ in the Helmholtz decomposition of ∇×w by its vertical derivative and a small multiple of the full curl, enabling recovery of w from its projection. Elliptic estimates for u are derived through the vorticity equation of w rather than the Stokes system, breaking a circular dependency. The two-tier energy method
Load-bearing premise
The background magnetic field must be exactly perpendicular to the horizontal boundary (b₀ = e₃). This specific orientation is what allows the reduction b = e₃ + ∂₃η and generates the ∂²₃η wave structure in the momentum equation. The authors note results may hold for non-perpendicular fields but the proof does not extend to that case. Additionally, the entire argument requires an a priori smallness assumption on the solution norm to control the Lagrangian coordinate perturb
What would settle it
If one could exhibit initial data arbitrarily close to the vertical steady field for which the solution develops a singularity in finite time, the theorem would fail. More concretely, if the derivative-loss mechanism from the antisotropic curl coupling between u and w cannot be compensated by the magnetic stabilizing term ∂²₃η at the energy level the authors claim, the energy closure would break down.
read the original abstract
This paper intends to understand the regularity and stability problem on the 3D incompressible magneto-micropolar equations with zero magnetic and angular viscosities in a strip domain. The magneto-micropolar system models the electrically conducting micropolar fluid in the presence of a magnetic field. The lack of magnetic diffusion and angular dissipation makes it impossible to prove even small data global well-posedness result, let alone general large data global regularity. This paper presents a steady-state setup around which any perturbations can be shown to be globally regular and stable. More precisely, any small perturbation near a steady magnetic field perpendicular to the horizontal boundary leads to a unique global classical solution. In addition, the solution is shown to converge to the steady state at an almost exponential rate as time goes to infinity. These appear to be the very first rigorous global results on the magneto-micropolar equations concerned here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies the global regularity and stability of the 3D incompressible magneto-micropolar equations with zero magnetic and angular viscosity in a strip domain. The authors consider perturbations near a steady magnetic field perpendicular to the horizontal boundary. By reformulating the system in Lagrangian coordinates, the background field b_0 = e_3 allows the reduction b = e_3 + ∂_3η, which introduces a wave-like structure ∂²₃η in the momentum equation. The proof employs the two-tier energy method of Guo-Tice, combined with a novel use of the Helmholtz projection to handle the coupling between the fluid velocity u and the angular velocity w. The main result (Theorem 2.1) establishes that sufficiently small perturbations yield a unique global classical solution with an almost exponential decay rate (1+t)^{4-2N} for N ≥ 4.
Significance. This appears to be the first rigorous global well-posedness and stability result for the 3D incompressible magneto-micropolar system without magnetic and angular viscosity. The key technical innovation is the use of the Helmholtz projection to decompose ∇×w and recover estimates on w via the anisotropic estimates of Lemma 3.2, which is necessitated by the incompressible pressure structure differing from the compressible case in [11]. The result provides a falsifiable prediction: stability holds specifically when the background field is perpendicular to the boundary. The energy scheme is extensive, spanning non-spatial, horizontal, and vertical derivative estimates, synthesized through an inductive argument in Section 6.
major comments (3)
- §6, Step 7 (transition from (6.45)–(6.46) to (6.47)): The inductive scheme in Steps 5–7 proceeds downward in k from k=2n−1 to k=0, absorbing lower-order W terms at each step. At the final step, (6.46) takes k=0 in (6.29) and absorbs ∥∂_t^l W∥²_{0,2n−1} from the right side of (6.45). However, (6.29) at k=0 contains the term σ∑_l ∥∂_t^l W∥²_{0,2n−1} on its right side (inherited from (6.23)–(6.24)), and the text states 'provided σ is suitably small' without explicitly verifying that a single σ > 0 works simultaneously across all 2n−1 inductive steps. Since the anisotropic estimate (3.11) produces a σ that depends on k and λ, one needs to confirm that the σ selected at each inductive step (k=2n−1,...,1) is compatible with the σ required for the final absorption at k=0. Please add a remark or a brief argument clarifying that a uniform σ > 0 suffices across all steps, or adjust the induction.
- §6, (6.29): The right-hand side of (6.29) contains the term ∑_l ∥∂_t^l w∥²_{0,2n}, which is of the same or higher order than the dissipation terms on the left. This term is carried through the inductive chain (6.38)–(6.45) and ultimately absorbed in (6.47) via reference to (5.6) and (5.25). The absorption of a 0,2n-norm quantity by ¯D*_n and ¯D♯_n should be made more explicit, as it is load-bearing for the closure of the dissipation estimate (6.47) that feeds into Lemma 7.2.
- §1.2, Remark 2.2: The reduction b = e_3 + ∂_3η and the resulting wave structure ∂²₃η in (1.7b) critically depend on b_0 = e_3 being perpendicular to the boundary. The remark that results 'may still hold true' for non-perpendicular fields is speculative. While this does not affect the correctness of the current result, the authors should clarify whether the perpendicularity is used only for the reduction (1.6)–(1.7b) or also more structurally in the anisotropic estimates of Lemma 3.2, which rely on the strip geometry.
minor comments (5)
- §2, (2.5): The definition of G_{2N}(t) involves the parameter θ > 0, but the constraint on θ (namely N−2 ≥ 1+θ) is only stated later in the proof of Lemma 7.4. It would help the reader to note this constraint when θ is first introduced.
- §3.2, Lemma 3.2: The proof uses frequency localization into |ξ| ≤ 1 and |ξ| ≫ 1 regimes, but the intermediate regime 1 < |ξ| < K̄ is not explicitly discussed. A brief remark that the low-frequency estimate (3.13) extends to any finite K would improve clarity.
- §5.1, Lemma 5.1: The proof is omitted with reference to 'similar arguments as in the proof of Lemma 4.2.' Given that this lemma is used at multiple critical points (n=2N and n=N+2), a slightly more detailed sketch, particularly for the highest-order terms in (5.3), would be beneficial.
- Typographical: In (6.3a), the equation reads '(µ+ζ)∂_t Q + Q + ζV = H'. The coefficient (µ+ζ) on ∂_t Q appears unusual given the structure of (6.3b); please verify this is correct and not a typo for the time derivative coefficient.
- §7, Lemma 7.1, (7.1): The term (E_n)² appears on the right side. The text explains this arises from X_n ≲ (E_n)² via (5.2), but the mechanism by which a quadratic term in the energy arises from the perturbation estimates could be stated more explicitly.
Simulated Author's Rebuttal
We thank the referee for a careful reading and for identifying three substantive points, all of which concern the closure of the energy scheme in Section 6. We address each point below. In brief: (1) we agree that the uniformity of σ across the inductive steps should be made explicit and will add a remark; (2) we agree that the absorption of the ∑_l ∥∂_t^l w∥²_{0,2n} term in (6.47) deserves a more transparent justification and will expand the text accordingly; (3) we will clarify the role of perpendicularity. No standing objections remain.
read point-by-point responses
-
Referee: §6, Step 7 (transition from (6.45)–(6.46) to (6.47)): The inductive scheme in Steps 5–7 proceeds downward in k from k=2n−1 to k=0, absorbing lower-order W terms at each step. At the final step, (6.46) takes k=0 in (6.29) and absorbs ∥∂_t^l W∥²_{0,2n−1} from the right side of (6.45). However, (6.29) at k=0 contains the term σ∑_l ∥∂_t^l W∥²_{0,2n−1} on its right side (inherited from (6.23)–(6.24)), and the text states 'provided σ is suitably small' without explicitly verifying that a single σ > 0 works simultaneously across all 2n−1 inductive steps. Since the anisotropic estimate (3.11) produces a σ that depends on k and λ, one needs to confirm that the σ selected at each inductive step (k=2n−1,...,1) is compatible with the σ required for the final absorption at k=0. Please add a remark or a brief argument clarifying that a uniform σ > 0 suffices across all steps, or adjust the induction.
Authors: We agree that this point deserves explicit clarification. The key observation is that the σ appearing in (6.23)–(6.24) and hence in (6.29) is the same parameter σ > 0 from Lemma 3.2, which is arbitrary (the lemma holds 'for any σ > 0'). At each inductive step k = 2n−1, 2n−2, ..., 1, the term σ∑_l ∥∂_t^l W∥²_{0,2n−1} on the right side of (6.29) is absorbed by the corresponding term ∑_l ∥∂_t^l W∥²_{k,2n−k−1} on the left side of the estimate at the next lower k, since ∥·∥_{k,2n−k−1} ≥ ∥·∥_{0,2n−1} for k ≥ 1. The crucial point is that the same σ > 0 appears at every step (it is not re-selected), and the absorption at each intermediate step is exact (no smallness of σ is needed—only the final absorption at k = 0 requires σ to be small relative to the dissipation on the left of (6.46)). More precisely: at each intermediate step, the term σ∑_l ∥∂_t^l W∥²_{0,2n−1} is bounded by σ∑_l ∥∂_t^l W∥²_{k−1,2n−k}, which is part of the dissipation on the left side of the estimate at level k−1, and is absorbed there regardless of the size of σ (it is a subset of the left-side dissipation). Only at the final step (k = 0, equation (6.46)) do we need σ sufficiently small so that σ∑_l ∥∂_t^l W∥²_{0,2n−1} can be absorbed by the full dissipation ∑_l ∥∂_t^l W∥²_{0,2n−1} on the left of (6.46). Since the lemma provides the estimate for any σ > 0, we simply choose σ small enough at the outset. We will add a remark to this effect in the revision. revision: yes
-
Referee: §6, (6.29): The right-hand side of (6.29) contains the term ∑_l ∥∂_t^l w∥²_{0,2n}, which is of the same or higher order than the dissipation terms on the left. This term is carried through the inductive chain (6.38)–(6.45) and ultimately absorbed in (6.47) via reference to (5.6) and (5.25). The absorption of a 0,2n-norm quantity by ¯D*_n and ¯D♯_n should be made more explicit, as it is load-bearing for the closure of the dissipation estimate (6.47) that feeds into Lemma 7.2.
Authors: The referee is correct that this absorption is load-bearing and should be made more explicit. We trace the argument as follows. The term ∑_l ∥∂_t^l w∥²_{0,2n} on the right side of (6.29) arises from the Hodge estimate (6.28), specifically from the term ∑_l ∥∂_t^l w∥²_{0,2n} in (6.28), which is the lower-order remainder in the div-curl decomposition. By the definition of the anisotropic norm (2.1), ∥∂_t^l w∥²_{0,2n} = ∑_{α₁+α₂≤2n} ∥∂₁^{α₁}∂₂^{α₂} ∂_t^l w∥²₀. For l = 0, this is ∥w∥²_{0,2n}, which is contained in ¯E*_n (see (5.5)) and its dissipation ∥w∥²_{1,2n} is contained in ¯D*_n (see (5.6)). For l ≥ 1, ∥∂_t^l w∥²_{0,2n} is contained in ¯E*_n and ∥∂_t^l w∥²_{1,2n} in ¯D*_n as well (the sum in (5.5)–(5.6) runs over j = 0, ..., n−1 with ∥∂_t^j w∥²_{0,2n−2j}, and for l ≤ n−1 we have 2n−2l ≥ 2n only when l = 0; for l ≥ 1, 2n−2l ≤ 2n−2, so ∥∂_t^l w∥²_{0,2n} is actually controlled by ∥∂_t^l w∥²_{0,2n−2l+2l} which is bounded by interpolation between the terms already in ¯D*_n). In the transition from (6.45) to (6.47), the term ∑_l ∥∂_t^l w∥²_{0,2n} is bounded by ¯D*_n + ¯D♯_n using (5.6) and (5.25), as stated. We will expand the text at (6.47) to spell out this decomposition explicitly, separating the l = 0 and l ≥ 1 cases and indicating which terms in ¯D*_n and ¯D♯_n absorb each piece. revision: yes
-
Referee: §1.2, Remark 2.2: The reduction b = e_3 + ∂_3η and the resulting wave structure ∂²₃η in (1.7b) critically depend on b_0 = e_3 being perpendicular to the boundary. The remark that results 'may still hold true' for non-perpendicular fields is speculative. While this does not affect the correctness of the current result, the authors should clarify whether the perpendicularity is used only for the reduction (1.6)–(1.7b) or also more structurally in the anisotropic estimates of Lemma 3.2, which rely on the strip geometry.
Authors: We thank the referee for this observation. To clarify: the perpendicularity of b_0 = e_3 is used in two distinct places. First, and most critically, it is used for the reduction (1.6)–(1.7b): the identity b·∇_A b = ∂²₃η holds specifically because b = e_3 + ∂₃η, which requires b_0 = e_3. This is what produces the wave-like structure ∂²₃η in the momentum equation and is essential to the entire energy scheme. Second, the anisotropic estimates of Lemma 3.2 rely on the strip geometry Ω = ℝ² × (0,1) and the explicit Fourier representation in the horizontal variables, but they do not themselves depend on the direction of b_0; they depend only on the domain geometry and the boundary conditions. So the perpendicularity enters structurally through the reduction (1.6)–(1.7b), not through Lemma 3.2. For a non-perpendicular constant background field b_0, the reduction b = (I + ∇η)b_0 still holds, but b·∇_A b would no longer simplify to a pure ∂²₃η term, and the wave structure would be lost or altered. Whether a suitable substitute structure exists is genuinely open. We will revise Remark 2.2 to make this distinction clear and to tone down the speculative language, explicitly stating that perpendicularity is used for the reduction (1.6)–(1.7b) and that the anisotropic estimates of Lemma 3.2 depend on the strip geometry but not on the direction of b_0. revision: yes
Circularity Check
No significant circularity found; the derivation is self-contained with standard a priori assumptions and independent external benchmarks.
full rationale
The paper proves global stability of a vertically charged steady magnetic field for 3D incompressible magneto-micropolar equations via a two-tier energy method in Lagrangian coordinates. The derivation chain is self-contained: (1) The a priori assumption G_{2N}(T) ≤ δ (eq. 4.1) is a standard bootstrap assumption for local solutions, not a definition of the output. The closure in Lemma 7.4 shows G_{2N}(t) ≲ E_{2N}(0) + F_{2N}(0), which is strictly smaller than δ when initial data is small — this is the standard bootstrap closure, not circular reasoning. (2) The self-citation to [11] (compressible version by overlapping authors Feng-Hong-Zhu) is explicitly distinguished: the paper states 'the crucial analysis in [11] cannot be accessed for the problem considered in this paper' because the incompressible pressure involves both u and w, unlike the compressible density. The new Helmholtz projection approach (Lemma 3.2) is developed precisely to overcome this, providing genuinely independent content. (3) The citation to [12] (Guo-Tice) provides the two-tier energy framework, which is an externally developed method, not a self-citation. (4) The reduction b = e_3 + ∂_3η (eq. 1.6) follows from det(I_3 + ∇η) = 1, which is derived from the PDE structure (∂_t J = J div_A u = 0), not assumed. (5) The inductive σ-absorption in Section 6 uses 'for any σ > 0' in Lemma 3.2, and the text states 'provided σ is suitably small' at each inductive step — this is a standard small-constant absorption argument where finitely many steps each require σ small, and taking the minimum works. This is a technical completeness concern (correctness risk), not circularity. No step in the derivation chain reduces to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- ε₀ (initial data smallness) =
sufficiently small (not numerically specified)
- δ (a priori bound) =
sufficiently small (not numerically specified)
- θ (decay parameter) =
any θ > 0 (chosen so N-2 ≥ 1+θ)
axioms (4)
- domain assumption Local well-posedness of the reduced system (1.7) in Lagrangian coordinates
- standard math Sobolev embedding H^s(Ω) ↪ L^∞(Ω) for s > 3/2 in the strip domain Ω = R² × (0,1)
- standard math Regularity theory for the Stokes system and elliptic equations on infinite strip domains
- domain assumption The initial magnetic field b₀ = e₃ (perpendicular to boundary)
Reference graph
Works this paper leans on
-
[1]
H. Abels and M. Wiegner, Resolvent estimates for the Stokes operator on an infinite layer, Differential Integral Equations,18(2005): 1081–1110
work page 2005
-
[2]
Ahmadi, Universal stability of thermo-magneto-micropolar fluid motions,Int
G. Ahmadi, Universal stability of thermo-magneto-micropolar fluid motions,Int. J. Engng Sci., 15(1976): 853–859
work page 1976
-
[3]
C. Cao, D. Regmi and J. Wu, The 2D MHD equations with horizontal dissipation and horizontal magnetic diffusion,J. Differential Equations,254(2013): 2661–2681
work page 2013
- [4]
-
[5]
J. Cheng and Y. Liu, Global regularity of the 2D magnetic micropolar fluid flows with mixed partial viscosity,Comput. Math. Appl.,70(2015): 66–72
work page 2015
-
[6]
B. Dong and Z. Zhang, Global regularity of the 2D micropolar fluid flows with zero angular viscosity,J. Differential Equations,249(2010): 200–213
work page 2010
-
[7]
G. Duvaut and J.-L. Lions, In´ equations en thermo´ elasticit´ e et magn´ etohydrodynamique,Arch. Ration. Mech. Anal.,46(1972): 241–279
work page 1972
-
[8]
Erdo˘ gan, Polar effects in the apparent viscosity of a suspension,Rheologica Acta,9(1970): 434–438
M. Erdo˘ gan, Polar effects in the apparent viscosity of a suspension,Rheologica Acta,9(1970): 434–438
work page 1970
-
[9]
Eringen, Theory of micropolar fluids,J
A.C. Eringen, Theory of micropolar fluids,J. Math. Mech.,16(1966): 1–18
work page 1966
-
[10]
J. Fan, S. Jiang and G. Nakamura, Stability of weak solutions to equations of magnetohydrody- namics with Lebesgue initial data,J. Differential Equations,251(2011): 2025–2036
work page 2011
-
[11]
Z. Feng, G. Hong and C. Zhu, Global classical solutions for 3d compressible magneto-micropolar fluids without resistivity and spin viscosity in a strip domain,Sci. China Math.,67(2024): 2485–2514
work page 2024
- [12]
- [13]
-
[14]
G. Hong, X. Hou, H. Peng and C. Zhu, Global existence for a class of large solutions to three- dimensional compressible magnetohydrodynamic equations with vacuum,SIAM J. Math. Anal., 49(2017): 2409–2441
work page 2017
-
[15]
X. Hu, Global existence for two dimensional compressible magnetohydrodynamic flows with zero magnetic diffusivity, preprint,arXiv:1405.0274, (2014)
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [16]
- [17]
-
[18]
Ladyzhenskaya,The Mathematical Theory of Viscous Incompressible Flow, Martino Fine Books, 2014
O.A. Ladyzhenskaya,The Mathematical Theory of Viscous Incompressible Flow, Martino Fine Books, 2014
work page 2014
-
[19]
H. Li, X. Xu and J. Zhang, Global classical solutions to 3D compressible magnetohydrodynamic equations with large oscillations and vacuum,SIAM J. Math. Anal.,45(2013): 1356–1387
work page 2013
-
[20]
F. Lin, L. Xu and P. Zhang, Global small solutions of 2-D incompressible MHD system,J. Differential Equations,259(2015): 5440–5485
work page 2015
- [21]
- [22]
-
[23]
L. Ma, Global existence of three-dimensional incompressible magneto-micropolar system with mixed partial dissipation, magnetic diffusion and angular viscosity,Comput. Math. Appl.75 (2018), 170-186
work page 2018
-
[24]
D. Regmi and J. Wu, Global regularity for the 2D magneto-micropolar equations with partial dissipation,J. Math. Study,49(2016): 169–194
work page 2016
-
[25]
X. Ren, J. Wu, Z. Xiang and Z. Zhang, Global existence and decay of smooth solution for the 2-D MHD equations without magnetic diffusion,J. Funct. Anal.,267(2014): 503–541
work page 2014
-
[26]
X. Ren, Z. Xiang and Z. Zhang, Global well-posedness for the 2D MHD equations without magnetic diffusion in a strip domain,Nonlinearity,29(2016): 1257–1291
work page 2016
-
[27]
M.A. Rojas-Medar and J.L. Boldrini, Magneto-micropolar fluid motion: existence of weak solu- tions,Rev. Mat. Complut.,11(1998): 443–460
work page 1998
-
[28]
M. Sermange and R. Temam, Some mathematical questions related to the MHD equations, Comm. Pure Appl. Math.,36(1983): 635–664
work page 1983
-
[29]
H. Shang and J. Wu, Global regularity for 2D fractional magneto-micropolar equations,Math. Z.,297(2021), 775–802
work page 2021
-
[30]
H. Shang and J. Zhao, Global regularity for 2D magneto-micropolar equations with only micro- rotational velocity dissipation and magnetic diffusion,Nonlinear Anal.150(2017), 194-209
work page 2017
-
[31]
R.M. Strain and Y. Guo, Almost exponential decay near Maxwellian,Comm. PDE,31(2006): 417–429
work page 2006
- [32]
-
[33]
Temam,Navier-Stokes Equations: Theory and Numerical Analysis, North-Holland, Amster- dam, 1977
R. Temam,Navier-Stokes Equations: Theory and Numerical Analysis, North-Holland, Amster- dam, 1977
work page 1977
-
[34]
Y. Wang and Y. Wang, Blow-up criterion for two-dimensional magneto-micropolar fluid equa- tions with partial viscosity,Math. Meth. Appl. Sci.,34(2011): 2125–2135
work page 2011
-
[35]
Z. Xiang and H. Yang, On the regularity criteria for the 3D magneto-micropolar fluids in terms of one directional derivative,Bound. Value Prob.2012, 139 (2012)
work page 2012
- [36]
-
[37]
K. Yamazaki, Global regularity of the two-dimensional magneto-micropolar fluid system with zero angular viscosity,Disc. Cont. Dyna. Syst.,35(2015): 2193–2207
work page 2015
-
[38]
Yuan, Regularity of weak solutions to magneto-micropolar fluid equations,Acta Math
B. Yuan, Regularity of weak solutions to magneto-micropolar fluid equations,Acta Math. Sci- entia,30(2010), 1469-1480. 42 (Z. Feng)School of Mathematical Sciences, Chongqing Normal University, China Email address:zefufeng@mails.ccnu.edu.cn (G. Hong)School of Mathematics, South China University of Technology, China Email address:magyhong@scut.edu.cn (J. Wu...
work page 2010
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.