SEARCH · Search NASA

Results for “classical 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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Continuous-Energy Verification of MCNP Calculations Using One-Group Spherical and Slab Criticality Benchmarks

This work presents a continuous-energy Monte Carlo verification study of one group spherical and slab criticality benchmarks using the MCNP ®1 code. Classical tabulations and newly generated benchmark solutions obtained by direct numerical evaluation by the author are considered. The benchmarks span weakly to strongly multiplying regimes and provide analytically defined critical radii as functions of a single parameter, c .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1973

Dissociation rate of bromine diatomics in an argon heat bath

The evolution of a collection of 300 K bromine diatomics embedded in a heat bath of argon atoms at 1800 K was studied by computer, and a dissociation-rate constant for the reaction Br2 + BR + Ar yields Br + Ar was determined. Previously published probability distributions for energy and angular momentum transfers in classical three-dimensional Br2-Ar collisions were used in conjunction with a newly developed Monte Carlo scheme for this purpose. Results are compared with experimental shock-tube data and the predictions of several other theoretical models. A departure from equilibrium is obtained which is significantly greater than that predicted by any of these other theories.

Razner, R.↗

DOE OSTI2025

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1975

Analysis of a semiclassical model for rotational transition probabilities

A semiclassical model proposed by Pearson and Hansen (1974) for computing collision-induced transition probabilities in diatomic molecules is tested by the direct-simulation Monte Carlo method. Specifically, this model is described by point centers of repulsion for collision dynamics, and the resulting classical trajectories are used in conjunction with the Schroedinger equation for a rigid-rotator harmonic oscillator to compute the rotational energy transition probabilities necessary to evaluate the rotation-translation exchange phenomena. It is assumed that a single, average energy spacing exists between the initial state and possible final states for a given collision.

Deiwert, G. S.↗

NASA NTRSDate not supplied

Analysis and Application of Novel and Heritage Acceleration Limiting Algorithms for SLS on Ascent

National Aeronautics and Space Administration (NASA) is currently building the Space Launch System (SLS) Block-1 launch vehicle, to be used as the crewed heavy-lift vehicle for the Artemis series of missions. The SLS Block-1 guidance subsystem utilizes a nonlinear algorithm derived from Shuttle heritage for limiting the vehicle’s induced maximum axial acceleration during ascent flight, known as g-limiting. Even though g-limiting has demonstrated stability and robustness through several design analysis cycles and the extensive Shuttle flight history, there are no available documents that demonstrate that this algorithm has been proven stable through conventional controls stability analysis. This paper highlights the non-linear nature of g-limit, presents an alternative methodology to employ a linear version of this algorithm, performs assessment of selected linear gains using classical stability analysis and, conducts a comparison of both approaches through Six Degrees-of-Freedom (6-DOF) Monte Carlo (MC) simulations.

Jason Everett↗

NASA NTRS2016

Spin Glass Patch Planting

In this paper, we propose a patch planting method for creating arbitrarily large spin glass instances with known ground states. The scaling of the computational complexity of these instances with various block numbers and sizes is investigated and compared with random instances using population annealing Monte Carlo and the quantum annealing DW2X machine. The method can be useful for benchmarking tests for future generation quantum annealing machines, classical and quantum mechanical optimization algorithms.

Quantum Annealing↗

DOE OSTI2026

Accelerating multicanonical sampling with irreversibility

Flat-histogram Monte Carlo simulations are well-established, robust methods to perform random walks in a physical observable or parameter space, making them suitable for finding ground states or studying phase transitions in complex systems in statistical physics. However, their efficiency can be limited by the time to attain the desired flat distribution, which is generally unknown prior to the simulations. In particular, they might suffer from slowing down towards the end of a simulation due to the diffusive nature of random walks. In this work we apply irreversibility to the multicanonical Monte Carlo method via the lifting approach to alleviate this behavior. We achieve a 2–4 times speedup in ground-state search for a two-dimensional (2D) Ising model, and up to an order of magnitude of speedup for finding the ground-state energy in an Edwards–Anderson spin glass, compared to traditional multicanonical sampling. In conclusion, the round-trip times between ground states show a narrower distribution and are significantly shorter compared to the reversible counterpart, suggesting that a lower convergence time with a smaller time variance is feasible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Bayesian optimization of PYTHIA 8 tunes

A new tune (set of model parameters) is found for the six most important parameters of the PYTHIA 8 final state parton shower and hadronization model using Bayesian optimization. The tune fits the Large Electron-Positron collider (LEPI) data from ALEPH better than the default tune in PYTHIA 8. To the best of our knowledge, we present the most comprehensive application of Bayesian optimization to the tuning of a parton shower and hadronization model using the LEPI data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Propagation of Noise Uncertainty Through Virtual Strain Gauge Formulations for 2D Digital Image Correlation

The effect of displacement uncertainty is examined on 2-dimensional strain, calculated using linear surfaces fitted to the displacement field. A classical engineering error propagation method is used to calculate uncertainty in Green-Lagrangian strain calculations. The derived uncertainty is compared to a Monte Carlo simulation and discrepancies under 2% are seen between these two methods. The effect of virtual strain gauge size, displacement uncertainty, and boundaries on the region of interest on the strain uncertainty are considered. Here, an exponential decay relationship is observed between strain uncertainty and virtual strain gauge size, while a linear relationship is seen between strain and displacement uncertainty. For boundaries in the region of interest, strain uncertainty is affected by the reduced number of points available to perform the regression.

42 ENGINEERING↗

NASA NTRS1994

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗

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 NTRS2018

A Blueprint for Demonstrating Quantum Supremacy with Superconducting Qubits

Long coherence times and high fidelity control recently achieved in scalable superconducting circuits paved the way for the growing number of experimental studies of many-qubit quantum coherent phenomena in these devices. Albeit full implementation of quantum error correction and fault tolerant quantum computation remains a challenge the near term pre-error correction devices could allow new fundamental experiments despite inevitable accumulation of errors. One such open question foundational for quantum computing is achieving the so called quantum supremacy, an experimental demonstration of a computational task that takes polynomial time on the quantum computer whereas the best classical algorithm would require exponential time and/or resources. It is possible to formulate such a task for a quantum computer consisting of less than a 100 qubits. The computational task we consider is to provide approximate samples from a non-trivial quantum distribution. This is a generalization for the case of superconducting circuits of ideas behind boson sampling protocol for quantum optics introduced by Arkhipov and Aaronson. In this presentation we discuss a proof-of-principle demonstration of such a sampling task on a 9-qubit chain of superconducting gmon qubits developed by Google. We discuss theoretical analysis of the driven evolution of the device resulting in output approximating samples from a uniform distribution in the Hilbert space, a quantum chaotic state. We analyze quantum chaotic characteristics of the output of the circuit and the time required to generate a sufficiently complex quantum distribution. We demonstrate that the classical simulation of the sampling output requires exponential resources by connecting the task of calculating the output amplitudes to the sign problem of the Quantum Monte Carlo method. We also discuss the detailed theoretical modeling required to achieve high fidelity control and calibration of the multi-qubit unitary evolution in the device. We use a novel cross-entropy statistical metric as a figure of merit to verify the output and calibrate the device controls. Finally, we demonstrate the statistics of the wave function amplitudes generated on the 9-gmon chain and verify the quantum chaotic nature of the generated quantum distribution. This verifies the implementation of the quantum supremacy protocol.

Kechedzhi, Kostyantyn↗

NASA NTRS1976

Monte Carlo calculations of diatomic molecule gas flows including rotational mode excitation

The direct simulation Monte Carlo method was used to solve the Boltzmann equation for flows of an internally excited nonequilibrium gas, namely, of rotationally excited homonuclear diatomic nitrogen. The semi-classical transition probability model of Itikawa was investigated for its ability to simulate flow fields far from equilibrium. The behavior of diatomic nitrogen was examined for several different nonequilibrium initial states that are subjected to uniform mean flow without boundary interactions. A sample of 1000 model molecules was observed as the gas relaxed to a steady state starting from three specified initial states. The initial states considered are: (1) complete equilibrium, (2) nonequilibrium, equipartition (all rotational energy states are assigned the mean energy level obtained at equilibrium with a Boltzmann distribution at the translational temperature), and (3) nonequipartition (the mean rotational energy is different from the equilibrium mean value with respect to the translational energy states). In all cases investigated the present model satisfactorily simulated the principal features of the relaxation effects in nonequilibrium flow of diatomic molecules.

Yoshikawa, K. K.↗

NASA NTRS1986

Monte Carlo models for the terrestrial exosphere over a solar cycle

Improved Monte Carlo exosphere simulations for the concentration and velocity distribution of hydrogen in the terrestrial exosphere were formulated for minimum, medium, and maximum solar cycle conditions. Both the classical exobase source of hot hydrogen and the plasmaspheric source were included, along with solar radiation pressure and photoionization. It was found that at solar minimum the hydrogen from the charge exchange of hot ions in the plasmasphere exceeds that from the exobase source not only for escape but for the population at geocentric distances greater than 2 earth radii. At about two earth radii, the equivalent temperature is about 50 percent greater than that of the exobase, a situation similar to that of the 'two-temperature' Venusian exosphere. Diurnal variations were calculated for altitudes of up to 30 earth radii; for all levels of solar activity, the inclusion of radiation pressure led to a semidiurnal variation at 3-8 earth radii and a high-altitude night/day asymmetry constituting the 'geotail'.

Tinsley, B. A.↗

DOE OSTI2024

Measuring signatures in photon angular spectra to distinguish nonlinear Compton scattering models

The collision of a high-energy electron beam with a laser pulse may be used to study radiation reaction and nonlinear Compton scattering among many other processes in strong-field quantum electrodynamics. Predictions from simulation and theory for these interactions rely on a number of approximations and assumptions that have not been experimentally tested. Here, experimentally measurable signatures are identified that might be able to distinguish between radiation reaction models, i.e., classical or quantum, or between the local constant field and local monochromatic approximations used to calculate the properties of the nonlinear Compton process. These signatures are considered through Monte Carlo simulations of various experimental conditions that are relevant to today's laser facilities. Potential detection schemes for measuring the signatures are proposed. We find that single-photon counting of keV photons to resolve harmonics and scintillator-based detection of MeV photons may allow us to validate nonlinear Compton scattering models and radiation reaction models respectively. This will require electron beams with divergence angles less than 2 mrad and less than 20% energy spread.

Russell, Brandon K↗

NASA NTRS2017

Evaluating the Stability of NASA's Space Launch System with Adaptive Augmenting Control

NASA's baseline Space Launch System (SLS) flight control system (FCS) design includes an adaptive augmenting control (AAC) component that modifies the attitude control system response to provide the classical gain-scheduled control architecture with additional performance and robustness. The NASA Engineering and Safety Center (NESC) teamed with the Space Launch System (SLS) Program to perform a comprehensive assessment of the stability and robustness of the FCS with AAC. This paper provides an overview of the approach, specific analysis techniques, and outcomes that were particularly relevant for the SLS Program. Multiple analysis techniques that specifically target the nonlinear AAC were commissioned as part of this assessment, which was completed outside of the Program's standard design analysis cycle. The following analyses were included, with each technique adding its own valuable insights: Lyapunov-based stability analysis, classical stability analysis with static AAC gain variations, circle criterion-based analysis of the FCS with a time-varying element, time-domain stability margin assessment, Monte Carlo simulations with expanded dispersions, and an extensive set of stressing cases. Several of the completed analyses focused on determining whether the inclusion of AAC introduced risk to the FCS, while others quantified the benefits of the adaptive augmentation.

Tannen S Vanzwieten↗

DOE OSTI2026

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

A Monte Carlo approach to the analysis of control system robustness

Stochastic robustness, a simple technique used to estimate the stability and performance robustness of linear, time-invariant systems, is described. The scalar probability of instability is introduced as a measure of stability robustness. Examples are given of stochastic performance robustness measures based on classical time-domain specifications. The relationship between stochastic robustness measures and control system design parameters is discussed. The technique is demonstrated by analyzing an LQG/LTR system designed for a flexible robot arm. It is concluded that the analysis of stochastic robustness offers a good alternative to existing robustness metrics.

Ray, Laura R.↗