Search NASA⌕ Search

SEARCH · Search NASA

Results for “Residual Monte Carlo”

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 OSTI2022

Transport error estimation using residual Monte Carlo

The residual Monte Carlo (RMC) method is also known in the literature as sequential Monte Carlo and reduced-source Monte Carlo. Given a Monte Carlo method for solving a linear equation and an approximate solution to that system, the residual method enables use of essentially the same Monte Carlo algorithm to directly compute the additive error or “defect” associated with the approximate solution. As the size of the defect decreases relative to the size of the solution, the residual Monte Carlo method becomes increasingly efficient relative to the standard Monte Carlo (SMC) method. Here we present a new RMC algorithm for evaluating the space-angle error in S n radiation transport solutions, and provide computational examples demonstrating that it can be far more efficient than SMC for this purpose. Herein we also describe a particular pitfall that must be avoided if RMC is to be efficient, and explain why the performance of RMC can significantly differ between different transport problems and different quantities of interest for the same problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Compositional inhomogeneities in a single Icelandic tholeiite flow

New trace element analyses of 25 basalt samples from a vertical traverse in an 11 m thick flow of Icelandic tholeiite show compositional differences which exceed analytical uncertainties and are random with respect to height. These variations suggest a short-range segregation model in which the amounts of phenocrysts, groundmass materials, and residual liquid vary randomly in different samples of the flow; a least-squares method is applied to determine if these variations reflect different proportions of crystallizing phases and residual liquid. Monte Carlo computations indicate that short-range segregation results in inhomogeneity of a magnitude of interflow compositional differences; this causes incoherence in elements with different geochemical behaviors while preserving coherence in elements of similar behavior.

Lindstrom, M. M.↗
NASA NTRS2012

Monte-Carlo Simulations of the Suzaku-XRS Residual Background Spectrum

Cryogenic micro-calorimeters are suitable to detect small amounts of energy deposited by electromagnetic and nuclear interactions, which makes them attractive in a variety of applications on ground and in space. The only X-ray microcalorimeter that operated in orbit to date is the X-Ray Spectrometer on-board of the Japanese Suzaku satellite. We discuss the analysis of the components of its residual background spectrum with the support of Monte-Carlo simulations.

X-ray↗
DOE OSTI2024

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Monte Carlo Simulation with CAD Interface for Calculation of 3D Maps of Residual Dose (CRADA)

Objective: To develop an easy-to-use software application to predict and mitigate radiation effects in research environment, space instruments, nuclear plants and medical facilities and help nonproliferation and national security efforts. Tech-X will develop standalone software libraries and command-line tools for ( 1) translating CAD into tessellated surfaces and tetrahedral meshes in GDML (for Geant4 and MARS 15), ROOT (for MARS 15) and HDF5 (for compact representation and for the visualization) formats, (2) healing CAD geometries to make them suitable for Monte Carlo simulations; (3) creating uniform and variable Cartesian and cylindrical meshes for detailed scoring; and ( 4) efficient Monte Carlo navigation in CAD geometries. JLAB will finish automation of simulations of residual dose in CAD geometries and integrate Tech-X software into Geant4 and MARS15. Finally, Tech-X will develop a Graphical User Interface to set up and heal CAD geometries, create input files, run and visualize simulations for residual dose. This application will run on local desktops, local and remote clusters and supercomputers and will be made available through public clouds, such as Amazon Web Services.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Upsampling Monte Carlo reactor simulation tallies in depleted LWR assemblies fueled with LEU and HALEU using a convolutional neural network

Simulating nuclear reactor cores at the highest achievable spatial and energy resolution is critical in modeling these systems accurately. Increasing the resolution, however, can dramatically increase the memory and central processing unit time required to run simulations. A convolutional neural network was shown previously to accurately upsample tally results of simulated light water reactor assemblies fueled with fresh, low enriched uranium. Here, we show that a convolutional neural network can be used to upsample tally results in assemblies containing fresh and depleted fuel enriched from 1.6 to 19.9 atom percent. The network was trained using neutron flux tallies from simulations of light water reactor assemblies with a range of fuel and coolant temperatures and a diverse selection of geometries. Accurate predictions of flux tallies are possible even on test assemblies with geometries and burnup levels well outside the range of those present in the training and validation data. The network improves the data density by a factor of 8 over a broad range of light water reactor assemblies while incurring insignificant additional computational cost to a Monte Carlo simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2024

Upsampling Monte Carlo Reactor Simulation Tallies in Depleted Sodium-Cooled Fast Reactor Assemblies Using a Convolutional Neural Network

The computational demand of neutron Monte Carlo transport simulations can increase rapidly with the spatial and energy resolution of tallied physical quantities. Convolutional neural networks have been used to increase the resolution of Monte Carlo simulations of light water reactor assemblies while preserving accuracy with negligible additional computational cost. Here, we show that a convolutional neural network can also be used to upsample tally results from Monte Carlo simulations of sodium-cooled fast reactor assemblies, thereby extending the applicability beyond thermal systems. The convolutional neural network model is trained using neutron flux tallies from 300 procedurally generated nuclear reactor assemblies simulated using OpenMC. Validation and test datasets included 16 simulations of procedurally generated assemblies, and a realistic simulation of a European sodium-cooled fast reactor assembly was included in the test dataset. We show the residuals between the high-resolution flux tallies predicted by the neural network and high-resolution Monte Carlo tallies on relative and absolute bases. The network can upsample tallies from simulations of fast reactor assemblies with diverse and heterogeneous materials and geometries by a factor of two in each spatial and energy dimension. The network’s predictions are within the statistical uncertainty of the Monte Carlo tallies in almost all cases. This includes test assemblies for which burnup values and geometric parameters were well outside the ranges of those in assemblies used to train the network.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Calibration and Performance of the ICARUS Detector at Fermilab

ICARUS is the largest Liquid Argon Time Projection Chamber (LArTPC) in operation and serves as the Far Detector of the Short Baseline Neutrino (SBN) program at Fermilab. Precise detector calibration is essential for reliable energy reconstruction and for maximizing the physics reach of the experiment. In this talk, I will describe the energy and timing calibration procedures developed at ICARUS. Thanks to the abundant flux of cosmic ray muons at the surface, the TPC response to ionization charge has been equalized across the detector, removing non-uniformities and in-transparency effects. A novel data-driven procedure has been used to tune the simulation of ionization signals and electronics noise yielding close agreement between Monte Carlo and data with minimal residual bias in reconstructed charge. In addition, the light detection system has been calibrated to sub-nanosecond timing resolution. This precision enables efficient cosmic background rejection and accurate synchronization of neutrino interactions with the time structure of both the Booster Neutrino Beam (BNB) and Neutrinos at the Main Injector (NuMI) beam. Together, these efforts demonstrate the robust performance of the detector and provide critical experience for future large LArTPCs, such as DUNE.

Vicenzi, Matteo [Brookhaven]↗
NASA NTRS2005

Proper Motions of Dwarf Spheroidal Galaxies from Hubble Space Telescope Imaging: Measurement for URSA Minor - 3

This article presents a measurement of the proper motion of the Ursa Minor dwarf spheroidal galaxy determined from images taken with the Hubble Space Telescope in two distinct fields. Each field contains a quasi-stellar object that serves as the "reference point". Integrating the motion of Ursa Minor in a realistic potential for the Milky Way produces orbital elements. The perigalacticon and apogalacticon are 40 (10, 76) and 89 (78, 160) kpc, respectively, where the values in the parentheses represent the 95% confidence intervals derived from Monte Carlo experiments. The eccentricity of the orbit is 0.39 (0.09, 0.79), and the orbital period is 1.5 (1.1, 2.7) Gyr. The orbit is retrograde and inclined by 124 degrees (94 deg, 36 deg ) to the Galactic plane. Ursa Minor is not a likely member of a proposed stream of galaxies on similar orbits around the Milky Way, nor is the plane of its orbit coincident with a recently proposed planar alignment of galaxies around the Milky Way. Comparing the orbits of Ursa Minor and Carina shows no reason for the different star formation histories of these two galaxies. Ursa Minor must contain dark matter to have a high probability of having survived disruption by the Galactic tidal force until the present.

FLUX RESIDUALS↗
DOE OSTI2022

Mirror nucleon removal reactions in p-shell nuclei

Nucleon removal reactions have been shown to be an effective tool for studying the single particle structure of nuclei. This work continues efforts to experimentally probe and benchmark the reaction and structure models used to calculate the removal reaction cross sections when using microscopic nuclear structure inputs. Here, three different single nucleon removal reactions were performed, from p -shell nuclei with masses A = 7, 9, and 10. The residual nuclei from the reactions were detected in coincidence with γ rays to determine partial cross sections to individual final states. The eikonal direct-reaction model is combined with overlap functions and residual nucleus densities from microscopic, variational Monte Carlo calculations to provide consistent nuclear structure input to the partial cross section calculations. Comparisons of measured and calculated cross sections, including for mirror reactions, are presented. The analysis of the partial cross sections leading to the ground states shows a similar behavior to the one observed from analyses of inclusive cross sections using shell model nuclear structure input: the theoretical description of the removal process is in better agreement with the data when removing weakly bound nucleons, than when removing well-bound ones. The two mirror reaction pairs presented here show consistent results between the respective members of the pairs. The results obtained for the population of the excited states, however, show a systematically different trend that appears connected to the structure part of the calculation. Additional cases are needed to better understand the respective roles of structure and dynamical effects in the deviations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2009

Lunar Outpost Life Support Architecture Study Based on a High Mobility Exploration Scenario

As scenarios for lunar surface exploration and habitation continue to evolve within NASA s Constellation program, so must studies of optimal life support system architectures and technologies. This paper presents results of a life support architecture study based on a 2009 NASA scenario known as Scenario 12. Scenario 12 represents a consolidation of ideas from earlier NASA scenarios and includes an outpost near the Lunar South Pole comprised of three larger fixed surface elements and four attached pressurized rovers. The scenario places a high emphasis on surface mobility, with planning assuming that all four crewmembers spend roughly 50% of the time away from the outpost on 3-14 day excursions in two of the pressurized rovers. Some of the larger elements can also be mobilized for longer duration excursions. This emphasis on mobility poses a significant challenge for a regenerative life support system in terms of cost-effective waste collection and resource recovery across multiple elements, including rovers with very constrained infrastructure resources. The current study considers pressurized rovers as part of a distributed outpost life support architecture in both stand-alone and integrated configurations. A range of architectures are examined reflecting different levels of closure and distributed functionality. Different lander propellant scavenging options are also considered involving either initial conversion of residual oxygen and hydrogen propellants to water or initial direct oxygen scavenging. Monte Carlo simulations are used to assess the sensitivity of results to volatile high-impact mission variables, including the quantity of residual lander propellants available for scavenging, the fraction of crew time away from the outpost on excursions, total extravehicular activity hours, and habitat leakage. Architectures are evaluated by estimating surpluses or deficits of water and oxygen per 180-day mission and differences in fixed and 10-year-total equivalent system mass (ESM) relative to a reference case. Results are presented based on current assumptions for Scenario 12 and based on Monte Carlo simulations with assumed probability distributions for the high-impact mission variables. The calculated probability of no water or oxygen resupply from Monte Carlo simulations provides a quantitative measure of system robustness that can be used for cost/benefit analyses to identify leading architecture candidates. Areas of technology improvement that are likely to have a significant impact are also suggested.

Lange, Kevin E.↗
NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗
NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗
DOE OSTI2025

Calibration and simulation of ionization signal and electronics noise in the ICARUS liquid argon time projection chamber

The ICARUS liquid argon time projection chamber (LArTPC) neutrino detector has been taking physics data since 2022 as part of the Short-Baseline Neutrino (SBN) Program. This paper details the equalization of the response to charge in the ICARUS time projection chamber (TPC), as well as data-driven tuning of the simulation of ionization charge signals and electronics noise. The equalization procedure removes non-uniformities in the ICARUS TPC response to charge in space and time. This work leverages the copious number of cosmic ray muons available to ICARUS at the surface. The ionization signal shape simulation applies a novel procedure that tunes the simulation to match what is measured in data. The end result of the equalization procedure and simulation tuning allows for a comparison of charge measurements in ICARUS between Monte Carlo simulation and data, showing good performance with minimal residual bias between the two.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1988

Observations of comets and asteroids

All available ground-based observational techniques are used to study the chemical and physical properties of the small bodies of the solar system, primarily comets and secondarily asteroids. The ultimate goal is to use these bodies to understand the formation and evolution of the solar system. The analysis of charge coupled device (CCD) images of comet P/Halley continued. This included a mojor revision of a paper on the periodicities of CN jets in Halley showing that their periodicity is 7.3 days, not 2.2 days. It also included an analysis of the continuum images which showed that the dust in jets is much redder than in the ambient coma. Calculations with Mie theory suggest that the particles that show the strongest effects of radiation pressure. Much effort was expended devising methods to estimate the level of sky background of the many images in which comet extends to the edge of the chip. Examined images of comet Wilson were taken over many months in 1987. Wilson, a dynamically new comet, shows none of the short-term variability or jet structure (either dusty jets or jets of radicals seen in Halley. This is presumably due to the lack of a mantle or the nucleus. Median-imaging of sky was as a technique for removing residual flat-fielding errors in CCD images. A Monte Carlo model and a convolution model were developed for studying temporal variability of gaseous daughter products.

Ahearn, Michael F.↗