Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion Monte Carlo”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Correlation effects on coupled electronic and structural properties of doped rare-earth trihydrides

Rare-earth trihydride (𝑅⁢H 3 ) compounds exhibit intriguing coupled electronic and structural properties as a function of doping, hydrogen vacancies, and thermodynamic conditions. Theoretical studies of these materials typically rely on density functional theory (DFT), including the use of small supercells that may underestimate strong correlation effects and structural distortions which in turn may influence their metallicity. Here, we elucidate the roles of lattice distortions and correlation effects on the electronic properties of pristine and doped 𝑅⁢H 3 compounds by adopting DFT +U and quantum Monte Carlo (QMC) methods. Linear-response constrained DFT (LR-cDFT) methods find Hubbard U ≈ 2 eV for 𝑅 𝑑 orbitals and U ≈ 6 eV for H 𝑠 ⁡/N 𝑝 /O 𝑝 orbitals. The small U on Lu 𝑑 orbitals is consistent with QMC calculations on LuH 3 and LuH 2.875 ⁢N 0.125 . In pure face-centered-cubic (FCC) 𝑅⁢H 3 (𝑅 = Lu,Y) compounds, neither DFT nor DFT +U with the self-consistently determined U is enough to create a band gap, however a supercell with hydrogen distortions creates a small gap whose magnitude increases when performing DFT +U with self-consistently determined U values. Furthermore, correlation effects, in turn, have a moderate influence on the coupled structural and electronic properties of doped 𝑅⁢H 3 compounds and may be important when considering the competition between structural distortions and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Advancing earth system model calibration: a diffusion-based method

Understanding of climate impact on ecosystems globally requires site-specific model calibration. Here we introduce a novel diffusion-based uncertainty quantification (DBUQ) method for efficient model calibration. DBUQ is a score-based diffusion model that leverages Monte Carlo simulation to estimate the score function and evaluates a simple neural network to quickly generate samples for approximating parameter posterior distributions. DBUQ is stable, efficient, and can effectively calibrate the model given diverse observations, thereby enabling rapid and site-specific model calibration on a global scale. This capability significantly advances Earth system modeling and our understanding of climate impacts on Earth systems. We demonstrate DBUQ's capability in E3SM land model calibration at the Missouri Ozark AmeriFlux forest site. Both synthetic and real-data applications indicate that DBUQ produces accurate parameter posterior distributions similar to those generated by Markov Chain Monte Carlo sampling but with 30X less computing time. This efficiency marks a significant stride in model calibration, paving the way for more effective and timely climate impact analyses.

Liu, Yanfang↗
DOE OSTI2021

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series

Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗