Convective flow with uniform vorticity in a vertical layer

The free convective circulation of liquid in plane vertical slits of circular and square cross section with a longitudinal horizontal temperature gradient at the boundaries was investigated experimentally. It was found that under such heating conditions there is a uniform-vorticity flow with a region of quasirigid rotation, which has the shape of a disk in a circular slit and the shape of a cross in a square slit; in each longitudinal section of this zone the liquid moves along concentric trajectories with constant angular velocity. Dimensionless numbers for the problem were established by tests with various liquids and cavities of different dimensions. In dimensionless form, the angular velocity of the vortex and the temperature gradient in it depend linearly on the temperature difference at the boundaries of the layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 160–165, May–June, 1984.     

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.     

Convective instability of a fluid in a vibration field under conditions of weightlessness

The problem of the convection and convective instability of a fluid in a high-frequency vibration field under conditions of weightlessness was formulated in an earlier paper of the authors [1]. In the present paper, the conditions of equilibrium are discussed and the boundaries of vibration instability are determined for some equilibrium states: a plane layer of fluid with transverse temperature gradient and arbitrary direction of the vibration, a cylindrical layer with radial gradient and longitudinal direction of the vibration, and an infinite circular cylinder with transverse and mutually perpendicular directions of the temperature gradient and the vibration axis.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 12–19, July–August, 1981.We thank G. I. Petrov for helpful discussions.     

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.     

Effect of thermal radiation exchange between boundaries on convection in a gas heated from below

Natural convection problems offer many examples of branching of the solutions [1]. Usually, such branching (from the standpoint of catastrophe theory) can be described by a Whitney fold or cusp. A characteristic feature of nontrivial branching is the presence of some small but finite disturbance of the convective equilibrium conditions. In this study the perturbation disturbing the convective equilibrium of a fluid heated from below is Stefan-law thermal radiation exchange between the boundaries of the enclosure. Natural convection with lateral heating and allowance for radiative heat transfer was previously investigated in [2].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 47–51, September–October, 1992.     

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.     

Thermal Convection in a Plane Layer Rotating About a Horizontal Axis

The thermal convection of a fluid in a plane vertical layer with a cylindrical lateral boundary, which rotates uniformly about a horizontal symmetry axis, is investigated experimentally. The structure and excitation limit of the convective flows are studied as functions of the rotation frequency, the temperature difference between the layer boundaries, and the layer thickness. The determining dimensionless parameters are found. It is shown that the period-average gravity action produces convection in the form of hexagonally distributed cells stationary in the reference system tied to the cavity.     

Numerical investigation of convective motion in a closed cavity

The investigation of thermal convection in a closed cavity is of considerable interest in connection with the problem of heat transfer. The problem may be solved comparatively simply in the case of small characteristic temperature difference with heating from the side, when equilibrium is not possible and when slow movement is initiated for an arbitrarily small horizontal temperature gradient. In this case the motion may be studied using the small parameter method, based on expanding the velocity, temperature, and pressure in series in powers of the Grashof number—the dimensionless parameter which characterizes the intensity of the convection [1–4]. In the problems considered it has been possible to find only two or three terms of these series. The solutions obtained in this approximation describe only weak nonlinear effects and the region of their applicability is limited, naturally, to small values of the Grashof number (no larger than 10³).With increase of the temperature difference the nature of the motion gradually changes—at the boundaries of the cavity a convective boundary layer is formed, in which the primary temperature and velocity gradients are concentrated; the remaining portion of the liquid forms the flow core. On the basis of an analysis of the equations of motion for the plane case, Batchelor [4] suggested that the core is isothermal and rotates with constant and uniform vorticity. The value of the vorticity in the core must be determined as the eigenvalue of the problem of a closed boundary layer. A closed convective boundary layer in a horizontal cylinder and in a plane vertical stratum was considered in [5, 6] using the Batchelor scheme. The boundary layer parameters and the vorticity in the core were determined with the aid of an integral method. An attempt to solve the boundary layer equations analytically for a horizontal cylinder using the Oseen linearization method was made in [7].However, the results of experiments in which a study was made of the structure of the convective motion of various liquids and gases in closed cavities of different shapes [8–13] definitely contradict the Batchelor hypothesis. The measurements show that the core is not isothermal; on the contrary, there is a constant vertical temperature gradient directed upward in the core. Further, the core is practically motionless. In the core there are found retrograde motions with velocities much smaller than the velocities in the boundary layer.The use of numerical methods may be of assistance in clarifying the laws governing the convective motion in a closed cavity with large temperature differences. In [14] the two-dimensional problem of steady air convection in a square cavity was solved by expansion in orthogonal polynomials. The author was able to progress in the calculation only to a value of the Grashof numberG=10⁴. At these values of the Grashof numberG the formation of the boundary layer and the core has really only started, therefore the author's conclusion on the agreement of the numerical results with the Batchelor hypothesis is not justified. In addition, the bifurcation of the central isotherm (Fig. 3 of [14]), on the basis of which the conclusion was drawn concerning the formation of the isothermal core, is apparently the result of a misunderstanding, since an isotherm of this form obviously contradicts the symmetry of the solution.In [5] the method of finite differences is used to obtain the solution of the problem of strong convection of a gas in a horizontal cylinder whose lateral sides have different temperatures. According to the results of the calculation and in accordance with the experimental data [9], in the cavity there is a practically stationary core. However, since the authors started from the convection equations in the boundary layer approximation they did not obtain any detailed information on the core structure, in particular on the distribution of the temperature in the core.In the following we present the results of a finite difference solution of the complete nonlinear problem of plane convective motion in a square cavity. The vertical boundaries of the cavity are held at constant temperatures; the temperature varies linearly on the horizontal boundaries. The velocity and temperature distributions are obtained for values of the Grashof number in the range 0<G4·10⁵ and for a value of the Prandtl number P=1. The results of the calculation permit following the formation of the closed boundary layer and the very slowly moving core with a constant vertical temperature gradient. The heat flux through the cavity is found as a function of the Grashof number.     

Secondary convective motions in a plane inclined layer

A linear theory of stability of a plane-parallel convective flow between infinite isothermal planes heated to different temperature was developed in [1–6]. At moderate Pr values the instability is monotonic and leads to the development of steady secondary motions. These motions for the case of a vertical layer have been investigated by the net [7, 8] and small-parameter [9] methods. In this paper steady secondary motions in an inclined layer are investigated. The small-parameter and net methods are used. The hard nature of excitation of secondary motions in a defined range of tilt angles is established. There are two types of secondary motions, whose regions of existence overlap — vortices at the boundary of countercurrent streams and convection rolls; the hard instability is due to the development of convection rolls. The analog of the Squire transformation obtained in [4] for infinitely small disturbances of a plane-parallel convective flow is extended to secondary motions of finite amplitude.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–9, May–June, 1977.I thank G. Z. Gershumi, E. M. Zhukhovitskii, and E. L. Tarunin for interest in the work and valuable discussion.     

Thermal structure of convective vortex lattices

The thermal structure of the convective motions of a rotating plane layer of fluid is experimentally investigated in the regular vortex structure regime. It is found that in such a system the intense vortex motion leads to a temperature distribution such that the mean fluid temperature falls linearly from the bottom of the layer to the surface, the temperature gradient being determined by the rate of rotation and depth of the fluid. By dimensional analysis it is shown that this gradient corresponds to heat transfer in which the Nusselt number isolines are parallel to the convection curve. The horizontal structure of the temperature field is investigated; it corresponds to motion in which the fluid descends within a narrow vortex-sink and rises along the edges of a cylinder which determines the characteristic dimension of the structure in rotating fluid convection.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 160–166, November–December, 1987.The author wishes to thank G. S. Golitsyn for his constant interest in the work.     

Non-Darcian Effects on the Onset of Convection in a Porous Layer with Throughflow

The Brinkman extended Darcy model including Lapwood and Forchheimer inertia terms with fluid viscosity being different from effective viscosity is employed to investigate the effect of vertical throughflow on thermal convective instabilities in a porous layer. Three different types of boundary conditions (free–free, rigid–rigid and rigid–free) are considered which are either conducting or insulating to temperature perturbations. The Galerkin method is used to calculate the critical Rayleigh numbers for conducting boundaries, while closed form solutions are achieved for insulating boundaries. The relative importance of inertial resistance on convective instabilities is investigated in detail. In the case of rigid–free boundaries, it is found that throughflow is destabilizing depending on the choice of physical parameters and the model used. Further, it is noted that an increase in viscosity ratio delays the onset of convection. Standard results are also obtained as particular cases from the general model presented here.     

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.     

Intrapore convection when the average temperature gradient is vertical

The heat conduction of a porous medium saturated with a fluid is usually regarded as being purely molecular [1]. The assumption here is that in the case of heating from below the local temperature gradient within each of the pores, like the averaged gradient in the complete layer, is strictly vertical, and, since the pores are as a rule small, this local gradient is less than the critical. It is therefore assumed that in the absence of large-scale convection the fluid in the pores is in equilibrium. However, for different thermal conductivities of the fluid and the porous skeleton surrounding it a vertical temperature gradient in the fluid and, accordingly, equilibrium of the fluid are possible only if a cavity is a sphere or an ellipsoid with a definite orientation [1]. Since the pores do not have such shapes, the convective motion that arises in each of the pores or in several communicating pores can lead to an increase in the effective thermal conductivity of the fluid and, accordingly, the effective thermal conductivity of the complete medium. The present paper is devoted to study of this effect.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 93–98, January–February, 1984.     

Convection and heat transfer in a vertical layer with allowance for radiation from nonisothermic walls

Many studies (for example, [1–5]) consider motion and heat transfer in closed vertical cavities with given different temperatures of the lateral boundaries. The majority of studies cover the case of convection, but of late studies have appeared (for example, [4]) in which joint radiative—convective heat transfer is taken into account. In the present study we consider motion and heat transfer in a rectangular cavity separating two media with given different temperatures. In contrast to [4], the temperature of the lateral boundaries is determined from the condition for interaction with the surrounding medium, and the air in the cavity is assumed to be transparent for the heat radiation of the walls. The problem considered is a mathematical model of the heat transfer through windows, and is necessary for the analysis of methods of improving the heat proofing of buildings.Translated-from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 25–30, 1987.     

Finite-amplitude regimes of mixed convection in a vertical layer with undulating boundaries

The method of finite differences is used to construct convective motions in a vertical layer with sinusoidally curved boundaries, fluid being pumped through longitudinally. Apart from steady and oscillation regimes, found earlier by analytical means for small amplitudes of undulation and slow pumping through [1, 2], new, essentially nonlinear, types of motion are discovered in the form of two-stroke cycles, and also of complex multi-revolution cycles which are two-dimensional resonance tori. The regions are determined in which regimes of various types exist.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 16–20, January–February, 1987.The author is grateful to E. M. Zhukhovitskii for constant interest in the study, and also to V. S. Anishchenko and A. A. Nepomnyashchii for useful discussions.     

Numerical modeling of thermocapillary convection in a fluid layer with local heating of the free surface

Thermocapillary convection in a plane horizontal fluid layer with concentrated heating of the free surface is modeled numerically using the Navier-Stokes equations and the heat transport equation. This makes it possible to examine the structure of the convection throughout the fluid volume, in particular in the region where the motion is weak. The deformation of the free surface is assumed to be negligibly small. In the case of a ponderable fluid this assumption is justified given certain upper and lower constraints on the temperature difference and the thickness of the layer, respectively, [9, 10]. Under conditions of weightlessness a fluid layer of constant thickness in a rectangular channel can be realized at a contact angle of 90° [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 108–113, July–August, 1987.     

The stability of a flow in a side-heated vertical layer with internal heat sources

The linear stability of a combined convective flow in a side-heated plane vertical layer of incompressible fluid is considered. The layer contains uniformly distributed heat sources. The steadystate flow is a superposition of the Gershuni convective flow and a flow induced by the action of the internal heat sources. The stability maps of the combined flow are constructed for different values of the Prandtl number. It is found that the internal sources damp the action of the mechanisms of the Gershuni flow crisis. The lateral heating may result in both the stabilization and the destabilization of the flow caused by the action of internal heat sources.     

Convective Instability in Superposed Fluid and Porous Layers with Vertical Throughflow

A closed form solution to the convective instability in a composite system of fluid and porous layers with vertical throughflow is presented. The boundaries are considered to be rigid-permeable and insulating to temperature perturbations. Flow in the porous layer is governed by Darcy–Forchheimer equation and the Beavers–Joseph condition is applied at the interface between the fluid and the porous layer. In contrast to the single-layer system, it is found that destabilization due to throughflow arises, and the ratio of fluid layer thickness to porous layer thickness, , too, plays a crucial role in deciding the stability of the system depending on the Prandtl number.     

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.     

Unsteady thermocapillary convection in a nonuniformly heated fluid layer

In many technological processes, thin extended layers of nonuniformly heated fluid are used [1–3]. If they are sufficiently thin, thermocapillary forces have a decisive influence on the occurrence and development of motion of the fluid [4–6]. Investigation of convective motion in such a layer is of great interest for estimating the intensity of heat and mass transfer in technological processes. This paper is a study of unsteady thermocapillary motion in a layer of viscous incompressible fluid with free surface in which a thermal inhomogeneity is created at the initial time. Approximate expressions are obtained for the fields of the velocity, temperature, and pressure in the fluid, and also for the shape of the free surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 17–25, May–June, 1991.