Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Analysis”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fine Temperature Grid Continuous Energy Cross Section Generation for Monte Carlo Analysis of Xe-100 Design

The standard “A Compact ENDF (ACE)” data libraries used by Monte Carlo based reactor physics codes calculations are provided by Los Alamos National Laboratory (LANL) with a temperature interval mostly of 300 K (e.g. 300 K, 600 K, 900 K) for the cross sections and between 100 K and 200 K for the thermal scattering libraries (TSL). However, some codes such as MCNP lack capability to perform on-the-fly temperature interpolation during simulation both for neutron and TSL cross-sections. To evaluate the impact related to Doppler broadening and spectrum shift associated with TSL changes, this paper explores the potential of adopting a temperature grid finer than the ones contained in the standard data libraries. A 50 K temperature grid was employed to quantify the error in neutronics calculations due to temperature grid resolution. This was achieved by comparing the results of this study (50 K temperature interval) against the results obtained with standard data libraries (>100 K temperature interval). While the adopted grid primarily relies on the ENDF/B-VII.1 library, for neutron cross-sections, it utilizes ENDF/B-VIII.0 library for TSL. The analyses confirmed that the accuracy of neutronics calculations is satisfactory when using a 50 K temperature grid. Notably, adopting a 50 K temperature grid, as opposed to standard libraries or coarser temperature grids, could lead to a difference of no more than a few hundred pcm in dk for both fresh fuel and burnt fuel. The most sensitive reaction type to the temperature grid was as expected identified as the capture cross-section of U-238.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Monte Carlo Analysis of Coolant Stream Impurity Gamma Emissions in Gas-Cooled Fast Reactors

It is well established that a rapid increase in the concentration of fission products in the reactor coolant stream can serve as an early indication of fuel failure. We use Monte Carlo simulations to investigate the feasibility of using several gamma detectors as diagnostic equipment to monitor the presence of major fission product isotopes in high-temperature gas-cooled reactor (HTGR) coolant streams for early detection of fuel failure and therefore the prevention of fuel failure conditions. Herein we model the response of high-purity germanium (HPGe), CdZnTe, NaI(Tl), and LaBr 3 (Ce) detectors of typical commercial sizes to the gamma emissions from nuclides expected to be found within the coolant stream of the Versatile Test Reactor (VTR) under development by the U.S. Department of Energy. The results indicate that for the 233- and 250-keV gamma rays from 133 Xe and 135 Xe, respectively, the 3σ detection criterion is met in under 1 min using a single HPGe detector. Changes in other spectral lines associated with Xe nuclides are detected within 1 h regardless of the choice of detector.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Monte Carlo Analysis of the Performance of the ITER Diagnostic Residual Gas Analyzer

The composition of exhausted gas is a key parameter in long-pulse plasma fusion experiments, and its evolution shall be monitored at timescales relevant to plasma dynamics and plasma-wall interactions. A diagnostic residual gas analyzer (DRGA) is a multisensor instrument particularly suited to these studies, and ITER will adopt DRGAs in the equatorial and in the divertor tokamak regions. In this work, we have revisited the design of the ITER divertor DRGA through simple vacuum analytical considerations supported by simulations conducted with Molflow+, a test particle Monte Carlo (TPMC) simulation code commonly used in the particle accelerator community. Here, starting with recommendations on the manufacturing of the vacuum piping of the DRGA, this work is followed by a complete vacuum characterization of the diagnostic vacuum setup (pressure profiles at base pressure and during sampling, orifice diameter, and length optimization), and finally, the in-vessel residence time of the most important gas species is simulated. These studies have allowed us to give insights into some experimental results recently found on the prototype DRGA installed in the Wendelstein W7-X stellarator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Monte Carlo simulation analysis for radiation damage on Glidcop Al-15 caused by 17-20 MeV/u heavy ions

The Facility for Rare Isotope Beams (FRIB) SRF linear heavy-ion accelerator is designed to accelerate all ions up to uranium to a maximum beam power of 400 kW. Several beam intercepting devices (BIDs) are essential to the successful operation of the accelerator, including a low power charge selector (LPCS) made of copper containing 0.15% precipitated aluminum oxide by weight called Glidcop Al-15. As FRIB ramps up the primary beam power beyond the current 20 kW level, the charge selector must withstand higher radiation damage rates, typically measured in displacements per atom (dpa). Significant beam induced radiation damage including significant deformation, like swelling, blistering, and cracking has been observed on a recently removed LPCS. We also observed physical features up to about 5.5 mm wide and appear to be deeper than the projected range of any ions. This paper presents the results of Monte Carlo simulations carried out using the Particle and Heavy Ion Transport code System (PHITS), quantifying the total damage dose and ion concentration, etc. Simulating an accurate irradiation history is essential to determining the scope of post irradiation examination (PIE) work.

Fletcher, Ethan [Michigan State U.]↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Energy Storage Analysis

This study presents a comprehensive techno-economic characterization of energy storage and exible low carbon power generation technologies that can shift energy across days, weeks, or months to balance daily, weekly, and seasonal disparities in supply and demand. Energy storage technologies evaluated here include pumped hydropower storage (PHS), adiabatic and diabatic compressed air energy storage (CAES), vanadium redox flow batteries (VRBs), pumped thermal energy storage (P-TES), and renewably produced hydrogen stored in either geologic formations or underground pipes with re-electrification via combustion turbines in combined cycles, stationary proton exchange membrane (PEM) fuel cells, or the novel use of PEM fuel cells designed for heavy-duty vehicles (HDVs) which are expected to have shorter lives but also lower capital costs. We also evaluate flexible low-carbon power generation systems, including ethanol combustion in gas turbines and natural gas combustion with carbon capture and sequestration (CCS). We estimate current costs with literature data, use learning rates to characterize future costs, and develop capacity factors calibrated to an 85% renewables grid to calculate the levelized cost of energy (LCOE) of each technology. Results illustrate that at the 12-hour storage duration, PHS and CAES have the lowest LCOE with current costs, and VRBs become competitive if future costs are achieved. At the 120-hour storage duration, hydrogen systems with geologic storage and natural gas with CCS achieve the lowest LCOE in both current and future capital cost scenarios. In particular, the new configuration of HDV-PEM fuel cells with hydrogen storage in geologic formations evaluated here could lower the LCOE by 22-27% compared to stationary fuel cell systems typically evaluated and might help enable very high (>80%) renewable energy electric power systems. P-TES and hydrogen stored in underground pipes are the least-cost options at the 120-hour storage duration rating that do not require some form of geologic storage. Sensitivity analysis and Monte Carlo analysis illustrate that the general trends seen in this study are valid for a wide range of capacity factors and future cost scenarios. The study also illustrates that coproducing and selling hydrogen to other markets could reduce the LCOE of hydrogen systems by up to 39%.

carbon capture and sequestration↗
DOE OSTI2024

Monte Carlo Perturbation Analysis of Fuel Temperature Variations in the MCNP Model of the Annular Core Research Reactor

The Annular Core Research Reactor (ACRR) Monte Carlo N-Particle (MCNP) model is used by ACRR reactor operators and experiment designers at Sandia National Laboratories for a variety of computational calculations ranging from reactor kinetics parameter estimates and safety analyses to experimental planning. To understand the dominant source of uncertainty within the MCNP model, perturbations in temperature were applied to individual ACRR MCNP fuel rods. Fuel rod temperatures were randomly sampled from a uniform distribution from operational temperatures to quantify temperature-related uncertainty effects. Stochastic mixing was used to blend the cross sections of the desired temperatures using the MCNP continuous and Thermal Neutron Scattering Treatment [S(α,β)] libraries in ENDF/B-VII.1. Furthermore, this uncertainty analysis produced a 640 row × 640 column correlation and covariance matrix of the neutron energy spectra. Positive covariance was produced around the 1-MeV region and the 0.2-eV region. Correlation was found in the thermal and fast energy regions, but no correlation was observed in the slowing-down energy region because interactions in this region are not dominated by fuel.

ACRR↗
DOE OSTI2021

Monte Carlo transport analysis to assess intensity dependent response of a carbon-doped GaN photoconductor

In this report evaluation of the photoresponse in wurtzite GaN photoconductive switches is presented based on kinetic Monte Carlo simulations. The focus is on electron transport physics and assessment of high frequency operation. The roles of GaN band structure, Pauli exclusion, and treatment of internal fields based on the fast multipole method are all comprehensively included. The implementation was validated through comparisons of velocity-field characteristics for GaN with computational results in the literature. Photocurrent widths of less than ~7 ps for the 1 μm device can be expected, which translates into a 100 GHz upper bound. Photocurrent pulse compression below the laser full width at half maxima at high applied fields are predicted based on the interplay of space-charge effects and the negative differential velocity characteristics of GaN.

36 MATERIALS SCIENCE↗
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 OSTI2026

Quantitative 14 N NMR with Monte Carlo Uncertainty Analysis of Nitrate/Nitrite in Alkaline Nuclear Waste

While monitoring of nitrate and nitrite concentrations is important for managing corrosion in nuclear waste systems, existing analytical methods are hindered by turbidity, spectral interference, and delays from sample handling. Here, we demonstrate quantitative 14 N nuclear magnetic resonance (qNMR) spectroscopy as a direct, matrix-tolerant approach for nitrate and nitrite detection at natural abundance. Monte Carlo resampling was integrated into the workflow to quantify random error, establish precision–time tradeoffs, and separate noise-limited uncertainty from systematic bias arising from shimming, transmitter offset, or excitation pulse conditions. Quantification of nitrate and nitrite were validated in controlled alkaline matrix challenges and in 18-component Hanford-type simulants. These results establish 14 N qNMR as a practical, uncertainty-bounded tool for monitoring redox-active nitrogen species in chemically complex environments and provide a generalizable framework for quantitative analysis of quadrupolar nuclei.

Graham, Trent R. [Pacific Northwest National Labor↗
DOE OSTI2023

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters

Accurately predicting the criticality of an experiment before interacting with the experimental components is very important for criticality safety. Radiation transport software can be utilized to calculate the effective neutron multiplication factor of a nuclear system. Because of the integral nature of the effective neutron multiplication factor, the value calculated contains various sources of nuclear-data induced uncertainty. The sensitivity analysis and data assimilation technique presented in this paper exhibit one possible method of identifying and reducing the effective neutron multiplication factor nuclear-data induced uncertainty. The results presented in this work show that it is possible to use relative sensitivity coefficients of the prompt neutron decay constant and the effective delayed neutron fraction to 239 Pu nuclear data to reduce nuclear-data induced uncertainties in the effective neutron multiplication factor. This work has been utilized by members of the Los Alamos National Laboratory project EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) for optimally designing a new experiment, which will be used to reduce compensating errors in 239 Pu nuclear data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗