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 271 records · Page 15

A first-order Lyapunov robustness method for linear systems with uncertain parameters

A method for stability-robustness analysis based on a quadratic Liapunov function that varies linearly with uncertainty parameters is derived. Linear time-invariant systems with structured uncertainties are discussed. The Liapunov function is optimized numerically to maximize the robustness region in parameter space. Numerical results are given for four examples in which the first-order method is compared to previous Liapunov methods. While the zero-order method is slightly better than the first-order method for one example, the first-order method is clearly superior in the other three (more realistic) examples. The first-order method is especially superior for the active control of flexible structures, where robustness with respect to (1) unmodeled coupling between modeled modes and (2) unmodeled modes is important. For such applications, the first-order method is much better at detecting the increased robustness associated with increased separation between frequencies.

Leal, M. A.↗

Linearized aerodynamic and control law models of the X-29A airplane and comparison with flight data

Flight control system design and analysis for aircraft rely on mathematical models of the vehicle dynamics. In addition to a six degree of freedom nonlinear simulation, the X-29A flight controls group developed a set of programs that calculate linear perturbation models throughout the X-29A flight envelope. The models include the aerodynamics as well as flight control system dynamics and were used for stability, controllability, and handling qualities analysis. These linear models were compared to flight test results to help provide a safe flight envelope expansion. A description is given of the linear models at three flight conditions and two flight control system modes. The models are presented with a level of detail that would allow the reader to reproduce the linear results if desired. Comparison between the response of the linear model and flight measured responses are presented to demonstrate the strengths and weaknesses of the linear models' ability to predict flight dynamics.

Bosworth, John T.↗

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

Research leading to the development of new nonlinear methodologies for the adaptive control and stability analysis of high angle of attack aircraft such as the F-18 is discussed. The emphasis has been on nonlinear adaptive control, but associated model development, system identification, stability analysis, and simulation were studied in some detail as well. Studies indicated that nonlinear adaptive control can outperform linear adaptive control for rapid maneuvers with large changes in angle of attack. Included here are studies on nonlinear model algorithmic controller design and an analysis of nonlinear system stability using robust stability analysis for linear systems.

Mohler, R. R.↗

High-Alpha Research Vehicle Lateral-Directional Control Law Description, Analyses, and Simulation Results

This report contains a description of a lateral-directional control law designed for the NASA High-Alpha Research Vehicle (HARV). The HARV is a F/A-18 aircraft modified to include a research flight computer, spin chute, and thrust-vectoring in the pitch and yaw axes. Two separate design tools, CRAFT and Pseudo Controls, were integrated to synthesize the lateral-directional control law. This report contains a description of the lateral-directional control law, analyses, and nonlinear simulation (batch and piloted) results. Linear analysis results include closed-loop eigenvalues, stability margins, robustness to changes in various plant parameters, and servo-elastic frequency responses. Step time responses from nonlinear batch simulation are presented and compared to design guidelines. Piloted simulation task scenarios, task guidelines, and pilot subjective ratings for the various maneuvers are discussed. Linear analysis shows that the control law meets the stability margin guidelines and is robust to stability and control parameter changes. Nonlinear batch simulation analysis shows the control law exhibits good performance and meets most of the design guidelines over the entire range of angle-of-attack. This control law (designated NASA-1A) was flight tested during the Summer of 1994 at NASA Dryden Flight Research Center.

Davidson, John B.↗

Large displacements and stability analysis of nonlinear propeller structures

The use of linear rigid formats in COSMIC NASTRAN without DMAP procedures for the analysis of nonlinear propeller structures is described. Approaches for updating geometry and applying follower forces for incremental loading are demonstrated. Comparisons are made with COSMIC NASTRAN rigid formats and other independent finite element programs. Specifically, the comparisons include results from the four approaches for updating the geometry using RIGID FORMAT 1, RIGID FORMATS 4 and 13, MARC and MSC/NASTRAN. It is shown that 'user friendly' updating approaches (without DMAPS) can be used to predict the large displacements and instability of these nonlinear structures. These user friendly approaches can be easily implemented by the user and predict conservative results.

Aiello, R. A.↗

Large displacements and stability analysis of nonlinear propeller structures

The use of linear rigid formats in COSMIC NASTRAN without DMAP procedures for the analysis of nonlinear propeller structures is described. Approaches for updating geometry and applying follower forces for incremental loading are demonstrated. The COSMIC NASTRAN rigid formats and other independent finite element programs are compared. The comparisons include results from the four approaches for updating the geometry using RIGID FORMAT 1, RIGID FORMATS 4 and 13, MARC and MSC/NASTRAN. It is shown that user friendly updating approaches can be used to predict the large displacements and instability of these nonlinear structures. The approaches are easily implemented by the user and predict conservative results.

Aiello, R. A.↗

Boundary-layer stability analysis of flight-measured transition data

Correlation of in-flight boundary-layer transition experiments with linear boundary-layer stability theory contributes both to the validation of the numerical methods as well as the analysis of the measured transition process. Transition results obtained in a recent flight experiment, in which the extent of laminar flow and the transition process on the wing of a business-jet fitted with an instrumented glove section were determined, are analyzed. The experiment was conducted at freestream Mach numbers from 0.55 to 0.82, chord Reynolds numbers from 10 to 20 x 10 to the 6th, and leading-edge sweep angles 17 deg to 20 deg. The growth of both Tollmien-Schlichting and crossflow instabilities are predicted using the e exp n method for several flight conditions and the calculated n-factors at transition onset are correlated. Comparison of the measured dominant boundary-layer disturbance frequencies and the predicted unstable frequencies shows fair agreement for several of the flight conditions studied.

Obara, Clifford J.↗

Computer Programs for Calculating and Plotting the Stability Characteristics of a Balloon Tethered in a Wind

Computer programs for calculating the stability characteristics of a balloon tethered in a steady wind are presented. Equilibrium conditions, characteristic roots, and modal ratios are calculated for a range of discrete values of velocity for a fixed tether-line length. Separate programs are used: (1) to calculate longitudinal stability characteristics, (2) to calculate lateral stability characteristics, (3) to plot the characteristic roots versus velocity, (4) to plot the characteristic roots in root-locus form, (5) to plot the longitudinal modes of motion, and (6) to plot the lateral modes for motion. The basic equations, program listings, and the input and output data for sample cases are presented, with a brief discussion of the overall operation and limitations. The programs are based on a linearized, stability-derivative type of analysis, including balloon aerodynamics, apparent mass, buoyancy effects, and static forces which result from the tether line.

Bennett, R. M.↗

Two-dimensional stability of laminar flames

The stability of laminar flames was studied numerically and the dependence of stability on finite rate chemistry with low activation energy and variable thermodynamic and transport properties was addressed. The calculations show that activation energy and details of chemistry play a minor role in altering the linear neutral stability results from asymptotic analysis. Variable specific heat makes a marginal change to the stability; variable transport properties, on the other hand, tend to substantially enhance the stability from a critical wave number of about 0.50 to 0.20. Also, the effects of variable properties tend to nullify the effects of nonunity Lewis number. When the Lewis number of a single species is different from unity, as is true in a hydrogen-air premixed flame, the stability results remain close to that of unity Lewis number.

Mukunda, H. S.↗

Non-linear controls influence functions in an aircraft dynamics simulator

In the development and testing of novel structural and controls concepts, such as morphing aircraft wings, appropriate models are needed for proper system characterization. In most instances, available system models do not provide the required additional degrees of freedom for morphing structures but may be modified to some extent to achieve a compatible system. The objective of this study is to apply wind tunnel data collected for an Unmanned Air Vehicle (UAV), that implements trailing edge morphing, to create a non-linear dynamics simulator, using well defined rigid body equations of motion, where the aircraft stability derivatives change with control deflection. An analysis of this wind tunnel data, using data extraction algorithms, was performed to determine the reference aerodynamic force and moment coefficients for the aircraft. Further, non-linear influence functions were obtained for each of the aircraft s control surfaces, including the sixteen trailing edge flap segments. These non-linear controls influence functions are applied to the aircraft dynamics to produce deflection-dependent aircraft stability derivatives in a non-linear dynamics simulator. Time domain analysis of the aircraft motion, trajectory, and state histories can be performed using these nonlinear dynamics and may be visualized using a 3-dimensional aircraft model. Linear system models can be extracted to facilitate frequency domain analysis of the system and for control law development. The results of this study are useful in similar projects where trailing edge morphing is employed and will be instrumental in the University of Maryland s continuing study of active wing load control.

Guerreiro, Nelson M.↗

Acceleration of convergence of vector sequences

A general approach to the construction of convergence acceleration methods for vector sequence is proposed. Using this approach, one can generate some known methods, such as the minimal polynomial extrapolation, the reduced rank extrapolation, and the topological epsilon algorithm, and also some new ones. Some of the new methods are easier to implement than the known methods and are observed to have similar numerical properties. The convergence analysis of these new methods is carried out, and it is shown that they are especially suitable for accelerating the convergence of vector sequences that are obtained when one solves linear systems of equations iteratively. A stability analysis is also given, and numerical examples are provided. The convergence and stability properties of the topological epsilon algorithm are likewise given.

Sidi, A.↗

Acceleration of convergence of vector sequences

A general approach to the construction of convergence acceleration methods for vector sequence is proposed. Using this approach, one can generate some known methods, such as the minimal polynomial extrapolation, the reduced rank extrapolation, and the topological epsilon algorithm, and also some new ones. Some of the new methods are easier to implement than the known methods and are observed to have similar numerical properties. The convergence analysis of these new methods is carried out, and it is shown that they are especially suitable for accelerating the convergence of vector sequences that are obtained when one solves linear systems of equations iterative. A stability analysis is also given, and numerical examples are provided. The convergence and stability properties of the topological epsilon algorithm are likewise given.

Sidi, A.↗

Transient Growth Analysis of Compressible Boundary Layers with Parabolized Stability Equations

The linear form of parabolized linear stability equations (PSE) is used in a variational approach to extend the previous body of results for the optimal, non-modal disturbance growth in boundary layer flows. This methodology includes the non-parallel effects associated with the spatial development of boundary layer flows. As noted in literature, the optimal initial disturbances correspond to steady counter-rotating stream-wise vortices, which subsequently lead to the formation of stream-wise-elongated structures, i.e., streaks, via a lift-up effect. The parameter space for optimal growth is extended to the hypersonic Mach number regime without any high enthalpy effects, and the effect of wall cooling is studied with particular emphasis on the role of the initial disturbance location and the value of the span-wise wavenumber that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary layer equations, mean flow solutions based on the full Navier-Stokes (NS) equations are used in select cases to help account for the viscous-inviscid interaction near the leading edge of the plate and also for the weak shock wave emanating from that region. These differences in the base flow lead to an increasing reduction with Mach number in the magnitude of optimal growth relative to the predictions based on self-similar mean-flow approximation. Finally, the maximum optimal energy gain for the favorable pressure gradient boundary layer near a planar stagnation point is found to be substantially weaker than that in a zero pressure gradient Blasius boundary layer.

Compressible boundary layer↗

Practical robustness measures in multivariable control system analysis

The robustness of the stability of multivariable linear time invariant feedback control systems with respect to model uncertainty is considered using frequency domain criteria. Available robustness tests are unified under a common framework based on the nature and structure of model errors. These results are derived using a multivariable version of Nyquist's stability theorem in which the minimum singular value of the return difference transfer matrix is shown to be the multivariable generalization of the distance to the critical point on a single input, single output Nyquist diagram. Using the return difference transfer matrix, a very general robustness theorem is presented from which all of the robustness tests dealing with specific model errors may be derived. The robustness tests that explicitly utilized model error structure are able to guarantee feedback system stability in the face of model errors of larger magnitude than those robustness tests that do not. The robustness of linear quadratic Gaussian control systems are analyzed.

Lehtomaki, N. A.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

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

Stabilization of large-scale systems - A spinning flexible spacecraft

Stabilization of high-order constant systems by multilevel control is proposed in the framework of the decomposition-aggregation method for stability analysis of large-scale systems. An 11-order linear model of a spinning flexible spacecraft is stabilized by the method using both the passive and the active stabilization devices. The method can take advantage of the special structural features of the model to provide the best estimate of the stability region for an important parameter which couples the wobble and the spin motion of the spacecraft.

Siljak, D. D.↗