Magnetohydrodynamic (MHD) turbulence is a cross-field process relevant to many systems. A prerequisite for understanding these systems is to constrain the role of MHD turbulence, and in particular, the energy exchange between kinetic and magnetic forms. The energetics of strongly magnetized and compressible turbulence has so far resisted attempts to understand them. Numerical simulations reveal that kinetic energy can be orders of magnitude higher than fluctuating magnetic energy. We solved this lack-of-balance puzzle by calculating the energetics of compressible and sub-Alfvénic turbulence based on the dynamics of coherent cylindrical fluid parcels. Using the MHD Lagrangian, we proved analytically that the bulk of the magnetic energy transferred to kinetic energy is the energy that is stored in the coupling between the ordered and fluctuating magnetic field. The analytical relations are in strikingly good agreement with numerical data, up to second-order terms.
more »
« less
On the properties and implications of collapse-driven MHD turbulence
ABSTRACT We investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, isothermal, magnetically supercritical molecular cloud core with transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along the gravitational acceleration vector, a proxy for the infall component of the velocity field. We find that (1) In spite of being supercritical, the core first contracts to a sheet perpendicular to the mean magnetic field, and the sheet itself collapses. (2) The solenoidal component of the turbulence remains at roughly its initial level throughout the simulation, while the compressible component increases continuously, implying that turbulence does not dissipate towards the centre of the core. (3) The distribution of simulation cells in the B–ρ plane occupies a wide triangular region at low densities, bounded below by the expected trend for fast MHD waves (B ∝ ρ, applicable for high-local Alfvénic Mach number MA) and above by the trend expected for slow waves (B ∼ constant, applicable for low local MA). At high densities, the distribution follows a single trend $$B \propto \rho ^{\gamma _{\rm eff}}$$, with 1/2 < γeff < 2/3, as expected for gravitational compression. (4) The mass-to-magnetic flux ratio λ increases with radius r due to the different scalings of the mass and magnetic flux with r. At a fixed radius, λ increases with time due to the accretion of material along field lines. (5) The solenoidal energy fraction is much smaller than the total turbulent component, indicating that the collapse drives the turbulence mainly compressibly, even in directions orthogonal to that of the collapse.
more »
« less
- Award ID(s):
- 2307840
- PAR ID:
- 10504396
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 530
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3431-3444
- Size(s):
- p. 3431-3444
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The transport of waves and turbulence beyond the photosphere is central to the coronal heating problem. Turbulence in the quiet solar corona has been modeled on the basis of the nearly incompressible magnetohydrodynamic (NI MHD) theory to describe the transport of low-frequency turbulence in open magnetic field regions. It describes the evolution of the coupled majority quasi-2D and minority slab component, driven by the magnetic carpet and advected by a subsonic, sub-Alfvénic flow from the lower corona. In this paper, we couple the NI MHD turbulence transport model with an MHD model of the solar corona to study the heating problem in a coronal loop. In a realistic benchmark coronal loop problem, we find that a loop can be heated to ∼1.5 million K by transport and dissipation of MHD turbulence described by the NI MHD model. We also find that the majority 2D component is as important as the minority slab component in the heating of the coronal loop. We compare our coupled MHD/NI MHD model results with a reduced MHD (RMHD) model. An important distinction between these models is that RMHD solves for small-scale velocity and magnetic field fluctuations and obtains the actual viscous/resistive dissipation associated with their evolution whereas NI MHD evolves scalar moments of the fluctuating velocity and magnetic fields and approximates dissipation using an MHD turbulence phenomenology. Despite the basic differences between the models, their simulation results match remarkably well, yielding almost identical heating rates inside the corona.more » « less
-
ABSTRACT Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical frameworks, we propose a novel Convolutional Neural Network (CNN) model to extract this embedded information, enabling the probe of 3D magnetic fields. This model examines the plane-of-the-sky magnetic field orientation (ϕ), the magnetic field’s inclination angle (γ) relative to the line-of-sight, and the total magnetization level (M$$_{\rm A}^{-1}$$) of the cloud. We train the model using synthetic emission lines of 13CO (J = 1–0) and C18O (J = 1–0), generated from 3D MHD simulations that span conditions from sub-Alfvénic to super-Alfvénic molecular clouds. Our tests confirm that the CNN model effectively reconstructs the 3D magnetic field topology and magnetization. The median uncertainties are under 5° for both ϕ and γ, and less than 0.2 for MA in sub-Alfvénic conditions (MA ≈ 0.5). In super-Alfvénic scenarios (MA ≈ 2.0), they are under 15° for ϕ and γ, and 1.5 for MA. We applied this trained CNN model to the L1478 molecular cloud. Results show a strong agreement between the CNN-predicted magnetic field orientation and that derived from Planck 353 GHz polarization. The CNN approach enabled us to construct the 3D magnetic field map for L1478, revealing a global inclination angle of ≈76° and a global MA of ≈1.07.more » « less
-
Abstract During its 10th orbit around the Sun, the Parker Solar Probe sampled two intervals where the local Alfvén speed exceeded the solar wind speed, lasting more than 10 hours in total. In this paper, we analyze the turbulence and wave properties during these periods. The turbulence is observed to be Alfvénic and unbalanced, dominated by outward-propagating modes. The power spectrum of the outward-propagating Elsässer z + mode steepens at high frequencies while that of the inward-propagating z − mode flattens. The observed Elsässer spectra can be explained by the nearly incompressible (NI) MHD turbulence model with both 2D and Alfvénic components. The modeling results show that the z + spectra are dominated by the NI/slab component, and the 2D component mainly affects the z − spectra at low frequencies. An MHD wave decomposition based on an isothermal closure suggests that outward-propagating Alfvén and fast magnetosonic wave modes are prevalent in the two sub-Alfvénic intervals, while the slow magnetosonic modes dominate the super-Alfvénic interval in between. The slow modes occur where the wavevector is nearly perpendicular to the local mean magnetic field, corresponding to nonpropagating pressure-balanced structures. The alternating forward and backward slow modes may also be features of magnetic reconnection in the near-Sun heliospheric current sheet.more » « less
-
Abstract We study anisotropic magnetohydrodynamic (MHD) turbulence in the slow solar wind measured by Parker Solar Probe (PSP) and Solar Orbiter (SolO) during its first orbit from the perspective of variance anisotropy and correlation anisotropy. We use the Belcher & Davis approach (M1) and a new method (M2) that decomposes a fluctuating vector into parallel and perpendicular fluctuating vectors. M1 and M2 calculate the transverse and parallel turbulence components relative to the mean magnetic field direction. The parallel turbulence component is regarded as compressible turbulence, and the transverse turbulence component as incompressible turbulence, which can be either Alfvénic or 2D. The transverse turbulence energy is calculated from M1 and M2, and the transverse correlation length from M2. We obtain the 2D and slab turbulence energy and the corresponding correlation lengths from those transverse turbulence components that satisfy an angle between the mean solar wind flow speed and mean magnetic field θ UB of either (i) 65° < θ UB < 115° or (ii) 0° < θ UB < 25° (155° < θ UB < 180°), respectively. We find that the 2D turbulence component is not typically observed by PSP near perihelion, but the 2D component dominates turbulence in the inner heliosphere. We compare the detailed theoretical results of a nearly incompressible MHD turbulence transport model with the observed results of PSP and SolO measurements, finding good agreement between them.more » « less