SEARCH · Search NASA

Results for “eigenvalue analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1996

The Effects of Acoustic Treatment on Pressure Disturbances From a Supersonic Jet in a Circular Duct

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined circular duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions at each axial location. These values are used to determine the growth rates and phase velocities of the instability waves and the near field pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, and wall impedance, and the effects of these changes on non-axisymmetric modes. The primary results indicate that the effects of the duct wall on stability characteristics diminish as the outer flow increases and as the jet azimuthal mode number increases. Also, wall reflections are reduced when using a finite impedance boundary condition at the wall; but in addition, reflections are reduced and growth rates diminished by keeping the imaginary part of the impedance negative when using the negative exponential for the harmonic dependence.

Dahl, Milo D.↗

NASA NTRS1997

Effects of Acoustically Lined Cylindrical Ducts on Instability Waves in Confined Supersonic Jets

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined cylindrical duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions. These values are used to determine the growth rates and phase velocities of the instability waves and the radial pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, wall impedance, and frequency. Results indicate that the effects of the duct wall on shear layer growth rates diminish as the outer flow increases. Also, wall reflections cause variations in growth rates depending on wall height and Strouhal number. These variations are due to the phase relationship between the outgoing and the reflected incoming pressure disturbances at the shear layer. The growth rate variations can be reduced and the maximum growth rate minimized by keeping the imaginary part of the impedance negative.

Dahl, Milo D.↗

DOE OSTI2016

A Methodology for Loading the Advanced Test Reactor Driver Core for Experiment Analysis

In support of experiments in the ATR, a new methodology was devised for loading the ATR Driver Core. This methodology will replace the existing methodology used by the INL Neutronic Analysis group to analyze experiments. Studied in this paper was the as-run analysis for ATR Cycle 152B, specifically comparing measured lobe powers and eigenvalue calculations.

42 ENGINEERING↗

DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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

NASA NTRS1996

Analysis of the Hessian for Aerodynamic Optimization: Inviscid Flow

In this paper we analyze inviscid aerodynamic shape optimization problems governed by the full potential and the Euler equations in two and three dimensions. The analysis indicates that minimization of pressure dependent cost functions results in Hessians whose eigenvalue distributions are identical for the full potential and the Euler equations. However the optimization problems in two and three dimensions are inherently different. While the two dimensional optimization problems are well-posed the three dimensional ones are ill-posed. Oscillations in the shape up to the smallest scale allowed by the design space can develop in the direction perpendicular to the flow, implying that a regularization is required. A natural choice of such a regularization is derived. The analysis also gives an estimate of the Hessian's condition number which implies that the problems at hand are ill-conditioned. Infinite dimensional approximations for the Hessians are constructed and preconditioners for gradient based methods are derived from these approximate Hessians.

Arian, Eyal↗

DOE OSTI2016

Hot zero power reactor calculations using the Insilico code

In this paper we describe the reactor physics simulation capabilities of the Insilico code. A description of the various capabilities of the code is provided, including detailed discussion of the geometry, meshing, cross section processing, and neutron transport options. Numerical results demonstrate that Insilico using an SP{sub N} solver with pin-homogenized cross section generation is capable of delivering highly accurate full-core simulation of various pressurized water reactor problems. Comparison to both Monte Carlo calculations and measured plant data is provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

Linear Static Structural and Vibration Analysis on High-Performance Computers

Parallel computers offer the opportunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on a massively-parallel computers hereafter referred to as Scalable High Performance Computers (SHPC) for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. Models of high speed civil transport). The goal of this research is to develop new efficient technique which extend structural analysis to SHPC and make large-scale structural analyses tractable.

Baddourah, Majdi↗

NASA NTRS1993

Linear static structural and vibration analysis on high-performance computers

Parallel computers offer the oppurtunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on massively-parallel computers hereafter referred to as Scalable High-Performance Computers (SHPC), for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. models for High-Speed Civil Transport). The goal of this research is to develop a new, efficient technique which extends structural analysis to SHPC and makes large-scale structural analyses tractable.

Baddourah, M. A.↗

DOE OSTI2024

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Testing the Activation Analysis for Fusion in OpenMC

OpenMC is a community-developed Monte Carlo neutron and photon transport simulation code. It can perform fission simulations such as fixed-source, k-eigenvalue, and subcritical multiplication calculations on models built using either a constructive solid geometry or CAD representation. To explore the use of OpenMC for fusion activation analysis, a detailed model of the Fusion Neutronics Science Facility (FNSF) was first developed for comparisons against an existing SERPENT model. A 90-degree model of FNSF in Standard-Triangle-Language (STL) CAD format was converted to Constructive Solid Geometry (CSG) using each code's built-in functions, and the geometries were validated by ensuring no cells overlapped and no particles were lost during simulations. The neutron fluxes were calculated and compared for multiple components close to the plasma. The results show differences mostly below 1% in fluxes and averaged 8% for activity and decay heat. Here, the work described in this study tests the CAD-based geometry using the DagMC toolkit in OpenMC and compares the activation analysis of OpenMC to SERPENT code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1988

A spectral collocation solution to the compressible stability eigenvalue problem

A newly developed spectral compressible linear stability code (SPECLS) (staggered pressure mesh) is presented for analysis of shear flow stability, and applied to high speed boundary layers and free shear flows. The formulation utilizes the first application of a staggered mesh for a compressible flow analysis by a spectral technique. An order of magnitude less number of points is needed for equivalent accuracy of growth rates compared to those calculated by a finite difference formulation. Supersonic disturbances which are found to have oscillatory structures were resolved by a spectral multi-domain discretization, which requires a factor of three fewer points than the single domain spectral stability code. It is indicated, as expected, that stability of mixing layers is enhanced by viscosity and increasing Mach number. The mean flow involves a jet being injected into a quiescent gas. Higher temperatures of the injected gas is also found to enhance stability characteristics of the free shear layer.

Macaraeg, Michele G.↗

DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗

NASA NTRS1975

Stability and transient motion of a plain journal mounted in flexible damped supports

Results are presented for an extensive study of the influence of damped flexible supports on the stability-threshold speed of a symmetric rotor mounted in nonlinear fluid-film bearings. A stability analysis is conducted by applying the Routh criteria and calculating the damped eigenvalues of the system. Stability maps are given in dimensionless form for a range of support stiffness, damping, and mass ratios. The effect of rotor imbalance on the stability of rotor systems is examined, and it is shown that time transient orbits of shaft and support motions clearly indicate their dynamic behavior for both stable and unstable operating speeds. Simple design criteria are proposed for choosing the support stiffness and damping on the basis of rotor weight and journal clearance to promote optimum stability. It is concluded that the use of a flexible damped support system may increase the stability threshold of a plain journal bearing and that damper supports should promote smoother and quieter machine operation by suppressing oil whirl and attenuating unbalanced rotational forces.

Kirk, R. G.↗

NASA NTRS1999

Aerodynamic and Aeroelastic Insights using Eigenanalysis

This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Heeg, Jennifer↗

NASA NTRS2004

Aerodynamic and Aeroelastic Insights using Eigenanalysis

Heeg, Jennifer↗

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↗