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

    基于Flow-3D的抗拖網(wǎng)海床基沖刷數(shù)值模擬

    2017-11-24 10:59:51于凱本林廣義孟慶健
    海洋學(xué)研究 2017年3期
    關(guān)鍵詞:拖網(wǎng)海床沙土

    于凱本,楊 濤,高 健,林廣義,孟慶健,宗 樂

    (1.國家深?;毓芾碇行?,山東 青島 266200;2.青島科技大學(xué) 機(jī)電工程學(xué)院,山東 青島 266100)

    基于Flow-3D的抗拖網(wǎng)海床基沖刷數(shù)值模擬

    于凱本1,楊 濤2,高 健2,林廣義2,孟慶健1,宗 樂1

    (1.國家深?;毓芾碇行?,山東 青島 266200;2.青島科技大學(xué) 機(jī)電工程學(xué)院,山東 青島 266100)

    基于Flow-3D軟件,建立抗拖網(wǎng)海床基沙土基床沖刷的三維數(shù)學(xué)模型,控制水流流速和泥沙粒徑,對(duì)不同條件下抗拖網(wǎng)海床基周圍沙土基床的沖刷進(jìn)行數(shù)值模擬,并分析了抗拖網(wǎng)海床基最大沖刷深度隨這兩種因素改變的趨勢。結(jié)果表明,抗拖網(wǎng)海床基沙土基床沖刷坑深度隨海流流速增大而增加,在設(shè)定其他條件不改變的情況下,沖刷坑的加深速度隨著流速增大會(huì)逐漸減緩,最終達(dá)到?jīng)_淤平衡。同時(shí),由于海底環(huán)境復(fù)雜多變,導(dǎo)致模型與實(shí)際情況有一定的差異,還需在后續(xù)的研究中進(jìn)一步進(jìn)行優(yōu)化。

    抗拖網(wǎng)海床基;泥沙沖刷;數(shù)值模擬

    0 引言

    海床基作為一種坐底式海洋監(jiān)測裝置,一般布放于海洋沿岸、近海潛岸以及河口等區(qū)域,對(duì)懸浮泥沙參數(shù)、流速剖面、水溫、潮汐及波浪等水流動(dòng)力學(xué)參數(shù)進(jìn)行長期、連續(xù)跟蹤監(jiān)測。由于布放區(qū)域是污染物、泥沙以及水流最大的區(qū)域,作為坐底長期環(huán)境監(jiān)測設(shè)備在大濃度泥沙水域極易被泥沙覆蓋,高的水流速度易使海床基平臺(tái)底部海底沉積物被掏空導(dǎo)致平臺(tái)傾斜倒塌。海床基的覆蓋和倒塌會(huì)降低海床基的回收效率,降低了海床基搭載設(shè)備的安全性和可靠性,增加海床基的應(yīng)用制造成本,導(dǎo)致不必要的經(jīng)濟(jì)損失[1]。因此,海床基坐底狀態(tài)下周圍沙土基床沖刷機(jī)理的研究對(duì)于提高海床基的回收安全性及對(duì)海床基外形優(yōu)化具有重要作用。

    本文運(yùn)用美國Flow science公司開發(fā)的計(jì)算流體動(dòng)力學(xué)(CFD)軟件Flow-3D中FAVOR技巧以及有關(guān)自由液面的VOF方法對(duì)海床基處于不同流速和不同泥沙粒徑的沖刷機(jī)理進(jìn)行了仿真模擬。以期為海床基的布放及整體結(jié)構(gòu)優(yōu)化提供一定的參考。

    1 控制方程

    針對(duì)本文擬計(jì)算算例,所建立的數(shù)值水池需模擬三維海床面上抗拖網(wǎng)海床基周圍泥沙在單向海流作用下的沖刷情況,即海流繞過抗拖網(wǎng)海床基與海床之間的相互作用。本文采用RNGk-ε模型進(jìn)行求解。RNGk-ε模型對(duì)湍動(dòng)粘度進(jìn)行修正,考慮了平均流動(dòng)中的旋轉(zhuǎn)及旋流流動(dòng)情況,可以更好地處理高應(yīng)變率及流線彎曲程度較大的流動(dòng)。

    RNGk-ε模型具體表示如下[2]:

    (1)

    (2)

    式中:xi為i方向坐標(biāo),t為時(shí)間(s),ρ和u分別代表流體的密度(kg/m3)和速度(m/s),gi為重力加速度(m/s2)。其中雷諾平均應(yīng)力定義為:

    (3)

    RNGk-ε湍流模型中的k方程和ε方程分別表示為:

    (4)

    (5)

    其中:

    2 沖刷沉積模型

    Flow-3D可以通過泥沙沖刷模塊對(duì)流體中泥沙運(yùn)動(dòng)狀況進(jìn)行分析。泥沙沖刷模塊是通過模擬在水流作用下的泥沙運(yùn)動(dòng)來預(yù)測泥沙侵蝕、平流擴(kuò)散和沉積。該模塊的主要應(yīng)用有:(1)計(jì)算懸移質(zhì)泥沙在平流水流中的擴(kuò)散輸移問題;(2)計(jì)算由于重力作用而產(chǎn)生的泥沙沉降問題;(3)計(jì)算由于剪切和水流擾動(dòng)而挾帶的泥沙問題;(4)計(jì)算由于泥沙顆粒隨著河床滾動(dòng)或者跳動(dòng)而引起的推移質(zhì)泥沙輸移問題。

    該模型的建立部分假定泥沙顆粒為球形并且他們的速度要小(從而使各個(gè)泥沙顆粒周圍的流體以粘性效應(yīng)占主導(dǎo)地位)。沉淀(漂移)系數(shù)Df通過式(6)自動(dòng)計(jì)算:

    (6)

    式中:SCRDIA是平均粒徑,RHOF是流體密度。沖刷的“起動(dòng)”組件是一個(gè)經(jīng)驗(yàn)?zāi)P?,該模型是以推移質(zhì)輸沙模型為基礎(chǔ)。在填充層的界面計(jì)算出的“起動(dòng)”速度是:

    (7)

    τc=SCRCRT×SCRDIA×g×(SCRRHO-RHOF)

    (8)

    式中:g為重力加速度;SCRRHO代表沉淀顆粒的密度,泥沙的典型密度值是2.65 g/cm3。

    當(dāng)沉積物含量大于等于沉積物臨界含量 (SCRFCR),沉積物即被認(rèn)為已經(jīng)堆積,堆積密度等于SCRFCR×SCRRHO。因此,在含沙量(CGS,單位為g/cm3)等于或大于SCRFCR×SCRRHO的地方有沉積物淤積,流體不會(huì)發(fā)生流動(dòng),施加在液體上的拖曳力被假定為無窮大(逆液拖曳數(shù)量,由Flow-3D計(jì)算的阻力系數(shù)為零)。在沉積物的體積分?jǐn)?shù)大于沉積物聚集含量(SCRFCO),但小于SCRFCR的區(qū)域,拖動(dòng)模型被激活;沉積物拖動(dòng)(SCRDRG)在這個(gè)模型中,相當(dāng)于TSDRG(默認(rèn)值是1.0)。當(dāng)沉積物的體積分?jǐn)?shù)小于SCRFCO,拖動(dòng)模型沒有被激活,但根據(jù)式(9)粘度會(huì)增強(qiáng):

    (9)

    式中:μ0是流體的粘度,μ是由于沉積物的存在引起的粘度增強(qiáng)。

    3 抗拖網(wǎng)海床基沙土基床沖刷模型

    3.1 沖刷模型

    本文對(duì)單向流作用下抗拖網(wǎng)海床基沙土基床沖刷進(jìn)行數(shù)值模擬分析,在Flow-3D軟件中建立與物理試驗(yàn)?zāi)P拖嗤娜S數(shù)值模型進(jìn)行模擬,通過對(duì)模擬結(jié)果分析,以驗(yàn)證沖刷模型條件參數(shù)及邊界條件的可靠性。

    抗拖網(wǎng)海床基模型圓臺(tái)形底面直徑φ=1.8 m,整體高度h=0.628 m。

    泥沙顆粒粒徑是仿真模擬的重要組成部分,粒徑的選取對(duì)模擬結(jié)果的準(zhǔn)確性有較大的影響。許炯心 等[3]利用1999年汛前深層鉆孔取樣資料和歷年淤積實(shí)測大斷面資料分析得到1960年前下游河道主槽淤積物中粒徑大于0.1 mm的泥沙占50.7%,這種泥沙在下游淤積比高達(dá)83.1%。所以在本次模擬中泥沙粒徑選取0.14 mm,泥沙休止角為31°,臨界希爾茲參數(shù)為0.05。數(shù)值水槽的長度為20 m,寬度為6 m,高度為4 m。模擬中以矩形模擬水域的下表面為基準(zhǔn)面,在此基準(zhǔn)面的基礎(chǔ)上建立高度為1 m的泥沙仿真區(qū)域。由于速度邊界提供的是沿整個(gè)水深均勻的流速,這樣就會(huì)在底部海床面泥沙范圍內(nèi)加入流速,不符合實(shí)際情況,且流速的存在可能會(huì)將入口附近的海床面泥沙帶走,為防止此情況發(fā)生,在左右邊界加入與海床面泥沙高度相同的實(shí)心擋板[4](圖1中模型的藍(lán)色和紅色部分)。

    圖1 抗拖網(wǎng)海床基沙土基床沖刷模型Fig.1 Scouring model of Trawl Resistant Seabed Basement on sand bed

    3.2 網(wǎng)格劃分

    利用CFD軟件中的網(wǎng)格劃分功能對(duì)抗拖網(wǎng)海床基沙土基床沖刷模型網(wǎng)格進(jìn)行了劃分,網(wǎng)格總數(shù)為100×60×50,形狀為均勻六面體網(wǎng)格。由于模型主要模擬抗拖網(wǎng)海床基周圍的流場變化和泥沙沖刷的地形變化,所以為了使模擬時(shí)計(jì)算能夠更精確,對(duì)抗拖網(wǎng)海床基周圍網(wǎng)格進(jìn)行了加密處理,如圖2和圖3所示。

    圖2 模型網(wǎng)格劃分平面圖Fig.2 Model meshing ichnography

    圖3 模型網(wǎng)格劃分立體圖Fig.3 Model meshing stereogram

    4 求解策略及模擬邊界條件設(shè)置

    本文采用廣義的GMRES算法,該方法的殘差極小,且收斂速度快,計(jì)算精度高,求解時(shí)不易產(chǎn)生發(fā)散,在求解N-S方程時(shí)與其他壓力速度的分離式解法相比有很高的效率。

    在沖刷模擬中需要考慮重力對(duì)結(jié)果的影響,因此本文設(shè)定重力矢量的大小為9.81 m/s。本文創(chuàng)建了兩相流(水和泥沙)模型,最底層為泥沙,高度為1 m;泥沙層上面的區(qū)域全部設(shè)置為海水。

    由于海水流速受季節(jié)、潮汐及深度等多方面的影響,參照文獻(xiàn)[5],對(duì)海水的模擬速度進(jìn)行選擇。

    表1 2000年南海東北部8—11月每月平均流速與方向[5]Tab.1 Monthly mean velocity and direction in the northeastern part of the South China Sea in Aug.-Nov.,2000[5]

    由表1可知,海水的速度隨著深度的增加而逐漸減小,并且因季節(jié)的不同而有所差異,為了模擬結(jié)果具有代表性,在海流速度的取值中取較大值。所以,進(jìn)流側(cè)設(shè)置恒定流速0.40 m/s,邊界條件的類型設(shè)置為Specified velocity。出口邊界:邊界條件的類型設(shè)置為Outflow。前后邊界和上下邊界:前后邊界采用默認(rèn)鏡像無通量邊界,其類型設(shè)置為Symmetry;下部邊界的類型設(shè)置為Wall。具體邊界類型如圖4所示。

    圖4 模型邊界條件設(shè)置Fig.4 Setting model boundary conditions

    5 沖刷的數(shù)值模擬與分析

    5.1 沖刷機(jī)理的數(shù)值模擬

    流速為0.40 m/s,泥沙粒徑為0.14 mm的條件下,抗拖網(wǎng)海床基沙土基床沖刷過程中不同時(shí)間下的沖刷坑深度值如表2所示。數(shù)值模擬過程中沖刷坑形成的速度比實(shí)際沖刷中沖刷坑形成的速度慢。當(dāng)沖刷開始時(shí),由于受到抗拖網(wǎng)海床基周圍水流和漩渦的作用,抗拖網(wǎng)海床基下面的泥沙被水流沖刷挾帶,在抗拖網(wǎng)海床基前端逐漸形成很小的沖刷坑,此過程中沖刷坑形成迅速,變化比較明顯。隨著沖刷的繼續(xù)進(jìn)行,抗拖網(wǎng)海床基沖刷坑的前后范圍逐漸增大,并向抗拖網(wǎng)海床基兩側(cè)延伸,沖刷坑的深度逐漸增加。由于抗拖網(wǎng)海床基屬于海底結(jié)構(gòu)物,整個(gè)平臺(tái)水平放置到海床面上,且抗拖網(wǎng)海床基對(duì)水流起阻礙作用,使得抗拖網(wǎng)海床基周圍的水流壓縮,導(dǎo)致局部的底流流速增大,從而抗拖網(wǎng)海床基底部周邊的泥沙被水流帶走。由于泥沙顆粒之間的相互摩擦力所起到阻礙作用,隨著水流的沖刷,從沖刷坑內(nèi)被水流沖走的泥沙愈來愈少,沖刷坑的加深速率逐漸變慢。最后,達(dá)到?jīng)_刷平衡,沖刷坑的深度不再增加,沖刷的結(jié)果如圖5所示。從圖5中可以清晰地看到從開始時(shí)刻到?jīng)_刷平衡時(shí)沖刷坑形成和床面高度變化的完整過程。

    表2 沖刷坑深度隨來流時(shí)間的變化Tab.2 Changes of scour depth with flow time

    圖5 抗拖網(wǎng)海床基沙土基床沖刷坑隨時(shí)間變化的情況Fig.5 Changes of scour with flow time around Trawl Resistant Seabed Basement

    水流遇到抗拖網(wǎng)海床基阻礙時(shí)形成了向下的沖擊水流,在水流的底部區(qū)域會(huì)形成與水流行進(jìn)方向相反的水平方向的反向漩渦(即馬蹄形漩渦),該漩渦對(duì)沙土基床沖刷起主要作用。同時(shí)因?yàn)闆_擊水流向下的沖擊作用,從而在抗拖網(wǎng)海床基迎流區(qū)域形成沖刷坑。而在抗拖網(wǎng)海床基背流區(qū)域由于向上水流和漩渦的共同作用,泥沙在抗拖網(wǎng)海床基背流區(qū)域內(nèi)被沖刷的效果同樣明顯。

    抗拖網(wǎng)海床基沙土基床沖刷不僅在前后端產(chǎn)生比較明顯的沖刷坑,抗拖網(wǎng)海床基左右兩側(cè)也會(huì)出現(xiàn)明顯的沖刷,從圖6可以看到,抗拖網(wǎng)海床基左右兩側(cè)形成了深度相當(dāng)?shù)臎_刷坑,沖刷坑范圍大概在1.8 m左右。在抗拖網(wǎng)海床基底部兩端由于馬蹄形漩渦和其他小漩渦的作用,抗拖網(wǎng)海床基周圍的沖刷坑都比較明顯。海床面也發(fā)生了相對(duì)比較大的變化,抗拖網(wǎng)海床基附近的床面沖刷變化比較大,在距離抗拖網(wǎng)海床基上游0.4 m和下游0.4 m處時(shí),床面基本沒變化。圖7為沖刷后的抗拖網(wǎng)海床基周圍地形的三維圖,其中左側(cè)為水流的來流方向。

    圖6 抗拖網(wǎng)海床基兩側(cè)的沖刷Fig.6 Scour on the sides of Trawl Resistant Seabed Basement

    圖7 沖刷后的抗拖網(wǎng)海床基周圍地形三維圖Fig.7 Three-dimensional map of terrain around Trawl Resistant Seabed Basement after scour

    5.2 不同流速下的泥沙沖刷分析

    為了更好地研究抗拖網(wǎng)海床基沙土基床的沖刷機(jī)理,在模擬中取泥沙粒徑d=0.5 mm,分別在行進(jìn)流速為0.20,0.25,0.30,0.40和0.50 m/s的條件下進(jìn)行模擬,具體的模擬結(jié)果如圖8所示。當(dāng)行進(jìn)流速為0.20 m/s時(shí),海床面上泥沙不起動(dòng),整個(gè)模擬過程中,海床面幾乎沒有變化,沖刷坑深度較小,為4.6 cm。當(dāng)流速為0.25 m/s時(shí),此時(shí)海床面上的泥沙也不起動(dòng),但是由于抗拖網(wǎng)海床基阻礙產(chǎn)生的向下阻流,使海床面產(chǎn)生泥沙沖刷,沖刷坑的深度相對(duì)較小,為7.52 cm。當(dāng)流速為0.50 m/s時(shí),海床面上的泥沙起動(dòng),最終達(dá)到?jīng)_刷平衡,其沖刷坑的深度為14.19 cm。

    圖8 沖刷坑深度隨流速的變化Fig.8 Changes of scour depth under different flow velocity

    由圖9可知,隨著流速的增大,抗拖網(wǎng)海床基沙土基床沖刷坑的深度也會(huì)相應(yīng)變大。這是抗拖網(wǎng)海床基沙土基床沖刷坑深度隨流速變化的總體趨勢。從圖中我們可以看到:

    (3)當(dāng)行進(jìn)流速V比海床面泥沙的起動(dòng)流速V0大并且逐漸增大時(shí),在海床面上的泥沙會(huì)發(fā)生起動(dòng),則沖刷坑內(nèi)會(huì)得到上游來沙的補(bǔ)充,沖刷坑的深度由于行進(jìn)流速的增大而增大的勢頭會(huì)減弱,此時(shí)為動(dòng)床沖刷。隨著行進(jìn)流速V的進(jìn)一步增大,沖刷坑深度也隨之在增大,但是沖刷速率明顯減小。圖中流速在0.40~0.50 m/s段可大體反映此種趨勢。

    圖9 不同流速下抗拖網(wǎng)海床基沖刷模擬圖Fig.9 Simulated diagram of scour around Trawl Resistant Seabed Basement under different flow velocity

    5.3 不同泥沙粒徑下的泥沙沖刷分析

    泥沙粒徑對(duì)抗拖網(wǎng)海床基沙土基床沖刷產(chǎn)生的影響,主要體現(xiàn)在泥沙顆粒的大小、形狀、阻力系數(shù)和組成結(jié)構(gòu)等四個(gè)方面。本文中,設(shè)定泥沙顆粒的中值粒徑D50分別為0.14,0.20,0.36,0.50和0.65 mm。行進(jìn)流速設(shè)為一個(gè)常數(shù),這里選取流速為0.50 m/s。沖刷結(jié)果具體見圖10和圖11。

    圖10 沖刷坑深度隨泥沙粒徑變化Fig.10 Changes of scour depth under different sand size

    但是在實(shí)際情況下,由于泥沙粒徑不同泥沙沖刷深度不同的現(xiàn)象并不顯著。在海底沉積物中,泥沙粒徑小則0.1 mm,大則1 000 mm,甚至更大,其比率相差有萬倍,有的甚至更多。而抗拖網(wǎng)海床基布放到海底后受海流作用而形成的泥沙沖刷坑的大小不過1.0~15 m,其比率相差僅15倍,前后兩者在數(shù)量級(jí)上差別很大。在數(shù)值模擬中,一般都是采用比尺模型進(jìn)行求解計(jì)算,所以沖刷坑大小范圍和泥沙顆粒大小范圍的變化數(shù)量級(jí)比較接近。由此可知,數(shù)值模擬的結(jié)果用于實(shí)際應(yīng)用需要慎重考慮。

    6 結(jié)論

    本文應(yīng)用Flow-3D軟件對(duì)抗拖網(wǎng)海床基沙土基床沖刷進(jìn)行數(shù)值模擬,并改變水流流速和泥沙粒徑,模擬不同條件下對(duì)沙土基床沖刷的影響。

    圖11 不同泥沙粒徑下抗拖網(wǎng)海床基沖刷模擬圖Fig.11 Simulated diagram of scour around Trawl Resistant Seabed Basement under different sand size

    模擬結(jié)果能夠直觀展現(xiàn)抗拖網(wǎng)海床基沙土基床沖刷坑的形成過程及沖刷床面的變換情況??雇暇W(wǎng)海床基沙土基床沖刷坑深度隨海流流速增大而增加,在設(shè)定其他條件不改變的情況下,沖刷坑的加深速度隨著流速增大會(huì)逐漸減緩,最終達(dá)到?jīng)_淤平衡。在其他條件不變的情況下,隨泥沙粒徑的增大,沖刷坑深度逐漸減小,而且減小的速度逐漸減緩并趨于不沖刷。

    本文的模擬結(jié)果還不能完全重現(xiàn)抗拖網(wǎng)海床基布放在海底后沖刷坑的實(shí)際情況,這主要是由于海底環(huán)境復(fù)雜多變,導(dǎo)致模型與實(shí)際情況有一定的差異,因此還需要使用更加準(zhǔn)確的海床面模型及初始條件對(duì)其進(jìn)行模擬。本研究希望能為抗拖網(wǎng)海床基結(jié)構(gòu)優(yōu)化設(shè)計(jì)和實(shí)際應(yīng)用提供理論參考。

    [1] LIU Jie, WU Cheng-xuan, Lü bin, et al. CFD based flow field analysis of a seabed monitoring platform[J]. Shandong Science,2012,25(1):65-68.

    劉杰,吳承璇,呂斌,等.基于CFD的海床基觀測平臺(tái)外部流場分析[J].山東科學(xué),2012,25(1):65-68.

    [2] WANG Fu-jun. Principle and application of CFD software for computational fluid dynamics analysis[M].Beijing:Tsinghua University Press,2004.

    王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004.

    [3] XU Jiong-xin, HU Chun-hong, CHEN Jian-guo. The influence of sediment of different size groups on the deposition of the lower reaches of the Yellow River and its significance in the management of the Yellow River[J]. Science in China: Series E,2009,39(2):310-317.

    許炯心,胡春宏,陳建國.不同粒徑組泥沙對(duì)黃河下游沉積的影響及其在黃河治理中的意義[J].中國科學(xué):E輯,2009,39(2):310-317.

    [4] ZHAO Yan-fei.Three-dimensional numerical simulation of wave force on the offshore wind turbine structure and scour around foundation[D]. Tianjing:Tianjin University,2010.

    趙雁飛.海上風(fēng)電支撐結(jié)構(gòu)波浪力及基礎(chǔ)沖刷的三維數(shù)值模擬研究[D].天津:天津大學(xué),2010.

    [5] YUAN Yao-chu, ZHAO Jin-ping, WANG Hui-qun, et al. Observations of 450 m currents in the shallow and deep layers of the northeastern South China Sea and their spectral analysis[J]. Science in China: Series D,2002,32(2):163-176.

    袁耀初,趙進(jìn)平,王惠群,等.南海東北部450 m以淺水層與深水層海流觀測結(jié)果及其譜分析[J].中國科學(xué):D輯,2002,32(2):163-176.

    NumericalsimulationanalysisofscouraroundTrawlResistantSeabedBasementbasedonFlow-3D

    YU Kai-ben1, YANG Tao2, GAO Jian2, LIN Guang-yi2, MENG Qing-jian1, ZONG Le1

    (1.NationalDeepSeaCenter,Qingdao266200,China; 2.CollegeofMechanicalandElectricalEngineering,QingdaoScience&TechnologyUniversity,Qingdao266100,China)

    Three-dimensional mathematical model of scour around Trawl Resistant Seabed Basement on sand bed was established using Flow-3D. By changing the flow velocity and sediment particle size, the scour on sand bed under different conditions was simulated, and the changes in the trend of maximum scour depth around Trawl Resistant Seabed Basement with the two factors was analyzed. The results show that the scour depth around Trawl Resistant Seabed Basement is increased with current velocity increases, under the other setting conditions unchanged, the scour pit deepened velocity with increasing flow velocity is gradually slowed down, and ultimately achieved sedimentation balance. At the same time, due to the complex and changeable seabed environment, the model and the actual situation are different, and therefore further research should be continued.

    Trawl Resistant Seabed Basement;sediment scour;numerical simulation

    于凱本,楊濤,高健,等.基于Flow-3D的抗拖網(wǎng)海床基沖刷數(shù)值模擬[J].海洋學(xué)研究,2017,35(3):91-98,

    10.3969/j.issn.1001-909X.2017.03.011.

    YU Kai-ben, YANG Tao, GAO Jian, et al. Numerical simulation analysis of scour around Trawl Resistant Seabed Basement based on Flow-3D[J].Journal of Marine Sciences,2017,35(3):91-98, doi:10.3969/j.issn.1001-909X.2017.03.011.

    2016-08-15

    2016-11-16

    “泰山學(xué)者工程”專項(xiàng)經(jīng)費(fèi)項(xiàng)目資助(tspd 20161007)

    于凱本(1977-),男,山東濟(jì)南市人,高級(jí)工程師,主要從事海洋檢測技術(shù)方面的研究。E-mail: kedayangtao@163.com

    P715.5;P737.2

    A

    1001-909X(2017)03-0091-08

    10.3969/j.issn.1001-909X.2017.03.011

    猜你喜歡
    拖網(wǎng)海床沙土
    拖網(wǎng)的人(外一首)
    拖網(wǎng)的人(外一首)
    開漁
    差別最大的字母
    山東90后小伙賣黃河沙土
    波流耦合作用下雙層砂質(zhì)海床累積液化特征數(shù)值分析?
    河北省框架拖網(wǎng)和桁桿拖網(wǎng)捕撈現(xiàn)狀調(diào)查
    沙土褲里的生命密碼
    波致砂土海床剪切與液化破壞特征對(duì)比研究?
    圍封治理對(duì)風(fēng)沙土養(yǎng)分含量的影響
    国产淫语在线视频| 两个人的视频大全免费| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 大片电影免费在线观看免费| 高清av免费在线| 亚洲一级一片aⅴ在线观看| 亚洲天堂av无毛| 男女边摸边吃奶| 免费高清在线观看视频在线观看| 免费观看a级毛片全部| 国产成人91sexporn| 一级片'在线观看视频| 国产精品精品国产色婷婷| 日日摸夜夜添夜夜爱| 毛片女人毛片| 亚洲国产色片| 日韩一区二区三区影片| 蜜桃久久精品国产亚洲av| 在线观看一区二区三区激情| 嫩草影院新地址| 噜噜噜噜噜久久久久久91| 亚洲国产色片| 精品久久久久久久人妻蜜臀av| 日本免费在线观看一区| 黑人高潮一二区| 精品久久久久久久末码| 国产亚洲av片在线观看秒播厂| 三级经典国产精品| 国内揄拍国产精品人妻在线| 久久精品久久精品一区二区三区| 一级黄片播放器| 狂野欧美白嫩少妇大欣赏| 老女人水多毛片| 观看免费一级毛片| 色5月婷婷丁香| 蜜臀久久99精品久久宅男| av一本久久久久| 乱系列少妇在线播放| 男人添女人高潮全过程视频| 亚洲精品亚洲一区二区| 噜噜噜噜噜久久久久久91| 国产淫片久久久久久久久| 观看免费一级毛片| 午夜福利网站1000一区二区三区| 国产av不卡久久| 最后的刺客免费高清国语| 尤物成人国产欧美一区二区三区| 国产免费福利视频在线观看| 久久人人爽人人爽人人片va| 国产精品久久久久久精品电影小说 | 国产老妇伦熟女老妇高清| 又爽又黄a免费视频| 三级经典国产精品| 亚洲精品色激情综合| 美女内射精品一级片tv| 又黄又爽又刺激的免费视频.| 亚洲成人av在线免费| 免费播放大片免费观看视频在线观看| videossex国产| 国产精品久久久久久精品古装| 最近2019中文字幕mv第一页| 2021天堂中文幕一二区在线观| av在线观看视频网站免费| 亚洲av男天堂| 麻豆国产97在线/欧美| 男女国产视频网站| 久久久精品欧美日韩精品| 国产亚洲最大av| 我要看日韩黄色一级片| 天美传媒精品一区二区| 亚洲av免费高清在线观看| 18禁在线无遮挡免费观看视频| 国产成人a∨麻豆精品| 国产黄片视频在线免费观看| 热99国产精品久久久久久7| 久久精品久久久久久久性| av在线老鸭窝| 亚洲国产日韩一区二区| 你懂的网址亚洲精品在线观看| 只有这里有精品99| 国产成人91sexporn| 精品一区在线观看国产| 好男人视频免费观看在线| 久久久久性生活片| 国产一区有黄有色的免费视频| 精品少妇黑人巨大在线播放| 国产精品国产三级专区第一集| 色哟哟·www| 熟女人妻精品中文字幕| 国产成人freesex在线| 欧美日韩精品成人综合77777| 日日摸夜夜添夜夜爱| 99热这里只有精品一区| .国产精品久久| 少妇人妻久久综合中文| 五月天丁香电影| 精品久久久噜噜| 亚洲欧美精品自产自拍| 日韩强制内射视频| 在线亚洲精品国产二区图片欧美 | 色视频在线一区二区三区| 中文字幕制服av| 九草在线视频观看| 免费不卡的大黄色大毛片视频在线观看| 免费观看的影片在线观看| 国产色婷婷99| 亚洲精品,欧美精品| 嫩草影院精品99| 青春草国产在线视频| 欧美人与善性xxx| 久久精品国产自在天天线| 美女内射精品一级片tv| 国产日韩欧美在线精品| 久久99精品国语久久久| 国产精品久久久久久av不卡| 天堂网av新在线| 国内精品美女久久久久久| 国产一区二区亚洲精品在线观看| 久久久久久久大尺度免费视频| 国产毛片a区久久久久| 亚洲伊人久久精品综合| 在线 av 中文字幕| 精品久久久久久久人妻蜜臀av| 国产精品成人在线| 欧美zozozo另类| 国产乱人偷精品视频| 国产成人精品久久久久久| 欧美少妇被猛烈插入视频| a级毛色黄片| 国产一区二区三区综合在线观看 | 18禁在线无遮挡免费观看视频| 男女那种视频在线观看| 成人毛片60女人毛片免费| 91久久精品国产一区二区三区| 国产69精品久久久久777片| 国产成人精品福利久久| 又粗又硬又长又爽又黄的视频| 午夜爱爱视频在线播放| 久久精品国产亚洲av涩爱| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 色视频在线一区二区三区| 成人亚洲精品av一区二区| h日本视频在线播放| 日韩三级伦理在线观看| 免费看日本二区| 99视频精品全部免费 在线| 又粗又硬又长又爽又黄的视频| 午夜激情福利司机影院| 白带黄色成豆腐渣| 国产国拍精品亚洲av在线观看| 国产黄色视频一区二区在线观看| 99久久九九国产精品国产免费| 亚洲一级一片aⅴ在线观看| 麻豆成人午夜福利视频| 日本熟妇午夜| 在线观看三级黄色| 干丝袜人妻中文字幕| 久久久久久久久久久免费av| 亚洲精品国产成人久久av| 嘟嘟电影网在线观看| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 欧美 日韩 精品 国产| 国产老妇女一区| 中文精品一卡2卡3卡4更新| 久久人人爽人人片av| 十八禁网站网址无遮挡 | 国产精品成人在线| 色播亚洲综合网| 极品少妇高潮喷水抽搐| freevideosex欧美| 国产免费视频播放在线视频| 日韩一区二区视频免费看| 亚洲熟女精品中文字幕| 国产老妇伦熟女老妇高清| 高清视频免费观看一区二区| 午夜激情久久久久久久| 狠狠精品人妻久久久久久综合| 欧美日韩一区二区视频在线观看视频在线 | 亚洲最大成人手机在线| 久久精品国产亚洲av涩爱| 色吧在线观看| 欧美3d第一页| 最近中文字幕2019免费版| 亚洲内射少妇av| 国产精品伦人一区二区| 亚洲第一区二区三区不卡| 欧美性感艳星| 免费人成在线观看视频色| 少妇人妻一区二区三区视频| 男插女下体视频免费在线播放| 免费看不卡的av| 精品国产三级普通话版| av女优亚洲男人天堂| 十八禁网站网址无遮挡 | 草草在线视频免费看| 亚洲av免费高清在线观看| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 国产精品一区www在线观看| 亚洲性久久影院| 国产成人免费观看mmmm| 国产亚洲精品久久久com| 日韩视频在线欧美| 91久久精品电影网| 国产一区亚洲一区在线观看| 欧美zozozo另类| 国产一级毛片在线| 人人妻人人爽人人添夜夜欢视频 | av免费在线看不卡| 成人漫画全彩无遮挡| 水蜜桃什么品种好| 国产欧美日韩精品一区二区| 永久网站在线| 中国三级夫妇交换| 国产亚洲精品久久久com| 99热这里只有是精品50| 亚洲电影在线观看av| 爱豆传媒免费全集在线观看| 男插女下体视频免费在线播放| 国产亚洲av片在线观看秒播厂| 日本免费在线观看一区| 国产精品一及| 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久精品电影小说 | 日日撸夜夜添| 国产精品久久久久久久久免| 日韩一本色道免费dvd| 一个人观看的视频www高清免费观看| 男女无遮挡免费网站观看| 3wmmmm亚洲av在线观看| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 偷拍熟女少妇极品色| 真实男女啪啪啪动态图| 男女下面进入的视频免费午夜| 国产精品99久久99久久久不卡 | 日韩 亚洲 欧美在线| 一个人看的www免费观看视频| 噜噜噜噜噜久久久久久91| 国产 一区 欧美 日韩| 国产精品成人在线| 老司机影院成人| 亚洲国产精品999| 久久久成人免费电影| 97精品久久久久久久久久精品| 网址你懂的国产日韩在线| 国产精品爽爽va在线观看网站| videos熟女内射| 亚洲精品,欧美精品| 国产精品嫩草影院av在线观看| 日韩在线高清观看一区二区三区| 午夜激情久久久久久久| 乱码一卡2卡4卡精品| 日产精品乱码卡一卡2卡三| 久久国内精品自在自线图片| 国产探花极品一区二区| 国产欧美另类精品又又久久亚洲欧美| 亚洲av不卡在线观看| 国产高潮美女av| 国产成人免费无遮挡视频| 午夜视频国产福利| 亚洲av福利一区| 18禁动态无遮挡网站| 22中文网久久字幕| 婷婷色综合www| 99re6热这里在线精品视频| 日本猛色少妇xxxxx猛交久久| av黄色大香蕉| 精品国产乱码久久久久久小说| 久久精品夜色国产| 交换朋友夫妻互换小说| 禁无遮挡网站| 99视频精品全部免费 在线| 久久精品久久精品一区二区三区| 伦精品一区二区三区| 久久久久久国产a免费观看| av在线app专区| av国产免费在线观看| 少妇人妻精品综合一区二区| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 插逼视频在线观看| 亚洲国产最新在线播放| 国产黄片美女视频| 丝袜美腿在线中文| 日韩,欧美,国产一区二区三区| 国产亚洲午夜精品一区二区久久 | 精品久久国产蜜桃| 亚洲人成网站高清观看| 亚洲综合色惰| 亚洲精品中文字幕在线视频 | 国产成人午夜福利电影在线观看| 欧美精品国产亚洲| 在线亚洲精品国产二区图片欧美 | av线在线观看网站| 久久国产乱子免费精品| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 亚洲综合精品二区| 美女脱内裤让男人舔精品视频| 精品少妇久久久久久888优播| 美女被艹到高潮喷水动态| 一本色道久久久久久精品综合| 欧美成人午夜免费资源| 免费播放大片免费观看视频在线观看| 亚洲精品国产av蜜桃| 欧美日韩在线观看h| 99热全是精品| 午夜日本视频在线| 在线免费十八禁| 99精国产麻豆久久婷婷| 六月丁香七月| 岛国毛片在线播放| 69av精品久久久久久| 国产伦在线观看视频一区| 97在线视频观看| 国产成人精品一,二区| 美女国产视频在线观看| 国产91av在线免费观看| 久久精品国产自在天天线| 女人被狂操c到高潮| 亚洲av中文av极速乱| 亚洲欧美中文字幕日韩二区| 两个人的视频大全免费| 简卡轻食公司| 日韩人妻高清精品专区| 水蜜桃什么品种好| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 校园人妻丝袜中文字幕| 欧美激情国产日韩精品一区| 免费大片18禁| 高清视频免费观看一区二区| 日本av手机在线免费观看| 亚洲精华国产精华液的使用体验| 在线观看国产h片| 亚洲丝袜综合中文字幕| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 免费看光身美女| 久久综合国产亚洲精品| 日韩国内少妇激情av| 一级二级三级毛片免费看| 国产日韩欧美在线精品| 国产精品人妻久久久影院| 老女人水多毛片| 久久久久久久久大av| 欧美极品一区二区三区四区| 亚洲美女视频黄频| 涩涩av久久男人的天堂| 亚州av有码| 国产高清三级在线| 国产成年人精品一区二区| 最近最新中文字幕大全电影3| 精品一区二区免费观看| 新久久久久国产一级毛片| 久久久久精品久久久久真实原创| 亚洲精品一区蜜桃| 亚洲av.av天堂| 精品人妻一区二区三区麻豆| 亚洲av免费在线观看| tube8黄色片| 色吧在线观看| 精品一区二区三卡| 亚洲内射少妇av| 可以在线观看毛片的网站| 2018国产大陆天天弄谢| 国产伦精品一区二区三区四那| 精品久久国产蜜桃| 91精品一卡2卡3卡4卡| 舔av片在线| 日韩成人伦理影院| 舔av片在线| 亚洲无线观看免费| 免费观看性生交大片5| 国产一区二区三区综合在线观看 | 高清毛片免费看| 国产精品秋霞免费鲁丝片| 99久久人妻综合| 国产精品av视频在线免费观看| 国产欧美日韩一区二区三区在线 | 国产乱来视频区| 人妻 亚洲 视频| av在线蜜桃| 亚洲色图综合在线观看| tube8黄色片| 2021天堂中文幕一二区在线观| 亚洲aⅴ乱码一区二区在线播放| 又黄又爽又刺激的免费视频.| 99热国产这里只有精品6| 蜜桃久久精品国产亚洲av| 亚洲怡红院男人天堂| 亚洲精品一区蜜桃| 精品午夜福利在线看| 在线观看一区二区三区激情| 国产成人福利小说| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 九草在线视频观看| 嫩草影院入口| 寂寞人妻少妇视频99o| 国产精品熟女久久久久浪| 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 久久99蜜桃精品久久| 欧美丝袜亚洲另类| 日本色播在线视频| 精品亚洲乱码少妇综合久久| 亚洲国产欧美在线一区| 小蜜桃在线观看免费完整版高清| 中国国产av一级| 日韩国内少妇激情av| 91精品伊人久久大香线蕉| 三级男女做爰猛烈吃奶摸视频| av国产久精品久网站免费入址| 97超碰精品成人国产| 亚洲内射少妇av| 亚洲精品日韩av片在线观看| 午夜免费鲁丝| 国产乱人偷精品视频| 久久国产乱子免费精品| 禁无遮挡网站| 成人高潮视频无遮挡免费网站| 亚洲精品国产成人久久av| 亚洲av成人精品一区久久| 久久久久久伊人网av| 蜜桃亚洲精品一区二区三区| 午夜福利在线在线| 在线免费观看不下载黄p国产| 国内少妇人妻偷人精品xxx网站| 久久国内精品自在自线图片| av天堂中文字幕网| 国产黄片视频在线免费观看| 91精品一卡2卡3卡4卡| 欧美国产精品一级二级三级 | 国产精品久久久久久久电影| 国产视频首页在线观看| 国产成人午夜福利电影在线观看| 女人十人毛片免费观看3o分钟| av天堂中文字幕网| 爱豆传媒免费全集在线观看| 人妻 亚洲 视频| 免费大片18禁| 久久久久久久久久人人人人人人| 亚洲精品日本国产第一区| 日日啪夜夜撸| 熟妇人妻不卡中文字幕| 欧美精品国产亚洲| 人人妻人人看人人澡| 一级a做视频免费观看| 乱系列少妇在线播放| 99热这里只有精品一区| 亚洲欧美精品专区久久| 久久精品夜色国产| 成人高潮视频无遮挡免费网站| 久久久国产一区二区| 天堂俺去俺来也www色官网| 永久免费av网站大全| 欧美激情在线99| 内射极品少妇av片p| 久久影院123| 欧美潮喷喷水| 色吧在线观看| 高清毛片免费看| 久久99蜜桃精品久久| 卡戴珊不雅视频在线播放| 亚洲欧美日韩东京热| 大片电影免费在线观看免费| 直男gayav资源| 亚洲av国产av综合av卡| 好男人视频免费观看在线| 国产精品99久久久久久久久| 日韩av免费高清视频| 男人添女人高潮全过程视频| 精品久久久久久电影网| 91精品国产九色| 国产精品成人在线| 国产亚洲一区二区精品| 午夜福利在线在线| 国产精品国产三级国产专区5o| 国产黄色免费在线视频| 亚洲欧美日韩另类电影网站 | 久久久精品免费免费高清| 国国产精品蜜臀av免费| 成年女人在线观看亚洲视频 | 乱码一卡2卡4卡精品| 九九在线视频观看精品| 成人黄色视频免费在线看| 人妻少妇偷人精品九色| 国产老妇女一区| 亚洲精品久久午夜乱码| 99久久精品一区二区三区| 国产成人精品婷婷| 免费黄频网站在线观看国产| 99久久精品热视频| 国产久久久一区二区三区| 又大又黄又爽视频免费| 精品99又大又爽又粗少妇毛片| 男女下面进入的视频免费午夜| 亚洲久久久久久中文字幕| 国产精品秋霞免费鲁丝片| 亚洲精品一区蜜桃| 视频区图区小说| 99久久人妻综合| 国产成人免费观看mmmm| 真实男女啪啪啪动态图| av又黄又爽大尺度在线免费看| 男女那种视频在线观看| 精品午夜福利在线看| 成人一区二区视频在线观看| 一级a做视频免费观看| 赤兔流量卡办理| 青春草亚洲视频在线观看| 欧美激情国产日韩精品一区| 久久97久久精品| 亚洲经典国产精华液单| 成人无遮挡网站| 一级毛片黄色毛片免费观看视频| 国内精品宾馆在线| 国产高潮美女av| 久久精品夜色国产| 黄片无遮挡物在线观看| 欧美日韩亚洲高清精品| 精品国产三级普通话版| 97在线人人人人妻| 各种免费的搞黄视频| 日韩视频在线欧美| 国产一区二区亚洲精品在线观看| 人人妻人人看人人澡| 超碰av人人做人人爽久久| 夜夜看夜夜爽夜夜摸| 啦啦啦啦在线视频资源| 18禁裸乳无遮挡动漫免费视频 | 国产高清国产精品国产三级 | 丰满人妻一区二区三区视频av| 久久久欧美国产精品| 成人无遮挡网站| 久久99精品国语久久久| tube8黄色片| 视频中文字幕在线观看| 亚洲欧美精品自产自拍| 久久久久久久大尺度免费视频| 国产黄a三级三级三级人| 女人久久www免费人成看片| 狂野欧美激情性bbbbbb| 亚洲伊人久久精品综合| 亚洲无线观看免费| 99久国产av精品国产电影| 少妇人妻一区二区三区视频| 青春草视频在线免费观看| 欧美激情在线99| 久久热精品热| 尤物成人国产欧美一区二区三区| 一级二级三级毛片免费看| 亚洲第一区二区三区不卡| 国产v大片淫在线免费观看| 国产黄a三级三级三级人| 国产黄片美女视频| 日韩一区二区三区影片| 亚洲自偷自拍三级| 国产毛片a区久久久久| 欧美xxxx性猛交bbbb| 美女国产视频在线观看| 高清视频免费观看一区二区| 在线 av 中文字幕| 亚洲欧美精品专区久久| 亚洲美女搞黄在线观看| 亚洲精品色激情综合| 特级一级黄色大片| 国产精品久久久久久精品电影小说 | 熟女电影av网| 嘟嘟电影网在线观看| 国产亚洲午夜精品一区二区久久 | 久久精品夜色国产| 免费电影在线观看免费观看| 国产探花极品一区二区| 熟女人妻精品中文字幕| 亚洲av.av天堂| 十八禁网站网址无遮挡 | 久久国内精品自在自线图片| 精品久久久久久久久亚洲| 久热这里只有精品99| 丰满少妇做爰视频| 日本色播在线视频| 午夜视频国产福利| 午夜福利视频1000在线观看| 欧美bdsm另类| 国产精品人妻久久久久久| 免费观看的影片在线观看| 欧美bdsm另类| 亚洲精品一二三| 99久国产av精品国产电影| 成人国产麻豆网| 国产伦理片在线播放av一区| 在线 av 中文字幕| 亚洲伊人久久精品综合| 十八禁网站网址无遮挡 | 国产老妇伦熟女老妇高清| 草草在线视频免费看| 国产伦在线观看视频一区| 日韩欧美精品免费久久| 神马国产精品三级电影在线观看| 亚洲人成网站在线播| 国产男女内射视频| 日日摸夜夜添夜夜添av毛片| 成人国产av品久久久| 高清视频免费观看一区二区| 在线观看免费高清a一片|