Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 343 records · Page 19

Image reconstruction from multiple 1-D scans using filtered localized projection

The spatial resolution that can be attained using scanning linear arrays (consisting of discrete IR solid-state detectors) for image acquisition is considered, and a filtered local projection (FLP) method is described which efficiently combines all available scan information into one rectangular grid without the need for explicit interpolation. Mathematically, the FLP algorithm consists of a localized summation followed by an inverse-filter operation, and it has application to nonlinear restoring techniques. The present method is applied, using a linear array, to simulated data for staggered parallel scans and to multiple scan directions. Noise effects and limitations of the technique are also considered.

Frieden, B. Roy↗

Perceptually Lossless Wavelet Compression

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp -1), where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We propose a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a 'perceptually lossless' quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗

Visibility of Wavelet Quantization Noise

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp)-L , where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We describe a mathematical model to predict DWT noise detection thresholds as a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗

Visibility of wavelet quantization noise

The discrete wavelet transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that we call DWT uniform quantization noise; it is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2-lambda, where r is display visual resolution in pixels/degree, and lambda is the wavelet level. Thresholds increase rapidly with wavelet spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from lowpass to horizontal/vertical to diagonal. We construct a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

NASA Center ARC↗

Context Modeler for Wavelet Compression of Spectral Hyperspectral Images

A context-modeling sub-algorithm has been developed as part of an algorithm that effects three-dimensional (3D) wavelet-based compression of hyperspectral image data. The context-modeling subalgorithm, hereafter denoted the context modeler, provides estimates of probability distributions of wavelet-transformed data being encoded. These estimates are utilized by an entropy coding subalgorithm that is another major component of the compression algorithm. The estimates make it possible to compress the image data more effectively than would otherwise be possible. The following background discussion is prerequisite to a meaningful summary of the context modeler. This discussion is presented relative to ICER-3D, which is the name attached to a particular compression algorithm and the software that implements it. The ICER-3D software is summarized briefly in the preceding article, ICER-3D Hyperspectral Image Compression Software (NPO-43238). Some aspects of this algorithm were previously described, in a slightly more general context than the ICER-3D software, in "Improving 3D Wavelet-Based Compression of Hyperspectral Images" (NPO-41381), NASA Tech Briefs, Vol. 33, No. 3 (March 2009), page 7a. In turn, ICER-3D is a product of generalization of ICER, another previously reported algorithm and computer program that can perform both lossless and lossy wavelet-based compression and decompression of gray-scale-image data. In ICER-3D, hyperspectral image data are decomposed using a 3D discrete wavelet transform (DWT). Following wavelet decomposition, mean values are subtracted from spatial planes of spatially low-pass subbands prior to encoding. The resulting data are converted to sign-magnitude form and compressed. In ICER-3D, compression is progressive, in that compressed information is ordered so that as more of the compressed data stream is received, successive reconstructions of the hyperspectral image data are of successively higher overall fidelity.

Kiely, Aaron↗

Pressure Sensitive Paint Applied to Flexible Models Project

One gap in current pressure-measurement technology is a high-spatial-resolution method for accurately measuring pressures on spatially and temporally varying wind-tunnel models such as Inflatable Aerodynamic Decelerators (IADs), parachutes, and sails. Conventional pressure taps only provide sparse measurements at discrete points and are difficult to integrate with the model structure without altering structural properties. Pressure Sensitive Paint (PSP) provides pressure measurements with high spatial resolution, but its use has been limited to rigid or semi-rigid models. Extending the use of PSP from rigid surfaces to flexible surfaces would allow direct, high-spatial-resolution measurements of the unsteady surface pressure distribution. Once developed, this new capability will be combined with existing stereo photogrammetry methods to simultaneously measure the shape of a dynamically deforming model in a wind tunnel. Presented here are the results and methodology for using PSP on flexible surfaces.

model deformation↗

Sampling Versus Filtering in Large-Eddy Simulations

A LES formalism in which the filter operator is replaced by a sampling operator is proposed. The unknown quantities that appear in the LES equations originate only from inadequate resolution (Discretization errors). The resulting viewpoint seems to make a link between finite difference approaches and finite element methods. Sampling operators are shown to commute with nonlinearities and to be purely projective. Moreover, their use allows an unambiguous definition of the LES numerical grid. The price to pay is that sampling never commutes with spatial derivatives and the commutation errors must be modeled. It is shown that models for the discretization errors may be treated using the dynamic procedure. Preliminary results, using the Smagorinsky model, are very encouraging.

Debliquy, O.↗

Dust from early galaxies and the X-ray background radiation

The effects of the scattering of cosmic X-ray background photons by dusty foreground galaxies are explored, assuming the X-rays themselves arise in discrete sources. The effective angle of scattering is about 1-10 arcmin in the 1-3 keV band, which can lead to a decrease in the spatial fluctuations of the background observed with an IPC-like imaging device with about 1 arcmin resolution. This can reduce the number of discrete sources required to explain the smoothness of the background to values greater than about 1000/sq deg. The effect may be stronger in the direction of distant clusters of galaxies, or if distant galaxies are much dustier than assumed here.

Rudak, B.↗

Evaluation of Galvanic Vestibular Stimulation System

Microgravity exposure results in an adaptive central reinterpretation of information from multiple sensory sources to produce a sensorimotor state appropriate for motor actions in this unique environment, but this new adaptive state is no longer appropriate for the 1-g gravitational environment on Earth. During these gravitational transitions, astronauts experience deficits in both perceptual and motor functions including impaired postural control, disruption in spatial orientation, impaired control of locomotion that include alterations in muscle activation variability, modified lower limb kinematics, alterations in head-trunk coordination as well as reduced dynamic visual acuity. Post-flight changes in postural and locomotor control might have adverse consequences if a rapid egress was required following a long-duration mission, where support personnel may not be available to aid crewmembers. The act of emergency egress includes, but is not limited to standing, walking, climbing a ladder, jumping down, monitoring displays, actuating discrete controls, operating auxiliary equipment, and communicating with Mission Control and recovery teams while maintaining spatial orientation, mobility and postural stability in order to escape safely. The average time to recover impaired postural control and functional mobility to preflight levels of performance has been shown to be approximately two weeks after long-duration spaceflight. The postflight alterations are due in part to central reinterpretation of vestibular information caused by exposure to microgravity. In this study we will use a commonly used technique of transcutaneous electrical stimulation applied across the vestibular end organs (galvanic vestibular stimulation, GVS) to disrupt vestibular function as a simulation of post-flight disturbances. The goal of this project is an engineering human-in-the-loop evaluation of a device that can degrade performance of functional tasks (e.g. to maintain upright balance) similar to what astronauts experience during transitions to new gravitational environments. Stochastic electrical stimulation can be applied to the vestibular system through electrodes placed over the mastoid process behind the ears in the binaural configuration resulting in stimulation in the mediolateral (side-to-side) plane. An additional electrode can be placed over the bony landmark of the tip of the c7 spinous process for the double monaural configuration, which will cause stimulation in the anteroposterior (forward-backward) plane. A portable constant current bipolar stimulator with subject isolation was designed and built to deliver the stimulus. The unit is powered using a 3.7 V battery pack and designed to produce currents up to 5 mA. The stimulator, controlled by a Raspberry Pi 3 computer, offers several stimulus signal generation options including a standalone mode, which uses onboard signal files stored on the flash memory card. Stochastic stimulation signals will be generated in 0-30 Hz frequency bandwidth. Stimulation amplitude can be increased incrementally to a maximum amplitude of 5.0 mA (e.g., 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0 mA). In control trials, subjects will be experiencing vestibular stimulation with 0-mA current applied through the electrodes. The system will be evaluated at various levels of stimulation and in both the binaural and double monaural electrode configurations. One of the objectives is to identify stimulation levels producing effects most comparable to the post-flight disturbances. This is a pilot study that will set the stage for a larger, more comprehensive study that will investigate wider aspects of post-flight sensorimotor dysfunction and set sensorimotor standards for crew health.

Kofman, I. S.↗

Application of the method of integral relations to unsteady fluid flow problems with shocks

A mixed method using both the method of integral relations and the finite difference technique is developed for the solution of unsteady flow problems. The integral relations method is based upon a chosen interpolating function dependent only on the time domain. The resulting local semidiscrete finite element equations obtained are assembled into a global form. The spatial derivatives at the nodes are replaced by finite difference operators and the discretized nonlinear algebraic system is solved by an iterative scheme. Solutions are obtained for the one-dimensional gasdynamics equation, the one-dimensional wave equation and Burger's model of turbulence. Agreement with other numerical and analytical solutions is excellent for the gasdynamics problem and satisfactory for the Burger equation in cases of small viscous effects.

Wadia, A. R.↗

Initialization of mesoscale models - The possible impact of remotely sensed data

Little or no improvement appears to have been achieved in mesoscale numerical prediction. This has been blamed on inherent unpredictability, insufficient spatial resolution, inadequate parameterization of important physical processes, and errors in the numerical discretization of the governing hydrodynamic equations. Attention is presently given to the more fundamental difficulty of forecasting when present conditions are not known with sufficient accuracy, and to the matching of remotely sensed sounding-capability satellite and all-weather Doppler radar data with mesoscale numerical models. This matching is not straightforward; simple and variational forms of four-dimensional assimilation, and Kalman filtering, are suggested approaches but their practical testing remains to be accomplished. The most severe problem faced is that of the temperature and humidity profiles' vertical resolution.

Gal-Chen, T.↗

Estimation of coefficients and boundary parameters in hyperbolic systems

Semi-discrete Galerkin approximation schemes in connection with inverse problems for the estimation of spatially varying coefficients and boundary condition parameters in second order hyperbolic systems typical of those arising in 1-D surface seismic problems are considered. Spline based algorithms are proposed for which theoretical convergence results along with a representative sample of numerical findings are given.

Banks, H. T.↗

Numerical methods for hypersonic boundary layer stability

Four different schemes for solving compressible boundary layer stability equations are developed and compared, considering both the temporal and spatial stability for a global eigenvalue spectrum and a local eigenvalue search. The discretizations considered encompass: (1) a second-order-staggered finite-difference scheme; (2) a fourth-order accurate, two-point compact scheme; (3) a single-domain Chebychev spectral collocation scheme; and (4) a multidomain spectral collocation scheme. As Mach number increases, the performance of the single-domain collocation scheme deteriorates due to the outward movement of the critical layer; a multidomain spectral method is accordingly designed to furnish superior resolution of the critical layer.

Malik, M. R.↗

Design of optimally smoothing multistage schemes for the Euler equations

A recently derived local preconditioning of the Euler equations is shown to be useful in developing multistage schemes suited for multigrid use. The effect of the preconditioning matrix on the spatial Euler operator is to equalize the characteristic speeds. When applied to the discretized Euler equations, the preconditioning has the effect of strongly clustering the operator's eigenvalues in the complex plane. This makes possible the development of explicit marching schemes that effectively damp most high-frequency Fourier modes, as desired in multigrid applications. The technique is the same as developed earlier for scalar convection schemes: placement of the zeros of the amplification factor of the multistage scheme in locations where eigenvalues corresponding to high-frequency modes abound.

Van Leer, Bram↗

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗

A Review of High-Order and Optimized Finite-Difference Methods for Simulating Linear Wave Phenomena

This paper presents a review of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, or elastic waves. The spatial operators reviewed include compact schemes, non-compact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods discussed include Runge-Kutta methods, Adams-Bashforth methods, and the leapfrog method. In addition, the following fourth-order fully-discrete finite-difference methods are considered: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method studied, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. Recommendations are made with respect to the suitability of the methods for specific problems and practical aspects of their use, such as appropriate Courant numbers and grid densities. Avenues for future research are suggested.

Zingg, David W.↗

Spatial resolution and cloud optical thickness retrievals

The impact of sensor spatial resolution and accurate cloud pixel identification on cloud property retrievals was investigated. Twelve fair weather cumulus (FWC) scenes of high spectral and spatial resolution Airborne Visible and Infrared Imaging Spectrometer (AVIRIS) data are analyzed. A variation of the 3-band ratio technique of Gao and Goetz is used to discriminate clouds from the background, and then a discrete ordinate radiative transfer model is used to obtain optical thickness of cloudy regions for each scene. To study the effect of spatial resolution upon retrieved optical thickness, the 20 m AVIRIS data was spatially degraded to spatial resolutions ranging from 40 to 960 m. Cloud area, scene average optical thickness, and distribution of retrieved optical thickness are determined at each spatial resolution. Finally, a comparison between the 3-band ratio technique and monospectral reflectance thresholding, using 20 m spatial resolution data, is presented.

Feind, Rand E.↗

Hybrid analog/digital, large format, photon counting detectors for astronomy

The development of a new microchannel plate photon-counting detector with an analog readout method based on a resistive anode is reported. This detector exhibits extremely high, stable electron gains of 10 to the 8th. At this gain, the spatial resolution is no longer primarily limited by the noise of the resistive anode, so that digital methods of readout, such as discrete conductors, lose their advantage. These detectors can be readily scaled to 40 mm and 70 mm formats to match plate scales of 2-m (and larger) telescopes. New, high speed digital electronics fully exploit the high spatial and time resolution made possible by gains of this level. Analysis of the theoretical performance of this detector shows that the major limitation to the spatial resolution is the proximity focus of the photocathode and the first microchannel plate. The detector has been mated to an echelle spectrograph developed.

Crocker, J.↗