魏 崢, 夏 超, 袁海東, 李啟良, 楊志剛,*

(1.上海市地面交通工具空氣動力與熱環(huán)境模擬重點實驗室, 上海 201804;2. 同濟大學 上海地面交通工具風洞中心, 上海 201804)

?

覆蓋多孔介質的圓柱尾跡實驗研究

魏 崢1,2, 夏 超1,2, 袁海東1,2, 李啟良1,2, 楊志剛1,2,*

(1.上海市地面交通工具空氣動力與熱環(huán)境模擬重點實驗室, 上海 201804;2. 同濟大學 上海地面交通工具風洞中心, 上海 201804)

利用風洞試驗研究了多孔介質對圓柱繞流的影響，通過熱線、煙線和粒子圖像測速對有/無覆蓋多孔介質柱體尾跡的頻率特征、平均流場和瞬時流場進行了測量和分析，進而探究多孔介質與流場之間的作用機理。結果表明：引入多孔介質后，尾跡的大尺度渦脫落結構得到有效抑制，圓柱渦脫落頻率降低；多孔介質減弱了近尾跡場的速度功率譜密度峰值，增大了遠場的速度功率譜密度峰值，渦脫落位置延后；多孔介質減小了尾跡區(qū)的速度，拓寬了尾跡的寬度，使得兩側剪切層變得細長，削弱了它們之間的相互作用，減緩了兩側渦相互作用的速率；同時，多孔介質減小了尾跡區(qū)的速度脈動以及雷諾應力，雷諾應力和平均渦量的峰值更加遠離圓柱，減弱了剪切層的不穩(wěn)定性，增大了渦形成長度。

圓柱；多孔介質；尾跡；粒子圖像測速

0 引 言

圓柱，作為典型的鈍體，在飛機起落架、高速列車受電弓、風力發(fā)電桅桿等，有著廣泛的應用。圓柱尾流區(qū)渦流的產生和變化對阻力、振動和噪聲有重要影響。因此，研究圓柱周圍流場及其控制方法具有重要的工程意義。

通常來說，按照是否需要外界提供能量，流動控制方法可分為主動控制和被動控制。主動控制需要提供復雜的機械設備，例如等離子體發(fā)生器，聲激勵等[1]。被動控制不需要能量，簡單、有效、更易實施，例如柔性壁[2]、改變表面形狀等[3]。另一種可行的措施是在鈍體表面覆蓋多孔介質[4-5]。多孔介質改變鈍體周圍流場的特性，如表面壓力分布、氣動力脈動、速度場、渦脫落頻率和尾跡結構。Sueki等[6-7]對包裹多孔介質的圓柱和高速列車受電弓進行風洞試驗測量其產生的噪聲，明確了多孔介質的降噪效果。Ruck等[8]通過阻力測量試驗發(fā)現(xiàn)，在圓柱背風側覆蓋多孔介質可以降低阻力達8%，而后又對其進行壓力測量得出造成上述結果的原因是背風側覆蓋的多孔介質使該處相比于光滑圓柱體壓力系數(shù)有所增加；同時利用PIV分析了流場，表明背風側覆蓋多孔介質的圓柱其分離渦位置延后，湍動能和速度脈動減小[9]。Liu等[10-11]對二維圓柱進行數(shù)值計算，結果表明多孔介質減小了氣動力脈動的幅值和頻率，降低渦脫落頻率，使得尾跡結構更加穩(wěn)定，具有降噪的效果。

目前國內外學者已對覆蓋多孔介質的圓柱的尾跡進行了一定的研究，并且明確了其對降噪的貢獻以及在實際工程中的作用。但是，在分析其對流場的影響機理上大多基于二維仿真結果，試驗方面還有所欠缺。基于此，本文通過一系列的風洞試驗研究圓柱覆蓋多孔介質前后尾跡流場的變化，并對多孔介質與流場之間的作用機理進行進一步的探究和分析。

本文主要利用熱線試驗測量圓柱尾跡設定測點的瞬時流向速度以獲得其頻率特征；煙線試驗捕捉流動細節(jié)，呈現(xiàn)瞬時渦結構；粒子圖像測速(PIV)技術測量圓柱水平面上的瞬時流場，進而可得到平均流場、速度脈動、渦形成長度和雷諾應力等。針對有/無覆蓋多孔介質的圓柱的流場結果對比分析，探究多孔介質對圓柱尾跡的控制機理。

1 多孔介質和試驗柱體

多孔材料根據結構可分為開孔式和閉孔式[12]。閉孔式各空隙之間相互獨立，對流場影響效果?。欢S蜂窩狀的開孔式結構對流場有明顯的干擾作用[7]。本文研究對象是覆蓋開孔式介質的圓柱，并與光滑圓柱進行對比。

試驗圓柱體長400 mm，外徑均為50 mm，多孔介質層厚10 mm，如圖1所示。代表多孔介質孔洞體積占比的孔隙率為97%；每英寸孔洞的數(shù)量，即PPI為13。

圖1 試驗柱體示意圖Fig.1 Outline of cylinder specimens

2 試驗設置和方法

試驗在開口回流式風洞內進行，噴口寬0.43 m、高0.28 m，試驗段長1.0 m。噴口最大風速45 m/s，湍流強度小于0.5%。圓柱垂直固定在中線上，距離噴口0.3 m。圓柱阻塞比為11.6%。

2.1 熱線試驗

利用兩個沿展向對稱布置的一維探針(Dantec 55P01)同時測量尾跡的瞬時流向速度。測點位置如圖2所示。探針連接在恒溫熱線風速儀上(Dantec streamline)，采樣頻率為2 kHz，采樣時間為16 s。試驗過程中，風洞測試段內的溫度變化小于0.5 ℃。試驗風速為30 m/s，其對應的雷諾數(shù)Re=1.0×105。

圖2 熱線探測點Fig.2 Schematic of hot-wire probe positions

2.2 煙線試驗

煙線流動顯示試驗通過加熱金屬絲上的油劑，產生清晰的煙絲，定性地分析圓柱尾部渦結構的變化。兩根金屬絲平行放置于柱體前后方，如圖3所示。100 V可調的直流電源用于加熱金屬絲。油劑采用潤滑油，其產生的煙霧被閃光燈照亮，由CANON 70D相機捕捉。閃光燈和拍攝由單片機控制協(xié)調。試驗風速為2.4 m/s，對應的雷諾數(shù)為8.0×103。

圖3 煙線試驗示意圖Fig.3 Schematic of smoke-wire visualization

2.3 PIV測量

本試驗采用的是TSI PowerView的PIV系統(tǒng)，主要包括CCD相機、激光器、控制器和同步器。Vlite-500脈沖激光器發(fā)射激光波長為532 nm，最大輸出能量可達到500 mJ/脈沖。CCD相機(PowerViewPlus 29M-HS)像素6600×4400，垂直于激光平面安裝。同步器控制片光的照射和圖片的捕捉，采樣頻率是1.5 Hz，樣本數(shù)量是500張。來流速度為30 m/s，對應的雷諾數(shù)是1.0×105。圖4為PIV試驗設置示意圖，所測平面為水平中截面，地面邊界層影響和端部效應可忽略。

圖4 PIV試驗設置Fig.4 Schematic of the PIV setup

3 結果討論

3.1 渦脫落頻率

圖5顯示了在圓柱一側y/D=-1.5、不同x/D處流向速度的功率譜密度(PSD)分布?？梢钥吹剑涸诓煌牧飨蛭恢?，覆蓋多孔介質的柱體St為0.162，小于光柱的St值(0.182)。當x/D≤1.5時，覆蓋多孔介質柱體的PSD峰值得到了抑制，尤其在x/D=1處，其峰值明顯小于光柱。而當x/D>1.5后，覆蓋多孔介質柱體的PSD峰值逐漸超過光柱。圖6和圖7展示了兩個展向對稱布置的探針測得的瞬時流向速度信號的相位特征，相位差接近±π表明渦從柱體兩側交替脫落。從圖可以觀察到：覆蓋多孔介質柱體的反相位出現(xiàn)的位置延后，而且其反相位的范圍縮小。綜上可知：多孔介質降低了柱體尾跡的渦脫落頻率，抑制了近尾跡區(qū)域的PSD峰值，使渦脫落位置延后。

圖5 y/D=-1.5處流向速度的PSD分布Fig.5 PSD distributions of U velocity at y/D=-1.5

圖7 多孔柱體尾跡兩個瞬時熱線信號間的相位譜Fig.7 Spectral phase between two hot-wire signals simultaneously measured in the wake of the cylinder with porous layer coating

3.2 平均流場

圖8呈現(xiàn)了覆蓋多孔介質前后圓柱尾跡的平均流場變化。相比較而言，多孔柱體尾跡區(qū)域速度明顯減小，尾跡被拓寬，并沿主流方向被拉長。通過圖9可以定量地比較不同位置流向平均速度的變化情況。如圖所示，覆蓋多孔介質后，柱體尾跡流向的平均速度減小，尤其在中線附近，但剪切層內速度梯度更大，多孔柱體尾流的低速區(qū)被拓寬。進一步觀察可以看出，平均流場的回流區(qū)(Umean<0)已從x/D=3延長至x/D=5。

圖10顯示了覆蓋多孔介質前后尾跡剪切層分布的變化。從圖中可以看到：多孔介質使剪切層變得更細長、彼此相距更遠，繼而有利于減弱兩側剪切層的相互作用，同時平均渦量的峰值更加遠離圓柱。

圖8 平均速度U分布Fig.8 Average velocity U distributions

圖9 流向平均速度Umean曲線Fig.9 Ensemble-averaged velocity profiles in the streamwise Umean

圖10 平均渦量圖Fig.10 Mean vorticity structure contours

3.3 瞬時流場

圖11利用煙線試驗展示了有無多孔介質的圓柱尾跡渦結構的變化。對比圖11(a)和(b)可以看出，光柱尾部流場呈現(xiàn)出大尺度的卡門渦脫落，尾跡展向擺動明顯；覆蓋多孔介質后流場中的大尺度渦結構被抑制，尾跡流場變得相對穩(wěn)定。PIV試驗得到的瞬時速度場(圖12)可以得出同樣的結論。同時，覆蓋多孔介質后，圓柱尾部的瞬時流場速度被顯著減小，尤其是近尾跡區(qū)域，有助于減緩兩側渦相互作用的速率，繼而會導致渦脫落頻率降低(如3.1節(jié)中所述)。

圖11 流動結構對比Fig.11 Comparison of detailed flow structures

通過圖13和圖14中流向和側向速度脈動曲線可以看出，覆蓋多孔介質后，平面上兩方向的速度脈動均明顯降低。多孔柱體Urms曲線雙峰間距增大，Vrms曲線單峰變雙峰，再次說明其尾跡被拓寬。并且，光柱流向和側向速度脈動沿主流方向先增大后減小，而多孔柱體在所測范圍內一直增大。這表明多孔介質使得渦脫落的位置延后，與3.1節(jié)的渦脫落頻率的分析結果相吻合。

圖13 光柱尾跡流向和側向速度脈動曲線(Urms和Vrms)Fig.13 Velocity fluctuations profiles in the streamwise Urms and spanwise Vrms in the wake of the bare cylinder

圖14 多孔柱體尾跡流向和側向速度脈動曲線(Urms和Vrms)Fig.14 Velocity fluctuations profiles in the streamwise Urms and spanwise Vrms in the wake of the cylinder with porous layer coating

用渦形成長度可以衡量渦在尾跡中的形成位置，該值影響尾跡的壓力分布和對柱體的作用力[13]。沿尾跡中心線，將柱體中心與流向速度脈動峰值之間的距離定義為渦形成長度[14-15]。本次試驗測得覆蓋多孔介質后渦形成長度由光柱的2.91D延長至6.58D。

圖15給出了有無覆蓋多孔介質的圓柱尾跡的雷諾應力分布，從圖中可以看出，覆蓋多孔介質的柱體的雷諾應力明顯減小，其峰值更加遠離圓柱，同時展向距離增大。多孔介質的引入削弱了剪切層的不穩(wěn)定性，使得渦脫落位置延后，增大了渦形成長度。

圖15 雷諾應力圖Fig.15 Reynolds stress contours

4 結 論

通過熱線、煙線和PIV，研究了多孔介質對圓柱尾跡的頻率特征、平均流場和瞬時流場的影響，得出如下結論：

1) 多孔介質可以降低渦脫落頻率，有效抑制尾跡的大尺度渦脫落結構，使得渦脫落位置延后；此外，多孔介質可以減弱近尾跡場的速度PSD峰值(x/D≤1.5)，但會增大遠場的速度PSD峰值(x/D>1.5)。

2) 多孔介質減小了尾跡區(qū)的速度，拓寬了圓柱的尾跡寬度，可以減緩兩側渦相互作用的速率；同時多孔介質使得兩側剪切層變得細長，減弱了兩者的相互作用，平均渦量的峰值更加遠離圓柱。

3) 多孔介質的引入減小了尾跡區(qū)的速度脈動以及雷諾應力，使得其峰值更加遠離圓柱，減弱了剪切層的不穩(wěn)定性，增大了渦形成長度。

[1]Fujisawa N, Takeda G. Flow control around a circular cylinder by internal acoustic excitation[J]. Journal of Fluids and Structures, 2003, 17(7): 903-913.

[2]Wiplier O, Ehrenstein U. Numerical simulation of linear and nonlinear disturbance evolution in a boundary layer with compliant walls[J]. Journal of Fluids and Structures, 2000, 14(2): 157-182.

[3]Lam K, Lin Y F. Large eddy simulation of flow around wavy cylinders at a subcritical Reynolds number[J]. International Journal of Heat and Fluid Flow, 2008, 29(4): 1071-1088.

[4]Heenan A F, Morrison J F. Passive control of pressure fluctuations generated by separated flow[J]. AIAA Journal, 1998, 36(6): 1014-1022.

[5]Jimenez J, Uhlmann M, Pinelli A, et al. Turbulent shear flow over active and passive porous surfaces[J]. Journal of Fluid Mechanics, 2001, 442: 89-117.

[6]Sueki T, Ikeda M, Takaishi T. Aerodynamic noise reduction using porous materials and their application to high-speed pantographs[J]. Quarterly Report of RTRI, 2009, 50(1): 26-31.

[7]Sueki T, Takaishi T, Ikeda M, et al. Application of porous material to reduce aerodynamic sound from bluff bodies[J]. Fluid Dynamics Research, 2010, 42(1): 015004.

[8]Ruck B, Klausmann K, Wacker T. The flow around circular cylinders partially coated with porous media[C]//Porous Media and its Applications in Science, Engineering, and Industry: Fourth International Conference. AIP Publishing, 2012, 1453(1): 49-54.

[9]Klausmann K, Ruck B. Pressure distribution around circular cylinders with porous layer coating in the lee and flow analysis[C]//Eighth International Colloquium on Bluff Body Aerodynamics and Applications, 2016.

[10]Liu H, Wei J, Qu Z. The interaction of porous material coating with the near wake of bluff body[J]. Journal of Fluids Engineering, 2014, 136(2): 021302.

[11]Liu H, Wei J, Qu Z. Prediction of aerodynamic noise reduction by using open-cell metal foam[J]. Journal of Sound and Vibration, 2012, 331(7): 1483-1497.

[12]Liu P S. Introduction to porous materials[M]. Beijing: Tsinghua University Press, 2004: 1-14. (in Chinese)劉培生. 多孔材料引論[M]. 北京: 清華大學出版社, 2004: 1-14.

[13]Jung J H, Yoon H S. Large eddy simulation of flow over atwisted cylinder at a subcritical Reynolds number[J]. Journal of Fluid Mechanics, 2014, 759: 579-611.

[14]Szepessy S, Bearman P W. Aspect ratio and end plate effects on vortex shedding from a circular cylinder[J]. Journal of Fluid Mechanics, 1992, 234: 191-217.

[15]Lam K, Wang F H, So R M C. Three-dimensional nature of vortices in the near wake of a wavy cylinder[J]. Journal of Fluids and Structures, 2004, 19(6): 815-833.

Experimental study on the wake of a circular cylinder with porous layer coating

Wei Zheng1,2, Xia Chao1,2, Yuan Haidong1,2, Li Qiliang1,2, Yang Zhigang1,2,*

(1.ShanghaiKeyLabofVehicleAerodynamicsandVehicleThermalManagementSystems,Shanghai201804,China;2.ShanghaiAutomotiveWindTunnelCenter,TongjiUniversity,Shanghai201804,China)

The effects of porous media on the flow around a circular cylinder were investigated experimentally using the synchronous measurement of hot-wire anemometry, smoke-wire visualization and particle image velocimetry (PIV). The frequency characteristics, time-average and instantaneous flow field of the wake were measured and analyzed for improving understanding on the mechanism of flow modification caused by porous layer coating. The results show that the large-scale vortex shedding structure is effectively inhibited with porous layer coating and the Strouhal number of vortex shedding is 0.162 for the porous one which is smaller than 0.182 for the bare one.Additionally, the peak power spectral density of velocity is decreased in the near wake (x/D≤1.5), while increased in the far wake(x/D>1.5) and the position of obvious vortex shedding shifts downstream.Porous media attenuate the velocity of the wake region and widen the wake, leading to the weaker interaction between the two slenderer shear layers and slower rate of vortex interaction. Meanwhile, the porous media significantly eliminate the velocity fluctuations in the streamwise and spanwise, especially in the centerline as well as the Reynolds stress in the wake. The positions of peak value of Reynolds stress and time-averaged vorticity are more far away from the cylinder, so that the porous media weaken the instability of the shear layer and increase the vortex formation length.

cylinder; porous coating; wake; PIV

0258-1825(2017)02-0265-06

2016-12-26;

2017-01-23

國家重點研發(fā)計劃課題(2016YFB1200503-04)

魏崢(1993-)，女，碩士研究生，研究方向：鈍體流動，流動控制. E-mail：xwz@#edu.cn

楊志剛*，男，教授，博導，研究方向：車輛空氣動力學，流動控制. E-mail：zhigang.yang@sawtc.com

魏崢, 夏超, 袁海東, 等. 覆蓋多孔介質的圓柱尾跡實驗研究[J]. 空氣動力學學報, 2017, 35(2): 265-270,289.

10.7638/kqdlxxb-2016.0169 Wei Z, Xia C, Yuan H D, et al. Experimental study on the wake of a circular cylinder with porous layer coating[J]. Acta Aerodynamica Sinica, 2017, 35(2): 265-270,289.

V211.7

A doi: 10.7638/kqdlxxb-2016.0169