Monte Carlo Renormalization Group for Classical Lattice Models with Quenched Disorder
Not provided.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SEARCH · Engineering Papers
Results for “Monte Carlo Modelling”
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.
Quantitative assessment of the universal thermopower in the Hubbard model
As primarily an electronic observable, the room-temperature thermopower S in cuprates provides possibilities for a quantitative assessment of the Hubbard model. Using determinant quantum Monte Carlo, we demonstrate agreement between Hubbard model calculations and experimentally measured room-temperature S across multiple cuprate families, both qualitatively in terms of the doping dependence and quantitatively in terms of magnitude. We observe an upturn in S with decreasing temperatures, which possesses a slope comparable to that observed experimentally in cuprates. From our calculations, the doping at which S changes sign occurs in close proximity to a vanishing temperature dependence of the chemical potential at fixed density. Our results emphasize the importance of interaction effects in the systematic assessment of the thermopower S in cuprates.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays
Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
The utility of Bayesian data reconciliation for separations
Data reconciliation methods for separation processes typically rely on classical statistical approaches to generate estimates of true mass flow rates from measurements. Knowledge regarding the uncertainty of these estimates has value in decision making, but is often not acquired. Bayesian approaches intrinsically quantify uncertainty; however, literature for Bayesian data reconciliation of separation processes is scarce. This publication outlines two Bayesian data reconciliation models and provides details for how the models were implemented for the BayesMassBal (V 1.0.0) software package written in R. To demonstrate the advantages of this approach for data reconciliation, the models were first applied to simulated data and then compared to a classical model through a Monte Carlo experiment. In this example, the Bayesian models were found to provide more accurate estimates of the simulated data, while also providing quantitative information on the estimate uncertainty. To demonstrate the use of the technique in a practical problem, the models were also applied to real data collected from a pilot-scale rare earth solvent extraction process. Here, this publication provides a small window into how Bayesian methods can be used for data reconciliation, but findings suggest Bayesian data reconciliation models for separation processes have distinct advantages over classical alternatives.
01 COAL, LIGNITE, AND PEAT↗
Context-aware learning of hierarchies of low-fidelity models for multi-fidelity uncertainty quantification
Multi-fidelity Monte Carlo methods leverage low-fidelity and surrogate models for variance reduction to make tractable uncertainty quantification even when numerically simulating the physical systems of interest with high-fidelity models is computationally expensive. This work proposes a context-aware multi-fidelity Monte Carlo method that optimally balances the costs of training low-fidelity models with the costs of Monte Carlo sampling. It generalizes the previously developed context-aware bi-fidelity Monte Carlo method to hierarchies of multiple models and to more general types of low-fidelity models. When training low-fidelity models, the proposed approach takes into account the context in which the learned low-fidelity models will be used, namely for variance reduction in Monte Carlo estimation, which allows it to find optimal trade-offs between training and sampling to minimize upper bounds of the mean-squared errors of the estimators for given computational budgets. This is in stark contrast to traditional surrogate modeling and model reduction techniques that construct low-fidelity models with the primary goal of approximating well the high-fidelity model outputs and typically ignore the context in which the learned models will be used in upstream tasks. Further, the proposed context-aware multi-fidelity Monte Carlo method applies to hierarchies of a wide range of types of low-fidelity models such as sparse-grid and deep-network models. Numerical experiments with the gyrokinetic simulation code Gene show speedups of up to two orders of magnitude compared to standard estimators when quantifying uncertainties in small-scale fluctuations in confined plasma in fusion reactors. This corresponds to a runtime reduction from 72 days to four hours on one node of the Lonestar6 supercomputer at the Texas Advanced Computing Center.
42 ENGINEERING↗
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↗
Biexciton Condensation in Electron-Hole-Doped Hubbard Bilayers: A Sign-Problem-Free Quantum Monte Carlo Study
The bilayer Hubbard model with electron-hole doping is an ideal platform to study excitonic orders due to suppressed recombination via spatial separation of electrons and holes. However, suffering from the sign problem, previous quantum Monte Carlo studies could not arrive at an unequivocal conclusion regarding the presence of phases with clear signatures of excitonic condensation in bilayer Hubbard models. Here, we develop a determinant quantum Monte Carlo algorithm for the bilayer Hubbard model that is sign-problem-free for equal and opposite doping in the two layers and study excitonic order and charge and spin density modulations as a function of chemical potential difference between the two layers, on-site Coulomb repulsion, and interlayer interaction. In the intermediate coupling regime and in proximity to the SU(4)-symmetric point, we find a biexcitonic condensate phase at finite electron-hole doping, as well as a competing (π, π) charge density wave state. In this study, we extract the Berezinskii-Kosterlitz-Thouless transition temperature from superfluid density and a finite-size scaling analysis of the correlation functions and explain our results in terms of an effective biexcitonic hard-core boson model.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Enhancing fatigue life by ductile-transformable multicomponent B2 precipitates in a high-entropy alloy
Catastrophic accidents caused by fatigue failures often occur in engineering structures. Thus, a fundamental understanding of cyclic-deformation and fatigue-failure mechanisms is critical for the development of fatigue-resistant structural materials. Here we report a high-entropy alloy with enhanced fatigue life by ductile-transformable multicomponent B2 precipitates. Its cyclic-deformation mechanisms are revealed by real-time in-situ neutron diffraction, transmission-electron microscopy, crystal-plasticity modeling, and Monte-Carlo simulations. Multiple cyclic-deformation mechanisms, including dislocation slips, precipitation strengthening, deformation twinning, and reversible martensitic phase transformation, are observed in the studied high-entropy alloy. Its improved fatigue performance at low strain amplitudes, i.e., the high fatigue-crack-initiation resistance, is attributed to the high elasticity, plastic deformability, and martensitic transformation of the B2-strengthening phase. This study shows that fatigue-resistant alloys can be developed by incorporating strengthening ductile-transformable multicomponent intermetallic phases.
36 MATERIALS SCIENCE↗
Off-Lattice Markov Chain Monte Carlo Simulations of Mechanically Driven Polymers
Here, we develop off-lattice simulations of semiflexible polymer chains subjected to applied mechanical forces by using Markov Chain Monte Carlo. Our approach models the polymer as a chain of fixed length bonds, with configurations updated through adaptive nonlocal Monte Carlo moves. This proposed method enables precise calculation of a polymer’s response to a wide range of mechanical forces, which traditional on-lattice models cannot achieve. Our approach has shown excellent agreement with theoretical predictions of persistence length and end-to-end distance in quiescent states as well as stretching distances under tension. Moreover, our model eliminates the orientational bias present in on-lattice models, which significantly impacts calculations such as the scattering function, a crucial technique for revealing the polymer conformation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo
Stochastic reaction network models are often used to explain and predict the dynamics of gene regulation in single cells.These models usually involve several parameters, such as the kinetic rates of chemical reactions, that are not directly measurable and must be inferred from experimental data. Bayesian inference provides a rigorous probabilistic frame-work for identifying these parameters by finding a posterior parameter distribution that captures their uncertainty.Traditional computational methods for solving inference problems such as Markov Chain Monte Carlo methods based on classical Metropolis-Hastings algorithm involve numerous serial evaluations of the likelihood function, which in turn requires expensive forward solutions of the chemical master equation (CME). We propose an alternate approach based on a multifidelity extension of the Sequential Tempered Markov Chain Monte Carlo (ST-MCMC) sampler. This algorithm is built upon Sequential Monte Carlo and solves the Bayesian inference problem by decomposing it into a sequence of efficiently solved subproblems that gradually increase both model fidelity and the influence of the observed data. We reformulate the finite state projection (FSP) algorithm, a well-known method for solving the CME, to produce a hierarchy of surrogate master equations to be used in this multifidelity scheme. To determine the appropriate fidelity, we introduce a novel information-theoretic criteria that seeks to extract the most information about the ultimate Bayesian posterior from each model in the hierarchy without inducing significant bias. This novel sampling scheme is tested with high performance computing resources using biologically relevant problems.
97 MATHEMATICS AND COMPUTING↗
Joint Bayesian Inference for Near-Surface Explosion Yield and Height-of-Burst
Forensic capabilities to understand chemical and nuclear explosions are greatly aided by an accurate estimate of explosive yield with uncertainty. The relationship between explosive size and geophysical observations of seismic, acoustic, and optical waves can be exploited to provide an estimate of yield. Any near-surface yield estimate is complicated by the surface interaction, so an estimate for the explosion height-of-burst is necessarily included in the relationship. Additionally, the relationship dictates a trade-off between estimates of yield and height-of-burst. Fortunately, the surface interaction for each type of observation is different, which breaks the trade-off, and the inclusion of height-of-burst with multiple data types improves yield estimation. We define simple parametric forward models to relate seismoäcoustoöptic observations from a data set of known explosive yields and height-of-bursts. The parameters of the models and a prediction for the yield and height-of-burst of a new event can then be estimated given new observations via Bayesian inference. We report posterior distribution estimates of the parametric models using a Markov chain Monte Carlo sampling technique. These models are then used to predict the yield and height-of-burst of SUGAR, a historical near-surface nuclear explosion, using its reported historical observations. The reported yield of 1.2 ktonne Trinitrotoluene (TNT)-equivalent (Department of Energy, 2015) is within the estimated posterior. Yield uncertainty can be estimated from the spread of the posterior, which is between 0.9 and 2.1 ktonne TNT-equivalent. The posterior for height-of-burst has a wider range between 10 m below and 8 m above ground that includes the true height-of-burst of 1 m.
58 GEOSCIENCES↗
Quantifying airborne fraction trends and the destination of anthropogenic CO 2 by tracking carbon flows in a simple climate model
Abstract Atmospheric carbon dioxide (CO 2 ) concentrations have increased as a direct result of human activity and are at their highest level over the last 2 million years, with profound impacts on the Earth system. However, the magnitude and future dynamics of land and ocean carbon sinks are not well understood; therefore, the amount of anthropogenic fossil fuel emissions that remain in the atmosphere (the airborne fraction) is poorly constrained. This work aims to quantify the sources and controls of atmospheric CO 2 , the fate of anthropogenic CO 2 over time, and the likelihood of a trend in the airborne fraction. We use Hector v3.0, a coupled simple climate and carbon cycle model with the novel ability to explicitly track carbon as it flows through the Earth system. We use key model parameters in a Monte Carlo analysis of 15 000 model runs from 1750 to 2300. Results are filtered for physical realism against historical observations and CMIP6 projection data, and we calculate the relative importance of parameters controlling how much anthropogenic carbon ends up in the atmosphere. Modeled airborne fraction was roughly 52%, consistent with observational studies. The overwhelming majority of model runs exhibited a negative trend in the airborne fraction from 1960–2020, implying that current-day land and ocean sinks are proportionally taking up more carbon than the atmosphere. However, the percentage of atmospheric CO 2 derived from anthropogenic origins can be much higher because of Earth system feedbacks. We find it peaks at over 90% between 2010–2050. Moreover, when looking at the destination of anthropogenic fossil fuel emissions, only a quarter ends up in the atmosphere while more than half of emissions are taken up by the land sink on centennial timescales. This study evaluates the likelihood of airborne fraction trends and provides insights into the dynamics of anthropogenic CO 2 in the Earth system.
54 ENVIRONMENTAL SCIENCES↗
Null Space Monte Carlo Evaluation of the Plateau to River Model
The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Deep Neural Network Informed Markov Chain Monte Carlo Methods
In subsurface flow modeling, quantifying the uncertainty of model parameters and the corresponding uncertainly on output quantities is a crucial task for groundwater management. Markov chain Monte Carlo (MCMC) methods can take advantage of observed data to estimate parameters in a Bayesian setting. However, MCMC can be slow to converge and produce highly correlated samples when the dimensions of the parameters is high. Using gradients for the posterior distribution can help samplers explore the parameter space more efficiently, but obtaining gradients can be computationally challenging.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo
We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.
08 HYDROGEN↗
Machine learning-assisted profiling of a kinked ladder polymer structure using scattering
Ladder polymers consisting of fused rings in the backbone have very limited conformational freedom, which results in very different properties from traditional linear polymers. However, accurately determining their size and chain conformations from solution scattering remains a challenge. Their chain conformations of kinked ladder polymers are largely governed by the structures and relative orientations or configurations of the repeat units, unlike conventional polymer chains whose bending angles between repeat units follow a unimodal Gaussian distribution. Meanwhile, traditional scattering models for polymer chains do not account for these unique structural features. This work introduces a novel approach that integrates machine learning with Monte Carlo simulations to construct a model that can describe the geometry of a type of kinked CANAL ladder polymers. We first develop a Monte Carlo simulation model for sampling the configuration space of CANAL ladder polymers, where each repeat unit is modeled as a biaxial segment. Then, we establish a machine learning-assisted scattering analysis framework based on Gaussian Process Regression. Finally, we conduct small-angle neutron scattering experiments on a CANAL ladder polymer solution to apply our approach. Our method uncovers structural features of such ladder polymers that conventional methods fail to capture.
Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes
Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.
97 MATHEMATICS AND COMPUTING↗
NETL Natural Gas Lifecycle Model
This is the excel-based life cycle model that contains all the parameters and Monte Carlo simulation capabilities to model the techno-regions contained in the report: Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile.
conventional gas↗