A physical conjecture for the dipolar–multipolar dynamo transition

Abstract

In numerical simulations of planetary dynamos there is an abrupt transition in the dynamics of both the velocity and magnetic fields at a ‘local’ Rossby number of 0.1. For smaller Rossby numbers there are helical columnar structures aligned with the rotation axis, which efficiently maintain a dipolar field. However, when the thermal forcing is increased, these columns break down and the field becomes multi-polar. Similarly, in rotating turbulence experiments and simulations there is a sharp transition at a Rossby number of ${\sim}0.4$. Again, helical axial columnar structures are found for lower Rossby numbers, and there is strong evidence that these columns are created by inertial waves, at least on short time scales. We perform direct numerical simulations of the flow induced by a layer of buoyant anomalies subject to strong rotation, inspired by the equatorially biased heat flux in convective planetary dynamos. We assess the role of inertial waves in generating columnar structures. At high rotation rates (or weak forcing) we find columnar flow structures that segregate helicity either side of the buoyant layer, whose axial length scale increases linearly, as predicted by the theory of low-frequency inertial waves. As the rotation rate is weakened and the magnitude of the buoyant perturbations is increased, we identify a portion of the flow which is more strongly three-dimensional. We show that the flow in this region is turbulent, and has a Rossby number above a critical value $Ro^{crit}\sim 0.4$, consistent with previous findings in rotating turbulence. We suggest that the discrepancy between the transition value found here (and in rotating turbulence experiments), and that seen in the numerical dynamos ($Ro^{crit}\sim 0.1$), is a result of a different choice of the length scale used to define the local $Ro$. We show that when a proxy for the flow length scale perpendicular to the rotation axis is used in this definition, the numerical dynamo transition lies at $Ro^{crit}\sim 0.5$. Based on this we hypothesise that inertial waves, continually launched by buoyant anomalies, sustain the columnar structures in dynamo simulations, and that the transition documented in these simulations is due to the inability of inertial waves to propagate for $Ro>Ro^{crit}$.