Search NASA⌕ Search

SEARCH · Search NASA

Results for “buoyancy-driven instability”

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.

Buoyancy-driven instability in a vertical cylinder: Binary fluids with Soret effect. I - General theory and stationary stability results

The buoyancy-driven instability of a monocomponent or binary fluid completely contained in a vertical circular cylinder is investigated, including the influence of the Soret effect for the binary mixture. The Boussinesq approximation is used, and the resulting linear stability problem is solved using a Galerkin technique. The analysis considers fluid mixtures ranging from gases to liquid metals. The flow structure is found to depend strongly on both the cylinder aspect ratio and the magnitude of the Soret effect. The predicted stability limits are shown to agree closely with experimental observations.

Hardin, G. R.↗

Buoyancy-driven instabilities and the nonlinear breakup of a sheared magnetic layer

Motivated by problems concerning the storage and subsequent escape of the solar magnetic field, a study has been made of how a magnetic layer embedded in a convectively stable atmosphere evolves due to axisymmetric instabilities driven by magnetic buoyancy. The initial equilibrium consists of a toroidal field sheared by a weaker poloidal component. The linear stability problem is investigated for both ideal and resistive MHD, and the nonlinear evolution is followed by numerical integration of the equations of motion. In all cases, the instability is greatly affected by the distribution and strength of the poloidal field. In particular, both the horizontal and vertical scales of the motions are controlled by the location of the surface on which the poloidal field vanishes: the resonant surface. In the nonlinear regime, a resonant surface close to the interface between the magnetized and field-free fluid leads to the localization of the instability, so that only a fraction of the magnetic region is disrupted by the motions. By contrast, a deeply seated resonant surface leads to the complete disruption of the layer and to the formation of large, helical magnetic fragments whose identity is preserved for the entire simulation.

Cattaneo, Fausto↗

Buoyancy-Driven Instabilities in Single-Bubble Sonoluminescence

The principal objectives of this study are to determine how gravity affects the emission of light from single-bubble sonoluminescence (SBSL), and whether or not the bubble extinction is directly related to gravity. Our experimental task involves designing glass or quartz spherical levitation cells that generate very stable SL bubbles. The cells must have minimized vibration, and some temperature control. The experimental system will reside in a light-tight enclosure. Aside from acceleration, the frequency, pressure amplitude, and light intensity must be measured. A computer program will be constructed to perform all aspects of the experiment.

Matula, Thomas J.↗

Asymmetric instabilities in buoyancy-driven flow in a tall vertical annulus

Linear stability of the one-dimensional flow between infinite vertical coaxial cylinders induced by heating the inner cylinder is considered for various ratios (kappa) of the inner radius to the outer radius, and for Prandtl numbers (P) appropriate to air and water. For air with P = 0.71 the least stable disturbance is nonaxisymmetric for kappa less than 0.44 and is axisymmetric for kappa greater than 0.44, and in either case the instability is due to the action of the shear forces. For P = 3.5, the situation is similar, except that the asymmetric shear mode is superseded by an axisymmetric instability driven by buoyancy forces for kappa = 0.03-0.16. Wave speeds, wavenumbers, and critical Grashof numbers for these cases and for the case of zero Prandtl number are given.

Mcfadden, G. B.↗

The nonlinear breakup of the sun's toroidal field

There are good reasons for believing that the sun has a strong toroidal magnetic field in the stably stratified region of convective overshoot sandwiched between the radiative zone and convective zone proper. The magnetic field in this region is modeled by studying the behavior of a layer of uniform field embedded in a subadiabatic atmosphere. Since the field can support extra mass, such a configuration is top-heavy, and instabilities of the Rayleigh-Taylor type can occur. Numerical integration of the two-dimensional compressible MHD equations makes it possible to follow the evolution of this instability into the nonlinear regime. The initial buoyancy-driven instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The somewhat surprising subsequent motions are determined primarily by the strong interactions between vortices.

Hughes, D. W.↗

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

Hydrodynamic interactions in directional solidification

A binary liquid undergoes unidirectional solidification. The one-dimensional steady state is susceptible to morphological instability that causes the solid/liquid interface to change from a planar state to a cellular pattern. This paper examines the effects on this transition of volume-change convection, buoyancy-driven convection or forced flows. It emphasizes how flows alter stability limits, create scale and pattern changes in morphology, and create, through coupling, new instabilities. Emphasis is placed on the physical mechanisms of the interactions.

Davis, Stephen H.↗

Frontal Polymerization in Microgravity

Frontal polymerization systems, with their inherent large thermal and compositional gradients, are greatly affected by buoyancy-driven convection. Sounding rocket experiments allowed the preparation of benchmark materials and demonstrated that methods to suppress the Rayleigh-Taylor instability in ground-based research did not significantly affect the molecular weight of the polymer. Experiments under weightlessness show clearly that bubbles produced during the reaction interact very differently than under 1 g.

Pojman, John A.↗

The effect of gravitational modulation on convection in vertical Bridgman growth

During vertical directional solidification of a binary alloy at constant velocity, buoyancy-driven solutal convection may occur due to the solute gradient associated with the solidification process. This problem is further complicated if time-periodic forcing is considered, which is relevant to materials processing in a microgravity environment or as a means of dynamic control of flow instabilities. The effect of time-periodic modulation is studied by introducing a gravitational acceleration which is a sinusoidal function of time. The onset of solutal convection is treated by a stability analysis of the linearized governing equations and boundary conditions. Solutions are obtained numerically employing two distinct computational implementations of Floquet theory. Results are presented for materials with large Schmidt number, and an analysis for large frequency and large Schmidt number yields a more complete description of the behavior in this relevant limit.

Murray, B. T.↗

Morphological stability in the presence of fluid flow in the melt

Recent experiments have shown that the presence of a vertical buoyancy-driven flow adjacent to an initially cylindrical crystal-melt interface may produce a time-dependent helical deformation of the interface, with a rotation period ranging from several minutes to many hours, depending upon the width of the melt. The temperature distribution is such that the interface is expected to be morphologically stable in the absence of fluid flow. A linear stability analysis reveals that the instability is due to a coupling between a basic hydrodynamic instability in the buoyant flow and the deformable boundary separating the two phases. The crystal-melt interface lowers the critical Grashof number of an analogous rigid-walled system by an order of magnitude for succinonitrile with a Prandtl number P = 22.8; furthermore, the hydrodynamic mode that is actually destabilized by the interface is not the least stable mode in the rigid-walled system for P = 22.8. The results show that the instability may be regarded either as a rather large alteration of a basic hydrodynamic instability by the crystal-melt interface, or as a significant modification of the morphological stability of the interface by the presence of the buoyant flow.

Mcfadden, G. B.↗

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗

Experimental Studies of the Interaction Between a Parallel Shear Flow and a Directionally-Solidifying Front

It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.

Zhang, Meng↗

Gravitational Effects on Flow Instability and Transition in Low Density Jets

Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Agrawal A. K.↗

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↗