SEARCH · Engineering Papers

Results for “Matrix 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Preconditioners for the spectral multigrid method

The systems of algebraic equations which arise from spectral discretizations of elliptic equations are full and direct solutions of them are rarely feasible. Iterative methods are an attractive alternative because Fourier transform techniques enable the discrete matrix-vector products to be computed with nearly the same efficiency as is possible for corresponding but sparse finite difference discretizations. For realistic Dirichlet problem preconditioning is essential for acceptable convergence rates. A brief description of Chebyshev spectral approximations and spectral multigrid methods for elliptic problems is given. A survey of preconditioners for Dirichlet problems based on second-order finite difference methods is made. New preconditioning techniques based on higher order finite differences and on the spectral matrix itself are presented. The preconditioners are analyzed in terms of their spectra and numerical examples are presented.

Phillips, T. N.↗

NASA NTRS1987

An optimal resolved rate law for kinematically redundant manipulators

The resolved rate law for a manipulator provides the instantaneous joint rates required to satisfy a given instantaneous hand motion. When the joint space has more degrees of freedom than the task space, the manipulator is kinematically redundant and the kinematic rate equations are underdetermined. These equations can be locally optimized, but the resulting pseudo-inverse solution has been found to cause large joint rates in some cases. A weighting matrix in the locally optimized (pseudo-inverse) solution is dynamically adjusted to control the joint motion as desired. Joint reach limit avoidance is demonstrated in a kinematically redundant planar arm model. The treatment is applicable to redundant manipulators with any number of revolute joints and to non-planar manipulators.

Bourgeois, B. J.↗

NASA NTRS1992

Non-LTE radiative transfer with lambda-acceleration - Convergence properties using exact full and diagonal lambda-operators

We investigate the convergence properties of Lambda-acceleration methods for non-LTE radiative transfer problems in planar and spherical geometry. Matrix elements of the 'exact' A-operator are used to accelerate convergence to a solution in which both the radiative transfer and atomic rate equations are simultaneously satisfied. Convergence properties of two-level and multilevel atomic systems are investigated for methods using: (1) the complete Lambda-operator, and (2) the diagonal of the Lambda-operator. We find that the convergence properties for the method utilizing the complete Lambda-operator are significantly better than those of the diagonal Lambda-operator method, often reducing the number of iterations needed for convergence by a factor of between two and seven. However, the overall computational time required for large scale calculations - that is, those with many atomic levels and spatial zones - is typically a factor of a few larger for the complete Lambda-operator method, suggesting that the approach should be best applied to problems in which convergence is especially difficult.

Macfarlane, J. J.↗

NASA NTRS1994

An efficient three-dimensional Poisson solver for SIMD high-performance-computing architectures

We present an algorithm that solves the three-dimensional Poisson equation on a cylindrical grid. The technique uses a finite-difference scheme with operator splitting. This splitting maps the banded structure of the operator matrix into a two-dimensional set of tridiagonal matrices, which are then solved in parallel. Our algorithm couples FFT techniques with the well-known ADI (Alternating Direction Implicit) method for solving Elliptic PDE's, and the implementation is extremely well suited for a massively parallel environment like the SIMD architecture of the MasPar MP-1. Due to the highly recursive nature of our problem, we believe that our method is highly efficient, as it avoids excessive interprocessor communication.

Cohl, H.↗

NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 2: Solutions and numerical results

This paper discusses an incremental-iterative nonlinear solution technique for solving the nonlinear finite element equations of the superconducting state of a superconductor. The untreated equations are highly ill-conditioned and are impossible to solve within the typical 16-place double precision supplied by most computers. A combination of matrix scaling and mesh grading techniques is used to reduce the condition number of the tangent stiffness matrix and increase the accuracy of the current carrying boundary layer representation. Numerical results for a one-dimensional model of a time-independent superconductor treated by the Ginzburg-Landau model are presented and discussed. The computed solutions clearly display the Meissner effect of magnetic field expulsion from the central region of the superconductor. These results are compared to the physics of a low-viscosity fluid problem. From this analogy, a physical argument is advanced about the macroscopic behavior of superconductors.

Schuler, James J.↗

NASA NTRS2000

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗

NASA NTRS2011

Iterative Strain-Gage Balance Calibration Data Analysis for Extended Independent Variable Sets

A new method was developed that makes it possible to use an extended set of independent calibration variables for an iterative analysis of wind tunnel strain gage balance calibration data. The new method permits the application of the iterative analysis method whenever the total number of balance loads and other independent calibration variables is greater than the total number of measured strain gage outputs. Iteration equations used by the iterative analysis method have the limitation that the number of independent and dependent variables must match. The new method circumvents this limitation. It simply adds a missing dependent variable to the original data set by using an additional independent variable also as an additional dependent variable. Then, the desired solution of the regression analysis problem can be obtained that fits each gage output as a function of both the original and additional independent calibration variables. The final regression coefficients can be converted to data reduction matrix coefficients because the missing dependent variables were added to the data set without changing the regression analysis result for each gage output. Therefore, the new method still supports the application of the two load iteration equation choices that the iterative method traditionally uses for the prediction of balance loads during a wind tunnel test. An example is discussed in the paper that illustrates the application of the new method to a realistic simulation of temperature dependent calibration data set of a six component balance.

Ulbrich, Norbert Manfred↗

NASA NTRS1986

Preconditioning matrices for Chebyshev derivative operators

The problem of preconditioning the matrices arising from pseudo-spectral Chebyshev approximations of first order operators is considered in both one and two dimensions. In one dimension a preconditioner represented by a full matrix which leads to preconditioned eigenvalues that are real, positive, and lie between 1 and pi/2, is already available. Since there are cases in which it is not computationally convenient to work with such a preconditioner, a large number of preconditioners were studied which were more sparse (in particular three and four diagonal matrices). The eigenvalues of such preconditioned matrices are compared. The results were applied to the problem of finding the steady state solution to an equation of the type u sub t = u sub x + f, where the Chebyshev collocation is used for the spatial variable and time discretization is performed by the Richardson method. In two dimensions different preconditioners are proposed for the matrix which arises from the pseudo-spectral discretization of the steady state problem. Results are given for the CPU time and the number of iterations using a Richardson iteration method for the unpreconditioned and preconditioned cases.

Rothman, Ernest E.↗

NASA NTRS1989

A new way to interpret the DIRAC equation in a non-Riemannian manifold

The idea of internal mass terms introduced in ref. (1), is shown not to be an appropriate hypothesis when it is placed in connection with the components of the generalized (matrix) vierbeins being proportional to the Riemannian (gravitational) vierbeins. It would result in an undesirable canceling of the Electromagnetic and the Yang-Mills components in the generalized metric. Another hypothesis is introduced where the wave function psi is Taylor expanded in a small parameter p.

Marques-Bonham, Sirley↗

NASA NTRS1973

Interior radiances in optically deep absorbing media. 1: Exact solutions for one-dimensional model

The exact solutions are obtained for a one dimensional model of a scattering and absorbing medium. The results are given for both the reflected and transmitted radiance for any arbitrary surface albedo as well as for the interior radiance. These same quantities are calculated by the matrix operator method. The relative error of the solutions is obtained by comparison with the exact solutions as well as by an error analysis of the equations. The importance of an accurate starting value for the reflection and transmission operators is shown. A fourth order Runge-Kutta method can be used to solve the differential equations satisfied by these operators in order to obtain such accurate starting values.

Kattawar, G. W.↗

NASA NTRS1983

Buckling of imperfect periodic lattice structures

A simplified buckling analysis is presented for a family of periodic lattice structures such as those proposed for large space structures. A transcendental 6 x 6 matrix of eigenvalues is shown to be sufficient for modeling buckling behavior because member stiffnesses are based on an exact solution of the beam-column equation. Exact stiffnesses are derived for a curved member, thus allowing modeling of imperfect lattice structures. Comparisons of predictions of the lattice model with those available from shell and beam theory underscore the inaccuracies introduced by treating the lattice structure as a continuum. Sample calculations are provided for an isogrid cylinder and a three element double-laced truss.

Anderson, M. S.↗

NASA NTRS1987

A coupled-mode theory for multiwaveguide systems satisfying the reciprocity theorem and power conservation

Two sets of coupled-mode equations for multiwaveguide systems are derived using a generalized reciprocity relation; one set for a lossless system, and the other for a general lossy or lossless system. The second set of equations also reduces to those of the first set in the lossless case under the condition that the transverse field components are chosen to be real. Analytical relations between the coupling coefficients are shown and applied to the coupling of mode equations. It is shown analytically that these results satisfy exactly both the reciprocity theorem and power conservation. New orthogonal relations between the supermodes are derived in matrix form, with the overlap integrals taken into account.

Chuang, Shun-Lien↗

NASA NTRS1988

On discrete inner-outer and spectral factorizations

Reliable algorithms are developed to perform inner-outer, coprime, and spectral factorizations for discrete FDLTI systems. It is shown that the discrete algebraic Riccati equation plays an important role in obtaining state-space representations for all key factorizations. The implementation of algorithms can be carried out efficiently using real matrix operations.

Chu, Cheng-Chih↗

NASA NTRS1989

A Transient Response Method for Linear Coupled Substructures

A method is presented for determining the transient response of a discrete coordinate model of a linear structural system composed of substructures. The method is applicable to systems consisting of any number of substructures, both determinate and indeterminate interface boundaries, and any topological arrangement of the substructures. The method is simple to implement from a computational point of view because the equations of motion of each of the substructures are solved independently, and the interface boundary compatibility conditions are enforced at each integration time step by a matrix multiplication. The method is demonstrated for a structural system consisting of two beam segments and acted upon by a time dependent force. The numerical results from the demonstration problem validates the accuracy of the method. The application of this method to structural systems with changing interface boundary conditions between substructures is discussed.

Admire, J. R.↗

NASA NTRS1991

Comparing direct and iterative equation solvers in a large structural analysis software system

Two direct Choleski equation solvers and two iterative preconditioned conjugate gradient (PCG) equation solvers used in a large structural analysis software system are described. The two direct solvers are implementations of the Choleski method for variable-band matrix storage and sparse matrix storage. The two iterative PCG solvers include the Jacobi conjugate gradient method and an incomplete Choleski conjugate gradient method. The performance of the direct and iterative solvers is compared by solving several representative structural analysis problems. Some key factors affecting the performance of the iterative solvers relative to the direct solvers are identified.

Poole, E. L.↗

NASA NTRS1992

An independent refinement and integration procedure in multiregion finite element analysis

An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.

Krishnamurthy, T.↗

NASA NTRS1993

On the cross-stream spectral method for the Orr-Sommerfeld equation

Cross-stream models are defined as solutions to the Orr-Sommerfeld equation which are propagating normal to the flow direction. These models are utilized as a basis for a Hilbert space to approximate the spectrum of the Orr-Sommerfeld equation with plane Poiseuille flow. The cross-stream basis leads to a standard eigenvalue problem for the frequencies of Poiseuille flow instability waves. The coefficient matrix in the eigenvalue problem is shown to be the sum of a real matrix and a negative-imaginary diagonal matrix which represents the frequencies of the cross-stream modes. The real coefficient matrix is shown to approach a Toeplitz matrix when the row and column indices are large. The Toeplitz matrix is diagonally dominant, and the diagonal elements vary inversely in magnitude with diagonal position. The Poiseuille flow eigenvalues are shown to lie within Gersgorin disks with radii bounded by the product of the average flow speed and the axial wavenumber. It is shown that the eigenvalues approach the Gersgorin disk centers when the mode index is large, so that the method may be used to compute spectra with an essentially unlimited number of elements. When the mode index is large, the real part of the eigenvalue is the product of the axial wavenumber and the average flow speed, and the imaginary part of the eigen value is identical to the corresponding cross-stream mode frequency. The cross-stream method is numerically well-conditioned in comparison to Chebyshev based methods, providing equivalent accuracy for small mode indices and superior accuracy for large indices.

Zorumski, William E.↗

NASA NTRS1996

MSC/NASTRAN DMAP Alter Used for Closed-Form Static Analysis With Inertia Relief and Displacement-Dependent Loads

Solving for the displacements of free-free coupled systems acted upon by static loads is a common task in the aerospace industry. Often, these problems are solved by static analysis with inertia relief. This technique allows for a free-free static analysis by balancing the applied loads with the inertia loads generated by the applied loads. For some engineering applications, the displacements of the free-free coupled system induce additional static loads. Hence, the applied loads are equal to the original loads plus the displacement-dependent loads. A launch vehicle being acted upon by an aerodynamic loading can have such applied loads. The final displacements of such systems are commonly determined with iterative solution techniques. Unfortunately, these techniques can be time consuming and labor intensive. Because the coupled system equations for free-free systems with displacement-dependent loads can be written in closed form, it is advantageous to solve for the displacements in this manner. Implementing closed-form equations in static analysis with inertia relief is analogous to implementing transfer functions in dynamic analysis. An MSC/NASTRAN (MacNeal-Schwendler Corporation/NASA Structural Analysis) DMAP (Direct Matrix Abstraction Program) Alter was used to include displacement-dependent loads in static analysis with inertia relief. It efficiently solved a common aerospace problem that typically has been solved with an iterative technique.

Source record↗