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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Application of an upwind Navier-Stokes code to two-dimensional transonic airfoil flow

An upwind-biased implicit approximate factorization Navier-Stokes algorithm is applied to a variety of steady transonic airfoil cases, using the NACA 0012, RAE 2822, and Jones supercritical airfoils. The thin-layer form of the compressible Navier-Stokes equations is used. Both the CYBER 205 and CRAY 2 supercomputers are utilized, with average computational speeds of about 18 and 16 microsec/gridpoint/iteration, respectively. Lift curves, drag polars, and variations in drag coefficient with Mach number are determined for the NACA 0012 and Jones supercritical airfoils. Also, several cases are computed for comparison with experiment. The effect of grid density and grid extent on a typical turbulent airfoil solution is shown. An algebraic eddy-viscosity turbulence model is used for all of the computations.

Rumsey, Christopher L.↗
NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1965

Threshold logic - a synthesis approach.

Book on threshold logic for logical design of digital systems including threshold element synthesis by iteration, tabulation and network synthesis

NETWORK SYNTHESIS↗
NASA NTRS1972

Spectral factorization in periodically time-varying systems and application to navigation problems.

Spectral factorization has been used previously to derive the steady-state solution of Kalman filtering equations without iteration for constant coefficient systems. The present work extends the spectral factorization algorithm to time-varying systems having periodic coefficient matrices for cases of both discrete and continuous systems. Time-consuming, expensive iterations of sequential covariance equations are not required to reach the final solution since this is an algebraic algorithm employing existing eigenvalue, eigenvector subroutines. The computer program incorporating the algorithm is suitable for sensitivity studies in formulating navigation and guidance strategies of low-thrust interplanetary missions. The determination of an optimum tracking pattern from an earth station is examined as an example.

Nishimura, T.↗
NASA NTRS1977

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗
NASA NTRS1979

Programs for high-speed Fourier, Mellin and Fourier-Bessel transforms

Several FORTRAN program modules for performing one-dimensional and two-dimensional discrete Fourier transforms, Mellin, and Fourier-Bessel transforms are described along with programs that realize the algebra of high speed Fourier transforms on a computer. The programs can perform numerical harmonic analysis of functions, synthesize complex optical filters on a computer, and model holographic image processing methods.

Ikhabisimov, D. K.↗
NASA NTRS1982

Development and application of the GIM code for the Cyber 203 computer

The GIM computer code for fluid dynamics research was developed. Enhancement of the computer code, implicit algorithm development, turbulence model implementation, chemistry model development, interactive input module coding and wing/body flowfield computation are described. The GIM quasi-parabolic code development was completed, and the code used to compute a number of example cases. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and implicit finite difference scheme were also added. Development was completed on the interactive module for generating the input data for GIM. Solutions for inviscid hypersonic flow over a wing/body configuration are also presented.

Stainaker, J. F.↗
NASA NTRS1982

Comparison of three-dimensional turbulent boundary-layer calculations with experiment

A numerical simulation is carried out of the incompressible three-dimensional turbulent boundary-layer experiments of van den Berg and Elsenaar (1972), Dechow (1977), and Muller (1982). The simulation is effected by integrating the boundary-layer equations together with an algebraic eddy-viscosity turbulence model. With this computational method, the downstream portion of the flow treated by van den Berg and Elsenaar, where the crossflow is large, cannot be predicted. It is believed that this flow is significantly influenced by elliptic flow-field effects. Even though Dechow's experiment also indicated departures from the boundary-layer concept, the calculations made here agree reasonably well with the mean-flow development up to separation. The normal pressure gradients in Muller's experiment are found to be negligible in regions with large skewing, making it possible to test turbulence models using the boundary-layer equations. The simulation of this flow compares favorably with the experimental data throughout the flow field.

Muller, U. R.↗
NASA NTRS1985

Simultaneous Earth observations from 2 satellites

Simultaneous co-located observations from two different orbits lead to several advantages (i.e., cross calibration of sensors and a wider range of solar-zenith and sensor look angles). The question was asked how many times per year (on the average) do the sub-satellite points of two satellites simultaneously come within D kilometers of each other?. For the Space Station (altitude: 500 km, inclination: 28 deg) and a Sun synchronous satellite (altitude 705 km, inclination 98.21 deg) the answers are 16, 41 and 82 times per year for encounter distances D of 20, 50, and 100 km espectively. The relationship between encounters per year and distance D is linear. The answers were obtained in two ways: (1) a closed form statistical approach which led to a simple algebraic expression, and (2) a Monte Carlo type computer solution. The largest difference between the two solutions was less than 12%.

Montgomery, H. E.↗
NASA NTRS1986

Euler solutions for high-speed flow about complex three-dimensional configurations

A numerical algorithm based on a finite-volume explicit scheme with Runge-Kutta time integration of the Euler equations is presented for calculating high-speed three-dimensional flow about complex aerospace configurations. The use of enhancing factors such as artificial dissipative terms, enthalpy damping and local time-stepping are described. An algebraic method for generating quasi-three-dimensional computational grids for realistic aerospace configurations is presented. Computed results for various three-dimensional bodies at different Mach numbers and angles of attack have been obtained using the methods for grid-generation and flow simulation. Comparison of computed and experimental data for an advanced tactical aircraft-like configuration is presented, and a reasonable agreement of the data is noticed.

Moitra, A.↗
NASA NTRS1986

Simultaneous earth observations from two satellites

Simultaneous co-located observations from two different orbits lead to several advantages (i.e., cross calibration of sensors and a wider range of solar-zenith and sensor look angles). The question was asked how many times per year (on the average) do the sub-satellite points of two satellites simultaneously come within D kilometers of each other? For the Space Station (altitude: 500 km, inclination: 28 deg) and a Sun synchronous satellite (altitude 705 km, inclination 98.21 deg) the answers are 16, 41 and 82 times per year for encounter distances D of 20, 50, and 100 km respectively. The relationship between encounters per year and distance D is linear. The answers were obtained in two ways: (1) a closed form statistical approach which led to a simple algebraic expression, and (2) a Monte Carlo type computer solution. The largest difference between the two solutions was less than 12 percent.

Montgomery, H. E.↗
NASA NTRS1991

Linear and circular digital spectral analysis of serial data

Two methods of digital spectral analysis of unevenly sampled data are developed and illustrated here. One method uses a linear function of time (or space), the other uses circular functions. The circular method turns out to be essentially equivalent to a least-squares sine-wave analysis. The linear, anharmonic method uses only the field of real numbers and elementary algebraic operations, and hence it can be made computationally very fast and accurate. Both methods are very general, properly handling all kinds of time series ranging from simple series consisting only of the times of events to complicated series consisting of pulses with long duty cycles. The two methods are here applied to the analysis of annual mean relative sunspot numbers.

Stothers, Richard B.↗
NASA NTRS1992

Algebraic surface design and finite element meshes

Some of the techniques are summarized which are used in constructing C sup 0 and C sup 1 continuous meshes of low degree, implicitly defined, algebraic surface patches in three dimensional space. These meshes of low degree algebraic surface patches are used to construct accurate computer models of physical objects. These meshes are also used in the finite element simulation of physical phenomena (e.g., heat dissipation, stress/strain distributions, fluid flow characteristics) required in the computer prototyping of both the manufacturability and functionality of the geometric design.

Bajaj, Chandrajit L.↗
NASA NTRS1996

GENIE(++): A Multi-Block Structured Grid System

The computer code GENIE++ is a continuously evolving grid system containing a multitude of proven geometry/grid techniques. The generation process in GENIE++ is based on an earlier version. The process uses several techniques either separately or in combination to quickly and economically generate sculptured geometry descriptions and grids for arbitrary geometries. The computational mesh is formed by using an appropriate algebraic method. Grid clustering is accomplished with either exponential or hyperbolic tangent routines which allow the user to specify a desired point distribution. Grid smoothing can be accomplished by using an elliptic solver with proper forcing functions. B-spline and Non-Uniform Rational B-splines (NURBS) algorithms are used for surface definition and redistribution. The built in sculptured geometry definition with desired distribution of points, automatic Bezier curve/surface generation for interior boundaries/surfaces, and surface redistribution is based on NURBS. Weighted Lagrance/Hermite transfinite interpolation methods, interactive geometry/grid manipulation modules, and on-line graphical visualization of the generation process are salient features of this system which result in a significant time savings for a given geometry/grid application.

Williams, Tonya↗
NASA NTRS1980

Investigation of a Three-Dimensional Shock Wave Separated Turbulent Boundary Layer

A detailed Investigation of a flow in which a three-dimensional shock wave separates a two-dimensional turbulent boundary layer is presented. The resulting flowfield is highly three dimensional with a significant portion of flow separation on the surface at the phi = 0 deg (windward) plane was well as a large zone of secondary surface flow off this plane. Mean and fluctuating experimental measurements were obtained throughout the entire flowfield. These measurements included mean pressures, flow angles and shear on the surface, as well as yaw angles, static pressures, turbulent shear stresses, and turbulent kinetic energies on selected planes throughout the flowfield. In addition, numerical predictions of this flow, obtained by solving the Navier-Stokes equations with an algebraic eddy viscosity turbulence model, are presented. These computations reasonably predict both the surface and flowfield quantities, despite the extremely complicated nature of the experimental flow.

Kussoy, M. I.↗
NASA NTRS2008

Application of Fast Multipole Methods to the NASA Fast Scattering Code

The NASA Fast Scattering Code (FSC) is a versatile noise prediction program designed to conduct aeroacoustic noise reduction studies. The equivalent source method is used to solve an exterior Helmholtz boundary value problem with an impedance type boundary condition. The solution process in FSC v2.0 requires direct manipulation of a large, dense system of linear equations, limiting the applicability of the code to small scales and/or moderate excitation frequencies. Recent advances in the use of Fast Multipole Methods (FMM) for solving scattering problems, coupled with sparse linear algebra techniques, suggest that a substantial reduction in computer resource utilization over conventional solution approaches can be obtained. Implementation of the single level FMM (SLFMM) and a variant of the Conjugate Gradient Method (CGM) into the FSC is discussed in this paper. The culmination of this effort, FSC v3.0, was used to generate solutions for three configurations of interest. Benchmarking against previously obtained simulations indicate that a twenty-fold reduction in computational memory and up to a four-fold reduction in computer time have been achieved on a single processor.

Dunn, Mark H.↗
NASA NTRS2008

Simulink Model of the Ares I Upper Stage Main Propulsion System

A numerical model of the Ares I upper stage main propulsion system is formulated based on first principles. Equation's are written as non-linear ordinary differential equations. The GASP fortran code is used to compute thermophysical properties of the working fluids. Complicated algebraic constraints are numerically solved. The model is implemented in Simulink and provides a rudimentary simulation of the time history of important pressures and temperatures during re-pressurization, boost and upper stage firing. The model is validated against an existing reliable code, and typical results are shown.

Burchett, Bradley T.↗