3.6 Wind tunnel
4.2.2 Mean velocities
Figure12showsthevelocityprolefortheaxialvelocity
U x
,forbothφ 1 = 90
andφ 1 = 180
.Thevelocityonthecenterlinewasmeasuredasecondtimeaftertheprolewas
taken, the resultwaswithin
1%
of therst measurement forbothproles. Both proleshavethesameshape.Ingure13thedatafor
φ 1 = 90
andφ 1 = 180
isplottedagainstthelogarithmic law. Torbergsen [9] didmeasurementsin the samepipe rig forRe = 75000
and obtained agood t with the logarithmic law usingκ = 0 . 41 andB = 5 . 5
. White [12] claimsthatκ = 0 . 41 andB = 5 . 0
giveabettert toexperimental data. Bothdatamatchesthechoiceof
κ = 0 . 41
. B=5.0assuggestedbyWhitegivesthebest tto themeasurementdatain thelog-lawregion.10 2 10 3
Inafullydevelopedturbulentpipeowonewouldexpectthetransversevelocity
components
U r
andU θ
to be zero. Figures 14 and 15 show that neither of thevelocitycomponentsareexactlyzeroacrossthepipe. Theradialpipevelocity
U r
,isfairlyconstantoverthecrosssectionofthepipe,butshowsomevariation,especially
close tothepipe walls. Therangeof variationin velocityis about
± 0 . 125 m/s
orapproximatelyanangle of
± 0 . 8 ◦
relativeto the averageaxial velocity. Both the proles forφ 1 = 90
andφ 1 = 180 ◦
showthe samekindof variationwith respect
to
y/R
but theyare osetrelativeto one another. Theoset equalsabout1 . 6 ◦
ofprobepitch,whichiswithintheerrorrangeonemustexpectwhentheprobeis
alignedwiththeowvisually.
Thecircumferentialvelocity
U θ
showapeculiar variationoverthecrossection of the pipe, varying overa rangeof± 0 . 45 m/s
or± 2 . 86 ◦
relativeto theaverage
velocity. If thedeviation of the circumferentialvelocity is compared to the local
axialvelocitytherangeofanglevariationexceeds
± 4 ◦
,this isshowningure16.Thevariationof
U θ
withrespecttoy/R
followsthesamepatternforbothseriesof measurementsand closely resembles atypicalinverse tangent function. If the
variation were to be explainedphysically it would implythat the owinside the
pipe wasrotating about the centerline. The velocity doeshowevernot decrease
close to the wall, but increases rapidly, this implies an unphysically large shear
stressonthewall.
Sincetheowisassumedtobeunphysicaltheradialvariationmustbecaused
byoneormoreerrorsin thesetup,datareductionorcausedbylimitationsofthe
−0.2 −0.15 −0.1 −0.05 0 0.05 0.1 0.15
Figure14: Radialmeanvelocity,
U r
−0.5 0 0.5
Figure15: Circumferentialmeanvelocity,
U θ
−5 0 5
Figure 16: Circumferentialmeanvelocity,
U θ
, expressedasangulardeviationrela-tiveto
U x
probe. Thedistinctshapeoftheprolecanmakeonewonderifthisvariationcould
be linked to an error in the arctanexpression in equation 18, but no such error
hasbeenfound. Aninterestingobservationcanhoweverbemadebyobservingthe
dierenceinthetwocenterlinemeasurements,bothfor
U r
andU θ
therelativeerrorbetweenthetwocenterlinepointsislarge. Whiletherepeatedmeasurementsgave
solutionsfortheaxialvelocity
U x
withinonepercentoftherstsolution,theradialand circumferential repeated solutions can vary up to
0 . 5 ◦
and
1 ◦
respectively.
Thisisalargedeviationcomparedtotherangeofthecalculatedvaluesfor
U r
andU θ
. By re-examiningthe measurementdata there was found to be asmall drift in wirevoltage, whichcould notbecorrected forby consideringthe temperaturechange. Ideally themeasurement seriesshould havebeenrepeated, but theerror
wasdiscoveredtolate. Basedonthisobservationpartsofthelargevariationfor
U r
and
U θ
mightbecausedbyvoltagedrift. TheshapeofthevariationofU r
andU θ
doeshoweverseemtobeafunctionofyorasomeotherpropertyrelatedtoy,not
onlyapossiblevoltagedrift. Butwhatpropertycouldthatbe?
U x
isafunctionofy,butissymmtricaboutthecentreline.Thegradientof
U x
alsovariesasafunctionof y but is notsymmtric. Insection 2.2 the possible error ofmeasuring in large
velocitygradientswas dicussed, this couldpossiblybethecause. This discussion
iscontinuedin thenextsectiononshearstresses.
4.2.3 Turbulentshear stresses
As discussed in section 2.3.3, the
u x u r
shear stress is expected to be thedomi-nantshearstressandbehavelinearlyacrossalargeportionofthepipecrossection
accordingtoequation37. Figure17showsthetheoreticalrelationandthe
experi-−1 −0.5 0 0.5 1
Theexperimental resultsfollowthelinear relationwell. The gradientis abit
smallerthan1inthecentreregionwhereviscousshearisassumedtobeneglectable
but thedeviationis small,averysimilar resultwasfoundbyTorbergesen[9]. At
thecentrelinethemeasuredshearstressisclosetozero,asexpected.
Closeto the pipe wallthe turbulent shear stressesare expected to decay and
droptozerointheviscoussublayer. Thespatialresolutionoftheprobeisnotlarge
enoughtomeasure closerto thewallthanabout
y + = 70
whichisfar outsidethe viscoussublayer. Themeasurementdatashowlittletendencytodropofclosetothewall. For
φ 1 = 90
there isalittledropfory/R < − 0 . 9
,whilefory/R > 0 . 9
there is actually anincreasein shear stress measuredfor both datasets. Theincreasedshearstresscouldjust beoutliers,butithappensforbothdatasets.
The two other shear stresses,
u x u θ
andu r u θ
should theoretically be zero as there is no meanvelocity gradientresultingin productionof neitherof them. Ingure18thenormalizedstressesareplotted. Inthecentreregionthemagnitudeof
u x u θ
andu r u θ
arerelativelysmallcomparedtothemaximumvalueofu x u r
,about3
%
, but not zero. Moving closer to the wall both shear stresses increase slowly until|y/R| ≥ 0 . 8
,where theshearstressesincreasemorerapidly. Themagnitude and variation of the shear stresses can not be explained physically, and must berelatedtothemeasurementprocess.
For
|y/R| ≥ 0 . 8
thevelocitygradientexperiencedbytheprobevolumeislarge.The exact eect of an exessively large velocity gradient compared to the probe
volume is unknown. But it will result in calculated velocities dierent from the
truevelocity,as the wallis approached. Theresult canbe agradientofboth
U r
and
U θ
withrespecttoyasobservedinsection4.2.2. Subsequentlythisislikelyto−0.4 −0.3 −0.2 −0.1 0 0.1 0.2 0.3
resultin unphysicalshear stresses. Insection4.2.2 thedeviation forboth
U r
andU θ
wasfoundtoincreaserapidlyfor|y/R| ≥ 0 . 8
,thismatchestheresultfoundforu x u θ
andu r u θ
and supports the theory that a largevelocity gradientbiasestheresult.
WhencomparingtheresultswiththatfoundbyAanesland[1],thedatascatter
fortheprolesobtainedin thisprojectislargeranddonotcollapseasneatlyinto
a line astheresults ofAanesland. What causes this is uncertain,but a to short
samplingtime couldpotentiallybethereason.
4.2.4 Turbulentnormal stresses
Figure 19 displaysthereduced normalstress forbothmeasurementseries, which
revealsthat
( u x ) +
isthelargestnormalstress,asonewould expect.Torbergsen [9]found
u x +
onthe centrelineto beapproximately0.85 for< = 0 . 75 e 5
. TheresultsshouldbecomparableastheReynoldsnumbersisof thesame order. According to the resultsof Torbergsen,u x +
is fairlyconstant onthecen-trelineforincreasingReynoldsnumbers,butincreasesclosertothewallduetothe
increasedvelocitygradient.
u x + cl. = 0 . 85
matchestheobtainedresultsfairlywell, there is however somescatter in the data as is already mentioned in the end ofsection4.2.3. Thedataseriesforthetwodierentprobesalsogiveslightlydierent
results. Movingcloserto thewall
u x +
isunderestimatedcompared tothe results ofTorbergsen,butmatchestheresultsofAanesland[1]better.u r +
andu θ +
showthe same variation asreported by Torbergsen for
|y/R| ≥ 0 . 6
but the scatter is largeforvariationofφ 1
.Valuesfor
u x +
,u r +
andu θ +
cannot beestimated without givinga relativlylargepotentialerror.
u θ +
fory/D = ± 0 . 6
,canforinstancebeestimatedas1 . 1
but−1 −0.5 0 0.5 1
Figure19: Reducednormalstress
thescatter isof theorder0.2, andcould potentiallybelargerifmeasurementsat
morevaluesof
φ 1
weretaken. Thelargevaritioncanbecausedbypooradjustmentof
φ 1
,itisasmentionedearliersetvisually.Ingure20theturbulenceintensityrelativetothelocalstreamwisevelocityis
plotted.
Thestreamwiseturbulence intensity onthe centreline,
u x + ≈ 3 . 3 − 3 . 6
. Tor-bergsenreporteda valueof approximately3 . 5%
. On thecentrelineu r +
andu θ +
should beequaldueto symmetry,there ishoweversomedierencewhich ismost
likelycausedbymisalignementoftheprobe.
4.3 Cylinder wake
Thecylinderwakecanbeanalyzedbothasameanowandasatimevaryingow.
Both approacheswillbetestedin thissection. Ideallyseveralmeasurementseries
for dierentvaluesof
φ 1
should havebeentakento gainmoreinformationabout theproberesponse. Onlyonemeasurement serieswastakenhowever. A velocityprole wastakenin the test sectionbefore thecylinder wasinserted,to mapthe
referencefreestreamconditions.
Vortexsheddingfromacircularcylinderisamuchstudiedow,itdidhowever
provediculttondnearwakeresultsinthesamerangeofReynoldsnumbers
( 10 4 )
, and downstreamdistancex/D.A surveyusingafourwirehotwireprobebyOngand Wallaceat
< = 3900
wastheclosestmatch found[7]. Theyusedafour wire hot wire probe with a crossection of 1 mm x 1 mm. For the ow investigatedin this project
< = 30717
, placing it in the subcritical range together with the resultsofOngandWallace. Theresultscanthereforebeassumedindependentof−1 −0.5 0 0.5 1
Figure20: Turbulenceintensity
Reynoldsnumberandcomparable.
Emphasisis placedonassesmentof thephysical validityof theresult, aswell
ascomparisonwiththeresultsofOngandWallace.