Search NASA⌕ Search

SEARCH · Search NASA

Results for “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

Stability of the stratifield cylindrical annulus flow

The linear stability analysis for the stratified flow between two rotating circular cylinders is formulated. Two approaches for the stability analysis are presented. The first approach results in an algebraic eigenvalue problem, while the second results in an initial value problem for the perturbation function. The advantages and disadvantages of both approaches are discussed and a preferable numerical solution technique is outlined.

Antar, B. N.↗

Aeroelastic effects in multirotor vehicles. Part 2: Methods of solution and results illustrating coupled rotor/body aeromechanical stability

This report is a sequel to the earlier report titled, Aeroelastic Effects in Multi-Rotor Vehicles with Application to Hybrid Heavy Lift System, Part 1: Formulation of Equations of Motion (NASA CR-3822). The trim and stability equations are presented for a twin rotor system with a buoyant envelope and an underslung load attached to a flexible supporting structure. These equations are specialized for the case of hovering flight. A stability analysis, for such a vehicle with 31 degrees of freedom, yields a total of 62 eigenvalues. A careful parametric study is performed to identify the various blade and vehicle modes, as well as the coupling between various modes. Finally, it is shown that the coupled rotor/vehicle stability analysis provides information on both the aeroelastic stability as well as complete vehicle dynamic stability. Also presented are the results of an analytical study aimed at predicting the aeromechanical stability of a single rotor helicopter in ground resonance. The theoretical results are found to be in good agreement with the experimental results, thereby validating the analytical model for the dynamics of the coupled rotor/support system.

Venkatesan, C.↗

Perturbation solutions of combustion instability problems

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Googerdy, A.↗

A perturbation solution for transverse pressure-sensitive combustion instability in an annulus

A method for analysis of combustion instability in rocket motors based on a combination of the Galerkin method and the two-variable (or multiple-scales) perturbation method is developed. The method is illustrated by applying it to the problem of pressure-sensitive combustion instability in a liquid-fuel annular combustion chamber. To the order of approximation inherent in the method it is found that the complete linear stability analysis and most of the nonlinear stability analysis can be carried out in closed form. The results thus obtained are used to clarify several aspects of the interpretation of numerical solutions obtained by previous investigators.

Googerdy, A.↗

The stability of tangential discontinuities

Intuition regarding the stability of MHD tangential discontinuities of solar wind has been borne out by a linear stability analysis of the MHD equations. In performing the calculation, the model used had no plasma flow and thus corresponded to the solar wind rest frame. The method calculated the change in total energy, plasma plus magnetic, due to the perturbation. This type of stability analysis is analogous to the test of the stability of an equilibrium point in a mechanical potential.

Birmingham, T.↗

Numerical Simulations of Unsteady Low-Reynolds-Number Flows Over the APEX Airfoil

Laminar and transitional separation bubbles are an important feature of low-Reynolds-number flows over airfoils. The separation bubbles are unsteady and have a significant impact on the aerodynamic properties of the airfoils. In this paper unsteady low-Reynolds-number separated flows over the APEX airfoil are calculated using a Navier-Stokes solver. The numerical results show the presence of unsteady separation bubbles in the flowfield. An analysis of the numerical results shows that flowfield disturbances are amplified significantly in the separation bubble, leading to periodic vortex shedding. A linear stability analysis of the separated boundary layer is performed and the results show that the dominant wavenumber and frequency in the numerical simulations agree with the most unstable wavenumber and frequency from the linear stability analysis. The numerical results also show the growth and interaction of disturbance waves in the separation bubble. For transonic flows over the APEX airfoil, the calculations show that the presence of shocks causes significant changes in the separation location and consequently, the overall flowfield.

Tatineni, Mahidhar↗

Robust controllers for the Middeck Active Control Experiment using Popov controller synthesis

Recent work in robust control with real parameter uncertainties has focused on absolute stability and its connections to real mu theory. In particular, the research has investigated the Popov stability criterion and its associated Lur'e-Postnikov Liapunov functions. State space representations of this Popov stability analysis tests are included in an H2 design formulation to provide a powerful technique for robust controller synthesis. This synthesis approach uses a state space optimization procedure to design controllers that minimize an overbound of an H2 cost functional and satisfy stability analysis tests based on the Popov multiplier. The controller and stability multiplier coefficients are optimized simultaneously, which avoids the iteration and curve-fitting procedures required by the D-K algorithm of mu synthesis. While previous work has demonstrated this synthesis approach on benchmark control problems, the purpose of this paper is to use Popov controller synthesis to design robust compensators for the Middeck Active Control Experiment (MACE).

How, Jonathan P.↗

Linear Stability Regime Analysis of the Compressible Reacting Mixing Layer

Previous investigations have shown that a compressible reacting mixing layer can develop two peaks in the mean density weighted vorticity profile. Linear stability analyses show that at these peaks two distinct 'outer' instability modes appear in addition to the more common central mode, which exists unaccompanied in incompressible nonreacting flows. The present study parametrically analyzes the effects of compressibility, heat release, stoichiometry, and density ratio on the amplification rate and obliquity of each stability mode. The mean profiles used in the spatial stability calculation are generated by self-similar solutions of the compressible boundary layer equations combined with the assumption of infinitely fast chemistry. It is shown that the influence of stoichiometry and density ratio on the peaks of the density weighted vorticity profile determines which modes will dominate. Of particular interest are the conditions where two modes are equally amplified, causing the mixing layer to develop into a 'colayer' structure.

Day, M. J.↗

Analysis of stability of a planar solid-liquid interface in a dilute binary alloy

This paper reconsiders the question of stability of a planar solid-liquid interface in an undercooled alloy melt without making the restrictive assumption of no heat flow in the solid (i.e., Gs = 0). The results of this analysis indicate that, provided the thermal gradient on the solid side of the interface, Gs, is positive, stability can be achieved in an undercooled alloy melt for growth rates R greater than Ra (where Ra is the absolute stability limit of Mullins and Sekerka, 1964). Thus, the absolute stability criterion for steady-state planar growth in an undercooled alloy melt is the same as derived earlier by Mullins and Sekerka for directional solidification. Relaxing the restrictive assumption of Gs = 0 also reveals that there is a regime of stability for low growth rates and low supercoolings.

Laxmanan, V.↗

Flutter Analysis with Stabilized Finite Elements Based on the Linearized Frequency-Domain Approach

When designing and certifying aircraft, engineers must take into consideration aeroelastic effects such as flutter. Design and certification of a vehicle may require analysis of thousands of aeroelastic responses. Standard tools in the aerospace industry are based on linear aerodynamic models such as the doublet-lattice method, but these methods can be nonconservative in certain situations such as in the transonic regime. While computational fluid dynamics (CFD) is a higher fidelity alternative, the time-marching approach has a drastically increased computational cost compared to the linear aerodynamic methods. By taking advantage of the periodic nature of flutter, frequency-domain methods offer a more efficient alternative to time-marching CFD. In this work, a linearized frequency-domain method is implemented and verified in the stabilized finite-element solver in FUN3D. The linearized frequency-domain method is demonstrated and compared to other methods for traditional benchmark cases for computational aeroelasticity: the AGARD 445.6 wing, the Benchmark Supercritical Wing, and the Benchmark NACA 0012Wing.

Kevin E Jacobson↗

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗

Nonlinear stability and control study of highly maneuverable high performance aircraft, phase 2

This research should lead to the development of new nonlinear methodologies for the adaptive control and stability analysis of high angle-of-attack aircraft such as the F18 (HARV). The emphasis has been on nonlinear adaptive control, but associated model development, system identification, stability analysis and simulation is performed in some detail as well. Various models under investigation for different purposes are summarized in tabular form. Models and simulation for the longitudinal dynamics have been developed for all types except the nonlinear ordinary differential equation model. Briefly, studies completed indicate that nonlinear adaptive control can outperform linear adaptive control for rapid maneuvers with large changes in alpha. The transient responses are compared where the desired alpha varies from 5 degrees to 60 degrees to 30 degrees and back to 5 degrees in all about 16 sec. Here, the horizontal stabilator is the only control used with an assumed first-order linear actuator with a 1/30 sec time constant.

Mohler, R. R.↗

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.

Antar, B. N.↗

New Directions in Asymptotically Stable Finite-dimensional Adaptive Control of Linear Distributed Parameter Systems

Distributed Parameter Systems (DPS), such as systems described by partial differential equations, require infinite-dimensional state space descriptions to correctly model their dynamical behavior. However, any adaptive control algorithm must be finite-dimensional in order to be implemented via on-line digital computers. Finite-dimensional adaptive control of linear DPS requires stability analysis of nonlinear, time-varying, infinite-dimensional systems. The structure of nonadaptive finite-dimensional control of linear DPS is summarized as it relates to the existence of limiting systems for adaptive control. Two candidate schemes for finite-dimensional adaptive control of DPS are described and critical issues in infinite-dimensional stability analysis are discussed, in particular, the invariance principle, center manifold theory, and relationships between input-output and internal stability.

Balas, M. J.↗

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗