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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

NASA NTRS1995

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

NASA NTRS1995

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↗

NASA NTRS1995

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↗

NASA NTRS1995

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

NASA NTRS1995

Estimation of polar stratospheric cloud infrared extinction climatology using visible satellite observations

Polar stratospheric clouds (PSC's) provide surfaces for heterogeneous processes which can dramatically alter the normal partitioning of odd nitrogen and chlorine families in the winter polar stratospheres, setting up conditions for significant ozone depletion as manifested in the springtime Antarctic ozone hole. The spatial and temporal distribution of PSC's is important for parameterizing PSC occurrence in multidimensional photochemical models whose use is essential for fully understanding observed Antarctic ozone losses as well as for accessing the possibility of a similar phemonenon occurring in the future in the Arctic. The Stratospheric Aerosol Measurement (SAM) 2 sensor, a single-channel (1mu m) photometer launched into a Sun-synchronous orbit aboard the Nimbus 7 satellite in October 1978, provided a unique database to establish the climatology of PSC's. Poole and Pitts (1994) used the record of high-latitude aerosol extinction obtained by SAM II from 1979-1989 to establish the climatology of PSC occurrences in the Arctic and Antarctic. Unfortunately, little information about PSC composition or type was detectable from the single-wavelength SAM II data.

Pitts, Michael C.↗

NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗

NASA NTRS1995

Visualization Techniques in Space and Atmospheric Sciences

Unprecedented volumes of data will be generated by research programs that investigate the Earth as a system and the origin of the universe, which will in turn require analysis and interpretation that will lead to meaningful scientific insight. Providing a widely distributed research community with the ability to access, manipulate, analyze, and visualize these complex, multidimensional data sets depends on a wide range of computer science and technology topics. Data storage and compression, data base management, computational methods and algorithms, artificial intelligence, telecommunications, and high-resolution display are just a few of the topics addressed. A unifying theme throughout the papers with regards to advanced data handling and visualization is the need for interactivity, speed, user-friendliness, and extensibility.

Szuszczewicz, E. P.↗

NASA NTRS1993

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗

NASA NTRS1996

The Kirchhoff Formulas for Moving Surfaces in Aeroacoustics - The Subsonic and Supersonic Cases

One of the active areas of computational aeroacoustics is the application of the Kirchhoff formulas to the problems of the rotating machinery noise predictions. The original Kirchhoff formula was derived for a stationary surface. In 1988, Farassat and Myers derived a Kirchhoff Formula obtained originally by Morgans using modem mathematics. These authors gave a formula particularly useful for applications in aeroacoustics. This formula is for a surface moving at subsonic speed. Later in 1995 these authors derived the Kirchhoff formula for a super-sonically moving surface. This technical memorandum presents the viewgraphs of a day long workshop by the author on the derivation of the Kirchhoff formulas. All necessary background mathematics such as differential geometry and multidimensional generalized function theory are discussed in these viewgraphs. Abstraction is kept at minimum level here. These viewgraphs are also suitable for understanding the derivation and obtaining the solutions of the Ffowcs Williams-Hawkings equation. In the first part of this memorandum, some introductory remarks are made on generalized functions, the derivation of the Kirchhoff formulas and the development and validation of Kirchhoff codes. Separate lists of references by Lyrintzis, Long, Strawn and their co-workers are given in this memorandum. This publication is aimed at graduate students, physicists and engineers who are in need of the understanding and applications of the Kirchhoff formulas in acoustics and electromagnetics.

Farassat, F.↗

NASA NTRS1996

A Genuinely Two-Dimensional Scheme for the Compressible Euler Equations

We present a new genuinely multidimensional discretization for the compressible Euler equations. It is the only high-resolution scheme known to us where Gauss-Seidel relaxation is stable when applied as a smoother directly to the resulting high-resolution scheme. This allows us to construct a very simple and highly efficient multigrid steady-state solver. The scheme is formulated on triangular (possibly unstructured) meshes.

Sidilkover, David↗

NASA NTRS1996

Accurate Finite Difference Algorithms

Two families of finite difference algorithms for computational aeroacoustics are presented and compared. All of the algorithms are single step explicit methods, they have the same order of accuracy in both space and time, with examples up to eleventh order, and they have multidimensional extensions. One of the algorithm families has spectral like high resolution. Propagation with high order and high resolution algorithms can produce accurate results after O(10(exp 6)) periods of propagation with eight grid points per wavelength.

Goodrich, John W.↗

NASA NTRS1997

A Model-Based Approach for the Measurement of Eye Movements Using Image Processing

This paper describes a video eye-tracking algorithm which searches for the best fit of the pupil modeled as a circular disk. The algorithm is robust to common image artifacts such as the droopy eyelids and light reflections while maintaining the measurement resolution available by the centroid algorithm. The presented algorithm is used to derive the pupil size and center coordinates, and can be combined with iris-tracking techniques to measure ocular torsion. A comparison search method of pupil candidates using pixel coordinate reference lookup tables optimizes the processing requirements for a least square fit of the circular disk model. This paper includes quantitative analyses and simulation results for the resolution and the robustness of the algorithm. The algorithm presented in this paper provides a platform for a noninvasive, multidimensional eye measurement system which can be used for clinical and research applications requiring the precise recording of eye movements in three-dimensional space.

Sung, Kwangjae↗

NASA NTRS1997

Verification and Improvement of Flamelet Approach for Non-Premixed Flames

Studies in the mathematical modeling of the high-speed turbulent combustion has received renewal attention in the recent years. The review of fundamentals, approaches and extensive bibliography was presented by Bray, Libbi and Williams. In order to obtain accurate predictions for turbulent combustible flows, the effects of turbulent fluctuations on the chemical source terms should be taken into account. The averaging of chemical source terms requires to utilize probability density function (PDF) model. There are two main approaches which are dominant in high-speed combustion modeling now. In the first approach, PDF form is assumed based on intuitia of modelliers (see, for example, Spiegler et.al.; Girimaji; Baurle et.al.). The second way is much more elaborate and it is based on the solution of evolution equation for PDF. This approach was proposed by S.Pope for incompressible flames. Recently, it was modified for modeling of compressible flames in studies of Farschi; Hsu; Hsu, Raji, Norris; Eifer, Kollman. But its realization in CFD is extremely expensive in computations due to large multidimensionality of PDF evolution equation (Baurle, Hsu, Hassan).

Zaitsev, S.↗

NASA NTRS1997

Dynamic Forms: Application to Aircraft Guidance - Part 2

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of waypoints through which the aircraft trajectory must pass. The waypoints typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory that satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multidimensional, multiaxis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions are required to be smooth. The proposed guidance algorithm is based on the inversion of the pure feedback approximation, followed by correction for the effects of zero dynamics. The paper describes the structure and major modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

NASA NTRS1998

Integrated Data Visualization and Virtual Reality Tool

The Integrated Data Visualization and Virtual Reality Tool (IDVVRT) Phase II effort was for the design and development of an innovative Data Visualization Environment Tool (DVET) for NASA engineers and scientists, enabling them to visualize complex multidimensional and multivariate data in a virtual environment. The objectives of the project were to: (1) demonstrate the transfer and manipulation of standard engineering data in a virtual world; (2) demonstrate the effects of design and changes using finite element analysis tools; and (3) determine the training and engineering design and analysis effectiveness of the visualization system.

Dryer, David A.↗

NASA NTRS1998

Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Frink, Neal T.↗

NASA NTRS1999

Performance and Application of Parallel OVERFLOW Codes on Distributed and Shared Memory Platforms

The presentation discusses recent studies on the performance of the two parallel versions of the aerodynamics CFD code, OVERFLOW_MPI and _MLP. Developed at NASA Ames, the serial version, OVERFLOW, is a multidimensional Navier-Stokes flow solver based on overset (Chimera) grid technology. The code has recently been parallelized in two ways. One is based on the explicit message-passing interface (MPI) across processors and uses the _MPI communication package. This approach is primarily suited for distributed memory systems and workstation clusters. The second, termed the multi-level parallel (MLP) method, is simple and uses shared memory for all communications. The _MLP code is suitable on distributed-shared memory systems. For both methods, the message passing takes place across the processors or processes at the advancement of each time step. This procedure is, in effect, the Chimera boundary conditions update, which is done in an explicit "Jacobi" style. In contrast, the update in the serial code is done in more of the "Gauss-Sidel" fashion. The programming efforts for the _MPI code is more complicated than for the _MLP code; the former requires modification of the outer and some inner shells of the serial code, whereas the latter focuses only on the outer shell of the code. The _MPI version offers a great deal of flexibility in distributing grid zones across a specified number of processors in order to achieve load balancing. The approach is capable of partitioning zones across multiple processors or sending each zone and/or cluster of several zones into a single processor. The message passing across the processors consists of Chimera boundary and/or an overlap of "halo" boundary points for each partitioned zone. The MLP version is a new coarse-grain parallel concept at the zonal and intra-zonal levels. A grouping strategy is used to distribute zones into several groups forming sub-processes which will run in parallel. The total volume of grid points in each group are approximately balanced. A proper number of threads are initially allocated to each group, and in subsequent iterations during the run-time, the number of threads are adjusted to achieve load balancing across the processes. Each process exploits the multitasking directives already established in Overflow.

Djomehri, M. Jahed↗

NASA NTRS1999

Inlet-Compressor Analysis using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. Isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (1-D) Euler inlet simulation showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition. To aid in this understanding and to help evaluate possible boundary conditions, an inlet-engine experiment was conducted at the University of Cincinnati. The interaction of acoustic pulses, generated in the inlet, with the engine were investigated. Because of the availability of experimental data for validation, it was decided to simulate the experiment using CFD. The philosophy here is that the inlet-engine system is best simulated by coupling (existing) specialized CFD component-codes. The objectives of this work were to aid in a better understanding of inlet-compressor interaction physics and the formulation of a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets. Previous simulations have used 1-D Euler engine simulations in conjunction with 1-D Euler and axisymmetric Euler inlet simulations. This effort is a first step toward CFD simulation of an entire engine by coupling multidimensional component codes.

Cole, Gary↗