Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Sampling”

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 289 records · Page 16

Spatially-dependent Hybridization of Potentials in LAMMPS

The purpose of the project was to facilitate joint Sandia/ExxonMobil collaborative research activities of mutual interest in the area of atomistic simulation of materials behavior using Sandia's LAMMPS molecular dynarnics software package. Work began in late 2013 and ended in early 2017. Initially, the main focus was on the development of a spatially-dependent hybridization capability for interatomic potentials in LAMMPS. Later, attention shifted to modeling and simulation of phase transformation kinetics in iron. Other research activities involved a variety of LAMMPS functionality, including interatomic potentials, Monte Carlo sampling, and energy minimization.

36 MATERIALS SCIENCE↗

The Quantum Alternating Operator Ansatz on Maximum k -Vertex Cover

The Quantum Alternating Operator Ansatz is a generalization of the Quantum Approximate Optimization Algorithm (QAOA) designed for finding approximate solutions to combinatorial optimization problems with hard constraints. In this paper, we study Maximum k-Vertex Cover under this ansatz due to its modest complexity, while still being more complex than the well studied problems of Max-Cut and Max E3-LIN2. Our approach includes (i) a performance comparison between easy-to-prepare classical states and Dicke states as starting states, (ii) a performance comparison between two XY -Hamiltonian mixing operators: the ring mixer and the complete graph mixer, (iii) an analysis of the distribution of solutions via Monte Carlo sampling, and (iv) the exploration of efficient angle selection strategies.

97 MATHEMATICS AND COMPUTING↗

Mapping Stochastic Devices to Probabilistic Algorithms

Probabilistic and Bayesian neural networks have long been proposed as a method to incorporate uncertainty about the world (both in training data and operation) into artificial intelligence applications. One approach to making a neural network probabilistic is to leverage a Monte Carlo sampling approach that samples a trained network while incorporating noise. Such sampling approaches for neural networks have not been extensively studied due to the prohibitive requirement of many computationally expensive samples. While the development of future microelectronics platforms that make this sampling more efficient is an attractive option, it has not been immediately clear how to sample a neural network and what the quality of random number generation should be. This research aimed to start addressing these two fundamental questions by examining basic “off the shelf” neural networks can be sampled through a few different mechanisms (including synapse “dropout” and neuron “dropout”) and examine how these sampling approaches can be evaluated both in terms of evaluating algorithm effectiveness and the required quality of random numbers.

97 MATHEMATICS AND COMPUTING↗

Posterior Regularized Bayesian Neural Network

Traditional NNs often lack the ability for uncertainty quantification. Bayesian NNs(BNNs) could help measure the confidence level by using distributions in NNs modeling. Besides, knowledge is commonly available and could improve the performance of BNNs if it can be properly incorporated. In this work, we propose a novel Posterior-Regularized BNN(PR-BNN) model by incorporating soft and hard constraints as a posterior regularization term. We also propose an augmented Lagrangian method and stochastic optimization algorithm for efficient updating via Monte Carlo sampling. The simulations and case studies for solar PV plants have shown the performance improvement of the proposed model over traditional BNNs.

97 MATHEMATICS AND COMPUTING↗

Improving the Parameterization of Cloud and Rain Microphysics in E3SM using Novel Observationally-Constrained Bayesian Approach (Final Technical Report)

In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.

54 ENVIRONMENTAL SCIENCES↗

LandScan Mosaic

The LandScan program at Oak Ridge National Laboratory (ORNL), in collaboration with the National Geospatial-Intelligence Agency (NGA), continues to deliver the most accurate and up to date global, high resolution gridded population data. Additionally, the latest advancements in the LandScan HD methodology led to reduced latency in development of rapid updates for geopolitical events. With momentum towards reporting more up to date population estimates, feedback from the user community expressed interest in reporting population estimates in ranges - whether to express a level of uncertainty or confirm to leadership and stakeholders the modeled data are estimates. Building upon the need to understand uncertainty or confidence in the modeled data and report ranges at the global scale, LandScan Mosaic was developed. LandScan Mosaic represents the next generation of high-resolution population modeling, building upon the established success of previous LandScan HD iterations. While LandScan HD employed a deterministic big data fusion approach, LandScan Mosaic enhances this methodology by integrating advanced machine learning techniques to impute missing, yet crucial, population model parameters. This advancement allows for probabilistic modeling of building occupancy and population distribution, incorporating uncertainty quantification through Monte Carlo sampling methods. By combining big data fusion with machine learning-driven imputation and stochastic modeling, LandScan Mosaic provides a more comprehensive and robust representation of population dynamics. LandScan Mosaic will be following the in the footsteps of its longstanding counterpart LandScan Global and releasing a global gridded population raster, at the 3-arcsecond resolution. This technical report documents the current stage of development of LandScan Mosaic, detailing the methodologies and data sources behind the modeling. Stakeholders are encouraged to use this document as an authoritative reference for insight into Mosaic’s data development processes. However, readers should note that LandScan Mosaic remains in a late-stage research and development phase, and methodologies and data presented here are subject to refinements ahead of the anticipated global release in Summer 2025. Feedback and inquiries from users and stakeholders are welcomed as we continue to refine and enhance this important population resource.

97 MATHEMATICS AND COMPUTING↗

Koopman-Based Approach to Nonintrusive Reduced Order Modeling: Application to Aerodynamic Shape Optimization and Uncertainty Propagation

A methodology for nonintrusive projection-based nonlinear model reduction originally presented by Renganathan et al. is further extended here toward parametric systems with a focus on application to aerospace design. Specifically, the method is extended to address static systems with parametric geometry (that deforms the mesh) in addition to parametric freestream boundary conditions. The main idea is to first perform a transformation on the governing equations such that it is lifted to a higher-dimensional but linear underdetermined system. This enables one to extract the system matrices easily as compared to that of the original nonlinear system. The underdetermined system is closed with a set of model-dependent nonlinear constraints upon which the model reduction is finally performed. The methodology is validated on the subsonic and transonic inviscid flows past the NACA0012 and the RAE2822 airfoils with parametrized shapes. The utility of the approach is further demonstrated by applying it to two common problems in aerospace design, namely, derivative-free global optimization and parametric uncertainty quantification with Monte Carlo sampling. Overall, the methodology is shown to achieve accuracy up to 5% and computational speedup of two to three orders of magnitude relative to the full order model. Comparison against another nonintrusive model reduction method revealed that the proposed approach is more robust, accurate, and retains the consistency between the state variables.

97 MATHEMATICS AND COMPUTING↗

Statistical aspects of carbon fiber risk assessment modeling

The probabilistic and statistical aspects of the carbon fiber risk assessment modeling of fire accidents involving commercial aircraft are examined. Three major sources of uncertainty in the modeling effort are identified. These are: (1) imprecise knowledge in establishing the model; (2) parameter estimation; and (3)Monte Carlo sampling error. All three sources of uncertainty are treated and statistical procedures are utilized and/or developed to control them wherever possible.

Gross, D.↗

Nonlinear filtering for spacecraft attitude estimation

Nonlinear filtering techniques are applied to spacecraft attitude estimation using quaternion parameterization for the attitude kinematics. By replacing the angular velocity vector by the gyro output vector, a state dependent noise vector is introduced in the seven-dimensional system equations. The resulting conditional probability density function from the Ito differential rule is governed by the Fokker Planck partial differential equation which is approximated by the second order mean and covariance differential equations. In order to minimize computer loading, the covariance propagation is carried out in six-dimensional state space using a matrix transformation. The star tracker data is used to update the covariance matrix in the seven-dimensional space. The algorithm is simulated for an earth pointing spacecraft mission, using Monte Carlo samples of gyro and star measurements. The performance of the second order filter is compared with the extended Kalman Filter through several simulation runs and drift rates have been identified.

Vathsal, S.↗

Monte Carlo models and analysis of galactic disk gamma-ray burst distributions

Gamma-ray bursts are transient astronomical phenomena which have no quiescent counterparts in any region of the electromagnetic spectrum. Although temporal and spectral properties indicate that these events are likely energetic, their unknown spatial distribution complicates astrophysical interpretation. Monte Carlo samples of gamma-ray burst sources are created which belong to Galactic disk populations. Spatial analysis techniques are used to compare these samples to the observed distribution. From this, both quantitative and qualitative conclusions are drawn concerning allowed luminosity and spatial distributions of the actual sample. Although the Burst and Transient Source Experiment (BATSE) experiment on Gamma Ray Observatory (GRO) will significantly improve knowledge of the gamma-ray burst source spatial characteristics within only a few months of launch, the analysis techniques described herein will not be superceded. Rather, they may be used with BATSE results to obtain detailed information about both the luminosity and spatial distributions of the sources.

Hakkila, Jon↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

Space Radiation Cancer Risk Projections for Exploration Missions: Uncertainty Reduction and Mitigation

In this paper we discuss expected lifetime excess cancer risks for astronauts returning from exploration class missions. For the first time we make a quantitative assessment of uncertainties in cancer risk projections for space radiation exposures. Late effects from the high charge and energy (HZE) ions present in the galactic cosmic rays including cancer and the poorly understood risks to the central nervous system constitute the major risks. Methods used to project risk in low Earth orbit are seen as highly uncertain for projecting risks on exploration missions because of the limited radiobiology data available for estimating HZE ion risks. Cancer risk projections are described as a product of many biological and physical factors, each of which has a differential range of uncertainty due to lack of data and knowledge. Monte-Carlo sampling from subjective error distributions represents the lack of knowledge in each factor to quantify risk projection overall uncertainty. Cancer risk analysis is applied to several exploration mission scenarios. At solar minimum, the number of days in space where career risk of less than the limiting 3% excess cancer mortality can be assured at a 95% confidence level is found to be only of the order of 100 days.

Cucinotta, Francis↗

Space Radiation Cancer Risks and Uncertainties for Mars Missions

Projecting cancer risks from exposure to space radiation is highly uncertain because of the absence of data for humans and because of the limited radiobiology data available for estimating late effects from the high-energy and charge (HZE) ions present in the galactic cosmic rays (GCR). Cancer risk projections involve many biological and physical factors, each of which has a differential range of uncertainty due to the lack of data and knowledge. We discuss an uncertainty assessment within the linear-additivity model using the approach of Monte Carlo sampling from subjective error distributions that represent the lack of knowledge in each factor to quantify the overall uncertainty in risk projections. Calculations are performed using the space radiation environment and transport codes for several Mars mission scenarios. This approach leads to estimates of the uncertainties in cancer risk projections of 400-600% for a Mars mission. The uncertainties in the quality factors are dominant. Using safety standards developed for low-Earth orbit, long-term space missions (>90 days) outside the Earth's magnetic field are currently unacceptable if the confidence levels in risk projections are considered. Because GCR exposures involve multiple particle or delta-ray tracks per cellular array, our results suggest that the shape of the dose response at low dose rates may be an additional uncertainty for estimating space radiation risks.

NASA Discipline Radiation Health↗

Uncertainties in Estimates of the Risks of Late Effects from Space Radiation

The health risks faced by astronauts from space radiation include cancer, cataracts, hereditary effects, and non-cancer morbidity and mortality risks related to the diseases of the old age. Methods used to project risks in low-Earth orbit are of questionable merit for exploration missions because of the limited radiobiology data and knowledge of galactic cosmic ray (GCR) heavy ions, which causes estimates of the risk of late effects to be highly uncertain. Risk projections involve a product of many biological and physical factors, each of which has a differential range of uncertainty due to lack of data and knowledge. Within the linear-additivity model, we use Monte-Carlo sampling from subjective uncertainty distributions in each factor to obtain a Maximum Likelihood estimate of the overall uncertainty in risk projections. The resulting methodology is applied to several human space exploration mission scenarios including ISS, lunar station, deep space outpost, and Mar's missions of duration of 360, 660, and 1000 days. The major results are the quantification of the uncertainties in current risk estimates, the identification of factors that dominate risk projection uncertainties, and the development of a method to quantify candidate approaches to reduce uncertainties or mitigate risks. The large uncertainties in GCR risk projections lead to probability distributions of risk that mask any potential risk reduction using the "optimization" of shielding materials or configurations. In contrast, the design of shielding optimization approaches for solar particle events and trapped protons can be made at this time, and promising technologies can be shown to have merit using our approach. The methods used also make it possible to express risk management objectives in terms of quantitative objective's, i.e., the number of days in space without exceeding a given risk level within well defined confidence limits.

Cucinotta, F. A.↗

Radiation Risk Projections for Space Travel

Space travelers are exposed to solar and galactic cosmic rays comprised of protons and heavy ions moving with velocities close to the speed of light. Cosmic ray heavy ions are known to produce more severe types of biomolecular damage in comparison to terrestrial forms of radiation, however the relationship between such damage and disease has not been fully elucidated. On Earth, we are protected from cosmic rays by atmospheric and magnetic shielding, and only the remnants of cosmic rays in the form of ground level muons and other secondary radiations are present. Because human epidemiology data is lacking for cosmic rays, risk projection must rely on theoretical understanding and data from experimental models exposed to space radiation using charged particle accelerators to simulate space radiation. Although the risks of cancer and other late effects from cosmic rays are currently believed to present a severe challenge to space travel, this challenge is centered on our lack of confidence in risk projections methodologies. We review biophysics and radiobiology data on the effects of the cosmic ray heavy ions, and the current methods used to project radiation risks . Cancer risk projections are described as a product of many biological and physical factors, each of which has a differential range of uncertainty due to lack of data and knowledge. Risk projections for space travel are described using Monte-Carlo sampling from subjective error di stributions that represent the lack of knowledge in each factor that contributes to the projection model in order to quantify the overall uncertainty in risk projections. This analysis is applied to space mi ssion scenarios including lunar colony, deep space outpost, and a Mars mission. Results suggest that the number of days in space where cancer mortality risks can be assured at a 95% confidence level to be below the maximum acceptable risk for radi ation workers on Earth or the International Space Station is only on the order of 100-200 days. Approaches to reduce these unceI1ainties and mitigate risks are described.

Cucinotta, Francis↗