SEARCH · Search NASA
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.
Monte carlo method of finding the radiative heat transfer and emissive power distribution between concentric cylinders enclosing a nonisothermal gas with a uniformly distributed heat source
Monte carlo method of finding radiative heat transfer and emissive power distribution between concentric cylinders enclosing nonisothermal gas with uniformly distributed heat source
Not Available
Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. Here, we investigate the compatibility, strengths, and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multideterminant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis-set full configuration interaction energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite–size solid-state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis-set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wave function than DMC while simultaneously exhibiting a much larger basis-set dependence. We outline some of the remaining challenges and opportunities for improving these approaches.
Monte Carlo transport methods are the most accurate schemes for solving problems with complex energy and spatial features, but they come with a high computational cost. Although hybrid methods have enabled the use of Monte Carlo transport for a large class of problems, they still require significant computing resources. Modern multicore CPUs with large numbers of compute cores and graphical processing units (GPUs) provide opportunities to optimize the memory and run-time costs of hybrid Monte Carlo methods. This paper documents the development and analysis of three Monte Carlo transport algorithms that support hybrid transport using the consistent adjoint-driven importance sampling (CADIS) and forward-weighted CADIS methods in the Shift Monte Carlo code: history-based transport using static and dynamic threading on multicore CPUs and event-based transport enabling weight window tracking on GPUs. The results are shown for two challenging hybrid problems on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility. The results show that all three methods yield good performance and enable solutions of difficult fixed-source transport problems in less than 2 min on 20 nodes of Frontier. Dynamic threading was observed to give up to 20% better scaling behavior than static threading. Moreover, the AMD Instinct 250X GPU was found to give 9 to 11 times greater throughput per graphics compute die than the best CPU performance. In conclusion, additional opportunities for optimization of hybrid transport on GPUs are discussed.
As an important Markov chain Monte Carlo (MCMC) method, the stochastic gradient Langevin dynamics (SGLD) algorithm has achieved great success in Bayesian learning and posterior sampling. Furthermore, SGLD typically suffers from a slow convergence rate due to its large variance caused by the stochastic gradient. In order to alleviate these drawbacks, we leverage the recently developed Laplacian smoothing technique and propose a Laplacian smoothing stochastic gradient Langevin dynamics (LS-SGLD) algorithm. We prove that for sampling from both log-concave and non-log-concave densities, LS-SGLD achieves strictly smaller discretization error in 2-Wasserstein distance, although its mixing rate can be slightly slower. Experiments on both synthetic and real datasets verify our theoretical results and demonstrate the superior performance of LS-SGLD on different machine learning tasks including posterior sampling, Bayesian logistic regression, and training Bayesian convolutional neural networks.
In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schrodinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, a Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.
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.
The Monte Carlo method delivers an unbiased estimate of the probability of failure. However, the variance of the estimate depends on the number of evaluated samples. This number must be very large for estimations of a low probability of failure. If the evaluation of each sample is computationally expensive, the crude Monte Carlo simulation strategy is impracticable. Therefore, subset simulations are used to reduce the required number of evaluations. Subset simulations require a Markov Chain Monte Carlo sampler, such as the random walk Metropolis-Hastings algorithm. The algorithm, however, struggles with sampling in low-probability regions, especially if they are narrow. As a consequence, advanced Markov Chain Monte Carlo simulations have been developed. In particular, the Hamiltonian Monte Carlo method explores the target distribution rapidly. Driven by the idea of Hamiltonian dynamics, this sampler provides a non-random walk through the target distribution. The incorporation of subset simulation and Hamiltonian Monte Carlo methods has shown promising results for reliability analysis. One downside of the Hamiltonian Monte Carlo method is that gradient evaluations are computationally expensive, especially when dealing with high-dimensional problems and evaluating long trajectories. We show that integrating Hamiltonian neural networks in Hamiltonian Monte Carlo simulations significantly speeds up the sampling task. Furthermore, the enhancement of adaptive trajectory length within the Hamiltonian Monte Carlo results in the efficient proposal of the following states. Based on this recent enhancement, we provide a fast sampling strategy for subset simulations using Hamiltonian neural networks to replace the evaluation of the gradient and significantly speed up the Hamiltonian Monte Carlo simulation.
Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.
Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy
Monte carlo method applied to hydrogen diffusion in thermosphere
First-Passage Kinetic Monte Carlo (FPKMC) simulations of species migration in duplex stainless steels were performed in order to establish relationships between alloy content, segregation behavior, and the formation of Ni-Si-Mn rich particles in cast duplex stainless steels during thermal aging. The Ni-Si-Mn-rich second phase forms after extended aging at reactor operating temperatures and degrades alloy properties. Simulations of Ni-Si-Mn cluster formation were validated through comparison with experimental results obtained through Atom Probe Tomography (APT) on similar alloy compositions, identifying several trends. First, Cu promotes the formation of Ni-Si-Mn clusters, but only when Ni or Mn prefer segregation to the surface of Cu particles; without the segregation of these species, the critical composition for the clusters to form was not achieved. Second, the width of precipitate-denuded zones near ?/d interfaces increases with decreasing Cu content. This finding was in strong agreement with APT results, further validating our model. Finally, our model predicts that Ni-Si-Mn cluster formation will be the most extensive when the Si:Mn ratio is approximately 1:1 and the least extensive when one of these key elements is less concentrated (Mn = 0.5 at.% or Si = 0.5 at.%). The implications of these results on how to improve alloy properties are discussed.
The NPDGamma experiment measures the parity-violating asymmetry in gamma-ray emission in the capture of polarized neutrons on liquid parahydrogen. The sensitivity to the asymmetry for each detector in the array is used as a parameter in the extraction of the physics asymmetry from the measured data. The detector array is approximately cylindrically symmetric around the target and a step-wise sinusoidal function has been used for the sensitivity in the previous iteration of the NPDGamma experiment, but deviations from cylindrical symmetry necessitate the use of a Monte Carlo model to determine corrections to the geometrical factors. For the calculations, source code modifications to MCNPX were done in order to calculate the sensitivity of each cesium iodide detector to the physics asymmetry. We describe the MCNPX model and results from calculations and how the results are validated through measurement of the parity violating asymmetry of gamma-rays from neutron capture on chlorine.
Monte Carlo analysis of rarefied-gas heat transfer between parallel plates under various flow conditions
Monte Carlo method analysis of rarefied gas heat transfer between parallel plates in terms of temperature, density and Knudsen number
Atomic-scale insights into the dissolution mechanisms of borosilicate glasses are critical to the development of the waste form corrosion models that inform performance assessments of nuclear waste disposal facilities. Monte Carlo (MC) simulations have been instrumental in advancing the state of knowledge of borosilicate glass dissolution; however, all such simulations to date have relied on lattices to represent the amorphous glass structure, thereby potentially introducing artifacts. Therefore, a new approach, referred to as the amorphous MC approach, has been developed whereby glass structures generated from molecular dynamics simulations are used as starting points for MC simulations of glass dissolution. MC simulations of sodium borosilicate glasses that covered a wide compositional range were performed and revealed correlations between dissolution rate and structural features of the glasses. Notably, the lattice MC approach predicted the fraction of tetrahedral boron to have a lesser influence on the dissolution rate than computed with the amorphous MC approach. The lattice MC simulations also overestimated the dissolution rate for a given glass composition. These findings highlight the importance of using truly amorphous structures in MC simulations of glass corrosion.