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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Eigenspace design techniques for active flutter suppression

The application of eigenspace design techniques to an active flutter suppression system for the DAST ARW-2 research drone is examined. Eigenspace design techniques allow the control system designer to determine feedback gains which place controllable eigenvalues in specified configurations and which shape eigenvectors to achieve desired dynamic response. Eigenspace techniques were applied to the control of lateral and longitudinal dynamic response of aircraft. However, little was published on the application of eigenspace techniques to aeroelastic control problems. This discussion will focus primarily on methodology for design of full-state and limited-state (output) feedback controllers. Most of the states in aeroelastic control problems are not directly measurable, and some type of dynamic compensator is necessary to convert sensor outputs to control inputs. Compensator design are accomplished by use of a Kalman filter modified if necessary by the Doyle-Stein procedure for full-state loop transfer function recovery, by some other type of observer, or by transfer function matching.

Garrard, W. L.↗
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 NTRS1985

The exact eigenfunctions and eigenvalues of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics

Exact eigenfunctions for a two-dimensional rigid rotor are obtained using Gaussian wave packet dynamics. The wave functions are obtained by propagating, without approximation, an infinite set of Gaussian wave packets that collectively have the correct periodicity, being coherent states appropriate to this rotational problem. This result leads to a numerical method for the semiclassical calculation of rovibrational, molecular eigenstates. Also, a simple, almost classical, approximation to full wave packet dynamics is shown to give exact results: this leads to an a posteriori justification of the De Leon-Heller spectral quantization method.

Reimers, J. R.↗
NASA NTRS1988

Development of an unsteady wake theory appropriate for aeroelastic analyses of rotors in hover and forward flight

The purpose of this research is the development of an unsteady aerodynamic model for rotors such that it can be used in conventional aeroelastic analysis (e.g., eigenvalue determination and control system design). For this to happen, the model must be in a state-space formulation such that the states of the flow can be defined, calculated and identified as part of the analysis. The fluid mechanics of the problem is given by a closed-form inversion of an acceleration potential. The result is a set of first-order differential equations in time for the unknown flow coefficients. These equations are hierarchical in the sense that they may be truncated at any number of radial or azimuthal terms.

Peters, David A.↗
NASA NTRS1989

Preconditioning matrices for the pseudospectral approximation of first-order operators

The behavior of the eigenvalues of preconditioning matrices for the pseudospectral approximation to the derivative operator has been analyzed in one and two dimensions. The one-dimensional analysis resulted in real and positive eigenvalues for the selected tridiagonal matrices. In the two-dimensional analysis, the eigenvalues of the selected block-diagonal matrices behaved well, but the preconditioner is full and therefore not suitable for applications. The Richardson scheme has been applied in the unpreconditioned as well as the preconditioned version to find the solution of the model problem.

Funaro, D.↗
DOE OSTI2018

On semi-convergence and inexact iteration of the GSS iteration method for nonsymmetric singular saddle point problems

In this paper, the generalized shift-splitting (GSS) iteration method is used to solve nonsymmetric singular saddle point systems with the positive real (1, 1)-block sub-matrix. By a careful analysis of the spectral properties of the corresponding iteration matrix, we investigate the semi-convergence property of the GSS iteration method, give the sharper bounds of the eigenvalues for the GSS iteration method and propose the inexact GSS iteration method, also. Numerical experiments are given to show the correctness of the theoretical analysis and the feasibility of the GSS preconditioner with appropriate parameters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

MPACT Subgroup Self-Shielding Efficiency Improvements

Recent developments to improve the efficiency of the MOC solvers in MPACT have yielded effective kernels that loop over several energy groups at once, rather that looping over one group at a time. These kernels have produced roughly a 2x speedup on the MOC sweeping time during eigenvalue calculation. However, the self-shielding subgroup calculation had not been reevaluated to take advantage of these new kernels, which typically requires substantial solve time. The improvements covered in this report start by integrating the multigroup kernel concepts into the subgroup calculation, which are then used as the basis for further extensions. The next improvement that is covered is what is currently being termed as “Lumped Parameter MOC”. Because the subgroup calculation is a purely fixed source problem and multiple sweeps are performed only to update the boundary angular fluxes, the sweep procedure can be condensed to allow for the instantaneous propagation of the flux across a spatial domain, without the need to sweep along all segments in a ray. Once the boundary angular fluxes are considered to be converged, an additional sweep that will tally the scalar flux is completed. The last improvement that is investigated is the possible reduction of the number of azimuthal angles per octant in the shielding sweep. Typically 16 azimuthal angles per octant are used for self-shielding and eigenvalue calculations, but it is possible that the self-shielding sweeps are less sensitive to the number of angles than the full eigenvalue calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Lieb-Liniger-like model of quantum solvation in CO-{sup 4}He{sub N} clusters

Small {sup 4}He clusters doped with various molecules allow for the study of “quantum solvation” as a function of cluster size. A peculiarity of quantum solvation is that, as the number of {sup 4}He atoms is increased from N = 1, the solvent appears to decouple from the molecule which, in turn, appears to undergo free rotation. This is generally taken to signify the onset of “microscopic superfluidity.” Currently, little is known about the quantum mechanics of the decoupling mechanism, mainly because the system is a quantum (N + 1)-body problem in three dimensions which makes computations difficult. Here, a one-dimensional model is studied in which the {sup 4}He atoms are confined to revolve on a ring and encircle a rotating CO molecule. The Lanczos algorithm is used to investigate the eigenvalue spectrum as the number of {sup 4}He atoms is varied. Substantial solvent decoupling is observed for as few as N = 5 {sup 4}He atoms. Examination of the Hamiltonian matrix, which has an almost block diagonal structure, reveals increasingly weak inter-block (solvent-molecule) coupling as the number of {sup 4}He atoms is increased. In the absence of a dopant molecule the system is similar to a Lieb-Liniger (LL) gas and we find a relatively rapid transition to the LL limit as N is increased. In essence, the molecule initially—for very small N—provides a central, if relatively weak, attraction to organize the cluster; as more {sup 4}He atoms are added, the repulsive interactions between the identical bosons start to dominate as the solvation ring (shell) becomes more crowded which causes the molecule to start to decouple. For low N, the molecule pins the atoms in place relative to itself; as N increases the atom-atom repulsion starts to dominate the Hamiltonian and the molecule decouples. We conclude that, while the notion of superfluidity is a useful and correct description of the decoupling process, a molecular viewpoint provides complementary insights into the quantum mechanism of the transition from a molecular cluster to a quantum solvated molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Vibration control of large linear quadratic symmetric systems

Some unique properties on a class of the second order lambda matrices were found and applied to determine a damping matrix of the decoupled subsystem in such a way that the damped system would have preassigned eigenvalues without disturbing the stiffness matrix. The resulting system was realized as a time invariant velocity only feedback control system with desired poles. Another approach using optimal control theory was also applied to the decoupled system in such a way that the mode spillover problem could be eliminated. The procedures were tested successfully by numerical examples.

Jeon, G. J.↗
NASA NTRS1993

Structural/control system optimization with variable actuator masses

A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.

Jin, Ik M.↗
DOE OSTI2019

Adaptive Reweighted Variance Estimation for Monte Carlo Eigenvalue Simulations

Monte Carlo (MC) simulation is used to solve the eigenvalue form of the Boltzmann transport equation to estimate various parameters such as fuel pin flux distributions that are crucial for the safe and efficient operation of nuclear systems (e.g., a nuclear reactor). Monte Carlo eigenvalue simulation uses a sample mean over many stationary cycles (iterations) to estimate various parameters important to nuclear systems. A variance estimate of the sample mean is often used for calculating the confidence intervals. However, MC eigenvalue simulation variance estimators that ignore the intercycle correlation underestimate the true variance of the estimated quantity. This paper presents novel data-adaptive approaches based on a simple autoregressive (AR) model and sigmoid functions to improve MC variance estimation. The standard MC sample-based variance estimator (or naïve estimator) and the spectral density–based MC variance estimator are enhanced by adding data-adaptive components that reduce their bias and improve performance. By investigating the frequency pattern of the AR(1) (order 1) model, two adaptive spectral estimators and one adaptive naïve estimator are proposed. The proposed estimators manifest superior performance when applied to three test problems compared to the standard spectral density–based estimator previously introduced by the authors. In conclusion, these new estimators are straightforward, as they use online algorithms and do not require storage of tallies from all active cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations

We study Chebyshev filter diagonalization as a tool for the computation of many interior eigenvalues of very large sparse symmetric matrices. In this technique the subspace projection onto the target space of wanted eigenvectors is approximated with filter polynomials obtained from Chebyshev expansions of window functions. After the discussion of the conceptual foundations of Chebyshev filter diagonalization we analyze the impact of the choice of the damping kernel, search space size, and filter polynomial degree on the computational accuracy and effort, before we describe the necessary steps towards a parallel high-performance implementation. Because Chebyshev filter diagonalization avoids the need for matrix inversion it can deal with matrices and problem sizes that are presently not accessible with rational function methods based on direct or iterative linear solvers. To demonstrate the potential of Chebyshev filter diagonalization for large-scale problems of this kind we include as an example the computation of the 10{sup 2} innermost eigenpairs of a topological insulator matrix with dimension 10{sup 9} derived from quantum physics applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

BWR Progression Problems

Under the Consortium for the Advanced Simulation of Light Water Reactors (CASL) Program, the Virtual Environment for Reactor Applications (VERA) was developed with the primary focus to model pressurized water reactors (PWRs). Recently, a new project was started to extend the modeling capability in VERA to model boiling water reactors (BWRs). The new project is called “Modeling and Analysis of Exelon BWRs for Eigenvalue and Thermal Limits Predictability,” and it is led by Oak Ridge National Laboratory (ORNL) and Exelon with participation from Global Nuclear Fuel and three universities: North Carolina State University (NCSU), the University of Michigan, and the University of Illinois.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Layered CAD/CSG geometry for spatially complex radiation transport scenarios

Many spatially complex fission, fusion, and national security Monte Carlo (MC) radiation transport scenarios involve combining computer-aided design (CAD) models with constructive solid geometry (CSG) models. A layered geometry method has been implemented in the Shift MC code to address this need. With layered geometry, multiple CAD and/or CSG models can be clipped, translated, rotated, and placed in overlapping layers to form transport-ready geometries. Here, the utility of this method is demonstrated with two problems: (1) a fixed-source simulation with a layered geometry consisting of a LiDAR-generated CAD model of the Combined Arms Collective Training Facility urban environment overlaid with CSG models of a mock hotel and a detector apparatus, and (2) a k-eigenvalue calculation using a layered geometry model of the Transformational Challenge Reactor consisting of CAD fuel elements placed in a CSG core. Tallied particle flux distributions match expectations, but tracking robustness must be improved prior to general-purpose use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 NTRS1980

A method of order reduction for structural dynamics

This paper presents a new method of order reduction based on recent work of a similar nature applicable to system dynamics and control. This method provides a practical computational procedure for producing a condensed model which exactly preserves the slowest n1 modes (where n1 is the number of 'slow' eigenvalues) of the total n modes of the original model for almost any set of n1 degrees of freedom retained in the condensed model. The method can also be used to compute the eigensolutions corresponding to the n1 slowest modes of the original structural dynamics problem. The accuracy of the condensed model and the speed/accuracy performance of the eigensolver are compared with standard methods for a 90 DOF cantilevered plate.

Anderson, L. R.↗
NASA NTRS1992

Adaptive control of Space Station with control moment gyros

An adaptive approach to Space Station attitude control is investigated. The main components of the controller are the parameter identification scheme, the control gain calculation, and the control law. The control law is a full-state feedback space station baseline control law. The control gain calculation is based on linear-quadratic regulator theory with eigenvalues placement in a vertical strip. The parameter identification scheme is a recursive extended Kalman filter that estimates the inertias and also provides an estimate of the unmodeled disturbances due to the aerodynamic torques and to the nonlinear effects. An analysis of the inertia estimation problem suggests that it is possible to estimate Space Station inertias accurately during nominal control moment gyro operations. The closed-loop adaptive control law is shown to be capable of stabilizing the Space Station after large inertia changes. Results are presented for the pitch axis.

Bishop, Robert H.↗
NASA NTRS2010

Broken Ergodicity in MHD Turbulence

Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.

Shebalin, John V.↗