No Arabic abstract
The $Omega$-phase of the liquid sodium $alpha$-$Omega$ dynamo experiment at NMIMT in cooperation with LANL has successfully demonstrated the production of a high toroidal field, $B_{phi} simeq 8times B_r$ from the radial component of an applied poloidal magnetic field, $B_r$. This enhanced toroidal field is produced by rotational shear in stable Couette flow within liquid sodium at $Rm simeq 120$. The small turbulence in stable Taylor-Couette flow is caused by Ekman flow where $ (delta v/v)^2 sim 10^{-3} $. This high $Omega$-gain in low turbulence flow contrasts with a smaller $Omega$-gain in higher turbulence, Helmholtz-unstable shear flows. This result supports the ansatz that large scale astrophysical magnetic fields are created within semi-coherent large scale motions in which turbulence plays only a smaller diffusive role that enables magnetic flux linkage.
In order to explore the magnetostrophic regime expected for planetary cores, experiments have been conducted in a rotating sphere filled with liquid sodium, with an imposed dipolar magnetic field (the DTS setup). The field is produced by a permanent magnet enclosed in an inner sphere, which can rotate at a separate rate, producing a spherical Couette flow. The flow properties are investigated by measuring electric potentials on the outer sphere, the induced magnetic field in the laboratory frame, and velocity profiles inside the liquid sodium using ultrasonic Doppler velocimetry. The present article focuses on the time-averaged axisymmetric part of the flow. The Doppler profiles show that the angular velocity of the fluid is relatively uniform in most of the fluid shell, but rises near the inner sphere, revealing the presence of a magnetic wind, and gently drops towards the outer sphere. The transition from a magnetostrophic flow near the inner sphere to a geostrophic flow near the outer sphere is controlled by the local Elsasser number. For Rossby numbers up to order 1, the observed velocity profiles all show a similar shape. Numerical simulations in the linear regime are computed, and synthetic velocity profiles are compared with the measured ones. In the geostrophic region, a torque-balance model provides very good predictions. We find that the induced magnetic field varies in a consistent fashion, and displays a peculiar peak in the counter-rotating regime. This happens when the fluid rotation rate is almost equal and opposite to the outer sphere rotation rate. The fluid is then almost at rest in the laboratory frame, and the Proudman-Taylor constraint vanishes, enabling a strong meridional flow. We suggest that dynamo action might be favored in such a situation.
A linearly unstable, sinusoidal $E times B$ shear flow is examined in the gyrokinetic framework in both the linear and nonlinear regimes. In the linear regime, it is shown that the eigenmode spectrum is nearly identical to hydrodynamic shear flows, with a conjugate stable mode found at every unstable wavenumber. In the nonlinear regime, turbulent saturation of the instability is examined with and without the inclusion of a driving term that prevents nonlinear flattening of the mean flow, and a scale-independent radiative damping term that suppresses the excitation of conjugate stable modes. A simple fluid model for how momentum transport and partial flattening of the mean flow scale with the driving term is constructed, from which it is shown that, except at high radiative damping, stable modes play an important role in the turbulent state and yield significantly improved quantitative predictions when compared with corresponding models neglecting stable modes.
The Kelvin-Helmholtz (KH) instability of a shear layer with an initially-uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with KH-unstable flows in hydrodynamics, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the magnetic field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited and thus transport less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear due to stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.
Magnetic fields pervade the entire Universe and affect the formation and evolution of astrophysical systems from cosmological to planetary scales. The generation and dynamical amplification of extragalactic magnetic fields through cosmic times, up to $mu$Gauss levels reported in nearby galaxy clusters, near equipartition with kinetic energy of plasma motions and on scales of at least tens of kiloparsecs, is a major puzzle largely unconstrained by observations. A dynamo effect converting kinetic flow energy into magnetic energy is often invoked in that context, however extragalactic plasmas are weakly collisional (as opposed to magnetohydrodynamic fluids), and whether magnetic-field growth and sustainment through an efficient turbulent dynamo instability is possible in such plasmas is not established. Fully kinetic numerical simulations of the Vlasov equation in a six-dimensional phase space necessary to answer this question have until recently remained beyond computational capabilities. Here, we show by means of such simulations that magnetic-field amplification via a dynamo instability does occur in a stochastically-driven, non-relativistic subsonic flow of initially unmagnetized collisionless plasma. We also find that the dynamo self-accelerates and becomes entangled with kinetic instabilities as magnetization increases. The results suggest that such a plasma dynamo may be realizable in laboratory experiments, support the idea that intracluster medium (ICM) turbulence may have significantly contributed to the amplification of cluster magnetic fields up to near-equipartition levels on a timescale shorter than the Hubble time, and emphasize the crucial role of multiscale kinetic physics in high-energy astrophysical plasmas.
Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport EMF averaging schemes that can meet this challenge, and using this strategy, we have developed a new Adaptive Mesh Refinement (AMR) MHD module of the ORION2 code. We investigate the effects of AMR on several statistical properties of a turbulent ideal MHD system with a thermal Mach number of 10 and a plasma $beta_0$ of 0.1 as initial conditions; our code is shown to be stable for simulations with higher Mach numbers ($M_rms = 17.3$) and smaller plasma beta ($beta_0 = 0.0067$) as well. Our results show that the quality of the turbulence simulation is generally related to the volume-averaged refinement. Our AMR simulations show that the turbulent dissipation coefficient for supersonic MHD turbulence is about 0.5, in agreement with unigrid simulations.