Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

Kinematic properties of rotary-wing and fixed-wing aircraft in steady coordinated high-g turns

An analytical approach to the study of flight dynamics of aircraft operating in a high-angle-of-attack flight regime and of helicopters operating in extreme thrust conditions is presented. Steady coordinated high-g turns are used to establish the initial equilibrium flight conditions near stall angles of attack. The kinematic properties of the aircraft in steady coordinated turns are examined: in high-g turns, pitch rate (independent of the angle of attack) is of a much larger magnitude than roll and yaw rate; a substantial roll rate is found to develop in steep turns for all angles of attack; the angle of attack also has a significant effect on the pitch attitude, with decreasing influence as the normal load factor increases. The exact small disturbance equations of motion of the aircraft in general steady turns are also developed for application to both rotary-wing and fixed-wing aircraft in extreme conditions. These equations are in a first-order, vector-matrix format, and are thus compatible with many efficient software packages developed in modern system theory.

Chen, R. T. N.↗

Determining 3-D motion and structure from image sequences

A method of determining three-dimensional motion and structure from two image frames is presented. The method requires eight point correspondences between the two frames, from which motion and structure parameters are determined by solving a set of eight linear equations and a singular value decomposition of a 3x3 matrix. It is shown that the solution thus obtained is unique.

Huang, T. S.↗

Kalman filter modeling

The formulation of appropriate state-space models for Kalman filtering applications is studied. The so-called model is completely specified by four matrix parameters and the initial conditions of the recursive equations. Once these are determined, the die is cast, and the way in which the measurements are weighted is determined foreverafter. Thus, finding a model that fits the physical situation at hand is all important. Also, it is often the most difficult aspect of designing a Kalman filter. Formulation of discrete state models from the spectral density and ARMA random process descriptions is discussed. Finally, it is pointed out that many common processes encountered in applied work (such as band-limited white noise) simply do not lend themselves very well to Kalman filter modeling.

Brown, R. G.↗

Calculation of compressible flow about three-dimensional inlets with auxiliary inlets, slats and vanes by means of a panel method

An efficient and user oriented method was constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: (1) the use of a panel method; (2) a technique of superposition for obtaining solutions at any inlet operating condition; and (3) employment of an advanced matrix iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that yield a more complete comprehension of the flowfield than was possible previously.

Hess, J. L.↗

Calculation of compressible flow about three-dimensional inlets with auxiliary inlets, slats and vanes by means of a panel method

An efficient and user-oriented method has been constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: the use of a panel method, a technique of superposition for obtaining solutions at any inlet operating condition, and employment of an advanced matrix-iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that, taken together, yield a more complete comprehension of the flowfield than had been possible previously. Examples of these features are presented for some complicated configurations, and where possible, comparisons are made between calculation and experiment.

Hess, J. L.↗

Development of improved gravity field estimates using existing altimeter data

The formulation of a Fourier series and geographic representation expressions for orbit error were reported. Addition activities include the study of undulation errors in the Lagrangian formulation, and the development and testing of all the necessary software for the above expression. Additionally, the matrix formulations that led to the formation of observed equations and covariance propagation were developed.

Source record↗

Compliance matrices for cracked bodies

An algorithm is developed to construct the compliance matrix for a cracked solid in the integral-equation formulation of two-dimensional linear-elastic fracture mechanics. The integral equation is reduced to a system of algebraic equations for unknown values of the dislocation-density function at discrete points on the interval from -1 to 1, using the numerical procedure described by Gerasoulis (1982). Sample numerical results are presented, and it is suggested that the algorithm is especially useful in cases where iterative solutions are required; e.g., models of fiber-reinforced concrete, rocks, or ceramics where microcracking, fiber bridging, and other nonlinear effects are treated as nonlinear springs along the crack surfaces (Ballarini et al., 1984).

Ballarini, R.↗

Determination of settings of a tilted head-cutter for generation of hypoid and spiral bevel gears

Kinematics of Gleason mechanisms of hypoid and spiral bevel cutting machines are considered. These mechanisms are designated to install the position and tilt of the head cutter. The tilt of the head cutter with standard blades provides the required pressure angle. The authors have developed the matrix presentation of kinematics of these meachanisms and basic equations for the required settings. An example is presented based on the developed computation procedure.

Litvin, F. L.↗

Krylov subspace methods on supercomputers

A short survey of recent research on Krylov subspace methods with emphasis on implementation on vector and parallel computers is presented. Conjugate gradient methods have proven very useful on traditional scalar computers, and their popularity is likely to increase as three-dimensional models gain importance. A conservative approach to derive effective iterative techniques for supercomputers has been to find efficient parallel/vector implementations of the standard algorithms. The main source of difficulty in the incomplete factorization preconditionings is in the solution of the triangular systems at each step. A few approaches consisting of implementing efficient forward and backward triangular solutions are described in detail. Polynomial preconditioning as an alternative to standard incomplete factorization techniques is also discussed. Another efficient approach is to reorder the equations so as to improve the structure of the matrix to achieve better parallelism or vectorization. An overview of these and other ideas and their effectiveness or potential for different types of architectures is given.

Saad, Youcef↗

Parallel processing techniques for finite element analysis of nonlinear large truss structures

Methods were developed for parallel processing of finite element solutions of large truss structures. The parallel processing techniques were implemented in two stages, i.e., the repeated forming of the nonlinear global stiffness matrix and the solving of the global system of equations. The Sequent Balance 21000 parallel computer was employed to demonstrate the procedures and the speed-up.

Chien, L. S.↗

A new technique for simulating composite material. Task 2: Analytical solutions with Generalized Impedance Boundary Conditions (GIBCs)

The diffraction problem associated with a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GSMF) in conjunction with the dual integral equation approach. The multilayer slab is replaced by a surface obeying a generalized impedance boundary condition (GIBC) to facilitate the computation of the pertinent Wiener Hopf split functions and their zeros. Both E(sub z) and H(sub z) polarizations are considered and a number of scattering patterns are presented, some of which are compared to exact results available for a homogeneous recessed slab.

Ricoy, M. A.↗

Electromagnetic scattering from two-dimensional thick material junctions

The problem of the plane wave diffraction is examined by an arbitrary symmetric two dimensional junction, where Generalized Impedance Boundary Conditions (GIBCs) and Generalized Sheet Transition Conditions (GSTCs) are employed to simulate the slabs. GIBCs and GSTCs are constructed for multilayer planar slabs of arbitrary thickness and the resulting GIBC/GSTC reflection coefficients are compared with exact counterparts to evaluate the GIBCs/GSTCs. The plane wave diffraction by a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GDMF) in conjunction with the dual integral equation approach. Various scattering patterns are computed and validated with exact results where possible. The diffraction by a material discontinuity in a thick dielectric/ferrite slab is considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of unknown constants is obtained, and these constants are evaluated for the recessed slab geometry by comparison with the solution obtained therein. Several other simplified cases are also presented and discussed. An eigenfunction expansion method is introduced to determine the unknown solution constants in the general case. This procedure is applied to the non-unique solution in terms of unknown constants; and scattering patterns are presented for various slab junctions and compared with alternative results where possible.

Ricoy, M. A.↗

Parallel-vector computation for CSI-design code

Computational aspects of Control-Structure Interaction (CSI) DESIGN code is reviewed. Numerical intensive computation portions of CSI-DESIGN code were identified. Improvements in computational speed for the CSI-DESIGN code can be achieved by exploiting parallel and vector capabilities offered by modern computers, such as the Alliant, Convex, Cray-2, and Cray-YMP. Four options to generate the coefficient stiffness matrix and to solve the system of linear, simultaneous equations are currently available in the CSI-DESIGN code. A preprocessor to use RCM (Reverse Cuthill-Mackee) algorithm for bandwidth minimization was also developed for the CSI-DESIGN code. Preliminary results obtained by solving a small-scale, 97 node CSI finite element model (for eigensolution) have indicated that this new CSI-DESIGN code is 5 to 6 times faster (using 1 Alliant processor) than the old version of CSI-DESIGN code. This speed-up was achieved due to the RCM algorithm and the use of a new skyline solver. Efforts are underway to further improve the vector speed for CSI-DESIGN code, to evaluate its performance on a larger scale CSI model (such as phase zero CSI model) to make the code run efficiently on multiprocessor, parallel computer environment, and to make the code portable among different parallel computers available at NASA LaRC, such as Alliant, Convex, and Cray computers.

Nguyen, Duc T.↗

A pressure flux-split technique for computation of inlet flow behavior

A method for calculating the flow field in aircraft engine inlets is presented. The phenomena of inlet unstart and restart are investigated. Solutions of the reduced Navier-Stokes (RNS) equations are obtained with a time consistent direct sparse matrix solver that computes the transient flow field both internal and external to the inlet. Time varying shocks and time varying recirculation regions can be 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.↗

Dynamic analysis of truss-beam system

A simple truss-beam method for determining the dynamic characteristics of the space structures intended to perform various tasks in orbit is presented. Algorithms are provided to determine the flexibility matrix of the truss beam for use in the equation of motion. The natural frequencies obtained through this method are compared with those obtained through the finite element method. An experimental procedure for verifying the theoretical results is described. It is concluded that the truss-beam method is a simple analysis technique that yields reasonably accurate results with a minimum of computational effort, especially for the critical lower bending modes of flexible structures. The method takes less computer time than more conventional methods and can be programed in both BASIC and FORTRAN for use in micro- and mainframe computers.

Abu-Saba, Elias G.↗

Eigenvalue calculation procedure for an Euler/Navier-Stokes solver with application to flows over airfoils

A Lanczos procedure is presently applied to a Navier-Stokes (N-S) solver for eigenvalues and eigenvectors associated with the small-perturbation analysis of the N-S equations' finite-difference representation for airfoil flows; the matrix used is very large, sparse, real, and nonsymmetric. The Lanczos procedure is shown to furnish complete spectral information for the eigenvalues, as required for transient-stability analysis of N-S solvers.

Mahajan, Aparajit J.↗

Computational micromechanics of woven composites

The bounds on the equivalent elastic material properties of a composite are presently addressed by a unified energy approach which is valid for both unidirectional and 2D and 3D woven composites. The unit cell considered is assumed to consist, first, of the actual composite arrangement of the fibers and matrix material, and then, of an equivalent pseudohomogeneous material. Equating the strain energies due to the two arrangements yields an estimate of the upper bound for the material equivalent properties; successive increases in the order of displacement field that is assumed in the composite arrangement will successively produce improved upper bound estimates.

Hopkins, Dale A.↗

Implementation of a multiblock sensitivity analysis method in numerical aerodynamic shape optimization

A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations. Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies. The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from the CFD discretized flow equations, recast in adjoint form, and solved using direct matrix inversion techniques. This investigation is performed to demonstrate an aerodynamic shape optimization technique on a multiblock domain and its applicability to complex geometries. The objectives are accomplished by shape optimizing two aerodynamic configurations. First, the shape optimization of a transonic airfoil is performed to investigate the behavior of the method in highly nonlinear flows and the effect of different grid blocking strategies on the procedure. Secondly, shape optimization of a two-element configuration in subsonic flow is completed. Cases are presented for this configuration to demonstrate the effect of simultaneously reshaping interfering elements. The aerodynamic shape optimization is shown to produce supercritical type airfoils in the transonic flow from an initially symmetric airfoil. Multiblocking effects the path of optimization while providing similar results at the conclusion. Simultaneous reshaping of elements is shown to be more effective than individual element reshaping due to the inclusion of mutual interference effects.

Lacasse, James M.↗