Search NASA⌕ Search

SEARCH · Search NASA

Results for “mcnp”

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 253 records · Page 14

Yield and PDV Measurements of the Godiva-IV Critical Assembly

A high dynamic range, scinillator based detector, the MHD-240, was used to measure the gamma and neutron flux of several prompt critical bursts of the Godiva-IV assembly. The measurements provided the alpha and FWHM for several different burst temperatures. Calibrations of the MHD240 gamma and neutron sensitivity along with MCNP based simulation of the Godiva-IV assembly allowed us to also determine the total fission yield of the burst integrated over approximately 5 ms. In addition a Photo Doppler Velocimetry measurement was conducted to determine the surface motion vibrational frequency of Godiva-IV assembly during the bursts. The combined data were collected to be used as input to multi-physics simulation of the Godiva-IV assembly. This is a preliminary report, there is an ongoing effort to update the MCNP simulations used in the analysis and develop a correction for the room return component of the measured signal.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Yield and PDV Measurements of the Godiva-IV Critical Assembly (Preliminary Report)

A high dynamic range, scinillator based detector, the MHD-240, was used to measure the gamma and neutron flux of several prompt critical bursts of the Godiva-IV assembly. The measurements provided the alpha and FWHM for several different burst temperatures. Calibrations of the MHD240 gamma and neutron sensitivity along with MCNP based simulation of the Godiva-IV assembly will allow us to also determine the total fission yield of the burst integrated over approximately 5 ms. Work is ongoing to check the calibrations of the MHD-240, and improve the MCNP simulations. This is a preliminary report, and the fission yields reported here are not final and are likely inaccurate. The measured alpha and FWHM are not impacted by the calibration issues and are considered final.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

KIPT ADS Target Cooling Analyses by MCNP/Fluent

This report presents multi-physics analyses of the target cooling system of the accelerator driven system (ADS) of the Kharkiv Institute of Physics and Technology (KIPT) using MCNP and Fluent computer programs. MCNP has been used to transport electrons, gammas, and neutrons, and to calculate the energy deposition in the target using a mesh-tally. The mesh-tally data have been imported in Fluent by a C subroutine that has been compiled and linked to Fluent as a user defined function. The KIPT ADS is in operation and Fluent used the CAD model from the manufacturing process of the target. The Fluent results for the reference case match very well the literature results obtained by STAR-CCM+ during the design phase. Other cases that differ from the reference one have been analyzed; in these cases, it is assumed a malfunction of the electron accelerator control or the primary water coolant pump. The target cooling system operates normally for all the analyzed cases except when the inlet water mass flow rate is decreased. If the water coolant mass flow decreases to 25%, the water maximum temperature exceeds the boiling point. However, the transient analysis showed that the target cooling system can operate for 180 seconds when the inlet water mass flow rate is decreased down by 50%.

42 ENGINEERING↗

Performance Analysis of the Expected Track Length Estimator

This work investigates whether performance of the MCNP code is improved by the implementation of the expected track length estimator. The expected track length estimator is known to decrease variance when compared to the track length estimator in many classes of Monte Carlo particle transport simulations. However, the computational cost of computing an exponent for every particle collision within a tallied cell has previously been considered prohibitively expensive. Therefore, the goal of this work is to determine whether the variance improvements seen by the expected track length estimator in most cases makes up for the increase in transport time when implemented in the MCNP code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

EUCLID: Experiments Underpinned by Computational Learning for Improvements in Nuclear Data [Slides]

EUCLID will design validation experiments optimized to resolve compensating errors and adjust nuclear data to experiments. A big part of the success of the EUCLID proposal was due to previous work supported by NCSP (MCNP, nuclear data, and NCERC capabilities). The work performed under EUCLID will similarly benefit the NCSP mission. It will lead to new MCNP capabilities, improved nuclear data and nuclear data capabilities, and new methodology and tools that will have large impact on future NCERC experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SCALE Calculations Replicating ARH-600 Data

This report documents Standardized Computer Analyses for Licensing Evaluation (SCALE) calculations replicating data from the ARH-600 criticality handbook. These calculations complement those reproduced with MCNP in PRC-NS-00009. Calculations were performed with SCALE version 6.1, using the v7-238 cross section library. The original ARH-600 data were calculated by the Atlantic Richfield Hanford (ARH) Company starting in the 1960s, and Monte Carlo N-Particle (MCNP) calculations were performed by the Hanford Plateau Remediation Company in 2010. This report is not intended to replace final analysis of fissile systems by qualified criticality personnel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Installing MCNP6.2 on Microsoft Windows

This video shows how to install the MCNP code, version 6.2, on a Microsoft Windows 10 computer. After installation, the Xming software is installed, a quick calculation is performed, and geometry is plotted. At this point, this setup is ready for use in an MCNP classroom setting and/or for performing calculations of the viewer’s own creation.

97 MATHEMATICS AND COMPUTING↗

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗

NCERC-CL TREAT Free Field Characterization Foil Reaction Rate Results

This report summarizes the analysis and preliminary MCNP reaction rate modeling of free field characterization measurements performed in the Transient Reaction Test Facility (TREAT) at Idaho National Laboratory (INL) between 12/18/2023 and 2/22/2024. A series of seven irradiations were performed using the Big-BUSTER (Broad Use Specimen Transient Experiment Rig) core configuration. Foil sets including 19.5% LEU-Zr alloy wire, S, Au, Fe, Ni, Co, Ti, and Zr were deployed in each irradiation. Identical foil sets were supplied to Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and INL. All foils were supplied by INL. Additional details on the experiment can be found in [2]. The foil sets were shipped to the National Criticality Experiments Research Center Counting Laboratory (NCERC-CL) branch at LANL (NCERC-CL NISC) and were received on 3/25/2024. Additional details on receipt and counting are reported in [3]. All reactor dosimetry measurements were performed in adherence to the ASTM standards applicable to reactor dosimetry. Measured reaction rates from Au, Co, Fe, Ni, and Ti are reported. Select reaction rate ratios and simulated reaction rates and ratios are discussed. Reaction rates for all measured reaction products were modeled with MCNP 6.3.1 and a TREAT input deck supplied by Edward Lum. Neutron emission estimates calculated using the measured and simulated reaction rates show excellent agreement between the three fast-threshold n-p reactions with the percent differences being < 9% between the three reactions. The capture reactions had poorer, but still reasonable, agreement with < 20% percent differences between the neutron emission estimates of the three capture reactions. The percent difference between the fast-threshold and capture reactions is ∼ 100% for each reaction. This supports the TREAT model is not accurately modeling the neutron spectrum and additional measurements are required to characterize TREAT and match measurements to simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gamma and Neutron Measurement and Modeling of Irradiated TRISO Fuel

Given the unique characteristics of the PBR fuel cycle, both gamma and neutron measurements are expected to play important roles in performing and maintaining nuclear material control and accounting for spent pebbles to safeguard the fuel cycle. Given the lack of irradiated pebbles in the US, a variety of irradiated TRISO fuel samples with wide ranges of burnups and cooling times available at ORNL were used in this work. A large number of gamma and neutron measurements have been performed on these samples to collect data to test the various detectors and to benchmark the computer models to simulate the depletion and decay of the fuel and the measurements themselves. Two neutron detectors, including a custom-made detector and the Very High-Performance Neutron Multiplicity Counting, were used to measure the neutrons emitted by these TRISO samples. Three gamma spectrometry detectors, including an HPGe and the M400 CZT detector, were used to measure gamma-ray emissions from these samples. The M400 was recently adopted by the IAEA for fresh uranium measurements, but it was tested for spent fuel measurements prior to this project. Detailed MCNP models were developed to simulate these neutron and gamma measurements. Some GADRAS models were also developed to cross check the MCNP models for the gamma measurements. It was found challenging to perform neutron measurements in the hot cell due to the high background counts. Close agreements were observed between the simulated and measured neutron count rates in both detectors’ measurements of californium calibration sources. Both the HPGe and M400 detectors were able to measure the 604 and 662 keV peaks from these samples, which are the two most important peaks used to infer fuel burnup. Although the M400 detector did not have nearly good energy resolution and did not detect some of the minor peaks as the HPGe detector, it was found to be capable of handling significantly higher dose rates than HPGe. Given the complexities in the TRISO samples (e.g., different samples sizes) and uncertainties in the alignments between the detector and the TRISO fuel inside the containers, large scatters were found between the peak area rates and the samples’ burnups. However, the 604/662 peak ratios were found to trend well with the samples’ burnups among most samples in both measured and simulated results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rapid Gamma Simulations of TRISO Fuel Elements

As energy demand rises, nuclear energy, particularly from reactors that use tristructural isotropic (TRISO) fuels, has gained attention due to the fuel’s enhanced resistance to radiation damage and high temperatures. This report investigates the modeling capabilities of the Gamma Detector Response and Analysis Software (GADRAS) for TRISO fuels, focusing on the gamma signatures of TRISO particles, which have not been extensively explored. Using the Monte Carlo N-Particle (MCNP) code as a benchmark, we developed both homogeneous and heterogeneous models of TRISO pebbles to analyze gamma spectra. Our findings reveal that the homogeneous and heterogeneous models produced different gamma signatures. Additionally, the GADRAS heterogeneous model significantly reduces computation times compared to MCNP, enabling effective modeling of gamma signatures for safeguards applications. This advancement is essential for the International Atomic Energy Agency (IAEA) in detecting anomalies and potential smuggling attempts in TRISO reactor fuel elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rates for neutron-capture reactions on tungsten isotopes in iron meteorites

High-precision W isotopic analyses by Harper and Jacobsen indicate the W-182/W-183 ratio in the Toluca iron meteorite is shifted by -(3.0 +/- 0.9) x 10(exp -4) relative to a terrestrial standard. Possible causes of this shift are neutron-capture reactions on W during Toluca's approximately 600-Ma exposure to cosmic ray particles or radiogenic growth of W-182 from 9-Ma Hf-182 in the silicate portion of the Earth after removal of W to the Earth's core. Calculations for the rates of neutron-capture reactions on W isotopes were done to study the first possibility. The LAHET Code System (LCS) which consists of the Los Alamos High Energy Transport (LAHET) code and the Monte Carlo N-Particle(MCNP) transport code was used to numerically simulate the irradiation of the Toluca iron meteorite by galactic-cosmic-ray (GCR) particles and to calculate the rates of W(n, gamma) reactions. Toluca was modeled as a 3.9-m-radius sphere with the composition of a typical IA iron meteorite. The incident GCR protons and their interactions were modeled with LAHET, which also handled the interactions of neutrons with energies above 20 MeV. The rates for the capture of neutrons by W-182, W-183, and W-186 were calculated using the detailed library of (n, gamma) cross sections in MCNP. For this study of the possible effect of W(n, gamma) reactions on W isotope systematics, we consider the peak rates. The calculated maximum change in the normalized W-182/W-183 ratio due to neutron-capture reactions cannot account for more than 25% of the mass 182 deficit observed in Toluca W.

Masarik, J.↗

Assay of the Martian Regolith with Neutrons

The purpose of the research is to combine experiments and Monte Carlo transport of neutrons through volume of soil in an attempt to model neutron leakage from planetary surfaces. Emphasis is given to the change of neutron spectra as a function of water content and location. During the first stage of effort, two experiments were conducted in which leakage of neutrons from a Pu-Be source through about 30 g/cm(exp 2) of soil were measured with several counters. A Monte Carlo code, MCNP, has been used to model many of the 100 individual runs of the experiment. Hydrogen is the element that has the most dramatic effect on the neutron spectrum and its effect on the neutron spectrum is almost the same whether it is in the form of water or polyethylene. In order to simulate various water configurations, sheets of polyethylene have been used between layers of soil as well as water in several concentrations up to 18%. Comparison of experimental results to theoretical predictions made with the MCNP code were disappointing for low concentrations of water. We have made extensive calculations to see if room return could be the cause of the discrepancies. Water concentrations of the 'dry' soil were measured by two different laboratories and differed only by 0.5%. We have made calculations to optimize the next experiment and are investigating other methods of determining the water content of 'dry' soil.

Drake, Darrell M.↗

Measurement of Insulation Compaction in the Cryogenic Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

The liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC) use expanded perlite as thermal insulation. Th ere is evidence that some of the perlite has compacted over time, com promising the thermal performance and possibly also structural integr ity of the tanks. Therefore an Non-destructive Testing (NDT) method for measuring the perlite density or void fraction is urgently needed. Methods based on neutrons are good candidates because they can readil y penetrate through the 1.75 cm outer steel shell and through the ent ire 120 cm thickness of the perlite zone. Neutrons interact with the nuclei of materials to produce characteristic gamma rays which are the n detected. The gamma ray signal strength is proportional to the atom ic number density. Consequently, if the perlite is compacted then the count rates in the individual peaks in the gamma ray spectrum will i ncrease. Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. With commercially available portable neutron generators it is possible to produce simul taneously fluxes of neutrons in two energy ranges: fast (14 MeV) and thermal (25 meV). Fast neutrons produce gamma rays by inelastic scatt ering which is sensitive to Fe and O. Thermal neutrons produce gamma rays by radiative capture in prompt gamma neutron activation (PGNA) and this is sensitive to Si, AI, Na, Kand H. Thus the two energy ranges produce complementary information. The R&D program has three phases: numerical simulations of neutron and gamma ray transport with MCNP s oftware, evaluation of the system in the laboratory on test articles and finally mapping of the perlite density in the cryogenic tanks at KSC. The preliminary MCNP calculations have shown that the fast/therma l neutron NDT method is capable of distinguishing between expanded an d compacted perlite with excellent statistics.

Livingston, R. A.↗

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ISOCS-based Extended In Situ Gamma Spectrometry Services Tool SuperISOCS for Waste Measurements - 20275

Canberra's In Situ Object Counting System (ISOCS) is an established tool used for gamma spectrometry. It allows physical representation of complex geometries and mathematical calculation of the calibration function while avoiding the need for radioisotope standards. ISOCS provides a large range of templates to generate the geometry for most nuclear waste components in nuclear facilities. Nonetheless, there are some complex geometries where the accuracy attainable is limited by the conventional ISOCS templates. An ISOCS-based extended in-situ Gamma Spectrometry service tool has been developed, based on the so-called SuperISOCS (SISOCS) software. This tool was developed to generate the most complex geometries. The principle of geometry construction is based on the composition of 8 different primitive geometries. The primitive object is able to be cut and to change a part of an object. An object is able to overlap to another one to form a more complex object. These functions refine the modelled geometry compared to the real objects and decrease the measurement systematic error because they allow the mathematical model to represent the real-life geometry more faithfully through reducing the approximations and simplifications that need to be made The algorithm of SISOCS is the same as that of standard ISOCS. A set of calculations was done in comparison with MCNP efficiency calculation and SISOCS results. 80 different models were compared between SISOCS and MCNP, with added 30 cases compared between ISOCS and SISOCS. The average differences were less than 2%, for all energy's efficiency calculations. In this paper we illustrate the benefits that this new tool can provide, in terms of improved accuracy for complex realistic geometries, through simulation of some example scenarios. Furthermore, we demonstrate, through SISOCS modelling of the real object data test, the comparison between SISOCS to ISOCS calculations for the same geometry. This paper describes the added functions that are easily employed using SISOCS. Our results show the potential application and benefit of this technique for real waste assay projects and the way to integrate the tool into the waste assay systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗