Search NASA⌕ Search

SEARCH · Search NASA

Results for “Safety factor gradient”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Thermoelastic analysis of solar window of a vacuum chamber

It has been noted during space simulation test programs that the chamber penetration window becomes the limiting factor as the solar intensity increases. The temperature at the center of the window normally attains a level of 550 C, while the edge temperature remains at about 35 C. The stresses produced by this thermal gradient combined with the atmospheric pressure are studied within the framework of classical theory of elasticity. The solution for the problem is obtained by solving the static uncoupled thermoelastic equations for a circular blank, subjected to nonuniform heat on one face, by introduction of two auxiliary variables. The displacements are expressed in terms of Bessel functions and parameters of interest are computed at various mesh points of the window. The maximum principal stresses thus obtained show that even with steep temperature gradient, the safety factor is acceptable.

Amba-Rao, C. L.↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Turbulence suppression at extreme plasma densities on DIII-D and EAST

Recent high-poloidal-beta (high-βP) experiments on DIII-D and EAST have made coordinated breakthroughs for high confinement quality at high density near the Greenwald limit. Density gradient amplification of turbulence suppression at high βP can explain both of these achievements. Experiments on DIII-D have achieved Greenwald fraction (fGr = line-averaged density/Greenwald density) above 1 simultaneously with normalized energy confinement (H98y2) around 1.5, as required in fusion reactor designs but never before verified in tokamak experiments with the divertor configuration. A synergy between increased H98y2 and fGr is observed with strong gas puffing, due to the build-up of an internal transport barrier at large radius in the temperature and density channels. Transport simulations reveal that the favorable trend of reduced turbulent energy transport at higher density is only expected when increasing the density gradient at high local safety factor and high β, thus at high βP to ensure strong α-stabilization. These conditions are crucial to many conceptual designs for steady-state reactors. New experiments on EAST have nearly doubled the ion temperature at fGr ∼ 0.9, consistent with predict-first modeling results based on the same physics revealed from the DIII-D analysis. All previous EAST long-pulse H-modes have Ti ≪ Te near plasma axis. Transport modeling indicates that the profiles are limited by ion-temperature-gradient modes at mid-radius. The modeling also suggested potential solutions, including reducing magnetic shear, enhancing density gradients, and higher impurity concentration. Following this guidance, EAST experiments directly show a strong enhancement of Ti achieved with a combination of a second plasma current ramp-up, a density gradient increase, and a Zeff perturbation by a short pulse (100 ms) of impurity injection, as predicted by the earlier modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of negative triangularity on microinstabilities in a low-recycling lithium-wall spherical tokamak

In this work, we present a linear gyrokinetic study of the impact of negative triangularity (NT) on microinstabilities in the Lithium Tokamak eXperiment-β (LTX-β), a low-recycling spherical tokamak with liquid lithium plasma-facing components that produce flat electron-temperature profiles [Elliott et al., IEEE Trans. Plasma Sci. 48, 1382 (2020)]. While NT is widely recognized as a stabilizing mechanism and often associated with improved confinement in conventional tokamaks, this study reveals that its effect is not universally stabilizing in the parameter regime of LTX-β and is shown to be highly sensitive to local equilibria. Using local linear simulations with the GS2 code [Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)] at ρ=0.3, 0.5, and 0.8 for two representative discharges (#103955 and #109355), and employing the Miller equilibrium model to isolate shaping effects, we find that NT can transition from stabilizing to destabilizing depending on magnetic shear, safety factor, and electron-temperature gradient. In shot #103955, NT reduces growth rates across radii, with strongest stabilization at the edge, whereas in shot #109355, it is stabilizing in the core but destabilizing at mid-radius and edge under experimental conditions. Parametric scans show that flattening the electron-temperature profile, increasing magnetic shear, and reducing the safety factor recover NT stabilization. These results demonstrate that NT stabilization is tunable rather than intrinsic and requires coordinated control of magnetic geometry and gradient drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗
DOE OSTI2023

The onset distribution of rotating m , n = 2 , 1 tearing modes and its consequences on the stability of high-confinement-mode plasmas in DIII-D

Abstract Analysis of a multi-scenario database of over 13 000 DIII-D H-mode discharges shows that the m , n = 2 , 1 magnetic islands are dominantly pressure gradient driven, stochastically triggered non-linear instabilities at all edge safety factor ( q 95 ) values. The instability onset time closely follows the exponential distribution in intermediate and high q 95 scenarios and is characterized by near constant onset rate ( λ ), in accordance with Poisson-point processes. This implies that the plasmas are operated in marginally stable conditions, characterized by a small threshold for instability growth and variations in the trigger amplitude and/or the stabilizing mechanisms with temporally uniform random distribution in this database. While the majority of the tearing modes occur in the first current-profile relaxation time of the β N flattop, constant λ throughout the β N flattop shows that the tearing onset is insensitive to the evolution of the equilibrium current profile. In low q 95 scenarios, where a large fraction of the plasmas are operated at low torque, λ increases over the course of the β N flattop, showing that these plasmas evolve toward more unstable conditions. The onset rate rapidly increases with β N , while it does not show a clear dependence on the current gradient at the mode rational surface. Overall, these observations support that the majority of the analyzed 2,1 tearing modes are non-linear, neoclassically driven instabilities and classical stability does not play a dominant role in their onset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Sizing and Topology Design of an Aeroelastic Wingbox Under Uncertainty

The goals of this work are to use a nested optimizer to conduct simultaneous sizing (inner level) and topology (outer level) design of a wingbox, considering uncertainties in the safety factors used to define the aeroelastic constraints. These uncertainties, propagated via sampling-driven polynomial chaos, are explicitly introduced at the inner level of the method, during gradient-based sizing optimization, resulting in a stochastic optimal sizing distribution. Measures of robustness in the total structural mass are then passed to the outer level, where a global optimizer evolves the topology parameters. The results demonstrate design choices needed to improve robustness in the face of uncertain safety factors, and the various physical mechanisms driving this process.

Stanford, Bret K.↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Isotope impact on Alfvén eigenmodes and fast ion transport in DIII-D

Abstract Measurements of beam driven Alfvén Eigenmode (AE) activity in matched deuterium (D) and hydrogen (H) DIII-D plasmas show a dramatic difference in unstable mode activity and fast ion transport for a given injected beam power. The dependence of the unstable AE spectrum in reversed magnetic shear plasmas on beam and thermal species is investigated in the current ramp by varying beam power in a sequence of discharges for fixed thermal and beam species at fixed density. In general, a spectrum of Reversed Shear Alfvén Eigenmodes (RSAEs) and Toroidal Alfvén Eigenmodes (TAEs) are driven unstable with sub-Alfvénic D beam injection while primarily only RSAEs are driven unstable for the H beam cases investigated. Further, for a given beam power, the driven AE amplitude is always reduced with H beams relative to D and for H thermal plasma relative to pure D or mixed D/H plasmas. Estimates of the fast ion stored energy combined with modeling using the hybrid kinetic-MHD code MEGA indicate that the dominant mechanism contributing to the difference between H and D beam drive is the faster classical slowing down of H beam ions relative to D and the resultant lower beam ion pressure. Calculations of the AE induced stored energy deficits using the reduced critical gradient model TGLFEP show quantitative agreement with the observed dependencies on injected power, isotope and minimum safety factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma

The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. Here, the analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP $β$ (NBI injection power), EP energy (NBI voltage) and the radial location of the EP density profile gradient (NBI radial deposition). The analysis indicates a transition from the soft (local plasma relaxation) to the hard MHD (global plasma relaxation) limit if the simulation EP $β$ $\unicode{x2A7E}$ $0.02$, leading to bursting MHD activity caused by radial AEs overlapping. MHD bursts cause an enhancement of the EP transport showing ballistic-like features as avalanche-like events. Simulations in the soft MHD limit show an increment of the EP density gradient as the EP $β$ increases. On the other hand, there is a gradient upper limit in the hard MHD limit, consistent with the critical-gradient behavior. AEs induce shear flows and zonal current leading to the deformation of the flux surfaces and the safety factor profile, respectively, particularly strong for the simulation in the hard MHD limit. Simulations in the hard MHD regime show a decrease of the AE frequency in the saturation phase; this is caused by the destabilization of a transitional mode between a $9/3-10/3$ TAE and a $9/3$ RSAE that may explain the AE frequency down-sweeping observed in some DIII-D discharges. Reducing the EP energy in the nonlinear simulations leads to a weakening of the plasma perturbation. On the other hand, increasing the EP energy causes the opposite effect. Nonlinear simulations of off-axis NBI profiles indicate a lower plasma perturbation as the EP density gradient is located further away from the magnetic axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

"BAAE" instabilities observed without fast ion drive

The instability that was previously identified as a fast-ion driven beta-induced Alfv´en-acoustic eigenmode (BAAE) in DIII-D was misidentified. In a dedicated experiment, low frequency modes (LFM) with characteristic “Christmas light” patterns of brief instability linked to the safety factor evolution occur in plasmas with electron temperature T e ≳ 2.1 keV but modest beta. To isolate the importance of different driving gradients on these modes, the electron cyclotron heating power and 80 keV, sub-Alfv´enic neutral beams are altered for 50-100 ms durations in reproducible discharges. Although beta-induced Alfv´en eigenmodes and reversed-shear Alfv´en eigenmodes stabilize when beam injection ceases (as expected for a fast-ion driven instability), the low frequency modes that were called BAAEs persist. Data mining reveals that characteristic LFM instabilities can occur in discharges with no beam heating but strong electron cyclotron heating. A large database of over 1000 discharges shows that LFMs are only unstable in plasmas with hot electrons but modest overall beta. The experimental LFMs have low frequencies (comparable to diamagnetic drift frequencies) in the plasma frame, occur near the minimum of the safety factor q min , and appear when q min is close to rational values. In conclusion, theoretical analysis suggests that the LFMs are a low frequency reactive instability of predominately Alfv´enic polarization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Advanced surrogate model for electron-scale turbulence in tokamak pedestals

We derive an advanced surrogate model for predicting turbulent transport at the edge of tokamaks driven by electron temperature gradient (ETG) modes. Our derivation is based on a recently developed sensitivity-driven sparse grid interpolation approach for uncertainty quantification and sensitivity analysis at scale, which informs the set of parameters that define the surrogate model as a scaling law. Our model reveals that ETG-driven electron heat flux is influenced by the safety factor q, electron beta β e and normalized electron Debye length λ D , in addition to well-established parameters such as the electron temperature and density gradients. To assess the trustworthiness of our model's predictions beyond training, we compute prediction intervals using bootstrapping. The surrogate model's predictive power is tested across a wide range of parameter values, including within-distribution testing parameters (to verify our model) as well as out-of-bounds and out-of-distribution testing (to validate the proposed model). Overall, validation efforts show that our model competes well with, or can even outperform, existing scaling laws in predicting ETG-driven transport.

fusion plasma↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗