Engineering Papers⌕ Search

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.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Understanding structure-processing relationships in metal additive manufacturing via featurization of microstructural images

Understanding and predicting accurate property-structure-processing relationships for additively manufactured components is important for both forward and inverse design of robust, reliable parts and assemblies. While direct mapping of process parameters to properties is sometimes plausible, it is often rendered difficult due to poor microstructural control. Exploring the direct relationship between processing conditions and microstructural features can thus provide significant physical insights and aid the overall design process. Here, in this study, we develop an automated high-throughput framework to simulate an uncertainty-aware additive manufacturing (AM) process, characterize microstructural images, and extract meaningful features/descriptors. A kinetic Monte Carlo (KMC) based model of the AM process is used to simulate microstructural evolution for a diverse set of experimentally relevant processing conditions. We perform a parametric study to explore the relationship between microstructural features and processing conditions. Our results indicate that a many-to-one mapping can exist between processing conditions and typical descriptors; therefore, multiple descriptors are thus necessary to unambiguously represent microstructural images. Our work provides crucial quantitative and qualitative in-formation that would aid in the selection of features for microstructural images. Featurized microstructures could then be utilized to build data-driven models for predictive control of microstructures and thereby properties of additively manufactured components.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Advanced manufacturing of 3D custom boron-carbide collimators designed for complex environments for neutron scattering

Scattered-beam collimation is a very useful method to reduce unwanted backgrounds and to boost the desired sample signal instead. This approach is of particular interest for samples contained within a complex environment that gives rise to much unwanted parasitic scatter. As neutron scattering instrument and techniques advances, small samples are becoming of more and more interest, which necessitates optimized collimation. Here, in this work, we describe a concept for the design and fabrication of advanced scattered-beam collimation 3D printed from B 4 C specifically tailored for samples contained within a complex environment. This concept is demonstrated through the use of a diamond anvil cell for high pressure experimentation, a technique that very typically requires small samples. The collimators here are designed through a modeling procedure via Monte Carlo neutron ray tracing that encompasses the entire system: the instrument, the complex environment and the collimator. Since the first approach of simply scaling up of the print-size was not successful, a novel concept of a multi-part alternate-blade collimator was developed. This approach addresses printing constraints but gives greater flexibility in design. Its performance is computationally compared against an unprintable progressively tighter blade collimator to assess the effect of alternating blades. No strong difference was observed. Its performance was validated through experimentation at the Spallation Neutron Source. The results emphasize the critical importance of ultra-high precision alignment while showing good overall agreement between simulation and experiment and underscore the feasibility of the method and its real-world application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

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

Non-Markovian Hole Excess Noise in Avalanche Amorphous Selenium Thin Films

Enhancing the signal-to-noise ratio in avalanche photodiodes by utilizing impact ionization gain requires materials exhibiting low excess noise factors. Amorphous selenium (a-Se) as a wide bandgap at ~2.1 eV, a solid-state avalanche layer, demonstrates single-carrier hole impact ionization gain and manifests ultralow thermal generation rates. A comprehensive study of the history dependent and non-Markovian nature of hot hole transport in a-Se was modeled using a Monte Carlo (MC) random walk of single hole free flights, interrupted by instantaneous phonon, disorder, hole–dipole, and impact-ionization scattering interactions. The hole excess noise factors were simulated for 0.1–15 μm a-Se thin-films as a function of mean avalanche gain. The hole excess noise factors in a-Se decreases with an increase in electric field, impact ionization gain, and device thickness. The history dependent nature of branching of holes is explained using a Gaussian avalanche threshold distance distribution and the dead space distance, which increases determinism in the stochastic impact ionization process. An ultralow non-Markovian excess noise factor of ~1 was simulated for 100 nm a-Se thin films corresponding to avalanche gains of 1000. Future detector designs can utilize the nonlocal/non-Markovian nature of the hole avalanche in a-Se, to enable a true solid-state photomultiplier with noiseless gain.

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

FORTE Measurements of Global Optical Lightning Waveforms and Implications for Optical Lightning Detection

Lightning processes generate a diverse collection of optical pulses from current traversing the lightning channels. These signals are then broadened spatially and temporally via scattering in the clouds. The resulting waveforms measured from space with instruments like the photodiode detector (PDD) on the Fast On-orbit Recording of Transient Events (FORTE) satellite have a variety of shapes. In this study, we use coincident optical and Radio-Frequency measurements to document the properties of optical PDD waveforms associated with different types of lightning, estimate delays from scattering in the clouds, and comment on how pulse shape impacts optical lightning detection. We find that the attributes of optical pulses recorded by the PDD are consistent with prior studies, but vary globally and with event amplitude. The most powerful lightning tends to be single-peaked with faster rise times (median: ~100 µs) and shorter effective widths (median: ~400 µs) than normal lightning. Particularly dim events, meanwhile, include cases of broad optical waveforms with sustained optical emission throughout the PDD record, which the pixelated FORTE Lightning Location System (LLS) instrument has difficulty detecting. We propose that this is due to the optical signal being divided between individual LLS pixels that are each, individually, not bright enough to trigger. We use PDD waveforms and Monte Carlo radiative transfer modeling to demonstrate that increasing the temporal and spatial resolution of a pixelated lightning imager will make it more difficult to detect these broad/dim pulses as their energy becomes divided between additional pixels/integration frames.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermonuclear neutron emission from a sheared-flow stabilized Z-pinch

We report the fusion Z-pinch experiment (FuZE) is a sheared-flow stabilized Z-pinch designed to study the effects of flow stabilization on deuterium plasmas with densities and temperatures high enough to drive nuclear fusion reactions. Results from FuZE show high pinch currents and neutron emission durations thousands of times longer than instability growth times. While these results are consistent with thermonuclear neutron emission, energetically resolved neutron measurements are a stronger constraint on the origin of the fusion production. This stems from the strong anisotropy in energy created in beam-target fusion, compared to the relatively isotropic emission in thermonuclear fusion. In dense Z-pinch plasmas, a potential and undesirable cause of beam-target fusion reactions is the presence of fast-growing, “sausage” instabilities. This work introduces a new method for characterizing beam instabilities by recording individual neutron interactions in plastic scintillator detectors positioned at two different angles around the device chamber. Histograms of the pulse-integral spectra from the two locations are compared using detailed Monte Carlo simulations. These models infer the deuteron beam energy based on differences in the measured neutron spectra at the two angles, thereby discriminating beam-target from thermonuclear production. An analysis of neutron emission profiles from FuZE precludes the presence of deuteron beams with energies greater than 4.65 keV with a statistical uncertainty of 4.15 keV and a systematic uncertainty of 0.53 keV. This analysis demonstrates that axial, beam-target fusion reactions are not the dominant source of neutron emission from FuZE. These data are promising for scaling FuZE up to fusion reactor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fixed Source Sensitivity Calculations for Inertial Confinement Fusion Applications

A numerical code library was developed for the radiation transport code MCNP6.3 to calculate generalized response sensitivity coefficients for fixed source neutron transport problems with applications to inertial confinement fusion (ICF) experiments. The new MCNP6.3 dependency is used to generate a novel time convolution response that represents a neutron time-of-flight (nToF) signal. The traditional suite of macroscopic cross-section sensitivities and constrained fixed source probability distribution sensitivities are available for both the standard and the new response tallies in this library. However, novel sensitivity coefficients for the constrained hyperparameters of analytic fixed source probability distributions are emphasized in this work for their connection to ICF neutron transport models. Particularly, advanced Monte Carlo methods are developed for calculating the sensitivity of a nToF signal to perturbations in an ICF plasma’s ion temperature and burn history as well as perturbations in the target liner mass density and the shape parameters of the nToF detector’s impulse response function. Together, these capabilities form an advanced suite of computational tools that can be used to analyze and extract information from any ICF experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solar reflection of dark matter with dark-photon mediators

Abstract We consider the scattering of low-mass halo dark-matter particles in the hot plasma of the Sun, focusing on dark matter that interact with ordinary matter through a dark-photon mediator. The resulting “solar-reflected” dark matter (SRDM) component contains high-velocity particles, which significantly extend the sensitivity of terrestrial direct-detection experiments to sub-MeV dark-matter masses. We use a detailed Monte Carlo simulation to model the propagation and scattering of dark-matter particles in the Sun, including thermal effects, with special emphasis on ultralight dark-photon mediators. We study the properties of the SRDM flux, obtain exclusion limits from various direct-detection experiments, and provide projections for future experiments, focusing especially on those with silicon and xenon targets. We find that proposed future experiments with xenon and silicon targets can probe the entire “freeze-in benchmark”, in which dark matter is coupled to an ultralight dark photon, including dark-matter masses as low as 𝒪(keV). Our simulations and SRDM fluxes are publicly available.

Astronomy & Astrophysics↗
DOE OSTI2021

Comparative study of He 6 β -decay based on different similarity-renormalization-group evolved chiral interactions

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multistep Coulomb excitation of Ni 64 : Shape coexistence and nature of low-spin excitations

Here, the structure of 64 Ni, the heaviest stable Ni isotope, has been investigated via high-statistics, multistep safe Coulomb excitation to search for shape coexistence, a phenomenon recently observed in neutron-rich 66 Ni and 70 Ni as well as in doubly magic, N = 40, 68 Ni. The study was motivated by recent, state-of-the-art Monte Carlo shell-model calculations (MCSM), where a Hamiltonian with effective interactions incorporating the monopole tensor force predicts the existence of shape coexistence, also in the lower-mass 62,64 Ni isotopes. A set of transition and static E2 matrix elements for both yrast and near-yrast structures was extracted from the differential Coulomb excitation cross sections. From comparisons between the new results and MCSM as well as other shell-model calculations, a clearer picture of the structure of 64 Ni emerges. Specifically, the low-spin states are shown to be dominated by proton and neutron excitations mainly within the fp shell, with minimal contribution from the g 9/2 shape-driving neutron orbital. The agreement between experimental data and MCSM results indicates a small oblate deformation for the $0$$^{+}_{2}$ level and a spherical shape for the $0$$^{+}_{3}$ state. In addition, the small upper limit determined for the B(E2) probability of a transition associated with the decay of the recently observed 3463-keV, $0$$^{+}_{4}$ state agrees with its proposed assignment to a prolate shape, herewith providing first evidence for triple shape coexistence in a stable Ni isotope.

59 ≤ A ≤ 89↗
DOE OSTI2023

First spectroscopic study of odd-odd Cu 78

Nuclei in the vicinity of 78 Ni are important benchmarks for nuclear structure, which can reveal changes in the shell structure far from stability. Spectroscopy of the odd-odd isotope 78 CU was performed for the first time in an experiment with the EURICA setup at the Radioactive Isotope Beam Factory at RIKEN Nishina Center. Excited states in the neutron-rich isotope were populated following the β decay of 78 Ni produced by in-flight fission and separated by the BigRIPS separator. A level scheme based on the analysis of $\gamma$-$\gamma$ coincidences is presented. Tentative spin and parity assignments were made when possible based on the β -decay feeding intensities and $\gamma$-decay properties of the excited states. Time correlations between β and $\gamma$ decay show clear indications of an isomeric state with a half-life of 3.8(4) ms. Large-scale Monte Carlo shell-model calculations were performed using the A3DA-m interaction and a valence space comprising the full f p shell and the 1⁢g 9/2 and 2⁢d 5/2 orbitals for both protons and neutrons. The comparison of the experimental results with the shell-model calculations allows interpreting the excited states in terms of spin multiplets arising from the proton-neutron interaction. In conclusion, the results provide further insight into the evolution of the proton single-particle orbitals as a function of neutron number, and quantitative information about the proton-neutron interaction outside the doubly magic 78 Ni core.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Determination of the spins and parities for the 0$^{+}_{4}$ and 0$^{+}_{5}$ states in 100 Zr

Here, two 0 + states at 1294.5 and 1774.0 keV, together with three 2 + and one 4 + levels, were identified or unambiguously spin-parity assigned for the first time in 100 Zr utilizing γ-ray spectroscopy and γ–γ angular correlation techniques with the Gammasphere spectrometer, following the β – decay of neutron-rich, mass-separated 100,100m Y isotopes. Comparisons with recent Monte Carlo shell-model calculations indicate that these two states are candidates for the bandhead of a sequence in a shape-coexisting spherical minimum predicted to be located around ≈1500 keV. According to the measured relative B(E2) relative transition probabilities, the 0$^+_5$ state exhibits decay properties which more closely align with those predicted for a spherical shape, while the 0$^+_4$ level is suggested to be associated with a weakly deformed shape similar to one related to the 0$^+_2$ state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Resilience of the surface code to error bursts

Quantum error correction works effectively only if the error rate of gate operations is sufficiently low. However, some rare physical mechanisms can cause a temporary increase in the error rate that affects many qubits; examples include ionizing radiation in superconducting hardware and large deviations in the global control of atomic systems. We refer to such rare transient spikes in the gate error rate as error bursts. In this work, we investigate the resilience of the rotated surface code to generic error bursts. We assume that, after appropriate mitigation strategies, the spike in the error rate lasts for only a single syndrome-extraction cycle; we also assume that the enhanced error rate is uniform across the code block. Under these assumptions, and for a circuit-level depolarizing noise model, we perform Monte Carlo simulations to determine the regime in burst error rate and background error rate for which the memory time becomes arbitrarily long as the code block size grows. Our results indicate that suitable hardware mitigation methods combined with standard decoding methods may suffice to protect against transient error bursts in the rotated surface code.

Quantum error correction↗
DOE OSTI2025

Estimation of backgrounds from jets misidentified as τ-leptons using the Universal Fake Factor method with the ATLAS detector

Processes with τ$$\tau $$-leptons in the final state are important for Standard Model measurements and searches for physics beyond the Standard Model. The ATLAS experiment at the Large Hadron Collider observes τ$$\tau $$-leptons produced in proton–proton collisions only through their decay products. Data analyses involving hadronically decaying τ$$\tau $$-leptons face challenges due to backgrounds from jets misidentified as τ$$\tau $$-leptons that are not modelled reliably by Monte Carlo simulations. Data-driven methods such as the fake-factor method allow such misidentified backgrounds to be predicted by measuring transfer factors, known as fake factors, in data from dedicated regions. This paper describes a refined technique for determining the fake factors, the Universal Fake Factor method. It evaluates the fake factors for a signal region by using fake factors from samples enriched in different sources of jets misidentified as τ$$\tau $$-leptons (light-quark, gluon, b-quark, and pile-up jets). Each fake factor is calculated as a linear combination of fake factors measured in these different enriched samples. For the full Run 2 data set, the systematic uncertainty of the calculated fake factors, evaluated using W(μν)$$W(\mu u )$$ enriched event sample, ranges from 15 to 35% depending on the τ$$\tau $$-lepton’s transverse momentum and charged-particle decay multiplicity.

Aad, G↗
DOE OSTI2026

Thermodynamic Modeling of Intrinsic Defects in MnBi₂Te₄

This repository contains the computational data supporting the manuscript titled “The critical role of intrinsic defects and many-body interactions on the stability of MnBi₂Te₄.” It includes: 1. DFT data generated using VASP, used for training and benchmarking electronic structure models. 2. Quantum Monte Carlo (QMC) data produced with QMCPACK, used to apply many-body corrections and validate the electronic and magnetic properties of MnBi₂Te₄. 3. Relevant scripts used to run, analyze, and process the calculations, enabling reproducibility and transparency of the workflows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Benchmark Evaluation of one Dimensional Array of HEU Moderated and Reflected by Lucite [Poster]

Two critical experiments with HEU-Lucite were performed using the Planet Universal Critical Assembly Machine at Los Alamos National Laboratory in 2019. HEU foils were interleaved with Lucite in a column stack for moderation and reflection from the square Lucite plates. The evaluation of the experiments is presented. A Monte-Carlo (MCNP6) calculation model was developed for the determination of uncertainties and their effect on the multiplication factor. Experimental uncertainties were found to be low with the overall uncertainty in k eff approximately 240 pcm. The comparison of calculations with experimental data is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The STAT7 Code for Statistical Propagation of Uncertainties in Steady State Thermal Hydraulics Analysis of Plate Fueled Reactors

The STAT7 software was developed to perform steady-state thermal hydraulic analyses. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as MITR-II. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model common fabrication parameters and other key reactor analysis uncertainties required for research and test reactor thermal hydraulic analyses. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one lateral stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each element in the core.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗