Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radionuclide transport”

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 307 records · Page 17

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

Infiltration Data Package for the E-Area Low-Level Waste Facility Performance Assessment

The infiltration data package contains the input parameters, the cap design and material properties assumptions, and the modeling results for the Hydrologic Evaluation of Landfill Performance (HELP) infiltration model simulations performed in support of the E-Area Low-Level Waste Facility Performance Assessment (PA). The infiltration estimates establish the upper boundary condition for the PORFLOW vadose-zone model and GoldSim model simulations for the following E-Area disposal unit types: Slit and Engineered Trenches, Low-Activity Waste Vault (LAWV), Intermediate-Level Vault (ILV), Componentin- Grout (CIG) special waste form trench segments, and the Naval Reactor Component Disposal Areas (NRCDAs). The infiltration data package builds upon relevant, foundational PA technical reports and memoranda from the past 15 years and is supported by three important components: the HELP model input parameter datasheets, HELP model input and output filenames and directory structure, and infiltration rates as a function of time for each scenario for each disposal unit type. For Slit and Engineered Trenches, a single best-estimate Bahia-grass case of two-percent slope and 585- foot slope length represents an upper bound on intact closure-cap infiltration rates over a 10,000-year period. Building upon the single, bounding, and intact case, two different approaches were considered for incorporating the effect of localized cap subsidence on infiltration rates in the PORFLOW flow and transport model simulations. In the first approach, a Monte Carlo probabilistic model was developed to calculate infiltration rates for 0.54%, 2%, 3.6%, and 4.9% subsidence scenariosa reflecting historical and future non-crushable content for the Slit and Engineered Trenches. The resulting blended (spatially and slope-length averaged) infiltration rates become an inflow upper boundary condition for the PORFLOW vadose zone model simulations. In the second and preferred approach, a Monte Carlo rectangle packing algorithm (Danielson, 2019) and a simplified equation for the total mass flux of water draining into a subsided compartment (Equation 1) will be employed to implement a weighted blending of radionuclide fluxes to the water table using the results of deterministic PORFLOW vadose zone simulations representative of specific subsidence cases. Finally, HELP-model results from sensitivity studies of the intact infiltration case were fit to a log-logistic function to generate infiltration profiles over a 10,000-year period for most-optimistic, more-optimistic, best-estimate, more-pessimistic, and most-pessimistic cases, for both the intact and four subsidence scenarios. For the LAWV and IL vaults, infiltration rate profiles for a 10,000-year period were developed for both an on-vault (above the concrete/metal vault roof) and an off-vault (10-foot soil zone adjacent to vault walls) scenario. The purpose of the off-vault simulations was to confirm that subsurface runoff from the concrete vault roof will adequately drain through the lowermost backfill layers adjacent to the vault walls. The LAWV and ILV infiltration profiles assume collapse of the concrete roof at relative Years 2,905 and 7,100, respectively. The identical Bahia-grass case with two-percent slope and 585-foot slope length assumed for the Slit and Engineered Trenches will represent an upper bound on intact infiltration rates for the NRCDAs and CIG special waste form as well. While the NRCDA casks are assumed to remain structurally intact (i.e., non-assumed to collapse 300 years after the end of operations. Because both rows of the existing CIG trench segments in ST23 (formerly CIG01) extend across the planned cap crest, the subsided CIG trench segments are treated conservatively in the HELP model as a 100%-subsidence case with zero subsurface run-on. The CIG special waste form trench segment infiltration profile, therefore, consists of an intact-cap period from relative Year 100 to Year 300, followed by a 100%-subsidence period from relative Year 300 to Year 10,100.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of flow on uranium speciation in soils impacted by acidic waste fluids

Radioactive acidic liquid waste is a common byproduct of uranium (U) and plutonium (Pu) enrichment and recycling processes whose accidental and planned release has led to a significant input of U into soils and sediments across the world, including at the U.S. DOE's Hanford site (WA, USA). Because of the particularly hazardous nature of U, it is important to predict its speciation when introduced into soils and sediments by acidic waste fluids. Of fundamental importance are the coupled effects of acid-driven mineral transformation and reactive transport on U speciation. To evaluate the effect of waste-fluid residence time and co-associated dissolved phosphate concentrations on U speciation in impacted soils and sediments, uncontaminated surface materials (from the Hanford Site) were reacted with U-containing synthetic acidic waste fluids (pH 2) amended with dissolved phosphate concentrations in both batch (no flow) and flow-through column systems for 7–365 days. By comparing dissolved U behavior and solid phase speciation as a function of flow regimen, we found that the availability of proton-promoted dissolution products (such as Si) to sequester U into uranyl silicates was dependent on waste fluid-sediment contact time as uranyl silicates were not detected in short contact time flow-through systems but were detected in no-flow, long contact time, reactors. Moreover, the dominance of uranyl phosphate as neoprecipitate U scavenger (principally in the form of meta-ankoleite) in phosphate amended systems confirmed the importance of phosphate amendments for an efficient sequestration of U in the soils and sediments. Overall, our experiments suggest that the formation of uranyl silicates in soils impacted by acidic waste fluids is likely to be limited unless reaction products are allowed to accumulate in soil pores, highlighting the importance of investigating soil U speciation in flow-through, transport-driven systems as opposed to no-flow, batch systems. Finally, this study provides insights into uranium speciation and its potential changes under acidic conditions for better prediction of risks and subsequent development of efficient remediation strategies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Can we safely go to 200 °C? An integrated approach to assessing impacts to the engineered barrier system in a high-temperature repository

Abstract. This presentation gives on overview of the complex thermo-hydro-mechanical-chemical (THMC) processes occurring during the disposal of heat-producing high-level radioactive waste in geologic repositories. A specific focus is on the role of compacted bentonite, which is commonly used as an engineered backfill material for emplacement tunnels because of its low permeability, high swelling pressure, and radionuclide retention capacity. Laboratory and field tests integrated with THMC modeling have provided an effective way to deepen our understanding of temperature-related perturbations in the engineered barrier system; however, most of this work has been conducted for maximum temperatures around 100 ∘C. In contrast, some international disposal programs have recently started investigations to understand whether local temperatures in the bentonite of up to 200 ∘C could be tolerated with no significant changes to safety relevant properties. Raising the maximum temperature is attractive for economical and safety reasons but faces the challenge of exposing the bentonite to significant temperature increases. Strong thermal gradients may induce complex moisture transport processes while geochemical processes, such as cementation and perhaps also illitization effects may occur, all of which could strongly affect the bentonite and near-field rock properties. Here, we present initial investigations of repository behavior exposed to strongly elevated temperatures. We will start discussing our current knowledge base for temperature effects in repositories exposed to a maximum temperature of 100 ∘C, based on data and related modeling analysis from a large heater experiment conducted for over 18 years in the Grimsel Test Site in Switzerland. We then show results from coupled THMC simulations of a nuclear waste repository in a clay formation exposed to a maximum temperature of 200 ∘C. We also explore preliminary data from a bench-scale laboratory mock-up experiment, which was designed to represent the strong THMC gradients occurring in a “hot” repository, and we finally touch on a full-scale field heater test to be conducted soon in the Grimsel Test Site underground research laboratory in Switzerland (referred to as HotBENT).

Birkholzer, Jens T. (ORCID:0000000279891912)↗

Codisposal Waste Package Loading Options for DOE SNF and HLW

The U.S. Department of Energy (DOE) is responsible for managing spent nuclear fuel (SNF) that is currently in, or will in the future come into, its possession. DOE must continue to safely store that SNF, transport it to an interim storage site or a repository, and dispose of it. These fuels come from a wide range of reactor types, such as light- and heavy-water-moderated reactors, graphite-moderated reactors, and breeder reactors, with various cladding materials and enrichments. Many of these reactors, now decommissioned, had unique design features, such as core configuration, fuel element and assembly geometry, moderator and coolant materials, operational characteristics, and neutron spatial and spectral properties resulting in a large diversity of reactor and fuel designs. Because of the wide variety and condition of SNF, a robust canister termed the DOE Standardized SNF Canister was proposed that would confine radionuclides and preclude moderator. The DOE Standard SNF Canister had four variations consisting of a 10- or 15-foot length and an 18- or 24-inch diameter. For ultimate disposal in the Yucca Mountain Repository these canisters were to be grouped with 24-inch diameter HLW canisters in a codisposal waste package. The 18-inch DOE Standardized SNF Canister could be placed in the middle of five HLW canisters. The 24-inch DOE Standardized SNF Canister would take the place of one the five HLW canister on the outer ring in the codisposal waste package. No DOE Standardized SNF Canisters have been loaded, and many fuel types can be loaded in multiple canisters. Depending on the selected diameter of the DOE Standardized SNF Canister for a particular fuel, the number of potential codisposal waste packages could vary. In addition, different types of DOE SNF may be able to be combined in a single canister. This could significantly change the number of DOE SNF canisters that are loaded and result in a decrease of the number of total potential codisposal waste packages. Ultimately, if too many HLW canisters need disposition, then a reduction in the number of DOE Standardized SNF Canisters will change the final total of codisposal waste packages very little. Likewise, if the ratio of DOE Standardized SNF Canisters to HLW canister produced is too high, then DOE Standardized SNF Canisters will not have enough HLW canisters to surround them in the codisposal waste package. This paper examines the results of differing loading strategies for a few DOE SNFs. It also compares the ranges of DOE Standardized SNF Canisters that may be produced to the ranges of HLW canisters that may be produced in order to determine potential scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Application of GOTHIC to Groundwater Transport Analysis - 20152

Migration of pollutants and hazardous wastes, potentially containing radioactive isotopes, via groundwater transport is a concern at most waste cleanup sites. Predictive analysis can be used to evaluate mitigating actions intended to minimize impact on the environment and public exposure. GOTHIC is a multipurpose thermal hydraulics code that is used extensively in the nuclear industry for design, licensing and operation evaluations. It combines the capabilities of typical one-dimensional system codes and the essential features of Computational Fluid Dynamics (CFD) codes for three-dimensional analysis. There are other codes that are specifically developed for groundwater transport analysis and the results presented here are consistent with prior analyses. However, GOTHIC has some unique features that offer advantages for applications related to nuclear waste. Most importantly, it has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B [2] and applicable portions of ASME NQA-1 [3] since 1995. Available GOTHIC capabilities that make the code especially useful for groundwater transport of nuclear materials include: - Tracking of any number of tracer elements for contaminants and other species of interest; - Radioactive decay and progeny of tracer elements; - Adsorption/desorption of tracer elements; - Tracking of any number of dissolved gases; - Release and absorption for dissolved gases; - Vapor phase tracking; - Non-Newtonian fluid modeling. The general porous body modeling approach makes GOTHIC well suited to groundwater transport analysis. The multi-region modeling approach used by GOTHIC simplifies model construction for regions of varying hydrologic characteristics and focuses the computational effort on regions of particular interest while simultaneously capturing the macroscopic response and any feedback effects across the larger domain. The applicability of GOTHIC to groundwater transport applications is demonstrated by comparing code results with available analytic or semi-analytic solutions for groundwater behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Electrokinetic In Situ Treatment of Metal-Contaminated Soil

An electrokinetic technique has been developed as a means of in situ remediation of soils, sludges, and sediments that are contaminated with heavy metals. Examples of common metal contaminants that can be removed by this technique include cadmium, chromium, zinc, lead, mercury, and radionuclides. Some organic contaminants can also be removed by this technique. In the electrokinetic technique, a low-intensity direct current is applied between electrodes that have been implanted in the ground on each side of a contaminated soil mass. The electric current causes electro-osmosis and migration of ions, thereby moving aqueous-phase subsurface contaminants from one electrode to the other. The half reaction at the anode yields H+, thereby generating an acid front that travels from the anode toward the cathode. As this acid front passes through a given location, the local increase in acidity increases the solubility of cations that were previously adsorbed on soil particles. Ions are transported towards one electrode or the other which one depending on their respective electric charges. Upon arrival at the electrodes, the ionic contaminants can be allowed to become deposited on the electrodes or can be extracted to a recovery system. Surfactants and other reagents can be introduced at the electrodes to enhance rates of removal of contaminants. Placements of electrodes and concentrations and rates of pumping of reagents can be adjusted to maximize efficiency. The basic concept of electrokinetic treatment of soil is not new. What is new here are some of the details of application and the utilization of this technique as an alternative to other techniques (e.g., flushing or bioremediation) that are not suitable for treating soils of low hydraulic conductivity. Another novel aspect is the use of this technique as a less expensive alternative to excavation: The cost advantage over excavation is especially large in settings in which contaminated soil lies near and/or under industrial buildings and therefore excavation would be made even more expensive by the need to prevent damage to numerous underground pipes and cables.

Quinn, Jacqueline↗

Evaluation of Observed and Modelled Aerosol Lifetimes Using Radioactive Tracers of Opportunity and an Ensemble of 19 Global Models

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (Cs-137) and xenon-133 (Xe-133) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. Cs-137 size distribution measurements taken close to the power plant suggested that accumulation mode (AM) sulfate aerosols were the main carriers of cesium. Hence, Cs-137 can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas Xe-133 behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of Cs-137that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and Xe-133 emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled Cs-137and Xe-133 concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime e, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95% confidence interval 13.1-15.7 days). The equivalent modelled e lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.42.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (Xe-133) concentrations in the Arctic as well but to a smaller extent than for the aerosol (Cs-137) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

Radioactive isotopes↗

A Rapid Sequential Separation Method for Determination of Actinides and TENORM in Fracking Wastes - 20223

With recent advances in unconventional drilling technology in the US and around the world, the risks of environmental contamination and exposure due to technologically enhanced naturally occurring radioactive materials (TENORM) from fracking wastes have also surfaced. Although the Permian Basin (southeastern New Mexico, western Texas) has long been a modest producer of oil, the advent of fracking technology a decade ago turned it into hot property for producers. A 2018 assessment of undiscovered, technically recoverable continuous oil and gas resources by the U.S. Geological Survey estimates an average 46.3 billion barrels of oil and 281 trillion cubic feet of gas. This is likely to increase the use of unconventional drilling and TENORM generation. TENORM concentrations in oil and gas exploration and production waste can be several hundred to several thousand times higher than background TENORM concentrations. Both the thorium and uranium decay chains contribute to airborne radionuclides arising from the radon gas escaping the ground and subsequently decaying as airborne particulates. Additionally, radon decay products {sup 210}Pb and {sup 210}Po can build up in scale on the internal surfaces of oil and gas handling pipes and in sludge in refineries, becoming potential inhalation and ingestion hazards for workers. Furthermore, southeastern NM is also home to world's only licensed and operating transuranic nuclear waste repository, the DoE's Waste Isolation Pilot Plant (WIPP). Plutonium isotopes ({sup 239+240}Pu) and {sup 241}Am, are expected to account for more than 99% of the total radioactivity scheduled for disposal in the WIPP repository. Thus, an improved understanding of the environmental fate and transport of TENORM, liberated by unconventional drilling, is essential to assess how best to protect individuals and the environment. In this context, accurate measurement of TENORM and actinides in environmental and biological samples is essential. In this presentation, a new sequential method for the separation and pre-concentration of actinides (Pu, Am, Np) and TENORM (Po, U and Th) derived from oil and gas exploration is proposed. The TEVA method involves a rapid co-precipitation step to remove matrix interferences from the samples, followed by plutonium oxidation state adjustment to Pu (IV) and an incubation period of ∼ 1 hour at 50-60 deg. C to allow the resultant Po (II) to oxidize into Po(IV). The polonium, neptunium, thorium and plutonium are then separated on a TEVA column, while americium and uranium are separated on a TRU column. After separation, the alpha counting source was prepared by micro-precipitation with copper sulfide (CuS) for polonium and neodymium fluoride (NdF{sub 3}) for actinides. The efficiency and reliability of the procedures were tested by analyzing filter, drinking water and frack sand samples. The method is simple, robust and can be performed quickly with excellent removal of interferences, high chemical recovery and very good alpha peak resolution. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

System Model Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

Two lined trenches (Trench 31 and 34) located within the 200 West Low-Level Burial Ground (LLBG) area of the U.S. Department of Energy’s (DOE) Hanford Site are designated for permanent disposal of low-level radioactive waste (LLW) and mixed low-level radioactive waste (MLLW). In accordance with DOE O 435.1, Radioactive Waste Management, radioactive waste shall be managed and disposed in a manner that is protective of worker, public health and safety, and the environment. DOE O 435.1 requires a site-specific radiological performance assessment (PA) that includes calculations of potential releases and subsequent doses to members of the public for a period of 1,000 years after closure of a low-level waste (LLW) disposal facility. The purpose of this environmental calculation file (ECF) is to document the various dose calculations performed to support DOE’s LLBG PA. The dose assessments are performed to evaluate the potential exposure of an all-pathway representative person to radionuclide contaminants of potential concern (COPCs) that may be released from the LLBG Trenches 31 and 34 to the point of calculation (POC) located at the outer edge of a 100 m buffer zone surrounding the trench boundary. This ECF uses inputs from other ECFs and their associated model package reports (MPRs) supporting the LLBG PA. The five major objectives of this ECF are to present the results of the following: (1) Deterministic dose assessments for both groundwater and atmospheric exposure scenarios; (2) Radon-222 (Rn-222) flux emanating at the surface from the disposed wastes; (3) Uncertainty analyses of the groundwater pathway annual dose due to uncertainty in input parameters; (4) Sensitivity analyses of the input parameters that could potentially impact the fate and transport of the contaminants and doses; (5) Doses arising from inadvertent intruder exposure scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A laboratory-scale process for producing dilithium beryllium tetrafluoride (FLiBe) with dissolved uranium tetrafluoride

Flibe Energy, Incorporated (FEI)'s conceptual Lithium Fluoride Thorium Reactor (LFTR) incorporates a chemical processing facility aimed at recovering uranium and other valuable volatile radionuclides while managing harmful radionuclides from the used fuel. The fuel utilized in this reactor is a combination of dilithium beryllium tetrafluoride (Li 2 BeF 4 or FLiBe) and uranium tetrafluoride (UF 4 ), (FLiBe/U). FEI's plan involves extracting the uranium and other valuable volatile fluoride-forming radionuclides using nitrogen trifluoride (NF 3 ). To facilitate laboratory-scale testing of uranium extraction using NF 3 and address the toxicity and physical hazards associated with beryllium and beryllium fluoride (BeF 2 ), we used a two-step process to prepare the simulated fuel salt. The first step entailed thermally decomposing ammonium beryllium tetrafluoride [(NH 4 ) 2 BeF 4 ] (ABeF) through a nominal 3-step process, combined with appropriate amounts of lithium fluoride (LiF) and UF 4 , resulting in the formation of beryllium fluoride (BeF 2 ). In the second step, the mixture was repeatedly melted and frozen at the melting point of FLiBe to prepare the eutectic FLiBe with dissolved UF 4 . Although the concept appears straightforward, the production of FLiBe/U involved various challenges. These challenges included transporting the gaseous decomposition products of ABeF, hydrogen fluoride (HF) and ammonia (NH 3 ), while preventing the formation of ammonium fluoride (NH 4 F). Additionally, it was necessary to control the reaction between the higher-than-anticipated water content in the commercial ABeF with NH 3 , HF, and the condensed NH 4 F, protect UF 4 from forming an unknown black compound, select suitable structural materials to mitigate fluoride corrosion, address the risks associated with beryllium toxicity through equipment design and operational protocols, and monitor process conditions. This article provides an account of the thermal decomposition chemistry observed in the commercial ABeF, describes the FLiBe/U production apparatus, describes the experiences and process refinements developed to prepare FLiBe/U, and presents our characterizations of prepared FLiBe/U.

Ammonium beryllium fluoride thermal decomposition↗

A Limited-In-Scope Comparison of Subsidence Scenarios For 3D Vadose Zone PORFLOW Trench Models

This technical memorandum summarizes the results from two separate sequences of 3D PORFLOW vadose zone flow and contaminant transport models that were used to test the impact of different spatial distributions of subsidence boundary conditions. More specifically, the two simulation setups used a general conceptual model for an engineered trench (ET): one with a subsided region specified in the center of the trench and one with a subsided region located at the end of the trench. Both subsided regions had the same geometric specifications (e.g., length, width, orientation), but used water infiltration rate boundary conditions that appropriately accounted for surface runoff and lateral drainage from the adjacent upslope intact regions of the final closure cap. Because of the substantially higher water infiltration rate associated with the subsided region at the end of the trench, this geometry produced higher absolute peak fluxes to the water table for all six radionuclides that were simulated. However, the slower release rate of the centrally located subsided hole geometry produced a higher flux to the water table in the time period subsequent to the absolute peak. A third model was implemented as a control, where the water infiltration rate that corresponds to the centrally located subsided region was applied at the end of the trench to highlight the purely geometric effects. The results show that a higher water infiltration rate produces a sharper peak, while a lower infiltration rate leads to a broader peak. A centrally located hole impacts more of the waste zone, producing a higher peak when an equal infiltration rate is applied. These generalized qualitative differences, and their impacts on the overall modeling workflow (i.e., from vadose zone models to limits and dose calculations), should receive appropriate consideration during the development of conceptual models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Calendar Year 2020 Post-Closure Monitoring Report for Corrective Action Unit 99: Rainier Mesa/Shoshone Mountain, Underground Test Area, Nevada National Security Site, Nevada (Rev. 1)

This report presents the results of monitoring conducted for water quality, water levels, and institutional controls in calendar year (CY) 2020 by the U.S. Department of Energy (DOE), Environmental Management (EM) Nevada Program’s Underground Test Area (UGTA) Activity at Corrective Action Unit (CAU) 99, Rainier Mesa/Shoshone Mountain (RM/SM), at the Nevada National Security Site (NNSS), Nevada. This report also presents analytical laboratory results (including results of quality assurance/quality control samples, such as field duplicates [FDs]), and verification of use restrictions (URs), institutional controls, and water usage. Groundwater samples were collected; water levels were measured; and well site surveillance was conducted in support of the Underground Test Area Closure Report (CR) for Corrective Action Unit 99: Rainier Mesa/Shoshone Mountain, Nevada National Security Site, Nevada (DOE/EMNV, 2020). Rainier Mesa (RM) was the site of 61 underground nuclear tests, and Shoshone Mountain (SM) was the site of 6 underground nuclear tests. As a result of these activities, which took place from 1957 to 1992, radionuclides (RNs) were released in the subsurface in the vicinity of the detonations. These 67 underground nuclear tests are associated with 66 specific corrective action site (CAS) numbers, as listed in the Federal Facility Agreement and Consent Order (FFACO) (1996, as amended). Two tests, HURON LANDING and DIAMOND ACE (conducted simultaneously), are included within one CAS. The CR (DOE/EMNV, 2020) establishes the regulatory boundaries and regulatory boundary objectives, monitoring program, UR boundaries and URs, and other institutional controls agreed to by EM Nevada Program and the Nevada Division of Environmental Protection (NDEP) for closure of CAU 99. Regulatory boundaries for the RM/SM CAU were established to protect receptors of groundwater from RN contamination within the three downgradient groundwater basins that receive recharge from RM. The SM regulatory boundary was established to verify that RN contamination does not reach the lower carbonate aquifer below SM. The UR boundaries were established based on contaminant simulation forecasts from the flow and transport simulations with consideration of site-specific topographical controls. The URs associated with the UR boundaries were identified to protect onsite workers and the public from inadvertent exposure to contaminated groundwater as forecasted to occur assuming current conditions. The UR is intended to restrict activities that might expose workers to contaminated groundwater within the UR area. In addition to the URs, other institutional controls are established to monitor and limit access to groundwater. These include federal ownership and management in perpetuity, controlled access of the NNSS and surrounding areas, and reporting water use on the NNSS and surrounding hydrographic basins.

54 ENVIRONMENTAL SCIENCES↗

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

61 RADIATION PROTECTION AND DOSIMETRY↗

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal electron sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comprehensive Analysis of Streaming and Shutdown Dose Rate Experiments at JET with ORNL Fusion Neutronics Workflows

Current experimental fusion systems and conceptual designs of fusion pilot plants (FPPs) are growing in complexity and size. Several radiation metrics are crucial to the safe operation of fusion machines, including neutron flux streaming through openings and the shutdown dose rate (SDDR). Most current designs of advanced experimental fusion systems—and the most probable candidates for FPPs—are based on the tokamak concept, which is prone to neutron streaming through the myriad openings needed for diagnostic and support systems. SDDR is caused by decay gamma rays from radionuclides that become activated by neutrons during the operation of a fusion system that use deuterium-deuterium (DD), tritium-tritium, or deuterium-tritium plasma. Because computational tools have become essential for determining these radiation metrics, they must be validated against reliable and applicable experimental data. Experiments at the Joint European Torus (JET) provide a unique source of experimental data for validating computational tools and nuclear data used to determine SDDR and neutron fluxes in streaming-dominated geometries. Here, this paper presents the comprehensive analysis of the high-performance DD JET SDDR, and streaming experiments performed using Oak Ridge National Laboratory (ORNL) fusion workflows. The computational results were compared with experimental results that consist of online SDDR measurements with ionization chambers and neutron fluence streaming measurements using thermoluminescent detectors. The ratio of calculated-to-experimental SDDR values ranges from 0.6 to 2.5, and the streaming results range from 0.5 to 8.0. Future work will include analyzing the JET 2021 DTE2 campaign alongside the integration of the Shift Monte Carlo transport code into all ORNL fusion neutronics workflows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the dissolution of a borosilicate glass with the use of isotopic tracing – Insights into the mechanism for the long-term dissolution rate

An understanding of the mechanisms responsible for controlling the long-term corrosion rate of nuclear waste glass is paramount if reliable glass dissolution models are to be used to calculate the controlled release of radionuclides from nuclear waste glass under geologic disposal conditions. Additionally, understanding silicate glass dissolution is also important for natural glasses to understand the role from the dissolution of these glasses has on the composition of natural aquatic systems. Two general mechanisms appear to be responsible for the elemental release from glass to the surrounding biosphere - ion exchange (release principally of alkalis) and matrix dissolution (release of the structural components of the glass). Key unknowns related to these mechanisms are the impact of surface layers on the altered glass and how these layers participate in corrosion. To better understand the impact of these layers on glass corrosion, we made two glasses of the same elemental composition, but with selected elements enriched in specific isotopes, e.g. more 10 B than 11 B. We reacted these glasses for 1 y then the solutions from each were exchanged and the chemistry of the solutions and the solids were followed for another 2 y and in some cases 3 y. This allowed us to follow changes in the glasses with respect to these elements and how these elements were transported through the alteration layers with time and, for at least one element, into the glass. Results from these experiments demonstrate that the release behavior of different elements is strongly dependent on their structural role in the glass (e.g. network formers and modifiers), but, more importantly, the role of water transport and subsequent ion exchange has in the long-term dissolution of the glass studied here. In this article, we highlight the behavior of four elements: lithium, sodium, silicon, and boron. Lithium and sodium, network modifying elements, have similar chemistries. The behavior of Li could be tracked in more detail by following both 6 Li and 7 Li; Li penetrates through the gel layer in both directions without hinderance and into and released from the glass deeper than sodium. Silicon, a network forming element, reacts with the silica-rich alteration layer. Boron, a network forming element, does not accumulate in the gel or the pristine glass and has a very sharp elemental profile between the pristine glass and alteration layer. Boron, lithium, and sodium elemental profiles suggest that there is little if any transport control within the alteration layer. The correlation between the Li and Na profiles and that of a water species suggest that the limiting release is a steady state between the diffusion of a water species and matrix dissolution that results in a low release rate for alkali and this provides a steady driving force for the long-term dissolution of glass.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Radiation Dose Modeling for Niowave’s Accelerator Driven Uranium Target Assembly 3

Molybdenum-99 is a high-value radionuclide commonly used for medical purposes within the United States. The National Nuclear Security Administration (NNSA) seeks to reliably produce the radioisotope 99 Mo without the use of highly enriched uranium. NNSA’s Office of Material Management and Minimization (M3) provides funding and government laboratory expertise to private companies to expedite the production process domestically and currently funds designs that use low-enriched uranium or other 99 Mo production pathways. Several production designs are being explored across the industry, including uranium fission and photonuclear conversion of 100 Mo targets. Niowave Inc. seeks to produce 99 Mo via a high-energy electron accelerator that strikes a lead-bismuth eutectic target that ultimately produces a consistent neutron flux. The neutron flux then interacts in a subcritical reactor core configuration to produce fission in low-enriched or natural uranium targets. These fissionable targets are then processed to extract 99 Mo. The purpose of this work is to estimate the neutron and photon dose response across Niowave’s proposed facility for worker safety during operation. Owing to the size of the proposed Niowave facility and necessary shielding, unbiased Monte Carlo radiation transport is impractical, and variance reduction methods are required. This work focuses on the weight window variance reduction method to produce high confidence dose response results within a Monte Carlo radiation transport code. Specifically, an adjoint-informed weight window methodology was created to improve the dose response estimates for accelerator-driven subcritical reactor designs. This adjoint-informed methodology was implemented for Niowave’s proposed design and improved dose results at far-field locations across the facility. Acceptable dose rate contours for the proposed facility were generated across the facility and are presented in this work.

07 ISOTOPE AND RADIATION SOURCES↗