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 703 records · Page 39

Stochastic estimation of organized structures in turbulent channel flow

Linear stochastic estimation is used to approximate the conditional vector fields associated with high Reynolds stress producing events in numerically simulated turbulent channel flow. Joint probability density distributions of u-v are presented for y(+) between 1.35 and 180, and at each position the u-v values that make the greatest contributions to the average uv3 are used to define conditional ejection- and sweep-type events. In the y-z plane the conditional fields appear to be pairs of counter-rotating vortices. When the conditional event defined by the weighted probability density function is specified at a point close to the wall, the spacing and size of these eddies are consistent with the accepted wall streak spacing. As this point moves outwards, the eddy size increases significantly. An abrupt change in the flow angle occurs in the buffer layer and may indicate transition from streamwise oriented wall layer structures to hairpin vortices characterizing the outer layer.

Moin, P.↗

Selecting Critical Scenarios of DER Adoption in Distribution Grids Using Bayesian Optimization

We develop a new methodology to select scenarios of DER adoption most critical for distribution grids. Anticipating risks of future voltage and line flow violations due to additional PV adopters is central for utility investment planning but continues to rely on deterministic or ad hoc scenario selection. We propose a highly efficient search framework based on multi-objective Bayesian Optimization. We treat underlying grid stress metrics as computationally expensive black-box functions, approximated via Gaussian Process surrogates and design an acquisition function based on probability of scenarios being Pareto-critical across a collection of line- and bus-based violation objectives. Our approach provides a statistical guarantee and offers an order of magnitude speed-up relative to a conservative exhaustive search. Case studies on realistic feeders with 200-400 buses demonstrate the effectiveness and accuracy of our approach.

Mulkin, Olivier↗

A statistical model for radar images of agricultural scenes

The presently derived and validated statistical model for radar images containing many different homogeneous fields predicts the probability density functions of radar images of entire agricultural scenes, thereby allowing histograms of large scenes composed of a variety of crops to be described. Seasat-A SAR images of agricultural scenes are accurately predicted by the model on the basis of three assumptions: each field has the same SNR, all target classes cover approximately the same area, and the true reflectivity characterizing each individual target class is a uniformly distributed random variable. The model is expected to be useful in the design of data processing algorithms and for scene analysis using radar images.

Frost, V. S.↗

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗

What is a Question?

A given question can be defined in terms of the set of statements or assertions that answer it. Application of the logic of inference to this set of assertions allows one to derive the logic of inquiry among questions. There are interesting symmetries between the logics of inference and inquiry; where probability describes the degree to which a premise implies an assertion, there exists an analogous quantity that describes the bearing or relevance that a question has on an outstanding issue. These have been extended to suggest that the logic of inquiry results in functional relationships analogous to, although more general than, those found in information theory. Employing lattice theory, I examine in greater detail the structure of the space of assertions and questions demonstrating that the symmetries between the logical relations in each of the spaces derive directly from the lattice structure. Furthermore, I show that while symmetries between the spaces exist, the two lattices are not isomorphic. The lattice of assertions is described by a Boolean lattice 2(sup N) whereas the lattice of real questions is shown to be a sublattice of the free distributive lattice FD(N) = 2(sup 2(sup N)). Thus there does not exist a one-to-one mapping of assertions to questions, there is no reflection symmetry between the two spaces, and questions in general do not possess unique complements. Last, with these lattice structures in mind, I discuss the relationship between probability, relevance and entropy.

Knuth, Kevin H.↗

A statistical model for interpreting computerized dynamic posturography data

Computerized dynamic posturography (CDP) is widely used for assessment of altered balance control. CDP trials are quantified using the equilibrium score (ES), which ranges from zero to 100, as a decreasing function of peak sway angle. The problem of how best to model and analyze ESs from a controlled study is considered. The ES often exhibits a skewed distribution in repeated trials, which can lead to incorrect inference when applying standard regression or analysis of variance models. Furthermore, CDP trials are terminated when a patient loses balance. In these situations, the ES is not observable, but is assigned the lowest possible score--zero. As a result, the response variable has a mixed discrete-continuous distribution, further compromising inference obtained by standard statistical methods. Here, we develop alternative methodology for analyzing ESs under a stochastic model extending the ES to a continuous latent random variable that always exists, but is unobserved in the event of a fall. Loss of balance occurs conditionally, with probability depending on the realized latent ES. After fitting the model by a form of quasi-maximum-likelihood, one may perform statistical inference to assess the effects of explanatory variables. An example is provided, using data from the NIH/NIA Baltimore Longitudinal Study on Aging.

NASA Discipline Neuroscience↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Computer simulation in template-directed oligonucleotide synthesis

It is commonly assumed that template-directed polymerizations have played a key role in prebiotic evolution. A computer simulation that models up to 33 competing reactions was used to investigate the product distribution in a template-directed oligonucleotide synthesis as a function of time and concentration of the reactants. The study focuses on the poly(C)-directed elongation reaction of oligoguanylates, and how it is affected by the competing processes of hydrolysis and dimerization of the activated monomer, which have the potential of severely curtailing the elongation and reducing the size and yield of the synthesized polymers. The simulations show that realistic and probably prebiotically plausible conditions can be found where hydrolysis and dimerization are either negligible or where a high degree of polymerization can be attained even in the face of substantial hydrolysis and/or dimerization.

Kanavarioti, Anastassia↗

Achieving unequal error protection with convolutional codes

This paper examines the unequal error protection capabilities of convolutional codes. Both time-invariant and periodically time-varying convolutional encoders are examined. The effective free distance vector is defined and is shown to be useful in determining the unequal error protection (UEP) capabilities of convolutional codes. A modified transfer function is used to determine an upper bound on the bit error probabilities for individual input bit positions in a convolutional encoder. The bound is heavily dependent on the individual effective free distance of the input bit position. A bound relating two individual effective free distances is presented. The bound is a useful tool in determining the maximum possible disparity in individual effective free distances of encoders of specified rate and memory distribution. The unequal error protection capabilities of convolutional encoders of several rates and memory distributions are determined and discussed.

Mills, D. G.↗

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗

A simulation framework for evaluating electronic order workflows in integrated health records

Electronic health record (EHR) systems are critical to modern healthcare delivery, yet the dynamic workflows that govern electronic order processing remain underexplored. Inefficiencies in these digital pathways can cause delays in care, repetitive workloads, and even patient harm. This study presents a discrete-event simulation framework used to reconstruct and evaluate EHR-based order workflows in a large integrated healthcare system. Using real-world data extracted from the Veterans Health Administration’s Corporate Data Warehouse, the authors mapped order events to standardized state transitions and modeled their progression across different facilities of varying complexity levels. After being calibrated with empirical distributions of transition times and validated against observed time-in-system metrics, the simulation demonstrates close alignment with historical performance. Scenario analyses reveal that resource capacity constraints significantly amplify the impact of electronic order surges, which are reflected in the disproportionate growth in backlogs and processing delays. Adjustments in transition probabilities further increased recirculation and extended workflow paths. Network-based analysis identified Reserved, InProgress, and Completed as structurally critical states that function as hubs within the process network but the transitions in-between also act as major bottlenecks. These results showcased the effectiveness of simulation-based approaches in monitoring EHR order processing performance and evaluating consequences of workflow changes on healthcare network resources planning. The proposed simulation framework provides a scalable data-driven tool to support operational decision-making and improve the efficiency of electronic order management in complex healthcare environments.

Engineering↗

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗

Split-Step Simulations to Assess the Effects of Atmospheric Boundary Layer Turbulence on the Dose Variability of N-Waves and Shaped Booms

The effects of atmospheric boundary layer turbulence on the loudness variability of a sonic boom N-wave and shaped boom are examined with split-step simulations. The shaped boom is representative of a design iteration of the NASA X-59 aircraft. Inhomogeneous atmospheric boundary layer turbulence is generated in the computational domain by a Fourier synthesis method. The N-wave and shaped boom are propagated through turbulent fields representing eight different convection levels measured at the NASA Kennedy Space Center and the NASA Armstrong Flight Research Center. Probability density functions of the formation of caustic regions along the propagation direction are computed from the N-wave results, and a parameter to collapse the caustic PDFs that accounts for both fluctuation intensities and length scales is proposed. Statistical results concerning loudness metric variability are presented, and the standard deviations of several metrics are shown to collapse across different convection levels of turbulence for small nondimensional propagation distances. The loudness metric distributions are observed to be well approximated by a normal distribution for a given range of propagation distances, and become increasingly skewed as distance increases. A model function for the dose variability is proposed, and the function parameters are found to be related to the convection level of the turbulence. The model for the dose variability distribution is compared to simulation data that were not used to find the regression parameters of the model. At several nondimensional propagation distances, agreement is observed between the model and the simulation data. These results indicate that the model may be suitable for providing quick estimates of noise dose variability in the primary carpet region across a wide range of atmospheric boundary layer conditions.

Sonic boom↗

Evaluating design safety margins in the American Society of Mechanical Engineers graphite core components design-by-analysis assessments

Graphite is an important material being used for core components in next-generation high-temperature gas-cooled nuclear reactors. The selection of graphite grade for a specific Designer is a complex task, dependent on reactor conditions, component functionality, and required reliability. The American Society of Mechanical Engineers (ASME) provides two semi-probabilistic design-by-analysis assessments to evaluate graphite core components against design reliability targets. The simplified assessment uses a 2-parameter Weibull distribution to describe the graphite grade’s tensile-strength distribution to establish component stress limits. The full assessment uses the 3-parameter Weibull distribution and a modified Weakest-Link Theory approach to calculate a component design probability of failure. The paper defines recommended assessment rules, which are the as-written simplified assessment and the full assessment with parameter lower bounds, the modulus update with threshold reduction, and the 2027 grouping rules. Code rules are applied to three grades: 2114, IG-110, and NBG-18. The baseline margin calculation is developed using the experimental tensile dogbone specimen. Percent margin is defined as the percent reduction in the median experimental load to obtain the allowable load per ASME assessments. Under the recommended rules, the SRC–1 margin in the simplified assessment ranged from 40.2 % to 52.7 % among the grades in this study and from 36.1 % to 49.8 % in the full assessment. The full assessment only decreases the margins by 2.5–4.5 % for the SRC-1 components and 0–1.5 % for the SRC-2 components for this baseline case. Margin is inversely related to material median strength (i.e., the strongest grade, 2114, has the lowest margin).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A mathematical examination of the press model for atmospheric turbulence

The random process used to model atmospheric turbulence in aircraft response problems is examined. The first, second, and higher order probability density and characteristic functions were developed. The concepts of the Press model lead to an approximate procedure for the analysis of the response of linear dynamic systems to a class of non-Gaussian random processes. The Press model accounts for both the Gaussian and non-Gaussian forms of measured turbulence data. The nonstationary aspects of measured data are explicitly described by the transition properties of the random process. The effects of the distribution of the intensity process upon calculated exceedances are examined. It is concluded that the press model with a Gaussian intensity distribution gives a conservative prediction of limit load values.

Sidwell, K.↗

Sharp edges to neutral hydrogen disks in galaxies and the extragalactic radiation field

It is shown that the very sharp truncation of the neutral hydrogen distribution seen in NGC 3198 (and probably M33) is well modeled as the result of ionization of the atomic gas by the extragalactic radiation field. Below a critical column density of about a few times 10 exp 19/sq cm the gas is dominantly ionized and undetectable in the 21-cm line. It is inferred from the photoionization models that the total disk gas distribution in NGC 3198 is actually fairly axisymmetric. The critical column density for ionization is not a strong function of galaxy mass or mass distribution; thus, all galaxies should show a cutoff at approximately the same column density. Specific models of 3198 suggest that the extragalactic ionizing photon flux is 5000-10,000 photons/sq cm s.

Maloney, Philip↗

Probabilistic Component Mode Synthesis of Nondeterministic Substructures

Standard methods of structural dynamic analysis assume that the structural characteristics are deterministic. Recognizing that these characteristics are actually statistical in nature researchers have recently developed a variety of methods that use this information to determine probabilities of a desired response characteristic, such as natural frequency, without using expensive Monte Carlo simulations. One of the problems in these methods is correctly identifying the statistical properties of primitive variables such as geometry, stiffness, and mass. We present a method where the measured dynamic properties of substructures are used instead as the random variables. The residual flexibility method of component mode synthesis is combined with the probabilistic methods to determine the cumulative distribution function of the system eigenvalues. A simple cantilever beam test problem is presented that illustrates the theory.

Brown, Andrew M.↗

Regions of constrained maximum likelihood parameter identifiability

This paper considers the parameter identification problem of general discrete-time, nonlinear, multiple-input/multiple-output dynamic systems with Gaussian-white distributed measurement errors. Knowledge of the system parameterization is assumed to be known. Regions of constrained maximum likelihood (CML) parameter identifiability are established. A computation procedure employing interval arithmetic is proposed for finding explicit regions of parameter identifiability for the case of linear systems. It is shown that if the vector of true parameters is locally CML identifiable, then with probability one, the vector of true parameters is a unique maximal point of the maximum likelihood function in the region of parameter identifiability and the CML estimation sequence will converge to the true parameters.

Lee, C.-H.↗