Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 487 records · Page 27

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗

Numerical Modeling of Nanoelectronic Devices

Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

Klimeck, Gerhard↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

Finite Element Analysis of Reverberation Chambers

The primary motivating factor behind the initiation of this work was to provide a deterministic means of establishing the validity of the statistical methods that are recommended for the determination of fields that interact in -an avionics system. The application of finite element analysis to reverberation chambers is the initial step required to establish a reasonable course of inquiry in this particularly data-intensive study. The use of computational electromagnetics provides a high degree of control of the "experimental" parameters that can be utilized in a simulation of reverberating structures. As the work evolved there were four primary focus areas they are: 1. The eigenvalue problem for the source free problem. 2. The development of a complex efficient eigensolver. 3. The application of a source for the TE and TM fields for statistical characterization. 4. The examination of shielding effectiveness in a reverberating environment. One early purpose of this work was to establish the utility of finite element techniques in the development of an extended low frequency statistical model for reverberation phenomena. By employing finite element techniques, structures of arbitrary complexity can be analyzed due to the use of triangular shape functions in the spatial discretization. The effects of both frequency stirring and mechanical stirring are presented. It is suggested that for the low frequency operation the typical tuner size is inadequate to provide a sufficiently random field and that frequency stirring should be used. The results of the finite element analysis of the reverberation chamber illustrate io-W the potential utility of a 2D representation for enhancing the basic statistical characteristics of the chamber when operating in a low frequency regime. The basic field statistics are verified for frequency stirring over a wide range of frequencies. Mechanical stirring is shown to provide an effective frequency deviation.

Bunting, Charles F.↗

Dynamic analysis of the joint dominated beam

A method for determining the vibrational modes of the joint dominated beam is presented with reference to space erectable structures. Mathematical models with and without joint imperfection are developed. Joint imperfection is represented by a simple flexibility, k, in the axial direction of the chord members only. An interactive computer program is used to obtain eigenvalue characteristics and mode shapes of the system. Results are presented on frequencies of the given truss beam with and without joint imperfections for three-panel and 10-panel models. Also presented are results of parametric studies of the beam truss models of three and 10 panels with k = 0.05, 0.10, 0.50, and 1.0 x 10 to the -7th.

Abu-Saba, Elias G.↗

Adaptive control of Space Station during nominal operations with CMGs

An adaptive control approach is investigated for the Space Station. The main components of the adaptive controller are the parameter identification scheme, the control gain calculation, and the control law. The control law is the Space Station baseline control law. The control gain calculation is based on linear quadratic regulator theory with eigenvalue placement in a vertical strip. The parameter identification scheme is a real-time recursive extended Kalman filter which 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 compute accurate estimates of the Space Station inertias during nominal CMG (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, R. H.↗

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗

Viscous and inviscid instabilities of a trailing vortex

The linear stability of the trailing line vortex model of Batchelor (1964) is studied using a spectral collocation and matrix eigenvalue method. The entire unstable region in the swirl/axial wavenumber parameter space is mapped out for various azimuthal wavenumbers for both the inviscid and viscous stability problem. The results of the study provide a direct numerical validation of the large-azimuthal-wavenumber asymptotic analysis of Leibovich and Stewartson (1983). It is shown that accurate results are obtained up to azimuthal wavenumbers of 10,000 and greater, and the agreement with the asymptotic theory is excellent.

Mayer, Ernst W.↗

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

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Linear beam stability in periodic focusing systems: Krein signature and band structure

The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characteristics of one-period maps, and discussed the stability properties of single-particle motions with skew quadrupoles and envelope perturbations in high-intensity beams on an equal footing. In particular, an application of the Krein theory to envelope instability analysis was newly attempted in this study. Finally, the appearance of instabilities is interpreted as the result of the collision between eigenmodes of opposite Krein signatures and the formation of a band gap.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effect of altered sensory conditions on multivariate descriptors of human postural sway

Multivariate descriptors of sway were used to test whether altered sensory conditions result not only in changes in amount of sway but also in postural coordination. Eigenvalues and directions of eigenvectors of the covariance of shnk and hip angles were used as a set of multivariate descriptors. These quantities were measured in 14 healthy adult subjects performing the Sensory Organization test, which disrupts visual and somatosensory information used for spatial orientation. Multivariate analysis of variance and discriminant analysis showed that resulting sway changes were at least bivariate in character, with visual and somatosensory conditions producing distinct changes in postural coordination. The most significant changes were found when somatosensory information was disrupted by sway-referencing of the support surface (P = 3.2 x 10(-10)). The resulting covariance measurements showed that subjects not only swayed more but also used increased hip motion analogous to the hip strategy. Disruption of vision, by either closing the eyes or sway-referencing the visual surround, also resulted in altered sway (P = 1.7 x 10(-10)), with proportionately more motion of the center of mass than with platform sway-referencing. As shown by discriminant analysis, an optimal univariate measure could explain at most 90% of the behavior due to altered sensory conditions. The remaining 10%, while smaller, are highly significant changes in posture control that depend on sensory conditions. The results imply that normal postural coordination of the trunk and legs requires both somatosensory and visual information and that each sensory modality makes a unique contribution to posture control. Descending postural commands are multivariate in nature, and the motion at each joint is affected uniquely by input from multiple sensors.

NASA Discipline Neuroscience↗

Instability analysis procedure for 3-level multi-bearing rotor-foundation systems

A procedure for the instability analysis of a three-level multispan rotor systems is described. This procedure is based on a distributed mass elastic representation of the rotor system in several eight-coefficient bearings. Each bearing is supported from an elastic foundation on damped, elastic pedestals. The foundation is represented as a general distributed mass elastic structure on discrete supports, which may have different stiffness and damping properties in the horizontal and vertical directions. This system model is suited to studies of instability threshold conditions for multirotor turbomachines on either massive or flexible foundations. The instability conditions is found by obtaining the eigenvalues of the system determinant, which is obtained by the transfer matrix method from the three-level system model. The stability determinant is solved for the lowest rotational speed at which the system damping becomes zero in the complex eigenvalue, and for the whirl frequency corresponding to the natural frequency of the unstable mode. An efficient algorithm for achieving this is described. Application of this procedure to a rigid rotor in two damped-elastic bearings and flexible supports is described. A second example discusses a flexible rotor with four damped-elastic bearings. The third case compares the stability of a six-bearing 300 Mw turbine generator unit, using two different bearing types. These applications validate the computer program and various aspects of the analysis.

Zhou, S.↗

Structural representation for analysis of a controlled structure

The purpose was to determine what reduced order structural representation is most appropriate for coupling with a control system. The goal was to choose a reduced order structural model which retains as closely as possible the characteristics of the closed loop model with a full order structural representation. By characteristics of the closed loop model, it is meant that the closed loop eigenvalues and the closed loop transfer functions from commands to loads and from commands to response. This process does not address the accuracy of the full order model (usually a finite element model) but only the loss of accuracy associated with reducing th model. For the purposes of this study, only collocated sensors and actuators are examined. The choice of a structural representation for noncollocated sensors and actuators is not so clear.

Blelloch, Paul A.↗

Flap-lag-torsion flutter analysis of a constant life rotor

The constant lift rotor (CLR) employs a control input of pitch moment to several airfoil sections which are free to pivot on a continuous spar, allowing them to change their pitch to obtain the desired lift. A flap-lag-torsion flutter analysis of a constant lift rotor blade in hover was developed. The blade model assumes rigid body flap and lead-lag motions at the root hinge and each strip undergoes an independent torsional motion. The results are presented in terms of root locus plots of complex eigenvalues as a function of thrust. The effects of several parameters (including structural damping, center of gravity and elastic axis offset from aerodynamic center, compressibility pitch-lag and pitch-flap coupling) on the blade dynamics are examined. With a suitable combination of lag damper and pitch-flap coupling, it is possible to design a constant lift rotor blade free from flutter instability.

Chopra, I.↗

SCALE Modeling of the Fast Spectrum Heat Pipe Reactor

As part of the severe accident analysis collaboration with Sandia National Laboratories (SNL) and the Nuclear Regulatory Commission (NRC), SCALE models were developed for a fast-spectrum heat pipe reactor. These models were based on the Idaho National Laboratory (INL) Design A concept, which is an alternative design to the Los Alamos National Laboratory (LANL) Special Purpose Reactor (SPR), also known as the Megapower reactor. The model contains 1,134 heat pipes, surrounded by hexagonal fuel elements, with a potassium working fluid; the fuel is UO 2 with 19.75 wt% 235 U enrichment. The model contains axial beryllium oxide (BeO) reflectors above and below the active fuel region along with a radial alumina reflector containing 12 B 4 C control drums. The center of the core is left unfueled to make room for two shutdown control rods, one annular and one solid. The active region of the core was discretized into twenty axial and five radial zones to analyze spatial variations in power and burnup. Infinite lattice unit cell sensitivity studies were used to perform verification between the SCALE and INL models. The eigenvalue results agreed well with the reported results to within roughly 50 percent mille (pcm). Full-core model verification was performed by analyzing system eigenvalues with differing configurations of control drum and shutdown rod positions. These full core results all had eigenvalue differences less than 310 pcm. Control drum and shutdown rod worths were also compared, with differences of 3.2% or less. Using the verified model, the isotopic inventory and decay heat, as well as temperature feedback coefficients, were calculated and provided to SNL as input to the MELCOR severe accident code to analyze potential releases from this class of reactor. The results of the MELCOR analysis are provided in a different report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Solving the homogeneous Bethe-Salpeter equation with a quantum annealer

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other nonsymmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the GEVP by means of a quantum annealer, i.e., to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the simulated-annealing package and the D-Wave . The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗