Search NASA⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An improved stochastic weighted particle method for boundary driven flows

Here, the stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Chemical Design of Spin Frustration to Realize Topological Spin Glasses

Patterning spins to generate collective behavior is at the core of condensed matter physics. Physicists develop techniques, including the fabrication of magnetic nanostructures and precision layering of materials specifically to engender frustrated lattices. As chemists, we can access such exotic materials through targeted chemical synthesis and create new lattice types by chemical design. Here, in this study, we introduce a new approach to induce magnetic frustration on a modified honeycomb lattice through a competition of alternating antiferromagnetic (AFM) and ferromagnetic (FM) nearest-neighbor interactions. By subtly modulating these two types of interactions through facile synthetic modifications, we created two systems: (1) a topological spin glass and (2) a frustrated spin-canted magnet with low-temperature exchange bias. To design this unconventional magnetic lattice, we used a metal–organic framework (MOF) platform, Ni 3 (pymca) 3 X 3 (NipymcaX where pymca = pyrimidine-2-carboxylato and X = Cl, Br). We isolated two MOFs, NipymcaCl and NipymcaBr, featuring canted Ni 2+ -based moments. Despite this similarity, differences in the single-ion anisotropies of the Ni 2+ spins result in distinct magnetic properties for each material. NipymcaCl is a topological spin glass, while NipymcaBr is a rare frustrated magnet with low-temperature exchange bias. Density functional theory calculations and Monte Carlo simulations on the NipymcaX lattice support the presence of magnetic frustration as a result of alternating AFM and FM interactions. Our calculations enabled us to determine the ground-state spin configuration and the distribution of spin–spin correlations relative to paradigmatic kagomé and triangular lattices. This modified honeycomb lattice is similar to the electronic Kekulé-O phase in graphene and provides a highly tunable platform to realize unconventional spin physics.

Exchange bias↗
DOE OSTI2025

Rate Limiting Regimes in Photochemical H2 Generation by Complexes of Colloidal CdS Nanorods and Hydrogenase

Driving redox enzyme catalysis with photoexcited semiconductor nanocrystals is a compelling approach for chemical conversion. We examined how the interplay of the many chemical steps involved determines the rates of photochemical H2 production with complexes of colloidal CdS nanorods and an [FeFe]-hydrogenase. We elucidated the roles of three critical and previously elusive processes-scavenging of photoexcited holes from nanorods, back-electron transfer, and H2 oxidation. Kinetic Monte Carlo simulations and fitting to experimental data revealed that hole transfer becomes the rate-limiting step at high illumination intensities. Comparisons of simulations to experimental H2 production showed that both back-electron transfer and H2 oxidation play an efficiency-limiting role at high catalyst loadings. This work provides guiding principles for tuning experimental parameters to minimize energy-wasting pathways and optimize photochemical product formation. More broadly, we demonstrate how critical but elusive chemical steps in photochemical reactions can be probed with a combination of experiments and simulations.

08 HYDROGEN↗
DOE OSTI2026

Techno-Economic Analysis of Recycling Strategies for Catalyst and Acid During Catalytic Graphitization

With the aim of meeting the urgent demand for active anode materials (AAM) in energy storage systems, bio-based graphite (biographite) emerges as an affordable solution to de-risk the turbulent supply chain of critical minerals. Anode grade biographite requires high crystallinity and purity, which can be achieved by catalytic graphitization with iron, followed by acid washing. Therefore, a well-conceived process integration that recycles catalyst can be the starting point to commercialization. This study evaluates closed-loop catalyst recovery, and byproducts valorization scenarios through a technoeconomic framework to help understand the scale-up potential of biographite. For the acid washing, three reactors in series meet the required biographite purity at 99.95%. Iron and acid recovery can reduce material consumption and waste generation by ~95%, albeit at the expense of ~80% increase in capital costs. Recovery scenarios present similar capital and operational expenses, yielding minimum selling prices (MSP) near $6 kg-1 of biographite. Monte Carlo methodology reveals that feedstock price accounts for ~60% of MSP variance, followed by plant capacity ~20%. The likelihood of reaching a competitive profit margin of 30% in the U.S. AAM market sits at 85% average for recovery scenarios, and 103% when iron oxide is sold as byproduct. Additionally, an IRR >= 15% can be achieved for half of Monte Carlo simulations, representing promising early-stage results. Biographite production offers a strategic pathway to stabilize the anode market beyond China by integrating established technologies for a scalable, economically viable, and sustainable process. The role of catalyst recovery and byproducts utilization is critical for advancing the biomaterials industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluation of a Variable-Aperture Full-Ring SPECT System using Large-Area Pixelated CZT Modules: A Simulation Study for Brain SPECT Applications

Single photon emission computed tomography (SPECT) scanner using cadmium zinc telluride (CZT) offers improved imaging capability over conventional NaI (Tl)-based SPECT scanner. We aim to demonstrate a full-ring SPECT system design with eight large-area CZT detectors that can be used for a broad spectrum of SPECT radiopharmaceuticals and provides higher sensitivity and better spatial resolution than those of conventional NaI (Tl)-based gamma cameras. A newly-designed full-ring SPECT system is composed of 8 large-area CZT cameras (128 mm × 179.2 mm effective area) that can be independently swiveled around their own axes of rotation independently and can have radial motion for varying aperture sizes that can be adapted to different sizes of imaging volume. Extended projection data were generated by conjoining projections of two adjacent detectors to overcome the limited field-of-view (FOV) by each CZT camera. Using Monte Carlo simulations, we evaluated this new system design with digital phantoms including a Derenzo hot-rod phantom and a Zubal brain phantom. Comparison of performance metrics such as spatial resolution, sensitivity, contrast-to-noise ratio (CNR), and contrast-recovery ratio were made between our design and a conventional SPECT scanner. The proposed scanner could result in up to about 3 times faster in acquisition time over conventional scan time at same acquisition time per step. The spatial resolution of our proposed scanner was similar or better to that of the conventional scanner, and there was significant performance improvement over the conventional scanner particularly in sensitivity (approximately 4 times). Overall, we successfully reconstructed the phantom image for both 99 mTc-based perfusion and 123 I-based dopamine transporter (DaT) brain studies simulated for our new design. In particular, the striatal/background contrast-recovery ratio in 3-to-1 reference ratio was over the 0.8 for the 123 I-based DaT study. In conclusion, we demonstrated the potential of our new full-ring CZT SPECT design, showing improved performance metrics such as improved system sensitivity and CNR while maintaining other important imaging parameters such as spatial resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Life-cycle cost-benefit analysis of a novel self-heating pavement made from coal-derived solid carbon

Winter weather events challenge the safety, efficiency, and sustainability of roadway networks, particularly road bridges vulnerable to climate impacts. Although effective, conventional de-icing methods incur high expenses and cause significant environmental contamination. To address these longstanding issues, a low-cost novel coal-derived carbon enabled smart pavement (CDC-SP) de-icing system was developed, aiming to alleviate traffic delays, severe corrosion, compromised safety, and environmental pollution. However, the economic performance of CDC-SP has not been quantified to guide decision-making. To this end, this paper presents a Life-Cycle Cost-Benefit Analysis (LCCBA) to assess the economic feasibility of the CDC-SP de-icing system for road bridges in five representative urban and rural areas across cool humid and cold humid climate zones. Given the variations and uncertainties of several factors that determine the costs and benefits, a sensitivity analysis was conducted to determine the system’s economic reliability. Additionally, a Monte Carlo simulation (MCS) was performed to quantify the impacts of important factors and identify the most beneficial scenarios. Furthermore, the CDC-SP de-icing systems have demonstrated substantial economic benefits and short payback periods, showing great applicability for rural and urban road bridges with different service levels across multiple climate zones as compared to conventional de-icing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coexistence and Interaction of Spinons and Magnons in an Antiferromagnet with Alternating Antiferromagnetic and Ferromagnetic Quantum Spin Chains

In conventional quasi-one-dimensional antiferromagnets with quantum spins, magnetic excitations are carried by either magnons or spinons in different energy regimes: they do not coexist independently, nor could they interact with each other. In this Letter, by combining inelastic neutron scattering, quantum Monte Carlo simulations and Random Phase Approximation calculations, we report the discovery and discuss the physics of the coexistence of magnons and spinons and their interactions in Botallackite-Cu 2 (OH) 3 Br. This is a unique quantum antiferromagnet consisting of alternating ferromagnetic and antiferromagnetic Spin-1/2 chains with weak interchain couplings. Furthermore, our study presents a new paradigm where one can study the interaction between two different types of magnetic quasiparticles: magnons and spinons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Edge-Mediated Annihilation of Vacancy Clusters in Monolayer Molybdenum Diselenide (MoSe 2 ) under Electron Beam Irradiation

Here annihilation of vacancy clusters in monolayer molybdenum diselenide (MoSe2) under electron beam irradiation is reported. In situ high-resolution transmission electron microscopy observation reveals that the annihilation is achieved by diffusion of vacancies to the free edge near the vacancy clusters. Monte Carlo simulations confirm that it is energetically favorable for the vacancies to locate at the free edge. By computing the minimum energy path for the annihilation of one vacancy cluster as a case study, it is further shown that electron beam irradiation and pre-stress in the suspended MoSe2 monolayer are necessary for the vacancies to overcome the energy barriers for diffusion. The findings suggest a new mechanism of vacancy healing in 2D materials and broaden the capability of electron beam for defect engineering of 2D materials, a promising way of tuning their properties for engineering applications.

2D materials↗
DOE OSTI2022

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , C 3 H 8 , and CO 2 across a temperature range of 250–400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure–temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption–desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300–550 K. Here, a temperature threshold of ~400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

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

Unravelling Magnetic Nanochain Formation in Dispersion for In Vivo Applications

Self-assembly of iron oxide nanoparticles (IONPs) into 1D chains is appealing, because of their biocompatibility and higher mobility compared to 2D/3D assemblies while traversing the circulatory passages and blood vessels for in vivo biomedical applications. In this work, parameters such as size, concentration, composition, and magnetic field, responsible for chain formation of IONPs in a dispersion as opposed to spatially confining substrates, are examined. In particular, the monodisperse 27 nm IONPs synthesized by an extended LaMer mechanism are shown to form chains at 4 mT, which are lengthened with applied field reaching 270 nm at 2.2 T. The chain lengths are completely reversible in field. Using a combination of scattering methods and reverse Monte Carlo simulations the formation of chains is directly visualized. The visualization of real-space IONPs assemblies formed in dispersions presents a novel tool for biomedical researchers. This allows for rapid exploration of the behavior of IONPs in solution in a broad parameter space and unambiguous extraction of the parameters of the equilibrium structures. Additionally, it can be extended to study novel assemblies formed by more complex geometries of IONPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗
DOE OSTI2021

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

We present the XFaster analysis package, XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-C ℓ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data (Netterfield et al. 2002; Jones et al. 2006), and also used as part of the Planck analysis (Rocha et al. 2011). Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument (SPIDER Collaboration 2021). The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Magnetic phase diagram of A 2 [FeCl 5 (H 2 O)] ( A = K, Rb, NH 4 )

Here, erythrosiderites with the formula A 2 FeX 5 ∙H 2 O, where A = Rb, K, and (NH 4 ) and X = Cl and Br are intriguing systems that possess various magnetic and electric phases, as well as multiferroic phases in which magnetism and ferroelectricity are coupled. In this report, we study the magnetic phase diagram of erythrosiderites as a function of superexchange interactions and magnetic anisotropies. To this end, we perform classical Monte Carlo simulations on magnetic Hamiltonians that contain five different superexchange interactions with single-ion anisotropies. Our phase diagram contains all magnetic ground states that have been experimentally observed in these materials. We argue that the ground states can be explained by varying the ratio of $\frac{J_4}{J_2}$. For $\frac{J_4}{J_2}$ > 0.95 a cycloidal spins structure is stabilized as observed in (NH 4 ) 2 FeCl 5 ∙H 2 O and otherwise a collinear spin structure is stabilized as observed in (K,Rb) 2 FeCl 5 ∙H 2 O. We also show that the difference in the single-ion anisotropy along a- and c- axes is essential to stabilize the intermediate state observed in (NH) 2 FeCl 5 ∙H 2 O.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Adjoint DSMC for nonlinear spatially-homogeneous Boltzmann equation with a general collision model

We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in Caflisch et al., which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Range cost-effectiveness of plug-in electric vehicle for heterogeneous consumers: An expanded total ownership cost approach

Plug-in electric vehicles (PEV) appears to have sales momentum in major personal vehicle markets but are still at the early market stage. Opportunities to accelerate PEV adoption can be discovered through comprehensive total cost of ownership (TCO) analysis. Understanding the cost-effective electric ranges of PEVs for consumers, manufacturers, and the society is critical for any discussion of PEV mass markets. This study expanded the traditional TCO approach by (1) fully considering heterogeneous consumer driving patterns, (2) quantifying the charging inconvenience and range anxiety cost of battery-electric vehicles (BEVs), and (3) monetizing both tangible and intangible PEV policies. Uncertainties were handled through Monte Carlo simulation. The results suggest that BEVs with an electric range of 250–350 km have the lowest TCO in cities with government-enacted purchase limitations, and internal combustion engine vehicles (ICEVs) have the lowest TCO in cities without purchase limitations, even when considering PEV subsidies. The lowest TCO for some consumer groups is obtained by BEVs with an electric range of 400–450 km, especially in northern China, where the weather is colder. The cost-effective all-electric range for BEVs in each city in 2025 will decrease due to improved battery performance in cold environments and an expanded charging infrastructure. Finally, based on TCO, plug-in hybrid electric vehicles (PHEVs) are currently more suitable for drivers with a high average daily mileage or a large mileage variance. However, by 2025, BEVs with a long driving range may become a more cost-effective choice for these drivers.

33 ADVANCED PROPULSION SYSTEMS↗