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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An approach to the constrained design of natural laminar flow airfoils

A design method has been developed by which an airfoil with a substantial amount of natural laminar flow can be designed, while maintaining other aerodynamic and geometric constraints. After obtaining the initial airfoil's pressure distribution at the design lift coefficient using an Euler solver coupled with an integml turbulent boundary layer method, the calculations from a laminar boundary layer solver are used by a stability analysis code to obtain estimates of the transition location (using N-Factors) for the starting airfoil. A new design method then calculates a target pressure distribution that will increase the larninar flow toward the desired amounl An airfoil design method is then iteratively used to design an airfoil that possesses that target pressure distribution. The new airfoil's boundary layer stability characteristics are determined, and this iterative process continues until an airfoil is designed that meets the laminar flow requirement and as many of the other constraints as possible.

Green, Bradford Earl↗
NASA NTRS1997

An Approach to the Constrained Design of Natural Laminar Flow Airfoils

A design method has been developed by which an airfoil with a substantial amount of natural laminar flow can be designed, while maintaining other aerodynamic and geometric constraints. After obtaining the initial airfoil's pressure distribution at the design lift coefficient using an Euler solver coupled with an integral turbulent boundary layer method, the calculations from a laminar boundary layer solver are used by a stability analysis code to obtain estimates of the transition location (using N-Factors) for the starting airfoil. A new design method then calculates a target pressure distribution that will increase the laminar flow toward the desired amount. An airfoil design method is then iteratively used to design an airfoil that possesses that target pressure distribution. The new airfoil's boundary layer stability characteristics are determined, and this iterative process continues until an airfoil is designed that meets the laminar flow requirement and as many of the other constraints as possible.

Green, Bradford E.↗
DOE OSTI2024

Wideband Impedance Passivation of MMCs for Suppressing Harmonic Oscillations

Harmonic instability and associated oscillation events have become one of the main concerns in MMC-based HVdc systems. These oscillations can appear in the range from a few hundred hertz to several kilohertz, and the root cause is identified as delay-induced negative damping of MMC impedance appearing at the system resonance frequency. Here, this paper introduces a wideband impedance reshaping method through MMC control to eliminate such negative damping for both grid-following (GFL) and grid-forming (GFM) MMCs, which makes the MMC impedance completely passive from the second harmonic frequency upward, thus preventing all harmonic oscillations. First, simplified impedance models of MMCs are derived for the harmonic stability analysis and the impedance reshaping control design. Next, a passivity-based impedance reshaping method is presented, as well as practical considerations for its implementation. In addition, to maintain the MMC's disturbance ride-through capability, an adaptive activation scheme is developed, which enables the wideband impedance reshaping control only in the presence of harmonic oscillation events. The effectiveness of the proposed method is validated by frequency domain analysis and by electromagnetic transient (EMT) simulations of two typical MMC-based power systems.

17 WIND ENERGY↗
NASA NTRS1986

Analysis and interpretation of the shocklike electrostatic noise observed during the AMPTE solar wind lithium releases

During the AMPTE (Active Magnetospheric Particle Tracer Explorers) solar wind lithium release on September 11, and again on September 20, 1984, an intense burst of electrostatic noise was observed near the upstream edge of the ion cloud. Comparisons with measurements by the IMP-6 and ISEE-1 spacecraft show that the spectrum and overall features of this noise are very similar to electrostatic noise observed at the Earth's bow shock. A stability analysis using realistic parameters shows that the electrostatic noise can be accounted for by an ion beam-plasma instability caused by the solar wind proton beam streaming through the nearly stationary lithium cloud. The growth rate of this instability is largest when the ion density and solar wind proton density are similar, which explains why the noise only occurs near the outer edge of the ion cloud. The similarity to the noise in the Earth's bow shock suggest that a shock may exist in the solar wind plasma flow upstream of the ion cloud. If the noise is associated with a shock, then it must be an electrostatic shock, since the ion cyclotron radii are too small for the existence of a MHD shock. Since the electrostatic instability occurs at phase velocities near the lithium thermal velocity, the electrostatic turbulence may play a role in heating the lithium ions and transferring momentum from the solar wind to the ion cloud.

Gurnett, D. A.↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Control-structure interaction in precision pointing servo loops

The control-structure interaction problem is addressed via stability analysis of a generic linear servo loop model. With the plant described by the rigid body mode and a single elastic mode, structural flexibility is categorized into one of three types: (1) appendage, (2) in-the-loop minimum phase, and (3) in-the-loop nonminimum phase. Closing the loop with proportional-derivative (PD) control action and introducing sensor roll-off dynamics in the feedback path, stability conditions are obtained. Trade studies are conducted with modal frequency, modal participation, modal damping, loop bandwidth, and sensor bandwidth treated as free parameters. Results indicate that appendage modes are most likely to produce instability if they are near the sensor rolloff, whereas in-the-loop modes are most dangerous near the loop bandwidth. The main goal of this paper is to provide a fundamental understanding of the control-structure interaction problem so that it may benefit the design of complex spacecraft and pointing system servo loops. In this framework, the JPL Pathfinder gimbal pointer is considered as an example.

Spanos, John T.↗
NASA NTRS1993

Thermal effects in rapid directional solidification - Linear theory

We study the morphological instability of the planar solid/liquid interface for a unidirectionally-solidified dilute binary mixture. We use a model developed by Boettinger et al. (1985, 1986), Aziz (1982), and Jackson et al. (1980), which allows for nonequilibrium effects on the interface through velocity-dependent segregation and attachment kinetics. Two types of instabilities are found in the linear stability analysis: (1) a cellular instability, and (2) an oscillatory instability driven by disequilibrium effects. Merchant and Davis (1990) characterized these instabilities subject to the frozen-temperature approximation (FTA). The present work relaxes the FTA by including the effects of latent heat and the full temperature distribution. Thermal effects slightly postpone the onset of the cellular instability but dramatically postpone the onset of the oscillatory instability; however, the absolute-stability conditions, at which at high speed the cellular and oscillatory instabilities are suppressed, remain unchanged from the FTA.

Huntley, D. A.↗
NASA NTRS1988

Stabilization and robustness of non-linear unity-feedback system - Factorization approach

The paper is a self-contained discussion of a right factorization approach in the stability analysis of the nonlinear continuous-time or discrete-time, time-invariant or time-varying, well-posed unity-feedback system S1(P, C). It is shown that a well-posed stable feedback system S1(P, C) implies that P and C have right factorizations. In the case where C is stable, P has a normalized right-coprime factorization. The factorization approach is used in stabilization and simultaneous stabilization results.

Desoer, C. A.↗
DOE OSTI2021

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Stability of sequences generated by nonlinear differential systems

A local stability analysis is presented for both the analytic and numerical solutions of the initial value problem for a system of ordinary differential equations. It is shown that, using a proper choice of Liapunov function, a connected region of stable initial values of both the analytic solution and the one-leg k-step numerical solution can be approximated. Attention is given to the example of the two-dimensional problem involving the stability of the longitudinal equations of motion of a gliding jet aircraft.

Brown, R. L.↗
NASA NTRS1992

Digital control algorithms for microgravity isolation systems

New digital control algorithms were developed to achieve the desired acceleration transmissibility function. The attractive electromagnets have been taken as actuators. The relative displacement and the acceleration of the mass were used as feedback signals. Two approaches were developed to find that controller transfer function in Z-domain, which yields the desired transmissibility at each frequency. In the first approach, the controller transfer function is obtained by assuming that the desired transmissibility is known in Z-domain. Since the desired transmissibility H sub d(S) = 1/(tauS+1)(exp 2) is given in S-domain, the first task is to obtain the desired transmissibility in Z-domain. There are three methods to perform this task: bilinear transformation, and backward and forward rectangular rules. The bilinear transformation and backward rectangular rule lead to improper controller transfer functions, which are physically not realizable. The forward rectangular rule does lead to a physically realizable controller. However, this controller is found to be marginally stable because of a pole at Z=1. In order to eliminate this pole, a hybrid control structure is proposed. Here the control input is composed of two parts: analog and digital. The analog input simply represents the velocity (or the integral of acceleration) feedback; and the digital controller which uses only relative displacement signal, is then obtained to achieve the desired closed-loop transfer function. The stability analysis indicates that the controller transfer function is stable for typical values of sampling period. In the second approach, the aforementioned hybrid control structure is again used. First, an analog controller transfer function corresponding to relative displacement feedback is obtained to achieve the transmissibility as 1/(tauS+1)(exp 2). Then the transfer function for the digital control input is obtained by discretizing this analog controller transfer function via bilinear transformation. The stability of the resulting Z-domain closed loop system is analyzed. Also, the frequency response of the Z-domain closed-loop transfer function is determined to evaluate the performance of the control system.

Sinha, Alok↗
NASA NTRS2012

Crossflow Instability on a Wedge-Cone at Mach 3.5

As a follow-on activity to the HyBoLT flight experiment, a six degree half angle wedge-cone model at zero angle of attack has been employed to experimentally and computationally study the boundary layer crossflow instability at Mach 3.5 under low disturbance freestream conditions. Computed meanflow and linear stability analysis results are presented along with corresponding experimental Pitot probe data. Using a model-mounted probe survey apparatus, data acquired to date show a well defined stationary crossflow vortex pattern on the flat wedge surface. This effort paves the way for additional detailed, calibrated flow field measurements of the crossflow instability, both stationary and traveling modes, and transition-to-turbulence under quiet flow conditions as a means of validating existing stability theory and providing a foundation for dynamic flight instrumentation development.

Beeler, George B.↗
NASA NTRSDate not supplied

Controls Algorithms, Models, and Stability Assessments of the RCS and TVC systems on Mars Ascent Vehicle

This paper presents a control systems overview of the Mars Ascent Vehicle (MAV) control design and stability analysis with two separate control methods: Thrust Vector Control (TVC) based control using main engine thrust and Reaction Control System (RCS) via thrusters. A pole-placement TVC controller is proposed to address quickly peaking thrust produced by the MAV solid propulsion engine. The controllerutilizes real-time parameter estimation to calculate control gains using an online pole-placement method. The pole-placement technique allows the controller to maintain constant stability margins throughout the flight under changing parameters and rapidly peaking thrust. A mathematical stability proof via Lyapunov inequality and Nichols method is demonstrated to support the stability of the method. In addition to the TVC control, a phase-plane based RCS control logic is applied to the vehicle for roll control during ascent and all axes during coast phases of flight. The logic implemented is an on-off type logic that applies pre-determined thrusters firings in reference to boundary parameters of pointing and rate error chosen to balance desired pointing with stability. A describing function methodology is applied to the non-linear phase plane loop response to determine linear stability across all flight phases. The controller additionally applies command alteration logic to output signals to account for RCS hardware limitations. The following paper will provide the methodology for controller implementation of the TVC and RCS, and provide a summary of results of the MAV performance applying the aforementioned control design.

Han Woong Bae↗
NASA NTRS1975

SSME turbopump technology improvements via transient rotordynamic analysis

The rotordynamic behavior of the high pressure oxygen turbopump and high pressure fuel pump was analyzed for the Space Shuttle Main Engine. The identification of potential rotordynamic problem areas which might arise during operation of these units prior to their testing was accomplished. Alternative procedures for correcting potential rotordynamic problems should they occur were investigated. An adequate analytic and physical understanding of the turbopump rotordynamics was developed to improve the probability of a correct diagnosis of rotordynamic problems from test data. Transient rotordynamic models were developed for both turbopumps. The transient models model the hydrodynamic forces of the turbopump seals. A linear stability analysis was performed for the turbopump rotordynamics models, which included gyroscopic effects, seal forces, speed-dependent bearing characteristics, and internal rotor damping. Results are presented and discussed.

Childs, D. W.↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Inviscid and viscous flow calculations for the F16XL configuration

The viewgraphs and discussion of the ongoing activity at NASA Langley Research Center (LaRC) in support of Supersonic Laminar FLow Control (SLFC) research, presented at the High-Speed Research Workshop, May 14-16, 1991, are provided. Details of the computation involved in obtaining the meanflow around bodies in high-speed flow and interfacing the results to a stablilty analysis is presented. Particular attention is given to the F-16XL configuration, which is the test-bed for the supersonic LFC experiment. Meanflow solutions for two geometries are discussed. The first one is for the F-16XL wing, with emphasis on the flow near the attachment line and the upper surface. Calculations were done with and without suction. The results were processed using an interface program and fed into a stability analysis program. The second geometry is a scale model of a swept wing leading edge at M = 3.5. Experimental measurements on transition on this model are planned at NASA LaRC. The computations are in support of this effort.

Iyer, Venkit↗
NASA NTRS1987

Development of a generic combustion stability code for liquid propellant rocket engines

The mathematical framework for a combustion stability analysis code is outlined. The goal for the code is to be general enough in problem treatment so that its validity and accuracy extend over a wide range of problem applications and that it lends the convenience for any future model improvement if necessary. An approach for modeling the combustion dynamics is devised to meet both requirements. An open-loop numerical procedure is also developed to mechanistically model various combustion processes for determining the stability parameters.

Fang, J. J.↗