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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

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

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

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↗

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

NASA NTRS1968

Analysis of Saturn booster performance by conventional control theory techniques Volume 1 of final report, 4 May 1967 - 3 Nov. 1968

Analysis of Saturn booster flight control system stability by linear feedback theory with emphasis on root locus technique

Smith, E. J.↗

NASA NTRS1974

Decomposition-aggregation stability analysis of the spinning Skylab

Stability of an 11-th order linear model of the spinning Skylab is determined by the decomposition-aggregation method based upon the comparison principle and vector Liapunov functions. To reduce the inherent conservativeness of the method an optimization problem is formulated and resolved producing the optimum comparison system. The system provides the best estimate of the stability region of the important structural parameter - asymmetry in the boom settings.

Cuk, S. M.↗

NASA NTRS1991

A geometric approach to regulator and tracker design for an aerospace plane

The paper presents a nonlinear design approach drawing from singular perturbations, feedback linearization, and variable structure control, that leads to regulators with automatic gain scheduling which exhibit similar dynamic behavior over the entire flight envelope of the aerospace plane. Additionally, design approach provides for a systematic way to counter disturbance effects as well as modeling uncertainties. The unifying feature of the three nonlinear feedback control methodologies is that they all have a geometric interpretation. First, the translational dynamics are decomposed into reduced-order slow and fast dynamics by way of a formal singular perturbation analysis. After feedback linearization the fast dynamics are robustly stabilized via a variable structure control approach. The slow dynamics are stabilized using conventional proportional-integral compensation based on the nominal slow dynamics. A number of sample command and disturbance responses at opposite ends of the flight envelope are presented for a nonlinear aerospace plane model.

Van Buren, Mark A.↗

NASA NTRS1992

Secondary instability in three-dimensional magnetic reconnection

We consider the transition to turbulence in three-dimensional reconnection of a magnetic neutral sheet. We find that the transition can occur via a three-step process. First, the sheet undergoes the usual tearing instability. Second, the tearing mode saturates to form a two-dimensional quasi-steady state. Third, this secondary equilibrium is itself unstable when it is perturbed by three-dimensional disturbances. Most of this paper is devoted to the analysis and simulation of the three-dimensional linear stability properties of the two-dimensional saturated tearing layer. The numerical simulations are performed with a semi-implicit, pseudospectral-Fourier collocation algorithm. We identify a three-dimensional secondary linear stability which grows on the ideal timescale. An examination of the modal energetics reveals that the largest energy transfer is from the mean field to the three-dimensional field, with the two-dimensional field acting as a catalyst.

Dahlburg, R. B.↗

NASA NTRS1996

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗

NASA NTRS1962

Use of an Inertia Sphere to Damp the Angular Motions of Spinning Space Vehicles

The general moment equations for a spin-stabilized vehicle with an inertia-reaction angular rate damper were considered, and it was noted that simplification would result if the damper had a spherical inertia distribution. A control system incorporating such a damper was postulated. The resulting equations were linearized, and conditions for stability were obtained from an analysis of the cubic characteristic equation. Two numerical examples were included.

Suddath, Jerrold H.↗

NASA NTRS2000

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗