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 91 records · Page 5

Extracting the Breakout Distance from the ECOT Trajectories: Gaussian Process Regression Approach

Enhanced Corner Turning (ECOT) experiments provide an important metric of performance of high explosive (HE) formulations. The breakout distance is a single scalar value that characterizes the corner turning efficiency of an HE. Extracting the breakout distance from the raw ECOT results, whether experimental or simulated, is a conceptually straightforward procedure which, however, is non-unique, especially in the presence of noise. More specifically, this procedure involves numerical smoothing and selecting particular values for parameters of this smoothing introduces human bias. In this work, we propose to use the Gaussian process regression to analyze ECOT results. This analysis involves the effective smoothing of the data, thus allowing for accurate extraction of the breakout distance. Most importantly, the parameters of this smoothing can be inferred from the ECOT data itself, rendering the approach effectively parameter-free and thus diminishing the human bias. An additional benefit of the Gaussian process regression, being a statistical inference method, is that not just the value of the breakout distance, but also its confidence interval can be extracted from the data. This report introduces the Gaussian process regression, as applied to ECOT, and demonstrates its usefulness by extracting the breakout distances for a selection of experimental and simulated data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

ARM-IRL: Adaptive Resilience Metric Quantification Using Inverse Reinforcement Learning

The resilience of safety-critical systems is gaining importance due to the rise in cyber and physical threats, especially within critical infrastructure. Traditional static resilience metrics may not capture dynamic system states, leading to inaccurate assessments and ineffective responses to cyber threats. This work aims to develop a data-driven, adaptive method for resilience metric learning. We propose a data-driven approach using inverse reinforcement learning (IRL) to learn a single, adaptive resilience metric. The method infers a reward function from expert control actions. Unlike previous approaches using static weights or fuzzy logic, this work applies adversarial inverse reinforcement learning (AIRL), training a generator and discriminator in parallel to learn the reward structure and derive an optimal policy. The proposed approach is evaluated on multiple scenarios: optimal communication network rerouting, power distribution network reconfiguration, and cyber–physical restoration of critical loads using the IEEE 123-bus system. The adaptive, learned resilience metric enables faster critical load restoration in comparison to conventional RL approaches.

97 MATHEMATICS AND COMPUTING↗

A Compound Poisson Generator Approach to Point-source Inference in Astrophysics

Abstract The identification and description of point sources is one of the oldest problems in astronomy, yet even today the correct statistical treatment for point sources remains one of the field’s hardest problems. For dim or crowded sources, likelihood-based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using compound Poisson generator (CPG) functionals that incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over non-Poissonian template fitting (NPTF)—an existing method—in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorized as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterization of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parameterization is presented for these priors that properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Inferring Type II-P Supernova Progenitor Masses from Plateau Luminosities

Abstract Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. presented bolometric light curves for a landscape of progenitors from spherically symmetric neutrino-driven core-collapse supernova (CCSN) simulations using an effective model. They find a tight relationship between the plateau luminosity of the Type II-P CCSN light curve and the terminal iron-core mass of the progenitor. Remarkably, this allows us to constrain progenitor properties with photometry alone. We analyze a large observational sample of Type II-P CCSN light curves and estimate a distribution of iron-core masses using the relationship of Barker et al. The inferred distribution matches extremely well with the distribution of iron-core masses from stellar evolutionary models and namely, contains high-mass iron cores that suggest contributions from very massive progenitors in the observational data. We use this distribution of iron-core masses to infer minimum and maximum masses of progenitors in the observational data. Using Bayesian inference methods to locate optimal initial mass function parameters, we find M min = 9.8 − 0.27 + 0.37 and M max = 24.0 − 1.9 + 3.9 solar masses for the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantile Data Analysis of Image Data

Quantile data analysis and functional statistical inference methods are introduced and applied to provide representations of spectral data which may lead to simple statistical discriminators effective for the estimation of ground truth from satellite spectral measurements. To estimate the ground truth of a pixel, the probability of each possible ground truth is estimated, given observed (estimated) quantile theoretic statistical characteristics of the multispectral image data corresponding to the pixel and its neighboring pixels. A strategy for determining which statistical characteristics discriminate best is described. Results are reported of quantile data analysis of an extensive collection of training files of image data.

Parzen, E.↗

Bayesian learning

In 1983 and 1984, the Infrared Astronomical Satellite (IRAS) detected 5,425 stellar objects and measured their infrared spectra. In 1987 a program called AUTOCLASS used Bayesian inference methods to discover the classes present in these data and determine the most probable class of each object, revealing unknown phenomena in astronomy. AUTOCLASS has rekindled the old debate on the suitability of Bayesian methods, which are computationally intensive, interpret probabilities as plausibility measures rather than frequencies, and appear to depend on a subjective assessment of the probability of a hypothesis before the data were collected. Modern statistical methods have, however, recently been shown to also depend on subjective elements. These debates bring into question the whole tradition of scientific objectivity and offer scientists a new way to take responsibility for their findings and conclusions.

Denning, Peter J.↗

Generic task problem solvers in Soar

Two trends can be discerned in research in problem solving architectures in the last few years. On one hand, interest in task-specific architectures has grown, wherein types of problems of general utility are identified, and special architectures that support the development of problem solving systems for those types of problems are proposed. These architectures help in the acquisition and specification of knowledge by providing inference methods that are appropriate for the type of problem. However, knowledge based systems which use only one type of problem solving method are very brittle, and adding more types of methods requires a principled approach to integrating them in a flexible way. Contrasting with this trend is the proposal for a flexible, general architecture contained in the work on Soar. Soar has features which make it attractive for flexible use of all potentially relevant knowledge or methods. But as the theory Soar does not make commitments to specific types of problem solvers or provide guidance for their construction. It was investigated how task-specific architectures can be constructed in Soar to retain as many of the advantages as possible of both approaches. Examples were used from the Generic Task approach for building knowledge based systems. Though this approach was developed and applied for a number of problems, the ideas are applicable to other task-specific approaches as well.

Johnson, Todd R.↗

Design and implementation of fuzzy logic controllers

The main objectives of our research are to present a self-contained overview of fuzzy sets and fuzzy logic, develop a methodology for control system design using fuzzy logic controllers, and to design and implement a fuzzy logic controller for a real system. We first present the fundamental concepts of fuzzy sets and fuzzy logic. Fuzzy sets and basic fuzzy operations are defined. In addition, for control systems, it is important to understand the concepts of linguistic values, term sets, fuzzy rule base, inference methods, and defuzzification methods. Second, we introduce a four-step fuzzy logic control system design procedure. The design procedure is illustrated via four examples, showing the capabilities and robustness of fuzzy logic control systems. This is followed by a tuning procedure that we developed from our design experience. Third, we present two Lyapunov based techniques for stability analysis. Finally, we present our design and implementation of a fuzzy logic controller for a linear actuator to be used to control the direction of the Free Flight Rotorcraft Research Vehicle at LaRC.

Abihana, Osama A.↗

Statistical studies of impulsive events at high latitudes

A statistical study has been made of the high-latitude impulsive events that were observed during the 1985-1986 South Pole Balloon Campaign. The events were selected by searching for unipolar pulses greater than or equal to 10 nT above background in the vertical component of the magnetic field on the ground and/or pedestal or 'W' shaped horizontal electric field perturbations greater than or equal to 10 mV/m in amplitude and accompanied by perturbations in the vertical electric field at balloon altitude. A main event list comprising 112 events was compiled from the 468 hours of data available. Three aspects of the events were examined: the solar wind conditions prior to the event, local time of observation, and intrinsic properties of the events. The local time distribution was obtained from the 112 entry main event list and was found to be nearly uniform across the dayside, with no midday gap. The event rate found using our low-amplitude selection criteria was 0.7 event/hr, comparable to expectations based on in situ studies of the magnetopause. A total of 42 events were found for which data were available from Interplanetary Monitoring Platform (IMP) 8. Of these events, 12 occurred when the Z(sub GSM) component (B(sub Z)) of the interplanetary magnetic field (IMF) was northward and 30 occurred when B(sub Z) was southward or fluctuating. Only three of the B(sub Z) northward cases and only five of the B(sub Z) southward cases were preceded by pressure pulses greater than 0.4 nPa in amplitude. Ten of the events were studied in detail by means of a model-fitting method discussed elsewhere. This method infers values of several parameters, including the total current flowing in a coaxial or monopole system and a two-dimensional dipole system. The intrinsic properties of the events showed that only approximately 10% of the total current contributed to momentum transfer to the high-latitude ionosphere, that the direction of the motion depended more on local time of observation than IMF B(sub y), and that events were usually several hundred kilometers in size. The observed B(sub z) control found in the 42 event list and the prevalence of coaxial current dominated events are inconsistent with the predictions of the pressure pulse model.

Lin, Z. M.↗

Web-based Weather Expert System (WES) for Space Shuttle Launch

The Web-based Weather Expert System (WES) is a critical module of the Virtual Test Bed development to support 'go/no go' decisions for Space Shuttle operations in the Intelligent Launch and Range Operations program of NASA. The weather rules characterize certain aspects of the environment related to the launching or landing site, the time of the day or night, the pad or runway conditions, the mission durations, the runway equipment and landing type. Expert system rules are derived from weather contingency rules, which were developed over years by NASA. Backward chaining, a goal-directed inference method is adopted, because a particular consequence or goal clause is evaluated first, and then chained backward through the rules. Once a rule is satisfied or true, then that particular rule is fired and the decision is expressed. The expert system is continuously verifying the rules against the past one-hour weather conditions and the decisions are made. The normal procedure of operations requires a formal pre-launch weather briefing held on Launch minus 1 day, which is a specific weather briefing for all areas of Space Shuttle launch operations. In this paper, the Web-based Weather Expert System of the Intelligent Launch and range Operations program is presented.

Bardina, Jorge E.↗

Variation in Global Chemical Composition of PM2.5: Emerging Results from SPARTAN

The Surface PARTiculate mAtter Network (SPARTAN) is a long-term project that includes characterization of chemical and physical attributes of aerosols from filter samples collected worldwide. This paper discusses the ongoing efforts of SPARTAN to define and quantify major ions and trace metals found in fine particulate matter (PM (sub 2.5). Our methods infer the spatial and temporal variability of PM (sub 2.5) in a cost-effective manner. Gravimetrically weighed filters represent multi-day averages of PM (sub 2.5), with a collocated nephelometer sampling air continuously. SPARTAN instruments are paired with AErosol RObotic NETwork (AERONET) sun photometers to better understand the relationship between ground-level PM (sub 2.5) and columnar aerosol optical depth (AOD). We have examined the chemical composition of PM (sub 2.5) at 12 globally dispersed, densely populated urban locations and a site at Mammoth Cave (US) National Park used as a background comparison. So far, each SPARTAN location has been active between the years 2013 and 2016 over periods of 2-26 months, with an average period of 12 months per site. These sites have collectively gathered over 10 years of quality aerosol data. The major PM (sub 2.5) constituents across all sites (relative contribution plus or minus Standard Deviation) are ammoniated sulfate (20 percent plus or minus 11 percent), crustal material (13.4 percent plus or minus 9.9 percent), equivalent black carbon (11.9 percent plus or minus 8.4 percent), ammonium nitrate (4.7 percent plus or minus 3.0 percent), sea salt (2.3 percent plus or minus 1.6 percent), trace element oxides (1.0 percent plus or minus 1.1 percent), water (7.2 percent plus or minus 3.3 percent) at 35 percent relative humidity, and residual matter (40 percent plus or minus 24 percent). Analysis of filter samples reveals that several PM (sub 2.5) chemical components varied by more than an order of magnitude between sites. Ammoniated sulfate ranges from 1.1 microns per cubic meter (Buenos Aires, Argentina) to 17 microns per cubic meter (Kanpur, India in the dry season). Ammonium nitrate ranged from 0.2 microns per cubic meter (Mammoth Cave, in summer) to 6.8 microns per cubic meter (Kanpur, dry season). Equivalent black carbon ranged from 0.7 microns per cubic meter (Mammoth Cave) to over 8 microns per cubic meter (Dhaka, Bangladesh and Kanpur, India). Comparison of SPARTAN vs. coincident measurements from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network at Mammoth Cave yielded a high degree of consistency for daily PM (sub 2.5) (r squared equals 0.76, slope equals 1.12), daily sulfate (r squared equals 0.86, slope equals 1.03), and mean fractions of all major PM (sub 2.5) components (within 6 percent). Major ions generally agree well with previous studies at the same urban locations (e.g. sulfate fractions agree within 4 percent for 8 out of 11 collocation comparisons). Enhanced anthropogenic dust fractions in large urban areas (e.g. Singapore, Kanpur, Hanoi, and Dhaka) are apparent from high Zn to Al ratios. The expected water contribution to aerosols is calculated via the hygroscopicity parameter kappa (sub v (volume)) for each filter. Mean aggregate values ranged from 0.15 (Ilorin) to 0.28 (Rehovot). The all-site parameter mean is 0.20 plus or minus 0.04. Chemical composition and water retention in each filter measurement allows inference of hourly PM (sub 2.5) at 35 percent relative humidity by merging with nephelometer measurements. These hourly PM (sub 2.5) estimates compare favourably with a beta attenuation monitor (MetOne) at the nearby US embassy in Beijing, with a coefficient of variation r squared equals 0.67 (number equals 3167), compared to r squared equals 0.62 when v (volume) was not considered. SPARTAN continues to provide an open-access database of PM (sub 2.5) compositional filter information and hourly mass collected from a global federation of instruments.

Snider, Graydon↗

Application of multi-method-multi-model inference to radiation related solid cancer excess risks models for astronaut risk assessment

The impact of including model-averaged excess radiation risks (ER) into a measure of radiation attributed decrease of survival (RADS) for the outcome all solid cancer incidence and the impact on the uncertainties is demonstrated. It is shown that RADS applying weighted model averaged ER based on AIC weights result in smaller risk estimates with narrower 95% CI than RADS using ER based on BIC weights. Further a multi-method-multi-model inference approach is introduced that allows calculating one general RADS estimate providing a weighted average risk estimate for a lunar and a Mars mission. For males the general RADS estimate is found to be 0.42% (95% CI: 0.38%; 0.45%) and for females 0.67% (95% CI: 0.59%; 0.75%) for a lunar mission and 2.45% (95% CI: 2.23%; 2.67%) for males and 3.91% (95% CI: 3.44%; 4.39%) for females for a Mars mission considering an age at exposure of 40 years and an attained age of 65 years. It is recommended to include these types of uncertainties and to include model-averaged excess risks in astronaut risk assessment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A Nonparametric Method for the Inference of Halo Occupation Distributions

The galaxy–halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a nonparametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past light cones constructed with the Santa Cruz semianalytic model. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the nonparametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modelling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

Kennedy, Jacob [Rutgers Univ., Piscataway, NJ (Uni↗

Improvement and generalization of ABCD method with Bayesian inference

To find New Physics or to refine our knowledge of the Standard Model at the LHC is an enterprise that involves many factors, such as the capabilities and the performance of the accelerator and detectors, the use and exploitation of the available information, the design of search strategies and observables, as well as the proposal of new models. We focus on the use of the information and pour our effort in re-thinking the usual data-driven ABCD method to improve it and to generalize it using Bayesian Machine Learning techniques and tools. We propose that a dataset consisting of a signal and many backgrounds is well described through a mixture model. Signal, backgrounds and their relative fractions in the sample can be well extracted by exploiting the prior knowledge and the dependence between the different observables at the event-by-event level with Bayesian tools. We show how, in contrast to the ABCD method, one can take advantage of understanding some properties of the different backgrounds and of having more than two independent observables to measure in each event. In addition, instead of regions defined through hard cuts, the Bayesian framework uses the information of continuous distribution to obtain soft-assignments of the events which are statistically more robust. To compare both methods we use a toy problem inspired by pp\to hh\to b\bar b b \bar b p p → h h → b b ‾ b b ‾ , selecting a reduced and simplified number of processes and analysing the flavor of the four jets and the invariant mass of the jet-pairs, modeled with simplified distributions. Taking advantage of all this information, and starting from a combination of biased and agnostic priors, leads us to a very good posterior once we use the Bayesian framework to exploit the data and the mutual information of the observables at the event-by-event level. We show how, in this simplified model, the Bayesian framework outperforms the ABCD method sensitivity in obtaining the signal fraction in scenarios with 1% and 0.5% true signal fractions in the dataset. We also show that the method is robust against the absence of signal. We discuss potential prospects for taking this Bayesian data-driven paradigm into more realistic scenarios.

Alvarez, Ezequiel↗

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗