Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Monte Carlo N-Particle”

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 37 records · Page 2

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN

MCNP® Code Version 6.3.2 Theory & User Manual (Rev. 1)

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.2, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

42 ENGINEERING

ORNLERASE - EPA COUNTS PER MINUTE CALCULATOR CONVERSION FACTORS

This database is used in the EPA's Superfund Counts Per Minute (CPM) Calculator (https://epa-cpm.ornl.gov/index.html). This database contains conversion factors used in the calculator that help users estimate the expected radiation detector reading, in counts per minute (CPM), corresponding to a measured radioactivity level reported in either pCi/cm² or pCi/g. Because surface contamination and volumetric contamination behave differently, the tool uses separate conversion processes, each based on Monte Carlo N-Particle (MCNP)–derived conversion factors. The CPM Calculator supports conversions for a range of common environmental media—including soil, steel, glass, drywall, concrete, and wood—to generate detector-ready CPM values. The primary goal of this tool is to enable more efficient, real-time field measurements, reducing reliance on laboratory analyses and ultimately saving both time and money during Superfund assessments.

Dolislager, Fred [ORNL] (ORCID:0009000325477921)

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.

Earth and Planetary Science Letters

Systematic measurements of the concentrations of cosmogen Ca-41 (half-life = 1.04 x 10(exp 5) yr) in the Apollo 15 long core 15001-15006 were performed by accelerator mass spectroscopy. Earlier measurements of cosmogenic Be-10, C-14, Al-26, Cl-36, and Mn-53 in the same core have provided confirmation and improvement of theoretical models for predicting production profiles of nuclides by cosmic ray induced spallation in the Moon and large meteorites. Unlike these nuclides, Ca-40 in the lunar surface is produced mainly by thermal neutron capture reactions on Ca-40. The maximum production of Ca-41, about 1 dpm/g Ca, was observed at a depth in the Moon of about 150 g/sq cm. For depths below about 300 g/sq cm, Ca-41 production falls off exponentially with an e-folding length of 175 g/sq cm. Neutron production in the Moon was modeled with the Los Alamos High Energy Transport Code System, and yields of nuclei produced by low-energy thermal and epithermal neutrons were calculated with the Monte Carlo N-Particle code. The new theoretical calculations using these codes are in good agreement with our measured Ca-41 concentrations as well as with Co-60 and direct neutron fluence measurements in the Moon.

K. Nishiizumi

Materials for Low-Energy Neutron Radiation Shielding

Various candidate aircraft and spacecraft materials were analyzed and compared in a low-energy neutron environment using the Monte Carlo N-Particle (MCNP) transport code with an energy range up to 20 MeV. Some candidate materials have been tested in particle beams, and others seemed reasonable to analyze in this manner before deciding to test them. The two metal alloys analyzed are actual materials being designed into or used in aircraft and spacecraft today. This analysis shows that hydrogen-bearing materials have the best shielding characteristics over the metal alloys. It also shows that neutrons above 1 MeV are reflected out of the face of the slab better by larger quantities of carbon in the material. If a low-energy absorber is added to the material, fewer neutrons are transmitted through the material. Future analyses should focus on combinations of scatterers and absorbers to optimize these reaction channels and on the higher energy neutron component (above 50 MeV).

Singleterry, Robert C., Jr.

Monte Carlo Analysis of Pion Contribution to Absorbed Dose from Galactic Cosmic Rays

Accurate knowledge of the physics of interaction, particle production and transport is necessary to estimate the radiation damage to equipment used on spacecraft and the biological effects of space radiation. For long duration astronaut missions, both on the International Space Station and the planned manned missions to Moon and Mars, the shielding strategy must include a comprehensive knowledge of the secondary radiation environment. The distribution of absorbed dose and dose equivalent is a function of the type, energy and population of these secondary products. Galactic cosmic rays (GCR) comprised of protons and heavier nuclei have energies from a few MeV per nucleon to the ZeV region, with the spectra reaching flux maxima in the hundreds of MeV range. Therefore, the MeV - GeV region is most important for space radiation. Coincidentally, the pion production energy threshold is about 280 MeV. The question naturally arises as to how important these particles are with respect to space radiation problems. The space radiation transport code, HZETRN (High charge (Z) and Energy TRaNsport), currently used by NASA, performs neutron, proton and heavy ion transport explicitly, but it does not take into account the production and transport of mesons, photons and leptons. In this paper, we present results from the Monte Carlo code MCNPX (Monte Carlo N-Particle eXtended), showing the effect of leptons and mesons when they are produced and transported in a GCR environment.

Aghara, S.K.

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic- metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic-metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry

Modeling Radiation Sources for Experimental Validation

To build on previously completed experimental work, I have been working to implement a model of our lab’s Americium-Beryllium neutron source and Ludlum detection instruments computationally to provide secondary verification of our experimental results. This has been done using the radiation transport code MCNP6 (Monte Carlo N-Particle). I have developed files to calculate the mass ratio of various composite materials that we have exposed to the radiation source, as well as the components of the detection instruments and converted all of these into geometric coordinates that define the experimental setup. Identifying the proper energy distribution and detector efficiencies as well as the correct interpretation of the detector physics have proven to be the most difficult and time-consuming challenges of the work. It is currently too early for results based on the most updated and accurate models. With this experience, the next step is to develop another source model of a proton beam to simulate exposures on other composite targets. This will be developed based on existing experimental datasets, but not used to validate them. Each of these projects will further define the radiation shielding needs for humans in space radiation environments that astronauts will someday be exposed to beyond Low Earth Orbit.

Juliana Simon

ANS Winter 2024 Summary: MCCAFE: The Monte Carlo Constructor for ATR Fuel Elements

The Irradiation Experiment Neutronics Analysis Department at Idaho National Laboratory (INL) has implemented a new analysis workflow for experiments in the Advanced Test Reactor (ATR). One key piece of this workflow is the Monte Carlo Constructor for ATR Fuel Elements, or MCCAFE. For each ATR operating cycle, the Reactor and Nuclear Safety Engineering (RNSE) Department first solves the core in eigenvalue mode and depletes the driver fuel materials. In a separate calculation, neutronics analysts model and deplete the materials of one or more irradiation experiments, usually in a series of fixed-source Monte Carlo N-Particle (MCNP) models of the ATR for neutron transport calculations. It was desirable to use the results of the former calculations to inform the models of the latter. MCCAFE is a Python program developed using American Society of Mechanical Engineers Nuclear Quality Assurance-1 procedures at INL. Its purpose is to take the calculated results from the RNSE depletion solutions and the measured or projected operating parameters from the Nuclear Data Management and Analysis System (NDMAS) to generate fixed-source models of the ATR core at given points in time across one or more cycles.

99 - GENERAL AND MISCELLANEOUS

Summary Report Of The FY25 Reactor Physics Verification And Validation Exercises In The Advanced Reactor Technologies - Gas-cooled Reactor Program

Valdiation and verification of numerical tools is critical for ensuring reasonable predictions for design scoping, licensing, and safety analsyis. In this report, two reactor physics verification and validation exercises are presented. The first of these exercises focuses on burnup analysis with data from the Advanced Gas Reactor (AGR) program. Simulations are performed with Monte Carlo N-Particle (MCNP) and are compared with the experimental measurements for the AGR 1 and 2 experiments that utilize both UCO and UO2 fuel. The second exercises utilizes data from the HTR-Proteus experiments to perform reactor physics validation. Specifications of the experimental facility are provdied, along with a demonstration of initial modeling efforts in Serpent for one of the determistic packing experiments. Both cases are part of the Generation-IV international forum (GIF) Very High-Temperature Reactor (VHTR) Computational Methods, Validation, and Benchmarking (CMVB) program, an international collaborative organization dedicated to the verification and validation of High-Temperature Gas-Cooled Reactor (HTGR) analysis. Participation in the CMVB allows the US Department of Energy (DOE) to leverage these existing validation activities to provide extra value through benchmarking activities with other CMVB members.

and Benchmarking (CMVB) program

Estimation of General Area Dose Rates During Closure Weld Operations on the DOE Standard Canister Pacific Basin Nuclear Conference 2024

The Department of Energy (DOE) Spent Nuclear Fuel (SNF) Road Ready Demonstration Project is an effort to package DOE-managed SNF for road-ready dry storage. A key operation in this demonstration is performing the closure weld on the DOE Standard Canister after it has been loaded with the SNF. Planning the closure weld operations is ongoing and must consider what parts of the operations (if any) can be done manually and what parts must be done remotely. A key parameter to inform the evolving planning of these operations is the dose rate in the vicinity of the canister and the dose hazard posed to the closure weld operators and collocated workers. This study developed a Monte Carlo N-Particle (MCNP) model of the closure weld operation, which was used to estimate general area dose rates as well as dose rates to the welding location, weld operators, and collocated workers. These estimates provide some insight into the potential radiation exposure hazard for workers during the operations and are a factor to be considered in the ongoing development of designs, technology, processes, and administrative controls for packaging DOE-managed SNF for road-ready dry storage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Comparing Control Performance Between Simulation and Experiment using the Microreactor Automated Control System Testbed

In the advanced reactor domain, a flexible and scalable software/hardware infrastructure is crucial for integrating and validating various control technologies. This study used the Microreactor Automated Control System (MACS) hardware platform as a testbed. MACS was originally designed to mirror Idaho National Laboratory (INL)'s Microreactor Applications Research Validation and Evaluation (MARVEL), a 85-kW thermal fission microreactor. It features control drums for simulated reactivity control; lights that function as a surrogate reactor core, with the brightness being proportional to the reactor power; and light sensors that emulate neutron detectors. To transform MACS into a physical twin of MARVEL for evaluating control methods, the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software was employed. This software integrated the hardware with two models of the MARVEL core, based on Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The study aimed to demonstrate the gap between control theory and actual practice—a gap that often necessitates empirical adjustments such as control gain retuning, filters, time discretization, and integrator anti-windup measures. Controllers were developed based on increasingly complex simulations without hardware, starting from the base MARVEL model and then introducing actuator saturation constraints and sensor noise. The final control strategy was then tested using MACS, and a comparative performance analysis was conducted.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN

Radiation Transport Tools for Space Applications: A Review

This slide presentation contains a brief discussion of nuclear transport codes widely used in the space radiation community for shielding and scientific analyses. Seven radiation transport codes that are addressed. The two general methods (i.e., Monte Carlo Method, and the Deterministic Method) are briefly reviewed.

Cosmic Ray Effects on Micro-Electronics (CREME96)