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

    波浪布拉格共振研究進(jìn)展

    2016-07-20 07:59:25鄭金海余豪豐陶愛(ài)峰海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室河海大學(xué)江蘇南京210098河海大學(xué)港口海岸與近海工程學(xué)院江蘇南京210098
    水利水電科技進(jìn)展 2016年3期
    關(guān)鍵詞:述評(píng)波浪

    鄭金海,余豪豐,陶愛(ài)峰,范 駿,王 懿(1.海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(河海大學(xué)),江蘇南京 210098;2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇南京 210098)

    ?

    波浪布拉格共振研究進(jìn)展

    鄭金海1,2,余豪豐1,2,陶愛(ài)峰1,2,范 駿1,2,王 懿1,2
    (1.海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(河海大學(xué)),江蘇南京 210098;2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇南京 210098)

    摘要:為給河口海岸防護(hù)工程設(shè)計(jì)和地形演變預(yù)測(cè)提供科學(xué)依據(jù),研究了波浪從外海向近岸傳播過(guò)程中在一定地形條件下產(chǎn)生的布拉格共振現(xiàn)象。介紹了波浪布拉格共振的定義、特征和研究意義,評(píng)述了現(xiàn)場(chǎng)觀測(cè)、理論研究、水槽試驗(yàn)、數(shù)值模擬等4種研究方法的優(yōu)缺點(diǎn),總結(jié)了波浪布拉格共振類型、低階共振反射系數(shù)公式和高階共振數(shù)值模擬的研究成果,提出了波浪布拉格共振主頻下移、不規(guī)則波布拉格共振特性、波流相互作用下布拉格共振機(jī)理等3個(gè)亟待深入研究的課題。

    關(guān)鍵詞:波浪;布拉格共振;非線性共振;主頻下移;波流相互作用;述評(píng)

    外海波浪傳入近岸淺水區(qū)域時(shí),一部分波浪能量會(huì)被海底地形反射,另一部分波浪能量則從地形上方繼續(xù)向岸傳播。1973年,Long[1]通過(guò)聯(lián)合北海波浪項(xiàng)目的實(shí)測(cè)數(shù)據(jù),分析了表面重力波通過(guò)不規(guī)則海底地形的情形,研究了當(dāng)?shù)貜?fù)雜地形對(duì)表面重力波的反射程度。1982年,Davies[2]發(fā)現(xiàn)許多海岸存在一系列的沿岸平行沙波,當(dāng)特定波長(zhǎng)的表面重力波通過(guò)這些周期性連續(xù)正弦型沙波時(shí),波浪與地形將產(chǎn)生共振現(xiàn)象。雖然沙波波高并不大,但由于兩者的共振作用會(huì)使得表面重力波大幅反射。Davies引入固體物理學(xué)中的布拉格反射原理對(duì)其進(jìn)行了解釋,將該現(xiàn)象稱為流體力學(xué)中的布拉格共振。由于周期性起伏地形引起的布拉格共振現(xiàn)象對(duì)于海岸防護(hù)工程具有重要價(jià)值,越來(lái)越多的科研人員開(kāi)始對(duì)這種共振現(xiàn)象產(chǎn)生興趣,并對(duì)其工程應(yīng)用展開(kāi)研究。

    1 波浪布拉格共振的定義和特征

    波浪布拉格共振是指在有限水深下,當(dāng)表面重力波通過(guò)周期性連續(xù)的正弦沙波地形時(shí),若沙波間距為表面重力波波長(zhǎng)一半的整數(shù)倍時(shí),波浪與地形將發(fā)生共振,此時(shí)表面重力波的反射最強(qiáng)。布拉格共振是1915年布拉格父子用X射線研究?jī)善叫芯w結(jié)構(gòu)時(shí)所發(fā)現(xiàn)的現(xiàn)象,當(dāng)晶體間距為X射線波長(zhǎng)一半的整數(shù)倍時(shí),X射線的反射最強(qiáng)。波浪和X射線的布拉格共振的本質(zhì)都是波列的相長(zhǎng)干涉。

    Davies等[3]基于正則攝動(dòng)法對(duì)周期性固定正弦地形上波浪布拉格共振的二維特征進(jìn)行了理論描述,指出在線性規(guī)則波和有限個(gè)沙波地形條件下,波浪的反射系數(shù)取決于波長(zhǎng)和沙波地形間隔的比例,當(dāng)兩者比例為2∶1時(shí)發(fā)生共振且反射系數(shù)最強(qiáng)。反射波的波幅隨沙波數(shù)量呈線性增長(zhǎng)趨勢(shì),且反射系數(shù)與沙波波幅和水深的比值呈正比。他們將整個(gè)共振波場(chǎng)劃分為兩塊,地形上方稱為“近域波場(chǎng)”,地形兩側(cè)稱為“遠(yuǎn)域波場(chǎng)”,在共振條件下,近域波場(chǎng)以及地形前方的遠(yuǎn)域波場(chǎng)會(huì)出現(xiàn)立波,將會(huì)影響地形前方的泥沙運(yùn)動(dòng),加速可侵蝕性海床的變化,這為沙波成長(zhǎng)機(jī)制的研究提供了科學(xué)依據(jù)。

    2 波浪布拉格共振的研究意義

    傳統(tǒng)的海岸防護(hù)工程以海堤直接阻擋波浪,消波作用明顯,但由于堤身直接承受全部的波浪荷載且堤前立波或不完全立波對(duì)堤腳沖刷劇烈,這些都可能導(dǎo)致工程有效使用壽命的縮短。傳統(tǒng)海堤在擋浪的同時(shí),還擋水擋沙,可能引起地形演變并阻礙人們的親水需求。20世紀(jì)80年代逐漸興起不露出水面的潛堤,既能反射部分波浪,還能允許部分水流和泥沙透過(guò)。然而由于其工作原理仍然是靠堤身反射和堤頂碎波,現(xiàn)有的潛堤堤頂高程仍然在平均低潮位附近,無(wú)法滿足通航要求,且堤腳沖刷仍然嚴(yán)重?;诓祭窆舱竦幕舅枷?可以通過(guò)設(shè)計(jì)堤頂高程低于平均低潮位的系列潛堤實(shí)現(xiàn)消浪效果。這種新型潛堤不影響海水交換和小型船舶通航,而且堤上容易長(zhǎng)出苔蘚,吸引浮游生物的聚集,對(duì)海岸生態(tài)環(huán)境有較好的維護(hù)作用。郭金棟等[4]探討了這類潛堤的應(yīng)用價(jià)值,認(rèn)為其不失為一種多效益兼顧的海岸防護(hù)工程。不少學(xué)者針對(duì)這類潛堤的結(jié)構(gòu)形式、尺寸和布設(shè)方案等開(kāi)展了研究,研究方法主要以波浪水槽試驗(yàn)和數(shù)值模擬方法為主,其中數(shù)值模擬方法以基于勢(shì)流理論的緩坡方程為主。Hsu等[5]、蔡立宏[6]、Wen等[7]、江鳴[8]分別研究了不同的潛堤形狀、海床坡度、堤高、堤寬、堤數(shù)等參數(shù)對(duì)共振反射系數(shù)的影響,Wen等[7]通過(guò)16000多組算例結(jié)果回歸出反射系數(shù)的經(jīng)驗(yàn)公式。羅恒[9]利用系列矩陣法改進(jìn)計(jì)算了三角形、梯形、半余弦形系列潛堤反射系數(shù)的解析解,并就長(zhǎng)波條件下潛堤的高度和寬度計(jì)算出最優(yōu)配置曲線。曾慧丹[10]在此基礎(chǔ)上補(bǔ)充了拋物線形潛堤的情形,在其研究中提及了全波譜的概念,但并非隨機(jī)波浪波譜,而是指頻率范圍較寬的單頻入射波。

    天然沙波上的波浪布拉格共振對(duì)河口海岸動(dòng)力地貌過(guò)程也有著重要影響。Elgar等[11]在Cape Cod Bay海域的實(shí)測(cè)資料顯示有約20%的入射波能被天然沿岸平行沙波反射,布拉格共振還會(huì)影響海岸及近海域的波譜演化,促進(jìn)穩(wěn)定的沿岸沙壩的形成。Renaud等[12]發(fā)現(xiàn)共振產(chǎn)生的立波會(huì)對(duì)船舶安全造成重要影響。Babcock等[13]認(rèn)為共振作用還可能引發(fā)背景水動(dòng)力場(chǎng)的脈動(dòng)噪音,可能引發(fā)深海區(qū)域的微震。我國(guó)各大江河口都分布著大面積的沙波群,楊世倫等[14]、李近元等[15]、孫杰等[16]分別對(duì)我國(guó)長(zhǎng)江口南支南港巷道、黃河口萊州灣以及珠江口內(nèi)伶仃島附近的沙波群的尺度和分布進(jìn)行了實(shí)測(cè)統(tǒng)計(jì),沙波波長(zhǎng)從微型沙波的2~3m到大型沙波的200~300m都有出現(xiàn),在浪、潮、流多動(dòng)力因素耦合影響下,復(fù)雜地形上的波浪布拉格共振對(duì)臨近水域的泥沙運(yùn)動(dòng)與地形演變都可能產(chǎn)生重要影響。

    3 波浪布拉格共振研究的主要進(jìn)展

    3.1 研究方法

    波浪布拉格共振現(xiàn)象的研究方法主要有現(xiàn)場(chǎng)觀測(cè)、理論解析、水槽試驗(yàn)和數(shù)值模擬等4種。

    a.現(xiàn)場(chǎng)觀測(cè)法。該方法是在大面積沿岸平行沙波群位置設(shè)立波浪測(cè)量點(diǎn),同時(shí)對(duì)海底沙波尺度進(jìn)行實(shí)測(cè),統(tǒng)計(jì)分析實(shí)測(cè)資料并計(jì)算波浪反射系數(shù)?,F(xiàn)場(chǎng)實(shí)測(cè)資料是理論解析、水槽試驗(yàn)和數(shù)值模擬方法得以順利開(kāi)展的必要基礎(chǔ)。Davies等[3]對(duì)布拉格共振的分析就是以1975年Short[17]針對(duì)沿岸短平行沙波的實(shí)測(cè)資料為基礎(chǔ)而展開(kāi)的。然而現(xiàn)場(chǎng)實(shí)測(cè)耗時(shí)耗力且具有一定的危險(xiǎn)性,同時(shí)需要大量資金或工程背景為基礎(chǔ)才能得以展開(kāi),我國(guó)尚缺乏對(duì)天然海域波浪布拉格共振的專門現(xiàn)場(chǎng)觀測(cè)。對(duì)于新型潛堤,工程建成后的波浪布拉格共振反射效果同樣需要現(xiàn)場(chǎng)實(shí)測(cè),并與相關(guān)遠(yuǎn)近海域的實(shí)測(cè)海浪譜進(jìn)行對(duì)比分析。

    b.理論解析法。該方法主要利用包括多重尺度展開(kāi)法在內(nèi)的攝動(dòng)方法對(duì)帶有特定邊界條件的非線性勢(shì)波方程進(jìn)行線性化,進(jìn)而推導(dǎo)主頻條件和非主頻條件下的反射系數(shù)公式,解釋波浪布拉格共振的現(xiàn)象及其原因,初步揭示基本物理機(jī)制,為思考并推導(dǎo)高階解,以及全面認(rèn)識(shí)布拉格共振現(xiàn)象做好鋪墊。理論解析方法不僅要求研究者具有嫻熟的數(shù)學(xué)推理能力,還要有清晰的物理概念,以便在推導(dǎo)過(guò)程中提出合理的假設(shè),進(jìn)行必要的簡(jiǎn)化?,F(xiàn)有針對(duì)波浪布拉格共振的理論解析研究,往往忽略了部分甚至全部的非線性項(xiàng),使得結(jié)論只能揭示較低階共振機(jī)制,難以解釋高階布拉格共振現(xiàn)象。

    c.水槽試驗(yàn)法。水槽試驗(yàn)法對(duì)波浪布拉格現(xiàn)象的實(shí)現(xiàn)和測(cè)量比現(xiàn)場(chǎng)觀測(cè)法經(jīng)濟(jì)方便,還可以排除真實(shí)海域里的干擾因素。由于布拉格共振研究方面實(shí)測(cè)資料較少,并且實(shí)測(cè)數(shù)據(jù)中非線性因素多而雜,并不適用于驗(yàn)證理論解析解的正確性,因而現(xiàn)有低階共振的理論解析結(jié)果一般都用水槽試驗(yàn)進(jìn)行驗(yàn)證。Heathershaw[18]通過(guò)水槽試驗(yàn)對(duì)Davies的第一類布拉格共振理論進(jìn)行了驗(yàn)證,并考察了共振作用對(duì)沙波成長(zhǎng)機(jī)制的影響。Guazzelli等[19]、Belzons 等[20]和Magne等[21]也分別對(duì)第二類和第三類布拉格共振的特性進(jìn)行過(guò)試驗(yàn)驗(yàn)證。然而,由于研究高階布拉格共振現(xiàn)象的物理試驗(yàn)對(duì)水槽和測(cè)量?jī)x器都有較高的要求,相關(guān)成果鮮有報(bào)道。

    d.數(shù)值模擬法。隨著現(xiàn)代計(jì)算機(jī)能力的提高和波浪模型數(shù)值算法的發(fā)展,利用數(shù)值模擬方法對(duì)各種波浪傳播過(guò)程進(jìn)行復(fù)演和預(yù)測(cè)已經(jīng)逐漸成為一種常用的研究手段,也成為高階非線性布拉格共振研究的主流?;谕卣沟木徠路匠踢\(yùn)算、連續(xù)應(yīng)用矩陣法、高階譜方法、VOF模型等都在布拉格共振的研究中得到應(yīng)用,模型的理論基礎(chǔ)及數(shù)值算法各有不同。以Dommermuth等[22-23]提出的高階譜方法為例,該方法不僅能夠分析大量子波的演變過(guò)程,還能考慮子波間任意高階的非線性相互作用,且能夠以指數(shù)收斂的效率模擬波浪和地形的相互作用。

    3.2 主要研究成果

    3.2.1 共振類型

    在1985年Heathershaw等[24]提出正向入射波通過(guò)單一周期系列正弦地形的布拉格共振理論解析公式的基礎(chǔ)上,Guazzelli等[19]、Liu等[25]都曾提出高階布拉格共振的共振條件。目前通常根據(jù)共振條件將波浪布拉格共振劃分為三類。

    第一類布拉格共振集中于三波二階相互作用。三波指入射表面波、反射表面波和一個(gè)單頻的地形波,二階是指分析結(jié)果限于表面波一階加地形波的一階,這種共振現(xiàn)象其實(shí)是線性共振。第二類布拉格共振的沙波地形由兩個(gè)不同頻率的正弦波組成,與表面入射波、表面反射波組成四波三階相互作用,三階是指自由表面波一階和地形波兩階。由于地形波的非線性影響,布拉格共振的頻帶會(huì)變寬,兩個(gè)地形波的差頻也會(huì)有相應(yīng)的次諧表面波反射。第三類布拉格共振同樣是四波三階相互作用,但四波是指三個(gè)表面波和一個(gè)地形波,三階是指表面波兩階加地形波一階。由于自由表面波成分的增加,將同時(shí)存在次諧波和超諧波,使布拉格共振特性更加復(fù)雜。三類布拉格共振的共振條件見(jiàn)表1。

    表1 三類布拉格共振的共振條件

    3.2.2 低階共振的反射系數(shù)公式

    Davies[2]利用正則攝動(dòng)法直接解析得到有限沙波數(shù)下的反射系數(shù)公式,但在無(wú)限個(gè)沙波時(shí)主頻的反射系數(shù)并不收斂。Mei[26]通過(guò)多重尺度展開(kāi)法對(duì)主頻反射系數(shù)公式進(jìn)行了修訂,Hara等[27]將理論拓展到二階,Mei等[28]利用該理論研究波浪斜向入射和斜坡底床情形下的低階共振機(jī)理。Miles[29]對(duì)拉普拉斯方程式積分,解析出等水深底床因微小高度變化所造成的反射系數(shù)公式。Rey等[30]基于拓展的多重尺度展開(kāi)法研究了第二類布拉格共振特性,就一階展開(kāi)和二階展開(kāi)的情形分別給出了解析公式,并討論了次諧波共振問(wèn)題。

    3.2.3 高階共振的數(shù)值模擬

    Kirby[31]基于緩坡方程模擬了“近域波場(chǎng)”中表面波的一般方程,同時(shí)考慮了波浪斜向入射、斜坡底床和水深變化等條件的影響,并將所得結(jié)論與Mei 等[28]的理論結(jié)果以及Davies等[3]的試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比,反射系數(shù)的結(jié)果基本一致,但模擬結(jié)果的主頻相對(duì)于理論值存在一定程度的下移量,并且與試驗(yàn)結(jié)果也不一致。Porter等[32]基于同樣的理論建立了三維模型,驗(yàn)證了布拉格共振在雙周期性起伏地形也會(huì)產(chǎn)生,發(fā)現(xiàn)一種給定模態(tài)的入射波會(huì)產(chǎn)生不符合Snell定律的不同模態(tài)的反射波,使得發(fā)生布拉格共振的條件更具多樣性。O’Hare等[33]將連續(xù)應(yīng)用矩陣法模型與Kirby模型進(jìn)行了對(duì)比,認(rèn)為連續(xù)應(yīng)用矩陣法模型運(yùn)算效率更高,但在主頻點(diǎn)和次主頻點(diǎn)的模擬上存在精度不高的問(wèn)題。Dalrymple 等[34]利用邊界積分法也做出模擬,所得結(jié)果和Mei 等[28]的理論結(jié)果以及Davies等[3]的試驗(yàn)數(shù)據(jù)相比同樣有較好的相似性。

    Liu等[35]將高階譜方法引入布拉格共振的數(shù)值計(jì)算中,得到反射系數(shù)的數(shù)值解。將波陡控制在線性理論可以接受的范圍內(nèi),通過(guò)改變波浪周期,以表面波波數(shù)與沙波波數(shù)的比值為變量,探討了布拉格共振的部分非線性特性。Alam等[36]以此為基礎(chǔ),探討了第三類布拉格共振特性,指出共振中將產(chǎn)生次諧波和超諧波,這些次生波被岸灘反射后加劇了近岸波浪的復(fù)雜程度。

    4 研究展望

    由于河口海岸水域波浪布拉格共振的復(fù)雜性及研究方法的局限性,仍然存在諸多亟待解決的問(wèn)題,主要包括共振主頻下移、不規(guī)則波布拉格共振特性、波流相互作用下布拉格共振機(jī)理等。

    4.1 波浪布拉格共振的共振主頻下移

    對(duì)于給定的周期性起伏地形,發(fā)生最強(qiáng)布拉格共振時(shí)相應(yīng)表面波的頻率稱為共振主頻,其精確計(jì)算是實(shí)際工程中應(yīng)用波浪布拉格共振達(dá)到最理想消浪效果的關(guān)鍵,所以理論研究主要集中在共振主頻的討論以及主頻反射系數(shù)的計(jì)算兩方面。迄今為止,研究只能得出線性規(guī)則波條件下波浪主頻反射系數(shù)的理論解析解。然而該解析理論依然離不開(kāi)Davies理論最初的假設(shè)[2],即線性規(guī)則波和有限個(gè)沙波地形滿足線性波浪理論的條件。大量的物理模型試驗(yàn)和數(shù)值模擬,如Heathershaw等[24]的物理模型試驗(yàn)結(jié)果與Liu等[25]基于高階譜方法的數(shù)值計(jì)算結(jié)果都表明,通過(guò)共振條件計(jì)算得出的布拉格共振主頻普遍存在著下移現(xiàn)象。Liu等[25]在數(shù)值模擬中增大沙波非線性后,發(fā)現(xiàn)主頻下移的現(xiàn)象更加明顯,并初步分析了該現(xiàn)象的產(chǎn)生機(jī)制,其將主頻下移的原因歸結(jié)于高階邊界效應(yīng)的影響,即諸如表面波波陡、沙波波陡、水深條件等因素不能嚴(yán)格滿足線性理論的假設(shè),使得共振波譜的主頻產(chǎn)生了下移。然而,這一結(jié)論還有待系統(tǒng)的試驗(yàn)數(shù)據(jù)的支持,且各種非線性因素的獨(dú)立影響乃至多因素的耦合作用將使共振主頻產(chǎn)生何種程度的下移尚無(wú)定論。

    4.2 不規(guī)則波的布拉格共振特性

    不規(guī)則波作用下波浪布拉格共振的復(fù)雜性,使其在理論解析和物理試驗(yàn)?zāi)M方面都難以開(kāi)展研究。Hsu等[37]提出了一套將Boussinesq方程應(yīng)用到布拉格共振中的數(shù)值模型。柯拓宇[38]以演進(jìn)型緩坡方程為基礎(chǔ),對(duì)不規(guī)則波作用下復(fù)合式矩形潛堤布拉格共振進(jìn)行了數(shù)值模擬,并提出了潛堤布置建議。林詹翰等[39]基于邊界積分法以JONSWAP譜為例數(shù)值模擬了梯形和矩形系列潛堤上的不規(guī)則波布拉格共振。然而,人們對(duì)不規(guī)則波布拉格共振機(jī)理的認(rèn)識(shí)尚不夠清晰,波譜特性、沙波形式等對(duì)布拉格共振的影響機(jī)制及影響程度都需要開(kāi)展大量的工作去研究解決。

    4.3 波流相互作用下的布拉格共振機(jī)理

    在河口海岸水域,存在著水流、岸灘反射波、海床蝕積等多種動(dòng)力因素,使得波浪布拉格共振的發(fā)生機(jī)理更加復(fù)雜。Kirby[40]基于多重尺度展開(kāi)理論,初步分析了同向水流對(duì)布拉格共振的影響,指出同向恒定流會(huì)使共振主頻產(chǎn)生下移,峰值變大。Magne等[21]進(jìn)行了波流共同作用的布拉格共振相關(guān)試驗(yàn),但試驗(yàn)僅為印證Kirby的理論結(jié)果而展開(kāi)。Yu等[41]認(rèn)為岸灘反射波反向通過(guò)沙波地形時(shí),若滿足一定條件,反而有可能使沙波地形前的波浪得到增強(qiáng)。Belibassakisa等[42-43]提出了模式耦合模型研究波浪、水流、地形三者相互作用,可以較好地模擬二維地形和簡(jiǎn)單沿岸剪切流條件下的波浪布拉格共振,并提出了緩坡剪切方程。但對(duì)于三維條件、復(fù)雜水流條件、非線性影響的波流相互作用下的布拉格共振數(shù)值模擬還有待進(jìn)一步研究,已有的模型模擬結(jié)果尚缺乏試驗(yàn)資料的支持。目前對(duì)于波浪、水流和沙波地形共存的布拉格共振機(jī)制的認(rèn)識(shí)還處于初級(jí)階段,系統(tǒng)性的試驗(yàn)證明和理論分析都亟待展開(kāi)。

    5 結(jié) 語(yǔ)

    本文評(píng)述了關(guān)于波浪布拉格共振研究的現(xiàn)場(chǎng)觀測(cè)、理論研究、水槽試驗(yàn)、數(shù)值模擬等4種方法的優(yōu)缺點(diǎn),總結(jié)了波浪布拉格共振類型、低階共振反射系數(shù)公式和高階共振數(shù)值模擬方面的主要研究成果。認(rèn)為現(xiàn)場(chǎng)實(shí)測(cè)資料是開(kāi)展理論研究、水槽試驗(yàn)和數(shù)值模擬研究的基礎(chǔ),由于現(xiàn)場(chǎng)觀測(cè)技術(shù)難度較大、成本昂貴且具有一定的危險(xiǎn)性,相關(guān)研究成果較少;理論解析法以多重尺度展開(kāi)法為基礎(chǔ),低階線性布拉格共振方面的理論解析研究已趨于完善;水槽試驗(yàn)為低階共振的理論解析研究提供驗(yàn)證數(shù)據(jù),但二者在高階共振方面的研究成果不多;數(shù)值模擬方法具有高效運(yùn)算能力,在高階布拉格共振研究中得到廣泛應(yīng)用,研究成果較多,然而對(duì)于布拉格共振主頻下移等非線性機(jī)制問(wèn)題的研究還有待深入細(xì)致地開(kāi)展。

    針對(duì)河口海岸水域波浪布拉格共振的復(fù)雜性,急需深入研究波浪布拉格共振主頻下移、不規(guī)則波布拉格共振特性、波流相互作用下布拉格共振機(jī)理等3個(gè)問(wèn)題,為其在河口海岸開(kāi)發(fā)與保護(hù)的工程應(yīng)用提供科學(xué)依據(jù)。

    參考文獻(xiàn):

    [1] LONG R R.Scattering of surface waves by an irregular bottom[J].Journal ofgeophysical Research, 1973, 78: 7861-7870.

    [2] DAVIES Ag.On the interaction between surface waves and undulations on the seabed [J].Journal ofmarine Research, 1982, 40(2): 331-368.

    [3] DAVIES Ag, HEATHERSHAW A D.Surface-wave propagation over sinusoidally varying topography [J].Journal of Fluidmechanics, 1984, 144: 419-443.

    [4]郭金棟,陳文俊,陳國(guó)書(shū).雙列潛堤對(duì)海灘防治效益之研究[C] //第二十一屆海洋工程研討會(huì)論文集.北京:海洋出版社, 1999: 307-313.

    [5] HSU T, CHANG H, TSAI L.Bragg reflection of waves by different shapes of artificial bars [J].China Ocean Engineering, 2002, 16: 343-358.

    [6]蔡立宏.波浪通過(guò)系列潛堤之布拉格反射研究[D].臺(tái)南:成功大學(xué), 2003.

    [7] WEN C, TSAI L.Numerical simulation of Bragg reflection based on linear waves propagation over a series of rectangular seabed[J].China Ocean Engineering, 2008, 22(1): 71-86.

    [8]江鳴.波浪通過(guò)系列矩形潛堤的數(shù)值模擬[D].天津:天津大學(xué), 2001.

    [9]羅恒.長(zhǎng)波條件下三角形、梯形和半余弦布拉格潛堤的優(yōu)化配置[D].南寧:廣西民族大學(xué), 2013.

    [10]曾慧丹.全波譜條件下三類布拉格防波堤的最優(yōu)配置[D].南寧:廣西民族大學(xué), 2014.

    [11] ELGAR S, RAUBENHEIMER B, HERBERS T H C.Bragg reflection of ocean waves from sandbars [J].Geophysical Research Letters, 2003, 30(1): 161-164.

    [12] RENAUDm, REZENDE F, WAALS O, et al.Secondorder wave loads on a LNG carrier inmulti-directional waves[C] //Proceedings of 27th International Conference on Offshoremechanics and Arctic Engineering.Estoril: ASME, 2008: 363-370.

    [13] BABCOCK Jm, KIRKENDALL B A, ORCUTT J A.Relationships between ocean bottom noise and the environment[J].Bulletin of the Seismological Society of America, 1994, 84: 1991-2007.

    [14]楊世倫,張正惕,謝文輝,等.長(zhǎng)江口南港航道沙波群研究[J].海洋工程, 1999, 17(2):79-88.(YANG Shilun, ZHANG Zhengti, XIE Wenhui, et al.A study of sandwaves in the south channel of the Yangtze Estuary [J].The Ocean Engineering, 1999, 17(2):79-88.(in Chinese))

    [15]李近元,范奉鑫,徐濤,等.萊州灣東部沙波地貌分布特征及其形成演化[J].海洋科學(xué), 2011,35(7): 51-54.(LI Jinyuan, FAN Fengxin, XU Tao, et al.The distribution characteristics and evolution of the sand wave in the east of Laizhou Bay[J].Marine Sciences, 2011, 35(7): 51-54.(in Chinese))

    [16]孫杰,詹文歡,賈建業(yè),等.珠江口海域?yàn)?zāi)害地質(zhì)因素及其與環(huán)境變化的關(guān)系[J].熱帶海洋學(xué)報(bào), 2010,29 (1): 104-110.(SUN Jie, ZHAN Wenhuan, JIA Jianye, et al.Hazardousgeology and its relationship with environmental evolution in the Pearl River Estuary[J].Journal of Tropical Oceanography, 2010,29(1): 104-110.(in Chinese))

    [17] SHORT A D.Multiple offshore bars and standing waves [J].Journal ofgeophysical Research, 1975, 80: 3838-3840.

    [18] HEATHERSHAW A D.Seabed-wave resonance and sand bargrowth[J].Nature, 1982, 296(25): 343-345.

    [19]gUAZZELLI E, REY V, BELZONSm.Higher-order Bragg reflection ofgravity surface waves by periodic beds [J].Journal of Fluidmechanics, 1992, 245: 301-317.

    [20] BELZONSm, REY V,gUAZZELLI E.Subharmonic Bragg resonance for surface water waves[J].Europhysics Letters, 1991, 16(2): 189-194.

    [21]mAGNE R, REY V.Measurement of wave scattering by topography in the presence of currents[J].Physics of Fluid, 2005, 17: 1-8.

    [22] DOMMERMUTH Dg, YUE D K P.A higher-order spectralmethod for the study of nonlineargravity waves [J].Journal of Fluidmechanics, 1987,184: 267-288.

    [23] DOMMERMUTH Dg, YUE D K P.The nonlinear threedimensional wavesgenerated by amoving surface disturbance[M].Washington D.C.: National Academy Press, 1988: 103-129.

    [24] HEATHERSHAW A D, DAVIES Ag.Resonant wave reflection by transverse bedforms and its relation to beaches and offshore bars[J].Marinegeology, 1985, 62(3): 321-338.

    [25] LIU Y, YUE D K P.Ongeneralized Bragg scattering of surface waves by bottom ripples [J].Journal of Fluidmechanics, 1998, 356: 297-326.

    [26]mEI C C.Resonant reflection of surface water waves by periodic sandbars [J].Journal of Fluidmechanics, 1985, 152: 315-335.

    [27] HARA T,mEI C C.Bragg scattering of surface waves by periodic bars: theory and experiment [J].Journal of Fluidmechanics, 1987, 178: 221-241.

    [28]mEI C C, HARA T, NACIRIm.Note on Bragg scattering of water waves by parallel bars on the seabed[J].Journal of Fluidmechanics, 1986, 186: 147-162.

    [29]mILES J W.Weakly nonlinear Kelvin-Helmholtz waves [J].Journal of Fluidmechanics, 1986, 172: 513-529.[30] REY V,gUAZZELLI E,mEI C C.Resonant reflection of surfacegravity waves by one dimensional doubly sinusoidal beds[J].Physics of Fluids, 1996, 8(6): 1525-1530.

    [31] KIRBY J T.Ageneral wave equation for waves over rippled beds [J].Journal of Fluidmechanics, 1986, 162: 171-186.

    [32] PORTER R, PORTER D.Interaction of water waves with three-dimensional periodic topography [J].Journal of Fluidmechanics, 2001, 434: 301-335.

    [33] O’HARE T J, DAVIES Ag.A comparison of twomodels for surface-wave propagation over rapidly varying [J].Applied Ocean Research, 1993, 15(1): 1-11.

    [34] DALRYMPLE R A, KIRBY J T.Water waves over ripples[J].Journal of Waterway, Port, Coastal, and Ocean Engineering, 1986, 112(2): 309-319.

    [35] LIU Y, DOMMERMUTH Dg, YUE D K P.A high-order spectralmethhod for nonlinear wave-body interactions [J].Journal of Fluidmechanics, 1992, 245: 115-136.

    [36] ALAMm R, LIU Y, YUE D K P.Oblique sub-and super-harmonic Bragg resonance of surface waves by bottom ripples[J].Journal of Fluidmechanics, 2010, 643: 437-447.

    [37] HSU T, HSIAO S, OU S, et al.An application of Boussinesq equation to Bragg reflection of irregular waves [J].China Ocean Engineering, 2007, 34: 870-883.

    [38]柯拓宇.不規(guī)則波浪作用下之布拉格共振[D].臺(tái)南:成功大學(xué), 2007.

    [39]林詹翰,唐宏結(jié),黃材成.不規(guī)則波通過(guò)系列潛堤之布拉格共振數(shù)值研究[C] //第十屆水下技術(shù)研討會(huì)暨國(guó)科會(huì)成果發(fā)表會(huì).高雄: [s.n.], 2008.

    [40] KIRBY J T.Current effects on resonant reflection of surface water waves by sand bars[J].Journal of Fluidmechanics,1988, 186: 501-520.

    [41] YU J,mEI C C.Do longshore bars shelter the shore[J].Journal of Fluidmechanics, 2000, 404: 251-268.

    [42] BELIBASSAKISA K A.A coupled-modemodel for the scattering of water waves by shearing currents in variable bathymetry[J].Journal of Fluidmechanics, 2007, 578: 413-434.

    [43] BELIBASSAKISA K A,gEROSTATHISB T P, ATHANASSOULISAg A.A coupled-modemodel for water wave scattering by horizontal, non-homogeneous current ingeneral bottom topography[J].Applied Ocean Research, 2011, 33: 384-397.

    中圖分類號(hào):P731.22

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1006- 7647(2016)03- 0083- 05

    DOI:10.3880/j.issn.1006- 7647.2016.03.017

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(51379071, 51137002),中央高??蒲袠I(yè)務(wù)費(fèi)項(xiàng)目(2014B04114)

    作者簡(jiǎn)介:鄭金海(1972—),男,教授,主要從事海岸動(dòng)力學(xué)與海岸工程研究。E-mail: jhzheng@ hhu.edu.cn

    收稿日期:(2015 08- 31 編輯:駱超)

    Research progress in wave Bragg resonance

    ZHENG Jinhai1,2, YU Haofeng1,2, TAO Aifeng1,2, FAN Jun1,2, WANGYi(1.Key Laboratory of Coastal Disaster and Defence(Hohai University),ministry of Education, Nanjing 210098, China;2.College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, China)

    1,2

    Abstract:Bragg resonance phenomena, occurring as waves propagating from the sea to shallow water under particular topographic conditions, were studied to provide a scientific basis for the design of estuarine and coastal protection engineering and the forecast ofmorphodynamic evolution.The definition and characteristics of wave Bragg resonance as well as the significance of research are introduced, and the pros and cons of four approaches for study of the phenomena are reviewed, including field observation, theoretical analysis, numerical simulation, and experimental investigation.Major achievements regarding the types of wave Bragg resonance, the formulas of the low-order resonance reflection coefficient, and high-order resonance numerical simulation are summarized.Three important topics for further intensive study are proposed, including the negative drift of dominant frequency, the characteristics of the irregular wave Bragg resonance, and thegenerationmechanism of Bragg resonance under wave-current interactions.

    Key words:wave;Bragg resonance;nonlinear resonance;negative drift of dominant frequency;wave-current interaction;review

    猜你喜歡
    述評(píng)波浪
    波浪谷和波浪巖
    波浪谷隨想
    去看神奇波浪谷
    走,到波浪巖去“沖浪”
    關(guān)于法學(xué)視角下的服務(wù)型政府建設(shè)探析
    東方教育(2016年3期)2016-12-14 20:36:15
    近十年來(lái)中國(guó)城鎮(zhèn)化發(fā)展的研究述評(píng)
    現(xiàn)代朱子哲學(xué)研究述評(píng)
    桂海論叢(2016年4期)2016-12-09 11:17:21
    國(guó)內(nèi)三十年語(yǔ)碼轉(zhuǎn)換研究述評(píng)
    文教資料(2016年22期)2016-11-28 13:32:22
    詩(shī)歌里的低訴,蒼涼中的守望
    人間(2016年28期)2016-11-10 21:51:46
    國(guó)內(nèi)科技投入研究現(xiàn)狀簡(jiǎn)評(píng)
    中文字幕av在线有码专区| 成人av一区二区三区在线看| 一区二区三区高清视频在线| 中文字幕av在线有码专区| 真实男女啪啪啪动态图| 麻豆av噜噜一区二区三区| 简卡轻食公司| 午夜免费激情av| 最近在线观看免费完整版| 老司机影院成人| 国产淫片久久久久久久久| 精品无人区乱码1区二区| 国产精品美女特级片免费视频播放器| 一a级毛片在线观看| 村上凉子中文字幕在线| 国产一区二区三区av在线 | 亚洲国产精品合色在线| 国产精品一区二区三区四区免费观看 | 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 精品人妻偷拍中文字幕| 波多野结衣高清无吗| 精品久久久久久久久亚洲| 免费av毛片视频| 色综合亚洲欧美另类图片| 一a级毛片在线观看| 自拍偷自拍亚洲精品老妇| 午夜影院日韩av| 天天躁夜夜躁狠狠久久av| 一进一出抽搐gif免费好疼| 18禁黄网站禁片免费观看直播| 欧美精品国产亚洲| 三级毛片av免费| 国产精品伦人一区二区| av中文乱码字幕在线| 午夜福利高清视频| 久久精品91蜜桃| 亚洲经典国产精华液单| 99视频精品全部免费 在线| 欧美一区二区精品小视频在线| 高清午夜精品一区二区三区 | 欧美bdsm另类| 丝袜美腿在线中文| 久久久国产成人精品二区| 日日摸夜夜添夜夜添av毛片| 少妇人妻一区二区三区视频| 毛片一级片免费看久久久久| 中文字幕人妻熟人妻熟丝袜美| 欧美人与善性xxx| 亚洲成人久久性| 亚洲av成人av| 成人性生交大片免费视频hd| 亚洲欧美成人精品一区二区| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 中文字幕免费在线视频6| 大香蕉久久网| 久久精品国产清高在天天线| 麻豆成人午夜福利视频| 国产伦在线观看视频一区| 自拍偷自拍亚洲精品老妇| 18禁在线播放成人免费| 亚洲美女搞黄在线观看 | 国产精品久久视频播放| 成人美女网站在线观看视频| 熟妇人妻久久中文字幕3abv| 国产精品99久久久久久久久| 丝袜美腿在线中文| 久久久久国内视频| 成人鲁丝片一二三区免费| 欧美成人免费av一区二区三区| 别揉我奶头~嗯~啊~动态视频| 午夜激情欧美在线| 免费看光身美女| 少妇被粗大猛烈的视频| 午夜免费激情av| 日本一本二区三区精品| 蜜臀久久99精品久久宅男| 亚洲最大成人av| 少妇熟女欧美另类| 免费无遮挡裸体视频| 久久亚洲精品不卡| 日本欧美国产在线视频| 精品国内亚洲2022精品成人| 少妇人妻一区二区三区视频| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区| 免费无遮挡裸体视频| 一进一出好大好爽视频| 久99久视频精品免费| 赤兔流量卡办理| av中文乱码字幕在线| 嫩草影院新地址| 欧美另类亚洲清纯唯美| 身体一侧抽搐| 嫩草影院精品99| 国产69精品久久久久777片| 中文资源天堂在线| 深爱激情五月婷婷| 亚洲av五月六月丁香网| 国产亚洲精品久久久久久毛片| 一个人看的www免费观看视频| 亚洲电影在线观看av| 国产一区二区在线观看日韩| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验 | 美女大奶头视频| 校园人妻丝袜中文字幕| 亚洲国产精品久久男人天堂| 你懂的网址亚洲精品在线观看 | 日韩精品有码人妻一区| 男插女下体视频免费在线播放| 国产精品永久免费网站| 色哟哟哟哟哟哟| 久久久久久国产a免费观看| 欧美色欧美亚洲另类二区| 搞女人的毛片| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 伦理电影大哥的女人| 免费电影在线观看免费观看| 久久精品国产亚洲网站| av天堂中文字幕网| 国产男人的电影天堂91| 久久草成人影院| 一个人看的www免费观看视频| 国产成人91sexporn| 免费观看的影片在线观看| 联通29元200g的流量卡| 国产一区二区在线观看日韩| av女优亚洲男人天堂| 亚洲成av人片在线播放无| 18禁黄网站禁片免费观看直播| 久久6这里有精品| 尤物成人国产欧美一区二区三区| 国产免费一级a男人的天堂| 91久久精品国产一区二区三区| 在线免费观看不下载黄p国产| 午夜日韩欧美国产| 婷婷六月久久综合丁香| av女优亚洲男人天堂| 久久久久免费精品人妻一区二区| 男女做爰动态图高潮gif福利片| 亚洲婷婷狠狠爱综合网| 国产精品不卡视频一区二区| 久久精品国产清高在天天线| 男人的好看免费观看在线视频| 欧美成人精品欧美一级黄| 免费在线观看成人毛片| 最近中文字幕高清免费大全6| 国产精品精品国产色婷婷| 久久午夜亚洲精品久久| 尤物成人国产欧美一区二区三区| 天天一区二区日本电影三级| 免费看a级黄色片| a级毛片免费高清观看在线播放| 日本一二三区视频观看| 一本精品99久久精品77| 禁无遮挡网站| 99久国产av精品| 在线播放无遮挡| 成人av一区二区三区在线看| 日韩三级伦理在线观看| 日韩欧美精品免费久久| 国产精华一区二区三区| 亚洲av电影不卡..在线观看| 亚洲精华国产精华液的使用体验 | 国内精品一区二区在线观看| 久久久久久久久久黄片| 女人被狂操c到高潮| .国产精品久久| 欧美精品国产亚洲| 亚洲美女搞黄在线观看 | a级毛色黄片| 波野结衣二区三区在线| 噜噜噜噜噜久久久久久91| 人妻久久中文字幕网| 亚洲国产高清在线一区二区三| 波多野结衣高清作品| 亚洲一区高清亚洲精品| 亚洲一区二区三区色噜噜| 国产精品精品国产色婷婷| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 免费看美女性在线毛片视频| 亚洲人成网站高清观看| 久久精品久久久久久噜噜老黄 | 91久久精品国产一区二区成人| 老师上课跳d突然被开到最大视频| 日本黄大片高清| 美女xxoo啪啪120秒动态图| 国产乱人偷精品视频| 欧美三级亚洲精品| 一级毛片久久久久久久久女| 中文字幕精品亚洲无线码一区| 午夜精品国产一区二区电影 | 国产真实伦视频高清在线观看| 在线观看av片永久免费下载| 国产av在哪里看| 1000部很黄的大片| 搡老岳熟女国产| 国产不卡一卡二| 在线观看一区二区三区| 久久精品综合一区二区三区| 一级毛片电影观看 | 不卡视频在线观看欧美| 成人永久免费在线观看视频| 91狼人影院| 啦啦啦观看免费观看视频高清| 久久久久久久久久成人| 少妇的逼水好多| 内射极品少妇av片p| 国产精品综合久久久久久久免费| 看免费成人av毛片| 少妇高潮的动态图| 亚洲经典国产精华液单| 亚洲国产精品sss在线观看| 成人特级av手机在线观看| 欧美成人a在线观看| 高清毛片免费看| 老熟妇乱子伦视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲18禁久久av| 男女视频在线观看网站免费| 热99re8久久精品国产| 真实男女啪啪啪动态图| 国产综合懂色| 99久久九九国产精品国产免费| 色综合色国产| 久久这里只有精品中国| 亚洲七黄色美女视频| 美女 人体艺术 gogo| 免费观看精品视频网站| 老熟妇乱子伦视频在线观看| av.在线天堂| 国产av在哪里看| 亚洲美女视频黄频| 一个人免费在线观看电影| av卡一久久| 午夜福利视频1000在线观看| 可以在线观看的亚洲视频| 国产精品美女特级片免费视频播放器| 久久久久精品国产欧美久久久| 中文资源天堂在线| 欧美又色又爽又黄视频| 91久久精品国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 如何舔出高潮| 精品福利观看| 99久久精品热视频| 给我免费播放毛片高清在线观看| 久久久欧美国产精品| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 成人av一区二区三区在线看| 99久久精品国产国产毛片| 人妻制服诱惑在线中文字幕| 最近视频中文字幕2019在线8| 尤物成人国产欧美一区二区三区| 国产精品无大码| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人中文| 99九九线精品视频在线观看视频| 亚洲欧美日韩高清在线视频| 久久久久久伊人网av| 麻豆av噜噜一区二区三区| 亚洲成人久久性| 97热精品久久久久久| 日本熟妇午夜| 日韩欧美免费精品| 在线国产一区二区在线| 九色成人免费人妻av| 国产av一区在线观看免费| av在线观看视频网站免费| 99久国产av精品国产电影| 网址你懂的国产日韩在线| 欧美性感艳星| 国产av麻豆久久久久久久| 日日撸夜夜添| 亚洲图色成人| 男人舔女人下体高潮全视频| 久久久国产成人精品二区| 久久热精品热| 99热网站在线观看| 午夜精品国产一区二区电影 | 中国美女看黄片| 久久精品夜色国产| 国产探花在线观看一区二区| 色5月婷婷丁香| 久久精品91蜜桃| 18禁黄网站禁片免费观看直播| 波多野结衣高清作品| 亚洲av一区综合| 最近手机中文字幕大全| 精品久久久久久久久久久久久| 中文字幕久久专区| 亚洲精品在线观看二区| 免费看美女性在线毛片视频| 老司机午夜福利在线观看视频| 国产伦在线观看视频一区| 男女下面进入的视频免费午夜| 18禁裸乳无遮挡免费网站照片| 天堂av国产一区二区熟女人妻| 色综合亚洲欧美另类图片| 晚上一个人看的免费电影| 自拍偷自拍亚洲精品老妇| 国产美女午夜福利| 成人毛片a级毛片在线播放| 99riav亚洲国产免费| 色综合色国产| 免费搜索国产男女视频| 伦精品一区二区三区| 天堂√8在线中文| av在线老鸭窝| 国产精品人妻久久久影院| 精品久久久久久久久久久久久| 国产真实伦视频高清在线观看| 亚洲国产精品成人久久小说 | 精品一区二区免费观看| 校园春色视频在线观看| 看非洲黑人一级黄片| 干丝袜人妻中文字幕| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 一区二区三区四区激情视频 | 一进一出抽搐gif免费好疼| 偷拍熟女少妇极品色| 成人鲁丝片一二三区免费| .国产精品久久| 色在线成人网| 少妇裸体淫交视频免费看高清| 亚洲国产精品成人综合色| 夜夜夜夜夜久久久久| 啦啦啦韩国在线观看视频| 亚洲第一区二区三区不卡| 亚洲自拍偷在线| 日韩欧美三级三区| 美女xxoo啪啪120秒动态图| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 国产精品人妻久久久影院| 夜夜夜夜夜久久久久| 国产老妇女一区| a级毛片a级免费在线| 人妻久久中文字幕网| 国产乱人视频| 男人的好看免费观看在线视频| 国内精品一区二区在线观看| 卡戴珊不雅视频在线播放| 黄色配什么色好看| 色播亚洲综合网| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频 | 精品人妻视频免费看| 国产中年淑女户外野战色| 亚洲精品亚洲一区二区| 97超级碰碰碰精品色视频在线观看| 国产高潮美女av| 国产片特级美女逼逼视频| 亚洲第一区二区三区不卡| 99久国产av精品国产电影| 精品一区二区三区视频在线观看免费| 欧美成人一区二区免费高清观看| 国产黄片美女视频| 插阴视频在线观看视频| 色哟哟哟哟哟哟| 一本一本综合久久| 一区二区三区高清视频在线| 非洲黑人性xxxx精品又粗又长| 国产高潮美女av| 成人性生交大片免费视频hd| 神马国产精品三级电影在线观看| 在线免费观看不下载黄p国产| 亚洲成人久久性| 99在线视频只有这里精品首页| av.在线天堂| 久久国内精品自在自线图片| 99久久久亚洲精品蜜臀av| 好男人在线观看高清免费视频| av.在线天堂| 亚洲av免费高清在线观看| 色哟哟哟哟哟哟| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 好男人在线观看高清免费视频| 在线天堂最新版资源| 禁无遮挡网站| 久久热精品热| 麻豆av噜噜一区二区三区| 亚洲欧美日韩无卡精品| 如何舔出高潮| 日韩欧美一区二区三区在线观看| 在线播放国产精品三级| 99热只有精品国产| 免费一级毛片在线播放高清视频| 麻豆国产97在线/欧美| 亚洲丝袜综合中文字幕| 狂野欧美激情性xxxx在线观看| 成人一区二区视频在线观看| 一级毛片电影观看 | 少妇熟女欧美另类| 成人三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 高清日韩中文字幕在线| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 久久久久久久久久成人| av卡一久久| 国产亚洲精品久久久久久毛片| 亚洲精品久久国产高清桃花| 精品一区二区三区视频在线观看免费| 国产精品嫩草影院av在线观看| 精品一区二区三区视频在线观看免费| 国内精品宾馆在线| 日韩av不卡免费在线播放| 国产日本99.免费观看| 成人三级黄色视频| 国产精品永久免费网站| 国产欧美日韩精品亚洲av| 精品久久久久久久人妻蜜臀av| 国国产精品蜜臀av免费| 亚洲,欧美,日韩| 亚洲七黄色美女视频| 97超碰精品成人国产| 亚洲熟妇中文字幕五十中出| 99热网站在线观看| 精品欧美国产一区二区三| 黄色一级大片看看| 久久久久久大精品| 极品教师在线视频| 日韩一本色道免费dvd| 一本久久中文字幕| 亚洲av五月六月丁香网| 欧美一区二区亚洲| 欧美日本视频| 亚洲aⅴ乱码一区二区在线播放| 91狼人影院| 欧美在线一区亚洲| 国产极品精品免费视频能看的| 午夜精品在线福利| 亚洲精品日韩在线中文字幕 | 在线播放国产精品三级| 久久人人爽人人爽人人片va| 国产美女午夜福利| 一进一出抽搐gif免费好疼| 中文字幕精品亚洲无线码一区| 亚洲激情五月婷婷啪啪| 啦啦啦啦在线视频资源| 午夜a级毛片| 中文在线观看免费www的网站| 精品乱码久久久久久99久播| 免费人成在线观看视频色| 中文字幕熟女人妻在线| 国产精品乱码一区二三区的特点| 免费人成视频x8x8入口观看| 亚洲av二区三区四区| 99久国产av精品国产电影| 此物有八面人人有两片| 晚上一个人看的免费电影| 香蕉av资源在线| 国产色婷婷99| 亚洲av一区综合| 亚洲专区国产一区二区| 日本欧美国产在线视频| av天堂在线播放| 欧美另类亚洲清纯唯美| 久久久久久久久久黄片| 久久久欧美国产精品| 久久6这里有精品| 亚洲电影在线观看av| 少妇人妻精品综合一区二区 | 波多野结衣高清无吗| 一进一出抽搐动态| 日韩 亚洲 欧美在线| 国产毛片a区久久久久| 俄罗斯特黄特色一大片| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 亚洲天堂国产精品一区在线| 亚洲av.av天堂| 午夜福利成人在线免费观看| 亚洲真实伦在线观看| 国产视频内射| 嫩草影视91久久| 亚洲欧美中文字幕日韩二区| 亚洲电影在线观看av| 精品欧美国产一区二区三| 亚洲电影在线观看av| 麻豆国产av国片精品| 身体一侧抽搐| 99久国产av精品国产电影| 日本黄色视频三级网站网址| 高清毛片免费观看视频网站| 欧美最新免费一区二区三区| 欧美zozozo另类| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 亚洲国产精品合色在线| 两个人视频免费观看高清| 国产免费男女视频| 久久久久久国产a免费观看| 日韩成人伦理影院| 我要搜黄色片| 亚洲av中文字字幕乱码综合| 久久久久久久午夜电影| 久久久a久久爽久久v久久| 香蕉av资源在线| 成人综合一区亚洲| 国产欧美日韩精品亚洲av| 国内精品久久久久精免费| 日韩欧美三级三区| 久久久a久久爽久久v久久| 日日撸夜夜添| 亚洲三级黄色毛片| 亚洲18禁久久av| 免费一级毛片在线播放高清视频| av天堂在线播放| 国产69精品久久久久777片| 亚洲成人精品中文字幕电影| 女人被狂操c到高潮| 免费看a级黄色片| 久久久国产成人精品二区| 村上凉子中文字幕在线| 久久久国产成人精品二区| 黄片wwwwww| 亚洲美女搞黄在线观看 | 99九九线精品视频在线观看视频| 欧美最黄视频在线播放免费| 亚洲av.av天堂| 国产色爽女视频免费观看| 亚洲国产欧洲综合997久久,| 久久久久久久久中文| 亚洲av免费在线观看| 国产精品伦人一区二区| 亚洲精品一区av在线观看| 国产精品一二三区在线看| 99视频精品全部免费 在线| 国产一区二区亚洲精品在线观看| 久久人妻av系列| 国产精品,欧美在线| 久久人妻av系列| 国产一区亚洲一区在线观看| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 搞女人的毛片| 亚洲av熟女| 国语自产精品视频在线第100页| 一个人免费在线观看电影| 白带黄色成豆腐渣| 禁无遮挡网站| 亚洲av中文av极速乱| 麻豆精品久久久久久蜜桃| 国产精品无大码| 精品久久久久久成人av| 在线观看一区二区三区| 欧美成人免费av一区二区三区| 久久精品国产亚洲av香蕉五月| 蜜桃亚洲精品一区二区三区| 国产精品av视频在线免费观看| 高清毛片免费看| av免费在线看不卡| 联通29元200g的流量卡| 国产一区二区亚洲精品在线观看| 日韩大尺度精品在线看网址| 精品乱码久久久久久99久播| 日韩欧美国产在线观看| 国产男人的电影天堂91| 亚洲av中文字字幕乱码综合| 91在线精品国自产拍蜜月| 99久久久亚洲精品蜜臀av| 蜜桃久久精品国产亚洲av| or卡值多少钱| 国产亚洲精品久久久久久毛片| 麻豆乱淫一区二区| 国产精品人妻久久久影院| 看免费成人av毛片| 三级毛片av免费| 久久久久久伊人网av| 亚洲,欧美,日韩| 国产淫片久久久久久久久| 国产精品av视频在线免费观看| 国产激情偷乱视频一区二区| 久久久久国内视频| 国产精品无大码| 国产精品电影一区二区三区| 一区二区三区免费毛片| av中文乱码字幕在线| 亚洲在线观看片| 久久韩国三级中文字幕| 婷婷精品国产亚洲av在线| 99热这里只有是精品50| 亚洲国产精品成人综合色| av天堂在线播放| 黄色欧美视频在线观看| 久久韩国三级中文字幕| av在线观看视频网站免费| 国产色婷婷99| 久久天躁狠狠躁夜夜2o2o| 亚洲成人精品中文字幕电影| 可以在线观看毛片的网站| 亚洲精品久久国产高清桃花| 成人三级黄色视频| 精品久久久久久久人妻蜜臀av| 亚洲婷婷狠狠爱综合网| 国产一区二区三区av在线 | 全区人妻精品视频| 亚洲av不卡在线观看| or卡值多少钱| 精品一区二区三区人妻视频| 亚洲精品乱码久久久v下载方式| 丰满乱子伦码专区| 精品一区二区三区视频在线观看免费| 国产单亲对白刺激| 女生性感内裤真人,穿戴方法视频|