Large Magnetic Prandtl Numbers Research Articles

Magnetic fields generated by hydromagnetic dynamos at the low Prandtl number in dependence on the Ekman and magnetic Prandtl numbers

This article investigates the dependence of hydromagnetic dynamos on the magnetic Prandtl number at low Prandtl number. In all the investigated cases, the generated magnetic fields are dipolar and neither transition to hemispherical dynamos nor weaker magnetic fields (which are less dipole dominated) were observed, although the inertia becomes important. The magnetic field becomes weak in the polar regions (is “convected out of polar regions”) only for low Prandtl numbers, when the inertia becomes important. It is a basic condition. However, whether the magnetic field gets weak in the polar regions (is “convected out of polar regions”) or not depends also on the magnetic Prandtl number. The magnetic Prandtl number has to exceed a minimum value in order to sustain dynamo action. If the magnetic diffusion is small (large magnetic Prandtl numbers) then this phenomenon does not exist but if it is large (small magnetic Prandtl numbers) it exists because the strong magnetic diffusion significantly weakens the magnetic field inside the tangent cylinder. The magnetic diffusion and inertia seem to act in the same direction as to weaken the magnetic field inside the tangent cylinder.

Plasmoid and Kelvin-Helmholtz instabilities in Sweet-Parker current sheets

A two-dimensional (2D) linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of reduced magnetohydrodynamics. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al., Phys. Plasmas 14, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids (k(max)L(CS)~S(3/8), where k(max) is the wave number of fastest growing mode, S=L(CS)V(A)/η is the Lundquist number, L(CS) is the length of the sheet, V(A) is the Alfvén speed, and η is the plasma resistivity), which grows super Alfvénically fast (γ(max)τ(A)~S(1/4), where γ(max) is the maximum growth rate, and τ(A)=L(CS)/V(A)). For typical background profiles, the growth rate and the wave number are found to increase in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfvén speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability γ(max)τ(A)~k(max)L(CS)~S(1/2). The effect of viscosity (ν) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers Pm=ν/η, it is found that γ(max)~S(1/4)Pm(-5/8) and k(max)L(CS)~S(3/8)Pm(-3/16), leading to the prediction that the critical Lundquist number for plasmoid instability in the Pm>>1 regime is S(crit)~10(4)Pm(1/2). These results are verified via direct numerical simulation of the linearized equations, using an analytical 2D SP equilibrium solution.

Dissipation in dynamos at low and high magnetic Prandtl numbers

AbstractUsing simulations of helically driven turbulence, it is shown that the ratio of kinetic to magnetic energy dissipation scales with the magnetic Prandtl number in power law fashion with an exponent of approximately 0.6. Over six orders of magnitude in the magnetic Prandtl number the magnetic field is found tobe sustained by large‐scale dynamo action of alphasquared type. This work extends a similar finding for small magnetic Prandtl numbers to the regime of large magnetic Prandtl numbers. At large magnetic Prandtl numbers, most of the energy is dissipated viscously, lowering thus the amount of magnetic energy dissipation, which means that simulations can be performed at magnetic Reynolds numbers that are large compared to the usual limits imposed by a given resolution. This is analogous to an earlier finding that at small magnetic Prandtl numbers, most of the energy is dissipated resistively, lowering the amount of kinetic energy dissipation, so simulations can then be performed at much larger fluid Reynolds numbers than otherwise. The decrease in magnetic energy dissipation at large magnetic Prandtl numbers is discussed in the context of underluminous accretion found in some quasars (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

The Tayler instability of toroidal magnetic fields in a columnar gallium experiment

AbstractThe nonaxisymmetric Tayler instability of toroidal magnetic fields due to axial electric currents is studied for conducting incompressible fluids between two coaxial cylinders without endplates. The inner cylinder is considered as so thin that the limit of Rin → 0 can be computed. The magnetic Prandtl number is varied over many orders of magnitudes but the azimuthal mode number of the perturbations is fixed to m = 1. In the linear approximation the critical magnetic field amplitudes and the growth rates of the instability are determined for both resting and rotating cylinders. Without rotation the critical Hartmann numbers do not depend on the magnetic Prandtl number but this is not true for the corresponding growth rates. For given product of viscosity and magnetic diffusivity the growth rates for small and large magnetic Prandtl number are much smaller than those for Pm = 1. For gallium under the influence of a magnetic field at the outer cylinder of 1 kG the resulting growth time is 5 s. The minimum electric current through a container of 10 cm diameter to excite the instability is 3.20 kA. For a rotating container both the critical magnetic field and the related growth times are larger than for the resting column (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Magnetic Dynamo action at Low Magnetic Prandtl Numbers

Amplification of magnetic field due to kinematic turbulent dynamo action is studied in the regime of small magnetic Prandtl numbers. Such a regime is relevant for planets and stars interiors, as well as for liquid-metal laboratory experiments. A comprehensive analysis based on the Kazantsev-Kraichnan model is reported, which establishes the dynamo threshold and the dynamo growth rates for varying kinetic helicity of turbulent fluctuations. It is proposed that in contrast with the case of large magnetic Prandtl numbers, the kinematic dynamo action at small magnetic Prandtl numbers is significantly affected by kinetic helicity, and it can be made quite efficient with an appropriate choice of the helicity spectrum.

Hall-magnetohydrodynamic small-scale dynamos

Magnetic field generation by dynamo action is often studied within the theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the Hall effect may become non-negligible. We present results from three-dimensional simulations of the Hall-MHD equations subjected to random nonhelical forcing. We study the role of the Hall effect in the dynamo efficiency for different values of the Hall parameter. For small values of the Hall parameter, the small-scale dynamo is more efficient, displaying faster growth and saturating at larger amplitudes of the magnetic field. For larger values of the Hall parameter, saturation of the magnetic field is reached at smaller amplitudes than in the MHD case. We also study energy transfer rates among spatial scales and show that the Hall effect produces a reduction of the direct energy cascade at scales larger than the Hall scale, therefore leading to smaller energy dissipation rates. Finally, we present results stemming from simulations at large magnetic Prandtl numbers, which is the relevant regime in the hot and diffuse interstellar medium. In the range of magnetic Prandtl numbers considered, the Hall effect moves the peak of the magnetic energy spectrum as well as other relevant magnetic length scales toward the Hall scale.

Islands in the stream: The effect of plasma flow on tearing stability

Reducing plasma flow clearly decreases the stability of tearing modes in multiple regimes (sawtooth, hybrid) in both high- and low-aspect-ratio tokamaks (DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)], Joint European Torus [M. Keilhacker and the JET Team, Plasma Phys. Controlled Fusion 41, 301 (1999)], National Spherical Torus Experiment [M. Ono, S. M. Kaye, Y.-K. M. Peng, Nucl. Fusion 40, 557 (2000)], each with distinct means of lessening rotation). Further, reducing flow makes pre-existing “saturated” islands larger at the same beta (β). Thus lower plasma flow impairs high-beta operation owing both to the destabilization and to the impact of tearing-mode islands. Experimental results suggest that flow shear (not flow) at the tearing rational surface is classically stabilizing, making the effective tearing stability index Δ′ of the total current density profile more negative (more stable). In this picture, with profiles and all else the same, the minimum metastable beta at which neoclassical tearing modes (NTMs) can destabilize is proportional to −Δ′ and hence lower flow and flow shear lead to possible destabilization (depending on seeding) at lower beta. Similarly, if destabilized, the saturated NTM island width is proportional to −β/Δ′ and thus increases as flow and flow shear are reduced. A working model gives a significant level of stabilizing shear if the plasma toroidal angular flow shear −dΩϕ/dr at a given rational surface is of order of the inverse of the product of the local values of the parallel magnetic shear length Ls and the Alfvén time τA. Experimental data are fitted for the effect of this normalization of flow shear in a simple empirical model for both onset and saturation of tearing modes. Most theoretical literature is on the consequence of flow shear on tearing stability at zero beta; tokamaks at high beta have large magnetic Prandtl number (an issue for the sign of the flow effect) and very large Lundquist number. It is in this regime that theory will be compared with experimentally based empirical models. The consequence for future tokamaks with low rotation may be lower tearing stability than now expected.

MAGNETIC DYNAMO ACTION IN ASTROPHYSICAL TURBULENCE

We investigate the structure of magnetic field amplified by turbulent velocity fluctuations, in the framework of the kinematic Kazantsev-Kraichnan model. We consider Kolmogorov distribution of velocity fluctuations, and assume that both Reynolds number and magnetic Reynolds number are very large. We present the full numerical solution of the model for the spectra and the growth rates of magnetic fluctuations. We consider astrophysically relevant limits of large and small magnetic Prandtl numbers, and address both helical and nonhelical cases.

Bifurcations and dynamo action in a Taylor–Green flow

We report successive bifurcations in direct numerical simulations (DNSs) of a Taylor–Green flow, in both a hydro- and a magneto-hydrodynamic case. Hydrodynamic bifurcations occur in between different metastable states with different dynamo action, and are triggered by the numerical noise. The various states encountered range from stationary to chaotic or turbulent through possible oscillatory states. The corresponding sequence of bifurcations is reminiscent of the sequence obtained in the von Karman (VK) flow, at aspect ratio Γ=2 (Nore et al 2003 J. Fluid Mech. 477 51). We then use kinematic simulations to compute the dynamo thresholds of the different metastable states. A more detailed study of the turbulent state reveals the existence of two windows of dynamo action. Stochastic numerical simulations are then used to mimic the influence of turbulence on the dynamo threshold of the turbulent state. We show that the dynamo threshold is increased (respectively decreased) by the presence of large scale (resp. small scale) turbulent velocity fluctuations. Finally, DNSs of the magneto-hydrodynamic equations are used to explore the linear and nonlinear stage of the dynamo instability. In the linear stage, we show that the magnetic field favours the bifurcation from the basic state directly towards the turbulent or chaotic stable state. The magnetic field can also temporarily stabilize a metastable state, resulting in cycles of dynamo action, with different Lyapunov exponents. The critical magnetic Reynolds number for dynamo action is found to increase strongly with the Reynolds number. Finally, we provide a preliminary study of the saturation regime above the dynamo threshold. At large magnetic Prandtl number, we have observed two main types of saturations, in agreement with an analytical prediction of Leprovost and Dubrulle (2005 Eur. Phys. J. B 44 395): (i) intermittent dynamo, with vanishing most probable value of the magnetic energy; (ii) dynamo with non vanishing mean value of the magnetic energy. We describe a sequence of bifurcation of the dynamo at low Reynolds number and at increasing magnetic Prandtl number, where the dynamo switches from stationary mode to chaotic modes in a complex manner, involving intermittency in a way reminiscent to what is observed in dynamical systems with low number of degrees of freedom.

The effects of different parameter regimes in geodynamo simulations

Over the past 10 years, geodynamo simulations have grown rapidly in sophistication. However, it is still necessary to make certain approximations in order to maintain numerical stability. In addition, models are forced to make assumptions about poorly known parameters for the Earth's core. Different magnetic Prandtl numbers have been used and different assumptions about the presence of radiogenic heating have been made. This study examines some of the consequences of different approximations and assumptions using the Glatzmaier–Roberts geodynamo model. Here, we show that the choice of magnetic Prandtl number has a greater influence on the character of the magnetic field produced than the addition of a plausible amount of radiogenic heating. In particular, we find that prescribing a magnetic Prandtl number of unity with Ekman number limited by current computing resources, results in magnetic fields with significantly smaller intensities and variabilities compared with the much more Earth-like results obtained from simulations with large magnetic Prandtl numbers. A magnetic Prandtl number of unity, with both the viscous and magnetic diffusivities set to the Earth's magnetic diffusivity, requires a rotation rate much smaller than that of the Earth for currently reachable Ekman numbers. This results in a reduced dominance of the Coriolis forces relative to the buoyancy forces, and therefore, a reduction in the magnetic field intensity and the variability compared to the large Prandtl number cases.

Dynamics of the Magnetoviscous Instability

In dilute astrophysical plasmas, the collisional mean free path of a particle can exceed its Larmor radius. Under these conditions, the thermal conductivity and viscosity of the plasma can be dramatically altered. This alteration allows outwardly decreasing angular velocity or temperature gradients to become strongly destabilizing. This paper generalizes an earlier, simple analysis of the viscosity instability by including the dynamical effects of magnetic field line tension. Such effects lower the growth rates found in the absence of such tension but still allow growth rates in excess of the maximum of the standard magnetorotational instability. We find very good quantitative agreement with more complex kinetic treatments of the same process. The combination of large growth rates and large magnetic Prandtl number suggests that protogalactic disks are powerful dynamos.

Simulations of the Small‐Scale Turbulent Dynamo

We report the results of an extensive numerical study of the small-scale turbulent dynamo. The primary focus is on the case of large magnetic Prandtl numbers Prm, which is relevant for hot low-density astrophysical plasmas. A Prm parameter scan is given for the model case of viscosity-dominated (low Reynolds number) turbulence. We concentrate on three topics: magnetic energy spectra and saturation levels, the structure of the magnetic field lines, and intermittency of the field strength distribution. The main results are as follows: (1) the folded structure of the field (direction reversals at the resistive scale, field lines curved at the scale of the flow) persists from the kinematic to the nonlinear regime; (2) the field distribution is self-similar and appears to be lognormal during the kinematic regime and exponential in the saturated state; and (3) the bulk of the magnetic energy is at the resistive scale in the kinematic regime and remains there after saturation, although the magnetic energy spectrum becomes much shallower. We propose an analytical model of saturation based on the idea of partial two-dimensionalization of the velocity gradients with respect to the local direction of the magnetic folds. The model-predicted saturated spectra are in excellent agreement with numerical results. Comparisons with large-Re, moderate-Prm runs are carried out to confirm the relevance of these results and to test heuristic scenarios of dynamo saturation. New features at large Re are elongation of the folds in the nonlinear regime from the viscous scale to the box scale and the presence of an intermediate nonlinear stage of slower than exponential magnetic energy growth accompanied by an increase of the resistive scale and partial suppression of the kinetic energy spectrum in the inertial range. Numerical results for the saturated state do not support scale-by-scale equipartition between magnetic and kinetic energies, with a definite excess of magnetic energy at small scales. A physical picture of the saturated state is proposed. Subject heading gs: magnetic fields — methods: numerical — MHD — plasmas — turbulence

Seed fields for galactic dynamos by the magnetorotational instability

A linear but global numerical model for the magnetorotational instability (MRI) in disk geometry is considered to estimate the instability parameters for galaxies. Similarity rules suggested by a local analysis are applied to reveal a universal behavior approached by the results of global calculations for large magnetic Prandtl numbers. The findings are used to estimate the MRI characteristics for galaxies with their very large magnetic Prandtl numbers which cannot be attained in any numerical simulations. The resulting minimum field for the instability, Bmin 10-25 G, is small compared to any seed fields currently discussed. The growth times of MRI are estimated to be on the order of the rotation period of the inner rigidly-rotating core, i.e. ~100 Myr. Global MRI excites preferentially the magnetic field modes of quadrupolar symmetry.

Saturated state of the nonlinear small-scale dynamo.

We consider the problem of incompressible, forced, nonhelical, homogeneous, and isotropic MHD turbulence with no mean magnetic field and large magnetic Prandtl number. This type of MHD turbulence is the end state of the turbulent dynamo, which generates folded fields with small-scale direction reversals. We propose a model in which saturation is achieved as a result of the velocity statistics becoming anisotropic with respect to the local direction of the magnetic folds. The model combines the effects of weakened stretching and quasi-two-dimensional mixing and produces magnetic-energy spectra in remarkable agreement with numerical results at least in the case of a one-scale flow. We conjecture that the statistics seen in numerical simulations could be explained as a superposition of these folded fields and Alfvén-like waves that propagate along the folds.

A numerical study of forced magnetic reconnection in the viscous Taylor problem

Two-dimensional, nonlinear magnetohydrodynamical simulations are used to investigate the so-called Taylor problem, in which a small amplitude boundary perturbation is suddenly applied to a tearing stable, slab plasma equilibrium—the perturbation being such as to drive magnetic reconnection within the plasma. This type of reconnection, which is not due to an intrinsic plasma instability, is generally known as “forced reconnection.” For numerical reasons, the investigation is restricted to the large magnetic Prandtl number limit. The simulation results are highly consistent with the analysis of Hahm and Kulsrud [Phys. Fluids 28, 2412 (1985)] (modified by strong plasma viscosity). At high perturbation amplitudes, the system exhibits a phase of Sweet–Parker reconnection, as predicted by Wang and Bhattacharjee [Phys. Fluids B 4, 1795 (1992)]. An expression for the threshold perturbation amplitude required to trigger Sweet–Parker reconnection is derived, and successfully benchmarked against numerical simulations. This expression suggests that a Sweet–Parker phase is only likely to occur during forced reconnection in tokamaks when the plasma is extremely hot and perturbation amplitude relatively large.

Differential rotation and meridional flow in the solar supergranulation layer: Measuring the eddy viscosity

We measure the eddy viscosity in the outermost layers of the solar convection zone by comparing the rotation law computed with the Reynolds stress resulting from f-plane simulations of the angular momentum transport in rotating convection with the observed differential rotation pattern. The simulations lead to a negative vertical and a positive horizontal angular momentum transport. The consequence is a subrotation of the outermost layers, as it is indeed indicated both by helioseismology and the observed rotation rates of sunspots. In order to reproduce the observed gradient of the rotation rate a value of about 1.5 x 10^{13} cm/s for the eddy viscosity is necessary. Comparison with the magnetic eddy diffusivity derived from the sunspot decay yields a surprisingly large magnetic Prandtl number of 150 for the supergranulation layer. The negative gradient of the rotation rate also drives a surface meridional flow towards the poles, in agreement with the results from Doppler measurements. The successful reproduction of the abnormally positive horizontal cross correlation (on the northern hemisphere) observed for bipolar groups then provides an independent test for the resulting eddy viscosity.

Imbalanced Weak Magnetohydrodynamic Turbulence

Weak MHD turbulence consists of waves that propagate along magnetic field lines, in both directions. When two oppositely directed waves collide, they distort each other, without changing their respective energies. Each wave suffers many collisions before cascading; by contrast, in strong MHD turbulence, waves cascade on the same timescale at which they collide. Imbalance means that more energy is going in one direction than the other. In general, MHD turbulence is imbalanced. Yet imbalanced MHD cascades are not understood. For example, turbulence in the solar wind is observed to be imbalanced, so solar wind turbulence will not be understood until a theory of the imbalanced cascade is developed. We solve MHD turbulence that is imbalanced. Of crucial importance is that the energies going in both directions are forced to equalize at the dissipation scale. This pinning of the energy spectra was discovered by Grappin and coworkers. It affects the entire inertial range. Weak MHD turbulence is particularly interesting because perturbation theory is applicable. Hence, it can be described with a simple kinetic equation. Galtier and coworkers derived this kinetic equation. We present a simpler, more physical derivation, based on the picture of colliding wavepackets. In the process, we clarify the role of the zero-frequency mode. We also explain why Goldreich & Sridhar claimed that perturbation theory is inapplicable, and why this claim is wrong. (Our weak is equivalent to Goldreich & Sridhar's intermediate.) We perform numerical simulations of the kinetic equation to verify our claims. We construct simplified model equations that illustrate the main effects. Finally, we show that a large magnetic Prandtl number does not have a significant effect, and that hyperviscosity leads to a pronounced bottleneck effect.

A model of nonlinear evolution and saturation of the turbulent MHD dynamo

The growth and saturation of magnetic field in conductingturbulent media with large magnetic Prandtl numbers areinvestigated. This regime is very common in low-density hotastrophysical plasmas. During the early (kinematic) stage, weakmagnetic fluctuations grow exponentially and concentrate at theresistive scale, which lies far below the hydrodynamic viscousscale. The evolution becomes nonlinear when the magnetic energyis comparable to the kinetic energy of the viscous-scale eddies. A physical picture of the ensuing nonlinear evolution of the MHD dynamo is proposed. Phenomenological considerations are supplemented with a simple Fokker-Planck model of thenonlinear evolution of the magnetic-energy spectrum. It isfound that, while the shift of the bulk of the magnetic energy from the subviscous scales to the velocity scales may bepossible, it occurs very slowly - at the resistive, ratherthan dynamical, timescale (for galaxies, this means thatthe generation of large-scale magnetic fields cannot be explainedby this mechanism). The role of Alfvénic motions and theimplications for the fully developed isotropic MHD turbulenceare discussed.

The Small‐Scale Structure of Magnetohydrodynamic Turbulence with Large Magnetic Prandtl Numbers

We study the intermittency and field-line structure of the MHD turbulence in plasmas with very large magnetic Prandtl numbers. In this regime, which is realized in the interstellar medium, some accretion disks, protogalaxies, galaxy-cluster gas, the early universe, etc., magnetic fluctuations can be excited at scales below the viscous cutoff. The salient feature of the resulting small-scale magnetic turbulence is the folded structure of the fields. It is characterized by very rapid transverse spatial oscillation of the field direction, while the field lines remain largely unbent up to the scale of the flow. Quantitatively, the fluctuation level and the field-line geometry can be studied in terms of the statistics of the field strength and of the field-line curvature. In the kinematic limit, the distribution of the field strength is an expanding lognormal, while that of the field-line curvature K is stationary and has a power tail ~K-13/7. The field strength and curvature are anticorrelated, i.e., the growing fields are mostly flat, while the sharply curved fields remain relatively weak. The field, therefore, settles into a reduced-tension state. Numerical simulations demonstrate three essential features of the nonlinear regime. First, the total magnetic energy is equal to the total kinetic energy. Second, the intermittency is partially suppressed compared to the kinematic case, as the fields become more volume-filling and their distribution develops an exponential tail. Third, the folding structure of the field is unchanged from the kinematic case: the anticorrelation between the field strength and the curvature persists, and the distribution of the latter retains the same power tail. We propose a model of back-reaction based on the folding picture that reproduces all of the above numerical results.

Magnetic instability in a sheared azimuthal flow

We study the magneto-rotational instability of an incompressible flow which rotates with angular velocity (r )= a + b=r 2 where r is the radius and a and b are constants. We nd that an applied magnetic eld destabilises the flow, in agreement with the results of Rudiger & Zhang (2001). We extend the investigation in the region of parameter space which is Rayleigh stable. We also study the instability at values of magnetic Prandtl number which are much larger and smaller than Rudiger & Zhang. Large magnetic Prandtl numbers are motivated by their possible relevance in the central region of galaxies (Kulsrud & Anderson 1992). In this regime we nd that increasing the magnetic Prandtl number greatly enhances the instability; the stability boundary drops below the Rayleigh line and tends toward the solid body rotation line. Very small magnetic Prandtl numbers are motivated by the current MHD dynamo experiments performed using liquid sodium and gallium. Our nding in this regime conrms Rudiger & Zhang's conjecture that the linear magneto-rotational instability and the nonlinear hydrodynamical instability (Richard & Zahn 1999) take place at Reynolds numbers of the same order of magnitude.