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

    激波/湍流邊界層干擾壓力脈動特性數(shù)值研究1)

    2021-11-09 06:26:08童福林段俊亦周桂宇李新亮
    力學(xué)學(xué)報 2021年7期
    關(guān)鍵詞:物面外層邊界層

    童福林 段俊亦 周桂宇 李新亮

    * (中國空氣動力研究與發(fā)展中心空氣動力學(xué)國家重點實驗室,四川綿陽 621000)

    ? (中國科學(xué)院力學(xué)研究所高溫氣體動力學(xué)國家重點實驗室,北京 100190)

    ** (中國空氣動力研究與發(fā)展中心計算空氣動力研究所,四川綿陽 621000)

    ?? (中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京 100049)

    引言

    激波與湍流邊界層的相互作用問題廣泛存在于各類高速飛行器中.已有研究表明[1-2],在強激波作用下,干擾區(qū)內(nèi)壓力脈動峰值急劇升高,分離激波的非定常運動導(dǎo)致物面脈動壓力低頻能量的急劇增強,這將使得飛行器結(jié)構(gòu)出現(xiàn)振動疲勞問題,進而嚴重影響飛行安全.因此,深入研究激波/湍流邊界層干擾區(qū)內(nèi)壓力脈動特性對認識和理解激波與湍流邊界層的相互作用機制非常重要.

    自20 世紀70 年代以來,國內(nèi)外學(xué)者一直致力于對該問題的研究.早期研究主要以風(fēng)洞試驗為主,受限于實驗測量技術(shù),風(fēng)洞試驗研究多局限于物面壓力脈動信號,并以此獲得干擾區(qū)內(nèi)典型流動結(jié)構(gòu)的特征信息,如分離激波非定常運動特性、湍流結(jié)構(gòu)、速度場等.Settle 等[3]研究了激波強度對壓縮拐角物面壓力脈動強度的影響規(guī)律.Ringuette 等[4]研究表明,低雷諾數(shù)下,干擾區(qū)內(nèi)物面壓力脈動強度降低而間歇性增強.在分離激波低頻振蕩方面,Erengil等[5]采用條件取樣和變窗口系綜平均技術(shù)分析了物面壓力分布與分離激波非定常運動的相關(guān)性.結(jié)果表明,分離泡內(nèi)的低頻壓力脈動與激波位置密切相關(guān):隨著激波往下游移動,分離泡內(nèi)壓力升高,反之則降低.Dolling 等[6]比較了不同構(gòu)型誘導(dǎo)的激波/湍流邊界層干擾區(qū)物面壓力脈動特性,研究發(fā)現(xiàn),分離激波的低頻振蕩機制主要來源于下游流場的局部或全局脈動,其過零頻率與上游湍流邊界層的特征頻率無關(guān).此外,Selig 等[7]基于物面壓力與邊界層質(zhì)量通量的時空關(guān)聯(lián)分析,給出了干擾區(qū)內(nèi)湍流結(jié)構(gòu)的演化規(guī)律.Piponniau 等[8]采用本征正交分解和線性隨機估計方法,通過實驗測得的物面壓力脈動信息對速度場進行了重構(gòu).

    由于風(fēng)洞試驗?zāi)軌虻玫降膲毫π畔O為有限,為此國內(nèi)外學(xué)者針對干擾區(qū)內(nèi)壓力脈動特性也開展了大量的高精度數(shù)值模擬研究.特別是近些年,直接數(shù)值模擬方法(direct numerical simulation,DNS)已經(jīng)在激波/湍流邊界層干擾問題方面取得了長足進展,其中Pirozzoli 等[9-11]的工作最具代表性.Pirozzoli等[9]分析比較了干擾區(qū)內(nèi)物面壓力脈動與剪切層壓力脈動的差異,并基于聲共振機制成功解釋了分離激波的非定常運動現(xiàn)象.Bernardini 等[10]著重探討了干擾區(qū)內(nèi)物面壓力脈動特性的演化規(guī)律,如脈動強度、時空關(guān)聯(lián)、對流速度和頻譜等,發(fā)現(xiàn)物面壓力脈動對流速度明顯降低而其展向尺度則顯著增強.隨后,Volpiani 等[11]進一步考察了壁面溫度對物面壓力脈動特性的影響規(guī)律.與此同時,國內(nèi)童福林等針對物面壓力脈動特性也開展了較為詳細的直接數(shù)值模擬研究,歸納了激波強度[12]、壁面溫度[13]、膨脹效應(yīng)[14-15]等因素的影響規(guī)律.

    盡管在激波/湍流邊界層干擾區(qū)壓力脈動特性研究方面,國內(nèi)外已開展了大量的風(fēng)洞試驗和數(shù)值研究工作,但是這些研究大多關(guān)注物面壓力脈動特性及其與分離激波非定常運動的相互關(guān)系方面,而關(guān)于邊界層內(nèi)壓力脈動在干擾區(qū)內(nèi)演化機制的認識仍然是很有限的.實際上,Duan 等[16]的研究結(jié)果表明,邊界層內(nèi)壓力脈動結(jié)構(gòu)特征與物面壓力脈動完全不同,其特征頻率沿法向?qū)⒔档投卣鞒叨葘⒓眲≡龃?然而,目前學(xué)術(shù)界對邊界層內(nèi)外層大尺度壓力脈動結(jié)構(gòu)與分離激波的相互作用機制仍然缺乏充分的認識,亟待開展相關(guān)研究工作.

    本文采用直接數(shù)值模擬方法研究M∞=2.25,Reθ=3567 來流下激波角β=33.2°的入射激波/平板湍流邊界層干擾區(qū)內(nèi)壓力脈動特性.為了便于比較和驗證結(jié)果,計算參數(shù)的選取與Dupont 等[17]的實驗和Fang 等[18]的DNS 相近.計算結(jié)果與前人的DNS 結(jié)果及風(fēng)洞實驗數(shù)據(jù)進行了仔細驗證,著重探討干擾區(qū)外層壓力脈動與物面壓力脈動演化機制的差異,特別是激波干擾對脈動強度、功率譜、特征尺度及時空關(guān)聯(lián)特性等的影響規(guī)律.

    1 DNS 計算

    控制方程為曲線坐標系(τ,ξ,η,ζ)下的三維可壓縮無量綱Navier-Stokes 方程組

    其中J?1為直角坐標系(x,y,z)變換為曲線坐標系(ξ,η,ζ) 的Jacobian 矩陣,U為守恒變量,F,G和H分別為(x,y,z)方向上對應(yīng)的無黏通量,Fv,Gv和Hv分別為(x,y,z)方向上對應(yīng)的黏性通量,具體表達式參見文獻[19].本文采用的解算器是高精度有限差分軟件OPENCFD-SC,該求解器目前已在壓縮拐角[19-20]、膨脹角入射激波/湍流邊界層干擾[14-15,21]等復(fù)雜流動問題中得到廣泛應(yīng)用和驗證確認,可以保證計算結(jié)果的可靠性.

    對流項計算時,采用Martin 等[22]優(yōu)化構(gòu)造的WENO_SYMBO_LMT 格式以及Steger-Warming流通量分裂方法,這樣能夠在極大抑制激波間斷區(qū)的強數(shù)值振蕩同時也保證了對邊界層內(nèi)不同尺度湍流脈動結(jié)構(gòu)的高分辨率捕捉.黏性項采用八階中心差分格式進行離散,時間方向上推進采用三階Runge-Kutta 方法.

    圖1 為計算模型示意圖,來流方向為圖中從左往右,來流參數(shù)如下:馬赫數(shù)M∞=2.25,基于單位長度來流雷諾數(shù)Re∞=2.5×104mm?1,靜溫T∞=169.44 K.模型長度Lx×Ly×Lz=137.6 mm×12.7 mm×4.4 mm,這里L(fēng)x,Ly和Lz分別為流向、法向和展向長度,坐標系原點取為計算域入口處.采用與文獻[21]類似的方法生成入射激波,在計算域上邊界xin=82.2 mm左右兩側(cè)分別設(shè)置為波前和波后參數(shù),即:在xin左側(cè)流動參數(shù)取為自由來流參數(shù),而在xin右側(cè)則按照激波關(guān)系式Rankine-Hugoniot 給出.這里激波角β取為33.2°,圖中入射激波在壁面上的名義入射點為xs=101.6 mm.

    圖1 計算模型Fig.1 Computational model

    計算網(wǎng)格點數(shù)為Nx×Ny×Nz=3700×300×250 (流向×法向×展向).圖2 為計算網(wǎng)格示意圖,為了便于展示,流向間隔了10 個點,法向間隔了5 個點.網(wǎng)格生成采用代數(shù)解析方法,這里流向網(wǎng)格點采用三段分布:上游轉(zhuǎn)捩區(qū)0 mm

    圖2 計算網(wǎng)格示意圖Fig.2 Sketch of the computational grid

    邊界條件設(shè)置如下:計算域入口處取為相同來流下的層流解(laminar inflow),具體分布如圖3 所示.壁面取為無滑移和等溫壁,壁溫Tw=321.9 K,上邊界采用簡單無反射邊界條件.出口處采用超聲速出口無反射條件及緩沖區(qū),以便消除擾動波反射帶來的影響.計算域展向兩側(cè)采用周期性邊界條件.

    圖3 入口處層流剖面Fig.3 Laminar profile at inlet

    計算時,通過在上游壁面添加多頻吹吸擾動從而轉(zhuǎn)捩生成充分發(fā)展湍流邊界層,進而與下游入射激波產(chǎn)生相互作用.這里,法向擾動速度Vbs的具體表達式如下[18]

    DNS 計算時取lmax=10,mmax=5,擾動幅值A(chǔ)=0.2,擾動頻率?=0.628U∞/δ,這里δ表示參考點處邊界層厚度,下文類似.擾動帶起始點分別為xa=7.62 mm 和xb=20.32 mm.式中相位差φl和φm取為0 到1 的隨機數(shù).

    表1 給出了參考點xref處湍流邊界層參數(shù)與Fang 等[18]DNS 結(jié)果比較情況,表中Reδ,Reδ*,Reθ分別為基于邊界層厚度、位移厚度和動量厚度的雷諾數(shù),Cf為摩阻系數(shù).兩者較為接近,這也表明本文計算結(jié)果的準確性.為了進一步驗證本文DNS 結(jié)果的可靠性,還與以往數(shù)值模擬結(jié)果和風(fēng)洞實驗數(shù)據(jù)進行了比較,包括湍流邊界層的速度剖面、脈動強度,物面壓力頻譜以及平均壓力分布.

    表1 參考點xref 處湍流邊界層參數(shù)Table 1 Turbulent boundary layer parameters at xref

    首先,圖4~ 圖6 分別給出了上游湍流邊界層在參考點xref處的統(tǒng)計特性.圖4 為平均流向速度剖面,定義為

    圖4 平均流向速度剖面Fig.4 Mean streamwise velocity profile at xref

    式中ρ和u分別為密度和流向速度,上標 ? 表示平均量.需要特別說明的是,本文中平均指的是時間和展向平均.速度分布在黏性底層y+<10 符合線性律,在30

    圖5 湍流脈動強度Fig.5 Turbulence intensities at xref

    圖6 物面壓力功率譜Fig.6 Power spectral density for wall pressure at xref

    此外,圖7 還給出了干擾區(qū)內(nèi)物面壓力分布情況,其中物面壓力P*和流向坐標x*分別定義為

    圖7 物面壓力分布Fig.7 Distribution of wall pressure

    式中,P1是反射激波的波后壓力值,xr*為轉(zhuǎn)換后的反射激波腳對應(yīng)的平均流向位置,這里采用與Fang等[18]類似的方法,即(P1+P∞)/2 對應(yīng)的流向位置.圖中紅色虛線代表壓力分布的無黏解.可見,計算得到的壓力分布與Fang 等[18]的DNS 結(jié)果及Dupont 等[17]的實驗結(jié)果相符較好.在黏性干擾作用下,干擾區(qū)內(nèi)壓力梯度變化較為緩和且在下游逐步趨近于無黏解.

    2 脈動強度

    圖8 給出了上游湍流邊界層(turbulent boundary layer,TBL)壓力脈動強度沿法向分布情況.可見,峰值位置出現(xiàn)在近壁區(qū),約為y/δ=0.04,其量值沿物面法向急劇降低,在邊界層外緣,約為近壁峰值的28%.采用當(dāng)?shù)匚锩婕羟袘?yīng)力τw進行無量綱化后,計算得到的脈動分布與Bernardini 等[28]和Duan 等[16]的DNS 結(jié)果符合較好.

    圖8 上游湍流邊界層壓力脈動強度分布Fig.8 Distribution of pressure fluctuation intensity for the incoming TBL

    從圖9 可以看出,與上游湍流邊界層相比,干擾區(qū)內(nèi)壓力脈動強度呈現(xiàn)截然不同的分布規(guī)律.從量值上來看,脈動峰值呈現(xiàn)先急劇升高,這主要是激波的增強作用,隨后在下游再附區(qū)伴隨著擾動邊界層的恢復(fù)則呈現(xiàn)緩慢衰減的流向演化趨勢.

    圖9 干擾區(qū)壓力脈動強度云圖Fig.9 Contour of pressure fluctuation intensity in the interaction region

    為了定量描述這一現(xiàn)象,圖10 給出了物面壓力脈動強度沿流向的分布情況.計算得到脈動強度與Fang 等[20]的DNS 結(jié)果符合良好,兩者在整體往下偏移0.02 的情況下與Dupont 等[17]的實驗數(shù)據(jù)也較為吻合.這主要是由于壓力傳感器截止頻率太小導(dǎo)致實驗數(shù)據(jù)明顯偏低[17].圖中符號S和R分別表示分離點和再附點位置,通過平均摩阻的過零點確定.可見,分離區(qū)內(nèi)物面壓力脈動急劇增強,其峰值約為上游的4.3 倍,隨后在下游再附區(qū)脈動峰值緩慢降低,在x*?xr*=12 時,仍約為上游的2.2 倍.另外,從峰值位置來看,其在分離區(qū)呈現(xiàn)逐漸遠離壁面的態(tài)勢,這與激波誘導(dǎo)生成分離區(qū)剪切層相關(guān),如圖9中平均音速線所示(粉色點劃線).強脈動主要出現(xiàn)在干擾區(qū)內(nèi)壓縮波系、剪切層附近等強剪切區(qū)域.此后,在下游再附邊界層,其峰值位置仍集中在外層區(qū)域.

    圖10 物面壓力脈動強度分布Fig.10 Distribution of wall pressure fluctuation intensity

    3 功率譜

    沿外層y/δ=0.8 和物面分別選取5 個流向站位進行對比分析,其中站位1 位于上游湍流邊界層內(nèi),站位2 位于分離區(qū)內(nèi),站位3~ 5 位于下游再附邊界層內(nèi),各站位具體位置見圖9.如無特別說明,下文中外層均為y/δ=0.8 處.

    圖11 比較了上游湍流邊界層站位1 壓力脈動的預(yù)乘功率譜ωφ(ω),其中ω和φ(ω)分別為角頻率和功率譜密度.功率譜密度的計算采用Pwelch 方法和Hamming 窗函數(shù),同時對分段數(shù)據(jù)按50%重疊處理.為了便于定量比較,圖中還采用當(dāng)?shù)孛}動強度的平方進行了歸一化處理.顯而易見,外層壓力脈動的時間尺度明顯要大于物面處,前者峰值頻率約為ωδ/U∞=3.0,后者約為ωδ/U∞=12.3,這與Duan 等[13]的研究結(jié)論是一致的.這是由于外層壓力脈動主要表征為大尺度脈動結(jié)構(gòu).

    圖11 站位1 壓力脈動預(yù)乘功率譜Fig.11 Pre-multiplied power spectral density for pressure fluctuations at station 1

    圖12 分別給出了干擾區(qū)外層和物面壓力脈動預(yù)乘功率譜的變化情況.從兩者的比較來看,以下幾個方面值得特別關(guān)注.首先,相較于上游站位1,分離區(qū)站位2 處兩者都出現(xiàn)了峰值頻率向更低頻區(qū)偏移的現(xiàn)象,這與激波的低頻振蕩運動密切相關(guān).

    圖12 干擾區(qū)內(nèi)不同站位壓力脈動預(yù)乘功率譜Fig.12 Pre-multiplied power spectral density for pressure fluctuations at various streamwise locations

    如圖12(a)和圖12(b)所示,外層和物面壓力峰值頻率分別約為ωδ/U∞=0.7 和ωδ/U∞=4.4.以往研究表明[29],在分離激波低頻振蕩的作用下,分離區(qū)內(nèi)物面壓力脈動的時間尺度通常比上游湍流邊界層低1~ 2 個數(shù)量級.本文計算結(jié)果進一步驗證了該結(jié)論,同時結(jié)果還表明,外層壓力脈動也同樣符合這一規(guī)律.在功率譜的定性分布方面,顯然,外層壓力脈動的變化較物面壓力更為劇烈,這很可能是由于外層站位非常靠近分離激波,受其非定常運動的影響更為直接,而物面壓力處于分離泡內(nèi)(見圖9),影響則要弱得多.

    從再附區(qū)站位3~ 5 的結(jié)果來看,外層壓力脈動呈現(xiàn)快速恢復(fù)的態(tài)勢,且峰值頻率向高頻區(qū)偏移,逐步趨近于上游站位1,這與物面壓力脈動的演化過程則完全不同.從圖12(a)中可以清楚看到,再附邊界層內(nèi)物面壓力脈動的低頻能量占比仍相對較高,峰值頻率依然維持在較低頻區(qū),并沒有出現(xiàn)往高頻區(qū)恢復(fù)的現(xiàn)象.造成兩者該差異的原因可以認為主要有兩方面:一方面,上游激波誘導(dǎo)分離形成了剪切層,分離泡上方剪切層內(nèi)不穩(wěn)定渦結(jié)構(gòu)沿流向急劇發(fā)展,并進入到下游再附邊界層內(nèi),這些不穩(wěn)定結(jié)構(gòu)極大地抑制了物面壓力脈動結(jié)構(gòu)的恢復(fù);另一方面,外層壓力脈動法向位置相對較高,盡管在站位2 處激波干擾對其影響顯著,但是過激波后外層壓力脈動受近壁區(qū)流動結(jié)構(gòu)的影響則要小得多,主要表征為外層大尺度脈動結(jié)構(gòu)的破碎過程,更為詳細的定量分析可見下一節(jié)的兩點相關(guān)統(tǒng)計結(jié)果.

    4 兩點相關(guān)

    通過兩點相關(guān)性考察干擾區(qū)外層壓力脈動擬序結(jié)構(gòu)長度尺度的變化規(guī)律,且定量分析比較其與物面壓力脈動結(jié)構(gòu)尺度的差異.兩點相關(guān)系數(shù)Cpp(Δx,Δz)定義為

    式中p′為雷諾平均后得到的壓力脈動,Δx和Δz分別代表流向間距和展向間距,x0為特征點流向位置,yref分別取物面和外層,上橫線代表時間和展向平均.

    圖13 分別給出了上游湍流邊界層外層和物面脈動壓力兩點相關(guān)系數(shù)的分布情況.以往的零壓力梯度平板結(jié)果表明[28],小間距時脈動壓力相關(guān)性曲線可近似為圓形,這是小尺度壓力脈動結(jié)構(gòu)的各向同性所決定的;大間距時,大尺度脈動結(jié)構(gòu)的各向異性導(dǎo)致相關(guān)性曲線在展向出現(xiàn)了急劇拉伸,近似為橢圓形分布,表征了展向拉長的大尺度結(jié)構(gòu).顯然,圖13(a)中的結(jié)果也完全符合該規(guī)律.然而,對于外層壓力脈動,如圖13(b)所示,整體來看,外層脈動結(jié)構(gòu)的流向尺度和展向尺度都明顯大于物面壓力脈動.這里以相關(guān)性系數(shù)0.3 為參考確定其特征長度,外層結(jié)構(gòu)的流向和展向尺度分別約為0.7δ和1.0δ,而物面壓力脈動則分別約為0.14δ和0.2δ.此外,盡管外層相關(guān)性等值線小間距時仍趨近于圓形,但隨著間距增大,等值線在流向和展向都被急劇拉伸,展向明顯強于流向,這表明外層壓力脈動擬序結(jié)構(gòu)仍以展向大尺度結(jié)構(gòu)為主,但其特征尺度要大得多.

    圖13 站位1 壓力脈動兩點相關(guān)系數(shù)Fig.13 Two-point correlation coefficient of pressure fluctuations at station 1

    圖14(a)和圖14(b)分別給出了站位2 物面和外層壓力脈動的兩點相關(guān)等值線,其中等值線從外到內(nèi)依次為0.3,0.5,0.7 和0.9.為了便于比較說明,圖中也給出了站位1 的結(jié)果.從分布規(guī)律上來看,等值線形態(tài)基本類似于上游湍流邊界層,仍然以展向拉伸的橢圓形分布為主,但整體呈現(xiàn)向外擴張的趨勢,這表明激波干擾作用增大了外層和物面壓力脈動結(jié)構(gòu)尺度,特別是在展向.需要特別關(guān)注的是,從圖14(b)中可以清楚看到,同樣以等值線0.3 為參考,站位1 處的展向尺度約為1.0δ,而站位2 處增大到約2.0δ,而流向尺度兩者基本維持在0.7δ.這很可能是由于站位2 外層位于入射激波附近(見圖9),激波的展向結(jié)構(gòu)增強了壓力脈動的展向相關(guān)性.更為細致的定量分析還有待深入研究.

    圖14 站位2 兩點相關(guān)系數(shù)(虛線:站位1)Fig.14 Two-point correlation coefficient of pressure fluctuations at station 2(dashed lines:station 1)

    如圖15 所示,再附區(qū)站位3~ 5 外層和物面壓力脈動的兩點相關(guān)等值線與上游湍流邊界層也基本類似.小間距時,曲線接近于圓形,這表明小尺度結(jié)構(gòu)并沒有發(fā)生實質(zhì)改變,仍以各向同性為主.從大尺度結(jié)構(gòu)的演化歷程來看,兩者差異顯著.研究發(fā)現(xiàn),大間距時,圖15(a)~ 圖15(c)中等值線以向外擴張為主,這表明再附區(qū)內(nèi)物面壓力大尺度脈動結(jié)構(gòu)的特征尺度仍在增大,而圖15(d)~ 圖15(f)等值線則以向內(nèi)收縮為主,對應(yīng)的是外層壓力大尺度脈動結(jié)構(gòu)特征尺度的減小.

    圖15 站位3~ 5 兩點相關(guān)系數(shù)(上:物面;下:y/δ=0.8;虛線:站位1)Fig.15 Two-point correlation coefficient at station 3~ 5 (top:wall;down:y/δ=0.8;dashed lines:station 1)

    為了定量描述該現(xiàn)象,圖16 進一步給出了干擾區(qū)內(nèi)物面及外層壓力脈動積分尺度的分布.積分尺度可以用來表征脈動場最大尺度結(jié)構(gòu)的特征尺度,這里流向和展向積分尺度分別定義為[10]

    可見,外層壓力脈動流向和展向積分尺度遠大于物面壓力脈動,這與Duan 等[16]的研究結(jié)果是一致的.從整體變化規(guī)律上來看,外層壓力Λx和Λz在站位2 處急劇增大,隨后在站位3~ 5 則又緩慢減小,而物面壓力的Λx和Λz整體上呈現(xiàn)逐步增大的趨勢,這與圖15 中相關(guān)性分析結(jié)果是相符的.從量值變化上來看,以圖16(b)中的展向積分尺度Λz為例,對于外層壓力,在激波的增強作用下,站位1 處Λz=0.73δ增大到站位2 處Λz=1.17δ,隨后在再附區(qū)又緩慢衰減到站位5 的Λz=1.02δ;而物面壓力,從站位1 處Λz=0.26δ逐步增大到站位5 處的Λz=0.82δ,后者約為前者的3.15 倍.可以看到,在站位5 處,此時再附區(qū)內(nèi)物面壓力脈動結(jié)構(gòu)的展向和流向積分尺度已非常接近于外層的積分尺度.綜上所述,可以認為這是由于外層壓力脈動大尺度結(jié)構(gòu)在過激波后呈現(xiàn)快速破碎的過程,而物面上大尺度結(jié)構(gòu)仍在緩慢增長,這也解釋了圖12 中物面壓力脈動低頻能量的相對升高而外層壓力脈動峰值頻率往高頻區(qū)的偏移現(xiàn)象.

    圖16 壓力脈動積分尺度分布Fig.16 Distribution of integral scale for pressure fluctuations

    5時空關(guān)聯(lián)

    進一步分析干擾區(qū)外層和物面壓力脈動的時空關(guān)聯(lián)特性,從而得到向下游傳播的壓力波的對流速度.

    圖17 分別給出了站位1 處物面和外層壓力的時空關(guān)聯(lián)云圖,圖中時空關(guān)聯(lián)系數(shù)Rpp(Δx,Δt)定義為

    式中Δt為延遲時間.以往研究結(jié)果表明[28,30],由于脈動結(jié)構(gòu)往下游傳播,物面壓力時空關(guān)聯(lián)云圖呈現(xiàn)傾斜的類橢圓形分布,如圖17(a)所示,本文計算結(jié)果也完全符合該結(jié)論.相較于物面壓力,從圖17(b)中可以看到,外層壓力時空關(guān)聯(lián)云圖的橢圓率和傾斜角都明顯增大.

    圖17 站位1 壓力脈動時空關(guān)聯(lián)系數(shù)Fig.17 Space-time correlation coefficient for pressure fluctuations at station 1

    圖18 和圖19 分別給出了分離區(qū)站位2 和再附區(qū)站位3~ 5 壓力脈動時空關(guān)聯(lián)系數(shù)與站位1 (黑色虛線)的比較情況.如圖18(a)所示,可以發(fā)現(xiàn),相較于上游湍流邊界層,分離區(qū)物面脈動壓力時空關(guān)聯(lián)曲線更為聚集,同時傾斜角也急劇減小,但整體形態(tài)仍滿足橢圓形(只是等值線0.1 較為雜亂).值得注意的是,對于站位2 外層脈動壓力,顯然,圖18(b)中等值線形態(tài)變化非常劇烈,已完全不符合橢圓形分布,特別是等值線0.1 和0.3,在時間方向上可近似為直線分布.造成這個現(xiàn)象的原因很可能是由于站位2 外層處于強激波區(qū)(見圖9),在分離激波的強干擾下,此時泰勒凍結(jié)假設(shè)完全失效,脈動壓力在時空上表征為強非線性.此外,從圖19 中還可以看到,無論是再附區(qū)外層還是物面,其脈動壓力時空關(guān)聯(lián)等值線整體形態(tài)與上游站位1 都較為類似.從站位3~站位5,強相關(guān)等值線被急劇拉長壓縮,且傾斜角增大,呈現(xiàn)逐步逼近上游站位1 的演化過程.從定性比較的角度來看,外層脈動壓力時空關(guān)聯(lián)曲線的恢復(fù)速度明顯要快于物面.

    圖18 站位2 壓力脈動時空關(guān)聯(lián)系數(shù)(虛線:站位1)Fig.18 Space-time correlation coefficient for pressure fluctuations at station 2 (dashed lines:station 1)

    圖19 站位3~ 5 壓力脈動時空關(guān)聯(lián)系數(shù)(上:物面;下:y/δ=0.8;虛線:站位1)Fig.19 Space-time correlation coefficient at station 3~ 5 (top:wall;down:y/δ=0.8;dashed lines:station 1)

    壓力脈動對流速度的計算方法如下[16,31]:首先給定延遲時間Δt,隨后依據(jù)時空關(guān)聯(lián)曲線通過

    獲得流向間距rx,進而得到對流速度Uc=rx/Δt.采用上述方法計算得到了上游湍流邊界層物面和外層脈動壓力對流速度隨延遲時間的分布情況.如圖20 和圖21 所示.總體來看,外層脈動壓力對流速度要遠大于物面脈動壓力,這與Duan 等[16]的研究結(jié)果也較為符合.另外,對于外層脈動壓力,可見,其對流速度對延遲時間的依賴性較小,基本維持在Uc=0.94U∞,與當(dāng)?shù)仄骄飨蛩俣冉咏?然而,物面脈動壓力對流速度隨延遲時間則變化劇烈,在小延遲時間(小尺度脈動)時,對流速度約為Uc=0.56U∞,而在大延遲時間(大尺度脈動) 時,其增大到Uc=0.73U∞.已有研究表明[13],這是由于物面大尺度脈動壓力與外層大尺度結(jié)構(gòu)相關(guān),而小尺度脈動壓力則主要受近壁區(qū)結(jié)構(gòu)影響,故前者對流速度明顯大于后者.

    從圖20 和圖21 中還可以看到,分離區(qū)物面脈動壓力對流速度降低,約為Uc=0.2U∞,隨后下游再附區(qū)則又逐漸增大,但仍小于上游湍流邊界層.Na等[32]在強逆壓力梯度平板的DNS 結(jié)果中也發(fā)現(xiàn)了類似規(guī)律,本文的計算結(jié)果進一步證實了該現(xiàn)象在入射激波/湍流邊界層干擾區(qū)內(nèi)的適用性.從圖20(b)中可以看到,采用當(dāng)?shù)赝鈱铀俣萓e對Uc進行歸一化后,分離區(qū)和再附區(qū)各站位處物面脈動壓力的對流速度分布均吻合較好,約為Uc=0.25Ue,與Bernardini 等[10]在跨聲速激波/湍流邊界層干擾下的研究結(jié)論也較為接近,但后者為初始分離(incipient separation)工況.如圖21 所示,干擾區(qū)外層不同站位處脈動壓力對流速度分布與上游結(jié)果較為類似,其對延遲時間的依賴性依然較弱,但其量值減小到約為Uc=0.4Ue.結(jié)合之前的功率譜和兩點相關(guān)分析,可以推測外層大尺度脈動結(jié)構(gòu)特征尺度的減小可能是造成該現(xiàn)象的重要因素.

    圖20 干擾區(qū)內(nèi)各站位物面壓力脈動對流速度Fig.20 Convection velocity for wall pressure fluctuations at various stations in the interaction region

    圖21 干擾區(qū)內(nèi)各站位外層壓力脈動對流速度Fig.21 Convection velocity for pressure fluctuations in the outer layer

    6 結(jié)論

    本文采用直接數(shù)值模擬方法研究了來流馬赫數(shù)2.25 下激波角33.2°的入射激波對平板湍流邊界層壓力脈動特性的影響規(guī)律,詳細地分析了外層和物面壓力脈動場的典型統(tǒng)計特征,如脈動強度、功率譜密度、兩點相關(guān)和時空關(guān)聯(lián)特性等,得到以下結(jié)論:

    (1) 數(shù)值模擬準確捕捉到了干擾區(qū)內(nèi)脈動壓力.研究發(fā)現(xiàn),分離區(qū)外層脈動壓力以低頻能量為主,高頻能量占比較小.與物面脈動壓力不同的是,再附區(qū)外層脈動壓力峰值頻率逐漸往高頻區(qū)偏移,而物面壓力脈動能量仍維持在相對較低頻率.

    (2) 激波干擾對壓力脈動結(jié)構(gòu)有顯著的增強作用,干擾區(qū)內(nèi)展向積分尺度明顯大于流向積分尺度.積分尺度在外層和物面的演化歷程差異明顯,前者經(jīng)歷了先急劇增長后緩慢衰減的過程,而后者整體上呈現(xiàn)增長趨勢.

    (3) 激波干擾極大地降低了壓力脈動場的對流速度.時空關(guān)聯(lián)結(jié)果表明,干擾區(qū)下游外層壓力波的對流速度仍明顯高于物面,前者約為0.4Ue,而后者約為0.25Ue.

    猜你喜歡
    物面外層邊界層
    一種溶液探測傳感器
    傳感器世界(2022年4期)2022-11-24 21:23:50
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    讓吸盤掛鉤更牢固
    一種購物袋
    科技資訊(2016年6期)2016-05-14 13:09:55
    新型單面陣自由曲面光學(xué)測量方法成像特性仿真
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    專題Ⅱ 物質(zhì)構(gòu)成的奧秘
    “人”字變身
    賽車空動優(yōu)化的秘密
    汽車之友(2014年21期)2014-11-03 17:33:45
    非特征邊界的MHD方程的邊界層
    一级毛片我不卡| 乱码一卡2卡4卡精品| 久久九九热精品免费| 亚洲七黄色美女视频| 亚洲av一区综合| 高清日韩中文字幕在线| 亚洲综合色惰| 此物有八面人人有两片| 看免费成人av毛片| 国产高清三级在线| av在线老鸭窝| 日韩欧美国产在线观看| 亚洲成a人片在线一区二区| 亚洲婷婷狠狠爱综合网| 欧美成人a在线观看| 国产成人a∨麻豆精品| 亚洲av第一区精品v没综合| 天堂动漫精品| 97碰自拍视频| 美女大奶头视频| 国产一区二区三区av在线 | 亚洲三级黄色毛片| 亚洲av成人精品一区久久| 免费黄网站久久成人精品| 精品99又大又爽又粗少妇毛片| 极品教师在线视频| 如何舔出高潮| 久久久国产成人精品二区| 日本黄色片子视频| 亚洲精品456在线播放app| 麻豆成人午夜福利视频| 日韩亚洲欧美综合| 国产亚洲av嫩草精品影院| 国产精华一区二区三区| 欧美成人a在线观看| 精品久久久久久久末码| 午夜福利视频1000在线观看| 久久国内精品自在自线图片| 久久午夜福利片| 一区二区三区免费毛片| 欧美激情在线99| 我要搜黄色片| 女人十人毛片免费观看3o分钟| 两个人视频免费观看高清| 91久久精品国产一区二区成人| 性插视频无遮挡在线免费观看| 色播亚洲综合网| 男女视频在线观看网站免费| videossex国产| 午夜精品在线福利| 麻豆成人午夜福利视频| 亚洲人与动物交配视频| 一级黄片播放器| 国产精品一二三区在线看| 如何舔出高潮| 久久久a久久爽久久v久久| 成年av动漫网址| 国产在视频线在精品| 成人综合一区亚洲| 看黄色毛片网站| 亚洲欧美成人精品一区二区| 内射极品少妇av片p| 在线观看66精品国产| 国内精品美女久久久久久| 国产乱人视频| 亚洲成人中文字幕在线播放| 草草在线视频免费看| 亚洲内射少妇av| 美女xxoo啪啪120秒动态图| 最近中文字幕高清免费大全6| 亚洲精品久久国产高清桃花| 最后的刺客免费高清国语| 一级毛片我不卡| 最近2019中文字幕mv第一页| 桃色一区二区三区在线观看| 国产精品一二三区在线看| 久久人人爽人人片av| 成年av动漫网址| 亚洲一区高清亚洲精品| 国产中年淑女户外野战色| 老司机午夜福利在线观看视频| 蜜桃亚洲精品一区二区三区| 久久久久久久亚洲中文字幕| 精品久久久久久久久久免费视频| 国产探花在线观看一区二区| 日韩中字成人| 一进一出好大好爽视频| 久久中文看片网| 少妇裸体淫交视频免费看高清| 国产成人影院久久av| 最近2019中文字幕mv第一页| 男人狂女人下面高潮的视频| 午夜久久久久精精品| 久久久久久久亚洲中文字幕| 人人妻人人看人人澡| 99久久九九国产精品国产免费| 精品不卡国产一区二区三区| 国产熟女欧美一区二区| 免费观看在线日韩| 村上凉子中文字幕在线| 国产精品av视频在线免费观看| 日韩欧美一区二区三区在线观看| 五月玫瑰六月丁香| 国产黄a三级三级三级人| 日韩强制内射视频| 观看美女的网站| 97热精品久久久久久| 麻豆国产av国片精品| 午夜精品国产一区二区电影 | 自拍偷自拍亚洲精品老妇| 一区福利在线观看| 长腿黑丝高跟| 91久久精品国产一区二区三区| 欧美日韩乱码在线| 一边摸一边抽搐一进一小说| 在线免费十八禁| 亚洲欧美精品综合久久99| 又黄又爽又刺激的免费视频.| 国产午夜精品论理片| 成人高潮视频无遮挡免费网站| 亚洲一区二区三区色噜噜| 亚洲精品亚洲一区二区| 国产伦精品一区二区三区视频9| 中国国产av一级| 亚洲精品影视一区二区三区av| 国产片特级美女逼逼视频| 少妇熟女欧美另类| 九九爱精品视频在线观看| 联通29元200g的流量卡| 国产免费男女视频| 国产精品亚洲一级av第二区| 亚洲av中文av极速乱| 草草在线视频免费看| 久久久久久大精品| 91在线精品国自产拍蜜月| 国产精品嫩草影院av在线观看| 男人狂女人下面高潮的视频| 亚洲国产精品国产精品| 啦啦啦韩国在线观看视频| 久久精品国产亚洲av涩爱 | 高清午夜精品一区二区三区 | 人人妻人人澡欧美一区二区| 久久婷婷人人爽人人干人人爱| 亚洲精品国产成人久久av| 色播亚洲综合网| 国产美女午夜福利| 校园春色视频在线观看| 国产综合懂色| 搞女人的毛片| 国产爱豆传媒在线观看| 国产精品av视频在线免费观看| 97在线视频观看| 18禁在线无遮挡免费观看视频 | 国产高清不卡午夜福利| 嫩草影视91久久| 精品乱码久久久久久99久播| 日本爱情动作片www.在线观看 | 国产精品福利在线免费观看| 九九热线精品视视频播放| 99riav亚洲国产免费| 久99久视频精品免费| 欧美3d第一页| 国产av一区在线观看免费| 伦精品一区二区三区| 熟女电影av网| 91av网一区二区| 床上黄色一级片| 午夜老司机福利剧场| 国产中年淑女户外野战色| 中文在线观看免费www的网站| 久久久久久久久久成人| 寂寞人妻少妇视频99o| 搡老岳熟女国产| 久久久精品94久久精品| 美女内射精品一级片tv| 99热6这里只有精品| 内地一区二区视频在线| 欧美日韩国产亚洲二区| 亚洲成人久久性| 99久久精品热视频| 在线观看美女被高潮喷水网站| 日韩高清综合在线| 春色校园在线视频观看| 欧美高清性xxxxhd video| 亚洲无线观看免费| 熟女电影av网| 日本色播在线视频| 床上黄色一级片| 男女视频在线观看网站免费| 国产日本99.免费观看| 亚洲人成网站高清观看| 白带黄色成豆腐渣| 两个人的视频大全免费| 亚洲美女视频黄频| 亚洲av第一区精品v没综合| 国内精品宾馆在线| 别揉我奶头 嗯啊视频| 久久鲁丝午夜福利片| 中文字幕av成人在线电影| 久久中文看片网| 卡戴珊不雅视频在线播放| 成人二区视频| 久久久久久久久中文| 精品熟女少妇av免费看| 男女下面进入的视频免费午夜| 国产国拍精品亚洲av在线观看| 99久国产av精品国产电影| 日本三级黄在线观看| 少妇熟女aⅴ在线视频| 99国产极品粉嫩在线观看| 变态另类丝袜制服| 午夜老司机福利剧场| 亚洲人成网站在线播| 亚洲欧美成人精品一区二区| 悠悠久久av| 国产精品,欧美在线| 一区二区三区四区激情视频 | 熟女人妻精品中文字幕| 国产在线精品亚洲第一网站| 夜夜看夜夜爽夜夜摸| 男人舔女人下体高潮全视频| 国产精品女同一区二区软件| 51国产日韩欧美| 又黄又爽又刺激的免费视频.| 国产精品亚洲一级av第二区| 丰满的人妻完整版| 性插视频无遮挡在线免费观看| 国产单亲对白刺激| 天堂av国产一区二区熟女人妻| 国产免费一级a男人的天堂| 久久中文看片网| 国模一区二区三区四区视频| 欧美成人精品欧美一级黄| 搞女人的毛片| 午夜福利18| 黄色视频,在线免费观看| 国产淫片久久久久久久久| 最近最新中文字幕大全电影3| 亚洲av五月六月丁香网| 美女免费视频网站| 免费观看人在逋| 国产亚洲欧美98| 免费av观看视频| 久久欧美精品欧美久久欧美| 日韩强制内射视频| 韩国av在线不卡| 久久久久久久久久成人| 精品一区二区三区人妻视频| ponron亚洲| 少妇的逼好多水| 插阴视频在线观看视频| 在线免费观看不下载黄p国产| 男人狂女人下面高潮的视频| 欧美潮喷喷水| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 国产一区二区在线av高清观看| 熟女人妻精品中文字幕| 中文字幕久久专区| 不卡一级毛片| 久久久久久久久久黄片| 欧美在线一区亚洲| 性欧美人与动物交配| 国产精品人妻久久久久久| 神马国产精品三级电影在线观看| 国产精品日韩av在线免费观看| 熟女人妻精品中文字幕| 国产在线精品亚洲第一网站| 精品一区二区三区人妻视频| 国产av不卡久久| 国产麻豆成人av免费视频| 日日摸夜夜添夜夜添小说| 成年女人看的毛片在线观看| 97超视频在线观看视频| 亚洲最大成人手机在线| 精品国内亚洲2022精品成人| 精品一区二区三区人妻视频| 99热只有精品国产| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 国产高清视频在线播放一区| 国产成年人精品一区二区| 一本一本综合久久| 99热精品在线国产| 中文字幕av成人在线电影| 91久久精品国产一区二区三区| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 一级a爱片免费观看的视频| 亚洲av第一区精品v没综合| 中出人妻视频一区二区| 国产精品乱码一区二三区的特点| 日韩精品有码人妻一区| 嫩草影院入口| 十八禁国产超污无遮挡网站| 女同久久另类99精品国产91| 丰满的人妻完整版| 成人无遮挡网站| 男女那种视频在线观看| 高清午夜精品一区二区三区 | 日韩欧美三级三区| 国产片特级美女逼逼视频| 国产精品电影一区二区三区| 在线观看av片永久免费下载| 一级毛片电影观看 | 中国国产av一级| 在线看三级毛片| 久久久久国产精品人妻aⅴ院| www.色视频.com| 在现免费观看毛片| 在线播放无遮挡| 美女免费视频网站| 午夜福利在线观看吧| 亚洲三级黄色毛片| 大型黄色视频在线免费观看| 国产精品精品国产色婷婷| 桃色一区二区三区在线观看| 久久精品国产亚洲av涩爱 | 久久这里只有精品中国| 国产精品久久久久久av不卡| 国产在视频线在精品| 欧美三级亚洲精品| 成年女人永久免费观看视频| 一进一出抽搐动态| 波多野结衣巨乳人妻| 无遮挡黄片免费观看| 天天躁日日操中文字幕| 精品久久久久久久久久免费视频| .国产精品久久| 人人妻人人澡欧美一区二区| 97碰自拍视频| 不卡一级毛片| 国产高清不卡午夜福利| 人妻制服诱惑在线中文字幕| 可以在线观看毛片的网站| 久久久久精品国产欧美久久久| 人人妻人人看人人澡| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品粉嫩美女一区| 午夜激情欧美在线| 国产精品99久久久久久久久| 亚洲欧美日韩无卡精品| 亚洲五月天丁香| 亚洲欧美日韩高清在线视频| 国内精品美女久久久久久| 俄罗斯特黄特色一大片| 欧美3d第一页| 亚洲av二区三区四区| 日本精品一区二区三区蜜桃| 热99re8久久精品国产| 人人妻人人看人人澡| 午夜福利18| 久久午夜福利片| 日韩大尺度精品在线看网址| 久久欧美精品欧美久久欧美| 午夜亚洲福利在线播放| 美女 人体艺术 gogo| 亚洲中文日韩欧美视频| 欧美bdsm另类| 日韩亚洲欧美综合| 18+在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 午夜免费男女啪啪视频观看 | 在线免费十八禁| 伊人久久精品亚洲午夜| 在线免费观看的www视频| 欧美丝袜亚洲另类| 亚洲欧美日韩东京热| 精品乱码久久久久久99久播| 欧美中文日本在线观看视频| 日本成人三级电影网站| 成人特级av手机在线观看| 免费黄网站久久成人精品| 又粗又爽又猛毛片免费看| 搡老妇女老女人老熟妇| 国产在线精品亚洲第一网站| 又粗又爽又猛毛片免费看| а√天堂www在线а√下载| 久久久精品大字幕| 国产午夜精品久久久久久一区二区三区 | 精品一区二区三区视频在线观看免费| 亚洲熟妇中文字幕五十中出| 久久久国产成人精品二区| 高清午夜精品一区二区三区 | 少妇裸体淫交视频免费看高清| av黄色大香蕉| 嫩草影视91久久| 亚洲精品影视一区二区三区av| 欧美最黄视频在线播放免费| 神马国产精品三级电影在线观看| 日韩av在线大香蕉| 又黄又爽又刺激的免费视频.| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 亚洲在线观看片| 最好的美女福利视频网| 性插视频无遮挡在线免费观看| 亚洲丝袜综合中文字幕| 欧美日韩精品成人综合77777| 99久国产av精品国产电影| 亚洲国产精品合色在线| 国产成人福利小说| 中文在线观看免费www的网站| 男女啪啪激烈高潮av片| 性插视频无遮挡在线免费观看| 中文字幕av在线有码专区| 色综合亚洲欧美另类图片| aaaaa片日本免费| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 亚洲第一电影网av| 此物有八面人人有两片| 男女视频在线观看网站免费| a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 欧美xxxx性猛交bbbb| 亚洲av五月六月丁香网| or卡值多少钱| 国产精品一区二区免费欧美| 91av网一区二区| 国产真实乱freesex| 精品不卡国产一区二区三区| 欧美日韩乱码在线| 青春草视频在线免费观看| 久久这里只有精品中国| 一级黄片播放器| 日本黄色视频三级网站网址| 天美传媒精品一区二区| 国产av麻豆久久久久久久| 午夜a级毛片| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 亚洲av第一区精品v没综合| 午夜视频国产福利| 精品一区二区免费观看| 麻豆av噜噜一区二区三区| 人妻夜夜爽99麻豆av| 性插视频无遮挡在线免费观看| 91狼人影院| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 变态另类成人亚洲欧美熟女| 激情 狠狠 欧美| 色尼玛亚洲综合影院| 亚洲成av人片在线播放无| 免费高清视频大片| 国产高清三级在线| 国产成人91sexporn| 久久午夜亚洲精品久久| 成人高潮视频无遮挡免费网站| 国产精品精品国产色婷婷| 一级毛片aaaaaa免费看小| 99热只有精品国产| 国产高清三级在线| 99久久中文字幕三级久久日本| 国产久久久一区二区三区| 国产精品综合久久久久久久免费| 日韩人妻高清精品专区| 99riav亚洲国产免费| 少妇熟女aⅴ在线视频| 一级毛片久久久久久久久女| 国产伦精品一区二区三区四那| 色吧在线观看| 精品久久久久久久久久久久久| 成年版毛片免费区| 国产一区二区三区在线臀色熟女| 国产av麻豆久久久久久久| 亚洲av美国av| 国产av一区在线观看免费| avwww免费| 欧美成人免费av一区二区三区| 最新中文字幕久久久久| 久久精品国产亚洲av天美| 日日撸夜夜添| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 可以在线观看毛片的网站| 国产一区二区在线观看日韩| 免费看日本二区| 亚洲精品日韩在线中文字幕 | 日韩av在线大香蕉| 18+在线观看网站| 成年女人永久免费观看视频| 尤物成人国产欧美一区二区三区| 欧美日韩精品成人综合77777| 国产精品伦人一区二区| 悠悠久久av| 可以在线观看的亚洲视频| 国产精品久久久久久亚洲av鲁大| 99热这里只有是精品在线观看| 日本色播在线视频| 中文字幕免费在线视频6| 日韩中字成人| 午夜老司机福利剧场| 不卡一级毛片| 精品人妻熟女av久视频| 国产精品一区二区性色av| 99在线视频只有这里精品首页| 日韩av不卡免费在线播放| 99热全是精品| 亚洲自拍偷在线| 狂野欧美白嫩少妇大欣赏| 亚洲精品在线观看二区| 国产中年淑女户外野战色| 亚洲精品日韩av片在线观看| 秋霞在线观看毛片| 亚洲一区二区三区色噜噜| 国产真实乱freesex| 成年女人看的毛片在线观看| 久久久色成人| 日韩欧美 国产精品| 我要搜黄色片| 天堂av国产一区二区熟女人妻| 国产精品精品国产色婷婷| 国产毛片a区久久久久| 久久6这里有精品| 国产精品1区2区在线观看.| 男插女下体视频免费在线播放| 悠悠久久av| 国产黄色视频一区二区在线观看 | 九九爱精品视频在线观看| 久久精品久久久久久噜噜老黄 | 国产精品爽爽va在线观看网站| 小蜜桃在线观看免费完整版高清| 精品久久国产蜜桃| 联通29元200g的流量卡| 少妇的逼好多水| 美女被艹到高潮喷水动态| 99在线人妻在线中文字幕| 国产欧美日韩精品一区二区| 久久久久久国产a免费观看| 亚洲中文日韩欧美视频| 乱人视频在线观看| 亚洲无线观看免费| 亚洲七黄色美女视频| 亚洲欧美日韩东京热| 美女xxoo啪啪120秒动态图| 欧美高清成人免费视频www| 狠狠狠狠99中文字幕| av在线蜜桃| 午夜福利成人在线免费观看| 床上黄色一级片| 成年版毛片免费区| 午夜激情福利司机影院| 久久精品国产清高在天天线| 亚洲性久久影院| 亚洲精品国产av成人精品 | av黄色大香蕉| 18禁在线播放成人免费| 国产一级毛片七仙女欲春2| 如何舔出高潮| 国产精品三级大全| .国产精品久久| 国产日本99.免费观看| 国产精品国产三级国产av玫瑰| 国产伦在线观看视频一区| 国内精品宾馆在线| 少妇被粗大猛烈的视频| 成人欧美大片| 亚洲婷婷狠狠爱综合网| av在线观看视频网站免费| 十八禁国产超污无遮挡网站| 91在线精品国自产拍蜜月| 国产高清视频在线播放一区| 亚洲精品乱码久久久v下载方式| 淫秽高清视频在线观看| 亚洲中文字幕日韩| 91精品国产九色| 麻豆国产97在线/欧美| 亚洲精品一卡2卡三卡4卡5卡| aaaaa片日本免费| 黑人高潮一二区| 国产精品乱码一区二三区的特点| 精品人妻熟女av久视频| 一进一出好大好爽视频| 欧美激情国产日韩精品一区| 最近在线观看免费完整版| 成人国产麻豆网| 亚洲欧美成人精品一区二区| 久久九九热精品免费| 少妇丰满av| 日韩欧美在线乱码| 久久九九热精品免费| 色噜噜av男人的天堂激情| 久久久久久久久中文| 美女cb高潮喷水在线观看| av国产免费在线观看| 99久久精品一区二区三区| 一级毛片我不卡| 18禁黄网站禁片免费观看直播| 亚洲av中文字字幕乱码综合| 高清毛片免费观看视频网站| av国产免费在线观看| 中文字幕av在线有码专区| 禁无遮挡网站| 国产不卡一卡二| 亚洲欧美成人精品一区二区| 男人和女人高潮做爰伦理| 一a级毛片在线观看| 亚洲欧美成人精品一区二区| 精品欧美国产一区二区三| 亚洲欧美成人综合另类久久久 | 国产av麻豆久久久久久久| 午夜老司机福利剧场| 97碰自拍视频| 免费黄网站久久成人精品| 天堂动漫精品| 看片在线看免费视频| 我要看日韩黄色一级片| 国产伦在线观看视频一区| 天美传媒精品一区二区| 国产一级毛片七仙女欲春2|