Effect of free-stream turbulence on surface friction in a turbulent boundary layer

The combined effect of the turbulence intensity , the turbulence scaleL, and the Reynolds number Re^** on the surface friction coefficientc _f in a turbulent boundary layer is studied. The dependence of the relative friction increment on the equivalent turbulence level _cq, which takes into account the simultaneous variation in ,L and Re^**, is determined. The threshold value _cq ^* below which the value ofc _f does not depend on _cq is found.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 65–75, March–April, 1995.     

Effective property models for homogeneous two-phase flows

Using an analogy between thermal conductivity of porous media and viscosity in two-phase flow, new definitions for two-phase viscosity are proposed. These new definitions satisfy the following two conditions: namely (i) the two-phase viscosity is equal to the liquid viscosity at the mass quality = 0% and (ii) the two-phase viscosity is equal to the gas viscosity at the mass quality = 100%. These new definitions can be used to compute the two-phase frictional pressure gradient using the homogeneous modeling approach. These new models are assessed using published experimental data of two-phase frictional pressure gradient in circular pipes, minichannels and microchannels in the form of Fanning friction factor (f_m) versus Reynolds number (Re_m). The published data include different working fluids such as R-12, R-22, argon (R740), R717, R134a, R410A and propane (R290) at different diameters and different saturation temperatures. Models are assessed on the basis minimizing the root mean square error (eRMS). It is shown that these new definitions of two-phase viscosity can be used to analyze the experimental data of two-phase frictional pressure gradient in circular pipes, minichannels and microchannels using simple friction models.     

Characterizing developing adverse pressure gradient flows subject to surface roughness

An experimental study was conducted to examine the effects of surface roughness and adverse pressure gradient (APG) on the development of a turbulent boundary layer. Hot-wire anemometry measurements were carried out using single and X-wire probes in all regions of a developing APG flow in an open return wind tunnel test section. The same experimental conditions (i.e., T _∞, U _ref, and C _p) were maintained for smooth, k ⁺ = 0, and rough, k ⁺ = 41–60, surfaces with Reynolds number based on momentum thickness, 3,000 < Re _θ < 40,000. The experiment was carefully designed such that the x-dependence in the flow field was known. Despite this fact, only a very small region of the boundary layer showed a balance of the various terms in the integrated boundary layer equation. The skin friction computed from this technique showed up to a 58% increase due to the surface roughness. Various equilibrium parameters were studied and the effect of roughness was investigated. The generated flow was not in equilibrium according to the Clauser (J Aero Sci 21:91–108, 1954) definition due to its developing nature. After a development region, the flow reached the equilibrium condition as defined by Castillo and George (2001), where Λ = const, is the pressure gradient parameter. Moreover, it was found that this equilibrium condition can be used to classify developing APG flows. Furthermore, the Zagarola and Smits (J Fluid Mech 373:33–79, 1998a) scaling of the mean velocity deficit, U _∞δ*/δ, can also be used as a criteria to classify developing APG flows which supports the equilibrium condition of Castillo and George (2001). With this information a ‘full APG region’ was defined.     

Experimental evaluation of pressure drop in round tubes provided with physically separated, multiple, short-length twisted tapes

Isothermal pressure drop tests were performed on horizontal round tubes each containing identical, physically separated, equally spaced, short-length twisted tapes (TTs). The tests investigated the dependence of the Darcy friction factor, f, on the empty tube Reynolds number Re, TT twist ratio y, and TT spacing s. The variation of f across TTs belonging to the same test section was also examined. Ranges of the experimental variables examined were: 10,000 ? Re ? 90,000, 1.5 ? y ? 6, and s = 30, 40 and 50. The number of 360° revolutions was held constant for all the tapes and equal to 1.5. Tap water at room temperature and nearly atmospheric pressure was used as the working fluid. A correlation for the Darcy friction factor, in the form f = f (Re, y, s), was developed from the collected experimental data, with excellent accuracy.     

Temporal stability of a particle-laden jet

The temporal instability of a particle-laden jet was investigated numerically which took into consideration the parametric effects of jet parameter, B, jet Reynolds number, Re_j, particle mass loading, Z and Stokes number, St. The linear stability theory was used to derive the instability equations of a viscous particle-laden jet flow. The single-phase instability of a top-hat jet was then calculated and compared with the available analytical theories. The numerical results agree well with the analytical results for both the axisymmetric (n = 0) and first azimuthal (n = 1) modes. The results show that the first azimuthal mode disturbance is usually more unstable than that of the axisymmetric mode. But the axisymmetric mode disturbance can be more unstable when Z is high enough (i.e., Z ? 0.1). The higher B and Re_j are, the more unstable the particle-laden jet will be. The existence of particles enhances the flow stability. With the increasing of Z, the jet flow will grow more stable. The inviscid single-phase jet is the most unstable. The wave amplification, c_i first decreases with the increasing of St and then increases afterwards. There exist certain values of St, at which the jet is the most stable.     

Laboratory investigation and commercial test for rotors-assembled strand applied in smooth tube

In the commercial test for smooth tube inserted with rotors-assembled strand comparing with non-inserted ones on condensers in electric power plant, using water as working fluid, the single-phase pressure drop and heat transfer were measured. It was found difficult to receive reliable and accurate enough data through commercial test. Meanwhile, the single-phase pressure drop and heat transfer in a rotors-assembled strand inserted tube were measured in laboratory, with the tube side Prandtl numbers varying from 5.67 to 5.80 and the tube side Reynolds numbers varying from 21,300 to 72,200. Before that, a validation experiment based on the same smooth tube was carried out to testify the experimental system and the data reduction method, in which fixed mounts were employed to eliminate entrance effects. The Prandtl numbers varied from 5.64 to 5.76 and the Reynolds numbers varied from 19,000 to 56,000 in the tube. The annular side Reynolds numbers remained nearly constant at the value of around 50,000 for all experiments, with the annular side Prandtl numbers varying from 8.02 to 8.22. The experimental results of smooth tube show that employment of fixed mounts leads to a visible bias of friction factor at relative low Reynolds numbers while it hardly affects the Nusselt numbers. On the other hand, experiment for the tube inserted with rotors-assembled strand show remarkable improvement for heat transfer with the Nusselt number increased by 9.764–11.87% and the overall heat transfer coefficient increased by 7.08–7.49% within the range of Reynolds number from about 21,300 to 55,500. Meanwhile, friction factor increases inevitably by 278.1–353.9% within the same range of Reynolds number. Based on through multivariant linear normal regression method, the Reynolds number and Prandtl number dependencies of the Nusselt number and friction factor were determined to be Nu = 0.0031Re^0.9Pr^1.0849 and f = 0.993Re^−0.22.     

Mixing characteristics of a film-exciting flapping jet

We have recently discovered a new type of self-excited flapping jets due to a flexible film whose leading edge is fixed at the nozzle exit [Exp Ther Fluid Sci, 106, 226-233]. This paper is to report the experimental investigation on mixing characteristics of the jet induced by a rectangular FEP film. Hot wire anemometry and flow visualization are used to examine the flapping jet flow versus the non-flapping counterpart. Experiments are conducted under the following conditions: i.e., L/D = 1.0 (fixed), W/D = 0.03 ~ 1.0 (varying) and Re = 10000 ~ 45000 (varying); where W and L are the film's width and length, D is the nozzle-exit diameter, and Re is the Reynolds number defined by Re U_oD/ν with U_o and ν being the jet-exit velocity and fluid viscosity.It is found that the jet-flapping frequency f_F varies with W in a complex fashion while it grows roughly linearly with increasing U_o for W/D ≥ 0.5. The flapping Strouhal number St_F f_FD/U_o ranges in 0.13 ≤ St_F ≤ 0.23 for Re = 15,000 ~ 45,000. These Strouhal numbers are substantially lower than that (≈ 0.45 ~ 0.7) for the primary vortex generation in the free jet, but one to two orders of magnitude higher than those from the conventional self-exciting fluidic devices. In general, the flapping jet decays and spreads more rapidly than does the free jet. As W increases, the decaying and spreading rates both grow. Of significance, the centerline evolutions of Taylor and Kolmogorov scales versus the integral scale are examined to characterize the small scales of turbulence against the large-scale motion.     

Delay in response of turbulent heat transfer against acceleration or deceleration of flow in a pipe

To predict the heat transfer enhancements that result from the application of a pulsating flow in a pipe, we experimentally investigated the turbulent heat transfer variations produced in response to sudden accelerations or decelerations to flows within a pipe. To accomplish this, the Reynolds numbers with the valve open (Re₁) and close (Re₀) were systematically varied in the range of 8,000 ≤ Re₁ ≤ 34,000 and 700 ≤ Re₀ ≤ 23,000, respectively, and in-pipe spatiotemporal heat transfer variations were measured using infrared thermography simultaneously with temporal variations to the in-pipe flow properties. Based on the experimental results, it was found that the heat transfer delays that occur in response to accelerations or decelerations can be characterized using the corresponding time lag Δt and first-order time constant τ. The values of Δt and τ can be expressed as non-dimensional forms of Δt/(ν/u_τ²) and τ/(R/u_τ), respectively, where u_τ is the pipe wall friction velocity, ν is the kinematic viscosity of the fluid, and R is the pipe radius.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.     

Polymer Induced Drag Reduction in a Turbulent Pipe Flow Subjected to a Coriolis Force

The skin friction factor f in a turbulent wall-bounded flow can be greatly reduced by using polymer solutions. In this paper we discuss experimental results on the effect of the Coriolis force on turbulent drag reduction. To study this, a horizontal smooth-walled pipe with internal diameter 25?mm is placed on a horizontal table rotating about its vertical axis. The rotation is made non-dimensional with friction velocity and pipe diameter, to form the Rotation number Ro. For a range of bulk Rotation number (Ro _b ) between 0 and 0.6 for two different Reynolds numbers (Re _b = 15 & 30 × 10³), the pressure drop is measured, from which the average friction factor f is obtained. Additionally the effect of four different polymer concentrations has been investigated. The single-phase results show that the friction factor increases monotonic but gradual with Rotation. With polymer additives a drag reduction is found that increases with concentration, but which is not affected by the rotation.     

A rotating elliptic airfoil in fluid at rest and in a parallel freestream

This paper reports results of DPIV measurements on a two-dimensional elliptic airfoil rotating about its own axis of symmetry in a fluid at rest and in a parallel freestream. In the former case, we examined three rotating speeds (Re _c,Ω = 400, 1,000 and 2,000), and in the later case, four rotating speeds (Ro _c,Ω = 2.4, 1.2, 0.6 and 0.4), together with two freestream velocities (Re _c,u  = 200 and 1,000) and two starting configurations of the airfoil (i.e., chord parallel to (α ₀ = 0°) or normal (α ₀ = 90°) to the freestream). Results show that a rotating airfoil in a stationary fluid produces two distinct types of vortex structures depending on the Reynolds number. The first type occurs at the lowest Reynolds number (Re _c,Ω = 400), where vortices shed from the two edges or tips of the airfoil dissipated quickly, resulting in the airfoil rotating in a layer of diffused vorticity. The second type occurs at higher Reynolds numbers (i.e., Re _c,Ω = 1,000 and 2,000), where the corresponding vortices rotated together with the airfoil. Due to the vortex suction effect, the torque characteristics are likely to be heavily damped for the first type because of the rapidly subsiding vortex shedding, and more oscillatory for the second type due to persistent presence of tip vortices. In a parallel freestream, increasing the tip-speed ratio (V/U) of the airfoil (i.e., decreasing the Rossby number, Ro _c,Ω) transformed the flow topology from periodic vortex shedding at Ro _c,Ω = 2.4 to the generation of a “hovering vortex” at Ro _c,Ω = 0.6 and 0.4. The presence of the hovering vortex, which has not been reported in literature before, is likely to enhance the lift characteristics of the airfoil. Freestream Reynolds number is found to have minimal effect on the vortex formation and shedding process, although it enhances shear layer instability and produces more small-scale flow structures that affect the dynamics of the hovering vortex. Likewise, initial starting configuration of the airfoil, while affecting the flow transient during the initial phase of rotation, has insignificant effect on the overall flow topology. Unfortunately, technical constraint of our apparatus prevented us from carrying out complimentary force measurements; nevertheless, the results presented herein, which are more extensive than those computed by Lugt and Ohring (1977), will provide useful benchmark data, from which more advanced numerical calculations can be carried out to ascertain the corresponding force characteristics, particularly for those conditions with the presence of hovering vortex.     

DNS and LES of turbulent flow in a closed channel featuring a pattern of hemispherical roughness elements

Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) were performed for fully-developed turbulent flow in channels with smooth walls and walls featuring hemispherical roughness elements at shear Reynolds numbers Re_τ = 180 and 400, with the goal of studying the effect of these roughness elements on the wall-layer structure and on the friction factor. The LES and DNS approaches were verified first by comparison with existing DNS databases for smooth walls. Then, a parametric study for the hemispherical roughness elements was conducted, including the effects of shear Reynolds number, normalized roughness height (k⁺ = 10–20) and relative roughness spacing (s⁺/k⁺ = 2–6). The sensitivity study also included the effect of distribution pattern (regular square lattice vs. random pattern) of the roughness elements on the walls. The hemispherical roughness elements generate turbulence, thus increasing the friction factor with respect to the smooth-wall case, and causing a downward shift in the mean velocity profiles. The simulations revealed that the friction factor decreases with increasing Reynolds number and roughness spacing, and increases strongly with increasing roughness height. The effect of random element distribution on friction factor and mean velocities is however weak. In all cases, there is a clear cut between the inner layer near the wall, which is affected by the presence of the roughness elements, and the outer layer, which remains relatively unaffected. The study reveals that the presence of roughness elements of this shape promotes locally the instantaneous flow motion in the lateral direction in the wall layer, causing a transfer of energy from the streamwise Reynolds stress to the lateral component. The study indicates also that the coherent structures developing in the wall layer are rather similar to the smooth case but are lifted up by almost a constant wall-unit shift y⁺ (∼10–15), which, interestingly, corresponds to the relative roughness k⁺ = 10.     

Fluid forces on a very low Reynolds number airfoil and their prediction

This paper presents the measurements of mean and fluctuating forces on an NACA0012 airfoil over a large range of angle (α) of attack (0-90°) and low to small chord Reynolds numbers (Re_c), 5.3 × 10³-5.1 × 10⁴, which is of both fundamental and practical importance. The forces, measured using a load cell, display good agreement with the estimate from the LDA-measured cross-flow distributions of velocities in the wake based on the momentum conservation. The dependence of the forces on both α and Re_c is determined and discussed in detail. It has been found that the stall of an airfoil, characterized by a drop in the lift force and a jump in the drag force, occurs at Re_c ? 1.05 × 10⁴ but is absent at Re_c = 5.3 × 10³. A theoretical analysis is developed to predict and explain the observed dependence of the mean lift and drag on α.     

Experimental studies on heat transfer and friction factor characteristics of Al2O3/water nanofluid in a circular pipe under laminar flow with wire coil inserts

In this paper, fully developed laminar flow convective heat transfer and friction factor characteristics of Al₂O₃/water nanofluid flowing through a uniformly heated horizontal tube with and without wire coil inserts is presented. For this purpose, Al₂O₃ nanoparticles of 43 nm size were synthesized, characterized and dispersed in distilled water to form stable suspension containing 0.1% volume concentration of nanoparticles. The Nusselt number in the fully developed region were measured and found to increase by 12.24% at Re = 2275 for plain tube with nanofluid compared to distilled water. Two wire coil inserts made of stainless steel with pitch ratios 2 and 3 were used which increased the Nusselt numbers by 15.91% and 21.53% respectively at Re = 2275 with nanofluid compared to distilled water. The better heat transfer performance of nanofluid with wire coil insert is attributed to the effects of dispersion or back-mixing which flattens the temperature distribution and make the temperature gradient between the fluid and wall steeper. The measured pressure loss with the use of nanofluids is almost equal to that of the distilled water. The empirical correlations developed for Nusselt number and friction factor in terms of Reynolds/Peclet number, pitch ratio and volume concentration fits with the experimental data within ±15%.     

Lock-on characteristics of a cavity shear layer

This investigation focuses on defining the lock-on regions of a cavity shear layer subject to local periodic excitations. A circular cylinder of small diameter (d=4 mm), located very close to the upstream edge of cavity, is used to generate the local periodic excitations in the form of oscillatory rotation about its center with various angular amplitudes (Δθ) and frequencies (f_e). All the experiments were conducted in a recirculating water channel at three different Reynolds numbers that are based on the momentum thickness at the upstream edge of cavity (Re_θ0=152, 216 and 278). The LDV system and the laser-sheet technique are employed to perform the quantitative velocity measurements and the qualitative flow visualization, respectively. For cavity flows at three Reynolds numbers studied, the resonant lock-on is found to be the primary lock-on region within the range of frequency ratio (f_e/f₀=0.28–2.0). Here f₀ denotes the natural instability frequency of an unexcited cavity shear layer. The frequency bandwidth of resonant lock-on region does increase with increasing excitation amplitudes (Δθ). While the excitation amplitudes are smaller than 5° (Δθ5°), the resonant lock-on region, at Reynolds numbers 216 and 278, distributes asymmetrically about f_e/f₀=1.0 and biases to the high frequency (or large f_e/f₀) side. However, the sidewise expansion of resonant lock-on region is enlarged and the degree of asymmetric distribution is alleviated at large excitation amplitudes (Δθ>5°). The amount of sidewise expansion of the resonant lock-on region biased toward the high-frequency side is more significant at the lowest Reynolds number (152) than those at two higher Reynolds numbers (216 and 278). Besides, there exists a sub-harmonic lock-on region only at the lowest Reynolds number 152. The existence of a sub-harmonic lock-on region clearly reveals that the differential equation governing the self-excited oscillation within a cavity contains the quadratic nonlinear term. Further, at the lowest Reynolds number (152), the sidewise expansion of the sub-harmonic lock-on region is much narrower than that of the resonant lock-on region.     

Visualization of a locally-forced separated flow over a backward-facing step

A laboratory water channel experiment was made of the separated flow over a backward-facing step. The flow was excited by a sinusoidally oscillating jet issuing from a separation line. The slit was connected to a cavity in which water was forced through a rigid pipe by a scotch-yoke system. The Reynolds number based on the step height (H) was fixed at Re _H =1200. The forcing frequency was varied in the range 0.305?St _H ?0.955 at the forcing amplitude A ₀=0.3. Time-averaged flow measurements were made by a LDV system, especially in the recirculating region behind the backward-facing step. To characterize the large-scale vortex evolution due to the local forcing, flow visualizations were performed by a dye tracer method with fluorescent ink. The vortex amalgamation process was captured at the effective forcing frequency (St _H =0.477) for laminar separation. This vortex merging process enhances flow mixing, which leads to the shortening of the reattachment length.