A&A 493, 375-383 (2009)
DOI: 10.1051/0004-6361:200809565
G. Rüdiger^{1} - D. A. Shalybkov^{2}
1 - Astrophysikalisches Institut Potsdam, An der Sternwarte 16,
14482 Potsdam, Germany
2 -
A.F. Ioffe Institute for Physics and Technology, 194021, St. Petersburg, Russia
Received 12 February 2008 / Accepted 10 September 2008
Abstract
Aims. The stability of the dissipative Taylor-Couette flow with a stable axial density stratification and a prescribed azimuthal magnetic field is considered.
Methods. Global nonaxisymmetric solutions of the linearized MHD equations with toroidal magnetic field, density stratification, and differential rotation are found for both insulating and conducting cylinders.
Results. Hydrodynamic calculations for various gap widths show that flat rotation laws such as the Kepler rotation are always unstable against SRI. Quasigalactic rotation laws, however, are stable for wide gaps. The influence of a current-free toroidal magnetic field on SRI strongly depends on the magnetic Prandtl number Pm: SRI is supported by
and it is suppressed by
.
For rotation laws that are too flat, when the hydrodynamic SRI ceases, a smooth transition exists to the instability that the toroidal magnetic field produces in combination with the differential rotation. For the first time this nonaxisymmetric azimuthal magnetorotational instability (AMRI) has been computed in the presence of an axial density gradient.
If the magnetic field between the cylinders is not current-free, then the Tayler instability occurs, too. The transition from the nonmagnetic centrifugal instability to the magnetic Tayler instability in the presence of differential rotation and density stratification proves to be complex. Most spectacular is the ``ballooning'' of the stability domain by the density stratification: already a small rotation stabilizes magnetic fields against the Tayler instability.
An azimuthal component of the electromotive force for the instability only exists for density-stratified flows. The related -effect for magnetic-influenced SRI with Kepler rotation appears to be positive for negative .
Key words: accretion, accretion disks - turbulence - instabilities - magnetohydrodynamics (MHD)
This work is motivated by the theoretical and experimental progress in studies of the stratorotational instability (SRI) and the magnetorotational instability (MRI) in MHD Taylor-Couette experiments. It has been shown theoretically (Molemaker et al. 2001; Yavneh et al. 2001; Dubrulle et al. 2005; Shalybkov & Rüdiger 2005; Umurhan 2006) and in the laboratory (Le Bars & Le Gal 2007) that a combination of a centrifugal-stable nonuniform rotation law and a stable axial density stratification leads to the SRI in the Taylor-Couette flow. This instability exists only for nonaxisymmetric disturbances. On the other hand, there are also nonaxisymmetric instabilities for a combination of Rayleigh-stable rotation laws and azimuthal magnetic fields (Pitts & Tayler 1985). The question is whether the combination of density stratification, differential rotation, and toroidal fields acts as stabilizing or destabilizing or whether even new instabilities arise.
Such a combination of axial density stratification, stable rotation law, and strong toroidal magnetic field is the typical constellation in accretion disks (Ogilvie & Pringle 1996; Curry & Pudritz 1996; Papaloizou & Terquem 1997). There the rotation is Keplerian with and the toroidal field is generated from weak large-scale poloidal fields by induction due to the differential rotation. The standard case is that the resulting toroidal field strongly exceeds the amplitude of the original poloidal field if the magnetic Reynolds number of the differential rotation is much larger than unity. If this is true, then the standard MRI, i.e. the influence of the poloidal field on thestability of the differential rotation, would be of minor importance^{}. The question is whether the toroidal magnetic field can reach such high amplitudes or whether it already becomes unstable for much lower values. It is known that toroidal fields with strong electric currents become unstable (Tayler instability, ``TI'') but it is also known that, in combination with differential rotation, toroidal fields also become unstable, which are current-free in the considered domain (azimuthal MRI, ``AMRI'', see Rüdiger et al. 2007a). In the latter case, with , the questions are whether the density stratification destabilizes AMRI and/or whether the toroidal field stabilizes the SRI too strongly, so that its existence becomes basically suppressed.
We consider the interaction of the differential rotation with both an axial density stratification and a toroidal magnetic field in the simplified Taylor-Couette geometry. The density stratification is always supposed to be stable, but the magnetic field and the rotation law between the cylinders can be both stable or unstable. If, in particular, the toroidal field is Tayler-unstable, then the interaction of differential rotation, density stratification and magnetic field becomes very complex. We find stabilization and destabilization to depend strongly on the magnetic Prandtl number. Again, in experiments, for small magnetic Prandtl number, the flows are predicted to be stabilized. For galaxies and protoneutron stars (PNS, with their high magnetic Prandtl numbers, we find the opposite. The magnetic influence supports the SRI leading to even smaller Reynolds numbers than in hydrodynamics.
Figure 1: The geometry of the problem - two concentric cylinders with radii and rotating with and . is the azimuthal magnetic field which, generally, is produced by both an axial current inside the inner cylinder and an axial current through the fluid. | |
Open with DEXTER |
(2) |
Similarly, the magnetic profiles are restricted by the azimuthal component of the induction equation to
(5) |
(7) |
In the presence of a vertical density gradient (
,
it has been shown for
that the system (8) allows the angular velocity profile (1) only in the limit of slow rotation and small stratification
The linear stability problem is considered for the perturbed state of , , , P_{0}+P_{1}+P'. Using the conditions (9), the linearized system (8) takes the Boussinesq form with the coefficients only depending on the radial coordinate. Hence, a normal mode expansion of the solution can be used with F as any of the disturbed quantities.
After a normalization we find
The dimensionless numbers of the problem are the magnetic
Prandtl number Pm, the Hartmann number Ha, the Reynolds number Re
According to the Rayleigh condition, the Taylor-Couette flow
with the rotation law (1) is stable for
Figure 2: The marginal-stability line for m=1 disturbances in a container with and . The instability is stabilized for stratifications that are both too weak and too strong. The minimum Reynolds number belongs to a Froude number of 1.4 depending on the rotation law. | |
Open with DEXTER |
As we know, a stable vertical density stratification destabilizes the Taylor-Couette flow with decreasing angular velocity even beyond the Rayleigh line. This SRI is nonaxisymmetric and the most unstable mode is the mode with m=1. Shalybkov & Rüdiger (2005) find that the instability condition is instead of , which has been confirmed experimentally by Le Bars & Le Gal (2007). This finding, however, based on a restricted number of calculations.
In Fig. 2 the marginal-stability line (the line that separates stable and unstable regions) is calculated as a function of for m=1 disturbances and for a container with and with . After (15) the last value lies beyond the Rayleigh line where nonstratified flows are linearly stable.
We are looking for the minimum critical Reynolds number that exists for a buoyancy frequency , or - which is the same - . Figure 2 reveals the SRI as a delicate balance of buoyancy and centrifugal force. Both the frequencies must be approximately equal, otherwise one of the two stabilities dominates the constellation and stabilizes the flow.
The question is how this balance depends on the rotation law and on the geometry of the container. In Fig. 3 three containers are used with different gap sizes. The computations also concern various density stratifications. Three vertical lines are given for any gap size. Left is the line for the Rayleigh limit where . The right line mimics the rotation law , which is typical of galaxies ( ), and the central line represents the Kepler rotation law , i.e. , which we call the quasi-Kepler limit. The main result of the Fig. 3 is that the quasi-Kepler rotation is always unstable. This, however, is not true for the flatter quasi-galactic rotation, which becomes stable for large gaps. From the bottom plot in Fig. 3, we have to realize that the quasi-galactic rotation law with is already too flat for the possible existence of nonmagnetic SRI in galaxies.
The present calculations do not confirm the line as the limit for the SRI. This limit can be higher for small-gap containers and can be lower for wide-gap containers. The stability limit seems to approach the Rayleigh line in the limit of very wide gaps and goes to unity for very narrow gaps. We ask in the following whether such results are influenced by the existence of a toroidal magnetic field.
Figure 3: The marginal-stability lines for a small gap ( top, ), modest gap ( , middle) and wide gap ( , bottom). The vertical dotted lines denote the Rayleigh limit, the quasi-Kepler rotation, and the quasi-galactic rotation ( ). The curves represent different density stratifications with (dot-dashed), (dashed), (dotted), (solid, maximum stratification). | |
Open with DEXTER |
Without density stratification the ideal Taylor-Couette flow with imposed azimuthal magnetic field
is stable against axisymmetric disturbances, if
The stability condition
for the azimuthal magnetic field against nonaxisymmetric disturbances after Tayler (1961) is
As , the stability interval (19) is smaller than the stability interval (17). In this sense the nonaxisymmetric (``kink'') disturbances are more unstable than the axisymmetric (``sausage'') disturbances. However, the situation is much more complex in the presence of (differential) rotation, and which instability really dominates depends on the parameters of the problem.
In the present paper the influence of stable vertical density stratifications on the Taylor-Couette flow stability with imposed azimuthal magnetic fields is considered for three different cases:
The azimuthal magnetic field with (i.e. for ) is considered which is current-free between the cylinders as the simplest example of an azimuthal magnetic field that is basically stable without rotation for both axisymmetric and asymmetric disturbances.
The typical stability diagram is presented by the Figs. 4 and 5 for containers with
conducting/insulating cylinders, medium gap width (
,
and a rather flat rotation law (quasi-Kepler
rotation with
and a slightly steeper rotation law with
),
which are both stable without density stratification and without magnetic field. We find that the results hardly depend
on the conducting properties of the cylinder material.
Figure 4: The marginal-stability lines for the flow with conducting cylinders, , current-free () azimuthal magnetic field, quasi-Kepler rotation ( ), for perturbation modes with m=1 (solid), and m=2 (dotted). The curves are labeled by their magnetic Prandtl numbers. For comparison the dashed line demonstrates AMRI (m=1) for homogeneous fluids (N=0). | |
Open with DEXTER |
Figure 5: The same as Fig. 4 but for insulating cylinders with (as in Fig. 2) but without m=2 mode. | |
Open with DEXTER |
Figures 4 and 5 indicate the (dashed) marginal-stability line of a homogeneous fluid with . There are no solutions at both the vertical axis and the horizontal axis. The MHD flow is only unstable as a combination of differential rotation and the azimuthal magnetic field. For any supercritical Hartmann number, there are two critical Reynolds numbers between which the fluid is unstable. Similarly, for any supercritical Reynolds number, there are two critical Hartmann numbers between which the fluid is unstable. We called this phenomenon the Azimuthal MagnetoRotational Instability (AMRI), which scales with the magnetic Reynolds number in the same way as for the standard MRI in Taylor-Couette experiments. However, magnetic Reynolds numbers around 100 are too high to be realized in the MHD laboratory.
The situation changes with a density stratification ( ). The following results can be interpreted as either AMRI for density-stratified fluids or as the influence of an azimuthal magnetic field on SRI. There is now an instability (SRI) without any magnetic field for at m=1 and for at m=2 (for conducting cylinders, quasi-Kepler rotation law, Fig. 4) and for at m=1 (insulating cylinders, , Fig. 5). For this case, the calculations are restricted to the kink (m=1) instability.
For conducting boundaries (see Fig. 4) and for a special case ( , m=1) the radial eigenfunctions are given for (Fig. 6), which demonstrate that the SRI modes do not form boundary phenomena. The profiles are normalized, since the amplitudes have no own physical meaning in linear theory.
For the linear growth time in units of the rotation period of the inner cylinder for the solution with ,
one obtains
Figure 6: The eigenfunctions for the radial components of flow (u_{R}, left) and field (b_{R}, right) for . Solid: real part, dashed: imaginary part; the (arbitrary) amplitude is normalized. , , m=1, see Fig. 4. | |
Open with DEXTER |
The influence of the azimuthal magnetic field on the SRI strongly depends on the magnetic Prandtl number Pm. For , the SRI is stabilized by the magnetic field, but it can be supported for (for fixed viscosity see Figs. 4, 5). For the critical Reynolds number is reduced (by a factor of two) for . In this case AMRI supports the SRI. This behavior is similar to that of the standard MRI with axial magnetic field and hydrodynamically unstable rotation where for the magnetic field also supports the centrifugal instability (Rüdiger & Shalybkov 2002).
Note also that the magnetic suppression of the SRI for small Pm is rather weak. Even for , a (large) Hartmann number around 100 only leads to the small critical Reynolds number of about 500. The lines for in Figs. 4 and 5 are identical to the lines for all smaller Pm; i.e., they also hold for gallium (10^{-6}) and mercury (10^{-7}). Here all the curves scale with Re rather than Rm so that the corresponding magnetic Reynolds numbers are quite small compared with those for AMRI. As the corresponding Hartmann numbers with are also not too high, laboratory experiments with fluid metals should be possible provided a sufficiently density stratification can be produced, e.g. by heating from above (see Gellert & Rüdiger 2008).
It is also shown in Fig. 4, that for a given Pm, the critical Reynolds numbers for m=2 exceed those for m=1. The differences, however, are not too big so that, for slightly supercritical Re, several modes should be excited. This is surprising insofar as the known smoothing action of differential rotation upon modes with high m here seems to be weak.
Figure 7: AMRI for stratified fluids with between insulating cylinders ( , ). The rotation law is so flat ( ) that nonmagnetic SRI does not exist (see Fig. 3, middle). Lines are labeled by their Pm numbers. | |
Open with DEXTER |
The small Reynolds numbers and Hartmann numbers lead to the impression that the instability should be observable in the laboratory. For hydrodynamically unstable situations, we only have
a suppression of the instability by the magnetic field for small Pm. For hydrodynamically stable situations (Figs. 4 and 5) it was shown for uniform fluids that the appropriate numbers of the problem are the magnetic
Reynolds number Rm and the Lundquist number S
Stratified flows without magnetic field become stable for -values slightly over . According to Fig. 7, the flow is stable for without a magnetic field, but it becomes unstable with magnetic fields. The critical Rm are not monotonous with Pm. The values of Rm decrease with Pm for small Pm and v.v. The result can be seen in Fig. 7, which differs in this respect from the AMRI without density stratification (see Rüdiger et al. 2007a).
The solid line in Fig. 8 demonstrates that the AMRI (like the standard MRI) exists for all rotation laws with decreasing angular velocity as a function of radius (i.e. for ).
Figure 8: The minimum Reynolds numbers for insulating cylinders, , and for the kink (m=1) mode. The vertical lines are the same as in Fig. 3. The numbers on the curves are Hartmann numbers which correspond to the minimum Reynolds numbers. The solid line is AMRI without a stratification, the dot-dashed line is AMRI with stratification ( , see Fig. 7) and the dashed line is SRI without magnetic field. | |
Open with DEXTER |
According to the Figs. 4 and 5, the marginal-stability lines have always a minimum for some Hartmann number including . These minimum Reynolds numbers are plotted in Fig. 8 as functions of . The dashed line is SRI without magnetic field, which disappears for rotation laws that are too flat ( ). There is a smooth transition, however, to the instability AMRI (current-free toroidal field plus differential rotation), which has been computed here for the first time with a density gradient. The comparison with standard AMRI without yields the expected suppression by density stratification. Note, however, that for quasi-galactic rotation ( in Fig. 8), the SRI with magnetic field can be excited more easily than without magnetic field. As the figure shows, the necessary magnetic field has a Hartmann number of about 137.
Because AMRI scales with Rm rather than Re, the ``smooth'' transition interval from the hydrodynamic solutions to the magnetohydrodynamic scales with Pm^{-1}. Hence, the dot-dashed line in Fig. 8 becomes steeper and steeper for decreasing Pm.
The influence of the vertical density stratification on the nonaxisymmetric
TI (i.e.
)
is considered for rotation laws that are steep enough to be Rayleigh unstable.
Again the gap width is
,
the outer cylinder is resting (
), and the toroidal field is fairly homogeneous
()
so that only the mode m=1 is unstable without rotation.
Calculations were made for both insulating and conducting cylinders
(Fig. 9). The difference between conducting
and insulating boundary conditions proves to be only quantitative.
Figure 9: The marginal-stability lines for a flow with insulating ( top) and conducting ( bottom) cylinders. It is , (resting outer cylinder) and . The curves are labeled by their magnetic Prandtl Pm. The m=0 (dotted) mode (no dependence on Pm) and m=1 (solid) mode for homogeneous fluid and m=0 (dot-dot-dot-dashed) and m=1 (dot-dashed) modes for density-stratified fluids with . | |
Open with DEXTER |
For m=1 the critical Hartmann numbers above which the flow is unstable for are 150 for conducting cylinders and 109 for insulating cylinders. These critical numbers are not influenced by the density stratification. No critical Hartmann numbers exist for m=0; the magnetic field stabilizes this Rayleigh-mode for both homogeneous and density-stratified fluids. For homogeneous fluids the marginal-stability line for axisymmetric disturbances does not depend on the magnetic Prandtl number Pm (Shalybkov 2006).
The critical Reynolds numbers for the Rayleigh instability for do not depend on the conducting properties of the cylinders. They are 68 (m=0) and 75 (m=1) without stratification. The density-stratification stabilizes the flow, and the critical Reynolds numbers are 294 for m=0and 226 for m=1 (for ). With density-stratification, the kink mode (m=1) proves to be the most unstable mode.
Moreover, for given Pm the kink mode is the most unstable one for all Ha for . For small Pm, this mode is stabilized by weak magnetic fields before a dramatic destabilization happens for larger magnetic fields. For increasing Pm the flow becomes unstable for decreasing Reynolds numbers at a fixed Hartmann number. The opposite is also true: the smaller Pm, the stronger is the magnetic field influence and the flow becomes unstable for smaller Hartmann numbers at a fixed Reynolds number.
Note that the competition between magnetic and centrifugal instabilities can lead to a rather complex behavior as illustrated by the marginal-stability line with . Generally, however, the stability region is much larger for density-stratified fluids so that we find that a stable stratification increases the flow stability. Of particular interest is the phenomenon that, for density-stratified fluids generally, the critical Hartmann numbers are larger than without rotation and the same is true for small Pm with respect to the critical Reynolds number (``ballooning''). For large magnetic Prandtl number, however, the magnetic field basically destabilizes the Rayleigh instability. Even for weak fields, the critical Reynolds number is smaller than the Reynolds number for the nonmagnetic case.
Figure 10: The marginal-stability lines for conducting cylinders with , , , ( top), and ( bottom). The lines connect the Rayleigh instability without magnetic field with the TI without rotation. The dotted lines are for m=0; the solid lines for m=1. The inner curves are for homogeneous fluids (N=0) and the outer curves for stratified fluids with . The density stratification stabilizes in both directions. | |
Open with DEXTER |
The most complex situation appears when both the modes with m=0 and m=1 are Tayler-unstable. Only conducting cylinders with medium-sized gap ( ), outer cylinder at rest and a toroidal magnetic field with strong currents () are analyzed. The differential rotation is Rayleigh unstable without magnetic field and density stratification and the magnetic field is Tayler-unstable without rotation (for m=0 and m=1). The critical Hartmann number for is nearly equal for all cases: m=0 gives and m=1 gives . For other parameters it can be more different. Either the m=0 or the m=1 mode can be the most unstable mode. Often the m=0 mode is also the most unstable for small Ha numbers and fast differential rotation (see Rüdiger et al. 2007b).
The critical Reynolds numbers above which the rotation becomes unstable without a magnetic field () for m=0 and m=1 are called Re_{0} and Re_{1}, resp., and Re_{0} is smaller than Re_{1} for homogeneous flows with . The density stratification increases Re_{0}and Re_{1} (Fig. 10). With density-stratification, the m=1 mode becomes the most unstable mode ( ). With a magnetic field, it is also the most unstable mode almost everywhere except for and .
``Ballooning'' of the stability region is produced by the density stratification. The density stratification in combination with the basic rotation stabilizes the flow for larger magnetic fields. Under the influence of differential rotation and density stratification stronger magnetic field amplitudes prove to be stable than they were without rotation. The effect already exists for rather slow rotation rates (see Fig. 10, ). This phenomenon is a consequence of the explicit inclusion of the stable density stratification into the MHD equations. One finds that both the Rayleigh instability and the TI are stabilized by the density stratification. Without rotation, however, the density influence of the Tayler instability would remain negligibly small.
Stratorotational instability under the influence of a toroidal field is tempting to apply to the concept of the mean-field electrodynamics in turbulent fields. The nonaxisymmetric components of flow and field can be used as fluctuations, while the axisymmetric components are considered as the mean quantities. Simply the averaging procedure is the integration over the azimuth .
It is standard to express the turbulence-induced electromotive force as
Figure 11: The correlation function (23) for AMRI with conducting cylinders and with , (quasi-Kepler), , . Density stratification is zero ( top) and finite ( bottom, ). The Reynolds numbers ( top) and ( bottom) were taken from Fig. 4. | |
Open with DEXTER |
In the present paper we only consider the -effect in the frame of a linear theory where all functions are free to one and the same (complex) arbitrary parameter. This is only possible if the second-order quantities such as
are normalized with a second-order quantity. In the considered case it makes sense to form the correlation function
One can easily show that second-order correlations of quantities running with after integration over do not depend on the coordinate z. The term kz only fixes the phase where the integration starts.
The correlation function
can be estimated for fast-rotating magnetoconvection in the high conductivity limit
and with
after Krause & Rädler (1980)
so that
This cannot be a boundary effect. A boundary effect is implausible because the preferred radial direction is perpendicular to the rotation axis for cylinders, which mimics the equator in rotating spheres and there the -effect must vanish. Indeed, the antisymmetry is reduced if a density stratification is allowed (Fig. 11, bottom). Then the density gradient is now the preferred direction in the system parallel to the rotation allowing the formation of large-scale helicity (as in rotating spheres at the poles). We are thus tempted to assume that the presented calculations for the presence of a density gradient have indeed demonstrated for the existence of an -effect for the magnetic SRI.
Figure 12: The correlation function (23) for , (quasi-Kepler), , . ). The Reynolds numbers are for and for and for . The -effect is quenched for stronger magnetic fields. | |
Open with DEXTER |
Figure 12 presents similar results but for the small magnetic Prandtl number . The correlation functions provide a dominance of the positive contributions, but for higher and higher magnetic fields, the difference becomes smaller. It seems that the numerical results reflect a magnetic quenching of the -effect here.
The stability of the dissipative Taylor-Couette flow under considering the joint influence of a stable vertical density stratification is considered and an azimuthal magnetic field is considered. The problem is of interest for future laboratory experiments but also within the frame of accretion disk physics. The Kepler rotation generates strong toroidal magnetic fields dominating the poloidal components. The standard MRI that works with only axial fields may not be very relevant for the stability of the Kepler rotation law compared with the azimuthal MRI in connection with the density stratification and Tayler instability. We mainly considered nonaxisymmetric ``kink'' modes (m=1), but there are also examples where the axisymmetric modes are the most unstable ones.
We started with a discussion of the SRI without magnetic fields. For a flat rotation law, a stratification value Fr exists for which the critical Reynolds number of the rotation has a minimum (see Fig. 2). The SRI is basically stabilized both by too weak or too strong stratification. The instability only appears if the characteristic buoyancy time approaches the rotation period. Also, if the rotation is too flat, the instability disappears.
The limiting ratio strongly depends on the gap width (see Fig. 3). For small gaps, rotation laws with even prove to be unstable while the condition results for SRI for wide gaps. Our previous finding that limits the SRI is only reproduced for medium-sized gaps. In all cases, however, the quasi-Kepler rotation proves to be unstable. New experiments with different gap sizes and larger Reynolds numbers could verify these results.
Figures 4 and 5 yield the basic results for SRIs subject to toroidal fields. The magnetic field is assumed to be as current-free in the fluid between the cylinders, i.e. excluding TI. Without any density stratification, no Rayleigh instability exists for quasi-Kepler rotation, but nonaxisymmetric modes with m=1 are unstable as a result of the interaction between differential rotation and magnetic field (``AMRI'').
The magnetic influence strongly depends on the magnetic Prandtl number. The instability needs a higher Reynolds number for , and it needs a lower Reynolds number for . After our experiences with MHD instabilities, this is not a surprise. It is a surprise, however, that the magnetic influence is always weak. Up to Hartmann numbers of , only a magnetic-induced factor of two plays a role^{}. Hence, our conclusion is that the SRI survives for rather high magnetic fields. For large Pm it is even supported by the toroidal magnetic field.
The combination of rotation, density stratification, and magnetic field leads to complex results. However, one basic finding is that the critical Reynolds numbers, if they exist, above which the flow becomes unstable without a magnetic field, are increased by the stable density stratification. In contrast, the critical Hartmann numbers, if any exist with the TI, above which the field becomes unstable without a rotation, do not depend on the stratification. Different transitions are possible between the two limits.
Generally speaking, the density stratification `balloons' the stability region, and in this sense it stabilizes the flow. For slow rotation the maximal stable magnetic field exceeds the critical magnetic field without rotation while the maximal stable magnetic field is smaller for faster rotation than this critical value. Even rather low values of the Reynolds numbers lead to a stabilization of those fields that are unstable for . The effect is strong for but becomes weaker for decreasing magnetic Prandtl numbers.
For steep radial profiles of the magnetic field (i.e. strong axial currents) and magnetic Prandtl numbers , one also finds the magnetic field destabilizing the Rayleigh instability; i.e., the critical Reynolds numbers with magnetic field are lower than the Reynolds numbers without magnetic field. This magnetic destabilization only exists for magnetic Prandtl numbers that are not too small. It exists for both uniform and density-stratified fluids (Fig. 10).
Finally, the -component of the electromotive force representing the -effect of the mean-field dynamo theory has been computed. The computations require an extreme degree of accuracy. The results demonstrate the importance of the density stratification for the existence of the -effect. Without density stratification, the correlations are vanishing in the radial average. With an axial density stratification, the calculations model the polar region of a rotating sphere or disk. With this interpretation accepted, we found the at the northern pole as positive for negative density gradients (Fig. 11).