Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 433 records · Page 24

A State-Space Model for Stability Boundary Analysis of Grid-Following Voltage Source Converters Considering Grid Conditions

With the growing significance of renewable energy resources and energy storage systems, the number of grid-connected inverters has been rising at an increasingly rapid pace. Generally, these inverters are directly integrated with the distribution network by synchronizing with the grid voltage at the point of common coupling. However, the low grid strength and varying R/X ratios, as the common characteristics of most distribution networks or weak grids, can lead to dynamic interactions that comprise stability and limit the power transfer capacity of grid-connected inverters. To ensure stable operation of the inverters, researchers must determine the stability boundary, described as the maximum power transfer capacity of grid-connected inverters under the premise of maintaining system small-signal stability. For this purpose, we propose to formulate a state-space model of the system in the synchronously rotating dq-frame of reference and perform eigenvalue analysis to determine the stability boundary. With a detailed model of the control structure and parameters of the grid-connected inverters, the stability boundary is identified as a surface with respect to different grid strengths and R/X ratios. Case study results of proposed eigenvalue analysis are compared with those of admittance model-based stability analysis as well as time-domain simulation using a switching model in Matlab/Simulink, validating the effectiveness and accuracy of the proposed eigenvalue analysis for stability boundary identification.

grid-connected inverters↗

Building a Practical Natural Laminar Flow Design Capability

A preliminary natural laminar flow (NLF) design method that has been developed and applied to supersonic and transonic wings with moderate-to-high leading-edge sweeps at flight Reynolds numbers is further extended and evaluated in this paper. The modular design approach uses a knowledge-based design module linked with different flow solvers and boundary layer stability analysis methods to provide a multifidelity capability for NLF analysis and design. An assessment of the effects of different options for stability analysis is included using pressures and geometry from an NLF wing designed for the Common Research Model (CRM). Several extensions to the design module are described, including multiple new approaches to design for controlling attachment line contamination and transition. Finally, a modification to the NLF design algorithm that allows independent control of Tollmien-Schlichting (TS) and cross flow (CF) modes is proposed. A preliminary evaluation of the TS-only option applied to the design of an NLF nacelle for the CRM is performed that includes the use of a low-fidelity stability analysis directly in the design module.

Campbell, Richard L.↗

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.

Lee-Rausch, Elizabeth M.↗

Comprehensive Analysis Aeroelastic Stability Predictions for the Mars Sample Recovery Helicopters and Ingenuity

A rotorcraft comprehensive analysis was used to generate aeroelastic stability predictions for both Ingenuity and the Mars Sample Recovery Helicopters. The analysis was performed multiple times for each rotor system with sequential increases to the level of modeling complexity. Initial cases used an independent blade analysis, trim inflow distribution, no airframe degrees of freedom, and a hover condition (a time-invariant problem). The final model used all blades, dynamic inflow, airframe degrees of freedom, and a Floquet theory approach to the flutter analysis. The approach of sequentially increasing modeling complexity was employed to ensure the causes of any observed instabilities could be identified. For each case, stability results were determined across a range of air densities, tip Mach numbers, and rotor collective pitch settings. For the Sample Recovery Helicopter, a forward flight condition analysis with the highest level of modeling complexity was also performed. The aeroelastic stability predictions showed both the Ingenuity and Sample Recovery Helicopter rotor systems to be stable for all analyzed conditions. Parameter excursions of structural properties were performed to explore conditions resulting in aeroelastic instabilities.

Aeroelastic↗

Unstalled flutter stability predictions and comparisons to test data for a composite prop-fan model

The aeroelastic stability analyses for three graphite/epoxy composite Prop-Fan designs and post-test stability analysis for one of the designs, the SR-3C-X2 are presented. It was shown that Prop-Fan stability can be effectively analyzed using the F203 modal aeroelastic stability analysis developed at Hamilton Standard and that first mode torsion-bending coupling has a direct effect on blade stability. Positive first mode torsion-bending coupling is a destabilizing factor and the minimization of this parameter will increase Prop-Fan stability. It was also shown that Prop-Fan stability analysis using F203 is sensitive to the blade modal data used as input. Calculated blade modal properties varied significantly with the structural analysis used, and these variations are reflected in the F203 calculations.

Turnberg, J. E.↗

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗

Effects of Mean Flow Profiles on Instability of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to investigate the effects of the mean flow profiles on the instability characteristics in the near-injector region of low-density gas jets injected into high-density ambient gas mediums. To achieve this, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round gas jet injected vertically upwards into a high-density ambient gas were performed by assuming three different sets of mean velocity and density profiles. The flow was assumed to be isothermal and locally parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The first set of mean velocity and density profiles assumed were those used by Monkewitz and Sohn for investigating absolute instability in hot jets. The second set of velocity and density profiles assumed for this study were the ones used by Lawson. And the third set of mean profiles included a parabolic velocity profile and a hyperbolic tangent density profile. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results for each set of mean profiles were delineated. Additional information is included in the original extended abstract.

Vedantam, Nanda Kishore↗

A sequential approach to multivariable stability robustness analysis

In sequential loop closure, the importance of evaluating the stability and stability robustness at the intermediate loop closures is well known, yet how the stability and stability robustness evaluated at the intermediate steps contribute to the stability and stability robustness of the overall feedback system must be developed. An analysis of the complete feedback system reveals the multivariable Nyquist contributions from the intermediate loop closures. It is also shown that the results greatly simplify if frequency separation exists between the intermediate loops. The analysis is presented with a two-step loop closure procedure using 'inner' and 'outer' loops which can be generalized to multi-step situations. The control of the longitudinal dynamics of an aircraft is addressed to further clarify and demonstrate the results.

Newman, Brett↗

Efficient and robust phase-split computations in the internal energy, volume, and moles ( UVN ) space

Phase-split computations in an isolated system, which in general may be defined as one where the total internal energy (U), volume (V), and the number of moles (N = N 1 , N 2 ,..., N n ) of the components are fixed at some specified set of values, involves the determination of the temperature (T), pressure (P), plus the amount and composition of the various phases that constitute the system. A simpler, but analogous problem is one where T, P, and N are specified instead. In TPN space, one may first perform a stability analysis to determine whether the system is stable, meaning whether at equilibrium it will split up into multiple phases or remain in single phase. If the single-phase state is unstable, the stability analysis reliably provides a good set of initial guesses in the subsequent phase-split computations. In UVN space, however, we demonstrate that the stability analysis (which is the main subject of our earlier study [1]) may not in general provide good enough initial guesses; we offer alternative strategies for setting up good initial guesses. Furthermore, we show that a combination of successive substitution iteration (SSI) and Newton's method---two iterative methods that are prevalent in the literature in TPN space---facilitates a robust and efficient algorithm for phase-split computations in isolated systems. Finally, this combination has so far not been applied in UVN space.

02 PETROLEUM↗

CFD validation of a supersonic laminar flow control concept

A three-dimensional Navier-Stokes code is used in conjunction with a linear compressible stability analysis code to develop a numerical procedure for prediction of laminar flow transition. The procedure is applied to a modified F-16XL fighter with a laminar flow control glove at supersonic speed. Details of boundary layer stability analysis indicate that, computationally, laminar flow could be realized on the highly swept wing in the absence of the leading edge attachment-line contamination. Effects of the three-dimensionality of the flow were shown to be important in the boundary layer stability analysis. The numerically predicted surface pressures compare favorably with the flight test data.

Woan, Chung-Jin↗

Ares I-X Flight Test Validation of Control Design Tools in the Frequency-Domain

A major motivation of the Ares I-X flight test program was to Design for Data, in order to maximize the usefulness of the data recorded in support of Ares I modeling and validation of design and analysis tools. The Design for Data effort was intended to enable good post-flight characterizations of the flight control system, the vehicle structural dynamics, and also the aerodynamic characteristics of the vehicle. To extract the necessary data from the system during flight, a set of small predetermined Programmed Test Inputs (PTIs) was injected directly into the TVC signal. These PTIs were designed to excite the necessary vehicle dynamics while exhibiting a minimal impact on loads. The method is similar to common approaches in aircraft flight test programs, but with unique launch vehicle challenges due to rapidly changing states, short duration of flight, a tight flight envelope, and an inability to repeat any test. This paper documents the validation effort of the stability analysis tools to the flight data which was performed by comparing the post-flight calculated frequency response of the vehicle to the frequency response calculated by the stability analysis tools used to design and analyze the preflight models during the control design effort. The comparison between flight day frequency response and stability tool analysis for flight of the simulated vehicle shows good agreement and provides a high level of confidence in the stability analysis tools for use in any future program. This is true for both a nominal model as well as for dispersed analysis, which shows that the flight day frequency response is enveloped by the vehicle s preflight uncertainty models.

Johnson, Matthew↗

Advanced grid-forming (GFM) inverter controls, modeling and system impact study for inverter dominated grids

Our main goals in this project are to understand: 1) the stability of the power system with very high penetration of PV or other inverter-based resources (IBR), 2) the interactions among IBRs and traditional generators, and 3) the impact of the GFM inverters on bulk power system. To achieve these goals, we took a holistic research and development approach. First, an impedance-based large system modeling and stability analysis tool was developed to aid the large IBR system small-signal stability analysis. Second, consensus-based control and transient overload ride through were developed to control and coordinate multiple GFM inverters in a large system to maintain transient stability during big grid events. Third, high fidelity electromagnetic transient (EMT) IBR models were developed and integrated into a large reduced WECC model to study the system impact from IBRs. Lastly, extensive validation tests of the GFM control which was developed in this project and implemented in two 1MW commercial solar inverters were performed using different validation platforms, e.g., control hardware-in-the-loop (CHIL) platform, power hardware- in-the-loop (PHIL) platform, and a real PV plant.

14 SOLAR ENERGY↗

A Reversed Impedance-Based Stability Criterion for IBR Grids

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗

Separated flows near the nose of a body of revolution

The solution of the Navier-Stokes equations for the problem of cross-flow separataion about a deforming cylinder was achieved by iteration. It was shown that the separation starts at the rear stagnation point and the point of primary separation moves upstram along the cylinder surface. A general method of linear stability analysis for nonparallel external flows was constructed, which consists of representing the eigenfunctions with complete orthogonal sets and forms characteristic equations with the Galerkin method. The method was applied to the Kovasznay flow which is an exact solution of the Navier-Stokes equation. The results show that when the critical parameter is exceeded, there are only a few isolated unstable eigen-frequencies. Another exact solution is shown to be absolutely and monotonically stable with respect to infinitesimal disturbances of all frequencies. The flow is also globally, asymptotically, and monotonically stable in the mean with respect o three-dimensional disturbances. This result forms the sound foundation of rigorous stability analysis for nonparallel flows, and provides an invaluable test ground for future studies of nonparallel flows in which the basic states do not posses exact solutions. The application of this method to the study of the formation of spiral vorticies near the nose of a rotating body of revolution is underway. The same method will be applied to the stability analysis of reversed flow over a plate with suction.

Lin, S. P.↗

Three-Dimensional Nonequilibrium Boundary Layer Stability & Transition Analysis For Martian Entry

In this work, we apply the boundary layer stability code LASTRAC, modified to include chemical and thermal nonequilibrium capabilities, to a simulated test case based on literature, with three-dimensional effects. LASTRAC provides Parabolized Stability Equations (PSE)as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semiempirical eN methods. Tests of the Mars Science Laboratory (MSL)indicated a range of laminar to transitional flow conditions. In order to reduce the weight of Thermal Protection Systems (TPS), reliable prediction of transition and an understanding of the environmental factors that contribute to transition are necessary. This work is a step towards addressing this need, including the effects of thermochemical nonequilibrium but neglecting radiation, ablation, and ionization effects. Our results indicate that an oblique first mode produces maximum N-Factors on the symmetry plane, however with N-Factors too small to explain transition. Although we are able to show nonequilibrium three-dimensional PSE results inCO2, further study with improved mean flow solutions, accounting for some of the uncertainties inherent in this analysis, or evaluation of alternative transition mechanisms is needed to explain transition in the wind tunnel cases examined here.

H L Kline↗

Resilient Cooperative Voltage Control for Distribution Network with High Penetration Distributed Energy Resources

Here, this paper considers the problem of designing a resilient distributed voltage control algorithm for distribution systems with high penetration of distributed energy resources in the presence of an unknown cyber-attack. The purpose of the attack is to force the system to violate the operating voltage limit by intercepting its communication channels and inserting exogenous signals to perturb and/or modify the information being exchanged. We first review the cooperative voltage control proposed in our previous work and provide a new stability analysis for it. Next, we present a resilient cooperative voltage control algorithm by introducing a virtual system interconnected with the original system such that the voltage can be maintained within the operational limit under unknown attacks. The resiliency of the proposed algorithm is demonstrated via simulations on the IEEE 8500-node system when subjected to an attack which consists of corrupting the data being exchanged in the communication network between two generation units.

cooperative control↗

A multiloop generalization of the circle criterion for stability margin analysis

In order to provide a theoretical tool suited for characterizing the stability margins of multiloop feedback systems, multiloop input-output stability results generalizing the circle stability criterion are considered. Generalized conic sectors with 'centers' and 'radii' determined by linear dynamical operators are employed to specify the stability margins as a frequency dependent convex set of modeling errors (including nonlinearities, gain variations and phase variations) which the system must be able to tolerate in each feedback loop without instability. The resulting stability criterion gives sufficient conditions for closed loop stability in the presence of frequency dependent modeling errors, even when the modeling errors occur simultaneously in all loops. The stability conditions yield an easily interpreted scalar measure of the amount by which a multiloop system exceeds, or falls short of, its stability margin specifications.

Safonov, M. G.↗