Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗
DOE OSTI2015

Molecular-Atomic Transition along the Deuterium Hugoniot Curve with Coupled Electron-Ion Monte Carlo Simulations

We have performed accurate simulations of the Deuterium Hugoniot using Coupled Electron Ion Monte Carlo (CEIMC). Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the principal Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we found a maximum compression of 4.85, roughly 10% larger than previous theoretical predictions and experimental data but still compatible with the latter because of their large uncertainty.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Landau-Level Mixing and Particle-Hole Symmetry Breaking for Spin Transitions in the Fractional Quantum Hall Effect

Here, the spin transitions in the fractional quantum Hall effect provide a direct measure of the tiny energy differences between differently spin-polarized states, and thereby serve as an extremely sen- sitive test of the quantitative accuracy of the theory of the fractional quantum Hall effect, and, in particular, of the role of Landau-level mixing in lifting the particle-hole symmetry. We report on an accurate quantitative study of this physics, evaluating the effect of Landau-level mixing in a nonperturbative manner using a fixed-phase diffusion Monte Carlo method. We find excellent agree- ment between our calculated critical Zeeman energies and the experimentally measured values. In particular, we find, as also do experiments, that the critical Zeeman energies for fractional quantum Hall states at filling factors ν = 2–n/(2n±1) are significantly higher than those for ν = n/(2n±1), a quantitative signature of the lifting of particle-hole symmetry due to Landau-level mixing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Variational Approach to Monte Carlo Renormalization Group

We present a Monte Carlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. Here, the two-dimensional Ising model is used to illustrate the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Crystallization in the Fractional Quantum Hall Regime Induced by Landau-Level Mixing

Here, the interplay between strongly correlated liquid and crystal phases for two-dimensional electrons exposed to a high transverse magnetic field is of fundamental interest. Through the non-perturbative fixed phase diffusion Monte Carlo method, we determine the phase diagram of the Wigner crystal in the ν – κ plane, where ν is the filling factor and κ is the strength of Landau level mixing. The phase boundary is seen to exhibit a striking ν dependence, with the states away from the magic filling factors ν = n/(2pn + 1) being much more susceptible to crystallization due to Landau level mixing than those at ν = n/(2pn + 1). Our results explain the qualitative difference between the experimental behaviors observed in n-doped and p-doped GaAs quantum wells, and, in particular, the existence of an insulating state for ν < 1/3 and also for 1/3 < ν < 2/5 in low density p-doped systems. We predict that in the vicinity of ν = 1/5 and ν = 2/9, increasing LL mixing causes a transition not into an ordinary electron Wigner crystal but rather into a strongly correlated crystal of composite fermions carrying two vortices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Variational Excitations in Real Solids: Optical Gaps and Insights into Many-Body Perturbation Theory

We present an approach to studying optical band gaps in real solids in which quantum Monte Carlo methods allow for the application of a rigorous variational principle to both ground and excited state wave functions. In tests that include small, medium, and large band gap materials, optical gaps are predicted with a mean absolute deviation of 3.5% against experiment, less than half the equivalent errors for typical many-body perturbation theories. The approach is designed to be insensitive to the choice of density functional, a property we exploit in order to provide insight into how far different functionals are from satisfying the assumptions of many-body perturbation theory. Here, we explore this question most deeply in the challenging case of ZnO, where we show that, although many commonly used functionals have shortcomings, there does exist a one-particle basis in which perturbation theory’s zeroth-order picture is sound. Insights of this nature should be useful in guiding the future application and improvement of these widely used techniques.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Here we present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Ab Initio Nuclear Thermodynamics

Here, we propose a new Monte Carlo method called the pinhole trace algorithm for ab initio calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first ab initio study of the density and temperature dependence of nuclear clustering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. In conclusion, our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Distilling the Essential Elements of Nuclear Binding via Neural-Network Quantum States

To distill the essential elements of nuclear binding, we seek the simplest Hamiltonian capable of modeling atomic nuclei with percent-level accuracy. A critical aspect of this endeavor consists of accurately solving the quantum many-body problem without incurring an exponential computing cost with the number of nucleons. Furthermore, we address this challenge by leveraging a variational Monte Carlo method based on a highly expressive neural-network quantum state ansatz. In addition to computing binding energies and charge radii of nuclei with up to 𝐴 = 20 nucleons, by evaluating their magnetic moments, we demonstrate that neural-network quantum states are able to correctly capture the self-emerging nuclear shell structure. To this end, we introduce a novel computational protocol based on adding an external magnetic field to the nuclear Hamiltonian, which allows the neural network to learn the preferred polarization of the nucleus within the given magnetic field.

Binding energy & masses↗
DOE OSTI2018

Dimensional control of defect dynamics in perovskite oxide superlattices

Point defects play a critical role in the structural, physical, and interfacial properties of perovskite oxide superlattices. However, understanding of the fundamental properties of point defects in superlattices, especially their transport properties, is rather limited. Here, we report predictions of the stability and dynamics of oxygen vacancies in SrTiO 3 /PbTiO 3 oxide superlattices using first-principles calculations in combination with the kinetic Monte Carlo method. By varying the stacking period, i.e., changing of n in n STO/ n PTO, we discover a crossover from three-dimensional diffusion to primarily two-dimensional planar diffusion. Such planar diffusion may lead to novel designs of ionic conductors. We show that the dominant vacancy position may vary in the superlattices, depending on the superlattice structure and stacking period, contradicting the common assumption that point defects reside at interfaces. Moreover, we predict a significant increase in room-temperature ionic conductivity for 3STO/3PTO relative to the bulk phases. As a result, considering the variety of cations that can be accommodated in perovskite superlattices and the potential mismatch of spin, charge, and orbitals at the interfaces, this paper identifies a pathway to control defect dynamics for technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion

Zirconium alloys are the most widely used nuclear fuel cladding materials for light water power reactors where irradiation damage causes solute redistribution, leading to degradation of alloy properties such as corrosion resistance. Designing radiation-tolerant zirconium alloys requires a thorough understanding of the atomic-scale transport behavior of the alloying elements in Zr. In this work, we perform density function theory calculations to investigate the diffusion of Sn, Cr, Fe, Be, Al, and Ni in the hexagonal close-packed (HCP) Zr matrix. We develop a methodology to accurately model the metastable vacancy states along the basal migration path, known to occur in group IV metals. We compute the vacancy-mediated solute diffusion coefficients and drag ratios using the kinetic Monte Carlo method and an analytic Green's function method—the agreement between the two validates our methodology. The computed diffusion coefficients of Sn and Al show good agreement with the experimental data and we expect these solutes to diffuse via the vacancy-mediated mechanism. We use a Green's function approach, parameterized with data from density functional theory calculations, to compute the interstitial diffusion coefficients of Cr, Fe, Be, and Ni in the HCP Zr lattice. The computed diffusion coefficients of Cr, Ni, and Be agree with the experimental measurements within one order of magnitude, while those of Fe are within two orders of magnitude of the experimental measurements. The drag ratios for Cr, Fe, Be, and Ni are positive up to 1235 K, which suggests that nonequilibrium vacancy fluxes could drag these solutes toward sinks such as dislocation loops and grain boundaries. We also compute the transport coefficients without including the metastable states, and using the eight- and thirteen-frequency model. Our results show significant differences in drag ratio for the eight- and thirteen-frequency model predictions compared with the Green's function methodology, but smaller errors in the solute diffusivity. Combining interstitial and vacancy-mediated diffusivities, we predict the unusual result that increased vacancy concentration slows down solute diffusivity, while a sufficiently high vacancy concentration can change the dominant mechanism to an accelerated vacancy-mediated diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Transport of hot carriers in plasmonic nanostructures

Plasmonic hot carrier devices extract excited carriers from metal nanostructures before equilibration and have the potential to surpass semiconductor light absorbers. However their efficiencies have so far remained well below theoretical limits, which necessitates quantitative prediction of carrier transport and energy loss in plasmonic structures to identify and overcome bottlenecks in carrier harvesting. Here, we present a theoretical and computational framework, nonequilibrium scattering in space and energy (NESSE), to predict the spatial evolution of carrier energy distributions that combines the best features of phase-space (Boltzmann) and particle-based (Monte Carlo) methods. Within the NESSE framework, we bridge first-principles electronic structure predictions of plasmon decay and carrier collision integrals at the atomic scale, with electromagnetic field simulations at the nano- to mesoscale. Finally, we apply NESSE to predict spatially-resolved energy distributions of photoexcited carriers that impact the surface of experimentally realizable plasmonic nanostructures at length scales ranging from tens to several hundreds of nanometers, enabling first-principles design of hot carrier devices.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Atomistic simulation of defect-dislocation interactions in concentrated solid-solution alloys

The interaction between point defects and dislocations plays a crucial role in governing material properties and microstructural evolutions under external stimuli, such as mechanical deformation and irradiation. Here, we present an atomistic study of the interactions between point defects and dislocations in concentrated solid-solution alloys (CSAs). Using molecular statics and kinetic Monte Carlo methods, we demonstrate that the strain energy and stress field distribution induced by a dislocation in CSAs are highly inhomogeneous along the dislocation line, which leads to heterogeneity of defect-dislocation interactions. Specifically, the interactions are spatially different and screened by the random arrangement of different elemental species. Furthermore, such localization of defect-dislocation interaction indicates that the “dislocation-bias” mechanism that is a driving force for radiation-induced void swelling can be suppressed in concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics

Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.

Bayesian methods↗
DOE OSTI2020

Ab initio study of ${(\nu_\ell,\ell^-)}$ and ${(\overline{\nu}_\ell,\ell^+)}$ inclusive scattering in $^{12}\hbox{C}$: confronting the MiniBooNE and $\hbox{T2K CCQE}$ data

We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections, measured on $^{12}$C by the MiniBooNE and T2K collaborations in the charged-current quasielastic regime. The calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one-and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with $\Lambda_A = 1.0$ or $1.15$ GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗