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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Analysis and Construction of Perfectly Matched Layers for the Linearized Euler Equations

We present a detailed analysis of a recently proposed perfectly matched layer (PML) method for the absorption of acoustic waves. The split set of equations is shown to be only weakly well-posed, and ill-posed under small low order perturbations. This analysis provides the explanation for the stability problems associated with the split field formulation and illustrates why applying a filter has a stabilizing effect. Utilizing recent results obtained within the context of electromagnetics, we develop strongly well-posed absorbing layers for the linearized Euler equations. The schemes are shown to be perfectly absorbing independent of frequency and angle of incidence of the wave in the case of a non-convecting mean flow. In the general case of a convecting mean flow, a number of techniques is combined to obtain a absorbing layers exhibiting PML-like behavior. The efficacy of the proposed absorbing layers is illustrated though computation of benchmark problems in aero-acoustics.

Hesthaven, J. S.↗
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Attitude stability of a spinning spacecraft during appendage deployment/retraction

The work presented is motivated by the need for a national satellite rescue policy, not the ad hoc policy now in place. In studying different approaches for a national policy, the issue of capture and stabilization of a tumbling spacecraft must be addressed. For a rescue mission involving a tumbling spacecraft, it may be advantageous to have a rescue vehicle which is compact and 'rigid' during the rendezvous/capture phase. After capture, passive stabilization techniques could be utilized as an efficient means of detumbling the resulting system (i.e., both the rescue vehicle and captures spacecraft). Since the rescue vehicle is initially compact and 'rigid,' significant passive stabilization through energy dissipation can only be achieved through the deployment of flexible appendages. Once stabilization is accomplished, retraction of the appendages before maneuvering the system to its final destination may also prove advantageous. It is therefore of paramount interest that we study the effect of appendage deployment/retraction on the attitude stability of a spacecraft. Particular interest should be paid to appendage retraction, since if this process is destabilizing, passive stabilization as proposed may not be useful. Over the past three decades, it has been an 'on-again-off-again affair' with the problem of spacecraft appendage deployment. In most instances, these studies have been numerical simulations of specific spacecraft configurations for which there were specific concerns. The primary focus of these studies was the behavior of the appendage during deployment; the effects of appendage retraction was considered only in one of these studies. What is missing in the literature is a thorough study of the effects of appendage deployment/retraction on the attitude stability of a spacecraft. This paper presents a rigorous analysis of the stability of a spinning spacecraft during the deployment or the retraction of an appendage. The analysis is simplified such that meaningful insights into the problem can be inferred; it is not overly simplified such that critical dynamical behavior is neglected. The system is analyzed assuming that the spacecraft hub is rigid. The appendage deployment mechanism is modeled as a point mass on a massless rod whose length undergoes prescribed changes. Simplified flexibility effects of the appendage are included. The system is examined for stability by linearizing the equations in terms of small deviations from steady, noninterfering coning motion. Routh's procedure for analyzing small deviations from steady motion in dynamical systems is utilized in the analysis. The system of equations are nondimensionalized to facilitate parametric studies. The results are presented in terms of a reduced number of nondimensional parameters so that some general conclusions may be drawn. Verification of the linear analysis is presented through numerical simulations of the complete nonlinear, nonautonomous, coupled equations.

Fitz-Coy, Norman↗
NASA NTRS1994

Cooperative control theory and integrated flight and propulsion control

This report documents the activities and research results obtained under a grant (NAG3-998) from the NASA Lewis Research Center. The focus of the research was the investigation of dynamic interactions between airframe and engines for advanced ASTOVL aircraft configurations, and the analysis of the implications of these interactions on the stability and performance of the airframe and engine control systems. In addition, the need for integrated flight and propulsion control for such aircraft was addressed. The major contribution of this research was the exposition of the fact that airframe and engine interactions could be present, and their effects could include loss of stability and performance of the control systems. Also, the significance of two directional, as opposed to one-directional, coupling was identified and explained. A multi variable stability and performance analysis methodology was developed, and applied to several candidate aircraft configurations. Also exposed was the fact that with interactions present along with some integrated control approaches, the engine command/limiting logic (which represents an important non-linear component of the engine control system) can impact closed-loop airframe/engine system stability. Finally, a brief investigation of control-law synthesis techniques appropriate for the class of systems was pursued, and it was determined that multi variable techniques, included model-following formulations of LQG and/or H (infinity) methods showed promise. However, for practical reasons, decentralized control architectures are preferred, which is an architecture incompatible with these synthesis methods.

Schmidt, David K.↗
NASA NTRS1972

A generalized method for the identification of aircraft stability and control derivatives from flight test data.

This paper discusses the application of a generalized identification method for flight test data analysis. The method is based on the maximum likelihood (ML) criterion and includes output error and equation error methods as special cases. Both the linear and nonlinear models with and without process noise are considered. The flight test data from lateral maneuvers of HL-10 and M2/F3 lifting bodies are processed to determine the lateral stability and control derivatives, instrumentation accuracies and biases. A comparison is made between the results of the output error method and the generalized ML method for M2/F3 data containing gusts. It is shown that better fits to time histories are obtained by using the generalized ML method.

Mehra, R. K.↗
NASA NTRS1996

Nonlinear stability of supersonic jets

The stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations are presented. In this analysis the large scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Tiwari, S. N.↗
NASA NTRS1995

Nonlinear Stability of Supersonic Jets

This paper presents stability calculations made for a shock-free supersonic jet using the model based on parabolized stability equations. In this analysis the large-scale structures, which play a dominant role in the mixing as well as the noise radiated, are modeled as instability waves. This model takes into consideration non-parallel flow effects and also nonlinear interaction of the instability waves. The stability calculations have been performed for different frequencies and mode numbers over a range of jet operating temperatures. Comparisons are made, where appropriate, with the solutions to Rayleigh's equation (linear, inviscid analysis with the assumption of parallel flow). The comparison of the solutions obtained using the two approaches show very good agreement.

Bhat, T. R. S.↗
NASA NTRS1973

Energy-like Liapunov functionals for linear elastic systems on a Hilbert space.

An approach is presented for generating energy-like functionals for linear elastic dynamic systems on a Hilbert space. The objective is to obtain a family of functionals which may be used for stability analysis of the equilibrium, i.e., Liapunov functionals. Although the energy functional, when one exists, is always a member of this family, the family is shown to exist even when an energy functional does not. Several discrete and distributed-parameter examples are presented, as are certain specific techniques for utilizing this approach.

Walker, J. A.↗
NASA NTRS1977

Stability and control of maneuvering high-performance aircraft

The stability and control of a high-performance aircraft was analyzed, and a design methodology for a departure prevention stability augmentation system (DPSAS) was developed. A general linear aircraft model was derived which includes maneuvering flight effects and trim calculation procedures for investigating highly dynamic trajectories. The stability and control analysis systematically explored the effects of flight condition and angular motion, as well as the stability of typical air combat trajectories. The effects of configuration variation also were examined.

Stengel, R. F.↗
NASA NTRS1992

Stochastic robustness synthesis for a benchmark problem

Stochastic robustness analysis guides the synthesis of robust linear quadratic Gaussian (LQG) regulators for a benchmark control problem. Probabilities of exceeding allowable design limits, including stability, setting time, and control usage, are estimated by Monte Carlo evaluation. Robust, low-gain compensators that fulfill objectives are designed by numerically minimizing quadratic functions of these probabilities. The method is straightforward and makes use of uncomplicated design principles.

Stengel, Robert F.↗
NASA NTRS1972

Longitudinal dynamic stability of a shuttle vehicle.

Analytical study of the longitudinal dynamic stability of a nonrolling, lifting vehicle gliding at hypersonic speeds. The analysis applies to shuttle vehicles designed for operating up to the rim of a planetary atmosphere. A general nondimensional time transformation is introduced to derive a unified second-order linear differential equation for the angle of attack, valid for all types of reentry of a general type of vehicle. The stability of motion is discussed for two fundamental regimes of flight that are based on widely different assumptions. For near ballistic entry along a straight line trajectory, the equation reduces to a confluent hypergeometric equation, the solution of which can be expressed in terms of Whittaker's function. Using a theorem in the theory of stability of differential equations, criteria for damped oscillations are derived. It is shown that the aerodynamic criteria for stability are the same as for the case of ballistic entry. In addition, for each vehicle configuration, and specified planetary atmosphere, there exists an altitude range where the angle of attack frequency is nearly equal to the orbital frequency causing instability in pitch. This resonance instability is due to the ellipticity of the orbit. Criteria for eccentricity instability are derived.

Vinh, N. X.↗
NASA NTRS1977

Finite Element Stability Analysis for Coupled Rotor and Support Systems

The effects of fuselage motions on stability and random response were analytically assessed. The feasibility of adequate perturbation models from non-linear trim conditions was studied by computer and hardware experiments. Rotor wake-blade interactions were assessed by using a 4-bladed rotor model with the capability of progressing and regressing blade pitch excitation (cyclic pitch stirring), by using a 4-bladed rotor model with hub tilt stirring, and by testing rotor models in sinusoidal up or side flow.

Hohenemser, K. H.↗
NASA NTRS1978

Some observations on four current subjects related to aeroelastic stability

After introductory comments on the literature and the purposes of this paper, a table is presented summarizing the author's views on some currently solved vs partially unsolved problems related to aeroelastic stability. The term 'solved' is used in the practical sense that engineers are able to cope confidently with that problem during the process of structural design. Selected entries in the table are reviewed, partially to motivate the topics in the rest of the paper. The 'four current subjects' are chosen both for timeliness and because they are among the ongoing interests of the Stanford group. The first involves the prediction of linearized unsteady aerodynamic loads due to arbitrary motions of streamlined shapes. Some contributions by Edwards are refined, which were motivated by the requirements of active control system design. The second subject is nonlinear unsteady aerodynamics for the transonic regime. After describing a few useful developments from locally-linear theory and computational fluid dynamics, there is suggested an empirical procedure for interim-analysis purposes. The third and fourth subjects concern recent discoveries regarding the aeroelastic stability of large-aspect-ratio wings and wind turbines. The former work is mainly that of Petre and Boyd. The latter includes some of the author's own preliminary discoveries about the performance and dynamics of vertical-axis machines.

Ashley, H.↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Analyzing Aeroelasticity in Turbomachines

ASTROP2-LE is a computer program that predicts flutter and forced responses of blades, vanes, and other components of such turbomachines as fans, compressors, and turbines. ASTROP2-LE is based on the ASTROP2 program, developed previously for analysis of stability of turbomachinery components. In developing ASTROP2- LE, ASTROP2 was modified to include a capability for modeling forced responses. The program was also modified to add a capability for analysis of aeroelasticity with mistuning and unsteady aerodynamic solutions from another program, LINFLX2D, that solves the linearized Euler equations of unsteady two-dimensional flow. Using LINFLX2D to calculate unsteady aerodynamic loads, it is possible to analyze effects of transonic flow on flutter and forced response. ASTROP2-LE can be used to analyze subsonic, transonic, and supersonic aerodynamics and structural mistuning for rotors with blades of differing structural properties. It calculates the aerodynamic damping of a blade system operating in airflow so that stability can be assessed. The code also predicts the magnitudes and frequencies of the unsteady aerodynamic forces on the airfoils of a blade row from incoming wakes. This information can be used in high-cycle fatigue analysis to predict the fatigue lives of the blades.

Reddy, T. S. R.↗
NASA NTRS2022

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow↗
DOE OSTI2022

Bifurcation Buckling Analysis and Non-Linear Collapse Analysis of Teardrop Shaped Vacuum Chamber

The Spallation Neutron Source (SNS) accelerator is being upgraded to increase the beam power from 1.4MW at 1GeV to 2.8MW at 1.3GeV. The currents in the middle two injection chicane magnets cannot simply be scaled up to accommodate the increased injection energy of 1.3GeV due to potential excessive H- stripping; the magnets must be replaced with longer, lower-field magnets and the associated vacuum chambers need to be redesigned. A new teardrop-shaped vacuum chamber was initially designed to accommodate the new magnets and the updated beam paths and instrumentation. This paper focuses on the structural stability study of the teardrop shape vacuum chamber based on buckling analysis. Protection against collapse from buckling according to the ASME BPVC requirement has been evaluated in depth. First, a Type-1 bifurcation buckling analysis using a linear eigenvalue solution to determine the critical load factor was performed. Subsequently, a Type-3 nonlinear collapse analysis was conducted using the static Riks method with elastic-plastic material properties and imperfections explicitly considered in the model geometry. The critical buckling load for the teardrop shape vacuum chamber was confidently estimated based upon this two-stage approach.

Jiang, Hao↗
NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗