Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)

This work describes the formalism for estimating thermonuclear reaction rates for astrophysical applications, emphasizing modern statistical approaches such as Monte Carlo sampling and Bayesian models. We discuss related topics including the calculation of resonance energies from nuclear Q values, indirect estimates of particle partial widths, and matching of reaction rates at elevated temperatures to statistical model results. We have evaluated available experimental data on cross sections, resonance energies and strengths, partial widths, lifetimes, spin-parities, and spectroscopic factors. Based on these results, we have estimated numerical values of 78 experimental charged-particle thermonuclear reaction rates for target nuclei in the A = 2–40 mass region, for temperatures ranging from 1 MK to 10 GK. For each reaction, three rate values are provided: low, median, and high, corresponding to the 16th, 50th, and 84th percentiles, respectively, of the cumulative reaction rate probability density distribution. Additionally, we present the factor uncertainty of each rate at each temperature grid point. These results enable users to sample the reaction rate probability density in nucleosynthesis calculations, facilitating uncertainty estimates of nuclidic abundances. The rates presented here refer to their laboratory values. For use in stellar model simulations, these values need to be corrected for the effects of thermal excitations of the interacting nuclei. For each reaction, we include graphs that illustrate the fractional contributions to the overall reaction rate along with the associated uncertainty. These visuals are designed to assist both stellar modelers and nuclear experimentalists by identifying the primary sources of rate uncertain=^texttx);ty at specific stellar temperatures. A graphical comparison with earlier Monte Carlo rates is also provided.

Nuclear astrophysics↗
DOE OSTI2024

A Diffusion‐Based Uncertainty Quantification Method to Advance E3SM Land Model Calibration

Abstract Calibrating land surface models and accurately quantifying their uncertainty are crucial for improving the reliability of simulations of complex environmental processes. This, in turn, advances our predictive understanding of ecosystems and supports climate‐resilient decision‐making. Traditional calibration methods, however, face challenges of high computational costs and difficulties in accurately quantifying parameter uncertainties. To address these issues, we develop a diffusion‐based uncertainty quantification (DBUQ) method. Unlike conventional generative diffusion methods, which are computationally expensive and memory‐intensive, DBUQ innovates by formulating a parameterized generative model and approximates this model through supervised learning, which enables quick generation of parameter posterior samples to quantify its uncertainty. DBUQ is effective, efficient, and general‐purpose, making it suitable for site‐specific ecosystem model calibration and broadly applicable for parameter uncertainty quantification across various earth system models. In this study, we applied DBUQ to calibrate the Energy Exascale Earth System Model land model at the Missouri Ozark AmeriFlux forest site. Results indicated that DBUQ produced accurate parameter posterior distributions similar to those from Markov Chain Monte Carlo sampling but with 30 times less computing time. This significant improvement in efficiency suggests that DBUQ can enable rapid, site‐level model calibration at a global scale, enhancing our predictive understanding of climate impacts on terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

VULCAN PLANETS: INSIDE-OUT FORMATION OF THE INNERMOST SUPER-EARTHS

The compact multi-transiting systems discovered by Kepler challenge traditional planet formation theories. These fall into two broad classes: (1) formation further out followed by migration and (2) formation in situ from a disk of gas and planetesimals. In the former, an abundance of resonant chains is expected, which the Kepler data do not support. In the latter, required disk mass surface densities may be too high. A recently proposed mechanism hypothesizes that planets form in situ at the pressure trap associated with the dead-zone inner boundary (DZIB) where radially drifting ''pebbles'' accumulate. This scenario predicts planet masses (M{sub p} ) are set by the gap-opening process that then leads to DZIB retreat, followed by sequential, inside-out planet formation (IOPF). For typical disk accretion rates, IOPF predictions for M{sub p} , M{sub p} versus orbital radius r, and planet-planet separations are consistent with observed systems. Here we investigate the IOPF prediction for how the masses, M{sub p,} {sub 1}, of the innermost (''Vulcan'') planets vary with r. We show that for fiducial parameters, M {sub p,} {sub 1} ≅ 5.0(r/0.1 AU) M {sub ⊕}, independent of the disk's accretion rate at time of planet formation. Then, using Monte Carlo sampling of a population of these innermost planets, we test this predicted scaling against observed planet properties, allowing for intrinsic dispersions in planetary densities and Kepler's observational biases. These effects lead to a slightly shallower relation M{sub p,} {sub 1}∝r {sup 0.9} {sup ±} {sup 0.2}, which is consistent with M{sub p,} 1∝r {sup 0.7} {sup ±} {sup 0.2} of the observed Vulcans. The normalization of the relation constrains the gap-opening process, favoring relatively low viscosities in the inner dead zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. Here, we perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Kekulé valence bond order in the honeycomb lattice optical Su-Schrieffer-Heeger model and its relevance to graphene

We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Forecasting Multi-Wave Epidemics Through Bayesian Inference

We present a simple, near-real-time Bayesian method to infer and forecast a multiwave outbreak, and demonstrate it on the COVID-19 pandemic. The approach uses timely epidemiological data that has been widely available for COVID-19. It provides short-term forecasts of the outbreak’s evolution, which can then be used for medical resource planning. The method postulates one- and multiwave infection models, which are convolved with the incubation-period distribution to yield competing disease models. The disease models’ parameters are estimated via Markov chain Monte Carlo sampling and information-theoretic criteria are used to select between them for use in forecasting. The method is demonstrated on two- and three-wave COVID-19 outbreaks in California, New Mexico and Florida, as observed during Summer-Winter 2020. We find that the method is robust to noise, provides useful forecasts (along with uncertainty bounds) and that it reliably detected when the initial single-wave COVID-19 outbreaks transformed into successive surges as containment efforts in these states failed by the end of Spring 2020.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

Abstract We describe an updated calibration and diagnostic framework, Balrog , used to directly sample the selection and photometric biases of the Dark Energy Survey (DES) Year 3 (Y3) data set. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant “source” galaxies and magnification biases of nearer “lens” galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ∼1–8 mmag, but for a small subset of objects, we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

No Galaxy Left Behind: Measuring the Transfer Function of the Dark Energy Survey with Balrog

In this dissertation, we describe a calibration and diagnostic framework called Balrog which was used to directly sample the selection and photometric biases of the Dark Energy Survey (DES) Year 3 (Y3) dataset. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with traditional generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant "source" galaxies and magnification biases of nearer "lens" galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the fiducial Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ~1-8 millimagnitudes, but for a small subset of objects we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. Finally, we discuss approaches to extend the current methodology to capture more aspects of the transfer function for future analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Use of a State-Level Integrated Assessment Model for Policy Analysis: An Uncertainty Analysis of Energy and Water

This presentation discusses some recent uncertainty assessment work being done on climate policy issues at the state level. The work is being done with the Platform for Regional Integrated Modeling and Analysis (PRIMA), a new climate change policy assessment framework developed at Pacific Northwest National Laboratory in Richland and Seattle, Washington and College Park, Maryland. The central computaional engine of PRIMA is the Global Change Assessment Model (GCAM), a dynamic-recursive integrated assessment model of the world economy, climate policy, and climate change that that has been regionalized to provide policy analysis down to the level of U.S. states. The PRIMA framework also includes detailed sectoral models that address issues in the energy-land-water interface geospatially at the sub-national level as the climate change and socioeconomic change occur. This presentation illustrates two uncertainty analyses conducted with GCAM and with water demand and supply components of PRIMA. The energy uncertainty analysis was conducted with GCAM alone and examines the impact of aggressive residential and commercial building equipment standards across the United States in the context of socioeconomic and climate uncertainty. The analysis focuses on the impact on future net energy service costs under aggressive standards in the residential and commercial building sectors from 2005 to 2050. The analysis uses a combination of discrete scenarios and Monte Carlo sampling approaches to address uncertainties. Results for Arizona and Florida are used to illustrate whether aggressive building equipment standards save energy and reduce net energy service costs when compared with current standards. Analysis of the output indicates a robust result under a broad range of future socioeconomic, policy, and climate scenarios. The water uncertainty analysis demonstrates a multi-model approach used to bound the prospects for future water shortages in Florida, given socioeconomic uncertainty and future water supply under a future climate. In this analysis, water demands from GCAM are coupled with a water management model that simulates regulated streamflows and water supply deficits given naturalized streamflows, water demands, and reservoir operating rules. The naturalized streamflows are simulated with a regional earth system model coupled to a river routing model. The climate is obtained from a regional climate projection of the Intergovernmental Panel on Climate Change (IPCC) Representative Concentration Pathway 8.5 (RCP 8.5) for all of the river basins in the United States. GCAM is used in Monte Carlo mode to estimate the outer bounds of water demand given socioeconomic uncertainty and RCP 8.5 climate. Future water demand is downscaled to a river basin level for Florida in high, low, and middle cases and is compared with available water supply for the same region for RCP 8.5 climate to identify the potential location and severity of water shortages.

Climate change, policy, mitigation, impacts, uncer↗
DOE OSTI2015

SU-E-T-144: Bayesian Inference of Local Relapse Data Using a Poisson-Based Tumour Control Probability Model

Purpose: The purpose of this project is to develop a robust method of parameter estimation for a Poisson-based TCP model using Bayesian inference. Methods: Bayesian inference was performed using the PyMC3 probabilistic programming framework written in Python. A Poisson-based TCP regression model that accounts for clonogen proliferation was fit to observed rates of local relapse as a function of equivalent dose in 2 Gy fractions for a population of 623 stage-I non-small-cell lung cancer patients. The Slice Markov Chain Monte Carlo sampling algorithm was used to sample the posterior distributions, and was initiated using the maximum of the posterior distributions found by optimization. The calculation of TCP with each sample step required integration over the free parameter α, which was performed using an adaptive 24-point Gauss-Legendre quadrature. Convergence was verified via inspection of the trace plot and posterior distribution for each of the fit parameters, as well as with comparisons of the most probable parameter values with their respective maximum likelihood estimates. Results: Posterior distributions for α, the standard deviation of α (σ), the average tumour cell-doubling time (Td), and the repopulation delay time (Tk), were generated assuming α/β = 10 Gy, and a fixed clonogen density of 10{sup 7} cm−{sup 3}. Posterior predictive plots generated from samples from these posterior distributions are in excellent agreement with the observed rates of local relapse used in the Bayesian inference. The most probable values of the model parameters also agree well with maximum likelihood estimates. Conclusion: A robust method of performing Bayesian inference of TCP data using a complex TCP model has been established.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Capabilities of multivariate Bayesian inference toward seismic hazard assessment

Multivariate Bayesian analysis can bring significant benefits to seismic hazard analysis: Its multivariate feature enables computing scalar and vector hazard without making any approximations; Correlations between intensity measures are implicitly modeled, permitting direct simulation of ground motion selection tools such as the conditional mean spectrum and the generalized conditioning intensity measure; and Its updating feature enables a seamless integration of new ground motion data into the hazard results. Here, we first develop a multivariate Bayesian ground motion model through the NGA-West2 database. The model functional form considers fault-type, magnitude, and distance dependencies, and also the linear and the rock intensity dependent site response. We use a hybrid Markov Chain Monte Carlo sampling to perform Bayesian inference consisting of Gibbs step and a multilevel Metropolis-Hastings step. We then perform several checks on the model and note that its performance is satisfactory. Finally, we illustrate the merits of this multivariate Bayesian analysis, which include: ground motion model updating with ground motion data recorded in the last four years not part of the NGA-West2 database; computation of scalar and vector seismic hazard using the un-updated and updated ground motion models for Los Angeles, CA; and simulation of the conditional mean spectrum under scalar and vector IM conditioning while accounting for different sources of aleatoric and epistemic uncertainties.

58 GEOSCIENCES↗
DOE OSTI2012

Mixed sneutrino dark matter in light of the 2011 XENON and LHC results

In the context of supersymmetric models in which small Dirac neutrino masses are generated by supersymmetry breaking, a mainly right-handed (RH) mixed sneutrino can be an excellent cold dark matter (DM) candidate. We perform a global analysis of the Minimal Supersymmetric Standard Model (MSSM)+RH neutrino parameter space by means of Markov Chain Monte Carlo sampling. We include all relevant constraints from collider and dark matter searches, paying particular attention to nuclear and astrophysical uncertainties. Two distinct cases can satisfy all constraints: heavy sneutrino DM with mass of order 100 GeV, as well as light sneutrino DM with mass of about 3–6 GeV. We discuss the implications for direct and indirect dark matter searches, as well as for SUSY and Higgs searches at the LHC for both, the light and the heavy sneutrino dark matter case. The light sneutrino case is excluded by the 125–126 GeV Higgs signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Multiscale Mathematics for Biomass Conversion to Renewable Hydrogen

Our two key accomplishments in the first three years were towards the development of, (1) a mathematically rigorous and at the same time computationally flexible framework for parallelization of Kinetic Monte Carlo methods, and its implementation on GPUs, and (2) spatial multilevel coarse-graining methods for Monte Carlo sampling and molecular simulation. A common underlying theme in both these lines of our work is the development of numerical methods which are at the same time both computationally efficient and reliable, the latter in the sense that they provide controlled-error approximations for coarse observables of the simulated molecular systems. Finally, our key accomplishment in the last year of the grant is that we started developing (3) pathwise information theory-based and goal-oriented sensitivity analysis and parameter identification methods for complex high-dimensional dynamics and in particular of nonequilibrium extended (high-dimensional) systems. We discuss these three research directions in some detail below, along with the related publications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

DISSECTING MAGNETAR VARIABILITY WITH BAYESIAN HIERARCHICAL MODELS

Neutron stars are a prime laboratory for testing physical processes under conditions of strong gravity, high density, and extreme magnetic fields. Among the zoo of neutron star phenomena, magnetars stand out for their bursting behavior, ranging from extremely bright, rare giant flares to numerous, less energetic recurrent bursts. The exact trigger and emission mechanisms for these bursts are not known; favored models involve either a crust fracture and subsequent energy release into the magnetosphere, or explosive reconnection of magnetic field lines. In the absence of a predictive model, understanding the physical processes responsible for magnetar burst variability is difficult. Here, we develop an empirical model that decomposes magnetar bursts into a superposition of small spike-like features with a simple functional form, where the number of model components is itself part of the inference problem. The cascades of spikes that we model might be formed by avalanches of reconnection, or crust rupture aftershocks. Using Markov Chain Monte Carlo sampling augmented with reversible jumps between models with different numbers of parameters, we characterize the posterior distributions of the model parameters and the number of components per burst. We relate these model parameters to physical quantities in the system, and show for the first time that the variability within a burst does not conform to predictions from ideas of self-organized criticality. We also examine how well the properties of the spikes fit the predictions of simplified cascade models for the different trigger mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Constraining white dwarf structure and neutrino physics in 47 Tucanae

We present a robust statistical analysis of the white dwarf cooling sequence in 47 Tucanae. We combine Hubble Space Telescope UV and optical data in the core of the cluster, Modules for Experiments in Stellar Evolution (MESA) white dwarf cooling models, white dwarf atmosphere models, artificial star tests, and a Markov Chain Monte Carlo sampling method to fit white dwarf cooling models to our data directly. We use a technique known as the unbinned maximum likelihood to fit these models to our data without binning. We use these data to constrain neutrino production and the thickness of the hydrogen layer in these white dwarfs. The data prefer thicker hydrogen layers (q{sub H}=3.2×10{sup −5}) and we can strongly rule out thin layers (q{sub H}=10{sup −6}). The neutrino rates currently in the models are consistent with the data. This analysis does not provide a constraint on the number of neutrino species.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Determination of Resonance Parameters and their Covariances from Neutron Induced Reaction Cross Section Data

Cross section data in the resolved and unresolved resonance region are represented by nuclear reaction formalisms using parameters which are determined by fitting them to experimental data. Therefore, the quality of evaluated cross sections in the resonance region strongly depends on the experimental data used in the adjustment process and an assessment of the experimental covariance data is of primary importance in determining the accuracy of evaluated cross section data. In this contribution, uncertainty components of experimental observables resulting from total and reaction cross section experiments are quantified by identifying the metrological parameters involved in the measurement, data reduction and analysis process. In addition, different methods that can be applied to propagate the covariance of the experimental observables (i.e. transmission and reaction yields) to the covariance of the resonance parameters are discussed and compared. The methods being discussed are: conventional uncertainty propagation, Monte Carlo sampling and marginalization. It is demonstrated that the final covariance matrix of the resonance parameters not only strongly depends on the type of experimental observables used in the adjustment process, the experimental conditions and the characteristics of the resonance structure, but also on the method that is used to propagate the covariances. Finally, a special data reduction concept and format is presented, which offers the possibility to store the full covariance information of experimental data in the EXFOR library and provides the information required to perform a full covariance evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Sensitivity and uncertainty analysis for a combined depletion and criticality calculation chain of a BWR FA

This paper examines the prerequisites and boundary conditions for extending burnup credit from pressurized water reactors to boiling water reactors with validated depletion codes by performing uncertainty and sensitivity analyses on a modern boiling water reactor fuel assembly design. We developed a computational chain consisting of an axial boiling water reactor axial irradiation model, followed by a series of depletion calculations and a three-dimensional criticality calculation. The irradiation model puts into relation the profiles of power, burnup and moderator density of the fuel assembly at 24 axial locations. In a simplified fashion it also includes the insertion of control blades. The depletion calculations of the fuel assembly were performed in axial slices with the 2D-depletion code HELIOS 1.12. The criticality calculations were done in 3D-geometry with KENO-Va from the SCALE 6.2.1 package, in which several fuel assemblies are placed into a generic transport cask. For the uncertainty and sensitivity analyses a Monte Carlo sampling based approach is used applying the GRS tool SUnCISTT. We performed several generic investigations testing the capabilities and the limitations of the whole computational chain. The impact of variations of the irradiation model on the resulting reactivity of the fuel assemblies inside a transport cask or storage pool was investigated using an uncertainty and sensitivity analysis on the whole computational chain. Geometrical uncertainties were also included in the evaluation and a parameter study on the use of control blades during irradiation was performed. The goal was to develop the procedures and tools to be able to obtain quantitative uncertainty data which can be used in the future to derive conservative irradiation profiles without unnecessary excess reactivity margins. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

A Flexible Transition State Searching Method for Atmospheric Reaction Systems

The precise and rapid exploration of transition states (TSs) is a major challenge when studying atmospheric reactions due to their complexity. In this work, a Monte Carlo Transition State Search Method (MCTSSM), which integrates Monte Carlo sampling technique with transition state optimization methods using an efficient computer script, has been developed for transition state searches. The efficiency and the potential application in atmospheric reactions of this method have been demonstrated by three types of test suits related to the reactions of atmospheric volatile organic compounds (VOCs): (1) OH addition, (2) OH hydrogen-abstraction, and (3) the other reactive group (e.g. Cl, O3, NO3), especially for the reaction of β-pinene-sCI (stabilized Criegee Intermediates) with water. It was shown that the application of this method with effective restricted parameters has greatly simplified the time-consuming and tedious manual search procedure for transition state (TS) of the bimolecular reaction systems.

Aerosol↗