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

    尾跡對(duì)低壓渦輪吸力面附面層流動(dòng)的影響

    2022-11-28 13:40:36曹惠玲吳興爽譚天榮宋立旗劉福林
    航空發(fā)動(dòng)機(jī) 2022年5期
    關(guān)鍵詞:附面層葉柵尾跡

    曹惠玲,吳興爽,譚天榮,宋立旗,劉福林

    (中國民航大學(xué)航空工程學(xué)院,天津 300300)

    0 引言

    低壓渦輪(Low Pressure Turbine,LPT)的級(jí)數(shù)和葉片較多,是影響航空發(fā)動(dòng)機(jī)效率和質(zhì)量的重要部件之一。采用高負(fù)荷LPT葉片設(shè)計(jì)可有效地減少葉片數(shù)量,研究并理解LPT內(nèi)部流動(dòng)特性及其機(jī)理有助于LPT設(shè)計(jì)。Curits等[1]發(fā)現(xiàn)吸力面產(chǎn)生的氣動(dòng)損失占葉型損失的60%;Mayle[2]發(fā)現(xiàn)LPT在低雷諾數(shù)下工作時(shí),其吸力面附面層可能會(huì)發(fā)生分離,造成嚴(yán)重的氣動(dòng)損失使發(fā)動(dòng)機(jī)效率降低。此外,在LPT工作時(shí),上游葉片產(chǎn)生周期性非定常尾跡將影響下游葉片的附面層流動(dòng)。

    一方面,大量學(xué)者針對(duì)低雷諾數(shù)下的低壓渦輪非定常流動(dòng)進(jìn)行了研究。Kaszeta等[3]研究了尾跡對(duì)低壓渦輪葉柵轉(zhuǎn)捩和分離的影響,指出尾跡產(chǎn)生的湍流使湍流水平總體上升;Zhang等[4]結(jié)合雷諾數(shù)和湍流度研究了非定常尾跡對(duì)超高負(fù)荷低壓渦輪附面層的影響,并指出在低湍流低雷諾數(shù)下,尾跡與分離附面層的相互作用主要是由尾跡負(fù)射流誘導(dǎo)的卷起渦引起的;Wolff等[5]使用熱線研究了超高負(fù)荷低壓渦輪葉片的非定常尾跡誘導(dǎo)附面層發(fā)展特性,并指出非定常尾跡在低雷諾數(shù)下對(duì)分離泡的抑制有助于減小葉型損失。

    另一方面,學(xué)術(shù)界對(duì)于附面層內(nèi)部流動(dòng)特性與轉(zhuǎn)捩機(jī)理的研究也較為深入。Jacobs等[6]發(fā)現(xiàn)Klebanoff條紋因葉片前緣的“剪切遮蔽”效應(yīng)產(chǎn)生;Liu等[7]發(fā)現(xiàn)自由流湍流會(huì)在附面層內(nèi)誘導(dǎo)出弱Klebanoff條紋;Coull[8]發(fā)現(xiàn)尾跡放大Klebanoff條紋的前緣、最強(qiáng)條紋、尾緣分別以88%、70%、50%的主流速度向下游運(yùn)動(dòng);Wissink等[9]在尾跡掃掠分離剪切層時(shí),觀測到因開爾文-亥姆霍茲(Kelvin-Helmholtz,KH)不穩(wěn)定性誘導(dǎo)形成了全展向卷起渦結(jié)構(gòu);McAuliffe等[10]發(fā)現(xiàn)Klebanoff條紋會(huì)扭曲分離剪切層,并抑制全展向KH卷起渦結(jié)構(gòu)的形成。

    中國針對(duì)非定常尾跡影響葉片流動(dòng)特性問題開展了大量研究。王英峰等[11]采用頻譜分析的方法進(jìn)行了上游尾跡對(duì)葉片分離的影響研究;羅華玲等[12]研究了尾跡在葉柵中的輸運(yùn)及葉片排之間的時(shí)序效應(yīng);李虹楊等[13]應(yīng)用自行開發(fā)的CFD程序?qū)χ芷谛晕槽E作用下的渦輪葉柵附面層轉(zhuǎn)捩進(jìn)行了研究;Sun等[14-15]研究了上游尾跡對(duì)低壓渦輪葉片附面層發(fā)展的影響,并對(duì)比了采用三角棒和圓棒模擬上游尾跡對(duì)吸力面附面層流動(dòng)與轉(zhuǎn)捩的影響。但中國學(xué)者對(duì)附面層瞬態(tài)分析的研究卻較為有限。Cui等[16-17]、Sun等[18]、Nagabhushana等[19-20]雖對(duì)葉柵通道內(nèi)的瞬態(tài)結(jié)構(gòu)進(jìn)行了分析,但主要集中于LES模型在主流與端區(qū)中的應(yīng)用,對(duì)于尾跡誘導(dǎo)轉(zhuǎn)捩和自然轉(zhuǎn)捩相互轉(zhuǎn)換主導(dǎo)附面層轉(zhuǎn)捩這一過程并無完整的時(shí)空演化分析。

    本文在1個(gè)完整尾跡掃掠周期內(nèi)對(duì)尾跡在葉片前緣形成的瞬態(tài)結(jié)構(gòu)進(jìn)行了研究,并分析了瞬態(tài)結(jié)構(gòu)對(duì)吸力面附面層流動(dòng)特性與多模態(tài)轉(zhuǎn)捩的影響。

    1 試驗(yàn)與數(shù)值模擬

    1.1 試驗(yàn)設(shè)置

    試驗(yàn)研究依托低速葉柵試驗(yàn)臺(tái)進(jìn)行,帶上游圓棒的葉柵通道模型及葉柵測試試驗(yàn)臺(tái)照片如圖1所示,葉柵和圓棒參數(shù)及來流條件見表1。葉柵試驗(yàn)臺(tái)使用運(yùn)動(dòng)圓棒模擬從上游葉片尾緣脫落的尾跡。圓棒安裝在由交流電機(jī)驅(qū)動(dòng)的輻條尾跡發(fā)生器上,并使用光電轉(zhuǎn)速傳感器測量輻條轉(zhuǎn)速。從圖1(a)中可見,氣流出口平面到葉片尾緣的距離為Cx,圓棒到葉片前緣的距離為0.6Cx,圓棒進(jìn)口平面到圓棒的距離為0.6Cx,柵距方向上的兩邊界間的距離為1個(gè)柵距P,圓棒間距為P。使用DSA 3217壓力傳感器測量吸力面壁面靜壓,表面壓力信號(hào)使用DSA 3217壓力傳感器進(jìn)行測量。DSA 3217壓力傳感器包含16個(gè)溫度補(bǔ)償壓阻式壓力傳感器和1個(gè)氣動(dòng)校準(zhǔn)閥、內(nèi)存、16位A/D轉(zhuǎn)換器,以及1個(gè)微型處理器。微型處理器用于補(bǔ)償溫度變化,控制校準(zhǔn)閥已執(zhí)行在線零位和多點(diǎn)校準(zhǔn),并能夠消除傳感器熱誤差,保證系統(tǒng)精度為0.05%。使用5孔探針和熱線分別測量出口動(dòng)壓與近壁面速度。

    圖1 帶上游圓棒的葉柵通道模型及葉柵測試試驗(yàn)臺(tái)

    表1 葉柵和圓棒參數(shù)及來流條件

    本文選取較低的雷諾數(shù)Re=50000和較低的自由來流湍流度Fsti=1.5%作為來流條件,根據(jù)真實(shí)LPT工況選擇流量系數(shù)Ф=1,折合頻率Fr=1.128。

    流量系數(shù)定義為

    式中:Ux為軸向速度;Ub為尾跡發(fā)生器移動(dòng)速度。

    折合頻率定義為

    式中:Cx為軸向弦長;P為柵距。

    1.2 數(shù)值模擬

    本文使用商用CFD軟件對(duì)流動(dòng)進(jìn)行仿真,采用LES湍流模型耦合Smagorinsky亞格子模型進(jìn)行數(shù)值模擬,并通過試驗(yàn)校核??臻g離散項(xiàng)設(shè)置為中心差分,時(shí)間離散項(xiàng)設(shè)置為2階精度向后歐拉積分。同時(shí)為彌補(bǔ)Smagorinsky亞格子模型帶來的近壁面渦黏度過大的問題,采用Van-Driest壁面函數(shù)[21]進(jìn)行修正,并將Smagorinsky系數(shù)設(shè)為0.1。Funazaki[22]指出,對(duì)于高質(zhì)量、網(wǎng)格點(diǎn)數(shù)多的計(jì)算網(wǎng)格,基于LES模型耦合Smagorinsky亞格子模型的數(shù)值仿真能夠準(zhǔn)確地模擬吸力面附面層特性。

    Packb網(wǎng)格拓?fù)浜陀?jì)算域如圖2所示。Wu等[23]的研究表明計(jì)算域展向?qū)挾炔粦?yīng)影響展向或流向的渦結(jié)構(gòu)發(fā)展。為保證能夠捕捉詳細(xì)的瞬態(tài)結(jié)構(gòu)進(jìn)行分析,本文仿真的計(jì)算域展向?qū)挾榷?.18Cx。環(huán)繞葉片近壁面區(qū)域1周為“O”型拓?fù)洌渫鈬鷧^(qū)域?yàn)椤癏”型拓?fù)?。?jì)算域網(wǎng)格在x方向上有560個(gè)節(jié)點(diǎn),y方向上有225個(gè)節(jié)點(diǎn),z方向上有51個(gè)節(jié)點(diǎn),x-y平面上環(huán)繞葉片近壁面1周的網(wǎng)格節(jié)點(diǎn)數(shù)為770。同時(shí)葉片壁面法向上第1層網(wǎng)格最大高度設(shè)置為0.003 mm,約為柵距的4.5×10-5倍,以保證在該雷諾數(shù)下y+<1。為捕捉附面層的發(fā)展,在葉片壁面法向上12%的柵距長度內(nèi),設(shè)置了100個(gè)節(jié)點(diǎn),并保證x+、y+小于40。在葉片計(jì)算域內(nèi)約有700萬個(gè)節(jié)點(diǎn),在尾跡發(fā)生器計(jì)算域內(nèi)則有約80萬個(gè)節(jié)點(diǎn),整個(gè)計(jì)算域約有800萬節(jié)點(diǎn)。

    圖2 Packb網(wǎng)格拓?fù)浜陀?jì)算域

    時(shí)間步長約為Δt=0.0256 ms,Courant數(shù)(UinΔt/Δxmin)約為0.176(Uin為進(jìn)口速度,Δxmin為網(wǎng)格在x方向的最小長度)。每個(gè)尾跡掃掠周期T=2000Δt。為保證數(shù)值模擬的穩(wěn)定性,采用收斂的雷諾平均(Reynolds Average Navier Stokes,RANS)結(jié)果作為初場,再采用SST-SAS模型計(jì)算1個(gè)尾跡掃掠周期,最后采用LES模型耦合Smagorinsky亞格子模型計(jì)算10個(gè)尾跡掃掠周期。統(tǒng)計(jì)時(shí)均結(jié)果時(shí),排除第1周期,對(duì)剩余的9個(gè)周期每隔3個(gè)周期進(jìn)行鎖相平均。使用1臺(tái)10核20線程的工作站進(jìn)行計(jì)算,每個(gè)工況耗時(shí)25天。

    1.3 試驗(yàn)驗(yàn)證

    本文吸力面時(shí)均靜壓系數(shù)Cp的定義為

    式中:Pt,in為進(jìn)口總壓;Ps為當(dāng)?shù)仂o壓;ρ為密度;Uout

    吸力面時(shí)均靜壓系數(shù)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比如圖3所示,其中吸力面流向位置s通過吸力面弧長S0進(jìn)行無量綱化。

    圖3 吸力面時(shí)均靜壓系數(shù)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比

    從圖中可見,數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果在速度峰值點(diǎn)、“平臺(tái)”起始點(diǎn)、“平臺(tái)”結(jié)束點(diǎn)和壓力恢復(fù)點(diǎn)等位置基本吻合。同時(shí)用黑點(diǎn)和數(shù)字分別將數(shù)值模擬結(jié)果的時(shí)均分離起始點(diǎn)、時(shí)均轉(zhuǎn)捩點(diǎn)和時(shí)均再附點(diǎn)標(biāo)出。一般認(rèn)為,在形狀因子H12>3.5時(shí)附面層發(fā)生分離[24],在間歇因子>0.25時(shí)開始發(fā)生轉(zhuǎn)捩。通過計(jì)算各時(shí)刻的形狀因子和間歇因子并進(jìn)行時(shí)均處理,得到時(shí)均分離起始點(diǎn)1位于60%S0,時(shí)均轉(zhuǎn)捩點(diǎn)2位于70%S0,時(shí)均再附點(diǎn)3位于78%S0。其中間歇因子采用Solomon等[25]給出的方法計(jì)算。

    Mayle[2]認(rèn)為LPT中的轉(zhuǎn)捩類型主要有旁路轉(zhuǎn)捩和分離流轉(zhuǎn)捩,分離流轉(zhuǎn)捩可能觀察不到T-S波的存在,但KH不穩(wěn)定性在這2種轉(zhuǎn)捩類型中均起到了重要的作用。T-S波通常出現(xiàn)在轉(zhuǎn)捩過程的初期,是最先出現(xiàn)在層流流動(dòng)中的不穩(wěn)定性表現(xiàn)形式,Walker[26]建立了T-S波頻率fTS和附面層參數(shù)間的經(jīng)驗(yàn)公式,并將fTS定義為

    式中:Ue為附面層邊緣速度;δ*為T-S波第1次出現(xiàn)在層流附面層內(nèi)時(shí)對(duì)應(yīng)位置的位移厚度;ν為動(dòng)力黏度。

    經(jīng)計(jì)算,在Re=50000工況下,fTS≈220 Hz。根據(jù)Hughes等[27]的研究結(jié)果,式(4)計(jì)算的fTS為最大的放大率,在試驗(yàn)中還能觀測到更低的T-S波頻率。Graveline等[28]計(jì)算的T-S波頻率與試驗(yàn)觀測也存在一定誤差,但仍處于同一量級(jí)。

    Volino[29]建立了KH不穩(wěn)定性主導(dǎo)頻率的經(jīng)驗(yàn)公式,并將其定義為

    式中:δs為自由剪切層厚度;k為系數(shù),Volino[29]認(rèn)為k≈1.0。

    經(jīng)計(jì)算,本文fKH≈2800 Hz。Graveline等[28]認(rèn)為頻譜圖中的峰值與渦掃過熱線探頭有關(guān),當(dāng)渦結(jié)構(gòu)通過熱線時(shí)會(huì)在頻譜圖上留下1個(gè)頻率范圍較寬的峰值,并認(rèn)為該峰值是由KH不穩(wěn)定性造成的。Tennekes等[30]指出渦的頻率約為0.6f~1.6f,其中f與渦通過的頻率有關(guān)。因此表征KH不穩(wěn)定性的fKH不僅隨流向位置的變化而變化,反映到頻譜圖中還應(yīng)當(dāng)是1個(gè)范圍值。

    此外,已有不少研究在理想試驗(yàn)條件[31-32]和CFD結(jié)果[33]中觀測到KH渦在向下游運(yùn)動(dòng)的過程中會(huì)發(fā)生渦配對(duì)(vortices pair-up)現(xiàn)象,意味著在渦向下游運(yùn)動(dòng)的同時(shí)fKH會(huì)減半。

    本文還使用熱線對(duì)吸力面表面的速度進(jìn)行測量,并對(duì)無量綱流向速度擾動(dòng)量的平方進(jìn)行頻譜分析得到能譜函數(shù)E1(f),其中為流向擾動(dòng)速度,Umain為主流速度,f為頻率。吸力面壁面速度頻譜試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比如圖4所示。對(duì)吸力面壁面50%S0(層流)、70%S0(轉(zhuǎn)捩過程中)、80%S0(完全湍流)位置處的頻譜進(jìn)行對(duì)比,并用黑色圓圈圈出尾跡掃掠擾動(dòng)(wake passing)、fTS和fKH。

    圖4 吸力面壁面速度頻譜試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比

    從圖4(a)中可見,約在125 Hz時(shí)出現(xiàn)的第1個(gè)明顯的波峰是由尾跡掃掠造成的,其對(duì)應(yīng)的頻率與尾跡掃掠的頻率基本符合;在170~190 Hz時(shí)出現(xiàn)了第2個(gè)明顯的波峰,這與本文計(jì)算的fTS≈220 Hz處于同一量級(jí),因此可以認(rèn)為第2個(gè)波峰是由T-S波造成的。在2400~3600 Hz內(nèi)可觀測到較為明顯的寬峰,根據(jù)Graveline等[28]和Tennekes等[30]的研究可知該寬峰是由KH不穩(wěn)定性引起的,同時(shí)在80%S0頻譜圖1000~1500 Hz內(nèi)可觀測到另1個(gè)寬峰,這是由KH渦向下游運(yùn)動(dòng)時(shí)發(fā)生的渦配對(duì)現(xiàn)象引起的頻率減半。此外,在80%S0時(shí)處于完全湍流狀態(tài)的頻譜線在103<f<104內(nèi)的衰減規(guī)律符合-5/3原則。

    對(duì)比圖4(a)、(b)可知,試驗(yàn)和數(shù)值模擬得到的尾跡掃掠、fTS、fKH的頻率基本吻合,在80%S0時(shí)頻譜在103<f<104內(nèi)衰減規(guī)律均符合-5/3原則,因此認(rèn)為數(shù)值仿真能夠準(zhǔn)確地模擬吸力面附面層流動(dòng)特性。

    2 結(jié)果分析

    主要對(duì)上游尾跡與吸力面附面層相互作用以及多種類型轉(zhuǎn)捩的瞬態(tài)流動(dòng)特性進(jìn)行了分析。在1個(gè)尾跡掃掠周期內(nèi)葉柵通道的瞬態(tài)流動(dòng)細(xì)節(jié)如圖5所示。圖中葉柵通道截面使用湍動(dòng)能(Turbulence Kinetic Energy,TKE)云圖標(biāo)注尾跡與附面層展向卷起渦,并用紅線標(biāo)出了近壁面速度流線。為觀測尾跡與附面層的作用以及附面層內(nèi)部的瞬態(tài)流動(dòng)結(jié)構(gòu),使用距離葉片吸力面的半附面層高度作為觀測面,在20張瞬態(tài)流動(dòng)圖中給出了觀測面上的湍動(dòng)能TKE、無量綱壁面法向渦量wy*和無量綱壁面法向速度uy*。根據(jù)Coull[8]的研究可知,一定強(qiáng)度的wy*和uy

    *能夠反映出附面層內(nèi)Klebanoff條紋和KH結(jié)構(gòu)。其中wy*越大,速度虧損越大,Klebanoff條紋的強(qiáng)度也越強(qiáng)。根據(jù)Jeong等[34]的研究可知,可以使用速度梯度張量特征方程的第2個(gè)不變量Q來標(biāo)識(shí)渦結(jié)構(gòu),Q=0.5(||Ω2||-||S2||),其中S為剪切應(yīng)變率,Ω為渦量。當(dāng)Q>0時(shí),存在渦結(jié)構(gòu)。因此在圖5(a)左下角使用Q等值面(Q=1.5×107)標(biāo)出了65%S0~100%S0的近壁面渦結(jié)構(gòu)和半附面層高度觀測面的wy*云圖。在圖5中還使用紅色實(shí)線和字母標(biāo)記出了尾跡前緣線L、尾跡中心線A、尾跡尾緣線T、尾跡誘導(dǎo)轉(zhuǎn)捩區(qū)尾緣線B和抑制區(qū)尾緣線C的位置。這些位置線通過計(jì)算各流向位置附面層邊緣速度擾動(dòng)在時(shí)間上的分布,并根據(jù)速度擾動(dòng)波峰波谷來確定。尾跡與附面層作用將誘導(dǎo)附面層發(fā)生轉(zhuǎn)捩形成湍斑,在向下游運(yùn)動(dòng)和生長的過程中,湍斑后部會(huì)形成1個(gè)抑制區(qū),湍斑尾緣和抑制區(qū)尾緣分別以50%和33%的主流速度向下游運(yùn)動(dòng)。一般認(rèn)為當(dāng)間歇因子大于0.25時(shí),轉(zhuǎn)捩開始發(fā)生。本文采用Solomon等[25]給出的方法計(jì)算間歇因子。

    圖5 在1個(gè)尾跡掃掠周期內(nèi)葉柵通道的瞬態(tài)流動(dòng)

    從圖5(a)中可見,在0.05T時(shí),上一條尾跡已經(jīng)離開葉柵通道,其誘導(dǎo)產(chǎn)生的轉(zhuǎn)捩區(qū)尾緣線B和抑制區(qū)尾緣線C尚未離開觀測面,表明上一條尾跡產(chǎn)生的尾跡誘導(dǎo)湍流區(qū)和抑制區(qū)在葉片尾緣與附面層發(fā)生相互作用;在70%S0~75%S0附近有強(qiáng)度較弱的流向結(jié)構(gòu)(黑色橢圓圈出),從左下角的渦結(jié)構(gòu)示意圖中可見相似的流向結(jié)構(gòu)(黑色橢圓圈出),二者可以對(duì)應(yīng),表明上一條尾跡產(chǎn)生的尾跡放大Klebanoff條紋仍與葉片尾緣附近的附面層發(fā)生相互作用;同時(shí),在該時(shí)刻下一條尾跡W1尚未與吸力面接觸,處于尾跡通過的間歇期,因此無法在葉片前部觀測到明顯的流向條紋結(jié)構(gòu);在C線上游、L線下游的65%S0~75%S0區(qū)域內(nèi)有渦卷起現(xiàn)象(箭頭標(biāo)出),以及渦卷起造成的觀測面的湍動(dòng)能增大。此時(shí)處于尾跡掃掠的間歇期,同時(shí)上一條尾跡產(chǎn)生的抑制區(qū)尾緣線C正在離開時(shí)均分離區(qū)域,表明吸力面附面層轉(zhuǎn)捩模式正在由尾跡誘導(dǎo)轉(zhuǎn)捩向自然轉(zhuǎn)捩轉(zhuǎn)變。

    從圖5(b)中可見,在0.10T時(shí),尾跡W1尚未與吸力面接觸,依然無法在葉片中部觀測到明顯的Klebanoff條紋,從流線和左下角的渦結(jié)構(gòu)示意圖中可見渦卷起并向下游運(yùn)動(dòng),造成觀測面上湍動(dòng)能TKE增大。此時(shí)處于尾跡掃掠的間歇期,同時(shí)抑制區(qū)尾緣線C已離開時(shí)均分離區(qū)域,吸力面附面層的轉(zhuǎn)捩模式為自然轉(zhuǎn)捩。

    從圖5(c)、(d)中可見,在0.15T~0.20T時(shí),隨著尾跡向下游運(yùn)動(dòng),在尾跡中心線A附近觀測到流向條紋結(jié)構(gòu)(黑色橢圓圈出),這是由于尾跡擾動(dòng)在滲透入吸力面前緣的附面層時(shí),因“剪切遮蔽”效應(yīng)形成尾跡放大Klebanoff條紋,尾跡放大條紋區(qū)前緣和尾緣分別以88%和50%的主流速度向下游運(yùn)動(dòng),其中最強(qiáng)條紋以70%的主流速度向下游運(yùn)動(dòng),因此條紋將始終落后于尾跡前緣線L。

    從圖5(e)中可見,在0.25T時(shí),尾跡W1開始與吸力面附面層接觸并發(fā)生相互作用,尾跡前緣線L抵達(dá)70%S0,并進(jìn)入分離剪切層中,尾跡擾動(dòng)開始滲透進(jìn)入分離剪切層并觸發(fā)KH不穩(wěn)定性,展向的KH渦不斷地沿分離泡邊緣形成、脫落并向下游運(yùn)動(dòng);從近壁面流線和左下角的渦結(jié)構(gòu)示意圖中可見65%S0~70%S0區(qū)域內(nèi)有全展向KH結(jié)構(gòu)(Full Span K-H Structure)(箭頭標(biāo)出,標(biāo)記為1),該KH結(jié)構(gòu)以37%的主流速度向下游運(yùn)動(dòng)并抬起,此時(shí)無法在觀測面uy*云圖上觀測到展向結(jié)構(gòu),這是由于該展向結(jié)構(gòu)尚未抵達(dá)半附面層高度。

    從圖5(f)中可見,在0.30T時(shí),標(biāo)記為1的全展向KH結(jié)構(gòu)在向下游運(yùn)動(dòng)的同時(shí)不斷抬起,并抵達(dá)半附面層高度,因此在觀測面uy*云圖72%S0附近可觀測到類似的展向結(jié)構(gòu)(黑色橢圓標(biāo)出)。

    從圖5(g)中可見,在0.35T時(shí),左下角渦結(jié)構(gòu)示意圖標(biāo)記為2的全展向KH渦緊隨全展向KH渦1于71%S0附近形成,并向下游運(yùn)動(dòng),在uy*云圖的對(duì)應(yīng)位置可見完整的全展向KH結(jié)構(gòu);同時(shí),從觀測面的wy*云圖中可見,尾跡放大Klebanoff條紋的發(fā)展呈現(xiàn)出楔形的展向不對(duì)稱性,這一楔形分布與觀測面的湍動(dòng)能分布形態(tài)吻合,與Nagabhushana[19]所觀測到的條紋展向不對(duì)稱性相似;從觀測面的wy*云圖中還可見,尾跡放大Klebanoff條紋最前緣抵達(dá)65%S0,尚未與位于71%S0的KH結(jié)構(gòu)接觸,因此可以觀測到較為完整的全展向KH結(jié)構(gòu)1、2。

    從圖5(h)中可見,在0.40T時(shí),尾跡放大Klebanoff條紋與全展向KH結(jié)構(gòu)接觸,同時(shí)從觀測面的uy*云圖和渦結(jié)構(gòu)示意圖中可見全展向KH結(jié)構(gòu)1、2開始發(fā)生扭曲,表明KH結(jié)構(gòu)在尾跡放大Klebanoff條紋的作用下開始崩潰。

    從圖5(i)中可見,在0.45T時(shí),楔形分布的尾跡放大Klebanoff條紋與全展向KH結(jié)構(gòu)1、2接觸,Klebanoff條紋會(huì)先沖擊全展向KH結(jié)構(gòu)的一側(cè),并使全展向KH結(jié)構(gòu)的一側(cè)發(fā)生劇烈扭曲(圓圈標(biāo)出)。

    從圖5(j)、(k)中可見,在0.50T~0.55T時(shí),已觀測不到全展向KH結(jié)構(gòu),但能在觀測面uy*云圖的73%S0~82%S0區(qū)域內(nèi)觀測到局部KH結(jié)構(gòu)(Part Span K-H Structure)(黑色方框標(biāo)出),這是由于尾跡放大Klebanoff條紋沖擊全展向KH結(jié)構(gòu),使全展向KH結(jié)構(gòu)破碎為局部KH結(jié)構(gòu),這一現(xiàn)象可在觀測面的wy*云圖和uy*云圖中觀測到;此外,從近壁面流線和渦結(jié)構(gòu)示意圖中已觀測不到新的卷起渦生成,表明在尾跡的作用下分離泡被完全抑制。

    從圖5(l)中可見,在0.60T時(shí),隨著尾跡放大,Klebanoff條紋不斷向下游運(yùn)動(dòng)并破壞KH結(jié)構(gòu),從渦結(jié)構(gòu)示意圖和uy*云圖中已看不到明顯的局部KH結(jié)構(gòu);同時(shí),從wy*云圖中可見尾跡放大Klebanoff條紋強(qiáng)度的衰減,這是由于尾跡不斷向下游運(yùn)動(dòng),最終不再與葉片前緣發(fā)生作用造成的,而已有條紋繼續(xù)向下游運(yùn)動(dòng)、拉伸,強(qiáng)度不斷衰減。

    從圖5(m)中可見,在0.65T時(shí),受尾跡放大Klebanoff條紋影響,已無法觀測到明顯的KH結(jié)構(gòu);同時(shí),尾跡誘導(dǎo)吸力面附面層發(fā)生尾跡誘導(dǎo)轉(zhuǎn)捩,在T線與B線之間形成尾跡誘導(dǎo)湍流區(qū),在B線和C線之間則形成抑制區(qū),湍流區(qū)和抑制區(qū)能夠有效地抑制分離,因此無法通過近壁面流線觀測到渦卷起現(xiàn)象。

    從圖5(n)中可見,在0.70T時(shí),隨著尾跡進(jìn)一步向下游運(yùn)動(dòng),尾跡尾緣線T抵達(dá)78%S0,通過近壁面流線觀測到1個(gè)較小的渦在70%S0附近卷起,分離泡開始重新生成,因此出現(xiàn)較弱的由于分離的形成而造成的渦卷起現(xiàn)象;此外,在渦結(jié)構(gòu)示意圖中的C線下游約75%S0處觀測到了Λ渦3的形成,并用箭頭標(biāo)出了Λ渦的2條“腿”,根據(jù)Durbin等[35]的研究可知,Λ渦是T-S波和Klebanoff條紋共同作用的結(jié)果,由于抑制區(qū)內(nèi)的T-S波擾動(dòng)較小,因此無法在B線和C線之間的抑制區(qū)內(nèi)觀測到明顯的Λ渦結(jié)構(gòu)。

    從圖5(n)、(o)左下角的渦結(jié)構(gòu)示意圖中可見,在0.70T~0.75T時(shí),Λ渦3拉伸并最終破裂為湍流。

    從圖5(p)中可見,在0.80T時(shí),新的Λ渦4在75%S0附近生成;由于尾跡放大Klebanoff條紋在葉片中后部強(qiáng)度衰減,渦的強(qiáng)度變?nèi)?,?條“腿”的不對(duì)稱性更為明顯;同時(shí),隨著尾跡向下游運(yùn)動(dòng),抑制區(qū)前緣線B不斷往下游運(yùn)動(dòng),Λ渦的強(qiáng)度將會(huì)變得更弱甚至被抑制。

    從圖5(q)、(r)中可見,在0.85T~0.90T時(shí),尾跡放大Klebanoff條紋繼續(xù)向下游運(yùn)動(dòng),其強(qiáng)度不斷衰減;Λ渦4在向下游運(yùn)動(dòng)的過程中強(qiáng)度不斷減弱,2條“腿”的不對(duì)稱性更為明顯;同時(shí),抑制區(qū)前緣線B和抑制區(qū)尾緣線C分別以50%和33%的主流速度向下游運(yùn)動(dòng),隨著時(shí)間的推移,抑制區(qū)在向下游運(yùn)動(dòng)的同時(shí),其所占面積也會(huì)增大,一方面抑制區(qū)會(huì)抑制渦,已無法在渦結(jié)構(gòu)示意圖中觀測到新Λ渦的生成,另一方面分離泡會(huì)在抑制區(qū)后重新生成,并能通過抑制區(qū)尾緣線C上游的近壁面流線觀測到渦卷起現(xiàn)象不斷增強(qiáng),分離泡尺寸也在不斷變大;此時(shí)尾跡W1即將離開葉柵通道,附面層轉(zhuǎn)捩模式由尾跡誘導(dǎo)轉(zhuǎn)捩向自然轉(zhuǎn)捩過渡。

    從圖5(r)、(s)、(t)中可見,在0.90T~1.00T時(shí),隨著尾跡向下游運(yùn)動(dòng)并最終離開葉柵通道,尾跡放大Klebanoff條紋在向下游運(yùn)動(dòng)的過程中,其在葉片中后部的強(qiáng)度不斷衰減;受抑制區(qū)影響,Λ渦在向下游運(yùn)動(dòng)的過程中強(qiáng)度也不斷減弱;同時(shí),下一條尾跡尚未與吸力面附面層接觸,既不會(huì)在葉片前部形成尾跡放大Klebanoff條紋,也不會(huì)觸發(fā)分離剪切層的全展向KH卷起渦,無法從wy*云圖中觀測到明顯的流向結(jié)構(gòu),但能通過流線觀測到分離泡以及不斷增強(qiáng)的渦卷起現(xiàn)象,并導(dǎo)致觀測面TKE增大。

    3 結(jié)論

    (1)在2條尾跡掃掠之間,上一條尾跡剛剛離開葉片尾緣,其尾跡誘導(dǎo)的湍流區(qū)和抑制區(qū)仍在影響葉片尾緣附面層的發(fā)展,同時(shí)下一條尾跡尚未與吸力面接觸,此時(shí)附面層轉(zhuǎn)捩模式由尾跡誘導(dǎo)轉(zhuǎn)捩向自然轉(zhuǎn)捩轉(zhuǎn)變,分離泡在抑制區(qū)上游重新生成。

    (2)下一條尾跡在葉片前緣與吸力面接觸,產(chǎn)生尾跡放大Klebanoff條紋,其前緣以88%的主流速度向下游運(yùn)動(dòng),隨后尾跡接觸附面層并與分離泡相互作用,觸發(fā)剪切層的KH不穩(wěn)定性形成全展向KH卷起渦,并抑制分離泡。尾跡放大Klebanoff條紋將追趕并沖擊全展向KH卷起渦,使其崩潰為局部KH渦并最終破裂形成Λ渦和全湍流。

    猜你喜歡
    附面層葉柵尾跡
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)
    亞聲速壓氣機(jī)平面葉柵及其改型的吹風(fēng)試驗(yàn)
    基于數(shù)值模擬的流場附面層邊緣識(shí)別方法
    基于EEMD-Hilbert譜的渦街流量計(jì)尾跡振蕩特性
    超聲壓氣機(jī)葉柵流場的數(shù)值模擬與試驗(yàn)驗(yàn)證
    超音速透平葉柵造型設(shè)計(jì)及分析
    針對(duì)軸流壓氣機(jī)的非軸對(duì)稱端壁造型優(yōu)化設(shè)計(jì)
    基于FABEMD和Goldstein濾波器的SAR艦船尾跡圖像增強(qiáng)方法
    發(fā)動(dòng)機(jī)進(jìn)口附面層測量試驗(yàn)與數(shù)值模擬
    成人亚洲精品一区在线观看| 亚洲国产看品久久| 国产老妇伦熟女老妇高清| e午夜精品久久久久久久| 亚洲九九香蕉| 天天影视国产精品| 九色亚洲精品在线播放| 精品第一国产精品| 国产精品国产三级专区第一集| 99国产精品免费福利视频| 免费看av在线观看网站| 啦啦啦视频在线资源免费观看| 极品人妻少妇av视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看一区二区三区激情| 水蜜桃什么品种好| 午夜免费鲁丝| 日本wwww免费看| 亚洲人成77777在线视频| 国产一区有黄有色的免费视频| 日韩一区二区三区影片| 天堂俺去俺来也www色官网| 黄色视频在线播放观看不卡| 国精品久久久久久国模美| 日韩av在线免费看完整版不卡| 天天操日日干夜夜撸| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 亚洲av美国av| 亚洲精品成人av观看孕妇| 久久久久久久精品精品| 亚洲人成电影免费在线| 欧美激情高清一区二区三区| 赤兔流量卡办理| 老司机在亚洲福利影院| 人人妻人人添人人爽欧美一区卜| 亚洲av成人精品一二三区| 日本a在线网址| 国产伦理片在线播放av一区| 成年av动漫网址| 亚洲一区中文字幕在线| 国产精品一区二区免费欧美 | 亚洲第一青青草原| 91字幕亚洲| 国产黄色视频一区二区在线观看| 看免费av毛片| 色综合欧美亚洲国产小说| 亚洲熟女精品中文字幕| 亚洲九九香蕉| 熟女av电影| 97在线人人人人妻| 亚洲欧美日韩另类电影网站| 高清不卡的av网站| 人人妻人人添人人爽欧美一区卜| 国产男女超爽视频在线观看| 美女高潮到喷水免费观看| 久久人妻福利社区极品人妻图片 | 男男h啪啪无遮挡| 中文字幕精品免费在线观看视频| 男女边摸边吃奶| 亚洲欧美精品综合一区二区三区| av国产精品久久久久影院| 狂野欧美激情性xxxx| 国产色视频综合| 51午夜福利影视在线观看| 尾随美女入室| 无遮挡黄片免费观看| 热re99久久精品国产66热6| 美国免费a级毛片| 欧美日韩视频高清一区二区三区二| 女人精品久久久久毛片| 国产精品国产三级专区第一集| 国产av精品麻豆| 久久av网站| 亚洲国产毛片av蜜桃av| 国产无遮挡羞羞视频在线观看| 在现免费观看毛片| 亚洲精品av麻豆狂野| 男女边吃奶边做爰视频| 亚洲国产av新网站| 亚洲五月婷婷丁香| 满18在线观看网站| 五月开心婷婷网| 亚洲一码二码三码区别大吗| 男女下面插进去视频免费观看| 肉色欧美久久久久久久蜜桃| 日韩熟女老妇一区二区性免费视频| 午夜福利,免费看| 精品免费久久久久久久清纯 | 在线观看免费午夜福利视频| 成人国语在线视频| 欧美日韩成人在线一区二区| 国产亚洲精品久久久久5区| 欧美黑人欧美精品刺激| 久久久国产精品麻豆| 黄片小视频在线播放| 国产精品麻豆人妻色哟哟久久| 老汉色∧v一级毛片| 黄色视频在线播放观看不卡| 欧美乱码精品一区二区三区| 久久人人爽av亚洲精品天堂| 日韩av免费高清视频| 亚洲国产av影院在线观看| 国产精品一国产av| 亚洲精品自拍成人| 黄频高清免费视频| 成人免费观看视频高清| kizo精华| h视频一区二区三区| 丁香六月天网| 久久久久久人人人人人| 亚洲成人免费电影在线观看 | 国产精品一区二区在线观看99| 久久亚洲国产成人精品v| 国产淫语在线视频| 久久久国产欧美日韩av| 久久久久久人人人人人| 天堂俺去俺来也www色官网| 老司机在亚洲福利影院| 免费在线观看视频国产中文字幕亚洲 | 高清视频免费观看一区二区| 午夜免费观看性视频| 欧美精品高潮呻吟av久久| 国产在线一区二区三区精| 老汉色∧v一级毛片| 日韩中文字幕欧美一区二区 | 国产有黄有色有爽视频| 女性生殖器流出的白浆| 熟女av电影| 1024香蕉在线观看| 精品一品国产午夜福利视频| 天天影视国产精品| 午夜福利视频在线观看免费| 亚洲国产精品国产精品| 三上悠亚av全集在线观看| 啦啦啦啦在线视频资源| 欧美日韩精品网址| 成年人午夜在线观看视频| 亚洲精品国产av成人精品| 好男人视频免费观看在线| 久久精品国产综合久久久| 美国免费a级毛片| 精品久久久精品久久久| 女人爽到高潮嗷嗷叫在线视频| 男人爽女人下面视频在线观看| 亚洲情色 制服丝袜| 男女高潮啪啪啪动态图| 在线av久久热| 国产97色在线日韩免费| √禁漫天堂资源中文www| 一级毛片 在线播放| 亚洲国产av影院在线观看| 可以免费在线观看a视频的电影网站| 亚洲综合色网址| 久久久国产一区二区| 欧美日韩福利视频一区二区| 亚洲精品自拍成人| 欧美在线一区亚洲| 亚洲国产精品一区二区三区在线| 免费观看人在逋| 中文字幕另类日韩欧美亚洲嫩草| 久久国产精品大桥未久av| 国产在线免费精品| 大香蕉久久网| 久久鲁丝午夜福利片| 日本黄色日本黄色录像| 欧美黑人欧美精品刺激| 99re6热这里在线精品视频| 欧美变态另类bdsm刘玥| 韩国精品一区二区三区| 国产片特级美女逼逼视频| 青春草亚洲视频在线观看| 日韩中文字幕视频在线看片| 各种免费的搞黄视频| 少妇的丰满在线观看| 可以免费在线观看a视频的电影网站| 亚洲精品乱久久久久久| 少妇的丰满在线观看| 大码成人一级视频| av又黄又爽大尺度在线免费看| 各种免费的搞黄视频| 国产男人的电影天堂91| 日韩av免费高清视频| svipshipincom国产片| 亚洲欧洲日产国产| 国产精品 欧美亚洲| 国产视频首页在线观看| 国产高清国产精品国产三级| 成人影院久久| 国产精品一区二区免费欧美 | 免费久久久久久久精品成人欧美视频| 啦啦啦在线观看免费高清www| 国产亚洲欧美精品永久| 熟女少妇亚洲综合色aaa.| 丰满人妻熟妇乱又伦精品不卡| 国产高清videossex| 国产一区二区 视频在线| 黄色a级毛片大全视频| 久久 成人 亚洲| 91九色精品人成在线观看| 国产在视频线精品| 亚洲午夜精品一区,二区,三区| 少妇 在线观看| 少妇 在线观看| 妹子高潮喷水视频| 天天躁夜夜躁狠狠躁躁| 精品人妻一区二区三区麻豆| 一区福利在线观看| 久久亚洲国产成人精品v| 看免费成人av毛片| 国产精品一区二区在线观看99| 国产高清不卡午夜福利| 国产免费福利视频在线观看| 老鸭窝网址在线观看| 视频在线观看一区二区三区| 香蕉国产在线看| bbb黄色大片| 国产精品人妻久久久影院| svipshipincom国产片| 男女国产视频网站| xxxhd国产人妻xxx| 天天躁日日躁夜夜躁夜夜| 曰老女人黄片| 秋霞在线观看毛片| 国产欧美日韩一区二区三区在线| 亚洲少妇的诱惑av| 精品视频人人做人人爽| 国产又色又爽无遮挡免| 国产伦理片在线播放av一区| 性高湖久久久久久久久免费观看| 亚洲人成77777在线视频| 免费观看av网站的网址| 99国产精品一区二区蜜桃av | 亚洲av欧美aⅴ国产| 久久精品亚洲av国产电影网| 成人黄色视频免费在线看| 黄色片一级片一级黄色片| 18禁裸乳无遮挡动漫免费视频| √禁漫天堂资源中文www| 一级毛片女人18水好多 | kizo精华| 国产精品欧美亚洲77777| 色视频在线一区二区三区| 9色porny在线观看| 好男人视频免费观看在线| 男人爽女人下面视频在线观看| 激情视频va一区二区三区| 乱人伦中国视频| 男人添女人高潮全过程视频| 热99久久久久精品小说推荐| 91精品三级在线观看| 欧美日韩综合久久久久久| 国产成人精品久久久久久| 中文字幕人妻丝袜制服| 1024视频免费在线观看| 国产极品粉嫩免费观看在线| 亚洲成国产人片在线观看| 亚洲精品国产色婷婷电影| 午夜免费鲁丝| 夜夜骑夜夜射夜夜干| 亚洲熟女毛片儿| av线在线观看网站| 18禁国产床啪视频网站| 91精品三级在线观看| 国产免费福利视频在线观看| 国产精品一国产av| 免费久久久久久久精品成人欧美视频| 少妇猛男粗大的猛烈进出视频| 97精品久久久久久久久久精品| 男女午夜视频在线观看| 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说| 黑人欧美特级aaaaaa片| 一区二区三区激情视频| 国产日韩欧美在线精品| 最新的欧美精品一区二区| 亚洲av电影在线观看一区二区三区| 777久久人妻少妇嫩草av网站| a级片在线免费高清观看视频| 午夜av观看不卡| 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 一本一本久久a久久精品综合妖精| 欧美精品一区二区大全| 午夜久久久在线观看| 亚洲国产精品成人久久小说| 精品福利观看| 国产在线观看jvid| 成人免费观看视频高清| 在线看a的网站| 亚洲熟女精品中文字幕| 欧美人与善性xxx| 黑丝袜美女国产一区| 伊人亚洲综合成人网| 亚洲天堂av无毛| 夫妻午夜视频| 欧美乱码精品一区二区三区| 亚洲中文字幕日韩| 亚洲少妇的诱惑av| 久久天堂一区二区三区四区| 亚洲av日韩在线播放| 一区二区av电影网| 精品一区在线观看国产| 亚洲成人手机| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 99精国产麻豆久久婷婷| 亚洲av国产av综合av卡| 性少妇av在线| 高清不卡的av网站| 亚洲成人免费电影在线观看 | a级毛片在线看网站| 少妇人妻久久综合中文| 国产av国产精品国产| 国产老妇伦熟女老妇高清| 晚上一个人看的免费电影| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久男人| 婷婷色综合大香蕉| 国产男女内射视频| 男男h啪啪无遮挡| 999精品在线视频| 国产欧美日韩一区二区三 | 免费在线观看影片大全网站 | 人人妻人人澡人人爽人人夜夜| 亚洲激情五月婷婷啪啪| 午夜激情久久久久久久| 不卡av一区二区三区| 亚洲专区中文字幕在线| 一区二区日韩欧美中文字幕| 大话2 男鬼变身卡| 久久久久久亚洲精品国产蜜桃av| 脱女人内裤的视频| 精品福利观看| 国产在线观看jvid| 丝袜美腿诱惑在线| 又大又黄又爽视频免费| 国产人伦9x9x在线观看| 超碰成人久久| 新久久久久国产一级毛片| 美女高潮到喷水免费观看| 最新的欧美精品一区二区| 精品久久蜜臀av无| 一本大道久久a久久精品| 成年av动漫网址| 欧美日本中文国产一区发布| 欧美精品人与动牲交sv欧美| 国产97色在线日韩免费| 精品久久久精品久久久| 日韩制服丝袜自拍偷拍| 在线精品无人区一区二区三| 电影成人av| 日韩免费高清中文字幕av| 久久国产精品影院| 男女无遮挡免费网站观看| 欧美激情高清一区二区三区| 观看av在线不卡| 搡老乐熟女国产| 亚洲欧美精品综合一区二区三区| 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲综合一区二区三区_| 日本色播在线视频| 国产99久久九九免费精品| 悠悠久久av| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠久久av| 国产亚洲午夜精品一区二区久久| 色婷婷久久久亚洲欧美| tube8黄色片| 黄色一级大片看看| 中文字幕人妻熟女乱码| 脱女人内裤的视频| 欧美成人精品欧美一级黄| 热99久久久久精品小说推荐| 永久免费av网站大全| 中文精品一卡2卡3卡4更新| 一本大道久久a久久精品| 乱人伦中国视频| 日韩熟女老妇一区二区性免费视频| 国产成人av教育| av国产久精品久网站免费入址| 美女视频免费永久观看网站| 国产不卡av网站在线观看| cao死你这个sao货| 欧美另类一区| 一二三四在线观看免费中文在| 国产又爽黄色视频| 亚洲黑人精品在线| 欧美国产精品一级二级三级| 下体分泌物呈黄色| 巨乳人妻的诱惑在线观看| 亚洲黑人精品在线| 搡老乐熟女国产| 亚洲精品第二区| 久久毛片免费看一区二区三区| 亚洲五月婷婷丁香| 校园人妻丝袜中文字幕| 一级黄色大片毛片| 久久久精品免费免费高清| 一区二区三区精品91| 日本av手机在线免费观看| 国产精品国产三级国产专区5o| 久久久欧美国产精品| 如日韩欧美国产精品一区二区三区| 黑人欧美特级aaaaaa片| 激情视频va一区二区三区| 熟女av电影| 欧美日韩福利视频一区二区| 大片免费播放器 马上看| 18禁观看日本| 国产精品一区二区在线观看99| 爱豆传媒免费全集在线观看| 五月天丁香电影| 日韩欧美一区视频在线观看| 亚洲人成电影免费在线| 亚洲精品中文字幕在线视频| 黄频高清免费视频| 免费在线观看日本一区| 国产伦理片在线播放av一区| 欧美国产精品一级二级三级| 亚洲伊人色综图| 在线观看免费午夜福利视频| 一区二区三区激情视频| av网站免费在线观看视频| av不卡在线播放| 成人亚洲精品一区在线观看| 黑人欧美特级aaaaaa片| 97在线人人人人妻| 欧美亚洲 丝袜 人妻 在线| 欧美日韩av久久| 在线看a的网站| 亚洲精品国产av成人精品| 波野结衣二区三区在线| 久久天躁狠狠躁夜夜2o2o | 国产精品一区二区免费欧美 | 精品国产一区二区三区四区第35| 女警被强在线播放| 日韩 亚洲 欧美在线| 国产福利在线免费观看视频| 91九色精品人成在线观看| 叶爱在线成人免费视频播放| 人妻人人澡人人爽人人| 久久久国产精品麻豆| 真人做人爱边吃奶动态| 深夜精品福利| 国产一区二区 视频在线| 热99久久久久精品小说推荐| 成人国产av品久久久| 最新的欧美精品一区二区| 777久久人妻少妇嫩草av网站| 狠狠婷婷综合久久久久久88av| 丁香六月欧美| 啦啦啦 在线观看视频| 精品视频人人做人人爽| 日本vs欧美在线观看视频| 亚洲精品第二区| 色网站视频免费| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 亚洲成国产人片在线观看| 亚洲精品在线美女| 肉色欧美久久久久久久蜜桃| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| 免费高清在线观看视频在线观看| 美国免费a级毛片| 国产精品av久久久久免费| av国产精品久久久久影院| 捣出白浆h1v1| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 亚洲精品国产一区二区精华液| 国产免费又黄又爽又色| 日韩大片免费观看网站| 亚洲伊人色综图| 精品熟女少妇八av免费久了| 久久久久久人人人人人| 又大又爽又粗| 下体分泌物呈黄色| www.999成人在线观看| 午夜久久久在线观看| 亚洲国产精品国产精品| 国产精品三级大全| 黄网站色视频无遮挡免费观看| 欧美日韩亚洲综合一区二区三区_| 精品少妇内射三级| 久久av网站| 极品人妻少妇av视频| 久久精品久久久久久久性| 久久人妻熟女aⅴ| av国产精品久久久久影院| 人妻 亚洲 视频| 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 一级,二级,三级黄色视频| 啦啦啦在线免费观看视频4| 欧美老熟妇乱子伦牲交| 成年女人毛片免费观看观看9 | 久久99热这里只频精品6学生| av欧美777| 最黄视频免费看| 黄色视频在线播放观看不卡| 视频区欧美日本亚洲| 久久影院123| 亚洲视频免费观看视频| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三区在线| 最新的欧美精品一区二区| 欧美国产精品va在线观看不卡| 欧美日韩综合久久久久久| 丰满少妇做爰视频| 91精品三级在线观看| 中国国产av一级| 精品人妻在线不人妻| 自拍欧美九色日韩亚洲蝌蚪91| 18禁观看日本| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区蜜桃| 日韩av在线免费看完整版不卡| 国产免费一区二区三区四区乱码| 尾随美女入室| 一边亲一边摸免费视频| 国产在线一区二区三区精| 两性夫妻黄色片| a 毛片基地| 久久中文字幕一级| 日日夜夜操网爽| 亚洲精品一二三| 久久久国产一区二区| 亚洲少妇的诱惑av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品一卡2卡三卡4卡5卡 | 国产成人欧美| 亚洲精品日韩在线中文字幕| 国产精品久久久久成人av| 女性生殖器流出的白浆| 亚洲情色 制服丝袜| 在线观看免费视频网站a站| 99热网站在线观看| 午夜视频精品福利| 丝袜人妻中文字幕| 一边摸一边抽搐一进一出视频| 色网站视频免费| 日日摸夜夜添夜夜爱| 免费在线观看完整版高清| 国产成人啪精品午夜网站| 免费黄频网站在线观看国产| 日本一区二区免费在线视频| 啦啦啦啦在线视频资源| 日韩av不卡免费在线播放| 免费在线观看视频国产中文字幕亚洲 | 女人爽到高潮嗷嗷叫在线视频| 在线天堂中文资源库| 久久热在线av| 国产午夜精品一二区理论片| 国产男人的电影天堂91| 亚洲少妇的诱惑av| 少妇人妻 视频| 亚洲国产成人一精品久久久| 91国产中文字幕| 亚洲五月色婷婷综合| 亚洲一码二码三码区别大吗| av电影中文网址| 最近最新中文字幕大全免费视频 | 成年人免费黄色播放视频| 欧美另类一区| 少妇 在线观看| 国产男女超爽视频在线观看| 久久 成人 亚洲| 一二三四社区在线视频社区8| 精品第一国产精品| kizo精华| 一级毛片电影观看| 午夜免费观看性视频| 久久久国产一区二区| av在线老鸭窝| 视频区图区小说| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 性色av一级| 亚洲美女黄色视频免费看| 国产亚洲av高清不卡| 色精品久久人妻99蜜桃| 国产精品一区二区精品视频观看| 高清黄色对白视频在线免费看| av不卡在线播放| 免费观看av网站的网址| 夜夜骑夜夜射夜夜干| 人体艺术视频欧美日本| 欧美日韩国产mv在线观看视频| videosex国产| 国产成人精品久久二区二区免费| 国产精品九九99| av国产精品久久久久影院| 成年人免费黄色播放视频| 看免费av毛片| 啦啦啦 在线观看视频| 黄色毛片三级朝国网站| av有码第一页| 99国产精品99久久久久| 色94色欧美一区二区| 91精品伊人久久大香线蕉| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 久久久国产精品麻豆| 精品少妇内射三级| 大陆偷拍与自拍| 9191精品国产免费久久| 大陆偷拍与自拍| 无限看片的www在线观看| 香蕉丝袜av| 成年人午夜在线观看视频| 大型av网站在线播放| 精品福利观看| 又大又爽又粗|