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 1,729 records · Page 96

Statistical Evaluation and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield/ultimate strength and high- cycle fatigue capability. This Technical Publication examines the cumulative distribution percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematica to calculate the combined value corresponding to any desired percentile is then presented along with a curve tit to this value. Another Excel macro that calculates the combination using Monte Carlo simulation is shown. Unlike the traditional techniques. these methods quantify the calculated load value with a consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Additionally, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can substantially lower the design loading without losing any of the identified structural reliability.

Brown, A. M.↗

Adaptive Optics Communications Performance Analysis

The performance improvement obtained through the use of adaptive optics for deep-space communications in the presence of atmospheric turbulence is analyzed. Using simulated focal-plane signal-intensity distributions, uncoded pulse-position modulation (PPM) bit-error probabilities are calculated assuming the use of an adaptive focal-plane detector array as well as an adaptively sized single detector. It is demonstrated that current practical adaptive optics systems can yield performance gains over an uncompensated system ranging from approximately 1 dB to 6 dB depending upon the PPM order and background radiation level.

Srinivasan, M.↗

Analysis of Terrestrial Interference Protection from UAS CNPC Satellite Transmitters

Unmanned aircraft (UA) are projected to have a major impact on future aviation. Larger UA operating at altitudes above 3000 feet will require at least occasional access to non-segregated, that is, controlled airspace. In order for unmanned aircraft to be integrated into the airspace and operate with other commercial aircraft, a very reliable command and control (a. k. a. control and non-payload communications, (CNPC)) link is required. For operations covering large distances or over remote locations, a beyond-line-of-sight (BLOS) CNPC link implemented through a satellite will almost always be required. Protected aviation spectrum (aeronautical mobile satellite (route) service, or AMS(R)S) would normally be used for such a safety-critical link, however studies have shown that currently available aviation safety satellite spectrum is inadequate to support the projected BLOS CNPC link bandwidth requirements. To address this inadequacy, the 2015 World Radio communication Conference studied the possible use of the Fixed Satellite Service (FSS) to provide CNPC, including possible allocations in Ku-Band and Ka-Band, under Agenda Item (AI) 1.5. Although UA CNPC satellite links in these bands were shown to meet operational availability and continuity requirements, a serious complication exists in that there are also terrestrial service allocations in these bands, in particular, Fixed Service (FS) point-to-point and point-to-multipoint microwave digital links. During the WRC-15 study cycle, much opposition to AI 1.5 was generated based on fears that UA CNPC satellite transmitters in these bands would impose unacceptable levels of interference to the FS receivers. NASA analyzed the possible interference from the UA transmitters based on probable UA transmission and FS receiver characteristics, and UA traffic distributions and densities to determine conditions under which UA could operate without imposing unacceptable interference levels to the FS. Ultimately, UA power flux density transmission limits were proposed as a way to insure protection of FS receivers and further studies were prepared on the various proposals. This paper presents the results of these studies and discusses possible implications on future UA BLOS operations.

aircraft communications↗

TURBOMAT: A Probabilistic Turbomachinery Aeroelastic Analysis Tool

An integration of aeroelastic analysis procedures with probabilistic analysis methods enables us to design safe reliable engines with quantified reliability. Towards this goal, a graphical user interface (GUI) based tool that integrates the codes Aeroelastic analysis of propfans (ASTROP2) and Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) is developed. The tool entitled TURBOMachinery Aeroelastic Analysis Tool (TURBOMAT), is developed utilizing the MATrix Laboratory (Matlab) Guide (Graphical User Interface Development) tool box. TURBOMAT provides a user friendly computational environment for rapid assessment of Turbomachinery blades flutter characteristics, subjected to uncertain loading conditions with variability in material and aerodynamic properties. The tool is seen as an education tool for new students and young engineers starting their careers in structural Aeroelasticity who want to learn and understand aeroelastic aspects of turbomachinery components, fans, compressors and turbines, including uncertainties in loading and material properties.A typical fan blade configuration geometry was chosen to demonstrate the tool. The results are presented in the form of probabilistic density function (PDF), the cumulative distribution function (CDF) and sensitivity factors. Both first order fast probability integration (FPI) and the Monto Carlo (MC) techniques are used in the analysis and compared. The tool enabled us to quantify blade flutter reliability as well as the ranking of uncertain variables and their importance to blade flutter response.

Aeroelastic Analysis↗

Collisional balance of the meteoritic complex

The present study has the objective to reevaluate the size distribution of interplanetary meteoroids on the basis of the most recent data, and to analyze the probable nature of the sinks and sources of meteoritic material. The flux of interplanetary meteorites at 1 AU is discussed, taking into account general characteristics, lunar crater distribution, flux curves, spatial densities, and cross-sectional distribution and light scattering. Collisional effects are examined, giving attention to catastrophic collisions, collision rate, and destroyed mass and generated fragments. The effect of radiation pressure on small particles is considered along with the difference between the lunar and interplanetary flux models, collisional evolution at 1 AU, potential sources for large meteoroids, and observational evidence of losses of small micrometeoroids.

Gruen, E.↗

Origin of the Valley Networks On Mars: A Hydrological Perspective

The geomorphology of the Martian valley networks is examined from a hydrological perspective for their compatibility with an origin by rainfall, globally higher heat flow, and localized hydrothermal systems. Comparison of morphology and spatial distribution of valleys on geologic surfaces with terrestrial fluvial valleys suggests that most Martian valleys are probably not indicative of a rainfall origin, nor are they indicative of formation by an early global uniformly higher heat flow. In general, valleys are not uniformly distributed within geologic surface units as are terrestrial fluvial valleys. Valleys tend to form either as isolated systems or in clusters on a geologic surface unit leaving large expanses of the unit virtually untouched by erosion. With the exception of fluvial valleys on some volcanoes, most Martian valleys exhibit a sapping morphology and do not appear to have formed along with those that exhibit a runoff morphology. In contrast, terrestrial sapping valleys form from and along with runoff valleys. The isolated or clustered distribution of valleys suggests localized water sources were important in drainage development. Persistent ground-water outflow driven by localized, but vigorous hydrothermal circulation associated with magmatism, volcanism, impacts, or tectonism is, however, consistent with valley morphology and distribution. Snowfall from sublimating ice-covered lakes or seas may have provided an atmospheric water source for the formation of some valleys in regions where the surface is easily eroded and where localized geothermal/hydrothermal activity is sufficient to melt accumulated snowpacks.

Gulick, Virginia C.↗

Relationship between radiation-induced aberrations in individual chromosomes and their DNA content: effects of interaction distance

PURPOSE: To study the effect of the interaction distance on the frequency of inter- and intrachromosome exchanges in individual chromosomes with respect to their DNA content. Assumptions: Chromosome exchanges are formed by misrejoining of two DNA double-strand breaks (DSB) induced within an interaction distance, d. It is assumed that chromosomes in G(0)/G(1) phase of the cell cycle occupy a spherical domain in a cell nucleus, with no spatial overlap between individual chromosome domains. RESULTS: Formulae are derived for the probability of formation of inter-, as well as intra-, chromosome exchanges relating to the DNA content of the chromosome for a given interaction distance. For interaction distances <1 microm, the relative frequency of interchromosome exchanges predicted by the present model is similar to that by Cigarran et al. (1998) based on the assumption that the probability of interchromosome exchanges is proportional to the "surface area" of the chromosome territory. The "surface area" assumption is shown to be a limiting case of d-->0 in the present model. The present model also predicts that the probability of intrachromosome exchanges occurring in individual chromosomes is proportional to their DNA content with correction terms. CONCLUSION: When the interaction distance is small, the "surface area" distribution for chromosome participation in interchromosome exchanges has been expected. However, the present model shows that for the interaction distance as large as 1 microm, the predicted probability of interchromosome exchange formation is still close to the surface area distribution. Therefore, this distribution does not necessarily rule out the formation of complex chromosomal aberrations by long-range misrejoining of DSB.

Non-programmatic↗

Probabilistic Assessment of Cancer Risk from Solar Particle Events

For long duration missions outside of the protection of the Earth s magnetic field, space radiation presents significant health risks including cancer mortality. Space radiation consists of solar particle events (SPEs), comprised largely of medium energy protons (less than several hundred MeV); and galactic cosmic ray (GCR), which include high energy protons and heavy ions. While the frequency distribution of SPEs depends strongly upon the phase within the solar activity cycle, the individual SPE occurrences themselves are random in nature. We estimated the probability of SPE occurrence using a non-homogeneous Poisson model to fit the historical database of proton measurements. Distributions of particle fluences of SPEs for a specified mission period were simulated ranging from its 5 th to 95th percentile to assess the cancer risk distribution. Spectral variability of SPEs was also examined, because the detailed energy spectra of protons are important especially at high energy levels for assessing the cancer risk associated with energetic particles for large events. We estimated the overall cumulative probability of GCR environment for a specified mission period using a solar modulation model for the temporal characterization of the GCR environment represented by the deceleration potential (^). Probabilistic assessment of cancer fatal risk was calculated for various periods of lunar and Mars missions. This probabilistic approach to risk assessment from space radiation is in support of mission design and operational planning for future manned space exploration missions. In future work, this probabilistic approach to the space radiation will be combined with a probabilistic approach to the radiobiological factors that contribute to the uncertainties in projecting cancer risks.

Kim, Myung-Hee Y.↗

2032 and 2036 risk enhancement from NEOs in the Taurid stream: Is there a significant coherent component to impact risk?

Impact risk is normally quantified by summing the product of the probability of an event and some measure of its consequences over the set of all possible events. The probability factor is considered to be more objective and is based on the size frequency distribution of Near Earth objects (NEOs) and an implicit assumption of randomness, which can be described as “stochastic risk”. Impact frequency does change with time, however, and there have been episodes in the deep geological past when the flux has been much higher. The hypothesis of “coherent catastrophism” suggests large variations on shorter timescales. It postulates the existence of a “Taurid resonant swarm” (TRS) of debris associated with Comet Encke that is stabilized by Jupiter and in a 7:2 resonance with it. The hypothetical cluster orbits in the broad Taurid stream, which crosses Earth's orbit twice a year at its nodes. Whereas the most extreme and fanciful versions of coherent catastrophism in recent geologic history have been comprehensively refuted, the possibility of a significant component of coherent (time-dependent) risk, associated with non-random correlated events, remains. This paper incorporates recently published data from observational campaigns associated with the Taurid stream. There is no evidence for objects in the Taurid stream that are above the global catastrophe threshold, but the possibility of a few large objects and a significant population of objects the size of the Chelyabinsk or Tunguska bodies has not been eliminated. Eyewitness accounts and comparison of airburst models to ground truth suggest that the Tunguska object was a Beta Taurid. The Tunguska event is probably an outlier when compared to the size frequency distributions under stochastic assumptions, but if there is a significant coherent component it may be representative of events that take place far more frequently, and risk assessments may have underestimated the contribution from airbursts. If so, the Earth will experience increased risk from objects this size, peaking in November 2032 and June 2036. We recommend targeted survey campaigns during these hypothetical TRS node crossings to quantify the population, search for potentially hazardous objects, and identify imminent impactors.

79 ASTRONOMY AND ASTROPHYSICS↗

Void statistics of the CfA redshift survey

Clustering properties of two samples from the CfA redshift survey, each containing about 2500 galaxies, are studied. A comparison of the velocity distributions via a K-S test reveals structure on scales comparable with the extent of the survey. The void probability function (VPF) is employed for these samples to examine the structure and to test for scaling relations in the galaxy distribution. The galaxy correlation function is calculated via moments of galaxy counts. The shape and amplitude of the correlation function roughly agree with previous determinations. The VPFs for distance-limited samples of the CfA survey do not match the scaling relation predicted by the hierarchical clustering models. On scales not greater than 10/h Mpc, the VPFs for these samples roughly follow the hierarchical pattern. A variant of the VPF which uses nearly all the data in magnitude-limited samples is introduced; it accounts for the variation of the sampling density with velocity in a magnitude-limited survey.

Vogeley, Michael S.↗

A Lumped Stochastic Model of Coupled Neutronic Assemblies

The objective of the present work is to apply the Harris ansatz to coupled assemblies in a lumped setting, i.e., ignoring the neutron phase dependence, which hitherto has not been considered. A Master equation for the single-time joint distribution of neutron number in each assembly, accounting for transfer of neutrons from one assembly to another via leakage transition probabilities and intercept or view factors is derived from which equations for the mean and variance in each assembly are developed. Because of the neutronic coupling, it is necessary to additionally compute the correlation function between pairs of assemblies which gives an expanded set of moment equations to solve. The reactivity of each assembly is allowed to vary with time, although a formal mechanism that is responsible for reactivity coupling between assemblies is not considered or proposed. Following Harris, it is then hypothesized that the number distribution in each assembly is a gamma distribution with assembly-specific mean and variance computed exactly from the moment equations. Illustrative numerical results are presented for the moments of a system of four coupled assemblies under various reactivity and neutronic coupling scenarios. The gamma distribution is then used to compute the probability that the neutron population exceeds a threshold value in any assembly. The report concludes with a discussion of recommended further work on this problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Logarithmic Bayesian Approach to Quantum Error Detection

We consider the problem of continuous quantum error correction from a Bayesian perspective, proposing a pair of digital filters using logarithmic probabilities that are able to achieve near-optimal performance on a three-qubit bit-flip code, while still being reasonable to implement on low-latency hardware. These practical filters are approximations of an optimal filter that we derive explicitly for finite time steps, in contrast with previous work that has relied on stochastic differential equations such as the Wonham filter. By utilizing logarithmic probabilities, we are able to eliminate the need for explicit normalization and can reduce the Gaussian noise distribution to a simple quadratic expression. The state transitions induced by the bit-flip errors are modeled using a Markov chain, which for log-probabilities must be evaluated using a LogSumExp function. We develop the two versions of our filter by constraining this LogSumExp to have either one or two inputs, which favors either simplicity or accuracy, respectively. Using simulated data, we demonstrate that the single-term and two-term filters are able to significantly outperform both a double threshold scheme and a linearized version of the Wonham filter in tests of error detection under a wide variety of error rates and time steps.

97 MATHEMATICS AND COMPUTING↗

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitary time, G(bar a, bar t), and that of the time when a crack reaches an arbitary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. TAking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitrary time, G(bar a, bar t), and that of the time when a crack reaches an arbitrary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. Taking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

Void probability as a function of the void's shape and scale-invariant models

The dependence of counts in cells on the shape of the cell for the large scale galaxy distribution is studied. A very concrete prediction can be done concerning the void distribution for scale invariant models. The prediction is tested on a sample of the CfA catalog, and good agreement is found. It is observed that the probability of a cell to be occupied is bigger for some elongated cells. A phenomenological scale invariant model for the observed distribution of the counts in cells, an extension of the negative binomial distribution, is presented in order to illustrate how this dependence can be quantitatively determined. An original, intuitive derivation of this model is presented.

Elizalde, E.↗

Void probability as a function of the void's shape and scale-invariant models

The dependence of counts in cells on the shape of the cell for the large scale galaxy distribution is studied. A very concrete prediction can be done concerning the void distribution for scale invariant models. The prediction is tested on a sample of the CfA catalog, and good agreement is found. It is observed that the probability of a cell to be occupied is bigger for some elongated cells. A phenomenological scale invariant model for the observed distribution of the counts in cells, an extension of the negative binomial distribution, is presented in order to illustrate how this dependence can be quantitatively determined. An original, intuitive derivation of this model is presented.

Elizalde, E.↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

Role of 'core' and 'halo' solar electrons in ionization of the interstellar medium

In the present consideration of the probability of interstellar wind H and He atoms' survival of ionization by solar wind electrons, a dual temperature electron distribution is for the first time used to model the effects of 'core' (10 eV) and 'halo' (60 eV) solar electrons on survival probabilities. The interstellar wind has a higher probability of surviving halo rather than core electron ionization only at heliocentric distances smaller than 0.5 AU. For larger distances, halo and core survival probabilities are similar, and similarly insensitive to changes in (1) interstellar wind velocity at infinity, (2) the solar ratio of radiation to gravitational force, and (3) solar electron temperatures. At distances smaller than 0.5 AU, however, sensitivity significantly increases.

Askew, S. D.↗