Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale analysis”

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 109 records · Page 6

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗

High-Fidelity Analysis of ULF Wave Mode Structure Following Interplanetary Shock Compression of the Dayside Magnetopause Using MMS Multi-Point Observations

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive

This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).

54 ENVIRONMENTAL SCIENCES↗

Gulf of Mexico Turbulence Intensity Model Data

The dataset is a subset of simulations carried out for the Gulf of Mexico region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing. Simulations are performed from January through June in 2020 and output every 10 minutes, which permits comparison with corresponding lidar observations.

17 WIND ENERGY↗

Gulf of Mexico Turbulence Intensity Model Data

The dataset is a subset of simulations carried out for the Gulf of Mexico region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from the NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing. Simulations are performed from January through June in 2020 and output every 10 minutes, which permits comparison with corresponding lidar observations.

17 WIND ENERGY↗

Gulf of Mexico Turbulence Intensity Model Data in Observational Space

The dataset archives model-simulated turbulence intensity and meteorological profiles and timeseries at the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars. The simulated data are interpolated in time and/or space according to observed quantities. The simulations were carried out for the Gulf of Mexico region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing.

17 WIND ENERGY↗

Mid-Atlantic Turbulence Intensity Model Data

The dataset is a subset of simulations carried out for the mid-Atlantic region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Air-Sea Interaction Tower (ASIT) of Woods Hole Oceanographic Institution’s Martha's Vineyard Coastal Observatory are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing. Simulations are performed from February through June in 2020 and output every 10 minutes, which permits comparison with corresponding lidar observations.

17 WIND ENERGY↗

Mid-Atlantic Turbulence Intensity Data in Observational Space

The dataset archives observed and model-simulated turbulence intensity and meteorological profiles and timeseries at the Air-Sea Interaction Tower (ASIT) of Woods Hole Oceanographic Institution’s Martha’s Vineyard Coastal Observatory (MVCO). The observational data were measured by lidar and buoy deployed at ASIT. The simulated profiles and timeseries data are interpolated in time and/or space according to observations. Simulations were carried out for the mid-Atlantic region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles at the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing.

17 WIND ENERGY↗

Northern Pacific Turbulence Intensity Model Data

The dataset is a subset of simulations carried out for the northern Pacific region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the lidar buoys deployed off the coast of California (Humbolt and Morro Bay) are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing. Simulations are performed from November to December in 2020 and May through June in 2021. Model outputs every 10 minutes, permitting comparison with corresponding lidar observations.

17 WIND ENERGY↗

Northern Pacific Turbulence Intensity Model Data in Observational Space

The dataset archives model-simulated turbulence intensity and meteorological profiles and timeseries at the lidar buoys sites off the coast of California (Humboldt and Morro Bay). The simulated data are interpolated in time and/or space according to observed quantities. The simulations were carried out for the north Pacific region using the revised Weather Research and Forecasting (WRF) model version 4.2 that incorporates the implementation of online turbulence intensity (TI) calculations (Tai et al. 2023). The simulated atmospheric profiles near the Shell Exploration and Production Corporation's Tension Leg Platforms Ursa and Mars are archived. Physics parameterizations chosen for the simulations include the Thompson microphysics parameterization, Mellor-Yamada-Nakanishi Niino (MYNN) boundary layer parameterization, Mellor-Yamada-Janjic surface layer parameterization, Unified Noah land-surface parameterization, and the RRTMG longwave and shortwave radiation parameterization. Initial and boundary conditions are taken from NOAA’s High-Resolution Rapid Refresh (HRRR) product. The JPL 0.01-degree Level 4 Multiscale Ultrahigh Resolution (MUR) Global Foundation Sea Surface Temperature (SST) Analysis (V4.1) data are used as the model’s SST forcing.

17 WIND ENERGY↗

Design and Manufacturing of High Performance, Reduced Charge Heat Exchangers

Air-to-fluid heat exchangers (HXs) play a critical role as the main heat transfer component in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) systems. However, their airside thermal resistance significantly inhibits their overall performance. Furthermore, these HXs must be continually more compact to meet the latest refrigerant charge limits to reduce emissions. Recent literature suggests that traditional HX geometries (e.g., round or flat tubes with fins), have reached their limits, and more sophisticated shape- and topology-optimized designs are required to achieve the next jump in performance. This research sheds light on the next generation of air-to-refrigerant HXs and aims to address several practical issues to commercialization such as novelty challenges (improved performance for significant charge reduction; modeling expertise & time investment), manufacturing challenges (non-round tube manufacturing; tube-header integration; product qualification, e.g., burst pressure testing, extreme operational environment, etc.), and operational challenges (flow maldistribution, fouling & wetting, noise & vibration). For example, a >20% improvement on one (or more) HX-level performance metrics (e.g., envelope volume, airside pressure drop, face area, capacity, refrigerant charge, weight, cost, etc.) must be achieved before a HX design is considered for commercialization. We present a new, comprehensive and experimentally validated air-to-refrigerant HX optimization framework with simultaneous thermal-hydraulic performance and mechanical strength considerations for novel, non-round, shape- and topology-optimized tubes capable of optimizing single and two-phase HX designs for any refrigerant choice and performance requirement with significant engineering time savings compared to conventional design practices. The framework was exercised for a wide range of applications and refrigerants, resulting in HXs which achieved greater than 20% improved performance, 20% reductions in size, and 25% reductions in refrigerant charge. To enable non-round tube bundle use in next generation HVAC&R equipment, novel manufacturing techniques were investigated, including the development of conventional manufacturing methods for small diameter, non-round tubes and novel tube-header integration strategies. In total, ten HX prototypes were manufactured, nine using conventional methods directly attributed to this project and one using advanced additive manufacturing methods. The five-year manufacturing feasibility of the proposed HXs was found to have a good outlook. The non-round tube HX simulated performance was validated through comprehensive experimental testing, including nine in-house component-level tests, one independent component-level test at an industry partner laboratory, and in-house system-level tests of using a commercially-available, residential packaged A/C unit which was retrofitted with a non-round tube prototype HX. It was found that HX designs proposed by the new framework can successfully predict experimental thermal-hydraulic performance within ±10-20% the first time with no manual design changes, eliminating the need for time-consuming and expensive prototyping efforts. This work will accelerate design and time to market for next generation HXs while simultaneously facilitating industry transition to new refrigerants at lower charge.

42 ENGINEERING↗

Multiscale Richtmyer-Meshkov instability experiments to isolate the strain rate dependence of strength

Theoretical analysis of Richtmyer-Meshkov instability (RMI) experiments for solid strength shows that the strain rate for a given shock should be inversely proportional to the length scale of the sine wave perturbations when η 0 k , the nondimensional amplitude to wavelength ratio, is held fixed. To isolate the effect of strain rate on strength, free-surface RMI specimens of annealed copper were prepared with three perturbation regions with the same η 0 k but different length scales, characterized by the wavelength λ varying by a factor of 4.9 from 65 to 130 to 320 µ m . Three such targets with different fixed η 0 k ′ s were impacted to a shock pressure of 25 GPa, and the instability evolution was measured with photon Doppler velocimetry. Strengths estimated by comparing hydrocode simulation to the data increased from 700 to 1200 MPa as λ decreased. The different η 0 k targets exercised increasing amounts of plastic strain yet showed no evidence of strain hardening. Physical regime sensitivity analysis determined that for 320 − 65 µ m wavelength perturbations, the effective strain rates increased from 8.7 × 10 6 to 3.3 × 10 7 s − 1 , a factor of 3.8. Thus, the predicted strain rate scaling was mostly achieved but slightly suppressed by increased strength at higher rates. The RMI strength estimates were plotted against constitutive testing data on copper from the literature to show striking evidence of the strength upturn at higher strain rates. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

An integrated coupling model for solving multiscale fluid-fluid coupling problems in SAM code

In this study, an integrated coupling method has been developed for solving multiscale fluid-fluid coupling problems in plant-scale safety analysis models in SAM (System Analysis Module) code. In this method, a higher-fidelity multi-dimensional (3D) flow module is used for reactor components of complex flow features (e.g., reactor core) and a lumped parameter one-dimensional (1D) flow module for plant-scale flow loops (e.g., primary loop pipe network), respectively. In this method, the 3D fluid equation/domain and 1D fluid equation/domain are tightly coupled at the residual level and solved simultaneously using the Newton’s method to overcome the convergence issues typically seen in existing approaches like separate domain approach, where the 3D fluid equation and 1D fluid equation are solved separately. Extensive and successful code verifications and demonstrations have been performed for this newly developed method. This new modeling approach significantly simplify the work flow in developing high-fidelity plant-scale safety analysis model, e.g. for pool-type reactors and pebble-bed reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Using In Situ High-Energy X-ray Diffraction to Quantify Electrode Behavior of Li-Ion Batteries from Extreme Fast Charging

Extreme fast charging (XFC, =15 min charging time) of Li-ion batteries (LIBs) has been proposed as an immediate target to increase the commercial appeal of electric vehicles. However, XFC of LIBs is associated with the degradation of battery performance and safety concerns. Quantitative and simultaneous characterization of various components during cell degradation represents a major experimental challenge. In this work, we outline a methodology for the use of spatially resolved, high-energy X-ray diffraction as a quantitative, in situ method of mapping the degradation of LIBs. We use this approach to study the battery cell capacity loss, both locally (mm scale) and globally over the entire cell (cm scale). Specifically, our workflow allows us to quantify the total amount of plated Li on the anode, as well as its spatial correlation to the structural properties of the anode and cathode. The method complements existing optical methods to resolve the spatial heterogeneity of local degradation mechanisms such as Li plating and provides simultaneous insights into concomitant anode state-of-charge variability. We apply it to commercially relevant single-layer pouch cells with the graphite anode and the LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode. Our results show that Li plating occurs heterogeneously on the graphite anode and that it is spatially correlated to the extent of anode lithiation. In this work, we anticipate that the described workflow will allow for understanding multiscale degradation in energy-storage devices beyond LIBs, where quantitative analysis at a local and global length scale can be performed without the necessity to tear down the device, due to the applicability of high-energy X-rays to probe in situ degradation.

25 ENERGY STORAGE↗

Surface and bulk characterization of reservoir and cap-rocks: Past, present, and future

Geologic formations provide potentially some of the largest volume capacities for CO 2 storage or sequestration. Potential storage sites can be deep saline aquifers, depleted oil reservoirs, and coal seams, surrounded by sealing layers to prevent CO 2 from leaking. It is therefore critical to understand mechanisms contributing to CO 2 trapping and CO 2 leaks. Both phenomena are governed by reactions at the interfaces of the reservoir and cap-rocks and are controlled by the complex chemistry and pore structures of rocks. Mechanisms at the macroscale are affected by the processes occurring at the nanoscale. This review highlights the necessity of multitechnique, multiscale characterization of rocks and points to the importance of surface analysis and surface science studies. Two shale rocks (seals) from Niobrara and Agardhfjellet formations with complex surface chemistry are used as examples throughout the paper. Typically, evaluation of rocks with x-ray diffraction, thermogravimetric analysis, Rock-Eval pyrolysis, gas adsorption, and electron microscopy combined with energy dispersive x-ray spectroscopy is conducted to provide valuable information about the bulk mineralogy, elemental composition, pore volume, and adsorbed species on the sample surface. These studies are necessary prior to designing surface sensitive experiments with x-ray photoelectron spectroscopy (XPS), guiding both sample preparation and sample analysis. XPS has been widely used to study the surface composition of rocks during the investigations of their fine-scale wettability, and the main findings are highlighted here. This paper also reviews the existing literature on ambient-pressure XPS, which provides new opportunities to study in situ chemical alteration due to interactions with CO 2 and offers recommendations for adapting this technique to study rock-fluid interactions, allowing for the identification of fundamental interactions during CO 2 sequestration and guide selection of formation sites for improved storage.

54 ENVIRONMENTAL SCIENCES↗

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗

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↗