Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gaussian Approximation”

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 343 records · Page 19

Cosmological constraints from the Subaru Hyper Suprime-Cam year 1 shear catalogue lensing convergence probability distribution function

Here we utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) year 1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude S 8 = σ 8 $\sqrt {Ω_{m}/0.3}$. The large-scale structure’s statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF—a complementary higher-order statistic. By defining the PDF over the standard-deviation-normalized convergence map, we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find S 8 = $0.860^{+0.066}_{–0.109}$ from the PDF alone and S 8 = 0.798$^{+0.029}_{–0.042}$ from the combination of the PDF and power spectrum (68% confidential level (CL)). The PDF improves the power-spectrum-only constraint by about 10%.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffuse reflectance of the ocean - The theory of its augmentation by chlorophyll a fluorescence at 685 nm

The radiative transfer equation is modified to include the effect of fluorescent substances and solved in the quasi-single scattering approximation for a homogeneous ocean containing fluorescent particles with wavelength independent quantum efficiency and a Gaussian shaped emission line. The results are applied to the in vivo fluorescence of chlorophyll a (in phytoplankton) in the ocean to determine if the observed quantum efficiencies are large enough to explain the enhancement of the ocean's diffuse reflectance near 685 nm in chlorophyll rich waters without resorting to anomalous dispersion. The computations indicate that the required efficiencies are sufficiently low to account completely for the enhanced reflectance. The validity of the theory is further demonstrated by deriving values for the upwelling irradiance attenuation coefficient at 685 nm which are in close agreement with the observations.

Gordon, H. R.↗

The possible modifications of the HISSE model for pure LANDSAT agricultural data

A method for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field) is discussed. Theoretical arguments are given which show that any significant refinement of the model beyond this proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Gaussian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗

Probability of Obstacle Collision for UAVs in Presence of Wind

For incorporation of unmanned aerial vehicles into the National Airspace, ensuring safety of the airspace including the vehicles, people, and property on the ground is of utmost importance. One of the safety-critical factors for unmanned aviation flights is the risk of deviating from a planned trajectory resulting in a variety of hazards, including potential loss of separation between vehicle and obstacles or unexpected battery energy consumption. Off-nominal conditions introduced by component failures, degraded controllability and environmental disturbances such as wind gusts can lead to unacceptable unexpected deviations from the flight trajectory. It is essential to accurately model such effects on the flight trajectory while computing safety thresholds such as minimum separation from surrounding obstacles, available battery resource to complete the mission or determining delay in the expected time of arrival of flights. In this paper, a tool is presented based on Gaussian Process Regression for wind representation over a pre-defined trajectory for fast, yet approximated, in-time evaluation of possible trajectory deviations caused by wind gusts. The deviation in the planned trajectory caused by wind is further simulated utilizing a 6 degrees-of-freedom (DOF) UAV trajectory simulator comprising of a rotorcraft lumped-mass model with LQRI controller. Both steady-state wind and wind gust effects are investigated. The probability of collision with obstacle is computed and demonstrated on real flight data from experimental flights of an octocopter at NASA Langley Research Center in the presence of simulated obstacles and wind conditions. Effect of varying wind conditions and varying UAV airspeed is further demonstrated on experimental flights in the presence of wind measured by ground based weather service stations. The proposed approach would eventually benefit timely mitigation of current and future safety-critical events in autonomous systems by enabling risk-informed decision making.

Portia Banerjee↗

Codiscovering graphical structure and functional relationships within data: A Gaussian Process framework for connecting the dots

Most problems within and beyond the scientific domain can be framed into one of the following three levels of complexity of function approximation. Type 1: Approximate an unknown function given input/output data. Type 2: Consider a collection of variables and functions, some of which are unknown, indexed by the nodes and hyperedges of a hypergraph (a generalized graph where edges can connect more than two vertices). Given partial observations of the variables of the hypergraph (satisfying the functional dependencies imposed by its structure), approximate all the unobserved variables and unknown functions. Type 3: Expanding on Type 2, if the hypergraph structure itself is unknown, use partial observations of the variables of the hypergraph to discover its structure and approximate its unknown functions. These hypergraphs offer a natural platform for organizing, communicating, and processing computational knowledge. While most scientific problems can be framed as the data-driven discovery of unknown functions in a computational hypergraph whose structure is known (Type 2), many require the data-driven discovery of the structure (connectivity) of the hypergraph itself (Type 3). We introduce an interpretable Gaussian Process (GP) framework for such (Type 3) problems that does not require randomization of the data, access to or control over its sampling, or sparsity of the unknown functions in a known or learned basis. Its polynomial complexity, which contrasts sharply with the super-exponential complexity of causal inference methods, is enabled by the nonlinear ANOVA capabilities of GPs used as a sensing mechanism.

Science & Technology - Other Topics↗

Analysis of low altitude atmospheric turbulence data measured in flight

All three components of turbulence were measured simultaneously in flight at each wing tip of a Beech D-18 aircraft. The flights were conducted at low altitude, 30.5 - 61.0 meters (100-200 ft.), over water in the presence of wind driven turbulence. Statistical properties of flight measured turbulence were compared with Gaussian and non-Gaussian turbulence models. Spatial characteristics of the turbulence were analyzed using the data from flight perpendicular and parallel to the wind. The probability density distributions of the vertical gusts show distinctly non-Gaussian characteristics. The distributions of the longitudinal and lateral gusts are generally Gaussian. The power spectra compare in the inertial subrange at some points better with the Dryden spectrum, while at other points the von Karman spectrum is a better approximation. In the low frequency range the data show peaks or dips in the power spectral density. The cross between vertical gusts in the direction of the mean wind were compared with a matched non-Gaussian model. The real component of the cross spectrum is in general close to the non-Gaussian model. The imaginary component, however, indicated a larger phase shift between these two gust components than was found in previous research.

Ganzer, V. M.↗

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗

Constraints on primordial non-Gaussianity from the cross-correlation of DESI luminous red galaxies and Planck CMB lensing

We use the angular cross-correlation between a luminous red galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey data release DR9 and the Planck cosmic microwave background (CMB) lensing maps to constrain the local primordial non-Gaussianity parameter, f NL , using the scale-dependent galaxy bias effect. The galaxy sample covers approximately 40% of the sky, contains galaxies up to redshift z ∼ 1.4, and is calibrated with the LRG spectra that have been observed for DESI Year 1 (Y1). We apply a nonlinear imaging systematics treatment based on neural networks to remove observational effects that could potentially bias the f NL measurement. Our measurement is performed without blinding, but the full analysis pipeline is tested with simulations including systematics. Using the two-point angular cross-correlation between LRG and CMB lensing only, we find f NL = 39 −38 +40 at the 68% confidence level, and our result is robust in terms of systematics and cosmological assumptions. If we combine this information with the autocorrelation of LRG, applying a scale cut to limit the impact of systematics, we find f NL = 24 −21 +20 at the 68% confidence level. Our results motivate the use of CMB lensing cross-correlations to measure f NL with future datasets, given its stability in terms of observational systematics compared to the angular autocorrelation. Furthermore, performing accurate systematics mitigation is crucially important in order to achieve competitive constraints on f NL from CMB lensing cross-correlation in combination with the tracers’ autocorrelation.

79 ASTRONOMY AND ASTROPHYSICS↗

Hierarchical ensemble Kalman methods with sparsity-promoting generalized gamma hyperpriors

This paper introduces a computational framework to incorporate flexible regularization techniques in ensemble Kalman methods, generalizing the iterative alternating scheme to nonlinear inverse problems. The proposed methodology approximates the maximum a posteriori (MAP) estimate of a hierarchical Bayesian model characterized by a conditionally Gaussian prior and generalized gamma hyperpriors. Suitable choices of hyperparameters yield sparsity-promoting regularization. We propose an iterative algorithm for MAP estimation, which alternates between updating the unknown with an ensemble Kalman method and updating the hyperparameters in the regularization to promote sparsity. Here, the effectiveness of our methodology is demonstrated in several computed examples, including compressed sensing and subsurface flow inverse problems.

Ensemble Kalman methods↗

Uncertainty Quantification via Stable Distribution Propagation

We propose a new approach for propagating stable probability distributions through neural networks. Our method is based on local linearization, which we show to be an optimal approximation in terms of total variation distance for the ReLU non-linearity. This allows propagating Gaussian and Cauchy input uncertainties through neural networks to quantify their output uncertainties. To demonstrate the utility of propagating distributions, we apply the proposed method to predicting calibrated confidence intervals and selective prediction on out-of-distribution data. The results demonstrate a broad applicability of propagating distributions and show the advantages of our method over other approaches such as moment matching.

Artificial Intelligence (cs.AI)↗

On the first-excursion probability in stationary narrow-band random vibration. II.

The first-excursion probability of a stationary narrow-band Gaussian process with mean zero has been studied. Within the framework of point process approach, series approximations derived from the theory of random points and approximations based on the maximum entropy principle have been developed. With the aid of numerical examples, merits of the approximations proposed previously as well as of those developed in this paper have been compared. The results indicate that the maximum entropy principle has not produced satisfactory approximations but the approximation based on nonapproaching random points is found to be the best among all the approximations proposed herein. A conclusion drawn from the present and the previous studies is that the point process approach produces a number of useful approximations for the first-excursion probability, particularly those based on the concepts of the Markov process, the clump-size, and the nonapproaching random points.

Yang, J.-N.↗

Soot and hydrocarbon formation in a turbulent diffusion flame

Formation and emission of soot and polycyclic aromatic hydrocarbons (PCAH) from a turbulent continuous flow combustor have been studied. Measurements included mass concentration of both soot and PCAH, composition of individual PCAH, and size distribution of soot particles, as a function of mixing intensity, fuel equivalence ratio and type of fuel (kerosene or benzene). Both soot and PCAH concentrations reach maxima early in the flame, after which PCAH decays rapidly and soot decays much slower. The maximum PCAH concentration always preceded that of soot in agreement with the concept that certain PCAH may serve as intermediates in soot formation. An approximate calculation based on the assumptions of local equilibrium with respect to soot formation and a Gaussian distribution of air-fuel mixedness gave satisfactory correlations of the data on soot formation at the higher cold gas velocity. At the lower velocity, the amount of soot and PCAH formed was drastically increased and strongly dependent on fuel atomization.

Prado, G. P.↗

Modern CACSD using the Robust-Control Toolbox

The Robust-Control Toolbox is a collection of 40 M-files which extend the capability of PC/PRO-MATLAB to do modern multivariable robust control system design. Included are robust analysis tools like singular values and structured singular values, robust synthesis tools like continuous/discrete H(exp 2)/H infinity synthesis and Linear Quadratic Gaussian Loop Transfer Recovery methods and a variety of robust model reduction tools such as Hankel approximation, balanced truncation and balanced stochastic truncation, etc. The capabilities of the toolbox are described and illustated with examples to show how easily they can be used in practice. Examples include structured singular value analysis, H infinity loop-shaping and large space structure model reduction.

Chiang, Richard Y.↗

A scattering model for ocean surface

A surface-scattering model based on an approximate solution of the integral equations for the surface tangential fields is developed for non-Gaussian distributed finitely conducting surfaces. It is a function of both the surface spectrum and the surface bispectrum which are defined to be the Fourier transforms of the surface-correlation function and the surface-skewness function, respectively. By employing a sea-surface spectrum proposed by Pierson (1964) and modified by Lee and Fung (1982) it is shown that good agreement is obtained between model and field measurements in polarization, incidence angle, and azimuth angle without adjusting surface parameters.

Chen, K. S.↗

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗

Characterization of the degradation of gamma-irradiated elastomers using Raman spectroscopy

This report presents key findings from Raman spectroscopic analysis of gamma-irradiated rubber samples extracted from a laminated lead-damped rubber (LDR) seismic isolation device. The samples were exposed to gamma radiation from a 60Co source in a Foss Therapy Services gamma irradiator, reaching absorbed doses up to 1600 kGy. A distinct threshold near 400 kGy was identified, beyond which significant spectral changes were observed. Two Raman peaks - at approximately 425 cm-1 and 2440 cm-1 - were tracked as a function of dose using Gaussian fitting. The 425 cm-1 peak, attributed to sulfur–sulfur (S–S) bond stretching (resulting from vulcanization of the rubber), exhibited a dose-dependent upshift, indicating radiation-induced crosslinking within the sulfur-based polymer network. Conversely, the 2440 cm-1 peak, likely associated with vibrational modes of additives or impurities, showed a downward shift with increasing dose, suggesting chain scission and degradation of non-rubber constituents. These results provide first-of-a-kind insights into the microstructural evolution of elastomers under high-dose gamma irradiation and establish a preliminary dose threshold for significant degradation. Future work will incorporate multi-modal characterization—including Fourier Transform Infrared (FTIR) spectroscopy, scanning electrom microscopy (SEM) of the rubber surface morphology, thermogravimetric analysis (TGA) to determine changes in thermal stability, and mechanical testing—to correlate molecular-level changes with macroscopic performance of these elastomers as damping media in seismic isolation devices. These findings are expected to provide regulatory guidance and design criteria for qualifying low-damping rubber seismic isolators in advanced nuclear reactor applications.

36 - MATERIALS SCIENCE↗

Probability density and exceedance rate functions of locally Gaussian turbulence

A locally Gaussian model of turbulence velocities is postulated which consists of the superposition of a slowly varying strictly Gaussian component representing slow temporal changes in the mean wind speed and a more rapidly varying locally Gaussian turbulence component possessing a temporally fluctuating local variance. Series expansions of the probability density and exceedance rate functions of the turbulence velocity model, based on Taylor's series, are derived. Comparisons of the resulting two-term approximations with measured probability density and exceedance rate functions of atmospheric turbulence velocity records show encouraging agreement, thereby confirming the consistency of the measured records with the locally Gaussian model. Explicit formulas are derived for computing all required expansion coefficients from measured turbulence records.

Mark, W. D.↗

The luminosity function of quasars

We propose a new evolutionary model for the optical luminosity function of quasars. Our analytical model is derived from fits to the empirical luminosity function estimated by Hartwick and Schade and Warren, Hewett, and Osmer on the basis of more than 1200 quasars over the range of redshifts 0 approximately less than z approximately less than 4.5. We find that the evolution of quasars over this entire redshift range can be well fitted by a Gaussian distribution, while the shape of the luminosity function can be well fitted by either a double power law or an exponential L(exp 1/4) law. The predicted number counts of quasars, as a function of either apparent magnitude or redshift, are fully consistent with the observed ones. Our model indicates that the evolution of quasars reaches its maximum at z approximately = 2.8 and declines at higher redshifts. An extrapolation of the evolution to z approximately greater than 4.5 implies that quasars may have started their cosmic fireworks at z(sub f) approximately = 5.2-5.5. Forthcoming surveys of quasars at these redshifts will be critical to constrain the epoch of quasar formation. All the results we derived are based on observed quasars and are therefore subject to the bias of obscuration by dust in damped Ly alpha systems. Future surveys of these absorption systems at z approximately greater than 3 will also be important if the formation epoch of quasars is to be known unambiguously.

Pei, Yichuan C.↗