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,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Computational Simplification of Robot-Arm Dynamics

Report presents two general methods for reducing mathematical complexities of state equations of robot-arm dynamics. Robot arms contain both rotary and linear joints. Both methods start with homogenous coordinates and Lagrangian formulation of mechanics briefly summarized in report. First method uses matrix-analysis techniques; second, vector-analysis techniques. Vector-analysis method includes new differential-vector representation of centripetal and Coriolis forces.

Bejczy, A. K.↗

NASA NTRS1985

Photoionization of sodium atoms and electron scattering from ionized sodium

The polarized-orbital method of Temkin (1957) is applied using polarized orbitals determined from Sternheimer's equation to compute the photoionization cross sections of Na atoms from threshold to about 60 eV. The approximations involved in the analysis are explained in detail; the explicit forms of the integrals and matrix expressions are given in appendices; and the results are presented in tables and graphs. Good agreement is found with the results of Chang and Kelly (1975), and the possibility that small amounts of molecular vapor in Na-photoionization experiments are responsible for the discrepancies between calculated and measured cross sections is considered.

Dasgupta, A.↗

NASA NTRS1990

Fundamental aspects of and failure modes in high-temperature composites

Fundamental aspects of and attendant failure mechanisms for high temperature composites are summarized. These include: (1) in-situ matrix behavior; (2) load transfer; (3) limits on matrix ductility to survive a given number of cyclic loadings; (4) fundamental parameters which govern thermal stresses; (5) vibration stresses; and (6) impact resistance. The resulting guidelines are presented in terms of simple equations which are suitable for the preliminary assessment of the merits of a particular high temperature composite in a specific application.

Chamis, Christos C.↗

NASA NTRS2020

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined.

Franciolini, Matteo↗

NASA NTRS2022

Multiple Scattering Effects of Trees on L-band Microwave Using a Hybrid Method

The multiple scattering of electromagnetic waves in forests is studied using a two-step hybrid method. First, the T-matrix of a single tree is calculated based on its far-field computed with the full-wave simulations of the FEKO software. The T-matrix captures the multiple scattering caused by the tree structure. Second, the the interactions among different trees are considered using the T-matrices of the individual trees and the Foldy-Lax equations. The result of the two-step hybrid method is validated with FEKO by solving scattering from two trees directly. The multiple scattering effects are illustrated by the field solutions.

Yueh, Simon↗

NASA NTRS1986

Analysis of Delamination Growth from Matrix Cracks in Laminates Subjected to Bending Loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg. plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg. plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg. ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg. ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg. ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, G. B.↗

NASA NTRS1988

Analysis of delamination growth from matrix cracks in laminates subjected to bending loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, Gretchen Bostaph↗

NASA NTRS1984

Preliminary results in implementing a model of the world economy on the CYBER 205: A case of large sparse nonsymmetric linear equations

A brief description of the Model of the World Economy implemented at the Institute for Economic Analysis is presented, together with our experience in converting the software to vector code. For each time period, the model is reduced to a linear system of over 2000 variables. The matrix of coefficients has a bordered block diagonal structure, and we show how some of the matrix operations can be carried out on all diagonal blocks at once.

Szyld, D. B.↗

NASA NTRS1974

On the theory of large amplitude Alfven waves

Large amplitude Alfvenic disturbances of arbitrary spatial shape and polarization are described by MHD equations, without resort to the usual assumption of planarity. However, because of their nonplanar nature, the direction of propagation of these disturbances cannot, in general, be determined by looking for minima in a variance matrix constructed from observed field fluctuations. When such minima exist, one is observing that subset of interplanetary Alfven waves that is essentially planar.

Goldstein, M. L.↗

NASA NTRS1979

Mechanical property characterization of intraply hybrid composites

An investigation of the mechanical properties of intraply hybrids made from graphite fiber/epoxy matrix hybridized with secondary S-glass or Kevlar 49 fiber composites is presented. The specimen stress-strain behavior was determined, showing that mechanical properties of intraply hybrid composites can be measured with available methods such as the ten-degree off-axis test for intralaminar shear, and conventional tests for tensile, flexure, and Izod impact properties. The results also showed that combinations of high modulus graphite/S-glass/epoxy matrix composites exist which yield intraply hybrid laminates with the best 'balanced' properties, and that the translation efficiency of mechanical properties from the constituent composites to intraply hybrids may be assessed with a simple equation.

Chamis, C. C.↗

NASA NTRS1991

A natural partitioning scheme for parallel simulation of multibody systems

A parallel partitioning scheme based on physical-coordinate variables is presented to systematically eliminate system constraint forces and yield the equations of motion of multibody dynamics systems in terms of their independent coordinates. Key features of the present scheme include an explicit determination of the independent coordinates, a parallel construction of the null space matrix of the constraint Jacobian matrix, an easy incorporation of the previously developed two-stage staggered solution procedure, and Schur complement based parallel preconditioned conjugate gradient numerical algorithm.

Chiou, J. C.↗

NASA NTRS1992

Fast methods to numerically integrate the Reynolds equation for gas fluid films

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid journal gas bearings by using Liebmann's iterative solution (LIS) for elliptic partial differential equations. The results are compared with known design criteria from experimental data. The developed methods show good accuracy and very short computer running time in comparison with methods based on an inverting of a matrix. The computer codes need a small amount of memory and can be run on either personal computers or on mainframe systems.

Dimofte, Florin↗

NASA NTRS1993

A natural partitioning scheme for parallel simulation of multibody systems

A parallel partitioning scheme based on physical-co-ordinate variables is presented to systematically eliminate system constraint forces and yield the equations of motion of multibody dynamics systems in terms of their independent coordinates. Key features of the present scheme include an explicit determination of the independent coordinates, a parallel construction of the null space matrix of the constraint Jacobian matrix, an easy incorporation of the previously developed two-stage staggered solution procedure and a Schur complement based parallel preconditioned conjugate gradient numerical algorithm.

Chiou, J. C.↗

NASA NTRS1995

Modeling of DC spacecraft power systems

Future spacecraft power systems must be capable of supplying power to various loads. This delivery of power may necessitate the use of high-voltage, high-power dc distribution systems to transmit power from the source to the loads. Using state-of-the-art power conditioning electronics such as dc-dc converters, complex series and parallel configurations may be required at the interface between the source and the distribution system and between the loads and the distribution system. This research will use state-variables to model and simulate a dc spacecraft power system. Each component of the dc power system will be treated as a multiport network, and a state model will be written with the port voltages as the inputs. The state model of a component will be solved independently from the other components using its state transition matrix. A state-space averaging method is developed first in general for any dc-dc switching converter, and then demonstrated in detail for the particular case of the boost power stage. General equations for both steady-state (dc) and dynamic effects (ac) are obtained, from which important transfer functions are derived and applied to a special case of the boost power stage.

Berry, F. C.↗

NASA NTRS2001

Inherent and Apparent Scattering Properties of Coated or Uncoated Spheres Embedded in an Absorbing Host Medium

The conventional Lorenz-Mie formalism is extended to the scattering process associated with a coated sphere embedded in an absorbing medium. It is shown that apparent and inherent scattering cross sections of a scattering particle, which are identical in the case of transparent host medium, are different if the host medium is absorptive. Here the inherent single-scattering properties are derived from the near-field information whereas the corresponding apparent counterparts are derived from the far-field asymptotic form of the scattered wave with scaling of host absorption that is assumed to be in an exponential form. The formality extinction and scattering efficiencies defined in the same manner as in the conventional sense can be unbounded. For a nonabsorptive particle embedded in an absorbing medium, the effect of host absorption on the phase matrix elements associated with polarization is significant. This effect, however, is largely reduced for strongly absorptive particles such as soot. For soot particles coated with water, the impurity can substantially reduce the single-scattering albedo of the particle if the size parameter is small. For water-coating soot and hollow ice spheres, it is shown that the phase matrix elements -P(sub 12)/P(sub 11) and P(sub 33)/P(sub 11) are unique if the shell is thin, as compared with the case for thick shell. Furthermore, the radiative transfer equation regarding a multidisperse particle system in an absorbing medium is discussed. It is illustrated that the conventional computation algorithms can be applied to solve the multiple scattering process if the scaled apparent single-scattering properties are applied.

Yang, P.↗

NASA NTRS2011

Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides

This work addresses the problem of excitation of a parallel plate waveguide by an array of rectangular waveguides that arises in applications such as the continuous transverse stub (CTS) antenna and dual-polarized parabolic cylindrical reflector antennas excited by a scanning line source. In order to design the junction region between the parallel plate waveguide and the linear array of rectangular waveguides, waveguide sizes have to be chosen so that the input match is adequate for the range of scan angles for both polarizations. Electromagnetic wave scattered by the junction of a parallel plate waveguide by an array of rectangular waveguides is analyzed by formulating coupled integral equations for the aperture electric field at the junction. The integral equations are solved by the method of moments. In order to make the computational process efficient and accurate, the method of weighted averaging was used to evaluate rapidly oscillating integrals encountered in the moment matrix. In addition, the real axis spectral integral is evaluated in a deformed contour for speed and accuracy. The MoM results for a large finite array have been validated by comparing its reflection coefficients with corresponding results for an infinite array generated by the commercial finite element code, HFSS. Once the aperture electric field is determined by MoM, the input reflection coefficients at each waveguide port, and coupling for each polarization over the range of useful scan angles, are easily obtained. Results for the input impedance and coupling characteristics for both the vertical and horizontal polarizations are presented over a range of scan angles. It is shown that the scan range is limited to about 35 for both polarizations and therefore the optimum waveguide is a square of size equal to about 0.62 free space wavelength.

Rengarajan, Sembiam↗

NASA NTRS1972

The general linear inverse problem - Implication of surface waves and free oscillations for earth structure.

The discrete general linear inverse problem reduces to a set of m equations in n unknowns. There is generally no unique solution, but we can find k linear combinations of parameters for which restraints are determined. The parameter combinations are given by the eigenvectors of the coefficient matrix. The number k is determined by the ratio of the standard deviations of the observations to the allowable standard deviations in the resulting solution. Various linear combinations of the eigenvectors can be used to determine parameter resolution and information distribution among the observations. Thus we can determine where information comes from among the observations and exactly how it constraints the set of possible models. The application of such analyses to surface-wave and free-oscillation observations indicates that (1) phase, group, and amplitude observations for any particular mode provide basically the same type of information about the model; (2) observations of overtones can enhance the resolution considerably; and (3) the degree of resolution has generally been overestimated for many model determinations made from surface waves.

Wiggins, R. A.↗

NASA NTRS1984

Kinematic control of robot with degenerate wrist

Kinematic resolved rate equations allow an operator with visual feedback to dynamically control a robot hand. When the robot wrist is degenerate, the computed joint angle rates exceed operational limits, and unwanted hand movements can result. The generalized matrix inverse solution can also produce unwanted responses. A method is introduced to control the robot hand in the region of the degenerate robot wrist. The method uses a coordinated movement of the first and third joints of the robot wrist to locate the second wrist joint axis for movement of the robot hand in the commanded direction. The method does not entail infinite joint angle rates.

Barker, L. K.↗