Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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

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

At least 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Studsvik R2 materials test reactor Ad Hoc depletion strategy for the derivation of the fuel isotopic composition of the MPCMIV benchmark

The Ad Hoc Depletion Strategy elaborated by the NINE company, developed in support of the organization of the MPCMIV (Multi-physics Pellet Cladding Mechanical Interaction Validation) benchmark input and output specifications, is presented. This work aims at illustrating the strategy itself and then showing the results obtained with its application over the Studsvik R2 Testing Reactor, which is analyzed in the benchmark. The objective of the application of the strategy is to compute the fuel elements isotopic compositions at the beginning of some core loadings of interest for the benchmark. To this objective, it is necessary to implement first the simulation model of the three single assembly types and perform the infinite lattice depletions, then, to build the full core model and to perform the simulation of the core cycle. All the models and simulations were carried out with the use of the Monte Carlo particle transport code Serpent 2. Finally, the simulations results are assessed against Studsvik isotopic compositions of the fuel elements discharged from the R2 Reactor at the end of the core loading. Several assumptions were necessary during all the steps of the strategy, to overcome the lack of information regarding the core management. For this reason, the solution found at the end of the current analysis may not be completely optimized and further improvements regarding the model assumptions will be tested in a future work. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Brazilian CBP - Technoeconomic analysis data

This data is related to the paper entitled "Techno-economic analysis of sugarcane bagasse and straw conversion into cellulosic ethanol via consolidated bioprocessing". That features the evaluation of sugarcane bagasse and straw conversion to ethanol at stand-alone facilities generating electricity from residues. The following scenarios were evaluated: Conventional, featuring hydrothermal pretreatment, fungal cellulase, and yeast fermentation (current commercial standard); Mid-term consolidated bioprocessing (CBP), relying on bagasse solubilization without pretreatment or cotreatment; and Mature CBP, incorporating cotreatment but no pretreatment and considering significant technological advance of the CBP. Available here are the spreadsheets used for Material and Energy balance calculation, Capital and Operational costs estimation and Cash flow analysis. Also available are the description and python code used for Monte Carlo analysis of the ethanol and capital investment variations. This data can be used as a source to implement other techno-economic analysis in the biorefinary context.

09 BIOMASS FUELS↗
DOE OSTI2011

Graphical User Interface for Simplified Neutron Transport Calculations

A number of codes perform simple photon physics calculations. The nuclear industry is lacking in similar tools to perform simplified neutron physics shielding calculations. With the increased importance of performing neutron calculations for homeland security applications and defense nuclear nonproliferation tasks, having an efficient method for performing simple neutron transport calculations becomes increasingly important. Codes such as Monte Carlo N-particle (MCNP) can perform the transport calculations; however, the technical details in setting up, running, and interpreting the required simulations are quite complex and typically go beyond the abilities of most users who need a simple answer to a neutron transport calculation. The work documented in this report resulted in the development of the NucWiz program, which can create an MCNP input file for a set of simple geometries, source, and detector configurations. The user selects source, shield, and tally configurations from a set of pre-defined lists, and the software creates a complete MCNP input file that can be optionally run and the results viewed inside NucWiz.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Description of Axial Detail for ROK Fuel

For the purpose of NDA simulations of the ROK fuel assemblies, we have developed an axial burnup distribution to represent the pins themselves based on gamma scans of rods in the G23 assembly. For the purpose of modeling the G23 assembly (both at ORNL and LANL), the pin-by-pin burnup map as simulated by ROK is being assumed to represent the radial burnup distribution. However, both DA and NDA results indicate that this simulated estimate is not 100% correct. In particular, the burnup obtained from the axial gamma scan of 7 pins does not represent exactly the same 'average' pin burnup as the ROK simulation. Correction for this discrepancy is a goal of the well-characterized assembly task but will take time. For now, I have come up with a correlation for 26 axial points of the burnup as obtained by gamma scans of 7 different rods (C13, G01, G02, J11, K10, L02, and M04, neglecting K02 at this time) to the average burnup given by the simulation for each of the rods individually. The resulting fraction in each axial zone is then averaged for the 7 different rods so that it can represent every fuel pin in the assembly. The burnup in each of the 26 axial zones of rods in all ROK assemblies will then be directly adjusted using this fraction, which is given in Table 1. Note that the gamma scan data given by ROK for assembly G23 included a length of ~3686 mm, so the first 12 mm and the last 14 mm were ignored to give an actual rod length of ~366 cm. To represent assembly F02 in which no pin-by-pin burnup distribution is given by ROK, we must model it using infinitely-reflected geometry but can look at the effects of measuring in different axial zones by using intermediate burnup files (i.e. smaller burnups than 28 GWd/MTU) and determining which axial zone(s) each burnup represents. Details for assembly F02 are then given in Tables 2 and 3, which is given in Table 1 and has 44 total axial zones to represent the top meter in explicit detail in addition to the other 26 zones. Note that the MCNP files for F02 were created using the Monte Carlo burnup linkage code Monteburns, which saves MCNP input files with detailed compositions as a function of burnup. The 'intermediate burnup files' produced for F02 include a cooling time of 27 years. The axial location of 5 spacers was also included in the ROK F02 assembly in which each spacer contained a length of 3.81 cm. Note that due to the nature of Monteburns, which was run in a special fashion for this problem, the step number increments after the 27 year decay, so the second column of Table 2 refers to the step number that should be used in the Monteburns files.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

DHS Summary Report -- Robert Weldon

This summer I worked on benchmarking the Lawrence Livermore National Laboratory fission multiplicity capability used in the Monte Carlo particle transport code MCNPX. This work involved running simulations and then comparing the simulation results with experimental experiments. Outlined in this paper is a brief description of the work completed this summer, skills and knowledge gained, and how the internship has impacted my planning for the future. Neutron multiplicity counting is a neutron detection technique that leverages the multiplicity emissions of neutrons from fission to identify various actinides in a lump of material. The identification of individual actinides in lumps of material crossing our boarders, especially U-235 and Pu-239, is a key component for maintaining the safety of the country from nuclear threats. Several multiplicity emission options from spontaneous and induced fission already existed in MCNPX 2.4.0. These options can be accessed through use of the 6th entry on the PHYS:N card. Lawrence Livermore National Laboratory (LLNL) developed a physics model for the simulation of neutron and gamma ray emission from fission and photofission that was included in MCNPX 2.7.B as an undocumented feature and then was documented in MCNPX 2.7.C. The LLNL multiplicity capability provided a different means for MCNPX to simulate neutron and gamma-ray distributions for neutron induced, spontaneous and photonuclear fission reactions. The original testing on the model for implementation into MCNPX was conducted by Gregg McKinney and John Hendricks. The model is an encapsulation of measured data of neutron multiplicity distributions from Gwin, Spencer, and Ingle, along with the data from Zucker and Holden. One of the founding principles of MCNPX was that it would have several redundant capabilities, providing the means of testing and including various physics packages. Though several multiplicity sampling methodologies already existed within MCNPX, the LLNL fission multiplicity was included to provide a separate capability for computing multiplicity as well as including several new features not already included in MCNPX. These new features include: (1) prompt gamma emission/multiplicity from neutron-induced fission; (2) neutron multiplicity and gamma emission/multiplicity from photofission; and (3) an option to enforce energy correlation for gamma neutron multiplicity emission. These new capabilities allow correlated signal detection for identifying presence of special nuclear material (SNM). Therefore, these new capabilities help meet the missions of the Domestic Nuclear Detection Office (DNDO), which is tasked with developing nuclear detection strategies for identifying potential radiological and nuclear threats, by providing new simulation capability for detection strategies that leverage the new available physics in the LLNL multiplicity capability. Two types of tests were accomplished this summer to test the default LLNL neutron multiplicity capability: neutron-induced fission tests and spontaneous fission tests. Both cases set the 6th entry on the PHYS:N card to 5 (i.e. use LLNL multiplicity). The neutron-induced fission tests utilized a simple 0.001 cm radius sphere where 0.0253 eV neutrons were released at the sphere center. Neutrons were forced to immediately collide in the sphere and release all progeny from the sphere, without further collision, using the LCA card, LCA 7j -2 (therefore density and size of the sphere were irrelevant). Enough particles were run to ensure that the average error of any specific multiplicity did not exceed 0.36%. Neutron-induced fission multiplicities were computed for U-233, U-235, Pu-239, and Pu-241. The spontaneous fission tests also used the same spherical geometry, except: (1) the LCA card was removed; (2) the density of the sphere was set to 0.001 g/cm3; and (3) instead of emitting a thermal neutron, the PAR keyword was set to PAR=SF. The purpose of the small density was to ensure that the spontaneous fission neutrons would not further interact and induce fissions (i.e. the mean free path greatly exceeded the size of the sphere). Enough particles were run to ensure that the average error of any specific spontaneous multiplicity did not exceed 0.23%. Spontaneous fission multiplicities were computed for U-238, Pu-238, Pu-240, Pu-242, Cm-242, and Cm-244. All of the computed results were compared against experimental results compiled by Holden at Brookhaven National Laboratory.

42 ENGINEERING↗
DOE OSTI2014

Monte Carlo Simulations of Neutron Resonance Transmission Analysis with the Dense Plasma Focus

Simulations were performed using the Monte Carlo radiation transport code MCNP in order to analyze the feasibility of using a dense plasma focus to perform neutron resonance transmission analysis. It was found that the dense plasma focus is well-suited for neutron resonance transmission analysis because of the relatively short pulse duration and large neutron yield. This makes lowenergy neutron spectroscopy using time of flight much more accurate, allowing for more detailed assays of smaller quantities of special nuclear material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Neutronics and Fuel Performance Evaluation of Accident Tolerant Fuel under Normal Operation Conditions

This report details the analysis of neutronics and fuel performance analysis for enhanced accident tolerance fuel, with Monte Carlo reactor physics code Serpent and INL’s fuel performance code BISON, respectively. The purpose is to evaluate two of the most promising candidate materials, FeCrAl and Silicon Carbide (SiC), as the fuel cladding under normal operating conditions. Substantial neutron penalty is identified when FeCrAl is used as monolithic cladding for current oxide fuel. From the reactor physics standpoint, application of the FeCrAl alloy as coating layer on surface of zircaloy cladding is possible without increasing fuel enrichment. Meanwhile, SiC brings extra reactivity and the neutron penalty is of no concern. Application of either FeCrAl or SiC could be favorable from the fuel performance standpoint. Detailed comparison between monolithic cladding and hybrid cladding (cladding + coating) is discussed. Hybrid cladding is more practical based on the economics evaluation during the transition from current UO 2 /zircaloy to Accident Tolerant Fuel (ATF) system. However, a few issues remain to be resolved, such as the creep behavior of FeCrAl, coating spallation, inter diffusion with zirconium, etc. For SiC, its high thermal conductivity, excellent creep resistance, low thermal neutron absorption cross section, irradiation stability (minimal swelling) make it an excellent candidate materials for future nuclear fuel/cladding system.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Surface Triangulation for CSG in Mercury

Visualization routines for rendering complicated geometries are very useful for engineers and scientists who are trying to build 3D prototypes of their designs. A common way to rapidly add interesting features to a 3D model is through the use of a concept called Constructive Solid Geometry. CSG uses compositions of the boolean set operations to manipulate basic geometric primitives to form more complicated objects. The most common boolean operations employed are union, intersection, and subtraction. Most computer-aided design software packages contain some sort of ability visualize CSG. The typical workflow for the user is as follows: The user specifies the individual primitive components, the user arbitrarily combines each of these primitives with boolean operations, the software generates a CSG tree structure which normally stores these solids implicitly with their defining equation, the tree is traversed and a general algorithm is applied to render the appropriate geometry onto the screen. Algorithms for visualizing CSG have been extensively developed for over a decade. Points sampled from the implicit solids are typically used as input by variations of algorithms like marching cubes and point-cloud surface reconstruction. Here, we explain a surface triangulation method from the graphics community that is being used for surface visualization in the framework of a Monte-Carlo neutron transport code called Mercury.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Incorporation of Spatial Stochastic Variability into Paracousti-UQ

Waves propagating through natural materials such as ocean water encounter spatial variations in material properties that cannot easily be predicted or known in advance. Deterministic wave simulation algorithms must assume that all properties throughout the model space are precisely known. However, a stochastic wave simulation tool can parameterize the material as a stochastic medium with a certain probability distribution and correlation length. This report documents the addition of spatial stochastic variability into Paracousti-UQ, Sandia Geophysics Department's 3-D full waveform acoustic algorithm within stochastic media. The ability of the code to replicate Monte Carlo solutions in 1-D spatially variable media is also evaluated.

58 GEOSCIENCES↗
DOE OSTI2018

Implementation and Validation of Photon Transport in OpenMC

OpenMC is a Monte Carlo particle transport code focused on reactor physics calculations. It stochastically simulates neutrons moving through 3D models using constructive solid geometry and continuous-energy cross sections. Recently, a photon transport capability was added to OpenMC. In reactor simulations, photon transport is needed, for example, to account for photon heating: although the majority of the energy deposited at a fission site comes from the kinetic energy of the fission fragments, a nontrivial portion can be carried away by photons. Therefore, photons produced in neutron reactions must be simulated in order to accurately model energy deposition. OpenMC can now simulate neutron-induced photon production as well as perform pure photon calculations. The four basic photon interactions with matter — coherent (Rayleigh) scattering, incoherent (Compton) scattering, the photoelectric effect, and pair production — are modeled. Secondary processes that can create new photons (i.e., atomic relaxation, electron-positron annihilation, and bremsstrahlung) are also included. The methods implemented here were validated against MCNP6 for a variety of materials and energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Activation and residual dose rates in the NSLS-II storage ring

It is well known that components in high energy accelerators may get activated and expose workers during maintenance activities. While the potential for activation is higher in positive ion accelerators due to the prolific production of neutrons, certain regions and components in high energy electron accelerators too can be activated enough to create residual dose rates that require careful attention for ALARA compliance. In this report, the activation and residual dose rates from a few materials are studied from various beam loss scenarios ranging from chronic point losses to intermittent real loss scenarios obtained from the operational parameters of the NSLS-II storage ring. FLUKA Monte Carlo radiation transport code is used to model the geometry, the time profile of the beam losses and score the resulting time evolution of the remnant activity and the residual dose rates.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling of Pyroprocessing Hot Cell for Process Monitoring Evaluation

In this report, we examine studies undertaken to determine the applicability of applying measurements with a neutron counter within a pyroprocessing hot cell. The studies focus on the modeling of the hot cell and neutron counter in Monte Carlo N- Particle Code (MCNP) to determine expected radiation effects and applicability of neutron counters. The neutron counter of focus was the High Dose Neutron Detector (HDND) developed at Los Alamos National Laboratory (LANL). We studied the HDND to determine what would be the detector response under normal background conditions with operating process equipment as well as the background’s potential effect on counting applications involving a U/TRU ingot produced by the pyroprocessing electrorefiner. To perform this task, we were required to develop models of the pyroprocessing process cell and unit operations and generate source terms for various different times in the process. The data for the pyroprocessing materials came from SNL’s Safeguards and Security Process Model (SSPM) as part of MPACT’s advanced integration mission. This report will discuss the development of the pyroprocessing process cell and unit operations models in MCNP as well as various simulations to determine the applicability of the HDND for process monitoring efforts in the hot cell.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Enhanced Visualization Capabilities in OpenMC: Report of Activities under a DOE GAIN Voucher

This work was performed under a Cooperative Research and Development Agreement (CRADA) through the U.S. Department of Energy’s Gateway for Accelerated Innovations in Nuclear (GAIN) program in partnership with Oklo Inc., the industry participant. The broad goal of this project was to improve the visualization capabilities of OpenMC related to model development and simulation results. OpenMC is a Monte Carlo particle transport code that is developed and maintained by a broad community of contributors, of which Argonne National Laboratory plays a leading role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Compendium of Material Composition Data for Radiation Transport Modeling

In 2011, Pacific Northwest National Laboratory (PNNL) produced a document known as the Materials Compendium, or Compendium of Material Composition Data for Radiation Transport Modeling, PNNL 15870, Rev. 1, that contains material information useable for modeling purposes for properties of 372 materials. This information is used in several modeling programs used by the radiological/nuclear community, though it is primarily tailored for the Monte-Carlo-N-Particle code produced by Los Alamos National Laboratory. This new document Revision 2 includes a complete review and update of all materials data and references, addressing discrepancies and changes in materials data or references that have occurred since the first revision, an additional 40 materials have been added, primarily newer detector materials developed since the last revision, and isotopic specificity.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Cost-optimized Automated Variance Reduction for Highly Angle-dependent Radiation Transport Analyses [Dissertation]

Monte Carlo variance-reduction techniques that directly bias particle direction have historically received limited attention despite being useful for many radiation transport applications such as those in which radiation travels through large regions with a low probability of colliding. One such technique is known as DXTRAN (short for deterministic transport) in the MCNP Monte Carlo radiation transport code. Until now, effectively applying DXTRAN in calculations is largely based on empirical observations of computational performance. Optimal DXTRAN parameters are identified through manual iteration. This work develops new mathematical descriptions of the DXTRAN variance-reduction process and demonstrates a new automated variance-reduction method that applies these mathematical formulations to determine the optimal application of DXTRAN in a given problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Modeling of the Advanced Test Reactor Using OpenMC, Cubit and Griffin

In the pursuit of the ability to perform multiphysics simulations of the Advanced Test Reactor, it is crucial to have a fast and highly accurate deterministic model. To achieve this, a contemporary two-step method is employed. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor multiphysics application based on the Multiphysics Object-Oriented Simulation Environment. To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material IDs are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Preliminary comparisons indicate good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 50 pcm in the two-dimensional geometry configuration. However, in three-dimensional calculations, an unacceptably large error is found in the Griffin solution. Future work is planned to resolve this discrepancy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Improved Fast Reactor Capability of Griffin in FY23

Griffin is a MOOSE based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE-NE NEAMS program. In FY23, we enhanced capabilities required for fast reactor analysis. This effort included primarily updating the cross-section generation workflow using MC2-3 for various reactor configurations, such as homogeneous, duct-heterogeneous, ring-heterogeneous, and fully-heterogeneous geometries. In addition, we initiated the implementation of a multi-cycle depletion and shuffling capability. To support fast reactor simulation capabilities, we significantly improved the performance of the DFEM-SN-based R-Z transport solver to efficiently solve ultrafine group (over 1000 groups) transport problems. Additionally, the performance of HFEM-PN was improved by introducing red-black iteration, the cmfd acceleration technique, and various optimizations. We also completed the pin power reconstruction capability to support multiphysics simulations while identifying and addressing issues associated with SPH equivalence parameter approach. These enhanced capabilities for fast reactor core simulations, specially HFEM-PN and pin power reconstruction features, were applied to benchmark problems involving ABTR and ABR-1000. These applications showcased excellent agreement with Monte Carlo and other code solutions in terms of eigenvalue, control rod worth, and assembly and pin powers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Radiation Dose Modeling for Niowave’s Accelerator Driven Uranium Target Assembly 3

Molybdenum-99 is a high-value radionuclide commonly used for medical purposes within the United States. The National Nuclear Security Administration (NNSA) seeks to reliably produce the radioisotope 99 Mo without the use of highly enriched uranium. NNSA’s Office of Material Management and Minimization (M3) provides funding and government laboratory expertise to private companies to expedite the production process domestically and currently funds designs that use low-enriched uranium or other 99 Mo production pathways. Several production designs are being explored across the industry, including uranium fission and photonuclear conversion of 100 Mo targets. Niowave Inc. seeks to produce 99 Mo via a high-energy electron accelerator that strikes a lead-bismuth eutectic target that ultimately produces a consistent neutron flux. The neutron flux then interacts in a subcritical reactor core configuration to produce fission in low-enriched or natural uranium targets. These fissionable targets are then processed to extract 99 Mo. The purpose of this work is to estimate the neutron and photon dose response across Niowave’s proposed facility for worker safety during operation. Owing to the size of the proposed Niowave facility and necessary shielding, unbiased Monte Carlo radiation transport is impractical, and variance reduction methods are required. This work focuses on the weight window variance reduction method to produce high confidence dose response results within a Monte Carlo radiation transport code. Specifically, an adjoint-informed weight window methodology was created to improve the dose response estimates for accelerator-driven subcritical reactor designs. This adjoint-informed methodology was implemented for Niowave’s proposed design and improved dose results at far-field locations across the facility. Acceptable dose rate contours for the proposed facility were generated across the facility and are presented in this work.

07 ISOTOPE AND RADIATION SOURCES↗