Three-dimensional waves in a boundary layer

As is known [1], two-dimensional waves develop in the boundary layer and then become three-dimensional waves with increase of the Reynolds number R. Since Squire [2] has shown that the linear growth of three-dimensional waves is more intense than that of the two-dimensional, it is natural that the behavior of three-dimensional waves in the boundary layer is explained by nonlinear intersection [3], However, Gaster [4] has noted that although disturbances which increase with time are usually considered, experimentally we observe disturbances which grow in space. (Squire's proof does not extend to this case.) It has been shown that the spatially growing disturbances cannot explain the occurrence of the three-dimensional waves (in the linear formulation).The author wishes to thank his scientific advisor G. I. Petrov and also A. A. Zaitsev for valuable discussions of the study.     

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.     

Waves on down-flowing fluid films: Calculation of resistance to arbitrary two-dimensional perturbations and optimal downflow conditions

Wavy downflow of viscous fluid films is studied. The full Navier-Stokes equations are used to calculate the hydrodynamic characteristics of the flow. The stability of calculated nonlinear waves to arbitrary two-dimensional perturbations is considered within the framework of the Floquet theory. It is shown that, for small values of the Kapitza number, the waves are stable over a wide range of wavelengths and values of the Reynolds number. It is found that, as the Kapitza number increases, the parameter range where nonlinear waves are calculated is divided into a series of alternating zones of stable and unstable solutions. A large number of narrow zones where the solutions are stable are revealed on the wavelength-Reynolds number parameter plane for large values of the Kapitza number. Optimal regimes of film downflow that correspond to the minimum value of average film thickness for nonlinear waves with different wavelengths are determined. The basic characteristics of these waves are calculated in a wide range of Reynolds and Kapitza numbers.     

Velocity-profile deformation in EHD flow in a two-dimensional channel

The study of unipolar-charged fluids in the presence of external and induced electric fields has recently taken on great importance. The characteristics of one-dimensional EGD flows [1, 2] and developed laminar flows of a viscous fluid [3] have been clarified in several studies made in this field. However, the study of three-dimensional flows of such media is actually just beginning. Here, along with the analysis of three-dimensional boundary layers and jets [4], there is considerable interest in the study of spatial (two-dimensional and three-dimensional) EHD flows of an inviscid fluid, since in many engineering devices the zone of interaction of the flow with the electric fields does not exceed a few channel diameters, which makes it possible to neglect viscous effects.In this paper we examine some aspects of two-dimensional EHD flows of a viscous incompressible medium for infinitely large electric Reynolds numbers. The perturbations of the hydrodynamic parameters of the flow downstream from the zone of action of the electrostatic forces are determined. It is shown that in many cases the flow parameters outside this zone may be determined without solving the complete system of EHD partial differential equations.     

Excitation of Tollmien-Schlichting waves py a vibrating section of a surface with flow round it

In the context of the problem of describing the transition of a laminar boundary layer to a turbulent, great interest attaches to the study of susceptibility, i.e., of the reaction of the flow to various external influences, such as acoustic perturbations, surface roughness, vibration of the wall, turbulence of the unperturbed flow, etc. A general property of the effect of the factors mentioned above on the flow in a laminar boundary layer was discovered in experimental and numerical studies and is noted in [1]: in all cases an external forcing perturbation leads to the excitation of normal modes of oscillation in the boundary layer which propagate downstream, namely, Tollmien-Schlichting waves. There is an analytical calculation in [2, 3] of the amplitude of a wave excited by harmonic oscillations of a narrow band on the surface of a plane plate, the Reynolds number having been assumed to be infinitely large, and the frequency of the vibrator corresponding to the neighborhood of the lower branch of the neutral cuirve [4], In [5] the amplitude of the wave of instability generated is calculated by the method of expansion of the solution in a biorthogonal system of eigenfunctions. The amplitudes of the Tollmien-Schlichting waves are calculated below by means of a generalization of the method of [2] for the whole range of Reynolds numbers and frequencies of the vibrator corresponding to the region of instability: for moderate Reynolds numbers the problem is solved numerically, while for large Reynolds numbers an asymptotic solution is constructed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 46–51, July–August, 1987.The author is grateful to M. N. Kogan and V. V. Mikhailov for useful discussions of the results of the study.     

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.     

Some geometric acoustic problems of one-dimensional unsteady and two-dimensional steady flows

The behavior of discontinuities (weak shocks) of the parameters of a disturbed flow and their interaction with the discontinuities of the basic flow in the geometric acoustics approximation, when the variation of the intensity of such shocks along the characteristics or the bicharacteristics is described by ordinary differential equations, has been investigated by many authors. Thus, Keller [1] considered the case when the undisturbed flow is three-dimensional and steady, and the external inputs do not depend on the flow parameters. An analogous study was made by Bazer and Fleischman for the MGD isentropic flow of an ideal conducting medium [2], while Lugovtsov [3] studied the three-dimensional steady flow of a gas of finite conductivity for small magnetic Reynolds numbers and no electric field. Several studies (for example, [4]) have considered the behavior of discontinuities of the solutions from the general positions of the theory of hyperbolic systems of quasilinear equations. Finally, the interaction of weak shocks (or the equivalent continuous disturbances) with shock waves was studied in [5–11].In what follows we consider one-dimensional (with plane, cylindrical, and spherical waves) and quasi-one-dimensional unsteady flows, and also plane and axisymmetric steady flows. Two problems are investigated: the variation of the intensity of weak shocks in the presence of inputs which depend on the stream parameters, and the interaction of weak shocks with strong discontinuities which differ from contact (tangential) discontinuities.The thermodynamic properties of the gas are considered arbitrary. We note that the resulting formulas for the interaction coefficients of the weak and strong discontinuities are also valid for nonequilibrium flow.     

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

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

Flow stability of a conducting liquid flowing down an inclined plane in the presence of a magnetic field

The stability of a steady flow of incompressible, conducting liquid down an inclined plane in the presence of longitudinal and transverse magnetic fields is studied. Solutions of the linearized magnetohydrodynamic equations with corresponding boundary conditions are found on the assumption that the Reynolds number R_g and the wave number are small. It is shown that the longitudinal magnetic field plays a stabilizing role. It is known [1] that the flow of a viscous liquid over a vertical wall is always unstable. In this article it is shown that the instability effect at small wave numbers may be eliminated if the longitudinal magnetic field satisfies the conditions found. The case when the Alfvén number and the wave number are small and the Reynolds number is finite is also examined.     

Bifurcation of the solution to the problem of steady traveling waves in a layer of viscous liquid on an inclined plane

Traveling waves in a viscous liquid flowing down an inclined plane can be described at small and moderate Reynolds numbers by an ordinary differential equation in the thickness of the layer [1, 2]. As the Reynolds number tends to zero, this equation goes over into an equation of third order with quadratic nonlinearity [3]. Periodic solutions of this last equation bifurcating from the plane-parallel solution have been investigated by Nepomnyashchii and Tsvelodub [3–6]. In the present paper, a study is made of the bifurcation of periodic solutions from periodic solutions, namely, an investigation is made of the values of the wave number for which a periodic solution is not unique; a bifurcation equation is derived, the number of bifurcating solutions is found, and their behavior near a bifurcation point is considered; and the bifurcating solutions are continued numerically with respect to a parameter (the wave number) from the neighborhoods of the bifurcation points.     

Double-diffusive convection in a cubical lid-driven cavity with opposing temperature and concentration gradients

A numerical study of three-dimensional incompressible viscous flow inside a cubical lid-driven cavity is presented. The flow is governed by two mechanisms: (1) the sliding of the upper surface of the cavity at a constant velocity and (2) the creation of an external gradient for temperature and solutal fields. Extensive numerical results of the three-dimensional flow field governed by the Navier-Stokes equations are obtained over a wide range of physical parameters, namely Reynolds number, Grashof number and the ratio of buoyancy forces. The preceding numerical results obtained have a good agreement with the available numerical results and the experimental observations. The deviation of the flow characteristics from its two-dimensional form is emphasized. The changes in main characteristics of the flow due to variation of Reynolds number are elaborated. The effective difference between the two-dimensional and three-dimensional results for average Nusselt number and Sherwood number at high Reynolds numbers along the heated wall is analyzed. It has been observed that the substantial transverse velocity that occurs at a higher range of Reynolds number disturbs the two-dimensional nature of the flow.     

Wavy flow of a liquid film in the presence of a cocurrent turbulent gas flow

Wavy downflow of viscous liquid films in the presence of a cocurrent turbulent gas flow is analyzed theoretically. The parameters of two-dimensional steady-state traveling waves are calculated for wide ranges of liquid Reynolds number and gas flow velocity. The hydrodynamic characteristics of the liquid flow are computed using the full Navier-Stokes equations. The wavy interface is regarded as a small perturbation, and the equations for the gas are linearized in the vicinity of the main turbulent flow. Various optimal film flow regimes are obtained for the calculated nonlinear waves branching from the plane-parallel flow. It is shown that for high velocities of the cocurrent gas flow, the calculated wave characteristics correspond to those of ripple waves observed in experiments.     

Linear stability of viscous incompressible flow in a circular viscoelastic tube

Flow stability in rigid tubes has been the subject of much research [1]. The overwhelming majority of authors of both theoretical and experimental studies now conclude that Poiseuille flow in a circular rigid tube is linearly stable. However, real tubes all possess elastic properties, the influence of which has not been investigated in such detail. For certain selected values of the parameters characterizing an elastic tube it has been shown that with respect to infinitesimal axisymmetric perturbations Poiseuille flow in the tube can be unstable [2]. In this case boundary conditions that did not take into account the fairly large velocity gradient of the undisturbed flow near the tube wall were used. The present paper reports the results of a numerical investigation of the linear stability of Poiseuille flow in a circular elastic tube with respect to three-dimensional perturbations in the form of traveling waves propagated along the system (azimuthal perturbation modes with numbers 0, 1, 2, 3, 4, and 5 are considered). It is shown that the elastic properties of the tube can have an important influence on the linear stability spectrum. In the case of axisymmetric perturbations it is possible to detect an instability which, at Reynolds numbers of more than 200, exists only for tubes whose modulus of elasticity is substantially less than that of materials in common use. The instability to perturbations of the second azimuthal mode is different in character, inasmuch as at Reynolds numbers greater than unity it occurs in stiffer tubes. Moreover, as the Reynolds number increases it can also occur in tubes of greater stiffness. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 126–134, November–December, 1986.     

Unsteady three-dimensional viscous shock layer in nonuniform gas flow in the neighborhood of a stagnation point

The combined influence of unsteady effects and free-stream nonuniformity on the variation of the flow structure near the stagnation line and the mechanical and thermal surface loads is investigated within the framework of the thin viscous shock layer model with reference to the example of the motion of blunt bodies at constant velocity through a plane temperature inhomogeneity. The dependence of the friction and heat transfer coefficients on the Reynolds number, the shape of the body and the parameters of the temperature inhomogeneity is analyzed. A number of properties of the flow are established on the basis of numerical solutions obtained over a broad range of variation of the governing parameters. By comparing the solutions obtained in the exact formulation with the calculations made in the quasisteady approximation the region of applicability of the latter is determined. In a number of cases of the motion of a body at supersonic speed in nonuniform media it is necessary to take into account the effect of the nonstationarity of the problem on the flow parameters. In particular, as the results of experiments [1] show, at Strouhal numbers of the order of unity the unsteady effects are important in the problem of the motion of a body through a temperature inhomogeneity. In a number of studies the nonstationary effect associated with supersonic motion in nonuniform media has already been investigated theoretically. In [2] the Euler equations were used, while in [3–5] the equations of a viscous shock layer were used; moreover, whereas in [3–4] the solution was limited to the neighborhood of the stagnation line, in [5] it was obtained for the entire forward surface of a sphere. The effect of free-stream nonuniformity on the structure of the viscous shock layer in steady flow past axisymmetric bodies was studied in [6, 7] and for certain particular cases of three-dimensional flow in [8–11].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 175–180, May–June, 1990.     

Flow reversal of an elastico-viscous liquid between a stationary and a torsionally oscillating disc

The flow of an elastico-viscous liquid contained between two infinite discs, when one is held at rest and the other performs small-amplitude torsional oscillations about their common axis, is considered. The liquid is characterized by equations of state more general than, and containing as special cases, the equations of state used by previous authors who have considered this problem.The phenomenon of flow reversal is examined for large values of the Reynolds number, and the apparently different conclusions of previous authors are explained in terms of their particular choices of material parameters.It is also shown that in general the flow at high frequency is dominated by a particular combination of material parameters.     

A spatial laminar boundary layer on a streamline in conical external flow of a uniform gas in the absence of sources and sinks

Theoretical study of a three-dimensional laminar boundary layer is a complex problem, but it can be substantially simplified in certain particular cases and even reduced to the solution of ordinary differential equations.One such particular case is the flow of a compressible gas on a streamline in conical external flow. The case is of considerable practical importance because the local heat fluxes may take extremal values on such lines.Such flow, except for the conical case, has been examined [1–4], and an approximate method has been given [1] on the basis of integral relationships and a special form for the approximating functions. A numerical solution has been given [2, 3] for such flow around an infinite cylinder. It was assumed in [1–3] that the Prandtl number and the specific heats were constant, and that the dynamic viscosity was proportional to temperature. Heat transfer has been examined [4] near a cylinder exposed to a flow of dissociated air.Here we give results from numerical solution of a system of ordinary differential equations for the flow of a compressible gas in a laminar boundary layer on streamlines in conical external flow, with or without influx or withdrawal of a homogeneous gas. It is assumed that the gas is perfect and that the dynamic viscosity has a power-law temperature dependence.     

Nonlinear equation for amplitude of Tollmien-Schlichting waves in a boundary layer

The nonlinear evolution of a Tollmien-Schlichting wave is analyzed with allowance for the flow being nonparallel in a boundary layer. In contrast to the early work of Zel'man [19], strict allowance is made for the fact that the extent to which the flow is nonparallel is not independent of the Reynolds number — the departure from parallel flow in a boundary layer is small only at large Reynolds numbers. Therefore, an asymptotic theory of Tollmien-Schlichting waves is constructed under the assumption that the Reynolds number tends to infinity.