SEARCH · Engineering Papers
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.
GAMBIT computer program provides multigroup flux and current averaged neutron cross sections for input to transport programs. The fast neutron compilation includes P/0/ and P/1/ elastic scattering transfer matrices, inelastic scattering matrices, n,2n cross section matrices, absorption and fission cross sections for all isotopes.
Our method of matrix synthesis of optical components and instruments is applied to the derivation of Jones's matrices appropriate for Fourier interferometers (spectrometers and spectropolarimeters). These matrices are obtained for both the source beam and the detector beam. In the course of synthesis, Jones's matrices of the various reflectors (plane mirrors; retroreflectors: roofed mirror, trihedral and prism cube corner, cat's eye) used by these interferometers are also obtained.
A technique is described which is utilized for modeling human operator dynamics. The technique is based upon a model-reference system configuration and a class of Liapunov functions formulated for this purpose which possess variable characteristics. It is shown how such a class of Liapunov functions and their time derivatives is formulated. The crux of the formulation lies in three variable positive definite matrices used for the construction of the class of Liapunov functions. The form and order of these matrices depend upon the form and order of the error differential equation of the human operator model-reference system. These matrices can be modified to include nonlinear functions of the human operator model.
A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.
The computer printout data generated during the Payload/Orbiter Contamination Control Requirement Study are presented. The computer listings of the input surface data matrices, the viewfactor data matrices, and the geometric relationship data matrices for the three orbiter/spacelab configurations analyzed in this study are given. These configurations have been broken up into the geometrical surfaces and nodes necessary to define the principal critical surfaces whether they are contaminant sources, experimental surfaces, or operational surfaces. A numbering scheme was established based upon nodal numbers that relates the various spacelab surfaces to a specific surface material or function. This numbering system was developed for the spacelab configurations such that future extension to a surface mapping capability could be developed as required.
The finite-element method has been applied to solve a combined radiative-conductive heat transfer problem for a large space telescope similar to those used in orbiting satellites. The derivation of the underlying matrices and associated solution algorithm for a 2-D triangular element is presented in detail. The resulting expressions for this triangular element typify such an analysis, which yields constitutive matrices when the heat equation is cast in the matrix form. The relevant matrices include those pertaining to thermal conductance, internal heat generation, radiative exchanges, and all possible external thermal loadings. Emphasis is placed on the treatment of non-linear radiative interchange between surfaces in an enclosure having mixed diffuse-specular surface characteristics. Essential differences in governing equations describing these distinctive surface characteristics are identified. Concluding remarks are drawn from an example simulating a Cassegrainian space telescope.
Radiative transfer of partially polarized radiation in an anisotropically scattering, inhomogeneous atmosphere containing an arbitrary polydispersion of particles is described using Jones's (1941) amplitude vectors and matrices. This approach exploits the close analogy between the quantum-mechanical states of spin 1/2 systems and the polarization states of electromagnetic radiation described by Jones's vector. The complete equivalence between the transport equation for Jones's vectors and the classical radiative-transfer equation for Stokes's intensity vectors is demonstrated in two independent ways after deriving the transport equations for the polarization coherency matrices and for the quaternions corresponding to the Jones's vectors. A compact-operator formulation of the theory is provided and used to derive the necessary equations for both a local and a global description of the transport of Jones's vectors. The integro-differential equations for the amplitude reflection and transmission matrices are derived and related to the usual corresponding equations.
The performance of Kalman-type, linear, discrete-time filters in the presence of modeling errors is considered. The discussion is limited to stationary performance, and bounds are obtained for the performance index, the mean-squared error of estimates for suboptimal and optimal (Kalman) filters. The computation of these bounds requires information on only the model matrices and the range of errors for these matrices. Consequently, a design can easily compare the performance of a suboptimal filter with that of the optimal filter, when only the range of errors in the elements of the model matrices is available.
This paper explores the concept of finite dynamic elements involving higher order dynamic correction terms in the associated stiffness and mass matrices. Such matrices are then developed for a rectangular prestressed membrane element. Next, efficient analysis techniques for the eigenproblem solution of the resulting quadratic matrix equations are described in detail. These are followed by suitable numerical examples which indicate that employment of such dynamic elements in conjunction with an efficient quadratic matric solution technique will result in a most significant economy in the free vibration analysis of structures.
A major difficulty in the practical application of linear-quadratic regulator theory is how to choose the weighting matrices in quadratic cost functions. The control system design with optimal weighting matrices was applied to a helicopter in the hover and approach phase. The weighting matrices were calculated to extremize the closed loop total system damping subject to constraints on the determinants. The extremization is really a minimization of the effects of disturbances, and interpreted as a compromise between the generalized system accuracy and the generalized system response speed. The trade-off between the accuracy and the response speed is adjusted by a single parameter, the ratio of determinants. By this approach an objective measure can be obtained for the design of a control system. The measure is to be determined by the system requirements.
The paper describes the split-Cholesky strategy for banded matrices arising from the large systems of equations in certain fluid mechanics problems. The basic idea is that for a banded matrix the computation can be carried out in pieces, with only a small portion of the matrix residing in core. Mesh considerations are discussed by demonstrating the manner in which the assembly of finite element equations proceeds for linear trial functions on a triangular mesh. The FORTRAN code which implements the out-of-core decomposition strategy for banded symmetric positive definite matrices (mass matrices) of a coupled initial value problem is given.
This paper describes experiments performed to measure the dynamic transfer matrices for cavitating (and noncavitating) pumps. These transfer matrices describe the relationship between small linear oscillatory perturbations in the pressures and mass flow rates at inlet and discharge from the hydraulic machine. The matrices were deduced from direct measurements of these fluctuating quantities for different modes of excitation of the machine. Results for a cavitating inducer are presented as functions of frequency and mean operating state. Though some of the trends in the data are consistent with existing theoretical models of inducer dynamics, others are not, indicating a need for further theoretical investigation of the dynamic characteristics of such flows. The results exhibit increasingly complex dynamics with increasing cavitation; it appears that the hydraulic machine deviates from an essentially passive response without cavitation to an increasingly active response as the cavitation number is reduced.
A finite-element thermal analysis procedure for elements with several temperature-dependent thermal parameters is presented. The procedure, based on an application of the Newton-Raphson iteration technique, is formulated by resolving element matrices into component matrices, one component for each thermal parameter. Component conductance matrices are evaluated by assuming constant thermal parameters within an element and are computed once per unit thermal parameter. Significant savings in computer time result from the unit thermal parameter concept. The solution procedure applied to a convectively cooled structure with significantly varying thermal parameters converged in four iterations.
The use of the SUDAN computer program for analyzing structural systems for their natural modes and frequencies of vibration is described. SUDAN is intended for structures which can be represented as an equivalent system of beam, spring, and rigid-body substructures. User-written constraint equations are used to analytically join the mass and stiffness matrices of the substructures to form the mass and stiffness matrices of the complete structure from which all the frequencies and modes of the system are determined. The SUDAN program can treat the case in which both the mass and stiffness matrices of the coupled system may be singular simultaneously. A general description of the FORTRAN IV program is given, the computer hardware and software specifications are indicated, and the input required by the program is described.
The method used to formulate the fluid/structure equations is described in the updates to the NASTRAN Theoretical Manual. The basis for defining the fluid are three dimensional finite elements connected to fluid grid points defining the Eulerian pressure at a point fixed in space. The use of a pressure single degree of freedom at each point rather than three displacements allows a finer mesh of elements with a reasonable matrix order. In the formulation of the fluid/structure system the interior fluid degrees of freedom are transformed and removed from the solution matrices. The eigenvalues of the combination are extracted from small, fully dense, symmetric mass and stiffness matrices, efficiently processed with the Givens method. The solution matrices are defined by the free surface displacements and the reduced structure coordinates.
Higher order correction terms for the stiffness and inertia matrices associated with a triangular plane stress-strain finite dynamic element are developed in detail. Numerical results are presented which indicate that the adoption of these matrices along with a suitable quadratic matrix eigenproblem solver effects a significant economy in the free vibration solution of structure when compared with the analysis based on the usual finite element procedure. Finally, a FORTRAN IV computer program listing of the various relevant matrices is given.
This report describes the third phase of work, the objective of which was to overcome the excessive brittleness of the previously developed UH-1 helicopter tail rotor drive shaft design which demonstrated a shaft train weight savings of 53.1% over the current 2024-T3 aluminum shaft train. A materials impact program demonstrated exceptionally noteworthy performance of two woven constructions containing E-glass and PRD 49-III (designation later changed to KEVLAR 49) fibers in an epoxy resin matrix. Thermoplastic matrices and PRD 49-III fiber provided impact resistance at low weight which was superior to composites having the same fiber in a thermoset resin matrix. A design, fabrication, and test program showed that shaft impact resistance could be improved over the previously developed graphite composite design at a cost in shaft train rate savings. The shaft train weight savings of the most impact tolerant construction was 4.0% over the current aluminum shaft train. Alternating plies of graphite and glass appear to provide substantially greater tube impact durability than that provided by hybridization of the two fibers into one tape wound to a ply design equivalent in strength and stiffness to that of the alternating ply design. Recommendations were made to continue research work to exploit the potential for more impact-durable structures through the use of KEVLAR 49 fiber, woven structures, thermoplastic matrices and THORNEL 50-S/KEVLAR 49 blends with thermoset matrices.
Petrographic analyses of CM matrices characterized four phyllosilicates in Murray and Murchison meteorites and Fe- and Mg-serpentines in Nogoya. All phyllosilicates and bulk matrices show enrichment of K relative to Na when compared with bulk meteorites; the loss of Na and some Cl, and the addition of H2O, CO2, and water-soluble organics during alteration indicates a partially open system. Synthesis of soluble organic materials may have occurred in CM matrices before aqueous alteration of the precursive phases. Nogoya was 95% altered and has a bulk C content of 5.2%, higher than any meteorite; also, it has the lowest measured C-13/C-12 ratio of any carbonaceous chondrite except for Karoonda.