王漢封

(中南大學(xué) 土木與建筑學(xué)院，湖南 長沙，410075)

有限長度鈍體尾流的三維特性

王漢封

(中南大學(xué) 土木與建筑學(xué)院，湖南 長沙，410075)

運(yùn)用X形熱線，對(duì)高寬比分別為3，5和7的有限長度方形截面柱體尾流進(jìn)行測(cè)量，并運(yùn)用相平均方法對(duì)其尾流場進(jìn)行研究。相關(guān)實(shí)驗(yàn)是在一低湍流度風(fēng)洞中，雷諾數(shù)Red為3 500的條件下進(jìn)行的。研究結(jié)果表明：壁面邊界層與柱體自由端面后的下降流造成了柱體尾流顯著的三維性；該三維性在柱體下游 10d(d為柱體的寬度)處仍能清晰的觀察到；對(duì)于高寬比為3的柱體，相平均的脈動(dòng)速度以及脈動(dòng)速度能譜在所有測(cè)量高度上均近似相同；而高寬比為5或7時(shí)，尾流中間高度附近流動(dòng)周期性最強(qiáng)，且相平均脈動(dòng)速度最大值為柱體頂部或底部附近對(duì)應(yīng)值的2倍以上。

尾流；渦街；熱線風(fēng)速儀；相平均

柱體的繞流和尾流問題普遍存在于工程中，如高層建筑物、水下結(jié)構(gòu)等受到的由流動(dòng)所誘發(fā)的力。針對(duì)二維鈍體尾流，已開展了大量的研究[1?5]。而工程中所遇到的往往是一端固定于壁面，另一端為自由端的有限長度鈍體的尾流。在平面邊界層以及鈍體自由端的影響下，其尾流與二維鈍體尾流有顯著區(qū)別[6?9]。Kawamura等[10]提出了有限長圓柱體的尾流模型，認(rèn)為當(dāng)柱體長徑比超過某一臨界值時(shí)，尾流大部分區(qū)域仍由 Kármán渦街控制，而尾流的上半部分則由一對(duì)產(chǎn)生于柱體自由端的沿流向的頂部渦所控制；當(dāng)長徑比小于該臨界值時(shí)，Kármán渦街消失，整個(gè)尾流由該頂部渦控制。Sakamoto等[11]研究有限長度鈍體尾流的斯坦頓數(shù)St與雷諾數(shù)Re的關(guān)系，發(fā)現(xiàn)當(dāng)長徑比為2～3時(shí)，St?Re曲線的斜率發(fā)生改變，并認(rèn)為此時(shí)在邊界層和頂部渦的作用下尾流結(jié)構(gòu)發(fā)生了改變。從文獻(xiàn)[6?7,11]可知，當(dāng)鈍體為有限長度時(shí)，平面邊界層與產(chǎn)生于自由端的頂部渦會(huì)影響柱體側(cè)面所脫落渦(即 Kármán渦街)的產(chǎn)生，使其明顯區(qū)別于二維鈍體的尾流。當(dāng)鈍體高寬比變化時(shí)，邊界層與頂部渦的相對(duì)作用范圍會(huì)有所變化，因此，尾流特性也會(huì)明顯受高寬比的影響[6?7]。為系統(tǒng)比較高寬比對(duì)正方形截面有限長度柱體尾流的影響，本文作者運(yùn)用了X形熱線對(duì)高寬比為3，5和7的柱體尾流進(jìn)行測(cè)量。運(yùn)用相平均方法對(duì)測(cè)量數(shù)據(jù)進(jìn)行處理，得到不同高度上時(shí)間平均速度、雷諾應(yīng)力以及相平均雷諾應(yīng)力等統(tǒng)計(jì)參數(shù)。

1 實(shí)驗(yàn)方法

實(shí)驗(yàn)在1個(gè)開式低湍流度風(fēng)洞中進(jìn)行，該風(fēng)洞本底湍流度為0.7%，工作段尺寸(長×寬×高)為1.2 m×0.3 m×0.3 m。一正方形截面柱體固定于1個(gè)水平平板上，柱體寬度d＝10 mm。實(shí)驗(yàn)裝置和坐標(biāo)系的定義如圖1所示。實(shí)驗(yàn)中，自由來流速度U∞＝6 m/s，對(duì)應(yīng)的基于柱體寬度的雷諾數(shù)Red＝3 500。水平板距風(fēng)洞底面約為0.1 m，平板前邊緣經(jīng)光滑處理以防止流動(dòng)分離。柱體中心距離水平板前邊緣距離為0.33 m。在柱體軸線位置處，水平板上的邊界層厚度約為 1 mm。實(shí)驗(yàn)中，分別測(cè)量了高寬比(H/d)為3，5和7的柱體尾流。X形熱線與相應(yīng)的參考熱線布置在柱體軸線下游10d的位置。對(duì)每一柱體，分別在不同高度上進(jìn)行測(cè)量?？紤]到流場的對(duì)稱性，實(shí)驗(yàn)僅在y≥0的范圍內(nèi)進(jìn)行。

圖1 實(shí)驗(yàn)裝置及坐標(biāo)系示意圖Fig.1 Experimental setup and coordinate system

為了對(duì)測(cè)量數(shù)據(jù)進(jìn)行相平均，在X熱線的另一側(cè)布置了1個(gè)單絲熱線以提供相位參考信號(hào)。參考熱線位于x*=10，y*≈4。本文中上標(biāo)“*”表示用U∞和d進(jìn)行無量綱化。在所有測(cè)量中，參考熱線高度與X熱線保持一致。參考熱線在流動(dòng)方向上稍為傾斜，如圖 1所示，以得到最好的周期性信號(hào)。實(shí)驗(yàn)中，每個(gè)通道的采樣頻率為3 kHz，每個(gè)測(cè)量點(diǎn)上的采樣時(shí)間為20 s。以二維方形柱體的St數(shù)可估算本實(shí)驗(yàn)中渦團(tuán)脫落頻率約為70 Hz。實(shí)驗(yàn)中熱線采樣頻率遠(yuǎn)高于該值。

圖2 柱體中間高度上v的瞬時(shí)值與過濾后的值Fig.2 Instantaneous and filtered v signal at mid-span

圖2所示為不同高寬比柱體中間高度上，y*=2處速度v以及其過濾后的值。圖2所示為各瞬時(shí)的相位關(guān)系(不同高寬比柱體對(duì)應(yīng)縱坐標(biāo)比例并不相同)。從圖2可以看出：當(dāng)H/d為5和7時(shí)，原始信號(hào)與過濾后的信號(hào)都具有顯著的周期性；當(dāng)H/d為3時(shí)，雖然過濾后信號(hào)仍具有周期性，但該周期性在原始信號(hào)中已經(jīng)相對(duì)較弱。由于所有過濾后信號(hào)均具有周期性，由此可確定任意時(shí)刻所對(duì)應(yīng)的相位。以式(1)確定相位參考點(diǎn)：0，π和2π(如圖2所示)，式中下標(biāo)“f ”表示過濾后的值。時(shí)間序列上任意其他點(diǎn)的相位可由式(2)確定。通過上述方法可確定所有測(cè)量數(shù)據(jù)的相位，具有某一相同相位的瞬時(shí)值可進(jìn)行相平均，具體方法參見文獻(xiàn)[12?13]。通過相平均方法，可以把任何瞬時(shí)量進(jìn)行三元分解[14]，如式(3)所示。

式中：為速度、渦量或雷諾應(yīng)力等瞬時(shí)量；為時(shí)間平均值；為周期性擬序結(jié)構(gòu)所帶來的脈動(dòng)部分；qr為非周期性運(yùn)動(dòng)帶來的隨機(jī)脈動(dòng)部分。

2 實(shí)驗(yàn)結(jié)果與分析

2.1 尾流周期性

圖3 尾流中不同高度上速度v的能譜Fig.3 Power spectral density function of v at various z*

圖3所示是H/d為3和5時(shí)，尾流中不同高度上y*=2處速度v的能譜。H/d=7的能譜與H/d=5時(shí)的能譜非常類似，故在此未給出。由圖 3(d)～(f)可知：當(dāng)H/d=5時(shí)，不同z*處能譜都有明顯的峰值，且其所對(duì)應(yīng)的頻率均相同；對(duì)應(yīng)的無量綱頻率即St數(shù)為0.12，略小于二維方形柱體尾流所對(duì)應(yīng)的值 0.13[2]；對(duì)于有限長柱體，自由端后的下降流使鈍體兩側(cè)所脫落的渦相互遠(yuǎn)離，使整個(gè)尾流變寬，從而降低了渦脫落頻率[15]。當(dāng)H/d≥5時(shí)，柱體中間高度附近能譜峰值最強(qiáng)，而靠近兩端時(shí)峰值減弱。這一現(xiàn)象是由于自由端后的下降流和底部平板邊界層的三維性削弱了其附近尾流的周期性。由圖 3(a)～(c)可知：當(dāng)H/d=3時(shí)，在所有高度上能譜都是類似的，且沒有明顯的峰值，這與H/d=5時(shí)的情況有著顯著的不同。表明當(dāng)H/d=3時(shí)，來自自由端的下降流和底部平板上的邊界層已經(jīng)影響了整個(gè)尾流，尾流中的Kármán渦街非常弱。

2.2 相平均脈動(dòng)速度

如圖3所示，H/d=3時(shí)尾流中周期性運(yùn)動(dòng)已很弱。因此可知由擬序結(jié)構(gòu)所帶來的相平均脈動(dòng)也很弱，故在此未給出。圖4所示為H/d=5時(shí)不同高度上的相平均脈動(dòng)速度。與尾流周期性類似，H/d=7時(shí)的相平均脈動(dòng)速度與H/d=5時(shí)的對(duì)應(yīng)值定性上是類似的，故在此也沒有給出。圖4中各圖長寬比例是由y方向的測(cè)量范圍與渦傳播速度以及渦脫落頻率來確定的[12]。從圖 4可知：在所有測(cè)量高度上都大于。在中間高度附近由擬序結(jié)構(gòu)所產(chǎn)生的和最大，說明該處尾流的周期性最強(qiáng)，這與圖3(b)所反映的信息是一致的。對(duì)比圖4(a)和(d)與圖4(c)和(f)可知：在靠近自由端時(shí)，和均弱于靠近底板時(shí)的對(duì)應(yīng)值。可知：在x*=10處，自由端后下降流的影響要大于底部邊界層的影響。值得注意的是，在z*=5時(shí)，的分布出現(xiàn)類似兩列渦街的分布，這應(yīng)該是自由端后下降流所帶來的速度較高的流體與由柱體兩側(cè)所脫落的準(zhǔn)二維的Kármán渦街相互作用的結(jié)果[15]。比較H/d分別為3, 5和7的結(jié)果可知：當(dāng)H/d≥5時(shí)，柱體中間高度處尾流具有顯著的周期性，而靠近兩端時(shí)相對(duì)較弱；而當(dāng)H/d≤3時(shí)，整個(gè)尾流場的周期性均較弱。Okajima[2]對(duì)不同長徑比的圓柱尾流進(jìn)行了研究，發(fā)現(xiàn)當(dāng)長徑比為2～3時(shí)，柱體尾流中的渦脫落頻率會(huì)發(fā)生突變。他認(rèn)為這是由于當(dāng)長徑比小于這一臨界值時(shí)流場的性質(zhì)發(fā)生了變化，即尾流從準(zhǔn)二維的 Kármán渦街結(jié)構(gòu)變成由下降流所控制的對(duì)稱的尾流結(jié)構(gòu)，周期性也變?nèi)?。這與實(shí)驗(yàn)觀測(cè)結(jié)果類似。

2.3 時(shí)間平均速度及雷諾應(yīng)力

圖4 H/d=5時(shí)的相平均脈動(dòng)速度Fig.4 Phase-average fluctuation velocity for H/d=5

圖5所示為不同高寬比柱體尾流中的沿流向的時(shí)間平均速度在不同高度上的分布情況。當(dāng)H/d為3, 5和7時(shí)，在= 1的高度上，的分布定性上都是類似的。這說明靠近柱體根部的流動(dòng)主要被邊界層所控制，而與柱體的高寬比沒有明顯的關(guān)系。在距離自由端為1d的高度上，3個(gè)高寬比柱體的尾流中均出現(xiàn)類似的分布，即在尾流中心附近(y*≈0)，出現(xiàn)了明顯的峰值，如圖5所示。這與二維鈍體尾流中對(duì)應(yīng)的速度分布有著顯著的區(qū)別。造成這種自由端附近呈現(xiàn)W形分布的原因是有限長鈍體自由端后形成的沿流向的頂部渦會(huì)將尾流區(qū)域外具有較高速度的流體帶入尾流中。這一規(guī)律與Wang等[6,9]所給出的有限長柱體尾流模型是吻合的。值得注意的是：在圖 5(c)中，的分布在=3和4處幾乎完全重合。由此可知：在H/d=7時(shí)，柱體中心高度附近的流動(dòng)已具有一定的二維性。

圖6所示為H/d為3和5時(shí)，尾流中時(shí)間平均的雷諾應(yīng)力和以及對(duì)應(yīng)的由擬序結(jié)構(gòu)所帶來的雷諾應(yīng)力，即和。如前所述，H/d=5和7時(shí)的結(jié)果是類似的，故H/d=7的結(jié)果在此沒有給出。由圖6可知：H/d=5時(shí)雷諾應(yīng)力遠(yuǎn)大于H/d=3的對(duì)應(yīng)結(jié)果，可見尾流中的周期性擬序結(jié)構(gòu)有增強(qiáng)雷諾應(yīng)力的趨勢(shì)，且對(duì)v2更為顯著。對(duì)于3個(gè)高寬比的柱體，在自由端高度上的雷諾應(yīng)力均小于其他高度上的對(duì)應(yīng)值。這是因?yàn)樽杂啥撕蟮南陆盗鲿?huì)將自由流中的流體引入尾流中，而這部分自由流具有較高的速度(如圖5所示)和較低的脈動(dòng)速度。如圖6(a)和(b)所示，在柱體自由端的高度上，的峰值出現(xiàn)在y*≈2附近，這與二維方形柱體尾流的情況完全不同。在二維方柱尾流中，當(dāng)x*＞4時(shí)，的峰值出現(xiàn)在尾流中心。這一現(xiàn)象同樣是下降流帶入尾流中的低湍流度的流體所造成的。由圖 6(a)和(b)可知：擬序結(jié)構(gòu)對(duì)時(shí)均雷諾應(yīng)力的貢獻(xiàn)非常小，且這一現(xiàn)象對(duì)于高寬比較小的柱體更為顯著。對(duì)于二維鈍體尾流，在x*≥40時(shí)，擬序結(jié)構(gòu)對(duì)雷諾應(yīng)力的貢獻(xiàn)近似為0[12]。而在本實(shí)驗(yàn)中，當(dāng)x*=10時(shí)，這一貢獻(xiàn)非常小。這說明在有限長度三維鈍體尾流中，擬序結(jié)構(gòu)在x*=10處已非常弱。

圖5 尾流中時(shí)間平均的流向速度的分布Fig.5 Lateral distribution of in wake

圖6 尾流中雷諾應(yīng)力的分布Fig.6 Distribution of Reynolds stresses in wake

當(dāng)H/d=5時(shí)，在中間高度上和明顯大于其他高度上的對(duì)應(yīng)值，這也驗(yàn)證了在此高度上擬序結(jié)構(gòu)是最強(qiáng)的。而當(dāng)H/d=3時(shí)，和在=1時(shí)最大。這可能是因?yàn)楫?dāng)z*=2時(shí)已經(jīng)受下降流的影響，抑制了和。擬序結(jié)構(gòu)對(duì)的貢獻(xiàn)遠(yuǎn)大于其對(duì)的貢獻(xiàn)，且這一現(xiàn)象隨柱體高寬比的增加更為顯著。此外，當(dāng)=1時(shí)，和對(duì)3個(gè)高寬比的柱體尾流都是近似的。也就是說，在這個(gè)高度上的流動(dòng)特性主要由平板邊界層所控制，而與柱體高寬比關(guān)系不大。

3 結(jié)論

(1) 對(duì)于有限長度方形截面柱體，當(dāng)高寬比H/d≥5時(shí)，尾流中間高度上的流動(dòng)具有明顯的周期性擬序結(jié)構(gòu)。這一擬序結(jié)構(gòu)隨著靠近柱體的兩端而逐漸減弱；當(dāng)H/d=3時(shí)，由于自由端后下降流以及底部邊界層已影響到整個(gè)尾流范圍，整個(gè)尾流中都未出現(xiàn)明顯的擬序結(jié)構(gòu)。

(2) 自由端后的下降流將一部分具有較高速度且較低湍流度的流體引入了尾流中，造成其與二維鈍體尾流有顯著區(qū)別。在靠近自由端的高度上，尾流中心部分速度較高，而雷諾應(yīng)力較小。這一規(guī)律隨高寬比的增大而更加明顯。

(3) 由于周期性擬序結(jié)構(gòu)較弱，有限長柱體尾流中雷諾應(yīng)力均小于二維鈍體尾流中的對(duì)應(yīng)值。這一規(guī)律隨高寬比的減小而更加顯著。擬序結(jié)構(gòu)對(duì)的貢獻(xiàn)遠(yuǎn)大于其對(duì)的貢獻(xiàn)。尾流中尤其是其中間高度上，擬序結(jié)構(gòu)對(duì)雷諾應(yīng)力的貢獻(xiàn)隨著高寬比的增大顯著增加。

[1]Saha A K, Muralidhar K, Biswas G. Experimental study of flow past a square cylinder at high Reynolds numbers[J]. Experiments in Fluids, 2000, 29(66): 553?563.

[2]Okajima A. Strouhal number of rectangular cylinders[J]. Journal of Fluid Mechanics, 1982, 123: 379?398.

[3]Durao D, Heitor M, Pereira J. Measurements of turbulent and periodic flows around a square cross-section cylinder[J].Experiments in Fluids, 1988, 6(5): 298?304.

[4]Zhou Y, Antonia R A. A study of turbulent vortices in the near wake of a cylinder[J]. Journal of Fluid Mechanics, 1993, 253:643?661.

[5]Zhang N, Li G Z, Huang J C. Numerical study of vortex shedding from a square cylinder and its effects on heat transfer of mixed convection[J]. Journal of Hydrodynamics B, 2003,15(3): 103?111.

[6]Wang H F, Zhou Y, Chan C, et al. Momentum and heat transport in a finite-length cylinder wake[J]. Experiments in Fluids, 2009,46(6): 1173?1185.

[7]Summer D, Heseltine J L, Dansereau O J P. Wake structure of a finite circular cylinder of small aspect ratio[J]. Experiments in Fluids, 2004, 37(5): 720?730.

[8]Adaramola M S, Akinlade O G, Sumner D, et al. Turbulent wake of a finite circular cylinder of small aspect ratio[J]. Journal of Fluids and Structures, 2006, 22(6): 919?928.

[9]Park C W, Lee S J. Free end effects on the near wake flow structure behind a finite circular cylinder[J]. International Journal of Wind Engineering and Industrial Aerodynamics, 2000,88(2): 231?246.

[10]Kawamura T, Hiwada M, Hibino T. Flow around a finite circular cylinder on a flat plate[J]. Bulletin of JSME, 1984, 27(232):2142?2151.

[11]Sakamoto H, Arie M. Vortex shedding from a rectangular prism and a circular cylinder placed vertically in a turbulent boundary layer[J]. Journal of Fluid Mechanics, 1983, 126: 147?165.

[12]Zhou Y, Zhang H J, Yiu M W. The turbulent wake of two side-by-side circular cylinder[J]. Journal of Fluid Mechanics,2002, 458: 303-332.

[13]Hu J C, Zhou Y. Flow structures behind two staggered circular cylinders. Part 2. Heat and momentum transport[J]. Journal of Fluid Mechanics, 2008, 607: 81?107.

[14]Reynolds W C, Hussain A K M F. The mechanism of an organized wake in turbulent shear flow. Part 3. Theoretical model and comparisons with experiments[J]. Journal of Fluid Mechanics, 1972, 54: 263-288.

[15]Zdravkovich M M. Flow around circular cylinder. vol 2:Applications[M]. New York: Oxford University Press, 2002:726?728.

(編輯 楊幼平)

Three dimension characteristics of finite-length bluff body wake

WANG Han-feng

(School of Civil Engineering and Architecture, Central South University, Changsha 410075, China)

The wake flow behind a finite-length cylinder with aspect ratio of 3, 5 and 7 was measured using X-wire.Phase-averaged technique was applied. All experiments were conducted in a low-speed wind-tunnel atRed=3 500. It is found that the boundary layer on the wall and the free-end downwash flow make the near wake highly three dimensional.This three-dimensionality is clearly observed at 10d(dis the cylinder width) downstream from the cylinder. ForH/d=3,the phase-averaged fluctuation velocity and the power spectra are similar in all spanwise measurement locations. AsH/d=5 and 7, the periodicity of the wake is the most pronounced at the mid-span, and the phase-averaged fluctuation velocity is twice bigger than that near the free-end and wall.

wake; vortex street; hot-wire anemometer; phase average

O351.3

A

1672?7207(2011)02?0495?06

2010?02?26；

2010?04?16

國家自然科學(xué)基金資助項(xiàng)目(50208019)；中南大學(xué)自由探索計(jì)劃項(xiàng)目(201012200215)

王漢封(1976?)，男，河南開封人，博士，副教授，從事流動(dòng)的測(cè)量與診斷以及流動(dòng)所誘發(fā)振動(dòng)的研究；電話：13875856504；E-mail:wanghfme@gmail.com