Search NASA⌕ Search

SEARCH · Search NASA

Results for “stability analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

Analysis of Wind Turbine Simulation Models: Assessment of Simplified versus Complete Methodologies: Preprint

This paper presents the current status of simplified wind turbine models used for power system stability analysis. This work is based on the ongoing work being developed in IEC 61400-27. This international standard, for which a technical committee was convened in October 2009, is focused on defining generic (also known as simplified) simulation models for both wind turbines and wind power plants. The results of the paper provide an improved understanding of the usability of generic models to conduct power system simulations.

17 WIND ENERGY↗

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD stability of negative triangularity DIII-D plasmas

Abstract Negative triangularity (NT) experiments in DIII-D point to an emergent reactor scenario free of sawteeth, endowed with benign, nondisruptive n = 2 tearing modes, which experience q min ⩾ 1 similar to the positive triangularity hybrid scenario. Plasmas exhibiting this behavior attain β N > 3 , high enough to reconsider long held views of NT stability. Ideal MHD and tearing stability analysis of well-diagnosed equilibrium reconstructions of experimental hybrid-like plasmas predict that among shape parameters, MHD stability limits are only sensitive to average triangularity. Operation is predicted to be possible at β N relevant to upcoming tokamaks and commercial NT reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic stability of unidirectional fiber-reinforced viscoelastic composite plates

This paper deals with a dynamic stability analysis of unidirectional fiber-reinforced composite viscoelastic plates subjected to compressive edge loads. The integrodifferential equations governing the stability problem are obtained by using, in conjunction with a Boltzmann hereditary constitutive law for a three-dimensional viscoelastic medium, a higher-order shear deformation theory of orthotropic plates. Such a theory incorporates transverse shear deformation, transverse normal stress, and rotatory inertia effects. The solution of the stability problem as considered within this paper concerns the determination of the critical in-plane edge loads yielding the asymptotic instability. Numerical applications, based on material properties derived within the framework of Aboudi's micromechanical model, are presented and pertinent conclusions concerning the nature of the loss of stability and the influence of various parameters are outlined.

Chandiramani, N. K.↗

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗

The effects of aeroelastic deformation on the unaugmented stopped-rotor dynamics of an X-Wing aircraft

A new design concept in the development of VTOL aircraft with high forward flight speed capability is that of the X-Wing, a stiff, bearingless helicopter rotor system which can be stopped in flight and the blades used as two forward-swept and two aft-swept wings. Because of the usual configuration in the fixed-wing mode, there is a high potential for aeroelastic divergence or flutter and coupling of blade vibration modes with rigid-body modes. An aeroelastic stability analysis of an X-Wing configuration aircraft was undertaken to determine if these problems could exist. This paper reports on the results of dynamic stability analyses in the lateral and longitudinal directions including the vehicle rigid-body and flexible modes. A static aeroelastic analysis using the normal vibration mode equations of motion was performed to determine the cause of a loss of longitudinal static margin with increasing airspeed. This loss of static margin was found to be due to aeroelastic washin of the forward-swept blades and washout of the aft-swept blades moving the aircraft aerodynamic center forward of the center of gravity. This phenomenon is likely to be generic to X-Wing aircraft.

Gilbert, Michael G.↗

The effects of aeroelastic deformation on the unaugmented stopped-rotor dynamics of an X-Wing aircraft

A new design concept in the development of vertical takeoff and landing aircraft with high forward flight speed capability is that of the X-Wing. The X-Wing is a stiff, bearingless helicopter rotor system which can be stopped in flight and the blades used as two forward-swept wings and two aft-swept wings. Because of the unusual configuration in the fixed-wing mode, there is a high potential for aeroelastic divergence or flutter and coupling of blade vibration modes with rigid-body modes. An aeroelastic stability analysis of an X-Wing configuration aircraft was undertaken to determine if these problems could exist. This paper reports on the results of dynamic stability analyses in the lateral and longitudinal directions including the vehicle rigid-body and flexible modes. A static aeroelastic analysis using the normal vibration mode equations of motion was performed to determine the cause of a loss of longitudinal static margin with increasing airspeed. This loss of static margin was found to be due to aeroelastic 'washin' of the forward-swept blades and 'washout' of the aft-swept blades moving the aircraft aerodynamic center forward of the center of gravity. This phenomenon is likely to be generic to X-Wing aircraft.

Gilbert, Michael G.↗

Conservation laws, bilinear forms and solitons for a fifth-order nonlinear Schrödinger equation for the attosecond pulses in an optical fiber

Under investigation in this paper is a fifth-order nonlinear Schrödinger equation, which describes the propagation of attosecond pulses in an optical fiber. Based on the Lax pair, infinitely-many conservation laws are derived. With the aid of auxiliary functions, bilinear forms, one-, two- and three-soliton solutions in analytic forms are generated via the Hirota method and symbolic computation. Soliton velocity varies linearly with the coefficients of the high-order terms. Head-on interaction between the bidirectional two solitons and overtaking interaction between the unidirectional two solitons as well as the bound state are depicted. For the interactions among the three solitons, two head-on and one overtaking interactions, three overtaking interactions, an interaction between a bound state and a single soliton and the bound state are displayed. Graphical analysis shows that the interactions between the two solitons are elastic, and interactions among the three solitons are pairwise elastic. Stability analysis yields the modulation instability condition for the soliton solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

Distributed non-negative matrix factorization with determination of the number of latent features

The holistic analysis and understanding of the latent (that is, not directly observable) variables and patterns buried in large datasets is crucial for data-driven science, decision making and emergency response. Such exploratory analyses require devising unsupervised learning methods for data mining and extraction of the latent features, and non-negative matrix factorization (NMF) is one of the prominent such methods. NMF is based on compute-intense non-convex constrained minimization, which, for large datasets requires fast and distributed algorithms. However, current parallel implementations of NMF fail to estimate the number of latent features. In practice, identifying these features is both difficult and significant for pattern recognition and latent feature analysis, especially for large dense matrices. Here, we introduce a distributed NMF algorithm coupled with distributed custom clustering followed by a stability analysis on dense data, which we call DnMFk, to determine the number of latent variables. The results on synthetic data and the classical Swimmer data set demonstrate the accuracy of model determination while scaling nearly linearly across multiple processors for large data. Further, we employ DnMFk to determine the number of hidden features from a terabyte matrix.

97 MATHEMATICS AND COMPUTING↗

PFPO plasma scenarios for exploration of long pulse operation in ITER

Long Pulse Scenarios (LPS) in ITER foreseen during the Pre-Fusion Power Operation (PFPO) phase of the ITER Research Plan (IRP) are assessed using 1.5D transport simulations within the ASTRA framework. Such assessment is required to predict the operational space for LPS operation in PFPO, as well as to evaluate which physics processes for LPS operation during Fusion Power Operation (FPO) could be studied during PFPO. An important aspect in the development of LPSs in PFPO is to minimize lifetime consumption of the Central Solenoid (CS) for these scenarios. The maximum pulse length achievable for LPSs in PFPO with no consumption of CS lifetime (currents in CS coils $⩽$30 kA per turn) has been assessed for a range of heating schemes and heating mixes, confinement regimes (L-mode and H-mode) and for helium and hydrogen plasmas. The operational space of LPS and pulse length has been explored through density scans with the Heating and Current Drive mix required for the FPO Q $⩾$ 5 steady-state plasma scenario (namely Neutral Beam Injection and Electron Cyclotron Heating) including acceptable shine through losses on the first wall for both helium and hydrogen plasmas. Fast particle physics aspects that are common between FPO plasmas and LPS PFPO H-mode plasmas at low densities are studied including MHD stability analysis with the KINX code and non-perturbative critical gradient model based on high-n Toroidal Alfven Eigenmodes (TAE) stability kinetic ballooning code HINST calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrochemical Decomposition of Primary Alcohol Groups in Deep Eutectic Solvents

This work investigates the electrochemical decomposition of an ethylene glycol/sodium chloride electrolyte and a deep eutectic solvent (DES) comprised of ethylene glycol/choline chloride (ethaline). Understanding these reactions can help to better understand the safety risks of overcharge when these electrolytes are used in a battery and to guide the design of DES electrolytes for improved electrochemical stability. Analysis of the decomposition products was performed using a combination of coulometry, gas chromatography, mass spectrometry, and Fourier transform infrared spectroscopy. It was found that even at relatively low (mM) concentrations of water, a stoichiometric rate of hydrogen is produced under cathodic conditions in both the ethaline and EG/NaCl electrolytes. This suggest an electrochemical hydrolysis cathodic decomposition reaction that produces both hydrogen gas and hydroxyl ions. This reaction was verified in both electrolytes by also identifying the liquid phase products. Only liquid phase products were observed under anodic conditions. A suspected alcohol oxidation to aldehyde anodic decomposition reaction was confirmed through a similar analysis of products. A tertiary alcohol only DES was investigated to eliminate these decomposition paths and was found to have increased stability window compared to the primary alcohol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Turbulent Spot Pressure Fluctuation Wave Packet Model

Wave packet analysis provides a connection between linear small disturbance theory and subsequent nonlinear turbulent spot flow behavior. The traditional association between linear stability analysis and nonlinear wave form is developed via the method of stationary phase whereby asymptotic (simplified) mean flow solutions are used to estimate dispersion behavior and stationary phase approximation are used to invert the associated Fourier transform. The resulting process typically requires nonlinear algebraic equations inversions that can be best performed numerically, which partially mitigates the value of the approximation as compared to a more complete, e.g. DNS or linear/nonlinear adjoint methods. To obtain a simpler, closed-form analytical result, the complete packet solution is modeled via approximate amplitude (linear convected kinematic wave initial value problem) and local sinusoidal (wave equation) expressions. Significantly, the initial value for the kinematic wave transport expression follows from a separable variable coefficient approximation to the linearized pressure fluctuation Poisson expression. The resulting amplitude solution, while approximate in nature, nonetheless, appears to mimic many of the global features, e.g. transitional flow intermittency and pressure fluctuation magnitude behavior. A low wave number wave packet models also recover meaningful auto-correlation and low frequency spectral behaviors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Utilizing Flight Data to Update Aeroelastic Stability Estimates

Stability analysis of high performance aircraft must account for errors in the system model. A method for computing flutter margins that incorporates flight data has been developed using robust stability theory. This paper considers applying this method to update flutter margins during a post-flight or on-line analysis. Areas of modeling uncertainty that arise when using flight data with this method are investigated. The amount of conservatism in the resulting flutter margins depends on the flight data sets used to update the model. Post-flight updates of flutter margins for an F/A-18 are presented along with a simulation of on-line updates during a flight test.

Lind, Rick↗

Role of anisotropy in determining stability of electrodeposition at solid-solid interfaces

Here, we investigate the stability of electrodeposition at solid-solid interfaces for materials exhibiting an anisotropic mechanical response. The stability of electrodeposition or resistance to the formation of dendrites is studied within a linear stability analysis. The deformation and stress equations are solved using the Stroh formalism and faithfully recover the boundary conditions at the interface. The stability parameter is used to quantify the stability of different solid-solid interfaces incorporating the full anisotropy of the elastic tensor of the two materials. Our results show a high degree of variability in the stability parameter depending on the crystallographic orientation of the solids in contact, and point to opportunities for exploiting this effect in developing Li metal anodes.

25 ENERGY STORAGE↗