Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct Simulation Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Design and simulation of a muon detector to characterize geological overburden

This study presents the design, construction, and simulation of a mobile muon detector tailored for geological overburden characterization. The detector employs plastic scintillator paddles with silicon photomultipliers (SiPMs) and a QuarkNet data acquisition system, offering a portable solution suitable for remote field deployment. The simulator’s modular aluminum frame allows for adjustable geometry and directional sensitivity, while its battery system supports over a week of autonomous operation. Preliminary experimental tests confirmed that its muon flux measurements were consistent with theoretical expectations. A comprehensive simulation framework using Geant4 and CORSIKA was developed to model detector response and overburden effects. Analytical and Monte Carlo methods were used to assess quadrant resolution and infer muon directionality. This work lays the foundation for future overburden mapping and supports the development of reconstruction algorithms for geological applications.

72 - PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Tetris-inspired detector with neural network for radiation mapping

Abstract Radiation mapping has attracted widespread research attention and increased public concerns on environmental monitoring. Regarding materials and their configurations, radiation detectors have been developed to identify the position and strength of the radioactive sources. However, due to the complex mechanisms of radiation-matter interaction and data limitation, high-performance and low-cost radiation mapping is still challenging. Here, we present a radiation mapping framework using Tetris-inspired detector pixels. Applying inter-pixel padding for enhancing contrast between pixels and neural networks trained with Monte Carlo (MC) simulation data, a detector with as few as four pixels can achieve high-resolution directional prediction. A moving detector with Maximum a Posteriori (MAP) further achieved radiation position localization. Field testing with a simple detector has verified the capability of the MAP method for source localization. Our framework offers an avenue for high-quality radiation mapping with simple detector configurations and is anticipated to be deployed for real-world radiation detection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Comparison of active interrogation methods for source location in a scattering and absorbing medium, consisting of PGNAA, and an AmBe quasi-forward biased directional source

Source location of Special Nuclear Material (SNM) encompassing 95% 235 U and 239 Pu is identified by utilizing a directional source from Patent No: US20190013109A1 using Prompt Gamma Neutron Activation Analysis (PGNAA) and neutron spectroscopy simulated with Monte-Carlo N-Particle transport 6.2 (MCNP). BC-408, HPGe, LaBr 3 detector arrays were used to identify the location of the SNM using total counts incident on each detector, and PGNAA photopeaks from HPGe and LaBr 3 detector arrays in a polyethylene shield. The conducted simulations varied the volume and location of the SNM in the MCNP input files to observe how the source location method behaved. PGNAA photopeaks used for source identification include 61 keV from fission, 2.223 MeV prompt gamma from hydrogen, 511 keV annihilation, and a single and double escape peaks from the prompt gamma interaction from hydrogen. The capabilities of each detector systems to acquire well resolved photopeaks with a 1% relative error or less, and total relative error for F4 and F8 tallies were less than 0.015% relative error. Source predictions of the SNM with uneven amounts of polyethylene shielding between the source and detectors was observed to overpredict and give invalid source location predictions. Source locations of the SNM with even amounts of polyethylene material between the source and each detector were found to be valid. With a 1 Ci 241 Am source activity, it was determined that 1630 s were needed to obtain the results for each detector system with the quasi-forward directional AmBe source. Coupling source and material identification together would increase acquisition time but would only require one system to determine.

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

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Feasibility of Observing Gamma-Ray Polarization from Cygnus X-1 Using a CubeSat

Instruments flown on CubeSats are small. Meaningful applications of CubeSats in astronomical observations rely on the choice of a particular subject that is feasible for CubeSats. Here we report the result of a feasibility study for observing gamma-ray polarization from Cygnus X-1 using a small Compton polarimeter on board a 3U CubeSat. Silicon detectors and cerium bromide scintillators were employed in the instrument models that we discuss in this study. Through Monte Carlo simulations with a Geant4-based MEGAlib package, we found that, with a 10 Ms on-axis, zenith-direction observation in a low-inclination, low-altitude, Earth-orbit radiation background environment, the minimum detectable polarization degree can be down to about 10% in 160–250 keV, 20% in 250–400 keV, and 65% in 400–2000 keV. A 3U CubeSat dedicated to observing Cygnus X-1 can therefore yield useful information on the polarization state of gamma-ray emissions from the brightest persistent X-ray black hole binary in the sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Framework to Analyze the Requirements of a Multiport Megawatt-Level Charging Station for Heavy-Duty Electric Vehicles

Widespread adoption of heavy-duty (HD) electric vehicles (EVs) will soon necessitate the use of megawatt (MW)-scale charging stations to charge high-capacity HD EV battery packs. Such a station design needs to anticipate possible station traffic, average and peak power demand, and charging/wait time targets to improve throughput and maximize revenue-generating operations. High-power direct current charging is an attractive candidate for MW-scale charging stations at the time of this study, but there are no precedents for such a station design for HD vehicles. We present a modeling and data analysis framework to elucidate the dependencies of a MW-scale station operation on vehicle traffic data and station design parameters and how that impacts vehicle electrification. This framework integrates an agent-based charging station model with vehicle schedules obtained through real-world vehicle telemetry data analysis to explore the station design and operation space. A case study applies this framework to a Class 8 vehicle telemetry dataset and uses Monte Carlo simulations to explore various design considerations for MW-scale charging stations and EV battery technologies. The results show a direct correlation between optimal charging station placement and major traffic corridors such as cities with ports, e.g., Los Angeles and Oakland. Corresponding parametric sweeps reveal that while good quality of service can be achieved with a mix of 1.2-megawatt and 100-kilowatt chargers, the resultant fast charging time of 35–40 min will need higher charging power to reach parity with refueling times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Effect of Confinement of SIA Cluster Diffusion by Impurities on Radiation Defect Accumulation due to 14MeV neutrons in Tungsten

Impurities acting as traps, reduce the self-interstitial atom (SIA) cluster mobility, thereby increasing recombination and in turn enhancing radiation-resistance of polycrystalline materials. However, at elevated temperatures, depending on SIA-trap (impurity) binding energy, SIA clusters may trap and detrap (detach) numerous times. In tungsten, SIA clusters of all sizes glide one-dimensionally (1D) along their Burgers vector direction. Hence, detrapped SIAs would retrace their path, confining their 1D-glide between traps. However, small SIA clusters can change the direction of their 1D-glide by overcoming a rotation barrier. In addition, impurities or solutes are known to reduce this barrier. A lower rotation barrier can result in detrapped SIA clusters that diffuse in a random 1D direction, effectively leading to 3D (net-3D) diffusion. Here we investigate the effect of two cases of SIA diffusion, namely confined-1D diffusion and net-3D diffusion, on the damage accumulation in polycrystalline tungsten under 14 MeV neutron irradiation using the object kinetic Monte Carlo method. Simulations were performed using cascade debris obtained from molecular dynamics simulations with primary knock-on atom (PKA) energies corresponding to 14 MeV neutrons. In the simulations with net-3D diffusion, SIA clusters of all sizes diffuse in a random 1D direction after the clusters detrap. While for the confined 1D diffusion case, SIA clusters larger than size 5 retain their original direction of 1D-glide. In both types of simulations, trapped vacancies are assumed to be permanently immobilized. Unexpectedly, the damage accumulation was lower with confined 1D diffusion than that with net-3D. We present a systematic comparison of the influence of the two cases of SIA diffusion on the radiation damage accumulation as a function of dose rate, detrapping barrier, and impurity concentration at 1025 K.

Accumulation and recovery of lattice damage, tungs↗
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 OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2026

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY↗
DOE OSTI2026

CADIS and FW-CADIS Variance Reduction in Gamma Transport for Predicting Prompt Forensics Signatures

The goal of prompt nuclear forensics is to determine the characteristics of a nuclear detonation based on the signatures available almost immediately after the explosion. An important characteristic is the reaction time history (RTH), a measure of the device’s rate of neutron multiplication. The RTH can be estimated by observation of the gamma radiation emitted from the detonation, which can be detected directly or observed indirectly as Teller light. Gamma transport simulations used to predict these radiation fields are often modeled stochastically using the Monte Carlo N-Particle (MCNP) code, which can be a computationally demanding task due to the number of particle histories needed to achieve statistical convergence. In an attempt to improve the efficiency of these calculations, we evaluate two variance reduction techniques: Consistent Adjoint-Driven Importance Sampling (CADIS) and Forward-Weighted Consistent Adjoint-Driven Importance Sampling (FW-CADIS). These methods use a deterministically calculated adjoint flux to create weight windows and source biasing that guide MCNP sampling. We study the utility of CADIS and FW-CADIS for their use in MCNP gamma transport for nuclear forensics prediction simulations. Furthermore, the results demonstrate that both CADIS and FW-CADIS improve the accuracy for forensics-focused simulations, with CADIS being most beneficial in direct detection and FW-CADIS being ideal for computing a global Teller light source.

CADIS↗
DOE OSTI2022

Induced seismicity potential based on probabilistic geomechanics for the San Juan Basin CarbonSAFE project

Although geologic carbon sequestration projects have yet to induce — or may never induce — a damaging earthquake, experiences from other deep injection industries such as hydraulic fracturing, enhanced geothermal systems, and saltwater disposal suggest that effective quantitative seismic risk assessment is necessary for deep saline carbon capture and sequestration (CCS) projects. One such imminent CCS project is the San Juan Basin CarbonSAFE Phase III program. We use Monte Carlo probabilistic geomechanical analyses combined with observations of the geologic and operational parameters of the San Juan Basin site and suggest that this project is of low-induced seismic risk. The primary analysis is split into four sections. First, we assess the literature for faults and past seismicity, and at least five faulting scenarios are directly relevant. Second, we develop and calibrate an integrated earth model for the project site. Third, we perform Monte Carlo simulations that consider reasonable uncertainties of the geomechanical parameters. Of five tested faulting scenarios, only one presents high Coulomb failure functions, but fourth, we determine the risk to be low based on the combined lack of historical seismicity, the geologic framework of the faulting scenario, and the presence of saltwater injection at the same depth as the proposed supercritical carbon dioxide injection. The most sensitive parameters in the geomechanical calculations are the fault dip and the coefficient of friction. The least sensitive are the fault strike and the orientation of the maximum horizontal principal stress.

58 GEOSCIENCES↗
DOE OSTI2023

High Intensity Gamma Ray Source Scintillation Attenuation Spectrometer and Filter Stack Monoenergetic Calibration

The Laser-Based X-ray Radiographic Imaging team at Los Alamos National Lab is looking for a monoenergetic MeV X-ray source to both verify our Monte Carlo N-particle simulations of detector performance and calibrate the instruments for future measurements. Funded by the Laboratory Directed Research and Development (LDRD) program, our overall goal is to improve the radiographic quality and reliability of laser-based X-ray sources for deployment at both dynamic and static radiography facilities. Laser-based X-ray sources have demonstrated smaller spot sizes to current electron accelerator based sources. The smaller laser spot size leads to significantly improved resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Rare Event Sampling Improves Mercury Instability Statistics

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare-event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately 10 –4 in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗