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

    湍流邊界層外區(qū)超大尺度相干結(jié)構(gòu)相位平均波形

    2016-07-25 08:33:23白建俠鄭小波姜楠
    實(shí)驗(yàn)流體力學(xué) 2016年5期
    關(guān)鍵詞:法向邊界層流向

    白建俠,鄭小波,姜楠,3,*

    (1.天津大學(xué)機(jī)械工程學(xué)院,天津300354;2.天津大學(xué)仁愛學(xué)院數(shù)學(xué)教學(xué)部,天津301636;3.天津市現(xiàn)代工程力學(xué)重點(diǎn)實(shí)驗(yàn)室,天津300350)

    湍流邊界層外區(qū)超大尺度相干結(jié)構(gòu)相位平均波形

    白建俠1,2,鄭小波1,姜楠1,3,*

    (1.天津大學(xué)機(jī)械工程學(xué)院,天津300354;2.天津大學(xué)仁愛學(xué)院數(shù)學(xué)教學(xué)部,天津301636;3.天津市現(xiàn)代工程力學(xué)重點(diǎn)實(shí)驗(yàn)室,天津300350)

    用IFA300恒溫式熱線風(fēng)速儀和X形二分量熱線探針精細(xì)測量了風(fēng)洞中平板湍流邊界層不同法向位置的瞬時流向、法向速度分量的時間序列信號。通過對不同速度分量進(jìn)行子波變換,得到了近壁區(qū)和外邊界區(qū)流向和法向脈動速度分量的分尺度湍動能隨尺度的分布,在外邊界區(qū)發(fā)現(xiàn)流向脈動速度和法向脈動速度的能量最大尺度隨法向坐標(biāo)遠(yuǎn)離壁面增大,相干結(jié)構(gòu)的平均猝發(fā)周期也相應(yīng)增加,說明在外區(qū)存在超大尺度的相干結(jié)構(gòu)。采用條件采樣和相位平均技術(shù),提取了近壁區(qū)和外邊界區(qū)能量最大尺度相干結(jié)構(gòu)流向及法向速度的噴射和掃掠的條件相位平均波形,發(fā)現(xiàn)流向和法向脈動速度及雷諾應(yīng)力噴射和掃掠事件的條件相位波形較近壁區(qū)都發(fā)生了改變。雷諾應(yīng)力負(fù)的幅值顯著降低,雷諾應(yīng)力表現(xiàn)為正負(fù)交替的波動現(xiàn)象。

    湍流邊界層;超大尺度結(jié)構(gòu);子波分析;噴射

    0 引 言

    湍流邊界層中存在著大尺度的相干結(jié)構(gòu)。通過實(shí)驗(yàn)發(fā)現(xiàn),大尺度的相干結(jié)構(gòu)可分為2類:一類存在于湍流邊界層的近壁區(qū),與近壁區(qū)的猝發(fā)現(xiàn)象和湍流產(chǎn)生有關(guān)[1],另一類存在于邊界層的外邊界區(qū)[2]。這種湍流邊界層外區(qū)的超大尺度條紋結(jié)構(gòu)由Marusic et al.[3]發(fā)現(xiàn)。Hutchins et al.[4]在湍流邊界層的對數(shù)律區(qū)和尾流區(qū)也發(fā)現(xiàn)了超大尺度低速和高速條紋結(jié)構(gòu),該條紋結(jié)構(gòu)流向長度可達(dá)到20倍的湍流邊界層厚度δ,而條紋最大寬度約為0.5δ。由于這種結(jié)構(gòu)具有很長的流向尺度,且這種結(jié)構(gòu)對湍流邊界層能量分布具有重要作用,因此將其稱為“Superstructures”。這種超大尺度結(jié)構(gòu)同時也在影響邊界層內(nèi)區(qū)的運(yùn)動。

    Adrian[5]認(rèn)為邊界層外區(qū)的超大尺度相干結(jié)構(gòu)通常以一組空間關(guān)聯(lián)程度較高的發(fā)卡渦(包)的形式存在,發(fā)卡渦(包)貫穿整個邊界層。Liu et al.[6]指出,超大尺度湍流結(jié)構(gòu)對湍流能量輸運(yùn)起主導(dǎo)作用,占有總湍動能一半以上,而且對湍流邊界層外區(qū)的雷諾應(yīng)力的貢獻(xiàn)高達(dá)2/3至3/4,對湍流的產(chǎn)生和維持有重要貢獻(xiàn)[7-11]。Guala et al.[12]通過對管道湍流中的“Superstructures”進(jìn)行實(shí)驗(yàn)研究發(fā)現(xiàn):除了發(fā)卡渦包,一定還存在其它原因?yàn)閷?shù)律區(qū)超大尺度結(jié)構(gòu)(>20δ)的形成起重要作用,而且湍流邊界層外區(qū)超大尺度結(jié)構(gòu)(>20δ)的形成機(jī)理與近壁區(qū)的大尺度結(jié)構(gòu)的形成機(jī)理不同。

    子波分析是近些年新發(fā)展起來的一種數(shù)學(xué)方法[1314],通過信號與一個被稱為子波的解析函數(shù)進(jìn)行卷積,將信號在時域與頻域空間同時進(jìn)行分解。姜楠等[15]提出了用子波分析的能量最大準(zhǔn)則檢測壁湍流相干結(jié)構(gòu)猝發(fā)事件的方法。

    一維速度信號u(t)在子波函數(shù)Wab(t)下的子波分析Wu(a,b)定義為:

    其中子波函數(shù)族Wab(t)是由子波母函數(shù)W(t)經(jīng)過平移(參數(shù)b)和伸縮(參數(shù)a)變換而來:

    子波分析又被稱為數(shù)學(xué)顯微鏡,其“放大倍數(shù)”為參數(shù)a,“焦點(diǎn)位置”為參數(shù)b,子波母函數(shù)W(t)的選擇則決定了它的“光學(xué)性質(zhì)”。

    根據(jù)子波系數(shù)Wu(a,b),速度信號u(t)的能量可以分解為

    姜楠等[15]提出了用子波分析的能量最大準(zhǔn)則確定壁湍流相干結(jié)構(gòu)猝發(fā)周期時間尺度的方法。在E(a)隨尺度a的分布中,每個尺度的湍流脈動速度所占有的湍流脈動動能是不同的,存在著一個峰值,對應(yīng)能量最大的尺度a*,該尺度就是相干結(jié)構(gòu)猝發(fā)周期對應(yīng)的時間尺度。

    本文中用到的小波母函數(shù)為二次B-樣條函數(shù)。

    零次B-樣條函數(shù)就是[0,1)區(qū)間的特征函數(shù),記作:

    各尺度信號所占有能量的總和等于信號的總能量。其中

    m次基數(shù)B-樣條函數(shù)Nm(x)由如下卷積迭代關(guān)系定義:

    顯然,Nm(x)是緊支的,其支撐寬度為

    其支撐中心為(m+1)/2,也是Nm(x)的對稱中心??梢缘玫?,二次B-樣條函數(shù)的數(shù)學(xué)表達(dá)式為:

    由于Nm(x)可以作為一個尺度函數(shù)而生成一個多分辨率分析,并且是一類分段光滑又在各段交界處具有一定光滑性的函數(shù)?;鶖?shù)B-樣條函數(shù)Nm(x)具有最小可能的支撐寬度,又便于算法的計(jì)算機(jī)處理和實(shí)時實(shí)現(xiàn)。不同子波母函數(shù)得到的分析結(jié)果有一定的差別,本研究選取的二次B-樣條小波和渦結(jié)構(gòu)具有較好的相似性,可以得到較好的結(jié)果。

    本文采用熱線測速實(shí)驗(yàn)技術(shù),在風(fēng)洞中開展湍流邊界層外區(qū)超大尺度相干結(jié)構(gòu)的實(shí)驗(yàn)研究。用子波變換對湍流邊界層中不同法向位置的流向和法向瞬時速度信號的時間序列進(jìn)行多尺度特征分析。目的是希望用子波分析的方法驗(yàn)證邊界層的外邊界區(qū)超大尺度相干結(jié)構(gòu)的存在性,提取湍流邊界層外區(qū)流向脈動速度和法向脈動速度的分尺度湍動能隨尺度和法向位置變化的規(guī)律和特征,提取超大尺度相干結(jié)構(gòu)的噴射和掃掠事件各物理量的條件相位平均波形,分析超大尺度相干結(jié)構(gòu)的相位平均波形與近壁湍流相干結(jié)構(gòu)的相位平均波形之間的差異。

    1 實(shí)驗(yàn)裝置

    實(shí)驗(yàn)在天津大學(xué)流體力學(xué)實(shí)驗(yàn)室的低速回流式木質(zhì)風(fēng)洞中進(jìn)行,實(shí)驗(yàn)段長1500mm,橫截面為切角的矩形,寬800mm,高600mm。實(shí)驗(yàn)段風(fēng)速在0.5~40m/s之間連續(xù)可調(diào),背景湍流度小于0.2%。在風(fēng)洞中沿中心線豎直放置一塊長1700mm、寬600mm、厚15mm的有機(jī)玻璃平板,平板前緣為8∶1橢圓。以平板前緣為原點(diǎn)建立空間直角坐標(biāo)系,ox軸正向表示流向,oy軸正向?yàn)榉ㄏ颍琽z軸正向?yàn)檎瓜颉T谄桨迩熬墄=8cm處安裝拌線和砂紙,用來加速轉(zhuǎn)捩并保證在其下游熱線探針測量點(diǎn)獲得充分發(fā)展的湍流邊界層(實(shí)驗(yàn)裝置見圖1)。

    圖1 實(shí)驗(yàn)裝置示意圖Fig.1 Schematic diagram of the experimental setup

    使用美國TSI公司生產(chǎn)的IFA300恒溫式熱線風(fēng)速儀以及TSI-1234-T1.5型雙絲邊界層熱線探針進(jìn)行測量。這種熱線探針的2根熱線在空間為互相垂直異面直線布置,熱線探針的敏感材料為長1.25mm,直徑4μm的鎢絲。熱線探針支桿固定在CCTS-1193E型三維自動控制坐標(biāo)架上。探針在距離平板前緣1090cm處的流向位置測量流向速度分量u和法向速度分量v。實(shí)驗(yàn)中選取130個不等間距的法向位置,距離壁面最近的位置4mm,探針移動的步長在近壁區(qū)為0.1mm,外區(qū)為2mm。

    由于在測量速度脈動信號時采樣頻率要高于湍流的最小時間尺度(耗散尺度),采樣時間要長于湍流的最大時間尺度(積分尺度),設(shè)定采樣頻率為100k Hz,低通截止頻率為50k Hz。每個測量點(diǎn)采集的樣本總量為4 194 304個,采樣時間為42s。實(shí)驗(yàn)中自由來流速度9.59m/s,邊界層的名義厚度為39.88mm,基于邊界層名義厚度的雷諾數(shù)為Reδ=24 092,基于動量損失厚度的雷諾數(shù)為Reθ=2183。

    2 實(shí)驗(yàn)結(jié)果分析

    圖2為由黏性內(nèi)尺度單位無量綱化的平板湍流邊界層的平均速度剖面。其中分別代表由摩擦速度和運(yùn)動黏性系數(shù)無量綱化的流向平均速度和法向位置,紅色線和藍(lán)色線分別表示單絲和雙絲測得的流向平均速度剖面。由于雙絲熱線探針的敏感元件與壁面呈45°傾斜放置,測點(diǎn)不能離壁面很近,從圖中可看出從緩沖層開始,單絲探針和雙絲探針的測量結(jié)果符合較好。本文主要分析雙絲測量的實(shí)驗(yàn)數(shù)據(jù)。

    圖2 平板湍流邊界層平均速度剖面Fig.2 Mean velocity profile of turbulent boundary layer

    用子波分析對湍流脈動速度信號進(jìn)行多尺度分解,把湍流流向和法向脈動速度信號分別分解成21個尺度。得到了子波系數(shù)在時間和尺度上的等值線圖,如圖3所示。不同顏色表示不同子波系數(shù)幅值。由圖3可知,不同尺度上存在著一個又一個擬序的湍渦結(jié)構(gòu),其占有較多的湍流脈動動能。圖中分別選取了近壁區(qū)y+=28.1和外區(qū)y+=679.3這2個不同的法向位置??梢钥闯?,圖中分布著不同尺度的渦結(jié)構(gòu)。在近壁區(qū)y+=28.1處,流向的能量最大尺度為第9尺度,在外區(qū)y+=679.3處,流向的能量最大尺度從第9尺度向第10尺度偏移,說明在外區(qū),相干結(jié)構(gòu)的流向尺度變大。

    圖4為不同法向位置多尺度流動結(jié)構(gòu)的能量分布線圖。圖4(a)為近壁區(qū)流向脈動速度能量隨尺度分布圖;圖4(b)為外區(qū)流向脈動速度能量隨尺度分布圖;圖4(c)為近壁區(qū)法向脈動速度能量隨尺度分布圖;圖4(d)為外區(qū)法向脈動速度能量隨尺度分布圖,圖4(e)為法向脈動速度近壁區(qū)過渡到外區(qū)能量隨尺度分布圖。由圖4可以看出,近壁區(qū)不同尺度的渦能量隨著法向位置遠(yuǎn)離壁面逐漸減少。表明緩沖層中的相干結(jié)構(gòu)最為活躍,占有特別重要的地位,是湍流產(chǎn)生的主要區(qū)域。

    圖3 近壁區(qū)和外區(qū)流向及法向速度信號的子波系數(shù)等值線圖Fig.3 Wavelet coefficient contour of streamwise velocity component and normal velocity component in the near-wall and outer region

    近壁區(qū)流向脈動速度的能量最大尺度都在第9尺度,外區(qū)y+=679.3處能量最大尺度仍保持在第9尺度,在y+=847.7處能量最大尺度變化到第10尺度。對于法向脈動速度,外區(qū)能量最大尺度在y+=679.3處由第8尺度變化到第9尺度,在y+=847.7處能量最大尺度從第9尺度變化到第10尺度。盡管遠(yuǎn)離壁面區(qū)域尺度增大,但是能量卻逐漸減弱,這種尺度和能量的變化趨勢與Adrian[5]和Ganapathisubramani B.[10]的描述是一致的。

    圖4 流法向脈動速度近壁區(qū)和外區(qū)能量隨尺度分布線圖Fig.4 Fluctuating kinetic energy distribution of streamwise velocity and normal velocity component in the nearwall and outer region

    用子波系數(shù)模極值法[16-17]提取多尺度相干結(jié)構(gòu)湍渦結(jié)構(gòu)。將流向脈動速度的子波系數(shù)是否為負(fù)的極小值作為判斷是否為噴射事件的標(biāo)準(zhǔn);子波系數(shù)是否為正的極大值作為判斷是否為掃掠事件的標(biāo)準(zhǔn)??梢酝ㄟ^測量湍流邊界層的近壁區(qū)域的脈動速度的時間序列信號,總結(jié)出相干結(jié)構(gòu)的基本特征,以此為根據(jù)制定相應(yīng)的判定準(zhǔn)則,依據(jù)這個判定準(zhǔn)則檢測相干結(jié)構(gòu),提取出相干結(jié)構(gòu)的條件相位平均波形。

    式中:De(aj,ti+t)和Ds(aj,ti+t)分別為噴射和掃掠事件的檢測函數(shù);〈f(aj,t)〉e和〈f(aj,t)〉s分別表示噴射和掃掠事件中各物理量對應(yīng)的條件相位波形,aj為能量最大尺度,E(aj)=max E(a)。

    圖5 近壁區(qū)不同法向位置流向及法向脈動速度噴射掃掠事件的條件相位平均波形Fig.5 Phase-averaging waveforms of streamwise velocity and normal velocity components in the near-wall region at ejection and sweep events

    圖5 給出了在近壁區(qū)不同法向位置的能量最大尺度對應(yīng)的噴射和掃掠事件相干結(jié)構(gòu)的流向和法向脈動速度的條件相位平均波形。從圖5中可以看出,在不同的法向位置,噴射和掃掠事件的條件相位波形幅值發(fā)生了改變。隨著法向位置增加,相干結(jié)構(gòu)的流向和法向脈動速度條件相位平均波形的幅值呈減少趨勢,表明噴射與掃掠事件的強(qiáng)度減弱。

    圖6 外區(qū)不同法向位置流向及法向脈動速度噴射掃掠事件的條件相位平均波形Fig.6 The phase-averaging waveforms of the streamwise velocity and normal velocity components in the outer region

    湍流脈動噴射和掃掠事件條件相位波形的周期也發(fā)生了變化。外區(qū)隨著能量最大尺度增加,噴射和掃掠事件的時間周期增加。流向脈動速度在外區(qū)y+=847.7處能量最大尺度由第9尺度變化到第10尺度,表明噴射和掃掠時間的尺度增加;法向脈動速度在外區(qū)y+=679.3和y+=847.7處能量最大尺度由第8尺度變化到第9和第10尺度,相應(yīng)的噴射和掃掠的時間周期也增加。說明在外區(qū)存在著尺度不斷增大的相干結(jié)構(gòu),這些超大尺度相干結(jié)構(gòu)的時間周期比近壁區(qū)大尺度結(jié)構(gòu)時間周期要大很多。隨著遠(yuǎn)離壁面,低速度條帶區(qū)域變寬、變長,這與Del Alamoa et al.[18]的發(fā)現(xiàn)是一致的。

    在噴射過程中,首先流向速度有一個緩慢的加速過程,然后是一個急劇的減速過程,減速過程非常短暫而強(qiáng)烈,在物理上表現(xiàn)為低速流體抬升遠(yuǎn)離壁面,低速流體到達(dá)測量點(diǎn)后,使當(dāng)?shù)販y量點(diǎn)速度減少。法向脈動速度主要表現(xiàn)為劇烈的加速過程,表明流體向上加速遠(yuǎn)離壁面。在掃掠過程中,首先是一個緩慢的減速過程,然后是一個急劇的加速過程,加速過程非常短暫而強(qiáng)烈,在物理上表現(xiàn)為高速流體沖向壁面的劇烈下掃,高速流體到達(dá)測量點(diǎn)后,使當(dāng)?shù)販y量點(diǎn)速度增加。而在法向是流體向下沖向壁面,法向速度由正變?yōu)樨?fù)。

    圖6給出了在外區(qū)不同法向位置的能量最大尺度對應(yīng)的噴射和掃掠事件相干結(jié)構(gòu)的流向和法向脈動速度的條件相位平均波形。在外區(qū),噴射事件(見圖6(a))表現(xiàn)出的負(fù)的流向脈動速度仍然比較明顯,隨著遠(yuǎn)離壁面其負(fù)向速度逐漸減弱。而法向脈動速度(見圖6(c))在y+=679.3處時間尺度變大,幅值減小較慢。在y+=847.7處,時間尺度再次變大,幅值沒有減小。表明在湍流邊界層的外邊界,確實(shí)存在著超大尺度相干結(jié)構(gòu)的摻混作用。

    在外區(qū),流向脈動速度噴射事件負(fù)的流向脈動減少,流向脈動速度轉(zhuǎn)為正值。掃掠事件也是負(fù)的流向脈動減少,流向脈動速度轉(zhuǎn)為正值,條件相位波形與近壁區(qū)比較發(fā)生了改變。法向脈動速度噴射事件也呈現(xiàn)出負(fù)的脈動減少、法向脈動速度轉(zhuǎn)為正值的趨勢。掃掠事件法向脈動速度轉(zhuǎn)為正值,條件相位波形與近壁區(qū)的波形比較發(fā)生了改變。這也與外區(qū)超大尺度相干結(jié)構(gòu)大量卷吸和摻混邊界層外的高速流體有關(guān)。Adrian[19]認(rèn)為,超大尺度相干結(jié)構(gòu)的形成和發(fā)卡渦(包)有關(guān),這些發(fā)卡渦(包)在流向連貫成一串,使得低動量流體依次通過發(fā)卡渦(包)向下傳遞,這些平行的發(fā)卡渦(包)誘導(dǎo)出趨于一致的低動量區(qū)域,最終形成了超大尺度運(yùn)動。在傳遞過程中,隨著流體向下游運(yùn)動,這些超大尺度結(jié)構(gòu)逐漸遠(yuǎn)離壁面,與邊界層外環(huán)境高速流體發(fā)生動量交換和摻混,其流向脈動速度逐漸增加,從負(fù)值逐漸轉(zhuǎn)為正值,法向脈動速度緩慢地小振幅正負(fù)波動。

    從圖7中可以看出,近壁區(qū)噴射和掃掠的過程中,雷諾切應(yīng)力〈u′v′〉的條件相位波形振幅較大,均為負(fù)值,說明流向和法向脈動速度異號特征明顯,雷諾切應(yīng)力〈u′v′〉在象限分布中主要集中在二四象限。由于噴射和掃掠過程中,流向和法向速度分量的條件相位波形總保持幾乎相反的相位,因此雷諾應(yīng)力〈u′v′〉保持為負(fù)號且具有很大的幅值。隨著法向位置的增加,雷諾應(yīng)力波形的幅值逐漸減小。

    圖7 近壁區(qū)不同法向位置雷諾應(yīng)力噴射掃掠事件的條件相位平均波形Fig.7 The phase-averaging waveforms of Reynolds shear stress in the near-wall region at ejection and sweep events

    從圖8可以看出,在外區(qū),流向脈動速度在y+=847.7處能量最大尺度由第9尺度變化到第10尺度,雷諾應(yīng)力〈u′v′〉的噴射和掃掠時間尺度相應(yīng)增加,負(fù)的幅值減少,雷諾應(yīng)力表現(xiàn)為正負(fù)交替的波動現(xiàn)象。法向脈動速度在y+=679.3和y+=847.7處能量最大尺度由第8尺度變化到第9和第10尺度,雷諾應(yīng)力〈u′v′〉的噴射和掃掠時間尺度也增加,雷諾應(yīng)力表現(xiàn)為正負(fù)交替的波動現(xiàn)象。

    圖8 外區(qū)不同法向位置雷諾應(yīng)力噴射掃掠事件的條件相位平均波形Fig.8 The phase-averaging waveforms of Reynolds shear stress in the outer region at ejection and sweep events

    3 結(jié) 論

    本文用熱線風(fēng)速儀和雙絲熱線探針精細(xì)測量了風(fēng)洞中湍流邊界層不同法向位置流向和法向速度分量的時間序列信號,得出了湍流邊界層近壁區(qū)及外區(qū)不同法向位置的流向和法向脈動速度的子波分尺度能量隨尺度分布以及能量最大尺度下噴射和掃掠事件的條件相位平均波形隨法向的變化規(guī)律,得到以下結(jié)論:

    (1)在近壁區(qū)由于噴射和掃掠事件導(dǎo)致脈動劇烈,隨著法向位置遠(yuǎn)離壁面,流向和法向脈動速度的能量隨尺度分布的強(qiáng)度逐漸減弱,表明隨著法向位置遠(yuǎn)離壁面,湍動能逐漸減少。

    (2)流向脈動速度在近壁區(qū)不同法向位置能量最大尺度不變,而到了外區(qū),能量最大尺度由第9增加到第10尺度;法向脈動速度在外區(qū)的能量最大尺度由由第8尺度到第9尺度再到第10尺度逐漸增大。證實(shí)了在湍流邊界層外區(qū)出現(xiàn)了超大尺度的相干結(jié)構(gòu)。

    (3)在不同的法向位置,噴射和掃掠事件的條件相位波形發(fā)生了改變。在外區(qū),流向脈動速度噴射和掃掠事件負(fù)的脈動減少,正的脈動增多;法向脈動速度噴射和掃掠事件負(fù)的脈動減少,正的脈動增多。雷諾應(yīng)力的條件相位波形發(fā)生了改變,雷諾應(yīng)力負(fù)值減少,正值增多,表現(xiàn)為正負(fù)交替的波動特征。

    [1]Robinson S K.Coherent motions in the turbulent boundary layer[J].Annu Rev Fluid Mech,1991,23:601-639.

    [2]Falco R E.Coherent motions in the outer region of turbulent boundary layers[J].Phys Fluids,1977,20:S124-S132.

    [3]Hutchins N,Marusic I.Evidence of very long meandering features in the logarithmic region of turbulent boundary layers[J].J Fluid Mech,2007,579:1-28.

    [4]Hutchins N,Hambleton W T,Marusic I.Inclined cross-stream stereo particle image velocimetry measurements in turbulent boundary layers[J].J Fluid Mech,2005,541:21-54.

    [5]Adrian R J,Meinhart C D,Tomkins C D.Vortex organization in the outer region of the turbulent boundary layer[J].J Fluid Mech,2000,422:1-54.

    [6]Liu Z,Adrian R J,Hanratty T J.Large-scale modes of turbulent channel flow:transport and structure[J].J Fluid Mech,2001,448:53-80.

    [7]Kim K C,Adrian R J.Very large-scale motion in the outer layer[J].Phys Fluids,1999,11(2):417-422.

    [8]Christensen K T,Adrian R J.Statistical evidence of hairpin vortex packets in wall turbulence[J].J Fluid Mech,2001,431:433-443.

    [9]Natrajan V K,Christensen K T.The role of coherent structures in subgrid-scale energy transfer within the log layer of wall turbulence[J].Phys Fluids,2006,18(6):065104.

    [10]Ganapathisubramani B,Hutchins N,Hambleton W T,et al.Investigation of large-scale coherence in a turbulent boundary layer using two-point correlations[J].J Fluid Mech,2005,524:57-80.

    [11]Wu Y,Christensen K T.Population trends of spanwise vortices in wall turbulence[J].J Fluid Mech,2006,568:55-76.

    [12]Guala M,Hommema S E,Adrian R J.Large-scale and verylarge-scale motions in turbulent pipe flow[J].J Fluid Mech,2006,554:521-542.

    [13]Farge M.Wavelet transforms and their applications to turbulence[J].Annu Rev Fluid Mech,1992,24:395-457.

    [14]Marie Farge,Nicholas Kevlahan,Valerie Perrier.Wavelets and turbulence[J].Proceedings of the IEEE,1996,84(4):639-669.

    [15]姜楠,王振東,舒瑋.子波分析辨識壁湍流猝發(fā)事件的能量最大準(zhǔn)則[J].力學(xué)學(xué)報(bào),1997,29(4):406-412.

    Jiang N,Wang Z D,Shu W.The maximum energy criterion for identifying burst events in wall turbulence using wavelet analysis[J].Chinese Journal of Theoretical and Applied Mechanics,1997,29(4):406-412.

    [16]Liu J H,Jiang N,Wang Z D,et al.Multi-scale coherent structures in turbulent boundary layer detected by locally averaged velocity structure functions[J].Applied Mathematics and Mechanics,2005,26(4):456-464.

    [17]Zheng X B,Jiang N.Experimental study on spectrum and multi-scale nature of wall pressure and velocity in turbulent boundary layer[J].Chinese Physics B,2015,24(6):064702.

    [18]Del Alamo J C,Jimenez J.Spectra of the very large anisotropic scales in turbulent channels[J].Phys Fluids,2003,15(6):L41-44.

    [19]Zhou J,Adrian R J,Balachandar S,et al.Mechanism for generating coherent packets of hairpin vortices in channel flow[J].J Fluid Mech,1999,387:353-396.

    Phase-averaging waveform s of superstructures in outer layer of turbulent boundary layer

    Bai Jianxia1,2,Zheng Xiaobo1,Jiang Nan1,3,*
    (1.School of Mechanical Engineering,Tianjin University,Tianjin 300354,China;2.Department of Mechanics,Renai College of Tianjin University,Tianjin 301636,China;3.Tianjin Key Laboratory of Modern Engineering Mechanics,Tianjin 300350,China)

    An experimental investigation of the turbulent boundary layer was conducted in a wind tunnel with IFA300 constant temperature anemometer.The simultaneous time series of the streamwise velocity component and the normal velocity component at different wall-normal positions in the turbulent boundary layer were finelymeasured by a double-sensor X hot-wire probe.The distributions of the turbulence kinetic energy alongwith the scalewere obtained in the near-wall region and the outer region of the turbulent boundary layer through wavelet transform on different velocity components.In the outer region of the turbulent boundary layer,the maximum energy scale and the average burst period of the coherent structure gradually increased.Large-scale vortical structure appeared in the outer region of the turbulent boundary layer.The conditional phase-averaging waveforms of the ejection and sweep events at themaximum energy scale of the coherent structure were extracted by the conditional phase-averaging technique in the near-wall region and the outer region of the turbulent boundary layer.The waveforms of ejection and sweep events in the outer region of the turbulent boundary layer were different from the waveforms in the near-wall region.The amplitude of negative Reynolds shear stresswas reduced.The conditional phase-averaging waveforms of Reynolds shear stress presented positive and negative fluctuations.

    turbulent boundary layer;superstructures;wavelet transform;eject

    O357.5

    A

    (編輯:楊 娟)

    1672-9897(2016)05-0001-08

    10.11729/syltlx20160064

    2016-04-18;

    2016-06-16

    國家自然科學(xué)基金項(xiàng)目(11272233,11332006,11471167,11572221);國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)課題(NO.2012CB720101,NO.2012CB720103)

    *通信作者E-mail:nanj@tju.edu.cn

    Bai J X,Zheng X B,Jiang N.Phase-averaging waveforms of superstructures in outer layer of turbulent boundary layer.Journal of Experiments in Fluid Mechanics,2016,30(5):1-8.白建俠,鄭小波,姜楠.湍流邊界層外區(qū)超大尺度相干結(jié)構(gòu)相位平均波形.實(shí)驗(yàn)流體力學(xué),2016,30(5):1-8.

    白建俠(1982-),女,河北滄州人,博士研究生。研究方向:實(shí)驗(yàn)流體力學(xué)。通信地址:天津市團(tuán)泊新城博學(xué)苑天津大學(xué)仁愛學(xué)院(301636)。E-mail:jianxiabai@yeah.net

    猜你喜歡
    法向邊界層流向
    落石法向恢復(fù)系數(shù)的多因素聯(lián)合影響研究
    小溪?。×飨蜻h(yuǎn)方
    井岡教育(2020年6期)2020-12-14 03:04:42
    基于HIFiRE-2超燃發(fā)動機(jī)內(nèi)流道的激波邊界層干擾分析
    低溫狀態(tài)下的材料法向發(fā)射率測量
    十大漲幅、換手、振副、資金流向
    流向逆轉(zhuǎn)的啟示
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    落石碰撞法向恢復(fù)系數(shù)的模型試驗(yàn)研究
    非特征邊界的MHD方程的邊界層
    不透明材料波段法向發(fā)射率在線測量方法
    小蜜桃在线观看免费完整版高清| 看非洲黑人一级黄片| 国产一区亚洲一区在线观看| 午夜老司机福利剧场| 极品教师在线视频| 不卡视频在线观看欧美| 亚洲熟女精品中文字幕| 欧美性感艳星| a级毛色黄片| 精品亚洲成a人片在线观看 | 成人高潮视频无遮挡免费网站| 亚洲av中文字字幕乱码综合| 欧美少妇被猛烈插入视频| 能在线免费看毛片的网站| 中文字幕亚洲精品专区| 网址你懂的国产日韩在线| 深夜a级毛片| 舔av片在线| 国产黄色视频一区二区在线观看| 国产精品人妻久久久久久| 在线天堂最新版资源| 99久久精品热视频| 午夜福利影视在线免费观看| 老女人水多毛片| 亚洲色图av天堂| 26uuu在线亚洲综合色| 天堂8中文在线网| 99热全是精品| 亚洲在久久综合| 亚洲精品,欧美精品| 男女国产视频网站| 99热这里只有是精品在线观看| 久久久精品免费免费高清| 午夜福利在线观看免费完整高清在| 少妇人妻久久综合中文| 2021少妇久久久久久久久久久| 久久精品国产鲁丝片午夜精品| 欧美激情极品国产一区二区三区 | 国产在线一区二区三区精| 看免费成人av毛片| 日韩精品有码人妻一区| 最后的刺客免费高清国语| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 精品人妻视频免费看| 欧美一区二区亚洲| 性色av一级| 国产精品不卡视频一区二区| 国产 一区精品| 综合色丁香网| 99久久人妻综合| 亚洲色图av天堂| 亚洲欧洲国产日韩| 亚洲精品国产av成人精品| 涩涩av久久男人的天堂| 午夜福利网站1000一区二区三区| 色哟哟·www| 黄色视频在线播放观看不卡| 午夜免费观看性视频| 久久久久久九九精品二区国产| 国产国拍精品亚洲av在线观看| 在线 av 中文字幕| 欧美变态另类bdsm刘玥| 欧美日韩视频精品一区| 美女内射精品一级片tv| 欧美一级a爱片免费观看看| 亚洲第一区二区三区不卡| 美女国产视频在线观看| 亚洲国产色片| 色5月婷婷丁香| 蜜臀久久99精品久久宅男| 搡老乐熟女国产| 三级国产精品片| 亚洲无线观看免费| 视频区图区小说| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| 1000部很黄的大片| 色哟哟·www| 午夜激情福利司机影院| 蜜桃在线观看..| 亚洲精品,欧美精品| 午夜免费鲁丝| 国产乱人偷精品视频| 美女脱内裤让男人舔精品视频| 十分钟在线观看高清视频www | 日韩 亚洲 欧美在线| 亚洲人成网站在线观看播放| 观看av在线不卡| 91精品伊人久久大香线蕉| 国产乱人偷精品视频| 国产av码专区亚洲av| av又黄又爽大尺度在线免费看| 久久久色成人| 国产成人一区二区在线| www.色视频.com| 日本一二三区视频观看| 91精品一卡2卡3卡4卡| 国产精品福利在线免费观看| 好男人视频免费观看在线| 在线观看人妻少妇| 国产乱人视频| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av天美| 色网站视频免费| 丰满人妻一区二区三区视频av| 2021少妇久久久久久久久久久| 成年女人在线观看亚洲视频| 成人黄色视频免费在线看| 国产免费视频播放在线视频| 中文字幕久久专区| 久久国产乱子免费精品| 亚洲欧美清纯卡通| 国产精品不卡视频一区二区| 少妇人妻久久综合中文| 少妇熟女欧美另类| 波野结衣二区三区在线| 精品久久国产蜜桃| 国产精品蜜桃在线观看| 免费观看在线日韩| 一级二级三级毛片免费看| 欧美 日韩 精品 国产| 成人亚洲精品一区在线观看 | 国产爽快片一区二区三区| 欧美区成人在线视频| 国产男人的电影天堂91| 嘟嘟电影网在线观看| .国产精品久久| 日韩三级伦理在线观看| 亚洲欧美日韩另类电影网站 | 夫妻午夜视频| 免费看av在线观看网站| 精品国产露脸久久av麻豆| 午夜激情久久久久久久| 国产精品国产三级国产专区5o| 99国产精品免费福利视频| 亚洲欧洲国产日韩| 少妇 在线观看| 99国产精品免费福利视频| 99久久精品热视频| 日韩一区二区三区影片| 国产av码专区亚洲av| 亚洲欧洲日产国产| av国产久精品久网站免费入址| 老司机影院成人| 建设人人有责人人尽责人人享有的 | 国产精品.久久久| .国产精品久久| 日韩中字成人| 欧美另类一区| 91精品国产国语对白视频| 少妇熟女欧美另类| 99热这里只有是精品50| 五月开心婷婷网| 最黄视频免费看| 精品久久久精品久久久| 日日摸夜夜添夜夜添av毛片| 亚洲av男天堂| 男人狂女人下面高潮的视频| 美女脱内裤让男人舔精品视频| 亚洲欧洲日产国产| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 午夜福利影视在线免费观看| 尾随美女入室| 黑人高潮一二区| 亚洲精品第二区| 人妻制服诱惑在线中文字幕| 最近最新中文字幕免费大全7| a 毛片基地| 国产伦在线观看视频一区| 亚洲伊人久久精品综合| 下体分泌物呈黄色| 国产黄片美女视频| 欧美国产精品一级二级三级 | 人妻制服诱惑在线中文字幕| 欧美国产精品一级二级三级 | 国产成人精品福利久久| 黑人猛操日本美女一级片| 中文字幕亚洲精品专区| av在线蜜桃| 一级毛片久久久久久久久女| 国产成人a区在线观看| 在线观看免费日韩欧美大片 | 日本爱情动作片www.在线观看| 尤物成人国产欧美一区二区三区| 欧美精品一区二区免费开放| 亚洲国产色片| 大陆偷拍与自拍| 色哟哟·www| 久久久久久久大尺度免费视频| 精品国产露脸久久av麻豆| 成人影院久久| 国产成人精品一,二区| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频 | 亚洲,一卡二卡三卡| 一级爰片在线观看| 又大又黄又爽视频免费| 精品人妻熟女av久视频| 熟女电影av网| 91久久精品国产一区二区成人| 亚洲第一av免费看| 少妇人妻 视频| 黄色视频在线播放观看不卡| 国产av码专区亚洲av| 成人18禁高潮啪啪吃奶动态图 | 夜夜骑夜夜射夜夜干| 欧美日本视频| 国产男人的电影天堂91| 国产在线一区二区三区精| 国产乱人视频| 亚洲av免费高清在线观看| 伦精品一区二区三区| 18禁裸乳无遮挡动漫免费视频| 韩国av在线不卡| 国产高清不卡午夜福利| 天美传媒精品一区二区| 人妻 亚洲 视频| 在线观看一区二区三区| 韩国高清视频一区二区三区| 中文天堂在线官网| 日本一二三区视频观看| 人人妻人人澡人人爽人人夜夜| 国国产精品蜜臀av免费| 国产日韩欧美在线精品| 草草在线视频免费看| 亚洲人与动物交配视频| 午夜福利在线观看免费完整高清在| 日韩成人av中文字幕在线观看| 久久久久久久亚洲中文字幕| 国产精品秋霞免费鲁丝片| 国产中年淑女户外野战色| 精品国产一区二区三区久久久樱花 | 国产精品一区www在线观看| 久久久久国产网址| 免费大片18禁| 久久毛片免费看一区二区三区| 精品久久久久久电影网| 国产69精品久久久久777片| 国产成人精品久久久久久| 日本av手机在线免费观看| 一级毛片aaaaaa免费看小| 久久久久精品久久久久真实原创| 超碰97精品在线观看| 亚洲第一区二区三区不卡| 91在线精品国自产拍蜜月| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人 | 免费黄网站久久成人精品| 国产乱人视频| 中文乱码字字幕精品一区二区三区| 女人久久www免费人成看片| 免费av中文字幕在线| 亚洲精品一区蜜桃| 日韩精品有码人妻一区| 极品教师在线视频| 欧美精品一区二区免费开放| 韩国av在线不卡| 国产成人精品久久久久久| 卡戴珊不雅视频在线播放| 少妇的逼水好多| 最近最新中文字幕免费大全7| 香蕉精品网在线| 日韩一本色道免费dvd| 午夜福利影视在线免费观看| 欧美最新免费一区二区三区| 国产成人免费无遮挡视频| 另类亚洲欧美激情| 九色成人免费人妻av| 亚洲精品一区蜜桃| 欧美日韩视频精品一区| 性色av一级| 亚洲av欧美aⅴ国产| 国模一区二区三区四区视频| 一个人看视频在线观看www免费| 又粗又硬又长又爽又黄的视频| 99久久人妻综合| 成人黄色视频免费在线看| 久久精品夜色国产| 22中文网久久字幕| 中文乱码字字幕精品一区二区三区| 黄色日韩在线| 久久韩国三级中文字幕| 岛国毛片在线播放| 91在线精品国自产拍蜜月| 亚洲欧美成人精品一区二区| 国产亚洲午夜精品一区二区久久| 一个人免费看片子| av免费在线看不卡| av线在线观看网站| 最新中文字幕久久久久| 夜夜骑夜夜射夜夜干| 国产精品麻豆人妻色哟哟久久| 久久亚洲国产成人精品v| 日韩,欧美,国产一区二区三区| 日韩精品有码人妻一区| 久久久色成人| 男女免费视频国产| 免费观看的影片在线观看| 网址你懂的国产日韩在线| 欧美日韩视频高清一区二区三区二| 99久久精品一区二区三区| 一本色道久久久久久精品综合| 熟妇人妻不卡中文字幕| 欧美bdsm另类| 国产精品蜜桃在线观看| 自拍偷自拍亚洲精品老妇| 国产精品一区www在线观看| 精品熟女少妇av免费看| 99久国产av精品国产电影| 久久国产精品大桥未久av | 天堂俺去俺来也www色官网| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 伦精品一区二区三区| 国国产精品蜜臀av免费| av专区在线播放| 亚洲国产欧美在线一区| 精品久久久噜噜| 久久久成人免费电影| 午夜精品国产一区二区电影| 黄片无遮挡物在线观看| 少妇熟女欧美另类| 亚洲真实伦在线观看| 国产美女午夜福利| 国产爱豆传媒在线观看| 女性被躁到高潮视频| 国产熟女欧美一区二区| 免费观看的影片在线观看| 色网站视频免费| 亚洲成人手机| av福利片在线观看| videossex国产| 国产黄片美女视频| 97热精品久久久久久| 国产黄片美女视频| 成人无遮挡网站| .国产精品久久| 国产成人a∨麻豆精品| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 女人久久www免费人成看片| av在线播放精品| 一个人看视频在线观看www免费| 成人18禁高潮啪啪吃奶动态图 | 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人久久小说| 永久免费av网站大全| 99久国产av精品国产电影| 免费观看av网站的网址| 国内精品宾馆在线| 欧美xxⅹ黑人| 简卡轻食公司| 亚洲av中文av极速乱| 久久久久国产精品人妻一区二区| 亚洲精品日韩在线中文字幕| 久久精品国产亚洲网站| 国产亚洲5aaaaa淫片| 97精品久久久久久久久久精品| 国产爱豆传媒在线观看| 久久女婷五月综合色啪小说| 欧美xxxx黑人xx丫x性爽| av福利片在线观看| 中文字幕精品免费在线观看视频 | 狂野欧美白嫩少妇大欣赏| 日韩人妻高清精品专区| 日韩一区二区三区影片| 国产 一区精品| 婷婷色av中文字幕| 久久国产乱子免费精品| 国产欧美日韩一区二区三区在线 | 国产精品成人在线| 香蕉精品网在线| 精品久久久久久久久亚洲| 久久久久久久大尺度免费视频| 3wmmmm亚洲av在线观看| 亚洲国产精品国产精品| 日韩大片免费观看网站| 欧美一级a爱片免费观看看| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻夜夜爽99麻豆av| 精品一区二区三区视频在线| 91久久精品电影网| 日韩一区二区三区影片| 一区在线观看完整版| 色视频www国产| 国产精品国产三级专区第一集| 人妻系列 视频| 免费人妻精品一区二区三区视频| 如何舔出高潮| 国产伦精品一区二区三区视频9| 久久精品国产a三级三级三级| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲 | 伊人久久精品亚洲午夜| 六月丁香七月| 国产免费福利视频在线观看| 人妻少妇偷人精品九色| 日日啪夜夜撸| 亚洲精品中文字幕在线视频 | 国产欧美另类精品又又久久亚洲欧美| 日韩电影二区| 日韩强制内射视频| 久久精品国产亚洲av涩爱| 亚洲内射少妇av| 亚洲美女搞黄在线观看| 久久久成人免费电影| 亚洲国产精品国产精品| 日韩中字成人| 国产69精品久久久久777片| 亚洲天堂av无毛| 纯流量卡能插随身wifi吗| 亚洲欧美成人综合另类久久久| 久久精品国产鲁丝片午夜精品| 十分钟在线观看高清视频www | 内射极品少妇av片p| 最近中文字幕高清免费大全6| 伦理电影大哥的女人| 午夜老司机福利剧场| 大话2 男鬼变身卡| 99热6这里只有精品| 2022亚洲国产成人精品| 在线亚洲精品国产二区图片欧美 | 哪个播放器可以免费观看大片| 精品99又大又爽又粗少妇毛片| 菩萨蛮人人尽说江南好唐韦庄| 久久久a久久爽久久v久久| 国产成人精品婷婷| 另类亚洲欧美激情| 久久97久久精品| 老司机影院毛片| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 欧美日韩亚洲高清精品| 久久久成人免费电影| 成人美女网站在线观看视频| 亚洲精品国产色婷婷电影| 欧美人与善性xxx| 午夜福利视频精品| 国产成人精品久久久久久| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 丝瓜视频免费看黄片| 爱豆传媒免费全集在线观看| 久久毛片免费看一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲色图av天堂| 国产亚洲精品久久久com| 久久久亚洲精品成人影院| 色视频www国产| 精品国产三级普通话版| 观看美女的网站| 黄色配什么色好看| 国产亚洲精品久久久com| 久久久亚洲精品成人影院| 超碰97精品在线观看| 草草在线视频免费看| 日本欧美视频一区| 深爱激情五月婷婷| 波野结衣二区三区在线| a 毛片基地| 日本av免费视频播放| 亚洲精品乱码久久久v下载方式| 日韩中文字幕视频在线看片 | 在线观看免费高清a一片| 欧美日韩一区二区视频在线观看视频在线| 日韩免费高清中文字幕av| 国产熟女欧美一区二区| 亚洲国产最新在线播放| 不卡视频在线观看欧美| 搡老乐熟女国产| 人妻 亚洲 视频| 亚洲精品乱久久久久久| 边亲边吃奶的免费视频| 最后的刺客免费高清国语| 少妇高潮的动态图| 亚洲av福利一区| 欧美日韩国产mv在线观看视频 | 日韩中文字幕视频在线看片 | 欧美日韩视频精品一区| 高清欧美精品videossex| 九九在线视频观看精品| 91精品伊人久久大香线蕉| 岛国毛片在线播放| 啦啦啦中文免费视频观看日本| 欧美激情国产日韩精品一区| 青春草亚洲视频在线观看| 亚洲精品一二三| av一本久久久久| 国产av一区二区精品久久 | 午夜激情福利司机影院| 国内少妇人妻偷人精品xxx网站| 日本欧美国产在线视频| 欧美三级亚洲精品| 亚洲综合精品二区| 人人妻人人澡人人爽人人夜夜| 国产av国产精品国产| 免费看光身美女| 午夜免费鲁丝| 国产精品不卡视频一区二区| 欧美少妇被猛烈插入视频| 久久婷婷青草| 国产亚洲91精品色在线| 在线观看免费视频网站a站| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 看免费成人av毛片| 99热网站在线观看| 在线免费观看不下载黄p国产| 高清黄色对白视频在线免费看 | 黑人猛操日本美女一级片| 日日啪夜夜爽| 最新中文字幕久久久久| 久久久久性生活片| 久久精品人妻少妇| av国产久精品久网站免费入址| 乱系列少妇在线播放| 国产精品久久久久久久电影| 欧美日韩国产mv在线观看视频 | 久久99热这里只有精品18| 纯流量卡能插随身wifi吗| 国产精品无大码| 99热6这里只有精品| 天天躁夜夜躁狠狠久久av| 大香蕉97超碰在线| 精品国产一区二区三区久久久樱花 | 一级毛片 在线播放| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区| 国产老妇伦熟女老妇高清| 国产精品一区二区性色av| 在现免费观看毛片| 成人毛片a级毛片在线播放| 国模一区二区三区四区视频| 欧美+日韩+精品| 日本猛色少妇xxxxx猛交久久| 老熟女久久久| 久久国产乱子免费精品| 国产精品一二三区在线看| 免费观看在线日韩| 欧美成人一区二区免费高清观看| 人人妻人人爽人人添夜夜欢视频 | 久久国产精品大桥未久av | 一本—道久久a久久精品蜜桃钙片| 热re99久久精品国产66热6| 天天躁夜夜躁狠狠久久av| 亚洲精品国产av成人精品| 九草在线视频观看| 亚洲国产精品一区三区| 精品久久久噜噜| 97超视频在线观看视频| 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 亚洲精品第二区| 久久人妻熟女aⅴ| 久久6这里有精品| 人人妻人人爽人人添夜夜欢视频 | 大陆偷拍与自拍| 国模一区二区三区四区视频| 插逼视频在线观看| 日韩免费高清中文字幕av| 免费av不卡在线播放| 嫩草影院入口| 黄色日韩在线| 免费观看性生交大片5| 尾随美女入室| 欧美少妇被猛烈插入视频| 亚洲丝袜综合中文字幕| 亚洲精品色激情综合| 国产一区有黄有色的免费视频| 97超视频在线观看视频| 中文字幕久久专区| 国产一区亚洲一区在线观看| 国产一区有黄有色的免费视频| 各种免费的搞黄视频| 精品久久久噜噜| 各种免费的搞黄视频| 亚洲三级黄色毛片| 国产精品免费大片| 秋霞伦理黄片| 国产国拍精品亚洲av在线观看| 日韩一区二区视频免费看| 精品久久久久久久久av| av视频免费观看在线观看| 成人影院久久| 国产午夜精品一二区理论片| 伦理电影免费视频| 亚洲性久久影院| 欧美一区二区亚洲| 熟女av电影| 亚洲欧美日韩卡通动漫| 国产精品.久久久| 能在线免费看毛片的网站| 色网站视频免费| 美女cb高潮喷水在线观看| 天堂中文最新版在线下载| 久久久久网色| 国产高潮美女av| 97精品久久久久久久久久精品| 亚洲精品久久久久久婷婷小说| 在线免费观看不下载黄p国产| 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 精品国产三级普通话版| 精品人妻一区二区三区麻豆| 亚洲精品乱久久久久久| 国产欧美日韩精品一区二区| 激情 狠狠 欧美| 国产探花极品一区二区| 亚洲aⅴ乱码一区二区在线播放| 国产白丝娇喘喷水9色精品|