Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “divertor geometry”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparison of private flux region instability in conventional and super-X divertor configurations

Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of lithium transport using the BOUT++ framework

A numerical model that calculates the collisional interactions between the lithium atoms from a lithium pellet and the background plasmas has been upgraded. The ion density ($N_t$), electron temperature ($T_e$), ion temperature ($T_i$) and parallel ion velocity ($V_{∥, i}$) are used to characterize the background plasmas. The lithium atom density ($N^{a}_{Li}$) and parallel velocity ($V_{∥,a}$) of lithium atoms evolve with time. For each lithium ion, the density ($N_{Li^{n+}}$), temperature ($T_{Li^{n+}}$) and parallel velocity ($V_{∥, Li^{n+}}$) are self-consistently calculated. A C-mod lower single null equilibrium is used to generate the grid for the BOUT++ simulation. The lithium atoms can be fully ionized to $Li^{3+}$ in ~2 μs. The rapid radial and poloidal expansion of the lithium ions are found in the simulation. After the collision interaction process, the electron temperature rapidly decreases at the pellet location; then, it rapidly poloidally expands, and the temperature at the pellet location starts to recover. The electron pressure increases at the pellet location despite the decrease in electron temperature because of the extra electrons from the lithium ionization. The ion pressure profile decreases in the pellet location due to the decrease in ion temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Statistical analysis of magnetic divertor configuration influence on H-mode transitions

DIII-D plasmas are compared for two upper divertor configurations: with the outer strike point on the small angle slot (SAS) divertor target and with the outer strike point on the horizontal divertor target (HT). Scanning the vertical distance between the magnetic null point and the divertor target over a range 0.10–0.16 m is shown to increase the threshold power, $P_\mathrm{th}$, and edge plasma power, $P_\mathrm{Loss},$ for the low-to-high confinement (L–H) and H–L transitions respectively, by up to a factor of 1.4. The X-point height scans were performed at three L-mode core plasma line average electron densities, $\bar{n}_\mathrm{e} = $ 1.2, 2.2 and 3.6 $\times 10^{19}\,\mathrm{m}^{-3}$, to investigate the density dependence of divertor magnetic configuration influence on $P_\mathrm{th}$. The X-point height, $Z_\textrm{x-pt}$, was further extended across the range 0.16–0.22 m with the more open HT divertor configuration, for which a clear decrease in $P_\mathrm{th}$ with increasing $Z_\textrm{x-pt}$ is observed. The dependence of $P_\textrm{th}$ on divertor magnetic geometry is further investigated using a time-dependent probability density function (PDF) model and information geometry to elucidate the roles played by pedestal plasma turbulence and perpendicular velocity flows. The degree of stochasticity of the plasma turbulence is observed to be sensitive to the plasma heating rate. The calculated square of the information rate shows changes in the relative density fluctuations and perpendicular velocity PDFs begin 2–5 ms prior to the L–H transition for three plasmas; providing a crucial measurement of the dynamic timescale of external transport barrier formation. Additionally, both information length and rate provide potential predictors of the L–H transition for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration

A comparative study of nitrogen versus neon has been carried out to analyze the impact of the two radiative species on power dissipation, SOL impurity distribution, divertor and pedestal characteristics. The experimental results show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. A significant amount of neon is measured in the plasma core leading to a steeper density gradient. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which for the first time at DIII-D includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. The impurity transport in the SOL is studied in terms of the parallel momentum balance showing that N is mostly retained in the divertor whereas Ne leaks out consistent with its higher ionization potential and longer mean free path. This is also in agreement with the enrichment factor calculations which indicate lower divertor enrichment for neon. The strong ionization source characterizing the SAS divertor causes a reversal of the main ions and impurity flows. The flow reversal together with plasma drifts and the effect of the thermal force contribute significantly in the shift of the impurity stagnation point affecting impurity leakage. This work provides a demonstration of the impurity leakage mechanism in a closed divertor structure and the consequent impact on pedestal. Since carbon is an intrinsic radiator at DIII-D, in this paper we have also demonstrated the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling. Here, carbon contributes more when neon seeding is injected compared to when nitrogen is used. Finally, the results highlight the importance of accompanying experimental studies with numerical modeling of plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry. In the cases presented here, plasma drifts and flow reversal caused by high level of closure in the slot upper divertor at DIII-D play an important role in the underlined mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improved core-edge compatibility using impurity seeding in the small angle slot (SAS) divertor at DIII-D

Impurity seeding studies in the slot divertor (SAS) at DIII-D have revealed a strong relationship between the detachment onset and pedestal characteristics with both target geometry and impurity species. N 2 seeding in the slot has led to the first simultaneous observation of detachment on the entire suite of boundary diagnostics viewing the SAS without degradation of core confinement. SOLPS-ITER simulations with D+C+N, full cross-field drifts and n-n collisions activated are performed for the first time in DIII-D to interpret the behavior. This illustrates a strong effect of divertor configuration and plasma drifts on the recycling source distribution with significant consequences on plasma flows. Flow reversal is found for both main ions and impurities affecting strongly the impurity transport and providing an explanation for the observed dependence on strike point location of the detachment onset and impurity leakage found in the experiments. Matched discharges with either nitrogen or neon injection show that while nitrogen does not significantly affect the pedestal, neon leads to increased pedestal pressure gradient and improved pedestal stability. Little nitrogen penetrates in the core, but a great amount of neon is found in the pedestal consistent with the different ionization potentials of the two impurities. This research demonstrates that neutral and impurity distributions in the divertor can be controlled through variations in strike point locations in a fixed baffle structure. Divertor geometry combined with impurity seeding enables mitigated divertor heat flux balancing core contamination and thus leading to enhanced divertor dissipation and improved core-edge compatibility which are essential for ITER and for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Divertor shaping with neutral baffling as a solution to the tokamak power exhaust challenge

Exhausting power from the hot fusion core to the plasma-facing components is one fusion energy’s biggest challenges. The MAST Upgrade tokamak uniquely integrates strong containment of neutrals within the exhaust area (divertor) with extreme divertor shaping capability. By systematically altering the divertor shape, this study shows the strongest evidence to date to our knowledge that long-legged divertors with a high magnetic field gradient (total flux expansion) deliver key power exhaust benefits without adversely impacting the hot fusion core. These benefits are already achieved with relatively modest geometry adjustments that are more feasible to integrate in reactor designs. Benefits include reduced target heat loads and improved access to, and stability of, a neutral gas buffer that ‘shields’ the target and enhances power exhaust (detachment). Analysis and model comparisons shows these benefits are obtained by combining multiple shaping aspects: long-legged divertors have expanded plasma-neutral interaction volume that drive reductions in particle and power loads, while total flux expansion enhances detachment access and stability. Containing the neutrals in the exhaust area with physical structures further augments these shaping benefits. These results demonstrate strategic variation in the divertor geometry and magnetic topology is a potential solution to one of fusion’s power exhaust challenge.

magnetically confined plasmas↗
DOE OSTI2022

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

Abstract MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Divertor characterization and access to dissipative divertor conditions in negative triangularity discharges in DIII-D

Reduced radial transport, short midplane-to-target parallel connection lengths, and a strong effect of cross-field drifts were responsible for the high densities required for detachment in strong negative triangularity (NT) configurations in DIII-D. Dissipative divertor conditions were achieved in NT discharges at different triangularity, injected power, plasma current, and toroidal field direction. Differences between negative and positive triangularity (PT) discharges are analyzed in this paper to understand the requirements for access to detached divertor conditions: power balance, geometry, radial transport and effect of cross-field particle drifts. Parametric dependencies of access to detachment on plasma current and power flowing into the scrape-off layer remained similar in negative and PT and impurity seeding was observed to reduce the density needed to detach by up to 30% at the expense of core impurity dilution. The impact of triangularity on core-edge integration was tested varying bottom triangularity at fixed top triangularity. The high density needed to detach was not intrinsic to the NT edge as shapes with positive lower triangularity and negative upper triangularity were able to detach at lower upstream densities while maintaining an ELM-free NT edge. Confinement degradation at deeper detachment levels was however observed in all NT shapes, often associated with radiation instabilities.

detachment↗
DOE OSTI2024

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

From MAST-U to STEP: Power exhaust control challenges and opportunities

This work explores the challenges and opportunities for power exhaust control in STEP, informed by dedicated MAST-U experiments. The STEP system primarily relies on hydrogenic and argon gas injection into the divertors, compensating for transients originating in the plasma core. Fast transients (>10 Hz) that evolve too quickly for actuators to compensate must be absorbed passively; the enhanced transient buffering of the foreseen long-legged divertor provides a clear advantage over conventional geometries, as demonstrated in the MAST-U Super-X divertor. Experiments further indicate that STEP’s tight divertor baffling enables near-independent control of the upper and lower divertors. The most challenging transients are expected from power-sharing fluctuations, as MAST-U experiments observed extremely fast dynamics. Fluctuations from core pellet fuelling remain relatively benign due to the small pellet size according to simulations. Exhaust control in a reactor like STEP requires an integrated approach befitting its machine-critical nature. A predictive control element, integrated with core plasma control, allows pre-emptive preparation of the divertor for incoming transients. An observer-based diagnostic approach is foreseen to monitor the divertor in the challenging reactor environment, supported by an extended diagnostic set in the non-nuclear phase to validate the required dynamic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neutral leakage, power dissipation and pedestal fueling in open vs closed divertors

The effect of divertor closure on neutral leakage, power dissipation and pedestal fueling is explored using the SOLPS code for different densities including attached and detached divertor conditions for DIII-D open and closed divertors. Here, a metric for divertor closure in terms of the percentage of neutrals escaping the divertor is presented. At the detachment onset the population of neutrals escaping the divertor is 12% for the closed divertor and 35% for the open divertor. This results in a lower contribution of the convective term to the total heat flux and increased power losses in the closed divertor. The comparison between the two geometries reveals that for the same upstream density, the upstream temperature is 23% higher in the closed divertor. The physical mechanisms for detachment at lower upstream density in the closed divertor are also identified with the reduction of the convective term directly resulting in 40% decrease and the increased power losses contributing with a factor of 2. The baffling also lowers the neutral flux refueling the pedestal resulting in further reduction of the core density for achieving detachment. The results introduced here show that the divertor recycling is the primary fueling mechanism in DIII-D. It is shown that fueling occurring in the region between the X-point and the outer midplane is at least four times lower in the closed divertor due the baffling structure which enhances target neutral trapping and consequently reduced upstream leakage. Due to the openness of the inner leg to the common flux and HFS regions both open and closed divertor experience almost identical fueling at the HFS highlighting the need to also consider the inboard divertor for shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Conceptual design of ELM control coils for the TCABR tokamak

An upgrade of the Tokamak à Chauffage Alfvén Brésilien (TCABR) is being designed to make it capable of creating a well controlled environment where the impact of resonant magnetic perturbation (RMP) fields on edge localised modes can be addressed over a wide range of (i) plasma shapes, (ii) divertor configurations, (iii) RMP coil geometries and (iv) perturbed magnetic field spectra. To this end, a unique set of in-vessel RMP coils is being designed and, in this work, their conceptual design is presented. This unique set of coils is composed of three toroidal arrays of coils on the low field side and three toroidal arrays of coils on the high field side. Each of these six toroidal arrays is composed of 18 coils thus allowing for the creation of RMP fields with toroidal mode numbers n ≤ 9 and with increased control of the poloidal mode number spectrum. To study dynamical effects of RMP fields of different toroidal mode numbers, all rotating simultaneously with different velocities, each of the 108 RMP coils will be powered independently by power supplies that can provide voltages of up to 4 kV and electric currents of up to 2 kA, with frequencies varying continuously from 0Hz up to 10kHz. A set of physical criteria were used to determine the optimal coil geometry and their respective number of turns to reduce the coil currents and voltages during operation with alternate current. Further, the conceptual design was carried out using both the vacuum approach (no plasma response) and the single-fluid response approach, which accounts for the response of a linear, single-fluid, visco-resistive plasma calculated using the M3D-C 1 code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Numerical analysis of gas exhaust in Wendelstein 7-X using the direct simulation Monte Carlo method

The present work is focused on a 3D numerical assessment of the Wendelstein 7-X (W7-X) particle exhaust. For all the numerical simulations the direct simulation Monte Carlo solver of the DIVGAS workflow, has been employed. The complex 3D geometry of the sub-divertor region includes the pumping gap panel, supporting structures, cooling pipes as well as the cryo-vacuum pump. All the considered flow simulations correspond to the Standard magnetic configuration of W7-X. The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; The coupling between EMC3-EIRENE and DIVGAS, which considers the fact that the incoming neutral particle flux at the sub-divertor is based on realistic plasma background, has been demonstrated. Three plasma scenarios have been considered, for which is clearly seen that by increasing the heating power, the neutral pressure as well as the resulting pumping efficiency is increased. The obtained numerical results of the neutral pressure in the sub-divertor lie within a more general scan matrix, which assumes a wider range of incoming particle flux, namely 10 19 –10 24 (s −1 ). It has been observed that, the sub-divertor neutral pressure is proportional to the incoming neutral particle flux, with the effective pumping speed to be a constant of proportionality. The influence of switching off the cryo-vacuum pump on the sub-divertor pressure is rather modest and a weak increase of the neutral pressure in the sub-divertor is expected. Correlations of the sub-divertor pressure with the total incoming particle flux as well as the individual pumped flux at each of the AEH and AEP sections have been deduced. Moreover, it has been demonstrated that the influence of the incoming neutral particle flux on the albedo coefficient at the AEH and AEP pumping gaps is rather weak. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of future experimental campaigns.

DSMC method↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2025

Recrystallization, cracking, and erosion of dispersoid-strengthened tungsten materials during exposure to divertor plasmas

In this study, we investigated the effects of combined intense particle and heat flux exposure on advanced tungsten plasma-facing materials within the DIII-D fusion facility. Our test matrix included two types of dispersoid-strengthened tungsten (containing either 100 nm diameter TiO 2 or Ni particles), along with high-purity polycrystalline tungsten as a reference. This experiment relied on a sample geometry angled at 15° relative to the divertor surface, thereby allowing the surfaces to intercept steady-state perpendicular heat fluxes (q ⟂ ) ranging from 10.1 to 19.6 MW/m 2 . During each shot, the samples were exposed to 42 Hz edge-localized modes (ELMs), allowing us to test the material response to transient heating. We correlated the exposure conditions with extensive post-test surface composition analysis and microscopy to determine how the plasma modified each surface. The angled specimens closest to the strike point received the highest combined heat and particle flux and melted midway through the experiment. EBSD analysis revealed they were completely recrystallized throughout, with an average grain size >100 µm. On the other hand, the specimens that received a lower steady state heat flux survived with more superficial surface damage. Whereas the high-purity polycrystalline tungsten exhibited a higher surface roughness, the dispersoid-strengthened material exhibited more extensive shallow inter-granular cracking. In addition, the surface was depleted of dispersoids following plasma exposure, possibly because of evaporation and/or sputtering. The results described here provide insights into the performance of these materials in a fusion environment which can guide further optimization for use in long-pulse devices.

Kolasinski, Robert D. [Sandia National Laboratorie↗
DOE OSTI2023

Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency

Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.

ICRH antenna↗