Search NASA⌕ Search

SEARCH · Search NASA

Results for “inference methods”

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,387 records · Page 77

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hail Storm Risk Assessment Using Space-Borne Remote Sensing Observations and Reanalyses

Much of the world is impacted by severe thunderstorms, but whether they become disasters depends upon resilience--our capacity to prepare, mitigate, respond, and recover. Hail is the costliest severe weather hazard for the insurance industry, generating ~70% of severe convective storm losses due to damage to assets such as homes, businesses, agriculture, and infrastructure. Most insurance companies do not reserve enough capital to cover catastrophes, so they acquire reinsurance. The reinsurance industry uses catastrophe models (CatModels) to statistically estimate risk to an insurer’s portfolio. Hail CatModels are developed with climatologies that define hailstorm frequency and severity. Hail-prone areas can be defined using hail reports from trained spotters, the media, and the general public. Extremely severe hail (2+ inch diameter) occurs nearly every day across the world. Weather radars can detect hail because hailstones strongly reflect microwave signals that they emit. However, hail climatologies are difficult to derive because hail covers small areas and there are neither hail reporting mechanisms (e.g. website or mobile app) nor radar networks in most places outside the US and Europe. This lack of ground truth on severe hail puts society and economies at risk. Hail is generated within storms by strong updrafts. These updrafts exhibit unique signatures in NASA and other agency satellite observations, offering new opportunities for hailstorm analysis. Geostationary (GEO) visible and infrared imagery has been collected for ~15-25 years across the world (region dependent) and methods have been developed at NASA Langley Research Center (LaRC) to detect hailstorm updrafts using GEO imagery. Climatological GEO updraft data has been used by Willis Towers Watson (WTW), a leader in catastrophe risk assessment for the insurance industry, and Karlsruhe Institute of Technology to develop CatModels over Europe and Australia. Hail can also be inferred with passive microwave imagery collected by low-Earth-orbiting sensors such as the GPM GMI, TRMM TMI, AMSR-E, AMSR-2, SSM/I, and SSMIS over the last 20+ years using methods developed at the Marshall Space Flight Center (MSFC). Hailstorms generate enhanced lightning flash rates that can be tracked using new GOES-R series GEO Lightning Mapping (GLM) imagery. Atmospheric reanalyses can be used to define favorable hailstorm environments for combination with the satellite-based storm detections. This presentation will describe a framework for developing continental to global hail climatologies and CatModels based on NASA satellite data and capabilities. This is a collaboration between LaRC and MSFC, WTW, and partners in Brazil, Argentina, and South Africa. This project seeks to mitigate hail disasters by aiding development of new satellite-based severe storm nowcasting tools by regional partners and developing climatologies to improve societal understanding of hail frequency. GEOO visible and infrared metrics of storm intensity, environmental conditions based on reanalyses, spotter hail reports and radar MESH observations are intercompared to quantify the detectability of hailstorms, and our ability to discriminate hailstorms from other severe storms. We are also maturing methods using land surface imaging satellite data (e.g. MODIS, Landsat, Sentinel 1 and 2) to identify hail damage to agriculture. Work with WTW will improve socioeconomic resilience through development of new CatModels. Southern Brazil, Uruguay, Paraguay, and Argentina feature some of the most intense thunderstorms on Earth. South America and South Africa are developing insurance markets of interest to WTW clients, and is similar to other regions routinely impacted by hail that do not have comprehensive hail reporting or radars to assess hailstorm frequency. Project datasets will be made available via online GIS-enabled tools developed at the LaRC Atmospheric Science Data Center (ASDC) which will visualize data and provide it in multiple formats for use in a wide range of open source and commercial tools.

Kristopher Michael Bedka↗

Fast Ray Tracing of Lunar Digital Elevation Models

Ray-tracing (RT) of Lunar Digital Elevation Models (DEM)'s is performed to virtually derive the degree of radiation incident to terrain as a function of time, orbital and ephemeris constraints [I- 4]. This process is an integral modeling process in lunar polar research and exploration due to the present paucity of terrain information at the poles and mission planning activities for the anticipated spring 2009 launch of the Lunar Reconnaissance Orbiter (LRO). As part of the Lunar Exploration Neutron Detector (LEND) and Lunar Crater Observation and Sensing Satellite (LCROSS) preparations RI methods are used to estimate the critical conditions presented by the combined effects of high latitude, terrain and the moons low obliquity [5-7]. These factors yield low incident solar illumination and subsequently extreme thermal, and radiation conditions. The presented research uses RT methods both for radiation transport modeling in space and regolith related research as well as to derive permanently shadowed regions (PSR)'s in high latitude topographic minima, e.g craters. These regions are of scientific and human exploration interest due to the near constant low temperatures in PSRs, inferred to be < 100 K. Hydrogen is thought to have accumulated in PSR's through the combined effects of periodic cometary bombardment and/or solar wind processes, and the extreme cold which minimizes hydrogen sublimation [8-9]. RT methods are also of use in surface position optimization for future illumination dependent on surface resources e.g. power and communications equipment.

McClanahan, Timothy P.↗

Shepherding Metadata Through the Building Lifecycle

Many different digital representations of a building are produced over the course of its lifecycle. These representations contain the metadata required to support different stages of the building, from initial planning and design, to construction and commissioning, through operations, audits, retrofits and maintenance. However, because of differences in the semantics, structure and syntax of these representations, the metadata they contain is not interoperable. We present a novel method for leveraging these representations to create a unified, authoritative Brick metadata model for a building that can be continually maintained over the course of the building lifecycle. A simple synchronization protocol relays inferred Brick metadata from existing metadata sources such as gbXML, BuildingSync, Project Haystack and Modelica to a central integration server, which merges the metadata into a valid Brick model.

Fierro, Gabe↗

Remote Sensing of Aerosol Water Fraction, Dry Size Distribution and Soluble Fraction Using Multi-Angle, Multi-Spectral Polarimetry

A framework to infer volume water fraction, soluble fraction and dry size distributions of fine-mode aerosol from multi-angle, multi-spectral polarimetry retrievals of column-averaged ambient aerosol properties is presented. The method is applied to observations of the Research Scanning Polarimeter (RSP) obtained during two NASA aircraft campaigns, namely the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) and the Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP2Ex). All aerosol retrievals are statistically evaluated using in situ data. Volume water fraction is inferred from the retrieved ambient real part of the refractive index, assuming a dry refractive index of 1.54 and by applying a volume mixing rule to obtain the effective ambient refractive index. The uncertainties in inferred volume water fraction resulting from this simplified model are discussed and estimated to be lower than 0.2 and decreasing with increasing volume water fraction. The daily mean retrieved volume water fractions correlate well with the in situ values with a mean absolute difference of 0.09. Polarimeter-retrieved ambient effective radius for daily data is shown to increase as a function of volume water fraction as expected. Furthermore, the effective variance of the size distributions also increases with increasing effective radius, which we show is consistent with an external mixture of soluble and insoluble aerosol. The relative variations of effective radius and variance over an observation period are then used to estimate the soluble fraction of the aerosol. Daily results of soluble fraction correlate well with in situ-observed sulfate mass fraction with a correlation coefficient of 0.79. Subsequently, inferred water and soluble fractions are used to derive dry fine-mode size distributions from their ambient counterparts. While dry effective radii obtained in situ and from RSP show similar ranges, in situ values are generally substantially smaller during the ACTIVATE deployments, which may be due to biases in RSP retrievals or in the in situ observations, or both. Both RSP and in situ observations indicate the dominance of aerosol with low hygroscopicity during the ACTIVATE and CAMP2Ex campaigns. Furthermore, RSP indicates a high degree of external mixing of particles with low and high hygroscopicity. These retrievals of fine-mode water volume fraction and soluble fraction may be used for the evaluation of water uptake in atmospheric models. Furthermore, the framework allows us to estimate the variation in the concentration of fine-mode aerosol larger than a specific dry radius limit, which can be used as a proxy for the variation in cloud condensation nucleus concentrations. This framework may be applied to multi-angle, multi-spectral satellite data expected to be available in the near future.

Amount of water in particulate matter↗

Predicting Fracture Porosity Evolution in Sandstone (Final Report)

To better understand porosity, strength, chemical reactivity, and patterns of fractures, we developed methods using mineral deposits to unravel how fracture growth and diagenesis interact to create and destroy fracture porosity. Quartz and dolomite cement textures and associated fluid inclusion assemblage sequences and thermal histories provided data to infer the duration and rates at which fractures open. We developed and tested a model that accounts for how fracture porosity and rock properties change as a function of thermal exposure and how fracture size and arrangement evolve. Effects of cement deposits on growth were investigated with new geomechanical models and new methods to quantify fracture spatial arrangement. We tested these concepts in a wide range of rock types using outcrop and core-based datasets. We documented opening histories of fractures and fracture arrays in contrasting tectonic and thermal regimes. Our methods set a new standard for fracture characterization and our diagenetic-geomechanical model accurately predicts fluid flow characteristics of fractured rocks based on coupled effects of diagenesis and deformation. Results represent a significant advance in fundamental understanding of fractures in the Earth and have application to improving subsurface engineering operations in a wide range of contexts including geothermal and sequestration applications. The breakthrough of successfully reconstructing the timing and opening rates of parts of single fractures and portions of fracture arrays to document the growth and linkage of fracture arrays was an accomplishment of this grant that is being followed up in our current BES grant DE-SC0022968, ‘Reconstructing and Predicting Fracture Pattern Evolution‘

58 GEOSCIENCES↗

Intrepid MCMC: Metropolis-Hastings with exploration

In engineering examples, one often encounters the need to sample from unnormalized distributions with complex shapes that may also be implicitly defined through a physical or numerical simulation model, making it computationally expensive to evaluate the associated density function. For such cases, MCMC has proven to be an invaluable tool. Random-walk Metropolis Methods (also known as Metropolis-Hastings (MH)), in particular, are highly popular for their simplicity, flexibility, and ease of implementation. However, most MH algorithms suffer from significant limitations when attempting to sample from distributions with multiple modes (particularly disconnected ones). Here, in this paper, we present Intrepid MCMC - a novel MH scheme that utilizes a simple coordinate transformation to significantly improve the mode-finding ability and convergence rate to the target distribution of random-walk Markov chains while retaining most of the simplicity of the vanilla MH paradigm. Through multiple examples, we showcase the improvement in the performance of Intrepid MCMC over vanilla MH for a wide variety of target distribution shapes. We also provide an analysis of the mixing behavior of the Intrepid Markov chain, as well as the efficiency of our algorithm for increasing dimensions. A thorough discussion is presented on the practical implementation of the Intrepid MCMC algorithm. Finally, its utility is highlighted through a Bayesian parameter inference problem for a two-degree-of-freedom oscillator under free vibration.

97 - MATHEMATICS AND COMPUTING↗

A Probabilistic Model-Based Diagnostic Framework for Nuclear Engineering Systems

A fault diagnostic framework was investigated in this study for applications in thermal–hydraulic systems of nuclear power plants. The proposed framework consists of quantitative model-based diagnosis, statistical change detection and probabilistic reasoning. The use of physics-based diagnostic models provides high detection sensitivity and allows noise and measurement uncertainty to be incorporated robustly. Performance-related parametric models for each component are constructed based on first principles. Numerical model residuals are generated using the concept of analytical redundancy. Statistical change detection methods are employed to detect non-zero residuals in the presence of uncertainty. The diagnosis task is performed using Bayesian inference to detect and localize possible faults. Application to a single-phase heat exchanger for demonstration showed that the proposed probabilistic framework can provide improved results in comparison with traditional approaches while remaining less sensitive to false alarms in the presence of measurement and modeling uncertainty.

Bayesian network↗

Automatic and rapid calibration of urban building energy models by learning from energy performance database

Urban building energy modeling (UBEM) is attracting increasing attention in the energy modeling filed. Unlike modeling a single building using detailed building systems information, UBEM generally uses limited high-level building stock data to infer default assumptions about building characteristics and operations. Additionally, this practice inherently brings uncertainty to UBEM. This study introduced a novel method of automatic and rapid calibration of UBEM based on the annual electricity and natural gas energy use data by learning the correlations between crucial model input parameters and the building energy use from the reference building models. A case study was presented to calibrate 72 large office buildings built before 1978 in San Francisco. Seventeen model parameters were selected and Monte Carlo sampling was used to create 1000 samples that reasonably represent the parameter space. Then 1000 simulations were performed for the reference building model to create an energy performance database. The results showed that by learning from the energy performance database, it took less than four simulation runs on average to calibrate a building model. After the calibration, the distributions of each parameter were obtained to replace their single predefined default values. For example, the default lighting power density of 21.39 W/m 2 was calibrated to be 7.50 W/m 2 on average. The case study successfully demonstrated the effectiveness of the novel calibration method for UBEM in the mild climate. The method will be further tested in future for other climate zones and other building types.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Three-dimensional realizations of flood flow in large-scale rivers using the neural fuzzy-based machine-learning algorithms

Machine learning methods have been extensively used to study the dynamics of complex fluid flows. One such algorithm, known as adaptive neural fuzzy inference system (ANFIS), can generate data-driven predictions for flow fields, but has not been applied to natural geophysical flows in large-scale rivers. Herein, we demonstrate the potential of ANFIS to produce three-dimensional (3D) realizations of the instantaneous flood flow field in several large-scale, virtual meandering rivers. The 3D dynamics of flood flow in large-scale rivers were obtained using large-eddy simulation (LES). The LES results, i.e., the 3D velocity components, were employed to train the learnable coefficients of an ANFIS. Further, the trained ANFIS, along with a few time-steps of LES results (precursor data) were then used to produce 3D realizations of flood flow fields in large-scale rivers with geometries other than the one the ANFIS was trained with. We also used the trained ANFIS to generate 3D realizations of river flow at a discharge other than that the ANFIS was trained with. The flow field results obtained from ANFIS were validated using separate LES runs to assess the accuracy of the 3D instantaneous realizations of the machine learning algorithm. An error analysis was conducted to quantify the discrepancies among the ANFIS and LES results for various flood flow predictions in large-scale rivers.

54 ENVIRONMENTAL SCIENCES↗

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Continuous-time probabilistic models for longitudinal electronic health records

Analysis of longitudinal Electronic Health Record (EHR) data is an important goal for precision medicine. Difficulty in applying Machine Learning (ML) methods, either predictive or unsupervised, stems in part from the heterogeneity and irregular sampling of EHR data. Here, we present an unsupervised probabilistic model that captures nonlinear relationships between variables over continuous-time. This method works with arbitrary sampling patterns and captures the joint probability distribution between variable measurements and the time intervals between them. Inference algorithms are derived that can be used to evaluate the likelihood of future using under a trained model. As an example, we consider data from the United States Veterans Health Administration (VHA) in the areas of diabetes and depression. Likelihood ratio maps are produced showing the likelihood of risk for moderate-severe vs minimal depression as measured by the Patient Health Questionnaire-9 (PHQ-9).

59 BASIC BIOLOGICAL SCIENCES↗

Assessing Metal Ion Assignment Accuracy in Protein Data Bank Models via Elemental Spectroscopy

Accurate representation of metal ions in macromolecular structures is critical for chemical interpretation, computational modeling, and machine-learning methods that rely on Protein Data Bank (PDB) entries. However, the elemental identity of metals modeled in crystallographic structures is often inferred indirectly and rarely validated experimentally. Here, we combine Particle Induced X-ray Emission (PIXE) and X-ray Fluorescence Spectroscopy (XRFS) to determine the elemental composition of protein samples used to generate 70 deposited metalloprotein crystal structures. By analyzing the original protein material employed for crystallization, but before the addition of crystallization buffer solutions, we assess whether the modeled metal ions in deposited structures are consistent with experimentally detectable elemental content. We find that in a majority of cases, the metals modeled in the corresponding PDB entries are inconsistent with the metals present in the protein samples before crystallization, or that additional metals are present but not represented in the structural models. Spectroscopic results were integrated with automated crystallographic validation metrics, including real-space Z-difference (RSZD) analysis and systematic rerefinement, to evaluate atomic-number mismatch at metal sites. PIXE and XRFS show strong agreement for dominant elemental signals and provide complementary, scalable approaches for identifying suspect metal assignments. This work does not address physiological or functional metalation but instead highlights a widespread data integrity issue in deposited macromolecular structures, PDB-wide. These results establish an experimentally corroborated link between elemental identity and crystallographic validation metrics, enabling the large-scale detection of chemically inconsistent annotations in structural databases used for computational modeling and machine learning.

Crystallization↗

Cosmological histories in bimetric gravity: a graphical approach

The bimetric generalization of general relativity has been proven to be able to give an accelerated background expansion consistent with observations. Apart from the energy densities coupling to one or both of the metrics, the expansion will depend on the cosmological constant contribution to each of them, as well as the three parameters describing the interaction between the two metrics. Even for fixed values of these parameters can several possible solutions, so called branches, exist. Different branches can give similar background expansion histories for the observable metric, but may have different properties regarding, for example, the existence of ghosts and the rate of structure growth. In this paper, we outline a method to find viable solution branches for arbitrary parameter values. We show how possible expansion histories in bimetric gravity can be inferred qualitatively, by picturing the ratio of the scale factors of the two metrics as the spatial coordinate of a particle rolling along a frictionless track. A particularly interesting example discussed is a specific set of parameter values, where a cosmological dark matter background is mimicked without introducing ghost modes into the theory.

79 ASTRONOMY AND ASTROPHYSICS↗

LARGER PLANET RADII INFERRED FROM STELLAR ''FLICKER'' BRIGHTNESS VARIATIONS OF BRIGHT PLANET-HOST STARS

Most extrasolar planets have been detected by their influence on their parent star, typically either gravitationally (the Doppler method) or by the small dip in brightness as the planet blocks a portion of the star (the transit method). Therefore, the accuracy with which we know the masses and radii of extrasolar planets depends directly on how well we know those of the stars, the latter usually determined from the measured stellar surface gravity, log g. Recent work has demonstrated that the short-timescale brightness variations ({sup f}licker{sup )} of stars can be used to measure log g to a high accuracy of ∼0.1-0.2 dex. Here, we use flicker measurements of 289 bright (Kepmag < 13) candidate planet-hosting stars with T {sub eff} = 4500-6650 K to re-assess the stellar parameters and determine the resulting impact on derived planet properties. This re-assessment reveals that for the brightest planet-host stars, Malmquist bias contaminates the stellar sample with evolved stars: nearly 50% of the bright planet-host stars are subgiants. As a result, the stellar radii, and hence the radii of the planets orbiting these stars, are on average 20%-30% larger than previous measurements had suggested.

79 ASTRONOMY AND ASTROPHYSICS↗

β decay of Ba 141

The β-decay strength function of nuclides produced in fission is important as it dictates the distribution of decay energy between electrons, neutrinos, and γ rays and so is critical for calculating decay heat in reactors and for estimating the reactor antineutrino spectrum. Several experimental techniques are available to determine this strength function, including electron spectroscopy, γ-ray calorimetry (TAGS spectroscopy), and detailed, high-resolution spectroscopy with modern large high-purity germanium arrays. This work investigates the decay of the well-known and strongly produced fission fragment 141 Ba. A beam of 141 Cs was implanted at the target position of the Gammasphere and the subsequent decay of the daughter 141 Ba was studied. Extensive decay spectroscopy was possible up to the decay Q value of 3.197(7) MeV, including a significant extension of the level scheme and detailed angular correlation measurements for all levels with greater than 0.25% β feeding. The distribution of the β-decay strength was then inferred and compared to previous calorimetric studies. In conclusion, the agreement was excellent and provides a benchmark for comparing strength function methods and data for a more detailed understanding of the structure of 141 La.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Probing Dark Matter Using Precision Measurements of Stellar Accelerations

Dark matter comprises the bulk of the matter in the Universe, but its particle nature and cosmological origin remain mysterious. Knowledge of the dark matter density distribution in the Milky Way Galaxy is crucial both to our understanding of the standard cosmological model and for grounding direct and indirect searches for the particles comprising dark matter. Current measurements of Galactic dark matter content rely on model assumptions to infer the forces acting upon stars from the distribution of observed velocities. Here, we propose to apply the precision radial velocity method, optimized in recent years for exoplanet astronomy, to measure the change in the velocity of stars over time, thereby providing a direct probe of the local gravitational potential in the Galaxy. Using numerical simulations, we develop a realistic strategy to observe the differential accelerations of stars in our Galactic neighborhood with next-generation telescopes, at the level of 10 –8 cm/s 2 . Our simulations show that detecting accelerations at this level with an ensemble of 10 3 stars requires the effect of stellar noise on radial velocity measurements to be reduced to < 10 cm/s. Furthermore, the measured stellar accelerations may then be used to extract the local dark matter density and morphological parameters of the density profile.

79 ASTRONOMY AND ASTROPHYSICS↗