Tungsten Erosion Research Articles

Simulations of tungsten fuzz growth and erosion under He/Ar mixed plasma irradiation on LP-MIES

Extrinsic noble impurities seeding has been commonly used to reduce the heat flux on divertor target. However, the physical sputtering yield on the target could be obviously enhanced due to larger mass of noble impurities. On the other hand, it has been reported that tungsten fuzz can suppress the physical sputtering due to the trapping effect on sputtered particles. Hence, it is interesting to check the erosion behaviors of tungsten fuzz under the irradiation of noble impurities. Recently, experiments of He/Ar mixed plasma irradiation on W target were conducted on LP-MIES device. Dedicated simulations have been performed by SURO-FUZZ code to study the impact of Ar flux on fuzz growth and erosion, which can well reproduce the experimental measurements on LP-MIES. For the case of pure He incidence, the growth of porous nanostructure leads to an obvious reduction in the fraction of the implanted He particles, which stabilizes to 15 % when the fuzz layer thickness exceeds 1500 nm. For He/Ar mixed plasma irradiation, the dependence of the position of the net sputtering on the local porosity has been investigated in detail. It is found that the net physical sputtering mainly occurs in the regions with the local porosity larger than 0.6.

Survey of tungsten gross erosion from main plasma facing components in WEST during a L-mode high fluence campaign

An initial high fluence campaign was performed in WEST, in 2023, on the newly installed actively cooled tungsten divertor composed of ITER-grade monoblocks. The campaign consisted in the repetition of a 60 s long Deuterium L-mode pulse in attached divertor conditions, cumulating over 10000s of plasma exposure. A maximum deuterium fluence of approximately 5⋅1026m−2 was reached in the outer strike point region, representative of a few high performance ITER pulses. Gross tungsten erosion inferred from visible spectroscopy shows that the most eroded plasma facing component is the inner divertor target with rates ten times larger than on the outer divertor target. The outer midplane tungsten bumpers, located a few centimeters from the plasma, show gross erosion rates two times lower than at the outer divertor. We conclude that the outer midplane bumpers have a negligible contribution to the long range tungsten migration and deposition onto the lower divertor. The cumulated gross erosion rate on the inner divertor translates in an effective gross erosion thickness of about 20μm, while it is about 2μm for the outer divertor. Strikingly, these orderings coincide with the thickness of deposits found locally on the divertor : the exposed surfaces of high field side monoblocks are covered with several tens of μm tungsten deposits, while on the lower field side, few μm thin tungsten deposits are only found on the magnetically shadowed parts of monoblocks. The strong impact of those deposits on WEST operation, namely perturbation of surface temperature measurement with infra-red thermography, and the emission of flakes causing radiative perturbation of the confined plasma, calls for anticipating similar issues in ITER. In particular, the start of research operation shall consider the definition of a divertor erosion budget in order to anticipate the formation of deleterious deposits.

An integral approach to plasma-wall interaction modelling for EU-DEMO

Abstract The integral approach to plasma-wall interaction (PWI) modelling for DEMO is presented, which is a part of EUROfusion Theory and Advanced Simulation Coordination (E-TASC) activities that were established to advance the understanding and predictive capabilities for modelling of existing and future fusion devices using a modern advanced computing approach. In view of DEMO design, the aim of PWI modelling activities is to assess safety-relevant information regarding erosion of plasma-facing components (PFC), including its impact on plasma contamination, dust production, fuel inventory, and material response to transient events. This is achieved using a set of powerful and validated computer codes that deal with particular PWI aspects and interact with each other by means of relevant data exchange. Steady state erosion of tungsten PFC and subsequent transport and re-deposition of eroded material are simulated with the ERO2.0 code using a DEMO plasma background produced by dedicated SOLPS-ITER simulations. Dust transport simulations in steady state plasma also rely on the respective SOLPS-ITER solution and are performed with the MIGRAINe code. On the way to improved simulations of tungsten erosion in the divertor of DEMO, relevant high density sheath models are being developed based on Particle-in-Cell (PIC) simulations with the state-of-the-art BIT code family. PIC codes of the SPICE code family, in turn, provide relevant information on multi-emissive sheath physics, such as semi-empirical scaling laws for field-assisted thermionic emission. These scalings laws are essential for simulations of material melting under transient heat loads that are performed with the recently developed MEMENTO code, the successor of MEMOS-U. Fuel retention simulations assess tritium retention in tungsten and structural materials and permeation to coolant accounting for neutron damage, in particular for divertor monoblocks of different sizes, using the FESTIM code, and for the first wall, in particular in view of tritium self-sufficiency, using the TESSIM code. Respective code-code dependencies and interactions, as well as modelling results achieved to date are discussed in the contribution.

Interaction of iron melt with tungsten and WFe composite structure evolution

Development of new plasma facing components (PFC) draw attention of numerous scientific groups in fusion energy field. Tungsten as a main PFC material main have poor machinability, thus various designs were proposed to overcome this problem. However, such technologies as functional graded layers, sintering and additive technologies are limited in production size and have low cost-efficiency. This research considers the alternative approach of steel melt injection or melt infiltration of tungsten mesh. The results obtained establish the mechanism of phase evolution during interaction of solid W with liquid steel and iron. Energy dispersive spectroscopy, X-ray diffraction and Thermocalc calculations used to analyze the phase composition of tungsten wetted with various steels and pure iron. Results show that interaction rate significantly depends on melt temperature and overheating above its liquidus. Then the overheating exceeds 150 °C erosion of tungsten occur.

Validated edge and core predictions of tungsten erosion and transport in JET ELMy H-mode plasmas

Predictive edge and core simulations of tungsten (W) erosion and transport in JET ITER-like wall plasmas are shown to be consistent with the experimentally inferred W density in the main plasma, within the uncertainty inherited from the measurements of the deuterium plasma conditions and from the W density measurements. The ERO2.0 code is applied to predicting the W erosion and edge transport, whereas JINTRAC predicts W transport from the pedestal top to the core plasma. The studied plasma scenarios range from L-mode to the highest-performance deuterium ELMy H-mode in JET.

Impact of helium and hydrogen plasma exposure on surface damage and erosion of tungsten

The impact of helium plasma exposure on the tungsten surface damage structure development and erosion has been investigated by comparing the impact of hydrogen plasma exposure. Crystal orientation dependence of the undulating surface structure formation and erosion rate is observed on the plasma-exposed tungsten surface independently from the plasma species. The top surface of the plasma exposed tungsten has a tendency to {100} plane independently from the initial surface orientation. Although hydrogen and/or helium cause no erosion in tungsten under incident ion energy exposure conditions below the sputtering threshold, inevitable minute impurities, like oxygen, play an essential role in erosion, and significant erosion can be observed even at 30 eV.

Modeling of tungsten divertor target erosion induced by impurity during edge‐localized modes by using a kinetic model

AbstractThe burst of edge‐localized modes (ELMs) leads to an increase in the energy and particle fluxes to the divertor target. Tungsten (W) is chosen as the primary candidate material for plasma‐facing components (PFCs) in the future fusion devices, and EAST has already upgraded all divertors to use W. Therefore, understanding tungsten target erosion during ELMs and finding the correlation between erosion rate and key ELM parameters are crucial for steady‐state operation. In this work, based on the Vlasov–Poisson model (VPM), we develop a one‐dimensional kinetic parallel transport code to investigate the parallel transport of particles in the EAST device during ELMs and the resulting target erosion. The EAST experiment (#102182) is simulated by VPM code. The simulation results are compared with experimental data as well as free‐stream model (FSM) calculation, showing the accuracy of the code. Considering the presence of lithium (Li) impurities in EAST discharge, the erosion of the W target is simulated. The results indicate that during the burst of ELM, the total average tungsten erosion rate, , is determined by both deuterium (D) and Li ions. D ions dominate the erosion when the ELM frequency (fELM) is low (ranging from 50 to 175 Hz), while Li impurities become more important than D+ in high‐frequency ELMs (fELM > 175 Hz). As fELM increases, the time‐averaged erosion of the W target first increases and then decreases. Therefore, the reduction of W erosion benefits from high‐frequency ELMs, with impurity ions being the primary contributor to the erosion.

Recent DIII-D progress toward validating models of tungsten erosion, re-deposition, and migration for application to next-step fusion devices

Fundamental mechanisms governing the erosion and prompt re-deposition of tungsten impurities in tokamak divertors are identified and analyzed to inform the lifetime of tungsten plasma-facing components in ITER and other future devices. Various experiments conducted at DIII-D to benchmark predictive models are presented, leveraging the DiMES removable sample exposure probe capability and the Metal Rings Campaign, in which toroidally symmetric rows of tungsten-coated tiles were installed in the DIII-D divertor. In tokamak divertors, the width of the electric sheath is of the order of the main ion Larmor radius, and a vast majority of sputtered tungsten impurities are typically ionized within the sheath. Therefore, W prompt redeposition is mainly governed by the ratio of the characteristic ionization mean-free path of neutral tungsten to the width of the sheath. In-situ monitoring of the prompt redeposition of tungsten impurities in divertors is demonstrated via the use of WII/WI line ratios and the ionizations/photon (S/XB) method in L-mode discharges. Even with this relatively limited set of emission measurements, net erosion measurements were found to be a consistent upper bound to an analytic scaling based on the ratio of the W ionization length, and the width of the magnetic sheath rather than the ratio of and the W+ gyro-radius. In the far-scrape-off layer (SOL) of the ITER divertor, however, it is calculated that the measurement of photon emissions associated with the ionization of tungsten impurities up to may be required. Finally, W deposition patterns on DiMES collector probes, interpreted via DIVIMP-WallDYN modelling, reveal the key roles of progressive W erosion/re-deposition staps and E × B drifts in regulating long-range high-Z material migration.

Erosion estimates for the divertor and main wall components from STEP

The tungsten erosion within Spherical Tokamak for Energy Production (STEP) assuming tungsten main wall and tungsten divertor has been estimated with ERO at the inner and outer divertor, at the inner and outer midplane and at the outboard baffle entrance. Plasma parameters are based on SOLPS simulations applying argon puffing for edge cooling. The plasma parameter range covers peak electron temperatures T e between 3 and 25 eV in the divertor. At the inner midplane T e ∼ 13 eV, at the outer midplane ∼7 eV and at the outboard baffle entrance between 1 eV and 4 eV. The modelled peak gross erosion is highest in the divertor with up to 1E19 W m−2 s−1 within the inner and 7E19 W m−2 s−1 in the outer one for the plasma parameter range studied. At the main wall the gross erosion is about 2E18 W m−2 s−1 at the inner midplane and 1.3E17 W m−2 s−1 at the outer one. However, tungsten deposition within the divertor is much larger with amounts between 88% and 98% and only between 10% and 60% at the midplane. At all locations studied, tungsten erosion due to deuterium ions is negligibly small compared to the erosion by argon ions. Erosion due to deuterium atoms has been studied for the outer midplane and is there at least four times smaller than the erosion due to argon ions. The simulations have been performed considering singly ionised Ar. However, according to the SOLPS runs the mean charge of Ar impinging the surfaces is about two at the locations of largest erosion, which leads to an increase of the gross erosion by a factor between 1.5 and 5 with the largest increase occurring at the outer divertor target.

Model validation of tungsten erosion and redeposition properties using biased tungsten samples on DiMES

An experiment was performed in the DIII-D lower divertor to validate numerical SOL tungsten (W) impurity erosion and redeposition simulations against experimental data. The net and gross erosion of W were calculated as a function of the voltage (or bias) applied to the exposed material. Five samples were inserted into the DIII-D lower divertor using the Divertor Material Evaluation System (DiMES) manipulator and exposed to constant L-mode attached plasma conditions. Each sample was partially coated with W. During plasma shots, samples were biased with respect to the machine vessel ground, ranging from −60 V to 25 V. The ERO2.0 code was used to numerically simulate the experiment aiming to compare the numerical results with experimental measures. A good agreement is found between estimated and measured tungsten erosion at least for negative biases.

Modelling of physical sputtering properties during tungsten fuzz growth under various impact energies on NAGDIS-II

Fibre-like nanostructure named fuzz can form on tungsten (W) surface under the irradiation of helium (He) plasma. While the W fuzz growth can be suppressed by the bombardment of energetic He particles. The experiments of tungsten fuzz growth and erosion exposed by He plasma with impact energies in the range of 200 ∼ 500 eV have been conducted on the linear divertor simulator NAGDIS-II. To study W physical sputtering properties during fuzz growth on NAGDIS-II, three-dimensional kinetic Monte Carlo code SURO-FUZZ has been upgraded by including the migration of W adatoms and the deposition of sputtered W particles. By means of the renewed SURO-FUZZ, dedicated simulations of W fuzz growth and erosion under various He impact energies have been performed to reproduce the NAGDIS-II experiments. The simulated fuzz layer thicknesses under erosive regimes show a good agreement with the measurements on NAGDIS-II. The W net physical sputtering yield on W fuzzy surfaces increases with He impact energy according to SURO-FUZZ modelling. Further, the relative W net physical sputtering yields of the fuzzy surface to the smooth surface (Yfuzz/Ysmooth) have been studied under different He impact energies, which reveals that the decay trend of Yfuzz/Ysmooth is mainly dependant on the fuzz layer thickness.

Morphological damage to tungsten eroded by abrasive water jet

Abrasive water jet is a promising method for processing rolled tungsten plates, which are considered the most promising plasma-facing material in fusion devices. To investigate the erosion of tungsten due to the abrasive water jet, we conducted erosion tests and numerical simulations. We examined the morphological characteristics of the tungsten erosion damage. An erosion model was developed by coupling smoothed-particle hydrodynamics with the finite-element method to simulate tungsten erosion. The mechanism behind the morphological damage characteristics of tungsten was analyzed. Our findings reveal that as erosion depth increases, the kinetic energy of the abrasive particles decreases while the erosion angle increases. The erosion damage to tungsten transitions from brittle fracture induced by microcutting and microplowing to impact deformation. As a result, the erosion section's morphology deteriorates, exhibiting strips from the cutting process, increasing both the trailing angle and surface roughness. The morphological damage characteristics of tungsten are intricately linked to the abrasive particles' motion. Therefore, by altering the motion characteristics of these particles, the erosion section's quality can be enhanced. This research offers valuable insights for enhancing the quality of tungsten sections processed by abrasive water jets in fusion devices.

Tungsten erosion and divertor leakage from the DIII-D SAS-VW tungsten-coated divertor in experiments with neon gas seeding

Collector probes have been used to examine tungsten divertor leakage in a variety of scenarios with low-Z impurity seeding during operation with the new tungsten-coated SAS-VW divertor in DIII-D. Measurements of tungsten deposition on collector probes inserted into the far Scrape-off-Layer (SOL) are used to deduce how efficiently tungsten leaks out of the closed, V-shaped divertor after it is eroded from the target surfaces. Qualitative differences in the tungsten deposition patterns across the collector probes provide clear experimental evidence that the SOL conditions depend on the low-Z impurity seeding conditions. These measurements show that in scenarios where neon gas is injected into the plasma, the tungsten divertor leakage and SOL transport depend on the poloidal location from which the neon is injected. In particular, neon injection from the Inner Midplane and Outer Midplane appear to each result in higher divertor leakage by a factor of 2 to 3 compared to cases with neon injection from either the SOL Crown or from the SAS-VW divertor itself.

Time-dependent collisional radiative modeling of tungsten in the magnetic sheath for erosion diagnosis

Abstract Tungsten is the material of choice for the divertors in ITER, SPARC and future fusion reactors. Accurate diagnosis of tungsten erosion and migration is important for first wall life time, slag production and core performance. The addition of a magnetic presheath requires time-dependent collisional radiative effects to be included for accurate neutral tungsten collisional radiative modeling. Gross erosion measurements could be modified by a factor of 10 due to the inclusion of time-dependent effects for ITER relevant divertor conditions. A simple sputtering model and sheath density model are developed to investigate time-dependent collisional radiative effects. Neutral tungsten spectral lines populated from different metastable levels depend on model parameters leading to potential spectroscopic diagnostics of plasma parameters. Electron temperatures inferred from spectroscopic line ratios are in agreement with Langmuir probe measurements in the Compact Toroidal Hybrid.

Study of the erosion and redeposition of W considering the kinetic energy distribution of incident ions through a semi-analytical model

In today’s nuclear fusion devices, erosion of high-Z metallic plasma-facing materials (PFMs) is mainly caused by physical sputtering. That is, by the exchange of energy between plasma ions and the atoms in the walls. In most of the numerical codes currently in use impinging plasma is approximated as a fluid. By averaging the incident particles’ energy distribution the high-energy population of the eroded material is underestimated. For heavy materials such as W, high-energy eroded particles tend to ionize far from the wall and they are less affected by the sheath electric field hence, not being attracted back to the wall, they have a higher chance to contaminate the core plasma. This could in turn result in an underestimation of the net erosion sources. In this work, a semi-analytical model was developed to include the energy distribution of the incident particles. Then, by Monte Carlo method, the net erosion of tungsten from a smooth PFM was calculated. The results show that the kinetic description in energy is important only for incident particles ionized once. For instance, it is particularly important for plasma ions such as Deuterium. It is seen that Deuterium contribution to the W net sources is not always negligible if compared to light impurities or to tungsten self-sputtering in the range of plasma parameters tested. Finally, results show that the difference between the fluid and kinetic models becomes more pronounced for high-screening plasma conditions.

KTIG–MIG hybrid arc welding process

Metal-inert gas welding is an important welding process but the process instability easily induces spatter and undercut defects. To improve the welding process stability, a new hybrid arc welding process, named KTIG–MIG, was proposed in this paper. Arc characteristics were preliminarily studied. It is found that in KTIG–MIG hybrid welding, MIG arc is stable even it is shielded with pure argon gas; stable spray transfer is easily achieved in a much wider process window, especially the lowest voltage to achieve spray transfer is greatly reduced; KTIG arc voltage decreases than that in single KTIG arc, and MIG arc current is higher than the pre-set value. Cooling state of the KTIG tungsten plays a key role in slowing tungsten erosion.

Absence of synergistic effects in quasi-simultaneous sputtering of tungsten by Ar and D ions

A quartz crystal microbalance was used to experimentally study the erosion of tungsten during rapidly alternating bombardment with 2 keV argon and deuterium projectiles. A key goal was to investigate whether the mean sputtering yield of the alternating irradiation can be predicted from data for sputtering yields of single ion species. In addition, influences by residual gas pressure in the UHV experiment and variable ion fluxes have been studied. Our results show that the mean sputtering yield of irradiations with alternating ion species can be well predicted for a range of different fluence ratios as a simple superposition of individual sputtering yields, weighted by the respective relative fluences. This finding supports that no synergistic sputtering effects were relevant in the investigated low-flux regime.

Multiphysics Modeling of Impurity Transport for FNSF Startup Scenario with ERO2.0

This study focuses on performing a multiphysics study using the ERO2.0 and UEDGE codes for two standard double null configurations for the Fusion Nuclear Science Facility: (a) 100% recycling and (b) 99% recycling. Results show that the main contributor to tungsten erosion along the divertor plates is impurities from the midplane waveguides. In addition, the standard high-recycling case (100% recycling) shows a significantly higher buildup of impurities along the divertor tiles during the startup phase, which can lead to a higher increase of energy loss in the plasma during steady-state operation. Last, for high recycling, anomalous diffusion can dominate over parallel field diffusion. The work performed in this study can be iteratively applied to a full operation scenario with additional physics such as those from neutrals, wall shaping, and additional external fields.

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C+6 and D+1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM’ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. The ELMs result in a ∼25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

Erosion of Fuzz Layers Formed in Steady-State Plasma Discharge

The erosion of nanostructured tungsten and titanium by high-heat plasma flux, laser, and arcing is investigated. To fabricate nanostructural fuzz layers and hierarchical granularity on the surfaces, samples were exposed to helium plasma in the steady-state plasma device PLM-M, which is a linear plasma trap of an eight-pole multicusp magnetic field with parameters similar to the scrape-off layer and divertor plasma in a tokamak. Arcing ignited with a Nd:YAG laser pulse on the target fuzzy surface in the helium plasma resulted in the melting of fibers and the creation of craters of several microns in depth and several tens of microns in diameter.