Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability and statistical methods”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling Buried Object Brightness and Visibility for Ground Penetrating Radar

Comparing the observed brightness of various buried objects is a straightforward way to characterize the performance of a ground penetrating radar (GPR) system. However, a limitation arises. A simple comparison of buried object brightness values does not disentangle the effects of the GPR system itself from the system's operating environment and the objects being observed. Therefore, with brightness values exhibiting an unknown synthesis of systemic, environmental, and object factors, GPR system analysis becomes a convoluted affair. In this work, we use an experimentally collected dataset of over 25,000 object observations from five different multistatic radar arrays to develop models of buried object brightness and control for these various effects. Additionally, our modeling efforts provide a means for quantifying the relative brightness of GPR systems, the objects they detect, and the physical properties of those objects that influence observed brightness. To evaluate the models' performance on new object observations, we repeatedly simulate fitting them to half the dataset and predicting the observed brightness values of the unseen half. In addition, we introduce a method for estimating the probability that individual observations constitute a visible object, which aids in failure analysis, performance characterization, and dataset cleaning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗
DOE OSTI2024

CMB-PAInT: An inpainting tool for the cosmic microwave background

Abstract The presence of astrophysical emissions in microwave observations forces us to perform component separation to extract the Cosmic Microwave Background (CMB) signal. However, even in the most optimistic cases, there are still strongly contaminated regions, such as the Galactic plane or those with emission from extragalactic point sources, which require the use of a mask. Since many CMB analyses, especially the ones working in harmonic space, need the whole sky map, it is crucial to develop a reliable inpainting algorithm that replaces the values of the excluded pixels by others statistically compatible with the rest of the sky. This is especially important when working withQandUsky maps in order to obtainE- andB-mode maps which are free fromE-to-Bleakage. In this work we study a method based on Gaussian Constrained Realizations (GCR), that can deal with both intensity and polarization. Several tests have been performed to asses the validation of the method, including the study of the one-dimensional probability distribution function (1-PDF),E- andB-mode map reconstruction, and power spectra estimation. We have considered two scenarios for the input simulation: one case with only CMB signal and a second one including also Planck PR4 semi-realistic noise. Even if we are limited to low resolution maps, N side = 64 ifT,QandUare considered, we believe that this is a useful approach to be applied to future missions such as LiteBIRD, where the target are the largest scales.

Astronomy & Astrophysics↗
NASA NTRS1975

A method for the evaluation of damage on the surface of carbon-carbon ablation samples

The effect of weave geometry on the onset of transition from laminar to turbulent flow around a sample of an orthogonally woven fine weave carbon-carbon material tested in an arc-jet is investigated. Seventeen post-test models which had undergone only partial transition around the circumference were examined as to location of turbulent regions relative to the weave axes. It was statistically shown that the weave axes had the lowest probability of being a position of initial transition onset. The probability generally increased with distance from a major axis, with some rays, notably the 45 deg rays, showing better transition resistance than the adjacent rays. The incidence of initial transition onset correlated well with the number of X or Y fibers that the flow passes over from stagnation point to the shoulder.

Heinonen, E. W.↗
NASA NTRS1992

Probabilistic finite elements for fatigue and fracture analysis

Attenuation is focused on the development of Probabilistic Finite Element Method (PFEM), which combines the finite element method with statistics and reliability methods, and its application to linear, nonlinear structural mechanics problems and fracture mechanics problems. The computational tool based on the Stochastic Boundary Element Method is also given for the reliability analysis of a curvilinear fatigue crack growth. The existing PFEM's have been applied to solve for two types of problems: (1) determination of the response uncertainty in terms of the means, variance and correlation coefficients; and (2) determination the probability of failure associated with prescribed limit states.

Belytschko, Ted↗
NASA NTRS2023

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↗
NASA NTRS1978

Cumulative fade statistics and prediction methods associated with the COMSTAR beacon signal at 28.56 GHz

The paper presents cumulative fade statistics for the period April 1977 to March 1978 and reviews the rain attenuation predictive efforts using radar and disdrometer data for the summer of 1977. Attention is given to the receiving and radar systems as well as the disdrometer-raingage system. Discussion of the fade statistics covers the yearly, monthly and time of day statistics. Consideration is given to measured and predicted attenuation events, measured and predicted probability distributions as well as an empirical calibration adjustment. It is noted that in spite of many reasons for noncorrelation, the radar results do correlate well. Finally, it is concluded that the results demonstrate the utility of using radar coupled with disdrometer measurements for predicting individual fade events as well as long term fade distributions associated with satellite communications through rain.

Goldhirsh, J.↗
NASA NTRS1982

Optimal weighting of data to detect climatic change - Application to the carbon dioxide problem

It is suggested that a weighting of surface temperature data, using information about the expected level of warming in different seasons and geographical regions and statistical information about the amount of natural variability in surface temperature, can improve the chances of early detection of carbon dioxide concentration-induced climatic warming. A preliminary analysis of the optimal weighting method presented suggests that it is 25 per cent more effective in revealing surface warming than the conventional method, in virtue of the fact that 25 per cent more data must conventionally be analyzed in order to arrive at a similar probability of detection. An approximate calculation suggests that the warming ought to have already been detected, if the only sources of significant surface temperature variability had time scales of less than one year.

Bell, T. L.↗
NASA NTRS2001

Normal and abnormal tissue identification system and method for medical images such as digital mammograms

A system and method for analyzing a medical image to determine whether an abnormality is present, for example, in digital mammograms, includes the application of a wavelet expansion to a raw image to obtain subspace images of varying resolution. At least one subspace image is selected that has a resolution commensurate with a desired predetermined detection resolution range. A functional form of a probability distribution function is determined for each selected subspace image, and an optimal statistical normal image region test is determined for each selected subspace image. A threshold level for the probability distribution function is established from the optimal statistical normal image region test for each selected subspace image. A region size comprising at least one sector is defined, and an output image is created that includes a combination of all regions for each selected subspace image. Each region has a first value when the region intensity level is above the threshold and a second value when the region intensity level is below the threshold. This permits the localization of a potential abnormality within the image.

Heine, John J.↗
NASA NTRS1974

An empirical analysis of the distribution of the duration of overshoots in a stationary gaussian stochastic process

This analysis utilizes computer simulation and statistical estimation. Realizations of stationary gaussian stochastic processes with selected autocorrelation functions are computer simulated. Analysis of the simulated data revealed that the mean and the variance of a process were functionally dependent upon the autocorrelation parameter and crossing level. Using predicted values for the mean and standard deviation, by the method of moments, the distribution parameters was estimated. Thus, given the autocorrelation parameter, crossing level, mean, and standard deviation of a process, the probability of exceeding the crossing level for a particular length of time was calculated.

Parrish, R. S.↗
DOE OSTI2024

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

New interplanetary proton fluence model

A new predictive engineering model for the interplanetary fluence of protons with above 10 MeV and above 30 MeV is described. The data set used is a combination of observations made from the earth's surface and from above the atmosphere between 1956 and 1963 and observations made from spacecraft in the vicinity of earth between 1963 and 1985. The data cover a time period three times as long as the period used in earlier models. With the use of this data set the distinction between 'ordinary proton events' and 'anomalously large events' made in earlier work disappears. This permitted the use of statistical analysis methods developed for 'ordinary events' on the entire data set. The greater than 10 MeV fluences at 1 AU calculated with the new model are about twice those expected on the basis of models now in use. At energies above 30 MeV, the old and new models agree. In contrast to earlier models, the results do not depend critically on the fluence from any one event and are independent of sunspot number. Mission probability curves derived from the fluence distribution are presented.

Feynman, Joan↗
NASA NTRS2013

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An algorithm to locate the centers of baryon acoustic oscillations

The cosmic structure formed from baryon acoustic oscillations (BAO) in the early universe is imprinted in the galaxy distribution observable in large-scale surveys and is used as a standard ruler in contemporary cosmology. Typically, BAOs are detected as a preferential length scale in two-point statistics, which gives little information about the location of the BAO structures in real space. The aim of the algorithm described in this paper is to find probable centers of BAOs in the cosmic matter distribution. The algorithm convolves the three-dimensional distribution of matter density with a spherical shell kernel of variable radius placed at different locations. The locations that correspond to the highest values of the convolution correspond to the probable centers of BAOs. This method is realized in an open-source, computationally efficient algorithm. Here, we describe the algorithm and present the results of applying it to the SDSS DR9 CMASS survey and associated mock catalogs. A detailed performance study demonstrates the ability of the algorithm to locate BAO centers and in doing so presents a novel detection of the BAO scale in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Some questions of an optimum probability synthesis of dynamic metalworking machine systems

Attempts were made to develop a dynamic metal working machine tool system using probability information on external factors and independent disturbances of the parameters. The problem was optimized in two stages: (1) creation of mathematical models under conditions determined by input influences, and (2) processing the results of calculations obtained on such models. Data were obtained on the patterns of distribution of dynamic quality characteristics of the tool, vibration frequencies and amplitudes, attenuation coefficients, and duration of a transition process. These data were obtained using a digital computer, statistical test methods, or logical possibility trees.

Dobrynin, S. A.↗
NASA NTRS2011

Method for Identifying Probable Archaeological Sites from Remotely Sensed Data

Archaeological sites are being compromised or destroyed at a catastrophic rate in most regions of the world. The best solution to this problem is for archaeologists to find and study these sites before they are compromised or destroyed. One way to facilitate the necessary rapid, wide area surveys needed to find these archaeological sites is through the generation of maps of probable archaeological sites from remotely sensed data. We describe an approach for identifying probable locations of archaeological sites over a wide area based on detecting subtle anomalies in vegetative cover through a statistically based analysis of remotely sensed data from multiple sources. We further developed this approach under a recent NASA ROSES Space Archaeology Program project. Under this project we refined and elaborated this statistical analysis to compensate for potential slight miss-registrations between the remote sensing data sources and the archaeological site location data. We also explored data quantization approaches (required by the statistical analysis approach), and we identified a superior data quantization approached based on a unique image segmentation approach. In our presentation we will summarize our refined approach and demonstrate the effectiveness of the overall approach with test data from Santa Catalina Island off the southern California coast. Finally, we discuss our future plans for further improving our approach.

Tilton, James C.↗