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 469 records · Page 26

On the application of Rice's exceedance statistics to atmospheric turbulence.

Discrepancies produced by the application of Rice's exceedance statistics to atmospheric turbulence are examined. First- and second-order densities from several data sources have been measured for this purpose. Particular care was paid to each selection of turbulence that provides stationary mean and variance over the entire segment. Results show that even for a stationary segment of turbulence, the process is still highly non-Gaussian, in spite of a Gaussian appearance for its first-order distribution. Data also indicate strongly non-Gaussian second-order distributions. It is therefore concluded that even stationary atmospheric turbulence with a normal first-order distribution cannot be considered a Gaussian process, and consequently the application of Rice's exceedance statistics should be approached with caution.

Chen, W. Y.↗

Kurtosis, skewness, and non-Gaussian cosmological density perturbations

Cosmological topological defects as well as some nonstandard inflation models can give rise to non-Gaussian density perturbations. Skewness and kurtosis are the third and fourth moments that measure the deviation of a distribution from a Gaussian. Measurement of these moments for the cosmological density field and for the microwave background temperature anisotropy can provide a test of the Gaussian nature of the primordial fluctuation spectrum. In the case of the density field, the importance of measuring the kurtosis is stressed since it will be preserved through the weakly nonlinear gravitational evolution epoch. Current constraints on skewness and kurtosis of primeval perturbations are obtained from the observed density contrast on small scales and from recent COBE observations of temperature anisotropies on large scales. It is also shown how, in principle, future microwave anisotropy experiments might be able to reveal the initial skewness and kurtosis. It is shown that present data argue that if the initial spectrum is adiabatic, then it is probably Gaussian, but non-Gaussian isocurvature fluctuations are still allowed, and these are what topological defects provide.

Luo, Xiaochun↗

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytic Neural Network Gaussian Process Enabled Chance-Constrained Voltage Regulation for Active Distribution Systems with PVs, Batteries and EVs

This paper proposes an analytic neural network Gaussian process (NNGP)-based chance-constrained real-time voltage regulation method for active distribution systems with photovoltaics (PVs), batteries, and electric vehicles (EVs). NNGP can utilize historical measurement data to achieve real-time probabilistic node voltage estimation through Bayesian inference. Then, NNGP is fully analytically embedded into the optimal power flow model to perform voltage regulation and adapt to various topological changes. The uncertainties of voltage estimations are easily considered via the chance constraint, and it has been shown that the adoption of this chance constraint can significantly improve the reliability of voltage regulation under various scenarios. The comparison results with other methods, carried out on a real 759-node distribution system located in western Colorado, U.S., show that the proposed method can achieve accurate voltage estimation across different topologies and reliably perform voltage regulation considering PVs, batteries, and EVs.

active distribution systems↗

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗

Approximation to the Probability Density at the Output of a Photmultiplier Tube

The probability density of the integrated output of a photomultiplier tube (PMT) is approximated by the Gaussian, Rayleigh, and Gamma probability densities. The accuracy of the approximations depends on the signal energy alpha: the Gamma distribution is accurate for all alpha, the Raleigh distribution is accurate for small alpha (approximate or less than 1 photon) and the Gaussian distribution is accurate for large alpha (approximate or greater than 10 photons).

Stokey, R. J.↗

Skewness in large-scale structure and non-Gaussian initial conditions

We compute the skewness of the galaxy distribution arising from the nonlinear evolution of arbitrary non-Gaussian intial conditions to second order in perturbation theory including the effects of nonlinear biasing. The result contains a term identical to that for a Gaussian initial distribution plus terms which depend on the skewness and kurtosis of the initial conditions. The results are model dependent; we present calculations for several toy models. At late times, the leading contribution from the initial skewness decays away relative to the other terms and becomes increasingly unimportant, but the contribution from initial kurtosis, previously overlooked, has the same time dependence as the Gaussian terms. Observations of a linear dependence of the normalized skewness on the rms density fluctuation therefore do not necessarily rule out initially non-Gaussian models. We also show that with non-Gaussian initial conditions the first correction to linear theory for the mean square density fluctuation is larger than for Gaussian models.

Fry, J. N.↗

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. Here in this work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗

Efficiency degradation due to tracking errors for point focusing solar collectors

An important parameter in the design of point focusing solar collectors is the intercept factor which is a measure of efficiency and of energy available for use in the receiver. Using statistical methods, an expression of the expected value of the intercept factor is derived for various configurations and control law implementations. The analysis assumes that a radially symmetric flux distribution (not necessarily Gaussian) is generated at the focal plane due to the sun's finite image and various reflector errors. The time-varying tracking errors are assumed to be uniformly distributed within the threshold limits and allows the expected value calculation.

Hughes, R. O.↗

Statistics of density maxima and the large-scale matter distribution

High peaks in Gaussian noise display enhanced clustering. The enhancement takes two forms: on large scales one obtains a linear amplification of the correlation function which is independent of scale. On smaller scales, but larger than the mass scale of the peaks themselves, a nonlinear (exponential) enhancement of the number density of high peaks in overdense regions arises. The large-scale correlations of Abell's rich clusters can be understood as a manifestation of this phenomenon. If the formation of bright galaxies favors the high overdensity peaks then the number of galaxies (per unit mass) in clusters and groups may be considerably enhanced. Consequences of these ideas for the density parameter and the large-scale matter distribution are discussed.

Kaiser, N.↗

Emperical Tests of Acceptance Sampling Plans

Acceptance sampling is a quality control procedure applied as an alternative to 100% inspection. A random sample of items is drawn from a lot to determine the fraction of items which have a required quality characteristic. Both the number of items to be inspected and the criterion for determining conformance of the lot to the requirement are given by an appropriate sampling plan with specified risks of Type I and Type II sampling errors. In this paper, we present the results of empirical tests of the accuracy of selected sampling plans reported in the literature. These plans are for measureable quality characteristics which are known have either binomial, exponential, normal, gamma, Weibull, inverse Gaussian, or Poisson distributions. In the main, results support the accepted wisdom that variables acceptance plans are superior to attributes (binomial) acceptance plans, in the sense that these provide comparable protection against risks at reduced sampling cost. For the Gaussian and Weibull plans, however, there are ranges of the shape parameters for which the required sample sizes are in fact larger than the corresponding attributes plans, dramatically so for instances of large skew. Tests further confirm that the published inverse-Gaussian (IG) plan is flawed, as reported by White and Johnson (2011).

White, K. Preston, Jr.↗

On the Response of a Nonlinear Structure to High Kurtosis Non-Gaussian Random Loadings

This paper is a follow-on to recent work by the authors in which the response and high-cycle fatigue of a nonlinear structure subject to non-Gaussian loadings was found to vary markedly depending on the nature of the loading. There it was found that a non-Gaussian loading having a steady rate of short-duration, high-excursion peaks produced essentially the same response as would have been incurred by a Gaussian loading. In contrast, a non-Gaussian loading having the same kurtosis, but with bursts of high-excursion peaks was found to elicit a much greater response. This work is meant to answer the question of when consideration of a loading probability distribution other than Gaussian is important. The approach entailed nonlinear numerical simulation of a beam structure under Gaussian and non-Gaussian random excitations. Whether the structure responded in a Gaussian or non-Gaussian manner was determined by adherence to, or violations of, the Central Limit Theorem. Over a practical range of damping, it was found that the linear response to a non-Gaussian loading was Gaussian when the period of the system impulse response is much greater than the rate of peaks in the loading. Lower damping reduced the kurtosis, but only when the linear response was non-Gaussian. In the nonlinear regime, the response was found to be non-Gaussian for all loadings. The effect of a spring-hardening type of nonlinearity was found to limit extreme values and thereby lower the kurtosis relative to the linear response regime. In this case, lower damping gave rise to greater nonlinearity, resulting in lower kurtosis than a higher level of damping.

Rizzi, Stephen A.↗

Physics of debris clouds from hypervelocity impacts

The protection scheme developed for long duration space platforms relies primarily upon placing thin metal plates or 'bumpers' around flight critical components. The effectiveness of this system is highly dependent upon its ability to break up and redistribute the momentum of any particle which might otherwise strike the outer surface of the spacecraft. Therefore it is of critical importance to design the bumpers such that maximum dispersion of momentum is achieved. This report is devoted to an in-depth study into the design and development of a laboratory instrument which would permit the in-situ monitoring of the momentum distribution as the impact event occurs. A series of four designs were developed, constructed and tested culminating with the working instrument which is currently in use. Each design was individually tested using the Space Environmental Effects Facility (SEEF) at the Marshall Space Flight Center in Huntsville, Alabama. Along with the development of the device, an experimental procedure was developed to assist in the investigation of various bumper materials and designs at the SEEF. Preliminary results were used to compute data which otherwise were not experimentally obtainable. These results were shown to be in relative agreement with previously obtained values derived through other methods. The results of this investigation indicated that momentum distribution could in fact be measured in-situ as the impact event occurred thus giving a more accurate determination of the effects of experimental parameters on the momentum spread. Data produced by the instrument indicated a Gaussian-type momentum distribution. A second apparatus was developed and it was placed before the shield in the line of travel utilized a plate to collect impact debris scattered backwards. This plate had a passage hole in the center to allow the particle to travel through it and impact the proposed shield material. Applying the law of conservation of angular momentum a backward momentum vector was determined from the angular velocity of the plate. The forward scattered and backward scattered momentum values were then analyzed to judge the distribution of debris. Loss of momentum was attributed to the inaccuracies of the means of measurement. Assumptions of symmetrical debris for the forward and backward scattered directions also contributed to this loss.

Zee, Ralph↗

Spectral characteristics of the return

Doppler spectra of VHF radar returns typically indicate a Gaussian back-ground shape with superimposed spikes. An average of 10 Doppler spectra are shown which are calculated from a time series of seven min of complex data. A proper Gaussian fit to the background distribution is possible by neglecting the strong amplitude spikes. If this background distribution is due to beam width broadening, either diffuse reflection or rather isotropic scattering is required. If beam width broadening is neglected, the width of the distribution is given by turbulent velocity fluctuations. The reflected component is about 1/3 of the scattered component. It is assumed that the amplitude spikes due to diffuse reflection indicate a Gaussian frequency distribution such as for the amplitudes due to scattering.

Rottger, J.↗

Successful Cleanroom Installation of PIP-II SSR2 Coupler Using Robotic Arm

The Fermilab Side-Coupled Linac accelerates H- beam from 116 MeV to 400 MeV through seven 805 MHz modules. Twelve wire scanners are present in the Side Coupled Linac and four are present in the transfer line between the Linac and the Booster synchrotron ring. These wire scanners act as important diagnostic instruments to directly collect information on the beam s transverse distribution. The manipulation of the conditions of wire scanner data collection enables further characterization of the beamline, such as calculating emittance and the Twiss parameters of the beam at select regions. Here we present the results of these studies and characterization of the non-Gaussian transverse beam distribution observed.

Narug, C.↗

Emittance Measurements with Wire Scanners in the Fermilab Side-Coupled Linac

The Fermilab Side-Coupled Linac accelerates H- beam from 116 MeV to 400 MeV through seven 805 MHz modules. Twelve wire scanners are present in the Side Coupled Linac and four are present in the transfer line between the Linac and the Booster synchrotron ring. These wire scanners act as important diagnostic instru-ments to directly collect information on the beam s transverse distribution. The manipulation of the condi-tions of wire scanner data collection enables further characterization of the beamline, such as calculating emittance and the Twiss parameters of the beam at select regions. Here we present the results of these studies and characterization of the non-Gaussian transverse beam distribution observed.

43 PARTICLE ACCELERATORS↗