Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Representatives”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electrochemical corrosion under controlled redox conditions

A method for predicting corrosion rates of a material during service conditions is provided, the method having the steps of determining a first phase composition of the material; exposing the material to service conditions chemical environment; applying an electrical potential to the exposed material to represent the solution redox; identifying ranges of the applied potential that correspond to different corrosion behaviors of the material; quantifying current and surface electrical properties during corrosion; and determining a second phase composition of the material to identify corroded phases. Also provided is a method for determining radionuclide source terms, the method having the steps of supplying a multiphase metallic waste containing the radionuclides; immersing the waste in a solution representing repository chemistry conditions; and oxidizing the immersed waste for a period of time and at particular imposed voltages representing solution redox values to establish a steady current representing corrosion rate of the waste.

Gattu, Vineeth Kumar↗
DOE OSTI2017

Installing Deep Groundwater Monitoring Wells in a Volcanic Dominated Rift Basin: Los Alamos National Laboratory, New Mexico - 17594

Los Alamos National Laboratory (LANL, or the Laboratory) is located in northern New Mexico near the west margin of the Espanola Basin. LANL is underlain by a complex basin-fill sequence of alluvial fan deposits, volcanic tuff, basaltic and dacitic lavas, and riverine deposits. On-going environmental remediation and characterization activities have evolved to improve the performance of groundwater monitoring wells, ensuring they yield representative groundwater and core samples. Groundwater occurs as canyon-fill alluvial water, perched intermediate-depth water associated with the basalts and other local-scale perching horizons, and as a regional-scale aquifer that is the main water supply for the area. Anthropogenic contaminants are found in all three groundwater zones. The regional aquifer below LANL is nominally 1000-1400 feet below ground surface (bgs). Monitoring wells at the Laboratory, are therefore, among the deepest routinely installed in the US. Logistical challenges at LANL include steep terrain, widespread cultural sites, and threatened and endangered species habitat restrictions. Additionally, the cost of installing monitoring wells is high because of groundwater depths and the variable drilling techniques that are required to characterize and install wells in a complex geologic environment. A major challenge to the characterization approach at LANL has been to meet sample quality requirements while minimizing drilling costs. In the late 1990's threaded drill-casing was advanced by rotary drilling methods using compressed air as the only circulation fluid. This resulted in relatively low impacts to the geochemical integrity of rock and groundwater samples, but led to frequent episodes of stuck casing. Additional fluids were introduced to the air-rotary program, improving penetration rates, but also resulting in non-representative groundwater samples due to long-term impacts to the aquifer in the vicinity of the well screen. Mud-rotary drilling was incorporated into the program, but it also had significant long-term impacts on groundwater quality, and often inhibited identification of perched groundwater zones and accurate identification of the water table in the regional aquifer. Current drilling practice makes limited use of fluid additives (e.g. foam) to supplement air circulation methods in order to advance drill casing through the vadose zone. Starting about 100 feet above the water table, drill casing is advanced to the target horizon in the regional aquifer using only air and municipal water (when needed) for circulation. Additionally, the threaded drill casing couplings have been replaced with welded connections, significantly reducing the occurrence of stuck casing. Careful observations of water production and water levels during drilling, coupled with the ability to retract casing for video and geophysical logs, now allows for robust characterization of perched groundwater systems and accurate definition of the top of regional saturation. Recent advances in the drilling program include combined dual-rotary/sonic drilling to depths up to 1150 feet and dual-rotary casing advance at angles up to 25 deg. from vertical. The former allows for core collection with minimal fluid addition, thereby preserving pore water chemistry. The latter method allows for well installation into portions of the aquifer otherwise inaccessible due to constraints at the surface (e.g., cultural sites or topography). Similar to the development of drilling techniques, both well design and groundwater sampling system design have evolved throughout the program. Current well design emphasizes a minimal annulus of filter pack, maximum screen slot size based on formation sieve analyses, and thorough well development to mitigate formation damage due to drilling. To sample groundwater at multiple levels in the same well, earlier sampling systems included Westbay{sup R} multi-zone systems and dual-screen single-pump and dual-screen dual-pump systems designed in collaboration with Baski, Inc. These multi-zone sampling systems proved difficult to develop and maintain given the depth to groundwater, and most current well installations are configured with a single screen and a submersible pump. The ability to purge from a well-developed single-screen well has substantially improved the collection of representative groundwater samples. A carefully crafted groundwater characterization program in target aquifers at depths ranging from 1,000 to 1,400 feet bgs within a complex geologic setting, while also retaining the ability to collect representative core and groundwater samples, may have application at a number of sites throughout the environmental industry. Optimizing each well to meet program objectives while reducing total project costs benefits all environmental investigations. (authors)

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

Scaling Input Distributions for Probabilistic Models - 19472

It is important to address scaling issues when building models to support risk-based decisions for radioactive waste disposal sites. Currently, probabilistic risk assessment and probabilistic performance assessment models are run for multiple realizations where each realization uses one value from each input distribution over the entire modeled time. The input distributions represent the current state of knowledge in the mean of the input variable over the spatial region of interest and over the time period of the model run. This paper focuses on the time factor in input distribution scaling because it presents more serious conceptual and distributional challenges than the spatial factor. One of the challenges in temporal scaling is the effect of non-linearity on propagation of distribution uncertainty through to uncertainty in results. Radioactive waste disposal models tend to have both linear and non-linear components in their models. Dose assessment calculations are typically represented by linear models, while contaminant transport models have linear and non-linear components. In principle, temporal scaling of input parameter distributions for a linear model is straightforward, but scaling of input distributions for non-linear models is not. This paper uses examples to investigate the effects of input distribution temporal scaling on simple linear and non-linear models. Different scaling scenarios were investigated by comparing three analysis tools: analytical mathematical representations, statistical software programming, and GoldSim Monte Carlo simulation software. R statistical software is an open source statistical programming language often used for statistical analyses and graphical representation of data. GoldSim is a robust probabilistic dynamic modeling platform used in performance assessment modeling of radioactive waste sites. A simple linear model, a quadratic model and a multiplicative model were the cases explored for each of these tools. For the three analysis tools considered and with each of the scaling cases, 1,000 realizations were implemented from distributions of the random variables of interest. This is a sufficient number of realizations for simple models to capture the model space. An annual temporal resolution of a 100-year simulation period was chosen for all models and analysis tools. The effects of temporal scaling were demonstrated by investigating two different options for time scaling. Option 1 was set up so that the input distributions represented the expected value on an annual time scale, and a new value from that distribution was chosen every year (at every time step which equaled a year) for 100 years. In Option 2, the input distributions were scaled to represent the expected mean value of the input variable over the 100 years of simulation, and one value from that distribution was chosen for the entire elapsed time. Results show agreement across the three analysis tools. As expected, the linear models have identical results for the two different scaling Options, and the non-linear models diverge in the results for different time scale implementation. The implications of this work point to the need for careful consideration of input distribution temporal scaling. (authors)

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

Nevada Higher Education Benefits from US Department of Energy (DOE) Environmental Management (EM) Nevada Program's Transfer of Geologic Samples - 20442

The Nevada National Security Site (NNSS), formerly the Nevada Test Site, was the location of 100 historic atmospheric and 828 historic underground nuclear tests from 1951 to 1992. Related to this historic nuclear testing, geologic and hydrologic studies of the site were conducted utilizing the skills and the expertise of the U.S. Geological Survey (USGS) and national laboratories (Los Alamos National Laboratory, Lawrence Livermore National Laboratory), who were principal in leading the development of weapons and were responsible for specific underground testing programs. The subsurface samples and technical/scientific data associated with NNSS geologic studies (past and present) are preserved and stored at the USGS Mercury Core Library and Data Center located at the NNSS. Currently, the facility stores over 2,000,000 linear feet of cores and cuttings from more than 2,600 drill holes that can be accessed for study. Most of the samples and historic geologic work was focused directly on subsurface geologic settings that relate directly to historic underground nuclear testing. Since 1992, the United States has observed a unilateral moratorium on full-scale nuclear testing, and the U.S. Department of Energy (DOE) Environmental Management (EM) Nevada Program is now responsible for hydrogeologic characterization of the potential impacts to the natural groundwater systems that may have resulted from historic underground nuclear tests. As part of this effort, the EM Nevada Program Underground Test Area (UGTA) Activity has drilled and completed over 58 deep (2,000 - 7,000 ft.) characterization wells, totaling in excess of 170,000 linear feet of cuttings and core samples. Cutting samples were collected as triplicate samples for each respective depth interval, to account for potential later nondestructive/destructive analysis and to preserve samples for regulatory purposes. Recently, it was recognized by the EM Nevada Program that opportunities may exist to reduce the cost and floor space required for the storage of geologic samples at the Mercury Core Library without impacting the integrity and representative nature of the samples necessary for project execution. An initiative was sponsored by EM Nevada Program to evaluate several options: 1) disposal of a portion of cuttings and cores in a land fill setting; 2) reduction in sample volume through skeletonizing core and cuttings from wells to a representative but much smaller number of samples; and 3) solicit potential interest in the academic community where other geoscientists could freely access the samples for studies. The Nevada state university system through the Nevada Bureau of Mines and Geology responded positively to the opportunity to receive these NNSS samples from the EM Nevada Program. In June 2019, approximately 17,000 geologic samples, representing greater than 170,000 linear feet of drilling and weighing over 20,000 pounds, were shipped from the NNSS to the Great Basin Science Sample and Records Library located in Reno, Nevada. The benefits resulting from this transfer were realized by the EM Nevada Program in terms of cost and space savings for geologic sample storage at the NNSS. The transfer did not impact the EM Nevada Program mission as access to representative geologic samples for regulatory and scientific purposes was preserved. Moreover, the Nevada Bureau of Mines and Geology acquired a significant resource of subsurface geologic samples and supporting technical data to support academic and scientific studies in a complex volcanic setting in southwest Nevada. (authors)

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

Development of a Conditioned System-Level Groundwater Model to Evaluate Long-Term Groundwater Impacts within a Performance Assessment - 20115

A preliminary performance assessment (PA) of single-shell tank Waste Management Area (WMA) A-AX located at the U.S. Department of Energy (DOE)'s Hanford Site in southeastern Washington is being conducted to satisfy the requirements of DOE Order 435.1 [1] as it relates to closure of the single-shell radioactive waste tanks in the WMA A-AX tank farms. A PA assesses the fate, transport, and impacts of radionuclides within a low-level radioactive waste disposal facility in its assumed closure configuration and the subsequent potential doses to humans over a 1,000-year compliance period and a 10,000-year performance evaluation period. The WMA A-AX preliminary PA evaluation is structured around the complementary use of process-level and system-level models to calculate the facility performance against established DOE Order 435.1 [1] performance objectives. Process-level models are those that represent a detailed phenomenological representation of processes of concern in the PA. Process models typically only represent one or a few of the components of the PA, such as groundwater flow and transport, and must be integrated with other modeling elements to perform PA calculations. System-level models are those that are abstracted from the process models, retaining the essential features of the process model, while allowing integration of all aspects of the PA in a single modeling framework. System-level models are often characterized by coarser numerical discretization, lower dimensionality, or other similar simplifications compared to the process-level model. Traditionally, a three-dimensional (3-D) process-level model is utilized primarily to evaluate the long-term impact on groundwater and the potential doses to individuals who consume contaminated groundwater. System-level models are also utilized to evaluate the groundwater pathway in PAs. These models typically have reduced dimensionality (1-D) and are conditioned utilizing flow fields (Darcy fluxes) and moisture content distributions that are abstracted from the process level models. The abstraction approach assures that the flow field in both models is consistent for a specific set of input parameters for flow, differing only in the discretization and dimensionality of the two models. The preliminary WMA A-AX PA incorporates a detailed representation of the geological system and hydraulic properties within the 3-D model STOMP{sup C} numerical code so that the effects of relevant features and processes on water flow and radionuclide transport in the subsurface can be evaluated. The complementary system-level model is developed utilizing the GoldSim{sup C} code to implement a simplified, 1-D equivalent model to represent the groundwater pathway. Contaminant transport through the vadose zone and unconfined aquifer for Tc-99 and I-129 were evaluated in each of the models. Adjustments to the saturated portion of the GoldSim{sup C} 1-D model were required to mimic the dispersion effect captured with the 3-D STOMP{sup C} model. Once this conditioning was conducted, highly similar results for the transport of Tc-99 and I-129 were achieved at the point of calculation, located 100 m downgradient from the WMA A-AX fenceline. These radionuclides represent elements that are regarded as primary dose drivers in the PA groundwater pathway analysis. The high degree of conformance between the two models suggests that the use of the equivalent 1-D system model is suitable for evaluating the full suite of radionuclides that will be released from the tank sources within WMA A-AX over the 1,000-year compliance period and over the 10,000-year evaluation period. (authors)

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

A qualitative assessment of current CCF guidance based on a review of safety system digital implementation changes with evolving technology - 148

This paper is a summary of a study performed in support of the US Nuclear Regulatory Commission (NRC) to evaluate current policy on software common-cause failure (CCF). The study reported in this paper is one of several that contributed to the technical basis to inform NRC staff. The study first reviewed policies and assessment guidance as discussed in the Staff Requirements Memorandum to the Secretary of the Commission, Office of the NRC (SECY) 93-087, Policy, Technical, and Licensing Issues Pertaining to Evolutionary and Advanced Light Water Reactor (ALWR) Designs, and Branch Technical Position (BTP) 7-19, Guidance on Evaluation of Defense-in-Depth and Diversity in Digital Computer-Based Instrumentation and Control Systems, as well as Chapter 7, 'Instrumentation and Controls', in NUREG-0800, Standard Review Plan for Review of Safety Analysis Reports for Nuclear Power Plants. The study then examined instrumentation and controls (I and C) technology implementations in nuclear power plants in the light of current CCF guidance. The intent was to assess whether the current position on CCF remains appropriate given the evolutions in digital safety system implementations and, if gaps in the guidance were found, to provide recommendations as to how these gaps could be closed. The methodology adopted was to review the vendors' technology and software implementation processes for digital safety systems as technology evolved. The following three representative safety systems were selected to provide illustrative examples: - Eagle 21 was selected to represent vintage microprocessor-based technology used from the 1980's to mid-1990's, - TELEPERM XS (TXS) was selected to represent second generation microprocessor-based technology used from the mid-1990's to the 2000's, and - The advanced logic system (ALS) was selected to represent the latest trend of using technology based on field programmable gate arrays (FPGAs). Technology implementations were reviewed in light of the basic premise of the NRC in BTP 7-19 (Revision 6), that software-based or software-logic-based digital system development errors are a credible source of CCF and therefore are susceptible to CCF because identical copies of the software-based logic and architecture are present in redundant divisions of safety-related systems. BTP 7-19 categorizes firmware and logic developed from software-based development systems all under software. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Assessment of Multi-Step CADIS Parameter Sensitivity to Neutron Flux Spectrum in Shutdown Dose Rate Calculations

Shutdown dose rate (SDDR) analysis requires (a) a neutron transport calculation to estimate neutron flux fields, (b) an activation calculation to compute radionuclide inventories and associated photon sources, and (c) a photon transport calculation to estimate final dose rate. In some applications, accurate full-scale Monte Carlo (MC) SDDR simulations are needed for very large systems with massive amounts of shielding materials. For example, SDDR assessments are required throughout the biological shield (bio-shield) of the ITER experimental facility to evaluate the required waiting period after the shutdown of ITER and to identify the locations where human accessibility should be prohibited. The bio-shield is a large cylindrical concrete structure (30 m tall and 34 m in diameter) surrounding the very complex tokamak machine. Determining the effects on SDDR of important factors, such as the cross talk (interactions) between the different ports of ITER, is possible only through full-scale simulations that involve all the complex inner details of the ITER tokamak machine. These simulations are impractical because calculation of space- and energy-dependent neutron fluxes throughout the structural materials is needed to estimate the distribution of radioisotopes causing the SDDR. Determining a good importance function that can be used in biasing a neutron MC calculation for SDDR is complicated because it is difficult to explicitly express the SDDR response function, which depends on subsequent computational steps. The Multi-Step Consistent Adjoint Driven Importance Sampling (MS-CADIS) hybrid MC/deterministic method was developed to speed up the SDDR MC neutron transport calculation using an importance function that represents the neutron importance to the final SDDR. The MS-CADIS method uses the CADIS method, which has been Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for the United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan). successfully used for more than a decade in radiation shielding applications to generate consistent source biasing and weight window (WW) variance reduction parameters to accelerate continuous-energy MC simulations. However, because MS-CADIS focuses on multi-step shielding calculations such as SDDR calculations, it develops an importance function for the initial radiation transport calculation (e.g., the neutron calculations in SDDR simulations) that represents the importance of particles to the final response of the overall simulation. To develop this importance function, the MS-CADIS method uses an adjoint calculation with a specific neutron adjoint source that represents the neutron's contribution to the final SDDR. Brute force calculation of the MS-CADIS neutron adjoint source requires a deterministic neutron transport calculation coupled with activation calculations in each element of the deterministic mesh. These calculations are needed to determine the MS-CADIS parameter for the neutron interactions that produce radioisotopes contributing to the SDDR. The MS-CADIS parameter is the microscopic cross section of the radioisotope-producing interaction, multiplied by the mass of this radioisotope existing at the end of the irradiation and decay scenario, and divided by the steady-state rate of the interaction. Because the activation calculations are performed for each radionuclide in each mesh element, implementation of the MS-CADIS method can be computationally expensive. However, these neutron transport and activation calculations are not needed if we assume that the radioisotopes are created due to single neutron interactions and their concentrations do not change as a result of burnup by neutron irradiation. The use of this approximation, which renders the MS-CADIS parameter independent of the neutron flux, seems reasonable in calculating the MS-CADIS importance function because this function is needed only to speed up the MC calculation, not to calculate the final SDDR. In this work, we investigate the validity of this approximation by investigating the behavior of the MS-CADIS parameter as a function of the neutron flux spectra in the Fusion Evaluated Nuclear Data Library (FENDL) Neutronics and Shielding Benchmark problem that contains representative isotopes and flux environments of fusion energy systems. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
NASA NTRS2020

An Algorithm for Characterization of Fiber Aggregation in Composite Microstructures

Composite structures are susceptible to localized flaws, or variability, that drive global failure. To capture this variance, a multiscale model needs to be introduced that not only accurately represents the statistical nature of the composite microstructure but is also efficient. At the microscale, fiber aggregation creates local stress concentrations, where failure is likely to occur sooner than expected. A method of rapid microstructure generation, in conjunction with cluster characterization, is examined to develop accurate reproductions of 2D composite cross-sections. Parameterization of shape and size of fiber clusters is used to characterize representative volume elements. The discrete element method is used to generate pseudo-microstructures. The clustering parameters from the pseudo- and actual microstructure arrangements (obtained from micrographs) can be compared to determine the validity of the representative volume element generated with the discrete element method. Results show a promising approach to evaluating randomness of fiber distributions and how to accurately recreate microstructures for strength analysis.

Carbon fiber↗
DOE OSTI2016

The impact of equilibrating hemispheric albedos on tropical performance in the HadGEM2-ES coupled climate model

Despite the fact that the southern hemisphere contains a far greater proportion of dark ocean than the northern hemisphere, the total amount of sunlight reflected from the hemispheres is equal. However, the majority of climate models do not adequately represent this equivalence. Here we examine the impact of equilibrating hemispheric albedos by various idealised methods in a comprehensive coupled climate model and find significant improvements in what have been considered longstanding and apparently intractable model biases. Monsoon precipitation biases almost vanish over all continental land areas, the penetration of monsoon rainfall across the Sahel and the west African monsoon “jump” become well represented, and indicators of hurricane frequency are significantly improved. The results appear not to be model specific, implying that hemispheric albedo equivalence may provide a fundamental constraint for climate models that must be satisfied if the dynamics driving these processes, in particular the strength of the Hadley cell, are to be adequately represented. Cross-equatorial energy transport is implicated as a crucial component that must be accurately modelled in coupled general circulation models. The results also suggest that the commonly used practice of prescribing sea-surface temperatures in models provides a less accurate represention of precipitation than constraining the hemispheric albedos.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Modeling stream temperature in the Anthropocene: An earth system modeling approach

A new large-scale stream temperature model has been developed within the Community Earth System Model (CESM) framework. The model is coupled with the Model for Scale Adaptive River Transport (MOSART) that represents river routing and a water management model (WM) that represents the effects of reservoir operations and water withdrawals on flow regulation. The coupled models allow the impacts of reservoir operations and withdrawals on stream temperature to be explicitly represented in a physically based and consistent way. The models have been applied to the Contiguous United States driven by observed meteorological forcing. It is shown that the model is capable of reproducing stream temperature spatiotemporal variation satisfactorily by comparison against the observed streamflow from over 320 USGS stations. Including water management in the models improves the agreement between the simulated and observed streamflow at a large number of stream gauge stations. Both climate and water management are found to have important influence on the spatiotemporal patterns of stream temperature. More interestingly, it is quantitatively estimated that reservoir operation could cool down stream temperature in the summer low-flow season (August – October) by as much as 1~2oC over many places, as water management generally mitigates low flow, which has important implications to aquatic ecosystems. In conclusion, sensitivity of the simulated stream temperature to input data and reservoir operation rules used in the WM model motivates future directions to address some limitations in the current modeling framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Influence of small-scale fluvial architecture on CO 2 trapping processes in deep brine reservoirs

A number of important candidate CO 2 reservoirs exhibit sedimentary architecture reflecting fluvial deposition. Recent studies have led to new conceptual and quantitative models for sedimentary architecture in fluvial deposits over a range of scales that are relevant to CO 2 injection and storage. We used a geocellular modeling approach to represent this multiscaled and hierarchical sedimentary architecture. With this model, we investigated the dynamics of CO 2 plumes, during and after injection, in such reservoirs. The physical mechanism of CO 2 trapping by capillary trapping incorporates a number of related processes, i.e., residual trapping, trapping due to hysteresis of the relative permeability, and trapping due to hysteresis of the capillary pressure. Additionally, CO 2 may be trapped due to differences in capillary entry pressure for different textural sedimentary facies (e.g., coarser-grained versus finer-grained cross sets). The amount of CO 2 trapped by these processes depends upon a complex system of nonlinear and hysteretic characteristic relationships including how relative permeability and capillary pressure vary with brine and CO 2 saturation. The results strongly suggest that representing small-scale features (decimeter to meter), including their organization within a hierarchy of larger-scale features, and representing their differences in characteristic relationships can all be critical to understanding trapping processes in some important candidate CO 2 reservoirs.

58 GEOSCIENCES↗
DOE OSTI2017

The coastal streamflow flux in the Regional Arctic System Model

The coastal streamflow flux from the Arctic drainage basin is an important driver of dynamics in the coupled ice-ocean system. Comprising more than one third of the total freshwater flux into the Arctic Ocean, streamflow is a key component of the regional and global freshwater cycle. To better represent the coupling of the streamflow flux to the ocean, we have developed and applied the RVIC streamflow routing model within the Regional Arctic System Model (RASM). The RASM is a high-resolution regional Earth System Model whose domain includes all of the Arctic drainage basin. Here in this paper, we introduce the RVIC streamflow routing model, detailing its application within RASM and its advancements in terms of representing high-resolution streamflow processes. We evaluate model simulated streamflow relative to in situ observations and demonstrate a method for improving model performance using a simple optimization procedure. We also present a new, spatially and temporally consistent, high-resolution data set of coastal freshwater fluxes for the Arctic drainage basin and surrounding areas that is based on a fully coupled RASM simulation and intended for use in Arctic Ocean modeling applications. This data set is evaluated relative to other coastal streamflow data sets commonly used by the ocean modeling community. We demonstrate that the RASM-simulated streamflow flux better represents the annual cycle than existing data sets, especially in ungauged areas. Finally, we assess the impact that streamflow has on the coupled ice-ocean system, finding that the presence of streamflow leads to reduced sea surface salinity, increased sea surface temperatures, and decreased sea ice thickness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Observational needs for estimating Alaskan soil carbon stocks under current and future climate

Abstract Representing land surface spatial heterogeneity when designing observation networks is a critical scientific challenge. Here we present a geospatial approach that utilizes the multivariate spatial heterogeneity of soil‐forming factors—namely, climate, topography, land cover types, and surficial geology—to identify observation sites to improve soil organic carbon (SOC) stock estimates across the State of Alaska, USA. Standard deviations in existing SOC samples indicated that 657, 870, and 906 randomly distributed pedons would be required to quantify the average SOC stocks for 0–1 m, 0–2 m, and whole‐profile depths, respectively, at a confidence interval of 5 kg C m −2 . Using the spatial correlation range of existing SOC samples, we identified that 309, 446, and 484 new observation sites are needed to estimate current SOC stocks to 1 m, 2 m, and whole‐profile depths, respectively. We also investigated whether the identified sites might change under future climate by using eight decadal (2020–2099) projections of precipitation, temperature, and length of growing season for three representative concentration pathway (RCP 4.5, 6.0, and 8.5) scenarios of the Intergovernmental Panel on Climate Change. These analyses determined that 12 to 41 additional sites (906 + 12 to 41; depending upon the emission scenarios) would be needed to capture the impact of future climate on Alaskan whole‐profile SOC stocks by 2100. The identified observation sites represent spatially distributed locations across Alaska that captures the multivariate heterogeneity of soil‐forming factors under current and future climatic conditions. This information is needed for designing monitoring networks and benchmarking of Earth system model results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Closing the scale gap between land surface parameterizations and GCMs with a new scheme, SiB3-Bins: SOIL MOISTURE SCALE GAP

The interaction of land with the atmosphere is sensitive to soil moisture (W). Evapotranspiration (ET) reacts to soil moisture in a nonlinear way, f(W), as soils dry from saturation to wilt point. This nonlinear behavior and the fact that soil moisture varies on scales as small as 1–10 m in nature, while numerical general circulation models (GCMs) have grid cell sizes on the order of 1 to 100s of kilometers, makes the calculation of grid cell-average ET problematic. It is impractical to simulate the land in GCMs on the small scales seen in nature, so techniques have been developed to represent subgrid scale heterogeneity, including: (1) statistical-dynamical representations of grid subelements of varying wetness, (2) relaxation of f(W), (3) moderating f(W) with approximations of catchment hydrology, (4) “tiling” the landscape into vegetation types, and (5) hyperresolution. Here we present an alternative method for representing subgrid variability in W, one proven in a conceptual framework where landscape-scale W is represented as a series of “Bins” of increasing wetness from dry to saturated. The grid cell-level f(W) is defined by the integral of the fractional area of the wetness bins and the value of f(W) associated with each. This approach accounts for the spatiotemporal dynamics of W. We implemented this approach in the SiB3 land surface parameterization and then evaluated its performance against a control, which assumes a horizontally uniform field of W. We demonstrate that the Bins method, with a physical basis, attenuates unrealistic jumps in model state and ET seen in the control runs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A new process sensitivity index to identify important system processes under process model and parametric uncertainty

Hydrological models are always composed of multiple components that represent processes key to intended model applications. When a process can be simulated by multiple conceptual-mathematical models (process models), model uncertainty in representing the process arises. While global sensitivity analysis methods have been widely used for identifying important processes in hydrologic modeling, the existing methods consider only parametric uncertainty but ignore the model uncertainty for process representation. To address this problem, this study develops a new method to probe multimodel process sensitivity by integrating the model averaging methods into the framework of variance-based global sensitivity analysis, given that the model averaging methods quantify both parametric and model uncertainty. A new process sensitivity index is derived as a metric of relative process importance, and the index includes variance in model outputs caused by uncertainty in both process models and model parameters. For demonstration, the new index is used to evaluate the processes of recharge and geology in a synthetic study of groundwater reactive transport modeling. The recharge process is simulated by two models that converting precipitation to recharge, and the geology process is also simulated by two models of different parameterizations of hydraulic conductivity; each process model has its own random parameters. The new process sensitivity index is mathematically general, and can be applied to a wide range of problems in hydrology and beyond.

Dai, Heng↗
DOE OSTI2018

Investigation of Turbulent Entrainment-Mixing Processes With a New Particle-Resolved Direct Numerical Simulation Model

Here, a new particle-resolved three dimensional direct numerical simulation (DNS) model is developed that combines Lagrangian droplet tracking with the Eulerian field representation of turbulence near the Kolmogorov microscale. Six numerical experiments are performed to investigate the processes of entrainment of clear air and subsequent mixing with cloudy air and their interactions with cloud microphysics. The experiments are designed to represent different combinations of three configurations of initial cloudy area and two turbulence modes (decaying and forced turbulence). Five existing measures of microphysical homogeneous mixing degree are examined, modified, and compared in terms of their ability as a unifying measure to represent the effect of various entrainment-mixing mechanisms on cloud microphysics. Also examined and compared are the conventional Damköhler number and transition scale number as a dynamical measure of different mixing mechanisms. Relationships between the various microphysical measures and dynamical measures are investigated in search for a unified parameterization of entrainment-mixing processes. The results show that even with the same cloud water fraction, the thermodynamic and microphysical properties are different, especially for the decaying cases. Further analysis confirms that despite the detailed differences in cloud properties among the six simulation scenarios, the variety of turbulent entrainment-mixing mechanisms can be reasonably represented with power-law relationships between the microphysical homogeneous mixing degrees and the dynamical measures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Potassium Metal Batteries: Stable Potassium Metal Anodes with an All-Aluminum Current Collector through Improved Electrolyte Wetting (Adv. Mater. 49/2020)

An explorer is kayaking in uncharted territory, representing the exciting scientific pursuit of emerging metal batteries (MBs). Water is flowing smoothly over the jagged rocks, representing dendrite-free metal anodes as well as enhanced electrolyte wetting, both due to geometrical roughness. The electric eels underneath represent potassium-based energy storage, as they utilize potassium and sodium ions for electricity. In article number 2002908, Pengcheng Liu, David Mitlin, and co-workers present a route via improved electrolyte wetting toward stable potassium metal batteries with an all-aluminum current collector.

25 ENERGY STORAGE↗