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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Molecular coupling competing with defects within insulator of the magnetic tunnel junction-based molecular spintronics devices

Nearly 70 years old dream of incorporating molecule as the device element is still challenged by competing defects in almost every experimentally tested molecular device approach. This paper focuses on the magnetic tunnel junction (MTJ) based molecular spintronics device (MTJMSD) method. An MTJMSD utilizes a tunnel barrier to ensure a robust and mass-producible physical gap between two ferromagnetic electrodes. MTJMSD approach may benefit from MTJ's industrial practices; however, the MTJMSD approach still needs to overcome additional challenges arising from the inclusion of magnetic molecules in conjunction with competing defects. Molecular device channels are covalently bonded between two ferromagnets across the insulating barrier. An insulating barrier may possess a variety of potential defects arising during the fabrication or operational phase. This paper describes an experimental and theoretical study of molecular coupling between ferromagnets in the presence of the competing coupling via an insulating tunnel barrier. We discuss the experimental observations of hillocks and pinhole-type defects producing inter-layer coupling that compete with molecular device elements. We performed theoretical simulations to encompass a wide range of competition between molecules and defects. Monte Carlo Simulation (MCS) was used for investigating the defect-induced inter-layer coupling on MTJMSD. Our research may help understand and design molecular spintronics devices utilizing various insulating spacers such as aluminum oxide (AlOx) and magnesium oxide (MgO) on a wide range of metal electrodes. This paper intends to provide practical insights for researchers intending to investigate the molecular device properties via the MTJMSD approach and do not have a background in magnetic tunnel junction fabrication.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Kinetic simulations of ignited mode cesium vapor thermionic converters

Cesium vapor thermionic converters are an attractive method of converting high-temperature heat directly to electricity, but theoretical descriptions of the systems have been difficult due to the multi-step ionization of Cs through inelastic electron–neutral collisions. This work presents particle-in-cell simulations of these converters, using a direct simulation Monte Carlo collision model to track 52 excited states of Cs. Here, these simulations show the dominant role of multi-step ionization, which also varies significantly based on both the applied voltage bias and pressure. The electron energy distribution functions are shown to be highly non-Maxwellian in the cases analyzed here. A comparison with previous approaches is presented, and large differences are found in ionization rates due especially to the fact that previous approaches have assumed Maxwellian electron distributions. Finally, an open question regarding the nature of the plasma sheaths in the obstructed regime is discussed. The one-dimensional simulations did not produce stable obstructed regime operation and thereby do not support the double-sheath hypothesis.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Calibration and simulation of ionization signal and electronics noise in the ICARUS liquid argon time projection chamber

The ICARUS liquid argon time projection chamber (LArTPC) neutrino detector has been taking physics data since 2022 as part of the Short-Baseline Neutrino (SBN) Program. This paper details the equalization of the response to charge in the ICARUS time projection chamber (TPC), as well as data-driven tuning of the simulation of ionization charge signals and electronics noise. The equalization procedure removes non-uniformities in the ICARUS TPC response to charge in space and time. This work leverages the copious number of cosmic ray muons available to ICARUS at the surface. The ionization signal shape simulation applies a novel procedure that tunes the simulation to match what is measured in data. The end result of the equalization procedure and simulation tuning allows for a comparison of charge measurements in ICARUS between Monte Carlo simulation and data, showing good performance with minimal residual bias between the two.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Short-range magnetic correlations in quasicrystalline i -Tb-Cd

Here, we report on elastic and inelastic neutron scattering from single-grain isotopically enriched samples to elucidate the local magnetic correlations between Tb 3+ moments in quasicrystalline i -Tb-Cd. The inelastic neutron scattering measurements of the crystalline electric field excitations demonstrated that the Tb 3+ moments are directed primarily along the local fivefold axes of the Tsai-type cluster as was found for the TbCd 6 approximant phase. Based on the inelastic measurements, we consider a simple Ising-type model for the moment configurations on a single Tb 3+ icosahedron and enumerate the lowest energy moment configurations. We then calculate the diffuse scattering from these configurations and compare with the experimental magnetic diffuse scattering measurements to identify the most likely single cluster moment configurations and find reasonable agreement between the broad features observed in our scattering simulations. We further use a heuristic glass model as well as large-scale Monte-Carlo simulations of a multicluster model to consider the role of higher-order (longer-range) intercluster correlations for magnetic frustration and the magnetic scattering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Chirality-induced zigzag domain wall in in-plane magnetized ultrathin films

In this report the domain structure in in-plane magnetized Fe/Ni/W(110) films is investigated using spin-polarized low-energy electron microscopy. A novel transition of the domain wall shape from a zigzaglike pattern to straight is observed as a function of the film thickness, which is triggered by the transition of the domain wall type from the out-of-plane chiral wall to the in-plane Néel wall. The contribution of the Dzyaloshinskii-Moriya interaction to the wall energy is proposed to explain the transition of the domain wall shape, which is supported by Monte Carlo simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Semiclassical treatment of photon cascades in nuclei

Here, we present a simple semiclassical treatment of photon cascades, suitable for use in nuclear fission simulation codes. The approximation is here developed for E⁢1 and E⁢2 transitions and its quality is illustrated for a variety of two-photon cascades. Implementation of the treatment into Monte Carlo simulations would make it possible to address photon correlation observables quantitatively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Monte Carlo N-Particle Transport Performance of Predicting Digital Radiographic IQI Inspection

The identification of porosity, geometric noncompliance, and other defect types are critical to the qualification of materials and components. X-ray radiographic nondestructive testing is a common industrial inspection method for process quality control and component qualification and certification. Digital radiography provides a quick and efficient alternative when compared to traditional film-based inspection. The quality of radiographic inspection is dependent on equipment specifications, such as the source spot size and detector pixel size, and the specific parameters selected for use for the radiographic technique. To evaluate if an x-ray system and technique is sufficient for a given requirement, a radiographic image quality indicator (IQI) can be used. Radiographic IQIs in hard to machine materials or hard to manufacture defects can be time consuming and expensive to manufacture. This study was conducted to evaluate current Savannah River National Laboratory (SRNL) x-ray imaging systems with a custom tantalum IQI and using Monte Carlo simulations to predict the performance of future systems. The tantalum IQI was tested using a Siefert Isovolt 420 keV x-ray tube with a Perkin Elmer XRD 1611 flat panel with 100-micron pixels. Using the Monte Carlo N-Particle transport software, the radiographic tally was used to simulate the photon flux through an identical tantalum IQI. These simulations provided a benchmark as to the best theoretical identification on a given system using our tantalum IQI. The simulations were refined to match SRNL’s current systems’ noise levels, leading to confidence in their ability to predict the performance of other systems that may be purchased and deployed in the future at the Savannah River Site. Future studies will be conducted to prove this research can be extended to artificially evaluate the ability for systems to identify critical defect sizes through x-ray radiographic inspection, drastically reducing the cost and time burdens of producing high-fidelity radiographic test articles.

digital X-ray radiography↗

DOE OSTI2024

Wavefunction matching for solving quantum many-body problems

Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions to quantum chemistry and from atomic and molecular systems to nuclear physics. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

A unifying equation for fermentation sustainability across the titer-rate-yield landscape

Industrial fermentation is central to the sustainable production of fuels and chemicals, yet commercial viability of emerging technologies hinges on improving fermentation titer, rate, and yield (TRY). How these metrics shape system cost remains difficult to generalize due to complex interactions among feedstocks, fermentation, separations, catalytic upgrading, waste management, and facility design. Here, we systematically map theoretical fermentation performance spaces (formed by all potential TRY combinations) for 32 representative biomanufacturing facilities—spanning distinct choices for feedstocks, fermentation regimes and products, separations, and catalytic upgrading—by simulating and evaluating them (via techno-economic analysis, TEA) under uncertainty (600,000 Monte Carlo simulations) and across TRY combinations (7500 TRY combinations for each of 32 configurations). Across this wide design and thermodynamic simulation space, we find the relationship between fermentation TRY and system cost is captured by a simple, generalizable mathematical equation (R 2 of 0.992 − 1.000 across our simulations; 0.954 − 1.000 when validated against prior studies that used different tools). We use this equation to elucidate key drivers that shape cost sensitivity to fermentation performance, generating widely applicable insights. By demonstrating a unifying relationship governs the impact of fermentation on biomanufacturing economics, this work establishes a foundation for agile, holistically predictive, resource-efficient strategies to prioritize fermentation research and development needs and accelerate commercialization of emerging biomanufacturing technologies.

applied mathematics↗

DOE OSTI2021

Cherenkov and Transition Radiation as Low-Energy Background Sources in SuperCDMS Detectors

The Weakly Interacting Massive Particle (WIMP) has historically been a prime candidate for dark matter due to its elegant compatibility with the Minimally Supersymmetric Standard Model (MSSM). Recent dark matter experiments have ruled out much of the GeV-TeV mass range predicted by the MSSM, however, and the newest generation of direct detection experiments, such as SuperCDMS SNOLAB, have begunto explore sub-GeV dark matter. As experiments push towards lower masses, sensitivity to eV-scale electron recoil events has become increasingly important. A variety of unexplained excesses at energy deposits of 1-100 eV have been found inmany such low-mass experiments, across different detection techniques and at different excess rates. These low-energy events are not thought to be dark matter, and they must be understood and effectively removed to further improve detector energy resolution. This thesis will outline the simulation of one possible source of low-energy excess at SuperCDMS: optical photons produced from charges in uniform motion, specifically from Cherenkov radiation (CR), transition radiation (TR), and the intermediate hybrid transition-Cherenkov radiation (HR). In the latter case, the standard formulae for HR contain divergences at certain angles in transparent media, which are an artifact of their derivation and become problematic to simulate. To rsolve this, we present a novel approach to normalize the divergent HR peaks, which allows HR to transition smoothly between TR and CR in a numerical simulation. We also add TR and CR as a physics process to the SuperCDMS Monte Carlo simulation package SuperSim, based on the Geant4 simulation toolkit. To verify the physics of our addition, we show some test simulations in comparison to theoretical predictions of optical radiation intensity. We also use the simulation to make some preliminary predictions on the contribution of TR and CR to the low energy background, with the expectation that more thorough analyses will be conducted in the future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

The Impact of Behind-the-Meter Heterogeneous Distributed Energy Resources on Distribution Grids: Preprint

The increasing integration of distributed energy resources (DERs) on the electric grid brings new challenges and opportunities for utility grid operations. With the rapid deployment of DERs, there is emerging interest in integrating these controllable devices with utility operations at all levels for monitoring and management. Different types of DERs such as Photovoltaics (PV), energy storage, electric vehicles (EVs), etc. have varying effects on the grid based on how and where they are deployed. Again, the DER deployment and its impact are network-dependent, while the traditional electric grids were not designed to host DERs. To understand the challenges with increasing behind-the-meter (BTM) DERs and to identify the needs in deploying advanced controls, a comprehensive grid impact study is indispensable. This paper investigates the impacts of integrating a mix of DERs in the utility distribution network in Colorado, U.S. Firstly, BTM DERs at the residential scale including distributed PV; battery energy storage systems (BESS); heating, ventilation, and air-conditioning (HVAC) load; electric water heater (EWH) load; and EVs are modeled. The impacts of integrating these resources on the distribution network are evaluated by conducting time-series simulations for different scenarios considering different days to capture the worst-case conditions. Monte Carlo simulations are conducted to generate the realistic EV charging profile. The voltage issues, substation transformer loadings, and critical nodes in the network are identified with the incorporation of uncoordinated EV charging loads and other DER in the network. This analysis helps visualize the DER impact on the grid, identify the challenges, and provides an insight into the new distribution management and control needs to enable reliable and resilient distribution grid operations. Additionally, further analysis is performed to estimate the daily residential electricity cost with the inclusion of DERs under time-of-use tariffs. The result shows the daily residential electricity cost reduced by 20:7% on average with DERs compared to the case without DERs.

27 ARPA - Advanced Research Projects Agency-Energy↗

DOE OSTI2022

Uncertainty Quantification of Bifacial Performance Modeling

Analysis on uncertainty in the annual energy of PV systems that can be attributed to parameters of particular importance to bifacial PV modules is presented. Monte Carlo simulations are used to evaluate the effect of uncertain module bifaciality factors, module transmission fractions, albedo values, and ground clearance. The analyses cover a wide spectrum of potential PV array archetypes through variation of installation parameters. The results of the Monte Carlo analysis reveal that the uncertainty is largely dependent on albedo uncertainty, but more simulations are needed to identify trends across system archetypes. The simulations are aimed at attributing an annual energy uncertainty factor for bifacial considerations that can be applied in post-processing of project probability of exceedance analysis.

energy modeling↗

DOE OSTI2022

Uncertainty Quantification of Bifacial Performance Modeling

bifacial↗

DOE OSTI2022

Uncertainty Quantification of Bifacial Performance Modeling: Preprint

energy modeling↗

DOE OSTI2022

Quantum simulations of dark sector showers

We consider dark sector scenarios where dark matter is accompanied by a dark photon and multiple-flavor dark fermions charged under the dark gauge group. We study quantum interference effects in dark sector jets, where multiple dark photons are emitted from high-energy dark fermions. We perform fully quantum simulations of dark sector showers and compare the results against those of the classical Monte-Carlo simulations. We find important differences in probability distributions of dark photon countings between quantum and classical computations. When the number of dark-fermion flavors is large, we find significant enhancements in large numbers of dark photon emissions. Such enhancements can provide distinguishing signals for our scenarios at particle colliders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Development of Accelerated Kinetic Monte Carlo Code for Simulation of Helium Bubble Evolution

A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Adjoint Klein-Nishina Sampling Methods: Efficiency, Speed and Applications

Three new rejection sampling methods for generating samples from the adjoint Klein-Nishina cross section are discussed: the two-branch rejection sampling procedure, the three-branch linear rejection sampling procedure and the three-branch inverse rejection sampling procedure. These methods have all been implemented in the Framework for REsearch in Nuclear ScIence and Engineering (FRENSIE). The efficiency and sample generation rate of each of these methods are evaluated to characterize the methods and to make recommendations regarding their use. The use of these methods in realistic transport simulations is also evaluated by incorporating a scattering function into the sampling process. Furthermore, the results of an infinite medium problem are presented to verify that the sampling procedure can be used in an adjoint Monte Carlo simulation to generate results that are in agreement with an equivalent forward simulation.

97 MATHEMATICS AND COMPUTING↗