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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The generalized perturbed equilibrium code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. Lastly, this work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Error field identification through torque balance on a saturated island in DIII-D

Measurement of the electromagnetic torque on a magnetic island could be an attractive method for error field identification in the early phase of ITER operation. Previous DIII-D experiments (Strait 2014 Nucl. Fusion 54 073004; Shiraki et al 2015 Plasma Phys. Control. Fusion 57 025016) have demonstrated the principle of this approach using a stationary or slowly rotating island, while recent developments in magnetic data analysis (Sweeney and Strait 2019 Phys. Plasmas 26 012509) allow the field of a rapidly rotating island to be readily distinguished from that of the wall currents induced by its rotation. In a recent experiment, a rotating n = 1 magnetic perturbation forced a saturated magnetic island to rotate, thus sampling all toroidal phases periodically in a single discharge. The phase and amplitude of the error field are inferred from analysis of the time-dependent torque balance on the island, including torques from the error field, the applied magnetic perturbation, and the wall currents induced by rotation of the applied perturbation and the island. Furthermore, results agree well with those from more conventional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Conceptual design of ELM control coils for the TCABR tokamak

An upgrade of the Tokamak à Chauffage Alfvén Brésilien (TCABR) is being designed to make it capable of creating a well controlled environment where the impact of resonant magnetic perturbation (RMP) fields on edge localised modes can be addressed over a wide range of (i) plasma shapes, (ii) divertor configurations, (iii) RMP coil geometries and (iv) perturbed magnetic field spectra. To this end, a unique set of in-vessel RMP coils is being designed and, in this work, their conceptual design is presented. This unique set of coils is composed of three toroidal arrays of coils on the low field side and three toroidal arrays of coils on the high field side. Each of these six toroidal arrays is composed of 18 coils thus allowing for the creation of RMP fields with toroidal mode numbers n ≤ 9 and with increased control of the poloidal mode number spectrum. To study dynamical effects of RMP fields of different toroidal mode numbers, all rotating simultaneously with different velocities, each of the 108 RMP coils will be powered independently by power supplies that can provide voltages of up to 4 kV and electric currents of up to 2 kA, with frequencies varying continuously from 0Hz up to 10kHz. A set of physical criteria were used to determine the optimal coil geometry and their respective number of turns to reduce the coil currents and voltages during operation with alternate current. Further, the conceptual design was carried out using both the vacuum approach (no plasma response) and the single-fluid response approach, which accounts for the response of a linear, single-fluid, visco-resistive plasma calculated using the M3D-C 1 code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D Toroidal Field Multipoles for Curved Accelerator Magnets

Curved magnets producing continuously rotating field multipoles along the length of the bend can provide strong and continuous transverse focusing, making them of interest for accelerator systems such as compact ion beam gantries and synchrotron light sources. Evaluating the utility of such rotating multipole systems requires an accurate description of field behavior for beam physics calculations. This paper presents a helical scalar potential solution in 3D toroidal harmonics relevant to boundary conditions on the surface of a torus. The resulting fields are evaluated for a curved helical quadrupole channel to illustrate field rotation and the effect of magnet curvature.

Brouwer, L.↗
DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Integrated control of edge localized modes and divertor flux using mixed toroidal harmonic resonant magnetic perturbations in EAST

Mixed harmonic resonant magnetic perturbations (RMPs) for integrated edge localized modes (ELMs) and divertor flux control are demonstrated in EAST target plasmas of low input torque and normalized beta β N ~1.7–1.9, which are close to the equivalent value in ITER high Q operation. The applied RMPs are designed to combine a static harmonic of the toroidal mode number n = 3 with a static or rotating harmonic of n = 2. ELM suppression is achieved without a drop of plasma energy confinement, and tungsten concentration is effectively reduced during the application of RMPs. With mixed harmonics, the toroidal varying steady state heat and particle fluxes on the divertor target can be modified with the rotating n = 2 harmonic, which agrees with the numerical modeling of three-dimensional magnetic topology, with plasma responses being taken into account. ELM suppression correlates with the times of larger n = 3 response with mixed n = 2 and n = 3 RMPs. The mixture of harmonics and the rotating n = 2 harmonic does not require additional coil current because the variation is only in the upper-lower coil current phase space. Furthermore, these results further affirm the effectiveness of integrated ELM and divertor flux control using RMPs with mixed harmonics and improve the understanding of the role of plasma responses in ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MARS-Q modeling of low- n resistive kink-peeling modes and edge harmonic oscillations in DIII-D

Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low-n (n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low-n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

EHO↗
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2024

Magnetic toroidicity

Abstract Directional non-reciprocity refers to the phenomenon where the motion in one direction differs from the motion in the opposite direction. This behavior is observed across various systems, such as one-way traffic and materials displaying electronic/optical directional dichroism, characterized by the symmetry of velocity vectors. Magnetic toroidal moments (MTMs), which typically arise from rotational spin arrangements, also possess the symmetry of velocity vectors, making them inherently directionally non-reciprocal. In this paper, we examine magnetic point groups (MPGs) that exhibit MTMs, subsequently leading to off-diagonal linear magnetoelectricity. Our focus is on the induction of MTMs through electric fields, magnetic fields, or shear stress, while enumerating the relevant MPGs. The findings of our study will serve as valuable guidance for future investigations on directional non-reciprocity, MTMs, and off-diagonal linear magnetoelectric effects.

Physics↗
DOE OSTI2024

Modeling of frequency-sweeping Alfvén modes in the TJ-II stellarator

Alfvénic activity has been observed in the TJ-II stellarator which resembles the frequency sweeping demonstrated by Alfvén cascade modes in tokamaks. A numerical validation study was conducted using a reduced magnetohydrodynamic (MHD) model to show that such modes could only have been observed in discharges where the rotational transform profile was non-monotonic. During experiments, coil current was varied which resulted in shifting of the minimum value of the rotational transform profile. To mimic this effect, we study the Alfvénic activity predicted by the reduced MHD model for a set of input rotational transform profiles with varying minima. A mode is found whose toroidal and poloidal mode numbers match those predicted in experiments which sweeps downward/upward in frequency as the minimum value of the rotational transform profile is increased/decreased. The results serve as a demonstration of the validity and utility of MHD spectroscopy.

Alfvén cascades↗
DOE OSTI2021

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗
DOE OSTI2026

Predicting core transport in ITER baseline discharges with neon injections

Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from a systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \unicode{x2A7D} Z_\textrm{eff} \unicode{x2A7D} 2.5$, and the corresponding power crossing the separatrix $P_\textrm{sep}$ is evaluated. We find that $P_\textrm{sep}$ varies by more than a factor of 1.7 across this scan and matches the ${\sim}100$ MW SOLPS-ITER prediction when $Z_\textrm{eff} \simeq 1.6$ or when auxiliary heating is reduced to ${\sim}75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_\textrm{sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_\textrm{eff} \approx 1.6$ –1.75 and $0.75 \lesssim f_{P_\textrm{aux}} \unicode{x2A7D} 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.

ITER↗
DOE OSTI2021

Optimizing the differential connection schemes for detecting 3D magnetic perturbations in DIII-D

Optimizing the differential pair connection scheme (i.e. the set of pairs) of a toroidal array of magnetic sensors dedicated to measuring slowly rotating asymmetric fields can enhance the mode number detection capability and failure-resilience. In this work the condition number obtained from singular value decomposition (SVD) of the design matrix is used as a metric to evaluate the quality of a connection scheme. A large number of possible pair connections are usually available, so evaluating all of them may require extensive use of computational resources and can be very time consuming. Alternative methods to reduce the number of pairs evaluated without losing the capabilities of toroidal mode detection are presented in this paper. Here, three examples of the applications of such analysis for the 3D magnetic diagnostic system of DIII-D are also presented: the addition of two new toroidal arrays with n > 3 detection capabilities, the modification of an existing toroidal array in the low field side of the machine to accommodate the addition of a helicon antenna, and the design of changes in several toroidal arrays in the high field side to accommodate the addition of a lower hybrid current drive antenna on the center post.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Halo current rotation scaling in post-disruption plasmas

Abstract Halo current (HC) rotation during disruptions can be potentially dangerous if resonant with the structures surrounding a tokamak plasma. We propose a drift-frequency-based scaling law for the rotation frequency of the asymmetric component of the HC as a function of toroidal field strength and plasma minor radius ( f rot ∝ 1/ B T a 2 ). This scaling law is consistent with results reported for many tokamaks and is motivated by the faster HC rotation observed in the HBT-EP tokamak. Projection of the rotation frequency to ITER and SPARC parameters suggest the asymmetric HC rotation will be on the order of 10 Hz and 60 Hz, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An extended variational method for the resistive wall mode in toroidal plasma confinement devices

The external-kink stability of a toroidal plasma surrounded by a rigid resistive wall is investigated. The well-known analysis of Haney and Freidberg is rigorously extended to allow for a wall that is sufficiently thick that the thin-shell approximation does not necessarily hold. A generalized Haney–Freidberg formula for the growth-rate of the resistive wall mode is obtained. Thick-wall effects do not change the marginal stability point of the mode but introduce an interesting asymmetry between growing and decaying modes. Growing modes have growth-rates that exceed those predicted by the original Haney–Freidberg formula. On the other hand, decaying modes have decay-rates that are less than those predicted by the original formula. The well-known Hu–Betti formula for the rotational stabilization of the resistive wall mode is also generalized to take thick-wall effects into account. Increasing wall thickness facilitates the rotational stabilization of the mode, because it decreases the critical toroidal electromagnetic torque that the wall must exert on the plasma. On the other hand, the real frequency of the mode at the marginal stability point increases with increasing wall thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗