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.     

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.     

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.     

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.     

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.     

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 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.     

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.     

Fluid flow in a rotating channel of square section

Laminar flow in a channel rotating about a transverse axis has been studied numerically [1–3] and analytically [4–7] at small Reynolds numbers. The drag coefficient of rotating channels with straight and curvilinear axes has been measured [4, 8, 9]. The present paper gives the results of an experimental investigation into the kinematics of water flow in a channel rotating with different intensities. The flow was visualized by means of hydrogen bubbles and a dye. A study was made of the process of flow separation in a rapidly rotating channel into a core with homogeneous velocity distribution in the direction parallel to the rotation axis and thin shear layers on the walls normal to this axis. The values of the dimensionless numbers were found that correspond to flow rearrangement accompanied by formation of longitudinally oriented vortex structures in the region of higher pressure, and also the values of the rotation parameter needed for the almost complete suppression of turbulence in the region of lower pressure. A general analysis is made of the forms of instability in the different regions of the flow and of the possible flow regimes in a rotating channel.     

Stability of convective flow in a vibration field with respect to three-dimensional perturbations

A plane-parallel convective flow in a vertical layer between boundaries maintained at different temperatures becomes unstable when the Grashof number reaches a critical value (see [1]). In [2, 3] the effect of high-frequency harmonic vibration in the vertical direction on the stability of this flow was investigated. The presence of vibration in a nonisothermal fluid leads to the appearance of a new instability mechanism which operates even under conditions of total weightlessness [4]. As shown in [2, 3], the interaction of the usual instability mechanisms in a static gravity field and the vibration mechanism has an important influence on the stability of convective flow. In this paper the flow stability is considered for an arbitrary direction of the vibration axis in the plane of the layer and the stability characteristics with respect to three-dimensional normal perturbations are determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 116–122, March–April, 1988.     

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.     

Three-dimensional instability of an elliptic Kirchhoff vortex

The instability of a Kirchhoff vortex [1–3] with respect to three-dimensional perturbations is considered in the linear approximation. The method of successive approximations is applied in the form described in [4–6]. The eccentricity of the core is used as a small parameter. The analysis is restricted to the calculation of the first two approximations. It is shown that exponentially increasing perturbations of the same type as previously predicted and observed in rotating flows in vessels of elliptic cross section [4–9] appear even in the first approximation. As distinct from the case of plane perturbations [1-3], where there is a critical value of the core eccentricity separating the stable and unstable flow regimes, instability is predicted for arbitrarily small eccentricity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 40–45, May–June, 1988.     

Circulation flow at the front surface of a sphere in a supersonic wake

Several theoretical and experimental studies of supersonic flow past a blunt body located in the wake behind another body have been made [1–7]. It has been shown that a reverse-circulation flow can occur in the shock layer at the front surface. The possibility of such a flow forming depends on the nonuniformity of the freestream flow and the Reynolds number. This paper presents new results of the theoretical study of the structure of the shock wave at the front surface of such a sphere, obtained on the basis of numerical solution of Navier-Stokes equations. It is shown that for a fixed nonuniformity of the freestream flow, an increase in the Reynolds number and cooling of the surface of the body lead to the formation of a secondary vortex in the region where the contour of the body intersects the axis of symmetry. A study is made of the variations of the drag and heat transfer parameters over the front surface of a cooled and thermally insulated sphere. The possibility of numerical simulation of the flow on the basis of the Euler equations is discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 143–148, May–June, 1985.     

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.     

Numerical study of eccentric Couette–Taylor flows and effect of eccentricity on flow patterns

In this study, the differential quadrature (DQ) method was used to simulate the eccentric Couette–Taylor vortex flow in an annulus between two eccentric cylinders with rotating inner cylinder and stationary outer cylinder. An approach combining the SIMPLE (semi-implicit method for pressure-linked equations) and DQ discretization on a non-staggered mesh was proposed to solve the time-dependent, three-dimensional incompressible Navier–Stokes equations in the primitive variable form. The eccentric steady Couette–Taylor flow patterns were obtained from the solution of three-dimensional Navier–Stokes equations. The reported numerical results for steady Couette flow were compared with those from Chou [1], and San and Szeri [2]. Very good agreement was achieved. For steady eccentric Taylor vortex flow, detailed flow patterns were obtained and analyzed. The effect of eccentricity on the eccentric Taylor vortex flow pattern was also studied.     

Theory of nonequilibrium inductive high-frequency discharge in a gas flow

The global nonequilibrium flow in the discharge chamber of an induction plasma generator is modeled. The problem for an equilibrium discharge was considered in [2, 3]. Here, on the basis of a numerical solution of the combined system of Navier-Stokes, Maxwell, energy, ionization kinetics and electron-gas energy balance equations, the structure of the nonequilibrium discharge is analyzed and the results obtained within the framework of the local one-dimensional approach [1] and on the basis of global numerical modeling of the flow are compared. As distinct from [2, 3], in finding the electromagnetic field distribution in the discharge chamber the boundary-value problem for the two-dimensional Maxwell equations is solved.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 153–160, March–April, 1991.The author is grateful to V. V. Lunev and G. N. Zalogin for their constant interest and useful discussions.