Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dose Rate Calculations”

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 91 records · Page 5

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Air Exchange Rate Impact on Activity Equilibrium Factors and Inhalation Fractional Equilibrium Factors for Rn, Xe, Kr, Ar, Ne, and Their Progeny in Vapor Intrusion, Risk, and Dose Models

Exposure to the radioactive noble gasses, especially radon, is of high concern and poses a significant risk to humans in an indoor air environment as the second leading cause of lung cancer in the United States. To evaluate and minimize the risks posed by these gasses, it is important to understand their radiological and physical properties. The EPA’s Radon Vapor Intrusion Screening Level (RVISL) calculator calculates indoor air RVISLs based on target working levels (WLs), target excess lifetime cancer risk (ELCR), and annual dose limits for the actinon (Rn-219), thoron (Rn-220), and radon (Rn-222) decay series. The RVISLs are based on inhalation and submersion in gas cloud exposure routes for residential and commercial settings. The RVISLs are analogous to preliminary remediation goals (PRGs) and dose compliance concentrations (DCCs), where the isotope-specific values are in units of activity concentration (activity per unit volume). If the concentration of a parent isotope of radon or its progeny is found to exceed the RVISL, then further action to ensure cleanup of the contaminant may be necessary. In residential and commercial settings, the RVISLs will vary based on the air exchange rate present. The EPA’s Radionuclide PRG and DCC Calculators also assess the risk/dose from noble gases in the air due to household use of water like showering. In this study, a computational method in MATLAB was developed to determine the impact of the air exchange rate on the activity equilibrium factor (A eq ) and the inhalation fractional equilibrium factor (F eq ). Both factors are values that reflect the equilibrium concentrations of progeny to their parent in the air. These factors have a direct impact on the RVISL, PRG, and DCC calculations of WL, ELCR, and annual dose, respectively. This study builds on a previous report that only focused on actinon, thoron, and radon by revisiting the original A eq and F eq calculation methods, as well as including the values for the Rn-207, Rn-209, Rn-210, Rn-211, Rn-215, Rn-216, Rn-217, Rn-218, Rn-223, Ne-24, Ar-42, Ar-43, Ar-44, Kr-74, Kr-75, Kr-76, Kr-77, Kr-88, Kr-89, Xe-120, Xe-121, Xe-122, Xe-123, Xe-135m, and Xe-138 decay chains, which are not currently available in literature. The EPA’s RVISL calculator will be updated to include the new A eq and F eq values for the actinon, thoron, and radon decay chains, while the rest of the calculators will incorporate all the new A eq and F eq values as appropriate.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Total Radioisotope Yield Calculator

An initial number of fissions on a fissile/fissionable material, defined by ENDF library availability, is used to produce an initial inventory of fission products (FPs). The FPs are decayed using full, analytical solutions to the Bateman equations to calculate the FP inventory at a user-specified time. Output can contain information about activity, mass, dose rate, and gamma-ray intensities based on user-defined parameters, given as .txt or .xlsx file formats. Full execution and display of results takes 5-10 seconds, depending on the complexity of the input parameters.

Holschuh, ThomasV [Idaho National Laboratory (INL)↗

Dose Summation Strategies for External Beam Radiation Therapy and Brachytherapy in Gynecologic Malignancy: A Review from the NRG Oncology and NCTN Medical Physics Subcommittees

Definitive, nonsurgical management of gynecologic malignancies involves external beam radiation therapy (EBRT) and/or brachytherapy (BT). Summation of the cumulative dose is critical to assess the total biologic effective dose to targets and organs at risk. Cumulative dose calculation from EBRT and BT can be performed with or without image registration (IR) and biologic dose summation. Among these dose summation strategies, linear addition of dose-volume histogram (DVH) parameters without IR is the global standard for composite dose reporting. This approach stems from an era without image guidance and simple external beam and brachytherapy treatment approaches. With technological advances, EBRT and high-dose-rate BT have evolved to allow for volume-based treatment planning and delivery. Modern conformal therapeutic radiation involves volumetric or intensity modulated EBRT, capable of simultaneously treating multiple targets at different specified dose levels. Therefore, given the complexity of modern radiation treatment, the linear addition of DVH parameters from EBRT and high-dose-rate BT is challenging to represent the combined dose distribution. Deformable image registration (DIR) between EBRT and image guided brachytherapy (IGBT) data sets may provide a more nuanced calculation of multimodal dose accumulation. However, DIR is still nascent in this regard, and needs further development for accuracy and efficiency for clinical use. Biologic dose summation can combine physical dose maps from EBRT and each IGBT fraction, thereby generating a composite DVH from the biologic effective dose. However, accurate radiobiologic parameters are tissue-dependent and not well characterized. A combination of voxel-based DIR and biologic weighted dose maps may be the best approximation of dose accumulation but remains invalidated. The purpose of this report is to review dose summation strategies for EBRT and BT, including conventional equivalent dose in 2-Gy fractions dose summation without image registration, physical dose summation using 3-dimensional rigid IR and DIR, and biologic dose summation. We also provide general clinical workflows for IGBT with a focus on cervical cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Review of Experimental Data for Validating Computer Codes Used in Shielding Calculations for Spent Fuel Storage and Transportation Systems

This report presents a review of available radiochemical assay data and shielding benchmarks applicable to spent nuclear fuel (SNF) shielding calculations. The relevant information reviewed herein includes the Spent Fuel Composition (SFCOMPO) database, the Shielding Integral Benchmark Archive and Database (SINBAD), the International Handbook of Evaluated Criticality Safety Benchmark Experiments, and published measurements of external dose rates of casks loaded with SNF. The relevant experimental data identified in this report may be used to support verification and validation of computer codes used in SNF cask/transport shielding applications, as well as development of calculation uncertainties. It should be noted that a relatively small subset of the identified experimental data (e.g., criticality alarm experiments) is available in a standard format established by the international community participating in experimental isotopic and shielding data evaluations. An effort of the SFCOMPO Technical Review Group (TRG) is underway to publish first isotopic evaluations of individual assay data using a standard data evaluation format. The SINBAD TRG has recently initiated benchmark evaluations and modernization of the database. Therefore, more relevant information is expected in the future that will enable users to select quality experimental data in depletion code and shielding code validations for SNF applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Air Pathway Dose Modeling for the E-Area Low-Level Waste Facility

The US Department of Energy (DOE) Order 435.1 performance assessment (PA) process (USDOE 1999) prescribes a performance objective (10 mrem y -1 ) for evaluating atmospheric releases of radionuclides from DOE low-level waste (LLW) disposal facilities. The potential dose to an individual from exposure to radionuclides released into the atmosphere from LLW disposals can be estimated by application of radionuclide -specific dose-release factors (DRFs) (mrem Ci -1 ) to estimated flux rates (Ci y -1 ) at a particular time and location. In Revision 0 of this report, Lee (2006) calculated DRFs for potential atmospheric releases of 15 volatile radionuclides from seven ELLWF disposal units in support of the 2008 ELLWF PA (WSRC 2008). The mainframe version of the US Environmental Protection Agency’s (EPA) dose model CAP88 was used in the Revision 0 assessment. In Revision 1 of this report (Dixon and Minter 2017), the site-specific Savannah River National Laboratory (SRNL) atmospheric dose models MAXDOSE-SR version 2013 (Stone and Jannik 2013a) and MAXINE version 2017 (Bell 2017) were used to calculate DRFs for ten radionuclides. This work (Revision 2) supersedes Revision 0 (Lee 2006) and Revision 1 (Dixon and Minter 2017) of this report. The Revision 2 method for estimating new DRFs for the ELLWF disposal units is described in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Air Pathway Dose Modeling for the E-Area Low-Level Waste Facility

The US Department of Energy (DOE) Order 435.1 performance assessment (PA) process (USDOE 1999) prescribes a performance objective (10 mrem y -1 ) for evaluating atmospheric releases of radionuclides from DOE low-level waste (LLW) disposal facilities. The potential dose to an individual from exposure to radionuclides released into the atmosphere from LLW disposals can be estimated by application of radionuclide-specific dose-release factors (DRFs) (mrem Ci -1 ) to estimated flux rates (Ci y -1 ) at a particular time and location. In Revision 0 of this report, Lee (2006) calculated DRFs for potential atmospheric releases of 15 volatile radionuclides from seven ELLWF disposal units in support of the 2008 ELLWF PA (WSRC 2008). The mainframe version of the US Environmental Protection Agency’s (EPA) dose model CAP88 was used in the Revision 0 assessment. In Revision 1 of this report (Dixon and Minter 2017), the site-specific Savannah River National Laboratory (SRNL) atmospheric dose models MAXDOSE-SR version 2013 (Stone and Jannik 2013a) and MAXINE version 2017 (Bell 2017) were used to calculate DRFs for ten radionuclides. Revision 2 provided updated DRFs using the current version of the US Environmental Protection Agency’s (EPA) dose model CAP88-PC Version 4.1. After Revision 2 was finalized, DRFs were requested for Ar-37, Ar-39, Kr-83m, and Hg-206. Revision 3 provides the additional DRFs calculated using the same methods as used in Revision 2 of this document. The method for estimating new DRFs for the ELLWF disposal units is described in this report.

61 RADIATION PROTECTION AND DOSIMETRY↗

Evaluation of Oxygen Consumption as a Sensitive Measure of Electrical Cable Polymer Insulation Degradation

Polymers used as dielectric materials in nuclear grade electrical cables, such as cross-linked polyethylene (XLPE) and ethylene-propylene rubber (EPR), are long-lived and durable. Estimating the decades-long useful service life of polymers experimentally in a feasible timeframe of weeks or months requires greatly accelerated aging in the laboratory at elevated temperatures and/or gamma dose rates. Many limitations of accelerated aging, such as due to artefacts of the process including diffusion limited oxidation and dose rate effects, have been identified that can make accurate lifetime prediction at service conditions challenging. This can be due to differences in degradation mechanisms between rapid and extended aging or to the presence of induction periods during which early aging is undetectable using common metrics of polymer insulation characterization such as tensile elongation at break. Conceptually, an ideal measure of aging would track clearly and predictably with material changes throughout the lifecycle of the evolving material state from early aging through end of useful life. With such a metric accelerated aging could be performed at much milder conditions, conditions closer to those found in long term service, and end of life calculated from the early data. Oxygen consumption has previously been proposed as a sensitive measure of insulation polymer oxidation that might meet such desirable characteristics. In this work, nuclear cable insulation samples of the two most common types, XLPE and EPR, from two of the most sourced manufacturers, RSCC and The Okonite Company, were subjected to gamma irradiation at dose rates of 100, 200, and 300 Gy/hr for up to 42 days at 26°C. The aged samples were used to evaluate oxygen consumption as a sensitive measure of aging for these materials and these conditions. For comparison and validation, the aged samples were also characterized using the more conventional tests of mass change, density, carbonyl index, and yellowness index. For both materials, oxygen consumption during irradiation was found to track consistently with applied dose—oxygen consumption rate was linear with time. Mass change and carbonyl index were also clearly proportional to exposure dose. Yellowness index, a surface sensitive measure, was observed to increase linearly with dose at low exposure levels before plateauing. Sample density was not observed to change under the conditions explored. Oxygen consumption was found to be a particularly sensitive metric to measure the degradation of polymeric cable insulation exposed to radiation. Oxygen consumption was compared to Arrhenius and time-temperature-superposition methods for activation energy calculation using literature date and found to be a useful and viable alternative for lifetime prediction. The techniques developed in this work are anticipated to be valuable in the estimation of insulation aging rates under mild conditions approaching service conditions. Direction comparison of useful life prediction by oxygen consumption with prediction via Arrhenius methodology for thermal aging and existing methods for radiation and combined thermal/radiation aging would serve to further validate this promising approach and refine our understanding of material-specific aging rates.

36 MATERIALS SCIENCE↗

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

Here, the target segments of the Oak Ridge National Laboratory Second Target Station (STS) neutron production facility become highly activated due to spallation reactions or nuclei transmutation by primary protons and emitted neutrons. Once the target segments are removed from their location within the core vessel, decay dose rates must be accurately quantified to determine the shielding configurations of remote-handling tools and transport casks and to aid in planning maintenance activities. For this analysis, we utilized a hybrid unstructured mesh (UM)/constructive solid geometry approach for calculating spallation products and neutron fluxes, activation calculations using the AARE package that includes the CINDER2008 activation code to calculate the decay photon source at different cooling times, and the ADVANTG code to accelerate the final decay photon transport calculation. Both Type 316 stainless steel (SS-316) and lead were investigated as candidates for shielding materials. The decay photon transport calculation through the thick SS-316 or lead shields exhibited between 25 and 30 orders-of-magnitude attenuations in the radial direction, depending on the shield. Such a difficult shielding calculation required advanced variance reduction. ADVANTG has some missing features, which limits its usability in spallation neutron source applications. It does not support volumetric sources created for MCNP6.2 UM capability. An approximate source was created for this problem. Not only was this approximate source needed for running the ADVANTG calculation to generate the weight windows, but also it was essential to develop source biasing (SB) parameters that were crucial for dramatically accelerating the decay photon transport in this problem. With this approximate source, the analysis was completed in a very reasonable computational time, and the design of the STS remote-handling equipment was finalized. This paper compares the efficiency of Monte Carlo simulations with different weight window and SB parameters calculated using different approximate ADVANTG calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extension of OpenMC for Fixed Source Transmutation Calculations

This report documents work performed under a Strategic Partnership Project between Argonne National Laboratory (ANL) and the United Kingdom Atomic Energy Authority (UKAEA). The overall goal of this project is to extend OpenMC [1], a community-developed Monte Carlo particle transport code, to be able to perform fixed-source transmutation calculations. In fission and fusion reactors, the high flux of energetic neutrons causes materials within the reactor to “transmute,” or undergo a nuclear reaction that results in the addition/removal of neutrons and protons from the nucleus of an atom. If subjected to these reactions for long enough, the overall composition and physical properties of the material itself begin to change as a result of transmutation. Such a feature is vital for predicting the decrease in tritium production rate within a breeder blanket during the lifetime of a fusion reactor. OpenMC is capable of simulating neutron transport in fission/fusion systems, thereby allowing it to estimate the flux that causes transmutation. It is also capable of solving the transmutation equations, which determine how the composition of a material changes over time due to neutron irradiation and radioactive decay. However, solving the transmutation equations was previously only possible when the source of neutrons came from a fission system. In a fusion system, the source of neutrons is typically determined by a separate code and then given as an input to the particle transport simulation. This is known as a fixed source calculation. Through this project, we have extended OpenMC to solve the transmutation equations for a fixed source calculation. Evaluating the change in material compositions due to transmutation and its effect on physical properties is of key importance to a range of engineering analyses for fission and fusion systems. For example, in a fusion reactor, estimating the dose rate at different physical locations resulting from irradiated materials in the reactor allows designers to ensure that workers are not exposed to doses beyond applicable regulations. In order to properly dispose of irradiated materials, designers also need to estimate the radiotoxicity, which again relies on knowledge of the material composition at some future time. The specific tasks for this project that were agreed to between ANL and UKAEA were as follows: 1. Make changes and additions in the openmc.deplete and related modules in OpenMC to support transmutation calculations following a fixed source transport simulation. 2. Make necessary changes to OpenMC to model transmutation due to an arbitrary set of reactions needed for fusion applications. Use this new capability to generate a depletion chain file based on the TENDL nuclear data library. 3. Improve the openmc.deplete module in OpenMC to keep track of gases produced as a by-product of nuclear reactions during transmutation calculations. 4. Validate the new capabilities by carrying out fixed-source transmutation calculations on a suitable benchmark problem using OpenMC and a comparable Monte Carlo neutron transport code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification and Validation Tests of the Gamma Library of the ARC Software Package (Rev.1)

As part of the verification and validation (V&V) efforts of the Argonne Reactor Code (ARC) software system to support the Versatile Test Reactor (VTR) project of the U.S. Department of Energy (DOE), work has been performed to verify and validate the gamma (or photon) library of the multigroup cross section generation code MC 2 -3 of the ARC system, which is the essential part of the coupled neutron and gamma heating calculations to determine the heat generation rate or power distribution in the core. This V&V study was carried out in three steps: 1) review of the procedures and utility programs for generating the MC 2 -3 gamma library, 2) verification tests of cross sections and KERMA factors by comparing the total heat generation rate of coupled neutron and gamma calculation with that obtained with reaction Q values under the assumption of local gamma energy deposition and by comparing the prompt heat generation rate with that calculated with the MCNP6.2 code, and 3) validation tests of cross sections and KERMA factors by analyzing the ZPPR-15D gamma dose data measured with thermoluminescent dosimeters (TLD).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

The Oak Ridge National Laboratory (ORNL) Second Target Station (STS) neutron production facility is an accelerator driven pulsed neutron source that is currently being actively developed at ORNL. The neutrons are produced by proton-induced spallation reactions. A proton beam of 700 kW power is delivered to a spallation target in short, less than 1 µs long pulses, with 15 Hz repetition rate. The spallation target of ORNL STS is a rotating water-cooled tungsten target with tantalum cladding housed in a stainless-steel shroud. It is divided into 21 segments. These segments become highly activated due to spallation reactions or nuclei transmutation by the emitted neutrons. The radioactive nuclides continue to decay after ceasing operation. The decay dose rates generated from the target segments once they are removed from their operational location within the core vessel must be accurately quantified to determine the shielding configurations of remote handling tools and transport casks and to aid in planning maintenance events. To determine the shielding configurations needed for an activated target segment after ceasing operation, both the hybrid unstructured mesh (UM)/constructive solid geometry (CSG) approach that was previously utilized for STS analyses [1] and the ADVANTG code [2] were used. Even though the ADVANTG code does not include UM capability, the utilization of its advanced variance reduction technique was crucial to accelerate the extremely difficult final photon transport calculation in this analysis. This paper also describes the procedures taken to mitigate the convergence issues that occur when ADVANTG uses a source definition that does not match the source of the final Monte Carlo (MC) calculation. These convergence issues often occur because of inconsistencies between source and transport biasing parameters.

Ibrahim, Ahmad↗

Tank 11H Low Temperature Aluminum Dissolution and Inhalation Dose Potential Analyses at Savannah River Site – 26018

Currently, there is approximately 34 million gallons of high-level radioactive tank waste in the Tank Farm at the Savannah River Site (SRS). The ultimate goal of operations at the Tank Farm is to remove the high level waste (HLW) from the tanks followed by stabilization of the waste through vitrification of the HLW into glass or grouting the decontaminated waste into saltstone. After bulk removal of the HLW consisting of sludge, saltcake, and supernatant, further efforts are made to reduce the residual waste present in the tank in order to declare preliminary cease waste removal (PCWR) signifying completion of HLW removal. These reduction efforts can include tank washing to remove soluble salts and radioisotopes and dissolution of solids including aluminum. Aluminum in the form of gibbsite and boehmite is relatively insoluble in water. Through addition of aqueous sodium hydroxide, the aluminum can be dissolved at mild temperatures. In order for the waste tank to meet closure mode requirements of the Concentration, Storage, and Transfer Facilities (CSTF), which includes the Tank Farm, Documented Safety Analysis (DSA), a component of the safety basis, the inhalation dose potential (IDP) and the radiolytic hydrogen generation rate of the stored waste must be demonstrated to be lower than their respective designated limits. These parameters are calculated from measured radiochemical analyses of isotopes that emit a high amount of radioactivity including Cs-137, Sr-90, Pu-238, Pu-239, Pu-240, Pu-241, Am-241, and Cm-244. Following the low temperature aluminum dissolution (LTAD) process, Tank 11H slurry samples were pulled from the tank and sent to Savannah River National Laboratory (SRNL) to measure the extent of aluminum dissolution, hydroxide concentration, densities of slurry and supernatant, weight percent solids analyses, and radionuclide activities. The analyses of the composite sample found that approximately 90% of the total aluminum in the slurry was dissolved, indicating successful reduction of the insoluble aluminum in the waste tank. Additionally, the weight percent insoluble solids (slurry basis) measurement of the composite sample was found to be approximately 1%, demonstrating that minimal solids still remain in the tank. Finally, the radiochemical analyses of the composite sample determined that the waste contents of the tank met the IDP and radiolytic hydrogen generation rate requirements of the CSTF DSA. These measurements have shown that the LTAD process in Tank 11H was successful in waste reduction efforts and a positive step towards declaring PCWR and tank closure at SRS.

Dekarske, John [Savannah River National Laboratory↗

High-throughput screening tools facilitate calculation of a combined exposure-bioactivity index for chemicals with endocrine activity

High-throughput and computational tools provide a new opportunity to calculate combined bioactivity of exposure to diverse chemicals acting through a common mechanism. We used high throughput in vitro bioactivity data and exposure predictions from the U.S. EPA’s Toxicity and Exposure Forecaster (ToxCast and ExpoCast) to estimate combined estrogen receptor (ER) agonist activity of non-pharmaceutical chemical exposures for the general U.S. population. High-throughput toxicokinetic (HTTK) data provide conversion factors that relate bioactive concentrations measured in vitro (µM), to predicted population geometric mean exposure rates (mg/kg/day). These data were available for 22 chemicals with ER agonist activity and were estimated for other ER bioactive chemicals based on the geometric mean of HTTK values across chemicals. For each chemical, ER bioactivity across ToxCast assays was compared to predicted population geometric mean exposure at different levels of in vitro potency and model certainty. Dose additivity was assumed in calculating a Combined Exposure-Bioactivity Index (CEBI), the sum of exposure/bioactivity ratios. Combined estrogen bioactivity was also calculated in terms of the percent maximum bioactivity of chemical mixtures in human plasma using a concentration-addition model. Estimated CEBIs vary greatly depending on assumptions used for exposure and bioactivity. In general, CEBI values were <1 when using median of the estimated general population chemical intake rates, while CEBI were ≥1 when using the upper 95th confidence bound for those same intake rates for all chemicals. Concentration-addition model predictions of mixture bioactivity yield comparable results. Based on current in vitro bioactivity data, HTTK methods, and exposure models, combined exposure scenarios sufficient to influence estrogen bioactivity in the general population cannot be ruled out. Future improvements in screening methods and computational models could reduce uncertainty and better inform the potential combined effects of estrogenic chemicals.

54 ENVIRONMENTAL SCIENCES↗

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculator For Radioisotope Yields, Targets, And Products

An initial number of atoms in a material (by element or by isotope), defined by ENDF library availability, is used to produce an initial reaction products (RPs). The RPs are decayed using full, analytical solutions to the Bateman equations to calculate the RP inventory at a user-specified time. Since each radioisotope's inventory is known, the proportion to the total may be shown for each isotope. The total dose rate and activity is the main objective for experimenters. Built-in information about many reactors and irradiation positions is available. Full execution and display of results takes 5-10 seconds, depending on the complexity of the input parameters.

Holschuh, ThomasV [Idaho National Laboratory (INL)↗