Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral decomposition”

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 235 records · Page 13

ICER-3D Hyperspectral Image Compression Software

Software has been developed to implement the ICER-3D algorithm. ICER-3D effects progressive, three-dimensional (3D), wavelet-based compression of hyperspectral images. If a compressed data stream is truncated, the progressive nature of the algorithm enables reconstruction of hyperspectral data at fidelity commensurate with the given data volume. The ICER-3D software is capable of providing either lossless or lossy compression, and incorporates an error-containment scheme to limit the effects of data loss during transmission. The compression algorithm, which was derived from the ICER image compression algorithm, includes wavelet-transform, context-modeling, and entropy coding subalgorithms. The 3D wavelet decomposition structure used by ICER-3D exploits correlations in all three dimensions of sets of hyperspectral image data, while facilitating elimination of spectral ringing artifacts, using a technique summarized in "Improving 3D Wavelet-Based Compression of Spectral Images" (NPO-41381), NASA Tech Briefs, Vol. 33, No. 3 (March 2009), page 7a. Correlation is further exploited by a context-modeling subalgorithm, which exploits spectral dependencies in the wavelet-transformed hyperspectral data, using an algorithm that is summarized in "Context Modeler for Wavelet Compression of Hyperspectral Images" (NPO-43239), which follows this article. An important feature of ICER-3D is a scheme for limiting the adverse effects of loss of data during transmission. In this scheme, as in the similar scheme used by ICER, the spatial-frequency domain is partitioned into rectangular error-containment regions. In ICER-3D, the partitions extend through all the wavelength bands. The data in each partition are compressed independently of those in the other partitions, so that loss or corruption of data from any partition does not affect the other partitions. Furthermore, because compression is progressive within each partition, when data are lost, any data from that partition received prior to the loss can be used to reconstruct that partition at lower fidelity. By virtue of the compression improvement it achieves relative to previous means of onboard data compression, this software enables (1) increased return of hyperspectral scientific data in the presence of limits on the rates of transmission of data from spacecraft to Earth via radio communication links and/or (2) reduction in spacecraft radio-communication power and/or cost through reduction in the amounts of data required to be downlinked and stored onboard prior to downlink. The software is also suitable for compressing hyperspectral images for ground storage or archival purposes.

Xie, Hua↗

Progress on Associate-Particle Imaging Algorithms, 2020

The present work describes progress on the development of imaging algorithms that use fast neutron signatures acquired using the associated-particle imaging (API) method. The present work complements ongoing work to develop neutron source and detector hardware to enable field inspection by investigating algorithms that are capable of discriminating between critical materials or extracting three-dimensional geometrical information from single-sided or transmission measurements. The present work is divided into three approaches:(1)Iterative reconstruction of inelastic gamma-ray emissions to perform three-dimensional time-of-flight imaging in a single view in either transmission or backscatter configurations. Iterative reconstruction enables image resolution better than the inherent TOF resolution.(2)Decomposition of registered neutron and x ray radiographs into an assumed material list for each pixel in the image.(3)Material identification using full spectral analysis that includes the emergent neutron and gamma ray energies, times, and angles.Progress for each approach is summarized for fiscal year 2020.

97 MATHEMATICS AND COMPUTING↗

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.↗

Benzoic acid degradation of polyacrylonitrile fibers

The reactions of polyacrylonitrile (PAN) fibers in the presence of benzoic acid have been studied. Polyacrylonitrile fibers oxidize more readily in the presence of benzoic acid than in air at temperatures in the range of 170 C. The product decreased in solubility with extent of reaction. Gel permeation chromatography of the soluble fraction showed change in polydispersity. The insoluble product exhibited differences in weight loss as a function of decomposition temperature compared to PAN fibers. Infrared analyses of the fiber product showed absorption peaks similar to air-oxidized PAN. High-energy photoelectron spectral analysis showed a carbon-rich surface which contained oxygen and nitrogen. An air oxidized sample of fiber contained more oxygen at the surface than a fiber treated first with benzoic acid and then air oxidized.

Varma, D. S.↗

Benzoic acid degradation of polyacrylonitrile fibers

The reactions of polyacrylonitrile (PAN) fibers in the presence of benzoic acid have been studied. Polyacrylonitrile fibers oxidize more readily in the presence of benzoic acid than in air at temperatures in the range of 170 C. The product decreased in solubility with extent of reaction. Gel permeation chromatography of the soluble fraction showed change in polydispersity. The insoluble product exhibited differences in weight loss as a function of decomposition temperature compared to PAN fibers. Infrared analyses of the fiber product showed absorption peaks similar to air-oxidized PAN. High-energy photoelectron spectral analysis showed a carbon-rich surface which contained oxygen and nitrogen. An air oxidized sample of fiber contained more oxygen at the surface than a fiber treated first with benzoic acid and then air oxidized.

Varma, D. S.↗

Automatic partitioning of unstructured grids into connected components

This paper presents two partitioning schemes that guarantee connected components given a connected initial grid. Connected components are important for convergence of methods such as domain decomposition or multigrid. For many of the grids tested, the schemes produce partitions as good (in terms of number of cut edges) or better than spectral partitioning and require only modest computational resources. This paper describes the two schemes in detail and presents comparison results from a number of two and three dimensional unstructured grids.

Dagum, Leonardo↗

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single molecule kinetics. Using a stochastic Markov state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated accurately for model problems for a physically realistic mechanism for an alkene hydrogenation reaction on a single atom catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quality indicators for trace gas concentrations determined from earthlimb emission spectra

Different techniques which can be used to indicate the quality of the mixing ratios retrieved from spectral radiance measurements made by the CLAES instrument for 11 different atmospheric gases are discussed. A key to retrieval from noisy data is the Singular Value Decomposition (SVD), which is described. The use of the SVD with simulated data and radiance measurements from a balloon floating at 30 km is discussed.

Uplinger, W. G.↗

Simulations of bypass transition for spatially evolving disturbances

The spatial evolution of disturbances in plane Poiseuille flow and zero pressure gradient boundary layer flow is considered. For disturbances governed by the linearized equations, potential for significant transient growth of the amplitude is demonstrated. The maximum amplification occurs for disturbances with zero or near zero frequencies. Spatial numerical simulations of the transition scenario involving a pair of oblique waves has been conducted for both flows. A fully spectral solver using a simple but efficient fringe region technique allowed the flows to be computed with high resolution into the fully turbulent domain. A modal decomposition of the simulation results indicates that non-linear excitation of the transient growth is responsible for the rapid emergence of low-frequency structures. Physically, this corresponds to streaky flow structures, as seen from the results of a numerical amplitude expansion. Thus, this spatial transition scenario has been found to be similar to the corresponding temporal one. In the boundary layer simulations the streaks are seen to break down from what appears to be a secondary instability.

Lundbladh, A.↗

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

Implementation of an Unsteady PSP System in the NASA Transonic Dynamics Tunnel

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment.

unsteady pressure-sensitive paint↗

Implementation of an Unsteady PSP System in the NASA TDT

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment. Note: this presentation is an MP4 video with sound, color with a run time of 10 minutes 37 seconds.

unsteady pressure-sensitive paint↗

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗

On the spectral stability of time integration algorithms for a class of constrained dynamics problems

Incomplete field formulations have recently been the subject of intense research because of their potential in coupled analysis of independently modeled substructures, adaptive refinement, domain decomposition, and parallel processing. This paper discusses the design and analysis of time-integration algorithms for these formulations and emphasizes the treatment of their inter-subdomain constraint equations. These constraints are shown to introduce a destabilizing effect in the dynamic system that can be analyzed by investigating the behavior of the time-integration algorithm at infinite and zero frequencies. Three different approaches for constructing penalty-free unconditionally stable second-order accurate solution procedures for this class of hybrid formulations are presented, discussed and illustrated with numerical examples. The theoretical results presented in this paper also apply to a large family of nonlinear multibody dynamics formulations. Some of the algorithms outlined herein are important alternatives to the popular technique consisting of transforming differential/algebraic equations into ordinary differential equations via the introduction of a stabilization term that depends on arbitrary constants and that influences the computed so1ution.

Farhat, Charbel↗

An infrared spectral database for gas-phase quantitation of volatile per- and polyfluoroalkyl substances (PFAS)

We report the construction of a database of vetted infrared spectra specifically targeting volatile fluorocarbon gases that may be emitted during thermal treatment of per- and polyfluoroalkyl substances (PFAS) to assist understanding of treatment processes and improve quantification. To populate this database, protocols derived from the Pacific Northwest National Laboratory (PNNL) infrared spectral database are used, curtailing the species selection for this data set. Each spectrum in the database is a weighted average derived from 10 or more individual measurements at different partial pressures (static method) or flow rates (dissemination method) to yield good fidelity of both strong and weak infrared signatures, with each composite spectrum ranging from = 6500 cm-1 to = 600 cm-1 with an apodized resolution of 0.112 cm-1. This resolution was chosen to fully resolve all spectral features, recognizing that atmospheric pressure broadening results in nearly all ro-vibrational lines having linewidths = 0.1 cm-1. As an example case, application of the database is demonstrated via identification and quantification of dominant 1H-perfluoroheptane and perfluorohept-1-ene fluorocarbon products resulting from thermal decomposition of perfluorooctanoate (PFOA) below 450 °C.

Infrared, Gas-phase spectra, FTIR, Spectral databa↗

On the Confidence Limit of Hilbert Spectrum

Confidence limit is a routine requirement for Fourier spectral analysis. But this confidence limit is established based on ergodic theory: For stationary process, temporal average equals the ensemble average. Therefore, one can divide the data into n-sections and treat each section as independent realization. Most natural processes in general, and climate data in particular, are not stationary; therefore, there is a need for the Hilbert Spectral analysis for such processes. Here ergodic theory is no longer applicable. We propose to use various adjustable parameters in the shifting processes of the Empirical Mode Decomposition (EMD) method to obtain an ensemble of Intrinsic Mode Function 0 sets. Based on such an ensemble, we introduce a statistical measure in. a form of confidence limits for the Intrinsic Mode Functions, and consequently, the Hilbert spectra. The criterion of selecting the various adjustable parameters is based on the orthogonality test of the resulting M F sets. Length-of-day data from 1962 to 2001 will be used to illustrate this new approach. Its implication in climate data analysis will also be discussed.

Huang, Norden↗

The Application of a Boundary Integral Equation Method to the Prediction of Ducted Fan Engine Noise

The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) is considered. Governing equations for the BIEM are based on linearized acoustics and describe the scattering of incident sound by a thin, finite-length cylindrical duct in the presence of a uniform axial inflow. A classical boundary value problem (BVP) is derived that includes an axisymmetric, locally reacting liner on the duct interior. Using potential theory, the BVP is recast as a system of hypersingular boundary integral equations with subsidiary conditions. We describe the integral equation derivation and solution procedure in detail. The development of the computationally efficient ducted fan noise prediction program TBIEM3D, which implements the BIEM, and its utility in conducting parametric noise reduction studies are discussed. Unlike prediction methods based on spinning mode eigenfunction expansions, the BIEM does not require the decomposition of the interior acoustic field into its radial and axial components which, for the liner case, avoids the solution of a difficult complex eigenvalue problem. Numerical spectral studies are presented to illustrate the nexus between the eigenfunction expansion representation and BIEM results. We demonstrate BIEM liner capability by examining radiation patterns for several cases of practical interest.

Dunn, M. H.↗