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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Contamination of argon x-ray spectra by tungsten and other elements commonly found in tokamaks

Here, emission lines which appear in the spectral ranges of ground state transitions from n> = 2 levels in He- and H-like argon ions are discussed. X-ray transitions from elements commonly found in tokamaks (tungsten, molybdenum, iron and sulphur) which radiate in the wavelength range from 3700–4000 mÅ are identified by comparison with atomic structure calculations. Individual lines from tungsten charge states in the vicinity of Zn-like W 44+ are documented, along with B-like Mo 37+ . The behaviour of line ratios as a function of electron temperature is examined, in support of the identifications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The very high n Rydberg series of Ar 16+ in Alcator C-Mod tokamak plasmas

We report x-ray transitions of the very high-n Rydberg series in Ar 16+ have been observed from Alcator C-Mod tokamak plasmas. Individual emission lines up to 1s16p-1s 2 have been resolved and the central chord line brightnesses with principal quantum number n between 7 and 16 are generally found to decay as 1/n α , with α slightly larger than 3. In the plasma periphery, emission from 1s9p-1s 2 and 1s10p-1s 2 are found to be significantly enhanced relative to this decrease, indicative of selected population of these levels through charge exchange between background neutral deuterium in the ground state and Ar 17+ . An unresolved feature between the wavelengths of 1s27p-1s 2 and 1s30p-1s 2 is also present, which arises through charge exchange with neutral deuterium in the n* = 3 excited state. The brightnesses of transitions populated by charge exchange are spatially up/down asymmetric, with an excess on the side of the magnetic surface X-point. The relative brightness of the unresolved very high-n feature compared to 1s7p-1s 2 is found to increase with electron temperature and decrease with electron density. Simulations of line emission just on the long wavelength side of the Ar 16+ ionization limit indicate that the principal quantum number decay exponent is closer to α = 4 at very high n. The brightness dependence on n below 16 is in excellent agreement with calculations from the flexible atomic code package.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Oxidation of lithium plasma facing components and its effect on plasma performance in the lithium tokamak experiment-β

The characteristics of lithium-coated plasma-facing components (PFCs) have been correlated with plasma performance on the lithium tokamak experiment-β (LTX-β). Previous experiments on LTX showed that the application of lithium to PFCs was needed to achieve higher performance discharges with flat electron temperature profiles and high edge temperatures. Samples that match the LTX-β PFCs were exposed to plasmas before and after PFCs were coated with lithium and transferred under vacuum to a surface analysis station. Measurements using x-ray photoelectron spectroscopy (XPS) revealed that the primary surface constituent was lithium oxide. Earlier XPS analysis of lithium-coated PFCs on LTX was only able to show the presence of surface oxygen. The new XPS data from LTX-β have sufficient resolution to clearly identify lithium compounds for the first time, and enable them to be correlated with how lithium-coated PFCs can reduce impurities and retain hydrogen to reduce recycling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

Here, this publication presents high spatio-temporal resolution (mm µs -1 ) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the tokamak à configuration variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the D α trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ~25–35 kHz quasi-coherent density fluctuation near the separatrix ρ ψ 0.993–1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick-up probes displaying a ballooning character and n = +1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microturbulence in edge of a tokamak plasma with medium density and steep temperature gradient

Gyrokinetic simulations of electrostatic micro-turbulence have been carried out for transport barriers (TBs) in tokamak plasmas. It is found that the ion temperature gradient (ITG) mode is dominant in the pedestal with medium density gradient and steep temperature gradient. The mode width shrinks with increase of the ratio of the density and temperature gradients, which is in good agreement with the analytic theory. Unstable mode with herringbone-like structure is excited in the simulation annulus. The simulation results also indicate that multiple ITG modes are induced by a high temperature gradient. Radial electric field shear decreases the mode growth rate while increases the real frequency. Lastly, the generation of the GAM and its interaction with the turbulence are also observed in the nonlinear simulation.

tokamak↗
DOE OSTI2021

Development of a reduced model for energetic particle transport by sawteeth in tokamaks

We report that the sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Single Gaussian process method for arbitrary tokamak regimes with a statistical analysis

Abstract Gaussian process regression is a Bayesian method for inferring profiles based on input data. The technique is increasing in popularity in the fusion community due to its many advantages over traditional fitting techniques including intrinsic uncertainty quantification and robustness to over-fitting. This work investigates the use of a new method, the change-point method, for handling the varying length scales found in different tokamak regimes. The use of the Student’s t-distribution for the Bayesian likelihood probability is also investigated and shown to be advantageous in providing good fits in profiles with many outliers. To compare different methods, synthetic data generated from analytic profiles is used to create a database enabling a quantitative statistical comparison of which methods perform the best. Using a full Bayesian approach with the change-point method, Matérn kernel for the prior probability, and Student’s t-distribution for the likelihood is shown to give the best results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Spectral transition of multiscale turbulence in the tokamak pedestal

The transition in the turbulence spectrum from ion-scale dominated regimes to multiscale transport regimes that couple ion and electron scales is studied with gyrokinetic simulations of turbulent transport. The simulations are based on DIII-D high-confinement mode (H-mode) plasma parameters in the tokamak pedestal. The transition is initiated by varying the ion temperature gradient. To our knowledge, no full multiscale simulations of pedestal-like transport have been done previously. The experimental parameters lie in a bifurcation region between the two regimes. At long wavelengths, a complex, ion-direction hybrid mode is the dominant linearly unstable drift wave, while an electron temperature gradient-driven mode is unstable at short wavelengths. In the transition from the multiscale branch to the ion-scale branch, the magnitude of the ion-scale poloidal wavenumber spectrum of the nonlinear turbulent energy flux increases and the magnitude of the high-wavenumber spectrum decreases. Here, the decrease in the electron-scale transport is due to nonlinear mixing with ion-scale fluctuations and the ion-scale-driven zonal flows. A shift in the total energy associated with the fluctuating electrostatic potential intensity from dominantly drift kinetic energy in the multiscale regime to dominantly potential intensity in the ion-scale regime is well-correlated with the trend in the total energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Linear instabilities in the hot-ion regime in a high-field spherical tokamak

Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.

ubiquitous mode↗
DOE OSTI2025

Bifurcated merging operations of two spherical tokamak plasmas for reconnection heating and magnetic helicity injection

In tokamak plasma merging experiments, high-power ion heating is achieved only when the current sheet thickness is compressed to ion gyroradius ρ i . The magnetic reconnection converts 40%–50% of the reconnecting magnetic field energy (B rec 2 /2μ 0 ) into the ion heating energy (W i ) with the scaling law of W i proportional to B rec 2 . If the compressed current sheet thickness is larger than ρ i the magnetic reconnection converts only 5%–10% of B rec 2 /2μ 0 into W i . Finally, the high-power ion heating can be understood from the experimental finding that the poloidal electric field due to electrostatic quadrupolar potential in the downstream region increases linearly with both B rec and the guide field B t .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2020

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Follow the power—pathways to steady-state tokamak reactors

Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The advanced tokamak path to a compact net electric fusion pilot plant

Abstract Physics-based simulations project a compact net electric fusion pilot plant with a nuclear testing mission is possible at modest scale based on the advanced tokamak concept, and identify key parameters for its optimization. These utilize a new integrated 1.5D core-edge approach for whole device modeling to predict performance by self-consistently applying transport, pedestal and current drive models to converge fully non-inductive stationary solutions, predicting profiles and energy confinement for a given density. This physics-based approach leads to new insights and understanding of reactor optimization. In particular, the levering role of high plasma density is identified, which raises fusion performance and self-driven ‘bootstrap currents’, to reduce current drive demands and enable high pressure with net electricity at a compact scale. Solutions at 6–7 T, ∼4 m radius and 200 MW net electricity are identified with margins and trade-offs possible between parameters. Current drive comes from neutral beam and ultra-high harmonic (helicon) fast wave, though other advanced approaches are not ruled out. The resulting low recirculating power in a double null configuration leads to a divertor heat flux challenge that is comparable to ITER, though reactor solutions may require more dissipation. Strong H-mode access (x2 margin over L–H transition scalings) and ITER-like heat fluxes are maintained with ∼20%–60% core radiation, though effects on confinement need further analysis. Neutron wall loadings appear tolerable. The approach would benefit from high temperature superconductors, as higher fields would increase performance margins while potential for demountability may facilitate nuclear testing. However, solutions are possible with conventional superconductors. An advanced load sharing and reactive bucking approach in the device centerpost region provides improved mechanical stress handling. The prospect of an affordable test device which could close the loop on net-electric production and conduct essential nuclear materials and breeding research is compelling, motivating research to validate the techniques and models employed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental evidence of runaway electron tail generation via localized helical structure in pellet-triggered tokamak disruptions

A novel detector, using stacked BGO crystals, is developed for runaway electron (REs) studies in the DIII-D tokamak. It is able to resolve fast dynamics of high-energy tail formation of REs with an ultra-high time resolution of ~1 μs. As a cost, the detector estimates the `effective' energy of a given shape of γ-ray spectra and sacrifices the energy resolution. In aid of the new measurement capability, a rapid, inhomogeneous growth of RE tail is observed in detail during a major disruption triggered by an argon pellet. It is found that both the population and energy of a well-confined RE tail significantly oscillate at the early period of the growth. The oscillation phase is locked to a slow rotating magnetohydrodynamic instability, which is briefly destabilized for only ~1 ms at the early period of the current quench. The oscillation ceases promptly, when the mode disappears. As a result, the data suggests that the high-energy RE tail is well-confined and accelerated via a localized helical structure in the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global ion heating/transport during merging spherical tokamak formation

Here we report global ion heating/transport characteristics of magnetic reconnection during merging spherical tokamak formation experiment on TS-6 (TS-3U). Using the 96CH/320CH ultra high resolution ion Doppler tomography diagnostics, the full- 2 D imaging measurement clearly revealed that magnetic reconnection initially forms localized hot spots in the downstream region of outflow jet with inboard/outboard asymmetry (more deposition in the high field side) but the continuous accumulation of the heating coupled with transport process expands the high temperature region globally and forms characteristic poloidally ring-like structure aligned with field lines. The dynamic ion heating/transport process is also affected by the polarity of toroidal field and poloidally tilted/rotating global structure has experimentally been found both during and after merging. The characteristic poloidal asymmetry gets flipped when toroidal field direction is reversed and it was found that higher temperature appears in the positive potential side, which is opposite to the conventional understanding/prediction of guide field reconnection. Through the parallel acceleration process coupled with global heat transport, poloidally asymmetric non-classical feature has experimentally been found for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Operating a full tungsten actively cooled tokamak: overview of WEST first phase of operation

WEST is an MA class superconducting, actively cooled, full tungsten (W) tokamak, designed to operate in long pulses up to 1000 s. In support of ITER operation and DEMO conceptual activities, key missions of WEST are: (i) qualification of high heat flux plasma-facing components in integrating both technological and physics aspects in relevant heat and particle exhaust conditions, particularly for the tungsten monoblocks foreseen in ITER divertor; (ii) integrated steady-state operation at high confinement, with a focus on power exhaust issues. During the phase 1 of operation (2017–2020), a set of actively cooled ITER-grade plasma facing unit prototypes was integrated into the inertially cooled W coated startup lower divertor. Up to 8.8 MW of RF power has been coupled to the plasma and divertor heat flux of up to 6 MW m –2 were reached. Long pulse operation was started, using the upper actively cooled divertor, with a discharge of about 1 min achieved. This paper gives an overview of the results achieved in phase 1. Perspectives for phase 2, operating with the full capability of the device with the complete ITER-grade actively cooled lower divertor, are also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗