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.     

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.     

Small vibrations of a sphere in a spherical volume of viscous fluid

Problems of the vibration of bodies in confined viscous fluids have been solved to determine the added masses and damping coefficients of rods [1–3] and floats [4–5]. The solutions of these problems, based on the use of simplifications of the boundary-layer method [4–6], are obtained analytically in general form and are in good agreement with the experimental data. However, in each specific case the possibility of using such solutions for given values of the fluid viscosity and vibration frequency must be justified either experimentally [2, 4, 5] or theoretically as, for example, in [1], where an analytic solution was obtained for concentric cylinders. The present paper offers a general solution of the problem of the small vibrations of a sphere in a spherical volume of fluid valid over a broad range of variation of the dimensionless kinematic viscosity. The limiting cases of this solution for both high and low viscosity are considered. The asymptotic expressions obtained are compared with calculations based on the analytic solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 29–34, March–April, 1986.     

Effect of plastic properties of fluids on phase permeabilities

The linearity of Darcy's law is known to be disturbed at both high and low flow velocities [1–3]. In the first case, this is associated with the increase in the inertial component of the hydraulic losses in the presence of large pressure gradients. The effect was theoretically investigated, for example, in [4]. In the second case, the nonlinearity is associated with the interphase interaction of the fluids and the skeleton of the porous material on the contact surface [5]. Here, within the context of the percolation approach [6, 7], the behavior of the phase permeabilities is analyzed for low flow velocities, when on the microlevel (flow in an individual pore channel) the fluids display plasticity properties [8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 110–115, March–April, 1991.     

General stability criterion for laminar tube flow

Initiation of turbelence is associated with disturbances of finite intensity [1]. Some attempts, as in [2], have been made to treat this region analytically. The nonuniformity of the local stability * of laminar fluid flow over the tube cross section has been established experimentally [3,4, 1]. On this finding is based the interpretation of a number of turbulent transition phenomena, including a characteristic singularity of the relationship between the resistance coefficient of rough tubes and the Reynolds number under transient conditions [5].Expression (1.1) introduced as a measure of stability, and the criterion q_* yield satisfactory quantitative results.     

Axially symmetric vortical flow of a nonviscous liquid in the semiinfinite gap between two coaxial circular cylinders

In the framework of the Hromek-Lamb equations we investigate the axially symmetric vortical flow of a nonviscous incompressible liquid in both semiinfinite and infinite gaps between two coaxial circular cylinders. The investigation is carried out for two circulation and flow functions and two different Bernoulli constants which are chosen in the form of a third-order polynomial in the flow function. This makes it possible to determine the effect of the azimuthal velocity component on the flow in an axial plane with radial and axial components of the velocity. It is shown that under certain circumstances wave oscillations in the flow are possible, in agreement with the results of [1–3] which investigated the flow in an infinite tube [1], in a semiinfinite tube with simpler circulation functions and Bernoulli constants [2], and in the two-dimensional case [3]. We determine the dependence of the formation of wave perturbations on the third term of the Bernoulli constant and on the azimuthal velocity component. The results of this work agree with investigations by other authors [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 38–45, September–October, 1977.The author thanks Yu. P. Gupalo and Yu. S. Ryazantsev for suggesting this problem and for their interest in the work. Thanks are also due to G. Yu. Stepanov for discussions and valuable comments.     

Motion of a spherical solid particle in a nonuniform flow of a viscous incompressible liquid

The effect of a particle on the basic flow is studied, and the equations of motion of the particle are formulated. The problem is solved in the Stokes approximation with an accuracy up to the cube of the ratio of the radius of the sphere to the distance from the center of the sphere to peculiarities in the basic flow. An analogous problem concerning the motion of a sphere in a nonuniform flow of an ideal liquid has been discussed in [1]. We note that the solution is known in the case of flow around two spheres by a uniform flow of a viscous incompressible liquid [2], and we also note the papers [3, 4] on the motion of a small particle in a cylindrical tube.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 71–74, July–August, 1976.     

Stationary flow past the lower surface of a piecewise planar delta wing with supersonic leading edges

The considered wing has any finite number of inflections in its plane with lines of inflection intersecting at the point of inflection of the leading edge. In the present paper, this generalizes the author's earlier work [1] on flow past the undersurface of a flat wing at unite angle of attack with finite angle of slip and supersonic leading edges. In [1], calculations were not given. The special case of flow without slip in the same situation was considered later in [2], However, this paper contains errors, indicated at the end of the present paper. The calculations given in [2] are not correct. In the quoted papers, the gas flow is assumed to be a perturbation of a homogeneous flow behind a plane oblique shock wave. Such flows are treated systematically in [3]. Here and in [1], we use and generalize the representation of the linearized conservation laws across the shock front as the conditions of a boundary-value problem for an analytic function of a complex variable as obtained in [4, 5]. Calculations are given of the pressure distribution over the span for a number of different flow regimes and the pressure coefficients in the middle of the wing are compared with a numerical solution presented partly in [6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 80–90, September–October, 1979.I am very grateful to V. I. Lapygin for making available a large number of variants of his numerical solution, and to L. E. Pekurovskii for assistance in the calculations.     

Flow of a viscous fluid in a closed cavity in the presence of slip effects

Slip at the wall is observed in the flow of non-Newtonian fluids [1–4] and rarefied gases [5]. The most complete information on the phenomenon is obtained in capillary viscosimetry. For small radii of the capillaries and in porous media the slip effect is manifested even for Newtonian fluids (water, kerosene, for example) [6]. Experiments [2, 4] show that the influence of the entrance section can be ignored if the length of the capillary exceeds its radius by about 100 times. For the measurement of the rheological characteristics of high-viscosity fluids the use of long capillaries is difficult, and it is necessary to calculate the two-dimensional flow at the entrance section with allowance for slip. The need for such calculations also arises, for example, when one is choosing the optimal parameters of the screw devices employed in the processing of polymers [7]. Two-dimensional flows of a viscous incompressible fluid are frequently calculated with the flow function and vorticity =– used as variables [8–14]. The expressions for the vorticity on the boundary are usually obtained from the viscous no-slip condition [8, 9]. In the present paper, expressions are obtained for the vorticity on a wall in the presence of slip. The obtained expressions are used to solve a test problem on the flow of a viscous incompressible fluid in a cavity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 10–16, January–February, 1980.     

Existence and uniqueness of a solution to a problem of steady composite waves on the surface of a liquid flowing over a wavy bottom

Composite waves on the surface of the stationary flow of a heavy ideal incompressible liquid are steady forced waves of finite amplitude which do not disappear when the pressure on the free surface becomes constant but rather are transformed into free nonlinear waves [1]. It will be shown that such waves correspond to the case of nonlinear resonance, and mathematically to the bifurcation of the solution of the fundamental integral equation describing these waves. In [2], a study is made of the problem of composite waves in a flow of finite depth generated by a variable pressure with periodic distribution along the surface of the flow. In [3], such waves are considered for a flow with a wavy bottom. In this case, composite waves are defined as steady forced waves of finite amplitude that, when the pressure becomes constant and the bottom is straightened, do not disappear but are transformed into free nonlinear waves over a flat horizontal bottom. However, an existence and uniqueness theorem was not proved for this case. The aim of the present paper is to fill this gap and investigate the conditions under which such waves can arise.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 88–98, July–August, 1980.     

Drag associated with transitional flow regime of a fiber suspension in a tube

A semlempirical model is constructed of the flow of a fiber suspension of low and medium concentration in regimes that are usually called mixed and undeveloped turbulent regimes [1–4]. It is shown that although the flow of fiber suspensions in these regimes has features similar to those of the turbulent flow of a Newtonian fluid, for example, a logarithmic velocity profile, the characteristic features of the flow in both regimes can be better explained, not by turbulence of the flow, but by orientation of the fibers in it and by plastic flow of the fiber continuum. For this reason, to distinguish the mixed and undeveloped turbulent regimes from a truly turbulent regime it is proposed here to describe them by a general name — transitional flow. The obtained expressions agree qualitatively and quantitatively with the experimental results of Lee and Duffy [2], Sanders and Meyer [3], and Mih and Parker [4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 40–45, May–June, 1984.I thank V. N. Nikolaevskii and A. N. Golubyatnikov for interest in the work and helpful comments.     

Flow of a conductive gas in a circular tube in a strong axiosymmetric magnetic field

The flow of a conductive gas along a channel in an external axiosymmetric magnetic field with a finite value of the magnetogasodynamic parameter N is examined. Numerical flow calculations are performed for a circular tube in such a field. Gas dynamic parameter fields, total pressure losses, and electric current intensities with the presence of transsonic zones and highly compressed regions are determined. Through comparison of the results obtained with linear theory data, the range of applicability of the latter is determined. Of the works dedicated to study of flow in external magnetic fields with N1, we should take note of [1], in which the process of entry of the gas into a transverse magnetic field was examined; [2], which studied one-dimensional transient motion with shock waves; and [3], where mixed flow in a Laval nozzle with an axiosymmetric homogeneous magnetic field was studied. Flow in a circular tube was examined in [4]; but the analysis performed by the characteristic method permitted calculation of only the initial supersonic flow zone. Motion in circular tubes in the presence of an axiosymmetric, magnetic field was studied in the linear formulation in [4, 5].Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–155, September–October, 1972.     

Effect of surface tension on the free outflow of a wetting fluid from a horizontal tube

The flow structure for which the Kármán hypothesis is valid is explored experimentally and theoretically. It is established that there exists not a point but a region of finite size, adjacent to the upper generator of the tube and including the moving line of contact between the phase interface and the wall, in which the fluid is stationary relative to the line of contact, i.e., in which the no-slip condition is not satisfied. The dimensions of the region depend on the surface tension. The action of this stagnant zone on the flow is fully explained by the effect of the surface tension in experiments [1]. It is established that depending on the ratio of the tube diameter to the dimension of the stagnant zone two flow regimes are possible: in sufficiently wide tubes, an inertial regime for which Kármán's hypothesis holds and the no-slip condition is not satisfied, and, in sufficiently narrow tubes, a creep regime in which the no-slip condition continues to apply. The values of the determining dimensionless parameter corresponding to the change of regime and the cessation of flow are calculated. They are similar to the experimental values.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 85–93, March–April, 1988.     

Electrization of dielectric liquids flowing in tubes

The methods of the mechanics of continuous media [1] are used to consider the problem of electrization of dielectric liquids flowing in tubes [2–6]. According to modern ideas [2–6], there is always dissolved in such liquids a slight admixture of an electrolyte, whose molecules in such a dilute solution dissociate to a certain extent into positively and negatively charged ions. On the walls, oxidizing and reducing reactions take place, as a result of which the negative and positive ions, respectively, give up to the wall surplus electrons or take missing electrons from it. Thus, a positive (respectively, negative) total electric charge is induced in the liquid by the flow. We consider in this paper the electrization of a dielectric liquid in laminar flow in a circular cylindrical tube. We find the distribution of the electric charge in the liquid, the maximal electric current, and the dependence of the length over which the distribution of the electric charge in the tube is established on the tube radius, the Debye radius of the liquid, and the Péclet diffusion number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 41–47, November–December, 1979.We thank V. V. Gogosov for helpful comments made in a discussion of thwe work.     

Stability of frontal displacement of fluids in a porous medium in the presence of interphase mass transfer and phase transitions

To preserve the stability of the front relative to small perturbations when one fluid is displaced by another the pressure gradient must decrease on crossing the front in the direction of displacement. Initially, this criterion was established for the piston displacement of fluids [1, 2], and later in the case of two-phase flow of immiscible fluids in porous media for the displacement front corresponding to the saturation jump in the Buckley—Leverett problem [3, 4]. Below it is shown that the same stability criterion remains valid for flows in porous media accompanied by interphase mass transfer and phase transitions [5, 6]. Processes of these kinds are encountered in displacing oil from beds using active physicochemical or thermal methods [7] and usually reduce to pumping into the bed a slug (finite quantity) of reagent after which a displacing agent (water or gas) is forced in. The slug volume may be fairly small, especially when expensive reagents are employed, and, accordingly, in these cases the question of the stability of displacement is one of primary importance. These active processes are characterized by the formation in the displacement zone of multiwave structures which, in the large-scale approximation (i.e., with capillary, diffusion and nonequilibrium effects neglected), correspond to discontinuous distributions of the phase saturations and component concentrations [5–10]. It is shown that the stability condition for a plane front, corresponding to a certain jump, does not depend on the type of jump [11, 12] and for a constant total flow is determined, as in simpler cases, by the relation between the total phase mobilities at the jump. An increase in total flow in the direction of displacement is destabilizing, while a decrease has a stabilizing influence on the stability of the front. Other trends in the investigation of the stability of flows in porous media are reviewed in [13].Translated fron Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 98–103, March–April, 1986.     

Forced oscillations of pressure in tubes of active material

A study is made in the quasione-dimensional inertialess approximation of the axisymmetric flow of a Newtonian fluid in a tube of finite length made of a nonlinear active material with the capability of reducing deformations in response to an increase in tensile stresses [1, 2]. A study is made of the influence of the frequency and amplitude of forced oscillations of pressure at the entrance of the tube on its flow rate characteristics and on the behavior of the tube, depending on its length and certain rheological parameters. The first attempts at a study within the framework of this model of flow for unsteady conditions at the ends of the tube and in the ambient medium are described in [3, 4]. A general solution of this problem for external periodic disturbances of low amplitude is constructed in [5]. The present study gives an analysis of certain results of the numerical solution of an analogous problem for a wide range of variations in the frequency and amplitude of the pressure oscillations at the entrance to the tube.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 88–90, March–April, 1985.     

Motion of a fluid between rotating cones

A study is made of the steady axisymmetric flow of a viscous fluid between two cones rotating in opposite ways round a common axis. It is shown that as in the case of the flow of fluid swirled by plane disks rotating at different speeds [1], there can be two regimes of motion in the system: a Batchelor regime with quasirigid rotation of the fluid outside the boundary layers [2] and a Stewartson regime in which the azimuthal flow is concentrated only in the boundary layers [3]. In the Stewartson regime, a boundary layer analogous to that in the single disk problem (see, for example, [4–6]) forms in the region of each cone far from the apex. For the flows outside the boundary layers, simple expressions are found which make it possible to obtain a conception of the circulation of the fluid as a whole. With minor alterations, the results can be applied to the case of the rotation of other curved surfaces.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 58–64, March–April, 1985.The author thanks A. M. Obukhov for suggesting the subject and for his interest in the work, and A. V. Danilov and S. V. Nesterov for useful discussions.     

Flow processes with change of regime in fractured reservoirs

Experimental and industrial observations indicate a strong nonlinear dependence of the parameters of the flow processes in a fractured reservoir on its state of stress. Two problems with change of boundary condition at the well — pressure recovery and transition from constant flow to fixed bottom pressure — are analyzed for such a reservoir. The latter problem may be formulated, for example, so as not to permit closure of the fractures in the bottom zone. For comparison, the cases of linear [1] and nonlinear [2] fractured porous media and a fractured medium [3] are considered, and solutions are obtained in a unified manner using the integral method described in [1]. Nonlinear elastic flow regimes were previously considered in [3–6], where the pressure recovery process was investigated in the linearized formulation. Problems involving a change of well operating regime were examined for a porous reservoir in [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 67–73, May–June, 1991.     

A method for calculating rheological and morphological properties of constant-volume polymer blend models in inhomogeneous shear fields

We show how to formulate two-point boundary-value problems in order to compute fully-developed laminar channel and tube flow profiles for viscoelastic fluid models. The formulation is applied to Couette and pressure-driven flows separately, or a combination of both. The application of this methodology is illustrated analytically for the Upper-Convected Maxwell Model, and it is applied computationally for the Phan-Thien/Tanner and Giesekus Models. Numerical solutions exist for the last two models [J.Y. Yoo, H.C. Choi, On the steady simple shear flows of the one-mode Giesekus fluid, Rheol. Acta 28 (1989) 13–24; P.J. Oliveira, F.T. Pinho, Analytical solution for fully developed channel and pipe flow of Phan-Thien–Tanner fluids, J. Fluid Mech. 387 (1999) 271–280; M.A. Alves, F.T. Pinho, P.J. Oliveira, Study of steady pipe and channel flows of a single-mode Phan-Thien–Tanner fluid, J. Non-Newtonian Fluid Mech. 101 (2001) 55–76], allowing verification of the computational technique. Subsequently, the computational algorithm is applied to the constant-volume polymer blend models of Maffettone and Minale [P.L. Maffettone, M. Minale, Equation of change for ellipsoidal drops in viscous flow, J. Non-Newtonian Fluid Mech. 84 (1999) 105–106 (Erratum), J. Non-Newtonian Fluid Mech. 78 (1998) 227–241] and Dressler and Edwards [M. Dressler, B.J. Edwards, The influence of matrix viscoelasticity on the rheology of polymer blends, Rheol. Acta 43 (2004) 257–282; M. Dressler, B.J. Edwards, Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution, J. Non-Newtonian Fluid Mech. 120 (2004) 189–205]. Rheological and morphological properties of the model blends are thus obtained as functions of the spatial position within the channel, applied pressure drop, and shear rate at the wall.     

Flow of a viscous heat-conducting gas and binary mixtures in a sink

The class of exact solutions of the one-dimensional Navier-Stokes equations corresponding to gas flows from a spherical source or sink has been investigated analytically and numerically on a number of occasions (see, for example, [1, 2]). Here, the solution for a sink is considered in the presence of heat transfer from the ambient medium. Apart from seeking the solution itself, the object of the investigation was to establish the conditions of transi tion from viscous to inviscid flow in the sink as the Reynolds number tends to infinity. As shown in [3], for zero heat flux at an infinitely remote point there is no such transition for flow in a sink. The sink flow characteristics of a binary gas mixture are investigated in detail. In the transonic flow region an asymptotic solution is obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 56–62, January–February, 1989.