Search NASA⌕ Search

SEARCH · Search NASA

Results for “toroidal magnetic confinement”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Magnetic fields with general omnigenity

Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has been either axisymmetric, quasi-symmetric or has poloidally closed contours of magnetic field strength $B$ . However, general omnigenous equilibria are a much larger design space than these subsets. A new model is presented and employed in the DESC stellarator optimization suite to represent and discover the full parameter space of omnigenous equilibria. Although exact omnigenity aside from quasi-symmetry is impossible, these results reveal that excellent particle confinement can be achieved in practice. Examples far from quasi-symmetry with poloidally, helically and toroidally closed $B$ contours are attained with DESC and shown to have low neoclassical collisional transport and fast particle losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Lyapunov-based nonlinear control of nonautonomous systems with individual input constraints

A control algorithm that can locally stabilize a specific class of multi-input multi-output nonautonomous nonlinear dynamical systems while satisfying individual input constraints is developed. The proposed Lyapunov-based state-feedback control law inherently accounts for the actuator amplitude saturation limits without the need for computationally expensive real-time optimization techniques. In addition to the control law, a formal definition for the local “controllable region” within which the controller can asymptotically drive the system states to the origin and satisfy the input saturation limits is also presented. The nonautonomous nature of the system dynamics implies that the “controllable region” continuously evolves with time. Therefore, a sufficient condition to maintain the system states within the “controllable region” is proposed in this work to make practical implementation feasible. The effectiveness of the controller is tested for a specific control problem arising in tokamaks, which are toroidal devices that use strong magnetic fields to confine a plasma (hot ionized gas). Here, the primary emphasis of tokamak research is to regulate the plasma properties around predetermined values to achieve stable plasma confinement. Nonlinear simulations show that the proposed controller can achieve the desired plasma control objectives in a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle resonances in toroidal fusion devices

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through avalanche. The existence of a mode particle resonance depends on properties of the equilibrium, particle trajectories, and perturbation mode harmonic content. Resonance location is a function of particle energy and equilibrium field line helicity. Different methods for finding resonance location and energy dependence are developed. The properties of resonances are discussed using examples from magnetic fusion devices. Here, we show that if mode resonances exist at low particle energy, they very likely also exist at high energy, thus modifying high energy beam particles and fusion products. It is possible for a resonance to appear due to mode induced orbit helicity modification when it is forbidden at low mode amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Avalanche transport of energetic-ions in magnetic confinement plasmas: nonlinear multiple wave-number simulation

Large burst activity, identified as toroidal Alfvén eigenmode (TAE) avalanche, occurs frequently in neutral-beam heated plasmas in National Spherical Torus Experiment (NSTX). Based on the typical experimental observation of TAE avalanche on NSTX, a self-consistent nonlinear multiple wave-number (k ∥ ≃ n/R, where n toroidal mode-number and R major radius) simulation associated with TAE avalanches is performed using the experimental parameters and profiles before the occurrence of TAE avalanche as the M3D-K input. The wave–wave nonlinear coupling among different modes and the resonant interaction between different modes and energetic-ions during TAE avalanches are identified in the nonlinear multiple wave-number simulations. The resonance overlap during the TAE avalanche is clearly observed in the simulation. It is found that the effective wave–wave coupling and a sufficiently strong drive are two important ingredients for the onset of TAE avalanches. TAE avalanche is considered to be a strongly nonlinear process and it is always accompanied by the simultaneous rapid frequency-chirping and large amplitude bursting of multiple modes and significant energetic-ion losses. The experimental phenomenon is observed on NSTX and is qualitatively reproduced by the simulation results in this work. These findings indicate that the onset of avalanche is triggered by nonlinearity of the system, and are also conducive to understanding the underlying mechanism of avalanche transport of energetic particles in the future burning plasmas, such as International Thermonuclear Experiment Reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fast neural Poincaré maps for toroidal magnetic fields

Poincaré maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincaré map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We propose a novel neural network architecture, the HénonNet, and show that it is capable of accurately learning realistic Poincaré maps from observations of a conventional field-line-following algorithm. After training, such learned Poincaré maps evaluate much faster than the field-line integration method. Moreover, the HénonNet architecture exactly reproduces the primary physics constraint imposed on field-line Poincaré maps: flux preservation. Furthermore, this structure-preserving property is the consequence of each layer in a HénonNet being a symplectic map. We demonstrate empirically that a HénonNet can learn to mock the confinement properties of a large magnetic island by using coiled hyperbolic invariant manifolds to produce a sticky chaotic region at the desired island location. This suggests a novel approach to designing magnetic fields with good confinement properties that may be more flexible than ensuring confinement using KAM tori.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental inference of flux tunneling between magnetic island chains in tokamaks

We report confinement bifurcations due to flux tunneling in magnetic islands in toroidal plasmas. This occurs in 2/1 islands, triggered by decoupling of non-overlapping adjacent 5/2 islands. While subject to electron cyclotron current drive (ECCD), the temperature at the 2/1 O-point is flat before but peaks after decoupling. Here, this is shown to be due to a bifurcation from stochastic to nested magnetic structure, caused by the removal of flux tunneling through intersecting manifolds of overlapping heteroclinic tangles. 2/1 island stability correlates with coupling events in an ITER baseline scenario plasma where ECCD can not stabilize the coupled islands, showing the critical impact of flux tunneling on disruption avoidance in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

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 OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
DOE OSTI2021

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurement of Zero-Frequency Fluctuations Generated by Coupling between Alfvén Modes in the JET Tokamak

We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with theoretical and numerical predictions of zonal modes pumped by Alfvén modes, and is correlated with an increase in the deep core ion temperature, temperature gradient, confinement factor H 89,P , and a reduction in the main ion heat diffusivity. Despite the energetic particle transport induced by the Alfvén eigenmodes, the generation of a zero-frequency fluctuation that can suppress the turbulence leads to an overall improvement of confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Piecewise Omnigenous Stellarators

In omnigenous magnetic fields, charged particles are perfectly confined in the absence of collisions and turbulence. For this reason, the magnetic configuration is optimized to be close to omnigenity in any candidate for a stellarator fusion reactor. However, approaching omnigenity imposes severe constraints on the spatial variation of the magnetic field. In particular, the topology of the contours of constant magnetic field strength on each magnetic surface must be such that there are no particles transitioning between different types of wells. This, in turn, usually leads to complicated plasma shapes and coils. This Letter presents a new family of optimized fields that display tokamak-like collisional energy transport while having transitioning particles. This result radically broadens the space of accessible reactor-relevant configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Periodic Korteweg–de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma- β equilibrium

Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B , that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equations. The nonlinear, overdetermined nature of the QS MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B . Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that B on a flux surface is determined by three or at most four flux functions, each of which determines a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

Differential geometry↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗