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 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Calculation of vaporization rates assuming various rate determining steps: Application to the resistojet

The various steps that could control the vaporization rate of a material are discussed. These steps include the actual vaporization, flow rate of matrix gas, chemical reaction, gas diffusion, and solid state diffusion. The applicable equations have been collected from diverse appropriate sources, and their use is explained. Rate equations are derived for conditions where more than one step is rate controlling. Calculations are made for two model materials: rhenium which vaporizes congruently, and tantalum carbide which vaporizes incongruently. The case of vaporization under thermal gradient conditions is also treated. The existence of a thermal gradient in the resistojet means that the vaporization rate of a material may be only one thousandth of that predicted under isothermal conditions. Calculations show that rhenium might have a 100,000 hr lifetime at temperature in a 2500 C resistojet. Tantalum carbide would have a life of only 660 sec under similar conditions.

May, C. E.↗

NASA NTRS1985

Parallel solution of closely coupled systems

An odd-even permutation and a nested dissection technique were used to circumvent the strong seriality of a system of closely coupled equations. The effect of transforming the n x n Hermitian definite positive matrix coefficient on the topology of Cholesky factors is discussed. A series of directed graphs is constructed in order to show the computational steps required for the odd-even permutation. Numerical expressions for the speed-up and efficiency of parallel N-processing techniques and sequential processing by a single computer are derived. Similar expressions are derived for the case of insufficient processing capacity. The application of the odd-even permutation to the ensemble class of computer architectures is demonstrated.

Utku, S.↗

NASA NTRS1986

Multi-color incomplete Cholesky conjugate gradient methods for vector computers

In this research, we are concerned with the solution on vector computers of linear systems of equations, Ax = b, where A is a larger, sparse symmetric positive definite matrix. We solve the system using an iterative method, the incomplete Cholesky conjugate gradient method (ICCG). We apply a multi-color strategy to obtain p-color matrices for which a block-oriented ICCG method is implemented on the CYBER 205. (A p-colored matrix is a matrix which can be partitioned into a pXp block matrix where the diagonal blocks are diagonal matrices). This algorithm, which is based on a no-fill strategy, achieves O(N/p) length vector operations in both the decomposition of A and in the forward and back solves necessary at each iteration of the method. We discuss the natural ordering of the unknowns as an ordering that minimizes the number of diagonals in the matrix and define multi-color orderings in terms of disjoint sets of the unknowns. We give necessary and sufficient conditions to determine which multi-color orderings of the unknowns correpond to p-color matrices. A performance model is given which is used both to predict execution time for ICCG methods and also to compare an ICCG method to conjugate gradient without preconditioning or another ICCG method. Results are given from runs on the CYBER 205 at NASA's Langley Research Center for four model problems.

Poole, E. L.↗

NASA NTRS1986

Differentiating matrices for arbitrarily spaced grid points

Differentiating matrices allow the numerical differentiation of functions defined at points of a discrete grid. Previous derivations of these matrices have been restricted to grids with uniformly spaced points, and the resulting derivative approximations have lacked precision, especially at endpoints. The present work derives differentiating matrices on grids with arbitrarily spaced points. It is shown that high accuracy can be achieved through use of differentiating matrices on non-uniform grids through the expedient of including 'near boundary' points. Use of the differentiating matrix as an operator in the solution of problems involving ordinary differential equations is also considered.

Lakin, W. D.↗

NASA NTRS1988

Large planar maneuvers for articulated flexible manipulators

An articulated flexible manipulator carried on a translational cart is maneuvered by an active controller to perform certain position control tasks. The nonlinear dynamics of the articulated flexible manipulator are derived and a transformation matrix is formulated to localize the nonlinearities within the inertia matrix. Then a feedback linearization scheme is introduced to linearize the dynamic equations for controller design. Through a pole placement technique, a robust controller design is obtained by properly assigning a set of closed-loop desired eigenvalues to meet performance requirements. Numerical simulations for the articulated flexible manipulators are given to demonstrate the feasibility and effectiveness of the proposed position control algorithms.

Huang, Jen-Kuang↗

NASA NTRS1988

Large planar maneuvers for articulated flexible manipulators

An articulated flexible manipulator carried on a translational cart is maneuvered by an active controller to perform certain position control tasks. The nonlinear dynamics of the articulated flexible manipulator are derived and a transformation matrix is formulated to localize the nonlinearities in the inertia matrix. Then a feeback linearization scheme is introduced to linearize the dynamic equations for controller design. Through a pole placement technique, a robust controller design is obtained by properly assigning a set of closed-loop desired eigenvalues to meet performance requirements. Numerical simulations for the articulated flexible manipulators are given to demonstrate the feasibility and effectiveness of the proposed position control algorithms.

Huang, Jen-Kuang↗

NASA NTRS1988

Active control of transient rotordynamic vibration by optimal control methods

Although considerable effort has been put into the study of steady state vibration control, there are few methods applicable to transient vibration control of rotorbearing systems. In this paper optimal control theory has been adopted to minimize rotor vibration due to sudden imbalance, e.g., blade loss. The system gain matrix is obtained by choosing the weighting matrices and solving the Riccati equation. Control forces are applied to the system via a feedback loop. A seven mass rotor system is simulated for illustration. A relationship between the number of sensors and the number of modes used in the optimal control model is investigated. Comparisons of responses are made for various configurations of modes, sensors, and actuators. Furthermore, spillover effect is examined by comparing results from collocated and noncollocated sensor configurations. Results show that shaft vibration is significantly attenuated in the closed loop system.

Palazzolo, A. B.↗

NASA NTRS1988

Symbolic generation of constitutive equations

The use of the symbolic manipulation program MACSYMA for the automatic generation of constitutive equations describing the thermomechanical behavior of complex materials such as metal-matrix composites is demonstrated. Problem-specific algorithms are developed; e.g., for partial differentiation by the chain rule, automatic generation of tree structures, special tensor analysis, and the simplification of expressions. The resulting subprogram package, SDICE, is shown to be fully applicable to both isotropic and anisotropic materials; its potential usefulness as a basis for FEM analyses is indicated.

Tan, H. Q.↗

NASA NTRS1989

A finite element approach for large motion dynamic analysis of multibody structures in space

A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Chang, Che-Wei↗

NASA NTRS1992

Symmetry properties in polarimetric remote sensing

This paper presents the relations among polarimetric backscattering coefficients from the viewpoint of symmetry groups. Symmetry of geophysical media encountered in remote sensing due to reflection, rotation, azimuthal, and centrical symmetry groups is considered for both reciprocal and nonreciprocal cases. On the basis of the invariance under symmetry transformations in the linear polarization basis, the scattering coefficients are related by a set of equations which restrict the number of independent parameters in the polarimetric covariance matrix. The properties derived under these transformations are general and apply to all scattering mechanisms in a given symmetrical configuration. The scattering coefficients calculated from theoretical models for layer random media and rough surfaces are shown to obey the derived symmetry relations. Use of symmetry properties in remote sensing of structural and environmental responses of scattering media is discussed. As a practical application, the results from this paper provide new methods for the external calibration of polarimetric radars without the deployment of man-made calibration targets.

Nghiem, S. V.↗

NASA NTRS1993

The use of linked lists in the simulation of controller-structure interaction

An algorithm for the computer simulation of large space structures under active control is considered. Linked lists are used in a matrix data structure to implement the trapezoidal rule on the system differential equations. The use of the trapezoidal rule ensures that the numerical stability is equivalent to the system stability, which is essential for this type of simulation. The sparsity of the system matrices is exploited by the linked lists, and the algorithm efficiently steps through the lists in an orderly fashion. Results of simulations on a NASA large space structure experiment are reported.

Quan, Ralph↗

NASA NTRS2012

OVERSMART Reporting Tool for Flow Computations Over Large Grid Systems

Structured grid solvers such as NASA's OVERFLOW compressible Navier-Stokes flow solver can generate large data files that contain convergence histories for flow equation residuals, turbulence model equation residuals, component forces and moments, and component relative motion dynamics variables. Most of today's large-scale problems can extend to hundreds of grids, and over 100 million grid points. However, due to the lack of efficient tools, only a small fraction of information contained in these files is analyzed. OVERSMART (OVERFLOW Solution Monitoring And Reporting Tool) provides a comprehensive report of solution convergence of flow computations over large, complex grid systems. It produces a one-page executive summary of the behavior of flow equation residuals, turbulence model equation residuals, and component forces and moments. Under the automatic option, a matrix of commonly viewed plots such as residual histograms, composite residuals, sub-iteration bar graphs, and component forces and moments is automatically generated. Specific plots required by the user can also be prescribed via a command file or a graphical user interface. Output is directed to the user s computer screen and/or to an html file for archival purposes. The current implementation has been targeted for the OVERFLOW flow solver, which is used to obtain a flow solution on structured overset grids. The OVERSMART framework allows easy extension to other flow solvers.

Kao, David L.↗

NASA NTRS1985

Short-wavelength buckling and shear failures for compression-loaded composite laminates

The short-wavelength buckling (or the microbuckling) and the interlaminar and inplane shear failures of multi-directional composite laminates loaded in uniaxial compression are investigated. A laminate model is presented that idealizes each lamina. The fibers in the lamina are modeled as a plate, and the matrix in the lamina is modeled as an elastic foundation. The out-of-plane w displacement for each plate is expressed as a trigonometric series in the half-wavelength of the mode shape for laminate short-wavelength buckling. Nonlinear strain-displacement relations are used. The model is applied to symmetric laminates having linear material behavior. The laminates are loaded in uniform end shortening and are simply supported. A linear analysis is used to determine the laminate stress, strain, and mode shape when short-wavelength buckling occurs. The equations for the laminate compressive stress at short-wavelength buckling are dominated by matrix contributions.

Shuart, M. J.↗

NASA NTRS1988

An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations are solved in a time-accurate manner in using the method of pseudocompressibility. Using this method, subiterations in pseudotime are required to satisfy the continuity equation at each time step. An upwind differencing scheme based on flux-difference splitting is used to compute the convective terms. The upwind differencing is biased based on the sign of the local eigenvalue of the Jacobian matrix. Third-order or fifth-order spatial accuracy is maintained throughout the interior grid points. The equations are solved using an implicit line-relaxation scheme. This solution scheme is stable and is capable of running at large time steps in pseudotime, leading to fast convergence for each physical time step. A variety of computed results are presented to validate the present scheme. Results for the flow over an oscillating plate are compared with the exact analytic solution, good agreement is seen. Excellent comparison is obtained between the computed solution and the analytical results for inviscid channel flow with an oscillating back pressure. Flow solutions over a circular cylinder with vortex shedding are also presented. Finally, the flow past an airfoil at -90 deg angle-of-attack is also computed.

Rogers, Stuart E.↗

NASA NTRS1949

Recommendations for numerical solution of reinforced-panel and fuselage-ring problems

Procedures are recommended for solving the equations of equilibrium of reinforced panels and isolated fuselage rings as represented by the external loads and the operations table established according to Southwell's method. From the solution of these equations the stress distribution can be easily determined. The method of systematic relaxations, the matrix-calculus method, and several other methods applicable in special cases are discussed. Definite recommendations are made for obtaining the solution of reinforced-panel problems which are generally designated as shear lag problems. The procedures recommended are demonstrated in the analysis of a number of panels. In the case of fuselage rings it is not possible to make definite recommendations for the solution of the equilibrium equations for all rings and loadings. However, suggestions based on the latest experience are made and demonstrated on several rings.

Hoff, N J↗

NASA NTRS1992

ICAN/DAMP-integrated composite analyzer with damping analysis capabilities: User's manual

This manual describes the use of the computer code ICAN/DAMP (Integrated Composite Analyzer with Damping Analysis Capabilities) for the prediction of damping in polymer-matrix composites. The code is written in FORTRAN 77 and is a version of the ICAN (Integrated Composite ANalyzer) computer program. The code incorporates a new module for synthesizing the material damping from micromechanics to laminate level. Explicit micromechanics equations based on hysteretic damping are programmed relating the on-axis damping capacities to the fiber and matrix properties and fiber volume ratio. The damping capacities of unidirectional composites subjected to off-axis loading are synthesized from on-axis damping values. The hygrothermal effect on the damping performance of unidirectional composites caused by temperature and moisture variation is modeled along with the damping contributions from interfacial friction between broken fibers and matrix. The temperature rise is continuously vibrating composite plies and composite laminates is also estimated. The ICAN/DAMP user's manual provides descriptions of the damping analysis module's functions, structure, input requirements, output interpretation, and execution requirements. It only addresses the changes required to conduct the damping analysis and is used in conjunction with the 'Second Generation Integrated Composite Analyzer (ICAN) Computer Code' user's manual (NASA TP-3290).

Saravanos, Dimitrious A.↗

DOE OSTI2025

Wormholes, branes and finite matrices in sine dilaton gravity

We compute the double trumpet in sine dilaton gravity via WdW quantization. The wormhole size is discretized. The wormhole amplitude matches the spectral correlation of a finite-cut matrix integral, where matrices have large but finite dimensions. This strongly suggests an identification of the sine dilaton gravity theory with the q-deformed JT gravity matrix integral. At the very least, it captures all universal content of that matrix model. The disk decomposes into the physical (gauge invariant) solutions of the WdW equation, which are trumpets with discrete sizes. This decomposition modifies the usual no-boundary wavefunction to a normalizable one in sine dilaton gravity.

2D Gravity↗

NASA NTRS1979

Analysis of a unidirectional composite containing broken fibers and matrix damage

An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a materials-modeling approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e. epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength depends on crack length while for large splitting it is crack length independent.

George, J. G.↗