Engineering Papers⌕ Search

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.

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

At least 775 records · Page 43

Scaling Laws Applied to a Modal Formulation of the Aeroservoelastic Equations

A method of scaling is described that easily converts the aeroelastic equations of motion of a full-sized aircraft into ones of a wind-tunnel model. To implement the method, a set of rules is provided for the conversion process involving matrix operations with scale factors. In addition, a technique for analytically incorporating a spring mounting system into the aeroelastic equations is also presented. As an example problem, a finite element model of a full-sized aircraft is introduced from the High Speed Research (HSR) program to exercise the scaling method. With a set of scale factor values, a brief outline is given of a procedure to generate the first-order aeroservoelastic analytical model representing the wind-tunnel model. To verify the scaling process as applied to the example problem, the root-locus patterns from the full-sized vehicle and the wind-tunnel model are compared to see if the root magnitudes scale with the frequency scale factor value. Selected time-history results are given from a numerical simulation of an active-controlled wind-tunnel model to demonstrate the utility of the scaling process.

Pototzky, Anthony S.↗

GENERAL INSTABILITY OF ORTHOGONALLY STIFFENED CYLINDRICAL SHELLS

Earlier research at the National Aeronautics Research Institute (N.L.R.), Amsterdam, which forms the basis of recent work is reviewed. This early work refers to 2 schemes: the orthotropic shell and, in view of buckling modes where the half wave length is of the order of the ring distance, the shell with continuously distributed stringers and discrete rings. Linear theory is considered to be adequate for these structures, where the imperfections are small in comparison to the height of the ring sections. Recent developments account for pressure difference in addition to axial compression, for the correct stiffness matrix of skin panels in the post-buckling stage and for stringer bending due to hoop stresses in the skin, which are of importance as has been shown by the investigation of the post-buckling behaviour. Numerical data for the stiffness matrix of skin panels have been established. Numerical evaluation of the stability equation has not been performed as yet.

SHELL STABILITY↗

Analysis of heat transfer in a porous cooled wall with variable pressure and temperature along the coolant exit boundary

Fluid from a reservior at constant pressure and temperature is forced through a porous wall of uniform thickness. The boundary through which the fluid exits has specified variations in pressure and temperature along it in one direction so that the flow and heat transfer are two-dimensional. The local fluid and matrix temperatures are assumed to be equal and therefore a single energy equation governs the temperature distribution within the wall. The solution is obtained by transforming this energy equation into potential plane coordinates, which results in a separable equation. A technique yielding an integral equation is used to adapt the general solution so that it satisfies the variable-pressure boundary condition. Analytical expressions are given for the normal exit velocity and heat flux along the exit boundary. Illustrative examples are carried out which indicate to what extent the solution is locally one-dimensional.

Siegel, R.↗

Two-dimensional finite-element temperature variance analysis

The finite element method is extended to thermal analysis by forming a variance analysis of temperature results so that the sensitivity of predicted temperatures to uncertainties in input variables is determined. The temperature fields within a finite number of elements are described in terms of the temperatures of vertices and the variational principle is used to minimize the integral equation describing thermal potential energy. A computer calculation yields the desired solution matrix of predicted temperatures and provides information about initial thermal parameters and their associated errors. Sample calculations show that all predicted temperatures are most effected by temperature values along fixed boundaries; more accurate specifications of these temperatures reduce errors in thermal calculations.

Heuser, J. S.↗

A numerical algorithm suggested by problems of transport in periodic media - The matrix case.

Extension of Allen and Wing's (1970) previous work on problems of transport in periodic media to the matrix case. A method in the form of a complete set of equations is presented that may be used without any further analytical work by investigators interested in computing solutions to problems of the type the method is designed to handle. All the formulas have been checked out numerically, and their effectiveness is demonstrated by numerical examples.

Allen, R. C., Jr.↗

An improved general algorithm for arbitrary pole assignment.

A new pole-placement algorithm for multi-input controllable systems is presented. The algorithm can be used to obtain real-valued feedback gains for arbitrarily specified real and complex, distinct or multiple eigenvalues. The multi-input controllable time-invariant system is converted to a single-input system. Pole assignment is achieved by making use of certain properties of the equation TF - AT = D, where A is the open-loop system matrix and F represents a matrix having the desired closed-loop eigenvalues.

Willner, L.↗

Diffusion in a unidirectional filament reinforced metal composite.

A theoretical analysis has been developed to describe composition changes resulting from exposure of a unidirectional filament reinforced metal composite to an elevated temperature environment. The analysis was based on a simple superposition of finite-difference solutions of the diffusion equation, and was applicable to the case for which the filament and matrix metals comprised an isomorphous binary alloy system. The possibility of chemical reaction between filaments and matrix was excluded. Excellent agreement was obtained between calculated concentration profiles and results of quantitative electron microprobe analyses of nickel filament reinforced copper samples.

Herring, H. W.↗

Aeroelastic stability of coupled flap-lag motion of hingeless helicopter blades at arbitrary advance ratios

Equations for large amplitude coupled flap-lag motion of a hingeless elastic helicopter blade in forward flight are derived. Only a torsionally rigid blade excited by quasi-steady aerodynamic loads is considered. The effects of reversed flow together with some new terms due to radial flow are included. Using Galerkin's method the spatial dependence is eliminated and the equations are linearized about a suitable equilibrium position. The resulting system of homogeneous periodic equations is solved using multivariable Floquet-Liapunov theory, and the transition matrix at the end of the period is evaluated by two separate methods. Computational efficiency of the two numerical methods is compared. Results illustrating the effects of forward flight and various important blade parameters on the stability boundaries are presented.

Friedmann, P.↗

Flap-lag dynamics of hingeless helicopter blades at moderate and high advance ratios

Equations for large amplitude coupled flaplag motion of a hingeless elastic helicopter blade in forward flight are derived. Only a torsionally rigid blade exicted by quasi-steady aerodynamic loads is considered. The effects of reversed flow together with some new terms due to forward flight are included. Using Galerkin's method the spatial dependence is eliminated and the equations are linearized about a suitable equilibrium position. The resulting system of equations is solved using multivariable Floquet-Liapunov theory, and the transition matrix at the end of the period is evaluated by two separate methods. Results illustrating the effects of forward flight and various important blade parameters on the stability boundaries are presented.

Friedman, P.↗

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗

Triangular system solutions on an optical systolic processor

It is noted that the simplified system of equations presented by Casasent and Ghosh (1983) can also be solved optically. One realization of this is described, with attention given to data flow and speed for the case of a lower triangular system of equations. The solution of an upper triangular system of equations can be handled in analogous fashion. The problem requires a matrix-vector processor. Only one carrier frequency input to the acoustooptic cell and only one output detector are used. To multiply a matrix by a vector in this system, the time-multiplexed vector elements are fed into the acoustooptic cell.

Ghosh, A.↗

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

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

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

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

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↗

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↗