Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal rotation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhancement of edge turbulence concomitant with ELM suppression during boron powder injection in EAST

A reproducible, quasi-stationary edge localized mode (ELM)-suppressed scenario was obtained over a wide range of plasma parameters by continuous injection of boron (B) powder into an upper-single null discharge in the experimental advanced superconducting tokamak [Sun et al., Nucl. Fusion 61, 014002 (2021)]. This powder-induced ELM-absent regime is associated with an edge harmonic mode (EHM) that provides continuous particle exhaust to maintain constant density without confinement degradation and impurity accumulation, the latter of which is often observed in ELM-free regimes. A flow rate threshold of B powder injection, leading to a threshold intensity of the EHM, is necessary for full ELM suppression. The fundamental harmonic of the EHM exhibits a toroidal mode number n = 1. The mode is observable in the entire poloidal cross section with a peak near the upper X-point in an upper-single null configuration. Here, the EHM spans radially across the pedestal and scrape-off layer, peaking inside the separatrix. The EHM appears to be insensitive to q 95 , heating power, plasma toroidal rotation, and pedestal collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurement of near- and far-field impurity flows during pellet-induced rapid shutdown in DIII-D

Both near-field (<1 m away toroidally from the pellet) and far-field (>1 m away toroidally from the pellet) poloidal and toroidal impurity (carbon ion) flows are measured using visible imaging and fast bolometry during a polypropylene pellet-induced rapid plasma shutdown in the DIII-D tokamak. In the near field, the pellet appears to increase poloidal flow in the ion diamagnetic direction, possibly due to the strong radial temperature gradient caused by the pellet ablation. In the far field, the poloidal impurity flow typically appears slower and in the opposite direction. Toroidal impurity flow appears to be strongly influenced by the plasma initial toroidal rotation, especially in the far field. These results demonstrate that rapid shutdown impurity flows are not necessarily global in structure but can be quite different close to and far from the injected pellet.

Electromagnetic radiation detectors↗
DOE OSTI2020

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of mean E × B shear in a stochastic magnetic field: ambipolarity breaking and radial current

The mean E $ \times $ B shear in a stochastic magnetic field is calculated, using the radial force balance relation and transport equations. This analysis is relevant to the L → H transition with resonant magnetic perturbations, and special focus is placed upon the physics of non-ambipolar transport and radial current. The key physical process is the flow of fluctuating currents along wandering magnetic fields. The increments in poloidal and toroidal rotation, density and ion pressure are calculated. The radial envelope of the magnetic perturbations inside the plasma defines a new scale ${\ell _{{\text{env}}}}$, which is the characteristic scale of the magnetic fluctuation intensity profile. Here, the net particle outflow due to stochastic magnetic fields is calculated and is determined by the net radial current through the separatrix. Implications for the L → H transition are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗
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 OSTI2025

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗
DOE OSTI2021

Visible spectroscopy of highly charged tungsten ions with the JET charge exchange diagnostic

Since the installation of the ITER-like wall on the JET tokamak, visible emission from a broad range of tungsten charge states has complicated plasma ion temperature and toroidal rotation measurements. A plethora of charge exchange transitions, from ions up to W 56+ , and 21 suspected magnetic dipole transitions have been observed. In particular W 39+ and W 46+ charge exchange emission contaminates the long-established carbon measurement at 529 nm. Fortunately the wavelength and relative intensity of the tungsten charge exchange lines is predictable and their influence can be mitigated when they are included in the spectral fit. Neon based charge exchange measurements at 525 nm are now preferred on JET since the ITER-like wall, however in this case an unidentified tungsten magnetic dipole transition contaminates the spectrum when the electron temperature is low.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Role of 3D neoclassical particle flux in density pump-out during ELM control by RMP in DIII-D

Application of three-dimensional (3-D) resonant magnetic perturbations, for the purpose of controlling edge localized modes (ELMs) in high confinement (H-mode) tokamak fusion operations, routinely leads to significantly reduction of the plasma density during the discharge, a phenomenon termed 'density pump-out'. Understanding the density pump-out phenomenon due to 3-D fields is a longstanding challenging issue. This work discusses a self-consistent toroidal model that allows initial value simulation of the density pump-out effect. The model successfully reproduces pump-out phenomenology in DIII-D ELM suppression experiments. The radial particle flux, associated with the 3-D neoclassical effects (the neoclassical toroidal viscosity, NTV), is discovered to play an important role in producing the observed large density pump- out. The key physics is associated with the drift kinetic Landau resonances between the toroidal precession of trapped thermal ions and electrons from one side, and the magnetic field perturbation from the other side. Such resonances, occurring near the top of the pedestal region, substantially enhances the radial particle flux associated with NTV. As a general consequence of the resonant NTV flux, larger density pump- out occurs in plasmas with lower collisionality and slower toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with ${\beta _{\rm{N}}}$ $ \geqslant $ 2.4 and ${H_{89}}$ $ \geqslant $ 2.0 at ${q_{95}}$ < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of ${\beta _{\rm{N}}} \lesssim 3.0$, ${H_{89}} \lesssim 2.4$ and G-factor ($ \equiv {\beta _{\rm{N}}}{H_{89}}/q_{95}^2$) $ \lesssim 0.4$6 has been achieved and sustained for $ \gtrsim 40{{\tau }_{\text{E}}}$ at ne/nGW ~ 0.7 with heating power of $ \lesssim $5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ${\it{\omega }_{E \times B}}$, and the magnetic shear effect. The electromagnetic finite $\beta $ stabilisation plays a role in the inner core region at ${\rho _{\rm{tor}}}\sim0.35$ together with the fast ion effect. The alpha stabilisation effect is also found at ${\rho _{\rm{tor}}}\sim0.5$. ${\omega _{E \times B}}$, which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ${\rho _{\rm{tor}}}\sim0.5 - 0.7$ with the highest contribution from the toroidal rotation. Furthermore, regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Divertor detachment in the pre-fusion power operation phase in ITER during application of resonant magnetic perturbations

Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗