Search NASA⌕ Search

SEARCH · Search NASA

Results for “toroidal flow”

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

Global gyrokinetic simulation with kinetic electron for collisionless damping of zonal flow in stellarators

Global gyrokinetic simulations with kinetic electrons for collisionless damping of zonal flows in LHD and W7-X stellarators show that the helical components of the equilibrium magnetic field responsible for helically trapped particles have significant impacts on zonal flow. Kinetic electrons reduce zonal flow residue and increase the frequency of low frequency oscillation (LFO). Here, the LFO is induced by dominant helical harmonics of magnetic field strength. Furthermore, linear toroidal coupling of multiple toroidal n-harmonics barely affects the zonal flows, but can generate long wavelength toroidal harmonics with the same toroidal number as the helical magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2022

The role of B T -dependent flows on W accumulation at the edge of the confined plasma

Abstract Near-separatrix impurity accumulation between the crown and the outer midplane of tokamaks is a common feature in results from codes such as SOLPS-ITER and DIVIMP; however, experimental evidence of accumulation has only recently been obtained and is reported here. The codes find that the poloidal distribution of impurity ions in the scrape-off layer (SOL) depends primarily on toroidal field ( B T )-dependent parallel flow patterns of the background plasma and the parallel ion temperature gradient (∇ ‖ T ion ) force. Experimentally, Mach probes used in L-mode plasmas with favorable (for H-mode access) B T measure fast ( M ∼ 0.3–0.5) inner-target-directed (ITD) background plasma flows at the crown of single-null discharges. This study reports a set of DIVIMP simulations for two similar H-mode discharges from the DIII-D W metal rings campaign differing primarily in B T -direction to assess the effect that fast ITD flows have on the distribution of W ions in the SOL. It is found that for imposed ITD flows of M = 0.3, W ions that otherwise accumulate due to the ∇ ‖ T ion -force are largely flushed out. It is also found that doubling the radial diffusion coefficient from 0.3 to 0.6 m 2 s −1 prevents accumulation due to rapid cross-field transport into the far-SOL, where background plasma flows drain W ions to the divertors. Far-SOL W distributions from DIVIMP are then used to specify input to the impurity transport code 3DLIM, which is used to interpretively model collector probe (CP) deposition patterns measured in the ‘wall-SOL’. It is demonstrated that the deposition patterns are consistent with the DIVIMP predictions of near-SOL accumulation for the unfavorable- B T direction, and little/no accumulation for the favorable- B T direction. The wall-SOL CPs have thus provided the first experimental evidence, albeit indirect, of near-SOL W accumulation—finding it occurs for the unfavorable- B T direction only. For the favorable- B T direction, fast flows can largely prevent accumulation from occurring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

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↗
DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Stability of a tokamak plasma with diffuse toroidal rotation

The present work considers the stability of a high- $\beta$ , large aspect ratio, circular plasma with diffuse profiles for the safety factor and the angular toroidal frequency (López & Guazzotto, Phys. Plasmas , vol. 24, 032501). An application of the Frieman–Rotenberg formalism results in a system of scalar eigenmode equations whose coupling is retained at the plasma–vacuum transition but is disregarded across the plasma column, which is a standard practice. The solution technique consists of a multidimensional shooting method for the poloidal harmonics; robust initial guesses are constructed by solving the dispersion relation in the static scenario with vanishing magnetic shear. Flow shear appears as a high- $\beta$ toroidal contribution, and we illustrate its destabilizing influence on $n=1$ external kink modes in the presence of ideal and resistive walls. Internal resonances are avoided by means of the selection of appropriate equilibrium parameters. The stabilizing influence of a finite positive average magnetic shear is also exemplified.

Physics↗
DOE OSTI2024

Transformer coupled toroidal wave-heated remote plasma sources operating in Ar/NF 3 mixtures

Remote plasmas are used in semiconductor device manufacturing as sources of radicals for chamber cleaning and isotropic etching. In these applications, large fluxes of neutral radicals (e.g. F, O, Cl, H) are desired with there being negligible fluxes of potentially damaging ions and photons. One remote plasma source (RPS) design employs toroidal, transformer coupling using ferrite cores to dissociate high flows of moderately high pressure (up to several Torr) electronegative gases. In this paper, results are discussed from a computational investigation of moderate pressure, toroidal transformer coupled RPS sustained in Ar and Ar/NF 3 mixtures. Operation of the RPS in 1 Torr (133 Pa) of argon with a power of 1.0 kW at 0.5 MHz and a single core produces a continuous toroidal plasma loop with current continuity being maintained dominantly by conduction current. Operation with dual cores introduces azimuthal asymmetries with local maxima in plasma density. Current continuity is maintained by a mix of conduction and displacement current. Operation in NF 3 for the same conditions produces essentially complete NF 3 dissociation. Electron depletion as a result of dissociative attachment of NF 3 and NF x fragments significantly alters the discharge topology, confining the electron density to the downstream portion of the source where the NFx density has been lowered by this dissociation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigations of plasma response associated with resonant magnetic perturbation fields using perturbation method in KSTAR H-mode plasmas

Abstract The plasma response associated with the resonant magnetic perturbation (RMP) field was investigated using the small edge perturbations induced by a modulated supersonic molecular beam injection (SMBI) in KSTAR. The modulated SMBI provides a time-varying perturbation of the plasma density source in the region just inside the last closed flux surface and a modulated flow damping rate. Radial propagation of the toroidal rotation perturbation induced by SMBI from the q = 3 surface to the q = 2 surface was observed. Theoretical analysis using the general perturbed equilibrium code of the RMP intensity profiles of the RMP field is consistent with the phase profile of the toroidal rotation perturbation.

Physics↗
DOE OSTI2022

Global gyrokinetic simulations of electrostatic microturbulent transport using kinetic electrons in LHD stellarator

Global gyrokinetic simulations of ion temperature gradient (ITG) and trapped electron mode (TEM) in the LHD stellarator are carried out using the gyrokinetic toroidal code (GTC) with kinetic electrons. ITG simulations show that kinetic electron effects increase the growth rate by more than 50% and more than double the turbulent transport levels compared with simulations using adiabatic electrons. Zonal flow dominates the saturation mechanism in the ITG turbulence. Nonlinear simulations of the TEM turbulence show that the main saturation mechanism is not the zonal flow but the inverse cascade of high to low toroidal harmonics. Further nonlinear simulations with various pressure profiles indicate that the ITG turbulence is more effective in driving heat conductivity whereas the TEM turbulence is more effective for particle diffusivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal modeling of 3D perturbations generated by current filaments in scrape-off layer in tokamak with biased divertor targets

The divertor biasing technique offers a promising alternative to control the edge localized mode (ELM) as well as the divertor heat load in tokamaks, as compared to the resonant magnetic perturbation (RMP) generated by magnetic coils. The linear resistive MHD code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is employed to study the plasma response to the n = 2 (n is the toroidal mode number) field perturbation in HL-2A, produced by the divertor biasing current filaments flowing in the scrape-off layer region. The response field due to biasing currents is compared with the RMP field utilized for ELM control in HL-2A. The strength of the radial resonant field perturbation, produced by the biasing currents of 100 A level, is found to be comparable with RMP generated by several kA of ELM control coil currents for the reference plasma. The plasma normal displacement near the X-point and the associated neoclassical toroidal viscosity torque are also computed to be similar between these two techniques. The modeling results thus strongly suggest that the biasing technique can be applied to control ELMs. Moreover, the biasing currents produce field perturbations, including the plasma response, that are localized more near the plasma edge (compared to the RMP counterpart), thus reducing the chance of mode locking associated with core perturbations. Particle orbit tracing also reveals that the biasing current produced magnetic perturbation tends to widen the heat deposition region and induce the strike point splitting of the ion saturation flow on the outer divertor surface, consistent with experimental observations in HL-2A. These toroidal modeling results confirm the possibility of ELM control and plasma exhaust solution by the divertor biasing technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Global gyrokinetic particle simulations of microturbulence in W7-X and LHD stellarators

Global gyrokinetic particle simulations of electrostatic ion temperature gradient (ITG) instability show that the most unstable eigenmode is localized to some magnetic fieldlines or discrete locations on the poloidal plane in the Wendelstein 7-X (W7-X) stellarator due to its mirror-like magnetic fields, which vary strongly in the toroidal direction and induce coupling of more toroidal harmonics (n) to form the linear eigenmode than in the Large Helical Device (LHD) stellarator. Nonlinear electrostatic simulation results show that self-generated zonal flows are the dominant saturation mechanism for the ITG instabilities in both the LHD and W7-X. Additionally, radial widths of the fluctuation intensity in both the LHD and W7-X are significantly broadened from the linear phase to the nonlinear phase due to turbulence spreading. Finally, nonlinear spectra in the W7-X are dominated by low-n harmonics, which can be generated both by nonlinear toroidal coupling of high-n harmonics and by linear toroidal coupling with large amplitude zonal flows due to the 3D equilibrium magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗