Search NASASearch

SEARCH · Search NASA

Results for “Matrix differential equation”

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 91 records · Page 5

A manual for PARTI runtime primitives

Primitives are presented that are designed to help users efficiently program irregular problems (e.g., unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial differential equations solvers) on distributed memory machines. These primitives are also designed for use in compilers for distributed memory multiprocessors. Communications patterns are captured at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry

A manual for PARTI runtime primitives, revision 1

Primitives are presented that are designed to help users efficiently program irregular problems (e.g., unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial differential equations solvers) on distributed memory machines. These primitives are also designed for use in compilers for distributed memory multiprocessors. Communications patterns are captured at runtime, and the appropriate send and receive messages are automatically generated.

Das, Raja

A computer program to generate equations of motion matrices, L217 (EOM). Volume 2: Supplemental system design and maintenance document

The equations of motion program L217 (EOM) is described. The program formulates the matrix coefficients for a second order linear differential equation which describes the motion of an airplane relative to its level equilibrium flight condition. Aerodynamic data from FLEXSTAB or Doublet Lattice (L216) programs are used to derive the equations for quasi-steady or complete unsteady aerodynamics.

Clemmons, R. E.

Orbit improvement from satellite imaging data obtainable from the Mariner/Jupiter/Saturn 1977 missions

Equations of motion are established for a dynamic system in which a spacecraft flies close to and interacts with an outer planet and one or more of its satellites. For the computation of the state and mass partials needed in the orbit corrections, a set of variational equations is derived. These assume a notably compact form through the introduction of a matrix operator D. The above system of differential equations is integrated numerically on a computer.

Aksnes, K.

Modeling the missile-launch tube problem in DYSCO

DYSCO is a versatile, general purpose dynamic analysis program which assembles equations and solves dynamics problems. The executive manages a library of technology modules which contain routines that compute the matrix coefficients of the second order ordinary differential equations of the components. The executive performs the coupling of the equations of the components and manages the solution of the coupled equations. Any new component representation may be added to the library if, given the state vector, a FORTRAN program can be written to compute M, C, K, and F. The problem described demonstrates the generality of this statement.

Berman, Alex

An invariant asymptotic formula for solutions of second-order linear ODE's

An invariant-matrix technique for the approximate solution of second-order ordinary differential equations (ODEs) of form y-double-prime = phi(x)y is developed analytically and demonstrated. A set of linear transformations for the companion matrix differential system is proposed; the diagonalization procedure employed in the final stage of the asymptotic decomposition is explained; and a scalar formulation of solutions for the ODEs is obtained. Several typical ODEs are analyzed, and it is shown that the Liouville-Green or WKB approximation is a special case of the present formula, which provides an approximation which is valid for the entire interval (0, infinity).

Gingold, H.

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.

A revised version of the transfer matrix method to analyze one-dimensional structures

A new and general method to analyze both free and forced vibration characteristics of one-dimensional structures is discussed in this paper. This scheme links for the first time the classical transfer matrix method with the recently developed integrating matrix technique to integrate systems of differential equations. Two alternative approaches to the problem are presented. The first is based upon the lumped parameter model to account for the inertia properties of the structure. The second releases that constraint allowing a more precise description of the physical system. The free vibration of a straight uniform beam under different support conditions is analyzed to test the accuracy of the two models. Finally some results for the free vibration of a 12th order system representing a curved, rotating beam prove that the present method is conveniently extended to more complicated structural dynamics problems.

Nitzsche, F.

A Continuous Square Root in Formation Filter-Swoother with Discrete Data Update

A differential equation for the square root information matrix is derived and adapted to the problems of filtering and smoothing. The resulting continuous square root information filter (SRIF) performs the mapping of state and process noise by numerical integration of the SRIF matrix and admits data via a discrete least square update.

SRIF mapping

Two-demensional analysis of heat and mass transfer in porous media using the strongly implicit procedure

Numerical results of the heat and mass transfer in a porous matrix are presented. The coupled, nonlinear partial differential equations describing this physical phenomenon are solved in finite difference form for two dimensions, using a new iterative technique (the strongly implicit procedure). The influence of the external environment conditions (heating and pressure) is shown to produce two-dimensional flow in the porous matrix. Typical fluid and solid temperature distributions in the porous matrix and internal pressure distributions are presented.

Curry, D. M.

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.

Unitary-matrix models as exactly solvable string theories

Exact differential equations are presently found for the scaling functions of models of unitary matrices which are solved in a double-scaling limit, using orthogonal polynomials on a circle. For the case of the simplest, k = 1 model, the Painleve II equation with constant 0 is obtained; possible nonperturbative phase transitions exist for these models. Equations are presented for k = 2 and 3, and discussed with a view to asymptotic behavior.

Periwal, Vipul

An intrinsically n-dimensional generalized flux vector splitting implicit finite element Euler algorithm

A generalized flux-vector splitting implicit Galerkin finite-element algorithm for the Euler equations in curvilinear coordinates for ideal and reacting gases is derived. For an arbitrary equation of state, the curvilinear-coordinate flux vector is split in kinematic and kinetic components, and the associated jacobian matrix eigenvalues explicitly depend on the metric data. After directional semidiscretization, the terminal ordinary differential-equation system is solved via a nonlinearly stable implicit Runge-Kutta scheme in concert with an accurate tensor matrix product factorization. The results for selected two-dimensional supersonic and axisymmetric hypersonic flows validate the algorithm and verify its robustness for curvilinear-coordinate computations. The evolution towards a steady state is achieved for large Courant numbers without indication of numerical instabilities.

Iannelli, G. S.

Integrating Matrix Method for Determining the Natural Vibrations of a Rotating, Unsymmetrical Beam with Application to Twisted Propeller Blades

A numerical method is Presented for determining the natural vibration frequencies, and the corresponding mode shapes, of a rotating cantilever beam which has a nonuniform, unsymmetrical cross section. Two coupled fourth-order differential equations of motion with variable coefficients are derived which govern the motion of such a beam having deformations in two directions. Through the development and utilization of the integrating matrix, the solution of the differential equations is obtained in the form of an eigenvalue problem. The solutions to the eigenvalue problem are determined by an iteration method based upon a special orthogonality relationship which is derived. Numerical examples, including an application to a twisted propeller blade, are presented with the results of the integrating matrix solutions being compared to exact solutions and experimental data.

Hunter, William F.