Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron temperature profile”

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 55 records · Page 3

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

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↗

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↗

Initial Results From the Newly Upgraded LTX- β

LTX-β, the upgraded Lithium Tokamak Experiment, recently began its first campaign with a goal to study the transport properties of gradient-free electron temperature profile equilibria with increased toroidal field and neutral beam injection. The temperature-gradient-free equilibria in LTX were enabled by the lithium plasma-facing surface. To make the lithium surface conditions more consistent for LTX-β, a new, faster lithium evaporator system has been developed. The toroidal magnetic field has been nearly doubled to ≥0.3 T. Auxiliary heating and fueling with a high flux (35 A, 20 kV) neutral beam has begun. The beam provides core fueling without cold edge neutrals, to test whether flat temperature profiles can be sustained, and makes possible core ion temperature measurements via charge exchange recombination spectroscopy. Neutral beam fueling has been demonstrated with increases in line-integrated density. The coupling of the beam matches well with NUBEAM predictions for shine through. Improved performance has been observed following lithium evaporative coatings and an increase in field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybrid simulations of sub-cyclotron compressional and global Alfven 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 (CAE) and global (GAE) Alfven Eigenmodes 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 power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, 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 simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

Compressional Alfven Eigenmodes↗

Real-time Electron Temperature and Density Profile Measurements for NSTX-U

This project set out to design, build, and deploy a real-time multi-point Thomson scattering (rt MPTS) system for NSTX-U, capable of providing electron temperature and density profiles—and their uncertainties—in real time for use in equilibrium reconstruction and plasma control. The long-term goal was to move Thomson scattering from a purely post-shot diagnostic to a central element of real-time situation awareness and feedback control on NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Public Data Set: Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection

This public data set contains openly-documented, machine readable digital research data corresponding to figures published in G.M. Bodner et al., ‘Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection,’ Physics of Plasmas 28, 102504 (2021) and its erratum in Physics of Plasmas 31, 129904 (2024).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Public Data Set: Erratum: “Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection” [Phys. Plasmas 28, 102504 (2021)]

This public data set contains openly-documented, machine readable digital research data corresponding to figures published in G.M. Bodner et al., ‘Erratum: “Initial Characterization of Electron Temperature and Density Profiles in PEGASUS Spherical Tokamak Discharges Driven Solely by Local Helicity Injection” [Phys. Plasmas 28, 102504 (2021)],’ Physics of Plasmas 31, 129904 (2024).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence for electron energization accompanying spontaneous formation of ion acceleration regions in expanding plasmas

We report experiments conducted in an expanding argon plasma generated in the inductive mode of a helicon source in the Hot hELIcon eXperiment–Large Experiment on Instabilities and Anisotropies facility. As the neutral gas pressure increases, the supersonic ion acceleration weakens. Increasing neutral pressure also alters the radial profile of electron temperature, density, and plasma potential upstream of the plasma expansion region. Langmuir probe measurements of the electron energy probability function (EEPF) show that heating of electrons at the plasma edge by RF fields diminishes with increasing gas pressure, yielding a plasma with a centrally peaked electron temperature, and flat potential profiles at higher neutral pressures. For neutral pressures at which ion acceleration regions develop in the expanding plasma plume, EEPFs reveal electrons with two temperature components.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simultaneous control of the electron temperature and safety factor profiles in DIII-D using model-based optimal control techniques

Future tokamak power plants will likely operate using a single, well-defined plasma scenario, either in steady state or for very long pulse lengths. In order to enhance the robustness of the scenario, feedback controllers for a variety of plasma properties will be necessary to counteract any disturbances and ensure safe operation. However, only a limited set of actuators will be available to control many different quantities. Because of this, it is necessary to develop controllers that are able to regulate multiple plasma properties using a limited set of actuators. To this end, a controller has been developed for the simultaneous regulation of both the electron temperature and safety factor profiles in DIII-D. This algorithm uses a linear quadratic integral control synthesis approach based on a linearized model of the dynamics of the two profiles. Two neural network surrogate models, NubeamNet and MMMnet, are included to improve the fidelity of the model. Furthermore, the controller has been tested in simulation using COTSIM, and has demonstrated the ability to simultaneously track changes in both the electron temperature and safety factor targets, including changes in both the magnitude and the shape of the profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian inference of electron density and ion temperature profiles from neutral beam and halo Balmer- α emission at Wendelstein 7-X

By employing Bayesian inference techniques, the full electron density profile from the plasma core to the edge of Wendelstein 7-X (W7-X) is inferred solely from neutral hydrogen beam and halo Balmer- α (H α ) emission data. The halo is a cloud of neutrals forming in the vicinity of the injected neutral beam due to multiple charge exchange reactions. W7-X is equipped with several neutral hydrogen beam heating sources and an H α spectroscopy system that views these sources from different angles and penetration depths in the plasma. As the beam and halo emission form complex spectra for each spatial point that are non-linearly dependent on the plasma density profile and other parameters, a complete model from the neutral beam injection and halo formation through to the spectroscopic measurements is required. The model is used here to infer electron density profiles for a range of common W7-X plasma scenarios. The inferred profiles show good agreement with profiles determined by the Thomson scattering and interferometry diagnostics across a broad range of absolute densities without any changes to the input or fitting parameters. The time evolution of the density profile in a discharge with continuous core density peaking is successfully reconstructed, demonstrating sufficient spatial resolution to infer strongly shaped profiles. Furthermore, it is shown as a proof of concept that the model is also able to infer the main ion temperature profile using the same data set.

forward modeling↗

Prediction of transport in the JET DTE2 discharges with TGLF and NEO models using the TGYRO transport code

Abstract The JET Deuterium-Tritium-Experiment Campaign 2 (DTE2) has demonstrated the highest-ever fusion energy production. To forecast the transport dynamics within these discharges, the TGLF and NEO models within the TGYRO transport code were employed. A critical development in this study is the new quasilinear transport model, TGLF-SAT2, specifically designed to resolve discrepancies identified in JET deuterium discharges. This model accurately describes the saturated three-dimensional (3D) fluctuation spectrum, aligning closely with a database of nonlinear CGYRO turbulence simulations, thereby enhancing the predictive accuracy of TGYRO simulations. In validating against the JET DTE2 discharges across two primary operating scenarios, TGYRO effectively predicted the temperature profiles within a broad radial window ( ρ ∼ 0.2–0.85), though with minor ion temperature discrepancies near the core. However, a consistent underprediction of electron density profiles by 20% across the simulation domain was noted, indicating areas for future refinement. To achieve a self-consistent steady-state solution based on the JET DTE2 discharges, an integrated modeling workflow TGYRO-STEP within the OMFIT framework was introduced. This workflow iterates among the core transport, the pedestal pressure and the MHD equilibrium, ultimately yielding a converged solution that significantly reduces dependence on experimental boundary conditions for temperature and density profiles. The integrated simulation results show negligible differences in electron density and temperature profiles compared to standalone TGYRO modeling, while the ion temperature profile is lower due to the updated boundary condition in TGYRO-STEP. The application of the TGYRO-STEP workflow to JET DTE2 discharges serves as a crucial test to validate its robustness and highlights its limitations, providing valuable insights for its potential future application in ITER and Fusion Power Plant deuterium and tritium prediction modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental study of high-k turbulence during an energy confinement degradation phase in EAST ohmic plasmas

In this paper, we present experimental studies of both high-k(r) ( cm(-1), and , respectively; ) and high- ( cm(-1), and , respectively; ) turbulence behavior during an ohmic energy confinement degradation phase in experimental advanced superconducting tokamak (EAST). High-k(r) turbulence from density fluctuation at and high- turbulence from density fluctuation at -0.97 were measured by tangential and poloidal CO2 laser collective scattering diagnostics, respectively. Note that k(r), , and are radial wavenumber, poloidal wavenumber, perpendicular wavenumber and ion gyroradius at electron temperature, respectively. Both high-k(r)/ turbulence power and energy confinement time are found to be temporally correlated to line-averaged electron density n(e) in the plasma current flat-top phase (I-p = 0.4 MA): when the n(e) shows continuous increase/decrease, the high-k(r)/ turbulence power increases/decreases correspondingly and the shows corresponding decrease/increase; the stable n(e) is related to stable both and high-k(r)/ turbulence power. Statistical results of high-k(r)/ turbulence power versus further imply that high-k(r)/ turbulence shows a strong correlation with plasma energy confinement degradation in high line-averaged n(e) condition, but high-k(r) and high- turbulence have relatively weak and no dependence on the transition between linear range in low-n(e) condition and energy confinement degradation in the high-n(e) condition, respectively. Moreover, profiles of electron temperature T-e and n(e) as well as their normalized gradients in a part of high- density fluctuation measurement region also have been given to qualitatively explain the enhancement of turbulence power with the increase of line-averaged n(e).

EAST tokamak↗

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗

Analysis of synthetic electron cyclotron emission from the high field side of HL-2M tokamak plasmas

A synthetic electron cyclotron emission (ECE) diagnostic is used to interpret ECE signals from preset plasma equilibrium profiles, including magnetic field, electron density, and electron temperature. According to the simulation results, the electron temperature (T e ) profile covering the harmonic overlap region can be obtained by receiving ECE signals at the high field side (HFS) of the HL-2M plasma. The third harmonic ECE at the low field side (LFS) cannot pass through the second harmonic resonance layer at the HFS unless the optical thickness (τ) of the second harmonic becomes gray (τ ≤ 2). In addition, the impact of the relativistic frequency down-shift has been evaluated and corrected. The measurable range of the HFS ECE has been calculated by scanning different parameters (electron density, temperature, and magnetic field). Higher plasma parameters allow a wider radial range of electron temperature measurements. The minimum inner measurable position can reach R = 120 cm (r/a = -0.89) when the product of core temperature (T e0 ) and density (n e0 ) is greater than 35 × 10 19 keV m -3 , which is extended by more than 30 cm inward compared with that of the LFS measurement. The HFS ECE will greatly improve the diagnostic ability of ECE systems on the HL-2M tokamak.

47 OTHER INSTRUMENTATION↗

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗