波流耦合作用下高樁擋板透空式防波堤水動力特性研究

張 彪，王心玉，宋軍港，趙 洋，劉 勇

(1.中國海洋大學 工程學院，青島 266100；2.中交天津港灣工程設計院有限公司，天津 300457)

高樁擋板透空式防波堤是一種環境友好的樁基海岸結構[1-5]，從高樁碼頭演化而來，在迎浪側設置擋浪板、胸墻進行擋浪，胸墻上方設置反浪弧來減小越浪量，在背浪側設置靠船構件兼顧碼頭功能。該結構輕便，成本較低，具有良好的透水性，使得港池內外的水交換比較暢通，尤其適用于深水地區和軟基地區。當擋浪板下方水流流速過大，還需要設置護底塊石來防止樁基沖刷。

許多學者研究了波浪作用下高樁擋板透空式防波堤的水動力特性。一些學者[6-9]基于勢流理論建立了波浪對高樁擋板透空式防波堤作用的二維理論模型，研究波浪參數和結構參數的變化對高樁擋板透空式防波堤透射特性和反射特性的影響，發現隨著波長的增大，透射系數逐漸增大，反射系數逐漸減小；隨著擋浪板入水深度的增加，透射系數逐漸減小，反射系數逐漸增大。也有許多學者[1, 10-11]通過物理模型試驗研究了波浪作用下防波堤的透射特性，發現隨著相對堤寬的增大，透射系數逐漸減小。HUANG和CHEN[2]試驗研究了不規則波作用下高樁擋板透空式防波堤的透射特性，發現透射系數在0.2～0.4變化，且隨波陡的增加而減小。許忠厚等[12]試驗研究了不規則波作用下擋浪板相對入水深度對高樁雙擋板防波堤透射系數的影響，發現透射系數隨著擋浪板相對入水深度的增大而減小。相比物理模型試驗，數值模擬成本更低，而且可以更好地揭示流場變化。ZHAI等[4]模擬了孤立波對高樁擋板透空式防波堤上部結構的作用，并對其流場進行了分析。王國玉等[5]基于OpenFOAM?研究了規則波對高樁擋板透空式防波堤上部結構作用，發現擋浪板入水深度越大，波浪透射越小，擋浪板下方附近流速越大。桂勁松等[3]考慮了上部結構擋板開孔的情況，發現隨著開孔率的增大，透射系數會先減小后增大，擋板開孔率為4%時，波浪透射最小。

在實際工程中，除了波浪的作用，海流、潮流等也會對結構物產生影響。當波流耦合作用時，其流場結構、剪應力特性和能量分布與單純波浪或水流作用相比發生了很大變化，因此有必要研究波流耦合作用下結構物的水動力特性。XIAO等[13]基于 RANS(Reynolds-Averaged Navier-Stokes)方程在數值波浪水槽底部設置入流速度邊界和出流速度邊界造流，實現了波流的非線性耦合。ZHANG等[14]同樣基于 RANS 方程建立了數值波流水槽，并利用 UMEYAMA[15]的試驗數據進行了驗證。SINGH等[16]利用快速傅立葉變換方法對試驗數據進行處理，發現波流耦合作用時，能量集中分布在入射波頻附近，分布規律與波浪單獨作用時相近。HUANG[17]試驗研究了波流耦合作用下開孔板的透射特性和反射特性，發現波流逆向時，水流的存在可以有效地增大波浪的能量耗散，降低波浪透射。HUANG[18]進一步試驗研究了波流逆向作用時雙排開孔板的透反射特性，同樣發現與波浪傳播方向相反的水流可以明顯降低波浪透射。

與以往研究不同，本文將開展物理模型試驗，并基于開源程序OpenFOAM?，數值模擬分析波流耦合作用下高樁擋板透空式防波堤的水動力特性。通過引入界面幾何重構方法 IsoAdvector，開發可計算帶多孔介質的透空式結構物的求解器 porousWaveIsoFoam，建立還原結構全模型的數值波流水槽，新建立的數值波流水槽可解決原有界面壓縮方法捕捉界面褶皺較為嚴重的問題。基于建立的數值模型研究波流耦合作用下結構的透射特性、群樁間的流速變化、結構物前后的流速變化、結構物前后波浪能量分布以及結構物周圍動壓場和速度場的變化，為樁基透空式結構物的工程設計和優化提供科學指導。

1 數學模型

1.1 控制方程

將流體假設為不可壓縮粘性流體，滿足連續方程和動量方程

(1)

(2)

式中：u是流體的速度；p是動水壓力；ρ是密度；g是重力加速度；μ是動力粘性系數。

將方程中的變量體積平均，并且考慮多孔介質的影響，連續方程和動量方程變為

(3)

(4)

(5)

(6)

式中：d50是多孔介質的中值粒徑，Keulegan-Carpenter 數KC=ucT/nd50，根據JENSEN等[20]的研究，uc是水質點最大振蕩速度，T是振蕩周期。αp和βp是經驗系數，本文采用 JENSEN等[20]的推薦值βp=2.0，αp=500。本文中的數值算例均采用LARSEN和FUHRMAN[21]提出的修正后的k-ω模型進行模擬，可有效解決波高衰減問題，在長時間模擬中得到穩定的波列。OpenFOAM?中數值水槽的自由液面通常采用界面壓縮方法來計算確定，但是界面壓縮方法捕捉到的界面是褶皺的。本文采用 ROENBY等[22]提出的幾何VOF方法，界面幾何重構方法IsoAdvector，捕捉到的水氣界面更加平滑，更加符合物理實際。

1.2 數值造波、造流和消波方法

采用速度入口邊界方法造波、造流，在側面入口邊界上實時給定波面位置和流速，實現水槽的數值造波。在數值水槽底部設置入流速度邊界條件和出流速度邊界條件，使水槽中形成穩定的剪切流，在傳播過程中受到邊界條件的影響，可以完成波流的非線性耦合，這與實驗室中物理水槽造流的方法相似。

數值模型采用JACOBSEN等[23]提出的松弛區域方法消波，在數值波流水槽的兩側設置松弛區域，可以有效防止入口邊界處波浪二次反射和出口邊界處波浪的反射，在松弛區域內每一時刻進行如下修正

φ=(1-αr)φt+αrφc

(7)

式中：φt為期望得到的目標速度、壓力、體積分數，φc為通過數值模擬得到的速度、壓力、體積分數，αr為與空間位置有關的加權函數，αr∈[0,1]，且滿足

(8)

式中：χ代表松弛區域內的局部坐標系，其值范圍在0～1。

2 數值模型驗證

2.1 網格收斂性檢驗

首先考慮了波流耦合后空水槽中波浪傳播的情況，數值波流水槽的布置如圖 1 所示，水槽的長度、寬度和高度分別為 23 m、0.4 m、1.0 m，兩端松弛區長度皆為 3 m。采用與波流對高樁擋板透空式防波堤作用數值模擬相同的波要素(h=0.4 m，H=0.05 m，T=1.2 s，Uc=0.08 m/s )進行網格收斂性檢驗。采用 blockMesh 工具建立三維數值水槽的背景網格△x=0.04 m、△y=0.012 m、△z=0.04 m，然后進一步加密區域 0 m

圖1 數值波流水槽布置示意圖

表1 加密區域的網格基本信息

在x=11 m 處設置波高儀，得到三種網格下波面的計算結果。如圖 2 所示，Mesh II 和 Mesh III 的波面歷時曲線吻合良好，而較為粗糙的網格 Mesh I 會一定程度上低估波峰、高估波谷的結果，原因是網格密度過小會造成比較嚴重的數值耗散。圖 2 中的對比結果表明 Mesh II 和 Mesh III 兩種網格密度都能得到可靠的數值計算結果，綜合考慮計算結果的準確性和計算效率，本文決定采用Mesh II來離散計算域。根據VANNESTE和TROCH[24]的研究，網格的劃分與波浪特性有關。在本算例中，在自由面附近的加密區域中，x方向的網格尺寸 Δx與波長L有關，這里將其設置為L/97，y方向的網格尺寸 Δy與波高H有關，這里將其設置為H/9，令z方向的網格尺寸Δz=Δx。在后文波浪、波流對高樁擋板透空式防波堤作用的數值模擬中，在自由面附近的加密區域，Δx=Lmin/97，Δy=Hmin/9，Δz=Δx，其中，Lmin和Hmin分別是模擬工況中波長和波高的最小值。對于波陡較大的工況，這種網格劃分方法可以獲得更加精細的網格，從而保證計算結果的準確性。

圖2 x=11 m 處不同網格所得數值結果對比

2.2 數值波流水槽驗證

如圖3 所示，二維數值波流水槽長25 m、高0.6 m，水槽兩端松弛區長度為3 m，波浪向右傳播，波高H=0.023 4 m，周期T=1 s。底部設置入流速度邊界條件和出流速度邊界條件，速度大小分別為0.05 m/s和-0.05 m/s ，水槽內形成平均流速為Uc=0.08 m/s 的穩定水流。水平方向設置5 000個網格，網格長度Δx=0.005 m，豎直方向設置400 個網格，網格高度 Δy=0.001 5 m。數值模型模擬時間采用ZHANG 等[14]的推薦，運行25 s(25個波浪周期)。監測水槽x=15 m處的波面和流速，取20 s之后5個波浪周期的數值結果的平均值與 UMEYAMA[15]的試驗數據對比。

圖3 數值波流水槽布置示意圖

圖4中給出了波流耦合作用下x=15 m 處的波面歷時曲線圖，其中η為自由液面。圖5中給出了同一波浪周期不同時刻(t/T=0、0.25、0.50、0.75)水平速度沿水深分布圖，其中，t為每個周期內的時間。數值模型所得波面歷時曲線和水平速度沿水深分布曲線均與UMEYAMA[15]的試驗數據符合良好，說明本文數值模型能夠有效捕捉波流耦合作用下的自由液面特征，合理模擬波流之間的非線性作用。

圖4 波流耦合作用下x=15 m 處的波面歷時曲線圖

5-a t/T=05-b t/T=0.255-c t/T=0.505-d t/T=0.75

2.3 波浪對高樁擋板透空式防波堤的作用

波浪對高樁擋板透空式防波堤作用的物理模型試驗在山東省海洋工程重點實驗室完成，實驗室水槽長60 m、寬0.8 m、高1.5 m。水槽的一側安裝了推板式造波機，由電腦控制可以產生目標波浪，另一側設置了斜坡型的消波網來消減波浪。高樁擋板透空式防波堤結構的幾何尺寸以及群樁分布情況見圖6和圖7，試驗比尺為1：25。擋浪板底部設置護底塊石，護底塊石的高度為0.064 m，中值粒徑d50=0.055 m，孔隙率為0.42。

圖6 高樁擋板透空式防波堤結構示意圖(單位：mm)

圖7 群樁布置示意圖(單位：mm)

三維數值水槽長23 m、高1.0 m、寬0.4 m(高樁擋板透空式防波堤結構具有對稱性，導致流場分布具有對稱性，因此在數值模型中引入symmetry對稱邊界條件，將數值水槽寬度簡化為0.4 m，可以大大提高三維數值模擬的計算效率)，兩端松弛區長度均為3 m。圖7中防波堤的寬度Lw=2.0 m，放置在數值水槽11 m

8-a 結構物前8-b 結構物后

本節模擬了波浪對高樁擋板透空式防波堤的作用，并利用開展的物理模型試驗數據進行驗證，透射波面的對比如圖9所示。結構物前水深為0.400 4 m，結構物后水深為 0.480 4 m。圖9-a中，入射波的波高H=0.188 0 m，周期T=1.952 s，圖9-b中，入射波的波高H=0.113 2 m，周期T=1.232 s。從圖9中可以看到，波浪作用下高樁擋板透空式防波堤的透射波面數值結果與試驗結果吻合良好，驗證了波浪對高樁擋板透空式防波堤作用的三維數值模型的合理性。

9-a H=0.188 0 m，T=1.952 s 9-b H=0.113 2 m，T=1.232 s

2.4 水流對高樁擋板透空式防波堤的作用

結構物前水深為0.4 m，結構物后水深為0.48 m。在3.5 m

圖10 流速0.08 m/s的水流對結構物作用后結構物前后的水位變化

3 結果討論與分析

3.1 透射波面與群樁間流速變化

波流對高樁擋板透空式防波堤作用的數值模擬示意圖如圖11所示，對入流速度邊界條件和出流速度邊界條件附近的局部流場圖進行了放大，其中，結構物前水深為0.4 m，結構物后水深為0.48 m，水槽中水流的平均流速Uc=0.08 m/s，規則波波高H=0.05 m，周期T=1.2 s，圖例中利用重力加速度g和結構物前的水深h對速度進行了無量綱化處理。數值模擬時間為40 s，取25～40 s的數據進行研究，從圖10中可以看到，此時，數值水槽已經達到穩定狀態。

注：左側箭頭所指為入流速度邊界條件附近流場的局部放大圖，右側箭頭所指為出流速度邊界條件附近流場的局部放大圖。

圖12給出波流耦合作用與波浪單獨作用下高樁擋板透空式防波堤透射波面的對比圖，可以看到，波流耦合作用下結構物的透射波高明顯大于波浪單獨作用時的透射波高，原因是水流的存在可以使波長變長[14]，對于高樁擋板透空式防波堤，波長越大，波浪透射越大。

圖12 波流耦合作用和波浪單獨作用下透射波面對比

圖13給出一個波浪周期內不同時刻(t/T=0、0.25、0.50、0.75)群樁中部位置(x=11.898 m處)群樁間水平速度沿水深分布圖。波浪單獨作用時，由于高樁擋板透空式防波堤擋浪效果良好，波浪透射小，群樁間的流速小，且波峰作用時(圖13-b)流速方向與波浪傳播方向一致，波谷作用時(圖13-d)流速方向與波浪傳播方向相反。波流耦合作用時，高樁擋板透空式防波堤對水流的阻擋效果較弱，群樁間的流速較大，因此，當波流耦合作用時，需要進一步考慮樁基沖刷問題。

13-a t/T=013-b t/T=0.2513-c t/T=0.5013-d t/T=0.75

3.2 高樁擋板透空式防波堤前后的流速變化

圖14和圖15對比了波流耦合作用下和波浪單獨作用下一個波浪周期內不同時刻(t/T=0、0.25、0.50、0.75)高樁擋板透空式防波堤前(x=10.97 m 處)、后(x=13.20 m 處)的水平流速沿水深分布情況。從圖14中可以看到，相比于波浪單獨作用，波流耦合作用下結構物前流速明顯增大，水平流速沿水深方向先增大后減小，擋浪板下方海床附近由于多孔介質的存在，流速較小，護底塊石和擋浪板之間的縫隙(0.064 m

14-a t/T=014-b t/T=0.2514-c t/T=0.5014-d t/T=0.75

圖15 波流耦合作用和波浪單獨作用下結構物后水平速度沿水深分布(x=13.20 m 處)

3.3 高樁擋板透空式防波堤前后的波浪能量分布變化

波流耦合作用和波浪單獨作用下結構物前后能量分布對比如圖16所示，圖中，頻率f用譜峰頻率fp無量綱化，對應頻率間隔△f內組成波的平均能量S(f)用波浪單獨作用時的譜峰頻率組成波的平均能量S(fp)w無量綱化。與波浪單獨作用相比，波流耦合作用時結構物前后能量在各個頻率上的分布規律并不發生改變，能量的大小發生變化。波浪總能量通過在整個頻率內對譜積分算得。波流耦合作用時結構物前的波浪總能量略大于波浪單獨作用的情況。如圖16-b中所示，透射波浪中出現了能量占比較大的低頻組分。結合圖12，可以發現，由于水流的存在，波浪透射明顯增大，透射波的總能量顯著增大。

16-a 結構物前x=9 m 處16-b 結構物后x=15 m 處

3.4 流場特性分析

在圖18-a和圖18-b中可以發現，波浪單獨作用時，在結構的前方出現了順時針的渦旋，在結構的后方出現了逆時針的渦旋(見圖18-a、圖18-b中的局部流場放大圖)，這也可能會導致動壓的減小。結合圖17-c、圖17-d以及圖18-c、圖18-d發現，波流耦合作用時，由于結構物后方流速較大，出現了較大的負壓。

17-a 波浪單獨作用

18-a 波浪單獨作用18-b 波浪單獨作用18-c 波流耦合作用18-d 波流耦合作用注：箭頭所指部分為流場的局部放大圖。

4 結論

基于 OpenFOAM?建立了波流對高樁擋板透空式防波堤作用的三維數值模型，引入界面幾何重構方法，可以捕捉到更加平滑、更加符合物理實際的自由液面，采用修正后的k-ω湍流模型，可以有效解決波高衰減問題，在長時間模擬中得到穩定的波列。基于建立的數值模型，分析了波流耦合作用下高樁擋板透空式防波堤的水動力特性。研究發現：

(1)單獨水流對高樁擋板透空式防波堤作用時，結構物前水位會升高，結構物后水位會降低；波流耦合作用時，結構物前后波能隨波頻的分布規律與波浪單獨作用下相似，但波浪透射顯著增大，透射能量顯著增大。

(2)波浪單獨作用時，擋浪板前方和靠船構件后方附近渦脫落現象明顯，會引起動壓的減小，波流耦合作用時，水質點運動速度明顯增大，渦脫落現象不再出現，波峰作用于擋浪板和胸墻，動壓顯著增加。

(3)波流耦合作用時，高樁擋板透空式防波堤前后的水平流速沿水深方向先增大后減小，靠船構件后方的流速、擋浪板底部和護底塊石頂部縫隙之間的流速明顯增大，護底塊石可以有效減小擋浪板下方的流速；波流耦合作用下，群樁間的流速顯著增大，且越靠近自由表面附近流速越大。