Incompressible Navier–Stokes–Maxwell System

• The incompressible Navier–Stokes–Maxwell system is a coupled set of partial differential equations that integrates incompressible fluid dynamics with electromagnetic fields through Ohm’s law.
• It leverages normalized formulations and energy methods to establish well-posedness and convergence, highlighting the interplay between parabolic (fluid) and hyperbolic-damped (electromagnetic) components.
• Applications include analysis of singular limits, kinetic derivations, and transitions to magnetohydrodynamic systems, providing deep insights into fluid-electromagnetic interactions.

Searching arXiv for recent and foundational papers on the incompressible Navier–Stokes–Maxwell system and closely related limits. arXiv Search Tool unavailable in this environment. Proceeding with the supplied arXiv records and citing them directly. The incompressible Navier–Stokes–Maxwell system denotes a class of coupled PDEs in which an incompressible viscous charged fluid is linked to Maxwell’s equations through Ohm’s law. In the classical formulation, the unknowns are the fluid velocity $u$, pressure $p$, electric field $E$, magnetic field $B$, and current $j$, and the model combines a parabolic incompressible Navier–Stokes equation with a hyperbolic–damped Maxwell subsystem. Across the literature, the term also borders on closely related systems: inhomogeneous variable-density versions, two-fluid formulations, thermally extended Navier–Stokes–Fourier–Maxwell models, and singular limits toward magnetohydrodynamics or toward solenoidal-Ohm-law reductions (Arsénio et al., 2018, Germain et al., 2012).

1. Canonical equations and constitutive structure

A classical incompressible Navier–Stokes–Maxwell system in space dimension $d=2$ or $3$ is written as

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$

Here $u$ is the fluid velocity, $E$ the electric field, $p$0 the magnetic field, $p$1 the pressure, $p$2 the current, $p$3 the speed of light, $p$4 the viscosity, and $p$5 the conductivity. In this formulation the current is not an independent unknown: it is determined by Ohm’s law. In two dimensions, the fields are embedded into $p$6 with

$p$7

“Incompressible” means $p$8, and the magnetic field is divergence-free as well, $p$9 (Arsénio et al., 2018).

A normalized version, used in wellposedness theory, sets the viscosity and conductivity equal to $E$0. The system then takes the form

$E$1

This representation makes explicit that the Lorentz forcing in the fluid equation is $E$2, while the damping term $E$3 belongs to the Maxwell operator itself (Germain et al., 2012).

A distinct limiting constitutive law appears in the solenoidal-Ohm-law reduction. In that regime the current is given by

$E$4

with $E$5, so that the scalar $E$6 enforces the solenoidal constraint on the current rather than on the velocity field alone (Guo et al., 16 Jul 2025).

2. Energy structure and analytical difficulties

The natural energy level for the classical system is

$E$7

together with the energy inequality

$E$8

where

$E$9

This energy level is natural but not by itself sufficient for general weak compactness of the nonlinear Lorentz term $B$0: a central obstacle is that the magnetic field lacks compactness, so finite-energy weak approximations do not permit straightforward passage to the limit in the coupling (Arsénio et al., 2018).

This hyperbolic–parabolic interaction dictates much of the modern analysis. One strategy is to propagate extra regularity of $B$1 in $B$2, which yields compactness for $B$3 and thereby permits passage to $B$4. In the large-energy weak-solution theory, this is coupled to refined energy estimates for Maxwell’s system, a Grönwall-like bootstrap, and a maximal parabolic estimate for the heat equation in Besov spaces strong enough to avoid Chemin–Lerner spaces altogether (Arsénio et al., 2018).

The damped character of the Maxwell subsystem is equally important in semigroup-based wellposedness theory. In the normalized formulation, the Maxwell equations imply dissipation for the electric field, and the magnetic field can be rewritten as a damped wave equation,

$B$5

which yields additional decay estimates for $B$6. This structure complements the parabolic smoothing of the Navier–Stokes component and underlies mild-solution frameworks based on the coupled linear semigroup (Germain et al., 2012).

3. Well-posedness theory and solution classes

The well-posedness theory separates sharply by dimension, regularity class, and whether the density is constant or variable. For the homogeneous incompressible system, a mild-solution theory establishes local existence for arbitrarily large initial data in spaces similar to the scale-invariant spaces classically used for Navier–Stokes, and global existence for sufficiently small data. In dimension $B$7, the theorem is formulated for

$B$8

while in dimension $B$9 it is stated for

$j$0

The proof uses semigroup bounds, Littlewood–Paley theory, Besov and logarithmically weighted spaces, product estimates based on Bony’s decomposition, and a contraction argument for the mild formulation

$j$1

with $j$2 (Germain et al., 2012).

A second line of results works at the energy level. In three dimensions, global weak solutions exist under the assumption that the initial velocity and electromagnetic fields have finite energy and that the initial electromagnetic field is small in $j$3 for $j$4. More precisely, if

$j$5

and

$j$6

then there exists a global weak solution satisfying the energy inequality and

$j$7

In two dimensions, the same analysis yields global weak solutions with

$j$8

for $j$9, together with estimates uniform in the speed of light $d=2$0 (Arsénio et al., 2018).

For the inhomogeneous planar system, the unknown density $d=2$1 satisfies the transport equation

$d=2$2

and the momentum equation becomes

$d=2$3

Global energy solutions are obtained in two dimensions when the initial density is bounded pointwise, bounded away from $d=2$4, and close to the constant state $d=2$5; the velocity satisfies

$d=2$6

and the electromagnetic field obeys

$d=2$7

The resulting solution is global and uniformly bounded with respect to the speed of light $d=2$8; if $d=2$9, then

$3$0

and the energy solution is unique in the class of all energy solutions (Arsénio et al., 2024).

4. Singular limits and reduced equations

Several asymptotic regimes reorganize the incompressible Navier–Stokes–Maxwell dynamics into reduced systems. One such limit starts from a two-fluid incompressible Navier–Stokes–Maxwell system for cation and anion velocities $3$1, and introduces the bulk velocity and current variables

$3$2

The reformulated system contains the equations

$3$3

$3$4

coupled to Maxwell’s equations and the divergence constraints. As $3$5, the $3$6-terms vanish and one obtains the incompressible Navier–Stokes–Maxwell system with solenoidal Ohm’s law,

$3$7

together with

$3$8

The notable feature is that the convergence is strong and occurs without loss of regularity, through a frequency-envelope method combined with a low–high frequency decomposition (Guo et al., 16 Jul 2025).

A second family of limits concerns the speed of light. In the two-dimensional homogeneous case, the global weak-solution theory yields estimates uniform in $3$9, and along subsequences the solutions converge as $$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$0 to the two-dimensional magnetohydrodynamic system

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$1

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$2

The proof uses compactness for both $$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$3 and $$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$4, with the magnetic compactness relying on the stronger uniform control obtained for the electromagnetic field (Arsénio et al., 2018).

The same nonrelativistic passage appears in the inhomogeneous planar system, where uniform-in-$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$5 bounds imply convergence toward the inhomogeneous MHD equations

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$6

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$7

$$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$8

Relative compactness follows from the energy and regularity bounds together with Aubin–Lions compactness (Arsénio et al., 2024).

A different singular regime begins with the full compressible Navier–Stokes–Maxwell system in a bounded smooth domain $$\left\{ \begin{aligned} &\partial_t u + u\cdot\nabla u - \mu\Delta u = -\nabla p + j\times B,\qquad \operatorname{div}u=0,\ &\frac1c\partial_t E - \nabla\times B = -j,\qquad j=\sigma(cE+u\times B),\ &\frac1c\partial_t B + \nabla\times E =0,\qquad \operatorname{div}B=0. \end{aligned} \right.$$9, with small Mach number and dielectric constant taken equal,

$u$0

Under well-prepared initial data and slip/electromagnetic perfect-conductor boundary conditions, uniform-in-$u$1 strong-solution estimates imply

$u$2

and

$u$3

The limit $u$4 solves the incompressible MHD system

$u$5

$u$6

$u$7

This shows that simultaneous removal of the acoustic scaling and the displacement current leads from compressible Navier–Stokes–Maxwell to incompressible MHD in a bounded domain (Fan et al., 2015).

5. Kinetic derivations and the Fourier–Maxwell extension

A substantial part of the recent theory derives incompressible fluid–Maxwell systems from kinetic equations. The macroscopic limit obtained in these works is usually the two-fluid incompressible Navier–Stokes–Fourier–Maxwell system with Ohm’s law rather than the bare incompressible Navier–Stokes–Maxwell system. Its core structure contains the incompressible momentum equation, Maxwell equations, and Ohm-type closure, together with a temperature equation and Boussinesq-type constraints (Jiang et al., 2019, Jiang et al., 2020).

For the two-species Vlasov–Maxwell–Boltzmann system under diffusive scaling, one rigorous limit yields

$u$8

$u$9

and an Ohm law of the form

$E$0

supplemented by

$E$1

One classical-solution result establishes global-in-time solutions uniform in the Knudsen number $E$2, proves convergence of the fluctuations, and shows that the microscopic part $E$3 vanishes strongly in $E$4 (Jiang et al., 2019). A companion Hilbert-expansion result constructs

$E$5

and derives a global-in-time remainder estimate uniform in $E$6 (Jiang et al., 2020).

These kinetic derivations have been extended beyond hard-sphere Boltzmann kernels. For soft potentials, uniform-in-$E$7 global classical estimates are proved for both non-cutoff and cutoff two-species Vlasov–Maxwell–Boltzmann systems. In the non-cutoff case the analysis assumes

$E$8

while the cutoff case considers

$E$9

The hydrodynamic limit again produces an incompressible Navier–Stokes–Fourier–Maxwell system with Ohm’s law. The same source explicitly notes that, for a Navier–Stokes–Maxwell-only discussion, the key equations are the momentum equation, incompressibility, Maxwell equations, and Ohm’s law, whereas the temperature equation is an extra feature of the Fourier extension (Jiang et al., 2023).

An analogous program has also been carried out for the two-species Vlasov–Maxwell–Landau system with Coulomb potential. There, a global-in-time uniform-in-Knudsen-number estimate shows that the microscopic part is dissipated at order $p$00, and the limit is once again the two-fluid incompressible Navier–Stokes–Fourier–Maxwell system with Ohm’s law, with transport coefficients $p$01, $p$02, and $p$03 determined by Landau-collision correctors (Lei et al., 2023).

6. Terminological boundaries and adjacent formulations

A recurrent source of confusion is the phrase “Navier–Stokes–Maxwell.” One important counterexample is the incompressible Navier–Stokes–Maxwell–Stefan system for multicomponent gaseous mixtures. That model consists of the species continuity equations,

$p$04

the Maxwell–Stefan diffusion laws,

$p$05

the Navier–Stokes momentum equation,

$p$06

and the incompressibility constraint $p$07. The corresponding 2018 paper develops an artificial compressibility approximation for this coupled fluid–diffusion model and proves existence of weak solutions for the approximating system together with convergence back to the incompressible Maxwell–Stefan limit as $p$08. It is explicitly not about the electromagnetic Navier–Stokes–Maxwell system (Dolce et al., 2018).

A second adjacent body of work comes from gauge/gravity duality. In Einstein–Maxwell and Gauss–Bonnet–Maxwell theories, a non-relativistic long-wavelength expansion around a charged AdS black brane and a timelike cutoff surface $p$09 yields the incompressible Navier–Stokes equation with external force density,

$p$10

while the Maxwell equations enforce

$p$11

This construction produces a forced incompressible Navier–Stokes equation on the cutoff surface rather than the full electromagnetic Navier–Stokes–Maxwell system with dynamical $p$12 and $p$13 fields (Niu et al., 2011).

Within these boundaries, the incompressible Navier–Stokes–Maxwell system is best viewed as a family of incompressible fluid–electromagnetic models organized by three structural choices: the Ohmic closure for $p$14, the regularity class in which the Lorentz term $p$15 can be controlled, and the asymptotic regime under consideration. The current literature therefore treats the system not as a single fixed PDE, but as a central electromagnetic fluid framework linked upward to two-fluid and kinetic descriptions and downward to MHD-type limits.