Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Aeroelastic Analysis of Highly Flexible Wings with Linearized Frequency-Domain Aerodynamics

Flutter analysis of configurations with geometric structural nonlinearities typically is done with time-domain analysis. The results from this process are computationally expensive and can yield cumbersome results that may be difficult to manage/interpret. Compared to time-domain methods, frequency-domain flutter analysis can provide additional insight into the characteristics of a flutter stability problem. By linearizing the aeroelastic problem about the nonlinear equilibrium state, this work applies frequency-domain aeroelastic analysis to the Pazy wing, the subject of the Large Deformation Working Group in the Aeroelastic Prediction Workshop. Generalized aerodynamic forces (GAFs) are computed with both a doublet-lattice method and a computational fluid dynamics solver at a range of reduced frequencies as well as a range of dynamic pressures to account for the dependence of the mode shapes on the nonlinear equilibrium state. These GAFs are used in a p-k flutter solver, which is modified to handle the dependence of the stiffness matrix and GAFs on the dynamic pressure.

Bret K Stanford↗
NASA NTRS1987

Combustion stability analysis of preburners in liquid propellant rocket engines during shutdown

A linearized one-dimensional lumped-parameter model capable of predicting the occurrence of the low frequency combustion instability (chugging) experienced during preburner shutdown in the Space Shuttle Main Engines is discussed, and predictions are compared with NASA experimental results. Results from a parametric study of parameters including chamber pressure, fuel and oxygen temperatures, and the effective bulk modulus of the liquid oxidizer suggest that chugging is probably affected by conditions at shutdown through the fuel and oxidizer temperatures. It is suggested that chugging is initiated when the fuel, oxidizer, and helium temperature and flow rates pass into an unstable region, and that chugging may be terminated by decaying pressures.

Lim, Kair-Chuan↗
NASA NTRS1990

Computer-aided-analysis of linear control system robustness

Stochastic robustness is a simple technique used to estimate the stability and performance robustness of linear, time-invariant systems. The use of high-speed graphics workstations and control system design software in stochastic robustness analysis is discussed and demonstrated. It is shown that stochastic robustness makes good use of modern computational and graphic tools, and it is easily implemented using commercial control system design and analysis software.

Stengel, Robert F.↗
NASA NTRS1986

Incompressible boundary-layer stability analysis of LFC experimental data for sub-critical Mach numbers

An incompressible boundary-layer stability analysis of Laminar Flow Control (LFC) experimental data was completed and the results are presented. This analysis was undertaken for three reasons: to study laminar boundary-layer stability on a modern swept LFC airfoil; to calculate incompressible design limits of linear stability theory as applied to a modern airfoil at high subsonic speeds; and to verify the use of linear stability theory as a design tool. The experimental data were taken from the slotted LFC experiment recently completed in the NASA Langley 8-Foot Transonic Pressure Tunnel. Linear stability theory was applied and the results were compared with transition data to arrive at correlated n-factors. Results of the analysis showed that for the configuration and cases studied, Tollmien-Schlichting (TS) amplification was the dominating disturbance influencing transition. For these cases, incompressible linear stability theory correlated with an n-factor for TS waves of approximately 10 at transition. The n-factor method correlated rather consistently to this value despite a number of non-ideal conditions which indicates the method is useful as a design tool for advanced laminar flow airfoils.

Berry, S. A.↗
NASA NTRS1993

Analysis of implicit second-order upwind-biased stencils

Truncation error and stability properties of several implicit upwind schemes for the two-dimensional Euler equations are examined. The schemes use linear data reconstruction methods to achieve second-order flux integrations where the implicit Jacobian operators are first order. The stability properties of the schemes are examined by a Von Neumann analysis of the linearized, constant-coefficient Euler equations. The choice of the data reconstruction method used to evaluate the flux integral has a dramatic effect on the convergence properties of the implicit solution method. In particular, the typical one-dimensional data reconstruction methods used with structured grids exhibit poor convergence properties compared to the unstructured grid method considered. Of the schemes examined, the one with the superior convergence properties is well-suited for both unstructured and structured grids, which has important implications for the design of implicit methods.

Roberts, Thomas W.↗
NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stability of superthermal strahl electrons in the solar wind

We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Digital signal processing and control and estimation theory -- Points of tangency, area of intersection, and parallel directions

A number of current research directions in the fields of digital signal processing and modern control and estimation theory were studied. Topics such as stability theory, linear prediction and parameter identification, system analysis and implementation, two-dimensional filtering, decentralized control and estimation, image processing, and nonlinear system theory were examined in order to uncover some of the basic similarities and differences in the goals, techniques, and philosophy of the two disciplines. An extensive bibliography is included.

Willsky, A. S.↗
NASA NTRS1983

Morphological and convective instabilities during solidification

Linear stability theory is employed in the present analysis of flow stability between two vertical, infinite, rigid coaxial cylinders at different temperatures. These calculations have been prompted by, and are found to be in general agreement with, experiments on succinonitrile. A long, vertical cylinder sample of this material was heated so that a vertical melt annulus formed between the coaxial heater and the surrounding crystal/melt interface. Above a critical Grashof number of about 200, a helical crystal/melt interface formed which steadily rotated about the cylinder axis and whose wave speed was several orders of magnitude lower than the base flow velocity.

Coriell, S. R.↗
NASA NTRS1992

A comparison of the impact of two time-differencing schemes on the NASA-GLAS climate model

Evidence is presented of the sensitivity of a GCM to the time-differencing scheme employed when the physical parameterizations and space discretization are not changed. Two-time marching schemes, the leapfrog and the Matsuno schemes, are analyzed and tested on the NASA-Goddard Laboratory for Atmospheric Studies fourth-order GCM in terms of the behavior and stability of two-month-averaged fields. Linear analysis indicates that Rossby waves are slightly accelerated and slightly damped when the Matsuno scheme is utilized and that these effects are scale selective, being smallest for the longest waves.

Pfeffer, Richard L.↗
NASA NTRS1989

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

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

Dynamic response and stability analysis of flexible, multibody systems

A general version of Lagrange's equations, including auxiliary nonholonomic, rheonomic conditions of constraint, is used in the dynamic simulation and stability analysis of interconnected flexible bodies. Modeling of the nonlinear flexible/rigid dynamic coupling effects, the interaction forces/torques, and the elastic deformation effects is discussed. A digital computer program is developed to obtain time-domain solution for the nonlinear response of systems represented as a collection of individual bodies, numerical linearization of system-governing equations, time-domain solution for the perturbation response about a nominal state, and a frequency-domain stability analysis corresponding to the linearization. The digital simulation code is employed to study the dynamic behavior of a typical satellite and a spacecraft with deployable experiment booms.

Bodley, C. S.↗
NASA NTRSDate not supplied

Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations

Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.

Boundary Layer Transition↗
DOE OSTI2023

Stability and Performance of an Undersea Kite Operating in a Turbulent Flow Field

In this paper, we examine the effects of flow disturbances resulting from turbulence on the dynamic behavior of an underwater energy-harvesting kite system that executes periodic figure-8 flight. Due to the periodic nature of the kite’s operation, we begin by assessing orbital stability using the Floquet analysis and stroboscopic intersection analysis of a Poincaré section, with the former analysis performed on a simplified “unifoil” model and the latter performed on a six-degree-of-freedom (6-DOF)/flexible tether model. With periodic stability established, a frequency-domain analysis based on a linearization about the kite’s path is used to predict the quality of flight path tracking as a function of the turbulence frequency. To validate the accuracy of these simulation-based predictions under flow disturbances, we compare the predictions of the kite’s behavior against the results of small-scale tow testing experiments performed in a controlled pool environment.

13 HYDRO ENERGY↗
NASA NTRS1988

An averaging analysis of discrete-time indirect adaptive control

An averaging analysis of indirect, discrete-time, adaptive control systems is presented. The analysis results in a signal-dependent stability condition and accounts for unmodeled plant dynamics as well as exogenous disturbances. This analysis is applied to two discrete-time adaptive algorithms: an unnormalized gradient algorithm and a recursive least-squares (RLS) algorithm with resetting. Since linearization and averaging are used for the gradient analysis, a local stability result valid for small adaptation gains is found. For RLS with resetting, the assumption is that there is a long time between resets. The results for the two algorithms are virtually identical, emphasizing their similarities in adaptive control.

Phillips, Stephen M.↗