1 Overview and Definition

1.1 Dimensionless form and standard expressions

The Lundquist number \(S\) is a dimensionless parameter used in magnetohydrodynamics (MHD) to compare how quickly magnetic structures diffuse relative to how quickly Alfvénic disturbances propagate. A common scaling form is \[ S \sim \frac{\mu_0 L V_A}{\eta}, \] where \(L\) is a characteristic length scale, \(V_A\) is the Alfvén speed, and \(\eta\) is the magnetic diffusivity. In many texts, \(\eta\) is linked to resistivity through \(\eta \propto 1/\sigma\), with \(\sigma\) the electrical conductivity, so that higher conductivity corresponds to smaller \(\eta\) and therefore larger \(S\).

1.2 Physical interpretation as a timescale ratio

Physically, \(S\) compares two characteristic timescales. One is an Alfvénic timescale, often written as \(\tau_A \sim L/V_A\), describing the time for information to travel along magnetic-field-aligned paths at roughly the Alfvén speed. The other is a diffusion timescale for magnetic fields, \(\tau_\eta \sim L^2/\eta\), describing how resistivity smooths magnetic gradients over the system size. With these definitions, \[ S \sim \frac{\tau_\eta}{\tau_A}, \] so large \(S\) indicates diffusion is slow compared with wave propagation, while small \(S\) indicates resistive effects act on comparable or faster timescales.

1.3 Relation to magnetic diffusivity and resistivity

Magnetic diffusivity \(\eta\) encodes how resistive effects enter the induction equation. In resistive MHD, the magnetic field evolves according to a balance between advection by bulk plasma motion and diffusion driven by finite resistivity. When \(\eta\) is small, magnetic flux is approximately “frozen” into the plasma over dynamical times, supporting an ideal-MHD description. When \(\eta\) is larger, diffusion competes more effectively with advection, breaking strict flux-freezing behavior.

1.4 Typical magnitudes and qualitative regimes

In many laboratory and computational settings, the Lundquist number can vary widely depending on geometry, plasma parameters, and effective transport coefficients. Qualitatively:

  • Large \(S\) (slow diffusion): ideal-MHD behavior is more accurate for a broader range of timescales, and current layers can become thinner before resistive effects dominate.
  • Small \(S\) (fast diffusion): resistive smoothing occurs more readily, weakening the ideal-MHD approximation and promoting earlier influence of diffusion on magnetic-field evolution.

These regimes are not absolute; they depend on the particular problem setup and on which length scale actually governs diffusion in the dynamics under study.

2 Derivation and Context in MHD

2.1 Governing equations and the role of resistivity

2.1.1 Induction equation and diffusive term

The resistive MHD induction equation can be written schematically as \[ \frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{u}\times \mathbf{B}) + \eta \nabla^2 \mathbf{B}, \] where \(\mathbf{B}\) is the magnetic field and \(\mathbf{u}\) is the bulk plasma velocity. The first term on the right represents field advection and stretching by flow, while the second term represents diffusion of magnetic gradients. The competition between these terms under nondimensionalization is what generates the Lundquist number.

2.1.2 Alfvén timescale identification

The emergence of \(V_A\) follows from the momentum equation and magnetic forces. Using a characteristic mass density \(\rho\), the Alfvén speed is \[ V_A \sim \frac{B}{\sqrt{\mu_0 \rho}}, \] so that the timescale for Alfvénic response over distance \(L\) is \(\tau_A \sim L/V_A\). In typical dynamical problems, this is the relevant propagation time against which diffusion is compared.

2.2 Nondimensionalization of MHD

2.2.1 Choosing characteristic length and velocity scales

Nondimensionalization begins by introducing dimensionless variables such as \(\tilde{\mathbf{x}}=\mathbf{x}/L\) and \(\tilde{t}=t/\tau\), along with a characteristic velocity \(\,U\) that can be taken as \(V_A\) in magnetic-dynamical regimes. Inserting these scalings into the induction equation makes the relative size of the diffusive term explicit. The parameter multiplying diffusion is proportional to \(\eta/(L U)\) and thus becomes the inverse of a diffusion-to-advection ratio.

2.2.2 Recovering \(S\) in reduced models

When one chooses \(U\sim V_A\) and \(\tau\sim \tau_A\), the dimensionless induction equation contains a factor that can be identified as \(1/S\). Since \(\tau_A \sim L/V_A\) and \(\tau_\eta \sim L^2/\eta\), the ratio \(\tau_\eta/\tau_A\) yields the standard scaling \[ S \sim \frac{L V_A}{\eta}, \] up to unit conventions. Including \(\mu_0\) in the definition depends on how \(\eta\) is expressed (e.g., via resistivity versus conductivity), but the physical meaning remains the same: the Lundquist number is the diffusion-versus-Alfvén competition.

3 Connection to Other Dimensionless Numbers

3.1 Comparison with magnetic Reynolds number

The magnetic Reynolds number \(R_m\) is a dimensionless quantity that also compares advection of magnetic fields to resistive diffusion. In common formulations, \[ R_m \sim \frac{U L}{\eta}, \] which, with \(U\sim V_A\), is essentially the same scaling as \(S\). Differences arise from how each parameter is defined with regard to geometry, characteristic velocity, or the chosen form of \(\eta\). Thus, \(S\) can be viewed as the magnetic Reynolds number specialized to MHD contexts where the Alfvén speed sets the dynamical velocity scale.

3.2 Interplay with Lundquist–Reynolds–Prandtl-type parameters

MHD problems often introduce additional dimensionless groups, including the (fluid) Reynolds number and the Prandtl number (and their magnetic analogs). These control how momentum transport (viscous effects) and thermal transport (conductive effects) interplay with magnetic diffusion. In resistive MHD, the Lundquist number sets the magnetic diffusion strength, while other parameters determine whether velocity gradients, temperature gradients, or both are sharply maintained. Consequently, the observed dynamics can depend on multiple competing transport processes rather than on \(S\) alone.

Because \(\eta\) is tied to conductivity \(\sigma\), and conductivity depends on plasma temperature and other microscopic effects, \(S\) inherits parameter dependence from those underlying scalings. In broad terms, higher conductivity reduces \(\eta\), increasing \(S\), while stronger resistive transport increases \(\eta\) and reduces \(S\). In practice, estimating \(S\) requires careful attention to what effective diffusivity is relevant: microscopic resistivity, collisionally determined transport, or enhanced (turbulent or numerical) diffusion.

4 Implications for Magnetic Field Dynamics

4.1 Ideal-MHD limit vs resistive effects

The ideal-MHD approximation corresponds to neglecting resistive diffusion over the timescale of interest. Since diffusion enters with a factor effectively proportional to \(1/S\), large Lundquist number supports a regime in which magnetic fields evolve primarily through advection and stretching by the flow. However, ideal behavior can break down locally when magnetic structures become sufficiently thin, so the relevant length scale for diffusion can be much smaller than the system size \(L\). The implications are therefore spatial as well as temporal.

4.2 Current sheets and thinning processes

Many MHD scenarios produce current layers where magnetic gradients concentrate. If the system allows such layers to thin, the local diffusion time over the current-sheet thickness \(\delta\) becomes \(\tau_\eta(\delta)\sim \delta^2/\eta\). Even when global \(S\) is large, diffusion can become significant once \(\delta\) is small enough that \(\tau_\eta(\delta)\) approaches the local dynamical time. This is a key reason \(S\) is often used to predict when and where resistive effects will matter.

4.3 Stability and onset of resistive phenomena

4.3.1 Criteria often expressed using \(S\)

Stability boundaries and the onset of resistive behaviors are frequently expressed in terms of \(S\) because it controls the effective strength of the diffusive term in the governing equations. In many model systems, thresholds scale as power laws in \(S\), with higher \(S\) shifting certain transitions to smaller spatial scales or to different growth-rate regimes. The specific exponents and conditions vary by instability type, equilibrium profile, and boundary conditions, but the general logic remains: increasing \(S\) weakens diffusion and can delay or reshape resistive dynamics.

5 Lundquist Number in Magnetic Reconnection Studies

5.1 Resistive reconnection basics

Magnetic reconnection is a process in which magnetic field lines change connectivity, enabled in resistive MHD by non-ideal behavior within localized regions where the induction equation cannot be treated as purely ideal. In a resistive framework, reconnection requires breaking the ideal constraint through effects controlled by \(\eta\). The Lundquist number organizes how rapidly this non-ideal region forms and how efficiently it allows topological change.

5.2 Role of \(S\) in reconnection rate scaling

In reconnection studies, the dependence of reconnection rate on \(S\) is central because it indicates whether reconnection becomes slower as diffusion weakens or whether additional physics (e.g., geometry-dependent thinning, kinetic corrections, or turbulent transport) can modify the scaling. In simple resistive models, reconnection often shows a trend in which increasing \(S\) reduces the effective diffusion influence and can lower the rate. Nonetheless, reconnection dynamics can be regime-dependent, and “effective” scaling may incorporate assumptions about the current-sheet structure and how it evolves with system parameters.

5.3 Boundary between fast and slow reconnection regimes

A commonly discussed conceptual distinction is between reconnection regimes that are comparatively fast (insensitive to very small resistivity, producing substantial rates even at large \(S\)) and slow (rates that diminish markedly as \(S\) grows). While the exact boundary is model- and geometry-dependent, \(S\) remains a key control parameter that indicates how strong resistive effects are relative to ideal dynamics. Observed or simulated reconnection behavior near the boundary often reflects how current sheets thin and how non-ideal regions respond to increasing \(S\).

6 Scaling Laws and Model Dependence

6.1 Dependence on geometry and boundary conditions

6.1.1 System size vs local diffusion length scales

Scaling arguments frequently use a global system size \(L\) to define \(S\), but many processes are controlled by local length scales such as the current-sheet half-thickness or diffusion-region size. If these local scales scale with \(S\) itself, then the effective importance of resistivity can differ from naive expectations based solely on \(L\). As a result, scaling laws for timescales or rates often involve combinations of \(S\) with empirically or theoretically determined relations for the local thickness \(\delta(S)\).

6.2 Self-similar scalings in simplified settings

In idealized geometries, such as planar or axisymmetric setups with simplified flow patterns, one can derive self-similar behaviors for how magnetic gradients evolve. Self-similar analyses often express the growth of thinning or the evolution of diffusion-region thickness as power laws in time, which then translate into scaling laws involving \(S\). These simplified results help clarify which assumptions are responsible for particular exponents and how robust a predicted dependence on \(S\) might be.

6.3 Numerical and theoretical considerations

Because \(S\) can be very large in real systems, numerical simulations may rely on lower effective \(S\) than those found experimentally or astrophysically. This raises questions of convergence and whether asymptotic scaling at extremely large \(S\) is already approached in accessible computations. Additionally, discretization, choice of resistivity model, and numerical diffusion can all affect the effective Lundquist number and thus the observed dynamics. The interpretation of results therefore requires attention to how \(S\) is implemented and measured in the simulation.

7 Measurement and Estimation in Practice

7.1 Determining characteristic scales for \(L\) and \(V_A\)

Estimating \(S\) requires selecting appropriate characteristic values:

  • \(L\) may be taken as a system size, an equilibrium gradient scale, or a more local length relevant to the process under consideration.
  • \(V_A\) depends on magnetic field strength and plasma density, typically using representative values such as peak magnetic field or density in the region of interest.

Different reasonable choices of \(L\) and \(V_A\) can lead to different numerical values for \(S\), so the definition should align with the physical question being asked.

7.2 Estimating \(\eta\) (effective vs microscopic)

The parameter \(\eta\) can correspond to microscopic resistivity determined by collisions, but in many practical contexts an effective diffusivity is more relevant. Effective diffusion may include contributions from turbulence, anomalous transport, or parameterizations used in models. Since \(S\propto 1/\eta\), uncertainties in \(\eta\) can strongly propagate into uncertainty in \(S\), affecting comparisons between theory, simulation, and experiment.

7.3 Uncertainty sources and sensitivity

Uncertainty in \(S\) typically comes from measurement errors in density and magnetic field (affecting \(V_A\)), ambiguity in the appropriate \(L\), and model-form uncertainty for \(\eta\). Sensitivity studies often show that the largest contributions to uncertainty stem from \(\eta\) or from which spatial scale is selected for diffusion. Consequently, reporting a Lundquist number alongside its underlying definitions (what \(L\) and \(\eta\) represent) is essential for meaningful interpretation.

8 Applications and Use Cases

8.1 Laboratory plasmas (e.g., confinement experiments)

In controlled fusion and related plasma experiments, Lundquist number is used to gauge how resistive effects compare with ideal MHD dynamics. It informs expectations about current-layer formation, stability behavior, and how quickly resistive processes may develop. Because experimental plasmas can be complex, researchers often compute \(S\) using equilibrium reconstruction and then interpret deviations from ideal predictions in light of finite resistivity.

8.2 Space and astrophysical plasmas (general resistive MHD context)

In space physics, reconnection and other resistive processes are frequently analyzed using resistive MHD concepts, with \(S\) serving as an indicator of the relative strength of diffusion. Even when fully kinetic effects may be important at small scales, resistive MHD can still provide a useful macroscopic baseline. In such applications, \(S\) helps translate measured or inferred plasma conditions into expectations about where resistive dynamics may become significant.

8.3 Engineering and technological contexts where applicable (e.g., liquid-metal MHD)

Beyond plasmas in traditional settings, MHD ideas appear in engineering contexts such as flows of conductive fluids and liquid-metal systems. While the microscopic basis differs from hot plasma resistivity, the mathematical structure of diffusion-versus-advection competition persists. In these contexts, a Lundquist-type number can help characterize whether magnetic fields remain closely coupled to the moving conductor or whether resistive diffusion significantly alters the field distribution.

9 Limitations and Extensions

9.1 Assumptions behind resistive MHD and constant \(\eta\)

Standard definitions of \(S\) usually assume a resistive MHD model with a diffusion coefficient \(\eta\) that is either constant or slowly varying relative to the dominant gradients. Real systems may exhibit strong spatial and temporal variations in conductivity or diffusivity. When \(\eta\) varies, a single global Lundquist number may be insufficient, and a more local or profile-aware definition becomes necessary.

9.2 When Hall/kinetic effects matter (conceptual extension)

At sufficiently small length scales or high frequencies, effects beyond resistive MHD—such as Hall physics or kinetic corrections—can become relevant. In that case, the Lundquist number remains a useful indicator of resistive coupling at a macroscopic level, but it no longer uniquely controls the dynamics. Additional scale lengths and non-ideal terms can set the structure and speed of reconnection or wave behavior, so interpreting results purely in terms of \(S\) can be misleading.

9.3 Interpretation cautions for inhomogeneous plasmas

In inhomogeneous plasmas, the characteristic \(L\) and effective \(V_A\) may change across the domain, and current sheets can form preferentially in certain regions. A single value of \(S\) may therefore mask the fact that local conditions correspond to different effective diffusion strengths. Practical use of the Lundquist number in such contexts often involves specifying where it is evaluated and whether it is intended to describe global dynamics, local current-sheet evolution, or a transition between behaviors.