Wall shear stress measurement in a turbulent pipe flow using ultrasound Doppler velocimetry

A turbulent boundary layer of a water flow is investigated by means of pulsed ultrasound Doppler velocimetry. The advantage of this method is the acquisition of complete velocity profiles along the sound propagation line within very short time intervals. The shear stress velocity, used for normalizing the velocity profiles, was determined by fitting the profiles to the universal profiles in a turbulent boundary layer obtained from Prandtl's mixing length theory. A coordinate transformation in the near-wall region is proposed to allocate the velocity data to "true" wall distances. From the experimental values of the wall shear stress velocity, the friction factors for a turbulent pipe flow are calculated and compared to the Blasius law. The overall error in measurement was estimated to NJ.4%.     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

Application of photobleaching molecular tagging velocimetry to turbulent bubbly flow in a square duct

To evaluate turbulence energy budget in bubbly flows, an image processing method in a photobleaching molecular tagging velocimetry is improved for accurate evaluation of velocity gradients. Turbulence properties in single-phase and two-phase dilute-bubbly flows in a square duct are measured using the improved method. As a result, the following conclusions are obtained: (1) The axial velocity and axial turbulent intensity measured by the present method agree well with those measured by laser Doppler velocimetry not only for the single-phase flow but also for the dilute-bubbly flow. (2) The present method can measure velocity components and velocity gradients in the vicinity of the wall, and therefore the present method is of great use in understanding the mechanism of turbulence generation and dissipation near the wall. (3) The present method can provide detailed information on turbulence structure such as turbulence kinetic energy budget. (4) Bubbles tend to increase not only the turbulence production but also the turbulence dissipation.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

High resolution measurement of turbulent structure in a channel with particle image velocimetry

High resolution particle image velocimetry is used to measure the turbulent velocity field for fully developed flow (Re = 2,872) in an enclosed channel. Photographs of particle displacement are obtained in a plane that is parallel to the flow and perpendicular to the walls. These are analyzed to give simultaneous measurements of two components of the velocity at more than 10,000 points. Maps of velocity vectors, spanwise vorticity and Reynolds stress reveal structural aspects of the turbulence. In particular, internal shear layers are observed, in agreement with predictions of direct numerical simulation. Ensemble-averaging of a number of photographs yields statistical properties of the velocity in good agreement with laser-Doppler velocimeter measurements, and with direct numerical simulations.     

Effect of measurement volume size on turbulent flow measurement using ultrasonic Doppler method

This paper presents the results of an investigation on the effects of measurement volume size on the mean velocity profile and the Reynolds stress for fully developed turbulent pipe flows. The study employs the ultrasonic velocity profile method, which is based on the ultrasonic Doppler method. The ultrasonic Doppler method offers many advantages over conventional methods for flow rate measurement in the nuclear power plant piping system. This method is capable of measuring the instantaneous velocity profile along the measuring line and is applicable for opaque liquids and opaque pipe wall materials. Furthermore, the method has the characteristic of being non-intrusive. Although it is applicable to various flow conditions, it requires a relatively large measurement volume. The measurement volume of the present method has a disk-shape determined by the effective diameter of the piezoelectric element and the number of the wave cycles of the ultrasonic pulse. Considering this disk-shaped measurement volume and expressing the time-averaged velocity in a truncated Taylor series expansion around the value at the center of the measuring control volume, the value of the velocity can be obtained. The results are then compared with the data obtained from DNS and LDA measurements. The result shows that the effect of the measurement volume size appears in the buffer region and viscous sublayer.     

Assessment of turbulent spectral bias in laser Doppler velocimetry

A critical evaluation is made of the spectral bias which occurs in the use of a laser Doppler velocimeter (LDV). In order to accommodate the randomly sampled LDV data, statistical treatments of particle arrival times are needed. This is modeled as a doubly stochastic Poisson process which includes the intensity function of the velocity field. Three processing algorithms are considered for spectral estimates: the sample and hold method (SH), the modified Shannon sampling technique (SR), and the direct transform (RG). Assessment is made of these for varying data densities (0.05 ≤ d.d ≤ 5) and turbulence levels (t.i.=30%, 100%). The effects of the values of the Reynolds stress coefficients and the transversal standard deviation on the spectral contents were examined. As an improved version of the spectral estimator, the utility of POCS (the projection onto convex sets) has been tested in the present study. This algorithm is found useful to be in the region when d.d. ? 3.     

Full-field bubbly flow velocity measurements by digital image pulsed laser velocimetry

Digital Pulsed Laser Velocimetry (DPLV) is a full-field, two dimensional, noninvasive, quantitative flow visualization technique. The technique described here includes the novel use of direct digitization of two-phase bubbly flow images using a high resolution imaging system. The image data is stored for further analysis by new image processing and analysis software developed for flow experiments.In the technique, ten consecutive frames of data separated by a time increment of 150 ms, are recorded. Each of these ten frames contains the images of bubbles at that one instant of time. A program smooths the instantaneous image and calculates bubble parameters. Another program matches the bubbles from each of the frames into tracks of bubbles through time. This program uses a statistical technique to determine the best possible path of the bubbles.The ability of pulsed laser velocimetry to capture simultaneous and quantitative rather than qualitative information along with these image processing techniques gives the experimentalist a powerful tool to perform flow visualization and analysis.     

Two-component planar Doppler velocimetry in the compressible turbulent boundary layer

 Non-intrusive Planar Doppler Velocimetry (PDV) has been extended to two-component measurements in a Mach 3, flat plate turbulent boundary layer and the distorted boundary layers downstream of 7° and 14° expansions. PDV results were compared to redundant Laser Doppler Velocimetry (LDV) measurements. Mean velocity results obtained with the two techniques agree to within ±5%. PDV measurements were obtained within 0.4 mm of the surface while LDV could be employed only to within approximately 2.0 mm, highlighting a near-wall resolution advantage for PDV. Effects including those associated with separate filtered and reference cameras led to PDV uncertainties of the same order as the encountered moderate turbulence intensities, precluding an investigation of instantaneous turbulence fields. Despite these difficulties, the current multi-component measurements in distorted, compressible boundary layers highlight the potential of PDV and represent progress in its ongoing evolution. Sources of error and improvements required for quantitative turbulence measurements are discussed. Further advances can be expected from ongoing development efforts. Received: 14 February 1997 / Accepted: 26 August 1997     

Development of a robust image processing technique for bubbly flow measurement in a narrow rectangular channel

This paper presents a robust image processing technique for bubbly flow measurement over a wide range of void fractions. The proposed algorithm combines geometrical, optical and topological information recorded with high speed cameras to separate and reconstruct the overlapping bubbles. The common difficulties such as overlapping, irregular bubble shape, surface deformation and large clustering in digital image processing are solved by combining different information based on a preset decision table and flow chart. Test with synthetic bubble images is performed to evaluate the reliability of the algorithm and quantify the uncertainty of the data. The result shows that the proposed algorithm can accurately measure bubbly flows with void fraction up to 18% for large bubbles. Four runs of bubbly flow images in a 30 mm  ×  10 mm rectangular channel are then recorded by three high speed cameras. The area-averaged void fraction of these test runs range from 2.4% to 9.1%. The axial and lateral distributions of bubble number density are obtained by the present algorithm for studying the characteristics of these flows.     

Taylor length scale measurement by laser Doppler velocimetry

 A new method for making direct measurements of the spatial velocity correlation coefficient, based on two-point laser Doppler velocimetry (LDV), has been developed. In this paper, the effects of control parameters on the correlation coefficient are being investigated. The main sources of experimental error have been identified and analysed. It appears that the probe volume length has a key effect on the accuracy of Taylor micro-scale measurement. A data processing procedure has been established and validated for the determination of this scale. The procedure shows that the portion of the correlation curve used to determine Taylor scale is a function of the integral scale to Taylor micro-scale ratio. Received: 7 June 1995/Accepted: 8 December 1997     

Velocity field in isothermal turbulent bubbly gas-liquid flow through a pipe

Velocity field was measured by laser Doppler velocimetry in isothermal, turbulent bubbly gas-liquid flow through a 26.6 mm inner diameter vertical pipe. The measurements were made about 33 diameters downstream from the pipe entrance, gas injection being just upstream of the entrance. The gas phase radial distribution at the measurement plane exhibited influence of the injection device in that higher gas fraction existed in the central region of the pipe. For comparison, velocity field was also measured in isothermal, turbulent single-phase liquid flow through the same pipe at the same axial plane. Measured were the radial distributions of liquid mean axial and radial velocities, axial and radial turbulent intensities, and axial Reynolds shear stress. The radial distributions of gas bubble mean axial velocity and axial velocity fluctuation intensity were also measured by LDV. A dualsensor fiberoptic probe was used at the same time to measure the radial distributions of gas fraction, bubble mean axial velocity and size slightly downstream of the LDV measurement plane.List of Symbols an average gas bubble diameter - f, f _TP friction factor, friction factor for gas-liquid flow - k _L liquid turbulent kinetic energy - , gas, liquid mass flow rate - R inner radius of pipe - r, {sitR}^* radial coordinate; nondimensional radial coordinate (=r/R) - Re _L liquid Reynolds number - U _G mean axial velocity of gas bubble - U _L mean axial velocity of liquid - U _LO mean axial velocity for flow at the total mass velocity with properties of the liquid phase - u _L ⁺ nondimensional mean axial velocity of liquid in wall coordinate - friction velocity - axial velocity fluctuation intensity of liquid - axial velocity fluctuation intensity of gas bubbles - V_L mean radial velocity of liquid - v _L radial velocity fluctuation intensity of liquid - (uv)_L single-point cross-correlation between axial and radial velocity fluctuations of liquid ( axial Reynolds shear stress) - T _in mean liquid temperature at test section inlet - x flow quality - y normal distance from wall - y ⁺ nondimensional normal distance from wall in wall coordinate (=yu/v_L) - _G gas phase residence time fraction - _L rate of dissipation in the liquid - _L Kolmogorov length scale in the liquid - _L liquid kinematic viscosity - _L characteristic turbulence length scale in the liquid - _G, _L density of gas, liquid - _m gas-liquid mixture density This work was partly supported by National Science Foundation, Thermal Transport and Thermal Processing Program, Chemical and Thermal Systems Division, under Grant No. CTS-9411898.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

Laser Doppler velocimetry measurements in an interior flow test facility: A database for CFD-code evaluation

Measurements in a test facility for a complex interior flow are provided as a database for CFD code evaluation. For pure forced convection as well as to a minor extent also for mixed convection flows, the three time-averaged components of the velocity vector and all components of the Reynolds stress tensor are measured in selected cross sections. Special attention is given to the inflow into the main flow section, since it is the important part of the boundary condition for a numerical solution.     

Blade manipulators in turbulent channel flow

We report here the results of a series of careful experiments in turbulent channel flow, using various configurations of blade manipulators suggested as optimal in earlier boundary layer studies. The mass flow in the channel could be held constant to better than 0.1%, and the uncertainties in pressure loss measurements were less than 0.1 mm of water; it was therefore possible to make accurate estimates of the global effects of blade manipulation of a kind that are difficult in boundary layer flows. The flow was fully developed at the station where the blades were mounted, and always relaxed to the same state sufficiently far downstream. It is found that, for a given mass flow, the pressure drop to any station downstream is always higher in the manipulated than in the unmanipulated flow, demonstrating that none of the blade manipulators tried reduces net duct losses. However the net increase in duct losses is less than the drag of the blade even in laminar flow, showing that there is a net reduction in the total skin friction drag experienced by the duct, but this relief is only about 20% of the manipulator drag at most.List of symbols A, A log law constants - c chord length of manipulator - D drag of the manipulator - dp/dx pressure gradient in the channel - h half height of the channel - H height of the channel (2h) - K log law constant - L length of the channel - L.E. leading edge of the manipulator - P static pressure - P _x static pressure at a location x on the channel - P _xm static pressure at the location x in the presence of manipulator - p _ref static pressure at any reference location x upstream of the manipulator - Re Reynolds number - t thickness of the manipulator - T.E. trailing edge of the manipulator - u velocity in the channel - U friction velocity - U ^* average velocity in the channel - u _c centre-line velocity in the channel - U ₊ U/U _* - u _m velocities in the channel downstream of the manipulators - u _ref velocities in the channel at reference location upstream of the manipulators - w Coles's wake function - W width of channel Also National Aeronautical Laboratory, Bangalore 560 017, India     

Unsteady flow in rotating drums using laser Doppler velocimetry

Non-destructive measurements by laser Doppler velocimetry is employed to study unsteady flow in a hollow drum filled with liquid. The drum is suddenly accelerated from rest or is suddenly decelerated from a steady rotation to rest. Pure water and glycerin-water mixtures are used as the test liquid in which polyethylenelatex particles are mixed as the light scattering tracer. The boundary layer formation, the time history of velocity, momentum and kinetic energy of the liquid, the wall-to-fluid force transfer, and the transient response time are determined. Also determined are the effects of side walls and fluid viscosity on the transient flow response. Of importance is the disclosure of Ekman layer instability near the inner radial wall of the test drum. It is actuated by the centripetal acceleration-induced buoyancy force.List of symbols A wetted surface area of test drum, cm² - a reciprocal of characteristic velocity, = t _sH, s/cm - B width of test drum, cm - b axial coordinate of test drum, cm - D diameter of test drum, cm; D ₁, inner diameter; D ₂, outer diameter - d diameter of laser beam, mm - d _p particle diameter, mm - E kinetic energy of liquid, kg · cm²/s²; E _s, steady value - F force transferred from drum walls to liquid, N - f focal length of lens, mm - G one-half of spacing between two parallel split beams, Fig. 1 - H characteristic length of test drum, cm; = V/A - M momentum of liquid, kg·cm/s; M _s, steady value - m mass of control volume, kg - r radial coordinate of test drum, cm - S fringe spacing, mm - t time, s - t _p time for particle to travel through fringe spacing, s - t _s transient time, s - u liquid velocity, cm/s - V liquid volume in test drum, cm³ - V _s effective volume of sample volume, mm³ - v velocity of tracer particle, cm/s; = S/t - W waist diameter of parabola in Fig. 2, mm - (x, y, z) coordinates for paraboloid in Fig. 2, mm - crossing angle of splitting beams, degrees - wavelength of laser length, cm - v kinematic viscosity, cm²/s - liquid density, kg/cm³ - Doppler frequency, l/s - s at steady state - 1 outer - 2 inner On leave from the Dept. of Mechanical Engineering, Musashi Institute of Technology, Tokyo, Japan     

Measurement of fully-developed turbulent pipe flow with digital particle image velocimetry

A new and unique high-resolution image acquisition system for digital particle image velocimetry (DPIV) in turbulent flows is used for the measurement of fully-developed turbulent pipe flow at a Reynolds number of 5300. The flow conditions of the pipe flow match those of a direct numerical simulation (DNS) and of measurements with conventional (viz., photographic) PIV and with laser-Doppler velocimetry (LDV). This experiment allows a direct and detailed comparison of the conventional and digital implementations of the PIV method for a non-trivial unsteady flow. The results for the turbulence statistics and power spectra show that the level of accuracy for DPIV is comparable to that of conventional PIV, despite a considerable difference in the interrogation pixel resolution, i.e. 32 × 32 (DPIV) versus 256 × 256 (PIV). This result is in agreement with an earlier analytical prediction for the measurement accuracy. One of the advantages of DPIV over conventional PIV is that the interrogation of the DPIV images takes only a fraction of the time needed for the interrogation of the PIV photographs.     

Dynamical eigenfunction decomposition of turbulent channel flow

The results of an analysis of low-Reynolds-number turbulent channel flow based on the Karhunen-Loéve(K-L) expansion are presented. The turbulent flow field is generated by a direct numerical simulation of the Navier-Stokes equations at a Reynolds number Re,= 80 (based on the wall shear velocity and channel half-width). The K-L procedure is then applied to determine the eigenvalues and eigenfunctions for this flow. The random coefficients of the K-L expansion are subsequently found by projecting the numerical flow field onto these eigenfunctions. The resulting expansion captures 90% of the turbulent energy with significantly fewer modes than the original trigonometric expansion. The eigenfunctions, which appear either as rolls or shearing motions, posses viscous boundary layers at the walls and are much richer in harmonics than the original basis functions. Chaotic temporal behaviour is observed in all modes and increases for higher-order eigenfunctions. The structure and dynamical behaviour of the eigenmodes are discussed as well as their use in the representation of the turbulent flow.     

Laminar bubbly flow in an open capillary channel in microgravity

The study of a bubbly laminar two-phase flow in an open capillary channel under microgravity conditions was conducted aboard the sounding rocket, Texus-45. The channel consists of two parallel plates of width b = 25 mm and distance a = 10 mm. The flow along the length l = 80 mm is confined by a free surface on one side and a plate on the opposite side. The bubbles are injected at the nozzle of the capillary channel via six capillary tubes of 100 μm in inner diameter. Different liquid and gas flow rates were tested leading to different liquid free surface shape and bubble size.