Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETOHYDRODYNAMIC TURBULENCE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

MHD Modeling of the Solar Wind with Turbulence Transport and Heating

We have developed a magnetohydrodynamic model that describes the global axisymmetric steady-state structure of the solar wind near solar minimum with account for transport of small-scale turbulence associated heating. The Reynolds-averaged mass, momentum, induction, and energy equations for the large-scale solar wind flow are solved simultaneously with the turbulence transport equations in the region from 0.3 to 100 AU. The large-scale equations include subgrid-scale terms due to turbulence and the turbulence (small-scale) equations describe the effects of transport and (phenomenologically) dissipation of the MHD turbulence based on a few statistical parameters (turbulence energy, normalized cross-helicity, and correlation scale). The coupled set of equations is integrated numerically for a source dipole field on the Sun by a time-relaxation method in the corotating frame of reference. We present results on the plasma, magnetic field, and turbulence distributions throughout the heliosphere and on the role of the turbulence in the large-scale structure and temperature distribution in the solar wind.

Goldstein, M. L.↗

The Transport of Density Fluctuations Throughout the Heliosphere

The solar wind is recognized as a turbulent magnetofluid, for which the properties of the turbulent velocity and magnetic field fluctuations are often described by the equations of incompressible magnetohydrodynamics (MHD). However, low-frequency density turbulence is also ubiquitous. On the basis of a nearly incompressible formulation of MHD in the expanding inhomogeneous solar wind, we derive the transport equation for the variance of the density fluctuations (Rho(exp 2)). The transport equation shows that density fluctuations behave as a passive scalar in the supersonic solar wind. In the absence of sources of density turbulence, such as within 1AU, the variance (Rho(exp 2)) approximates r(exp -4). In the outer heliosphere beyond 1 AU, the shear between fast and slow streams, the propagation of shocks, and the creation of interstellar pickup ions all act as sources of density turbulence. The model density fluctuation variance evolves with heliocentric distance within approximately 300 AU as (Rho(exp 2)) approximates r(exp -3.3) after which it flattens and then slowly increases. This is precisely the radial profile for the density fluctuation variance observed by Voyager 2. Using a different analysis technique, we confirm the radial profile for Rho(exp 2) of Bellamy, Cairns, & Smith using Voyager 2 data. We conclude that a passive scalar description for density fluctuations in the supersonic solar wind can explain the density fluctuation variance observed in both the inner and the outer heliosphere.

Zank, G. P.↗

Statistical Theory of the Ideal MHD Geodynamo

A statistical theory of geodynamo action is developed, using a mathematical model of the geodynamo as a rotating outer core containing an ideal (i.e., no dissipation), incompressible, turbulent, convecting magnetofluid. On the concentric inner and outer spherical bounding surfaces the normal components of the velocity, magnetic field, vorticity and electric current are zero, as is the temperature fluctuation. This allows the use of a set of Galerkin expansion functions that are common to both velocity and magnetic field, as well as vorticity, current and the temperature fluctuation. The resulting dynamical system, based on the Boussinesq form of the magnetohydrodynamic (MHD) equations, represents MHD turbulence in a spherical domain. These basic equations (minus the temperature equation) and boundary conditions have been used previously in numerical simulations of forced, decaying MHD turbulence inside a sphere [1,2]. Here, the ideal case is studied through statistical analysis and leads to a prediction that an ideal coherent structure will be found in the form of a large-scale quasistationary magnetic field that results from broken ergodicity, an effect that has been previously studied both analytically and numerically for homogeneous MHD turbulence [3,4]. The axial dipole component becomes prominent when there is a relatively large magnetic helicity (proportional to the global correlation of magnetic vector potential and magnetic field) and a stationary, nonzero cross helicity (proportional to the global correlation of velocity and magnetic field). The expected angle of the dipole moment vector with respect to the rotation axis is found to decrease to a minimum as the average cross helicity increases for a fixed value of magnetic helicity and then to increase again when average cross helicity approaches its maximum possible value. Only a relatively small value of cross helicity is needed to produce a dipole moment vector that is aligned at approx.10deg with the rotation axis.

Shebalin, J. V.↗

Turbulence in Core-collapse Supernovae

It is understood in a general sense that turbulent fluid motion below the shock front in a core-collapse supernova stiffens the effective equation of state of the fluid and aids in the revival of the explosion. However, when one wishes to be precise and quantify the amount of turbulence in a supernova simulation, one immediately encounters the problem that turbulence is difficult to define and measure. Using the 3D magnetohydrodynamic code ELEPHANT, we study how different definitions of turbulence change one’s conclusions about the amount of turbulence in a supernova and the extent to which it helps the explosion. We find that, while all the definitions of turbulence we use lead to a qualitatively similar growth pattern over time of the turbulent kinetic energy in the gain region, the total amount of turbulent kinetic energy, and especially the ratios of turbulent to total kinetic energy, distinguish them. Some of the definitions appear to indicate turbulence is a necessary contributor to the explosion, and others indicate it is not. The different definitions also produce turbulence maps with different correlations with maps of the enstrophy, a quantity widely regarded as also indicating the presence of turbulence. We also compute the turbulent adiabatic index and observe that, in regions of low enstrophy, this quantity is sensitive to the definition used. As a consequence, the effective adiabatic index depends upon the method used to measure the turbulence, and thus it alters one’s conclusions regarding the impact of turbulence within the supernova.

computational astronomy↗

Fading Coronal Structure and the Onset of Turbulence in the Young Solar Wind

Above the top of the solar corona, the young, slow solar wind transitions from low-beta, magnetically structured flow dominated by radial structures to high-beta, less structured flow dominated by hydrodynamics. This transition, long inferred via theory, is readily apparent in the sky region close to 10deg from the Sun in processed, background-subtracted solar wind images. We present image sequences collected by the inner Heliospheric Imager instrument on board the Solar-Terrestrial Relations Observatory (STEREO/HI1) in 2008 December, covering apparent distances from approximately 4deg to 24deg from the center of the Sun and spanning this transition in the large-scale morphology of the wind. We describe the observation and novel techniques to extract evolving image structure from the images, and we use those data and techniques to present and quantify the clear textural shift in the apparent structure of the corona and solar wind in this altitude range. We demonstrate that the change in apparent texture is due both to anomalous fading of the radial striae that characterize the corona and to anomalous relative brightening of locally dense puffs of solar wind that we term "flocculae." We show that these phenomena are inconsistent with smooth radial flow, but consistent with the onset of hydrodynamic or magnetohydrodynamic instabilities leading to a turbulent cascade in the young solar wind.

DeForest, C. E.↗

Renormalization group analysis of reduced magnetohydrodynamics with application to subgrid modeling

The technique for obtaining a subgrid model for Navier-Stokes turbulence, based on renormalization group analysis (RNG), is extended to the reduced magnetohydrodynamic (RMND) equations. It is shown that a RNG treatment of the Alfven turbulence supported by the RMHD equations leads to effective values of the viscosity and resistivity at large scales, k yields 0, dependent on the amplitude of turbulence. The effective viscosity and resistivity become independent of the molecular quantities when the RNG analysis is augmented by the Kolmogorov argument for energy cascade. A self-contained system of equations is derived for the range of scales, k = 0-K, where K = pi/Delta is the maximum wave number for a grid size Delta. Differential operators, whose coefficients depend upon the amplitudes of the large-scale quantities, represent in this system the resistive and viscous dissipation.

Longcope, D. W.↗

Cascade and Dissipation of Solar Wind Turbulence at Electron Scales: Whistlers or Kinetic Alfv\'en Waves?

Over the past few decades, large-scales solar wind (SW) turbulence has been studied extensively, both theoretically and observationally. Observed power spectra of the low frequency turbulence, which can be described in the magnetohydrodynamic (MHD) limit, are shown to obey the Kolmogorov scaling, $k"{ -5/3 }$, down the local proton gyrofrequency ($C{ci} \sim O.l$-Hz). Turbulence at frequencies above $C{ci}$ has not been thoroughly investigated and remains far less well understood. Above $C{ ci}$ the spectrum steepens to $\sim f"{ -2.5}$ and a debate exists as to whether the turbulence has become dominated by dispersive kinetic Alfven waves (KA W) or by whistler waves, before it is dissipated at small scales, In a case study Sahraoui et al., PRL (2009) have reported the first direct determination of the dissipation range of solar wind turbulence near the electron gyroscale using the high resolution Cluster magnetic and electric field data (up to $10"2$-Hz in the spacecraft reference frame). Above the Doppler-shifted proton scale $C{\rho i}$ a new inertial range with a scaling $\sim f"{ -2.3}$ has been evidenced and shown to remarkably agree with theoretical predictions of a quasi-two-dimensional cascade into KA W turbulence. Here, we use a wider sample of data sets of small scale SW turbulence under different plasma conditions, and investigate under which physical criteria the KA W (or the whistler) turbulence may be observed to carry out the cascade at small scales, These new observations/criteria are compared to the predictions on the cascade and the (kinetic) dissipation from the Vlasov theory. Implications of the results on the heating problem of the solar wind will be discussed.

Sahraoui, Fouad↗

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Here, we use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number S limit (above ~104), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate ~0.01–0.13 is weakly dependent on S. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a superdiffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗

Layered patterns of active scalar fields in a two-dimensional magnetohydrodynamic system

Here, we observe the formation of staircase patterns in the magnetic potential (𝐴) in a weakly magnetized two-dimensional magnetohydrodynamic system driven by a forced, fluctuating vortex array. Layering occurs due to inhomogeneous mixing of 𝐴 by vortex cells. Magnetic Reynolds number (𝑅 𝑚 )–dependent quenching of the turbulent diffusion of 𝐴 by weak magnetic fields increases the disparity between the (short) cell circulation time and the (long) time for intercell transport of magnetic potential. Thus, magnetic fields strengthen transport barriers between cells and reinforce the staircase, relative to its passive scalar counterpart. The analysis reveals a feedback mechanism, which promotes staircase formation. Magnetic staircases persist in both the flux expulsion (𝑅 𝑚 ⁢𝑣$^2_𝐴$/𝑈$^2_0$<1) and vortex disruption (𝑅 𝑚 ⁢𝑣$^2_𝐴$/𝑈$^2_0$≥1) limits. In the latter case, residual vortex cells homogenize 𝐴. Global layering morphology is shown to be well characterized by staircase curvature. Stochastic forcing of magnetic potential can support magnetic staircases against resistive decay.

magnetohydrodynamics↗

Turbulence link to L-mode, I-mode, and H-mode confinement in the DIII-D tokamak

Understanding the physics of low-confinement (L-), improved-confinement (I-), and high-confinement (H-) modes is critical for fusion reactors. The finding herein reports observations of two types of turbulence coexisting near the L-mode edge, one magnetohydrodynamic (MHD)-like and another micro-tearing mode (MTM)-like, linked to the H-mode and I-mode confinement in the DIII-D tokamak. Ion-scale magnetic and density turbulence is measured using a Faraday-effect radial-interferometer-polarimeter and beam-emission-spectroscopy (BES). Broadband turbulence spectra of up to ~600 kHz are observed in two discharges where transitions between L-mode, I-mode, and H-mode occurs. Turbulence is found to be inversely correlated with confinement, meaning lower turbulence power at higher confinement. Distinctively, the high-frequency (HF, >~100 kHz) magnetic turbulence power changes by the most (55%) during transitions primarily involving energy confinement change, whereas the low-frequency (LF, <~100 kHz) magnetic and density turbulence power changes by the most (80%) during transitions primarily involving particle confinement change. The LF turbulence amplitude oscillates with and leads to deuterium-alpha emission oscillations before an H-mode. These results imply that HF turbulence mainly affects energy confinement whereas LF turbulence can affect particle confinement. The magnetic and density turbulence exhibits coherence up to 0.6 and cross-phase magnitude close to π/2 in most cases, suggesting they have a common origin in both the LF and HF ranges. BES suggests that LF turbulence resides at the edge (ρ = 0.95) and HF turbulence can be at the outer core (ρ = 0.8) or edge (ρ = 0.95). Comparisons of measurements, theory, and gyrokinetic simulations suggest that HF turbulence is MTM-like in all cases, whereas LF turbulence is more consistent with MHD-like modes and the exact instability might change during transitions—except that a drift-wave origin is possible in a low collisionality H-mode. These results suggest that the H-mode involves suppressed MHD-like turbulence, whereas the I-mode mitigates MTM-like turbulence along with largely unchanged MHD-like turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetofluid Turbulence in the Solar Wind

The solar wind shows striking characteristics that suggest that it is a turbulent magnetofluid, but the picture is not altogether simple. From the earliest observations, a strong correlation between magnetic fluctuations and plasma velocity fluctuations was noted. The high corrections suggest that the fluctuations are Alfven waves. In addition, the power spectrum of the magnetic fluctuation showed evidence of an inertial range that resembled that seen in fully-developed fluid turbulence. Alfven waves, however, are exact solutions of the equations of incompressible magnetohydrodynamics. Thus, there was a puzzle: how can a magnetofluid consisting of Alfven waves be turbulent? The answer lay in the role of velocity shears in the solar wind that could drive turbulent evolution. Puzzles remain: for example, the power spectrum of the velocity fluctuations is less steep than the slope of the magnetic fluctuations, nor do we understand even now why the solar wind appears to be nearly incompressible with a -5/3 power-spectral index.

Goldstein, Melvyn L.↗

Turbulence and Global Properties of the Solar Wind

The solar wind shows striking characteristics that suggest that it is a turbulent magnetofluid, but the picture is not altogether simple. From the earliest observations, a strong correlation between magnetic fluctuations and plasma velocity fluctuations was noted. The high corrections suggest that the fluctuations are Alfven waves. In addition, the power spectrum of the magnetic fluctuation showed evidence of an inertial range that resembled that seen in fully-developed fluid turbulence. Alfven waves, however, are exact solutions of the equations of incompressible magnetohydrodynamics. Thus, there was a puzzle: how can a magnetofluid consisting of Alfven waves be turbulent? The answer lay in the role of velocity shears in the solar wind that could drive turbulent evolution. Puzzles remain: for example, the power spectrum of the velocity fluctuations is less steep than the slope of the magnetic fluctuations. The plasma in the magnetic tail of Earth's magnetosphere also shows aspects of turbulence, as does the plasma in the dayside magnetosphere near the poles the dayside cusps. Recently, new analyses of high time resolution magnetic field data from Cluster have offered a glimpse of how turbulence is dissipated, thus heating the ambient plasma.

Goldstein, Melvyn L.↗

Shell Models of RMHD Turbulence and the Heating of Solar Coronal Loops

A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B|| and stratification, nicknamed 'Shell-Atm,' is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfve'n waves propagating along B||. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and haracteristic properties of heating events are distributed as power laws. Crosscorrelations show a delay of dissipation compared with energy content.

corona↗