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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Entanglement entropy across the lattice-continuum correspondence

This paper revisits the standard calculations of free field entanglement entropy in light of the newly developed lattice-continuum correspondence. This correspondence prescribes an explicit method to extract an approximately continuum quantum field theory out of a fully regularized lattice theory. This prescription will here be extended to subregion algebras, and it will be shown how entropies of continuum boson and fermion theories can be computed by working purely with lattice quantities. This gives a clear picture of the origin of divergences in entanglement entropy while also presenting a concise and detailed recipe for calculating this important quantity in continuum theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

A linear programming manual

Computer solutions of linear programming problems are outlined. Information covers vector spaces, convex sets, and matrix algebra elements for solving simultaneous linear equations. Dual problems, reduced cost analysis, ranges, and error analysis are illustrated.

Tuey, R. C.↗
NASA NTRS1974

Discretized partial differential equations - Examples of control systems defined on modules

The purpose of this paper is to show how the important problems of linear system theory can be solved concisely for a particular class of linear systems, namely block circulant systems, by exploiting the algebraic structure. This type of system arises in lumped approximations to linear partial differential equations. The computation of the transition matrix, the variation of constants formula, observability, controllability, pole allocation, realization theory, stability and quadratic optimal control are discussed. In principle, all questions which are solved here could also be solved by standard methods; the present paper clearly exposes the structure of the solution, and thus permits various savings in computational effort.

Brockett, R. W.↗
NASA NTRS1974

Nonlinear finite-element analysis of laminated composite shells

Two aspects of the geometrically nonlinear analysis of laminated composite shells are considered in this paper. The first pertains to the simplifications and computational advantages resulting from the use of group-theoretic methods in conjunction with computerized symbolic (algebraic) manipulation in the element development. The second aspect pertains to the combined effects of shear deformation and nonlinearity on the accuracy and convergence of two-dimensional, shear-flexible stiffness finite element models.

Noor, A. K.↗
NASA NTRS1979

Orbital theories of outer planet satellites

An analytical theory of Mimas and Tethys was developed, retaining terms which produce perturbations of the order of + or - 10 km. The theory uses a novel set of variables, and was developed with the Hori-Lie algorithm, using the algebraic manipulation language TRIGMAN. the perturbations were implemented by means of FORTRAN subroutines produced by the computer.

Jefferys, W. H.↗
NASA NTRS1984

Pitch angle diffusion in the Jovian magnetodisc

Using a linear approximation to the field variation and a mathematical technique for distorting the orbit integration into the complex gyrophase plane, a formula is derived for the nonadiabatic change in magnetic moment. This change varies as the cosine of the particle gyrophase at the equator, so that if a particle traverses the equator several times with different values of this gyrophase, it suffers randomly correlated values for the nonadiabatic change in magnetic moment. The algebraic formula for this nonadiabatic change is found to be in excellent agreement with previous numerical computations. The phase space density, averaged over several bounce periods, is found to satisfy a diffusion equation in pitch angle. The diffusion coefficient is dependent on particle energy, pitch angle, and the field line along which the particles are moving.

Birmingham, T. J.↗
NASA NTRS1986

Improvement of the second- and third-moment modeling of turbulence: A study of Reynolds-stress closure model

Four parts of the Reynolds-stress closure modeling are reported: (1) improvement of the k and epsilon equaitons; (2) development of the third-moment transport equation; (3) formulation of the diffusion coefficient of the momentum equation by using the algebraic-stress model of turbulence; and (4) the application of the Reynolds-stress model to a heat exchanger problem. It was demonstrated that the third-moment transport model improved the prediction of the triple-velocity products in the recirculating and reattaching flow regions in comparison with the existing algebraic models for the triple-velocity products. Optimum values for empirical coefficients are obtained for the prediction of the backward-facing step flows. A functional expression is derived for the coefficient of the momentum diffusion by employing the algebraic-stress model. The second-moment closure is applied to a heat transfer problem. The computations for the flow in a corrugated-wall channel show that the second-moment closure improves the prediction of the heat transfer rates by 30% over the k - epsilon model.

Amano, R. S.↗
NASA NTRS1989

Image understanding and the man-machine interface II; Proceedings of the Meeting, Los Angeles, CA, Jan. 17, 18, 1989

Image understanding concepts and models, image understanding systems and applications, advanced digital processors and software tools, and advanced man-machine interfaces are among the topics discussed. Particular papers are presented on such topics as neural networks for computer vision, object-based segmentation and color recognition in multispectral images, the application of image algebra to image measurement and feature extraction, and the integration of modeling and graphics to create an infrared signal processing test bed.

Barrett, Eamon B.↗
NASA NTRS1992

A comparative study of turbulence models for overset grids

The implementation of two different types of turbulence models for a flow solver using the Chimera overset grid method is examined. Various turbulence model characteristics, such as length scale determination and transition modeling, are found to have a significant impact on the computed pressure distribution for a multielement airfoil case. No inherent problem is found with using either algebraic or one-equation turbulence models with an overset grid scheme, but simulation of turbulence for multiple-body or complex geometry flows is very difficult regardless of the gridding method. For complex geometry flowfields, modification of the Baldwin-Lomax turbulence model is necessary to select the appropriate length scale in wall-bounded regions. The overset grid approach presents no obstacle to use of a one- or two-equation turbulence model. Both Baldwin-Lomax and Baldwin-Barth models have problems providing accurate eddy viscosity levels for complex multiple-body flowfields such as those involving the Space Shuttle.

Renze, Kevin J.↗
NASA NTRS1994

Program Generates Two-Dimensional Computational Grids

TDIGG is fast and versatile computer program for generating two-dimensional computational grids for use in programs solving equations of flow by finite-difference methods. Both algebraic and elliptic grid-generation systems included. Enables user to view results of each iteration. Written in FORTRAN 77.

Vu, Bruce T.↗
NASA NTRS1998

User Manual for Beta Version of TURBO-GRD: A Software System for Interactive Two-Dimensional Boundary/ Field Grid Generation, Modification, and Refinement

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.

Choo, Yung K.↗
NASA NTRS1991

Viscous three-dimensional calculations of transonic fan performance

A 3-D flow analysis code was used to compute the design speed operating line of a transonic fan rotor, and the results were compared with experimental data. The code is an explicit finite difference code with an algebraic turbulence model. The transonic fan, called rotor 67, was tested experimentally at NASA-Lewis with conventional aerodynamic probes and with user anemometry and was included as one of the AGARD test cases for the computation of internal flows. The experimental data are described. Maps of total pressure ratio and adiabatic efficiency versus mass flow were computed and are compared with the experimental maps, with good agreement. Detailed comparisons between calculations and experiment are made at two operating points, one near peak efficiency and the other near stall. Blade-to-blade contour plots are used to show the shock structure. Comparisons of circumferentially integrated flow quantities downstream of the rotor show spanwise distributions of several aerodynamic parameters. Calculated Mach number distributions are compared with laser anemometer data within the blade row and the wake to quantify the accuracy of the calculations. Particle traces are used to show the nature of secondary flow.

Chima, Rodrick V.↗
NASA NTRS1978

Numerical solution of two-dimensional turbulent blunt body flows with an impinging shock

An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.

Tannehill, J. C.↗
NASA NTRS1991

An exact algebraic solution of the infimum in H-infinity optimization with output feedback

This paper presents a simple and noniterative procedure for the computation of the exact value of the infimum in the standard H-infinity-optimal control with output feedback. The problem formulation is general and does not place any restrictions on the direct feedthrough terms between the control input and the controlled output variables, and between the disturbance input and the measurement output variables. The method is applicable to systems that satisfy (1) the transfer function from the control input to the controlled output is right-invertible and has no invariant zeros on the j(w) axis and, (2) the transfer function from the disturbance to the measurement output is left-invertible and has no invariant zeros on the j(w) axis. A set of necessary and sufficient conditions for the solvability of H-infinity-almost disturbance decoupling problem via measurement feedback with internal stability is also given.

Chen, Ben M.↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Navier-Stokes computations of symmetric and asymmetric vortex shedding around slender bodies

The unsteady, compressible Reynolds-averaged thin-layer Navier-Stokes equations are used to solve for symmetric and asymmetric vortical flows around slender pointed bodies of revolution. The modified Baldwin and Lomax algebraic two-layer turbulent model is used for the eddy viscosity calculation in prescribed turbulent regions. The implicit, upwind flux-difference splitting finite-volume scheme is used to obtain the solutions. Computational results are presented for a low-speed symmetrical vortical flow around a 3.5-caliber tangent-ogive cylinder for fully laminar flow. Computed results are presented for a high-speed asymmetric vortex-shedding flow around a cone at 3.286 relative incidence, using a fine computational grid. The asymmetric vortex-shedding flow has been produced through a small asymmetric transitional perturbation in the cross-flow plane.

Kandil, Osama A.↗
NASA NTRS1993

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.

Hodges, Dewey H.↗