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 127 records · Page 7

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗

Superconducting air core inductor systems and methods

Provided is a low-weight, high-efficiency inductor design for use with or in electrical power equipment, such as inverters. A toroidal power inductor includes a support structure comprising an outer shell, an inner shell, and one or more coolant channels formed therebetween, a plurality of conductors wrapped around and supported by an exterior surface of the outer shell, and an interior cavity substantially enclosed by the inner shell of the toroidal support structure. The plurality of conductors are configured to provide an inductance for the toroidal power inductor, and the one or more coolant channels are distributed beneath the exterior surface of the outer shell to cool the plurality of conductors. An air-core power inductor may implement the conductors using high-temperature superconducting (HTS) tapes cooled by cryogenic fluid flowing within the coolant channels.

42 ENGINEERING↗

Quasi-linear toroidal simulations of resonant magnetic perturbations in eight ITER H-mode scenarios

Abstract Both linear and quasi-linear aspects of the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for various H-mode scenarios in ITER, covering the pre-fusion power operation and the fusion power operation phases. Linear response computations for eight ITER scenarios, with varying plasma current and toroidal magnetic field, reveal that the best coil current phasing for controlling the type-I edge localized modes (ELMs) scales roughly linearly with the edge safety factor. The coil phasing is defined as the relative toroidal phase of the coil currents between different rows, for a given toroidal harmonic. Quasi-linear initial value simulation, which is the focus of the present study, shows that application of the n = 3 ( n is the toroidal mode number) RMP field has a minimum side effect on the plasma core momentum confinement but potentially a large effect on the global particle transport. Generally, the RMP field with the best (worst) coil phasing for ELM control produces the strongest (weakest) effect on the plasma edge flow and the overall density. This robustly holds for all eight ITER scenarios. Consequently, in order to minimize the RMP induced side effects while achieving ELM control (suppression) in ITER, a compromise is necessary in choosing the coil current configuration.

Physics↗

Divertor characterization and access to dissipative divertor conditions in negative triangularity discharges in DIII-D

Reduced radial transport, short midplane-to-target parallel connection lengths, and a strong effect of cross-field drifts were responsible for the high densities required for detachment in strong negative triangularity (NT) configurations in DIII-D. Dissipative divertor conditions were achieved in NT discharges at different triangularity, injected power, plasma current, and toroidal field direction. Differences between negative and positive triangularity (PT) discharges are analyzed in this paper to understand the requirements for access to detached divertor conditions: power balance, geometry, radial transport and effect of cross-field particle drifts. Parametric dependencies of access to detachment on plasma current and power flowing into the scrape-off layer remained similar in negative and PT and impurity seeding was observed to reduce the density needed to detach by up to 30% at the expense of core impurity dilution. The impact of triangularity on core-edge integration was tested varying bottom triangularity at fixed top triangularity. The high density needed to detach was not intrinsic to the NT edge as shapes with positive lower triangularity and negative upper triangularity were able to detach at lower upstream densities while maintaining an ELM-free NT edge. Confinement degradation at deeper detachment levels was however observed in all NT shapes, often associated with radiation instabilities.

detachment↗

Calculation of tearing mode stability in an inverse aspect-ratio expanded tokamak plasma equilibrium

The tearing mode stability of an inverse aspect-ratio expanded tokamak plasma equilibrium of general shape is investigated using asymptotic matching techniques. Particular emphasis is placed on the conservation of toroidal electromagnetic angular momentum. The TJ code, which is a specific implementation of the results of the investigation, is described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗

Analytic investigation of combined effects of anisotropic thermal transport and energetic particles on stability of resistive plasma resistive wall mode

The combined effects of anisotropic thermal transport and trapped energetic particles (EPs) on the stability of the resistive plasma resistive wall mode (RPRWM) are investigated by an energy-principle based analytical model. The results qualitatively confirm that of a recent toroidal modeling study [Bai et al., Phys. Plasmas 27, 072502 (2020)], in which the thermal transport can stabilize the RPPWM depending on the parameters of both the plasma equilibrium and energetic particles. The analytical model predicts a complete stabilization of the RPRWM in highly resistive plasmas, at sufficiently high EPs’ pressure, and finite plasma flow. The stabilizing effect of thermal transport originates from its enhancement of energy dissipations associated with both the resistive layer and the trapped energetic particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microturbulence Suppression by Alfvén Eigenmodes in the DIII-D Tokamak

Mitigation and suppression of low-𝑘 turbulence are observed during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in DIII-D experiments. Turbulence mitigation begins when the dominant TAE starts to depart from the typical shear Alfvén wave polarization. TAE activity then evolves into a new state, characterized by more discrete coherent modes, increased total amplitude, and more localized radial structures. During this transition, a narrow shear flow layer forms, driven by an enhanced Reynolds stress force, with a shearing rate that exceeds the local turbulence decorrelation rate, leading to a full turbulence suppression. Furthermore, these observations indicate the imbalance between Reynolds and Maxwell stress forces during the nonlinear evolution of the TAE and their important roles in shear flow generation and turbulence reduction.

Alfvén waves↗

Empirical probability and machine learning analysis of m , n = 2, 1 tearing mode onset parameter dependence in DIII-D H-mode scenarios

m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗

Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium

Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of divertor electrical drifts on particle distribution and detachment near the divertor target plate in DIII-D

Strong impacts of drifts on the divertor plasma in–out asymmetry and detachment are demonstrated in DIII-D with an open divertor configuration. For forward toroidal field, B T , i.e., with the ion B × ∇B drift toward the divertor, the particle flux to the inner divertor, as represented by the Langmuir probe measured ion saturation current (J sat ), exhibits a double peak structure, with electron temperature, lower at the inner target. Reversing the B T direction reverses both the radial and poloidal E × B flows, leading to a broad particle flux profile in the outboard scrape-off layer (SOL) with a similar double-peak structure to that observed at the inner target with forward B T . The correlation of a double peak structure with divertor temperature profiles confirms physical coupling between the drift flow and sheath boundary condition and their strong impact on divertor profiles. In addition, under reversed B T conditions, increasing the density flattens the target temperature profile. However, J sat remains high away from the strike point, rendering it difficult to achieve an “effective” detached plasma, i.e., with effective reduction in both peak heat flux and peak temperature (in the far SOL). In contrast, divertor detachment with a cold and flat temperature profile can be achieved at both target plates with the forward B T .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessing physics of ion temperature gradient turbulence via hierarchical reduced-model representations

In this work, the saturation physics of ion temperature gradient (ITG) turbulence is probed by studying how amplitudes and scalings with key parameters vary across a hierarchy of reduced models. The models derive from nonlinear fluid equations for toroidal ITG turbulence under approximations to the mode coupling interactions in wavenumber space and the representation of turbulent decorrelation. Mode coupling approximations include local-in-wavenumber treatments like the spectral density of flux in quasilinear theory, a truncation to three nonlinearly interacting waves, and the interactions in a cascade to high radial wavenumber mediated by a single zonal flow. Turbulent decorrelation treatments are based on the triplet correlation time with and without eddy damping. Model fidelity is assessed by the scalings and magnitudes of the squared amplitudes of unstable mode, stable mode, and zonal flow with respect to the flow-damping rate and temperature gradient. It is shown that all models reproduce fundamental scalings, provided they incorporate the coupling of unstable mode, stable mode, and zonal flow. Accurate amplitude prediction requires eddy damping in the triplet correlation time and proper representation of the zonal-flow drive by interactions associated with the radial wavenumber cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗

Modeling study of divertor particle and heat flux asymmetries for EAST H-mode discharges

The BOUT++ transport code is run to study the effects of plasma drifts on the divertor out-in asymmetries (DOIAs) of particle and heat fluxes and their decay widths for EAST lower single null H-mode discharges. The diamagnetic drift seems to have no effects on the DOIAs of total particle and heat fluxes due to its divergence-free nature. However, it could significantly increase the DOIAs of peak particle and heat fluxes and the flux decay widths. The E × B drift is found to induce a large plasma flow to the divertor region, enhancing the DOIAs of both total and peak particle and heat fluxes and the flux decay widths. Both the radial and poloidal components of the E × B drift are necessary in increasing the DOIAs, however, the radial E × B drift seems to play a more important role. The effects on the DOIAs caused by both diamagnetic and E × B drifts are reversed with the reverse of toroidal magnetic field. The heat flux decay width λ q and spreading width S q are important physical and engineering parameters for the divertors and could be obtained by fitting the heat flux profiles at divertor targets. The λ q at the outer target from the simulation case with all drifts could well match with the multi-machine scaling proposed by Eich and the DOIA of λ q is in reasonable agreement with the scaling proposed by Goldston.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗