基于測波雷達(dá)現(xiàn)實(shí)方向譜的三維波浪數(shù)值模擬研究

王作超，石愛國，吳 明，劉 猛

(海軍大連艦艇學(xué)院，遼寧 大連116018)

0 引 言

傳統(tǒng)海浪環(huán)境的描述，一般是從頻率和波浪能量關(guān)系的角度，采用海浪頻譜來進(jìn)行研究。這樣的海浪頻譜是假設(shè)海浪能量只沿一個方向傳播，所以一般也被認(rèn)作單方向海浪譜。

實(shí)際上，波浪沿多個方向傳播。因此要更加準(zhǔn)確地描述波浪環(huán)境，必須把波浪能量的方向分布情況描述出來，由此產(chǎn)生了方向譜的概念。海浪的方向譜描述了波浪能量隨頻率和方向的分布情況。顯然，在海運(yùn)工程、船舶設(shè)計、船舶性能計算等領(lǐng)域，與海浪頻譜相比，方向譜可以生成更為接近真實(shí)的海洋環(huán)境，具有更大的優(yōu)勢。因此，自20世紀(jì)90年代起，有越來越多的學(xué)者從物理水池、數(shù)值水池(NWD)等方面開展了大量方向譜的研究工作［1］。近年來，隨著計算機(jī)計算能力的大幅提高，基于方向譜的數(shù)值水池研究更是發(fā)展迅速［2－3］。目前，方向譜研究領(lǐng)域存在的一個普遍問題是將方向譜一般看作頻譜和方向擴(kuò)散函數(shù)乘積的統(tǒng)計譜［4］，同真實(shí)的海浪環(huán)境仍有差別，因此利用真實(shí)波浪的觀測譜代替統(tǒng)計譜就成為深化數(shù)值波浪模擬研究中的必然課題。

在現(xiàn)有的各種方向譜觀測手段中，諸如遙感、水面、水下測波儀器，雖然類型不少，也各有特點(diǎn)，但可裝載于艦艇，且簡便易行，能實(shí)時獲取波浪信息，并具有較高觀測精度，目前唯有X 波段雷達(dá)測波系統(tǒng)。

自20世紀(jì)60年代，Oudshoorn 等［5］開始探索利用雷達(dá)獲取海浪信息以來，經(jīng)過近50年的發(fā)展，這一技術(shù)已趨于成熟。其中，德國GKSS 研發(fā)的WaMoSⅡ［6］，通過了德國船級社(GL)與挪威船級社(DNV)的認(rèn)證，并得到了廣泛安裝和應(yīng)用，具有代表性。為此，本文將WaMoSⅡ所采集的真實(shí)海浪信息代替?zhèn)鹘y(tǒng)的統(tǒng)計譜用于三維波浪的數(shù)值模擬，其成果將提高數(shù)值水池造波的仿真度，有助于提高船舶耐波性的預(yù)報水平。

1 測波雷達(dá)

雷達(dá)測波的系統(tǒng)組成、測量流程可以參考文獻(xiàn)［5，7－9］。測波雷達(dá)系統(tǒng)會受到障礙物、雨水、降雪等外界環(huán)境的影響，影響測量的準(zhǔn)確性。為了數(shù)據(jù)的可靠性，WaMoSⅡ在系統(tǒng)中引入數(shù)據(jù)質(zhì)量因子IQ(Quality Index)，以此判別所獲信息的準(zhǔn)確性，只有當(dāng)IQ 為000 時，才表示所獲數(shù)據(jù)可靠［10］。

本文所用的WaMoSⅡ測波設(shè)備，能測量海區(qū)波浪的頻譜和方向譜，主要測量參數(shù)如表1所示。

表1 WaMoSⅡ測量浪參數(shù)Tab.1 Measurement wave parameter of WaMoSⅡ

2 數(shù)值水池

短峰不規(guī)則波可視為無限多個頻率不等、方向不同、振幅變化且相位隨機(jī)的微幅簡諧波疊加而成。在空間位置(x，y)處，t 時刻的波面高度可表示為:

式中:aij為第i 個頻率、第j 個方向的成分波的波幅;ωi、ki分別為第i 個頻率的成分波的圓頻率和波數(shù);θj為第j 個方向的成分波的方向;εij為第i 個頻率、第j 個方向的成分波的隨機(jī)相位角，是［0，2π］區(qū)間內(nèi)均勻分布的隨機(jī)數(shù)。

其速度場的表達(dá)式為:

方向譜密度可分解為:

式中:S(ω)為頻譜;D(ω，θ)為歸一化的方向分布函數(shù)。一般的海浪研究和海浪數(shù)值模擬，S(ω)取統(tǒng)計譜，它通過長時間的觀測得到，反映了海浪譜的一般規(guī)律，同真實(shí)海浪差別較大。圖1 為WaMoSⅡ?qū)崪y的頻譜和以相同參數(shù)條件計算的ITTC 雙參數(shù)譜對比圖。從中可以看出二者的明顯差別。就這一具體海區(qū)的海浪狀況而言，實(shí)測譜的能量更高，也更為集中。

圖1 WaMoSⅡ頻譜和ITTC 譜頻對比圖Fig.1 The comparison of spectrum between WaMoSⅡand ITTC

據(jù)此，本文直接采用WaMoSⅡ的實(shí)測數(shù)據(jù)，從其所給的方向譜數(shù)據(jù)文件中提取方向、頻率和對應(yīng)的譜密度值，以求營造更真實(shí)的海浪環(huán)境。

本文數(shù)值模擬中的控制方程及其離散、湍流模型可參考文獻(xiàn)［2］。采用邊界造波法生成波浪［11］，采用阻尼消波法實(shí)現(xiàn)消波［12］。

3 數(shù)值計算

3.1 數(shù)據(jù)描述

根據(jù)上文提到的WaMoSⅡ有關(guān)數(shù)據(jù)的質(zhì)量要求，本文選取了“海洋22”號于2009年9月13日1 000 時刻觀測的方向譜信息。其有義波高為3.8 m，平均周期9.6 s，主波向67 °，IQ 為000。WamosⅡ提供的方向目標(biāo)譜如圖2所示。

圖2 WaMoSⅡ提供的方向譜(Date:20090913 Time:1000)Fig.2 The directional spectrum given by WaMoSⅡ(Date:20090913 Time:1000)

從圖2 中可以看出，所選樣本數(shù)據(jù)的能量主要集中在40 ° ～110 °之間，其他角度的能量可忽略不計。

3.2 求解絕對譜密度函數(shù)

當(dāng)利用海浪譜進(jìn)行數(shù)值波浪模擬時，其需要具有絕對譜密度值的方向譜。而WaMoSⅡ給出的是幅值歸一化的相對方向譜，不能滿足要求。本文利用譜形相似條件，將相對方向譜密度解算為絕對方向譜密度。求解過程如下:

1)采用辛普森法，求得相對方向譜的譜體積

2)根據(jù)測得的H1/3，由求得真實(shí)的海浪方差mt。

3)計算譜值比例系數(shù):

4)求解絕對譜密度函數(shù):

根據(jù)上述步驟，解得所采數(shù)據(jù)的mf=41.137 6，k=0.021 8，完成相對譜密度值向絕對譜密度值的轉(zhuǎn)換。

3.3 計算設(shè)置

各項(xiàng)設(shè)置如表2所示。

表2 計算設(shè)置Tab.2 The setting up of numerical simulation

3.4 邊界條件設(shè)置

1)造波區(qū)入口邊界——速度入口，給定波浪沿3 個方向的分速度，根據(jù)波高給出水的體積分?jǐn)?shù)及k，ε 值;

2)消波區(qū)入口邊界——速度入口，給定流速沿3 個方向的分速度，根據(jù)靜水面位置給出水的體積分?jǐn)?shù)及k，ε 值;

3)上下邊界——速度入口，給定3 個方向流速及水的體積分?jǐn)?shù);

4)出口邊界——壓力出口，設(shè)定靜壓力、底面位置和自由面高度，回流的湍動能k，耗散率ε。

圖3 計算域及邊界設(shè)置Fig.3 Grid domain and setting up of boundary layer

4 數(shù)值模擬結(jié)果及分析

4.1 頻譜分析

在短峰不規(guī)則波的模擬中，采用虛擬5 點(diǎn)浪高儀陣列(見圖4)監(jiān)測波面時歷。圖5 中給出的是所模擬的短峰波瞬時波面圖。圖6 是5 點(diǎn)的部分波高時歷。

根據(jù)所采集的波浪時歷，經(jīng)FFT 得到頻譜，同時與WaMoSⅡ給出的目標(biāo)頻譜和同參數(shù)條件下求解的ITTC 譜進(jìn)行比較(見圖7)。由圖中可以看出，相較于目標(biāo)頻譜，通過數(shù)值模擬所獲的頻譜能量雖有所衰減，但譜形吻合程度較高。為抑制衰減，可采取網(wǎng)格加密、改變網(wǎng)格分布及尺寸等措施。

圖4 浪高儀陣列Fig.4 The probe array of surface elevation

圖5 短峰波的瞬時波面圖Fig.5 The animation of surface of short crested wave

圖6 浪高儀陣列波面時歷Fig.6 The time history of surface elevation of probe array

圖7 頻譜對比圖Fig.7 The comparison of spectrum of NWD，ITTC and WaMosⅡdata

4.2 方向譜分析

方向譜主要的分析方法有:直接傅里葉變換方法(DFTM)、參數(shù)法(PM)、最大似然法(MLM)、擴(kuò)展的最大似然法(EMLM)、最大熵法(MEP)和貝葉斯方法(BDM)［13－14］。上述方法中，BDM 法具有較高的分辨率，當(dāng)陣列儀器數(shù)不少于4 時就能給出符合實(shí)際的方向分布，受儀器間距、方向分布寬度以及互譜誤差影響最小［15］。故本文采用BDM 法對所獲的波面時歷進(jìn)行方向譜分析，給出波能隨方向的分布情況(見圖8和圖9)，獲得的三一波高、平均周期、主浪向等參數(shù)如表3所示。

圖8 數(shù)值模擬的三維波方向譜Fig.8 The 3－directional spectrum of numerical simulation

圖9 數(shù)值模擬的方向譜等值線圖Fig.9 The contour of directional spectrum of numerical simulation

表3 所獲參數(shù)與目標(biāo)值對比Tab.3 The comparison of parameters between numerical simulation and WaMosⅡ

從模擬結(jié)果看，與WaMoSⅡ的目標(biāo)方向譜譜形吻合較好，除有義波高外，波浪的平均周期和主浪向與目標(biāo)值均較為接近，模擬精度較高。考慮到計算效率和模擬精度，在造波過程中，對WaMoSⅡ給出的原始數(shù)據(jù)，刪除了部分對譜密度值影響較低的角度和頻率，這是導(dǎo)致方向譜、頻譜波動的總能量同目標(biāo)能量有一定偏差的原因之一。

5 結(jié) 語

本文基于WaMoSⅡ給出的方向譜數(shù)據(jù)開展的三維波浪數(shù)值模擬，達(dá)到了較高模擬精度，與傳統(tǒng)的統(tǒng)計譜相比，能更好地模擬真實(shí)海浪，為后續(xù)船舶耐波性的研究提供了更加真實(shí)的海浪環(huán)境。

本文的三維波浪數(shù)值模擬，受限于現(xiàn)有的計算資源和計算能力，在權(quán)衡計算效率與計算精度后，在網(wǎng)格、頻率、角度的選取上做了一定的舍取，對波浪模擬的準(zhǔn)確性有一定影響。隨著計算機(jī)硬件技術(shù)的發(fā)展與波浪模擬算法的改進(jìn)，對真實(shí)海浪環(huán)境的數(shù)值模擬必將精度更高、耗費(fèi)時間更短，更加實(shí)用化，用于海上作業(yè)實(shí)踐。

需要指出的是:本文的意義，在于為在航船獲取實(shí)時的真實(shí)海浪環(huán)境，實(shí)現(xiàn)風(fēng)浪航行優(yōu)化決策的創(chuàng)新作業(yè)模式，提供了技術(shù)支撐。

［1］ALLENDER J，et al．The WADIC Project:a comprehensive field evaluation of directional wave instrumentation［J］．Ocean Enginering，1989，16:505－536．

［2］吳乘勝，朱德祥，顧民．?dāng)?shù)值波浪水池及頂浪中船舶水動力計算［J］．船舶力學(xué)，2008，12(2):171－179．

WU Cheng-sheng，ZHU De-xiang，GU Min．Computation of hydrodynamic forces for a ship in regular heading waves by a viscous numerical wave tank［J］．Journal of Ship Mechanics，2008，12(2):171－179．

［3］馮光，吳乘勝，鄭文濤，等．?dāng)?shù)值水池短峰不規(guī)則波模擬研究［J］．船舶力學(xué)，2010，14(4):347－354．

FENG Guang，WU Cheng-sheng，ZHENG Wen-tao，et al．Reproduction of short-crested irregular waves by a viscous numerical wave tank［J］．Journal of Ship Mechanics，2010，14(4):347－354．

［4］PARK J C，et al．Numerical reproduction of fully nonlinear multi－directional waves by a viscous 3D numerical wave tank［J］．Ocean Engineering，2004，31(11):1549－1565．

［5］OUDSHOORN H M．The use of radar in hydro-dynamic surveying［J］．Coastal Engineering，1960(7):59－76．

［6］Reference list［EB/OL］．http://www．oceanwaves．org /start．html．

［7］VALENZUELA G R．Theories for the interaction of electromagnetic and oceanic waves a review［J］．Boundary-Layer Meteorology，1978(19):87－105．

［8］YOUNG I R，ROSENTHAL W，ZIEMER F．A threedimensional analysis of marine radar images for the determination of ocean wave directionality and surface currents［J］．Journal of Geophysical Research，1985，90:1049－1059．

［9］ZIEMER F，GUNTHER H．A system to monitor ocean wave fields［C］．Proceedings of the Second International Conference on Air－Sea Interaction and Meteorology and Oceanography of the Coastal Zone，Lisbon，1994．

［10］WaMos Ⅱ Wave Monitoring System Operating Manual Installation Guide，2006，44．

［11］劉霞，譚國煥，王大國．基于邊界造波法的二階Stokes 波的數(shù)值生成［J］．遼寧工程技術(shù)大學(xué)學(xué)報(自然科學(xué)版)，2010，29(1):107－111．

LIU Xia，TAN Guo-huan，WANG Da-guo．Numerical simulation of second-order Stokes based on wave-generation method of defining inlet boundary conditions［J］．Journal of Liaoning Technical University(Natural Science)，2010，29(1):107－111．

［12］李凌，林兆偉，尤云祥，等．基于動量源方法的黏性流數(shù)值波浪水槽［J］．水動力學(xué)研究與進(jìn)展(A 輯)，2007，22(1):76－82．

LI Ling，LIN Zhao-wei，YOU Yun-xiang，et al．The numerical wave flume of the viscous fluid based on the momentum source method［J］．Journal of Hydrodynamics(Ser．A)，2007，22(1):76－82．

［13］HASHIMOTO N，NAGAI T，ASAI T．Extension of maximum entropy principle method(MEP)for estimating directional wave spectrum，24thInt．Conf．Coastal Eng，1994，11:232－246．

［14］HASHIMOTO N，KDBUNE K，KAMBYAMA Y．Estimation of directional spectrum using the Bayesian approach，and its application to field data analysis［R］．Report of the Port and Harbor Research Institute，1987(26):57－100．

［15］柳淑學(xué)，俞聿修．方向譜的分析方法［J］．港口工程，1994(1):25－32．

LIU Shu-xue，YU Yu-xiu．The analytical methods of directional spectrum．Port Engineering［J］．1994(1):25－32．