Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge modes”

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 379 records · Page 21

Explainable tokamak-agnostic forecasting of fusion plasma instability via megahertz turbulent fluctuations

Scientific applications of artificial intelligence (AI) often remain limited by device-specific training and unexplained “black-box” approaches, creating fundamental barriers to cross-system generalization. This challenge is critical for nuclear fusion, where future reactors will have limited operational data for AI training. Here, we demonstrate that our neural network, trained solely on megahertz-scale turbulence measurements from one machine (DIII-D), forecasts Type-I edge localized mode (ELM) onsets in a different tokamak (KSTAR) through zero-shot weight transfer following physics-consistent preprocessing without device-specific retraining. Through an explainable AI framework combining gradient-weighted class activation mapping with physics validation, we reveal that our network can internalize physics relationships governing the ELM instabilities rather than memorizing device-specific patterns. The network perceives spatiotemporal features that correlate consistently with independently calculated instability growth rates, magnetohydrodynamic stability limits, and pedestal structure dynamics. Statistical analyses of dimensionally-reduced saliency features reveal the identical triangular features between the saliency representations, instability growth rates, and prediction probability across tokamaks, providing evidence that our forecasting system can show tokamak-agnostic generalization. This work contributes to a foundation for explainable scientific AI systems, where cross-system developments are essential for transcending traditional domain-specific constraints.

AI↗

The effect of high-power transient events on tungsten and tungsten coatings used for radio frequency launcher applications

High-temperature plasma-facing material coatings used for radio frequency (RF) launchers need to be robust enough to survive RF breakdown arcing or other transient events from the plasma (e.g., an edge localized mode) without causing a catastrophic failure of the coating. High-power transient effects are being explored by using an RF-induced vacuum arc to determine the robustness of tungsten coatings made by a variety of manufacturing methods. A 1/4-wavelength resonant section of vacuum transmission line terminated with an open circuit electrode structure with a well-defined electric field (30-60 kV/mm) produces repeatable arcing conditions. The initial focus is on tungsten as a plasma-facing material, including sintered tungsten, tungsten coatings on steel produced via physical vapor deposition (PVD), and functionally graded tungsten/steel coatings deposited by low-pressure plasma-spraying (LPPS). Thin PVD coatings (1-2 microns) fail catastrophically from an arc and result in severe delamination of the coating. The arc-induced damage of thicker coatings, such as those made via LPPS, tend to be restricted to the top few microns of the surface. Arcing often initiates on sharp surface microstructures and causes localized melting of tungsten at the surface of all the materials studied and results in resolidified melt pools with surface cracks. The resolidified surface results in a reduction in overall deuterium retention when exposed to typical RF plasma sheath conditions.

Caughman, John [ORNL] (ORCID:0000000206091164)↗

Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.

3-D field↗

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Volumetric recombination in EMC3-EIRENE: Implementation and first application to the pre-fusion power operation phase in ITER

An upgrade to the 3D plasma boundary model EMC3-EIRENE for detached divertor plasmas is introduced and evaluated. Stabilization of the iterative approximation of a self-consistent plasma–neutral gas solution is required at low divertor plasma temperatures, and this has been achieved by linearization of the energy loss term associated with the electron–neutral gas interaction. An implicit method for volumetric recombination has been implemented, which improves agreement with SOLPS-ITER results for the pre-fusion power operation phase in ITER. Resonant magnetic perturbations for control of edge localized modes are found to result in a partially detached divertor plasma with a significantly broader reduction of particle loads in the traditional strike zone, which is attributed to volumetric recombination. The additional non-axisymmetric far scrape-off layer strike points, on the other hand, remain attached at relatively high temperature so that volumetric recombination is not effective there.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Internal rotation of ELM filaments on NSTX

Edge localized modes (ELMs) are a threat to tokamaks due to their high heat and particle loads on plasma facing components. A significant portion of this energy is carried and deposited by the emerging ELM filaments, whose dynamics are directly connected to their impact. Therefore, understanding their underlying physics is important for the operation of future fusion reactors. Our paper extends our knowledge of ELM filaments by reporting on their internal rotation (spinning) around the magnetic field lines along which they are extended. Our analysis of gas-puff imaging data on National Spherical Torus Experiment shows that ELM filaments are characterized by internal rotation in the direction of the ion-gyromotion with ω = 15:2 krad/s median angular velocity, which is approximately three times faster than the blob rotation in the background turbulence. The characteristic size of the ELM filament was also assessed and found to be similar to the blobs. A nearly linear trend was found between the angular velocity and the radial velocity of the ELM filament. The angular velocity was found to be linearly dependent on the distance of the filament from the separatrix, as well. An analytical model called the shear-induced rotation model was identified as a candidate for explaining the physics of the observations. Our results show that the modeled mechanism could significantly influence the rotation of the ELM filament; however, it cannot be a sole contributor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements of multiple heat flux components at the divertor target by using surface eroding thermocouples (invited)

The Surface Eroding Thermocouple (SETC) is a robust diagnostic utilized in DIII-D to provide fast, edge-localized modes (ELMs) resolved heat flux measurements, in particular in geometric regions that are too shadowed for traditional infrared thermography. In order to further investigate the power dissipation in the divertor region, a combination of flush-mounted and recessed SETCs was developed to assess the effect on surface heating from non-charged particles at the divertor target. Utilizing the Divertor Materials Evaluation System sample exposure platform, the first demonstration of the feasibility of using this new method to distinguish between the heat flux from charged particles and that from neutrals and radiative heating was achieved. This paper details the process of using the combination of flush SETCs and recessed SETCs to measure the multiple heat flux components at the divertor target and further discusses how to determine two important ratios, α (ratio of heat flux from charged particles deposit on recessed SETC to that deposit on flush SETC) and β (ratio of heat flux from non-charged particles deposit on recessed SETC to that deposit on flush SETC), in the estimation of the heat flux from non-charged particle sources. Using a time dependent ratio α, it was found that ~50% of the total incident heat flux is attributable to the non-charged particles in the fully detached open divertor in DIII-D. Finally, the new application of similar SETC diagnostics in the Small Angle Slot divertor with a V-like configuration and partial tungsten coated surface (SAS-VW) is also introduced.

47 OTHER INSTRUMENTATION↗

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

Here we present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatial heterodyne spectroscopy for fast local magnetic field measurements of magnetized fusion plasmas

A novel spectroscopy diagnostic for measuring internal magnetic fields in high temperature magnetized plasmas has been developed. It involves spectrally resolving the Balmer- α (656 nm) neutral beam radiation split by the motional Stark effect with a spatial heterodyne spectrometer (SHS). The unique combination of high optical throughput (3.7 mm 2 sr) and spectral resolution (δλ ~ 0.1 nm) allows these measurements to be made with time resolution $\ll$1 ms. The high throughput is effectively utilized by incorporating a novel geometric Doppler broadening compensation technique in the spectrometer. The technique significantly reduces the spectral resolution penalty inherent to using large area, high-throughput optics while still collecting the large photon flux provided by such optics. In this work, fluxes of order 10 10 s –1 support the measurement of deviations of <5 mT (Δ λ Stark ~ 10 –4 nm) in the local magnetic field with 50 µs time resolution. Example high time resolution measurements of the pedestal magnetic field throughout the ELM cycle of a DIII-D tokamak plasma are presented. Local magnetic field measurements give access to the dynamics of the edge current density, which is essential to understanding stability limits, edge localized mode generation and suppression, and predicting performance of H-mode tokamaks.

47 OTHER INSTRUMENTATION↗

Novel angular velocity estimation technique for plasma filaments

Magnetic field aligned filaments such as blobs and edge localized mode filaments carry significant amounts of heat and particles to the plasma facing components and they decrease their lifetime. The dynamics of these filaments determine at least a part of the heat and particle loads. These dynamics can be characterized by their translation and rotation. In this paper, we present an analysis method novel for fusion plasmas, which can estimate the angular velocity of the filaments on frame-by-frame time resolution. After pre-processing, the frames are two-dimensional (2D) Fourier-transformed, then the resulting 2D Fourier magnitude spectra are transformed to log-polar coordinates, and finally the 2D cross-correlation coefficient function (CCCF) is calculated between the consecutive frames. The displacement of the CCCF’s peak along the angular coordinate estimates the angle of rotation of the most intense structure in the frame. Further, the proposed angular velocity estimation method is tested and validated for its accuracy and robustness by applying it to rotating Gaussian-structures. The method is also applied to gas-puff imaging measurements of filaments in National Spherical Torus Experiment plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calibration improvements expand filterscope diagnostic use

The filterscope diagnostic on DIII-D utilizes photomultiplier tubes to measure visible light emission from the plasma. The system has undergone a substantial upgrade since previous attempts to cross-calibrate the filterscope with other spectroscopic diagnostics were unsuccessful. The optics now utilize a dichroic mirror to initially split the light at nearly 99% transmission or reflectance for light below or above 550 nm. This allows the system to measure D α emission without degrading visible light emission from the plasma for wavelengths below 550 nm (to measure D β , D γ , W–I, C-III, etc.). Additional optimization of the optical components and calibration techniques reduce the error in the signal up to 10% in some channels compared to previous methods. Cross-calibration measurements with two other high resolution spectroscopic diagnostics now show excellent agreement for the first time. This expands the capabilities of the filterscope system allowing measurement of divertor detachment, emission profiles, edge-localized mode behavior, and plasma–wall interactions. Additionally, it enables direct comparisons against calculations from boundary plasma simulations. These were not possible before.

47 OTHER INSTRUMENTATION↗

Experimental evidence of enhanced radial transport in small ELM regimes at DIII-D

Small/type-II edge-localized-modes (ELMs), carrying 1% of the plasma stored energy, are found to deposit only 45 ± 5% of the ELM power near the strike point, and the remaining 55 ± 5% to the far scrape-off-layer (SOL). Small ELMs spread their power over a larger area compared to type-I ELMs, where such a ratio is about 60% and 40% for near- and far-SOL regions, respectively. The larger spread is reflected in the heat flux width (λq) in the SOL for the intra-small ELMs profile of 6.0 mm, almost a factor 2 larger than that of type-I ELMs of 3.15 mm, for similar plasma conditions and magnetic configuration. At the ELM peak, the small ELMs λq is found to be up to 4 times larger than for the type-I ELMs, going from 2 to 7.9 mm, indicating enhanced radial transport in the neon-seeded small ELM scenario. Inter-ELM λqs have been also calculated at the secondary outer divertor in quasi-double-null (QDN) discharges. It is found that, on average, λq is 2.2 times larger in the high-separatrix-density small ELM regime, compared to a reference type-I ELM one. These findings are supported by small ELMs radial velocity profiles, measured at the outer midplane with a fast reciprocating probe, showing a decay length (λvr) in the SOL of 12.8 cm, which is 3.3 times larger than that for the type-I ELMs of 3.9 cm. This analysis shows that small ELMs, although attractive for future machines due to low peak heat flux and large λq, might be of concern for the larger flux to the outer wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement of Stark-split beam and carbon charge exchange emissions for simultaneous B-field and temperature/rotation analysis at DIII-D

A set of two newly designed, single-channel Czerny–Turner spectrometers has been deployed at the DIII-D tokamak for measurements of the motional Stark effect (MSE) split beam emission and the C6+ (CVI) carbon charge exchange recombination (CER) emission at high spectral (δλ = 0.13 nm) and temporal (1–5 kHz) resolution. High throughput optics (f/# = 2.8) allow for good signal-to-noise at high time resolution using fast EMCCD detectors. The MSE emission allows for spectral fitting of the magnitude and direction of the local B-field, while the carbon emission yields local ion temperature and toroidal rotation information. To reduce so-called Doppler broadening of the MSE emission, a new channel-specific variable lens-masking approach has been developed. Experimental data collected from the 2023 DIII-D experimental campaign demonstrate the signal quality and instrument fidelity for both diagnostic measurements. Moreover, initial CER data analysis shows a clear evolution of the toroidal rotation during edge localized modes. Initial progress on the advanced MSE model, including a new validated ray-trace model of the DIII-D collection optics, is shown via sensitivity analysis.

Instruments & Instrumentation↗

First-wall fluxes from large and small ELMs

Experimental analysis at DIII-D shows that small edge localized modes (ELMs) deposit a larger fraction of their energy to the first wall, compared to type-I ELMs in similar magnetic configuration and input power. The energy (λQ) and particle (λΓ) flux decay lengths in the scrape-off layer (SOL) are up to 3 and 5 times larger, respectively, for small ELMs than for larger type-I ELMs. Transport dynamics of ELM filaments in the SOL are found to be related to divertor conditions, where high divertor collisionality, typical for partially detached plasmas, is associated with increased cross field ELM radial fluxes. Results show that a sufficiently large outer wall gap and/or limiters might be needed in future scenarios to protect the first wall, if operating with small ELMs and a cold divertor. This might also have implications for RF heating in future devices, where the coupling efficiency is dependent to some degree to the outer-wall gap.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗

Reconfiguration of an Electrothermal-Arc Plasma Source for In Situ PMI Studies

An electrothermal-arc plasma source (ET-Arc) has been developed to produce transient plasma heat and particle fluxes similar to those produced by edge localized modes onto divertor plasma-facing components in tokamaks. The ET-Arc utilizes a capacitive discharge to send current through a 4-mm-diameter, 9-cm-long capillary source liner. The liner material is ablated to form a high-velocity plasma jet that impacts the target downstream. With the current discharge circuit configuration, pulse lengths are 1 to 2 ms in duration and deliver heat fluxes of 0.25 to 2.1 GW m-2. The plasma was previously characterized with optical emission spectroscopy (OES) on helium emission lines. The He I line ratios were interpreted with collisional radiative analysis to calculate ne and Te. The electron temperature and electron density ranged from Te = 1 to 5 eV and ne = 1022 to 1028 electrons/m3, respectively.

Lindquist, Elizabeth G.↗

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗