No Arabic abstract
Via amplification by turbulent dynamo, magnetic fields can be potentially important for the formation of the first stars. To examine the dynamo behavior during the gravitational collapse of primordial gas, we extend the theory of nonlinear turbulent dynamo to include the effect of gravitational compression. The relative importance between dynamo and compression varies during contraction, with the transition from dynamo- to compression-dominated amplification of magnetic fields with the increase of density. In the nonlinear stage of magnetic field amplification with the scale-by-scale energy equipartition between turbulence and magnetic fields, reconnection diffusion of magnetic fields in ideal magnetohydrodynamic (MHD) turbulence becomes important. It causes the violation of flux-freezing condition and accounts for (a) the small growth rate of nonlinear dynamo, (b) the weak dependence of magnetic energy on density during contraction, (c) the saturated magnetic energy, and (d) the large correlation length of magnetic fields. The resulting magnetic field structure and the scaling of magnetic field strength with density are radically different from the expectations of flux-freezing.
The excitation and further sustenance of large-scale magnetic fields in rotating astrophysical systems, including planets, stars and galaxies, is generally thought to involve a fluid magnetic dynamo effect driven by helical magnetohydrodynamic turbulence. While this scenario is appealing on general grounds, it however currently remains largely unconstrained, notably because a fundamental understanding of the nonlinear asymptotic behaviour of large-scale fluid magnetism in the astrophysically-relevant but treacherous regime of large magnetic Reynolds number $Rm$ is still lacking. We explore this problem using local high-resolution simulations of turbulent magnetohydrodynamics driven by an inhomogeneous helical forcing generating a sinusoidal profile of kinetic helicity, mimicking the hemispheric distribution of kinetic helicity in rotating turbulent fluid bodies. We identify a transition at large $Rm$ to an asymptotic nonlinear state, followed up to $Rmsimeq 3times 10^3$, characterized by an asymptotically small resistive dissipation of magnetic helicity, by its efficient spatial redistribution across the equator through turbulent fluxes driven by the hemispheric distribution of kinetic helicity, and by the presence in the tangled dynamical magnetic field of plasmoids typical of reconnection at large $Rm$.
The growth and saturation of magnetic field in conducting turbulent media with large magnetic Prandtl numbers are investigated. This regime is very common in low-density hot astrophysical plasmas. During the early (kinematic) stage, weak magnetic fluctuations grow exponentially and concentrate at the resistive scale, which lies far below the hydrodynamic viscous scale. The evolution becomes nonlinear when the magnetic energy is 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 the nonlinear evolution of the magnetic-energy spectrum. It is found that, while the shift of the bulk of the magnetic energy from the subviscous scales to the velocity scales may be possible, it occurs very slowly -- at the resistive, rather than dynamical, time scale (for galaxies, this means that generation of large-scale magnetic fields cannot be explained by this mechanism). The role of Alfvenic motions and the implications for the fully developed isotropic MHD turbulence are discussed.
We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.
Small-scale turbulent dynamo is responsible for the amplification of magnetic fields on scales smaller than the driving scale of turbulence in diverse astrophysical media. Most earlier dynamo theories concern the kinematic regime and small-scale magnetic field amplification. Here we review our recent progress in developing the theories for the nonlinear dynamo and the dynamo regime in a partially ionized plasma. The importance of reconnection diffusion of magnetic fields is identified for both the nonlinear dynamo and magnetic field amplification during gravitational contraction. For the dynamo in a partially ionized plasma, the coupling state between neutrals and ions and the ion-neutral collisional damping can significantly affect the dynamo behavior and the resulting magnetic field structure. We present both our analytical predictions and numerical tests with a two-fluid dynamo simulation on the dynamo features in this regime. In addition, to illustrate the astrophysical implications, we discuss several examples for the applications of the dynamo theory to studying magnetic field evolution in both preshock and postshock regions of supernova remnants, in weakly magnetized molecular clouds, during the (primordial) star formation, and during the first galaxy formation.
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $ abla cdot {bf B}=0$ to machine precision, albeit at significant computational expense.