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 469 records · Page 26

Extending the predictive power of perturbative QCD using the principle of maximum conformality and the Bayesian analysis

In addition to the evaluation of high-order loop contributions, the precision and predictive power of perturbative QCD (pQCD) predictions depends on two important issues: (1) how to achieve a reliable, convergent fixed-order series, and (2) how to reliably estimate the contributions of unknown higher-order terms. The recursive use of renormalization group equation, together with the Principle of Maximum Conformality (PMC), eliminates the renormalization scheme-and-scale ambiguities of the conventional pQCD series. The result is a conformal, scale-invariant series of finite order which also satisfies all of the principles of the renormalization group. In this paper we propose a novel Bayesian-based approach to estimate the size of the unknown higher order contributions based on an optimized analysis of probability distributions. We show that by using the PMC conformal series, in combination with the Bayesian analysis, one can consistently achieve high degree of reliability estimates for the unknown high order terms. Thus the predictive power of pQCD can be greatly improved. We illustrate this procedure for two pQCD observables: R e+e– and R τ , which are each known up to four loops in pQCD. Numerical analyses confirm that by using the scale-independent and more convergent PMC conformal series, one can achieve reliable Bayesian probability estimates for the unknown higher-order contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Potential health risks from postulated accidents involving the Pu-238 RTG on the Ulysses solar exploration mission

Potential radiation impacts from launch of the Ulysses solar exploration experiment were evaluated using eight postulated accident scenarios. Lifetime individual dose estimates rarely exceeded 1 mrem. Most of the potential health effects would come from inhalation exposures immediately after an accident, rather than from ingestion of contaminated food or water, or from inhalation of resuspended plutonium from contaminated ground. For local Florida accidents (that is, during the first minute after launch), an average source term accident was estimated to cause a total added cancer risk of up to 0.2 deaths. For accidents at later time after launch, a worldwide cancer risk of up to three cases was calculated (with a four in a million probability). Upper bound estimates were calculated to be about 10 times higher.

Goldman, Marvin↗

Game Strategies for Entanglement Paths in Quantum Network Infrastructure

Entanglement distribution is a core function of quantum networks, and the paths used for this purpose are composed of quantum and conventional network components and are routed through physical facility sites. A game theoretic model is formulated for the defense of entanglement paths in a quantum network infrastructure by modeling the correlations and probabilities of reinforcement and failure of its components. A sum-form utility function is used to capture the cost-benefit trade-offs in reinforcing the entanglement path components to defend against their failures and attacks. Under Nash Equilibrium criteria, estimates of survival probabilities of entanglement paths are derived using the parameters and correlations of quantum, conventional, hybrid, and facility components. They provide insights into the dependencies of entanglement paths on its components, including cross-boundary effects of conventional, quantum, and facility components.

Rao, Nageswara [ORNL] (ORCID:0000000234085941)↗

Probabilistic Hazard Assessment for Tornadoes, Straight-Line Wind, and Extreme Precipitation at the Savannah River Site

Recent data sets for three meteorological phenomena with the potential to inflict damage on SRS facilities – tornadoes, straight-line winds, and heavy precipitation – are analyzed using appropriate statistical techniques to estimate the occurrence probabilities for these events in the future. Summaries of the results for DOE-mandated return periods and comparisons to similar calculations performed in 2013 by Werth et al. (W2013) are given. Using tornado statistics for i) the combined states of Georgia and South Carolina, and ii) a 2⁰ square area surrounding SRS, we calculated the probability per year of any location at SRS being struck by a tornado (the ‘strike’ probability) and the probability that any point will experience winds above set thresholds. The strike probability was calculated to be 7.04E-4 (1 chance in 1420) per year and tornadic wind speeds for DOE mandated return periods of 50,000 years (corresponding to wind design category 3 (WDC-3), and 125,000 years (meeting WDC-4) (USDOE, 2016) were estimated to be 132 mph and 147 mph, respectively. By contrast, default tornado wind speeds taken from ANSI/ANS-2.3-2011 are somewhat higher: 161 mph for return periods of 50,000 years and 173 mph every 125,000 years (ANS, 2011). Although the ANS and the SRS evaluation used the same basic model (Ramsdell and Rishel, 2007), the region defined in ANS 2.3 that encompasses the SRS also includes areas of the Great Plains and lower Midwest, regions with much higher occurrence frequencies of strong tornadoes. The SRS straight-line wind values associated with various return periods were calculated by fitting existing wind data to a GEV1 distribution and extrapolating the values for any return period from the tail of that function. For the DOE mandated return periods, we expect straight-line winds of 117 mph every 2500 years (the required WDC-3 standard) and 125 mph every 6250 years (WDC-4) at any point within the SRS. These values are similar to those from the ANS-2.3-2011 report, which has wind speeds of 125mph and 133 mph for return periods of 2500 years and 6250 years, respectively. For extreme precipitation, we compared the fits of two different theoretical extreme-value distributions and applied the one that fit the data best for each of several accumulation periods. The DOE mandated 6-hr accumulated rainfall for return periods of 10,000 years (corresponding to precipitation design category 3 (PDC-3) and 25,000 years (PDC-4) were estimated as 9.1 inches and 10.1 inches, respectively. For the 24-hr rainfall return periods of 10,000 years and 25,000 years, total rainfall estimates were 12.02 inches and 13.17 inches, respectively, higher than comparable values provided in the W2013 report.

54 ENVIRONMENTAL SCIENCES↗

Rare Events via Cross-Entropy Population Monte Carlo

Rare events are events that happen with very low frequency. Estimating rare event probabilities using Monte Carlo techniques is computationally expensive, often to the point of intractability, and special methods are required. Importance sampling (IS) is a well known technique that uses a proposal distribution in place of a target distribution to lower the variance of the estimator. Key to the success of IS methods is the choice of a proposal distribution, or the parameters governing the distribution. Adaptive importance sampling improves the parameters of a family or population of proposal distributions iteratively through trials. We present a novel cross-entropy population Monte Carlo algorithm, which adapts the parameters of proposals through the cross-entropy method. The proposed method stands apart from previous work in that we are not optimizing a mixture distribution. Instead, we leverage deterministic mixture weights and optimize the distributions individually through a reinterpretation of the typical derivation of the cross-entropy method. Demonstrations on rare event examples show that the algorithm can outperform existing resampling based population Monte Carlo methods, especially for higher-dimensional problems. Finally, we also demonstrate efficacy on a conjunction analysis problem.

97 MATHEMATICS AND COMPUTING↗

Estimation and classification by sigmoids based on mutual information

An estimate of the probability density function of a random vector is obtained by maximizing the mutual information between the input and the output of a feedforward network of sigmoidal units with respect to the input weights. Classification problems can be solved by selecting the class associated with the maximal estimated density. Newton's s method, applied to an estimated density, yields a recursive maximum likelihood estimator, consisting of a single internal layer of sigmoids, for a random variable or a random sequence. Applications to the diamond classification and to the prediction of a sun-spot process are demonstrated.

Baram, Yoram↗

Acreage estimation, feature selection, and signature extension dependent upon the maximum likelihood decision rule

A maximum likelihood estimation technique is used for the analysis of agricultural remote sensor data. The m-class probability of misclassification is estimated using unlabeled test samples and labeled training samples. A bound on the variance of a proposed unbiased estimator of the m-class probability of error is derived. The particular case in which each class density is assumed to be a mixture of multivariate normal densities is considered. The extension of spectral signatures in space and time is discussed.

Quirein, J. A.↗

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

Model ensembles are effective tools for estimating prediction uncertainty in deep learning atomistic force fields. However, their widespread adoption is hindered by high computational costs and overconfident error estimates. In this work, we address these challenges by leveraging distributions of per-sample errors obtained during training and employing a distance-based similarity search in the model latent space. Our method, which we call LTAU (Loss Trajectory Analysis for Uncertainty), efficiently estimates the full probability distribution function of errors for any test point using the logged training errors, achieving speeds that are 2–3 orders of magnitudes faster than typical ensemble methods and allowing it to be used for tasks where training or evaluating multiple models would be infeasible. We apply LTAU towards estimating parametric uncertainty in atomistic force fields (LTAU-FF), demonstrating that it produces well-calibrated confidence intervals and predicts errors that correlate strongly with the true errors for data near the training domain. Furthermore, we show that the errors predicted by LTAU-FF can be used in practical applications for detecting out-of-domain data, tuning model performance, and predicting failure during simulations. We believe that LTAU will be a valuable tool for uncertainty quantification in atomistic force fields and is a promising method that should be further explored in other domains of machine learning.

97 MATHEMATICS AND COMPUTING↗

Probabilistic risk analysis of flying the space shuttle with and without fuel turbine discharge temperature redline protection

An exact mathematical model is presented of the Space Shuttle Main Engine computer voting logic in the presence of High Pressure Fuel Turbine (HPFT) overtemperature events and fuel turbine temperature sensor failures. The model provides estimates of the probability of erroneous engine shutdown and the probability of not detecting a HPFT overtemperature event. Because it is believed that the likelihood of sensor failures and overtemperature events in the HPFT greatly overshadows those in the High Pressure Oxygen Turbine (HPOT), this modeling effort focused on the HPFT. However, because the redline protection logic is the same for both turbines, estimation of the model parameters using relevant HPOT data would provide estimates of the risk of erroneous engine shutdown and an undetected overtemperature event in the HPOT. Because of the complexity of the model, it was necessary to program the solution which thus makes it feasible to accommodate a changing data base. This is considered to be of great interest because of the subjective nature in determining the relevancy of certain failures and the fact that the data base is constantly changing as a result of the frequent engine tests.

Howell, Leonard↗

Hidden Markov models for fault detection in dynamic systems

The invention is a system failure monitoring method and apparatus which learns the symptom-fault mapping directly from training data. The invention first estimates the state of the system at discrete intervals in time. A feature vector x of dimension k is estimated from sets of successive windows of sensor data. A pattern recognition component then models the instantaneous estimate of the posterior class probability given the features, p(w(sub i) (vertical bar)/x), 1 less than or equal to i isless than or equal to m. Finally, a hidden Markov model is used to take advantage of temporal context and estimate class probabilities conditioned on recent past history. In this hierarchical pattern of information flow, the time series data is transformed and mapped into a categorical representation (the fault classes) and integrated over time to enable robust decision-making.

Smyth, Padhraic J.↗

Hidden Markov models for fault detection in dynamic systems

The invention is a system failure monitoring method and apparatus which learns the symptom-fault mapping directly from training data. The invention first estimates the state of the system at discrete intervals in time. A feature vector x of dimension k is estimated from sets of successive windows of sensor data. A pattern recognition component then models the instantaneous estimate of the posterior class probability given the features, p(w(sub i) perpendicular to x), 1 less than or equal to i is less than or equal to m. Finally, a hidden Markov model is used to take advantage of temporal context and estimate class probabilities conditioned on recent past history. In this hierarchical pattern of information flow, the time series data is transformed and mapped into a categorical representation (the fault classes) and integrated over time to enable robust decision-making.

Smyth, Padhraic J.↗

Studies on possible propagation of microbial contamination in planetary clouds

One of the key parameters in estimation of the probability of contamintion of the outer planets (Jupiter, Saturn, Uranus, etc.) is the probability of growth (Pg) of terrestrial microorganisms on or near these planets. For example, Jupiter appears to have an atmosphere in which some microbial species could metabolize and propagate. This study includes investigation of the likelihood of metabolism and propagation of microbes suspended in dynamic atmospheres. It is directed toward providing experimental information needed to aid in rational estimation of Pg for these outer planets. Current work is directed at demonstration of aerial metabolism under near optimal conditions and tests of propagation in simulated Jovian atmospheres.

Dimmick, R. L.↗

Sea-Level Rise: Projections for Maryland 2018

In fulfillment of requirements of the Maryland Commission on Climate Change Act of 2015, this report provides updated projections of the amount of sea-level rise relative to Maryland coastal lands that is expected into the next century. These projections represent the consensus of an Expert Group drawn from the Mid-Atlantic region. The framework for these projections is explicitly tied to the projections of global sea-level rise included in the Intergovernmental Panel on Climate Change Fifth Assessment (2014) and incorporates regional factors such as subsidence, distance from melting glaciers and polar ice sheets, and ocean currents. The probability distribution of estimates of relative sea-level rise from the baseline year of 2000 are provided over time and, after 2050, for three different greenhouse gas emissions pathways: Growing Emissions (RCP8.5), Stabilized Emissions (RCP4.5), and meeting the Paris Agreement (RCP2.6). This framework has been recently used in developing relative sea-level rise projections for California, Oregon, Washington, New Jersey, and Delaware as well as several metropolitan areas. The Likely range (66% probability) of the relative rise of mean sea level expected in Maryland between 2000 and 2050 is 0.8 to 1.6 feet, with about a one-in-twenty chance it could exceed 2.0 feet and about a one-in-one hundred chance it could exceed 2.3 feet. Later this century, rates of sea-level rise increasingly depend on the future pathway of global emissions of greenhouse gases during the next sixty years. If emissions continue to grow well into the second half of the 21st century, the Likely range of sea-level rise experienced in Maryland is 2.0 to 4.2 feet over this century, two to four times the sea-level rise experienced during the 20th century. Moreover, there is a one-in-twenty chance that it could exceed 5.2 feet. If, on the other hand, global society were able to bring net greenhouse gas emissions to zero in time to meet the goals of the Paris Climate Agreement and reduce emissions sufficient to limit the increase in global mean temperature to less than 2°Celsius over pre-industrial levels, the Likely range for 2100 is 1.2 to 3.0 feet, with a one-in-twenty chance that it would exceed 3.7 feet. The difference in sea-level rise between these contrasting scenarios would diverge even more during the next century, with the failure to reduce emissions in the near term resulting in much greater sea-level rise 100 years from now. Moreover, recent research suggests that, without imminent and substantial reductions in greenhouse gas emissions, the loss of polar ice sheets-and thus the rate of sea-level rise-may be more rapid than assumed in these projections, particularly under the Growing Emissions scenario. These probabilistic sea-level rise projections can and should be used in planning and regulation, infrastructure siting and design, estimation of changes in tidal range and storm surge, developing inundation mapping tools, and adaptation strategies for high-tide flooding and saltwater intrusion.

Boesch, D. F.↗

Radio emission from the outer parts of stellar atmospheres (mantles) at centimeter wavelengths and the winds of early-type stars

It is noted that when magnetic lines of force are present in the mantle of an early-type star, gyroresonance radiation at radio wavelengths may occur at detectable levels in addition to bremsstrahlung. An expression is derived for the effective value of the ratio of the emissivity in gyroresonance radiation to that in bremsstrahlung, and some numerical values are estimated using values of the parameters which may be typical for conditions in the mantles of early-type stars. The volume emissivity in gyroresonance radiation may exceed that for bremsstrahlung by a large factor. The implications of this result for the interpretation of the radio fluxes from early-type stars are discussed. It is shown that the radio-frequency radiation emerging from a hot star probably consists of gyroresonance radiation from an array of magnetic loops as well as bremsstrahlung from a low-density wind. It is argued that it is probable that the estimates of rate of mass loss made ignoring the possibility of gyroresonance radiation may be too large by a factor of at least 10.

Underhill, A. B.↗

GPS Spoofing Mitigation and Timing Risk Analysis in Networked Phasor Measurement Units via Stochastic Reachability

To address phasor measurement unit (PMU) vulnerability to spoofing, we propose the use of a set-valued state estimation technique known as stochastic reachability (SR)-based distributed Kalman filter (DKF) that computes secure global positioning system (GPS) timing across a network of receivers. Utilizing SR, we estimate not only GPS time but also its stochastic reachable set, which is parameterized by probabilistic zonotope (p-Zonotope). While requiring known measurement error bounds in only non-spoofed conditions, we designed a two-tiered approach. We first performed measurement-level spoofing mitigation via deviation of a measurement innovation from its expected p-Zonotope. We then performed state-level timing risk analysis via a determination of the intersection probability of the estimated p-Zonotope with an unsafe set that violates IEEE C37.118.1a-2014 standards. Finally, we validated our SR-DKF algorithm by subjecting it to a simulated receiver network to coordinate signal-level spoofing. We demonstrate improved timing accuracy and successful spoofing mitigation via the use of our SR-DKF algorithm. We also validated the robustness of the estimated timing risk as the number of receivers were varied.

47 OTHER INSTRUMENTATION↗

Stochastic robustness

To carry out stochastic robustness analysis, an expected probability distribution is assigned to each uncertain parameter in the system. The Monte Carlo analysis proceeds by repeatedly assigning shaped random values to each plant parameter, evaluating the stability of performance metric, and performing the binary classification (stable/unstable, etc.). If the system is stable, the state response to a unit disturbance impulse can be propagated to establish whether the response would violate settling time envelopes and whether peak actuator use would violate predetermined maximums. The final estimates of the probability of each form of unacceptable behavior are found by dividing the number of cases in which the overall system had that form of unacceptability by the number of cases run. Stability robustness can be portrayed graphically using the stochastic root locus and by using histograms of parameter values found in the unacceptable cases.

Marrison, C.↗

Chance and necessity in the selection of nucleic acid catalysts

In Tom Stoppard's famous play [Rosencrantz and Guildenstern are Dead], the ill-fated heroes toss a coin 101 times. The first 100 times they do so the coin lands heads up. The chance of this happening is approximately 1 in 10(30), a sequence of events so rare that one might argue that it could only happen in such a delightful fiction. Similarly rare events, however, may underlie the origins of biological catalysis. What is the probability that an RNA, DNA, or protein molecule of a given random sequence will display a particular catalytic activity? The answer to this question determines whether a collection of such sequences, such as might result from prebiotic chemistry on the early earth, is extremely likely or unlikely to contain catalytically active molecules, and hence whether the origin of life itself is a virtually inevitable consequence of chemical laws or merely a bizarre fluke. The fact that a priori estimates of this probability, given by otherwise informed chemists and biologists, ranged from 10(-5) to 10(-50), inspired us to begin to address the question experimentally. As it turns out, the chance that a given random sequence RNA molecule will be able to catalyze an RNA polymerase-like phosphoryl transfer reaction is close to 1 in 10(13), rare enough, to be sure, but nevertheless in a range that is comfortably accessible by experiment. It is the purpose of this Account to describe the recent advances in combinatorial biochemistry that have made it possible for us to explore the abundance and diversity of catalysts existing in nucleic acid sequence space.

Review, Tutorial↗

Joint Modeling of Wind Speed and Wind Direction Through a Conditional Approach

Atmospheric near surface wind speed and wind direction play an important role in many applications, ranging from air quality modeling, building design, wind turbine placement to climate change research. It is therefore crucial to accurately estimate the joint probability distribution of wind speed and direction. In this work, we develop a conditional approach to model these two variables, where the joint distribution is decomposed into the product of the marginal distribution of wind direction and the conditional distribution of wind speed given wind direction. To accommodate the circular nature of wind direction, a von Mises mixture model is used; the conditional wind speed distribution is modeled as a directional dependent Weibull distribution via a two-stage estimation procedure, consisting of a directional binned Weibull parameter estimation, followed by a harmonic regression to estimate the dependence of the Weibull parameters on wind direction. A Monte Carlo simulation study indicates that our method outperforms two other approaches in estimation efficiency: one that utilizes periodic spline quantile regression and another that generates data from the commonly used Abe-Ley distribution for cylindrical data. We illustrate our method by using the output from a regional climate model to investigate how the joint distribution of wind speed and direction may change under some future climate scenarios. Our method indicates significant changes in the variation of wind speed with respect to some directions.

17 WIND ENERGY↗