Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Investigation of radiative interaction in laminar flows using Monte Carlo simulation

The Monte Carlo method (MCM) is employed to study the radiative interactions in fully developed laminar flow between two parallel plates. Taking advantage of the characteristics of easy mathematical treatment of the MCM, a general numerical procedure is developed for nongray radiative interaction. The nongray model is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. To validate the Monte Carlo simulation for nongray radiation problems, the results of radiative dissipation from the MCM are compared with two available solutions for a given temperature profile between two plates. After this validation, the MCM is employed to solve the present physical problem and results for the bulk temperature are compared with available solutions. In general, good agreement is noted and reasons for some discrepancies in certain ranges of parameters are explained.

Liu, Jiwen↗
DOE OSTI2023

Improved thermonuclear rate of 42 Ti( p , γ ) 43 V and its astrophysical implication in the rp process

Context. Accurate 42 Ti(p,γ) 43 V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rp process) that occurs in X-ray bursts. Aims. We aim to improve the thermonuclear rates of 42 Ti(p,γ) 43 V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on the rp process. Methods. We reevaluated the reaction rate of 42 Ti(p,γ) 43 V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn. Results. The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus 43 V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region of rp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 10 8 K. Using a trajectory with a peak temperature of 1.95×10 9 K, we performed the rp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that the rp-process path does not bypass the isotope 43 V. Conclusions. Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for 42 Ti(p,γ) 43 V. The adoption of the new rates confirms that the reaction path of 42 Ti(p,γ) 43 V(p,γ) 44 Cr(β + ) 44 V is a key branch of the rp process in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Uncertainty-quantification-enabled inversion of nuclear responses

Nuclear quantum many-body methods rely on integral transform techniques to infer properties of electroweak response functions from ground-state expectation values. Retrieving the energy dependence of these responses is highly nontrivial, especially for quantum Monte Carlo methods, as it requires inverting the Laplace transform, a notoriously ill-posed problem. Here, in this work, we propose an artificial neural network architecture suitable for accurate response function reconstruction with precise estimation of the uncertainty of the inversion. We demonstrate the capabilities of this new architecture benchmarking it against maximum entropy and previously developed neural network methods designed for a similar task, paying particular attention to its robustness noise in the Euclidean

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2008

Ares I Flight Control System Overview

This paper describes the control challenges posed by the Ares I vehicle, the flight control system design and performance analyses used to test and verify the design. The major challenges in developing the control system are structural dynamics, dynamic effects from the powerful first stage booster, aerodynamics, first stage separation and large uncertainties in the dynamic models for all these. Classical control techniques were employed using innovative methods for structural mode filter design and an anti-drift feature to compensate for translational and rotational disturbances. This design was coded into an integrated vehicle flight simulation and tested by Monte Carlo methods. The product of this effort is a linear, robust controller design that is easy to implement, verify and test.

Hall, Charles↗
NASA NTRS2008

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗
DOE OSTI2022

Estimating Cosmological Constraints from Galaxy Cluster Abundance using Simulation-Based Inference

Inferring the values and uncertainties of cosmological parameters in a cosmology model is of paramount importance for modern cosmic observations. In this paper, we use the simulation-based inference (SBI) approach to estimate cosmological constraints from a simplified galaxy cluster observation analysis. Using data generated from the Quijote simulation suite and analytical models, we train a machine learning algorithm to learn the probability function between cosmological parameters and the possible galaxy cluster observables. The posterior distribution of the cosmological parameters at a given observation is then obtained by sampling the predictions from the trained algorithm. Our results show that the SBI method can successfully recover the truth values of the cosmological parameters within the 2σ limit for this simplified galaxy cluster analysis, and acquires similar posterior constraints obtained with a likelihood-based Markov Chain Monte Carlo method, the current state-of the-art method used in similar cosmological studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1985

Advanced reliability methods for structural evaluation

Fast probability integration (FPI) methods, which can yield approximate solutions to such general structural reliability problems as the computation of the probabilities of complicated functions of random variables, are known to require one-tenth the computer time of Monte Carlo methods for a probability level of 0.001; lower probabilities yield even more dramatic differences. A strategy is presented in which a computer routine is run k times with selected perturbed values of the variables to obtain k solutions for a response variable Y. An approximating polynomial is fit to the k 'data' sets, and FPI methods are employed for this explicit form.

Wirsching, P. H.↗
DOE OSTI2021

25 Siβ + -decay spectroscopy

Background: β-decay spectroscopy provides valuable information on exotic nuclei and a stringent test for nuclear theories beyond the stability line. Purpose: To search for new β-delayed protons and γ rays of 25 Si to investigate the properties of 25 Al excited states. Method: 25 Si β decays were measured by using the Gaseous Detector with Germanium Tagging system at the National Superconducting Cyclotron Laboratory. The protons and γ rays emitted in the decay were detected simultaneously. Here, a Monte Carlo method was used to model the Doppler broadening of 24 Mg γ-ray lines caused by nuclear recoil from proton emission. Shell-model calculations using two newly developed universal sd-shell Hamiltonians were performed. Results: The most precise 25Si half-life to date has been determined. A new proton branch at 724(4) keV and new proton-γ-ray coincidences have been identified. Three 24 Mg γ-ray lines and eight 25 Al γ-ray lines are observed for the first time in 25 Si decay. The first measurement of the 25 Si β-delayed γ-ray intensities through the 25 Al unbound states is reported. All the bound states of 25 Al are observed to be populated in the β decay of 25 Si. Several inconsistencies between the previous measurements have been resolved, and new information on the 25 Al level scheme is provided. An enhanced decay scheme has been constructed and compared to the mirror decay of 25 Na and the shell-model calculations. Conclusions: The measured excitation energies, γ-ray and proton branchings, log $ft$ values, and Gamow-Teller transition strengths for the states of 25 Al populated in the β decay of 25 Si are in good agreement with the shell model calculations, offering gratifyingly consistent insights into the fine nuclear structure of 25 Al.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neural simulation-based inference of the neutron star equation of state directly from telescope spectra

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A continuous symmetry breaking measure for finite clusters using Jensen-Shannon divergence

A quantitative measure of symmetry breaking is introduced that allows the quantification of which symmetries are most strongly broken due to the introduction of some kind of defect in a perfect structure. The method uses a statistical approach based on the Jensen-Shannon divergence. The measure is calculated by comparing the transformed atomic density function with its original. Software code is presented that carries the calculations out numerically using Monte Carlo methods. The behavior of this symmetry breaking measure is tested for various cases including finite size crystallites (where the surfaces break the crystallographic symmetry), atomic displacements from high symmetry positions, and collective motions of atoms due to rotations of rigid octahedra. Finally, the approach provides a powerful tool for assessing local symmetry breaking and offers new insights that can help researchers understand how different structural distortions affect different symmetry operations.

atomic & molecular structure↗
NASA NTRSDate not supplied

Monocular Horizon Navigation

The Artemis program advances towards operations on the lunar surface, where precise surface localization is a driving need for safety and science. Using the observable horizon as an image landmark allows for estimating the photographer’s position. This work analyzes the application of Perspective-n-Point (PnP) algorithms to this lunar localization problem. Batch simulations using lunar topography display the accuracy and drawbacks of these methods. Monte Carlo techniques show the pipeline’s solutions as measurements provided to a navigation filter. When filtered, these solutions have position errors of 50 meters or better, which is similar to or better than the performance of other methods of surface localization.

optical navigation↗
DOE OSTI2021

Probabilistic Assessment and Uncertainty Analysis of CO2 Storage Capacity of the Morrow B Sandstone—Farnsworth Field Unit

This paper presents probabilistic methods to estimate the quantity of carbon dioxide (CO2) that can be stored in a mature oil reservoir and analyzes the uncertainties associated with the estimation. This work uses data from the Farnsworth Field Unit (FWU), Ochiltree County, Texas, which is currently undergoing a tertiary recovery process. The input parameters are determined from seismic, core, and fluid analyses. The results of the estimation of the CO2 storage capacity of the reservoir are presented with both expectation curve and log probability plot. The expectation curve provides a range of possible outcomes such as the P90, P50, and P10. The deterministic value is calculated as the statistical mean of the storage capacity. The coefficient of variation and the uncertainty index, P10/P90, is used to analyze the overall uncertainty of the estimations. A relative impact plot is developed to analyze the sensitivity of the input parameters towards the total uncertainty and compared with Monte Carlo. In comparison to the Monte Carlo method, the results are practically the same. The probabilistic technique presented in this paper can be applied in different geological settings as well as other engineering applications.

03 NATURAL GAS↗
NASA NTRS1972

A modified Monte Carlo model for the ionospheric heating rates

A Monte Carlo method is adopted as a basis for the derivation of the photoelectron heat input into the ionospheric plasma. This approach is modified in an attempt to minimize the computation time. The heat input distributions are computed for arbitrarily small source elements that are spaced at distances apart corresponding to the photoelectron dissipation range. By means of a nonlinear interpolation procedure their individual heating rate distributions are utilized to produce synthetic ones that fill the gaps between the Monte Carlo generated distributions. By varying these gaps and the corresponding number of Monte Carlo runs the accuracy of the results is tested to verify the validity of this procedure. It is concluded that this model can reduce the computation time by more than a factor of three, thus improving the feasibility of including Monte Carlo calculations in self-consistent ionosphere models.

Mayr, H. G.↗
DOE OSTI2024

CORRLA-RS

The CORRLA-RS package provides a suite of statistical methods for sampling multidimensional distributions and to conduct sensitivity and correlation analysis of large scale data in the Rust programming language. The software provides a unique solution to multidimensional constrained sampling problems utilizing a combination of parallelized Markov Chain Monte Carlo methods and traditional rejection sampling. The sensitivity and correlation analysis methods are backed by a high performance randomized singular value decomposition implementation which enables datasets larger than the random access memory (RAM) size to be analyzed. Additionally, CORRLA-RS implements the active subspace identification method using a KD-Tree and the randomized singular value decomposition acting in concert.

Gurecky, William [Oak Ridge National Laboratory (O↗
NASA NTRS1985

Comments on the measurements of multiple muon phenomena

The extensive air showers in the energy around 10 to the 15th power eV include those initiated by astrophysical primary gamma-rays. The observations need a precise measurement on the directions of primary particles. It is one of the methods to measure the directions of high-energy muons in air showers. The accuracy in measuring the direction, by calculating the cosmic-ray phenomena in the atmosphere at very high energy was investgated. The results calculated by Monte Carlo method suggest that one may determine the direction of primary cosmic-rays within errors of 10/3 rad in observing muons of above 100 GeV at sea level.

Sato, T.↗
DOE OSTI2024

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗