Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tokamak plasma physics”

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 469 records · Page 26

Comparative studies of cross-phase dynamics in turbulent momentum flux and particle flux at the tokamak edge

Turbulent transport events, including turbulent transport flux of momentum (i.e., turbulent momentum flux or Reynolds stress) and turbulent transport flux of particle (i.e., turbulent particle flux), have important effects on the confinement performance of magnetic confinement fusion devices. Poloidal Reynolds stress is the ensemble average of the product of radial velocity fluctuations and poloidal velocity fluctuations, i.e., $\langle {\widetilde{v}}_{r}{\widetilde{v}}_{\theta }\rangle$. Turbulent particle flux is the ensemble average of the product of radial velocity fluctuations and density fluctuations, i.e., $\langle \widetilde{n}{\widetilde{v}}_{r}\rangle$. Changes in either amplitude of fluctuations or cross phase between fluctuations can cause changes in turbulent transport. In this paper, cross-phase dynamics in the Reynolds stress and turbulent particle flux at the tokamak edge are studied in detail. Reynolds stress and turbulent particle flux are, respectively, written as the product of fluctuation amplitudes and an average cross-phase factor. The mathematical expressions of the average cross-phase factors are derived. The average cross-phase factors and the power spectra of cross phase are obtained using experimental measurement data. It is found that the cross-phase dynamics in Reynolds stress and particle flux are very different. Reynolds stress is found to be more sensitive to cross phase than particle flux is. In the strong $E\times B$ shear layer, spatial slips of cross phase lead to the obvious radial gradient of Reynolds stress. In the no/weak $E\times B$ shear region, the cross phase in Reynolds stress tends to lock. Here, phase locking refers to that the power spectra of phase tend to distribute around a fixed phase which does not change with radial position, while phase slip means that the power spectra of cross phase tend to distribute around a phase that varies with radial position. Phase slip or locking mainly describes the central phase weighted by the power spectra, while the phase scattering mainly describes the dispersion of the power spectrum distribution of the phase. The increased scattering of cross phase, which indicates the power spectra distribution of the phase is more dispersed, contributes to the decreased Reynolds stress for higher collisionality. The cross phase in particle flux tends to lock in both strong and no/weak shear regions. The degree of scattering of cross phase in the particle flux does not change obviously as collisionality increases. For higher collisionality, it is the increased density fluctuation amplitude rather than cross-phase dynamics that leads to the increased particle flux. The underlying physical mechanism that causes Reynolds stress and particle flux to exhibit different phase dynamics is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal Response of a Lithium Vapor Divertor to Cyclical Operation

The lithium vapor divertor concept is being developed as a method to achieve detached divertor conditions in a tokamak while minimizing impurity radiation losses from the core plasma. SOLPS-ITER modeling has previously been used to identify some of the geometric constraints and required lithium evaporation rate of a lithium vapor divertor in a medium-sized tokamak during steady-state operation. Here an updated conceptual design based on these operating requirements is introduced and the thermal response of the system is modeled during cyclical operation, consistent with operation in a short-pulse tokamak. Controllability of the temperature of the lithium capillary porous system (CPS) is achieved by adopting a design where there is no line-of-sight for radiation from the plasma to reach the heated CPS surface. Operational strategies to minimize the amount of lithium evaporated between plasma discharges while achieving steady evaporation rates during plasma discharges are discussed and modeled here. The optimal feedforward control strategy demonstrated in this work is to ramp up the temperature of the evaporator as quickly as possible immediately before a plasma discharge and then reduce the heating to match the desired steady-state net evaporation rate just before the plasma discharge begins, allowing the thermal inertia of the system to stabilize the evaporation rate during the first second of the plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Optical design and synthetic analysis of the electron cyclotron emission imaging diagnostic of HL-2M tokamak

The electron cyclotron emission imaging (ECEI) diagnostic is a powerful tool to study the MHD and turbulent transport in magnetically confined fusion plasmas. In this work, the optical system including the local oscillator (LO) coupling and radio frequency (RF) receiving optics, has been designed and analyzed for the HL-2M ECEI diagnostic. The LO optics can illuminate the antenna array and drive the mixer diode to work efficiently, with more than 36% of the beam intensity for the channels which probe the plasma edge relative to those of channels which probe the plasma core. The RF optics aims at guiding the plasma emission signal to the antenna array. To meet different physical requirements, three types of field of view have been achieved for the plasma imaging, with zoom factors of around 1, 1.5 and 2, respectively. The focal surfaces are almost flat, with a maximum off-mapping (defined as the radial distance between the beam waists of the lower-/uppermost antennas and electron cyclotron emission layer) less than 2.5 cm, which can match the electron cyclotron emission layer quite well and suppress the image distortion in the plasma edge region. Furthermore, the impact of possible vibrations and installation error on the focal plane has been evaluated. In addition, the predefined MHD and turbulence perturbations are well reproduced by combining the optical simulation results and synthetic ECEI modeling, which further verifies the good performance of the RF optics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c , that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c , due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic-ballooning-limited pedestals in spherical tokamak plasmas

Abstract A theoretical model is presented that for the first time matches experimental measurements of the pedestal width-height Diallo scaling in the low-aspect-ratio high- β tokamak NSTX. Combining linear gyrokinetics with self-consistent pedestal equilibrium variation, kinetic-ballooning, rather than ideal-ballooning plasma instability, is shown to limit achievable confinement in spherical tokamak pedestals. Simulations are used to find the novel Gyrokinetic Critical Pedestal constraint, which determines the steepest pressure profile a pedestal can sustain subject to gyrokinetic instability. Gyrokinetic width-height scaling expressions for NSTX pedestals with varying density and temperature profiles are obtained. These scalings for STs depart significantly from that of conventional aspect ratio tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anisotropic E × B shearing rate in a magnetic island

We derive the E × B shearing rate associated with vortex flow inside a macroscopic magnetic island (MI) in axisymmetric toroidal geometry. Due to the elongation of the MI and incompressibility of the E × B flow, the shearing rate near X-points is much lower than that near the mid-plane (x-axis of the local Cartesian coordinate) of the MI on the same flux surface. Furthermore, the rate formally vanishes at the X-points where the local poloidal magnetic field associated with the MI stagnates. Here, this calculation of E × B shearing profile and, in particular, minimal E × B shear near the X-points is consistent with the recent experimental finding that turbulence tends to spread into an MI through regions around the X-points and can contribute to more thorough quantitative interpretation of the results from experiments and simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of triangularity on ion-temperature-gradient-driven turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. Here, the effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ|≤0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulated effect of edge plasma density parameters on lower hybrid wave scattering in EAST

The incorporation of lower hybrid (LH) wave spectrum broadening in the poloidal wavenumber (⁠k θ ) space at the last close field surface (LCFS) is reported to lead to better agreement of the modeled LH wave current profile with the experimental results [Baek et al., Nucl. Fusion 61, 106034 (2021)]. To further understand its underlying mechanism and find the possible influence factors, effects of wave scattering caused by drift-wave type density fluctuation on the probability distribution of the LH wave polar refractive index (⁠N θ ) at the LCFS are studied under density parameters in the scrape-off-layer. According to a scattering model [P. T. Bonoli and E. Ott, Phys Fluids 25(2), 359–375 (1982)], scattering probability and scattering angle distribution are two main factors that determine the degree of spectral broadening. Studies presented here show that the total scattering probability increases first and then decreases as the wave propagates toward a smaller normalized radius of poloidal magnetic flux (⁠ρ⁠). The degree of spectral broadening is found to depend on the density and density fluctuation together by changing the intensity and a proportion of the geometrical optics approximation term and the E×B drift term in the scattering model. Additionally, the fluctuation correlation length can significantly modify the probability distribution of N θ at the LCFS, which is found to significantly change the LH wave current profile.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evidence on the effects of main-chamber neutrals on density shoulder broadening

Evidence that density shoulder broadening is dependent on high main-chamber neutral density is presented. Shoulder broadening does not occur when the sources for main-chamber neutrals are minimized using divertor baffles and wide gaps to the first wall (~3× the density decay length). Removing the baffles or reducing the gap to the inner wall both act to increase the density shoulder amplitude in otherwise identical TCV discharges. Radial turbulent transport is correlated with shoulder amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement and simulation of small cryogenic neon pellet Ne-I 640 nm photon efficiency during ablation in DIII-D plasma

Small (~1 mm) neon pellet fragments are fired into DIII-D H-mode plasmas, and resulting trajectory-averaged photon efficiencies (neutral neon ionization events for every photon emitted) of S/XB ≈ 85 are estimated for Ne-I 640 nm by dividing the estimated initial pellet fragment mass by the measured number of emitted Ne-I photons. The experiments are modeled by running the Lagrangian particle (LP) fluid/magneto-hydrodynamic pellet code to estimate axial ablation plume neon density profiles and temperature profiles at each pellet position. These solutions are then fed into the PrismSPECT collisional-radiative code, which calculates resulting neon charge states and photon emission rates, giving a profile-average of S/XB ≈ 109. The burnthrough plasma minor radius predicted by LP (ρ ≈ 0.63) is reasonably close to the experimental observation ρ ≈ 0.6. Furthermore, the modeling indicates that local S/XB is not constant along the pellet trajectory but tends to increase with increasing ablation rate. Non-equilibrium kinetics are predicted to be very important, while line trapping is predicted to be relatively unimportant (for Ne-I 640 nm S/XB).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Close coupling of the peeling ballooning and external kink modes in high poloidal-beta discharges with strong shape and large pedestal width

In plasmas with strong shape and large pedestal width, like DIII-D high poloidal-beta discharges, the eigenvalue spectrum of the peeling ballooning mode (PBM) is shown to shift to the very low-n regime near one, where n is the toroidal mode number. This peeling-type eigenvalue spectrum has a further shift to n = 1, thus being smoothly connected to the n = 1 external kink mode (EKM), as the normalized beta (β N ) increases. Once this connection occurs, the mode takes a mixture form of the PBM and EKM with its mode structure varying from the PBM-like to the EKM-like one as β N increases. The mode stability also becomes sensitive to both the local pedestal gradient and global β N , thus allowing an anti-correlation between the two driving forces. Further, these results appear to provide a qualitative explanation of the two unexpected features observed in the DIII-D high poloidal-beta discharges, that is, the dominance of the n = 1 mode in the edge-localized-modes and the negative correlation between the pedestal height and the internal transport barrier strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗