Transition flow of a fibrous suspension in a pipe as the flow of an anisotropic fluid

The transition flow is considered of a fibrous suspension in a pipe. The flow region consists of two subregions: at the center of the flow a plug formed by interwoven fibers and fluid moves as a rigid body; between the solid wall and the plug is a boundary layer in which the suspension is a mixture of the liquid phase and fibers separated from the plug [1–3]. In the boundary region the suspension is simulated as an anisotropic Ericksen—Leslie fluid [4, 5] which satisfies certain additional conditions. Equations are obtained for the velocity profile and drag coefficient of the pipe, which are both qualitatively and quantitatively in good agreement with the experimental results [6–8]. Within the framework of the model, a mechanism is found for reducing the drag in the flow of a fibrous suspension as compared to the drag of its liquid phase.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 91–98, September–October, 1985.     

Influence of the viscous boundary layer on the deposition of particles from a gas suspension flowing past a sphere

The flow of incompressible gas containing particles past bodies of simple shapes at moderate and high velocities is investigated in [1–5], in which the flow of the carrier medium is assumed to be irrotational. The estimates made in [3] for the neighborhood of the stagnation point show that it is necessary to take into account the viscous boundary layer in the case of fine particles. In the present paper, the viscous flow of a gas suspension over the front surface of a sphere at Reynolds numbers R = 10³–10⁷ is considered. It is assumed that the carrier gas is incompressible and the particle concentra ion negligibly small. The influence of the boundary layer on the particle trajectories and the deposition of the disperse phase on the surface of the sphere is investigated. It is shown that there is a wide range of flow parameters for the gas suspension in which the influence of the boundary layer is important. The limits of this range are established.Translated from Izvestiya Akademli Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–66, January–February, 1982.I thank Yu. P. Savel'ev for a helpful discussion of the work.     

Simulation of ordered structures in open turbulent flows

Extensive experimental material [1–4] indicates that ordered (coherent) structures play an important part in determining the nature of the flow, the generation of Reynolds stresses and turbulence energy, and the transport of heat, momentum, and passive admixtures in a turbulent flow. In the present paper, a model is constructed for describing coherent structures in which, given the profile of the mean velocity, one can determine the characteristic sizes, the propagation velocities, and also the frequency and amplitude characteristics of these ordered motions. The model is based on the analogy between the ordered formations and secondary flows in a subsidiary laminar flow whose velocity profile is the same as the turbulent profile of the mean velocity. The influence of small-scale pulsations is described by the introduction of the coefficient of turbulent viscosity. In the framework of the model, numerical calculations are made for two-dimensional turbulent flows in a mixing layer, a jet, and a wake behind a cylinder. The results of the calculations are compared with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 45–52, July–August, 1981.     

Theoretical investigation of coherent structures in stratified turbulent channel flow

A theoretical model is constructed of turbulent stratified flow in a flat horizontal channel with allowance for coherent structures that arise in it. The ordered part is separated from the turbulence of the flow and to describe the Reynolds-type equations are derived. The remaining part of the turbulence is taken into account parametrically in the form of an effective exchange coefficient. The flow is divided into a core, in which the ordered structures are manifested quite clearly, and wall regions, in which ordered large-scale structures are weakly manifested. To study the coherent structures in the core of the flow, an approach analogous to one already used to model ordered structures in open flows [4] is used. Monin-Obukhov scaling theory is used to describe the turbulence in the wall region.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 32–38, May–June, 1991.     

Estimation of the hate of flow in porous media

In a number of problems of the theory of flow in porous media it is particularly important to find the integral characteristics of the flow for regions that constitute extended stream tubes. (Such a region can be imagined as the result of the continuous deformation of a cylinder whose lateral surfaces are impermeable while the bases are surfaces of constant pressure, the inlet and outlet of the flow. In the plane case the cylinder becomes a rectangle.) Usually, the flow rate Q is to be found from the difference of head H. In some cases it possible to obtain upper and lower bounds for the flow rate by varying the flow region, the flow resistance field or the form of the flow law [1–4]. The aim of the present study is to find the shape of the region of fixed volume (in the plane case area) which for given constraints realizes an extremum of the steady-state flow rate. It is shown that the extremality requirement leads to an additional local condition on the unknown part of the boundary. A class of plane problems for which the resulting boundary-value problem with an unknown boundary is effectively solved is identified. Examples are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 80–87, March–April, 1986.     

Contribution to the theory of thin-profile trailing edge separation

The conditions of nonsymmetric trailing edge flow with separation are investigated. Solutions of the equations for the interaction zone in the neighborhood of the trailing edge of a thin profile at an angle of attack of the order O(Re^–1/16) in the separated flow regime are constructed numerically. It is shown that for this zone a solution exists up to a certain angle of attack. In all the regimes the value of the friction on the upper surface at the very end of the trailing edge remains a positive quantity. The solution of the equations in the separated flow regimes is found to be nonunique. The flow over the leading edge is assumed to be unseparated, and the separation at the trailing edge, if present, is assumed to be localized in the interior of the boundary layer. The flow over a Kutta profile at zero angle of attack is taken as an example. In this case the satisfaction of the Chaplygin-Joukowsky condition at the trailing edge ensures smooth flow over both the trailing and leading edges.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 55–59, July–August, 1989.     

Hydraulic modeling of the flow in the delta of the River Neva and in the eastern part of the Gulf of Finland

In connection with the construction of Leningrad's flood defenses the flow in the mouth of River Neva has been investigated on hydraulic models of different scales both unpressurized (see, for example, [1, 2]) and pressurized (see, for example, [3]). Since the results of these investigations are being intensively debated in connection with the discussion of the effect of defensive structures on the ecology of the Leningrad region, it is important to analyze the scientific basis for the hydraulic modeling of such flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 57–65, September–October, 1990.     

Role of localized streamwise structures in the process of transition to turbulence in boundary layers and jets (review)

Results of the analysis of specific features of the laminar–turbulent transition in various subsonic shear flows, which are caused by localized stationary and nonstationary streamwise structures, are presented. One mechanism of flow turbulization is considered, which involves the origination and development of secondary highfrequency disturbances in regions of flow instability generated by its modulation by streamwise structures. It is shown that this process is identical in different types of shear flows (boundary layers and jets) and in flows of the type of localized streamwise structures (stationary or nonstationary).     

Wave propagation and resonance phenomena in inhomogeneous media

The waveguide and resonance properties of inhomogeneous penetrable one–dimensional–periodical structures that consist of two different media are studied within the framework of a one–dimensional approximation. The pass and stop bands are determined. A dispersion relation for all the waveguide modes is obtained. Explicit expressions for low waveguide frequencies and corresponding phase velocities of waveguide modes for mono– and polydisperse media are found. The influence of the polydispersity of the sizes of heterogeneities on the low frequencies of a pass band is considered. A pass band in the range of low frequencies is detected. It is shown that the polydispersity does not affect the waveguide properties of a medium at low frequencies of the first pass band. The resonance phenomena in periodical media and structures are investigated. The resonance phenomena are shown to occur for an unlimited discrete set of frequencies if the group velocity of the waveguide mode for them is zero; in this case, the growth of the oscillation amplitude is localized in the neighborhood of a source (localization of the resonance). A synchrophasotron resonance at which the infinite chain of oscillation sources has the oscillations phase of a corresponding traveling wave from the pass band is detected.     

A Numerical Study of Micro-Heterogeneity Effects on Upscaled Properties of Two-Phase Flow in Porous Media

Commonly, capillary pressure–saturation–relative permeability (P ^c–S–K _r) relationships are obtained by means of laboratory experiments carried out on soil samples that are up to 10–12 cm long. In obtaining these relationships, it is implicitly assumed that the soil sample is homogeneous. However, it is well known that even at such scales, some micro-heterogeneities may exist. These heterogeneous regions will have distinct multiphase flow properties and will affect saturation and distribution of wetting and non-wetting phases within the soil sample. This, in turn, may affect the measured two-phase flow relationships. In the present work, numerical simulations have been carried out to investigate how the variations in nature, amount, and distribution of sub-sample scale heterogeneities affect P ^c–S–K _r relationships for dense non-aqueous phase liquid (DNAPL) and water flow. Fourteen combinations of sand types and heterogeneous patterns have been defined. These include binary combinations of coarse sand imbedded in fine sand and vice versa. The domains size is chosen so that it represents typical laboratory samples used in the measurements of P ^c–S–K _r curves. Upscaled drainage and imbibition P ^c–S–K _r relationships for various heterogeneity patterns have been obtained and compared in order to determine the relative significance of the heterogeneity patterns. Our results show that for micro-heterogeneities of the type shown here, the upscaled P ^c–S curve mainly follows the corresponding curve for the background sand. Only irreducible water saturation (in drainage) and residual DNAPL saturation (in imbibition) are affected by the presence and intensity of heterogeneities.     

Micro-scale instability in a continuously stratified fluid

Experimental observations of small-scale structures caused by flow instabilities at layers of high density gradient in the wake behind a cylinder in a fluid with a continuous salt concentration stratification are reported. In the density wake it is possible to discern a number of structures such as wedge-like structures or cusps; small-scale instabilities (breakers) in the zones of interaction of attached internal waves and the high-gradient wake envelope; small-scale instability of the density boundary layer with a complicated density gradient pattern superimposed on a smooth velocity profile, and small-scale wake structures behind attached vortices in the case of a closed (central) wake envelope.Translated from Izvestiya RossiiskoiAkademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–10, May–June, 1995.     

Interference of steady shock waves traveling in the same direction

A classification of the possible types of shock-wave structures formed as a result of the interference between overtaking shocks in a homogeneous flow is developed on the basis of a previous study [1]. A series of analytic and numerical interaction type criteria is obtained, which makes it possible to justify and supplement the analysis, carried out in [2], of the regions of the governing flow parameters in which steady-state solutions for shock-wave structures of different types exist. The calculations are found to be consistent with the known experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 143–152, July–August, 1987.     

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.     

Two-phase three-component flow through a porous medium with variable total flow

In analyzing the processes of the displacement of oil, in which intensive interphase mass transfer takes place, it is normally assumed that the partial volumes of the components as they mix are additive (Amagat's Law) [1, 2]. Then the equations of motion have an integral, which is the total volume flow rate through the porous medium, and the basic problems of frontal displacement, if there are not too many components in the system, permit an exact analytical study to be made [3–5]. If this assumption is rejected, the total flow becomes variable [3, 6, 7]. It appears that the consequences of this as applied to the processes of the displacement of oil by high pressure gases have not previously been considered. The results of such a study, developing the approach outlined in [4], are given below. The initial multicomponent system is simulated by a three-component one which contains oil (the component being displaced), gas (the neutral or main displacing component), and intermediate hydrocarbon fractions or solvent (the active component). It is shown that instead of the triangular phase diagram (TPD) normally used where the partial volumes of the components are additive, in this case it is convenient to use a special spatial phase diagram (SPD) of the apparent volume concentrations of the components to construct the solutions and to interpret them graphically. The method of constructing the SPD and its main properties are explained. A corresponding graphoanalytical technique is developed for constructing the solutions of the basic problems of frontal displacement which correspond to motions with variable total flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 113–120, May–June, 1985.     

Flow and Mixing Fields for Transported Scalar PDF Simulations of a Piloted Jet Diffusion Flame (‘Delft Flame III’)

Numerical simulation results are presented for ‘Delft Flame III’, a piloted jet diffusion flame with strong turbulence–chemistry interaction. While pilot flames emerge from 12 separate holes in the experiments, the simulations are performed on a rectangular grid, under the assumption of axisymmetry. In the first part of the paper, flow and mixing field results are presented with a non-linear first order k–ε model, with the transport equation for ε based on a modeled enstrophy transport equation, for cold and reactive flows. For the latter, the turbulence model is applied in combination with pre-assumed β-PDF modeling for the turbulence–chemistry interaction. The mixture fraction serves as conserved scalar. Two chemistry models are considered: chemical equilibrium and a steady laminar flamelet model. The importance of the turbulence model is highlighted. The influence of the chemistry model is noticeable too. A procedure is described to construct appropriate inlet boundary conditions. Still, the generation of accurate inlet boundary conditions is shown to be far less important, their effect being local, close to the nozzle exit. In the second part of the paper, results are presented with the transported scalar PDF approach as turbulence–chemistry interaction model. A C₁ skeletal scheme serves as chemistry model, while the EMST method is applied as micro-mixing model. For the transported PDF simulations, the model for the pilot flames, as an energy source term in the mean enthalpy transport equation, is important with respect to the accuracy of the flow field predictions. It is explained that the strong influence on the flow and mixing field is through the turbulent shear stress force in the region, close to the nozzle exit.     

Characteristics of an unstably localized disturbance in a compressed boundary layer

Experiments have demonstrated [1] that the transition of streamline-type flow into turbulent flow in a boundary layer occurs as a result of the formation and development of turbulent spots apparently arising from small natural disturbances. A study of the nonlinear evolution and interaction of localized disturbances requires knowledge of their characteristics to a linear approximation [2]. In the current work, results are presented of calculations of such characteristics for the first two unstable modes in a supersonic boundary layer on a two-dimensional plate (M = 4.5, T_w = 4.44).Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 50–53, January–February, 1976.     

Determination of relative permeabilities in a three-phase flow in a porous medium

A study is made of the problem of determining the parameters of flow described by the Buckley-Leverett system of equations by using functions that admit direct measurement. The well-known solution to the analogous problem for two-phase flow [1–3] is generalized. In contrast to [4], the general case is considered, when the fractions of the phases in the flow and the phase permeabilities depend on two variables.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 187–189, September–October, 1984.The author wishes to thank B. V. Shalimov for his helpful advice.     

Longitudinal-transverse interaction in a three-dimensional boundary layer

Three-dimensional flow is considered for an incompressible fluid in a boundary layer developing along a curved solid surface during interaction between it and a small uneven area (projection or depression) on the surface. It is shown that an important part in the formation of the flow round the uneven area may be played by the drop in the pressure across the boundary layer. Conditions are formulated under which this effect, which is connected with the action of centrifugal forces, is realized. On the assumption that the longitudinal dimension of the uneven area is of the order of Re^–3/14, its width of the order of O(Re^–3/7), and its height O(Re^–4/7), where Re is the Reynolds number, asymptotic equations are derived which describe the motion of the fluid in the neighborhood of the uneven area.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–50, May–June, 1987.     

Development and instability of steady convective flows in a square cavity heated from below and a field of vertically directed vibrational forces

The present paper is devoted to numerical investigation of the spatial structure and stability of secondary vibrational convective flows resulting from instability of the equilibrium of a fluid heated from below. Vibrations parallel to the vector of the gravitational force (vertical vibrations) are considered. As in earlier work [7–9], a region of finite size is used — a square cavity heated from below. It is shown that enhancement of the vibrational disturbance of the natural convective flow may either stabilize or destabilize flows with different spatial structures; it may also stabilize certain solutions of the system of convection equations that are unstable in the absence of vibrational forces. In addition, increase of the vibrational Rayleigh number can lead to a change of the mechanisms responsible for equilibrium instability and oscillatory instability of the secondary steady flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–18, March–April, 1991.I thank G. Z. Gershuni for assistance and extremely fruitful discussions of the results of the paper.     

Effect of Evaporation of Liquid Droplets on the Distribution of Parameters in a Two–Species Laminar Flow

A calculation model was developed, and the heat– and mass–transfer characteristics in a laminar air—vapor—droplet flow moving in a round tube were studied numerically. The distributions of parameters of the two–phase flow over the tube radius were obtained for varied initial concentrations of the gas phase. The calculated heat and mass transfer is compared to experimental data and calculations of other authors. It is shown that evaporation of droplets in a vapor—gas flow leads to a more intense heat release as compared to a one–species vapor—droplet flow and one–phase vapor flow