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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Mercury's Na Exosphere from MESSENGER Data

MESSENGER entered orbit about Mercury on March 18, 2011. Since then, the Ultraviolet and Visible Spectrometer (UWS) channel of MESSENGER's Mercury Atmospheric and Surface Composition Spectrometer (MASCS) has been observing Mercury's exosphere nearly continuously. Daily measurements of Na brightness were fitted with non-uniform exospheric models. With Monte Carlo sampling we traced the trajectories of a representative number of test particles, generally one million per run per source process, until photoionization, escape from the gravitational well, or permanent sticking at the surface removed the atom from the simulation. Atoms were assumed to partially thermally accommodate on each encounter with the surface with accommodation coefficient 0.25. Runs for different assumed source processes are run separately, scaled and co-added. Once these model results were saved onto a 3D grid, we ran lines of sight from the MESSENGER spacecraft :0 infinity using the SPICE kernels and we computed brightness integrals. Note that only particles that contribute to the measurement can be constrained with our method. Atoms and molecules produced on the nightside must escape the shadow in order to scatter light if the excitation process is resonant-light scattering, as assumed here. The aggregate distribution of Na atoms fits a 1200 K gas, with a PSD distribution, along with a hotter component. Our models constrain the hot component, assumed to be impact vaporization, to be emitted with a 2500 K Maxwellian. Most orbits show a dawnside enhancement in the hot component broadly spread over the leading hemisphere. However, on some dates there is no dawn/dusk asymmetry. The portion of the hot/cold source appears to be highly variable.

Killen, Rosemary M.↗
NASA NTRS2013

Uncertainty Quantification in Remaining Useful Life of Aerospace Components using State Space Models and Inverse FORM

This paper investigates the use of the inverse first-order reliability method (inverse- FORM) to quantify the uncertainty in the remaining useful life (RUL) of aerospace components. The prediction of remaining useful life is an integral part of system health prognosis, and directly helps in online health monitoring and decision-making. However, the prediction of remaining useful life is affected by several sources of uncertainty, and therefore it is necessary to quantify the uncertainty in the remaining useful life prediction. While system parameter uncertainty and physical variability can be easily included in inverse-FORM, this paper extends the methodology to include: (1) future loading uncertainty, (2) process noise; and (3) uncertainty in the state estimate. The inverse-FORM method has been used in this paper to (1) quickly obtain probability bounds on the remaining useful life prediction; and (2) calculate the entire probability distribution of remaining useful life prediction, and the results are verified against Monte Carlo sampling. The proposed methodology is illustrated using a numerical example.

Sankararaman, Shankar↗
NASA NTRS2013

A Green's Function Approach to Simulate DNA Damage by the Indirect Effect

The DNA damage is of fundamental importance in the understanding of the effects of ionizing radiation. DNA is damaged by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research, many questions on the DNA damage by ionizing radiation remains. In the recent years, the Green's functions of the diffusion equation (GFDE) have been used extensively in biochemistry [1], notably to simulate biochemical networks in time and space [2]. In our future work on DNA damage, we wish to use an approach based on the GFDE to refine existing models on the indirect effect of ionizing radiation on DNA. To do so, we will use the code RITRACKS [3] developed at the NASA Johnson Space Center to simulate the radiation track structure and calculate the position of radiolytic species after irradiation. We have also recently developed an efficient Monte‐Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state [4], which can be modified and adapted to simulate DNA damage by free radicals. To do so, we will use the known reaction rate constants between radicals (OH, eaq, H,...) and the DNA bases, sugars and phosphates and use the sampling algorithms to simulate the diffusion of free radicals and chemical reactions with DNA. These techniques should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double‐strand breaks.

Plante, Ianik↗
NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗
NASA NTRS2010

Failure Bounding And Sensitivity Analysis Applied To Monte Carlo Entry, Descent, And Landing Simulations

In the study of entry, descent, and landing, Monte Carlo sampling methods are often employed to study the uncertainty in the designed trajectory. The large number of uncertain inputs and outputs, coupled with complicated non-linear models, can make interpretation of the results difficult. Three methods that provide statistical insights are applied to an entry, descent, and landing simulation. The advantages and disadvantages of each method are discussed in terms of the insights gained versus the computational cost. The first method investigated was failure domain bounding which aims to reduce the computational cost of assessing the failure probability. Next a variance-based sensitivity analysis was studied for the ability to identify which input variable uncertainty has the greatest impact on the uncertainty of an output. Finally, probabilistic sensitivity analysis is used to calculate certain sensitivities at a reduced computational cost. These methods produce valuable information that identifies critical mission parameters and needs for new technology, but generally at a significant computational cost.

Gaebler, John A.↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗
NASA NTRSDate not supplied

Application of a Bayesian Framework for Plasticity Model Selection

Interpretable Machine Learning (IML) has performed well when tasked with deriving constitutive material models. However, IML has been shown to prefer models that overfit noise in data, which tends to lead to bloat and a decrease in interpretability. Due to these issues, the ability of IML to reliably derive models that fit the data and are both interpretable and generalizable is limited. A method developed recently has shown promise to improve upon traditional IML by using a Bayesian fitness definition for the evolution of free-form models with non-deterministic parameters. This framework was developed for genetic-programming-based symbolic regression(GPSR) and involves model parameter estimation using Sequential Monte Carlo sampling (SMC).The method has demonstrated a reduction in bloat when dealing with noisy data in comparison to conventional GPSR. The results of this framework applied to stress-strain data for copper show models that more effectively predict the experimental data better than was previously shown with GPSR.

plasticity↗
NASA NTRS2023

Land Surfaces at the Tipping-Point for Water and Energy Balance Coupling

The surface water and energy balances can be coupled or uncoupled depending on whether the evaporation regime is water-limited or energy-limited. As the landscape loses soil moisture during drydowns, a transition between the regimes may occur, which signifies a nonlinear change in water-energy-carbon coupling. Regions that switch often between these two regimes, that is, are dominated by neither regime, are particularly vulnerable to climate variability and change. To robustly identify these tipping points, we identify drydown events based on global soil moisture data sets from remote sensing. The event identification does not rely on precipitation information and is robust with respect to measurement noise. Then, the soil moisture thresholds delineating the evaporation regime transitions are determined by Sequential Monte Carlo Sampling and a two-stage parametrization strategy. Based on the estimated soil moisture thresholds across the globe, we estimate observation-based water availability indices which quantify the nonlinear controls of soil moisture on evaporation. This framework is tested and applied globally using Soil Moisture Active Passive soil moisture retrievals. Combined with a new tippling-point metric that describes the frequency of evaporation regime transitions, we identify regions that switch often between different evaporation regimes at the global scale. Given unit shifts in soil moisture, these regions will experience the most change in how their surface water and energy are coupled.

Jianzhi Dong↗
NASA NTRSDate not supplied

L2-Charged Particle Environment (L2-CPE)Low Energy Radiation Fluence Model

The L2 Charged Particle Environment (L2-CPE) model provides estimates of number flux and fluence for the low energy electron, proton, and alpha particle populations in the near Earth solar wind and the Earth’s distant magnetosheath and magnetotail. The model is an engineering tool for specifying radiation environments over an energy range from a few eV to a few MeV of importance to surface dose and radiation damage to thin space exposed materials and is intended for use in space system design applications. Mission fluences are obtained by simulating a spacecraft flight trajectory through time-dependent bow shock and magnetopause boundaries with dimensions and orientations driven by solar wind parameters. Monte Carlo sampling of flux environments within individual plasma regimes and/or fluence accumulated along a fight trajectory through multiple regimes is used to determine statistical variations (means and extremes) of the differential number flux and fluence environments for each of the three charged particle species. Model output is differential (in energy) number fluence and flux for the three particle species in units of particles/cm2-keV and particles/cm2-sec-keV, respectively. Users can select the energy range of interest but the current version of the code (L2-CPE Version 1.4.2d) is limited to an energy range of 1 eV to 10 MeV. Flux is integrated over angle to give flux and fluence to surfaces in the ±XGSE, ±YGSE, and/or ±ZGSE directions. Figure 1 shows the opening screen from the L2-CPE graphical user interface (GUI) with options for flux, fluence, plotting model output. The flux scene generate is not currently implemented in the GUI.

Joseph I Minow↗
NASA NTRS2025

Assessment of Model Outcomes Between the Integrated Medical Model (IMM) and the Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT)

The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) is a computational model that provides human health and medical risk predictions for crewed spaceflight missions. MEDPRAT utilizes discrete event modeling and dynamic probabilistic simulation to predict critical mission outcomes (total medical events, crew health index, quality time lost, loss of crew life, removal to definitive care), condition occurrences, and resource consumption. Input parameters for MEDPRAT include crew attributes (e.g., sex), types of mission activities (e.g., whether and where crew members perform an extravehicular activity (EVA)), available resources, treatment information, and probability distributions for medical conditions. As an evolution of the Integrated Medical Model (IMM), MEDPRAT provides enhanced capabilities and higher fidelity, and incorporates more appropriate assumptions for long-duration spaceflight. IMM is the currently accepted standard for quantifying spaceflight mission medical risk in NASA operations that uses a probabilistic risk assessment (PRA) approach. MEDPRAT builds on the same logical foundation as IMM but implements the model architecture through highly optimized Monte Carlo sampling methods. An analysis is performed comparing the outputs from IMM with those from MEDPRAT V1.0 and V2.0 for the same reference missions in order to quantify similarities and differences in the model outcomes. The juxtaposition between IMM and MEDPRAT V1.0 and 2.0 shown in this report demonstrates that these two models generate very similar results; where differences in outcomes are shown, these are in accordance with the underlying assumptions and differences in the model architectures. This validation effort further establishes the credibility and reliability of the MEDPRAT software.

Matthew T Prelich↗
NASA NTRSDate not supplied

Quantifying the Sensitivity of Condition Incidence Parameters in the Evidence Library

One approach to quantifying spaceflight risk at NASA makes use event driven probabilistic techniques. The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) is such a tool that estimates medical risk metrics via simulation and enables optimization of medical resources subject to mission constraints [1]. Previous analyses have informed medical set composition, exercise countermeasures, and water intake, where each analysis quantifies the risk associated with proposed variations in system design. As future mission profiles extend beyond Low-Earth Orbit (LEO) and lengthen in duration, understanding these risks and contributing factors is critical. MEDPRAT employs Monte Carlo sampling techniques to simulate missions and track the occurrence of medical events. These events follow fault-tree-like progressions through levels of severity and mitigation via medical treatment to many possible outcomes and these are reported throughout the mission. Making this possible, are the medical databases that contain evidence gathered by the Human Research Program (HRP). Quantifying the impact of uncertainty or variability in the input data is an important step in evaluating the credibility of modeling and simulation results. In this work, we investigate the sensitivity of medical risk metrics with respect to the condition incidence parameters within the Evidence Library (EL) [2] as the medical database input for MEDPRAT. The medical conditions, contained in the EL, are equipped with incidence rates that describe the likelihood that the condition will occur. These incidence rates reflect historical spaceflight data or when appropriate, terrestrial data. In this presentation, we will explore how uncertainty in these rates propagate to the medical risk described by MEDPRAT. These results identify the conditions and parameters with the largest contribution to medical risks.

Ian Lim↗
DOE OSTI2012

NUCLEAR ENERGY SYSTEM COST MODELING

The U.S. Department of Energy’s Fuel Cycle Technologies (FCT) Program is preparing to perform an evaluation of the full range of possible Nuclear Energy Systems (NES) in 2013. These include all practical combinations of fuels and transmuters (reactors and sub-critical systems) in single and multi-tier combinations of burners and breeders with no, partial, and full recycle. As part of this evaluation, Levelized Cost of Electricity at Equilibrium (LCAE) ranges for each representative system will be calculated. To facilitate the cost analyses, the 2009 Advanced Fuel Cycle Cost Basis Report is being amended to provide up-to-date cost data for each step in the fuel cycle, and a new analysis tool, NE-COST, has been developed. This paper explains the innovative “Island” approach used by NE-COST to streamline and simplify the economic analysis effort and provides examples of LCAE costs generated. The Island approach treats each transmuter (or target burner) and the associated fuel cycle facilities as a separate analysis module, allowing reuse of modules that appear frequently in the NES options list. For example, a number of options to be screened will include a once-through uranium oxide (UOX) fueled light water reactor (LWR). The UOX LWR may be standalone, or may be the first stage in a multi-stage system. Using the Island approach, the UOX LWR only needs to be modeled once and the module can then be reused on subsequent fuel cycles. NE-COST models the unit operations and life cycle costs associated with each step of the fuel cycle on each island. This includes three front-end options for supplying feedstock to fuel fabrication (mining/enrichment, reprocessing of used fuel from another island, and/or reprocessing of this island’s used fuel), along with the transmuter and back-end storage/disposal. Results of each island are combined based on the fractional energy generated by each islands in an equilibrium system. The cost analyses use the probability distributions of key parameters and employs Monte Carlo sampling to arrive at an island’s cost probability density function (PDF). When comparing two NES to determine delta cost, strongly correlated parameters can be cancelled out so that only the differences in the systems contribute to the relative cost PDFs. For example, one comparative analysis presented in the paper is a single stage LWR-UOX system versus a two-stage LWR-UOX to LWR-MOX system. In this case, the first stage of both systems is the same (but with different fractional energy generation), while the second stage of the UOX to MOX system uses the same type transmuter but the fuel type and feedstock sources are different. In this case, the cost difference between systems is driven by only the fuel cycle differences of the MOX stage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Mads.jl

Mads.jl (Model analysis and decision support in Julia) is a code that streamlines the process of using data and models for analysis and decision support. It is based on another open-source code developed at LANL and written in C/C++ (MADS; http://mads.lanl.gov; LA-CC-11- 035). Mads.jl can work with external models of arbitrary complexity as well as built-in models of flow and transport in porous media. It enables a number of data- and model-based analyses including model calibration, sensitivity analysis, uncertainty quantification, and decision analysis. The code also can use a series of alternative adaptive computational techniques for Bayesian sampling, Monte Carlo, and Bayesian Information-Gap Decision Theory. The code is implemented in the Julia programming language, and has high-performance (parallel) and memory management capabilities. The code uses a series of third party modules developed by others. The code development will also include contributions to the existing third party modules written in Julia; this contributions will be important for the efficient implementation of the algorithm used by Mads.jl. The code also uses a series of LANL developed modules that are developed by Dan O'Malley; these modules will be also a part of the Mads.jl release. Mads.jl will be released under GPL V3 license. The code will be distributed as a Git repo at gitlab.com and github.com. Mads.jl manual and documentation will be posted at madsjulia.lanl.gov.

Vesselinov, Velimir↗
DOE OSTI2014

ANALYSIS OF SAFETY IMPACTS FROM EXTERNAL FLOODING USING THE RISK-INFORMED SAFETY MARGIN CHARACTERIZATION (RISMC) TOOLKIT

The existing fleet of U.S. nuclear power plants is in the process of extending its lifetime and increasing the power generated from these plants via power uprates. In order to evaluate the impact of these factors on the safety of the plant, the Risk-Informed Safety Margin Characterization (RISMC) project aims to provide insight to decision makers through a series of simulations of the plant dynamics for different initial conditions (e.g., probabilistic analysis and uncertainty quantification). This paper demonstrates how Idaho National Laboratory (INL) researchers use the RISMC Toolkit to investigate complex nuclear plant phenomena using RAVEN and RELAP-7. The analysis focused on a highly relevant topic currently facing some nuclear power plants – specifically flooding issues. This research and development looked at challenges to a hypothetical pressurized water reactor, including: (1) a potential loss of off-site power followed by the possible loss of all diesel generators (i.e., a station black-out event), (2) earthquake induced station-blackout, and (3) a potential earthquake induced tsunami flood. The analysis is performed by using a set of codes: a thermal-hydraulic code (RELAP-7), a flooding simulation tool (NEUTRINO) and a stochastic analysis tool (RAVEN) – these are currently under development at INL. Using RAVEN, we were able to perform multiple RELAP-7 simulation runs by changing specific parts of the model in order to reflect specific aspects of different scenarios, including both the failure and recovery of critical components. The simulation employed traditional statistical tools (such as Monte-Carlo sampling) and more advanced machine-learning based algorithms to perform uncertainty quantification in order to understand changes in system performance and limitations as a consequence of power uprate. Qualitative and quantitative results obtained gave a detailed picture of the issues associated with potential accident scenarios. These types of insights can provide useful material for decision makers to perform risk-informed margins management.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

On the short-term uncertainty in performance f a point absorber wave energy converter

Of interest, in this study, is the quantification of uncertainty in the performance of a two-body wave point absorber (Reference Model 3 or RM3), which serves as a wave energy converter (WEC). We demonstrate how simulation tools may be used to establish short-term relationships between any performance parameter of the WEC device and wave height in individual sea states. We demonstrate this methodology for two sea states. Efficient structural reliability methods, validated using more expensive Monte Carlo sampling, allow the estimation of uncertainty in performance of the device. Such methods, when combined with metocean data quantifying the likelihood of different sea states, can be useful in long-term studies and in reliability-based design.

16 TIDAL AND WAVE POWER↗
DOE OSTI2017

Probabilistic Wind Power Ramp Forecasting Based on a Scenario Generation Method: Preprint

Wind power ramps (WPRs) are particularly important in the management and dispatch of wind power, and they are currently drawing the attention of balancing authorities. With the aim to reduce the impact of WPRs for power system operations, this paper develops a probabilistic ramp forecasting method based on a large number of simulated scenarios. An ensemble machine learning technique is first adopted to forecast the basic wind power forecasting scenario and calculate the historical forecasting errors. A continuous Gaussian mixture model (GMM) is used to fit the probability distribution function (PDF) of forecasting errors. The cumulative distribution function (CDF) is analytically deduced. The inverse transform method based on Monte Carlo sampling and the CDF is used to generate a massive number of forecasting error scenarios. An optimized swinging door algorithm is adopted to extract all the WPRs from the complete set of wind power forecasting scenarios. The probabilistic forecasting results of ramp duration and start time are generated based on all scenarios. Numerical simulations on publicly available wind power data show that within a predefined tolerance level, the developed probabilistic wind power ramp forecasting method is able to predict WPRs with a high level of sharpness and accuracy.

17 WIND ENERGY↗
DOE OSTI2016

Joint Image Framework

The software implements probabilistic modeling of astronomical images of stars and galaxies as seen by multiple telescopes and/or at multiple observation times. The image models are rendered with the external software package GalSim (https://github.com GalSim- developers GalSim ). JIF primarily handles the parsing of image model parameters, Monte Carlo sampling, and diagnostic plots.

Dawson, William↗
DOE OSTI2017

AN INVESTIGATION OF CONSISTENT UNCERTAINTY PROPAGATION AND SENSITIVITY ANALYSIS IN A MULTIPHYSICS ENVIRONMENT

The framework for a LOCA toolkit for the U.S. light water reactors (LOTUS) is outlined. LOTUS allows for the coupling of the core design code VERA-CS with the fuel performance code FRAPCON such that uncertainties can be propagated consistently. An analysis of inherent model discrepancies, confirms that all differences are either negligible or the focus on immediate future work. Monte Carlo sampling is performed on 15 inputs for a depletion case of one fuel assembly. Both uncertainty quantification (UQ) and sensitivity analysis (SA) are performed on the figures of merit (FOM) of minimum departure from nucleate boiling ratio (MDNBR), maximum fuel centerline temperature (MFCT), and gap conductance at peak power (GCPP). Pearson correlation coefficient, Sobol sensitivity indices, and the density-based, delta moment independent measure are used as SA metrics. All measures yield consistent rankings of inputs matching past work.

10CFR50.46c↗