Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space plasma physics”

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

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗
DOE OSTI2026

Multi‐Species Energy‐Banded Ions in the Ionosphere During the 21 January 2005 Magnetic Storm: Low‐Altitude Edge of the Warm Plasma Cloak

During the 21–22 January 2005 magnetic storm, the FAST satellite observed warm (< few keV) ions in discrete energy bands on the dayside at ∼3,000 km altitude for more than 6.5 hr. We suggest that the ionospheric energy-banded ions represent the low-altitude edge of the warm plasma cloak observed simultaneously by magnetospheric satellites. This is a clear example of the multi-species ion energy bands (10 eV to several keV) observed during strong magnetic storms by the FAST satellite, stretching from the diffuse auroral region to the plasmapause with lifetimes up to 12 hr. The close association of these energy-banded ions with magnetic storms, their broad latitudinal extent, and the presence of multiple ion species in the same energy band, rather than at the same velocity, indicate that this is a distinct phenomenon from other types of energy-banded ions. During the 21–22 January 2005 magnetic storm, the dayside ion energy band structures, centered at 10 eV (H + ), 40 eV (H + and He + ), and 160 eV (H + , He + , and O + ), were consistent with a “time-of-flight and velocity filter” formation process acting on a near-cusp, impulsive outflow of a <200 eV multi-species ion-source population, poleward and in the same hemisphere as FAST. Understanding the sources and dynamics of warm energy-banded ions and their linkage to the warm plasma cloak is important because during superstorms these ions are transported to L values as low as L ∼ 1.2 in the dawn sector, significantly altering the energetics of the mid-latitude ionosphere.

58 GEOSCIENCES↗
DOE OSTI2025

Zero curvature is a necessary and sufficient condition for a spin-orbital decomposition

There has been an extended debate regarding the existence of a spin-orbital decomposition of the angular momentum of photons and other massless particles. It was recently shown that there are both geometric and topological obstructions preventing any such decomposition. Here we show that any geometric connection on a particle’s state space induces a splitting of the angular momentum into two operators. These operators are well-defined angular momentum operators if and only if the connection has zero curvature. Massive particles have two canonical curved connections corresponding to boosts and rotations, respectively. Furthermore, these can be uniquely combined to produce a flat connection, and this gives a novel derivation of the Newton-Wigner position operator and the corresponding spin and orbital angular momenta for relativistic massive particles. When the mass is taken to zero, transverse boosts and rotations degenerate, leaving only a single connection for massless particles. This connection produces a commonly proposed splitting of the massless angular momentum into two operators. However, the connection is not flat, explaining why these operators do not satisfy the angular momentum commutation relations and are thus not true spin and orbital angular momentum operators.

Angular momentum↗
NASA NTRSDate not supplied

Next Generation Big Data Storage for Long Space Missions

This paper presents the results of the HELIOS (Hardened Extremely Long Life In-formation Optical Storage) mission on the International Space Station (ISS) which tested a unique solution for the long-term storage and retrieval of data in space. For this mission Creative Technology (CTech) developed test media—termed WORF (Write Once, Read Forever)—to validate whether this patented technology will survive all critical parameters for harsh space-based environments including microgravity and ionizing radiation. The HELIOS experiment confirmed that the WORF media is impervious to ionizing radiation, microgravity, solar (plasma) eruptions, and the stress from 8 Gs of the launch including extreme temperature expo-sure. The principal results indicate that there has been no discernible degradation of the media after 8 months on the ISS as compared to a control set of media stored on the ground. This data validated the media’s survivability for harsh space environments for long-term and deep space missions. In addition to the space environment, we are confident that WORF technology can be used for data storage where space-related and other long-term or archival integrity is critical such as: geospatial collections from satellites; space weather archives; past, ongoing, and future space mission media and documentation files; the deep space Gate-way program; as well as Big Data applications such as the Vera C. Rubin astronomical observatory (formerly the LSST). WORF technology for the HELIOS experiment uses a proven archival media, redesigned, re-purposed and patented by CTech to store digital data for long periods, measured in decades and possibly centuries. The media stores standing waves embedded in a substrate that capture the precise col-ors or wavelengths projected onto the media. The colors represent numerical data, with each data location storing multiple superimposed wavelengths, which facilitate the storage of multiple data bytes (rather than just zeros and ones); advanced mathematical permutations allow for extremely large data density equal to or greater than contemporary data storage de-vices. These colors cannot fade or degrade over time since the standing waves are physically stabilized (fully oxidized) metallic silver; no dyes are embedded for this storage system, and silver ions resist micro-bacterial and fungal contamination

Rodney Grubbs↗
NASA NTRS2019

Highly Structured Slow Solar Wind Emerging From an Equatorial Coronal Hole

At solar minimum, the solar wind is observed at high solar latitudes as a predominantly fast (> 500 km/s), highly Alfvenic, rarefied stream of plasma originating deep within coronal holes, while near the ecliptic plane it is interspersed with a more variable slow (< 500 kms) wind. The precise origins of the slow wind streams are less certain, with theories and observations supporting sources from the tips of helmet streamers, interchange reconnection near coronal hole boundaries, and origins within coronal holes with highly diverging magnetic fields. The heating mechanism required to drive the solar wind is also an open question and candidate mechanisms include Alfven wave turbulence, heating by reconnection in nanoflares, ion cyclotron wave heating and acceleration by thermal gradients1. At 1 au, the wind is mixed and evolved and much of the diagnostic structure of these sources and processes has been lost. Here we present new measurements from Parker Solar Probe at 36 to 54 solar radii that show clear evidence of slow, Alfvenic solar wind emerging from a small equatorial coronal hole. The measured magnetic field exhibits patches of large, intermittent reversals associated with jets of plasma and enhanced Poynting flux and interspersed in a smoother and less turbulent flow with near-radial magnetic field. Furthermore, plasma wave measurements suggest electron and ion velocity-space micro-instabilities that have been identified with plasma heating and thermalization processes. Our measurements suggest an impulsive mechanism associated with solar wind energization and a heating role for micro-instabilities and provide strong evidence for low latitude coronal holes as a significant contribution to the source of the slow solar wind.

Solar physics↗
DOE OSTI2025

Modeling fast ion losses due to tearing and internal kink perturbations in MAST-U

Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of fishbone instability through zonal flows driven by energetic particle transport in tokamak plasmas

Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the&#xD;fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Elevated stress hormone levels relate to Epstein-Barr virus reactivation in astronauts

OBJECTIVE: The objective of this study was to determine the effects of stress and spaceflight on levels of neuroendocrine hormones and Epstein-Barr virus (EBV)-specific antibodies in astronauts. METHODS: Antiviral antibody titers and stress hormones were measured in plasma samples collected from 28 astronauts at their annual medical exam (baseline), 10 days before launch (L-10), landing day (R+0), and 3 days after landing (R+3). Urinary stress hormones were also measured at L-10 and R+0. RESULTS: Significant increases (p <.01) in EBV virus capsid antigen antibodies were found at all three time points (L-10, R+0, and R+3) as compared with baseline samples. Anti-EBV nuclear antigen antibodies were significantly decreased at L-10 (p <.05) and continued to decrease after spaceflight (R+0 and R+3, p <.01). No changes were found in antibodies to the nonlatent measles virus. The 11 astronauts who showed evidence of EBV reactivation had significant increases in urinary epinephrine and norepinephrine as compared with astronauts without EBV reactivation. CONCLUSION: These findings indicate that physical and psychological stresses associated with spaceflight resulted in decreased virus-specific T-cell immunity and reactivation of EBV.

NASA Center JSC↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2025

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First fluctuation measurements using an Imaging Neutral Particle Analyzer on DIII-D

A recent upgrade to the Imaging Neutral Particle Analyzer (INPA) on DIII-D has allowed for the first fluctuation measurements using an INPA to be taken during a neoclassical tearing mode (NTM). The INPA signal tracked the mode over 150 ms as the mode frequency dropped to zero, capturing both the NTM with poloidal and toroidal mode numbers $m/n=2/1$ and a $3/2$ mode at double the frequency. Analysis shows that the signals originate from charge exchange events near the edge of the plasma, and relative fluctuation amplitudes are greater than 25% for the duration of the NTM. Filtered signals show frequency beating patterns that are phase-space dependent. Simulated signal is dominated by prompt transport from the neutral beams to the INPA sightline, while the contribution from the slowing down distribution is significantly lower. Simulated measurements in the range of pitches that the diagnostic is sensitive to (0.5≤|v∥/v|≤ 0.75) show the signal is dominated by trapped orbits that pass through magnetic islands near the edge of the plasma. Calculations of expected fluctuation levels show that only direct interaction with the NTM can provide the strong relative fluctuation levels seen experimentally. The prompt nature of the orbits and thin radial layer found to contribute to synthetic signals suggest INPA passive data may be used to measure the perturbations of confined orbits on a single pass through a plasma instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Low Loss, Fast Ramp HTS Solenoid Prototype for Compact Spherical Tokamaks

To explore AC loss characteristics of commercial Bi-2212 conductors made in the USA, a small ohmic heating (OH) prototype coil is being developed for compact tokamaks using high performance Bi-2212 wires newly manufactured by Bruker EST. Here, the goal is to build and test Bi-2212 prototype coils for fast ramp central solenoid (CS) operation of spherical tokamaks (ST). For a ST pilot plant, one of the key challenges is the tight space available for CS in a center column where in-board legs of the toroidal field (TF) and OH coils are assembled into a cryostat with a few millimeters clearance. A high current density Rutherford cable consisting of 17 legacy Bi-2212 wires, a contribution of the US Magnet Development Program (MDP), was fabricated at the Lawrence Berkeley National Laboratory (LBNL) for the development of fast ramp CS magnets in the ST Advanced Reactor (STAR) designed by PPPL. Characterization of the new and legacy HTS wires is presented, along with the coil design wound directly with new wire, and a Rutherford cable-wound solenoid using the legacy wire. A comparison of AC losses is made based on expected performance for a fast ramp ST operation. The HTS cable and cabled coil fabrication issues such as the heat treatment, coil mandrel design, assembly and 2212 quench characteristics are discussed, along with a test plan for CS design validation.

Zhai, Yuhu [Princeton Plasma Physics Laboratory (P↗
DOE OSTI2025

Turbulent Heating in Collisionless Low-beta Plasmas: Imbalance, Landau Damping, and Electron–Ion Energy Partition

An understanding of how turbulent energy is partitioned between ions and electrons in weakly collisional plasmas is crucial for modeling many astrophysical systems. Using theory and simulations of a four-dimensional reduced model of low-beta gyrokinetics (the “Kinetic Reduced Electron Heating Model”), we investigate the dependence of collisionless heating processes on plasma beta and imbalance (normalized cross-helicity). These parameters are important because they control the helicity barrier, the formation of which divides the parameter space into two distinct regimes with remarkably different properties. In the first, at lower beta and/or imbalance, the absence of a helicity barrier allows the cascade of injected power to proceed to small (perpendicular) scales, but its slow cascade rate makes it susceptible to significant electron Landau damping, in some cases leading to a marked steepening of the magnetic spectra on scales above the ion Larmor radius. In the second, at higher beta and/or imbalance, the helicity barrier halts the cascade, confining electron Landau damping to scales above the steep “transition-range” spectral break, resulting in dominant ion heating. We formulate quantitative models of these processes that compare well to simulations in each regime, and combine them with results of previous studies to construct a simple formula for the electron–ion heating ratio as a function of beta and imbalance. This model predicts a “winner takes all” picture of low-beta plasma heating, where a small change in the fluctuations' properties at large scales (the imbalance) can cause a sudden switch between electron and ion heating.

Interplanetary turbulence↗
NASA NTRS2011

Fast Plasma Investigation for MMS: Simulation of the Burst Triggering System

The Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. To meet these requirements, the Fast Plasma Instrument (FPI) consists of eight (8) identical half top-hat electron sensors and eight (8) identical ion sensors and an Instrument Data Processing Unit (IDPU). The sensors (electron or ion) are grouped into pairs whose 6 degree x 180 degree fields-of-view (FOV) are set 90 degrees apart. Each sensor is equipped with electrostatic aperture steering to allow the sensor to scan a 45 degree x 180 degree fan about the its nominal viewing (0 deflection) direction. Each pair of sensors, known as the Dual Electron Spectrometer (DES) and the Dual Ion Spectrometer (DIS), occupies a quadrant on the MMS spacecraft and the combination of the eight electron/ion sensors, employing aperture steering, image the full-sky every 30-ms (electrons) and 150-ms (ions), respectively. To probe the diffusion regions of reconnection, the highest temporal/spatial resolution mode of FPI results in the DES complement of a given spacecraft generating 6.5-Mb (raised dot) per second of electron data while the DIS generates 1.1-Mb (raised dot) per second of ion data yielding an FPI total data rate of 6.6-Mb (raised dot) per second. The FPI electron/ion data is collected by the IDPU then transmitted to the Central Data Instrument Processor (CIDP) on the spacecraft for science interest ranking. Only data sequences that contain the greatest amount of temporal/spatial structure will be intelligently down-linked by the spacecraft. This requires a data ranking process known as the burst trigger system. The burst trigger system uses pseudo physical quantities to approximate the local plasma environments. As each pseudo quantity will have a different value, a set of two scaling factors is employed for each pseudo term. These pseudo quantities are then combined at the instrument, spacecraft, and observatory level leading to a final ranking of data based on expected scientific interest. Here, we present simulations of the fixed point burst trigger system for the FPI. A variety of data sets based on previous mission data as well as analytical formulations are tested. Comparisons of floating point calculations versus the fixed point hardware simulation are shown. Analysis of the potential sources of error from overflows, quantization, etc. are examined and mitigation methods are presented. Finally a series of calibration curves are presented, showing the expected error in pseudo quantities based solely on the scale parameters chosen and the expected data range. We conclude with a presentation of the current base-lined FPI burst trigger approach.

Barrie, A. C.↗
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↗