Engineering Papers⌕ Search

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Some rules for polydimensional squeezing

The review of the following results is presented: For mixed state light of N-mode electromagnetic field described by Wigner function which has generic Gaussian form, the photon distribution function is obtained and expressed explicitly in terms of Hermite polynomials of 2N-variables. The momenta of this distribution are calculated and expressed as functions of matrix invariants of the dispersion matrix. The role of new uncertainty relation depending on photon state mixing parameter is elucidated. New sum rules for Hermite polynomials of several variables are found. The photon statistics of polymode even and odd coherent light and squeezed polymode Schroedinger cat light is given explicitly. Photon distribution for polymode squeezed number states expressed in terms of multivariable Hermite polynomials is discussed.

Manko, Vladimir I.↗
NASA NTRS1999

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗
NASA NTRS2004

Network Theory: A Primer and Questions for Air Transportation Systems Applications

A new understanding (with potential applications to air transportation systems) has emerged in the past five years in the scientific field of networks. This development emerges in large part because we now have a new laboratory for developing theories about complex networks: The Internet. The premise of this new understanding is that most complex networks of interest, both of nature and of human contrivance, exhibit a fundamentally different behavior than thought for over two hundred years under classical graph theory. Classical theory held that networks exhibited random behavior, characterized by normal, (e.g., Gaussian or Poisson) degree distributions of the connectivity between nodes by links. The new understanding turns this idea on its head: networks of interest exhibit scale-free (or small world) degree distributions of connectivity, characterized by power law distributions. The implications of scale-free behavior for air transportation systems include the potential that some behaviors of complex system architectures might be analyzed through relatively simple approximations of local elements of the system. For air transportation applications, this presentation proposes a framework for constructing topologies (architectures) that represent the relationships between mobility, flight operations, aircraft requirements, and airspace capacity, and the related externalities in airspace procedures and architectures. The proposed architectures or topologies may serve as a framework for posing comparative and combinative analyses of performance, cost, security, environmental, and related metrics.

Holmes, Bruce J.↗
NASA NTRS2009

Gaussian and Lognormal Models of Hurricane Gust Factors

A document describes a tool that predicts the likelihood of land-falling tropical storms and hurricanes exceeding specified peak speeds, given the mean wind speed at various heights of up to 500 feet (150 meters) above ground level. Empirical models to calculate mean and standard deviation of the gust factor as a function of height and mean wind speed were developed in Excel based on data from previous hurricanes. Separate models were developed for Gaussian and offset lognormal distributions for the gust factor. Rather than forecasting a single, specific peak wind speed, this tool provides a probability of exceeding a specified value. This probability is provided as a function of height, allowing it to be applied at a height appropriate for tall structures. The user inputs the mean wind speed, height, and operational threshold. The tool produces the probability from each model that the given threshold will be exceeded. This application does have its limits. They were tested only in tropical storm conditions associated with the periphery of hurricanes. Winds of similar speed produced by non-tropical system may have different turbulence dynamics and stability, which may change those winds statistical characteristics. These models were developed along the Central Florida seacoast, and their results may not accurately extrapolate to inland areas, or even to coastal sites that are different from those used to build the models. Although this tool cannot be generalized for use in different environments, its methodology could be applied to those locations to develop a similar tool tuned to local conditions.

Merceret, Frank↗
DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Probing the large-scale velocity field with clusters of galaxies

What is the role of clusters of galaxies in probing the large-scale velocity field of the universe? We investigate the distribution of peculiar velocities of clusters of galaxies in the popular low-density (omega = 0.3) flat cold dark matter (CDM) cosmological model, which best fits many large-scale structure observations. An omega = 1 CDM model is also studied for comparison. We find that clusters of galaxies are efficient tracers of the large-scale velocity field. The clusters exhibit a Maxwellian distribution of peculiar velocities, as expected from Gaussian initial density fluctuations. The cluster three-dimensional velocity distribution for the omega = 0.3 model peaks at nu approximately greater than 400 km/s and extends to high velocities of nu approximately 1200 km/s. The rms peculiar velocity of the clusters is 440 km/s. Approximately 10% of all model clusters move with high peculiar velocities nu greater or equal to 700 km/s. The observed velocity distribution of clusters of galaxies is compared with the predictions from cosmological models. The observed data exhibit a larger velocity tail than seen in the model simulations; however, due to the large observational uncertainties, the data are consistent at approximately equal to 3 sigma level with the odel predictions, and with a Gaussian initial density field. The large peculiar velocities reported for some clusters of galaxies (nu approximately greater than 3000 km/s) are likely to be overestimated, if the current model is viable.

Bahcall, Neta A.↗
DOE OSTI2020

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Capacity of PPM on Gaussian and Webb Channels

This paper computes and compares the capacities of M-ary PPM on various idealized channels that approximate the optical communication channel: (1) the standard additive white Gaussian noise (AWGN) channel;(2) a more general AWGN channel (AWGN2) allowing different variances in signal and noise slots;(3) a Webb-distributed channel (Webb2);(4) a Webb+Gaussian channel, modeling Gaussian thermal noise added to Webb-distributed channel outputs.

Capacity PPM Webb Optical↗
NASA NTRS1995

Inverse sequential procedures for the monitoring of time series

When one or more new values are added to a developing time series, they change its descriptive parameters (mean, variance, trend, coherence). A 'change index (CI)' is developed as a quantitative indicator that the changed parameters remain compatible with the existing 'base' data. CI formulate are derived, in terms of normalized likelihood ratios, for small samples from Poisson, Gaussian, and Chi-Square distributions, and for regression coefficients measuring linear or exponential trends. A substantial parameter change creates a rapid or abrupt CI decrease which persists when the length of the bases is changed. Except for a special Gaussian case, the CI has no simple explicit regions for tests of hypotheses. However, its design ensures that the series sampled need not conform strictly to the distribution form assumed for the parameter estimates. The use of the CI is illustrated with both constructed and observed data samples, processed with a Fortran code 'Sequitor'.

Radok, Uwe↗
NASA NTRS1992

Interpretation of lightcurves of atmosphereless bodies. I - General theory and new inversion schemes

Theoretical models of light-curve inversion that can be used in photomorphography, i.e., determination of the 3D shape and/or the light-scattering behavior of the surface of a body from disk-integrated photometry. These methods can be applied to atmosphereless bodies, mainly asteroids and planetary satellites, in the solar system. The inversion comprises three steps. First, a function containing data on the shape and the albedo variegation of an asteroid is determined. In the second step, separate expressions for the inverse of the Gaussian curvature and the albedo distribution are derived from the information obtained in the first step. In the third step, the nontrivial problem of determining the radius vector of the surface from the Gaussian curvature is solved by using iterative optimization procedures.

Kaasalainen, M.↗
NASA NTRS1994

1919+479: Big WAT in a poor cluster

New x-ray, optical, and redshift data are presented for the cluster of galaxies associated with the giant, 1 Mpc diameter, wide-angle tailed (WAT) radio galaxy 1919+479. The ROSAT Position-Sensitive Proportional Counter (PSPC) pointed observation shows an x-ray peak on the WAT and elongated diffuse emission tracing the galaxy distribution. In addition, an asymmetric extension of emission exists between the tails of the WAT. The fitting of a Raymond-Smith thermal model to the x-ray spectra suggests an approximately = 2 keV temperature intracluster medium (ICM). The cooling time and irregular morphology rule out a cluster-wide cooling flow. The x-ray luminosity and temperature are consistent with the velocity dispersion, 480 km/s, estimated from 31 galaxy velocities. However, this velocity distribution is significantly non-Gaussian, which along with the x-ray morphology, suggests incomplete virialization in the cluster. Substructure analysis does not reveal significant clumping in the velocities/positions; but, the spatial distribution of galaxies is very elongated. Also, the cD galaxy producing the WAT does not have a significant radial peculiar velocity with respect to the cluster centroid. These characteristics are consistent with a merger scenario in which a subcluster has crossed the cluster core in the plane of the sky and has dispersed. We compare this cluster with the post-merger cluster Abell 2634 (Pinkney et al., 1993), containing the prototype WAT, and with a recent N-Body/Hydro simulation of merging clusters. The similarities indicate that the cluster 1919+479 may be the poor extreme of a class of clusters in which the bulk motion in the ICM, caused by a subcluster merger, in shaping the central WAT.

Pinkney, Jason↗
DOE OSTI2020

Mass-asymmetric fission of 205,207,209 Bi at energies close to the fission barrier using proton bombardment of 204,206,208 Pb

Background: Recent observation of mass-asymmetric fission in neutron-deficient Hg and Pt nuclei has reignited interest in fission fragment mass distributions close to Pb. Investigations at energies close to the fission barrier, where mass-asymmetric fission is expected to be most obvious and the sensitivity to shell effects is maximized, are limited in this mass region. Purpose: To measure fission mass distributions for 205,207,209 Bi nuclei at the lowest possible excitation energies to determine how the mass distributions change with excitation energy and the neutron number of the compound nucleus. Method: Proton beams bombarding targets of 204,206,208 Pb were used to study the fission of 205,207,209 Bi at energies from just above to 10 MeV above their fission barriers. Fission fragments were measured using the CUBE fission spectrometer. Fission fragment mass distributions were determined using a newly developed time difference analysis method. Mass distributions were characterized by triple-Gaussian fits to determine the systematic trends across each isotope with excitation energy. Results: Measured mass distributions of all three Bi isotopes exhibit a component of mass-asymmetric fission at all energies studied. The probability of mass-asymmetric fission decreases significantly with increasing excitation energy, from ≈70 to ≈40% over a 10-MeV range. Comparisons between the three Bi isotopes hint at an increase in the mass-symmetric fission yield with increasing neutron number, which could be due to a decrease in the difference between the symmetric and asymmetric fission barriers. The centroids of the mass-asymmetric peaks suggest that several deformed shell gaps in the fission fragments could be contributing to the presence of the mass-asymmetric fission mode with Z light ≃ 38, Z heavy ≃ 45, and N light ≃ 56 all present in the fission fragments. Conclusions: Measurements of fission mass distributions at the lowest possible excitation energies above the fission barrier provide an excellent platform to investigate the origins of the mass-asymmetric fission mode. Finally, further systematic measurements at these energies offer an opportunity to rigorously test new models of fission in this mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Theoretical investigation of gas-surface interactions

A Dirac-Hartree-Fock code was developed for polyatomic molecules. The program uses integrals over symmetry-adapted real spherical harmonic Gaussian basis functions generated by a modification of the MOLECULE integrals program. A single Gaussian function is used for the nuclear charge distribution, to ensure proper boundary conditions at the nuclei. The Gaussian primitive functions are chosen to satisfy the kinetic balance condition. However, contracted functions which do not necessarily satisfy this condition may be used. The Fock matrix is constructed in the scalar basis and transformed to a jj-coupled 2-spinor basis before diagonalization. The program was tested against numerical results for atoms with a Gaussian nucleus and diatomic molecules with point nuclei. The energies converge on the numerical values as the basis set size is increased. Full use of molecular symmetry (restricted to D sub 2h and subgroups) is yet to be implemented.

Dyall, Kenneth G.↗
DOE OSTI2023

Robust neural network-enhanced estimation of local primordial non-Gaussianity

When applied to the nonlinear matter distribution of the universe, neural networks have been shown to be very statistically sensitive probes of cosmological parameters, such as the linear perturbation amplitude σ 8 . However, when used as a “black box,” neural networks are not robust to baryonic uncertainty. We propose a robust architecture for constraining primordial non-Gaussianity f NL , by training a neural network to locally estimate σ 8 , and correlating these local estimates with the large-scale density field. We apply our method to N-body simulations, and show that σ⁡(f NL ) is 3.5 times better than the constraint obtained from a standard halo-based approach. Finally, we show that our method has the same robustness property as large-scale halo bias: baryonic physics can change the normalization of the estimated f NL , but cannot change whether f NL is detected.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Electron velocity distributions near interplantary shocks

Previous work on electron heating mechanisms operating at the supercritical earth's bow shock is extended to the generally lower Mach number interplanetary shocks. This is done by studying electron velocity distributions observed both upstream and downstream of interplanetary shocks observed by ISEE 3 between August 1978 and December 1979. It is found that perpendicular heating accompanied by little or no parallel heating occurs at the weaker shocks but that parallel heating is dominant for the stronger shocks. When perpendicular heating is dominant, downstream velocity distributions are for the most part Gaussian at low energies (core) and Lorentzian at high energies (halo). When parallel heating is dominant, however, these distributions are flat-topped at low energies, having power law extensions to higher energies. At the weaker shocks, the observations are consistent with electron heating resulting from conservation of the magnetic moment, whereas at the stronger shocks the heating is quite similar to that observed at the earth's bow shock.

Feldman, W. C.↗
NASA NTRS2016

Erasing the Milky Way: New Cleaning Technique Applied to GBT Intensity Mapping Data

We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15- and 1-h field data of the GBT observations previously presented in Masui et al. (2013) and Switzer et al. (2013), covering about 41 square degrees at 0.6 less than z is less than 1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative sub-traction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.

Wolz, L.↗
NASA NTRS2012

Effects of Tropospheric Spatio-Temporal Correlated Noise on the Analysis of Space Geodetic Data

The standard VLBI analysis models the distribution of measurement noise as Gaussian. Because the price of recording bits is steadily decreasing, thermal errors will soon no longer dominate. As a result, it is expected that troposphere and instrumentation/clock errors will increasingly become more dominant. Given that both of these errors have correlated spectra, properly modeling the error distributions will become increasingly relevant for optimal analysis. We discuss the advantages of modeling the correlations between tropospheric delays using a Kolmogorov spectrum and the frozen flow assumption pioneered by Treuhaft and Lanyi. We then apply these correlated noise spectra to the weighting of VLBI data analysis for two case studies: X/Ka-band global astrometry and Earth orientation. In both cases we see improved results when the analyses are weighted with correlated noise models vs. the standard uncorrelated models. The X/Ka astrometric scatter improved by approx.10% and the systematic Delta delta vs. delta slope decreased by approx. 50%. The TEMPO Earth orientation results improved by 17% in baseline transverse and 27% in baseline vertical.

instrumentation errors↗
DOE OSTI2024

Statistical Behavior of Low-Amplitude Power System Point-on-Wave Measurements

The power grid is undergoing massive changes to ensure resiliency and reliability in a more decentralized world. Distributed energy resources are becoming a prominent source of generation, potentially leading to a lack of centralized generation sources. Due to these new behaviors and system topologies, it is important to install measurement devices that are 1) accurate and 2) self-aware of their measurement quality. In this paper, a residential-scale microgrid is used to generate voltage and current waveforms, captured by Verivolt and National Instruments measurement equipment. A least-squares approach is used to separate the “clean” signals from the noise. Finally, Gaussian mixture modeling is used to approximate noise distributions, and it is shown these higher-order distribution estimates are a better fit to voltage and current noise profiles than single-mode Gaussian estimates.

24 POWER TRANSMISSION AND DISTRIBUTION↗