Search NASA⌕ Search

SEARCH · Search NASA

Results for “nstx”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. II. Interaction via Landau resonance

We report that conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes driven by the Landau resonance with super-Alfvénic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite ω/ω ci and |k ∥ /k ⊥ | is retained and found to be responsible for the existence of the GAE instability via this resonance. Finally, encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

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

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Filament-assisted mode conversion in magnetized plasmas

At density filaments in magnetized plasmas, electromagnetic waves are guided along the filament ("wave-filament bound states"). Several wave-filament bound states exist beyond those predicted. Here, the new bound states occur under experimentally relevant conditions, and are especially common under typical NSTX conditions, where they are a plausible mechanism for the experimentally observed power losses in the edge plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design of microwave broadband CMOS transmitter and receiver circuits for MIR and ECEI plasma diagnostics

To efficiently determine the plasma electron density fluctuations using the MIR diagnostic technique, a 55–75 GHz 65 nm-CMOS transmitter has been developed where four separate intermediate frequency (IF) signals are up-converted, amplified, and then combined to generate an 8-tone RF output; a broadband 90 nm-CMOS receiver has also been constructed, which consists of an RF-low noise amplifier (LNA), mixer, and IF amplifier. The circuits and their corresponding modules will soon be deployed on the DIII-D and NSTX-U fusion devices. A 110–140 GHz 65 nm-CMOS receiver has also been designed, which is suitable for measuring the deep-core temperature fluctuations in the DIII-D tokamak using the electron cyclotron emission imaging diagnostic system. In addition to the RF-LNA/balun, mixer, and IF amplifier, an LO balun/tripler and driving amplifier are now included in this highly integrated circuit chip. By adopting the microwave and millimeterwave system-on-chip concept in the front-end system design, this paper demonstrates that compact transmitter and receiver modules can be easily built, which, in turn, facilitates array implementation and maintenance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Initial operation and data processing on a system for real-time evaluation of Thomson scattering signals on the Large Helical Device

A scalable system for real-time analysis of electron temperature and density based on signals from the Thomson scattering diagnostic, initially developed for and installed on the NSTX-U experiment, was recently adapted for the Large Helical Device and operated for the first time during plasma discharges. During its initial operation run, it routinely recorded and processed signals for four spatial points at the laser repetition rate of 30 Hz, well within the system’s rated capability for 60 Hz. We present examples of data collected from this initial run and describe subsequent adaptations to the analysis code to improve the fidelity of the temperature calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Installation of a solid state neutral particle analyzer array on mega ampere spherical tokamak upgrade

A compact solid state neutral particle analyzer (SSNPA) diagnostic, previously installed at NSTX-U, has been moved to MAST-U and successfully operated in the first physics campaign (MU01). The SSNPA operates by detecting the flux of fast neutral particles produced by charge exchange (CX) reactions to diagnose the fast ion distribution. The diagnostic consists of three 16-channel sensors, which provide a radial view of the plasma and have a sightline intersection with the South–South neutral beam line. From this radial geometry, active CX signals from mostly trapped particles are observed. Each channel has a spatial resolution of 3–4 cm, a temporal resolution of 200 kHz, and an average pitch angle resolution of a few degrees. The three-sensor configuration allows for coarse energy resolution of the CX signals; each sensor sees similar sightlines but different filter thicknesses alter the energy cutoffs by known amounts. Experimental data show that all channels are collecting data as intended. The signal to noise ratio is typically around 15. Preliminary data analysis shows a correlation between the SSNPA signal and magnetohydrodynamic activity in the plasma as expected.

Instruments & Instrumentation↗
DOE OSTI2023

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations

The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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

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

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

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Eigenmode twist by energy flux during fast-ion driven instabilities

Twisting of the spatial structure of eigenmodes by energy fluxes generated due to the spatial channeling (SC)—the fast-ion energy and momentum transfer across the magnetic field by destabilized modes—is considered. It is revealed that there is a correlation between the direction of the radial wave energy flux, orientation of the mode twist (MT), and the direction of mode rotation. The mode twist parameter is introduced and relations connecting it with the energy flux are established. It is shown the energy flux transforms zeros of the radial profile of the mode amplitude into minima (i.e., zeros disappear). It is found that the radial group velocity and phase velocity of reversed shear Alfvén eigenmodes (RSAEs) have opposite directions. These findings can be used for diagnostics of the energy fluxes during fast-ion driven instabilities and understanding whether the SC degrades or improves plasma performance. Specific calculations are carried out for fast magnetoacoustic modes, FMMs (known also as compressional Alfvén eigenmodes, CAEs), and Alfvén eigenmodes, AEs. A DIII-D experiment where RSAEs were observed is considered. It is concluded that an outward SC took place in this experiment. Peculiarities of various modes are discussed, which may explain why twisting of AEs, but not of FMMs, was observed experimentally in the DIII-D and NSTX tokamaks.

Kolesnichenko, Ya. I. (ORCID:0000000323219681)↗
DOE OSTI2025

Machine learning enhanced predictions of ICRF heating: Overcoming numerical limitations via data curation

In this work, we present the development of robust surrogate models for Ion Cyclotron Range of Frequencies (ICRF) and High-Harmonic Fast Wave (HHFW) heating predictions in fusion plasmas. Building upon our previous efforts to achieve real-time capable models, we identify the cause of the outliers found using TORIC in certain HHFW heating scenarios. The outliers are observed to be spurious ion Bernstein wave (IBW)-like modes caused by a wavelength control algorithm designed to address challenging scenarios with high perpendicular wavenumbers. The effect arises from the modulation in the perpendicular susceptibility, which can induce sign reversal and IBW-like propagation for scenarios featuring normalized ion Larmor radius λ i ≫ 1. We use TORIC with this algorithm disabled to generate a novel HHFW-NSTX database that is free of outliers. Surrogate models trained on this database, including Random Forest Regressor (RFR), Multi-Layer Perceptrons, and Gaussian Process Regressors (GPR), demonstrate the ability to accurately predict HHFW heating profiles, with regression scores of R 2 ∈[0.93−0.99]. Additionally we demonstrate that it is possible to generalize predictions beyond training data by the use of both RFR and GPR models, enabling the prediction of scenarios previously limited to the original model. GPR models also provide uncertainty quantification, offering insights into model confidence. This work introduces a comprehensive Verification, Validation, and Uncertainty Quantification methodology for surrogate modeling, applicable not only to ICRF heating but also to other RF heating challenges and fusion physics problems. Beyond accelerated inference, these models show effective extrapolation capabilities, providing an alternative for addressing numerical challenges.

Artificial neural networks↗
DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2022

$\mathrm{T}$ransient $\mathrm{CHI}$ System Design Studies for $\mathrm{P}$ EGASUS -$\mathrm{III}$

We report transient coaxial helicity injection (transient CHI), first developed on the Helicity Injected Torus-II (HIT-II) and later on the National Spherical Torus Experiment (NSTX) for implementing solenoid-free plasma current startup capability in a spherical tokamak (ST), is now planned to be tested on the Pegasus-III ST using a novel double-biased configuration. Such a configuration is likely needed for transient CHI deployment in a reactor. The transient CHI system optimization will be studied on Pegasus-III to enable startup toroidal persisting currents at the limits permitted by the external poloidal field coils. A transient CHI discharge is generated by driving injector current along magnetic field lines that connect the inner and outer divertor plates on one end of the ST. Simulations using the Tokamak Simulation Code are used to assess the transient CHI toroidal current generation potential and electrode gap location on the Pegasus-III. While past transient CHI systems have used high-voltage, oil-filled capacitors for driving the injector current, for improved safety, Pegasus-III will use a high-current capacitor bank based on low-voltage electrolytic capacitors. The designed and fabricated system is capable of over 32 kA. The modular design features permit the system to be upgraded to higher currents, as needed, to meet the future needs of the Pegasus-III facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resonant wave–filament interactions as a loss mechanism for HHFW heating and current drive

Abstract Perkins et al (2012 Phys. Rev. Lett. 109 045001) reported unexpected power losses during high harmonic fast wave (HHFW) heating and current drive in the National Spherical Torus Experiment (NSTX). Recently, Tierens et al (2020 Phys. Plasmas 27 010702) proposed that these losses may be attributable to surface waves on field-aligned plasma filaments, which carry power along the filaments, to be lost at the endpoints where the filaments intersect the limiters. In this work, we show that there is indeed a resonant loss mechanism associated with the excitation of these surface waves, and derive an analytic expression for the power lost to surface wave modes at each filament.

Tierens, W. (ORCID:0000000269798140)↗
DOE OSTI2024

Vertical instability forecasting and controllability assessment of multi-device tokamak plasmas in DECAF with data-driven optimization

Abstract Reliable vertical position control will be an essential element of any future tokamak-based fusion power plant in order to reduce disruptions and maximize performance. We investigate methods to improve vertical controllability boundary determination in plasma operational space and demonstrate a data-driven approach based on direct pseudoinversion of operational space data that is rigorously quantitative, applicable in real-time plasma control systems, and physically intuitive to interpret. Applied to historical shot data from entire run campaigns on the MAST-U, KSTAR, and NSTX tokamaks, this approach, implemented in DECAF, improves vertical displacement event identification accuracy to 98.9%–100%. Further, we explore the application of a physics-based vertical stability metric as an early warning forecaster for vertical displacement events. The development of a linear surrogate model for the plasma current density profile, with a coefficient of determination of 0.992 on the training dataset, enables potential employment of this forecaster in real-time. The application of this approach on historical data from the MAST-U MU02 campaign yields a forecaster with 62.6% accuracy, indicating promise for this method when further refined and potentially coupled with other stability metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗