3D Stationary Micropolar Fluid Equations

• 3D stationary micropolar fluid equations are a continuum model coupling velocity and angular momentum via microrotation fields, generalizing the Navier–Stokes system.
• Analytic techniques such as localized energy estimates, Sobolev regularity, and annular interpolation are key to establishing existence, regularity, and Liouville-type uniqueness results.
• The derivation of explicit Green-type solutions, like micropolar Stokeslets and Oseenlets, underscores the distinctions from classical Newtonian flows.

Searching arXiv for papers on 3D stationary micropolar fluid equations and related Liouville, existence, and fundamental-solution results. The three-dimensional stationary micropolar fluid equations are a steady incompressible continuum model in which the usual macroscopic velocity field is coupled to a microrotation, or angular-velocity, field. In the cited literature, the system is treated on the whole space $\mathbb{R}^3$, with unknowns typically denoted by $u$ for velocity, $\omega$ or $w$ for microrotation, and $p$ or $\pi$ for pressure. It is studied both as a particular generalization of the usual Navier–Stokes system where the microrotations of the fluid particles must be taken into account, and as the zero-magnetic-field reduction of the stationary magneto-micropolar system (Chamorro et al., 2023, Kim et al., 2022).

1. Governing formulations

A forced stationary micropolar system used in the whole-space analysis is

$$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$

where $A=-\Delta$, $V=\nabla\times$, $u$0 is a given viscosity parameter, $u$1 are distributional forces, $u$2 is the fluid velocity, $u$3 is the pressure, and $u$4 is the microrotation field (Chamorro et al., 2023).

A normalized unforced form considered in recent Liouville theory is

$u$5

with $u$6 the vortex-viscosity parameter (Zhang et al., 4 Aug 2025).

The stationary micropolar system also appears as the pure-fluid reduction of the stationary three-dimensional magneto-micropolar system. If one formally sets the magnetic field $u$7, then

$u$8

which is precisely the classical stationary micropolar fluid system (Kim et al., 2022).

Across these formulations, the structural additions relative to stationary Navier–Stokes are the second vector field $u$9, the coupling terms $\omega$0 in the momentum equation and $\omega$1 in the microrotation equation, and the terms $\omega$2 and $\omega$3, or their equivalent normalized versions, in the angular-momentum balance (Chamorro et al., 2023).

2. Whole-space setting and solution classes

The stationary theory summarized in the cited works is formulated on $\omega$4, with no boundary. Decay at infinity is encoded either by Sobolev and Lebesgue membership or by explicit vanishing assumptions (Chamorro et al., 2023).

For weak-solution existence, one seeks

$\omega$5

under

$\omega$6

while strong solutions occur when $\omega$7 are smooth or zero, and one shows ultimately $\omega$8 (Chamorro et al., 2023).

For Liouville-type analysis, the principal setting is that of smooth entire solutions. One recent framework imposes

$\omega$9

and measures decay or growth at infinity by Lebesgue norms on dyadic annuli

$w$0

through the quantities

$w$1

This annular formulation is central in the 2025 Liouville results (Zhang et al., 4 Aug 2025).

A second whole-space framework uses potentials for the fields. In the magneto-micropolar setting, smooth solutions $w$2 are assumed to admit matrix-valued potentials $w$3 with

$w$4

together with a mean-oscillation bound on balls $w$5. This provides a BMO-type route to rigidity (Kim et al., 2022).

3. Existence and regularity

For the forced stationary problem in $w$6, there exists at least one weak solution $w$7 whenever $w$8 and $w$9, with

$p$0

(Chamorro et al., 2023).

The construction proceeds by mollifying the system with a cutoff $p$1 localizing to $p$2 and an extra viscous term $p$3, rewriting the mollified problem as a fixed-point equation

$p$4

in $p$5, where $p$6 with the divergence-free constraint on $p$7. The cited proof checks continuity, compactness, and a uniform bound for all $p$8-scale-invariant solutions $p$9, applies Schaefer’s theorem, and then passes $\pi$0 and $\pi$1 by weak and strong compactness (Chamorro et al., 2023).

The same work gives a regularity bootstrap. From the equation

$\pi$2

one recovers $\pi$3 by product laws in Sobolev spaces. Applying the Leray projector to the momentum equation yields

$\pi$4

so $\pi$5; then the $\pi$6-equation yields $\pi$7, and iteration gives arbitrarily high Sobolev regularity, whence smoothness (Chamorro et al., 2023).

Within the summarized literature, these existence and regularity results are distinct from the Liouville theorems. The former guarantee solvability for forced data in natural Sobolev classes, while the latter identify conditions under which unforced entire stationary solutions must collapse to the trivial state.

4. Liouville-type rigidity

A central topic is the absence of nontrivial entire stationary solutions under suitable integrability, decay, or growth restrictions.

For the stationary magneto-micropolar system in $\pi$8, smooth solutions are trivial under three classes of assumptions: first, a mean-oscillation growth condition on potentials $\pi$9 for some exponent $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$0; second, global $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$1 membership of $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$2 for $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$3; and third, global $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$4 membership for $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$5 together with

$$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$6

In each case one obtains

$$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$7

By setting $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$8, the same statements apply to the pure stationary micropolar system (Kim et al., 2022).

For the stationary micropolar system itself, Chamorro–Llerena–Vergara-Hermosilla prove a Liouville-type uniqueness statement for weak solutions: if $$A\,u-(u\cdot\nabla)u-\nabla p+V\,\omega+f=0,\qquad \operatorname{div}u=0,$$9 and

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$0

then the only solution is

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$1

Remarkably, no smallness or extra condition on $$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$2 is required—only a single Lebesgue-decay condition on $$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$3 (Chamorro et al., 2023).

A later annular-norm formulation sharpens the picture. For $$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$4, if

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$5

and there exists a sequence $$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$6 such that

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$7

then the only smooth solution is trivial. If

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$8

and

$$A\,\omega+\nabla\operatorname{div}\omega-(u\cdot\nabla)\omega-K\,\omega+V\,u+g=0,$$9

with the endpoint condition specified in the statement, then again $A=-\Delta$0. The corresponding corollary states that if

$A=-\Delta$1

then any smooth solution of the stationary micropolar system must be trivial (Zhang et al., 4 Aug 2025).

A recurrent issue in this literature is whether the additional microrotation field permits new entire steady states in $A=-\Delta$2. The cited Liouville theorems show that, under the stated whole-space decay, integrability, or mean-oscillation hypotheses, it does not.

5. Analytical methods

The proofs of the rigidity results follow a localized energy strategy. In the magneto-micropolar setting, one chooses a large ball $A=-\Delta$3, a smooth cutoff $A=-\Delta$4, and tests the equations by

$A=-\Delta$5

where $A=-\Delta$6 solves

$A=-\Delta$7

in an annular shell. Integration by parts cancels the pressure term through the divergence correction, and the nonlinear and coupling terms are estimated by Hölder and Young inequalities. One then applies an iteration lemma of Giaquinta type to an inequality for the localized Dirichlet integral, and finally lets $A=-\Delta$8 (Kim et al., 2022).

Two auxiliary tools are singled out in that analysis. One is a solenoidal extension lemma: if $A=-\Delta$9 has zero mean, then there exists $V=\nabla\times$0 with $V=\nabla\times$1 and $V=\nabla\times$2. The other is a mean-oscillation estimate for vector fields possessing a potential $V=\nabla\times$3 whose oscillation satisfies the ball-growth condition

$V=\nabla\times$4

These tools convert potential growth assumptions into $V=\nabla\times$5 and $V=\nabla\times$6 control for the fields (Kim et al., 2022).

In the 2025 micropolar analysis, the method is again a localized energy-iteration scheme, but with a Bogovskiĭ corrector on thin annuli and a specifically defined energy profile

$V=\nabla\times$7

From testing by $V=\nabla\times$8 and $V=\nabla\times$9, one obtains

$u$00

The proof then relies on a novel combination of two different interpolation inequalities to estimate the microrotation and mixed terms. In particular, the paper states that it allows the $u$01-norm of the angular velocity on the annuli to grow polynomially at any degree (Zhang et al., 4 Aug 2025).

These methods collectively illustrate a characteristic analytic pattern for stationary micropolar systems: pressure elimination by divergence correction, localized energy testing, annular interpolation, and an iteration lemma in the spirit of Giaquinta.

6. Linearized steady singular solutions and broader significance

The steady theory also includes explicit fundamental solutions for linearized three-dimensional incompressible micropolar flow. In the Oseen setting with a uniform background velocity $u$02, the linearized equations are

$u$03

$u$04

$u$05

with $u$06 a point force and $u$07 a point couple. In the pure Stokes limit $u$08, the convective terms drop out (Shu et al., 2014).

Explicit three-dimensional steady Green-type singularities derived in this framework include the micropolar Stokeslet, the micropolar Stokes couplet, the micropolar Oseenlet, and the micropolar Oseen couplet. These fundamental solutions do not exist in Newtonian flow due to the absence of microrotation velocity field. They satisfy the equations away from the origin, reproduce the required delta-singularity, and in the Stokes case decay at least as fast as $u$09 at infinity (Shu et al., 2014).

This linearized theory complements the nonlinear Liouville results. The Green-function constructions demonstrate that steady micropolar flow admits distinctive singular responses to point forcing and point couples, while the whole-space Liouville theorems show that smooth entire solutions with suitable decay or integrability are forced to be trivial. A plausible implication is that the stationary three-dimensional micropolar theory separates sharply between singularly forced linear responses and globally decaying nonlinear equilibria.

The cited works also make the comparison with classical stationary Navier–Stokes explicit. The micropolar system augments Navier–Stokes by the angular-velocity equation and coupling terms, yet the whole-space rigidity theorems show that the microrotation equation does not introduce new non-trivial entire-space steady solutions under the same or closely related decay ranges (Chamorro et al., 2023). In that sense, the stationary micropolar equations retain the classical Liouville theme while introducing substantially richer PDE structure.