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 973 records · Page 54

A split delta-V technique for drift control of geosynchronous spacecraft

A split-delta-V technique is presented for the east-west stationkeeping of geosynchronous spacecraft using the linearized optimal transfer strategies for near-circular coplanar orbits. The maneuvers are carried out by taking advantage of the dynamics of the spacecraft in order to cancel out the effect of the drift accumulated between successive impulsive delta-V's. A stability analysis is then carried out in order to determine the effect of implementing a series of precanned maneuvers at the end of each drift cycle in an autonomous mode by taking into account initial orbit determination uncertainties and execution errors emanating from each delta-V. It is shown that small tolerance deadbands are violated after only a few cycles and the system is unstable.

Kechichian, F. A.↗

On the OLP prediction of the unstable modes of the flat plate turbulent boundary layer

The Orr et al. (OLP) stability analysis is reformulated to yield profiles from hot-wire velocity covariance data in the boundary layer of a flow over a flat plate. The proper orthogonal decomposition theorem (PODT) is used to maximize the mean square of the inner product of the velocity vector with a deterministic candidate vector function in Hilbert space. The OLP formulation couples the PODT for the eddy viscosity with a global extremum principle for the disturbance kinetic energy. Derivation of the governing equations is described analytically, as are the extraction of the mean velocity profile and the eddy viscosity. Finally, techniques for characterizing the dominant mode of the velocity field are introduced and applied to three-dimensional flow over a flat plate. Good agreement is found between the predictions, using laboratory data, and the measured energy and three-dimensional structures of the dominant eddies.

Wadia, A. R.↗

Flare loop radiative hydrodynamics. IV - Dynamic evolution of unstable semiempirical loop models

The evolution of the unstable solar atmosphere into the nonlinear phase, in response to various perturbations, is followed. The initial dynamic evolution of the atmosphere follows the predictions of linear stability analysis. In the nonlinear phase, rapid changes are confined to the transition region; these changes are manifested as a propagation of the transition region through the plasma, i.e., chromospheric evaporation or condensation. Global evolution therefore proceeds on the coronal conductive time scale. The rate of propagation of the transition region is determined by the imbalance between the energy supplied by thermal conduction from the corona and radiative cooling within the transition region itself. Flow velocities in the lower corona during evaporation or condensation are, in the cases studied, of order 3 km/s. The observed dynamic evolution is consistent with the existence of relatively long-lived coronal loops whose brightnesses vary on the evaporative time scale.

An, C.-H.↗

Unsymmetrically laminated composites

Classical lamination theory predicts that the room-temperature shapes of all elevated-temperature cure, unsymmetrically laminated composites are saddles. However, experimental observation indicates that the shapes are often cylindrical. In addition, a second cylindrical can sometimes be obtained from the first by a simple snap-through action. A geometrically nonlinear extension to classical lamination theory is used to explain this behavior. Approximate solutions to the nonlinear extension are obtained by using a Rayleigh-Ritz minimization of the laminate's total potential energy. A stability analysis explains the dual cylindrical shapes.

Hyer, M. W.↗

The Riemann disks. II - Stability

The Riemann disks are the two-dimensional analogs of the Riemann S-type ellipsoids. They form a simple yet rich family of galactic bar models. A method is developed to find their normal modes to arbitrary order, and the stability analysis is carried out for all modes with order not greater than 5.

Weinberg, M. D.↗

Studies of convection in a solidifying binary mixture at reduced gravity

A great deal of interest was generated recently in the possibility of producing new materials in the reduced gravity environment provided during the forthcoming missions of Spacelab. The range of possibilities extend from producing large crystals of uniform properties to manufacturing materials with unique properties. Most of these processes involve the solidification of materials from the liquid state. Convective motions within the liquid during solidification can influence the local material composite and the shape of the solid-liquid interface which may result in solids with non-uniform properties and crystal defects. The microgravity environment of Spacelab is being viewed as one in which the buoyancy forces are eliminated so that convection driven by thermal gradients does occur, resulting in an improved solidification process. However, convection may occur for other reasons and whether convection is negligible or not during solidification constitutes processing in low-gravity environment. Little information exists presently on convection during solidification under such circumstances. A continuation of an analytical investigation into the nature of convective motion in a binary liquid layer due to surface tension forces during its solidification is reported. The onset of convection will be determined through a stability analysis which is described.

Antar, B. N.↗

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗

A selection principle in Benard-type convection

In a Benard-type convection problem, the stationary flows of an infinite layer of fluid lying between two rigid horizontal walls and heated uniformly from below are determined. As the temperature difference across the layer increases beyond a certain value, other convective motions appear. These motions areoften cellular in character in that their streamlines are confined to certain well-defined cells having, for example, the shape of rolls or hexagons. A selection principle that explains why hexagonal cells seem to be preferred for certain ranges of the parameters is formulated. An operator-theoretical formulation of one generalized Bernard problem is given. The infinite dimensional problem is reduced to one of solving a finite dimensional system of equations, namely, the selection equations. These equations are solved and a linearized stability analysis of the resultant stationary flows is presented.

Knightly, G. H.↗

Electrostatic Kelvin-Helmholtz instability in a radially injected plasma cloud

An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.

Sydora, R. D.↗

Natural laminar flow flight experiments on a swept-wing business jet

Flight test experiments have been conducted to measure the extent and nature of laminar flow on a smoothed test region of a swept-wing business jet. Surface hot film anemometers and sublimating chemicals were used for transition detection. Surface pressure distributions were measured using pressure belts. Instrumentation techniques and problems are presented. Correlation was obtained between the hot film and sublimating chemicals for transition detection. Taylor-Goertler vortices were observed for some flight conditions. Boundary layer stability analysis is being conducted using measured pressure distributions to determine the correlation possible for transition prediction on swept wings.

Wentz, W. H., Jr.↗

Effect of Scattering on Instabilities in Line-driven Stellar Winds

Line driven O-B stellar winds are unstable to perturbations of short spatial wavelength and the growth rates for such instabilities are very rapid. The nature of the nonlinear development of this instability is unknown, but might possibly be one of blobs of gas driven through ambient gas or a quasiregular train of outward moving shocks. In either case the resulting dissipation of mechanical energy might explain the observed anomalous heating n O-B stars as evidenced by their X-ray emission and high ionization state. It also might explain the observed fine structure of the absorption lines and their time variability. The driving due to the absorption of the stellar continuum flux was considered and the effects of the diffuse, scattered radiation field were neglected. It was shown that under certain special conditions the effect of scattering could reduce the instability growth rate to zero. A stability analysis that includes scattering, but that uses the more physically realistic assumption of complete redistribution instead of coherent scattering, and that includes the effects of transverse velocity gradients, which become important as the flow moves away from the stellar surface is presented. It is found that the instability is eliminated right at the base of the wind, but that as the flow moves outward the instability rate rapidly becomes equal to a substantial fraction of the calculated value, the fraction asymptotically reaching 80% at large radii. Since this still implies many e folds in a characteristic outflow time, the primary conclusion that these winds are highly unstable is unchanged.

Owocki, S. P.↗

Thermosolutal convection during directional solidification

During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly - a phenomenon also indicated by linear stability analysis. The lateral inhomogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time.

Mcfadden, G. B.↗

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗

Convective influence on the stability of a cylindrical solid-liquid interface

Experiments in which a long vertical, heated wire is surrounded by concentric annuli of a melt and its crystalline solid show that the convection state changes from a stable unicell surrounded by a stationary cylindrical solid-liquid interface, to a complex time-dependent flow surrounded by a rotating, helical solid-liquid interface. This transition occurs at a Grashof number of approximately 150, which is an order of magnitude less than the critical Grashof number calculated for a liquid annulus surrounded by rigid walls. A linear stability analysis has been carried out for an infinitely tall vertical annulus. When the deformable nature of the crystal-melt interface is taken into account in the boundary conditions, two new modes of instability arise. The most dangerous mode is asymmetrical and corresponds to helical waves travelling vertically upwards. The critical Grashof number and the scaling properties of the eigenstate agree with experiments. The results clearly demonstrate the coupling of convection with crystal-melt interfacial instabilities.

Fang, Q. T.↗

Magnetohydrodynamic thermal instabilities in cool inhomogeneous atmospheres

The stability of magnetic loops to current-driven filamentation instabilities is investigated. The unperturbed atmosphere is assumed to be composed of an (upper) isothermal optically thin low-density portion and a (lower) higher-density portion which is in radiative equilibrium; in both cases, the atmosphere is in hydrostatic equilibrium, so that gravitational stratification is taken into account. In order to provide specific equilibrium conditions for evaluation of the dispersion relation, conditions appropriate for the surface of a solar-type star are adopted; i.e., a fairly low temperature (T = 5000 K) appropriate for a 'precoronal' state associated, for example, with magnetic flux emerging from photospheric levels under the action of magnetic buoyancy. A linear stability analysis is performed, and numerical results show that physically plausible current densities, which would be generated by typical loop-footpoint motions, are effective in driving MHD instabilities in such a plasma. The instability growth rates are strongly dependent on the assumed current density distribution and on the density scale height.

Bodo, G.↗

A selection principle for Benard-type convection

In a Benard-type convection problem, the stationary flows of an infinite layer of fluid lying between two rigid horizontal walls and heated uniformly from below are determined. As the temperature difference across the layer increases beyond a certain value, other convective motions appear. These motions are often cellular in character in that their streamlines are confined to certain well-defined cells having, for example, the shape of rolls or hexagons. A selection principle that explains why hexagonal cells seem to be preferred for certain ranges of the parameters is formulated. An operator-theoretical formulation of one generalized Bernard problem is given. The infinite dimensional problem is reduced to one of solving a finite dimensional system of equations, namely, the selection equations. These equations are solved and a linearized stability analysis of the resultant stationary flows is presented.

Knightly, G. H.↗

Instabilities and transition in the wall boundary layers of low-disturbance supersonic nozzles

Linear stability analysis and experimental nozzles were employed in studying means to control Tollmein-Schlichting (TS) waves and Taylor-Goertler (TG) vortices at transition points along walls in a supersonic wind tunnel. It was hoped that control of the TS waves and TG vortices would cause the wall flows to remain laminar, thereby reducing the turbulence noise in the tunnel. The test nozzles injecting the waves and vertices had Mach numbers from 3 to 5. TG vortices triggered transition, while axisymmetric nozzle flows suppressed it. TS waves were controllable with favorable pressure graients at Mach 3.5. Finally, advances were made in the capability of predicting the locations of the transition points.

Chen, F.-J.↗

The effect of topography on the evolution of unstable disturbances in a baroclinic atmosphere

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, their stability, and the long term evolution of incipient unstable waves. The flow is forced by latitudinally dependent radiative heating. Dissipation is in the form of Rayleigh friction. An analytical solution is found for the propagating finite amplitude waves which result from baroclinic instability of the zonal winds when topography is absent. The appearance of this solution for wavelengths just longer than the Rossby radius of deformation and disappearance of ultra-long wavelengths is interpreted in terms of the Hopf bifurcation theory. Simple dynamic and thermodynamic criteria for the existence of periodic Rossby solutions are presented. A Floquet stability analysis shows that the waves are neutral. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance.

Clark, J. H. E.↗