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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Data-Driven Global Sensitivity Analysis of Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution system with respect to stochastic PV variations plays an important role in designing optimal voltage control schemes. This paper proposes a Kriging, i.e., Gaussian process modeling enabled data-driven GSA method. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol index can be conveniently calculated to assess the global sensitivity of voltage to various power injection variations. Comparison results with other model-based GSA methods on the IEEE 37-bus feeder, such as the polynomial chaos expansion and the Monte Carlo approaches demonstrate that the proposed method can achieve accurate GSA outcomes while maintaining high computational efficiency.

distribution system analysis↗
NASA NTRS2004

Characteristics of Electromagnetic Pulse Propagation in Metal

It is well known that the solution of the diffusion equation for an electromagnetic field with a time harmonic term, e(sup iwt), is in the form of a traveling wave whose amplitude attenuates over distance into a conducting medium. As the attenuation is an increasing function of frequency, the high frequency components attenuate more rapidly than those of low ones upon entering a well conducting object. At the same time, the phase velocity of an individual component is also an increasing function of frequency causing a broadening of the pulse traveling inside a conductor. In the results of our previous study of numerical simulations, the problem of using a gaussian input pulse was immediately clear. First, having the dominant frequency components distributed around zero, the movement of the peak was not well defined. Second, with the amplitude of fourier components varying slowly over a wide range, the dispersion-induced blurring of the peak position was seen to be severe. For the present study, we have used a gaussian modulated single frequency sinusoidal wave, i. e., the carrier, as an input pulse in an effort to improve the issues related to the unclear movement of peak and dispersion as described above. This was based on the following two anticipated advantages: First, the packet moves in a conductor at the group velocity calculated at the carrier frequency, which means it is well controllable. Second, the amplitude of frequency components other than that of the carrier can be almost negligible, such that the effect of dispersion can be significantly reduced. A series of experiments of transmitting electromagnetic pulses through aluminum plates of various thickness was performed to test the validity of the above points. The results of numerical simulation based on wave propagation are discussed with respect to the experimental results. Finally, a simple simulation was performed based on diffusion of a continuous sine wave input and the results are compared with those of a single frequency sinusoidal wave observed over time at difference locations inside a conductor.

Namkung, M.↗
DOE OSTI2022

Quantifying radial migration in the Milky Way: inefficient over short time-scales but essential to the very outer disc beyond ∼15 kpc

ABSTRACT Stellar radial migration plays an important role in reshaping a galaxy’s structure and the radial distribution of stellar population properties. In this work, we revisit reported observational evidence for radial migration and quantify its strength using the age–[Fe/H] distribution of stars across the Milky Way with APOGEE data. We find a broken age–[Fe/H] relation in the Galactic disc at r > 6 kpc, with a more pronounced break at larger radii. To quantify the strength of radial migration, we assume stars born at each radius have a unique age and metallicity, and then decompose the metallicity distribution function (MDF) of mono-age young populations into different Gaussian components that originated from various birth radii at rbirth < 13 kpc. We find that, at ages of 2 and 3 Gyr, roughly half the stars were formed within 1 kpc of their present radius, and very few stars (<5 per cent) were formed more than 4 kpc away from their present radius. These results suggest limited short-distance radial migration and inefficient long-distance migration in the Milky Way during the last 3 Gyr. In the very outer disc beyond 15 kpc, the observed age–[Fe/H] distribution is consistent with the prediction of pure radial migration from smaller radii, suggesting a migration origin of the very outer disc. We also estimate intrinsic metallicity gradients at ages of 2 and 3 Gyr of −0.061 and −0.063 dex kpc−1, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Heavy-ion fusion cross section formula and barrier height distribution

Methods for obtaining fusion cross section formulas are discussed, especially for those that take on an empirical form. A new expression starting with σ( E) = $\frac{1}{E}$ $\int ^{E}_{E_0}$ $(\int^{E'}_{E_0} B (E^\shortparallel) dE^\shortparallel)dE'$ has been explored. Here, B(E) is a reasonable, assumed function of the barrier height distribution, $\frac{d^2 (\sigma E)}{dE^2}$. Further, the resulting analytic cross section formula reproduces very well the excitation functions for many light and heavy fusion systems across wide energy ranges, when B(E) is assumed to be a multi-Gaussian function. This study offers an improved determination of the fusion barrier height distribution over other numerical techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗
DOE OSTI2020

A discrete integral transform for rapid spectral synthesis

Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed. A Python implementation of the method is compared against state-of-the-art spectral code RADIS, and is since adopted as RADIS's default synthesis method. The synthesis of a benchmark CO2 spectrum consisting of 1.8 M spectral lines and 200k spectral points took only 3.1 s using the proposed method (1011 lines × spectral points/s), a factor ~300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. Finally, an experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2~3 orders performance increase, achieving up to 5 ∙ 10 14 lines × spectral points/s.

42 ENGINEERING↗
DOE OSTI2022

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Nonlinear models for coupling the effects of stimulated Raman scattering to inertial confinement fusion codes

Laser plasma instabilities (LPI) reduce driver-target coupling, alter implosion symmetry, and therefore can fundamentally limit fusion performance in inertial confinement fusion (ICF). Developing a predictive modeling capability for LPI effects can critically advance the success of the field. We perform vector particle-in-cell simulations of multi-speckled laser beams undergoing stimulated Raman scattering (SRS) at various densities and intensities relevant to mainly indirectly driven and a subset of parameter space for directly driven ICF systems, focusing on the regimes with intensities above the onset of electron trapping. Based on the wavenumber of the SRS daughter electron plasma wave, we identify several regions with underpinning SRS saturation physics: the electron-trapping dominated region with intermediate kλD values, the strong Landau damping region at larger kλD values, and the region where the Langmuir decay instability arises at lower kλD values. We develop a nonlinear SRS reflectivity model that features the base trapping-dominated scaling of (kλD)−4 and its modifications. Electron trapping deforms the initialized electron distribution functions, and we have developed a new δf-Gaussian-mixture algorithm for an accurate characterization of the trapped hot electron population. With this SRS hot electron description, we construct a nonlinear energy deposition model and a hot electron source model—based on a modified Manley–Rowe relation—suitable for including SRS effects as a sub-grid module in a high-fidelity ICF design code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spot size measurement of a deuterium–tritium dense plasma focus using neutron radiography

Neutron radiography is a technique uniquely suited to applications in nuclear diagnostics, non-destructive testing, and subcritical experiments. The spatial resolution of neutron radiographs is degraded by optical blur in the imaging system and the neutron source size, where the ideal source is point-like to optimize the point-spread function. A potential neutron source for radiography is the dense plasma focus (DPF), a coaxial Z-pinch that produces thermonuclear and beam-target neutrons. To assess if the source size is suitable for radiography, a neutron imaging system was used to measure the source size of the 4 MA Sodium DPF at the Nevada National Security Site operating with deuterium–tritium gas-fill. The source size was measured using the edge-spread function of tungsten objects, each having a rolled (convex) edge. The spot size was found to be 7–12 mm full-width at half-max (FWHM) assuming a Gaussian source, though comparison is presented for Lorentzian and Bennett distributions. The average FWHM was found to be 8.6 ± 1.2 mm vertically and 10.8 ± 1.2 mm horizontally with respect to the image plane, averaging over varied edges and alignments. The results were sensitive to source alignment and edge metrology, which introduced notable uncertainties. These results are consistent with separate experimental measurements as well as magnetohydrodynamics simulations of this DPF, which suggest that neutron production can originate from pinches ∼5–7 mm off-axis. These results suggest that the DPF should be used for radiography at low magnification (M < 1) where spot size does not dominate spatial blur.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

QSO photometric redshifts from SDSS, WISE, and GALEX colours

ABSTRACT Machine learning techniques, specifically the k-nearest neighbour algorithm applied to optical band colours, have had some success in predicting photometric redshifts of quasi-stellar objects (QSOs): Although the mean of differences between the spectroscopic and photometric redshifts, Δ$z$, is close to zero, the distribution of these differences remains wide and distinctly non-Gaussian. As per our previous empirical estimate of photometric redshifts, we find that the predictions can be significantly improved by adding colours from other wavebands, namely the near-infrared and ultraviolet. Self-testing this, by using half of the 33 643 strong QSO sample to train the algorithm, results in a significantly narrower spread in Δ$z$ for the remaining half of the sample. Using the whole QSO sample to train the algorithm, the same set of magnitudes return a similar spread in Δ$z$ for a sample of radio sources (quasars). Although the matching coincidence is relatively low (739 of the 3663 sources having photometry in the relevant bands), this is still significantly larger than from the empirical method (2 per cent) and thus may provide a method with which to obtain redshifts for the vast number of continuum radio sources expected to be detected with the next generation of large radio telescopes.

Curran, S. J.↗
DOE OSTI2022

Fiber Uncertainty Visualization for Bivariate Data With Parametric and Nonparametric Noise Models

Visualization and analysis of multivariate data and their uncertainty are top research challenges in data visualization. Constructing fiber surfaces is a popular technique for multivariate data visualization that generalizes the idea of level-set visualization for univariate data to multivariate data. Here, in this paper, we present a statistical framework to quantify positional probabilities of fibers extracted from uncertain bivariate fields. Specifically, we extend the state-of-the-art Gaussian models of uncertainty for bivariate data to other parametric distributions (e.g., uniform and Epanechnikov) and more general nonparametric probability distributions (e.g., histograms and kernel density estimation) and derive corresponding spatial probabilities of fibers. In our proposed framework, we leverage Green's theorem for closed-form computation of fiber probabilities when bivariate data are assumed to have independent parametric and nonparametric noise. Additionally, we present a nonparametric approach combined with numerical integration to study the positional probability of fibers when bivariate data are assumed to have correlated noise. For uncertainty analysis, we visualize the derived probability volumes for fibers via volume rendering and extracting level sets based on probability thresholds. We present the utility of our proposed techniques via experiments on synthetic and simulation datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Measures of illuminant-induced metameric mismatch: theory, comparative analysis, and implications for application

Illuminant-induced metameric mismatch is an important consideration in the specification of light sources for some architectural environments, yet there is currently no standardized performance measure. The goal of this work was to evaluate two recent research proposals: the metameric uncertainty index ( R t ) and the metamer mismatching color rendering index (MMCRI). To compare the relative performance of these two measures, 100,000 spectral power distributions were generated with 3, 4, 5, 6, and 7 Gaussian spectral components and spectral widths varying from 1 nm (monochromatic) to 100 nm. Both measures generally agree with the theory that broadband radiation should cause less metameric mismatch than narrowband radiation. The two measures have relatively better agreement for broadband SPDs and relatively worse agreement for narrower spectra. Despite some similarities, non-parametric statistical tests suggest that R t and MMCRI are significantly different quantifications of illuminant-induced metameric mismatch ( p < 0.0001 for all comparisons). Characteristics of the MMCRI computation that are potentially problematic for applied lighting were observed.

42 ENGINEERING↗
DOE OSTI2020

We Know Where We Don't Know: 3D Bayesian CNNs for Credible Geometric Uncertainty

Deep learning has been successfully applied to the segmentation of 3D Computed Tomography (CT) scans. Establishing the credibility of these segmentations requires uncertainty quantification (UQ) to identify untrustworthy predictions. Recent UQ architectures include Monte Carlo dropout networks (MCDNs), which approximate deep Gaussian processes, and Bayesian neural networks (BNNs), which learn the distribution of the weight space. BNNs are advantageous over MCDNs for UQ but are thought to be computationally infeasible in high dimension, and neither architecture has produced interpretable geometric uncertainty maps. We propose a novel 3D Bayesian convolutional neural network (BCNN), the first deep learning method which generates statistically credible geometric uncertainty maps and scales for application to 3D data. We present experimental results on CT scans of graphite electrodes and laser-welded metals and show that our BCNN outperforms an MCDN in recent uncertainty metrics. The geometric uncertainty maps generated by our BCNN capture distributions of sigmoid values that are interpretable as confidence intervals, critical for applications that rely on deep learning for high-consequence decisions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A New Probability Table Methodology with Covariances on Probability Tables

Through this work, we test a new methodology of evaluating probability tables informed by a combination of physics models and experimental data. This methodology uses generalized least squares to produce a joint mean and covariance for average cross sections and probability tables. This methodology has been tested and validated using total Monte Carlo. Attempts were made to parameterize the probability table using a formulation provided in. Although the normal inverse Gaussian fit was successful for the total cross section probability distribution, the rational-function fit to the capture, inelastic, and fission cross sections was unsuccessful due to time constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

A comprehensive model to determine the effects of temperature and species fluctuations on reactions in turbulent reacting flows

The joint probability distribution function (pdf), which is a modification of the bivariate Gaussian pdf, is discussed and results are presented for a global reaction model using the joint pdf. An alternative joint pdf is discussed. A criterion which permits the selection of temperature pdf's in different regions of turbulent, reacting flow fields is developed. Two principal approaches to the determination of reaction rates in computer programs containing detailed chemical kinetics are outlined. These models represent a practical solution to the modeling of species reaction rates in turbulent, reacting flows.

Antaki, P. J.↗
NASA NTRS1989

A finite element formulation for the large deflection random response of thermally buckled beams

The effects of temperature and acoustic loading are included in a theoretical finite element large deflection formulation for thin, isotropic beams. Thermal loads are applied as steady-state temperature distributions, and acoustic loads are taken to be ergodic and Gaussian with zero mean and uniform magnitude and phase along the length of the beam. Material properties are considered presently to be independent of temperature. Also, inplane and rotary inertia terms are assumed to be negligible, and all inplane edge conditions are taken to be immovable. For the random response analysis, both auto- and cross-correlation terms are included. The nature of the loads leads to the solution of two separate problems. First, the problem of thermal postbuckling is solved to determine the deflections and stresses due to the thermal load only. These deflections and stresses are then used as initial deflections and stresses for the random vibration analysis. Root-mean-square (RMS) maximum deflections and strains are obtained and compared with previous classical equivalent linearization results.

Locke, James↗
NASA NTRS1990

Distributed memory approaches for robotic neural controllers

The suitability is explored of two varieties of distributed memory neutral networks as trainable controllers for a simulated robotics task. The task requires that two cameras observe an arbitrary target point in space. Coordinates of the target on the camera image planes are passed to a neural controller which must learn to solve the inverse kinematics of a manipulator with one revolute and two prismatic joints. Two new network designs are evaluated. The first, radial basis sparse distributed memory (RBSDM), approximates functional mappings as sums of multivariate gaussians centered around previously learned patterns. The second network types involved variations of Adaptive Vector Quantizers or Self Organizing Maps. In these networks, random N dimensional points are given local connectivities. They are then exposed to training patterns and readjust their locations based on a nearest neighbor rule. Both approaches are tested based on their ability to interpolate manipulator joint coordinates for simulated arm movement while simultaneously performing stereo fusion of the camera data. Comparisons are made with classical k-nearest neighbor pattern recognition techniques.

Jorgensen, Charles C.↗
NASA NTRS1993

Disk plus halo models of gamma-ray burst sources

The present consideration of Gaussian halo, exponential halo, and standard 'dark matter' halo forms for the distribution of Galactic halo sources indicates that about 70 percent of the gamma-ray bursts observed by the BATSE program can come from local Galactic disk neutron stars, with the rest lying in an extended Galactic halo. Two specific disk-halo combinations are considered as examples which use an exponential disk; in the first case, this is combined with a Gaussian halo, and in the second, with a dark matter halo.

Smith, I. A.↗