An experimental investigation of the flow about a sphere rotating in a rivlin-ericksen fluid

Tangential and radial velocity profiles were measured for the flow about a sphere rotating slowly in a Newtonian fluid, contained in a rectangular tank. Velocities were determined from enlarged streak photographs of aluminium particles moving in a collimated “sheet” of light, at several planes through the flow field. Similar velocity profiles were measured for the flow of a 1.50% Natrosol 250 H solution about two spheres of different diameters rotating in two different sized rectangular tanks. A set of velocity distributions were also measured for a sphere rotating in a 0.9% Natrosol 250 H solution. A dye tracer study of the flow about a sphere rotating in this liquid is presented as well. Both Natrosol solutions exhibited viscoelastic behaviour. The Newtonian fluid study was carried out at a Reynolds number of 1.2 and the viscoelastic fluid studies were within the Reynolds number range of 0.05–1.24.The zero shear viscosities of the Natrosol solutions were measured using the falling-sphere method. The non-Newtonian material parameters were obtained by fitting the theoretical curves to the measured velocity data. The values of the elastic and shear thinning parameters for the two liquids obtained in the different geometrical and dynamical situations are compared.     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

New measurements of the flow-curves for Carbopol dispersions without slip artefacts

The thickening properties of many commercial thickeners are difficult to measure because of wall slip artefacts. Here we report a series of experiments on a typical thickener where these artefacts have been successfully eliminated. As a result, complete, steady-state flow-curves of aqueous Carbopol 980 (the toxicologically preferred version of the older and more well-known Carbopol 940) dispersions are reported for a range of concentrations of 0.045–1.0 wt%. The vane-and-basket flow geometry was used to avoid slip problems at low shear stress, with the geometry housed in a TA AR1000-N controlled-stress rheometer, whilst a Haake RV2 viscometer with an SV2P and MV2P concentric-cylinder geometries were used at higher shear rates. The flow-curves obtained show a smooth but steep transition from a very high Newtonian viscosity at low shear stress to a much lower viscosity at high shear stress. No real yield stresses were detected, but the higher shear rate results can be fitted to the Herschel-Bulkley model, which assumes an apparent yield stress. The various model parameters are displayed as a function of Carbopol concentration. Received: 29 November 2000/Accepted: 26 February 2001     

Averaged fluid velocities near a cylinder in turbulent flow in an open channel: Experiment

The results of experiments in which a circular cylinder located near the bottom of a rectangular channel was exposed to transverse statistically stationary turbulent subcritical flow with free surface are presented. The particle image velocimetry (PIV) was used to obtain data on the averaged velocity field near the cylinder. The gradients of the longitudinal velocity component were used to determine the shear stresses on the bottom of the channel. It is shown that the presence of the cylinder in the flow causes considerable averaged vertical velocities and a significant change in the shear stresses on the bottom of the channel.     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

Nonlinear three-dimensional stretched flow of an Oldroyd-B fluid with convective condition,thermal radiation,and mixed convection

The effect of non-linear convection in a laminar three-dimensional Oldroyd-B fluid flow is addressed. The heat transfer phenomenon is explored by considering the non-linear thermal radiation and heat generation/absorption. The boundary layer assumptions are taken into account to govern the mathematical model of the flow analysis. Some suitable similarity variables are introduced to transform the partial differential equations into ordinary differential systems. The Runge-Kutta-Fehlberg fourth-and fifth-order techniques with the shooting method are used to obtain the solutions of the dimensionless velocities and temperature. The effects of various physical parameters on the fluid velocities and temperature are plotted and examined. A comparison with the exact and homotopy perturbation solutions is made for the viscous fluid case, and an excellent match is noted. The numerical values of the wall shear stresses and the heat transfer rate at the wall are tabulated and investigated. The enhancement in the values of the Deborah number shows a reverse behavior on the liquid velocities. The results show that the temperature and the thermal boundary layer are reduced when the nonlinear convection parameter increases. The values of the Nusselt number are higher in the non-linear radiation situation than those in the linear radiation situation.     

Study of laminar–turbulent flow transition under pulsatile conditions in a constricted channel

In this paper, direct numerical simulation is performed to investigate a pulsatile flow in a constricted channel to gain physical insights into laminar–turbulent–laminar flow transitions. An in-house computer code is used to conduct numerical simulations based on available high-performance shared memory parallel computing facilities. The Womersley number tested is fixed to 10.5 and the Reynolds number varies from 500 to 2000. The influences of the degree of stenosis and pulsatile conditions on flow transitions and structures are investigated. In the region upstream of the stenosis, the flow pattern is primarily laminar. Immediately after the stenosis, the flow recirculates under an adverse streamwise pressure gradient, and the flow pattern transitions from laminar to turbulent. In the region far downstream of the stenosis, the flow becomes re-laminarised. The physical characteristics of the flow field have been thoroughly analysed in terms of the mean streamwise velocity, turbulence kinetic energy, viscous wall shear stresses, wall pressure and turbulence kinetic energy spectra.     

Two-phase flow modelling: the closure issue for a two-layer flow

The fully-developed, laminar flow of two fluid layers in a horizontal channel is studied by means of the height-averaged balance equations. The closure issue is addressed and the closure relations for the wall and interfacial shear stresses are given for some particular cases.     

Measurements of velocity,temperature and velocity fluctuation distributions in falling liquid films

The velocity, temperature and velocity fluctuation distributions within falling spindle oil films in an inclined rectangular channel were measured using hot-wire techniques and thin thermocouples. The interfacial shear was caused by cocurrent air flow.The results indicate that the liquid films are as a whole much more laminar-like than turbulent in a range of Reynolds numbers (4γ/μ) up to the experimental limit of 6000. Mixing motion occurs in the vicinity of the interface; however, the flow near the wall surface exhibits no sign of such eddy motions, as predicted by the wall law for single phase turbulent flow. Although velocity fluctuation is observed within films with interfacial shear, mean velocity profiles are approximately the same as those obtained by the laminar film prediction.     

Wall functions for numerical modeling of laminar MHD flows

A general wall function treatment is presented for the numerical modeling of laminar magnetohydrodynamic (MHD) flows. The wall function expressions are derived analytically from the steady-state momentum and electric potential equations, making use only of local variables of the numerical solution. No assumptions are made regarding the orientation of the magnetic field relative to the wall, nor of the magnitude of the Hartmann number, or the wall conductivity. The wall functions are used for defining implicit boundary conditions for velocity and electric potential, and for computing mass flow and electrical currents in near wall-cells. The wall function treatment was validated in a finite volume formulation, and compared with an analytic solution for a fully developed channel flow in a transverse magnetic field. For the case with insulating walls, a uniform 20×20 grid, and Hartmann numbers Ha={10,30,100}, the accuracy of pressure drop and wall shear stress predictions was {1.1%,1.6%,0.5%}, respectively. Comparable results were obtained also with conducting Hartmann walls. The accuracy of predicted pressure drop and wall shear stress was essentially independent of the resolution of the Hartmann layers. When applied also to the parallel walls, the wall functions reduced the errors by a factor two to three. The wall functions can be implemented in any general flow solver, to allow accurate predictions at reasonable cost even for complex geometries and nonuniform magnetic fields.     

Radial flow velocity profiles of a yield stress fluid between smooth parallel disks

In rock grouting, idealized 2D-radial laminar flow of yield stress fluids (YSF) is a fundamental flow configuration that is used for cement grout spread estimation. A limited amount of works have presented analytical and numerical solutions on the radial velocity profiles between parallel disks. However, to the best of our knowledge, there has been no experimental work that has presented measured velocity profiles for this geometry. In this paper, we present velocity profiles of Carbopol (a simple YSF), measured by pulsed ultrasound velocimetry within a radial flow model. We describe the design of the physical model and then present the measured velocity profiles while highlighting the plug-flow region and slip effects observed for three different apertures and volumetric flow rates. Although the measured velocity profiles exhibited wall slip, there was a reasonably good agreement with the analytical solution. We then discuss the major implications of our work on radial flow.

     

Unsteady laminar MHD flow and heat transfer in the stagnation region of an impulsively spinning and translating sphere in the presence of buoyancy forces

 An analysis has been carried out to determine the development of momentum and heat transfer occurring in the laminar boundary layer of an incompressible viscous electrically conducting fluid in the stagnation region of a rotating sphere caused by the impulsive motion of the free stream velocity and the angular velocity of the sphere. At the same time the wall temperature is also suddenly increased. This analysis includes both short and long-time solutions. The partial differential equations governing the flow are solved numerically using an implicit finite-difference scheme. There is a smooth transition from the short-time solution to the long-time solution. The surface shear stresses in the longitudinal and rotating directions and the heat transfer are found to increase with time, magnetic field, buoyancy parameter and the rotation parameter. Received on 27 January 2000     

Three-dimensional particle-tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow

Using a three-dimensional (3-D) particle-tracking velocimeter, detailed turbulent flow measurements were made in a plane channel with a one-sided 50% abrupt expansion, which acted as a backward-facing step. The turbulent channel flow reached a fully developed state well upstream of the step. The Reynolds number based on the upstream centerline velocity and the step height H was 5540. With the mean reattachment point located at 6.51H downstream of the step, the measurement region ranged from −2H upstream to 12H downstream of the step. Various turbulent statistics and the energy budget were calculated from numerous instantaneous vector distributions. As in previous experimental investigations, the Reynolds normal and shear stresses had maximum values upstream of the reattachment. The stress anisotropy tensor revealed a peculiar phenomenon near the reattachment wall, wherein the spanwise normal stress was the largest among the three normal stresses. The triple velocity correlations indicated large values in the separating shear layer, and hence the turbulent diffusion was a major term in the energy budget. Comparison was made between the present results and those of the direct numerical simulation (DNS) of Le et al. (1993), and it was found that the mean and fluctuating velocities, the Reynolds shear stress, and the turbulent energy budget were in excellent agreement, although there was a considerable difference in the inflow conditions.     

Particle dynamics calculations of wall stresses and slip velocities for couette flow of smooth inelastic spheres

The objective of this work is to characterize the effects of boundary geometry on the flow of dry granular materials composed of smooth, inelastic spheres between parallel, bumpy walls in the absence of gravity. Particle dynamic simulations are done in which wall stresses and slip velocities are computed over a wide range of parameters, including shear gap height, geometry of the wall particles, wall to free particle diameter ratio, and normal restitution coefficient. Calculated wall stresses and slip velocities are found to be highly sensitive to boundary roughness, which is characterized in terms of the mean spacing value for the wall particles. Most noticeable is a pronounced stress drop for dense flows with an associated large slip velocity in a system having a narrow shear gap height.     

A velocity dependent effective angle method for calibration of X-probes at low velocities

A velocity dependent effective angle (VDEA) method for the calibration of yaw response of hot-wire X-probes at low flow velocities (0.5–6 m/s) is presented. Comparisons with a full velocity vs. yaw-angle method (Österlund 1999) in a smooth wall channel flow indicate that there is only moderate advantage in using the latter method, which is considerably more laborious. Comparisons with direct numerical simulations (DNS) (Moser et al. 1999) and the more common fixed effective angle method (FEA) show that the VDEA method significantly improves estimates of Reynolds stresses compared to the FEA method.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

Turbulent flow in a ribbed channel: Flow structures in the vicinity of a rib

PIV measurements are performed in a channel with periodic ribs on one wall. The emphasis of this study is to investigate the flow structures in the vicinity of a rib in terms of mean velocities, Reynolds stresses, probability density functions (PDF), and two-point correlations. The PDF distribution of u′ is bimodal in the separated shear layer downstream of the rib. The maximum Reynolds shear stresses occur at the leading edge of the rib. Based on quadrant analysis, it is found that ejection motions make a dominant contribution to the Reynolds shear stress in this region. Moreover, topology-based visualization is applied to the separation bubble upstream of the rib. Salient critical points and limit cycles are extracted, which gives clues to the physical processes occurring in the flow.     

Three‐dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model

Turbulent flow in a compound meandering open channel with seminatural cross sections is one of the most complicated turbulent flows as the flow pattern is influenced by the combined action of various forces, such as centrifugal force, pressure, and shear stresses. In this paper, a three‐dimensional (3D) Reynolds stress model (RSM) is adopted to simulate the compound meandering channel flows. Governing equations of the flow are solved numerically with finite‐volume method. The velocity fields, wall shear stresses, and Reynolds stresses are calculated for a range of input conditions. Good agreement between the simulated results and measurements indicates that RSM can successfully predict the complicated flow phenomenon. Copyright © 2008 John Wiley & Sons, Ltd.     

Channel, tube, and Taylor–Couette flow of complex viscoelastic fluid models

We show how to formulate two-point boundary value problems to compute laminar channel, tube, and Taylor–Couette flow profiles for some complex viscoelastic fluid models of differential type. The models examined herein are the Pom-Pom Model [McLeish and Larson 42:81–110, (1998)] the Pompon Model [Öttinger 40:317–321, (2001)] and the Two Coupled Maxwell Modes Model (Beris and Edwards 1994). For the two-mode Upper-Convected Maxwell Model, we calculate analytical solutions for the three flow geometries and use the solutions to validate the numerical methodology. We illustrate how to calculate the velocity, pressure, conformation tensor, backbone orientation tensor, backbone stretch, and extra stress profiles for various models. For the Pom-Pom Model, we find that the two-point boundary value problem is numerically unstable, which is due to the aphysical non-monotonic shear stress vs shear rate prediction of the model. For the other two models, we compute laminar flow profiles over a wide range of pressure drops and inner cylinder velocities. The volumetric flow rate and the nonlinear viscoelastic material properties on the boundaries of the flow geometries are determined as functions of the applied pressure drop, allowing easy analysis of experimentally measurable quantities.     

Effect of system rotating on turbulent boundary layer flow

This paper experimentally investigated the effect of rotating on the turbulent boundary layer flow using hot-wire. The experiments were completed in a rotating rig with a vertical axis and four measured positions along the streamwise direction in channel, which focuses on the flow flied in the rotating channel. The rotating effects on velocity profile, wall shear stress and semi-logarithmic mean velocity profile are discussed in this paper. The results indicated that: due to the Coriolis force induced by rotating, the phenomenon of velocity deficit happens near the leading side. The velocity deficit near the leading side, do not increase monotonically with the increase of Ro. The trend of the velocity deficit near the leading side is also affected by the normal component of pressure gradient, which is another important force in the cross-section of the rotating channel. The wall shear stress near the trailing side is larger than that on the leading side, and the semi-logarithmic mean velocity profile is also different under rotating effects. The phenomenon reveals that the effect of rotation penetrates into the logarithm region, and the flow near the leading side tends to turn into laminar under the effect of rotation. The rotation correction of logarithmic law is performed in current work, which can be used in the wall function of CFD to increase the simulating accuracy at rotating conditions.