劉師輝　胡金鵬

摘 要：該文基于不可壓縮黏性流體的N-S方程和RNG 湍流模型，利用Fluent軟件的二次開發(fā)功能，采用VOF方法建立了二維波浪數(shù)值水槽。通過UDF定義動邊界的運動規(guī)律以及在水槽尾端添加消波區(qū)的方式，實現(xiàn)了二維線性規(guī)則波和二階stokes波的造波和消波，得到的數(shù)值波浪形狀和理論波形擬合良好。參照已有物理模型試驗資料，數(shù)值模擬了波浪對水平板的沖擊過程，對比物理試驗和數(shù)值試驗板下11個測壓點峰值壓強數(shù)據(jù)，發(fā)現(xiàn)數(shù)值試驗數(shù)據(jù)較好的擬合了物理試驗數(shù)據(jù)，表明Fluent在模擬波浪與結(jié)構(gòu)物的相互作用方面具有較高的可靠性，為進一步探究影響波浪對水平板沖擊壓強大小的因素奠定了基礎(chǔ)。

關(guān)鍵詞：Fluent；數(shù)值造波；水平板；沖擊壓強

中圖分類號：P731.22 文獻標識碼：A

Abstract： Based on the N-S equations for viscous， incompressible fluid and RNG turbulence model， a two-dimensional numerical wave tank is established by utilizing the secondary development function of the software Fluent and the VOF method. A two-dimensional linear regular wave and a second-order stokes wave are generated by defining the motion of moving boundary and adding wave absorption zone at the end of the wave tank. The obtained numerical wave shape agreed well with the theoretical shape. The impact of waves on the horizontal plate is simulated referring to the existing physical model test. By Comparison of the physical and the numerical pressure data of 11 pressure points under the horizontal plate， it is found that they fit well， indicating the effectiveness of Fluent in simulating the interaction between waves and structures， which lays a foundation for further exploring the factors that affect the magnitude of the impact pressure on the horizontal plate.

Key words： Fluent； Numerical waves； Horizontal plate； Impact pressure

1 前言

隨著海洋經(jīng)濟的飛速發(fā)展，探究波浪的生成、演變及其對海工結(jié)構(gòu)物作用的規(guī)律成為了一個熱門學科。為了更好的研究波浪對海工結(jié)構(gòu)物的作用的規(guī)律性，在研究中常采用價格低廉且可重復的人工模擬造波的方式進行。人工模擬造波分為物理模型試驗和計算機數(shù)值模擬試驗兩種方法：相較于物理模型試驗，計算機數(shù)值模擬試驗的優(yōu)點主要是具有較高的經(jīng)濟性和較小的比尺效應，并且可以有效避免物理模型試驗中傳感器對流場的影響，因此計算機數(shù)值模擬造波得到越來越多學者的青睞。

目前關(guān)于數(shù)值造波的研究多集中于造波和消波方法上。

梁修峰[1]仿物理造波方法中的揺板式造波，結(jié)合Fluent軟件中用戶自定義UDF及動網(wǎng)格功能，較好的模擬了不規(guī)則波的運動；Kim[2]等利用改進的MAC方法，實現(xiàn)了三維數(shù)值水槽中不規(guī)則波浪的模擬；劉加海[3]通過設(shè)定造波板邊界的運動規(guī)律來模擬物理造波板的運動過程，實現(xiàn)了二維規(guī)則波的數(shù)值模擬，并進一步分析了模擬波高隨造波板運動周期、振幅和坐標的關(guān)系；李宏偉[4]同時采用揺板式造波和源造波兩種方法，模擬生成了線性規(guī)則波、不規(guī)則波、stokes波以及孤立波，并對比了兩種方法的優(yōu)缺點；辛穎[5]利用UDF宏功能定義了板式造波機的運動規(guī)律并實現(xiàn)了消波區(qū)的消波功能，所模擬的數(shù)值波浪與理論解及邊界元結(jié)果吻合良好。

在消波技術(shù)方面，董志[6]提出了多孔介質(zhì)消波方法，認為多孔介質(zhì)消波效果良好并給出了最佳的造波-消波組合；韓鵬[7]基于VOF方法，在水槽右端設(shè)置一定長度的海綿層吸收波浪，對比了五種不同海綿層衰減系數(shù)的水槽反射率，認為線性衰減系數(shù)的消波效果最好。

在數(shù)值水槽中加入結(jié)構(gòu)物以模擬現(xiàn)實海洋環(huán)境下波浪對海洋建筑物的作用是當下十分流行的趨勢。鄭艷娜[8]利用Fluent軟件的二次開發(fā)功能，建立了二維波浪數(shù)值水槽，通過在波浪場設(shè)置潛堤，探究了Fluent在模擬波浪與結(jié)構(gòu)物相互作用方面的有效性；Rolf Baarholm[9]采用非線性邊界元方法模擬了波浪沖擊水平板的過程；孫家文[10]應用改進的不可壓縮流體動力學（ISPH）方法，成功建立了波浪對水平板沖擊作用的數(shù)值模型，模型結(jié)果與物理試驗結(jié)果吻合良好。

本文基于Fluent軟件用戶自定義函數(shù)（UDF）功能，采用推板式造波方法模擬生成二維線性規(guī)則波及二階stokes波。在數(shù)值水槽添加水平板，模擬規(guī)則波對水平板的沖擊作用，通過對比數(shù)值試驗數(shù)據(jù)和物理試驗數(shù)據(jù)，以驗證Fluent在模擬波浪與結(jié)構(gòu)物作用上的可靠性。endprint

2 數(shù)學模型

流體連續(xù)性方程：

3 數(shù)值模擬

3.1 幾何模型

利用gambit軟件建立二維數(shù)值水槽（見圖1）：水槽長16 m、高1.2 m、水深0.8 m，消波區(qū)位于水槽最右端12～16 m處，造波板設(shè)定在左邊界；

流場網(wǎng)格劃分（見圖2）：網(wǎng)格總數(shù)16 800個；在非消波區(qū)域網(wǎng)格沿波浪傳播方向尺寸Δx=0.05 m；消波區(qū)采用漸變疏方式劃分，Δx=0.05～0.15 m；縱向網(wǎng)格采用指數(shù)網(wǎng)格劃分方式，在自由液面處加密處理以提高計算精度。

3.2 波浪的模擬

所模擬的規(guī)則波和二階stokes波要素相同：波高H=0.1 m、周期T=1.6 s、水深h=0.8 m、波長L=3.55 m；推板的運動規(guī)律通過UDF中的DEFINE_CG_MOTION宏來定義；消波源項則通過DEFINE_SOURCE宏來編寫，并最終在 Define命令窗口導入Fluent。

3.3 Fluent求解

模型采用推板式造波方法，故定義左邊界為動邊界、頂部邊界為壓力入口邊界，其余為固壁邊界；自由液面捕捉采用VOF方法；湍流模型使用RNG k-ε湍流模型，壓力速度項采用PISO算法進行迭代計算；動網(wǎng)格模塊選擇鋪層進行網(wǎng)格更新以實現(xiàn)造波板的運動，時間步長取0.005 s，總計算時長50 s。

3.4 計算結(jié)果分析

圖3 （a）、（b）分別為規(guī)則波1倍和3倍波長處虛擬浪高儀所測得的波浪位移隨時間變化歷程曲線與理論曲線的對比。從圖中可以發(fā)現(xiàn)：數(shù)值模擬波浪與理論解擬合的較好，各位置處模擬波浪周期與理論周期基本相同，相位也幾乎不存在偏差；但波高沿程有衰減，波高衰減程度隨波浪傳播距離增加而增加，這主要是因為基于勢流理論的規(guī)則波是在無旋無粘的理想狀態(tài)下提出的，而本文模擬的數(shù)值波浪考慮了水的粘性效應和紊動效應，因此波浪能沿程會損耗，而這也更符合現(xiàn)實波浪的傳播情況。

圖4 （a）、（b）分別為二階stokes波數(shù)值解與理論解對比曲線以及數(shù)值解與理論規(guī)則波波面對比曲線：從圖（a）中可以發(fā)現(xiàn)，該數(shù)值水槽模擬的二階stokes波波形曲線與理論波形曲線十分吻合，說明該模型模擬生成的stokes波效果良好；圖（b）中的波形曲線對比則很好的體現(xiàn)了stokes波“波峰尖陡，波谷平坦”的非線性特征。

4 波浪與水平板的相互作用

4.1 水平板及測壓點的布置

大連理工大學任冰[11]博士做過系列波浪對水平板的沖擊物理模型試驗，參照其規(guī)則波對水平板的沖擊作用物理模型試驗，波浪數(shù)值水槽幾何尺寸及平板布置見圖5：水槽長16 m、高1.0 m、水深0.6 m，消波區(qū)位于12～16 m處；水平板中心位于x=8 m處，板長1.0 m、板厚0.02 m，板底距自由水面0.02 m。

板底等間距布置11個虛擬測壓點，測壓點之間間距0.09 m，最左端和最右端測壓點距離板左右兩端點間距為0.05 m。

4.2 波浪對水平板的沖擊作用

根據(jù)相對板長的不同進行了兩組數(shù)值試驗：（1）H=0.1 m、T=1.5 s、Δh=0.03 m、L1=2.99 m、Lm/ L1 =0.33；（2）H=0.1 m、T=1.2 s、Δh=0.03 m、L2=2.12 m、Lm/ L2 =0.47。每組試驗模擬時間為30個波浪周期，每一個周期將有一個峰值壓強，將每個測點30個峰值壓強的前1/3大值的平均值PC，1/3作為該測點沖擊壓強特征值，并與物理模型試驗的壓強數(shù)據(jù)對比，其結(jié)果見圖6。

從圖6可以看出，不同工況條件下，數(shù)值模擬峰值壓強沿板底分布規(guī)律與物理模型試驗所得的規(guī)律大致相同，說明Fluent在模擬波浪對平板結(jié)構(gòu)物的沖擊作用方面具有較高的準確性。

5 結(jié)論

本文基于Fluent軟件平臺，利用其用戶自定義函數(shù)功能（UDF）仿物理推板式造波方法，成功模擬了二維線性規(guī)則波及二階stokes波，并通過在數(shù)值水槽中加入水平板式結(jié)構(gòu)物，模擬了波浪對水平板的沖擊過程?？偨Y(jié)如下：（1）該文采用RNG k-ε模型結(jié)合VOF自由液面捕捉方法模擬出的二維規(guī)則波及二階stokes波均很好的擬合了理論波形，stokes波較好的體現(xiàn)了“波峰尖陡，波谷平坦”的非線性特征，波高的沿程衰減在一定程度上反映了真實粘性流波浪的傳播狀態(tài)；（2）對照已有物理模型試驗，數(shù)值模擬規(guī)則波對水平板的沖擊過程，統(tǒng)計板下11個測點的特征峰值壓強并與物模數(shù)據(jù)對比，發(fā)現(xiàn)數(shù)值結(jié)果具有較高精度。不足之處在于數(shù)值試驗未考慮波浪的損耗以及水平板對波浪場的干擾等因素對計算值的影響，需后續(xù)研究進一步完善。

本文建立的波浪數(shù)值水槽為后續(xù)進一步研究波浪對板式結(jié)構(gòu)物沖擊作用的影響奠定了基礎(chǔ)。

參考文獻

[1] 梁修鋒， 楊建民， 李欣，等.不規(guī)則波浪的數(shù)值模擬[J].船舶力

學， 2010，14（5）.

[2] Kim M H， Celebi M S， Kim D J. Fully nonlinear interactions of waves with a

three-dimensional body in uniform currents[J]. Applied Ocean Research， 1998，

20（5）.

[3] 劉加海， 楊永全， 張洪雨，等. 二維數(shù)值水槽波浪生成過程及波浪形

態(tài)分析[J]. 四川大學學報（工程科學版）， 2004， 36（6）.

[4] 李宏偉. 數(shù)值水池造波方法研究[D]. 哈爾濱工程大學， 2009.

[5] 辛穎. Fluent UDF方法在數(shù)值波浪水槽中的應用研究[D]. 大連理工大

學， 2013.

[6] 董志， 詹杰民. 基于VOF方法的數(shù)值波浪水槽以及造波、消波方法研

究[J]. 水動力學研究與進展， 2009， 24（1）.

[7] 韓朋. 基于VOF方法的不規(guī)則波阻尼消波研究[D]. 大連理工大學，

2009.

[8] 鄭艷娜， 劉卓， 陳昌平，等. 基于Fluent軟件二維數(shù)值波浪水槽的研

究[J]. 中國海洋平臺， 2015， 30（6）.

[9] Baarholm R， Faltinsen O M. Wave impact underneath horizontal decks[J].

Journal of Marine Science & Technology， 2004， 9（1）.

[10] 孫家文. 波浪對水平板沖擊作用的實驗研究與數(shù)值模擬[D]. 大連理

工大學， 2011.

[11] 任冰. 隨機波浪對不同接岸型式碼頭上部結(jié)構(gòu)的沖擊作用研究[D].

大連理工大學， 2003.endprint