Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability distribution”

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 451 records · Page 25

How Many Elements Matter?

Some studies of stars' multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and $\mathrm{log}\,g$ minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [X/H] ≲ 0.02 dex, consistent with APOGEE's reported statistical uncertainties of ~0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

79 ASTRONOMY AND ASTROPHYSICS↗

A Simple Model of Cirrus Horizontal Inhomogeneity and Cloud Fraction

A simple model of horizontal inhomogeneity and cloud fraction in cirrus clouds has been formulated on the basis that all internal horizontal inhomogeneity in the ice mixing ratio is due to variations in the cloud depth, which are assumed to be Gaussian. The use of such a model was justified by the observed relationship between the normalized variability of the ice water mixing ratio (and extinction) and the normalized variability of cloud depth. Using radar cloud depth data as input, the model reproduced well the in-cloud ice water mixing ratio histograms obtained from horizontal runs during the FIRE2 cirrus campaign. For totally overcast cases the histograms were almost Gaussian, but changed as cloud fraction decreased to exponential distributions which peaked at the lowest nonzero ice value for cloud fractions below 90%. Cloud fractions predicted by the model were always within 28% of the observed value. The predicted average ice water mixing ratios were within 34% of the observed values. This model could be used in a GCM to produce the ice mixing ratio probability distribution function and to estimate cloud fraction. It only requires basic meteorological parameters, the depth of the saturated layer and the standard deviation of cloud depth as input.

Smith, Samantha A.↗

A probabilistic fracture mechanics approach for structural reliability assessment of space flight systems

A probabilistic fracture mechanics approach for predicting the failure life distribution due to subcritical crack growth is presented. A state-of-the-art crack propagation method is used in a Monte Carlo simulation to generate a distribution of failure lives. The crack growth failure model expresses failure life as a function of stochastic parameters including environment, loads, material properties, geometry, and model specification errors. A stochastic crack growth rate model that considers the uncertainties due to scatter in the data and mode misspecification is proposed. The rationale for choosing a particular type of probability distribution for each stochastic input parameter and for specifying the distribution parameters is presented. The approach is demonstrated through a probabilistic crack growth failure analysis of a welded tube in the Space Shuttle Main Engine. A discussion of the results from this application of the methodology is given.

Sutharshana, S.↗

NO PLIF flow visualization and time-resolved temperature distributions in laser induced breakdown plumes

NO planar laser induced fluorescence (PLIF) is used to obtain images of laser-induced breakdown plasma plumes in NO-seeded nitrogen and dry air at near atmospheric pressure. Single-shot PLIF-images show that the plume development 5–50 μs after the breakdown pulse is fairly reproducible shot-to-shot, although the plume becomes increasingly stochastic on longer timescales, 100–500 μs. The stochastic behavior of the plume is quantified using probability distributions of the loci of the plume boundary. Analysis of the single-shot images indicates that the mixing of the plume with ambient gas on sub-ms time scale is insignificant. The induced flow velocity in the plume is fairly low, up to 30 m s –1 , suggesting that laser breakdowns are ineffective for mixing enhancement in high speed flows. The ensemble-averaged PLIF images indicate the evolution of the plume from an initially elongated shape to near-spherical to toroidal shape, with a subsequent radial expansion and formation of an axial jet in the center. Temperature distributions in the plume in air are obtained from the NO PLIF images, using two rotational transitions in the NO(X, v' = 0 → A, v'' = 0) band, J'' = 6.5 and 12.5 of the Q R 12 + Q 2 branch. The results indicate that the temperature in the plume remains high, above 1000 K, for approximately 100 μs, after which it decays gradually, to below 500 K at 500 μs. The residual NO fraction in the plume is ~0.1%, indicating that repetitive laser-assisted ignition may result in significant NO-generation. Furthermore, these measured temperature and velocity distributions can be used for detailed validation of kinetic models of laser-induced breakdown and assessment of their predictive capability.

42 ENGINEERING↗

Product Distribution Theory for Control of Multi-Agent Systems

Product Distribution (PD) theory is a new framework for controlling Multi-Agent Systems (MAS's). First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (probability distribution of) the joint stare of the agents. Accordingly we can consider a team game in which the shared utility is a performance measure of the behavior of the MAS. For such a scenario the game is at equilibrium - the Lagrangian is optimized - when the joint distribution of the agents optimizes the system's expected performance. One common way to find that equilibrium is to have each agent run a reinforcement learning algorithm. Here we investigate the alternative of exploiting PD theory to run gradient descent on the Lagrangian. We present computer experiments validating some of the predictions of PD theory for how best to do that gradient descent. We also demonstrate how PD theory can improve performance even when we are not allowed to rerun the MAS from different initial conditions, a requirement implicit in some previous work.

Lee, Chia Fan↗

The Sensitivity of GPz Estimates of Photo-z Posterior PDFs to Realistically Complex Training Set Imperfections

The accurate estimation of photometric redshifts is crucial to many upcoming galaxy surveys, for example, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Almost all Rubin extragalactic and cosmological science requires accurate and precise calculation of photometric redshifts; many diverse approaches to this problem are currently in the process of being developed, validated, and tested. In this work, we use the photometric redshift code GPz to examine two realistically complex training set imperfections scenarios for machine learning based photometric redshift calculation: (i) where the spectroscopic training set has a very different distribution in color–magnitude space to the test set, and (ii) where the effect of emission line confusion causes a fraction of the training spectroscopic sample to not have the true redshift. By evaluating the sensitivity of GPz to a range of increasingly severe imperfections, with a range of metrics (both of photo-z point estimates as well as posterior probability distribution functions, PDFs), we quantify the degree to which predictions get worse with higher degrees of degradation. In particular, we find that there is a substantial drop-off in photo-z quality when line-confusion goes above ~1%, and sample incompleteness below a redshift of 1.5, for an experimental setup using data from the Buzzard Flock synthetic sky catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗

Distributed Optimization

We demonstrate a new framework for analyzing and controlling distributed systems, by solving constrained optimization problems with an algorithm based on that framework. The framework is ar. information-theoretic extension of conventional full-rationality game theory to allow bounded rational agents. The associated optimization algorithm is a game in which agents control the variables of the optimization problem. They do this by jointly minimizing a Lagrangian of (the probability distribution of) their joint state. The updating of the Lagrange parameters in that Lagrangian is a form of automated annealing, one that focuses the multi-agent system on the optimal pure strategy. We present computer experiments for the k-sat constraint satisfaction problem and for unconstrained minimization of NK functions.

Macready, William↗

Space Radiation Risk Assessment for Future Lunar Missions

For lunar exploration mission design, radiation risk assessments require the understanding of future space radiation environments in support of resource management decisions, operational planning, and a go/no-go decision. The future GCR flux was estimated as a function of interplanetary deceleration potential, which was coupled with the estimated neutron monitor rate from the Climax monitor using a statistical model. A probability distribution function for solar particle event (SPE) occurrence was formed from proton fluence measurements of SPEs occurred during the past 5 solar cycles (19-23). Large proton SPEs identified from impulsive nitrate enhancements in polar ice for which the fluences are greater than 2 10(exp 9) protons/sq cm for energies greater than 30 MeV, were also combined to extend the probability calculation for high level of proton fluences. The probability with which any given proton fluence level of a SPE will be exceeded during a space mission of defined duration was then calculated. Analytic energy spectra of SPEs at different ranks of the integral fluences were constructed over broad energy ranges extending out to GeV, and representative exposure levels were analyzed at those fluences. For the development of an integrated strategy for radiation protection on lunar exploration missions, effective doses at various points inside a spacecraft were calculated with detailed geometry models representing proposed transfer vehicle and habitat concepts. Preliminary radiation risk assessments from SPE and GCR were compared for various configuration concepts of radiation shelter in exploratory-class spacecrafts.

Kim, Myung-Hee Y.↗

Effective diffusion equation in a random velocity field

The effects are studied of assumed random velocity fields on diffusion in a binary fluid. Random velocity fields can result, for example, from the high-frequency components of residual accelerations onboard spacecraft (often called g-jitter). An effective diffusion equation is derived for an average concentration which includes spatial and temporal correlations induced by the fluctuating velocity fields assumed to be Gaussianly distributed. The resulting equation becomes nonlocal, and if correlations between different components of the velocity field exist, it is also anisotropic. The simple limiting case of short correlation times is discussed and an effective diffusivity is obtained which reflects the enhanced mixing caused by the velocity fields. The results obtained in the limit of short correlation times are valid even if the probability distribution of the velocity field is not Gaussian.

Vinals, Jorge↗

Data Applicability of Heritage and New Hardware For Launch Vehicle Reliability Models

Bayesian reliability requires the development of a prior distribution to represent degree of belief about the value of a parameter (such as a component's failure rate) before system specific data become available from testing or operations. Generic failure data are often provided in reliability databases as point estimates (mean or median). A component's failure rate is considered a random variable where all possible values are represented by a probability distribution. The applicability of the generic data source is a significant source of uncertainty that affects the spread of the distribution. This presentation discusses heuristic guidelines for quantifying uncertainty due to generic data applicability when developing prior distributions mainly from reliability predictions.

Al Hassan, Mohammad↗

Reliability-Informed Life-Cycle Warranty Cost Analysis: A Case Study on a Transmission in Agricultural Equipment

In agricultural and industrial equipment, both new and remanufactured systems are often available for warranty coverage. In such cases, it may be challenging for equipment manufacturers to properly trade-off between the system reliability and the cost associated with a replacement option (e.g., replace with a new or remanufactured system). To address this problem, we present a reliability-informed life-cycle warranty cost (LCWC) analysis framework that enables equipment manufacturers to evaluate different warranty policies. These warranty policies differ in whether a new or remanufactured system is used for replacement in the case of product failure. The novelty of this LCWC analysis framework lies in its ability to incorporate real-world field reliability data into warranty policy assessment using probabilistic warranty cost models that consider multiple life cycles. First, the reliability functions for the new and remanufactured systems are built as the time-to-failure distributions that provide the best-fit to the field reliability data. Then, these reliability functions and their corresponding warranty policies are used to build the LCWC models according to the specific warranty terms. Finally, Monte Carlo simulation is used to propagate the time-to-failure uncertainty of each system, modeled by its reliability function, through each LCWC model to produce a probability distribution of the LCWC. The effectiveness of the proposed reliability-informed LCWC analysis framework is demonstrated with a real-world case study on a transmission used in some agricultural equipment.

agricultural equipment↗

Climate Impacts of Convective Cloud Microphysics in NCAR CAM5

Here we improved the treatments of convective cloud microphysics in the NCAR Community Atmosphere Model version 5.3 (CAM5.3) by 1) implementing new terminal velocity parameterizations for convective ice and snow particles, 2) adding graupel microphysics, 3) considering convective snow detrainment, and 4) enhancing rain initiation and generation rate in warm clouds. Furthermore, we evaluated the impacts of improved microphysics on simulated global climate, focusing on simulated cloud radiative forcing, graupel microphysics, convective cloud ice amount, and tropical precipitation. Compared to CAM5.3 with the default convective microphysics, the too-strong cloud shortwave radiative forcing due primarily to excessive convective cloud liquid is largely alleviated over the tropics and midlatitudes after rain initiation and generation rate is enhanced, in better agreement with the CERES-EBAF estimates. Geographic distributions of graupel occurrence are reasonably simulated over continents; whereas the graupel occurrence remains highly uncertain over the oceanic storm-track regions. When evaluated against the CloudSat–CALIPSO estimates, the overestimation of convective ice mass is alleviated with the improved convective ice microphysics, among which adding graupel microphysics and the accompanying increase in hydrometeor fall speed play the most important role. The probability distribution function (PDF) of rainfall intensity is sensitive to warm rain processes in convective clouds, and enhancement in warm rain production shifts the PDF toward heavier precipitation, which agrees better with the TRMM observations. Common biases of overestimating the light rain frequency and underestimating the heavy rain frequency in GCMs are mitigated.

54 ENVIRONMENTAL SCIENCES↗

Management of SSME hardware life utilization

Statistical and probabilistic reliability methodologies were developed for the determination of hardware life limits for the Space Shuttle Main Engine (SSME). Both methodologies require that a mathematical reliability model of the engine (system) performance be developed as a function of the reliabilities of the components and parts. The system reliability model should be developed from the Failute Modes and Effects Analysis/Critical Items List. The statistical reliability methodology establishes hardware life limits directly from the failure distributions of the components and parts obtained from statistically-designed testing. The probabilistic reliability methodology establishes hardware life limits from a decision analysis methodology which incorporates the component/part reliabilities obtained from a probabilistic structural analysis, a calibrated maintenance program, inspection techniques, and fabrication procedures. Probilistic structural analysis is recommended as a tool to prioritize upgrading of the components and parts. The Weibull probability distribution is presently being investigated by NASA/MSFC to characterize the failure distribution of the SSME hardware from a limited data base of failures.

Pauschke, J. M.↗

Design prediction for long term stress rupture service of composite pressure vessels

Extensive stress rupture studies on glass composites and Kevlar composites were conducted by the Lawrence Radiation Laboratory beginning in the late 1960's and extending to about 8 years in some cases. Some of the data from these studies published over the years were incomplete or were tainted by spurious failures, such as grip slippage. Updated data sets were defined for both fiberglass and Kevlar composite stand test specimens. These updated data are analyzed in this report by a convenient form of the bivariate Weibull distribution, to establish a consistent set of design prediction charts that may be used as a conservative basis for predicting the stress rupture life of composite pressure vessels. The updated glass composite data exhibit an invariant Weibull modulus with lifetime. The data are analyzed in terms of homologous service load (referenced to the observed median strength). The equations relating life, homologous load, and probability are given, and corresponding design prediction charts are presented. A similar approach is taken for Kevlar composites, where the updated stand data do show a turndown tendency at long life accompanied by a corresponding change (increase) of the Weibull modulus. The turndown characteristic is not present in stress rupture test data of Kevlar pressure vessels. A modification of the stress rupture equations is presented to incorporate a latent, but limited, strength drop, and design prediction charts are presented that incorporate such behavior. The methods presented utilize Cartesian plots of the probability distributions (which are a more natural display for the design engineer), based on median normalized data that are independent of statistical parameters and are readily defined for any set of test data.

Robinson, Ernest Y.↗

Probabilistic inference in very large universes

Our current favored cosmological theories allow for the striking and controversial possibility that the observable universe is just a small part of a much larger universe in which parameters that describe the effective, low-energy laws of physics vary from one region to another. The controversy is largely driven by the fact that such a “very large universe” is mostly observationally inaccessible to us, so the issue arises of how we can reasonably assess a theory that describes such a universe. In this paper, we propose a Bayesian method for theory assessment based on theory-generated probability distributions for our observations. We focus on the principles that define this method, leaving aside concerns about how, in practice, one would carry out the required calculations. (One important issue that we set aside is the measure problem.) We argue that cosmological theories can be tested by the standard method of Bayesian updating, but we need to use theoretical predictions for “first-person” probabilities—that is, probabilities that we should use for our observations, taking into account all relevant selection effects. These selection effects can vary from one observer to another and can vary with time, so, in principle, first-person probabilities are defined for each observer instant—an observer at a specific instant of time. Calculations of first-person probabilities should take into account everything that the observer believes about herself and her surroundings, which we refer to as her subjective state. If the universe is very large, a theory might predict that there are many observer instants in the same subjective state; we argue that first-person probabilities should be calculated using a principle of self-locating indifference (PSLI), the assumption that any real observer should make predictions for her future as if she were chosen randomly and uniformly from the theoretically predicted observer instants that share her subjective state. We believe the PSLI is intuitively very reasonable, but we also argue that, if the theory is correct, the use of this principle maximizes the expected fraction of observers who will make correct predictions. A further complication is that cosmological theories are not expected to fully predict the detailed properties of the universe, but rather will predict a set of possible universes, each with a probability. Different possible universes will generically have different numbers of observers. We argue that, in the calculation of first-person probabilities, the probability for each possible universe should be weighted by the number of observer instants in the specified subjective state that it contains. These issues have been controversial in the literature, so we also provide a rebuttal to the claim that principles like the PSLI involve a “selection fallacy”; a rebuttal to what we dub the principle of required certainty; an argument rejecting theories that predict a preponderance of Boltzmann brains; a rebuttal to a parable about humans and Jovians used by Hartle and Srednicki to argue that assumptions of typicality can lead to absurd consequences; and, finally, a discussion about how the use of “old evidence” can be fit into a Bayesian mold.

Azhar, Feraz [University of Notre Dame, IN (United↗

Simple Acceptance Decisions and Acceptance Risk

Acceptance risk may be estimated using cumulative probability distribution functions applied to populations of measurands and test measurement values. “Simple acceptance” is a decision rule that sets the acceptance range of a test result equal to the tolerance range specification. While acceptance risk is comprised of both consumer risk and producer risk, this paper compares the effects of simple acceptance decision rules and guardbanding decision rules on consumer risk. Consumer risk is also known as the probability of false acceptance. The terms describe the risk of accepting test results as passing when the actual values exceed specification limits. False acceptance is only possible when the true value of a measurand is out of tolerance and the test result indicates that the measurand is within tolerance. Metrologists generally have some information about the test uncertainty regarding a specific acceptance test result. Along with a general lack of knowledge regarding the measurand population parameters, the complicated interplay between risk, dispersion, and central values generally prevents calibration laboratories from fully characterizing acceptance probabilities. Organizations that model the measurand populations can reduce consumer risk by avoiding certification of items whose measurand populations are not well centered. The models in this paper present a recurring trend: guardbanding reduces nonnegligible risks of false acceptance when compared to simple acceptance. However, guardbanding is not the most effective means of acceptance risk mitigation. Systematic characterization of measurand populations can provide the information a calibration laboratory needs to reliably control acceptance risk.

42 ENGINEERING↗

A Probabilistic Mass Estimation Algorithm for a Novel 7- Channel Capacitive Sample Verification Sensor

A document describes an algorithm created to estimate the mass placed on a sample verification sensor (SVS) designed for lunar or planetary robotic sample return missions. A novel SVS measures the capacitance between a rigid bottom plate and an elastic top membrane in seven locations. As additional sample material (soil and/or small rocks) is placed on the top membrane, the deformation of the membrane increases the capacitance. The mass estimation algorithm addresses both the calibration of each SVS channel, and also addresses how to combine the capacitances read from each of the seven channels into a single mass estimate. The probabilistic approach combines the channels according to the variance observed during the training phase, and provides not only the mass estimate, but also a value for the certainty of the estimate. SVS capacitance data is collected for known masses under a wide variety of possible loading scenarios, though in all cases, the distribution of sample within the canister is expected to be approximately uniform. A capacitance-vs-mass curve is fitted to this data, and is subsequently used to determine the mass estimate for the single channel s capacitance reading during the measurement phase. This results in seven different mass estimates, one for each SVS channel. Moreover, the variance of the calibration data is used to place a Gaussian probability distribution function (pdf) around this mass estimate. To blend these seven estimates, the seven pdfs are combined into a single Gaussian distribution function, providing the final mean and variance of the estimate. This blending technique essentially takes the final estimate as an average of the estimates of the seven channels, weighted by the inverse of the channel s variance.

Wolf, Michael↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗