A&A 397, 393-399 (2003)
DOI: 10.1051/0004-6361:20021568
V. Archontis^{1} - S. B. F. Dorch^{2,3} - Å. Nordlund^{3}
1 - Instituto de Astrofisica de Canarias, Via Lactea s/n, 38200
La Laguna, Spain
2 -
The Institute for Solar Physics of the
Royal Swedish Academy of Sciences,
SCFAB, 10691 Stockholm, Sweden
3 -
The Niels Bohr Institute for Astronomy, Physics and Geophysics,
Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark
Received 2 May 2002 / Accepted 21 October 2002
Abstract
Numerical simulations of kinematic dynamo action in
steady and three-dimensional ABC flows are presented with special
focus on the difference in growth rates between cases with single
and multiple periods of the prescribed velocity field. It is found
that the difference in growth rate (apart from a trivial factor
stemming from a scaling of the rate of strain with the wavenumber
of the velocity field) is due to differences in the recycling of
the weakest part of the magnetic field. The single wavelength
classical ABC-flow experiments impose stronger symmetry
requirements, which results in a suppression of the growth rate.
The experiments with larger wave number achieve growth rates that
are more compatible with the turn-over time scale by breaking the
symmetry of the resulting dynamo-generated magnetic field.
Differences in topology in cases with and without stagnation
points in the imposed velocity field are also investigated, and it
is found that the cigar-like structures that develop in the
classical A=B=C dynamos are replaced by ribbon structures in
cases where the flow is without stagnation points.
Key words: diffusion - magnetic fields - magnetohydrodynamics (MHD) - turbulence - plasmas
In the simplest formulation of the kinematic dynamo problem, an initially weak magnetic field is amplified due to the motions of a prescribed, stationary velocity field. Flows that have chaotic transport properties are prime candidates for fast dynamos. They have the ability to amplify the magnetic field exponentially by stretching the magnetic field lines. One such class of flows that is often invoked in fast dynamos studies is the ABC flows (named after Arnold, Beltrami & Childress). A first detailed report of the kinematic dynamo action of ABC flows may be found in Dombre et al. (1986), who examined especially the internal symmetries of the velocity field. There is, however, no guarantee that such symmetries are preserved in the eigenmodes of the magnetic field. Some discussion of possible symmetry breaking in the eigenmodes may be found in Galloway et al. (1984) and in Galanti et al. (1992). They found no evidence of symmetry breaking when the magnetic Reynolds number is limited to very low values, although symmetries are broken for higher values of the magnetic Reynolds number. The presence of symmetry breaking was also shown by Galloway et al. (1984) when the wavenumber of the periodic ABC flow is equal to 3 and 4. The conclusion was that the physical reason for the symmetry breaking is unclear and requires further investigation.
Dynamo action in ABC-flows has also been studied to determine the growth rate of the magnetic field for different values of the magnetic Reynolds number (Arnold et al. 1983; Galloway et al. 1984; Lau & Finn 1993). Results for the case with unit wavenumber, for magnetic Reynolds number up to (Dorch 2000) showed that the growth rate saturates at a value close to 0.07. Galanti et al. (1992) found that doubling the wavenumber of the flow leads to a sharp increase of the growth rate, while for larger wavenumbers the growth rate scales with the strain rate of the flow (i.e., in proportion to the wave number).
In the current paper we address some of the issues mentioned above by means of numerical kinematic MHD simulations. The dynamics of the magnetic structures in physical space is examined in connection with the amplification mechanism which drives the dynamo action in ABC flows. These flows become hydrodynamically unstable at low kinetic Reynolds number, Re (Podvigina & Pouquet 1994). We briefly discuss the dependence of the growth rate on but defer a detailed discussion of this topic to a forthcoming paper, and concentrate on the dynamo properties of laminar ABC flows in the current paper. In particular, we investigate, the physical reason that causes the increase of the growth rate in an ABC flow for wavenumbers larger than unity. To achieve this we follow the motion of magnetic field lines in time and demonstrate that the transportation of the weak part of the field is a key factor of the dynamo process.
The numerical aspects of the simulations are described in the following section. The growth rates produced by fast kinematic dynamo action of ABC flows are presented in Sect. 3. Visualizations of the magnetic structures and field lines in physical space are shown in Sect. 4. Dynamo action produced by flows without stagnation points is studied in Sect. 5. Section 6 contains the overall conclusions of the numerical simulations presented here.
Our aim is to study kinematic dynamo action in ABC flows. The time
evolution of B is governed by the (dimensionless) induction
equation:
Equation (1) is solved numerically on a staggered mesh using array valued functions to evaluate the space and time derivatives. A third-order predictor-corrector method (Hyman 1979) is used for the time-stepping. The applied numerical scheme (by Galsgaard, Nordlund and others, see Galsgaard & Nordlund 1997; Nordlund et al. 1992) ensures Eq. (2). Periodicity of is assumed over a three-dimensional computational domain. The initial magnetic field is a random perturbation with an amplitude of 10^{-5}.
Previous numerical simulations (e.g. Arnold et al. 1983; Galloway et al. 1984; Galanti et al. 1992; Dorch 2000) have shown that the normal (A=B=C=k=1) ABC flow very likely acts as a fast dynamo with a growth rate close to 0.077 (corresponding to a time-scale of 0.25 turn-over times). When the wavenumber k is larger than unity (where k=1 corresponds to the largest scale of in the computational domain) the growth rate scales approximately with the rate of strain. These results are confirmed in our simulations and we continue the study of the growth rate of the magnetic field for higher values of .
We first discuss the case with A=B=C=1 and k=2 before turning our attention to higher wavenumbers: As shown in Fig. 1 the magnetic field is amplified exponentially without oscillations. Initially, several modes may be present, associated with oscillatory growth of the magnetic field (as in the case of a normal ABC flow, Dorch 2000), but eventually the mode with the highest growth rate is the only one that remains in the solution (i.e. the dominant mode).
Figure 1: Temporal evolution of the magnetic energy for the kinematic ABC dynamo with A=B=C=1 and k=2 for and . | |
Open with DEXTER |
For low
the growth rate increases with
and reaches its maximum value at
.
This value
of the magnetic Reynolds number belongs to the first dynamo window
for the normal ABC flow dynamo, and corresponds to a maximum
growth rate of the magnetic field in this window (Galanti et al. 1992). With k=2 the lowest
possible
for dynamo action is found to be approximately equal to 2. For
large
the growth rate decreases slightly (in
contrast to the low
results of Galanti et al. 1992) but always remains positive and bounded
away from zero (Fig. 2). The largest value of
considered here is
and the
corresponding numerical resolution is 128^{3}. In the high
regime the growth rate saturates at a value
close to 0.3.
Figure 2: Growth rate of the magnetic energy (normalized by the rate of strain) versus the magnetic Reynolds number for the kinematic ABCdynamo with A=B=C=1 and k=2. | |
Open with DEXTER |
The evidence indicates that the ABC flow with k=2 is a fast dynamo, since the growth rate seems to saturate at a non-zero value for high . The reason for the enhancement of the growth rate when k=2 is discussed in the following section.
The results for the growth rate of the magnetic energy at
different values of magnetic Reynolds number when k=3 and k=4are shown in Fig. 3.
Figure 3: Two panels showing the growth rate of the magnetic energy normalized by the rate of strain for the kinematic ABC dynamo with k=3 (left panel) and k=4 (right panel). | |
Open with DEXTER |
Figure 4: Temporal evolution of the magnetic energy for the normal ABC flow with (top). Comparison with over a long stretch of the kinematic regime (bottom). | |
Open with DEXTER |
The normal ABC flow loses its stability at moderate values of Re and becomes increasingly more complex as Re increases (Podvigina & Pouquet 1994). However, as we discuss in more detail in a forthcoming paper, the growth rate is not significantly different in the nonlinear flows that develop at high Re, provided that the forcing is controlled so as to keep the turnover times approximately the same. An experiment with at 40^{3} resolution shows that the magnetic field is amplified exponentially with an average growth rate close to 0.064 (Fig. 4, top). Doubling the values of Reynolds numbers and of the resolution, the temporal evolution of the magnetic energy is examined and it is found that, apart from a transient increase of the magnetic energy as the structures shrink in size, the growth rate is very similar (Fig. 4, bottom). This result indicates that the growth rate is mainly determined by the large scale properties of the flow (in particular the turn over time), and we expect that the only important change in the magnetic field structure with increasing Reynolds number is the development of more fine structure.
An experiment with
and a numerical
resolution of 64^{3} was performed in order to study the
dynamics of the magnetic structures in physical space when the
wavenumber k of the flow is equal to 2.
The computational domain
then contains eight periodic cells, where each one (initially)
reflects the geometrical properties of a standard ABC flow. The
magnetic field evolves according to the induction equation Eq. (1) and flux "cigars" are formed at the
so-called -type stagnation points (Dombre et al. 1986). These are points (stable manifolds of the
flow) with converging streamlines in the plane of the stagnation
point and the field is rapidly advected to these regions (Fig. 5). The flux cigars are aligned along the axis of
divergence through the -type stagnation points.
Figure 5: Magnetic field strength isosurfaces (dark) and stagnation points of the flow (light). The magnetic field is concentrated along flux cigars centered at the -type stagnation points of the flow. The stagnation points with no flux cigars centered at them are the -type points. In each cell three cigars form a triangle around a -type stagnation point and another one is pointing towards it. The magnetic field strength of the isosurfaces is of the peak value in the snapshot. | |
Open with DEXTER |
Figure 6: A snapshot showing the structure of the magnetic field when the fastest growing mode appears and the symmetry is not conserved. Each cell consists of flux cigars which have different size and different field strength. | |
Open with DEXTER |
At the -type stagnation points (unstable manifolds of the flow) the field is weaker, and takes the form of flux sheets on both sides of the plane of divergence. The flux sheets are formed because of the converging flow along the axis of divergence from the -type stagnation points. The flux cigars thus appear to be connected by flux sheets through heteroclinic stream lines.
When the symmetry breaking appears, and the fastest growing eigenmode takes over, the topology of the magnetic field changes. The fastest growing mode does not follow the symmetric configuration of a standard ABC flow - the strong magnetic field is still elongated into flux "cigars'' but their size and field strength vary from cell to cell (Fig. 6).
The rate of growth (and decay) of individual flux cigars is
different, rather than being uniform as in the normal k=1 case.
The magnetic field for the latter case is recycled inside its cell
because of the exact symmetry. The time evolution of the magnetic
energy follows a periodic (oscillatory) behavior which is
associated with the conservation of symmetry in the fastest
growing eigenmode (Dorch 2000) in contrast to the
dominant mode for the k=2 case which is a symmetry breaking mode
(Galloway et al. 1984).
Figure 7: Visualization of magnetic field lines and flux cigars when the dominant mode starts to take over. The stretching of the field lines at the corner of a triangle at the bottom-left cell is obvious. These field lines get pressed up against another cigar in a neighboring cell (up-right). | |
Open with DEXTER |
Figure 8: Folding of week flux sheets at the plane of a -type stagnation point. The field lines (top of the figure) are coming from a neighboring flux cigar and change direction across the plane. | |
Open with DEXTER |
In view of the role played by the stagnation points in the normal ABC flow dynamo, it is natural to ask whether they are a necessary ingredient for fast dynamo action. When the parameters A, B and C are varied it is possible for the stagnation points to be created or destroyed (Dombre et at. 1986). When A=5, B=2, C=2 or A=5, B=2, C=1 (abbreviated as the 5:2:2and 5:2:1 case, respectively) the flow has no stagnation points. Numerical simulations were performed for both cases and the results obtained are similar. Therefore, only the former case is discussed in detail.
The 5:2:2 flow has four points where the velocity is close to zero (but not zero, hence there are no stagnation points). The disappearance of the stagnation points means that there are no longer any heteroclinic stream lines that connect the points of minimum velocity. However, there are stream lines that converge to these points and stream lines that diverge from these points (Fig. 9). The stream lines of the 5:2:2flow do not follow any specific symmetry as in the normal ABC flow (which has a three-fold symmetry of stream lines around the stagnation points) and the topology of the flow becomes more complex.
The importance of these minimum velocity points is that they
coincide with local extrema of the "perpendicular convergence" of
the local motion to a magnetic field line
Figure 9: The bead-shaped isosurfaces show the position of the points with minimum velocity of the 5:2:2 flow. Visualization of stream lines of the flow shows that they diverge when they reach the neighborhood of other points with minimum velocity. | |
Open with DEXTER |
Numerical experiments were performed with up to 3200 and numerical resolution up to 240^{3} to test whether the 5:2:2 case has vigorous dynamo properties in the kinematic regime. The initial seed field was a random perturbation with an amplitude of 10^{-5}. It was found that the magnetic field is amplified exponentially and the growth rate increases with as illustrated in Fig. 10.
Figure 10: Growth rate of the magnetic energy versus the magnetic Reynolds number for the 5:2:2 flow. | |
Open with DEXTER |
The 5:2:1 case was also examined at high magnetic Reynolds number. The highest value of considered here is . The growth rate seems to saturate at a value of about 0.67 as illustrated in Fig. 11.
As mentioned above, the 5:2:2 flow has four points with minimum velocity that coincide with local maxima of the function in Eq. (4).
When the induction equation is evolved from a weak random field, flux ribbons arise close to the minimum velocity points. The thickness of the ribbons decreases as increases, presumably scaling as , similarly to the flux cigars for the standard ABC flow. The sign (polarity) of their magnetic field changes from ribbon to ribbon. Thus, field lines with opposite polarity may reconnect when two ribbons are close together. The magnetic flux is released through the reconnected field lines (Fig. 12). However, the magnetic ribbons are not confined at certain points but are advected by the flow, and field lines with the same polarity could also be found to approach each other.
Figure 11: Growth rate of the magnetic energy versus the magnetic Reynolds number for the 5:2:1 flow. | |
Open with DEXTER |
Figure 12: The magnetic field takes the form of flux ribbons when A:B:C=5:2:2. The reconnection of field lines is shown when ribbons with opposite polarity field are found close together. | |
Open with DEXTER |
The process responsible for driving the amplification mechanism of the mode with the largest growth is the constructive folding of the ribbons.
Two of the most important properties of the 5:2:2 flow is the exponential stretching of the field lines and the way in which the magnetic field is advected and the field lines are packed. Neighboring regions with magnetic field of opposite polarity are destroyed whereas regions with like polarity flux are conserved. Therefore the rate of field generation depends on the global folding of the field lines.
Figure 13: A snapshot showing the merge of the ribbons (dark) with like polarity in the neighborhood of the minimum velocity points (light). | |
Open with DEXTER |
The flow folds the field so as to bring together flux in the same direction. Thus, stretched flux ribbons with the same polarity come into alignment and merge. The constructive folding of the sheets occurs in the vicinity of the minimum velocity points (Fig. 13). Due to the stretching ability of the flow, which is maximum close to these points, the field lines are stretched in these regions.
The flow advects the field so as to bring together stretched field lines with the same polarity. They come into alignment and are folded constructively. These field lines are continuously stretched and also advected by the flow, and move to the vicinity of another neighboring point with minimum velocity, where they fold with new stretched field lines of the same polarity. The folding is again constructive.
Due to the stretching and constructive reinforcement of the magnetic field discussed above, the magnetic energy increases each time the process is repeated.
The amplification of the dynamo-generated field is associated with a combination of stretching, twisting, folding and reconnection of magnetic field lines. The importance of the weak part of the field in the dynamo process is emphasized when the wavenumber of the flow increases. Symmetry breaking allows the transportation of stretched weak magnetic field from cell to cell, where it gets packed against the local flux structures ("cigars"). Most of the space in the bulk part of the flow has passive advection, which does not change at increasing magnetic Reynolds number. The evidence indicates that the invariant properties of the flows (i.e. stretching ability) do not change in the limit of arbitrarily small magnetic diffusivity and the growth rate of the normal ABC flow reaches an asymptotic value (for high ) even after the onset of its hydrodynamic instability.
The study of kinematic dynamo action in ABC flows is, from a physical point of view, simple enough since it does not contain the complexity of non-linear back-reaction of the magnetic field. A fully non-linear study of the dynamo problem connects the initial phase of the growth of a weak magnetization state with the saturation of the dynamo generated field.
We have demonstrated that the fast dynamo action reflects the same mechanisms at work whatever the values of the parameters (such as and k) are. The main mechanism is the stretching of the weak magnetic field. This is a first step towards the understanding of the processes at work in a turbulent environment and during the saturation state of a dynamo where the flow has still a good grip on the weakest part of the magnetic field.
It is of course also of great interest to perform dynamo simulations in more realistic (e.g. spherical) geometries, and to include more physical processes such as differential rotation and convection. The advantage of the simulations presented in the current paper is that the individual mechanisms are easily identified and can be studied in great detail. Thus, kinematic dynamo experiments such as these can provide very useful clues for the general theory of astrophysical dynamos.
Acknowledgements
VA thanks the EU-TMR for support through a Marie Curie Fellowship. SBFD was supported through an EU-TMR grant to the European Solar Magnetometry Network. The work of ÅN was supported in part by the Danish National Research Foundation, through its establishment of the Theoretical Astrophysics Center. Computing time at the UNIC computing center was provided by the Danish Natural Science Research Council.