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 577 records · Page 32

Simplified Analysis Of Vehicle/Payload Vibrations

Simplified mathematical model of coupled vibrations of vehicle and its payloads used in stability analysis for control purposes. Multiple-payload stability criterion extension of single-payload criterion based on measurements and detailed calculations of fundamental vibrational modes of vehicle and payload. Represented in simplified form by coupled masses and springs vibrating in one dimension, each spring having stiffness that gives rise to previously determined fundamental vibrational frequency associated mass in absence of coupling with other mass. Developed to avoid complexity and cost of full vibrations/stability calculations for each different combination of Space Shuttle payloads, simplified analysis technique also useful in assessing stability in loaded airplanes, ships, trucks, cranes, and conveyor systems.

Yu, Reginald R. J.↗

Pit rim decomposition into multiple quantum dots on surfaces of epitaxial thin films grown on pit-patterned substrates

Here, we report results of dynamical simulations according to an experimentally validated surface morphological evolution model on the formation of regular arrays of quantum dot molecules (QDMs) consisting of 1D arrays of smaller interacting quantum dots (QDs). These QD arrays form along the sides of each pit rim on the surface of a coherently strained thin film epitaxially deposited on a semiconductor substrate, the surface of which consists of a periodic pattern of inverted prismatic pits with rectangular pit openings. We find that this complex QDM pattern results from the decomposition of the pit rim from a “quantum fortress” with four elongated QDs into four 1D arrays of multiple smaller QDs arranged along each side of the pit rim. Systematic parametric analysis indicates that varying the pit opening dimensions and the pit wall inclination directly impacts the number of QDs in the resulting QDM pattern, while varying the pit depth only affects the dimensions of the QDs in the nanostructure pattern. Therefore, the number, arrangement, and sizes of QDs in the resulting pattern of QDMs on the film surface can be engineered precisely by proper tuning of the pit design parameters. Our simulation results are supported by predictions of morphological stability analysis, which explains the pit rim decomposition into multiple QDs as the outcome of a tip-splitting instability and provides a fundamental characterization of the post-instability nanostructure pattern. Our theoretical findings can play a vital role in designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

36 MATERIALS SCIENCE↗

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temporal model of an optically pumped co-doped solid state laser

Currently, research is being conducted on the optical properties of materials associated with the development of solid state lasers in the two micron region. In support of this effort, a mathematical model describing the energy transfer in a holmium laser sensitized with thulium is developed. In this paper, we establish some qualitative properties of the solution of the model, such as non-negativity, boundedness, and integrability. A local stability analysis is then performed from which conditions for asymptotic stability are attained. Finally, we report on our numerical analysis of the system and how it compares with experimental results.

Wangler, T. G.↗

Control Parameter Sensitivity Study for Inverter-Based-Resource Dominated Grids: A Small Signal Stability Approach and Framework

The growing adoption of renewable energy is driving the prevalence of inverter-based resources (IBRs) within power grids. Future power grids will integrate both grid-following IBRs (GFM-IBRs) and grid-forming IBRs (GFL-IBRs) alongside synchronous generators. Therefore, it is crucial to perform stability studies that account for all components and especially control interactions related to IBRs. Extensive research has performed to study the IBR-related stability, however, the sensitivity study of IBRs' control parameters on system stability has not been adequately studied yet, especially from a systematic way. Therefore, this paper conducts a small signal stability analysis for a generic grid with multiple types of resources and develops an analytical framework for assessing the sensitivity of control parameters affecting stability margins. To achieve that, the non-autonomous reduced-order non-linear dynamic model is developed for a generic power system with multiple synchronous generator-based resources (SGBRs), GFM-IBRs, and GFL-IBRs. Based on the analytic model, a systematic framework for parametric sensitivity on systems' asymptotic stability is developed. A parameter sensitivity analysis based on eigenvalue methods is proposed. The impact of the droop controllers of GFM-IBRs, PQ-dispatch and the PLL controller of GFL-IBR on the system asymptotic stability is discussed. This sensitivity study is aiming to provide deep insights on control parameters' impact on system stability, and gives direction for parameter tuning in case of instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Grid Integration of Offshore Wind Power: Standards, Control, Power Quality and Transmission

Offshore wind is expected to be a major player in the global efforts toward decarbonization, leading to exceptional changes in modern power systems. Understanding the impacts and capabilities of the relatively new and uniquely positioned assets in grids with high integration levels of inverter-based resources, however, is lacking, raising concerns about grid reliability, stability, power quality, and resilience, with the absence of updated grid codes to guide the massive deployment of offshore wind. To help fill the gap, this paper presents an overview of the state-of-the-art technologies of offshore wind power grid integration. First, the paper investigates the most current grid requirements for wind power plant integration, based on a harmonized European Network of Transmission System Operators (ENTSO-E) framework and notable international standards, and it illuminates future directions. The paper discusses the wind turbine and wind power plant control strategies, and new control approaches, such as grid-forming control, are presented in detail. The paper reviews recent research on the ancillary services that offshore wind power plants can potentially provide, which, when harmonized, will not only comply with regulations but also improve the value of the asset. The paper explores topics of wind power plant harmonics, reviewing the latest standards in detail and outlining mitigation methods. The paper also presents stability analysis methods for wind power plants, with discussions centered on validity and computational efficiency. Finally, the paper discusses wind power plant transmission solutions, with a focus on high-voltage direct-current topologies and controls.

17 WIND ENERGY↗

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗

A theoretical approach for analyzing the restabilization of wakes

A linear stability analysis provides the basis for an analysis of the effects of prescribed changes upon the global stability characteristics of a bluff body wake. The analysis is applied in particular to the problem of placing a small cylinder in the wake of a larger one in order to suppress the vortex street. A prediction of the alteration in the critical Reynolds number of the system as a function of the location of the second cylinder relative to the first is presented and is shown to compare favorably with experimental results. The effect upon shedding frequency is also considered. The approach relies on the solution of an adjoint eigenvalue problem and an inhomogeneous adjoint boundary value problem. These adjoint solutions are shown to be a key component in the analysis of a variety of wake control strategies.

Hill, D. C.↗

Long-Term Stability of Radio Sources in VLBI Analysis

Positional stability of radio sources is an important requirement for modeling of only one source position for the complete length of VLBI data of presently more than 20 years. The stability of radio sources can be verified by analyzing time series of radio source coordinates. One approach is a statistical test for normal distribution of residuals to the weighted mean for each radio source component of the time series. Systematic phenomena in the time series can thus be detected. Nevertheless, an inspection of rate estimation and weighted root-mean-square (WRMS) variations about the mean is also necessary. On the basis of the time series computed by the BKG group in the frame of the ICRF2 working group, 226 stable radio sources with an axis stability of 10 as could be identified. They include 100 ICRF2 axes-defining sources which are determined independently of the method applied in the ICRF2 working group. 29 stable radio sources with a source structure index of less than 3.0 can also be used to increase the number of 295 ICRF2 defining sources.

Engelhardt, Gerald↗

Human Stability While Exercising on a VIS Device in Zero Gravity

BACKGROUND: A Vibration Isolation and Stabilization (VIS) system is being designed for use with the European Enhanced Exploration Exercise Device (E4D) [1]. The exercise bar is attached to cables that extend and retract through openings in the E4D platform, which the subject stands on while exercising. Since the E4D is mounted on a moving VIS device, the whole system moves under the subject’s feet. While tension in the cables does help stabilize a crew member in 0g by bracing him or her against the moving platform, flexible cables do not offer full support against falling. This raises the question of whether various exercises that are successfully performed in 1g on a stationary device can be performed without losing balance in 0g on a moving device. Quantifying stability conditions and establishing stability requirements for exercise countermeasure systems is a long-standing challenge and this work helped to inform this specific need for integrated E4D/VIS flight project development. METHODS AND RESULTS: A simplified stability analysis can be attempted based on platform accelerations generated by the VIS device simulation. In this type of analysis, the subject is conceived as standing on the platform in 0g while being held down to it by the tension in the cables. If the sideways acceleration of the platform is such that the cables cannot generate a sufficient moment relative to, e.g., the heels or the toes of the subject to overcome the tipping moment from the inertial forces on the accelerated subject’s body, the subject becomes unstable. This approach, which we refer to as the ‘static’ approximation, indicated failure of all the exercises that were considered in the comprehensive VIS analysis of the E4D. It was realized, however, that platform accelerations are not externally imposed but are themselves induced by the motion of the subject’s body during exercise, and a model calculation confirmed that this drastically altered the tipping moment on the subject, with the potential to even switch the direction in which the body would tip over. This made it imperative to consider in a coupled manner both the motion of the exercising subject’s body and the induced VIS platform motion while considering stability. The coupled dynamics requirement was met by the VIS simulation incorporating time-dependent human mass properties atop the platform and driven by the inertial forces from the prescribed subject motion relative to the platform based on motion-capture recorded human trajectories in 1g on a stationary E4D [2]. To analyze simulation results, a stability criterion was also needed. We do not know how to account for the human ‘control system’ that would take visual and vestibular cues as inputs as the platform moves in 0g. We thus chose to follow an approach we had previously used to analyze the dynamic feasibility of performing a task in lunar gravity along the human body trajectory recorded in 1g [3]. In this approach, one analyzes the center of pressure (COP) between the shoes and the ground (or platform) and checks if the COP remains within the convex hull of the footprints on the ground (or platform), commonly referred to as base of support (BOS). Since pressure is vertical and, in the absence of foot restraints, positive, the COP going outside the BOS indicates that the trajectory recorded in 1g is not dynamically feasible in microgravity and/or on the moving platform of the device. The outcome of the COP-based analysis was that some of the key exercises, especially the deadlift and the back squat exercises, were found to be dynamically feasible for at least some number of the exercise cycles, this number increasing with the cable tension. While this type of stability analysis does not account for the likely alteration of the exercise trajectory from its 1g form in space, it does show that there exist at least some realistic trajectories that pass the dynamic feasibility criterion. This adds confidence that, with further adjustment by the subject of the exercise form in 0g on a moving device, a number of exercises can be performed without losing balance.

D Frenkel↗

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron-field instability: Excitation of electron plasma waves by an electric field

Electric fields are commonplace in plasmas and affect transport by driving currents and, in some cases, instabilities. The necessary condition for instability in collisionless plasmas is commonly understood to be described by the Penrose criterion, which quantifies a sufficient relative drift between different populations of particles that must be present for wave amplification via inverse Landau damping. For example, electric fields generate drifts between electrons and ions that can excite the ion-acoustic instability. Here, we use particle-in-cell simulations and linear stability analysis to show that the electric field can drive a fundamentally different type of kinetic instability, named the electron-field instability. This instability excites electron plasma waves with wavelengths ≳30λ De , has a growth rate that is proportional to the electric field strength, and does not require a relative drift between electrons and ions. The Penrose criterion does not apply when accounting for the electric field. Furthermore, the large value of the observed frequency, near the electron plasma frequency, further distinguishes it from the standard ion-acoustic instability, which oscillates near the ion plasma frequency. The ubiquity of macroscopic electric fields in quasineutral plasmas suggests that this instability is possible in a host of systems, including low-temperature and space plasmas. In fact, damping from neutral collisions in such systems is often not enough to completely damp the instability, adding to the robustness of the instability across plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sequence Impedance Measurement of Utility-Scale Wind Turbines and Inverters – Reference Frame, Frequency Coupling, and MIMO/SISO Forms

Sequence impedance responses with or without considering frequency coupling in both MIMO and SISO forms are increasingly used for the stability analysis of three-phase power electronic systems; however, many aspects of sequence impedance measurement are not fully explored. It is not clear if the sequence impedance has a reference frame similar to the dq impedance. If so, the role of the grid voltage angle estimation in aligning the sequence impedance reference frame has not been discussed. Additionally, existing methods for measuring the sequence impedance with frequency coupling are complicated, are not feasible for large wind turbines and inverters, and provide the sequence impedance responses in only either MIMO or SISO form. This paper presents a sequence impedance measurement method that considers the frequency coupling, performs reference frame alignment, demonstrates the impact of the grid voltage angle estimation, and obtains the sequence impedance response in both MIMO and SISO forms. This paper demonstrates the proposed method and practical problems associated with the sequence impedance measurement of utility-scale wind turbines and inverters on a 1.9-MW Type III wind turbine and a 2.2-MVA inverter using an impedance measurement system built around a 7-MW/13.8-kV grid simulator and a 5-MW dynamometer.

17 WIND ENERGY↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

H-mode pedestal improvements with neon injection in DIII-D

Improvements in the H-mode pedestal pressure and global energy confinement were determined on DIII-D with both 3 and 6 MW of neutral beam (NBI) heating with neon injection. The pedestal pressure increased, primarily associated with a significant increase in the pedestal density with a relatively smaller decrease in the pedestal temperature. At the same neon injection rate, a rapid decrease in edge localized mode (ELM) frequency and ~35% increase in plasma stored energy were observed in the 3 MW discharge, while in the 6 MW discharge, the ELM frequency was gradually decreased and the plasma stored energy gradually increased by up to 20%. The measured pedestal widths from both 3 and 6 MW discharges matched the EPED1.0 width scaling, ${\Delta _{{\psi _N}}} = 0.089\beta _{p,ped}^{1/2}$ ,data-, within 20% deviations before and after neon injection. The peeling-ballooning mode (PBM) stability analysis showed that after neon injection, the ballooning stability boundary was increased, while the peeling stability boundary only dropped once the neon level was relatively high. The ballooning stability boundary increase with neon injection is consistent with increase of the ion diamagnetic stabilization (i.e. including the contribution from carbon, neon as well as deuterium ions) which is used in the criterion to determine when the calculated linear growth rates of PBM are unstable. Consequently, the new ELITE computed PBM boundaries were more consistent with the operating points than when only deuterium was included. This in turn suggests that the impurity may help to improve pedestal pressure through affecting ion diamagnetic frequency in the pedestal region, which could positively lead to the ballooning stability boundary extension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Definition and application of longitudinal stability derivatives for elastic airplanes

A set of longitudinal stability derivatives for elastic airplanes is defined from fundamental principles allowing perturbations in forward speed. Application of these derivatives to longitudinal stability analysis by use of approximate expressions for static stability and control parameters as well as the dynamic equations of motion is illustrated. One commonly used alternative formulation for elastic airplanes is shown to yield significant inaccuracies because of inappropriate interpretation of inertial effects.

Kemp, W. B., Jr.↗

Combustion instability with partial length acoustic liners

The development is discussed of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors with partial length acoustic liners, injector face baffles, and nonlinear combustion responses. Three primary objectives included: (1) extension of the Green's Function integral-iteration technique previously developed so that nonlinear three dimensional wave propagation effects could be included in stability predictions; (2) development of an analytical technique and a computer program for the prediction of the damping effects caused by injector face baffles; and (3) inclusion of a nonlinear, wave distortation dependent, combustion response model in the global stability analysis and determination of its effect on combustor stability.

Mitchell, C. E.↗

STARS: An integrated general-purpose finite element structural, aeroelastic, and aeroservoelastic analysis computer program

The details of an integrated general-purpose finite element structural analysis computer program which is also capable of solving complex multidisciplinary problems is presented. Thus, the SOLIDS module of the program possesses an extensive finite element library suitable for modeling most practical problems and is capable of solving statics, vibration, buckling, and dynamic response problems of complex structures, including spinning ones. The aerodynamic module, AERO, enables computation of unsteady aerodynamic forces for both subsonic and supersonic flow for subsequent flutter and divergence analysis of the structure. The associated aeroservoelastic analysis module, ASE, effects aero-structural-control stability analysis yielding frequency responses as well as damping characteristics of the structure. The program is written in standard FORTRAN to run on a wide variety of computers. Extensive graphics, preprocessing, and postprocessing routines are also available pertaining to a number of terminals.

Gupta, Kajal K.↗