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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Pinhole occulter experiment

Viewgraphs and charts from a briefing summarize the accomplishments, results, conclusions, and recommendations of a feasibility study using the Pinhole Occulter Facility (POF). Accomplishments for 1986 include: (1) improved IPS Gimbal Model; (2) improved Crew Motion Disturbance Model; (3) use of existing shuttle on-orbit simulation to study the effects of orbiter attitude deadband size on POF performance; (4) increased understanding of maximum performance expected from current actuator/sensor set; (5) use of TREETOPS nonlinear time domain program to obtain system dynamics describing the complex multibody flexible structures; (6) use of HONEY-X design tool to design and evaluate multivariable compensator for stability, robustness, and performance; (7) application of state-of-the-art compensator design methodology Linear Quadratic Gaussian/Loop Transfer Recovery (LQG/LTR); and (8) examination of tolerance required on knowledge of the POF boom flexible mode frequencies to insure stability, using structure uncertainty analysis.

Ring, Jeff↗
NASA NTRS1987

Stability of a rigid rotor supported on oil-film journal bearings under dynamic load

Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.

Majumdar, B. C.↗
NASA NTRS1994

Nonlinear instability of accelerating shock waves with application to supernovae

We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.

Luo, Ding↗
DOE OSTI2020

Impacts of Experimentally Obtained Harmonic Spectrums of Residential Appliances on Distribution Feeder: Preprint

This paper focuses on the harmonic load flow study for a set of residential loads and the impacts of the harmonic contents on the distribution feeder. Owing to the increased use of the power electronic based appliances and energy efficient home equipment like the lighting loads, the percentage of non-linear loads has increased in the distribution system which can impact the system performance, loss, and stability. Thus a comprehensive knowledge of their power quality and harmonic analysis is essential to improve the load models, voltage stability assessment, determining possible interactions at harmonic frequencies, protection planning, and the effect of the system impedance. For this study, the household appliances include the lighting, power electronics, resistive, and motor loads with a nominal supply voltage of 120 V single-phase or 240 V split-phase at a nominal frequency of 60 Hz. The harmonic spectrum as obtained from the experimental data is used to solve the harmonic load flow for a detailed secondary network including the distribution transformer extracted from a real distribution feeder's model. Additionally, the impact of the voltage harmonic on a three phase test load is presented through experimental results. The harmonic data of the transformer terminal voltage as obtained from the former load flow study is scaled to generate the source voltage for the three phase configuration of the grid simulator to make the test setup very similar to a real scenario.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1996

A microstructural approach to cytoskeletal mechanics based on tensegrity

Mechanical properties of living cells are commonly described in terms of the laws of continuum mechanics. The purpose of this report is to consider the implications of an alternative approach that emphasizes the discrete nature of stress bearing elements in the cell and is based on the known structural properties of the cytoskeleton. We have noted previously that tensegrity architecture seems to capture essential qualitative features of cytoskeletal shape distortion in adherent cells (Ingber, 1993a; Wang et al., 1993). Here we extend those qualitative notions into a formal microstructural analysis. On the basis of that analysis we attempt to identify unifying principles that might underlie the shape stability of the cytoskeleton. For simplicity, we focus on a tensegrity structure containing six rigid struts interconnected by 24 linearly elastic cables. Cables carry initial tension ("prestress") counterbalanced by compression of struts. Two cases of interconnectedness between cables and struts are considered: one where they are connected by pin-joints, and the other where the cables run through frictionless loops at the junctions. At the molecular level, the pinned structure may represent the case in which different cytoskeletal filaments are cross-linked whereas the looped structure represents the case where they are free to slip past one another. The system is then subjected to uniaxial stretching. Using the principal of virtual work, stretching force vs. extension and structural stiffness vs. stretching force relationships are calculated for different prestresses. The stiffness is found to increase with increasing prestress and, at a given prestress, to increase approximately linearly with increasing stretching force. This behavior is consistent with observations in living endothelial cells exposed to shear stresses (Wang & Ingber, 1994). At a given prestress, the pinned structure is found to be stiffer than the looped one, a result consistent with data on mechanical behavior of isolated, cross-linked and uncross-linked actin networks (Wachsstock et al., 1993). On the basis of our analysis we concluded that architecture and the prestress of the cytoskeleton might be key features that underlie a cell's ability to regulate its shape.

NASA Discipline Cell Biology↗
NASA NTRS2018

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that stationary disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip. Furthermore, the present analysis shows significant nonmodal growth of traveling disturbances that peak within the entropy layer and could also play a role in the transition reversal phenomenon.

Freestream Mach Number↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Stability analysis of flexible wind turbine blades using finite element method

Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.

Kamoulakos, A.↗
DOE OSTI2021

Degradation Comparison of Cyclic and Linear Siloxane Contamination on Solid Oxide Fuel Cells Ni-YSZ Anode

The solid oxide fuel cell (SOFC) nickel-yttria stabilized zirconia (Ni-YSZ) anode degradation due to different types of siloxane contamination is investigated. A cyclic structure siloxane, octamethylcyclotetrasiloxane (D4), and a linear structure siloxane, decamethyltetrasiloxane (L4), are mixed with H 2 +N 2 as the fuel for SOFCs at 750°C. The electrochemical characterization results after stability experiments suggest that the SOFC contaminated with cyclic siloxane, D4, had higher degradation. Pure YSZ pellets with different surface hydroxylation extents were also tested to investigate the D4/L4 adsorption and deposition process. Postmortem SEM/WDS, XRD and Raman analysis all indicate that cyclic siloxane has more deposition than linear siloxane on the anode. Further analysis demonstrates that high adsorption and low desorption rates of cyclic siloxane on YSZ are linked to the degradation. Besides the silicon deposition, SiC and amorphous carbon deposition were also observed from the XRD and Raman analysis.

25 ENERGY STORAGE↗
NASA NTRS2018

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis Under NATO STO AVT-240

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics-based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip.

Numerical analysis↗
DOE OSTI2021

Kinklike structures in an arcsin real scalar dynamics

Highlights: • Kinklike analyses in a real scalar theory with an arcsin dynamic. • First-order formalism implementation and deformation procedure in an arcsin scalar dynamics. In this paper, we analyze kink-like analytical solutions in a real scalar theory with an arcsin dynamics inspired by the arcsin electrodynamics presented in Kruglov (2015). This analysis is done by means of the first-order formalism. This formalism provides a framework where the equations of motion can be simplified by preserving the linear stability of the theory. In this work, the deformation procedure is implemented with the aim of finding exact solutions in systems with generalized dynamics. Along the paper, we explore how the first-order formalism is implemented in the arcsin kinetics and how such a term influences the kink-like solutions. As a part of the result of our paper, we show that the kink-like solutions are similar to the ones obtained in the standard scalar kinetic theory. We also show that the extra parameter, that controls the non-linearities of the model, plays an essential role in the energy densities and stability potentials. These quantities vary according to this parameter. The goals here are to show how the first-order framework is implemented in this arcsin scenario and to present the analytical kink-like solutions that can be found by means of the first-order framework and deformation method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Relativistic stellar stability - Preferred-frame effects

In a previous paper, the PPN (parametrized post-Newtonian) formalism was used to analyze relativistic influences on stellar stability in nearly all metric theories of gravity. That analysis omitted all preferred-frame terms. In this paper, possible preferred-frame effects on stellar stability are examined and no new instabilities are found. Although terms linear in the preferred-frame velocity w (time-odd terms, analogous to viscosity and energy generation) change the shapes of the normal modes, their symmetry properties prevent them from changing the characteristic frequencies. Thus, no new vibrational or secular instabilities can occur. Terms quadratic in w do not change either the shapes of the normal modes or the characteristic frequencies for radial pulsations (except for the effects due to the renormalization of the gravitation constant which does not affect stability). Thus, they have no influence on radial stability. Terms quadratic in w do change both the normal modes and the characteristic frequencies of nonradial pulsations; but in the limit of a neutral mode these changes vanish.

Ni, W.-T.↗
NASA NTRS1995

Static stability of a three-dimensional space truss

In order to deploy large flexible space structures it is necessary to develop support systems that are strong and lightweight. The most recent example of this aerospace design need is vividly evident in the space station solar array assembly. In order to accommodate both weight limitations and strength performance criteria, ABLE Engineering has developed the Folding Articulating Square Truss (FASTMast) support structure. The FASTMast is a space truss/mechanism hybrid that can provide system support while adhering to stringent packaging demands. However, due to its slender nature and anticipated loading, stability characterization is a critical part of the design process. Furthermore, the dire consequences surely to result from a catastrophic instability quickly provide the motivation for careful examination of this problem. The fundamental components of the space station solar array system are the (1) solar array blanket system, (2) FASTMast support structure, and (3) mast canister assembly. The FASTMast once fully deployed from the canister will provide support to the solar array blankets. A unique feature of this structure is that the system responds linearly within a certain range of operating loads and nonlinearly when that range is exceeded. The source of nonlinear behavior in this case is due to a changing stiffness state resulting from an inability of diagonal members to resist applied loads. The principal objective of this study was to establish the failure modes involving instability of the FASTMast structure. Also of great interest during this effort was to establish a reliable analytical approach capable of effectively predicting critical values at which the mast becomes unstable. Due to the dual nature of structural response inherent to this problem, both linear and nonlinear analyses are required to characterize the mast in terms of stability. The approach employed herein is one that can be considered systematic in nature. The analysis begins with one and two-dimensional failure models of the system and its important components. From knowledge gained through preliminary analyses a foundation is developed for three-dimensional analyses of the FASTMast structure. The three-dimensional finite element (FE) analysis presented here involves a FASTMast system one-tenth the size of the actual flight unit. Although this study does not yield failure analysis results that apply directly to the flight article, it does establish a method by which the full-scale mast can be evaluated.

Shaker, John F.↗
NASA NTRS2012

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗
NASA NTRS2014

Impact of Cross-Axis Structural Dynamics on Validation of Linear Models for Space Launch System

A feasibility study was performed to examine the advisability of incorporating a set of Programmed Test Inputs (PTIs) during the Space Launch System (SLS) vehicle flight. The intent of these inputs is to provide validation to the preflight models for control system stability margins, aerodynamics, and structural dynamics. During October 2009, Ares I-X program was successful in carrying out a series of PTI maneuvers which provided a significant amount of valuable data for post-flight analysis. The resulting data comparisons showed excellent agreement with the preflight linear models across the frequency spectrum of interest. However unlike Ares I-X, the structural dynamics associated with the SLS boost phase configuration are far more complex and highly coupled in all three axes. This presents a challenge when implementing this similar system identification technique to SLS. Preliminary simulation results show noticeable mismatches between PTI validation and analytical linear models in the frequency range of the structural dynamics. An alternate approach was examined which demonstrates the potential for better overall characterization of the system frequency response as well as robustness of the control design.

Pei, Jing↗
DOE OSTI2023

The onset distribution of rotating m , n = 2 , 1 tearing modes and its consequences on the stability of high-confinement-mode plasmas in DIII-D

Abstract Analysis of a multi-scenario database of over 13 000 DIII-D H-mode discharges shows that the m , n = 2 , 1 magnetic islands are dominantly pressure gradient driven, stochastically triggered non-linear instabilities at all edge safety factor ( q 95 ) values. The instability onset time closely follows the exponential distribution in intermediate and high q 95 scenarios and is characterized by near constant onset rate ( λ ), in accordance with Poisson-point processes. This implies that the plasmas are operated in marginally stable conditions, characterized by a small threshold for instability growth and variations in the trigger amplitude and/or the stabilizing mechanisms with temporally uniform random distribution in this database. While the majority of the tearing modes occur in the first current-profile relaxation time of the β N flattop, constant λ throughout the β N flattop shows that the tearing onset is insensitive to the evolution of the equilibrium current profile. In low q 95 scenarios, where a large fraction of the plasmas are operated at low torque, λ increases over the course of the β N flattop, showing that these plasmas evolve toward more unstable conditions. The onset rate rapidly increases with β N , while it does not show a clear dependence on the current gradient at the mode rational surface. Overall, these observations support that the majority of the analyzed 2,1 tearing modes are non-linear, neoclassically driven instabilities and classical stability does not play a dominant role in their onset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Stability of flow over axisymmetric bodies with porous suction strips

Linear triple deck, closed form solutions for mean-flow quantities are developed for axisymmetric incompressible flow past a body with porous strips. The solutions account for upstream influence and are linear superpositions of the flow past the body without suction plus the perturbations due to the suction strips. Flow past the suctionless body is calculated using the Transition Analysis Program System, and a simple linear optimization scheme to determine number, spacing, and mass flow rate through the strips on an axisymmetric body is developed using the linear, triple-deck, closed-form solutions. The theory is demonstrated by predicting optimal strip distributions, and the effect of various adverse pressure-gradient situations on stability is studied.

Nayfeh, A. H.↗