Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eigenvalue problems”

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 595 records · Page 33

(U) LNK3DNT Capability in SENSMG

The SENSMG multigroup neutron sensitivity code has been extended to read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Not all of the capability of SENSMG is available for redoin/lnk3dnt inputs. The only responses allowed are k eff or α eigenvalues. Fixed-source problems are not allowed, and there is no mechanism to input reaction rates or reaction-rate ratios. The next section of this report describes the user interface. Section III describes the verification of the capability. Section IV is a summary and conclusions. The input files used for verification are listed in the appendix.

97 MATHEMATICS AND COMPUTING↗

Analytical and experimental study of the vibration of bonded beams with a lap joint

A theoretical model to study the flexural vibration of a bonded lap joint system is described in this paper. First, equations of motion at the joint region are derived using a differential element approach. The transverse displacements of the upper and lower beam are considered to be different. The adhesive is assumed to be linearly viscoelastic and the widely used Kelvin-Voight model is used to represent the viscoelastic behavior of the adhesive. The shear force at the interface between the adhesive and the beam is obtained from the simple bending motion equations of the two beams. The resulting equations of motion are combined with the equations of transverse vibration of the beams in the unjointed regions. These are later solved as a boundary value problem to obtain the eigenvalues and eigenvectors of the system. The model can be used to predict the natural frequencies, modal damping ratios, and mode shapes of the system for free vibration. Good agreement between numerical and experimental results was obtained for a system of graphite epoxy beams lap-jointed by an epoxy adhesive.

Rao, M. D.↗

A new direction in hydrodynamic stability: Beyond eigenvalues

Fluid flows that are smooth at low speeds become unstable and then turbulent at higher speeds. This phenomenon has traditionally been investigated by linearizing the equations of flow and looking for unstable eigenvalues of the linearized problem, but the results agree poorly in many cases with experiments. Nevertheless, it has become clear in recent years that linear effects play a central role in hydrodynamic instability. A reconciliation of these findings with the traditional analysis can be obtained by considering the 'pseudospectra' of the linearized problem, which reveals that small perturbations to the smooth flow in the form of streamwise vortices may be amplified by factors on the order of 10(exp 5) by a linear mechanism, even though all the eigenmodes are stable. The same principles apply also to other problems in the mathematical sciences that involve non-orthogonal eigenfunctions.

Trefethen, Lloyd N.↗

Frequency domain FIR and IIR adaptive filters

A discussion of the LMS adaptive filter relating to its convergence characteristics and the problems associated with disparate eigenvalues is presented. This is used to introduce the concept of proportional convergence. An approach is used to analyze the convergence characteristics of block frequency-domain adaptive filters. This leads to a development showing how the frequency-domain FIR adaptive filter is easily modified to provide proportional convergence. These ideas are extended to a block frequency-domain IIR adaptive filter and the idea of proportional convergence is applied. Experimental results illustrating proportional convergence in both FIR and IIR frequency-domain block adaptive filters is presented.

Lynn, D. W.↗

Application of a systematic finite-element model modification technique to dynamic analysis of structures

A systematic finite-element model modification technique has been applied to two small problems and a model of the main wing box of a research drone aircraft. The procedure determines the sensitivity of the eigenvalues and eigenvector components to specific structural changes, calculates the required changes and modifies the finite-element model. Good results were obtained where large stiffness modifications were required to satisfy large eigenvalue changes. Sensitivity matrix conditioning problems required the development of techniques to insure existence of a solution and accelerate its convergence. A method is proposed to assist the analyst in selecting stiffness parameters for modification.

Robinson, J. C.↗

Design of bearings for rotor systems based on stability

Design of rotor systems incorporating stable behavior is of great importance to manufacturers of high speed centrifugal machinery since destabilizing mechanisms (from bearings, seals, aerodynamic cross coupling, noncolocation effects from magnetic bearings, etc.) increase with machine efficiency and power density. A new method of designing bearing parameters (stiffness and damping coefficients or coefficients of the controller transfer function) is proposed, based on a numerical search in the parameter space. The feedback control law is based on a decentralized low order controller structure, and the various design requirements are specified as constraints in the specification and parameter spaces. An algorithm is proposed for solving the problem as a sequence of constrained 'minimax' problems, with more and more eigenvalues into an acceptable region in the complex plane. The algorithm uses the method of feasible directions to solve the nonlinear constrained minimization problem at each stage. This methodology emphasizes the designer's interaction with the algorithm to generate acceptable designs by relaxing various constraints and changing initial guesses interactively. A design oriented user interface is proposed to facilitate the interaction.

Dhar, D.↗

An automated procedure for computing flutter eigenvalues.

A new, fast and economical automated procedure for implementing the traditional V-g method of flutter solution is described. The procedure requires as input the generalized aerodynamic forces for a range of reduced frequencies obtained from an aerodynamic program. These aerodynamic forces are interpolated with respect to reduced frequency using a newly developed, partially tabulated cubic spline that is both fast in execution and economical in storage. The flutter solution is then obtained using an eigenvalue routine that has been developed to take advantage of the parametric nature of the V-g type of solution. Furthermore, the routine takes care of the fundamental and troublesome problem of properly sorting the output eigenvalues. By solving the root-sorting problem, the interpolation for flutter crossings and automatic plotting are accomplished efficiently. The computational techniques used in this new program are described and some sample results are given.

Desmarais, R. N.↗

Eigenvalues of singular differential operators by finite difference methods. I.

Approximation of the eigenvalues of certain self-adjoint operators defined by a formal differential operator in a Hilbert space. In general, two problems are studied. The first is the problem of defining a suitable Hilbert space operator that has eigenvalues. The second problem concerns the finite difference operators to be used.

Baxley, J. V.↗

An Eigenvalue Analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L(sub 2) stability on a finite domain.

Warming, Robert F.↗

An eigenvalue analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L (sub 2) stability on a finite domain.

Warming, Robert F.↗

Sturm-Liouville eigenproblems with an interior pole

The eigenvalues and eigenfunctions of self-adjoint Sturm-Liouville problems with a simple pole on the interior of an interval are investigated. Three general theorems are proved, and it is shown that as n approaches infinity, the eigenfunctions more and more closely resemble those of an ordinary Sturm-Liouville problem. The low-order modes differ significantly from those of a nonsingular eigenproblem in that both eigenvalues and eigenfunctions are complex, and the eigenvalues for all small n may cluster about a common value in contrast to the widely separated eigenvalues of the corresponding nonsingular problem. In addition, the WKB is shown to be accurate for all n, and all eigenvalues of a normal one-dimensional Sturm-Liouville equation with nonperiodic boundary conditions are well separated.

Boyd, J. P.↗

Impact of Radial Reflector Fidelity on Neutronics and Vessel Fluence Simulations

The Consortium for Advanced Simulation of Light Water Reactors is developing the Virtual Environment for Reactor Applications (VERA), and the MPACT code, which is the primary deterministic neutron transport solver in VERA, provides sub-pin level flux and power distributions as part of full-scale cycle depletion and analysis. In such calculations, an important aspect is the radial reflector treatment. To improve the fidelity of the radial reflector treatment, MPACT was extended to approximate the modeling of the reactor’s structural components such as the core shroud, barrel, neutron pads, and vessel. This work explores several modeling configurations with varying levels of fidelity and computational burden and assesses the importance of modeling fidelity on the eigenvalue and pin power distribution. Two two-dimensional (2-D) problems were analyzed to assess the impact on eigenvalue and pin power distributions with low-fidelity, coarse square cell reflector representations: (1) a Watts Bar Nuclear Plant Unit 1 (WBN1) quarter-core slice with depletion and (2) an AP1000 quarter-core slice. In this work, the analyses showed that the effect on eigenvalue is fairly small, but the effect on pin power is more pronounced, especially locally in the assemblies closest to the periphery, where the maximum pin power difference is nearly 3.5% in the AP1000 case. Two additional 2-D problems were used to assess the comparison between the low-fidelity coarse square cell treatment and a high-fidelity geometric representation that uses subpin material specification: (1) the same WBN1 quarter-core slice and (2) a representative model of the NuScale small modular reactor (SMR), which features a solid reflector design with moderator holes. These results demonstrate that even a coarse, low-fidelity representation adequately captures the necessary simulation characteristics. Last, these capabilities were applied to the 2-D WBN1 quarter-core depletion to assess the impact on vessel fluence using VeraShift. From adjoint calculations, pins along the periphery were observed to be of highest importance for fluence calculation, so the impact of the reflector representation in MPACT could theoretically substantially affect the predicted result. However, it was observed that the change in pin powers along the periphery minimally impacts the maximum vessel fluence with a difference within the statistical uncertainty but provides terrific insight on the sensitivity of the peripheral pins.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Relativistic levels of mesic atoms

We revisit the derivation of the covariant two-body scalar-fermion equation with a Coulomb interaction, presented in a previous paper. We show that it can be given the formal aspect of a Dirac equation, but for the fact that the eigenvalue is also contained in one of the coefficients and thus it is not linearly included. The discussion of the boundary value problem is therefore different, although some properties of the Dirac equation can be recovered in an approximation bringing back to the concept of reduced mass. We discuss a mixed analytic-numerical method of solution which allows to obtain very accurate results and we calculate the lowest levels, states and QED corrections for pionic and kaonic atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sampling electronic structure quadratic unconstrained binary optimization problems (QUBOs) with Ocean and Mukai solvers

The most advanced D-Wave Advantage quantum annealer has 5000+ qubits, however, every qubit is connected to a small number of neighbors. As such, implementation of a fully-connected graph results in an order of magnitude reduction in qubit count. To compensate for the reduced number of qubits, one has to rely on special heuristic software such as qbsolv, the purpose of which is to decompose a large quadratic unconstrained binary optimization (QUBO) problem into smaller pieces that fit onto a quantum annealer. In this work, we compare the performance of the open-source qbsolv which is a part of the D-Wave Ocean tools and a new Mukai QUBO solver from Quantum Computing Inc. (QCI). The comparison is done for solving the electronic structure problem and is implemented in a classical mode (Tabu search techniques). The Quantum Annealer Eigensolver is used to map the electronic structure eigenvalue-eigenvector equation to a QUBO problem, solvable on a D-Wave annealer. We find that the Mukai QUBO solver outperforms the Ocean qbsolv with one to two orders of magnitude more accurate energies for all calculations done in the present work, both the ground and excited state calculations. This work stimulates the further development of software to assist in the utilization of modern quantum annealers.

97 MATHEMATICS AND COMPUTING↗

On eigenvectors of multiple eigenvalues obtained in NASTRAN

In the case of nonmultiple eigenvalues, each of the three real eigenvalue extraction methods available in NASTRAN will, for a given type of normalization, give essentially the same eigenvectors, but this is not so in the case of multiple eigenvalues. This discrepancy is explained and illustrated by considering the example of a NASTRAN demonstration problem that has both multiple and nonmultiple eigenvalues.

Pamidi, P. R.↗

Verification of the pin power reconstruction method in NESTLE

NESTLE is a 3D nodal core simulator that solves the few group neutron diffusion equation using the nodal expansion method (NEM) for the steady state (eigenvalue or fixed-source), adjoint (eigenvalue), and transient (eigenvalue or fixed-source) problems. NESTLE has been in use for over 25 years and has been continuously improved and developed. NESTLE can be used to model either Cartesian geometries (for use with typical PWR, BWR fuel assemblies) or hexagonal geometries (for use with VVER and other triangularly pitched fuel assemblies). As part of NESTLE's ongoing verification and validation process, several computational comparisons have been made to quantify the discrepancies in NESTLE's hexagonal pin power reconstruction method. In this paper, we present the results of these code-to-code comparisons and find that the root mean squared (RMS) error in pin power is comparable to other nodal codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical pole assignment by eigenvalue Jacobian inversion

A numerical procedure for solving the linear pole placement problem is developed which operates by the inversion of an analytically determined eigenvalue Jacobian matrix. Attention is given to convergence characteristics and pathological situations. It is not concluded that the algorithm developed is suitable for computer-aided control system design with particular reference to the scan platform pointing control system for the Galileo spacecraft.

Sevaston, George E.↗

Stability and bifurcation of quasiparallel Alfven solitons

The inverse scattering transformation (IST) is used to study the one-parameter and two-parameter soliton families of the derivative nonlinear Schroedinger (DNLS) equation. The two-parameter soliton family is determined by the discrete complex eigenvalue spectrum of the Kaup-Newell scattering problem and the one-parameter soliton family corresponds to the discrete real eigenvalue spectrum. The structure of the IST is exploited to discuss the existence of discrete real eigenvalues and to prove their structural stability to perturbations of the initial conditions. Also, though the two-parameter soliton is structurally stable in general, it is shown that a perturbation of the initial conditions may change the two-parameter soliton into a degenerate soliton which, in turn, is structurally unstable. This degenerate, or double pole, soliton may bifurcate due to a perturbation of the initial conditions into a pair of one-parameter solitons. If the initial profile is on compact support, then this pair of one-parameter solitons must be compressive and rarefactive respectively. Finally, the Gelfand-Levitan equations appropriate for the double pole soliton are solved.

Hamilton, R. L.↗