Engineering PapersSearch

SEARCH · Engineering Papers

Results for “GAUSSIAN DISTRIBUTION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

One dimensional representations in quantum optics

The possibility of representing the quantum states of a harmonic oscillator not on the whole alpha-plane but on its one dimensional manifolds is considered. It is shown that a simple Gaussian distribution along a straight line describes a quadrature squeezed state while a similar Gaussian distribution along a circle leads to the amplitude squeezed state. The connection between the one dimensional representations and the usual Glauber representation is discussed.

Janszky, J.

Plasma sheet turbulence observed by Cluster II

Cluster fluxgate magnetometer (FGM) and ion spectrometer (CIS) data are employed to analyze magnetic field fluctuations within the plasma sheet during passages through the magnetotail region in the summers of 2001 and 2002 and, in particular, to look for characteristics of magnetohydrodynamic (MHD) turbulence. Power spectral indices determined from power spectral density functions are on average larger than Kolmogorov's theoretical value for fluid turbulence as well as Kraichnan's theoretical value for MHD plasma turbulence. Probability distribution functions of the magnetic fluctuations show a scaling law over a large range of temporal scales with non-Gaussian distributions at small dissipative scales and inertial scales and more Gaussian distribution at large driving scales. Furthermore, a multifractal analysis of the magnetic field components shows scaling behavior in the inertial range of the fluctuations from about 20 s to 13 min for moments through the fifth order. Both the scaling behavior of the probability distribution functions and the multifractal structure function suggest that intermittent turbulence is present within the plasma sheet. The unique multispacecraft aspect and fortuitous spacecraft spacing allow us to examine the turbulent eddy scale sizes. Dynamic autocorrelation and cross correlation analysis of the magnetic field components allow us to determine that eddy scale sizes fit within the plasma sheet. These results suggest that magnetic field turbulence is occurring within the plasma sheet resulting in turbulent energy dissipation.

Weygand, James M.

Chance-Constrained Guidance With Non-Convex Constraints

Missions to small bodies, such as comets or asteroids, require autonomous guidance for descent to these small bodies. Such guidance is made challenging by uncertainty in the position and velocity of the spacecraft, as well as the uncertainty in the gravitational field around the small body. In addition, the requirement to avoid collision with the asteroid represents a non-convex constraint that means finding the optimal guidance trajectory, in general, is intractable. In this innovation, a new approach is proposed for chance-constrained optimal guidance with non-convex constraints. Chance-constrained guidance takes into account uncertainty so that the probability of collision is below a specified threshold. In this approach, a new bounding method has been developed to obtain a set of decomposed chance constraints that is a sufficient condition of the original chance constraint. The decomposition of the chance constraint enables its efficient evaluation, as well as the application of the branch and bound method. Branch and bound enables non-convex problems to be solved efficiently to global optimality. Considering the problem of finite-horizon robust optimal control of dynamic systems under Gaussian-distributed stochastic uncertainty, with state and control constraints, a discrete-time, continuous-state linear dynamics model is assumed. Gaussian-distributed stochastic uncertainty is a more natural model for exogenous disturbances such as wind gusts and turbulence than the previously studied set-bounded models. However, with stochastic uncertainty, it is often impossible to guarantee that state constraints are satisfied, because there is typically a non-zero probability of having a disturbance that is large enough to push the state out of the feasible region. An effective framework to address robustness with stochastic uncertainty is optimization with chance constraints. These require that the probability of violating the state constraints (i.e., the probability of failure) is below a user-specified bound known as the risk bound. An example problem is to drive a car to a destination as fast as possible while limiting the probability of an accident to 10(exp -7). This framework allows users to trade conservatism against performance by choosing the risk bound. The more risk the user accepts, the better performance they can expect.

FROM

Characterizing the System Impulse Response Function from Photon-Counting LiDAR Data

NASA's Multiple Altimeter Beam Experimental LiDAR (MABEL) is an aircraft-based photon-counting laser altimeter designed as a simulator to test measurement techniques and algorithms for Advanced Topographic Laser Altimeter System (ATLAS), the sole instrument on NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission. By measuring the time of flight, pointing angle, and absolute position for individual photons, ICESat-2 provides detailed elevation measurements of earth's surface. Calculating accurate and precise elevations requires an understanding of how photons interact with surfaces, and characterization of the photon distribution after returning from surfaces. Neither MABEL nor ATLAS records the transmitted laser pulse shape, relying instead on aggregating several pulses worth of photons, often using histograms, to characterize the pulse shape. In this paper, we assess the limitations of using histograms and propose a more robust method to describe MABEL's system impulse-response function using an exponentially modified Gaussian distribution. We also provide standard error estimates for the arithmetic mean and standard deviation calculations, and for exponentially modified Gaussian parameters using a Monte Carlo sensitivity analysis. We apply this method to photon returns from a sea ice lead and from a dry salt lake bed as case studies for estimating the standard error associated with sample size for the arithmetic mean and standard deviation, and for the exponentially modified Gaussian parameters. We use these standard errors to calculate the minimum number of photons required to find both Gaussian and exponentially modified Gaussian distribution parameters within 3 cm of their parent population values.

photoncounting

A comparison of theoretical and experimental video compression designs

This paper compares theoretical and experimental picture compression designs, for images processed in 8 x 8 blocks using the Walsh-Hadamard transform (WHT). The optimum picture compression design is well known, if the mean-square error (mse) is used as the measure of distortion, and if it is assumed that the video process is a stationary first-order Markov process with a Gaussian distribution. This theoretical design gives useful results when the transform processing is done on full pictures, but gives inferior results (relative to empirical design) when transform processing is done on small 8 x 8 blocks. The use of non-Gaussian distributions for the transform components fails to improve this poor performance, which is due to the nonstationary nature of the video process. An experimentally based design procedure, which considers nonstationarity, yields significantly improved mse and subjective performance.

Jones, H. W., Jr.

Cosmic ray isotope resolution obtained using a new Cerenkov-total energy telescope

The details of a recent measurement of the cosmic ray isotopic composition of Z equals or exceeds 10 nuclei are given, using a new technique combining a Cerenkov and total energy counter. In terms of sigma, the characteristic width of a Gaussian distribution of mass measurement, values of mass resolution have been obtained ranging from about 0.4 AMU at a Z about 14, to about 0.8 AMU at a Z about 26, using the Cerenkov-total energy (C x E) technique. It is customary to claim that two equal Gaussian distributions may be resolved for values of sigma less than 0.4 x the separation of their means. With these criteria isotopes separated by delta A = 2 are being resolved up to a Z about 26 using the C x E technique. An important source of error in these measurements is the determination of the mass of each of the isotopes on an absolute scale.

Webber, W. R.

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian

Performance Metrics for the Assessment of Satellite Data Products: An Ocean Color Case Study

Performance assessment of ocean color satellite data has generally relied on statistical metrics chosen for their common usage and the rationale for selecting certain metrics is infrequently explained. Commonly reported statistics based on mean squared errors, such as the coefficient of determination (r2), root mean square error, and regression slopes, are most appropriate for Gaussian distributions without outliers and, therefore, are often not ideal for ocean color algorithm performance assessment, which is often limited by sample availability. In contrast, metrics based on simple deviations, such as bias and mean absolute error, as well as pair-wise comparisons, often provide more robust and straightforward quantities for evaluating ocean color algorithms with non-Gaussian distributions and outliers. This study uses a SeaWiFS chlorophyll-a validation data set to demonstrate a framework for satellite data product assessment and recommends a multimetric and user-dependent approach that can be applied within science, modeling, and resource management communities.

remote sensing

Bayes classification of terrain cover using normalized polarimetric data

The normalized polarimetric classifier (NPC) which uses only the relative magnitudes and phases of the polarimetric data is proposed for discrimination of terrain elements. The probability density functions (PDFs) of polarimetric data are assumed to have a complex Gaussian distribution, and the marginal PDF of the normalized polarimetric data is derived by adopting the Euclidean norm as the normalization function. The general form of the distance measure for the NPC is also obtained. It is demonstrated that for polarimetric data with an arbitrary PDF, the distance measure of NPC will be independent of the normalization function selected even when the classifier is mistrained. A complex Gaussian distribution is assumed for the polarimetric data consisting of grass and tree regions. The probability of error for the NPC is compared with those of several other single-feature classifiers. The classification error of NPCs is shown to be independent of the normalization function.

Yueh, H. A.

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois

Normal probabilities for Vandenberg AFB wind components - monthly reference periods for all flight azimuths, 0- to 70-km altitudes

Vandenberg Air Force Base (AFB), California, wind component statistics are presented to be used for aerospace engineering applications that require component wind probabilities for various flight azimuths and selected altitudes. The normal (Gaussian) distribution is presented as a statistical model to represent component winds at Vandenberg AFB. Head tail, and crosswind components are tabulated for all flight azimuths for altitudes from 0 to 70 km by monthly reference periods. Wind components are given for 11 selected percentiles ranging from 0.135 percent to 99.865 percent for each month. The results of statistical goodness-of-fit tests are presented to verify the use of the Gaussian distribution as an adequate model to represent component winds at Vandenberg AFB.

Falls, L. W.

Normal probabilities for Cape Kennedy wind components: Monthly reference periods for all flight azimuths. Altitudes 0 to 70 kilometers

This document replaces Cape Kennedy empirical wind component statistics which are presently being used for aerospace engineering applications that require component wind probabilities for various flight azimuths and selected altitudes. The normal (Gaussian) distribution is presented as an adequate statistical model to represent component winds at Cape Kennedy. Head-, tail-, and crosswind components are tabulated for all flight azimuths for altitudes from 0 to 70 km by monthly reference periods. Wind components are given for 11 selected percentiles ranging from 0.135 percent to 99,865 percent for each month. Results of statistical goodness-of-fit tests are presented to verify the use of the Gaussian distribution as an adequate model to represent component winds at Cape Kennedy, Florida.

Falls, L. W.

Intermittent turbulence in solar wind from the south polar hole

The magnetic fields measured by the Ulysses spacecraft are used to study solar wind turbulence in the fast solar wind from the south polar hole. The spacecraft was at about 46 deg south latitude and 3.9 AU. For a magnetic field with a Gaussian distribution the power spectrum (second-order structure function) is sufficient to completely characterize the turbulence. However, the actual distribution is non-Gaussian so that the effects of intermittency must be taken into account. The observed spectral exponents include effects of intermittency and cannot be directly compared with the standard second-order spectral theories such as the Kolmogorov and Kraichnan theories. To permit a better comparison of the observations with the theoretical models, we study the structure characteristics of the data. We find the exponents of the second-order structure functions (power spectra) and the higher-order normalized structure functions for the components of the magnetic fields. We show that these sets of exponents can be approximately described by two basic numbers: the spectral exponent and the intermittency exponent. The intermittency exponent characterizes correlation properties of the energy cascade from large to small scales. Before comparing the observations to the theoretically expected values, a reduction must be made to the observed spectral exponent. The amount of the reduction depends on both the intermittency exponent and the model of the energy cascade assumed in the turbulence theory. We reduce the measured spectral indices according to a simple model for Alfven turbulence that is described here. We then compare our reduced spectral indices with second-order spectral theory. The reduced spectral indices for the period range of 1 min to about a half hour are remarkably constant and in good agreement with the value of 3/2. Thus our treatment is self-consistent. Our tentative conclusion is that the high-frequency turbulence appears to agree with the model of random-phased Alfven waves. This tentative conclusion must be tested by further theoretical and observational work.

Ruzmaikin, A. A.

Distributions of the interplanetary magnetic field revisited

The adequacy of the power spectrum to characterize the variations of a parameter depends on whether or not the parameter has a Gaussian distribution. We here perform very simple tests of Gaussianity on the distribution. We here perform very simple tests of Gaussianity on the distributions of the magnitudes of the interplanetary magnetic field, and on the distributions of the components; that is, we find the first four cumulants of the distributions (mean, variance, skewness, and kurtosis) and their solar cycle variations. We find, consistent with other recent analyses, that the traditional distributions of the 1-hour averaged magnitude are not distributed normally or lognomally as has often been assumed and the 1-hour averaged z component is found to have a nonzero kurtosis. Thus the power spectrum is insufficient to completely characterize these variations and polyspectra are needed. We have isolated variations in the 1/f frequency region of the spectrum and show that the distributions of the magnitudes have nonzero skewness and kurtosis, the magnitudes are not distributed lognormally, and the distributions of the components have nonzero kurtosis. Thus higher-order spectra are again needed for a full characterization.

Feynman, Joan

A Convex Optimization Approach to Improving Suboptimal Hyperparameters of Sliced Normal Distributions

Sliced Normal (SN) distributions are a generalization of Gaussian distributions where the quadratic argument of the exponential is replaced with a sum of squares polynomial. SNs may be used to represent the distribution of a diverse set of random variables including multi-modal, non-symmetric, and skewed distributions. Unfortunately, the likelihood function of a SN includes a normalization constant and the inclusion of this normalization constant makes the likelihood a non-convex function of the hyperparameters which define the SN. In previous work, suboptimal fitting of the hyperparameters was performed by transforming the given data into a higher dimensional monomial basis and selecting the optimal hyperparameters of a Gaussian fit in this space. However, this approach did not account for the effect of lifting on the normalization constant. Indeed, it was observed that as the number of monomials is increased the likelihood of the Sliced Normal can decrease. In this paper, we increase the likelihood of Sliced Normals found using the previous method by developing a convex formulation which scales the covariance matrix of the Gaussian fit such that the likelihood of the Sliced Normal is maximized. The result is significant improvements of the log likelihood of fitted SN distributions, including a significant increase, especially for problems with 500+ monomials.

Convex optimization approach to improving suboptim

Photogrammetric and photometric investigation of a smoke plume viewed from space.

Use of detailed analyses of an Apollo 6 stereographic photograph of a smoke plume which originated in southern Arizona and crossed over into Mexico to illustrate how high-resolution photography can aid meteorologists in evaluating specific air pollution events. Photogrammetric analysis of the visible smoke plume revealed that the plume was 8.06 miles long and attained a maximum width of 4000 ft, 3.0 miles from the 570-ft chimney emitting the effluent. Stereometric analysis showed that the visible top of the plume rose nearly 2400 ft above stack top, attaining 90% of this total rise 1.75 miles downwind from the source. Photometric analysis of the plume revealed a field of plume optical density that portrayed leptokurtic and bimodal distributions rather than a true Gaussian distribution. A horizontal eddy diffusivity of about 650,000 sq cm/sec and a vertical eddy diffusivity of 230,000 sq cm/sec were determined from the plume dimensions. Neutron activation analysis of plume samples revealed the elemental composition of the smoke to be copper, arsenic, selenium, indium and antimony, with trace amounts of vanadium and scandium.

Randerson, D.