---
title: Discrete Magnetic Advection Operator
url: https://www.emergentmind.com/topics/discrete-magnetic-advection-operator
type: topic
---

# Discrete Magnetic Advection Operator

Searching arXiv for the cited papers to ground the article in current records.
arxiv_search.query({"search_query":"id:1309.4398 OR id:1102.0494 OR id:2412.20106 OR id:1207.2955 OR id:1408.1346 OR id:1707.03227 OR id:1708.05057 OR id:2412.01927 OR id:2509.09913","start":0,"max_results":10})
arxiv_search.query({"query":"id:1309.4398 OR id:1102.0494 OR id:2412.20106 OR id:1207.2955 OR id:1408.1346 OR id:1707.03227 OR id:1708.05057 OR id:2412.01927 OR id:2509.09913","max_results":10})
arxiv_search({"query":"1309.4398","max_results":5})
A discrete magnetic advection operator is a numerical representation of the transport part of magnetic evolution equations. In the literature surveyed here, that role appears in several mathematically distinct but structurally related forms: as the transport term $(\mathbf{u}\cdot\nabla)\mathbf{B}$ in magnetic induction equations, as the curl form $\nabla\times(\mathbf{u}\times\mathbf{B})$ advanced by constrained transport, as the field-aligned operator $\nabla_\parallel=\mathbf{b}\cdot\nabla$ in magnetized plasma models, as the Lie derivative $-\boldsymbol{\beta}\times(\nabla\times\mathbf{u})+\nabla(\boldsymbol{\beta}\cdot\mathbf{u})$ in $\mathbf{H}(\mathrm{curl})$ formulations, and, in one analytic reformulation, as an effective advection generated by gradients of magnetic diffusivity [1309.4398]. The resulting discrete operators differ in placement of unknowns, conservation principles, and stabilization strategy, but they are unified by the requirement that magnetic transport be represented without destroying the geometric and energetic structure of the underlying PDEs.

## 1. Continuous operators and analytical origins

In resistive MHD with spatially varying magnetic diffusivity $\eta(\mathbf{x})$, the induction equation is
\[
\frac{d\mathbf{B}}{dt}
=
\nabla\times(\mathbf{v}\times\mathbf{B})
-
\nabla\times\bigl(\eta\,\nabla\times\mathbf{B}\bigr).
\]
Zita showed that, under $\nabla\cdot\mathbf{B}=0$, isotropic scalar $\eta$, and a first-order approximation in $\nabla\eta$, the diffusivity term can be decomposed as
\[
\nabla\times\left[\eta\,\nabla\times\mathbf{B}\right]
=
\nabla'\times\nabla'\times\mathbf{B}
-
(\mathbf{U}_\eta\cdot\nabla')\mathbf{B},
\]
with
\[
\nabla'=\sqrt{\eta}\,\nabla,
\qquad
\mathbf{U}_\eta=\nabla\sqrt{\eta},
\]
so that the induction equation becomes
\[
\frac{d\mathbf{B}}{dt}
=
\nabla\times(\mathbf{v}\times\mathbf{B})
+
(\mathbf{U}_\eta\cdot\nabla')\mathbf{B}
-
\nabla'\times\nabla'\times\mathbf{B}.
\]
Equivalently,
\[
(\mathbf{U}_\eta\cdot\nabla')\mathbf{B}
=
\frac{1}{2}(\nabla\eta\cdot\nabla)\mathbf{B}.
\]
This identifies an advection-like operator induced purely by diffusivity gradients, even when the physical flow $\mathbf{v}$ is zero [1309.4398].

In magnetic induction equations with prescribed velocity field, another standard continuous form is the non-conservative symmetric representation
\[
\partial_t \mathbf{B}+(\mathbf{u}\cdot\nabla)\mathbf{B}=M(D)\mathbf{B},
\]
or, in two dimensions,
\[
\mathbf{B}_t+\Lambda_1\mathbf{B}_x+\Lambda_2\mathbf{B}_y-C\mathbf{B}=0.
\]
Here the magnetic advection part is
\[
(\mathbf{u}\cdot\nabla)\mathbf{B}=u^1\mathbf{B}_x+u^2\mathbf{B}_y,
\]
while $C\mathbf{B}$ accounts for stretching and rotation [1102.0494].

In magnetized plasma models, the central transport operator is often the parallel gradient
\[
\nabla_\parallel f=\mathbf{b}\cdot\nabla f,
\qquad
\nabla\cdot\mathbf{b}=0.
\]
Because $\nabla\cdot\mathbf{b}=0$, one has
\[
\nabla\cdot(\mathbf{b}f)=\mathbf{b}\cdot\nabla f,
\]
and therefore
\[
\int_\Omega p\,\nabla_\parallel q\,d\Omega+\int_\Omega q\,\nabla_\parallel p\,d\Omega
=
\int_{\partial\Omega}pq\,\mathbf{b}\cdot\mathbf{n}\,ds.
\]
Under homogeneous Dirichlet conditions, the right-hand side vanishes, giving the skew-symmetry that many discrete magnetic advection operators seek to mimic [2412.20106].

A further continuous representation arises in formulations using the magnetic vector potential. Expanding $\mathbf{u}\times\mathbf{B}$ with $\mathbf{B}=\nabla\times\mathbf{A}$ gives
\[
\partial_t A_i
=
-
u_j\partial_j A_i
+
u_j\partial_i A_j
-
\eta\mu_0 J_i.
\]
The first term is a pure advection term on $\mathbf{A}$; the second is a shear term. In orbital-advection algorithms this splitting is exploited directly [1708.05057].

## 2. Principal discrete formulations

The surveyed literature does not present a single canonical discrete magnetic advection operator. Instead, several operator classes appear, each tied to a specific state variable and discretization framework.

| Formulation | Primary transported quantity | Discrete magnetic advection operator |
|---|---|---|
| SBP-SAT finite difference | $\mathbf{B}$ on tensor grids | $\mathcal{A}_h(\mathbf{u})V=u^1\circ d_xV+u^2\circ d_yV$ |
| Mimetic staggered-grid transport | Scalars on staggered $p$/$q$ grids | $C_{pq}=\alpha A_{pq}+(1-\alpha)B_{pq}$, $C_{qp}=\beta A_{qp}+(1-\beta)B_{qp}$ |
| CT orbital advection | Face-centered $\mathbf{B}$ | EMF-based linear transport step combined with CT update |
| Vector-potential orbital advection | Point-collocated $\mathbf{A}$ | Explicit residual advection plus implicit spectral shift by $\mathcal{L}=\bar u_\phi r^{-1}\partial_\phi$ |
| Variational DEC / Lagrangian labeling | Discrete magnetic flux cochains | Built-in advection: $\langle B,\sigma^2\rangle=\langle B_0,\sigma_0^2\rangle$ |
| $\mathbf{H}(\mathrm{curl})$ SUPG | $\mathbf{u}$, later identified with $\mathbf{A}$ | $\tilde L_{\boldsymbol{\beta},h}\mathbf{u}_h=L_{\boldsymbol{\beta},h}\mathbf{u}_h-\mathbf{r}_\alpha(\boldsymbol{\beta}\cdot\mathbf{n}\,\llbracket\mathbf{u}_h\rrbracket)$ |

In the SBP-SAT framework for the magnetic induction equations, the discrete magnetic advection operator is defined componentwise by
\[
\boxed{\mathcal{A}_h(\mathbf{u})V=u^1\circ d_xV+u^2\circ d_yV,}
\]
with $d_x,d_y$ the 2D SBP derivative matrices and $\circ$ the Hadamard product. If both field components are stacked, $\mathcal{A}_h(\mathbf{u})$ becomes block diagonal, with identical scalar advection operators on each block. The operator is central in the interior; boundary stability is supplied by SAT penalties rather than by upwind bias [1102.0494].

For divergence-free field-aligned transport, the mimetic finite difference construction uses two staggered grids. The discrete parallel gradient from the $q$-grid to the $p$-grid is
\[
\nabla_\parallel|_{pq}f
=
\alpha\,[\mathbf{b}\cdot\nabla|_{pq}f]
+
(1-\alpha)\,[\nabla|_{pq}\cdot(\mathbf{b}f)],
\]
and analogously from $p$ to $q$,
\[
\nabla_\parallel|_{qp}f
=
\beta\,[\mathbf{b}\cdot\nabla|_{qp}f]
+
(1-\beta)\,[\nabla|_{qp}\cdot(\mathbf{b}f)].
\]
In matrix form these are $C_{pq}$ and $C_{qp}$. When $\alpha+\beta=1$, they satisfy
\[
C_{pq}=-C_{qp}^T,
\]
which is the discrete skew-adjointness condition underpinning energy preservation [2412.20106].

In CT orbital-advection schemes for ideal MHD, the discrete magnetic advection operator is not written as a single matrix acting on $\mathbf{B}$. Instead it is assembled from time-integrated EMFs and discrete Stokes updates on face-centered magnetic fields. During the orbital step with $\mathbf{w}=w\,\mathbf{e}_y$,
\[
\partial_t\mathbf{B}=\nabla\times(\mathbf{w}\times\mathbf{B}),
\]
and EMFs such as
\[
\mathcal{E}^x=-\int wB_z\,dt,
\qquad
\mathcal{E}^z=\int wB_x\,dt
\]
are computed by an integer shift plus fractional flux construction. These EMFs are then inserted into CT face updates. This defines a linear mapping
\[
\mathbf{B}^{n+1}=\mathcal{A}_{\mathrm{orb}}^B(\Delta t)\,\mathbf{B}^n
\]
that is divergence-preserving by construction [1207.2955].

In $\mathbf{H}(\mathrm{curl})$ finite elements for the magnetic advection-diffusion problem, the elementwise Lie derivative
\[
L_{\boldsymbol{\beta}}\mathbf{u}
=
-\boldsymbol{\beta}\times(\nabla\times\mathbf{u})
+
\nabla(\boldsymbol{\beta}\cdot\mathbf{u})
\]
is modified by a lifting of interelement jumps. The discrete magnetic advection operator is
\[
\tilde L_{\boldsymbol{\beta},h}\mathbf{u}_h|_T
=
L_{\boldsymbol{\beta},h}\mathbf{u}_h|_T
-
\mathbf{r}_\alpha|_T(\boldsymbol{\beta}\cdot\mathbf{n}\,\llbracket\mathbf{u}_h\rrbracket),
\]
where the lifting operator is defined locally by
\[
\int_T \mathbf{r}_\alpha(\mathbf{v})\cdot\mathbf{w}_h\,dx
=
\int_{\partial T}\alpha|_{\partial T}\,\mathbf{v}\cdot\mathbf{w}_h\,ds.
\]
This turns facet fluxes into a volume operator and makes the advection term usable inside SUPG residual stabilization [2509.09913].

## 3. Geometric and structure-preserving realizations

A central division in the literature is between operators that approximate a transport PDE explicitly and operators in which magnetic advection is built into the geometric representation itself.

In the DEC variational integrator based on Newcomb’s Lagrangian, the magnetic field is represented as a discrete 2-form on mesh faces, with flux dofs $\langle B,\sigma^2\rangle$. Magnetic advection is not implemented through a stencil for $\nabla\times(\mathbf{v}\times\mathbf{B})$; rather, it is hard-wired through the algebraic identity
\[
\langle B,\sigma^2\rangle=\langle B_0,\sigma_0^2\rangle.
\]
The material faces move with the mesh, while the face flux coefficients remain constant. In that sense, the discrete magnetic advection operator is the identity on flux cochains combined with time-dependent mesh geometry. If $\mathrm{d}B_0=0$ initially, then $\mathrm{d}B=0$ for all time because the cochain values and incidence structure remain unchanged [1408.1346].

A related ideal-MHD variational discretization in 2D represents both velocity and magnetic field as primal 1-forms on a staggered Cartesian grid. The discrete induction equation is
\[
(\Delta_t B^x)_{i,j+1/2,n+1/2}+\phi^x_{i,j+1/2,n+1/2}(V,B)=0,
\]
\[
(\Delta_t B^y)_{i+1/2,j,n+1/2}+\phi^y_{i+1/2,j,n+1/2}(V,B)=0,
\]
where the discrete magnetic advection operator is the pair $(\phi^x,\phi^y)$ built from centered differences of the cell-centered flux
\[
F_{i,j,n+1/2}
=
\bar V^x_{i,j,n+1/2}\,\bar B^y_{i,j,n+1/2}
-
\bar V^y_{i,j,n+1/2}\,\bar B^x_{i,j,n+1/2}.
\]
This is the discrete analogue of Lie advection of the magnetic flux 2-form, expressed in edge-based components [1707.03227].

In vector-potential orbital advection, the large mean orbital flow is split as
\[
\mathbf{u}=\bar u_\phi(r,z)\,\hat{\boldsymbol{\phi}}+\mathbf{u}',
\]
and the operator
\[
\mathcal{L}=\frac{\bar u_\phi(r,z)}{r}\,\frac{\partial}{\partial\phi}
\]
is treated implicitly through a Fourier shift. The discrete magnetic advection operator then consists of two pieces: explicit high-order finite-difference discretization of
\[
-u'_j\partial_j A_i + u_j\partial_i A_j
\]
and exact orbital translation by $\mathcal{Q}(t)=e^{-t\mathcal{L}}$ applied at every RK stage. This suggests a split operator in which the stiff, nearly uniform magnetic advection is handled as translation rather than by local differencing [1708.05057].

These constructions show that “discrete magnetic advection operator” may refer either to an explicit discrete differential operator or to a transport map induced by mesh motion, edge EMFs, or spectral translation. This suggests a broad operator class rather than a single algebraic template.

## 4. Stability, conservation, and divergence control

Across formulations, the decisive design criterion is usually not pointwise monotonicity but preservation of the invariants associated with magnetic transport.

In SBP-SAT finite differences, stability is derived in the SBP norm
\[
\|V^n\|_P^2=(V^n,V^n)_P.
\]
If the SAT penalty parameters satisfy the inflow-sign conditions
\[
(\sigma_R)_{N-1,j}\le \tfrac12 u^{1,-}(1,y_j),\quad
(\sigma_L)_{0,j}\le -\tfrac12 u^{1,+}(0,y_j),
\]
\[
(\sigma_U)_{i,M-1}\le \tfrac12 u^{2,-}(x_i,1),\quad
(\sigma_D)_{i,0}\le -\tfrac12 u^{2,+}(x_i,0),
\]
then the fully discrete scheme satisfies
\[
\|V^n\|_P^2\le e^{KT}\|V^0\|_P^2.
\]
The essential point is that the central discrete advection operator becomes energy-stable only after boundary fluxes are controlled by SBP-SAT [1102.0494].

In the mimetic staggered-grid construction, stability is stronger: when $\alpha+\beta=1$,
\[
\mathbf{p}^T X_p C_{pq}\mathbf{q}+\mathbf{q}^T X_q C_{qp}\mathbf{p}=0,
\]
so the semi-discrete wave system preserves
\[
\|\mathbf{x}\|^2=\|\mathbf{p}\|^2+\|\mathbf{q}\|^2.
\]
If instead $\alpha=\beta=1$ or $\alpha=\beta=0$, the discrete operator is not skew-adjoint, the spectrum acquires eigenvalues with positive real parts, and the solution grows exponentially in time [2412.20106]. A common misconception is that centered transport operators are automatically neutral; in these plasma transport discretizations, neutral stability depends on an exact discrete adjoint pairing, not merely on centered differencing.

In CT orbital advection, divergence control is topological. The face-centered updates are written as edge-EMF circulations, so the discrete divergence change is identically zero: each edge EMF appears twice with opposite sign when the divergence of the updated field is taken. The method therefore preserves the discrete divergence-free condition of the magnetic field to machine precision during the linear transport step [1207.2955].

In the variational ideal-MHD schemes, conservation is more restrictive still. Because the advection equations are built into the discrete action, magnetic flux through each material face is preserved exactly, and there is no separate induction solve that could introduce numerical reconnection. One formulation states that the method is “free of numerical resistivity” and that magnetic field line topology is preserved, with total energy, magnetic helicity, and cross helicity satisfied within machine accuracy in 2D numerical examples [1707.03227]. Another reports that numerical reconnection does not take place when singular current sheets are present [1408.1346]. A plausible implication is that these methods realize magnetic advection as exact flux freezing at the discrete level rather than as approximate Eulerian transport.

For the $\mathbf{H}(\mathrm{curl})$ SUPG method, coercivity is obtained in the norm
\[
\|\mathbf{u}_h\|^2
=
\varepsilon\|\nabla\times\mathbf{u}_h\|_0^2
+
\rho_0\|\mathbf{u}_h\|_0^2
+
\sum_T\delta_T\|\tilde L_{\boldsymbol{\beta},h}\mathbf{u}_h\|_{0,T}^2
+\text{jump and boundary terms},
\]
and one proves
\[
a_h(\mathbf{v}_h,\mathbf{v}_h)\ge \frac12 \|\mathbf{v}_h\|^2
\]
under the upwind condition on $\alpha$ and the bound
\[
\delta_T\le
\min\left\{
\frac{h_T^2}{2C_{inv}^2\varepsilon},
\frac{\rho_0}{2\|\gamma\|_{L^\infty(T)}^2}
\right\}.
\]
Here the discrete magnetic advection operator itself enters the norm, so stability is directly tied to controlling $\tilde L_{\boldsymbol{\beta},h}\mathbf{u}_h$ [2509.09913].

## 5. Accuracy, spectral behavior, and dissipation

The accuracy of a discrete magnetic advection operator depends strongly on how transport is split, where unknowns are placed, and whether stabilization acts in the same direction as the physical transport.

For the SBP-SAT induction scheme, the advection operator inherits the underlying derivative orders. SBP2 gives second-order interior accuracy and overall second-order behavior; SBP4 gives fourth-order interior and second-order boundary accuracy, with overall third-order spatial behavior. In the rotating-pulse test at $t=2\pi$, the reported discrete $l^2$ errors were
\[
6.9\times10^{-1},\ 2.1\times10^{-1},\ 5.5\times10^{-2},\ 1.3\times10^{-2},\ 3.3\times10^{-3}
\]
for SBP2 on $40\times40$ through $640\times640$ grids, with rates approaching $2.0$, and
\[
8.0\times10^{0},\ 5.0\times10^{-1},\ 4.5\times10^{-2},\ 5.1\times10^{-3},\ 6.4\times10^{-4}
\]
for SBP4 with rates approaching $3.0$ [1102.0494].

In orbital-advection MHD, the dominant accuracy issue is not solely formal order but dissipation under large background flow. In the PLUTO FARGO-MHD scheme, the orbital step is a conservative shift-plus-interpolation operator; with PPM reconstruction it is third-order accurate in that step. Reported MRI growth rates in a global disk show that at $16$ cells per wavelength, FARGO-MHD with third-order orbital advection reaches $\gamma^{MRI}\approx0.74$, compared with $0.70$ for the standard scheme, close to the ideal value $0.75$ [1207.2955]. In Pencil Code orbital advection, the combined RK3-2N plus spectral-shift scheme is proven third order in time, and a 1D time-accuracy test gives a slope of $3$ in the log-log error plot [1708.05057].

For field-aligned transport, spectral fidelity becomes the decisive metric. In the anti-symmetry representation of parallel diffusion, the 1D anti-symmetric discrete response is
\[
\mathcal{R}_{\mathrm{AS}}(k)
=
\frac{2}{\Delta x^2}
\left[
2\left(\frac{k\Delta x}{2}\right)\sin\left(\frac{k\Delta x}{2}\right)
\right],
\]
whereas the conventional finite-difference response is
\[
\mathcal{R}_{\mathrm{FD}}(k)
=
\frac{2}{\Delta x^2}\bigl[1-\cos(k\Delta x)\bigr].
\]
The paper states that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart, and that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods [2412.01927]. Although that work concerns parallel diffusion rather than magnetic induction directly, the operator structure is explicitly advection-like, with diffusion recast as a flow operator aligned with $\mathbf{b}$. A plausible implication is that the same skew-symmetric staggered structure is attractive for discrete magnetic advection along field lines.

In the $\mathbf{H}(\mathrm{curl})$ SUPG method, the stabilization direction is the discrete magnetic advection operator itself. The analysis gives energy-norm convergence of order $k+1/2$ and $L^2$-norm convergence of order $k+1$ for smooth advection-dominated cases, while numerical experiments indicate that using only the jump term $S_h^1$ improves stability but does not recover full accuracy, whereas the SUPG residual term built on $\tilde L_{\boldsymbol{\beta},h}$ does [2509.09913].

## 6. Applications, distinctions, and recurrent misunderstandings

The application range of discrete magnetic advection operators is unusually broad. Zita’s diffusivity-gradient analysis was motivated by solar dynamo studies, and the paper discusses steep turbulent magnetic diffusivity gradients near the tachocline and in the upper convection zone, as well as space and laboratory plasmas including the Reversed-Field Pinch [1309.4398]. SBP-SAT induction operators were developed for prescribed-velocity magnetic induction problems and tested on rotating divergence-free fields [1102.0494]. Mimetic staggered-grid field-line transport was applied to electrostatic shear Alfvén waves in tokamak-like geometry [2412.20106]. Orbital-advection operators were designed for accretion and proto-planetary disks in Cartesian, cylindrical, and spherical coordinates [1207.2955; 1708.05057]. Variational flux-freezing schemes were used to study singular current sheets and ideal coalescence instability without numerical reconnection [1408.1346; 1707.03227]. The $\mathbf{H}(\mathrm{curl})$ SUPG formulation targets the magnetic advection-diffusion problem for vector potentials and reduced MHD models [2509.09913].

Several distinctions recur across this literature.

First, “magnetic advection” does not always mean advection of $\mathbf{B}$ by a prescribed fluid velocity. In some papers it means advection of the magnetic field itself by $\mathbf{u}$; in others it means advection along magnetic field lines, $\mathbf{b}\cdot\nabla$; in the variational DEC setting it means exact transport of magnetic flux cochains; and in Zita’s derivation it means an effective transport velocity proportional to $\nabla\sqrt{\eta}$ [1309.4398; 2412.20106; 1408.1346].

Second, advection by diffusivity gradients should not be conflated with mean-field diamagnetic pumping. Zita explicitly contrasts the MHD result
\[
\mathbf{U}_\eta=\nabla\sqrt{\eta}
\]
with the mean-field weak-field diamagnetic pumping velocity
\[
\mathbf{U}^{\mathrm{dia}}=-\frac12\nabla\eta_T,
\]
emphasizing that the derivations, sign, and dependence differ [1309.4398]. This is a recurrent source of confusion because both appear as transport terms driven by diffusivity variation.

Third, centered schemes are not necessarily unstable, and upwinding is not the only route to robust magnetic transport. SBP-SAT uses central high-order interior stencils with weak boundary penalties and still proves stability [1102.0494]. The mimetic plasma operator is energy-preserving only when the advective and divergence forms are combined so that $\alpha+\beta=1$ [2412.20106]. The $\mathbf{H}(\mathrm{curl})$ SUPG method uses a lifted Lie derivative plus residual stabilization rather than purely one-sided stencils [2509.09913].

Fourth, the advective gauge in vector-potential orbital advection is analytically natural but numerically problematic. The Pencil Code work states that the advective gauge is unstable in practice, and argues that implicit treatment of the orbital advection term makes the pseudo-advective gauge sufficient because it precludes the shear term from canceling the advection term [1708.05057].

Taken together, these results define the discrete magnetic advection operator not as a single formula but as a research domain organized around a common problem: how to represent magnetic transport discretely while retaining, as appropriate to the formulation, divergence preservation, flux freezing, skew-adjointness, exact or weak conservation laws, controlled interface fluxes, and high spectral fidelity.

Source: https://www.emergentmind.com/topics/discrete-magnetic-advection-operator