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

    深水聚焦波群的雙層Boussinesq數(shù)值模擬

    2022-09-21 03:40:58劉必勁張振偉劉忠波
    船舶力學(xué) 2022年9期
    關(guān)鍵詞:波面波峰雙層

    劉必勁,張振偉,劉忠波

    (1.浙江大學(xué)海南研究院,海南三亞 572025;2.廈門理工學(xué)院土木工程與建筑學(xué)院,福建廈門 361024;3.大連海事大學(xué)交通運(yùn)輸工程學(xué)院,遼寧大連 116024)

    0 引 言

    畸形波是非線性極強(qiáng)的瞬態(tài)大浪,受實(shí)際海洋中各種不確定的眾多環(huán)境因素影響,導(dǎo)致人類對(duì)這種畸形波的產(chǎn)生機(jī)制并沒有完全掌握。當(dāng)這種瞬態(tài)大浪作用在過往的船只或海上建筑物時(shí),將會(huì)嚴(yán)重影響它們的安全。為了探求畸形波的生成機(jī)制,國內(nèi)外學(xué)者嘗試從相位聚焦、波-流非線性相互作用、地形聚焦作用、表面風(fēng)的作用、波浪調(diào)制不穩(wěn)定性等方面開展畸形波的研究。

    對(duì)畸形波的研究主要方法有物理模型實(shí)驗(yàn)室模擬[1-5]、數(shù)值模擬[6]和理論分析[7-8]等。近年來,伴隨計(jì)算機(jī)性能的大幅度提升和各類數(shù)值模型的不斷完善,采用數(shù)值模型開展畸形波的研究成為比較常見且主要的手段,其中時(shí)域計(jì)算模型多采用勢流理論、歐拉方程或粘性流理論等。在勢流理論方面,有高階邊界元求解Laplace 方程的數(shù)值模型[6]、高階譜數(shù)值模型[9-11]和Boussinesq 型水波數(shù)值模型[12-13]。歐拉方程或粘性Navier-Stokes方程基礎(chǔ)上建立起來的非靜壓模型[14],將內(nèi)部造波法引入CIP求解對(duì)流項(xiàng)的Navier-Stokes 方程的數(shù)值模型[15]。相關(guān)研究表明,采用相位聚焦方式可以形成畸形波,由于畸形波的強(qiáng)非線性特征,線性理論和二階Stokes波理論均無法精準(zhǔn)預(yù)報(bào)聚焦波面和波峰面下的速度場[1]。

    聚焦波演化形成畸形波的過程中,多種頻率組成成分的波浪非線性相互作用引起的多波共振或調(diào)制不穩(wěn)定等是生成畸形波的最有可能的原因之一。精確預(yù)報(bào)畸形波的生成與發(fā)展過程勢必要求數(shù)值模型具有足夠精準(zhǔn)的色散性和非線性特征。一方面,聚焦波中波浪頻率組成成分越多,在同一水深下,這些波浪具有不同的無因次水深,要求模型在更大水深范圍內(nèi)具備優(yōu)良的線性色散特性,才能保證每一種頻率波浪的計(jì)算結(jié)果與線性解析解具有良好的相位匹配關(guān)系[16];另一方面,不同頻率的波浪之間存在較為明顯的非線性相互作用,這將產(chǎn)生低頻、高頻等其他頻率的差頻與和頻波浪,精準(zhǔn)的非線性特征(二階和差頻、三階非線性波幅、波幅色散)則可以保證數(shù)值模型在模擬非線性波浪傳播時(shí)具有較好的適用性[16-19]。當(dāng)方程在描述速度分布的精度不準(zhǔn)確時(shí),數(shù)值計(jì)算得到的速度場與相關(guān)實(shí)驗(yàn)結(jié)果會(huì)出現(xiàn)一定的差異[13]。劉必勁等[13]采用Liu 和Fang[20]的單層Boussinesq 型水波數(shù)值模型,在開展聚焦波的數(shù)值模擬時(shí)發(fā)現(xiàn),聚焦波面與試驗(yàn)結(jié)果的吻合程度良好,但聚焦波峰值下的水平速度場計(jì)算值呈S 形,與Baldock 等[1]試驗(yàn)的速度分布對(duì)比時(shí)僅基本吻合,這表明受限于單層Boussinesq 水波型方程在速度分布方面的精度,該數(shù)值模型無法勝任深水聚焦波速度場的模擬。為此,本文選用了Liu 和Fang[20]的雙層Boussinesq 型方程,建立相應(yīng)的數(shù)值模型,研究該數(shù)值模型在應(yīng)用于模擬深水聚焦波問題時(shí),波面位移特別是聚焦波峰下水平速度的準(zhǔn)確性。

    1 雙層Boussinesq數(shù)值模型

    1.1 控制方程

    在無旋、無粘的假設(shè)下,Liu 和Fang[20]保留了垂向速度,經(jīng)推導(dǎo)給出了適合緩變地形的雙層Boussinesq 方程,立面二維模型由9 個(gè)方程組成,含有9 個(gè)獨(dú)立待求變量:波面位移η、自由表面處水平和垂向速度uη和wη、靜止水位處的速度u0和w0,以及計(jì)算速度(u*i,w*i)(i=1,2)。方程的具體表達(dá)形式如下:

    在自由面處,滿足精確的運(yùn)動(dòng)學(xué)邊界條件和動(dòng)力學(xué)邊界條件:

    自由面處速度與靜止水位處速度之間的表達(dá)式為

    靜止水位處速度與第一層計(jì)算速度的關(guān)系式為

    在連接水深位置處,滿足速度連續(xù),則有

    在水底處,滿足運(yùn)動(dòng)學(xué)邊界條件:

    式中,下標(biāo)x表示對(duì)空間x求導(dǎo),g是重力加速度。

    方程(1)~(9)構(gòu)成了本文數(shù)值模型對(duì)應(yīng)的雙層Boussinesq水波方程,該方程具有較好的線性與非線性性能。對(duì)波浪沿水深分布的速度進(jìn)行優(yōu)化,Liu等[18]給出了方程中的系數(shù)(α1,α2,β12,β13,β22,β23)=(0.129,0.371,0.941,0.896,0.948,0.608)。

    1.2 數(shù)值模型的求解

    雙層Boussinesq 水波方程的最高空間導(dǎo)數(shù)是3,其與目前常用于海岸波浪模擬的開源程序FUNWAVE[21]對(duì)應(yīng)的Boussinesq 水波[22]方程最高空間導(dǎo)數(shù)的階數(shù)上是一致的。為此,該方程在時(shí)空差分格式的離散上均采用與該開源代碼一樣的數(shù)值格式。在時(shí)間差分格式上,采用混合4 階Adams-Bashforth-Moulton 格式。其中,預(yù)報(bào)時(shí)利用3 階Adams-Bashforth 格式求解方程(1)和(2),可以得到波面位移以及波面處水平速度的預(yù)報(bào)值;分別求解方程(3)、(5)和(7)可得到3個(gè)水平速度的預(yù)報(bào)值;進(jìn)而求解方程(9)、(8)、(6)和(4),可依次得到4 個(gè)垂向速度的預(yù)報(bào)值。校正時(shí)則利用4 階Adams-Moulton 時(shí)間步進(jìn)格式,得到波面位移以及波面處水平速度的校正值,其他速度變量校正值的求解過程類似于預(yù)報(bào)階段,不同點(diǎn)在于所有方程右端項(xiàng)均是當(dāng)前預(yù)報(bào)層的值。當(dāng)9 個(gè)變量的校正值與預(yù)報(bào)值之間的誤差控制在0.000 1 內(nèi),當(dāng)前迭代結(jié)束,否則更新變量的預(yù)報(bào)值,其采用校正值和預(yù)報(bào)值的加權(quán)值,繼續(xù)進(jìn)行當(dāng)前步的迭代過程。校正值的加權(quán)系數(shù)R為0.05~0.1,而預(yù)報(bào)值加權(quán)系數(shù)為1-R。

    2 數(shù)值模型的驗(yàn)證

    本章將重點(diǎn)研究聚焦波的一些特性,并與相關(guān)理論解析解和實(shí)驗(yàn)結(jié)果進(jìn)行比較,驗(yàn)證本文數(shù)值模型的適用性。首先開展線性聚焦波的模擬,考察數(shù)值模型在具有多個(gè)不同周期組成下線性色散的適用性問題,其最重要的表現(xiàn)是疊加后波面的相位問題,當(dāng)多波組成的波面與解析波面一致時(shí),則表明數(shù)值模型在更大無因次水深范圍內(nèi)有較好的色散精度;其次,開展非線性深水聚焦波的模擬,考察數(shù)值模型波-波非線性相互作用和非線性色散的精確性問題。

    2.1 線性聚焦波的模擬

    雙層Boussinesq 方程具有確切的色散關(guān)系表達(dá)式,其與Stokes 波色散關(guān)系表達(dá)式的對(duì)比在Liu 和Fang[20]以及Liu 等[18]中均明確給出,本文的雙層Boussinesq 數(shù)值模型中采用了Liu 等[18]給出的參數(shù)值,原圖可詳見Liu 等[18]中的圖3,這里為了展現(xiàn)其色散性能對(duì)圖進(jìn)行了放大,重新繪制得到的無因次相速度見圖1。由圖1 可見,方程的相速度在無因次水深0<kh<20范圍內(nèi)與解析解的誤差僅為0.003 68%,這表明雙層Boussinesq 水波方程具有極為精確的色散特征。盡管方程具有精確的色散關(guān)系式,但數(shù)值模型因離散誤差以及數(shù)值積累誤差等,一般來說,其計(jì)算出來的色散方面的性能低于理論分析得到的色散性能。

    當(dāng)忽略非線性項(xiàng)后,雙層Boussinesq 模型退化為線性模型,利用該模型對(duì)線性聚焦波進(jìn)行模擬,目的是考察建立的雙層Boussinesq 數(shù)值模型的精準(zhǔn)度,這里選擇Baldock 等[1]給出的波浪頻率進(jìn)行數(shù)值計(jì)算。Baldock等[1]的聚焦波深水實(shí)驗(yàn)含有不同波浪頻率段范圍,在不同周期段內(nèi)等分周期成29份,每個(gè)頻率波幅取值相等,其中的B 組(T=0.6~1.4 s,寬譜)和D 組(T=0.8~1.2 s,窄譜)對(duì)應(yīng)的線性無因次水深kh分別為1.568~7.825 和2.026~4.403。在數(shù)值模擬中,水深為0.7 m,計(jì)算域設(shè)置為24 m,末端設(shè)置4 m長的海綿邊界層吸收波浪,時(shí)間步長為0.01 s,空間步長為0.04 m。設(shè)置線性聚焦位置在8 m,線性聚焦時(shí)間為60 s,聚焦波幅為0.055 m,選擇這一聚焦時(shí)間點(diǎn)的原因是避開數(shù)值模型在模擬波群演化中波前(wave front)的影響。在聚焦時(shí)刻t=60 s,B 組和D 組的計(jì)算結(jié)果與一階Stokes 波解析解的比較見圖2 和圖3。由圖可見,數(shù)值模擬結(jié)果與解析解吻合程度高,這一方面表明數(shù)值模型可以精準(zhǔn)捕捉不同頻率波浪的相位,另一方面也從數(shù)值角度佐證了如圖1展示的雙層Boussinesq水波方程在較大頻率范圍內(nèi)精準(zhǔn)的相速度。此外,在聚焦位置處,兩組工況的計(jì)算波峰面均是0.054 7 m,與解析波峰面0.055 m的誤差是0.5%。

    圖1 無因次相速度Fig.1 Non-dimensional phase velocity

    圖2 B55計(jì)算結(jié)果與Stokes波解析解的比較Fig.2 Comparison of the simulated surface elevation with the analytical solution for Case B55

    圖3 D55計(jì)算結(jié)果與Stokes波解析解的比較Fig.3 Comparison of the simulated surface elevation with the analytical solution for Case D55

    2.2 非線性聚焦波的模擬

    在2.1 節(jié)線性聚焦波的模擬中,不同成分的線性波浪滿足疊加關(guān)系,以上僅是表征方程在線性波群中相位捕捉的能力。實(shí)際中,非線性波浪并不滿足線性疊加原理。一方面,聚焦波組成成分之間存在明確的非線性相互作用,這會(huì)導(dǎo)致波群演化過程中出現(xiàn)多個(gè)頻率之間產(chǎn)生和頻項(xiàng)與差頻項(xiàng),因此,要求方程必須具有足夠精確的和頻與差頻性能;另一方面,根據(jù)三階Stokes波理論可知,波幅越大,波長將越長(同等條件下),在相位聚焦生成大波的過程中,波峰面越大,將會(huì)導(dǎo)致其與設(shè)置的線性聚焦位置發(fā)生更多偏離[1,14]。

    在進(jìn)行數(shù)值模擬之前,首先確認(rèn)方程的非線性和速度分布特征方面的適用范疇。在1%誤差下,雙層Boussinesq 水波方程在二階非線性、三階非線性和波幅色散三個(gè)特性方面的最大適用水深kh是35.7、33.1和36.7,而線性水平速度和垂向速度的最大適用水深kh均是17.6,關(guān)于這些特性對(duì)應(yīng)的詳圖可見參考文獻(xiàn)[18]。這表明雙層Boussinesq 數(shù)值模型可以勝任深水聚焦波演化的模擬。其次,對(duì)所采用的聚焦波物理試驗(yàn)進(jìn)行說明。第一,水槽中開展的波浪試驗(yàn),波高在沿程方向存在一定程度的衰減現(xiàn)象;第二,造波機(jī)驅(qū)動(dòng)信號(hào)產(chǎn)生的波浪與理論解存在差異;第三,實(shí)際水槽無法克服邊界的反射問題,這種反射波會(huì)“污染”試驗(yàn)結(jié)果。

    針對(duì)Baldock 等[1]試驗(yàn)中B 組和D 組的系列工況開展數(shù)值模擬研究,在計(jì)算中,數(shù)值水槽長24 m,末端設(shè)置4 m 的海綿邊界層消波。空間步長和時(shí)間步長分別是0.04 m 和0.01 s,計(jì)算時(shí)間為40 s。對(duì)B 組和D 組分別含有三組線性聚焦波峰面Af(分別是22 mm、38 mm 和55 mm)工況進(jìn)行數(shù)值模擬。在入射邊界,造波采用線性波浪條件,設(shè)置的線性聚焦位置是8 m,通過反算可以確定入射波浪的情況,聚焦時(shí)間設(shè)置為30 s,其他條件均與Baldock 等[1]實(shí)驗(yàn)描述的一致。實(shí)際計(jì)算得到的最大聚焦位置處,波面位移與實(shí)驗(yàn)結(jié)果[1]的對(duì)比見圖4。相應(yīng)地,最大聚焦波峰面下的水平速度與實(shí)驗(yàn)結(jié)果的對(duì)比見圖5。從圖4中可見:非線性較強(qiáng)的工況B55和D55的計(jì)算結(jié)果比實(shí)驗(yàn)結(jié)果大;作為非線性相對(duì)弱一些的D38 工況,計(jì)算波面比實(shí)驗(yàn)結(jié)果小,其他非線性工況計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果相差不大。從圖5 對(duì)比的水平速度場來分析,除D22 外,其他模擬結(jié)果與實(shí)驗(yàn)結(jié)果的吻合程度較好,同時(shí)雙層模型計(jì)算的速度剖面比較光滑,不會(huì)出現(xiàn)單層模型給出的S 形[13],這表明雙層Boussinesq 水波方程具有良好的速度分布特性。

    圖4 計(jì)算波面位移與Baldock等[1](1996)實(shí)驗(yàn)結(jié)果的對(duì)比Fig.4 Comparison of the calculated surface elevation with the experimental results of Baldock et al[1](1996)

    圖5 水平速度與Baldock等[1](1996)實(shí)驗(yàn)結(jié)果的對(duì)比Fig.5 Comparison of the calculated horizontal velocity with the experimental results of Baldock et al[1](1996)

    此外,實(shí)驗(yàn)中流速儀無法捕捉到波峰面處的速度,而數(shù)值模型可以提供波峰面速度,因此,速度場的最大值高于試驗(yàn)值是可以接受的。所有工況算例表明,雙層Boussinesq數(shù)值模型可以較為準(zhǔn)確地模擬深水聚焦波演化和垂向分布的速度場。

    根據(jù)前面提及水槽實(shí)驗(yàn)中的波高衰減等問題,對(duì)兩組非線性最強(qiáng)的兩個(gè)算例(B55和D55)減少入射波高進(jìn)行模擬,其中,B55工況對(duì)應(yīng)的線性聚焦波峰面Af設(shè)置為0.053 m,而D55工況對(duì)應(yīng)的Af設(shè)置為0.053 5 m,具體計(jì)算得到的波面位移和水平速度結(jié)果見圖6。由該圖可見,計(jì)算的波面位移和水平速度均與實(shí)驗(yàn)結(jié)果吻合。

    圖6 計(jì)算的波面位移和水平速度與Baldock等(1996)實(shí)驗(yàn)結(jié)果的對(duì)比Fig.6 Comparison of the calculated surface elevation and horizontal velocity with the experimental results of Baldock et al(1996)

    3 數(shù)值模擬結(jié)果的進(jìn)一步討論與分析

    從已有的實(shí)驗(yàn)研究[1]和相關(guān)數(shù)值模擬研究[6,13]可確知,保持中心頻率不變,在同樣線性聚焦峰面情況下,周期范圍越寬,則非線性越弱。因此本文以D組為例,首先與Baldock等[1]的實(shí)驗(yàn)條件類似,保持線性聚焦位置與聚焦時(shí)間不變化,討論本文建立的雙層Boussinesq數(shù)值模型計(jì)算得到的最大波面位移ηf和最大波面位移處的水平速度uf隨線性聚焦波峰幅值A(chǔ)f的變化情況;其次,以Baldock 等[1]的D55 為例,討論線性聚焦位置(xf=6 m,7 m,8 m,9 m,10 m)對(duì)ηf和uf的影響,同時(shí)為避免短波沒有傳播到聚焦位置等對(duì)數(shù)值結(jié)果的影響,數(shù)值模型中將線性聚焦時(shí)間tf設(shè)置為35 s。最后,仍以Baldock 等[1]的D55為例,討論了計(jì)算域內(nèi)每個(gè)空間點(diǎn)的最大波峰面位移η和對(duì)應(yīng)η處的水平速度uη。

    在數(shù)值模擬中,保持水深為0.7 m,計(jì)算域仍然設(shè)置為24 m,末端設(shè)置4 m 長的海綿吸收邊界,時(shí)間步長為0.01 s,空間步長為0.04 m,計(jì)算時(shí)間設(shè)置為40 s。

    3.1 線性聚焦波峰幅值A(chǔ)f對(duì)ηf和uf的影響

    圖7和圖8分別給出了D組系列工況中的最大波面位移ηf和最大波面位移處的水平速度uf隨線性聚焦波峰幅值A(chǔ)f(0.01~0.055 m,每個(gè)工況間隔0.005 m)的變化情況。由圖7 可見,數(shù)值模型計(jì)算的最大波峰面在不同工況時(shí)與實(shí)驗(yàn)結(jié)果[1]趨勢基本一致且比較接近;進(jìn)一步分析可知,伴隨線性聚焦波峰面(入射波高)的增大,最大波峰面呈現(xiàn)非線性變化。類似地,圖8給出的最大波峰面處的水平速度也具有同樣變化趨勢。

    圖7 最大波面位移隨線性聚焦波峰幅值的變化Fig.7 Variation of the maximum surface elevation with the amplitude of the linearly focused wave crest

    圖8 最大波峰面處水平速度隨線性聚焦波峰幅值的變化Fig.8 Variation of horizontal velocity at the maximum surface elevation with the amplitude of the linearly focused wave crest

    3.2 線性聚焦位置xf對(duì)ηf和uf的影響

    表1 給出了Af=55 mm 工況(tf=35 s)最大波面位移ηf和最大波面位移處的水平速度uf隨線性聚焦位置(xf=6 m,7 m,8 m,9 m,10 m)的變化情況。在表中也給出了相鄰兩次大波峰面之間的傳播速度cf,其計(jì)算公式可寫為

    表1 ηf和uf隨xf的變化情況Tab.1 Variation of ηf and uf with xf

    式中:xf1和xf2分別為相鄰兩次的最大波峰面位置,對(duì)應(yīng)著時(shí)間tf1和tf2;xf1為設(shè)置線性聚焦位置后的第一個(gè)聚焦波峰面位置,xf2為相鄰的最大波峰面位置。需要強(qiáng)調(diào)的是,從數(shù)值計(jì)算結(jié)果來看,xf2可以在xf1之前,也可以在xf1之后,根據(jù)具體的數(shù)值模擬情況確定。

    由表1 可知:一方面,伴隨線性聚焦位置的增大,線性聚焦位置后的第一個(gè)聚焦波峰面位置xf1處的波面和對(duì)應(yīng)的水平速度值均在增大,且變化范圍較大;另一方面,較短的距離范圍內(nèi),不同頻率之間的非線性相互作用還未充分發(fā)展。因此,在線性聚焦位置為xf(5m,6m 和7m)的第1次聚焦之后,又出現(xiàn)了更大的聚焦波峰面,在本文計(jì)算中并未發(fā)現(xiàn)其后有第3 次聚焦大峰。而在其他算例(xf=8m,9m,10 m 和12 m)中,均在最大聚焦位置前,也出現(xiàn)了較大的聚焦波峰面(可以近似認(rèn)為是局部最大波峰)。綜合來看,所有計(jì)算工況全域內(nèi)的最大波峰面在74.42~81.27 mm 之間,這說明線性聚焦位置的設(shè)置對(duì)計(jì)算的最大波峰面產(chǎn)生了較大影響;類似地,所有計(jì)算工況中最大波峰面處的水平速度變化范圍為0.596~0.702 m/s之間,說明線性聚焦位置的設(shè)置對(duì)uf也產(chǎn)生了較大影響。此外,有趣的是式(10)給出的相鄰兩次最大波峰面間的傳播速度cf卻變化很小,其變化范圍在0.79~0.81 m/s之間。所有工況中,波峰面的最大速度為0.702 m/s,其小于最大波峰間的傳播速度0.81 m/s,這間接表明本文計(jì)算工況中波浪均未發(fā)生破碎。

    3.3 波峰面最大值及水平速度的沿程變化

    針對(duì)Baldock等[1]的D55工況,計(jì)算得到的波峰面及水平速度的沿程變化見圖9和圖10。從圖9中可知,以70 mm為界,存在較多的空間點(diǎn)波峰面超過該值,且存在一個(gè)較大的共振區(qū)域。類似地,圖10也表明,波峰面處的水平速度也存在一個(gè)較大區(qū)域內(nèi)的共振現(xiàn)象。

    圖9 波峰面最大值的沿程變化Fig.9 Maximum crest elevations along the flume

    圖10 波峰面最大值處水平速度的沿程變化Fig.10 Horizontal velocity of maximum crest elevations along the flume

    4 結(jié) 論

    本文基于雙層Boussinesq 水波方程建立了聚焦波立面二維數(shù)值模型,將其應(yīng)用于深水聚焦波的數(shù)值模擬中,得出的主要結(jié)論如下:

    (1)數(shù)值模型計(jì)算線性聚焦波面處的波面位移與一階Stokes波解析解吻合程度很高,表明本文建立的雙層Boussinesq數(shù)值模型能夠精確捕捉線性相位問題。

    (2)非線性聚焦演化中,波浪非線性起重要作用,數(shù)值模擬波面與實(shí)驗(yàn)結(jié)果吻合較好,表明雙層Boussinesq數(shù)值模型能勝任聚焦波波面演化的模擬,也印證了數(shù)值模型對(duì)應(yīng)的方程具有良好的非線性特性。此外,對(duì)波峰面下的水平速度場,本文建立的雙層Boussinesq數(shù)值模型也能較好地模擬。

    (3)在強(qiáng)非線性波浪聚焦過程中,將會(huì)出現(xiàn)一個(gè)較大范圍的波浪聚焦區(qū)。以Baldock 等(1996)的D55 為例,在設(shè)置不同線性聚焦位置的情況下,最大波峰面的傳播速度保持在0.79~0.81 m/s 之間,這表明不同線性聚焦位置的設(shè)置對(duì)最大波峰面的傳播速度影響不大,但其對(duì)波峰面下的速度影響較大。

    猜你喜歡
    波面波峰雙層
    作用于直立堤墻與樁柱的波峰高度分析計(jì)算
    基于恒定陡度聚焦波模型的分析與討論
    水道港口(2020年6期)2020-02-22 11:33:50
    墨爾本Fitzroy雙層住宅
    多普勒效應(yīng)中觀察者接收頻率的計(jì)算
    淺談光的干涉和衍射的區(qū)別和聯(lián)系
    中文信息(2018年2期)2018-05-30 11:45:10
    兒童標(biāo)準(zhǔn)12導(dǎo)聯(lián)T波峰末間期的分析
    次級(jí)通道在線辨識(shí)的雙層隔振系統(tǒng)振動(dòng)主動(dòng)控制
    Dynamic Loads and Wake Prediction for Large Wind Turbines Based on Free Wake Method
    波面位移非線性特征數(shù)值研究
    傳統(tǒng)Halbach列和雙層Halbach列的比較
    五月玫瑰六月丁香| 亚洲国产精品国产精品| 免费看光身美女| 91久久精品电影网| 国产成人精品婷婷| 午夜老司机福利剧场| 国产无遮挡羞羞视频在线观看| 久久这里有精品视频免费| 色网站视频免费| 少妇被粗大猛烈的视频| 欧美zozozo另类| a级毛色黄片| 久久婷婷青草| 男女边摸边吃奶| 黑人高潮一二区| 欧美成人a在线观看| 国产成人aa在线观看| 夜夜爽夜夜爽视频| 天堂中文最新版在线下载| 成人国产麻豆网| 欧美日韩视频高清一区二区三区二| 狂野欧美激情性bbbbbb| av视频免费观看在线观看| 久久久久久久久久人人人人人人| 18禁动态无遮挡网站| 99热这里只有是精品在线观看| 51国产日韩欧美| 最近的中文字幕免费完整| a 毛片基地| 精品一区二区三区视频在线| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说 | 最近中文字幕高清免费大全6| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜爱| 在线观看美女被高潮喷水网站| 亚洲婷婷狠狠爱综合网| 午夜精品国产一区二区电影| 国产成人一区二区在线| 街头女战士在线观看网站| 国产视频首页在线观看| 黄色视频在线播放观看不卡| 亚洲精品一二三| 久久久亚洲精品成人影院| 亚洲自偷自拍三级| 精品一区二区三区视频在线| 免费人成在线观看视频色| 国产免费福利视频在线观看| 日韩三级伦理在线观看| 亚洲欧美一区二区三区国产| 中国美白少妇内射xxxbb| 亚洲精品亚洲一区二区| 国产精品精品国产色婷婷| 99热6这里只有精品| h视频一区二区三区| 欧美精品人与动牲交sv欧美| 欧美性感艳星| 热re99久久精品国产66热6| 欧美老熟妇乱子伦牲交| 国产黄频视频在线观看| av又黄又爽大尺度在线免费看| 伦理电影免费视频| 黄色怎么调成土黄色| 久久国产亚洲av麻豆专区| 少妇熟女欧美另类| 成人一区二区视频在线观看| 免费人成在线观看视频色| 亚洲成色77777| 国内少妇人妻偷人精品xxx网站| 午夜福利在线观看免费完整高清在| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 日韩av在线免费看完整版不卡| 国产一区二区三区av在线| 精品久久久久久电影网| 日韩av免费高清视频| 久久人人爽av亚洲精品天堂 | 内射极品少妇av片p| 五月开心婷婷网| 亚洲欧美精品自产自拍| av国产免费在线观看| 一级片'在线观看视频| 好男人视频免费观看在线| 久久韩国三级中文字幕| 两个人的视频大全免费| 99热国产这里只有精品6| 一级a做视频免费观看| 久久综合国产亚洲精品| 一区二区三区免费毛片| 女性被躁到高潮视频| 麻豆乱淫一区二区| 久久精品久久久久久久性| 美女视频免费永久观看网站| 欧美成人午夜免费资源| 精品一区二区三卡| 亚洲无线观看免费| 免费看光身美女| 国产免费福利视频在线观看| 久久久午夜欧美精品| 久久久久精品久久久久真实原创| 亚洲一级一片aⅴ在线观看| 久久97久久精品| 一区二区三区四区激情视频| 日本av免费视频播放| av网站免费在线观看视频| 这个男人来自地球电影免费观看 | 草草在线视频免费看| 日韩三级伦理在线观看| 激情 狠狠 欧美| 久久午夜福利片| 亚洲国产精品成人久久小说| 亚洲精品国产av蜜桃| 精品亚洲乱码少妇综合久久| a级一级毛片免费在线观看| 伦理电影大哥的女人| 国产高潮美女av| www.色视频.com| 春色校园在线视频观看| 亚洲性久久影院| 五月伊人婷婷丁香| 最后的刺客免费高清国语| 夜夜爽夜夜爽视频| 久久久亚洲精品成人影院| 欧美一区二区亚洲| 亚洲精品456在线播放app| 女性生殖器流出的白浆| 亚洲av福利一区| 久久久久精品性色| 国产成人精品福利久久| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| 国产白丝娇喘喷水9色精品| 久久99热这里只有精品18| 久久国内精品自在自线图片| 精品一品国产午夜福利视频| 亚洲美女黄色视频免费看| 免费观看无遮挡的男女| 国产极品天堂在线| 国模一区二区三区四区视频| 人人妻人人看人人澡| 精品人妻视频免费看| 一二三四中文在线观看免费高清| 色视频www国产| 联通29元200g的流量卡| av免费观看日本| 王馨瑶露胸无遮挡在线观看| 国产免费又黄又爽又色| 精品熟女少妇av免费看| 丰满乱子伦码专区| 亚洲人与动物交配视频| 欧美97在线视频| 天天躁夜夜躁狠狠久久av| 久久ye,这里只有精品| 蜜桃亚洲精品一区二区三区| 国产成人精品久久久久久| 狂野欧美激情性xxxx在线观看| 亚洲精品中文字幕在线视频 | 肉色欧美久久久久久久蜜桃| 夫妻午夜视频| 免费看日本二区| 亚洲经典国产精华液单| 麻豆国产97在线/欧美| 欧美xxxx性猛交bbbb| 99久久精品一区二区三区| freevideosex欧美| 国产黄片视频在线免费观看| 免费少妇av软件| 天天躁夜夜躁狠狠久久av| 美女xxoo啪啪120秒动态图| 国产欧美日韩一区二区三区在线 | 亚洲国产日韩一区二区| 高清日韩中文字幕在线| 亚洲欧美日韩另类电影网站 | 夫妻午夜视频| 欧美成人午夜免费资源| 欧美精品一区二区大全| 国产黄片视频在线免费观看| 日韩电影二区| 日本wwww免费看| 97热精品久久久久久| 日韩视频在线欧美| 欧美精品一区二区免费开放| 国产片特级美女逼逼视频| 欧美高清成人免费视频www| 日韩伦理黄色片| 男人和女人高潮做爰伦理| 久久av网站| 久久久成人免费电影| 女性被躁到高潮视频| 国产毛片在线视频| 亚洲国产欧美在线一区| 欧美日韩亚洲高清精品| 日韩av在线免费看完整版不卡| 久久 成人 亚洲| 亚洲第一av免费看| 欧美bdsm另类| a 毛片基地| 黄色配什么色好看| 色哟哟·www| 国产精品蜜桃在线观看| 午夜日本视频在线| 成人无遮挡网站| 国产伦精品一区二区三区四那| 3wmmmm亚洲av在线观看| 最近的中文字幕免费完整| 国产白丝娇喘喷水9色精品| 日韩av免费高清视频| 欧美 日韩 精品 国产| 久久人妻熟女aⅴ| 欧美日韩视频高清一区二区三区二| videos熟女内射| 欧美亚洲 丝袜 人妻 在线| 91精品伊人久久大香线蕉| 精品视频人人做人人爽| 亚洲国产精品专区欧美| 国产精品久久久久久精品古装| 国产精品久久久久久久电影| 国精品久久久久久国模美| 一级爰片在线观看| 亚洲欧洲国产日韩| 亚洲av综合色区一区| 夜夜看夜夜爽夜夜摸| 免费大片18禁| a级一级毛片免费在线观看| 亚洲欧美成人综合另类久久久| 男人爽女人下面视频在线观看| 免费看光身美女| 大片电影免费在线观看免费| 伊人久久国产一区二区| 观看av在线不卡| 久久久色成人| 99精国产麻豆久久婷婷| 高清不卡的av网站| 欧美高清成人免费视频www| 国产成人免费观看mmmm| 午夜精品国产一区二区电影| 只有这里有精品99| 亚洲欧美中文字幕日韩二区| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 精品亚洲成国产av| 国产精品成人在线| 热99国产精品久久久久久7| 伊人久久精品亚洲午夜| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 九九在线视频观看精品| 亚洲国产精品专区欧美| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 午夜免费观看性视频| 老司机影院毛片| 日韩大片免费观看网站| 精品久久久精品久久久| av黄色大香蕉| 熟女电影av网| 国产一区亚洲一区在线观看| 在线免费十八禁| 久久婷婷青草| 婷婷色麻豆天堂久久| 久久久久人妻精品一区果冻| 国产视频内射| 精品少妇黑人巨大在线播放| 3wmmmm亚洲av在线观看| 好男人视频免费观看在线| 久久久a久久爽久久v久久| 男女免费视频国产| 激情 狠狠 欧美| 成人二区视频| 成人亚洲精品一区在线观看 | 久久久久久久大尺度免费视频| 在线观看av片永久免费下载| av不卡在线播放| 人妻制服诱惑在线中文字幕| 激情 狠狠 欧美| tube8黄色片| 亚洲成人中文字幕在线播放| 51国产日韩欧美| 欧美日本视频| 日韩中文字幕视频在线看片 | 2018国产大陆天天弄谢| 丝袜喷水一区| 亚洲av.av天堂| 久热久热在线精品观看| freevideosex欧美| 国产美女午夜福利| 亚洲精品aⅴ在线观看| 欧美成人一区二区免费高清观看| 一个人看的www免费观看视频| 丰满迷人的少妇在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品秋霞免费鲁丝片| 国产视频内射| 只有这里有精品99| 国产精品99久久99久久久不卡 | 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 色5月婷婷丁香| 男女边吃奶边做爰视频| 成人综合一区亚洲| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产鲁丝片午夜精品| 亚洲熟女精品中文字幕| 性色av一级| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 亚洲国产高清在线一区二区三| 91精品伊人久久大香线蕉| 久久精品国产亚洲av涩爱| 国产精品爽爽va在线观看网站| 美女中出高潮动态图| 七月丁香在线播放| 中文字幕免费在线视频6| 自拍偷自拍亚洲精品老妇| 亚洲国产精品成人久久小说| 男女免费视频国产| 久久国内精品自在自线图片| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 日韩国内少妇激情av| 免费av中文字幕在线| 乱码一卡2卡4卡精品| 街头女战士在线观看网站| 亚洲欧美一区二区三区国产| 中文字幕av成人在线电影| 一本一本综合久久| 亚洲av日韩在线播放| 欧美日韩国产mv在线观看视频 | 久久国内精品自在自线图片| 成年美女黄网站色视频大全免费 | 久久精品久久精品一区二区三区| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 一级毛片aaaaaa免费看小| 国产欧美另类精品又又久久亚洲欧美| 国产成人精品一,二区| www.色视频.com| 亚洲图色成人| 精品久久久久久电影网| 新久久久久国产一级毛片| 一区二区av电影网| 老熟女久久久| 亚洲伊人久久精品综合| 久久精品国产鲁丝片午夜精品| 成人国产麻豆网| 亚洲av不卡在线观看| av播播在线观看一区| 狂野欧美激情性bbbbbb| 中文字幕av成人在线电影| 久久久久网色| 两个人的视频大全免费| 国产成人aa在线观看| 国产爱豆传媒在线观看| 久久精品夜色国产| 亚洲精品自拍成人| 国产中年淑女户外野战色| 亚洲精品自拍成人| 成年免费大片在线观看| 日本与韩国留学比较| 亚洲精品成人av观看孕妇| 亚洲精品亚洲一区二区| 99热这里只有是精品在线观看| 蜜臀久久99精品久久宅男| 免费黄色在线免费观看| 国产精品偷伦视频观看了| 少妇熟女欧美另类| 国产午夜精品久久久久久一区二区三区| 一级黄片播放器| 亚洲一区二区三区欧美精品| 日本黄色片子视频| 国产 一区精品| 热re99久久精品国产66热6| 国产精品成人在线| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 国产91av在线免费观看| 99热网站在线观看| 国产成人a区在线观看| 国产一区亚洲一区在线观看| 精品久久久噜噜| 亚洲国产欧美在线一区| 国产伦精品一区二区三区四那| 多毛熟女@视频| 秋霞伦理黄片| 国产亚洲最大av| 日韩一区二区视频免费看| 日韩av免费高清视频| 日韩一区二区视频免费看| 日本猛色少妇xxxxx猛交久久| 一级毛片aaaaaa免费看小| 精品久久国产蜜桃| 久热这里只有精品99| 国产精品欧美亚洲77777| 国产精品国产av在线观看| 伦理电影免费视频| 日韩国内少妇激情av| 国精品久久久久久国模美| 亚洲中文av在线| 91精品伊人久久大香线蕉| 看免费成人av毛片| 毛片一级片免费看久久久久| 精品国产一区二区三区久久久樱花 | 免费看av在线观看网站| 久久精品久久久久久噜噜老黄| 亚洲精品乱码久久久v下载方式| 精品国产三级普通话版| 新久久久久国产一级毛片| 亚州av有码| 久久国产精品男人的天堂亚洲 | 十八禁网站网址无遮挡 | 久久精品国产亚洲av天美| 一区二区av电影网| 亚洲第一av免费看| 狂野欧美激情性xxxx在线观看| 91狼人影院| 啦啦啦中文免费视频观看日本| 伦理电影大哥的女人| 久久午夜福利片| 欧美三级亚洲精品| 狂野欧美白嫩少妇大欣赏| 国产黄频视频在线观看| 蜜桃久久精品国产亚洲av| 国产成人91sexporn| 亚洲av成人精品一区久久| 精品久久久久久久久av| 日日摸夜夜添夜夜添av毛片| 亚洲aⅴ乱码一区二区在线播放| av福利片在线观看| 男人爽女人下面视频在线观看| 亚洲精品乱码久久久久久按摩| 色视频www国产| 美女主播在线视频| av黄色大香蕉| 麻豆国产97在线/欧美| 91精品一卡2卡3卡4卡| 久久青草综合色| 日韩国内少妇激情av| 亚洲av电影在线观看一区二区三区| 日韩电影二区| 日日撸夜夜添| 国产v大片淫在线免费观看| 男的添女的下面高潮视频| 热99国产精品久久久久久7| 大陆偷拍与自拍| 嫩草影院新地址| 高清视频免费观看一区二区| 国产成人免费观看mmmm| 麻豆乱淫一区二区| 国产高潮美女av| 性色avwww在线观看| 看免费成人av毛片| 丰满迷人的少妇在线观看| 日韩一区二区视频免费看| av又黄又爽大尺度在线免费看| 99久久人妻综合| 黄片wwwwww| 国产无遮挡羞羞视频在线观看| 日韩人妻高清精品专区| av不卡在线播放| 九九在线视频观看精品| 亚洲aⅴ乱码一区二区在线播放| 国产成人一区二区在线| 男女国产视频网站| 免费人成在线观看视频色| 18禁裸乳无遮挡动漫免费视频| 蜜桃久久精品国产亚洲av| 精品久久久久久久久亚洲| 国产69精品久久久久777片| av黄色大香蕉| 日日撸夜夜添| 亚洲av中文字字幕乱码综合| 综合色丁香网| 插逼视频在线观看| 欧美精品亚洲一区二区| 亚洲天堂av无毛| 国产日韩欧美亚洲二区| 丰满乱子伦码专区| 五月开心婷婷网| 毛片女人毛片| 汤姆久久久久久久影院中文字幕| 中文字幕免费在线视频6| 免费不卡的大黄色大毛片视频在线观看| a级一级毛片免费在线观看| 免费大片18禁| 人体艺术视频欧美日本| 最近的中文字幕免费完整| 日本午夜av视频| 国产日韩欧美亚洲二区| 亚洲精品国产色婷婷电影| 黄色一级大片看看| 国产老妇伦熟女老妇高清| 国产黄频视频在线观看| 极品少妇高潮喷水抽搐| 亚洲精品国产色婷婷电影| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 亚洲一区二区三区欧美精品| 内地一区二区视频在线| 麻豆国产97在线/欧美| 国产探花极品一区二区| 久久午夜福利片| 人体艺术视频欧美日本| 99热全是精品| 亚洲国产精品专区欧美| 亚洲人成网站在线观看播放| 毛片一级片免费看久久久久| 1000部很黄的大片| 欧美 日韩 精品 国产| 亚洲精品色激情综合| 91狼人影院| 一级毛片我不卡| 国内少妇人妻偷人精品xxx网站| 亚洲av综合色区一区| 丰满人妻一区二区三区视频av| h视频一区二区三区| 国产高清有码在线观看视频| 亚洲av日韩在线播放| 性色av一级| 我的女老师完整版在线观看| 国产又色又爽无遮挡免| 一区二区三区免费毛片| 亚洲av国产av综合av卡| 狠狠精品人妻久久久久久综合| 日韩中文字幕视频在线看片 | 91在线精品国自产拍蜜月| 在线天堂最新版资源| 青春草国产在线视频| 一区二区av电影网| 久久毛片免费看一区二区三区| 日日摸夜夜添夜夜爱| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| 在线观看三级黄色| 国产伦在线观看视频一区| 日本一二三区视频观看| 在线观看免费日韩欧美大片 | 蜜桃亚洲精品一区二区三区| 日韩三级伦理在线观看| 欧美人与善性xxx| 国产亚洲精品久久久com| 亚洲精品,欧美精品| 新久久久久国产一级毛片| 伦理电影免费视频| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 亚洲国产av新网站| 国产av一区二区精品久久 | 高清视频免费观看一区二区| 亚洲精品国产av成人精品| 这个男人来自地球电影免费观看 | 国产v大片淫在线免费观看| 成人国产av品久久久| 嘟嘟电影网在线观看| 亚州av有码| 欧美人与善性xxx| 久久国内精品自在自线图片| av免费观看日本| 日韩视频在线欧美| 国产免费视频播放在线视频| 欧美精品一区二区大全| 亚洲精品乱码久久久v下载方式| 麻豆国产97在线/欧美| 一个人看的www免费观看视频| 亚洲欧美一区二区三区国产| 少妇人妻久久综合中文| 美女中出高潮动态图| 亚洲aⅴ乱码一区二区在线播放| 交换朋友夫妻互换小说| 久久精品国产亚洲av天美| 国产精品欧美亚洲77777| 久久国产精品大桥未久av | 另类亚洲欧美激情| av国产久精品久网站免费入址| 亚洲精品自拍成人| 最近最新中文字幕大全电影3| 成人18禁高潮啪啪吃奶动态图 | 夜夜爽夜夜爽视频| 国精品久久久久久国模美| 午夜福利网站1000一区二区三区| 亚洲人成网站在线播| 欧美另类一区| 全区人妻精品视频| 干丝袜人妻中文字幕| 在线观看三级黄色| 岛国毛片在线播放| 亚洲精品成人av观看孕妇| 七月丁香在线播放| 国产在线视频一区二区| 丝袜脚勾引网站| 亚洲综合色惰| 在线精品无人区一区二区三 | 91aial.com中文字幕在线观看| 国产高清三级在线| 啦啦啦中文免费视频观看日本| 国产爱豆传媒在线观看| 久久精品国产a三级三级三级| 深爱激情五月婷婷| 国产免费福利视频在线观看| 欧美日韩亚洲高清精品| 91久久精品国产一区二区三区| 成人免费观看视频高清| 青春草国产在线视频| 五月玫瑰六月丁香| 国产黄片美女视频| 99热全是精品| 中文精品一卡2卡3卡4更新| 日本黄色日本黄色录像| 少妇 在线观看| 国产精品国产三级国产av玫瑰| 欧美激情极品国产一区二区三区 | 少妇被粗大猛烈的视频| 赤兔流量卡办理| 2021少妇久久久久久久久久久| 亚洲成人中文字幕在线播放| 久久毛片免费看一区二区三区|