丁宏卿 張 翔
(河北大學建筑工程學院,河北保定 071002)
圓弧面防波堤是在半圓形防波堤應用經(jīng)驗基礎上開發(fā)的一種新型防波堤。
目前研究得比較多的是計算波浪力的經(jīng)驗公式,數(shù)值模擬方面也做了一定的工作。本文用NEWFLUME數(shù)值水槽模擬波浪作用于防波堤的全過程。
由于破碎波復雜的流動和紊流結構,現(xiàn)在這一現(xiàn)象的模擬已經(jīng)成為眾多海岸研究者的熱點研究課題。Lin和Liu(1998)提出了一個二維數(shù)值模型,基于雷諾平均的NS方程求解平均流動,采用改進的k—ε方程模擬紊流輸移,從而封閉NS方程。模型已經(jīng)通過實驗進行了驗證。模型有了多方面改進,如增加了多孔介質(zhì)內(nèi)流動的模擬功能,增加了內(nèi)部源函數(shù)造波功能等。
對于波源,有邊界造波和內(nèi)部源函數(shù)造波兩種方法。
邊界造波法:由于大多數(shù)周期波浪理論都會出現(xiàn)質(zhì)量凈輸移現(xiàn)象,而邊界造波法指定自由表面位移和速度值,因此會有凈質(zhì)量流進入計算域,累計后會導致計算溢出。
源函數(shù)造波法:采用源函數(shù)域內(nèi)造波法的優(yōu)勢在于不用指定邊界處的波浪條件,避免邊界波浪反射問題,并且由于波浪在域內(nèi)產(chǎn)生,因此左右邊界都可以設置為輻射邊界條件,很容易就可以吸收外傳的波列。
本文采用內(nèi)部源函數(shù)造波,為了保證波浪作用防波堤的穩(wěn)定,源函數(shù)造波區(qū)域離防波堤的距離為3倍~5倍的波長。
實驗采用1/4圓弧沉箱防波堤進行不規(guī)則同步壓力試驗。模型按重力相似準則,幾何比尺1∶36。
波浪要素分組如表1所示。

表1 波浪要素
實驗模型是按重力相似準則設計,幾何比尺1∶36。根據(jù)重力相似準則,試驗采用的波高為實際波高的1/36,而試驗周期為實際周期的1/6。得到的試驗波高和周期見表2。

表2 試驗的波高和周期
斷面圖見圖1。

本研究重在計算堤體上的波浪力。圖2~圖4是數(shù)值模型指定位置的位移、速度、壓力的時間歷程曲線圖。



由圖2~圖4可以得出,水平速度與自由表面位移具有相同的相位,垂向速度則有90°的相位差。
通過NEWFLUME數(shù)值水槽模擬波浪作用于防波堤的全過程,可以得出以下結論:
1)NEWFLUME數(shù)值水槽??梢詫Ψ啦ǖ踢M行詳細分析,可以得到質(zhì)量、動能、勢能等一系列數(shù)據(jù),可應用到實際工程中去。2)NEWFLUME數(shù)值水槽模在利用輸出文件生成圖形過程中,會出現(xiàn)一些額外的誤差。
[1]Tanimoto K,Takahashi S.Japanese Experiences on Composite Breakwaters[C].Proc.Intern.Workshop on Wave Barriers in Deepwaters.Yokosuda,1994:12-21.
[2]Xie Shileng.Wave Forces on Submerged Semi-Circular Breakwater and Similar Structures[J].China Ocean Engineering,1999,13(1):63-72.
[3]俞聿修,張寧川,饒永紅.半圓型防波堤的水力特性研究[J].海洋工程,1999,17(4):39-48.
[4]袁德奎,陶建華.半圓型防波堤波浪力的計算方法[J].中國港灣建設,2002(2):11-15.
[5]李炎保,馬強華,谷漢斌.半圓型防波堤波浪力計算方法的討論[J].港工技術,2003(6):2.
[6]謝世楞,李炎保,吳永強,等.圓弧面防波堤波浪力初步研究[J].海洋工程,2006,24(1):15-16.
[7]Jiang Xue-lian,Gu Han-bin,Li Yan-bao.Numberical Simulation on Hydraulic Performances of Quarter Circular Breakwater[J].China Ocean Engineering,2008,22(4):585-594.
[8]郄祿文,秦一楠.1/4圓弧面胸墻沉箱防波堤波浪力實驗研究[J].港工技術,2010,12(6):23-24.