Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Molecular dynamics simulation of thermal transport in UO 2 containing uranium, oxygen, and fission-product defects

Uranium dioxide (UO 2 ) is the most commonly used fuel in light-water nuclear reactors and thermal conductivity controls the removal of heat produced by fission, thereby governing fuel temperature during normal and accident conditions. The use of fuel performance codes by the industry to predict operational behavior is widespread. A primary source of uncertainty in these codes is thermal conductivity, and optimized fuel utilization may be possible if existing empirical models are replaced with models that incorporate explicit thermal-conductivity-degradation mechanisms during fuel burn up. This approach is able to represent the degradation of thermal conductivity due to each individual defect type, rather than the overall burn-up measure typically used, which is not an accurate representation of the chemical or microstructure state of the fuel that actually governs thermal conductivity and other properties. To generate a mechanistic thermal conductivity model, molecular dynamics (MD) simulations of UO 2 thermal conductivity including representative uranium and oxygen defects and fission products are carried out. These calculations employ a standard Buckingham-type interatomic potential and a potential that combines the many-body embedded-atom-method potential with Morse-Buckingham pair potentials. Potential parameters for UO 2+x and ZrO 2 are developed for the latter potential. Physical insights from the resonant phonon-spin-scattering mechanism due to spins on the magnetic uranium ions are introduced into the treatment of the MD results, with the corresponding relaxation time derived from existing experimental data. High defect scattering is predicted for Xe atoms compared to that of La and Zr ions. Uranium defects reduce the thermal conductivity more than oxygen defects. For each defect and fission product, scattering parameters are derived for application in both a Callaway model and the corresponding high-temperature model typically used in fuel-performance codes. The model is validated by comparison to low-temperature experimental measurements on single-crystal hyperstoichiometric UO 2+x samples and high-temperature literature data. Furthermore, this work will enable more accurate fuel-performance simulations and will extend to new fuel types and operating conditions, all of which improve the fuel economics of nuclear energy and maintain high fuel reliability and safety.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Summary report on UO 2 thermal conductivity model refinement and assessment studies

Uranium dioxide (UO 2 ) is the most commonly used fuel in light water nuclear reactors and thermal conductivity controls the removal of heat produced by fission, therefore, governing fuel temperature during normal and accident conditions. The use of fuel performance codes by the industry to predict operational behavior is widespread. A primary source of uncertainty in these codes is thermal conductivity, and optimized fuel utilization may be possible if existing empirical models were replaced with models that incorporate explicit thermal conductivity degradation mechanisms during fuel burn-up. This approach is able to represent the degradation of thermal conductivity due to each individual defect type, rather than the overall burn-up measure typically used which is not an accurate representation of the chemical or microstructure state of the fuel that actually governs thermal conductivity and other properties. To generate a mechanistic thermal conductivity model, molecular dynamics (MD) simulations of UO 2 thermal conductivity including representative uranium and oxygen defects and fission products are carried out. These calculations employ a standard Buckingham type interatomic potential and a potential that combines the many-body embedded atom method potential with Morse-Buckingham pair potentials. Potential parameters for UO 2+x and ZrO 2 are developed for the latter potential. Physical insights from the resonant phonon-spin scattering mechanism due to spins on the magnetic uranium ions have been introduced into the treatment of the MD results, with the corresponding relaxation time derived from existing experimental data. High defect scattering is predicted for Xe atoms compared to that of La and Zr ions. Uranium defects reduce the thermal conductivity more than oxygen defects. For each defect and fission product, scattering parameters are derived for application in both a Callaway model and the corresponding high-temperature model typically used in fuel performance codes. The model is validated by comparison to low-temperature experimental measurements on single crystal hyper-stoichiometric UO 2+x samples and high-temperature literature data. Ongoing works include investigation of the effect of phase separation to UO 2 +U 4 O 9 on the low temperature thermal conductivity of UO 2+x , and modeling of thermal conductivity using the Green-Kubo method. Ultimately, this work will enable more accurate fuel performance simulations as well as extension to new fuel types and operating conditions, all of which improve the fuel economics of nuclear energy and maintain high fuel reliability and safety.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Umpire: Integration with Other Software Technologies

Umpire is an open-source C++ library that provides memory and resource management functions for use with heterogeneous memory spaces. The Umpire library is being developed to simplify coordination and sharing of limited memory resources for complex integrated applications developed in the ECP/ATDM and ASC programs, and other multi-physics codes that combine many physics packages and libraries. The overarching goal of Umpire is to provide a single application-facing API for memory management on all platforms by decoupling memory operations from platform-specific memory spaces. It provides operations for: memory space discovery, memory provisioning, allocation and transfers, and introspection. Umpire development began in the Center for Applied Scientific Computing at Lawrence Livermore National Laboratory (LLNL). Umpire is funded primarily by LLNL ASC and ATMD portfolios and receives a small amount of support as part of ECP ST Project 2.3.1.06 STMP08-ISC. The source code is hosted and developed on GitHub. In this report, we explore possible integrations with other memory hierarchy management software. We examine both projects supported by ECP, as well as third-party and vendor provided tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

A systematic study of magnetic braking in low-mass binaries

A short summary is given of Rappaport et al. (1983) which described results of extending a simplified stellar evolution code covering the evolution of low-mass compact binaries. Magnetic braking is probably an important process in the evolution of such binaries (such as cataclysmic variables and low-mass X-ray sources). The initial simplified code describes the mass-losing star as an n = 3/2 polytrope and was developed to study the evolution of binaries with a secondary of low mass (between 0.01 and 0.4 solar mass) when the angular momentum losses are due to gravitational radiation. In the extended code, a composite polytrope model is used for the secondary, wherein the structure of the radiative core is described by an n = 3 polytrope and the convective envelope by an n = 3/2 polytrope.

Verbunt, F.↗
DOE OSTI2016

SU-F-T-54: Determination of the AAPM TG-43 Brachytherapy Dosimetry Parameters for A New Titanium-Encapsulated Yb-169 Source by Monte Carlo Calculations

Purpose: To determine the AAPM TG-43 brachytherapy dosimetry parameters of a new titanium-encapsulated Yb-169 source designed to maximize the dose enhancement during gold nanoparticle-aided radiation therapy (GNRT). Methods: An existing Monte Carlo (MC) model of the titanium-encapsulated Yb-169 source, which was described in the current investigators’ published MC optimization study, was modified based on the source manufacturer’s detailed specifications, resulting in an accurate model of the titanium-encapsulated Yb-169 source that was actually manufactured. MC calculations were then performed using the MCNP5 code system and the modified source model, in order to obtain a complete set of the AAPM TG-43 parameters for the new Yb-169 source. Results: The MC-calculated dose rate constant for the new titanium-encapsulated Yb-169 source was 1.05 ± 0.03 cGy per hr U, indicating about 10% decrease from the values reported for the conventional stainless steel-encapsulated Yb-169 sources. The source anisotropy and radial dose function for the new source were found similar to those reported for the conventional Yb-169 sources. Conclusion: In this study, the AAPM TG-43 brachytherapy dosimetry parameters of a new titanium-encapsulated Yb-169 source were determined by MC calculations. The current results suggested that the use of titanium, instead of stainless steel, to encapsulate the Yb-169 core would not lead to any major change in the dosimetric characteristics of the Yb-169 source, while it would allow more low energy photons being transmitted through the source filter thereby leading to an increased dose enhancement during GNRT. Supported by DOD/PCRP grant W81XWH-12-1-0198 This investigation was supported by DOD/PCRP grant W81XWH-12-1- 0198.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

MCNP6 Elevation Scaling of Cosmic Ray Backgrounds - Paper 29

The galactic cosmic-ray (GCR) source option was implemented in the all-energy, all-particle transport code MCNP62 in 2010. This source option along with other newly added features to MCNP6 have been used to produce increasingly accurate global cosmic and terrestrial background flux data files (background.dat).3 This data file is produced on a grid and is read and sampled by MCNP6 whenever a user invokes the background source option. The most recent adaptation of this data file, Version 4, will be completed in 2014. This paper reports on current progress toward enhancements that have been made to the modeling and simulation of these background spectra, in particular the elevation scaling used to produce more accurate spectra at off-grid locations. Cosmic radiation continuously bombards Earth with various particles, such as protons, alpha particles, and heavier nuclei, some of which are deflected by the Earth's shielding magnetic field. Particles that carry sufficient momentum can overcome the deflection and penetrate into the atmosphere. The sufficient momentum is dependent on the terrestrial coordinates due to the shape of the Earth's magnetic field and the Lorentz force's proportionality to the sine of the angle between the velocity vector of the incoming particle and the magnetic field direction. As the particles propagate through the atmosphere, collisions with atmospheric molecules generate secondary particles such as neutrons, protons, photons, muons, pions, and other exotic particles. These secondary particles often have sufficient energy to undergo additional nuclear interactions, and so on, forming what is known as a cascade shower. The particles in the cascade shower are constantly undergoing energy losses due to interactions, especially charged particles whose paths are bent into helical trajectories around magnetic field lines. These energy losses make the integral flux of the cosmic background at the ground level highly dependent on the atmospheric depth traveled by the incident particle and its progeny, and therefore dependent on the elevation of the ground where the cascade deposits its energy. For relatively small depth changes of approximately 100 g/cm{sup 2} or less, the ground-level cosmic flux follows an exponential model. The tabulation of background particle fluxes on the surface of the earth is important for a variety of reasons, one of which is the design of nuclear material detection systems. This method is also useful in producing user-specified effective attenuation lengths for particles. Attenuation lengths are an essential parameter used to compare particle spectra from different locations, elevations, and energy ranges. Traditionally these attenuation lengths are inferred from measured data. Now with the ability to simulate an accurate galactic cosmic ray source and transport it and its shower through a modeled atmosphere, these attenuation lengths can be calculated for any set of conditions. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRSDate not supplied

An Open-Source Numerical Model for Mitigating Refractory Alloy Hot Cracking Susceptibility

Refractory alloys are susceptible to solidification cracking during welding and 3D printing. Composition control is an effective method of controlling solidification cracking. This work evaluates the effect of compositional variation in refractory metal systems on a computed solidification cracking susceptibility. A numerical model has been developed using Python code and open-source CALPHAD software to calculate Kou’s crack susceptibility index. The model is validated against past weldability studies performed on several refractory alloy systems. The approach is extended towards the development of new alloys with improved 3D printability and weldability and is shown to have utility in defining compositional limits for existing alloys and feedstocks. Furthermore, the model will aid in determining process controls for powder reuse and recycling.

pycalphad↗
NASA NTRS2012

SE-FIT

The mathematical theory of capillary surfaces has developed steadily over the centuries, but it was not until the last few decades that new technologies have put a more urgent demand on a substantially more qualitative and quantitative understanding of phenomena relating to capillarity in general. So far, the new theory development successfully predicts the behavior of capillary surfaces for special cases. However, an efficient quantitative mathematical prediction of capillary phenomena related to the shape and stability of geometrically complex equilibrium capillary surfaces remains a significant challenge. As one of many numerical tools, the open-source Surface Evolver (SE) algorithm has played an important role over the last two decades. The current effort was undertaken to provide a front-end to enhance the accessibility of SE for the purposes of design and analysis. Like SE, the new code is open-source and will remain under development for the foreseeable future. The ultimate goal of the current Surface Evolver Fluid Interface Tool (SEFIT) development is to build a fully integrated front-end with a set of graphical user interface (GUI) elements. Such a front-end enables the access to functionalities that are developed along with the GUIs to deal with pre-processing, convergence computation operation, and post-processing. In other words, SE-FIT is not just a GUI front-end, but an integrated environment that can perform sophisticated computational tasks, e.g. importing industry standard file formats and employing parameter sweep functions, which are both lacking in SE, and require minimal interaction by the user. These functions are created using a mixture of Visual Basic and the SE script language. These form the foundation for a high-performance front-end that substantially simplifies use without sacrificing the proven capabilities of SE. The real power of SE-FIT lies in its automated pre-processing, pre-defined geometries, convergence computation operation, computational diagnostic tools, and crash-handling capabilities to sustain extensive computations. SE-FIT performance is enabled by its so-called file-layer mechanism. During the early stages of SE-FIT development, it became necessary to modify the original SE code to enable capabilities required for an enhanced and synchronized communication. To this end, a file-layer was created that serves as a command buffer to ensure a continuous and sequential execution of commands sent from the front-end to SE. It also establishes a proper means for handling crashes. The file layer logs input commands and SE output; it also supports user interruption requests, back and forward operation (i.e. undo and redo), and others. It especially enables the batch mode computation of a series of equilibrium surfaces and the searching of critical parameter values in studying the stability of capillary surfaces. In this way, the modified SE significantly extends the capabilities of the original SE.

Chen, Yongkang↗
NASA NTRS2008

Initial Integration of Noise Prediction Tools for Acoustic Scattering Effects

This effort provides an initial glimpse at NASA capabilities available in predicting the scattering of fan noise from a non-conventional aircraft configuration. The Aircraft NOise Prediction Program, Fast Scattering Code, and the Rotorcraft Noise Model were coupled to provide increased fidelity models of scattering effects on engine fan noise sources. The integration of these codes led to the identification of several keys issues entailed in applying such multi-fidelity approaches. In particular, for prediction at noise certification points, the inclusion of distributed sources leads to complications with the source semi-sphere approach. Computational resource requirements limit the use of the higher fidelity scattering code to predict radiated sound pressure levels for full scale configurations at relevant frequencies. And, the ability to more accurately represent complex shielding surfaces in current lower fidelity models is necessary for general application to scattering predictions. This initial step in determining the potential benefits/costs of these new methods over the existing capabilities illustrates a number of the issues that must be addressed in the development of next generation aircraft system noise prediction tools.

Nark, Douglas M.↗
DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

SU-E-T-380: Evaluation of BEBIG HDR 60Co System for AccuBoost Therapy

Purpose: In this project, the possibility of utilizing the BEBIG 60Co HDR system for AccuBoostTM treatment has been evaluated. Methods: Dose distributions in various breast sizes have been calculated for both Co-60 and Ir-192 sources using the MCNP5 code. These calculations were performed in breast tissues with thicknesses of 4cm, 6cm, and 8cm. The initial calculations were performed with the same applicator dimensions as the existing applicators used with the HDR Ir-192 system. The activity of the Co-60 source was selected such that the dose at the breast center was the same as the values from 192Ir. Then, the applicator thicknesses were increased to twice of those used with HDR Ir-192 system, for reducing skin and chest doses by Co-60 system. Dose to breast skin and chest wall were compared for both applicators types, with and without inclusion of a focusing cone at the applicator center. Results: The results showed that loading HDR Co-60 source inside the thin applicators impose higher doses to breast skin and chest wall compared to the 192Ir source. The area of the chest wall covered by 10Gy when treated by Co-60 with the thin and thick applicators, or treated by Ir-192 with thin applicator are 79cm2, 39cm2, and 3.8cm2, respectively. These values are reduced to 34cm2, 0cm2, and 0cm2 by using the focusing cone. It is worth noting that the breast skin areas covered by the 60Gy isodose line are 9.9cm2 and 7.8cm2 for Co-60 with the thin and thick applicators, respectively, while it is 20cm2 for Ir-192 when no focusing cone is present. These values are 0cm2, 0cm2, and 11cm2 in the presence of the focusing cone. Conclusion: The results indicate that using Co-60 with the thicker applicators is beneficial because of the higher half-life of Co-60, and the reduced maximum skin dose when compared with Ir-192.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Current State of Neutron Source Term Calculations

SOURCES-4C has some significant deficiencies that are mostly due to problems with the default data libraries that it uses. Some research groups have modified these libraries to fix most of the problems, but not all of them. Adding modern data libraries to SOURCES-4C is one of the most promising approaches to resolving these issues. Other promising options include using either NEDIS or a GEANT4 tool. Both of these codes approximately match the accuracy of the modified SOURCES-4C code. The NEDIS code is Russian. The GEANT4 tool might be computationally slow and it might also no longer be actively maintained. There are a handful of other (α,n) codes, but they are not accurate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Fringe biasing: A variance reduction technique for optically thick meshes

Fringe biasing is a stratified sampling scheme applicable to Monte Carlo thermal radiation transport codes. The thermal emission source in optically thick cells is partitioned into separate contributions from the cell interiors (where the likelihood of the particles escaping the cells is virtually zero) and the 'fringe' regions close to the cell boundaries. Thermal emission in the cell interiors can now be modelled with fewer particles, the remaining particles being concentrated in the fringes so that they are more likely to contribute to the energy exchange between cells. Unlike other techniques for improving the efficiency in optically thick regions (such as random walk and discrete diffusion treatments), fringe biasing has the benefit of simplicity, as the associated changes are restricted to the sourcing routines with the particle tracking routines being unaffected. This paper presents an analysis of the potential for variance reduction achieved from employing the fringe biasing technique. The aim of this analysis is to guide the implementation of this technique in Monte Carlo thermal radiation codes, specifically in order to aid the choice of the fringe width and the proportion of particles allocated to the fringe (which are interrelated) in multi-dimensional simulations, and to confirm that the significant levels of variance reduction achieved in simulations can be understood by studying the behaviour for simple test cases. The variance reduction properties are studied for a single cell in a slab geometry purely absorbing medium, investigating the accuracy of the scalar flux and current tallies on one of the interfaces with the surrounding medium. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data

Background: Microbiomes are now recognized as the main drivers of ecosystem function ranging from the oceans and soils to humans and bioreactors. However, a grand challenge in microbiome science is to characterize and quantify the chemical currencies of organic matter (i.e., metabolites) that microbes respond to and alter. Critical to this has been the development of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS), which has drastically increased molecular characterization of complex organic matter samples, but challenges users with hundreds of millions of data points where readily available, user-friendly, and customizable software tools are lacking. Results: Here, we build on years of analytical experience with diverse sample types to develop MetaboDirect, an open-source, command-line-based pipeline for the analysis (e.g., chemodiversity analysis, multivariate statistics), visualization (e.g., Van Krevelen diagrams, elemental and molecular class composition plots), and presentation of direct injection high-resolution FT-ICR MS data sets after molecular formula assignment has been performed. When compared to other available FT-ICR MS software, MetaboDirect is superior in that it requires a single line of code to launch a fully automated framework for the generation and visualization of a wide range of plots, with minimal coding experience required. Among the tools evaluated, MetaboDirect is also uniquely able to automatically generate biochemical transformation networks (ab initio) based on mass differences (mass difference network-based approach) that provide an experimental assessment of metabolite connections within a given sample or a complex metabolic system, thereby providing important information about the nature of the samples and the set of microbial reactions or pathways that gave rise to them. Finally, for more experienced users, MetaboDirect allows users to customize plots, outputs, and analyses. Conclusion: Application of MetaboDirect to FT-ICR MS-based metabolomic data sets from a marine phage-bacterial infection experiment and a Sphagnum leachate microbiome incubation experiment showcase the exploration capabilities of the pipeline that will enable the research community to evaluate and interpret their data in greater depth and in less time. It will further advance our knowledge of how microbial communities influence and are influenced by the chemical makeup of the surrounding system. The source code and User’s guide of MetaboDirect are freely available through (https://github.com/Coayala/MetaboDirect) and (https://metabodirect.readthedocs.io/en/latest/), respectively.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Gamma-ray imaging with a rotating hexagonal uniformly redundant array

Laboratory experiments have been performed to demonstrate the capabilities of a gamma-ray imaging system employing a NaI Anger camera and a rotating coded aperture mask. The mask incorporates in its design a new type of hexagonal uniformly redundant array (HURA) which is essentially antisymmetric under 60 deg rotation. The image formation techniques are described and results are presented that demonstrate the imaging capability of the system for individual and multiple point sources of gamma-ray emission. The results are compared to analytical predictions for the imaging and point source localization capabilities of coded aperture systems using continuous detectors.

Cook, W. R.↗
DOE OSTI2019

Exploiting fission chain reaction dynamics to image fissile materials

Radiation imaging is one potential method to verify nuclear weapons dismantlement. In this work, we present a method to discriminate between non-multiplying and multiplying neutron sources using a neutron coded aperture imaging system. This method applies time-correlated pulse-height (TCPH) analysis to identify neutrons in fission chain reactions and recreate images using only these events. This analysis was applied to measurements of weapons-grade plutonium (WGPu) metal and 252 Cf performed at the Nevada National Security Site (NNSS) Device Assembly Facility in July 2015. The results demonstrate it is possible to eliminate the non-fissile 252 Cf source from the image while preserving the fissile WGPu source. TCPH analysis was also applied to additional scenes in which the WGPu and 252 Cf sources were measured individually. The results of these separate measurements further demonstrate the ability to remove the non-fissile 252 Cf source and retain the fissile WGPu source. Simulations performed using MCNPX-PoliMi indicate that in a one hour measurement, hollow WGPu spheres are retained at a 1σ level for neutron multiplications M ≃ 2:7 and above. This work has potential application to dismantlement verification where it may be valuable to verify that all neutron sources in an image are multiplying.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗