Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte-Carlo methods”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A comprehensive framework for validating simulation models of power system equipment using terminal measurements

Accurate simulation of power-plants is essential to the planning and operation of modern power grids. The current methods used to periodically check power-plant simulation models have many open questions about their limitations and accuracy. The research in this project explored using Monte-Carlo Experimentation (MCE) as a means for answering these important questions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Advancing X-ray quantum imaging through Monte-Carlo simulations

Imaging with X-rays poses fundamental limits due to radiation damage of the highly energetic photons. This becomes problematic for sensitive biological systems such as subcellular structures. Lowering the radiation dose, without sacrificing the signal-to-noise ratio, would be desirable for any kind of imaging modalities involving X-rays. To achieve this goal, quantum imaging with entangled X-ray photons constitutes a promising route. Production of biphotons have been demonstrated in the X-ray regime by the process of Spontaneous Parametric Down-Conversion (SPDC). However, compared to SPDC in the regime of visible light, the production rate for X-ray biphotons is extremely low. With the introduction of new high average brightness X-ray sources, such as 4th generation synchrotrons and high repetition rate Free-Electron X-ray Lasers (XFEL), quantum imaging may become practical. We introduce a ray tracing approach using Monte-Carlo sampling, specifically designed for quantum imaging with entangled X-ray photons generated by SPDC. By simulation, the superior image quality of quantum over classical imaging methods is demonstrated using realistic experimental conditions available at high repetition rate XFELs. With these simulations, we can efficiently assist the design of future experiments at beam lines, which can substantially accelerate the advancement of X-ray quantum imaging and reduce costs.

Entangled Photons↗
NASA NTRSDate not supplied

Design, Development, and Test of the Advanced Apollo Orbital Assembly System

As spaceflight moves toward commercial solutions for Crew Vehicles and Space Stations, opportunity exists to lower costs with novel designs. Probe and cone docking systems provide a lightweight, low cost, and high-performance docking solution. This work revisits the Apollo probe and cone design and modifies it for the requirements of today’s computer-controlled spacecraft. This new system is called the Advanced Apollo Orbital Assembly (APOA) system, and is intended to support satellite, cargo, and space station assembly missions. A simulation of the APOA was programmed, components were sized, and a design was created by scaling the probe to the size of an EELV Secondary Payload Adapter (ESPA) tunnel. Prototype test hardware was 3D printed using Fused Deposition Modeling (FDM) methods with Polylactic Acid (PLA) material. Testing of the APOA-ESPA was conducted at Marshall Space Flight Center’s (MSFC) Flat Floor, and a test-correlated simulation is used to evaluate a Monte-Carlo of Initial Contact Conditions (ICC’s) to establish baseline performance. The successful development, test, and correlation of the APOA-ESPA proves the design validity and increases Technology Readiness Level (TRL) from 2 to 4. This work opens the door to construction of an APOA-ESPA from flight like materials, and to develop a larger scale prototype APOA. When the full scale APOA is incorporated with the Common Berthing Mechanism (CBM), becomes the Hybrid Berthing System (HBS), which allows for berthing without a robotic arm.

Berthing↗
DOE OSTI2024

G EANT 4 atomic relaxation data for transfermium nuclei (Z = 101–104)

Advanced theoretical methods can accurately calculate various atomic observables and predict electronic structure. Still, systematic computations of the radiative and non-radiative transition probabilities and energies are missing for the actinides and all the transfermium elements. However, these compilations are needed for comprehensive Monte-Carlo simulations (such as GEANT4) of the radioactive decay of transfermium nuclei. These simulations can forma basis for data analysis of experiments, especially with complex detection setups. Investigation of the transfermium nuclei is crucial for understanding the nature of the nuclear force. In this study, simulations based on data from the Jena Atomic Calculator (JAC) and the data from the Evaluated Atomic Data Library (EADL) present in GEANT4 were found compatible for the three elements Ba(Z = 56), U(Z = 92), and Fm(Z = 100), thus, validating the JAC calculations. For Z> 100, we also found sound agreement between simulations that used data generated with JAC and experimental results involving No(Z = 102) and Rf(Z = 104) isotopes. In conclusion, these results demonstrate that JAC can produce reliable atomic data sets for transfermium elements, which will assist in analyzing nuclear-decay-spectroscopy experiments.

GEANT4↗
NASA NTRS1991

Translating expert system rules into Ada code with validation and verification

The purpose of this ongoing research and development program is to develop software tools which enable the rapid development, upgrading, and maintenance of embedded real-time artificial intelligence systems. The goals of this phase of the research were to investigate the feasibility of developing software tools which automatically translate expert system rules into Ada code and develop methods for performing validation and verification testing of the resultant expert system. A prototype system was demonstrated which automatically translated rules from an Air Force expert system was demonstrated which detected errors in the execution of the resultant system. The method and prototype tools for converting AI representations into Ada code by converting the rules into Ada code modules and then linking them with an Activation Framework based run-time environment to form an executable load module are discussed. This method is based upon the use of Evidence Flow Graphs which are a data flow representation for intelligent systems. The development of prototype test generation and evaluation software which was used to test the resultant code is discussed. This testing was performed automatically using Monte-Carlo techniques based upon a constraint based description of the required performance for the system.

Becker, Lee↗
DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗
DOE OSTI2026

Double-real corrections to color singlet decay in a parton-shower inspired scheme

We introduce a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, with counterterms based on scalar radiators and pure splitting functions. Overlapping singularities in the multipole radiation pattern are disentangled by partial fractioning, and the kinematics mapping corresponds to iterated next-to-leading order kinematics. We verify that the double-real remainder to $e^+e^-\to\;q\bar{q}$ is rendered finite in the single and double unresolved limits and investigate the numerical convergence of the Monte-Carlo integral. We compute the phase-space integrals of the scalar counterterms in the back-to-back configuration, both analytically and with the help of numerical techniques based on sector decomposition.

Campbell, John M. [Fermilab]↗
NASA NTRS1978

Performance evaluation of integrated guidance and navigation systems

The method of linear covariance analysis is applied to the study of an aerospace vehicle integrated guidance-and-navigation system. This requires that the rms trajectory deviations of the vehicle be computed along with the rms errors of the navigation system. The technique is used to examine the performance of the shuttle orbiter vehicle's guidance-and-navigation system during the terminal phase of the landing. Essentially the same rms trajectory deviations were obtained from a single covariance run as from a 30-sample Monte-Carlo run set. Of particular interest in the study was the result that rms trajectory deviations for the integrated-system were generally significantly larger than the rms navigation errors.

Tao, Y. C.↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2025

Particle-in-cell simulations of burning inertial confinement fusion capsule implosions

Anomalies observed in the neutron spectral shift of high-yield shots at the National Ignition Facility (NIF) suggest the presence of suprathermal ions [E. P. Hartouni et al., Nat. Phys. 19, 72 (2023)], implying that kinetic effects play a significant role in burning inertial confinement fusion (ICF) plasmas. Furthermore, recent measurements of reaction-in-flight (RIF) neutrons offer a direct probe of the stopping power in the burning fuel region of high-energy alpha particles and up-scattered fuel ions. We have developed the particle-in-cell code PICNIC, an exactly energy-conserving particle-in-cell Monte-Carlo collision (PIC-MCC) code to simulate the burn stage in ICF. We present results from 1D spherical simulations of NIF shot N210808, which was the first to exceed the Lawson criterion for laser fusion. We find that the suprathermal ions generated by large-angle Rutherford and nuclear elastic scattering (NES) with fusion alphas produce an alpha knock-on neutron (AKN) signal consistent with the extent and relative yield of the AKN spectrum identified in ignition experiments at the NIF. Furthermore, we also find that the inclusion of large-angle scattering physics does not explain the anomalously large spectral shift observed in the experiment.

High-energy-density plasmas↗
NASA NTRS1989

Monte-Carlo modeling of polar wind photoelectron distributions with anomalous heat flux

In situ measurements above the polar cap by the DE satellites show asymmetric field-aligned electron velocity distributions in the photoelectron energy range 5-60 eV. The formation of these anisotropic distributions is described by a collisional kinetic equation that includes the effects of Coulomb collisions (via a Fokker-Planck operator), those of a parallel electric field, and the magnetic mirror effects of the geomagnetic field. This equation is solved using a Monte Carlo simulation method. The simulation illustrates the formation of the portion of the observed electron distributions originating in the ionosphere and makes it possible to evaluate the significance of the suprathermal electrons at higher altitude. Both observations and calculations indicate that these electrons possess an energy flux which, by its effect on the polar-wind electric field, may influence the dynamics of the polar-wind outflow.

Yasseen, F.↗