Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty propagation”

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 559 records · Page 31

Solar activity and coronal mass ejections on the western hemisphere of the Sun in mid-August 1989: Association with interplanetary observations at the ICE and IMP 8 spacecraft

During the 10-day period from August 12 to 21, 1989, a sequence of coronal mass ejections (CMEs) was observed above the west limb of the Sun by the Solar Maximum Mission (SMM) coronagraph. Most of these CMEs apparently originated in the vicinity of one particularly active region during its passage from near central meridian to behind the west limb of the Sun. We present observations made at 1 AU during this period by the International Cometary Explorer (ICE) (formerly International Sun Earth Explorer-3 (ISEE 3)) and Interplanetary Monitoring Platform (IMP 8) spacecraft which were separated by approximately 75 deg in heliolongitude. Following CMEs on August 12 associated with solar events at approximately W40 deg, IMP 8 (in Earth orbit) detected a strong shock followed by signatures in magnetic field, solar wind plasma, and energetic ion data which suggest that CME-related material ('ejecta') forming the shock driver engulfed the spacecraft. This spacecraft only observed weak shocks, and no ejecta, from later CMEs originating further west of the spacecraft. In contrast ICE, off the west limb at approximately W75 deg, observed the shock from the W40 deg event but failed to encounter the shock driver, whereas clear ejecta signatures were observed following events further west, closer to the spacecraft heliolongitude. The disappearance of these signatures (which include bidirectional energetic ion flows, bidirectional solar wind heat fluxes, quiet, enhanced magnetic fields and anomalously cool plasma) at IMP 8 and their emergence at ICE as the solar source region moved westward supports the association of such signatures with ejecta related to CMEs. The dual-spacecraft observations are also consistent with the conclusion of Richardson and Cane (1993) that ejecta at 1 AU typically extend approximately 50 degs in longitude from the solar source. Some plausible associations between particular intervals of ejecta signatures at ICE and individual CMEs are made. However, these associations are complicated by the large number of CMEs present, by intermittent ICE data coverage, and by uncertainties in the ejecta propagation speeds to the spacecraft.

Richardson, I. G.↗

Formal results regarding metric space techniques for the study of astrophysical maps

We extend a newly developed formal system for the description of astrophyscial maps. In this formalism, we consider the difference between maps to be the distance between elements of a pseudometric space (the space of all such maps). This Ansatz allows us to measure quantitatively the difference between any two maps and to order the space of all maps. For each physical characteristic of interest, this technique assigns an 'output' function to each map; the difference between the maps is then determined from the difference between their corresponding output functions. In this present study, we show that the results of this procedure are invariant under a class of transformations of the maps and the domains of the maps. In addition, we study the propagation of errors (observational uncertainties) through this formalism. We show that the uncertainties in the output functions can be controlled provided that the signal to noise ratios in the original astrophysical maps are sufficiently high. The results of this paper thus increase the effectiveness of this formal system for the description, classification, and analysis of astrophysical maps.

Adams, Fred C.↗

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

Spacecraft Charging Sensitivity to Material Properties

Evaluating spacecraft charging behavior of a vehicle in the space environment requires knowledge of the material properties relevant to the charging process. Implementing surface and internal charging models requires a user to specify a number of material electrical properties including electrical resistivity parameters (dark and radiation induced), dielectric constant, secondary electron yields, photoemission yields, and breakdown strength in order to correctly evaluate the electric discharge threat posed by the increasing electric fields generated by the accumulating charge density. In addition, bulk material mass density and/or chemical composition must be known in order to analyze radiation shielding properties when evaluating internal charging. We will first describe the physics of spacecraft charging and show how uncertainties in material properties propagate through spacecraft charging algorithms to impact the results obtained from charging models. We then provide examples using spacecraft charging codes to demonstrate their sensitivity to material properties. The goal of this presentation is to emphasize the importance in having good information on relevant material properties in order to best characterize on orbit charging threats.

Minow, Joseph I.↗

Stardust Interstellar Preliminary Examination X: Impact Speeds and Directions of Interstellar Grains on the Stardust Dust Collector

On the basis of an interstellar dust model compatible with Ulysses and Galileo observations, we calculate and predict the trajectories of interstellar dust (ISD) in the solar system and the distribution of the impact speeds, directions, and flux of ISD particles on the Stardust Interstellar Dust Collector during the two collection periods of the mission. We find that the expected impact velocities are generally low (less than 10 km per second) for particles with the ratio of the solar radiation pressure force to the solar gravitational force beta greater than 1, and that some of the particles will impact on the cometary side of the collector. If we assume astronomical silicates for particle material and a density of 2 grams per cubic centimeter, and use the Ulysses measurements and the ISD trajectory simulations, we conclude that the total number of (detectable) captured ISD particles may be on the order of 50. In companion papers in this volume, we report the discovery of three interstellar dust candidates in the Stardust aerogel tiles. The impact directions and speeds of these candidates are consistent with those calculated from our ISD propagation model, within the uncertainties of the model and of the observations.

interstellar↗

Model-Based Fatigue Prognosis of Fiber-Reinforced Laminates Exhibiting Concurrent Damage Mechanisms

Prognostics of large composite structures is a topic of increasing interest in the field of structural health monitoring for aerospace, civil, and mechanical systems. Along with recent advancements in real-time structural health data acquisition and processing for damage detection and characterization, model-based stochastic methods for life prediction are showing promising results in the literature. Among various model-based approaches, particle-filtering algorithms are particularly capable in coping with uncertainties associated with the process. These include uncertainties about information on the damage extent and the inherent uncertainties of the damage propagation process. Some efforts have shown successful applications of particle filtering-based frameworks for predicting the matrix crack evolution and structural stiffness degradation caused by repetitive fatigue loads. Effects of other damage modes such as delamination, however, are not incorporated in these works. It is well established that delamination and matrix cracks not only co-exist in most laminate structures during the fatigue degradation process but also affect each other's progression. Furthermore, delamination significantly alters the stress-state in the laminates and accelerates the material degradation leading to catastrophic failure. Therefore, the work presented herein proposes a particle filtering-based framework for predicting a structure's remaining useful life with consideration of multiple co-existing damage-mechanisms. The framework uses an energy-based model from the composite modeling literature. The multiple damage-mode model has been shown to suitably estimate the energy release rate of cross-ply laminates as affected by matrix cracks and delamination modes. The model is also able to estimate the reduction in stiffness of the damaged laminate. This information is then used in the algorithms for life prediction capabilities. First, a brief summary of the energy-based damage model is provided. Then, the paper describes how the model is embedded within the prognostic framework and how the prognostics performance is assessed using observations from run-to-failure experiments

Prognostics↗

Precomputing Process Noise Covariance for Onboard Sequential Filters

Process noise is often used in estimation filters to account for unmodeled and mismodeled accelerations in the dynamics. The process noise covariance acts to inflate the state covariance over propagation intervals, increasing the uncertainty in the state. In scenarios where the acceleration errors change significantly over time, the standard process noise covariance approach can fail to provide effective representation of the state and its uncertainty. Consider covariance analysis techniques provide a method to precompute a process noise covariance profile along a reference trajectory using known model parameter uncertainties. The process noise covariance profile allows significantly improved state estimation and uncertainty representation over the traditional formulation. As a result, estimation performance on par with the consider filter is achieved for trajectories near the reference trajectory without the additional computational cost of the consider filter. The new formulation also has the potential to significantly reduce the trial-and-error tuning currently required of navigation analysts. A linear estimation problem as described in several previous consider covariance analysis studies is used to demonstrate the effectiveness of the precomputed process noise covariance, as well as a nonlinear descent scenario at the asteroid Bennu with optical navigation.

onboard↗

Real-Time Prediction of Safety Margins in the National Airspace

Underlying all operations in the National Airspace System (NAS) is the concept of safety. Safety, as defined by acceptable levels of risk, is to be maintained at all times. The real-time safety monitoring (RTSM) framework is under development to provide an automated system to quantify safety in the NAS, estimate the current level of safety, and predict the future evolution of safety and the occurrence of events that pose an increased risk to flights so that these occurrences can be managed strategically rather than mitigated reactively. This paper presents the mathematical framework, the models, and the monitoring and prediction algorithms used to achieve this. RTSM computes safety as expressed through a set of safety margins based on user-defined safety metrics, thresolds, and events. Sources of uncertainty are modeled and propagated through the predictions in order to compute the probabilistic evolution of safety and the probability of events that introduce increased risk to operations. A prototype implementation is discussed and results demonstrating feasibility are presented. The results highlight the kinds of predictions that can be computed and the fidelity that is currently achieved.

national airspace system↗

Monte Carlo simulations of neutrino and charged lepton propagation in the Earth with nuPyProp

An accurate modeling of neutrino flux attenuation and the distribution of leptons they produce in transit through the Earth is an essential component to determine neutrino flux sensitivities of underground, sub-orbital and space-based detectors. Through neutrino oscillations over cosmic distances, astrophysical neutrino sources are expected to produce nearly equal fluxes of electron, muon and tau neutrinos. Of particular interest are tau neutrinos that interact in the Earth at modest slant depths to produce leptons. Some leptons emerge from the Earth and decay in the atmosphere to produce extensive air showers. Future balloon-borne and satellite-based opticalCherenkov neutrino telescopes will be sensitive to upward air showers from tau neutrino induced lepton decays. We present nuPyProp, a python code that is part of the nuSpaceSim package. nuPyProp generates look-up tables for exit probabilities and energy distributions for (see paper for proper equation) propagation in the Earth. This flexible code runs with either stochastic or continuous electromagnetic energy losses for the lepton transit through the Earth. Current neutrino cross section models and energy loss models are included along with templates for user input of other models. Results from nuPyProp are compared with other recent simulation packages for neutrino and charged lepton propagation. Sources of modeling uncertainties are described and quantified.

Yosui Akaike↗

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↗

Effect of Correlated Precision Errors on Uncertainty of a Subsonic Venturi Calibration

An uncertainty analysis performed in conjunction with the calibration of a subsonic venturi for use in a turbine test facility produced some unanticipated results that may have a significant impact in a variety of test situations. Precision uncertainty estimates using the preferred propagation techniques in the applicable American National Standards Institute/American Society of Mechanical Engineers standards were an order of magnitude larger than precision uncertainty estimates calculated directly from a sample of results (discharge coefficient) obtained at the same experimental set point. The differences were attributable to the effect of correlated precision errors, which previously have been considered negligible. An analysis explaining this phenomenon is presented. The article is not meant to document the venturi calibration, but rather to give a real example of results where correlated precision terms are important. The significance of the correlated precision terms could apply to many test situations.

Hudson, S. T.↗

Charge and energy dependence of the residence time of cosmic ray nuclei below 15 GeV/nucleon

The relative abundance of nuclear species measured in cosmic rays at Earth has often been interpreted with the simple leaky box model. For this model to be consistent an essential requirement is that the escape length does not depend on the nuclear species. The discrepancy between escape length values derived from iron secondaries and from the B/C ratio was identified by Garcia-Munoz and his co-workers using a large amount of experimental data. Ormes and Protheroe found a similar trend in the HEAO data although they questioned its significance against uncertainties. They also showed that the change in the B/C ratio values implies a decrease of the residence time of cosmic rays at low energies in conflict with the diffusive convective picture. These conclusions crucially depend on the partial cross section values and their uncertainties. Recently new accurate cross sections of key importance for propagation calculations have been measured. Their statistical uncertainties are often better than 4% and their values significantly different from those previously accepted. Here, these new cross sections are used to compare the observed B/C+O and (Sc to Cr)/Fe ratio to those predicted with the simple leaky box model.

Soutoul, A.↗

Correlation Calculations for the Russian Pu Metal Fast Experiments

Nuclear criticality experiments are often conducted in campaigns with multiple variations. These experiments reuse the same basic components, like the fuel, moderator, or positioning machines. The components have uncertainties in their geometry and composition that propagate to models of the experiments. Shared components create shared uncertainty between the $k_{eff}$ of benchmarks. The shared uncertainty is commonly quantified with a covariance, or correlation coefficient. These covariances can impact criticality safety and nuclear data validation applications. While benchmark evaluations tabulate an experiment’s uncertainty, they often lack a detailed calculation of correlations between experiments. Even some very commonly used benchmarks, like the Russian Pu Metal Fast (PMF) experiments, have missing correlations. This paper presents our approach to calculate the correlations for five of the Russian PMF experiments. The experiments share hemispherical Pu shells that induce a correlation between modeled $k_{eff}$ values. We estimated the correlations with simplified and detailed models of the experiments through linear-perturbation theory. The correlations between the experiments vary significantly between the detailed vs. simplified models. We also investigate how the correlations affect validation metrics of the experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generic Post-Closure Safety Assessment for Deep Borehole Disposal of Heat-Generating Waste - 20451

The deep borehole disposal (DBD) concept for high-activity waste consists of drilling a large-diameter borehole into crystalline basement rock to a depth of about 5,000 m, placing waste packages in the lower portion of the borehole - the waste emplacement zone - and sealing and plugging the upper portion of the borehole with a combination of bentonite, cement plugs, and sand or crushed rock ballast. The waste emplacement zone is several times deeper than typical mined repositories and is well below the typical maximum depth of fresh groundwater resources. Recent work describes a generic reference case for deep borehole disposal of Cs and Sr capsules stored at the Hanford Site. Two scenarios are described: an undisturbed (nominal) scenario, and disturbed scenario, in which a waste package is stuck adjacent to a transmissive fracture zone in crystalline basement. Simulations of the processes affecting radionuclide release and transport in the two scenarios provide the basis for quantitative post-closure performance assessment (PA) of the reference case. This paper presents a probabilistic PA of the stuck waste package scenario, i.e. uncertainties affecting radionuclide mobilization and transport are propagated through 200 realizations of a predictive model to quantify uncertainties in predicted radionuclide concentrations within the borehole and the fracture zone. Simulations predict no radionuclide transport in the borehole to distances greater than 100 m above the stuck waste package. An assumed lateral pressure gradient results in predictions of limited radionuclide transport in the fracture zone at a lateral distance of 200 m from the stuck waste package. These results complement previous results for the nominal scenario and contribute to the development of a generic safety case for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

GOES dynamic propagation of attitude

The spacecraft in the next series of Geostationary Operational Environmental Satellites (GOES-Next) are Earth pointing and have 5-year mission lifetimes. Because gyros can be depended on only for a few years of continuous use, they will be turned off during routine operations. This means attitude must, at times, be determined without benefit of gyros and, often, using only Earth sensor data. To minimize the interruption caused by dumping angular momentum, these spacecraft have been designed to reduce the environmental torque acting on them and incorporate an adjustable solar trim tab for fine adjustment. A new support requirement for GOES-Next is that of setting the solar trim tab. Optimizing its setting requires an estimate of the unbalanced torque on the spacecraft. These two requirements, determining attitude without gyros and estimating the external torque, are addressed by replacing or supplementing the gyro propagation with a dynamic one, that is, one that integrates the rigid body equations of motion. By processing quarter-orbit or longer batches, this approach takes advantage of roll-yaw coupling to observe attitude completely without Sun sensor data. Telemetered momentum wheel speeds are used as observations of the unbalanced external torques. GOES-Next provides a unique opportunity to study dynamic attitude propagation. The geosynchronous altitude and adjustable trim tab minimize the external torque and its uncertainty, making long-term dynamic propagation feasible. This paper presents the equations for dynamic propagation, an analysis of the environmental torques, and an estimate of the accuracies obtainable with the proposed method.

Markley, F. Landis↗

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-Fidelity Computational Studies of Intermediate Pressure Capacitively Coupled Plasmas

This report discussed progress made in the development of a computational model for intermediate pressure capacitively coupled plasma discharge. A new full fidelity argon chemistry model has been developed that includes electron, argon ion, three lumped states for electronically excited species, and the background ground state argon. As such the mechanism has been validated over a wide range of intermediate pressure conditions from 100 mTorr to a few Torr. Uncertainties in the electron impact cross sections and the rate coefficients for the reactions are propagated through the model to establish uncertainties that emerge in the model due the reaction mechanism and its impact on the discharge properties. The model is therefore validated in face of uncertainties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncertainty methodology for in-flight thrust determination

A methodology is proposed for the evaluation of uncertainty in the in-flight determination of aircraft thrust, which provides error traceability to a national standards laboratory, and is independent of the procedure used to calculate or measure thrust in flight, thereby yielding a consistent means for the evaluation of measurement capabilities. Attention is given to the factors of measurement error, precision, bias, uncertainty, error estimation and classification, error propagation, ground testing, and the related problems of model bias error, model precision error, and the uncertainty limit.

Adams, G. R.↗