Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A MONTE CARLO CODE FOR RELATIVISTIC RADIATION TRANSPORT AROUND KERR BLACK HOLES

We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption mechanisms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code and also present results from a number of numerical tests to assess its accuracy and convergence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Treatments for Neutron Resonance Elastic Scattering Using the Multipole Formalism in Monte Carlo Codes

The elastic scattering cross sections of resonant nuclides can vary drastically within small energy intervals. In reactor physics simulations, accurate treatment of neutron resonance elastic scattering is vital, as the phenomenon can have a strong effect on scattering kinematics which, in turn, affect macroscopic values, such as the fuel temperature reactivity coefficient. To reduce memory requirements and facilitate coupling with heat transfer, a new representation for cross sections was recently adopted in OpenMC based on the multipole formalism. However, this new approach invalidates the previous implementations of correct resonance treatments, including the Doppler Broadening Rejection Correction (DBRC) method. This work thus developed a modified DBRC algorithm, called DBRC-MP, compatible with the new data representation. This new method is also verified against the previous DBRC method. While slightly more computationally costly than the traditional DBRC approach using ACE format cross sections, the use of the multipole representation in treating resonance scattering enables the reduction of memory requirements by a hundred-fold and facilitates the representation of temperature dependent cross sections. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2022

Expanded Validation of Uranium Systems with the KENO Monte Carlo Codes and SCALE 6.2.4

The development of accident-tolerant fuel and advanced reactors fueled with high-assay low-enriched uranium has created an increased focus on systems with uranium enrichments between 5 and 20 wt % 235 U. SCALE criticality safety validation reports generated at Oak Ridge National Laboratory contain very few benchmark experiments in this range. The work presented herein is an initial investigation of expanded validation of the entire range of enrichments for uranium-fueled systems, including consideration of over 300 newly developed KENO models. The vast majority of the new models are between 4.4 wt % and 37.5 wt % enrichment and include a variety of thermal, intermediate, and mixed neutron spectra. This work finds no strong evidence of an enrichment-dependent bias, but there are indications of a spectrum-dependent bias. There are, however, only a limited number of independent experiments included at this time to distinguish intermediate enrichment effects from intermediate spectrum effects. Further work in this area is ongoing for intermediate enrichment systems. In the long run, the addition of these models to the Verified, Archived Library and Inputs and Data is planned, as funding and staff availability allow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Construction of an Object-Oriented Interface to an Enhanced EGS Monte Carlo Code

The version of SLAC’s EGS Code System (known as EGS4) was clearly recognized throughout the world as a “gold standard” for solving physics and engineering problems related to the transport of electrons and photons in any element, compound, or mixture. However, EGS4 was written in a computer language called Mortran3 and the learning curve for this rather obscure language was difficult and steep. One major objective of this project was to convert the EGS4 Mortran3 coding into Fortran77, making it much simpler for first-time users to understand and appreciate. In addition, since the introduction of EGS4 in 1985, a number of enhancements in electron-photon transport, as well as a growing need for additional physics processes, had been developed, and the new code (called EGS5) would incorporate more than 20 of these improvements. Of course, many key benchmark simulations had to be made during its development in order to prove that EGS5 was just as credible as its predecessor, EGS4. All of the above was tasked to SLAC, with outside support coming from colleagues at the University of Michigan (Nuclear Engineering Department) and at the High Energy Accelerator Research Organization (KEK) in Japan.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Expanded Validation of Uranium Systems with the KENO Monte Carlo Codes and SCALE 6.2.4 [Slides]

There is likely a trend in bias for uranium-fueled systems as a function of spectrum, but probably not as a function of enrichment. Some enrichment trends appear, but are potentially the result of outlier experimental results; further investigation is needed where possible. There are good examples in this data of markedly different biases from different systems. Future work includes review and vetting for inclusion in the VALID library and analysis with TSUNAMI.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗