Search NASA⌕ Search

SEARCH · Search NASA

Results for “model structure uncertainty”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification (LDRD 218317)

This report summarizes the work done under the Laboratory Directed Research and Development (LDRD) project entitled "Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification". In this project we investigated multi-fidelity strategies for fusing data from information sources of varying cost and accuracy. Most existing strategies exploit hierarchical relationships between models, for example that occur when different models are generated by refining a numerical discretization parameter. In this work we focused on encoding the relationships between information sources using directed acyclic graphs. The multi-fidelity networks can have general structure and represent a significantly greater variety of modeling relationships than recursive networks used in the current state literature. Numerical results show that a non-hierarchical multi-fidelity Monte Carlo strategy can reduce the cost of estimating uncertainty in predictions of a model of plasma expanding in a vacuum by almost two orders of magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Quantifying uncertainties in the structural response of SSME blades

To quantify the uncertainties associated with the geometry and material properties of a Space Shuttle Main Engine (SSME) turbopump blade, a computer code known as STAEBL was used. A finite element model of the blade used 80 triangular shell elements with 55 nodes and five degrees of freedom per node. The whole study was simulated on the computer and no real experiments were conducted. The structural response has been evaluated in terms of three variables which are natural frequencies, root (maximum) stress, and blade tip displacements. The results of the study indicate that only the geometric uncertainties have significant effects on the response. Uncertainties in material properties have insignificant effects.

Nagpal, Vinod K.↗
DOE OSTI2021

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

A new interpretation of luminous blue stars

A major revision of current theoretical ideas about the brightest blue stars must be made if Carson's (1976) radiative opacities are adopted in stellar models. Unlike earlier opacities, these exhibit a large 'bump' due to CNO ionization, which leads to very strong central condensation, convective instability, and pulsational instability in hot diffuse stellar envelopes. Despite a number of theoretical uncertainties, the new picture of the structure of very luminous stars is reasonably successful in accounting for a variety of previously unexplained observations. The stellar models for the phase of core hydrogen burning predict large radii and rather cool effective temperatures for O stars and a spreading out of the main-sequence band in the H-R diagram toward luminous cool supergiants for masses higher than about 20 solar masses. In massive X-ray binary systems, circular orbits and supergiant-like visual companions are expected to be quite common. Long-period variability is predicted to exist for massive blue supergiants of luminosity class Ia. The models for helium stars predict large radii and rather cool effective temperatures for Wolf-Rayet stars, as well as multimodal pulsational instability and, possibly, surface turbulence for these stars.

Stothers, R.↗
NASA NTRS1983

Sensitivity of model calculations to uncertain inputs, with an application to neutron star envelopes

A method is given for determining the sensitivity of certain types of calculations to the uncertainties in the input physics or model parameters; this method is applicable to problems that involve solutions to coupled, ordinary differential equations. In particular the sensitivity of calculations of the thermal structure of neutron star envelopes to uncertainties in the opacity and equation of state is examined. It is found that the uncertainties in the relationship between the surface and interior temperatures of a neutron star are due almost entirely to the imprecision in the values of the conductive opacity in the region where the ions form a liquid; here the conductive opacity is, for the most part, due to the scattering of electrons from ions.

Epstein, R. I.↗
NASA NTRS1996

RSRM 10 percent Scale Model Drilled Hole Plate Tests

The RSRM 10% Scaled Model under design will make use of drilled hole liners to provide mass addition along the axial length of the model. The model will have two sets of liners in use at a time. The outer most liner is a flow distribution tube, the purpose of which is to help distribute the flow evenly over each model segment. The inner most liner will simulate the propellant burning surface at a burn time of 80 seconds. This liner will replicate as closely as possible the actual geometry of the full scale RSRM at the 80 second burn time. In order to obtain the correct mass flow rate for the burn time selected, it is necessary to determine the porosity of the holes drilled in each liner and the performance of those holes. The pressure drop across the liners directly effects the uniformity of the flow in the axial direction for a given model section. It is desired to have a pressure drop across the liners which is greater than the axial pressure drop in a given section. However, the pressure drop across the liner also has a bearing on the structural soundness of the model. The performance of the model was determined analytically, but there was some uncertainty as to the value of the discharge coefficient used. This uncertainty was the impetus for these drilled hole plate tests. Experimentally obtaining the discharge coefficients for sample plates of the porosity to be used in the model would increase the fidelity of the model design. These tests were developed in order to provide the required information with the least amount of testing time and hardware.

Purinton, David C.↗
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) for matching single-input single-output (SISO) transfer function models to empirical frequency response data is described. The program was written in MATLAB and has a graphical user interface (GUI). It is interactive in that the user manually builds the transfer function model using ``elementary factors'' (gain, delay, differentiators and integrators, and first- and second-order poles and zeros) and adjusts their values with sliders or entry fields. A nonlinear optimization can also be used to determine maximum-likelihood estimates of the transfer function parameters and their associated uncertainties. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and other system identification problems. FRFit is demonstrated using example problems, including the identification of aircraft transfer functions and rational function approximations of Theodorsen's function.

Frequency response↗
NASA NTRS1993

Linear, multivariable robust control with a mu perspective

The structured singular value is a linear algebra tool developed to study a particular class of matrix perturbation problems arising in robust feedback control of multivariable systems. These perturbations are called linear fractional, and are a natural way to model many types of uncertainty in linear systems, including state-space parameter uncertainty, multiplicative and additive unmodeled dynamics uncertainty, and coprime factor and gap metric uncertainty. The structured singular value theory provides a natural extension of classical SISO robustness measures and concepts to MIMO systems. The structured singular value analysis, coupled with approximate synthesis methods, make it possible to study the tradeoff between performance and uncertainty that occurs in all feedback systems. In MIMO systems, the complexity of the spatial interactions in the loop gains make it difficult to heuristically quantify the tradeoffs that must occur. This paper examines the role played by the structured singular value (and its computable bounds) in answering these questions, as well as its role in the general robust, multivariable control analysis and design problem.

Packard, Andy↗
DOE OSTI2020

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Massively Parallel Bayesian Model Calibration and Uncertainty Quantification with Applications to Nuclear Fuels and Materials

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Fuel behavior in a nuclear reactor is governed by a complex network of mechanisms interacting with various other physics aspects in the reactor system. Any model developed to represent the fuel behavior will likely be idealized resulting in uncertainties in their predictions compared to the observed data. As such, this report was motivated by the need to identify the sources of uncertainties and quantify and propagate them through the fuel model outputs. Such quantification of uncertainties will establish a level of model trustworthiness, identify approaches to improve the model trustworthiness, and even guide optimal experiment design for maximal information gain. To accomplish the uncertainty quantification for computational models, this report has relied on the Bayesian framework which provides probabilistic treatment of models their inputs and outputs. The current state-of-the-art on performing Bayesian Uncertainty Quantification (UQ) for nuclear engineering models using High Performance Computing (HPC) resources have been reviewed. Implementation of capabilities for massively parallel Bayesian UQ in Multiphysics Object-Oriented Simulation Environment (MOOSE) is discussed. Several verification cases are discussed to verify the accuracy of the quantified uncertainties using the developed computational capabilities in MOOSE. Then, the problem of quantifying the uncertainties in TRI-Structural isOtropic (TRISO) fuel silver release is addressed. For the first time, the uncertainties arising from the TRISO Fission Gas Release (FGR) model due to model inadequacy and experimental noise are quantified. Also, the Bayesian capabilities are applied to the calibration of the MATPRO creep model, a widely used model in several fuel assessment cases. The impact of the prediction uncertainties in the MATPRO model on the fuel cladding behavior as part of the TRIBULATION assessment case (which is an integral effects case) is investigated. This report concludes with a discussion on the future work for the UQ for computational models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2021

Satellite Constraints on the Latitudinal Distribution and Temperature Sensitivity of Wetland Methane Emissions

Wetland methane (CH4) emissions comprise about one-third of the global CH4 source. The latitudinal distribution and climate sensitivity of wetland CH4 fluxes are the key determinants of the global CH4-climate feedback. However, large differences exist between bottom-up estimates, informed by ground-based flux measurements, and top-down estimates derived from spaceborne total column CH4. Despite the extensive coverage of satellite CH4 concentration observations, challenges remain with using top-down estimates to test bottom-up models, mainly because of the uncertainties in the satellite retrievals, the model representation errors, the variable prior emissions, and the confounding role of the posterior error covariance structures. Here, we use satellite-based top-down CH4 flux estimates (2010–2012) to test and refine 42 bottom-up estimates of wetland emissions that use a range of hypothesized wetland extents and process controls. Our comparison between bottom-up models and satellite-based fluxes innovatively accounts for cross-correlations and spatial uncertainties typically found in top-down inverse estimates, such that only the information from satellite observations and the atmospheric transport model is kept as a constraint. We present a satellite-constrained wetland CH4 ensemble product derived from assembling the highest-performance bottom-up models, which estimates global wetland CH4 emissions of 148 (117–189, 5th–95th percentile) Tg CH4 yr−1. We find that tropical wetland emissions contribute 72% (63%–85%) to the global wetland total. We also find that a lower-than-expected temperature sensitivity agrees better with atmospheric CH4 measurements. Overall, our approach demonstrates the potential for using satellites to quantitatively refine bottom-up wetland CH4 emission estimates, their latitudinal distributions, and their sensitivity to climate.

wetland ch4↗
DOE OSTI2026

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗
DOE OSTI2021

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to 89,90 Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of 90 Sr and 89 Sr intake to humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Application of Non-Deterministic Methods to Assess Modeling Uncertainties for Reinforced Carbon-Carbon Debris Impacts

The Space Shuttle Columbia Accident Investigation Board (CAIB) made several recommendations for improving the NASA Space Shuttle Program. An extensive experimental and analytical program has been developed to address two recommendations related to structural impact analysis. The objective of the present work is to demonstrate the application of probabilistic analysis to assess the effect of uncertainties on debris impacts on Space Shuttle Reinforced Carbon-Carbon (RCC) panels. The probabilistic analysis is used to identify the material modeling parameters controlling the uncertainty. A comparison of the finite element results with limited experimental data provided confidence that the simulations were adequately representing the global response of the material. Five input parameters were identified as significantly controlling the response.

Lyle, Karen H.↗
NASA NTRS1994

A Comparison of Metallographic Cooling Rate Methods Used in Meteorites

The primary objective of this study was to test the postulate that cooling rates acquired from metal grains in chondrites are consistent with those from iron meteorites. Both types of metal occur in some Group IAB meteorites, which are mixtures of massive metal with well-developed Widmanstatten structures and chondritic inclusions with dispersed metal grains. The grains have textures and compositions similar to chondritic metal, including negligible P. The meteorites studied show little or no sign of shock reheating and textural evidence indicates that silicates and metal were mixed before Widmanstatten patterns formed during cooling. Cooling rates were obtained by comparing measured to modeled taenite grain or lamellae dimensions and central Ni contents. Modeling entails solving diffusion equations using experimental diffusion coefficients, phase relations, and bulk or local Ni and P contents, taking into account geometry, undercooling, and impingement. There is one set of parameters for grains and another, quite different set for Widmanstatten lamellae, including a factor of 30 difference in diffusion coefficients. Yet cooling rates obtained from Widmanstatten structures and metal grains in chondritic inclusions of the same meteorite are consistent; uncertainties in the best data are +/- 10 K/Ma, equivalent to a factor of 1 +/- 0.25. This agreement implies that the data and models are correct or contain fortuitously offsetting errors, which is quite unlikely. Cooling rates range from 40 K/Ma to 70 K/Ma in IAB meteorites that contain both grains and Widmanstatten structures. Rates based on grains in Ni-poor and Ni-rich meteorites lacking Widmansatten patterns expand the range from 30 K/Ma to perhaps 200 K/Ma. Cooling rates correlate with Ni content; Ni-poor meteorites have slower rates than Ni-rich ones. Evidently, IAB meteorites were radially distributed over greater than 30km in a body with a radius less than 50km. A comparison of the available Ar ages with cooling times inferred from the cooling rates suggests that the parent body cooled more slowly after the metallographic cooling rates were established.

Herpfer, Marc A.↗
NASA NTRS2011

Control System Modeling for the Thirty Meter Telescope Primary Mirror

The Thirty Meter Telescope primary mirror is composed of 492 segments that are controlled to high precision in the presence of wind and vibration disturbances, despite the interaction with structural dynamics. The higher bandwidth and larger number of segments compared with the Keck telescopes requires greater attention to modeling to ensure success. We focus here on the development and validation of a suite of quasi-static and dynamic modeling tools required to support the design process, including robustness verification, performance estimation, and requirements flowdown. Models are used to predict the dynamic response due to wind and vibration disturbances, estimate achievable bandwidth in the presence of control-structure-interaction (CSI) and uncertainty in the interaction matrix, and simulate and analyze control algorithms and strategies, e.g. for control of focus-mode, and sensor calibration. Representative results illustrate TMT performance scaling with parameters, but the emphasis is on the modeling framework itself.

uncetainty↗