Convection in a closed cavity heated from below when equilibrium conditions are disrupted

The equilibrium of a fluid is possible in a closed cavity in the presence of a strictly vertical temperature gradient (heating from below) [1]. There is a distinct sequence of critical Rayleigh numbers R_i at which this equilibrium loses its stability relative to low characteristic perturbations. The presence of different finite perturbations, unavoidable in an experiment, leads to the absence of a strict equilibrium when R < R₁. The problem of the influence of the perturbation on the convection conditions near the critical points arises in this context [2, 3]. The case in which the cavity is heated not strictly from below is investigated in [2] and the case in which the perturbation of the equilibrium is due to the slow movement of the upper boundary of the region is investigated in [3]. In [2, 3] the perturbation has the structure of a first critical motion and thus the results of these papers coincide qualitatively. The perturbation of the temperature in the horizontal sections of the boundary, which creates a perturbation with a two-vortex structure corresponding to the second critical point R₂, is examined in this paper. A similar type of perturbation is characteristic for experiments in which the thermal conductivity properties of the fluid and the cavity walls are different. The nonlinear convection conditions are investigated numerically by the net-point method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 203–207, March–April, 1977.The author wishes to thank D. B. Lyubimova, V. I. Chernatynskii, and A. A, Nepomnyashchii for their helpful comments.     

Nonlinear capillary waves on a viscous fluid jet surface

Capillary instability of a fluid jet is one of the classical problems of hydrodynamics [1]. Studying it is of practical interest, particularly for the optimization of the ignition of a liquid propellant and the development of granulating apparatus in the chemical industry [2]. Until recently, the main attention has been paid to analyzing linear problems. Dispersion equations have been obtained for small perturbations of a jet surface with the viscosity of the external medium taken into account [3]. The construction of a theory of finite-amplitude waves on an ideal fluid jet surface was started in [4, 5]. Up to now this theory has achieved substantial results, as can be assessed by the successful numerical modeling of the dissociation of an inviscid fluid jet into drops [6] (see [7, 8] also). This paper is devoted to a discussion of the nonlinear development stage of viscous fluid jet instability under conditions allowing the influence of the surrounding medium and the gravity field to be neglected.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 179–182, March–April, 1977.The author is grateful to B. M. Konyukhov and G. D. Kuvatov for suggesting this problem and performing the experiment and to M. I. Rabinovich for useful discussions.     

Thermomechanical behavior of a rectangular viscoelastic prism exposed to repeated stretching and contraction

Vibrational heat production is an important problem in studying the efficiency of viscoelastic elements of structures experiencing cyclic loads. The design of heat regimes constitutes one of the fundamental problems in the construction of such types of vibrational-proof systems as laminated rods, plates, and shells [1] and fiberglass and rubber-metal products, in particular, shock absorbers [2, 3]. Calculation of the critical parameters beyond which a rapid growth in temperature occurs (the phenomenon of thermal explosion), which leads to partial or complete loss of the supporting power of the product as a result of softening of the material, is of particular interest. A variational method has been used [4] to calculate heat production in a twodimensional shock absorber. The boundary conditions are satisfied on the basis of the St. Venant principle. In the current work, the stress-strain state, self-heating temperature field, and thermal instability of a long rectangular prism being periodically loaded (plane deformation) are investigated.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1 pp, 149–156, January–February, 1976.     

Operating regimes of a chemical reactor with longitudinal and transverse mixing in the case of large Péclet numbers

A two-dimensional model of a chemical reactor with longitudinal and transverse mixing is investigated in the case of large Péclet numbers calculated from the effective thermal conductivity in the transverse direction. For this model the existence of at least one steady-state regime has been demonstrated [1], sufficient criteria of its uniqueness have been determined, an asymptotic expansion of the solution has been constructed in the case of small Péclet numbers, and the critical ignition and quenching parameters have been found. In this paper the other limiting case of the model, in which heat is propagated in the transverse direction much more slowly than it is transported by the flow along the reactor (large Péclet numbers), is analyzed in detail. An asymptotic expansion of the solution which closely coincides with the data of numerical calculations is constructed. The critical quenching and ignition conditions of the process are determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 120–127, January–February, 1987.     

Effect of viscosity on the autoignition of a reacting moving mixture

Ya. B. Zel'dovich has established [1] that in a continuous-flow reactor two ignition regimes are possible: forced ignition and autoignition.It is important to consider the special properties of the autoignition regime associated with the hydromechanics of laminar flow and heat transfer through the pipe wall. In [2, 3] it was shown that the effect of heat of friction on heat transfer in long pipes is qualitative in character. Moreover, according to Schlichting [4], in certain cases the temperature gradient for such flows due to the heat of friction may reach 10–30°, which is comparable with the preexplosion temperature rise in the stationary theory of thermal explosion [5]. In this connection it is clear that under certain conditions the heat of friction may considerably reduce the explosion limit.This paper is devoted to a study of the effect of heat of friction on the explosion limit of a reacting fluid in a long cylindrical pipe. The dynamic autoignition regime due to heat of friction is examined. In particular, it is established that, other things being equal, by increasing the pressure drop it is possible to obtain explosion of the reacting system.     

Investigation of the hypersonic viscous shock layer around blunt bodies in a nonuniform flow

Unseparated viscous gas flow past a body is numerically investigated within the framework of the theory of a thin viscous shock layer [13–15]. The equations of the hypersonic viscous shock layer with generalized Rankine-Hugoniot conditions at the shock wave are solved by a finite-difference method [16] over a broad interval of Reynolds numbers and values of the temperature factor and nonuniformity parameters. Calculation results characterizing the effect of free-stream nonuniformity on the velocity and temperature profiles across the shock layer, the friction and heat transfer coefficients and the shock wave standoff distance are presented. The unseparated flow conditions are investigated and the critical values of the nonuniformity parameter a_k [10] at which reverse-circulatory zones develop on the front of the body are obtained as a function of the Reynolds number. The calculations are compared with the asymptotic solutions [10, 12].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 154–159, May–June, 1987.     

Dynamic stability analysis and DQM for beams with variable cross-section

The present paper deals with the dynamic behaviour of a clamped beam subjected to a sub-tangential follower force at the free end. The aim of this work is to obtain the frequency–axial load relationship for a beam with a variable circular cross-section. In this way, one can identify both divergence critical loads – where the frequency goes to zero – and the flutter critical load – in correspondence with two frequencies coalescence. The numerical approach adopted for solving the partial differential equation of motion is the differential quadrature method (henceforth DQM). This method was proposed by Bellmann and Casti [Bellmann, R.E., Casti, J., 1971. Differential quadrature and long-term integration. J. Math. Anal. 34, 235–238] and has been employed recently in the solution of solid mechanics problems by Bert and Malik [Bert, C.W., Malik, M., 1996. Differential quadrature method in computational mechanics: a review. Appl. Mech. Rev., ASME, 49 (1), 1–28] and Chen et al. [Chen, W., Stritz, A.G., Bert, C.W., 1997. A new approach to the differential quadrature method for fourth-order equations. Int. J. Numer. Method Eng. 40, 1941–1956]. More precisely, a modified version of this method has been used, as proposed by De Rosa and Franciosi [De Rosa, M.A., Franciosi, C., 1998a. On natural boundary conditions and DQM. Mech. Res. Commun. 25 (3), 279–286; De Rosa, M.A., Franciosi, C., 1998b. Non classical boundary conditions and DQM. J. Sound Vibrat. 212(4), 743–748] to satisfy all the boundary conditions.Some frequencies–axial loads relationships are reported in order to show the influence of tapering on the critical loads.     

Influence of anisotropy of the permeability on the convection and heat transfer in a porous annular layer

Investigations into the convective transport of heat in porous materials are of interest for many applications in connection with the problem of increasing the efficiency of thermal insulation. In [1–5], convection in Isotropic porous media was considered. However, in many cases porous materials have an essential anisotropy of their permeability. Convective heat transfer has been inadequately studied for this case. In [6], the linearized equations were used to study the convection between infinite horizontal planes with a filling of an anisotropic material; the value of the critical Rayleigh number was found, and this agreed satisfactorily with experimental data. In the present paper, we investigate numerically convection between two infinite coaxial cylinders with an anisotropic porous filling, using the equations of convection in the Darcy—Boussinesq approximation [1–3]. The permeability tensor in the annular region is constructed from its principal values, which can be found experimentally. A method of calculation is developed and a parametric study made of the structure of the flow and of the local and averaged characteristics of the heat transfer, which are of interest for the design of thermal insulation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–64, January–February, 1980.     

Nonisothermal instability of high-velocity elastoplastic flows

An important feature of the high-velocity deformation of solids is the localization of deformation, one of the causes of which may be the nonisothermal instability of plastic flow [1–6]. In connection with the intensive development of high-velocity technology in the treatment of materials, the investigation of the criteria for nonisothermal stability of processes of plastic deformation is of fundamental interest, since in certain cases they determine the optimum technological regimes [5]. The critical values of deformation velocities, above which the effects of thermal instability becomes decisive in the process of deformation of solids, are estimated by semiempirical methods in [1]. The non-boundary-value problem of the criteria for nonisothermal instability is analyzed in [2] for the point of view of flow stability in the so-called coupled formulation. The latter means that the heat-conduction equation is added to the basic equations determining the dynamics of an elastoplastic medium. The problem is solved in [6] in an analogous formulation, but for flow averaged over the spatial coordinate. The solution of the boundary-value problem for one-dimensional flow in this formulation is given in the present paper.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 133–138, May–June, 1986.     

Convective stability of a weakly conducting liquid in an electric field

In an inhomogeneously heated weakly conductive liquid (electrical conductivity 10^–12–1 cm^–1) located in a constant electric field a volume charge is induced because of thermal inhomogeneity of electrical conductivity and dielectric permittivity. The ponderomotive forces which develop set the liquid into intense motion [1–6]. However, under certain conditions equilibrium proves possible, and in that case the question of its stability may be considered. A theoretical analysis of liquid equilibrium stability in a planar horizontal condenser was performed in [2, 4]. Critical problem parameters were found for the case where Archimedean forces are absent [2]. Charge perturbation relaxation was considered instantaneous. It was shown that instability is of an oscillatory character. In [4] only heating from above was considered. Basic results were obtained in the limiting case of disappearingly small thermal diffusivity in the liquid (infinitely high Prandtl numbers). In the present study a more general formulation will be used to examine convective stability of equilibrium of a vertical liquid layer heated from above or below and located in an electric field. For the case of a layer with free thermally insulated boundaries, an exact solution is obtained. Values of critical Rayleigh number and neutral oscillation frequency for heating from above and below are found Neutral curves are constructed. It is demonstrated that with heating from below instability of both the oscillatory and monotonic types is possible, while with heating from above the instability has an oscillatory character. Values are found for the dimensionless field parameter at which the form of instability changes for heating from below and at which instability becomes possible for heating from above.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 16–23, September–October, 1976.In conclusion, the author thanks E. M. Zhukhovitskii for this interest in the study and valuable advice.     

Optimization of body shape at small reynolds numbers

The problem of the optimization of the shape of a body in a stream of viscous liquid or gas was treated in [1–5]. The necessary conditions for a body to offer minimum resistance to the flow of a viscous gas past it were derived in [1], The necessary optimality conditions when the motion of the fluid is described by the approximate Stokes equations were derived in [2], The shape of a body of minimum resistance was found numerically in [3] in the Stokes approximation. The optimality conditions when the motion of the fluid is described by the Navier—Stokes equations were derived in [4, 5], and in [4] these conditions were extended to the case of a fluid whose motion is described in the boundary-layer approximation. The necessary optimality conditions when the motion of the fluid is described by the approximate Oseen equations were derived in [5] and an asymptotic analysis of the behavior of the optimum shape near the critical points was performed for arbitrary Reynolds numbers.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp, 87–93, January–February, 1978.     

Determining the form of a body resulting in minimum heat flux,with laminar flow in the boundary layer

The exact solution of the problem of determining the optimal body shape for which the total thermal flux will be minimal for high supersonic flow about the body involves both computational and theoretical difficulties. Therefore, at the present time wide use is made of the inverse method, based on comparing the thermal fluxes for bodies of various specified form [1, 2]. The results of such calculations cannot always replace the solution of the direct variational problem. Therefore it is advisable to consider the direct variational problem of determining the form of a body with minimal thermal flux by using the approximate Newton formula for finding the gasdynamic parameters at the edge of the boundary layer. This approach has been used in finding the form of the body of minimal drag in an ideal fluid [3–5] arid with account for friction [6], and also for determining the form of a thin two-dimensional profile with minimal thermal flux for given aerodynamic characteristics [7].     

Nonisothermal Couette flow of a non-Newtonian fluid under the action of a pressure gradient

Nonisothermal Couette flow has been studied in a number of papers [1–11] for various laws of the temperature dependence of viscosity. In [1] the viscosity of the medium was assumed constant; in [2–5] a hyperbolic law of variation of viscosity with temperature was used; in [6–8] the Reynolds relation was assumed; in [9] the investigation was performed for an arbitrary temperature dependence of viscosity. Flows of media with an exponential temperature dependence of viscosity are characterized by large temperature gradients in the flow. This permits the treatment of the temperature variation in the flow of the fluid as a hydrodynamic thermal explosion [8, 10, 11]. The conditions of the formulation of the problem of the articles mentioned were limited by the possibility of obtaining an analytic solution. In the present article we consider nonisothermal Couette flows of a non-Newtonian fluid under the action of a pressure gradient along the plates. The equations for this case do not have an analytic solution. Methods developed in [12–14] for the qualitative study of differential equations in three-dimensional phase spaces were used in the analysis. The calculations were performed by computer.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 26–30, May–June, 1981.     

Ignition of a gas in a boundary layer at a heated plate

The problem of the ignition of a moving homogeneous gaseous combustible mixture in a boundary layer along a heated flat semiinfinite plate is one of the main problems of the ignition of a combustible mixture in a flow (for example, [1]). The formulation of the problem includes the two-dimensional equations of motion and the equations of the transfer of heat and of the reacting substance, written taking a chemical reaction into consideration, as well as boundary conditions, and should lead to determination of the steady-state fields of the concentration and the temperature and, by the same token, of the position of the combustion zone. Different approximate numerical solutions of the problem were analyzed in [1–5]. One of the most important characteristics of the process is the length of the ignition, i.e., the distance from the edge of the plate to the point at which, thanks to the intrinsic chemical heat evolution in the gas, the heat flux from the plate to the gas becomes equal to zero. In the present work, for the case of large values of the activation energy of the chemical reaction and a sufficiently great temperature difference between the wall and the flow, an approximate expression is obtained for the length of the ignition.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 142–148, September–October, 1977.The authors thank V. M. Shevtsov for his aid in making the calculations.     

Transition flow of a fibrous suspension in a pipe as the flow of an anisotropic fluid

The transition flow is considered of a fibrous suspension in a pipe. The flow region consists of two subregions: at the center of the flow a plug formed by interwoven fibers and fluid moves as a rigid body; between the solid wall and the plug is a boundary layer in which the suspension is a mixture of the liquid phase and fibers separated from the plug [1–3]. In the boundary region the suspension is simulated as an anisotropic Ericksen—Leslie fluid [4, 5] which satisfies certain additional conditions. Equations are obtained for the velocity profile and drag coefficient of the pipe, which are both qualitatively and quantitatively in good agreement with the experimental results [6–8]. Within the framework of the model, a mechanism is found for reducing the drag in the flow of a fibrous suspension as compared to the drag of its liquid phase.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 91–98, September–October, 1985.     

Nonstatiqnary inner heating of glaciers

The problem of stationary inner heating of cold glaciers is considered in [1, 2]. It is shown that under certain conditions in the stationary state the heat released inside the glacier due to viscous dissipation is not carried off through the boundary surfaces, and that there is an analogy with the stationary theory of thermal explosions [3]. It was assumed that the temperature at the bed of the glacier could reach the melting point due to nonstationary heating and that the conditions would be produced for motion of the glacier. Here we examine the nonstationary response of a glacier to small deviations of the system parameters from critical values corresponding to the stationary state.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 169–172, September–October, 1975.     

The uniqueness of the solution of boundary-value problems of a thermal model of discharge

The paper is devoted to a nonlinear analysis of superheating [1, 2] instability of an electric discharge stabilized by electrodes [3] in the framework of a thermal model [4] where the stability of the discharge relative to the long-wave and short-wave perturbations is proved in a linear approximation. Similar boundary-value problems arise in the theories of chemically and biologically reacting mixtures [5–7], thermal breakdown of dielectrics [8], thermal explosion [9], in the investigation of nonlinear waves in semiconductors and superconductors [10, 11], and in the investigation of Couette flow with variable viscosity [12]. The uniqueness of the one-dimensional steady solutions of the thermal model of discharge and the stability relative to the small spatial perturbations, respectively, for the exponential and step dependence of the electrical conductivity on the temperature are proved in [3, 13]. The uniqueness of the solutions in the one-dimensional case for the same electrode temperature and arbitrary dependences of the electrical and thermal conductivity on the temperature is established in paper [14]. In the present paper, the existence and uniqueness of steady solutions of the thermal model of discharge in a three-dimensional formulation for arbitrary fairly smooth electrical and thermal conductivity functions of the temperature in the case of isothermal isopotential electrodes are proved analytically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 140–145, January–February, 1986.The author expresses his gratitude to A. G. Kulikovskii and A. A. Barmin for the formulation of the problem and their discussions.     

On nonisothermal electroconvection

The flow from the tip of a needle electrode is caused by the Coulomb force acting on the space charge [1–3]. This charge is formed because of the dependence of the conductivity on the temperature, nonuniformity of which is due to Joule heating [1] and the electric field intensity [2] or processes near the electrode [3–5]. The present paper considers the stability of a dielectric liquid between spherical electrodes in order to elucidate the possibility of a thermoelectrohydrodynainic flow due to Joule heating. In the presence of external heating, the possibility of such a flow has been demonstrated both experimentally and theoretically [6–8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 133–137, March–April, 1980.     

Theory of vortical helical ideal gas flows in laval nozzles

An asymptotic solution is found for the direct problem of the motion of an arbitrarily vortical helical ideal gas flow in a nozzle. The solution is constructed in the form of double series in powers of parameters characterizing the curvature of the nozzle wall at the critical section and the intensity of stream vorticity. The solution obtained is compared with available theoretical results of other authors. In particular, it is shown that it permits extension of the known Hall result for the untwisted flow in the transonic domain [1]. The behavior of the sonic line as a function of the vorticity distribution and the radius of curvature of the nozzle wall is analyzed. Spiral flows in nozzles have been investigated by analytic methods in [2–5] in a one-dimensional formulation and under the assumption of weak vorticity. Such flows have been studied by numerical methods in a quasi-one-dimensional approximation in [6, 7]. An analogous problem has recently been solved in an exact formulation by the relaxation method [8, 9]. A number of important nonuniform effects for practice have consequently been clarified and the boundedness of the analytical approach used in [2–7] is shown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–137, March–April, 1978.The authors are grateful to A. N. Kraiko for discussing the research and for valuable remarks.     

Boundary layer on a blunt body in a flow of dusty gas

The problem of interaction of gas-dust flows with solid surfaces arose in connection with the study of the motion of aircraft in a dusty atmosphere [1–2], the motion of a gas suspension in power generators, and in a number of other applications [3]. The presence of a disperse admixture may lead to a significant increase in the heat fluxes [4] and to erosion of the surface [5]. These phenomena are due to the joint influence of several factors — the change in the structure of the carrier-phase boundary layer due to the presence of the particles, collisions of the particles with the surface, roughness of the ablating surface, and so forth. This paper continues an investigation begun earlier [6–7] into the influence of particles on the structure of the dynamical and thermal two-phase boundary layer formed around a blunt body in a flow. The model of the dusty gas [8] has an incompressible carrier phase. The method of matched asymptotic expansions [9] is used to obtain the equations of the two-phase boundary layer. In the frame-work of the refined classification made by Stulov [6], it is shown that the form of the boundary layer equations is different in the presence and absence of inertial precipitation of the particles. The equations are solved numerically in the neighborhood of the stagnation point of the blunt body. The temperature and phase velocity distributions in the boundary layer, and also the friction coefficients and the heat transfer of the carrier phase are found for a wide range of the determining parameters. In the case of an admixture of low-inertia particles that are not precipitated on the body, it is shown that even when the mass concentration of the particles in the undisturbed flow is small their accumulation in the boundary layer can lead to a sharp increase in the thermal fluxes at the stagnation point.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 99–107, September–October, 1985.I thank V. P. Strulov for a discussion.