Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability estimates”

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 289 records · Page 16

Machine Learning Predictions of Transition Probabilities in Atomic Spectra

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

74 ATOMIC AND MOLECULAR PHYSICS↗

Probabilistic Assessment of Radiation Risk for Astronauts in Space Missions

Accurate predictions of the health risks to astronauts from space radiation exposure are necessary for enabling future lunar and Mars missions. Space radiation consists of solar particle events (SPEs), comprised largely of medium energy protons, (less than 100 MeV); and galactic cosmic rays (GCR), which include protons and heavy ions of higher energies. While the expected frequency of SPEs is strongly influenced by the solar activity cycle, SPE occurrences themselves are random in nature. A solar modulation model has been developed for the temporal characterization of the GCR environment, which is represented by the deceleration potential, phi. The risk of radiation exposure from SPEs during extra-vehicular activities (EVAs) or in lightly shielded vehicles is a major concern for radiation protection, including determining the shielding and operational requirements for astronauts and hardware. To support the probabilistic risk assessment for EVAs, which would be up to 15% of crew time on lunar missions, we estimated the probability of SPE occurrence as a function of time within a solar cycle using a nonhomogeneous Poisson model to fit the historical database of measurements of protons with energy > 30 MeV, (phi)30. The resultant organ doses and dose equivalents, as well as effective whole body doses for acute and cancer risk estimations are analyzed for a conceptual habitat module and a lunar rover during defined space mission periods. This probabilistic approach to radiation risk assessment from SPE and GCR is in support of mission design and operational planning to manage radiation risks for space exploration.

Kim, Myung-Hee↗

Tin Whisker Electrical Short Circuit Characteristics

Existing risk simulations make the assumption that when a free tin whisker has bridged two adjacent exposed electrical conductors, the result is an electrical short circuit. This conservative assumption is made because shorting is a random event that has an unknown probability associated with it. Note however that due to contact resistance electrical shorts may not occur at lower voltage levels. In our first article we developed an empirical probability model for tin whisker shorting. In this paper, we develop a more comprehensive empirical model using a refined experiment with a larger sample size, in which we studied the effect of varying voltage on the breakdown of the contact resistance which leads to a short circuit. From the resulting data we estimated the probability distribution of an electrical short, as a function of voltage. In addition, the unexpected polycrystalline structure seen in the focused ion beam (FIB) cross section in the first experiment was confirmed in this experiment using transmission electron microscopy (TEM). The FIB was also used to cross section two card guides to facilitate the measurement of the grain size of each card guide's tin plating to determine its finish.

Courey, Karim J.↗

Efficient Calibration of Expensive Computational Models

Accounting for uncertainty when calibrating expensive computational models is a common challenge faced by scientists and engineers. Often Bayesian techniques are adopted to estimate a probability density function over the model parameters given noisy empirical data. The methods used to perform this type of probabilistic calibration are computationally prohibitive in that they require a large number of evaluations of the expensive model. In these cases, surrogate modeling -- that is, using a fast-to-evaluate, lower fidelity stand-in for the original computational model -- may be the only option to alleviate this computational burden. However, the upfront cost of generating training data to build a surrogate model can itself be expensive. As such, it is important to be judicious when selecting training points at which the full-fidelity model is evaluated. Here, an active learning approach is proposed that enables efficient selection of training points using approximate samples of the calibrated parameter probability density function. In this way, the training points can be concentrated in regions where the calibration algorithm requires high model accuracy.

active learning↗

A study of solar flare soft X rays and their relation to particle events observed with IMP 8

Data from the Johns Hopkins University Applied Physics Laboratory charged particle measurement experiment aboard IMP 8 was used to study solar flare events and find correlations between soft X rays (4-17 A) and particle events. The results show that a greater probability exists of observing a proton event if the observed soft X ray flare has a relatively large decay time. It is also found that the H-alpha brightening area is a good indicator of the strength of the soft X ray burst. Results also show that the strength of the soft X ray burst can be used to estimate a probability that a proton event will occur, along with the strength of that particle event.

Nonnast, J. H.↗

Ultimate limits to error probabilities for ionospheric models based on solar geophysical indices and how these compare with the state of the art

An ideal model based on a given set of geophysical indices is defined as a model that provides a least squares fit to the data set as a function of the indices considered. Satellite measurements of electron content for three stations at different magnetic latitudes were used to provide such data sets which were each fitted to the geophysical indices. The magnitude of the difference between the measured value and the derived equation for the data set was used to estimate the probability of making an error greater than a given magnitude for such an ideal model. Atmospheric Explorer C data is used to examine the causes of the fluctuations and suggestions are made about how real improvements can be made in ionospheric forecasting ability. Joule heating inputs in the auroral electrojets are related to the AL and AU magnetic indices. Magnetic indices based on the time integral of the energy deposited in the electrojets are proposed for modeling processes affected by auroral zone heating.

Nisbet, J. S.↗

Motivating the sure bounds

Motivation is provided for a theorem that provides upper and lower bounds for the reliability of reconfigurable digital control systems. The reliability goals for these systems are too high to be established by natural life testing, which means the probability of system failure must be computed from mathematical models that capture the essential elements of fault occurence and system fault recovery. The upper and lower bound theorem shows that system recovery can be adequately described by its first two moments, provided component failure rate is low and system recovery is fast. This result greatly simplifies both the fault injection experiments that study system recovery and the numerical computations that estimate the probability of system failure from a mathematical model.

White, Allan L.↗

A methodology for probabilistic prediction of structural failures of launch vehicle propulsion systems

A quantitative methodology for estimating the probability of occurrence of critical structural failure modes of launch vehicle propulsion systems during mission operation is presented. In this approach, engineering analysis models that characterize specific failure modes based on the physics and mechanics of the failure phenomenon are used in a prescribed probabalistic structure in which quantified uncertainty in the information used in the analysis and quantified uncertainty of the engineering analysis models themselves are both considered in estimating failure risk. The methodology is applicable to high cycle fatigue, low cycle fatigue, flaw propagation stress rupture, and other failure modes of launch vehicle propulsion systems and other space-flight systems. The probabalistic failure assessment presented provides the capability of performing quantitative analyses to rank options to improve reliability, thereby enabling limited financial resources of development and improvement programs to be more effectively allocated.

Moore, N.↗

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗

Turbine Damage Probability

Damage probability maps for offshore wind turbines exposed to tropical cyclones (TCs) under both historical and future climate scenarios along the U.S. Atlantic and Gulf Coasts are presented in this dataset. TCs are generated using the Risk Analysis Framework for Tropical Cyclones (RAFT), forced by CMIP6 historical and future global climate simulations. Maximum wind speeds for 20- and 50-year TCs are processed through a fragility function specific to offshore wind (OSW) turbines in order to estimate the probability of damage – specifically yielding and buckling – based on wind speed intensity.

17 WIND ENERGY↗

Statistical Study of Aircraft Icing Probabilities at the 700- and 500- Millibar Levels over Ocean Areas in the Northern Hemisphere

A statistical study is made of icing data reported from weather reconnaissance aircraft flown by Air Weather Service (USAF). The weather missions studied were flown at fixed flight levels of 500 millibars (18,000 ft) and 700 millibars (10,000 ft) over wide areas of the Pacific, Atlantic, and Arctic Oceans. This report is presented as part of a program conducted by the NACA to obtain extensive icing statistics relevant to aircraft design and operation. The thousands of in-flight observations recorded over a 2- to 4-year period provide reliable statistics on icing encounters for the specific areas, altitudes, and seasons included in the data. The relative frequencies of icing occurrence are presented, together with the estimated icing probabilities and the relation of these probabilities to the frequencies of flight in clouds and cloud temperatures. The results show that aircraft operators can expect icing probabilities to vary widely throughout the year from near zero in the cold Arctic areas in winter up to 7 percent in areas where greater cloudiness and warmer temperatures prevail. The data also reveal a general tendency of colder cloud temperatures to reduce the probability of icing in equally cloudy conditions.

Perkins, Porter J.↗

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗

J-PLUS: Spectral evolution of white dwarfs by PDF analysis

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), which provides 12 photometric optical passbands over 2176 deg 2 . We analyzed 5926 white dwarfs with r ≤ 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen- and helium-dominated atmospheres. We estimated the probability distribution functions for effective temperature (T eff ), surface gravity, parallax, and composition; and the probability of having a H-dominated atmosphere (p H ) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as a reference, and derived a self-consistent prior for the atmospheric composition as a function of T eff . We described the fraction of white dwarfs with a He-dominated atmosphere (f He ) with a linear function of the effective temperature at 5000 < T eff < 30 000 K. We find f He = 0.24 ± 0.01 at T eff = 10 000 K, a change rate along the cooling sequence of 0.14 ± 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 ± 0.02 at the high-temperature end. We tested the obtained p H by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 M ⊙ , and T eff > 6000 K. The result for H-dominated white dwarfs agrees with previous studies, with a dominant M = 0.59 M ⊙ peak and the presence of an excess at M ~ 0.8 M ⊙ . This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provide a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 ± 3% increase in the fraction of He-dominated white dwarfs from T eff = 20 000 K to T eff = 5000 K.

79 ASTRONOMY AND ASTROPHYSICS↗

Geoeffectiveness (D (sub st) and K (sub p)) of Interplanetary Coronal Mass Ejections During 1995-2009 and Implications for Storm Forecasting

We summarize the geoeffectiveness (based on the Dst and Kp indices) of the more than 300 interplanetary coronal mass ejections (ICMEs) that passed the Earth during 1996-2009, encompassing solar cycle 23. We subsequently estimate the probability that an ICME will generate geomagnetic activity that exceeds certain thresholds of Dst or Kp, including the NOAA "G" storm scale, based on maximum values of the southward magnetic field component (Bs), the solar wind speed (V), and the y component (Ey) of the solar wind convective electric field E = -V x B, in the ICME or sheath ahead of the ICME. Consistent with previous studies, the geoeffectiveness of an ICME is correlated with Bs or Ey approx.= VBs in the ICME or sheath, indicating that observations from a solar wind monitor upstream of the Earth are likely to provide the most reliable forecasts of the activity associated with an approaching ICME. There is also a general increase in geoeffectiveness with ICME speed, though the overall event-to-event correlation is weaker than for Bs and Ey. Nevertheless, using these results, we suggest that the speed of an ICME approaching the Earth inferred, for example, from routine remote sensing by coronagraphs on spacecraft well separated from the Earth or by all-sky imagers, could be used to estimate the likely geoeffectiveness of the ICME (our "comprehensive" ICME database provides a proxy for ICMEs identified in this way) with a longer lead time than may be possible using an upstream monitor

Richardson, I. G.↗

An Information Theory Approach on Deciding Spectroscopic Follow-ups

Classification and characterization of variable phenomena and transient phenomena are critical for astrophysics and cosmology. These objects are commonly studied using photometric time series or spectroscopic data. Given that many ongoing and future surveys are conducted in a time domain, and given that adding spectra provides further insights but requires more observational resources, it would be valuable to know which objects we should prioritize to have a spectrum in addition to a time series. We propose a methodology in a probabilistic setting that determines a priori which objects are worth taking a spectrum of to obtain better insights, where we focus on the insight of the type of the object (classification). Objects for which we query their spectrum are reclassified using their full spectral information. We first train two classifiers, one that uses photometric data and another that uses photometric and spectroscopic data together. Then for each photometric object we estimate the probability of each possible spectrum outcome. We combine these models in various probabilistic frameworks (strategies), which are used to guide the selection of follow-up observations. The best strategy depends on the intended use, whether it is obtaining more confidence or accuracy. For a given number of candidate objects (127, equal to 5% of the data set) for taking spectra, we improve the class prediction accuracy by 37% as opposed to 20% of a non-naive (non-random) best-baseline strategy. Our approach provides a general framework for follow-up strategies and can be extended beyond classification to include other forms of follow-ups beyond spectroscopy.

79 ASTRONOMY AND ASTROPHYSICS↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

A dataset of eco-evidence tools to inform early-stage environmental impact assessments of hydropower development

The datasets described herein provide the foundation for a decision support prototype (DSP) toolkit aimed at assisting stakeholders in determining evidence of which aspects of river ecosystems have been impacted by hydropower. The DSP toolkit and its application are presented and described in the article “Evidence-based indicator approach to guide preliminary environmental impact assessments of hydropower development” [1]. Development of the DSP and the output for decision support centralize around 42 river function indicators describing the dimensionality of river ecosystems through six main categories: biota and biodiversity, water quality, hydrology, geomorphology, land cover, and river connectivity. Three main tools are represented in the DSP: A science-based questionnaire (SBQ), an environmental envelope model (EEM), and a river function linkage assessment tool (RFLAT). The SBQ is a structured survey-style questionnaire whose objective is to provide evidence of which indicators have been impacted by hydropower. Based on a global literature review, 140 questions were developed from general hypotheses regarding the impacts of dams on rivers. The EEM is a model to predict the likelihood of hydropower impacting indicators based on a several variables. The intended use of the EEM is for situations of new hydropower development where results of the SBQ are incomplete or highly uncertain. The EEM was developed through the compilation of a dataset containing attributes of dams, reservoirs, and geospatial information on environmental concerns, which was combined with data on ecological indicators documented at those sites through literature review. The model operates through 247 “envelopes” and weighting factors, representing the individual effect of each variable on each indicator, all available through spreadsheets. Finally, the RFLAT is a tool to examine causal relationships amongst indicators. Inter-indicator relationships were hypothesized based on literature review and summarized into node and edge datasets to represent the structure of a graphical network. Bayes theorem was used estimate conditional probabilities of inter-indicator relationships based on the output of the SBQ. Nodes and edges were imported into R programming environment to visualize ecological indicator networks. The datasets can be expanded upon and enriched with more detailed questions for the SBQ, building upon the EEM with to develop more sophisticated models, and identifying new relationships for the RFALT. Additionally, once the tools are applied to numerous hydropower developments, the output of the tools (e.g. evidence of impacted indicators) becomes a very useful dataset for meta-analyses of hydropower impacts.

13 HYDRO ENERGY↗

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗