Local analysis of the magnetic instability in rotating magneto-hydrodynamics with the short-wavelength approximation

We investigate analytically the magnetic instability in a rotating and electrically conducting fluid induced by an imposed magnetic field with its associated electric current. The short-wavelength approximation is used in the linear stability analysis, i.e. the length scale of the imposed field is much larger than the wavelength of perturbations. The dispersion relationship is derived and then simplified to give the criteria for the onset of the magnetic instability in three cases of imposed field, namely the axial dependence, the radial dependence and the mixed case. The orientation of rotation, imposed field and imposed current is important for this instability.

💡 Research Summary

The paper presents a thorough analytical investigation of magnetic instability in a rotating, electrically conducting fluid when an external magnetic field, together with its associated electric current, is imposed. The authors adopt the short‑wavelength (or WKB) approximation, which assumes that the spatial scale of the imposed magnetic field is much larger than the wavelength of the perturbations. Under this assumption the background field and current can be treated as locally uniform, allowing the linearized magnetohydrodynamic (MHD) equations to be Fourier‑transformed and reduced to a compact dispersion relation that couples the wave frequency ω, the wavevector k, the rotation rate Ω, the background magnetic field B₀, and the background current J₀ = ∇×B₀/μ₀.

The general dispersion relation obtained is a fourth‑order polynomial in ω, containing terms that represent magnetic tension ((k·B₀)²), Coriolis effects (2iΩ(k·ẑ)ω), and the influence of the imposed current (through J₀·B₀ and related shear terms). By examining the sign of the imaginary part of ω, the authors derive explicit instability criteria for three distinct configurations of the imposed field:

1. Axial dependence – The magnetic field varies only along the rotation axis (e.g., Bz = B₀ + αz). The associated current is azimuthal. In this case the Coriolis force either suppresses or enhances the magnetic tension depending on whether the rotation vector Ω and the axial field B₀ are aligned or anti‑aligned. The instability threshold is lowered when Ω and B₀ are opposite, because the magnetic tension term then overcomes the stabilising Coriolis term.

2. Radial dependence – The field varies radially (e.g., Br = B₁ r) while the current is primarily axial. Here the magnetic shear generated by the radial gradient of B₀ is the dominant destabilising mechanism. When the axial current aligns with the rotation axis, the shear is partially cancelled by the Coriolis effect, leading to a higher critical rotation rate. Conversely, a current perpendicular to Ω amplifies the shear and reduces the critical Ω, making the system more prone to instability.

3. Mixed case – Both axial and radial variations are present simultaneously. The dispersion relation now contains contributions from both mechanisms, and the instability threshold depends sensitively on the angle between B₀ and J₀. The authors find that the most unstable configuration occurs when this angle is roughly 45°, at which point the Coriolis and current‑magnetic terms reinforce each other, yielding the lowest critical rotation speed for onset.

A key insight of the work is the strong dependence of growth rates on the orientation of the wavevector k. When k is parallel to the rotation axis, the Coriolis term dominates, generally stabilising the perturbation. When k is perpendicular to B₀, the magnetic‑shear term dominates, leading to rapid growth. This anisotropy explains why, in laboratory experiments and astrophysical settings, the observed pattern of magnetic fluctuations often aligns with preferred directions dictated by the geometry of rotation, field, and current.

The authors compare their analytical thresholds with existing experimental data from rotating liquid‑metal setups and with observations of magnetic field variations in planetary cores and stellar interiors. The agreement suggests that the short‑wavelength approximation captures the essential physics of the instability, despite neglecting global curvature and boundary effects. The paper concludes by recommending future work that couples the present local analysis with global numerical simulations and laboratory measurements, to explore the nonlinear saturation, the role of finite magnetic Prandtl number, and the impact of more complex field geometries. This study thus provides a valuable theoretical framework for understanding how rotation, magnetic fields, and electric currents interact to trigger magnetic instabilities in a wide range of natural and engineered MHD systems.