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 1,405 records · Page 78

Ares-I-X Stability and Control Flight Test: Analysis and Plans

The flight test of the Ares I-X vehicle provides a unique opportunity to reduce risk of the design of the Ares I vehicle and test out design, math modeling, and analysis methods. One of the key features of the Ares I design is the significant static aerodynamic instability coupled with the relatively flexible vehicle - potentially resulting in a challenging controls problem to provide adequate flight path performance while also providing adequate structural mode damping and preventing adverse control coupling to the flexible structural modes. Another challenge is to obtain enough data from the single flight to be able to conduct analysis showing the effectiveness of the controls solutions and have data to inform design decisions for Ares I. This paper will outline the modeling approaches and control system design to conduct this flight test, and also the system identification techniques developed to extract key information such as control system performance (gain/phase margins, for example), structural dynamics responses, and aerodynamic model estimations.

Brandon, Jay M.↗

Thermal Design of the Instrument for the Transiting Exoplanet Survey Satellite

TESS observatory is a two year NASA Explorer mission which will use a set of four cameras to discover exoplanets. It will be placed in a high-earth orbit with a period of 13.7 days and will be unaffected by temperature disturbances caused by environmental heating from the Earth. The cameras use their stray-light baffles to passively cool the cameras and in turn the CCD's in order to maintain operational temperatures. The design has been well thought out and analyzed to maximize temperature stability. The analysis shows that the design keeps the cameras and their components within their temperature ranges which will help make it a successful mission. It will also meet its survival requirement of sustaining exposure to a five hour eclipse. Official validation and verification planning is underway and will be performed as the system is built up. It is slated for launch in 2017.

Thermal Engineering↗

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗

Structural Architectures for Self-Erecting Lunar Towers

The load carrying performance for deployable lunar tower architectures utilizing thin-ply composite booms are parametrically evaluated in this paper. These self-erecting towers support a communications and sensor platform on the lunar surface. The primary design consists of the Corrugated Rollable Tubular Boom (COROTUB) acting as the primary mast of the tower which supports the payload tip mass. Guy wires connecting the boom tip to spreader bars attached at the base deployer serve to correct any lateral eccentricities and provide dimensional stability. The analysis compares the maximum tip mass that this primary design can carry relative to a reference design that does not contain any guy wires or spreader bars. The results provide insights into whether the additional system mass and complexity associated with the cable-stayed design is worth the increased payload mass that can be supported. By setting the design limiting condition to a nominal maximum allowable deflection, the analysis reveals that primary design outperforms reference design for the majority of the design space. This outcome is found to be the case when the spreader bar length relative to the tower height exceeds7 a critical value. The analysis is extended to lunar towers with Collapsible Tubular Mast (CTM) booms, and similar trends are observed.

Jacob G Daye↗

Structural Architectures for Self-Erecting Lunar Towers

The load carrying performance for deployable lunar tower architectures utilizing thin-ply composite booms are parametrically evaluated in this paper. These self-erecting towers support a communications and sensor platform on the lunar surface. The primary design consists of the Corrugated Rollable Tubular Boom (COROTUB) acting as the primary mast of the tower which supports the payload tip mass. Guy wires connecting the boom tip to spreader bars attached at the base deployer serve to correct any lateral eccentricities and provide dimensional stability. The analysis compares the maximum tip mass that this primary design can carry relative to a reference design that does not contain any guy wires or spreader bars. The results provide insights into whether the additional system mass and complexity associated with the cable-stayed design is worth the increased payload mass that can be supported. By setting the design limiting condition to a nominal maximum allowable deflection, the analysis reveals that primary design outperforms reference design for the majority of the design space. This outcome is found to be the case when the spreader bar length relative to the tower height exceeds a critical value. The analysis is extended to lunar towers with Collapsible Tubular Mast (CTM) booms, and similar trends are observed.

Jacob G Daye↗

Dual Quaternion Visual Servo Control

This paper focuses on a dual quaternion-based estimation and control approach for position-based visual servoing (PBVS). The pose estimation of the camera is achieved using a dual quaternion-based Extended Kalman Filter (EKF), which estimates the position and orientation of the camera based on feature points acquired through a sequence of camera images. Based on the estimation, a dual quaternion control law is developed to regulate the camera to the desired pose. Leveraging the local exponential stability of the EKF and the global exponential stability of the designed controller, a nonlinear separation principle is used to prove the stability of the joint estimation and control for PBVS. The method is distinguished from other PBVS methods in the sense that a compact representation of dual quaternion is used to represent the pose, and a joint stability of estimator and controller for PBVS in dual quaternion space is presented. The proposed dual quaternion PBVS method is validated using a simulation.

Quaternions↗

Development of the X-57 Mod III/IV Piloted Simulator and Discussion of the Resultant Flighting Qualities Predictions

The all-electric X-57 Mod III aircraft was designed to demonstrate wingtip propulsion drag reduction benefits on a high-aspect-ratio wing. The X-57 Mod IV aircraft was designed to expand on the Mod III version and demonstrate the benefits of distributed electric propulsion. This paper discusses the development of theall-electricX-57 Mod III and Mod IV fixed-base pilot-in-the-loop nonlinear simulator. The paper describes the model development and the simulator cockpit construction for uses such as flight training and flying qualities analysis. Results from a stability and flying qualities analysis for both Mod III and Mod IV are presented. These results predict that both Mod III and Mod IV are stable throughout their flight envelopes. Mod III is predicted to have satisfactory flying qualities while Mod IV is predicted to have areas of adequate flying qualities.

Ryan Wallace↗

Dynamic Phasor Modeling of Three Phase Voltage Source Inverters

With the increase in the development and implementation of distributed energy resources, application of parallel connected inverters is increasing and as a result having accurate modeling and simulation tools that can help in the design, analysis and stability assessment of the grid is of great importance. Development of fundamental methods that can achieve accurate, reliable, and computationally efficient results can be very beneficial. Application of Dynamic phasor (DP) modeling method has been limited to study of limited harmonics and small combination of interconnected converters due to the complexity associated with developing models that describe larger systems. In this paper, application of DP modeling method is expanded to model any number of parallel connected three phase voltage source inverter (VSI) with inclusion of a wider harmonic content including fundamental, subharmonics, inter-harmonics, switching frequency and their sidebands. Results achieved from this modeling method is compared with conventional average model as well as detailed switching model and the effect of inclusion of wider harmonic content on accuracy of DP modeling method is demonstrated.

Xue, Yaosuo↗

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗

A new finite element formulation for computational fluid dynamics. IX - Fourier analysis of space-time Galerkin/least-squares algorithms

A Fourier stability and accuracy analysis of the space-time Galerkin/least-squares method as applied to a time-dependent advective-diffusive model problem is presented. Two time discretizations are studied: a constant-in-time approximation and a linear-in-time approximation. Corresponding space-time predictor multi-corrector algorithms are also derived and studied. The behavior of the space-time algorithms is compared to algorithms based on semidiscrete formulations.

Shakib, Farzin↗

Stability derivatives of triangular wings at supersonic speeds

The analysis of the stability derivatives of low-aspect-ratio triangular wings at subsonic and supersonic speeds, given in NACA TN no. 1423, is extended to apply to triangular wings having large vertex angles and traveling at supersonic speeds. The lift, rolling moment due to sideslip, and damping in roll and pitch for this more general case have been treated elsewhere on the basis of the theory of small disturbances. The surface potentials for angle of attack and rolling taken therefrom are used to obtain the several side-force and yawing-moment derivatives that depend on leading-edge suction, and a tentative value for the rolling moment due to yawing. The lift and moment due to downward acceleration are obtained on the basis of an unpublished unsteady-flow solution. All the known stability derivatives of the triangular wing at supersonic speeds, regardless of source, are summarized for convenience and presented with respect to both body axes and stability axes. The results are limited to Mach numbers for which the triangular wing is contained within the Mach cone for its vertex. The spanwise variation of Mach number in the case of yawing is neglected, although the effect must be of importance.

Ribner, Herbert S↗

Passive Roll Stabilization of the Near Earth Asteroid Scout Solar Sail Mission

The Near Earth Asteroid (NEA) Scout is a small satellite, solar sail mission set to launch on Artemis I. Analysis of the NEA Scout solar radiation pressure model determined that, for certain solar incidence angles, there exists at least one locally stable equilibrium point about the “roll” axis, or the axis normal to the solar sail plane. This analysis is extended to three other geometrically similar solar sail models, and a Lyapunov stability analysis is conducted demonstrating that the stable equilibria are locally stable for the majority of roll angles in a wide range of solar incidence angles. Under certain assumptions, it is possible that only two axis control is required for non-spinning solar sails. If the roll axis is controlled, an active rate damping controller can be designed to asymptotically stabilize the system. Simulation results are presented that follow the stability analysis, where in the homogenous case, the system enters a limit cycle and, in the actively damped case, the system asymptotically converges to the equilibrium point.

Ivan Rodrigues Bertaska↗

Probing Thermal Stability of Li-Ion Battery Ni-Rich Layered Oxide Cathodes by means of Operando Gas Analysis and Neutron Diffraction

Thermal stability is a crucial parameter that must be considered within the overall performance metrics of Ni-rich layered oxide cathodes. While the intrinsic structural stability of the cathode materials under thermal conditions is important, it is also critical to consider their reaction with electrolytes. In this study, operando gas analysis, ex situ neutron diffraction, and differential scanning calorimetry were combined to give a broader picture of the thermal stability of Ni-rich NMC cathodes. Li 1– x Ni 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) composite materials with four different states of charge were investigated with and without the presence of an electrolyte. It has been found that electrolyte can greatly accelerate both the structural and thermal decomposition of the cathode materials. A higher state of charge will also make cathode materials more susceptible to thermal shock. Without an electrolyte, O 2 release inducing a structural change from layered to rock salt was the major observation during thermal runaway. However, all samples retained some levels of layered structure after annealing up to 300 °C without the electrolyte. In comparison, almost all the O3-type layered phase transformed to the rock-salt phase for the cathode materials heated with the electrolyte at the same experimental conditions. Moreover, the amount of CO 2 increased by nearly two orders of magnitude when annealing cathode materials with electrolyte in operando gas analysis experiments. More delithiated NMC811 samples released larger amounts of CO 2 at earlier onset temperatures, resonating well with the differential scanning calorimetry (DSC) results showing that more delithiated samples release more heat during thermal shock.

25 ENERGY STORAGE↗