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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Models of resistive wall tearing mode disruptions

Disruptions are a serious issue in tokamaks. In a disruption, the thermal energy is lost by means of an instability which could be a resistive wall tearing mode (RWTM). During precursors to a disruption, the plasma edge region cools, causing the current to contract. Model sequences of contracted current equilibria are given, and their stability is calculated. A linear stability study shows that there is a maximum value of edge q a ≈ 3 for RWTMs to occur. Further, this also implies a minimum rational surface radius normalized to plasma radius from RWTMs to be unstable. Nonlinear simulations are performed using a similar model sequence derived from an equilibrium reconstruction. There is a striking difference in the results, depending on whether the wall is ideal or resistive. With an ideal wall, the perturbations saturate at moderate amplitude, causing a minor disruption without a thermal quench. With a resistive wall, there is a major disruption with a thermal quench, if the edge q a ≤ 3. There is a sharp transition in nonlinear behavior at q a = 3. This is consistent with the linear model and with experiments. If disruptions are caused by RWTMs, then devices with highly conducting walls, such as the International Tokamak Experimental Reactor will experience much milder, tolerable, disruptions than presently predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pedestal stability analysis on MAST in preparation for MAST-U

In preparation for the upcoming MAST-U campaign, pedestal stability of spherical tokamaks is revisited by investigating standard H-mode discharges on MAST. As a step beyond previous studies, both ion and electron profiles are used for obtaining equilibria and a diverse set of pedestals is evaluated. Stability analysis with the ELITE and CGYRO codes shows that MAST pedestals are constrained by kinetic ballooning modes and medium toroidal mode number peeling-ballooning modes, with most unstable modes ranging from n = 25 to n = 45. In discharges with a steep q profile at the edge a larger number of poloidal harmonics is excited for each toroidal mode. In conclusion, a comparison with discharges on DIII-D with matched shape and similar non-dimensional parameters indicates that the increased shear at lower aspect ratio stabilizes low n peeling modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally Launched Helicon Waves

Helicon waves (a.k.a whistler waves) satisfying the normal wave-particle cyclotron resonance are observed to limit the growth and maximum energy of relativistic electrons (REs) in low-density Ohmic DIII-D tokamak plasmas. Following the application of helicon waves, pitch-angle scattering of high-energy REs causes an increase in both synchrotron and electron-cyclotron emissions. The hard x-ray emission, a proxy for the RE population, ceases to grow. Energy-resolved hard x-ray measurements also show a striking decrease in the number of high-energy REs (above the resonance at approximately 8MeV) to below the noise floor and an increase in low-energy (∼ 4 MeV) REs. This occurs despite the toroidal electric field remaining high enough to drive exponential RE growth in the absence of helicon waves. Furthermore, these results open new directions for limiting the maximum energy of RE populations in laboratory and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Enhancement of edge turbulence concomitant with ELM suppression during boron powder injection in EAST

A reproducible, quasi-stationary edge localized mode (ELM)-suppressed scenario was obtained over a wide range of plasma parameters by continuous injection of boron (B) powder into an upper-single null discharge in the experimental advanced superconducting tokamak [Sun et al., Nucl. Fusion 61, 014002 (2021)]. This powder-induced ELM-absent regime is associated with an edge harmonic mode (EHM) that provides continuous particle exhaust to maintain constant density without confinement degradation and impurity accumulation, the latter of which is often observed in ELM-free regimes. A flow rate threshold of B powder injection, leading to a threshold intensity of the EHM, is necessary for full ELM suppression. The fundamental harmonic of the EHM exhibits a toroidal mode number n = 1. The mode is observable in the entire poloidal cross section with a peak near the upper X-point in an upper-single null configuration. Here, the EHM spans radially across the pedestal and scrape-off layer, peaking inside the separatrix. The EHM appears to be insensitive to q 95 , heating power, plasma toroidal rotation, and pedestal collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

DECAF Cross-device characterization of tokamak disruptions indicated by abnormalities in plasma vertical position and current

Abnormal (deviating from target) variations in the plasma vertical position Z and current I p (such as vertical displacements, transient I p 'spikes' and quenches) constitute common elements of a disruption, a phenomenon that is to be mitigated, or ultimately avoided in future reactor-relevant tokamaks. While those abnormalities are generally recognized cross-shot and cross-device, details in terms of appearance (or not) and order of those abnormalities in disruption event chains is bound with the plasma state at the time of the chain initiation. Detection of those abnormalities is thus indicative not only of the onset of the plasma collapse itself, but also of the disruption driving cause that is promoted at the particular plasma state. Here, occurrence of disruptions, explored via detection of a Ip quench, and analysis of disruption event chains constituted by Ip and Z abnormalities, is reported for in total 7 full years of operation of 3 devices (KSTAR, MAST-U and NSTX-U) using the DECAFTM code expanded tools and capabilities. It is shown that the disruption occurrence depends not only on details of the plasma state, but also on (device-dependent) technical elements of the shot exit scenario. A year-to-year change in main disruption causes and a reduction of the disruptivity rate, bound with device and operation upgrades, are reported. Particular trigger instances of disruption event chains (and the full chains, when applicable) are shown to occupy different parts of the operation space diagrams, in accordance with prior expectation. Plasma elongation is identified as an important factor influencing details of the chains and its role will be further explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D

Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. Cryogenically frozen deuterium pellets are injected from the low-field side of the DIII-D tokamak at frequencies lower than the natural ELM frequency, and heat flux is measured by infrared cameras. Ideal MHD pedestal stability calculations show that without pellet injection, these low collisionality pedestals were limited by their current density (peeling-limited) rather than their pressure gradient (ballooning-limited). ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. Cameras observing both heat flux and D-alpha emission often find significant toroidally asymmetric striations in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets. Toroidal asymmetries at the inner strike point are similar between natural and pellet-triggered ELMs, indicating that the reduction in peak heat flux and total fluence at that location is robust for the conditions reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Projected global stability of high beta MAST-U spherical tokamak plasmas

Assessment of the limits of stability of tokamak plasmas is essential to operation in high fusion performance ranges without disruption of the plasma current. Projected equilibria have been generated for the MAST-U spherical tokamak experiment, an upgrade of the previous MAST device, in order to prepare for operation. These equilibria are scanned in pressure and current profiles, and assessed with the DCON and MARS-F stability codes to find the so-called 'no-wall' beta limit, above which resistive wall mode instabilities can be expected in the absence of other stabilising effects. The no-wall limit was generally found to increase as plasma internal inductance increased. The equilibria are also assessed for the 'with-wall' limit, theoretically the highest achievable performance point, again with the DCON and MARS-F codes, including different approximate axisymmetric walls, and with the VALEN code which includes a 3D model of the surrounding conducting structure. Similar limits were found, despite the difference between the 2D and 3D codes in the treatment of the wall. Conducting passive stabilisation plates, which were newly installed in MAST-U, are in a region of significant mode perturbation when the plasma β N is sufficiently high and eddy currents are driven in these structures. Due to the increased stabilising effect of the wall in MAST-U vs. MAST, a significant gap exists between the approximate no-wall limits of β N / l i = 6.71 and 7.13, found from DCON and MARS-F respectively, and the with-wall limits of β N / l i = 8.23 and 8.53 for the equilibrium profiles analysed in this study. This opens a region of high beta operating space in MAST-U for potentially stable operation if non-ideal effects or active control can stabilise the resistive wall mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental inference of flux tunneling between magnetic island chains in tokamaks

We report confinement bifurcations due to flux tunneling in magnetic islands in toroidal plasmas. This occurs in 2/1 islands, triggered by decoupling of non-overlapping adjacent 5/2 islands. While subject to electron cyclotron current drive (ECCD), the temperature at the 2/1 O-point is flat before but peaks after decoupling. Here, this is shown to be due to a bifurcation from stochastic to nested magnetic structure, caused by the removal of flux tunneling through intersecting manifolds of overlapping heteroclinic tangles. 2/1 island stability correlates with coupling events in an ITER baseline scenario plasma where ECCD can not stabilize the coupled islands, showing the critical impact of flux tunneling on disruption avoidance in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental Observation of Magnetic Island Heteroclinic Bifurcation in Tokamaks

We report empirical observations of magnetic island heteroclinic bifurcation for the first time. This behavior is observed in interacting coupled 2/1 tearing modes in the core of a DIII-D tokamak plasma. Poincaré maps constrained by measured magnetic amplitudes and phasing show bifurcation from heteroclinic to homoclinic topology in the 2/1 island as the 4/2 relative amplitude (R 4/2 ) decreases. Initially, the local electron temperature peak in the 2/1 island splits, consistent with two O points. As R 4/2 decreases a single peak forms, consistent with one O point. Lastly, these call for developing tearing stability theory and control solutions for heteroclinic islands in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quantifying experimental edge plasma evolution via multidimensional adaptive Gaussian process regression

The edge density and temperature of tokamak plasmas are strongly correlated with energy and particle confinement and their quantification is fundamental to understanding edge dynamics. These quantities exhibit behaviours ranging from sharp plasma gradients and fast transient phenomena (e.g. transitions between low and high confinement regimes) to nominal stationary phases. Analysis of experimental edge measurements therefore require robust fitting techniques to capture potentially stiff spatiotemporal evolution. Additionally, fusion plasma diagnostics inevitably involve measurement errors and data analysis requires a statistical framework to accurately quantify uncertainties. This paper outlines a generalized multidimensional adaptive Gaussian process routine capable of automatically handling noisy data and spatiotemporal correlations. We focus on the edge-pedestal region in order to underline advancements in quantifying time-dependent plasma profiles including transport barrier formation on the Alcator C-Mod tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Visualization of fast ion phase-space flow in plasmas well-below, near and well-above Alfvén eigenmode stability threshold in tokamak

Abstract Transport of fast ions along certain local phase space paths, referred to as fast ion phase-space flow, has been systematically measured by an imaging neutral particle analyzer (INPA) in three plasma regimes, which are well-below, near, and well-above the Alfvén eigenmode stability threshold (Du et al 2021 Phys. Rev. Lett. 127 235002). (1) In plasmas well-below the Alfvén eignenmodes (AE) stability threshold, fast ions are well-confined on passing particle orbits without noticeable transport over the phase space. The observed INPA images agree well with the synthetic INPA images, using the fast ion distribution predicted by neoclassical theory. (2) In plasmas near the AE stability threshold, INPA images in the presence of AE activity moderately deviate from those without AE activity. The image difference can be well interpreted by AE-driven, phase-space fast ion flow. Paths of this flow over the velocity space (or streamlines) are reconstructed by the intersection lines of curved E ʹ and µ surfaces, referred to as E ʹ and µ line (where E ′ ≡ E − ( ω / n ) P ζ ; E , P ζ and µ are the energy, canonical toroidal momentum and magnetic moment of ions; ω and n are the angular frequency and toroidal mode number of AEs, respectively). Resonant fast ions move radially inward by gaining energy and move radially outward by losing energy and the trajectory well aligns with the E ʹ and µ lines that pass through the mode resonances near the injection energy of neutral beams. (3) In plasmas well-above the AE stability threshold, fast ion phase-space dynamics shows additional features. Fast ions are transported out of the birth positions so promptly along the streamlines that the slowing-down process from the injection energy is not observable, exhibiting strong critical gradient behavior at local phase space. As a result, the increase of electron temperature is very small, in spite of an increase of beam power by ∼ 45 % . It should be emphasized that the directions of phase-space transport, induced by AEs with different frequencies, structures and mode numbers, do not largely differ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Power balance analysis at the L-H transition in JET-ILW NBI-heated deuterium plasmas

The understanding of the physics underlying the L-H transition has strong implications for ITER experimental reactor and demonstration power plant (DEMO). In many tokamaks, including JET, it has been observed that, at a particular plasma density, n e,min , the power necessary to access H-mode P L-H is minimum. In the present work, L-H transitions of JET deuterium plasmas heated by neutral beam injection (NBI) are studied for the first time by means of a power balance analysis to characterize the main contributions in the transition, through integrated transport modelling. In the pulses analysed, we do observe a minimum of the L-H power threshold in density, indicating the presence of density branches and of n e,min . Electron and ion heat fluxes at the transition are estimated separately. The electron/ion equipartition power results in favour of the ions, as shown by QuaLiKiz quasilinear gyrokinetic simulations, which predict a larger ion transport that causes T e > T i . The resulting edge ion heat flux also shows a clear change of slope below n e,min , similarly to ASDEX-Upgrade (AUG) NBI pulses (Ryter et al 2014 Nucl. Fusion 54 083003). JET NBI data are compared to radio-frequency heated AUG and Alcator C-mod pulses (Schmidtmayr et al 2018 Nucl. Fusion 58 056003), showing a different trend of the power, coupled to ions at the L-H transition with respect to the linearity observed in the radio-frequency heated plasmas. Finally, the presence of n e,min and the role of the ion heat flux is discussed in the paper, although it seems it is not possible to explain the presence of a P L-H minimum in density by a critical ion heat flux and by the equipartition power for the JET NBI-heated plasmas analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scenario adaptive disruption prediction study for next generation burning-plasma tokamaks

Next generation High Performance (HP) tokamaks risk damage from unmitigated disruptions at high current and power. Achieving reliable disruption prediction for a device's HP operation based on its Low Performance (LP) data is a key to its success. In this letter, through explorative data analysis and dedicated numerical experiments on multiple existing tokamaks, we demonstrate how the operational regimes of tokamaks can affect the power of a trained disruption predictor. First, our results suggest data-driven disruption predictors trained on abundant LP discharges work poorly on the HP regime of the same tokamak, which is a consequence of the distinct distributions of the tightly correlated signals related to disruptions in these two regimes. Second, we find that matching operational parameters among tokamaks strongly improves cross-machine accuracy which implies our model learns from the underlying scalings of dimensionless physics parameters like q 95 , β p and confirms the importance of these parameters in disruption physics and cross machine domain matching from the data-driven perspective. Finally, our results show in the absence of HP data from the target devices, the best predictivity of the HP regime for the target machine can be achieved by combining LP data from the target with HP data from other machines. Furthermore, these results provide a possible disruption predictor development strategy for next generation tokamaks, such as ITER and SPARC, and highlight the importance of developing baseline scenario discharges of future tokamaks on existing machines to collect more relevant disruptive data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented in this work is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameter-file driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High precision strike point control to support experiments in the DIII-D small angle slot divertor

Improvements in strike point control enabled tests with the new Small Angle Slot (SAS) divertor installed at DIII-D with precision as tight as ±1 mm, barring brief perturbations from ELMs. Aside from being narrower than previous baffle structures in DIII-D, simulations indicate that performance of the SAS divertor should be very sensitive to strike point position, making the improved strike point control essential for experiments with the SAS divertor. Improved control entails regulating more divertor geometry target locations than before and setting higher gain on the strike point target. This is achieved by tasking more poloidal field (PF) shaping coils with divertor topology control, fine tuning of the coil configuration and algorithm gains, and adding regulation of the power supply common rail. Here, the result is that although more variables are controlled, precision is better: high frequency jitter in the strike point position is lower on average compared to related discharges with standard divertor shape control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗