Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MCNP modeling”

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 217 records · Page 12

Project Description for Student Symposium 2022: Modeling Approach to Critical in the Upcoming CERBERUS Experiment

The Critical Experiment Reflected By copper to BettEr Understand Scattering [CERBERUS] is an integral experiment that will be executed at the end of FY23. The experiment is designed to study the intermediate energy range, in particular elastic and inelastic scattering of neutrons in copper. Approaching criticality safely is of utmost importance, and having a computer mode is a useful tool when determining changes that can safely be made to the experiment. Multiple simulations are needed to analyze the approach to critical, and manually creating MCNP input files is both tedious and error prone. To solve this problem, a python script was created to develop and run these input files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Polyethylene density uncertainty in the simulation of neutron detectors

Simulations with the Monte Carlo N-Particle® (MCNP®) code are frequently used to study neutron coincidence counters for international safeguards. Reducing sources of uncertainty enables more accurate simulations, which reduces the need for expensive or impossible measurements. Accurately specifying the density of polyethylene is crucial for accurate simulations. The density of polyethylene ranges from 0.88 to 0.97 g/cm 3 and high-density polyethylene (HDPE) ranges from about 0.944 to 0.965 g/cm 3 . The density of a specific detector should be measured exactly for the best simulations. Here, the density of the blue/white Active Well Coincidence Counter (AWCC) bulk polyethylene was measured to be 0.9594 +- 0.0013 g/cm 3 at 1-sigma confidence. The density was calculated from repeated measurements of the various polyethylene dimensions for volume, and weighing the polyethylene for the mass. The AWCC neutron detection efficiency was simulated for each polyethylene density. The relative range of efficiencies (3σ) was 3.4% for unspecified polyethylene, 0.9% for HDPE, and 0.3% for the measured density. The impact of modeling small features such as screws was also studied and was found to be a negligible 0.2%. Measuring polyethylene density can reduce its role in uncertainty from one of the largest to one of the smallest contributors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Photoneutron Production Using an Electron Linear Accelerator for Applications in Neutron Imaging

Photoneutron production is possible using an electron linear accelerator and a target capable of generating photonuclear reactions. A short pulse neutron source can be useful for neutron imaging dynamic experiments. This study is aimed at the feasibility of photoneutron production using a 20 MeV electron linear accelerator and tungsten and depleted uranium targets of various thicknesses. MCNP6 (Monte Carlo N-Particle) code will be used to develop a computational model to estimate total neutron yield, and this will be verified at the Idaho State University’s Accelerator Center. After verification of the neutron yield and energy spectra, an additional MCNP6 model will be developed to analyze the neutron imaging processes. This study will potentially prove it is possible to conduct multi-mode imaging experiments on the anticipated Scorpius electron linear accelerator at the Nevada National Security Sites.

43 PARTICLE ACCELERATORS↗

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than \building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of MCNP, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Novel Algorithm for CAD to CSG Conversion in McCAD

Modeling and simulation lie at the heart of the design process of any nuclear application. An accurate representation of the radiation environment ensures not only the feasibility of new technologies, but it also aids in operation, maintenance, and even decommissioning. With increasingly complex designs, high-fidelity models have become a necessity for design maturity. McCAD has been under development for many years at Karlsruhe Institute of Technology (KIT) to facilitate the process of generating suitable models for nuclear analyses. In this paper, an overview of the major advances in the new version of the code is presented. A novel conversion algorithm has proven to be robust in significantly reducing the processing time to generate radiation transport models, making it easier to iterate on design details. A first-of-a-kind capability to generate hierarchical void cells is also discussed with preliminary analysis showing performance gains for particle tracking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Total Measurement Uncertainty in Neutron Coincidence Multiplicity Analysis

Neutron multiplicity counting is the most commonly used nondestructive assay technique for determining the plutonium mass within containers of scrap PuO 2 or mixed oxide (MOX). In multiplicity analysis, the 240 Pu eff mass, leakage multiplication, and alpha ratio (the ratio of [α, n]-to-spontaneous fission neutron production) are the three primary unknown sample properties. They must be determined simultaneously. To solve for these three unknowns in a multiplicity assay, three measured values are needed: the singles, doubles, and triples neutron count rates. While the analysis is limited to solving for three unknowns, there are many additional factors that impact the observed count rates and contribute to the measurement uncertainty. In this study we investigate the various uncertainty contributors for the multiplicity analysis through a combination of traditional uncertainty propagation techniques supplemented by Monte Carlo simulations to address the dependences not explicitly expressed by the point source model. Uncertainties arising from counting statistics, calibration parameters, calibration method, nuclear data, and various material characteristics (isotopic abundances, chemical form, density, and impurities) are considered. A Total Measurement Uncertainty (TMU) estimate is then developed from these uncertainty contributors. This study is confined to multiplicity analysis of items commonly encountered in international safeguards applications. That is, the study focused on Pu oxides and MOX materials for the masses ranging up to 4000 grams total Pu. Multiplicity measurements were simulated using MCNP V6 based on the Plutonium Scrap Multiplicity Counter (PSMC), Epithermal Multiplicity Counter (ENMC), Pyrochemical Multiplicity Counter, and Large Epithermal Multiplicity Counter (LEMC) for this study; however, this report focuses on the parameterization of the uncertainties for the PSMC. The performance differences between the PSMC and the other multiplicity counting systems are relatively small, primarily manifesting in the impact on measurement precision so that the evaluation developed for the PSMC can be applied to the other multiplicity counting systems. Finally an analysis tool, the Multiplicity TMU Estimator, was developed from this study to serve as an aid for evaluation of the total measurement uncertainty of multiplicity assay results obtained from the commonly used INCC acquisition and analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Characterization of the Fast-Neutron Irradiator and the Fast-Flux Tube Irradiation Fixtures at the Pennsylvania State Breazeale Reactor

Accurate knowledge of the neutron spectrum at a nuclear research reactor is a prerequisite for planning irradiation experiments, as well as for evaluating irradiation exposure results. The neutron-flux spectrum in the fast-neutron irradiator (FNI) and the fast-flux tube (FFT) irradiation fixtures at the Pennsylvania State Breazeale Reactor (PSBR) were characterized using the multi-foil neutron activation method. These irradiation fixtures make use of graded shielding to produce unique neutron fields. Multiple foil sets were irradiated in the fixtures with different exposure times and reactor powers to understand the stability over a wide range of operating conditions. Measured results were evaluated against a MCNP6 simulation to produce a measurement-informed neutron flux-energy spectrum for each fixture using STAYSL_PNNL. Simulated estimates of the FNI fixture, a newer fixture (~25 years old), demonstrated excellent agreement with measured results; the FFT did not. Thermal neutron measurements from the FFT suggest there is additional thermal leakage not captured in the simulation model. In conclusion, possible explanations for the discrepancy include burn-out or degradation (i.e., micro-cracking, gaps, etc.) in the boral and cadmium liners over the lifetime of the fixture (~40 years old).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of material accountancy techniques for 233 Pa from thorium nuclear fuels

Thorium is a promising alternative to uranium as nuclear fuel with advantages such as higher abundance, lower production of long-lived transuranic elements, and potentially better proliferation resistance. However, thorium presents a potential pathway for proliferation where produced 233 Pa can be diverted for the clandestine production of safeguarded 233 U. To prevent this, the ability to detect and measure 233 Pa must be assessed. This paper reviews several nuclear material accountancy techniques to determine their suitability for detecting 233 Pa extracted from irradiated thorium fuel. Hybrid K-edge densitometry and passive gamma spectroscopy have been found to be the best options based on technology maturity, cost, accuracy, and acquisition time. Thorium can be used in various reactor designs such as pressurized water reactors (PWRs), Canada deuterium uranium (CANDU) reactors, and molten salt reactors (MSRs). Therefore, thorium-uranium oxide fueling was modeled for three representative reactors (PWR, CANDU, MSR), burning the fuel to 47 GWd/MTHM for PWR, 19 GWd/MTHM for CANDU, and at a steady power of 52.711 MW/MTHM for MSR. Within each model, the protactinium element in the used fuel was extracted and its isotopic content analyzed. Simulated results indicated that 233 Pa can be detected using passive gamma spectroscopy in each fuel type at all decay times (0–300 days) following separation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DRiFT Current Mode, Trigger Settings and Flexible Detector Specifications Applied to Scintillator Arrays

MCNP radiation transport output is post-processed by DRiFT, a Detector Response Function Toolkit to simulate detailed nuclear instrumentation response. DRiFT can be used to assess the performance and potential limitations of scintillator, gas, and semiconductor detectors under a variety of simulated conditions not easily achievable in a laboratory setting. This work describes new updates in DRiFT for scintillator simulations which focus on the capability to simulate scintillators in current mode, an expansion of trigger options, and the ability to customize individual detector properties in a simulation. These improvements are designed to facilitate the ability to model large arrays of scintillator detectors with higher fidelity than was previously possible and are demonstrated in three examples. The first shows the difference between operating DRiFT in current and pulse mode. In the second example, which is intended to demonstrate deviations in individual detector performance, each detector has properties (PMT gain, optical transport, scintillation yield, etc.) that vary between detectors and are specified in DRiFT. A final example examines how DRiFT could be used to optimize digitizer settings in high rate measurements with split signals using the new common trigger option.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Initial Neutronics Investigation of a Chlorine Salt-Based Breeder Blanket

Tritium breeding blankets within D-T-fueled fusion reactors contain lithium compounds and typically require neutron multiplier materials to achieve a tritium breeding ratio (TBR) consistent with self-sustaining operation. Liquid breeder blankets have some advantages over solid blankets, and previous blanket studies have investigated liquid metal as well as liquid salt–based blankets. Liquid salts have reduced magnetohydrodynamic effects as compared to liquid metals, but typically have a lower TBR. Recently, advanced fission reactor concepts have considered chloride-based salts in their design, and there is a significant amount of research work occurring to study these salts. Chloride salts have previously been considered for fusion reactors, but studies have typically found lower breeding ratios than for fluoride salts, such as 2(LiF)-BeF 2 (flibe) so they have not been further developed. In this work, we use a one-dimensional cylindrical radiation transport model of a conceptual tokamak fusion reactor to investigate the neutronics feasibility of using a chloride salt–based blanket that uses chlorine enriched in 37 Cl, which has both a low capture cross section and a substantial (n,2n) cross section. It is found that chloride salts (LiCl mixed with BeCl 2 and/or PbCl 2 ) can potentially achieve a ~3% to 5% higher TBR than fluoride molten salts, notably flibe, in the absence of a solid multiplier. Including a solid multiplier, however, does narrow this advantage, with TBRs estimated within ~1% of flibe with a 2-cm Be multiplier. Chloride salts can also reach lower melting points than flibe, potentially improving the scope for the use of reduced activation ferritic-martensitic steel as a structural material. There is substantial uncertainty in the calculations driven by limited thermochemical data for the Cl salts, plus cross-section uncertainties. The production of 36 Cl through 35 Cl(n,g) and 37 Cl(n,2n) has the potential to challenge the waste disposal rating of the blanket. Calculations indicate that, while this is not an immediate showstopper, this case depends upon the exact waste disposal rating criteria used for 36 Cl. Further work could reduce these uncertainties with improved thermochemical data, higher-fidelity modeling for downselected salts, and more refined waste disposal calculations and regulatory guidance. Lastly, it must be recognized that, as for all molten salts, corrosion and chemistry can present appreciable technical challenges that require further assessment in developing a practical blanket concept, and also that the enrichment of chlorine presents an additional technical and supply chain challenge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiation Damage Analysis of FNSF Components Using McCad and MCNP

The Fusion Energy System Studies Fusion Nuclear Science Facility (FESS-FNSF) concept represents a transitional step between ITER and a commercial fusion power plant. The FNSF is a conceptualized D-T fueled tokamak with 518 MW of fusion power that has been extensively used to explore and optimize design features. The energetic 14.1-MeV neutrons can produce significant localized heating and activations, and can cause damage to plasma-facing components, which can determine maintenance/outage scheduling needs and also impact the lifetime of the device as a whole. This study illustrates a neutronics analysis that was conducted on a 22.5-degree symmetric sector of the FNSF with the goal of understanding the neutron heating and radiation damage that can be characterized by quantifying the displacements per atom (dpa). Concurrently, this study also focused on the development of analysis capabilities by converting a three-dimensional computer-aided design model of the FNSF into MCNP6.2 input using the McCad code. Accordingly, some confirmatory results on tritium production and the tritium breeding ratio (TBR) are provided to support model validation. The results produced by MCNP6.2 simulations showed that the highest heating and damage occurred in the outboard region, which concentrated approximately 290 MW of the total nuclear heating, in contrast to 97 MW within the inboard region. These results are consistent with previous studies that employed earlier versions of the FNSF concept and different modeling approaches. Furthermore, this study also provides additional details on neutron wall loading, as well as total heating from neutrons and gammas, results which show the total heating of the device (16 sectors) is approximately 477.83 ± 0.80% MW, indicating a neutron energy multiplication factor of 1.15. Additionally, the capability to calculate hydrogen and helium production, as well as dpa, is illustrated. Finally, the neutronics effects of using alternative materials to tungsten carbide were evaluated for the vacuum vessel, low-temperature shield, and structural ring components, which showed that compounds like YH 2 , Mg(BH 4 ) 2 , and ZrH 2 could reduce the total heating on the magnet and also reduce the TBR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle Swarm Optimization Algorithm for Critical Experiment Design

Nuclear criticality experiments are used to validate nuclear cross section data used by simulation software. This is typically achieved by designing a critical system with a high sensitivity to a certain material’s cross section. Once the experiment has been carried out, a high fidelity model of the system is developed into a benchmark. When this benchmark model is simulated by a transport code, some of the difference between the experimental and computational effective neutron multiplication factor can be attributed to inaccurate nuclear data. Nuclear data evaluators then can make adjustments accordingly to improve cross section data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗

Impact of fusion reactor neutronics modeling for transmutation and thermal feedback *

Abstract Fusion neutronics calculations provide important metrics pertinent to fusion device operations, such as tritium breeding ratios (TBRs) and data on heat deposition, material activation, and damage. Because of the high computational burden required to generate a high-fidelity Monte Carlo simulation of a 3D fusion device, various assumptions are made to reduce computational time by simplifying the reactor model or the calculation iteration. This paper explores the impact of fusion neutronics metrics such as the TBR and decay heat of structural materials based on assumptions of material composition in the fusion reactor and temperature modeling of materials. Results show that for compact tokamaks with high power and long operational cycles, the transmutation of structural materials is significant enough to cause a substantial change in the flux spectrum and decrease the TBR by 1.68% after 2 years of full power operation. Additionally, assuming a constant temperature and material density can impact the TBR calculations up to 3%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗