SEARCH · Search NASA

Results for “Monte Carlo simulations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Deep Direct-Use Feasibility Study Numerical Modeling and Uncertainty Analysis using iTOUGH2 for West Virginia University

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. As a result, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided to use single permeability model for further analysis; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the perdition, therefore, is considered in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis are used for economic analysis. The dataset submitted here support the described study. Manuscript is submitted to Geothermics, will be linked once paper is accepted.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Evaluation of light collection from highly scattering media using wavelength-shifting fibers

Opaque scintillators are designed to have a short scattering length such that scintillation photons are probabilistically confined to a small region of space about their origin. The benefit of this feature is that information on the interaction event topology can be recorded with greater fidelity than in traditional highly transparent media with sensors at large distances from the light production region. Opaque scintillator detectors rely on wavelength-shifting fibers to extract the scintillation light; however, the efficiency of light collection has not yet been directly measured in experiment. We measured the efficiency of light collection as a function of the optical parameters of an opaque liquid and the distance from the origin of the light to the fiber. We use the experimental data to validate a Monte Carlo model of light transport and collection and discuss a simple diffusion model that reproduces the results of Monte Carlo simulation with high fidelity. Furthermore, this combination of validated models has the potential for use in predictions of performance in various designs of future opaque scintillator detectors such as LiquidO.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Uncertainty evaluation for twist drilling stability model

This paper describes the first uncertainty analysis for drilling stability using a frequency-domain drilling stability model. The stability model inputs include: the modal parameters for the torsional-axial vibration mode from the twist drill-holder-spindle axial frequency response function; and the mechanistic coefficients that relate the torque and thrust force to chip area for the selected drill-workpiece material combination. Furthermore, Monte Carlo simulation is applied to propagate the input uncertainties to output uncertainty in the predicted stability map, which separates stable from unstable (chatter) zones in the spindle speed-chip width parameter space. Additionally, the mean stability boundary and its 95% confidence intervals are determined for five cases: varying all four inputs simultaneously and varying them individually. This enables the individual sensitivities to be compared. Experimental results from drilling tests are included for comparison to the prediction. Additionally, Matlab code is provided to implement the stability model and Monte Carlo uncertainty analysis.

42 ENGINEERING↗

DOE OSTI2020

E-Area Low-Level Waste Facility GoldSim System Model

This report documents the development of the E-Area Low-Level Waste Facility (ELLWF) trench system model. The GoldSim® Monte Carlo simulation software is utilized to model the release and transport of radiological inventory disposed (both currently and in the future) within Engineered and Slit Trenches. This model is in support of the sensitivity and uncertainty analysis for the ELLWF Performance Assessment. The ELLWF system model utilizes a hybrid-approach to accurately describe the disposal system. The Hydrologic Evaluation of Landfill Performance model provides the infiltration data to both PORFLOW and GoldSim. PORFLOW is used to calibrate the GoldSim model to ensure confidence in the stochastic results. Finally, the concentrations from GoldSim transport simulations are fed into the SRNL Dose Toolkit to calculate dose impacts and assess plume interaction

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗

DOE OSTI2023

MOFX-DB: An Online Database of Computational Adsorption Data for Nanoporous Materials

Machine learning and data mining coupled with molecular modeling have become powerful tools for materials discovery. Metal-organic frameworks (MOFs) are a rich area for this due to their modular construction and numerous applications. Here, we make data from several previous large-scale studies in MOFs and zeolites from our groups (and new data for N 2 and Ar adsorption in MOFs) easily accessible in one place. The database includes over 3 million simulated adsorption data points for H 2 , CH 4 , CO 2 , Xe, Kr, Ar, and N 2 in over 160 000 MOFs and zeolites, textural properties like pore sizes and surface areas, and the structure file for each material. We include metadata about the Monte Carlo simulations to enable reproducibility. The database is searchable by MOF properties, and the data are stored in a standardized JSON format that that is interoperable with the NIST adsorption database. We also identify several MOFs that meet high performance targets for multiple applications, such as high storage capacity for both hydrogen and methane or high CO 2 capacity plus good Xe/Kr selectivity. Here, by providing this data publicly, we hope to facilitate machine learning studies on these materials, leading to new insights on adsorption in MOFs and zeolites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Modulation of intensity emerging from zero effort (MIEZE) with extended Fourier time at large scattering angle

Modulation of Intensity Emerging from Zero Effort (MIEZE) is a neutron resonant spin echo technique that allows one to measure time correlation scattering functions in materials by implementing radio-frequency (RF) intensity modulation at the sample and the detector. The technique avoids neutron spin manipulation between the sample and the detector and, thus, could find applications in cases where the sample depolarizes the neutron beam. However, the finite sample size creates a variance in the path length between the locations where scattering and detection happen, which limits the contrast in intensity modulation that one can detect, in particular, toward long correlation times or large scattering angles. We propose a modification to the MIEZE setup that will enable one to extend those detection limits to longer times and larger angles. Furthermore, we use Monte Carlo simulations of a neutron scattering beamline to show that by tilting the RF flippers in the primary spectrometer with respect to the beam direction, one can shape the wave front of the intensity modulation at the sample to compensate for the path variance from the sample and the detector. The simulation results indicate that this change enables one to operate a MIEZE instrument at much increased RF frequencies, thus improving the effective energy resolution of the technique. For the MIEZE instrument simulated, it shows that for an incident beam with the maximum divergence of 0.33°, the maximum Fourier time can be increased by a factor of 3.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Sputter-deposited low-stress boron carbide films

Significant challenges have been faced in the manufacturing of low-stress B 4 C films. Here, we demonstrate a set of process parameters for direct-current magnetron sputter deposition of smooth, high-purity, amorphous B 4 C films with near-zero total residual stress and with thicknesses up to 10 μm. Films are characterized by a combination of high-energy ion scattering, x-ray diffraction, electron microscopy, and stress measurements based on substrate curvature. In order to facilitate the process transfer to other sputtering geometries, the favorable deposition conditions are correlated with distributions of landing energies and incident angles of depositing species estimated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. Based on such simulations, a decrease in compressive stress with increasing Ar working gas pressure is attributed to the corresponding broadening of the angular distribution of depositing species and associated atomic shadowing effects.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

Here, the target segments of the Oak Ridge National Laboratory Second Target Station (STS) neutron production facility become highly activated due to spallation reactions or nuclei transmutation by primary protons and emitted neutrons. Once the target segments are removed from their location within the core vessel, decay dose rates must be accurately quantified to determine the shielding configurations of remote-handling tools and transport casks and to aid in planning maintenance activities. For this analysis, we utilized a hybrid unstructured mesh (UM)/constructive solid geometry approach for calculating spallation products and neutron fluxes, activation calculations using the AARE package that includes the CINDER2008 activation code to calculate the decay photon source at different cooling times, and the ADVANTG code to accelerate the final decay photon transport calculation. Both Type 316 stainless steel (SS-316) and lead were investigated as candidates for shielding materials. The decay photon transport calculation through the thick SS-316 or lead shields exhibited between 25 and 30 orders-of-magnitude attenuations in the radial direction, depending on the shield. Such a difficult shielding calculation required advanced variance reduction. ADVANTG has some missing features, which limits its usability in spallation neutron source applications. It does not support volumetric sources created for MCNP6.2 UM capability. An approximate source was created for this problem. Not only was this approximate source needed for running the ADVANTG calculation to generate the weight windows, but also it was essential to develop source biasing (SB) parameters that were crucial for dramatically accelerating the decay photon transport in this problem. With this approximate source, the analysis was completed in a very reasonable computational time, and the design of the STS remote-handling equipment was finalized. This paper compares the efficiency of Monte Carlo simulations with different weight window and SB parameters calculated using different approximate ADVANTG calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Zero-point motion of liquid and solid hydrogen

We present an inelastic neutron scattering study of liquid and solid hydrogen carried out using the wide Angular Range Chopper Spectrometer at Oak Ridge National Laboratory. From the observed dynamic structure factor, we obtained empirical estimates of the molecular mean-squared displacement and average translational kinetic energy. Here we find that the former quantity increases with temperature, indicating that a combination of thermal and quantum effects is important near the liquid-solid phase transition, contrary to previous measurements. We also find that the kinetic energy drops dramatically on melting of the crystals, a consequence of the large increase in molar volume together with the Heisenberg indeterminacy principle. Our results are compared with quantum Monte Carlo simulations based on different model potentials. In general, there is good agreement between our findings and theoretical predictions based on the Silvera-Goldman and Buck potentials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

First-principles analysis of the Al-rich corner of Al-Li-Cu phase diagram

The phase diagram of Al-Li-Cu system in the Al-rich region was determined by means of first-principles calculations and statistical mechanics. The mixing enthalpies of many configurations for different lattices in the whole Al-Li-Cu system were determined by density functional theory simulations to find the stable phases in the convex hull. They were fitted with a cluster expansion to calculate the free energy of the configurations with different compositions as a function of temperature in the Al-rich region (Al content > 40 at.%) by means of Monte Carlo simulations. It was found that the ground state phases in the Al-rich part of the Al-Li-Cu phase diagram were α-Al, θ' (Al 2 Cu), $δ$' (Al 3 Li), $δ$ (AlLi) and T 1 (Al 6 Cu 4 Li 3 ), while θ'' (Al 3 Cu), T 1' (Al 2 CuLi) and Al 3 Cu 2 Li were found on the lowest mixing enthalpy surfaces of their lattices and were metastable. α-Al, $δ$ and T 1 are stable phases in the whole temperature range while $δ$' becomes metastable at very low temperature and θ (Al 2 Cu) replaces θ' as the stable phase at approximately 550 K due to the vibrational entropic contribution. In addition, the phase diagram in the Al-rich region was built and it was shown in isothermal sections from 100 to 900 K. They were in good agreement with the limited experimental data in the literature and provided new information regarding the stability, solubility and stoichiometry of the different phases. This information is important to understand the precipitation mechanisms during high temperature aging.

36 MATERIALS SCIENCE↗

DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗

DOE OSTI2022

The throughput calibration of the VERITAS telescopes

Context. The response of imaging atmospheric Cherenkov telescopes to incident γ -ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ -ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time. Methods. The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy. Results. Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired γ -ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗