Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multidimensional”

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 1,369 records · Page 76

Optimal band selection for dimensionality reduction of hyperspectral imagery

Hyperspectral images have many bands requiring significant computational power for machine interpretation. During image pre-processing, regions of interest that warrant full examination need to be identified quickly. One technique for speeding up the processing is to use only a small subset of bands to determine the 'interesting' regions. The problem addressed here is how to determine the fewest bands required to achieve a specified performance goal for pixel classification. The band selection problem has been addressed previously Chen et al., Ghassemian et al., Henderson et al., and Kim et al.. Some popular techniques for reducing the dimensionality of a feature space, such as principal components analysis, reduce dimensionality by computing new features that are linear combinations of the original features. However, such approaches require measuring and processing all the available bands before the dimensionality is reduced. Our approach, adapted from previous multidimensional signal analysis research, is simpler and achieves dimensionality reduction by selecting bands. Feature selection algorithms are used to determine which combination of bands has the lowest probability of pixel misclassification. Two elements required by this approach are a choice of objective function and a choice of search strategy.

Stearns, Stephen D.↗

Optimizing tertiary storage organization and access for spatio-temporal datasets

We address in this paper data management techniques for efficiently retrieving requested subsets of large datasets stored on mass storage devices. This problem represents a major bottleneck that can negate the benefits of fast networks, because the time to access a subset from a large dataset stored on a mass storage system is much greater that the time to transmit that subset over a network. This paper focuses on very large spatial and temporal datasets generated by simulation programs in the area of climate modeling, but the techniques developed can be applied to other applications that deal with large multidimensional datasets. The main requirement we have addressed is the efficient access of subsets of information contained within much larger datasets, for the purpose of analysis and interactive visualization. We have developed data partitioning techniques that partition datasets into 'clusters' based on analysis of data access patterns and storage device characteristics. The goal is to minimize the number of clusters read from mass storage systems when subsets are requested. We emphasize in this paper proposed enhancements to current storage server protocols to permit control over physical placement of data on storage devices. We also discuss in some detail the aspects of the interface between the application programs and the mass storage system, as well as a workbench to help scientists to design the best reorganization of a dataset for anticipated access patterns.

Chen, Ling Tony↗

VLSI neuroprocessors

Electronic and optoelectronic hardware implementations of highly parallel computing architectures address several ill-defined and/or computation-intensive problems not easily solved by conventional computing techniques. The concurrent processing architectures developed are derived from a variety of advanced computing paradigms including neural network models, fuzzy logic, and cellular automata. Hardware implementation technologies range from state-of-the-art digital/analog custom-VLSI to advanced optoelectronic devices such as computer-generated holograms and e-beam fabricated Dammann gratings. JPL's concurrent processing devices group has developed a broad technology base in hardware implementable parallel algorithms, low-power and high-speed VLSI designs and building block VLSI chips, leading to application-specific high-performance embeddable processors. Application areas include high throughput map-data classification using feedforward neural networks, terrain based tactical movement planner using cellular automata, resource optimization (weapon-target assignment) using a multidimensional feedback network with lateral inhibition, and classification of rocks using an inner-product scheme on thematic mapper data. In addition to addressing specific functional needs of DOD and NASA, the JPL-developed concurrent processing device technology is also being customized for a variety of commercial applications (in collaboration with industrial partners), and is being transferred to U.S. industries. This viewgraph p resentation focuses on two application-specific processors which solve the computation intensive tasks of resource allocation (weapon-target assignment) and terrain based tactical movement planning using two extremely different topologies. Resource allocation is implemented as an asynchronous analog competitive assignment architecture inspired by the Hopfield network. Hardware realization leads to a two to four order of magnitude speed-up over conventional techniques and enables multiple assignments, (many to many), not achievable with standard statistical approaches. Tactical movement planning (finding the best path from A to B) is accomplished with a digital two-dimensional concurrent processor array. By exploiting the natural parallel decomposition of the problem in silicon, a four order of magnitude speed-up over optimized software approaches has been demonstrated.

Kemeny, Sabrina E.↗

Application of a new finite difference algorithm for computational aeroacoustics

Acoustic problems have become extremely important in recent years because of research efforts such as the High Speed Civil Transport program. Computational aeroacoustics (CAA) requires a faithful representation of wave propagation over long distances, and needs algorithms that are accurate and boundary conditions that are unobtrusive. This paper applies a new finite difference method and boundary algorithm to the Linearized Euler Equations (LEE). The results demonstrate the ability of a new fourth order propagation algorithm to accurately simulate the genuinely multidimensional wave dynamics of acoustic propagation in two space dimensions with the LEE. The results also show the ability of a new outflow boundary condition and fourth order algorithm to pass the evolving solution from the computational domain with no perceptible degradation of the solution remaining within the domain.

Goodrich, John W.↗

Bayes classification of interferometric TOPSAR data

We report the Bayes classification of terrain types at different sites using airborne interferometric synthetic aperture radar (INSAR) data. A Gaussian maximum likelihood classifier was applied on multidimensional observations derived from the SAR intensity, the terrain elevation model, and the magnitude of the interferometric correlation. Training sets for forested, urban, agricultural, or bare areas were obtained either by selecting samples with known ground truth, or by k-means clustering of random sets of samples uniformly distributed across all sites, and subsequent assignments of these clusters using ground truth. The accuracy of the classifier was used to optimize the discriminating efficiency of the set of features that was chosen. The most important features include the SAR intensity, a canopy penetration depth model, and the terrain slope. We demonstrate the classifier's performance across sites using a unique set of training classes for the four main terrain categories. The scenes examined include San Francisco (CA) (predominantly urban and water), Mount Adams (WA) (forested with clear cuts), Pasadena (CA) (urban with mountains), and Antioch Hills (CA) (water, swamps, fields). Issues related to the effects of image calibration and the robustness of the classification to calibration errors are explored. The relative performance of single polarization Interferometric data classification is contrasted against classification schemes based on polarimetric SAR data.

Michel, T. R.↗

Nonlinear hydrodynamics of cosmological sheets. 1: Numerical techniques and tests

We present the numerical techniques and tests used to construct and validate a computer code designed to study the multidimensional nonlinear hydrodynamics of large-scale sheet structures in the universe, especially the fragmentation of such structures under various instabilities. This code is composed of two codes, the hydrodynamical code ZEUS-2D and a particle-mesh code. The ZEUS-2D code solves the hydrodynamical equations in two dimensions using explicit Eulerian finite-difference techniques, with modifications made to incorporate the expansion of the universe and the gas cooling due to Compton scattering, bremsstrahlung, and hydrogen and helium cooling. The particle-mesh code solves the equation of motion for the collisionless dark matter. The code uses two-dimensional Cartesian coordinates with a nonuniform grid in one direction to provide high resolution for the sheet structures. A series of one-dimensional and two-dimensional linear perturbation tests are presented which are designed to test the hydro solver and the Poisson solver with and without the expansion of the universe. We also present a radiative shock wave test which is designed to ensure the code's capability to handle radiative cooling properly. And finally a series of one-dimensional Zel'dovich pancake tests used to test the dark matter code and the hydro solver in the nonlinear regime are discussed and compared with the results of Bond et al. (1984) and Shapiro & Struck-Marcell (1985). Overall, the code is shown to produce accurate and stable results, which provide us a powerful tool to further our studies.

Anninos, Wenbo Y.↗

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.↗

Microphysical modeling of cirrus. 1: Comparison with 1986 FIRE IFO measurements

We have used a one-dimensional model of cirrus formation to study the development of cirrus clouds during the 1986 First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) intensive field observations (IFO). The cirrus model includes microphysical, dynamical, and radiative processes. Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing nucleation, growth by water vapor deposition, evaporation, coagulation, and vertical transport. We have focused on the cirrus observed on November 1, 1986. Vertical wind speed for the one-dimensional simulation is taken from a mesoscale model simulation for the appropriate time period. The mesoscale model simulation suggested that strong upward motions over Wyoming and subsequent horizontal transport of upper level moisture were responsible for the cirrus observed over Wisconsin on this date. We assumed that our one-dimensional model could be used to represent a vertical column moving from Wyoming to Wisconsin over a period of several hours. Ice crystal nucleation occurs in our model in the 8 to 10-km region as a result of the strong updrafts (and cooling) early in the simulation. Growth, coagulation, and sedimentation of these ice crystals result in a broad cloud region (5-10 km thick) with an optical depth of 1-2 after a few hours, in agreement with the FIRE measurements. Comparison with aircraft microphysical measurements made over Wisconsin indicates that the simulation generated reasonable ice water content, but the predicted ice number densities are too low, especially for radii less than about 50 microns. Sensitivity tests suggest that better agreement between simulated and observed microphysical properties is achieved if the nucleation rate is higher or stronger vertical mixing (perhaps associated with multidimensional motions) is present.

Jensen, Eric J.↗

Model study of atmospheric transport using carbon 14 and strontium 90 as inert tracers

The observed excess carbon 14 in the atmosphere from 1963 to 1970 provides unique, but limited, data up to an altitude of about 35 km for testing the air motions calculated by 11 multidimensional atmospheric models. Strontium 90 measurements in the atmosphere from 1964 to mid-1967 provide data that have more latitude coverage than those of carbon 14 and are useful for testing combined models of air motions and aerosol settling. Model calculations for carbon 14 begin at October 1963, 9 months after the conclusion of the nuclear bomb tests; the initial conditions for the calculations are derived by three methods, each of which agrees fairly well with measured carbon 14 in October 1963 and each of which has widely different values in regions of the stratosphere where there were no carbon 14 measurements. The model results are compared to the stratospheric measurements, not as if the observed data were absolute standards, but in an effort to obtain new insight about the models and about the atmosphere. The measured carbon 14 vertical profiles at 31 deg N are qualitatively different from all of the models; the measured vertical profiles show a maximum mixing ratio in the altitude range of 20 to 25 km from October 1963 through July 1966, but all modeled profiles show mixing ratio maxima that increase in altitude from 20 km in October 1963 to greater than 40 km by April 1966. Both carbon 14 and strontium 90 data indicate that the models differ substantially among themselves with respect to stratosphere-troposphere exchange rate, but the modeled carbon 14 stratospheric residence times indicate that differences among the models are small with respect to transport rate between the middle stratosphere and the lower stratosphere. Strontium 90 data indicate that aerosol settling is important up to at least 35 km altitude. Relative to the measurements, about three quarters of the models transport carbon 14 from the lower stratosphere to the troposphere too rapidly, and all models with all three sets of initial conditions appear to sweep carbon 14 out of the midstratosphere (above 28 km) much more slowly than was observed from 1963 to 1970. Atmospheric dynamicists should seriously consider that the carbon 14 measurements imply with respect to advection in the midstratosphere.

Kinnison, D. E.↗

Equilibrium structure of solar magnetic flux tubes: Energy transport with multistream radiative transfer

We examine the equilibrium structure of vertical intense magnetic flux tubes on the Sun. Assuming cylindrical geometry, we solve the magnetohydrostatic equations in the thin flux-tube approximation, allowing for energy transport by radiation and convection. The radiative transfer equation is solved in the six-stream approximation, assuming gray opacity and local thermodynamic equilibrium. This constitutes a significant improvement over a previous study, in which the transfer was solved using the multidimensional generalization of the Eddington approximation. Convection in the flux tube is treated using mixing-length theory, with an additional parameter alpha, characterizing the suppression of convective energy transport in the tube by the strong magnetic field. The equations are solved using the method of partial linearization. We present results for tubes with different values of the magnetic field strength and radius at a fixed depth in the atmosphere. In general, we find that, at equal geometric heights, the temperature on the tube axis, compared to the ambient medium, is higher in the photosphere and lower in the convection zone, with the difference becoming larger for thicker tubes. At equal optical depths the tubes are generally hotter than their surroundings. The results are comparatively insensitive to alpha but depend upon whether radiative and convective energy transport operate simultaneously or in separate layers. A comparison of our results with semiempirical models shows that the temperature and intensity contrast are in broad agreement. However, the field strengths of the flux-tube models are somewhat lower than the values inferred from observations.

Hasan, S. S.↗

Target compressions by working fluids driven with solid liner implosions

Compression by a spherical solid liner of a gold target surrounded by a hydrogen plasma is simulated. Two-dimensional simulations that treat only a subset of the physics included in the one-dimensional code were performed in an attempt to assess multidimensional effects. A one-dimensional numerical code has been developed to study the effects of thermal radiation and conduction. Results of pressure, density, and energy deposited for different initial plasma conditions are presented and discussed. Results from both one- and two-dimensional codes show that the average target density at peak compression is 39-43 g/cu cm, using the SHIVA Star facility at 90 kV discharge.

Chiang, P.-R.↗

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

Parallel computing using a Lagrangian formulation

This paper adopts a new Lagrangian formulation of the Euler equation for the calculation of two dimensional supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, we have achieved better than six times speed-up on a 8192-processor CM-2 over a single processor of a CRAY-2.

Liou, May-Fun↗

Algorithmic trends in computational fluid dynamics; The Institute for Computer Applications in Science and Engineering (ICASE)/LaRC Workshop, NASA Langley Research Center, Hampton, VA, US, Sep. 15-17, 1991

The purpose here is to assess the state of the art in the areas of numerical analysis that are particularly relevant to computational fluid dynamics (CFD), to identify promising new developments in various areas of numerical analysis that will impact CFD, and to establish a long-term perspective focusing on opportunities and needs. Overviews are given of discretization schemes, computational fluid dynamics, algorithmic trends in CFD for aerospace flow field calculations, simulation of compressible viscous flow, and massively parallel computation. Also discussed are accerelation methods, spectral and high-order methods, multi-resolution and subcell resolution schemes, and inherently multidimensional schemes.

Hussaini, M. Y.↗

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗

Advanced compilation techniques in the PARADIGM compiler for distributed-memory multicomputers

The PARADIGM compiler project provides an automated means to parallelize programs, written in a serial programming model, for efficient execution on distributed-memory multicomputers. .A previous implementation of the compiler based on the PTD representation allowed symbolic array sizes, affine loop bounds and array subscripts, and variable number of processors, provided that arrays were single or multi-dimensionally block distributed. The techniques presented here extend the compiler to also accept multidimensional cyclic and block-cyclic distributions within a uniform symbolic framework. These extensions demand more sophisticated symbolic manipulation capabilities. A novel aspect of our approach is to meet this demand by interfacing PARADIGM with a powerful off-the-shelf symbolic package, Mathematica. This paper describes some of the Mathematica routines that performs various transformations, shows how they are invoked and used by the compiler to overcome the new challenges, and presents experimental results for code involving cyclic and block-cyclic arrays as evidence of the feasibility of the approach.

Su, Ernesto↗

Selecting a general-purpose data compression algorithm

The National Space Science Data Center's Common Data Formate (CDF) is capable of storing many types of data such as scalar data items, vectors, and multidimensional arrays of bytes, integers, or floating point values. However, regardless of the dimensionality and data type, the data break down into a sequence of bytes that can be fed into a data compression function to reduce the amount of data without losing data integrity and thus remaining fully reconstructible. Because of the diversity of data types and high performance speed requirements, a general-purpose, fast, simple data compression algorithm is required to incorporate data compression into CDF. The questions to ask are how to evaluate and compare compression algorithms, and what compression algorithm meets all requirements. The object of this paper is to address these questions and determine the most appropriate compression algorithm to use within the CDF data management package that would be applicable to other software packages with similar data compression needs.

Mathews, Gary Jason↗

Aerodynamic, unsteady, kinetic and heat loss effects on the dynamics and structure of weakly-burning flames

The first objective of the program is to introduce the meritorious counterflow methodology in microgravity in order to quantify the steady and unsteady characteristics of weakly-burning premixed and diffusion flames for a wide variety of conditions including elevated pressures. Subsequently, through detailed modeling and comparisons with the experimental data, to provide physical insight into the elementary mechanisms controlling the flame response. The configuration offers good control over the parameters of interest and can be modelled closely. The knowledge which will be gained from the counterflow flames will be subsequently used to analyze near-limit phenomena related to other configurations by conducting detailed numerical simulations including multidimensional ones. Among the problems to be analyzed are the downward and upward propagation of near-limit flames in tubes and phenomena observed in spherical and cylindrical geometries.

Egolfopoulos, Fokion N.↗