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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Validation of realistic Monte Carlo plasma gamma-ray source on JET discharges

A novel modelling methodology has been developed for the creation of a realistic plasma gamma-ray source for Monte Carlo transport simulations in the tokamak JET. The methodology couples the TRANSP code for plasma transport calculations with the MCNP Monte Carlo particle transport code, thus connecting plasma physics with gamma-ray transport. Here, we present the validation of the developed source methodology by comparing calculated gamma-ray spectra with measurements performed at JET. The validation focuses on gamma-ray spectra measured by the tangential gamma-ray spectrometer during two JET three ion RF scenario discharges, performed in the JET 2019 deuterium experimental campaign. For validation the calculated plasma gamma-ray spectrum was combined with the neutron induced prompt gamma-ray background, originating in the vacuum vessel, and scaled to absolute values calculating the total number of plasma gamma-ray and neutron emitting reactions. The comparison between calculated and measured gamma-ray spectra shows good agreement with the shape of the calculated gamma-ray spectra matching that of measurements for both studied discharges. Moreover, the calculated absolute values of the gamma-ray spectra were of the same order of magnitude at the position of the gamma-ray detector located at the end of a long line-of-sight in a biological shield. Finally, the comparison has validated the developed plasma gamma-ray source methodology for MCNP photon transport calculations at JET. The validation provides a basis for the developed plasma gamma-ray source to be used as a support for the development of future tokamaks such as DEMO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The HTR-Proteus Benchmark: Analysis and Use as a Verification and Validation Case

This presentation outlines the evaluation and application of the HTR-Proteus benchmark as a verification and validation (V&V) case for advanced reactor modeling tools. The work supports the U.S. Department of Energy’s HALEU Availability Program (HAP) and the joint DOE/NRC DNCSH project, which aims to reduce criticality safety uncertainties in commercial-scale HALEU fuel cycle and transportation systems. The HTR-Proteus experiments, conducted at the Paul Scherrer Institute, provide high-fidelity data for TRISO-fueled, graphite-moderated pebble bed reactors with high neutron leakage—conditions relevant to HALEU transport scenarios. This study focuses on Cores 4.2 and 4.3 of the HTR-Proteus benchmark, analyzing key sources of uncertainty including pebble packing, TRISO particle positioning, and core height. Using Project Chrono for realistic pebble geometries and Serpent, SHIFT, and MCNP for neutronics simulations, the study quantifies the impact of these uncertainties on the effective multiplication factor (keff). Results show that a sample size of 110 pebble configurations is sufficient to converge keff, with ±30 pcm uncertainty due to packing randomness. TRISO positioning and core height variations also significantly influence keff, highlighting the importance of detailed modeling in V&V efforts. The benchmark serves as a valuable test case for validating the Griffin reactor physics code and improving confidence in HALEU system simulations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Updates to Relevance Vector Machine: Multiclass Classification, Variable Selection, and Proof-of-Concept Application to Safeguards Fresh Fuel Verification using List-Mode Neutron Collar Data

To expand the capabilities of safeguards authorities to verify the integrity of fresh fuel assemblies, Oak Ridge National Laboratory has retrofit the existing electronics of the JCC-71 uranium neutron coincidence collar, which contains 18 3 He neutron detectors and an external 241 AmLi(α, n) neutron interrogation source arranged to surround a fresh nuclear fuel assembly. The new electronics system allows analysts to record list-mode neutron multiplicity data in addition to the singles and doubles rates that are currently measured. Based on previous proof-of-concept research, analysis of these new data will identify off-normal fuel configurations in an assembly and characterize or localize the specific partial fuel defects. The purpose of this report it to document the analysis algorithm development and then to demonstrate its capability for the safeguards verification of fresh fuel assemblies using list mode neutron collar data. To analyze the complex list-mode data collected with the upgraded uranium neutron collar, multivariate classification algorithms are being developed using a novel classification method, the relevance vector machine. This approach may be applied to multiclass problems to estimate the probability that test data belongs to one of many possible classes of data. In addition, our method identifies the most useful variables/channels for making predictions, which illuminates the basis for the model’s predictions, and this interpretability is largely unique among data analytics methods. Variable selection occurs during model training and parameter tuning and does not need any external hyperparameter tuning routines. Finally, we apply the modified relevance vector machine to a simulated dataset of list-mode neutron collar data generated with the radiation transport code MCNP. The method can correctly identify off-normal fuel configurations, categorize the data according to four fuel defect scenarios, and rank the channels in the data according to prediction utility. For nuclear safeguards applications, it is concluded that this method has the potential to increase the sensitivity and reliability to detect missing fuel rods from a standard 17 x 17 Pressurized Water Reactor (PWR) fresh fuel assembly. Within this analysis, “off-normal” (i.e., missing fuel rods) were correctly classified in 17 simulated test scenarios with one quarter (25%) of the fresh fuel rods missing using a training data set of 58 simulated measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling and Simulation of MC-Small Design

In order to design a lightweight, highly portable neutron multiplicity counter for use in higher risk regions by mobile field teams, a redesign of the MC-15 is being performed, with the goal of producing a detector below 30 lbs for ease of transport. To achieve this, multiple MCNP® simulations were run, evaluating the number and positioning of the He-3 tubes in the detector as well as the amount of high-density polyethylene (HDPE) surrounding the tubes. Based on these simulations, the proposed redesign of the MC-15 (aka MC-Small) possesses 11 tubes and is 13 in wide, 3.21 in thick and 13 in tall. Including the changes to the positioning of the tubes and the removal of HDPE from less sensitive regions of the detectors (i.e. bases of tubes, back of detector), this resulted in an efficiency of 0.56% when measuring a bare Cf-252 source 19.69 in away, only 28% less efficient (relatively) than the MC-15 while weighing 18 lbs, 12 lbs less than the weight of the MC-15.

42 ENGINEERING↗
DOE OSTI2021

Radiation Transport Simulation of the Sulfur Counting System in the Radiation Metrology Laboratory (RML)

At the request of staff members from the Radiation Metrology Laboratory (RML), a series of Monte Carlo radiation transport calculations were performed using two different models of the detector geometry of the RMLs sulfur counting system. The fraction of electrons from each β-decay of 32 P in the sulfur pellet that enter the window of the sulfur counting system was calculated with both MCNP and ITS. In addition, the differential energy distributions of the electrons entering the counting system window were computed. There was significant agreement between the integral and differential quantities calculated by the two transport codes. Summary tables are for the integral efficiency values are found in the body of the report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

LANL Critical Benchmark Comparison Study and Subsequent Revision for Cases Involving LEU and MIX

As part of an international collaboration within the DOE Nuclear Criticality Safety Program (NCSP), LANL is involved in a comparison study to quantify differences in k-effective results from neutron transport simulations of critical benchmark experiments. The DOE NCSP Mission and Vision details the activity in which the French Institut de Radioprotection et de Sûreté Nucléaire (IRSN) leads the study with LANL, ORNL, and LLNL to compare results of various neutron transport codes and nuclear data libraries for ICSBEP benchmarks held in common by the entities. The task statement from the DOE NCSP Five-Year Execution Plan is given as: "The proposal is for IRSN to lead a new intercomparison based on the MORET code with the latest JEFF-3.3 data and ENDF/B-VIII.0 data, when available, using their existing comprehensive selection of 2,714 benchmarks and collate their results together with those from LLNL (COG), LANL (MCNP) and ORNL (SCALE). Due to the large number of benchmarks involved, this effort is envisioned to take three years with an additional year for IRSN to complete a summary report. The benchmark development will be performed independently to minimize modeling errors through discovery and resolution of discrepant results. A summary report will be generated (led by IRSN) to document the results of this study." This report documents results obtained for revisions made to cases involving Intermediate Enriched Uranium (IEU) and a mixture of Pu and Uranium (MIX), with a focus on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of particular benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

RadSim: Particle Transport Wrapper

RadSim is being developed to provide the capability to simulate radiation source emissions, interpolate results from radiation transport tools into a common format to prepare incident flux, and model radiation detector response. This software package includes wrapper API code which is utilized in the development of the second (transport) task of the RadSim project. This software package only contains tools to interface RadSim with radiation transport tools (GEANT4 and MCNP) and does not include any part of the transport source code itself. The users will need prior access to the transport tools (GEANT4 and MCNP), and can use this software package to interface with the transport tool of their choice. This software package takes the geometry file for a radiation transport problem and converts it into an input deck for the transport tools (GEANT4 and MCNP). Additionally, it also contains filters that convert the output of the transport tools to a common flux format.

Cheung, Hoi Sing↗
DOE OSTI2023

Assessing Several Modeling Approaches in Depletion Studies of a TRISO-Fueled Microreactor

The impact of MCNP6 depletion resolution on core lifetime is examined in the context of the Snowflake microreactor with explicit TRISO fuel. The change in core lifetime and isotope mass as a result of different tracked isotopes, timesteps, and spatial regions is discussed. Calculation speed of a prototype MCNP delta tracking module is compared to a reactivity equivalent physical transformation (RPT). Using a single depletion material underpredicts core lifetime by 15%, and the resolution necessary to converge isotope mass greatly depends on the specific isotope, in addition to the size (and location) of the depletion region. The prototype delta tracking module decreases the CPU time of explicit TRISO criticality calculations by 30%, but does not always result in a speedup when used with depletion. Significant depletion speedup is obtained using RPT (50% faster), and all isotope masses agreed within three percent.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

ACRRF High-Bay Dose Calculations using MCNP (Part B)

ACRR radiation outputs through vertical cavities have been documented in two reports. In Part B, the maximum dose from the unshielded central cavity, FREC-II, and NRS is calculated to inform the safety basis. The maximum dose is ≈1430 mrem per 300 MJ, or ≈40,930 rem/hr for full-power operation. A verification study completes the V&V of the modeling. Supplementary studies of variance-reduction techniques, model sensitivities, and aircraft dose above the ACRRF are included.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

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↗
DOE OSTI2023

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Equivalent Point Sources for Simulating Nuclear/Radiological Search Operations

Modeling and simulation can be used to estimate the detection of different types and strengths of radioactive sources. Simulations of shielded gamma-ray sources in complex environments, such as searches in urban areas using vehicle-mounted detectors, can demonstrate what source strengths and shielding combinations should be able to be detected. For faster calculations, gamma-emitting sources that could have different shielding layers can be simulated ahead of time, resulting in equivalent point source emission distributions that can be used in larger simulation models. Larger-scale models will then run faster using equivalent point sources because the shielding layers around the sources have already been accounted for. Accompanying this paper are 260 equivalent point source emission files for 72 sources with different shielding layers, formatted for three different Monte Carlo codes. Users wanting to request versions of the SCALE and/or MCNP software can find details here: https://rsicc.ornl.gov/

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Monte Carlo Simulations of the Water Draining Experiment of Giacint Critical Assembly

The MCNP6 computer program has been successfully extended to simulate reactor dynamics problems with moving parts of the geometries. Different from the dynamic method developed in other Monte Carlo codes, a movement scheme has been developed to account for the geometrical parts motion during the particle random walk. The MCNP6 computer program has been used to simulate two transient experiments of the Giacint critical assembly. The MCNP6 calculated total neutron flux was compared with that from the Serpent simulation. An excellent agreement was obtained between the results of the two Monte Caro computer programs. The MCNP6 calculated total neutron flux was also compared with the two measured transients. The MCNP6 results predicted a faster transient than the experimental data. The MCNP6 transient simulation was improved with an adjusted geometrical model which shifts the fuel rods slightly to match the measured reactivity worth due to the drained water. (c) 2021 Elsevier Ltd. All rights reserved.

MCNP, Moving Geometries, Monte Carlo Reactor Trans↗
DOE OSTI2022

Scoping Study of Effects on the Advanced Test Reactor’s Center Flux Trap Due to Various Loadings in the H and Inner A Test Positions

This study determined the effects experienced by the In-Pile Tube, of the Advanced Test Reactor’s center flux trap, due to various loading configurations of the H and Inner A positions. The effects of concern included: the axial profile of the linear heat generation rate of each fuel rodlet, the axial flux shape, and spectrum, in the fuel, and cladding, of each rodlet, and the total fission heating rate of the entire In-Pile Tube. Additionally, the reactivity worth of each configuration, relative to an ATR configuration with water filling all concerning test positions, was determined. The test positions were filled with commonly used test specimens, backup test specimens, and homogeneous substances. The study was performed using the modeling software MCNP6.2.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

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.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

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↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

AWCC Simulations

Neutron well counters provide a means to measure fissile mass by detecting coincident neutrons, which is a unique signature to fission. For coincidence counting, mass is determined using a calibration curve established with items that have the same material characteristics. Multiplicity analysis solves the point model equations for mass using known detector parameters. In both cases, mass measurements may exhibit a bias due to item characteristics and should be corrected for. Experiments can be done to estimate bias and material effects, but they require having both the well counter and the material well defined. In pursuit of estimating measurement bias and uncertainty in a High Efficiency Neutron Counter (HENC), an MCNP simulation was matched to measurements performed with an Active Well Coincidence Counter (AWCC) to de ne a neutron source term. The neutron source term will be used to predict measurement performance of similar material types in a HENC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗