Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal magnetic fields”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Islands and current singularities in quasisymmetric toroidal plasmas

The presence of current singularities in a quasisymmetric magnetic field is explored. Quasisymmetry is shown effective in isolating Pfirsch–Schlüter singularities, to leading order, to a single magnetic surface resonant with the helicity of the symmetry. The effects of pressure driven currents are analyzed, indicating that exclusion of this surface from the plasma volume reduces the potential opening of islands but does not generally eliminate them completely due to higher order asymmetric geometric effects. These three-dimensional effects are contained in quasisymmetry and indicate the complexity of finding consistent solutions and their potential sensitivity. The δ-function current singularities show a distinct quasisymmetric behavior only when the higher-order Fourier content of B is relevant for the resistive stability parameter D R (not included in leading-order near-axis expansions). In such scenarios, quasisymmetry proves advantageous, both in simplicity and avoidance of amplification by low-order rational surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high- β are reduced to the level of that of collisionless trapped electron modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dependence of alpha-particle-driven Alfven eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conven- tional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfven eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently- developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvenic or nearly sub-Alfvenic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Neutronics analysis of spin-polarized fuel in spherical tokamaks

The effect of spin-polarized deuterium–tritium fuel on the blanket, magnets, and surrounding structures is analyzed for a spherical tokamak fusion power plant. The focus of this paper is the neutron flux distribution in the fusion device and the corresponding engineering considerations, such as tritium breeding ratio (TBR) and magnet fluence. Significant benefits were found for anti-aligned polarization: a 2.7% increase in the total TBR (which could increase significantly with further optimization) and an approximately 68% increase in the magnet lifetime. If the inboard breeding blanket can be replaced with extra shielding and magnets, then it could significantly increase fusion power with no loss of magnet lifetime. This would be accomplished by increasing the toroidal field strength.

fusion↗
DOE OSTI2023

New linear stability parameter to describe low-β electromagnetic microinstabilities driven by passing electrons in axisymmetric toroidal geometry

In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Verification of the Shafranov shift in free-boundary VMEC, DESC, and SPEC calculations with elliptical geometry

Here, given a set of current-carrying filaments that create a magnetic field with rotating elliptical flux surfaces, we present verification calculations using the VMEC, DESC, and SPEC free-boundary magnetohydrostatic equilibrium codes for both the “vacuum” equilibrium, for which there is no plasma pressure and no plasma currents, and for a non-zero pressure, zero toroidal current equilibrium. For the vacuum case, these codes are quantitatively compared to the magnetic field produced by the coils. For the non-zero pressure case, they are compared to each other. As the stepped pressure profile used in SPEC approaches the continuous profile used in VMEC and DESC, the SPEC magnetic axis approaches the VMEC and DESC axes. To within an error given by the difference between the pressure profiles, the VMEC, DESC, and SPEC calculations give the same pressure-induced geometric shift of the magnetic axis, known as the Shafranov shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High-Field Ignition-Capable Stellarator i-ASTER: Initial Structural Evaluation

i-ASTER is a preliminary conceptual design of a high-field pulsed experimental stellarator aimed essentially at achieving and understanding plasma ignition. The averaged magnetic field on the axis of i-ASTER.v1 is B ~ 10 T and the plasma volume is 30 m 3 . A toroidal monolithic support structure, which is located externally to the modular coils, is defined to accurately arrange the coils and withstand the magnetic Lorentz forces from the high magnetic fields. The magnetic field in the coils and the resulting Lorentz forces on the coils were calculated and utilized for the stress calculation in monolithic support of constant thickness. Subsequently, a monolithic support of variable thickness is proposed and studied. The thickness of the structure is variable and adjusted to the magnitude of the mechanical stress in each area of the toroidal structure so as to smooth the stresses while minimizing the quantity of material used. A fabrication method for the monolithic structure has been tested. For that, a scaled-down mock-up of one period of stellarator has been additively manufactured and strengthened with fiber-reinforced resin. Unsurmountable difficulties have not been found for this fabrication method.

Queral, Vicente↗
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 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 OSTI2020

Geometric concepts for stellarator permanent magnet arrays

The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.

stellarator, permanent magnet, optimization↗
DOE OSTI2020

Geometric concepts for stellarator permanent magnet arrays

The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Here, two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Magnetic geometry effects on turbulent density fluctuations in Wendelstein 7-X

Line-integrated density fluctuation measurements from the core of the optimised stellarator Wendelstein 7-X (W7-X) are studied in the context of magnetic field geometry changes. The amplitude as well as the spectral distribution of fluctuations is affected by varying the magnetic field configuration. Weaker turbulent fluctuations are observed in configurations with higher rotational transform and lower magnetic mirror depth. Simulations with the gyrokinetic codes stella and GENE-3D are used to investigate the role of magnetic geometry for the experimental observations. The stabilising effect of the high rational transform configuration is reproduced by simulations and can be associated with stronger shaping of the magnetic field. This leads to increased finite Larmor radius and local shear effects, which limit the toroidal extent of ion temperature gradient (ITG) modes. For W7-X configurations with higher magnetic mirror depth, gyrokinetic theory and simulations predict weaker fluctuations, since they are closer to having the maximum-J property, which has a stabilising effect on trapped-electron mode and ITG. The opposite trend is observed experimentally, which can therefore not be explained by geometry effects on turbulence alone.

Wendelstein-7X↗
DOE OSTI2021

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations

The structure of the non-axisymmetric heat load distribution at the divertor plates is determined not only by the toroidal but also from the poloidal spectrum of non-axisymmetric field perturbations. Whether they are intrinsic, like error fields, or they are applied through 3D coils, the non-axisymmetric fields produce complex 3D edge magnetic topologies (footprints) that alter the properties of the heat and particle flux distributions on the divertor target plates. In this manuscript, a study of the impact of applied 3D field poloidal spectrum on the footprint size and structure is done for the DIII-D tokamak using the resistive MHD code M3D-C1 coupled with the field line tracing code TRIP3D. To resolve the impact of the poloidal spectrum of the magnetic perturbation, the relative phase of the two rows of in-vessel 3D coils used to produce both a n = 2 and a n = 3 perturbation is varied, where n is the toroidal harmonic of the magnetic perturbation. This shows that the largest footprint is predicted when the relative phase of the two rows is close to zero, which is also where the resonant coupling with the plasma is maximized. These results suggest that it will be challenging to decouple the footprint size from the requisite resonant coupling for RMP–ELM control. The correlation between the measured heat load and particle flux distributions at the outer divertor plates in DIII-D and the magnetic measurements is in good agreement with the predicted dependence of the magnetic footprint size on the amplitude of the resonant component of the plasma response.

3D fields↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗