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 649 records · Page 36

A Plasma Drag Hypervelocity Particle Accelerator (HYPER)

Current debris models are able to predict the growth of the space debris problem and suggest that spacecraft must employ armor or bumper shields for some orbital altitudes now and that the problem will become worse as a function of time. The practical upper limit to the velocity distribution is on the order of 40 km/s and is associated with the natural environment. The velocity distribution of the man-made component peaks at 9-10 km/s with maximum velocity in the 14-16 km/s range. Experience in space has verified that the "high probability of impact" particles are in the microgram to milligram range. These particles can have very significant effects on coatings, insulators, and thin metallic layers. The surface of thick materials becomes pitted and the local debris component is enhanced by ejecta from the impact events. In this paper, the HYPER facility is described which produces a reasonable simulation of the man-made space debris spectrum in a controlled environment. The facility capability is discussed in terms of drive geometry, energetics, velocity distribution, diagnostics, and projectile/debris loading. The facility has been used to study impact phenomena on Space Station Freedom's solar array structure, the calibration of space debris collectors, other solar array materials, potential structural materials for use in space, electrical breakdown in the space environment, and as a means of clarifying or duplicating the impact phenomena on surfaces which have been exposed in space.

Best, Steve R.↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Observation of a driver gas-tangential discontinuity

A complete analysis of an interplanetary disturbance of Nov. 19, 1970 using the Apollo 12-SIDE (Suprathermal Ion Detector Experiment) is presented. The SIDE detectors were pointing at 26.3 degrees from the normal solar-wind direction during the observations. The data were least-squares fitted (using a parabolic hypersurface approximation) to a convected Maxwell-Boltzmann distribution function. The results of the fit combined with two other experiments showed a drastic change in the wind speed (from an steady 352 km/sec down to 219 km/sec), direction, and temperature. Except for a delta-function increase at the onset, the density remained constant. There was a considerable enhancement in the abundance of He and probably of heavier elements. The interplanetary magnetic field exhibited a jump of 21 gamma with a change in latitude from -56 to -76 degrees in solar ecliptic coordinates. It is concluded that the disturbance was due to the driver gas-tangential discontinuity of a solar flare-induced shock wave. The characteristic of the tangential discontinuity fit well with theoretical prediction.

Medrano, R. A.↗

Improved techniques for outgoing wave variational principle calculations of converged state-to-state transition probabilities for chemical reactions

Improved techniques and well-optimized basis sets are presented for application of the outgoing wave variational principle to calculate converged quantum mechanical reaction probabilities. They are illustrated with calculations for the reactions D + H2 yields HD + H with total angular momentum J = 3 and F + H2 yields HF + H with J = 0 and 3. The optimization involves the choice of distortion potential, the grid for calculating half-integrated Green's functions, the placement, width, and number of primitive distributed Gaussians, and the computationally most efficient partition between dynamically adapted and primitive basis functions. Benchmark calculations with 224-1064 channels are presented.

Mielke, Steven L.↗

Spatial analysis of cell patterning to aid genetic and phenotypic understanding of grass stomatal density: A case study in maize

Biological processes involve complex hierarchies where composite traits result from multiple component traits. However, holistically understanding of how sets of component traits interact to underpin genotype-to-phenotype relationships is generally lacking. Stomatal density (SD) is a tractable model system for exploring how high-throughput phenotyping (HTP) data could be exploited by a new spatial analysis approach to better understand a developmentally and functionally important trait. SD is a composite trait, resulting from various components related to cell identity and size, which are themselves governed by a series of spatio-developmental processes. Data from 192 recombinant inbred lines of maize [Zea mays (L.)] were analyzed by a new stomatal patterning phenotype (SPP) to (1) describe the average spatial probability distribution of the nearest neighboring stomata; (2) derive a core set of component traits related to cell size, cell packing, and positional probabilities; (3) build a structural equation model of component traits underlying SD; and (4) identify stomatal patterning quantitative trait loci (QTL). The core set of SPP-derived traits explained 74% of the variation in SD. Analyzing SPP component traits allowed some loci previously identified as generic SD QTL to be recognized as specific to lateral versus longitudinal elements of stomatal patterning. Therefore, this study highlights how novel insights can be gained by decomposing a composite trait (e.g., SD) into a set of component traits that were present in HTP data but not previously exploited.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗

Spacecraft dynamics characterization and control system failure detection, volume 1

The work under this grant has been directed to two aspects of the control of flexible spacecraft: (1) the modeling of deployed or erected structures including nonlinear joint characteristics; and (2) the detection and isolation of failures of the components of control systems for large space structures. The motivation for the first of these research tasks is the fact that very large assemblies in space will have to be built or deployed in situ. A likely scenario is, in fact, a combination of these wherein modules which are folded for transportation into orbit are erected to their final configuration and then jointed with other such erected modules to form the full assembly. Any such erectable modules will have joints. It remains to be seen whether or not joints designed for operational assemblies will have nonlinear properties, but it seems prudent to develop a methodology for dealing with that possibility. The motivation for the second of these research tasks is the fact that we foresee large assemblies in space which will require active control to damp vibrations and/or hold a desired shape. Lightweight structures will be very flexible, with many elastic modes having very low frequencies. In order to control these modes well, the control system will likely require many sensors and many actuators, probably distributed over much of the structure. The combination of a large number of control system components with long operational periods virtually guarantees that these systems will suffer control system component failures during operation. The control system must be designed to tolerate failures of some sensors and actuators, and still be able to continue to perform its function.

Source record↗

Evolved stars in Omega Centauri. I - Radial distribution of blue subdwarfs

A U - V color-magnitude diagram containing over 5000 stars from 10 fields in Omega Centauri at or brighter than the main-sequence turnoff is presented. Completeness corrections are obtained as a continuous function of magnitude and radial distance are obtained and used to show that the blue subdwarfs are centrally concentrated with respect to subgiants and stars on the horizontal branch proper. The chance probability of this result is less than 1 percent. Since the blue subdwarfs probably consist of helium-burning cores of about 0.5 solar mass surrounded by a thin hydrogen envelope, mass segregation could not produce this result if these stars had evolved singly. Two kinds of possible precursor systems are considered: pairs of degenerate dwarfs which subsequently merge to form helium-burning stars, and moderately wide binaries in which mass transfer is initiated shortly before the helium flash is ignited.

Bailyn, Charles D.↗

Interactive Soil Dust Aerosol Model in the GISS GCM: Sensitivity of the Soil Dust Cycle to Radiative Properties of Soil Dust Aerosols - Part 1

The sensitivity of the soil dust aerosol cycle to the radiative forcing by soil dust aerosols is studied. Four experiments with the NASA/GISS atmospheric general circulation model, which includes a soil dust aerosol model, are compared, all using a prescribed climatological sea surface temperature as lower boundary condition. In one experiment, dust is included as dynamic tracer only (without interacting with radiation), whereas dust interacts with radiation in the other simulations. Although the single scattering albedo of dust particles is prescribed to be globally uniform in the experiments with radiatively active dust, a different single scattering albedo is used in those experiments to estimate whether regional variations in dust optical properties, corresponding to variations in mineralogical composition among different source regions, are important for the soil dust cycle and the climate state. On a global scale, the radiative forcing by dust generally causes a reduction in the atmospheric dust load corresponding to a decreased dust source flux. That is, there is a negative feedback in the climate system due to the radiative effect of dust. The dust source flux and its changes were analyzed in more detail for the main dust source regions. This analysis shows that the reduction varies both with the season and with the single scattering albedo of the dust particles. By examining the correlation with the surface wind, it was found that the dust emission from the Saharan/Sahelian source region and from the Arabian peninsula, along with the sensitivity of the emission to the single scattering albedo of dust particles, are related to large scale circulation patterns, in particular to the trade winds during Northern Hemisphere winter and to the Indian monsoon circulation during summer. In the other regions, such relations to the large scale circulation were not found. There, the dependence of dust deflation to radiative forcing by dust particles is probably dominated by physical processes with short time scales. The experiments show that dust radiative forcing can lead to significant changes both in the soil dust cycle and in the climate state. To estimate dust concentration and radiative forcing by dust more accurately, dust size distributions and dust single scattering albedo in the model should be a function of the source region, because dust concentration and climate response to dust radiative forcing are sensitive to dust radiative parameters.

Perlwitz, Jan↗

The calculation of radial dose from heavy ions: predictions of biological action cross sections

The track structure model of heavy ion cross sections was developed by Katz and co-workers in the 1960s. In this model the action cross section is evaluated by mapping the dose-response of a detector to gamma rays (modeled from biological target theory) onto the radial dose distribution from delta rays about the path of the ion. This is taken to yield the radial distribution of probability for a "hit" (an interaction leading to an observable end-point). Radial integration of the probability yields the cross section. When different response from ions of different Z having the same stopping power is observed this model may be indicated. Since the 1960s there have been several developments in the computation of the radial dose distribution, in the measurement of these distributions, and in new radiobiological data against which to test the model. The earliest model, by Butts and Katz made use of simplified delta ray distribution functions, of simplified electron range-energy relations, and neglected angular distributions. Nevertheless it made possible the calculation of cross sections for the inactivation of enzymes and viruses, and allowed extension to tracks in nuclear emulsions and other detectors and to biological cells. It set the pattern for models of observable effects in the matter through which the ion passed. Here we outline subsequent calculations of radial dose which make use of improved knowledge of the electron emission spectrum, the electron range-energy relation, the angular distribution, and some considerations of molecular excitation, of particular interest both close to the path of the ion and the outer limits of electron penetration. These are applied to the modeling of action cross sections for the inactivation of several strains of E-coli and B. subtilis spores where extensive measurements in the "thin-down" region have been made with heavy ion beams. Such calculations serve to test the radial dose calculations at the outer limit of electron penetration. We lack data from which to test these calculations in regions close to the path of the ion aside from our earliest work on latent tracks in plastics, though it appears that the criterion then suggested for the threshold of track formation, of a minimal dose at a minimal distance (of about 20 angstroms, in plastics), remains valid.

NASA Discipline Radiation Health↗

Use of collateral information to improve LANDSAT classification accuracies

Methods to improve LANDSAT classification accuracies were investigated including: (1) the use of prior probabilities in maximum likelihood classification as a methodology to integrate discrete collateral data with continuously measured image density variables; (2) the use of the logit classifier as an alternative to multivariate normal classification that permits mixing both continuous and categorical variables in a single model and fits empirical distributions of observations more closely than the multivariate normal density function; and (3) the use of collateral data in a geographic information system as exercised to model a desired output information layer as a function of input layers of raster format collateral and image data base layers.

Strahler, A. H.↗

Soft X-ray observations of pre-main sequence stars in the chamaeleon dark cloud

Einstein IPC observations of the nearby Chamaeleon I star forming cloud show 22 well-resolved soft X-ray sources in a 1x2 deg region. Twelve are associated with H-alpha emission line pre-main sequence (PMS) stars, and four with optically selected PMS stars. Several X-ray sources have two or more PMS stars in their error circles. Optical spectra were obtained at CTIO of possible stellar counterparts of the remaining X-ray sources. They reveal 5 probable new cloud members, K7-MO stars with weak or absent emission lines. These naked X-ray selected PMS stars are similar to those found in the Taurus-Auriga cloud. The spatial distributions and H-R diagrams of the X-ray and optically selected PMS stars in the cloud are very similar. Luminosity functions indicate the Chamaeleon stars are on average approximately 5 times more X-ray luminous than Pleiad dwarfs. A significant correlation between L sub x and optical magnitude suggests this trend may continue within the PMS phase of stellar evolution. The relation of increasing X-ray luminosity with decreasing stellar ages is thus extended to stellar ages as young as 1 million years.

Feigelson, Eric D.↗

Determining large-scale heliospheric structure using ultraviolet resonance line observations

Currently the Pioneer 10 and Voyager 1 and 2 spacecraft are beyond the orbit of Pluto, traveling outward from the Sun. Each is capable of detecting ultraviolet radiation resonantly scattered from hydrogen and helium atoms in the heliosphere and local interstellar medium. These observations are particularly well suited for the investigation of the large-scale heliospheric H and He distributions because the Voyager spacecraft are heading upstream, into the direction of local interstellar flow, whereas Pioneer 10 is heading downstream. Observations of the brightest resonance line, H Lyman-alpha, reveals that beyond about 20 AU from the Sun, upstream intensities decrease less quickly as a function of solar distance than downstream intensities. This implies that the heliospheric H distributions in the upstream and downstream directions are significantly different. Heliospheric H atoms originate in the local interstellar flow, and must penetrate through the heliospheric interface, where they are subject to charge exchange collisions with solar wind and interstellar protons. Models indicate that this process is probably responsible for the upstream/downstream difference in H Lyman-alpha. In addition, a recent spectroscopic determination of the H atom velocity distribution in the inner heliosphere implies a significant deceleration in the bulk flow speed of the heliospheric hydrogen gas relative to the helium flow, an effect that is also likely due to H-p charge exchange occurring in the upstream heliospheric interface region. In this presentation, recent heliospheric resonance line observations and their interpretations will be reviewed, focusing on their sensitivity to large-scale heliospheric structure.

Hall, Doyle T.↗

The SMM Model as a Boundary Value Problem Using the Discrete Diffusion Equation

A generalized single step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

Campbell, Joel↗

Probability distribution of vertical longitudinal shear fluctuations.

This paper discusses some recent measurements of third and fourth moments of vertical differences (shears) of longitudinal velocity fluctuations obtained in unstable air at the NASA 150 m meteorological tower site at Cape Kennedy, Fla. Each set of measurements consisted of longitudinal velocity fluctuation time histories obtained at the 18, 30, 60, 90, 120 and 150 m levels, so that 15 wind-shear time histories were obtained from each set of measurements. It appears that the distribution function of the longitudinal wind fluctuations at two levels is not bivariate Gaussian. The implications of the results relative to the design and operation of aerospace vehicles are discussed.-

Fichtl, G. H.↗

Extreme Mean and Its Applications

Extreme value statistics obtained from normally distributed data are considered. An extreme mean is defined as the mean of p-th probability truncated normal distribution. An unbiased estimate of this extreme mean and its large sample distribution are derived. The distribution of this estimate even for very large samples is found to be nonnormal. Further, as the sample size increases, the variance of the unbiased estimate converges to the Cramer-Rao lower bound. The computer program used to obtain the density and distribution functions of the standardized unbiased estimate, and the confidence intervals of the extreme mean for any data are included for ready application. An example is included to demonstrate the usefulness of extreme mean application.

Swaroop, R.↗

Impact of stratospheric aircraft on calculations of nitric acid trihydrate cloud surface area densities using NMC temperatures and 2D model constituent distributions

A parameterization of NAT (nitric acid trihydrate) clouds is developed for use in 2D models of the stratosphere. The parameterization uses model distributions of HNO3 and H2O to determine critical temperatures for NAT formation as a function of latitude and pressure. National Meteorological Center temperature fields are then used to determine monthly temperature frequency distributions, also as a function of latitude and pressure. The fractions of these distributions which fall below the critical temperatures for NAT formation are then used to determine the NAT cloud surface area density for each location in the model grid. By specifying heterogeneous reaction rates as functions of the surface area density, it is then possible to assess the effects of the NAT clouds on model constituent distributions. We also consider the increase in the NAT cloud formation in the presence of a fleet of stratospheric aircraft. The stratospheric aircraft NO(x) and H2O perturbations result in increased HNO3 as well as H2O. This increases the probability of NAT formation substantially, especially if it is assumed that the aircraft perturbations are confined to a corridor region.

Considine, David B.↗

Decomposition of Electron Ionization Mass Spectra for Space Application Using a Monte-Carlo approach

Rationale: Quadrupole mass spectrometers equipped with an electron ionization (EI) sources have been widely used in space exploration to investigate the composition of planetary surfaces and atmospheres. However, the complexity of the samples and the minimal calibration for the fragmentation of molecules in the ionization chambers have prevented the deconvolution of the majority of the mass spectra obtained at different targets, thus limiting the determination of the exact composition of the samples analyzed. We propose a Monte‐Carlo approach to solve this issue mathematically. Methods: We decomposed simulated mass spectra of mixtures acquired with unit resolving power mass spectrometers and EI sources into the sum of the single components fragmentation patterns weighted by their relative concentration using interior‐point least‐square fitting. To fit compounds with poorly known fragmentation patterns, we used a Monte‐Carlo method to vary the intensity of individual fragment ions. We then decomposed the spectrum thousands of times to obtain a statistical distribution. Results: By performing the deconvolution on a mixture of seven different molecules with interfering fragmentation patterns (H2O, O2, CH4, Ar, N2, C2H4, and C2H6) we show that this approach retrieves the mixing ratio of the individual components more accurately than regular mass spectra decomposition methods that rely on fragmentation patterns from general databases. It also provides the probability density function for each species's mixing ratio. Conclusions: By removing the solution degeneracy in the decomposition of mass spectra, the method described herein could significantly increase the scientific retrieval from archived space flight mass spectrometry data, where calibration of the ionization source is no longer an option.

Thomas Gautier↗