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 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

Chiou, Jin-Chern↗

NASA NTRS1990

Micro-mechanical analysis of damage growth and fracture in discontinuous fiber reinforced metal matrix composites

The near-crack-tip stresses in any planar coupon of arbitrary geometry subjected to mode 1 loading may be equated to those in an infinite center-cracked panel subjected to the appropriate equivalent remote biaxial stresses (ERBS). Since this process can be done for all such mode 1 coupons, attention may be focused on the behavior of the equivalent infinite cracked panel. To calculate the ERBS, the constant term in the series expansion of the crack-tip stress must be retained. It is proposed that the ERBS may be used quantitatively to explain different fracture phenomena such as crack branching.

Goree, James G.↗

NASA NTRS1981

Two-boundary grid generation for the solution of the three dimensional compressible Navier-Stokes equations

A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.

Smith, R. E.↗

NASA NTRS1977

An extension to the Chahine method of inverting the radiative transfer equation

An extension of the Chahine relaxation method (1970) for inverting the radiative transfer equation is presented. This method is superior to the original method in that it takes into account in a realistic manner the shape of the kernel function, and its extension to nonlinear systems is much more straightforward. A comparison of the new method with a matrix method due to Twomey (1965), in a problem involving inference of vertical distribution of ozone from spectroscopic measurements in the near ultraviolet, indicates that in this situation this method is stable with errors in the input data up to 4%, whereas the matrix method breaks down at these levels. The problem of non-uniqueness of the solution, which is a property of the system of equations rather than of any particular algorithm for solving them, remains, although it takes on slightly different forms for the two algorithms.

Twomey, S.↗

NASA NTRS1996

T-Matrix Computations of Light Scattering by Nonspherical Particles: A Review

We review the current status of Waterman's T-matrix approach which is one of the most powerful and widely used tools for accurately computing light scattering by nonspherical particles, both single and composite, based on directly solving Maxwell's equations. Specifically, we discuss the analytical method for computing orientationally-averaged light-scattering characteristics for ensembles of nonspherical particles, the methods for overcoming the numerical instability in calculating the T matrix for single nonspherical particles with large size parameters and/or extreme geometries, and the superposition approach for computing light scattering by composite/aggregated particles. Our discussion is accompanies by multiple numerical examples demonstrating the capabilities of the T-matrix approach and showing effects of nonsphericity of simple convex particles (spheroids) on light scattering.

Mischenko, Michael I.↗

NASA NTRS1972

Space shuttle GN and C equation document no. 20: Fine alignment gyro torquing angles

The procedure for aligning the inertial reference unit of the space shuttle by celestial references is discussed. The gyro torquing angles required to go from the present stable member to the desired stable member are determined. The inputs are the transformation matrix from the reference to the desired stable member coordinates and the unit vectors to the two target stars in both the reference and present stable member coordinates. The equations involved in establishing the reference are developed.

St.amand, J.↗

NASA NTRS2019

Comprehensive Thematic T-Matrix Reference Database: a 2017-2019 Update

Since its inception in the mid-1960s, the T-matrix method has been one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It has been widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This compilation is the ninth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2017. It also includes a few earlier publications that have so far been overlooked.

Michael I Mishchenko↗

NASA NTRS1973

The inference of atmospheric ozone using satellite nadir measurements in the 1042/cm band

A description and detailed analysis of a technique for inferring atmospheric ozone information from satellite nadir measurements in the 1042 cm band are presented. A method is formulated for computing the emission from the lower boundary under the satellite which circumvents the difficult analytical problems caused by the presence of atmospheric clouds and the watervapor continuum absorption. The inversion equations are expanded in terms of the eigenvectors and eigenvalues of a least-squares-solution matrix, and an analysis is performed to determine the information content of the radiance measurements. Under favorable conditions there are only two pieces of independent information available from the measurements: (1) the total ozone and (2) the altitude of the primary maximum in the ozone profile.

Russell, J. M., III↗

NASA NTRS1974

Crystallization of dense neutron matter

The equation of state for cold neutron matter at high density is studied in the t-matrix formulation, and it is shown that energetically it is convenient to have neutrons in a crystalline configuration rather than in a liquid state for values of the density exceeding 1600 Tg/cu cm. The study of the mechanical properties indicates that the system is stable against shearing stresses. A solid core in the deep interior of heavy neutron stars appears to offer the most plausible explanation of speed-ups observed in the Vela pulsar.

Canuto, V.↗

NASA NTRS1977

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

An accurate mathematical model for sound propagation in axisymmetric aircraft engine ducts with compressible mean flow is reported. The model is based on the usual perturbation of the basic fluid mechanics equations for small motions. Mean flow parameters are derived in the absence of fluctuating quantities and are then substituted into the equations for the acoustic quantities which were linearized by eliminating higher order terms. Mean swirl is assumed to be zero from the restriction of axisymmetry. A linear rectangular serendipity element is formulated from these equations using a Galerkin procedure and assembled in a special purpose computer program in which the matrix map for a rectangular mesh was specifically coded. Representations of the fluctuating quantities, mean quantities and coordinate transformations are isoparametric. The global matrix is solved by foreward and back substitution following an L-U decomposition with pivoting restricted internally to the blocks. Results from the model were compared with results from several alternative analyses and yielded satisfactory agreement.

Abrahamson, A. L.↗

NASA NTRS1978

Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics (SOUSSA) for complex aircraft configurations

The Green's function method and the computer program SOUSSA (Steady Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed. The Green's function method is applied to the fully unsteady potential equation yielding an integro-differential-delay equation. This equation is approximated by a set of differential-delay equations in time using the finite element method. The Laplace transform is used to yield a matrix relating the velocity potential to the normal wash. The matrix of the generalized aerodynamic forces is obtained by premultiplying and postmultiplying the matrices relating generalized forces to the potential and the normal wash by the generalized coordinates. The program SOUSSA is compared with existing numerical results. Results indicate that the program is not only general, flexible, and easy to use, but also accurate and fast.

Morino, L.↗

NASA NTRS1982

An exact zooming method for finite element analyses

An exact zooming technique which employs static condensation and exact structural reanalysis methods was developed. Successive application of static condensation reduces the system to one that is only associated with the degrees of freedom (dof) of the original model. Application of an exact static reanalysis technique permits the displacements at the dof of the original model that are contained in the zoomed portion of the structure to be obtained first. The response external to the zoom, as well as the response of additional dof within various levels of zooming, is computed. With the triangular factor of the stiffness matrix of the original system available, this approach involves only the solution of a system of equations of small order.

Hirai, I.↗

NASA NTRS1984

Predicting Thermal Conductivity

Empirical equation predicts thermal conductivity of composite insulators consisting of cellular, granular or fibrous material embedded in matrix of solid viscoelastic material. Application in designing custom insulators for particular environments.

Penn, B.↗

NASA NTRS1989

Atomic data for opacity calculations. X - Oscillator strengths and photoionisation cross sections for O III

Energy levels and electric dipole radiative transitions were determined for O III. The wavefunctions for the bound and continuum states were derived by solving coupled integrodifferential equations in the close-coupling approximation with the aid of two methods, the R-matrix and linear algebraic methods. Configuration interaction wavefunctions are presented for eight states of the O IV target with configurations of 2s(x)sp(y) (x + y = 3). Oscillator strengths are calculated for the transitions between, and photoionization cross sections from, bound states of O III with configurations of 2s(x)2p(y)nl (for n of not greater than 10 and l of not greater than 3).

Luo, D.↗

NASA NTRS1996

Research in Computational Aeroscience Applications Implemented on Advanced Parallel Computing Systems

Improving the numerical linear algebra routines for use in new Navier-Stokes codes, specifically Tim Barth's unstructured grid code, with spin-offs to TRANAIR is reported. A fast distance calculation routine for Navier-Stokes codes using the new one-equation turbulence models is written. The primary focus of this work was devoted to improving matrix-iterative methods. New algorithms have been developed which activate the full potential of classical Cray-class computers as well as distributed-memory parallel computers.

Wigton, Larry↗

NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗

NASA NTRS1992

Kinematic equations for control of the redundant eight-degree-of-freedom advanced research manipulator 2

The forward position and velocity kinematics for the redundant eight-degree-of-freedom Advanced Research Manipulator 2 (ARM2) are presented. Inverse position and velocity kinematic solutions are also presented. The approach in this paper is to specify two of the unknowns and solve for the remaining six unknowns. Two unknowns can be specified with two restrictions. First, the elbow joint angle and rate cannot be specified because they are known from the end-effector position and velocity. Second, one unknown must be specified from the four-jointed wrist, and the second from joints that translate the wrist, elbow joint excluded. There are eight solutions to the inverse position problem. The inverse velocity solution is unique, assuming the Jacobian matrix is not singular. A discussion of singularities is based on specifying two joint rates and analyzing the reduced Jacobian matrix. When this matrix is singular, the generalized inverse may be used as an alternate solution. Computer simulations were developed to verify the equations. Examples demonstrate agreement between forward and inverse solutions.

Williams, Robert L., II↗

NASA NTRS1975

Prolongation structures of nonlinear evolution equations

A technique is developed for systematically deriving a 'prolongation structure' - a set of interrelated potentials and pseudopotentials - for nonlinear partial differential equations in two independent variables. When this is applied to the Korteweg-de Vries equation, a new infinite set of conserved quantities is obtained. Known solution techniques are shown to result from the discovery of such a structure: related partial differential equations for the potential functions, linear 'inverse scattering' equations for auxiliary functions, Backlund transformations. Generalizations of these techniques will result from the use of irreducible matrix representations of the prolongation structure.

Wahlquist, H. D.↗