Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tokamak plasma physics”

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

Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC

Here, the potential formation of multi-mega-ampere beams of relativistic 'runaway' electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium–tritium plasmas.

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 OSTI2022

Turbulent Transport Studies at ASDEX Upgrade Enabled by Correlation Electron Cyclotron Emission and nT-phase Diagnostic. Final report

Work performed on project titled, "Turbulent Transport Studies at ASDEX Upgrade Enabled by Correlation Electron Cyclotron Emission and nT-phase Diagnostic". We conducted an integrated project of experiment and model validation at ASDEX Upgrade (AUG) to expand knowledge of the underlying physics of high-performance, ELM-free regimes suitable for long pulse operation. Both the I-mode and QH-mode high-confinement regimes have potential as operating scenarios in a future reactor. Additional fluctuation measurements provided critical constraints on transport models, and elucidate the transport mechanisms associated with benevolent quasi-coherent modes in the edge and pedestal region of QH-mode and I-mode. This research expanded knowledge of turbulent-transport in tokamaks in support of ITER and other burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A tutorial on radio frequency sheath physics for magnetically confined fusion devices

Radio frequency (RF) sheaths occur under a wide variety of conditions when RF waves, material surfaces and plasma coexist. RF sheaths are of special importance in describing the interaction of ion cyclotron range of frequency (ICRF) waves with the boundary plasma in tokamaks, stellarators and other magnetic confinement devices. In this article the basic physics of RF sheaths is discussed in the context of magnetic fusion research. Techniques for modelling RF sheaths, their interaction with RF wave fields and the resulting consequences are highlighted. The article is intended as a guide for the early-career ICRF researcher, but it may equally well serve to provide an overview of basic RF sheath concepts and modelling directions for any interested fusion scientist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local transport dynamics of cold pulses in tokamak plasmas

Abstract For over two decades, our fundamental understanding of energy transport dynamics in the core of tokamak plasmas had been challenged by the striking observation of temperature perturbation reversals following the injection of cold pulses at the plasma edge. These phenomena were first discovered by Gentle et al. (Phys. Rev. Lett. 74(18):3620–3623, 1995) in 1995 and had long been suggested to be evidence of nonlocal transport effects. In recent years, a new explanation to these phenomena has emerged, fully consistent with the theory of turbulent transport in magnetized plasmas and in remarkable agreement with experiment. This article reviews the experimental observation of temperature reversals in tokamak plasmas and presents the explanation based on local transport physics.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗
DOE OSTI2024

Public Data Set: Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection

This public data set contains openly-documented, machine readable digital research data corresponding to figures published in G.M. Bodner et al., ‘Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection,’ Physics of Plasmas 28, 102504 (2021) and its erratum in Physics of Plasmas 31, 129904 (2024).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Real-time confinement regime detection in fusion plasmas with convolutional neural networks and high-bandwidth edge fluctuation measurements

Abstract A real-time detection of the plasma confinement regime can enable new advanced plasma control capabilities for both the access to and sustainment of enhanced confinement regimes in fusion devices. For example, a real-time indication of the confinement regime can facilitate transition to the high-performing wide-pedestal (WP) quiescent H-mode, or avoid unwanted transitions to lower confinement regimes that may induce plasma termination. To demonstrate real-time confinement regime detection, we use the 2D beam emission spectroscopy (BES) diagnostic system to capture localized density fluctuations of long wavelength turbulent modes in the edge region at a 1 MHz sampling rate. BES data from 330 discharges in either L-mode, H-mode, quiescent H (QH)-mode, or WP QH-mode were collected from the DIII-D tokamak and curated to develop a high-quality database to train a deep-learning classification model for real-time confinement detection. We utilize the 6×8 spatial configuration with a time window of 1024 µ s and recast the input to obtain spectral-like features via fast Fourier transform preprocessing. We employ a shallow 3D convolutional neural network for the multivariate time-series classification task and utilize a softmax in the final dense layer to retrieve a probability distribution over the different confinement regimes. Our model classifies the global confinement state on 44 unseen test discharges with an average F 1 score of 0.94, using only ∼1 ms snippets of BES data at a time. This activity demonstrates the feasibility for real-time data analysis of fluctuation diagnostics in future devices such as ITER, where the need for reliable and advanced plasma control is urgent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Integrated package of electron cyclotron emission imaging data processing and forward modeling in OMFIT

An Electron Cyclotron Emission Imaging (ECEI) data analysis module has been developed for the OMFIT platform to accommodate the needs of users at the DIII-D tokamak for physics applications. The user can easily access the ECEI spatial observation windows in the plasma that are calculated based on the automatically retrieved hardware setup and available DIII-D equilibria, perform spectral analysis, and make 2D electron temperature fluctuation images. The module provides a powerful data post processing package for extracting important physics parameters from the 2D measurements, including the radial structure and poloidal mode number of Alfven Eigenmodes, as well as the frequency-vs-wavenumber dispersion relationship of broadband MHD. Lastly, the module propagates characterized synthetic fluctuations for the user, so one can perform forward modeling tasks with simple analytical fluctuations.

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 OSTI2024

Risk-Aware Framework Development for Disruption Prediction: Alcator C-Mod and DIII-D Survival Analysis

Abstract Survival regression models can achieve longer warning times at similar receiver operating characteristic performance than previously investigated models. Survival regression models are also shown to predict the time until a disruption will occur with lower error than other predictors. Time-to-event predictions from time-series data can be obtained with a survival analysis statistical framework, and there have been many tools developed for this task which we aim to apply to disruption prediction. Using the open-source Auton-Survival package we have implemented disruption predictors with the survival regression models Cox Proportional Hazards, Deep Cox Proportional Hazards, and Deep Survival Machines. To compare with previous work, we also include predictors using a Random Forest binary classifier, and a conditional Kaplan-Meier formalism. We benchmarked the performance of these five predictors using experimental data from the Alcator C-Mod and DIII-D tokamaks by simulating alarms on each individual shot. We find that developing machine-relevant metrics to evaluate models is an important area for future work. While this study finds cases where disruptive conditions are not predicted, there are instances where the desired outcome is produced. Giving the plasma control system the expected time-to-disruption will allow it to determine the optimal actuator response in real time to minimize risk of damage to the device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas

This research is a concentrated and coordinated experimental and physics analysis program by a small U.S. team, focused on the topical area of “Long Pulse Control” as stated in the proposal solicitation. The overall goal of the research is to produce world-class analysis regarding forecasting and avoidance of disruptions in a long-pulse, superconducting tokamak at high stability parameters (e.g. normalized beta) with direct, quantifiable figures of merit (i.e. plasma disruption rate - “disruptivity”) and its reduction based on improved stability physics and control. The research targets the backbone of the problem – the basic, main issues – to produce this quantitative analysis by: (1) Creating and applying a new physics-based paradigm to understand how to eliminate disruptions: disruption event characterization, and event chain analysis – a direct manifestation of the approach suggested in the DOE Transients Report, (2) Ensuring that the components required to produce high-quality measurements and analysis needed for disruption event characterization for the tokamak in question (here, KSTAR) exist (e.g. pressure profile peaking, high-accuracy pitch angle data, compensated low frequency toroidal mode number n = 1 mode measurement) to create a primarily physics-based disruption forecasting capability, and (3) Addressing disruption avoidance for the first time in KSTAR by attacking primary disruption ‘event chains’, utilizing the most important associated plasma measurements (dynamic n = 1 field). The focus of this work will be avoidance of macroscopic instabilities, especially global modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physical mechanisms for the transition from type-III to large ELMs induced by impurity injection on EAST

Transition from type-III to large-amplitude ELMs induced by neon injection has been observed in the EAST tokamak at overlapping q 95 space between large and small ELMs. With neon injection, pedestal density gradient shows a remarkable increase accompanied by some decrease of pedestal electron temperature, and consequently the pressure gradient increases moderately and edge bootstrap current has minimal change. Further experiment demonstrates that the occurrence of large ELMs after neon injection is highly correlated with the change in edge density. Linear peeling-ballooning stability analysis indicates that the large ELM case is more unstable than the type-III ELM case during the ELM transition. A scan of pedestal density gradient in linear stability analysis shows that the direct destabilizing effect of steep pedestal density gradient on peeling-ballooning instabilities via two-fluid effects could also facilitate the transition to large ELMs. In conclusion, these results could provide more insight into the role of pedestal density gradient on pedestal stability and ELM behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal magnetic field ripple in the presence of misaligned toroidal field coils on the COMPASS-U tokamak

Homogeneity of the toroidal magnetic field is a critical aspect of the tokamak fusion device. It is affected by the design of toroidal magnetic field coils and their manufacturing and assembly misalignments. The magnetic field of the COMPASS-U project is analysed using toroidal field ripple quantity. Here ANSYS modelling of bitter-plate-designed toroidal field coils is performed to prepare an accurate simplistic representation of the coils using a virtual current line. Toroidal field ripple is analysed using a Gaussian distribution of possible misalignments. The worst-case toroidal field ripple corresponding to the technical drawings only slightly exceeds the chosen limit of 0.5%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the KSTAR tokamak

The EPEC code is employed to model the q95 windows for n = 2 and n = 1 resonant magnetic perturbation (RMP)-induced edge localized mode (ELM) suppression in typical KSTAR H-mode discharges. The plasma equilibria used in the study are derived by rescaling the experimental plasma equilibrium in KSTAR discharge #18594 measured at time t = 6450 ms. The n = 2 and n = 1 ELM suppression windows predicted by EPEC are comparatively narrow in q95 (i.e., Δq 95 ≃ 0.1) and are such that (assuming a relative pedestal pressure reduction of 15% is required to trigger ELM suppression), even in the middle of the windows, there is a threshold RMP coil current of about 1–3 kA/turn that must be exceeded before ELM suppression occurs. Furthermore, the n = 2 and n = 1 ELM suppression windows calculated by the EPEC code are consistent with the experimental observations and are also broadly similar to those recently calculated by the TM1 code [Hu et al., Phys. Plasmas 28, 052505 (2021)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗