潛式雙層水平板型防波堤的消浪機理分析

李靖波，張寧川，劉愛珍

（1.大連元堃海洋科技有限公司，大連116000；2.大連理工大學海岸和近海工程國家重點實驗室，大連116024；3.交通運輸部天津水運工程科學研究所，天津300456）

潛式雙層水平板型防波堤的消浪機理分析

李靖波1，張寧川2，劉愛珍3

（1.大連元堃海洋科技有限公司，大連116000；2.大連理工大學海岸和近海工程國家重點實驗室，大連116024；3.交通運輸部天津水運工程科學研究所，天津300456）

雙層水平板型防波堤是一種新型結構形式的防波堤。該防波堤與波浪相互作用的流場和渦量場變化特性及消浪機理研究，文獻中還鮮有論述。文章用數值方法研究了波浪與潛式雙層水平板型防波堤的相互作用。以連續(xù)方程、雷諾時均方程和k-ε湍流模型作為控制方程，建立了波浪與防波堤相互作用的數值模型。數值研究了潛式雙層水平板型防波堤附近的流場、渦量場變化過程和不同位置的波能流。研究結果表明，上層板上方水域的渦動能量損耗及雙層板之間和下層板下方波能流的逆向傳遞是防波堤的主要消浪機理。

雙層水平板型防波堤；流場；渦；波能流；消浪機理

雙層水平板型防波堤是一種新型的環(huán)保型防波堤結構，與傳統防波堤相比有很多優(yōu)點：（1）允許水體在防波堤內外自由交換，不破壞海洋環(huán)境；（2）結構簡單，便于施工；（3）造價隨水深的變化不敏感。

近年來，水平板型防波堤越來越引起研究者們的興趣。在大家的努力下，已取得了許多有價值的成果。Heins［1］基于Wiener?Hopf技術較早研究了有限水深條件下，波浪與半無限長潛式水平板的作用問題，得到了反射系數和透射系數的解析表達式。Patarapanich［2］采用有限單元法數值研究了潛式水平板的反射和受力問題，分析了相對板寬、相對水深和相對潛深的變化對單層水平板的反射系數及波浪荷載的影響?？紤]到潮差較大時，水平板的淹沒深度將有較大的變化，此時單層板很難達到理想的消浪效果。鑒于此，許多學者對雙層水平板和多層水平板的消浪性能進行了研究。Patarapanich和Cheong［3］用試驗方法對雙層水平板型防波堤的反射系數和透射系數進行了研究。該防波堤結構的下層板潛入水下，上層板位于自由水面處并且放在下層板的下游。研究結果表明，相對板寬相同時，該雙層水平板結構的反射系數比任意淹沒深度的單層板要大得多，并且下層板的淹沒深度與水深的比值約在0.1～0.20時，雙層板的透射系數最小。Usha 和Gayathri［4］從線性勢流理論出發(fā)，研究了雙層水平板型防波堤對二維波浪的反射和透射情況。研究結果表明，相對板間距為0.22～0.23，板寬為波長的0.37～0.39時，透射系數最小。王國玉［5］用試驗方法對多層（11層）水平板防波堤的透射系數和反射系數的變化特征進行了研究。結果顯示，相對板寬比相對板間距對消浪效果的影響要明顯得多。Neelaman和Gayathri［6］通過試驗研究了單層板和雙層水平板型防波堤的水動力特性及波浪的透射系數和反射系數。研究結果表明，通常情況下，雙層板的水動力特性要優(yōu)于單層板，且能更有效的增加反射波浪和減少透射波浪。Li和Zhang［7］建立了波浪與潛式雙層水平板相互作用的數值模型，數值模型可以準確的模擬波面和波壓力，研究結果發(fā)現在堤后產生了高階諧波分量，而且堤后可見明顯的旋渦。李靖波和張寧川［8］通過數值計算研究了潛式雙層水平板防波堤的透射系數和波浪力，并擬合出了透射系數和波浪力的計算公式。

綜上所述，研究者關于波浪與水平板型防波堤的相互作用已經做了大量研究工作。研究內容涉及了反射、透射系數及波浪力的相關影響因素，但對水動力特性的研究較少，特別是關于水平板防波堤與波浪相互作用的流場和渦量場變化特性研究，文獻中還鮮有論述。研究流場和渦量場變化特性，可以更好的從能量耗散的角度，綜合解釋雙層水平板的消浪機理，從而推進該結構向實用化方向邁進。這正是本文的研究目的。本文數值研究了波浪與雙層水平板型防波堤的相互作用。數值研究了潛式雙層水平板防波堤周圍的流場和渦量場變化特性，并進一步探討了防波堤的消浪機理。

圖1 波浪與雙層水平板防波堤作用示意圖Fig.1 Sketch of waves interaction with a twin?plate breakwater

1 控制方程

當波浪與建筑物相互作用時，其經常破碎，水質點出現紊動摻混現象。紊動影響應該考慮進去。所以，本文采用連續(xù)方程，雷諾時均方程和k-ε湍流模型作為控制方程，并引入流體體積函數F（x，y，t）處理自由表面?？刂品匠碳癋方程的表達式可詳見文獻［8］。

2 數值方法

數值模型的計算網格劃分、方程離散形式、邊界條件及數值模型的驗證可詳見文獻［8］。

3 數值結果及分析

波浪與潛式雙層水平板型防波堤相互作用的示意圖可參見圖1。

圖2 防波堤附近的流場變化過程（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）Fig.2 Velocity field around the breakwater（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）

3.1 防波堤附近的流場變化過程

圖2給出了計算參數為周期T=1.4 s，波高H=0.06 m，水深D=0.48 m，板寬B=0.9 m，板間距S=10.5 cm條件下，潛式堤周圍的流場分布示例。圖中顯示，防波堤的存在對水質點流場影響很大，在t/T=0時刻，入射波浪在堤前與水平板上、下方和兩板中間產生的離岸回流相互作用，水質點速度方向垂直向上，板前端的波面出現雍高現象。從t/T=3/7時刻，波浪從水平板上方經過時，與板上方的水體發(fā)生了激烈的碰撞，導致波峰前側變陡，并在水平板的上方發(fā)生破碎。在t/T=1時刻，波浪傳至水平板的后方，在水面附近形成了強烈的向岸流，水平板后方產生了一個大的旋渦，同時有一部分水體斜向下流入兩板中間，傳向堤前方。

兩板之間及下層板到水底之間的水質點基本做水平方向的往復運動。具體而言：從t/T=0至t/T=2/7時刻，這兩個區(qū)域的水體基本水平流向堤前方；從t/T=3/7至t/T=5/7時刻，水質點改變流動方向，水體基本水平流向堤后方；從t/T=6/7至t/T=1時刻，水質點再次向堤前方做水平流動。

3.2 防波堤附近的渦量場變化過程

圖3給出了防波堤附近的渦量場分布在一個波浪周期內不同相位時刻的計算結果示例。其中渦量采用下式計算

圖4中可以看出，水平板周圍水域最大渦強度主要分布在上層板和自由水面之間。一個完整周期范圍內，最大渦強位置隨波浪傳播向前移動。

圖3 防波堤附近的渦量場變化過程（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）Fig.3 Vorticity field around the breakwater（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）

在t/T=0時刻，上層板上表面的前端一個負渦發(fā)生脫落，隨著時間的發(fā)展，這個負渦逐漸向尾端脫落；直至t/T=5/7時刻，負渦脫落至上層板上表面的尾端。而在t/T=6/7和t/T=1時刻，尾端的流場發(fā)生轉向（見圖2），一個正渦在尾端發(fā)生脫落。可見，渦隨著流場的變化而變化。在負渦從前端向尾端脫落的同時，還可以看出，從t/T=1/7到t/T=1時刻，在自由表面附近一個正向的渦從板首向板尾發(fā)生脫落，在這過程中渦的強度和面積越來越大，直到t/ T=1時刻，波浪傳至防波堤后方，渦的強度和面積開始減小。

圖4 剖面劃分示意圖Fig.4 Sketch of longitudinal section

以上可見，從板首到板尾，正渦和負渦一前一后交替發(fā)生脫落。

從圖2還可以看出，由于雙層板的存在，波浪對兩板之間和下層板下方的水體運動擾動很小，水質點基本在水平方向做往復運動，因此在這個區(qū)域內渦量很小。

3.3不同位置的波能流

為了探討消浪機理，進一步研究波浪能量在雙層水平板防波堤附近的傳遞過程是必要的。在此計算波能流沿程的變化，該變化對了解防波堤的消浪機理有重要的作用。

圖5 波能流沿程變化（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）Fig.5 Variations of the wave energy flux along the breakwater（H=0.06 m，D=0.48 m，T=1.4 s，B=0.9 m，S=10.5 cm）

首先沿水平板從水底到自由水面等間距做垂向剖面，剖面編號從1～36，每個剖面又分為三部分，表示為a-c。a部分從自由水面到上層板上表面，b部分為雙層板之間，c部分從水底到下層板下表面，詳細如圖4所示。

為方便表述，做以下定義：Pna，Pnb，Pnc分別表示通過每個剖面a部分，b部分，c部分的波能流，上標n為剖面編號。PnT表示Pna，Pnb，Pnc之和，也即通過整個剖面的總波能流。波能流P可由公式（2）計算

式中：p為波壓力；u為水質點的水平速度。

無雙層板時，從水底到自由水面的總波能流PwD可通過下式計算

式中：c為波浪傳播的相速度；k為波數。

無雙層板時，從水底到任意深度的波能流Ph可由下式計算

因此無板時通過每個剖面a部分，b部分，c部分的波能流顯然就都可求出，并分別用Pwa，Pwb，Pwc表示。

圖5給出了一個周期內，通過潛式堤的每個剖面a、b、c三部分的波能流以及總波能流。橫坐標為無量綱的距離x/B，即剖面到板首的距離與板寬的比值；縱坐標為波能流，并用無板時對應的波能流無量綱化。圖中可以發(fā)現，通過剖面b部分和c部分的波能流均為負值，也就是說在一個波浪周期內，波浪能量不但沒能從雙層板之間和下層板下方傳到堤后，而且還會有堤后能量從這兩個區(qū)域向堤前方傳遞，說明這兩個區(qū)域對波浪有阻塞作用。通過板首剖面a部分的波能流比無板時的波能流要大很多，這是因為b部分和c部分有阻塞作用，波能量被迫從a部分向堤后傳遞。對比圖2，從t/T=1/7到t/T=3/7時刻，由于上層板下方的回流作用，部分波浪能量從板首流入上層板的上方，導致板上方的波浪增大，因此會發(fā)生破碎，再加上紊動耗散以及a部分自身亦會反射部分波浪能量，所以板尾處的波能流比板首處的波能流減小了很多。通過板尾剖面的總波能流與無板時的總波能流相比已經明顯減小，因此防波堤起到了很好的消浪作用。

4結論

本文通過數值模擬的方法，研究了潛式雙層水平板型防波堤附件的流場、渦量場變化過程和不同位置的波能流。結果表明：（1）防波堤雙層板之間的流場呈水平往復流動狀態(tài)，下層板到水底之間的流動也接近水平往復流。（2）防波堤上層板上表面和自由水面之間，從板首到板尾，正渦和負渦一前一后交替發(fā)生脫落。渦的脫落隨著流場的變化而變化。最大渦強隨波浪傳播而移動。（3）防波堤上層板上方水域的渦動能量損耗及雙層板之間和下層板下方波能流的逆向傳遞是防波堤的主要消浪機理。

［1］Heins A E.Water waves over a channel of finite depth with a submerged plane barrier［J］.Canadian Journal of Mathematics?Journal Canadien De Mathematiques，1950，2（2）：210-222.

［2］Patarapanich M.Forces and moment on a horizontal plate due to wave scattering［J］.Coastal Engineering，1984，8（3）：279-301.

［3］Cheong H F，Patarapanich M.Reflection and transmission of random waves by a horizontal double?plate breakwater［J］.Ocean Engi?neering，1992，18（1-2）：63-82.

［4］Ushaa R，Gayathri T.Wave motion over a twin?plate breakwater［J］.Ocean Engineering，2005，32（8-9）：1 054-1 072.

［5］王國玉.多層水平板透空式防波堤消浪性能試驗研究［J］.大連理工大學學報，2005，45（6）：865-870. WANG G Y.Experimental study of wave?damping performance of multiple layer breakwater［J］.Journal of Dalian University of Technology，2005，45（6）：865-870.

［6］Neelamani S，Gayathri T.Wave interaction with twin plate wave barrier［J］.Ocean Engineering，2006，33（3-4）：495–516.

［7］LI Jing?bo，ZHANG Ning?chuan，GUO Chuan?sheng.Numerical simulation of waves interacting with a submerged horizontal twin?plate breakwater［J］.China Ocean Engineering，2010，24（4）：627-640.

［8］李靖波，張寧川，劉愛珍.潛式雙層水平板型防波堤的數值研究［J］.水道港口，2014，35（4）：317-324. LI J B，ZHANG N C，LIU A Z.Numerical investigation of the submerged horizontal twin?plate breakwater［J］.Journal of Waterway and Harbor 2014，35（4）：317-324.

Analysis of wave dissipation mechanism of submerged horizontal twin?plate breakwater

LI Jing?bo1，ZHANG Ning?chuan2，LIU Ai?zhen3
（1.Dalian Yuan?kun Marine Science and Technology Co.,Ltd.,Dalian 116000,China;2.State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian 116024,China;3.Tianjin Research Institute for Water Transport Engineering,Tianjin 300456,China）

Horizontal twin?plate breakwater is a relatively new type of breakwater.At the present,there is little literature available on the velocity field,vorticity field and wave dissipation mechanism of the breakwater.In this pa?per,wave interaction with a submerged horizontal twin?plate breakwater was investigated numerically.To address this issue,the continuity equation,Reynolds time?averaged equations and the two equationk-εmodel were used as the governing equations.Then,a numerical wave model was established.The velocity field and the vorticity field around the breakwater and the wave energy flux along the plate were numerically investigated.The computed re?sults reveal that the wave energy flux between the two plates can transfer to the opposite direction.And the wave en?ergy flux beneath the lower plate has the same transfer mode.The wave dissipation mechanism is the energy loss due to the vortex and the reversion transmission of the wave energy flux.

twin?plate breakwater;velocity field;vortex;wave energy flux;wave dissipation mechanism

TV 143；O 242.1

A

1005-8443（2016）03-0255-05

2015-07-03；

2015-08-03

國家自然科學基金支持（51221961）

李靖波（1979-），男，吉林洮南人，博士，工程師，主要從事波浪與建筑物的相互作用研究工作。

Biography：LI Jing?bo（1979-），male，engineer.