Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elevation models”

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 577 records · Page 32

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗

Microstructure and properties of additively manufactured Al–Ce–Mg alloys

Additive manufacturing of aluminum alloys is largely dominated by a near-eutectic Al-Si compositions, which are highly weldable, but have mechanical properties that are not competitive with conventional wrought Al alloys. In addition, there is a need for new Al alloys with improved high temperature properties and thermal stability for applications in the automotive and aerospace fields. In this work, we considered laser powder bed fusion additive manufacturing of two alloys in the Al–Ce–Mg system, designed as near-eutectic (Al–11Ce–7Mg) and hyper-eutectic (Al–15Ce–9Mg) compositions with respect to the binary L → Al + Al 11 Ce eutectic reaction. The addition of magnesium is used to promote solid solution strengthening. A custom laser scan pattern was used to reduce the formation of keyhole porosity, which was caused by excessive vaporization due to the high vapor pressure of magnesium. The microstructure and tensile mechanical properties of the alloys were characterized in the as-fabricated condition and following hot isostatic pressing. The two alloys exhibit significant variations in solidification structure morphology. These variations in non-equilibrium solidification structure were rationalized using a combination of thermodynamic and thermal modeling. Both alloys showed higher yield strength than AM Al-10Si-Mg for temperatures up to 350 °C and better strength retention at elevated temperatures than additively manufactured Scalmaloy.

36 MATERIALS SCIENCE↗

Continued Investigations of Respirable Release Fractions for Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using relatively high backfill pressures (up to approximately 800 kPa) in the canister to enhance internal natural convection. This pressure differential offers a comparatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.89 mm (0.350 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to demonstrate a new capability to characterize SCCs under well-controlled boundary conditions. Modeling efforts were also initiated that evaluate the depletion of aerosols in a commercial dry storage canister. These preliminary modeling and ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation [Abstract]

The Oak Ridge National Laboratory (ORNL) Health Physics Research Reactor (HPRR) was a research reactor designed and built at ORNL in 1961. The critical assembly was using a highly enriched uranium and molybdenum alloy as the fuel, and it could be operated in steady-state or burst modes. The reactor was used for about 25 years to produce a lot of publications related to dosimetry, radiobiology and radiation detectors testing before its decommissioning in 1987. In recent years, the idea of using legacy operation data from the to create a valuable critical accident alarm system shielding benchmark arose. Such a benchmark has been submitted to the International Criticality Safety Benchmark Experiment Project (ICSBEP) Technical Review Group for a potential inclusion in the 2022 version of the handbook. Another way to use the valuable data from the operation of the HPRR is to evaluate the feasibility of the creation of a subcritical or prompt supercritical benchmark for inclusion in the ICSBEP or the International Reactor Physics Experiments Evaluation Project (IRPhEP) handbooks. To initiate a burst, the HPRR had to be operated in a slightly subcritical state for a few minutes. Then, the insertion of the burst control rod would greatly increase the reactivity of the system and start the burst. No critical configuration of the HPRR critical assembly could be located. The only information available concerns stable subcritical and prompt supercritical states, found in a burst experiments’ logbook. In the recovered logbook pages, information about 8 different bursts is available. The information includes the rods positions before and during a burst, the recorded subcritical reactor period and reactivity, and the burst fission yield derived from the temperature elevation sulfur pellet irradiation analysis. By using the HPRR logbook information and the as-built drawings of the critical assembly, a highly detailed model of the HPRR was created with SCALE 6.2.4/KENO-VI. Eight KENO-VI models were created to replicate the sub-critical assembly configurations described in the eight bursts from the recovered logbook pages. KENO-VI calculates k eff and it can be linked to a reactivity value in cents by using the delayed neutron fraction B eff , also calculated by KENO-VI. KENO-VI can also be used to model the prompt super-critical configurations of the HPRR and to assess the similarity with the burst measurements by comparing the calculated k eff and the measured fission yields between each burst. Unfortunately, high uncertainty exist and the obtained discrepancies between experiments and calculation results are high, compromising the creation of a valuable critical benchmark from HPRR operation data. The reasons of the discrepancies and potential ways to solve them are explored.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mechanical response and microstructure evolution from multiaxial stress relaxation in textured zircaloy nuclear cladding

Modern energy markets provide motivation for nuclear reactors to increase operational flexibility and load following capacity of existing reactor designs. The operational power changes result in fuel pellet expansion and contraction resulting in pellet-cladding mechanical interaction (PCMI) and high cladding stresses that are strongly related to stressed corrosion cracking (SCC) which can lead to fuel rod failure. This work aims to quantify the stress fields in cladding with a focus on multiaxial stress relaxation response from an imposed strain and the resulting microstructure changes in textured zircaloy claddings. Stress relaxation tests resulting in multiaxial stress fields are performed to measure the material response of Zircaloy-4 nuclear cladding. Multiple experiments using samples from the same piece of cladding are performed so samples with the same processing history may be removed at strategic points throughout mechanical testing for microstructure characterization. The microstructure evolution is analyzed with mechanical response to evaluate the static recovery impact on crystallographic texture and grain morphology. Experiments are performed at Idaho National Laboratory using an experimental device capable of applying internal pressure and independent axial load to commercial nuclear cladding at elevated temperatures with stress- or strain-controlled experimental capability. X-ray diffraction and electron backscatter diffraction are used for microstructure analysis. Rheological models are used to visualize the multiaxial stress-strain relationships and isolate elastic, plastic, and viscoplastic material properties for analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Estimating soil N 2 O emissions induced by organic and inorganic fertilizer inputs using a Tier-2, regression-based meta-analytic approach for U.S. agricultural lands

Consistent methods are essential for generating country and region-specific estimates of greenhouse gas (GHG) emissions used for reporting and policymaking. The estimates of direct N 2 O emissions from U.S. agricultural soils have primarily relied on the use of emission factors (EFs, Tier-1) and process-based models (Tier-3). However, Tier-1 estimates are relatively crude while Tier-3 calculations can be costly. This work addressed this gap by developing a Tier-2, regression-based approach by leveraging a meta-database containing 1883 field N 2 O observations together with environmental and management covariates from 139 studies. Our results estimated higher monthly soil N 2 O emissions (N 2 O m , kg N/ha) during the growing season (0.38) than the fallow period (0.15), highlighting the importance of considering measurement periods when utilizing meta-databases for analyzing N 2 O drivers. Significantly different N 2 O m were found for tillage practices (conventional > no-till: 0.42 > 0.27), fertilizer type (liquid > solid manure: 0.55 > 0.32), and soil texture (fine > coarse: 0.36 > 0.22). The comparisons of the influence of crop type and rotation, water management, and soil order on N 2 O emissions are complicated by regional data availability and interactions among different factors. Additionally, the finding that N 2 O emissions reported based on area (N 2 O m ), N input rate (EF), or yield can alter treatment rankings underscores the need to establish transparent criteria for rewarding or discouraging regionally-based management practices using N 2 O metrics. Finally, we show how General Linear Models (GLMs) can be used to estimate country and regional Tier-2 N 2 O m using a suite of covariates. Our GLMs identified tillage, water management, N input type and rate, soil properties, and elevation as the most influential covariates for the conterminous U.S. The limited accuracy of regional-scale GLMs, however, suggests the need to further improve the quality and availability of GHG and covariate data through concerted efforts in data collection.

54 ENVIRONMENTAL SCIENCES↗

Throughfall-Reduction Drying Effects on Soil Respiration, Moisture, and Temperature in Four Distinct Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that throughfall reduction (i.e. chronic drying) would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We conducted a field study in four distinct lowland Panamanian forests that vary in background rainfall and soil fertility. In each forest we constructed throughfall reduction structures that reduce moisture inputs to plots by ~50% (see methods). In drying and paired control plots we measured soil CO2 fluxes, soil moisture (volume/volume), soil temperature (°C), air temperature (°C), and forest floor biomass during wet-dry cycles (2015 – 2022). We also collected litterfall biomass and surveyed all tree species for diameter at breast height (dbh), and identified tree species in 2018 and 2019.Results:We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks. Datasets included: Datasets included here include .csv and .xls files for soil respiration (carbon dioxide flux) collected using the LiCor (see methods), air temperature °C , soil moisture (vl/vl), and soil temperature °C. Soil moisture and soil temperature were collected at 0-10cm depths using hand-held probes (see methods). Data are on ~3 week timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗

Throughfall-Reduction Drying Effects on Soil Moisture and Soil Temperature at 30-Minute Intervals Four Lowland Panamanian Forests

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (Cusack et al. 2023). Datasets included: Datasets here include .csv and .xls files for volumetric soil moisture (volume/volume) and soil temperature (°C). Data were taken at 30-minute intervals at 5cm and 20cm depths using permanent probes and data loggers (see methods). Data presented are from 2018-2019. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

Light and Elevated Temperature Induced Degradation (LeTID) in a Utility-Scale Photovoltaic System

Here, we present a detailed case study of degradation in monocrystalline silicon photovoltaic modules operating in a utility-scale power plant over the course of approximately three years. We present the results of degradation analysis on arrays within the site, and find that five of the six arrays degraded faster than the best performing array, even though the arrays consist of modules of the same manufacturer and model. We also describe the results of extensive laboratory characterization of modules returned from the field, including module- and cell-level current–voltage characterization, luminescence imaging, and accelerated testing. The laboratory test results and the field performance are consistent with light and elevated temperature induced degradation (LeTID). Notably, we observe differences in back contact technology between affected and unaffected modules. This article also demonstrates a method to identify possible LeTID degradation in the field and confirm the result with laboratory testing of a small number of modules.

14 SOLAR ENERGY↗

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗

Distribution Development for Residual Inventory at the New York West Valley Site - 20400

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351 ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Operations led to contamination in a number of facilities and locations. Some of that contamination has migrated from waste disposal zones to other layers, formations, and features on and off the WNYNSC. Phase I decommissioning activities are ongoing and involve the removal of a number of areas and structures that have been associated with contamination. The purpose of this work is to outline the approach for characterizing contamination not associated with disposed wastes, contaminated structures, or specific releases. In this work, the term, residual radiological activity, is used to describe environmental contamination that exists subsequent to the completion of Phase I decommissioning activities, that is not associated with disposed wastes, contaminated structures, or specific releases. Contamination from the Site was quantified relative to data that characterize the concentrations of radionuclides that exist in background. Background concentrations are those present in the area but having no influence from Site related activities. The existence of residual radiological activity that is elevated relative to background has the potential to contribute to future risks to human health and the environment. As a consequence, the residual inventory information is used to inform the West Valley Probabilistic Performance Assessment (PPA) model to characterize potential future risks to human health and the environment. The centralized West Valley Data Management System (DMS) was the source of information for the data assembled in this analysis. The DMS is a fairly large compilation consisting of thousands of records from investigation studies, with sample dates ranging from 1990 to present. Samples from monitoring wells, boreholes, geoprobe studies, surface water, surface soils, storm water outfalls, ventilation stack filters, plant and animal tissues, and more are included in the DMS. Results are typically reported in units of activity per unit volume. For the purpose of the analyses presented here, all results were converted into consistent units of pCi per unit volume. Since 1990, data have been collected from various locations across the WNYNSC at different times with varying frequency over the course of several decades. As a consequence, a number of potential issues can arise with respect to the assembly of a dataset that is deemed adequate for the characterization of residual radiological activity. These issues were assessed and resolved to the extent possible through careful consideration of the properties of the distributions. The intent was to use data which characterize the current state of the Site. Radionuclides can be designated to one of several groups depending on their origin. In this work the groups considered were 1) Naturally Occurring Radioactive Material (NORM), 2) fallout, and 3) Other (including power plant, medical research, etc). This grouping is a useful construct with respect to the interpretation of fixed laboratory results. For example, NORM radionuclides that exist within a decay chain should have approximately equivalent distributions of concentrations if they are representative of background conditions. Insights such as these can be used as a check to identify sample results that need to be further investigated or omitted due to issues associated with reported results from fixed laboratory analyses. This type of analysis provided a foundation for the assessment of the adequacy of sample results for use in subsequent components of an assessment. The general process for the assessment of residual radiological contamination at the Site consists of a sequence of several steps. First, for each analyte, several statistical tests were performed to assess the weight of evidence against the null hypothesis that the mean of the distribution of concentrations was equal to zero. If the mean of the distribution of concentrations for a given radionuclide was not found to be greater than zero, then it was removed from consideration as a component of the residual radiological contamination. If there was significant evidence to reject the hypothesis of the mean being equal to zero, the second step was to compare the distribution of the data from the Site to that of the corresponding background. A suite of tests was used to compare the distributions of the site and background data. The results of these tests were collectively used to determine if site data are elevated relative to background. The third step was to develop distributions using a Bayesian framework to characterize the distribution of mean of the increment present above background for each of the radionuclides. The Bayesian model implemented allowed for the comparison of site-specific records to background concentrations to better approximate contamination attributed to the Site. A final screening step was employed for radionuclides that exceed background. This screening step compared 95% upper confidence limits (UCLs) from the increment distribution developed in the previous step to the risk screening levels. This approach yields a list of analytes that were determined to be elevated relative to background.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Model 9975 O Ring Fixture LongTerm Leak Performance (FY2021 Status Report)

Leak testing experiments to monitor the aging performance of Viton® GLT and GLT-S O-rings used in the model 9975 shipping package has been ongoing since 2004 at Savannah River National Laboratory. Seventy tests using mock-up 9975 primary containment vessels (PCVs) with GLT O-rings were assembled and heated to temperatures ranging from 200 to 450 °F. Due to material substitution, fourteen tests with GLT-S O-rings were initiated in 2008 and heated to temperatures ranging from 200 to 400 °F. The conditioning temperatures are elevated compared to the calculated maximum O-ring temperature in a 9975 package in storage, 158 °F, to accelerate aging and enable observations of O-ring failures in a reasonable time frame. The fixtures are leak tested periodically and all GLT O-ring fixtures aging at 350 °F or above have failed to maintain a leak-tight seal. Eight GLT O-ring fixtures aging at 300 °F have failed after 2.8 to 5.7 years at temperature while the remaining fixtures at 300 °F were retired from testing following more than five years without failure. There has been one GLT O-ring fixture which failed after 13.4 years aging at 200 °F. However, 19 other GLT O-rings aging at 200 °F have remained leak-tight for over 14 years and remain in test. No failures have yet been observed in GLT O-ring fixtures aging at 270 °F for 9.5 years. All GLT-S O-ring fixtures aging at 300 °F or above have failed their leak test. No failures have yet been observed in GLT-S O-ring fixtures aging at 200 and 250 °F for 11.5 years. Data from the O-ring fixtures are generally consistent with results from compression stress-relaxation testing, and provide confidence in the predictive models based on those results. However, uncertainty exists in extrapolating these elevated temperature results to the lower temperatures of interest for normal storage in K-Area Complex (KAC). Oxygen consumption testing, which includes results from temperatures near KAC normal storage temperatures, is ongoing to provide further confidence and corroborate these extrapolations. The collective data from these test efforts suggest the minimum O-ring service life at KAC normal storage conditions should be at least 34 years for GLT and GLT-S O-rings. Measurement of compression set in O-rings removed from failed fixtures, compared to that from KAC surveillance O-rings, indicate significant margin remains for O-rings still in service in 9975 packages in KAC. Aging and periodic leak testing will continue for the remaining 26 mock-up PCV fixtures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

On the use of air temperature and precipitation as surrogate predictors in soil respiration modelling

Soil respiration (R S ), the soil-to-atmosphere CO 2 flux that is a major component of the global carbon cycle, is strongly influenced by local soil temperature (T soil ) and water content (SWC). Regional to global-scale R S modelling thus requires this information at local scales, but few high-quality, wall-to-wall (global) T soil and SWC data exist. As a result, such modelling efforts commonly use air temperature (T air ) and monthly precipitation (P m ) as surrogate predictors, but their site-scale accuracy and potential bias are unknown. In this report we used monthly data from 880 sites across a wide variety of different environmental conditions (i.e., climate, ecosystem type, elevation, vegetation leaf habit and drainage conditions) to determine the suitability of T air as a surrogate for T soil , and data from 507 sites to examine the suitability of P m as a surrogate for SWC. Site-specific linear and second-order exponential non-linear models were compared using model evaluation metrics (i.e., slope, p-value of slope, root mean square error [RMSE], index of agreement and model efficiency). We found that T soil and T air are highly correlated and explain similar R S variability. In contrast, P m is not a good surrogate for SWC, even though P m explains a similar amount of R S variability to SWC. The wide variability in the site-specific relationships between R S and SWC means that no single relationship can be used for large-scale modelling. The results from this study support the use of T air in continental-to-global scale R S models, and highlight the urgent need for continental-to-global scale SWC datasets for the modelling and evaluation of future soil carbon dynamics under global climate change.

54 ENVIRONMENTAL SCIENCES↗

Inferred inflow forecast horizons guiding reservoir release decisions across the United States

Abstract. Medium- to long-range forecasts often guide reservoir release decisions to support water management objectives, including mitigating flood and drought risks. While there is a burgeoning field of science targeted at improving forecast products and associated decision support models, data describing how and when forecasts are applied in practice remain undeveloped. This lack of knowledge may prevent hydrological modelers from developing accurate reservoir release schemes for large-scale, distributed hydrology models that are increasingly used to assess the vulnerabilities of large regions to hydrological stress. We address this issue by estimating seasonally varying, regulated inflow forecast horizons used in the operations of more than 300 dams throughout the conterminous United States (CONUS). For each dam, we take actual forward observed inflows (perfect foresight) as a proxy for forecasted flows available to the operator and then identify for each week of the year the forward horizon that best explains the release decisions taken. Resulting “horizon curves” specify for each dam the inferred inflow forecast horizon as a function of the week of the water year. These curves are analyzed for strength of evidence for contribution of medium- to long-range forecasts in decision making. We use random forest classification to estimate that approximately 80 % of large dams and reservoirs in the US (1553±50 out of 1927 dams with at least 10 Mm3 storage capacity) adopt medium- to long-range inflow forecasts to inform release decisions during at least part of the water year. Long-range forecast horizons (more than 6 weeks ahead) are detected in the operations of reservoirs located in high-elevation regions of the western US, where snowpack information likely guides the release. A simulation exercise conducted on four key western US reservoirs indicates that forecast-informed models of reservoir operations may outperform models that neglect the horizon curve – including during flood and drought conditions.

54 ENVIRONMENTAL SCIENCES↗

Model 9975 O-Ring Fixture Long-Term Leak Performance (FY2022 Status Report)

Leak testing experiments to monitor the aging performance of Viton® GLT and GLT-S O-rings used in the model 9975 shipping package has been ongoing since 2004 at Savannah River National Laboratory. Seventy tests using mock-up 9975 primary containment vessels (PCVs) with GLT O-rings were assembled and heated to temperatures ranging from 200 to 450 °F. Due to material substitution, fourteen tests with GLT-S O-rings were initiated in 2008 and heated to temperatures ranging from 200 to 400 °F. The conditioning temperatures are elevated compared to the calculated maximum O-ring temperature in a 9975 package in storage, 158 °F, to accelerate aging and enable observations of O-ring failures in a reasonable time frame. The fixtures are leak tested periodically and all GLT O-ring fixtures aging at 350 °F or above have failed to maintain a leak-tight seal. Eight GLT O-ring fixtures aging at 300 °F have failed after 2.8 to 5.7 years at temperature while the remaining fixtures at 300 °F were retired from testing following more than five years without failure. There has been one GLT O-ring fixture which failed after 13.4 years aging at 200 °F. However, 20 other GLT O-rings aging at 200 °F have remained leak-tight for over 14 years and remain in test. No failures have yet been observed in GLT O-ring fixtures aging at 270 °F for 10.4 years. All GLT-S O-ring fixtures aging at 300 °F or above have failed their leak test. No failures have yet been observed in GLT-S O-ring fixtures aging at 250 °F for 12.5 years. This year, one GLT-S O-ring fixture failed after 12.4 years aging at 200 °F. Data from the O-ring fixtures are generally consistent with results from compression stress-relaxation testing, and provide confidence in the predictive models based on those results. However, uncertainty exists in extrapolating these elevated temperature results to the lower temperatures of interest for normal storage in K-Area Complex (KAC). Oxygen consumption testing, which includes results from temperatures near KAC normal storage temperatures, is ongoing to provide further confidence and corroborate these extrapolations. The collective data from these test efforts suggest the minimum O-ring service life at KAC normal storage conditions should be at least 34 years for GLT and GLT-S O-rings. Measurement of compression set in O-rings removed from failed fixtures, compared to that from KAC surveillance O-rings, indicate significant margin remains for O-rings still in service in 9975 packages in KAC. Aging and periodic leak testing will continue for the remaining 25 mock-up PCV fixtures.

36 MATERIALS SCIENCE↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗