Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear 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 163 records · Page 9

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.↗

Separated flows near the nose of a body of revolution

The solution of the Navier-Stokes equations for the problem of cross-flow separataion about a deforming cylinder was achieved by iteration. It was shown that the separation starts at the rear stagnation point and the point of primary separation moves upstram along the cylinder surface. A general method of linear stability analysis for nonparallel external flows was constructed, which consists of representing the eigenfunctions with complete orthogonal sets and forms characteristic equations with the Galerkin method. The method was applied to the Kovasznay flow which is an exact solution of the Navier-Stokes equation. The results show that when the critical parameter is exceeded, there are only a few isolated unstable eigen-frequencies. Another exact solution is shown to be absolutely and monotonically stable with respect to infinitesimal disturbances of all frequencies. The flow is also globally, asymptotically, and monotonically stable in the mean with respect o three-dimensional disturbances. This result forms the sound foundation of rigorous stability analysis for nonparallel flows, and provides an invaluable test ground for future studies of nonparallel flows in which the basic states do not posses exact solutions. The application of this method to the study of the formation of spiral vorticies near the nose of a rotating body of revolution is underway. The same method will be applied to the stability analysis of reversed flow over a plate with suction.

Lin, S. P.↗

Mechanisms by which heat release affects the flow field in a chemically reacting, turbulent mixing layer

The mechanisms by which heat release affects the fluid dynamics in a turbulent reacting mixing layer are studied by direct numerical simulation. In agreement with previous laboratory experiments, the heat release is observed to lower the rate at which the mixing layer grows and to reduce the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer are shown to produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of fluid elements. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers are shown to result from vorticity generation by baroclinic torques. The density reductions also lower the generation rates of turbulent kinetic energy and the turbulent shear stresses, resulting in less turbulent mixing of fluid elements. Calculations of the energy in the various wave numbers show that the heat release has a stabilizing effect on the growth rates of individual modes. A linear stability analysis of a simlified model problem confirms this, showing that low density fluid in the mixing region will result in a shift of the frequency of the unstable modes to lower wave numbers (longer wavelengths). The growth rates of the unstable modes decrease, contributing to the slower growth of the mixing layer.

Riley, J. J.↗

Small-scale motions observed by aircraft in the tropical lower stratosphere - Evidence for mixing and its relationship to large-scale flows

Measurements of temperature and ozone from instrumented aircraft in the tropical lower stratosphere show the presence of small-scale disturbances generated by underlying cumulus convection and Kelvin-Helmholtz instability. The disturbances associated with underlying convection have peak-to-peak vertical parcel excursions of about 300 m. Flying conditions were smooth, suggesting an ensemble of gravity waves and little or no turbulent mixing. It is nevertheless possible that these waves break at other altitudes, leading to turbulent mixing and net fluxes of vertically stratified tracers. Disturbances attributed to KH instability implied vertical parcel excursions of 300-400 m. The disturbances coincided with rough flying conditions, suggesting turbulent mixing. A linear stability analysis of the atmospheric basic state defined by high-resolution radiosondes shows fastest growing waves with horizontal wavelengths of 1.4-1.8 km, consistent with the aircraft observations. The strong shears responsible for the KH instability are due to large-scale waves propagating into a region of small intrinsic frequency. Radiosonde observations show that the zonal length scale of these waves is about 1000 km.

Pfister, Leonhard↗

Effect of a forced Couette flow on coupled convective and morphological instabilities during unidirectional solidification

The effect of a forced Couette flow, parallel to a horizontal crystal-melt interface during directional solidification of an alloy of lead containing tin, on the onset of convective and morphological instabilities, is calculated numerically via a linear stability analysis. Such a flow does not affect perturbations with wave vectors perpendicular to the flow. For perturbations with wave vectors parallel to the flow, the onset of morphological instability is somewhat suppressed and thermosolutal convection is greatly suppressed. When instabilities occur, they are oscillatory and correspond to travelling waves. For values of the crystal growth velocity for which mixed morphological and convective modes occur, the presence of a forced flow produces sufficient decoupling to allow otherwise degenerate branches to be identified.

Coriell, S. R.↗

Boundary-layer stability and airfoil design

Several different natural laminar flow (NLF) airfoils have been analyzed for stability of the laminar boundary layer using linear stability codes. The NLF airfoils analyzed come from three different design conditions: incompressible; compressible with no sweep; and compressible with sweep. Some of the design problems are discussed, concentrating on those problems associated with keeping the boundary layer laminar. Also, there is a discussion on how a linear stability analysis was effectively used to improve the design for some of the airfoils.

Viken, Jeffrey K.↗

The effect of an electric field on the morphological stability of the crystal-melt interface of a binary alloy

A fully time-dependent linear stability analysis of the morphological stability of a planar interface during directional solidification of a binary alloy at constant velocity in the presence of an electric field, is performed. The electromigration of solute and the differing electrical conductivities of solid and liquid for a model in which the temperature gradient is constant are taken into account. The present results are compared with the constitutional supercooling criterion, and it is shown there may be substantial differences. A modified constitutional supercooling criterion which is valid over a large range of conditions is derived. It is also found under certain conditions that the onset of instability may be time dependent.

Wheeler, A. A.↗

The possible role of meridional flows in suppressing magnetic buoyancy

The equation of motion for a toroidal flux ring in a stellar convective envelope is derived, and the equilibrium of such a ring is considered. Necessary conditions for the stability of toroidal flux rings are derived, and results of stability calculations for a particular model of the meridional flow are presented. The motions of the flux rings when the rings are far from their equilibrium position or when equilibrium does not exist are considered. The results confirm the linear stability analysis, and show that in the absence of stable equilibrium, the rings move toward the solar surface along a trajectory which is parallel to the rotation axis. It is expected that viscosity will tend to reduce the rotational velocity difference between the flux ring and its surroundings, thus reducing the Coriolis force and altering the equilibrium. The storage time of toroidal flux rings is estimated, and some implications for the sun are discussed.

Van Ballegooijen, A. A.↗

A delayed action oscillator for ENSO

A simple nonlinear model is proposed for the El Nino/Southern Oscillation (ENSO) phenomenon. Its key feature is the inclusion of oceanic wave transit effects through a negative, delayed feedback. A linear stability analysis and numerical results are presented, showing that the period of the oscillation is typically several times the delay. It is suggested that such an effect can account for the long time scale of ENSO.

Suarez, Max J.↗

Effects of heat release on the structure and stability of a coflowing, chemically reacting jet

It is shown how linear stability analysis and full simulations can be used to study the stability and structure of reacting flows. The small effect of gravity on steady laminar flame heights can be modeled if allowance is made for the density variations due to heat release. It is found that cold transitional jet simulations yield perturbation growth rates which are consistent with the predictions of the inviscid stability theory.

Mahalingam, S.↗

Buoyancy effects on morphological instability during directional solidification

The onset of morphological instability during the directional solidification of a single-phase binary alloy at constant velocity vertically upwards is treated by a linear stability analysis. The case in which a heavier solute is rejected at the solidifying interface is considered, and the effect of natural convection on the critical concentration for the onset of instability is studied. For tin containing lead, a small destabilization of the system at low growth velocities, and a large increase in the wavelength of the instability at the onset are found. Calculations show that the destabilization is enhanced as the variation of density with solute concentration is reduced, and in the limit of neutrally-dense solute, there is a long wavelength instability for which the critical solute concentration is several orders of magnitude lower than that predicted by the Mullins and Sekerka (1964) analysis in the absence of convection. For the neutrally-dense solute, a simplified analysis indicates the roles played by the interface deformation and thermal convection in promoting the instability. In particular, the destabilization is very sensitive to the ratio of crystal and melt thermal conductivities.

Coriell, S. R.↗

Absolute and convective instability of a cylinder wake

The stability behavior of a circular cross-section cylinder's wake at Reynolds number values of up to 45 is presently investigated by means of local linear stability theory. The steady-wake profiles computed are Navier-Stokes solutions of a uniform, incompressible viscous flow around a cylinder obtained by a spectral method. An absolutely unstable region is found to begin to form at a Reynolds number of about 20, and grows with incresing Reynolds number. The onset of global instability response must be characterized by a critical length of an absolutely unstable region; a critical Reynolds number criterion and preferred frequency are proposed based on linear stability analysis.

Yang, X.↗

The evolution of cooling flows - Self-similar cooling waves

Similarity solutions are presented for the time-dependent evolution of cooling flows. A cooling flow expands when the cooling time computed from the initial gas distribution increases with radius. If certain simplifying conditions are met, the evolution of the resulting cooling wave becomes self-similar. The similarity solutions obtained here assumed subsonic flow, which is valid in the outer parts of observed cooling flows; they are matched to transonic accretion solutions valid in the central parts. The models are applied to the cooling flow onto M87, with unsatisfactory agreement at small radius. A local linear stability analysis shows that the similarity solution suffers isobaric thermal instability, suggesting that the neglect of star formation may be responsible for the poor agreement. The similarity solution suggests that the cooling flow around M87 was more vigorous in the past.

Bertschinger, Edmund↗

Curvature effects on the stability of three-dimensional laminar boundary layers

The linear stability equations for compressible, three-dimensional laminar boundary layer flow are derived in an orthogonal curvilinear coordinate system. The system of equations is solved using a finite difference scheme in order to study the effects of streamline and surface curvature and compressibility on the stability of the flow past a swept wing. It is known that convex surface curvature can have a stabilizing effect on the laminar boundary layer. Conversely, concave surface curvature can be destabilizing. The magnitude of these effects for swept wings is determined. Results indicate that amplification rates and hence, N-factors, for the flow over the convex upper surface of a swept wing can be reduced by about 15 to 45 percent when curvature effects are included in the linear stability analysis. The results of the calculations show that concave curvature destabilizes crossflow type disturbances with a significant increase in amplification rate. In addition, comparisons are made with some experimental results on a swept concave-convex surface. Calculated velocity vector plots show good agreement with observed disturbances in the laminar boundary layer over the concave surface.

Collier, F. S., Jr.↗

The role of gravity on macrosegregation in alloys

During dendritic solidification liquid flow is induced both by buoyancy forces and solidification shrinkage. There is strong evidence that the major reason for the liquid flow is the former, i.e., thermosolutal convection. In the microgravity environment, it is thought that the thermosolutal convection will be greatly diminished so that convection will be confined mainly to the flow of interdendritic liquid required to satisfy the solidification shrinkage. An attempt is made to provide improved models of dendritic solidification with emphasis on convection and macrosegregation. Macrosegregation is an extremely important subject to the commercial casting community. The simulation of thermosolutal convection in directionally solidified (DS) alloys is described. A linear stability analysis was used to predict marginal stability curves for a system that comprises a mushy zone underlying an all-liquid zone. The supercritical thermosolutal convection in directionally solidified dendritic alloys was also modeled. The model assumes a nonconvective initial state with planar and horizontal isotherms and isoconcentration that move upward at a constant solidification velocity. Results are presented for systems involving lead-tin alloys and show significant differences with results of plane-front solidification.

Poirier, D. R.↗