SEARCH · Engineering Papers

Results for “parallel time integration”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Assessing parallel path cooling tower performance via artificial neural networks

Real-time monitoring of a research nuclear reactor, a system in which all generated power is dissipated to the environment, can be performed via analysis of the heat rejection from the cooling system. Given an inlet water temperature and flow rate, the reactor power can be well-approximated from the outlet water temperature; however, the instrumentation to measure outlet conditions may not be robust or accurate. If we know how a cooling tower performs from historical data, but cannot measure the outlet temperature, a mathematical representation of the system can be inverted to obtain the outlet water temperature that describes the cooling capacity. Unfortunately, model inversion processes are computationally expensive. To address this, an artificial neural network (ANN) is implemented to assess the performance of a multi-cell cooling tower for a nuclear reactor. This approach leverages the Merkel model to obtain an extensive data set describing performance of the cooling tower cells throughout a wide array of potential operating conditions. The Merkel model is expressed as a function of four parameters: the inlet and outlet water temperatures, inlet air wet bulb temperature, and ratio of liquid-to-gas mass flow rates (L/G), which together provide a non-dimensional number indicative of cooling tower performance, called the Merkel integral. Computing a 4-dimensional data structure that describes finite combinations of the Merkel integral, an inverse model is then generated using an ANN to determine the cell outlet water temperature from the other three model parameters along with the computed Merkel integral. Compared to traditional model inversion methods, the ANN reduces the computational time by approximately 4 orders of magnitude, with effectively no sacrifice to solution accuracy, and could be applied for different cooling towers in the event the performance curve is known. Finally, three use cases of the ANN are then reviewed: (1) determining the cell outlet water temperatures when gas flow at rated conditions (GFRC) is known, (2) performing the prior case without knowledge of the GRFC, and (3) assessing performance differences between the individual tower cells.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Merged Observatory Data Files (MODFs): an integrated observational data product supporting process-oriented investigations and diagnostics

A large and ever-growing body of geophysical information is measured in campaigns and at specialized observatories as a part of scientific expeditions and experiments. These collections of observed data include many essential climate variables (as defined by the Global Climate Observing System) but are often distinguished by a wide range of additional non-routine measurements that are designed to not only document the state of the environment but also the drivers that contribute to that state. These field data are used not only to further understand environmental processes through observation-based studies but also to provide baseline data to test model performance and to codify understanding to improve predictive capabilities. To address the considerable barriers and difficulty in utilizing these diverse and complex data for observation–model research, the Merged Observatory Data File (MODF) concept has been developed. A MODF combines measurements from multiple instruments into a single file that complies with well-established data format and metadata practices and has been designed to parallel the development of corresponding Merged Model Data Files (MMDFs). Using the MODF and MMDF protocols will facilitate the evolution of model intercomparison projects into model intercomparison and improvement projects by putting observation and model data “on the same page” in a timely manner. The MODF concept was developed especially for weather forecast model studies in the Arctic. The surprisingly complex process of implementing MODFs in that context refined the concept itself. Thus, this article explains the concept of MODFs by providing details on the issues that were revealed and resolved during that first specific implementation. Detailed instructions are provided on how to make MODFs, and this article can be considered a MODF creation manual.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Investigation of pedestal parameters and divertor heat fluxes in small ELM regimes in DIII-D

Abstract Divertor heat flux and its correlation with pedestal parameters within various small edge localized mode (ELM) regimes, including high beta poloidal, type-II and ELMs with negative triangularity H-modes were investigated in DIII-D. The parallel energy fluences of type-II and high beta poloidal small ELM regimes fall below the linear scaling with pedestal electron pressure for type-I ELMs put forward in Eich et al 2017 ( Nucl. Mater. Energy 12 84–90). The negative triangularity of H-mode ELMs follow the Eich scaling for type-I ELMs. The parallel heat flux and total heat loads to the divertor were determined using high-time resolution infrared thermography, while pedestal parameters were obtained through self-consistent kinetic equilibrium reconstructions. Linear regressions for the type-II and high beta poloidal regimes demonstrate that an equivalent 7.5 MA small ELM scenario in ITER would fall below the ~5 MJ m − 2 leading edge melting limit for tungsten (Gunn et al 2017 Nucl. Fusion 57 046025). Utilizing fast thermography, the scrape-off layer power fall-off length for both inter-ELM and intra-ELM was determined and compared to the Eich scaling with poloidal magnetic field in Eich et al (ASDEX Upgrade Team and JET EFDA Contributors 2013 Nucl. Fusion 53 093031). Except for the high beta poloidal scenario, all the small ELM regimes during both inter- and intra-ELM periods had power fall-off lengths ( λ q ) larger then would be expected from the B pol , MP − 1 scaling associated with type-I ELMs, signifying their potential in managing heat loads and offering a solution for core–edge integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Regional-scale fault-to-structure earthquake simulations with the EQSIM framework: Workflow maturation and computational performance on GPU-accelerated exascale platforms

Continuous advancements in scientific and engineering understanding of earthquake phenomena, combined with the associated development of representative physics-based models, is providing a foundation for high-performance, fault-to-structure earthquake simulations. However, regional-scale applications of high-performance models have been challenged by the computational requirements at the resolutions required for engineering risk assessments. The EarthQuake SIMulation (EQSIM) framework, a software application development under the US Department of Energy (DOE) Exascale Computing Project, is focused on overcoming the existing computational barriers and enabling routine regional-scale simulations at resolutions relevant to a breadth of engineered systems. This multidisciplinary software development—drawing upon expertise in geophysics, engineering, applied math and computer science—is preparing the advanced computational workflow necessary to fully exploit the DOE’s exaflop computer platforms coming online in the 2023 to 2024 timeframe. Achievement of the computational performance required for high-resolution regional models containing upward of hundreds of billions to trillions of model grid points requires numerical efficiency in every phase of a regional simulation. This includes run time start-up and regional model generation, effective distribution of the computational workload across thousands of computer nodes, efficient coupling of regional geophysics and local engineering models, and application-tailored highly efficient transfer, storage, and interrogation of very large volumes of simulation data. This article summarizes the most recent advancements and refinements incorporated in the workflow design for the EQSIM integrated fault-to-structure framework, which are based on extensive numerical testing across multiple graphics processing unit (GPU)-accelerated platforms, and demonstrates the computational performance achieved on the world’s first exaflop computer platform through representative regional-scale earthquake simulations for the San Francisco Bay Area in California, USA.

58 GEOSCIENCES↗

DOE OSTI2022

Manufacturing Supply Chain Development for Modular Solar-Thermochemical Conversion Platform - CRADA 387 (Final Report)

Modular chemical process intensification (MCPI) is an emerging field where chemical processing is performed using small-scale modular equipment instead of conventional large centralized chemical plants. Conventional chemical plants benefit from economies of scale that encourage scale-up to ever larger plants. A goal of MCPI is to develop technology that intensifies processing so that equipment can be dramatically smaller and integrated into modular systems. Scale-up occurs by adding more modules in parallel rather than making the equipment larger. A key concept is that equipment and modules can ultimately be cheaper by leveraging economies of mass production, analogous to the automotive industry, in manufacturing the equipment. This project made significant progress toward this outcome by meeting the RAPID institute metric to reduce equipment cost by 20% for each doubling in manufacturing volume. The MCPI application was thermochemical technology that is being commercialized by STARS Technology Corporation, one of the CRADA partners. The technology converts solar and renewable power to chemical energy to produce renewable hydrogen, fuels, and chemicals. The benefit to the public is reduction in greenhouse gases that are contributing to climate change. The project transitioned the steam methane reforming (SMR) reactor from conventional fabrication methods to additive manufacturing (AM) direct metal laser sintering (DMLS) process. This is projected to reduce the cost of making a reactor by 58% when producing 100 reactors per year. Innovations in the DMLS process produced a patented design that reduces reactor weight by 60%. Reductions in material costs and processing time extend the DMLS advantage to higher production volumes. The new design promises to be 38% cheaper than the conventional processes at 1000 units per year. The resulting 87% reduction in the steam methane reforming (SMR) module cost in scaling from current costs meets the RAPID metric. The project was successful in producing and testing the first ever additively manufactured SMR reactors. A reactor achieved over 82% efficiency in converting electric power to chemical energy, which is a world record for an inductively heated SMR. The project has contributed to the design and assembly of a first demonstration plant that is headed to a hydrogen bus filling station in Thousand Palms, CA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Temporal coherent control of resonant two-photon double ionization of the hydrogen molecule via doubly excited states

Here, we use time-delayed, counter-rotating, circularly polarized few-cycle attosecond nonoverlapping pulses to study the temporal coherent control of the resonant process of two-photon double ionization (TPDI) of hydrogen molecule via doubly excited states for pulse propagation direction along $\hat k$ either parallel or perpendicular to the molecular axis $\hat R$. For $\hat k$∥ $\hat R$ and a pulse carrier frequency of 36 eV resonantly populating the Q 2 1 Π$^+_u$ (1) doubly excited state as well as other 1 Π$^+_u$ doubly excited states, we find that the indirect ionization pathway through these doubly excited states changes the character of the kinematical vortex-shaped momentum distribution produced by the two direct ionization pathways from fourfold to twofold rotational symmetry. This result is similar to what found in TPDI of the He atom involving 1 P$^ o_ {±1}$ doubly excited states; however, angular distributions exhibiting a quantum beat effect between the ground state and a doubly excited state seen for the He atom are observed here for its molecular counterpart with an anomaly in shape and magnitude, not in frequency. The sixfold differential probability integrated over the azimuthal angle of the photoelectron pair shows that this anomaly is due to autoionization decays and quantum beats between doubly excited states. For $\hat k$⟂$\hat R$ and a broadband pulse carrier frequency of 30 eV populating the Q 1 1 Π$^+_u$(1), Q 1 1 Σ$^+_u$(1), Q 2 1 Π$^+_u$(1) , and Q 1 1 Σ$^+_u$ (2) doubly excited states, the momentum distribution is shown to exhibit dynamical electron vortices with four spiral arms, which originates from the interplay between the 1 Δ$^+_ g$ , 1 Π$^+_ g$, and 1 Σ$^+_ g$ dynamical ionization amplitudes. Our treatment within either the adiabatic-nuclei approximation or fixed-nuclei approximation shows that the latter provides a very good account for this correlated process.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

National Ignition Facility Opacity Time Resolved Spectrometer Systems Engineering Final Project

The National Ignition Facility (NIF) is the world’s largest and most energetic laser facility. The NIF system is designed to produce high energy density (temperature and pressure) conditions through the application of its 192 laser beams. One of the users of NIF is the opacity platform developed to study the opacities at temperatures and densities relevant to the solar interior and stellar evolution. The platform was developed to study iron (Fe) opacity at temperatures relevant to the solar interior. The opacity campaign uses spectrometers to gather data. Spectrometers utilize crystals to produce x-ray spectra that are recorded on time-integrated and time-resolved detectors. The opacity spectrometer (OpSpec) currently fielded and in use at NIF uses a time integrated film channel to collect data. The opacity spectrometer time resolved (OpSpecTR) will utilize novel hCMOS detectors to capture time resolved images of spectra of interest. The key stakeholders identified for OpSpecTR included the physicists responsible for OpSpec and OpSpecTR, the Target Area Science and Engineering (TASE) department at NIF, the NIF and Photon Science (NIF & PS) Opacity program, the Nevada National Security Site (NNSS) Physics and Engineering program, the Sandia hCMOS manufacturing and testing program, and the Los Alamos National Laboratory (LANL) program sponsor. The Target and Experimental Operations (TEXOPS) was identified as a key stakeholder because the group includes the individuals that will physically interact with the OpSpecTR system as it participates in NIF experiments. The opacity platform collects data in a unique orientation relative to the existing diagnostics fielded at NIF. The existing infrastructure at NIF uses a diagnostic manipulator (DIM) to insert the diagnostic near the target chamber center to collect data during a NIF shot. Existing diagnostics collect data through the center line of the DIM axis and collect relevant data perpendicular to this axis. The opacity platform requires crystals mounted in a specific orientation which requires data collection parallel to the DIM axis. This deviation from standard NIF practices was a key factor in developing requirements.

42 ENGINEERING↗

DOE OSTI2020

Nitrogen-seeded divertor detachment in TCV L-mode plasmas

Most of the detachment experiments done to date on the tokamak à configuration variable (TCV), both in standard and alternative divertor geometries, focused on L-mode integrated core density ramps. In view of extending these studies to high-power, high-confinement regimes, where impurity seeding will be necessary for detachment, the properties of nitrogen seeded L-mode detachment in TCV are assessed here with the extensive set of edge and divertor diagnostics and similarities and differences with integrated core density ($\langle n_e \rangle$) ramp detachment experiments are elucidated. It is found that in high current, reversed field plasmas, detachment at the outer target is achieved with N 2 -seeding and density ramps, with target heat flux reductions of up to 90%, while the inner target only detaches with seeding. The Scrape-Off Layer radiation fraction reaches values of 60%–80% and in all situations, a stable radiator can form around the X-point. In this work, the most striking difference between seeding and density ramp is the behavior of the upstream quantities. During the $\langle n_e \rangle$-ramp, a broadening of the upstream density profile (density ‘shoulder’) occurs, concurrent with the outer target ion flux roll-over, while no such behavior occurs during nitrogen seeded detachment. Separatrix density, electron temperature and pressure also evolve strongly with increasing density, and are largely unaffected by the injection of nitrogen. Comparison of upstream and target pressures reveals that, in all cases, the outer target ion flux reduction coincides with the development of a parallel gradient of the total pressure. Common to all cases is also a reduction of energy confinement time with detachment, although this effect is weak for seeding at relatively high density. Studying the impact of the ∇B-drift direction in both nitrogen seeding and core density ramps reveals that drifts mainly affect the behavior at the inner strike point, highlighting the need to include them in edge transport simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Scalable Data-Intensive Geocomputation: A Design for Real-Time Continental Flood Inundation Mapping

The convergence of data-intensive and extreme-scale computing enables an integrated software and data ecosystem for scientific discovery. Developments in this realm will fuel transformative research in data-driven interdisciplinary domains. Geocomputation provides computing paradigms in Geographic Information Systems (GIS) for interactive computing of geographic data, processes, models, and maps. Because GIS is data-driven, the computational scalability of a geocomputation workflow is directly related to the scale of the GIS data layers, their resolution and extent, as well as the velocity of the geo-located data streams to be processed. Unique in high user interactivity and low end-to-end latency requirements, geocomputation applications will dramatically benefit from the convergence of high-end data analytics (HDA) and high-performance computing (HPC). The application level challenge, however, is to identify and eliminate computational bottlenecks that arise along a geocomputation workflow. Indeed, poor scalability at any of the workflow components is detrimental to the entire end-to-end pipeline. Here, we study a large geocomputation use case in flood inundation mapping that handles multiple national-scale geospatial datasets and targets low end-to-end latency. We discuss benefits and challenges for harnessing both HDA and HPC for data-intensive geospatial data processing and intensive numerical modeling of geographic processes. We propose an HDA+HPC geocomputation architecture design that couples HDA (e.g., Spark)-based spatial data handling and HPC-based parallel data modeling. Key techniques for coupling HDA and HPC to bridge the two different software stacks are reviewed and discussed.

Liu, Yan↗

DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

ACHILLES-GENIE Interface for Neutrino Simulations, ADRIANO2 Tile Prototype for High-Granularity Dual-Readout Calorimetry

ICARUS (Imaging Cosmic And Rare Underground Signals) is a liquid argon time projection chamber (LArTPC) detector that pursues the sterile neutrino, which relies on accurate simulations of neutrino-argon interactions. REDTOP (Rare Eta Decays To Observe new Physics) is a proposed low-energy, high-intensity meson factory designed to explore rare $\eta$/$\eta'$ meson decays and probe physics beyond the Standard Model. As a next-generation experiment, this requires both accurate simulations and innovative detector technologies. This project contributes to both ICARUS, from a simulation perspective, and REDTOP, from both a simulation and detection perspective, through the event generation of lepton-nucleon interactions and the physical enhancement of the calorimeter technology within the REDTOP detector. We developed an interface between ACHILLES (A CHIcago Land Lepton Event Simulator), a theory-driven lepton-level event generator, and GENIE, a robust event generator framework used for neutrino physics. By incorporating the precise theoretical cross-section calculations of ACHILLES into the experimental realism of GENIE, the interface allows for improved accuracy of neutrino-nucleon simulations, which can be adapted for the proton beam specifications of the REDTOP meson factory as well as for the ICARUS experiment. In parallel, we developed an improved prototype for the ADRIANO2 (A Dual Readout Integrally Active Non-segmented Option) dual-readout calorimeter tiles for the REDTOP detector. To improve the efficiency of the lead-glass tiles trapping Cherenkov light for energy reconstruction and particle identification, we optimized the application of a highly reflective coating. Through viscosity and thickness control, masking, and a custom spray technique, we refined the coating process to reduce surface defects and improve light yield. Together, these efforts strengthen the ICARUS neutrino program and REDTOP's capability of detecting rare decay events.

Visser, Erin [Michigan State U.]↗

DOE OSTI2025

Enhancing ICARUS and REDTOP Software and Hardware: Event Generator Interface Development and Calorimeter Tile Prototype

Visser, Erin [Michigan State U.] (ORCID:0009000184↗

DOE OSTI2024

Wire arc additive manufactured A36 steel performance for marine renewable energy systems

Additive manufacturing has established itself to be advantageous beyond small-scale prototyping, now supporting full-scale production of components for a variety of applications. Despite its integration across industries, marine renewable energy technology is one largely untapped application with potential to bolster clean energy production on the global scale. Wave energy converters (WEC) are one specific facet within this realm that could benefit from AM. As such, wire arc additive manufacturing (WAAM) has been identified as a practical method to produce larger scale marine energy components by leveraging cost-effective and readily available A36 steel feedstock material. The flexibility associated with WAAM can benefit production of WEC by producing more complex structural geometries that are challenging to produce traditionally. Additionally, for large components where fine details are less critical, the high deposition rate of WAAM in comparison to traditional wrought techniques could reduce build times by an order of magnitude. In this context of building and supporting WEC, which experience harsh marine environments, an understanding of performance under large loads and corrosive environments must be understood. Hence, WAAM and wrought A36 steel tensile samples were manufactured, and mechanical properties compared under both dry and corroded conditions. Here, the unique microstructure created via the WAAM process was found to directly correlate to the increased ultimate tensile and yield strength compared to the wrought condition. Static corrosion testing in a simulated saltwater environment in parallel with electrochemical testing highlighted an outperformance of corroded WAAM A36 steel than wrought, despite having a slighter higher corrosion rate. Ultimately, this study shows how marine energy systems may benefit from additive manufacturing components and provides a foundation for future applications of WAAM A36 steel.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Data Mining and Visualization of High-Dimensional ICME Data for Additive Manufacturing

Integrated computational materials engineering (ICME) methods combining CALPHAD with process-based simulations can produce rich, high-dimensional data for alloy and process design. In ICME methods for metallurgical applications, the visualization and interpretation of such high-dimensional data has previously been through heat maps represented in 2 or 3 dimensions. While such an approach is ideal when one variable is varied at a time, in the case of high-dimensional data with multiple variables varied simultaneously, as is the case in additive manufacturing, interpreting the trends through two- or three-dimensional heat maps becomes challenging. Here, we propose a strategy of mixed visual data mining and quantitative analysis for high-dimensional metallurgical and process data using high-throughput thermodynamic calculations. Two case studies show the application of the proposed approach. The first case study investigated the effects of feedstock chemistry on the δ ferrite formation in 316L stainless steel powders used for binder jet additive manufacturing. The second case study linked Scheil–Gulliver calculations to a process model for dissimilar joining of aluminum alloys 5356 and 6111 during laser hot-wire additive manufacturing. Both cases contained thousands of calculated data points, showcasing the utility of visual data analysis through parallel coordinate plotting, Pearson correlation coefficient matrices, and scatter matrices compared to traditional process maps. These visualization techniques can be extended to many additive manufacturing problems to capture process–structure–property relationships for additively manufactured components.

36 MATERIALS SCIENCE↗

DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

DOE OSTI2026

PRISMA: PARALLEL REFINEMENT AND INTEGRATION SYSTEM FOR MULTI-AZIMUTHAL ANALYSIS

The Parallel Refinement and Integration System for Multi-azimuthal Analysis (PRISMA, version 1.1.0) is a Python application for processing X-ray diffraction (XRD) image data. PRISMA wraps GSAS-II to perform azimuthally-binned peak refinement, computes per-frame strain and d-spacing from those fits, and provides three PyQt5 graphical interfaces: (1) a Recipe Builder for selecting GSAS-II control (.imctrl) files, optional mask (.immask) files or threshold-ased masking, reference and experiment image sets, peaks, zimuthal range and bin size, and an optional ceria-based auto-calibration; (2) a Batch Processor that uses Dask on local workstations and pure MPI (mpi4py.futures.MPICommExecutor) on HPC to distribute GSAS-II refinement across cores or compute nodes and write results to a 4-dimensional (peaks x frames x azimuths x measurements) Zarr dataset; and (3) a Data Analyzer that renders heatmaps of fit parameters, strain, frame-to-frame deltas, and percent-change-vs-reference, and exports user-defined subsections to CSV or Excel. The peak-refinement algorithm is deterministic. Benchmark on ALCF Crux: a 20,000-image set, single-peak fit in frame mode with 44 azimuthal bins on 128 nodes x 128 workers, 48 seconds total wall time.

Lorenzo Martin, Maria De La Cinta [Argonne Nationa↗

DOE OSTI2023

Microgrid energy scheduling under uncertain extreme weather: Adaptation from parallelized reinforcement learning agents

Microgrids are useful solutions for integrating renewable energy resources and providing seamless green electricity to minimize carbon footprint. In recent years, extreme weather events happened often worldwide and caused significant economic and societal losses. Such events bring uncertainties to the microgrid energy scheduling problems and increase the challenges of microgrid operation. Traditional optimization approaches suffer from the inaccuracy of the uncertain microgrid model and the unseen events. Existing reinforcement learning (RL) - based approaches are also hampered by the limited generalization and the increasing computational burden when stochastic formulations are required to accommodate the uncertainties. This paper proposes a new parallelized reinforcement learning (PRL) method based on the probabilistic events to handle the microgrid energy uncertainties. Specifically, several local learning agents are employed to interact with pertinent microgrid environments in a distributed manner and report outcomes to the global agent, which will optimize microgrid energy resources online during extreme events. The stochastic microgrid energy optimization problem is reformulated to include all possible scenarios with probabilities. The advantage estimate functions of learning agents are designed with a backward sweep to transfer the outcomes to the value function updating process. Two simulation studies, stochastic optimization and online testing, are performed to compare with several existing RL approaches. Results substantiate that the proposed PRL method can achieve up to 20% optimization performance improvement with 4 and 28 times less computation cost than Q-learning with experience replay and multi-agent Q-learning approaches, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Fork Experiments in the Hot Cell Using Spent Fuel Rods for International Nuclear Safeguards

This work leveraged the rare availability of 25 full-length pressurized water reactor spent fuel rods and 1 irradiated mixed-oxide rod at an Oak Ridge National Laboratory hot cell. This was done to collect measurement data with two Fork detectors to assess the detectors’ capability of verifying operator declaration data and detecting partial defects in spent fuel, which are the two primary goals of international safeguards on spent nuclear fuel. The data can also be used to benchmark the ORIGEN module, which has been adopted in the International Atomic Energy Agency’s (IAEA’s)/European Atomic Energy Community’s (Euratom’s) Integrated Review and Analysis Program to predict the Fork detector count rates in real time. In this project, the authors first calibrated two Fork detectors—a standard one and a modified one—by using known strong neutron and gamma sources. Then, the authors measured all 26 fuel rods at multiple locations along the length. The fuel rods were then assembled into three arrays—2 × 2, 3 × 3, and 5 × 5—by using specially designed support grids to mimic fuel assemblies and measure the arrays with both detectors. For the 5 × 5 array, 4 and 8 fuel rods of the array were replaced in two separate cases with short stainless-steel rods to mimic two partial defect scenarios, and the arrays were measured before and after the replacements. Polyethylene blocks were used in this experiment to mimic water. The results show that the Fork detectors were able to verify operator declarations and detect partial defects in spent fuel, and the authors were the first to demonstrate this through experiments. A discovery was also made that determined the root cause of the nonlinear response to gamma dose in the ion chambers used in IAEA and Euratom’s Fork detectors. After the experiments, both detectors were retrieved from the hot cell for future use. The data collected in this project will be used in a parallel International Nuclear Safeguards Engagement Program (INSEP) project to enhance the safeguards in the Finnish spent fuel encapsulation plant, and the data will be useful to other projects in the future given the increased safeguards needs due to spent fuel transfer and disposal activities worldwide.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗