A nonlinear dynamo wave riding on a spatially varying background

A systematic asymptotic investigation of a pair of coupled nonlinear one-dimensional amplitude equations, which provide a simplified model of solar and stellar magnetic activity cycles, is presented. Specifically, an α -dynamo in a thin shell of small gap-to-radius ratio &b.epsi; (< 1) is considered, in which the -effect (the differential rotation) is prescribed but the α-effect is quenched by the finite-amplitude magnetic field. The unquenched system is characterized by a latitudinally θ-dependent dynamo number D, with a symmetric single-hump profile, which vanishes at both the pole, θ = π/2, and the equator, θ = 0, and has a maximum, D, at mid-latitude, θM = π/4. The shape D(θ)/D is fixed, so that there is only a single driving parameter D. At onset of global instability, D = DL( ): = DT + O( ), a travelling wave, of frequency ω = ωL( ): = ωT+ O( ) and wavelength O( ), is localized at a low latitude θPT(< θM); DT and ωT are constants independent of . As a consequence of the spatial separation of θPT and θM, the squared field amplitude increases linearly with the excess dynamo number D - DL in the weakly nonlinear regime, as usual, but with a large constant of proportionality dependent on some numerically small power of exp(1/ ). Whether the bifurcation is sub- or supercritical is extremely sensitive to the value of . In the nonlinear regime, the travelling wave localized at θPT at global onset expands and lies under an asymmetric envelope that vanishes smoothly at a low latitude θP but terminates abruptly on a length O( ) - comparable to the wavelength - across a front at high latitude θF. The criterion of Dee and Langer, applied to the local linear evanescent disturbance ahead of the front, determines the lowest order value of the frequency close to the global onset value ωT. The global transition is characterized by the abrupt shift of θF from θPT to θM; during that passage, D executes O( -1) oscillations of increasing magnitude about DL. Fully developed nonlinearity occurs when θF > θM. In that regime, Meunier and co-workers showed that the O(1) quantities θF - θM and (ω - ωT)/ 2/3 increase together in concert with D - DT. By analysing the detailed structure of the front of width O( ), we obtain ω correct to the higher order O( ) and show improved agreement with numerical integrations performed by Meunier and co-workers of the complete governing equations at finite .