Saturation of Zeldovich stretch-twist-fold map dynamos

Zeldovich's stretch-twist-fold (STF) dynamo provided a breakthrough in conceptual understanding of fast dynamos, including the small-scale fluctuation dynamos. We study the evolution and saturation behaviour of two types of generalized Baker's map dynamos, which have been used to model Zeldovich's STF dynamo process. Using such maps allows one to analyse dynamos at much higher magnetic Reynolds numbers Re_M as compared to direct numerical simulations. In the two-strip map dynamo there is constant constructive folding, while the four-strip map dynamo also allows the possibility of a destructive reversal of the field. Incorporating a diffusive step parametrized by Re_M into the map, we find that the magnetic field B(x) is amplified only above a critical Re_M=R_crit ∼ 4 for both types of dynamos. The growing B(x) approaches a shape-invariant eigenfunction independent of initial conditions, whose fine structure increases with increasing Re_M. Its power spectrum M(k) displays sharp peaks reflecting the fractal nature of B(x) above the diffusive scale. We explore the saturation of these dynamos in three ways: via a renormalized reduced effective Re_M (case I) or due to a decrease in the efficiency of the field amplification by stretching, without changing the map (case IIa), or changing the map (case IIb), and a combination of both effects (case III). For case I, we show that B(x) in the saturated state, for both types of maps, approaches the marginal eigenfunction, which is obtained for Re_M=R_crit independent of the initial Re_M=R_M0. On the other hand, in case II, for the two-strip map, we show that B(x) saturates, preserving the structure of the kinematic eigenfunction. Thus the energy is transferred to larger scales in case I but remains at the smallest resistive scales in case II, as can be seen from both B(x) and M(k). For the four-strip map, B(x) oscillates with time, although with a structure similar to the kinematic eigenfunction. Interestingly, the saturated state in case III shows an intermediate behaviour, with B(x) similar to the kinematic eigenfunction at an intermediate Re_M=R_sat, with R_M0 R_sat R_crit. The R_sat value is determined by the relative importance of the increased diffusion versus the reduced stretching. These saturation properties are akin to the range of possibilities that have been discussed in the context of fluctuation dynamos.