Search NASA⌕ Search

SEARCH · Search NASA

Results for “RF 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurement and modeling of the radio frequency sheath impedance in a large magnetized plasma

The DC and RF properties of radio frequency (RF) driven sheaths were studied in the Large Plasma Device (LAPD) at the University of California, Los Angeles. The experiments diagnosed RF sheaths on field lines connected to a grounded plate at one end and an ion cyclotron range of frequencies (ICRF) antenna at the other end. The experimental setup permitted measurement of the RF sheath impedance at the plate as a function of DC sheath voltage, with the latter controlled by varying the RF current applied to the antenna. Here, the DC current-voltage characteristics of these sheaths and the RF sheath impedance measurements were compared with modeling. Hot electrons, present in the LAPD plasma, were inferred to contribute significantly to both the DC and RF currents and hence the RF impedance. It was postulated that at very low power hot electrons could not access the region of the plasma subject to RF waves, resulting in an increased RF impedance. Within some experimental limitations and significant assumptions, an RF sheath impedance model was verified by the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High-energy ballistic electrons in low-pressure radio-frequency plasmas

This study demonstrates the presence of a small number of high-energy ballistic electrons (HEBEs) that originate from secondary electrons in low-pressure radio-frequency (rf) plasmas. The kinetic behaviors of the HEBEs are illustrated through electron energy probability functions from the fully kinetic particle-in-cell simulations, showing two wavy high-energy tails and two bifurcations during one rf cycle. Test-particle simulations and a semi-analytical method associated with nonlocal electron kinetics are performed to characterize the HEBE trajectories, which reveal the ballistic nature of the HEBEs and their typical bouncing features between the rf sheaths. Parameter dependence of the HEBEs on the discharge conditions (e.g., gas pressure, gap distance, and rf frequency) are identified, which is relevant to the plasma collisionality. With a pronounced presence of HEBEs, the overall impacts of the secondary electron emission on discharge parameters, such as electron power absorption and ionization rate, are also illustrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The geometry of the ICRF-induced wave–SOL interaction. A multi-machine experimental review in view of the ITER operation

As a part of ITPA-Integrated Operational Scenario activities, this contribution reviews recent experimental characterizations of radio-frequency (RF)-induced scrape-off layer (SOL) modifications of various tokamaks worldwide and of the LArge Plasma Device at UCLA. The phenomenology, as observed using a large variety of measurement techniques, is consistent with expectations from RF sheath rectification. Emphasis is then put on the complex three-dimensional (3D) spatial patterns of the RF–SOL interaction, in relation to the magnetic topology and the spatial distribution of RF currents over the metallic structures surrounding the RF wave launchers. Dependence on the local plasma parameters in the antenna vicinity is also briefly addressed. The final part discusses implications for future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigation of two-dimensional radio-frequency sheath properties using a microscale fluid model

Abstract In previous work (Kohno and Myra 2023 Comput. Phys. Commun. 291 108841), we developed a numerical scheme based on a two-dimensional microscale radio-frequency (RF) sheath model with periodically curved wall boundaries. Here, we expand the capability of this scheme through modification of the boundary conditions (BCs) on the conducting walls, which allows the ion flow to turn back to the plasma at locations on the walls where the electromagnetic force on the ions is reversed from its usual direction. Numerical simulations are carried out to investigate the dependences of the surface-integrated admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle, and to visualize the sheath structures in several cases. One of the main results is the ion cyclotron admittance resonance observed under the condition of low ion mobility (high normalized frequency). It is shown that the amplitude of the resonance peak depends on the wall bump height and the ion velocity is reversed on the sides of the bump in an RF cycle for the resonance cases. Furthermore, the differences in the admittances between the one- and two-dimensional microscale models are assessed for the purpose of understanding non-locality of the sheath near the wall surface for the parameters considered in this study. This information will be essential for improving the sheath BC for macroscale calculations in the future.

Kohno, H. (ORCID:0000000332690010)↗
DOE OSTI2023

Benchmark between antenna code TOPICA, RAPLICASOL and Petra-M for the ICRH ITER antenna

ITER will be equipped with three plasma heating systems: neutral beam (NB), electron cyclotron (EC), and ion cy-clotron resonance heating (ICRH). The latter consists of two identical ICRH antennas to deliver 20 MW to the plasma (baseline, upgradable to 40 MW). ICRH will play a crucial role in the ignition and sustainment of burning plasmas in ITER. A high fidelity and robust modeling effort to understand the interaction of the IC waves with the scrape-off-layer (SOL) plasma is a very important aspect. Among the main important research topics, we have the assessment of the antenna loading for different plasma scenarios, the role of the lower hybrid resonance in front of the antenna and how to include it in our models, and the RF sheath boundary conditions to evaluate the antenna impurity generation. In this work, we tackle the first of these by reporting on ICRF simulations employing the Petra-M code, which is an electromagnetic simulation tool for modeling RF wave propagation based on MFEM [http://mfem.org] for the ITER ICRH antenna. Moreover, a benchmark between the well tested antenna codes TOPICA, RAPLI-CASOL, which is based on COMSOL [www.comsol.com], and the Petra-M code is also presented. In conclusion, S- and Z-matrices and wave electric field are compared showing an excellent agreement among these codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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↗
DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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↗
DOE OSTI2025

Erosion of high-Z refractory coatings for helicon plasma source window

Proto-MPEX (Materials Plasma Exposure eXperiment), a linear plasma device (LPD), using a high power radio frequency (RF) (⩾100 kW, 13.56 MHz) helicon plasma source, has suffered RF sheath-induced window erosion. The sputtered impurities from the window surface transport toward the downstream target affecting plasma-material interaction studies. The rectified sheath voltage formed was high enough to cause window erosion by light ions due to its low-Z components (e.g. Si 3 N 4 and AlN). Hence, we are proposing impurity mitigation strategies, which involve the application of a Faraday screen to lower the rectified sheath voltage and the application of high-Z refractory coatings on the existing window plasma-facing surface. In this work, we report the erosion of two different coatings, i.e., tantalum oxide (Ta 2 O 5 ) and hafnium oxide (HfO 2 ) on a silicon nitride (Si 3 N 4 ) window material under conditions of low-energy deuterium (D) ion impact, well below the Ta 2 O 5 sputtering energy threshold (i.e. 250 eV). The test samples were RF-biased and exposed to a high D-ion fluence (∼10 26 m −2 ) in Plasma Interaction with Surface Component Experimental Station (PISCES-A) LPD. The experimentally measured sputtering yields of the high-Z refractory coatings below the sputtering energy threshold of Ta 2 O 5 indicate an increased erosion due to the plasma impurities, especially due to oxygen. Improving the vacuum level could reduce the oxygen impurities and increase the lifespan of the coated helicon window for plasma operation. Post-surface analysis indicates no preferential erosion of oxygen or enrichment of the high-Z surface component. This is likely due to an effective energy transfer from the impurity ion for sputtering of the high-Z component in the coating. With the reduced rectified sheath voltage at the window surface and improved vacuum conditions, the proposed high-Z refractory coatings offer a promising solution for reducing window erosion in future MPEX plasma operations at ORNL.

RF bias↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Comparison between Alcator C-Mod ICRF experiments and 3D full wave simulations

Reliable modeling of ion cyclotron range of frequencies (ICRF) antenna performance is essential for interpreting present experiments and guiding the design of future reactors. In this work, a 3D model of the Alcator C-Mod field-aligned antenna is implemented in the Petra-M finite-element framework [S. Shiraiwa et al 2023 Nucl. Fusion 63 026024] and benchmarked against experimental results. Four experimental cases are examined. First, the simulated rectified sheath potentials on the antenna limiters are compared with measurements from a power tapering experiment. Second, proof-of-principle far-field sheath simulations are performed. In scenarios with low single-pass absorption, simulations predict enhanced sheath potentials on a distant poloidal limiter in the far field of the antenna, consistent with experimental observations. Third, the simulated antenna loading during edge localized modes follows experimental trends and appears to be dominated by the density gradient at the pedestal. Finally, Petra-M predicts the unintended excitation of high-$k$ ∥ modes during monopole phasing operation, in agreement with experimental evidence of poor wave coupling and accessibility to the plasma core. Overall, reasonable agreement is found between the simulations and experiments. At the same time, areas of imperfect agreement are identified. These provide important guidance on the limits of the current state-of-the-art modeling, which should be kept in mind when using it as a predictive tool for future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The effect of high-power transient events on tungsten and tungsten coatings used for radio frequency launcher applications

High-temperature plasma-facing material coatings used for radio frequency (RF) launchers need to be robust enough to survive RF breakdown arcing or other transient events from the plasma (e.g., an edge localized mode) without causing a catastrophic failure of the coating. High-power transient effects are being explored by using an RF-induced vacuum arc to determine the robustness of tungsten coatings made by a variety of manufacturing methods. A 1/4-wavelength resonant section of vacuum transmission line terminated with an open circuit electrode structure with a well-defined electric field (30-60 kV/mm) produces repeatable arcing conditions. The initial focus is on tungsten as a plasma-facing material, including sintered tungsten, tungsten coatings on steel produced via physical vapor deposition (PVD), and functionally graded tungsten/steel coatings deposited by low-pressure plasma-spraying (LPPS). Thin PVD coatings (1-2 microns) fail catastrophically from an arc and result in severe delamination of the coating. The arc-induced damage of thicker coatings, such as those made via LPPS, tend to be restricted to the top few microns of the surface. Arcing often initiates on sharp surface microstructures and causes localized melting of tungsten at the surface of all the materials studied and results in resolidified melt pools with surface cracks. The resolidified surface results in a reduction in overall deuterium retention when exposed to typical RF plasma sheath conditions.

Caughman, John [ORNL] (ORCID:0000000206091164)↗
DOE OSTI2025

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↗
DOE OSTI2023

Field reversal in low pressure, unmagnetized radio frequency capacitively coupled argon plasma discharges

In general, the radio frequency (rf) electric field within a sheath points toward the metal electrode in low pressure, unmagnetized rf electropositive capacitively coupled plasma (CCP) glow discharges. This is due to the large ratio of electron to ion mobility and the formation of an ion sheath. In this work, we studied, using fully kinetic particle-in-cell simulations, a reversed electric field induced by the strong secondary electron emission during the phase of sheath collapse in a high-voltage rf-driven low pressure CCP glow discharge. We explored the transition behavior of the formation of field reversal as a function of driving voltage amplitude and found that field reversal starts to form at around 750 V, for a discharge with an electrode spacing of 4 cm at 10 mTorr argon pressure driven at 13.56 MHz. Accordingly, the energy distribution function of electrons incident on the electrode shows peaks from around 3 to 10 eV while varying the driving voltage from 150 to 2000 V, showing potentially beneficial effects in plasma material processing where relatively directional electrons are preferred to solely thermal diffusion electrons.

Physics↗
DOE OSTI2026

Direct Measurement of ICRF-Enhanced Plasma Potentials on WEST Using Reciprocating Emissive Probes

An extensive documentation of ICRF-enhanced plasma potentials has been conducted over two experimental campaigns on the WEST tokamak using reciprocating emissive probes magnetically connected to two ICRF antennas. The collected data spans a wide range of antenna electrical settings (coupled power, toroidal phasing, left–right power balance) and plasma parameters (density at the antenna limiter above and below the lower hybrid resonance, plasma current, minority fraction). By scanning the edge safety factor across multiple probe plunges, the magnetic connection between the probe and the antenna varied, enabling the construction of a 2D map of the plasma and floating potentials around an active ICRF antenna. This dataset will be used to validate RF simulation tools equipped with the sheath boundary condition and used to predict RF rectified potentials and ICRF-induced impurity sputtering in future machines. This paper presents the diagnostic and some initial measurements, while the rest will be reported elsewhere.

Diab, Raymond [Massachusetts Inst. of Technology (↗