Search NASA⌕ Search

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Structural stability and artificial buckling modes in topology optimization

Abstract This paper demonstrates how a strain energy transition approach can be used to remove artificial buckling modes that often occur in stability constrained topology optimization problems. To simulate the structural response, a nonlinear large deformation hyperelastic simulation is performed, wherein the fundamental load path is traversed using Newton’s method and the critical buckling load levels are estimated by an eigenvalue analysis. The goal of the optimization is to minimize displacement, subject to constraints on the lowest critical buckling loads and maximum volume. The topology optimization problem is regularized via the Helmholtz PDE-filter and the method of moving asymptotes is used to update the design. The stability and sensitivity analyses are outlined in detail. The effectiveness of the energy transition scheme is demonstrated in numerical examples.

Dalklint, Anna (ORCID:0000000346195205)↗
DOE OSTI2016

Fast Multilevel Solvers for a Class of Discrete Fourth Order Parabolic Problems

In this paper, we study fast iterative solvers for the solution of fourth order parabolic equations discretized by mixed finite element methods. We propose to use consistent mass matrix in the discretization and use lumped mass matrix to construct efficient preconditioners. We provide eigenvalue analysis for the preconditioned system and estimate the convergence rate of the preconditioned GMRes method. Furthermore, we show that these preconditioners only need to be solved inexactly by optimal multigrid algorithms. Our numerical examples indicate that the proposed preconditioners are very efficient and robust with respect to both discretization parameters and diffusion coefficients. We also investigate the performance of multigrid algorithms with either collective smoothers or distributive smoothers when solving the preconditioner systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Synergistic data analytics for electromechanical oscillation in electric power systems

Accurate real-time estimation of the four electromechanical oscillation properties, i.e., dominant oscillation modes, mode shapes, participation factors, and coherent groups, is of great importance to assess and mitigate potential electromechanical oscillations in interconnected power systems. While eigenvalue analysis can realize such estimation, it requires precise linearized dynamic models and accurate parameters, which are highly difficult to obtain in practice. Data fusion-based modal estimation methods can extract the properties of electromechanical oscillations from measurement data without the power system model and parameters, but most of the time only one or two property assessments can be accomplished each time. To overcome this challenge, this paper presents a synergistic data analytics solution to characterize the dynamic behaviors of electromechanical oscillations from real-time measurement data. Here, the proposed method uses optimized variable projection, and it is capable of estimating all four electromechanical oscillation behavior properties from measured responses. Case studies are performed using the simulated measurement data of a 16-generator 68-bus test system and the field measurements collected by the PMUs deployed in the Yunnan Power Grid. The results demonstrate that the proposed synergistic data analytics solution can achieve satisfactory performance in capturing the properties of electromechanical oscillations from measurement data and exhibit strong robustness against measurement noise when compared with existing measurement-based methods.

42 ENGINEERING↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗
DOE OSTI2018

New aspects of longitudinal instabilities in electron storage rings

Here, novel features of the longitudinal instability of a single electron bunch circulating in a low-emittance electron storage ring are discussed. Measurements and numerical simulations, performed both in time and frequency domain, show a non-monotonic increase of the electron beam energy spread as a function of single bunch current, characterized by the presence of local minima and maxima, where a local minimum of the energy spread is interpreted as a higher-order microwave instability threshold. It is also shown that thresholds related to the same zero-intensity bunch length depend linearly on the accelerating radio frequency voltage. The observed intensity-dependent features of the energy spread, confirmed by measurements with two independent diagnostics methods, i.e. horizontal beam profile measurements by a synchrotron light monitor and photon energy spectrum measurements of undulator radiation, are given a theoretical interpretation by applying a novel eigenvalue analysis based on the linearized Vlasov equation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Convective mixing in vertically-layered porous media: The linear regime and the onset of convection

In this paper, we study the effect of permeability heterogeneity on the stability of gravitationally unstable, transient, diffusive boundary layers in porous media. Permeability is taken to vary periodically in the horizontal plane normal to the direction of gravity. In contrast to the situation for vertical permeability variation, the horizontal perturbation structures are multimodal. We therefore use a two-dimensional quasi-steady eigenvalue analysis as well as a complementary initial value problem to investigate the stability behavior in the linear regime, until the onset of convection. We find that thick permeability layers enhance instability compared with thin layers when heterogeneity is increased. On the contrary, for thin layers the instability is weakened progressively with increasing heterogeneity to the extent that the corresponding homogeneous case is more unstable. For high levels of heterogeneity, we find that a small change in the permeability field results in large variations in the onset time of convection, similar to the instability event in the linear regime. However, this trend does not persist unconditionally because of the reorientation of vorticity pairs due to the interaction of evolving perturbation structures with heterogeneity. Consequently, an earlier onset of instability does not necessarily imply an earlier onset of convection. A resonant amplification of instability is observed within the linear regime when the dominant perturbation mode is equal to half the wavenumber of permeability variation. On the other hand, a substantial damping occurs when the perturbation mode is equal to the harmonic and sub-harmonic components of the permeability wavenumber. Finally, the phenomenon of such harmonic interactions influences both the onset of instability as well as the onset of convection.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical study of transition between even and odd toroidal Alfvén eigenmodes on EAST

In this work, linear simulations of toroidal Alfvén eigenmodes (TAEs) driven by energetic particles (EPs) on EAST (Experimental Advanced Superconducting Tokamak) are performed using the hybrid-kinetic MHD model implemented in NIMROD code. Furthermore, the EAST equilibrium is reconstructed using the EFIT code based on the experimental measurement. The “slowing down” distribution is adopted for modeling the equilibrium distribution of the energetic ions from the deuterium neutral beam injection on EAST. The frequency, the dominant poloidal mode number, the radial location and the detailed 2D mode structure of the TAE/RSAE (Reverse Shear Alfvén Eigenmode)/EPM (energetic particle mode) modes are consistent between the eigenvalue analysis and the NIMROD simulation. As the β fraction of EP increases, a transition between even and odd TAEs occurs, along with that between the ballooning and antiballooning mode structures. It is demonstrated that when the β fraction of EP is close to the transition threshold, both types of TAEs coexist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2022

AlgoSCR: an algorithm for solar contamination removal from radio interferometric data

Hydrogen intensity mapping is a new field in astronomy that promises to make three-dimensional maps of the matter distribution of the Universe using the redshifted 21cm line of neutral hydrogen gas (HI). Several ongoing and upcoming radio interferometers, such as Tianlai, CHIME, HERA, HIRAX, etc., are using this technique. These instruments are designed to map large swaths of the sky by drift scanning over periods of many months. One of the challenges of the observations is that the daytime data are contaminated by strong radio signals from the Sun. In the case of Tianlai, this results in almost half of the measured data being unusable. We try to address this issue by developing an algorithm for solar contamination removal (AlgoSCR) from the radio data. The algorithm is based on an eigenvalue analysis of the visibility matrix and hence is applicable only to interferometers. We apply AlgoSCR to simulated visibilities, as well as real daytime data from the Tianlai dish array. The algorithm can reduce strong solar contamination by about 95 per cent without seriously affecting other weaker sky signals and thus makes the data usable for certain applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Impact of Frequency-Watt Control on the Dynamics of a High DER Penetration Power System

Inverter-coupled generation is slowly displacing synchronous generators in the power system, resulting in reduced system inertia and primary frequency response, which makes the system susceptible to larger frequency deviations. Such deviations may result in load shedding or tripping of inverters, thus reducing the reliability. Frequency-watt control of inverters is an autonomous approach to support grid frequency at a fast timescale by adjusting the inverter output with changes in frequency. This study investigates the dynamics of a power system with a high penetration of distributed PV generation with frequency-watt droop control enabled. A simple aggregated reduced-order two-source model of the O'ahu Island (Hawaii) power system is developed and validated against a full dynamic model, and the small-signal stability of the reduced model is then examined the under various PV penetrations. The impact of frequency-watt parameters on the system frequency dynamics is investigated by observing trajectories of small- signal eigenvalues. The results of the eigenvalue analysis are validated through time-domain simulation. Considerations for frequency-watt parameter selection are discussed.

14 SOLAR ENERGY↗
DOE OSTI2020

Isochronous Architecture-Based Voltage-Active Power Droop for Multi-Inverter Systems

This article proposes an isochronous architecture for parallel inverters with only voltage-active power droop (VP-D) control for improving active power sharing as well as plug-and-play of multi-inverter-based distributed energy resources (DERs). The isochronous framework obviates the need for an explicit regulation of the frequency while allowing for sharing of reactive power. The article shows that the detrimental effects of circulating currents between inverters can be addressed in the framework developed. The isochronous architecture is implemented by employing a GPS to disseminate clock timing signals that enable the microgrid to maintain the nominal system frequency. Small signal eigenvalue analysis of a Multi-inverter Microgrid system near the steady-state operating point is presented to evaluate the system stability. Moreover, unlike traditional droop-based methods, the isochronous strategy lends itself to analytical guarantees which are developed in the article. Furthermore, the effect of delays in receiving GPS clock signals by an inverter is experimentally characterized and an inverter isolation mechanism is proposed that disconnects inverters facing large GPS communication delays. Hardware experiments on an 1.2 kVA-prototype and controller-hardware-in-the-loop validation on the CIGRE distribution network are conducted to demonstrate the effectiveness of the proposed architecture towards active and reactive power sharing between inverters with load scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1973

A simplified modal plotting technique for the representation of complex structural models

A plotting method has been developed that has proven to be useful in rapidly defining the mode shapes of a complex spacecraft. The reduction of the complex model to a simple stick plot is accomplished by augmenting the existing NASTRAN plot package with constrained plot elements. Both the existing NASTRAN modal plot package and the stick plots of a representative mode are presented for comparison. The NASTRAN plot package provides the capability of generating modal deformations (mode shapes) resulting from real eigenvalue analysis. These modal deformations of the structural model may be displayed in the deformed shape either alone or superimposed on the undeformed shape. Another available method of displaying the modal deformations is by means of displacement vectors at the grid points

Hanlein, S. L.↗
NASA NTRS1974

Attitude stability of a flexible solar electric spacecraft - A parametric study

The influence of large area solar array flexibility on the attitude control of a solar electric spacecraft is analyzed. The spacecraft consists of a rigid central body, two flexible roll-out solar arrays, and a cluster of three electric propulsion engines. Each engine has gimballing or translational capacity, or both, for attitude control purposes. A parametric stability study of various spacecraft configurations was made, using hybrid coordinate techniques in an eigenvalue analysis. Parameters varied were solar array aspect ratio, solar array first natural frequency, solar array rotation angle about the yaw axis, and sensor gain factors. Only the first six modes of solar array vibration were included in the study. The attitude response characteristics of the various spacecraft configurations were obtained by means of a digital computer simulation program of the system's equations of motion. The feasibility of using the thrusters for attitude control of electric propulsion spacecraft is demonstrated. It is also shown that although the dynamics effects of large-area flexible solar arrays can be harmful, changes in the control system parameters can insure stable attitude control.

Marsh, E. L.↗
NASA NTRS1978

Tidal frequency estimation for closed basins

A method was developed for determining the fundamental tidal frequencies for closed basins of water, by means of an eigenvalue analysis. The mathematical model employed, was the Laplace tidal equations.

Eades, J. B., Jr.↗
NASA NTRS1983

Large space structure damping design

Several FORTRAN subroutines and programs were developed which compute complex eigenvalues of a damped system using different approaches, and which rescale mode shapes to unit generalized mass and make rigid bodies orthogonal to each other. An analytical proof of a Minimum Constrained Frequency Criterion (MCFC) for a single damper is presented. A method to minimize the effect of control spill-over for large space structures is proposed. The characteristic equation of an undamped system with a generalized control law is derived using reanalysis theory. This equation can be implemented in computer programs for efficient eigenvalue analysis or control quasi synthesis. Methods to control vibrations in large space structure are reviewed and analyzed. The resulting prototype, using electromagnetic actuator, is described.

Pilkey, W. D.↗
NASA NTRS1983

Large amplitude MHD waves upstream of the Jovian bow shock

Observations of large amplitude magnetohydrodynamics (MHD) waves upstream of Jupiter's bow shock are analyzed. The waves are found to be right circularly polarized in the solar wind frame which suggests that they are propagating in the fast magnetosonic mode. A complete spectral and minimum variance eigenvalue analysis of the data was performed. The power spectrum of the magnetic fluctuations contains several peaks. The fluctuations at 2.3 mHz have a direction of minimum variance along the direction of the average magnetic field. The direction of minimum variance of these fluctuations lies at approximately 40 deg. to the magnetic field and is parallel to the radial direction. We argue that these fluctuations are waves excited by protons reflected off the Jovian bow shock. The inferred speed of the reflected protons is about two times the solar wind speed in the plasma rest frame. A linear instability analysis is presented which suggests an explanation for many of the observed features of the observations.

Goldstein, M. L.↗