Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma facing components”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

LLNL FESP Theory Highlights: August 2024

The 2024 ABOUND SciDAC and BOUT++ combined workshop was held August 5-9 th 2024 at the University of California Livermore Collaboration Center (UCLC) in Livermore. Bringing together leading scientists and researchers from across the globe, this pivotal event focused on advancing plasma physics and boundary plasma dynamics within the context of fusion energy research. Key discussions throughout the meeting highlighted significant advancements in the BOUT++ framework, including enhanced simulations of small Edge Localized Modes (ELMs) and the initiation of integrating the integration of the 5D GEM gyrokinetic turbulence core code with the 2D SOLPS-ITER boundary transport code. These developments are crucial for managing heat loads in fusion reactors and supporting the longevity of plasma-facing components. The event also featured a session on Inter-SciDAC Collaborations, where principal investigators from multiple U.S. FES SciDAC tokamak projects explored opportunities for cross-collaboration. Additionally, the meeting showcased cutting-edge advancements in GPU acceleration and AI/ML technologies, poised to drive the next generation of fusion research. In his closing remarks, Dr. Xueqiao Xu emphasized the importance of the collaborative efforts and discussions that took place, noting their potential to shape future breakthroughs in fusion energy. The event underscored the global nature of the BOUT++ collaboration, with contributions from over 57 institutions worldwide. The 2024 BOUT++ and ABOUND Joint Hybrid Meeting continues to drive forward the research and innovations needed to achieve fusion energy, setting the stage for future collaboration and discovery.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTX-U National Research Program: White Paper in Response to Call from FESAC Sub-Committee

Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗
DOE OSTI2023

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗
DOE OSTI2021

Modeling of ExB effects on tungsten re-deposition and transport in the DIII-D divertor

Mixed-material DIVIMP-WallDYN modelling, now incorporating ExB drifts, is presented that simultaneously reproduces tungsten (W) erosion and deposition patterns observed during the DIII-D Metal Rings Campaign, in which a toroidally symmetric set of W-coated tiles were installed in the carbon (C) DIII-D divertor. Since most reactor plasma facing component (PFC) designs call for mixed-material environments, including ITER’s W/Be enviroment, the divertor targets will quickly evolve into reconstitued surfaces of multiple elements. This work identifies controlling physics that affects material migration patterns in the divertor, which impact PFC lifetimes and impurity leakage from the divertor to the core. These simulations indicate that radial and poloidal ExB transport dominates over parallel force balance for high-Z impurities such as W in the divertor region of DIII-D. It is demonstrated that ExB drifts are required to reproduce the experimental observation of non-local W and C co-accumulation in a band ~7-9 cm outboard of the outer-strike-point (OSP) W source, for attached Lmode conditions in the unfavorable ion grad-B drift direction. In addition, W gross erosion is localized to the region outboard of the OSP, as the formation of C co-deposits suppresses W erosion at the strike point. Time-dependent simulations with scaled ExB impurity drifts (60% of the OEDGE-calculated drift velocity) and W re-erosion quantitatively reproduce these features, including depth-resolved W/C ratios, within a factor of 2 over ~115 seconds of accumulated plasma exposure. The location of co-deposition regions is shown to be well represented by an analytic leakage model, driven largely by poloidal ExB drifts. Qualitative agreement is also found between campaign-integrated W deposition measurements and simulations for the favorable ion grad-B drift direction, the standard mode of operation for most tokamaks. Furthermore, these results imply that a longterm inward radial migration of material from the outer divertor through the private flux region may occur in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The effectiveness of D 2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign

Abstract Edge localized modes (ELMs) in H-mode plasmas erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes D 2 pellet injection on the DIII-D fusion experiment used as an ELM mitigation technique applied during the 2016 tungsten Metal Rings Campaign to reduce W erosion during ELMs. The 400.9 nm photon wavelength line emission intensity of tungsten atoms (WI) filterscope channels and Langmuir probes were used to infer the gross erosion rate of tungsten-coated tiles installed in the divertor of DIII-D. D 2 mass injection rates ranging from 34 to 41 arbitrary units (A.U.) and no D 2 injection resulted in a similar total W erosion rate during ELMs (intra-ELM). On average, results show a 29% increase in the total gross W erosion rate with intermediate mass injection rates (∼13–23 A.U.) compared to the no pellets and the highest injection rate cases. On average, the fast D 2 mass injection rate cases had 15% less erosion in the inter-ELM phase than the case with no pellets. Generally, higher D 2 mass injection rates increased the ELM frequency, and the highest injection rates reduced the average erosion per ELM and fractional carbon impurities at the top of the pedestal by nearly 40% when compared to the no-pellet case. As expected, a higher D 2 pellet injection rate led to a higher plasma density and lower plasma temperature in the divertor. Additionally, an increasing divertor inter-ELM plasma electron density directly correlated to more frequent pellet injection and a decrease in both the average gross intra-ELM W erosion and the total gross intra-ELM W erosion rate. Simulations of intra-ELM erosion using the ‘free-streaming plus recycling model’ (FSRM) underestimate W erosion during pellet injection by about 30% on average. The discrepancies between the experimental measurements and the FSRM intra-ELM W erosion predictions are postulated to be due to C/W material mixing. A simple analytic mixed-material model is presented and results in better agreement with the experimental data. These results highlight the importance of incorporating the effects of a mixed-material layer in the analysis of PFC erosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Lithium–Divertor Interactions and Helium/Hydrogen Trapping in Lithiated Metals (Final Technical Report)

The goal of this project was to develop a fundamental understanding of helium and hydrogen behavior in lithium, both in bulk and at interfaces with tungsten, in order to inform the design of lithium based plasma facing components for fusion devices. Over the course of the award, the project produced the first comprehensive, peer reviewed dataset describing helium energetics, migration behaviors, and defect interactions in lithium. The research demonstrated that helium behaves in ways not previously observed in any other body centered cubic (BCC) metal: its interstitial configurations are more stable than substitutional ones, and its migration barriers are extraordinarily low, in some cases more than an order of magnitude below those in tungsten or iron. These discoveries reveal that helium in lithium diffuses so rapidly that its transport may be dominated by translational motion rather than the vibrationally activated mechanisms that underpin conventional solid state diffusion. This work lays a scientific foundation for understanding gas retention, bubble formation, and wall evolution in lithium based fusion environments and provides new computational tools, most notably a newly developed Li–He interatomic potential, for advancing future modeling efforts.

36 MATERIALS SCIENCE↗
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 OSTI2023

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Internal rotation of ELM filaments on NSTX

Edge localized modes (ELMs) are a threat to tokamaks due to their high heat and particle loads on plasma facing components. A significant portion of this energy is carried and deposited by the emerging ELM filaments, whose dynamics are directly connected to their impact. Therefore, understanding their underlying physics is important for the operation of future fusion reactors. Our paper extends our knowledge of ELM filaments by reporting on their internal rotation (spinning) around the magnetic field lines along which they are extended. Our analysis of gas-puff imaging data on National Spherical Torus Experiment shows that ELM filaments are characterized by internal rotation in the direction of the ion-gyromotion with ω = 15:2 krad/s median angular velocity, which is approximately three times faster than the blob rotation in the background turbulence. The characteristic size of the ELM filament was also assessed and found to be similar to the blobs. A nearly linear trend was found between the angular velocity and the radial velocity of the ELM filament. The angular velocity was found to be linearly dependent on the distance of the filament from the separatrix, as well. An analytical model called the shear-induced rotation model was identified as a candidate for explaining the physics of the observations. Our results show that the modeled mechanism could significantly influence the rotation of the ELM filament; however, it cannot be a sole contributor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predictive outlook for experiments resolving prompt vs local redeposition of high- Z materials in tokamaks

High-Z plasma facing components redeposit within the sheath through a combination of two distinct mechanisms: prompt (or geometric-driven) and local (or sheath-driven) redeposition. Experimental efforts are needed to determine the leading-order parameters influencing prompt-vs-local trade-off, which sets the fraction of material entering the scrape-off layer. In preparation for such experiments, leading-order parameters are isolated within the PYEAD-RustBCA-GITR coupled net erosion code using Sobol’ sensitivity analysis. Then, experiments resolving prompt-vs-local trade-off under variation of these leading-order parameters are proposed using an isotopic coupon design with multifaceted diagnostic coverage. The measurability of these experiments is evaluated using synthetic diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Initial results from boron powder injection experiments in WEST lower single null L-mode plasmas

Using a recently installed impurity powder dropper (IPD), boron powder (<150 μm) was injected into lower single null (LSN) L-mode discharges in WEST. IPDs possibly enable real-time wall conditioning of the plasma-facing components and may help to facilitate H-mode access in the full-tungsten environment of WEST. The discharges in this experiment featured I p = 0.5 MA, B T = 3.7 T, q 95 = 4.3, t pulse = 12–30 s, n e,0 ~ 4 × 10 19 m –2 , and P LHCD ~ 4.5 MW. Estimates of the deuterium and impurity particle fluxes, derived from a combination of visible spectroscopy measurements and their corresponding S/XB coefficients, showed decreases of ~50% in O + , N + , and C + populations during powder injection and a moderate reduction of these low-Z impurities (~50%) and W (~10%) in the discharges that followed powder injection. Along with the improved wall conditions, WEST discharges with B powder injection observed improved confinement, as the stored energy W MHD , neutron rate, and electron temperature T e increased significantly (10%–25% for W MHD and 60%–200% for the neutron rate) at constant input power. Notably, these increases in confinement scale up with the powder drop rate and are likely due to the suppression of ion temperature gradient (ITG) turbulence from changes in Z eff and/or modifications to the electron density profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Developing the hybrid scenario in DIII-D and KSTAR for W-compatible long-pulse operation

Long-pulse, high-performance operation with tungsten plasma-facing components is a crucial challenge for achieving steady-state (SS) goals in ITER and future fusion devices. This paper presents a joint international effort between the DIII-D and KSTAR tokamaks aimed at demonstrating the compatibility of the high q 95 ~ 6 hybrid scenario with a tungsten divertor and superconducting coils. A DIII-D hybrid scenario recipe is successfully tailored to KSTAR’s operational constraints, which facilitates cross-device studies. The results show successful 30 s long-pulse operation in KSTAR at a normalized beta (β N = 2.4) with a benign m/n = 4/3 mode and effective mitigation of tungsten accumulation. Despite this success, the KSTAR scenario exhibits a greater than 25% reduction in thermal confinement compared to the DIII-D reference. To investigate this discrepancy, we performed a predictive TGYRO validation study, which suggests that the difference is likely caused by a weaker temperature pedestal in the KSTAR discharge due to a higher and narrower density pedestal. In addition, such lower thermal confinement leads to higher impurity accumulation, possibly due to a change in neoclassical transport with increased inward pinch, which further degrades global confinement. We also report that stationary, long-pulse operation with sustained thermal confinement requires consistent fueling and wall conditioning. These findings highlight the importance of optimized fueling and real-time wall conditioning in achieving a high temperature pedestal and strong global confinement during stable long-pulse operation. We propose potential approaches toward fully SS, long-pulse operation, including real-time wall conditioning using an impurity dropper in KSTAR and future tokamaks.

DIII-D↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

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

Abstract 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, λ iz , and the width of the magnetic sheath rather than the ratio of λ iz 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 W 5 + 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗