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 469 records · Page 26

Dynamic and thermal response finite element models of multi-body space structural configurations

Presented is structural dynamics modeling of two multibody space structural configurations. The first configuration is a generic space station model of a cylindrical habitation module, two solar array panels, radiator panel, and central connecting tube. The second is a 15-m hoop-column antenna. Discussed is the special joint elimination sequence used for these large finite element models, so that eigenvalues could be extracted. The generic space station model aided test configuration design and analysis/test data correlation. The model consisted of six finite element models, one of each substructure and one of all substructures as a system. Static analysis and tests at the substructure level fine-tuned the finite element models. The 15-m hoop-column antenna is a truss column and structural ring interconnected with tension stabilizing cables. To the cables, pretensioned mesh membrane elements were attached to form four parabolic shaped antennae, one per quadrant. Imposing thermal preloads in the cables and mesh elements produced pretension in the finite element model. Thermal preload variation in the 96 control cables was adjusted to maintain antenna shape within the required tolerance and to give pointing accuracy.

Edighoffer, Harold H.↗

Simultaneous iterations algorithm for general eigenvalue problems on parallel processors

The method of simultaneous iteration with shift is extended to extraction of m-eigenpairs of a general eigenvalue problem of large order n in a parallel processing environment. The algorithm combines the power method and the Jacobi technique, and reduces to performing four basic operations. Parallel implementation of the algorithm is discussed in detail. The analysis accounts for computation and communication costs, and utilizes a parallel processing architecture of the ensemble type. Expressions for the computational efficiency and speedup are defined as a function of the problem and hardware parameters. Selected representative problems exhibit efficiencies ranging from 60 to 98 percent.

Utku, S.↗

Digital controller design: Analysis of the annular suspension pointing system

The analog controllers of the annular suspension pointing system are designed for control of the chi, phi sub 1, and phi sub 2 bandwidth dynamics through decoupling and pole placement. Since it is virtually impossible to find an equivalent bandwidth of the overall system and establish a general eigenvalue requirement for the system, the subsystem dynamics are decoupled through state feedback and the poles are placed simultaneously to realize the desired bandwidths for the three system components. Decoupling and pole placement are also used to design the closed-loop digital system through approximation.

Kuo, B. C.↗

Quasi-deuteron model at low renormalization group resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photoabsorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteronlike quantum numbers. In Phys. Rev. C 104, 034311 (2021), we showed that SRC physics can be naturally accounted for at low RG resolution. We describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photoabsorption to that for photodisintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the polarization state of hydromagnetic fluctuations in the solar wind

From presently available observations one can infer that the Alfvenic turbulence measured in the solar wind, predominantly on trailing edges of high-speed streams, is a mixture of modes with two different polarizations, namely, Alfvenic modes and modes which are the incompressible limit of slow magnetosonic waves. Using Helios 2 magnetic data and a variance analysis, parallel (to the mean field) and perpendicular components of the fluctuations are separated, and the possible correlation between such components which would be predicted as a consequence of the incompressible character of the turbulence is studied. Correlations between eigenvalues of the variance matrix are also investigated and discussed.

Bavassano, B.↗

Robust Predictive Control for Modular Solid-State Transformer With Reduced DC Link and Parameter Mismatch

This paper presents the analysis and implementation of a predictive control method for dc-link regulation and voltage balance in a cascaded modular reduced dc-link solid-state transformer (SST). Passive components like bulky dc links limit the power density of power converters, especially medium-voltage (MV) SST. Reduced dc-link or low-inertia converters can dramatically reduce the size, cost, and weight by tolerating larger dc-link ripples and improve the reliability with electrolytic capacitor-less dc link. However, a small dc link leads to tight coupling between the input and the output stages, which is a challenge for control design. In stacked low-inertia converters (SLIC), the low-inertia converter modules are stacked for MV applications, resulting in coupling between the modules and making the control more challenging. A new model predictive control method which can achieve deadbeat regulation on the dc link without weighting factors has been proposed to address this novel problem. This paper focuses on analyzing the condition of the low-inertia dc link up to 80% ripple, the robustness of the control under parameter mismatches, high-order terms, and important implementation issues such as model-based sampling and computation delay compensation. Significantly, the high-order terms are introduced because of the large dc-link ripple. These high-order terms are unique to the SLIC and negligible in conventional high-inertia converters. A discrete-time large-signal model is built to capture the dc-link’s nonlinear dynamics, and the eigenvalues of a small-signal Jacobian matrix are analyzed with Floquet theory to evaluate stability, using the modular soft-switching solid-state transformer (M-S4T) as an example of the SLIC. Simulation and experimental results of an MVDC M-S4T verify the analysis and the predictive control method. Finally, the general application of the predictive control to low-inertia converters is compared against a conventional PI controller using a reduced dc-link active-front-end (AFE) rectifier as an example.

14 SOLAR ENERGY↗

Aeroelastic analysis of a troposkien-type wind turbine blade

The linear aeroelastic equations for one curved blade of a vertical axis wind turbine in state vector form are presented. The method is based on a simple integrating matrix scheme together with the transfer matrix idea. The method is proposed as a convenient way of solving the associated eigenvalue problem for general support conditions.

Nitzsche, F.↗

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

Wentzel-Kramers-Brillouin method in the Bargmann representation

It is demonstrated that the Bargmann representation of quantum mechanics is ideally suited for semiclassical analysis, using as an example the WKB method applied to the bound-state problem in a single well of one degree of freedom. For the harmonic oscillator, this WKB method trivially gives the exact eigenfunctions in addition to the exact eigenvalues. For an anharmonic well, a self-consistent variational choice of the representation greatly improves the accuracy of the semiclassical ground state. Also, a simple change of scale illuminates the relationship of semiclassical versus linear perturbative expansions, allowing a variety of multidimensional extensions.

Voros, A.↗

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

A study of delamination buckling of laminates

The subject of this paper is the buckling of laminated plates, with a preexisting delamination, subjected to in-plane loading. Each laminate is modelled as an orthotropic Mindlin plate. The analysis is carried out by a combination of the finite element and asymptotic expansion methods. By applying the finite element method, plates with general delamination regions can be studied. The asymptotic expansion method reduces the number of unknown variables of the eigenvalue equation to that of the equation for a single Kirchhoff plate. Numerical results are presented for several examples. The effects of the shape, size, and position of the delamination on the buckling load are studied through these examples.

Mukherjee, Yu-Xie↗

Adjoint-based Sensitivities of Flutter Predictions based on the Linearized Frequency-domain Approach

Flutter is a critical factor in designing and certifying aircraft. The linearized frequency-domain method offers a lower cost alternative to time-marching computational fluid dynamics for high-fidelity flutter analysis. In this work, adjoint-based sensitivities are added to a flutter analysis based on the linearized frequency-domain method to efficiently compute derivatives of flutter cost functions with respect to design variables or uncertain parameters. The derivation of the adjoint equations, which involve complications such as derivatives of a nonlinear generalized eigenvalue problem with complex-valued inputs and derivatives of the linearized Navier-Stokes equations, is provided. The implemented adjoint terms and derivatives are verified before demonstrating the approach for derivatives of flutter dynamic pressure with respect to Mach number for the AGARD 445.6 wing.

Aeroelasticity↗

Buckling analysis of a missile interstage

The NASTRAN program was used to find the buckling location and the degradation factor due to a specified temperature field for an interstage of a typical three stage missile. The model was made up of 6332 independent degrees of freedom. This large size required the use of an IBM 370/195 for evaluation of the eigenvalues. The effect of elevated temperature was to lower the critical buckling load by 3%.

Dudley, D. P.↗

Experiments at Sandia to Measure the Effect of Temperature on Critical Systems

Estimation of the uncertainty in a critical experiment attributable to uncertainties in the measured experiment temperature is done by calculating the variation of the eigenvalue of a benchmark configuration as a function of temperature. In the low-enriched water-moderated critical experiments performed at Sandia, this is done by 1) estimating the effects of changing the water temperature while holding the UO 2 fuel temperature constant, 2) estimating the effects of changing the UO 2 temperature while holding the water temperature constant, and 3) combining the two results. This assumes that the two effects are separable. The results of such an analysis are nonintuitive and need experimental verification. Critical experiments are being planned at Sandia National Laboratories (Sandia) to measure the effect of temperature on critical systems and will serve to test the methods used in estimating the temperature effects in critical experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experiments at Sandia to Measure the Effect of Temperature on Critical Systems [Slides]

Estimation of the uncertainty in a critical experiment attributable to uncertainties in the measured experiment temperature is done by calculating the variation of the eigenvalue of a benchmark configuration as a function of temperature. In the low-enriched water-moderated critical experiments performed at Sandia, this is done by 1) estimating the effects of changing the water temperature while holding the UO 2 fuel temperature constant, 2) estimating the effects of changing the UO 2 temperature while holding the water temperature constant, and 3) combining the two results. This assumes that the two effects are separable. The results of such an analysis are nonintuitive and need experimental verification. Critical experiments are being planned at Sandia National Laboratories (Sandia) to measure the effect of temperature on critical systems and will serve to test the methods used in estimating the temperature effects in critical experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗