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 343 records · Page 19

Assessment of Existing Transportation Packages for Use with HALEU

Commercial light water reactor operators and fuel vendors in the United States are pursuing changes to fuel that include increased 235 U enrichment. Economic studies generally anticipate maximum near-term fuel assembly designs with up to 8 wt.% 235 U. Many next-generation nuclear reactor designs require high-assay low-enriched uranium (HALEU) (19.75 wt.% > 235 U > 5 wt.%) fuel. One necessary element for the commercial-scale use of HALEU is the ability to safely transport large quantities of enriched fuel material in multiple forms. However, there is uncertainty as to whether subcriticality requirements can be satisfied with existing package designs and whether existing critical benchmark experiment data are sufficient to support criticality safety code validation for HALEU transportation applications. This study assesses the potential to use currently licensed transportation packages for the transportation of increased enrichment unirradiated U fuel forms. The assessment uses selected package designs that represent the five categories of fuel form-boiling water reactor pins and assemblies, pressurized water reactor pins and assemblies, UF 6 , U-metal and tristructural isotropic (TRISO) particles, and UO 2 pellets or powder-and focuses on demonstrating subcriticality and identifying benchmark critical experiments appropriate for use in criticality computer code validation. Key quantities of interest that relate to subcriticality are limiting conditions (e.g., optimum moderation), package or package array k eff , package capacity, and package transportation array size. The SCALE TSUNAMI-3D and TSUNAMI Indices and Parameters codes are used for sensitivity and uncertainty calculations and for identification of candidate critical benchmark experiments for code validation. The primary metric for identifying candidate benchmarks is the similarity coefficient, c k . For each fuel form category, a representative package is evaluated. Results provided for each package evaluation include enrichment and packaging limits (e.g., maximum transportation array size as a function of enrichment) and benchmark critical experiment similarity coefficients. Results indicate that there are viable means for increasing enrichments into the HALEU range across the spectrum of fuel forms with differing increase amounts available for different packages. Sources of subcriticality margin to offset increased enrichment reactivity include reduced transportation array size, reduced fissile mass, burnable absorber credit, and safety analysis margin harvesting. For all packages except the DN-30, numerous critical benchmark experiment candidates for validation were identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental and modeling studies of PR and ND oxalate solubility to high ionic strengths: Insight into actinide(III) oxalates

Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a series of various files organized by state sub-directories. The first two characters of each file name is the abbreviation for the state the in which the file contains data for. After the abbreviation is the text which describes the contents of the file. Here is each file type described in detail: XXFullHarmonizedRights.csv: A file of the combined groundwater and surface water records for each state. Essentially, this file is the merging of XXGroundwaterHarmonizedRights.csv and XXSurfaceWaterHarmonizedRights.csv by state. The column headers for each of this type of file are: state - The name of the state the data comes from. FIPS - The two-digit numeric state ID code. waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. priorityDate - The priority date associated with the right. origWaterUse - The original stated water use(s) from the state. waterUse - The water use category under the unified use categories established here. source - Whether the right is for surface water or groundwater. basinNum - The alpha-numeric identifier of the WMA the record belongs to. CFS - The maximum flow of the allocation in cubic feet per second (ft3s-1). Arizona is unique among the states, as its surface and groundwater resources are managed with two different sets of boundaries. So, for Arizona, the basinNum column is missing and instead there are two columns: surBasinNum - The alpha-numeric identifier of the surface water WMA the record belongs to. grdBasinNum - The alpha-numeric identifier of the groundwater WMA the record belongs to. XXStatePOD.shp: A shapefile which identifies the location of the Points of Diversion for the state's water rights. It should be noted that not all water right records in XXFullHarmonizedRights.csv have coordinates, and therefore may be missing from this file. XXStatePOU.shp: A shapefile which contains the area(s) in which each water right is claimed to be used. Currently, only Idaho and Washington provided valid data to include within this file. XXGroundwaterHarmonizedRights.csv: A file which contains only harmonized groundwater rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. XXSurfaceWaterHarmonizedRights.csv: A file which contains only harmonized surface water rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. Additionally, one file, stateWMALabels.csv, is not stored within a sub-directory. While we have referred to the spatial boundaries that each state uses to manage its water resources as WMAs, this term is not shared across all states. This file lists the proper name for each boundary set, by state. For those whom may be interested in exploring our code more in depth, we are also making available an internal data file for convenience. The file is in .RData format and contains everything described above as well as some minor additional objects used within the code calculating the cumulative curves. For completeness, here is a detailed description of the various objects which can be found within the .RData file: states: A character vector containing the state names for those states in which data was collected for. More importantly, the index of the state name is also the index in which that state's data can be found in the various following list objects. For example, if California is the third index in this object, the data for California will also be in the third index for each accompanying list. rightsByState_ground: A list of data frames with the cleaned ground water rights collected from each state. This object holds the the data that is exported to created the xxGroundwaterHarmonizedRights.csv files. rightsByState_surface: A list of data frames with the cleaned surface water rights collected from each state. This object holds the the data that is exported to created the xxSurfaceWaterHarmonizedRights.csv files. fullRightsRecs: A list of the combined groundwater and surface water records for each state. This object holds the the data that is exported to created the xxFullHarmonizedRights.csv files. projProj: The spatial projection used for map creation in the beginning of the project. Specifically, the World Geodetic System (WGS84) as a coordinate reference system (CRS) string in PROJ.4 format. wmaStateLabel: The name and/or abbreviation for what each state legally calls their WMAs. h2oUseByState: A list of spatial polygon data frames which contain the area(s) in which each water right is claimed to be used. It should be noted that not all water right records have a listed area(s) of use in this object. Currently, only Idaho and Washington provided valid data to be included in this object. h2oDivByState: A list of spatial points data frames which identifies the location of the Point of Diversion for the state's water rights. It should be noted that not all water right records have a listed Point of Diversion in this object. spatialWMAByState: A list of spatial polygon data frames which contain the spatial WMA boundaries for each state. The only data contained within the table are identifiers for each polygon. It is worth reiterating that Arizona is the only state in which the surface and groundwater WMA boundaries are not the same. wmaIDByState: A list which contains the unique ID values of the WMAs for each state. plottingDim: A character vector used to inform mapping functions for internal map making. Each state is classified as either "tall" or "wide", to maximize space on a typical 8x11 page. The code related to the creation of this dataset can be viewed within HarDWR GitHub Repository/dataHarmonization.

Economics↗

Recommendations to Improve the Nuclear Regulatory Commission Reactor Licensing and Approval Process

Due to the urgency around climate change and associated goals for clean energy transition, and Battelle Energy Alliance's (BEA) and Idaho National Laboratory's (INL) role as the nation's nuclear energy laboratory, numerous stakeholders have asked for BEA’s thoughts and recommendations to reduce the time and costs associated with licensing new reactors at the U.S. Nuclear Regulatory Commission (NRC). As an M&O contractor for an FFRDC, BEA is a long-term partner with the Government in seeking to achieve clean energy goals, yet has a level of independence needed to appropriately evaluate this topic. The views herein are informed by extensive BEA experience supporting nuclear energy endeavors including ongoing discussions with current and former regulators, applicants, and licensees. With this background in mind, the United States benefits from having an agency such as the NRC, which is viewed internationally as the leader in nuclear safety licensing and regulation. Nonetheless, while acknowledging the important nuclear safety role satisfied by the NRC, it is apparent that one of the most significant time and resource intensive activities for new reactor developers is the NRC licensing process. The time and cost to obtain NRC licenses add significant financial stress for new reactor projects and may result in abandonment of projects or failure to even begin new projects. The challenge is particularly acute for advanced reactors which may raise unique or new regulatory questions and may be smaller in size, resulting in a much higher proportional impact from regulatory and cost challenges. This situation presents a particularly troublesome risk for the nation given the urgency in which utilities are working to transition to clean, non-carbon-emitting energy sources like nuclear energy. Reforms to the NRC licensing process have the potential to greatly increase certainty and support the successful progress of new reactors. The NRC can retain its world-class nuclear safety reputation while becoming a world leader for regulatory efficiency and a critical enabler to the clean energy transition. This report describes potential NRC reforms, focusing on those with a statutory connection. Recognizing the potential tradeoffs with any proposed changes, the report attempts to highlight those considerations in the analysis of the reforms. The recommendations are presented as a set of options for consideration. Unless noted, they are independent options, offering stakeholders the option to select a subset for further consideration. Although difficult to calculate precise time improvements for some of the changes, the reforms have the potential for substantial improvements, perhaps even by a factor of two.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Modeling and Analysis of the Transport and Disposal of Beryllium Moderator Blocks and Greater than Class C (GTCC) Waste

Radioactive waste in the United states is categorized based on its source, radioactivity, and security risk. The categorization method employed by the Nuclear Regulatory Committee (NRC) for low-level radioactive waste is not currently applied to any waste generated by the Department of Energy (DOE) or disposed of in DOE facilities. As a result, there is a large volume of DOE-generated waste with characteristics similar to the NRC’s “Greater than Class C” (GTCC) waste category which presently do not have a path for long-term disposal. Much of this waste cannot undergo the standard commercial disposition process due to it being categorized as Transuranic (TRU) waste under DOE guidelines, due to elevated concentrations of key fission products. Possible solutions to this dilemma are explored in the Environmental Impact Statement for the “Disposal of Greater-Than-Class-C Low Level Radioactive Wave” (EIS-0375). This paper analyzes several EIS possible paths for disposal for this “orphaned” waste. One example that highlights this categorization issue is the spent beryllium cladding that has been extracted from the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) throughout its operational lifetime. These beryllium blocks surround the reactor’s main chambers and components and serve as neutron moderators. This paper will use this particular waste form to analyze the economic and technological viability of transporting and storing this material under the storage and transportation criteria of the Waste isolation Pilot Plant (WIPP), which was deemed the most feasible disposition path according to the EIS. The example of this waste form used in this paper that highlights this categorization issue is the spent beryllium cladding extracted from Idaho National Laboratory’s (INL) Advanced Test Reactor (ATR). These beryllium blocks surround the reactor’s components and serve as neutron moderators, and thus have accumulated high concentration of high-activity transuranic isotopes. This paper analyzed the viability of transporting and storing this material under the storage and transportation criteria of the five primary disposition paths covered in EIS-0375. An array of calculations was carried out to assess the viability of the various disposition paths for the beryllium shipments as well as other waste shipments that fall within the GTCC category. Modeling with MCNP 6.2 was conducted to determine whether the beryllium blocks could be safely stored and transported within a 72-B cask; the standard shipping container for remote-handled, transuranic waste, while also meeting regulatory limits at various disposition paths. A cost analysis of the transportation and long-term disposal paths mentioned in the EIS was also carried out using available data and information from similar waste shipments. Lastly, a geochemical analysis of the various geological repository discussed in the EIS was also carried out using the Geochemists Workbench Release 14. Long-term disposal of the beryllium blocks at the Waste Isolation Pilot Plant (WIPP) proved to be the most cost-effective long-term disposition option of the ones considered in the EIS. A dose rate calculation at both contact and remote-handling distances indicate that the analyzed beryllium shipments should not exceed the exposure limits at any of the considered locations. Greater-than-class-C waste has been left in a regulatory state of limbo for years, resulting in backlogs of inventory across several research sites in the United States. Due to its high activity and presence of transuranic isotopes, it is imperative to ensure that it remains inaccessible and sequestered both in its short-term interim as well as a long-term geologic time scale. The information and assessments done in the paper could potentially serve as a reference for any future shipments of this waste form at the WIPP facility as well as other disposition paths that may be considered in the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Position-specific kinetic isotope effects for nitrous oxide: a new expansion of the Rayleigh model

Nitrous oxide (N 2 O) is a potent greenhouse gas and the most significant anthropogenic ozone-depleting substance currently being emitted. A major source of anthropogenic N 2 O emissions is the microbial conversion of fixed nitrogen species from fertilizers in agricultural soils. Thus, understanding the enzymatic mechanisms by which microbes produce N 2 O has environmental significance. Measurement of the 15 N/ 14 N isotope ratios of N 2 O produced by purified enzymes or axenic microbial cultures is a promising technique for studying N 2 O biosynthesis. Typically, N 2 O-producing enzymes combine nitrogen atoms from two identical substrate molecules (NO or NH 2 OH). Position-specific isotope analysis of the central (N α ) and outer (N β ) nitrogen atoms in N 2 O enables the determination of the individual kinetic isotope effects (KIEs) for N α and N β , providing mechanistic insight into the incorporation of each nitrogen atom. Previously, position-specific KIEs (and fractionation factors) were quantified using the Rayleigh distillation equation, i.e., via linear regression of δ 15 N α or δ 15 N β against [–f In f / (1 – f)], where f is the fraction of substrate remaining in a closed system. This approach, however, is inaccurate for N α and N β because it does not account for fractionation at N α affecting the isotopic composition of substrate available for incorporation into the β position (and vice versa). Therefore, we developed a new expansion of the Rayleigh model that includes specific terms for fractionation at the individual N 2 O nitrogen atoms. By applying this Expanded Rayleigh model to a variety of simulated N 2 O synthesis reactions with different combinations of normal, inverse, and/or no KIEs at N α and N β , we demonstrate that our new model is both accurate and robust. We also applied this new model to two previously published datasets describing N 2 O production from NH 2 OH oxidation in a methanotroph culture (Methylosinus trichosporium) and N 2 O production from NO by a purified Histoplasma capsulatum (fungal) P450 NOR, demonstrating that the Expanded Rayleigh model is a useful tool in calculating position-specific fractionation for N 2 O synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

iSPECTRON: a simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software

We introduce iSPECTRON, an open source (under the Educational Community License version 2.0) program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained through quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The code capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Two levels of electronic structure theory, TDDFT (with NWChem) and RASSCF/RASPT2 (with OpenMolcas), are compared where possible. Acknowledgements: F.S., A.N., D.R.N., N.G., S.M, M.G. acknowledge support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division under Award Nos. DE-SC0019484, KC-030103172684. The Spectron code was developed with support from the National Science Foundation (Grant CHE- 1953045). This research benefited from computational resources provided by EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at PNNL. PNNL is operated by Battelle Memorial Institute for the United States Department of Energy under DOE Contract No. DE-AC05-76RL1830.

Segatta, Francesco↗

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗

Dosimetric and biological impact of activity extravasation of radiopharmaceuticals in PET imaging

The increasing use of nuclear medicine and PET imaging has intensified scrutiny of radiotracer extravasation. To our knowledge, this topic is understudied but holds great potential for enhancing our understanding of extravasation in clinical PET imaging. This work aims to (1) quantify the absorbed doses from radiotracer extravasation in PET imaging, both locally at the site of extravasation and with the extravasation location as a source of exposure to bodily organs and (2) assess the biological ramifications within the injection site at the cellular level. A radiation dosimetry simulation was performed using a whole-body 4D Extended Cardiac-Torso (XCAT) phantom embedded in the GATE Monte Carlo platform. A 10-mCi dose of 18 F-FDG was chosen to simulate a typical clinical PET scan scenario, with 10% of the activity extravasated in the antecubital fossa of the right arm of the phantom. The extravasation volume was modeled as a 5.5 mL rectangle in the hypodermal layer of skin. Absorbed dose contributions were calculated for the first two half-lives, assuming biological clearance thereafter. Dose calculations were performed as absorbed doses at the organ and skin levels. Energy deposition was simulated both at the local extravasation site and in multiple organs of interest and converted to absorbed doses based on their respective masses. Each simulation was repeated ten times to estimate Monte Carlo uncertainties. Biological impacts on cells within the extravasated volume were evaluated by randomizing cells and exposing them to a uniform radiation source of 18 F and 68 Ga. Particle types, their energies, and direction cosines were recorded in phase space files using a separate Geant4 simulation to characterize their entry into the nucleus of the cellular volume. Subsequently, the phase space files were imported into the TOPAS-nBio simulation to assess the extent of DNA damage, including double-strand breaks (DSBs) and single-strand breaks (SSBs). Organ-level dosimetric estimations are presented for 18 F and 68 Ga radionuclides in various organs of interest. With 10% extravasation, the hypodermal layer of the skin received the highest absorbed dose of 1.32 ± 0.01 Gy for 18 F and 0.99 ± 0.01 Gy for 68 Ga. The epidermal and dermal layers received absorbed doses of 0.07 ± 0.01 Gy and 0.13 ± 0.01 Gy for 18 F, and 0.14 ± 0.01 Gy and 0.29 ± 0.01 Gy for 68 Ga, respectively. In the extravasated volume, 18 F caused an average absorbed dose per nucleus of 0.17 ± 0.01 Gy, estimated to result in 10.58 ± 0.50 DSBs and 268.11 ± 12.43 SSBs per nucleus. For 68 Ga, the absorbed dose per nucleus was 0.11 ± 0.01 Gy, leading to an estimated 6.49 ± 0.34 DSBs and 161.24 ± 8.12 SSBs per nucleus. Absorbed doses in other organs were on the order of micro-gray (µGy). The likelihood of epidermal erythema resulting from extravasation during PET imaging is low, as the simulated absorbed doses to the epidermis remain below the thresholds that trigger such effects. Moreover, the organ-level absorbed doses were found to be clinically insignificant across various simulated organs. The minimal DNA damage at the extravasation site suggests that long-term harm, such as radiation-induced carcinogenesis, is highly unlikely.

DNA strand breaks↗

California Price Response Potential Study

California's energy landscape is undergoing a significant transformation, driven by the increasing integration of renewable energy sources, the increased adoption of distributed energy resources, the electrification of end-use loads, and the growing need for grid efficiency. To address these challenges, recent revisions to the State’s Load Management Standards (LMS) require all of California’s large utilities and community choice aggregators (CCAs) to offer dynamic electricity pricing options to customers by 2027. Dynamic pricing, which involves varying electricity rates based on real-time supply and demand conditions, offers a promising solution for optimizing grid operations, reducing costs, and incentivizing efficient use of grid capacity. Effective implementation of dynamic pricing requires understanding the potential impacts on customer bills, system load, and the cost-effectiveness of automation technologies. This study aims to evaluate the load response of various end-use devices to hourly dynamic prices. The end-uses studied here are space cooling, space heating, water heating, crop irrigation, pool and spa pumps, and electric vehicle (EV) charging, all for both residential and commercial applications, except for crop irrigation. In 2030, these end uses are forecasted to account for 18% of annual electricity demand in the state, but 40% of demand in the peak net load hour. By modeling possible price-responsive load dispatch algorithms and assessing the resulting impacts on both individual bills and the overall grid, we seek to inform policymakers and utilities about the potential benefits and challenges associated with dynamic pricing, and considerations for the design of dynamic pricing tariffs. Additionally, we will explore the cost effectiveness of adopting automation technologies to enable devices to respond more effectively to real-time price signals. This study considers a range of price profiles, accounting for differences across utilities and customer classes, and presents scenarios for dynamic price design via variation in the percentage of total customer electric costs that are allocated dynamically (versus constituting a fixed portion of the hourly volumetric price). We present results focused primarily on 2030, forecasting electricity prices under both low and high-cost scenarios, to inform longer-term tariff design considerations. We design tariffs by starting with 2019 prices that were calculated according to CalFUSE guidance (CPUC, 2022) and that have been used in recent studies; these prices are all-in volumetric rates that vary by utility and are revenue-neutral to each customer class. They are developed by considering six electricity cost components that are allocated hourly based on system load indicators (gross and net load, and wholesale prices). These prices are forecasted to 2030 for low and high cost scenarios, considering recent trends in total electricity costs with and without years of substantial wildfire mitigation investments. These tariffs, which allocate all costs on an hourly basis, are considered our “Full” dynamic tariff design scenario, while two additional scenarios explore allocating a portion of costs as a flat volumetric charge: the “Medium” scenario allocates 50% of revenue dynamically (and keeps 50% flat), while the “Mild” scenario allocates 20% of revenue dynamically. The 20% dynamic allocation on the Mild scenario aims to represent a case where only the marginal operating costs of the grid are included in the dynamic price.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Final Report for ARPA-E LOCOMOTIVES Advanced Locomotive Technology and Rail Infrastructure Optimization System (ALTRIOS) Project

The Advanced Locomotive Technology and Rail Infrastructure Optimization System (ALTRIOS) is a unique, fully integrated, open-source software tool used to evaluate strategies for cost-effectively deploying advanced locomotive technologies and associated infrastructure. ALTRIOS simulates freight-demand-driven train scheduling, mainline meet-pass planning, locomotive dynamics, train dynamics, energy conversion efficiencies, and energy storage dynamics of line-haul train operations. Because new locomotives represent a significant long-term capital investment and new technologies must be thoroughly demonstrated before deployment, this tool provides guidance on the risk/reward trade-offs and operation integration of different technology rollout strategies. An open, integrated simulation tool is valuable for identifying future research needs and making decisions on technology development, routes, and train selection. This final report details the ALTRIOS software architecture, major modules and components, and data validation process. It demonstrates the software's utility through a 30-year rollout case study targeting high penetration of advanced powertrain technologies by 2050 for two BNSF Railway routes: loaded taconite ore trains from Hibbing, Minnesota, to Superior, Wisconsin, and mixed-freight trains from Superior to Minneapolis, Minnesota.

33 ADVANCED PROPULSION SYSTEMS↗

Sneutrino tribrid inflation, metastable cosmic strings and gravitational waves

Here, we present a successful realization of sneutrino tribrid inflation model based on a gauged U(1)_B-L extension of Minimal Supersymmetric Standard Model (MSSM). A single interaction term involving the B-L Higgs field and the right-handed neutrinos serves multiple purposes. These include the generation of heavy Majorana masses for the right-handed neutrinos to provide an explanation for the tiny neutrino masses via the seesaw mechanism, a realistic scenario for reheating and non-thermal leptogenesis with a reheat temperature as low as 10 6 GeV, and a successful realization of inflation with right-handed sneutrino as the inflaton. The matter parity which helps avoid rapid proton decay survives as a Z 2 subgroup of a U(1) R-symmetry. Depending on the choice of model parameters yields the following predicted range of the tensor to scalar ratio, 3 x 10 -11 ≲ r≲ 7x 10 -4 (6 x 10 -7 ≲ r ≲ 0.01), and the running of the scalar spectral index, -0.00022 ≲ dn_s/dln k ≲ -0.0026 (-0.00014 ≲ dn_s/dln k ≲ 0.005), along with the B-L breaking scale, 3 x 10 14 ≲ M/ GeV≲ 5 x 10 15 (6 x 10 15 ≲ M/ GeV≲ 2 x 10 16 ), calculated at the central value of the scalar spectral index, n_s =0.966, reported by Planck 2018. The possibility of realizing metastable cosmic strings in a grand unified theory setup is briefly discussed. The metastable cosmic string network admits string tension values in the range 10 -8 ≲ Gμ s ≲ 10 -6 , and predicts a stochastic gravitational wave background lying within the 2-σ bounds of the recent NANOGrav 12.5-yr data.

supersymmetry and cosmology↗

Energy Northwest - Horn Rapids Solar and Storage: An Assessment of Battery Technical Performance

Chartered in 1957 as a joint action agency of the state, Energy Northwest (ENW) is a consortium of 27 public utility districts and municipalities across Washington state. ENW takes advantage of economies of scale and shared services to help utilities run their operations more efficiently and at lower cost, to the benefit of more than 1.5 million customers. ENW develops, owns, and operates a diverse mix of electricity generating resources, including hydro, solar, and wind projects – and the Northwest’s only active nuclear energy facility. These projects provide enough reliable, affordable, and environmentally responsible energy to power more than a million homes each year, and that carbon-free electricity is provided at the cost of generation. The agency continually explores new generation projects to meet its members’ needs. In 2017, as part of the second round of funding from the Washington state Clean Energy Fund, the Washington State Department of Commerce granted up to $3 million in matching funds to develop an estimated $6.5 million project that deployed a 4-MW, 20-acre solar generating array of photovoltaic (PV) panels coupled with a 1 MW/5.5 MWh lithium-iron-phosphate battery energy storage system (BESS) in Richland, Washington. The combination of PV and BESS will provide a predictable, renewable generating source and will also serve as a training ground for solar and battery technicians throughout the nation. The City of Richland will purchase the power from the project and utilize the benefits of the energy storage. The project provides Washington state with its first opportunity to integrate a large-scale solar and storage facility into its clean mix of hydro, nuclear, and wind resources. This first-of-its-kind facility combines solar generation with battery storage and technician training. In 2019, Pacific Northwest National Laboratory (PNNL) worked with ENW to assess the integrated PV and BESS in representative use cases that could benefit the City of Richland. Between March and May 2022, extensive testing was conducted, and the results were used to assess the technical performance of the BESS subjected to actual field operations. Both reference performance and use case tests were performed: (A) Reference performance tests assess the general technical capabilities of the BESS, such as energy capacity, round-trip efficiency (RTE), ramp rate, and signal tracking capability. These are the first tests performed (baseline), and they are repeated after use case tests (post cycle). A standardized U.S. Department of Energy (DOE) energy storage performance protocol was used to characterize the BESS, including representative duty cycle profiles, test procedure guidance, and calculation guidance for determining key characteristics. (B) Use case tests examine the performance of the BESS for specific use cases using duty cycles developed by PNNL in collaboration with ENW. Five use cases were selected for testing: 1) demand charge reduction, 2) load shaping, 3) transmission charge reduction, 4) Volt-VAR service, and 5) outage mitigation. The use case duty cycles were developed based on utility and site-specific characteristics in addition to the technical characteristics and physical capabilities of the BESS. Use case tests were performed between the baseline and post cycle tests. This report describes the BESS and its components, presents testing and performance analysis results, and shares key insights and lessons learned from this project. Outcomes of the tests and analyses will help ENW understand the performance of the Horn Rapids BESS in its current state and design appropriate operational strategies for this and other BESSs over the long term.

14 SOLAR ENERGY↗

Potential State Regulatory Pathways to Facilitate Low-Carbon Fuels

States and the federal government are increasingly engaged in the challenges around decarbonizing the electric grid. In particular, regulators, consumers, stakeholders, and utilities recognize the need to carefully consider the role natural gas will play in a decarbonized future. A variety of technology and policy options to reduce greenhouse gas emissions associated with natural gas use are available, including energy efficiency programs, demand reduction tools, strategic electrification, and strategies to reduce emissions from natural gas production, transportation, and consumption. Low-carbon fuels – mainly renewable natural gas (RNG) and clean hydrogen – are being considered an important component of decarbonization goals. RNG and hydrogen may be able to meaningfully reduce emissions from processes independent of geologic natural gas, displacing emissions of methane, a powerful greenhouse gas. Although RNG and hydrogen are not cost-competitive today with geologic natural gas and are smaller in scale and potential than other decarbonization options, they can be explored as potential critical tools to decarbonize sectors that are difficult to electrify or shift off of natural gas entirely, such as air travel, industrial processes, maritime transport, long-distance trucking, space heating on cold days, and railroads (Nadel, 2022). The role of this report is to provide informational context for state utility regulators to understand the impacts of and challenges associated with broader integration of low-carbon fuels, followed by examples of state regulatory actions taken to date to facilitate the development of low-carbon fuels. Setting clear guidance to calculate the environmental benefits of low-carbon fuels and continuing federal and state investments in research and development to reduce costs relative to fossil fuels will be important steps to take to signal the desire to grow the market for these fuels. State public utility commissions may play a key role in setting regulatory frameworks for low-carbon fuels and ensuring that ratepayer funds, if utilized, are done so to further the public interest. This report is intended to summarize decisions that states have made to date on low-carbon fuels. In the spirit of understanding the current market and sharing information, this report provides success stories, and lessons learned across states as regulators implement varying strategies to achieve decarbonization objectives while maintaining their focus on affordability, safety, and reliability of the energy system. The report begins with an introduction of the role of natural gas in the U.S. economy (Section I) and background information on natural gas use, decarbonization, and low-carbon fuels (Section II). Next, the report describes the current market by discussing the scale of current production, emissions intensity, resource potential, and costs of low-carbon fuels compared to geologic natural gas (Section III). Following these sections, the report describes four strategies states have employed to facilitate low-carbon fuels: opening exploratory dockets, approving voluntary tariffs for customers, approving interconnection tariffs for producers, and considering portfolio-wide procurement targets (Section IV). This section lists states that have taken actions in each category, citing utility filings, commission decisions, stakeholder comments, and other relevant sources. Finally, the report concludes with suggested questions regulators may wish to consider regarding low-carbon fuels, in the interest of preparing to make decisions in the future (Section V). These questions include: Are there existing regulatory or technical barriers to voluntary purchases of low-carbon fuels? Can customers work with utilities to procure low-carbon fuels; are producers able to interconnect projects without significant barriers to entry? Should the infrastructure and/or commodity costs of low-carbon fuels be socialized among all ratepayers, or borne solely by the large commercial and industrial (C&I) customers currently driving the market? Should regulated natural gas and/or electric utilities own and operate low-carbon fuel production? How should regulators consider the unique decarbonization potential of low-carbon fuels, particularly for hard-to-abate sectors, in decision-making? Is additional direction or clarity from state policymakers needed? What no-regrets approaches can help facilitate both near-term RNG development and long-term development of hydrogen and other zero-carbon fuels? We collectively wish to express our gratitude to the U.S. Department of Energy, Office of Fossil Energy and Carbon Management, for supporting this report and other technical assistance resources for state regulators on natural gas topics. State regulators operate under a variety of policy environments, and states have vastly different types of energy resources, infrastructure, and customers. While there is no optimal regulatory, policy, or technological solution that will be successful in every state, state regulators can benefit by exchanging lessons learned with their peers across the country. We look forward to continued engagement with our fellow commissioners, commission staff, NARUC, the U.S. Department of Energy, and other stakeholders to develop sound regulation in the public interest.

03 NATURAL GAS↗

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Collaboration with NSTX in Calculations of Radiofrequency and Neutral Beam Heating and Current Drive

The overall objective of this project was to advance the understanding of heating and current drive necessary for toroidal magnetic confinement fusion energy reactors, and to provide accurate models for these sources to aid in understanding magneto-hydrodynamic, transport and boundary physics issues. Increased understanding is provided by computer simulation using the underlying physics laws, and will lead to greater reliability and increased efficiency of a power producing fusion reactor. The proposed work entailed applications of existing CompX Fokker-Planck, ray tracing, full-wave microwave, and particle diffusion calculator codes to the National Spherical Torus Experiment at Princeton Plasma Physics Laboratory. Special features of NSTX include that large ion orbits are emphasized in its geometry and thus offer a major opportunity within fusion research for focused exploration of their ubiquitous effects, including comparison with simulation codes. The CompX CQL3D Fokker-Planck code was be further adapted to the finite-orbit-width physics. Comparison with experiment ensured that the computer codes have sufficient accuracy to accurately model the experiments. Nonthermal distributions of energetic ions are produced in experiments-- and modeled with the codes-- exhibiting the signatures of auxiliary heating, transport, and radiation processes within the plasma. Using physics-based modeling, in conjunction with experimental data, nonthermal distributions can provide keys to validation and quantification of auxiliary heating physics. Nonthermal distributions also hold particular information on radial transport processes, not obscured by the smoothing Maxwellization effects of Coulomb collisions. Generally, radial transport is not well understood but is a key to efficient fusion energy production. By comparison with Fokker-Planck modeling, the velocity dependence of radial transport coefficients becomes experimentally accessible, particularly in cases of well-modeled sources such as neutral beams, fast wave power, or the toroidal electric field. It is clear from modeling of fast-wave interactions with ions that finite-orbit-width effects play a significant role. For example, a different set of ions is computed to intersect a given cyclotron resonance depending on whether or not radial orbit drifts are taken into account. Also, radial transport, associated bootstrap current, and alpha channeling depend on radial drifts. Ion cyclotron heating above the first harmonic interacts preferentially with higher velocity ions, giving large radial drifts. The CQL3D Fokker-Planck code, maintained by the proposers, is a physics based model used extensively within the US fusion community to compare with and interpret experimental observations. At the beginning of this grant, bounce-averages of the rf quasilinear and other operators have been performed neglecting the radial drifts which give finite orbit widths to the particle orbits. The proposed work included a major upgrade to CQL3D to include finite-orbit-width effects be accounted for by performing the bounce-averages of rf and neutral beam fast ions along guiding center orbits. This CQL3D enhancement provided microwave and neutral beam contributions to neoclassical radial transport coefficients for nonthermal distributions, and substantially increased the comprehensiveness of the ion model. Additional work validated Fokker-Planck and ray tracing modeling of neutral beam and high harmonic fast wave absorption and current drive against the NSTX experiment, as had been very successful for electron cyclotron current drive in the DIII-D tokamak at General Atomics and lower hybrid waves on C-Mod at MIT. The validation work was in conjunction with staff and graduate students at Princeton Plasma Physics Laboratory. This DOE grant supported fusion energy research, a potential long-term solution to the world's energy needs. The research grant under consideration, DE-SC0006614, provided research support at the 0.3 Full-Time-Equivalent level for application of computer models to aid in understanding and projecting efficacy of heating and current drive sources in the National Spherical Torus Experiment, a tokamak variant, at the Princeton Plasma Physics Laboratory. The NSTX experiment explores the physics of very tight aspect ratio, almost spherical tokamaks, aiming at producing steady-state fusion plasmas. The current drive is an integral part of the steady-state concept, maintaining the magnetic geometry in the steady-state tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗