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 307 records · Page 17

Structure of the Compressible Reacting Mixing Layer: A Linear Stability Analysis

Previous investigations have demonstrated that a mixing layer in compressible reacting conditions can develop two instability modes in addition to the more common central mode, which exists unaccompanied in incompressible non-reacting flows. These 'outer' modes were shown in the investigations of Planche and Reynolds to have a detrimental effect on mixing and combustion efficiency and therefore their presence has important consequences for applications of supersonic combustion. Results from a parametric study of the compressible reacting mixing layer's regime space using a linear stability technique will be discussed. The purpose is to develop a more complete understanding behind the effects of compressibility, heat release and the density, equivalence and velocity ratios on the amplification rate and structure of each instability mode.

Day, M. J.↗

Dynamic Stability Analysis Using Free-Flight CFD

Blunt body probes are the primary design for atmospheric entry. The blunted capsule shape, while desirable for the decreased forebody heating, typically experiences dynamic instability at low supersonic and transonic flight. At moderate to low supersonic conditions, the unsteady wake begins to adversely interact with the vehicle in free-flight resulting in attitude oscillations which can grow in time and result in tumbling. Understanding vehicle dynamics is thus important when considering vehicle entry design and guidance navigation and control (GN&C) considerations. Historically, dynamic stability has been carried out exclusively using experimental techniques. The wide range of experimental approaches provide reasonable coverage in Mach and dynamic pressure space for a given trajectory, however achieving flight similitude is often difficult or impossible. Advancements in Computational Fluid Dynamics (CFD) techniques allow for dynamic motion of the vehicle surface during the simulation with either forced/prescribed motion or free motion in response to fluid forces. Several examples of CFD software with the ability to simulate fluid dynamics and vehicle motion are available in literature. The Free-Flight CFD (FF-CFD) capability developed by NASA Ames within the US3D flow solver has shown extensive validation and verification (V&V) efforts from moderate to low supersonic flow across a range of vehicle architectures. The solver is able to simulate full rotational and translation motion resulting in full 6 degrees-of-freedom (DoF) motion. The validation and verification efforts the solver's ability to predict the motion of a ballistic range model as compared to experimental results. Similar to other computational capabilities in literature, FF-CFD relied on heritage data reduction methodologies and processes for reducing dynamic data into aerodynamic coefficients intended for use in design and trajectory codes. These heritage approaches, initially developed with the assumptions (and restrictions) of ground test facilities, resulted in the significant culling of the full 6-DoF dynamic data set to discrete values which are then fit using an aerodynamic model. Recent development of simulation and data reduction methodologies has resulted in an improved approach to FF-CFD simulations of vehicle dynamics which leverage the control of dynamic constraints within a simulation set-up. The presentation of this work will cover the constrained dynamic simulation approach for deriving dynamic coefficients for a Genesis capsule. The methodology will be outlined and a discussion of implications in motion restriction will be discussed. Finally, reconstructed vehicle trajectories using FF-CFD derived aerodynamic coefficients will be presented.

Joseph M Brock↗

Computer program system for dynamic simulation and stability analysis of passive and actively controlled spacecraft. Volume 1. Theory

A theoretical development and associated digital computer program system is presented. The dynamic system (spacecraft) is modeled as an assembly of rigid and/or flexible bodies not necessarily in a topological tree configuration. The computer program system may be used to investigate total system dynamic characteristics including interaction effects between rigid and/or flexible bodies, control systems, and a wide range of environmental loadings. Additionally, the program system may be used for design of attitude control systems and for evaluation of total dynamic system performance including time domain response and frequency domain stability analyses. Volume 1 presents the theoretical developments including a description of the physical system, the equations of dynamic equilibrium, discussion of kinematics and system topology, a complete treatment of momentum wheel coupling, and a discussion of gravity gradient and environmental effects. Volume 2, is a program users' guide and includes a description of the overall digital program code, individual subroutines and a description of required program input and generated program output. Volume 3 presents the results of selected demonstration problems that illustrate all program system capabilities.

Bodley, C. S.↗

Correlation of boundary layer stability analysis with flight transition data

Recently, NASA completed a boundary-layer transition flight test on an F-14 aircraft which has variable-sweep capability. Transition data were acquired for a wide variety of sweep angles, pressure distributions, Mach numbers, and Reynolds numbers. In this paper, the F-14 flight test is briefly described and N-factor correlations with measured transition locations are presented for one of two gloves flown on the F-14 wing in the flight program; a thin foam and fiberglass glove which provided a smooth sailplane finish on the basic F-14, modified NACA 6-series airfoil. For these correlations, an improved linear boundary-layer stability theory was utilized that accounts for compressibility and surface and streamline curvature effects for the flow past swept wings.

Collier, F. S., Jr.↗

Co-Simulation of PSS/E, OpenDSS, and PSCAD for Power Systems Stability Analysis With Inverter-Based Resources: Preprint

The increasing penetration of inverter-based resources (IBRs) is reshaping the dynamic behavior of power systems. IEEE standard 1547-2018 suggests that distributed energy resources (DERs) should provide grid services such as voltage and frequency supports. On the other hand, dynamic events caused by IBRs such as sub-synchronous oscillation have been reported. These developments necessitate improvement in the current modeling capabilities to better understand the interdependencies within power systems. These include interactions between transmission and distribution systems, among IBRs themselves, and between IBRs and conventional resources. In this paper, we present a co-simulation model integrating PSS/E, OpenDSS, and PSCAD to analyze IBR impacts on the stability of transmission and distribution systems. A key challenge in developing a co-simulation model is ensuring interoperability among different simulators (interfacing and data flow) while maintaining accurate results. Using the developed model, we simulate the impact of IBRs on power systems in two test cases: 1) fault ride through (FRT) capability during a generation trip contingency; 2) IBR-induced sub-synchronous oscillation. The results show the effectiveness of the co-simulation model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Stability Analysis for a Hydrodynamic Three-Wave Journal Bearing

The influence of the wave amplitude and oil supply pressure on the dynamic behavior of a hydrodynamic three-wave journal bearing is presented. Both, a transient and a small perturbation technique, were used to predict the threshold to fractional frequency whirl (FFW). In addition, the behavior of the rotor after FFW appeared was determined from the transient analysis. The turbulent effects were also included in the computations. Bearings having a diameter of 30 mm, a length of 27.5 mm, and a clearance of 35 microns were analyzed. Numerical results were compared to experimental results obtained at the NASA GRC. Numerical and experimental results showed that the above-mentioned wave bearing with a wave amplitude ratio of 0.305 operates stably at rotational speeds up to 60,000 rpm, regardless of the oil supply pressure. For smaller wave amplitude ratios, a threshold of stability was found. It was observed that the threshold of stability for lower wave amplitude strongly depends on the oil supply pressure and on the wave amplitude. When the FFW occurs, the journal center maintains its trajectory inside the bearing clearance and therefore the rotor can be run safely without damaging the bearing surfaces.

Ene, Nicoleta M.↗

Stability Analysis for HIFiRE Experiments

The HIFiRE-1 flight experiment provided a valuable database pertaining to boundary layer transition over a 7-degree half-angle, circular cone model from supersonic to hypersonic Mach numbers, and a range of Reynolds numbers and angles of attack. This paper reports selected findings from the ongoing computational analysis of the measured in-flight transition behavior. Transition during the ascent phase at nearly zero degree angle of attack is dominated by second mode instabilities except in the vicinity of the cone meridian where a roughness element was placed midway along the length of the cone. The growth of first mode instabilities is found to be weak at all trajectory points analyzed from the ascent phase. For times less than approximately 18.5 seconds into the flight, the peak amplification ratio for second mode disturbances is sufficiently small because of the lower Mach numbers at earlier times, so that the transition behavior inferred from the measurements is attributed to an unknown physical mechanism, potentially related to step discontinuities in surface height near the locations of a change in the surface material. Based on the time histories of temperature and/or heat flux at transducer locations within the aft portion of the cone, the onset of transition correlated with a linear N-factor, based on parabolized stability equations, of approximately 13.5. Due to the large angles of attack during the re-entry phase, crossflow instability may play a significant role in transition. Computations also indicate the presence of pronounced crossflow separation over a significant portion of the trajectory segment that is relevant to transition analysis. The transition behavior during this re-entry segment of HIFiRE-1 flight shares some common features with the predicted transition front along the elliptic cone shaped HIFiRE-5 flight article, which was designed to provide hypersonic transition data for a fully 3D geometric configuration. To compare and contrast the crossflow dominated transition over the HIFiRE-1 and HIFiRE-5 configurations, this paper also analyzes boundary layer instabilities over a subscale model of the HIFiRE-5 flight configuration that was tested in the Mach 6 quiet tunnel facility at Purdue University.

Li, Fei↗

Orbital stability analysis of hypothetical Earth-mass and Luna-mass moons in the Sagarmatha (HD 100777) star system

Abstract Six of the solar system planets have 150 confirmed moons (Earth—1, Mars—2, Jupiter—53, Saturn—53, Uranus—27, Neptune—14) and seven of them (Ganymede, Titan, Callisto, Io, the Moon, Europa, Triton) have masses >0.001 $M_\oplus$. However, no exomoons have yet been discovered despite the successful detection of ∼5000 exoplanets. We can infer, based on the solar system planets, that these exoplanets are capable of hosting one or more exomoons. In this paper, we study the possible existence of hypothetical Earth-mass and Luna-mass moons orbiting the Jupiter-mass planet, Laligurans (HD 100777b) in the Sagarmatha (HD 100777) star system by means of orbital stability. We apply long-term orbital integrations and the MEGNO (mean exponential growth of nearby orbits) chaos indicator to study the orbital stability of the moons and predict a phase-space region comprising periodic, chaotic, and unstable orbits. The phase spaces primarily constitute the moon's semimajor axis, which extends from the host planet's Roche radius to the Hill radius, and full range of eccentricity. Specific points are picked from three different regions of the MEGNO map and run as single-orbit integration for up to 10 billion periods of the innermost orbit. Furthermore, the lifetime and maximum eccentricity maps are generated from the direct integration to inspect the stable and unstable orbital configurations. The analyses of these maps, with the aid of time-series plots, show that both moons maintain stable orbits in the low-eccentricity regime and semimajor axis between the Roche limit and 28.4% of the Hill radius of the planet.

Gautam, Aayush (ORCID:0000000247257589)↗

Unstructured Grid Euler Method Assessment for Longitudinal and Lateral/Directional Stability Analysis of the HSR Reference H Configuration at Transonic Speeds

Transonic Euler computations, based on unstructured grid methodology, are performed for a proposed High Speed Civil Transport (HSCT) configuration, designated as the Reference H configuration within the High Speed Research (HSR) Program. The predicted results are correlated with appropriate experimental wind-tunnel data for the baseline configuration with and without control surface deflections for a range of angle of attack at M(sub infinity) = 0.95. Good correlations between the predictions and measured data have been obtained for the longitudinal aerodynamic characteristics of the baseline configuration. The incremental effects in the longitudinal aerodynamic characteristics due to horizontal rail deflections as well as wing leading-edge and trailing-edge flap deflections have also been predicted reasonably well. Computational results and correlations with data are also presented for the lateral and directional stability characteristics for a range of angle of attack at a constant sideslip angle as well as a range of sideslip angles at a constant angle of attack. In addition, the results are presented to assess the computational method performance and convergence characteristics.

Ghaffari, Farhad↗

Fully Coupled Aeroelastic Stability Analysis of Adaptive Shape Memory Alloy Structural Technologies for Airframe Noise Reduction

The objective of this work is the development of computational models and analysis of the coupled fluid-structure response of a slat gap filler (SGF) noise treatment applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The representative airframe chosen for this work is NASA’s High-lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs must be assessed for static aeroelastic deflection as well as dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL CSD) tools. The technical approach consists of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design is assessed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability.

Fluid structure interaction↗

Stability analysis of multiple-robot control systems

In a space telerobotic service scenario, cooperative motion and force control of multiple robot arms are of fundamental importance. Three paradigms to study this problem are proposed. They are distinguished by the set of variables used for control design. They are joint torques, arm tip force vectors, and an accelerated generalized coordinate set. Control issues related to each case are discussed. The latter two choices require complete model information, which presents practical modeling, computational, and robustness problems. Therefore, focus is on the joint torque control case to develop relatively model independent motion and internal force control laws. The rigid body assumption allows the motion and force control problems to be independently addressed. By using an energy motivated Lyapunov function, a simple proportional derivative plus gravity compensation type of motion control law is always shown to be stabilizing. The asymptotic convergence of the tracing error to zero requires the use of a generalized coordinate with the contact constraints taken into account. If a non-generalized coordinate is used, only convergence to a steady state manifold can be concluded. For the force control, both feedforward and feedback schemes are analyzed. The feedback control, if proper care has been taken, exhibits better robustness and transient performance.

Wen, John T.↗

The Radiometric Scaling of the MODIS and VIIRS Imagers to a Common Reference and Stability Analysis for the Next Edition of CERES Products

The NASA CERES project has provided the climate quality observed TOA and computed surface fluxes to the scientific community. CERES uses MODIS and VIIRS imagers to retrieve cloud properties needed to convert Terra, Aqua, NPP, and NOAA-20 CERES footprint radiance observations into fluxes. The imagers are also used to radiometrically scale the geostationary sensors (GEO) radiances to the imager calibration reference to ensure that the GEO derived cloud properties and broadband TOA fluxes are consistent in both space and time. Both the imager and GEO retrieved cloud properties are used to compute the surface fluxes. The CERES imager and GEO calibration team uses ray-matched radiance pairs to radiometrically scale the SNPP and NOAA-20 VIIRS sensors to the Aqua-MODIS calibration reference. The radiometric scaling is further validated using geostationary imagers as transfer radiometers. The team will rely primarily on deep convective clouds to monitor the imager channel calibration stability and to correct for short term calibration drifts. The DSCOVR satellite was launched on February 25, 2015 and orbits around the L1 Lagrange point directly between the Earth and the sun. The EPIC sensor contains no onboard calibration systems. However, multiple inter-calibration studies have shown that the EPIC imager is very stable in time. The excellent radiometric stability of EPIC allows the use of EPIC observations as a stable reference to validate both the short term drift corrections of the imagers, as well as to validate the radiometric scaling factors between them. Examples of the imager relative calibration using EPIC before and after radiometric scaling will be shown along with the results from the use of invariant targets to remove imager calibration drifts.

DSCOVR-EPIC↗

Stability Analysis of Streaks Induced by Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗

Stability Analysis of Streaks Induced By Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗

Dynamic Stability Analysis of Blunt Body Entry Vehicles Using Time-Lagged Aftbody Pitching Moments

This analysis defines an analytic model for the pitching motion of blunt bodies during atmospheric entry. The proposed model is independent of the pitch damping sum coefficient present in the standard formulation of the equations of motion describing pitch oscillations of a decelerating blunt body, instead using the principle of a time-lagged aftbody moment as the forcing function for oscillation divergence. Four parameters, all with intuitive physical relevance, are introduced to fully define the aftbody moment and the associated time delay. It is shown that the dynamic oscillation responses typical to blunt bodies can be produced using hysteresis of the aftbody moment in place of the pitch damping coefficient. The approach used in this investigation is shown to be useful in understanding the governing physical mechanisms for blunt body dynamic stability and in guiding vehicle and mission design requirements. A validation case study using simulated ballistic range test data is conducted. From this, parameter identification is carried out through the use of a least squares optimizing routine. Results show good agreement with the limited existing literature for the parameters identified, suggesting that the model proposed could be validated by an experimental ballistic range test series. The trajectories produced by the identified parameters were found to match closely those from the MER ballistic range tests for a wide array of initial conditions and can be identified with a reasonable number of ballistic range shots and computational effort.

Supersonic↗

Stability Analysis of Parameter Varying Genetic Toggle Switches Using Koopman Operators

The genetic toggle switch is a well known model in synthetic biology that represents the dynamic interactions between two genes that repress each other. The mathematical models for the genetic toggle switch that currently exist have been useful in describing circuit dynamics in rapidly dividing cells, assuming fixed or time-invariant kinetic rates. There is a growing interest in being able to model and extend synthetic biological function for growth conditions such as stationary phase or during nutrient starvation. As cells transition from one growth phase to another, kinetic rates become time-varying parameters. In this paper, we propose a novel class of parameter varying nonlinear models that can be used to describe the dynamics of genetic circuits, including the toggle switch, as they transition from different phases of growth. We show that there exists unique solutions for this class of systems, as well as for a class of systems that incorporates the microbial phenomena of quorum sensing. Further, we show that the domain of these systems, which is the positive orthant, is positively invariant. We also showcase a theoretical control strategy for these systems that would grant asymptotic monostability of a desired fixed point. We then take the general form of these systems and analyze their stability properties through the framework of time-varying Koopman operator theory. A necessary condition for asymptotic stability is also provided as well as a sufficient condition for instability. A Koopman control strategy for the system is also proposed, as well as an analogous discrete time-varying Koopman framework for applications with regularly sampled measurements.

Harrison, Jamiree (ORCID:0000000265351579)↗

Stability Analysis of Intertank Formed Skin/Stringer Compression Panel with Simulated Damage

The External Tank (ET) is a component of the Space Shuttle launch vehicle that contains fuel and oxidizer. During launch, the ET supplies the space shuttle main engines with liquid hydrogen and liquid oxygen. In addition to supplying fuel and oxidizer, it is the backbone structural component of the Space Shuttle. It is comprised of a liquid hydrogen (LH2) tank and a liquid oxygen (LOX) tank, which are separated by an Intertank. The Intertank is a stringer-stiffened cylindrical structure with hat-section stringers that are roll formed from aluminum-lithium alloy Al-2090. Cracks in the Intertank stringers of the STS-133 ET were noticed after a November 5, 2010 launch attempt. The cracks were approximately nine inches long and occurred on the forward end of the Intertank (near the LOX tank), along the fastener line, and were believed to have occurred while loading the ET with the cryogenic propellants. These cracks generated questions about the structural integrity of the Intertank. In order to determine the structural capability of the Intertank with varying degrees of damage, a finite element model (FEM) simulating a 1995 compression panel test was analyzed and correlated to test data. Varying degrees of damage were simulated in the FEM, and non-linear stability analyses were performed. The high degree of similarity between the compression panel and the Intertank provided confidence that the ET Intertank would have similar capabilities.

Harper, David W.↗