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.

SourceResearch connectionFrom yearThrough year

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

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Relative Navigation for Formation Flying of Spacecraft

This paper presents a robust and efficient approach for relative navigation and attitude estimation of spacecraft flying in formation. This approach uses measurements from a new optical sensor that provides a line of sight vector from the master spacecraft to the secondary satellite. The overall system provides a novel, reliable, and autonomous relative navigation and attitude determination system, employing relatively simple electronic circuits with modest digital signal processing requirements and is fully independent of any external systems. Experimental calibration results are presented, which are used to achieve accurate line of sight measurements. State estimation for formation flying is achieved through an optimal observer design. Also, because the rotational and translational motions are coupled through the observation vectors, three approaches are suggested to separate both signals just for stability analysis. Simulation and experimental results indicate that the combined sensor/estimator approach provides accurate relative position and attitude estimates.

Alonso, Roberto↗

NASA NTRS2001

Triana Safehold: A New Gyroless, Sun-Pointing Attitude Controller

Triana is a single-string spacecraft to be placed in a halo orbit about the sun-earth Ll Lagrangian point. The Attitude Control Subsystem (ACS) hardware includes four reaction wheels, ten thrusters, six coarse sun sensors, a star tracker, and a three-axis Inertial Measuring Unit (IMU). The ACS Safehold design features a gyroless sun-pointing control scheme using only sun sensors and wheels. With this minimum hardware approach, Safehold increases mission reliability in the event of a gyroscope anomaly. In place of the gyroscope rate measurements, Triana Safehold uses wheel tachometers to help provide a scaled estimation of the spacecraft body rate about the sun vector. Since Triana nominally performs momentum management every three months, its accumulated system momentum can reach a significant fraction of the wheel capacity. It is therefore a requirement for Safehold to maintain a sun-pointing attitude even when the spacecraft system momentum is reasonably large. The tachometer sun-line rate estimation enables the controller to bring the spacecraft close to its desired sun-pointing attitude even with reasonably high system momentum and wheel drags. This paper presents the design rationale behind this gyroless controller, stability analysis, and some time-domain simulation results showing performances with various initial conditions. Finally, suggestions for future improvements are briefly discussed.

Chen, J.↗

NASA NTRS2001

Hypersonic Boundary Layer Instability Over a Corner

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Balakumar, Ponnampalam↗

NASA NTRS2001

Noise Sources in a Low-Reynolds-Number Turbulent Jet at Mach 0.9

The mechanisms of sound generation in a Mach 0.9, Reynolds number 3600 turbulent jet are investigated by direct numerical simulation. Details of the numerical method are briefly outlined and results are validated against an experiment at the same flow conditions. Lighthill's theory is used to define a nominal acoustic source in the jet, and a numerical solution of Lighthill's equation is compared to the simulation to verify the computational procedures. The acoustic source is Fourier transformed in the axial coordinate and time and then filtered in order to identify and separate components capable of radiating to the far field. This procedure indicates that the peak radiating component of the source is coincident with neither the peak of the full unfiltered source nor that of the turbulent kinetic energy. The phase velocities of significant components range from approximately 5% to 50% of the ambient sound speed which calls into question the commonly made assumption that the noise sources convect at a single velocity. Space-time correlations demonstrate that the sources are not acoustically compact in the streamwise direction and that the portion of the source that radiates at angles greater than 45 deg. is stationary. Filtering non-radiating wavenumber components of the source at single frequencies reveals that a simple modulated wave forms for the source, as might be predicted by linear stability analysis. At small angles from the jet axis the noise from these modes is highly directional, better described by an exponential than a standard Doppler factor.

Freund, Jonathan B.↗

NASA NTRS2002

Boundary Layer Transition on Slender Cones in Conventional and Low Disturbance Mach 6 Wind Tunnels

An experimental investigation was conducted on a 5-degree half-angle cone and a 5-degree half-angle flared cone in a conventional Mach 6 wind tunnel to examine the effects of facility noise on boundary layer transition. The influence of tunnel noise was inferred by comparing transition onset locations determined from the present test to that previously obtained in a Mach 6 low disturbance quiet tunnel. Together, the two sets of experiments are believed to represent the first direct comparison of transition onset between a conventional and a low disturbance wind tunnel using a common test model and transition detection technique. In the present conventional hypersonic tunnel experiment, separate measurements of heat transfer and adiabatic wall temperatures were obtained on the conical models at small angles of attack over a range of Reynolds numbers, which resulted in laminar, transitional, and turbulent flow. Smooth model turbulent heating distributions are compared to that obtained with transition forced via discrete surface roughness. The model nosetip radius was varied to examine the effects of bluntness on transition onset. Despite wall to total temperature differences between the transient heating measurements and the adiabatic wall temperature measurement, the two methods for determining sharp cone transition onset generally yielded equivalent locations. In the 'noisy' mode of the hypersonic low disturbance tunnel, transition onset occurred earlier than that measured in the conventional hypersonic tunnel, suggesting higher levels of freestream acoustic radiation relative to the conventional tunnel. At comparable freestream conditions, the transition onset Reynolds number under low disturbance conditions was a factor of 1.3 greater than that measured on flared cone in the LaRC conventional hypersonic tunnel and a factor of 1.6 greater that the flared cone run in the low disturbance tunnel run 'noisy'. Navier-Stokes mean flow computations and linear stability analysis were conducted to assess the experimental results and have indicated N factors associated with sharp flared cone transition onset to be approximately a factor of 2 lower than that inferred from the corresponding low disturbance tunnel measurements.

Horvath, Thomas J.↗

NASA NTRS2003

Effects of Traveling Magnetic Field on Dynamics of Solidification

The Lorentz body force induced in electrically conducting fluids can be utilized for a number of materials processing technologies. An application of strong static magnetic fields can be beneficial for damping convection present during solidification. On the other hand, alternating magnetic fields can be used to reduce as well as to enhance convection. However, only special types of time dependent magnetic fields can induce a non-zero time averaged Lorentz force needed for convection control. One example is the rotating magnetic field. This field configuration induces a swirling flow in circular containers. Another example of a magnetic field configuration is the traveling magnetic field (TMF). It utilizes axisymmetric magnetostatic waves. This type of field induces an axial recirculating flow that can be advantageous for controlling axial mass transport, such as during solidification in long cylindrical tubes. Incidentally, this is the common geometry for crystal growth research. The Lorentz force induced by TMF can potentially counter-balance the buoyancy force, diminishing natural convection, or even setting up the flow in reverse direction. Crystal growth process in presence of TMF can be then significantly modified. Such properties as the growth rate, interface shape and macro segregation can be affected and optimized. Melt homogenization is the other potential application of TMF. It is a necessary step prior to solidification. TMF can be attractive for this purpose, as it induces a basic flow along the axis of the ampoule. TMF can be a practical alloy mixing method especially suited for solidification research in space. In the theoretical part of this work, calculations of the induced Lorentz force in the whole frequency range have been completed. The basic flow characteristics for the finite cylinder geometry are completed and first results on stability analysis for higher Reynolds numbers are obtained. A theoretical model for TMF mixing is also developed. In the experimental part, measurements of flow induced by TMF in a column of mercury (Hg) are presented. Also, an alloy mixing of Bi-Sn of the eutectic composition is demonstrated. A traveling magnetic field of 4mT at 3kHz applied for 120 minutes is found to be sufficient to homogenize an alloy enclosed in a 1cm diameter and 12 cm long tube.

Source record↗

NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗

NASA NTRS2004

Cross Flow Effects on Glaze Ice Roughness Formation

The present study examines the impact of large-scale cross flow on the creation of ice roughness elements on the leading edge of a swept wing under glaze icing conditions. A three-dimensional triple-deck structure is developed to describe the local interaction of a 3 D air boundary layer with ice sheets and liquid films. A linear stability analysis is presented here. It is found that, as the sweep angle increases, the local icing instabilities enhance and the most linearly unstable modes are strictly three dimensional.

Tsao, Jen-Ching↗

NASA NTRS2003

A Comparison between Growth Morphology of "Eutectic" Cells/Dendrites and Single-Phase Cells/Dendrites

Directionally solidified (DS) intermetallic and ceramic-based eutectic alloys with an in-situ composite microstructure containing finely distributed, long aspect ratio, fiber, or plate reinforcements are being seriously examined for several advanced aero-propulsion applications. In designing these alloys, additional solutes need to be added to the base eutectic composition in order to improve heir high-temperature strength, and provide for adequate toughness and resistance to environmental degradation. Solute addition, however, promotes instability at the planar liquid-solid interface resulting in the formation of two-phase eutectic "colonies." Because morphology of eutectic colonies is very similar to the single-phase cells and dendrites, the stability analysis of Mullins and Sekerka has been extended to describe their formation. Onset of their formation shows a good agreement with this approach; however, unlike the single-phase cells and dendrites, there is limited examination of their growth speed dependence of spacing, morphology, and spatial distribution. The purpose of this study is to compare the growth speed dependence of the morphology, spacing, and spatial distribution of eutectic cells and dendrites with that for the single-phase cells and dendrites.

Tewari, S. N.↗

NASA NTRS2004

Effects of Gravity on Bubble Formation in an Annular Jet

The effects of gravity on the bubble formation in an annular jet were studied. The experiments were conducted in the 2.2-second drop tower at the NASA Glenn Research Center. Terrestrial gravity experiments were conducted at the Fluid Dynamics Research Laboratory at the University of Oklahoma. Stainless steel tubing with inner diameters of 1/8" (gas inner annulus) and 5/16" (liquid outer annulus) served as the injector. A rectangular test section, 6" x 6" x 14" tall, made out of half-inch thick Lexan was used. Images of the annular jet were acquired using a high-speed camera. The effects of gravity and varying liquid and gas flow rates on bubble size, wavelength, and breakup length were documented. In general, the bubble diameter was found to be larger in terrestrial gravity than in microgravity for varying Weber numbers (0.05 - 0.16 and 5 - 11) and liquid flow rates (1.5 ft/s - 3.0 ft/s). The wavelength was found to be larger in terrestrial gravity than in microgravity, but remained constant for varying Weber numbers. For low Weber numbers (0.05 - 0.16), the breakup length in microgravity was significantly higher than in terrestrial gravity. Comparison with linear stability analysis showed estimated bubble sizes within 9% of experimental bubble sizes. Bubble size compared to other terrestrial gravity experiments with same flow conditions showed distinct differences in bubble size, which displayed the importance of injector geometry on bubble formation.

Koepp, R. A.↗

NASA NTRS2004

Effects of Mean Flow Profiles on the Instability of a Low-Density Gas Jet Injected into a High-Density Gas

The effects of the mean velocity profiles on the instability characteristics in the near-injector region of axisymmetric low density gas jets injected vertically upwards into a high-density gas medium were investigated using linear inviscid stability analysis. The flow was assumed to be isothermal and locally parallel. Three velocity profiles, signifying different changes in the mean velocity in the shear layer, were used in the analysis. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the instability for each set of mean profiles were delineated. At a large Froude number (negligible gravity), a critical density ratio was found for the three profiles at which the jet became absolutely unstable. The critical density ratio for each velocity profile was increased as the Froude number was reduced. A critical Froude number was found for the three sets of profiles, below which the jet was absolutely unstable for all the density ratios less than unity, which demarcated the jet flow into the momentum-driven regime and the buoyancy-driven regime.

Vedantam, NandaKishore↗

NASA NTRS2005

The Use of a Gyroless Wheel-Tach Controller in SDO Safehold Mode

This paper describes the progression of the Safehold mode design on the Solar Dynamics Observatory satellite. Safehold uses coarse Sun sensors and reaction wheel tachometers to keep the spacecraft in a thermally safe and power-positive attitude. The control algorithm is described, and simulation results shown. Specific control issues arose when the spacecraft entered eclipse, and a description of the trade study which added gyroscopes to the mode is included. The paper concludes with the results from the linear and nonlinear stability analysis.

Bourkland, Kristin L.↗

NASA NTRS2006

A Measure-Theoretic Proof of the Markov Property for Hybrid Systems with Markovian Inputs

The behavior of a general hybrid system in discrete time can be represented by a non-linear difference equation x(k+1) = Fk(x(k), theta(k)), where theta(k) is assumed to be a finite state Markov chain. An important step in the stability analysis of these systems is to establish the Markov property of (x(k), theta(k)). There are, however, no complete proofs of this property which are simple to understand. This paper aims to correct this problem by presenting a complete and explicit proof, which uses only basic measure-theoretical concepts.

Tejada, Arturo↗

NASA NTRS2006

On the Markov Property for Nonlinear Discrete-Time Systems with Markovian Inputs

The behavior of a general hybrid system in discrete-time can be represented by a non-linear difference equation x(k + 1) = Fk(x(k), theta(k)), where theta(k) is assumed to be a finite-state Markov chain. An important step in the stability analysis of these systems is to establish the Markov property of (x(k), theta(k)). There are, however, no complete proofs of this property which are simple to understand. This paper aims to correct this problem by presenting a complete and explicit proof, which uses only fundamental measure-theoretical concepts.

Tejada, Arturo↗

NASA NTRS2008

Ares-I Bending Filter Design using a Constrained Optimization Approach

The Ares-I launch vehicle represents a challenging flex-body structural environment for control system design. Software filtering of the inertial sensor output is required to ensure adequate stable response to guidance commands while minimizing trajectory deviations. This paper presents a design methodology employing numerical optimization to develop the Ares-I bending filters. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares-I time-varying dynamics and control system can be frozen over a short period of time, the bending filters are designed to stabilize all the selected frozen-time launch control systems in the presence of parameter uncertainty. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constrains minimizes performance degradation caused by the addition of the bending filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The bending filter designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC MAVERIC 6DOF nonlinear time domain simulation.

Hall, Charles↗

NASA NTRS2009

Re-evaluation of an Optimized Second Order Backward Difference (BDF2OPT) Scheme for Unsteady Flow Applications

Recent experience in the application of an optimized, second-order, backward-difference (BDF2OPT) temporal scheme is reported. The primary focus of the work is on obtaining accurate solutions of the unsteady Reynolds-averaged Navier-Stokes equations over long periods of time for aerodynamic problems of interest. The baseline flow solver under consideration uses a particular BDF2OPT temporal scheme with a dual-time-stepping algorithm for advancing the flow solutions in time. Numerical difficulties are encountered with this scheme when the flow code is run for a large number of time steps, a behavior not seen with the standard second-order, backward-difference, temporal scheme. Based on a stability analysis, slight modifications to the BDF2OPT scheme are suggested. The performance and accuracy of this modified scheme is assessed by comparing the computational results with other numerical schemes and experimental data.

Vatsa, Veer N.↗

NASA NTRS2010

Ares I Flight Control System Design

The Ares I launch vehicle represents a challenging flex-body structural environment for flight control system design. This paper presents a design methodology for employing numerical optimization to develop the Ares I flight control system. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares I time-varying dynamics and control system can be frozen over a short period of time, the flight controllers are designed to stabilize all selected frozen-time launch control systems in the presence of parametric uncertainty. Flex filters in the flight control system are designed to minimize the flex components in the error signals before they are sent to the attitude controller. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constraints minimizes performance degradation caused by the addition of the flex filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The flight control system designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC 6DOF nonlinear time domain simulation.

Jang, Jiann-Woei↗

NASA NTRS2011

On Using Exponential Parameter Estimators with an Adaptive Controller

Typical adaptive controllers are restricted to using a specific update law to generate parameter estimates. This paper investigates the possibility of using any exponential parameter estimator with an adaptive controller such that the system tracks a desired trajectory. The goal is to provide flexibility in choosing any update law suitable for a given application. The development relies on a previously developed concept of controller/update law modularity in the adaptive control literature, and the use of a converse Lyapunov-like theorem. Stability analysis is presented to derive gain conditions under which this is possible, and inferences are made about the tracking error performance. The development is based on a class of Euler-Lagrange systems that are used to model various engineering systems including space robots and manipulators.

Patre, Parag↗