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

    展向離散抽吸法控制邊界層轉(zhuǎn)捩實驗研究

    2014-07-10 13:15:24李小寶王海文馮玉龍
    實驗流體力學 2014年6期
    關(guān)鍵詞:邊界層條帶振幅

    郭 輝,李小寶,王海文,馮玉龍

    (1.北京航空航天大學航空科學與工程學院,北京100191;2.流體力學教育部重點實驗室(北京航空航天大學),北京100191)

    展向離散抽吸法控制邊界層轉(zhuǎn)捩實驗研究

    郭 輝,李小寶,王海文,馮玉龍

    (1.北京航空航天大學航空科學與工程學院,北京100191;2.流體力學教育部重點實驗室(北京航空航天大學),北京100191)

    通過實驗成功發(fā)展了一種生成大振幅穩(wěn)定條帶的有效方法,即展向離散抽吸方法。在此基礎(chǔ)上,開展條帶控制邊界層轉(zhuǎn)捩的研究。實驗在水洞中進行,以零壓力梯度平板邊界層為研究對象,利用氫氣泡時間線法觀測引入條帶前后人工激發(fā)轉(zhuǎn)捩邊界層中擾動發(fā)展變化,分析條帶對轉(zhuǎn)捩的控制效果及參數(shù)影響。結(jié)果表明,展向離散抽吸方法生成的穩(wěn)定條帶最大振幅可達28.4%U(U為自由流速度);實驗中引入的寬度14和28mm的穩(wěn)定條帶都能起到抑制轉(zhuǎn)捩的作用;條帶振幅越大、寬度越窄,抑制效果越明顯。研究結(jié)果為探索降低水/空氣中航行器摩阻的新技術(shù)提供有價值的參考。

    條帶;邊界層;轉(zhuǎn)捩;流動控制;氫氣泡法

    0 引 言

    有效降低阻力對于提高空氣/水中航行器性能和經(jīng)濟性以及節(jié)能環(huán)保都是非常重要的。巡航狀態(tài)的空氣/水中航行器大都是表現(xiàn)為小迎角的流線型外形,流動分離不顯著或僅限于后緣很小的區(qū)域,以摩擦阻力為主。如在巡航狀態(tài)的飛行器摩阻占到總阻力的50%,水中航行器的摩阻可高達80%。通常邊界層的發(fā)展始于層流,經(jīng)轉(zhuǎn)捩發(fā)展為湍流。典型飛行雷諾數(shù)條件下,邊界層內(nèi)層流的摩阻比湍流小一個量級,雷諾數(shù)越大差別越顯著。于是減小摩阻發(fā)展出2條途徑:一是延遲邊界層轉(zhuǎn)捩,又稱層流控制(laminar flow control,LFC),目的在于保持更大的層流范圍[1-2];二是發(fā)展湍流邊界層減阻技術(shù)[3]。本項研究關(guān)注的是第一種技術(shù)途徑。

    上游引入穩(wěn)定條帶是近幾年提出的一種新的延遲轉(zhuǎn)捩方法[4-5]。已有研究表明,上游引入穩(wěn)定條帶可以有效抑制T-S波的增長,從而使轉(zhuǎn)捩得到延遲[6-8]。條帶對T-S波的抑制作用強弱受條帶振幅的影響,振幅越大,抑制效果越好。但以往采用的粗糙元方法能夠生成的穩(wěn)定條帶的最大振幅只有12%U(U為自由流速度),遠小于理論預測的條帶穩(wěn)定的振幅上限閾值26%U[9]。這種方法之所以無法生成振幅更高的穩(wěn)定條帶,原因在于生成振幅更高的條帶需要更高的粗糙元,而粗糙元的高度增加后,局部繞流引起的渦脫落形式的擾動將直接引起邊界層旁路轉(zhuǎn)捩,條帶延遲轉(zhuǎn)捩作用失效。因此尋求有效生成大振幅穩(wěn)定條帶來提高延遲轉(zhuǎn)捩效果成為該流動控制技術(shù)研究關(guān)注的一大問題[10-11],最近的渦流發(fā)生器法生成的穩(wěn)定條帶最大振幅甚至高于30%U[11]。另一個問題是,盡管條帶對T-S波增長的抑制作用很明確,也有實驗給出引入條帶能夠抑制轉(zhuǎn)捩的實例,但這并不代表引入任何條帶都會延遲轉(zhuǎn)捩。新近的數(shù)值模擬和穩(wěn)定性研究結(jié)果也表明,不同的條帶參數(shù)(如條帶寬度)可能會得到完全不同的轉(zhuǎn)捩控制效果,即轉(zhuǎn)捩有可能會被延遲,也有可能會得到促進而提早發(fā)生[12-14]。原因在于T-S波增長只是轉(zhuǎn)捩過程中的起始階段,條帶對后續(xù)的二次不穩(wěn)定等擾動發(fā)展的影響,也會直接關(guān)系到最終的控制效果。因此需要開展不同條帶參數(shù)影響的實驗驗證和影響機制研究,找出能夠延遲轉(zhuǎn)捩的條帶參數(shù)及最優(yōu)控制參數(shù),為該流動控制技術(shù)的工程應用奠定良好基礎(chǔ)。

    在我們嘗試利用渦流法生成大振幅穩(wěn)定條帶的實驗中,偶然發(fā)現(xiàn)展向離散抽吸更能有效生成大振幅穩(wěn)定條帶,以此為基礎(chǔ)開展條帶控制轉(zhuǎn)捩實驗研究。本文實驗工作首先對未加任何擾動的平板邊界層基本流動性態(tài)進行了檢測;接下來引入固定頻率T-S波擾動激發(fā)邊界層轉(zhuǎn)捩,然后通過展向離散抽吸法引入條帶,利用氫氣泡時間線方法觀測引入條帶前后轉(zhuǎn)捩過程中擾動發(fā)展的變化,分析條帶對轉(zhuǎn)捩的控制效果,即抑制還是促進轉(zhuǎn)捩,包括條帶振幅和寬度等參數(shù)變化的影響。

    1 實驗裝置和方法

    實驗在北京航空航天大學回流式多用途低速水洞中完成。水洞實驗段尺寸1.2m×1.0m×16m(高×寬×長),最高流速1.0m/s。本次實驗來流速度保持在U=0.2m/s,此流速下實驗段湍流度和截面速度不均勻度均低于1%。實驗裝置見圖1。實驗平板立于水洞底壁,尺寸2.5m×1.05m×20mm(長×寬×厚),平板前緣截面為1∶3半橢圓,后緣裝有偏角β可調(diào)的尾鰭,用以控制前緣駐點位置,保證平板工作面上基本邊界層流動為穩(wěn)定的層流。

    激發(fā)邊界層轉(zhuǎn)捩的T-S波擾動通過距平板前緣xs=1000mm處的展向狹縫引入(見圖1(a)),狹縫寬1mm,長300mm。擾動引入裝置如圖1(b),擾動信號來自信號發(fā)生器生成的正弦波,通過功率放大器驅(qū)動揚聲器,引起揚聲器封閉內(nèi)腔容積的周期性變化,導致與其連通的狹縫處形成周期性零質(zhì)量射流,在邊界層內(nèi)可發(fā)展成二維T-S波擾動[15]。信號發(fā)生器輸出的正弦信號的頻率和幅值可調(diào),可用來改變引入邊界層內(nèi)T-S波的頻率和振幅。定義坐標如圖1(a)所示:x沿流向,x=0位于前緣;y沿平板工作面外法向,y= 0位于壁面;z為展向,z=0位于狹縫展向中心處。

    條帶通過距平板前緣xst=218mm的展向等距布置的抽吸孔引入,如圖1(c)所示,抽吸孔直徑4mm,與平板傾角為15°。抽吸通過虹吸原理實現(xiàn),各孔抽吸量Q保持相等,Q的大小可以通過改變各虹吸管出口的高度進行調(diào)節(jié),用以改變生成條帶的振幅。本次實驗中抽吸間距Δz主要取14、20和28mm 3種情況,在下游生成相同寬度(波長)的條帶。

    圖1 實驗裝置Fig.1 Sketches of experimental setup

    實驗觀測采用氫氣泡時間線法,時間線頻率12 Hz。主要采用俯視觀測方式,即生成氫氣泡的鉑絲沿展向(z向)布置于邊界層內(nèi)某一高度,觀測到的是此高度平行平板的流動形態(tài),可顯示出轉(zhuǎn)捩過程擾動發(fā)展的各種典型形式如T-S波、Λ結(jié)構(gòu)(發(fā)卡渦)以及流動紊亂化現(xiàn)象等。通過氫氣泡時間線法,還可以得到部分定量信息,如邊界層速度分布和條帶振幅等等。本次實驗氫氣泡時間線法測量流速和條帶振幅相對誤差均小于1%。

    2 實驗結(jié)果和分析

    2.1 基本邊界層流動

    在關(guān)于邊界層轉(zhuǎn)捩及其控制的實驗中,首要工作是需先明確未引入任何擾動時基本邊界層的流動性態(tài)。圖2(a)和2(b)分別給出本次實驗中氫氣泡時間線法測得的基本邊界層中不同流向位置速度分布和邊界層位移厚度δ*,測量結(jié)果與Blasius理論解相吻合。需要說明的是圖2(b)中實驗測得的邊界層位移厚度結(jié)果與理論結(jié)果相比向下游偏移Δx=43mm,這是由于半橢圓前緣區(qū)域存在的局部壓強梯度引起的,對下游的流動無進一步影響[15]。

    圖2 基本邊界層流動特性及其與Blasius結(jié)果的比較Fig.2 Characteristics of basic boundary layer flow and comparison with those of Blasius results

    基本邊界層流動在實驗觀測范圍內(nèi)(x=500~2000mm)的氫氣泡時間線俯視圖中都顯示出穩(wěn)定的層流特征:相鄰時間線等距且相互平行,沿展向無明顯變化,如圖3在x=2000mm即x-xs=1000mm處的流動顯示結(jié)果。本文中給出的氫氣泡流動顯示圖片中流動方向都是從左向右,圖片說明中參數(shù)x-xs和y/δ分別代表氫氣泡發(fā)生線(鉑絲)的流向和法向位置。其中δ為當?shù)剡吔鐚映叨?,定義為δ=(xν/U)1/2,ν為水的運動粘性系數(shù)。

    2.2 展向離散抽吸生成條帶

    本部分探討在基本邊界層中利用展向離散抽吸方法生成穩(wěn)定條帶的現(xiàn)象、機制和振幅發(fā)展特性,條帶寬度(抽吸孔間距)Δz=20mm。

    展向離散抽吸生成條帶機制如圖4(a)所示。在流向位置1(x=250mm)設(shè)置橫截面激光片光“切割”上游展向布置的鉑絲發(fā)出的氫氣泡面,結(jié)果見圖1(b),可以清晰的看到在每個抽吸孔下游展向兩側(cè)會出現(xiàn)一對反向旋轉(zhuǎn)的流向渦對。由于誘導作用,流向渦對可以有效地引起近壁區(qū)低動量的流體抬升,同時外層高動量的流體下掃,這就是所謂的“抬升效應”,在下游形成穩(wěn)定的條帶,如圖4(c),表現(xiàn)為高/低速區(qū)沿展向規(guī)則交替排列。

    圖3 氫氣泡時間線法顯示的層流邊界層流動(xxs=1000mm,y/δ=1.6)Fig.3 The laminar boundary layer flow visualized by hydrogen bubble timelines(x-xs=1000mm,y/δ=1.6)

    圖4 展向離散抽吸法生成條帶機制Fig.4 Sketch of the mechanisms by which spanwise discrete suctions generate streaks

    條帶沿流向發(fā)展主要體現(xiàn)為振幅的變化,條帶寬度(展向波長)保持不變。某一流向位置條帶的振幅Ast定義為條帶的一個展向波長內(nèi)最高和最低速度差值的一半在法向的極大值:

    實驗測得該極大值發(fā)生的法向高度y/δ=1.63(略低于粗糙元情況y/δ=2.2),沿流向基本保持不變。從實驗測得的條帶振幅沿流向變化結(jié)果圖5可看到,條帶振幅首先經(jīng)過短暫的增長之后,達到最大值,表現(xiàn)出典型的瞬時增長過程(transient growth),緊接著便逐漸衰減,這與粗糙元法生成的條帶相同。條帶的最大振幅決定于抽吸量,抽吸量越大,條帶振幅越大,如圖5中Q=2.83和0.944ml/s兩種情況比較。當抽吸量Q=2.83ml/s時,生成的穩(wěn)定條帶最大振幅可高達28.4%U(見圖5),遠高于相同條件下粗糙元方法能夠生成的最大振幅12%U,甚至超過理論預測的條帶穩(wěn)定的最大振幅26%U。當抽吸量繼續(xù)增大,條帶失穩(wěn),表現(xiàn)為非對稱破裂(亦稱正弦式破裂)。圖5還反映出盡管抽吸量不同,但生成的穩(wěn)定條帶的最大振幅都出現(xiàn)在幾乎相同的流向位置,對應無量綱波數(shù)β=2πδ/Δz=0.45左右,與理論預測及粗糙元生成條帶結(jié)果相同[6]。在相同的平板模型和來流條件下,我們也進行了粗糙元生成條帶實驗,實驗測得的穩(wěn)定條帶最大振幅約為12.3%U,與Fransson實驗吻合[6]。一方面驗證了氫氣泡時間線法測量條帶振幅的可靠性,更證實了展向抽吸方法生成大振幅穩(wěn)定條帶的能力。本部分實驗結(jié)果對于提高條帶控制轉(zhuǎn)捩能力,發(fā)展高效轉(zhuǎn)捩控制技術(shù)有參考意義。

    圖5 條帶振幅沿流向變化(Δz=20mm)Fig.5 Streamwise variation of the streak amplitude for streaks ofΔz=20mm

    2.3 條帶對人工激發(fā)邊界層轉(zhuǎn)捩的控制效果

    2.3.1 人工激發(fā)邊界層轉(zhuǎn)捩

    人工激發(fā)轉(zhuǎn)捩通過狹縫(見圖1)引入二維T-S波擾動實現(xiàn)。擾動頻率f=0.83 Hz,同文獻[15]。狹縫處擾動條件如圖6中“*”所示,對應無量綱參數(shù)R=426,F(xiàn)=1.43×10-4。其中為當?shù)乩字Z數(shù)平方根,F(xiàn)=2πfν/U2為無量綱頻率。“*”下游水平直線表示實驗觀測的流向范圍(R=426~602),對應x-xs=0~1000mm。作為參照,圖中同時給出中性穩(wěn)定性曲線,包括理論和實驗結(jié)果。通過多次測試,激發(fā)T-S波的信號發(fā)生器輸出電壓幅值選定為US=2.5V,一方面防止擾動過小邊界層內(nèi)T-S波發(fā)展不夠充分無法實現(xiàn)激發(fā)轉(zhuǎn)捩,另一方面避免擾動過強引起旁路轉(zhuǎn)捩。

    圖6 人工激發(fā)邊界層轉(zhuǎn)捩的擾動條件及中性穩(wěn)定性曲線Fig.6 The perturbation to excite boundary layer transition together with the neutral stability curve

    圖7的流動顯示結(jié)果給出激發(fā)邊界層轉(zhuǎn)捩過程中擾動發(fā)展表現(xiàn)出的典型流動形態(tài)。首先是T-S波的出現(xiàn),對應圖7(a)中左半部分顯示出的流態(tài),突出特證是沿流向時間線間距大小出現(xiàn)周期性變化,展向保持很好的二維性;然后,擾動向下游發(fā)展生成Λ結(jié)構(gòu)(發(fā)卡渦),如圖7(a)中箭頭A和B所指;最后在遠下游出現(xiàn)流動紊亂化,如圖7(b)所示。更多人工激發(fā)轉(zhuǎn)捩流動細節(jié)參見文獻[15]。

    圖7 人工激發(fā)轉(zhuǎn)捩過程中典型流動結(jié)構(gòu)(y/δ=1.6)Fig.7 The typical flow structures in the excited boundary layer transition(y/δ=1.6)

    2.3.2 條帶對邊界層轉(zhuǎn)捩的控制和條帶參數(shù)影響

    本部分實驗在上述人工激發(fā)轉(zhuǎn)捩邊界層中引入不同參數(shù)條帶,觀測條帶對轉(zhuǎn)捩的控制作用和條帶參數(shù)影響。選取Δz=14和28mm分析條帶寬度影響;在Δz=14mm條件下選取3個抽吸量Q來分析條帶振幅變化影響,各測試條帶的振幅見圖8,圖中只給出狹縫下游即引入T-S波后轉(zhuǎn)捩發(fā)展區(qū)的條帶幅值。

    圖8 控制轉(zhuǎn)捩的條帶的振幅Fig.8 The amplitudes of the steaks used to control the transition

    上述引入的都是穩(wěn)定條帶,即在基本邊界層中引入這些條帶不會發(fā)生條帶失穩(wěn)引起轉(zhuǎn)捩。圖9給出這些條帶被引入基本邊界層流動后的典型流動形態(tài)。圖9(a)給出在Δz=14mm和Q=1.36ml/s的條帶參數(shù)條件下,鉑絲位于x-xs=500和1000mm處的觀測結(jié)果。條帶流動的典型特征表現(xiàn)得十分明顯,即高/低速區(qū)沿展向交替排列。相比之下,在x-xs= 500mm處條帶振幅更大,展向流動變化更規(guī)則;在下游x-xs=1000mm處,由于條帶衰減,振幅明顯減小,且流動沿展向變化的規(guī)律性減弱,但仍保持穩(wěn)定的層流狀態(tài)。圖9(b)的流動顯示結(jié)果對應的是Δz= 28mm的條帶,流動特征與Δz=14mm條帶情況類似,但條帶寬度增加,振幅增大,且在x-xs=1000mm處仍能夠保持規(guī)整的波形。

    在2.3.1節(jié)所述的人工激發(fā)轉(zhuǎn)捩邊界層中引入寬度Δz=14mm、振幅(對應Q= 1.36ml/s)的條帶后,流動顯示結(jié)果見圖10。這里為狹縫xs=1000mm處條帶振幅(參見圖8),為方便起見,在下面的結(jié)果對比中統(tǒng)一將作為引入的條帶振幅大小的度量。圖10(b)與圖7(b)在x-xs=1000mm處的流動比較可以明顯看出,該條帶的引入完全抑制了人工激發(fā)轉(zhuǎn)捩在當?shù)氐牧鲃游蓙y化,表現(xiàn)出與基本流動中引入條帶類似的形態(tài)(圖9(a)中x-xs=1000mm處流動)。進一步,圖10(a)與圖7(a)在x-xs=200mm處流動比較,表明條帶能夠抑制發(fā)卡渦的發(fā)展,此時擾動表現(xiàn)為與引入的條帶相同展向波長的三維擾動,如圖10(a)中箭頭所示,這些擾動向下游不會進一步發(fā)展成發(fā)卡渦,而是又恢復成穩(wěn)定的條帶流動特征(圖10(b))。

    為分析條帶振幅變化對轉(zhuǎn)捩的影響,圖11(a)和(b)分別給出Δz=14mm條帶的振幅分別降至=10%U和7%U(分別對應Q=0.901和0.497ml/s)后,x-xs=1000mm處的流動顯示結(jié)果。二者之間以及連同圖10(b)的比較表明振幅增大,條帶對流動紊亂化的抑制作用增強,亦即延遲邊界層轉(zhuǎn)捩效果增強。

    為分析條帶寬度影響,圖12(a)給出引入Δz= 28mm條帶時x-xs=1000mm處的氫氣泡顯示結(jié)果,與圖7(b)(未引入條帶)和圖10(b)(Δz=14mm)的比較可以判斷出:Δz=28mm條帶對流動紊亂化同樣有抑制作用;但條帶寬度增大后,對流動紊亂化的抑制作用減弱,盡管Δz=28mm條帶的振幅(23%U)比Δz=14mm條帶的振幅(13%U)要大。注意到與圖10(b)相比,圖12(a)中條帶寬度增大對抑制流動紊亂化作用的減弱表現(xiàn)的不是十分顯著,為此,我們還嘗試了引入更寬的Δz=42mm條帶,結(jié)果如圖12(b)所示,在x-xs=1000mm處Δz=42mm條帶對流動紊亂化的抑制作用明顯減弱,甚至在個別局部地區(qū)流動紊亂化現(xiàn)象仍然存在,如圖12(b)中的2個箭頭所示位置。需要說明的是,生成Δz=42mm條帶的抽吸量與Δz=28mm情況相同(Q=2.53ml/s),但由于寬度Δz過大,生成的條帶未能得到充分發(fā)展,高速區(qū)流速存在明顯的等值段,所以無法明確地給出條帶振幅。

    圖9 穩(wěn)定條帶流動形態(tài)(y/δ=1.6)Fig.9 Visualized flow of the steady streaks

    圖10 轉(zhuǎn)捩邊界層中引入Δz=14mm和=13%U條帶后的流動(y/δ=1.6)Fig.10 The visualized flow in the transitional boundary layer after the streaks ofΔz=14mm andAst*=13%Uwere introduced

    圖11 Δz=14mm條帶不同振幅對轉(zhuǎn)捩過程流動紊亂化的影響(x-xs=1000mm,y/δ=1.6)Fig.11 Influence of streak width variation on the flow breakdown during transition(x-xs=1000mm,y/δ=1.6)

    圖12 不同條帶寬度對轉(zhuǎn)捩過程流動紊亂化的影響(x-xs=1000mm,y/δ=1.6)Fig.12 Influence of streak spacing variation on the flow breakdown during transition(x-xs=1000mm,y/δ=1.6)

    以上實驗研究結(jié)果表明利用展向離散抽吸方法引入穩(wěn)定條帶能夠有效抑制邊界層轉(zhuǎn)捩,從而為降低摩阻提供了一種可能的新途徑。需要指出的是,條帶的引入也會引起當?shù)啬Σ翍Φ脑黾印ransson等和Shahinfar等分別對此進行了估算[8,16],結(jié)果表明引入12%U中等振幅穩(wěn)定條帶后引起的摩擦阻力系數(shù)增大量小于3.5%,引入25%U大振幅穩(wěn)定條帶引起的摩擦阻力系數(shù)增大約12%,遠小于相同條件下湍流的摩擦應力(通常相比層流高一個量級)。因此,合理地設(shè)計引入條帶的位置(在邊界層轉(zhuǎn)捩位置上游不遠處),完全能夠有效降低總的摩擦阻力。此外,這種離散抽吸方法的實際應用還要考慮很多具體問題,如實現(xiàn)吸氣需增加相應裝置和控制系統(tǒng),由此會增加航行器載荷以及可能帶來設(shè)計上的麻煩,需要對整體能效比進行細致評估,這是所有主動控制技術(shù)都必須考慮的問題。本文僅在探索生成大振幅穩(wěn)定條帶的新方法及其對邊界層轉(zhuǎn)捩的控制效果方面開展了初步的、基礎(chǔ)性的研究工作。

    3 結(jié) 論

    發(fā)展了一種生成穩(wěn)定條帶的新方法,在此基礎(chǔ)上開展了條帶控制人工激發(fā)邊界層轉(zhuǎn)捩和條帶參數(shù)影響研究,主要結(jié)論:

    (1)展向離散抽吸是生成大振幅穩(wěn)定條帶的有效方法,能夠生成穩(wěn)定條帶的最大振幅可達28.4%U,遠高于粗糙元法生成的穩(wěn)定條帶的最大振幅12%U,甚至超過穩(wěn)定條帶的理論預測振幅上限26%U。

    (2)本實驗引入的兩種寬度的穩(wěn)定條帶都對轉(zhuǎn)捩過程中流動紊亂化進程有抑制效果,起到抑制轉(zhuǎn)捩的作用。條帶振幅增大,抑制作用增強;條帶寬度(展向波長)減小,對轉(zhuǎn)捩抑制效果更好。

    [1]王菲,額日其太,王強,等.基于升華法的后掠翼混合層流控制研究[J].實驗流體力學,2010,24(3):54-58.Wang Fei,Eriqitai,Wang Qiang,et al.Investigation of HLFC on swept wing based on sublimation technique[J].Journal of Experiments in Fluid Mechanics,2010,24(3):54-58.

    [2]鄧雙國,額日其太,聶俊杰.后掠翼模型混合層流控制實驗研究[J].實驗流體力學,2011,25(3):30-33.Deng Shuangguo,Eriqitai,Nie Junjie.Hybrid laminar flow control experiment on swept wing model[J].Journal of Experiments in Fluid Mechanics,2011,25(3):30-33.

    [3]盧浩,王兵,張會強,等.展向不同間距壁面粗糙元對近壁湍流擬序結(jié)構(gòu)的影響[J].實驗流體力學,2010,24(5):83-87. Lu Hao,Wang Bing,Zhang Huiqiang,et al.Effects of roughness elements in different spanwise interval arrayal on near-wall turbulence coherent structures[J].Journal of Experiments in Fluid Mechanics,2010,24(5):83-87.

    [4]Cossu C,Brandt L.Stabilization of Tollmien-Schlichting waves by finite amplitude optimal streaks in the Blasius boundary layer[J].Phys.Fluids,2002,14:L57-L60.

    [5]Cossu C,Brandt L.On Tollmien-Schlichting-like waves in streaky boundary layers[J].Eur J Mech/B Fluids,2004,23:815-833.

    [6]Fransson JHM,Brandt L,Talamelli A,et al.Experimental and theoretical investigation of the non-modal growth of steady streaks in a flat plate boundary layer[J].Phys Fluids,2004,16:3627-3638.

    [7]Fransson J H M,Brandt L,Talamelli A,et al.Experimental study of the stabilization of Tollmien-Schlichting waves by finite amplitude streaks[J].Phys Fluids,2005,17:054110.

    [8]Fransson J H M,Talamelli A,Brandt L,et al.Delaying transition to turbulence by a passive mechanism[J].Phys Rev Lett,2006,96:064501.

    [9]Andersson P,Brandt L,Bottaro A,et al.On the breakdown of boundary layers streaks[J].J Fluid Mech,2001,428:29-60.

    [10]Hollands M,Cossu C.Adding streamwise streaks in the plane Poiseuille flow[J].Comptes Rendus Mecanique,2009,337(3):179-183.

    [11]Fransson J H M,Talamelli A.On the generation of steady streamwise streaks in flat-plate boundary layers[J].J Fluid Mech,2012,698:211-234.

    [12]Liu Y,Zaki T A,Durbin P A.Boundary layer transition by interaction of discrete and continuous modes[J].J Fluid Mech,2008,604:199-233.

    [13]Liu Y,Zaki T A,Durbin P A.Floquet analysis of secondary instability of boundary layers distorted by Klebanoff streaks and Tollmien-Schlichting waves[J].Phys Fluids,2008,20:124102.

    [14]Durbin P A,Zaki T A,Liu Y.Interaction of discrete and continuous boundary layer modes to cause transition[J].International Journal of Heat and Fluid Flow,2009,30:403-410.

    [15]Guo H,Lian Q X,Li Y,et al.A visual study on complex flow structures and flow breakdown in boundary layer transition[J].Experiments in Fluids,2004,37:311-322.

    [16]Shahinfar S,Sattarzadeh SS,F(xiàn)ransson J H M,et al.Revival of classical vortex generators now for transition delay[J].Phys Rev Lett,2012,109:074501.

    Experimental study on boundary-layer transition control by spanwise discrete suction

    Guo Hui,Li Xiaobao,Wang Haiwen,F(xiàn)eng Yulong
    (1.School of Aeronautic Science and Engineering,Beihang University,Beijing 100191,China;2.Key Laboratory of Fluid Mechanics(Beihang University),Ministry of Education,Beijing 100191,China)

    Delaying boundary layer transition is an important way to reduce friction drag of vehicles in air and water.Introducing steady streaks upstream is a new idea proposed recently to delay transition.Available studies demonstrated that the steady streaks of higher amplitude have higher ability to depress T-S wave propagation,and thus are more effective in delaying transition.However,only relatively low amplitude steady streaks(12%U,whereUis the freestream velocity)can be obtained by the roughness elements method commonly used previously.This is because the generation of higher amplitude steady streaks needs higher roughness elements,which will lead to the bypass transition due to the stronger perturbation of vortex shedding caused by the elements.Moreover,recent theoretical studies show that there are some kinds of streaks wich do not suppress but instead promote the transition because of their acceleration of secondary instability.But this has not been verified by experiments yet.Thus,at least two problems are valuable to be studied for the flow control technique.One is to develop an effective method to generate higher amplitude steady streaks;the other is to test the control effect of different streaks on transition.In our studies,a new method,called spanwise discrete suction,has been developed to generate high amplitude steady streaks.Based on this,the streak control on the artificially excited transition in flat plate boundary layer has been studied using hydrogen bubble timeline method,by detecting the variation of perturbation evolution before and after the steady streaks introduced.The influences of streak spacing(corresponding to the suction hole spacing)and streak amplitude(depending on the suction intensity)are tested.The results show that spanwise discrete suction is an effective method to generate high amplitude streaks.The maximum amplitude of steady streaks generated by this method reaches as high as 28.4%U,evenhigher than the upper amplitude threshold value of 26%Upredicted theoretically for steady streaks.Moreover,both of the two steady streaks with spacing of 14mm and 28 mm introduced in our study suppress the flow breakdown,and thus delay the transition;The higher the streak amplitude or the narrower the streak spacing is,the more effective the streaks are in depressing the transition.

    streaks;boundary layer;transition;flow control;hydrogen bubble method

    O357.4;V211.1+9

    :A

    1672-9897(2014)06-0013-07doi:10.11729/syltlx20140027

    (編輯:李金勇)

    2014-03-04;

    :2014-10-15

    自然基金資助項目“條帶對邊界層轉(zhuǎn)捩的影響和轉(zhuǎn)捩控制實驗研究”(11072017)

    GuoH,LiXB,WangHW,etal.Experimentalstudyonboundary-layertransitioncontrolbyspanwisediscretesuction.JournalofExperimentsinFluidMechanics,2014,28(6):13-19.郭 輝,李小寶,王海文,等.展向離散抽吸法控制邊界層轉(zhuǎn)捩實驗研究.實驗流體力學,2014,28(6):13-19.

    郭 輝(1967-),男,河南唐河人,副教授。研究方向:實驗流體力學。通信地址:北京市海淀區(qū)學院路37號,北京航空航天大學流體所(100191)。E-mail:guohui@buaa.edu.cn

    猜你喜歡
    邊界層條帶振幅
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    基于條帶模式GEOSAR-TOPS模式UAVSAR的雙基成像算法
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    基于 Savitzky-Golay 加權(quán)擬合的紅外圖像非均勻性條帶校正方法
    中國光學(2015年1期)2015-06-06 18:30:20
    非特征邊界的MHD方程的邊界層
    鄭州市春季邊界層風氣候變化研究
    河南科技(2014年23期)2014-02-27 14:19:08
    av片东京热男人的天堂| 久久久久久久大尺度免费视频| 国产成人91sexporn| 亚洲伊人久久精品综合| 久久久久久久久久久久大奶| 大话2 男鬼变身卡| 久久精品国产a三级三级三级| 99热全是精品| 美女福利国产在线| 丁香六月欧美| 一边亲一边摸免费视频| 又粗又硬又长又爽又黄的视频| 十分钟在线观看高清视频www| 日韩人妻精品一区2区三区| 国产精品久久久久久人妻精品电影 | 丝瓜视频免费看黄片| 午夜精品国产一区二区电影| 好男人电影高清在线观看| 手机成人av网站| 欧美日韩精品网址| 亚洲精品国产av成人精品| 欧美国产精品一级二级三级| 成人影院久久| 成人18禁高潮啪啪吃奶动态图| 我的亚洲天堂| 校园人妻丝袜中文字幕| 真人做人爱边吃奶动态| 一级毛片我不卡| 亚洲精品国产av成人精品| 午夜激情av网站| 秋霞在线观看毛片| 黄网站色视频无遮挡免费观看| 亚洲五月色婷婷综合| 午夜精品国产一区二区电影| www.自偷自拍.com| 极品少妇高潮喷水抽搐| 在线观看人妻少妇| 波多野结衣av一区二区av| 成人国产一区最新在线观看 | 九草在线视频观看| www.av在线官网国产| 国产亚洲av片在线观看秒播厂| 精品亚洲成a人片在线观看| 日本五十路高清| 18禁观看日本| 男人舔女人的私密视频| videos熟女内射| av在线老鸭窝| 精品人妻在线不人妻| 日韩电影二区| 久久精品国产a三级三级三级| av线在线观看网站| 9热在线视频观看99| 亚洲欧美清纯卡通| 男人爽女人下面视频在线观看| 一级黄片播放器| 狠狠精品人妻久久久久久综合| 亚洲精品中文字幕在线视频| 性色av一级| 久久久久久久大尺度免费视频| 国产精品麻豆人妻色哟哟久久| 蜜桃在线观看..| 亚洲av日韩精品久久久久久密 | 久久久国产一区二区| 99热网站在线观看| 中文精品一卡2卡3卡4更新| bbb黄色大片| 99热网站在线观看| 国产无遮挡羞羞视频在线观看| 欧美少妇被猛烈插入视频| 日韩av在线免费看完整版不卡| 在线观看一区二区三区激情| 日韩制服丝袜自拍偷拍| 午夜视频精品福利| 一本大道久久a久久精品| 色婷婷av一区二区三区视频| 99热国产这里只有精品6| 久久人妻熟女aⅴ| 精品福利永久在线观看| 一二三四社区在线视频社区8| 99re6热这里在线精品视频| 99精国产麻豆久久婷婷| 岛国毛片在线播放| 午夜福利在线免费观看网站| 亚洲人成77777在线视频| 免费在线观看完整版高清| 一区二区日韩欧美中文字幕| 欧美变态另类bdsm刘玥| 亚洲av美国av| 女人久久www免费人成看片| 欧美日本中文国产一区发布| 国产又爽黄色视频| 欧美精品高潮呻吟av久久| 黑人巨大精品欧美一区二区蜜桃| 丁香六月天网| 99久久人妻综合| 国产日韩欧美在线精品| 亚洲av美国av| 19禁男女啪啪无遮挡网站| 满18在线观看网站| 一区在线观看完整版| 中文字幕人妻熟女乱码| 亚洲人成电影观看| 成人黄色视频免费在线看| a级毛片在线看网站| 久久久亚洲精品成人影院| 伊人亚洲综合成人网| 国产欧美日韩综合在线一区二区| 一本—道久久a久久精品蜜桃钙片| 久久久精品94久久精品| 两个人看的免费小视频| 亚洲av片天天在线观看| 亚洲精品日本国产第一区| 国产成人免费观看mmmm| 午夜激情久久久久久久| 久久 成人 亚洲| 下体分泌物呈黄色| 黄色 视频免费看| 999精品在线视频| 美女国产高潮福利片在线看| 美女扒开内裤让男人捅视频| 亚洲一区二区三区欧美精品| 伊人亚洲综合成人网| 国产精品一国产av| 18禁裸乳无遮挡动漫免费视频| 久久亚洲精品不卡| 啦啦啦在线免费观看视频4| 欧美人与性动交α欧美精品济南到| 99香蕉大伊视频| 国产精品亚洲av一区麻豆| 国产一区二区 视频在线| 一区二区三区乱码不卡18| 午夜两性在线视频| 亚洲成人免费电影在线观看 | 男的添女的下面高潮视频| 欧美精品人与动牲交sv欧美| 一区在线观看完整版| 成年女人毛片免费观看观看9 | 老司机午夜十八禁免费视频| 亚洲国产欧美在线一区| 国产片内射在线| 成人18禁高潮啪啪吃奶动态图| 男女边吃奶边做爰视频| 久9热在线精品视频| 国产欧美日韩精品亚洲av| 亚洲精品国产区一区二| 在线看a的网站| 成人午夜精彩视频在线观看| 亚洲久久久国产精品| 一本一本久久a久久精品综合妖精| 黄色一级大片看看| 一区二区日韩欧美中文字幕| 久久久国产精品麻豆| 婷婷色麻豆天堂久久| 国产一区亚洲一区在线观看| 高清不卡的av网站| 一区二区三区激情视频| 久久久国产欧美日韩av| 精品人妻熟女毛片av久久网站| 亚洲av日韩在线播放| 国产成人一区二区三区免费视频网站 | 成人国产一区最新在线观看 | 日韩av免费高清视频| 亚洲国产最新在线播放| 久久久久久久大尺度免费视频| 国产淫语在线视频| 久久女婷五月综合色啪小说| 老鸭窝网址在线观看| 一级片免费观看大全| 大香蕉久久成人网| 69精品国产乱码久久久| 成人国产av品久久久| 亚洲成人免费电影在线观看 | 国产精品成人在线| 国产精品免费大片| 男男h啪啪无遮挡| 国产精品一区二区精品视频观看| 999久久久国产精品视频| 超碰成人久久| 精品亚洲成a人片在线观看| 在线观看免费视频网站a站| 十八禁网站网址无遮挡| 性色av乱码一区二区三区2| 久久精品亚洲av国产电影网| netflix在线观看网站| 国产真人三级小视频在线观看| 一区二区三区精品91| 91精品伊人久久大香线蕉| 18禁观看日本| 捣出白浆h1v1| 亚洲国产欧美一区二区综合| 777米奇影视久久| 精品久久久精品久久久| 中文字幕亚洲精品专区| 少妇被粗大的猛进出69影院| 亚洲黑人精品在线| 91麻豆av在线| a 毛片基地| 欧美精品高潮呻吟av久久| 国产午夜精品一二区理论片| 人妻一区二区av| 男女床上黄色一级片免费看| 久久中文字幕一级| 久久免费观看电影| 丝袜人妻中文字幕| av网站在线播放免费| 一区二区三区四区激情视频| 一边摸一边抽搐一进一出视频| 欧美日韩黄片免| 在线观看免费视频网站a站| 欧美激情极品国产一区二区三区| 男女边摸边吃奶| 亚洲成国产人片在线观看| 美女扒开内裤让男人捅视频| 久久青草综合色| 欧美日韩av久久| videos熟女内射| 黄色视频不卡| 国产精品成人在线| 亚洲欧美日韩高清在线视频 | 亚洲欧洲国产日韩| 婷婷丁香在线五月| 精品久久久久久电影网| 高清不卡的av网站| 999久久久国产精品视频| 五月开心婷婷网| 别揉我奶头~嗯~啊~动态视频 | 脱女人内裤的视频| 亚洲欧美日韩另类电影网站| 国产又爽黄色视频| 大型av网站在线播放| 日韩av免费高清视频| 中文乱码字字幕精品一区二区三区| 嫩草影视91久久| 91麻豆精品激情在线观看国产 | 免费在线观看影片大全网站 | 只有这里有精品99| 亚洲精品乱久久久久久| 久久久久久久精品精品| 日韩免费高清中文字幕av| 丝瓜视频免费看黄片| 80岁老熟妇乱子伦牲交| 成年美女黄网站色视频大全免费| 国产黄色免费在线视频| 一本色道久久久久久精品综合| 一级毛片女人18水好多 | www.av在线官网国产| 久久99一区二区三区| 精品国产乱码久久久久久小说| 亚洲精品乱久久久久久| 国产人伦9x9x在线观看| 91精品伊人久久大香线蕉| 亚洲专区国产一区二区| 国产成人欧美| 熟女av电影| 日本午夜av视频| 校园人妻丝袜中文字幕| 在线观看www视频免费| 在线亚洲精品国产二区图片欧美| 无遮挡黄片免费观看| 成年人免费黄色播放视频| 视频在线观看一区二区三区| 18禁黄网站禁片午夜丰满| 熟女少妇亚洲综合色aaa.| av国产精品久久久久影院| 亚洲精品一卡2卡三卡4卡5卡 | 只有这里有精品99| 黄色片一级片一级黄色片| 狂野欧美激情性xxxx| 在线观看免费高清a一片| 一级毛片我不卡| 国产亚洲一区二区精品| 精品一区二区三卡| 男女免费视频国产| svipshipincom国产片| 亚洲国产精品成人久久小说| 国产真人三级小视频在线观看| 国产精品久久久久久人妻精品电影 | 爱豆传媒免费全集在线观看| 我要看黄色一级片免费的| 国产精品香港三级国产av潘金莲 | 丰满少妇做爰视频| av线在线观看网站| www.999成人在线观看| 大香蕉久久网| 天天躁夜夜躁狠狠久久av| 1024视频免费在线观看| 中文字幕人妻熟女乱码| 国产亚洲午夜精品一区二区久久| 亚洲自偷自拍图片 自拍| 不卡av一区二区三区| 久热爱精品视频在线9| 97人妻天天添夜夜摸| 日韩电影二区| 香蕉丝袜av| 精品欧美一区二区三区在线| 亚洲自偷自拍图片 自拍| 午夜福利乱码中文字幕| 欧美在线一区亚洲| 99热网站在线观看| 在线观看免费高清a一片| 国产亚洲精品第一综合不卡| 美女大奶头黄色视频| 男女国产视频网站| 女性被躁到高潮视频| 嫩草影视91久久| 亚洲成av片中文字幕在线观看| 欧美人与善性xxx| 国产又爽黄色视频| 久热这里只有精品99| 老司机影院毛片| 欧美黄色淫秽网站| 欧美 亚洲 国产 日韩一| 蜜桃国产av成人99| 人妻 亚洲 视频| 亚洲一码二码三码区别大吗| 操美女的视频在线观看| 男人爽女人下面视频在线观看| 1024香蕉在线观看| 成人影院久久| 国产又色又爽无遮挡免| 男女边吃奶边做爰视频| 两个人免费观看高清视频| 一本色道久久久久久精品综合| 乱人伦中国视频| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 视频在线观看一区二区三区| 99热全是精品| 亚洲五月婷婷丁香| 亚洲精品在线美女| 波多野结衣一区麻豆| 丰满饥渴人妻一区二区三| 亚洲av国产av综合av卡| 亚洲三区欧美一区| 超色免费av| 国产黄色视频一区二区在线观看| 亚洲一区二区三区欧美精品| 国产日韩欧美视频二区| 中文欧美无线码| 一区二区三区激情视频| 精品亚洲成a人片在线观看| 美女大奶头黄色视频| 在线 av 中文字幕| 大话2 男鬼变身卡| 老司机在亚洲福利影院| 麻豆av在线久日| 亚洲精品一卡2卡三卡4卡5卡 | 99国产综合亚洲精品| 国产成人免费无遮挡视频| 日韩,欧美,国产一区二区三区| 丰满迷人的少妇在线观看| kizo精华| 精品一区在线观看国产| 尾随美女入室| 脱女人内裤的视频| 日韩熟女老妇一区二区性免费视频| 日日夜夜操网爽| 欧美日本中文国产一区发布| 在线观看免费高清a一片| 亚洲中文av在线| 欧美日韩视频高清一区二区三区二| 91麻豆精品激情在线观看国产 | 韩国高清视频一区二区三区| 在线观看免费日韩欧美大片| 欧美亚洲日本最大视频资源| 日本午夜av视频| 欧美成人精品欧美一级黄| 日本午夜av视频| 少妇粗大呻吟视频| 免费在线观看黄色视频的| 三上悠亚av全集在线观看| 日本色播在线视频| 午夜视频精品福利| 国产一区二区在线观看av| 成年av动漫网址| 亚洲欧洲国产日韩| 91成人精品电影| 18禁裸乳无遮挡动漫免费视频| 少妇精品久久久久久久| 精品欧美一区二区三区在线| 欧美黑人精品巨大| 男女高潮啪啪啪动态图| 亚洲精品一二三| 9色porny在线观看| 黄色一级大片看看| 丰满少妇做爰视频| 久久99一区二区三区| 亚洲国产精品成人久久小说| 中国美女看黄片| 亚洲专区国产一区二区| 国产色视频综合| 美女视频免费永久观看网站| 男人操女人黄网站| 亚洲精品国产av蜜桃| 黄色片一级片一级黄色片| 亚洲av在线观看美女高潮| 国产成人91sexporn| 女人高潮潮喷娇喘18禁视频| 在线av久久热| 精品一区二区三卡| av天堂久久9| 男女床上黄色一级片免费看| 国产1区2区3区精品| 看十八女毛片水多多多| 又黄又粗又硬又大视频| 电影成人av| 一边摸一边做爽爽视频免费| 亚洲美女黄色视频免费看| 亚洲av综合色区一区| 后天国语完整版免费观看| 国产成人欧美| 色婷婷av一区二区三区视频| 日韩 欧美 亚洲 中文字幕| 欧美中文综合在线视频| 人人妻人人爽人人添夜夜欢视频| 久久精品久久精品一区二区三区| 成在线人永久免费视频| 不卡av一区二区三区| 久久久久久久久久久久大奶| 热re99久久精品国产66热6| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产av蜜桃| 老司机午夜十八禁免费视频| 美女中出高潮动态图| 男女国产视频网站| 成年人免费黄色播放视频| 国产成人精品久久二区二区免费| 国产成人精品久久久久久| 九色亚洲精品在线播放| 91国产中文字幕| 在线亚洲精品国产二区图片欧美| 日韩大片免费观看网站| 精品人妻在线不人妻| 中文字幕最新亚洲高清| 黄色一级大片看看| 我要看黄色一级片免费的| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 日本av免费视频播放| 在线亚洲精品国产二区图片欧美| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三 | 欧美人与性动交α欧美软件| 亚洲激情五月婷婷啪啪| 欧美 日韩 精品 国产| 免费观看a级毛片全部| 九色亚洲精品在线播放| 国产精品偷伦视频观看了| 两个人看的免费小视频| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 欧美日韩亚洲高清精品| 亚洲av欧美aⅴ国产| 美国免费a级毛片| videosex国产| 国产成人av激情在线播放| 精品一品国产午夜福利视频| 中文字幕最新亚洲高清| 成人免费观看视频高清| 国产精品久久久久成人av| 日韩av不卡免费在线播放| 亚洲三区欧美一区| 精品一区二区三区四区五区乱码 | 免费在线观看黄色视频的| 无遮挡黄片免费观看| svipshipincom国产片| 一二三四社区在线视频社区8| 国产免费一区二区三区四区乱码| 精品免费久久久久久久清纯 | 久久久久久亚洲精品国产蜜桃av| 天天躁夜夜躁狠狠躁躁| 精品国产乱码久久久久久小说| 婷婷色av中文字幕| cao死你这个sao货| 亚洲国产最新在线播放| 国产男女超爽视频在线观看| 高清不卡的av网站| 十八禁网站网址无遮挡| 天堂8中文在线网| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 男女无遮挡免费网站观看| 国产精品亚洲av一区麻豆| 亚洲精品一区蜜桃| 亚洲av成人精品一二三区| 欧美成人精品欧美一级黄| 别揉我奶头~嗯~啊~动态视频 | 精品福利观看| 久久精品国产亚洲av高清一级| 99国产精品一区二区三区| 精品国产一区二区三区久久久樱花| 午夜日韩欧美国产| 岛国毛片在线播放| 我的亚洲天堂| 亚洲第一青青草原| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| av网站在线播放免费| 日韩视频在线欧美| 国产97色在线日韩免费| 国产精品久久久av美女十八| 国产精品人妻久久久影院| av欧美777| 一级毛片电影观看| 80岁老熟妇乱子伦牲交| 成人影院久久| 蜜桃在线观看..| 激情五月婷婷亚洲| 久久久亚洲精品成人影院| 亚洲午夜精品一区,二区,三区| 精品福利永久在线观看| 美女国产高潮福利片在线看| 大型av网站在线播放| 热99国产精品久久久久久7| 国产男女超爽视频在线观看| 99九九在线精品视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产av成人精品| 色94色欧美一区二区| 亚洲国产最新在线播放| 性色av乱码一区二区三区2| 免费在线观看视频国产中文字幕亚洲 | √禁漫天堂资源中文www| 国产日韩欧美亚洲二区| 国产男人的电影天堂91| 国产精品香港三级国产av潘金莲 | 在线观看免费视频网站a站| 一级毛片电影观看| 18禁国产床啪视频网站| 久久人人爽av亚洲精品天堂| 在线精品无人区一区二区三| 精品少妇久久久久久888优播| 看免费成人av毛片| 日日夜夜操网爽| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 欧美日韩综合久久久久久| 久久人人爽av亚洲精品天堂| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一卡2卡三卡4卡5卡 | 丝袜在线中文字幕| 欧美日韩一级在线毛片| 国产无遮挡羞羞视频在线观看| 成人黄色视频免费在线看| 欧美成人精品欧美一级黄| 久久九九热精品免费| 亚洲精品在线美女| 国产熟女午夜一区二区三区| 亚洲av综合色区一区| 一区在线观看完整版| 女人被躁到高潮嗷嗷叫费观| 制服诱惑二区| 热re99久久精品国产66热6| 后天国语完整版免费观看| 亚洲国产欧美日韩在线播放| 欧美精品人与动牲交sv欧美| 一边摸一边抽搐一进一出视频| 高清黄色对白视频在线免费看| 免费日韩欧美在线观看| 欧美在线一区亚洲| 老司机影院成人| 一级片免费观看大全| 一区二区三区激情视频| 午夜av观看不卡| 国产精品久久久av美女十八| av在线播放精品| 久久免费观看电影| 国产在线一区二区三区精| 国产高清videossex| 天天添夜夜摸| 亚洲成人免费电影在线观看 | 欧美激情 高清一区二区三区| 咕卡用的链子| 91国产中文字幕| 精品一区二区三卡| 国产成人精品久久二区二区免费| 嫁个100分男人电影在线观看 | 黄色a级毛片大全视频| 欧美xxⅹ黑人| 午夜视频精品福利| 男女国产视频网站| 免费女性裸体啪啪无遮挡网站| 少妇粗大呻吟视频| 亚洲精品国产色婷婷电影| 成人亚洲欧美一区二区av| 国产精品一二三区在线看| 中文精品一卡2卡3卡4更新| 国产免费现黄频在线看| 一级黄色大片毛片| 久久久久视频综合| 女人久久www免费人成看片| 亚洲熟女毛片儿| 爱豆传媒免费全集在线观看| 国产一区二区激情短视频 | 欧美黑人欧美精品刺激| www日本在线高清视频| 久久久精品国产亚洲av高清涩受| a级片在线免费高清观看视频| 高清欧美精品videossex| 99热国产这里只有精品6| 亚洲成人国产一区在线观看 | 后天国语完整版免费观看| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美网| 免费人妻精品一区二区三区视频| 男人添女人高潮全过程视频| 精品国产国语对白av| 一区在线观看完整版| 久久毛片免费看一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 精品福利永久在线观看| 国产深夜福利视频在线观看| 国产日韩一区二区三区精品不卡| 中文字幕另类日韩欧美亚洲嫩草|