王 勇, 郝南松, 耿子海, 王萬波

(1. 中國空氣動力研究與發(fā)展中心 空氣動力學(xué)國家重點實驗室, 四川 綿陽 621000; 2. 中國空氣動力研究與發(fā)展中心 氣動噪聲控制重點實驗室, 四川 綿陽 621000)

0 引 言

在實際工程應(yīng)用中，如航空中的飛機翼間支柱、海洋開發(fā)中的深海鉆井平臺及其鉆桿、風(fēng)工程中的樁基碼頭和大橋橋墩等，以圓柱繞流為代表的鈍物體繞流是一個廣泛存在的問題。由于圓柱繞流尾跡中常常伴隨著復(fù)雜的流動現(xiàn)象，如邊界層分離、旋渦形成與脫落及耗散等[1-2]，開展圓柱體繞流尾跡特性研究具有重要的實際意義。

對圓柱繞流尾跡的研究，有理論分析、實驗研究和數(shù)值模擬3種方法。由于理論簡化受限、實驗耗費較大，現(xiàn)階段的研究以數(shù)值模擬為主、實驗研究和理論分析為輔。熱線風(fēng)速儀[3-4](Hot-Wire Anemometer, HWA) 和粒子圖像測速[5-6](Particle Image Velocimetry, PIV)是常用于圓柱繞流研究的實驗設(shè)備。熱線風(fēng)速儀價格便宜、頻率響應(yīng)高(達20kHz以上)，但只能進行單點測量且與流場直接接觸，會對繞流尾跡產(chǎn)生顯著影響；另一方面，PIV能夠非接觸地獲取較大范圍內(nèi)整個流場的瞬時速度信息。因此，PIV是對圓柱繞流尾跡進行定性的流場顯示和定量的流場測量與分析的理想工具。

Julio[7]利用采樣頻率fPIV=25Hz的數(shù)字式PIV對雷諾數(shù)Re=875和769的圓柱繞流進行了測量，討論了垂直于圓柱軸線的平面和經(jīng)圓柱軸線沿流向平面上的瞬時速度場和渦量場。Krothapalli等[8]采用fPIV=10Hz的PIV對Re=(0.85～1.5)×105的圓柱近尾流區(qū)域進行了測量，實驗結(jié)果給出了瞬時速度場和瞬時渦量場。Braza等[9]和Perrin等[10]采用fPIV=4Hz的2D2C PIV和2D3C PIV研究了Re=1.4×105的圓柱繞流的湍流特性，并通過相位平均分析展示了繞流尾跡中旋渦交替脫落和發(fā)展的演化過程。Sung和Yoo[11]在水洞中采用fPIV=125Hz的PIV對Re=360的圓柱繞流進行了測量，提出了一種從二維PIV數(shù)據(jù)中獲取三維相位平均數(shù)據(jù)的方法。Konstantinidis等[12]采用PIV在水洞中對Re=2160的自然圓柱繞流和強迫圓柱繞流進行了測量，對比分析了2種情況下的平均速度場、雷諾應(yīng)力和相位平均速度場。張瑋等[13]采用fPIV=10Hz的低頻PIV對雷諾數(shù)Re=500的圓柱繞流尾跡進行了測量，給出了圓柱兩側(cè)各1倍直徑、圓柱后方2.2倍直徑所圍成區(qū)域內(nèi)的一系列瞬態(tài)速度矢量場、渦量場和流線圖，并著重分析了1個擬周期內(nèi)旋渦的發(fā)展演化過程。張孝棣等[14]采用fPIV=8Hz的PIV系統(tǒng)對Re=104～105范圍內(nèi)的圓柱繞流尾跡進行了測量，主要分析了繞流尾跡的瞬時速度場、渦量場和流線圖。涂程旭等[15]采用fPIV=15Hz的PIV對Re=630、800和950等3種雷諾數(shù)的圓柱繞流場進行了實驗研究，獲得了速度場、渦量場以及渦脫落的時空演化規(guī)律。

盡管基于PIV的圓柱繞流研究已經(jīng)取得了很多成果，并為理論分析和數(shù)值模擬提供了基礎(chǔ)數(shù)據(jù)，但現(xiàn)有的研究存在2個方面的不足。首先，PIV的采樣頻率普遍偏低，只有幾十赫茲，而高雷諾數(shù)下圓柱繞流的旋渦脫落頻率往往達成百上千赫茲，導(dǎo)致PIV采樣的奈奎斯特(Nyquist)頻率難以超過旋渦脫落頻率。這就使得現(xiàn)有研究通常只能給出流場的速度場、渦量場和流線圖，而難以對流場進行深入的統(tǒng)計分析(如頻譜分析等)，對圓柱繞流旋渦脫落頻率的分析需要借助于數(shù)值模擬或者升力系數(shù)、阻力系數(shù)等直接與壓力有關(guān)的測量。其次，圓柱繞流的相位平均分析中，需要引入1個額外的參考信號以區(qū)分不同時刻流場的相位，如Sung和Yoo[11]、Wlezien和Way[16]使用選定點或熱線測量的速度信號，Braza等[9]、Perrin等[10]、Cantwell和Coles[17]使用圓柱表面的壓強信號。但是，參考信號和尾流速度信號之間隨機產(chǎn)生的相位抖動(Phase jitter)會使得相位平均結(jié)果高估湍流運動并使卡曼渦的結(jié)構(gòu)更光滑[18]。

本文中使用采樣率達fPIV=1000Hz的時間解析PIV系統(tǒng)對圓柱后方7.5倍直徑、圓柱兩側(cè)各3.3倍直徑所圍成范圍內(nèi)的繞流尾跡在雷諾數(shù)Re=2.74×104下的非定常流場進行了測量。針對PIV獲得的速度場數(shù)據(jù)，進行流場和頻譜特性分析，探討了圓柱繞流尾跡中的渦量場、旋轉(zhuǎn)強度和旋渦脫落的頻率特性等。同時，提出了一種簡單有效的基于速度場之間相關(guān)性的相位平均分析方法，系統(tǒng)分析了圓柱上下兩側(cè)旋渦交替生成、脫落、發(fā)展并耗散的完整演化過程。

盡管關(guān)于圓柱尾流研究已經(jīng)有比較豐富的研究成果，但本文在實驗測量、數(shù)據(jù)處理與分析方面有一些新的創(chuàng)新點：(1) 借助時間解析PIV高時間分辨率和高空間分辨率的優(yōu)點，對圓柱繞流尾跡進行了精細(xì)測量。高時間分辨率測量使得可以直接對PIV速度場數(shù)據(jù)進行頻譜分析而得到旋渦脫落的頻率特性。另一方面，全流場、高空間分辨率測量使得我們可以得到一個較大平面上的非定常流場的頻率特性(如斯特勞哈爾數(shù)云圖)，而不只是某幾個點的頻率特性。(2) 本文提出了基于速度相關(guān)性的相位平均方法來處理非定常結(jié)果，該方法簡單有效，能有效地捕捉圓柱尾流中旋渦交替生成、脫落、發(fā)展并耗散的時空演化過程，在非定常流場測量(如柱翼組合模型和空腔模型的流場)方面具有普遍推廣意義。

1 實驗裝置及方法

實驗在中國空氣動力研究與發(fā)展中心的0.55m×0.4m聲學(xué)風(fēng)洞開口試驗段中進行。該風(fēng)洞是一座單回流式低速風(fēng)洞，具有開/閉口2個試驗段，其中開口試驗段長1.5m、寬0.55m、高0.4m，最大風(fēng)速達100m/s。

速度場測量采用TSI公司的時間解析PIV系統(tǒng)，照明光源為Photonics公司的DM30單腔激光器(能量2×15mJ@1000Hz、雙脈沖工作方式)，通過片光組件以水平方向照射圓柱軸向中央截面流場，片光厚度約1mm。使用Photran公司的SA4相機(配備Nikon 50/1.2鏡頭，分辨率1024pixel×1024pixel，最大滿幀率3600幀/s)以垂直片光方向?qū)ζ庹彰鲄^(qū)域進行拍攝并采集PIV圖像。實驗中設(shè)置粒子圖像對的采樣頻率為1000Hz，2束激光的時間間隔為40μs。激光與相機之間的同步控制采用BNC公司的555時序控制器。PIV試驗過程中采用DF-1500型煙霧發(fā)生器產(chǎn)生DEHS油霧粒子對流場進行示蹤。

對PIV采集的粒子圖像采用德國PIVTEC公司開發(fā)的PIVview2C軟件對任意相鄰的粒子圖像對做互相關(guān)分析以獲取所拍攝平面區(qū)域的速度場。處理過程中采用標(biāo)準(zhǔn)FFT相關(guān)算法、Nyquist頻率濾波和9點最小二乘高斯擬合峰值算法等確定相關(guān)性，采用多次相關(guān)(Multiple correlation)和多網(wǎng)格查詢(Multi-grid interrogation)方法提高分析精度，初始查詢窗為32×32、最終查詢窗為16×16，重疊因子50%。

實驗使用的實心鋁制圓柱模型(展長L=550mm，直徑D=20mm)一端固定連接于支撐系統(tǒng)的頂部連接板，另一端用鐵絲拉緊以消除模型振動。實驗中，圓柱表面噴涂黑色啞光漆以消除強激光在模型表面反光而對PIV圖像產(chǎn)生影響。整個實驗的設(shè)置示意圖如圖1所示。

實驗中，自由來流風(fēng)速U∞=20m/s，基于圓柱直徑的雷諾數(shù)為Re=ρU∞D(zhuǎn)/μ=2.74×104，此時圓柱繞流處于亞臨界狀態(tài)，邊界層為層流分離而尾流已轉(zhuǎn)變?yōu)橥牧鳒u街。PIV采樣時間為2.7s，共獲取2700幅粒子圖像對。

圖1 實驗設(shè)置示意圖(俯視圖)

(1)

其中，cov(Xk,Xl)為Xk和Xl的協(xié)方差，var(Xk)和var(Xl)分別為Xk和Xl的方差。

2 實驗結(jié)果及分析

2.1 平均流場特性分析

圖2給出了根據(jù)PIV測量結(jié)果計算得到的平均速度場的流線圖，從圖中可以看出，在圓柱壁面正后方約1.4D的區(qū)域內(nèi)存在2個與圓柱半徑同尺度的反向旋轉(zhuǎn)的旋渦，旋渦中心位于0.8D的位置附近，圓柱上下兩側(cè)交替脫落的旋渦進入尾流場后形成了低速回流區(qū)。

圖2 平均速度場的流線圖

圖3給出了平均速度場的歸一化渦量圖以及根據(jù)λci[19]準(zhǔn)則確定的旋轉(zhuǎn)強度(Swirling strength)圖。從圖中可以看出：(1) 渦量模值較大的區(qū)域主要位于圓柱的上、下兩側(cè)，“帶狀”分布的渦量沿圓柱后方回流區(qū)的邊緣向下游擴散、向旋渦區(qū)聚集并逐漸減弱。(2) 圓柱上、下兩側(cè)的渦量強度相同但符號相反，順時針旋轉(zhuǎn)的旋渦對應(yīng)于負(fù)的渦量，而逆時針旋轉(zhuǎn)的旋渦對應(yīng)于正的渦量。(3) 旋轉(zhuǎn)強度沿著圓柱后方回流區(qū)的邊緣逐漸增強，在0.8D的位置附近旋轉(zhuǎn)強度達最大值，隨后旋轉(zhuǎn)強度迅速減弱。(4) 渦量和旋轉(zhuǎn)強度均表明，回流中心區(qū)域(0.8D的位置附近)是流動結(jié)構(gòu)變化最劇烈的區(qū)域。

圖3 平均速度場的歸一化渦量圖和旋轉(zhuǎn)強度圖

2.2 脈動流場特性分析

(2)

圖4 脈動速度場的湍動能和歸一化雷諾剪切應(yīng)力

2.3 頻譜特性分析

由于實驗中采用的時間解析PIV采樣頻率達1000Hz，大于圓柱尾渦脫落頻率的2倍，因此可以對PIV速度場數(shù)據(jù)進行頻譜分析。采用pWelch方法計算每個PIV網(wǎng)格點上速度場的功率譜密度(Power Spectra Density, PSD)。pWelch方法中引入Hamming窗以減小信號分析的“頻譜泄露”，分塊大小為1024，重疊因子為50%，頻率分辨率為0.98Hz。

圖5給出了圓柱繞流速度場的斯特勞哈爾數(shù)St=fD/U∞云圖，其中f為功率譜密度的最大值對應(yīng)的頻率，即旋渦脫落頻率。從圖中可以看出，除了視場的上下邊緣和前緣外，圓柱繞流尾跡中的斯特勞哈爾數(shù)穩(wěn)定在0.2左右，這與文獻[20-26]中報道的亞臨界雷諾數(shù)條件下圓柱繞流的St≈0.20～0.21是一致的。

圖5 斯特勞哈爾數(shù)云圖

2.4 相位平均分析

圖7和8分別給出了根據(jù)速度場之間的相關(guān)性確定的相位平均流線圖和渦量圖，從圖中可以看到：(1) 圓柱上、下兩側(cè)不斷地有旋渦交替地生成，從圓柱表面脫落向下游發(fā)展并隨著能量的衰竭而耗散，形成較規(guī)則的Karman渦街脫落。(2) 在t=0T(T為圓柱尾渦脫落的周期)時刻，圓柱正后方約1.3D處存在1個順時針旋轉(zhuǎn)的大尺度旋渦，在旋渦附近存在負(fù)的渦量聚集區(qū)。在緊鄰圓柱表面的右下方，1個逆時針旋轉(zhuǎn)的小旋渦開始形成。(3) 在t=1/5T時刻，順時針旋轉(zhuǎn)的大尺度旋渦沿流向明顯拉伸并向四周擴散，其渦量降低并開始衰竭。圓柱下方的小旋渦進一步長大、渦量逐漸聚集并增大到最大值。(4) 在t=2/5T時刻，圓柱右下方逆時針旋轉(zhuǎn)的旋渦已經(jīng)發(fā)展得非常飽滿，成為大尺度旋渦，其中心位于圓柱正后方約1.0D的位置，正的渦量聚集區(qū)進一步發(fā)展。圓柱右上方t=0T時刻生成的大尺度旋渦已經(jīng)擴散并向下游移動，渦量進一步降低、聚集區(qū)逐漸消失。此時，與圓柱表面相鄰的右上方，1個新的旋渦結(jié)構(gòu)開始醞釀形成。(5) 在t=3/5T時刻，逆時針旋轉(zhuǎn)的旋渦隨著向下游流動，旋渦中心后移、渦量降低，開始進入被主流拉伸的初始階段。圓柱右上方形成的小尺度旋渦繼續(xù)長大。(6) 在t=4/5T時刻，圓柱上方順時針旋轉(zhuǎn)的旋渦已經(jīng)發(fā)展得比較飽滿、其中心位于圓柱正后方約0.9D的位置，負(fù)的渦量聚集區(qū)進一步擴大。t=2/5T時刻發(fā)展得非常飽滿的大尺度旋渦已經(jīng)擴散并向下游移動，一個新的逆時針旋轉(zhuǎn)的旋渦結(jié)構(gòu)開始在緊鄰圓柱的右下方醞釀形成。(7)t=T時刻的圓柱尾流特征，不論是流線還是渦量分布，都與t=0T時刻的尾流特征完全一致，圓柱繞流的旋渦演化過程完成1個完整的周期。

圖6 速度場之間的相關(guān)性分析

圖7 相位平均流線圖

圖8 相位平均渦量場

3 結(jié) 論

本文使用采樣頻率為fPIV=1000Hz的時間解析PIV系統(tǒng)對圓柱后方7.5倍直徑、圓柱兩側(cè)各3.3倍直徑所圍成范圍內(nèi)的繞流尾跡在雷諾數(shù)Re=2.74×104下的非定常流場進行了細(xì)致地測量。PIV采樣頻率高于旋渦的變化頻率，一方面使得可以直接利用速度場數(shù)據(jù)分析圓柱繞流尾跡中的旋渦脫落頻率，另一方面由于每個渦脫落周期能捕獲到若干幅PIV圖像從而使得基于速度場之間相關(guān)性的相位平均分析成為可能。實驗研究結(jié)果表明：

(1) 在圓柱后方1.4D的范圍內(nèi)存在1個低速回流區(qū)，其中心0.8D位置附近是流動結(jié)構(gòu)變化最劇烈的區(qū)域；

(2) 圓柱后方1.9D位置附近是上/下兩側(cè)脫落旋渦交匯、耦合的區(qū)域，湍流強度和雷諾剪切應(yīng)力較大，湍流脈動最強；

(3) 每個PIV網(wǎng)格點上速度場的頻譜分析表明，圓柱繞流尾跡中的旋渦脫落頻率對應(yīng)的斯特勞哈爾數(shù)穩(wěn)定在0.2左右；

(4) PIV速度場之間相關(guān)性系數(shù)的頻譜分析表明，功率譜密度的峰值頻率與圓柱后方的旋渦脫落頻率是完全一致的?；谒俣葓鲋g相關(guān)性的相位平均分析方法簡單有效，可以準(zhǔn)確地識別繞流尾跡中旋渦交替生成、脫落、發(fā)展并耗散的時空演化過程，并且該方法可以很容易地推廣應(yīng)用到其它非定常流場的測量和分析中。

[1]Von Karman T H, Rubach H. On the mechanism of resistance in fluid[J]. Physikalische Zeitschrift, 1912, 13(5): 351-358.

[2]Williamson C H K. Vortex shedding in the cylinder wake[J]. Annu Rev Fluid Mech, 1996, 28: 477-539.

[3]Raffel M, Willert C, Wereley S, et al. Particle image velocimetry: a practical guide[M]. 2nd ed. Berlin Heidelberg: Springer-Verlag, 2007.

[4]Adrian R J. Twenty years of particle image velocimetry[J]. Exp Fluids, 2005, 39(2): 159-169.

[5]Comte-Bellot G. Hot-wire anemometry[J]. Annu Rev Fluid Mech, 1976, 8: 209-231.

[6]Perry A E. Hot-wire anemometry[M]. Oxford: Clarendon Press, 1982.

[7]Julio S. An investigation of the near wake of a circular cylinder using a video-based digital cross-correlation particle image velocimetry technique[J]. Experimental Thermal and Fluid Science, 1996, 12(2): 221-233.

[8]Krothapalli A, Shih C, Lourenco L. The near wake of a circular cylinder at 0.3

[9]Braza M, Perrin R, Hoarau Y. Turbulence properties in the cylinder wake at high Reynolds numbers[J]. Journal of Fluids and Structures, 2006, (22): 757-771.

[10]Perrin R, Cid E, Cazin S, et al. Phase-averaged measurements of the turbulence properties in the near wake of a circular cylinder at high Reynolds number by 2C-PIV and 3C-PIV[J]. Exp Fluids, 2007, (42): 93-109.

[11]Sung J, Yoo J Y. Three-dimensional phase averaging of time-resolved PIV measurement data[J]. Meas Sci Technol, 2001, (12): 655-662.

[12]Konstantinidis E, Balabani S, Yianneskis M. Conditional averaging of PIV plane wake data using a cross-correlation approach[J]. Exp Fluids, 2005, (39): 38-47.

[13]張瑋, 王元, 徐忠, 等. 圓柱繞流渦系演變的DPIV測試[J]. 空氣動力學(xué)學(xué)報, 2002, 20(4): 379-387.

Zhang W, Wang Y, Xu Z, et al. Experimental investigation of vortex evlovement around a circular cylinder by digital particle image velocimetry (DPIV)[J]. Acta Aerodynamica Sinca, 2002, 20(4): 379-387.

[14]張孝棣, 蔣甲利, 賈元勝, 等. 圓柱體繞流尾跡的PIV測量[J]. 實驗流體力學(xué), 2005, 19(2): 74-78.

Zhang X D, Jiang J L, Jia Y S, et al. Measurement of cylinder’s wake by PIV[J]. Journal of Experiments in Fluid Mechanics, 2005, 19(2): 74-78.

[15]涂程旭, 王昊利, 林建忠. 圓柱繞流的流場特性及渦脫落規(guī)律研究[J]. 中國計量學(xué)院學(xué)報, 2008, 19(2): 98-102.

Tu C X, Wang H L, Lin J Z. Experimental research on the flow characteristics and vortex shedding in the flow around a circular cylinder[J]. Journal of China Jiliang University, 2008, 19(2): 98-102.

[16]Wlezien R W, Way J L. Techniques for the experimental investigation of the near wake of a circular cylinder[J]. AIAA J, 1979, 17(6): 563-570.

[17]Cantwell B, Coles D. An experimental study of entrainment and transport in the turbulent near wake of a circular cylinder[J]. J Fluid Mech, 1983, 136: 321-374.

[18] Perrin R, Braza M, Cid E, et al. Coherent and turbulent process analysis in the flow past a circular cylinder at high Reynolds number[J]. Journal of Fluids and Structures, 2008, 24(8): 1313-1325.

[19]Zhou J, Adrian R J, Balachandar S, et al. Mechanisms for generating coherent packets of hairpin vortices in channel flow[J]. J Fluid Mech, 1999, 387: 353-396.

[20]Kravchenko A G, Moin P. Numerical studies of flow over a circular cylinder atRe=3900[J]. Physics of Fluids, 2000, 12(2): 403-417.

[21] Parnaudeau P, Carlier J, Heitz D, et al. Experimental and numerical studies of the flow over a circular cylinder at Reynolds number 3900[J]. Physics of Fluids, 2008, 20(8): 261-271.

[22] Schlichting H. Boundary layer theory[M]. New York: McFraw-hill Book Company, 1968.

[23] Franke J, Frank W. Large eddy simulation of the flow past a circular cylinder atRe=3900[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(10): 1191-1206.

[24]Breuer M. A challenging test case for large eddy simulation: high Reynolds number circular cylinder flow[J]. Int J Heat Fluid Flow, 2000, 21: 648-654.

[25] Dong S, Karniadakis E, Ekmekci A, et al. A combined direct numberical simulation-partical image velocimetry study of the turbulent near wake[J]. Journal of Fluid Mechanics, 2006, 569(12): 185-207.

[26]Fey U, Konig M, Eckelmann H. A new Strouhal-Reynolds number relationship for the circular cylinder in the range 47