Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular value 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Spherical Harmonic Analysis of Particle Velocity Distribution Function: Comparison of Moments and Anisotropies using Cluster Data

This paper presents a spherical harmonic analysis of the plasma velocity distribution function using high-angular, energy, and time resolution Cluster data obtained from the PEACE spectrometer instrument to demonstrate how this analysis models the particle distribution function and its moments and anisotropies. The results show that spherical harmonic analysis produced a robust physical representation model of the velocity distribution function, resolving the main features of the measured distributions. From the spherical harmonic analysis, a minimum set of nine spectral coefficients was obtained from which the moment (up to the heat flux), anisotropy, and asymmetry calculations of the velocity distribution function were obtained. The spherical harmonic method provides a potentially effective "compression" technique that can be easily carried out onboard a spacecraft to determine the moments and anisotropies of the particle velocity distribution function for any species. These calculations were implemented using three different approaches, namely, the standard traditional integration, the spherical harmonic (SPH) spectral coefficients integration, and the singular value decomposition (SVD) on the spherical harmonic methods. A comparison among the various methods shows that both SPH and SVD approaches provide remarkable agreement with the standard moment integration method.

Gurgiolo, Chris↗
NASA NTRS2010

Equivalent Dipole Vector Analysis for Detecting Pulmonary Hypertension

Various 12-lead ECG criteria have been established to detect right ventricular hypertrophy as a marker of pulmonary hypertension (PH). While some criteria offer good specificity they lack sensitivity because of a low prevalence of positive findings in the PH population. We hypothesized that three-dimensional equivalent dipole (ED) model could serve as a better detection tool of PH. We enrolled: 1) 17 patients (12 female, 5 male, mean age 57 years, range 19-79 years) with echocardiographically detected PH (systolic pulmonary arterial pressure greater than 35 mmHg) and no significant left ventricular disease; and 2) 19 healthy controls (7 female, 12 male, mean age 44, range 31-53 years) with no known heart disease. In each subject we recorded a 5-minute high-resolution 12-lead conventional ECG and constructed principal signals using singular value decomposition. Assuming a standard thorax dimension of an adult person with homogenous and isotropic distribution of thorax conductance, we determined moving equivalent dipoles (ED), characterized by the 3D location in the thorax, dipolar strength and the spatial orientation, in time intervals of 5 ms. We used the sum of all ED vectors in the second half of the QRS complex to derive the amplitude of the right-sided ED vector (RV), if the orientation of ED was to the right side of the thorax, and in the first half the QRS to derive the amplitude of the left-sided vector (LV), if the orientation was leftward. Finally, the parameter RV/LV ratio was determined over an average of 256 complexes. The groups differed in age and gender to some extent. There was a non-significant trend toward higher RV in patients with PH (438 units 284) than in controls (280 plus or minus 140) (p = 0.066) but the overlap was such that RV alone was not a good predictor of PH. On the other hand, the RV/LV ratio was a better predictor of PH, with 11/17 (64.7%) of PH patients but only in 1/19 (5.3%) control subjects having RV/LV ratio greater than or equal to 0.70 (p less than 0.001). The use of ED for evaluating PH shows good specificity at a reasonable sensitivity. The results are limited due to the small study groups and differences in age and gender, but further investigations are warranted, including of ED's diagnostic accuracy for PH versus that of other proposed ECG and VCG criteria.

Harlander, Matevz↗
NASA NTRS2010

A Framework for Propagation of Uncertainties in the Kepler Data Analysis Pipeline

The Kepler space telescope is designed to detect Earth-like planets around Sun-like stars using transit photometry by simultaneously observing 100,000 stellar targets nearly continuously over a three and a half year period. The 96-megapixel focal plane consists of 42 charge-coupled devices (CCD) each containing two 1024 x 1100 pixel arrays. Cross-correlations between calibrated pixels are introduced by common calibrations performed on each CCD requiring downstream data products access to the calibrated pixel covariance matrix in order to properly estimate uncertainties. The prohibitively large covariance matrices corresponding to the ~75,000 calibrated pixels per CCD preclude calculating and storing the covariance in standard lock-step fashion. We present a novel framework used to implement standard propagation of uncertainties (POU) in the Kepler Science Operations Center (SOC) data processing pipeline. The POU framework captures the variance of the raw pixel data and the kernel of each subsequent calibration transformation allowing the full covariance matrix of any subset of calibrated pixels to be recalled on-the-fly at any step in the calibration process. Singular value decomposition (SVD) is used to compress and low-pass filter the raw uncertainty data as well as any data dependent kernels. The combination of POU framework and SVD compression provide downstream consumers of the calibrated pixel data access to the full covariance matrix of any subset of the calibrated pixels traceable to pixel level measurement uncertainties without having to store, retrieve and operate on prohibitively large covariance matrices. We describe the POU Framework and SVD compression scheme and its implementation in the Kepler SOC pipeline.

Clarke, Bruce D.↗
NASA NTRS2012

Thermographic Imaging of the Space Shuttle During Re-Entry Using a Near Infrared Sensor

High resolution calibrated near infrared (NIR) imagery of the Space Shuttle Orbiter was obtained during hypervelocity atmospheric re-entry of the STS-119, STS-125, STS-128, STS-131, STS-132, STS-133, and STS-134 missions. This data has provided information on the distribution of surface temperature and the state of the airflow over the windward surface of the Orbiter during descent. The thermal imagery complemented data collected with onboard surface thermocouple instrumentation. The spatially resolved global thermal measurements made during the Orbiter s hypersonic re-entry will provide critical flight data for reducing the uncertainty associated with present day ground-to-flight extrapolation techniques and current state-of-the-art empirical boundary-layer transition or turbulent heating prediction methods. Laminar and turbulent flight data is critical for the validation of physics-based, semi-empirical boundary-layer transition prediction methods as well as stimulating the validation of laminar numerical chemistry models and the development of turbulence models supporting NASA s next-generation spacecraft. In this paper we provide details of the NIR imaging system used on both air and land-based imaging assets. The paper will discuss calibrations performed on the NIR imaging systems that permitted conversion of captured radiant intensity (counts) to temperature values. Image processing techniques are presented to analyze the NIR data for vignetting distortion, best resolution, and image sharpness. Keywords: HYTHIRM, Space Shuttle thermography, hypersonic imaging, near infrared imaging, histogram analysis, singular value decomposition, eigenvalue image sharpness

Zalameda, Joseph N.↗
NASA NTRS2013

Developing an Accurate CFD Based Gust Model for the Truss Braced Wing Aircraft

The increased flexibility of long endurance aircraft having high aspect ratio wings necessitates attention to gust response and perhaps the incorporation of gust load alleviation. The design of civil transport aircraft with a strut or truss-braced high aspect ratio wing furthermore requires gust response analysis in the transonic cruise range. This requirement motivates the use of high fidelity nonlinear computational fluid dynamics (CFD) for gust response analysis. This paper presents the development of a CFD based gust model for the truss braced wing aircraft. A sharp-edged gust provides the gust system identification. The result of the system identification is several thousand time steps of instantaneous pressure coefficients over the entire vehicle. This data is filtered and downsampled to provide the snapshot data set from which a reduced order model is developed. A stochastic singular value decomposition algorithm is used to obtain a proper orthogonal decomposition (POD). The POD model is combined with a convolution integral to predict the time varying pressure coefficient distribution due to a novel gust profile. Finally the unsteady surface pressure response of the truss braced wing vehicle to a one-minus-cosine gust, simulated using the reduced order model, is compared with the full CFD.

Bartels, Robert E.↗
NASA NTRS2014

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data. Part II: Using Maximum Covariance Analysis to Effectively Compare Spatiotemporal Variability of Satellite and AERONET Measured Aerosol Optical Depth

Moderate Resolution Imaging SpectroRadiometer (MODIS) and Multi-angle Imaging Spectroradiomater (MISR) provide regular aerosol observations with global coverage. It is essential to examine the coherency between space- and ground-measured aerosol parameters in representing aerosol spatial and temporal variability, especially in the climate forcing and model validation context. In this paper, we introduce Maximum Covariance Analysis (MCA), also known as Singular Value Decomposition analysis as an effective way to compare correlated aerosol spatial and temporal patterns between satellite measurements and AERONET data. This technique not only successfully extracts the variability of major aerosol regimes but also allows the simultaneous examination of the aerosol variability both spatially and temporally. More importantly, it well accommodates the sparsely distributed AERONET data, for which other spectral decomposition methods, such as Principal Component Analysis, do not yield satisfactory results. The comparison shows overall good agreement between MODIS/MISR and AERONET AOD variability. The correlations between the first three modes of MCA results for both MODIS/AERONET and MISR/ AERONET are above 0.8 for the full data set and above 0.75 for the AOD anomaly data. The correlations between MODIS and MISR modes are also quite high (greater than 0.9). We also examine the extent of spatial agreement between satellite and AERONET AOD data at the selected stations. Some sites with disagreements in the MCA results, such as Kanpur, also have low spatial coherency. This should be associated partly with high AOD spatial variability and partly with uncertainties in satellite retrievals due to the seasonally varying aerosol types and surface properties.

aerosols↗
NASA NTRS2014

Characterization of the Inner Knot of the Crab: The Site of the Gamma-Ray Flares?

Subsequent to the detections AGILE and Fermi/LAT of the gamma-ray flares from the Crab Nebula in the fall of 2010, this team has been monitoring the X-Ray emission from the Crab on a regular basis. Initially X-Ray observations took place once per month when viewing constraints allowed. More recently observations with Chandra and HST have taken place four times per year. There have been notable exceptions, e.g. in April of 2011 and March 2013 when we initiated a set of Chandra Target of opportunity observations in conjunction with bright gamma-ray flares. Often Keck observations were obtained. The aim of this program to characterize, in depth, the X-ray, optical, and infrared variations that take place in the nebula, and, by so doing, determine the regions which contribute to the harder X-ray variations and, if possible, determine the precise location within the Nebula of the origin of the gamma-ray flares. As part of this project members of the team have applied Singular Value Decomposition techniques to sequences of images in order to more accurately characterize features and their behavior. The current status of the project will be discussed highlighting studies of the inner knot and possible correlations with the gamma-ray flares.

Chandra↗
NASA NTRS2015

Characterization of the Inner Knot of the Crab: the Site of the Gamma-ray Flares?

One of the most intriguing recent discoveries has been the detection of powerful gamma-ray flares from the Crab Nebula. Such events, with a recurrence time of about once per year, can be so dramatic to make the system the brightest source in the gamma‐ray sky as occurred, e.g. in April 2011. These flares challenge our understanding of how pulsar wind nebulae work and defy current astrophysical models for particle acceleration. We present here our study of the inner knot located within a fraction of an arcsecond from the pulsar with the aim of characterizing the feature and asking if this might be the site of the origin of the gamma-ray flares. We took data using Keck, HST, and Chandra obtained as part of our multi‐wavelength campaign to identify the source of the enigmatic flares. We set an upper limit as to the gamma-ray flux from the knot. We also find that the dimensions, surface brightness, flux, etc. of the optical and infrared knot are all correlated with distance from the pulsar. This distance, in turn, varies with time. In addition to this most thorough characterization of the inner knot's properties, we examine the hypothesis that the knot may be the site of the flares by examining the knot separation versus the Fermi/LAT gamma-ray flux. Finally, as part of this research, we make use of a new approach employing singular value decomposition (SVD) for analyzing time series of images and compare the approach to more traditional methods. Our conclusions are only refined but not impacted by using the new approach.

Gamma-ray Flares↗
NASA NTRS2014

Monitoring of the Crab Nebula with Chandra and Other Observatories Including HST

Subsequent to the detections AGILE and Fermi/LAT of the gamma‐ray flares from the Crab Nebula in the fall of 2010, this team has been monitoring the X‐Ray emission from the Crab on a regular basis. X‐Ray observations have taken place typically once per month when viewing constraints allow and more recently four times per year. There have been notable exceptions, e.g. in April of 2011 and March 2013 when we initiated a set of Chandra Target of opportunity observations in conjunction with bright gamma‐ray flares. For much of the time regular HST observations were made in conjunction with the Chandra observations. The aim of this program to further characterize, in depth, the X‐Ray and optical variations that take place in the nebula, and by so doing determine the regions which contribute to the harder X‐ray variations and, if possible, determine the precise location within the Nebula of the origin of the gamma‐ray flares. As part of this project members of the team have developed Singular Value Decomposition techniques to sequences of images in order to more accurately characterize features. The current status of the project will be presented highlighting studies of the inner knot and possible correlations with the flares.

Weisskopf, Martin C.↗
NASA NTRS2016

Erasing the Milky Way: New Cleaning Technique Applied to GBT Intensity Mapping Data

We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15- and 1-h field data of the GBT observations previously presented in Masui et al. (2013) and Switzer et al. (2013), covering about 41 square degrees at 0.6 less than z is less than 1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative sub-traction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.

Wolz, L.↗
NASA NTRS2018

Decadal Seasonal Shifts of Precipitation and Temperature in TRMM and AIRS Data

We present results from an analysis of seasonal phase shifts in the global precipitation and surface temperatures. We use data from the TRMM (Tropical Rainfall Measuring Mission) Multi-satellite Precipitation Algorithm (TMPA), and the Atmospheric Infrared Sounder (AIRS) on Aqua satellite, all hosted at NASA Goddard Earth Science Data and Information Services Center (GES DISC). We explore the information content and data usability by first aggregating daily grids from the entire records of both missions to pentad (5-day) series which are then processed using Singular Value Decomposition approach. A strength of this approach is the normalized principal components that can then be easily converted from real to complex time series. Thus, we can separate the most informative, the seasonal, components and analyze unambiguously for potential seasonal phase drifts. TMPA and AIRS records represent correspondingly 20 and 15 years of data, which allows us to run simple “phase learning†from the first 5 years of records and use it as reference. The most recent 5 years are then phase-compared with the reference. We demonstrate that the seasonal phase of global precipitation and surface temperatures has been stable in the past two decades. However, a small global trend of delayed precipitation, and earlier arrival of surface temperatures seasons, are detectable at 95% confidence level. Larger phase shifts are detectable at regional level, in regions recognizable from the Eigen vectors to having strong seasonal patterns. For instance, in Central North America, including the North American Monsoon region, confident phase shifts of 1-2 days per decade are detected at 95% confidence level. While seemingly symbolic, these shifts are indicative of larger changes in the Earth Climate System. We thus also demonstrate a potential usability scenario of Earth Science Data Records curated at the NASA GES DISC in partnership with Earth Science Missions.

surface temperatures↗
NASA NTRSDate not supplied

Reduced Dimensionality Analysis of TEMPO Ozone Profile Retrievals Using the Compact Phase Space (CPSR) Algorithm

TEMPO ozone (O 3 ) profile retrievals are expected to have fidelity in the troposphere due the sensitivities of the associated averaging kernels. However, those averaging kernels are severely rank deficiency meaning that a visual inspection of the vertical structure of the averaging kernel profile sensitivities is misleading due linear dependencies in the profile. The Compact Phase Space Retrieval (CPSR) algorithm use singular value decompositions of the averaging kernels and the ‘compressed’ retrieval solution error covariance to project the transformed averaging kernels into a space that removes the linear dependencies and accounts for the solution error uncertainties. In this oral presentation and poster, we apply the CPSR dimensional reduction analysis to TEMPO and TROPOMI O 3 profile retrievals for 13:45 UTC March 29, 2024 to study the phase space characteristics of the transformed averaging kernels as a function of latitude for North America. Our results show that TEMPO generally has more phase space vertical structure in the troposphere than TROPOMI. TEMPO has four to five dominant modes, and TROPOMI has five to six dominant modes. That means that dimensional reduction can reduce the TEMPO resource requirements by ~77% and the TROPOMI requirements by ~81%. Finally, we found that after removing linear dependences and after accounting for solution uncertainties TEMPO still has sensitivities throughout the troposphere.

TEMPO↗
NASA NTRS1996

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗
NASA NTRS1996

Balanced Placement Of Vibration Sensors And Actuators

Mathematical approach that involves Hankel singular values leads to method for estimation of optimal placement of relatively few vibration sensors and actuators in multivariable flexible structure. Approach extends and complements previous developments regarding approximate decomposition, to Hankel singular values, of singular values of controllability and observability grammian matrices of multivariable flexible structure. Mathematical constructs also described in "Designing An Approximately Balanced LQG Compensator" (NPO-19000).

Gawronski, Wodek↗
NASA NTRS1987

Robust eigensystem assignment for flexible structures

An improved method is developed for eigenvalues and eigenvectors placement of a closed-loop control system using either state or output feedback. The method basically consists of three steps. First, the singular value of QR decomposition is used to generate an orthonormal basis that spans admissible eigenvector space corresponding to each assigned eigenvalue. Secondly, given a unitary matrix, the eigenvector set which best approximates the given matrix in the least-square sense and still satisfy eigenvalue cosntraints is determined. Thirdly, a unitary matrix is sought to minimize the error between the unitary matrix and the assignable eigenvector matrix. For use as the desired eigenvector set, two matrices, namely, the open-loop eigenvector matrix and its closest unitary matrix are proposed. The latter matrix generally encourages both minimum conditioning and control gains. In addition, the algorithm is formulated in real arithmetic for efficient implementation. To illustrate the basic concepts, numerical examples are included.

Juang, Jer-Nan↗
NASA NTRS1990

Robust eigensystem assignment for second-order estimators

An approach for the robust eigensystem assignment of flexible structures using full state or output feedback is developed. Using the second-order dynamic equations, the approach can assign the eigenvalues of the system via velocity and displacement feedbacks, or acceleration and velocity feedbacks. The eigenvalues and eigenvectors of the system are assigned, via the second-order eigenvalue problem for the structural system, in two steps. First, an orthonormal basis spanning the attainable closed-loop eigenvector space corresponding to each desired closed-loop eigenvalue is generated using the Singular Value or QR decompositions. Second, a sequential procedure is used to choose a set of closed-loop eigenvectors that are as close as possible to the column space of a well-conditioned target matrix. Among the possible choices of the target matrix, the closest unitary matrix to the open-loop eigenvector matrix appears to be a suitable choice. A numerical example is given to illustrate the proposed algorithm.

Juang, Jer-Nan↗
NASA NTRS1990

Robust eigensystem assignment for second-order dynamic systems

A novel approach for the robust eigensystem assignment of flexible structures using full state or output feedback is developed. Using the second-order dynamic equations, the approach can assign the eigenvalues of the system via velocity and displacement feedbacks, or acceleration and velocity feedbacks. The eigenvalues and eigenvectors of the system are assigned, via the second-order eigenvalue problem for the structural system, in two steps. First, an orthonormal basis spanning the attainable closed-loop eigenvector space corresponding to each desired closed-loop eigenvalue is generated using the Singular Value or QR decompositions. Second, a sequential procedure is used to choose a set of closed-loop eigenvectors that are as close as possible to the column space of a well-conditioned target matrix. Among the possible choices of the target matrix, the closest unitary matrix to the open-loop eigenvector matrix appears to be a suitable choice. A numerical example is given to illustrate the proposed algorithm.

Juang, Jer-Nan↗