Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 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↗
DOE OSTI2024

Dosimetric and biological impact of activity extravasation of radiopharmaceuticals in PET imaging

The increasing use of nuclear medicine and PET imaging has intensified scrutiny of radiotracer extravasation. To our knowledge, this topic is understudied but holds great potential for enhancing our understanding of extravasation in clinical PET imaging. This work aims to (1) quantify the absorbed doses from radiotracer extravasation in PET imaging, both locally at the site of extravasation and with the extravasation location as a source of exposure to bodily organs and (2) assess the biological ramifications within the injection site at the cellular level. A radiation dosimetry simulation was performed using a whole-body 4D Extended Cardiac-Torso (XCAT) phantom embedded in the GATE Monte Carlo platform. A 10-mCi dose of 18 F-FDG was chosen to simulate a typical clinical PET scan scenario, with 10% of the activity extravasated in the antecubital fossa of the right arm of the phantom. The extravasation volume was modeled as a 5.5 mL rectangle in the hypodermal layer of skin. Absorbed dose contributions were calculated for the first two half-lives, assuming biological clearance thereafter. Dose calculations were performed as absorbed doses at the organ and skin levels. Energy deposition was simulated both at the local extravasation site and in multiple organs of interest and converted to absorbed doses based on their respective masses. Each simulation was repeated ten times to estimate Monte Carlo uncertainties. Biological impacts on cells within the extravasated volume were evaluated by randomizing cells and exposing them to a uniform radiation source of 18 F and 68 Ga. Particle types, their energies, and direction cosines were recorded in phase space files using a separate Geant4 simulation to characterize their entry into the nucleus of the cellular volume. Subsequently, the phase space files were imported into the TOPAS-nBio simulation to assess the extent of DNA damage, including double-strand breaks (DSBs) and single-strand breaks (SSBs). Organ-level dosimetric estimations are presented for 18 F and 68 Ga radionuclides in various organs of interest. With 10% extravasation, the hypodermal layer of the skin received the highest absorbed dose of 1.32 ± 0.01 Gy for 18 F and 0.99 ± 0.01 Gy for 68 Ga. The epidermal and dermal layers received absorbed doses of 0.07 ± 0.01 Gy and 0.13 ± 0.01 Gy for 18 F, and 0.14 ± 0.01 Gy and 0.29 ± 0.01 Gy for 68 Ga, respectively. In the extravasated volume, 18 F caused an average absorbed dose per nucleus of 0.17 ± 0.01 Gy, estimated to result in 10.58 ± 0.50 DSBs and 268.11 ± 12.43 SSBs per nucleus. For 68 Ga, the absorbed dose per nucleus was 0.11 ± 0.01 Gy, leading to an estimated 6.49 ± 0.34 DSBs and 161.24 ± 8.12 SSBs per nucleus. Absorbed doses in other organs were on the order of micro-gray (µGy). The likelihood of epidermal erythema resulting from extravasation during PET imaging is low, as the simulated absorbed doses to the epidermis remain below the thresholds that trigger such effects. Moreover, the organ-level absorbed doses were found to be clinically insignificant across various simulated organs. The minimal DNA damage at the extravasation site suggests that long-term harm, such as radiation-induced carcinogenesis, is highly unlikely.

DNA strand breaks↗
DOE OSTI2021

Nematic order driven by superconducting correlations

The interplay of nematicity and superconductivity has been observed in a wide variety of quantum materials. To explore this interplay, we consider a two-dimensional (2D) array of nematogens, local droplets with Z 3 nematicity, coupled to a network of Josephson junction wires. Using finite temperature classical Monte Carlo simulations, we elucidate the phase diagram of this model and show that the development of superconducting correlations and the directional delocalization of Cooper pairs can promote nematogen ordering, resulting in long-range nematic order. We obtain the transport properties of our model within an effective resistor network picture. We discuss these ideas in the context of the 2D electron gas at the (111) KTaO 3 interface and the doped topological insulators Nb x Bi 2 Se 3 and Cu x Bi 2 Se 3 . Overall, our work makes contact with Phil Anderson’s numerous contributions to broken symmetries driven by the saving of kinetic energy, including double exchange ferromagnetism and the interlayer tunneling theory of high Tc csuperconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

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

Designs and Modeling of Architect Multi-material Scintillator System

We present a new conceptual radiation detector, the Architected Multimaterial Scintillator System, that utilizes a scintillator composed of multiple materials arranged in architected structures to enable new capabilities. By structuring differently-dyed materials, the wavelength of the scintillation light encodes additional information in radiation measurements. These structures can be realized through additive manufacture (3D-printing). Two classes of this concept are described and modelled using Monte Carlo simulations to evaluate their performance. The first detector design uses dye microstructures to encode particle tracking information, allowing for directional neutron detection and gamma/neutron discrimination. The second design uses a dye gradient to indicate the position of radiation along the gradient. The simulation results indicate this new concept in radiation detection can achieve strong performance in a variety of capabilities including particle identification, directionality and spectroscopy measurements, and particle position reconstruction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

On the Relationship Between Shallow Cumulus Cloud Field Properties and Surface Solar Irradiance

Abstract Shallow cumulus clouds exhibit highly three‐dimensional (3‐D) spatial structure leading to complex variability in the surface solar irradiance (SSI) beneath. This variability is captured by the typically bimodal shape of the SSI probability density function (PDF). Using large eddy simulation to generate well‐resolved cloud fields and Monte Carlo 3‐D radiative transfer to reproduce realistic associated SSI PDFs, we seek direct relationships between the cloud field properties and the SSI PDF shape. Applying both random forest and artificial neural network algorithms, we find variations in the two modes of the SSI PDF are well predicted by just a handful of cloud field properties. The two algorithms utilize cloud properties similarly, with indistinguishable performance despite their different architectures. These results offer a marked improvement in realism relative to one‐dimensional radiative transfer while bypassing computationally expensive 3‐D radiative transfer, with immediate application to renewable energy assessments, and potential for several other geophysical applications.

58 GEOSCIENCES↗
DOE OSTI2021

High-resolution planar electron beam induced current in bulk diodes using high-energy electrons

Understanding the impact of high-energy electron radiation on device characteristics remains critical for the expanding use of semiconductor electronics in space-borne applications and other radiation harsh environments. Here in this paper, we report on in situ measurements of high-energy electron radiation effects on the hole diffusion length in low threading dislocation density homoepitaxial bulk n-GaN Schottky diodes using electron beam induced current (EBIC) in high-voltage scanning electron microscopy mode. Despite the large interaction volume in this system, quantitative EBIC imaging is possible due to the sustained collimation of the incident electron beam. This approach enables direct measurement of electron radiation effects without having to thin the specimen. Using a combination of experimental EBIC measurements and Monte Carlo simulations of electron trajectories, we determine a hole diffusion length of 264 ± 11 nm for n-GaN. Irradiation with 200 kV electron beam with an accumulated dose of 24 × 10 16 electrons cm -2 led to an approximate 35% decrease in the minority carrier diffusion length.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Monte Carlo Modeling and Design of Photon Energy Attenuation Layers for >10× Quantum Yield Enhancement in Si-Based Hard X-ray Detectors

High-energy (>20 keV) X-ray photon detection at high quantum yield, high spatial resolution, and short response time has long been an important area of study in physics. Scintillation is a prevalent method but limited in various ways. Directly detecting high-energy X-ray photons has been a challenge to this day, mainly due to low photon-to-photoelectron conversion efficiencies. Commercially available state-of-the-art Si direct detection products such as the Si charge-coupled device (CCD) are inefficient for >10 keV photons. Here, we present Monte Carlo simulation results and analyses to introduce a highly effective yet simple high-energy X-ray detection concept with significantly enhanced photon-to-electron conversion efficiencies composed of two layers: a top high-Z photon energy attenuation layer (PAL) and a bottom Si detector. We use the principle of photon energy down conversion, where high-energy X-ray photon energies are attenuated down to ≤10 keV via inelastic scattering suitable for efficient photoelectric absorption by Si. Our Monte Carlo simulation results demonstrate that a 10–30× increase in quantum yield can be achieved using PbTe PAL on Si, potentially advancing high-resolution, high-efficiency X-ray detection using PAL-enhanced Si CMOS image sensors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗