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.     

Oscillations of a vessel containing a viscous fluid

The oscillations of a rigid body having a cavity partially filled with an ideal fluid have been studied in numerous reports, for example, [1–6]. Certain analogous problems in the case of a viscous fluid for particular shapes of the cavity were considered in [6, 7]. The general equations of motion of a rigid body having a cavity partially filled with a viscous liquid were derived in [8]. These equations were obtained for a cavity of arbitrary form under the following assumptions: 1) the body and the liquid perform small oscillations (linear approximation applicable); 2) the Reynolds number is large (viscosity is small). In the case of an ideal liquid the equations of [8] become the previously known equations of [2–6]. In the present paper, on the basis of the equations of [8], we study the free and the forced oscillations of a body with a cavity (vessel) which is partially filled with a viscous liquid. For simplicity we consider translational oscillations of a body with a liquid, since even in this case the characteristic mechanical properties of the system resulting from the viscosity of the liquid and the presence of a free surface manifest themselves.The solutions are obtained for a cavity of arbitrary shape. We then consider some specific cavity shapes.     

Solution of a problem of flow in a porous medium with limiting gradient

For the law of flow in a porous medium with limiting gradient studied previously in [1], an exact solution is found for the problem formulated in [2] of the plane steady motion of an incompressible fluid in a channel with a rectangular step. Particular cases of the solution obtained are given; these represent the solutions of the problem of flow past a broken wall and of motion from a point source in a strip.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 76–78, January–February, 1985.     

Asymptotic theory of the flow around an obstacle by a sonic flow

The first investigation of the problem of the flow around an obstacle by a gas flow whose velocity is equal to the speed of sound at infinity was carried out in [1, 2], where it is shown in particular that the principal term of the appropriate asymptotic expansion is a self-similar solution of Tricomi's equation, to which the problem reduces in the first approximation upon a hodographic investigation. The requirement that the stream function be analytic as a function of the hodographic variables on the limiting characteristic was an important condition determining the selection of the self-similarity exponent n (xy^–n is an invariant of the self-similar solution). The analytic nature of the velocity field everywhere in the flow above the shock waves, which arise from necessity upon flow around an obstacle, follows from this condition. The latter was found in [3], where one of the branches of the solution obtained in [1] was used in the region behind the shock waves. The principal and subsequent terms of the asymptotic expansion describing a sonic flow far from an obstacle were discussed in [4], where the author restricted himself to Tricomi's equation. Each term of the series constructed in [4] contains an arbitrary coefficient (we will call it a shape parameter) which is not determined within the framework of a local investigation, and consideration of the problem of flow around a given obstacle as a whole is necessary in order to determine these shape parameters. It follows from the results of [4] that the problem of higher approximations to the solution of [1] coincides with the problem, of constructing a flow in the neighborhood of the center of a Laval nozzle with an analytic velocity distribution along the longitudinal axis (a Meyer-type flow). Along with the Meyer-type flow in the vicinity of the nozzle center, which corresponds to a self-similarity exponent n=2, two other types of flow are asymptotically possible with n=3 and 11, given in [5]. The appropriate solutions are written out in algebraic functions in [6]. The results of [5] show that the condition that the velocity vector be analytic on the limiting characteristic in the flow plane is broader than the condition that the stream function be analytic as a function of the hodographic variables, which is employed in [1, 2, 4]. Therefore, the necessity has arisen of reconsidering the problem of higher approximations for the obstacle solution of F. I. Frankl'. It has proved possible for the region in front of the shock waves to use a series which is more general than in [4], which implies the inclusion of an additional set of shape parameters. The solution is given in the hodograph plane in the form of the sum of two terms; the series discussed in [4] corresponds to the first one, and the series generated by the self-similar solution with n=3 or with n=11 corresponds to the second one.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 99–107, May–June, 1979.The authors thank S. V. Fal'kovich for a useful discussion.     

Motion of inhomogeneities of a developed fluidized bed at small reynolds numbers

The instability of a fluidized system in which the particles are uniformly distributed in space [1–3] leads to the development of local inhomogeneities in the internal structure, these taking the form of more or less stable formations of packets of particles [4]. In accordance with the existing ideas based on experimental data [5–8, 13], the particle concentration within a packet may vary in a wide range from very small values (10^–2–10^–3 [8]) for bubbles to the concentration of the unfluidized bed for bunches of particles in a nearly closely packed state. The paper considers the steady disturbed motion of the fluid and solid phases near an ascending or descending packet of particles in a developed fluidized bed. It is assumed that the motion of the solid phase corresponds to a creeping flow of viscous fluid, and the viscosity of the fluidizing agent is taken into account only in the terms that describe the interphase interaction. The velocity fields and pressure distributions of the phases inside and outside a packet are determined. If the particle concentration within a packet tends to zero, the solution describes the slow motion of a bubble in a fluidized bed. The results of the paper are compared with results obtained earlier for the model of ideal fluids [9] and Batchelor's model [10], in which the fluidized bed is treated in a simplified form as a viscous quasihomogeneous continuum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 57–65, July–August, 1984.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

Invariant solutions of the equations of a multicomponent charged-particle beam

The concept of the invariant-group solution (H-solution) was introduced and a general method for obtaining it was developed in [1–3]. The group properties of the equations of a monoenergetic charged-particle beam with the same value and sign of the specific charge, assuming univalency of the velocity vector V, were studied in [4–6], where all essentially different H-solutions were also constructed. Below, the results of [4–6] are extended to the case of a beam in the presence of a fixed background of density ₀ (§1), and also to the case of multivelocity (V is an s-valued function) and multicomponent beams (i.e., beams formed by particles of several kinds) (§2). A number of analytic solutions that describe some nonstationary processes in devices with plane, cylindrical, and spherical geometry —among them a continuous periodic solution for a plane diode with a period determined by the background density -are obtained in §1. A transformation that contains arbitrary functions of time and preserves Vlasov's equations is given (§2). The equations studied can be treated as the equations of a rarefied plasma in the magnetohydrodynamic approximation, when the pressure gradients are negligible as compared with forces of electromagnetic origin.     

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.     

Motion of inviscid fluid tongues and bubbles in a Hele-Shaw cell occupied by a viscoelastic fluid

Problems of steady-state motion of an inviscid fluid tongue or bubble in a Hele-Shaw cell occupied by a viscoelastic fluid are considered. The problems generalize the well-known Saffman-Taylor problems [1, 2] in the approximation in which the pressure jump on the interface between the two fluids depends on the normal velocity of the interface [4]. In the zero approximation the results obtained by means of an iteration technique are in good agreement with the analytic solutions obtained by Saffman and Taylor [1, 2].__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 117–123. Original Russian Text Copyright © 2005 by Entov, Kolganov, and Kolganova.     

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.     

Self-similar problems of propagation of an extended hydrofracture in a permeable medium

The propagation of an extended hydrofracture in a permeable elastic medium under the influence of an injected viscous fluid is considered within the framework of the model proposed in [1, 2]. It is assumed that the motion of the fluid in the fracture is turbulent. The flow of the fluid in the porous medium is described by the filtration equation. In the quasisteady approximation and for locally one-dimensional leakage [3] new self-similarity solutions of the problem of the hydraulic fracture of a permeable reservoir with an exponential self-similar variable are obtained for plane and axial symmetry. The solution of this two-dimensional evolution problem is reduced to the integration of a one-dimensional integral equation. The asymptotic behavior of the solution near the well and the tip of the fracture is analyzed. The difficulties of using the quasisteady approximation for solving problems of the hydraulic fracture of permeable reservoirs are discussed. Other similarity solutions of the problem of the propagation of plane hydrofractures in the locally one-dimensional leakage approximation were considered in [3, 4] and for leakage constant along the surface of the fracture in [5–7].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 91–101, March–April, 1992.     

Some problems in the theory of plane jets of conducting fluid

Using the boundary-layer equations as a basis, the author considers the propagation of plane jets of conducting fluid in a transverse magnetic field (noninductive approximation).The propagation of plane jets of conducting fluid is considered in several studies [1–12]. In the first few studies jet flow in a nonuniform magnetic field is considered; here the field strength distribution along the jet axis was chosen in order to obtain self-similar solutions. The solution to such a problem given a constant conductivity of the medium is given in [1–3] for a free jet and in [4] for a semibounded jet; reference [5] contains a solution to the problem of a free jet allowing for the dependence of conductivity on temperature. References [6–8] attempt an exact solution to the problem of jet propagation in any magnetic field. An approximate solution to problems of this type can be obtained by using the integral method. References [9–10] contain the solution obtained by this method for a free jet propagating in a uniform magnetic field.The last study [10] also gives a comparison of the exact solution obtained in [3] with the solution obtained by the integral method using as an example the propagation of a jet in a nonuniform magnetic field. It is shown that for scale values of the jet velocity and thickness the integral method yields almost-exact values. In this study [10], the propagation of a free jet is considered allowing for conduction anisotropy. The solution to the problem of a free jet within the asymptotic boundary layer is obtained in [1] by applying the expansion method to the small magnetic-interaction parameter. With this method, the problem of a turbulent jet is considered in terms of the Prandtl scheme. The Boussinesq formula for the turbulent-viscosity coefficient is used in [12].This study considers the dynamic and thermal problems involved with a laminar free and semibounded jet within the asymptotic boundary layer, propagating in a magnetic field with any distribution. A system of ordinary differential equations and the integral condition are obtained from the initial partial differential equations. The solution of the derived equations is illustrated by the example of jet propagation in a uniform magnetic field. A similar solution is obtained for a turbulent free jet with the turbulent-exchange coefficient defined by the Prandtl scheme.     

Equations of laminar boundary layer in a two-phase medium

The motion of a two-phase medium in which the carrier component has low viscosity is considered. The equations obtained in [1], to which the viscous stress tensor in the fluid is added, are used. The boundary layer method [2] makes it possible to obtain asymptotic equations for the wall region. These equations have different forms depending on the characteristic values of the dimensionless determining parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 51–60, January–February, 1979.I thank A. N. Kraiko for discussing the work.     

Choice of a self-similar solution in boundary-layer theory

If the speed of the outer flow at the edge of the boundary layer does not depend on the time and is specified in the form of a power-law function of the longitudinal coordinate, then a self-similar solution of the boundary-layer equations can be found by integrating a third-order ordinary differential equation (see [1–3]). When the exponent of the power in the outerflow velocity distribution is negative, a self-similar solution satisfying the equations and the usually posed boundary conditions is not uniquely determinable [4], A similar result was obtained in [5] for flows of a conducting fluid in a magnetic field. In the present paper we study the behavior of non-self-similar perturbations of a self-similar solution, enabling us to provide a basis for the choice of a self-similar solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 42–46, July–August, 1974.     

Waves from an underground explosion

The problem of the propagation of a spherical detonation wave in water-saturated soil was solved in [1, 2] by using a model of a liquid porous multicomponent medium with bulk viscosity. Experiments show that soils which are not water saturated are solid porous multicomponent media having a viscosity, nonlinear bulk compression limit diagrams, and irreversible deformations. Taking account of these properties, and using the model in [2], we have solved the problem of the propagation of a spherical detonation wave from an underground explosion. The solution was obtained by computer, using the finite difference method [3]. The basic wave parameters were determined at various distances from the site of the explosion. The values obtained are in good agreement with experiment. Models of soils as viscous media which take account of the dependence of deformations on the rate of loading were proposed in [4–7] also. In [8] a model was proposed corresponding to a liquid multicomponent medium with a variable viscosity.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 34–41, May–June, 1984.     

Dynamics of a cavitational cavity in a non-newtonian fluid

The dynamics of a spherical cavity in a non-Newtonian fluid, described by the Reiner-Rivlin rheological equation [1], is investigated. The equation of radial cavity motion is obtained, where the gas in the cavity is subject to a polytropic law and surface tension is taken into account. The equation of cavity motion is solved numerically for a number of values of the transverse viscosity coefficient. The influence of the transverse viscosity on the collapse process of vapor and gas-filled cavities is shown. Numerical computations are also carried out for the rate of energy dissipation and the pressure distribution in the fluid.Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 170–173, July–August, 1973.The authors are grateful to A. T. Listrove for attention to the research.     

Computation of the pressure on the surface of a fuselage with a wing in an ideal fluid

The flows around complex three-dimensional bodies by an ideal fluid are computed by methods [1–3] using approximation of the surface by a set of plane elements. A layer of surface singularities, whose intensity is found by solving a system of linear algebraic equations of very high order, is distributed continuously over each element. Evaluation of the system coefficients and its solution require significant machine time expenditures on powerful electronic computers. If in the method of [2] the total system of equations is separated successfully into several subsystems by simplifications and an approximate solution of the problem is obtained more rapidly than by the method in [1], then the same author practically used the method in [1] to design specific fuselages in [3]. A method [4] developed for a fuselage is expanded in this paper to design a wing-fuselage combination. This method turns out to be less tedious, without being inferior in accuracy, by being different from the method in [1] in the means of solving the fundamental integral equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–115, May–June, 1977.     

Strong mass injection at the surface of a blunt body in a supersonic flow

The case of supersonic flow over a blunt body when another gas is injected through the surface of the body in accordance with a given law is theoretically investigated. If molecular transport processes are neglected, the flow between the shock wave and the surface of the body should be regarded as two-layer, that is, as consisting of the flow in the shock layer between the shock wave and the contact surface and the flow in the layer of injected gas. A numerical solution of the problem is obtained near the front of the body and its accuracy is estimated. Approximate analytic solutions are obtained in the injected-gas layer: a constant-density solution and a solution of the boundary-layer type in the local similarity approximation. Near the flow axis the numerical and analytic solutions are fairly close, but at a distance from the axis the assumptions made reduce the accuracy of the approximate solutions. The flow in question can serve as a gas-dynamic model of a series of problems describing the radiant heating of blunt bodies in a hypersonic flow. In the presence of intense radiative heat transfer, vaporization is so great that the thickness of the vapor layer is comparable with the thickness of the shock layer. Moreover, the thermal shielding of various kinds of obstacles in channels through which a radiating plasma flows can be organized by means of the forced injection of a strong absorber. The formulation of a similar problem was reported in [1], but the results of the solution were not given. A two-layer model of the flow of an ideal gas over a blunt body was used in [2, 3] for the analysis of radiative heat transfer. In [2] the neighborhood of the stagnation point is considered. In [3] preliminary results relating to two-layer flow over blunt cones are presented. The solution is obtained by Maslen's approximate method.Moscow. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 89–97, March–April, 1972.     

Surface waves and stability of tangential velocity discontinuity on a solid-fluid boundary

Solutions of the Rayleigh-wave type on the boundary of an elastic half-space and a moving layer of ideal fluid are obtained. The limiting cases of zero flow velocity and a tangential velocity discontinuity in the fluid were investigated in [1–3]. In [4] the order of magnitude of the critical flow velocity was estimated. An increase in the velocity scales used in engineering and experimental practice (see [5], for instance) has aroused interest in a more thorough analysis of the effect.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 43–46, May–June, 1981.     

Dynamics of a flat cavity in a low-viscosity liquid

The problem of the motion of a cavity in a plane-parallel flow of an ideal liquid, taking account of surface tension, was first discussed in [1], in which an exact equation was obtained describing the equilibrium form of the cavity. In [2] an analysis was made of this equation, and, in a particular case, the existence of an analytical solution was demonstrated. Articles [3, 4] give the results of numerical solutions. In the present article, the cavity is defined by an infinite set of generalized coordinates, and Lagrange equations determining the dynamics of the cavity are given in explicit form. The problem discussed in [1–4] is reduced to the problem of seeking a minimum of a function of an infinite number of variables. The explicit form of this function is found. In distinction from [1–4], on the basis of the Lagrauge equations, a study is also made of the unsteady-state motion of the cavity. The dynamic equations are generalized for the case of a cavity moving in a heavy viscous liquid with surface tension at large Reynolds numbers. Under these circumstances, the steady-state motion of the cavity is determined from an infinite system of algebraic equations written in explicit form. An exact solution of the dynamic equations is obtained for an elliptical cavity in the case of an ideal liquid. An approximation of the cavity by an ellipse is used to find the approximate analytical dependence of the Weber number on the deformation, and a comparison is made with numerical calculations [3, 4]. The problem of the motion of an elliptical cavity is considered in a manner analogous to the problem of an ellipsoidal cavity for an axisymmetric flow [5, 6]. In distinction from [6], the equilibrium form of a flat cavity in a heavy viscous liquid becomes unstable if the ratio of the axes of the cavity is greater than 2.06.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 15–23, September–October, 1973.The author thanks G. Yu. Stepanov for his useful observations.