Search NASA⌕ Search

SEARCH · Search NASA

Results for “spherical tokamak plasma”

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

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ideal MHD Limited Electron Temperature in Spherical Tokamaks

It is well documented that the central electron temperature in the national spherical torus experiment (NSTX) remains largely unchanged as the external heating power, and hence the normalized volume averaged plasma pressure β increases [Stutman, Phys. Rev. Lett. 102, 115002 (2009)]. Herein we present a hypothesis that low n, pressure driven ideal magnetohydrodynamic (MHD) instabilities that are nondisruptive, can break magnetic surfaces in the central region and thereby flatten the electron temperature profiles. We demonstrate this mechanism in a 3D resistive MHD simulation of a NSTX discharge. By varying the toroidal magnetic field strength, and/or the heating power, we show that there is a critical value of β, above which the central temperature profile no longer peaks on axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient ECCD non-inductive plasma current start-up, ramp-up, and sustainment for an ST fusion reactor

The elimination of the need for an Ohmic heating solenoid may be the most impactful design driver for the realization of economical compact fusion tokamak reactor systems. However, this would require fully non-inductive start-up and current ramp-up from zero plasma current and low electron temperature of sub-keV to the full plasma current of ~10–15 MA at 20–30 keV electron temperature. To address this challenge, an efficient solenoid-free start-up and ramp-up scenario utilizing a low-field-side-launched extraordinary mode at the fundamental electron cyclotron harmonic frequency (X–I) is proposed, which has more than two orders of magnitude higher electron cyclotron current drive (ECCD) efficiency than the conventional ECCD for the sub-keV start-up regime. A time dependent model was developed to simulate the start-up scenarios. For the Spherical Tokamak Advanced Reactor (STAR) (Menard et al 2023 Next-Step Low-Aspect-Ratio Tokamak Design Studies (IAEA)), it was found that to fully non-inductively ramp-up to 15 MA, it would take about 25 MW of EC power at 170 GHz. Because of the relatively large plasma volume of STAR, radiation losses must be considered. It is important to make sure that high Z impurities are kept sufficiently low during the early current start-up phase where the temperature is sub-keV range. Since the initial current ramp up takes place at a factor of ten lower density compared to the sustained regimes, it is important to transition into a higher bootstrap fraction discharge at lower density to minimize the ECCD power requirement during the densification. For the sustainment phase an array of eight gyrotron launchers with a total of about 60 MW of fundamental O-mode was found to be sufficient to provide the required axis-peaked external current drive. High efficiencies between 19–57 kA MW –1 were found with optimal aiming, and these were resilient to small changes in aiming angles and density and temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗
DOE OSTI2022

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high- β are reduced to the level of that of collisionless trapped electron modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2023

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MHD-induced SOL filaments and divertor heat flux striations in NSTX

Magnetohydrodynamic (MHD) activity in NSTX has been found to induce striated heat flux patterns at the divertor. These striated heat fluxes were previously interpreted as strike-point splitting caused by stochastic magnetic lobes near the X-point. Recent analysis indicates that MHD can directly generate filaments in the scrape-off layer (SOL), and the MHD-related filaments are consistent with the rotation and the frequency of the MHD. Magnetic field line tracing from the mid-plane of the MHD-related filaments to the divertor is consistent with the location of striated heat fluxes on the divertor. Additionally, gas puff imaging (GPI) diagnostic data shows the generation of regular blobs following each edge harmonic oscillation (EHO) toroidal mode number n = 1 perturbation, which indicates that the MHD perturbation directly drives the blobs. These blobs are responsible for the divertor heat flux striations. Here, this finding provides a new framework for understanding MHD-induced edge transport and has important implications for heat-load control in NSTX-U and future spherical tokamaks.

MHD↗
DOE OSTI2021

Scoping study of detecting high harmonic fast waves in NSTX-U hot core plasma directly using beam emission spectroscopy

In this paper, we used radio frequency (RF) full wave simulation to explore the feasibility of generating a density beat oscillation formed by high harmonic fast waves (HHFWs) in the National Spherical Tokamak Experiment-Upgrade (NSTX-U) and of measuring it numerically with 2D beam emission spectroscopy (BES). This oscillation is of interest as it can be used to determine the injected RF fields. We began by computing HHFW fields in NSTX-U under a variety of experimental conditions and antenna configurations with Petra-M, a 3D RF wavefield solver. These fields were then used to evaluate the amplitude of δne,beat/ne0 in the volume of measurement, which we predict to be between 5 × 10−10 and 10−9, although it is highly dependent on plasma parameters. Trends in magnetic field, core temperature, core density, and antenna phasing are presented, and the parameters most suitable to this measurement are identified. Finally, we developed a synthetic BES diagnostic to compare the experimental BES signal against.

Rutherford, Grant (ORCID:0000000163620037)↗
DOE OSTI2024

High Current Density Cables for Simpler HTS Magnets in Fusion Energy Systems

High current density cables are needed for the engineering design of potentially low cost, simpler geometry high temperature superconducting (HTS) magnets in the promising magnetic configurations as a fusion pilot plant (FPP) option. Significant technology maturation efforts are underway by privately funded startups with the goal to demonstrate mature HTS magnet technology. Test results, however, indicate critical engineering issues remain to be addressed to meet performance goals, and demonstrate HTS coil operation repeatability and reliability. To this end, exploring and enabling multiple viable conductor and cable options is vital. Partnering with a private fusion startup and manufacturers of superconducting strands and cables, Princeton Plasma Physics Laboratory (PPPL) is exploring and seeking to de-risk the aggressive high field approach presently targeted by others. Our main objective is to develop, test and calibrate novel high current density cables for a broad deployment of affordable and reliable coils using Bi-2212 conductors. If successful, such a project will provide technical feasibility for promising FPP configurations including spherical tokamaks (ST) and compact stellarators. Here, we aim at the state-of-the-art Bi-2212 cable technologies toward a current density of 100 A/mm 2 at 16 T and 4.2 K–10 K operation for low cost, simpler geometry toroidal field (TF) coils for compact stellarators developed by fusion startup companies on a timeline consistent with the FPP initiatives and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Concept Design and Analysis of the Magnet System of the Sustained High Power Density Tokamak Facility

A recent National Academy study recommended a next-step sustained high power density (SHPD) facility for the United States that can be a bridge to a compact fusion pilot plant. A design has been initiated to investigate a possible SHPD non-deuterium-tritium device that builds upon recent low aspect ratio tokamak studies. Here, a 1.2-m device with a 2.4 aspect ratio has been chosen to evaluate physics performance goals within a double-null plasma machine arrangement centered on three component features: high current density toroidal field (TF) and central ohmic heating/poloidal field (PF) coils, liquid metal divertor/first wall systems, and the integration of a limited set of outboard dual-coolant lead lithium test blankets that can be maintained within a vertical maintenance approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Error field correction strategies in preparation to MAST-U operation

In magnetic fusion devices, unwanted non-axisymmetric magnetic field perturbations, known as error fields (EF), can have detrimental effects on plasma stability and confinement. Such EFs may originate from several sources, such as axi-symmetric coil misalignments, coil feeds, 3D structures in the wall surrounding the plasma, presence of ferromagnetic materials near the plasma surface, blanket materials in future devices. To minimize their impact on plasma performance and on the available operational space, it is important to identify the EF sources and develop EF control strategies. MAST Upgrade (MAST-U) is a spherical tokamak which will operate in the near future after a series of enhancements from the previous MAST experiment. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal eld (P) and divertor (D) coils. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal field (P) and divertor (D) coils, whose magnetic eld measurements were available when writing this paper. Such measurements reveal that P and D coils are 3D deformed and thus are responsible for intrinsic EFs, with mainly n=1 and n=2 toroidal mode numbers. In preparation to MAST-U operation, both passive and active EF control strategies have been adopted for n=1 and n=2 EFs compensation. Passive EF control consisted of installing finely each D and P coil within MAST-U device so as to minimize the intrinsic n=1 EF amplitude. The optimal coil alignment has been determined based on magnetic eld measurements and the corresponding 3D electro-magnetic modelling, and envisaged coil shift and tilt of 3.2 mm and 0.7 mrad, respectively, in the case of P coil named P4. Conversely, active control will be use during MAST-U operation to reduce the n=2 EF which is associated mainly with P4 and P5, as well. Since these coils have been re-used from the MAST device, studies attempting n=2 EF control, based on MAST plasmas, have been modelled utilizing the MARS-F code to interpret experimental results in MAST and to give hints for future n=2 EF control studies in MAST-U. Here, a model-based control set for n=2 EF control has been identified which would allow for the minimization of rotation braking, of the resonant magnetic field at the q=2 and of the plasma displacement in MAST-U, simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Multi-harmonic electron cyclotron heating and current drive scenarios for non-inductive start-up and ramp-up in high field ST-40 spherical tokamak

We report non-inductive start-up and ramp-up is an important topic for spherical tokamak reactor design as the central solenoid implementation is highly restrictive particularly for the low-aspect-ratio tokamak configuration. In the high field spherical tokamak (ST), ST-40 with B T0 ≤ 3 T, a preparation is underway for high power ECH and ECCD current start-up/ramp-up experiments utilizing two MW-class 140/105 GHz gyrotrons. Here, we explored various ECH/ECCD scenarios for a low-field-side (LFS) launch-angle steerable waveguide launcher placed near the mid-plane region. Due to the large toroidal field variation of ST configuration, multiple cyclotron harmonic resonance layers could exist within the plasma. In this start-up and ramp-up regime, both fundamental and second harmonic ECH resonances must be considered. We find that even with the presence of X-II resonance layer in the plasma, an efficient X-I ECH and ECCD regime can be accessed for the low electron temperature T e0 as low as 200 eV which is a typical starting temperature of ECH heated plasmas in an open-field-line configuration. The presence of X-II resonance could become significant at higher T e0 as X-II absorption increases with T e0 which could reduce the current ramp-up efficiency as the power reaching X-I is reduced. Finally for the pure X-I regime where the 2Ω e resonance is moved outside the plasma with B T0 ~ 3.4 T, we find that it is possible to reach the full current of I p ~ 1 MA fully non-inductively with the ECH power of ~1 MW at n e0 ~ 1.0 × 10 19 m -3 using 105 GHz frequency gyrotron. By reducing the outer limiter position R L ~ 78 cm to 70 cm, the pure X-I regime is recovered at the rated ST-40 magnetic field of B T0 ~ 3.0 T. This X-I regime is accessible with a relatively broad range of launched n II or the launching angles. A survey of X-mode X-II ECH and ECCD at higher density regimes is also shown for completeness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plasma shape and position control development for NSTX-U using the GSEvolve plasma simulator

NSTX-U discharges have been reproduced using the GSEvolve plasma simulator in order to improve their control and reduce oscillations during the plasma current (Ip) ramp-up phase, and with the ultimate goal of enabling reliable H-mode access with a double-null shape at high elongation, low internal inductance, and high plasma performance. The control objective is to eliminate undesired oscillations in the plasma shape and vertical position that start when the plasma diverts during the Ip ramp-up phase, and that happened recurrently during the first NSTX-U experimental campaign. Model-based control solutions have been developed, including alternative H-mode access recipes and a Multi-Input Multi-Output (MIMO) decoupling controller. These solutions have demonstrated their capability to reduce or eliminate the aforementioned oscillations in GSEvolve simulations, enabling reproducible H-mode access. The results of this work aim to inform NSTX-U experimental operations on potential control developments to maximize its scientific output, and highlight the powerful capabilities of high-fidelity control-simulation frameworks like GSEvolve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The impact of kinetic and global effects on ideal ballooning 2nd stable pedestals of conventional and low aspect-ratio tokamaks

The EPED model (Snyder et al 2011 Nucl. Fusion 51 103016) has had success in describing the pedestal structure of Type-I ELM and QH-mode plasmas in tokamaks, by combining kinetic ballooning mode (KBM) and peeling–ballooning (PB) constraints. Within EPED, the KBM constraint is often approximated by a technique using calculated ideal ballooning mode (IBM) thresholds, fit to a functional form designed to capture non-local effects. It has been noted that quantitative differences between local ideal MHD and gyro-kinetic (GK) ballooning stability can be larger at low aspect ratio. KBM critical pedestals are consistent with observations in initial studies on conventional and spherical tokamaks. In this work, the application of a reduced model for the calculation of the kinetic ballooning stability boundary is presented based on a novel and newly developed gyro-fluid system (GFS) code (Staebler et al 2023 Phys. Plasmas 30 102501). GFS is observed to capture KBMs in DIII-D as well as the NSTX pedestals, opening a route to integrating this model into EPED. Finally, high but finite $n$ global ballooning modes are observed to limit the access to the local 2nd stability and thus provide a transport mechanism that constrains the pedestal evolution with $β_{p,\textrm{ped}}$. The high $n$ global ballooning stability is approximated by its ideal MHD analog using ELITE. It is shown that nearly-local high $n$ modes with $k_yρ_s$ ~ $0.25 - 0.5$ can provide a proxy for the critical $β_{p,\textrm{ped}}$ when a 2nd stable access exists on DIII-D plasmas. The use of GFS and ELITE scaling in EPED provides improved agreement to EPED1 in an initial comparison with DIII-D pedestal data.

Anastopoulos Tzanis, Michail Savvas [Oak Ridge Nat↗