Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Sampling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Multilevel ensemble Kalman filtering

This study embeds a multilevel Monte Carlo sampling strategy into the Monte Carlo step of the ensemble Kalman filter (EnKF) in the setting of finite dimensional signal evolution and noisy discrete-time observations. The signal dynamics is assumed to be governed by a stochastic differential equation (SDE), and a hierarchy of time grids is introduced for multilevel numerical integration of that SDE. Finally, the resulting multilevel EnKF is proved to asymptotically outperform EnKF in terms of computational cost versus approximation accuracy. The theoretical results are illustrated numerically.

ltering↗
DOE OSTI2024

Correlation function distributions for O ( N ) lattice field theories in the disordered phase

Numerical computations in strongly interacting quantum field theories are often performed using Monte Carlo sampling methods. A key task in these calculations is to estimate the value of a given physical quantity from the distribution of stochastic samples that are generated using the Monte Carlo method. Typically, the sample mean and sample variance are used to define the expectation values and uncertainties of computed quantities. However, the Monte Carlo sample distribution contains more information than these basic properties, and it is useful to investigate it more generally. In this work, the exact form of the probability distributions of two-point correlation functions at zero momentum in O ( N ) lattice field theories in the disordered phase and in infinite volume are determined. These distributions allow for a robust investigation of the efficacy of the Monte Carlo sampling procedure and are shown also to allow for improved estimators of the target physical quantity to be constructed. The theoretical expectations are shown to agree with numerical calculations in the O ( 2 ) model. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sensitivity Analyses for Monte Carlo Sampling-Based Particle Simulations

Computational design-based optimization is a well-used tool in science and engineering. Our report documents the successful use of a particle sensitivity analysis for design-based optimization within Monte Carlo sampling-based particle simulation—a currently unavailable capability. Such a capability enables the particle simulation communities to go beyond forward simulation and promises to reduce the burden on overworked analysts by getting more done with less computation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

MCCCS-MN

The MCCCS‒MN (Monte Carlo for Complex Chemical Systems‒Minnesota) software is developed by the Siepmann research group at the University of Minnesota. MCCCS‒MN allows for the simulation of multi-component molecular systems in the canonical, isobaric-isothermal (including constant stress for solids), grand-canonical, semi-grand, and Gibbs (NVT, NPT, and more than two simulation boxes) ensembles. It uses the configurational-bias Monte Carlo method to efficiently sample phase space for linear, branched and cyclic chain molecules, the adiabatic nuclear and electronic sampling Monte Carlo method to treat many-body polarization effects, and the aggregation-volume-bias Monte Carlo algorithm to efficiently sample the spatial distribution of associating molecules. MCCCS-MN employs a molecular representation of the system where force fields contain bonded and non-bonded terms. Funding for the development of MCCCS-MN through grants from the National Science Foundation (simulation of fluid phase equilibria and chromatography) and the Department of Energy (simulation of adsorption equilibria) is gratefully acknowledged.

Siepmann, J.Ilja [University of Minnesota - Twin C↗
DOE OSTI2026

Renovating Monte Carlo Methods and Codebases with Generative Models

The code will implement a standardized interface for Monte Carlo sampling methods, including conventional techniques, and going beyond current available packages to also incorporate generative model-enabled Monte Carlo sampling to provide a unified framewor

Garcia-Cardona, Cristina↗
DOE OSTI2019

Analytical Approach to Evaluate Uncertainty in Discharged Inventory Estimates at the Hanford Site - 19012

Many waste sites (>400) located within the Central Plateau portion of the Hanford Site received liquid discharges during the nuclear weapons material production mission (1944-1988). The Hanford Soil Inventory Model (SIM) is developed to quantitatively evaluate the discharged inventory for important radionuclides and chemicals at a given waste site to support fate and transport modeling of contaminants. An analytic approach was developed to estimate inventory and associated uncertainty. The approach was implemented for SIM and compared to a prior stochastic implementation applying a Monte-Carlo sampling approach. The analytic approach provides exact mean and variance results and the run times are much faster compared to the Monte-Carlo sampling approach. With the improved solution method using the analytical approach, it is feasible to run inventory updates and reduce associated uncertainty ranges. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

RESOLVE AND ECO: THE HALO MASS-DEPENDENT SHAPE OF GALAXY STELLAR AND BARYONIC MASS FUNCTIONS

In this work, we present galaxy stellar and baryonic (stars plus cold gas) mass functions (SMF and BMF) and their halo mass dependence for two volume-limited data sets. The first, RESOLVE-B, coincides with the Stripe 82 footprint and is extremely complete down to baryonic mass M {sub bary} ∼ 10{sup 9.1} M {sub ⊙}, probing the gas-rich dwarf regime below M {sub bary} ∼ 10{sup 10} M {sub ⊙}. The second, ECO, covers a ∼40× larger volume (containing RESOLVE-A) and is complete to M {sub bary} ∼ 10{sup 9.4} M {sub ⊙}. To construct the SMF and BMF we implement a new “cross-bin sampling” technique with Monte Carlo sampling from the full likelihood distributions of stellar or baryonic mass. Our SMFs exhibit the “plateau” feature starting below M {sub star} ∼ 10{sup 10} M {sub ⊙} that has been described in prior work. However, the BMF fills in this feature and rises as a straight power law below ∼10{sup 10} M {sub ⊙}, as gas-dominated galaxies become the majority of the population. Nonetheless, the low-mass slope of the BMF is not as steep as that of the theoretical dark matter halo MF. Moreover, we assign group halo masses by abundance matching, finding that the SMF and BMF, separated into four physically motivated halo mass regimes, reveal complex structure underlying the simple shape of the overall MFs. In particular, the satellite MFs are depressed below the central galaxy MF “humps” in groups with mass <10{sup 13.5} M {sub ⊙} yet rise steeply in clusters. Our results suggest that satellite destruction and stripping are active from the point of nascent group formation. We show that the key role of groups in shaping MFs enables reconstruction of a given survey’s SMF or BMF based on its group halo mass distribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Denoising Autoencoder for Reconstructing Sensor Observation Data and Predicting Evapotranspiration: Noisy and Missing Values Repair and Uncertainty Quantification

Abstract Machine learning (ML) methods applied in scientific research often deal with interrelated features in high‐dimensional data. Reducing data noise and redundancy is needed to increase prediction accuracy and efficiency especially when dealing with data from field sensors. We explored an unsupervised learning method, the denoising autoencoder (DAE), to extract the underlying data structure from noisy raw data in the context of predicting hydrologic quantities from multiple field sensors. These sensors have intrinsic instrumental noise and occasional malfunctions that cause missing values. Our DAE neural network reconstructed meteorological sensor data containing noise and missing values to predict evapotranspiration in a mountainous watershed. The DAE reconstructed the sensor variables with a mean coefficient of determination value of 0.77 across 15 dimensions representing individual sensors. It reduced variance and bias uncertainties compared to a classical autoencoder model. The reconstruction quality varied across dimensions depending on their cross‐correlation and alignment with the underlying data structure. Uncertainties arising from the model structure were overall higher than those resulting from data corruption. We attached the DAE structure to a downstream ET‐prediction neural network in three formats and achieved reasonably accurate ET predictions . The use of the DAE notably reduced variance uncertainty in ET prediction. However, excessive variance reduction may be accompanied by an increase in bias due to the intrinsic bias‐variance tradeoff. Our method of evaluating and reducing uncertainties in aggregated data from different sources can be used to improve predictive models, process understanding, and uncertainty quantification for better water resource management. Plain Language Summary We present a machine learning method, namely the denoising autoencoder, which reduces the effects of data noise and missing values typically present in scientific data sets collected through sensor measurements. This method selects the most relevant information from noisy raw data collected by the instruments and fills in missing values. To demonstrate the effectiveness of our method, we applied it to predict evapotranspiration, a hydrologic variable that represents the water moved from the land surface to the atmosphere through a combination of evaporation and plant water use (transpiration). We also used a random sampling technique (the Monte Carlo method) to compare the uncertainty in the predictions when using the raw and noisy data versus the reconstructed data. The denoising process produced more accurate predictions of evapotranspiration with less uncertainty. Improved predictions of evapotranspiration can lead to a better understanding and accounting of water budgets. This ML approach is broadly suitable for a wide variety of applications that involve noisy sensor data with missing values. Key Points We used a denoising autoencoder (DAE) neural network to reduce noise in meteorological and soil sensor observations by on average We used Monte Carlo sampling to estimate the bias and variance of all model outputs, including uncertainty sources from data and the model We attached the DAE component to a downstream neural network to predict ET with the variance reduced by , compared to that without the DAE

denoising autoencoder↗
DOE OSTI2016

An Atomic Mix Closure for Stochastic Media Transport Problems

The accurate modeling of stochastic media in radiation transport calculations remains on ongoing research challenge. The Levermore-Pomraning (LP) closure for the statistical transport equation is simple and inexpensive, but it is often inaccurate. More accurate closures have been proposed, but they can be prohibitively expensive. In this work we present a relatively inexpensive closure that we have observed to be more accurate than LP in almost all cases studied. Obtaining solutions to stochastic transport problems can be quite difficult. Transport calculations on an ensemble of explicit realizations generated by Monte Carlo sampling can be prohibitively expensive. Atomic mix or LP calculations are relatively inexpensive, but they can be inaccurate. The present work illustrates an approach that is less expensive than Monte Carlo sampling yet more accurate than atomic mix or LP approximations. We still need to analyze the reasons why the atomic mix closures derived for thicker problems lead to source iteration instabilities. It is not clear if the closure itself induces an ill-posed problem or if it is merely the iterative process that is problematic. We hope to examine alternative source terms to drive the subsidiary calculations. We also want to extend the work to multigroup and multidimensional problems. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Conditional Point Sampling: A Monte Carlo Method for Radiation Transport in Stochastic Media.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Progress in Normalizing Flows for 4d Gauge Theories

Normalizing flows have arisen as a tool to accelerate Monte Carlo sampling for lattice field theories. This work reviews recent progress in applying normalizing flows to 4-dimensional nonabelian gauge theories, focusing on two advancements: an architectural improvement referred to as learned active loops, and the application of correlated ensemble methods to QCD with N f = 2 dynamical fermions.

Abbott, Ryan [Massachusetts Institute of Technolog↗
DOE OSTI2017

Determination of critical experiment correlations for experiments involving arrays of low-enriched fuel rods

Quantifying correlations among critical experiments has become a topic of interest over the last 10-15 years. Correlations arise from the shared use of fissile material, non-fissile materials, fixtures, and procedures in multiple experiments or multiple experiment series. A Monte-Carlo sampling technique proposed as a viable technical method to evaluate these correlations has been implemented in the SCALE 6.2 Sampler sequence. This paper presents an investigation of the use of this new tool and technique for low-enriched uranium pin array experiments with a thermal neutron energy spectrum. The Nuclear Energy Agency's Expert Group on Uncertainty Analyses for Criticality Safety Assessments has been examining methods of calculating the correlation coefficients among experiments from the LCT-007 and LCT-039 series. Results for two scenarios are presented to demonstrate the impact of pin pitch uncertainty modeling on the correlation coefficients. The convergence of the correlation coefficients is also investigated. Additional studies examine the impact of the stochastic uncertainty and the use of continuous-energy Monte-Carlo for individual realizations on the correlation coefficients. Correlation coefficients are also determined for the LCT-042 experiment series. This experiment series differs from other experiments because unique poison panels and fuel array separation distances may reduce correlation of the individual experiments within the series. The effects of pin pitch uncertainty and stochastic uncertainty are also examined for this series. Comparing the results of the two sets of experiments shows potential conclusions that can be drawn regarding viability of the Monte-Carlo sampling technique for determining critical experiment correlations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Determination of critical experiment correlations using the sampler sequence within scale 6.2 - 14602

The validation of neutron transport methods used in nuclear criticality safety analyses is required by consensus American National Standards Institute/American Nuclear Society (ANSI/ANS) standards. In the last decade, there has been an increased interest in correlations among critical experiments used in validation that have shared physical attributes and which impact the independence of each measurement. The statistical methods included in many of the frequently cited guidance documents on performing validation calculations incorporate the assumption that all individual measurements are independent, so little guidance is available to practitioners on the topic. Typical guidance includes recommendations to select experiments from multiple facilities and experiment series in an attempt to minimize the impact of correlations or common-cause errors in experiments. Recent efforts have been made both to determine the magnitude of such correlations between experiments and to develop and apply methods for adjusting the bias and bias uncertainty to account for the correlations. This paper describes recent work performed at Oak Ridge National Laboratory using the Sampler sequence from the SCALE code system to develop experimental correlations using a Monte Carlo sampling technique. Sampler will be available for the first time with the release of SCALE 6.2, and a brief introduction to the methods used to calculate experiment correlations within this new sequence is presented in this paper. Techniques to utilize these correlations in the establishment of upper subcritical limits are the subject of a companion paper and will not be discussed here. Example experimental uncertainties and correlation coefficients are presented for a variety of low-enriched uranium water-moderated lattice experiments selected for use in a benchmark exercise by the Working Party on Nuclear Criticality Safety Subgroup on Uncertainty Analysis in Criticality Safety Analyses. The results include studies on the effect of fuel rod pitch on the correlations, and some observations are also made regarding difficulties in determining experimental correlations using the Monte Carlo sampling technique. (authors)

42 ENGINEERING↗
DOE OSTI2017

Determination of Critical Experiment Correlations for Experiments Involving Arrays of Low-Enriched Fuel Rods

Quantifying correlations among critical experiments has become a topic of interest over the last 10–15 years. Correlations arise from the shared use of fissile material, nonfissile materials, fixtures, and procedures in multiple experiments or multiple experiment series. A Monte Carlo sampling technique proposed as a viable technical method to evaluate these correlations has been implemented in the SCALE 6.2 Sampler sequence. This paper presents an investigation of the use of this new tool and technique for low-enriched uranium pin array experiments with a thermal neutron energy spectrum. The Nuclear Energy Agency’s Expert Group on Uncertainty Analyses for Criticality Safety Assessments has been examining methods of calculating the correlation coefficients among experiments from the LCT-007 and LCT-039 series. Results for two scenarios are presented to demonstrate the impact of pin pitch uncertainty modeling on the correlation coefficients. The convergence of the correlation coefficients is also investigated. Additional studies examine the impact of the stochastic uncertainty and the use of continuous-energy Monte Carlo for individual realizations on the correlation coefficients. Correlation coefficients are also determined for the LCT-042 experiment series. This experiment series differs from other experiments because unique poison panels and fuel array separation distances may reduce correlation of the individual experiments within the series. The effects of pin pitch uncertainty and stochastic uncertainty are also examined for this series. Comparing the results of the two sets of experiments shows potential conclusions that can be drawn regarding viability of the Monte Carlo sampling technique for determining critical experiment correlations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

A Bayesian approach to calibrating hydrogen flame kinetics using many experiments and parameters

First-principles Markov Chain Monte Carlo sampling is used to investigate uncertainty quantification and uncertainty propagation in parameters describing hydrogen kinetics. Specifically, we sample the posterior distribution for thirty-one parameters focusing on the H 2 O 2 and HO 2 reactions resulting from conditioning on ninety-one experiments. Established literature values are used for the remaining parameters in the mechanism as well as other thermodynamic and transport data needed to specify fluid properties. The samples are computed using an affine invariant sampler starting with broad, noninformative priors. Autocorrelation analysis shows that O(1M) samples are sufficient to obtain a reasonable sampling of the posterior. The resulting distribution identifies strong positive and negative correlations and several non-Gaussian characteristics. Using samples drawn from the posterior, we investigate the impact of parameter uncertainty on the prediction of two more complex flames: a 2D premixed flame kernel and the ignition of a hydrogen jet issuing into a heated chamber. The former represents a combustion regime similar to the target experiments used to calibrate the mechanism and the latter represents a different combustion regime. For the premixed flame, the net amount of product after a given time interval has a standard deviation of less than 2% whereas the standard deviation of the ignition time for the jet is more than 10%. The samples used for these studies are posted online. These results indicate the degree to which parameters consistent with the target experiments constrain predicted behavior in different combustion regimes. Finally, this process provides a framework for both identifying reactions for further study from candidate mechanisms as well as combining uncertainty quantification and propagation to, ultimately, tie uncertainty in laboratory flame experiments to uncertainty in end-use numerical predictions of more complicated scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Adversarial autoencoder ensemble for fast and probabilistic reconstructions of few-shot photon correlation functions for solid-state quantum emitters

Second-order photon correlation measurements [g (2) (τ) functions] are widely used to classify single-photon emission purity in quantum emitters or to measure the multiexciton quantum yield of emitters that can simultaneously host multiple excitations – such as quantum dots – by evaluating the value of g (2) (τ = 0). Accumulating enough photons to accurately calculate this value is time consuming and could be accelerated by fitting of few-shot photon correlations. Here, we develop an uncertainty-aware, deep adversarial autoencoder ensemble (AAE) that reconstructs noise-free g (2) (τ) functions from noise-dominated, few-shot inputs. The model is trained with simulated g (2) (τ) functions that are facilely generated by Poisson sampling time bins. The AAE reconstructions are performed orders-of-magnitude faster, with reconstruction errors and estimates of g (2) (τ = 0) that are lower in variance and similar in accuracy compared to Maximum likelihood estimation and Levenberg-Marquardt least-squares fitting approaches, for simulated and experimentally measured few-shot g (2) (τ) functions (~100 two-photon events) of InP/ZnS/ZnSe and CdS/CdSe/CdS quantum dots. The deep-ensemble model comprises eight individual autoencoders, allowing for probabilistic reconstructions of noise-free g (2) (τ) functions, and we show that the predicted variance scales inversely with number of shots, with comparable uncertainties to computationally intensive Markov chain Monte Carlo sampling. Furthermore, this work demonstrates the advantage of machine learning models to perform uncertainty-aware, fast, and accurate reconstructions of simple Poisson-distributed photon correlation functions, allowing for on-the-fly reconstructions and accelerated materials characterization of solid-state quantum emitters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗