Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Eigensensitivity analysis for space structures with applications in design optimization

An approximate structural eigenvalue/vector analysis technique which uses eigensensitivities in a truncated Taylor series expression is presented. It is shown that this technique can provide the computational efficiency urgently needed for large scale control-structure optimization problems. In addition, a unified formulation and solution approach for calculating eigenvalue and eigenvector derivatives of the real symmetric structural eigenvalue problem is presented.

Kenny, Sean P.↗
NASA NTRS1987

Stability of streamwise vortices

A brief overview of some theoretical and computational studies of the stability of streamwise vortices is given. The local induction model and classical hydrodynamic vortex stability theories are discussed in some detail. The importance of the three-dimensionality of the mean velocity profile to the results of stability calculations is discussed briefly. The mean velocity profile is provided by employing the similarity solution of Donaldson and Sullivan. The global method of Bridges and Morris was chosen for the spatial stability calculations for the nonlinear eigenvalue problem. In order to test the numerical method, a second order accurate central difference scheme was used to obtain the coefficient matrices. It was shown that a second order finite difference method lacks the required accuracy for global eigenvalue calculations. Finally the problem was formulated using spectral methods and a truncated Chebyshev series.

Khorrami, M. K.↗
DOE OSTI2018

Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method

For a nuclear system in which the entire α-eigenvalue spectrum is known, eigenfunction expansion yields the time-dependent flux response to any arbitrary source. Applications in which this response is of interest include pulsed-neutron experiments, accelerator-driven subcritical systems, and fast burst reactors, where a steady-state assumption used in neutron transport is invalid for characterizing the time-dependent flux. To obtain the α-eigenvalue spectrum, the transition rate matrix method (TRMM) tallies transition rates describing neutron behavior in a discretized position-direction-energy phase space using Monte Carlo. Interpretation of the resulting Markov process transition rate matrix as the operator in the adjoint α-eigenvalue problem provides an avenue for determining a large finite set of α eigenvalues and eigenfunctions of a nuclear system. Results from the TRMM are verified using analytic solutions, time-dependent Monte Carlo simulations, and modal expansion from diffusion theory. For simplified infinite-medium and one-dimensional geometries, the TRMM accurately calculates eigenvalues, eigenfunctions, and eigenfunction expansion solutions. Applications and comparisons to measurements are made for the small fast burst reactor CALIBAN and the Fort St. Vrain high-temperature gas-cooled reactor. For large three-dimensional geometries, discretization of the large position-energy-direction phase space limits the accuracy of eigenfunction expansion solutions using the TRMM, but it can still generate a fair estimate of the fundamental eigenvalue and eigenfunction. These results show that the TRMM generates an accurate estimate of a large number of α eigenvalues. This is not possible with existing Monte Carlo–based methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Investigating a hybrid perturbation-Galerkin technique using computer algebra

A two-step hybrid perturbation-Galerkin method is presented for the solution of a variety of differential equations type problems which involve a scalar parameter. The resulting (approximate) solution has the form of a sum where each term consists of the product of two functions. The first function is a function of the independent field variable(s) x, and the second is a function of the parameter lambda. In step one the functions of x are determined by forming a perturbation expansion in lambda. In step two the functions of lambda are determined through the use of the classical Bubnov-Gelerkin method. The resulting hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Bubnov-Galerkin methods applied separately, while combining some of the good features of each. In particular, the results can be useful well beyond the radius of convergence associated with the perturbation expansion. The hybrid method is applied with the aid of computer algebra to a simple two-point boundary value problem where the radius of convergence is finite and to a quantum eigenvalue problem where the radius of convergence is zero. For both problems the hybrid method apparently converges for an infinite range of the parameter lambda. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Andersen, Carl M.↗
NASA NTRS1985

On the design derivatives of eigenvalues and eigenvectors for distributed parameter systems

In this paper, analytic expressions are obtained for the design derivatives of eigenvalues and eigenfunctions of self-adjoint linear distributed parameter systems. Explicit treatment of boundary conditions is avoided by casting the eigenvalue equation into integral form. Results are expressed in terms of the linear operators defining the eigenvalue problem, and are therefore quite general. Sufficiency conditions appropriate to structural optimization of eigenvalues are obtained.

Reiss, R.↗
NASA NTRS1992

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗
NASA NTRS1966

What is H sub 0/ques/

Zero order Hamiltonian in perturbation theory for eigenvalue problems of atom in magnetic field

EIGENVALUE↗
DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗
NASA NTRS1981

Eigenvalues of the Rayleigh-Benard and Marangoni problems

The eigenvalues of the linear Benard-Marangoni stability problem are discussed. Pearson and Nield boundary conditions, which correspond to a rigid, isothermal lower boundary and a stress-free conducting upper boundary are considered. It is shown that although a critical value of the Marangoni number can be determined, the number is not, strictly speaking, an eigenvalue and cannot be used as an eigenvalue parameter for the determination of an eigenvector set.

Rosenblat, S.↗
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 OSTI2019

Problem Description

Gobstopper.inp is a test problem that is designed to be a nice middle-of-the-road problem that balances work done in tracking and collisions. It is domain decomposed to stress particle communication between domains while in the tracking phase. This is an eigenvalue problem that is forced to run for 25 cycles with 10,000 particles per cycle. The final eigenvalue (keff) should be around 0.998 with some uncertainty reported. The cross sections are stored in a 230 energy multigroup fashion, but the collisions use continuous energy physics. The geometry is 1,000 nested spheres made of random materials from the list of materials found in a reactor. The shells are uniformly spaced up to 10 centimeters. However, 80% of the shells are made of highly enriched uranium found in the Godiva criticality experiment. This was to push the criticality eigenvalue (keff) closer to 1.0. The material compositions may be found in the material section. The densities are all the same for each material and can be found in the cell description. The geometry is further decomposed into 64 (4 X 4 X 4) cube domains. The problem is initially sourced in the first cycle with particles uniformly over a 10 cm sphere. Their energy spectrum is represented by the outgoing fission spectrum from U-235 nucleus hit by a 1.5 MeV neutron. The tallies are mostly turned off except for the balance, eigenvalue and spectrum tallies. The spectrum and eigenvalue tallies are written out to the .plt file. The balance and eigenvalue tallies are written to the .out file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

The two-dimensional kinetic ballooning theory for trapped electron mode in tokamak

The two-dimensional (2D) kinetic theory for a collisionless trapped electron mode is developed based on the Fourier-ballooning transform in an up-down symmetric equilibrium (illustrated via concentric circular magnetic surfaces). The system consists of two equations: the ballooning (integral) equation with a parameterized Floquet phase and a second order differential equation for the distribution of the Floquet phase. The coupled equations are, then, numerically solved as an eigenvalue problem yielding the 2D mode structure (in real space) as well as the global (phase-independent) eigenvalue for an L-mode parameter set. The 2D mode structure exhibits apparent radial-poloidal asymmetry; due to the poloidal coupling, the radial correlation length is found to be, at least, twice as large as the poloidal one. The global (phase-independent) eigenvalue of the mode differs considerably from the conventional local (phase-dependent) estimate. Finally, this paper shares many technical aspects with a published paper that works out the 2D kinetic theory for the ion temperature gradient mode [Xie et al., Phys. Plasmas 24, 102506 (2017)]

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Sensitivity analysis for large-scale problems

The development of efficient techniques for calculating sensitivity derivatives is studied. The objective is to present a computational procedure for calculating sensitivity derivatives as part of performing structural reanalysis for large-scale problems. The scope is limited to framed type structures. Both linear static analysis and free-vibration eigenvalue problems are considered.

Noor, Ahmed K.↗
NASA NTRS1994

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗
NASA NTRS1980

Linear state feedback, quadratic weights, and closed loop eigenstructures

Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalue and the directional derivatives of closed loop eigenvectors. An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties. An algorithm is presented that can be used to select a feedback gain matrix for the linear state feedback problem which produces a specified asymptotic eigenstructure. Another algorithm is given to compute the asymptotic eigenstructure properties inherent in a given set of quadratic weights. Finally, it is shown that optimal root loci for nongeneric problems can be approximated by generic ones in the nonasymptotic region.

Thompson, P. M.↗
NASA NTRS1972

The calculation of the eigenvalues and eigenfunctions of Mathieu's equation

The eigenfunctions of Mathieu's equation are expanded in trigonometric series, and the resulting eigenvalue problem is cast in matrix form. This matrix is found to be a symmetric, triagonal matrix, and the eigenvalues are computed using the bisection method. The eigenfunction expansion coefficients are obtained by the standard recursion method. This computational technique for the eigenvalues and eigenfunctions of Mathieu's equation is both rapid and accurate.

Hodge, D. B.↗
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.↗
NASA NTRS1987

Removal of spurious modes encountered in solving stability problems by spectral methods

A technique based on the Galerkin approximation is developed to remove spurious roots arising when Chebyshev spectral methods are used to solve eigenvalue problems in hydrodynamic stability. The derivation of Galerkin-Chebyshev approximations is explained, and numerical results for the Orr-Sommerfeld equations of plane Poiseuille flow and a Blasius profile are presented in tables and compared with those obtained by the method of Zebib (1984). It is pointed out that the present method does not increase the size of the algebraic system to be solved.

Zebib, Abdelfattah↗