Stability of steady plane-parallel convective motion of a chemically active medium

A study is made of a vertical plane layer of reacting fluid whose boundaries are maintained at constant equal temperatures. As a result of heating due to a chemical reaction of zeroth order taking place in the fluid a steady plane-parallel convective flow develops in the layer, and if the internal heat release is sufficiently intense this can become unstable. The linear stability of this motion has hitherto been considered only in the hydro-dynamic formulation [1], in which one can ignore the thermal perturbations and their influence on the development of the hydrodynamic perturbations (the region of small Prandtl numbers). In the present paper, the stability boundary is determined for arbitrary values of the Prandtl number and the Frank-Kamenetskii parameter FK characterizing the steady plane-parallel regime. An important difference between this flow and the types of convective motion hitherto studied [2] is that the basic planeparallel flow of the reacting medium is possible only in a definite range of the parameter FK: At values of the parameter exceeding a critical value, there is a thermal explosion — abrupt strong heating of the fluid. This is due to the essentially nonlinear dependence of the heat release of a chemical reaction on the temperature.     

Secondary vibrational convective motions in a flat vertical layer of liquid

Investigations of the stability of steady-state plane-parallel convective motion between vertical planes heated to different temperatures [1–5] have shown that this motion, depending on the value of the Prandtl number P, exhibits instability of two types. With small and moderate Prandtl numbers, the instability is of a hydrodynamic nature. It is brought about by monotonic perturbations which, in the supercritical region, develop into a periodic, with respect to the vertical, system of steady-state vortices at the interface between the opposing convective flows. Articles [6, 7] are devoted to the numerical investigation of nonlinear secondary steady-state flows. If P>11.4, there appears a new mode of instability, i.e., running thermal waves increasing in the flow; with P>12, this mode becomes more dangerous [4]. This instability is connected with the development of vibrational perturbations, and it can be considered that in the supercritical region the perturbations lead to the establishment of steady-state vibrations. Linear theory has made it possible to determine the boundaries of the regions of stability. In the present article a numerical investigation is made of nonlinear supercritical conditions developing as a result of a loss of stability of the steady-state flow with respect to vibrational perturbations.     

Stability of convective flow in a vibration field with respect to three-dimensional perturbations

A plane-parallel convective flow in a vertical layer between boundaries maintained at different temperatures becomes unstable when the Grashof number reaches a critical value (see [1]). In [2, 3] the effect of high-frequency harmonic vibration in the vertical direction on the stability of this flow was investigated. The presence of vibration in a nonisothermal fluid leads to the appearance of a new instability mechanism which operates even under conditions of total weightlessness [4]. As shown in [2, 3], the interaction of the usual instability mechanisms in a static gravity field and the vibration mechanism has an important influence on the stability of convective flow. In this paper the flow stability is considered for an arbitrary direction of the vibration axis in the plane of the layer and the stability characteristics with respect to three-dimensional normal perturbations are determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 116–122, March–April, 1988.     

Vibrations of a body in a bounded volume of viscous fluid

The problem of the vibrations of a body in a bounded volume of viscous fluid has been studied on a number of occasions [1–4]. The main attention has been devoted to determining the hydrodynamic characteristics of elements in the form of rods. Analytic solution of the problem is possible only in the simplest cases [2]. In the present paper, in which large Reynolds numbers are considered, the asymptotic method of Vishik and Lyusternik [5] and Chernous' ko [6] is used to consider the general problem of translational vibrations of an axisymmetric body in an axisymmetric volume of fluid. Equations of motion of the body and expressions for the coefficients due to the viscosity of the fluid are obtained. It is shown that in the first approximation these coefficients differ only by a constant factor and are completely determined if the solution to the problem for an ideal fluid is known. Examples are given of the determination of the “viscous” added mass and the damping coefficient for some bodies and cavities. In the case of an ideal fluid, general estimates are obtained for the added mass and also for the influence of nonlinearity. Ritz's method is used to solve the problem of longitudinal vibrations of an ellipsoid of revolution in a circular cylinder. The hydrodynamic coefficients have been determined numerically on a computer. The theoretical results agree well with the results of experimental investigations.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

Motion of a solidcontaininga spherical damper

For the purpose of modeling the motion of a solid with a cavity filled with a viscous fluid, M. A. Lavrent'ev [1] has proposed a model in the form of a solid with a spherical cavity in which another solid, spherical in shape, is enclosed. The sphere is separated from the cavity walls by a small, clearance in which viscous forces act (a lubricating film). This simple model with a finite number of degrees of freedom possesses certain mechanical properties of a solid with a cavity containing a viscous fluid. Study of this model is therefore of interest.The present paper examines certain properties of the model, which will be termed a solid with a damper. It is shown that for a highviscosity lubricant the motion of a solid with a damper can be described by the same equations which pertain to the motion of a solid with a spherical cavity filled with a high-viscosity fluid. Expressions relating the parameters of the systems are obtained. If these relations are fulfilled, the systems become mechanically equivalent.The steady motions of a free solid with a damper and their stability conditions are determined.These motions and stability conditions hold for a body with a cavity filled with a viscous fluid [2].     

Angular vibrations of an ellipsoid of revolution in a confined viscous fluid

Vibrations of bodies in confined viscous fluids have been studied on a number of occasions, transverse vibrations of rods being the main subject of investigation [1–3]. The present authors [4] have considered the general problem of translational vibrations of an axisymmetric body in an axisymmetric region containing a low-viscosity fluid. The present paper follows the same approach and considers the problem of small angular vibrations of an ellipsoid of revolution in a circular cylinder with flat ends. In the general case, the hydrodynamic coefficients in the equation of motion of the ellipsoid are determined numerically for different values of the dimensionless geometrical parameters using Ritz's method. In the case of an unconfined fluid, analytic dependences in terms of elementary functions are obtained for the hydrodynamic coefficients. The theoretical results agree well with experimental investigations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 34–39, March–April, 1984.     

Control of Thermo- and Solutocapillary Flows in FZ Crystal Growth by High-Frequency Vibrations

The paper deals with the numerical investigation of the possibilities to control convective flows in the liquid bridge in zero gravity conditions applying axial vibrations. The surface tension is assumed to be dependent both on the temperature and on the solute concentration. The free surface deformations and the curvature of the phase change surfaces are neglected but pulsational deformations of the free surface are accounted for. The first part of the paper concerns axisymmetric steady flows. The calculations show that the evolution of convective flow with the variation of thermal Marangoni number at a fixed value of the solutal Marangoni number is accompanied by the hysteresis phenomenon, which is related to the existence of two stable steady regimes in a certain parameter range. One of these regimes is thermocapillary dominated, it corresponds to the two-vortex flow, and the other is solutocapillary dominated, it corresponds to the single-vortex flow. Under vibrations, the range of the Marangoni numbers where hysteresis is observed becomes narrower and is shifted to the area of larger values. The second part of the paper concerns the stability of axisymmetric thermo-and solutocapillary flows and the transition to three-dimensional regimes. Significant mutual influence of flows generated by each process on the stability of the other is discovered. Stability maps in the parametric plane for the thermal Marangoni number, the solutal Marangoni number, are obtained for different values of vibration parameters. It is shown, that vibrations exert a stabilizing effect, increasing critical Marangoni numbers for all modes of instability. However, this effect is different for different modes and at high vibration intensity destabilization is possible. Consequently, vibrations can modify the scenario of the transition to the three-dimensional mode.     

Development of a steady relief at the interface of fluids in a vibrational field

The vibrations of a vessel strongly influence the behavior of the interface of the fluids in it. Thus, vertical vibrations can lead both to the parametric excitation of waves (Faraday ripples) and to the suppression of the Rayleigh-Taylor instability [1–2]. At the present time, the influence of vertical vibrations on the behavior of a fluid surface have been studied in sufficient detail (see, for example, review [3]). The behavior of an interface of fluids in the case of horizontal vibrations has been studied less. An interesting phenomenon has been revealed in the experimental papers [4, 5]: in the case of fairly strong horizontal vibrations of a vessel containing a fluid with a free surface, the fluid collects near one of the vertical vessel walls, the free surface being practically plane and stationary with respect to the vessel, while its angle of inclination to the horizon depends on the vibration rate. But if there is a system of immiscible fluids with comparable but different densities in the vessel, horizontal vibrations lead to the formation of a steady wave relief at the interface. An explanation of the behavior of a fluid with a free boundary was given in [6] on the basis of averaged equations of fluid motion in a vibrational field. The present paper is devoted to an analysis of the behavior of the interface of fluids with comparable densities in a high-frequency vibrational field. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 8–13, November–December, 1986.     

Development and instability of steady convective flows in a square cavity heated from below and a field of vertically directed vibrational forces

The present paper is devoted to numerical investigation of the spatial structure and stability of secondary vibrational convective flows resulting from instability of the equilibrium of a fluid heated from below. Vibrations parallel to the vector of the gravitational force (vertical vibrations) are considered. As in earlier work [7–9], a region of finite size is used — a square cavity heated from below. It is shown that enhancement of the vibrational disturbance of the natural convective flow may either stabilize or destabilize flows with different spatial structures; it may also stabilize certain solutions of the system of convection equations that are unstable in the absence of vibrational forces. In addition, increase of the vibrational Rayleigh number can lead to a change of the mechanisms responsible for equilibrium instability and oscillatory instability of the secondary steady flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–18, March–April, 1991.I thank G. Z. Gershuni for assistance and extremely fruitful discussions of the results of the paper.     

Vibrational-convective instability of a horizontal fluid layer with internal heat sources

The convective motion of a nonisothermal fluid in a gravity field in a vibrating cavity is caused by two mechanisms: the usual static mechanism and a vibrational mechanism. The same mechanisms are also responsible for mechanical equilibrium crisis under the conditions in which such equilibrium is possible. The research on these questions is reviewed in [1]. The problems of vibrational-convective stability examined so far relate to cases in which the nonisothermicity was created by specifying the temperature at the boundaries of the region. The present study is concerned with the vibrational-convective stability of a fluid in which the temperature nonuniformity is created by internal heat generation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 3–7, September–October, 1985.     

Convective Instability in Slip Flow in a Vertical Circular Porous Microchannel

The paper represents an analysis of convective instability in a vertical cylindrical porous microchannel performed using the Galerkin method. The dependence of the critical Rayleigh number on the Darcy, Knudsen, and Prandtl numbers, as well as on the ratio of the thermal conductivities of the fluid and the wall, was obtained. It was shown that a decrease in permeability of the porous medium (in other words, increase in its porosity) causes an increase in flow stability. This effect is substantially nonlinear. Under the condition Da?>?0.1, the effect of the porosity on the critical Rayleigh number practically vanishes. Strengthening of the slippage effects leads to an increase in the instability of the entire system. The slippage effect on the critical Rayleigh number is nonlinear. The level of nonlinearity depends on the Prandtl number. With an increase in the Prandtl number, the effect of slippage on the onset of convection weakens. With an increase in the ratio of the thermal conductivities of the fluid and the wall, the influence of the Prandtl number decreases. At high values of the Prandtl numbers (Pr?>?10), its influence practically vanishes.

     

The effect of vertical vibrations on the onset of convection in a planar rotating layer

The effect of vertical vibrations on the convection in a rotating planar fluid layer heated from below was studied. In this case a modulation parameter, the acceleration due to gravity, appears in the problem. The modulation of the parameter may have a significant effect on the onset of convective instability. Parameter modulation in nonrotating layers has been investigated in earlier work [1–3]. The presence of rotation significantly increases the complexity of the mathematical problem, introducing an additional dependence of the solution on the Taylor number Ta and the Prandtl number Pr. Furthermore, an oscillatory convection regime can occur at the stability limit in rotating fluids with Pr < 1. Parameter modulation in the rotating fluid may not only lead to a change in the stability limit and critical wavelength but also to a change in the eigenfrequency of the oscillatory convection. Rauscher and Kelly [4] examined the effect of parameter modulation on the convective stability of a rotating fluid only for the particular case of a sinusoidal variation in the temperature gradient with a small amplitude for Pr = 1, i.e., the effect of modulation was studied on only a steady convection regime.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 12–22, July–August, 1984.     

Stability of Thermosolutal Natural Convection in Superposed Fluid and Porous Layers

This article deals with the onset of thermosolutal natural convection in horizontal superposed fluid and porous layers. A linear stability analysis is performed using the one-domain approach. As in the thermal convection case, the results show a bimodal nature of the marginal stability curves where each mode corresponds to a different convective instability. At small wave numbers, the convective flow occurs in the whole cavity (“porous mode”) while perturbations of large wave numbers lead to a convective flow mainly confined in the fluid layer (“fluid mode”). Furthermore, it is shown that the onset of thermosolutal natural convection is characterized by a multi-cellular flow in the fluid region for negative thermal Rayleigh numbers. For positive thermal Rayleigh numbers, the convective flow takes place both in the fluid and porous regions. The influence of the depth ratio and thermal diffusivity ratio is also investigated for a wide range of the thermal Rayleigh numbers.     

Oscillations of a body filled with a viscous fluid

The oscillations of a physical pendulum containing a spherical cavity filled with an incompressible viscous liquid were discussed in [1]. In this paper we consider the mote general problem of the motion of an axially symmetric solid with a spherical cavity filled with an incompressible viscous fluid and moving about a fixed point. It is assumed that the center of the cavity and the fixed point lie on the axis of symmetry of the body.     

Steady convective motions in a plane horizontal fluid layer with permeable boundaries

In a plane horizontal fluid layer bounded by permeable plane surfaces which are heated to different temperatures and between which transverse flow takes place with uniform velocity, convection occurs at a definite critical Rayieigh number. The study of the disturbance spectrum and the convective stability, made within the framework of linear theory in [1], showed that convective instability in the layer with permeable boundaries, just as in the case of the Rayieigh problem, is associated with the development of monotonie disturbances. It turns out that the transverse motion in the layer leads to a considerable increase of the Rayieigh number. Linear theory does not permit analysis of the development of the disturbances in the supercritical region. Analysis of the developed nonlinear motion can be made only on the basis of the complete nonlinear convection equations.In this investigation we made a numerical study of nonlinear motions in the supercritical region. Calculations were made on a computer via the grid method. Solutions are obtained for the nonlinear equations of motion over a wide range of Rayieigh numbers for different values of the Peclet number, whichdefines the intensity of the transverse motion in the layer.The author wishes to thank E. M. Zhukovitskii for his guidance, and G. Z. Gershuni and E. L. Tarunin for their interest and assistance in the study.     

Convection in rotating spherical layers with allowance for centrifugal forces

The investigation of convection in rotating spherical layers with a central gravitational field g(r) is very important for the study of the global motions in the atmospheres of large planets and the convective zones of stars. In recent years, many studies of these questions have been made (they have been reviewed, for example, by Yavorskaya and Belyaev [1]), but the centrifugal convective force has been ignored in all the numerical and analytic investigations. In some cases, for example, for large planets, the centrifugal force may reach an appreciable value, O.1g, and have a strong influence on the convective motion. The present paper studies the occurrence of convection in slowly rotating spherical layers with allowance for centrifugal forces. It is shown that the centrifugal force leads to the appearance in a layer of an axisymmetric flow, at the stability limit of which convective cells of banana or toroidal shape can develop. The latter are possible only in layers with undeformable boundaries at sufficiently large values of the Froude number. Irrespective of the form of the layer and the magnitude of the centrifugal force, the banana-shaped cells propagate in a wavelike manner in the direction opposite to the rotation. In the case of undeformable boundaries, the centrifugal force stabilizes the motion of the fluid as compared with the case of a layer at rest. Deformation of one or both of the boundaries under the influence of the centrifugal force leads to destabilization of the basic flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 14–21, March–April, 1984.     

Influence of temperature dependent conductivity of a nanofluid in a vertical rectangular duct

Natural convective heat transfer and fluid flow in a vertical rectangular duct filled with a nanofluid is studied numerically assuming the thermal conductivity to be dependent on the fluid temperature. The transport equations for mass, momentum and energy formulated in dimensionless form are solved numerically using finite difference method. Particular efforts have been focused on the effects of the thermal conductivity variation parameter, Grashof number, Brinkman number, nanoparticles volume fraction, aspect ratio and type of nanoparticles on the fluid flow and heat transfer inside the cavity. It is found that the flow was enhanced for the increase in Grashof number, Brinkman number and aspect ratio for any values of conductivity variation parameter and for regular fluid and nanofluid. The heat transfer rate for regular fluid is less than that for the nanofluid for all governing parameters.     

Stability of stationary convective flow of a mixture in a vertical porous layer

In contrast to the corresponding viscous flow, the convective flow of a homogeneous liquid in a planar vertical layer whose boundaries are maintained at different temperatures is stable [1]. When a porous layer is saturated with a binary mixture, in the presence of potentially stable stratification one must expect an instability of thermal-concentration nature to be manifested. This instability mechanism is associated with the difference between the temperature and concentration relaxation times, which leads to a buoyancy force when an element of the fluid is displaced horizontally. In viscous binary mixtures, the thermal-concentration instability is the origin of the formation of layered flows, which have been studied in detail in recent years [2–4]. The convective instability of the equilibrium of a binary mixture in a porous medium was considered earlier by the present authors in [5]. In the present paper, the stability of stationary convective flow of a binary mixture in a planar vertical porous layer is studied. It is shown that in the presence of sufficient longitudinal stratification the flow becomes unstable against thermal-concentration perturbations; the stability boundary is determined as a function of the parameters of the problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 150–157, January–February, 1980.     

Stability of plane poiseuille flow of a viscous anisotropic fluid

Experimental investigations show that the presence in a fluid of fibers and rigid asymmetric particles leads to a greater stability of flow in tubes and lowers the turbulent frictional resistance in a certain range of Reynolds numbers [1]. In the present paper, the anisotropic structure of a fluid with additives is described by Ericksen's rheological model [2]. The parameters of the model are particularized in accordance with the paper [3] of Pilipenko, Kalinichenko, and Lemak, and in the limiting case of weak Brownian motion allowance is made for the effect of the predominant orientation of the particles and the influence of additives on the longitudinal and shear viscosity. The stability of the Poiseuille flow is considered in the linear formulation. In an anisotropic viscous fluid, an equation of Orr-Sommerfeld type has a singular point. A rule for choosing the path of integration avoiding the singular point is obtained on the basis of a generalization of the method of Dikii [4] proposed in an investigation of the stability of the flow of an ideal fluid. The results of numerical calculations of the neutral stability curve for two-dimensional perturbations are given.