RESUMEN
Alfvénic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons (REs), and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compressional Alfvén eigenmodes (CAEs) in DIII-D disruption scenarios, providing an explanation of the CQ modes. Simulation results reveal that high energy trapped REs can have resonance with the Alfvén mode through their toroidal precession motion, and the resonance frequency is proportional to the energy of REs. The mode frequencies and their relationship with the RE energy are consistent with experimental observations. The perturbed magnetic fields from the modes can lead to spatial diffusion of REs including the nonresonant passing ones, thus providing the theoretical basis for a potential approach for RE mitigation.
RESUMEN
Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a previously undiscovered range of energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain a spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, for example, solar coronal heating, can be further explored with current and future spacecraft and telescopes.
RESUMEN
As noted in Wikipedia, skin in the game refers to having 'incurred risk by being involved in achieving a goal', where 'skin is a synecdoche for the person involved, and game is the metaphor for actions on the field of play under discussion'. For exascale applications under development in the US Department of Energy Exascale Computing Project, nothing could be more apt, with the skin being exascale applications and the game being delivering comprehensive science-based computational applications that effectively exploit exascale high-performance computing technologies to provide breakthrough modelling and simulation and data science solutions. These solutions will yield high-confidence insights and answers to the most critical problems and challenges for the USA in scientific discovery, national security, energy assurance, economic competitiveness and advanced healthcare. This article is part of a discussion meeting issue 'Numerical algorithms for high-performance computational science'.
RESUMEN
The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At a sufficiently large magnetic Reynolds number (R_{m}=10^{6}), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index -2.2 as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.
RESUMEN
The effects of kinetic whistler wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from experiments at the DIII-D National Fusion Facility, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results. We find an enhancement of the RE avalanche for low density and high electric field, but for high density and low electric field the scattering can suppress the avalanche and raise the threshold electric field, bringing the present model much closer to observations. The excitation of kinetic instabilities and the scattering of resonant electrons are calculated self-consistently using a quasilinear model and local approximation. We also explain the observed fast growth of electron cyclotron emission signals and excitation of very low-frequency whistler modes observed in the quiescent RE experiments at DIII-D tokamak. Simulations using ITER parameters show that by controlling the background thermal plasma density and temperature, the plasma waves can also be excited spontaneously in tokamak disruptions and the avalanche generation of runaway electrons may be suppressed.
RESUMEN
Recent theory has demonstrated a novel physics regime for magnetic reconnection in high-energy-density plasmas where the magnetic field is advected by heat flux via the Nernst effect. Here we elucidate the physics of the electron dissipation layer in this regime. Through fully kinetic simulation and a generalized Ohm's law derived from first principles, we show that momentum transport due to a nonlocal effect, the heat-flux-viscosity, provides the dissipation mechanism for magnetic reconnection. Scaling analysis, and simulations show that the reconnection process comprises a magnetic field compression stage and quasisteady reconnection stage, and the characteristic width of the current sheet in this regime is several electron mean-free paths. These results show the important interplay between nonlocal transport effects and generation of anisotropic components to the distribution function.