Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Analytic prediction of airplane equilibrium spin characteristics

The nonlinear equations of motion are solved algebraically for conditions for which an airplane is in an equilibrium spin. Constrained minimization techniques are employed in obtaining the solution. Linear characteristics of the airplane about the equilibrium points are also presented and their significance in identifying the stability characteristics of the equilibrium points is discussed. Computer time requirements are small making the method appear potentially applicable in airplane design. Results are obtained for several configurations and are compared with other analytic-numerical methods employed in spin prediction. Correlation with experimental results is discussed for one configuration for which a rather extensive data base was available. A need is indicated for higher Reynolds number data taken under conditions which more accurately simulate a spin.

Adams, W. M., Jr.↗
NASA NTRS1975

Some properties of extremal solutions of zero-sum differential games

Necessary conditions for infinite-time linear-quadratic two-player differential games with perfect state information are examined in the light of Willems' work on the algebraic Riccati equation. Necessary conditions for the nonzero-sum game do not always reduce precisely to those of the zero-sum game, when the weighting parameters of the cost index are taken as those of a zero-sum game. Formulation of zero-sum game problems in extended L2-spaces reveals the existence of nonstabilizing but infimal controls, and explicit frequency-domain conditions for existence of solutions and stabilizing solutions are established through the algebraic Riccati equation (ARE).

Johnson, T. L.↗
NASA NTRS1983

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity is studied in detail via the eigenstrain approach. A complete formulation and a treatment of both types of eigenstrains for equivalence between the inhomogeneity problem and the inclusion problem are given. This approach is shown to be mathematically identical to other approaches such as the direct volume integral formulation. Expanding the eigenstrains and applied strains in the polynomial form in the position vector and satisfying the equivalence conditions at every point, the governing simultaneous algebraic equations for the unknown coefficients in the eigenstrain expansion are derived. The elastodynamic field outside an ellipsoidal inhomogeneity in a linear elastic isotropic medium is given as an example. The angular and frequency dependence of the induced displacement field, as well as the differential and total cross sections are formally given in series expansion form for the case of uniformly distributed eigenstrains.

Fu, L. S.↗
NASA NTRS1984

Nonlinear evolution of magnetopause tearing modes

Since the magnetosheath plasma is highly turbulent, reconnection at the dayside magnetopause is likely to be temporally unsteady. The tearing mode can be viewed as a model for the unsteady development of a reconnecting magnetic topology. Magnetopause tearing occurs in the guide-field limit and has a wave packet spatial structure in the east-west direction. This paper solves for the nonlinear evolution of a single wavelength guide-field tearing mode including the effects of finite transit time on the Landau resonant electrons. Short wavelength modes evolve algebraically in time with perturbation amplitudes proportional to t-squared. Long wavelength modes are fully nonlinear, and the amplitude grows linearly in time.

Coroniti, F. V.↗
NASA NTRS1988

Recursive mass matrix factorization and inversion: An operator approach to open- and closed-chain multibody dynamics

This report advances a linear operator approach for analyzing the dynamics of systems of joint-connected rigid bodies.It is established that the mass matrix M for such a system can be factored as M=(I+H phi L)D(I+H phi L) sup T. This yields an immediate inversion M sup -1=(I-H psi L) sup T D sup -1 (I-H psi L), where H and phi are given by known link geometric parameters, and L, psi and D are obtained recursively by a spatial discrete-step Kalman filter and by the corresponding Riccati equation associated with this filter. The factors (I+H phi L) and (I-H psi L) are lower triangular matrices which are inverses of each other, and D is a diagonal matrix. This factorization and inversion of the mass matrix leads to recursive algortihms for forward dynamics based on spatially recursive filtering and smoothing. The primary motivation for advancing the operator approach is to provide a better means to formulate, analyze and understand spatial recursions in multibody dynamics. This is achieved because the linear operator notation allows manipulation of the equations of motion using a very high-level analytical framework (a spatial operator algebra) that is easy to understand and use. Detailed lower-level recursive algorithms can readily be obtained for inspection from the expressions involving spatial operators. The report consists of two main sections. In Part 1, the problem of serial chain manipulators is analyzed and solved. Extensions to a closed-chain system formed by multiple manipulators moving a common task object are contained in Part 2. To retain ease of exposition in the report, only these two types of multibody systems are considered. However, the same methods can be easily applied to arbitrary multibody systems formed by a collection of joint-connected regid bodies.

Rodriguez, G.↗
NASA NTRS1991

Effect of periodic accelerations on interface stability in a multilayered fluid configuration

The increasing number of research opportunities in a microgravity environment will benefit not only fundamental studies in fluid dynamics, but also technological applications such as those involving materials processing. In particular, fluid configurations that involve fluid-fluid interfaces would occur in a variety of experimental investigations. This work investigates the stability of a configuration involving fluid-fluid interfaces in the presence of a time-dependent (periodic) forcing. The fluid configuration is multilayered and infinite in extent. The analysis is linear and inviscid, and the acceleration vector is oriented perpendicular to each interface. A Floquent analysis is employed, and the resulting algebraic eigensystem is truncated. Nondimensional parameters appear in the algebraic system. A numerical study is performed to elucidate the regions of instability and the effect of parameter variation on the fluid configuration stability.

Lyell, M. J.↗
NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗
NASA NTRS1997

Complementary Reliability-Based Decodings of Binary Linear Block Codes

This correspondence presents a hybrid reliability-based decoding algorithm which combines the reprocessing method based on the most reliable basis and a generalized Chase-type algebraic decoder based on the least reliable positions. It is shown that reprocessing with a simple additional algebraic decoding effort achieves significant coding gain. For long codes, the order of reprocessing required to achieve asymptotic optimum error performance is reduced by approximately 1/3. This significantly reduces the computational complexity, especially for long codes. Also, a more efficient criterion for stopping the decoding process is derived based on the knowledge of the algebraic decoding solution.

Fossorier, Marc P. C.↗
NASA NTRS1975

Controllability in nonlinear systems

An explicit expression for the reachable set is obtained for a class of nonlinear systems. This class is described by a chain condition on the Lie algebra of vector fields associated with each nonlinear system. These ideas are used to obtain a generalization of a controllability result for linear systems in the case where multiplicative controls are present.

Hirschorn, R. M.↗
NASA NTRS1996

Generalized Bezout's Theorem and its applications in coding theory

This paper presents a generalized Bezout theorem which can be used to determine a tighter lower bound of the number of distinct points of intersection of two or more curves for a large class of plane curves. A new approach to determine a lower bound on the minimum distance (and also the generalized Hamming weights) for algebraic-geometric codes defined from a class of plane curves is introduced, based on the generalized Bezout theorem. Examples of more efficient linear codes are constructed using the generalized Bezout theorem and the new approach. For d = 4, the linear codes constructed by the new construction are better than or equal to the known linear codes. For d greater than 5, these new codes are better than the known codes. The Klein code over GF(2(sup 3)) is also constructed.

Berg, Gene A.↗
NASA NTRSDate not supplied

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗
NASA NTRS2001

Glovebox Integrated Microgravity Isolation Technology (g-LIMIT): A Linearized State-Space Model

Vibration acceleration levels on large space platforms exceed the requirements of many space experiments. The Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) is being built by the NASA Marshall Space Flight Center to attenuate these disturbances to acceptable levels. G-LIMIT uses Lorentz (voice-coil) magnetic actuators to levitate and isolate payloads at the individual experiment/sub-experiment (versus rack) level. Payload acceleration, relative position, and relative orientation measurements are fed to a state-space controller. The controller, in turn, determines the actuator Currents needed for effective experiment isolation. This paper presents the development of an algebraic, state-space model of g-LIMIT, in a form suitable for optimal controller design. The equations are first derived using Newton's Second Law directly, then simplified to a linear form for the purpose of controller design.

Hampton, R. David↗
NASA NTRS1973

Suboptimal design of a class of nonlinear controllers.

A systematic design procedure is presented for suboptimal nonlinear controller synthesis. Such a control yields faster response characteristics than are normally possible with a linear controller. The technique is geared to single-input linear systems where excursions of one state variable are of more importance than others. Design steps are formulated in terms of readily applicable algebraic equations. The proposed method is applied to several design examples and the advantages are shown by comparison with the results obtained from existing techniques.

Rouse, W. B.↗
NASA NTRS1989

Algebraic multigrid methods applied to problems in computational structural mechanics

The development of algebraic multigrid (AMG) methods and their application to certain problems in structural mechanics are described with emphasis on two- and three-dimensional linear elasticity equations and the 'jacket problems' (three-dimensional beam structures). Various possible extensions of AMG are also described. The basic idea of AMG is to develop the discretization sequence based on the target matrix and not the differential equation. Therefore, the matrix is analyzed for certain dependencies that permit the proper construction of coarser matrices and attendant transfer operators. In this manner, AMG appears to be adaptable to structural analysis applications.

Mccormick, Steve↗
NASA NTRS1976

Computation and sensitivity considerations of the Ricatti equation in spacecraft CMG steering laws

The steady-state matrix Ricatti equation is studied with respect to sensitivity and computational burden for a linear regulator controller in spacecraft. The class of systems studied use control moment gyros (CMGs) as the torque driver, assuming that the CMGs are free of gimbal stops. Different solution techniques are evaluated relative to computation time and accuracy to determine which solution method is most acceptable with respect to the prescribed requirements. It is shown that it is not feasible to eliminate off-diagonal Ricatti equation terms in the relevent updata equation due to their high sensitivity, that Potter's algebraic technique requires the least computation time of the methods studied, and that the simple Euler method yields the best results when accuracy and computation time are equally weighted. If an optimal linear regulator formulation is employed, a full Ricatti equation solution of n(n + 1)2 elements of the symmetric positive definite gain matrix which is the solution to the matrix Ricatti equation must be computed at each update cycle.

Colburn, B. K.↗
NASA NTRS1982

Controller scheduling - A possible algebraic viewpoint

In the applications, one common way to design a control system for a nonlinear plant is to localize its behavior along lines of operation specified by the plant manufacturer, to develop linear multi-variable controls for these localizations, and to schedule those controls with key plant variables which vary smoothly along operating lines. An important part of practical design lore, the art of controller scheduling has received little modern attention from the conceptual point of view. This paper describes four basic types of scheduling questions and outlines some of the theoretical issues associated with them. Schedules are considered in terms of state equations; however, some relations with the input/output description are discussed, together with an analysis of the effects on the overall configuration of approximations made to the individual subsystems.

Sain, M.↗
NASA NTRS1990

GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 2: User's manual and program listing

An efficient computer program, called GRID2D/3D, was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. The theory and method used in GRID2D/3D is described.

Bailey, R. T.↗