SEARCH · Engineering Papers
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.
A new end-plugging method for rotating plasmas is identified and analyzed. It uses the ponderomotive potential associated with an azimuthal magnetostatic wiggler. Studied both analytically and numerically, this process compares favorably to other end-plugging methods in open field line magnetized plasma devices.
How much of the density pedestal structure profile is determined by fueling versus transport is still a valid question. In most current tokamaks, the Scrape-Off Layer is not opaque enough to screen the neutrals and limit their penetration past the separatrix. As such recycled neutrals as well as additional gas puffing can directly affect the pedestal structure through ionization. Yet, in ITER and other future burning plasma magnetic confinement devices, this ionization source will be pushed further out into the Scrape-Off Layer. This poses the question, how much of the pedestal density structure is governed by this edge ionization source versus plasma transport effects. Theoretically, several turbulent modes have been identified which could provide a 'pinch' like up-gradient transport mechanism. Up-gradient transport is necessary in a source-free region to obtain a peaked density profile, which is often observed in the core of tokamaks. In an opaque Scrape-Off Layer, without the existence of such a pinch, the pedestal density structure would eventually disappear. Not only would this limit our ability to fuel a fusion reactor, it would also affect the pedestal stability. In this paper we will give an overview of the state of the art on what sets the pedestal density structure based on experimental observations as well as theoretical models. We will complement these studies with new results from Alcator C-Mod and DIII-D where the opaqueness is increased to reach values similar to those in ITER.
Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those conditions are met, a new turbulence regime is accessed where reconnection-driven energy dissipation becomes common, rather than the rare feature that it must be when they are not.
Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.
Magnetic reconnection (MR) is a fundamental process in space and laboratory plasmas. The appearance of high power lasers opens a new way to investigate MR under the relativistic condition. In this paper, relativistic collisionless MR driven by two ultra-intense lasers and a pair of asymmetric targets is studied numerically via the kinetic simulations. The static magnetic fields produced by the electron vortex structures with opposite magnetic polarities approach each other driven by the magnetic pressure and the density gradient. The antiparallel magnetic fields annihilate accompanied with the topological variation and the corresponding magnetic field energy is being dissipated to the kinetic energy of the nonthermal charged particles. Besides the outflows along the current sheet, a fast particle bunch is accelerated perpendicularly contributed by the displacement current.
Rotating spokes are observed in a partially magnetized plasma using a two-dimensional full fluid moment (FFM) model. In the present setup, where the radial electric field and plasma density gradient exist in opposite directions, it is observed that the spokes propagate in the direction of the diamagnetic drift and not the E × B drift. This is contrary to the modified Simon–Hoh instability, and the results suggest that the spokes can be driven by a strong diamagnetic drift. Different parameters, including magnetic field amplitude and physical domain size, influence the growth of the rotational instability as well as the dominant wave modes that arise. Furthermore, the propagation speed of the rotating spokes obtained from the FFM simulation are in good agreement with the observations in experimental and other computational work.
Most of the light from blazars, active galactic nuclei with jets of magnetized plasma that point nearly along the line of sight, is produced by high-energy particles, up to around 1 TeV. Although the jets are known to be ultimately powered by a supermassive black hole, how the particles are accelerated to such high energies has been an unanswered question. The process must be related to the magnetic field, which can be probed by observations of the polarization of light from the jets. Measurements of the radio to optical polarization—the only range available until now—probe extended regions of the jet containing particles that left the acceleration site days to years earlier and hence do not directly explore the acceleration mechanism, as could X-ray measurements. Here we report the detection of X-ray polarization from the blazar Markarian 501 (Mrk 501). We measure an X-ray linear polarization degree Π X of around 10%, which is a factor of around 2 higher than the value at optical wavelengths, with a polarization angle parallel to the radio jet. This points to a shock front as the source of particle acceleration and also implies that the plasma becomes increasingly turbulent with distance from the shock.
The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.
Abstract For over two decades, our fundamental understanding of energy transport dynamics in the core of tokamak plasmas had been challenged by the striking observation of temperature perturbation reversals following the injection of cold pulses at the plasma edge. These phenomena were first discovered by Gentle et al. (Phys. Rev. Lett. 74(18):3620–3623, 1995) in 1995 and had long been suggested to be evidence of nonlocal transport effects. In recent years, a new explanation to these phenomena has emerged, fully consistent with the theory of turbulent transport in magnetized plasmas and in remarkable agreement with experiment. This article reviews the experimental observation of temperature reversals in tokamak plasmas and presents the explanation based on local transport physics.
A common approach to assess the nature of energy conversion in a classical fluid or plasma is to compare power densities of the various possible energy conversion mechanisms. A leading research area is quantifying energy conversion for systems that are not in local thermodynamic equilibrium (LTE), as is common in a number of fluid and plasma systems. Here we introduce the “higher-order nonequilibrium term” (HORNET) effective power density, which quantifies the rate of change of departure of a phase space density from LTE. It has dimensions of power density, which allows for quantitative comparisons with standard power densities. We employ particle-in-cell simulations to calculate HORNET during two processes, magnetic reconnection and decaying kinetic turbulence in collisionless magnetized plasmas, that inherently produce non-LTE effects. We investigate the spatial variation of HORNET and the time evolution of its spatial average. By comparing HORNET with power densities describing changes to the internal energy (pressure dilatation, Pi − D , and divergence of the vector heat flux density), we find that HORNET can be a significant fraction of these other measures (8% and 35% for electrons and ions, respectively, for reconnection; up to 67% for both electrons and ions for turbulence), meaning evolution of the system towards or away from LTE can be dynamically important. Applications to numerous plasma phenomena are discussed. Published by the American Physical Society 2024
Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫ σ , and form a hard power-law energy distribution, f i ∝ γ −1 . This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to either σ in the strong cooling regime or the radiation burnoff limit, γ syn , when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit, ε c ≈ 16 MeV, is only efficiently produced in the strong cooling regime, γ syn < σ . In this regime, we find that the spectral cutoff scales as ε cut ≈ ε c ( σ / γ syn ) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.
Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for non-thermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic non-thermal particle acceleration in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically observed dependencies of the power-law index $p$ and high-energy cutoff $\gamma _c$ of the resulting non-thermal particle energy spectrum $f(\gamma )$ on the ambient plasma magnetization $\sigma$ , and (for $\gamma _c$ ) on the system size $L$ . In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{\rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution.
Here we study within a fully kinetic framework the generation of “seed” magnetic fields through the Weibel instability, driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and the saturation of the Weibel instability when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated fields on key parameters (e.g., ratio of system scale to electron skin depth and forcing amplitude) are confirmed by two-dimensional and three-dimensional particle-in-cell simulations of an electron–positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow and therefore has important implications for magnetogenesis in dilute astrophysical systems.
This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.
The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.
This paper investigates energetic electron transport in magnetized toroidal plasmas with magnetic fields characterized by island chains and regions of stochastic field lines produced by coil perturbations. We report on experiments performed in the DIII-D tokamak, which utilize electron cyclotron heating and current drive pulses to ‘tag’ electron populations within different locations across the discharge. The cross-field transport of these populations is then inferred from electron cyclotron emission measurements and gamma emission signals from scintillator detectors. Two types of energetic particles are distinguished and discussed: non-relativistic suprathermal electrons and relativistic runaway electrons. The magnetic field topology in each discharge is reconstructed with field-line tracing codes, which are also used to determine the location and scale of magnetic islands and stochastic regions. Comparison of simulations and experiments suggests that suprathermal transport is suppressed when the tagging is performed at a smaller radial location than the location of the $q = 1$ island chain and enhanced otherwise. Here q is the safety factor. We further demonstrate that increasing the width of the stochastic region within the edge plasma yields enhancement of the suprathermal electron transport.
Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.
This paper discusses thermal equilibrium states of single-species plasmas, such as pure electron plasmas and pure positron plasmas, that are confined in a dipole trap. Thermal equilibrium states for such plasmas are routinely realized in the homogeneous magnetic field of Penning–Malmberg traps. We generalize the theory of these states to include inhomogeneous magnetic dipole fields. The approach to thermal equilibrium takes place in two stages with well separated time scales. On the collision time scale, thermal equilibrium is established along each magnetic field line. On the much longer transport time scale, heat conduction and viscosity bring the plasmas on different flux contours into thermal equilibrium, we call this a state of global thermal equilibrium. We present numerical results for local and global thermal equilibria. These results agree with the analytic predictions for charge collections that are large compared with the Debye length. There is, in principle, no limit to the confinement time of a single-species plasma in a global thermal equilibrium state. Experiments with hot electron–ion plasmas performed in the LDX and RT1 devices give us confidence that, in contrast to a Penning–Malmberg trap, a magnetic dipole field can also confine cold quasi-neutral electron–positron pair plasmas on the time scale of the phenomena of interest. Such pair plasmas are assumed to form in the magnetosphere of neutron stars but have so far not been realized in a laboratory. The creation of an electron–positron pair plasma is the main goal of the APEX collaboration.