Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Algebra”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A general multiblock Euler code for propulsion integration. Volume 1: Theory document

A general multiblock Euler solver was developed for the analysis of flow fields over geometrically complex configurations either in free air or in a wind tunnel. In this approach, the external space around a complex configuration was divided into a number of topologically simple blocks, so that surface-fitted grids and an efficient flow solution algorithm could be easily applied in each block. The computational grid in each block is generated using a combination of algebraic and elliptic methods. A grid generation/flow solver interface program was developed to facilitate the establishment of block-to-block relations and the boundary conditions for each block. The flow solver utilizes a finite volume formulation and an explicit time stepping scheme to solve the Euler equations. A multiblock version of the multigrid method was developed to accelerate the convergence of the calculations. The generality of the method was demonstrated through the analysis of two complex configurations at various flow conditions. Results were compared to available test data. Two accompanying volumes, user manuals for the preparation of multi-block grids (vol. 2) and for the Euler flow solver (vol. 3), provide information on input data format and program execution.

Chen, H. C.↗
NASA NTRS1990

Time-accurate Navier-Stokes computations of self-excited two-dimensional unsteady cavity flows

Time-accurate Navier Stokes computations were performed to study self-excited oscillatory shear layer behavior over a two-dimensional rectangular cavity and to gain insight into the vortex-acoustic coupling that sustains these oscillations. Results were obtained solving the full, unsteady compressible Navier-Stokes equations for Reynolds-averaged turbulent air flow past a cavity with L/d = 2. Conditions duplicated test conditions in experiments by Rossiter at Mach numbers varied from 0.6 to 1.2. The computational scheme was the USA 2D series code developed by Rockwell's Science Center with an algebraic turbulence model for the shear layer and recirculating region. Two physical conditions were simulated: an isolated cavity situated (1) in an infinite plane wall and (2) in an infinite channel with an opposing parallel wall. Computed results for (1) showed excellent agreement with Rossiter's measured amplitudes and frequencies. A tuned acoustic duct resonance effect was observed in (2), demonstrating a controlled change in the acoustic feedback mechanism and resultant amplification in computed oscillations.

Dougherty, N. Sam↗
NASA NTRS1982

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗
NASA NTRS1994

Boundary element analysis on vector and parallel computers

Boundary element analysis (BEA) can be characterized as a numerical technique that generally shifts the computational burden in the analysis toward numerical integration and the solution of nonsymmetric and either dense or blocked sparse systems of algebraic equations. Researchers have explored the concept that the fundamental characteristics of BEA can be exploited to generate effective implementations on vector and parallel computers. In this paper, the results of some of these investigations are discussed. The performance of overall algorithms for BEA on vector supercomputers, massively data parallel single instruction multiple data (SIMD), and relatively fine grained distributed memory multiple instruction multiple data (MIMD) computer systems is described. Some general trends and conclusions are discussed, along with indications of future developments that may prove fruitful in this regard.

Kane, J. H.↗
NASA NTRS1981

Two-boundary grid generation for the solution of the three dimensional compressible Navier-Stokes equations

A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.

Smith, R. E.↗
NASA NTRS1973

A rational approach to the use of Prandtl's mixing length model in free turbulent shear flow calculations

Prandtl's basic mixing length model was used to compute 22 test cases on free turbulent shear flows. The calculations employed appropriate algebraic length scale equations and single values of mixing length constant for planar and axisymmetric flows, respectively. Good agreement with data was obtained except for flows, such as supersonic free shear layers, where large sustained sensitivity changes occur. The inability to predict the more gradual mixing in these flows is tentatively ascribed to the presence of a significant turbulence-induced transverse static pressure gradient which is neglected in conventional solution procedures. Some type of an equation for length scale development was found to be necessary for successful computation of highly nonsimilar flow regions such as jet or wake development from thick wall flows.

Rudy, D. H.↗
NASA NTRS2004

Simplified computational approach for dual-probe heat-pulse method

Two equations are currently available for estimating soil volumetric heat capacity (pc) with the dual-probe heat-pulse (DPHP) method. One is simple but gives only approximate results because it assumes that the DPHP sensor releases an impulse of heat instantaneously. The other explicitly accounts for the finite duration of heating and gives exact results. Unfortunately, the equation that gives exact results involves the exponential integral function, which is not available in most computer spreadsheet software packages or data logger function libraries. In this note we introduce an approximation of the exact equation that contains only simple algebraic functions. The approximation consists of the first five terms of a Taylor series, which are written as a telescoped polynomial for computational purposes. For most applications of the DPHP method, the polynomial approximation gives estimates of pc that are at least an order of magnitude more accurate than estimates obtained from the simple equation based on instantaneous heating.

NASA Center JSC↗
NASA NTRS1992

Computation of vortex wake flows and control of their effects on trailing wings

The near-vortex-wake flow of a large aspect-ratio rectangular wing is accurately computed by using the thin-layer and full Navier-Stokes (NS) equations. The chordwise section of the wing is a NACA-0012 airfoil and its tip is round. The computations have been carried out on a fine C-O grid using an implicit, upwind, flux-difference splitting, finite-volume scheme. The thin-layer NS results have been obtained with and without flux limiters, and the full NS results have been obtained without flux limiters. Flow transition from laminar to turbulent is mimicked by turning-on the Baldwin-Lomax algebraic model at an experimentally prescribed chord-station location of 0.05. Comparison of computed results and experimental data shows that the full NS results give the best resolution of the near-vortex-wake flow. Next, the strength of the wing-tip vortex has been reduced substantially without reducing the lift coefficient by using flow-injection from a slot along a portion of the wing tip. The flow injection is directed in the wing plane at 45 deg with the wing-tip chord.

Wong, Tin-Chee↗
NASA NTRS1986

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗
NASA NTRS1994

Implementation and Assessment of Advanced Analog Vector-Matrix Processor

This paper discusses the design and implementation of an analog optical vecto-rmatrix coprocessor with a throughput of 128 Mops for a personal computer. Vector matrix calculations are inherently parallel, providing a promising domain for the use of optical calculators. However, to date, digital optical systems have proven too cumbersome to replace electronics, and analog processors have not demonstrated sufficient accuracy in large scale systems. The goal of the work described in this paper is to demonstrate a viable optical coprocessor for linear operations. The analog optical processor presented has been integrated with a personal computer to provide full functionality and is the first demonstration of an optical linear algebra processor with a throughput greater than 100 Mops. The optical vector matrix processor consists of a laser diode source, an acoustooptical modulator array to input the vector information, a liquid crystal spatial light modulator to input the matrix information, an avalanche photodiode array to read out the result vector of the vector matrix multiplication, as well as transport optics and the electronics necessary to drive the optical modulators and interface to the computer. The intent of this research is to provide a low cost, highly energy efficient coprocessor for linear operations. Measurements of the analog accuracy of the processor performing 128 Mops are presented along with an assessment of the implications for future systems. A range of noise sources, including cross-talk, source amplitude fluctuations, shot noise at the detector, and non-linearities of the optoelectronic components are measured and compared to determine the most significant source of error. The possibilities for reducing these sources of error are discussed. Also, the total error is compared with that expected from a statistical analysis of the individual components and their relation to the vector-matrix operation. The sufficiency of the measured accuracy of the processor is compared with that required for a range of typical problems. Calculations resolving alloy concentrations from spectral plume data of rocket engines are implemented on the optical processor, demonstrating its sufficiency for this problem. We also show how this technology can be easily extended to a 100 x 100 10 MHz (200 Cops) processor.

Gary, Charles K.↗
NASA NTRS1981

Interruption as a test of the user-computer interface

In order to study the effects different logic systems might have on interrupted operation, an algebraic calculator and a reverse polish notation calculator were compared when trained users were interrupted during problem entry. The RPN calculator showed markedly superior resistance to interruption effects compared to the AN calculator although no significant differences were found when the users were not interrupted. Causes and possible remedies for interruption effects are speculated. It is proposed that because interruption is such a common occurrence, it be incorporated into comparative evaluation tests of different logic system and control/display system and that interruption resistance be adopted as a specific design criteria for such design.

Kreifeldt, J. G.↗
NASA NTRS1984

Luneburg lens and optical matrix algebra research

Planar, as opposed to channelized, integrated optical circuits (IOCs) were stressed as the basis for computational devices. Both fully-parallel and systolic architectures are considered and the tradeoffs between the two device types are discussed. The Kalman filter approach is a most important computational method for many NASA problems. This approach to deriving a best-fit estimate for the state vector describing a large system leads to matrix sizes which are beyond the predicted capacities of planar IOCs. This problem is overcome by matrix partitioning, and several architectures for accomplishing this are described. The Luneburg lens work has involved development of lens design techniques, design of mask arrangements for producing lenses of desired shape, investigation of optical and chemical properties of arsenic trisulfide films, deposition of lenses both by thermal evaporation and by RF sputtering, optical testing of these lenses, modification of lens properties through ultraviolet irradiation, and comparison of measured lens properties with those expected from ray trace analyses.

Wood, V. E.↗
NASA NTRS1966

Matrix methods and automation in structural engineering

The partial differential equations of motion of shell structures subject to arbitrary loads may be extremely difficult or even impossible to solve. Since the digital computer is now an available tool for the structural engineer, current research has been directed toward methods which involve matrix formulation of large systems of algebraic equations and matrix solutions for discrete elements rather than the solution of the partial differential equations. These methods require rapid and accurate computer solutions. Recognizing the accuracy problems inherent in working with large matrices, a comprehensive survey of available computer programs was performed for large matrix inversion and for eigenvalue and eigenvector solutions under a Research Grant from the National Aeronautics and Space Administration and the techniques are developed in this paper for using these programs most efficiently for structural applications. The contents of the paper include: automation of matrix compilation, methods of very large matrix inversion and solution of simultaneous equation, and techniques for finding eigenvalues and vectors. In addition, a finite element stiffness matrix approach developed at the Denver Research Institute for both plates and shells subject to arbitrary dynamic loads is described as it was instrumented with complete automation on the digital computer.

Anita S West↗
NASA NTRS1984

Flow simulations for nacelle-propeller configurations using Euler equations

Euler codes for both axisymmetric and general three-dimensional nacelle-propeller flow analysis have been developed. Surface-fitted grids are generated either by a finite difference method or by an algebraic method. The propeller is represented by an actuator disk along a computational plane where proper boundary conditions are assigned to simulate the propeller power loading. Computed results for a NASA SR3 propeller as well as a NASA turboprop configuration are compared with test data. Good agreement has been achieved through the present simulation method.

Yu, N. J.↗
NASA NTRS1987

Viscous calculations of supersonic flows past cylinders at angles of attack

Two- and three-dimensional separated flows around cylindrical sections and a body-of-revolution at angles of attack from -10 deg to +20 deg are computed. The blunt-nose-cylinders chosen for the analyses have different configurations, and one of them has a boat-tail. The freestream Mach numbers are 1.5 and 1.6, the unit Reynolds number and the total temperature are 2 million/ft, and 585 R, respectively. Among the possible areas of application for this study are the aerodynamics of separating stores and missiles. The Reynolds-averaged Navier-Stokes equations are solved by an approximately-factored, implicit, upwind-biased, finite-volume algorithm. The convergence to steady state was accelerated using a multigrid strategy. The results helped understanding the structure of the flowfield with varying incidence and configuration, and the suitability of the computational algorithm used for this class of flows. A two-layer algebraic eddy viscosity model is used for the two-dimensional turbulent flows. Laminar three-dimensional calculations are benchmarked through comparisons with wind tunnel test results.

Baysal, O.↗
NASA NTRS1988

Calculations of transonic flows with shocks using Newton's method and direct solver. II - Solution of Euler equations

Transonic flows with shocks are simulated using steady Euler equations and by simultaneously solving the resulting nonlinear algebraic equations using Newton's method. At each iteration, a direct solver computes the corrections and the process is repeated until convergence is achieved. The corrections and errors are reduced quadratically with the present method, allowing solutions of machine accuracy to be obtained in a few steps. Nonunique inviscid solutions and nonunique solutions of the Navier Stokes equations for quasi-one-dimensional flows in nozzles are presented. Calculations are also presented for steady two-dimensional inviscid flows around a cylinder in the transonic regime.

Hafez, M.↗
NASA NTRS1991

Experiments with conjugate gradient algorithms for homotopy curve tracking

There are algorithms for finding zeros or fixed points of nonlinear systems of equations that are globally convergent for almost all starting points, i.e., with probability one. The essence of all such algorithms is the construction of an appropriate homotopy map and then tracking some smooth curve in the zero set of this homotopy map. HOMPACK is a mathematical software package implementing globally convergent homotopy algorithms with three different techniques for tracking a homotopy zero curve, and has separate routines for dense and sparse Jacobian matrices. The HOMPACK algorithms for sparse Jacobian matrices use a preconditioned conjugate gradient algorithm for the computation of the kernel of the homotopy Jacobian matrix, a required linear algebra step for homotopy curve tracking. Here, variants of the conjugate gradient algorithm are implemented in the context of homotopy curve tracking and compared with Craig's preconditioned conjugate gradient method used in HOMPACK. The test problems used include actual large scale, sparse structural mechanics problems.

Irani, Kashmira M.↗