Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A general method for dynamic analysis of structures overview

The presented research deals with the development of a dynamic analysis method for structural systems. The modeling approach is essentially a finite element method in the sense that the structure is divided into n elements. An element is defined as any structural unit whose degree of freedom (dofs) can be categorized as either interface or non-interface dofs. An element could be a fundamental unit such as a rod, a beam, a plate etc., or it could be an entire structural component. Furthermore, the parameters for the element could be distributed or lumped. The choice of elements is totally arbitrary and is a matter of user convenience. In particular, issues of accuracy and convergence do not enter on the level of example that bookkeeping is reduced to a minimum. Each element is modeled using a set of interface constraint modes (ICM) combined with a set of interface restrained normal models (IRNM). The next step is the solution of the system eigenvalue problem. The procedure calls for the sequential solution of a number of small eigenvalue problems based on a truncation principle for IRNM. In addition, the form of these eigenvalue problems is very simple such that an escalator type of eigenvalue problem solver can be used which is extremely cost-effective and fast.

Engels, Remi C.↗
NASA NTRS2002

Linear Temporal Stability Analysis of a Low-Density Round Gas Jet Injected into a High-Density Gas

It has been observed in previous experimental studies that round helium jets injected into air display a repetitive structure for a long distance, somewhat similar to the buoyancy-induced flickering observed in diffusion flames. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis of a round helium jet injected into air was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. The temporal growth rates of the disturbances increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity.

Lawson, Anthony L.↗
NASA NTRS1992

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗
NASA NTRS1981

Graph theory approach to the eigenvalue problem of large space structures

Graph theory is used to obtain numerical solutions to eigenvalue problems of large space structures (LSS) characterized by a state vector of large dimensions. The LSS are considered as large, flexible systems requiring both orientation and surface shape control. Graphic interpretation of the determinant of a matrix is employed to reduce a higher dimensional matrix into combinations of smaller dimensional sub-matrices. The reduction is implemented by means of a Boolean equivalent of the original matrices formulated to obtain smaller dimensional equivalents of the original numerical matrix. Computation time becomes less and more accurate solutions are possible. An example is provided in the form of a free-free square plate. Linearized system equations and numerical values of a stiffness matrix are presented, featuring a state vector with 16 components.

Reddy, A. S. S. R.↗
NASA NTRS1984

A solution of Mx(double dot) + Cx(dot) + Kx = 0 applicable to the design of active dampers

A solution is presented for the equations of motion for the damped linear oscillator, Mx(double dot) + Cx(dot) + Kx = 0. The algorithm solves a transformed set of equations in terms of the modal variables of the undamped system and, at the same time, solves the adjoint equation of the transformed problem. The adjoint solution is normalized to give the inverse of the solution matrix of the transformed problem. The normalized inverse is useful in design for direct computation of sensitivity derivatives of damping ratios with respect to damping rates. The algorithm is programmed to reduce storage requirements by a factor of three-fourths compared to standard complex eigenvalue subroutines. A numerical example is included.

Thurston, G. A.↗
NASA NTRS1992

Minimization of the vibration energy of thin-plate structure

An optimization method is proposed to reduce the vibration of thin plate structures. The method is based on a finite element shell analysis, a modal analysis, and a structural optimization method. In the finite element analysis, a triangular shell element with 18 dof is used. In the optimization, the overall vibration energy of the structure is adopted as the objective function, and it is minimized at the given exciting frequency by varying the thickness of the elements. The technique of modal analysis is used to derive the sensitivity of the vibration energy with respect to the design variables. The sensitivity is represented by the sensitivities of both eigenvalues and eigenvectors. The optimum value is computed by the gradient projection method and a unidimensional search procedure under the constraint condition of constant weight. A computer code, based on the proposed method, is developed and is applied to design problems using a beam and a plate as test cases. It is confirmed that the vibration energy is reduced at the given exciting frequency. For the beam excited by a frequency slightly less than the fundamental natural frequency, the optimized shape is close to the beam of uniform strength.

Inoue, Katsumi↗
NASA NTRS1972

A fast numerical method for determining the optimum SNR of an array subject to a Q factor constraint.

The problem of finding the optimum current excitations for a point-detector array which extracts a signal from background noise using the criterion of maximizing the SNR subject to a constraint on the supergain ratio (or Q factor) has recently been solved. The formulation reduced the problem to one of solving for the roots of a polynomial whose coefficients are, in general, complex. We demonstrate that by making a certain transformation, it becomes possible to employ a state variable technique which changes the problem into one of finding the eigenvalues of a real matrix which is considerably faster to do.

Winkler, L. P.↗
NASA NTRS1986

Symmetry and stability in Taylor-Couette flow

The flow of a fluid between concentric rotating cylinders (the Taylor problem) is studied by exploiting the symmetries of the system. The Navier-Stokes equations, linearized about Couette flow, possess two zero and four purely imaginary eigenvalues at a suitable value of the speed of rotation of the outer cylinder. There is thus a reduced bifurcation equation on a six-dimensonal space which can be shown to commute with an action of the symmetry group 0(2) x S0(2). The group structure is used to analyze this bifurcation equation in the simplest (nondegenerate) case, and to compute the stabilities of solutions. In particular, when the outer cylinder is counterrotated, transitions which seem to agree with recent experiments of Andereck, Liu, and Swinney (1984) are obtained. It is also possible to obtain the 'main sequence' in this model. This sequence is normally observed in experiments when the outer cylinder is held fixed.

Golubitsky, M.↗
NASA NTRSDate not supplied

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗
DOE OSTI2025

Entanglement Witness for Indistinguishable Electrons Using Solid-State Spectroscopy

Characterizing entanglement in quantum materials is crucial for advancing next-generation quantum technologies. Despite recent strides in witnessing entanglement in magnetic materials with distinguishable spin modes, quantifying entanglement in systems formed by indistinguishable electrons remains a formidable challenge. To solve this problem, we introduce a method to extract various four-fermion correlations by analyzing the nonlinearity in resonant inelastic x-ray scattering spectra. These correlations constitute the primary components of the cumulant two-particle reduced density matrix. We further derive bounds for its eigenvalues and demonstrate the linear scaling with fermionic entanglement depth, providing a reliable witness for entanglement. Using the material-relevant strongly correlated models as examples, we show how this entanglement witness can efficiently quantify multipartite entanglement across different phase regions, highlighting its advantage over quantum Fisher information. Published by the American Physical Society 2025

Liu, Tongtong (ORCID:0000000295324061)↗
NASA NTRS1973

Optical potential approach to the electron-atom impact ionization threshold problem

The problem of the threshold law for electron-atom impact ionization is reconsidered as an extrapolation of inelastic cross sections through the ionization threshold. The cross sections are evaluated from a distorted wave matrix element, the final state of which describes the scattering from the Nth excited state of the target atom. The actual calculation is carried for the e-H system, and a model is introduced which is shown to preserve the essential properties of the problem while at the same time reducing the dimensionability of the Schrodinger equation. Nevertheless, the scattering equation is still very complex. It is dominated by the optical potential which is expanded in terms of eigen-spectrum of QHQ. It is shown by actual calculation that the lower eigenvalues of this spectrum descend below the relevant inelastic thresholds; it follows rigorously that the optical potential contains repulsive terms. Analytical solutions of the final state wave function are obtained with several approximations of the optical potential.

Temkin, A.↗
NASA NTRS1973

Divergence and Necessary Conditions for Extremums

The problem is considered of finding a dimension reducing transformation matrix B that maximizes the divergence in the reduced dimension for multi-class cases. A comparitively simple expression for the gradient of the average divergence with respect to B is developed. The developed expression for the gradient contains no eigenvectors or eigenvalues; also, all matrix inversions necessary to evaluate the gradient are available from computing the average divergence.

Quirein, J. A.↗
NASA NTRS1987

Improvements to the fastex flutter analysis computer code

Modifications to the FASTEX flutter analysis computer code (UDFASTEX) are described. The objectives were to increase the problem size capacity of FASTEX, reduce run times by modification of the modal interpolation procedure, and to add new user features. All modifications to the program are operable on the VAX 11/700 series computers under the VAX operating system. Interfaces were provided to aid in the inclusion of alternate aerodynamic and flutter eigenvalue calculations. Plots can be made of the flutter velocity, display and frequency data. A preliminary capability was also developed to plot contours of unsteady pressure amplitude and phase. The relevant equations of motion, modal interpolation procedures, and control system considerations are described and software developments are summarized. Additional information documenting input instructions, procedures, and details of the plate spline algorithm is found in the appendices.

Taylor, Ronald F.↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2020

Survey of Dynamic Mode Decomposition Methods

Dynamic mode decomposition (DMD) is a data-driven reduced order modeling (ROM) technique used for dynamic systems. The widely adopted algorithm was first introduced and demonstrated on fluid flow data by Schmid. In recent years, various other fields, such as nuclear engineering, have begun to adopt this method. For example, DMD has been used for estimating α-eigenvalues, as an ROM for pulsed neutron problems, for predicting isotopic composition in burnup calculations, as acceleration techniques for iterative methods, and in capturing dynamic behaviors in molten salt reactor transients. This report seeks to demonstrate the capabilities and limits of the standard DMD algorithm, and identify problem spaces where variants may be better suited. The primary variant this report considers is Multi-Resolution DMD (mrDMD). Because this serves as a survey, synthetically produced data is used in lieu of simulation results. The remainder of this report will go into detail on the DMD theory, outline the standard DMD and mrDMD algorithms, present test cases highlighting the applicability of each, and finally present a discussion on how to determine the best suited algorithm for a given problem. All calculations performed in this report are carried out using the open source DMD library, PyDMD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary layer flow (i.e., the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so called doubly generalized inflexion condition is derived, which is a condition for the existence of so called subsonic neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two freestream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to rapidly disappear as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗
NASA NTRS1973

An automated procedure for computing flutter eigenvalues.

A new, fast and economical automated procedure for implementing the traditional V-g method of flutter solution is described. The procedure requires as input the generalized aerodynamic forces for a range of reduced frequencies obtained from an aerodynamic program. These aerodynamic forces are interpolated with respect to reduced frequency using a newly developed, partially tabulated cubic spline that is both fast in execution and economical in storage. The flutter solution is then obtained using an eigenvalue routine that has been developed to take advantage of the parametric nature of the V-g type of solution. Furthermore, the routine takes care of the fundamental and troublesome problem of properly sorting the output eigenvalues. By solving the root-sorting problem, the interpolation for flutter crossings and automatic plotting are accomplished efficiently. The computational techniques used in this new program are described and some sample results are given.

Desmarais, R. N.↗
NASA NTRS1990

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary-layer flow (i.e. the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so-called 'doubly generalized' inflexion condition is derived, which is a condition for the existence of so-called 'subsonic' neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two free-stream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to disappear rapidly as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗