Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tokamak devices”

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

Design and analysis of actively-cooled, edge-transport diagnostic for long-pulsed operation in WEST

Next step fusion devices that will operate in steady-state will require complex plasma-facing components (PFCs) that can survive the harsh environment over long timescales not common in current devices. This will require robust plasma facing surfaces that are integrated with active cooling systems. In a collaboration between CEA and ORNL, a plasma-interacting diagnostic is being designed for the W Environment in Steady-state Tokamak (WEST) in Cadarache, France which requires plasma-facing protection like those needed for steady-state PFCs. This integrated diagnostic studies edge transport and impurity migration within WEST, and is planned to include imbedded temperature Langmuir probe sensors, as well as removeable sample slots for ex-situ surface analysis of plasma-material interactions. The entire assembly is expected to move into the plasma edge for periods up to 1000 s having an energy removeable capability of ∼6 kW and seeing a peak heat flux of more than 7 MW/m2. The assembly compliments WEST high fluence campaigns that plan for multiple 1000 s pulses. Within these specifications, the assembly will require a refractory metal plasma facing surface and integral cooling in order to function within the limited space allotted for such diagnostics. Because of the limited space and linear actuator needs, additive manufacturing of the high heat flux working end of the assembly is being considered which could allow for precision cooling-channels and lighter weight designs. Conceptual design along with simulation and analysis results will be presented for this complex diagnostic with novel PFCs.

Lumsdaine, Arnold↗
DOE OSTI2026

Multimachine analysis of intermittent fluctuations in the scrape-off layer of magnetically confined fusion plasmas

Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) [and subsequently titled] Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022) (Final Report)

This is the Final Technical Report for the DOE Measurement Innovation grant DE-SC0017998 titled Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) and subsequently titled Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 September 2017 – 31 May 2022. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that could improve understanding of magnetic equilibrium, edge localized modes in tokamaks, and bootstrap and Pfirsch-Schlüter currents in stellarators. We have performed simulations of a novel non-hydrogen beam-based diagnostic that would determine a local value (in the plasma) of the poloidal flux function by measuring the toroidal angle of the beam particles' velocity at a detector external to the plasma. The approach is similar to that of a Heavy Ion Beam Probe (HIBP), but would be able to operate with a smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key results include: Identification of the neutral and singly charged atoms that are most appropriate for use in the beam based diagnostic and the ionization cross-sections that are important for studying the plasma edge of a large device; Development of a promising geometry for operation in which the beam is injected horizontally near the torus midplane; Implementation of simulations using banana orbit trajectories of secondary beam ions to access plasma regions that are commonly in-accessible to an HIBP diagnostic, including sample volumes near the X-point of a diverted tokamak; Assessment of the ability of the diagnostic concept (originally developed under the assumption of, and reliant on, symmetry) to work in stellarators by performing simulations with both quasi-symmetric fields and non-quasi-symmetric fields; Measurement of noise levels with a prototype detector on a stellarator and exploration of techniques to minimize effects of noise on a diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

First Evidence of Local E × B Drift in the Divertor Influencing the Structure and Stability of Confined Plasma near the Edge of Fusion Devices

The structure of the edge plasma in a magnetic confinement system has a strong impact on the overall plasma performance. Here we uncover for the first time a magnetic-field-direction dependent density shelf, i.e., local flattening of density radial profile near the magnetic separatrix, in high confinement plasmas with low edge collisionality in the DIII-D tokamak. The density shelf is correlated with a doubly peaked density profile near the divertor target plate, which tends to occur for operation with the ion BX∇B drift direction away from the X-point, as currently employed for DIII-D advanced tokamak scenarios. This double-peaked divertor plasma profile is connected via the EXB drifts, arising from a strong radial electric field induced by the radial electron temperature gradient near the divertor target. The drifts lead to the reversal of the poloidal flow above the divertor target, resulting in the formation of the density shelf. The edge density shelf can be also enhanced at higher heating power, preventing large, periodic bursts of the plasma, i.e. Edge-localized Modes, in the edge region, consistent with ideal magnetohydrodynamics calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗
DOE OSTI2020

Anomalous losses of energetic particles in the presence of an oscillating radial electric field in fusion plasmas

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. Here, it is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measuring the electron temperature and identifying plasma detachment using machine learning and spectroscopy

A machine learning approach has been implemented to measure the electron temperature directly from the emission spectra of a tokamak plasma. This approach utilized a neural network (NN) trained on a dataset of 1865 time slices from operation of the DIII-D tokamak using extreme ultraviolet/vacuum ultraviolet emission spectroscopy matched with high-accuracy divertor Thomson scattering measurements of the electron temperature, T e . This NN is shown to be particularly good at predicting T e at low temperatures (T e < 10 eV) where the NN demonstrated a mean average error of less than 1 eV. Trained to detect plasma detachment in the tokamak divertor, a NN classifier was able to correctly identify detached states (T e < 5 eV) with a 99% accuracy (an F 1 score of 0.96) at an acquisition rate 10× faster than the Thomson scattering measurement. The performance of the model is understood by examining a set of 4800 theoretical spectra generated using collisional radiative modeling that was also used to predict the performance of a low-cost spectrometer viewing nitrogen emission in the visible wavelengths. Furthermore, these results provide a proof-of-principle that low-cost spectrometers leveraged with machine learning can be used to boost the performance of more expensive diagnostics on fusion devices and be used independently as a fast and accurate T e measurement and detachment classifier.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of shape on pedestal characteristics in the wide pedestal quiescent H-mode in the DIII-D tokamak

Experiments on the DIII-D tokamak have advanced the operational limits of wide pedestal quiescent H-mode (WPQH) plasmas towards increased ITER relevance by simultaneously demonstrating well-matched plasma shape and net zero injected torque. Wide pedestal QHmodes are a compelling candidate regime for a future power producing device because they maintain a stationary pedestal without ELMs via additional edge transport. The pedestal is wider than what would be predicted from kinetic ballooning mode physics due to enhanced edge transport generated by broadband turbulence, a limit cycle oscillation, or some combination thereof. Here, compared to the double null shape, the lower single null shape is observed to have a lower density, narrower pedestal width, larger density fluctuations over a broad range of wavenumber, and operates closer to the peeling-ballooning instability boundary calculated from the simple pedestal scaling, Δ ψΝ , which is still observed to be wider than the EPED prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

I-mode pedestal relaxation events in the Alcator C-Mod and ASDEX Upgrade tokamaks

Abstract In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light on their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks — Alcator C-Mod and ASDEX Upgrade — is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, a growing oscillating precursor before the PRE onset is observed in the region close to the separatrix in both devices, and a discussion on a possible triggering mechanism is outlined. The PRE relative energy loss from the confined region is found to increase with decreasing pedestal top collisionality ν ped * . Similarly, also the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing ν ped * . Based on these relations, the PRE relative energy loss in future devices such as DEMO and ARC is estimated. Finally, the divertor peak energy fluence due to the PRE is measured on each device. Those values are then compared to the model introduced in Eich et al (2017 Nucl. Mater. Energy 12 84–90) for type-I edge localized modes. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections of the PRE divertor energy fluence to DEMO and ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Performance assessment of long-legged tightly-baffled divertor geometries in the ARC reactor concept

Extremely intense power exhaust channels are projected for tokamak-based fusion power reactors; a means to handle them remains to be demonstrated. Advanced divertor configurations have been proposed as potential solutions. Recent modelling of tightly baffled, long-legged divertor geometries for the divertor test tokamak concept, ADX, has shown that these concepts may access passively stable, fully detached regimes over a broad range of parameters. The question remains as to how such divertors may perform in a reactor setting. To explore this, numerical simulations are performed with UEDGE for the long-legged divertor geometry proposed for the ARC pilot plant conceptual design - a device with projected heat flux power width (λq||) of 0.4 mm and power exhaust of 93 MW - first for a simplified Super-X divertor configuration (SXD) and then for the actual X-point target divertor (XPTD) being proposed. It is found that the SXD, combined with 0.5% fixed-fraction neon impurity concentration, can produce passively stable, detached divertor regimes for power exhausts in the range of 80-108 MW - fully accommodating ARC's power exhaust. The XPTD configuration is found to reduce the strike-point temperature by a factor of ~10 compared to the SXD for small separations (~1.4λq||) between main and divertor X-point magnetic flux surfaces. Even greater potential reductions are identified for reducing separations to ~1λq|| or less. The power handling response is found to be insensitive to the level of cross-field convective or diffusive transport assumed in the divertor leg. By raising the separatrix density by a factor of 1.5, stable fully detached divertor solutions are obtained that fully accommodate the ARC exhaust power without impurity seeding. To our knowledge, this is the first time an impurity-free divertor power handling scenario has been obtained in edge modelling for a tokamak fusion power reactor with λq|| of 0.4 mm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First coupled GENE–XGC microturbulence simulations

Covering the core and the edge region of a tokamak, respectively, the two gyrokinetic turbulence codes Gyrokinetic Electromagnetic Numerical Experiment (GENE) and X-point Gyrokinetic Code (XGC) have been successfully coupled by exchanging three-dimensional charge density data needed to solve the gyrokinetic Poisson equation over the entire spatial domain. Certain challenges for the coupling procedure arise from the fact that the two codes employ completely different numerical methods. This includes, in particular, the necessity to introduce mapping procedures for the transfer of data between the unstructured triangular mesh of XGC and the logically rectangular grid (in a combination of real and Fourier space) used by GENE. Constraints on the coupling scheme are also imposed by the use of different time integrators. First, coupled simulations are presented. We have considered collisionless ion temperature gradient turbulence, in both circular and fully shaped plasmas. Coupled simulations successfully reproduce both GENE and XGC reference results, confirming the validity of the code coupling approach toward a whole device model. Here, many lessons learned in the present context, in particular, the need for a coupling procedure as flexible as possible, should be valuable to our and other efforts to couple different kinds of codes in pursuit of a more comprehensive description of complex real-world systems and will drive our further developments of a whole device model for fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

X-ray sources for in situ wavelength calibration of x-ray imaging crystal spectrometers

X-ray sources for a range of wavelengths are being considered for in situ calibration of X-ray Imaging Crystal Spectrometers (XICSs) and for monitoring line shifts due to changes in the crystal temperature, which can vary during experimental operation over a day [A. Ince-Cushman et al., Rev. Sci. Instrum. 79, 10E302 (2008), L. Delgado-Aparicio et al., Plasma Phys. Control. Fusion 55, 125011 (2013)]. Such crystal temperature dependent shifts, if not accounted for, could be erroneously interpreted as Doppler shifts leading to errors in plasma flow-velocity measurements. The x-ray sources encompass characteristic x-ray lines falling within the wavelength range of 0.9–4.0 Å, relevant for the XICSs on present and future fusion devices. Several technological challenges associated with the development of x-ray sources for in situ calibration are identified and are being addressed in the design of multiple x-ray tubes, which will be installed inside the spectrometer housing of the XICS for the JT-60SA tokamak. These x-ray sources will be especially useful for in situ calibration between plasma discharges. Here, in this paper, laboratory experiments are described that were conducted with a Cu x-ray source, a heated quartz (102) crystal, and a pixelated Pilatus detector to measure the temperature dependent shifts of the Cu Kα 1 and Kα 2 lines at 1.5405 and 1.5443 Å, respectively, and to evaluate the 2d-lattice constant for the Bragg reflecting crystal planes as a function of temperature, which, in the case of in situ wavelength calibration, would have to be used for numerical analysis of the x-ray spectra from the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and preliminary results of 270 GHz microwave forward scattering diagnostic system on the experimental advanced superconducting tokamak (EAST)

To measure localized (intermediate, high) poloidal wavenumber electron density fluctuations, a tangential millimeter-wave collective scattering system has been designed and successfully developed on the experimental advanced superconducting tokamak (EAST). This innovative system employs a 270 GHz mm-wave probe beam, emitted from the mid-plane of Port K and directed towards Port P after passing through a combination of two in-vessel mirrors. Here these two ports are located on the EAST device’s midplane and are 110° apart. The downward scattered signals pass through two in-vessel mirrors, and exit through the P window. The received multi-channel separated scattering signals are arranged along the poloidal direction. This optical arrangement makes the monitored turbulent wave number is poloidal dominant and up to 40 cm −1 . The tangential microwave scattering scheme enables this diagnostic to monitor local density fluctuations, with typically a radial spatial resolution △ R = 5 cm, and a wavenumber resolution of △ k = 0.4 cm −1 . The scattering system underwent comprehensive laboratory testing in 2023, with installation finalized in 2024. The system features steerable launch and receiver optics, allowing remote control to position the scattering volume from normalized radius ρ –0 to the pedestal region on a shot-by-shot basis. Effective experimental data, characterized by distinct off-center spectral peaks, have been successfully obtained in neutral beam-heated plasmas on EAST. The beam tracing forward modeling has been employed for experimental data interpretation. The advanced millimeter-wave scattering system provides fluctuation measurement capability from ion temperature gradient to electron temperature gradient scale on low field side. Preliminary experimental density fluctuation data have been successfully obtained in neutral beam-heated plasmas on EAST. The localized measurement will be used for kinetic turbulence transport numerical simulation validation.

high-k collective scattering↗
DOE OSTI2026

Predictions of m/n = 2/1 neoclassical tearing mode stabilization via high field side lower hybrid current drive on the DIII-D tokamak

Neoclassical tearing modes (NTMs) are a class of resistive instabilities that arise in tokamaks at rational surfaces and form magnetic islands. These islands flatten the pressure gradient, reducing plasma performance and may lead to disruptions if they grow large enough. Driving current within the island can stabilize the NTM, which has been achieved with electron cyclotron current drive (ECCD) on multiple devices. An alternative to ECCD is lower hybrid current drive (LHCD), which offers the advantages of increased current drive efficiency and reduced system cost. LHCD has been viewed as poorly suited for NTM suppression due to the large spatial extent of the driven current when in the multi-pass absorption regime (as has been the case in all past LHCD experiments). However, the driven current is more localized when in the single pass absorption regime, as is predicted for the DIII-D high field side (HFS) LHCD experiment. This work evaluates the feasibility of NTM suppression with HFS LHCD on DIII-D by predicting the island growth rate for a set of representative DIII-D plasmas via the modified Rutherford equation with and without the application of LHCD. In these plasmas, NTM suppression is achieved at reasonable power levels, even with finite misalignment between LH current and the island. The effect of current condensation was included and found to be most significant at smaller island sizes, assuming an experimentally typical temperature perturbation amplitude of 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a 1 mm radius vitreous carbon pellet (~3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗