Search NASA⌕ Search

SEARCH · Search NASA

Results for “tungsten erosion”

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

Multi-physics modeling of tungsten collector probe samples during the WEST C4 He campaign

Abstract We describe the results of a multi-scale, multi-physics modeling assessment of SOLPS-ITER, hPIC2, RustBCA and Xolotl, in which five single-crystal tungsten (W) samples were placed in a reciprocating collector probe and exposed to helium (He) plasma in the WEST fusion device. In our models, we considered a pure (100 %) He plasma, as well as one with oxygen (O) present (95% He 5% O) corresponding to the impurity concentration estimated during the C4 He campaign in WEST. Our SOLPS simulations approximately match experimental reciprocating Langmuir probe plasma measurements of plasma density and temperature. Using these plasma parameters as input, hPIC2 and RustBCA predict that the presence of oxygen impurities lead to a 15%–20% decrease in ion and heat fluxes to the surface, and an order of magnitude higher sputtering yields (compared with a pure He plasma). Xolotl predictions for the response of tungsten to plasma surface interactions (PSIs) agree with experimental LAMS analysis, and indicate large near-surface He concentrations, which quickly decay with depth. Our model also shows an increasing role of erosion—in removing the near-surface He—with time. Overall, slightly higher retention is predicted for tungsten exposed to a pure He plasma, with the largest differences in the near-surface gas content caused by the large oxygen-induced erosion. This highlights the important role that impurities play in PSI. Therefore, future work will focus on providing a fully self-consistent description of oxygen (and oxides, etc.) in our models, through multi-species implementation in GITR and inclusion of oxygen and tungsten oxide formation in Xolotl.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

This project focused on the development and integration of high-performance simulation tools to predict the operating behavior of Plasma-Facing Components in magnetic confinement fusion systems. A key objective at Illinois was to assess the impact of the dynamic interplay between the evolving material surface and the magnetized plasma sheath, and characterize the impact of tungsten-based PFCs on plasma contamination, including phenomena such as surface erosion, dynamic recycling of fuel species, and tritium retention, which are critical for the success of future magnetic fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

Abstract All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The dependence of tungsten fuzz layer thickness and porosity on tungsten deposition rate and helium ion fluence

Abstract Fuzz formation on a heated tungsten surface in the presence of a helium-containing plasma and tungsten deposition source was investigated. Tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, W/He ion flux ratio of ∼ 0.4 × 10 − 4 , and varied helium ion fluence from 0.18 to 3.4 × 10 26 m −2 . Fuzz thickness was measured by cross-sectional scanning electron microscopy to increase from 0.22 to 15 µ m with increasing helium ion fluence. No indication of saturation in fuzz thickness at high fluence was observed, in contrast to fuzz produced on a tungsten surface without tungsten deposition. Additional tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, helium ion fluence of ∼ 3.4 × 10 26 m −2 , and varied W/He ion flux ratio from 0.26 to 3.0 × 10 − 4 . Fuzz thickness increased from 7.5 to 120 µ m with increasing W/He ion ratio. A final sample exposed at 1123 K to a mixed helium-deuterium plasma with ion incident energy of 76 eV, helium ion fluence of 0.18 × 10 26 m −2 , and W/He ion flux ratio of 2.2 × 10 − 4 developed nearly identical fuzz structures to that developed in a pure He plasma. As a function of deposited tungsten fluence, all results were found to trace out a single layer-growth curve given by a power law relation, indicating that fuzz thickness is independent of the W/He ion flux ratio in the range investigated and independent of any deuterium present in the plasma. As a result, for tungsten plasma facing walls in magnetic fusion devices at 1000–2000 K with 10 −4 W/He ion flux ratio, fuzz with thicknesses greater than hundreds of microns may form in as little as 10 4 s (in the absence of ELM-induced erosion or annealing), and may more significantly affect its thermophysical properties than fuzz generated without a tungsten deposition source.

Physics↗

Electrode Erosion and Prefire Studies Towards Fusion Scale Pulsed Power

This study presents a comprehensive investigation of electrode erosion and discharge behavior in spark gap switches over long switching cycle lifetimes. Brass, copper–tungsten (CuW), and stainless steel electrodes are tested under controlled conditions to quantify material degradation, debris accumulation, and changes in breakdown voltage. High-resolution imaging and statistical analysis of spark channel locations and gap breakdown voltages reveal how surface evolution influences long-term performance and reliability. These results provide essential data for lifetime modeling and inform design strategies for pulsed power systems in emerging applications such as private sector fusion energy and large-scale facilities like Sandia’s Z Machine and proposed ZX upgrades, where high repetition reliability and predictable behavior are critical.

electrical breakdown↗

Design and physics basis for the upcoming DIII-D SAS-VW campaign to quantify tungsten leakage and transport in a new slot divertor geometry

A set of experiments are planned to exploit the high SOL collisionality enabled by a tightly baffled slot divertor geometry to suppress tungsten leakage in DIII-D. A toroidal row of graphite tiles from the Small Angle Slot (SAS) divertor is being coated with 10–15 μ m of tungsten. New spectroscopic viewing chords with in-vacuo optics will measure the W gross erosion source from the divertor surface with high spatial and temporal resolution. In parallel, the bottom of the SAS divertor is changed from a flat to a ‘V’ shape. New SOLPS-ITER/DIVIMP simulations conducted with drifts using the planned ‘V’ shape predict a substantial reduction in W sourcing and SOL accumulation in either B × ∇B direction relative to either the old SAS divertor shape or the open, lower divertor. Dedicated studies are planned to carefully characterize the level of W sourcing, leakage, and scrape-off-layer (SOL) accumulation in DIII-D over a wide range of plasma scenarios. Various actuators will be assessed for their efficacy in further reducing high-Z impurity sources and leakage from the slot divertor geometry. Finally, this coupled code-experiment validation effort will be used to stress-test physics models and build confidence in extrapolations to advanced, high-Z divertor geometries for next-step devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liquid lithium wetting and percolation in a porous tungsten/liquid Li plasma facing component (PFC)

This work explores the concept of a hybrid material system designed to protect the plasma from high-Z material emission while also protecting an underlying high-Z wall material from erosion by integrating a low-Z component in the liquid phase. Plasma-surface interaction properties of a porous tungsten (W)-liquid metal hybrid system, having the favorable bulk thermomechanical properties of W while serving as a scaffold for a liquid metal with self-healing and radiative vapor shielding characteristics, are examined. W-substrates with 70% density of bulk W and 1–5 μm sized pores have been fabricated with 50-nm W powders using spark plasma sintering. Enhanced lithium (Li) wettability, driven by percolation through the porous W architecture, is demonstrated with in-situ liquid Li drop measurements. Finally, results show complete wetting of liquid Li at 250 °C in the porous W, 100 °C lower temperature than what has been observed on traditional W surfaces. In operando ion beam analysis of Li films deposited on porous W substrates provides microscopic evidence of Li percolation into the substrate and a first look at D retention up to a depth of 100 nm after plasma exposure at a flux of 10 20 m –2 s –1 , 250 eV/amu D + plasma to a fluence of 1.16 × 10 23 m –2 at a temperature of 200 °C.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beryllium-driven structural evolution at the divertor surface

Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational exploration of DIII-D QH-mode with extended MHD

This work will enable improved understanding of the scientific basis for DIII-D Quiescent H-mode (QH-mode) by leveraging recent computational developments of the NIMROD code. Experimentally, it is established that the QH-mode tokamak operation regime is free of edge-localized modes (ELMs) with good energy confinement. ELMs produce an impulsive heat-flux on the divertor that is unacceptable under burning-plasma conditions. Beyond DIII-D in future high-power-density tokamaks, accommodation of a high heat flux at the divertor represents a key challenge to overcome. Proposed solutions involve operation at high divertor density to achieve divertor detachment, operation with a liquid lithium wall and/or introduction of impurities to enhance radiative losses. Additionally, the first wall material is likely to be a high-Z material such as Tungsten in future devices as materials such as graphite tiles restrict the lifetime of the first wall due to a high erosion rate. While the atomic species composition of the plasma is known to impact the edge-pedestal structure and stability, the specifics of the underlying physics basis for this interaction is not firm. Studying this interaction is the goal of this proposal without focusing specifically on the fine details of the first wall and divertor solution. It is the goal is to lay the ground work to establish the compatibility of the divertor solutions and wall material with pedestal solutions that eliminate ELMs by leveraging present DIII-D experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control of the Plasma-Material Interface for Long Pulse Optimization in EAST and KSTAR

This project's major goal was to understand and control the plasma-material interface to improve long pulse discharge and performance in the EAST and KSTAR devices, emphasizing the unique strengths of each machine. In particular, MIT has developed implanted depth markers to diagnose the net erosion of high-Z materials. The goals included the characterization and utilization of these markers in various high-Z materials of interest, including TZM, tungsten, and refractory metal base high entropy alloys.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the effect of B T direction on W source-to-core pathways during the SAS-VW campaign on DIII-D

Experiments using the V-shaped closed slot tungsten (W) coated SAS-VW divertor in DIII-D studied the effects of the BT direction on core contamination of eroded tungsten from a closed slot divertor configuration. Core W content is inferred using soft-X ray tomography (SXR) and vacuum ultraviolet spectroscopy (SPRED), while W divertor erosion is inferred from visible spectroscopy of W emission (400.9 nm) measured by in-slot filterscopes (filtered photo-multipliers). Post-mortem analysis from the campaign discovered tile misalignment leading to suspected pronounced leading-edge erosion in the unfavorable BT direction (ion B ⃗×∇B ⃗ drift away from divertor) likely not captured by diagnostics. However, empirical findings show up to ~2-3x larger core contamination in the favorable BT direction even considering no additional W erosion from leading edges. A “source-to-core efficiency factor” is derived to estimate the effects of leading-edge erosion and compare W contamination for two pairs of H-mode discharges in opposite BT directions. While having differing absolute parameters, similar core impurity density gradients suggest comparable core impurity transport. These results show that favorable BT may have stronger source-to-core pathways for W impurities sourced from the outer divertor region. Possible explanations could include the effects of E ⃗×B ⃗ drifts on W transport in the scrape-off-layer (SOL) as well as previously determined fast SOL inner target directed flows in favorable BT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor Plasma Detachment Control Neural Network

DivControlNN is a state-of-the-art software tool that leverages advanced machine learning techniques to predict and control divertor plasma behavior in fusion reactors. Plasma, a highly energetic and electrically charged gas, requires meticulous management to protect reactor components and maintain optimal energy production. Conventional simulation methods, although extremely detailed, typically demand extensive computational time-making them unsuitable for real-time control scenarios. DivControlNN addresses this challenge by learning from tens of thousands of high-fidelity simulations, thereby creating a rapid surrogate model that can deliver near-instantaneous predictions. At the core of its functionality is a sophisticated technique known as latent space mapping, which condenses complex, high-dimensional plasma data into a compact, lower-dimensional representation. This streamlined representation enables the system to quickly forecast essential plasma properties and determine the precise conditions required for effective detachment. Detachment is a crucial process in which the plasma is cooled before reaching the divertor plates, thereby reducing heat loads and mitigating material erosion. In recent experiments conducted on the KSTAR tokamak in South Korea, DivControlNN successfully guided the detachment process without any fine-tuning-even when applied to a new tungsten divertor configuration. By achieving a computational speed-up of over one hundred million times compared to traditional simulation methods while maintaining low prediction errors, DivControlNN stands to significantly enhance real-time control and diagnostic capabilities in future fusion reactors. This breakthrough paves the way for safer, more reliable reactor operation and represents a major advancement toward realizing fusion energy as a practical, sustainable, and clean power source.

Xu, Xueqiao [Lawrence Livermore National Laborator↗

Synergistic effects of W–Ta alloys exposed to various ELM-like fusion events

Abstract Recently tungsten alloys have shown capabilities in suppressing morphological changes when exposed to harsh fusion conditions. Morphological changes due to He + and D + ions at high temperatures (1000–2000 K) can lead to easier erosion and sputtering of tungten (W) particles towards quenching the fusion plasma. This study investigated W–3Ta alloy as potential plasma facing materials (PFMs) exposed to different edge-localized mode (ELM)-like energies and ion loadings. Performance of W–3Ta has been tested as a function of ELM-like energies (0.6–1.5 GW m −2 ) and varying He + and D + ion concentrations (100% He + , 50% He + :50% D + ,10% He + :90% D + (similar ELM-like reactor condition), and 100% D + ions). Early-stage fuzz formation and surface pores were observed during 100% He + ion irradiation and mixed ion loadings but not observed during 100% D + ion loadings. Pores decreased in quantity and increased in size with increasing transient heat loading. Grain boundaries acted as trapping sites of helium and deuterium atoms during dual beam ion irradiations, as observed through increased pore size along the grain boundary. Qualitative erosion estimates were measured using witness plates. Increased erosion was observed with increasing D + ion concentration, which is attributed to D supersaturation in the near surface region affecting thermal and mechanical properties. More erosion was observed for W–3Ta samples compared to pure W, evident through large (>10 µ m) cracks on the surface, and an increase in W metal x-ray photoelectron spectroscopy peaks under similar loading conditions. The slight enhancement in erosion is attributed to the slightly decreased thermal properties, ∼10% thermal diffusivity drop, of W–3Ta at high temperatures. The damage threshold was impacted by ion loading, where distinct cracking was observed for 100% He + at a lower heat loading. This matrix study using reactor-like loading supports the need for extensive studies on future W alloy PFMs before selecting one as an alternative.

Physics↗

Be Surface Structures on W(110) and W(211): A DFT Study

Beryllium and Tungsten are promising candidates for use as plasma facing materials (PFMs) in upcoming fusion reactors. Many complex and competing phenomena however complicate the understanding and development of materials that operate in the harsh conditions of the reactor. In particular, redeposition of Be onto the W divertor must be considered due to the expected erosion of the Be first wall under exposure to the plasma. It is known that a build up of Be on W allows for the formation of BeW alloys which can harm the longevity and performance of the fusion divertor. In an effort to understand the interaction of Be with W surfaces, a study of Be structures on W(110) and W(211) as a function of Be coverage has been carried out using Density Functional Theory. We have found that both W surfaces develop monolayers of Be characterized by a densely packed hexagonal structure. Below this monolayer coverage, Be structures are found to have two motifs of coverage, a densely packed hexagonal pattern punctuated by areas of low-density coverage. The structures found here produce work function values and trends in good agreement with experimental measurements.

36 MATERIALS SCIENCE↗

Modelling the limiter ramp-up of WEST for addressing the future challenges of ITER

This paper presents a joint experimental and numerical investigation into the physics of long limited plasma ramp-up in tokamaks with tungsten (W) first walls, a critical phase for ITER operation. The comparison between the average plasma quantities simulated using the SolEdge-HDG code and the measurements taken during three successive WEST discharges after boronisation shows how challenging it is to predict this phase. While simulations reproduce the general trends at the midplane, with reasonable match in density profiles, they consistently underestimate core temperatures possibly due to too large perpendicular heat conductivity. On the contrary, at the high field side (HFS) limiter, simulations overestimate the measured quantities, and highlights the limitation of using Bohm boundary conditions at grazing magnetic angles. Experimental measurements reveal that the boron layer is rapidly eroded, on a timescale comparable to a single ITER discharge. The subsequent transition from a boron-coated to a tungsten wall increases recycling and significantly degrades the core plasma, reducing the electron temperature by nearly half due to W contamination, despite wall parameters remaining stable. Furthermore, comparisons with Langmuir probes, bolometry, reflectometry, and spectroscopy indicate that the experimental far scrape-Off layer (SOL) is significantly wider than simulated. This wide SOL implies that boron erosion extends along the entire HFS limiter rather than being confined to the contact point. This work highlights some characteristics of the plasma during this phase of the discharge and emphasizes the current modelling issues that need to be resolved in order to obtain reliable predictions concerning the ITER ramp-up.

ITER↗

Liquid metal walls

Here, the plasma performance of fusion devices depends strongly on the chosen wall materials. Solid plasma-facing components (PFCs) are predominantly used in present devices, and are the most investigated candidates for fusion designs. High-Z materials such as tungsten (W) are the leading solid PFC material candidates. To date, a material choice that scales to steady-state reactor conditions has not been identified. Moreover, if plasma transient events such as edge-localized modes and disruptions cannot be altogether avoided or sufficiently mitigated, the projected peak heat and particle loads far exceed the power exhaust capabilities of solids. Liquid metal (LM) PFCs represent an intrinsic self-healing boundary that are both resilient to surface damage from transients, and that could handle high steady-state heat and particle fluxes. Flowing LM PFCs can be designed to remove “slag,” the buildup of material erosion due to plasma-material interactions, including charge exchange sputtering in the main chamber. Further, LM offer the prospect to manage hydrogenic species otherwise retained in the PFCs, which is important from a safety and inventory standpoint. The two most promising LMs are lithium (Li) and tin (Sn), although Sn–Li eutectics may be considered. While Sn offers a higher temperature window with low vapor pressure and low hydrogen retention, Li offers the prospect of enhanced energy confinement and higher acceptable core contamination limits, and this section focuses on Li PFCs. An LM PFC development research program developed LM PFC concepts for a nuclear fusion device via engineering design calculations, single-effect experiments, and staged prototypical experiments. A self-consistent design window was identified with liquid Li flow speeds ~5–10 m/s; plasma contamination was negligible for predicted Li evolution rates. While these preconceptual designs hold promise, there is substantial R&D needed to advance the technical readiness levels of LM PFCs for application to fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗