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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Adaptive Key Component Control of Nonlinear Evolving Flexible Structures

An optimal alarm system and its approximations may use Kalman filtering for univariate linear dynamic systems driven by Gaussian noise to provide a layer of predictive capability. Predicted Kalman filter future process values and a fixed critical threshold can be used to construct a candidate level-crossing event over a predetermined prediction window. An optimal alarm system can be designed to elicit the fewest false alarms for a fixed detection probability in this particular scenario.

Balas, Mark↗
NASA NTRS1988

Focal region fields of distorted reflectors

The problem of the focal region fields scattered by an arbitrary surface reflector under uniform plane wave illumination is solved. The physical optics (PO) approximation is used to calculate the current induced on the reflector. The surface of the reflector is described by a number of triangular domain-wise 5th degree bivariate polynomials. A 2-dimensional Gaussian quadrature is employed to numerically evaluate the integral expressions of the scattered fields. No Freshnel or Fraunhofer zone approximations are made. The relation of the focal fields problem to surface compensation techniques and other applications are mentioned. Several examples of distorted parabolic reflectors are presented. The computer code developed is included, together with instructions on its usage.

Buris, N. E.↗
NASA NTRS1993

A function approximation approach to anomaly detection in propulsion system test data

Ground test data from propulsion systems such as the Space Shuttle Main Engine (SSME) can be automatically screened for anomalies by a neural network. The neural network screens data after being trained with nominal data only. Given the values of 14 measurements reflecting external influences on the SSME at a given time, the neural network predicts the expected nominal value of a desired engine parameter at that time. We compared the ability of three different function-approximation techniques to perform this nominal value prediction: a novel neural network architecture based on Gaussian bar basis functions, a conventional back propagation neural network, and linear regression. These three techniques were tested with real data from six SSME ground tests containing two anomalies. The basis function network trained more rapidly than back propagation. It yielded nominal predictions with, a tight enough confidence interval to distinguish anomalous deviations from the nominal fluctuations in an engine parameter. Since the function-approximation approach requires nominal training data only, it is capable of detecting unknown classes of anomalies for which training data is not available.

Whitehead, Bruce A.↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Small-Scale Gravity Waves in ER-2 MMS/MTP Wind and Temperature Measurements during CRYSTAL-FACE

Lower stratospheric wind and temperature measurements made from NASA's high-altitude ER-2 research aircraft during the CRYSTAL-FACE campaign in July 2002 were analyzed to retrieve information on small scale gravity waves (GWs) at the aircraft's flight level (typically approximately 20 km altitude). For a given flight segment, the S-transform (a Gaussian wavelet transform) was used to search for and identify small horizontal scale GW events, and to estimate their apparent horizontal wavelengths. The horizontal propagation directions of the events were determined using the Stokes parameter method combined with the cross S-transform analysis. The vertical temperature gradient was used to determine the vertical wavelengths of the events. GW momentum fluxes were calculated from the cross S-transform. Other wave parameters such as intrinsic frequencies were calculated using the GW dispersion relation. More than 100GW events were identified. They were generally high frequency waves with vertical wavelength of approximately 5 km and horizontal wavelength generally shorter than 20 km. Their intrinsic propagation directions were predominantly toward the east, whereas their ground-based propagation directions were primarily toward the west. Among the events, approximately 20% of them had very short horizontal wavelength, very high intrinsic frequency, and relatively small momentum fluxes, and thus they were likely trapped in the lower stratosphere. Using the estimated GW parameters and the background winds and stabilities from the NCAR/NCEP reanalysis data, we were able to trace the sources of the events using a simple reverse ray-tracing. More than 70% of the events were traced back to convective sources in the troposphere, and the sources were generally located upstream of the locations of the events observed at the aircraft level. Finally, a probability density function of the reversible cooling rate due to GWs was obtained in this study, which may be useful for cirrus cloud models.

Wang, L.↗
NASA NTRS1976

Development and application of an atmospheric turbulence model for use in flight simulators in flight simulators

The influence of simulated turbulence on aircraft handling qualities was investigated. Pilot opinion of the handling qualities of a light general aviation aircraft were evaluated in a motion-base simulator using a simulated turbulence environment. A realistic representation of turbulence disturbances is described in terms of rms intensity and scale length and their random variations with time. The time histories generated by the proposed turbulence models showed characteristics which appear to be more similar to real turbulence than the frequently-used Gaussian turbulence model. In addition, the proposed turbulence models can flexibly accommodate changes in atmospheric conditions and be easily implemented in flight simulator studies. Six turbulence time histories, including the conventional Gaussian model, were used in an IFR-tracking task. The realism of each of the turbulence models and the handling qualities of the simulated airplane were evaluated. Analysis of pilot opinions shows that at approximately the same rms intensities of turbulence, the handling quality ratings transit from the satisfactory level, for the simple Gaussian model, to an unacceptable level for more realistic and compositely structured turbulence models.

Jacobson, I. D.↗
NASA NTRS2011

Deep UV Luminosity Functions at the Infall Region of the Coma Cluster

We have used deep GALEX observations at the infall region of the Coma cluster to measure the faintest UV luminosity functions (LFs) presented for a rich galaxy cluster thus far. The Coma UV LFs are measured to M(sub uv) = -10.5 in the GALEX FUV and NUV bands, or 3.5 mag fainter than previous studies, and reach the dwarf early-type galaxy population in Coma for the first time. The Schechter faint-end slopes (alpha approximately equal to -1.39 in both GALEX bands) are shallower than reported in previous Coma UV LF studies owing to a flatter LF at faint magnitudes. A Gaussian-plus-Schechter model provides a slightly better parametrization of the UV LFs resulting in a faint-end slope of alpha approximately equal to -1.15 in both GALEX bands. The two-component model gives faint-end slopes shallower than alpha = -1 (a turnover) for the LFs constructed separately for passive and star forming galaxies. The UV LFs for star forming galaxies show a turnover at M(sub UV) approximately equal to -14 owing to a deficit of dwarf star forming galaxies in Coma with stellar masses below M(sub *) = 10(sup 8) solar mass. A similar turnover is identified in recent UV LFs measured for the Virgo cluster suggesting this may be a common feature of local galaxy clusters, whereas the field UV LFs continue to rise at faint magnitudes. We did not identify an excess of passive galaxies as would be expected if the missing dwarf star forming galaxies were quenched inside the cluster. In fact, the LFs for both dwarf passive and star forming galaxies show the same turnover at faint magnitudes. We discuss the possible origin of the missing dwarf star forming galaxies in Coma and their expected properties based on comparisons to local field galaxies.

Hammer, D. M.↗
DOE OSTI2024

Projective Integral Updates for High-Dimensional Variational Inference

Variational inference is an approximation framework for Bayesian inference that seeks to improve quantified uncertainty in predictions by optimizing a simplified distribution over parameters to stand in for the full posterior. Capturing model variations that remain consistent with training data enables more robust predictions by reducing parameter sensitivity. This work introduces a fixed-point optimization for variational inference that is applicable when every feasible log density can be expressed as a linear combination of functions from a given basis. In such cases, the optimizer becomes a fixed-point of projective integral updates. When the basis spans univariate quadratics in each parameter, the feasible distributions are Gaussian mean-fields and the projective integral updates yield quasi-Newton variational Bayes (QNVB). Other bases and updates are also possible. Since these updates require high-dimensional integration, this work begins by proposing an efficient quasirandom sequence of quadratures for mean-field distributions. Each iterate of the sequence contains two evaluation points that combine to correctly integrate all univariate quadratic functions and, if the mean-field factors are symmetric, all univariate cubics. More importantly, averaging results over short subsequences achieves periodic exactness on a much larger space of multivariate polynomials of quadratic total degree. The corresponding variational updates require four loss evaluations with standard (not second-order) backpropagation to eliminate error terms from over half of all multivariate quadratic basis functions. Furthermore, this integration technique is motivated by first proposing stochastic blocked mean-field quadratures, which may be useful in other contexts. A PyTorch implementation of QNVB allows for better control over model uncertainty during training than competing methods. Experiments demonstrate superior generalizability for multiple learning problems and architectures.

Gaussian mean-field↗
DOE OSTI2022

Diffraction and thermal effect of a Bessel-Gaussian laser for Ag nanoparticle deposition

Nanoparticles are known to sinter at much lower temperatures than the corresponding bulk or micro size particles. A laser-assisted sintering process is considered in this study to sinter Ag nanoparticles by dispensing Ag paste onto an indium tin oxide-coated Si substrate. The Gaussian beam of a CO 2 laser source is propagated through axicon and biconvex lenses, and the resulting hollow beam is focused on the Ag paste with a hollow parabolic mirror. A Bessel-Gaussian irradiance distribution is obtained at the focal plane of the parabolic mirror due to the interference of the hollow laser cone. The Fresnel diffraction approximation is considered to determine the phasor of the laser and an analytical approach is implemented to calculate the irradiance distribution of the Bessel-Gaussian beam. This irradiance distribution is utilized as a heat source in a heat conduction model and the temperature distribution is analyzed for thin Ag films formed during the laser sintering of Ag nanoparticles. An analytical expression is obtained for the temperature distribution by solving the heat conduction equation using Fourier transform for finite media. The widths of the deposited Ag lines are predicted from the temperature profiles and the model predictions compare well with the experimental results. The isotherms are found to be geometrically noncongruent with convex and concave tips depending on the locally maximum and minimum irradiances of the Bessel-Gaussian beam, respectively. The convex and concave tips, however, appear in the same isotherm for sufficiently high substrate speed relative to the laser beam.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1992

Topology of large-scale structure in seeded hot dark matter models

The topology of the isodensity surfaces in seeded hot dark matter models, in which static seed masses provide the density perturbations in a universe dominated by massive neutrinos is examined. When smoothed with a Gaussian window, the linear initial conditions in these models show no trace of non-Gaussian behavior for r0 equal to or greater than 5 Mpc (h = 1/2), except for very low seed densities, which show a shift toward isolated peaks. An approximate analytic expression is given for the genus curve expected in linear density fields from randomly distributed seed masses. The evolved models have a Gaussian topology for r0 = 10 Mpc, but show a shift toward a cellular topology with r0 = 5 Mpc; Gaussian models with an identical power spectrum show the same behavior.

Beaky, Matthew M.↗
NASA NTRS1994

New description of charged particle propagation in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by random components that are superposed on the guiding field. In the simplest analysis of this situation, scattering causes the particles to diffuse parallel to the guiding field. At the next level of approximation, moving pulses that correspond to a coherent mode of propagation are present, but they are represented by delta-functions whose infinitely narrow width makes no sense physically and is inconsistent with the finite duration of coherent pulses observed in solar energetic particle events. To derive a more realistic description, the transport problem is formulated in terms of 4 x 4 matrices, which derive from a representation of the particle distribution function in terms of eigenfunctions of the scattering operator, and which lead to useful approximations that give explicit predictions of the detailed evolution not only of the coherent pulses, but also of the diffusive wake. More specifically, the new description embodies a simple convolution of a narrow Gaussian with the solutions above that involve delta-functions, but with a slightly reduced coherent velocity. The validity of these approximations, which can easily be calculated on a desktop computer, has been exhaustively confirmed by comparison with results of Monte Carlo simulations which kept track of 50 million particles and which were carried out on the Maspar computer at Goddard Space Flight Center.

Earl, James A.↗
NASA NTRS2020

Electromagnetic scattering by discrete random media illuminated by a Gaussian beam I: Derivation of the radiative transfer equation

In this paper we present the vector radiative transfer theory for a discrete random medium illuminated by a Gaussian beam. The analysis is based on a plane wave spectrum representation for a Gaussian beam and uses an approach developed previously for a discrete random medium illuminated by a plane electromagnetic wave. Specifically, we establish an integral representation for the coherent field, define an approximate coherent field that satisfies the differential equation fulfilled by the coherent field corresponding to a plane electromagnetic wave and matches the Gaussian beam at the interface of the particulate medium, and finally, derive the vector radiative transfer equation. For weakly focused Gaussian beams, the resulting equation is the traditional radiative transfer equation

Gaussian beam↗
DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Improving the electrostatic design of the Jefferson Lab 300 kV DC photogun

The 300 kV DC high voltage photogun at Jefferson Lab was redesigned to deliver electron beams with a much higher bunch charge and improved beam properties. The original design provided only a modest longitudinal electric field ( E z ) at the photocathode, which limited the achievable extracted bunch charge. To reach the bunch charge goal of approximately few nC with 75 ps full-width at half-maximum Gaussian laser pulse width, the existing DC high voltage photogun electrodes and anode–cathode gap were modified to increase E z at the photocathode. In addition, the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the non-symmetric nature of the inverted insulator photogun design. We present the electrostatic design of the original photogun and the modified photogun and beam dynamics simulations that predict vastly improved performance. We also quantify the impact of the photocathode recess on beam quality, where recess describes the actual location of the photocathode inside the photogun cathode electrode relative to the intended location. A photocathode unintentionally recessed/misplaced by sub-millimeter distance can significantly impact the downstream beam size.

47 OTHER INSTRUMENTATION↗
NASA NTRS1980

An experimental study of the surface elevation probability distribution and statistics of wind-generated waves

Laboratory experiments were performed to measure the surface elevation probability density function and associated statistical properties for a wind-generated wave field. The laboratory data along with some limited field data were compared. The statistical properties of the surface elevation were processed for comparison with the results derived from the Longuet-Higgins (1963) theory. It is found that, even for the highly non-Gaussian cases, the distribution function proposed by Longuet-Higgins still gives good approximations.

Huang, N. E.↗
NASA NTRS1988

Range data description based on multiple characteristics

An algorithm for describing range images based on Mean curvature (H) and Gaussian curvature (K) is presented. Range images are unique in that they directly approximate the physical surfaces of a real world 3-D scene. The curvature parameters are derived from the fundamental theorems of differential geometry and provides visible invariant pixel labels that can be used to characterize the scene. The sign of H and K can be used to classify each pixel into one of eight possible surface types. Due to the sensitivity of these parameters to noise the resulting HK-sing map does not directly identify surfaces in the range images and must be further processed. A region growing algorithm based on modeling the scene points with a Markov Random Field (MRF) of variable neighborhood size and edge models is suggested. This approach allows the integration of information from multiple characteristics in an efficient way. The performance of the proposed algorithm on a number of synthetic and real range images is discussed.

Al-Hujazi, Ezzet↗
NASA NTRS1985

Gaussian beams for surface waves in laterally slowly-varying media

Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. The elastodynamic equations of motion in ray-centered coordinates are derived, and a laterally slowly-varying approximation for elastodynamic equations is obtained. Parabolic equations for Love and Rayleigh waves are studied and solved, and the properties of Gaussian beams of seismic surface waves are examined.

Yomogida, K.↗
DOE OSTI2024

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗