Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eigenvector”

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 289 records · Page 16

Eigenfunction methods in magnetospheric radial-diffusion theory

Complete sets of orthonormal basis functions constructed according to a generalization of the quantum-mechanical WKB approximation can be used to generate a nearly-diagonal matrix representation of the radial-transport operator for ring-current ions in the presence of radial diffusion and charge exchange. The resulting eigenfunctions (constructed by weighting the basis functions in proportion to the respective components of the eigenvectors of the matrix representation) and eigenvalues provide a spatial and temporal description of the evolving phase-space density during and following a magnetospheric disturbance (e.g., a magnetic storm). A linear superposition of the basis functions can also be used to eliminate any discrepancy between the steady-state solution of the transport equation and the appropriate WKB approximation of this steady-state solution.

Schulz, Michael↗

Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogeneous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the source terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

Yee, H. C.↗

Improvement of modal scaling factors using mass additive technique

A general investigation into the improvement of modal scaling factors of an experimental modal model using additive technique is discussed. Data base required by the proposed method consists of an experimental modal model (a set of complex eigenvalues and eigenvectors) of the original structure and a corresponding set of complex eigenvalues of the mass-added structure. Three analytical methods,i.e., first order and second order perturbation methods, and local eigenvalue modification technique, are proposed to predict the improved modal scaling factors. Difficulties encountered in scaling closely spaced modes are discussed. Methods to compute the necessary rotational modal vectors at the mass additive points are also proposed to increase the accuracy of the analytical prediction.

Zhang, Qiang↗

Vibration analysis of quasi-symmetric structures

A computational technique is developed to reduce the size of the dynamical model of a large complex quasi-symmetric structure prior to free vibration analysis. The structure is discretized using a two-field mixed FEM model; the response vectors are decomposed into symmetric and antisymmetric components; the eigenvectors of the quasi-symmetric structure are approximated by a linear combination of a few preselected sets of symmetric and antisymmetric modes generated by the reduced-size model; global approximation vectors are selected and generated; and the resulting eigenvalue problem is significantly reduced via a classical Bubnov-Galerkin procedure. Numerical results for a shallow laminated anisotropic panel of quadrilateral planform are presented in extensive tables and graphs, demonstrating the effectiveness of the method.

Noor, Ahmed K.↗

Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media

The transfer of monochromatic radiation in a scattering, absorbing, and emitting plane-parallel medium with a specified bidirectional reflectivity at the lower boundary is considered. The equations and boundary conditions are summarized. The numerical implementation of the theory is discussed with attention given to the reliable and efficient computation of eigenvalues and eigenvectors. Ways of avoiding fatal overflows and ill-conditioning in the matrix inversion needed to determine the integration constants are also presented.

Stamnes, Knut↗

Sequential design of linear quadratic state regulators via the optimal root-locus techniques

The use of well-known root-locus techniques for sequentially finding the weighting matrices and the linear quadratic state regulators of multivariable control systems in the frequency domain is considered. This sequential design method permits the retention of some stable open-loop poles and the associated eigenvectors in the closed-loop system; it also allows some optimal closed-loop poles to be placed in a specific region of the complex plane. In addition, it provides a design procedure for determining the weighting matrices and linear quadratic state regulators for the optimal control of multivariable systems in the frequency domain.

Shieh, L. S.↗

Attitude determination and parameter estimation using vector observations

Procedures for attitude determination based on Wahba's loss function are generalized to include the estimation of parameters other than the attitude, such as sensor biases. Optimization with respect to the attitude requires either the singular value decomposition of a 3x3 matrix or finding the maximum eigenvalue and corresponding eigenvector of a 4x4 symmetric matrix, but does not require an a priori estimate of the attitude. Optimization with respect to the other parameters employs an iterative approach, which does require an a priori estimate of these parameters. Conventional state estimation methods require a priori estimates of both the parameters and the attitude, while the algorithms presented in this paper always compute the exact optimal attitude for given values of the parameters. The proposed method is shown to give the correct solution of an example problem. An expression for the covariance of the attitude and parameter estimates is derived.

Markley, F. Landis↗

Experimental verification of a large flexible manipulator

A large experimental lightweight manipulator would be useful for material handling, for welding, or for ultrasonic inspection of a large structure, such as an airframe. The flexible parallel link mechanism is designed for high rigidity without increasing weight. This constrained system is analyzed by singular value decomposition of the constraint Jacobian matrix. A verification of the modeling using the assumed mode method is presented. Eigenvalues and eigenvectors of the linearized model are compared to the measured system natural frequencies and their associated mode shapes. The modeling results for large motions are compared to the time response data from the experiments. The hydraulic actuator is verified.

Lee, Jac Won↗

Analytical approach to peel stresses in bonded composite stiffened panels

A closed-form solution was obtained for the stresses and displacements of two bonded beams. A system of two fourth-order and two second-order differential equations with the associated boundary equations was determined using a variational work approach. A FORTRAN computer program was devised to solve for the eigenvalues and eigenvectors of this system and to calculate the coefficients from the boundary conditions. The results were then compared with NASTRAN finite-element solutions and shown to agree closely.

Barkey, Derek A.↗

The optimal polarizations for achieving maximum contrast in radar images

There is considerable interest in determining the optimal polarizations that maximize contrast between two scattering classes in polarimetric radar images. A systematic approach is presented for obtaining the optimal polarimetric matched filter, i.e., that filter which produces maximum contrast between two scattering classes. The maximization procedure involves solving an eigenvalue problem where the eigenvector corresponding to the maximum contrast ratio is an optimal polarimetric matched filter. To exhibit the physical significance of this filter, it is transformed into its associated transmitting and receiving polarization states, written in terms of horizontal and vertical vector components. For the special case where the transmitting polarization is fixed, the receiving polarization which maximizes the contrast ratio is also obtained. Polarimetric filtering is then applies to synthetic aperture radar images obtained from the Jet Propulsion Laboratory. It is shown, both numerically and through the use of radar imagery, that maximum image contrast can be realized when data is processed with the optimal polarimeter matched filter.

Swartz, A. A.↗

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented. These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models. The expressions are formulated in the context of either a two or three temperature model. Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v. The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set.

Peter A. Gnoffo↗

Robust eigenstructure assignment by a projection method: Application using multiple optimization criteria

New ideas which lead to feedback control laws for large flexible structures which are insensitive to model uncertainty are presented. A pole placement method is presented which leads to near-unitary closed loop eigenvectors, and a new method is introduced to design the control while simultaneously considering three competing measures of optimality. Robustness versus Integral algorithms are applicable to at least moderately high-dimensioned systems. In the present discussion, controls for two coupled flexible bodies are considered. A 6x24 gain matrix is designed to control a 12 modes system using 6 actuators. Researchers also developed control laws for the R2P2 simulator at Martin Marietta; in this case 3 actuators are used to control a 12th order system. Simulation studies indicate that researchers indeed achieved robust designs without significant difficulties associated with spillover into the uncontrolled modes. Here, several key ideas and numerical results are given. In the references, details of the formulation, discussions of salient features, and connection to the availiable literature are given.

Junkins, J. L.↗

The V3, V4 and V6 bands of formaldehyde: A spectral catalog from 900 cm(-1) to 1580 cm(-1)

The results of a complete high resolution study of the three vibration-rotation bands v sub 3, v sub 4, and V sub 6 using both TDLs and FT-IR spectroscopy are presented. The reults are given in terms of a table of over 8000 predicted transition frequencies and strengths. A plot of the predicted and calculated spectra is shown. Over 3000 transitions were assigned and used in the simultaneous analysis of the three bands. The simultaneous fit permitted a rigorous study of Coriolis and other type iterations among bands yielding improved molecular constants. Line intensities of 28 transitions measured by a TDL and 20 transitions from FTS data were used, along with the eigenvectors from the frequency fitting, in a least squares analysis to evaluate the band strengths.

Nadler, Shachar↗

Analysis of modified SMI method for adaptive array weight control

An adaptive array is applied to the problem of receiving a desired signal in the presence of weak interference signals which need to be suppressed. A modification, suggested by Gupta, of the sample matrix inversion (SMI) algorithm controls the array weights. In the modified SMI algorithm, interference suppression is increased by subtracting a fraction F of the noise power from the diagonal elements of the estimated covariance matrix. Given the true covariance matrix and the desired signal direction, the modified algorithm is shown to maximize a well-defined, intuitive output power ratio criterion. Expressions are derived for the expected value and variance of the array weights and output powers as a function of the fraction F and the number of snapshots used in the covariance matrix estimate. These expressions are compared with computer simulation and good agreement is found. A trade-off is found to exist between the desired level of interference suppression and the number of snapshots required in order to achieve that level with some certainty. The removal of noise eigenvectors from the covariance matrix inverse is also discussed with respect to this application. Finally, the type and severity of errors which occur in the covariance matrix estimate are characterized through simulation.

Dilsavor, R. L.↗

Iterative methods for design sensitivity analysis

A numerical method is presented for design sensitivity analysis, using an iterative-method reanalysis of the structure generated by a small perturbation in the design variable; a forward-difference scheme is then employed to obtain the approximate sensitivity. Algorithms are developed for displacement and stress sensitivity, as well as for eignevalues and eigenvector sensitivity, and the iterative schemes are modified so that the coefficient matrices are constant and therefore decomposed only once.

Belegundu, A. D.↗

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the Goddard Earth Model (GEM)-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment. These calibrations (utilizing each of the major data subsets within the solution) yield very stable calibration factors which vary by approximately 10 percent over the range of tests employed. Measurements of gravity anomalies obtained from altimetry were also used directly as observations to show that GEM-T1 is calibrated. The mathematical representation of the covariance error in the presence of unmodeled systematic error effects in the data is analyzed and an optimum weighting technique is developed for these conditions. This technique yields an internal self-calibration of the error model, a process which GEM-T1 is shown to approximate.

Lerch, F. J.↗

Solving large sparse eigenvalue problems on supercomputers

An important problem in scientific computing consists in finding a few eigenvalues and corresponding eigenvectors of a very large and sparse matrix. The most popular methods to solve these problems are based on projection techniques on appropriate subspaces. The main attraction of these methods is that they only require the use of the matrix in the form of matrix by vector multiplications. The implementations on supercomputers of two such methods for symmetric matrices, namely Lanczos' method and Davidson's method are compared. Since one of the most important operations in these two methods is the multiplication of vectors by the sparse matrix, methods of performing this operation efficiently are discussed. The advantages and the disadvantages of each method are compared and implementation aspects are discussed. Numerical experiments on a one processor CRAY 2 and CRAY X-MP are reported. Possible parallel implementations are also discussed.

Philippe, Bernard↗

Numerical solution of large nonsymmetric eigenvalue problems

Several methods are discribed for combinations of Krylov subspace techniques, deflation procedures and preconditionings, for computing a small number of eigenvalues and eigenvectors or Schur vectors of large sparse matrices. The most effective techniques for solving realistic problems from applications are those methods based on some form of preconditioning and one of several Krylov subspace techniques, such as Arnoldi's method or Lanczos procedure. Two forms of preconditioning are considered: shift-and-invert and polynomial acceleration. The latter presents some advantages for parallel/vector processing but may be ineffective if eigenvalues inside the spectrum are sought. Some algorithmic details are provided that improve the reliability and effectiveness of these techniques.

Saad, Youcef↗