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 685 records · Page 38

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

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

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

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

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↗

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

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

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

Numerical and Analytical Assessment of a Coupled Rotating Detonation Engine and Turbine Experiment

An analysis is presented of an experimental rig comprising a rotating detonation engine (RDE) with bypass ejector flow coupled to a downstream turbine. The analysis used a validated computational fluid dynamics RDE simulation combined with straightforward algebraic mixing equations for the bypass flow. The objectives of the analysis were to supplement and interpret the necessarily sparse measurements from the rig, and to assess the performance of the RDE itself, which was not instrumented in this installation. The analysis is seen to agree reasonably well with available data. It shows that the RDE is operating in an unusual fashion, with subsonic flow throughout the exhaust plane. The detonation event itself is producing a total pressure rise relative to the pre-detonative pressure; however, the length of the device and the substantial flow restriction at the inlet yield an overall pressure loss. This is not surprising since the objective of the rig test was primarily aimed at investigating RDE/turbine interactions, and not on performance optimization. Furthermore, the RDE was designed for fundamental detonation studies and not performance. Nevertheless, the analysis indicates that with some small alterations to the design, an RDE with an overall pressure rise is possible.

detonation↗

Numerical and Analytical Assessment of a Coupled Rotating Detonation Engine and Turbine Experiment

An analysis is presented of an experimental rig comprising a rotating detonation engine (RDE) with bypass ejector flow coupled to a downstream turbine. The analysis used a validated computational fluid dynamics RDE simulation combined with straightforward algebraic mixing equations for the bypass flow. The objectives of the analysis were to supplement and interpret the necessarily sparse measurements from the rig, and to assess the performance of the RDE itself, which was not instrumented in this installation. The analysis is seen to agree reasonably well with available data. It shows that the RDE is operating in an unusual fashion, with subsonic flow throughout the exhaust plane. The detonation event itself is producing a total pressure rise relative to the pre-detonative pressure; however, the length of the device and the substantial flow restriction at the inlet yield an overall pressure loss. This is not surprising since the objective of the rig test was primarily aimed at investigating RDE turbine interactions, and not on performance optimization. Furthermore, the RDE was designed for fundamental detonation studies and not performance. Nevertheless, the analysis indicates that with some small alterations to the design, an RDE with an overall pressure rise is possible.

propulsion↗

Numerical and Analytical Assessment of a Coupled Rotating Detonation Engine and Turbine Experiment

An analysis is presented of an experimental rig comprising a rotating detonation engine (RDE) with bypass flow coupled to a downstream turbine. The analysis used a validated computational fluid dynamics RDE simulation combined with straightforward algebraic mixing equations for the bypass flow. The objectives of the analysis were to supplement and interpret the necessarily sparse measurements from the rig, and to assess the performance of the RDE itself (which was not instrumented in this installation). The analysis is seen to agree reasonably well with available data. It shows that the RDE is operating in an unusual fashion, with subsonic flow throughout the exhaust plane. The detonation event itself is producing a total pressure rise relative to the pre-detonative pressure; however, the length of the device and the substantial flow restriction at the inlet yield an overall pressure loss. This is not surprising since the objective of the rig test was primarily aimed at investigating RDEturbine interactions, and not on performance optimization. Furthermore, the RDE was designed for fundamental detonation studies and not performance. Nevertheless, the analysis indicates that with some small alterations to the design, an RDE with an overall pressure rise is possible.

detonation↗

Heat Transfer on a Flat Plate with Uniform and Step Temperature Distributions

Heat transfer associated with turbulent flow on a step-heated or cooled section of a flat plate at zero angle of attack with an insulated starting section was computationally modeled using the GASP Navier-Stokes code. The algebraic eddy viscosity model of Baldwin-Lomax and the turbulent two-equation models, the K- model and the Shear Stress Turbulent model (SST), were employed. The variations from uniformity of the imposed experimental temperature profile were incorporated in the computations. The computations yielded satisfactory agreement with the experimental results for all three models. The Baldwin- Lomax model showed the closest agreement in heat transfer, whereas the SST model was higher and the K-omega model was yet higher than the experiments. In addition to the step temperature distribution case, computations were also carried out for a uniformly heated or cooled plate. The SST model showed the closest agreement with the Von Karman analogy, whereas the K-omega model was higher and the Baldwin-Lomax was lower.

Bahrami, Parviz A.↗

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗

Algebraic discrete quantum harmonic oscillator with dynamic resolution scaling

We develop an algebraic formulation for the discrete quantum harmonic oscillator (DQHO) from the Hamiltonian for two, coupled QHOs and provide a physical picture for the Kravchuk function eigenstates of the oscillator. The familiar $\mathfrak{su}(2)$ structure of the coupled QHO Hamiltonian divides its spectrum into sets corresponding to the DQHO at different resolutions. In addition to energy ladder operators, the formulation allows for the introduction of resolution ladder operators connecting all DQHOs with different resolutions, thus enabling the dynamic scaling of the resolution of finite degree-of-freedom quantum simulations. The coherent state of the DQHO is constructed, and its expected position is proven to oscillate as a classical harmonic oscillator. The DQHO coherent state recovers that of the quantum harmonic oscillator at large resolution.

97 MATHEMATICS AND COMPUTING↗

State-Space System Realization with Input- and Output-Data Correlation

This paper introduces a general version of the information matrix consisting of the autocorrelation and cross-correlation matrices of the shifted input and output data. Based on the concept of data correlation, a new system realization algorithm is developed to create a model directly from input and output data. The algorithm starts by computing a special type of correlation matrix derived from the information matrix. The special correlation matrix provides information on the system-observability matrix and the state-vector correlation. A system model is then developed from the observability matrix in conjunction with other algebraic manipulations. This approach leads to several different algorithms for computing system matrices for use in representing the system model. The relationship of the new algorithms with other realization algorithms in the time and frequency domains is established with matrix factorization of the information matrix. Several examples are given to illustrate the validity and usefulness of these new algorithms.

Juang, Jer-Nan↗

CFD applications in pump flows

The objective of the proposed paper is to develop a computational procedure that solves incompressible Navier-Stokes equations for pump flows. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. As a benchmark problem, the flow through the Rocketdyne inducer is numerically simulated. A coarse grid solution is obtained with a single zone by using an algebraic turbulence model. In multi-zone fine grid computation, a one-equation Baldwin-Barth turbulence model is utilized. Numerical results are compared with experimental measurements and a good agreement is found between the two. The resulting computer code is then applied to the flow analysis inside a two-stage fuel pump impeller operating at 80 percent, 100 percent, and 120 percent of design flow.

Kiris, Cetin↗

An Algebra of Machine Learners with Applications

Machine learning (ML) methods are increasingly being applied to solve complex, data-driven problems in diverse areas, by exploiting the physical laws derived from first principles such as thermal hydraulics and the abstract laws developed recently for data and computing infrastructures. These physical and abstract laws encapsulate, typically in compact algebraic forms, the critical knowledge that complements data-driven ML models. We present a unified perspective of these laws and ML methods using an abstract algebra (A;⊕,Ⓧ), wherein the performance estimation and classification tasks are characterized by the additive ⊕ operations, and the diagnosis, reconstruction, and optimization tasks are characterized by the difference Ⓧ operations. This abstraction provides ML codes and their performance characterizations that are transferable across different areas. We describe practical applications of these abstract operations using examples of throughput profile estimation tasks in data transport infrastructures, and power-level and sensor error estimation tasks in nuclear reactor systems.

Rao, Nageswara↗