Search NASA⌕ Search

Engineering topics

Peplow, Douglas E.

Publications and source records attributed to Peplow, Douglas E..

At least 19 records

Fractionation of UF6 and daughter progeny in storage cylinders from external heating

Abstract Uranium hexafluoride (UF 6 ) is a significant concern for material accountancy and verification in the international safeguards community. Verification of the contents of UF 6 cylinders is generally attempted with gamma spectroscopy but the current methods assume a uniform, homogeneous UF 6 mass distribution within the cylinder. In this work, it was found experimentally and confirmed via modeling, that under an external heat load (the sun), the UF 6 and its daughter products undergo fractionation in the cylinder. This fractionation of the UF 6 and daughter products can cause an errant measurement of the enrichment of the cylinder when using the current verification methods.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

SCALE 6.3 Validation: Radiation Shielding

Safe and reliable use of scientific and engineering computer codes requires validation for the types of applications in which they will be used. An example in the nuclear reactor engineering and licensing field is radiation transport employed in shielding analyses. The validity of computer codes for shielding applications is demonstrated in this report for SCALE version 6.3.0. Representative benchmarks corresponding to shielding analyses are selected for the validation study. Typical measurement results analyzed from these benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma dose rates, neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between measurement and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, the Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) radiation transport computer code with built-in variance reduction methods distributed with the SCALE computer code system agrees well with the measurement results. In selected benchmarks, MAVRIC is also compared to Monte Carlo N- Particle® (MCNP® ) 1 calculations. Both computer codes generally agree well within the estimated uncertainties. With the release of SCALE 6.3.0, Shift was integrated as an alternative transport solver in MAVRIC, denoted MAVRIC-Shift. Although the traditional MAVRIC using Monaco was used primarily in this validation study, many results have also been generated using MAVRIC-Shift. Agreement between MAVRIC-Monaco and MAVRIC-Shift is generally very good. The benchmarks presented in this report were obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook, the Shielding Integral Benchmark Archive & Database, and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to expand the shielding validation suite.

61 RADIATION PROTECTION AND DOSIMETRY↗

Threat Sources for Creating Synthetic Urban Search Data

Equivalent point source energy emission distributions were computed for various threat sources for use in simulating the detector responses for urban search scenarios. The sources include standard isotopic sources used in detector testing, medical and industrial sources occasionally encountered in urban searches, and several types of special nuclear materials. Most of the equivalent point source distributions represent small sources inside some amount of shielding, but the special nuclear material sources represent volumetrically distributed sources in spheres of metal. Text-based inputs for emission distributions are available for the Monte Carlo transport codes Monte Carlo N-Particle, SCALE/MAVRIC, and Omnibus/Shift, any of which can easily be converted to other formats. These sources were developed for use in the Radiological Anomaly Detection and Identification (RADAI) project and the follow-on project, the RADAI-Extended project, sponsored by the National Nuclear Security Administration Office of Defense Nuclear Nonproliferation Research and Development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SCALE 6.2.4 Validation: Radiation Shielding

For safe and reliable use of computer codes by the community, accuracy must be clearly evaluated. In particular, the nuclear reactor engineering and licensing field needs accurate tools for radiation shielding modeling. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is one such tool, with built-in variance reduction methods distributed within the SCALE code, and its validity is demonstrated in this report for the released version 6.2.4. Representative benchmarks corresponding to shielding analysis are selected for the validation study. Typical experimental results analyzed from those benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma doses, foil neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between experiment and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, MAVRIC agrees well with the experiment results. MAVRIC is also compared to Monte Carlo N-Particle (MCNP) calculations when available, and both codes generally produce good agreements within estimated uncertainties. The selected benchmarks are obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook (ICSBEP Handbook), the Shielding Integral Benchmark Archive & Database (SINBAD), and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to incrementally expand the shielding validation suite incrementally.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Monte Carlo Simulation of Background and Source Measurements with CSG and CAD Geometries

Detecting radioactive sources in an urban environment is difficult due to the large magnitude and variability of the background radiation. To support the search mission of the National Nuclear Security Administration, a project was undertaken to determine if first-principles modeling and simulation could be used to accurately predict the response of radiation detectors in urban environments. This study demonstrates that the simulated detector responses to photon radiation from both an isotopic source and the background compare well to benchmark-quality measurements in a large urban-like test environment. Simulation results using a traditional combinatorial solid geometry and a CAD model generated from LiDAR measurements gave very similar detector responses. With this validation, further simulations could be used to predict detector responses for various situations in real urban settings.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Skyshine Calculations for a Large Spent Nuclear Fuel Storage Facility with SCALE 6.2.3

The SCALE code system developed at Oak Ridge National Laboratory includes state-of-the-art capabilities for radiation source term and radiation transport simulations that can be used in numerous applications, including dose rate analyses of complex consolidated interim storage facilities (CISFs). A licensed CISF could be used to store tens of thousands of tonnes of spent nuclear fuel discharged from commercial power reactors using various cask and storage pad designs. A CISF design must comply with the regulatory requirements provided in 10 CFR Part 72, including requirements related to annual dose limits applicable to real individuals located beyond the area controlled by the licensee. Therefore, calculating a dose to the public is a necessary part of the licensing process for the construction of a CISF. These calculations are very challenging because of the complexity of the CISF design and the low magnitude of dose rate at large distances from the facility. This paper describes detailed far-field dose rate calculations performed for a proposed CISF using MAVRIC, the Monte Carlo radiation shielding sequence in SCALE 6.2.3, with automated variance reduction based on discrete ordinates calculations. The method presented in this paper uses a detailed Monte Carlo radiation transport simulation in one step from source to dose rate. A series of independent simulations was made using the complete site geometry (all casks present), but with only one cask containing radiation sources to obtain the dose rate maps produced by each storage cask. The CISF dose rate map was obtained by adding the dose rate maps produced by the independent individual cask simulations. Ample volumes of air and soil extending beyond the location of interest for dose rate calculation were included in the calculation model to properly simulate important radiation attenuation and scattering events that affect far-field dose rates. Furthermore, a comprehensive sensitivity study is included in this paper to illustrate the importance of selecting appropriate air volume, mass density, and composition for CISF skyshine dose rate calculations. Dry soil and soil containing water were analyzed to determine their effects on groundshine radiation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data for training and testing radiation detection algorithms in an urban environment

The detection, identification, and localization of illicit nuclear materials in urban environments is of utmost importance for national security. Most often, the process of performing these operations consists of a team of trained individuals equipped with radiation detection devices that have built-in algorithms to alert the user to the presence nuclear material and, if possible, to identify the type of nuclear material present. To encourage the development of new detection, radioisotope identification, and source localization algorithms, a dataset consisting of realistic Monte Carlo-simulated radiation detection data from a 2 in. × 4 in. × 16 in. NaI(Tl) scintillation detector moving through a simulated urban environment based on Knoxville, Tennessee, was developed and made public in the form of a Topcoder competition. The methodology used to create this dataset has been verified using experimental data collected at the Fort Indiantown Gap National Guard facility. Realistic signals from special nuclear material and industrial and medical sources are included in the data for developing and testing algorithms in a dynamic real-world background.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Generation of Synthetic Data for a Radiation Detection Algorithm Competition

This paper details the generation of synthetic radiation data using large-scale Monte Carlo transport models to evaluate radiation search detection algorithms. Modular 3-D Monte Carlo models spanning multiple city blocks were constructed, loosely based on downtown Knoxville, TN, containing buildings composed of multiple materials (brick, granite, and concrete), sidewalks, a four-lane road, side streets, parking lots, and grassy fields. Background and simulated source detector response calculations from these models were used to create synthetic list mode data sets for a 2"×4"×16" NaI(Tl) detector moving through a city street at a constant speed. For the background simulations, major isotopes were computed individually so that background composition and variability could be computed efficiently outside Monte Carlo. The source detector response included six simulated sources placed at 15 source locations. Detector response was developed to be periodic through the city street so that a detector path could begin at one end of the model and wrap around to the other. This framework allowed for the creation of diverse data sets, each with its own unique background and simulated source detector response. Synthetic data allows for high-quality labels, which are useful in developing data-driven radiation detection algorithms. This methodology was used to create synthetic data sets which were released as part of a public data competition to spur the development of new radiation detection algorithms for radiological search applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Specific Gamma-Ray Dose Constants with Current Emission Data

The specific gamma-ray dose constant represents the gamma effective dose rate due to a point source of unit activity of a given nuclide at 1 m. New tabulations of specific gamma-ray dose constants have been made using current gamma emission data from the SCALE 6.2.3 software package and International Commission on Radiological Protection Publication 107, combined with the effective dose per fluence conversion coefficients (antero-posterior orientation) of International Commission on Radiological Protection Publication 116. SCALE data cover 1,264 nuclides, and International Commission on Radiological Protection Publication 107 data include 1,192 nuclides, with only 777 nuclides in common between the two sets.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Radioisotope Power System Dose Estimation Tool (RPS-DET) User Manual

Radioisotope power systems (RPSs) convert decay heat from radioactive materials to electrical power. The National Aeronautics and Space Administration (NASA) uses RPSs for certain orbiters, landers, rovers, and probes that cannot rely on solar or chemical power. Most of the emitted charged particles in RPS fuels are locally absorbed, and the particles’ kinetic energy is converted into thermal energy as they slow down in the surrounding materials. This heat is subsequently converted into electricity through some form of energy conversion process (e.g., thermoelectric, heat engine, thermionic). While nearly all the charged particles are stopped locally, a non-trivial portion of the emitted neutral particles (neutrons and photons) escape the fuel and are free to interact with surrounding materials. This escaping radiation field can be operationally problematic for certain instruments, materials, electronics, etc., and can cause problems during terrestrial shipping, integration, testing for personnel, instrumentation, etc. Furthermore, the fuel’s emission spectrum evolves over time as the parent isotopes decay away and daughter products accumulate. This changing fuel composition alters characteristics of the ambient radiation field over time, which is further affected by RPS orientation and the surrounding environment. Tools that can rapidly estimate the radiation characteristics of RPSs with various fuel compositions, fuel ages, various RPS orientations, and in/on different spaceflight and planetary scenarios provide researchers with an opportunity to quickly estimate the radiation implications of a given RPS scenario. The RPS dose estimation tool (RPS-DET) provides engineers, scientists and analysts who are not trained in the art of particle transport software with a solution for quickly building, simulating, and analyzing the radiation impacts of various RPS designs in relevant spaceflight and terrestrial environments.

30 DIRECT ENERGY CONVERSION↗

A Directional Detector Response Function for Anisotropic Detectors

For a radiation detector that is not isotropic, a directional detector response is needed to accurately account for the variation in a detector’s behavior depending on the incoming particle direction. The concept of the detector response function has been extended to include particle direction using a set of pregenerated detector responses based on the orientation of the incoming radiation and the detector. This directional detector response function (DDRF) then can be applied to the flux and current tallies computed by a Monte Carlo simulation. Validation of the new approach has been done by comparing simulated count rates processed with the DDRF to measured count rates taken with a 5.08 × 10.16 × 40.64-cm NaI(Tl) detector. The comparisons show that the applied method produces good agreement with both background and source measurements with a 137 Cs source. Furthermore, separation of the detector response generation from Monte Carlo particle transport calculations provides greater flexibility in locating single or multiple detectors without any interference in the model and also enables simulation of various models using the same detector response without the need for generating additional detector responses if the same detector is being used.

42 ENGINEERING↗

Preliminary Validation of the Shift Monte Carlo Code for Fixed-Source Radiation Transport Problems

Simulations of neutron and photon transport using the Shift Monte Carlo radiation transport code are compared with experimental measurements and their corresponding benchmark simulations from several sets of experiments. Overall, Shift results match the calculations made by the benchmark teams quite well and match the measured values, which typically have large uncertainties, fairly well. A variety of attenuation/scattering problems are examined, as well as a streaming problem and a skyshine problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Attribution of gamma-ray background collected by a mobile detector system to its surroundings using panoramic video

Mobile radiation detector systems are important tools for detecting radiological and nuclear sources outside of regulatory control, but due to their mobility, they are subject to complex and varying backgrounds in most realistic operational scenarios. Recent work has found correlations between non-radiological contextual information and gamma-ray spectral features that can be used to decrease false alarm rates, however a more complete understanding of background source terms has been elusive. A measurement campaign developed a full radiological characterization of a controlled facility that roughly corresponded to two city blocks. As part of the campaign, the Radiological Multi-sensor Analysis Platform (RadMAP) collected extensive multi-sensor data. In this work, RadMAP’s panoramic video data were used to visually identify several different materials and quantify the detector response to each material. A linear model was fit between the material responses and radiological features, and the results demonstrate reasonable agreement with ground truth flux measurements for each material. This paper will describe the data collection, processing, and analysis of the gamma-ray and video data from RadMAP. The paper will conclude with perspectives on the applicability of such a method to less controlled environments both with respect to achieving better understanding the sources of variability of background radiation in urban environments and whether such methods could be leveraged in operational scenarios.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detailed SCALE Dose Rate Evaluations for a Consolidated Interim Spent Nuclear Fuel Storage Facility

This paper presents a confirmatory dose rate calculation for a proposed consolidated interim spent nuclear fuel (SNF) storage facility (CISF). The proposed CISF was modeled in detail using SCALE 6.2.1 to determine the radiation dose map of the surrounding areas of the facility and determine the location of the controlled area boundary based on 10 CFR 72.104 requirements. The spatial dose rate distribution was determined as the total dose rate contributions by each of the 467 storage casks comprising the CISF. A rectangular boundary outside of which no dose rate exceeds 25 mrem/yr has the dimensions of 1.25 km x 1.25 km. It denotes the smallest rectangle that could serve as the controlled area boundary. The results (dose map) show the unique capability of the SCALE code system for a high-resolution radiation dose simulation of a large SNF storage facility.

Radulescu, Georgeta↗

A Methodology for Determining the Concentration of Naturally Occurring Radioactive Materials in an Urban Environment

The detection of radioactive sources in an urban setting is greatly complicated by natural background radiation, which emanates from various materials including roadways, sidewalks, soil, and building exteriors. The method presented and demonstrated here represents an effort to characterize the concentration of naturally occurring radioactive material (NORM) in these types of materials. The location surveyed in this work was the Fort Indiantown Gap Combined Arms Collective Training Facility in Lebanon County, Pennsylvania. Over 70 measurements with a high-purity germanium detector were performed to ascertain the NORM concentrations present in the soil, asphalt, gravel, concrete, and walls found throughout the site. Monte Carlo radiation transport calculations were used to obtain detector responses for these various geometries and materials to convert these measurements into NORM concentration estimates. Lastly, synthetic spectra were simulated using the predicted source terms and compared to actual measurements, showing acceptable agreement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling Urban Scenarios & Experiments: Fort Indiantown Gap Data Collections Summary and Analysis

This report summarizes experimental radiation detector, contextual sensor, weather, and global positioning system (GPS) data collected to inform and validate a comprehensive, operational radiation transport modeling framework to evaluate radiation detector system and algorithm performance. This framework will be used to study the influence of systematic effects (such as geometry, background activity, background variability, environmental shielding, etc.) on detector responses and algorithm performance using synthetic time series data. This work consists of performing data collection campaigns at a canonical, controlled environment for complete radiological characterization to help construct and benchmark a high-fidelity model with quantified system geometries, detector response functions, and source terms for background and threat objects. This data also provides an archival, benchmark dataset that can be used by the radiation detection community. The data reported here spans four data collection campaigns conducted between May 2015 and September 2016.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Comparison of the FW-CADIS and MR-CADIS Variance Reduction Methods

The Consistent Adjoint Driven Importance Sampling (CADIS) method, which uses an adjoint deterministic estimate for making variance reduction parameters that optimize Monte Carlo (MC) shielding calculations, was introduced in 1998. This method was designed to optimize a single tally, typically an integrated response such as a detector count rate or a dose rate. Because both deterministic and stochastic methods are employed, CADIS is referred to as a hybrid method. In 2007, the forward-weighted CADIS (FW-CADIS) method was introduced to optimize several tallies at once, giving more uniform relative uncertainties among the tallies for both high-flux and low-flux areas of the problem. This method can be used for several small detectors, for semi-global problems (mesh tally over a large fraction of the problem), or for global problems (mesh tally over the entire problem). Optimization can be done for an integrated response such as a dose rate D(r) as a function of position, the total flux {phi}(r) as a function of position, or the space/energy flux {phi}(r,E). Both the CADIS and FW-CADIS methods have been in the MAVRIC sequence of SCALE for three versions (6.0 in 2009, 6.1 in 2011, and 6.2 in 2016). The ADVANTG package released in 2015 uses the CADIS and FW-CADIS methods to accelerate MCNP neutron/photon shielding calculations. In 2014, the multi-response CADIS (MR-CADIS) method was presented as a way to obtain even more uniform relative uncertainties for problems with multiple tallies than those obtained with FW-CADIS. The method does not use a forward deterministic calculation but instead uses one adjoint calculation for every tally to be optimized in the Monte Carlo calculation. This work compares the above methods.

42 ENGINEERING↗