Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

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 487 records · Page 27

Theory of the large-amplitude plane magnetoacoustic wave propagating transverse to the magnetic field in a hot collisionless plasma

An exact solution of the kinetic and electromagnetic equations for a large-amplitude plane magnetoacoustic wave propagating transverse to the magnetic field in a hot collisionless plasma is presented. The solution gives simple relations among the magnetic-field strength, density, stress tensor, and plasma velocity, all of which are measurable in the interplanetary plasma. These relations are independent of the electron and ion velocity distributions, subject to certain restrictions on 'high-velocity tails.' The magnetic field of the wave is linearly polarized. The wave steepens to form a shock much as the analogous waves of MHD theory do.

Barnes, A.↗

Generation of plasma waves by thick-target electron beams, and the expected radiation signature

The production of plasma waves by a nonthermal beam of high-energy electrons injected into a background thermal plasma was investigated. The coupled kinetic equations for nonthermal electrons were solved, and the level of the plasma waves produced by the bump-on-tail instability arising as a result of the Coulomb collisions of nonthermal electrons with background plasma was estimated. In addition, the radiation signatures of these waves were examined. It is shown that the wave-particle interactions have a significant effect on the overall distribution of the electrons and that it is unlikely that such effects can be discerned in the observed bremsstrahlung or synchrotron radiation of the nonthermal electrons. However, a significant fraction of the energy lost by the beam is transferred to plasma waves and is then dissipated by the thermal plasma.

Hamilton, Russell J.↗

Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes

Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.

Zimmerman, Michael I.↗

Review of the 1st EUV Light Sources Code Comparison Workshop

We review the results of the 1 st Extreme Ultraviolet (EUV) Light Sources Code Comparison Workshop. The goal of this workshop was to provide a platform for specialists in EUV light source plasma modeling to benchmark and validate their numerical codes using well-defined case studies. Detailed consideration of a plethora of atomic collisional and radiative processes is required for modeling EUV light source plasmas. Eight institutions spanning four countries contributed data to the workshop. Two topics were addressed, namely (i) the atomic kinetics and radiative properties of tin plasmas under EUV-generating conditions and (ii) laser absorption in a fully ionized, one-dimensional hydrogen plasma. In this paper, we summarize the key findings of the workshop and outline plans for future iterations of the code comparison activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗

Kinetic simulations underestimate the effects of waves during magnetic reconnection

Collisionless plasma systems are often studied using fully kinetic simulations, where protons and electrons are treated as particles. Due to their computational expense, it is necessary to reduce the ion-to-electron mass ratio $m_i/m_e$ or the ratio between plasma and cyclotron frequencies in simulations of large systems. In this Letter we show that when electron-scale waves are present in larger-scale systems, numerical parameters affect their amplitudes and effects on the larger system. Using lower-hybrid drift waves during magnetic reconnection as an example, we find that the ratio between the wave electric field and the reconnection electric field scales as $\sqrt{m_i/m_e}$, while the phase relationship is also affected. The combination of these effects means that the anomalous drag that contributes to momentum balance in the reconnection region can be underestimated by an order of magnitude. The results are relevant to the coupling of electron-scale waves to ion-scale reconnection regions, and other systems such as collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microphysically modified magnetosonic modes in collisionless, high-β plasmas

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-β collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma β is increased. For mode amplitudes |δB/B 0 |≳0.3, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies |δB/B 0 |≳2β -1 . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-β plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

To Boldly Go: America's Next Era in Space. The Plasma Universe

Dr. France Cordova, NASA's Chief Scientist, chaired this, the eighth seminar in the Administrator's Seminar Series. She introduced the NASA Administrator, Daniel S. Goldin, who, in turn, introduced the subject of plasma. Plasma, an ionized gas, is a function of temperature and density. We ve learned that, at Jupiter, the radiation is dense. But, Goldin asked, what else do we know? Dr. Cordova then introduced Dr. James Van Allen, for whom the Van Allen radiation belt was named. Dr. Van Allen, a member of the University of Iowa faculty, discussed the growing interest in practical applications of space physics, including radiation fields and particles, plasmas and ionospheres. He listed a hierarchy of magnetic fields, beginning at the top, as pulsars, the Sun, planets, interplanetary medium, and interstellar medium. He pointed out that we have investigated eight of the nine known planets,. He listed three basic energy sources as 1) kinetic energy from flowing plasma such as constitutional solar wind or interstellar wind; 2) rotational energy of the planet, and 3) orbital energy of satellites. He believes there are seven sources of energetic particles and five potential places where particles may go. The next speaker, Dr. Ian Axford of New Zealand, has been associated with the Max Planck Institut fuer Aeronomie and plasma physics. He has studied solar and galactic winds and clusters of galaxies of which there are several thousand. He believes that the solar wind temperature is in the millions of degrees. The final speaker was Dr. Roger Blanford of the California Institute of Technology. He classified extreme plasmas as lab plasmas and cosmic plasmas. Cosmic plasmas are from supernovae remnants. These have supplied us with heavy elements and may come via a shock front of 10(sup 15) electron volts. To understand the physics of plasma, one must learn about x-rays, the maximum energy of acceleration by supernova remnants, particle acceleration and composition of cosmic rays, maximum acceleration, and how fast protons are heated by ions. He asked questions about where high energy cosmic rays are made, what accelerates electrons, radiates gamma rays, makes electronpositron plasma, and finally noted that pulsars are good time keepers, but we need a better understanding of their mechanism and of plasmas, both cosmic and ground-based. In the discussion period, Goldin asked if NASA should put up an x-ray interferometer. The answer was no; gamma rays are of greater interest just now. Goldin also asked what the assembled scientists would like to see for a future mission? They expressed an interest in learning more about the origin of galaxies, cosmic rays, solar systems, planets, the existence of life "out there", gamma ray sources, the nature of gamma ray bursts, and the flow of gases around black holes. The discussion concluded with a suggestion that NASA should communicate to the general public more information regarding actual technological trials and tribulations involved in getting an experiment to work. The speakers thought that this would help non-scientists to better appreciate what it is that NASA does in connection with the benefits that are achieved.

Source record↗

On the violation of the zeroth law of turbulence in space plasmas

The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. This behaviour, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation – no matter how small – can thermalise the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetised plasma turbulence such as that recently observed by the Parker Solar Probe spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is ‘imbalanced’ – when the large-scale energy input is dominated by Alfvénic perturbations propagating in one direction (the most common situation in space plasmas) – nonlinear conservation laws imply the existence of a ‘barrier’ at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong resemblance to those observed in situ , exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales. The effect thus offers a possible solution to the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The existence of the ‘barrier’ also suggests that, how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.

Physics↗

The effect of gas injection location on a lithium vapor box divertor in NSTX-U

Abstract The lithium vapor box divertor is a proposed divertor design that will minimize contamination of the upstream plasma in a fusion device, while also ensuring protection of the target. In this design lithium is evaporated near the target by high temperature lithium surfaces, dissipating the plasma heat flux. The lithium vapor box has been predicted via the fluid-kinetic code scrape off layer plasma simulator (SOLPS-ITER) to achieve low ( n L i / n e ∼ 0.05) upstream concentrations of lithium and low target heat fluxes. Here we compare two choices of deuterium gas puff location using SOLPS-ITER, the private flux region (PFR) and the common flux region (CFR), and find significant differences in the contamination level required to reach an acceptable target heat flux (defined here as q tar max ⩽ 10 MW m −2 ). Deuterium gas puffing from the PFR is seen to better reduce upstream lithium contamination. The difference in puffing location is seen to cause changes in the upstream flow of lithium ions. The PFR puff, having better access to the separatrix, can better reduce the upstream-directed flow of lithium near the separatrix which is the primary source of contamination due to a large thermal force in this region. Puffing from the CFR, partially due to inefficacy at reducing separatrix lithium flow, has higher lithium concentration within the plasma. Solutions that reduce the heat flux to below 10 MW m −2 have a range of lithium concentrations between n L i / n e ∼ 0.01–0.12 depending on puff intensity, location, evaporator temperature and recycling at the various plasma facing components. The efficacy of the puffs is tested for sensitivity to deuterium recycling coefficient at the target, evaporators, and main chamber walls. A CFR located puff is found to be more dependent on the recycling coefficients used than a PFR located puff. regardless of the set of recycling coefficients chosen, PFR puffing achieves lower lithium contamination than CFR puffing for a given heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗