SEARCH · Engineering Papers
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.
Through Small Business Innovation Research funds from the Stennis Space Center, MathSoft, Inc., developed a system that can provide the building blocks for signal analysis and rapid prototyping. The product is the result of work to help NASA develop a complete understanding propulsion test data by using time frequency displays, automatic estimation and denoising, and data analysis plots for wavelet decomposition.
We present the high time resolution observations of one of the Langmuir wave packets obtained in the source region of a solar type III radio burst. This wave packet satisfies the threshold condition of the supersonic modulational instability, as well as the criterion of a collapsing Langmuir soliton, i.e., the spatial scale derived from its peak intensity is less than that derived from its short time scale. The spectrum of t his wave packet contains an intense spectral peak at local electron plasma frequency, f(sub pe) and relatively weaker peaks at 2f(sub pe) and 3f(sub pe). We apply the wavelet based bispectral analysis technique on this wave packet and compute the bicoherence between its spectral components. It is found that the bicoherence exhibits two peaks at (approximately f(sub pe), approximately f(sub pe)) and (approximately f(sub pe) approximately 2f(sub pe)), which strongly suggest that the spectral peak at 2f(sub pe) probably corresponds to the second harmonic radio emission, generated as a result of the merging of antiparallel propagating Langmuir waves trapped in the collapsing Langmuir soliton, and, the spectral peak at 3f(sub pe) probably corresponds to the third harmonic radio emission, generated as a result of merging of a trapped Langmuir wave and a second harmonic electromagnetic wave.
GRACE satellites have detected regional‐scale preseismic, coseismic, and postseismic gravity changes associated with great earthquakes during the GRACE era (2002‐2017). Earthquakes also excite global‐scale transient gravity changes associated with free oscillations that may be discerned for a few days. In this study, we examine such global gravity changes due to Earth's free oscillations and quantify how they affect GRACE measurements. We employ the normal mode formalism to synthesize the global gravity changes after the 2004 Sumatra earthquake and simulate the (gravitational) free oscillation signals manifested in the GRACE K‐band ranging (KBR) measurements. Using the Kaula orbit perturbation theory, we show how GRACE inter‐satellite distances are perturbed through a complex coupling of eigenfrequencies of the normal modes with the Earth's rotation rate and the GRACE satellites' orbital frequency. It is found that a few gravest normal modes can generate range‐rate perturbations as large as 0.2 μm/s, which are comparable to actual errors of GRACE KBR ranging and accelerometer instruments. Wavelet time‐frequency analysis of the GRACE KBR residual data in December 2004 reveals the existence of a significant transient signal after the 2004 Sumatra earthquake. This transient signal is characterized by a frequency of ~0.022 mHz that could be potentially associated with the largest excitation due to the “football” mode of the Earth's free oscillation. However, the results are also affected by low‐frequency noise of the GRACE accelerometers. Improved space‐borne gravitational instrumentation may open new opportunities to study the Earth's interior and earthquakes independently from global seismological analysis.
Clouds play an important role in the radiation budget of the atmosphere. A good understanding of how clouds interact with solar radiation is necessary when considering their effects in both general circulation models and climate models. This study examined the radiative properties of clouds in both an inhomogeneous cloud system, and a simplified cloud system through the use of a Monte Carlo model. The purpose was to become more familiar with the radiative properties of clouds, especially absorption, and to investigate the excess absorption of solar radiation from observations over that calculated from theory. The first cloud system indicated that the absorptance actually decreased as the cloud's inhomogeneity increased, and that cloud forcing does not indicate any changes. The simplified cloud system looked at two different cases of absorption of solar radiation in the cloud. The absorptances calculated from the Monte Carlo is compared to a correction method for calculating absorptances and found that the method can over or underestimate absorptances at cloud edges. Also the cloud edge effects due to solar radiation points to a possibility of overestimating the retrieved optical depth at the edge, and indicates a possible way to correct for it. The effective cloud fraction (Ne) for a long time has been calculated from a cloud's reflectance. From the reflectance it has been observed that the N, for most cloud geometries is greater than the actual cloud fraction (Nc) making a cloud appear wider than it is optically. Recent studies we have performed used a Monte Carlo model to calculate the N, of a cloud using not only the reflectance but also the absorptance. The derived Ne's from the absorptance in some of the Monte Carlo runs did not give the same results as derived from the reflectance. This study also examined the inhomogeneity of clouds to find a relationship between larger and smaller scales, or wavelengths, of the cloud. Both Fourier transforms and wavelet transforms were used to analyze the liquid water content of marine stratocumulus clouds taken during the ASTEX project. From the analysis it was found that the energy in the cloud is not uniformly distributed but is greater at the larger scales than at the smaller scales. This was determined by examining the slope of the power spectrum, and by comparing the variability at two scales from a wavelet analysis.
Explore the source record for details and available documents.
Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.
Wavelets are the central idea of a broad framework for thinking about displaying, and analyzing data - not just a new specialized technique. They are localized basis functions that provide an extremely flexible and efficient representation of time series. Current developments in the theory and application of wavelets are yielding powerful multiresolution techniques for the analysis of time series and image data. Wavelet methods can effect noise removal. characterization of stochastic behavior, depiction of power spectrum time-variations, and data compression - for time series data, image data, and data of higher dimensions. I review published applications of wavelets in astronomy.
We present evidence that the quasi-periodic oscillations (QPO) and very low frequency noise (VLFN) characteristic of many accretion sources are different aspects of the same physical process. We analyzed a long, high time resolution EXOSAT observation of the low-mass X-ray binary (LMXB) Sco X-1. The X-ray luminosity varies stochastically on time scales from milliseconds to hours. The nature of this variability - as quantified with both power spectrum analysis and a new wavelet technique, the scalegram - agrees well with the dripping handrail accretion model, a simple dynamical system which exhibits transient chaos. In this model both the QPO and VLFN are produced by radiation from blobs with a wide size distribution, resulting from accretion and subsequent diffusion of hot gas, the density of which is limited by an unspecified instability to lie below a threshold.
Traditionally, Fourier Transforms have been utilized for performing signal analysis and representation. But although it is straightforward to reconstruct a signal from its Fourier transform, no local description of the signal is included in its Fourier representation. To alleviate this problem, Windowed Fourier transforms and then wavelet transforms have been introduced, and it has been proven that wavelets give a better localization than traditional Fourier transforms, as well as a better division of the time- or space-frequency plane than Windowed Fourier transforms. Because of these properties and after the development of several fast algorithms for computing the wavelet representation of any signal, in particular the Multi-Resolution Analysis (MRA) developed by Mallat, wavelet transforms have increasingly been applied to signal analysis problems, especially real-life problems, in which speed is critical. In this paper we present and compare efficient wavelet decomposition algorithms on different parallel architectures. We report and analyze experimental measurements, using NASA remotely sensed images. Results show that our algorithms achieve significant performance gains on current high performance parallel systems, and meet scientific applications and multimedia requirements. The extensive performance measurements collected over a number of high-performance computer systems have revealed important architectural characteristics of these systems, in relation to the processing demands of the wavelet decomposition of digital images.
Estimation of the state of the atmosphere with the Kalman filter remains a distant goal because of high computational cost of evolving the error covariance for both linear and nonlinear systems. Wavelet approximation is presented here as a possible solution that efficiently compresses both global and local covariance information. We demonstrate the compression characteristics on the the error correlation field from a global two-dimensional chemical constituent assimilation, and implement an adaptive wavelet approximation scheme on the assimilation of the one-dimensional Burger's equation. In the former problem, we show that 99%, of the error correlation can be represented by just 3% of the wavelet coefficients, with good representation of localized features. In the Burger's equation assimilation, the discrete linearized equations (tangent linear model) and analysis covariance are projected onto a wavelet basis and truncated to just 6%, of the coefficients. A nearly optimal forecast is achieved and we show that errors due to truncation of the dynamics are no greater than the errors due to covariance truncation.
This paper describes the application of neural network adaptive wavelets for fault diagnosis of space station power system. The method combines wavelet transform with neural network by incorporating daughter wavelets into weights. Therefore, the wavelet transform and neural network training procedure become one stage, which avoids the complex computation of wavelet parameters and makes the procedure more straightforward. The simulation results show that the proposed method is very efficient for the identification of fault locations.
This paper presents a physics-aligned framework for fault diagnosis in multi-phase power-electronic systems using cycle-synchronous windowing and multi-domain features derived from Fourier, wavelet, and Hilbert–Huang representations. While both logistic regression and multilayer perceptron (MLP) models achieve perfect performance under standard evaluation, blind unseen testing reveals a critical failure in a baseline MLP. This is shown to arise from model selection based on validation accuracy. Using validation-loss-based selection restores correct unseen performance and improves confidence. Feature ablation shows that Fourier and wavelet features dominate, while computational analysis indicates that feature extraction, particularly HHT, governs runtime.
Image registration is one of the basic image processing operations in remote sensing. With the increase in the number of images collected every day from different sensors, automated registration of multi-sensor/multi-spectral images has become an important issue. A wide range of registration techniques has been developed for many different types of applications and data. The objective of this paper is to present an automatic registration algorithm which uses a multiresolution analysis procedure based upon the wavelet transform. The procedure is completely automatic and relies on the grey level information content of the images and their local wavelet transform modulus maxima. The registration algorithm is very simple and easy to apply because it needs basically one parameter. We have obtained very encouraging results on test data sets from the TM and SPOT sensor images of forest, urban and agricultural areas.
"Multi-Scale Extended Kalman" (MSEK) is a novel model-based control approach recently found to be effective for suppressing combustion instabilities in gas turbines. A control law formulated in this approach for fuel modulation demonstrated steady suppression of a high-frequency combustion instability (less than 500Hz) in a liquid-fuel combustion test rig under engine-realistic conditions. To make-up for severe transport-delays on control effect, the MSEK controller combines a wavelet -like Multi-Scale analysis and an Extended Kalman Observer to predict the thermo-acoustic states of combustion pressure perturbations. The commanded fuel modulation is composed of a damper action based on the predicted states, and a tones suppression action based on the Multi-Scale estimation of thermal excitations and other transient disturbances. The controller performs automatic adjustments of the gain and phase of these actions to minimize the Time-Scale Averaged Variances of the pressures inside the combustion zone and upstream of the injector. The successful demonstration of Active Combustion Control with this MSEK controller completed an important NASA milestone for the current research in advanced combustion technologies.
With the rapid expansion of large-scale photovoltaic (PV) plants, it is paramount for solar stakeholders to understand the reliability and efficiency of their plants to inform maintenance decisions, increase production, and understand the design factors that impact performance. Diagnosing underperformance in PV plants is challenging due to the relatively few monitoring points with respect to the large geographic footprint of the plant. This work introduces a cutting-edge method that transforms the analysis and management of key factors influencing PV plant performance, including performance loss rate (PLR), recoverable soiling, and major system changes. Identifying these factors is critical for deriving actionable insights. Leveraging advanced analytical techniques such as wavelet transformation, robust regression, and extreme point analysis, this approach provides a nuanced understanding of these factors. This method has been tested across two synthetic datasets and one real dataset, consistently surpassing existing benchmarks by achieving a lower median mean absolute error and reduced error variability across all comparable components.
New challenges concerning system health-monitoring and life-extending robust controls for the Ultra-Efficient Engine Technology Project, as well as other advanced engine and power system concepts at NASA and elsewhere, have renewed the control community s interest in smart, model-based methods. In particular, these challenges have further motivated efforts at the NASA Glenn Research Center to exploit the versatility and superiority of the dynamic features extraction of multiscale analysis for controls--such as with "wavelets" and "wavelet filter-banks.' The accomplishments reported herein pertain to the active suppression of combustion instabilities in liquid-fuel combustors via fuel modulation. The fundamentals and initial success of this innovation were reported for a unique demonstration of active combustion control (a research collaboration of NASA Glenn with Pratt & Whitney and the United Technologies Research Center, UTRC). This demonstration, conducted in 2002 at UTRC on the NASA single nozzle rig (SNR) combustor, was the first known suppression of high-frequency instability with a liquid-fueled combustor. The SNR is based on a high-powered military engine combustor that exhibited well-known instabilities.
While wavelet transforms have been used in the analysis of seismic signals, the methods used for the compression of seismic signals have been limited to Linear Predictive Coding and Transform Coding.
Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤 CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ CDM) and 90% (𝑤 CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 (Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ CDM) and Ω m = 0.273 ±0.029 (𝑤 CDM). In the context of the 𝑤 CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2𝜎.