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

    旋渦聲散射特性的尺度效應(yīng)數(shù)值研究*

    2021-06-01 08:31:58馬瑞軒王益民張樹海武從海王勛年
    物理學(xué)報(bào) 2021年10期
    關(guān)鍵詞:散射截面旋渦聲壓

    馬瑞軒 王益民 張樹海 武從海 王勛年

    1) (中國空氣動力研究與發(fā)展中心, 空氣動力學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 綿陽 621000)

    2) (中國空氣動力研究與發(fā)展中心, 氣動噪聲控制重點(diǎn)實(shí)驗(yàn)室, 綿陽 621000)

    以聲波為主要表現(xiàn)形式的膨脹過程和以旋渦為主要表現(xiàn)形式的剪切過程之間的非線性耦合問題一直以來都是流體力學(xué)的研究熱點(diǎn).尤其是旋渦對聲波的散射問題, 具有重要的科學(xué)意義與工程應(yīng)用背景.本文通過線性緊致格式直接數(shù)值求解二維歐拉方程, 獲得了平面聲波穿過均熵Taylor渦的散射特性.與之前經(jīng)典文獻(xiàn)中的標(biāo)準(zhǔn)算例比較, 結(jié)果極其吻合, 直接驗(yàn)證了研究所采用的高精度高分辨率空間差分和時間推進(jìn)格式以及遠(yuǎn)場無反射邊界條件(緩沖區(qū))的計(jì)算方法在時域同時解析動力學(xué)量和聲學(xué)量(量級遠(yuǎn)遠(yuǎn)小于動力學(xué)量)的有效性.通過引入散射截面, 將全區(qū)域的散射分為長波近似區(qū)、共振散射區(qū)和幾何聲學(xué)區(qū).針對每個子區(qū)域,重點(diǎn)分析了無量綱尺度量旋渦強(qiáng)度和長度尺度比對散射聲場的影響, 給出了散射聲場關(guān)于上述兩個關(guān)鍵無量綱參數(shù)的尺度律關(guān)系, 并且得到了極低馬赫數(shù)極大波長時散射聲場的分布函數(shù).在此基礎(chǔ)上給出了關(guān)于旋渦聲散射物理機(jī)制的一種解釋.

    1 引 言

    旋渦流動中的聲傳播問題作為氣動聲學(xué)的經(jīng)典問題, 不僅可以直接應(yīng)用于聲目標(biāo)識別與探測,旋渦主動降噪等實(shí)際工程問題[1?5], 而且對于發(fā)展氣動聲學(xué)計(jì)算方法[6]以及認(rèn)識復(fù)雜流動(如剪切層、湍流)與聲波相互作用及其發(fā)聲機(jī)理有著重要意義[7?12], 極具應(yīng)用與科研價(jià)值.

    聲波穿過旋渦流動時會產(chǎn)生強(qiáng)烈的非線性散射, 頻率、幅值和相位會發(fā)生顯著的變化, 而且聲波的傳播也會影響流場本身隨時間發(fā)展的運(yùn)動規(guī)律.旋渦流動對聲波的散射問題有著悠久的研究歷史, 其中聲波與湍流的相互作用問題的研究歷史甚至超過氣動聲學(xué)本身.早在1941年, 前蘇聯(lián)學(xué)者Obukhov[13]就開展了湍流對聲散射的研究.Kraichnan[10]和Lighthill[11]針對湍流聲散射問題采用Lighthill聲比擬方法發(fā)展了最早的理論預(yù)估模型.Howe[14]發(fā)展了用來描述包含多尺度連續(xù)旋渦對聲波散射的動力學(xué)方程.Clifford和Brown[15]基于含源項(xiàng)的亥姆霍茲方程發(fā)展了描述散射聲場強(qiáng)度的理論模型, 并且發(fā)現(xiàn)流場的平均輸運(yùn)速度會帶來散射聲波的多普勒頻移現(xiàn)象.采用拋物型的近似方程, Ostashev等[16]開展了理論研究, Dallois等[17]開展了數(shù)值研究.

    目前絕大多數(shù)針對湍流聲散射的研究主要采用半經(jīng)驗(yàn)的理論預(yù)估模型和數(shù)值計(jì)算方法, 缺乏描述湍流的精確物理模型.而且研究中發(fā)現(xiàn)對聲散射起決定作用的流動結(jié)構(gòu)是一系列具有特定尺度的旋渦結(jié)構(gòu), 因此從單個旋渦出發(fā)研究聲散射具有重要的意義.在早期, Fetter[18]從線性歐拉方程,Oshea[19]采用聲比擬理論都獲得了經(jīng)典點(diǎn)渦對聲波的散射特性.Colonius[20,21]等選取多個基本旋渦流動結(jié)構(gòu), 分別采用直接數(shù)值模擬和理論研究的方法開展了聲散射研究.其中針對Rankine渦, 基于聲比擬理論分析了低馬赫數(shù)下聲散射特性, 采用聲線法獲取了高頻近似解.Ford等[22?24]和Hattori等[25]采用匹配漸近展開的方法研究了球?qū)ΨQ旋渦聲散射的特性, 并且采用直接數(shù)值計(jì)算的方法驗(yàn)證了其理論結(jié)果.Howe[26], Kopiev和Belyaev[27]分別采用渦聲理論和諧波展開的方法研究了Rankine渦的聲散射特性.

    現(xiàn)有的數(shù)學(xué)理論方法直接嚴(yán)格求解旋渦聲散射的控制方程時會遇到一系列技術(shù)上的難題, 發(fā)展理論模型時必須要進(jìn)行數(shù)學(xué)上的近似, 比如漸近展開時要求旋渦流動是低馬赫數(shù), 玻恩近似時要求聲波波長遠(yuǎn)遠(yuǎn)大于渦核半徑, 幾何聲學(xué)理論中的射線追蹤法只適用于高頻聲散射.但是在實(shí)際推導(dǎo)過程中, 大多時候這些近似所需要的數(shù)學(xué)條件無法同時滿足, 而且近似條件也大大地限制了結(jié)果的適用范圍.目前還不存在普遍適用且準(zhǔn)確性高的理論模型, 絕大多數(shù)理論模型缺乏對散射聲場細(xì)節(jié)的解析, 比如相位信息的缺失.近年來, 隨著計(jì)算機(jī)軟硬件的提升, 開展精細(xì)化的數(shù)值模擬成為了一種常用的研究旋渦聲散射的手段.其中, Candel[28]通過數(shù)值求解近似的拋物型方程獲得了聲波穿過Rankine渦的畸變特性, Colonius等[20,21]通過直接數(shù)值求解Navier-Stokes方程獲得了典型緊致渦與非緊致渦的聲散射特性, Clair和Gabara[29]通過數(shù)值求解線性歐拉方程獲得了旋渦運(yùn)動對散射聲場的影響特性.但是這些研究中, 大多是將數(shù)值計(jì)算作為驗(yàn)證理論結(jié)果的手段, 沒有對計(jì)算結(jié)果系統(tǒng)分析進(jìn)而給出散射聲場的定量描述.本文針對均熵Taylor渦, 通過高精度高分辨率方法直接數(shù)值求解二維歐拉方程, 獲得了散射聲場隨時間演化的全部信息.通過分析散射截面與旋渦強(qiáng)度和長度尺度比的尺度律關(guān)系, 將旋渦聲散射進(jìn)行了分類.針對每一類型的散射, 重點(diǎn)分析了散射聲場關(guān)于上述兩個無量綱量的尺度效應(yīng).

    2 物理模型

    2.1 均熵Taylor渦

    首先, 給出本文所用旋渦模型的詳細(xì)數(shù)學(xué)表示.考慮二維圓對稱均熵流動, 在柱坐標(biāo) (r,θ) 下其法向和徑向速度表示如下:

    這就是經(jīng)典Taylor渦 (Taylor, 1918) 在無粘情況下的表述.其中, 速度分量和空間坐標(biāo)分別采用無窮遠(yuǎn)處的聲速c∞和渦核半徑R進(jìn)行無量綱化處理.旋渦的馬赫數(shù)(旋渦強(qiáng)度)定義為

    這里Uv是旋渦的最大速度.再由均熵關(guān)系:

    和聲速關(guān)系:

    可以得到密度和壓強(qiáng)的表達(dá)式:

    其中密度ρ和壓強(qiáng)p分別采用無窮遠(yuǎn)處密度ρ∞和無量綱化;γ表示絕熱指數(shù), 空氣中一般為1.4.

    2.2 平面聲波

    選取無窮遠(yuǎn)處沿x正方向傳播的平面波作為入射波, 具體表達(dá)式如下:

    這里聲學(xué)量的密度ρ′;x方向速度分量u′;y方向速度分量v′和壓強(qiáng)p′采用和 旋 渦動力學(xué)量相同參考量進(jìn)行無量綱化;ω為平面波的角頻率;ε=10?5?1, 以保證聲學(xué)量遠(yuǎn)遠(yuǎn)小于動力學(xué)量.

    3 計(jì)算方法

    與之前的研究類似[25], 我們認(rèn)為聲傳播的過程可以忽略粘性和熱傳導(dǎo)效應(yīng).因此在本文的研究中, 通過數(shù)值求解二維歐拉方程獲取平面聲波穿過旋渦后的散射特性.其守恒形式如下:

    其中U是守恒變量;F和G分別是x,y方向的流通矢量, 具體形式為

    ρ,u,v,p分別是密度、x和y方向的速度、壓強(qiáng);總能量密度

    為了更好地從背景流場(動力學(xué)量)中解析出聲學(xué)量, 空間離散采用經(jīng)典的6階中心緊致格式(邊界處為3階精度), 并且采用4階龍格庫塔格式進(jìn)行時間推進(jìn).為了保持計(jì)算的穩(wěn)定性, 采用8階中心緊致濾波技術(shù)抑制計(jì)算中產(chǎn)生的非物理高頻振蕩.

    計(jì)算在直角坐標(biāo)系(x,y)下進(jìn)行, 矩形計(jì)算域如圖1所示.上節(jié)給出的均熵Taylor渦是非定常歐拉方程的精確解, 在計(jì)算中保持定常狀態(tài), 可作為理想的背景流動.旋渦中心位于計(jì)算域的幾何中心.本文只考慮順時針旋轉(zhuǎn)的旋渦, 逆時針對應(yīng)的結(jié)果只需沿x軸對稱變換即可.左邊界為入口邊界, 按照(8)式給定沿x軸正方向傳播的平面聲波.右邊界為出口邊界, 為了減少邊界反射對計(jì)算的影響, 在右邊界的外設(shè)置緩沖區(qū)(buffer zone)[30,31].在上下邊界處, 我們將物理量分成定常部分和非定常部分, 定常部分按照(1)式, (2)式, (6)式和(7)式給定, 非定常部分采取高階單邊插值的方法處理.

    圖1 計(jì)算域示意圖Fig.1.Schematic diagram of computation configuration.

    計(jì)算中采用等間距的笛卡爾網(wǎng)格, 網(wǎng)格間距滿足能同時解析旋渦動力學(xué)量和聲學(xué)小擾動量的要求, 取 ?x=?y=min{λ/8,R/16}, 其中λ為入射聲波的波長, 且有λ=2πc∞/ω.計(jì)算域設(shè)置為緩沖區(qū)的厚度為 2λ.柯朗數(shù)c∞?t/?x=0.4 , 總計(jì)算時長為20(R+λ)/c∞.

    跟隨文獻(xiàn)[21]中的思想, 散射聲場按如下方式計(jì)算:

    其中psc,p,pvor和pinc分別表示散射聲壓、總壓強(qiáng)、旋渦的動力學(xué)壓強(qiáng)和入射波的聲壓.等到計(jì)算穩(wěn)定后, 開始統(tǒng)計(jì)散射聲壓計(jì)算得到散射聲壓的均方根prms.

    4 方法驗(yàn)證

    在時域直接數(shù)值模擬聲傳播的方法主要有三種: 求解Navier-Stokes方程(文獻(xiàn)[21])、求解線性歐拉方程(文獻(xiàn)[29])、還有本文所采用直接數(shù)值求解歐拉方程[32].

    為了考察本文方法的可行性和計(jì)算的準(zhǔn)確性,選取文獻(xiàn)[21, 29]中的標(biāo)準(zhǔn)算例計(jì)算并進(jìn)行對比.算例中入射聲波波長λ=4R, 均熵Taylor渦的馬赫數(shù)Mv=0.125.為了更好的與文獻(xiàn)結(jié)果對比, 對所得到的散射聲壓均方根prms(r= 8R)使用入射聲壓的幅值pi進(jìn)行歸一化處理, 從而得到無量綱量prms/pi, 其中pi=ε=10?5.下文中如果不作特殊說明, 在指向性分布計(jì)算中, 散射聲壓均方根prms均取r= 8R處的值.

    由圖2可知, 采用本文所建方法得到的結(jié)果與Colonius等在文獻(xiàn)[21]和Clair等在文獻(xiàn)[29]中的計(jì)算結(jié)果吻合度很高, 尤其與文獻(xiàn)[29]中結(jié)果極度吻合, 充分證明了本文所建方法的有效性.相比直接求解Navier-Stokes方程, 歐拉方程不考慮粘性和熱傳導(dǎo)效應(yīng), 保證了均熵Taylor渦在計(jì)算過程中保持定常狀態(tài)(不考慮數(shù)值耗散), 有利于有效地提取出散射聲場, 因此歐拉方程更加適合具有解析解的旋渦對聲波的散射問題.而且相比線性歐拉方程的人為固定背景流動, 歐拉方程的計(jì)算是將動力學(xué)量和聲學(xué)量耦合在一起, 考慮了兩者之間的非線性相互作用, 更加接近真實(shí)物理情況.但是,采用歐拉方程計(jì)算的結(jié)果中同時也包含了聲學(xué)量和動力學(xué)量, 而將這兩部分物理量分離一直以來都是氣動聲學(xué)的難題.因此, 之前大多數(shù)聲散射研究采用將完整的歐拉方程近似為線性歐拉方程的方法達(dá)到分離動力學(xué)量和聲學(xué)量的目的.而本研究采用歐拉方程進(jìn)行計(jì)算, 并將動力學(xué)量與聲學(xué)量成功分離是建立在等熵Taylor渦隨距離指數(shù)衰減的特性.在渦核區(qū)域外, 流體動力學(xué)量快速衰減, 背景流動接近靜止空氣介質(zhì).那么在遠(yuǎn)場存在的擾動成分就只有聲波.

    圖2 驗(yàn)證算例Fig.2.Comparison with previous studies.

    5 計(jì)算結(jié)果與分析

    由第2節(jié)中旋渦和聲波的數(shù)學(xué)表達(dá)式中容易看出: 關(guān)于旋渦對聲波的散射問題中, 有四個重要的物理量—旋渦的渦核半徑R、旋渦最大速度Uv、聲波的波長λ和聲波的傳播速度c∞.前兩個物理量決定了渦的運(yùn)動, 后兩個決了聲傳播的特性.為了衡量兩者的相互作用, 我們互相進(jìn)行無量綱化, 就可以得到兩個重要的無量綱參數(shù): 速度尺度比(這里我們稱之為旋渦強(qiáng)度, 即旋渦馬赫數(shù))Mv和長度尺度比(聲波波長與旋渦渦核半徑之比)λ/R.在本文, 針對旋渦強(qiáng)度Mv從0.015625到0.25和長度尺度比λ/R從0.125到16范圍內(nèi)旋渦聲散射特性進(jìn)行計(jì)算并分析研究.這里旋渦強(qiáng)度和長度尺度比的具體選擇采用二分法, 從需要模擬的最大值開始, 每次縮小1/2, 直至需要模擬的最小值.前期經(jīng)過大量數(shù)值模擬結(jié)果摸索出了這種取法, 可以用盡可能少的算例, 得到比較完整的規(guī)律(包含了聲散射的三個子區(qū)域).

    5.1 散射截面

    為了更加直觀地評價(jià)兩個無量綱尺度參數(shù)Mv和λ/R對旋渦聲散射特性的影響(尺度效應(yīng)), 引進(jìn)散射截面Σ表征旋渦對聲波散射的強(qiáng)弱, 其中:

    可以看出散射截面依賴于積分圓的半徑r.考慮到散射波以球面波的形式向外輻射, 在遠(yuǎn)場(r?R)基本上以衰 減[21,29].因 此 以 πr進(jìn) 行 無 量 綱 化,既使得在遠(yuǎn)場時散射截面保持與半徑r無關(guān), 又可以保證入射波與散射波強(qiáng)度相當(dāng)時, 散射截面的大小與1接近.不失一般性, 計(jì)算散射截面時, 當(dāng)時,r取8R; 當(dāng)λ>πR時,r取16R.

    圖3給出了當(dāng)長度尺度比λ/R為0.125, 0.25,0.5, 1, 2, 4, 8, 16時, 散射截面 Σ 與旋渦強(qiáng)度Mv之間的關(guān)系.對于所計(jì)算范圍的所有λ/R, Σ 都會隨著Mv的增加而增加, 直至達(dá)到1的量級時, 不再會有明顯的變化.一旦散射截面超過100, 其與旋渦強(qiáng)度的平方關(guān)系將不再成立.散射截面與旋渦強(qiáng)度的定量關(guān)系依賴于長度尺度比的大小, 當(dāng)λ/R>1 時,Σ與Mv的平方呈正比; 當(dāng)λ/R<1 時,只有對旋渦強(qiáng)度較小時成立.顯然, 散射波的強(qiáng)弱與旋渦的強(qiáng)度呈正相關(guān), 且最大只能達(dá)到與入射波相同的量級.

    圖3 不同長度尺度比下散射截面與旋渦強(qiáng)度的關(guān)系(對數(shù)坐標(biāo)系)Fig.3.Scattering cross-section Σ potted against M v at different λ /R (Logarithmic coordinate system).

    圖4 給出了當(dāng)旋渦強(qiáng)度Mv為0.015625, 0.3125,0.0625, 0.125, 0.25時, 散射截面Σ與長度尺度比λ/R之間的關(guān)系.對于所計(jì)算范圍的所有Mv, 散射截面Σ都會隨著λ/R的增加而減小, 且減小的速度隨著λ/R的增加而越來越快, 在長度尺度比較大時, 衰減速度的接近 (λ/R)?4.經(jīng)典的理論分析證明在玻恩近似(即長波近似,λ/R→+∞時)下,衰減的速度可達(dá)到最大值 (λ/R)?4.按照文獻(xiàn)[22]針對希爾球渦聲散射的理論分析方法, 對于Taylor渦的聲散射結(jié)果, 在低馬赫數(shù)大波長的理想情況下, 同樣可以得到:

    圖4 不同旋渦強(qiáng)度下散射截面與長度尺度比的關(guān)系(對數(shù)坐標(biāo)系)Fig.4.Scattering cross-section Σ plotted against λ /R at different M v (Logarithmic coordinate system).

    下面考慮旋渦強(qiáng)度Mv和長度尺度比λ/R兩者共同對散射截面 Σ 的影響.將整個旋渦聲散射定義在關(guān)于Mv和λ/R的平面上, 并且將該平面劃分為:長波近似區(qū), 共振散射區(qū)和幾何聲學(xué)區(qū).具體劃分如圖5所示, 以為特征變量.在長波近似區(qū)近似有這是由經(jīng)典的長波近似理論, 即玻恩近似所得到的.在共振散射區(qū)有而λ/R的–4次方律不再成立.在該區(qū)域里入射聲波波長與渦核尺度相當(dāng), 散射聲場強(qiáng)度與入射聲場也相當(dāng), 產(chǎn)生了類似“共振”的散射效應(yīng); 在幾何聲學(xué)區(qū),Mv的2次方律和λ/R的–4次方律都不再成立, 該區(qū)域的聲散射現(xiàn)象可以通過幾何聲學(xué)理論進(jìn)行分析.其中, 長波近似區(qū)和共振散射區(qū)的分界線大致在λ/R=2π , 正好對應(yīng)于亥姆霍茲數(shù)kR=2πR/λ=1 , 其中k為波數(shù).而共振散射區(qū)和幾何聲學(xué)區(qū)通常以S= 1為分界線.下面分別對這三個區(qū)域的散射特性的尺度效應(yīng)做詳細(xì)分析.

    圖5 旋渦強(qiáng)度和長度尺度比共同對散射截面的影響(對數(shù)坐標(biāo)系)Fig.5.Scattering cross-section as a function of M v and λ/R(Logarithmic coordinate system).

    5.2 長波近似區(qū)

    圖6 給出了無量綱時間t= 280時的無量綱化的散射聲壓psc/pi.在長波近似區(qū), 旋渦對聲波的散射很弱, 散射聲場的強(qiáng)度遠(yuǎn)遠(yuǎn)小于入射聲場.散射聲場的指向性類似于四極子聲場, 出現(xiàn)了四道波束(beam), 分別位于θ=45?,?45?,?135?,135?附近, 依次記作第一波束(first beam), 第二波束(second beam), 第三波束(third beam), 第四波束(fourth beam).靠近x正方向的兩道波束強(qiáng)于靠近x軸負(fù)方向的兩道波束, 且任意相鄰的兩道波束之間相位相反.

    圖6 t = 280時的散射聲壓, 其中 M v=0.125 ,λ/R=10Fig.6.Snapshot of scattered pressure withMv=0.125 and λ /R=10 at t = 280.

    圖7 給出了長波近似區(qū)五個典型狀態(tài)的聲壓均方根 (r= 16R) 的指向性特征, 采用無量綱量作為特征量.和瞬時散射聲壓類似, 均方根值也出現(xiàn)了四道波束, 且每道波束的形狀類似一個半正弦函數(shù).散射聲壓均方根prms與旋渦強(qiáng)度Mv呈正比.但是隨著Mv增加, 散射聲場的指向性分布表現(xiàn)出更多的不對稱性, 第一波束的強(qiáng)度逐漸增強(qiáng), 相應(yīng)地, 第二波束減弱.這與文獻(xiàn)[29]中的完全對稱的描述不符.這種差異主要來源于計(jì)算所使用的的控制方程不同.本文所采用歐拉方程模型是比文獻(xiàn)[29]中忽略高階小項(xiàng)的線性歐拉方程更加貼合真實(shí)的物理模型, 而且在旋渦的馬赫數(shù)較高時, 這種非線性的高階項(xiàng)愈發(fā)重要, 因此本文所得到的結(jié)果更為可信.在極低馬赫數(shù)時,指向性分布關(guān)于聲波入射方向完全對稱, 這與基于低馬赫數(shù)大波長的理論分析結(jié)果(玻恩近似[18,19]和匹配漸近展開[22?25])類似.同時, 隨著Mv減小,在θ=0?處的散射逐漸趨于0.散射聲場關(guān)于長度尺度比λ/R的“–2次方律”基本符合.同樣隨著λ/R增加, 不同波束的強(qiáng)度逐步接近, 散射聲場表現(xiàn)了良好的對稱性.

    圖7 長波近似區(qū)的散射聲壓均方根指向性分布Fig.7.Directivity for root-mean-square pressure of the scattered fields in long-wavelength domain.

    綜合以上結(jié)果, 可以得到在長波近似區(qū)散射聲場強(qiáng)度關(guān)于旋渦強(qiáng)度和長度尺度比的尺度律關(guān)系:

    圖8給出了低馬赫數(shù) (Mv=0.015625) 、大波長 (λ/R8) 時的散射聲場指向性分布(r= 16R).從圖8可以看出, 隨著波長的增加, 指向性分布逐步接近四個等強(qiáng)度的波束(參考曲線).我們可以得到在極低馬赫數(shù)、極大波長時散射聲場的分布函數(shù)有如下形式:

    其中AT是一個常數(shù), 針對等熵Taylor渦,AT=36.2038671967512.

    圖8 低馬赫數(shù)大波長下的散射聲壓均方根指向性分布Fig.8.Directivity for root-mean-square pressure of the scattered fields with low Mach number and long wavelength.

    5.3 共振散射區(qū)

    隨著長度尺度比λ/R的減小, 旋渦對聲波的散射逐漸過渡到共振散射區(qū).與長波散射相比, 共振散射區(qū)的散射聲場中, 靠近x軸負(fù)方向的兩道波束持續(xù)衰減, 直至消失.如圖9所示, 采用無量綱量作為特征量, 會發(fā)現(xiàn)靠近x軸正方向的兩道波束其強(qiáng)度仍然滿足λ/R的“–2次方律”,并逐漸向聲波入射方向靠攏.這是長波近似區(qū)向共振散射區(qū)過渡的顯著特征.從圖10可以看出, 當(dāng)λ/R=4時, 散射聲場已經(jīng)完全表現(xiàn)共振散射的特性, 只存在靠近x軸方向的兩道反相的波束.

    圖9 共振散射區(qū)較大波長散射聲場指向性分布Fig.9.Directivity for root-mean-square pressure of the scattered fields in resonance domain with relatively long wavelength.

    圖10 t = 140時的散射聲壓, 其中 M v=0.125 ,λ/R=4Fig.10.Snapshot of scattered pressure withMv=0.125 and λ /R=4 at t = 140.

    圖11 旋渦強(qiáng)度對共振散射區(qū)散射聲場指向性的影響Fig.11.Directivity for root-mean-square pressure of the scattered fields in resonance domain at different vortex strength.

    下面采用無量綱量prms/(piMv) 作為特征量考察散射聲場與旋渦強(qiáng)度的關(guān)系.從圖11可以看出,在共振散射區(qū), 散射聲壓均方根值prms與旋渦強(qiáng)度Mv呈正比.和長波散射類似, 隨著Mv的增加, 散射聲場逐漸表現(xiàn)出不對稱性.結(jié)合以上結(jié)果, 可以得到在共振散射區(qū)散射聲場的強(qiáng)度關(guān)于旋渦強(qiáng)度有如下關(guān)系式:

    針對共振散射區(qū)的短波散射, 如圖12所示,兩道主要的波束逐漸向聲波入射方向靠攏,x軸負(fù)方向的散射徹底消失.從圖13可以看出隨著旋渦強(qiáng)度的增加, 兩道波束發(fā)生干涉, 標(biāo)志著散射由共振散射區(qū)轉(zhuǎn)向幾何聲學(xué)區(qū).

    在共振散射區(qū)域, 散射聲波的強(qiáng)度相比于長波近似區(qū)有明顯的增強(qiáng), 接近入射聲波的強(qiáng)度.這個區(qū)域的聲散射特性對于采用聲波探測旋渦結(jié)構(gòu)、基于旋渦陣列主動降噪等工程實(shí)用的聲學(xué)技術(shù)具有很強(qiáng)的現(xiàn)實(shí)意義.

    圖12 t = 60時的散射聲壓, 其中 M v=0.125 ,λ/R=1Fig.12.Snapshot of scattered pressure withMv=0.125 and λ /R=1 at t = 60.

    圖13 旋渦強(qiáng)度對共振散射區(qū)小波長散射聲場指向性的影響Fig.13.Directivity for root-mean-square pressure of the scattered fields in resonance domain with relatively small wavelength at different vortex strength.

    5.4 幾何聲學(xué)區(qū)

    隨著旋渦強(qiáng)度Mv增加和長度尺度比λ/R的減小, 聲散射逐漸進(jìn)入到幾何聲學(xué)區(qū), 其散射聲場強(qiáng)度高于入射聲場, 具體參考圖14和圖15.由于主波束之間的相互干涉, 產(chǎn)生許多強(qiáng)度較弱的次級波束, 且隨著λ/R減小和Mv的增加, 次級波束的數(shù)量逐漸增多, 散射聲場的指向性更加不規(guī)則.在幾何聲學(xué)區(qū), 散射聲場的強(qiáng)度不再與旋渦強(qiáng)度呈正比.與長波近似區(qū)和共振散射區(qū)不同, 對于低馬赫數(shù)的旋渦, 幾何聲學(xué)區(qū)的散射聲場也表現(xiàn)出了不對稱性.

    5.5 旋渦聲散機(jī)制的初步討論

    圖14 散射聲壓, 其中旋渦強(qiáng)度 M v=0.25 (a) λ /R=1 ; (b) λ /R=0.5 ; (c) λ /R=0.25 ; (d)λ/R=0.125Fig.14.Snapshot of scattered pressure with M v=0.25 : (a) λ /R=1 ; (b) λ /R=0.5 ; (c) λ /R=0.25 ; (d) λ /R=0.125.

    圖15 不同旋渦強(qiáng)度對幾何聲學(xué)區(qū)散射聲場指向性的影響Fig.15.Directivity for root-mean-square pressure of the scattered field in geometrical acoustics domain at different vortex strength.

    旋渦對聲波的散射包含了兩種聲散射的機(jī)制.首先是非線性散射的效應(yīng), 即旋渦與聲波非線性耦在渦核內(nèi).另外一種是長程折射效應(yīng), 主要發(fā)生在渦核外部.表現(xiàn)為背景流動對渦核內(nèi)所產(chǎn)生的散射聲波傳播路徑長距離的偏轉(zhuǎn)作用, 是造成非對稱性的原因, 稱之為旋轉(zhuǎn)模態(tài).

    采用以上分析結(jié)果可以很好地解釋: 為什么隨著波長的減小, 反向散射越來越不明顯.原本由渦核所產(chǎn)生的散射波在前向和反向的強(qiáng)度是一致的,但是與入射波疊加后各個方向出現(xiàn)了差異.而這種差異隨著波長的減小越發(fā)明顯.使得反向的折射效應(yīng)強(qiáng)于正向, 原本反向的散射波便被偏轉(zhuǎn)為接近入射波的方向.合發(fā)聲機(jī)制.由于入射聲波的存在, 旋渦不再保持原來定常狀態(tài), 而是受到聲波的激發(fā)向外輻射出散射聲波以抵消入射聲波的影響.散射聲波的強(qiáng)度與旋渦強(qiáng)度成正比, 此外還與聲波和旋渦的尺度比值密切相關(guān), 當(dāng)兩者接近時, 非線性發(fā)聲達(dá)到最大值即通常意義上的共振.旋渦受激輻射聲波有良好的對稱性, 可以稱之為對稱模態(tài), 且該過程主要發(fā)生

    6 結(jié) 論

    本文通過數(shù)值求解二維歐拉方程, 研究了無量綱尺度參數(shù)(旋渦強(qiáng)度Mv和長度尺度比λ/R)對均熵Taylor渦聲散射特性的影響.隨著旋渦強(qiáng)度的增加與長度尺度比的減小, 旋渦聲散射現(xiàn)象依次經(jīng)歷長波近似區(qū)、共振散射區(qū)以及幾何聲學(xué)區(qū), 并且散射聲場強(qiáng)度逐漸增加, 指向性表現(xiàn)出更多的非對稱性.

    1) 在長波近似區(qū), 散射聲場很微弱, 由四道波束組成, 相鄰兩道波束反相.散射強(qiáng)度與旋渦強(qiáng)度呈正比, 與長度尺度比的平方呈反比; 旋渦強(qiáng)度較弱時, 散射聲場關(guān)于入射方向表現(xiàn)出良好的對稱性;

    2) 在共振散射區(qū), 散射主要發(fā)生在靠近入射波的方向, 包含兩道反相的波束, 其強(qiáng)度與旋渦強(qiáng)度呈正比;

    3) 在幾何聲學(xué)區(qū), 散射聲場強(qiáng)度和入射聲場相當(dāng), 主要集中在旋渦后方區(qū)域, 沒有表現(xiàn)出明確的指向性;

    4) 旋渦對聲波的散射包含了兩種不同的機(jī)制,即非線性的散射效應(yīng)與線性的長程折射效應(yīng).

    需要注意的是本文只考慮了相對靜止的旋渦對聲波的散射.因此聲波穿過旋渦后, 頻率并沒有改變.下一步研究將考慮運(yùn)動旋渦對散射聲場的影響, 尤其是不同形式的運(yùn)動(平動、振動等)給散射聲場帶來的多普勒頻移效應(yīng).

    猜你喜歡
    散射截面旋渦聲壓
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    小心,旋渦來啦
    大班科學(xué)活動:神秘的旋渦
    LHCb =8 TeV的Drell-Yan-Z→e+e-數(shù)據(jù)對部分子分布函數(shù)的影響
    旋渦笑臉
    山間湖
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    基于微波倍頻源太赫茲頻段雷達(dá)散射截面測量
    115In中子非彈性散射截面的實(shí)驗(yàn)測量及蒙特卡羅修正
    核技術(shù)(2016年4期)2016-08-22 09:05:22
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    電測與儀表(2016年9期)2016-04-12 00:30:02
    床上黄色一级片| 99久久精品国产国产毛片| 一级毛片电影观看| av免费观看日本| 欧美区成人在线视频| 精品熟女少妇av免费看| 免费观看在线日韩| av在线蜜桃| 欧美精品一区二区大全| 精品熟女少妇av免费看| 人妻夜夜爽99麻豆av| 日韩欧美 国产精品| 亚洲av中文av极速乱| 亚洲av中文av极速乱| av免费观看日本| 精品一区二区免费观看| 六月丁香七月| 天堂影院成人在线观看| 欧美日韩综合久久久久久| 91av网一区二区| 狂野欧美激情性xxxx在线观看| 亚洲av免费高清在线观看| 亚洲图色成人| 亚洲色图av天堂| 国产一区二区三区av在线| 国产亚洲最大av| 亚洲欧美清纯卡通| 老师上课跳d突然被开到最大视频| 十八禁网站网址无遮挡 | 熟女人妻精品中文字幕| 少妇的逼水好多| 两个人的视频大全免费| 亚洲18禁久久av| 精品久久久精品久久久| 真实男女啪啪啪动态图| 一区二区三区高清视频在线| 精品久久国产蜜桃| 三级毛片av免费| 可以在线观看毛片的网站| 欧美日韩一区二区视频在线观看视频在线 | 久久久久网色| 天堂网av新在线| 国产精品久久久久久久电影| 在线免费观看的www视频| 天天躁日日操中文字幕| 国产久久久一区二区三区| 国产一区二区三区综合在线观看 | 久久精品人妻少妇| 深爱激情五月婷婷| 色视频www国产| 久久亚洲国产成人精品v| 日韩av在线免费看完整版不卡| 精品久久久久久久久av| 91在线精品国自产拍蜜月| 高清日韩中文字幕在线| 看非洲黑人一级黄片| 伊人久久精品亚洲午夜| 看黄色毛片网站| 精品人妻一区二区三区麻豆| 丝瓜视频免费看黄片| 真实男女啪啪啪动态图| 人人妻人人看人人澡| 亚洲成色77777| 91久久精品电影网| av播播在线观看一区| 男女下面进入的视频免费午夜| 色吧在线观看| 午夜免费激情av| 久久久a久久爽久久v久久| 肉色欧美久久久久久久蜜桃 | 久久午夜福利片| 又爽又黄无遮挡网站| 特级一级黄色大片| 亚洲精品成人久久久久久| or卡值多少钱| 建设人人有责人人尽责人人享有的 | 搡老妇女老女人老熟妇| 欧美成人精品欧美一级黄| 国产探花极品一区二区| 99久久中文字幕三级久久日本| 精品久久久精品久久久| 赤兔流量卡办理| 男人舔奶头视频| 舔av片在线| 女人被狂操c到高潮| 欧美+日韩+精品| 中文天堂在线官网| 自拍偷自拍亚洲精品老妇| 天天躁夜夜躁狠狠久久av| 神马国产精品三级电影在线观看| 爱豆传媒免费全集在线观看| 51国产日韩欧美| 80岁老熟妇乱子伦牲交| 国产亚洲5aaaaa淫片| 一级二级三级毛片免费看| 啦啦啦啦在线视频资源| 欧美最新免费一区二区三区| 亚洲国产欧美在线一区| 老司机影院成人| 免费人成在线观看视频色| 男的添女的下面高潮视频| 91久久精品国产一区二区三区| 少妇人妻精品综合一区二区| 不卡视频在线观看欧美| av免费在线看不卡| 最近的中文字幕免费完整| 国产午夜精品一二区理论片| 亚洲成人中文字幕在线播放| 精品久久久精品久久久| 黄色欧美视频在线观看| 只有这里有精品99| 欧美成人a在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲av一区综合| 美女高潮的动态| 国产成人freesex在线| 日韩精品有码人妻一区| 一二三四中文在线观看免费高清| 午夜福利成人在线免费观看| 亚洲经典国产精华液单| 亚洲怡红院男人天堂| 18禁在线无遮挡免费观看视频| 日韩欧美国产在线观看| 肉色欧美久久久久久久蜜桃 | 国产精品美女特级片免费视频播放器| 18+在线观看网站| 街头女战士在线观看网站| 校园人妻丝袜中文字幕| 中文字幕免费在线视频6| av福利片在线观看| 黑人高潮一二区| 在线免费十八禁| 免费av观看视频| 亚洲精品日韩av片在线观看| 在线观看免费高清a一片| 中文字幕人妻熟人妻熟丝袜美| 日本-黄色视频高清免费观看| 国产国拍精品亚洲av在线观看| 国产成人aa在线观看| 91狼人影院| eeuss影院久久| 午夜福利视频精品| 老女人水多毛片| 午夜福利高清视频| 免费大片18禁| 精品久久久久久成人av| 你懂的网址亚洲精品在线观看| 亚洲av福利一区| 最近的中文字幕免费完整| 搞女人的毛片| 亚洲欧洲国产日韩| 大香蕉97超碰在线| 免费看不卡的av| 哪个播放器可以免费观看大片| 亚洲精品影视一区二区三区av| av在线老鸭窝| 国产乱人视频| 一夜夜www| 国产在线一区二区三区精| 成年女人看的毛片在线观看| 一级毛片电影观看| 欧美日韩亚洲高清精品| 国产精品久久久久久久电影| 色综合站精品国产| 男人狂女人下面高潮的视频| 国产精品.久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产精品一区二区在线观看99 | 极品少妇高潮喷水抽搐| 成人毛片a级毛片在线播放| a级毛片免费高清观看在线播放| 欧美丝袜亚洲另类| 亚洲,欧美,日韩| 亚洲欧洲国产日韩| 在线 av 中文字幕| kizo精华| 又大又黄又爽视频免费| 高清欧美精品videossex| kizo精华| 婷婷六月久久综合丁香| 国内精品宾馆在线| 中文资源天堂在线| 在线观看一区二区三区| 午夜福利在线观看吧| 99久久精品一区二区三区| 亚洲自拍偷在线| 午夜老司机福利剧场| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 最后的刺客免费高清国语| 青春草视频在线免费观看| 国产黄频视频在线观看| 好男人视频免费观看在线| 欧美潮喷喷水| 3wmmmm亚洲av在线观看| 午夜免费激情av| 在线播放无遮挡| 搡女人真爽免费视频火全软件| 男人和女人高潮做爰伦理| 亚洲精品一二三| 亚洲国产av新网站| 久久久久精品性色| 日韩国内少妇激情av| 亚洲自偷自拍三级| 麻豆av噜噜一区二区三区| 美女脱内裤让男人舔精品视频| 国产片特级美女逼逼视频| 欧美97在线视频| 哪个播放器可以免费观看大片| 国内精品宾馆在线| 欧美人与善性xxx| 久久久久久伊人网av| 免费黄色在线免费观看| 国产精品一区二区三区四区免费观看| 一级片'在线观看视频| 如何舔出高潮| 国产女主播在线喷水免费视频网站 | 欧美潮喷喷水| 日韩中字成人| 大香蕉97超碰在线| 国产成人freesex在线| 美女国产视频在线观看| 午夜激情福利司机影院| 免费无遮挡裸体视频| 国产黄a三级三级三级人| 韩国高清视频一区二区三区| 国产av码专区亚洲av| 搡老乐熟女国产| 国产在线一区二区三区精| 日韩精品青青久久久久久| 日本av手机在线免费观看| 天堂网av新在线| 五月伊人婷婷丁香| av女优亚洲男人天堂| 精品酒店卫生间| 国产精品三级大全| 欧美成人午夜免费资源| 91午夜精品亚洲一区二区三区| 男女视频在线观看网站免费| 久久久精品免费免费高清| 亚洲怡红院男人天堂| 久久久久性生活片| 精品国产露脸久久av麻豆 | 精品一区二区三卡| 色综合色国产| 七月丁香在线播放| 天天躁日日操中文字幕| 亚洲成人久久爱视频| 欧美高清成人免费视频www| 又大又黄又爽视频免费| 大香蕉97超碰在线| 国产免费又黄又爽又色| 三级国产精品欧美在线观看| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 九九在线视频观看精品| 国产精品日韩av在线免费观看| 一级av片app| 成年女人看的毛片在线观看| 天堂√8在线中文| 精品人妻熟女av久视频| 免费观看精品视频网站| 国产爱豆传媒在线观看| 欧美精品国产亚洲| 国产人妻一区二区三区在| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看| 国产一级毛片在线| 日本-黄色视频高清免费观看| 综合色av麻豆| 18禁动态无遮挡网站| 肉色欧美久久久久久久蜜桃 | 日韩伦理黄色片| 色播亚洲综合网| 秋霞伦理黄片| 国产午夜精品久久久久久一区二区三区| 精品久久久久久久末码| 赤兔流量卡办理| 日韩精品有码人妻一区| 插阴视频在线观看视频| 午夜视频国产福利| 亚洲精品,欧美精品| 一区二区三区乱码不卡18| 国产综合懂色| 成人毛片60女人毛片免费| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 国产乱人视频| 青春草亚洲视频在线观看| 亚洲精品一区蜜桃| 国产精品久久久久久久电影| av又黄又爽大尺度在线免费看| 午夜激情欧美在线| 十八禁网站网址无遮挡 | 听说在线观看完整版免费高清| 成年免费大片在线观看| 日日摸夜夜添夜夜添av毛片| 久久这里只有精品中国| 国产淫片久久久久久久久| 午夜视频国产福利| 中文欧美无线码| 啦啦啦中文免费视频观看日本| 欧美xxxx性猛交bbbb| 内射极品少妇av片p| 成人国产麻豆网| 免费观看无遮挡的男女| 中文字幕亚洲精品专区| 人体艺术视频欧美日本| 午夜日本视频在线| 男女下面进入的视频免费午夜| 日韩强制内射视频| 好男人视频免费观看在线| 久久久久网色| 精品99又大又爽又粗少妇毛片| 赤兔流量卡办理| 亚洲电影在线观看av| 亚洲精品日韩av片在线观看| 欧美精品国产亚洲| 午夜福利视频精品| 国产乱来视频区| 大又大粗又爽又黄少妇毛片口| 岛国毛片在线播放| 亚洲激情五月婷婷啪啪| 亚洲伊人久久精品综合| 一个人免费在线观看电影| 国产精品久久久久久精品电影小说 | ponron亚洲| 99久国产av精品国产电影| 精品不卡国产一区二区三区| 一个人观看的视频www高清免费观看| 特级一级黄色大片| 一级二级三级毛片免费看| 小蜜桃在线观看免费完整版高清| 高清av免费在线| 免费少妇av软件| 久久久久久久久久黄片| 日日啪夜夜爽| 美女大奶头视频| 好男人视频免费观看在线| 蜜桃久久精品国产亚洲av| 成年人午夜在线观看视频 | 18禁在线播放成人免费| av国产久精品久网站免费入址| 九九爱精品视频在线观看| 亚洲精品自拍成人| 97热精品久久久久久| 国产又色又爽无遮挡免| 亚洲乱码一区二区免费版| av在线播放精品| 国产成人午夜福利电影在线观看| 久久久久精品久久久久真实原创| 日韩欧美精品免费久久| 欧美变态另类bdsm刘玥| 婷婷六月久久综合丁香| 婷婷色av中文字幕| 日韩欧美 国产精品| 日本爱情动作片www.在线观看| av在线老鸭窝| 精品久久久久久久久亚洲| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线 | 亚洲丝袜综合中文字幕| 久久国产乱子免费精品| 日韩欧美精品免费久久| 搞女人的毛片| a级毛片免费高清观看在线播放| 最近最新中文字幕大全电影3| eeuss影院久久| 国产精品美女特级片免费视频播放器| 国产精品1区2区在线观看.| 美女xxoo啪啪120秒动态图| 一级爰片在线观看| 少妇被粗大猛烈的视频| 国产69精品久久久久777片| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 亚洲av中文字字幕乱码综合| 国产 亚洲一区二区三区 | 麻豆乱淫一区二区| 在线天堂最新版资源| 国产精品女同一区二区软件| 99热这里只有是精品50| 国产精品无大码| 成人av在线播放网站| 成年版毛片免费区| 成人特级av手机在线观看| 亚洲在线观看片| 日韩精品有码人妻一区| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站 | 国产v大片淫在线免费观看| 国产爱豆传媒在线观看| 亚洲综合精品二区| 可以在线观看毛片的网站| 亚洲av男天堂| 国产v大片淫在线免费观看| 国产亚洲5aaaaa淫片| 性色avwww在线观看| 高清视频免费观看一区二区 | 欧美成人午夜免费资源| 国产精品一区二区三区四区免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲一区高清亚洲精品| 一级毛片黄色毛片免费观看视频| 国精品久久久久久国模美| 少妇的逼水好多| 国产不卡一卡二| 亚洲人与动物交配视频| 大陆偷拍与自拍| 一区二区三区四区激情视频| 精品久久久久久久人妻蜜臀av| 高清日韩中文字幕在线| 18禁动态无遮挡网站| 男女啪啪激烈高潮av片| 一级毛片我不卡| 日韩,欧美,国产一区二区三区| 九色成人免费人妻av| 亚洲精品一区蜜桃| 亚洲熟女精品中文字幕| 51国产日韩欧美| 国产女主播在线喷水免费视频网站 | 老司机影院毛片| 禁无遮挡网站| 午夜福利网站1000一区二区三区| 最近的中文字幕免费完整| 只有这里有精品99| 女的被弄到高潮叫床怎么办| 99久久精品热视频| 国产成人精品一,二区| 亚洲国产欧美人成| 91狼人影院| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花 | 日韩精品青青久久久久久| 麻豆av噜噜一区二区三区| 黄片wwwwww| 亚洲国产精品国产精品| 亚洲18禁久久av| 精品久久久久久久人妻蜜臀av| 免费大片18禁| 成人亚洲精品一区在线观看 | xxx大片免费视频| 午夜视频国产福利| 一夜夜www| av福利片在线观看| 久久久久网色| 人人妻人人澡欧美一区二区| 国产精品国产三级国产专区5o| 男女边摸边吃奶| 亚洲国产色片| 97超碰精品成人国产| 免费观看在线日韩| 亚洲高清免费不卡视频| 色网站视频免费| 亚洲va在线va天堂va国产| 亚洲最大成人中文| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 综合色丁香网| 免费av观看视频| 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添av毛片| 国产女主播在线喷水免费视频网站 | 国产综合精华液| 精品亚洲乱码少妇综合久久| 日日啪夜夜撸| 亚洲欧美清纯卡通| 国产白丝娇喘喷水9色精品| 国产精品女同一区二区软件| 久久久久久久大尺度免费视频| a级一级毛片免费在线观看| 黄色欧美视频在线观看| 亚洲精品久久久久久婷婷小说| 人妻少妇偷人精品九色| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 国产有黄有色有爽视频| 一本久久精品| 岛国毛片在线播放| 亚洲av一区综合| 国产成人freesex在线| 91在线精品国自产拍蜜月| 看免费成人av毛片| 69av精品久久久久久| 看免费成人av毛片| 久久6这里有精品| 一个人看视频在线观看www免费| 嫩草影院新地址| 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 欧美成人午夜免费资源| 国产熟女欧美一区二区| 一级黄片播放器| 又大又黄又爽视频免费| 啦啦啦韩国在线观看视频| 欧美性感艳星| 国产熟女欧美一区二区| 人妻系列 视频| 国产三级在线视频| 高清在线视频一区二区三区| av一本久久久久| www.av在线官网国产| 2021少妇久久久久久久久久久| 国产极品天堂在线| 精品亚洲乱码少妇综合久久| a级毛片免费高清观看在线播放| 岛国毛片在线播放| 亚洲国产色片| 日韩中字成人| 免费少妇av软件| 国产综合精华液| 熟女电影av网| 韩国高清视频一区二区三区| 我的女老师完整版在线观看| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 国内精品美女久久久久久| 青春草国产在线视频| www.色视频.com| 国产激情偷乱视频一区二区| 两个人视频免费观看高清| av.在线天堂| 久久99热6这里只有精品| 久久久久久久久大av| 国产成人午夜福利电影在线观看| 好男人视频免费观看在线| 久久精品人妻少妇| 直男gayav资源| 亚洲精品国产成人久久av| 亚洲av二区三区四区| 最近2019中文字幕mv第一页| 热99在线观看视频| 看黄色毛片网站| 卡戴珊不雅视频在线播放| 国内少妇人妻偷人精品xxx网站| 国产精品久久久久久久久免| 国产av不卡久久| 亚洲最大成人中文| 亚洲国产精品成人久久小说| 男人舔女人下体高潮全视频| 99久久精品一区二区三区| 卡戴珊不雅视频在线播放| 久久久久久九九精品二区国产| 欧美极品一区二区三区四区| 国产亚洲91精品色在线| 欧美激情久久久久久爽电影| 熟女人妻精品中文字幕| 精品国产三级普通话版| 少妇的逼水好多| 午夜福利在线观看吧| 汤姆久久久久久久影院中文字幕 | 在线播放无遮挡| 亚洲综合精品二区| 日韩欧美精品v在线| 99热这里只有是精品50| 免费观看性生交大片5| 日韩欧美 国产精品| 99热这里只有精品一区| 国产亚洲午夜精品一区二区久久 | www.色视频.com| 亚洲国产日韩欧美精品在线观看| 最后的刺客免费高清国语| 久久99蜜桃精品久久| 国产成人精品一,二区| 免费观看性生交大片5| 成人毛片a级毛片在线播放| 九九爱精品视频在线观看| 91久久精品国产一区二区三区| 最近中文字幕2019免费版| 色视频www国产| 国产男人的电影天堂91| 久久鲁丝午夜福利片| av女优亚洲男人天堂| 搡女人真爽免费视频火全软件| kizo精华| 日本欧美国产在线视频| 少妇丰满av| 欧美高清成人免费视频www| 春色校园在线视频观看| 婷婷色综合大香蕉| 亚洲av中文字字幕乱码综合| 一区二区三区乱码不卡18| 国产 一区精品| 日韩一区二区视频免费看| 嫩草影院入口| 久久99蜜桃精品久久| 联通29元200g的流量卡| 国产精品久久视频播放| 日韩欧美精品v在线| 成人欧美大片| 亚洲婷婷狠狠爱综合网| 中文欧美无线码| 久久精品熟女亚洲av麻豆精品 | 一个人看视频在线观看www免费| 在线免费十八禁| 国产精品久久久久久精品电影| 欧美+日韩+精品| 成人性生交大片免费视频hd| 免费av毛片视频| 不卡视频在线观看欧美| 狂野欧美白嫩少妇大欣赏| 在线免费观看的www视频| 九草在线视频观看| 亚洲人与动物交配视频| 亚洲精品乱久久久久久| 色综合色国产| 建设人人有责人人尽责人人享有的 | 久久久精品免费免费高清| 亚洲va在线va天堂va国产| 欧美日韩精品成人综合77777| 久久久久久久久中文| 午夜精品在线福利| 亚洲av成人精品一二三区| 欧美日本视频| 亚洲国产精品成人久久小说|