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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Determination of the spins and parities for the $0^+_4$ and $0^+_5$ states in $^{100}Zr$

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 (MCSM) 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 OSTI2024

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

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 (MCSM) 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 OSTI2022

First Measurement of Differential Charged Current Quasielastic-like $\nu_\mu$–Argon Scattering Cross Sections in Kinematic Imbalance Variables with the MicroBooNE Detector

We report the first measurement of flux-integrated multi-differential cross sections for charged-current events with muon neutrinos scattering on argon with solely a muon and a single proton in the final state as a function of kinematic imbalance variables. The measurement was carried out using the Booster Neutrino Beam at Fermi National Accelerator Laboratory within the MicroBooNE Liquid Argon Time Projection Chamber detector with an exposure of 6.79 × 10 20 protons on target. Events were selected to enhance the contribution of charged-current mesonless interactions with one proton detected in the final state. The data discussed here are reported in terms of multi-differential cross sections in kinematic imbalance variables, which are generally sensitive to nuclear effects. The double-differential results in these variables can provide an excellent handle to disentagle specific nuclear aspects not easily isolated via single differential cross sections. Our results pave a path towards identifying regions of the phase-space where future interaction modeling development and Monte Carlo neutrino generator tuning efforts should concentrate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

STAT7 v1.2 User Guide: 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, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. 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 Massachusetts Institute of Technology Research Reactor. 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 uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for thermal hydraulics analyses of research and test reactors. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to the fit functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one latera stripe of a fuel element. 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 fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Simulation Study of Monoenergetic Photon Source LiDAR System Performance for Emergency Response Applications

The following report summarizes the results of a Monte Carlo simulation and modeling study of the Lawrence Berkeley National Laboratory’s (LBL) Monoenergetic Photon Source (MPS) as part of a high-energy LiDAR system for use in emergency response applications. A set of image quality indicators (IQI) and nonsensitive phantoms or generalized surrogates for search targets have been identified for use in the planned measurement campaign at the Berkeley Lab Laser Accelerator (BELLA) Center toward the end FY2023. These items were modeled and the time-of-flight (TOF) spectral responses following irradiation were calculated for a set of varying operating conditions. The goal of this study was to both predict the LiDAR system performance as well as determine if minor modifications to the IQI and phantom designs were required before fabrication. It was found that the IQI and target designs both overlap and extend towards the limit of performance of the system and are thus suitable for the final project measurement campaign.

43 - PARTICLE ACCELERATORS↗