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 577 records · Page 32

Theories of resonant satellite pairs in Saturn's system

Analytical theories have been devised for the three resonant pairs in Saturn's satellite system: Enceladus-Dione, Titan-Hyperion, and Mimas-Tethys. An algebraic manipulation language was adopted to calculate the theories; use of the canonically invariant Hori-Lie (1966) theory in noncanonical variables made it possible to perform direct computations of perturbations in such quantities as longitude, latitude and radius vector. Important terms at the 1-100 km level have been found second-order terms of the resonant satellite are calculated.

Jefferys, W. H.↗

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

Nonlinear analysis of compressively loaded linked-plate structures

A method is presented for the geometrically nonlinear analysis of the equilibrium behavior of a general class of compressively loaded, prismatic, linked-plate structures. The semi-analytical method accounts for geometric shape imperfections, and allows the use of laminated composite plate properties. Buckling eigenfunctions are used as displacement shape functions, and the associated second-order displacement fields are incorporated. The nonlinear plate equilibrium equations are expressed in terms of the assumed displacement form, and the Galerkin method is used to obtain a set of algebraic equations which are nonlinear in the modal amplitudes and the load parameter. Buckling eigensolutions for the class of structures under consideration can be obtained from the VIPASA computer code; thus, subject to some additional restrictions in loading and plate anisotropy, structural configurations which can be analyzed with regard to buckling response by VIPASA can be analyzed with regard to geometrically nonlinear response using the new method.

Stoll, Frederick↗

Two Axis Pointing System (TAPS) attitude acquisition, determination, and control

The Two Axis Pointing System (TAPS) is a 2 axis gimbal system designed to provide fine pointing of Space Transportation System (STS) borne instruments. It features center-of-mass instrument mounting and will accommodate instruments of up to 1134 kg (2500 pounds) which fit within a 1.0 by 1.0 by 4.2 meter (40 by 40 by 166 inch) envelope. The TAPS system is controlled by a microcomputer based Control Electronics Assembly (CEA), a Power Distribution Unit (PDU), and a Servo Control Unit (SCU). A DRIRU-II inertial reference unit is used to provide incremental angles for attitude propagation. A Ball Brothers STRAP star tracker is used for attitude acquisition and update. The theory of the TAPS attitude determination and error computation for the Broad Band X-ray Telescope (BBXRT) are described. The attitude acquisition is based upon a 2 star geometric solution. The acquisition theory and quaternion algebra are presented. The attitude control combines classical position, integral and derivative (PID) control with techniques to compensate for coulomb friction (bias torque) and the cable harness crossing the gimbals (spring torque). Also presented is a technique for an adaptive bias torque compensation which adjusts to an ever changing frictional torque environment. The control stability margins are detailed, with the predicted pointing performance, based upon simulation studies. The TAPS user interface, which provides high level operations commands to facilitate science observations, is outlined.

Azzolini, John D.↗

Techniques for grid manipulation and adaptation

Two approaches have been taken to provide systematic grid manipulation for improved grid quality. One is the control point form (CPF) of algebraic grid generation. It provides explicit control of the physical grid shape and grid spacing through the movement of the control points. It works well in the interactive computer graphics environment and hence can be a good candidate for integration with other emerging technologies. The other approach is grid adaptation using a numerical mapping between the physical space and a parametric space. Grid adaptation is achieved by modifying the mapping functions through the effects of grid control sources. The adaptation process can be repeated in a cyclic manner if satisfactory results are not achieved after a single application.

Choo, Yung K.↗

A recurrence matrix solution for the dynamic response of aircraft in gusts

A systematic procedure developed for the calculation of the structural response of aircraft flying through a gust by use of difference equations in the solution of dynamic problems is first illustrated by means of a simple-damped-oscillator example. A detailed analysis is then given which leads to a recurrence matrix equation for the determination of the response of an airplane in a gust. The method takes into account wing bending and twisting deformations, fuselage deflection, vertical and pitching motion of the airplane, and some tail forces. The method is based on aerodynamic strip theory, but compressibility and three-dimensional aerodynamic effects can be taken into account approximately by means of over-all corrections. Either a sharp-edge gust or a gust of arbitrary shape in the spanwise or flight directions may be treated. In order to aid in the application of the method to any specific case, a suggested computational procedure is included. The possibilities of applying the method to a variety of transient aircraft problems, such as landing, are brought out. A brief review of matrix algebra, covering the extent to which it is used in the analysis, is also included. (author)

Houbolt, John C↗

ADPAC v1.0: User's Manual

The overall objective of this study was to evaluate the effects of turbulence models in a 3-D numerical analysis on the wake prediction capability. The current version of the computer code resulting from this study is referred to as ADPAC v7 (Advanced Ducted Propfan Analysis Codes -Version 7). This report is intended to serve as a computer program user's manual for the ADPAC code used and modified under Task 15 of NASA Contract NAS3-27394. The ADPAC program is based on a flexible multiple-block and discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Turbulence models now available in the ADPAC code are: a simple mixing-length model, the algebraic Baldwin-Lomax model with user defined coefficients, the one-equation Spalart-Allmaras model, and a two-equation k-R model. The consolidated ADPAC code is capable of executing in either a serial or parallel computing mode from a single source code.

Hall, Edward J.↗

Numerical study of scramjet and ramjet flow fields

Two computer programs have been developed to numerically calculate complex, two-dimensional flow fields in scramjets. The first program is written for inlet analysis whereas the second program is written primarily for combustor analysis. Both programs solve the full two-dimensional Navier-Stokes equations by a well-known explicit, predictor-corrector technique. Turbulence is modeled by an algebraic eddy-viscosity model. The combustor program also includes one or more species conservation equations to calculate mixing and reacting flows. The hydrogen/air chemistry in this program is modeled by a complete reaction model. The combustor program has been recently modified to analyze axisymmetric ramjet dump combustor flow field. Results from these computer programs are presented that predict the flow in several scramjet inlet configurations, two model scramjet engine configurations, and in a dump combustor simulator. Computed results are also compared with available experimental data to allow assessment of the programs.

Kumar, A.↗

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

Data Retrieval Algorithms for Validating the Optical Transient Detector and the Lightning Imaging Sensor

A linear algebraic solution is provided for the problem of retrieving the location and time of occurrence of lightning ground strikes from an Advanced Lightning Direction Finder (ALDF) network. The ALDF network measures field strength, magnetic bearing, and arrival time of lightning radio emissions. Solutions for the plane (i.e., no earth curvature) are provided that implement all of these measurements. The accuracy of the retrieval method is tested using computer-simulated datasets, and the relative influence of bearing and arrival time data an the outcome of the final solution is formally demonstrated. The algorithm is sufficiently accurate to validate NASA:s Optical Transient Detector and Lightning Imaging Sensor. A quadratic planar solution that is useful when only three arrival time measurements are available is also introduced. The algebra of the quadratic root results are examined in detail to clarify what portions of the analysis region lead to fundamental ambiguities in sc)iirce location, Complex root results are shown to be associated with the presence of measurement errors when the lightning source lies near an outer sensor baseline of the ALDF network. For arbitrary noncollinear network geometries and in the absence of measurement errors, it is shown that the two quadratic roots are equivalent (no source location ambiguity) on the outer sensor baselines. The accuracy of the quadratic planar method is tested with computer-generated datasets, and the results are generally better than those obtained from the three-station linear planar method when bearing errors are about 2 deg.

Koshak, W. J.↗

On the representation matrices of the spin permutation group

A method is presented for the determination of the representation matrices of the spin permutation group (symmetric group), a detailed knowledge of these matrices being required in the study of the electronic structure of atoms and molecules. The method is characterized by the use of two different coupling schemes. Unlike the Yamanouchi spin algebraic scheme, the method is not recursive. The matrices for the fundamental transpositions can be written down directly in one of the two bases. The method results in a computationally significant reduction in the number of matrix elements that have to be stored when compared with, say, the standard Young tableaux group theoretical approach.

Wilson, S.↗

Application of the implicit MacCormack scheme to the PNS equations

The two-dimensional parabolized Navier-Stokes equations are solved using MacCormack's (1981) implicit finite-difference scheme. It is shown that this method for solving the parabolized Navier-Stokes equations does not require the inversion of block tridiagonal systems of algebraic equations and allows the original explicit scheme to be employed in those regions where implicit treatment is not needed. The finite-difference algorithm is discussed and the computational results for two laminar test cases are presented. Results obtained using this method for the case of a flat plate boundary layer are compared with those obtained using the conventional Beam-Warming scheme, as well as those obtained from a boundary layer code. The computed results for a more severe test of the method, the hypersonic flow past a 15 deg compression corner, are found to compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.

Lawrence, S. L.↗

Application of Artificial Intelligence technology to the analysis and synthesis of reliable software systems

The development of a methodology for the production of highly reliable software is one of the greatest challenges facing the computer industry. Meeting this challenge will undoubtably involve the integration of many technologies. This paper describes the use of Artificial Intelligence technologies in the automated analysis of the formal algebraic specifications of abstract data types. These technologies include symbolic execution of specifications using techniques of automated deduction and machine learning through the use of examples. On-going research into the role of knowledge representation and problem solving in the process of developing software is also discussed.

Wild, Christian↗

Investigation of third-order closure model of turbulence for the computation of incompressible flows in a channel with a backward-facing step

To predict the diffusion process of the Reynolds stresses in reattaching shear flows, the transport model for the triple-velocity products has been developed and tested for the computation of the flow in a channel with a backward-facing step. Upon comparison of the results of uuv, uvv, and vvv with those obtained by using existing algebraic correlations, it was shown that the present model improved the prediction of the triple-velocity products.

Amano, R. S.↗

Navier-Stokes calculations of scramjet-afterbody flowfields

Scramjet simulant-gas exhaust flows are computed using two-dimensional Navier-Stokes (NS) equations. The NS equations are solved using an implicit, upwind, finite-volume scheme, and the Reynolds stresses are modeled algebraically. The solutions are compared with experimental data where the freestream and the nozzle exhaust fluids are air. The NS equations for the flow, where the scramjet exhaust simulant gas is a mixture of Freon-12 and argon, are solved by an explicit, finite difference scheme. The results are compared with experimental data, and certain discrepancies are explained.

Baysal, O.↗

The symmetric turbulent plane wake downstream of a sharp trailing edge

The analysis and modeling of the symmetric turbulent plane wake downstream of a sharp trailing edge is addressed. A compact description of the flow near the trailing edge is formulated using the results of a previous asymptotic analysis. The new description retains the two-layered structure identified in the previous work, and it clarifies the principal dynamics of the flow in the near-wake outer layer, away from the wake centerline. For zero-pressure-gradient flow, the near-wake outer layer is shown to be represented to leading order by the similarity solution that governs the outer region of the surface boundary layer. The leading perturbation in the outer layer due to the developing near-wake inner-layer flow is identified, and this is shown to be asymptotically smaller than undetermined higher-order terms associated with the surface boundary-layer flow. Results of the new near-wake analysis are used to formulate an algebraic eddy viscosity model for wake flow predictions at arbitrary distances from the trailing edge. The model is used in a numerical solution of the boundary layer equations, and computed velocity and Reynolds stress profiles are shown to compare well with experimental data.

Bogucz, E. A.↗

Viscous Analysis of Ref. H Based Wing/Bodies

Computations have been performed on the baseline Reference H wing/body configuration, as well as the Wing 704 configuration, an optimized wing and fuselage combination derived from Ref. H through automated optimization. The parabolized Navier-Stokes solver UPS was employed with viscous terms in two directions in an effort to understand the source and level of potential viscous/inviscid interactions. The paper briefly describes the UPS code and the grids used to obtain the solutions before the discussion of results. Results of these computations indicate that viscous/inviscid interaction can contribute increments to both the pressure- and friction-related drag. Computations were performed for wind tunnel conditions-1.675% scale models at a Reynolds number of 4 million per foot. Turbulent flow results were obtained using the Baldwin-Lomax algebraic turbulence model and were compared with laminar flow results. The laminar flow fields were used to obtain upper bounds on potential interaction effects.

Lawrence, Scott L.↗