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 199 records · Page 11

Features of wavy vortices in a curved channel from experimental and numerical studies

Results are reported from an experimental study obtaining evidence of time-dependent, wavy vortex motions associated with undulating and twisting Dean vortices in a curved channel with 40-to-1 aspect ratio, and mild curvature (radius ratio = 0.979). The results are compared with direct numerical simulations of time-dependent 3D Navier-Stokes equations using boundary conditions in the spanwise and streamwise directions. When viewed in cross section, experimental visualizations of undulating and twisting vortex flows show rocking motion and changes in the direction of the flow between vortices that are like those observed in the simulations. Experimental spectra show that undulating vortices are replaced by the higher-frequency, shorter streamwise wavelength twisting vortices at higher Reynolds numbers. When undulating vortices are present, experimental power spectra and visualizations give frequencies that are somewhat lower than the most unstable frequencies predicted by linear stability analysis.

Ligrani, P. M.↗

The effect of gravity modulation on thermosolutal convection in an infinite layer of fluid

The effect of time-periodic vertical gravity modulation on the onset of thermosolutal convection in an infinite horizontal layer with stress-free boundaries is investigated using Floquet theory for the linear stability analysis. Situations for which the fluid layer is stably stratified in either the fingering or diffusive regimes of double-diffusive convection are considered. Results are presented both with and without steady background acceleration. Modulation may stabilize an unstable base solution or destabilize a stable base solution. In addition to synchronous and subharmonic response to the modulation frequency, instability in the double diffusive system can occur via a complex conjugate mode. In the diffusive regime, where oscillatory onset occurs in the unmodulated system, regions of resonant instability occur and exhibit strong coupling with the unmodulated oscillatory frequency. The response to modulation of the fundamental instability of the unmodulated system is described both analytically and numerically; in the double-diffusive system this mode persists under subcritical conditions as a high-frequency lobe.

Saunders, B. V.↗

Cellular instability in rapid directional solidification - Bifurcation theory

Merchant and Davis performed a linear stability analysis on a model for the directional solidification of a dilute binary alloy valid for all speeds. The analysis revealed that nonequilibrium segregation effects modify the Mullins and Sekerka cellular mode, whereas attachment kinetics has no effect on these cells. In this paper, the nonlinear stability of the steady cellular mode is analyzed. A Landau equation is obtained that determines the amplitude of the cells. The Landau coefficient here depends on both nonequilibrium segregation effects and attachment kinetics. This equation gives the ranges of parameters for subcritical bifurcation (jump transition) or supercritical bifurcation (smooth transition) to cells.

Braun, R. J.↗

Thermal effects in rapid directional solidification - Linear theory

We study the morphological instability of the planar solid/liquid interface for a unidirectionally-solidified dilute binary mixture. We use a model developed by Boettinger et al. (1985, 1986), Aziz (1982), and Jackson et al. (1980), which allows for nonequilibrium effects on the interface through velocity-dependent segregation and attachment kinetics. Two types of instabilities are found in the linear stability analysis: (1) a cellular instability, and (2) an oscillatory instability driven by disequilibrium effects. Merchant and Davis (1990) characterized these instabilities subject to the frozen-temperature approximation (FTA). The present work relaxes the FTA by including the effects of latent heat and the full temperature distribution. Thermal effects slightly postpone the onset of the cellular instability but dramatically postpone the onset of the oscillatory instability; however, the absolute-stability conditions, at which at high speed the cellular and oscillatory instabilities are suppressed, remain unchanged from the FTA.

Huntley, D. A.↗

Marangoni instability in a liquid layer with two free surfaces

The onset of the Marangoni instability in a liquid layer with two free nearly insulating surfaces heated from below is studied. Linear stability analysis yields a condition for the emergence of a longwave or a finite wavelength instability from the quiescent equilibrium state. Using the method of asymptotic expansions, a weakly nonlinear evolution equation describing the spatiotemporal behavior of the velocity and temperature fields at the onset of the longwave instability is derived. The latter is given by delta(M) = 24, delta(M) being the difference between the upper and the lower Marangoni numbers. It is shown that in some parametric range one convective cell forms across the layer, while in other parametric domains two convective cells emerge between the two free surfaces.

Deissler, Robert J.↗

The effect of crossflow on Taylor vortices: A model problem

A number of practically relevant problems involving the impulsive motion or the rapid rotation of bodies immersed in fluid are susceptible to vortex-like instability modes. Depending upon the configuration of any particular problem the stability properties of any high-wavenumber vortices can take on one of two distinct forms. One of these is akin to the structure of Gortler vortices in boundary layer flows while the other is similar to the situation for classical Taylor vortices. Both the Gortler and Taylor problems have been extensively studied when crossflow effects are excluded from the underlying base flows. Recently, studies were made concerning the influence of crossflow on Gortler modes and a linearized stability analysis is used to examine crossflow properties for the Taylor mode. This work allows us to identify the most unstable vortex as the crossflow component increases and it is shown how, like the Gortler case, only a very small crossflow component is required in order to completely stabilize the flow. Our investigation forms the basis for an extension to the nonlinear problem and is of potential applicability to a range of pertinent flows.

Otto, S. R.↗

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗

Overview of the Miniature Sensor Technology Integration (MSTI) spacecraft attitude control system

Msti2 is a small, 164 kg (362 lb), 3-axis stabilized, low-Earth-orbiting satellite whose mission is missile booster tracking. The spacecraft is actuated by 3 reaction wheels and 12 hot gas thrusters. It carries enough fuel for a projected life of 6 months. The sensor complement consists of a Horizon Sensor, a Sun Sensor, low-rate gyros, and a high rate gyro for despin. The total pointing control error allocation is 6 mRad (.34 Deg), and this is while tracking a target on the Earth's surface. This paper describes the Attitude Control System (ACS) algorithms which include the following: attitude acquisition (despin, Sun and Earth acquisition), attitude determination, attitude control, and linear stability analysis.

Mcewen, Rob↗

Use of rotation to suppress thermosolutal convection in directionally solidified binary alloys

Effects of rotation on onset of convection during plane-front directional solidification of Pb-Sn and the pseudobinary system mercury cadmium telluride (Hg(1-x)Cd(x)Te), and on dendritic solidification of Pb-Sn have been studied by means of linear stability analysis. Incorporating Coriolis and centrifugal accelerations into the momentum equation of Coriell et al., we find that under realistic processing conditions, a large degree of stabilization can be achieved using modest rotation rates for both Pb-Sn and mercury cadmium telluride. At a growth velocity of 5 micron/sec and nominal liquid-side temperature gradient of 200 K/cm in Pb-Sn, rotation at 500 rpm results in a hundredfold increase in the critical Sn concentration. Large increases in the maximum allowable growth velocity at fixed melt composition are also attainable with modest rotation rates. The effect is amplified under conditions of reduced gravitational acceleration. For Hg(1-x)Cd(x)Te, we have also studied the nonrotating case. The key differences are due to the existence of a composition range for Hg(1-x)Cd(x)Te in which the melt density has a local maximum as a function of temperature. When the melt solidifies by cooling from below, the liquid density may initially increase with distance above the interface, before ultimately decreasing as the melt temperature increases above the value at which the local density maximum occurs. In contrast to the Pb-Sn case where density depends monotonically on temperature and composition, for Hg(1-x)Cd(x)Te there exists a critical value of the growth velocity above which plane-front solidification is unstable for all bulk CdTe mole fractions. Again, rotation leads to significant inhibition of onset. We identify the predicted stabilization with the Taylor-Proudman mechanism by which rotation inhibits thermal convection in a single-component fluid heated from below. In a binary liquid undergoing solidification, rotation inhibits the onset of buoyancy-driven convection, and has no effect on the short-wavelength morphological instability. At large growth velocities, the plane-front interface between liquid and solid becomes unstable with respect to a morphological instability and solidification occurs dendritically, with a mushy zone of dendrites and interdendritic fluid separating the solid from the melt. For the Pb-Sn system, rotation substantially suppresses the onset of convection in the mushy zone and in the overlying liquid, holding open the promise that rotation can suppress freckling and other macrosegregation defects.

Pearlstein, Arne J.↗

The behavior of unsteady thermocapillary flows

The thermocapillary-driven flow of a liquid in a shallow slot provides a simple model to explore the effects of free-surface deflection and end walls on the stability of the flow. Such an investigation may be useful in understanding the complex flow seen in processes such as the float-zone refining of single crystals. A linear stability analysis of the viscously-dominated slot flow indicates stable basic states to both tow- and three-dimensional infinitesimal disturbances for capillary numbers below 43.3. Above this capillary number steady solutions do not exist. We also discuss current work involving the development of an inertially-dominated slot-flow model using singular perturbation methods. If these flows are unstable, inertial effects would be the only possible cause of the instability.

Vrane, David R.↗

Shear stabilization of solidification fronts

A linear-stability analysis is performed on the interface formed during the directional solidification of a dilute binary alloy in the presence of a time-periodic flow. In general, the flow is generated by the elliptical motion of the crystal parallel to the interface. The presence of the flow can either stabilize or destabilize the system relative to the case without flow with the result depending on the frequency and amplitude of the oscillations as well as the properties of the material. In the particular, we find that the morphological instability present in the absence of flow can be entirely suppressed with respect to disturbances of the same frequency as the oscillation.

Davis, Stephen H.↗

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗

Role of cross-field current instability in substorm onsets and intensifications

A cross field current instability is investigated as a potential mechanism for current reduction/disruption during substorm onsets and intensifications. Linear stability analysis shows that, for sufficiently strong current density, the instability can occur in the plasma sheet with a growth time comparable to the substorm onset time and excite waves with a significant electromagnetic component. Nonlinear analysis shows that the wave growth reaches the nonlinear stage in less than one ion gyroperiod and can reduce the cross tail current by approximately 15 to 28 of its initial value at saturation. The resulting anomalous resistivity is 11 to 12 orders of magnitude above the classical value. For a typical current reduction, the plasma in the disruption region is subjected to an earthward force. The substorm development scenario constructed based on this instability can readily account for a large number of substorm features.

Lui, Anthony T. Y.↗

Intraseasonal variability in a barotropic model with seasonal forcing

It has recently been suggested that oscillatory topographic instability could contribute to low-frequency variability over the Northern Hemisphere midlatitudes. A barotropic potential vorticity model, with a hierarchy of forcing and topography configurations on the sphere, is used to investigate the nature of low-frequency oscillations induced by such instabilities. Steady-state solutions of the model include multiple unstable equilibria that sustain oscillatory instabilities with periods of 10 - 15 days, and 150 - 180 days, for a realistic forcing pattern, Time-dependent solutions exhibit chaotic behavior with episodic oscillations, featuring both the intraseasonal (35 - 50 day) and biweekly (10 - 15 day) modes. The former is dominated by standing spatial patterns, the latter by traveling wave patterns. The phases of the intraseasonal oscillation are robust for all cases, exhibiting a clear oscillatory exchange of atmospheric angular momentum with the solid earth via mountain troque. It is demonstrated, through linear stability analysis on the sphere, that the intraseasonal oscillations are induced by topographic instabilities. The role of the seasonal cycle is studied by prescribing an annual cycle in the forcing. In this case, the winter forcing is more favorable than the summer for the occurrence of episodic intraseasonal oscillations. Recent observations are consistent with this model result.

Strong, Christopher M.↗

An instability in planetary rings due to ballistic transport

Ballistic transport in planetary rings is the net radial transport of mass and angular momentum due to exchanges of meteoroid impact ejecta between neighboring ring regions. The detailed linear stability analysis in this paper demonstrates that ballistic transport causes wavelike disturbances to grow and propagate in an otherwise uniform ring. The growth is strongest for intermediate values of the normal ring optical depth tau = 0.1 to 1.0 and goes to zero as tau approaches 0 and tau approaches infinity. For nominal values of various physical parameters, the minimum e-folding time is approximately 10(exp 5) years for tau approximately 0.4. The direction of propagation is opposite to the sense of any asymmetry that may exist in the ejecta direction distribution (inward for prograde ejecta and outward for retrograde ejecta). The additional effect of viscous transport tends to damp wavelike perturbations strongly at short wavelengths and at high values of tau. The quantitative agreement between this analytic work and numerical simulations reported elsewhere is generally quite good. As applied to Saturn's rings, the results in this paper strengthen the earlier conclusion from numerical calculations that the 100-km structure in the inner B Ring is caused by ballistic transport. However, it is also clear that ballistic transport cannot produce the complex structure seen in the outer two-thirds of the B Ring where tau greater than or approximately 1.5. Wavelike structures in the C Ring might also be attributed to ballistic transport; but this requires further study.

Durisen, Richard H.↗

Effects of Meteoroid Erosion in Planetary Rings

This grant supported continuing studies of the effects of ballistic transport on the evolution of Saturn's rings. Ballistic transport, as used in this context, refers to the net transport of mass and angular momentum caused by the exchange of meteoroid impact ejecta between neighboring ring regions (Ip 1983, 1984, Morfill et al. 1983, Lissauer 1984, Durisen 1984a,b). The characteristic time scale associated with this process is the gross erosion time t(sub g) the time it would take a ring region to be completed eroded if all impact ejecta were lost. This time scale is estimated to be about 10(exp 5) to 10(exp 6) years for a ring region with normal optical depth tau approximately 1. Earlier work by myself and collaborators developed the physical theory and simulation techniques to model this process (Durisen et al. 1989, Cuzzi and Durisen 1990). Detailed simulations, supported in part by this grant, have demonstrated that ballistic transport can produce observed structures in Saturn's rings, especially at and near the inner edges of the A and B Rings (Durisen et al. 1992, 1996). The structures of interest in the real rings are illustrated in Figures 1 and 2. Most of these structures were previously unexplained. The computational results plus analytic treatments place useful constraints on fundamental ring properties, including the indication of a relatively young ring age less than or equal 10(exp 8) years (see reviews by Nicholson and Dones 1991, Esposito 1993, Cuzzi 1995, and Porco 1995). This grant also supported development of the faster computational algorithms necessary to permit longer evolutions. Resulting simplications in the ballistic transport equations permitted an analytic linear stability analysis (Durisen 1995), which has provided considerable insight into ballistic transport processes and applications to Saturn's rings. All these accomplishments are described in more detail below.

Durisen, Richard H.↗

Hydrodynamic Stability of Multicomponent Droplet Gasification in Reduced Gravity

This investigation addresses the problem of hydrodynamic stability of a two-component droplet undergoing spherically-symmetrical gasification. The droplet components are assumed to have characteristic liquid species diffusion times that are large relative to characteristic droplet surface regression times. The problem is formulated as a linear stability analysis, with a goal of predicting when spherically-symmetric droplet gasification can be expected to be hydrodynamically unstable from surface-tension gradients acting along the surface of a droplet which result from perturbations. It is found that for the conditions assumed in this paper (quasisteady gas phase, no initial droplet temperature gradients, diffusion-dominated gasification), surface tension gradients do not play a role in the stability characteristics. In addition, all perturbations are predicted to decay such that droplets were hydrodynamically stable. Conditions are identified, however, that deserve more analysis as they may lead to hydrodynamic instabilities driven by capillary effects.

Aharon, I.↗