Search NASA⌕ Search

SEARCH · Search NASA

Results for “tokamak”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Analysis of RF Sheath-Driven Tungsten Erosion at RF Antenna in the WEST Tokamak

This study applies the newly developed STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to analyze tungsten (W) erosion at RF antenna structures in the WEST tokamak. STRIPE integrates SolEdge3x for edge plasma backgrounds, COMSOL for 3D RF sheath potentials, RustBCA for sputtering yields, and GITR for impurity transport and ion energy–angle distributions. Building on prior work by Kumar et al. (2025) Nuclear Fusion, 65, 076039, which validated STRIPE for WEST ICRH discharge #57877, the present study provides a spatially resolved assessment of gross W erosion at both Q2 antenna limiters under ohmic and ICRH conditions. Simulations using 2D SolEdge3x profiles in COMSOL capture rectified sheath potentials exceeding 300 V, leading to strong upper-limiter localization. Both poloidal and toroidal asymmetries are observed and attributed to RF sheath effects, with modeled erosion patterns deviating from experiment—highlighting sensitivity to sheath geometry and plasma resolution. Erosion is driven primarily by high-charge-state oxygen ions (O6+–O8+), while D+ plays a negligible role. Assuming a plasma composition of 1% oxygen and 98% deuterium, STRIPE predicts a 30-fold increase in gross W erosion from ohmic to ICRH phases, consistent with a >25-fold rise in W-I (400.9 nm) brightness. Quantitative agreement is within 5% in the ohmic phase and 30% under ICRH, demonstrating predictive capability. Importantly, the study shows that the magnitude of ICRH-driven W erosion depends strongly on the concentration of light impurities (O, B, N, C), which drive sputtering through high charge states. Cleaner plasma conditions with reduced impurity content are therefore expected to substantially mitigate antenna W sources in WEST and other toroidal fusion devices. These findings establish STRIPE as a predictive framework for RF-induced plasma–material interactions and support its application to reactor-scale antenna design.

Kumar, Atul [ORNL] (ORCID:0000000277475210)↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Initial electron cyclotron heating/current drive scoping study for the SMall Aspect Ratio Tokamak (SMART)

In this work an initial electron cyclotron heating/current drive (EC H&CD) scoping study is presented for the SMall Aspect Ratio Tokamak (SMART), constructed and operated by the Plasma Science and Fusion Technology (PSFT) Laboratory of the University of Seville. We consider two planned phases of operation with different magnitude of magnetic field (B 0.4, and 1 T). This EC H&CD scoping study is carried out by using the ray tracing code TRAVIS, which allows to quickly perform a series of scans on the main SMART plasma scenarios. A scan in the poloidal and toroidal launching angles is performed for the B 0.4 T plasma scenario assuming three launching locations (mid-plane, off-mid-plane, and top-launcher) and the extraordinary (X-) mode polarization. The absorption occurs at the second harmonic and, overall, the highest power absorption is obtained launching the EC beam from the top of the machine. A numerical analysis of the fundamental X-mode EC H&CD start-up regime for the B 1.0 T plasma scenario with two different wave frequencies (28 and 35 GHz) is presented showing a very high current drive efficiency.

Bertelli, Nicola [Princeton Plasma Physics Laborat↗
DOE OSTI2020

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gyrokinetic simulation of turbulence and transport in the SPARC tokamak

The turbulence and transport expected in the SPARC tokamak Primary Reference Discharge (PRD) [P. Rodriguez-Fernandez et al., J. Plasma Phys. 86, 865860503 (2020)] have been investigated with the gyrokinetic code CGYRO [J. Candy et al., J. Comput. Phys. 324, 73–93 (2016)]. Linear and nonlinear simulations that focus on ion (kθρs<1.0) and electron-scale (kθρs>1.0) turbulence were used to probe the nature of the turbulence and the resulting transport in the fusion core. It is found that in the SPARC PRD, ion temperature gradient (ITG) turbulence is expected to dominate transport over most of the profile with some potential trapped electron mode impact in the near edge. Stiff turbulence is observed over a part of the plasma core such that SPARC's ion temperature profile will likely be pinned to just above the critical gradient for ITG. The role of electromagnetic turbulence, rotation, and electron-scale turbulence was investigated to provide some insight into the physics required to accurately predict SPARC performance via gyrokinetics. Additionally, predictions of impurity peaking for potential low- and high-Z SPARC first-wall materials are probed using ion-scale simulation. The dominance of low-k turbulence in SPARC provides a potential opportunity for more tractable prediction of plasma profiles using nonlinear gyrokinetics. This work is the first step toward full gyrokinetic profile prediction of SPARC kinetic profiles and the resulting fusion power and plasma gain.

Howard, N. T. (ORCID:0000000287876309)↗
DOE OSTI2022

Design elements and first data from a new Doppler backscattering system on the MAST-U spherical tokamak

A new Doppler backscattering (DBS) system has been installed and tested on the MAST-U spherical tokamak. It utilizes eight simultaneous fixed frequency probe beams (32.5, 35, 37.5, 40, 42.5, 45, 47.5, and 50 GHz). These frequencies provide a range of radial positions from the edge plasma to the core depending on plasma conditions. The system utilizes a combination of novel features to provide remote control of the probed density wavenumber, the launched polarization (X vs O-mode), and the angle of the launched DBS to match the magnetic field pitch angle. The range of accessible density turbulence wavenumbers (kθ) is reasonably large with normalized wavenumbers kθρs ranging from ≤0.5 to 9 (ion sound gyroradius ρs = 1 cm). This wavenumber range is relevant to a variety of instabilities believed to be important in establishing plasma transport (e.g., ion temperature gradient, trapped electron, electron temperature gradient, micro-tearing, kinetic ballooning modes). The system is specifically designed to address the requirement of density fluctuation wavevector alignment which can significantly reduce the SNR if not accounted for.

Instruments & Instrumentation↗
DOE OSTI2022

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of aspect ratio on plasma response to resonant magnetic perturbations in tokamak devices

A systematic numerical study is carried out, computing and comparing the plasma response to the resonant magnetic perturbation (RMP) field, applied for controlling edge localized modes (ELMs), in a series of tokamak plasmas with varying aspect ratio and utilizing the MARS-F code. The aspect ratio is scanned either by varying the plasma major radius at a fixed minor radius or by varying the latter while fixing the former. Both approaches yield similar results when compared in terms of quantities with proper normalizations. In general, a non-monotonic dependence of the resonant response field (normalized by the vacuum counterpart) near the plasma edge is found with varying aspect ratio, indicating that a given ELM control coil current configuration strongly favors plasmas with a certain aspect ratio. This optimal aspect ratio, on the other hand, depends on the toroidal as well as poloidal (i.e., coil phasing) spectra of the applied RMP field. The equilibrium (edge) safety factor, the plasma shape, and the plasma toroidal flow are all fixed to ensure that the effects identified here are predominantly due to the plasma aspect ratio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Overview of the Langmuir probe system on the Mega Ampere Spherical Tokamak (MAST) Upgrade

A detailed description of the Langmuir probe system on Mega Ampere Spherical Tokamak Upgrade is presented. The system features 850 tile-embedded probes and 40 bespoke electronic modules that each have the capability to drive and acquire data from up to 16 probes in a time-multiplexed manner. The system provides spatiotemporal-resolved measurements (1 cm and ~1 ms, respectively) in the divertor region of ion saturation current, electron temperature, and floating potential. The standard interpretation of current-voltage (IV) characteristics is to apply a four-parameter fit, based on unmagnetized probe theory, which includes a linear model for the ion saturation region. In conclusion, to mitigate the effect of the magnetic field, analysis is restricted to the region of the IV characteristic, which is sensitive to only the tail of the electron energy distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Divertor turbulent transport in the single null and snowflake in the TCV tokamak

The relative importance of divertor radial turbulent particle fluxes is considered by comparing it against the fluxes in the main-chamber outer midplane (OMP) in a variety of conditions and divertor geometries in the tokamak à configuration variable. Within the first power falloff length, the radial turbulent fluxes in the leg of the outer divertor are consistently found to be small, and about 1/5th the magnitude measured at the OMP. In a low-density single null divertor, magnetic shear is found to play a strong role in isolating the main-chamber turbulence from the divertor. The snowflake divertor is purported to have turbulence-enhancing properties in the volume between the two X-points but was instead found to further reduce the divertor turbulent fluxes compared to the single null. Depending on the collisionality, the electric field fluctuations and radial turbulent fluxes were higher near the X-point than at the outer midplane, which is likely due to the binormal compression of the flux bundles consistent with analytical models of the resistive X-point mode. Density and potential fluctuation amplitudes decrease monotonically with distance from the OMP with a slope that depends on collisionality.

Physics↗
DOE OSTI2024

Probabilistic locked mode predictor in the presence of a resistive wall and finite island saturation in tokamaks

We present a framework for estimating the probability of locking to an error field in a rotating tokamak plasma. This leverages machine learning methods trained on data from a mode-locking model, including an error field, resistive magnetohydrodynamics modeling of the plasma, a resistive wall, and an external vacuum region, leading to a fifth-order ordinary differential equation (ODE) system. It is an extension of the model without a resistive wall introduced by Akçay et al. [Phys. Plasmas 28, 082106 (2021)]. Tearing mode saturation by a finite island width is also modeled. We vary three pairs of control parameters in our studies: the momentum source plus either the error field, the tearing stability index, or the island saturation term. The order parameters are the time-asymptotic values of the five ODE variables. Normalization of them reduces the system to 2D and facilitates the classification into locked (L) or unlocked (U) states, as illustrated by Akçay et al., [Phys. Plasmas 28, 082106 (2021)]. This classification splits the control space into three regions: L̂, with only L states; Û, with only U states; and a hysteresis (hysteretic) region Ĥ, with both L and U states. In regions L̂ and Û, the cubic equation of torque balance yields one real root. Region Ĥ has three roots, allowing bifurcations between the L and U states. The classification of the ODE solutions into L/U is used to estimate the locking probability, conditional on the pair of the control parameters, using a neural network. We also explore estimating the locking probability for a sparse dataset, using a transfer learning method based on a dense model dataset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Resistive hose modes in tokamak runaway electron beams

Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear dynamics of multiple Alfvén modes driven by trapped energetic ions in tokamaks: A triplet paradigm

Non-linear dynamics of multiple infernal Alfvén eigenmodes—a subset of global Alfvén eigenmodes in tokamak plasmas with extended low-shear central core [Marchenko et al., Phys. Plasmas 16, 092502 (2009)]—is studied. The analysis is carried out for a mode triplet with toroidal mode-numbers n=1, 2, 3. It was assumed that the n = 1 mode was linearly unstable due to precession resonance with trapped fast ions, whereas the other modes were linearly damped. The modes were coupled due to a non-linearity in a bounce-averaged drift kinetic equation for the distribution function of fast ions. Nonlinear equations for the mode amplitudes and phases are derived and solved numerically. It is found that the temporal evolution of the amplitudes and the phase (responsible for the frequency chirping) of the modes exhibit Hopf bifurcations to stable limit cycles. This can explain a synchronous cyclic destabilization of multiple modes in Alfvén avalanches (sudden growth of amplitudes of the mode cluster with different n and approximately equal frequency spacing) in NSTX and bursting modes in MAST—events, which resulted in enhanced loss of fast ions.

Marchenko, V. S. (ORCID:0000000261480533)↗
DOE OSTI2024

Design and commissioning of resistive foil bolometer diagnostics on the ST40 tokamak

Following successful campaigns on the compact high field spherical tokamak ST40, resistive gold foil bolometers have now been installed to measure the radiated power profile. Positioned on the midplane, two bolometer cameras offer perspectives of the horizontal and vertical planes, while the third camera, situated above the midplane, provides a vertical view of the top divertor plates, including the X-point and the strike points. Further, these cameras use a commercial off-the-shelf data acquisition platform with modules specifically designed for bolometers in a Wheatstone bridge configuration, providing both phase sensitive detection and in situ calibration. This paper presents an overview of the design and functionality of all three cameras and the commissioning of the horizontal plane camera. The commissioning of the horizontal camera shows data from a range of plasma pulses, presenting a positive correlation between the relative profiles measured from the bolometer camera and an absolute extreme ultraviolet diode camera viewing along similar lines of sight.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First measurements of energetic protons in Mega Amp Spherical Tokamak Upgrade (MAST-U)

First proton production rates from the d(d,p)t reaction in the Mega Amp Spherical Tokamak Upgrade (MAST-U) are measured. The data were taken during the MAST-U experimental campaign with an upgraded version of the proton detector (PD) previously used in MAST. The new detector array consists of three collimated silicon surface barrier detectors with a depletion depth of 300 μm and a collimated 120 μm thick diamond detector, mounted on the MAST-U reciprocating probe arm. This array measures the energies of unconfined energetic 3 MeV protons and 1 MeV tritons mainly produced by beam-thermal DD reactions during neutral beam injection heating. Diamond detectors have the potential to be uniquely suited to detect charged fusion products as they promise to be much more radiation resistant and much less sensitive to temperature variations compared to silicon-based detectors. Using silicon and diamond-based detectors simultaneously allowed us to directly compare the performance of these two detector types. PD particle rates measured during different plasma scenarios are presented and compared to neutron rates measured using the neutron camera upgrade and TRANSP predictions.

Instruments & Instrumentation↗
DOE OSTI2024

Design and development of the magnetic diagnostic systems for the first operational phase of the SMART tokamak

A set of magnetic diagnostics has been designed, manufactured, and calibrated for the first operational phase of the small aspect ratio tokamak. The sensor suite comprises of Rogowski coils; 2D magnetic probes; and poloidal, saddle, and diamagnetic flux loops. Here, a set of continuous Rogowski coils has been manufactured for the measurement of plasma current and induced eddy currents in conductive elements. A set of flux loops and magnetic probes will be used as input for the reconstruction of the magnetohydrodynamic equilibrium. The quantity and position of these sensors have been verified to be sufficient with synthetic equilibrium reconstructions using the equilibrium fitting code and baseline scenarios computed with the Fiesta code. These sensors will also be used as input for the real-time control system, and magnetic probes will be used for the detection of plasma instabilities. The calibration procedure for the magnetic probes is described, and the results are shown. The signal conditioning and data acquisition systems are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗