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.     

Turbulent friction in the flow of liquid in tubes and channels

A calculating relationship is presented for turbulent flow; it takes a unique form over the whole cross section of the flow. A relationship is also derived between turbulent friction and the mean velocity profile on the basis of the equation for the maximum turbulent friction, which follows directly from the equation of motion. The proportionality factor in this relationship is obtained with due allowance for twelve boundary conditions relating to the turbulent flow, the mean velocity, and their derivatives. The resultant turbulent-friction profiles agree with the experimental data of Laufer. The profile parameters may be related to the Reynolds number.Leningrad. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 140–145, March–April, 1972.     

Influence of the initial conditions on axisymmetric jets in a parallel flow

An experimental investigation was made of the initial-section flow of axisymmetric helium, air, and freon-12 jets in a parallel air flow for two different velocity profiles at the nozzle exit near the boundary of the jet. In one case, the velocity profile was determined by boundary layers on the nozzle walls; in the other case, it was produced artificially by means of a honeycomb of tubes of variable length. Measurements were made of the profiles of the mean and the pulsation velocity and the temperature. The flow was also photographed. The investigations showed that, depending on the initial conditions, the intensity of mixing of the jets in the initial section at Reynolds numbers Re 10⁴ (calculated using the jet diameter) can change from the level determined by molecular diffusion to the level characteristic of developed turbulent flow. The flow structure in the annular mixing layer also depends strongly on the initial conditions. The observed ordered structures in the mixing layer are related to a section of development of perturbations near the nozzle. The ordered structures are strongly influenced by the effect on the jet of acoustic vibrations from an external source. When the initial velocity profile is produced by the honeycomb, the transition to developed turbulence may be due to the development of long-wavelength perturbations or to the development of small-scale turbulence generated by the flow over the end of the honeycomb.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 18–24, July–August, 1980.We thank V. M. levlev and K. I. Artamonov for assistance and for discussing the work.     

Calculation of three-dimensional weakly expanding flow in jets and channels

The results are given of a calculation of laminar flow in a channel of square section and the motion of a turbulent jet from a cruciform nozzle in an ambient flow. To calculate the secondary flows, the field of the transverse velocity is decomposed into irrotational and solenoidal components. The results of the calculation of the flow in the channel are compared with the calculations of other authors and experimental data. To calculate the flow in the turbulent jet, a one-parameter turbulence model is used, and the influence of the inhomogeneity of the distribution of the longitudinal component of the velocity on the components of the Reynolds stress tensor is taken into account. The results of calculation of the flow in the jet behind a cruciform nozzle are compared with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 36–44, July–August, 1984.     

Symmetries, Invariance and Scaling-Laws in Inhomogeneous Turbulent Shear Flows

An approach to derive turbulent scaling laws based on symmetry analysis is presented. It unifies a large set of scaling laws for the mean velocity of stationary parallel turbulent shear flows. The approach is derived from the Reynolds averaged Navier–Stokes equations, the fluctuation equations, and the velocity product equations, which are the dyad product of the velocity fluctuations with the equations for the velocity fluctuations. For the plane case the results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the centre of a Couette flow and in the centre of a rotating channel flow, and a new exponential mean velocity profile that is found in the mid-wake region of high Reynolds number flat-plate boundary layers. The algebraic scaling law is confirmed in both the centre and the near wall regions in both experimental and DNS data of turbulent channel flows. For a non-rotating and a moderately rotating pipe about its axis an algebraic law was found for the axial and the azimuthal velocity near the pipe-axis with both laws having equal scaling exponents. In case of a rapidly rotating pipe, a new logarithmic scaling law for the axial velocity is developed. The key elements of the entire analysis are two scaling symmetries and Galilean invariance. Combining the scaling symmetries leads to the variety of different scaling laws. Galilean invariance is crucial for all of them. It has been demonstrated that two-equation models such as the k– model are not consistent with most of the new turbulent scaling laws.     

Experimental investigation of flow in a trench

An experimental investigation was made of the flow of a viscous incompressible liquid in a trench of square transverse cross section, using a laser Doppler velocimeter. The investigation was made with two values of the Reynolds number Re, corresponding to laminar and turbulent flow conditions in the channel. The experimental data show that a core with a constant vorticity is formed in the trench, that a jet propagates near the walls of the trench, and that there are secondary eddies in the corners of the trench. The motion of a viscous liquid in a trench of rectangular cross section is part of a broad class of breakaway flows. Experimental data on the investigation of flow in trenches are extremely few. A majority of the existing information is limited to visual observations [1–4]. In [2, 5, 6] the question of the unstable character of flow in trenches was discussed. Quantitative measurements of stable eddy flows in trenches were made in [7–9] using a thermoanemometer, and in [7] measurements were made of the pressure at the bottom and walls of trenches; there are data on the distribution of the velocity in the middle sections of trenches. In [8] the mean velocity, the intensity of the turbulence, and the stress of the turbulent flow were obtained in several sections parallel to the side walls of the trench, In [9] a measurement was made of the velocities also in two cross sections of a trench in which one component of the velocity prevails. A brief analysis of the existing experimental results shows that these data are insufficient to form a detailed representation of the character of flow in a trench.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 76–86, March–April, 1976.     

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

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

Universal Velocity Defect Law for the Turbulent Boundary Layer

Turbulent plane boundary layer flows of an incompressible fluid are considered. A refinement of the known Coles wake law is proposed. This refinement makes it possible to ensure the smooth matching of the turbulent boundary layer velocity profile with the outer flow and to extend the range of validity of the law to the case of large positive pressure gradients. The accuracy of the analytical approximation obtained is verified by comparison with the known experimental equilibrium velocity profiles. Using the approximation proposed, a relation for calculating the cross-sectional distribution of the Reynolds stress in the equilibrium boundary layer is derived. The pressure distributions for which the equilibrium turbulent boundary layer flows are single- and two-valued are distinguished.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, 2005, pp. 89–101.Original Russian Text Copyright © 2005 by Mikhailov.     

4D Magnetic resonance velocimetry for mean velocity measurements in complex turbulent flows

An adaptation of a medical magnetic resonance imaging system to the noninvasive measurement of three-component mean velocity fields in complex turbulent engineering flows is described. The aim of this paper is to evaluate the capabilities of the technique with respect to its accuracy, time efficiency and applicability as a design tool for complex turbulent internal geometries. The technique, called 4D magnetic resonance velocimetry (4D-MRV), is used to measure the mean flow in fully developed low-Reynolds number turbulent pipe flow, Re=6400 based on bulk mean velocity and diameter, and in a model of a gas turbine blade internal cooling geometry with four serpentine passages, Re=10,000 and 15,000 based on bulk mean velocity and hydraulic diameter. 4D-MRV is capable of completing full-field measurements in three-dimensional volumes with sizes on the order of the magnet bore diameter in less than one hour. Such measurements can include over 2 million independent mean velocity vectors. Velocities measured in round pipe flow agreed with previous experimental results to within 10%. In the turbulent cooling passage flow, the average flow rates calculated from the 4D-MRV velocity profiles agreed with ultrasonic flowmeter measurements to within 7%. The measurements lend excellent qualitative insight into flow structures even in the highly complex 180° bends. Accurate quantitative measurements were obtained throughout the Re=10,000 flow and in the Re=15,000 flow except in the most complex regions, areas just downstream of high-speed bends, where velocities and velocity fluctuations exceeded MRV capabilities for the chosen set of scan parameters. General guidelines for choosing scanning parameters and suggestions for future development are presented.     

Turbulent air flow over a moving curved surface

A numerical model of turbulent air flow over a curved surface is described. The model is based on two-dimensional nonlinear Reynolds equations and continuity equations written in a coordinate system moving with the profile of the curved surface. The Reynolds stresses are represented in the form of the product of the isotropic turbulent viscosity coefficient, which increases linearly with height, and the deformation tensor of the mean velocity field. Flow over a stationary sinusoidal surface and a sinusoidal gravity wave on water is simulated. The structure of the velocity and pressure wave fields is obtained. The differences in flow over stationary and moving surfaces are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 20–24, September–October, 1986.     

入口速度比对射流拟序结构的影响

为了深入了解湍流流动机理以及湍流拟序结构发现过程的影响因素，本文采用大涡模拟方法对不同入口射流伴流速度比的平面湍射流流动进行了数值模拟。采用分步投影法求解动量方程，亚格子项采用标准Smagorinsky亚格子模式模拟，压力泊松方程采用修正的循环消去法快速求解，空间方程采用二阶精度的差分格式，在时间方向上采用二阶精度的显式差分格式。模拟结果给出了平面射流中湍流拟序结构的瞬态发展演变过程，分析了入口速度比对射流拟序结构发展演化过程及宏观流场形态的影响。为进一步研究射流拟序结构及其在湍流流动中的作用提供了基础。     

Quasihomogeneous model of cavitation flows in diffuser channels

Starting with the experiments carried out by Reynolds in 1894, the flow in Venturi tubes has traditionally been used to study and demonstrate various forms of cavitation. Numerous authors have carried out experimental research on the various flow regimes in diffuser channels [1–7] or have investigated theoretical models of such flows [6, 8]. The occurrence and development of cavitation is closely associated with the phenomenon of turbulent separation complicated by the presence of two-phase flow in the dissipation zone. For a long time these effects were considered separately, until Gogish and Stepanov [9] proposed a single model of cavitation and separation based on the theory of intense interaction of an incompressible potential flow and a turbulent cavitation layer of variable density and embracing the various stages of cavitation. The object of this study is to demonstrate the possibilities of this model with reference to the simple example of flows accompanied by cavitation and separation in plane and axisymmetric diffuser channels of the Venturi tube type with straight and curved walls. The dissipative flow near the walls is described by a quasihomogeneous model of turbulent two-phase flow, in which the presence of two phases is taken into account only by varying the mean density. The potential core of the flow is considered in the one-dimensional formulation. The displacement thickness serves as the flow interaction parameter. The conditions of ocurrence and development of circulatory flows are determined. Examples of symmetrical and nonsymmetrical flows are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 47–54, September–October, 1986.     

Calculation of the momentum and heat transfer in turbulent gas-particle jet flows

The equations for the second moments of the dispersed-phase velocity and temperature fluctuations are used for calculating gas-suspension jet flows within the framework of the Euler approach. The advantages of introducing the equations for the second moments of the particle velocity fluctuations has previously been quite convincingly demonstrated with reference to the calculation of two-phase channel boundary flows [9–11]. The flows considered below have a low solid particle volume concentration, so that interparticle collisions can be neglected and, consequently, the stochastic motion of the particles is determined exclusively by their involvement in the fluctuating motion of the carrier flow. In addition to the equations for the turbulent energy of the gas and its dissipation, the calculation scheme includes the equations for the turbulent energy and turbulent heat transfer of the solid phase; however, the model constructed does not contain additional empirical constants associated with the presence of the particles in the flow.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.3, pp. 69–80, May–June, 1992.     

Turbulent convection in a plane horizontal layer

The problem of convection in an incompressible fluid between two horizontal planes maintained at a constant temperature without friction on the boundaries is considered. The medium is assumed to be turbulent. A theoretical model is constructed using mathematical modeling of the coherent structure in the turbulent flow. This turbulent convection-model has one empirical constant in the relations closing the generalized Reynolds equations. The problem formulated is solved analytically by means of the Stuart-Landau method. The main characteristics of the finite-amplitude ordered convection are obtained and their dependence on the empirical constant is studied.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 49–56, November–December, 1993.     

Heat transfer enhancement with deformation of the velocity profile

The problem of enhancing the heat transfer in channels and boundary layers by the appropriate deformation of the fluid velocity profile is considered. The resulting additional hydraulic losses, the price of heat transfer enhancement, are determined. The possibilities of controlling heat transfer by redistributing the fluid velocity in channels are demonstrated with reference to flows at low Prandtl numbers. Laminar and turbulent liquid and gas flows with heat transfer in channels and boundary layers are numerically modeled on the basis of modern models of turbulence (flow development in channels with different initial velocity profiles, flows with wall roughness and boundary layer flows with forces acting on the flow to cause deformation of the velocity profile). In all cases it is found that the heat transfer can be enhanced only at the expense of a considerable increase in the hydaulic losses. A class of self-similar thermal problems for flows in plane diffusers is formulated. The eigenfunctions — temperature modes — for various velocity profiles are determined with allowance for the nonuniqueness of the solution of the classical dynamical problem for a plane diffuser and the corresponding heat transfer coefficients are found.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 94–105, May–June, 1993.The authors are grateful to A. Yu. Klimenko for useful discussions.     

Intense Reynolds-stress events in turbulent ducts

The aim of the present work is to investigate the role of intense Reynolds shear-stress events in the generation of the secondary flow in turbulent ducts. We consider the connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensional quadrant analysis proposed by Lozano-Durán et al. (J. Fluid Mech., vol. 694, 2012, pp. 100–130). We examine both the geometrical properties of these structures and their contribution to the mean in-plane velocity components, and we perfom a comparison with turbulent channel flow at similar Reynolds number. The contribution to a certain mean quantity is defined as the ensemble average over the detected coherent structures, weighted with their own occupied volume fraction. In the core region of the duct, the contribution of intense events to the wall-normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three-dimensional objects do not differ significantly in both flows. In the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity. However, such contribution is less relevant than that of the complementary portion of the flow not included in such objects. Our results show that strong Reynolds shear-stress events are affected by the presence of a corner but, despite the important role of these structures in the dynamics of wall-bounded turbulent flows, their contribution to the secondary flow is relatively low, both in the core and in the corner.     

Turbulence in a three‐dimensional wall‐bounded shear flow

A new turbulent flow with distinct three‐dimensional characteristics has been designed in order to study the impact of mean‐flow skewing on the turbulent coherent vortices and Reynolds‐averaged statistics. The skewing of a unidirectional plane Couette flow was achieved by means of a spanwise pressure gradient. Direct numerical simulations of the statistically steady Couette–Poiseuille flow enabled in‐depth explorations of the turbulence field in the skewed flow. The imposition of a modest spanwise gradient turned the mean flow about 8° away from the original Couette flow direction and this turning angle remained nearly the same over the entire cross section. Nevertheless, a substantial non‐alignment between the turbulent shear stress angle and the mean velocity gradient angle was observed. The structure parameter turned out to slightly exceed that in the pure Couette flow, contrary to the observations made in some other three‐dimensional shear flows. Coherent flow structures, which are known to be associated with the Reynolds shear stress in near‐wall regions, were identified by the λ₂‐criterion. Instantaneous and ensemble‐averaged vortices resembled those found in the unidirectional Couette flow. In the skewed flow, however, the vortex structures were turned to align with the local mean‐flow direction. The conventional symmetry between Case 1 and Case 2 vortices was broken due to the mean‐flow three‐dimensionality. The turning of the coherent vortices and the accompanying symmetry‐breaking gave rise to secondary and tertiary turbulent shear stress components. By averaging the already ensemble‐averaged shear stresses associated with Case 1 and Case 2 vortices in the homogeneous directions, a direct link between the educed near‐wall structures and the Reynolds‐averaged turbulent stresses was established. These observations provide evidence in support of the hypothesis that the structural model proposed for two‐dimensional turbulent boundary layers remains valid also in flows with moderate mean three‐dimensionality. Copyright © 2009 John Wiley & Sons, Ltd.     

Order and disorder in particle–liquid flows in a pipe

The spherical expanded polystyrene particle–oil two-phase flow in a vertical pipe was used to simulate the dispersed phase distribution in laminar bubbly flows. A three-dimensional particle image tracking technique was used to track the particles in the flow to study the ordered structure of dispersed phase distribution and its transition to disorder. The ordered structures behaved as particle strings aligned in the flow direction as induced by the flow shear. The structures were quite durable in high liquid velocity flows and dispersed gradually as the liquid velocity decreased. In lower velocity flows, the particles tended to form clusters in the horizontal direction, as predicted by potential theory for spherical bubbles rising in a quiescent inviscid liquid and as observed in experiments on non-shear bubbly water flows.     

Two-layer turbulent flow over a rough rotating surface

A model of turbulent incompressible fluid flow over a rough surface under the action of the Coriolis force with a turbulent transfer coefficient corresponding to the Prandtl mixing length is proposed. A solution of the problem, asymptotic in the small Coriolis parameter, is presented for horizontally uniform steady-state flow. It is shown that for a small Coriolis parameter the velocity profile and the turbulent transfer coefficient can differ substantially from the limiting expressions known from Prandtl theory. The smaller the roughness coefficient, the greater the difference.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 61–67, July–August, 1995.     

Coherent Vortex Simulation (CVS), A Semi-Deterministic Turbulence Model Using Wavelets

In the spirit of Ha Minh's semi-deterministic model, we propose a new method for computing fully-developed turbulent flows, called Coherent Vortex Simulation (CVS). It is based on the observation that turbulent flows contain both an organized part, the coherent vortices, and a random part, the incoherent background flow. The separation into coherent and incoherent contributions is done using the wavelet coefficients of the vorticity field and the Biot–Savart kernel to reconstruct the coherent and incoherent velocity fields. The evolution of the coherent part is computed using a wavelet basis, adapted at each time step to resolve the regions of strong gradients, while the incoherent part is discarded during the flow evolution, which models turbulent dissipation. The CVS method is similar to LES, but it uses nonlinear multiscale band-pass filters, which depend on the instantaneous flow realization, while LES uses linear low-pass filters, which do not adapt to the flow evolution. As example, we apply the CVS method to compute a time developing two-dimensional mixing layer and a wavelet forced two-dimensional homogeneous isotropic flow. We also demonstrate how walls or obstacles can be taken into account using penalization and compute a two-dimensional flow past an array of cylinders. Finally, we perform the same segmentation into coherent and incoherent components in a three-dimensional homogeneous isotropic turbulent flow. We show that the coherent components correspond to vortex tubes, which exhibit non-Gaussian statistics and long-range correlation, with the same k ^−5/3power-law energy spectrum as the total flow. In contrast, the incoherent components correspond to an homogeneous random background flow which does not contain organized structures and presents an energy equipartition together with a Gaussian PDF of velocity. This justifies their elimination during the CVS computation to model turbulent dissipation. This revised version was published online in August 2006 with corrections to the Cover Date.