Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bending Magnet Synchrotron Radiation Imaging with Large Orbital Collection Angles

Synchrotron radiation (SR) from bending magnets, wigglers, and undulators is now extensively produced for users at storage ring based light sources, with unique properties in terms of average brightness and stability. Here, we present a profound study of bending magnet SR intensity distribution in the image plane of a focusing optical system. Measurements of this intensity distribution at the MAX-IV low emittance storage ring are compared to theoretical predictions, and found to be in excellent agreement. This work shows upon the possibility of performing high resolution emittance diagnostics with visible or near-visible SR on upcoming low-emittance storage ring based light sources. As a byproduct of our study, we derive a closed analytical expression for the intensity distribution from a zero-emittance beam, in the limiting case of wide orbital collection angles. This expression finally allows us to demonstrate the meeting between classical electrodynamics applied to SR emission and focusing, and the Landau and Lifshitz prediction of radiation intensity distribution nearby a caustic.

36 MATERIALS SCIENCE↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2026

Elucidating interfacial active sites in ruthenium–boron nitride nanotube catalysts for efficient low-temperature ammonia-to-hydrogen conversion

Tailoring the interaction between metal nanoparticles and catalyst support presents a prominent strategy to enhance both the activity and durability in hydrogen (H 2 ) production catalysts. In this work, ruthenium nanoparticles (NPs) supported on boron nitride nanotubes (Ru/BNNT) are introduced as efficient and thermally robust catalysts for low-temperature ammonia (NH 3 ) decomposition. The unique curvature and ionic nature of BNNTs enable uniform Ru dispersion and metal-support interactions (MSIs), resulting in exceptional H 2 generation efficiency and long-term operational stability. In-situ transmission electron microscopy (TEM) reveals remarkable thermal resistance of Ru/BNNT with minimal nanoparticle sintering, while density functional theory (DFT) calculations uncover a dual-site mechanism in which interfacial Ru atoms promote NH 3 dissociation and adjacent Ru sites facilitate 2H* recombination and H 2 desorption. This cooperative interaction between metal NPs and the BNNT support underpins the outstanding catalytic performance and durability observed. In conclusion, the findings highlight the strategic potential of BNNTs as versatile supports for high-performance and stable catalysts in sustainable H 2 energy conversion and related catalytic processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data for Microbial-Explicit Processes and Refined Perennial Plant Traits Improve Modeled Ecosystem Carbon Dynamics

Globally, soils hold approximately half of ecosystem carbon and can serve as a source or sink depending on climate, vegetation, management, and disturbance regimes. Understanding how soil carbon dynamics are influenced by these factors is essential to evaluate proposed natural climate solutions and policy regarding net ecosystem carbon balance. Soil microbes play a key role in both carbon fluxes and stabilization. However, biogeochemical models often do not specifically address microbial-explicit processes. Here, we incorporated microbial-explicit processes into the DayCent biogeochemical model to better represent large perennial grasses and mechanisms of soil carbon formation and stabilization. We also take advantage of recent model improvements to better represent perennial grass structural complexity and life-history traits. Specifically, this study focuses on: 1) a plant sub-model that represents perennial phenology and more refined plant chemistry with downstream implications for soil organic matter (SOM) cycling though litter inputs, 2) live and dead soil microbe pools that influence routing of carbon to physically protected and unprotected pools, 3) Michaelis-Menten kinetics rather than first-order kinetics in the soil decomposition calculations, and 4) feedbacks between decomposition and live microbial pools. We evaluated the performance of the plant sub-model and two SOM cycling sub-models, Michaelis-Menten (MM) and first-order (FO), using observations of net ecosystem production, ecosystem respiration, soil respiration, microbial biomass, and soil carbon from long-term bioenergy research plots in the mid-western United States. The MM sub-model represented seasonal dynamics of soil carbon fluxes better than the FO sub-model which consistently overestimated winter soil respiration. While both SOM sub-models were similarly calibrated to total, physically protected, and physically unprotected soil carbon measurements, the models differed in future soil carbon response to disturbance and climate, most notably in the protected pools. Adding microbial-explicit mechanisms of soil processes to ecosystem models will improve model predictions of ecosystem carbon balances but more data and research are necessary to validate disturbance and climate change responses and soil pool allocation.

Field Data↗
DOE OSTI2023

Pt Nanoparticles on Atomic-Metal-Rich Carbon for Heavy-Duty Fuel Cell Catalysts: Durability Enhancement and Degradation Behavior in Membrane Electrode Assemblies

Proton exchange membrane fuel cells (PEMFCs) are a promising zero-emission power source for heavy-duty vehicles (HDVs). However, long-term durability of up to 25,000 h is challenging because current carbon support, catalyst, membrane, and ionomer developed for traditional light-duty vehicles cannot meet the stringent requirement. Therefore, understanding catalyst degradation mechanisms under the HDV condition is crucial for rationally designing highly active and durable platinum group metal (PGM) catalysts for high-performance membrane electrode assemblies (MEAs). Herein, we report a PGM catalyst consisting of platinum nanoparticles with a high content (40 wt %) on atomic-metal-site (e.g., MnN 4 )-rich carbon support. MEAs with the Pt (40 wt %)/Mn–N–C cathode catalyst achieved significantly enhanced performance and durability, generating 1.41 A cm –2 at 0.7 V under HDV conditions (0.25 mgPt cm –2 and 250 kPa abs pressure) and retaining 1.20 A cm –2 after an extended and accelerated stress test up to 150,000 voltage cycles. Electron microscopy studies indicate that most fine Pt nanoparticles are retained on or/and in the carbon support covered with the ionomer throughout the catalyst layer at the end of life. During the long-term stability test, the observed electrochemical active surface area reduction and performance loss primarily result from Pt depletion in the catalyst layer due to Pt dissolution and redeposition at the interface of the cathode and membrane. Importantly, the first-principle density functional theory calculations further reveal a support entrapment effect of the Mn–N–C, in which the MnN 4 site can specifically adsorb the Pt atom and further retard the Pt dissolution and migration, therefore enhancing long-term MEA durability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Interfacial Reactivity-Triggered Oscillatory Lattice Strains of Nanoalloys

Understanding the structure evolution of nanoalloys under reaction conditions is vital to the design of active and durable catalysts. Herein, we report an operando measurement of the dynamic lattice strains of dual-noble-metal alloyed with an earth-abundant metal as a model electrocatalyst in a working proton-exchange membrane fuel cell using synchrotron high-energy X-ray diffraction coupled with pair distribution function analysis. Additionally, the results reveal an interfacial reaction-triggered oscillatory lattice strain in the alloy nanoparticles upon surface dealloying. Analysis of the lattice strains with an apparent oscillatory irregularity in terms of frequency and amplitude using time-frequency domain transformation and theoretical calculation reveals its origin from a metal atom vacancy diffusion pathway to facilitate realloying upon dealloying. This process, coupled with surface metal partial oxidation, constitutes a key factor for the nanoalloy's durability under the electrocatalytic oxygen reduction reaction condition, which serves as a new guiding principle for engineering durable or self-healable electrocatalysts for sustainable fuel cell energy conversion.

Wu, Zhi-Peng [State University of New York (SUNY),↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Agent-Based Coordination Scheme for PV Integration (ABC4PV)

Renewables and especially photovoltaics (PV) have benefitted significantly from a host of incentives and policies targeted toward enhanced integration and adoption of specific energy technologies. However, with the push to move forward into a subsidy-free market framework, behind-the-meter residential PV applications have generally struggled to retain their value (unlike utility scale and commercial projects) [1]. This project focused on developing control-theoretic solutions aimed at improving the integration and interaction of behind-the-meter residential PV with other distribution system assets (controllable and non-controllable) to enhance the integrated value of residential PV. To this end, a suite of decentralized control methodologies have been developed to enable effective coordination and control of behind-the-meter residential load customers’ PV, battery storage systems (BSS), controllable loads and other similar assets within a distribution feeder. This interaction aims at procuring energy savings and, thus, energy bill savings. The main source of savings is drawn from reducing the effect of demand charge pricing and is realized at the feeder level, assuming community level interaction and management among the aforementioned assets. Optimal control of the assets is implemented with a distributed optimization methodology, leveraging consensus-based algorithms. The results gathered from the optimal control simulations demonstrates that the savings can be duly achieved and the algorithm decision times (to dynamically control asset set points, for example) are fast. As for the overall efficiency of PV+BSS systems, to procure energy savings from curtailment of the demand charge pricing effects, the optimal control is set up so as to minimize the variance of the load for all customers, throughout a feeder and throughout time in a rolling horizon scheduling with model predictive control. The control takes into account inter-temporal electrochemical storage (battery) degradation costs: specifically, we have developed a long-term lifetime model for the BSS that weighs in the effect of the degradation factor in the dispatch formulations, thus, a considerable operating cost that affects energy decision making. The levelized cost of energy (LCOE – redefined for the purpose of quantifying asset integration effectiveness through the customers’ energy cost) is shown to be below the threshold set for the combined PV+BSS topology of $ 0.14/kWh for multiple cases of PV penetration all the way up to 50%, provided that a policy of shared ownership of and savings is in place. Further, the LCOE calculated for the case before the deployment PV+BSS systems is also achievable, i.e. the deployment of PV+BSS, if planned and scheduled optimally. will have no effect on customers’ energy costs. From the control methodology viewpoint, the developed consensus-based algorithms are shown to converge for a wide range of problem cases (spanning normal operating scenarios and contingencies), guaranteeing dispatch solutions under forecasting errors, communication break-downs and cyber-security attacks. The proposed control solutions are scalable and real-time implementable, with dispatch computations and device set-point updates converging in less than 2s in most practical instances of the above events.

14 SOLAR ENERGY↗
DOE OSTI2024

Toward a Machine Learning Approach to Interpreting X-ray Spectra of Trace Impurities by Converting XANES to EXAFS

The fact that the photoabsorption spectrum of a material contains information about the atomic structure, commonly understood in terms of multiple scattering theory, is the basis of the popular extended X-ray absorption spectroscopy (EXAFS) technique. How much of the same structural information is present in other complementary spectroscopic signals is not obvious. Here we use a machine learning approach to demonstrate that within theoretical models that accurately predict the EXAFS signal, the extended near-edge region does indeed contain the EXAFS-accessible structural information. We do this by exhibiting deep operator neural networks (DeepONets) that have learned the relationship between the extended and near edge portions of the X-ray absorption spectrum to predict the former from the latter. We find that we can accurately predict the EXAFS spectrum between 6 and 14 Å –1 from the first 6 Å –1 (≈100 eV) of the absorption spectrum of Cu 2 + substitutional defects in the Fe 3+ mineral hematite (α-Fe 2 O 3 ). This surprising finding implies that theoretical analyses of X-ray absorption spectra could be implemented that extract the same conclusions as high-quality EXAFS studies from spectra collected over a much smaller range of photon energies. This relaxes a host of experimental limitations related to the X-ray source and measurement sample, including collection time, minimum dopant concentration, source brilliance, and energy range. We describe the theoretical data sets and DeepONet construction and show that the resulting DeepONets produce EXAFS that recovers linear combination fits to experimental data with accuracy approaching the original ab initio calculations. We discuss the implications of our findings for minor constituent characterization and for understanding the information content of spectroscopic data more broadly, including how this approach might be applied to measured experimental spectra. In conclusion, to encourage similar efforts, the simulated X-ray spectra, machine learning, and fitting code are publicly available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Model and Standard Operating Procedures Supporting Signal Variation Flow Graph Analysis

This document will describe the principles of the Signal Variation Flow Technique, and the uncertainty models generated using it. We focus on the capture of the variation between the ideal signal and the measured signal. The ideal signal is defined to represent the signal output of a system whose full state behavior is known, with no variation in environment or during operation, and whose photons trajectories are not modified by the object. A CT uncertainty analysis is the result of two main steps. First, the radiography regime extends from the source to the collected image, which is 2D in the case of standard CT. This maps all upstream uncertainties into the variation observed on the radiograph and captures all variation in the physical domain. Second, the reconstruction regime extends from the captured images to the reconstructed 3D image. This regime is purely in the mathematical domain and corrects reconstruction algorithm artifacts/anomalies. The present work focuses on building the model through the radiography regime. The reconstruction regime is expected to be largely a study of algorithmic sensitivity, requiring the definition of a range of standardized tests through which the algorithms would be run. Radiographic variation would then be mapped through reconstruction sensitivities to predict the signal variation in the final image. An incomplete list for the reconstructed image variation output basis includes voxel density, edge blur and length variation, in analogous fashion to the basis functions presented in this work. The signal variation occurs in several forms at the radiograph. These forms are gathered into a complete basis set of functions describing all variation on the radiograph. The scale of each basis function is calculated independently via a specific SVFG. This includes 0D pixel noise (0DI), 0D energy noise (0DE), 1D blur (1DB), 1D length (1DL) or 2D position (2DP). These models are orthogonal in that they each explore a space in the signal variation domain that cannot be reached by the other basis functions. The variation basis functions, and associated SVFG models are split by output dimensionality, a term loosely used for categorization, and explained further below.

42 ENGINEERING↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2022

Production of Fischer-Tropsch Synfuels at Nuclear Plants

A case study analysis was performed to evaluate nuclear-powered synthetic fuel production in the midwestern United States (U.S.). A Fischer-Tropsch (FT) fuel synthesis plant design was used as the basis for the analysis. The FT plant design was configured to produce a product slate consisting of diesel fuel, jet fuel, and motor gasoline blend stocks from carbon dioxide (CO 2 ) and hydrogen (H 2 ) feedstocks. The CO 2 feedstock for the FT plant was assumed to be sourced from biorefineries in the region around a Midwest light water reactor (LWR) nuclear power plant (NPP). The analysis specifies that power from the LWR is used to produce H 2 via high-temperature steam electrolysis and to operate the FT synfuel production plant. Capital costs were estimated for the FT plant while capital costs for the electrolysis plant were based on previous Idaho National Laboratory (INL) studies. In addition to labor and maintenance costs for the FT and electrolysis plants, operating costs also include the costs for CO 2 feedstock transport. An analysis was performed to determine the cost of transporting CO 2 from the distributed biorefinery sources to the centralized fuel synthesis plant as a function of the synfuel plant capacity and corresponding CO 2 demand. The primary revenue streams are associated with sales of the synthetic fuel products. The synthetic fuel products will likely follow the same market trends as the conventional fuel products. The synfuel price data was thus based on projections made by the U.S. Energy Information Administration (EIA) 2021 Annual Energy Outlook (AEO) for conventional fuel products minus federal and state taxes, as well as marketing and distribution costs. The economic analysis also considered cases that included and excluded revenues from the 2022 Inflation Reduction Act (IRA) clean hydrogen production tax credit (PTC) of $\$ $3.00/kg for the first ten years of operation. The economic analysis calculated the net present value (NPV) for cases involving steady-state synfuel production for comparison with the NPV for a business-as-usual case in which NPP continues to sell only electric power to the grid. A synfuel production “Reference Case” was considered in addition to sensitivity cases in which the plant capacity, electricity price, and synthetic fuel product prices were perturbed. The synfuel production Reference Case considered a scenario in which the electrolysis and synfuel plants utilized a combined electrical load of 1000 megawatt electrical (MWe) from the LWR with the balance of the LWR power output being sold to the electric grid. The economic analysis suggests that the synfuel production Reference Case evaluated in this analysis would lead to considerable economic potential for near-term deployment of a nuclear-based synfuel production plant. Specifically, the economic analysis suggests that the deployment of a 1000 megawatt (MW) nuclear-powered synfuel plant could result in a NPV increase of approximately $\$ $1.7 billion for a case with no clean synfuel price premium relative to conventional petroleum fuels when accounting for the additional revenues from the 2022 IRA clean hydrogen PTCs of $\$ $3/kg. Sensitivity analysis was performed to evaluate the effect of perturbation of selected model input parameters on the NPV for the synfuel production Reference Case. The sensitivity analysis indicates that the plant capacity has the largest impact on the differential NPV, with a smaller synfuel production capacity resulting in a decrease in revenue when a larger fraction of the power from the NPP is sold to the grid and a smaller fraction of the power is used to produce synthetic fuel products. The synfuel product pricing has the next largest impact on the differential NPV, with lower synfuel prices resulting in decreased NPV from decreased synfuel sales revenue while higher synfuel prices result in increased NPV from increased synfuel sales revenue. Electricity pricing has a smaller effect on the NPV than the fuel sales price since, in the Reference Case, most of the energy from the NPP is used for synfuel production and a smaller amount of the system revenues are associated with electrical power sales. However, the electricity price sensitivity does indicate that the Synfuel Integrated Energy System (IES) would have a greater NPV than the business-as-usual case (e.g., grid power sales only) when electricity market prices are low, suggesting that synfuel production could provide a strategy for decreasing the economic risks to NPPs posed by a loss of revenues attributed to falling electricity market prices.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Interplay Between Metastability and Mechanically Induced Structural Instability in CsPbBr 3 Photovoltaic Perovskite

An atomic-level understanding of the underlying structural metastability is still absent in halide perovskite photovoltaic systems. Focusing on the model material CsPbBr 3 , the impact of mechanically induced atomic structure alteration is elucidated through structural modeling and X-ray diffraction measurements. Sudden Cs–Br bond breaking drives the system metastability, where the first-order transition arises from cation-halide bond strain relief by severing the corner connectivity of the PbBr 6 units. The pressure–volume and entropy terms govern the Gibbs free-energy landscape in halide perovskites. In conclusion, metastability is revealed as a critical factor limiting the performance of lead halide perovskites under extreme but natural environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Near-Term Effects of Perennial Grasses on Soil Carbon and Nitrogen in Eastern Nebraska

Incorporating native perennial grasses adjacent to annual row crop systems managed on marginal lands can increase system resiliency by diversifying food and energy production. This study evaluated (1) soil organic C (SOC) and total N stocks (TN) under warm-season grass (WSG) monocultures and a low diversity mixture compared to an adjacent no-till continuous-corn system, and (2) WSG total above-ground biomass (AGB) in response to two levels of N fertilization from 2012 to 2017 in eastern Nebraska, USA. The WSG treatments consisted of (1) switchgrass (SWG), (2) big bluestem (BGB), and (3) low-diversity grass mixture (LDM; big bluestem, Indiangrass, and sideoat grama). Soils were sampled at fixed depth increments (0–120 cm) in the WSG plots and in the adjacent corn experiment in 2012 and 2017. Soil stocks (Mg ha−1) of SOC and TN were calculated on an equivalent soil mass (ESM) basis and compared within the three WSG treatments as well as between experiments (corn compared to the mean of all WSGs). Soil organic C and TN stocks within soil layers and cumulative stocks responded to the main effect of WSG (PWSG < 0.05) but were no different when comparing the WSGs to corn (Pexpt = NS). Both SOC/TN stocks and cumulative stocks were generally greater in the LDM compared to the BGB. Neither SOC nor TN changed over time under either the WSGs or corn. Warm-season grass AGB responded to a three-way interaction of year, N rate, and WSG (p = 0.0007). Decreases in AGB over time were significant across WSGs and N levels except for SWG at 56 kg N ha−1 and LDM at 112 kg N ha−1. Above-ground biomass was generally greater in the LDM after the first harvest year (2013). Results suggest that incorporating WSGs into marginal cropland can maintain SOC and TN stocks while providing a significant source of biomass to be used in energy production or in integrated livestock systems.

Ramirez, II, Salvador (ORCID:000000030152153X)↗
DOE OSTI2022

Experimental study of the Si 30 ( He 3 , d ) P 31 reaction and thermonuclear reaction rate of Si 30 ( p , γ ) P 31

Background: Abundance anomalies in some globular clusters, such as the enhancement of potassium and the depletion of magnesium, can be explained in terms of an earlier generation of stars polluting the presently observed ones. It was shown that the potential range of temperatures and densities of the polluting sites depends on the strength of a few critical reaction rates. The 30 Si(p,γ) 31 P reaction has been identified as one of these important reactions. Purpose: The key ingredient for evaluating the thermonuclear 30 Si(p,γ) 31 P reaction rate is the strength of the resonances which, at low energy, are proportional to their proton width. Therefore, the goal of this work is to determine the proton widths of unbound 31 P states. Method: States in 31 P were studied at the Maier-Leibnitz-Laboratorium using the one-proton 30 Si( 3 He,d) 31 P transfer reaction. Deuterons were detected with the high resolution one quadrupole three dipoles (Q3D) magnetic spectrometer. Angular distribution and spectroscopic factors were extracted for 27 states, and proton widths and resonance strengths were calculated for the unbound states. Results: Several 31 P unbound states have been observed for the first time in a one-proton transfer reaction. Above 20 MK, the 30 Si(p,γ) 31 P reaction rate is now entirely estimated from the observed properties of 31 P states. The reaction rate uncertainty from all resonances other than the E$^{c.m.}_{r}$ = 149 keV resonance has been reduced down to less than a factor of 2 above that temperature. The unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance dominates the uncertainty in the rate in the relevant temperature range. Conclusion: Here, the remaining source of uncertainty on the 30 Si(p,γ) 31 P reaction rate comes from the unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance which can change the reaction rate by a factor of 10 in the temperature range of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Electronic and magnetic properties of hole-doped topological kagome Fe 1−𝑥 ⁢Mn 𝑥 ⁢Sn thin films

We have investigated the electronic and magnetic structures of topological kagome Fe 1−𝑥 ⁢Mn 𝑥 ⁢Sn (0 ≤ 𝑥 ≤ 0.3) thin films via neutron diffraction, electronic transport measurements, and ab initio density functional theory (DFT) to understand the interplay between hole doping, magnetism, and the electronic structures. Temperature-dependent neutron diffraction measurements on parent FeSn reveal the Néel temperature to be 𝑇 N ∼ 355 K and the underlying A-type antiferromagnetic ordering is associated with a wave vector 𝒒 = (001/2). Upon Mn doping to 𝑥 = 0.15, 𝑇 N decreases slightly while the magnetic ordering vector remains the same. Resistivity measurements show metallic characteristics and in-plane anisotropy down to 10 K for all the investigated samples. The effects of hole doping are mapped in terms of electronic ground state calculations via DFT which show that the Dirac point is moved closer to the Fermi level (𝐸 F ) and the flat bands get pushed away from 𝐸 F upon hole doping. However, a comparison between hole-doped Fe 1−𝑥⁢ Mn 𝑥⁢ Sn and electron-doped Fe 1−𝑥 ⁢Co 𝑥 ⁢Sn indicates that the Néel temperature does not scale with the position of 𝐸 F relative to the flat band. Furthermore, our results establish the antiferromagnetic state of FeSn and Fe 1−𝑥 ⁢Mn 𝑥⁢ Sn films at room temperature, laying the groundwork for future studies of magnetism in kagome heterostructures.

36 MATERIALS SCIENCE↗