Results for “eigenvalue analysis” — Page 15

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

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 NTRS1983

A modal analysis of flexible aircraft dynamics with handling qualities implications

A multivariable modal analysis technique is presented for evaluating flexible aircraft dynamics, focusing on meaningful vehicle responses to pilot inputs and atmospheric turbulence. Although modal analysis is the tool, vehicle time response is emphasized, and the analysis is performed on the linear, time-domain vehicle model. In evaluating previously obtained experimental pitch tracking data for a family of vehicle dynamic models, it is shown that flexible aeroelastic effects can significantly affect pitch attitude handling qualities. Consideration of the eigenvalues alone, of both rigid-body and aeroelastic modes, does not explain the simulation results. Modal analysis revealed, however, that although the lowest aeroelastic mode frequency was still three times greater than the short-period frequency, the rigid-body attitude response was dominated by this aeroelastic mode. This dominance was defined in terms of the relative magnitudes of the modal residues in selected vehicle responses.

Schmidt, D. K.↗

NASA NTRS1995

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗

NASA NTRS1989

The application of TSIM software to act design and analysis on flexible aircraft

The TSIM software is described. This is a package which uses an interactive FORTRAN-like simulation language for the simulation on nonlinear dynamic systems and offers facilities which include: mixed continuous and discrete time systems, time response calculations, numerical optimization, automatic trimming of nonlinear aircraft systems, and linearization of nonlinear equations for eigenvalues, frequency responses and power spectral response evaluation. Details are given of the application of TSIM to the analysis of aeroelastic systems under the RAE Farborough extension FLEX-SIM. The aerodynamic and structural data for the equations of motion of a flexible aircraft are prepared by a preprocessor program for incorporation in TSIM simulations. Within the simulation, the flexible aircraft model may then be selected interactively for different flight conditions and modal reduction techniques applied. The use of FLEX-SIM is demonstrated by an example of the flutter prediction for a simple aeroelastic model. By utilizing the numerical optimization facility of TSIM, it is possible to undertake identification of required parameters in the TSIM model within the simulation.

Kaynes, Ian W.↗

DOE OSTI2025

Multi-Objective Optimization of Uranium Target Assembly–3: A Comparison of Genetic and Traditional Methods

Commonly produced as a byproduct of uranium fission, 99 Mo is a key medical isotope that is in high demand in the United States. An international goal is to switch from medical isotope production technologies that require highly enriched uranium to medical isotope production technologies that require only low-enriched uranium. Niowave Inc. is contributing to this goal by developing an accelerator-driven subcritical assembly called the Uranium Target Assembly (UTA). This work compares the performance of Dakota’s Multi-Objective Genetic Algorithm (MOGA) against traditional sensitivity analysis in the neutronic optimization of the UTA-3 system. The design objectives are k-eigenvalue (k eff ) and natural uranium fission power, which are directly correlated with the amount of 99 Mo produced. Dakota:MOGA did not perform as well as human engineering ingenuity in optimization studies with high numbers of input parameters, such as fuel rod type selection and fuel rod placement. However, Dakota:MOGA did outperform traditional sensitivity analysis in optimization studies with fewer than 20 parameters and revealed the degree to which each parameter influences the optimal design space for k eff and natural uranium fission power (to a lesser extent). As the design model became more complex in the final stage of design, the computational resources required to calculate the design objective values in the Monte Carlo N-Particle transport code from selected input parameter combinations limited Dakota:MOGA’s performance, and, unfortunately, human intervention was required to discern the optimal design space. In conclusion, future work will attempt to reduce computational resource constraints by incorporating areduced-order neutronics model into the optimization cycle.

Accelerator-driven systems↗

NASA NTRS1976

Nonlinear periodic waves

Systematic perturbation procedures for the analysis of nonlinear problems are reviewed. The cases when the multiplicity of an eigenvalue is finite or infinite are treated for self-sustained and forced oscillations. The possibility of the formation of shock waves is discussed. Applications to acoustic problems are presented.

Ting, L.↗

NASA NTRS1979

Analysis of integral controls in linear quadratic regulator design

The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.

Slater, G. L.↗

NASA NTRS1984

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗

NASA NTRS1984

Perturbation analysis of optimal integral controls

The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.

Slater, G. L.↗

NASA NTRS1985

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗

NASA NTRS1992

Negative-energy states in the Dirac-Hartree-Fock problem - The effect of omission of two-electron integrals involving the small component

The effect of omission of two-electron integrals involving basis functions for the small component of the wavefunction on the eigenvalue spectrum in the Dirac-Hartree-Fock problem is studied. From an analysis of the Fock matrix it is shown that omission of these integrals moves the negative-energy states down, not up. Their complete omission does not give rise to intruder states. The appearance of intruder states occurs when only some of the core integrals are omitted, due to the nature of particular contraction schemes used for the core basis functions. Use of radially localized functions rather than atomic functions alleviates the intruder state problem.

Dyall, Kenneth G.↗

NASA NTRS1993

Sensitivity of discrete systems

Sensitivity techniques are examined for structures discretized by finite element analysis. The discussion covers finite difference sensitivities, sensitivity of static response, sensitivity of eigenvalues and eigenvectors, sensitivity of linear transient response, and an example of use of transient sensitivity.

Haftka, Raphael T.↗

NASA NTRS1993

Viscous analysis of three-dimensional rotor flows using a multigrid method

A three-dimensional code for rotating blade-row flow analysis was developed. The space discretization uses a cell-centered scheme with eigenvalues scaling for the artificial dissipation. The computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. An application is presented for the NASA rotor 67 transonic fan. Due to the blade stagger and twist, a zonal, non-periodic H-type grid is used to minimize the mesh skewness. The calculation is validated by comparing it with experiments in the range from the maximum flow rate to a near-stall condition. A detailed study of the flow structure near peak efficiency and near stall is presented by means of pressure distribution and particle traces inside boundary layers.

Arnone, A.↗

NASA NTRS2005

Sunspot Pattern Classification using PCA and Neural Networks (Poster)

The sunspot classification scheme presented in this paper is considered as a 2-D classification problem on archived datasets, and is not a real-time system. As a first step, it mirrors the Zuerich/McIntosh historical classification system and reproduces classification of sunspot patterns based on preprocessing and neural net training datasets. Ultimately, the project intends to move from more rudimentary schemes, to develop spatial-temporal-spectral classes derived by correlating spatial and temporal variations in various wavelengths to the brightness fluctuation spectrum of the sun in those wavelengths. Once the approach is generalized, then the focus will naturally move from a 2-D to an n-D classification, where "n" includes time and frequency. Here, the 2-D perspective refers both to the actual SOH0 Michelson Doppler Imager (MDI) images that are processed, but also refers to the fact that a 2-D matrix is created from each image during preprocessing. The 2-D matrix is the result of running Principal Component Analysis (PCA) over the selected dataset images, and the resulting matrices and their eigenvalues are the objects that are stored in a database, classified, and compared. These matrices are indexed according to the standard McIntosh classification scheme.

Rajkumar, T.↗

NASA NTRS1984

Sensitivity analysis of discrete structural systems: A survey

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗