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 811 records · Page 45

Extended Lyapunov stability criterion using a nonlinear algebraic relation with application to adaptive control

A new nonlinear stability criterion is developed by use of a class of Lyapunov functionals for model-reference adaptive systems (MRAS). Results are compared with traditional results, and a comparative design technique is used to illustrate its function in improving the transient response of an MRAS controller. For a particular system structure and class of input signals, the new stability criterion is shown to include traditional sufficiency stability conditions as a special case. An example is cited to illustrate the use of the nonlinear criterion and its definite advantages in helping improve the adaptive error transient response of a system. Analysis of results is effected by use of a linearization technique on the resulting adaptive equations.

Colburn, B. K.↗

Unsteady flow in a supersonic cascade with strong in-passage shocks

Linearized theory is used to study the unsteady flow in a supersonic cascade with in-passage shock waves. We use the Wiener-Hopf technique to obtain a closed-form analytical solution for the supersonic region. To obtain a solution for the rotational flow in the subsonic region we must solve an infinite set of linear algebraic equations. The analysis shows that it is possible to correlate quantitatively the oscillatory shock motion with the Kutta condition at the trailing edges of the blades. This feature allows us to account for the effect of shock motion on the stability of the cascade. Unlike the theory for a completely supersonic flow, the present study predicts the occurrence of supersonic bending flutter. It therefore provides a possible explanation for the bending flutter that has recently been detected in aircraft-engine compressors at higher blade loadings.

Goldstein, M. E.↗

Toward a unified theory of atmospheric convective instability

A nonlinear three-dimensional truncated spectral model of shallow and moist Boussinesq convection indicates that parallel instability and thermal forcing are linked, in view of the fact that only one convective mode exists in which either or both mechanisms are operating to generate convection in the planetary boundary layer. It is also established that the wind field causes two-dimensional roll convection formation, an alignment of the convection with the wind in a preferred manner, and a propagation speed that is related to the wind component perpendicular to the roll axis. Latent heating is responsible for the decrease of the critical value of the environmental lapse rate in accordance with the slice method stability criterion. When only the upper part of the upward branch is moist and all of the downward branch is dry, latent heating also causes a finite-amplitude convective solution for Rayleigh number values lower than the critical value of linear analysis.

Shirer, H. N.↗

Issues in modeling and controlling the SCOLE configuration

The parametric study of the in-plane Spacecraft Control Laboratory Experiment (SCOLE) system, the Floquet Stability Analysis, and three dimensional formulations of the SCOLE system dynamics are examined. Control issues are discussed, such as: control of large structures with delayed input in continuous time; control with delayed input in discrete time; control law design for SCOLE using Linear Quadratic Gaussian (LQC)/TRR technique; and optimal torque control for SCOLE slewing maneuvers.

Bainum, Peter M.↗

Dynamic evaluation of a traction-drive joint for space telerobots

Presented is an analysis and evaluation of a prototype traction-drive joint for robotic manipulators, developed under NASA sponsorship. A dynamic model is developed using the Lagrange formulation. Controllability, observability, dynamic stability, and response characteristics of the joint to test inputs are studied. A linear quadratic regulator (LQR) is implemented on the joint model to determine a basis for evaluating the performance of the traction-drive joint under servo control. An evaluation of the results and directions for future investigations are presented.

Desilva, Clarence W.↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Practical Methodology for the Inclusion of Nonlinear Slosh Damping in the Stability Analysis of Liquid-Propelled Space Vehicles

One of the challenges of developing flight control systems for liquid-propelled space vehicles is ensuring stability and performance in the presence of parasitic minimally damped slosh dynamics in the liquid propellants. This can be especially difficult when the fundamental frequencies of the slosh motions are in proximity to the frequency used for vehicle control. The challenge is partially alleviated since the energy dissipation and effective damping in the slosh modes increases with amplitude. However, traditional launch vehicle control design methodology is performed with linearized systems using a fixed slosh damping corresponding to a slosh motion amplitude based on heritage values. This papers presents a method for performing the control design and analysis using damping at slosh amplitudes chosen based on the resulting limit cycle amplitude of the vehicle thrust vector system due to a control-slosh interaction under degraded phase and gain margin conditions.

Ottander, John A.↗

Practical Methodology for the Inclusion of Nonlinear Slosh Damping in the Stability Analysis of Liquid-propelled Space Vehicles

One of the challenges of developing flight control systems for liquid-propelled space vehicles is ensuring stability and performance in the presence of parasitic minimally damped slosh dynamics in the liquid propellants. This can be especially difficult when the fundamental frequencies of the slosh motions are in proximity to the frequency used for vehicle control. The challenge is partially alleviated since the energy dissipation and effective damping in the slosh modes increases with amplitude. However, traditional launch vehicle control design methodology is performed with linearized systems using a fixed slosh damping corresponding to a slosh motion amplitude based on heritage values. This papers presents a method for performing the control design and analysis using damping at slosh amplitudes chosen based on the resulting limit cycle amplitude of the vehicle thrust vector system due to a control-slosh interaction under degraded phase and gain margin conditions.

Ottander, John A.↗

Hypersonic Transition Analysis for HIFiRE Experiments

The HIFiRE-1 flight experiment provided a valuable database pertaining to boundary layer transition over a 7-degree half-angle, circular cone model from supersonic to hypersonic Mach numbers, and a range of Reynolds numbers and angles of attack. This paper reports selected findings from the ongoing computational analysis of the measured in-flight transition behavior. Transition during the ascent phase at nearly zero degree angle of attack is dominated by second mode instabilities except in the vicinity of the cone meridian where a roughness element was placed midway along the length of the cone. The growth of first mode instabilities is found to be weak at all trajectory points analyzed from the ascent phase. For times less than approximately 18.5 seconds into the flight, the peak amplification ratio for second mode disturbances is sufficiently small because of the lower Mach numbers at earlier times, so that the transition behavior inferred from the measurements is attributed to an unknown physical mechanism, potentially related to step discontinuities in surface height near the locations of a change in the surface material. Based on the time histories of temperature and/or heat flux at transducer locations within the aft portion of the cone, the onset of transition correlated with a linear N-factor, based on parabolized stability equations, of approximately 13. Due to the large angles of attack during the re-entry phase, crossflow instability may play a significant role in transition. Computations also indicate the presence of pronounced crossflow separation over a significant portion of the trajectory segment that is relevant to transition analysis.

Li, Fei↗

Boundary-layer stability analysis of NLF and LFC experimental data at subsonic and transonic speeds

NASA-Lewis has conducted wind tunnel experiments to ascertain the effectiveness of state-of-the-art in natural laminar flow (NLF) and LFC airfoils for subsonic and transonic speeds, such as the NLF(1)-0414F and the SCLFC(1)-0513F. Attention is given to the effects of Tollmien-Schlichting (TS) and/or crossflow linear mechanisms amplifying small disturbances to generate turbulence. It is found that the incompressible TS transitional n-factors were generally in the 9-12 range, in agreement with earlier correlation studies; the TS instability was the dominant instability mode on a swept-planform LFC airfoil over the entire range of test conditions.

Berry, Scott A.↗

Stability analysis of double-peaked proton distribution functions in the solar wind

Twelve proton velocity distribution functions as measured for different solar wind speeds by the plasma analyzer at the Helios spacecraft are modeled analytically by a superposition of two bi-Maxwellians. The stability of these distributions is investigated by solving the linear dispersion equation for parallel propagation of both left- and right-hand polarized modes. Due to the temperature anisotropy of the main proton component the presence of the resonant left-hand ion cyclotron instability is found to be typical for high-speed streams and is independent of an ion beam. The results indicate that the ratio of the beam density to total particle density and the beam drift velocity are generally close to the threshold values necessary to drive a highly unstable nonresonant right-hand mode. This instability is found for one proton distribution and is a consequence of the beam particles. It is concluded that both the ion cyclotron instability and the ion beam instability act as a regulating mechanism on the shape of solar wind double humped proton distribution functions.

Leubner, Manfred P.↗

Analysis of Alfven eigenmode destabilization in ITER using a Landau closure model

Alfvén eigenmodes (AE) can be destabilized during ITER discharges driven by neutral beam injection (NBI) energetic particles (EP) and alpha particles. The aim of the present study is to analyze the AE stability of different ITER operation scenarios considering multiple energetic particle species. We use the reduced magneto-hydrodynamic (MHD) equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the EP species including the effect of the acoustic modes. Furthermore, the AEs driven by the NBI EP and alpha particles are stable in the configurations analyzed, only MHD-like modes with large toroidal couplings are unstable, although both can be destabilized if the EP increases above a threshold. The threshold is two times the model value for the NBI EP and alpha particles in the reverse shear (RS) case, leading to the destabilization of Beta induced AE (BAE) near the magnetic axis with a frequency of kHz and toroidal or elliptical AE (TAE/EAE) in the RS region with a frequency of kHz, respectively. On the other hand, the hybrid and steady state configurations show a threshold 3 times larger with respect to the model for the alpha particle and 40 times for the NBI EP, also destabilizing BAE and TAE between the inner and middle plasma region. In addition, a extended analysis of the RS scenario where the of both alpha particles and NBI EP are above the AE threshold, multiple EP damping effects are also identified as well as optimization trends regarding the resonance properties of the alpha particle and NBI EP with the bulk plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analyzing Impact of Distributed PV Generation on Integrated Transmission & Distribution System Voltage Stability—A Graph Trace Analysis Based Approach

The use of a Graph Trace Analysis (GTA)-based power flow for analyzing the voltage stability of integrated Transmission and Distribution (T&D) networks is discussed in the context of distributed Photovoltaic (PV) generation. The voltage stability of lines and the load carrying capability of buses is analyzed at various PV penetration levels. It is shown that as the PV generation levels increase, an increase in the steady state voltage stability of the system is observed. Moreover, within certain regions of stability margin changes, changes in voltage stability margins of transmission lines are shown to be linearly related to changes in the loading of the lines. Two case studies are presented, where one case study involves a model with eight voltage levels and 784,000 nodes. In one case study, a voltage-stability heat map is used to demonstrate the identification of weak lines and buses.

14 SOLAR ENERGY↗

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

Formation Control of the MAXIM L2 Libration Orbit Mission

The Micro-Arcsecond X-ray Imaging Mission (MAXIM), a proposed concept for the Structure and Evolution of the Universe (SEU) Black Hole Imager mission, is designed to make a ten million-fold improvement in X-ray image clarity of celestial objects by providing better than 0.1 micro-arcsecond imaging. Currently the mission architecture comprises 25 spacecraft, 24 as optics modules and one as the detector, which will form sparse sub-apertures of a grazing incidence X-ray interferometer covering the 0.3-10 keV bandpass. This formation must allow for long duration continuous science observations and also for reconfiguration that permits re-pointing of the formation. To achieve these mission goals, the formation is required to cooperatively point at desired targets. Once pointed, the individual elements of the MAXIM formation must remain stable, maintaining their relative positions and attitudes below a critical threshold. These pointing and formation stability requirements impact the control and design of the formation. In this paper, we provide analysis of control efforts that are dependent upon the stability and the configuration and dimensions of the MAXIM formation. We emphasize the utilization of natural motions in the Lagrangian regions to minimize the control efforts and we address continuous control via input feedback linearization (IFL). Results provide control cost, configuration options, and capabilities as guidelines for the development of this complex mission.

Folta, David↗

Whirl Motion of a Seal Test Rig with Squeeze-Film Dampers

This paper presents the experimental behavior and dynamic analysis of a high speed test rig with rolling element bearings mounted in squeeze film oil damper bearings. The test rotor is a double overhung configuration with rolling element ball bearings mounted in uncentered squeeze-film oil dampers. The damper design is similar to that employed with various high-speed aircraft HP gas turbines. The dynamic performance of the test rig with the originally installed dampers with an effective damper length of length 0.23-inch was unacceptable. The design speed of 40,000 RPM could not be safely achieved as nonsynchronous whirling at the overhung seal test disk and high amplitude critical speed response at the drive spline section occurred at 32,000 RPM. In addition to the self excited stability and critical speed problems, it was later seen from FFT data analysis, that a region of supersynchronous dead band whirling occurs between 10,000 to 15,000 RPM which can lead to bearing distress and wear. The system was analyzed using both linear and nonlinear techniques. The extended length damper design resulting from the analysis eliminated the rotor subsynchronous whirling, high amplitude critical speed, and the dead band whirling region allowing the system to achieve a speed of 45,000 RPM. However, nonlinear analysis shows that damper lockup could occur with high rotor unbalance at 33,000 RPM, even with the extended squeeze-film dampers. The control of damper lockup will be addressed in a future paper.

Proctor, Margaret P.↗