Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma material interaction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX. I. Sheath equivalent dielectric layer for modeling the RF sheath

The pulsed linear plasma device Prototype Material Plasma Exposure eXperiment (Proto-MPEX) uses a radio frequency (RF) helicon antenna with an aluminum nitride ceramic window for plasma production. The RF sheath created under the helicon antenna is sufficient to cause ion impact energies to be greater than the sputtering threshold of the AlN helicon window material and for impurities to be created. In this study, we investigate the RF sheath on the inner diameter of the helicon window and its impact on the impurity production rates in Proto-MPEX. Three models—a 3D COMSOL finite element RF model of the Proto-MPEX helicon region, a rectified DC sheath potential model, and the 3D Global Impurity TRansport code—are coupled together to study impurity production and transportation. This novel method of impurity generation and transport modeling spans length scales ranging from the sheath (millimeters or less) up to the full device (meters) and can be applied to other radio frequency sources and antennas in a wide range of plasma physics studies, including basic plasmas, low-temperature processing plasmas, plasma thrusters, and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evidence on the effects of main-chamber neutrals on density shoulder broadening

Evidence that density shoulder broadening is dependent on high main-chamber neutral density is presented. Shoulder broadening does not occur when the sources for main-chamber neutrals are minimized using divertor baffles and wide gaps to the first wall (~3× the density decay length). Removing the baffles or reducing the gap to the inner wall both act to increase the density shoulder amplitude in otherwise identical TCV discharges. Radial turbulent transport is correlated with shoulder amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Continuum kinetic investigation of the impact of bias potentials in the current saturation regime on sheath formation

In this work, we examine sheath formation in the presence of bias potentials in the current saturation regime for pulsed power fusion experiments. It is important to understand how the particle and heat fluxes at the wall may impact the wall material and affect electrode degradation. Here, simulations are performed using the 1X-1V Boltzmann–Poisson system for a proton–electron plasma in the presence of bias potentials ranging from 0 to 10 kV. The results indicate that the sheath near the anode remains generally the same as that of a classical sheath without the presence of a bias potential. However, the sheath near the cathode becomes more prominent with a larger potential drop, a significant decrease in the electron density, and larger sheath lengths. The spatially constant current density increases to a saturation value with the increasing bias potential. For high bias potentials, the material choice needs to consider that the anode faces significantly larger particle and heat fluxes compared to the cathode. In general, the results trend with theory with differences attributed to the simplified assumptions in the theory and the kinetic effects considered in the simulations. Due to the significant computational cost of a well resolved 1X-2V simulation, only one such simulation is performed for the 5 kV case showing higher current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Estimates of global recycling coefficients for LTX- β discharges

We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dependence of high- Z redeposition on the field-to-surface pitch angle and other sheath parameters in tokamaks

Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Voltage waveform tailoring for high aspect ratio plasma etching of SiO 2 using Ar/CF 4 /O 2 mixtures: Consequences of low fundamental frequency biases

The use of non-sinusoidal waveforms in low pressure capacitively coupled plasmas intended for microelectronics fabrication has the goal of customizing ion and electron energy and angular distributions to the wafer. One such non-sinusoidal waveform uses the sum of consecutive harmonics of a fundamental sinusoidal frequency, f 0 , having a variable phase offset between the fundamental and even harmonics. In this paper, we discuss results from a computational investigation of the relation between ion energy and DC self-bias when varying the fundamental frequency f 0 for capacitively coupled plasmas sustained in Ar/CF 4 /O 2 and how those trends translate to a high aspect ratio etching of trenches in SiO 2 . The fundamental frequency, f 0 , was varied from 1 to 10 MHz and the relative phase from 0° to 180°. Two distinct regimes were identified. Average ion energy onto the wafer is strongly correlated with the DC self-bias at high f 0 , with there being a maximum at φ = 0° and minimum at φ = 180°. In the low frequency regime, this correlation is weak. Average ion energy onto the wafer is instead dominated by dynamic transients in the applied voltage waveforms, with a maximum at φ = 180° and minimum at φ = 0°. The trends in ion energy translate to etch properties. In both, the high and low frequency regimes, higher ion energies translate to higher etch rates and generally preferable final features, though behaving differently with phase angle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2021

Ion heating in the PISCES-RF liquid-cooled high-power, steady-state, helicon plasma device

Radio frequency (RF) driven helicon plasma sources are commonly used for their ability to produce high-density argon plasmas ( n > 10 19 m -3 ) at relatively moderate powers (typical RF power < 2 kW). Typical electron temperatures are <10 eV and typical ion temperatures are <0.6 eV. A newly designed helicon antenna assembly (with concentric, double-layered, fully liquid-cooled RF-transparent windows) operates in steady-state at RF powers up to 10 kW. We report on the dependence of argon plasma density, electron temperature and ion temperature on RF power. At 10 kW, ion temperatures >2 eV in argon plasmas are measured with laser induced fluorescence, which is consistent with a simple volume averaged 0D power balance model. 1D Monte Carlo simulations of the neutral density profile for these plasma conditions show strong neutral depletion near the core and predict neutral temperatures well above room temperatures. The plasmas created in this high-power helicon source (when light ions are employed) are ideally suited for fusion divertor plasma-material interaction studies and negative ion production for neutral beams.

Physics↗
DOE OSTI2020

Deuterium retention in Be-D co-deposits formed over an ITER relevant parameter space

Be-D co-deposit samples are produced over an extended range of plasma-material interaction parameters reflective of the conditions expected in ITER deposition dominant locations. Deposition temperature, D impact energy, atomic Be deposition rate, and D2 neutral pressure were varied in the ranges 373–623 K, ~2–100 eV, 1–10 × 10 19 m 2 s –1 , and 0.27–13.3 Pa respectively. Thermal desorption spectrometry was carried out on samples, subsequently, to measure the D release profile and retention. In this work, it was found that the D retention in samples agrees well with the existing GDT D/Be scaling over the extended parameter space, and that increased D impact energy and D2 pressure can lead to the observation of a sharp low-temperature release peak at ~440 K in desorption data. Desorption data without the presence of the sharp release mechanism are modeled with TMAP-7 utilizing three traps of de-trapping energies, 0.8, 0.9 and 1.0 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Large radiation back-flux from Monte Carlo simulations of fusion neutron–material interactions

Abstract Fusion power reactors will generate intense neutron fluxes into plasma-facing and structural materials (SMs). Radiation back-fluxes, generated from neutron–material interactions under these fluxes, can dramatically impact the plasma dynamics, e.g. by seeding runaway electrons during disruptions via Compton scattering of background electrons by wall-emitted gamma radiation. Here, we quantify these back-fluxes, including neutrons, gamma rays, and electrons, using Monte Carlo calculations for a range of SM candidates and first wall (FW) thicknesses. The radiation back-flux magnitudes are remarkably large, with neutron and gamma radiation back-fluxes on the same order of magnitude as the incident fusion neutron flux. Electron back-fluxes are two orders of magnitudes lower, but are emitted at sufficiently high energies to impact the sheath and boundary plasma dynamics. Material configuration plays a key role in determining back-flux magnitudes. The SM chiefly determines the neutron back-flux magnitude, while the FW thickness principally attenuates the gamma ray and electron back-fluxes. In addition to prompt back-fluxes, which are emitted immediately after fusion neutrons impact the surface, significant delayed gamma ray and electron back-fluxes arise from nuclear decay processes in the activated materials. These delayed back-flux magnitudes range from 2% to 7% of the prompt back-fluxes, and remain present during transients when fusion no longer occurs. During disruptions, build-up of delayed gamma radiation back-flux represents potential runaway electron seeding mechanisms, posing additional challenges for disruption mitigation in a power reactor compared with non-nuclear plasma operations. This work highlights the impact of these radiation back-fluxes plasma performance and demonstrates the importance of considering back-flux generation in materials selection for fusion power reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantitative Measurement of Positive and Negative Ion Species Ejected from a Li–O–H Surface by Hydrogen and Noble Gas Ion Irradiation

We report sputtering yields of Li + , H - , O - , and OH x - ion species from an Li–O–H surface for H, D, He, Ne, and Ar ion irradiation at 45° incidence in the energy range of 30–2000 eV. A Li film was deposited on a stainless steel target using Li evaporators in the LTX-β vessel, using the LTX-β Sample Exposure Probe (SEP), which includes an ultrahigh vacuum suitcase for transferring targets without significant contamination from air exposure. The SEP was used to transfer the Li-coated target from LTX-β to a separate Sample Exposure Station (SES) to perform ion exposure measurements. The SEP was also used for characterization of the Li-coated target utilizing X-ray photoelectron spectroscopy in a different chamber, showing that the lithium film surface was oxidized. Ion exposures were performed using an electron cyclotron resonance plasma source in the SES. Sputtered/ejected species were sampled by a quadrupole mass spectrometer with capabilities for detecting positive and negative ions, and an energy filter for determining the mean kinetic energy of the ejected ion species. All ion irradiations caused Li + ions to be ejected, while causing impurity ions such as H + , H - , O - and OH - to be ejected. Measured ion energies of Li + ions from a Li–O–H surface suggested that the typical sheath potential on the divertor surface can trap sputtered Li + ions, which were previously reported as ~ 60% of total sputtered Li species from Li targets (Allain and Ruzic in Nucl Fusion 42:202, 2002). Hence, our results for the sputtering yields of ejected ion species and their associated ion energies from a Li–O–H surface indicates that lithium sputtering is suppressed and impurity removal is enhanced due to the sheath potential at the divertor surface for fusion reactor applications.

Plasma-material interaction↗