Search NASA⌕ Search

SEARCH · Search NASA

Results for “source zone”

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 235 records · Page 13

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Water Coolant Pool

Mechanistic assessments of radionuclide release during postulated accidents are expected to be included in advanced reactor license applications. The mechanistic source term (MST) provides an opportunity for vendors to realistically evaluate the radiological consequences of an incident, and may aid in justifying reduced emergency planning zones and plant sites. However, the development of MSTs for advanced nuclear reactors is challenging because there are numerous phenomena that can affect the transport and retention of radionuclides. As part of a trial MST assessment for a metal-fueled, pool-type sodium cooled fast reactor (SFR), led by Argonne National Laboratory, a simplified radionuclide transport code (SRT code) was developed, which includes models to estimate the quantity of fission product aerosols scrubbed in the sodium pool during postulated accident scenarios. In a pool-type SFR, when fission products are released into the coolant pool due to failure of fuel pins, most of the radionuclides are scrubbed by the coolant pool, but some have the potential to migrate to the cover gas region through entrainment within gas bubbles. The SRT code contains a model that evaluates this scrubbing behavior and calculates the fraction of fission product aerosols that reach the cover gas. Due to a lack of available validation data for sodium pool scrubbing, the U.S. Department of Energy funded an experiment at the University of Wisconsin-Madison to measure aerosol scrubbing by injecting air bubbles containing aerosol into a coolant pool. Prior to performing an experiment with liquid sodium, a water loop experiment was performed. Their experiment evaluated the effect of changing the aerosol size, aerosol density, aerosol concentration, bubble size, and pool depth on the aerosol scrubbing efficiency of the pool. In this benchmark experiment, the base tests were conducted by repeated tests of isolated bubbles. Afterwards, more prototypic tests with bubble swarms were performed to evaluate the interactions between the bubbles. The bubble swarm test was able to confirm that a larger amount of aerosol scrubbing occurred than the single bubble test. It was also confirmed that as the bubble size, aerosol density, and pool height increase, the extent of pool scrubbing also increases and does not change with the aerosol concentration. In addition, since the degree of scrubbing is the lowest at aerosol sizes between 0.01 and 1 μm, that is, the largest amount of aerosol is emitted, it was confirmed that the analysis of this size in MST is the most important. This benchmark experiment informs the direction of future sodium experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Local Resource Assessment: Beetle-Killed Spruce in Alaska

This report synthesizes forest inventory and spatial data from federal and state sources to present a snapshot of beetle-killed spruce across Southcentral and Interior Alaska. This natural resource is under investigation for use in cellulose-based products, particularly building insulation. Research funded by the Department of Energy's ARPA-E program, has successfully developed insulation - Celium - using beetle-killed spruce combined with mycelium, the root network of fungi. The resulting insulation product can be manufactured locally and used across all Alaskan climate zones for residential, commercial, and shipping applications. Removing standing dead wood also contributes to wildfire risk reduction and mitigates associated public health impacts. This Local Resource Assessment, funded by the Denali Commission, contextualizes and reconciles unaligned forestry data sources, establishing a pathway to refined, replicable analysis and, ultimately, presenting estimates of beetle-killed spruce biomass to inform the potential scale and development of a Celium industry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Multiwavelength Observations of the Blazar VER J0521+211 during an Elevated TeV Gamma-Ray State

We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10 –7 photons m –2 s –1 , or ~37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗

Radiological Releases from Novel Fuel Forms in Advanced Reactors During Severe Accidents for Consequence Analyses

Various advanced reactor developers are exploring the potential for reductions in the size of physical security forces and emergency planning zones. These reductions are based on robust fuel forms and inherently safe reactor designs. However, such reductions in physical protection measures could increase the risk of sabotage. To assess the possibility of reducing these measures, sabotage-induced radiological consequence analyses were carried out. These analyses considered accident scenarios that were beyond design basis accidents and overly conservative (Shah, 2025a; Shah, 2025b; Shah and Hartanto, 2026), yielding very large release fractions. These fractions, which can be used to evaluate physical protection and emergency planning requirements, have been crudely determined and applied as demonstrations for a sodium-cooled fast reactor (SFR) (Shah and Hartanto, 2025a), a high-temperature gas-cooled reactor (HTGR) (Shah and Hartanto, 2025b), a heat pipe–cooled reactor (HPR) (Shah and Hartanto, 2025c), and a molten salt–cooled reactor (MSR) (Shah et al., 2026). A Sandia National Laboratories (SNL) team used MELCOR—a fully integrated severe accident analysis code—to demonstrate the code’s capability to analyze advanced (i.e., not light water–cooled) reactors (including a fluoride salt–cooled high-temperature reactor [FHR]) and calculate radiological releases to the environment during severe accidents (Wagner et al., 2022a, 2022b, 2022c, 2023a, and 2023b). Although the analyses were carried out to demonstrate MELCOR’s growing capability, the release source terms were estimated for advanced reactors, providing valuable insights into the accident progression and radiological releases. These findings from prior SNL studies, including estimated source terms and related sensitivity studies, were leveraged to derive source terms for postulated sabotage-induced accidents. Insights from these sensitivity studies informed the scaling of SNL’s estimated source terms for the defined accident scenarios. The derived release fractions for the severe accident scenarios for the respective reactor designs can be used to perform more nuanced dose consequence analyses to evaluate the reactors’ physical protection and emergency planning zone requirements. These analyses are in accordance with the risk-informed, performance-based approach proposed under 10 CFR Part 53. This study builds on the prior source term analyses and associated sensitivity studies by SNL to derive time-dependent and design-informed release fractions. Section 2 describes the diverse advanced reactor designs analyzed by the SNL team. Section 3 discusses the severe accident analyses, the release fractions calculated, and the limitations and assumptions of the demonstration project. Section 4 presents the release percentages derived for the hypothetical sabotage-induced severe accidents at the advanced reactors. Section 5 summarizes the study’s findings and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

In-Situ Nanoscale Focusing of Extreme Ultraviolet Solid-State High Harmonics

Extreme ultraviolet light delivering radiation with a wavelength shorter than approximately 100 nm is now available from solid-state sources. However, despite exceptional progress, efficient focusing of extreme ultraviolet photons to their ultimate diffraction limit remains a formidable challenge because of the precision of the focusing by curved, optical surfaces. Here we integrate coherent short-wavelength high-order harmonics from a MgO crystal, with a high-numerical-aperture, nanostructured, focusing element etched onto the surface of the crystal itself. We focus extreme ultraviolet light, seventh harmonic of an 800-nm laser, with a zone plate of numerical aperture 0.35, down to a waist radius of 150 nm, with 18% focusing efficiency. The estimated intensity approaches 10 7 W/cm 2 . Future developments may demonstrate nanoscale laser ablation and miniaturization of extreme ultraviolet coherent sources on a chip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of the Table of Initial Isolation and Protective Action Distances for the 2024 Emergency Response Guidebook

The transportation of hazardous materials creates numerous opportunities for the release of toxic substances into the environment, whether caused by traffic accidents, train derailments, equipment failures, or human error. Such releases can pose acute hazards to the general public and to emergency response personnel who are the first to arrive at the scene. To help first responders determine whether a shipment is potentially hazardous and decide what actions should be taken if a toxic spill does occur, the Emergency Response Guidebook (ERG) is published by the U.S. Department of Transportation (DOT), Transport Canada, and the Secretariat of Transport and Communications of Mexico; with contributions from Centro de Informaciòn Quìmica para Emergencias of Argentina. The most recent version is the 2024 edition of the ERG (ERG 2024), titled 2024 Emergency Response Guidebook (ERG2024). The ERG provides essential information about firefighting, spill response, and potential public health effects. For chemicals that are toxic by inhalation (TIH) and chemicals that produce TIH gases upon reaction with water (TIH by water reactivity or TIHWR), the ERG provides initial isolation distances (IIDs) and protective action distances (PADs). The IID defines the radius of the zone around the spill that should be accessed solely by people who are directly involved in emergency response. The PAD is the distance downwind of the source of the release within which persons should be either evacuated or sheltered in place, depending on the severity of the incident and the nature of the population (e.g., density, age, health).

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

In Situ Laser Induced Florescence Measurements of Fuel Dilution from Low Load to Stochastic Pre Ignition Prone Conditions

This work employs a novel laser induced fluorescence (LIF) diagnostic to measure fuel dilution in a running single cylinder research engine operated at stochastic pre ignition (SPI) and non-SPI prone conditions. Measurements of LIF based fuel dilution are quantified over a range of engine loads and fuel injection timings for two fuels. The in situ LIF measurements of fuel/lubricant interactions illustrate regions of increased fuel dilution from fuel-wall interactions and is believed to be a fundamental underpinning to generating top ring zone liquid conditions conducive to SPI. Furthermore, a novel level of dye doped in the fuel, between 50 to 500 ppm was used as the fluorescence source, at engine operating speed of 2000r/min from 5 to 18 bar of IMEPg injection timings was swept for two fuels of varying volatility. The direct real time LIF measurements highlight that there are non-linear trends in fuel dilution beyond simple dependencies of fuel volatility, injection duration or injection timing, suggesting that further understanding of spray interaction with engine surfaces and the turbulent field are needed to quantify fuel dilution effects that are conducive to SPI. Moreover, results show the potential of this diagnostic technique as an additional tool for quantifying spray and fuel mixing in fundamental studies deployable across a variety of engine loads from law to full load.

47 OTHER INSTRUMENTATION↗

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗

Characterizing major and trace element compositions in fallout melt glass from a near-surface nuclear test

We report the chemical and isotopic compositions of fallout melt glasses from nuclear tests contain a range of information constraining the physical conditions within the fireball and the mechanisms of fallout formation but historic studies tended to exclude the behavior of stable major and trace elements. Here, we present a large study specifically focused on major and trace element relationships within a population of macroscale fallout samples from a single event. We interpret these data to better constrain how fallout melt glass formation in near surface environments is influenced by that environment and demonstrate how major and trace element abundances can provide useful insights into chemical processes within the fireball. Data confirm that the uranium in the fallout glass population derives from two isotopically distinct endmembers: isotopically enriched uranium (presumably from the weapon), and natural composition uranium that may be a combination of anthropogenic and environmental materials from within the blast zone. The similarity between major and trace element concentrations in fallout and corresponding local soils from the event site confirm the local soils as the most probable source of entrained material into the fireball and the source of carrier material into which the bomb vapor was incorporated. The lack of correlation between major and trace element abundances with size indicates that volatility driven processes, such as condensation from the fireball, do not control the composition of macroscale fallout melt glass. Although the fallout has major and trace element chemical characteristics broadly similar to those of the local, associated soils, some systematic differences are observed between the two populations. Fallout melt glass is depleted in volatile elements such as K, Na, Tl and Pb, consistent with heating to temperatures above ~1000 °C for 3-10 s. This is supported by the results of laser heating experiments performed on rhyolitic soil at temperatures (1600-2200 °C) and timescales (1-120 s) that are broadly relevant to fallout formation conditions. Relative enrichments of metals such as Cu and Co do not correlate with the abundance of uranium, suggesting that fallout also records input of near field anthropogenic materials. Our observations suggest that major chemical features can be related to processing in the fireball and used to inform the thermal-chemical evolution of the system. Ultimately, these data are consistent with a fallout formation mechanism that involves rapid melting of surface materials to form carrier material melts with minor incorporation of bomb vapor and a degree of volumetric volatile loss due to heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗

Variations of Stress Parameters in the Southern California Plate Boundary Around the South Central Transverse Ranges

We examine the stress field in Southern California with a focus on the region near the South Central Transverse Ranges (SCTR), using a refined stress inversion methodology. We independently invert declustered and aftershocks focal mechanisms from the years 1981 to 2017. Comparisons between the stress field orientations and the stress ratios from these inversions provide information on local sources of the stress field. Regionally, the SHmax tends to be toward NNE and the stress ratios vary from transtensional stress regime near the Eastern California Shear Zone (ECSZ), to strike-slip faulting near the SCTR, and toward transpression near the Western Transverse Ranges. Detailed analysis of stress parameters near the SCTR indicates deviations from the regional strike-slip faulting including the following: (1) The San Bernardino Mountain area shows S Hmax direction toward NNW and transpressional stress components likely associated with the motions of the San Andreas Fault and ECSZ; (2) the Cajon Pass and San Gorgonio Pass show transpressional stress regime near the bottom of the seismogenic zone likely associated with the elevated topography; and (3) in Crafton Hills, rotation of the principal stress plunges and SHmax direction and transtensional stress regime below ~10 km, along with lower estimated apparent friction coefficient, suggest a weak fault possibly associated with deep creep. The performed multiscale analysis resolves effects of regional and local loadings. The average stress parameters in the analyzed 37 years do not show significant temporal variations of stress field near the SCTR region.

58 GEOSCIENCES↗

A Study of Cost-Saving Potential of Load Flexibility Measures in Grid-Interactive Multifamily Buildings

With recent advances in smart technologies, more and more smart devices are penetrating the residential and commercial buildings market. The introduction of these smart devices is also helping IoT companies emerge with load aggregator roles in the sector. With more utility companies on the track of supporting OpenADR protocols, the aggregators could play a significant role in providing load flexibilities by automatically responding to demand response (DR) events and coordinating load flexibility measures between customers. This would benefit utility companies by reducing stress on the grid during critical peak demand hours as well as customers by allowing them to utilize utility rate structures advantageous to those able to reduce electric usage during high-demand hours. This study evaluates cost and energy savings from adopting multiple load flexibility measures in multifamily buildings. Combinations of different load flexibility measures, including space temperature floating, light dimming, automatic window shading, and water heater temperature floating, are considered. The simulations are performed using OpenStudio®, an open-source U.S. Department of Energy (DOE) simulation platform. For the case study, we used a midrise apartment building with weather conditions from Denver, Colorado. To compare climate zone differences, simulations were also performed for Los Angeles, California, and Chicago, Illinois. Initial results indicate that the application of automated load flexibility measures without careful consideration of dispatching strategies and DR program enrollments could significantly affect the savings. To get meaningful cost savings, aggregators need to encourage tenant awareness to curtail energy usage through occupant behavior in addition to dispatching automatic load flexibility measures. The outcomes from this study are believed to help load aggregators understand the risks and benefits of load flexibility opportunities.

building energy modeling↗

OB200-DV-1 Treatability Testing: Final Results

The Hanford Site in Washington state previously generated plutonium for nuclear weapons. During operations, radionuclide byproducts and chemical process fluids were intentionally and/or unintentionally released to the subsurface, resulting in more than 800 contaminated waste sites across the Central Plateau, where historical chemical separations and waste management activities took place. As the Hanford Site mission transitioned from operations to site cleanup, remediation of the vadose zone and groundwater became a priority. However, given the depth of the unsaturated zone contamination above the groundwater, the unique nature of the waste, and the continuing impacts on groundwater quality, technologies needed to be identified and evaluated for in situ remediation in the deep vadose zone (DVZ). A laboratory treatability study has been completed to evaluate site-relevant effectiveness for nine in situ technologies that may be used to treat continuing sources of contaminants in specific areas of the Central Plateau waste sites that are grouped into the 200-DV-1 Operable Unit (OU). The 200-DV-1 OU was established in 2010 to address 43 Central Plateau waste sites with complex DVZ remediation challenges. Eight of these technologies were identified through a prescreening effort that evaluated remedial technologies potentially applicable to DVZ contamination in the Central Plateau . These eight technologies were selected for further study based on site specific knowledge gaps about their effectiveness. A ninth technology was added to the treatability study based on new information from separate laboratory investigations (conducted following the prescreening effort) demonstrating the technology’s potential effectiveness (see Section 1.2 for more information) and value for inclusion in the treatability study.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗

Investigating the Role of Accident Tolerant Cladding on Source Term Reduction for High-Burnup PWRs Using MELCOR

The use of accident tolerant fuel (ATF) cladding can increase coping times during and beyond design basis accidents. While such gains may be incremental, they provide a margin that can potentially be recovered to enable high-burnup (HBU) operation. Realizing such a margin requires demonstrating that the combination of HBU and ATF has not led to an overall increase in source term. This study investigates the influence of cladding technology (Zr-based, Cr-coated Zr, and FeCrAl) and fuel cycle length (18 and 24 months) on radiological dose at the boundary of the exclusion zone for a four-loop pressurized water reactor to investigate whether ATF claddings can provide such benefits. We analyze a recovered large break loss-of-coolant accident scenario to investigate the impact of transient timescale on the benefits of such coping time increases. The simulations have been performed using the MELCOR and MELCOR Accident Consequence Code System codes. For the cases analyzed, increased fuel cycle length did not necessarily increase radionuclide release and hydrogen generation, as these were found to be sensitive to the core power distribution. Similarly, off-site dose consequence is dominated by short-lived radionuclides that tend to saturate earlier in the burnup, so higher burnup operation did not necessarily increase the source term for the phenomena and transients analyzed here. Delays in recovery of the lowpressure safety injection system increase hydrogen production and radionuclide release, especially between 780 s and 1620 s, due to the nonlinear oxidation and core degradation behavior. Results show that Cr-coated Zr enhances safety by delaying heatup and gap release. Here, when uncertainty propagation on oxidation properties is considered, FeCrAl exhibits the lowest overall radionuclide release and off-site dose throughout the spectrum. However, while the considered “base model” performance is superior under delayed injection scenarios, upper-bound cases display hydrogen generation risk comparable to the Zr-based cladding.

Accident Tolerant Fuel↗