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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1984

An urnful of blinding functions

There is a fundamental connection between the mathematical theory of discrete probability distributions and the parametric curves and surfaces of computer-aided geometric design. It is no accident that the blending functions of Bezier curves and surfaces have an obvious probabilistic interpretation, nor is it a coincidence that the normalized uniform B-spline basis functions also model a simple stochastic process. The link between probability and geometry, and how to exploit simple probabilistic arguments to derive many of the classical geometric properties of the parametric curves and surfaces currently in vogue in computer-aided geometric design is discussed. This probabilistic approach is also used to introduce many new types of curves and surfaces, and it is demonstrated how probability theory can be used to simplify, unify, and generalize many well-known results.

Goldman, R. N.↗

NASA NTRS1998

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.↗

NASA NTRS1992

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.↗

NASA NTRS2004

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↗

NASA NTRS2005

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↗

NASA NTRS2007

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.↗

NASA NTRS1992

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↗

NASA NTRS2015

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↗

NASA NTRS1986

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.↗

NASA NTRS1992

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.↗

NASA NTRS2012

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↗

NASA NTRS1997

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗

NASA NTRS2004

Improving Project Management with Simulation and Completion Distribution Functions

Despite the critical importance of project completion timeliness, management practices in place today remain inadequate for addressing the persistent problem of project completion tardiness. A major culprit in late projects is uncertainty, which most, if not all, projects are inherently subject to. This uncertainty resides in the estimates for activity durations, the occurrence of unplanned and unforeseen events, and the availability of critical resources. In response to this problem, this research developed a comprehensive simulation based methodology for conducting quantitative project completion time risk analysis. It is called the Project Assessment by Simulation Technique (PAST). This new tool enables project stakeholders to visualize uncertainty or risk, i.e. the likelihood of their project completing late and the magnitude of the lateness, by providing them with a completion time distribution function of their projects. Discrete event simulation is used within PAST to determine the completion distribution function for the project of interest. The simulation is populated with both deterministic and stochastic elements. The deterministic inputs include planned project activities, precedence requirements, and resource requirements. The stochastic inputs include activity duration growth distributions, probabilities for events that can impact the project, and other dynamic constraints that may be placed upon project activities and milestones. These stochastic inputs are based upon past data from similar projects. The time for an entity to complete the simulation network, subject to both the deterministic and stochastic factors, represents the time to complete the project. Repeating the simulation hundreds or thousands of times allows one to create the project completion distribution function. The Project Assessment by Simulation Technique was demonstrated to be effective for the on-going NASA project to assemble the International Space Station. Approximately $500 million per month is being spent on this project, which is scheduled to complete by 2010. NASA project stakeholders participated in determining and managing completion distribution functions produced from PAST. The first result was that project stakeholders improved project completion risk awareness. Secondly, using PAST, mitigation options were analyzed to improve project completion performance and reduce total project cost.

Cates, Grant R.↗

NASA NTRS2011

The Impacts of an Observationally-Based Cloud Fraction and Condensate Overlap Parameterization on a GCM's Cloud Radiative Effect

It has been shown that the details of how cloud fraction overlap is treated in GCMs has substantial impact on shortwave and longwave fluxes. Because cloud condensate is also horizontally heterogeneous at GCM grid scales, another aspect of cloud overlap should in principle also be assessed, namely the vertical overlap of hydrometeor distributions. This type of overlap is usually examined in terms of rank correlations, i.e., linear correlations between hydrometeor amount ranks of the overlapping parts of cloud layers at specific separation distances. The cloud fraction overlap parameter and the rank correlation of hydrometeor amounts can be both expressed as inverse exponential functions of separation distance characterized by their respective decorrelation lengths (e-folding distances). Larger decorrelation lengths mean that hydrometeor fractions and probability distribution functions have high levels of vertical alignment. An analysis of CloudSat and CALIPSO data reveals that the two aspects of cloud overlap are related and their respective decorrelation lengths have a distinct dependence on latitude that can be parameterized and included in a GCM. In our presentation we will contrast the Cloud Radiative Effect (CRE) of the GEOS-5 atmospheric GCM (AGCM) when the observationally-based parameterization of decorrelation lengths is used to represent overlap versus the simpler cases of maximum-random overlap and globally constant decorrelation lengths. The effects of specific overlap representations will be examined for both diagnostic and interactive radiation runs in GEOS-5 and comparisons will be made with observed CREs from CERES and CloudSat (2B-FLXHR product). Since the radiative effects of overlap depend on the cloud property distributions of the AGCM, the availability of two different cloud schemes in GEOS-5 will give us the opportunity to assess a wide range of potential cloud overlap consequences on the model's climate.

Oreopoulos, Lazaros↗

NASA NTRS2018

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.↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

Gyekenyesi, J. P.↗

NASA NTRS1991

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.↗