The next generation dynamo experiment currently under development at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) will consist of a precessing cylindrical container filled with liquid sodium. We perform numerical simulations of kinematic dynamo action applying a velocity field that is obtained from hydrodynamic models of a precession driven flow. So far, the resulting magnetic field growth-rates remain below the dynamo threshold for magnetic Reynolds numbers up to Rm=2000.

The azimuthal version of the magnetorotational instability (MRI) is a non-axisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid metal Taylor-Couette experiment that shows the occurrence of an azimuthal MRI in the expected range of Hartmann numbers.

In the last decades, the experimental study of dynamo action has made great progress. However, after the dynamo experiments in Karlsruhe and Riga, the von-Karman-Sodium (VKS) dynamo is only the third facility that has been able to demonstrate fluid flow driven self-generation of magnetic fields in a laboratory experiment. Further progress in the experimental examination of dynamo action is expected from the planned precession driven dynamo experiment that will be designed in the framework of the liquid sodium facility DRESDYN (DREsden Sodium facility for DYNamo and thermohydraulic studies). In this paper, we briefly present numerical models of the VKS dynamo that demonstrate the close relation between the axisymmetric field observed in that experiment and the soft iron material used for the flow driving impellers. We further show recent results of preparatory water experiments and design studies related to the precession dynamo and delineate the scientific prospects for the final set-up.

Many astrophysical phenomena (such as the slow rotation of neutron stars or the rigid rotation of the solar core) can be explained by the action of the Tayler instability of toroidal magnetic fields in the radiative zones of stars. In order to place the theory of this instability on a safe fundament it has been realized in a laboratory experiment measuring the critical field strength, the growth rates as well as the shape of the supercritical modes. A strong electrical current flows through a liquid-metal confined in a resting columnar container with an insulating outer cylinder. As the very small magnetic Prandtl number of the gallium-indium-tin alloy does not influence the critical Hartmann number of the field amplitudes, the electric currents for marginal instability can also be computed with direct numerical simulations. The results of this theoretical concept are confirmed by the experiment. Also the predicted growth rates of the order of minutes for the nonaxisymmetric perturbations are certified by the measurements. That they do not directly depend on the size of the experiment is shown as a consequence of the weakness of the applied fields and the absence of rotation.

The eigenvalues and eigenfunctions of a linear {\alpha}^{2}-dynamo have been computed for different spatial distributions of an isotropic \alpha-effect. Oscillatory solutions are obtained when \alpha exhibits a sign change in the radial direction. The time-dependent solutions arise at so called exceptional points where two stationary modes merge and continue as an oscillatory eigenfunction with conjugate complex eigenvalues. The close proximity of oscillatory and non-oscillatory solutions may serve as the basic ingredient for reversal models that describe abrupt polarity switches of a dipole induced by noise. Whereas the presence of an inner core with different magnetic diffusivity has remarkable little impact on the character of the dominating dynamo eigenmodes, the introduction of equatorial symmetry breaking considerably changes the geometric character of the solutions. Around the dynamo threshold the leading modes correspond to hemispherical dynamos even when the symmetry breaking is small. This behavior can be explained by the approximate dipole-quadrupole degeneration for the unperturbed problem. More complicated scenarios may occur in case of more realistic anisotropies of \alpha- and \beta-effect or through non-linearities caused by the back-reaction of the magnetic field (magnetic quenching).

In the current-driven, kink-type Tayler instability (TI) a sufficiently strong azimuthal magnetic field becomes unstable against non-axisymmetric perturbations. The TI has been discussed as a possible ingredient of the solar dynamo mechanism and a source of the helical structures in cosmic jets. It is also considered as a size limiting factor for liquid metal batteries. We report on a liquid metal TI experiment using a cylindrical column of the eutectic alloy GaInSn to which electrical currents of up to 8 kA are applied. We present results of external magnetic field measurements that indicate the occurrence of the TI in good agreement with numerical predictions. The interference of TI with the competing large scale convection, resulting from Joule heating, is also discussed.

Starting from the present version of the Riga dynamo experiment with its rotating magnetic eigenfield dominated by a single frequency we ask for those modifications of this set-up that would allow for a non-trivial magnetic field behaviour in the saturation regime. Assuming an increased ratio of azimuthal to axial flow velocity, we obtain energy oscillations with a frequency below the eigenfrequency of the magnetic field. These new oscillations are identified as magneto-inertial waves that result from a slight imbalance of Lorentz and inertial forces. Increasing the azimuthal velocity further, or increasing the total magnetic Reynolds number, we find transitions to a chaotic behaviour of the dynamo.

The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end-caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end-caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one.

We analyze numerically the magnetorotational instability of a Taylor-Couette flow in a helical magnetic field (HMRI) using the inductionless approximation defined by a zero magnetic Prandtl number (Pm=0). The Chebyshev collocation method is used to calculate the eigenvalue spectrum for small amplitude perturbations. First, we carry out a detailed conventional linear stability analysis with respect to perturbations in the form of Fourier modes that corresponds to the convective instability which is not in general self-sustained. The helical magnetic field is found to extend the instability to a relatively narrow range beyond its purely hydrodynamic limit defined by the Rayleigh line. There is not only a lower critical threshold at which HMRI appears but also an upper one at which it disappears again. The latter distinguishes the HMRI from a magnetically-modified Taylor vortex flow. Second, we find an absolute instability threshold as well. In the hydrodynamically unstable regime before the Rayleigh line, the threshold of absolute instability is just slightly above the convective one although the critical wave length of the former is noticeably shorter than that of the latter. Beyond the Rayleigh line the lower threshold of absolute instability rises significantly above the corresponding convective one while the upper one descends significantly below its convective counterpart. As a result, the extension of the absolute HMRI beyond the Rayleigh line is considerably shorter than that of the convective instability. The absolute HMRI is supposed to be self-sustained and, thus, experimentally observable without any external excitation in a system of sufficiently large axial extension.

The magnetorotational instability (MRI) plays an essential role in the formation of stars and black holes. By destabilizing hydrodynamically stable Keplerian flows, the MRI triggers turbulence and enables outward transport of angular momentum in accretion discs. We present the results of a liquid metal Taylor-Couette experiment under the influence of helical magnetic fields that show typical features of MRI at Reynolds numbers of the order 1000 and Hartmann numbers of the order 10. Particular focus is laid on an improved experiment in which split end caps are used to minimize the Ekman pumping.

It is widely known that cosmic magnetic fields, i.e. the fields of planets, stars, and galaxies, are produced by the hydromagnetic dynamo effect in moving electrically conducting fluids. It is less well known that cosmic magnetic fields play also an active role in cosmic structure formation by enabling outward transport of angular momentum in accretion disks via the magnetorotational instability (MRI). Considerable theoretical and computational progress has been made in understanding both processes. In addition to this, the last ten years have seen tremendous efforts in studying both effects in liquid metal experiments. In 1999, magnetic field self-excitation was observed in the large scale liquid sodium facilities in Riga and Karlsruhe. Recently, self-excitation was also obtained in the French "von Karman sodium" (VKS) experiment. An MRI-like mode was found on the background of a turbulent spherical Couette flow at the University of Maryland. Evidence for MRI as the first instability of an hydrodynamically stable flow was obtained in the "Potsdam Rossendorf Magnetic Instability Experiment" (PROMISE). In this review, the history of dynamo and MRI related experiments is delineated, and some directions of future work are discussed.

The asymmetric time dependence and various statistical properties of polarity reversals of the Earth's magnetic field are utilized to infer some of the most essential parameters of the geodynamo, among them the effective (turbulent) magnetic diffusivity, the degree of supercriticality, and the relative strength of the periodic forcing which is believed to result from the Milankovic cycle of the Earth's orbit eccentricity. A time-stepped spherically symmetric alpha^2-dynamo model is used as the kernel of an inverse problem solver in form of a downhill simplex method which converges to solutions that yield a stunning correspondence with paleomagnetic data.

On 11 November 1999, a self-exciting magnetic eigenfield was detected for the first time in the Riga liquid sodium dynamo experiment. We report on the long history leading to this event, and on the subsequent experimental campaigns which provided a wealth of data on the kinematic and the saturated regime of this dynamo. The present state of the theoretical understanding of both regimes is delineated, and some comparisons with other laboratory dynamo experiments are made.

Using a spherical symmetric mean field alpha^2-dynamo model for Earth's magnetic field reversals, we show the coexistence of the noise-induced phenomena coherence resonance and stochastic resonance. Stochastic resonance has been recently invoked to explain the 100 kyr periodicity in the distribution of the residence time between reversals. The comparison of the resulting residence time distribution with the paleomagnetic one allows for some estimate of the effective diffusion time of the Earth's core which may be 100 kyr or slightly below rather than 200 kyr as it would result from the molecular resistivity.

The conventional magnetic induction equation that governs hydromagnetic dynamo action is transformed into an equivalent integral equation system. An advantage of this approach is that the computational domain is restricted to the region occupied by the electrically conducting fluid and to its boundary. This integral equation approach is first employed to simulate kinematic dynamos excited by Beltrami-like flows in a finite cylinder. The impact of externally added layers around the cylinder on the onset of dynamo actions is investigated. Then it is applied to simulate dynamo experiments within cylindrical geometry including the von Karman sodium (VKS) experiment and the Riga dynamo experiment. A modified version of this approach is utilized to investigate magnetic induction effects under the influence of externally applied magnetic fields which is also important to measure the proximity of a given dynamo facility to the self-excitation threshold.

The asymmetric shape of reversals of the Earth's magnetic field indicates a possible connection with relaxation oscillations as they were early discussed by van der Pol. A simple mean-field dynamo model with a spherically symmetric $\alpha$ coefficient is analysed with view on this similarity, and a comparison of the time series and the phase space trajectories with those of paleomagnetic measurements is carried out. For highly supercritical dynamos a very good agreement with the data is achieved. Deviations of numerical reversal sequences from Poisson statistics are analysed and compared with paleomagnetic data. The role of the inner core is discussed in a spectral theoretical context and arguments and numerical evidence is compiled that the growth of the inner core might be important for the long term changes of the reversal rate and the occurrence of superchrons.

The magnetorotational instability (MRI) plays a key role in the formation of stars and black holes, by enabling outward angular momentum transport in accretion disks. The use of combined axial and azimuthal magnetic fields allows the investigation of this effect in liquid metal flows at moderate Reynolds and Hartmann numbers. A variety of experimental results is presented showing evidence for the occurrence of the MRI in a Taylor-Couette flow using the liquid metal alloy GaInSn.

We study experimentally the flow of a liquid metal confined between differentially rotating cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For increasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the cylinders. The wavelengths and speeds of these structures, as well as the field strengths and rotation rates at which they arise, are broadly consistent with theoretical predictions of such a traveling wave magnetorotational instability.

A recent paper [R. Hollerbach and G. Rudiger, Phys. Rev. Lett. 95, 124501 (2005)] has shown that the threshold for the onset of the magnetorotational instability (MRI) in a Taylor-Couette flow is dramatically reduced if both axial and azimuthal magnetic fields are imposed. In agreement with this prediction, we present results of a Taylor-Couette experiment with the liquid metal alloy GaInSn, showing evidence for the existence of the MRI at Reynolds numbers of order 1000 and Hartmann numbers of order 10.

We consider the magnetorotational instability (MRI) of a hydrodynamically stable Taylor-Couette flow with a helical external magnetic field in the inductionless approximation defined by a zero magnetic Prandtl number ($\Pm=0)$. This leads to a considerable simplification of the problem eventually containing only hydrodynamic variables. First, we point out that the energy of any perturbation growing in the presence of magnetic field has to grow faster without the field. This is a paradox because the base flow is stable without the magnetic while it is unstable in the presence of a helical magnetic field without being modified by the latter as it has been found recently by Hollerbach and Rudiger [Phys. Rev. Lett. 95, 124501 (2005)]. We revisit this problem by using a Chebyshev collocation method to calculate the eigenvalue spectrum of the linearized problem. In this way, we confirm that MRI with helical magnetic field indeed works in the inductionless limit where the destabilization effect appears as an effective shift of the Rayleigh line. Second, we integrate the linearized equations in time to study the transient behavior of small amplitude perturbations, thus showing that the energy arguments are correct as well. However, there is no real contradiction between both facts. The linear stability theory predicts the asymptotic development of an arbitrary small-amplitude perturbation, while the energy stability theory yields the instant growth rate of any particular perturbation, but it does not account for the evolution of this perturbation.

The intention of the ''von Karman sodium'' (VKS) experiment is to study the hydromagnetic dynamo effect in a highly turbulent and unconstrained flow. Much effort has been devoted to the optimization of the mean flow and the lateral boundary conditions in order to minimize the critical magnetic Reynolds number and hence the necessary motor power. The main focus of this paper lies on the role of ''lid layers'', i.e. layers of liquid sodium between the impellers and the end walls of the cylinder. First, we study an analytical test flow to show that lid layers can have an ambivalent effect on the efficiency of the dynamo. The critical magnetic Reynolds number shows a flat minimum for a small lid layer thickness, but increases for thicker layers. For the actual VKS geometry it is shown that static lid layers yield a moderate increase of the critical magnetic Reynolds number by approximately 12 per cent. A more dramatic increase by 100 until 150 per cent can occur when some rotational flow is taken into account in those layers. Possible solutions of this problem are discussed for the real dynamo facility.

The irregular polarity reversals of the Earth's magnetic field have attracted much interest during the last decades. Despite the fact that recent numerical simulations of the geodynamo have shown nice polarity transitions, the very reason and the basic mechanism of reversals are far from being understood. Using a paradigmatic mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha we attribute the essential features of reversals to the magnetic field dynamics in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. At such exceptional (branch) points of square root type two real eigenvalues coalesce and continue as a complex conjugated pair of eigenvalues. Special focus is laid on the comparison of numerically computed time series with paleomagnetic observations. It is shown that the considered dynamo model with high supercriticality can explain the observed time scale and the asymmetric shape of reversals with a slow decay and a fast field recovery.

In a recent paper (Phys. Rev. Lett. 94 (2005), 184506; physics/0411050) it was shown that a simple mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha can exhibit a number of features which are typical for Earth's magnetic field reversals. In particular, the model produces asymmetric reversals, a positive correlation of field strength and interval length, and a bimodal field distribution. All these features are attributable to the magnetic field dynamics in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. The negative slope of the growth rate curve between the nearby local maximum and the exceptional point makes the system unstable and drives it to the exceptional point and beyond into the oscillatory branch where the sign change happens. A weakness of this reversal model is the apparent necessity to fine-tune the magnetic Reynolds number and/or the radial profile of alpha. In the present paper, it is shown that this fine-tuning is not necessary in the case of higher supercriticality of the dynamo. Numerical examples and physical arguments are compiled to show that, with increasing magnetic Reynolds number, there is strong tendency for the exceptional point and the associated local maximum to move close to the zero growth rate line. Although exemplified again by the spherically symmetric alpha^2 dynamo model, the main idea of this ''self-tuning'' mechanism of saturated dynamos into a reversal-prone state seems well transferable to other dynamos. As a consequence, reversing dynamos might be much more typical and may occur much more frequently in nature than what could be expected from a purely kinematic perspective.

Using a mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha which is dynamically quenched and disturbed by additional noise, the basic features of geomagnetic polarity reversals are shown to be generic consequences of the dynamo action in the vicinity of exceptional points of the spectrum. This simple paradigmatic model yields long periods of constant polarity which are interrupted by self-accelerating field decays leading to asymmetric polarity reversals. It shows the recently discovered bimodal field distribution, and it gives a natural explanation of the correlation between polarity persistence time and field strength. In addition, we find typical features of coherence resonance in the dependence of the persistence time on the noise.

The spectral branching behavior of the 2x2 operator matrix of the magneto-hydrodynamic alpha^2-dynamo is analyzed numerically. Some qualitative aspects of level crossings are briefly discussed with the help of a simple toy model which is based on a Z_2-graded-pseudo-Hermitian 2x2 matrix. The considered issues comprise: the underlying SU(1,1) symmetry and the Krein space structure of the system, exceptional points of branching type and diabolic points, as well as the algebraic and geometric multiplicity of corresponding degenerate eigenvalues.

With only a few exceptions, the numerical simulation of cosmic and laboratory hydromagnetic dynamos has been carried out in the framework of the differential equation method. However, the integral equation method is known to provide robust and accurate tools for the numerical solution of many problems in other fields of physics. The paper is intended to facilitate the use of integral equation solvers in dynamo theory. In concrete, the integral equation method is employed to solve the eigenvalue problem for a hydromagnetic dynamo model with a spherically symmetric, isotropic helical turbulence parameter alpha. Three examples of the function alpha(r) with steady and oscillatory solutions are considered. A convergence rate proportional to the inverse squared of the number of grid points is achieved. Based on this method, a convergence accelerating strategy is developed and the convergence rate is improved remarkably. Typically, quite accurate results can be obtained with a few tens of grid points. In order to demonstrate its suitability for the treatment of dynamos in other than spherical domains, the method is also applied to alpha^2 dynamos in rectangular boxes. The magnetic fields and the electric potentials for the first eigenvalues are visualized.

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of Biot-Savart's law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos in finite domains. For the spherically symmetric alpha^2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and the toroidal field components. The integral equation formulation for spherical dynamos with general velocity fields is also derived. Two numerical examples, the alpha^2 dynamo model with radially varying alpha, and the Bullard-Gellman model illustrate the equivalence of the approach with the usual differential equation method.

Until recently, the existence of oscillatory mean-field dynamos of the alpha^2-type with spherically symmetric and isotropic alpha was an open question. We find such dynamos by means of an evolutionary strategy, and we illustrate the spectral properties of the corresponding dynamo operators.

After the dynamo experiment in November 1999 had shown magnetic field self-excitation in a spiraling liquid metal flow, in a second series of experiments emphasis was placed on the magnetic field saturation regime as the next principal step in the dynamo process. The dependence of the strength of the magnetic field on the rotation rate is studied. Various features of the saturated magnetic field are outlined and possible saturation mechanisms are discussed.

For a simple spherically symmetric mean-field dynamo model we investigate the possibility of determining the radial dependence of the coefficient $\alpha$. Growth rates for different magnetic field modes are assumed to be known by measurement. An evolutionary strategy (ES) is used for the solution of the inverse problem. Numerically, we find quite different $\alpha$-profiles giving nearly the same eigenvalues. The ES is also applied to find functions $\alpha(r)$ yielding zero growth rates for the lowest four magnetic field modes. Additionally, a slight modification of the ES is utilized for an ''energetic'' optimization of $\alpha^2$-dynamos. The consequences of our findings for inverse dynamo theory and for the design of future dynamo experiments are discussed.

The paper deals with the integral equation approach to steady kinematic dynamo models in finite domains based on Biot-Savart's law. The role of the electric potential at the boundary is worked out explicitly. As an example, a modified version of the simple spherical $\alpha$-effect dynamo model proposed by Krause and Steenbeck is considered in which the $\alpha$-coefficient is no longer constant but may vary with the radial coordinate. In particular, the results for the original model are re-derived. Possible applications of this integral equation approach for numerical simulations of dynamos in arbitrary geometry and for an ''inverse dynamo theory'' are sketched.

In an experiment at the Riga sodium dynamo facility, a slowly growing magnetic field eigenmode has been detected over a period of about 15 seconds. For a slightly decreased propeller rotation rate, additional measurements showed a slow decay of this mode. The measured results correspond satisfactory with numerical predictions for the growth rates and frequencies.