Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA SHEATH”

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 55 records · Page 3

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LLNL FESP Theory Highlights: November 2024

This talk presented a refined understanding of plasma-facing thermionic emitter cathodes by examining the limitations of the conventional space charge-limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge-limited regimes based on electron current constraints, recent findings reveal that trapped ions can form a quasi-neutral "second plasma," expanding upstream and interacting with the primary plasma. These effects are amplified in non-planar geometries, such as cylindrical or spherical sheaths, where trapped ions create a positive feedback loop with the upstream plasma. To address these phenomena, a new "Aid-and-Compete" model was developed, offering insights into current enhancement and sheath mode transitions relevant for plasma devices, including emissive cathodes in electric propulsion, fusion applications, thin film deposition systems, and plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiofrequency sheath rectification on WEST: application of the sheath-equivalent dielectric layer technique in tokamak geometry *

Radiofrequency sheath rectification is a phenomenon relevant to the operation of Ion Cyclotron Range of Frequencies (ICRFs) actuators in tokamaks. Techniques to model the sheath rectification on 3D ICRF antenna geometries have only recently become available (Shiraiw et al 2023 Nucl. Fusion 63 026024; Beers et al 2021 Phys. Plasmas 28 093503). In this work, we apply the 'sheath-equivalent dielectric layer' technique, used previously only on linear devices (Beers et al 2021 Phys. Plasmas 28 103508), in tokamak geometry, computing rectified sheath potentials on the WEST ICRF antenna. Advancing the state of the art in sheath rectification modeling, we compute the sheath potentials not just on the limiters, but also on the Faraday Screen bars. The calculations show a peak rectified DC potential of 300 V on the limiters and 500 V on the Faraday screen. Assuming a typical sputtering yield curve, the RF sheath rectification increases the sputtering yield from the limiters by a factor of 2.6 w.r.t. the sputtering due to the non-rectified thermal sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A variational method for the sheath potential of hypersonic leading edges with space-charge limitations

Electron transpiration cooling for the leading edges (LE) of hypersonic aircraft utilizes thermionic emission; however, space-charge effects limit the electron emission rate, potentially diminishing the efficiency of this cooling mechanism. We develop a variational weak form of the Poisson equation that describes the sheath potential and then numerically solve it using the finite element method. This formulation has two main benefits: (1) the space-charge limit condition can be incorporated as a constraint and (2) it allows for the analysis of three-dimensional geometries with complex boundary conditions. We demonstrate that the current emitted from the surface of an LE is generally a small fraction of the Child–Langmuir limit due to space charge. We then propose several methods to enhance the emitted current from the surface and to boost the cooling effect of thermionic emission. These include increasing the plasma density, applying a negative surface potential, and using fringe fields under suitable geometric conditions. For a LaB6 emitting LE, the total emitted current is shown to be minimal and independent of the temperature of a surface with floating potential. However, when a negative potential is applied and the surface is heated, the emitted current follows the Richardson–Dushman relationship up to a critical temperature, beyond which it remains constant. At an applied surface potential of −5 V, the critical temperature is around 1700 K.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictive outlook for experiments resolving prompt vs local redeposition of high- Z materials in tokamaks

High-Z plasma facing components redeposit within the sheath through a combination of two distinct mechanisms: prompt (or geometric-driven) and local (or sheath-driven) redeposition. Experimental efforts are needed to determine the leading-order parameters influencing prompt-vs-local trade-off, which sets the fraction of material entering the scrape-off layer. In preparation for such experiments, leading-order parameters are isolated within the PYEAD-RustBCA-GITR coupled net erosion code using Sobol’ sensitivity analysis. Then, experiments resolving prompt-vs-local trade-off under variation of these leading-order parameters are proposed using an isotopic coupon design with multifaceted diagnostic coverage. The measurability of these experiments is evaluated using synthetic diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A tutorial on the micro-trench technique for incident ion angle, material erosion, and impurity deposition measurements at plasma-facing surfaces

The micro-trench method is a new technique to measure polar and azimuthal incident ion angles, material erosion, and impurity deposition at plasma-facing surfaces, and has been applied in the DIII-D divertor. This article gives a tutorial of the micro-trench technique consisting of six steps: (1) micro-trench fabrication, (2) pre-exposure measurement of the fabricated micro-trench geometry, (3) tracer material deposition, (4) plasma exposure, (5) post-exposure observation of the tracer material, and (6) post-exposure measurement of the micro-trench geometry. Two criteria need to be satisfied to apply the micro-trench method successfully: (i) uniform impurity deposition on the micro-trench floor, and (ii) erosion dominated by the physical sputtering induced by the impinging ions. Furthermore, when those two criteria are satisfied, post-exposure analysis of the impurity deposition patterns on the micro-trench floor may be used to determine the polar and azimuthal incident ion directions (mean values of the ion angle directions), erosion rate, and impurity deposition rate during plasma exposure without computational interpretation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mitigation of ICRF—edge plasma interaction in Alcator C-Mod

The parasitic interaction between ion cyclotron range of frequencies (ICRF) waves and the edge plasma has long limited high-power ICRF utilization in magnetically confined fusion plasmas. In this paper, we report on the successful mitigation of this interaction by power tapering a four-strap field-aligned (FA) antenna in the Alcator C-Mod tokamak. With the antenna operated in dipole phasing ($0\pi0\pi$), the ratio of the power coupled by the central two straps ($P_\mathrm {cent}$) to the power coupled by the outer two straps was varied at fixed total coupled ICRF power ($P_\mathrm {tot}$). With $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8-0.9$, no enhancement of the plasma potential was measured despite 1 MW of coupled ICRF power. When the antenna is powered, energy deposition increases on its four corners and on neighboring limiter tiles magnetically connected to some of those corners. The energy deposited on the antenna and limiter is minimized together with near-field plasma potentials. Additionally, we show that deviation from $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8$ enhances wave-plasma interactions in places not magnetically connected to the antenna. These far-field effects are correlated with changes in the antenna spectrum. Hence, on C-Mod, near-field and far-field interactions were minimized for about the same range of $P_\mathrm {cent}$/$P_\mathrm {tot}$. Optimal power tapering was also found to improve ICRF heating efficiency and facilitate H-mode access. Finally, we show that the heat flux pattern on the toroidally aligned antennas reverses with the direction of the RF-induced E×B flow in front of the antenna, a correlation not observed for the FA antenna. We argue that the RF-induced E×B flow, which carries high-density plasma from the near-SOL, is aligned with the FA antenna and does not intercept its corners, where the plasma potential is large. This could explain the previously reported smaller impurity generation by the FA antenna (Wukitch et al 2013 Phys. Plasmas 20 056117).

Alcator C-Mod↗

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗

Impurity transport in PISCES-RF

Linear plasma devices (LPD) utilizing a helicon plasma source, a high density light ion source, can generate impurities due to progressive erosion of the radio frequency (RF) transmission window caused by rectified sheath voltage. These source-born impurities can entrain and be transported by the plasma toward a target, affecting plasma-material interaction studies. Earlier work on material testing in Prototype-Materials Plasma Exposure eXperiment at ORNL revealed significant source impurity deposition on downstream targets. However, using a similar RF source, no target impurity deposition is observed in Plasma Interaction Surface Component Experimental Station (PISCES)-RF despite evidence of RF window erosion in the source region, thereby motivating the present work. Experimentally, using various magnetic field configurations upstream of the PISCES-RF plasma source and seeding titanium (Ti) impurities at various axial locations, impurity transport and deposition along the machine axis were investigated. It was found that Ti deposition was localized to the side of the plasma source where the Ti impurity was seeded. In contrast, aluminum (Al) deposition, originating from the sputtering of the helicon window, occurred predominantly upstream of the plasma source, suggesting an asymmetry in the axial transport of eroded RF window material. These observations suggest a stagnation of the parallel plasma flow immediately downstream of the plasma source, with impurity ions remaining unmagnetized near the source upstream. Al deposition in magnetic field-free regions in PISCES-RF indicates that sputtered Al impurities likely remained neutral due to their large ionization mean-free path under PISCES-RF conditions. Plasma modeling and simulation supported this, indicating that Al-neutrals transport toward the helicon source upstream for low electron density cases. It was found that the Larmor radius of the Al ions was greater than the plasma radius towards the source upstream and remained weakly magnetized in PISCES-RF, meaning that plasma source-born impurities are not efficiently entrained in the plasma flow. These findings provide critical insights into impurity transport in helicon plasma-based LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory of the thermionic current beyond the traditional space charge limit enabled by trapped ions in the virtual cathode

We show that trapped ions in virtual cathode potential wells can raise the transmitted current of emitted electrons into a plasma much closer to the full emission than is predicted by cathode sheath theories without trapped ions. The transmitted current is controlled by the well barrier voltage, which must adjust to balance the creation of low-energy ions within the well, and their loss. Our model considers the case of a plasma-facing cathode where trapped ions are created passively via charge-exchange collisions and lost passively via thermal leakage over the well. We quantify these rates and estimate the current in terms of system parameters for thermionic emission into a plasma with several cathode geometries. A general prediction is that the current as a function of emitted flux does not saturate at the traditional space charge limit (the onset of a well) but can reach far higher values until the trapped ion balance breaks down, causing instability. The maximum stable current depends on parameters but in principle can be arbitrarily high if active techniques are used to manipulate the trapped ion balance. We conclude that major improvements in plasma technologies with hot cathodes might be achieved by optimizing the current enhancement enabled by trapped ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Proton beam enhancements from ultrafast laser interactions with compound parabolic concentrators

Target normal sheath acceleration (TNSA) is among the most widely studied laser–plasma ion acceleration mechanisms. In this Letter, we report on studies of proton acceleration from flat Cu targets with cone-like focusing structures called compound parabolic concentrators (CPCs) coupled to their front surface. The CPC acts as a non-imaging focusing optic that enhances the laser intensity at the target's front surface by ∼2×, from 9 × 10 18 to 2 × 10 19 W/cm 2 . This effect drives a reduction in the particle source size from the typical laser spot size of 120 μm to the CPC tip size of 65 μm and significant enhancements in the peak energy and temperature (2.6 × increase) of the resulting TNSA proton beam. We also observe increased opening angle of the beam. 2D PIC simulations have been conducted and replicate the experimental behavior. These increases have implications for long focal length facilities that require higher energy proton beams, which now can be achieved without significant infrastructure changes.

3D printing↗

Experimental parameters for plasma wakefield acceleration in a narrow plasma column

Recent theoretical advancements propose multiple positron acceleration schemes in plasma wakefield acceleration (PWFA). One of the most promising ideas involves the creation of an electron-driven blowout wake within a finite-radius pre-ionized plasma column. This leads to the formation of an elongated region of sheath electrons at the closing of the first wake period capable of accelerating positrons while simultaneously providing a transverse focusing force. Additionally, the proposed scheme has shown to be a potential means of suppressing instabilities. We present an experimental opportunity to explore the narrow column PWFA at the Facility for Advanced Accelerator Experimental Tests II (FACET-II): the E333 experiment. As a pivotal first step towards achieving positron acceleration in PWFA, we have planned a precursor experiment utilizing the currently available single-bunch electron beam to study the physics of the narrow plasma PWFA scheme. We outline the feasible experimental parameters, including beam and ionization laser parameters, along with the required laser optics for the experiment. The primary observable signatures for this stage of the experiment are the final energy spectrum and transverse position of the electron bunch, anticipating reduced energy loss and enhanced beam guidance in the narrow plasma column compared to the nominal PWFA. Comprehensive simulations are used to detail our experimental plan.

Lee, Valentina [University of Colorado, Boulder, C↗

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗