Search NASA⌕ Search

SEARCH · Search NASA

Results for “System of Linear Equations”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Newton-like minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗
NASA NTRS1989

The FORCE: A portable parallel programming language supporting computational structural mechanics

This project supports the conversion of codes in Computational Structural Mechanics (CSM) to a parallel form which will efficiently exploit the computational power available from multiprocessors. The work is a part of a comprehensive, FORTRAN-based system to form a basis for a parallel version of the NICE/SPAR combination which will form the CSM Testbed. The software is macro-based and rests on the force methodology developed by the principal investigator in connection with an early scientific multiprocessor. Machine independence is an important characteristic of the system so that retargeting it to the Flex/32, or any other multiprocessor on which NICE/SPAR might be imnplemented, is well supported. The principal investigator has experience in producing parallel software for both full and sparse systems of linear equations using the force macros. Other researchers have used the Force in finite element programs. It has been possible to rapidly develop software which performs at maximum efficiency on a multiprocessor. The inherent machine independence of the system also means that the parallelization will not be limited to a specific multiprocessor.

Jordan, Harry F.↗
NASA NTRS1995

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗
NASA NTRS1973

An improved method for the aerodynamic analysis of wing-body-tail configurations in subsonic and supersonic flow. Part 1: Theory and application

A new method has been developed for calculating the pressure distribution and aerodynamic characteristics of wing-body-tail combinations in subsonic and supersonic potential flow. A computer program has been developed to perform the numerical calculations. The configuration surface is subdivided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. The normal components of velocity induced at specified control points by each singularity distribution are calculated and make up the coefficients of a system of linear equations relating the strengths of the singularities to the magnitude of the normal velocities. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations using an interactive procedure. Once the singularity strengths are known, the pressure coefficients are calculated, and the forces and moments acting on the configuration determined by numerical integration.

Woodward, F. A.↗
NASA NTRS1990

Doppler radar imaging of spherical planetary surfaces

This paper describes a technique for using echo power-spectra for making radar images of spherical planetary targets. In developing the Doppler-radar imaging system, the target's reflectivity distribution is expanded in a truncated spherical harmonic series; the distribution of echo power in rotational phase and the Doppler frequency are obtained in form of a system of linear equations which express the echo spectra as functions of the series coefficients. To estimate the coefficients from an observed phase-Doppler distribution, the inversion is cast as a least-squares problem, and is solved using singular value decompositions, yielding a linear imaging system which can be fully characterized by its impulse response as a function of latitude. Simulations were designed and used to explore the capabilities of Doppler radar imaging, and the sensitivity of imaging to such factors as subradar altitude coverage and signal-to-noise ratio was investigated.

Hudson, R. Scott↗
NASA NTRS1992

Minimization of deviations of gear real tooth surfaces determined by coordinate measurements

The deviations of a gear's real tooth surface from the theoretical surface are determined by coordinate measurements at the grid of the surface. A method was developed to transform the deviations from Cartesian coordinates to those along the normal at the measurement locations. Equations are derived that relate the first order deviations with the adjustment to the manufacturing machine-tool settings. The deviations of the entire surface are minimized. The minimization is achieved by application of the least-square method for an overdetermined system of linear equations. The proposed method is illustrated with a numerical example for hypoid gear and pinion.

Litvin, F. L.↗
NASA NTRS1984

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner.

Bayliss, A.↗
NASA NTRS1990

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗
NASA NTRS1992

Mapping unstructured grid computations to massively parallel computers

Investigated here is this mapping problem: assign the tasks of a parallel program to the processors of a parallel computer such that the execution time is minimized. First, a taxonomy of objective functions and heuristics used to solve the mapping problem is presented. Next, we develop a highly parallel heuristic mapping algorithm, called Cyclic Pairwise Exchange (CPE), and discuss its place in the taxonomy. CPE uses local pairwise exchanges of processor assignments to iteratively improve an initial mapping. A variety of initial mapping schemes are tested and recursive spectral bipartitioning (RSB) followed by CPE is shown to result in the best mappings. For the test cases studied here, problems arising in computational fluid dynamics and structural mechanics on unstructured triangular and tetrahedral meshes, RSB and CPE outperform methods based on simulated annealing. Much less time is required to do the mapping and the results obtained are better. Compared with random and naive mappings, RSB and CPE reduce the communication time two fold for the test problems used. Finally, we use CPE in two applications on a CM-2. The first application is a data parallel mesh-vertex upwind finite volume scheme for solving the Euler equations on 2-D triangular unstructured meshes. CPE is used to map grid points to processors. The performance of this code is compared with a similar code on a Cray-YMP and an Intel iPSC/860. The second application is parallel sparse matrix-vector multiplication used in the iterative solution of large sparse linear systems of equations. We map rows of the matrix to processors and use an inner-product based matrix-vector multiplication. We demonstrate that this method is an order of magnitude faster than methods based on scan operations for our test cases.

Hammond, Steven Warren↗
NASA NTRS1996

Asymptotic solution of the turbulent mixing layer for velocity ratio close to unity

The equations describing the first two terms of an asymptotic expansion of the solution of the planar turbulent mixing layer for values of the velocity ratio close to one are obtained. The first term of this expansion is the solution of the well-known time-evolving problem and the second, which includes the effects of the increase of the turbulence scales in the stream-wise direction, obeys a linear system of equations. Numerical solutions of these equations for a two-dimensional reacting mixing layer show that the correction to the time-evolving solution may explain the asymmetry of the entrainment and the differences in product generation observed in flip experiments.

Higuera, F. J.↗
NASA NTRS1992

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗
NASA NTRS1971

The problem of an elastic stiffener bonded to a half plane.

The contact problem of an elastic stiffener bonded to an elastic half plane with different mechanical properties is considered. The governing integral equation is reduced to an infinite system of linear algebraic equations. It is shown that, depending on the value of a parameter which is a function of the elastic constants and the thickness of the stiffener, the system is either regular or quasi-regular. A complete numerical example is given for which the strength of the stress singularity and the contact stresses are tabulated.

Erdogan, F.↗
NASA NTRS1982

Evaluation of automated decisionmaking methodologies and development of an integrated robotic system simulation. Volume 2, Part 2: Appendixes B, C, D and E

The derivation of the equations is presented, the rate control algorithm described, and simulation methodologies summarized. A set of dynamics equations that can be used recursively to calculate forces and torques acting at the joints of an n link manipulator given the manipulator joint rates are derived. The equations are valid for any n link manipulator system with any kind of joints connected in any sequence. The equations of motion for the class of manipulators consisting of n rigid links interconnected by rotary joints are derived. A technique is outlined for reducing the system of equations to eliminate contraint torques. The linearized dynamics equations for an n link manipulator system are derived. The general n link linearized equations are then applied to a two link configuration. The coordinated rate control algorithm used to compute individual joint rates when given end effector rates is described. A short discussion of simulation methodologies is presented.

Lowrie, J. W.↗
NASA NTRS1996

On the Retrieval of Lightning Radio Sources from Time-of-Arrival Data

We examine the problem of retrieving three-dimensional lightning locations from radio frequency Time-Of-Arrival (TOA) measurements. Arbitrary antenna locations are considered. By judiciously differencing measurements that are related to the location of the antennas and their excitation times, the problem is converted from the initial spherical nonlinear form to a system of linear equations. In the linear formalism, the source location and time-of-occurrence is viewed geometrically as an intersection of hyperplanes in the four-dimensional Minkowski space (x,y,z,t). The linear equations are solved to obtain explicit analytic expressions for the location and time variables. Retrieval errors are not interpreted with conventional Geometrical Dilution of Precision (GDOP) arguments as discussed by Holmes and Reedy (1951), but with more recent inversion analyses considered by Twomey (1977). Measurement errors are propagated analytically so that the specific effect of these errors on the solution is clarified. The sensitivity of the solution on the number of antennas used, antenna network geometry, source position, and measurement differencing schemes are discussed in terms of the eigenvalues of the linear system.

Koshak, William J.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the loop unrolling technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗