Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Distribution”

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,315 records · Page 73

Combined, time-resolved, in situ neutron reflectometry and X-ray diffraction analysis of dynamic SEI formation during electrochemical N 2 reduction

One means of improving performance for electrochemical ammonia production through the Li-mediated N 2 reduction reaction (Li-NRR) is by cycling the current driving the reaction between open-circuit conditions and periods of applied current density. Herein, we have investigated the dynamics of the electrode–electrolyte interface under Li-NRR conditions during current cycling using in situ time-resolved neutron reflectometry and grazing-incidence synchrotron X-ray diffraction. During cycling, measured neutron reflectivity curves indicated bilayer formation in which Li-containing species such as LiOH, Li 2 O, and small quantities of Li 3 N and metallic Li primarily appeared in a thin layer at the cathode surface, above which formed a much larger, porous, ‘solid–electrolyte interface’ (SEI) layer. Upon return to open-circuit conditions, Li-containing species quickly moved out of the thin layer, leaving a compact, stable layer of decomposition products underneath the SEI layer. This SEI layer concomitantly filled with electrolyte or dissolved, becoming indistinguishable from the electrolyte via contrast in scattering-length density (SLD). During the second current cycle, Li-containing species again preferentially deposited directly atop the cathode, with the thick SEI-like layer again appearing within a minute. This SEI layer exhibited a lower SLD more quickly than in the first cycle, which might suggest that Li-containing species become distributed within the porous SEI layer. Furthermore, these time-resolved observations of SEI and plated layers during current cycling suggest that benefits associated with return to open-circuit conditions between periods of applied current density may be related to the concomitant loss of Li-containing species from a thin layer at the cathode surface into a porous SEI layer that becomes filled with electrolyte or dissolves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence for anti-synergism between ion-assisted etching and in-plasma photoassisted etching of silicon in a high-density chlorine plasma

Etching of p-Si in 60 mTorr 10%Cl2/90%Ar Faraday-shielded inductively coupled high density plasmas was investigated under both ionassisted etching (IAE) and photoassisted etching (PAE) conditions. Real-time etching rates and after-etching Si surface chemical compositions were obtained by laser interferometry and vacuum-transfer x-ray photoelectron spectroscopy (XPS), respectively. Precisely controlled ion energy distributions (IEDs) were generated by applying pulsed negative DC bias on the conductive sample stage. Above a 36 eV threshold at a total flow rate of 250 SCCM, the IAE rate increased with the square root of the ion energy. In contrast to the DC bias, etching under RF bias did not exhibit a threshold ion energy because of the wide IED. XPS spectra revealed that the surface layer under PAE conditions had a significantly lower chlorine content, composed of only SiCl. Under IAE conditions, however, silicon dangling bonds (Si•), SiCl 2 , and SiCl 3 were found on the surface, in addition to SiCl, with a relative abundance of SiCl > SiCl 2 , and SiCl 3 . The absence of higher chlorides and Si• under PAE conditions suggested that vacuum ultraviolet photons and above threshold-energy ions interact with the surface very differently. By varying the duty cycle of the pulsed DC bias, it was found that the IAE rate scaled with the energetic ion dose, but only for low duty cycles. For higher duty cycles, the apparent IAE yield fell off with an increasing Cl coverage on the surface, as the duty cycle went up, which pointed to a negative synergy (antisynergism) between PAE and IAE as the explanation. This antisynergism was further supported by the observed decrease of the total etching rate with an increasing period of the pulsed DC bias. A plausible mechanism is that increasing the pulsing period causes more near-surface damage, creating more recombination centers that lead to a higher loss rate of electron-hole pairs through recombination, thereby reducing the PAE rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Post synthetic amine functionalization of MOF-808 for CO 2 sorption via ligand exchange and Michael addition

Post synthetic modification of metal organic frameworks presents a viable route for amine functionalization which can significantly enhance CO 2 sorption capacity. We present a facile means of amine incorporation using limited synthetic steps and low-cost reagents that results in a high density of primary alkyl amines distributed through a zirconium-based metal organic framework (MOF). Both the MOF synthesis and the post synthetic modification take place under aqueous conditions and result in strongly bound molecular amines available for sorbate interaction throughout the MOF pores. Furthermore, this amine incorporation protocol results in a significantly increased CO 2 capacity compared with the unmodified MOF-808. Specifically, CO 2 isotherms collected at 298 K for the unmodified MOF-808 show uptakes of 0.06, 0.2, and 1.2 mmol/g at 4, 15, and 100 kPa, respectively, which can be compared with 0.3, 0.7, and 2.5 mmol/g for the glycine grafted MOF-808 and 0.5, 0.9, and 2.3 mmol/g for ethylenediamine grafted MOF-808.

Amine sorbent↗

A neutron tomography study to visualize fluoride salt (FLiNaK) intrusion in nuclear-grade graphite

Manufactured graphite is a preferred material for in-core components of molten salt reactors and fluoride salt-cooled high-temperature reactors, which are in permanent contact with liquid salts. However, owing to the porous nature of nuclear graphite, under certain conditions, molten salts may intrude graphite's pores and affect graphite's properties and functionality. Therefore, a better understanding of molten salt intrusion (distribution across sample cross section and penetration depth) is needed to assess its effects. Here, in this work, we have demonstrated the use of neutron imaging (computed tomography) in the evaluation of salt penetration and distribution of a wide range of graphite grades with diverse microstructures that have been subjected to FLiNaK (LiF–NaF–KF) intrusion at 750 °C and 5 bar pressure for 12 h. Because of the great neutron attenuation contrast from scattering and adsorption between Li (from FLiNaK) and the graphite matrix, we have obtained direct visualization of FLiNaK salt distribution in the salt-impregnated graphites for the first time. Three-dimensional reconstructed images and cross-sectional concentration profiles demonstrate that salt penetration and density distribution are greatly dependent on the microstructural properties of the graphite grade.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗

Resource assessment of ocean thermal energy conversion in Puerto Rico and U.S. Virgin Islands

Island communities often struggle to establish and maintain traditional electric grids and are therefore heavily reliant on costly imported fossil fuels. In the case of Puerto Rico, these challenges are enhanced by extreme weather and other natural hazards that threaten the local electricity generation and transmission infrastructure. Ocean thermal energy conversion (OTEC) could play an important role in establishing a more resilient electrical grid in the region. Here, in this study, a detailed analysis is conducted to characterize the ocean thermal resource and power potential of OTEC in Puerto Rico based on a 14-year dataset of modeled ocean temperature. The assessment considers seasonal and interannual variability in the region's thermal resource and examines the operational limitations associated with minimal thermal gradients required to run a typical OTEC heat engine. Notably, the local thermal resource is found to be sensitive to El Niño-Southern Oscillation (ENSO) climate patterns, with La Niña conditions linked to greater OTEC power availability. Seven areas of opportunity are identified based on their resource potential and proximity to existing electrical distribution lines, including two that could benefit the nearby U.S. Virgin Islands. The greatest OTEC power potential is observed to the south of the main island of Puerto Rico in the Caribbean Sea with an estimated capacity of 138 MW for a plant pumping cold water from a depth of 1,000 m, or the equivalent amount of electricity required to power 219,000 households.

OTEC↗

Characterization of the In Vivo Deuteration of Native Phospholipids by Mass Spectrometry Yields Guidelines for Their Regiospecific Customization

Customization of deuterated biomolecules is vital for many advanced biological experiments including neutron scattering. However, because it is challenging to control the proportion and regiospecificity of deuterium incorporation in live systems, often only two or three synthetic lipids are mixed together to form simplistic model membranes. This limits the applicability and biological accuracy of the results generated with these synthetic membranes. Despite some limited prior examination of deuterating Escherichia coli lipids in vivo, this approach has not been widely implemented. In this report an extensive mass spectrometry-based profiling of E. coli phospholipid deuteration states with several different growth media was performed, and a computational method to describe deuterium distributions with a one-number summary is introduced. The deuteration states of 36 lipid species were quantitatively profiled in 15 different growth conditions, and tandem mass spectrometry was used to reveal deuterium localization. Regressions were employed to enable the prediction of lipid deuteration for untested conditions. Small-angle neutron scattering was performed on select deuterated lipid samples, which validated the deuteration states calculated from the mass spectral data. Based on these experiments, guidelines for the design of specifically deuterated phospholipids are described. This unlocks even greater capabilities from neutron-based techniques, enabling experiments that were formerly impossible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray Diffraction of Water in Polyvinylpyrrolidone

PVP is a hydrophilic polymer commonly used as an excipient in pharmaceutical formulations. Here we have performed time-resolved high-energy X-ray scattering experiments on pellets of PVP at different humidity conditions for 1-2 days. A two-phase exponential decay in water sorption is found with a peak in the differential pair distribution function at 2.85 Å, which is attributed to the average (hydrogen bonded) carbonyl oxygen-water oxygen distance. Additional scattering measurements on powders with fixed compositions ranging from 2 to 12.3 wt % H 2 O were modeled with Empirical Potential Structure Refinement (EPSR). Further, the models reveal approximately linear relations between the carbonyl oxygen-water oxygen coordination number (nO C - O W ) and the water oxygen-water oxygen coordination number (nO W - O W ) versus water content in PVP. A stronger preference for water-water hydrogen bonding over carbonyl-water bonding is found. At all the concentrations studied the majority of water molecules were found to be randomly isolated, but a wide distribution of coordination environments of water molecules is found within the PVP polymer strands at the highest concentrations. Overall, the EPSR models indicate a continuous evolution in structure versus water content with nO W - O W =1 occurring at similar to 12 wt % H 2 O, i.e., the composition where, on average, each watermolecule is surrounded by one other water molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of Mixed-Wettability on Brine Drainage and Supercritical CO 2 Storage Efficiency in a 2.5-D Heterogeneous Micromodel

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO 2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO 2 (scCO 2 ) is used in enhanced oil recovery. Here, we performed scCO 2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from –8.1 to –4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO 2 saturation increase at both capillary fingering and crossover flow regimes (–8.1 ≤ logCa ≤ –6.1), the one heterogeneous wetting to scCO 2 results in ~0.09 saturation increase only at the crossover flow regime (–7.1 ≤ logCa ≤ –6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO 2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > –6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO 2 saturation and distribution suggest the importance of wettability on CO 2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

58 GEOSCIENCES↗

Optimal deformation procedure for hybrid adaptive x-ray mirror based on mechanical and piezo-driven bending system

Here, a hybrid deformable x-ray mirror consisting of a mechanical bender and a bimorph deformable mirror has been developed to realize adaptive optical systems, such as zoom condenser optics, for synchrotron-radiation-based x-ray microscopy. In the developed system, both bending mechanisms comprehensively contribute to the formation of the target mirror shape and can narrow the role of piezoelectric actuators, thereby enabling a more stable operation. In this study, the behavior of the bimorph mirror under the clamped condition was investigated, and the sharing of the deformation amount for each bending mechanism was optimized to minimize the amplitude of the voltage distribution of the bimorph mirror.

47 OTHER INSTRUMENTATION↗

Magneto-optical study of Nb thin films for superconducting qubits

Abstract Among the recognized sources of decoherence in superconducting qubits, the spatial inhomogeneity of the superconducting state and the possible presence of magnetic-flux vortices remain comparatively underexplored. Niobium is commonly used as a structural material in transmon qubits that host Josephson junctions, and excess dissipation anywhere in the transmon can become a bottleneck that limits overall quantum performance. The metal/substrate interfacial layer may simultaneously host pair-breaking loss channels (e.g. two-level systems) and control thermal transport, thereby affecting dissipation and temperature stability. Here, we use quantitative magneto-optical imaging of the magnetic-flux distribution to characterize the homogeneity of the superconducting state and the critical current density, j c , in niobium films fabricated under different sputtering conditions. The imaging reveals distinct flux-penetration regimes, ranging from a nearly ideal Bean critical state to strongly nonuniform thermo-magnetic dendritic avalanches. By fitting the measured magnetic-induction profiles, we extract j c and try to correlate it with film physical properties and with measured qubit internal quality factors. Our results indicate that the Nb/Si interlayer can be a significant contributor to decoherence and should be considered an important factor that must be optimized.

Datta, Amlan [Ames National Laboratory; Iowa State↗

Impact of gas injection location and divertor surface material on ITER fusion power operation phase divertor performance assessed with SOLPS-ITER *

Abstract The ITER divertor design and performance assessment, primarily based on the SOLPS-4.3 burning plasma database (Pitts R. et al 2019 Nucl. Mater. Energy 20 100696), assumes the use of beryllium (Be) as the divertor surface material and the injection of gas from the main chamber top. However, the current ITER baseline favors gas injection from the more toroidally symmetric sub-divertor region. This paper evaluates the implications of these assumptions for divertor performance in the ITER fusion power operation phase. The impact of the divertor surface material and the gas injection location on the main ions mirrors the hydrogen only low power phase scenario shown in Park J.-S. et al (2020 Nucl. Fusion 61 016021). However, during burning plasma operation, extrinsic impurity seeding will be required. In the case of neon (Ne), studied here, impurity retention is influenced by both the divertor surface material and the fueling location. Neon leakage increases due to more energetic reflection from tungsten than beryllium, but equivalent divertor performance can be achieved by adjusting the neon seeding rate. While the impurity seeding location does not affect the distributions of impurity or radiation, the fueling location does. Top fueling provides local ionization sources mainly in the mid-SOL under detached conditions, enhancing divergences of the flux there (source-driven flow), bringing stagnation points close to the fueling location, and equilibrating flows towards both targets. In contrast, the global flow pattern (in the absence of fluid drifts) in the case of sub-divertor fueling is biased towards the inner target. Impurity flows, driven by force balance, largely mirror those of the main ion flow, including the stagnation point. The case with top fueling enhances Ne retention and corresponding radiation in the outer divertor, effectively reducing the total and peak target heat fluxes by 20%–40%, compared to the case with divertor fueling. Meanwhile, the case with outer target fueling also achieves similar reductions by enhancing plasma-neutral interactions. These results suggest the possibility that the selection of the fueling location and throughput can be used as an actuator to control impurity divertor retention and divertor radiation asymmetry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Counter Data Paucity through Adversarial Invariance Encoding: A Case Study on Modeling Battery Thermal Runaway

Lithium-ion batteries, widely used for their durability and high energy storage, face the risk of internal short circuits leading to catastrophic thermal runaway events. These events, triggered by external stimuli like mechanical loads, pose safety concerns in applications such as electric vehicles. Detecting and understanding thermal runaway events is crucial, but physics-driven models struggle to explain the non-linear evolution of battery temperature during these events, considering factors like material composition and state-of-charge. Due to the rarity of these events and the cost of data collection, we propose a deep learning (DL) model to predict battery temperature responses during thermal runaway. The challenge lies in the scarcity of data, making traditional DL models prone to overfitting and learning low-quality representations of the complex process.Our approach introduces a novel few-shot architecture that incorporates an adversarially governed invariant encoding process. This architecture aims to distill "invariant" relationships by addressing distributional shifts in data across various battery properties, facilitating the detection of thermal runaway events. Specifically, our results demonstrate that deep learning models conditioned on these "invariant" representations outperform state-of-the-art baselines, achieving a remarkable 96.8% performance improvement in terms of the popular metric MAPE. This framework presents a promising direction for enhancing battery safety modeling, particularly in the context of rare and complex events like thermal runaway. Our code and code and dataset used for the paper are public1.

Tabassum, Anika [ORNL] (ORCID:0000000254600955)↗

Diverse Super-Resolution (diversity_SR) [SWR-21-60]

Deep learning tools for enhancing the spatial resolution of wind data. The software is developed in Python using the TensorFlow deep learning package. Models for diversity super-resolution is provided. Included in the package are pretrained models with example code/data to perform the super-resolution as well as tools of training models for different enhancement- or data-types. The super-resolution is an inherently ill-conditioned problem, with multiple high-resolution fields plausibly mapping to the same coarse field. Considitional GANs provide a framework for generating a distribution of high-resolution realizations from a given low-resolution input. Stochastic estimation is used to inform the network of the expected degree and location of sub-grid diversity. The package includes a pretrained network to generate distributions of 10x-enhanced fields of wind data.

Glaws, Andrew↗

Dynamic pulse propagation modelling for predictive femtosecond-laser-microbonding of transparent materials

A dynamic pulse propagation modeling for femtosecond laser bonding of Borofloat glass is presented. The temperature evolution and internal modifications are predicted by incorporating the nonlinear electron dynamics along with temperature dependent thermal properties. The modelling predicts the spatial and temporal distribution of absorption coefficient and plasma density that gives quantitative estimations of the heat affected zone and weld geometry. The impact of focusing condition on heat affected zone and weld geometry is investigated, which for the first time to our knowledge allows to numerically determine the desired relative position between the geometrical focus of a femtosecond-laser-pulse and the interface of the two substrates to be welded. The prediction of the modelling on the offset distance is applied to weld Borofloat glass plates having optical contact and can be applied to other dielectric solids.

Sahoo, Pankaj K. (ORCID:0000000209437319)↗

Multi-Scale Lithofacies and Chemostratigraphic Analysis of Two Middle Devonian Marcellus Shale Wells in Northern West Virginia, USA

The Marcellus Shale Energy and Environmental Laboratory (MSEEL) consists of two project areas within the dry gas producing region of the Marcellus shale play in Monongalia County, West Virginia. MSEEL is a collaborative field project led by West Virginia University, with Northeast Natural Energy LLC, several industrial partners, and sponsored by the US Department of Energy National Energy Technology Laboratory. The study areas are drilled approximately 8.5 miles apart to better understand the vertical and lateral changes in stratigraphy over a short distance. Two vertical pilot wells, MIP-3H and Boggess 17H were drilled in the fall of 2015 and spring of 2019, respectively. Core was recovered from the MIP-3H (API: 47-061-01707-00-00) 112 feet (34m) between depths of 7445 to 7557 feet, and from the Boggess 17H (API: 47-061-01812-00-00) 139 feet (42m) between depths of 7908 and 8012 ft. A full suite of triple combo (gamma ray, neutron, density logs), image logs, and advanced logging tools were run in both wells and calibrated to core analysis. Core analysis includes medical computed tomography (CT) scans, mineralogy and chemostratigraphy determined from handheld X-Ray fluorescence (hhXRF) and X-Ray powder diffraction (XRD) measurements, and determination of total organic content (TOC). Lithofacies were determined at core-scale using traditional core description techniques and medical CT-scan images. Log-scale facies are based on mineralogy and TOC data and developed using petrophysical logging data calibrated to core data (XRD and pyrolysis data). Chemostratigraphic analysis utilized hhXRF data to determine the major and trace element trends in the cores. In the two wells six shale lithofacies were recognized at the core and log scale. Both wells show organic-rich facies (TOC > 6.5%) primarily in the middle and lower Marcellus, with a slight decrease in thickness of this interval in the Boggess 17H. This interval is interpreted as an increase in paleo-productivity (increased Ni, Zn, and V), decreased sedimentation (decreased detrital proxies), and anoxic to euxinic conditions (increased Mo and chalcophile elements). Paleo-redox conditions in both wells are dynamic throughout deposition transitioning between euxinic/anoxic to dysoxic/oxic. This trend is seen through elemental proxies and calcite/pyrite concretion distributions.

03 NATURAL GAS↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K2O, LaBr3, KBr, Cu, Zr, F, Ga, Ga2O3, NaNO2, Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated assuming 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradioactive chemical

54 ENVIRONMENTAL SCIENCES↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradradioactive chemicals modeled. And finally, the calculations assume all 12 tests will be performed at the same place at both locations, and all 11 radioactive material types will be used for each test. This is conservative because it is anticipated that no more than two material types will be used per test.

99 GENERAL AND MISCELLANEOUS↗