SEARCH · Search NASA

Results for “Monte Carlo Simulation”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2026

IK-Frag: Frag data generator for the PHITS simulation with the inverse kinematic reaction producing a focused neutron beam

IK-Frag has been developed for the creation of the nuclear cross-section data format, which is named frag data and can be used in PHITS, a Monte Carlo simulation code. IK-Frag focuses on the inverse kinematic reactions between a lithium or beryllium ion and a proton target. These reactions achieve naturally collimated neutrons and potentially reduce the necessity of radiation shielding. IK-Frag enables PHITS users to conduct simulations for the inverse kinematic reactions. The present software aims to contribute to future development of the neutron source system using the inverse kinematic reactions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

NETL Natural Gas Lifecycle Model

conventional gas↗

DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗

DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗

DOE OSTI2025

QCD Theory Meets Information Theory

We present a novel technique to incorporate precision calculations from quantum chromodynamics into fully differential particle-level Monte Carlo simulations. By minimizing an information-theoretic quantity subject to constraints, our reweighted Monte Carlo incorporates systematic uncertainties absent in individual Monte Carlo predictions, achieving consistency with the theory input in precision and its estimated systematic uncertainties. Our method can be applied to arbitrary observables known from precision calculations, including multiple observables simultaneously. It generates strictly positive weights, thus offering a clear path to statistically powerful and theoretically precise computations for current and future collider experiments. As a proof of concept, we apply our technique to event-shape observables at electron-positron colliders, leveraging existing precision calculations of thrust. Our analysis highlights the importance of logarithmic moments of event shapes, which have not been previously studied in the collider physics literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Investigating the ecological fallacy through sampling distributions constructed from finite populations

Correlation coefficients and linear regression values computed from group averages can differ from correlation coefficients and linear regression values computed using individual scores. This observation known as the ecological fallacy often assumes that all the individual scores are available from a population. In many situations, one must use a sample from the larger population. In such cases, the computed correlation coefficient and linear regression values will depend on the sample that is chosen and the underlying sampling distribution. The sampling distribution of correlation coefficients and linear regression values for group averages will be identical to the sampling distribution for individuals for normally distributed variables for random samples drawn from infinitely large continuous distributions. However, data that is acquired in practice is often acquired when sampling without replacement from a finite population. Our objective is to demonstrate through Monte Carlo simulations that the sampling distributions for correlation and linear regression will also be similar for individuals and group averages when sampling without replacement from normally distributed variables. These simulations suggest that when a random sample from a population is selected, the correlation coefficients and linear regression values computed from individual scores will not be more accurate in estimating the entire population values compared to samples when group averages are used as long as the sample size is the same.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

NETL Natural Gas Lifecycle Model 2020 Rev1

This is the excel-based life cycle model that contains all the parameters and Monte Carlo simulation capabilities to model the techno-regions contained in the report: Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile. This is an update to a previous version (https://edx.netl.doe.gov/dataset/netl-natural-gas-lifecycle-model) and corrects a modeling error in the post-processing natural gas composition.

Calcium Chloride↗

DOE OSTI2025

Are better combinations of DERs more profitable?: Combinatorial optimization for aggregation of DERs in wholesale electricity markets

Recently, regulatory changes in various countries have enabled the participation of small-scale distributed energy resources (DERs) aggregated in virtual power plants (VPPs) in wholesale electricity markets. The inherent uncertainty and variability of resources comprising VPPs can lead to imbalances between forecasted and metered outputs, potentially resulting in the deficient settlement of generation under imbalance settlement rules. To address this challenge, it is essential to manage variability in the planning phase and uncertainty in the operation phase. Most current research focuses on managing forecasting errors in the operational phase, with insufficient attention given to the planning phase. Here, to bridge this gap, this paper proposes an optimal combination strategy for DERs to maximize the market participation revenue of VPPs by proactively managing variability in the planning phase. To estimate the expected revenue, we conducted analyses for homogeneous and heterogeneous DERs using Monte Carlo simulations and genetic algorithms. Remarkably, the proposed method demonstrated approximately 8 % higher revenue compared to the neighboring group case when considering diversity in DER set configuration with equal proportions of photovoltaics and wind.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

An Uncertainty-Informed and High-Fidelity Performance Forecasting Framework for Heliostat Fields

Concentrating Solar Thermal (CST) tower systems employ heliostat fields to direct solar energy to a central receiver, which then transfers the heat either directly to a thermal process (e.g., steam production) or to a thermal energy storage system for future use. Heliostat fields compose a significant proportion of the project costs of a CST tower system and the performance of the heliostats determines a plant's productivity at a given location. While CST characterization tools such as SolarPILOT and System Advisor Model (SAM) include a large collection of inputs that influence the performance of a CST tower system, many are uncertain prior to the development of the project and may have a significant impact on the overall energy delivery and profitability of a project; moreover, the fidelity of these models under default conditions may be insufficient to determine the value of component improvements such as those under development in the Heliostat Consortium. This work introduces a Monte Carlo simulation framework that incorporates uncertainty in key performance parameters to generate confidence intervals and percentile estimates for a CST solar field's energy delivery.

14 SOLAR ENERGY↗

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 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 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 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 OSTI2025

The role of excess vacancies in stabilizing solute clusters in low-alloy steels

The formation of solute clusters in irradiated low-alloy steels, such as reactor pressure vessel steels, is a critical cause of radiation hardening and embrittlement. However, the chemical interactions among solute elements and excess vacancies remain to be fully understood. This study employs density functional theory, cluster expansion, and lattice-based Monte Carlo simulations to examine the stability and morphology of nano-size coherent solute clusters with excess vacancies. The findings reveal that excess vacancies are crucial in stabilizing and promoting the growth of Mn-Ni-Si-vacancy clusters. A minimum of seven vacancies is required for stable nucleation and growth of these clusters. These Mn-Ni-Si clusters act as defect sinks, effectively trapping and absorbing mobile vacancies generated under irradiation. Furthermore, phosphorus (P) preferentially dissolves in Mn-Ni-Si clusters due to chemical coupling with vacancies. In conclusion, this study offers new insights into solute-vacancy interactions, enhancing our understanding of solute-defect cluster formation and embrittlement in low-alloy steels.

36 - MATERIALS SCIENCE↗

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 OSTI2025

Experimental study of energy-dependent angular broadening of MeV electron beams for high-resolution imaging in thick samples

In scanning transmission electron microscopy (STEM), spatial resolution is primarily influenced by the projected size of the electron probe within the specimen. In thin samples, a large semi-convergence angle enables a tightly focused beam and sub-nanometer resolution. However, in thick specimens, resolution is fundamentally limited by transverse beam broadening from multiple large-angle scattering events—for example, a probe with 10 mrad angular divergence can broaden by ∼100 nm over a 10 μm path. Since this broadening scales inversely with beam energy, MeV-STEM offers a promising route for high-resolution imaging in thick materials. To quantitatively assess this effect, we performed high-precision measurements at UCLA’s PEGASUS beamline, characterizing beam divergence and intensity profiles for 3–8 MeV electrons transmitted through a wedged-silicon sample of varying thickness. Our results reconcile discrepancies among analytical models and validate Monte Carlo simulations. Here, we find that increasing beam energy from 3.0 to 5.8 MeV reduces angular broadening by a factor of 2.6, with diminishing returns observed at 7.6 MeV. These findings provide a quantitative framework for optimizing MeV-STEM parameters in high-resolution imaging of thick biological and microelectronic specimens, and for guiding beam energy selection in other advanced imaging modes beyond STEM.

36 MATERIALS SCIENCE↗