Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIII-D, JET”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Trends in runaway electron plateau partial recombination by massive H 2 or D 2 injection in DIII-D and JET and first extrapolations to ITER and SPARC

Abstract Experimental trends in thermal plasma partial recombination resulting from massive D 2 injection into high-Z (Ar) containing runaway electron (RE) plateaus in DIII-D and JET are studied for the purpose of achieving sufficiently low electron density ( n e ≈ 10 18 m − 3 ) to increase RE final loss MHD levels. In both DIII-D and JET, thermal electron density n e is found to drop by ∼100 × when the thermal plasma partially recombines, with a minimum at a vacuum vessel-averaged D 2 density in the range 10 20 − 10 21 m − 3 . RE effective resistivity also drops after partial recombination, indicating expulsion of the Ar content. The n e level after partial recombination is found to increase as RE current is increased. The amount of initial Ar in the RE plateau is not observed to have a strong effect on partial recombination. Partial recombination timescales of order 5 ms in DIII-D and 15 ms in JET are observed. These basic trends and timescales are matched with a 1D diffusion model, which is then used to extrapolate to ITER and SPARC tokamaks. Within the approximations of this model, it is predicted that ITER will be able to achieve sufficiently low n e values on time scales faster than expected RE plateau vertical drift timescales (of order 100 ms), provided sufficient D 2 or H 2 is injected. In SPARC, it is predicted that achieving significant n e recombination will be challenging, due to the very high RE current density. In both ITER and SPARC, it is predicted that achieving low n e will be easier with Ar as a background impurity (rather than Ne).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
DOE OSTI2025

Runaway electron-induced plasma facing component damagein tokamaks

This Roadmap article addresses the critical and multifaceted challenge of plasma facing component (PFC) damage caused by runaway electrons (REs) in tokamaks, a phenomenon that poses a significant threat to the viability and longevity of future fusion reactors such as ITER and DEMO. The dramatically increased RE production expected in future high-current tokamaks makes it very difficult to avoid or mitigate REs in such devices when a plasma discharge terminates abnormally. Preventing damage from the intense localised heat loads they can cause requires a holistic approach that considers plasma, REs and PFC damage. Despite decades of progress in understanding the physics of REs and the thermomechanical response of PFCs separately, their complex interplay remains poorly understood. This document aims to initiate a coordinated, interdisciplinary approach to bridge this gap by reviewing experimental evidence, advancing diagnostic capabilities, and improving modelling tools across different scales, dimensionalities, and fidelities. Key topics include RE beam formation and transport, damage mechanisms in both brittle and metallic PFCs, and observed effects in major facilities such as JET, DIII-D, WEST and EAST. The Roadmap emphasises the urgency of predictive, high-fidelity modelling validated against well-diagnosed controlled experiments, particularly in the light of recent changes in ITER’s wall material strategy and the growing importance of private sector fusion initiatives. Each section of the Roadmap article is written to provide a concise overview of one area of this multidisciplinary subject, with an assessment of the status, a look at current and future challenges, and a brief summary. The ultimate goal of this initiative is to guide future mitigation strategies and design resilient components that can withstand the intense localised loads imposed by REs, thus ensuring the safe and sustainable operation of the next generation of fusion power plants.

Ratynskaia, Svetlana [KTH Royal Inst. of Technolog↗
DOE OSTI2021

Recent progress in shattered pellet injection technology in support of the ITER disruption mitigation system*

Shattered pellet injection (SPI) has been selected as the baseline technology for the disruption mitigation (DM) system for ITER. Typical SPI utilizes cryogenic cooling to desublimate low pressure (<100 mbar) gases onto a cold zone within a pipe gun barrel, forming a cylindrical pellet. Pellets are dislodged from the barrel and accelerated using either a gas driven mechanical punch or high-pressure light-gas delivered by a fast-opening valve. SPI technology developed at Oak Ridge National Laboratory is currently deployed and operational on DIII-D, JET, and KSTAR. These SPI systems are used in experiments for physics scaling to ITER thermal mitigation and runaway electron dissipation/avoidance. The pellet sizes used for these machines are in the range of 4 to 12.5 mm in diameter with length to diameter ratios (L/D) of ~1.5. The current plan for ITER SPI is to utilize pellets that are 28.5 mm in diameter with an L/D of ~2. Furthermore, the large pellet sizes, high steady-state magnetic fields, and limitations of operating in a radiation environment render much of the current technology unusable. In addition to technology improvements, a deeper understanding of pellet material properties, formation, and release is being developed for implementation in future SPI designs, specifically ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A novel path to runaway electron mitigation via deuterium injection and current-driven MHD instability

Relativistic electron (RE) beams at high current density (low safety factor, qa) yet very low free-electron density accessed with D2 secondary injection in the DIII-D and JET tokamak are found to exhibit large-scale MHD instabilities that benignly terminate the RE beam. In JET, this technique has enabled termination of MA-level RE currents without measurable first-wall heating. This scenario thus offers an unexpected alternate pathway to achieve RE mitigation without collisional dissipation. Benign termination is explained by two synergistic effects. First, during the MHD-driven RE loss events both experiment and MHD orbit-loss modeling supports a significant increase in the wetted area of the RE loss. Second, as previously identified at JET and DIII-D, the fast kink loss timescale precludes RE beam regeneration and the resulting dangerous conversion of magnetic to RE kinetic energy. During the termination, the RE kinetic energy is lost to the wall, but the current fully transfers to the cold bulk thus enabling benign Ohmic dissipation of the magnetic energy on longer timescales via a conventional current quench. Hydrogenic (D2) secondary injection is found to be the only injected species that enables access to the benign termination. D2 injection: 1) facilitates access to low qa in existing devices (via reduced collisionality & resistivity), 2) minimizes the RE avalanche by ‘purging’ the high-Z atoms from the RE beam, 3) drives recombination of the background plasma, reducing the density and Alfven time, thus accelerating the MHD growth. Furthermore, this phenomenon is found to be accessible when crossing the low qa stability boundary with rising current, falling toroidal field, or contracting minor radius - the latter being the expected scenario for vertically unstable RE beams in ITER. While unexpected, this path scales favorably to fusion-grade tokamaks and offers a novel RE mitigation scenario in principle accessible with the day-one disruption mitigation system (DMS) of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗
DOE OSTI2023

MST resistive wall tearing mode simulations

Here, the Madison Symmetric Torus (MST) is a toroidal device that, when operated as a stan dard tokamak, does not have major disruptions. Unlike most tokamaks, the MST plasma is surrounded by a close fitting highly conducting wall, with a resistive wall penetration time two orders of magnitude longer than in JET or DIII-D, and three times longer than in ITER. The MST can operate with edge q a &LessSlantEqual; 2, unlike standard tokamaks. Simulations presented here indicate that the MST is unstable to resistive wall tearing modes (RWTMs) and resistive wall modes (RWMs). They could in principle cause disruptions, but the predicted thermal quench time is much longer than the experimental pulse time. If the MST thermal quench time were comparable to measurements in JET and DIII-D, theory and simulations predict that disruptions would have been observed in MST. This is consistent with the modeling herein, predicting that disruptions are caused by RWTMs and RWMs. In the low q a ~ 2 regime of MST, the RWTM asymptotically satisfies the RWM dispersion relation. The transition from RWTM to RWM occurs smoothly at q a = m/n, where m, n are poloidal and toroidal mode numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Systematic multi-machine analysis of the exhaust time-dependent behavior in tokamaks

The understanding of the time-scales and associated transient behavior of fusion exhaust plasmas plays a crucial role in its dynamic modeling and its control. This work presents an overview of experimental investigations of the exhaust dynamics in TCV, MAST-U, ASDEX-Upgrade, WEST, DIII-D, and JET. From the presented experiments, a clear picture arises on properties of the exhaust dynamics across machines. Particularly, we observe that the scrape-off layer equilibrates on fast time-scales ($>$ 70 Hz) and that exhaust dynamics measured in response to gas valve modulations mostly behave smoothly and linearly, with similarities across devices, across scenarios (H-mode, L-mode), injected species, and injection locations. The measurements presented have formed the basis for systematic exhaust control on the considered devices. We now present this database for the essential validation of dynamic exhaust models for reactor design and control.

control↗
DOE OSTI2022

A survey of pedestal magnetic fluctuations using gyrokinetics and a global reduced model for microtearing stability

This article presents a global reduced model for slab-like microtearing modes (MTMs) in the H-mode pedestal, which reproduces distinctive features of experimentally observed magnetic fluctuations, such as chirping and discrete frequency bands at noncontiguous mode numbers. Our model, importantly, includes the global variation of the diamagnetic frequencies, which is necessary to reproduce the experimental observations. The key insight underlying this model is that MTM instability is enabled by the alignment of a rational surface with the peak in the profile of the diamagnetic frequency. Conversely, MTMs are strongly stabilized for toroidal mode numbers for which these quantities are misaligned. This property explains the discrete fluctuation bands in several DIII-D and JET discharges, which we survey using our reduced model in conjunction with global gyrokinetic simulations. A fast yet accurate reduced model for MTMs enables rapid interpretation of magnetic fluctuation data from a wide range of experimental conditions to help assess the role of MTM in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multi-device analysis of energy loss duration and pellet penetration with implications for shattered pellet injection in ITER

A robust disruption mitigation system (DMS) requires accurate characterization of key disruption timescales, one of the most notable being the thermal quench (TQ). Recent modeling of shattered pellet injection (SPI) into ITER plasmas, using JOREK and INDEX, suggests long TQ durations (6–10 ms) and slow cold front propagation due to the large plasma size. If validated, these predictions would have an impact on the desired pellet parameters and mitigation strategies for the ITER DMS. To resolve these questions, a database of SPI experiments from several small-to-large sized devices (J-TEXT, KSTAR, AUG, DIII-D, and JET) has been compiled under the auspices of the International Tokamak Physics Activity MHD, disruptions, and control topical group. Analysis of the energy loss duration (proxy for the TQ duration) with machine size is presented for both mixed neon/deuterium (Ne/D) SPI and pure deuterium (D) SPI. Several metrics for the energy loss onset (e.g. soft x-ray signal drop, I p dip, and radiation flash) were considered as the conventional metric, electron cyclotron emission, is often cut-off during SPI. Several scalings with different onset metrics showed an increase in energy loss duration with machine size. The energy loss duration was additionally shown to be a function of the ratio between the number of SPI neon atoms injected and the stored energy. Analysis of the pellet shard position relative to the cold front found that in larger devices, pellets are typically found inboard of the q = 2 surface at the energy loss onset. Lastly, the delay between the pellet shards hitting the q = 2 surface and the energy loss onset was additionally found to increase with machine size. This suggests that the pellet shards in large devices will penetrate faster and further than the cooling front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Runaway electron current reconstitution after a nonaxisymmetric magnetohydrodynamic flush

Benign termination of mega-ampere (MA) level runaway current has been convincingly demonstrated in recent JET and DIII-D experiments, establishing it as a leading candidate for runaway mitigation on ITER. This comes in the form of a runaway flush by parallel streaming loss along stochastic magnetic field lines formed by global magnetohydrodynamic instabilities, which are found to correlate with a low-Z injection that purges the high-Z impurities from a post-thermal-quench plasma. Here, in this study, we show the competing physics that govern the postflush reconstitution of the runaway current in an ITER-like reactor where significantly higher current is expected. The trapped “runaways” are found to dominate the seeding for runaway reconstitution, and the incomplete purge of high-Z impurities helps drain the seed but produces a more efficient avalanche, two of which compete to produce a 2–3 MA step in current drop before runaway reconstitution of the plasma current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗
DOE OSTI2023

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Disruption Mitigation Solutions for Long Pulse Tokamaks (Final Technical Report)

The project "Disruption Mitigation Solutions for Long Pulse Tokamaks" was a collaborative effort between ORNL, GA, and UCSD. The main goal of the project was to study the injection of shattered cryogenic pellets (SPI) as a method of reducing wall damage from global tokamak plasma instabilities called disruptions. This work was primarily carried out in the JET tokamak in the UK. The team already had extensive experience studying PSI in the medium sized DIII-D tokamak in the US, but adding data from the larger JET tokamak enabled more confident extrapolation to the even larger planned ITER tokamak in France. Each of the collaborating institutions focused on different aspects of SPI disruption mitigation. The UCSD collaboration focused on the volume recombination of post-disruption runaway electron (RE) beams by massive hydrogen injection. Disruptions can sometimes result in the formation of toroidal (ring shaped) beams of highly energetic (relativistic) electron beams, which can cause serious local wall damage, even in present-day tokamaks. A promising method to reduce wall damage from the RE beams is to inject massive amounts of H 2 or D 2 into the beams using SPI. although the precise method by which the injected H 2 reducses RE wall damage is still being investigated, it seems clear that a necessary and required step for the reduction to occur is for the H 2 injection to cause volume recombination (free electron density reduction) of the background cold plasma surrounding the fast electron beam. The UCSD collaboration focused on answering the question of if this recombination could be achievable in the future reactor-size tokamak ITER using realistic levels of H 2 injection. A model capturing the most essential physics was developed and was constrained using both DIII-D and JET data. This model was then used to extrapolate to ITER. It was predicted that RE beam recombination should be achievable in ITER, although in the worst-case scenari (a 10 MA Re beam in ITER), recombination is barely achieved within the maximum allowed levels of H 2 injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of breech geometry and propellant flow on the release of large pellets for the ITER disruption mitigation system

Studies have been performed on the release mechanism for large pellets using high pressure gas in a shattered pellet injector. Typically, pellets are dislodged from the cryogenic surface and accelerated down a barrel using high pressure gas delivered by a fast-acting propellant valve. The pellets impact an angled surface which shatters the pellet into many small fragments before entering the plasma. This technique was initially demonstrated on DIII-D (Commaux et al 2016 Nucl. Fusion 56 046007) and is now deployed on JET, KSTAR, ASDEX-Upgrade, and other tokamaks around the world in support of ITER's disruption mitigation system design and physics basis. The large hydrogen, 28.5 mm diameter, 2 length-to-diameter ratio, pellets foreseen for ITER SPI operation have low material strength and low heat of sublimation, which cause the pellets to be fragile and highly reactive to the impact of warm propellant gas. Due to the size of the pellets, significantly more propellant gas is required to dislodge and accelerate them. This creates a potentially significant propellant gas removal issue as 2–6 bar-L of gas is expected to be required for release and speed control. The research presented in this paper is an in-depth exploration of the parameters that are keys to reliable pellet release and speed control. Computational fluid dynamics (CFD) modeling of propellant flows through various breech designs was conducted to determine the force generated on the back surface of a pellet. These simulations assumed the use of the ORNL designed flyer plate valve. CFD modeling combined with experimental measurements provide adequate insight to determine a path to an optimal valve and breech design for ITER SPI pellet release and speed control while minimizing propellant gas usage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Basic Physical Processes Involving Dust in Fusion Plasmas

This report presents the main outcomes of our research focused on understanding the behavior and effects of dust particles in fusion plasmas, particularly in the edge regions of tokamaks and stellarators. We conducted computer modeling studies of burst injections of carbon and tungsten dust particles in DIII-D like divertor plasmas. The studies investigated effects of transient influx of the low- and high-Z material plasma contaminants in form of dust on the edge plasma dynamics in a modern mid-size tokamak. We also performed theoretical and computational studies of the forces acting on non-spherical dust grains in magnetized and nonmagnetized plasmas. In addition, we cooperated with General Atomics slag management group on development of DIII-D dust injection experiments to measure trajectories of dust grains in the divertor plasma and compare them with DUSTT code predictions for validation of the dust modeling capabilities. We collaborated with JET experimentalists on evaluation of radiative losses induced by injection of mixed neon-deuterium ice pellets for disruption and run-away electron generation mitigation during thermal and current quench phases. We also cooperated with experimentalists at LHD fusion device (Japan) on modeling support of experiments on injection of boron granules in fusion plasma discharges to assess their dynamics and effectiveness for in situ dynamic wall conditioning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predictive turbulence-driven flux model of scrape-off layer widths across confinement regimes in tokamaks

Reliable scrape-off layer (SOL) profile decay lengths predictions are needed to design and operate future tokamaks. The present manuscript describes a new model based on turbulent transport that is able to predict SOL widths for both L-mode and H-mode plasmas. The model is based upon the sheared-spectral filament paradigm (Peret et al (WEST Team) 2022 Phys. Plasmas 29 072306), however, incorporating the effects of thermal transport in order to calculate the parallel heat fluxes. The effects of magnetic shear and ExB shear on the cross-field transport are crucial to explain the shorter SOL decay lengths found in H-mode. The model is validated against a database of thousands of DIII-D L-mode and H-mode SOL profiles. We also calculate SOL decay length predictions in terms of plasma and engineer control parameters, which are in agreement with the multi-machine empirical H-mode scaling (Eich et al (ASDEX Upgrade Team and JET EFDA Contributors) 2013 Nucl. Fusion 53 093031), however, with an additional device geometry dependence. ITER SOL width predictions by the model are 3 times higher than the empirical scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗