Flow of a viscous fluid in an annular gap with intensive blowing from the walls

Self-similar solutions are obtained in [1, 2] to the Navier-Stokes equations in gaps with completely porous boundaries and with Reynolds number tending to infinity. Approximate asymptotic solutions are also known for the Navier-Stokes equations for plane and annular gaps in the neighborhood of the line of spreading of the flow [3, 4]. A number of authors [5–8] have discovered and studied the effect of increase in the stability of a laminar flow regime in channels of the type considered and a significant increase in the Reynolds number of the transition from the laminar regime to the turbulent in comparison with the flow in a pipe with impermeable walls. In the present study a numerical solution is given to the system of Navier-Stokes equations for plane and annular gaps with a single porous boundary in the neighborhood of the line of spreading of the flow on a section in which the values of the local Reynolds number definitely do not exceed the critical values [5–8]. Generalized dependences are obtained for the coefficients of friction and heat transfer on the impermeable boundary. A comparison is made between the solutions so obtained and the exact solutions to the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 21–24, January–February, 1987.     

Initial section in a square duct rotating about the transverse axis

Laminar flow in a duct uniformly rotating about the transverse axis has been the subject of a number of studies. These are reviewed, for example, in [1]. However, in all these studies only the case of stablized flow was subjected to theoretical investigation. In this paper the development of laminar flow in a prismatic duct uniformly rotating about an axis perpendicular to one of the walls is analyzed on the basis of the results of a numerical integration of the parabolized system of Navier-Stokes equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 41–46, September–October, 1985.     

Calculation of three-dimensional flow in stationary and rotating channels with permeable walls

The results of numerical integration of the complete Navier-Stokes equations are presented for steady-state laminar incompressible viscous flow in a rotating channel of square cross section with a blind end and two permeable walls perpendicular to the axis of rotation. The formation of reverse flow zones is noted.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 179–181, September–October, 1989.The author is grateful to E. M. Smirnov for discussing his results.     

Example of parametric resonance in a rotating flow

Parametric resonance is one of the common types of instability of mechanical systems [1]. A standard example of the equations describing parametric oscillations is the Mathieu equation and its generalizations. In hydrodynamics these oscillations have been closely studied in connection with the problem of the vertical oscillations of a vessel containing an incompressible fluid in a uniform gravity field [1–5]. In this paper a new example of a flow whose stability problem reduces to the Mathieu equation is given. This is a flow of special type in a rotating cylindrical channel. The direction of the angular velocity is perpendicular to the channel axis, and its magnitude varies periodically with time. Flows with this geometry are of potential interest in technical applications [6, 7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 175–177, March–April, 1987.     

Analysis of flow in an annular duct with large injection at the walls

Viscous heat-conducting compressible fluid flow in an annular duct formed by two coaxial cylinders with large injection at the walls is investigated. An asymptotic solution exhibiting the influence of the axial symmetry of the duct is obtained in the vicinity of the y axis and is compared with the results of exact numerical calculations. Asymptotic solutions of the Navier-Stokes equations have been obtained earlier for flows in a plane channel with various rates of wall injection in the case of an incompressible gas [1, 2] and a compressible gas [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 135–139, May–June, 1976.     

Calculation of autooscillations of the type of azimuthal waves occurring at loss of stability of flow of a viscous fluid between concentric cylinders rotating in opposite directions

Starting with the Navier-Stokes equation we use the Lyapunov-Schmidt method to investigate the nature of the loss of stability of Couette flow between cylinders as the Reynolds number passes through its critical value. We consider the rotation of the cylinders in opposite directions with the ratio of the angular velocities such that the role of the most dangerous disturbances passes over from rotationally symmetric to nonrotationally symmetric disturbances. Branching nonstationary secondary flows (autooscillations) are found in the form of azimuthal waves; the longitudinal wave number and the azimuthal wave number m are assumed given. The amplitude of autooscillations and the wave velocity are calculated for m = 1, and it is shown that depending on the value of both weak excitation of stable and strong excitation of unstable autooscillations are possible and the wave number for which the critical Reynolds number is a minimum corresponds to a stable wave regime in the supercritical region. The linear problem of the stability of the circular flow of a viscous fluid with respect to nonrotationally symmetric disturbances is discussed in [1–3]. Di Prima [1] solved the problem numerically by the Galerkin method when the gap is small and the cylinders rotate in the same direction. Di Prima's analysis is extended in [2] to cylinders rotating in opposite directions, and in [3] it is extended to gaps which are not small. The nonlinear stability problem is treated in [4], where for fixed = 3 and cylinders rotating in opposite directions the axisymmetric stationary secondary flow the Taylor vortex is calculated. The formation of azimuthal waves in the fluid between the cylinders was studied experimentally in detail by Coles [5].Translated from Zhurnal Prikladnoi Mekhanika i Tekhnicheskoi Fiziki, No. 2, pp. 68–75, March–April, 1976.     

Rayleigh stability of shear flow in relaxing media

The stability of steady-state flow is considered in a medium with a nonlocal coupling between pressure and density. The equations for perturbations in such a medium are derived in the linear approximation. The results of numerical integration are given for shear motion. The stability of parallel layered flow in an inviscid homogeneous fluid has been studied for a hundred years. The mathematics for investigating an inviscid instability has been developed, and it has been given a physical interpretation. The first important results in flow stability of an incompressible fluid were obtained in the papers of Helmholtz, Rayleigh, and Kelvin [1] in the last century. Heisenberg [2] worked on this problem in the 1920's, and a series of interesting papers by Tollmien [3] appeared subsequently. Apparently one of the first problems in the stability of a compressible fluid was solved by Landau [4]. The first investigations on the boundary-layer stability of an ideal gas were carried out by Lees and Lin [5], and Dunn and Lin [6]. Mention should be made of a series of papers which have appeared quite recently [7–9]. In all the papers mentioned flow stability is investigated in the framework of classical single-phase hydrodynamics. Meanwhile, in recent years, the processes by which perturbations propagate in media with relaxation have been intensively studied [10–12].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 87–93, May–June, 1976.     

Calculation of unsteady supersonic flow past a two-dimensional cascade of plates ween vorticity inhomogeneities of the flow act on it

In the framework of the linear theory of small perturbations the problem of unsteady subsonic flow past a two-dimensional cascade of plates has been considered in a number of papers. Thus, the unsteady aerodynamic characteristics of a cascade of vibrating plates were calculated in [1] by the method of integral equations, while the same method was used in [2, 3] to calculate the sound fields that are excited when sound waves Coming from outside or vorticity inhomogeneities of the oncoming flow act on the cascade. The problem of a two-dimensional cascade of vibrating plates in a supersonic flow was solved in [4, 5]. In [4] the solution was constructed on the basis of the well-known solution of the problem of vibrations of a single plate, while in [5] a variant of the method of integral equations was used which differed slightly from the usual formulation of this method [1–3]. The approach proposed in [5] is used below to calculate the unsteady flow past a two-dimensional cascade of plates in the case when vorticity inhomogeneities of a supersonic oncoming flow act on it. Equations are obtained for the strength of the unsteady pressure jumps arising in such a flow and the vortex wakes shed from the trailing edges of the plates. Examples of the calculations illustrating the accuracy of the method and its possibilities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp, 152–160, May–June, 1986.     

Investigation of hypersonic three-dimensional flow past blunt bodies within the framework of the parabolized Navier-Stokes equations

The flow of a homogeneous gas in a three-dimensional hypersonic viscous shock layer, which includes the shock wave structure, is examined within the framework of the parabolic approximation of the Navier-Stokes equations. The Navier-Stokes equations are simplified on the basis of the asymptotic analysis carried out in [1], are written in variables of the Dorodnitsyn type [2] and are solved by the method proposed in [3, 4] extended to the case of three-dimensional flows. The flow at zero angle of attack past elliptic paraboloids, two-sheeted hyperboloids and triaxial ellipsoids is calculated. Some results of investigating the flow past such bodies are presented. Flow past a sphere in the analogous approximation was considered in [5], where a comparison was also made with the solution of the complete Navier-Stokes equations [6, 7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 134–142, July–August, 1987.In conclusion, the authors wish to express their warm thanks to V. V. Lunev and G. A. Tirskii for useful discussion and valuable comments.     

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.     

Hydrodynamic flows and heat transfer in slightly noncylindrical channels

Viscous incompressible laminar flow and heat transfer in channels with a small arbitrary deviation from a cylindrical surface are examined. A linear system of equations and boundary conditions for the disturbed dynamic and thermal fields, obtained by linearizing the complete system of Navier-Stokes equations with respect to the solution for developed flows in cylindrical tubes of arbitrary cross section, is presented. In the important practical case in which the perturbations of the channel surface are concentrated on an interval of finite length it is shown that the integral dynamic and thermal characteristics of the channel can be found without solving the three-dimensional equations by going over to effective two-dimensional boundary-value problems which are fundamentally no more difficult to solve than those for developed flows. Extensions of the theory to flows with low-efficiency power sources are given. Applications to plane channels and circular tubes with deformed surfaces are considered. Among the numerous applications requiring information about the integral characteristics of flows in channels whose initially cylindrical surface is slighty deformed, we note the problem of heat transfer intensification by slightly deforming the tube surface with careful estimation of the accompanying increase in resistance [1] and the calculation of the resistance of capillaries and biological transport systems in the form of tubes and channels when the walls are deformed [2]. Below we consider laminar flow in channels with deformed walls. Whereas for the first problem this class of flows is only one of those possible (in general it is necessary to analyze the transition, turbulence and flow separation effects), in the second case, which is characterized by low Reynolds numbers, the laminar flow model is perfectly adequate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 21–30, March–April, 1990.The authors are grateful to A. Yu. Klimenko for useful discussions.     

Direct numerical simulation of the laminar-turbulent transition in a thick spherical layer

The results of a numerical study of the laminar-turbulent transition in unsteady isothermal three-dimensional flows of viscous incompressible fluid in a thick spherical layer between counter-rotating spherical boundaries are presented. The calculations are performed for the governing parameters corresponding to the experimental data [1, 2]. The numerical investigations include both solving the complete system of Navier-Stokes equations and analyzing the linear stability of steady-state axisymmetric flows with respect to three-dimensional disturbances. A stochastic flow regime is calculated for the first time. The limits of existence of different flow regimes and the hysteresis regions are found. The spatial flow patterns and frequency characteristics are obtained, which makes it possible to extend and refine the existing experimental data.     

Numerical simulation of instabilities and secondary regimes in spherical couette flow

In all previous numerical investigations of spherical Couette flow only axisymmetric regimes were considered. At the same time, in experiments [1–4] it was found that when both spheres rotate and the layer is thin centrifugal instability of the main flow leads to the appearance of nonaxisymmetric secondary flows of the azimuthal traveling wave type. The results of an initial numerical investigation of these flows are presented below. Solving the linear problem of the stability of the main flow and simulating the secondary flows on the basis of the complete nonlinear Navier-Stokes equations has made it possible to supplement and explain many of the results obtained experimentally. The type of bifurcation and the structure of the disturbances whose growth leads to the appearance of three-dimensional nonstationary flows are determined, and the transitions between different secondary regimes in the region of weak supercriticality are described.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 3–15, January–February, 1995.     

Development and stability of swirling flows

The purpose of this research is to investigate steady axisymmetric swirling flows in channels and free vortices and also to establish the role of hydrodynamic instability in the formation of the sharp changes in flow structure associated with an increased rate of rotation. On the basis of numerical solutions of the complete Navier-Stokes equations obtained by a finite-difference method swirling flows in pipes with impermeable and permeable walls and in a free vortex are investigated. The stability of the swirling axisymmetric flows is considered on the assumption of local parallelism: the problem of the normal modes developing against the background of the axisymmetric flow determined by the velocity profiles in local cross sections of the flow is solved. Attention is mainly concentrated on free vortex flows with reverse current zones, their structure and stability.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–11, July–August, 1988.     

Numerical investigation of viscous gas secondary flows in a rotating cylinder with sources and sinks

A source-sink model of secondary flow excitation in a rotating cylinder, which describes the interaction between a circulator and a rotating gas, is proposed for a nonlinear system of Navier-Stokes equations, and the results of a numerical calculation of the resulting circulating flows are presented. The modified Newton's method employed in the numerical solution makes use of regularizing perturbations to ensure its stability and convergence at low Ekman numbers and high rates of rotation of the cylinder. The combined effect of mechanical and thermal means of flow excitation and the influence of viscous energy dissipation are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 39–44, July–August, 1989.     

Stability of flows of the Hamel type in channels with permeable walls

The stability of nonparallel flows of a viscous incompressible fluid in an expanding channel with permeable walls is studied. The fluid is supplied to the channel through the walls with a constant velocity v₀ and through the entrance cross section, where a Hamel velocity profile is assigned. The resulting flow in the channel depends on the ratio of flow rates of the mixing streams. This flow was studied through the solution of the Navier—Stokes equations by the finite-difference method. It is shown that for strong enough injection of fluid through the permeable walls and at a distance from the initial cross section of the channel the flow approaches the vortical flow of an ideal fluid studied in [1]. The steady-state solutions obtained were studied for stability in a linear approximation using a modified Orr—Sommerfeld equation in which the nonparallel nature of the flow and of the channel walls were taken into account. Such an approach to the study of the stability of nonparallel flows was used in [2] for self-similar Berman flow in a channel and in [3] for non-self-similar flows obtained through a numerical solution of the Navier—Stokes equations. The critical parameters *, R*, and Cr* at the point of loss of stability are presented as functions of the Reynolds number R₀, characterizing the injection of fluid through the walls, and the parameter , characterizing the type of Hamel flow. A comparison is made with the results of [4] on the stability of Hamel flows with R₀ = 0.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 125–129, November–December, 1977.The author thanks G.I. Petrov for a discussion of the results of the work at a seminar at the Institute of Mechanics of Moscow State University.     

The stability of steady regimes in an isothermal chemical flow reactor

One of the fundamental problems in the theory of chemical reactors is the determination of the number of steady regimes and their stability. The problem of the number of steady regimes has been considered in many studies, for example, in [1–4]. The stability of a steady regime is usually established from an analysis of the behavior of small perturbations. The corresponding linear boundary-value problem for perturbations has been studied mainly in the limiting cases of ideal mixing and ideal displacement. When account was taken of longitudinal mixing, the only criteria obtained were ones which imposed fairly severe restrictions on the parameters [5]. In the present study numerical analysis is used in order to investigate the stability of steady concentration distributions in an isothermal chemical flow reactor with longitudinal mixing in the case of a single chemical reaction. The eigenvalues were obtained for the Sturm-Liouville problem, which fully characterize the stability for several laws of variation of the chemical reaction rate as a function of the concentration. A knowledge of the eigenvalues is essential, for example, in order to construct the stabilization system proposed in [6] for the unsteady regime.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 179–182, March–April, 1985.     

Existence and nonuniqueness of local separation zones in viscous jets

The plane steady problem of the flow of a viscous wall jet past a smoothed break in the contour of a body is considered. For convenience, the flow in the neighborhood of the junction between two flat plates inclined at an angle to each other is chosen for study. As a result of the small extent of the region investigated the flow field is divided into two layers: the main part of the jet, which undergoes inviscid rotation, and a thin sublayer at the wall, which ensures the satisfaction of the no-slip condition. Particular interest attaches to the flow regime in which the solution in the sublayer satisfies the Prandtl boundary layer equations with a given pressure gradient. A similar problem was studied in [1–4]. The present case is distinguished by the structure of the free interaction region in a small neighborhood of the point of zero surface friction stress. By means of the method of matched asymptotic expansions, applied to the analysis of the Navier-Stokes equations, it is established that the interaction mechanism is that described in [5–7]. As a result, an integrodifferential equation describing the behavior of the surface friction stress function is obtained. A numerical solution of this equation is presented. The range of plate angles on which solutions of the equation obtained exist and, therefore, flows of this general type are realized is determined. The essential nonuniqueness of the possible solutions is established, and in particular attention is drawn to the possible existence of six permissible friction distributions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 38–45, January–February, 1986.The author wishes to thank V. V. Sychev and A. I. Ruban for their useful advice and discussion of the results.     

Effect of supersonic stream parameters on the characteristic time of transient step flow

There have been many studies of viscous compressible gas flow in wakes and behind steps [1–6] in which attention has been focused on the steady-state flow regime. The problem of the supersonic flow of a viscous compressible heat-conducting gas past a plain backward-facing step is considered. The problem is solved numerically within the framework of the complete system of Navier-Stokes equations. The passage of the solution from the initial data to the steady-state regime and the effect of the gas dynamic parameters of the external flow on the characteristic flow stabilization time are investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 137–140, July–August, 1989.     

Stability of unsteady viscous flows

A study is made of the linear stability of plane-parallel unsteady flows of a viscous incompressible fluid: in the mixing layer of two flows, in a jet with constant flow rate, and near a wall suddenly set in motion [1]. The slow variation of these flows in time compared with the rate of change of the perturbations makes it possible to use the method of two-scale expansions [2]. The stability of nonparallel flows with allowance for their slow variation with respect to the longitudinal coordinate was investigated, for example, in [3–6]. The unsteady flows considered in the present paper have a number of characteristic properties of non-parallel flows [1], but in contrast to them are described by exact solutions of the Navier-Stokes equations. In addition, for unsteady planeparallel flows a criterion of neutral stability can be uniquely established by means of the energy balance equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, 138–142, July–August, 1981.I thank G. I. Petrov for suggesting the problem, and also S. Ya. Gertsenshtein and A. V. Latyshev for assisting in the work.