基于PIV技術不同形式防波堤周圍渦旋場特性研究

郄祿文，杜闖，張翔，蔣學煉

（1.河北大學 建筑工程學院，河北 保定 071002；2.天津城市建設學院 天津 300384）

PIV（Particle Image Velocimetry），又稱粒子圖像測速法，是20世紀70年代末發展起來的一種瞬態、多點、無接觸式的基于流場圖像分析的非接觸式二維流場測量技術。PIV技術的特點是超出了單點測速技術（如LDV）的局限性，能在同一瞬態記錄下大量空間點上的速度分布信息，并可提供豐富的流場空間結構以及流動特性。PIV技術除向流場散布示蹤粒子外，所有測量裝置并不介入流場，有能夠實現無擾動精確有效測量二維流場的流速分布形式的特點。粒子成像測速技術是流場可視化研究中的一個最新發展，它是伴隨著光纖、微電子、激光、數字信號處理和圖形圖像處理等高新技術的發展而發展起來的，目前在各領域應用得非常廣泛。

多年來，渦旋波流動的理論和應用研究一直受到很多學者的關注。但由于早期實驗流體力學理論和測試技術的限制，以前的研究者大都采用相機對流場進行拍照的實驗手段，對渦旋波流場進行了一些定性的分析和比較，還未能獲得渦旋波形成機理的精確描述。Pedley等（1985）用平行光對上壁面周期性壓縮的槽道內流體進行拍照，觀察并得到了渦流波流場的靜態照片。Sobey（1985）用配有250W閃光燈的尼康相機在周期震蕩的流體經過的槽道中，用熒光珍珠粉作為示蹤粒子，在柯達膠卷上拍攝了渦旋波流場的靜態照片，定性的分析了渦旋波流場的特性。PIV技術在不干擾槽道內部原始流場分布的狀態下進行流速測定，具有較高的準確度和精確性。劉鳳霞等（2006）運用PIV技術分析出渦旋波流場的速度和剪力分布。劉愛華等（2006）運用PIV技術分析了渦流中間包流場的特性，分析得出了一些重要的結果。董敏（2007）利用PIV技術研究破碎波對直墻建筑物作用過程，對破碎波作用力特點進行研究。龍曉警（2009）研究了PIV技術在水槽波浪中的研究，分析了整個測試流場的瞬時速度場。袁德奎（2004）對半圓型防波堤水動力學特性進行了數值模擬的方法和應用研究。劉長根等人研究了用時變雷諾方程模擬孤立波與有臺階的半圓型防波堤的相互作用，并研究了SPH方法在水波與半圓型防波堤之間作用。

深水防波堤周圍流場存在復雜的紊動變化，防波堤較容易發生破壞。本文采用PIV流場測試技術，對周期性規則波流經防波堤槽道所對應的流場進行了實驗測試研究，獲得了規則波流場的瞬時速度分布。通過對實驗數據的分析和整理，定量地研究了防波堤堤身前后渦旋波流場的特性。

1 實驗設計與參數

1.1 實驗布置與流程

實驗在大連理工大學海岸和近海工程國家重點實驗室進行。實驗流程見圖1、圖2。實驗裝置主要由四部分組成：槽道、往復泵，PIV系統和外觸發裝置。槽道由有機玻璃制成，以實現流動可視（圖 3）。

實驗模型擬采用三種堤型：直立堤、1/4圓弧堤、半圓堤斷面型式（圖4-圖6）。

圖1 PIV系統示意圖

圖2 PIV實驗流程圖

圖3 實驗現場圖

圖4 直立堤斷面示意圖

圖5 1/4圓弧堤斷面示意圖

圖6 半圓堤斷面示意圖

表1 防波堤實驗波浪水力條件

如圖7試驗時的采集點位于No.4和No.5，當圖左側造波機開始造波，波浪經過No.5號波高儀，并在終端上反應出3個完整波浪周期圖像時，開始采集；當終端屏幕上全部呈現周期規則波浪圖像時，停止采集。波高儀采集間隔0.02 s，數據截選時間16 s，即800次的樣本采集所得到的數據區間。圖7右側采用碎石搭建斜坡，用于消減回浪，避免二次浪干擾實驗數據。

1.2 渦量的定義

渦量是一個矢量，在數值上等于速度矢量的旋度，其定義如下：

式中：Ω為流場中量測點上的渦量值，U、W分別為水質點的水平速度與垂直速度，而其計算的方法則是利用中心差分而得：

圖7 實驗布置示意圖

其中：Ωp為位置在（i，k）時之渦度值，i，k分別代表x，z方向的陣列碼。Ωp為正值時，表示渦流逆時針旋轉；Ωp為負值時，表示渦流順時針旋轉。

2 實驗數據處理及分析

利用PIV技術方法，本實驗對6種工況進行斷面波浪實驗，每種工況采集45個瞬時（間隔0.15 S）時刻各測點的速度值，采樣區測點數目分布為48×35=1 680個測點。對測點速度值和渦量值進行分析，得到該時刻速度場圖和渦量場等值線圖。

2.1 速度與渦量的相對關系

通過對計算機采集的測點實驗數據，可得每一時刻各測點的速度值，并進而可計算觀測區各點渦量及渦量場的分布。圖8是半圓堤在第一種水力狀態下的同一幀的速度場，渦量場云圖和渦量場的對比。

由圖8可以觀察得出速度集中和速度變化率較高的地方，容易形成渦量值高的區域。同時也印證了前面的公式，速度與渦量之間的內在聯系。

2.2 波面過程與渦量分布

試驗時PIV都采集3 S，為得到大約3個周期的過程，對應45幀則是每幀間隔1/15 S。波面和PIV圖像是同時采集的，但位置并不完全重疊。波高儀可監測反射系數和透射系數，以便比較不同渦流情況下的反透射變化。波高儀采集次數為800次，間隔0.02 S。如圖9所示。

實驗測得半圓堤波面過程曲線和渦量值代數均值曲線，如圖10、圖11、圖12所示。

以上結果看出，渦量代數均值曲線和波面過程曲線并不重合，并得到如下規律：

圖9 采集數據區域圖

圖10 水力條件1下半圓堤迎、背浪面波況與渦量代數均值曲線

圖11 水力條件2下半圓堤迎、背浪面波況與渦量代數均值曲線

圖12 水力條件4下半圓堤迎、背浪面波況與渦量代數均值曲線

（1）堤前波面過程曲線與堤后曲線呈相反分布。

（2）波浪的波高變化并不會改變渦量代數均值曲線和波面過程曲線的相關性，但相關性變得更加明顯。

（3）波浪的不同周期影響了原有渦量代數均值曲線和波面過程曲線的相關性。

2.3 渦量數值特性分析

2.3.1 隨機測點渦量代數值分析

為了探討渦量值的形成規律，分析半圓堤背浪面隨機挑選不同行不同列測點（同列測點1，測點2，測點3的水深是依次加大；同行測點1，測點2，測點3依次遠離防波堤），在采樣時間段內的堤前與堤后渦量值曲線圖，采樣時間內的渦量代數值曲線具有相同的規律，同列測點間的曲線滿足同一規律，曲線的走向近似相同。1/4圓弧堤、直立堤情形呈現相同規律。

圖13 半圓堤同行、同列隨機測點渦量代數值曲線圖

2.3.2 渦量代數值最值分析

隨后對相同水力條件下三種防波堤迎浪面和背浪面的渦量最值進行分析。在采樣區間內，分別得出每一幀的最大值和最小值曲線圖，如圖14、圖15所示。

三種防波堤的實驗結果定量分析表明，迎浪面區域，半圓堤產生的渦量最大值和最小值分別是1/4圓弧堤的86%和117%。背浪面區域，半圓堤產生的渦量最大值和最小值分別是1/4圓弧堤的40%和30%，直立堤產生的渦量最大值和最小值分別是1/4圓弧堤的31%和33%。

2.3.3 渦量均值分析

對水力條件1下，三種防波堤在采樣區間內每一幀上所對應的采樣區間上所產生的渦量均值進行分析。將45幀平均值繪成曲線圖，如下：

圖14 半圓堤、1/4圓弧堤、直立堤迎浪面最值曲線圖

圖15 半圓堤、1/4圓弧堤、直立堤背浪面最值曲線圖

由圖16可以看出，均值曲線均呈現出正弦曲線的特性，堤前較堤后特征更加明顯。可見，防波堤背浪面貼近防波堤表面的地方，渦旋運動非常微弱，遠離防波堤處，底部渦旋運動變得劇烈；迎浪面情況正好與之相反，貼近海面渦量波動劇烈。其他部位測點渦量波動基本一致，并且變化平緩。

2.3.4 不同波高下渦量代數值分析

對于同種防波堤，不同水力條件下，渦量的數值差異研究，本次試驗選取了半圓堤，在僅改變波高的情況下，進行數值差異研究（圖17、18）。

圖16 防波堤迎浪面、背浪面渦量均值曲線

由圖17、圖18可以看出，波高和渦量代數極值成正相關關系。波高越高，產生的極大值越大（極小值越小）。均值曲線呈正弦曲線特征，迎浪面數值曲線相對背浪面數值曲線較平穩，特征較明顯。

2.3.5 不同周期下，三種防波堤的渦量特性分析

分別選取水力條件1和4的實驗數據進行渦量數值分析（圖19）。

由圖19可以看出，當波浪的周期減小時，半圓堤和1/4圓弧堤所產生的渦量積累相應的增加，而直立堤的曲線反映出來的情況正好相反，會隨著周期的減小渦量的積累相應的減少。所有的代數均值曲線均反映正弦曲線特征。隨著波浪周期降低，在防波堤背面會生成更多的渦量，1/4圓弧堤和直立堤曲線特征表現的比較明顯，呈現出正弦曲線的特征。周期的改變對正弦曲線的波動情況影響很大，兩條曲線之間近似差了半個周期的距離。

隨著波浪周期降低，堤前區域產生了較少的渦旋運動，而堤后正好相反。波浪周期的改變對渦量曲線的周期影響不明顯。由圖20可以得出，當周期從1.223 s降低為0.996 s（約減少18.6%），渦量代數均值曲線的振幅，半圓堤平均約增加的16%，直立堤平均約減少4.5%。

圖17 半圓堤迎浪面三種波高下渦量最小值、最大值、均值曲線

2.4 不同防波堤背浪面渦流形成與分布特性分析

對于三種防波堤而言，渦量的形成的阻礙效果各不相同。而在同種波況下，渦量值偏高的區域也各不相同。

圖18 半圓堤背浪面三種波高下渦量最小值、最大值、均值曲線

對于半圓防波堤來說，正渦量區域多集中在采樣區的左上部，即貼近防波堤背浪面的區域（圖21-1）。而負渦量區域多集中在采樣區的中場地帶（圖21-2）。即遠離防波堤的背浪面處多產生順時針渦旋，中場地帶多產生逆時針渦旋。而這些渦旋分布隨著時間的推移也會產生相應的移動，正渦量區會向近堤方向移動，原區域變為負渦量區（圖21-3）。一段時間后正渦量區又會回到原位，負渦量區回到中場。對于半圓防波堤來說，整個采樣時間段內出現0渦量區域的時間較少（圖21-4），采樣區間大多時間內受到渦旋的影響。

圖19 半圓堤、1/4圓弧堤、直立堤背浪面渦量對比

圖20 不同水力條件半圓堤、直立堤渦量均值對比

圖21-1 半圓堤1幀渦量場

圖21-2 半圓堤6幀渦量場

圖21-3 半圓堤15幀渦量場

圖21-4 半圓堤31幀渦量場

對于1/4圓堤來說，正負渦量區間常伴隨出現，并沒有明顯的正渦量區和負渦量區。開始時正負渦量區間共同出現在采樣區的左上側，即遠離防波堤背浪面（圖22-1）。隨后共同向采樣區中上部移動（圖22-2），到逐漸消失，最后又回到初始位置（圖22-3）。整個周期相對于半圓堤較長。另外相對于半圓防波堤，在整個采樣時間段內，采樣區域中部為0渦量區域的時刻增多（圖22-4），在采樣區域內渦旋生成相對于半圓防波堤較弱。

對于直立堤來說，正負渦量區域分布在采樣區上部（圖23-1），范圍較小，正負渦量區伴隨出現，隨著時間的推移正負渦量區域逐漸變小（圖23-2），經過一段時間后恢復初始狀態（圖23-3）。相對于半圓堤和1/4圓堤，直立堤的采樣區域內出現了最多的0渦量區，即采樣區間內形成渦旋的區域也相對較少。

圖22-1 1/4圓堤1幀渦量場

圖22-2 1/4圓堤10幀渦量場

圖22-3 1/4圓堤42幀渦量場

圖22-4 1/4圓堤13幀渦量場

圖23-1 直立堤1幀渦量場

圖23-2 直立堤5幀渦量場

圖23-3 直立堤26幀渦量場

3 結論

本文利用PIV技術，基于模型實驗，對規則波作用下防波堤表面的速度流場及其渦流特性進行了實驗研究，通過對不同型式防波堤在不同的水力條件下的波浪流場實驗數據的分析，定量地研究了防波堤堤身水域的渦旋流場的水力特性。

（1）區域內受到的波浪運動與區域內產生的渦量成相關性，這種相關性受到波浪周期和頻率的影響；防波堤會改變這種相關性。

（2）區域內同一時刻產生的渦量正負極值是對稱出現的。

（3）防波堤迎浪面貼近防波堤表面的地方，渦旋運動非常劇烈，遠離防波堤處，底部渦旋運動變得微弱，貼近海面渦量波動劇烈。迎浪面情況正好與之相反。

（4）波高和渦量代數極值成正相關關系。周期的增加會影響渦量的積累和產生。周期越小渦旋運動愈多，但是積累增加不大。渦量的均值曲線也滿足正弦性質，并且堤前部分特性更加明顯。

（5）防波堤底部因渦旋作用對堤底造成滑移影響，1/4圓弧堤較為劇烈，半圓堤次之，直立堤最弱。

（6）三種防波堤在同種水力條件下，在背浪面采集區域內，渦旋運動的遷徙方向和聚集分布基本相似，但是渦旋形成的劇烈程度和渦旋的大小卻各不相同。直立堤相對其他兩種堤型，背浪面采樣局域內渦旋運動相對較弱，遷徙周期相對較短。

Okajina A，1982.Strouhal Nunmbers of Retangular Cylinder.J Fluid Mech,123：379-398.

Sobey IJ,1985.Observation of Wavesduring Oscillatory Channel Flow.J.Fluid Mech,151：395-426.

Westerwell J.1997,Fundamentals of digital particle image velocimetery.Measurement Scienceand Technology,8(12):1379-1392.

董敏,2007.破碎波對直墻建筑物作用機理的試驗研.大連理工大學碩士論文,

段俐,康琦,申功忻,2000.PIV技術的粒子圖像處理方法,北京航空航天大學學報,2(1):79-82.

蔣學煉,李珊珊,李炎保,2012.防波堤滑移運動解析解及其應用研究,海洋通報,31(4):441-448.

李炎保,孔暉,蔣學煉,2007.沉箱式防波堤蛇行破壞的期望滑移量模型,海洋通報,26(2）78-87．

栗鴻飛,宋文武,2009.PIV技術在流動測試與研究中的應用,西華大學學報(自然科學學報）,9(5):27-31.

劉愛華,候勤福,王煥洋,等,2006.PIV測試技術在渦流中間包流場研究中的應用,材料與冶金學報,5(1):26-29.

劉鳳霞,劉志軍,王琳,等,2006.用PIV技術測試渦旋波流場的速度和剪應力分布實驗力學,3(6):378-284.

劉鳳霞,史啟財,劉志軍,等,2006.應用PIV技術測試渦旋波流場,試驗力學,21(5):611-616.

龍曉警,2009.粒子圖像測速(PIV）技術在水槽波浪中的研究,天津大學碩士論文．

齊鄂榮,盧煒娟,邱蘭,等,2007.應用PIV系統研究橫流中水平力柱繞流旋渦特征.武漢大學學報(工學版),40(1):19-24.

申功忻,1997.全場觀測技術概念、進展與展望．北京航天航空大學學報,23(3):332-340.

史小軍,王樹立,李田恩,等,2008.利用PIV測試水槽流暢的實驗研究.安徽工業大學學報,25(2):151-154.

孫鶴泉,沈永明,王永學,等,2004.PIV技術的幾種實現方法.水科學進展,15(1):105-108.

王軍,李國元,溫國珍,等,2004.水流空心渦帶流場的PIV測試與分析,華中科技大學學報(自然科學版）,12(12):66-68.