Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Inference”

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 1,495 records · Page 83

An evidential reasoning extension to quantitative model-based failure diagnosis

The detection and diagnosis of failures in physical systems characterized by continuous-time operation are studied. A quantitative diagnostic methodology has been developed that utilizes the mathematical model of the physical system. On the basis of the latter, diagnostic models are derived each of which comprises a set of orthogonal parity equations. To improve the robustness of the algorithm, several models may be used in parallel, providing potentially incomplete and/or conflicting inferences. Dempster's rule of combination is used to integrate evidence from the different models. The basic probability measures are assigned utilizing quantitative information extracted from the mathematical model and from online computation performed therewith.

Gertler, Janos J.↗

An experimental and theoretical study of the ice accretion process during artificial and natural icing conditions

Real-time measurements of ice growth during artificial and natural icing conditions were conducted using an ultrasonic pulse-echo technique. This technique allows ice thickness to be measured with an accuracy of + or - 0.5 mm; in addition, the ultrasonic signal characteristics may be used to detect the presence of liquid on the ice surface and hence discern wet and dry ice growth behavior. Ice growth was measured on the stagnation line of a cylinder exposed to artificial icing conditions in the NASA Lewis Icing Research Tunnel (IRT), and similarly for a cylinder exposed in flight to natural icing conditions. Ice thickness was observed to increase approximately linearly with exposure time during the initial icing period. The ice accretion rate was found to vary with cloud temperature during wet ice growth, and liquid runback from the stagnation region was inferred. A steady-state energy balance model for the icing surface was used to compare heat transfer characteristics for IRT and natural icing conditions. Ultrasonic measurements of wet and dry ice growth observed in the IRT and in flight were compared with icing regimes predicted by a series of heat transfer coefficients. The heat transfer magnitude was generally inferred to be higher for the IRT than for the natural icing conditions encountered in flight. An apparent variation in the heat transfer magnitude was also observed for flights conducted through different natural icing-cloud formations.

Kirby, Mark S.↗

Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors

Cloud drop number concentrations (N d ) over the western North Atlantic Ocean (WNAO) are generally highest during the winter (DJF) and lowest in summer (JJA), in contrast to aerosol proxy variables (aerosoloptical depth, aerosol index, surface aerosol mass concentrations, surface cloud condensation nuclei (CCN) concentrations) that generally peak inspring (MAM) and JJA with minima in DJF. Using aircraft, satellite remote sensing, ground-based in situ measurement data, and reanalysis data, we characterize factors explaining the divergent seasonal cycles and furthermore probe into factors influencing N d on seasonal timescales. The results can be summarized well by features most pronounced in DJF, including features associated with cold-air outbreak (CAO) conditions such as enhanced values of CAO index, planetary boundary layer height (PBLH),low-level liquid cloud fraction, and cloud-top height, in addition to winds aligned with continental outflow. Data sorted into high- and low-N d days in each season, especially in DJF, revealed that all of these conditions were enhanced on the high-N d days, including reduced sea level pressure and stronger wind speeds. Although aerosols may be more abundant in MAM and JJA, the conditions needed to activate those particles into cloud droplets are weaker than in colder months, which is demonstrated by calculations of the strongest (weakest) aerosol indirect effects in DJF (JJA) based on comparing N d to perturbations in four different aerosol proxy variables (total and sulfate aerosol optical depth, aerosol index, surface mass concentration of sulfate). We used three machine learning models and up to 14 input variables to infer about most influential factors related to N d for DJF and JJA, with the best performance obtained with gradient-boosted regression tree (GBRT) analysis. The model results indicated that cloud fraction was the most important input variable, followed by some combination (depending on season) of CAO index and surface mass concentrations of sulfate and organic carbon. Future work is recommended to further understand aspects uncovered here such as impacts of free tropospheric aerosol entrainment on clouds, degree of boundary layer coupling, wet scavenging, and giant CCN effects on aerosol–N d relationships, updraft velocity, and vertical structure of cloud properties such as adiabaticity that impact the satellite estimation of N d .

54 ENVIRONMENTAL SCIENCES↗

Detailed Modeling and Analysis of the CPFM Dataset

A quantitative understanding of photolysis rate coefficients (or "j-values") is essential to determining the photochemical reaction rates that define ozone loss and other crucial processes in the atmosphere. j-Values can be calculated with radiative transfer models, derived from actinic flux observations, or inferred from trace gas measurements. The principal objective of this study is to cross-validate j-values from the Composition and Photodissociative Flux Measurement (CPFM) instrument during the Photochemistry of Ozone Loss in the Arctic Region In Summer (POLARIS) and SAGE I11 Ozone Loss and Validation Experiment (SOLVE) field campaigns with model calculations and other measurements and to use this detailed analysis to improve our ability to determine j-values. Another objective is to analyze the spectral flux from the CPFM (not just the j-values) and, using a multi-wavelength/multi-species spectral fitting technique, determine atmospheric composition.

Swartz, William H.↗

Intra-seasonal Oscillations Inferred from SABER (TIMED) and MLS (UARS) Temperature Measurements

In the zonal mean meridional winds of the upper mesosphere, intra-seasonal oscillations with periods between 1 and 4 months have been inferred from UARS measurements and independently predicted with the Numerical Spectral Model WSM). The wind oscillations tend to be confined to low latitudes and appear to be driven, at least in part, by small-scale gravity waves propagating in the meridional direction. Winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. Investigating this phenomenon, we have analyzed SABER temperatures from TIMED in the altitude range between 55 and 95 km to delineate with an empirical model, the year-long variability of the migrating tides and zonal mean components. The inferred seasonal variations of the diurnal tide, characterized by amplitude maxima near equinox, are in substantial agreement with UARS observations and results from the NSM. For the zonal mean, the dominant seasonal variations in the SABER temperatures, with annual (12 months) and semiannual (6 months) periodicities, agree well with those derived from UARS measurements. The intra-seasonal variations with periods between 2 and 4 months have amplitudes close to 2 K, almost half as large as those for the dominant seasonal variations. Their amplitudes are in qualitative agreement with the corresponding values inferred from UARS during different years. The SABER and UARS temperature variations reveal pronounced hemispherical asymmetries, consistent with meridional wind oscillations across the equator. The phase of the semi-annual temperature oscillations from the NSM agrees with the observations from UARS and SABER. But the amplitudes are systematically smaller, which may indicate that planetary waves are more important than is allowed for in the model. For the shorter-period intra-seasonal variations, which can be generated by gravity wave drag, the model results are generally in better agreement with the observations.

Huang, F. T.↗

Spike-and-Slab Shrinkage Priors for Structurally Sparse Bayesian Neural Networks

Network complexity and computational efficiency have become increasingly significant aspects of deep learning. Sparse deep learning addresses these challenges by recovering a sparse representation of the underlying target function by reducing heavily overparameterized deep neural networks. Specifically, deep neural architectures compressed via structured sparsity (e.g., node sparsity) provide low-latency inference, higher data throughput, and reduced energy consumption. In this article, we explore two well-established shrinkage techniques, Lasso and Horseshoe, for model compression in Bayesian neural networks (BNNs). To this end, we propose structurally sparse BNNs, which systematically prune excessive nodes with the following: 1) spike-and-slab group Lasso (SS-GL) and 2) SS group Horseshoe (SS-GHS) priors, and develop computationally tractable variational inference, including continuous relaxation of Bernoulli variables. We establish the contraction rates of the variational posterior of our proposed models as a function of the network topology, layerwise node cardinalities, and bounds on the network weights. Furthermore, we empirically demonstrate the competitive performance of our models compared with the baseline models in prediction accuracy, model compression, and inference latency.

97 MATHEMATICS AND COMPUTING↗

Data Analysis with Graphical Models: Software Tools

Probabilistic graphical models (directed and undirected Markov fields, and combined in chain graphs) are used widely in expert systems, image processing and other areas as a framework for representing and reasoning with probabilities. They come with corresponding algorithms for performing probabilistic inference. This paper discusses an extension to these models by Spiegelhalter and Gilks, plates, used to graphically model the notion of a sample. This offers a graphical specification language for representing data analysis problems. When combined with general methods for statistical inference, this also offers a unifying framework for prototyping and/or generating data analysis algorithms from graphical specifications. This paper outlines the framework and then presents some basic tools for the task: a graphical version of the Pitman-Koopman Theorem for the exponential family, problem decomposition, and the calculation of exact Bayes factors. Other tools already developed, such as automatic differentiation, Gibbs sampling, and use of the EM algorithm, make this a broad basis for the generation of data analysis software.

Buntine, Wray L.↗

The Livingstone Model of a Main Propulsion System

Livingstone is a discrete, propositional logic-based inference engine that has been used for diagnosis of physical systems. We present a component-based model of a Main Propulsion System (MPS) and say how it is used with Livingstone (L2) in order to implement a diagnostic system for integrated vehicle health management (IVHM) for the Propulsion IVHM Technology Experiment (PITEX). We start by discussing the process of conceptualizing such a model. We describe graphical tools that facilitated the generation of the model. The model is composed of components (which map onto physical components), connections between components and constraints. A component is specified by variables, with a set of discrete, qualitative values for each variable in its local nominal and failure modes. For each mode, the model specifies the component's behavior and transitions. We describe the MPS components' nominal and fault modes and associated Livingstone variables and data structures. Given this model, and observed external commands and observations from the system, Livingstone tracks the state of the MPS over discrete time-steps by choosing trajectories that are consistent with observations. We briefly discuss how the compiled model fits into the overall PITEX architecture. Finally we summarize our modeling experience, discuss advantages and disadvantages of our approach, and suggest enhancements to the modeling process.

Bajwa, Anupa↗

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

SBbadger: biochemical reaction networks with definable degree distributions

Abstract Motivation An essential step in developing computational tools for the inference, optimization and simulation of biochemical reaction networks is gauging tool performance against earlier efforts using an appropriate set of benchmarks. General strategies for the assembly of benchmark models include collection from the literature, creation via subnetwork extraction and de novo generation. However, with respect to biochemical reaction networks, these approaches and their associated tools are either poorly suited to generate models that reflect the wide range of properties found in natural biochemical networks or to do so in numbers that enable rigorous statistical analysis. Results In this work, we present SBbadger, a python-based software tool for the generation of synthetic biochemical reaction or metabolic networks with user-defined degree distributions, multiple available kinetic formalisms and a host of other definable properties. SBbadger thus enables the creation of benchmark model sets that reflect properties of biological systems and generate the kinetics and model structures typically targeted by computational analysis and inference software. Here, we detail the computational and algorithmic workflow of SBbadger, demonstrate its performance under various settings, provide sample outputs and compare it to currently available biochemical reaction network generation software. Availability and implementation SBbadger is implemented in Python and is freely available at https://github.com/sys-bio/SBbadger and via PyPI at https://pypi.org/project/SBbadger/. Documentation can be found at https://SBbadger.readthedocs.io. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Dynamics and Power Balance of Near Unity Target Gain Inertial Confinement Fusion Implosions

The change in the power balance, temporal dynamics, emission weighted size, temperature, mass, and areal density of inertially confined fusion plasmas have been quantified for experiments that reach target gains up to 0.72. It is observed that as the target gain rises, increased rates of self-heating initially overcome expansion power losses. This leads to reacting plasmas that reach peak fusion production at later times with increased size, temperature, mass and with lower emission weighted areal densities. Analytic models are consistent with the observations and inferences for how these quantities evolve as the rate of fusion self-heating, fusion yield, and target gain increase. At peak fusion production, it is found that as temperatures and target gains rise, the expansion power loss increases to a near constant ratio of the fusion self-heating power. Finally, this is consistent with models that indicate that the expansion losses dominate the dynamics in this regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An initial framework for the rapid qualification of long-term creep rupture strength via microstructural modeling

This report describes the development and testing of a new method for extrapolating short-term creep rupture test data to predict long-term rupture strength. The goal of this work is to reduce the time required to qualify new materials for nuclear service by reducing the lead time required for dedicated, long-term material testing to establish key long-term material properties. The new approach described here uses a physics-based model to predict the long-term creep rupture strength of 316H stainless steel using only short-term test data. The key idea is to use Bayesian inference to find the statistical distribution of the model parameters that best explain the short-term rupture data. Because the model is physics-based these parameters are all microstructural quantities that can be measured through detailed material characterization experiments. The Bayesian prior distributions provide a means for incorporating this characterization data into the final model to improve the accuracy of the long-term model predictions. However, where such data is not available the process still produces an accurate model based on an uniformed prior. Our hypothesis is that this approach more accurately extrapolates the short-term test data when compared to current, empirical methods. The report proves this hypothesis using actual long-term rupture data available for 316H, including tests with rupture times greater than 200,000 hours. The general approach developed here could be applied to other materials and other time-dependent material properties. Applying this new technique to develop long-term qualified material properties, potentially in conjunction with other accelerated qualification approaches like staggered qualification test programs, could greatly reduce the time required to qualify new materials for nuclear service.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ocean Heat and Carbon Uptake in Transient Climate Change: Identifying Model Uncertainty

Global warming on decadal and centennial timescales is mediated and ameliorated by the oceansequestering heat and carbon into its interior. Transient climate change is a function of the efficiency by whichanthropogenic heat and carbon are transported away from the surface into the ocean interior (Hansen et al. 1985).Gregory and Mitchell (1997) and Raper et al. (2002) were the first to identify the importance of the ocean heat uptakeefficiency in transient climate change. Observational estimates (Schwartz 2012) and inferences from coupledatmosphere-ocean general circulation models (AOGCMs; Gregory and Forster 2008; Marotzke et al. 2015), suggest thatocean heat uptake efficiency on decadal timescales lies in the range 0.5-1.5 W/sq m/K and is thus comparable to theclimate feedback parameter (Murphy et al. 2009). Moreover, the ocean not only plays a key role in setting the timing ofwarming but also its regional patterns (Marshall et al. 2014), which is crucial to our understanding of regional climate,carbon and heat uptake, and sea-level change. This short communication is based on a presentation given by A.Romanou at a recent workshop, Oceans Carbon and Heat Uptake: Uncertainties and Metrics, co-hosted by US CLIVARand OCB. As briefly reviewed below, we have incomplete but growing knowledge of how ocean models used in climatechange projections sequester heat and carbon into the interior. To understand and thence reduce errors and biases inthe ocean component of coupled models, as well as elucidate the key mechanisms at work, in the final section we outlinea proposed model intercomparison project named FAFMIP. In FAFMIP, coupled integrations would be carried out withprescribed overrides of wind stress and freshwater and heat fluxes acting at the sea surface.

periodic variations↗

Differentially Private Adaptive Noise Injection (DP-ANI) v1.0

Location data is collected from users continuously to understand their mobility patterns. Releasing the user trajectories may compromise user privacy. Therefore, the general practice is to release aggregated location datasets. However, private information may still be inferred from an aggregated version of location trajectories. Differential privacy (DP) protects the query output against inference attacks regardless of background knowledge. This software implements a differential privacy-based privacy model that protects the user's origins and destinations from being inferred from aggregated mobility datasets. This is achieved by injecting Planar Laplace noise to the user origin and destination GPS points. The noisy GPS points are then transformed into a link representation using a link-matching algorithm. Finally, the link trajectories form an aggregated mobility network. The injected noise level is selected using the Sparse Vector Mechanism. This DP selection mechanism considers the link density of the location and the functional category of the localized links. Compared to the different baseline models, including a k-anonymity method, our differential privacy-based aggregation model offers query responses that are close to the raw data in terms of aggregate statistics at both the network and trajectory-levels with maximum 9% deviation from the baseline in terms of network length.

Peisert, Sean [Lawrence Berkeley National Laborato↗

Applications of emulation and Bayesian methods in heavy-ion physics

Abstract Heavy-ion collisions provide a window into the properties of many-body systems of deconfined quarks and gluons. Understanding the collective properties of quarks and gluons is possible by comparing models of heavy-ion collisions to measurements of the distribution of particles produced at the end of the collisions. These model-to-data comparisons are extremely challenging, however, because of the complexity of the models, the large amount of experimental data, and their uncertainties. Bayesian inference provides a rigorous statistical framework to constrain the properties of nuclear matter by systematically comparing models and measurements. This review covers model emulation and Bayesian methods as applied to model-to-data comparisons in heavy-ion collisions. Replacing the model outputs (observables) with Gaussian process emulators is key to the Bayesian approach currently used in the field, and both current uses of emulators and related recent developments are reviewed. The general principles of Bayesian inference are then discussed along with other Bayesian methods, followed by a systematic comparison of seven recent Bayesian analyses that studied quark-gluon plasma properties, such as the shear and bulk viscosities. The latter comparison is used to illustrate sources of differences in analyses, and what it can teach us for future studies.

Paquet, Jean-François (ORCID:0000000187368171)↗

The influence of polar heterogeneous processes on reactive chlorine at middle latitudes - Three dimensional model implications

Three-dimensional model calculations with the NASA/GSFC chemistry and transport model have been designed to consider the impact of heterogeneous processes occurring on polar stratospheric clouds (PSCs) in the Arctic vortex on the HCl distribution. By examining the HCl concentration for a calculation with PSCs relative to a calculation with gas phase chemistry only, the impact of polar processing on reactive chlorine species at middle latitudes is inferred. Results from the chemistry and transport model reproduce basic features of the ClO measurements (Toohey et al., 1991), which were made on the ferry flights of the ER-2 from Stavanger, Norway to Moffett Field, California via Wallops Island, Virginia on February 20 and 21, 1989. The model indicates that perturbed air which is contained within the polar vortex during winter is not homogeneously mixed, and that the ferry flights were made through air with the largest conversion of HCl to reactive chlorine that is seen at middle latitudes.

Douglass, Anne R.↗

Pluto’s Sputnik Planitia: Composition of geological units from infrared spectroscopy

We have compared spectroscopic data of Sputnik Planitia on Pluto, as acquired by New Horizons’ Linear Etalon Imaging Spectral Array (LEISA) instrument, to the geomorphology as mapped by White et al. (2017) using visible and panchromatic imaging acquired by the LOng-Range Reconnaissance Imager (LORRI) and the Multi-spectral Visible Imaging Camera (MVIC). We have focused on 13 of the geologic units identified by White et al. (2017), which include the plains and mountain units contained within the Sputnik basin. We divided the map of Sputnik Planitia into 15 provinces, each containing one or more geologic units, and we use LEISA to calculate the average spectra of the units inside the 15 provinces. Hapke-based modeling was then applied to the average spectra of the units to infer their surface composition, and to determine if the composition resulting from the modeling of LEISA spectra reflects the geomorphologic analyses of LORRI data, and if areas classified as being the same geologically, but which are geographically separated, share a similar composition. We investigated the spatial distribution of the most abundant ices on Pluto’s surface - CH4, N2, CO, H2O, and a non-ice component presumed to be a macromolecular carbon-rich material, termed a tholin, that imparts a positive spectral slope in the visible spectral region and a negative spectral slope longward of ~1.1 μm. Because the exact nature of the non-ice component is still debated and because the negative spectral slope of the available tholins in the near infrared does not perfectly match the Pluto data, for spectral modeling purposes we reference it generically as the negative spectral slope endmember (NSS endmember). We created maps of variations in the integrated band depth (from LEISA data) and areal mass fraction (from the modeling) of the components. The analysis of correlations between the occurrences of the endmembers in the geologic units led to the observation of an anomalous suppression of the strong CH4 absorption bands in units with compositions that are dominated by H2O ice and the NSS endmember. Exploring the mutual variation of the CH4 and N2 integrated band depths with the abundance of crystalline H2O and NSS endmember revealed that the NSS endmember is primarily responsible for the suppression of CH4 absorptions in mountainous units located along the western edge of Sputnik Planitia. Our spectroscopic analyses have provided additional insight into the geological processes that have shaped Sputnik Planitia. A general increase in volatile abundance from the north to the south of Sputnik Planitia is observed. Such an increase first observed and interpreted by Protopapa et al., 2017 and later confirmed by climate modeling (Bertrand et al., 2018) is expressed geomorphologically in the form of preferential deposition of N2 ice in the upland and mountainous regions bordering the plains of southern Sputnik Planitia. Relatively high amounts of pure CH4 are seen at the southern Tenzing Montes, which are a natural site for CH4 deposition owing to their great elevation and the lower insolation they are presently receiving. The NSS endmember correlates the existence of tholins within certain units, mostly those coating the low-latitude mountain ranges that are co-latitudinal with the tholin-covered Cthulhu Macula. The spectral analysis has also revealed compositional differences between the handful of occurrences of northern non-cellular plains and the surrounding cellular plains, all of which are located within the portion of Sputnik Planitia that is presently experiencing net sublimation of volatiles, and which do not therefore exhibit a surface layer of bright, freshly-deposited N2 ice. The compositional differences between the cellular and non-cellular plains here hint at the effectiveness of convection in entraining and trapping tholins within the body of the cellular plains, while preventing the spread of such tholins to abutting non-cellular plains.

Pluto's Sputnik Planitia↗

Modeling the Swift BAT trigger algorithm with machine learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift/BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of ≳97% (≲3% error), which is a significant improvement on a cut in GRB flux, which has an accuracy of 89.6% (10.4% error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of n{sub 0}∼0.48{sub −0.23}{sup +0.41} Gpc{sup −3} yr{sup −1} with power-law indices of n{sub 1}∼1.7{sub −0.5}{sup +0.6} and n{sub 2}∼−5.9{sub −0.1}{sup +5.7} for GRBs above and below a break point of z{sub 1}∼6.8{sub −3.2}{sup +2.8}. This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting.

79 ASTRONOMY AND ASTROPHYSICS↗