Search NASA⌕ Search

SEARCH · Search NASA

Results for “neoclassical”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Early internal detection of magnetic tearing and implications for tokamak magnetohydrodynamic stability

Internal fluctuation measurements with Faraday-effect polarimetry in the DIII-D tokamak reveal the onset of a tearing mode with toroidal mode number n = 3 well before it is detected by the sensing coils external to the plasma. This mode appears before the n = 2, 1 modes and is first detected with internal measurements at a lower value of the ideal-wall kink beta limit than is indicated at the time of first detection by the sensing coils. When the mode is first detected, the linear resistive stability parameter, Δ′, indicates marginal stability and continues to do so until later when the mode amplitude begins increasing linearly with time—together suggesting a neoclassical origin for this mode.

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 OSTI2026

Enhancement of mid-/high-Z impurity transport by continuous Li-granule dropping in a stellarator plasma

An enhancement of core impurity transport is observed in high-density plasmas of the stellarator large helical device, heated by neutral beam injection, when continuous lithium (Li) granule dropping is performed. In the reported experiments, in which the tracer-encapsulated solid pellet is employed to inject trace amounts of titanium (Ti) and molybdenum (Mo) into the plasma core, confinement times for these impurities are seen to reduce significantly when Li dropping is applied, this reduction being more notable for Mo. To gain some initial insight into these observations, simulations are performed using the drift-kinetic transport code SFINCS for the Mo case. These simulations indicate that, while neoclassical transport prevails for the main plasma components (electrons, majority ions, and low-Z impurities), the classical contribution appears dominant for transporting Mo impurities. In summary, this work reports the first experimental observation of the degradation of mid-Z and high-Z impurity confinement induced by the continuous dropping of Li granules into a high-density stellarator plasma. In the case of the Mo impurity, simulations suggest that classical transport is the key mechanism underlying the enhanced impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of a synthetic phase contrast imaging diagnostic for turbulence studies at Wendelstein 7-X

Abstract We present a synthetic phase contrast imaging (PCI) diagnostic for studying turbulence at the Wendelstein 7-X (W7-X) stellarator. We first describe the implemented instrument response model, which captures diffraction effects, detector noise, and the long-wavelength cutoff due to the phase plate of the PCI system. To verify the instrument response model, we show that it is capable of reproducing the PCI signal generated by the sound wave speaker used for calibration at W7-X. Next, we discuss the calculation of synthetic PCI signals based on the global, nonlinear gyrokinetic codes GENE-3D and EUTERPE, including results from some of the first stellarator simulations of this type with kinetic electrons (KEs) in GENE-3D. While the simulations used in this work lack a neoclassical radial electric field, which is crucial for reproducing experimental PCI signals, they do indicate that the dominant rotation direction and velocities of the turbulent fluctuations can be inferred from the wave number-frequency spectra of the PCI signals, as expected. The synthetic PCI wave number spectra are further shown to be similar to those of the line-integrated fluctuating electron density, with distinct differences between adiabatic and KE simulations, explainable by previously published turbulence models. For example, the wave number spectra of all adiabatic electron simulations analyzed here follow a power law with an exponent close to −5 for sufficiently large wave numbers. This indicates that universal features of electron density turbulence at W7-X may be studied using the PCI system.

EUTERPE↗
DOE OSTI2022

Comparison of Doppler back-scattering and charge exchange measurements of E × B plasma rotation in the DIII-D tokamak under varying torque conditions

Measurements of the E × B toroidal angular velocity, $\omega_{E \times B}\,=\,E_\mathrm{r}/RB_\theta$ ($E_\mathrm{r}$ is the radial electric field, B θ is the poloidal magnetic field), are made using the Doppler back-scattering (DBS) and charge-exchange recombination (CER) spectroscopy diagnostics. DBS uses the Doppler shift of wavenumber-resolved density fluctuations while CER uses the Doppler shift of impurity emission lines to independently measure plasma parameters for calculating the local radial electric field. DBS and CER profiles of $\omega_{E \times B}$ as a function of normalized toroidal flux (ρ) are compared at various levels of neutral beam applied torque on the plasma. Under standard neoclassical theory $\omega_{E\times B}$ is a flux surface quantity, making it appropriate to compare across diagnostics. DBS and CER generally show good agreement when comparing $\omega_{E \times B}$ profiles at different levels of neutral beam injection-applied torque. Furthermore, the DBS values have close to the same precision as CER values when averaged over a similar time-scale and effects, such as prompt-torque are considered. DBS is able to observe the rapid ($\lt$10 ms) modification of the $E_\mathrm{r}$ profile by the diagnostic neutral beam 'blips'. This modification is most pronounced when the blip applies a large relative change in torque on the plasma. Overall, these results could have implications on transport analysis and suggests using DBS and CER in conjunction to constrain values of the E × B-shear (sometimes called $\gamma_{E \times B}$).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗
DOE OSTI2023

Pedestal particle balance studies in JET-ILW H-mode plasmas

JET-ILW type I ELMy H-modes at 2.5 MA/2.8 T with constant NBI heating (23 MW) and gas fuelling rate were performed, utilising edge localised mode (ELM) pacing by vertical kicks and plasma shaping (triangularity, δ) as tools to disentangle the effects of ELMs, inter-ELM transport and edge stability on the pedestal particle balance. In agreement with previous studies, the pedestal confinement improves with increasing δ, mostly due to a significant increase in pedestal density while the ELM frequency ($f_{\mathrm{ELM}}$) is decreased. Improved pedestal confinement with increasing δ was observed even when the pedestal MHD stability was degraded artificially by vertical kicks, implying that increased triangularity may favourably affect the inter-ELM pedestal recovery. The workflow developed to quantify the pedestal particle balance uses high time-resolution profile reflectometry to characterise the inter-ELM evolution of the plasma particle content ($\mathrm{d} N/\mathrm{d}t$), the NEO drift-kinetic solver to evaluate the neoclassical fluxes and interpretative EDGE2D-EIRENE simulations to estimate the edge particle source. The edge particle source is then constrained by deuterium Balmer-α line intensity measurements in the main chamber, which are, however, strongly affected by reflections from the metal walls. The reflections are accounted for by the CHERAB code taking the divertor emission (the brightest light source in the torus) distribution from imaging spectroscopy measurements as input. Our analysis shows that in the second half of the ELM cycle, the volume-integrated particle source is larger than $\mathrm{d} N/\mathrm{d}t$, indicating that transport plays a key role in the inter-ELM pedestal recovery.

fuelling↗
DOE OSTI2023

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Macroscopic trends of linear tearing stability in cylindrical current profiles

Abstract The likelihood of realising tokamak power-plants will be greatly improved by the discovery of high-gain equilibria that resist the formation of small islands and hence avoid the disruptive neoclassical tearing mode. We propose a series of studies to understand how simple tokamak design can leverage aspects of tearing onset physics to maximise passive resistance to island formation. Here we investigate the variation that current profiles can bring about in preventing tearing onset through the cylindrical linear tearing stability parameter Δ ′ . A database of 159148 realistic pilot-plant current profiles was generated with Monte Carlo sampling, and the distribution of Δ ′ values was linked with interpretable profile characteristics. In agreement with prior theoretical and experimental studies, Δ ′ was found to be strongly correlated with the existence and steepness of a local toroidal current well or hill, with the former destabilising and the latter stabilising. In the absence of these two cases, the remaining Δ ′ values were linearly bounded by the toroidal current gradient at the rational surface.

Physics↗
DOE OSTI2025

Technology readiness assessment of magnetohydrodynamic stability control

The technology readiness level scale is revised with clear definitions for magnetohydrodynamic stability control technologies for magnetic confinement fusion. The definitions are in the form of a sentence with three dimensions, which are the system complexity, the test environment, and the cost. A clear and simple definition is sought after to make an assessment based on the evidence which is agreeable to the fusion energy development community. A sample technology readiness assessment (TRA) for the electron cyclotron current drive-based control of neoclassical tearing mode for ST Advanced Reactor, using the proposed scale, reveals that the risk of a TRA being misinterpreted is high, which can only be avoided with an understanding of the procedure described in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗
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 OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2026

First use of an InSb crystal for x-ray imaging spectroscopy of highly ionized tungsten in the Wendelstein 7-X stellarator

Advances in time and space resolved measurements of highly charged states of tungsten (W) through x-ray imaging spectroscopy have enabled investigation of impurity transport in the Wendelstein 7-X (W7-X) stellarator. The high-resolution x-ray imaging spectrometer (HR-XIS) system on W7-X utilizes the Bragg diffraction properties of a set of multiple crystals to measure a range of impurity emission lines within sections of the 1–7 Å wavelength range, including transitions of W. A new indium–antimonide crystal has been installed on the HR-XIS system to allow viewing of the 5.5–6.2 Å region focusing on emissions of W. Consequently, many bright W emission lines from W 40+ to W 47+ were observed in this wavelength range, both in plasmas with injected W and in those with only intrinsic W impurity sources, showing the high sensitivity of the diagnostic. Three W 46+ emissions correspond in wavelength and intensity with calculated photon emissivity coefficients and can be exploited for W transport and concentration applications in plasmas with T e ≳ 2.1 keV. Here, an estimate of the core W density behavior in two separate turbulence-reduced ‘high-performance’ (HP) discharges on W7-X is done using the 5.6893 Å W 46+ line. The n W behavior in HP scenarios can be explained by previous experimental results and neoclassical predictions.

XICS↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
DOE OSTI2020

Cause and impact of low-frequency chirping modes in DIII-D hybrid discharges

Significant variations in MHD activity and fast-ion transport are observed in the DIIID high-beta, steady-state hybrid discharges with a mixture of Electron Cyclotron (EC) waves and Neutral Beam Injection (NBI). When Electron Cyclotron Heating (ECH) or Current Drive (ECCD) is applied, Alfven Eigenmodes (AEs) are usually suppressed and replaced by low-frequency bursting modes. The analysis of a recently compiled database of hybrid discharges suggests that the change of the fast-ion pressure especially the perpendicular pressure is the main factor responsible for the instability transition. The lower ratio of fast-ion injection speed v inj to Alfven speed v alfven and slight drop of safety factor q min during ECCD may also facilitate the transition. The database shows that AEs mainly occur when the fast-ion fraction P f /P total is less than 0.53 and v inj /v alfven is greater than 0.50, while low-frequency bursting modes appear in the opposite regime. In this research, P f and P total are the central fast-ion pressure from classical prediction and total plasma pressure, respectively. The reason why the instability transition correlates with P f /P total and v inj /v alfven is that they can significantly modify the drive of low-frequency bursting modes and AEs. The explanation is supported by the observation that low-frequency bursting modes are rarely seen in the hybrids with NBI only, with EC waves and counter-NBI, or with high plasma density. A careful check of the low-frequency bursting modes suggests that they are mainly chirping (neoclassical) tearing modes [referred to as chirping (N)TM], i.e. the mode frequency firstly jumps up from the steady (N)TM frequency, then chirps down, and finally returns to the steady (N)TM frequency. Sometimes, the (N)TMs are fully stabilized and replaced with pure fishbones. Lastly, the resonance condition calculation and ’Kick’ model simulations suggest that (N)TMs and fishbones can interact through modification of the fast ion distribution in phase space, which influences the drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗