3D Magneto-Micropolar Equations

• 3D magneto-micropolar equations are coupled nonlinear PDEs modeling electrically conducting micropolar fluids with internal spin and magnetic stretching effects.
• They feature antisymmetric velocity-spin coupling and linear damping in the micro-rotation, enhancing dissipation and stabilizing fluid dynamics.
• The framework spans various domains and settings, integrating decay theory, Liouville theorems, and energy-stable numerical schemes to support robust analysis.

The 3D magneto-micropolar equations are coupled nonlinear partial differential equations for electrically conducting micropolar fluids in three spatial dimensions. In their standard incompressible form, the unknowns are the velocity $u$, the micro-rotation or angular-velocity field $\omega$, the magnetic field $b$, and the pressure $p$. Relative to classical magnetohydrodynamics, the system adds an internal rotational degree of freedom and an antisymmetric $u$–$\omega$ coupling, together with linear damping in the spin equation. Across recent work, the subject includes whole-space Cauchy problems, periodic and strip-domain formulations, stationary Liouville theory, partial regularity, pressure-based continuation criteria, energy-stable numerical schemes, compressible extensions, and application-driven models with explicit magnetization dynamics (Niche et al., 2020).

1. Canonical equations and principal variants

A common incompressible formulation on $\mathbb R^3$ is

$\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$

where $\mu>0$ is the kinematic viscosity, $\gamma>0$ the micro-rotation viscosity, $\omega$0 the magnetic diffusivity, and $\omega$1 the vortex-viscosity. The term $\omega$2 is a linear damping term and plays a decisive role in decay and stability theory (Niche et al., 2020).

The stationary equations are obtained by removing the time derivatives. In the normalization used by Zhang–Zu, one seeks smooth fields $\omega$3 on $\omega$4 satisfying

$\omega$5

with $\omega$6 (Zhang et al., 4 Aug 2025).

Generalized models replace the Laplacians in the $\omega$7, $\omega$8, and $\omega$9 equations by fractional powers $b$0 and $b$1, while retaining the divergence-free constraints and the damping term $b$2. Kim’s pressure-criterion study treats

$b$3

with corresponding fractional equations for $b$4 and $b$5, for $b$6 (Kim, 2023).

Other normalizations differ by constant factors but preserve the same structural triad: incompressible transport, magnetic stretching, and micro-rotation coupling. On $b$7, for example, Shang studies

$b$8

$b$9

$p$0

again with $p$1 (Shang, 3 Mar 2026). In strip-domain work without resistivity and spin viscosity, the $p$2-equation may reduce to an ODE-type damping-coupling law $p$3 (Zhao, 21 May 2026).

A standard structural observation is that suppressing the micro-rotation field recovers a classical MHD-type system. Several papers use this comparison to isolate what the microstructure contributes analytically: faster decay, effective dissipation, and in some regimes a stabilizing mechanism (Shang, 3 Mar 2026).

2. Geometric settings and solution frameworks

The theory is distributed across several domains and function-space regimes. For the whole-space nonstationary problem on $p$4, Niche–Perusato work with Leray-type weak solutions satisfying

$p$5

with initial data $p$6, $p$7 (Niche et al., 2020).

For stationary problems on $p$8, the usual hypothesis is full smoothness, $p$9, with decay or growth encoded by Lebesgue norms over expanding annuli $u$0. Zhang–Zu introduce the scale-invariant quantities

$u$1

which organize the Liouville theory in a form sensitive to scaling and annular growth (Zhang et al., 4 Aug 2025).

Periodic problems on $u$2 are typically posed for zero-mean Sobolev data. Shang assumes $u$3 or $u$4, with $u$5, and analyzes global strong solutions near the zero state or near a constant background magnetic field $u$6 (Shang, 3 Mar 2026).

Bounded-domain formulations split into two major classes. In numerical analysis, Qiu considers a bounded Lipschitz domain $u$7, homogeneous Dirichlet conditions for $u$8 and $u$9, and perfect-conductor conditions $\omega$0, $\omega$1 on $\omega$2 (Qiu, 2023). In the compressible Coulomb-force problem of Liu–Zhong, $\omega$3 is a bounded simply connected $\omega$4 domain, and slip conditions are imposed: $\omega$5 on $\omega$6 (Liu et al., 2022).

The strip domain $\omega$7 has become central in the non-resistive, non-spin-viscous theory. There the velocity satisfies the no-slip condition $\omega$8, while the magnetic and micro-rotation equations are handled through Lagrangian reformulations, div–curl structure, and refined trace estimates rather than standard dissipative boundary conditions (Zhao, 21 May 2026, Feng et al., 6 Jul 2026).

3. Sharp decay and asymptotic behavior

The sharp $\omega$9-decay theory for the 3D whole-space incompressible system is organized by the decay character $\mathbb R^3$0, which measures the lowest-frequency behavior of the initial datum through

$\mathbb R^3$1

For $\mathbb R^3$2, Niche–Perusato prove that if $\mathbb R^3$3 and $\mathbb R^3$4, then

$\mathbb R^3$5

and for $\mathbb R^3$6 there is a matching lower bound

$\mathbb R^3$7

The micro-rotation field decays faster: $\mathbb R^3$8 so one always has at least $\mathbb R^3$9-decay for $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$0, regardless of $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$1 (Niche et al., 2020).

That faster decay is not a secondary refinement but a direct consequence of the linear damping. Niche–Perusato also compare the nonlinear solution with the linearized flow

$\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$2

showing

$\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$3

and establishing first-order asymptotic estimates for $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$4 and $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$5. Their proof combines decay character, Fourier splitting, and Duhamel estimates exploiting the factor $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$6 (Niche et al., 2020).

Anisotropic dissipation produces a different but related decay theory. For the 3D system with horizontal dissipation, Zhang and collaborators establish global well-posedness when the initial data and their vertical derivatives are sufficiently small in $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$7. Under the additional assumption of anisotropic negative regularity $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$8, with $\begin{cases} \partial_t u + (u\cdot\nabla)u + \nabla p = (\mu+\chi)\,\Delta u + \chi\,\nabla\times\omega + (b\cdot\nabla)b,\[4pt] \partial_t \omega + (u\cdot\nabla)\omega = \gamma\,\Delta\omega + \nabla(\nabla\!\cdot\omega) + \chi\,\nabla\times u - 2\chi\,\omega,\[4pt] \partial_t b + (u\cdot\nabla)b = \nu\,\Delta b + (b\cdot\nabla)u,\[4pt] \nabla\!\cdot u=\nabla\!\cdot b=0, \end{cases}$9, the global solution satisfies

$\mu>0$0

$\mu>0$1

which are stated as the optimal large-time decay rates for the $\mu>0$2-norm of the solution (Lu et al., 23 Sep 2025).

4. Stability with partial or degenerate dissipation

Recent work has emphasized regimes in which classical viscous or magnetic smoothing is missing, and where the microstructure compensates for that loss. On $\mu>0$3, Shang studies two such cases. In the first, $\mu>0$4 but $\mu>0$5, $\mu>0$6, and $\mu>0$7; for small $\mu>0$8 data, there exists a unique global strong solution with

$\mu>0$9

In the second, $\gamma>0$0 and $\gamma>0$1, one perturbs around a constant background magnetic field $\gamma>0$2, assumes $\gamma>0$3, a Diophantine condition

$\gamma>0$4

and obtains a unique global solution with algebraic decay in lower norms, for example

$\gamma>0$5

The central conclusion is that the microstructure has the effect of enhancing dissipation and contributes to stabilize the fluid (Shang, 3 Mar 2026).

This point is especially striking because the corresponding MHD limits remain difficult. Shang explicitly notes that when $\gamma>0$6, the first regime becomes the 3D inviscid and resistive MHD equations, whose stability problem is still a challenging open problem, and the second becomes 3D ideal MHD, for which weighted spaces are usually required in global results (Shang, 3 Mar 2026). A plausible implication is that the $\gamma>0$7–$\gamma>0$8 coupling acts as an effective dissipation mechanism even when direct velocity viscosity or magnetic diffusion is absent.

The strip-domain theory develops the same theme with boundary effects. Zhao studies the incompressible system in $\gamma>0$9 without resistivity and spin viscosity,

$\omega$00

$\omega$01

and proves global well-posedness of small classical solutions by a two-layer energy method adapted from Guo–Tice. The resulting lower-order energy and dissipation satisfy

$\omega$02

while higher-order norms remain uniformly controlled (Zhao, 21 May 2026).

Luo and collaborators analyze a related strip problem with zero magnetic and angular viscosities near the steady state

$\omega$03

For sufficiently small perturbations in high Sobolev norms, they obtain a unique global classical solution and establish decay

$\omega$04

For $\omega$05, the decay is faster than any polynomial and is described as “almost exponential” (Feng et al., 6 Jul 2026).

5. Stationary equations and Liouville theory

The stationary 3D magneto-micropolar equations on $\omega$06 support a substantial Liouville theory, whose aim is to identify conditions under which every smooth stationary solution must be trivial. In the 2022 study of Cho, Neustupa, and Yang, the system is considered with $\omega$07, and three distinct criteria are proved. First, if there exist smooth matrix-valued potentials $\omega$08 with

$\omega$09

and their mean oscillations satisfy a BMO-type growth bound for some $\omega$10, then any smooth solution is trivial. Second, if $\omega$11 for some $\omega$12, then $\omega$13. Third, if $\omega$14, $\omega$15, and $\omega$16 as $\omega$17, then again the solution is trivial (Kim et al., 2022).

These theorems are based on a cut-off energy inequality, an auxiliary divergence-correction field $\omega$18 solving a Galdi-type divergence problem, and estimates of the form

$\omega$19

followed by an iteration lemma (Kim et al., 2022).

Zhang–Zu substantially sharpen this picture. Assuming $\omega$20, they prove Liouville theorems under annular growth conditions encoded by $\omega$21, $\omega$22, and $\omega$23. In the “large-$\omega$24” and “small-$\omega$25” regimes they obtain triviality under explicit algebraic relations among $\omega$26, $\omega$27, $\omega$28, and $\omega$29, and in the Lebesgue-space corollary they show $\omega$30 whenever $\omega$31 with one of the index ranges

$\omega$32

or

$\omega$33

or

$\omega$34

They also prove purely micropolar analogues for $\omega$35, with the sharp threshold $\omega$36 (Zhang et al., 4 Aug 2025).

The most distinctive feature of Zhang–Zu’s result is the treatment of the angular velocity. Compared with the velocity and magnetic fields, they raise the most relaxed restriction for the angular velocity, allowing the $\omega$37-norm on annuli to grow polynomially at any degree. Their proof combines a Bogovskii-map correction, an iterative energy-cutoff inequality, Gagliardo–Nirenberg interpolation for $\omega$38 and $\omega$39, and a novel two-case interpolation for $\omega$40, followed by Giaquinta’s lemma (Zhang et al., 4 Aug 2025). This directly contradicts the common intuition that $\omega$41 must satisfy decay assumptions comparable to those imposed on $\omega$42 and $\omega$43.

6. Regularity theory and numerical approximation

Local regularity theory for the nonstationary 3D equations has been developed in the suitable-weak-solution framework. In the perturbed system studied by Campanato-type methods, one introduces Elsässer variables $\omega$44, $\omega$45, a divergence-free perturbation field $\omega$46, external forces $\omega$47, and a local energy inequality for

$\omega$48

There exists an absolute $\omega$49 such that if

$\omega$50

then $\omega$51, $\omega$52, and $\omega$53 belong to $\omega$54 for some $\omega$55. An equivalent cubic criterion using

$\omega$56

is also available (Chamorro et al., 2021).

Kim’s pressure-based regularity criteria place continuation conditions directly on the pressure in Lorentz spaces. For Leray–Hopf weak solutions of the generalized fractional system, if $\omega$57 and either

$\omega$58

with the norm small, or

$\omega$59

then the weak solution is smooth on $\omega$60. In the MHD-type case $\omega$61, an additional criterion involving small Lorentz norms of $\omega$62 and $\omega$63 also yields smoothness up to time $\omega$64 (Kim, 2023).

The numerical analysis of the full 3D nonstationary equations has concentrated on preserving the underlying energy law. In a bounded Lipschitz domain, Qiu derives the continuous identity

$\omega$65

and constructs two families of unconditionally energy stable schemes: a first-order semi-implicit Euler discretization and a second-order Crank–Nicolson scheme with extrapolated nonlinear terms preserving skew-symmetry (Qiu, 2023).

Under standard regularity assumptions, the semi-implicit Euler method yields first-order accuracy $\omega$66 for $\omega$67, $\omega$68, $\omega$69, and the pressure in the stated norms, while the Crank–Nicolson–extrapolation scheme improves the time order to $\omega$70 with spatial order $\omega$71. Qiu also proposes fully decoupled first-order schemes based on Gauss–Seidel-type splitting, and reports numerical tests with manufactured solutions, monotone decay of the discrete energy $\omega$72, and a 3D lid-driven cavity computation at high Reynolds number (Qiu, 2023).

7. Compressible and application-driven extensions

The compressible theory couples magneto-micropolar dynamics with density evolution and, in some models, electrostatic effects. Jia–Tan–Zhou study the compressible magneto-micropolar Navier–Stokes–Poisson system on $\omega$73, with unknowns $\omega$74, and show global existence near a constant state under the smallness condition $\omega$75, even when higher Sobolev norms may be large. They further prove optimal decay rates in negative Sobolev or Besov scales,

$\omega$76

and obtain the usual $\omega$77-$\omega$78 decay without requiring the $\omega$79 norm of the initial data to be small (Jia et al., 2019).

Liu–Zhong treat a bounded-domain compressible isentropic system with Coulomb force and slip boundary conditions. If the initial total energy

$\omega$80

is sufficiently small, then there exists a unique global classical solution satisfying

$\omega$81

for every $\omega$82. Their proof relies on piecewise estimates, slip-boundary identities, and effective viscous fluxes involving both velocity and micro-rotational velocity (Liu et al., 2022).

A separate application-driven line introduces explicit magnetization and micromagnetorotation (MMR). Sabeel and Hameed formulate a 3D incompressible magneto-micropolar model with unknowns $\omega$83, including the MMR torque

$\omega$84

in the angular-momentum equation and an MMR constitutive law with new nondimensional parameters

$\omega$85

They then derive a boundary-layer reduction and study its ODE form numerically (Khan et al., 2023).

Aslani and collaborators use a related micropolar MHD blood-flow model with magnetization law

$\omega$86

and micromagnetic torque $\omega$87, and solve the 3D stenosed-artery problem numerically with the OpenFOAM solvers epotMicropolarFoam and epotMMRFoam. Their computations indicate that when MMR is disregarded, the magnetic field does not significantly alter blood flow regardless of its intensity, whereas with MMR the flow exhibits reductions of up to $\omega$88 in velocity and vorticity and up to $\omega$89 in microrotation, with damping of vortices and disturbances (Aslani et al., 18 Apr 2025).

Taken together, these extensions show that the 3D magneto-micropolar framework is not a single equation set but a family of closely related models. The core analytic features—antisymmetric velocity–spin coupling, magnetic stretching, and often a damping term in the micro-rotation equation—persist across incompressible, compressible, periodic, bounded, strip-domain, and magnetization-augmented formulations. The current literature suggests that these features are not merely constitutive embellishments: they materially alter decay, Liouville rigidity, regularity thresholds, and global stability mechanisms (Niche et al., 2020, Shang, 3 Mar 2026).