• <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
    俄罗斯特黄特色一大片| 亚洲色图 男人天堂 中文字幕| 亚洲免费av在线视频| av天堂久久9| 男男h啪啪无遮挡| 欧美亚洲日本最大视频资源| 日韩精品免费视频一区二区三区| 成在线人永久免费视频| 男人操女人黄网站| 国产亚洲午夜精品一区二区久久| 男男h啪啪无遮挡| 成人av一区二区三区在线看 | h视频一区二区三区| 叶爱在线成人免费视频播放| 日韩中文字幕视频在线看片| 国产免费一区二区三区四区乱码| 日韩制服丝袜自拍偷拍| 午夜福利免费观看在线| 日本一区二区免费在线视频| 免费一级毛片在线播放高清视频 | 精品少妇黑人巨大在线播放| 久久香蕉激情| netflix在线观看网站| 亚洲全国av大片| 99国产精品免费福利视频| 这个男人来自地球电影免费观看| 爱豆传媒免费全集在线观看| 欧美 亚洲 国产 日韩一| 97在线人人人人妻| 免费女性裸体啪啪无遮挡网站| 日本wwww免费看| 丁香六月欧美| 自线自在国产av| 女人久久www免费人成看片| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久精品精品| 亚洲 国产 在线| 国产欧美亚洲国产| 日本wwww免费看| 国产成人a∨麻豆精品| 久久久久久人人人人人| 老司机影院毛片| 69精品国产乱码久久久| 国产av国产精品国产| 黄色 视频免费看| 亚洲成av片中文字幕在线观看| 免费人妻精品一区二区三区视频| 窝窝影院91人妻| 每晚都被弄得嗷嗷叫到高潮| 十八禁网站网址无遮挡| 少妇人妻久久综合中文| 欧美精品啪啪一区二区三区 | 久久这里只有精品19| 丁香六月欧美| 一个人免费看片子| 一级片免费观看大全| svipshipincom国产片| 精品福利观看| tube8黄色片| 国产亚洲精品久久久久5区| 一本—道久久a久久精品蜜桃钙片| 久久精品国产综合久久久| 久久久久久久大尺度免费视频| 捣出白浆h1v1| av片东京热男人的天堂| 亚洲一区中文字幕在线| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 久久国产亚洲av麻豆专区| 三上悠亚av全集在线观看| 国产欧美亚洲国产| 国产成人欧美在线观看 | 亚洲国产欧美网| 亚洲伊人色综图| 亚洲性夜色夜夜综合| 亚洲免费av在线视频| 欧美人与性动交α欧美软件| 久久毛片免费看一区二区三区| 这个男人来自地球电影免费观看| 中文字幕av电影在线播放| 18禁观看日本| 成年美女黄网站色视频大全免费| 成人影院久久| 天天躁夜夜躁狠狠躁躁| 精品久久久久久久毛片微露脸 | 12—13女人毛片做爰片一| 啦啦啦 在线观看视频| 99久久国产精品久久久| 欧美精品啪啪一区二区三区 | 老司机午夜福利在线观看视频 | 老汉色∧v一级毛片| 欧美黑人欧美精品刺激| 三级毛片av免费| 黑人操中国人逼视频| 国产日韩欧美亚洲二区| 最近最新免费中文字幕在线| 免费在线观看日本一区| 欧美精品高潮呻吟av久久| 国产免费福利视频在线观看| 国产亚洲欧美精品永久| 欧美日韩黄片免| 19禁男女啪啪无遮挡网站| 久久精品久久久久久噜噜老黄| 少妇猛男粗大的猛烈进出视频| 黄色 视频免费看| 亚洲第一青青草原| 亚洲欧美日韩另类电影网站| 咕卡用的链子| 美国免费a级毛片| 午夜91福利影院| 国产精品99久久99久久久不卡| 飞空精品影院首页| 国产日韩欧美在线精品| 国产伦理片在线播放av一区| 国产福利在线免费观看视频| 亚洲欧美成人综合另类久久久| 欧美大码av| 一本久久精品| 啦啦啦啦在线视频资源| 精品国产一区二区三区四区第35| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品 国内视频| 18禁裸乳无遮挡动漫免费视频| 下体分泌物呈黄色| 国产精品九九99| 色婷婷av一区二区三区视频| 午夜激情av网站| 亚洲,欧美精品.| 久9热在线精品视频| 亚洲成人手机| 亚洲第一av免费看| 国产成人影院久久av| 人人澡人人妻人| 纯流量卡能插随身wifi吗| 天天添夜夜摸| 亚洲精品国产色婷婷电影| 亚洲中文字幕日韩| 最近中文字幕2019免费版| 人妻久久中文字幕网| 久久人人爽人人片av| 老汉色av国产亚洲站长工具| 国产男人的电影天堂91| 蜜桃国产av成人99| 国产福利在线免费观看视频| 精品第一国产精品| 无遮挡黄片免费观看| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 精品亚洲成国产av| 国产黄频视频在线观看| 9191精品国产免费久久| 亚洲精品一二三| 午夜精品久久久久久毛片777| 热re99久久国产66热| 日韩 欧美 亚洲 中文字幕| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 午夜影院在线不卡| 成人黄色视频免费在线看| 天天躁日日躁夜夜躁夜夜| 久热这里只有精品99| 桃红色精品国产亚洲av| 精品国产国语对白av| avwww免费| 超色免费av| 天堂8中文在线网| 日韩制服丝袜自拍偷拍| 电影成人av| 精品国内亚洲2022精品成人 | 久久久久网色| 大香蕉久久网| 欧美少妇被猛烈插入视频| 久久久久久人人人人人| 欧美精品人与动牲交sv欧美| 午夜日韩欧美国产| h视频一区二区三区| 国产熟女午夜一区二区三区| 水蜜桃什么品种好| 久久狼人影院| 性少妇av在线| 另类精品久久| 亚洲黑人精品在线| 亚洲中文日韩欧美视频| 18禁国产床啪视频网站| 在线av久久热| 性色av一级| 亚洲人成77777在线视频| 欧美另类亚洲清纯唯美| 亚洲av成人不卡在线观看播放网 | 国产一级毛片在线| 啪啪无遮挡十八禁网站| 日日摸夜夜添夜夜添小说| 亚洲欧美色中文字幕在线| 久久性视频一级片| 色老头精品视频在线观看| 亚洲精品美女久久久久99蜜臀| 在线av久久热| 午夜精品国产一区二区电影| 日韩大码丰满熟妇| 不卡av一区二区三区| 99re6热这里在线精品视频| 亚洲免费av在线视频| 亚洲三区欧美一区| 日韩欧美免费精品| av国产精品久久久久影院| 亚洲精品中文字幕在线视频| 亚洲av欧美aⅴ国产| 啦啦啦中文免费视频观看日本| tube8黄色片| 欧美精品一区二区免费开放| 老熟妇乱子伦视频在线观看 | 久久久久久久久免费视频了| 美女午夜性视频免费| 一区在线观看完整版| 免费在线观看视频国产中文字幕亚洲 | 在线av久久热| 不卡一级毛片| 久久久久国内视频| 久久久久网色| 美女脱内裤让男人舔精品视频| 久久影院123| av在线老鸭窝| 欧美人与性动交α欧美精品济南到| av不卡在线播放| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 99国产极品粉嫩在线观看| 国产成人影院久久av| 国产在视频线精品| 欧美激情极品国产一区二区三区| 亚洲精品国产色婷婷电影| 日本五十路高清| 亚洲av国产av综合av卡| 女性被躁到高潮视频| 国产亚洲欧美在线一区二区| 黄片播放在线免费| 考比视频在线观看| 精品欧美一区二区三区在线| 亚洲少妇的诱惑av| 在线观看免费午夜福利视频| 亚洲精品在线美女| 亚洲av国产av综合av卡| 9热在线视频观看99| 女人久久www免费人成看片| 精品一区在线观看国产| avwww免费| 久久久国产一区二区| 国产无遮挡羞羞视频在线观看| 狠狠狠狠99中文字幕| 精品第一国产精品| 日韩欧美一区二区三区在线观看 | 亚洲性夜色夜夜综合| 亚洲情色 制服丝袜| 国产日韩欧美在线精品| 视频区图区小说| 黑人巨大精品欧美一区二区mp4| 91av网站免费观看| 日日夜夜操网爽| 日韩电影二区| 大码成人一级视频| 男人舔女人的私密视频| 一区在线观看完整版| 日韩 亚洲 欧美在线| 国内毛片毛片毛片毛片毛片| 精品高清国产在线一区| 国产在线免费精品| 久久青草综合色| 美女中出高潮动态图| 日本wwww免费看| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久人妻精品电影 | 王馨瑶露胸无遮挡在线观看| 亚洲精品粉嫩美女一区| 亚洲情色 制服丝袜| 国产精品香港三级国产av潘金莲| 精品欧美一区二区三区在线| 乱人伦中国视频| 国产欧美日韩一区二区精品| 亚洲欧美激情在线| 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 叶爱在线成人免费视频播放| 黑人操中国人逼视频| 在线天堂中文资源库| 久久久精品94久久精品| 午夜免费观看性视频| 久久久久久亚洲精品国产蜜桃av| 欧美+亚洲+日韩+国产| 久久性视频一级片| 国产精品一二三区在线看| 永久免费av网站大全| 久久精品国产综合久久久| 成年人免费黄色播放视频| 伦理电影免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 两人在一起打扑克的视频| 纯流量卡能插随身wifi吗| 大香蕉久久成人网| 日韩一区二区三区影片| 最新的欧美精品一区二区| 亚洲午夜精品一区,二区,三区| 丝袜美腿诱惑在线| 黑人操中国人逼视频| 在线天堂中文资源库| 丰满迷人的少妇在线观看| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 欧美 日韩 精品 国产| a在线观看视频网站| 黄色视频在线播放观看不卡| 电影成人av| 色婷婷久久久亚洲欧美| 亚洲一码二码三码区别大吗| 咕卡用的链子| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 国产av精品麻豆| 热re99久久国产66热| 成年人免费黄色播放视频| 国产一区二区 视频在线| 国产免费av片在线观看野外av| 欧美日韩亚洲国产一区二区在线观看 | www.精华液| 国产精品麻豆人妻色哟哟久久| 精品一区二区三区四区五区乱码| 亚洲精品一区蜜桃| 桃花免费在线播放| 亚洲成人免费av在线播放| 午夜老司机福利片| 久久精品国产a三级三级三级| 日韩有码中文字幕| 欧美+亚洲+日韩+国产| av福利片在线| 亚洲精品国产av蜜桃| 成人黄色视频免费在线看| 97精品久久久久久久久久精品| 99国产精品一区二区蜜桃av | 咕卡用的链子| 日韩欧美免费精品| 亚洲va日本ⅴa欧美va伊人久久 | 91九色精品人成在线观看| 国产免费现黄频在线看| 精品国产国语对白av| 午夜福利免费观看在线| 欧美日韩国产mv在线观看视频| 久久久久精品人妻al黑| 日本猛色少妇xxxxx猛交久久| 97在线人人人人妻| 国产色视频综合| 99国产精品99久久久久| 欧美人与性动交α欧美软件| 人妻一区二区av| 亚洲免费av在线视频| 成年av动漫网址| 午夜两性在线视频| 黄色怎么调成土黄色| 一区二区日韩欧美中文字幕| 在线av久久热| 真人做人爱边吃奶动态| 午夜精品国产一区二区电影| 免费看十八禁软件| 欧美日韩亚洲国产一区二区在线观看 | 美女午夜性视频免费| 国产极品粉嫩免费观看在线| netflix在线观看网站| 亚洲,欧美精品.| 精品久久蜜臀av无| 三上悠亚av全集在线观看| 久久免费观看电影| 欧美日韩国产mv在线观看视频| 日韩三级视频一区二区三区| 一级片'在线观看视频| 国产视频一区二区在线看| 91av网站免费观看| 美女福利国产在线| 日本黄色日本黄色录像| 欧美国产精品一级二级三级| 午夜福利视频在线观看免费| 久久久精品国产亚洲av高清涩受| 美女视频免费永久观看网站| 亚洲人成电影观看| 国产老妇伦熟女老妇高清| 黄色视频,在线免费观看| 亚洲精品久久午夜乱码| 国产不卡av网站在线观看| 多毛熟女@视频| 亚洲中文字幕日韩| 51午夜福利影视在线观看| 制服人妻中文乱码| 成年动漫av网址| 老司机午夜十八禁免费视频| 久久久久久免费高清国产稀缺| 中亚洲国语对白在线视频| 亚洲精品国产av成人精品| 久久精品熟女亚洲av麻豆精品| 黑丝袜美女国产一区| 国产精品一区二区免费欧美 | a级片在线免费高清观看视频| 桃花免费在线播放| 国产一级毛片在线| 男女边摸边吃奶| 中文字幕精品免费在线观看视频| 国产欧美日韩精品亚洲av| 男女国产视频网站| 久久国产精品男人的天堂亚洲| 亚洲欧美成人综合另类久久久| 午夜两性在线视频| 纵有疾风起免费观看全集完整版| 咕卡用的链子| 国产欧美亚洲国产| 亚洲国产精品一区三区| 国产主播在线观看一区二区| 一本色道久久久久久精品综合| 丝袜在线中文字幕| 18在线观看网站| 成年人午夜在线观看视频| 精品福利永久在线观看| 国产精品一二三区在线看| 免费观看a级毛片全部| 嫩草影视91久久| 桃花免费在线播放| 午夜免费观看性视频| av天堂久久9| 十八禁网站网址无遮挡| 亚洲欧美一区二区三区久久| 99国产精品一区二区蜜桃av | 一级毛片电影观看| 国产淫语在线视频| 无遮挡黄片免费观看| 亚洲黑人精品在线| 真人做人爱边吃奶动态| 日韩,欧美,国产一区二区三区| 国产精品一区二区免费欧美 | 在线av久久热| 三上悠亚av全集在线观看| 91国产中文字幕| 韩国高清视频一区二区三区| av超薄肉色丝袜交足视频| 搡老熟女国产l中国老女人| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 两个人免费观看高清视频| 午夜久久久在线观看| 老熟妇乱子伦视频在线观看 | 中文字幕人妻丝袜制服| 黄色怎么调成土黄色| 久久中文看片网| 国产精品 欧美亚洲| h视频一区二区三区| 人妻久久中文字幕网| 国产亚洲av片在线观看秒播厂| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 最近最新中文字幕大全免费视频| 男男h啪啪无遮挡| 成年人免费黄色播放视频| 国产日韩一区二区三区精品不卡| 日韩电影二区| 可以免费在线观看a视频的电影网站| av线在线观看网站| 最黄视频免费看| 热re99久久精品国产66热6| 亚洲精品久久久久久婷婷小说| 欧美在线黄色| 午夜福利视频在线观看免费| svipshipincom国产片| 国产99久久九九免费精品| 国产老妇伦熟女老妇高清| 秋霞在线观看毛片| 宅男免费午夜| 欧美日韩一级在线毛片| 老司机亚洲免费影院| 国产亚洲av片在线观看秒播厂| 最近最新中文字幕大全免费视频| 老司机午夜十八禁免费视频| 老司机亚洲免费影院| 亚洲成人国产一区在线观看| 最黄视频免费看| 丝袜脚勾引网站| 国产男女超爽视频在线观看| 少妇的丰满在线观看| 丰满饥渴人妻一区二区三| 国产成人精品在线电影| 视频区图区小说| 亚洲精品国产av成人精品| 少妇被粗大的猛进出69影院| 亚洲精华国产精华精| 啦啦啦在线免费观看视频4| 欧美在线黄色| 两性夫妻黄色片| 丰满少妇做爰视频| 亚洲成人国产一区在线观看| 日韩大片免费观看网站| 亚洲精品中文字幕在线视频| 99国产精品一区二区蜜桃av | 国产精品二区激情视频| 欧美 亚洲 国产 日韩一| 十八禁高潮呻吟视频| 免费看十八禁软件| 国产成人a∨麻豆精品| 午夜福利视频在线观看免费| 国产精品偷伦视频观看了| 男男h啪啪无遮挡| 国产主播在线观看一区二区| 亚洲视频免费观看视频| 亚洲天堂av无毛| 美国免费a级毛片| 午夜福利乱码中文字幕| 日本一区二区免费在线视频| 国产精品久久久久久精品电影小说| 亚洲成人国产一区在线观看| 十分钟在线观看高清视频www| 亚洲欧美激情在线| 久久中文字幕一级| 淫妇啪啪啪对白视频 | 成人亚洲精品一区在线观看| 高清在线国产一区| 国产深夜福利视频在线观看| 黄频高清免费视频| 丰满迷人的少妇在线观看| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕| 无限看片的www在线观看| 亚洲男人天堂网一区| 国产极品粉嫩免费观看在线| 亚洲精品国产av成人精品| 男女床上黄色一级片免费看| 大码成人一级视频| 亚洲精品国产区一区二| 国产成人精品无人区| 精品高清国产在线一区| 免费少妇av软件| 一区二区av电影网| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜制服| 两个人看的免费小视频| 搡老岳熟女国产| 欧美日韩福利视频一区二区| 动漫黄色视频在线观看| 久久久久久久久久久久大奶| 国产亚洲精品一区二区www | 国产男人的电影天堂91| 久久久久久久大尺度免费视频| 欧美精品人与动牲交sv欧美| 亚洲av成人一区二区三| 国产精品二区激情视频| 国产成人啪精品午夜网站| 精品人妻1区二区| 国产欧美日韩一区二区精品| 下体分泌物呈黄色| 日韩欧美免费精品| 久久 成人 亚洲| 日韩欧美一区二区三区在线观看 | 脱女人内裤的视频| 精品一品国产午夜福利视频| 亚洲欧美成人综合另类久久久| 韩国高清视频一区二区三区| 久久99一区二区三区| 午夜91福利影院| 99国产精品免费福利视频| 亚洲欧洲日产国产| 成人国产一区最新在线观看| 18在线观看网站| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 国产欧美日韩一区二区三 | 久久精品国产综合久久久| 美女高潮喷水抽搐中文字幕| 久久久精品94久久精品| 欧美97在线视频| 日本vs欧美在线观看视频| 蜜桃国产av成人99| a级毛片在线看网站| 国产精品久久久久成人av| 国产一区二区激情短视频 | 欧美av亚洲av综合av国产av| 亚洲国产av新网站| 亚洲精品国产区一区二| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说| 久久国产精品大桥未久av| 精品少妇内射三级| 最近中文字幕2019免费版| 淫妇啪啪啪对白视频 | 亚洲精品日韩在线中文字幕| 色老头精品视频在线观看| 久久中文字幕一级| 9191精品国产免费久久| 一个人免费在线观看的高清视频 | 啦啦啦免费观看视频1| av网站在线播放免费| 欧美在线黄色| 搡老乐熟女国产| 99re6热这里在线精品视频| 久久国产亚洲av麻豆专区| 欧美激情高清一区二区三区| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 狂野欧美激情性bbbbbb| 18在线观看网站| 欧美日韩亚洲高清精品| 人人妻人人澡人人看| 97在线人人人人妻| 一本—道久久a久久精品蜜桃钙片| 国产真人三级小视频在线观看| 色综合欧美亚洲国产小说| 在线永久观看黄色视频| 青青草视频在线视频观看| 性高湖久久久久久久久免费观看| 亚洲精品国产av成人精品| 母亲3免费完整高清在线观看| 亚洲成人免费电影在线观看| 国产一区二区三区av在线| 国产日韩欧美亚洲二区|