Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Ridge Regression Signal Processing

The introduction of the Global Positioning System (GPS) into the National Airspace System (NAS) necessitates the development of Receiver Autonomous Integrity Monitoring (RAIM) techniques. In order to guarantee a certain level of integrity, a thorough understanding of modern estimation techniques applied to navigational problems is required. The extended Kalman filter (EKF) is derived and analyzed under poor geometry conditions. It was found that the performance of the EKF is difficult to predict, since the EKF is designed for a Gaussian environment. A novel approach is implemented which incorporates ridge regression to explain the behavior of an EKF in the presence of dynamics under poor geometry conditions. The basic principles of ridge regression theory are presented, followed by the derivation of a linearized recursive ridge estimator. Computer simulations are performed to confirm the underlying theory and to provide a comparative analysis of the EKF and the recursive ridge estimator.

Kuhl, Mark R.↗
NASA NTRS1998

Pursuit Latency for Chromatic Targets

The temporal dynamics of eye movement response to a change in direction of stimulus motion has been used to compare the processing speeds of different types of stimuli (Mulligan, ARVO '97). In this study, the pursuit response to colored targets was measured to test the hypothesis that the slow response of the chromatic system (as measured using traditional temporal sensitivity measures such as contrast sensitivity) results in increased eye movement latencies. Subjects viewed a small (0.4 deg) Gaussian spot which moved downward at a speed of 6.6 deg/sec. At a variable time during the trajectory, the dot's direction of motion changed by 30 degrees, either to the right or left. Subjects were instructed to pursue the spot. Eye movements were measured using a video ophthalmoscope with an angular resolution of approximately 1 arc min and a temporal sampling rate of 60 Hz. Stimuli were modulated in chrominance for a variety of hue directions, combined with a range of small luminance increments and decrements, to insure that some of the stimuli fell in the subjects' equiluminance planes. The smooth portions of the resulting eye movement traces were fit by convolving the stimulus velocity with an exponential having variable onset latency, time constant and amplitude. Smooth eye movements with few saccades were observed for all stimuli. Pursuit responses to stimuli having a significant luminance component are well-fit by exponentials having latencies and time constants on the order of 100 msec. Increases in pursuit response latency on the order of 100-200 msec are observed in response to certain stimuli, which occur in pairs of complementary hues, corresponding to the intersection of the stimulus section with the subjects' equiluminant plane. Smooth eye movements can be made in response to purely chromatic stimuli, but are slower than responses to stimuli with a luminance component.

Mulligan, Jeffrey B.↗
NASA NTRS1999

Validity of the Relation Between Spontaneous and Stimulated Emissions in Semiconductors

The Einstein relation between spontaneous emission and absorption was originally derived for a system consists of a two-state subsystem representing matter and harmonic fields representing radiation. The derivation is based on the detailed balance between these two subsystems under thermal equilibrium. The relationship was later investigated in connection with the interactions between radiation field and solids or semiconductors. The simple derivation dose not hold for semiconductors in general. In certain limiting cases, simple relation was obtained. The validity of this relation is important not only because of its fundamental role connecting two of the most fundamental optical processes in semiconductors, but mostly also because of its wide use as a practical method to measure the optical gain of a semiconductor. The validity of this relation for semiconductors has been an issue of controversial for some time. In this paper we numerically examine the validity of this relationship for several different lineshapes including Lorentzian, Gaussian, Sech, and a convoluted double Lorentzians (CDL). We find out that at relatively low density above transparency level, all first three lineshapes violate the Einstein relation. The relation is approximately valid at high density. At very high density, the validity of the Einstein relation holds well for all three lineshapes. The reason behind this observation is explained. The CDL lineshape has been shown analytically to obey the Einstein relationship previously. We show that for a 2D semiconductor with parabolic bands, the CDL lineshape can be integrated analytically. This analytic lineshape is compared with a simple Lorentzian lineshape.

Ning, Cun-Zheng↗
DOE OSTI2016

Conservative tightly-coupled simulations of stochastic multiscale systems

Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domain-decomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domain-decomposition algorithm in which tight coupling is achieved by employing either Picard's or Newton's iterative method. Coupled diffusion equations, one of which has a Gaussian white-noise source term, provide a computational testbed for analysis of these two coupling strategies. Fully-converged (“implicit”) coupling with Newton's method typically outperforms its Picard counterpart, especially at high noise levels. This is because the number of Newton iterations scales linearly with the amplitude of the Gaussian noise, while the number of Picard iterations can scale superlinearly. At large time intervals between two subsequent inter-solver communications, the solution error for single-iteration (“explicit”) Picard's coupling can be several orders of magnitude higher than that for implicit coupling. Increasing the explicit coupling's communication frequency reduces this difference, but the resulting increase in computational cost can make it less efficient than implicit coupling at similar levels of solution error, depending on the communication frequency of the latter and the noise strength. This trend carries over into higher dimensions, although at high noise strength explicit coupling may be the only computationally viable option.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Uncertainty quantification of microstructural properties due to variability in measured pole figures

Experimental pole figures are an important input for microstructure homogenization models. In this paper, we derive an exact analytical formulation to quantify the uncertainties in homogenized properties due to uncertainty in the experimentally measured pole figures. The pole figures are acquired from a set of Ti-7Al alloy samples. These samples were obtained from the same process: by compressing a beta forged ingot at room temperature followed by annealing. The samples were taken from different regions of the original ingot, and this created variability in the resulting pole figures. The joint multivariate probability distributions of the computed orientation distribution function (ODF) is then found using the method of characteristic functions based on a Gaussian model of the pole figures. Engineering properties such as elastic modulus can be obtained by volume averaging over the ODF. We also show that uncertainty in elastic properties can be analytically obtained using direct transformation of random variables in the homogenization approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Apodization Specific Fitting for Improved Resolution, Charge Measurement, and Data Analysis Speed in Charge Detection Mass Spectrometry

Short-time Fourier transforms with short segment lengths are typically used to analyze single ion charge detection mass spectrometry (CDMS) data either to overcome effects of frequency shifts that may occur during the trapping period or to more precisely determine the time at which an ion changes mass or charge, or enters an unstable orbit. The short segment lengths can lead to scalloping loss unless a large number of zero-fills are used, making computational time a significant factor in real-time analysis of data. Apodization specific fitting leads to a 9-fold reduction in computation time compared to zero-filling to a similar extent of accuracy. This makes possible real-time data analysis using a standard desktop computer. Rectangular apodization leads to higher resolution than the more commonly used Gaussian or Hann apodization and makes it possible to separate ions with similar frequencies, a significant advantage for experiments in which the masses of many individual ions are measured simultaneously. Equally important is a >20% increase in S/N obtained with rectangular apodization compared to Gaussian or Hann, which directly translates to a corresponding improvement in accuracy of both charge measurements and ion energy measurements that rely on the amplitudes of the fundamental and harmonic frequencies. Finally, combined with computing the fast Fourier transform in a lower-level language, this fitting procedure eliminates computational barriers and should enable real-time processing of CDMS data on a laptop computer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

The Bogolyubov-Born-Green-Kirkwood-Yvon hierarchy and Fokker-Planck equation for many-body dissipative randomly driven systems

By generalizing Bogolyubov’s reduced description method, we suggest a formalism to derive kinetic equations for many-body dissipative systems in external stochastic field. As a starting point, we use a stochastic Liouville equation obtained from Hamilton’s equations taking dissipation and stochastic perturbations into account. The Liouville equation is then averaged over realizations of the stochastic field by an extension of the Furutsu-Novikov formula to the case of a non-Gaussian field. As the result, a generalization of the classical Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy is derived. In order to get a kinetic equation for the single-particle distribution function, we use a regular cutoff procedure of the BBGKY hierarchy by assuming weak interaction between the particles and weak intensity of the field. Within this approximation, we get the corresponding Fokker-Planck equation for the system in a non-Gaussian stochastic field. Two particular cases are discussed by assuming either Gaussian statistics of external perturbation or homogeneity of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Observations of CO in the Magellanic irregular galaxy NGC 55

The content of molecular gas in galaxies, mainly H2, is one of the key observations necessary for the understanding of star formation processes and history. As the CO molecule is the most widely distributed molecule after H2 and has easily observable mm lines, it is used as a tracer for the molecular gas. CO was detected towards the direction where the H alpha and 6 cm radio continuum emission is strongest (Hummel et al. 1986). Here, researchers present the Gaussian line parameters in tabular form. The distribution of CO corresponds well with the intense HI cloud near the bar of NGC 55. The extent of the CO cloud is about 975 pc perpendicular to the major axis. As the radio continuum and the H alpha emission also peaks in this region, it is most probably associated with the star forming region in NGC 55. Assuming that the molecular gas is in virial equilibrium, researchers derive a mass of about 8 times 10(exp 7) solar magnitude. The molecular mass found indicates that the conversion factor for the molecular mass in Irr galaxies as inferred from CO line emission is indeed higher by up to a factor of 20 compared to the canonical value for the Galaxy.

Heithausen, Andreas↗
DOE OSTI2016

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗
DOE OSTI2015

Theoretical detection threshold of the proton-acoustic range verification technique

Purpose: Range verification in proton therapy using the proton-acoustic signal induced in the Bragg peak was investigated for typical clinical scenarios. The signal generation and detection processes were simulated in order to determine the signal-to-noise limits. Methods: An analytical model was used to calculate the dose distribution and local pressure rise (per proton) for beams of different energy (100 and 160 MeV) and spot widths (1, 5, and 10 mm) in a water phantom. In this method, the acoustic waves propagating from the Bragg peak were generated by the general 3D pressure wave equation implemented using a finite element method. Various beam pulse widths (0.1–10 μs) were simulated by convolving the acoustic waves with Gaussian kernels. A realistic PZT ultrasound transducer (5 cm diameter) was simulated with a Butterworth bandpass filter with consideration of random noise based on a model of thermal noise in the transducer. The signal-to-noise ratio on a per-proton basis was calculated, determining the minimum number of protons required to generate a detectable pulse. The maximum spatial resolution of the proton-acoustic imaging modality was also estimated from the signal spectrum. Results: The calculated noise in the transducer was 12–28 mPa, depending on the transducer central frequency (70–380 kHz). The minimum number of protons detectable by the technique was on the order of 3–30 × 10{sup 6} per pulse, with 30–800 mGy dose per pulse at the Bragg peak. Wider pulses produced signal with lower acoustic frequencies, with 10 μs pulses producing signals with frequency less than 100 kHz. Conclusions: The proton-acoustic process was simulated using a realistic model and the minimal detection limit was established for proton-acoustic range validation. These limits correspond to a best case scenario with a single large detector with no losses and detector thermal noise as the sensitivity limiting factor. Our study indicated practical proton-acoustic range verification may be feasible with approximately 5 × 10{sup 6} protons/pulse and beam current.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Robust negativity in the quantum-to-classical transition of Kerr dynamics

Here, we quantify the quantum-to-classical transition of the single-mode Kerr nonlinear dynamics in the presence of loss. We establish three timescales that govern the dynamics, each with distinct characteristics. For times short compared with the Ehrenfest time, the evolution is classical, characterized by Gaussian dynamics. For sufficiently long times, as we increase the initial photon number, unitary Kerr evolution would generate macroscopic superpositions of coherent states (so-called kitten states). However, this is severely restricted in the presence of small photon loss, and the expectation values of observables coincide with their classical values. The intermediate timescale, however, shows resilient quantum behavior in the macroscopic limit. We show that in the mean-field non-Gaussian regime, the Kerr Hamiltonian (with small photon loss) generates a significant amount of Wigner-negativity, and classical flow is recovered only if the loss rate grows with system size. Our results broaden the usual understanding of quantum-to-classical transitions and demonstrate the potential for creating robust nonclassical resources for continuous-variable quantum information processing in the presence of loss.

Raza, Mohsin [University of New Mexico, Albuquerqu↗
DOE OSTI2017

HALOGAS Observations of NGC 4559: Anomalous and Extraplanar H i and its Relation to Star Formation

We use new deep 21 cm H i observations of the moderately inclined galaxy NGC 4559 in the HALOGAS survey to investigate the properties of extraplanar gas. We use TiRiFiC to construct simulated data cubes to match the H i observations. We find that a thick-disk component of scale height ∼2 kpc, characterized by a negative vertical gradient in its rotation velocity (lag) of ∼13 ± 5 km s{sup −1} kpc{sup −1} is an adequate fit to extraplanar gas features. The tilted ring models also present evidence for a decrease in the magnitude of the lag outside R {sub 25}, and a radial inflow of ∼10 km s{sup −1}. We extracted lagging extraplanar gas through Gaussian velocity profile fitting. From both the 3D models and extraction analyses we conclude that ∼10%–20% of the total H i mass is extraplanar. Most of the extraplanar gas is spatially coincident with regions of star formation in spiral arms, as traced by H α and GALEX FUV images, so it is likely due to star formation processes driving a galactic fountain. We also find the signature of a filament of a kinematically “forbidden” H i feature, containing ∼1.4 × 10{sup 6} M {sub ⊙} of H i, and discuss its potential relationship to a nearby H i hole. We discover a previously undetected dwarf galaxy in H i located ∼0.°4 (∼58 kpc) from the center of NGC 4559, containing ∼4 × 10{sup 5} M {sub ⊙}. This dwarf has counterpart sources in SDSS with spectra typical of H ii regions, and we conclude that it is two merging blue compact dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Statistics of multi-look AIRSAR imagery: A comparison of theory with measurements

The intensity and amplitude statistics of SAR images, such as L-Band HH for SEASAT and SIR-B, and C-Band VV for ERS-1 have been extensively investigated for various terrain, ground cover and ocean surfaces. Less well-known are the statistics between multiple channels of polarimetric of interferometric SAR's, especially for the multi-look processed data. In this paper, we investigate the probability density functions (PDF's) of phase differences, the magnitude of complex products and the amplitude ratios, between polarization channels (i.e. HH, HV, and VV) using 1-look and 4-look AIRSAR polarimetric data. Measured histograms are compared with theoretical PDF's which were recently derived based on a complex Gaussian model.

Lee, J. S.↗
NASA NTRS2000

The High Mosaicity Illusion: Revealing the True Physical Characteristics of Macromolecular Crystals

An experimental system and software have been developed for simultaneously measuring the diffraction resolution and mosaic spread of macromolecular crystals. Hundreds of reflection profiles over a wide resolution range were rapidly measured by using a charge coupled device (CCD) area detector in combination with superfine phi slicing data collection. The contributions of the X-ray beam to the reflection widths were minimized by using a highly-parallel, highly-monochromatic synchrotron source. These contributions and Lorentz effects were evaluated and deconvoluted from the recorded data. Data collection and processing is described. From one degree of superfine phi slice data collected on a crystal of manganese superoxide dismutase the mosaicity of 261 reflections were measured. The average mosaicity was 0.0101 degrees (0.0035) at the full-width-at-half-maximum (FWHM) and ranged from 0.0011 degrees to 0.0188 degrees. Each reflection profile was individually fit with two gaussian profiles with the first gaussian contributing 55% and the second contributing 35% of the reflection. On average, the mosaicity of the first gaussian was 0.0054 degrees (0.0015) and the second was 0.0061 degrees (0.0023). The mosaicity of the crystal was anisotropic with fh, f k, and fl values of 0.0068 degrees, 0.0140 degrees and 0.0046 degrees, respectively at the FWHM. The anisotropic mosaicity analysis indicates that the crystal is the most perfect in the I direction which corresponds to the favored growth direction of the crystal.

Bellamy, Henry↗
NASA NTRS2022

Modeling the East-West Asymmetry of Energetic Particle Fluence in Large Solar Energetic Particle Events Using the iPATH Model

It has been noted that in large solar energetic particle (SEP) events, the peak intensities show an East-West asymmetry with respect to the source flare locations. Using the 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we investigate the origin of this longitudinal trend. We consider multiple cases with different solar wind speeds and eruption speeds of the coronal mass ejections (CMEs) and fit the longitudinal distributions of time-averaged fluence by symmetric/asymmetric Gaussian functions with three time intervals of 8, 24 and 48 hr after the flare onset time respectively. The simulation results are compared with a statistical study of three-spacecraft events. We suggest that the East-West asymmetry of SEP fluence and peak intensity can be primarily caused the combined effect of an extended shock acceleration process and the evolution of magnetic field connection to the shock front. Our simulations show that the solar wind speed and the CME speed are important factors determining the East-West fluence asymmetry.

Zheyi Ding↗
NASA NTRS1976

Digital processing of satellite imagery application to jungle areas of Peru

The author has identified the following significant results. The use of clustering methods permits the development of relatively fast classification algorithms that could be implemented in an inexpensive computer system with limited amount of memory. Analysis of CCTs using these techniques can provide a great deal of detail permitting the use of the maximum resolution of LANDSAT imagery. Potential cases were detected in which the use of other techniques for classification using a Gaussian approximation for the distribution functions can be used with advantage. For jungle areas, channels 5 and 7 can provide enough information to delineate drainage patterns, swamp and wet areas, and make a reasonable broad classification of forest types.

Pomalaza, J. C.↗
DOE OSTI2021

Estimating basis functions in massive fields under the spatial mixed effects model

Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.

Pazdernik, Karl↗
DOE OSTI2024

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗