Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wavelet functions”

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 55 records · Page 3

Complexity Induced Anisotropic Bimodal Intermittent Turbulence in Space Plasmas

The "physics of complexity" in space plasmas is the central theme of this exposition. It is demonstrated that the sporadic and localized interactions of magnetic coherent structures arising from the plasma resonances can be the source for the coexistence of nonpropagating spatiotemporal fluctuations and propagating modes. Non-Gaussian probability distribution functions of the intermittent fluctuations from direct numerical simulations are obtained and discussed. Power spectra and local intermittency measures using the wavelet analyses are presented to display the spottiness of the small-scale turbulent fluctuations and the non-uniformity of coarse-grained dissipation that can lead to magnetic topological reconfigurations. The technique of the dynamic renormalization group is applied to the study of the scaling properties of such type of multiscale fluctuations. Charged particle interactions with both the propagating and nonpropagating portions of the intermittent turbulence are also described.

Chang, Tom↗

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI↗

Wavelet Analyses of F/A-18 Aeroelastic and Aeroservoelastic Flight Test Data

Time-frequency signal representations combined with subspace identification methods were used to analyze aeroelastic flight data from the F/A-18 Systems Research Aircraft (SRA) and aeroservoelastic data from the F/A-18 High Alpha Research Vehicle (HARV). The F/A-18 SRA data were produced from a wingtip excitation system that generated linear frequency chirps and logarithmic sweeps. HARV data were acquired from digital Schroeder-phased and sinc pulse excitation signals to actuator commands. Nondilated continuous Morlet wavelets implemented as a filter bank were chosen for the time-frequency analysis to eliminate phase distortion as it occurs with sliding window discrete Fourier transform techniques. Wavelet coefficients were filtered to reduce effects of noise and nonlinear distortions identically in all inputs and outputs. Cleaned reconstructed time domain signals were used to compute improved transfer functions. Time and frequency domain subspace identification methods were applied to enhanced reconstructed time domain data and improved transfer functions, respectively. Time domain subspace performed poorly, even with the enhanced data, compared with frequency domain techniques. A frequency domain subspace method is shown to produce better results with the data processed using the Morlet time-frequency technique.

Martin J Brenner↗

Image Compression Algorithm Altered to Improve Stereo Ranging

A report discusses a modification of the ICER image-data-compression algorithm to increase the accuracy of ranging computations performed on compressed stereoscopic image pairs captured by cameras aboard the Mars Exploration Rovers. (ICER and variants thereof were discussed in several prior NASA Tech Briefs articles.) Like many image compressors, ICER was designed to minimize a mean-square-error measure of distortion in reconstructed images as a function of the compressed data volume. The present modification of ICER was preceded by formulation of an alternative error measure, an image-quality metric that focuses on stereoscopic-ranging quality and takes account of image-processing steps in the stereoscopic-ranging process. This metric was used in empirical evaluation of bit planes of wavelet-transform subbands that are generated in ICER. The present modification, which is a change in a bit-plane prioritization rule in ICER, was adopted on the basis of this evaluation. This modification changes the order in which image data are encoded, such that when ICER is used for lossy compression, better stereoscopic-ranging results are obtained as a function of the compressed data volume.

Kiely, Aaron↗

Specification of Surface Roughness Over Oceans in the GLAS Fourth Order GCM

The surface roughness height z sub 0 is a parameter which measures the effect of surface irregularities on the mean profiles which occur in the atmospheric surface layer. Over land surfaces, z sub 0 is a function of soil composition, plant canopy type and distribution and other factors such as buildings and structures. Over open water, z sub 0 depends only on the distribution of waves, wavelets and other surface irregularities. These, in turn, are thought to be determined by the mean wind within the surface layer. Large and Pond have derived an empirical relationship between the mean wind speed U sub 10 at the ten meter level over open water and the neutral drag coefficient C sub DN. Their relationship covers the range U sub 10 or = M/S. Earlier, Kondo had computed a relationship between U sub 10 and C sub DN on the basis of laboratory observations which extended down to the range of infinitely small wind speeds. An interpretation between the relationships of Large and Pond for large wind speed and of Kondo for small wind speed is presented.

Helfand, H. M.↗

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗

The Impact of Stellar Variability on the Detection of Transiting Earth-like Planets

Stellar variability represents a key limitation on the detectability of weak transit signatures of small, Earth-size planets with space-based missions such as NASA's Kepler Mission or ESA's PLATO Mission. The expectations for the performance of the Kepler Mission in the face of solar-like variability were predicated on observations of the Sun with the Variability of Solar Irradiance and Gravity Oscillations (VIRGO) instrument aboard NASA and ESA's joint Solar and Heliospheric Observatory (SOHO) mission, which was launched in December 1995. Analyses of SOHO data indicated that solar-like variability would contribute approx.10 ppm to Kepler's noise budget at 6.5 ours, a typical grazing transit timescale. This proved to be optimistic as the typical 12th mag G2V star in the Kepler field of view exhibited 20 ppm of stellar variability-related noise at 6.5 hours, making the detection of true Earth-Sun analogs more difficult. While the consequent increase in the total noise at 6.5 hours from 20 ppm to 30 ppm could have been overcome by extending the mission to a total duration of 7 to 8 years. Unfortunately, a second reaction wheel failed in May 2013, ending the Kepler primary mission. Nevertheless, Kepler proved to be a veritable cornucopia of science results, both for exoplanets and for astrophysics. The phenomenal photometric precision and continuous observations required in order to identify small, rocky transiting planets enables the study of a large range of phenomena contributing to stellar variability for many thousands of solar-like stars in Kepler's field of view in exquisite detail. These effects range from less than 1 ppm acoustic oscillations on timescales from a few minutes and longward, to flares on timescales of hours, to spot-induced modulation on timescales of days to weeks to activity cycles on timescales of months to years. Kepler discovered over 2600 validated and confirmed exoplanets, and measured the pressure-mode oscillations of over 15000 stars over the course of its mission and the reconstituted K2 mission. We present the adaptive, wavelet-based matched filter used for both the Kepler and the TESS missions, and discuss the connections between this detection algorithm and fundamental detection theory, which allows the detector to function as a noise characterization engine. This provides a dynamical measurement of the photometric noise at transit timescales and thus, a key performance metric for missions like PLATO and Kepler. We give an overview of the stellar variability we see across the full range of spectral types observed by Kepler, from the cool, small red M stars to the hot, large late A stars, both in terms of amplitude as well as timescale.

transit surveys↗

The Impact of Stellar Variability on the Detection of Transiting Earth-Like Planets

Stellar variability represents a key limitation on the detectability of weak transit signatures of small, Earth-size planets with space-based missions such as NASA's Kepler Mission or ESA's PLATO Mission. The expectations for the performance of the Kepler Mission in the face of solar-like variability were predicated on observations of the Sun with the Variability of Solar Irradiance and Gravity Oscillations (VIRGO) instrument aboard NASA and ESA's joint Solar and Heliospheric Observatory (SOHO) mission, which was launched in December 1995. Analyses of SOHO data indicated that solar-like variability would contribute ~10 ppm to Kepler's noise budget at 6.5 ours, a typical grazing transit timescale. This proved to be optimistic as the typical 12th mag G2V star in the Kepler field of view exhibited 20 ppm of stellar variability-related noise at 6.5 hours, making the detection of true Earth-Sun analogs more difficult. While the consequent increase in the total noise at 6.5 hours from 20 ppm to 30 ppm could have been overcome by extending the mission to a total duration of 7 to 8 years. Unfortunately, a second reaction wheel failed in May 2013, ending the Kepler primary mission. Nevertheless, Kepler proved to be a veritable cornucopia of science results, both for exoplanets and for astrophysics. The phenomenal photometric precision and continuous observations required in order to identify small, rocky transiting planets enables the study of a large range of phenomena contributing to stellar variability for many thousands of solar-like stars in Kepler's field of view in exquisite detail. These effects range from <1 ppm acoustic oscillations on timescales from a few minutes and longward, to flares on timescales of hours, to spot-induced modulation on timescales of days to weeks to activity cycles on timescales of months to years. Kepler discovered over 2600 validated and confirmed exoplanets, and measured the pressure-mode oscillations of over 15000 stars over the course of its mission and the reconstituted K2 mission. We present the adaptive, wavelet-based matched filter used for both the Kepler and the TESS missions, and discuss the connections between this detection algorithm and fundamental detection theory, which allows the detector to function as a noise characterization engine. This provides a dynamical measurement of the photometric noise at transit timescales and thus, a key performance metric for missions like PLATO and Kepler. We give an overview of the stellar variability we see across the full range of spectral types observed by Kepler, from the cool, small red M stars to the hot, large late A stars, both in terms of amplitude as well as timescale.

exoplanets↗

Normal and abnormal tissue identification system and method for medical images such as digital mammograms

A system and method for analyzing a medical image to determine whether an abnormality is present, for example, in digital mammograms, includes the application of a wavelet expansion to a raw image to obtain subspace images of varying resolution. At least one subspace image is selected that has a resolution commensurate with a desired predetermined detection resolution range. A functional form of a probability distribution function is determined for each selected subspace image, and an optimal statistical normal image region test is determined for each selected subspace image. A threshold level for the probability distribution function is established from the optimal statistical normal image region test for each selected subspace image. A region size comprising at least one sector is defined, and an output image is created that includes a combination of all regions for each selected subspace image. Each region has a first value when the region intensity level is above the threshold and a second value when the region intensity level is below the threshold. This permits the localization of a potential abnormality within the image.

Heine, John J.↗

Unsteady Heat-Flux Measurements of Second-Mode Instability Waves in a Hypersonic Boundary Layer

In this paper we report on the application of the atomic layer thermopile (ALTP) heat- flux sensor to the measurement of laminar-to-turbulent transition in a hypersonic flat plate boundary layer. The centerline of the flat-plate model was instrumented with a streamwise array of ALTP sensors and the flat-plate model was exposed to a Mach 6 freestream over a range of unit Reynolds numbers. Here, we observed an unstable band of frequencies that are associated with second-mode instability waves in the laminar boundary layer that forms on the flat-plate surface. The measured frequencies, group velocities, phase speeds, and wavelengths of these instability waves are in agreement with data previously reported in the literature. Heat flux time series, and the Morlet-wavelet transforms of them, revealed the wave-packet nature of the second-mode instability waves. In addition, a laser-based radiative heating system was developed to measure the frequency response functions (FRF) of the ALTP sensors used in the wind tunnel test. These measurements were used to assess the stability of the sensor FRFs over time and to correct spectral estimates for any attenuation caused by the finite sensor bandwidth.

Kergerise, Michael A.↗

Wavelets and spacetime squeeze

It is shown that the wavelet is the natural language for the Lorentz covariant description of localized light waves. A model for covariant superposition is constructed for light waves with different frequencies. It is therefore possible to construct a wave function for light waves carrying a covariant probability interpretation. It is shown that the time-energy uncertainty relation (Delta(t))(Delta(w)) is approximately 1 for light waves is a Lorentz-invariant relation. The connection between photons and localized light waves is examined critically.

Han, D.↗

Development of a Mars Surface Imager

The Mars Surface Imager (MSI) is a multispectral, stereoscopic, panoramic imager that allows imaging of the full scene around a Mars lander from the lander body to the zenith. It has two functional components: panoramic imaging and sky imaging. In the most recent version of the MSI, called PIDDP-cam, a very long multi-line color CCD, an innovative high-performance drive system, and a state-of-the-art wavelet image compression code have been integrated into a single package. The requirements for the flight version of the MSI and the current design are presented.

Squyres, Steve W.↗

Profiles of Gamma-Ray Bursts and Their Component Pulses

One physically informative regularity of their otherwise heterogeneous ensemble, is that many Gamma-Ray Bursts consist of well defined pulses. To objectively quantify the temporal structure of BATSE bursts, we have developed an automatic modeling procedure that separates overlapping pulses and determines the energy-dependence of the pulse-shape parameters. No binning of photon arrival times is needed, so when applied to time-tagged events (TTE) the procedure captures variability information down to the shortest time scales present in the raw data. Maximizing the Bayesian likelihood function Pr(data/model) yields estimates of the model parameters, including the number of pulses present, and allows intercomparison of models of different forms. As with any nonlinear optimization, good initial guesses are crucial to avoid convergence to undesirable local minima. We find excellent initial pulse decompositions by wavelet-denoising a cumulative distribution of the raw photon arrival data; differentiation then gives a time profile mostly free of the systematic effects of degraded resolution (as in ordinary Fourier smoothing) and binning. We present statistical information on pulse rise-time, decay-time, peakedness, and amplitudes, plus their energy dependences - both within a single burst and for a large ensemble of bursts.

Scargle, Jeff D.↗

Distant Cluster Hunting: A Comparison of X-Ray and Optical Cluster Detection Techniques and Catalogs from the ROSAT Optical X-Ray Survey - II

We present and analyze the optical and X-ray catalogs of moderate-redshift cluster candidates from the ROSA TOptical X-Ray Survey, or ROXS. The survey covers the sky area contained in the fields of view of 23 deep archival ROSA T PSPC pointings, 4.8 square degrees. The cross-correlated cluster catalogs were con- structed by comparing two independent catalogs extracted from the optical and X-ray bandpasses, using a matched-filter technique for the optical data and a wavelet technique for the X-ray data. We cross-identified cluster candidates in each catalog. As reported in Paper 1, the matched-filter technique found optical counter- parts for at least 60% (26 out of 43) of the X-ray cluster candidates; the estimated redshifts from the matched filter algorithm agree with at least 7 of 1 1 spectroscopic confirmations (Az 5 0.10). The matched filter technique. with an imaging sensitivity of ml N 23, identified approximately 3 times the number of candidates (155 candidates, 142 with a detection confidence >3 u) found in the X-ray survey of nearly the same area. There are 57 X-ray candidates, 43 of which are unobscured by scattered light or bright stars in the optical images. Twenty-six of these have fairly secure optical counterparts. We find that the matched filter algorithm, when applied to images with galaxy flux sensitivities of mI N 23, is fairly well-matched to discovering z 5 1 clusters detected by wavelets in ROSAT PSPC exposures of 8000-60,000 s. The difference in the spurious fractions between the optical and X-ray (30%) and IO%, respectively) cannot account for the difference in source number. In Paper I, we compared the optical and X-ray cluster luminosity functions and we found that the luminosity functions are consistent if the relationship between X-ray and optical luminosities is steep (Lx o( L&f). Here, in Paper 11, we present the cluster catalogs and a numerical simulation of the ROXS. We also present color-magnitude plots for several of the cluster candidates, and examine the prominence of the red sequence in each. We find that the X-ray clusters in our survey do not all have a prominent red sequence. We conclude that while the red sequence may be a distinct feature in the color-magnitude plots for virialized massive clusters, it may be less distinct in lower mass clusters of galaxies at even moderate redshifts. Multiple, complementary methods of selecting and defining clusters may be essential, particularly at high redshift where all methods start to run into completeness limits, incomplete understanding of physical evolution, and projection effects.

Donahue, Megan↗

Locally-Based Kernal PLS Smoothing to Non-Parametric Regression Curve Fitting

We present a novel smoothing approach to non-parametric regression curve fitting. This is based on kernel partial least squares (PLS) regression in reproducing kernel Hilbert space. It is our concern to apply the methodology for smoothing experimental data where some level of knowledge about the approximate shape, local inhomogeneities or points where the desired function changes its curvature is known a priori or can be derived based on the observed noisy data. We propose locally-based kernel PLS regression that extends the previous kernel PLS methodology by incorporating this knowledge. We compare our approach with existing smoothing splines, hybrid adaptive splines and wavelet shrinkage techniques on two generated data sets.

Rosipal, Roman↗

Low-complexity wavelet filter design for image compression

Image compression algorithms based on the wavelet transform are an increasingly attractive and flexible alternative to other algorithms based on block orthogonal transforms. While the design of orthogonal wavelet filters has been studied in significant depth, the design of nonorthogonal wavelet filters, such as linear-phase (LP) filters, has not yet reached that point. Of particular interest are wavelet transforms with low complexity at the encoder. In this article, we present known and new parameterizations of the two families of LP perfect reconstruction (PR) filters. The first family is that of all PR LP filters with finite impulse response (FIR), with equal complexity at the encoder and decoder. The second family is one of LP PR filters, which are FIR at the encoder and infinite impulse response (IIR) at the decoder, i.e., with controllable encoder complexity. These parameterizations are used to optimize the subband/wavelet transform coding gain, as defined for nonorthogonal wavelet transforms. Optimal LP wavelet filters are given for low levels of encoder complexity, as well as their corresponding integer approximations, to allow for applications limited to using integer arithmetic. These optimal LP filters yield larger coding gains than orthogonal filters with an equivalent complexity. The parameterizations described in this article can be used for the optimization of any other appropriate objective function.

Majani, E.↗

A New Method for Nonlinear and Nonstationary Time Series Analysis and Its Application to the Earthquake and Building Response Records

A new method for analyzing nonlinear and nonstationary data has been developed. The key part of the method is the Empirical Mode Decomposition method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF). An IMF is defined as any function having the same numbers of zero-crossing and extrema, and also having symmetric envelopes defined by the local maxima and minima respectively. The IMF also admits well-behaved Hilbert transform. This decomposition method is adaptive, and, therefore, highly efficient. Since the decomposition is based on the local characteristic time scale of the data, it is applicable to nonlinear and nonstationary processes. With the Hilbert transform, the Intrinsic Mode Functions yield instantaneous frequencies as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum, Example of application of this method to earthquake and building response will be given. The results indicate those low frequency components, totally missed by the Fourier analysis, are clearly identified by the new method. Comparisons with Wavelet and window Fourier analysis show the new method offers much better temporal and frequency resolutions.

Huang, Norden E.↗