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 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Interactive application of quadratic expansion of chi-square statistic to nonlinear curve fitting

This report contains a detailed theoretical description of an all-purpose, interactive curve-fitting routine that is based on P. R. Bevington's description of the quadratic expansion of the Chi-Square statistic. The method is implemented in the associated interactive, graphics-based computer program. Taylor's expansion of Chi-Square is first introduced, and justifications for retaining only the first term are presented. From the expansion, a set of n simultaneous linear equations is derived, then solved by matrix algebra. A brief description of the code is presented along with a limited number of changes that are required to customize the program of a particular task. To evaluate the performance of the method and the goodness of nonlinear curve fitting, two typical engineering problems are examined and the graphical and tabular output of each is discussed. A complete listing of the entire package is included as an appendix.

Badavi, F. F.↗

NASA NTRS1986

Vector solution for the mean electromagnetic fields in a layer of random particles

The mean electromagnetic fields are found in a layer of randomly oriented particles lying over a half space. A matrix-dyadic formulation of Maxwell's equations is employed in conjunction with the Foldy-Lax approximation to obtain equations for the mean fields. A two variable perturbation procedure, valid in the limit of small fractional volume, is then used to derive uncoupled equations for the slowly varying amplitudes of the mean wave. These equations are solved to obtain explicit expressions for the mean electromagnetic fields in the slab region in the general case of arbitrarily oriented particles and arbitrary polarization of the incident radiation. Numerical examples are given for the application to remote sensing of vegetation.

Lang, R. H.↗

NASA NTRS1989

Control of a slow moving space crane as an adaptive structure

Assuming that the space crane is an adaptive structure with length-adjustable bars and taking as controls the length-adjustments of these bars, the computation of the incremental controls corresponding to the motion of a payload along its minimum-energy trajectory is given in terms of the inverse-transpose of matrix B of the joint equilibrium equations Bs = p, where s lists the bar forces and p lists the nodal loads. The compensation of the controls for elastic deformations and support movements are shown. It is also shown that the computations may be done automatically and in real time by an attached processor once the characteristics of the crane's maneuver are keyed in.

Utku, S.↗

NASA NTRS1989

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗

NASA NTRS1990

A flux-split solution procedure for unsteady inlet flows

The unstart and restart of an axisymmetric inlet is investigated using a flux-split procedure applied to the Euler and Reduced Navier Stokes (RNS) equations. A time consistent direct sparse matrix solver is applied to compute the transient flow field both internal and external to the inlet. Time varying oblique and normal shocks are captured. The code is quite general and is applicable for subsonic, transonic and supersonic free streams. The current analysis is concerned with supersonic flight conditions.

Pordal, H. S.↗

NASA NTRS1990

A diagonal implicit scheme for computing flows with finite-rate chemistry

A new algorithm for solving steady, finite-rate chemistry, flow problems is presented. The new scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The source Jacobian matrix is replaced by a diagonal matrix which is tailored to account for the fastest reactions in the chemical system. A point-implicit procedure is discussed and then the algorithm is included into the LU-SGS scheme. Solutions are presented for hypervelocity reentry and Hydrogen-Oxygen combustion. For the LU-SGS scheme a CFL number in excess of 10,000 has been achieved.

Eberhardt, Scott↗

NASA NTRS1990

Transient behavior of supersonic flow through inlets

A solution technique to compute inlet flow behavior is presented. The phenomena of inlet unstart and restart are investigated using a flux-split procedure applied to the Euler and Reduced Navier Stokes (RNS) equations. A time consistent direct sparse matrix solver in conjunction with a domain decomposition strategy is applied to compute the transient flow behavior both internal and external to the inlet. Time varying shocks and time varying recirculation regions are efficiently analyzed. The code is quite general and is suitable for the computation of flow for a wide variety of geometries and over a wide range of Mach and Reynolds numbers.

Pordal, H. S.↗

NASA NTRS1990

A manual for PARTI runtime primitives

Primitives are presented that are designed to help users efficiently program irregular problems (e.g., unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial differential equations solvers) on distributed memory machines. These primitives are also designed for use in compilers for distributed memory multiprocessors. Communications patterns are captured at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗

NASA NTRS1991

Coupled Riccati equations for complex plane constraint

A new Linear Quadratic Gaussian design method is presented which provides prescribed imaginary axis pole placement for optimal control and estimation systems. This procedure contributes another degree of design freedom to flexible spacecraft control. Current design methods which interject modal damping into the system tend to have little affect on modal frequencies, i.e., they predictably shift open plant poles horizontally in the complex plane to form the closed loop controller or estimator pole constellation, but make little provision for vertical (imaginary axis) pole shifts. Imaginary axis shifts which reduce the closed loop model frequencies (the bandwidths) are desirable since they reduce the sensitivity of the system to noise disturbances. The new method drives the closed loop modal frequencies to predictable (specified) levels, frequencies as low as zero rad/sec (real axis pole placement) can be achieved. The design procedure works through rotational and translational destabilizations of the plant, and a coupling of two independently solved algebraic Riccati equations through a structured state weighting matrix. Two new concepts, gain transference and Q equivalency, are introduced and their use shown.

Strong, Kristin M.↗

NASA NTRS1991

A manual for PARTI runtime primitives, revision 1

Das, Raja↗

NASA NTRS1990

Unsteady airloading panel method for prop-fans

An unsteady, linear, compressible, three dimensional lifting surface panel method is presented for use with propellers and prop-fans. Both gust and blade vibration problems can be treated. Derivation of the theory is based on the acceleration potential method and can be considered a generalization of the analogous wing theory to include rotation and multiple blades. As with the wing methods, the integral equation is discretized and solved by matrix inversion. Sample calculations are use to explore the difference between 2D and 3D results and the significance of unsteadiness in angular inflow calculations for propellers. It is shown at the once-per-revolution frequency of propellers, that quasi-steady methods should be adequate. At the high frequencies of interest in noise analysis, 2D strip methods give adequate results away from blade tips. For flutter, frequencies lie in an intermediate range where a full 3D unsteady technique seems to be required.

Hanson, D. B.↗

NASA NTRS1990

Time-dependent treatment of scattering. II - Novel integral equation approach to quantum wave packets

The novel wave-packet propagation scheme presented is based on the time-dependent form of the Lippman-Schwinger integral equation and does not require extensive matrix inversions, thereby facilitating application to systems in which some degrees of freedom express the potential in a basis expansion. The matrix to be inverted is a function of the kinetic energy operator, and is accordingly diagonal in a Bessel function basis set. Transition amplitudes for various orbital angular momentum quantum numbers are obtainable via either Fourier transform of the amplitude density from the time to the energy domain, or the direct analysis of the scattered wave packet.

Sharafeddin, Omar A.↗

NASA NTRS1991

Variational finite element-tensor formulation for the large deflection random vibration of composite plates

A numerical integration routine is derived from a set of unified single step integration algorithms using a weighted satisfaction of the equilibrium equation governing the large deflection random response of laminated composite plates. The equilibrium equation is derived using a constant matrix large deflection finite-element formulation. In-plane inertia terms are considered in the formulation, however, rotary inertia terms are assumed negligible. Probability density, spectral density and autocorrelation functions of the maximum displacement and strain responses are presented for three acoustic excitation levels. Classical thin plate boundary conditions and pseudo white noise excitation are used in this investigation.

Robinson, Jay H.↗

NASA NTRS1946

Numerical procedures for the calculation of the stresses in monocoques III : calculation of the bending moments in fuselage frames

This report deals with the calculation of the bending moments in and the distortions of fuselage rings upon which known concentrated and distributed loads are acting. In the procedure suggested, the ring is divided into a number of beams each having a constant radius of curvature. The forces and moments caused in the end sections of the beams by individual unit displacements of the end sections are listed in a table designated as the operations table in conformity with Southwell's nomenclature. The operations table and the external loads are equivalent to a set of linear equations. For their solution the following three procedures are presented: 1) Southwell's method of systematic relaxations. This is a step-by-step approximation procedure guided by the physical interpretation of the changes in the values of the unknown. 2) The growing unit procedure in which the individual beams are combined successively into beams of increasing length until finally the entire ring becomes a single beam. In each step of the procedure a set of not more than three simultaneous linear equations is solved. 3) Solution of the entire set of simultaneous equations by the methods of the matrix calculus. In order to demonstrate the manner in which the calculations may be carried out, the following numerical examples are worked out: 1) Curved beam with both its end sections rigidly fixed. The load is a concentrated force. 2) Egg-shape ring with symmetric concentrated loads. 3) Circular ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 4) Same with V-braces incorporated in the ring. 5) Egg-shape ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 6) Same with V-braces incorporated in the ring. The results of these calculations are checked, whenever possible, by calculations carried out according to known methods of analysis. The agreement is found to be good. The amount of work necessary for the solution of ring problems by the methods described in the present report is practically independent of the degree of redundancy of the structure. For this reason the methods are recommended for use particularly in problems of rings having one or more internal bracing elements.

LOADS - FUSELAGES↗

NASA NTRS1993

Preconditioning and the limit to the incompressible flow equations

The use of preconditioning methods to accelerate the convergence to a steady state for both the incompressible and compressible fluid dynamic equations are considered. The relation between them for both the continuous problem and the finite difference approximation is also considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Hence, the steady state of the preconditioned system is the same as the steady state of the original system. For finite difference methods the preconditioning can change and improve the steady state solutions. An application to flow around an airfoil is presented.

Turkel, E.↗

NASA NTRS1993

Pseudo-compressibility methods for the incompressible flow equations

Preconditioning methods to accelerate convergence to a steady state for the incompressible fluid dynamics equations are considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Thus the steady state of the preconditioned system is the same as the steady state of the original system. The method is compared to other types of pseudo-compressibility. For finite difference methods preconditioning can change and improve the steady state solutions. An application to viscous flow around a cascade with a non-periodic mesh is presented.

Turkel, Eli↗

NASA NTRS1994

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗

NASA NTRS1996

Robust Stabilization of Uncertain Systems Based on Energy Dissipation Concepts

Robust stability conditions obtained through generalization of the notion of energy dissipation in physical systems are discussed in this report. Linear time-invariant (LTI) systems which dissipate energy corresponding to quadratic power functions are characterized in the time-domain and the frequency-domain, in terms of linear matrix inequalities (LMls) and algebraic Riccati equations (ARE's). A novel characterization of strictly dissipative LTI systems is introduced in this report. Sufficient conditions in terms of dissipativity and strict dissipativity are presented for (1) stability of the feedback interconnection of dissipative LTI systems, (2) stability of dissipative LTI systems with memoryless feedback nonlinearities, and (3) quadratic stability of uncertain linear systems. It is demonstrated that the framework of dissipative LTI systems investigated in this report unifies and extends small gain, passivity, and sector conditions for stability. Techniques for selecting power functions for characterization of uncertain plants and robust controller synthesis based on these stability results are introduced. A spring-mass-damper example is used to illustrate the application of these methods for robust controller synthesis.

Gupta, Sandeep↗