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 811 records · Page 45

Solar imaging vector magnetograph

This report describes an instrument which has been constructed at the University of Hawaii to make observations of the magnetic field in solar active regions. Detailed knowledge of active region magnetic structures is crucial to understanding many solar phenomena, because the magnetic field both defines the morphology of structures seen in the solar atmosphere and is the apparent energy source for solar flares. The new vector magnetograph was conceived in response to a perceived discrepancy between the capabilities of X ray imaging telescopes to be operating during the current solar maximum and those of existing magnetographs. There were no space-based magnetographs planned for this period; the existing ground-based instruments variously suffered from lack of sensitivity, poor time resolution, inadequate spatial resolution or unreliable sites. Yet the studies of flares and their relationship to the solar corona planned for the 1991-1994 maximum absolutely required high quality vector magnetic field measurements. By 'vector' measurements we mean that the observation attempts to deduce the complete strength and direction of the field at the measurement site, rather than just the line of sight component as obtained by a traditional longitudinal magnetograph. Knowledge of the vector field permits one to calculate photospheric electric currents, which might play a part in heating the corona, and to calculate energy stored in coronal magnetic fields as the result of such currents. Information about the strength and direction of magnetic fields in the solar atmosphere can be obtained in a number of ways, but quantitative data is best obtained by observing Zeeman-effect polarization in solar spectral lines. The technique requires measuring the complete state of polarization at one or more wavelengths within a magnetically sensitive line of the solar spectrum. This measurement must be done for each independent spatial point for which one wants magnetic field data. All the measurements need to be done in a time short compared to the time scale for changes of the solar features being observed. Were it possible, one would want to record all the needed data simultaneously, since temporal variation of atmospheric seeing degrades both the image and the polarization sensitivity. Since the measurements must span four dimensions, two spatial plus polarization and wavelength, we had some freedom to design the instrument to favor some dimensions over others in terms of simultaneity. Our earlier instrument, the Haleakala Stokes Polarimeter, records a range of wavelengths spanning two spectral lines in each reading, but requires two seconds to determine the polarization state and obtains spatial information only by assembling a long sequence of measurements at single locations on the sun. The new instrument sacrifices spectral detail and accuracy in favor of greatly improved imaging characteristics. The scientific goals for this instrument were to measure surface magnetic fields with enough accuracy to permit calculations of photospheric currents, but with a field of view covering an entire typical active region, high spatial resolution, and a fast enough temporal cadence for detecting flare-associated changes in magnetic structures.

Canfield, Richard C.↗

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↗

Hydropower Operations in the Colorado River Basin: Institutional Analysis of Opportunities and Constraints

The Colorado River Basin is facing an unprecedented drought. In ongoing drought management efforts, limited attention has been paid to hydropower generation. While some studies do exist on hydropower, they are quantitative in nature and focus on calculating the reduction in megawatts generated at dams in the Basin with declining water availability. These studies simplify the complex process of hydropower generation; water availability is but one factor that impacts hydropower generation. At a more fundamental level, formal institutional arrangements, that is, laws, policies, and rules create the framework within which dams are operated and hydropower is generated. This paper conducts a comparative institutional analysis of water, environment, and energy laws and policies and changes therein to understand the constraints and opportunities faced by hydropower generation in the Colorado River Basin. To tease out the nuances in how institutional arrangements affect dam operations and hydropower generation, the comparative analysis focuses on the two largest and strategically important dams in the Basin: Hoover and Glen Canyon. This paper uses Elinor Ostrom’s Institutional Analysis and Development Framework to analyze laws and policies at three levels: constitutional-choice, collective-choice, and operational levels. Constitutional-choice level laws and policies apply to the entire Basin, whereas collective-choice level and operational level laws and policies are dam specific. Hoover and Glen Canyon Dams face similar biophysical challenges by the virtue of their location in the same river basin. Yet, despite the similarity in the biophysical setting, the analysis in this study finds that the differences in the applicability of constitutional-choice level laws along with the differences in dam specific collective-choice and operational level institutional arrangements produce a distinct set of constraints for hydropower generation at Hoover and Glen Canyon Dams. Even without a drought, water and environmental laws at both the constitutionalchoice and collective-choice levels as well as power contracts constrain hydropower generation and limit the flexibility with which Glen Canyon Dam can be operated. Water and environmental laws also impose specific water release requirements that, at times, require off-peak power generation at Glen Canyon Dam. On the other hand, even with a drought, Hoover Dam faces limited hydropower generation constraints and can operate flexibly. This is because constitutional-choice level laws and dam-specific collective-choice and operational level laws pose limited constraints for flexible daily operations at Hoover. The result is that Hoover Dam can generate hydropower at the same level as it did three decades ago and operate flexibly to provide ancillary services and peaking generation. While water and environmental laws and policies pose constraints for hydropower generation, the analysis in this study further finds that specific historic provisions within energyrelated institutional arrangements and recent changes within power contracts have maintained and even enhanced the value of hydropower to power customers. Historic institutional provisions ensure that hydropower is sold ‘at cost’ making this resource economically competitive with wholesale electricity market rates. Recent power contract modifications have resulted in the amendment of an older resale prohibition clause to expand the flexibility available to power customers in using their capacity and/or energy allocation in RTOs, ISOs, and bulk power markets. This amendment has opened up an opportunity for customers, especially Hoover power customers, to use flexible generation and ancillary services in a market environment. In addition, the extension of power contract duration to the legally maximum term has enhanced the reliability and stability of this resource for customers. In the Colorado River Basin, despite the enduring economic responsibility of power customers—where laws require customers to pay for a large portion of construction and O&M costs whether or not they actually receive hydropower— the persistent threat of a drought-induced water shortage, and constraints imposed by water and environmental laws and policies, power customers continue to invest in this resource as energy-related institutional arrangements and power contract provisions protect the reasons why they value hydropower. Lastly, the analysis in this study finds that the consequences of changes in hydropower generation for energy users, irrigators, and environmental programs in the Basin depend on how specific institutional arrangements tie electricity revenues to irrigation aid and environmental programs, and how the power contracts themselves are set up. Collective-choice level institutional arrangements create a higher level of financial dependency of irrigation aid and environmental programs on electricity revenues in the Upper Basin—the legal subdivision of Colorado River where Glen Canyon Dam is located—compared to the Lower Basin—the legal subdivision of Colorado River where Hoover Dam is located. Therefore, changes in hydropower generation or the way its revenue is collected and used will have far reacting detrimental consequences for the Upper Basin. Likewise, differences in the nature of power contracts for Glen Canyon and Hoover Dams also creates differences in the financial impact incurred by energy users when there is a reduction in hydropower generation. While this study identifies the types of impacts on resource users as a result of specific institutional arrangements, the calculation of extent of impact warrants further attention. Hydropower in the United States is in a unique position today. The strategic importance of this resource for the nation’s electricity sector is rapidly growing even as its contribution to overall electricity generation remains fairly small. This strategic importance, however, is built hydropower’s ability to operate flexibly in order to support the integration of intermittent renewable generating sources and the expansion of electricity markets. As this study shows, such flexibility may not be available at certain plants not due to the lack of water availability but because of institutional constraints. Institutional arrangements may also require dam operators to first consider high priority water uses (such as irrigation or environmental needs), which in turn may limit the ability to generate hydropower when it is most valuable or useful. Engineering and quantitative models, such as production cost models, recognize policy constraints for hydropower operations but often inadequately capture or assume away such constraints in the models. A failure to account for policy constraints in these models runs the risk of inaccurate representation of the operational flexibility and capacity available at specific hydropower plants, which can result in over/underestimation of hydropower’s ability to support the integration of variable renewable resources and address grid reliability concerns. Against this background, this paper and the analysis herein serves as an example of how we can systematically identify institutional constraints (and opportunities) that influence the flexibility in not only generating electricity at specific dams but also using this hydropower once it is generated.

13 HYDRO ENERGY↗

Estimating Future Changes of Energy Demand for Heating and Cooling Buildings at NASA Centers GC23J-1198

With its unique and trusted earth observations, NASA is a critical source in informing decisions that will help achieve the U.S. goal of Net-Zero Greenhouse Gas (GHG) Emissions by 2050. NASA’s Prediction of Worldwide Energy Resource (POWER) project facilitates the use of NASA Earth Science data holdings within the energy, agricultural, and building heating/cooling design industries. POWER packages solar and meteorological data from several NASA projects in a user friendly GIS-enabled web services system (https://power.larc.nasa.gov). As part of the development of new data products to support the energy and building heating/cooling design communities, we estimate the changes in energy required to heat and cool buildings in the future climate at 14 different NASA site locations spread throughout the continental United States, as projected by CMIP6 climate models under different emissions scenarios. Bias-corrected downscaled time series of meteorological variables are taken from NASA Earth Exchange (NEX) Global Daily Downscaled Projections (GDDP-CMIP6) downscaled climate model data. The spread between the different model projections is accounted for by analyzing both the ensemble average of 22 CMIP6 models and 6 representative models with different climate sensitivities and different interannual variability. Changes in energy use are estimated in two ways. First, changes in the total annual heating and cooling degree days (HDD and CDD, respectively) are calculated relative to the current climate. This is done at all 14 sites. Second, the downscaled time series are used as inputs into RETScreen(R), a clean energy management decision tool, to give an estimate of heating/cooling energy use for a typical office building. We use this estimation method with model data at Langley Research Center. In the next 50 years, the annual total of HDD (CDD) is projected to decrease by 8-38% (increase by 5-28%) at all sites, with the increase in CDD typically a larger magnitude the decrease in HDD. For a typical small office building at Langley Research Center, the amount of energy needed to cool increases by 33-54% and the amount of energy to heat decreases by 29-40%. POWER is working to develop long term climate data services based on these results to include in future data products to provide to users.

Bradley M. Hegyi↗

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↗

Enhancement of the Earth Science and Remote Sensing Group's Website and Related Projects

The major problem addressed throughout the term was the need to update the group's current website, as it was outdated and required streamlining and modernization. The old Gateway to Astronaut Photography of the Earth website had multiple components, many of which involved searches through expansive databases. The amount of work required to update the website was large and due to a desired release date, assistance was needed to help build new pages and to transfer old information. Additionally, one of the tools listed on the website called Image Detective had been underutilized in the past. It was important to address why the public was not using the tool and how it could potentially become more of a resource for the team. In order to help with updating the website, it was necessary to first learn HTML. After assisting with small edits, I began creating new pages. I utilized the "view page source" and "developer" tools in the internet browser to observe how other websites created their features and to test changes without editing the code. I then edited the code to create an interactive feature on the new page. For the Image Detective Page I began an evaluation of the current page. I also asked my fellow interns and friends at my University to offer their input. I took all of the opinions into account and wrote up a document regarding my recommendations. The recommendations will be considered as I help to improve the Image Detective page for the updated website. In addition to the website, other projects included the need for additional, and updated image collections, along with various project requests. The image collections have been used by educators in the classroom and the impact crater collection was highly requested. The glaciers collection focused mostly on South American glaciers and needed to include more of the earth's many glaciers. The collections had not been updated or created due to the fact that related imagery had not been catalogued. The process of cataloging involves identifying the center point location of the image and feature identification. Other project needs included collecting night images of India in for publishing. Again, many of the images were not catalogued and the database was lacking in night time imagery for that region. The last project was to calculate the size of mega fans in South Africa. Calculating the fan sizes involved several steps. To expedite the study, calculations needed to be made after the base maps had been created. Using data files that included an outline of the mega fans on a topographic map, I opened the file in Photoshop, determined the number of pixels within the outlined area, created a one degree squared box, determined the pixels within the box, converted the pixels within the box to kilometers, and then calculated the fan size using this information. Overall, the internship has been a learning experience for me. I have learned how to use new programs and I developed new skills. These These skills can help me as I enter into the next phase of my career. Learning Photoshop and HTML in addition to coding in Dreamweaver are highly sought after skills that are used in a variety of fields. Additionally, the exposure to different aspects of the team and working with different people helped me to gain a broader set of skills and allowed me to work with people with different experiences. The various projects I have worked on this summer have directly benefitted the team whether it was completing projects they did not have the time to do, or by helping the team reach deadlines sooner. The new website will be the best place to see all of my work as it will include the newly designed pages and will feature my updates to collections.

Coffin, Ashley↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗

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↗

The Nature and Evolutionary History of GRO J1744-28

GRO J1744-28 is the first known X-ray source to display bursts, periodic pulsations, and quasi-periodic oscillations. This source may thus provide crucial clues that will lead to an understanding of the differences in the nature of the X-ray variability from various accreting neutron stars. The orbital period is 11.8 days, and the measured mass function of 1.31 x 10(exp -4) solar mass is one of the smallest among all known binaries. If we assume that the donor star is a low-mass giant transferring matter through the inner Lagrange point, then we can show that its mass is lower than approximately 0.7 solar mass and probably closer to 0.25 solar mass. Higher mass, but unevolved, donor stars are shown to be implausible. We also demonstrate that the current He core mass of the donor star lies in the range of 0.20-0.25 solar mass. Thus, this system is most likely in the final stages of losing its hydrogen-rich envelope, with only a small amount of mass remaining in the envelope. If this picture is correct, then GRO J1744-28 may well represent the closest observational link that we have between the low-mass X-ray binaries and recycled binary pulsars in wide orbits. We have carried out a series of binary evolution calculations and explored, both systematically and via a novel Monte Carlo approach, the range of initial system parameters and input physics that can lead to the binary parameters of the present-day GRO J1744-28 system. The input parameters include both the initial total mass and the core mass of the donor star, the neutron-star mass, the strength of the magnetic braking, the mass-capture fraction, and the specifics of the core mass/radius relation for giants. Through these evolution calculations, we compute probability distributions for the current binary system parameters (i.e., the total mass, core mass, radius, luminosity, and K-band magnitude of the donor star, the neutron star mass, the orbital inclination angle, and the semimajor axis of the binary). Our calculations yield the following values for the GRO J1744-28 system parameters (with 95% confidence limits in parentheses): donor star mass: 0.24 solar mass (0.2-0.7 solar mass); He core mass of the donor star: 0.22 solar mass (0.20-0.25 solar mass); neutron-star mass: 1.7 solar mass (1.39-1.96 solar mass); orbital inclination angle: 18deg (7deg-22deg); semi- major axis: 64 lt-s (60-67 lt-s); radius of the donor star: 6.2 solar radius(6-9 solar radius); luminosity of donor star: 23 solar luminosity (15-49 solar luminosity), and long-term mass transfer rate at the current epoch: 5 x 10(exp -10)solar mass/yr (2 x 10(exp -10) to 5 x 10(exp -9)solar mass/yr). We deduce that the magnetic field of the underlying neutron star lies in the range of approximately 1.8 x 10(exp 11)G to approximately 7 x 10(exp 11)G, with a most probable value of 2.7 x 10(exp 11)G. This is evidently sufficiently strong to funnel the accretion flow onto the magnetic polar caps and suppress the thermonuclear flashes that would otherwise give rise to the type 1 X-ray bursts observed in most X-ray bursters. We present a simple paradigm for magnetic accreting neutron stars where X-ray pulsars, GRO J1744-28, the Rapid Burster, and the type 1 X-ray bursters may form a continuum of possible behaviors among accreting neutron stars, with the strength of the neutron-star magnetic field serving as a crucial parameter that determines the mode of X-ray variability from a given object.

Rappaport, S.↗

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↗