Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Harmonic analysis of spacecraft power systems using a personal computer

The effects that nonlinear devices such as ac/dc converters, HVDC transmission links, and motor drives have on spacecraft power systems are discussed. The nonsinusoidal currents, along with the corresponding voltages, are calculated by a harmonic power flow which decouples and solves for each harmonic component individually using an iterative Newton-Raphson algorithm. The sparsity of the harmonic equations and the overall Jacobian matrix is used to an advantage in terms of saving computer memory space and in terms of reducing computation time. The algorithm could also be modified to analyze each harmonic separately instead of all at the same time.

Williamson, Frank↗
NASA NTRS1991

Vibration suppression using a constrained rate-feedback-threshold control strategy

A finite time, minimum force rate-feedback-threshold controller is developed to bring a system with or without known external disturbances back into an 'allowable' state bound in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The optimal control is defined by a two-point boundary value problem coupled to a set of definite integral constraints. Quasi-closed form solutions are derived which replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples are provided which demonstrate the threshold control technique and compare the quasi-closed form solutions with numerical and MACSYMA generated exact solutions.

Zimmerman, D. C.↗
NASA NTRS1991

Predictions of the critical strain for matrix cracking of ceramic matrix composites

The critical strain for matrix cracking of ceramic matrix composites has been studied; emphasis is placed on the effects of the fiber/matrix debonding length and interfacial debonding energy. Based on a modified shear-lag model, stresses in the fiber and matrix have been found for both the bonded and debonded regions. An energy-balance approach, based upon the stress field, is then adopted to evaluate the critical strain for matrix cracking. From the general equation for the critical strain, close form solutions have been deduced for two limiting cases: complete debonding and perfect bonding. Numerical solutions are given for the cases of nonzero debonding energy and partial fiber debonding. The results show that the interfacial debonding energy, which has been ignored by most of the investigators, is an important factor in determining both the critical strain and the debonding length.

Kuo, Wen-Shyong↗
NASA NTRS1999

Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation

In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.

Kouatchou, Jules↗
NASA NTRS2006

Analyzing Aeroelastic Stability of a Tilt-Rotor Aircraft

Proprotor Aeroelastic Stability Analysis, now at version 4.5 (PASTA 4.5), is a FORTRAN computer program for analyzing the aeroelastic stability of a tiltrotor aircraft in the airplane mode of flight. The program employs a 10-degree- of-freedom (DOF), discrete-coordinate, linear mathematical model of a rotor with three or more blades and its drive system coupled to a 10-DOF modal model of an airframe. The user can select which DOFs are included in the analysis. Quasi-steady strip-theory aerodynamics is employed for the aerodynamic loads on the blades, a quasi-steady representation is employed for the aerodynamic loads acting on the vibrational modes of the airframe, and a stability-derivative approach is used for the aerodynamics associated with the rigid-body DOFs of the airframe. Blade parameters that vary with the blade collective pitch can be obtained by interpolation from a user-defined table. Stability is determined by examining the eigenvalues that are obtained by solving the coupled equations of motions as a matrix eigenvalue problem. Notwithstanding the relative simplicity of its mathematical foundation, PASTA 4.5 and its predecessors have played key roles in a number of engineering investigations over the years.

Kvaternil, Raymond G.↗
NASA NTRS2016

Scattering of Gaussian Beams by Disordered Particulate Media

A frequently observed characteristic of electromagnetic scattering by a disordered particulate medium is the absence of pronounced speckles in angular patterns of the scattered light. It is known that such diffuse speckle-free scattering patterns can be caused by averaging over randomly changing particle positions and/or over a finite spectral range. To get further insight into the possible physical causes of the absence of speckles, we use the numerically exact superposition T-matrix solver of the Maxwell equations and analyze the scattering of plane-wave and Gaussian beams by representative multi-sphere groups. We show that phase and amplitude variations across an incident Gaussian beam do not serve to extinguish the pronounced speckle pattern typical of plane-wave illumination of a fixed multi-particle group. Averaging over random particle positions and/or over a finite spectral range is still required to generate the classical diffuse speckle-free regime.

radiative transfer theory↗
NASA NTRS1983

Riccati group invariants of linear hamiltonian systems

The action of the Riccati group on the Riccati differential equation is associated with the action of a subgroup of the symplectic group on a set of hamiltonian matrices. Within this framework various sets of canonical forms are developed for the matrix coefficients of the Riccati differential equation. The canonical forms presented are valid for arbitrary Kronecker indices, and it is shown that the Kronecker indices are invariants for this group action. These canonical forms are useful for studying problems arising in the areas of optimal decentralized control and the spectral theory of optimal control problems.

Garzia, M. R.↗
NASA NTRS1987

Strategies for vectorizing the sparse matrix vector product on the CRAY XMP, CRAY 2, and CYBER 205

Large, randomly sparse matrix vector products are important in a number of applications in computational chemistry, such as matrix diagonalization and the solution of simultaneous equations. Vectorization of this process is considered for the CRAY XMP, CRAY 2, and CYBER 205, using a matrix of dimension of 20,000 with from 1 percent to 6 percent nonzeros. Efficient scatter/gather capabilities add coding flexibility and yield significant improvements in performance. For the CYBER 205, it is shown that minor changes in the IO can reduce the CPU time by a factor of 50. Similar changes in the CRAY codes make a far smaller improvement.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

Diagonalized Lagrangian Robot Dynamics

A diagonal equation for robot dynamics is developed by combining recent mass matrix factorization results with classical Lagrangian mechanics.

working joint moments total joint rotational rates↗
NASA NTRS1972

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.↗
NASA NTRS1999

Strain Rate Dependent Modeling of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Strain rate dependent inelastic constitutive equations have been developed to model the polymer matrix, and have been incorporated into a micromechanics approach to analyze polymer matrix composites. The Hashin failure criterion has been implemented within the micromechanics results to predict ply failure strengths. The deformation model has been implemented within LS-DYNA, a commercially available transient dynamic finite element code. The deformation response and ply failure stresses for the representative polymer matrix composite AS4/PEEK have been predicted for a variety of fiber orientations and strain rates. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗
DOE OSTI2024

Impact of Reordering on the LU Factorization Performance of Bordered Block-Diagonal Sparse Matrix

Power engineers rely on computer-based simulation tools to assess grid performance and ensure security. At the core of these tools are solvers for sparse linear equations. When transformed into a bordered block-diagonal (BBD) structure, part of the sparse linear equation solving can be parallelized. This work focuses on using the Schur-complement-based method for LU factorization on BBD matrices, specifically, Jacobian matrices from large-scale systems. Our findings show that the natural ordering method outperforms the default ordering method in computational performance for each block of the BBD matrix. This observation is validated using synthetic 25k-bus and 70k-bus cases, showing a speedup of up to 38% when using natural ordering without permutation. Additionally, the impact of the number of partitions is studied, and the result shows that computational performance improves with more, smaller partitions in the BBD matrices.

BBD matrix↗
NASA NTRS1981

Research on numerical algorithms for large space structures

Numerical algorithms for analysis and design of large space structures are investigated. The sign algorithm and its application to decoupling of differential equations are presented. The generalized sign algorithm is given and its application to several problems discussed. The Laplace transforms of matrix functions and the diagonalization procedure for a finite element equation are discussed. The diagonalization of matrix polynomials is considered. The quadrature method and Laplace transforms is discussed and the identification of linear systems by the quadrature method investigated.

Denman, E. D.↗
NASA NTRS1995

Finite-element time evolution operator for the anharmonic oscillator

The finite-element approach to lattice field theory is both highly accurate (relative errors approximately 1/N(exp 2), where N is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly preserved at the lattice sites). In this talk I construct matrix elements for dynamical variables and for the time evolution operator for the anharmonic oscillator, for which the continuum Hamiltonian is H = p(exp 2)/2 + lambda q(exp 4)/4. Construction of such matrix elements does not require solving the implicit equations of motion. Low order approximations turn out to be extremely accurate. For example, the matrix element of the time evolution operator in the harmonic oscillator ground state gives a results for the anharmonic oscillator ground state energy accurate to better than 1 percent, while a two-state approximation reduces the error to less than 0.1 percent.

Milton, Kimball A.↗
NASA NTRS2007

Corrigendum: New Form of Kane's Equations of Motion for Constrained Systems

A correction to the previously published article "New Form of Kane's Equations of Motion for Constrained Systems" is presented. Misuse of the transformation matrix between time rates of change of the generalized coordinates and generalized speeds (sometimes called motion variables) resulted in a false conclusion concerning the symmetry of the generalized inertia matrix. The generalized inertia matrix (sometimes referred to as the mass matrix) is in fact symmetric and usually positive definite when one forms nonminimal Kane's equations for holonomic or simple nonholonomic systems, systems subject to nonlinear nonholonomic constraints, and holonomic or simple nonholonomic systems subject to impulsive constraints according to Refs. 1, 2, and 3, respectively. The mass matrix is of course symmetric when one forms minimal equations for holonomic or simple nonholonomic systems using Kane s method as set forth in Ref. 4.

Roithmayr, Carlos M.↗
NASA NTRS2016

Time-Dependent Stress Rupture Strength Degradation of Hi-Nicalon Fiber-Reinforced Silicon Carbide Composites at Intermediate Temperatures

The stress rupture strength of silicon carbide fiber-reinforced silicon carbide composites with a boron nitride fiber coating decreases with time within the intermediate temperature range of 700 to 950 degree Celsius. Various theories have been proposed to explain the cause of the time-dependent stress rupture strength. The objective of this paper is to investigate the relative significance of the various theories for the time-dependent strength of silicon carbide fiber-reinforced silicon carbide composites. This is achieved through the development of a numerically based progressive failure analysis routine and through the application of the routine to simulate the composite stress rupture tests. The progressive failure routine is a time-marching routine with an iterative loop between a probability of fiber survival equation and a force equilibrium equation within each time step. Failure of the composite is assumed to initiate near a matrix crack and the progression of fiber failures occurs by global load sharing. The probability of survival equation is derived from consideration of the strength of ceramic fibers with randomly occurring and slow growing flaws as well as the mechanical interaction between the fibers and matrix near a matrix crack. The force equilibrium equation follows from the global load sharing presumption. The results of progressive failure analyses of the composite tests suggest that the relationship between time and stress-rupture strength is attributed almost entirely to the slow flaw growth within the fibers. Although other mechanisms may be present, they appear to have only a minor influence on the observed time-dependent behavior.

Progressive Failure Analysis↗
NASA NTRS2015

A Numerical Solution Routine for Investigating Oxidation-Induced Strength Degradation Mechanisms in SiC/SiC Composites

The stress rupture strength of silicon carbide fiber-reinforced silicon carbide (SiCSiC) composites with a boron nitride (BN) fiber coating decreases with time within the intermediate temperature range of 700-950 C. Various theories have been proposed to explain the cause of the time dependent stress rupture strength. Some previous authors have suggested that the observed composite strength behavior is due to the inherent time dependent strength of the fibers, which is caused by the slow growth of flaws within the fibers. Flaw growth is supposedly enabled by oxidation of free carbon at the grain boundaries. The objective of this paper is to investigate the relative significance of the various theories for the time-dependent strength of SiCSiC composites. This is achieved through the development of a numerically-based progressive failure analysis routine and through the application of the routine to simulate the composite stress rupture tests. The progressive failure routine is a time marching routine with an iterative loop between a probability of fiber survival equation and a force equilibrium equation within each time step. Failure of the composite is assumed to initiate near a matrix crack and the progression of fiber failures occurs by global load sharing. The probability of survival equation is derived from consideration of the strength of ceramic fibers with randomly occurring and slow growing flaws as well as the mechanical interaction between the fibers and matrix near a matrix crack. The force equilibrium equation follows from the global load sharing presumption. The results of progressive failure analyses of the composite tests suggest that the relationship between time and stress-rupture strength is attributed almost entirely to the slow flaw growth within the fibers. Although other mechanisms may be present, they appear to have only a minor influence on the observed time dependent behavior.

degradation↗