SEARCH · Search NASA

Results for “toroidal”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Overview of recent experimental results on the EAST Tokamak

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (~1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high β P has been significantly extended (H 98y2 > 1.3, β P ~ 4.0, β N ~ 2.4 and n e /n GW ~ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q 95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m -1 ) with EC pre-ionization is developed. A beneficial role on the lower hybrid wave injection to control the tungsten concentration in the NBI discharge is observed for the first time in EAST suggesting a potential way toward steady-state H-mode NBI operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Role of electrostatic perturbation on kinetic resistive wall mode with application to spherical tokamak

Abstract A more complete non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid formulation is developed by including the perturbed electrostatic potential δφ in the particle Lagrangian. The fluid-like counter-parts of the hybrid equations, in the Chew-Goldberger-Low high-frequency limit, are also derived and utilized to test the new toroidal implementation in the MARS-K code. Application of the updated non-perturbative hybrid model for a high- β spherical tokamak plasma in MAST finds that the perturbed electrostatic potential generally plays a minor role in the n = 1 ( n is the toroidal mode number) resistive wall mode instability. The effect of δφ is largely destabilizing, with the growth rate of the instability increased by several (up to 20) percent as compared to the case without including δφ . A similar relative change is also obtained for the kinetic-induced resonant field amplification effect at high- β in the MAST plasma considered. The updated capability of the MARS-K code allows quantitative exploration of drift kinetic effects on various MHD instabilities and the antenna-driven plasma response where the electrostatic perturbation, coupled to magnetic perturbations, may play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Active control of Alfvén eigenmodes by external magnetic perturbations with different spatial spectra

Alfvén eigenmodes have been suppressed and excited in tokamak plasmas by (just) modifying the poloidal spectra of externally applied static magnetic perturbations. This effect is observed experimentally when toroidal spectra of n = 2, n = 4 as well as a mixed spectrum of n = 2 and n = 4 is applied. Under the n = 2 magnetic perturbations, the modes are excited or suppressed by modifying the coil phasing between the upper and the lower set of coils. Regardless of the absolute rotation, an even parity for the n = 4 perturbation is observed to reduce the amplitude of the Alfvénic instabilities, while an odd parity amplifies it. To combine the stabilizing (and destabilizing) effect of n = 2 and n = 4, a mixed spectrum is applied, finding similar reduction (and amplification) trends. However, the impact on the mode amplitude is more subtle, due to the reduced coil current required for a mixed spectrum. The signal level on the fast-ion loss detector is sensitive to the applied poloidal spectrum, which is consistent with Hamiltonian full-orbit modelling of an edge resonant transport layer activated by the 3D perturbative fields. An internal redistribution of the fast-ion population is induced, modifying the phase-space gradients driving the Alfvénic instabilities, and ultimately determining their existence. The calculated edge resonant layers for both n = 2 and n = 4 toroidal spectra are consistent with the observed suppressed and excited phases. Moreover, hybrid kinetic-magnetohydrodynamic (MHD) simulations reveal that this edge resonant transport layer overlaps in phase-space with the population responsible for the fast-ion drive. The results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

Alfvén waves↗

DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Plasma performance and operational space with an RMP-ELM suppressed edge

Abstract The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3 × 10 19 m −3 . This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Velocity profile shapes in Alcator C-Mod plasmas

Abstract Toroidal rotation velocity spatial profiles ( r / a < 0.8) have been obtained from C-Mod over a wide range of operational conditions, including H-mode, I-mode, ICRF-heated L-mode and Ohmic L-mode (LOC and SOC), and in plasmas with ITBs, LH wave injection and MCFD. Peaked, flat and hollow rotation profiles have been observed. In H- and I-mode plasmas, generally with co-current peaked profiles, the peaking is correlated with temperature profile peaking, and both increase with toroidal magnetic field (decrease with ρ ∗ ). Any dependence on density peaking is unclear. For Ohmic L-mode discharges, with LOC, the velocity profiles are usually flat and most often directed co-current, while with SOC the profiles are hollow, mostly co-current at the edge and counter-current in the core. Both of these Ohmic rotation states exist with matched density and temperature profiles (and gradients), indicating that neither gradient is relevant during rotation reversals. For plasmas with LH wave injection and discharges with ITBs, the velocity profiles are hollow while the density and temperature profiles exhibit substantial peaking. Broadly speaking for all operational regimes, there is no unifying ordering of the velocity gradient with plasma parameters.

Rice, J. E. (ORCID:0000000183195971)↗

DOE OSTI2024

Invariant regimes of Spencer scaling law for magnetic compression of rotating FRC plasma

Abstract The scaling laws for the magnetic compression of a toroidally rotating field reversed configuration (FRC) have been investigated in this work. The magnetohydrodynamics (MHD) simulations of the magnetic compression on rotating FRCs employing the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355), are compared with the Spencer’s one-dimensional (1D) theory (Spencer et al 1983 Phys. Fluids 26 1564) for a wide range of initial flow speeds and profiles. The toroidal flow can influence the scalings directly through the alteration of the compressional work as also evidenced in the 1D adiabatic model, and indirectly by reshaping the initial equilibrium. However, in comparison to the static initial FRC equilibrium cases, the pressure and the radius scalings remain invariant for the magnetic compression ratio B w 2 / B w 1 up to 6 in presence of the initial equilibrium flow, suggesting a broader applicable regime of the Spencer scaling law for FRC magnetic compression. The invariant scaling has been proven a natural consequence of the conservation of angular momentum of both fluid and magnetic field during the dynamic compression process.

Ma, Yiming↗

DOE OSTI2024

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗

DOE OSTI2024

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Gyrokinetic simulations of the effects of magnetic islands on microturbulence in KSTAR

Gyrokinetic simulations are utilized to study effects of magnetic islands (MIs) on the ion temperature gradient (ITG) turbulence in the KSTAR tokamak with resonant magnetic perturbations. Simulations show that the transport is controlled by the nonlinear interactions between the ITG turbulence and self-generated vortex flows and zonal flows, leading to an anisotropic structure of fluctuation and transport on the poloidal plane and in the toroidal direction. MIs greatly enhance turbulent transport of both particle and heat. The turbulent transport exhibits variations in the toroidal direction, with transport through the resonant layer near the island X-point being enhanced when the X-point is located at the outer mid-plane. A quantitative agreement is shown between simulations and KSTAR experiments in terms of time frequency and perpendicular wavevector spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Unified ELM Suppression on KSTAR and DIII-D via Adaptive Feedback Control Strategies

This paper reports on the extension of our amplitude-based resonant magnetic perturbation (RMP) edge localized mode (ELM) controller to support phasing control (relative toroidal phases of RMP waveforms between rows of coils), multiple toroidal mode numbers, and new ‘jump’ and ‘probing’ strategies, all deployed on KSTAR and DIII-D. By treating the control algorithm as device-independent and adjusting only the real-time interfaces to sensors and power supplies, we have confirmed that the same finite state machine—based feedback logic can be ported between machines with minor modifications. In experiments using n = 2 RMPs on KSTAR and n = 3 on DIII-D, the controller successfully modulated RMP amplitudes in real time to sustain ELM suppression while minimizing confinement degradation. Phasing control broadened the suppression window, as it permitted the system to avoid locked-mode regions and safely access ELM-free conditions. A rotating RMP phasing scheme, integrated into the same framework, distributes divertor heat loads more uniformly, making it a promising strategy for protecting plasma-facing components during long discharges. New ‘jump’ and ‘probing’ techniques demonstrate the possibility for the controller to preempt imminent ELMs and refine the minimum required RMP amplitude without returning to ELMy conditions. Taken together, these upgrades enable extended ELM-free operation while mitigating confinement degradation, providing a practical framework for real-time ELM control in future high-performance tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Observation of stationary filaments with resonant magnetic perturbations in NSTX

The application of resonant magnetic perturbations (RMPs) induces striated heat fluxes on the divertor plate in tokamaks. These striated heat fluxes were previously considered strike point splitting by the intersection between the divertor and predicted magnetic lobe structures. Recently, the existence of stationary 3D filaments has been observed for the first time in National Spherical Torus Experiment (NSTX) plasmas with toroidal mode number n = 3 RMPs. Magnetic field line tracing of these filaments from the midplane to the divertor is consistent with the quantity and the radial and toroidal location of striated heat fluxes observed on the divertor. This study suggests that divertor striated heat flux patterns are linked to these stationary filaments rather than to strike point splitting or magnetic lobe structures for NSTX. This finding will help guide future efforts to understand the heat flux striations in NSTX with RMPs in order to control and reduce the divertor peak heat flux.

RMPs↗

DOE OSTI2025

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗

DOE OSTI2025

Nonlinear Alfvén instability simulation and EP transport for ITER reversed shear (steady-state) and monotonic q-profile regimes *

The nonlinear evolution and dynamics of energetic particle (EP) driven instabilities in both steady-state (reversed shear q-profile) and monotonic q-profile regimes of ITER are examined using the global gyro-Landau closure model FAR3d. Both neutral beam and alpha components are included, leading to synergistic effects between the two populations. In order to maintain computational feasibility, the present nonlinear simulation includes all toroidal mode numbers ranging from n = 0 to 15. It cannot be excluded that n’s above this range could introduce larger linear growth rates; however, as indicated from the results, the n = 0 to 15 range is sufficient to provide a strong linear drive and significant nonlinear transport effects. While the toroidal mode numbers n = 6,13 dominate the linear instability growth phase, nonlinear energy transfers in the saturated phase reverse this trend, leading to the dominance of lower n’s in the saturated phase of the simulation. Zonal flow structure generation and EP density flattening are further consequences of the nonlinear phase. The density profile flattening is caused by collectively driven EP radial transport fluxes, which have non-local characteristics. Instantaneous alpha particle transport fluxes are significant (Γ α ~ 2 x 10 20 m –2 s –1 ) for the reversed shear regime and lower (Γ α ~ 1.1 x 10 20 m –2 s –1 ) for the monotonic q-profile case.

AE↗

DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: III. Novel sub-cyclotron modes

Dedicated DIII-D experiments studied instabilities above and below the ion cyclotron frequency by injecting thirteen different neutral beam populations into L-mode plasmas with different thermal compositions of H, D, and 3 He at five different values of toroidal field. Two companion papers describe the results of these experiments in detail. This paper focuses on data and analysis of two instabilities that occur at ${\sim}0.6$ of the deuterium and hydrogen cyclotron frequencies $f_\mathrm {cD}$ and $f_\mathrm {cH}$ in mixed hydrogen and deuterium plasmas. The modes below $f_\mathrm {cD}$ are often driven unstable by deuterium beams but are rarely destabilized by hydrogen beams. The modes below $f_\mathrm {CD}$ are reported here for the first time in a fusion plasma. These $\mathrm{sub-}f_\mathrm {cH}$ modes are only driven by deuterium beams but, unlike the $\mathrm{sub-}f_\mathrm {cD}$ modes, they are only unstable in plasmas with deuterium concentrations that are ${\lesssim}20\%$; they also have a smaller amplitude than the sub-deuterium modes. Addition of fully-ionized 3 He to the plasma has little effect on the sub-cyclotron deuterium modes but is stabilizing for the sub-cyclotron hydrogen modes. The modes are excited in the core and have toroidal mode numbers around 20. These trends are compared with calculations of the growth rate in a simplified model. The theory successfully predicts global Alfvén eigenmode instability with frequencies and mode numbers close to experimental values, the stronger amplitude of the sub-deuterium modes, and the species-mix dependence of the sub-hydrogen modes. Increased damping associated with the ion-ion hybrid, cutoff, and crossover frequencies plausibly accounts for the stabilization of the $\mathrm{sub-}f_\mathrm {cH}$ modes when 3 He is puffed. However, the relative stability of $\mathrm{sub-}f_\mathrm {cH}$ modes during hydrogen injection remains unexplained.

electromagnetic ion cyclotron instability↗

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 OSTI2026

Validation of Hermes-3 turbulence simulations against the TCV-X21 diverted L-mode reference case

Electrostatic flux-driven turbulence simulations with the Hermes-3 code are performed in TCV L-mode conditions in forward and reversed toroidal field configurations, and compared to the TCV-X21 reference dataset (Oliveira et al 2022 Nucl. Fusion 62 096001) qualitatively and with a quantitative methodology. Using only the magnetic equilibrium, total power across the separatrix (120 kW) and total particle flux to the targets (3 x 10 21 s −1 ) as inputs, the simulations produce time-averaged plasma profiles in good agreement with experiment. Shifts in the target peak location when the toroidal field direction is reversed are reproduced in simulation, including the experimentally observed splitting of the outer strike point into two density peaks. The overall normalized discrepancy between simulation and observation is better than any previously reported in the reversed field configuration, and matches the best previously reported in forward field configuration. Differences between simulation and experiment include density profiles inside the separatrix and at the inner target in forward (favorable $\bigtriangledown B$) field configuration. These differences in target temperature in forward field configuration lead to differences in the balance of current to the inner and outer divertor in the private flux region. The cause of these differences is most likely the lack of neutral gas in these simulations, indicating that even in low recycling regimes neutral gas plays an important role in determining edge plasma profiles. These conclusions are consistent with findings in Oliveira et al (2022 https://github.com/SPCData/TCV-X21).

Physics - Plasma physics↗

DOE OSTI2026

Macroscopic trends of neoclassical tearing stability in high-field H-mode tokamak pilot plants

The neoclassical tearing mode (NTM) stability metric—minimum marginally stable island width $w$$^{*}_{m}$—was compared across 14651 inductive high-field tokamak pilot plant equilibria. Larger devices with reduced elongation and/or increased minor radius demonstrated an order-of-magnitude increase in $w$$^{*}_{m}$, primarily due to a reduction in bootstrap drive. This work is part of an ongoing effort to ensure passive NTM-stability in the ARC tokamak, in which the technology to achieve active tearing-suppression with localised electron cyclotron current drive does not yet exist. The equilibrium scenarios in the database were Monte Carlo generated and normalised to the same >400MW fusion power, minimum pressure scenario at a range of plasma currents, before tearing analysis using the modified Rutherford equation was applied for all resonant poloidal and toroidal m, n modes up to n = 4. Single-helicity toroidal Δ' calculations in resistive DCON set the minimum marginally stable island width, and a simple modal scaling proportional to –m 2 n –1 was identified for high-m Δ' values. The dominant correlates of $w$$^{*}_{m}$ and Δ' across the database were analysed using interpretable machine learning techniques.

NTM seeding↗