3D Stationary Magneto-Micropolar Fluid Equations

• 3D stationary magneto-micropolar fluid equations are steady PDE models that couple velocity, micro-rotation, and magnetic fields, extending Navier–Stokes and MHD theories.
• Recent studies use Liouville-type theorems with decay, integrability, and oscillation conditions to prove that smooth solutions must be trivial under certain regimes.
• Perturbative analyses in strip and periodic geometries demonstrate nonlinear stability around steady magnetic configurations, highlighting the role of microstructure effects.

Searching arXiv for recent and foundational papers on 3D stationary magneto-micropolar fluid equations. arXiv search query: "3D stationary magneto-micropolar fluid equations Liouville type theorems". The 3D stationary magneto-micropolar fluid equations are the steady incompressible PDEs that couple a velocity field, a micro-rotational or angular-velocity field, a magnetic field, and a pressure, thereby extending both stationary Navier–Stokes and stationary MHD models to fluids with microstructure. In the recent literature, the subject is organized around two principal themes: Liouville-type rigidity on $\mathbb{R}^3$, which identifies regimes in which every smooth stationary solution must be trivial, and nonlinear stability of distinguished steady magnetic configurations in strip or periodic geometries (Kim et al., 2022, Zhang et al., 4 Aug 2025, Feng et al., 6 Jul 2026, Shang, 3 Mar 2026).

1. Governing systems and model variants

A standard 3D stationary incompressible magneto-micropolar system on $\mathbb{R}^3$ is

$$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$

where $u$ is the fluid velocity, $w$ is the micro-rotational velocity, $b$ is the magnetic field, and $\Pi$ is the pressure; in one formulation the positive parameters $\gamma,\nu,\chi$ are set to $1$ for simplicity (Kim et al., 2022).

A closely related formulation retains an explicit micro-rotational viscosity parameter $\chi>0$ and writes the angular variable as $\mathbb{R}^3$0: $\mathbb{R}^3$1 Setting $\mathbb{R}^3$2 reduces this system to the stationary micropolar fluid equations (Zhang et al., 4 Aug 2025).

These equations encode three coupled mechanisms. The velocity equation combines viscous diffusion, inertial transport, magnetic self-interaction, and coupling to microrotation. The angular equation contains diffusion, a $\mathbb{R}^3$3 term, linear damping such as $\mathbb{R}^3$4 or $\mathbb{R}^3$5, and curl-coupling back to the fluid velocity. The magnetic equation is an elliptic stationary induction equation constrained by $\mathbb{R}^3$6. The stationary setting removes time evolution and therefore shifts analysis toward rigidity, decay, and structure at spatial infinity.

2. Liouville rigidity on $\mathbb{R}^3$7

A central question is whether smooth stationary solutions on all of $\mathbb{R}^3$8 can exist under natural decay, integrability, or oscillation hypotheses. The paper "Some Liouville-type theorems for the stationary 3D magneto-micropolar fluids" proves that, under several such hypotheses, the only smooth solution is

$\mathbb{R}^3$9

(Kim et al., 2022).

Three representative regimes are as follows:

Regime	Assumptions	Conclusion
Potential-growth criterion	Potentials $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$0 with $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$1, $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$2, $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$3 and mean oscillation growth bounded by $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$4 for $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$5	$$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$6
$$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$7 criterion	$$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$8 with $$\left\{ \begin{aligned} -\Delta u+(u \cdot \nabla) u + \nabla \Pi &= \nabla \times w + (b \cdot \nabla) b, \ -\Delta w+(u \cdot \nabla) w &= \nabla(\nabla \cdot w) + \nabla \times u - 2w, \ -\Delta b + (u\cdot \nabla)b &= (b\cdot \nabla)u, \ \nabla \cdot u &= 0,\quad \nabla \cdot b = 0, \end{aligned} \right.$$9	$u$0
Low-$u$1 plus decay	$u$2 with $u$3 and $u$4 at infinity	$u$5

The first result uses tensor potentials $u$6 satisfying

$u$7

together with the oscillation bound

$u$8

for all $u$9 and some $w$0 (Kim et al., 2022). This extends Liouville theory beyond plain Lebesgue assumptions to a mean-oscillation framework.

The same paper identifies $w$1 as the critical exponent coming from the scale-invariance of the system, and states that the inclusion of the micro-rotational vector $w$2 does not fundamentally change the scaling or the range of critical exponents for Liouville-type theorems (Kim et al., 2022). This directly addresses a common misconception that the micropolar coupling necessarily shifts the rigidity threshold.

The proofs combine cut-off function arguments, interpolation inequalities, and exploitation of the divergence-free structure and potential theory (Kim et al., 2022). In this framework, Liouville theorems are not merely uniqueness statements; they also exclude globally defined steady states with insufficient decay or overly controlled oscillation.

3. Annular growth criteria and relaxed control of angular velocity

A later development replaces global $w$3 assumptions by scale-sensitive conditions on annuli at infinity. The paper "Some new Liouville type theorems for the 3D stationary magneto-micropolar fluid equations" introduces

$w$4

and the scaled quantities

$w$5

Its main point is that the angular velocity can be treated much more flexibly than the velocity and magnetic fields (Zhang et al., 4 Aug 2025).

For $w$6, $w$7, $w$8, $w$9, $b$0, $b$1, and $b$2, if there exists $b$3 such that

$b$4

and either

$b$5

or

$b$6

then $b$7 (Zhang et al., 4 Aug 2025).

For $b$8, a parallel result assumes

$b$9

with either (A1) or

$\Pi$0

and again concludes $\Pi$1 (Zhang et al., 4 Aug 2025).

The paper emphasizes two improvements. First, it allows the $\Pi$2-norm of the angular velocity on the annuli to grow polynomially at any degree. Second, it improves the admissible integrability range for the magnetic field to $\Pi$3, from earlier $\Pi$4 (Zhang et al., 4 Aug 2025). The associated corollary gives Lebesgue-space vanishing criteria under the cases (B1), (B2), and (B3) stated there.

Methodologically, these results rely on an iterative procedure, the Bogovskii operator to handle pressure terms while preserving divergence-free constraints, a Caccioppoli-type inequality, Giaquinta's iteration lemma, and a novel combination of two different interpolation inequalities (Zhang et al., 4 Aug 2025). The flexible treatment of $\Pi$5 is the distinctive feature of this annular approach.

4. Steady states and nonlinear stability in bounded geometries

The stationary theory also includes perturbative analysis around explicit steady magnetic configurations. In a strip domain $\Pi$6, one paper studies the 3D incompressible magneto-micropolar equations without magnetic diffusivity and angular viscosity and focuses on the steady state

$\Pi$7

It proves that any small perturbation near this steady magnetic field perpendicular to the horizontal boundary leads to a unique global classical solution, and that the solution converges to the steady state at an almost exponential rate as time goes to infinity (Feng et al., 6 Jul 2026).

In the torus $\Pi$8, another line of work studies perturbations near a constant background magnetic field $\Pi$9. With zero kinematic viscosity but positive microstructure effects, it proves global stability and exponential decay in one partially dissipative regime, and global stability with algebraic decay in the case of zero kinematic viscosity and zero magnetic diffusion, provided $\gamma,\nu,\chi$0 and the initial perturbation is sufficiently small in high Sobolev norms (Shang, 3 Mar 2026). The same paper states that, if the micro-rotation effect is neglected, the system reduces to 3D inviscid and resistive MHD or to 3D ideal MHD, for which the corresponding stability problems remain open or require weighted spaces (Shang, 3 Mar 2026).

A related strip-domain result establishes global well-posedness of classical solutions for the 3D incompressible magneto-micropolar fluid system without resistivity and spin viscosity and shows algebraic time-decay of solutions toward the equilibrium (Zhao, 21 May 2026). That work identifies three major analytical obstacles: degeneracy from the lack of magnetic diffusion and spin viscosity, the non-dissipative anti-symmetric coupling between velocity and micro-rotation, and the interaction among velocity, magnetic field, and pressure when the pressure acts as a non-state variable (Zhao, 21 May 2026).

These results do not replace whole-space Liouville theorems. Rather, they address a different question: whether a particular stationary configuration is a nonlinear attractor under small perturbations. This distinction is essential in the modern theory.

5. Analytical structures and proof technology

The stationary Liouville literature is built around localization and iteration. The 2022 rigidity paper uses cut-off arguments, interpolation inequalities, divergence-free identities, and potential theory (Kim et al., 2022). The 2025 paper sharpens this by combining two interpolation inequalities within an iterative annular scheme, using the Bogovskii operator for pressure elimination under divergence constraints and an iterative Caccioppoli-type inequality closed by Giaquinta's lemma (Zhang et al., 4 Aug 2025).

The steady-state stability literature in domains with boundary uses a different toolkit. The strip-domain stability paper without magnetic and angular viscosities adopts a Lagrangian reformulation, a two-tier energy method inspired by Guo and Tice, Helmholtz decomposition for projected curls, and anisotropic div-curl and Stokes estimates adapted to strip geometry (Feng et al., 6 Jul 2026). A closely related well-posedness paper for the non-resistive, no-spin-viscosity case uses a multi-layer energy framework, refined trace inequalities, Stokes regularity, div-curl decomposition, and time-weighted energy inequalities yielding algebraic decay (Zhao, 21 May 2026).

A broader regularity backdrop is provided by the interior $\gamma,\nu,\chi$1-regularity theory for incompressible magneto-micropolar equations with a perturbation term. There, the key criterion is the smallness of the localized dissipation

$\gamma,\nu,\chi$2

which implies Hölder continuity near $\gamma,\nu,\chi$3 for suitable weak solutions (Chamorro et al., 2021). Although this is a time-dependent theory, it clarifies the role of localized dissipation in controlling singular behavior in magneto-micropolar systems.

6. Position within the broader magneto-micropolar literature

The stationary equations sit at the intersection of several PDE traditions. If the magnetic field is removed by setting $\gamma,\nu,\chi$4, the stationary magneto-micropolar system reduces to the stationary micropolar fluid equations (Zhang et al., 4 Aug 2025). If the micro-rotation field is neglected, the problem reduces toward classical MHD or Navier–Stokes settings, and several papers explicitly compare the resulting stability questions with open problems in inviscid resistive MHD or ideal MHD (Shang, 3 Mar 2026).

A second recurrent theme is the role of microstructure. In the Liouville setting, the 2022 whole-space results state that the inclusion of the micro-rotational vector does not fundamentally change the scaling or the range of critical exponents (Kim et al., 2022). In contrast, the perturbative stability literature emphasizes that microstructure can enhance dissipation and stabilize electrically conducting fluids, even in regimes where the corresponding MHD problems are open (Shang, 3 Mar 2026). This suggests that scaling-critical rigidity and nonlinear stability probe different aspects of the coupling.

The dynamical literature also shows why stationary states matter. Sharp decay estimates for the 3D magneto-micropolar system establish that the micro-rotational field decays faster because of a linear damping term, and that the nonlinear solution is asymptotically governed by the linearized dynamics (Niche et al., 2020). A plausible implication is that stationary and near-stationary configurations are natural asymptotic reference states, even when the primary analysis is time-dependent.

From the standpoint of mathematical scope, the currently documented stationary theory is strongest in two directions: global Liouville rigidity on $\gamma,\nu,\chi$5 under quantitative growth or integrability hypotheses, and small-perturbation stability of explicit steady magnetic fields in periodic or strip geometries. What remains difficult is precisely what several of these papers identify as difficult: weak dissipation, anti-symmetric velocity–microrotation coupling, boundary pressure effects, and the lack of magnetic diffusion or angular viscosity in genuinely three-dimensional settings (Feng et al., 6 Jul 2026, Zhao, 21 May 2026).