Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core tools”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Sustainable production of plastic degrading enzymes

The 21st century faces the challenge of transitioning to a climate-neutral economy. This hurdle will need to be addressed in a myriad of approaches, including all societal spheres, as highlighted in several UN and EU reports, including new and accessed technologies. Among them, biotechnology can play a core role by exploiting its environmental benefits (e.g., substituting fossil-based products and expected greenhouse emission reductions). In some cases, biotech can even outperform fossil counterparts' features (e.g., the biodegradability of some bioplastics). Microalga is an attractive bioplastic feedstock in this scenario since it directly substitutes petroleum and captures CO2. However, albeit its environmental benefits, the microalgal process still needs to be tipped economically and the net-positive ecological impact improved to drive adoption. To that end, we envisioned the use of a green alga as a source of bioplastic feedstock and industrial enzymes, combining natural traits of microalgae with engineered properties in a resourceful process to petroleum-based products replacement. As proof of concept, we demonstrated that we could use C. reinhardtii, a green algae, to produce plastic-degrading enzymes. We assemble a vector to secrete PHL7, a recently described PETase, capable of degrading post-consumer PET plastic. We successfully engineered the cell line to secrete the plastic-degrading enzyme. We confirmed the secretion of the enzyme by observing halos in agar plates supplemented with plastic dispersions, enzymatic assays, zymograms, and protein sequencing of supernatant samples. This approach allows using green algae to replace petrochemical-derived plastic while producing the tool necessary to recycle existing plastics sustainably.

Dutra Molino, Joao Vitor↗
DOE OSTI2023

Colorimetric Chemosensor for Cu2+ and Fe3+ Based on a meso-Triphenylamine-BODIPY Derivative

Optical chemosensors are a practical tool for the detection and quantification of important analytes in biological and environmental fields, such as Cu2+ and Fe3+. To the best of our knowledge, a BODIPY derivative capable of detecting Cu2+ and Fe3+ simultaneously through a colorimetric response has not yet been described in the literature. In this work, a meso-triphenylamine-BODIPY derivative is reported for the highly selective detection of Cu2+ and Fe3+. In the preliminary chemosensing study, this compound showed a significant color change from yellow to blue–green in the presence of Cu2+ and Fe3+. With only one equivalent of cation, a change in the absorption band of the compound and the appearance of a new band around 700 nm were observed. Furthermore, only 10 equivalents of Cu2+/Fe3+ were needed to reach the absorption plateau in the UV-visible titrations. Compound 1 showed excellent sensitivity toward Cu2+ and Fe3+ detection, with LODs of 0.63 µM and 1.06 µM, respectively. The binding constant calculation indicated a strong complexation between compound 1 and Cu2+/Fe3+ ions. The 1H and 19F NMR titrations showed that an increasing concentration of cations induced a broadening and shifting of the aromatic region peaks, as well as the disappearance of the original fluorine peaks of the BODIPY core, which suggests that the ligand–metal (1:2) interaction may occur through the triphenylamino group and the BODIPY core.

Pinto, Sónia C. S.↗
DOE OSTI2020

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Porous Colloidal Nanoparticles as Injectable Multimodal Contrast Agents for Enhanced Geophysical Sensing

Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.

NMR↗
DOE OSTI2025

A model to assess Zircaloy’s mechanical property changes following a transient beyond critical heat flux

Maintaining the integrity of nuclear fuel rods is essential for ensuring public health and safety in nuclear power generation. During reactor operation, this integrity is confirmed by demonstrating compliance with established regulatory acceptance criteria. For moderate-frequency events, such as limiting transients and anticipated operational occurrences (AOOs), the current fuel integrity criterion is based on preventing boiling transition. This criterion assumes that prevention of boiling transition will prevent excessive cladding heating and, thus, fuel failure during normal operations. While conservative, this approach places significant constraints on core design, fuel cycle economics, and a plant’s ability to perform major power uprates, leading to suboptimal fuel utilization and inefficient carbon-free energy production. A more efficient approach could be achieved by revising the failure criterion to a material-specific limit rather than strictly preventing the boiling transition, since boiling transition per se is not a cause of fuel cladding failure. Here, as a result, a new licensing framework based on material properties, termed time-at-temperature (t@T), is needed. This approach would allow for brief periods of post–critical heat flux operation during an AOO without compromising safety. Implementing the t@T licensing strategy requires a robust technical foundation in material properties, which must be established through comprehensive data collection on both unirradiated and irradiated fuel and cladding materials. This foundation would enable the development of a safety basis that ensures safe operation while providing greater flexibility and efficiency for reactor operation. This paper documents a thorough review of the available data to establish a baseline knowledge that can inform the development of cladding mechanical models, as well as identify experimental data gaps that need to be addressed in future research. Machine learning and data informatics were utilized to extract the importance of parameters on the t@T parameter. Industry tools were used to perform baseline analyses to define the relevant transient conditions for data analysis. The subsequent review successfully identified applicable experimental data, as well as sufficient data to evaluate changes in cladding mechanical properties following an AOO transient. Rather than developing new models, this work coupled existing irradiation annealing and recrystallization models to calculate changes in hardness, yield stress, and ultimate tensile stress following an AOO event. The findings from this review were summarized to highlight the experimental data needs required to fill remaining gaps and support the development of future t@T licensing methodologies.

Cladding performance↗
DOE OSTI2021

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Progress toward Accelogic compression in ROOT

For the last 7 years, Accelogic pioneered and perfected a radically new theory of numerical computing codenamed “Compressive Computing”, which has an extremely profound impact on real-world computer science [1]. At the core of this new theory is the discovery of one of its fundamental theorems which states that, under very general conditions, the vast majority (typically between 70% and 80%) of the bits used in modern large-scale numerical computations are absolutely irrelevant for the accuracy of the end result. This theory of Compressive Computing provides mechanisms able to identify (with high intelligence and surgical accuracy) the number of bits (i.e., the precision) that can be used to represent numbers without affecting the substance of the end results, as they are computed and vary in real time. The bottom-line outcome will be to provide state-of-the-art compression algorithms --and accompanying software libraries-- able to surpass the performance of the compression engines currently available in the ROOT [7] framework. The resulting technology has the capability to enable substantial economic and operational gains (including speedup) for High Energy and Nuclear Physics data storage/analysis. In our initial studies, a factor of nearly x4 (3.9) compression was achieved with RHIC/STAR data where ROOT compression managed only x1.4 [6].As a collaboration of experimental scientists, private industry, and the ROOT Team, our aim is to capitalize on the substantial success delivered by the initial effort and produce a robust technology properly packaged as an open-source tool that could be used by virtually every experiment around the world as means for improving data management and accessibility.In this contribution, we will present our efforts integrating our concepts of “functionally lossless compression” within the ROOT framework implementation, with the purpose of producing a basic solution readily integrated into HENP applications. We will also present our progress applying this compression through realistic examples of analysis from both the STAR and CMS experiments.

Canal, Ph.↗
DOE OSTI2025

Novel application of neutrinos to evaluate U.S. nuclear weapons performance

There is a growing realization that neutrinos can be used as a diagnostic tool to better understand the inner workings of a nuclear weapon. Robust estimates demonstrate that an Inverse Beta Decay (IBD) neutrino scintillation detector built at the Nevada Test Site with a 1000-ton active target mass at a standoff distance of 500 m would detect thousands of antineutrino events per nuclear test. This would provide less than 4% statistical error on the measured antineutrino rate and 5% error on antineutrino energy. Extrapolating this to an error on the test device explosive yield requires knowledge from evaluated nuclear databases, non-equilibrium fission rates, and assumptions on internal neutron fluxes. Initial calculations demonstrate that the total number of neutrinos emitted per fission in the first 10 3 s after a short pulse of 239 Pu fission is about a factor of two less than that from Pu fissioning under steady state conditions. Furthermore, there are significant energy spectral differences as a function of time after the pulse that must be considered. These and other model dependencies will be discussed in the paper. In the absence of nuclear weapons testing, many of the technical and theoretical challenges of a full nuclear test could be mitigated with a low cost smaller scale 20 ton fiducial mass IBD demonstration detector placed near a pulsed reactor. Potential reactors include the Texas A&M University TRIGA 1 GW–10 ms pulsed facility or the Sandia Annular Core Research Reactor. The short duty cycle and repeatability of pulses would provide critical real environment testing and measurements, which would be valuable for planning a possible real test shot in the future. Furthermore, the antineutrino rate as a function of time data would provide unique constraints on fission databases and model assumptions. Finally, there are impactful science drivers such as sensitive searches for ∼1 eV 2 sterile neutrinos and ∼MeV scale axions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Explosive Soot Challenge (Final Report)

This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Automated and Accelerated Continuum Model Development for Electrochemical Systems (Abbreviated Report)

Despite the availability of computational resources and advancements in numerical computing capabilities, the multiscale models core to understanding, predicting the behaviors of, and designing energy and environmental systems involving porous media are still 1.) developed through by-hand derivations and 2.) limited by many methodological assumptions employed during model derivation. As a result, the advancement of effective media models for engineering DOE mission-critical systems (e.g., batteries, flow batteries, electrolyzers, geothermal systems, subsurface chemical storage systems, etc.) is slow (i.e., it takes years for models to traverse from stages of “development” to “practical utilization”), hindering our ability to effectively optimize such systems and stay at the cutting-edge of the energy frontier. In this work, we aimed to address these limitations by 1.) automating and accelerating multiscale model derivation via symbolic computing and 2.) develop a novel multiscale modeling methodology for flow and transport through porous media that avoids the typical assumptions hindering previous models. As a result of our efforts, we 1.) developed a hybrid symbolic-numeric code called Fouriera for fully-automating the implementation of multiphysical and phase-field models via the Fourier spectral method for materials science research, and 2.) advanced a multiscale modeling methodology called The Method of Finite Averages that rigorously predicts the behaviors of flow and transport through heterogeneous porous media under the influence of non-local effects and strong advection. Ultimately, these deliverables provide strong foundations from which further efforts can advance multiscale modeling tools and capabilities that do not intrinsically rely on 1.) the speed and mathematical capabilities of humans, nor 2.) the methodological assumptions limiting current models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quarkyonic matter equation of state in beta-equilibrium

Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasiparticles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass (≳ 2 $M_⊙$) and the radius of a 1.4 $M_⊙$ star ($R_{1.4}$ ≲ 13.5 km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. Here, we propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Refining Processing Engines from SAPHIRE: Initialization of Fault Tree/Event Tree Solver

SAPHIRE has been extensively employed for over 35 years to model risk-important systems and quantify risk models. As a well-established and thoroughly documented Probabilistic Risk Assessment (PRA) tool, SAPHIRE has continuously tracked computational trends and received regular updates. Despite its ongoing evolution, there remains a need for further enhancements, particularly in dealing with the quantification of exceptionally large models. These improvements could take the form of algorithmic advancements, harnessing the power of parallel computing, and exploring the potential benefits of cloud computing solutions. Considering these aspirations, the notion of a remote solve option was introduced and subsequently evolved into a dedicated project within the SAPHIRE development team. A significant outcome of this initiative is SAPHSOLVE, an engine extracted from the SageRisk API designed specifically for remote solving capabilities. The ongoing project is nearing its culmination, marked by a series of discoveries that have brought undocumented aspects to light. Among these revelations is the intricacy of the input and output format for the SAPHSOLVE engine. This document serves the crucial purpose of meticulously delineating the precise formats for both input and output, as they form an indispensable foundation. The importance of documenting these formats cannot be overstated, as it is a pivotal step in facilitating rigorous testing and comparison. Whether it involves scrutinizing SAPHSOLVE results against those of the internal integrated solver or other external solvers, the ability to construct models or transform existing ones into a compatible SAPHSOLVE format is imperative. Chapter 1 offers a succinct introduction to both SAPHIRE and SAPHSOLVE, followed by Chapter 2 which outlines the roadmap for enhancing SAPHSOLVE. The core of this report is Chapter 3, which intricately elucidates the intricacies of the input and output file formats. To provide a tangible illustration of these formats, a rudimentary example has been compiled and is available in Appendix. SAPHSOLVE represents a novel external solving mechanism developed by the SAPHIRE team, although it has not yet reached the full spectrum of capabilities possessed by SAPHIRE's internal solver. However, the SAPHIRE team has set a comprehensive course for incorporating the functionalities of SAPHSOLVE. A comprehensive outlook on the future of SAPHSOLVE is expounded upon in Chapter 4.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Advanced Transmission Technologies – GETs and HPCs Session 2: Advanced Power Flow Control and Transmission Topology Optimization

The INL TADA GETs Cohort Session 2, held on November 7, 2025, conducted in collaboration with ScottMadden, focused on two core Advanced Transmission Technologies (ATTs): Advanced Power Flow Control (APFC) and Transmission Topology Optimization (TTO). These technologies are pivotal in enhancing grid flexibility, reliability, and cybersecurity resilience. APFC, particularly through modular FACTS devices like Modular Static Synchronous Series Compensators (M-SSSCs), enables dynamic voltage injection to reroute power flows. The session highlighted the deployment benefits of APFC, such as rapid installation, minimal civil works, and re-deployability. Regulatory drivers like FERC Order 2023 mandate the inclusion of Grid-Enhancing Technologies (GETs) in interconnection studies. Case studies from Central Hudson, CAISO, and National Grid (UK) demonstrated APFC’s effectiveness in congestion relief and cost savings. The session also addressed cybersecurity concerns, including firmware vulnerabilities, SCADA integration risks, and supply chain dependencies. Participants engaged in interactive exercises to rank cybersecurity and supply chain risks, emphasizing the need for robust digital assurance strategies. TTO involves software-based reconfiguration of transmission networks to optimize power flow without new infrastructure. The session showcased its operational value, with examples from SPP, PJM, and MISO showing significant congestion cost reductions. Cybersecurity vulnerabilities were discussed, particularly in API security and software supply chains, referencing incidents like SolarWinds and attacks on Danish utilities. Digital assurance exercises explored worst-case scenarios, attack paths, and mitigation responsibilities between vendors and utilities. Reliability challenges such as algorithm stability, vendor dependency, and operator trust were also examined. Cross-cutting themes emphasized the importance of digital assurance tools, including Software Bills of Materials (SBOMs) and hardware-in-loop testing. Human performance, training, and operational confidence were identified as critical enablers of technology adoption. The session concluded with a preview of Session 3, which will focus on High Performance Conductors (HPCs) and risk-based cybersecurity tools. Session 2 of 3.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Environmental Protection Coating System for Refractory Metal Alloys (EPCS for RMAs)

As part of ARPA-E’s Ultra High Temperature Impervious Materials Advancing Turbine Efficiency (UITIMATE) program, this project aims to develop an environmental protection coating system (EPCS) for refractory metal alloys (RMAs) – which will provide radical improvement in long term protection for ultra-high temperature refractory metal components in harsh gas-turbine environments. The drive for higher fuel efficiency and higher core power of gas turbines used in electric power generation and in aircraft propulsion requires higher peak operation temperatures. The temperature capability of state-of-the-art materials cannot meet the new requirements. Refractory metal alloys have desirable capabilities for achieving the required peak temperatures, but most of them are subject to oxidation in the gas turbine environment. Thus, there is a strong need to develop novel environmental protective coatings that enable the refractory metals to operate in an ultra-high temperature engine environment. This project is composed of four major technical innovations that can provide a potentially transformational solution in providing environmental protection for refractory metal alloys over a wide range of temperature. 1) Multi-layer self-healing environmental barrier coating (EBC) that provides protection against combustion gas environment. 2) Oxidation resistant diffusion barrier (DB) on the alloy surface to enhance coating stability. 3) Synergistic integration of the multi-layer self-healing EBC with the diffusion barrier to extend coating life. 4) Guidance of physics-based modeling to accelerate coating development and optimization. Our team has successfully developed a multi-layer self-healing environmental protection coating system (EPCS) that can provide oxidation protection for commercially available molybdenum based refractory metal alloys over a wide temperature range of up to 1600 °C for fuel-efficient gas turbine applications. The team also developed a multi-layer high temperature coating design tool to accelerate coating development and optimization based on physics-based modeling. The team also engaged RTX business unite, notable Pratt and Whitney, for future commercialization of this technology. The successful completion of this Phase I effort led to the selection of performing Phase II of the program. In Phase II, the team plans to apply the EPCS in the protection of a newly discovered molybdenum-based metal alloy which has the potential of superior high temperature properties to the commercial products. Th coating will be evaluated under test conditions and environment that mimic gas turbine operation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multi-Scale Characterization of Porosity and Cracks in Silicon Carbide Cladding after Transient Reactor Test Facility Irradiation

Silicon carbide (SiC) ceramic matrix composite (CMC) cladding is currently being pursued as one of the leading candidates for accident-tolerant fuel (ATF) cladding for light water reactor applications. The morphology of fabrication defects, including the size and shape of voids, is one of the key challenges that impacts cladding performance and guarantees reactor safety. Therefore, quantification of defects’ size, location, distribution, and leak paths is critical to determining SiC CMC in-core performance. This research aims to provide quantitative insight into the defect’s distribution under multi-scale characterization at different length scales before and after different Transient Reactor Test Facility (TREAT) irradiation tests. A non-destructive multi-scale evaluation of irradiated SiC will help to assess critical microstructural defects from production and/or experimental testing to better understand and predict overall cladding performance. X-ray computed tomography (XCT), a non-destructive, data-rich characterization technique, is combined with lower length scale electronic microscopic characterization, which provides microscale morphology and structural characterization. This paper discusses a fully automatic workflow to detect and analyze SiC-SiC defects using image processing techniques on 3D X-ray images. Following the XCT data analysis, advanced characterizations from focused ion beam (FIB) and transmission electron microscopy (TEM) were conducted to verify the findings from the XCT data, especially quantitative results from local nano-scale TEM 3D tomography data, which were utilized to complement the 3D XCT results. In this work, three SiC samples (two irradiated and one unirradiated) provided by General Atomics are investigated. The irradiated samples were irradiated in a way that was expected to induce cracking, and indeed, the automated workflow developed in this work was able to successfully identify and characterize the defects formation in the irradiated samples while detecting no observed cracking in the unirradiated sample. These results demonstrate the value of automated XCT tools to better understand the damage and damage propagation in SiC-SiC structures for nuclear applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗