We examine the relationship between magnetic flux generation, taken as an indicator of large-scale dynamo action, and magnetic helicity, computed as an integral over the dynamo volume, in a simple dynamo. We consider dynamo action driven by Magneto-Rotational Turbulence (MRT) within the shearing-box approximation. We consider magnetically open boundary conditions that allow a flux of helicity in or out of the computational domain. We circumvent the problem of the lack of gauge invariance in open domains by choosing a particular gauge -- the winding gauge -- that provides a natural interpretation in terms of average winding number of pairwise field lines. We use this gauge precisely to define and measure the helicity and helicity flux for several realizations of dynamo action. We find in these cases, that the system as a whole does not break reflectional symmetry and the total helicity remains small even in cases when substantial magnetic flux is generated. We find no particular connection between the generation of magnetic flux and the helicity or the helicity flux through the boundaries. We suggest that this result may be due to the essentially nonlinear nature of the dynamo processes in MRT.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

The Earth's magnetic field is generated by dynamo action driven by convection in the outer core. For numerical reasons, inertial and viscous forces play an important role in geodynamo models; however, the primary dynamical balance in the Earth's core is believed to be between buoyancy, Coriolis and magnetic forces. The hope has been that by setting the Ekman number to be as small as computationally feasible, an asymptotic regime would be reached in which the correct force balance is achieved. However, recent analyses of geodynamo models suggest that the desired balance has still not yet been attained. Here we adopt a complementary approach consisting of a model of rapidly rotating convection in which inertial forces are neglected from the outset. Within this framework we are able to construct a new branch of solutions in which the dynamo generates a strong magnetic field that satisfies the expected force balance. The resulting strongly magnetized convection is dramatically different to the corresponding solutions in which the field is weak.

In this work we study the process of energy dissipation triggered by a slow large scale motion of a magnetized conducting fluid. Our consideration is motivated by the problem of heating the solar corona, which is believed to be governed by fast reconnection events set off by the slow motion of magnetic field lines anchored in the photospheric plasma. To elucidate the physics governing the disruption of the imposed laminar motion and the energy transfer to small scales, we propose a simplified model where the large-scale motion of magnetic field lines is prescribed not at the footpoints but rather imposed volumetrically. As a result, the problem can be treated numerically with an efficient, highly-accurate spectral method, allowing us to use a resolution and statistical ensemble exceeding those of the previous work. We find that, even though the large-scale deformations are slow, they eventually lead to reconnection events that drive a turbulent state at smaller scales. The small-scale turbulence displays many of the universal features of field-guided MHD turbulence like a well developed inertial range spectrum. Based on these observations, we construct a phenomenological model that gives the scalings of the amplitude of the fluctuations and the energy dissipation rate as functions of the input parameters. We find a good agreement between the numerical results and the predictions of the model.

We present an approach to deriving global properties of accretion disks from the knowledge of local solutions derived from numerical simulations based on the shearing box approximation. The approach consists of a two-step procedure. First a local solution valid for all values of the disk height is constructed by piecing together an interior solution obtained numerically with an analytical exterior radiative solution. The matching is obtained by assuming hydrostatic balance and radiative equilibrium. Although in principle the procedure can be carried out in general, it simplifies considerably when the interior solution is fully convective. In these cases, the construction is analogous to the derivation of the Hayashi tracks for protostars. The second step consists of piecing together the local solutions at different radii to obtain a global solution. Here we use the symmetry of the solutions with respect to the defining dimensionless numbers--in a way similar to the use of homology relations in stellar structure theory--to obtain the scaling properties of the various disk quantities with radius.

Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{\perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{\perp}^{-1.5}$ for smaller simulation domain, and $k_{\perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.

The inertial-interval energy spectrum of strong magnetohydrodynamic (MHD) turbulence with a uniform background magnetic field was observed numerically to be close to $k^{-3/2}$ by a number of independent groups. A dissenting opinion has been voiced by Beresnyak, A. 2011, PRL, 106, 075001-. 2012, MNRAS, 422, 3495-. 2014, ApJ, 784, L20 that the spectral scaling is close to $k^{-5/3}$. The conclusions of these papers are however incorrect as they are based on numerical simulations that are drastically unresolved, so that the discrete numerical scheme does not approximate the physical solution at the scales where the measurements are performed. These results have been rebutted in our more detailed papers Perez, J. C., Mason, J., Boldyrev, S., & Cattaneo, F. 2012, PRX, 2, 041005-. 2014, ApJL, 793, L13; here, by popular demand, we present a brief and simple explanation of our major criticism of Beresnyak's work.

A variety of studies of magnetised plasma turbulence invoke theories for the advection of a passive scalar by turbulent fluctuations. Examples include modelling the electron density fluctuations in the interstellar medium, understanding the chemical composition of galaxy clusters and the intergalactic medium, and testing the prevailing phenomenological theories of magnetohydrodynamic turbulence. While passive scalar turbulence has been extensively studied in the hydrodynamic case, its counterpart in MHD turbulence is significantly less well understood. Herein we conduct a series of high-resolution direct numerical simulations of incompressible, field-guided, MHD turbulence in order to establish the fundamental properties of passive scalar evolution. We study the scalar anisotropy, establish the scaling relation analogous to Yaglom's law, and measure the intermittency of the passive scalar statistics. We also assess to what extent the pseudo Alfven fluctuations in strong MHD turbulence can be modelled as a passive scalar. The results suggest that the dynamics of a passive scalar in MHD turbulence is considerably more complicated than in the hydrodynamic case.

Magnetohydrodynamic (MHD) turbulence in the majority of natural systems, including the interstellar medium, the solar corona, and the solar wind, has Reynolds numbers far exceeding the Reynolds numbers achievable in numerical experiments. Much attention is therefore drawn to the universal scaling properties of small-scale fluctuations, which can be reliably measured in the simulations and then extrapolated to astrophysical scales. However, in contrast with hydrodynamic turbulence, where the universal structure of the inertial and dissipation intervals is described by the Kolmogorov self-similarity, the scaling for MHD turbulence cannot be established based solely on dimensional arguments due to the presence of an intrinsic velocity scale -- the Alfven velocity. In this Letter, we demonstrate that the Kolmogorov first self-similarity hypothesis cannot be formulated for MHD turbulence in the same way it is formulated for the hydrodynamic case. Besides profound consequences for the analytical consideration, this also imposes stringent conditions on numerical studies of MHD turbulence. In contrast with the hydrodynamic case, the discretization scale in numerical simulations of MHD turbulence should decrease faster than the dissipation scale, in order for the simulations to remain resolved as the Reynolds number increases.

We consider the generation of magnetic activity --- dynamo waves --- in the astrophysical limit of very large magnetic Reynolds number. We consider kinematic dynamo action for a system consisting of helical flow and large-scale shear. We demonstrate that large-scale dynamo waves persist at high $Rm$ if the helical flow is characterised by a narrow band of spatial scales and the shear is large enough. However for a wide band of scales the dynamo becomes small-scale with a further increase of $Rm$, with dynamo waves re-emerging only if the shear is then increased. We show that at high $Rm$ the key effect of the shear is to suppress small-scale dynamo action, allowing large-scale dynamo action to be observed. We conjecture that this supports a general "suppression principle" --- large-scale dynamo action can only be observed if there is a mechanism that suppresses the small-scale fluctuations.

We consider the problem of convergence in stratified isothermal shearing boxes with zero net magnetic flux. We present results with the highest resolution to-date--up to 200 grid-point per pressure scale height--that show no clear evidence of convergence. Rather, the Maxwell stresses continue to decrease with increasing resolution. We propose some possible scenarios to explain the lack of convergence based on multi-layer dynamo systems.

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. We choose radiative boundary conditions at the vertical boundaries in which the heat flux is propor- tional to the fourth power of the temperature. We compare the results with the corresponding cases in which fixed temperature boundary conditions are applied. The most notable result is that the formation of a fully convective state in which the density is nearly constant as a function of height and the heat is transported to the upper and lower boundaries by overturning motions is robust and persists even in cases with radiative boundary conditions. Interestingly, in the convective regime, although the diffusive transport is negligible the mean stratification does not relax to an adiabatic state.

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero mean magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. The calculations begin from an isothermal state spanning three scale heights above and below the mid-plane. After a long transient the layers settle to a stationary state in which thermal losses out of the boundaries are balanced by dissipative heating. We identify two regimes. A conductive regime in which the heat is transported mostly by conduction and the density decreases with height. In the limit of large thermal diffusivity this regime resembles the more familiar isothermal case. Another, the convective regime, observed at smaller values of the thermal diffusivity, in which the layer becomes unstable to overturning motions, the heat is carried mostly by advection and the density becomes nearly constant throughout the layer. In this latter constant-density regime we observe evidence for large-scale dynamo action leading to a substantial increase in transport efficiency relative to the conductive cases.

The energy spectrum of magnetohydrodynamic turbulence attracts interest due to its fundamental importance and its relevance for interpreting astrophysical data. Here we present measurements of the energy spectra from a series of high-resolution direct numerical simulations of MHD turbulence with a strong guide field and for increasing Reynolds number. The presented simulations, with numerical resolutions up to 2048^3 mesh points and statistics accumulated over 30 to 150 eddy turnover times, constitute, to the best of our knowledge, the largest statistical sample of steady state MHD turbulence to date. We study both the balanced case, where the energies associated with Alfv\'en modes propagating in opposite directions along the guide field, E^+ and $E^-, are equal, and the imbalanced case where the energies are different. In the balanced case, we find that the energy spectrum converges to a power law with exponent -3/2 as the Reynolds number is increased, consistent with phenomenological models that include scale-dependent dynamic alignment. For the imbalanced case, with E^+>E^-, the simulations show that E^- ~ k_{\perp}^{-3/2} for all Reynolds numbers considered, while E^+ has a slightly steeper spectrum at small Re. As the Reynolds number increases, E^+ flattens. Since both E^+ and E^- are pinned at the dissipation scale and anchored at the driving scales, we postulate that at sufficiently high Re the spectra will become parallel in the inertial range and scale as E^+ ~ E^- ~ k_{\perp}^{-3/2}. Questions regarding the universality of the spectrum and the value of the "Kolmogorov constant" are discussed.

Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for the turbulent dynamics can be built. Numerical simulations of MHD turbulence are widely used to guide and test the theoretical predictions; however, simulating MHD turbulence and accurately measuring its scaling properties is far from straightforward. Computational power limits the calculations to moderate Reynolds numbers and often simplifying assumptions are made in order that a wider range of scales can be accessed. After describing the theoretical predictions and the numerical approaches that are often employed in studying strong incompressible MHD turbulence, we present the findings of a series of high-resolution direct numerical simulations. We discuss the effects that insufficiencies in the computational approach can have on the solution and its physical interpretation.

We consider the problem of convergence in homogeneous shearing box simula- tions of magneto-rotationally driven turbulence. When there is no mean magnetic flux, if the equations are non dimensionalized with respect to the diffusive scale, the only free parameter in the problem is the size of the computational domain. The problem of convergence then relates to the asymptotic form of the solutions as the computational box size becomes large. By using a numerical code with a high order of accuracy we show that the solutions become asymptotically inde- pendent of domain size. We also show that cases with weak magnetic flux join smoothly to the zero flux cases as the flux vanishes. These results are consistent with the operation of a subcritical small-scale dynamo driving the turbulence. We conclude that for this type of turbulence the angular momentum transport is a proportional to the diffusive flux and therefore has limited relevance in as- trophysical situations.

Magnetised turbulence is ubiquitous in astrophysical systems, where it notoriously spans a broad range of spatial scales. Phenomenological theories of MHD turbulence describe the self-similar dynamics of turbulent fluctuations in the inertial range of scales. Numerical simulations serve to guide and test these theories. However, the computational power that is currently available restricts the simulations to Reynolds numbers that are significantly smaller than those in astrophysical settings. In order to increase computational efficiency and, therefore, probe a larger range of scales, one often takes into account the fundamental anisotropy of field-guided MHD turbulence, with gradients being much slower in the field-parallel direction. The simulations are then optimised by employing the reduced MHD equations and relaxing the field-parallel numerical resolution. In this work we explore a different possibility. We propose that there exist certain quantities that are remarkably stable with respect to the Reynolds number. As an illustration, we study the alignment angle between the magnetic and velocity fluctuations in MHD turbulence, measured as the ratio of two specially constructed structure functions. We find that the scaling of this ratio can be extended surprisingly well into the regime of relatively low Reynolds number. However, the extended scaling becomes easily spoiled when the dissipation range in the simulations is under-resolved. Thus, taking the numerical optimisation methods too far can lead to spurious numerical effects and erroneous representation of the physics of MHD turbulence, which in turn can affect our ability to correctly identify the physical mechanisms that are operating astrophysical systems.

A review of MHD dynamos and turbulence.

Hydromagnetic dynamo theory provides the prevailing theoretical description for the origin of magnetic fields in the universe. Here we consider the problem of kinematic, small-scale dynamo action driven by a random, incompressible, non-helical, homogeneous and isotropic flow. In the Kazantsev dynamo model the statistics of the driving flow are assumed to be instantaneously correlated in time. Here we compare the results of the model with the dynamo properties of a simulated flow that has equivalent spatial characteristics as the Kazantsev flow but different temporal statistics. In particular, the simulated flow is a solution of the forced Navier-Stokes equations and hence has a finite correlation time. We find that the Kazantsev model typically predicts a larger magnetic growth rate and a magnetic spectrum that peaks at smaller scales. However, we show that by filtering the diffusivity spectrum at small scales it is possible to bring the growth rates into agreement and simultaneously align the magnetic spectra.

Numerical simulations are an important tool in furthering our understanding of turbulent dynamo action, a process that occurs in a vast range of astrophysical bodies. It is important in all computational work that comparisons are made between different codes and, if non-trivial differences arise, that these are explained. Kapyla et al (2010: MNRAS 402, 1458) describe an attempt to reproduce the results of Hughes & Proctor (2009: PRL 102, 044501) and, by employing a different methodology, they arrive at very different conclusions concerning the mean electromotive force and the generation of large-scale fields. Here we describe why the simulations of Kapyla et al (2010) are simply not suitable for a meaningful comparison, since they solve different equations, at different parameter values and with different boundary conditions. Furthermore we describe why the interpretation of Kapyla et al (2010) of the calculation of the alpha-effect is inappropriate and argue that the generation of large-scale magnetic fields by turbulent convection remains a problematic issue.

Plasma experiments in laboratory settings offer unique opportunities to address fundamental aspects of the solar dynamo and magnetism in the solar atmosphere. We argue here that ground-based laboratory experiments have direct connections to NASA based missions and NSF programs, and that a small investment in laboratory heliophysics may have a high payoff. We advocate for broad involvement in community-scale plasma experiments.

Magnetohydrodynamic (MHD) turbulence is pervasive in astrophysical systems. Recent high-resolution numerical simulations suggest that the energy spectrum of strong incompressible MHD turbulence is $E(k_{\perp})\propto k_{\perp}^{-3/2}$. So far, there has been no phenomenological theory that simultaneously explains this spectrum and satisfies the exact analytic relations for MHD turbulence due to Politano & Pouquet. Indeed, the Politano-Pouquet relations are often invoked to suggest that the spectrum of MHD turbulence instead has the Kolmogorov scaling -5/3. Using geometrical arguments and numerical tests, here we analyze this seeming contradiction and demonstrate that the -3/2 scaling and the Politano-Pouquet relations are reconciled by the phenomenon of scale-dependent dynamic alignment that was recently discovered in MHD turbulence.

We present numerical simulations of circular Couette flow in axisymmetric and fully three-dimensional geometry of a cylindrical annulus inspired by Princeton MRI liquid gallium experiment. The incompressible Navier-Stokes equations are solved with the spectral element code Nek5000 incorporating realistic horizontal boundary conditions of differentially rotating rings. We investigate the effect of changing rotation rates (Reynolds number) and of the horizontal boundary conditions on flow structure, Ekman circulation and associated transport of angular momentum through the onset of unsteadiness and three-dimensionality. A mechanism for the explanation of the dependence of the Ekman flows and circulation on horizontal boundary conditions is proposed.

We discuss the applicability of the kinematic $\alpha$-effect formalism at high magnetic Reynolds numbers. In this regime the underlying flow is likely to be a small-scale dynamo, leading to the exponential growth of fluctuations. Difficulties arise with both the actual calculation of the $\alpha$ coefficients and with its interpretation. We argue that although the former may be circumvented -- and we outline several procedures by which the the $\alpha$ coefficients can be computed in principle -- the interpretation of these quantities in terms of the evolution of the large-scale field may be fundamentally flawed.

Aims: We study the changes in the properties of turbulence driven by the magnetorotational instability in a shearing box, as the computational domain size in the radial direction is varied relative to the height Methods: We perform 3D simulations in the shearing box approximation, with a net magnetic flux, and we consider computational domains with different aspect ratios Results: We find that in boxes of aspect ratio unity the transport of angular momentum is strongly intermittent and dominated by channel solutions in agreement with previous work. In contrast, in boxes with larger aspect ratio, the channel solutions and the associated intermittent behavior disappear. Conclusions: There is strong evidence that, as the aspect ratio becomes larger, the characteristics of the solution become aspect ratio independent. We conclude that shearing box calculations with aspect ratio unity or near unity may introduce spurious effects.

We report the results of an extensive set of direct numerical simulations of forced, incompressible, magnetohydrodynamic turbulence with a strong guide field. The aim is to resolve the controversy regarding the power law exponent (\alpha, say) of the field perpendicular energy spectrum E(k_\perp) \propto k_\perp ^ {\alpha}. The two main theoretical predictions, \alpha=-3/2 and \alpha=-5/3, have both received some support from numerical simulations carried out by different groups. Our simulations have a resolution of 512^3 mesh points, a strong guide field, an anisotropic simulation domain, and implement a broad range of large-scale forcing routines, including those previously reported in the literature. Our findings indicate that the spectrum of well developed, strong incompressible MHD turbulence with a strong guide field is E(k_{\perp})\propto k_{\perp}^{-3/2}.

Motivated by recent analytic predictions, we report numerical evidence showing that in driven incompressible magnetohydrodynamic turbulence the magnetic- and velocity-field fluctuations locally tend to align the directions of their polarizations. This dynamic alignment is stronger at smaller scales with the angular mismatch between the polarizations decreasing with the scale \lambda approximately as \theta_\lambda ~ \lambda^{1/4}. This can naturally lead to a weakening of the nonlinear interactions and provide an explanation for the energy spectrum E(k) ~ k^{-3/2} that is observed in numerical experiments of strongly magnetized turbulence.

We investigate analytically the amplification of a weak magnetic field in a homogeneous and isotropic turbulent flow lacking reflectional symmetry (helical turbulence). We propose that the spectral distributions of magnetic energy and magnetic helicity can be found as eigenmodes of a self-adjoint, Schr\"odinger-type system of evolution equations. We argue that large-scale and small-scale magnetic fluctuations cannot be effectively separated, and that the conventional alpha-model is, in general, not an adequate description of the large-scale dynamo mechanism. As a consequence, the correct numerical modeling of such processes should resolve magnetic fluctuations down to the very small, resistive scales.

We analyze the initial, kinematic stage of magnetic field evolution in an isotropic and homogeneous turbulent conducting fluid with a rough velocity field, v(l) ~ l^alpha, alpha<1. We propose that in the limit of small magnetic Prandtl number, i.e. when ohmic resistivity is much larger than viscosity, the smaller the roughness exponent, alpha, the larger the magnetic Reynolds number that is needed to excite magnetic fluctuations. This implies that numerical or experimental investigations of magnetohydrodynamic turbulence with small Prandtl numbers need to achieve extremely high resolution in order to describe magnetic phenomena adequately.

We employ the magnetic and velocity fields from turbulent dynamo simulations to synthesize the polarization of a typical photospheric line. The synthetic Stokes profiles have properties in common with those observed in the quiet Sun. The simulated magnetograms present a level of signal similar to that of the Inter-Network regions. Asymmetric Stokes V profiles with two, three and more lobes appear in a natural way. The intensity profiles are broadened by the magnetic fields in fair agreement with observational limits. Furthermore, the Hanle depolarization signals of the Sr I 4607 A line turn out to be within the solar values. Differences between synthetic and observed polarized spectra can also be found. There is a shortage of Stokes V asymmetries, that we attribute to a deficit of structuring in the magnetic and velocity fields from the simulations as compared to the Sun This deficit may reflect the fact that the Reynolds numbers of the numerical data are still far from solar values. We consider the possibility that intense and tangled magnetic fields, like those in the simulations, exist in the Sun. This scenario has several important consequences. For example, less than 10% of the existing unsigned magnetic flux would be detected in present magnetograms. The existing flux would exceed by far that carried by active regions during the maximum of the solar cycle. Detecting these magnetic fields would involve improving the angular resolution, the techniques to interpret the polarization signals, and to a less extent, the polarimetric sensitivity.