SEARCH · Search NASA

Results for “sampling variability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2016

The Two-Column Aerosol Project: Phase I-Overview and impact of elevated aerosol layers on aerosol optical depth

The Two-Column Aerosol Project (TCAP), conducted from June 2012 through June 2013, was a unique study designed to provide a comprehensive data set that can be used to investigate a number of important climate science questions, including those related to aerosol mixing state and aerosol radiative forcing. The study was designed to sample the atmosphere between and within two atmospheric columns; one fixed near the coast of North America (over Cape Cod, MA) and a second moveable column over the Atlantic Ocean several hundred kilometers from the coast. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) was deployed at the base of the Cape Cod column, and the ARM Aerial Facility was utilized for the summer and winter intensive observation periods. One important finding from TCAP is that four of six nearly cloud-free flight days had aerosol layers aloft in both the Cape Cod and maritime columns that were detected using the nadir pointing second-generation NASA high-spectral resolution lidar (HSRL-2). In addition, these layers contributed up to 60% of the total observed aerosol optical depth (AOD). Many of these layers were also intercepted by the aircraft configured for in situ sampling, and the aerosol in the layers was found to have increased amounts of biomass burning material and nitrate compared to aerosol found near the surface. Lastly, while there was a great deal of spatial and day-to-day variability in the aerosol chemical composition and optical properties, no systematic differences between the two columns were observed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Towards retrieving critical relative humidity from ground-based remote sensing observations

Nearly all parameterisations of large-scale cloud require the specification of the critical relative humidity (RHcrit). This is the gridbox-mean relative humidity at which the subgrid fluctuations in temperature and water vapour become so large that part of a subsaturated gridbox becomes saturated and cloud starts to form. Until recently, the lack of high-resolution observations of temperature and moisture variability has hindered a reasonable estimate of the RHcrit from observations. However, with the advent of ground-based measurements from Raman lidar, it becomes possible to obtain long records of temperature and moisture (co-)variances with sub-minute sample rates. Lidar observations are inherently noisy and any analysis of higher-order moments will be very dependent on the ability to quantify and remove this noise. We present an exporatory study aimed at understanding whether current noise levels of lidar-retrieved temperature and water vapour are sufficient to obtain a reasonable estimate of the RHcrit. We show that vertical profiles of RHcrit can be derived for a gridbox length of up to about 30 km (120) with an uncertainty of about 4 % (2 %). RHcrit tends to be smallest near the scale height and seems to be fairly insensitive to the horizontal grid spacing at the scales investigated here (30 - 120 km). However, larger sensitivity was found to the vertical grid spacing. As the grid spacing decreases from 400 to 100 m, RHcrit is observed to increase by about 6 %, which is more than the uncertainty in the RHcrit retrievals.

Van Weverberg, Kwinten↗

DOE OSTI2021

Properties and reaction mechanism of phosphoric acid activated metakaolin geopolymer at varied curing temperatures

Reaction mechanism of phosphoric activated metakaolin is studied by techniques of calorimetry, ICP-OES, FTIR and NMR taking curing temperature and [H{sub 3}PO{sub 4}] as variables. A two-stage curing method including pre-curing at 40 °C for 24 h and second curing at 60/80 °C for another 24 h effectively mitigates thermal-cracking of the obtained silico-aluminate phosphate (SAP) samples, whose maximum compressive strength reaches 120 MPa. Dealumination of metakaolin in H{sub 3}PO{sub 4} solution produces low aluminate-containing silicate units and subsequently condensation of silicate tetrahedrons occurs. Higher temperature and H{sub 3}PO{sub 4} concentration accelerate the two processes. After complete reaction, A{sub VI} structure from metakaolin disappears and the dissolved Al preferably forms Al{sub VI}-O-P units. During reaction, metastable P-O-P intermediate appears in SAPs with higher P/Al ratio and finally transforms into Al-O-P structure upon PO{sub 4}{sup 3−} consumption. Phosphate units as P(OAl){sub x}(H{sub 2}O){sub 4-x} are connected with Si or Al structure units. Major structural units of the SAPs include Al-O-P, Si-O-P, Si-O-Si, and Si-O-Al.

36 MATERIALS SCIENCE↗

DOE OSTI2019

A comprehensive modeling procedure for estimating statistical properties of forced ignition

In this work, a comprehensive modeling procedure for estimating the probability of ignition with application to high altitude relights of aircraft combustors is developed. In these configurations, an ignitor is used to introduce high-enthalpy discharge into a fuel-laden but stratified flow. Because of the inherent variabilities inflow conditions, kernel discharge process, and the chaotic turbulent flow, ignition cannot be described deterministically, but only as a probabilistic measure. The proposed modeling framework consists of three components. The turbulent flow is represented using the large eddy simulation (LES) framework. The ignition process is modeled using a manifold approach, where the initial kernel evolution is represented using a homogeneous reactor while the latter part of the evolution is represented as a competition between diffusion and chemical reactions using a flamelet-type mapping. A combined lookup table that can track the evolution of the kernel through these two distinct reaction stages is developed. The table lookup variables are solved within the LES framework. The resulting simulation tool is then embedded within an uncertainty quantification approach, where variations in the turbulent flow, as well as operating and kernel properties, are simulated using a Monte-Carlo-based sampling approach. Techniques to reduce computational cost are used to obtain a robust, numerically accurate, and physically representative model for engine relight. The method is validated using experimental data for ignition of methane/air mixtures. Due to the comprehensive nature of the modeling procedure, it is found that the simulation tool reproduces experimentally observed ignition probabilities over a wide range of operating conditions.

42 ENGINEERING↗

DOE OSTI2020

A Satellite-Based Estimate of Convective Vertical Velocity and Convective Mass Flux: Global Survey and Comparison with Radar Wind Profiler Observations

Convective vertical velocity (w c ) and convective mass flux (M c ) lie at the heart of GCM cumulus parameterizations, but few observations of these critical parameters are available. In this paper, we develop and evaluate a novel, satellite-based method for estimating profiles of w c and M c . Here, comparisons with collocated ground-based radar wind profiler (RWP) observations show that satellite estimated median w c is slightly greater than the RWP estimates, but they show solid agreement when compared at the 95th percentiles (intense updrafts). RWP-derived and satellite estimated M c are broadly comparable in the lower and middle troposphere, with some differences in the upper troposphere due to differences in convective core sampling. A k-means cluster analysis of multiple years of w c data shows that convective characteristics are distinctly different among extratropical convection, tropical land convection, and tropical oceanic convection. Tropical land convection is significantly more intense and more variable than the oceanic counterpart.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Deep Learning Reconstruction of Daily Soil CO 2 Efflux Reveals Biogeochemical Insights and Reduces Annual Estimate Uncertainty Despite Limited Daily Predictability

Soil CO 2 efflux is commonly measured monthly or seasonally, leaving daily dynamics poorly resolved and contributing to global estimation uncertainty. We trained a single Long Short-Term Memory (LSTM) model to predict daily soil CO 2 efflux across 82 globally distributed sites in COSORE, with 0.2%–46.9% daily data coverage from 2003 to 2020. Despite using far fewer sites than are typically used to train a single deep learning model, with observations biased toward temperate mesic sites, the LSTM model performed well at approximately one-third of sites, reconstructed nearly 2 decades of daily efflux, and outperformed commonly used approaches for estimating daily efflux when applied to the same data set. Performance was weakest at pronounced peaks and troughs and at non-temperate sites with <1.5 years of observations and irregular data patterns. Nevertheless, annual efflux from reconstructed daily data had <40% error even at underperforming sites, substantially improving estimates derived from monthly and seasonal sampling (maximum errors of 95% and 136%, respectively). Temperature sensitivity (Q 10 ) estimated from reconstructed daily predictions closely matched estimates from daily observations, whereas Q 10 values derived from monthly or seasonal observations deviated substantially, suggesting that coarse temporal sampling may contribute to uncertainty in reported Q 10 values. Consistent daily reconstructions further enabled trend analyses for well-performing, predominantly temperate sites and showed increasing soil CO 2 efflux at most sites from 2003 to 2020, with more variable summer trends. Despite limitations, these results demonstrate the potential of LSTM models to reconstruct daily soil CO 2 efflux and reduce estimation uncertainties from sparse observations.

Smykalov, Valerie [Pennsylvania State University, ↗

DOE OSTI2018

SSAGES: Software Suite for Advanced General Ensemble Simulations

Molecular simulation has emerged as an essential tool for modern-day research, but obtaining proper results and making reliable conclusions from simulations requires adequate sampling of the system under consideration. To this end, a variety of methods exist in the literature that can enhance sampling considerably, and increasingly sophisticated, effective algorithms continue to be developed at a rapid pace. Implementation of these techniques, however, can be challenging for experts and non-experts alike. There is a clear need for software that provides rapid, reliable, and easy access to a wide range of advanced sampling methods, and that facilitates implementation of new techniques as they emerge. Here we present SSAGES, a publicly available Software Suite for Advanced General Ensemble Simulations designed to interface with multiple widely used molecular dynamics simulations packages. SSAGES allows facile application of a variety of enhanced sampling techniques—including adaptive biasing force, string methods, and forward flux sampling—that extract meaningful free energy and transition path data from all-atom and coarse grained simulations. A noteworthy feature of SSAGES is a user-friendly framework that facilitates further development and implementation of new methods and collective variables. In this study, the use of SSAGES is illustrated in the context of simple representative applications involving distinct methods and different collective variables that are available in the current release of the suite.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Precise measurement of the tt¯ production cross-section and lepton differential distributions in eμ dilepton events from s=13TeVpp collisions with the ATLAS detector

The inclusive top quark pair (tt¯$$t\bar{t}$$) cross-section σtt¯$$\sigma _{t\bar{t}}$$ has been measured in proton–proton collisions at s=13TeV$$\sqrt{s}=13\,\text {TeV}$$, using 140fb-1$$140\,{\text {fb}^{-1}} $$ of data collected by the ATLAS experiment at the Large Hadron Collider. Using events with an opposite-charge eμ$$e\mu $$ pair and b-tagged jets, the cross-section is measured to be: σtt¯=829.3±1.3(stat)±8.0(syst)±7.3(lumi)±1.9(beam)pb,$$\begin{aligned} \sigma _{t\bar{t}} & = 829.3\pm 1.3\,\mathrm {(stat)}\ \pm 8.0\,\mathrm {(syst)}\ \pm 7.3\,\mathrm {(lumi)}\ \\ & \quad \pm 1.9\,\mathrm {(beam)}\,\textrm{pb}, \end{aligned}$$where the uncertainties reflect the limited size of the data sample, experimental and theoretical systematic effects, the integrated luminosity, and the proton beam energy, giving a total uncertainty of 1.3%. The result is used to determine the top quark pole mass via the dependence of the predicted cross-section on mtpole$${m_{t}^\textrm{pole}}$$, giving mtpole=172.8-1.7+1.5$${m_{t}^\textrm{pole}}=172.8^{+1.5}_{-1.7}$$ GeV$$\text {GeV}$$. The same event sample is used to measure absolute and normalised differential cross-sections for the tt¯→eμνν¯bb¯$$t\bar{t} \rightarrow e\mu u \bar{ u }b\bar{b} $$ process as a function of single-lepton and dilepton kinematic variables. Complementary measurements of eμbb¯$$e\mu b\bar{b} $$ production, treating both tt¯$$t\bar{t}$$ and Wt events as signal, are also provided. Both sets of differential cross-sections are compared to the predictions of various Monte Carlo event generators, demonstrating that the state-of-the-art generators Powheg MiNNLO and Powhegbb4l describe the data better than Powheghvq.

Aad, G↗

DOE OSTI2014

Estimation of k-ε parameters using surrogate models and jet-in-crossflow data

We demonstrate a Bayesian method that can be used to calibrate computationally expensive 3D RANS (Reynolds Av- eraged Navier Stokes) models with complex response surfaces. Such calibrations, conditioned on experimental data, can yield turbulence model parameters as probability density functions (PDF), concisely capturing the uncertainty in the parameter estimates. Methods such as Markov chain Monte Carlo (MCMC) estimate the PDF by sampling, with each sample requiring a run of the RANS model. Consequently a quick-running surrogate is used instead to the RANS simulator. The surrogate can be very difficult to design if the model's response i.e., the dependence of the calibration variable (the observable) on the parameter being estimated is complex. We show how the training data used to construct the surrogate can be employed to isolate a promising and physically realistic part of the parameter space, within which the response is well-behaved and easily modeled. We design a classifier, based on treed linear models, to model the "well-behaved region". This classifier serves as a prior in a Bayesian calibration study aimed at estimating 3 k - ε parameters ( C μ , C ε2 , C ε1 ) from experimental data of a transonic jet-in-crossflow interaction. The robustness of the calibration is investigated by checking its predictions of variables not included in the cal- ibration data. We also check the limit of applicability of the calibration by testing at off-calibration flow regimes. We find that calibration yield turbulence model parameters which predict the flowfield far better than when the nomi- nal values of the parameters are used. Substantial improvements are still obtained when we use the calibrated RANS model to predict jet-in-crossflow at Mach numbers and jet strengths quite different from those used to generate the ex- perimental (calibration) data. Thus the primary reason for poor predictive skill of RANS, when using nominal values of the turbulence model parameters, was parametric uncertainty, which was rectified by calibration. Post-calibration, the dominant contribution to model inaccuraries are due to the structural errors in RANS.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Thermal Characterization and Exhumation of Northwest San Juan Basin Area, NM

The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.

02 PETROLEUM↗

DOE OSTI2017

Testing the Accuracy of Data-driven MHD Simulations of Active Region Evolution

Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of the magnetic field lie at the photosphere, necessitating the development of coronal models which are “data-driven” at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth data set, to supply both the photospheric information and to perform the validation of the data-driven method. We focus our investigation on how the accuracy of the data-driven model depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA’s Solar Dynamics Observatory produces full-disk vector magnetic field measurements at a 12-minute cadence. Using our framework we show that ARs that emerge over 25 hr can be modeled by the data-driving method with only ∼1% error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors, such as reduced resolution, instrument bias, and noise.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Application of High Performance Computing for Simulating Cycle-to-Cycle Variation in Dual-Fuel Combustion Engines

Interest in operational cost reduction is driving engine manufacturers to consider lower-cost fuel substitution in heavy-duty diesel engines. These dual-fuel (DF) engines could be operated either in diesel-only mode or operated with premixed natural gas (NG) ignited by a pilot flame of compression-ignited direct-injected diesel fuel. One promising application is that of large-bore, medium-speed engines such as those used in locomotives. With realistic natural gas substitution levels in the fleet of locomotives currently in service, such fuel substitution could result in billions of dollars of savings annually in the US alone. However, under certain conditions, dual-fuel operation can result in increased cycle-to-cycle variability (CCV) during combustion, resulting in variations in cylinder pressure and work extraction. In certain situations, the CCV of dual-fuel operation can be notably higher than that of diesel-only combustion under similar operating conditions. Excessive CCV can limit the NG substitution rate and operating range of a dual-fuel engine by increasing emissions and reducing engine stability, reliability and fuel efficiency via incomplete natural-gas combustion. Running multiple engine cycles in series to simulate CCV can be quite time consuming. Hence innovative modelling techniques and large computing resources are needed to investigate the factors affecting CCV in dual-fuel engines. This paper discusses the use of the High Performance Computing resource Titan, at the Oak Ridge Leadership Computing Facility at Oak Ridge National Laboratory, to investigate cycle-to-cycle combustion variability of a dual-fuel engine. The CONVERGE CFD software was used to simulate multiple, parallel single cycles of dual-fuel combustion with perturbed operating parameters and boundary conditions. These perturbations are imposed according to a sparse grids sampling of the parameter space. The sampling scheme chosen is similar to a design of experiments method but uses functions designed to minimize the number of samples required to achieve a desired degree of accuracy. The perturbed input parameters were connected to engine performance output over a large parameter space. Additionally, the same methodology was used to perform a set of diesel-only simulations at similar operating conditions to assess the numerical capability of capturing relatively low combustion variation. Lastly, this technique is expected to be useful to understand and predict combustion stability limits (knock and misfire) which are important for designing dual-fuel engines with improved operability and reliability.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Measurement of charged-current $\nu_\mu$ and $\bar{\nu}_\mu$ cross sections on hydrocarbon in a shallow inelastic scattering region

This MINERvA analysis is the first neutrino and antineutrino study of a shallow inelastic scattering region, which is the transition region between resonant production and deep inelastic scattering processes. This transition is explicitly included in this study by expanding the scope of shallow inelastic scattering to include not only the mainly lower-$Q^2$ nonresonant pion production but also the kinematic region where pion production off quarks within the nucleon becomes significant with $Q^2$ below the onset of the deep-inelastic scattering region defined in this analysis. To reduce the resonance background the kinematic region 1.5 $<$$W$$<$ 2 GeV/$c^2$ was chosen. In addition to the inclusive differential cross section measurements, to emphasize SIS interactions off quarks within the nucleon a sample with $Q^2$$\geq$ 1 GeV/$c^2$ was also analyzed. The measurements of one-dimensional cross sections at $\left\langle E_\nu \right\rangle \sim 6$ GeV on hydrocarbon of $Q^2$, Bjorken x and muon momentum variables are compared with modified predictions from the GENIE 2 neutrino generator as well as predictions of other neutrino simulators GiBUU, NEUT, NuWro and an alternative GENIE 3 version. Significant discrepancies both in shape and magnitude between measurements and neutrino simulator predictions of all variables have been observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Marine and continental stratocumulus cloud microphysical properties obtained from routine ARM Cimel sunphotometer observations

This study investigates marine and continental stratocumulus (Sc) cloud properties obtained from an automated implementation of a multispectral photometer retrieval. Photometer methods simultaneously retrieve cloud optical depth (τ) and cloud droplet effective radius (r e ), with estimates for liquid water path (LWP) calculated on the availability of those quantities. These applied methods evaluate retrieved cloud properties for Sc identified during a recent 6 year period over the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program sites in Oklahoma, USA (SGP) and in the Azores, Portugal (ENA). Modest agreement in key quantity retrievals is found between the routine photometer products and multisensor collocated profiling references. Cumulative breakdowns contingent on cloud thickness indicate increases in all retrieved quantities in thicker clouds, with larger discrepancies in the relative performance between the retrievals collected in the presence of drizzle. Under continental cloud conditions, the clouds of a similar thickness and r e to those sampled under marine conditions report a factor of 1.5 larger τ and LWP. An r 2 ≅0.65 is found between photometer τ retrievals and shadowband radiometer measurements, with photometer retrievals reporting a high (relative) bias. The τ intercomparisons indicate that variability between retrievals is a factor of three larger than errors reported from individual retrieval input perturbation tests. Photometer r e retrievals suggest a low r 2 (< 0.1) having a standard deviation ≅ 3 µm when compared to ARM baseline multi-sensor radar/radiometer references (accounting for offsets in the cloud droplet number concentration assumptions of the latter). However, photometer LWP calculations remain relatively unbiased in non-drizzling conditions, with errors O (50 g m −2 ) and r 2 ≅0.5 to collocated radiometer and interferometer references. Additional sensitivity tests for island influences on marine Sc properties suggest that while island-influenced winds may promote larger cloud LWP or thickness, the influence could be within retrieval method uncertainty and/or collocated instrument variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Target allocation yields for massively multiplexed spectroscopic surveys with fibers

We present Simulated Annealing fiber-to-target allocation simulations for the proposed DESI and 4MOST massively multiplexed spectroscopic surveys, and for both Poisson and realistically clustered mock target samples. We simulate both Echidna and theta-phi actuator designs, including the restrictions caused by the physical actuator characteristics during repositioning. For DESI, with theta-phi actuators, used in 5 passes over the sky for a mock ELG/LRG/QSO sample, with matched fiber and target densities, a total target allocation yield of 89.3% was achieved, but only 83.7% for the high-priority Ly-alpha QSOs. If Echidna actuators are used with the same pitch and number of passes, the yield increases by 5.7% and 16% respectively. Echidna also allows a factor-of-two increase in the number of close Ly-alpha QSO pairs that can be observed. Echidna spine tilt causes a variable loss of throughput, with average loss being the same as the loss at the rms tilt. With a natural tilt minimization scheme, we find an rms tilt always close to 0.58 x maximum. There is an additional but much smaller defocus loss, equivalent to an average defocus of 30microns. These tilt losses offset the gains in yield for Echidna, but because the survey strategy is driven by the higher priority targets, a clear survey speed advantage remains. For 4MOST, high and low latitude sample mock catalogs were supplied by the 4MOST team, and allocations were carried out with the proposed Echidna-based positioner geometry. At high latitudes, the resulting target completeness was 85.3% for LR targets and 78.9% for HR targets. At low latitude, the target completeness was 93.9% for LR targets and 71.2% for HR targets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

Simulating field conditions for the curing of salt-stone - 15501

salt-stone is a cementitious, low-level waste (LLW) material processed and stored at the Savannah River Site (SRS). It consists of fly ash (FA), blast furnace slag (BFS), and Portland cement that are blended with LLW salt solution. The dry feed materials and the salt solution are mixed at the salt-stone Production Facility (SPF), and the fresh grout subsequently pumped into a salt-stone Disposal Unit (SDU) in the salt-stone Disposal Facility (SDF) where it cures to form a rigid, cementitious monolith. In order to demonstrate the long-term performance of salt-stone, non-radioactive simulant samples are processed and cured in a laboratory environment and subsequently tested with respect to key performance properties. Of particular relevance to assessing the long-term stabilization and immobilization of radionuclides in salt-stone is the permeability, or saturated hydraulic conductivity, which determines the rate of infiltration for oxygen-bearing ground waters that can affect redox sensitive contaminants, such as {sup 99}Tc. Environmental conditions in the SDU, in particular temperature and humidity, are believed to affect the cured properties of salt-stone. With respect to humidity it has been determined that the average vapor space humidity in an SDU is at least 95%, and in addition the assumption can be made that underlying grout is subjected to a moisture-saturated environment by the grout poured over it. With respect to temperature, the SDUs contain thermo-well arrays that provide detailed vertical, and to a lesser degree radial, temperature profiles. Temperature data from the SDUs have subsequently been utilized to better simulate curing temperatures for simulant salt-stone in the laboratory. However, no two thermal profiles are the same, varying with respect to the following parameters: initial temperature, rate of temperature increase, maximum temperature attained, time at maximum temperature, and cooling rate. These temperature variables are dependent on the depth of grout (due to insulating effects of underlying and overlying grout), frequency of grout addition to the SDU (with more frequent additions causing grout temperatures to rise due to the exothermic hydration reactions), and, to a lesser degree, the ambient temperature of the external environment. Due to the variability in each thermal profile recorded in an SDU, the question arises as to which profile (or profiles) should be used in the simulated laboratory environment to best mimic the curing of field emplaced salt-stone, such that there is confidence that the properties measured for laboratory-prepared samples are representative of the grout in the SDU. Since the current SDU design contains sixty-six thermo-wells it is obviously not practical to process simulant samples according to each profile. Rather the aim would be to choose a reduced number of profiles that adequately encompass the range of recorded profiles, or to choose a single profile that is considered most conservative with respect to the resulting properties of interest. This is typically thought to be the profile exhibiting the highest peak temperature since for cement-based materials higher curing temperatures (> 60 deg. C), while resulting in enhanced reaction kinetics and strength development for short curing durations (∼28 days), have been shown to adversely affect the physical properties of the material after extended cure times. In choosing representative profiles it is important to consider not only the maximum temperature, but also the influence of rate of temperature rise to the maximum temperature, and time at or close to the maximum temperature on the post-cured properties. The intent of this paper is to consider the varied curing temperature profiles to which the salt-stone is exposed in an SDU, and to analyze the currently available data for elevated temperature curing of salt-stone to determine the potential thermal effects on the saturated hydraulic conductivity of salt-stone emplaced in an SDU. This may subsequently facilitate a more informed approach to simulating the curing of salt-stone in a laboratory environment. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Environmental fate of radiocesium in biota inhabiting a contaminated ecosystem on the U.S. Department of Energy's Savannah River Site

Although biomagnification of radiocesium ( 137 Cs) has been reported in food webs, most previous research has been limited to select trophic linkages. Few studies have included a comprehensive survey of fauna associated with aquatic, semi-aquatic, and terrestrial habitats within a single study framework. The objectives of this study were to advance our understanding of the dynamics of 137 Cs accumulation within food webs by quantifying 137 Cs activity across a wide range of biota found within a contaminated canal, as well as test the hypothesis that life-stage and body size influence 137 Cs bioaccumulation in select herpetofauna. With extensive sampling across multiple taxa collected from a contaminated canal system and associated floodplain on the Savannah River Site, we assessed 137 Cs activity and stable nitrogen isotopes for both aquatic organisms that were restricted to the contaminated effluent canal, and semi-aquatic organisms able to move freely between the contaminated canal and the adjacent uncontaminated terrestrial habitat. We found 137 Cs activity to be highly variable among species, with evidence for and against biomagnification in semi-aquatic and aquatic organisms, respectively. Furthermore, 137 Cs activity decreased with life stage and body size in bullfrogs (Lithobates catesbeianus), despite post-metamorphic bullfrogs having a more carnivorous diet compared to tadpoles, while cottonmouths (Agkistrodon piscivorus) retained similar 137 Cs activity regardless of their age and size. While evidence of biomagnification has been observed in some contaminated systems, results of our study suggest the extent to which 137 Cs biomagnifies within food webs is context-dependent and likely influenced by a suite of biotic and abiotic factors. Further, our data indicate sampling of a broad suite of species and environmental attributes are needed to elucidate the fate and dynamics of anthropogenic pollutants within contaminated ecosystems.

54 ENVIRONMENTAL SCIENCES↗