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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
NASA NTRS1997

Effects of Neutral Density on Energetic Ions Produced Near High-Current Hollow Cathodes

Energy distributions of ion current from high-current, xenon hollow cathodes, which are essential information to understand erosion phenomena observed in high-power ion thrusters, were obtained using an electrostatic energy analyzer (ESA). The effects of ambient pressure and external flow rate introduced immediately downstream of hollow cathode on ion currents with energies greater than that associated with the cathode-to-anode potential difference were investigated. The results were analyzed to determine the changes in the magnitudes of ion currents to the ESA at various energies. Either increasing the ambient pressure or adding external flow induces an increase in the distribution of ion currents with moderate energies (epsilon less than 25 to 35 eV) and a decrease in the distribution for high energies (epsilon greater than 25 to 35 eV). The magnitude of the current distribution increase in the moderate energy range is greater for a cathode equipped with a toroidal keeper than for one without a keeper, but the distribution in the high energy range does not seem to be affected by a keeper. An MHD model, which has been proposed to describe energetic-ion production mechanism in hollow cathode at high discharge currents, was developed to describe these effects. The results show, however, that this model involves no mechanism by which a significant increase of ion current could occur at any energy. It was found, on the other hand, that the potential-hill model of energetic ion production, which assumes existence of a local maximum of plasma potential, could explain combined increases in the currents of ions with moderate energies and decreases in high energy ions due to increased neutral atom density using a charge-exchange mechanism. The existing, simplified version of the potential-hill model, however, shows poor quantitative agreement with measured ion-current-energy-distribution changes induced by neutral density changes.

Kameyama, Ikuya↗
DOE OSTI2025

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

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

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

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↗
NASA NTRS1995

Low-state disks and low-beta disks

Stellar black hole candidates (BHCs) exhibit bimodal spectral states. We calculate nonthermal disk spectra, demonstrating that a large photon index (alpha (sub x) approximately 2-3) observed in the soft (high) state is due to a copious soft photon supply, whereas soft photon starvation leads to a smaller index (alpha (sub x) approximately 1.5-2) in the hard (low) state. Thus, the absence of the soft component flux in the low state cannot be due to obscuration. A possible disk configuration during the low state is discussed. We proposed that a low-state disk may be a low-beta disk in which magnetic pressure may exceed gas pressure becuase of the suppression of field escape by a strong shear. As a result, disk material will take the form of blobs constricted by mainly toroidal magnetic fields. Fields are dissipated mainly by occasional reconnection events with a huge energy release. This will account for large-amplitude, aperiodic X-ray variations (flickering) and high-energy radiation with small alpha(sub x) from hard state BHCs and possibly from active galactic nuclei. Further, we propose a hysteretic relation between the mass-flow rate and plasma-beta, a ratio of gas pressure to magnetic pressure, for the spectral evolution of transient BHCs. The disk is in the low-beta state in quiescence and early rise. The low-beta disk is optically thin and affected by advection. A hard-to-soft transition occurs before the peak luminosity, since there is no advection-dominated branch at higher luminosities. An optically thick, high-beta disk appears at small radii. In the decay phase of the light curve, the standard-type disk becomes effectively optically thin, when a soft-hard transition is triggered. High-beta plasmas in the main body shrink to form minute blobs, and low-beta coronal plasma fills interblob space.

Mineshige, Shin↗
DOE OSTI2025

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

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

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↗
NASA NTRS1984

Turbulent swirling combustion

The non-reacting flow produced by two confined, co- and/or counter-swirling jets wre measured by means of a two-color laser Doppler velocimeter. These results are compared with earlier experiments and with numerical predictions. It is shown that under both swirl conditions, a closed recirculation zone is created at the combustor center line. This zone is characterized by the presence of a one-cell toroidal vortex, low tangential velocities, high turbulent intensity, and large dissipation rates of the kinetic energy of turbulence.

Sommer, H. T.↗
DOE OSTI2026

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

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↗
NASA NTRS2000

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗
NASA NTRS1979

Fine particulate capture device

To capture fine particulate matter in a gas such as air, a dielectric fluid is directed to the center of whichever face of a rotating disc is exposed to the air flow. The disc is comprised of two or more segments which bear opposite electrostatic potentials. As the dielectric fluid is centrifuged towards the periphery of the rotating disc, the fluid becomes charged to the same potential as the segment over which it is passing. Particulate matter is attracted to the charged segment and is captured by the fluid. The fluid then carries the captured particulate matter to a collection device such as a toroidal container disposed around the periphery of the disc. A grounded electrically-conductive ring may be disposed at the outer periphery of the disc to neutralize the captured particles and the fluid before they enter the container.

Peterson, V. S.↗
DOE OSTI2021

Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B → . The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e. the temperatures are assumed to be isotropic. Additionally, in certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the two-dimensional (2D) transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗