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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Estimating coil features from an equilibrium

We present an explicit theoretical framework for constructing artificial modular coils for vacuum stellarator fields based solely on equilibrium properties, achieved through the formulation of a current potential defined on flux surfaces. Contours of constant Boozer toroidal angle can be directly interpreted as proxy coils, and so we demonstrate that key measures of coil complexity—particularly coil non-planarity—are strongly governed by local magnetic-field properties. This approach shows promise as predictor for more realistic coil configurations, providing both a pathway toward deeper understanding of equilibrium–coil relationships and a potential practical proxy for coil design.

Rodríguez, Eduardo [Max Planck Institute for Plasm↗
DOE OSTI2024

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Poor confinement in stellarators at high energy

High-energy particle resonances can modify particle distributions and even cause significant particle loss. Resonances can be present in any toroidal confinement device and can easily be found numerically. Many stellarators have weak magnetic shear so that large islands and large chaotic regions can be produced by resonant perturbations with small amplitudes. Additionally, while the choice of the field line helicity profile in the plasma can limit the presence of resonances at low particle energy, the resonance location is energy-dependent, and they can move into the plasma at higher energy. If resonances match the toroidal variation of the equilibrium, they can produce wide islands in the phase space of orbits even in the absence of perturbations due to instabilities. These islands increase in size with particle energy and can seriously affect the confinement of high-energy ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗
DOE OSTI2025

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗
NASA NTRS1977

Inward transport of a toroidally confined plasma subject to strong radial electric fields

The paper aims at showing that the density and confinement time of a toroidal plasma can be enhanced by radial electric fields far stronger than the ambipolar values, and that, if such electric fields point into the plasma, radially inward transport can result. The investigation deals with low-frequency fluctuation-induced transport using digitally implemented spectral analysis techniques and with the role of strong applied radial electric fields and weak vertical magnetic fields on plasma density and particle confinement times in a Bumpy Torus geometry. Results indicate that application of sufficiently strong radially inward electric fields results in radially inward fluctuation-induced transport into the toroidal electrostatic potential well; this inward transport gives rise to higher average electron densities and longer particle confinement times in the toroidal plasma.

Roth, J. R.↗
DOE OSTI2024

Diocotron modes in pure electron plasmas in the APEX levitating dipole trap

Abstract Observations of persistent toroidal modes in dipole-trapped pure electron plasmas are presented. Non-neutral electron plasmas were confined by the poloidal magnetic field of the APEX (A Positron Electron eXperiment) levitating dipole trap. The negative bias of the electron emitter was gated to control the filling time. Large amplitude, narrow-band oscillations were consistently observed in wall probe measurements after the fill. The oscillatory image-charge currents induced by the trapped electron plasma typically persisted for many seconds with fundamental frequencies, f 0 = 50 to 500 kHz, that scale with the E → × B → drift velocity. We attribute this feature to the dipole-equivalent of the diocotron mode that is commonly observed in non-neutral plasmas in linear traps.

Deller, A. (ORCID:0000000234301501)↗
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of MARS-F modeling of plasma response to RMPs using internal measurements on DIII-D

The linear resistive plasma response model is validated against the plasma internal measurement data from DIII-D edge-localized mode (ELM) control experiments with applied resonant magnetic perturbation (RMP). Considered are DIII-D discharges where the n = 1, 2, and 3 (n is the toroidal mode number) RMP field was applied. Experimental data for the plasma boundary displacement, as well as the three-dimensional (3D) pressure perturbation in the edge pedestal region, are deduced from the vertical Thomson scattering (TS) system and the horizontal charge exchange recombination (CER) system on DIII-D. The linear response model produces results that are in reasonable quantitative agreement with the DIII-D internal measurements. The plasma boundary displacement of up to 15 mm is modeled, with the pedestal pressure perturbation reaching 3 kPa. As an important insight, the larger plasma displacement measured by the vertical TS system, as compared to that measured by the horizontal CER system, is due to the contribution from the tangential component of the plasma displacement to the former. This mixing of displacement components is also found to influence the sensitivity of the CER measurement comparisons. The results of this study provide further confidence in the linear resistive plasma response model for analyzing ELM control experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Characteristics of the NASA Lewis bumpy-torus plasma generated with positive applied potentials

Experimental observations were made during steady-state operation of a bumpy-torus plasma at input powers up to 150 kW in deuterium and helium gas and with positive potentials applied to the midplane electrodes. In this steady-state ion heating method a modified Penning discharge is operated such that the plasma is acted upon by a combination of strong electric and magnetic fields. Experimental investigation of a deuterium plasma revealed electron temperatures from 14 to 140 eV and ion kinetic temperatures from 160 to 1785 eV. At least two distinct modes of operation exist. Experimental data shows that the average ion residence time in the plasma is virtually independent of the magnetic field strength. Data was taken when all 12 anode rings were at high voltage, and in other symmetric configurations in which the toroidal plasma was generated by applying positive potentials to six anode rings, three anode rings, and a single anode ring.

Roth, J. R.↗
DOE OSTI2023

A divergence-free constrained magnetic field interpolation method for scattered data

An interpolation method to evaluate magnetic fields, given its unstructured and scattered magnetic data, is presented. The method is based on the reconstruction of the global magnetic field using a superposition of orthogonal functions. The coefficients of the expansion are obtained by minimizing a cost function defined as the L2 norm of the difference between the ground truth and the reconstructed magnetic field evaluated on the training data. The divergence-free condition is incorporated as a constraint in the cost function, allowing the method to achieve arbitrarily small errors in the magnetic field divergence. An exponential decay of the approximation error is observed and compared with the less favorable algebraic decay of local splines. Compared to local methods involving computationally expensive search algorithms, the proposed method exhibits a significant reduction of the computational complexity of the field evaluation, while maintaining a small error in the divergence even in the presence of magnetic islands and stochasticity. Finally, applications to the computation of Poincaré sections using data obtained from numerical solutions of the magnetohydrodynamic equations in toroidal geometry are presented and compared with local methods currently in use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

HTS to NbTi Joints for Fusion Magnet Cost Savings

This paper describes methods of employing high field High Temperature Superconductors (HTS), only where high field capability and high current density is needed in superconducting coils. The approach lends itself, in particular, to Toroidal Field (TF) coils used in fusion reactors. The initial idea was to develop a conductor joint design or designs that would allow joints between expensive high field High Temperature Superconductors (HTS) Rare Earth Barium Copper Oxide (REBCO) conductors to low cost NbTi lower field conductors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recent advance progress of HL-3 experiments

Since the first plasma realized in 2020, a series of key systems on HL-3 (known as HL-2M before) tokamak have been equipped/upgraded, including in-vessel components (the first wall, lower divertor, and toroidal cryogenic/water-cooling/baking/glow discharge systems, etc.), auxiliary heating system of 11 MW, and 28 diagnostic systems (to measure the plasma density, electron temperature, radiation, magnetic field, etc.). Magnet field systems were commissioned firstly for divertor plasma discharges. During the 2nd experimental campaign of HL-3 tokamak, several great progresses have been achieved. Firstly, the successful operation with plasma current larger than 1 MA was achieved under a divertor configuration. Secondly, the advanced divertor concept with two distinct snowflake configurations was realized. It is found that the distribution of ion saturation current and heat flux on bottom plate becomes wide due to magnetic surface expansion, demonstrating the advantage of such configuration in the heat flux mitigation. In addition, using the combination of NBI, ECRH and LHCD, the standard sawtoothing high confinement mode of megampere plasma was firstly accessed on the HL-3. The successful commissioning of HL-3 is beneficial for the initial operation of ITER.

1 MA plasma↗
NASA NTRS2019

Reassessment of Electromagnetic Core-Mantle Coupling and Its Implications to the Earth's Decadal Polar Motion

The observed Earth's polar motion on decadal time scales has long been conjectured to be excited by the exchange of equatorial angular momentum between the solid mantle and the fluid outer core, via the mechanism of electromagnetic (EM) core-mantle coupling. However, past estimations of the EM coupling torque from surface geomagnetic observations is too weak to account for the observed decadal polar motion. Our recent estimations from numerical geodynamo simulations have shown the opposite. In this paper, we re-examine in detail the EM coupling mechanism and the properties of the magnetic field in the electrically conducting lower mantle (characterized by a thin D″-layer at the base of the mantle). Our simulations find that the toroidal field in the D″-layer from the induction and convection of the toroidal field in the outer core could be potentially much stronger than that from the advection of the poloidal field in the outer core. The former, however, cannot be inferred from geomagnetic observations at the Earth's surface, and is missing in previous EM torque estimated from geomagnetic observations. Our deduction suggests further that this field could make the actual EM coupling torque sufficiently strong, at approximately 5 × 10(exp 19) Nm, to excite, and hence explain, the decadal polar motion to magnitude of approximately 10 mas.

Kuang, Weijia↗