Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigenvalues”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

FOS: Physical sciences↗
NASA NTRS1966

Matrix methods and automation in structural engineering

The partial differential equations of motion of shell structures subject to arbitrary loads may be extremely difficult or even impossible to solve. Since the digital computer is now an available tool for the structural engineer, current research has been directed toward methods which involve matrix formulation of large systems of algebraic equations and matrix solutions for discrete elements rather than the solution of the partial differential equations. These methods require rapid and accurate computer solutions. Recognizing the accuracy problems inherent in working with large matrices, a comprehensive survey of available computer programs was performed for large matrix inversion and for eigenvalue and eigenvector solutions under a Research Grant from the National Aeronautics and Space Administration and the techniques are developed in this paper for using these programs most efficiently for structural applications. The contents of the paper include: automation of matrix compilation, methods of very large matrix inversion and solution of simultaneous equation, and techniques for finding eigenvalues and vectors. In addition, a finite element stiffness matrix approach developed at the Denver Research Institute for both plates and shells subject to arbitrary dynamic loads is described as it was instrumented with complete automation on the digital computer.

Anita S West↗
NASA NTRS1968

Internal low Reynolds number turbulent heat transfer

The results of a semi-analytical and experimental study of adiabatic tube flow and an analytical and experimental investigation of the thermal entrance region for gas flowing through electrically heated circular tubes are presented. Emphasis is placed on the low Reynolds number turbulent flow regime--defined as fully turbulent flow at bulk Reynolds numbers from 3,000 to about 15,000. Adiabatic air velocity and friction data and localized heat transfer measurements for air and helium, at low heating rates, are presented for this range. The adiabatic data were obtained in a 1.61 inch 1D tube for flow at bulk Reynolds numbers from 3,000 to 15,000. A continuous, Reynolds number-dependent, profile is developed from the data by using a modification of Reichardts' wall and middle law eddy diffusivity expressions. The velocity profile satisfies continuity. It is valid for all Reynolds numbers in excess of 3,000 for which the flow is fully turbulent and the Blasius friction factor expression is valid. The thermal entrance problem for a fully developed velocity profile is solved analytically by the method of Sparrow, Hallman, and Siegel. The solution is based on the profile developed from the velocity study. Tabular values of the eigenvalues and normalized Nusselt numbers for gases are presented for it range of Reynolds numbers from 3,000 to 50,000. The axial variation of Nusselt number is found to be correlated by [Nu/Nu_(∞)] = 1 + 0.8 (1+ 70,000 Re^(-3/2)) ((x/D)^(-1)) to within ± 5 per cent for x/D ≥2. The fully developed value agrees with the Dittus - Boelter correlation, Nu_(∞) = 0.021 Re^(0.8) Pr^(0.4) For the eigenvalues, λ^(2)_(n), and the associated constants, A_(n), correlations of the form λ^(2)_(n) =A_(1,n) Re^(-b_(1,n) + C_(1,n) Re^(-d_(1,n)) A_(n) = -A_(2,n) Re^(-b_(2,n) + C_(2,n) Re^(-d_(2,n)) are obtained. The coefficients and powers are presented in tabular form. Heat transfer data are presented, primarily for helium, for the low Reynolds number turbulent range. A one-quarter inch, resistively heated, vertical, circular tube was used for the study. The data cover an axial range from 1.2 to 96 diameters; wall-to-bulk temperature ratios vary from 1 to 1.4. In the low Reynolds number turbulent regime, these data clearly support the present analytical solution rather than the prediction obtained by applying the eddy diffusivity distribution used by Sparrow, Hallman, and Siegel.

Harold C. Reynolds, Jr.↗
NASA NTRS1971

A contracting-interval program for the Danilewski method

The concept of contracting-interval programs is applied to finding the eigenvalues of a matrix. The development is a three-step process in which (1) a program is developed for the reduction of a matrix to Hessenberg form, (2) a program is developed for the reduction of a Hessenberg matrix to colleague form, and (3) the characteristic polynomial with interval coefficients is readily obtained from the interval of colleague matrices. This interval polynomial is then factored into quadratic factors so that the eigenvalues may be obtained. To develop a contracting-interval program for factoring this polynomial with interval coefficients it is necessary to have an iteration method which converges even in the presence of controlled rounding errors. A theorem is stated giving sufficient conditions for the convergence of Newton's method when both the function and its Jacobian cannot be evaluated exactly but errors can be made proportional to the square of the norm of the difference between the previous two iterates. This theorem is applied to prove the convergence of the generalization of the Newton-Bairstow method that is used to obtain quadratic factors of the characteristic polynomial.

Harris, J. D.↗
NASA NTRS1972

Synthesis of shuttle vehicle damping using substructure test results

An empirical method is developed for predicting the modal damping of a combined parallel-stage shuttle model by means of damping measurements performed on the individual substructures. Correlations are first determined for each component in terms of damping energy as a function of peak kinetic energy and modal amplitude. The results are then used to predict component damping energies corresponding to the respective kinetic energies and amplitudes that occur for the new modes of the combined system. Modal characteristics for the system, other than damping, are obtained by a real eigenvalue solution of dynamic equations developed by Hurty's procedure of substructures. System equations, which include component modal damping, are also solved by a complex eigenvalue approach for comparison with results of the empirical method.

Kana, D. D.↗
NASA NTRS1972

On the calculation of panel flutter boundaries.

Methods are described for the complete automation of flutter boundary calculations when the aerodynamic forces are derived from linear three-dimensional unsteady potential flow theory. The usual process of visual curve fairing in the mass ratio versus structural damping plane is replaced by numerical procedures for ordering the eigenvalues in such a way that the n-th eigenvalue is always associated with the same flutter boundary. The mass ratio versus structural damping curves are interpolated via parametric cubic spline functions to produce the desired plots in the stiffness-parameter/mass-ratio plane. The entire process is accomplished in a single computer run.

Gaspers, P. A., Jr.↗
NASA NTRS1972

An improved general algorithm for arbitrary pole assignment.

A new pole-placement algorithm for multi-input controllable systems is presented. The algorithm can be used to obtain real-valued feedback gains for arbitrarily specified real and complex, distinct or multiple eigenvalues. The multi-input controllable time-invariant system is converted to a single-input system. Pole assignment is achieved by making use of certain properties of the equation TF - AT = D, where A is the open-loop system matrix and F represents a matrix having the desired closed-loop eigenvalues.

Willner, L.↗
NASA NTRS1973

Multidimensional analysis of fast-spectrum material replacement measurements for systematic estimation of cross section uncertainties

A series of central core and core-reflector interface sample replacement experiments for 16 materials performed in the NASA heavy-metal-reflected, fast spectrum critical assembly (NCA) were analyzed in four and 13 groups using the GAM 2 cross-section set. The individual worths obtained by TDSN and DOT multidimensional transport theory calculations showed significant differences from the experimental results. These were attributed to cross-section uncertainties in the GAM 2 cross sections. Simultaneous analysis of the measured and calculated sample worths permitted separation of the worths into capture and scattering components which systematically provided fast spectrum averaged correction factors to the magnitudes of the GAM 2 absorption and scattering cross sections. Several Los Alamos clean critical assemblies containing Oy, Ta, and Mo as well as one of the NCA compositions were reanalyzed using the corrected cross sections. In all cases the eigenvalues were significantly improved and were recomputed to within 1 percent of the experimental eigenvalue. A comparable procedure may be used for ENDF cross sections when these are available.

Klann, P. G.↗
NASA NTRS1972

SYMAT, COVAR: Test Procedures for Matrix Calculations

The FORTRAN subroutine SYMAT and related subroutines are described. In essence SYMAT is an iterative algorithm in which the problem of finding eigenvalues and eigenvectors of a real symmetric matrix is transformed into an equivalent problem of finding eigenvalues and eigenvectors of an infinite sequence of matrices of order two. A DEMO PROGRAM contains a subroutine COVAR which is used to compute the covariance matrix (denoted by A) of a data matrix (denoted by X). Since a covariance matrix is symmetric it can be analyzed by using subroutine SYMAT.

Morris, W. L.↗
NASA NTRS1972

Diffraction of a plane pulse by a three-dimensional corner.

The conical solutions for the incidence of a plane pulse on a three-dimensional corner are presented. The corner is represented by a trihedron with one edge perpendicular to the other two. Both the boundary condition of the first kind and that of the second kind are considered. Outside the characteristic sphere of the vertex of the corner, the solution is represented by the well-known conical solutions in two variables. Inside the characteristic sphere, the problem involves three conical variables. By the separation of variables, the problem is reduced to that of an eigenvalue problem with an irregular boundary which is in turn reduced to a system of homogeneous algebraic equations. The eigenvalues are then determined numerically. By the superposition of the conical solutions for plane pulses, the solution for the incidence of a plane wave is obtained. Numerical examples simulating the incidence of a sonic boom on the corner of a structure are presented.

Ting, L.↗
NASA NTRS1973

Transient analysis using conical shell elements

The use of the NASTRAN conical shell element in static, eigenvalue, and direct transient analyses is demonstrated. The results of a NASTRAN static solution of an externally pressurized ring-stiffened cylinder agree well with a theoretical discontinuity analysis. Good agreement is also obtained between the NASTRAN direct transient response of a uniform cylinder to a dynamic end load and one-dimensional solutions obtained using a method of characteristics stress wave code and a standing wave solution. Finally, a NASTRAN eigenvalue analysis is performed on a hydroballistic model idealized with conical shell elements.

Yang, J. C. S.↗
NASA NTRS1973

Fast modal extraction in NASTRAN via the FEER computer program

A new eigensolution routine, FEER (Fast Eigensolution Extraction Routine), used in conjunction with NASTRAN at Israel Aircraft Industries is described. The FEER program is based on an automatic matrix reduction scheme whereby the lower modes of structures with many degrees of freedom can be accurately extracted from a tridiagonal eigenvalue problem whose size is of the same order of magnitude as the number of required modes. The process is effected without arbitrary lumping of masses at selected node points or selection of nodes to be retained in the analysis set. The results of computational efficiency studies are presented, showing major arithmetic operation counts and actual computer run times of FEER as compared to other methods of eigenvalue extraction, including those available in the NASTRAN READ module. It is concluded that the tridiagonal reduction method used in FEER would serve as a valuable addition to NASTRAN for highly increased efficiency in obtaining structural vibration modes.

Newman, M. B.↗
NASA NTRS1974

Results on the two population feature selection problem using probability of correct classification as a criterion

Variational equations are presented for maximizing the probability of correct classification as a function of a 1xn feature selection matrix B for the two-population problem. For the special case of equal covariance matrices the optimal B is unique up to scalar multiples and rank one sufficient. For equal population means, the best 1xn B is an eigenvector corresponding either to the largest or smallest eigenvalue of sigma sub 2 to the minus 1 power sigma sub 1 where sigma sub 1 and sigma sub 2 are the nxn covariance matrices of the two populations. The transformed probability of correct classification depends only on the eigenvalue. Finally, a procedure is proposed for constructing an optimal or nearly optimal kxn matrix of rank k without solving the k-dimensional variational equation.

Peters, B. C.↗
NASA NTRS1974

Perturbation solutions for the influence of forward flight on helicopter rotor flapping stability

The stability of the flapping motion of a helicopter rotor blade in forward flight is investigated, using a perturbation technique which gives analytic expressions for the eigenvalues, including the influence of the periodic aerodynamic forces in forward flight. The perturbation solutions are based on small advance ratio (the ratio of the helicopter forward speed to the rotor tip speed). The rotor configurations considered are a single, independent blade; a teetering rotor; a gimballed rotor with three, four, and five or more blades; and a rotor with N independent blades. The constant coefficient approximation with the equations and degrees of freedom in the nonrotating frame represents the flap dynamic reasonably well for the lower frequency modes, although it cannot, of course, be completely correct. The transfer function of the rotor flap response to sinusoidal pitch input is examined, as an alternative to the eigenvalues as a representation of the dynamic characteristics of the flap motion.

Johnson, W.↗
NASA NTRS1974

Application of antiresonance theory to helicopters

Antiresonance theory is the principle underlying nonresonant nodes in a structure, and covers both nonresonant nodes occurring naturally and those introduced by devices such as dynamic absorbers and antiresonant isolators. The dynamic antiresonant vibration isolator (DAVI) and the nodale module are examples of the applications of transfer antiresonances. It is shown that antiresonances are eigenvalues, and that they can be determined by matrix iteration. Applications of antiresonance theory to helicopter engineering problems, using the antiresonant eigenvalue equation are suggested.

Bartlett, F. D., Jr.↗
NASA NTRS1973

Theoretical study of H2/+/ spectroscopic properties. II, III

Description of the theoretical spectroscopic properties of the 2p pi/sub u/ and 3d sigma/sub g/ excited states of the H2/+/ hydrogen molecular ion. Numerical integration of the Schrodinger equation is used to determine vibration-rotation eigenvalues. Dunham power series expansions are used to determine the equilibrium separation, potential coefficients, and spectroscopic constants. The eigenvalues are used to determine delta-G, Bv, Dv, and Hv.

Beckel, C. L.↗
NASA NTRS1974

Analytical and experimental investigation of a 1/8-scale dynamic model of the shuttle orbiter. Volume 1: Summary report

A 1/8-scale structural dynamics model of the space shuttle orbiter was analyzed using the NASA Structural Analysis System (NASTRAN). Comparison of the calculated eigenvalues with preliminary test data for the unrestrained condition indicate that the analytical model was consistently stiffer, being about 20% higher in the first mode. The eigenvectors show reasonably good agreement with test data. A series of analytical and experimental investigations undertaken to resolve the discrepancy are described. Modifications in the NASTRAN model based upon these investigations resulted in close agreement for both eigenvalues and eigenvectors.

Mason, P. W.↗