Numerical Simulations of Dynamos Generated in Spherical Couette Flows

We numerically investigate the efficiency of a spherical Couette flow at generating a self-sustained magnetic field. No dynamo action occurs for axisymmetric flow while we always found a dynamo when non-axisymmetric hydrodynamical instabilities are excited. Without rotation of the outer sphere, typical critical magnetic Reynolds numbers $Rm_c$ are of the order of a few thousands. They increase as the mechanical forcing imposed by the inner core on the flow increases (Reynolds number $Re$). Namely, no dynamo is found if the magnetic Prandtl number $Pm=Rm/Re$ is less than a critical value $Pm_c\sim 1$. Oscillating quadrupolar dynamos are present in the vicinity of the dynamo onset. Saturated magnetic fields obtained in supercritical regimes (either $Re>2 Re_c$ or $Pm>2Pm_c$) correspond to the equipartition between magnetic and kinetic energies. A global rotation of the system (Ekman numbers $E=10^{-3}, 10^{-4}$) yields to a slight decrease (factor 2) of the critical magnetic Prandtl number, but we find a peculiar regime where dynamo action may be obtained for relatively low magnetic Reynolds numbers ($Rm_c\sim 300$). In this dynamical regime (Rossby number $Ro\sim -1$, spheres in opposite direction) at a moderate Ekman number ($E=10^{-3}$), a enhanced shear layer around the inner core might explain the decrease of the dynamo threshold. For lower $E$ ($E=10^{-4}$) this internal shear layer becomes unstable, leading to small scales fluctuations, and the favorable dynamo regime is lost. We also model the effect of ferromagnetic boundary conditions. Their presence have only a small impact on the dynamo onset but clearly enhance the saturated magnetic field in the ferromagnetic parts. Implications for experimental studies are discussed.

💡 Research Summary

The paper presents a comprehensive numerical investigation of dynamo action in spherical Couette flow, i.e., the flow generated between a rotating inner sphere and a surrounding outer sphere that may be stationary or rotating. The authors explore a wide parameter space defined by the Reynolds number (Re), magnetic Reynolds number (Rm), magnetic Prandtl number (Pm = Rm/Re), Ekman number (E = ν/(ΩL²)), and Rossby number (Ro = ΔΩ/Ω). Their main findings can be summarized as follows.

1. Axisymmetric versus non‑axisymmetric flow – When the flow remains axisymmetric (inner sphere rotating, outer sphere fixed), no dynamo is observed regardless of how large Re becomes. The magnetic field decays because the induced currents are confined to toroidal loops that do not generate a net poloidal component. By contrast, once the Reynolds number exceeds a hydrodynamic critical value (Re > Re_c) and non‑axisymmetric instabilities (e.g., wavy vortices, shear‑layer instabilities) develop, the flow acquires sufficient three‑dimensional structure to sustain a magnetic field. In this regime the critical magnetic Reynolds number is of order a few thousand (Rm_c ≈ 2–5 × 10³).

2. Role of magnetic Prandtl number – Because Rm = Pm · Re, the existence of a dynamo depends critically on the ratio Pm. The simulations show a sharp threshold Pm_c ≈ 1: for Pm < Pm_c the magnetic field always decays, even when Re is very large. When Pm ≥ Pm_c, the dynamo can be excited and, once supercritical (either Re > 2 Re_c or Pm > 2 Pm_c), the system settles into a saturated state where magnetic and kinetic energies are in approximate equipartition.

3. Effect of global rotation – Introducing rotation of the outer sphere adds Coriolis forces, characterized by the Ekman number. Two values are examined, E = 10⁻³ and 10⁻⁴. Global rotation modestly lowers the critical Pm by roughly a factor of two, but a particularly favorable regime appears when the inner and outer spheres rotate in opposite directions (Rossby number Ro ≈ ‑1). At E = 10⁻³ this configuration creates a strong shear layer around the inner sphere, dramatically reducing the dynamo threshold to Rm_c ≈ 300. The enhanced shear concentrates electric currents and amplifies the inductive coupling between toroidal and poloidal fields. However, when the Ekman number is reduced further to E = 10⁻⁴ the shear layer itself becomes unstable, generating small‑scale turbulence that destroys the coherent current sheet, and the low‑Rm dynamo regime disappears.

4. Boundary magnetic properties – The authors also model ferromagnetic boundary conditions (high relative permeability μ_r ≫ 1) on the inner and/or outer sphere. These boundaries have only a minor influence on the onset of dynamo action (the critical Rm changes by less than 10 %). Nevertheless, in the saturated regime the magnetic energy stored in the ferromagnetic parts is significantly larger—up to about 20 % higher—because the high‑μ material concentrates magnetic flux.

5. Temporal behavior near onset – Close to the dynamo threshold the magnetic field often takes the form of an oscillating quadrupole. This time‑dependent solution reflects a nonlinear interaction between the dominant non‑axisymmetric flow mode and the magnetic field, leading to a periodic exchange of energy between toroidal and poloidal components.

6. Implications for laboratory experiments – The study provides clear guidance for experimental dynamo setups. To achieve self‑sustained magnetic fields in a spherical Couette device, one must (i) drive the flow into a non‑axisymmetric regime by exceeding the hydrodynamic instability threshold, (ii) operate with a fluid whose magnetic Prandtl number is of order unity (e.g., liquid sodium at temperatures that give Pm ≈ 10⁻⁵ is far below the threshold, so a higher‑Pm alloy or a different working fluid would be required), (iii) consider rotating the outer sphere in the opposite direction to the inner sphere to exploit the low‑Rm shear‑layer dynamo, and (iv) recognize that ferromagnetic liners can boost the saturated field but will not substantially lower the critical rotation rates.

In summary, the paper demonstrates that dynamo action in spherical Couette flow is fundamentally linked to the presence of non‑axisymmetric hydrodynamic instabilities, a magnetic Prandtl number above a critical value of order one, and, when global rotation is present, to the formation of a strong shear layer at moderate Ekman numbers. The findings bridge the gap between idealized theoretical models and realistic laboratory experiments, offering a roadmap for achieving dynamo action in laboratory‑scale spherical Couette devices.