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 253 records · Page 14

Physics-informed machine learning modeling for predictive control using noisy data

Due to the occurrence of over-fitting at the learning phase, the modeling of chemical processes via artificial neural networks (ANN) by using corrupted data (i.e., noisy data) is an ongoing challenge. Therefore, this work investigates the effect of both Gaussian and non-Gaussian noise on the performance of process-structure based recurrent neural networks (RNN) models, which take the form of partially-connected RNN models in this work, that are used to approximate a class of multi-input-multi-outputs nonlinear systems. Furthermore, two different techniques, specifically Monte Carlo dropout and co-teaching, are utilized in the development of partially-connected RNN models. Here, these two techniques are employed to reduce the over-fitting in ANNs when noisy data is used in the training process and, hence, to improve the open-loop accuracy as well as the closed-loop performance under a Lyapunov-based model predictive controller (MPC). Aspen Plus Dynamics, a well-known high-fidelity process simulator, is used to simulate a large-scale chemical process application in order to demonstrate the anticipated improvements in both open-loop approximation and closed-loop controller performance in the presence of Gaussian and non-Gaussian noise in the data set using physics-informed RNNs.

97 MATHEMATICS AND COMPUTING↗

Statistical characteristics of laser interference and its effect on IR optoelectronic observation systems

We consider statistical characteristics of laser interference recorded in the focal region of an IR optoelectronic observation system under field and out-of-field illumination by laser radiation. The applicability of the analytical description of experimental histograms of the laser interference signal distribution at the photodetector array output using the probability density approximations in the form of gamma distribution and the Gaussian function is substantiated. It is shown that the specific form of these functions is determined both by the mean value (mathematical expectation) of the interference signal and by the characteristic parameter M, which depends on the statistical properties of the laser interference. (laser applications and other issues in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗

Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments

The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hubble space telescope observations of young star clusters in NGC-4038/4039, 'the antennae' galaxies

New, high-resolution images of the disks of NGC 4038/4039 obtained with the Wide Field Camera of the Hubble Space Telescope (HST) are presented. NGC 4038/4039, nicknamed The Antennae, is a prototypical example of a pair of colliding galaxies believed to be at an early stage of a merger. Down to the limiting magnitude of V approximately 23 mag, the HST images reveal a population of over 700 blue pointlike objects within the disks. The mean absolute magnitude of these objects is M(sub V) = -11 mag, with the brightest objects reaching M(sub V) approximately -15. Their mean apparent color indices ar U - V = -0.7 mag and V - 1 = 0.8 mag on the Johnson UVI passband system, while their mean indices corrected for internal reddening are (u - v)(sub 0) = -1.0 mag and (V - I(sub 0) = 0.5. Their mean effective radius, determined from slightly resolved images, is 18 pc (for H(sub 0) = 50 km/s /Mpc). Based on their luminosities and resolution, most of these objects cannot be individual stars, but are likely young compact star clusters. The brighter ones are similar to the objects found in NGC 1275 and NGC 7252, which appear to be young globular clusters formed during recent galazy mergers. Based on their U - V and V - I colors, the brightest, bluest clusters of NGC 4038/4039 appear to be less than 10 Myr old. Most of these bright clusters are relatively tightly clustered themselves, with typically a dozen individual clusters belonging to a complex identified as a giant H II region from ground-based observations. The cluster luminosity function (LF) is approximately a power law, phi(L)dL proportional to L(exp -1.78+/-0.05)dL, with no hint of a turnover at fainter magnitudes. This power-law shape agrees with the LF of Magellanic Cloud clusters and Galactic open clusters, but differs from the LF of old globular cluster systems that is typically Gaussian with a Full Width at Half Maximum (FWHM) of approximately 3 mag. Besides the blue clusters, we also find about a dozen extremely red objects with V - I greater than 3.0. The highest number density of these red objects is found in the SE quadrant, where star formation appears to be most recent. We propose that these objects may be very young star clusters still embedded in their placental dust cocoons.

Whitmore, Bradley C.↗

Solution of the stochastic control problem in unbounded domains.

Bellman's dynamic programming equation for the optimal index and control law for stochastic control problems is a parabolic or elliptic partial differential equation frequently defined in an unbounded domain. Existing methods of solution require bounded domain approximations, the application of singular perturbation techniques or Monte Carlo simulation procedures. In this paper, using the fact that Poisson impulse noise tends to a Gaussian process under certain limiting conditions, a method which achieves an arbitrarily good approximate solution to the stochastic control problem is given. The method uses the two iterative techniques of successive approximation and quasi-linearization and is inherently more efficient than existing methods of solution.

Robinson, P.↗

The motions of clusters and group of galaxies

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard cold dark matter (CDM) (omega = 1); low-density CDM (omega = 0.3); hot dark matter (HDM) (omega = 1); and primeval baryonic isocurvature (PBI) (omega = 0.3). All models are normalized to the microwave background fluctuations observed by Cosmic Background Explorer (COBE). We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models, as expected from a Gaussian initial density field. The clusters appear to be fundamental and efficient tracers of the large-scale velocity field. The cluster three-dimensional velocity distribution typically peaks at v approximately 600 km/s and extends to high cluster velocities of v approximately 2000 km/s. The low-density CDM model exhibits somewhat lower velocities: it peaks at approximately 400 km/s and extends to approximately 1200 km/s. Approximately 10% (approximately 1% for low-density CDM) of all model rich clusters move with high peculiar velocities of V greater than or = 10(exp 3) km/s. The highest velocity clusters frequently originate in dense superclusters. The model velocity distributions of rich clusters are compared with the model velocity distributions of small groups of galaxies, and of the total matter. The group velocity distribution is, in general, similar to the velocity distribution of the rich clusters. The matter velocity distribution is similar to that of the rich clusters for the omega = 0.3 models; these models exhibit Maxwellian velocity distributions for clusters, for groups, and for matter that are all similar to one another. The mass distribution in the omega = 1 models, however, exhibits a longer tail of high velocities than do the clusters. This high-velocity tail originates mostly from the high velocities that exist within rich clusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data exhibit a larger high-velocity tail, to v(sub r) greater than or = 2000 km/s, than is seen in the model simulations (except for HDM). Because of the large observational uncertainties, however, the data are consistent at a approximately 1-3 sigma level with the model predictions and with a Gaussian initial density field. Accurate observations of cluster peculiar velocities, especially at the high-velocity tail, shoud provide powerful constraints on the cosmological models.

Bahcall, Neta A.↗

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Detection of a sub-milliarcsecond radio component in the RS CVn system HR 1099

The RS CVn system HR 1099 was observed with a five station VLBI array at a frequency of 8.4 GHz during a strong radio outburst of approximately 400 mJy. The data are consistent with a circular Gaussian source of 0.8 + or - 0.12 milli-arcsec (FWHM), corresponding to a linear size of 4 + or - 0.6 x 10 to the 11th cm. This is comparable to the distance between the surfaces of the two stars, or to 75 percent of the diameter of the chromospherically active K star. Extrapolation of published photometric data shows that the starspot formation of HR 1099 was facing toward the observers at the times of observations. The high equivalent brightness temperature, approximately 10 billion K, is consistent with gyrosynchrotron emission from a power-law energy distribution of electrons in a magnetic field of strength about 30 gauss.

Mutel, R. L.↗

Note on the modified two-stream approximation of Sagan and Pollack.

The modified two-stream approximation is derived from the radiative transfer equation using a two-point gaussian quadrature and neglecting terms of order l equal to or greater than 3 in the Legendre expansion of the phase function. The results are compared with the Eddington approximation and with exact results for the special cases of perfect absorption, perfect scattering, and a semi-infinite layer with isotropic scattering.

Lyzenga, D. R.↗

Learning likelihood ratios with neural network classifiers

The likelihood ratio is a crucial quantity for statistical inference in science that enables hypothesis testing, construction of confidence intervals, reweighting of distributions, and more. Many modern scientific applications, however, make use of data- or simulation-driven models for which computing the likelihood ratio can be very difficult or even impossible. By applying the so-called “likelihood ratio trick,” approximations of the likelihood ratio may be computed using clever parametrizations of neural network-based classifiers. A number of different neural network setups can be defined to satisfy this procedure, each with varying performance in approximating the likelihood ratio when using finite training data. We present a series of empirical studies detailing the performance of several common loss functionals and parametrizations of the classifier output in approximating the likelihood ratio of two univariate and multivariate Gaussian distributions as well as simulated high-energy particle physics datasets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Radial density profile and stability of capillary discharge plasma waveguides of lengths up to 40 cm

We measured the parameter reproducibility and radial electron density profile of capillary discharge waveguides with diameters of 650 μm to 2 mm and lengths of 9 to 40 cm. To the best of the authors’ knowledge, 40 cm is the longest discharge capillary plasma waveguide to date. This length is important for ≥10 GeV electron energy gain in a single laser-driven plasma wakefield acceleration stage. Evaluation of waveguide parameter variations showed that their focusing strength was stable and reproducible to <0.2% and their average on-axis plasma electron density to <1%. These variations explain only a small fraction of laser-driven plasma wakefield acceleration electron bunch variations observed in experiments to date. Measurements of laser pulse centroid oscillations revealed that the radial channel profile rises faster than parabolic and is in excellent agreement with magnetohydrodynamic simulation results. We show that the effects of non-parabolic contributions on Gaussian pulse propagation were negligible when the pulse was approximately matched to the channel. However, they affected pulse propagation for a non-matched configuration in which the waveguide was used as a plasma telescope to change the focused laser pulse spot size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A spatially regularized detector for emergent/re-emergent disease outbreaks

Early detection of outbreaks caused by emergent pathogens, using epidemiological surveillance data i.e., daily case counts, is difficult. This is because the data tend to be noisy during the early epoch of the outbreak. In contrast, the spread-rate of the disease tends to be well-behaved, as it depends only on the mixing patterns of the population and the characteristics of the pathogen, neither of which behave erratically in space-time. In this report, we explore whether the spread-rate can be used for epidemiological surveillance, conditional on case count data. Estimating the spread-rate from case count data allows us to exploit exogenous information, e.g., incubation period distributions etc., which can considerably smooth out any erratic temporal behavior. Further, epidemiological dynamics are spatially correlated, and if case counts are available for multiple areal units e.g., counties, these correlations could potentially be used to suppress noise in the early epoch data. These exogenous information and structure are not exploited by conventional syndromic surveillance detectors to extract a well-behaved latent variable for monitoring purposes. The technical challenge lies in the estimation of the spread-rate field jointly over a collection of areal units; further, the spread-rate varies over time. We develop a method based on mean-field variational inference to approximately estimate the spread-rate field, using a Gaussian Random Field Model for spatial regularization. The method is tested on the estimation of spread-rate in the thirty-three counties of New Mexico and detect the arrival of the Fall 2020 COVID-19 wave in September 2020. We find that the method is scalable, but underestimates the uncertainty in the estimated spread-rate field. We detect the arrival of the Fall 2020 wave a week ahead of conventional syndromic surveillance algorithms, but our simplistic detection algorithm, based on simple anomaly detection, suffers from a high false positive rate, similar to conventional detectors.

59 BASIC BIOLOGICAL SCIENCES↗

Scattering of electromagnetic waves from a randomly perturbed quasiperiodic surface

Electromagnetic-wave scattering by a quasi-periodic surface with random perturbations (as in the remote sensing of plowed fields) is investigated analytically, applying the Kirchhoff approximation and modeling the plowed fields by means of Gaussian random variation, sinusoidal variation, and Gaussian random variation about the spatial frequency. Coherent and incoherent bistatic scattering coefficients are derived in closed form by evaluating the physical-optics integral and shown to be proportional, in the geometric-optics limit, to the occurrence probability of slopes which reflect the incident wave specularly in the direction of the scattered wave. Backscattering cross sections are plotted as functions of incidence angle for a number of cases, demonstrating the strong effect of row direction.

Shin, R. T.↗

Hard X-Ray Flare Source Sizes Measured with the Ramaty High Energy Solar Spectroscopic Imager

Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations of 18 double hard X-ray sources seen at energies above 25 keV are analyzed to determine the spatial extent of the most compact structures evident in each case. The following four image reconstruction algorithms were used: Clean, Pixon, and two routines using visibilities maximum entropy and forward fit (VFF). All have been adapted for this study to optimize their ability to provide reliable estimates of the sizes of the more compact sources. The source fluxes, sizes, and morphologies obtained with each method are cross-correlated and the similarities and disagreements are discussed. The full width at half-maximum (FWHM) of the major axes of the sources with assumed elliptical Gaussian shapes are generally well correlated between the four image reconstruction routines and vary between the RHESSI resolution limit of approximately 2" up to approximately 20" with most below 10". The FWHM of the minor axes are generally at or just above the RHESSI limit and hence should be considered as unresolved in most cases. The orientation angles of the elliptical sources are also well correlated. These results suggest that the elongated sources are generally aligned along a flare ribbon with the minor axis perpendicular to the ribbon. This is verified for the one flare in our list with coincident Transition Region and Coronal Explorer (TRACE) images. There is evidence for significant extra flux in many of the flares in addition to the two identified compact sources, thus rendering the VFF assumption of just two Gaussians inadequate. A more realistic approximation in many cases would be of two line sources with unresolved widths. Recommendations are given for optimizing the RHESSI imaging reconstruction process to ensure that the finest possible details of the source morphology become evident and that reliable estimates can be made of the source dimensions.

Dennis, Brian R.↗

Equivalent linearization for fatigue life estimates of a nonlinear structure

An adaptation of the method of equivalent linearization is developed in order to estimate the fatigue life of a nonlinear structure excited by Gaussian white noise. While conventional equivalent linearization has been successfully applied to approximate the mean square responses of the variety of complex nonlinear systems, it has not been possible to apply the method for fatigue predictions. This limitation results from the assumption that the probability density of the response of the nonlinear system is Gaussian, which is not generally true. In the present effort, a modification of equivalent linearization, developed by the author, is employed which has been shown in recent studies to provide a much more accurate description of the random response of a nonlinear structure than provided by conventional equivalent linearization. A method is presented for estimating fatigue life based on this approach and results are compared with those obtained using both numerical simulation and the classical method. The fatigue life is estimated for a nonlinear plate vibrating in a single resonant mode with random white noise excitation. Excellent agreement is found between all three methods. These results indicate that equivalent linearization may be extended to predict the fatigue life of more complex nonlinear systems than can be analyzed at present.

Miles, R. N.↗