Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GAUSSIAN PROCESSES”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Image gathering and digital restoration for fidelity and visual quality

The fidelity and resolution of the traditional Wiener restorations given in the prevalent digital processing literature can be significantly improved when the transformations between the continuous and discrete representations in image gathering and display are accounted for. However, the visual quality of these improved restorations also is more sensitive to the defects caused by aliasing artifacts, colored noise, and ringing near sharp edges. In this paper, these visual defects are characterized, and methods for suppressing them are presented. It is demonstrated how the visual quality of fidelity-maximized images can be improved when (1) the image-gathering system is specifically designed to enhance the performance of the image-restoration algorithm, and (2) the Wiener filter is combined with interactive Gaussian smoothing, synthetic high edge enhancement, and nonlinear tone-scale transformation. The nonlinear transformation is used primarily to enhance the spatial details that are often obscurred when the normally wide dynamic range of natural radiance fields is compressed into the relatively narrow dynamic range of film and other displays.

Huck, Friedrich O.↗
NASA NTRSDate not supplied

Classifying Agnostic Biosignatures using Raman, VNIR, and Elemental Data

How can we use our current wealth of terrestrial data, encompassing biogenic and abiogenic systems, to determine the distinguishing properties of life? SCOBI (Statistical Classification of Biosignature Information) uses machine learning techniques to algorithmically identify combinations of measurements that are “indicative of life”. A set of ~1000 observations, comprising elemental abundance, isotopic fractionation, VNIR reflectance, and (in progress) Raman spectra, have been assembled from existing literature and databases. The observations cover systems classified as “indicative alive” (e.g., cells, vegetation), “indicative non-alive” (e.g., fossils, teeth), “mixed indicative” (e.g., soil, pond water), or “non-indicative” (e.g., rocks, meteorites). VNIR data was preprocessed by linear interpolation from 400-2100 nm and smoothed with a Savitzky-Golay filter. To limit the amount of Earth-biochemistry-specific (non-agnostic) information included, the first five spectral features extracted were number of peaks, number of troughs, mean reflectance, mean peak width, and broadest peak width. To help further emphasize agnostic biosignatures, Earth-specific features such as chlorophylls have been manually flagged so that feature importance with and without them can be compared. Classifiers including k-nearest neighbors (KNN), Gaussian Naïve Bayes (GNB), logistic regression (LR), random forest (RF), and support vector machine (SVM) were implemented, as was a combination voting classifier. Performance metrics included false positive rates, false negative rates, and AUC with 50-50 test/train splits (Monte Carlo simulations). Key takeaways from this stage, prior to the inclusion of Raman spectra, are (1) the overall success rate of 0.933 AUC was most heavily influenced by the elemental abundance data; and (2) VNIR reflectance had the lowest classification performance with 0.52 AUC (58% of objects correctly classified). The next steps are to complete integration of Raman spectral data and to improve the approach to pre-processing and feature extraction for both types of spectral data, such as automated baseline removal, whole spectrum matching, and dimensionality reduction.

Biosignatures↗
NASA NTRS1988

The lowest ionization potentials of Al2

Potential curves for the lowest two electronic states (X 2 sigma g + and A 2 pi u) of Al2(+) were computed using complete active space SCF/multireference CI wave functions and large Gaussian basis sets. The lowest observable vertical ionization potential (to Al2(+) X 2 sigma g +) of the Al2 X 3 pi u ground state is calculated to occur around 6.1 eV, in excellent agreement with the experimental range of 6.0 to 6.42 eV obtained in recent cluster ionization studies by Cox and co-workers. The second vertical ionization potential (to Al2(+) A 2 pi u) occurs near 6.4 eV, also within the experimental range. The adiabatic IP of 5.90 eV is in good agreement with the value of 5.8 to 6.1 eV deduced by Hanley and co-workers from the difference in thresholds between collision induced dissociation processes of Al3(+). The computed IP values are somewhat larger than those deduced from branching ratios in cluster fragmentation experiments by Jarrold and co-workers. The observation of an ionization threshold below 6.42 eV is shown to be incompatible with an Al2 ground electronic state assignment of 3 sigma g -, but the separation between the two lowest states of Al2 is so small that it is likely that both are populated in the experiments, so that this does not provide unambiguous support for the recent theoretical assignment of the ground state as 3 pi u.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2021

The Dynamic of Complex Two-Phase Mixtures During Coarsening: From Dendritic to Bicontinuous Mixtures

Dendrites are tree-like crystals that frequently form during the solidification of castings. In systems with nearly isotropic interfacial energy, these dendrites possess secondary and sometimes even tertiary arms, and lead to morphologically complex two-phase mixtures. Dealloying, in which one component is selectively leached, leads to nanoporous structures that have proven useful for many technological applications. These bicontinuous two-phase mixtures are topologically and morphologically complex with spatially varying mean and Gaussian curvatures. The evolution of these structures during coarsening is of considerable technological importance since the morphology of the dendrites and nanoporous structure sets the properties of the resulting material. During coarsening, the total surface area per volume of the mixture decreases, thus resulting in an increase in the size-scale of the two-phase mixture. The evolution of interfacial morphology of these interconnected two-phase mixtures is being studied using large-scale computer simulations. The simulations are used to develop a theory of coarsening in these systems that will elucidate the impact of the complicated morphology found in dendritic and nanoporous systems on the coarsening process. The experiments employ four-dimensional X-ray microtomography, which enables the morphology of the two-phase mixture to be followed in three dimensions and time (an additional dimension). The experiments also serve as a test of simulations and theory of the coarsening process, as well as to provide initial conditions for the simulations. Through these experiments and theory, we aim to develop a comprehensive description of coarsening in these complex and technologically important two-phase mixtures.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Intensity Conserving Spectral Fitting

The detailed shapes of spectral line profiles provide valuable information about the emitting plasma, especially when the plasma contains an unresolved mixture of velocities, temperatures, and densities. As a result of finite spectral resolution, the intensity measured by a spectrometer is the average intensity across a wavelength bin of non-zero size. It is assigned to the wavelength position at the center of the bin. However, the actual intensity at that discrete position will be different if the profile is curved, as it invariably is. Standard fitting routines (spline, Gaussian, etc.) do not account for this difference, and this can result in significant errors when making sensitive measurements. Detection of asymmetries in solar coronal emission lines is one example. Removal of line blends is another. We have developed an iterative procedure that corrects for this effect. It can be used with any fitting function, but we employ a cubic spline in a new analysis routine called Intensity Conserving Spline Interpolation (ICSI). As the name implies, it conserves the observed intensity within each wavelength bin, which ordinary fits do not. Given the rapid convergence, speed of computation, and ease of use, we suggest that ICSI be made a standard component of the processing pipeline for spectroscopic data.

Klimchuk, J. A.↗
DOE OSTI2021

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Iterative LQG Controller Design Through Closed-Loop Identification

This paper presents an iterative Linear Quadratic Gaussian (LQG) controller design approach for a linear stochastic system with an uncertain open-loop model and unknown noise statistics. This approach consists of closed-loop identification and controller redesign cycles. In each cycle, the closed-loop identification method is used to identify an open-loop model and a steady-state Kalman filter gain from closed-loop input/output test data obtained by using a feedback LQG controller designed from the previous cycle. Then the identified open-loop model is used to redesign the state feedback. The state feedback and the identified Kalman filter gain are used to form an updated LQC controller for the next cycle. This iterative process continues until the updated controller converges. The proposed controller design is demonstrated by numerical simulations and experiments on a highly unstable large-gap magnetic suspension system.

Hsiao, Min-Hung↗
DOE OSTI2021

Spatially resolved x-ray detection with photonic crystal scintillators

We study the self-collimation phenomenon in photonic crystals (PhC) of wide bandgap materials for ultra-fast and high spatial resolution x-ray detection. We work on various heavy inorganic scintillators: BaF 2 , GaN, ZnO, CsI:Tl, NaI:Tl, LYSO, WO 4 compounds, and plastic scintillators. Conventional scintillator detectors do not rely on a direct detection mechanism; hence, they require intricate design and fabrication processes. We offer a PhC design to observe self-collimation phenomena and overcome the ongoing spatial resolution challenges with these types of materials. We investigate the photonic band diagrams and iso-frequency contours. Fourier transforms based on finite-difference time-domain and frequency domain simulations are done for verifying and analyzing the self-collimation with the selected material. Light extraction efficiency at the PhC–air interface, depending on the truncation distance from the excitation point, is measured. Beam divergence values are calculated at 1 mm propagation distance. The vertical field profiles are obtained to observe the confinement. For the spatial resolution analysis, cross-sectional beam profiles have been examined. Furthermore, Gaussian envelopes are fitted to beam profiles for a consistent data analysis, and full-width-at-half-maximum values are considered. As a result, we theoretically prove and demonstrate the spatially resolved x-ray detection at the sub-micrometer level for a wide range of scintillator materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fast ground-state-to-ground-state separation of small ion crystals

Rapid separation of linear crystals of trapped ions into different subsets is critical for realizing trapped ion quantum computing architectures where ions are rearranged in trap arrays to achieve all-to-all connectivity between qubits. Here we introduce a general theoretical framework that can be used to describe the separation of same-species and mixed-species crystals into smaller subsets. The framework relies on an efficient description of the evolution of Gaussian motional states under quadratic Hamiltonians that only requires a special solution of the classical equations of motion of the ions to describe their quantum evolution under the influence of a time-dependent applied potential and the ions' mutual Coulomb repulsion. We provide time-dependent applied potentials suitable for separation of a mixed-species three-ion crystal on timescales similar to that of free expansion driven by Coulomb repulsion, with all modes along the crystal axis starting and ending close to their ground states. Three separately confined mixed-species ions can be combined into a crystal held in a single well without energy gain by time-reversal of this separation process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

The LQR/LTR procedure for multivariable feedback control design

This paper provides a tutorial overview of the linear quadratic Gaussian with loop transfer recovery (LQG/LTR) design procedure for linear multivariable feedback systems. LQR/LTR is interpreted as the solution of a specific weighted H-squared tradeoff between transfer functions in the frequency domain. Properties of this solution are examined for both minimum-phase and nonminimum-phase systems. This leads to a formal weight augmentation procedure for the minimum-phase case which permits essentially arbitrary specification of system sensitivity functions in terms of the weights. While such arbitrary specifications are not possible for nonminimum-phase problems, a direct relationship between weights and sensitivities is developed for nonminimum-phase SISO and certain nonminimum-phase MIMO cases which guides the weight selection process.

Stein, Gunter↗
DOE OSTI2021

Influence of Disorder and State Filling on Charge-Transfer-State Absorption and Emission Spectra

We conduct comprehensive temperature-dependent measurements of the charge-transfer-(CT) state photocurrent and emission spectra for two organic small molecule donor:fullerene (C 60 ) acceptor bulk heterojunction solar cells. We reveal that the CT spectral width and position are affected by static energetic disorder in the blend, especially evident at low temperatures. The relative contributions of the static and dynamic disorder broadening in the CT spectra are effectively extracted through consideration of a Gaussian CT energetic distribution. However, electroluminescence (EL) spectra can only be interpreted when injected carriers reach thermal equilibrium sites within the disordered density of states and emission occurs from the lowest possible CT energy. For the blend with the smaller energetic disorder, this is the case near room temperature; for the other blend with larger static disorder, carriers fail to reach thermal equilibrium sites even at room temperature and EL spectra need to be interpreted with care. For example, in the latter case, the effect of energetic disorder might not be apparent from EL spectra because the lowest energy sites are not participating. Nonetheless, these states contribute to the photocurrent generation-recombination and energy-loss processes and thus demand accurate characterization, which we show is feasible through temperature-dependent external quantum-efficiency measurements.

14 SOLAR ENERGY↗
DOE OSTI2021

A framework for data-driven solution and parameter estimation of PDEs using conditional generative adversarial networks

We employ and adapt the image-to-image translation concept based on conditional generative adversarial networks (cGAN) for learning a forward and an inverse solution operator of partial differential equations (PDEs). We focus on steady-state solutions of coupled hydromechanical processes in heterogeneous porous media and present the parameterization of the spatially heterogeneous coefficients, which is exceedingly difficult using standard reduced-order modeling techniques. We show that our framework provides a speed-up of at least 2,000 times compared to a finite-element solver and achieves a relative root-mean-square error (r.m.s.e.) of less than 2% for forward modeling. For inverse modeling, the framework estimates the heterogeneous coefficients, given an input of pressure and/or displacement fields, with a relative r.m.s.e. of less than 7%, even for cases where the input data are incomplete and contaminated by noise. The framework also provides a speed-up of 120,000 times compared to a Gaussian prior-based inverse modeling approach while also delivering more accurate results.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

The reflection of impulses from a nonlinear random sea

A model of the reflection of radar impulses from the sea at near-vertical incidence is used to account for non-Gaussian ocean waves statistics. The joint probability density function (pdf), of wave height and slope, is calculated according to the theory of Longuet-Higgins (1963) on the distribution of variables in a 'weakly nonlinear' random era. The long-crested approximation is made, a Phillips wave spectrum is assumed, and the Gram-Charlier series is truncated after skewness terms. It is found that the height and height-slope skewness coefficients bear the ratio 1:2 and that the derived impulse response and conditional cross section versus wave height are in excellent agreement with previous observations. Finally, it is suggested that the empirically determined and theoretically predicted sea state bias be corrected for in the routine processing of satellite radar altimeter data.

Jackson, F. C.↗
NASA NTRS2008

A New Bond Albedo for Performing Orbital Debris Brightness to Size Transformations

We have developed a technique for estimating the intrinsic size distribution of orbital debris objects via optical measurements alone. The process is predicated on the empirically observed power-law size distribution of debris (as indicated by radar RCS measurements) and the log-normal probability distribution of optical albedos as ascertained from phase (Lambertian) and range-corrected telescopic brightness measurements. Since the observed distribution of optical brightness is the product integral of the size distribution of the parent [debris] population with the albedo probability distribution, it is a straightforward matter to transform a given distribution of optical brightness back to a size distribution by the appropriate choice of a single albedo value. This is true because the integration of a powerlaw with a log-normal distribution (Fredholm Integral of the First Kind) yields a Gaussian-blurred power-law distribution with identical power-law exponent. Application of a single albedo to this distribution recovers a simple power-law [in size] which is linearly offset from the original distribution by a constant whose value depends on the choice of the albedo. Significantly, there exists a unique Bond albedo which, when applied to an observed brightness distribution, yields zero offset and therefore recovers the original size distribution. For physically realistic powerlaws of negative slope, the proper choice of albedo recovers the parent size distribution by compensating for the observational bias caused by the large number of small objects that appear anomalously large (bright) - and thereby skew the small population upward by rising above the detection threshold - and the lower number of large objects that appear anomalously small (dim). Based on this comprehensive analysis, a value of 0.13 should be applied to all orbital debris albedo-based brightness-to-size transformations regardless of data source. Its prima fascia genesis, derived and constructed from the current RCS to size conversion methodology (SiBAM Size-Based Estimation Model) and optical data reduction standards, assures consistency in application with the prior canonical value of 0.1. Herein we present the empirical and mathematical arguments for this approach and by example apply it to a comprehensive set of photometric data acquired via NASA's Liquid Mirror Telescopes during the 2000-2001 observing season.

Mulrooney, Mark K.↗
DOE OSTI2023

Improving Dark Energy Constraints Using Low-Redshift Large-Scale Structures

The primary goal of this project was to improve constraints on dark energy measurements by improving our ability to extract cosmological information from low redshift large-scale structures. PI Clowe's project's primary aim was to reduce the bias in measurements of the masses of clusters of galaxies to a level where the evolution of the cluster mass function can be used in the Vera Rubin Observatory's Legacy Survey of Space and Time Dark Energy Science Collaboration survey to improve the accuracy of the measurement of dark energy and other cosmological parameters. Co-PI Seo's project studied observational systematics affecting large-scale clustering of galaxies, which will be used to improve dark energy constraints from the Dark Energy Spectroscopic Instrument (DESI). The cluster lensing project employed a series of simulations and observations of clusters of galaxies to test numerous potential systematic errors in cluster mass measurements using weak gravitational lensing as the accuracy of current weak lensing measurements are more than order of magnitude worse than what is required to use clusters of galaxies for accurate determination of dark energy parameters. PI Clowe and group developed and analyzed simulations to test for, and correct biases introduced in, the weak lensing measurement process. Finally, PI Clowe and group developed a method of detecting clusters using galaxy overdensities and applied the method to the BLISS and DES surveys. The success of spectroscopic dark energy mission such as the extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI) will depend on a thorough understanding of various observational systematics in the target density fluctuations that would give rise to spurious, non-cosmological signals. PI Seo and group developed a deep learning, artificial neural network (ANN) technique that modeled and mitigated such effects, aimed at deriving more robust galaxy clustering signals not only for the baryon acoustic oscillation feature and redshift-space distortions but also for primordial non-Gaussianity constraint.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

FY2021 Improvements to the New CTH Code Verification & Validation Test Suite

Over the past few years, the CTH multiphysics hydrocode has overhauled its software quality and testing processes, implementing current best practices in software quality and building a robust V&V test suite comprised of traditional hydrocode verification problems, including ASC Tri-Lab Test Suite and Enhanced Tri-Lab Test Suite problems, as well as validation problems for some of CTH’s most frequently used equations of state, materials models, and other key capabilities. Substantial progress towards building this new test suite was made in FY19 and FY20. In FY21, the test suite has been expanded to include verification and validation tests of the Steinberg-Guinan-Lund (ST) viscoplastic model and the Johnson Cook (JFRAC) fracture model. Additionally, two new verification tests were added, covering hydrodynamics and high explosive (HE) modeling capabilities: the Kidder Gaussian density problem and the Escape of HE Products (EHEP) problem from the Tri-Lab Test Suite. This report discusses each of these test problems in detail. Verification test results are compared to analytic solutions. Validation test results are compared to experimental data. Wherever possible, convergence or mesh refinement studies are included. Additionally, while implementing the Kidder verification problem, a bug was identified that affects the use of tables to initialize pressure or density in 1D or 2D calculations. A brief discussion of the bug and its fix is included. CTH demonstrates good performance overall on the new test suite problems. Simulation results showed good agreement with analytic solutions for the Kidder problem, with convergence rates ranging between 1.8 and sub-linear, and relatively good agreement for the EHEP problem, though convergence rates for pressure and density were nearly 0. The ST and JFRAC strain rate loading verification tests show good agreement with analytic solutions. Likewise, CTH simulation results show good agreement with experimental validation data, including Taylor rod impact testing, for the materials tested. Future V&V work will focus on adding 2D and 3D versions of existing verification tests as well as adding validation tests of other frequently used capabilities such as other fracture models.

42 ENGINEERING↗
DOE OSTI2020

Improving qubit readout with hidden Markov models

We demonstrate the application of pattern recognition algorithms via hidden Markov models (HMM) for qubit readout. This scheme provides a state-path trajectory approach capable of detecting qubit-state transitions and makes for a robust classification scheme with higher starting-state assignment fidelity than when compared to a multivariate Gaussian or a support vector machine scheme. Therefore, the method also eliminates the qubit-dependent readout time optimization requirement in current schemes. Using a HMM state discriminator we estimate fidelities reaching the ideal limit. Unsupervised learning gives access to transition matrix, priors, and IQ distributions, providing a toolbox for studying qubit-state dynamics during strong projective readout.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗