Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space plasma”

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 325 records · Page 18

An overview of 3D field optimization for control of transport and edge instabilities on KSTAR

An international team from several laboratories and universities has made key advances over the last few years in the control of plasma transport and edge instabilities with applied 3D fields in the KSTAR tokamak to optimize long pulse operation scenarios. This overview begins with the optimization of both core and edge resonant magnetic perturbations (RMP) to improve fast ion confinement to avoid excessive limiter heat loads due to fast ion losses and successful modeling of the experimental results. Integrated and advanced plasma control techniques with machine learning (ML) and adaptive control were then used to optimize the 3D field spectrum in real-time to control edge localized modes (ELMs) while avoiding core locked modes that could disrupt the plasma. Accelerating the offline model of 3D fields with a surrogate ML model can optimize ELM suppression in the edge while limiting the impact of the applied RMP fields deeper in the plasma core in real-time. In addition, the impact of the 3D fields on the divertor heat load has been modeled and compared with experimental measurements. An analysis of a multi-machine database including KSTAR has been performed to better understand the metrics for the observed RMP thresholds for ELM suppression and the resulting plasma performance. Predictive modeling of the operational space for ELM suppression and density pumpout due to RMP has shown the importance of magnetic islands in the plasma edge and their impact on plasma turbulence. This research has culminated in the development of successful long pulse operational scenarios on KSTAR while attempting to overcome challenges of the new tungsten divertor.

3D fields↗

Operando XPS for Plasma Process Monitoring: A Case Study on the Hydrogenation of Copper Oxide Confined under h-BN

Here, we demonstrate that ambient pressure x-ray photoelectron spectroscopy (APXPS) can be used for in situ studies of dynamic changes in surface chemistry in a plasma environment. This opens a new and vast application space for XPS and greatly complements modern spectroscopy techniques to probe plasma-solid/liquid interactions relevant to process monitoring in the semiconductor industry, bio-medical plasma applications and plasma remediation technologies. Hexagonal boron nitride (h-BN) grown on Cu was used in this study as a well-defined model system for plasma process monitoring and because of its unique chemical, optical and electrical properties that make it a prospective material for advanced electronics. To better understand the stability and surface chemistry of h-BN during plasma assisted processing, we track in real time the plasma-induced chemical state changes of B, N and the underlying Cu substrate using APXPS equipped with an AC discharge plasma source operating at 13 Pa. Residual gas analysis (RGA) mass-spectra were concurrently collected during plasma-XPS to track reaction products formed during plasma exposure. A clear reduction of Cu x O is seen, while an h-BN layer remains intact, suggesting hydrogen radical (H • ) species can attack the exposed and h-BN covered Cu oxide patches and partially reduce the underlying substrate without significantly damaging the overlaying h-BN, which is of practical importance for development of h-BN encapsulated devices and interfaces. In addition to demonstration of plasma-XPS capabilities we discuss the observed challenges (e.g., parasitic plasma-chamber walls reactions and charging effects) and propose potential solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Variable-spectrum mode control of high poloidal beta discharges

DIII-D experiments demonstrate that high pressure, broad current profile equilibria can be accessed in the high poloidal beta regime by optimizing the MHD mode control poloidal spectrum. A novel, variable spectrum (VS) magnetic feedback scheme implemented using the DIII-D internal non-axisymmetric coils (I-coils) facilitated access to reduced internal inductance $l$ i operation above the no-wall beta limit compared with both no feedback and fixed spectrum feedback. In addition, the VS feedback helped avoid beta collapses caused by marginally unstable resistive wall mode activity. The lower and upper I-coil rows were configured in two independent feedback loops, allowing the feedback field's poloidal spectrum to vary and track changes in the plasma mode structure as the edge safety factor q 95 varied from 11 to 6 during the discharges. The q 95 dependence of the measured phase difference between the lower and upper I-coil rows during VS feedback is qualitatively compatible with ideal MHD simulations of the least-stable plasma kink mode and with plasma response simulations that included kinetic modifications to ideal MHD. The VS feedback approach is a straightforward way to improve resilience to variations in mode structure that occur as plasma parameters change. The demonstrated expansion of the operating space to lower $l$ i is expected to improve the coupling of the plasma kink mode to external fields and beneficial wall eddy currents, and is compatible with high bootstrap fraction operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulence driven widening of the near-SOL power width in ASDEX Upgrade H-Mode discharges

Divertor heat flux data from infra red (IR) from various tokamaks in H-Mode regime show that the power width, $\lambda_{q}$, scales approximately like 1.6 $(R/a) ρ_{s,pol}$. However, the IR based scaling comes with the restriction that only discharges were considered with attached divertor conditions. Here we set up a new data base covering low edge density data but also including high density plasmas reaching the H-mode density limit. We use Thomson-scattering to measure the electron temperature decay length which will set the near-SOL power width through parallel heat conduction. We present scaling results for the density, temperature and pressure separatrix fall-off lengths. As the principal result we derive a generalized power width scaling which reads as $\lambda_{q}\propto\,\rho_{s,pol}\cdot(1+2.1\alpha_{t}^{1.7})$ where $\alpha_{t}$ describes a normalized collisionality ($\alpha_{t} = 3\cdot\,10^{-18}\,R\,q^2\,n\,Z_{eff}T^{-2}$). The parameter $\alpha_{t}$ describes the relative importance of the interchange effect on drift-wave turbulence as proposed by Scott and is found to be closely linked to the diamagnetic parameter αd introduced by Rogers, Drake and Zeiler (RDZ). The edge operational space for ASDEX Upgrade displays a boundary qualitatively in agreement with the region labelled as density limit in the edge plasma phase space proposed by RDZ. The data base shows in the limit of low edge densities ($\alpha_{t}\simeq\,0.2$) accurate agreement to the IR based power width scaling. At elevated separatrix densities ($\alpha_{t}\simeq\,1$) the power width is broadened by a factor of up to about three accompanied by a reduction of the $\text{H}_{98,y2}$ confinement factor. Finally, plasmas with higher triangularity show higher values for $\text{H}_{98,y2}$ in the same range of normalized collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simultaneous access to high normalized density, current, pressure, and confinement in strongly-shaped diverted negative triangularity plasmas

Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized density, current, pressure, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the key parameters, expanded operating spaces featuring edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are observed, even simultaneously, and all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗

Radio Frequency Transients Correlated with Electron Flux Measured On-Board the STP-Sat6

We analyzed transient electromagnetic signatures, referred to as spacecraft environment discharges (SEDs) or on-board discharges (OBDs), and the in-situ particle environment simultaneously measured by the Radio Frequency Sensor (RFS) and the Space and Atmospheric Burst Reporting System Z-plasma spectrometer (SABRS-ZPS), respectively. These sensors are on-board the U.S. Department of Defense Space Test Program Satellite 6 (STP-Sat6) in geostationary orbit. Radio frequency signatures of SEDs on STP-Sat6 have relatively high peak amplitudes; they often caused amplitude-saturation in the so-called high band (120–140 MHz frequency range) of the RFS likely due to the proximity of their source to the RFS’s dipole antennas. We examined the occurrence rate of SEDs reported by the RFS and its temporal relationship to electron flux in the 7.9–12.2 keV range. Peaks in SED count rates were time-correlated with peaks in electron flux. The electron-flux peaks appeared to lead the SED count rate peaks by 24–45 min for the 7.9–12.2 keV energies. This indicates that incidence of electrons at this relatively low-energy range likely resulted in ramping up of the charging of the spacecraft leading to periods of high occurrence rates of discharges in and around the spacecraft.

54 ENVIRONMENTAL SCIENCES↗

Plasma Imaging, LOcal Measurement, and Tomographic Experiment (PILOT): A Mission Concept for Transformational Multi-Scale Observations of Mass and Energy Flow Dynamics in Earth’s Magnetosphere

We currently do not understand the fundamental physical processes that govern mass and energy flow through the Earth’s magnetosphere. Knowledge of these processes is critical to understanding the mass loss rate of Earth’s atmosphere, as well as for determining the role that a planetary magnetic field plays in atmospheric retention, and therefore habitability, for Earth-like planets beyond the solar system. Mass and energy flow processes are challenging to determine at Earth in part because Earth’s planetary magnetic field creates a complex “system of systems” composed of interdependent plasma populations and overlapping spatial regions that perpetually exchange mass and energy across a broad range of temporal and spatial scales. Further, the primary mass carrier in the magnetosphere is cold plasma (as cold as ~0.1 eV), which is invisible to many space-borne instruments that operate in the inner magnetosphere. The Plasma Imaging LOcal and Tomographic experiment (PILOT) mission concept, described here, provides the transformational multi-scale observations required to answer fundamental open questions about mass and energy flow dynamics in the Earth’s magnetosphere. PILOT uses a constellation of spacecraft to make radio tomographic, remote sensing, and in-situ measurements simultaneously, fully capturing cold plasma mass dynamics and its impact on magnetospheric systems over an unprecedented range of spatial and temporal scales. This article details the scientific motivation for the PILOT mission concept as well as a potential mission implementation.

79 ASTRONOMY AND ASTROPHYSICS↗

A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.

divertor physics↗

The evolution of curvature in planar, photoionization-driven heat fronts

Photoionized plasmas are common in astrophysics and cosmology, especially in space near compact objects, and there are effects from photoionization in high-energy-density plasmas due to the large radiation fields present. Photoionized plasmas are an active area of laboratory research and there are currently experiments to study photoionization-supported heat fronts. These photoionization fronts differ from the physics of diffusive radiation waves, commonly called Marshak waves, that are also an active area of research. This work uses a geometric argument to describe the expected evolution of the photoionization front curvature, in a planar geometry. It then compares this curvature to that of a Marshak wave as a method of diagnosing a heat front experiment. It is found that while the curvature of a planar Marshak wave increases in time, it decreases for a photoionization front. A comparison of radiation energy and electron heat fluxes through the container for the heat front propagating medium demonstrates that the geometric argument for the photoionization front curvature is sufficient. Furthermore, this comparison also demonstrates that wall losses are not significant in a photoionization front because the post-front region is very optically thin. A discussion of the implication this work has on material choice in the targets for an experiment follows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Reversed magnetic shear scenario development in NSTX-U using TRANSP

Abstract Understanding and control of electron thermal transport is a critical point of research in magnetic fusion experiments. Previous experiments have shown that operation with reversed magnetic shear (RMS) can suppress electron thermal transport, resulting in the generation of internal transport barriers (ITBs), with the location of the ITB correlated with the location of minimum magnetic shear. The recent upgrades to NSTX—increased magnetic field up to 1 T, increased plasma current up to 2 MA, 2nd neutral beam—present an increased operating space in which to explore electron thermal transport in RMS plasmas. Utilizing TRANSP, we have developed operating scenarios by which to generate RMS in NSTX-U. The results suggest that RMS in NSTX-U can be generated through fast current ramp and early beam injection into a large plasma volume. This is very similar to the procedure that was followed in both TFTR and NSTX to generate RMS. Sustainment of RMS, disregarding non-( q min = 1) MHD events, requires maintaining a large plasma volume, and increasing the core T e , either via increased plasma current and/or adding heating power. Using this procedure, RMS was sustained for ∼1 s, with q min > 1 for that period.

Galante, M. E. (ORCID:0009000098149425)↗

Topological phases and bulk-edge correspondence of magnetized cold plasmas

Abstract Plasmas have been recently studied as topological materials. However, a comprehensive picture of topological phases and topological phase transitions in cold magnetized plasmas is still missing. Here we systematically map out all the topological phases and establish the bulk-edge correspondence in cold magnetized plasmas. We find that for the linear eigenmodes, there are 10 topological phases in the parameter space of density n , magnetic field B , and parallel wavenumber k z , separated by the surfaces of Langmuir wave-L wave resonance, Langmuir wave-cyclotron wave resonance, and zero magnetic field. For fixed B and k z , only the phase transition at the Langmuir wave-cyclotron wave resonance corresponds to edge modes. A sufficient and necessary condition for the existence of this type of edge modes is given and verified by numerical solutions. We demonstrate that edge modes exist not only on a plasma-vacuum interface but also on more general plasma-plasma interfaces. This finding broadens the possible applications of these exotic excitations in space and laboratory plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗

Calibration and Initial Results of Space Radiation Dosimetry Using the iMESA-R

The integrated Miniaturized Electrostatic Analyzer-Reflight (iMESA-R) is a space weather instrument designed to measure plasma density, temperature, and spacecraft charging, along with total ionizing dose and dose rate. A constellation consisting of five nearly identical instruments has been designed and developed to act as science payloads hosted on-board Department of Defense satellites operating in low-Earth orbit. To validate the dosimetry component of the iMESA-R, a radiation test study was performed to calibrate the dosimeter and ascertain the attenuation due to the instrument aluminum housing. An 80 curie Cobalt-60 radioisotope source emitting beta, X-ray, and gamma ray radiation decay products was used to calibrate the dosimeter response as a function of distance and instrument shielding. We present results of the calibration study along with initial on-orbit data presented from the first operational iMESA-R hosted on-board Space Test Program Satellite 5 in a polar Earth orbit. The initial on-orbit data demonstrate the ability to map the radiation environment, particularly the South Atlantic Anomaly and the auroral regions using total ionizing dose rate due to electrons of E > 1.5 MeV and protons of E > 25 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of temporal variations in plasma conditions on the electric potential near self-organized dust chains

Self-organization of dust grains into stable filamentary dust structures (or “chains”) largely depends on dynamic interactions between individual charged dust grains and complex electric potential arising from the distribution of charges within a local plasma environment. Recent studies have shown that the positive column of the gas discharge plasma in the Plasmakristall-4 (PK-4) experiment at the International Space Station supports the presence of fast-moving ionization waves, which lead to variations of plasma parameters by up to an order of magnitude from the average background values. The highly variable environment resulting from ionization waves may have interesting implications for the dynamics and self-organization of dust particles, particularly concerning the formation and stability of dust chains. Here, we investigate the electric potential surrounding dust chains in the PK-4 experiment by employing a molecular dynamics model of the dust and ions with boundary conditions supplied by a particle-in-cell with Monte Carlo collision simulation of the ionization waves. The model is used to examine the effects of the plasma conditions within different regions of the ionization wave and compare the resulting dust structure to that obtained by employing the time-averaged plasma conditions. Here, the comparison between simulated dust chains and experimental data from the PK-4 experiment shows that the time-averaged plasma conditions do not accurately reproduce observed results for dust behavior, indicating that more careful treatment of plasma conditions in the presence of ionization waves is required. It is further shown that commonly used analytic forms of the electric potential do not accurately describe the electric potential near charged dust grains under these plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conservative velocity mappings for discontinuous Galerkin kinetics

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.

Discontinuous Galerkin↗