Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical algorithm computing”

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 865 records · Page 48

Flightspeed Integral Image Analysis Toolkit

The Flightspeed Integral Image Analysis Toolkit (FIIAT) is a C library that provides image analysis functions in a single, portable package. It provides basic low-level filtering, texture analysis, and subwindow descriptor for applications dealing with image interpretation and object recognition. Designed with spaceflight in mind, it addresses: Ease of integration (minimal external dependencies) Fast, real-time operation using integer arithmetic where possible (useful for platforms lacking a dedicated floatingpoint processor) Written entirely in C (easily modified) Mostly static memory allocation 8-bit image data The basic goal of the FIIAT library is to compute meaningful numerical descriptors for images or rectangular image regions. These n-vectors can then be used directly for novelty detection or pattern recognition, or as a feature space for higher-level pattern recognition tasks. The library provides routines for leveraging training data to derive descriptors that are most useful for a specific data set. Its runtime algorithms exploit a structure known as the "integral image." This is a caching method that permits fast summation of values within rectangular regions of an image. This integral frame facilitates a wide range of fast image-processing functions. This toolkit has applicability to a wide range of autonomous image analysis tasks in the space-flight domain, including novelty detection, object and scene classification, target detection for autonomous instrument placement, and science analysis of geomorphology. It makes real-time texture and pattern recognition possible for platforms with severe computational restraints. The software provides an order of magnitude speed increase over alternative software libraries currently in use by the research community. FIIAT can commercially support intelligent video cameras used in intelligent surveillance. It is also useful for object recognition by robots or other autonomous vehicles

Thompson, David R.↗

Computation of viscous transonic flow over porous airfoils

The viscous effects on transonic flow past an airfoil which contains a shallow cavity beneath a porous surface are studied numerically. The porous region occupies a small portion of the total airfoil surface, and is located near the shock. Both an interactive boundary layer (IBL) algorithm and a thin-layer Navier-Stokes (TLNS) algorithm have been modified for use in studying the outer flow, whereas a stream-function formulation has been used to model the inner flow in the small cavity. The coupling procedure at the porous surface is based on Darcy's law and on the assumption of a constant total presusre in the cavity. In addition, a modified Baldwin-Lomax turbulence model is used to consider the transpired turbulent boundary layer in the TLNS approach, and the Cebeci-Smith turbulence model is used in the IBL approach. According to the present analysis, a porous surface can reduce the wave drag appreciably, but it can also increase viscous losses. As has been observed experimentally, the numerical results indicate that the total drag is reduced at higher Mach numbers and increased at lower Mach numbers when the angles of attack are small. Furthermore, the streamline patterns of passive-shock and boundary-layer interaction are revealed in this study.

Chow, Chuen-Yen↗

A numerical investigation of the influence of heating on the excitation of high and low Mach number jet flows

The influence of heating and velocity on the unsteady nature of jet mixing layers is investigated by solving the time-dependent compressible Navier-Stokes equations using MacCormack's explicit finite difference algorithm. The computations are performed for jet Mach numbers of 0.3 and 0.8 with flow total temperatures up to 800 K. The Reynolds number ranges from 3 x 10 to the 5th to 1.3 x 10 to the 6th. Excitation is accomplished by imposing an acoustic pressure signal inside the jet duct. The objective of this effort is to compare the response of heated and unheated jets with acoustic excitation at high and low subsonic Mach numbers. The preliminary results indicate that without acoustic excitation the mixing in the heated jet is greater than in the unheated case. It has also been found that for the low speed heated jets (M = 0.3, Tt = 672K) acoustic excitation causes the production of large scale vortex structures to occur in a regular periodic manner with organized periodic pairing at some Strouhal numbers. These results are in agreement with experimental data.

Scott, J. N.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool.

Chima, R. V.↗

Numerical predictions of residence times behind a rearward facing step with transverse injection

A numerical investigation of the two-dimensional recirculation zone formed between a rearward facing step and a transverse hydrogen jet in a supersonic flow is presented. Calculations were made using an elliptic flow computer code with a time-relaxation finite difference integration algorithm. Cases were examined for jet static pressures between 2 and 5 times mainstream with jets located 3 to 7 step heights downstream of the step. Nominal mainstream flow conditions were Mach 2.4 air at 1 atm pressure and 1000 K (1800 R) static temperature. Results in the form of mean properties in the recirculation zone and residence time are presented. Calculations of the decay of hydrogen mass from the steady state are presented for selected cases. The combustion potential for each case is assessed by comparison of computed conditions in the recirculation zone with the results from a well-stirred reactor theory.

Rogers, R. C.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool. Previously announced in STAR as N84-22527

Chima, R. V.↗

The effect of Mach number on the stability of a plane supersonic wave

The influence of compressibility on the mechanisms governing the various stages of transition in a supersonic wake is investigated. Results from linear stability theory are used to provide physical insights into the observed reduction in growth rate at high Mach numbers. A newly developed hybrid algorithm is used to solve the compressible inviscid linear disturbance equations. Growth rates for both antisymmetric and symmetric modes of two-dimensional and oblique waves are computed for a wide range of Mach numbers. Results from two and three-dimensional direct numerical simulations of a forced compressible time-developing wake are presented in order to understand the nonlinear stages of transition at high Mach numbers. Observed nonlinear growth rate comparisons are made for wakes at two different Mach numbers. The reduction in growth rate at high Mach numbers is explained by examining contour plots of baroclinic torques and the product of dilatation and vorticity.

Chen, Jacqueline H.↗

Computational fluid dynamics - A personal view

This paper provides a personal view of computational fluid dynamics. The main theme is divided into two categories - one dealing with algorithms and engineering applications and the other with scientific investigations. The former category may be termed computational aerodynamics, with the objective of providing reliable aerodynamic or engineering predictions. The latter category is essentially basic research, where the algorithmic tools are used to unravel and elucidate fluid-dynamic phenomena hard to obtain in a laboratory. A critique of the numerical solution techniques for both compressible and incompressible flows is included. The discussion on scientific investigations deals in particular with transition and turbulence.

Hussaini, M. Y.↗

Incompressible viscous flow computations for the pump components and the artificial heart

A finite-difference, three-dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. In this work, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two. Included in the appendix is a paper on incompressible viscous flow through artificial heart devices with moving boundaries. Time-accurate calculations, such as impeller and diffusor interaction, will be reported in future work.

Kiris, Cetin↗

Use of non-spherical gravity harmonics for relative motion GN/C (category 2)

Detailed analysis of the Automatic Rendezvous and Capture problem indicate a need for three different regions of mathematical description for the GN&C algorithms: (1) multi-vehicle orbital mechanics to the rendezvous interface point, i.e., within 100 nm; (2) relative motion solutions (such as Clohessy-Wiltshire type) from the far-field to the near-field interface, i.e., within 1 nm and; (3) close proximity motion - the near-field motion where the relative differences in the gravitational and orbit inertial accelerations can be neglected from the equations of motion. Limit boundaries to these regions can be precisely defined by further analysis and will be functions of the tracking measurement accuracies and the computer resources available for the solution of the algorithms. This paper analyzes the relative motion in Regions 2 and 3 above and present the derivation and discussion of the general case of non-spherical gravitational perturbed relative motion. Mathematical deviations from the numerically integrated spherical gravity case and solutions from the Clohessy-Wiltshire equations are presented in the analysis. Based upon this preliminary analysis, it is recommended that further efforts be used to assess the relative position and velocity differences in Region 2 due to non-spherical gravity harmonics and that viable GN&C algorithms be developed to include these gravity perturbations (especially the effects of the first gravity harmonic, J2).

Henderson, David M.↗

Parallel methods for dynamic simulation of multiple manipulator systems

In this paper, efficient dynamic simulation algorithms for a system of m manipulators, cooperating to manipulate a large load, are developed; their performance, using two possible forms of parallelism on a general-purpose parallel computer, is investigated. One form, temporal parallelism, is obtained with the use of parallel numerical integration methods. A speedup of 3.78 on four processors of CRAY Y-MP8 was achieved with a parallel four-point block predictor-corrector method for the simulation of a four manipulator system. These multi-point methods suffer from reduced accuracy, and when comparing these runs with a serial integration method, the speedup can be as low as 1.83 for simulations with the same accuracy. To regain the performance lost due to accuracy problems, a second form of parallelism is employed. Spatial parallelism allows most of the dynamics of each manipulator chain to be computed simultaneously. Used exclusively in the four processor case, this form of parallelism in conjunction with a serial integration method results in a speedup of 3.1 on four processors over the best serial method. In cases where there are either more processors available or fewer chains in the system, the multi-point parallel integration methods are still advantageous despite the reduced accuracy because both forms of parallelism can then combine to generate more parallel tasks and achieve greater effective speedups. This paper also includes results for these cases.

Mcmillan, Scott↗

Eigenvalue routines in NASTRAN: A comparison with the Block Lanczos method

The NASA STRuctural ANalysis (NASTRAN) program is one of the most extensively used engineering applications software in the world. It contains a wealth of matrix operations and numerical solution techniques, and they were used to construct efficient eigenvalue routines. The purpose of this paper is to examine the current eigenvalue routines in NASTRAN and to make efficiency comparisons with a more recent implementation of the Block Lanczos algorithm by Boeing Computer Services (BCS). This eigenvalue routine is now available in the BCS mathematics library as well as in several commercial versions of NASTRAN. In addition, CRAY maintains a modified version of this routine on their network. Several example problems, with a varying number of degrees of freedom, were selected primarily for efficiency bench-marking. Accuracy is not an issue, because they all gave comparable results. The Block Lanczos algorithm was found to be extremely efficient, in particular, for very large size problems.

Tischler, V. A.↗

A generic efficient adaptive grid scheme for rocket propulsion modeling

The objective of this research is to develop an efficient, time-accurate numerical algorithm to discretize the Navier-Stokes equations for the predictions of internal one-, two-dimensional and axisymmetric flows. A generic, efficient, elliptic adaptive grid generator is implicitly coupled with the Lower-Upper factorization scheme in the development of ALUNS computer code. The calculations of one-dimensional shock tube wave propagation and two-dimensional shock wave capture, wave-wave interactions, shock wave-boundary interactions show that the developed scheme is stable, accurate and extremely robust. The adaptive grid generator produced a very favorable grid network by a grid speed technique. This generic adaptive grid generator is also applied in the PARC and FDNS codes and the computational results for solid rocket nozzle flowfield and crystal growth modeling by those codes will be presented in the conference, too. This research work is being supported by NASA/MSFC.

Mo, J. D.↗

Issues in the analysis of oligonucleotide tiling microarrays for transcript mapping

Traditional microarrays use probes complementary to known genes to quantitate the differential gene expression between two or more conditions. Genomic tiling microarray experiments differ in that probes that span a genomic region at regular intervals are used to detect the presence or absence of transcription. This difference means the same sets of biases and the methods for addressing them are unlikely to be relevant to both types of experiment. We introduce the informatics challenges arising in the analysis of tiling microarray experiments as open problems to the scientific community and present initial approaches for the analysis of this nascent technology.

Review↗

Optimum propellant usage for reaction jet systems of space vehicles

The on-off type control for reaction jet systems is proven to be the optimal fuel scheme. However, due to the nonlinear characteristics of this type control, no direct method of solution is known for this optimal process being applied to the attitude control of space vehicles. This paper will discuss the optimization of fuel usage for an attitude control of space vehicles. A computation technique is developed for the calculation of optimal control law. The method of calculus of variations is applied to the estimation of the changes of performance index as well as terminal constraints. Thus an algorithm is obtained by the steepest descent method. A numerical example is given in the paper.

Liu, T. C.↗

Inventory estimation on the massively parallel processor

This paper describes algorithms for efficiently computing inventory estimates from satellite based images. The algorithms incorporate a one dimensional feature extraction which optimizes the pairwise sum of Fisher distances. Biases are eliminated with a premultiplication by the inverse of the analytically derived error matrix. The technique is demonstrated with a numerical example using statistics obtained from an actual Landsat scene. Attention was given to implementation of the Massively Parallel processor (MPP). A timing analysis demonstrates that the inventory estimation can be performed an order of magnitude faster on the MPP than on a conventional serial machine.

Argentiero, P. D.↗

Numerical simulation of excited jet mixing layers

A numerical simulation of unsteady flow in jet mixing layers, both with and without external excitation, has been performed by solving the time-dependent compressible Navier-Stokes equations. Computations were performed on a CRAY X-MP computer using MacCormick's explicit finite difference algorithm. Different excitation methods were investigated and were shown to be very effective in controlling the well organized periodic production, shedding and pairing of large scale vortex structures. It is found that pressure excitation was generally more effective than temperature excitation, and that grid refinement results in substantial improvement in the resolution of unsteady features. The location and orientation, in addition to the frequency, of the excitation source are shown to have a significant influence on the production and interaction of large scale vortex structures in the jet mixing layer.

Scott, J. N.↗

A new approach to the solution of boundary value problems involving complex configurations

A new approach for solving certain types of boundary value problems about complex configurations is presented. Numerical algorithms from such diverse fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the memory, speed and architecture of current and emerging supercomputers. Although the approach has application to many branches of computational physics, the present effort is concentrated in areas of Computational Fluid Dynamics (CFD) such as steady nonlinear aerodynamics, time harmonic unsteady aerodynamics, and aeroacoustics. The most significant attribute of the approach is that it can handle truly arbitrary boundary geometries and eliminates the difficult task of generating surface fitted grids.

Rubbert, P. E.↗