• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    可滲透島礁環(huán)境下的波浪運(yùn)動(dòng)研究

    2016-05-04 06:12:47陳彧超施鵬飛
    船舶力學(xué) 2016年8期
    關(guān)鍵詞:消波環(huán)礁瀉湖

    王 科,陳彧超,施鵬飛,邊 疆

    (大連理工大學(xué) 工業(yè)裝備與結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,工程力學(xué)系,遼寧 大連 116024)

    可滲透島礁環(huán)境下的波浪運(yùn)動(dòng)研究

    王 科,陳彧超,施鵬飛,邊 疆

    (大連理工大學(xué) 工業(yè)裝備與結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,工程力學(xué)系,遼寧 大連 116024)

    文章應(yīng)用邊界單元法,依據(jù)波浪繞射理論研究了潛淹沒(méi)島礁上的波浪運(yùn)動(dòng)問(wèn)題。研究中分析了島礁的滲透性、排布的數(shù)量、環(huán)礁瀉湖的深度等因素對(duì)于島礁消波作用的影響,獲得了不同島礁環(huán)境下波浪的透射系數(shù)、反射系數(shù)。研究發(fā)現(xiàn),可滲透島礁的消波效果比不可滲透島礁差,且孔隙率越大,島礁對(duì)于波浪的消波效果越??;單個(gè)島礁對(duì)長(zhǎng)波有明顯的消波作用,而對(duì)中短波幾乎沒(méi)有影響;但隨著島礁排布數(shù)量的增加,島礁群對(duì)于中波的消波效果增加;環(huán)礁的消波效果受到瀉湖深度的影響,在一定范圍內(nèi),瀉湖的深度越大,環(huán)礁對(duì)于中短波的消波效果越好,但當(dāng)瀉湖深度達(dá)到一定數(shù)值后,繼續(xù)增大深度,消波效果基本沒(méi)有提高。

    可滲透島礁;邊界單元法;波浪透射及反射系數(shù);瀉湖

    0 引 言

    我國(guó)的海域面積超過(guò)400萬(wàn)平方公里,南海地區(qū)有許多寬度達(dá)數(shù)公里、且潛深較小的珊瑚礁,這些自然島礁與潛式防波堤類似,有一定的消波能力。對(duì)于這些島礁,在開(kāi)發(fā)其漁業(yè)及旅游價(jià)值的同時(shí),還可以充分利用其消波和護(hù)島功能,在周邊設(shè)置浮式結(jié)構(gòu)物,甚至構(gòu)建機(jī)場(chǎng)、補(bǔ)給基地等大型設(shè)施,來(lái)輔助礦產(chǎn)開(kāi)采、漁業(yè)生產(chǎn)、科研觀測(cè)等海洋開(kāi)發(fā)活動(dòng)[1]。要實(shí)現(xiàn)這些構(gòu)想,對(duì)于島礁環(huán)境下的波浪運(yùn)動(dòng)研究顯得非常重要。

    早期主要通過(guò)監(jiān)測(cè)及試驗(yàn)手段來(lái)研究島礁環(huán)境下的波浪運(yùn)動(dòng)問(wèn)題,Hardy等人[2]在1990年對(duì)澳大利亞大堡礁進(jìn)行可觀測(cè)試驗(yàn)來(lái)研究珊瑚礁與波浪的相互作用,結(jié)果表明島礁的消波作用很大程度上受其礁頂深度的影響。Gourlay(1994)[3],用模型實(shí)驗(yàn)研究了島礁環(huán)境下不同波浪條件和不同水深條件下的波浪透射特性。Lowe等人(2005)[4]在夏威夷卡內(nèi)奧赫灣的堡礁開(kāi)展了兩星期的實(shí)驗(yàn),布置多個(gè)測(cè)點(diǎn)來(lái)監(jiān)測(cè)波高、波速,研究了島礁對(duì)于波浪的消波作用。

    近年來(lái),考慮到島礁有著與潛式防波堤類似的特點(diǎn),一些水動(dòng)力計(jì)算方法開(kāi)始被用于島礁波浪的相互作用研究。Hearn等人(1999)[5]引入了深度系數(shù),用一維模型研究了島礁及瀉湖的水動(dòng)力特性,分析了深度系數(shù)對(duì)島礁水動(dòng)力特性的影響,并用澳大利亞寧格羅暗礁(Ningaloo reef)和美國(guó)卡內(nèi)奧赫灣(Kaneohe Bay)的相關(guān)數(shù)據(jù)檢驗(yàn)了模型的準(zhǔn)確性。丁軍等人(2014)[6],用Morison公式和Green函數(shù)法研究了島礁環(huán)境下半潛式平臺(tái)的運(yùn)動(dòng)響應(yīng)。田超等人(2014)[7-8]通過(guò)模型實(shí)驗(yàn),研究了島礁環(huán)境下波浪變形及浮式結(jié)構(gòu)的運(yùn)動(dòng)響應(yīng),并通過(guò)緩坡方程引入了半潛式平臺(tái)在島礁環(huán)境下的島礁影響系數(shù),提出了一種島礁環(huán)境下浮式結(jié)構(gòu)運(yùn)動(dòng)響應(yīng)的預(yù)報(bào)方法。Fang等人(2014)[9]用Boussinesq波浪模型并考慮波浪破碎及底部摩擦因素,對(duì)不同波浪條件及島礁環(huán)境下的波浪運(yùn)動(dòng)特性進(jìn)行了數(shù)值模擬,結(jié)果與監(jiān)測(cè)數(shù)據(jù)吻合較好。

    本文依據(jù)波浪繞射理論,發(fā)展了能求解任意復(fù)雜邊界的邊界元算法。通過(guò)算例驗(yàn)證,證明了本文采用的邊界元算法對(duì)于潛淹沒(méi)結(jié)構(gòu)上波浪運(yùn)動(dòng)問(wèn)題的精確性,并對(duì)潛淹沒(méi)島礁上的波浪運(yùn)動(dòng)問(wèn)題進(jìn)行了研究。本文計(jì)算了不同島礁環(huán)境下波浪的透射系數(shù)、反射系數(shù),分析了島礁的滲透性、排布的數(shù)量、環(huán)礁瀉湖的深度這幾個(gè)因素對(duì)于島礁消波作用的影響。

    1 數(shù)學(xué)理論及基本計(jì)算公式

    圖1 計(jì)算示意圖Fig.1 Calculation sketch

    1.1 邊界值問(wèn)題基本公式

    如圖1所示,一個(gè)拋物線形的島礁位于水底,水深H,島礁寬度為B=2a,其頂部距水面距離為HS。建立如圖所示的笛卡爾坐標(biāo)系,取原點(diǎn)o位于水面,波浪沿著x軸正方向入射??紤]不可壓縮、無(wú)粘性理想流體,則流速可以用速度勢(shì)函數(shù)的梯度表示,流體速度勢(shì)可以表示為:

    在自由表面SF上,速度勢(shì)滿足線性自由表面條件:

    若水底為不可滲透邊界,則在水底滿足:

    若物面S0可滲透,則可將流場(chǎng)分成兩個(gè)區(qū)域:外流域Φ1和可滲透結(jié)構(gòu)內(nèi)流域Φ2,兩個(gè)流域的速度勢(shì)都滿足拉普拉斯方程??紤]到相鄰流域必須在交界面上滿足連續(xù)性方程,則速度勢(shì)函數(shù)滿足:

    其中:ε表示可滲透結(jié)構(gòu)的孔隙率,S為慣性系數(shù),f為線性阻尼系數(shù)(Sollitt and Cross,1972),f可以通過(guò)下列公式求得[10]:

    其中:ν為運(yùn)動(dòng)粘度,Cf為湍流阻力系數(shù),KP為滲透介質(zhì)的固有滲透系數(shù),q表示滲流速度的實(shí)部。Cf和KP是滲透介質(zhì)的本身特性,需要由實(shí)驗(yàn)測(cè)得,f則需要迭代求解得到。

    計(jì)算邊界離島礁足夠遠(yuǎn)時(shí)(大于5倍水深),可認(rèn)為反射和透射波浪達(dá)到定常狀態(tài),速度勢(shì)滿足:

    基于線性假設(shè),可以將速度勢(shì)分解為入射勢(shì)φI與繞射勢(shì)φD:

    其中:φI為微幅波入射勢(shì),可表示為:

    其中:A為入射波振幅。

    1.2 邊界積分方程

    對(duì)于波浪繞射問(wèn)題,根據(jù)格林定理得到如下邊界積分方程:

    如取格林函數(shù)基本解為下式,則可采用鏡像方法消除水底邊界積分:

    當(dāng)物面邊界不可滲透時(shí),將邊界條件(4)、(5)、(6)和(10)代入(15)式,得到:

    1.3 透射系數(shù)與反射系數(shù)

    透、反射系數(shù)可以計(jì)算如下,定義:

    其中:eHR,eHI分別為eH的實(shí)部和虛部,相位α1為:

    類似地,定義:

    其中:eSR,eSI分別為eS的實(shí)部和虛部,相位α2為:

    反射系數(shù)RF和透射系數(shù)RT可表示為:

    1.4 可滲透島礁波浪運(yùn)動(dòng)研究

    圖2 可滲透結(jié)構(gòu)計(jì)算示意圖Fig.2 Sketch of permeable structure

    對(duì)于物面邊界不可滲透的情況,只存在一個(gè)流域,只需對(duì)外域的積分方程進(jìn)行求解,計(jì)算速度勢(shì)。而當(dāng)物面邊界為可滲透邊界時(shí),流場(chǎng)分為外部流域和滲透結(jié)構(gòu)的內(nèi)部流域,見(jiàn)圖2,物面邊界條件變?yōu)閮?nèi)外速度勢(shì)φ1,φ2及其導(dǎo)數(shù)的相互關(guān)系,此時(shí)需對(duì)內(nèi)外流域進(jìn)行迭代求解,來(lái)計(jì)算內(nèi)外速度勢(shì)。

    一般采用如下迭代步驟:

    (5)具體程序執(zhí)行框圖見(jiàn)圖3。

    圖3 迭代計(jì)算示意圖Fig.3 Sketch of iterative computation

    2 結(jié)果與分析

    2.1 算例驗(yàn)證

    為了驗(yàn)證本文采用邊界元算法的精確性,本文對(duì)潛淹沒(méi)矩形防波堤的波浪運(yùn)動(dòng)問(wèn)題進(jìn)行了研究,計(jì)算了其反射系數(shù),并將結(jié)果與Cho和Lee(2004)[11]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了比較。如圖4所示,考慮一組不可滲透的矩形結(jié)構(gòu),在水下依次排布,其尺寸為寬度la=高度lb=0.4 m,潛深Hs=0.4 m,矩形結(jié)構(gòu)的間距d為2.0 m,水深H為0.8 m,波高h(yuǎn)為0.04 m。

    圖4 潛沒(méi)式矩形防波堤的示意圖Fig.4 Sketch of submerged rectangular breakwater

    圖5 N=1時(shí)反射系數(shù)計(jì)算結(jié)果比較Fig.5 Comparison of reflection coefficient when N=1

    Cho和Lee(2004)[11]對(duì)此問(wèn)題進(jìn)行了實(shí)驗(yàn)研究,分別測(cè)得矩形個(gè)數(shù)N=1~3情況下的反射系數(shù),并與特征函數(shù)展開(kāi)法的計(jì)算結(jié)果進(jìn)行了比較。本文通過(guò)邊界元方法對(duì)此算例進(jìn)行計(jì)算,對(duì)于每個(gè)潛淹沒(méi)矩形結(jié)構(gòu)劃分180個(gè)邊界單元,得到反射系數(shù)結(jié)果如圖5~7所示,圖中橫坐標(biāo)為無(wú)因次化波數(shù)KH??梢钥吹剑吔缭ǖ挠?jì)算結(jié)果與Cho和Lee(2004)[11]的實(shí)驗(yàn)數(shù)據(jù)吻合較好,與特征函數(shù)展開(kāi)法的結(jié)果接近但也存在一定差異:邊界元法的結(jié)果曲線在前半段略小于特征函數(shù)展開(kāi)法,而在后半段則略大于特征函數(shù)展開(kāi)法,但是總體上兩者的計(jì)算結(jié)果比較接近。這說(shuō)明對(duì)于波浪與結(jié)構(gòu)物相互作用的問(wèn)題,邊界元方法有較高的計(jì)算精度。

    2.2 島礁的滲透性對(duì)于波浪運(yùn)動(dòng)的影響分析

    對(duì)于島礁計(jì)算模型,經(jīng)反復(fù)試算,將圖2中流域邊界沿水平方向劃分成840個(gè)邊界單元,垂直方向劃分140個(gè)邊界單元,將每個(gè)島礁邊界劃分為180個(gè)邊界單元進(jìn)行計(jì)算??蓾B透島礁潛沒(méi)于水面下,物面邊界S0為可滲透邊界,計(jì)算中具體取值如下:?jiǎn)蝹€(gè)島礁寬度B=2a=2.0 m,水深H=4.0 m,潛深HS=0.5 m,取湍流阻力系數(shù)Cf=0.295,固有滲透系數(shù)KP=1.057 2×10-7m2(Hsu和Wu,1999)[12],運(yùn)動(dòng)粘度ν=1.0126×10-6m2/s[13]。取平均孔隙率ε分別為0.1,0.2,0.3,0.4,0.5。下圖中橫坐標(biāo)為無(wú)因次化波數(shù)Ka,入射波浪振幅均為單位值。

    圖6 N=2時(shí)反射系數(shù)計(jì)算結(jié)果比較Fig.6 Comparison of reflection coefficient when N=2

    圖7 N=3時(shí)反射系數(shù)計(jì)算結(jié)果比較Fig.7 Comparison of reflection coefficient when N=3

    圖8 可滲透島礁的透射系數(shù)Fig.8 Transmission coefficient of a permeable reef with different porosity

    圖9 可滲透島礁的反射系數(shù)Fig.9 Reflection coefficient of a permeable reef with different porosity

    圖8~9中ε=0的透、反射系數(shù)曲線為不可滲透島礁的計(jì)算結(jié)果??梢钥吹剑瑵B透性并不影響曲線隨著相對(duì)波數(shù)Ka的增減變化趨勢(shì),以及峰值出現(xiàn)的相應(yīng)位置;但隨著平均孔隙率ε的增大,反射系數(shù)逐漸減小,透射系數(shù)逐漸增大。當(dāng)島礁不可滲透時(shí),反射系數(shù)在Ka=0.17處達(dá)到最大值,此時(shí)反射系數(shù)為0.63;而當(dāng)孔隙率ε達(dá)到0.5時(shí),反射系數(shù)的最大值僅為0.50。

    根據(jù)王新志等人(2008)[14]的實(shí)驗(yàn)結(jié)果,我國(guó)南沙群島的礁灰?guī)r孔隙率為45%~55%;本文的計(jì)算結(jié)果顯示,ε=0.5時(shí)的反射系數(shù)峰值約為不可滲透情況下80% ,可見(jiàn),礁灰?guī)r島礁的防浪作用要明顯差于不可滲透島礁。

    2.3 島礁排布數(shù)量對(duì)于波浪運(yùn)動(dòng)的影響

    考慮一組不可滲透的島礁,在水下連續(xù)排布,單個(gè)島礁寬度仍然取為B=2a=2.0 m,水深為H=4.0 m,潛深HS=1.0 m,取島礁間距d=B=2.0 m,計(jì)算島礁排布數(shù)量N分別為1,2,3時(shí)的相關(guān)參數(shù),圖10為多島礁的計(jì)算示意圖,計(jì)算結(jié)果如圖11~12所示。

    單個(gè)島礁的透、反射系數(shù)曲線呈拋物線狀,反射曲線的最大值約為0.45,透射曲線的最小值約為0.9;當(dāng)N=2時(shí),透、反射系數(shù)曲線出現(xiàn)兩處峰值:在Ka=0.05處,反射系數(shù)約為0.45,透射系數(shù)為0.89,在Ka=0.45處,反射系數(shù)達(dá)到0.5,透射系數(shù)減小到0.86;當(dāng)N=3時(shí),結(jié)果曲線的變化更為劇烈,有3處比較明顯的波峰:在Ka=0.02和Ka=0.2處,透射、反射系數(shù)達(dá)到峰值,反射系數(shù)為0.45,透射系數(shù)為0.9,在Ka=0.55處,反射系數(shù)達(dá)到最大值,約為0.62,透射系數(shù)則減小為0.78。

    隨著島礁排布數(shù)量N的增大,透射、反射系數(shù)曲線的變化愈發(fā)劇烈。以反射系數(shù)為例,隨著N的增大,出現(xiàn)了更多的波峰。雙島礁情況下產(chǎn)生的第一個(gè)波峰,和三個(gè)島礁情況下產(chǎn)生的前兩個(gè)波峰,其峰值和單個(gè)島礁情況下的最大值基本相同,約為0.45,且都出現(xiàn)在0<Ka<0.3的長(zhǎng)波范圍內(nèi),這表明島礁數(shù)量的增加不影響島礁對(duì)于長(zhǎng)波的作用,其對(duì)于長(zhǎng)波的消波效果和單個(gè)島礁基本相同;而在0.3<Ka<0.8的中波范圍內(nèi)產(chǎn)生的波峰,其峰值明顯大于其他波峰,且隨著島礁數(shù)量的增加而增大:N=2時(shí)反射系數(shù)為0.5,N=3時(shí),反射系數(shù)達(dá)到0.62,這表明島礁數(shù)量的增加能加強(qiáng)島礁對(duì)中波的消波作用。當(dāng)Ka>0.8時(shí),不同島礁對(duì)波浪的反射系數(shù)在0.2以下。

    圖10 多島礁的計(jì)算示意圖Fig.10 Sketch of multiple reefs

    圖11 不同數(shù)量島礁的透射系數(shù)Fig.11 Wave transmission coefficients of multiple reefs

    圖12 不同數(shù)量島礁的反射系數(shù)Fig.12 Wave reflection coefficients of multiple reefs

    圖13 環(huán)礁的計(jì)算示意圖Fig.13 Schematic sketch of lagoon reef

    綜上所述,當(dāng)潛深較大時(shí),單個(gè)島礁僅對(duì)長(zhǎng)波有一定的消波作用,但是對(duì)中短波幾乎沒(méi)有效果;島礁排布數(shù)量的增加,并不會(huì)改善島礁對(duì)于長(zhǎng)波的效果,但能夠顯著增強(qiáng)對(duì)于中波的阻礙作用。

    2.4 環(huán)礁的波浪運(yùn)動(dòng)分析

    我國(guó)南海的地區(qū)有許多環(huán)礁,例如華光礁、美濟(jì)礁、仙娥礁等。珊瑚蟲(chóng)圍繞島嶼繁殖形成環(huán)狀珊瑚,當(dāng)發(fā)生地質(zhì)變化,島嶼下沉至海面以下,就形成了環(huán)礁及圍繞其中的瀉湖。

    環(huán)礁的計(jì)算模型如圖13所示,環(huán)礁潛淹沒(méi)于水面以下,水深H=4.0 m,島礁寬度B=2a=4.0 m,瀉湖寬度為b=2.0 m,礁頂潛深HS=0.5 m,分別取瀉湖相對(duì)礁頂?shù)纳疃菻C為0.5 m,1.0 m,1.5 m,2.0 m,2.5 m和3.0 m,計(jì)算其透射系數(shù)和反射系數(shù)。

    從圖14~15可以看到和雙島礁排布的情況類似,環(huán)礁的透射系數(shù)和反射系數(shù)隨著相對(duì)寬度Ka的變化出現(xiàn)兩個(gè)明顯的波峰:第一個(gè)波峰出現(xiàn)在Ka=0.1處,第二個(gè)波峰出現(xiàn)在Ka=1.2處。以反射系數(shù)為例,瀉湖相對(duì)礁頂?shù)纳疃菻C變化對(duì)第一個(gè)波峰的峰值影響很小,隨著HC增大,第一個(gè)波峰的峰值幾乎不變,約為0.65。HC由0.5增大到1.5的過(guò)程中,第二個(gè)波峰的峰值開(kāi)始增大:HC=0.5時(shí),為0.47,HC=1.0時(shí)為0.51,HC=1.5時(shí),為0.52,而當(dāng)HC達(dá)到1.5后,瀉湖深度繼續(xù)增大對(duì)反射系數(shù)的影響很小,這時(shí)環(huán)礁的消波作用不再隨深度變大而增強(qiáng)。

    圖14 瀉湖深度不同時(shí)環(huán)礁的透射系數(shù)Fig.14 Transmission coefficient of lagoon reef with different lagoon depth

    圖15 瀉湖深度不同時(shí)環(huán)礁的反射系數(shù)Fig.15 Reflection coefficient of lagoon reef with different lagoon depth

    3 結(jié) 論

    本文根據(jù)波浪繞射理論,應(yīng)用邊界元方法,對(duì)于島礁環(huán)境下的波浪運(yùn)動(dòng)進(jìn)行了計(jì)算,分析了島礁滲透性,排布數(shù)量,環(huán)礁瀉湖的潛深等因素對(duì)于波浪運(yùn)動(dòng)的影響,得到以下結(jié)果:

    (1)可滲透島礁的消波效果差于不可滲透島礁,且孔隙率越大,島礁對(duì)于波浪的阻礙越小。

    (2)單個(gè)島礁的消波作用集中在長(zhǎng)波,而對(duì)中短波幾乎沒(méi)有效果;島礁排布數(shù)量的增加,會(huì)增強(qiáng)島礁對(duì)于中波的阻礙作用。

    (3)環(huán)礁的消波效果受到其瀉湖深度的影響,在一定范圍內(nèi),瀉湖深度越大,環(huán)礁對(duì)于中短波的消波效果越好;但當(dāng)深度達(dá)到一定程度后,繼續(xù)增大,則基本不影響其消波效果。

    [1]吳有生,田 超,宗 智,等.波浪環(huán)境下超大型浮式結(jié)構(gòu)物的水彈性響應(yīng)研究[C].第二十五屆全國(guó)水動(dòng)力學(xué)研討會(huì)暨第十二屆全國(guó)水動(dòng)力學(xué)學(xué)術(shù)會(huì)議文集(上冊(cè)),2013.

    [2]Hardy T A,Young I R,Nelson R C,et al.Wave attenuation on an offshore coral reef[J].Coastal Engineering Proceedings, 1990,1(22):330-344.

    [3]Gourlay M R.Wave transformation on a coral reef[J].Coastal Engineering,1994,23(1):17-42.

    [4]Fernandes L,Day J O N,Lewis A,et al.Establishing representative no-take areas in the great barrier reef:Large-scale implementation of theory on marine protected areas[J].Conservation Biology,2005,19(6):1733-1744.

    [5]Hearn C J.Wave-breaking hydrodynamics within coral reef systems and the effect of changing relative sea level[J].Journal of Geophysical Research:Oceans,1999,104(C12):30007-30019.

    [6]丁 軍,田 超,張 凱,等.近島礁半潛式平臺(tái)的運(yùn)動(dòng)及載荷響應(yīng)研究[C].第十三屆全國(guó)水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨第二十六屆全國(guó)水動(dòng)力學(xué)研討會(huì)論文集—F船舶與海洋工程流體力學(xué),2014.

    [7]田 超,丁 軍,吳有生.近島礁浮式平臺(tái)水動(dòng)力模型試驗(yàn)研究[C].第十三屆全國(guó)水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨第二十六屆全國(guó)水動(dòng)力學(xué)研討會(huì)論文集—D水動(dòng)力學(xué)實(shí)驗(yàn)和測(cè)試技術(shù),2014.

    [8]田 超,丁 軍,楊 鵬.島礁波浪環(huán)境下浮式結(jié)構(gòu)物的動(dòng)響應(yīng)預(yù)報(bào)[J].船舶力學(xué),2014,18(11):1284-1291. Tian Chao,Ding Jun,Yang Peng.Prediction of dynamic responses of floating structures under wave environment near islands and reefs[J].Journal of Ship Mechanics,2014,18(11):1284-1291.

    [9]Fang K,Yin J,Liu Z,et al.Revisiting study on Boussinesq modeling of wave transformation over various reef profiles[J]. Water Science and Engineering,2014,7(3):306-318.

    [10]Sollitt C K,Cross R H.Wave transmission through permeable breakwaters[J].Coastal Engineering Proceedings,1972,1(13): 1827-1846.

    [11]Cho Y,Lee J,Kim Y.Experimental study of strong reflection of regular water waves over submerged breakwaters in tandem [J].Ocean Engineering,2004,31(10):1325-1335.

    [12]Hsu H,Wu Y.Scattering of water wave by a submerged horizontal plate and a submerged permeable breakwater[J].Ocean Engineering,1998,26(4):325-341.

    [13]Wu Y T,Hsiao S C.Propagation of solitary waves over a submerged permeable breakwater[J].Coastal Engineering,2013, 81:1-18.

    [14]王新志,汪 稔,孟慶山,等.南沙群島珊瑚礁礁灰?guī)r力學(xué)特性研究[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(11):2221-2226. Wang Xinzhi,Wang Ren,Meng Qingshan,et al.Research on characteristics of coral reef calcareous rock in Nansha islands [J].Chinese Journal of Rock Mechanics and Engineering,2008,27(11):2221-2226.

    Study on wave motion over submerged reefs

    WANG Ke,CHEN Yu-chao,SHI Peng-fei,BIAN Jiang
    (State Key Laboratory of Structural Analysis for Industrial Equipment,Department of Engineering Mechanics, Dalian University of Technology,Dalian 116024,China)

    Based on wave diffraction theory,boundary element method is applied to study wave motion over submerged reefs.The transmission coefficient,reflection coefficient under different conditions such as permeability,number of arrays and the lagoon’s depth are obtained.It is found that,porous reefs have less wave elimination effect than impermeable reefs.When the porosity of reefs is bigger,their wave elimination ability will become worse.A single reef is only efficient on the wave elimination of long waves,while an array of reefs have an efficient wave elimination effect on medium waves.Lagoon reef’s wave elimination ability is determined by lagoon’s depth.The wave elimination effect of lagoon reef on medium and short waves will get better when lagoon’s depth increases.But the elimination effect will show little change when lagoon’s depth reaches a certain value.

    permeable reefs;boundary element method;wave transmission and reflection coefficient; lagoon

    U661.1

    A

    10.3969/j.issn.1007-7294.2016.08.005

    1007-7294(2016)08-0964-09

    2016-05-09

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2013CB036101);國(guó)家自然科學(xué)基金(51379037)

    王 科(1970-),男,副教授,E-mail:kwang@dlut.edu.cn;陳彧超(1990-),男,碩士研究生。

    猜你喜歡
    消波環(huán)礁瀉湖
    斜坡非淹沒(méi)柔性植物對(duì)船行波消減效能的試驗(yàn)研究
    循環(huán)水槽多層孔板消波裝置開(kāi)發(fā)及消波特性數(shù)值模擬
    密克羅尼西亞
    華夏地理(2018年3期)2018-08-06 03:59:24
    大洋洲的朗伊羅阿環(huán)礁
    印焦慮馬爾代夫向沙特賣(mài)礁
    圣約瑟夫環(huán)礁
    中外文摘(2016年11期)2016-08-31 19:52:50
    多層變孔徑傾斜孔板式消波裝置的試驗(yàn)研究
    英國(guó)潮汐瀉湖發(fā)電潛力評(píng)估
    斯旺西灣潮汐瀉湖項(xiàng)目取得進(jìn)展
    T型透空式防波堤消波性能的理論分析
    亚洲va在线va天堂va国产| 亚洲成人精品中文字幕电影| 给我免费播放毛片高清在线观看| 2021天堂中文幕一二区在线观| 国产精品永久免费网站| 嫩草影院新地址| 日本色播在线视频| 99热这里只有是精品50| 亚洲七黄色美女视频| 又粗又爽又猛毛片免费看| 中文字幕人妻熟人妻熟丝袜美| 免费人成视频x8x8入口观看| 国产熟女欧美一区二区| 91精品国产九色| 久久影院123| 亚洲欧美日韩无卡精品| 国产精品人妻久久久久久| 人妻系列 视频| 18禁裸乳无遮挡动漫免费视频| 亚洲在久久综合| 这个男人来自地球电影免费观看 | 美女主播在线视频| 亚洲精品,欧美精品| 久久影院123| 18禁裸乳无遮挡免费网站照片| 五月玫瑰六月丁香| 日日摸夜夜添夜夜爱| 久久99热这里只有精品18| 在线观看一区二区三区| 国产高清有码在线观看视频| 成人免费观看视频高清| 在线免费观看不下载黄p国产| 欧美高清性xxxxhd video| tube8黄色片| kizo精华| 午夜福利在线在线| 伦理电影免费视频| 狂野欧美激情性bbbbbb| 午夜免费观看性视频| 欧美最新免费一区二区三区| 成人国产麻豆网| 亚洲国产精品一区三区| 最近2019中文字幕mv第一页| 老女人水多毛片| 十八禁网站网址无遮挡 | 亚洲一区二区三区欧美精品| 午夜激情久久久久久久| 精品视频人人做人人爽| 中文字幕人妻熟人妻熟丝袜美| 亚洲久久久国产精品| 亚洲精品久久午夜乱码| 亚洲电影在线观看av| 亚洲欧美精品专区久久| 久久久久人妻精品一区果冻| 少妇 在线观看| 久久精品国产鲁丝片午夜精品| 久久99精品国语久久久| 亚洲精品中文字幕在线视频 | 美女视频免费永久观看网站| 久久韩国三级中文字幕| 美女cb高潮喷水在线观看| 国精品久久久久久国模美| 国国产精品蜜臀av免费| 亚洲电影在线观看av| h视频一区二区三区| 欧美另类一区| 特大巨黑吊av在线直播| 午夜老司机福利剧场| a 毛片基地| 欧美精品亚洲一区二区| 日本欧美视频一区| 九草在线视频观看| 亚洲av国产av综合av卡| 欧美精品国产亚洲| 欧美成人午夜免费资源| 亚洲,欧美,日韩| 亚洲精品亚洲一区二区| 97超视频在线观看视频| 亚洲欧美日韩卡通动漫| 亚洲精品成人av观看孕妇| 狂野欧美激情性bbbbbb| 亚洲国产精品国产精品| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 国产精品福利在线免费观看| 人人妻人人添人人爽欧美一区卜 | 亚洲色图av天堂| 日韩强制内射视频| 亚洲av免费高清在线观看| .国产精品久久| 久久精品久久精品一区二区三区| 男女国产视频网站| 春色校园在线视频观看| 日韩一区二区三区影片| 美女视频免费永久观看网站| 久久女婷五月综合色啪小说| 亚洲熟女精品中文字幕| 男女边吃奶边做爰视频| 春色校园在线视频观看| 亚洲欧美日韩另类电影网站 | 老女人水多毛片| 午夜免费男女啪啪视频观看| 亚洲欧美精品自产自拍| 精品一区二区免费观看| 一二三四中文在线观看免费高清| 韩国av在线不卡| 中文在线观看免费www的网站| 国产伦精品一区二区三区四那| 国产在线视频一区二区| 国产一区二区三区av在线| 日本vs欧美在线观看视频 | 少妇的逼好多水| 美女cb高潮喷水在线观看| 在线免费观看不下载黄p国产| 老师上课跳d突然被开到最大视频| 在线精品无人区一区二区三 | 一区二区三区乱码不卡18| 麻豆乱淫一区二区| 亚洲在久久综合| 免费高清在线观看视频在线观看| 国产亚洲最大av| 国产精品国产三级国产专区5o| 久久久久久久国产电影| 欧美丝袜亚洲另类| 亚洲av日韩在线播放| 久久久久久久亚洲中文字幕| 久久精品国产鲁丝片午夜精品| 91精品国产九色| 国产一级毛片在线| 夜夜骑夜夜射夜夜干| 亚洲在久久综合| 大片免费播放器 马上看| 日日撸夜夜添| 亚洲精品日韩在线中文字幕| 欧美日韩综合久久久久久| 黄色日韩在线| 午夜福利网站1000一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 免费观看av网站的网址| 亚洲欧美清纯卡通| 只有这里有精品99| 久久毛片免费看一区二区三区| 亚洲av综合色区一区| 黄色怎么调成土黄色| 免费大片18禁| av在线播放精品| 色哟哟·www| 国产精品国产三级国产专区5o| 久久人人爽人人爽人人片va| 99久久人妻综合| 国产av国产精品国产| 青春草亚洲视频在线观看| 亚洲成人中文字幕在线播放| 99国产精品免费福利视频| 中文在线观看免费www的网站| 国产一区有黄有色的免费视频| 赤兔流量卡办理| 亚洲无线观看免费| 国产精品99久久久久久久久| 成年女人在线观看亚洲视频| 亚洲图色成人| 在线观看国产h片| 3wmmmm亚洲av在线观看| 久久国产乱子免费精品| 最新中文字幕久久久久| 亚洲熟女精品中文字幕| 一级毛片久久久久久久久女| 精品亚洲成国产av| 国产永久视频网站| 日本爱情动作片www.在线观看| 丰满少妇做爰视频| 久久久久视频综合| kizo精华| 男人添女人高潮全过程视频| 天堂俺去俺来也www色官网| 久久久久国产精品人妻一区二区| 丝袜脚勾引网站| 特大巨黑吊av在线直播| 在线精品无人区一区二区三 | 人妻少妇偷人精品九色| 99九九线精品视频在线观看视频| 亚洲欧洲国产日韩| 欧美亚洲 丝袜 人妻 在线| 亚洲av中文av极速乱| 寂寞人妻少妇视频99o| 高清黄色对白视频在线免费看 | 久久久a久久爽久久v久久| 免费黄网站久久成人精品| 国产91av在线免费观看| 99热6这里只有精品| 亚洲精品第二区| a级毛色黄片| 欧美极品一区二区三区四区| 一个人免费看片子| 欧美日本视频| 色网站视频免费| 欧美另类一区| 晚上一个人看的免费电影| 国产探花极品一区二区| 精品久久久久久久久av| 亚洲精品成人av观看孕妇| 久久婷婷青草| 欧美日韩在线观看h| 啦啦啦中文免费视频观看日本| 国产精品一二三区在线看| 欧美老熟妇乱子伦牲交| 免费播放大片免费观看视频在线观看| 国产淫片久久久久久久久| 久热这里只有精品99| 久久久亚洲精品成人影院| 嫩草影院新地址| 国产爽快片一区二区三区| 成人免费观看视频高清| 亚洲av成人精品一二三区| av女优亚洲男人天堂| 国产欧美日韩精品一区二区| 亚洲av中文av极速乱| 亚洲婷婷狠狠爱综合网| 男女无遮挡免费网站观看| 91久久精品电影网| 亚洲成人手机| 熟妇人妻不卡中文字幕| 欧美成人午夜免费资源| 97在线视频观看| 日本猛色少妇xxxxx猛交久久| av卡一久久| 国产精品99久久99久久久不卡 | 青青草视频在线视频观看| 国产免费一区二区三区四区乱码| 超碰97精品在线观看| 国产色爽女视频免费观看| 国产成人午夜福利电影在线观看| 欧美xxⅹ黑人| 街头女战士在线观看网站| 高清av免费在线| 日本与韩国留学比较| av视频免费观看在线观看| 色视频www国产| 久久久久久久久久久丰满| 日日啪夜夜爽| 中文字幕人妻熟人妻熟丝袜美| 久久国产亚洲av麻豆专区| 五月天丁香电影| 乱系列少妇在线播放| 男人和女人高潮做爰伦理| 亚洲内射少妇av| 丰满人妻一区二区三区视频av| 免费观看无遮挡的男女| 亚洲国产精品国产精品| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单| 亚洲激情五月婷婷啪啪| 成人国产av品久久久| 777米奇影视久久| 免费大片18禁| 日韩一区二区三区影片| 日韩视频在线欧美| 免费av不卡在线播放| 精品久久久久久久久亚洲| 在线观看三级黄色| 王馨瑶露胸无遮挡在线观看| 亚洲图色成人| 最近最新中文字幕大全电影3| 在线看a的网站| 精品亚洲成a人片在线观看 | 国产男女超爽视频在线观看| 精品久久久久久久久亚洲| 卡戴珊不雅视频在线播放| 日韩一区二区三区影片| 欧美97在线视频| 国产成人精品福利久久| 91久久精品电影网| av免费在线看不卡| 免费观看性生交大片5| 午夜激情福利司机影院| 国产久久久一区二区三区| 熟女av电影| 国产精品偷伦视频观看了| 国产 一区精品| 国产伦精品一区二区三区视频9| 视频中文字幕在线观看| 舔av片在线| 精品久久久久久久末码| 久久久久网色| 久久久久久人妻| 人人妻人人看人人澡| 在线看a的网站| 一区二区av电影网| 国产男女内射视频| 高清av免费在线| 美女xxoo啪啪120秒动态图| 国产av国产精品国产| 久久国产乱子免费精品| 国产亚洲午夜精品一区二区久久| 日韩免费高清中文字幕av| 免费看av在线观看网站| 欧美日韩亚洲高清精品| 成人一区二区视频在线观看| 色综合色国产| 九草在线视频观看| 国产亚洲5aaaaa淫片| 久久久久久久久大av| 五月天丁香电影| 国产免费又黄又爽又色| 免费观看在线日韩| a级毛片免费高清观看在线播放| 中文字幕制服av| 在线天堂最新版资源| 我的女老师完整版在线观看| 国语对白做爰xxxⅹ性视频网站| 中文在线观看免费www的网站| 三级经典国产精品| 一本—道久久a久久精品蜜桃钙片| 亚洲四区av| 少妇的逼水好多| 成人黄色视频免费在线看| 久久午夜福利片| 亚洲人成网站在线观看播放| 精品一区二区免费观看| 欧美xxxx性猛交bbbb| 人妻制服诱惑在线中文字幕| 国产精品国产三级国产av玫瑰| 在线免费观看不下载黄p国产| 国产白丝娇喘喷水9色精品| 各种免费的搞黄视频| 麻豆乱淫一区二区| 国产白丝娇喘喷水9色精品| 国产成人精品婷婷| 国产精品伦人一区二区| 久久久久久久久久成人| 我要看日韩黄色一级片| 97在线人人人人妻| 这个男人来自地球电影免费观看 | 婷婷色综合www| 成人午夜精彩视频在线观看| 久久精品人妻少妇| 人妻少妇偷人精品九色| 亚洲av二区三区四区| 色哟哟·www| 中文字幕亚洲精品专区| 成人毛片60女人毛片免费| 国产综合精华液| 韩国av在线不卡| 老熟女久久久| 国产亚洲午夜精品一区二区久久| 在线观看国产h片| 偷拍熟女少妇极品色| 18禁裸乳无遮挡免费网站照片| 成年av动漫网址| 99国产精品免费福利视频| 国内少妇人妻偷人精品xxx网站| 97在线视频观看| 一边亲一边摸免费视频| 少妇人妻一区二区三区视频| 国产毛片在线视频| 国产成人freesex在线| 欧美最新免费一区二区三区| 久久久午夜欧美精品| 亚洲自偷自拍三级| 亚洲综合精品二区| 亚洲人成网站高清观看| 中文字幕av成人在线电影| 国产精品欧美亚洲77777| 国产高清有码在线观看视频| 男人添女人高潮全过程视频| 亚洲伊人久久精品综合| 免费黄色在线免费观看| 国产精品一区二区在线观看99| 亚洲av福利一区| 建设人人有责人人尽责人人享有的 | 色吧在线观看| 妹子高潮喷水视频| 人妻夜夜爽99麻豆av| 国产精品女同一区二区软件| 中文字幕人妻熟人妻熟丝袜美| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 久久久欧美国产精品| 色婷婷av一区二区三区视频| 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 亚洲色图综合在线观看| av福利片在线观看| 99九九线精品视频在线观看视频| 高清视频免费观看一区二区| 日韩中文字幕视频在线看片 | 色视频www国产| 亚洲精品乱码久久久久久按摩| av在线观看视频网站免费| 我要看黄色一级片免费的| 亚洲国产高清在线一区二区三| 午夜免费鲁丝| 亚洲三级黄色毛片| 多毛熟女@视频| 国产日韩欧美在线精品| 亚洲欧美成人综合另类久久久| 三级国产精品片| 久久精品人妻少妇| 免费观看的影片在线观看| 久久久久久久久久成人| 在线观看美女被高潮喷水网站| 夜夜爽夜夜爽视频| av又黄又爽大尺度在线免费看| 日本黄色片子视频| 日本午夜av视频| 国产精品一二三区在线看| 永久免费av网站大全| 亚洲精品久久午夜乱码| 成人免费观看视频高清| 久久 成人 亚洲| 欧美高清成人免费视频www| 日韩中字成人| 纯流量卡能插随身wifi吗| 午夜福利高清视频| 国产精品久久久久成人av| 国产精品人妻久久久久久| 少妇的逼好多水| 高清在线视频一区二区三区| 在线看a的网站| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 黄色怎么调成土黄色| 五月开心婷婷网| av一本久久久久| 青春草视频在线免费观看| 中文字幕av成人在线电影| 国产淫片久久久久久久久| 日本免费在线观看一区| 最近中文字幕高清免费大全6| 人妻制服诱惑在线中文字幕| 视频中文字幕在线观看| 赤兔流量卡办理| 综合色丁香网| 另类亚洲欧美激情| 国产一区亚洲一区在线观看| 久热这里只有精品99| 丰满人妻一区二区三区视频av| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 国产精品99久久99久久久不卡 | 18禁在线播放成人免费| 国产伦在线观看视频一区| 亚洲精品国产成人久久av| 最近最新中文字幕大全电影3| 大码成人一级视频| 国产高清三级在线| 丰满人妻一区二区三区视频av| 精品视频人人做人人爽| 久久人人爽人人爽人人片va| 亚洲人成网站在线播| 精品国产三级普通话版| 干丝袜人妻中文字幕| av网站免费在线观看视频| 久久人人爽人人片av| 国产免费一级a男人的天堂| 久久午夜福利片| 国产精品一区二区性色av| 美女cb高潮喷水在线观看| 大陆偷拍与自拍| 99久久综合免费| 欧美+日韩+精品| 91久久精品国产一区二区成人| 亚洲欧美日韩东京热| 国产精品精品国产色婷婷| 久久久精品免费免费高清| xxx大片免费视频| 久久精品国产a三级三级三级| 久久国产乱子免费精品| 亚洲av免费高清在线观看| 国产一区二区三区综合在线观看 | 欧美bdsm另类| 久久精品久久久久久久性| 在线精品无人区一区二区三 | 少妇人妻精品综合一区二区| 国产精品福利在线免费观看| 女性被躁到高潮视频| 黑丝袜美女国产一区| 18禁在线无遮挡免费观看视频| 伦理电影大哥的女人| 青青草视频在线视频观看| 欧美高清成人免费视频www| 精品久久国产蜜桃| 久久久精品免费免费高清| 一边亲一边摸免费视频| 18禁动态无遮挡网站| 一级毛片黄色毛片免费观看视频| 夫妻午夜视频| 嘟嘟电影网在线观看| 国产男女超爽视频在线观看| 九色成人免费人妻av| 婷婷色av中文字幕| 日本免费在线观看一区| 激情 狠狠 欧美| 在线免费观看不下载黄p国产| 久久av网站| 日韩不卡一区二区三区视频在线| 嫩草影院新地址| av在线老鸭窝| 国内少妇人妻偷人精品xxx网站| 亚洲精品久久午夜乱码| 性色av一级| 国产男女超爽视频在线观看| 亚洲国产欧美人成| 欧美成人精品欧美一级黄| 国产av精品麻豆| 欧美性感艳星| 国产日韩欧美亚洲二区| 91久久精品电影网| 国产亚洲91精品色在线| 日韩 亚洲 欧美在线| 久久av网站| 日日撸夜夜添| 国产亚洲一区二区精品| 欧美精品一区二区免费开放| 亚洲成人一二三区av| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 街头女战士在线观看网站| 久久精品久久久久久久性| 成人黄色视频免费在线看| 国产69精品久久久久777片| 亚洲av中文av极速乱| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 超碰av人人做人人爽久久| 国产高清三级在线| 国产精品偷伦视频观看了| 高清欧美精品videossex| 色综合色国产| 日韩欧美 国产精品| 97在线视频观看| 国内精品宾馆在线| 亚洲中文av在线| 一级a做视频免费观看| 一本—道久久a久久精品蜜桃钙片| 热re99久久精品国产66热6| 亚洲欧美日韩东京热| 啦啦啦在线观看免费高清www| 免费观看无遮挡的男女| 在线观看一区二区三区激情| av在线观看视频网站免费| 日本一二三区视频观看| 尤物成人国产欧美一区二区三区| av免费观看日本| 久久人人爽人人片av| 最近中文字幕2019免费版| 天美传媒精品一区二区| 国产在线一区二区三区精| 亚洲成人av在线免费| 国产精品三级大全| 国产高清三级在线| 伦理电影大哥的女人| 男女边吃奶边做爰视频| 亚洲精品中文字幕在线视频 | 久久久久视频综合| 亚洲图色成人| 国产亚洲5aaaaa淫片| 日本av手机在线免费观看| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 欧美极品一区二区三区四区| 免费观看在线日韩| 国产综合精华液| 日本欧美视频一区| 视频区图区小说| 国产精品伦人一区二区| 国产精品国产三级专区第一集| 国语对白做爰xxxⅹ性视频网站| 精品酒店卫生间| 国国产精品蜜臀av免费| 国产乱来视频区| 熟女电影av网| 国产黄片视频在线免费观看| 在线精品无人区一区二区三 | 午夜免费观看性视频| 永久免费av网站大全| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| 搡老乐熟女国产| 中文字幕制服av| 成年美女黄网站色视频大全免费 | 爱豆传媒免费全集在线观看| 久热久热在线精品观看| 肉色欧美久久久久久久蜜桃| 久久人人爽av亚洲精品天堂 | 中文字幕人妻熟人妻熟丝袜美| 日韩成人伦理影院| 亚洲美女视频黄频| 成人二区视频| 久久国产乱子免费精品| 国产老妇伦熟女老妇高清| 国产毛片在线视频| 久久精品国产自在天天线| 亚洲精品第二区| 看免费成人av毛片| 97超视频在线观看视频| 久久精品国产亚洲av天美| 国产亚洲欧美精品永久| 51国产日韩欧美| 精品国产一区二区三区久久久樱花 | 亚洲欧美日韩另类电影网站 | 在线观看av片永久免费下载| 久久韩国三级中文字幕| 午夜免费男女啪啪视频观看| 有码 亚洲区| 国产乱人偷精品视频| 久久99精品国语久久久| 一级黄片播放器| 国产有黄有色有爽视频| 成人影院久久| 一个人免费看片子| 日本wwww免费看| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 简卡轻食公司|