Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability Distribution Function”

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 595 records · Page 33

Microstructural characterization of random packings of cubic particles

Understanding the properties of random packings of solid objects is of critical importance to a wide variety of fundamental scientific and practical problems. The great majority of the previous works focused, however, on packings of spherical and sphere-like particles. We report the first detailed simulation and characterization of packings of non-overlapping cubic particles. Such packings arise in a variety of problems, ranging from biological materials, to colloids and fabrication of porous scaffolds using salt powders. In addition, packing of cubic salt crystals arise in various problems involving preservation of pavements, paintings, and historical monuments, mineral-fluid interactions, CO 2 sequestration in rock, and intrusion of groundwater aquifers by saline water. Not much is known, however, about the structure and statistical descriptors of such packings. We have developed a version of the random sequential addition algorithm to generate such packings, and have computed a variety of microstructural descriptors, including the radial distribution function, two-point probability function, orientational correlation function, specific surface, and mean chord length, and have studied the effect of finite system size and porosity on such characteristics. Here, the results indicate the existence of both spatial and orientational long-range order in the packing, which is more distinctive for higher packing densities.

36 MATERIALS SCIENCE↗

Extreme_SeaState_Contour_v1

This software generates environmental contours of extreme sea states using buoy observations of significant wave height and energy period or peak period. The code transforms these observations using principal component analysis (PCA) to create an uncorrelated representation of the data. The subsequent components are modeled using probability distributions and parameter fitting functions. The inverse first-order reliability method (I-FORM) is then applied to these models in order to generate an extreme event contour based on a given return period (i.e., 100 years).The subsequent contour is then transformed back into the original input space defined by the variables of interest in order to create an environmental contour of extreme sea states.

Eckert, Aubrey↗

A New Probability Table Methodology with Covariances on Probability Tables

Through this work, we test a new methodology of evaluating probability tables informed by a combination of physics models and experimental data. This methodology uses generalized least squares to produce a joint mean and covariance for average cross sections and probability tables. This methodology has been tested and validated using total Monte Carlo. Attempts were made to parameterize the probability table using a formulation provided in. Although the normal inverse Gaussian fit was successful for the total cross section probability distribution, the rational-function fit to the capture, inelastic, and fission cross sections was unsuccessful due to time constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An alternative to Monte Carlo.

Histogram construction via algebraic Monte Carlo method, determining simultaneous effect of several input statistical variables on output variable

PROBABILITY THEORY↗

On the probability of cycle-slipping in first-order phase-locked loops.

The first-passage time boundary value problem for first-order phase-locked loops (PLL) is analyzed, and spectral representations are developed for the probability density function (pdf), the distribution function, and the moments of the first time to passage (or cycle-slip). For the sinusoidal PLL, an asymptotic formula, that is surprisingly accurate even at low loop SNR's and large frequency offsets, is obtained for the pdf of the time to cycle-slip, in terms of the mean time to slip.

La Frieda, J. R.↗

Turbulence control on an airborne laser platform

An active flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time and then release all of the flow instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. A flat plate towed in a water channel is used as a test bed. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, it is shown that the cyclic jet will produce periodic structures that are similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall of the plate. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.

Gad-El-hak, Mohamed↗

Turbulence control on airborne laser platform

An avctive flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A flat plate towed in a water channel is used as a test bed. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time during its 'on' period. When the jet is turned 'off', all of the turbulent fluid is released instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, the cyclic jet produces periodic structures that are qualitatively similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.

Gad-El-hak, Mohamed↗

Probabilistic material degradation model for aerospace materials subjected to high temperature, mechanical and thermal fatigue, and creep

A probabilistic general material strength degradation model has been developed for structural components of aerospace propulsion systems subjected to diverse random effects. The model has been implemented in two FORTRAN programs, PROMISS (Probabilistic Material Strength Simulator) and PROMISC (Probabilistic Material Strength Calibrator). PROMISS calculates the random lifetime strength of an aerospace propulsion component due to as many as eighteen diverse random effects. Results are presented in the form of probability density functions and cumulative distribution functions of lifetime strength. PROMISC calibrates the model by calculating the values of empirical material constants.

Boyce, L.↗

Potential Sources of Polarized Light from a Plant Canopy

Field measurements have demonstrated that sunlight polarized during a first surface reflection by shiny leaves dominates the optical polarization of the light reflected by shiny-leafed plant canopies having approximately spherical leaf angle probability density functions ("Leaf Angle Distributions" - LAD). Yet for other canopies - specifically those without shiny leaves and/or spherical LADs - potential sources of optically polarized light may not always be obvious. Here we identify possible sources of polarized light within those other canopies and speculate on the ecologically important information polarization measurements of those sources might contain.

Polarized Light↗

Application of the Global Precipitation Forecast dataset for Global Landslide Forecasting System

Landslides are extremely damaging, pervasive and cause fatalities and economic impacts. Extreme rainfall events, coupled with inopportune surface conditions, are the primary triggers of landslides around the world. Forecasting landslide events represent an area of open research. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model has been created that provides near real-time dynamic landslide characterization using Integrated Multi-Satellite Retrievals for Global Precipitation Mission (IMERG). However, it does not provide information on prediction of landslides into future. This study considers how global precipitation forecast data compares to satellite rainfall at different spatiotemporal scales and outlines the potential for its use in landslide hazard prediction/forecasting system. NASA Goddard Earth Observing System (GEOS)-Forecast model assimilates new observations every 6 hours, at 00, 06, 12, and 18 UTC. The GEOS-Forecast model is initialized at 00 UTC and is evaluated with IMERG Early, using Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) precipitation product as a reference over contiguous United States (CONUS). Categorical and continuous statistics along with probability density functions and cumulative distribution functions are considered to assess the performance of the precipitation products, with a focus on landslide hotspots. Seasonality appears to influence the performance of both the GEOS-Forecast and IMERG Early product. Moreover, the global comparison between the GEOS-Forecast and IMERG Early is carried out in terms of percentile difference, correlation, and bias maps. For extreme rainfall events in regions such as Mekong, Colombia, and Tajikistan, GEOS-Forecast appears to resolve high rainfall relative to IMERG Early more frequently. Validation over landslide points reveal that for 24hr rainfall accumulations > 100mm, GEOS-Forecast and IMERG coincide. GEOS-Forecast and IMERG Early precipitation matches more closely for tropical cyclones than other types of storms. Overall, the performance varies with respect to location, rainfall intensities, and type of precipitation events.

Sana Khan↗

Extreme_SeaState_Contour_v1

This software generates environmental contours of extreme sea states using buoy observations of significant wave height and energy period or peak period. The code transforms these observations using principal component analysis (PCA) to create an uncorrelated representation of the data. The subsequent components are modeled using probability distributions and parameter fitting functions. The inverse first-order reliability method (I-FORM) is then applied to these models in order to generate an extreme event contour based on a given return period (i.e., 100 years).The subsequent contour is then transformed back into the original input space defined by the variables of interest in order to create an environmental contour of extreme sea states.

Source record↗

The velocity distribution function of galaxy clusters as a cosmological probe

We present a new approach for quantifying the abundance of galaxy clusters and constraining cosmological parameters using dynamical measurements. In the standard method, galaxy line-of-sight velocities, v, or velocity dispersions are used to infer cluster masses, M, to quantify the halo mass function (HMF), dn(M)/dlog(M), which is strongly affected by mass measurement errors. In our new method, the probability distributions of velocities for each cluster in the sample are summed to create a new statistic called the velocity distribution function (VDF), dn(v)/dv. The VDF can be measured more directly and precisely than the HMF and can be robustly predicted with cosmological simulations that capture the dynamics of subhalos or galaxies. We apply these two methods to realistic (ideal) mock cluster catalogs with (without) interlopers and forecast the bias and constraints on the matter density parameter Ω{sub m} and the amplitude of matter fluctuations σ{sub 8} in flat ΛCDM cosmologies. For an example observation of 200 massive clusters, the VDF with (without) interloping galaxies constrains the parameter combination σ{sub 8} Ω{sub m}{sup 0.29(0.29)}=0.589±0.014 (0.584±0.011) and shows only minor bias. However, the HMF with interlopers is biased to low Ω{sub m} and high σ{sub 8} and the fiducial model lies well outside of the forecast constraints, prior to accounting for Eddington bias. When the VDF is combined with constraints from the cosmic microwave background, the degeneracy between cosmological parameters can be significantly reduced. Upcoming spectroscopic surveys that probe larger volumes and fainter magnitudes will provide clusters for applying the VDF as a cosmological probe.

79 ASTRONOMY AND ASTROPHYSICS↗

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

A family of random walks with generalized Dirichlet steps

We analyze a class of continuous time random walks in R{sup d},d≥2, with uniformly distributed directions. The steps performed by these processes are distributed according to a generalized Dirichlet law. Given the number of changes of orientation, we provide the analytic form of the probability density function of the position (X{sub {sub d}}(t),t>0) reached, at time t > 0, by the random motion. In particular, we analyze the case of random walks with two steps. In general, it is a hard task to obtain the explicit probability distributions for the process (X{sub {sub d}}(t),t>0). Nevertheless, for suitable values for the basic parameters of the generalized Dirichlet probability distribution, we are able to derive the explicit conditional density functions of (X{sub {sub d}}(t),t>0). Furthermore, in some cases, by exploiting the fractional Poisson process, the unconditional probability distributions of the random walk are obtained. This paper extends in a more general setting, the random walks with Dirichlet displacements introduced in some previous papers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Breast dose in mammography is about 30% lower when realistic heterogeneous glandular distributions are considered

Purpose: Current dosimetry methods in mammography assume that the breast is comprised of a homogeneous mixture of glandular and adipose tissues. Three-dimensional (3D) dedicated breast CT (bCT) data sets were used previously to assess the complex anatomical structure within the breast, characterizing the statistical distribution of glandular tissue in the breast. The purpose of this work was to investigate the effect of bCT-derived heterogeneous glandular distributions on dosimetry in mammography. Methods: bCT-derived breast diameters, volumes, and 3D fibroglandular distributions were used to design realistic compressed breast models comprised of heterogeneous distributions of glandular tissue. The bCT-derived glandular distributions were fit to biGaussian functions and used as probability density maps to assign the density distributions within compressed breast models. The MCNPX 2.6.0 Monte Carlo code was used to estimate monoenergetic normalized mean glandular dose “DgN(E)” values in mammography geometry. The DgN(E) values were then weighted by typical mammography x-ray spectra to determine polyenergetic DgN (pDgN) coefficients for heterogeneous (pDgN{sub hetero}) and homogeneous (pDgN{sub homo}) cases. The dependence of estimated pDgN values on phantom size, volumetric glandular fraction (VGF), x-ray technique factors, and location of the heterogeneous glandular distributions was investigated. Results: The pDgN{sub hetero} coefficients were on average 35.3% (SD, 4.1) and 24.2% (SD, 3.0) lower than the pDgN{sub homo} coefficients for the Mo–Mo and W–Rh x-ray spectra, respectively, across all phantom sizes and VGFs when the glandular distributions were centered within the breast phantom in the coronal plane. At constant breast size, increasing VGF from 7.3% to 19.1% lead to a reduction in pDgN{sub hetero} relative to pDgN{sub homo} of 23.6%–27.4% for a W–Rh spectrum. Displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width in the superior and inferior directions, resulted in a 37.3% and a −26.6% change in the pDgN{sub hetero} coefficient, respectively, relative to the centered distribution for the Mo–Mo spectrum. Lateral displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width, resulted in a 1.5% change in the pDgN{sub hetero} coefficient relative to the centered distribution for the W–Rh spectrum. Conclusions: Introducing bCT-derived heterogeneous glandular distributions into mammography phantom design resulted in decreased glandular dose relative to the widely used homogeneous assumption. A homogeneous distribution overestimates the amount of glandular tissue near the entrant surface of the breast, where dose deposition is exponentially higher. While these findings are based on clinically measured distributions of glandular tissue using a large cohort of women, future work is required to improve the classification of glandular distributions based on breast size and overall glandular fraction.

60 APPLIED LIFE SCIENCES↗

Identification of topology changes in power grids using phasor measurements

Phasor measurement units (PMUs) are increasingly important for monitoring the state of an electrical power grid and quickly detecting topology changes caused by events such as lines going down or large loads being dropped. Phasors are complex‐valued measurements of voltage and current at various points of generation and consumption. If a line goes down or a load is removed, power flows change throughout the grid according to known physical laws, and the probability distribution of phasor measurements changes accordingly. This paper develops a method to estimate the current topology of a power grid from phasor measurements and considers the design goal of placing PMUs at strategic points in a distribution system to achieve good sensitivity to single‐line outages. From a vector of phasor measurements, probabilities are computed corresponding to the scenario that all power lines are operational and to alternate scenarios in which each line goes down individually. These probabilities are functions of the joint distributions of phasor measurements under each possible scenario, obtained through Monte Carlo simulations with random load profiles. We use log‐spline densities to estimate marginal distributions of phasor measurements and fold these into a multivariate Gaussian copula to capture important correlations. Sensitivity to outages varies according to which line goes down and where PMUs are placed on the grid. A greedy search algorithm is demonstrated for placing PMUs at locations that provide good sensitivity to single‐line outages. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.

Wiel, Scott Vander↗