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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1989

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

NASA NTRS1990

Time-dependent wave-packet method for the complete determination of S-matrix elements for reactive molecular collisions in three dimensions

An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.

Judson, Richard S.↗

NASA NTRS1981

Vibration-rotation intensities of SiO

Dipole moment matrix elements have been computed for a large number of transitions of astrophysical interest for the more abundant isotopes of SiO. The wave functions utilized were obtained from a direct solution of the Schroedinger equation with an accurate RKR potential. The dipole moment function, in the form of a Pade approximant, was chosen to reproduce the experimental measurements near equilibrium, to have the proper united and separated atom limits, and to have the correct long-range asymptotic dependence on internuclear separation. Because of the large number of transitions involved, and to facilitate applications, the squares of the dipole moment matrix elements were fitted by a least-squares procedure to polynomials in v and J. In addition, Einstein A coefficients are given for observed maser transitions and for selected vibrational bands. These latter are compared with previous calculations, and it is concluded that for the higher Delta v transitions, the present results represent a significant improvement.

Tipping, R. H.↗

NASA NTRS2024

MSFC-TVC | TB-03: Derivation of Thrust Vector Control (TVC) Actuator-Force / Gimbal-Torque Transformation Matrix

Thrust vector control (TVC) systems for rocket engine propulsion traditionally use a simple linear relationship to convert between actuator forces and torques about the engine gimbal’s center-of-rotation (COR). As shown in Equation (1), the torque about the gimbal COR is proportional to the applied actuator force and the TVC moment arm (MA)—the perpendicular distance between the TVC actuator’s line-of-action (LOA) and the engine gimbal’s COR. While this fundamental relationship remains valid and accurate in a two-dimensional (2D), one-degree-of-freedom (1-DOF) context—particularly in its non-linear formulation as described in the ER63 Technical Bulletin TB-02 (Derivation of Thrust Vector Control (TVC) Engine-Gimbal / Actuator Moment-Arm Geometry)—it becomes limited when extended to three-dimensional (3D), two-degree-of-freedom (2-DOF) analyses. In 3D, out-of-plane angular displacements can arise, causing the gimbaled engine plane-of-motion to become non-coplanar with the respective actuator plane-of-motion. Such misalignments occur due to the rod-end (RE) and/or tail-stock (TE) mounting geometry of individual TVC actuators. These geometric complexities lead to inaccuracies in calculating engine torques and corresponding actuator forces if only the traditional TVC MA relationship is employed. For higher gimbal angular displacements (e.g., greater than 5 degrees), these inaccuracies become more pronounced, necessitating a more robust mathematical framework. This document presents a comprehensive 3D (2-DOF) derivation of a transformation matrix that accurately converts between gimbaled engine torques and TVC actuator forces. By incorporating the necessary geometric and rotational considerations, this new approach corrects the limitations of the traditional TVC MA method. Subsequent sections compare the newly formulated approach to the traditional method, demonstrating its enhanced accuracy and reliability for a broad range of gimbaled engine conditions.

Thrust Vector Control↗

NASA NTRS1972

Corrections to star catalogues from satellite observations.

Theoretical investigation of the possibility of obtaining systematic corrections to star catalogs from observations of artificial satellites. A model is established to represent the system of equations formed when such corrections are determined simultaneously with corrections to the geocentric position of the observer. An analytical expression for the covariance matrix is set up for the two-dimensional (planar) case. The resulting correlations and weights are discussed. Since this method presents an independent way of obtaining the systematic errors in star catalogs, apart from the orientation to the equatorial system, it is related to the measurement of all quantities which depend on star positions. This institutes a completely new principle of obtaining systematic corrections to star catalogs. It is based not on rotation, but on comparing the angle subtended at the satellite by two stars to the angle subtended at the satellite by two points on the surface of the earth.

Williams, C. A.↗

NASA NTRS1974

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗

NASA NTRS1974

Prediction of compressor stall for distorted and undistorted flow by use of a multistage compressor simulation on the digital computer

Daniele, C. J.↗

NASA NTRS1975

Prediction of compressor stall for distorted and undistorted flow by use of a multistage compressor simulation on the digital computer

A simulation technique is presented for the prediction of compressor stall for axial-flow compressors for clean and distorted inlet flow. The simulation is implemented on the digital computer and uses stage stacking and lumped-volume gas dynamics. The resulting nonlinear differential equations are linearized about a steady-state operating point, and a Routh-Hurwitz stability test is performed on the linear system matrix. Parallel compressor theory is utilized to extend the technique to the distorted inlet flow problem. The method is applied to the eight-stage J85-13 compressor. Analytical stall prediction for the undistorted stall line shows good agreement with experimental results. The predicted stall line for distorted inlet flow is in agreement with experimental results only for large distortion extents and/or low distortion levels. Results for low distortion extents and high distortion level do not agree with experimental results.

Daniele, C. J.↗

NASA NTRS1977

Dynamic characteristics of rotor blades with pendulum absorbers

The point transmission matrix for a vertical plane pendulum on a rotating blade undergoing combined flapwise bending, and chordwise bending and torsion is derived. The equilibrium equation of the pendulum is linearized for small oscillations about the steady state. A FORTRAN program was written for the case of a vertical plane pendulum attached to a uniform blade with flapwise bending degree of freedom for cantilever boundary conditions. The frequency has a singular value right at the uncoupled pendulum natural frequency and thus introduces two frequencies corresponding to the nearest natural frequency of the blade without pendulum. In both of these modes it was observed that the pendulum deflection is large. One frequency can be thought of as a coupled pendulum frequency and the other as a coupled bending and pendulum frequency.

Murthy, V. R.↗

NASA NTRS1977

An efficient numerical technique for calculating thermal spreading resistance

An efficient numerical technique for solving the equations resulting from finite difference analyses of fields governed by Poisson's equation is presented. The method is direct (noniterative)and the computer work required varies with the square of the order of the coefficient matrix. The computational work required varies with the cube of this order for standard inversion techniques, e.g., Gaussian elimination, Jordan, Doolittle, etc.

Gale, E. H., Jr.↗

NASA NTRS1979

A high order accurate numerical solution algorithm for turbulent boundary layer flow

A fourth-order accurate numerical solution algorithm is derived using finite element interpolation theory for the non-linear parabolic equations governing turbulent boundary layer flow including a two-equation turbulence closure model. The results of carefully controlled numerical experiments firmly quantize for the first time performance differences between finite element and finite difference solution methodology for this type of equation. The developed algorithm takes advantage of the apparent semi-analytical formulational procedure, in establishment of a single, retarded-evaluation Jacobian matrix iterative solution algorithm. Numerical results document performance of solution economy features in terms of computer requirements and solution accuracy. The developed algorithm should find wide application in aerodynamics analysis.

Soliman, M. O.↗

NASA NTRS1981

A new global operator for two-particle delta functions

A new type of global operator to be used in evaluating matrix elements of two-particle delta functions is introduced. It is based, like the Trivedi one-particle operator, on the Poisson equation and is easier to apply than the method of Hiller, Sucher and Feinberg. After a test in the helium isoelectronic sequence, the new method is applied successfully to the interesting problem of hyperfine structure in muonic helium.

Drachman, R. J.↗

NASA NTRS1983

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.

Zebker, H. A.↗

NASA NTRS1989

Free-vibration characteristics and correlation of a space station split-blanket solar array

Two methods for studying the free-vibration characteristics of a large split-blanket solar array in a zero-g cantilevered configuration are presented. The zero-g configuration corrresponds to an on-orbit configuration of the Space Station solar array. The first method applies the equations of continuum mechanics to determine the natural frequencies of the array; the second uses the finite element method program, MSC/NASTRAN. The stiffness matrix from the NASTRAN solution was found to be erroneously grounded. The results from the two methods are compared. It is concluded that the grounding does not seriously compromise the solution to the elastic modes of the solar array. However, the correct rigid body modes need to be included to obtain the correct dynamic model.

Carney, Kelly S.↗

NASA NTRS1989

Free-vibration characteristics and correlation of a Space Station split-blanket solar array

Two methods for studying the free-vibration characteristics of a large split-blanket solar array in a zero-g cantilevered configuration are presented. The zero-g configuration corresponds to an on-orbit configuration of the Space Station solar array. The first method applies the equations of continuum mechanics to determine the natural frequencies of the array; the second uses the finite element method program, MSC/NASTRAN. The stiffness matrix from the NASTRAN solution was found to be erroneously grounded. The results from the two methods are compared. It is concluded that the grounding does not seriously compromise the solution to the elastic modes of the solar array. However, the correct rigid body modes need to be icluded to obtain the correct dynamic model.

Carney, Kelly S.↗

NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗

NASA NTRS1997

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗

NASA NTRS2005

System Study: Technology Assessment and Prioritizing

The objective of this NASA funded project is to assess and prioritize advanced technologies required to achieve the goals for an "Intelligent Propulsion System" through collaboration among GEAE, NASA, and Georgia Tech. Key GEAE deliverables are parametric response surface equations (RSE's) relating technology features to system benefits (sfc, weight, fuel burn, design range, acoustics, emission, etc...) and listings of Technology Impact Matrix (TIM) with benefits, debits, and approximate readiness status. TIM has been completed for GEAE and NASA proposed technologies. The combined GEAE and NASA TIM input requirement is shown in Table.1. In the course of building the RSE's and TIM, significant parametric technology modeling and RSE accuracy improvements were accomplished. GEAE has also done preliminary ranking of the technologies using Georgia Tech/GEAE USA developed technology evaluation tools. System level impact was performed by combining beneficial technologies with minimum conflict among various system figures of merits to assess their overall benefits to the system. The shortfalls and issues with modeling the proposed technologies are identified, and recommendations for future work are also proposed.

Source record↗