SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2000

SUPREM-DSMC: A New Scalable, Parallel, Reacting, Multidimensional Direct Simulation Monte Carlo Flow Code

An AFRL/NRL team has recently been selected to develop a scalable, parallel, reacting, multidimensional (SUPREM) Direct Simulation Monte Carlo (DSMC) code for the DoD user community under the High Performance Computing Modernization Office (HPCMO) Common High Performance Computing Software Support Initiative (CHSSI). This paper will introduce the JANNAF Exhaust Plume community to this three-year development effort and present the overall goals, schedule, and current status of this new code.

Campbell, David↗

NASA NTRS1994

Application of electron paramagnetic resonance imaging to the characterization of the Ultem(R) exposed to 1 MeV electrons. Correlation of radical density data to tiger code calculations

A major long-term goal of the Materials Division at the NASA Langley Research Center is the characterization of new high-performance materials that have potential applications in the aircraft industry, and in space. The materials used for space applications are often subjected to a harsh and potentially damaging radiation environment. The present study constitutes the application of a novel technique to obtain reliable data for ascertaining the molecular basis for the resilience and durability of materials that have been exposed to simulated space radiations. The radiations of greatest concern are energetic electrons and protons, as well as galactic cosmic rays. Presently, the effects of such radiation on matter are not understood in their entirety. It is clear however, that electron radiation causes ionization and homolytic bond rupture, resulting in the formation of paramagnetic spin centers in the polymer matrices of the structural materials. Since the detection and structure elucidation of paramagnetic species are most readily accomplished using Electron Paramagnetic Resonance (EPR) Spectroscopy, the NASA LaRC EPR system was brought back on-line during the 1991 ASEE term. The subsequent 1992 ASEE term was devoted to the adaptation of the EPR core system to meet the requirements for EPR Imaging (EPRI), which provides detailed information on the spatial distribution of paramagnetic species in bulk media. The present (1994) ASEE term was devoted to the calibration of this EPR Imaging system, as well as to the application of this technology to study the effects of electron irradiation on Ultem(exp R), a high performance polymer which is a candidate for applications in aerospace. The Ultem was exposed to a dose of 2.4 x 10(exp 9) Rads (1-MeV energy/electron) at the LaRC electron accelerator facility. Subsequently, the exposed specimens were stored in liquid nitrogen, until immediately prior to analyses by EPRI. The intensity and dimensions of the EPR Images that were generated for the irradiated specimens showed that the electrons penetrated the material to a depth of approximately 0.125 inch. These data show a very high degree of correlation to the energy deposition profile as predicted by the Tiger Code, a Monte Carlo code that provides guidelines for the transport of electrons in matter. Subsequent efforts will focus on delineating the transport properties of energetic protons in Ultem(R).

Suleman, Naushadalli K.↗

DOE OSTI2023

ITS Version 6.7: The Integrated TIGER Series of Coupled Electron/Photon Monte Carlo Transport Codes User's Manual

ITS is a powerful software package permitting state-of-the-art Monte Carlo solution of linear time-independent coupled electron/photon radiation transport problems, with or without the presence of macroscopic electric and magnetic fields of arbitrary spatial dependence. Our goal has been to simultaneously maximize operational simplicity and physical accuracy. Through a set of preprocessor directives, the user selects one of the many ITS codes. The ease with which the make system is applied combines with an input scheme based on order-independent descriptive keywords that makes maximum use of defaults and internal error checking to provide experimentalists and theorists alike with a method for the routine but rigorous solution of sophisticated radiation transport problems. Physical rigor is provided by employing accurate cross sections, sampling distributions, and physical models for describing the production and transport of the electron/photon cascade from 1.0 GeV down to 1.0 keV. The availability of source code permits the more sophisticated user to tailor the codes to specific applications and to extend the capabilities of the codes to more complex applications. Version 6, the latest version of ITS, contains (1) improvements to the ITS 5.0 codes, and (2) conversion to Fortran 95. The general user friendliness of the software has been enhanced through memory allocation to reduce the need for users to modify and recompile the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗

DOE OSTI2023

Monte Carlo Analyses of the FOEHN Experiment

The FOEHN critical experiment, which has been carried out at Cadarache (France) in the reactor EOLE in the early 70s, was designed to verify RHF design analyses. This study focuses on the impact of the FOEHN experiment uncertainties on the calculation of the energy deposition. The latter has been obtained by the Monte Carlo codes Serpent and MCNP. In the calculation of the energy deposition, the first code relies on KERMA factors whereas the latter code relies on Q-values. The two Monte Carlo codes share the same geometry and material specifications. In addition to the design and experiment parameters uncertainties, nuclear data also impacts the calculation of the energy deposition. In this study, the MCNP simulations use three different nuclear data library sets: ENDF/B versions VII.0, VII.1, and VIII.0. Finally, the ksens card of MCNP, using perturbation theory, has been used to identify the cross sections and isotopes with largest impact on the effective multiplication factor. By showing that the use of different software and computational methodologies and design parameters uncertainties do not significantly impact the core power distribution, this work demonstrate the appropriateness of the tools and methods employed for the conversion of involute-plate reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Reassessing the MCNP Random Number Generator

Random number generators are integral components to Monte Carlo codes. They provide the pseudorandom number sequence used to actually sample the distributions of interest. As a result, they are one of the most important components to the software. The current recommended MCNP random number generator is a 63-bit linear congruential generator (LCG). This generator is quite fast, but it has some drawbacks. First, it only has a period of 2 63 . Due to the necessarily non-optimal usage of random numbers to ensure parallel reproducibility, this amount is too few to guarantee random number sequences are not reused in all configurations the code runs under. As simulation size increases, users will need to be aware of the limitations of the generator and tune configuration variables to best suit their simulations, or they will need to assume that reuse is not negatively affecting their answers. Neither of these are optimal. Second, small LCGs are fairly weak in bit generation quality, and this can have an unknown impact on the quality of the simulation. This paper is an investigation into whether or not more modern random number generators can supersede the current ones. The goal is to find a generator that is similar or superior in speed to the LCGs, has a state space large enough to make strong guarantees about random number reuse, and passes all modern random number test suites. If such a generator is found, it would eliminate the need for the user to even be aware of the limitations of the random number generator and would simplify the use of the code. This paper will be broken into several parts. Sec. 2 will discuss the evolution of the random number generator within the MCNP code. Sec. 3 will go over what a Monte Carlo code needs from a generator to be reproducible and portable and how the current generator behaves in that light. Sec. 4 goes through how each generator was tested. Finally, Sec. 5 will discuss improvements that could be made to the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2004

Infrared Spectroscopy of Star Formation in Galactic and Extragalactic Regions

Last year we submitted and had accepted a paper entitled "The Far-Infrared Emission Line and Continuum Spectrum of the Seyfert Galaxy NGC 1068," by Spinoglio, L., Malkan, M., Smith. HA, Gonzalez-Alfonso, E., and Fischer, J. This analysis was based on the SWAS Monte Carlo code modeling of the OH lines in galaxies observed by ISO. Since that meeting last spring considerable effort has been put into improving the Monte Carlo code. A group of European astronomers, including Prof. Eduardo Gonzalez-Alfonso, had been performing Monte Carlo modeling of other molecules seen in ISO galaxies. We used portions of this grant to bring Prof. Gonzalez-Alfonso to Cambridge for an intensive working visit. A second major paper on the ISO IR spectroscopy of galaxies, "The Far Infrared Spectrum of Arp 220," Gonzalez-Alfonso, E., Smith. H., Fischer, J., and Cernicharo, J., is in press. Spitzer science development was the major component of this past year;s research. This program supported the development of five Early Release Objects for Spitzer observations on which Dr. Smith was Principal Investigator or Co-Investigator, and another five proposals for GO time. The early release program is designed to rapidly present to the public and the scientific community some exciting results from Spitzer in the first months of its operation. The Spitzer instrument and science teams submitted proposals for ERO objects, and a competitive selection process narrowed these down to a small group with exciting science and realistic observational parameters. This grant supported Dr. Smith's participation in the ERO process, including developing science goals, identifying key objects for observation, and developing the detailed AOR (observing formulae) to be use by the instruments for mapping, integrating, etc.). During this year Dr. Smith worked on writing up and publishing these early results. The attached bibliography includes six of Dr. Smith's articles. During this past year Dr. Smith also led or helped to develop proposals for ten Spitzer GO Programs, and three others. Appendix B lists the programs involved.

Smith, Howard A.↗

DOE OSTI2021

Reactor cell neutron dose for the molten salt breeder reactor conceptual design

The private sector’s interest in the active development of molten salt reactors has led to the need to develop and test advanced modeling and simulation tools to analyze various advanced reactor types under numerous conditions. This paper discusses the effort undertaken to model the Oak Ridge National Laboratory (ORNL) Molten Salt Breeder Reactor (MSBR) design using ORNL’s Shift Monte Carlo code. The MSBR model integrates a Monte Carlo N-Particle (MCNP) MSBR core model with an MCNP model that was generated from a CAD model of the external components and the reactor building, which was subsequently run in Shift. This paper focuses on development of the fully integrated model and its use in performing neutron transport calculations in the reactor cell area. This model is intended to aid in understanding radiological dose conditions during operation, as well as the iron dpa rates in the reactor vessel. The neutron biological dose rates and flux calculated in the reactor cell are much higher in the MSBR than in typical light-water reactors. The implications of these results and future work are also discussed in this paper.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group

Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing. In this work, to respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics. G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes. This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2012

Use of Fluka to Create Dose Calculations

Monte Carlo codes provide an effective means of modeling three dimensional radiation transport; however, their use is both time- and resource-intensive. The creation of a lookup table or parameterization from Monte Carlo simulation allows users to perform calculations with Monte Carlo results without replicating lengthy calculations. FLUKA Monte Carlo transport code was used to develop lookup tables and parameterizations for data resulting from the penetration of layers of aluminum, polyethylene, and water with areal densities ranging from 0 to 100 g/cm^2. Heavy charged ion radiation including ions from Z=1 to Z=26 and from 0.1 to 10 GeV/nucleon were simulated. Dose, dose equivalent, and fluence as a function of particle identity, energy, and scattering angle were examined at various depths. Calculations were compared against well-known results and against the results of other deterministic and Monte Carlo codes. Results will be presented.

Lee, Kerry T.↗

DOE OSTI2022

Comparative analysis of energy deposition modes available in Serpent 2 within the framework of the supercritical water reactor - Fuel qualification test reactor physics benchmark

A joint European Canadian Chinese development of a supercritical water-cooled small modular reactor (SCW-SMR) technology is in progress since September 2020 in the framework of a Horizon 2020 project called ECC-SMART. As a main purpose of the project, proper estimates of energy deposition and its spatial distribution are prerequisites for the accurate analysis of safety related parameters of the SCW-SMR concept under development. A supercritical water reactor fuel computational benchmark model, provided by Canadian Nuclear Laboratories, was applied for detailed comparison of different energy deposition calculation options available in the Serpent 2 Monte Carlo code. The effect of energy deposition options on the normalization of the results as well as on the spatial distribution of the energy deposition are discussed. Consistent energy deposition calculation methods are presented between three Monte Carlo codes, viz., Serpent 2, MCNP6 and OpenMC. Although resource-intensive, the use of the coupled neutron-photon transport mode of Serpent 2 is recommended for accurate spatial and quantitative characterization of energy deposition in the SCW-SMR fuel assemblies, accounting for both neutron and photon heating of all the materials. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2018

Universal and Automated Monte Carlo Method Code for Uncertainty Propagation in Metrology Databases

A method of providing transformed target points for integrating a component into an assembly includes collecting a set of component target points, collecting a set of assembly target points, identifying target points common to the set of component target points and the set of assembly target points; performing a specified number of Monte Carlo transformations of selected ones of the common target points to yield a set of transformed target points and vectors and an associated uncertainty value for each transformed target point and vector, and using certain ones of the transformed target points for integrating the component into the assembly based on the associated uncertainty value for each of the transformed target points.

Hayden, Joseph↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗