Search NASA⌕ Search

SEARCH · Search NASA

Results for “aliasing”

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 163 records · Page 9

Investigation of spectral analysis techniques for randomly sampled velocimetry data

It is well known that velocimetry (LV) generates individual realization velocity data that are randomly or unevenly sampled in time. Spectral analysis of such data to obtain the turbulence spectra, and hence turbulence scales information, requires special techniques. The 'slotting' technique of Mayo et al, also described by Roberts and Ajmani, and the 'Direct Transform' method of Gaster and Roberts are well known in the LV community. The slotting technique is faster than the direct transform method in computation. There are practical limitations, however, as to how a high frequency and accurate estimate can be made for a given mean sampling rate. These high frequency estimates are important in obtaining the microscale information of turbulence structure. It was found from previous studies that reliable spectral estimates can be made up to about the mean sampling frequency (mean data rate) or less. If the data were evenly samples, the frequency range would be half the sampling frequency (i.e. up to Nyquist frequency); otherwise, aliasing problem would occur. The mean data rate and the sample size (total number of points) basically limit the frequency range. Also, there are large variabilities or errors associated with the high frequency estimates from randomly sampled signals. Roberts and Ajmani proposed certain pre-filtering techniques to reduce these variabilities, but at the cost of low frequency estimates. The prefiltering acts as a high-pass filter. Further, Shapiro and Silverman showed theoretically that, for Poisson sampled signals, it is possible to obtain alias-free spectral estimates far beyond the mean sampling frequency. But the question is, how far? During his tenure under 1993 NASA-ASEE Summer Faculty Fellowship Program, the author investigated from his studies on the spectral analysis techniques for randomly sampled signals that the spectral estimates can be enhanced or improved up to about 4-5 times the mean sampling frequency by using a suitable prefiltering technique. But, this increased bandwidth comes at the cost of the lower frequency estimates. The studies further showed that large data sets of the order of 100,000 points, or more, high data rates, and Poisson sampling are very crucial for obtaining reliable spectral estimates from randomly sampled data, such as LV data. Some of the results of the current study are presented.

Sree, Dave↗

Preliminary solutions for the lunar gravity field from analysis of lunar orbiter tracking data

Knowledge of the gravitation field, in combination with surface topography, provides one of the principal means of inferring the internal structure of a planetary body. Previous analyses of the lunar gravitational field have been based on data from the Lunar Orbiters, the Apollo subsatellites, and the low altitude passes of the Apollo spacecraft. Recently, Konopliv et al. have reanalyzed all available Lunar Orbiter and Apollo subsatellite tracking data, producing a 60th degree and order solution. In preparation for the Clementine Mission to the Moon, we have also initiated a reanalysis of the Lunar Orbiter and Apollo subsatellite data. Our reanalysis takes advantage of advanced force and measurement modeling techniques as well as modern computational facilities. We applied the least squares collocation technique which stabilizes the behavior of the solution and high degree and order. The extension of the size of the field reduces the aliasing coming from the omitted portion of the gravitational field. This is especially important for the analysis of the tracking data from the Lunar Orbiters, as the periapse heights frequently ranged from 50 to 100 km.

Lemoine, F. G.↗

Grayscale/resolution trade-off for text: Model predictions and psychophysical results for letter confusion and letter discrimination

In a series of papers presented in 1994, we examined the grayscale/resolution trade-off for natural images displayed on devices with discrete pixellation, such as AMLCD's. In the present paper we extend this study by examining the grayscale/resolution trade-off for text images on discrete-pixel displays. Halftoning in printing is an example of the grayscale/resolution trade-off. In printing, spatial resolution is sacrificed to produce grayscale. Another example of this trade-off is the inherent low-pass spatial filter of a CRT, caused by the point-spread function of the electron beam in the phosphor layer. On a CRT, sharp image edges are blurred by this inherent low-pass filtering, and the block noise created by spatial quantization is greatly reduced. A third example of this trade-off is text anti-aliasing, where grayscale is used to improve letter shape, size and location when rendered at a low spatial resolution. There are additional implications for display system design from the grayscale/resolution trade-off. For example, reduced grayscale can reduce system costs by requiring less complexity in the framestore, allowing the use of lower cost drivers, potentially increasing data transfer rates in the image subsystem, and simplifying the manufacturing processes that are used to construct the active matrix for AMLCD (active-matrix liquid-crystal display) or AMTFEL (active-matrix thin-film electroluminescent) devices. Therefore, the study of these trade-offs is important for display designers and manufacturing and systems engineers who wish to create the highest performance, lowest cost device possible. Our strategy for investigating this trade-off is to generate a set of simple test images, manipulate grayscale and resolution, predict discrimination performance using the ViDEOS(Sarnoff) Human Vision Model, conduct an empirical study of discrimination using psychophysical procedures, and verify the computational results using the psychophysical results.

Gille, Jennifer↗

Bicycle ergometer instrumentation to determine muscle and bone forces during exercise

It is hypothesized that bone loss experienced by astronauts in zero gravity conditions may be curtailed by appropriate exercise. According to Wolf's law, bone regenerates when muscles produce stresses by pulling on the bone during daily activity and/or exercise on Earth. to use this theory to prevent or decrease bone loss, one needs to quantify musculoskeletal loads and relate them to bone density changes. In the context of the space program, it is desirable to determine musculoskeletal loads during exercise (using the bicycle ergometer in this case) so that one may make similar measurements on Earth and in space. In this manner, load measurements on Earth may be used as reference to generate similar loads during exercise in space. The work reported in this document entails a musculoskeletal load measurement system that, when complete, will provide forces at muscle insertion points and other contact points, on bone. This data will be used by Dr. Beth A. Todd, who is also a SSF working with Dr. Shackelford, as input to a finite element model of bone sections to determine stress distributions. A bicycle ergometer has been instrumented to measure parameters needed to determine musculoskeletal forces during exercise. A primary feature of the system is its compactness. It uses small/light sensors without line-of-sight requirements. The system developed includes sensors, signal processing, a data acquisition system, and software to collect the data. The sensors used include optical encoders to measure position and orientation of the pedal (foot), accelerometers to determine kinematic parameters of the shank and thigh, load cells to measure pedal forces on the sagittal plane, and EMG probes to measure muscle activity. The signals are processed using anti-aliasing filters and amplifiers. The sensors' output is digitized using 30 channels of a board mounted inside a 486 class PC. A program sets the data acquisition parameters and collects data during a time period specified by the user. The data is put directly into a file on the hard disk in binary form. The 30 channels are sampled at 200 KHz, and each 30 channel scan is done at a rate of 1000 Hz. The instrumented ergometer has been flown in the KC-135 zero-gravity (zero-g) flight to collect information needed to determine musculoskeletal forces under these conditions. Similar information has been collected in 1-g conditions for comparision with the results from the zero-g case. At this time, the sets of data from both experiments are being processed. An existing methodology will be used to determine the kinematic parameters of the shank and thigh using accelerometer and encoder data. This methodology was developed during the fellow's previous NASA/ASEE fellowship and thanks to a Director's Grant. In the future, a methodology to determine the musculoskeletal forces using Newton's Law of Motion and optimization techniques will be developed to determine forces exerted by particular muscles.

Figueroa, Fernando↗

Evidence for eddy formation in the eastern Arabian Sea during the northeast monsoon

The seasonal formation of a large (500-800 km diameter) anticyclonic eddy in the upper 300-400 m of the eastern Arabian Sea during the northest monsoon period (December-April) is indicated fom hydrographic and satellite altimetry sea level observations, as well as from numerical model experiments. The center of the eddy circulation is approximately 10 deg N, 70 deg E, just to the west of the north-south Laccadive Island chain. In this paper the eddy is called the Laccadive High (LH). In some ways it is like a mirrorlike counterpart to the Great Whirl that develops during the southwest monsoon of the Somali coast (western Arabian Sea). The LH occurs at the same latitude but on the opposite side of the basin during the reversed monsoon. It is different from the Great Whirl, however, in its formation process, its intensity, and its decay. The hydrographic data obtained from surveys all during a single season give sufficiently close station spacing to allow reasonable contouring of the geopotential surfaces and of the properties within and around the LH region with minimum time aliasing. The Geostat altimeter record extends over 4 years, during which the seasonal variability of the LH indicates a dynamic relief of approximately 15-20 cm, which is in good agreement with the hydrographics observations. The altimetry time series also suggests a westward translation of the LH by January with a subsequent dissipation in midbasin. The model used is a wind-forced three-layer primitive equation model which depicts a LH agreement with the timing, position, and amplitude of both the hydrographic and altimetric measurements. The numerical simulation includes a passive tracer located in the Western Bay of Bengal; the western advection of the tracer around the south coasts of Sri Lanka and India in December and January is consistent with the appearance of low-salinity water observed to extend into the Arabian Sea during this period. The modeling studies suggest that both local and remote forcing are important in formation of the LH.

Bruce, John G.↗

The resolution capability of an irregularly sampled dataset: With application to Geosat altimeter data

A formalism is presented for determining the wavenumber-frequency transfer function associated with an irregularly sampled multidimensional dataset. This transfer function reveals the filtering characteristics and aliasing patterns inherent in the sample design. In combination with information about the spectral characteristics of the signal, the transfer function can be used to quantify the spatial and temporal resolution capability of the dataset. Application of the method to idealized Geosat altimeter data (i.e., neglecting measurement errors and data dropouts) concludes that the Geosat orbit configuration is capable of resolving scales of about 3 deg in latitude and longitude by about 30 days.

Chelton, Dudley B.↗

Comparison of TOPEX sea surface heights and tide gauge sea levels

TOPEX sea surface height data from the first 300 days of the mission are compared to sea level data from 71 tide gauges. The initial comparison uses sea surface height data processed according to standard procedures as defined in the users handbook. It is found that the median correlations for island and for coastal tide gauges are 0.53 and 0.42, respectively. The analogous root mean square (RMS) differences between the two data sets are 7.9 and 10.4 cm. The comparisons improve significantly when a 60-day harmonic is fit to the differences and removed. This period captures aliased M(sub 2) and S(sub 2) tidal energy that is not removed by the tide model. Making this correction and smoothing the sea surface height data over 25-km along-track segments results in median correlations of 0.58 and 0.46 for the islands and coastal stations, and median RMS differences of 5.8 and 7.7 cm, respectively. Removing once per revolution signals from the sea surface heights results in degraded comparisons with the sea levels. It is also found that a number of stations have poor comparisons due to propagating signals that introduce temporal lags between the altimeter and tide gauge time series. A final comparison is made by eliminating stations where this propagation effect is large, discarding two stations that are suspected to have problems with the sea level data, smoothing over 10-day intervals, and restricting attention to islands gauges. This results in a set of 552 data pairs that have a correlation of 0.66 and a RMS difference of 4.3 cm. The conclusion is that on timescales longer than about 10 days the RMS sea surface height errors are less than or of the order of several centimeters.

Mitchum, Gary T.↗

A preliminary tidal analysis of TOPEX/POSEIDON altimetry

Approximately 12 months of data from the TOPEX/POSEIDON satellite altimeter mission are analyzed for the major short-period oceanic tides. A harmonic analysis is performed on data captured within bins defined on a deep-ocean grid, which, owing to tidal aliasing considerations, must have a relatively coarse spatial resolution. Our analysis is in terms of corrections to the Schwiderski and Cartwright-Ray models, and it confirms many of the Schwiderski differences previously reported by Cartwright and Ray. Our differences with respect to the Geosat-based Cartwright and Ray model form a sectorial pattern in M(sub 2) with high/low differences separated roughly 180 deg in longitude. We suggest that these sectorial errors were most likely induced by Geosat's relatively large orbit error. Comparisons to independent data validate the improved TOPEX/POSEIDON solutions; in-situ 'ground truth' shows M(sub 2) root mean square (RMS) differences of 4.10 cm (Schwiderski), 3.86 cm (cartwright and Ray), 2.63 cm (this paper). Global rates of energy dissipation confirm earlier estimates for M(sub 2), and show improved agreement with satellite tracking studies for K(sub 1) and S(sub 2). These preliminary exercises confirm that TOPEX/POSEIDON should results in a new generation of improved global tidal models.

Schrama, E. J. O.↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

Diurnal and semidiurnal variations in Earth orientation determined from LAGEOS laser ranging

Variations in universal time and polar motion due to ocean tides at nearly diurnal and nearly semidiurnal frequencies are determined from analysis of laser ranging to the LAGEOS satellite over the period from 1987 to 1992. The adjusted diurnal tides were K(sub 1), S(sub 1), P(sub 1), O(sub 1), and Q(sub 1), while the semidiurnal tides were K(sub 2), S(sub 2), M(sub 2), and N(sub 2). A formulation was used that explicitly separated prograde and retrograde terms in the polar motion in order to eliminate aliasing from the singularity of retrograde wobble with long period orbit error and nutation. The results are compared to other experimentally derived observations from very long baseline interferometry (VLBI) and with predictions from an ocean tide model. The results of this study were well with those from the VLBI studies, typically at the 2-3 microsecond level in universal time and 30-50 microarc sec (muas) level in polar motion. The agreement with predictions from the ocean tide model were roughly a factor of 3 worse for UT although better for polar motion, particularly in the diurnal band.

Watkins, Michael M.↗

Artifacts Of Spectral Analysis Of Instrument Readings

Report presents experimental and theoretical study of some of artifacts introduced by processing outputs of two nominally identical low-frequency-reading instruments; high-sensitivity servo-accelerometers mounted together and operating, in conjunction with signal-conditioning circuits, as seismometers. Processing involved analog-to-digital conversion with anti-aliasing filtering, followed by digital processing including frequency weighting and computation of different measures of power spectral density (PSD).

Wise, James H.↗

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim↗

Rapid calibration of seven-hole probes

This paper summarizes the major conclusions and some of the key supporting analyses resulting from the calibration and application of two small seven hole probes at NASA Lewis Research Center. These probes can produce reasonably accurate and rapid surveys of unknown steady flow fields which may include flow angles up to 70 degrees and Mach numbers up to 0.8. The probes were calibrated with both 'complete' and 'reduced' test matrices. Both types of test matrices produced similar results suggesting the the reduced matrices are adequate for most purposes. The average accuracy fo the calibration was about the same as that achieved in previous seven hole probe calibrations. At the higher Mach numbers, the calibration was sensitive to the diameter of the free jet in the calibration facility. Over a narrow angular range at the higher Mach numbers, the system had serious repeatability problems. This lack or repeatability apparently results from aliasing of high frequency (20 to 40 Hz) noise with the data acquisition system sampling frequency of 10 Hz. Analyses show that these noise frequencies are probably not related to airflow dynamics in the connecting tubing.

Ericksen, A. L.↗

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

Blaisdell, Gregory A.↗

Plasma Heating and Flow in an Auroral Arc

We report direct observations of the three-dimensional velocity distribution of selected topside ionospheric ion species in an auroral context between 500 and 550 km altitude. We find heating transverse to the local magnetic field in the core plasma, with significant heating of 0(+), He(+), and H(+), as well as tail heating events that occur independently of the core heating. The 0(+) velocity distribution departs from bi-Maxwellian, at one point exhibiting an apparent ring-like shape. However, these observations are shown to be aliased within the auroral arc by temporal variations that arc not well-resolved by the core plasma instrument. The dc electric field measurements reveal superthermal plasma drifts that are consistent with passage of the payload through a series of vortex structures or a larger scale circularly polarized hydromagnetic wave structure within the auroral arc. The dc electric field also shows that impulsive solitary structures, with a frequency spectrum in the ion cyclotron frequency range, occur in close correlation with the tail heating events. The drift and core heating observations lend support to the idea that core ion heating is driven at low altitudes by rapid convective motions imposed by the magnetosphere. Plasma wave emissions at ion frequencies and parallel heating of the low-energy electron plasma are observed in conjunction with this auroral form; however, the conditions are much more complex than those typically invoked in previous theoretical treatments of superthermal frictional heating. The observed ion heating within the arc clearly exceeds that expected from frictional heating for the light ion species H(+) and He(+), and the core distributions also contain hot transverse tails, indicating an anomalous transverse heat source.

Moore, T. E.↗

Structure and Rotation of the Solar Interior: Initial Results from the MDI Medium-L Program

The medium-l program of the Michelson Doppler Imager instrument on board SOHO provides continuous observations of oscillation modes of angular degree, l, from 0 to approximately 300. The data for the program are partly processed on board because only about 3% of MDI observations can be transmitted continuously to the ground. The on-board data processing, the main component of which is Gaussian-weighted binning, has been optimized to reduce the negative influence of spatial aliasing of the high-degree oscillation modes. The data processing is completed in a data analysis pipeline at the SOI Stanford Support Center to determine the mean multiplet frequencies and splitting coefficients. The initial results show that the noise in the medium-l oscillation power spectrum is substantially lower than in ground-based measurements. This enables us to detect lower amplitude modes and, thus, to extend the range of measured mode frequencies. This is important for inferring the Sun's internal structure and rotation. The MDI observations also reveal the asymmetry of oscillation spectral lines. The line asymmetries agree with the theory of mode excitation by acoustic sources localized in the upper convective boundary layer. The sound-speed profile inferred from the mean frequencies gives evidence for a sharp variation at the edge of the energy-generating core. The results also confirm the previous finding by the GONG (Gough et al., 1996) that, in a thin layer just beneath the convection zone, helium appears to be less abundant than predicted by theory. Inverting the multiplet frequency splittings from MDI, we detect significant rotational shear in this thin layer. This layer is likely to be the place where the solar dynamo operates. In order to understand how the Sun works, it is extremely important to observe the evolution of this transition layer throughout the 11-year activity cycle.

Kosovichev, A. G.↗

Accelerometer Data Analysis and Presentation Techniques

The NASA Lewis Research Center's Principal Investigator Microgravity Services project analyzes Orbital Acceleration Research Experiment and Space Acceleration Measurement System data for principal investigators of microgravity experiments. Principal investigators need a thorough understanding of data analysis techniques so that they can request appropriate analyses to best interpret accelerometer data. Accelerometer data sampling and filtering is introduced along with the related topics of resolution and aliasing. Specific information about the Orbital Acceleration Research Experiment and Space Acceleration Measurement System data sampling and filtering is given. Time domain data analysis techniques are discussed and example environment interpretations are made using plots of acceleration versus time, interval average acceleration versus time, interval root-mean-square acceleration versus time, trimmean acceleration versus time, quasi-steady three dimensional histograms, and prediction of quasi-steady levels at different locations. An introduction to Fourier transform theory and windowing is provided along with specific analysis techniques and data interpretations. The frequency domain analyses discussed are power spectral density versus frequency, cumulative root-mean-square acceleration versus frequency, root-mean-square acceleration versus frequency, one-third octave band root-mean-square acceleration versus frequency, and power spectral density versus frequency versus time (spectrogram). Instructions for accessing NASA Lewis Research Center accelerometer data and related information using the internet are provided.

Rogers, Melissa J. B.↗

Sarnoff JND Vision Model for Flat-Panel Design

This document describes adaptation of the basic Sarnoff JND Vision Model created in response to the NASA/ARPA need for a general-purpose model to predict the perceived image quality attained by flat-panel displays. The JND model predicts the perceptual ratings that humans will assign to a degraded color-image sequence relative to its nondegraded counterpart. Substantial flexibility is incorporated into this version of the model so it may be used to model displays at the sub-pixel and sub-frame level. To model a display (e.g., an LCD), the input-image data can be sampled at many times the pixel resolution and at many times the digital frame rate. The first stage of the model downsamples each sequence in time and in space to physiologically reasonable rates, but with minimum interpolative artifacts and aliasing. Luma and chroma parts of the model generate (through multi-resolution pyramid representation) a map of differences-between test and reference called the JND map, from which a summary rating predictor is derived. The latest model extensions have done well in calibration against psychophysical data and against image-rating data given a CRT-based front-end. THe software was delivered to NASA Ames and is being integrated with LCD display models at that facility,

Brill, Michael H.↗