Issue 
A&A
Volume 519, September 2010



Article Number  A80  
Number of page(s)  7  
Section  Numerical methods and codes  
DOI  https://doi.org/10.1051/00046361/200913702  
Published online  16 September 2010 
An efficient method for computing the eigenfunctions of the dynamo equation
M. Schrinner^{1,2}  D. Schmitt^{1}  J. Jiang^{1}  P. Hoyng^{3}
1  MaxPlanckInstitut für Sonnensystemforschung, MaxPlanckStr. 2,
37191 KatlenburgLindau, Germany
2  MAG (ENS/IPGP), LRA, École Normale Supérieure, 24 rue Lhomond,
75252 Paris Cedex 05, France
3  SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584
CA Utrecht, The Netherlands
Received 19 November 2009 / Accepted 19 May 2010
Abstract
Aims. We present an elegant method of determining
the eigensolutions of the induction and dynamo equations in a fluid
embedded in a vacuum.
Methods. The magnetic field is expanded in a
complete set of functions. The new method is based on the
biorthogonality of the adjoint electric current and the vector
potential with an inner product defined by a volume integral over the
fluid domain. The advantage of this method is that the velocity and the
dynamo coefficients of the induction and the dynamo equation do not
have to be differentiated and thus even numerically determined
tabulated values of the coefficients produce reasonable results.
Results. We provide test calculations and compare
with published results obtained by the classical treatment based on the
biorthogonality of the magnetic field and its adjoint. We especially
consider dynamos with meanfield coefficients determined from direct
numerical simulations of the geodynamo and compare with initial value
calculations and the full MHD simulations.
Key words: magnetohydrodynamics (MHD)  magnetic fields  methods: numerical
1 Introduction
The generation and evolution of magnetic fields in cosmic bodies like
the
planets and stars is generally thought to be governed by induction
processes
due to motions in their electrically conducting fluid interior. The
magnetic
field
is described by the induction equation
where
Here represents the velocity and the magnetic diffusivity.
In the framework of meanfield theory (e.g. Krause &
Rädler 1980; Moffatt
1978),
and
are considered as means, e.g. ensemble averaged, quantities, whereas
the action of the smallscale turbulent flow on the mean magnetic field
is parametrised by the socalled dynamo coefficients,
and .
They are, in general, tensors of second and third ranks, respectively.
We use the following compact notation of the mean field coefficients,
which include the socalled ,
,
and effects,
see e.g. Rädler (1980). Then,
the operator D reads
instead of (2), and acts on the mean magnetic field. Except for the additional and terms in the D operator, the induction and the dynamo equation are formally equivalent. Thus, the new method presented here equally applies to both.
The dynamo region is located in a flow domain V with exterior vacuum E. In this work we assume V to be either a sphere or a spherical shell. The magnetic field is continuous through the boundary and potential in E.
In kinematic dynamo theory, all coefficients (,
,
,
and
)
are assumed given and independent of the magnetic field. Thus the
dynamo equation is linear in the magnetic field and can be solved by
considering an eigenvalue problem
with eigenvalues describing the time evolution proportional to of the magnetic field .
Many studies have been made of the eigenvalues of the dynamo operator for various celestial bodies and with many forms of the dynamo coefficients (e.g.Bullard & Gellman 1954; Steenbeck & Krause 1969; Kumar & Roberts 1975; Gubbins et al. 2000; Deinzer et al. 1993; Gubbins 1973; Jiang & Wang 2007; Schubert & Zhang 2001; Deinzer & Stix 1971; Livermore & Jackson 2004; Dudley & James 1989; Roberts 1960; Jiang & Wang 2006; Roberts & Stix 1972; Roberts 1972; Schmitt & Schüssler 1989; Livermore & Jackson 2005). Often the coefficients are approximated by simple analytical functions of position, and their tensorial character is disregarded. Recently, the testfield method, developed by Schrinner et al. (2005,2007) (see also Ossendrijver et al. 2002,2001), allows one to determine all tensorial components of and directly from selfconsistent numerical simulations (Käpylä et al. 2009; Brandenburg et al. 2008). These coefficients are sometimes strongly varying functions of position. This may introduce large errors because the dynamo operator involves differentiation of the dynamo coefficients, and these are only available as numerically determined tabulated values.
In this paper we present a new method that does not require differentiation, so it is also applicable to numerically determined dynamo coefficients. The method is based on the biorthogonality of the electric current and the vector potential with an inner product defined by a volume integral over the fluid. This property has already been noted by Rädler & Bräuer (1987), Hoyng (1988), Fuchs et al. (1993), Hoyng & Schutgens (1995), and Rädler et al. (2002).
The method is described in detail in Sect. 2. Extensive test calculations have been performed and compared with published results by other eigenvalue methods. Some of these tests are presented in Sect. 3. In Sect. 4 we apply the new method to eigenmodes of the dynamo operator with coefficients obtained from geodynamo models. Our conclusion are drawn in Sect. 5.
2 Eigenvalue problem
We expand the field
of the dynamo in a complete set of functions
:
Here and in the following we make use of the summation convention for two identical indices. The expansion functions are often eigenfunctions of some differential operator. Since this operator is, in general, not selfadjoint, the functions are not orthogonal. This problem is handled by using the adjoint set , with the following inner product
The integration volume X can be either the whole space V+E or the fluid domain V alone. The base functions and constitute a biorthogonal set. For a given set of functions , the adjoint set depends on the choice of the integration domain X, so in principle we have two different sets , one for X=V and one for X=V+E.
Later we adopt the free magnetic decay modes, for which the base functions and their adjoints are known. But at this point there is no need to specify which set we actually use.
2.1 Biorthogonal sets
Starting from a set
and
that is biorthogonal on V+E, a
very useful biorthogonal set on V is provided by
the associated electric current
and the vector potential
where ,
with the inner product
Here we have absorbed a factor of in the definition of the current . The relation (7) is derived with the help of the vector identity and a volume integration over V+E; and go fast enough to zero at infinity. Surface integrals vanish because the field and the vector potential are continuous through . Volume integrals containing currents are restricted to V since in E. The inner product (7) is invariant under a gauge transformation because , as currents and their adjoints run parallel to the boundary.
Electric currents and vector potentials thus form a biorthogonal set on V. This is essential for the new eigenvalue method presented in Sect. 2.3.
2.2 Classical eigenvalue method
Inserting the expansion (5)
in the dynamo eigenvalue Eq. (4)
yields
Subsequently, we take the inner product (6) based on V with the adjoint magnetic field. This leads to
A partial integration to shift the curl from the second to the first term, as done in (7) above and used in the new method below, is not possible because the surface term need not vanish here.
We mention as an aside that the magnetic field is often decomposed in its poloidal and toroidal components (see Appendix A) after which the dynamo equation is formulated in terms of the defining scalars P and T. If the dynamo coefficients possess certain symmetry properties, the solutions can be split into two independent subsets, describing magnetic fields symmetric and antisymmetric with respect to the equator.
2.3 New eigenvalue method
We start again with (8),
which we uncurl to obtain
Taking now the inner product (7) with the adjoint current results in
The gradient term drops out as discussed in Sect. 2.1 above. The corresponding adjoint functions here are different from those in Sect. 2.2 as they pertain to a different inner product.
The matrices M_{ki} and N_{ki} have the same eigenvalues . The advantage of the new method using N_{ki} in (11) instead of M_{ki} in (9) is that no differentiation of the operator D is required, so even numerically computed or tabulated values of , , and produce accurate results.
2.4 Choice of b and numerical handling of Eq. (11)
For the set of base functions, we adopt the free magnetic decay modes whose magnetic fields are known analytically in V+E in terms of the defining scalars P and T as described in Appendix A. The decay modes are continuous through and potential in E, so they satisfy the boundary conditions. They are characterised by three numbers, the radial order n, the latitudinal degree l, and the azimuthal order m.
Another advantage of the decay modes is that they are selfadjoint on V+E so that the adjoint functions are the complex conjugates and likewise . Normalisation on V+E, i.e. , is thus readily achieved, see Appendix A.
The computation of the matrix elements N_{ki}
is now straightforward. Once we know the matrix elements, the
eigenvalue problem (11)
is solved
numerically using LAPACK routines (http://www.netlib.org/lapack), and
we obtain
the eigenvalues
and eigenvectors ,
such that
(12) 
is eigenfunction of with eigenvalue . Each mode k contains, in general, a mixture of n, l, and mvalues.
Table 1: Eigenvalues of the fundamental dipolar mode and the fifth and tenth overtones of the sphere.
In the following we consider only velocities and dynamo coefficients that are independent of azimuth , but this is not a necessary constraint. Thus each value of m can be treated separately. Although we present only results for m=0 here, we have tested and applied other values of m as well. We employ the robust GaussLegendre quadrature in r and to compute the matrix elements since the basis functions are heavily oscillatory in r for high values of n and in for high degree l. For the GaussLegendre integration we used 66 quadrature points here in the radial and 80 in the latitudinal direction, respectively. In general, this depends of course on the required resolution.
2.5 Adjoint eigenfunctions
We now show how one may construct the adjoint set of eigenfunctions of a set of eigenfunctions of the dynamo operator . Although these adjoints are not needed in the present paper, they appear in applications. For example, let be the actual magnetic field of the dynamo, then it is often advantageous to expand in dynamo eigenfunctions, i.e., . To find the coefficients c_{i}(t) we use the adjoint set, to find .
This illustrates that we need the adjoints ,
and these may be
constructed as follows. Let
be the representation of
in terms of
the selfadjoint magnetic decay modes as above. Then we
write
and ,
and we require
(13) 
Therefore in matrix notation, and we find a unique biorthogonal set. Here, indicates the Hermitean adjoint, , where T indicates the transposed and * complex conjugation.
Three important messages follow from this construction: (i) the adjoint of is , that is, the adjoint operation commutes with ; (ii) to obtain the adjoint eigenfunctions, it is not necessary to know the explicit form of the adjoint dynamo operator ; and (iii) the eigenfunctions and their adjoints have the same boundary conditions because they are a linear combination of the decay modes and their complex conjugates, respectively.
Table 2: Eigenvalues of the fundamental mode for the flow for two magnetic Reynolds numbers .
3 Test results
3.1 sphere
We first consider the socalled sphere of unit radius r_{0}=1, represented by , , and , which also can be treated analytically (Krause & Rädler 1980, Chap. 14). The eigenvalues are independent of azimuth m, and the eigenfunctions decouple in latitudinal quantum number l. For , the modes couple in radial number n, as they do between the poloidal and toroidal components.
For ^{}, the first mode is a stationary dipole, while the overtones decay with the rates given by (Hoyng & van Geffen 1993, HvG93). We successfully reproduced the fundamental mode and the overtones. In Table 1 we consider the convergence in the eigenvalues as a function of the maximum radial number for some dipolar (l=1) modes. Higher l modes behave similarly. We also reproduced the eigenfunction plots as provided by Krause & Rädler (1980).
Table 3: Eigenvalues of the first and fifth modes of the dynamo (see Sect. 3.3).
3.2 Spherical flows
As a next test, we apply the spherical stationary
(MDJ) flow of Livermore & Jackson
(2004) which is given by
=  
(14) 
with and such that the rms poloidal to toroidal energy ratio is 0.5, and the flow has an rms value of u_{0}. P_{l}^{m} and are defined in Appendix A.
Like Livermore & Jackson (2005) we consider the axisymmetric (m=0) and equatorially antisymmetric magnetic field solution for a unit sphere (r_{0}=1) embedded in a vacuum. Table 2 shows the convergence in the eigenvalue with the largest real part as a function of truncation and for two magnetic Reynolds numbers and , together with the converged values given by Livermore & Jackson (2005) (LJ05). There is a difference of about one permille between their value and ours for .
3.3 and dynamos
We reproduced the critical dynamo numbers further for the dipolar (l=1) mode of an isotropic dynamo with and N=1,2 in a spherical shell of inner and outer radius r_{i} and r_{0} with r_{0}r_{i}=1 and r_{i}/r_{0}=0.35 and 0.8 surrounded by a vacuum and either an insulating or a conducting inner core, as reported in Table 2 of Schubert & Zhang (2001). With this test we treated in particular two different aspect ratios of a spherical shell (a thick and a thin one) and two different molecular diffusivities (insulating or conducting) of the inner core.
Finally we applied our method to an dynamo of Jiang & Wang (2006) who employ the classical eigenvalue treatment for the poloidal and toroidal scalars P and T expanded in spherical harmonics in the angular coordinates and in Chebychev polynomials in rdirection. We set , , , except with r_{i}=0.5, r_{0}=1, embedded in a vacuum inside and outside, and and . Some results obtained by the new and the classical methods are compiled in Table 3. Numbers in parentheses are the real and imaginary parts of complex eigenvalues. The real part denotes the growth rate, the imaginary part the frequency of the mode in units of . The modes are axisymmetric (m=0), the fundamental mode is monotonously growing and symmetric (indicated by S) with respect to the equator, and the fourth overtone is damped, oscillatory and antisymmetric (indicated by A). Modes with higher m are more strongly damped. refers to the maximum radial number of the decay modes (spherical Bessel functions) for the new method and to the maximum degree of the Chebychev polynomials for the code of Jiang & Wang (2006), respectively. Since the modes have smaller length scales in latitudinal than in radial direction, higher values of l than of n are required for convergence. We find remarkably similar convergence of the eigenvalues for both methods. This also applies to modes with higher m. Of course we have also verified that the eigenfunctions obtained with the two methods are identical.
4 Geodynamo models
Having proven that the new method works correctly and efficiently, we now apply it to determine the eigensolutions of the dynamo operator with meanfield coefficients obtained from selfconsistent numerical simulations of the geodynamo. For a recent review of numerical geodynamo simulations, see Christensen & Wicht (2007). Schrinner et al. (2007) developed an efficient method of calculating all tensorial meanfield coefficients and and compared the results of meanfield and direct numerical simulations of the geodynamo. We plan to use the eigenmodes of the dynamo equation to decompose the magnetic field of the numerical simulations and to determine the statistical properties of the mode coefficients (Hoyng 2009) to analyse the working of the geodynamo.
4.1 Benchmark dynamo
We examine a quasisteady geodynamo model which has been used before as a numerical benchmark dynamo (Christensen et al. 2001, case 1). The governing parameters are Ekman number E=10^{3}, Rayleigh number Ra = 100, Prandtl number Pr=1, and magnetic Prandtl number Pm = 5. The convection pattern is columnar with a natural 4fold azimuthal symmetry and is stationary except for an azimuthal drift. The intensity of the fluid motion is characterised by a magnetic Reynolds number of , defined with a characteristic flow velocity, the thickness of the convecting shell, and the molecular magnetic diffusivity. The magnetic energy density exceeds the kinetic one by a factor of 20.
In Schrinner et al. (2007), the meanfield coefficients are derived from the numerical simulation. We solved the dynamo equation with these meanfield coefficients by the new method and obtained the eigenvalues and eigenfunctions. Since the coefficients are spatially variable to a considerable degree, converged solutions require high truncation levels in n and l. The eigenvalues of the first two modes are shown in Table 4. Beyond and , the eigensolution of the first mode does not change significantly and is displayed in Fig. 1. The convergence of the second mode requires a larger of about 32. The results for high values of may be affected by the spatial variation of the meanfield coefficients and would require more than 66 radial quadrature points to compute the matrix elements.
A comparison of Fig. 1 with its counterpart Fig. 10 of Schrinner et al. (2007) shows that the field of the antisymmetric fundamental mode resembles the field of an initialvalue meanfield dynamo calculation remarkably well as it does the axisymmetric component of the direct numerical simulation. The mode here grows slightly with a rate around , the field of the initial value calculation decays slightly with a rate of approximately  , while the solution of the direct numerical simulation is stationary^{}. Here L=r_{0}r_{i}=1 is the thickness of the spherical shell. The difference in these rates between the eigenvalue and initial value calculation comes from the higher numerical diffusivity of the latter at the chosen resolution of 33 radial and 80 latitudinal grid points. The difference is actually small, much less than one effective decay rate, because the relevant turbulent diffusivity, described by the coefficient with values up to , is much higher than the molecular one.
Besides the true physical eigenmodes, we find growing unphysical spurious eigenmodes. Their eigenvalues depend strongly on the resolution, and their eigenfunctions are highly structured. We attribute their appearance to a locally confined inappropriate parametrisation of the mean electromotive force by the meanfield coefficients and (Schrinner et al. 2007). The spurious modes are present neither in the initial value calculation nor in the following example of a timedependent dynamo, because of a higher numerical and molecular diffusivity, respectively.
4.2 A timedependent dynamo in the columnar regime
The next example has stronger forcing with parameters E=10^{4}, Ra=334, Pr=1, and Pm=2. The numerical simulation by Olson et al. (1999, case 2) shows a highly timedependent, but still dominantly columnar convection characterised by a magnetic Reynolds number of . The magnetic energy exceeds the kinetic energy by a factor of three. The magnetic field has a strong axial dipole contribution. Although chaotically timedependent, the velocity field is symmetric and the magnetic field antisymmetric with respect to the equatorial plane.
Table 4: Eigenvalues in units of of the first two eigenmodes of the benchmark dynamo.
Figure 1: Magnetic field structure of the fundamental antisymmetric eigenmode for the benchmark dynamo. Compare with Fig. 10 of Schrinner et al. (2007). For each plot the grey scale is separately adjusted to its maximum modulus with white as negative and black as positive. The contour lines correspond to 0.1, 0.3, 0.5, 0.7, and 0.9 of the maximum modulus. 

Open with DEXTER 
The meanfield coefficients are obtained as before by the testfield method of Schrinner et al. (2007). The coefficients are now of course also highly timedependent. A time average yields coefficients that roughly resemble those for the benchmark dynamo, although there are differences in some profiles and amplitudes.
For the timeaveraged dynamo operator the eigenvalues of the first two antisymmetric eigenmodes for various values of and are shown in Table 5. It seems that a value of is sufficient for convergence, while is needed. Figure 2 shows the eigenfunctions of these modes. The eigensolutions for and , values for and , are symmetric with respect to the equator.
Table 5: Eigenvalues in units of of the first two antisymmetric eigenmodes of the temporally averaged dynamo operator obtained from the timedependent dynamo (case 2, Sect. 4.2).
Figure 2: Radial component of the first two antisymmetric eigenmodes for the case2 dynamo. Left: fundamental mode; middle: real part of the first overtone; right: imaginary part of the first overtone. Grey scales and contours as in Fig. 1. 

Open with DEXTER 
An initialvalue, meanfield dynamo calculation with the same mean velocity and dynamo coefficients shows a slighly decaying solution with a decay rate of approximately which is to be compared with the eigenvalue of the fundamental mode. Again, the turbulent diffusivity exceeds the molecular one by a factor of up to 23 in this case. The difference in the decay rates is therefore much less than one effective decay rate. As for the benchmark dynamo, the profile of the antisymmetric fundamental mode is again remarkably similar to the solution of the initial value calculation and to the axisymmetric component of the direct numerical simulation.
A decomposition of the actual magnetic field of the simulation by Olson et al. (1999, case 2) in eigenfunctions of the timeaveraged dynamo operator, i.e., , shows that the antisymmetric fundamental mode contributes to about 75 percent and, together with the first antisymmetric overtone (see Table 5 and Fig. 2), to about 85 percent of the total magnetic energy. The variability in time of the magnetic field of the direct numerical simulation is reflected in the variability of the expansion coefficients. More details are presented in Schrinner et al. (2010).
5 Conclusions and outlook
We presented a new method for computing the eigenvalues and eigenfunctions of the induction and the dynamo equation. The method is based on the biorthogonality of the adjoint electric current and the vector potential with an inner product defined by a volume integral over the fluid domain. The advantage of the method is that the velocity and dynamo coefficients do not have to be differentiated. The method is therefore wellsuited for spatially strongly variable dynamo coefficients.
We tested the new method against the classical treatment and proved that it works correctly and efficiently. We applied it to two cases with dynamo coefficients derived from direct numerical simulations of the geodynamo. The obtained dynamo eigenmodes are promising candidates for decomposing the magnetic field of the numerical simulations and for analysing the statistical properties of the mode coefficients as proposed by Hoyng (2009).
AcknowledgementsWe thank Ulrich Christensen, Johannes Wicht, and Robert Cameron for many useful discussions and support. We further thank the referee, Matthias Rheinhardt, for his detailed comments that helped to improve the paper.
Appendix A: Free magnetic decay modes in a sphere or spherical shell embedded in vacuum
We decompose the magnetic field in its poloidal and toroidal
components
(A.1) 
with defining scalars and and unit vector in radial direction in spherical coordinates . The equation of free magnetic decay
(A.2) 
with constant magnetic diffusivity then reads
(A.3) 
where is the horizontal Laplacian
(A.4) 
The solutions are the free magnetic decay modes
(A.5)  
(A.6) 
with x=r/r_{0} where r_{0} is the radius of the sphere. The growth rates are given by
(A.7) 
and are independent of the azimuthal degree m. The constants p_{ln} and t_{ln} are
(A.8) 
where j_{l,n} is the nth zero of j_{l}. The Y_{l}^{m} are the spherical harmonics and normalised to unity by taking
(A.9) 
using Ferrer's definition of the Legendre functions of first kind P_{l}^{m} with degree l and order m.
For a sphere embedded in vacuum the radial functions are given by
(A.10)  
(A.11) 
with the spherical Bessel functions of first kind j_{l}. This ensures regularity in the origin of the sphere, vanishing toroidal component at its outer boundary and smooth transition of the poloidal component to a potential field in the vacuum outside.
For a spherical shell with inner radius r_{i}
(x_{i}=r_{i}/r_{0})
and outer radius
r_{0} (x_{0}=1)
embedded in vacuum the radial functions inside the shell are
given by
f_{l}(p_{ln}x)=j_{l}(p_{ln}x)y_{l}(p_{ln}x)j_{l+1}(p_{ln}x_{i})/y_{l+1}(p_{ln}x_{i})  (A.12) 
and
g_{l}(t_{ln}x)=j_{l}(t_{ln}x)y_{l}(t_{ln}x)j_{l}(t_{ln}x_{i})/y_{l}(t_{ln}x_{i}),  (A.13) 
and the constants in the arguments are the roots of
j_{l+1}(p_{ln}x_{i})y_{l1}(p_{ln})j_{l1}(p_{ln})y_{l+1}(p_{ln}x_{i})=0  (A.14) 
for p_{ln} and of
j_{l}(t_{ln})y_{l}(t_{ln}x_{i})j_{l}(t_{ln}x_{i})y_{l}(t_{ln})=0  (A.15) 
for t_{ln}. Here y_{l} are the spherical Bessel functions of second kind.
The magnetic field of the decay modes
is obtained by inserting the
spatial parts of the defining scalars P_{lmn}
and T_{lmn},
respectively,
into Eq. (A.1). Here we have comprised the three indices into
one. The decay modes are selfadjoint on V+E,
so that the adjoint functions are obtained simply by complex
conjugation:
and likewise .
Normalisation on V+E, i.e.,
,
is thus straightforward.
For a unit sphere the radial functions are normalised to unity by
scaling the f_{l}
with
(A.16) 
and the g_{l} with
(A.17) 
For a spherical shell the normalisation constants are more lengthy expressions, which we suppress here.
The free magnetic decay modes form a complete and orthogonal set of functions, and they obey the boundary conditions of the magnetic field between the dynamo volume V and the exterior vacuum E.
We mention for completeness that the poloidal decay modes are not selfadjoint on V, i.e., . If we like to work with an inner product defined on V, the adjoint functions can be constructed by requiring , similar to the one described in Sect. 2.5.
References
 Brandenburg, A., Rädler, K.H., Rheinhardt, M., & Subramanian, K. 2008, ApJ, 687, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Bullard, E., & Gellman, H. 1954, Phil. Trans. R. Soc. Lond. A, 247, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Christensen, U. R., & Wicht, J. 2007, in Treatise of Geophysics, Core Dynamics, ed. G. Schubert, 8, 245 (Amsterdam: Elsevier) [Google Scholar]
 Christensen, U. R., Aubert, J., Cardin, P., et al. 2001, Phys. Earth Planet. Inter., 128, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Deinzer, W., & Stix, M. 1971, A&A, 12, 11 [NASA ADS] [Google Scholar]
 Deinzer, W., Grosser, H., & Schmitt, D. 1993, A&A, 273, 405 [NASA ADS] [Google Scholar]
 Dudley, M., & James, R. 1989, Proc. R. Soc. Lond. A, 425, 407 [NASA ADS] [CrossRef] [Google Scholar]
 Fuchs, H., Rädler, K.H., & Schüler, M. 1993, IAU Symp., 157, 129 [NASA ADS] [Google Scholar]
 Gubbins, D. 1973, Phil. Trans. R. Soc. Lond. A, 274, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Gubbins, D., Barber, C., Gibbons, S., & Love, J. 2000, Proc. R. Soc. Lond. A, 456, 1333 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyng, P. 1988, ApJ, 332, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyng, P. 2009, Phys. Rev. E, 79, 046320 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyng, P., & van Geffen, J. H. G. M. 1993, Geophys. Astrophys. Fluid Dynamics, 68, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyng, P., & Schutgens, N. A. J. 1995, A&A, 293, 777 [NASA ADS] [Google Scholar]
 Jiang, J., & Wang, J. X. 2006, Chin. J. Astron. Astrophys., 6, 227 [NASA ADS] [CrossRef] [Google Scholar]
 Jiang, J., & Wang, J. X. 2007, MNRAS, 377, 711 [NASA ADS] [CrossRef] [Google Scholar]
 Käpylä, P. J., Korpi, M. J., & Brandenburg, A. 2009, A&A, 500, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Krause, F., & Rädler, K.H. 1980, Meanfield magnetophydrodynamics and dynamo theory, (Oxford: Pergamon Press) [Google Scholar]
 Kumar, S., & Roberts, P. H. 1975, Proc. R. Soc. Lond. A, 344, 235 [NASA ADS] [CrossRef] [Google Scholar]
 Livermore, P. W., & Jackson, A. 2004, Proc. R. Soc. Lond. A, 460, 1453 [NASA ADS] [CrossRef] [Google Scholar]
 Livermore, P. W., & Jackson, A. 2005, Geophys. Astrophys. Fluid Dynamics, 99, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Moffatt, H. K. 1978, Magnetic field generation in electrically conducting fluids (Cambridge: Cambridge University Press) [Google Scholar]
 Olson, P., Christensen, U., & Glatzmaier, G. A. 1999, J. Geophys. Res., 104, 10383 [NASA ADS] [CrossRef] [Google Scholar]
 Ossendrijver, M., Stix, M., & Brandenburg, A. 2001, A&A, 376, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ossendrijver, M., Stix, M., Brandenburg, A., & Rüdiger, G. 2002, A&A, 394, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rädler, K.H. 1980, Astron. Nachr., 301, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Rädler, K.H., & Bräuer, H.J. 1987, Astron. Nachr., 308, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Rädler, K.H., Rheinhardt, M., Apstein, E., & Fuchs, H. 2002, Magnetohydrodynamics, 38, 41 [NASA ADS] [EDP Sciences] [Google Scholar]
 Roberts, P. H. 1960, J. Math. Analysis Applic., 1, 195 [CrossRef] [Google Scholar]
 Roberts, P. H. 1972, Phil. Trans. R. Soc. Lond. A, 272, 663 [NASA ADS] [CrossRef] [Google Scholar]
 Roberts, P. H., & Stix, M. 1972, A&A, 18, 453 [NASA ADS] [Google Scholar]
 Schmitt, D., & Schüssler, M. 1989, A&A, 223, 343 [NASA ADS] [Google Scholar]
 Schrinner, M., Rädler, K.H., Schmitt, D., Rheinhardt, M., & Christensen, U. R. 2005, Astron. Nachr., 326, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Schrinner, M., Rädler, K.H., Schmitt, D., Rheinhardt, M., & Christensen, U. R. 2007, Geophys. Astrophys. Fluid Dynamics, 101, 81 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Schrinner, M., Schmitt, D., Cameron, R., & Hoyng, P. 2010, Geophys. J. Int., 182, 675 [NASA ADS] [CrossRef] [Google Scholar]
 Schubert, G., & Zhang, K. 2001, ApJ, 557, 930 [NASA ADS] [CrossRef] [Google Scholar]
 Steenbeck, M., & Krause, F. 1969, Astron. Nachr., 291, 49 and 271 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
 ...^{}
 The numerical value of is equal to the first zero of the spherical Bessel function j_{1}.
 ... stationary^{}
 In Schrinner et al. (2007) the mean flow entered with a sign error into the initial value calculation, leading to a stronger decay of  . We apologise and correct this value here.
All Tables
Table 1: Eigenvalues of the fundamental dipolar mode and the fifth and tenth overtones of the sphere.
Table 2: Eigenvalues of the fundamental mode for the flow for two magnetic Reynolds numbers .
Table 3: Eigenvalues of the first and fifth modes of the dynamo (see Sect. 3.3).
Table 4: Eigenvalues in units of of the first two eigenmodes of the benchmark dynamo.
Table 5: Eigenvalues in units of of the first two antisymmetric eigenmodes of the temporally averaged dynamo operator obtained from the timedependent dynamo (case 2, Sect. 4.2).
All Figures
Figure 1: Magnetic field structure of the fundamental antisymmetric eigenmode for the benchmark dynamo. Compare with Fig. 10 of Schrinner et al. (2007). For each plot the grey scale is separately adjusted to its maximum modulus with white as negative and black as positive. The contour lines correspond to 0.1, 0.3, 0.5, 0.7, and 0.9 of the maximum modulus. 

Open with DEXTER  
In the text 
Figure 2: Radial component of the first two antisymmetric eigenmodes for the case2 dynamo. Left: fundamental mode; middle: real part of the first overtone; right: imaginary part of the first overtone. Grey scales and contours as in Fig. 1. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.