Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced eigenvalue problem”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Reanalysis and design in structural dynamics

A unified efficient formulation for the static and dynamic reanalysis of locally modified structures is presented. The reanalysis problem is to find the structural responses when some of the element properties are adjusted as a result of design modifications or when additional structures are appended to the original system. The reanalysis problem is formulated as a problem of much lower order than the original system. This is achieved by utilizing the linearity property of the structure using the pseudoload concept together with the solution of the original system. The modifications to the structure are treated as displacement-dependent pseudoloads of the system. By expressing the modified system response as linear combinations of the response of the original system and a term depending on the pseudoload, a reduced set of response equations can be obtained. In the static and sinusoidal steady state analysis, this leads to a set of linear algebraic equations. For a free vibration analysis, this results in a reduced eigenvalue problem. Using these reduced equations, a design problem can be solved; that is, the magnitude of a specified local modification for given dynamic response characteristics can be calculated. Several examples are presented to illustrate the general formulations.

Wang, B. P.↗
DOE OSTI2020

Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals

The normal modes of a trapped ion crystal are derived using an approach based on the Hermitian properties of the system’s dynamical matrix. This method is equivalent to the standard Bogoliubov method, but for classical systems it is arguably simpler and more general in that canonical coordinates are not necessary. The theory is developed for stable, unstable, and neutrally-stable systems. The method is then applied to ion crystals in a Penning trap. Reduced eigenvalue problems for the case of large applied magnetic field are developed, for which the spectrum breaks into E x B drift modes, axial modes, and cyclotron modes. Thermal fluctuation levels in these modes are analyzed and shown to be consistent with the Bohr-van-Leeuwen theorem, provided that neutrally-stable modes associated with crystal rotations are included in the analysis. An expression for the rotational inertia of the crystal is derived, and a magnetic contribution to this inertia, which dominates in large magnetic fields, is described. An unusual limit is discovered for the special case of spherically-symmetric confinement, in which the rotational inertia does not exist and changes in angular momentum leave the rotation frequency unaffected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Evaluation and reduction of errors induced by the Guyan transformation

The Guyan reduction refers to a method used to reduce the number of degrees of freedom in a structural model for dynamic analysis. Experience shows that, if the method is properly employed, then this reduction method does in fact provide a "reasonably" accurate approximation of the dynamic characteristics of the unreduced model. A method that provides an accurate estimate of the error in each mode of the reduced eigenvalue problem and, in addition, gives first order correction terms that greatly improve the accuracy of the eigenvectors and eigenvalues is presented. The method is demonstrated by standard COSMIC NASTRAN DAMP alters to rigid format 3, real normal mode analysis. Comparative computer time of modal extraction versus error analysis is discussed for the VAX-11/780 version of COSMIC NASTRAN.

Fox, G. L.↗
DOE OSTI2022

Solving generalized eigenvalue problems by ordinary differential equations on a quantum computer

Many eigenvalue problems arising in practice are often of the generalized form A x = λ B x . One particularly important case is symmetric, namely A , B are Hermitian and B is positive definite. The standard algorithm for solving this class of eigenvalue problems is to reduce them to Hermitian eigenvalue problems. For a quantum computer, quantum phase estimation is a useful technique to solve Hermitian eigenvalue problems. In this work, we propose a new quantum algorithm for symmetric generalized eigenvalue problems using ordinary differential equations. The algorithm has lower complexity than the standard one based on quantum phase estimation. Moreover, it works for a wider case than symmetric: B is invertible, B - 1 A is diagonalizable and all the eigenvalues are real.

Science & Technology - Other Topics↗
DOE OSTI2019

Reduced-order modeling of parameterized multi-group diffusion k-eigenvalue problems

In this paper, a group-wise Reduced-Order Model, based on a Proper Orthogonal Decomposition (POD) technique, is presented for parameterized multi-group diffusion k-eigenvalue problems. The group-wise approach is shown to be more robust than model-order reduction based on monolithic multi-group POD modes. In addition, fewer POD modes are required. The method of snapshots is employed to generate the offline training set. Two Naïve and a Greedy snapshot-generation strategies are compared. The input parameter space is surveyed using Latin Hypercube Sampling techniques. The effectiveness of the group-wise model-order reduction method is demonstrated on two reactor-physics benchmarks, a two-group 3D Pressurized Water Reactor problem with 26 uncertain input parameters and a seven-group UO 2 -MOX fuel 3D mini-core with 287 uncertain input parameters. In both examples the Reduced-Order Models built using the snapshots generated by the Greedy algorithm proved to be slightly better. Model-order reduction yields results that are within 1 pcm of the full-order model, with speed-up factors of about 1000–3000, depending upon the test case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
NASA NTRS1972

A NASTRAN correlation study for vibrations of a cross-stiffened ship's deck

To evaluate the effectiveness of NASTRAN for predicting the vibration modes of panels with bending-membrane coupling, a cross-stiffened ship deck was analyzed. In correlation with experimental data, one NASTRAN finite element representation gave results slightly more accurate than a previous analytical solution. Computational time was excessively long, due to the Guyan method of reducing the eigenvalue problem. It is recommended that a more efficient method of matrix reduction be implemented for the lumped mass formulation.

Thornton, E. A.↗
NASA NTRS1975

Ellipsometer measurements of epitaxial GaAs layers: A concept

Report discusses analysis involving Maxwell's equations formed in a 6-by-6 matrix. By applying boundary conditions at proper points in sample, equation for propagation of light through anisotropic medium is reduced to eigenvalue problem resulting in 4-by-4 matrix.

Desmet, D. J.↗
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 NTRS1987

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.↗
NASA NTRS1988

A comparative study of theoretical methods on Goertler instability

Goertler vortices arise in boundary layers along concave surfaces due to centrifugal effects and these vortices in combination with Tollmien-Schlichting (TS) waves and crossflow vortices may play an important role in triggering early transition. There have been two distinct theoretical approaches to solve the problem. In the classical, normal-mode approach (NMA) the linear development of Goertler vortices is reduced to an eigenvalue problem. In the other method, initiated by Hall, the governing partial differential equations for the perturbations which are parabolic in the streamwise direction, are solved as an initial value problem. This method predicts multiple neutral curves depending on how and where the mean flow is perturbed, unlike the NMA where a unique but different neutral curve is predicted by each investigator. The present paper attempts to compare the two techniques and shows that mutually compatible results can be obtained when physically consistent assumptions are made.

Kalburgi, Vijay↗
NASA NTRSDate not supplied

Towards Formalization of Advanced Linear Algebra with Applications to Dynamical Systems using PVS

Linear Algebra is essential for numerous aerospace problems of interest. Formal reasoning about hybrid systems that contain variables modeled by differential equations rely on concepts from Linear Algebra such as eigenvalues, matrix decompositions, and matrix valued functions. For example, the long-term dynamics of a system of differential equations depend on the stability/instability of its equilibrium points, which often reduces to an eigenvalue problem. This talk will embark on a quest to formalize theorems and results about eigenvalues and eigenvectors using PVS. We shall start our journey with 2 x 2 complex matrices, where we will apply our PVS code to a simple example of a dynamical system. Since it can be difficult or impossible to give simple expressions of eigenvalues for larger matrices (i.e. 5 x 5 or higher), we then move towards specifying the power method for verified computation of eigenvalue approximations in PVS. This effort requires development of multivariate complex arithmetic. At the end of the day, having such additions to the PVS NASA libraries will help move towards the use of formal methods to verify concepts of control theory and system level verification.

Linear Algebra↗
DOE OSTI2020

A Shift Selection Strategy for Parallel Shift-invert Spectrum Slicing in Symmetric Self-consistent Eigenvalue Computation

The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A Shift Selection Strategy for Parallel Shift-Invert Spectrum Slicing in Symmetric Self-Consistent Eigenvalue Computation

The central importance of large scale eigenvalue problems in scientific computation necessitates the development massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world's largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing excessive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions which will be explored in future work.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Application of the Finite Element Method to Rotary Wing Aeroelasticity

A finite element method for the spatial discretization of the dynamic equations of equilibrium governing rotary-wing aeroelastic problems is presented. Formulation of the finite element equations is based on weighted Galerkin residuals. This Galerkin finite element method reduces algebraic manipulative labor significantly, when compared to the application of the global Galerkin method in similar problems. The coupled flap-lag aeroelastic stability boundaries of hingeless helicopter rotor blades in hover are calculated. The linearized dynamic equations are reduced to the standard eigenvalue problem from which the aeroelastic stability boundaries are obtained. The convergence properties of the Galerkin finite element method are studied numerically by refining the discretization process. Results indicate that four or five elements suffice to capture the dynamics of the blade with the same accuracy as the global Galerkin method.

Straub, F. K.↗
NASA NTRS1989

Computation of canonical correlation and best predictable aspect of future for time series

The canonical correlation between the (infinite) past and future of a stationary time series is shown to be the limit of the canonical correlation between the (infinite) past and (finite) future, and computation of the latter is reduced to a (generalized) eigenvalue problem involving (finite) matrices. This provides a convenient and essentially, finite-dimensional algorithm for computing canonical correlations and components of a time series. An upper bound is conjectured for the largest canonical correlation.

Pourahmadi, Mohsen↗
NASA NTRS1976

Determination of poles and zeros of transfer functions for flexible spacecraft attitude control

The transfer function matrix is obtained for a three-input and three-output model of minimum sensors and actuators for the attitude control system of flexible spacecraft, and a method is described for determining the poles and zeros of this transfer function. Three cases are considered: (1) the actuators and the sensors are all attached to the primary body, (2) the actuators are on the primary body and the sensors are on the sub-body, and (3) the actuators are on the sub-body and the sensors are on the primary body. The zero-determination problem is shown to reduce to eigenvalue calculations of a matrix which is constructed from the inertial and modal matrices in a simple fashion.

Ohkami, Y.↗
NASA NTRS1981

Instability of non-zonal baroclinic flows - Multiple-scale analysis

The linear instability of a non-zonal flow can be reduced to an eigenvalue-eigenfunction problem, governed by a nonseparable partial differential equation (Niehaus, 1980). Approximate solutions, found by the method of multiple scales, are derived here and compared with earlier results found using a spectral method. The amplitude maxima are correctly located. The zonal variations of local wavenumber and of amplitude are qualitatively correct, but not sufficiently extreme. Because the method is oversensitive to local conditions, and less sensitive to global constraints, this comparison provides theoretical limits to the possibility of parameterizing transient eddies in terms of the local time mean state of the atmosphere. The method can be extended easily to flows with more realistic vertical structure.

Niehaus, M. C. W.↗