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

    運(yùn)動(dòng)單幀圖像法湍流邊界層相干結(jié)構(gòu)的實(shí)驗(yàn)測(cè)量

    2017-11-22 10:20:07曲寧寧蔡小舒周騖劉超群周雷
    化工學(xué)報(bào) 2017年11期
    關(guān)鍵詞:測(cè)量結(jié)構(gòu)

    曲寧寧,蔡小舒,周騖,劉超群,周雷

    (1上海理工大學(xué)顆粒與兩相流測(cè)量研究所,上海 200093;2 上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室(上海理工大學(xué)),上海 200093;3德州大學(xué)阿靈頓分校,美國(guó))

    運(yùn)動(dòng)單幀圖像法湍流邊界層相干結(jié)構(gòu)的實(shí)驗(yàn)測(cè)量

    曲寧寧1,2,蔡小舒1,2,周騖1,2,劉超群3,周雷1,2

    (1上海理工大學(xué)顆粒與兩相流測(cè)量研究所,上海 200093;2上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室(上海理工大學(xué)),上海 200093;3德州大學(xué)阿靈頓分校,美國(guó))

    在研究單幀單曝光圖像法(SFSE)流場(chǎng)測(cè)量的基礎(chǔ)上,提出了測(cè)量湍流邊界層相干結(jié)構(gòu)的運(yùn)動(dòng)單幀長(zhǎng)曝光圖像法(MSFLE),該方法是具有拉格朗日性質(zhì)的測(cè)量方法,測(cè)量時(shí)相機(jī)以被拍攝相干結(jié)構(gòu)的遷移速度基本相同的速度移動(dòng),同時(shí)采用長(zhǎng)曝光的方法記錄粒子的運(yùn)動(dòng)軌跡。應(yīng)用該方法,得到了充分發(fā)展湍流邊界層流向-展向平面內(nèi)相干結(jié)構(gòu)的直觀圖像,清晰地展示了相干結(jié)構(gòu)隨時(shí)間與空間的發(fā)展演變過程,并對(duì)其中的運(yùn)動(dòng)特征進(jìn)行了分析研究。實(shí)驗(yàn)結(jié)果表明,高低速條帶之間的相互作用可以導(dǎo)致發(fā)夾渦的位移、合并和耗散,鞍點(diǎn)及附近強(qiáng)剪切層的存在是相干結(jié)構(gòu)發(fā)展演化的關(guān)鍵。

    湍流;邊界層;相干結(jié)構(gòu);流動(dòng);圖像法;測(cè)量

    引 言

    湍流邊界層中相干結(jié)構(gòu)的發(fā)現(xiàn)是湍流研究領(lǐng)域的一次重大突破,相干結(jié)構(gòu)的發(fā)展和演變對(duì)湍流邊界層的壁面摩阻、傳熱傳質(zhì)和湍動(dòng)能的產(chǎn)生有重要影響[1]。20世紀(jì)60年代,Kline等[2]利用氫氣泡流動(dòng)顯示技術(shù)最先觀測(cè)到了湍流邊界層內(nèi)的條帶結(jié)構(gòu)和猝發(fā)事件,改變了人們傳統(tǒng)地認(rèn)為湍流是完全隨機(jī)流動(dòng)的觀點(diǎn)。Smits等[3]將湍流邊界層內(nèi)的主要相干結(jié)構(gòu)分為4類,包括高低速條帶[4]、發(fā)夾渦[5]、大尺度結(jié)構(gòu)和超大尺度結(jié)構(gòu)[6],其中高低速條帶與發(fā)夾渦對(duì)湍流發(fā)生、維持和發(fā)展起到關(guān)鍵作用。Adrian等[7-14]通過PIV實(shí)驗(yàn)發(fā)現(xiàn)湍流邊界層外區(qū)存在著沿流向排列的發(fā)夾渦,且具有相同的遷移速度,形成了大尺度發(fā)夾渦包結(jié)構(gòu),每個(gè)渦包由5~10個(gè)發(fā)夾渦組成,共同誘導(dǎo)產(chǎn)生了較長(zhǎng)的低速條帶。Kim等[15]認(rèn)為流向排列的大尺度結(jié)構(gòu)相互合并形成超大尺度結(jié)構(gòu),而Mckeon等[16]則認(rèn)為超大尺度結(jié)構(gòu)是由流動(dòng)的線性不穩(wěn)定性引起的。Tomkins等[17]及Ganapathisubramani等[18]對(duì)湍流邊界層流向-展向平面進(jìn)行了實(shí)驗(yàn)研究,發(fā)現(xiàn)低速條帶周圍總是出現(xiàn)反向旋轉(zhuǎn)的旋渦對(duì), 這些旋渦對(duì)被認(rèn)為是發(fā)夾渦的渦腿。王洪平等[19]利用層析PIV技術(shù)對(duì)湍流邊界層內(nèi)的相干結(jié)構(gòu)進(jìn)行了統(tǒng)計(jì)分析,發(fā)現(xiàn)在強(qiáng)剪切區(qū)域存在著周期性分布的渦結(jié)構(gòu)。潘翀等[20]通過對(duì)相干結(jié)構(gòu)的時(shí)間-空間相關(guān)分析,認(rèn)為發(fā)夾渦的對(duì)流決定了邊界層內(nèi)的流體輸運(yùn)特性。

    目前人們對(duì)于湍流邊界層內(nèi)相干結(jié)構(gòu)有兩種模型來描述,一種模型是Adrian等[7]提出的由發(fā)夾渦組成的發(fā)夾渦包結(jié)構(gòu)通過自組織及合并產(chǎn)生新的流動(dòng)結(jié)構(gòu)。另一種[21-25]是通過湍流邊界層內(nèi)的上拋、下掃和旋渦來解釋相干結(jié)構(gòu),高低速條帶的產(chǎn)生與上拋和下掃有關(guān),旋渦與剪切層的產(chǎn)生有關(guān),然而對(duì)于上拋和下掃的產(chǎn)生機(jī)理并沒有做出解釋,只是認(rèn)為這是湍流的基本屬性。

    綜上所述,相干結(jié)構(gòu)在湍流邊界層中起重要作用,人們對(duì)其特征形態(tài)有了一定的了解,但對(duì)其產(chǎn)生和發(fā)展的機(jī)理還不是很清楚,想要深入理解湍流邊界層相干結(jié)構(gòu)的本質(zhì),就必須得到相干結(jié)構(gòu)隨時(shí)間與空間的動(dòng)力學(xué)演化過程。本文利用基于單幀單曝光圖像法(single frame and single exposure imaging,SFSE)[26-29]發(fā)展而來的運(yùn)動(dòng)單幀長(zhǎng)曝光圖像法(motion single frame and long exposure imaging,MSFLE),得到了湍流邊界層相干結(jié)構(gòu)在流向-展向平面內(nèi)直觀的運(yùn)動(dòng)發(fā)展過程,并對(duì)條帶和渦結(jié)構(gòu)的相互作用關(guān)系進(jìn)行了分析研究,希望據(jù)此對(duì)湍流邊界層內(nèi)的相干結(jié)構(gòu)有更深入的理解。

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

    1.1 SFSE與MSFLE方法介紹

    在流場(chǎng)中布撒跟隨性較好的示蹤粒子,通過觀測(cè)這些粒子的運(yùn)動(dòng)軌跡,可以得到流場(chǎng)信息。單幀單曝光圖像法(SFSE)是通過設(shè)定合適的曝光時(shí)間,將流場(chǎng)中示蹤粒子的運(yùn)動(dòng)軌跡清晰地記錄在單幀圖像中,運(yùn)動(dòng)軌跡的路徑代表了流體的運(yùn)動(dòng),其長(zhǎng)度代表了曝光時(shí)間內(nèi)示蹤粒子的運(yùn)動(dòng)距離。圖1是采用SFSE圖像法得到的典型的示蹤粒子運(yùn)動(dòng)軌跡圖像及計(jì)算模型。

    圖1 典型的運(yùn)動(dòng)軌跡圖像及計(jì)算模型Fig.1 Typical motion trajectory image and calculation model

    示蹤粒子的運(yùn)動(dòng)速度可以由運(yùn)動(dòng)距離與相機(jī)曝光時(shí)間相比得出,計(jì)算公式如下

    式中,V為示蹤粒子的運(yùn)動(dòng)速度,m·s-1,S為運(yùn)動(dòng)軌跡總長(zhǎng)度,m;D為運(yùn)動(dòng)軌跡的寬度,m;M為相機(jī)鏡頭的放大倍率;Δt為相機(jī)的曝光時(shí)間,s。

    與PIV方法相比,SFSE方法有以下幾個(gè)優(yōu)點(diǎn):① PIV是通過兩次短時(shí)間間隔曝光得到的兩幅圖像來確定粒子運(yùn)動(dòng)的開始和結(jié)束位置,忽略了粒子的運(yùn)動(dòng)過程,而SFSE通過單次長(zhǎng)曝光方法從單幀圖像中直觀顯示了粒子的運(yùn)動(dòng)過程,得到的粒子運(yùn)動(dòng)信息更全面;② 在SFSE獲得的單幀圖像中,流場(chǎng)信息可以利用粒子的運(yùn)動(dòng)軌跡和曝光時(shí)間得到,避免了PIV復(fù)雜的互相關(guān)算法;③ SFSE的實(shí)驗(yàn)環(huán)境要求低,僅需較小體積的連續(xù)激光器和工業(yè)相機(jī),適于現(xiàn)場(chǎng)測(cè)量。

    在圖像法測(cè)量中,相機(jī)的視場(chǎng)和圖像的分辨率是一對(duì)矛盾。要得到高分辨率渦運(yùn)動(dòng)的精細(xì)結(jié)構(gòu),就要求相機(jī)鏡頭有較大的放大倍率,而這使得相機(jī)的視場(chǎng)較小,無法捕捉快速運(yùn)動(dòng)的流動(dòng)結(jié)構(gòu)。為解決視場(chǎng)與高分辨率間的矛盾,在SFSE方法的基礎(chǔ)上,提出了運(yùn)動(dòng)單幀長(zhǎng)曝光圖像法(MSFLE)。在該方法中,相機(jī)以被拍攝的渦的遷移速度基本相同的速度移動(dòng),采用長(zhǎng)曝光的方法記錄粒子的運(yùn)動(dòng)軌跡。

    湍流邊界層內(nèi)的相干結(jié)構(gòu)是跟隨流體運(yùn)動(dòng)的,具有拉格朗日性質(zhì),存在著沿流向的遷移速度。使用歐拉法定點(diǎn)測(cè)量時(shí),只能得到在某一時(shí)刻流體在該位置的運(yùn)動(dòng)信息,受遷移速度的影響導(dǎo)致無法觀測(cè)到相干結(jié)構(gòu),而采用拉格朗日性質(zhì)的測(cè)量方法可以跟蹤流體相干結(jié)構(gòu)。在測(cè)量時(shí)相機(jī)以一定平移速度沿流向運(yùn)動(dòng)并同時(shí)進(jìn)行實(shí)驗(yàn)記錄,當(dāng)相機(jī)運(yùn)動(dòng)速度與視野中相干結(jié)構(gòu)的遷移速度相同時(shí),不僅可以得到相干結(jié)構(gòu)的直觀圖像,而且可以顯示其隨時(shí)間與空間的運(yùn)動(dòng)發(fā)展過程,無須旋渦判別準(zhǔn)則和伽利略速度分解,可以直接分析研究。

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

    實(shí)驗(yàn)在低速循環(huán)水洞中進(jìn)行,實(shí)驗(yàn)段長(zhǎng) 250 cm、寬16 cm、高20 cm。在水洞中心水平安裝了一塊有機(jī)玻璃平板,平板尺寸1700 mm×150 mm×4 mm,正對(duì)來流方向的平板前緣設(shè)計(jì)為半橢圓形,長(zhǎng)短軸比例為 4:1,以提高平板來流的穩(wěn)定性。示蹤粒子選用平均粒徑為20 μm的空心玻璃微珠,密度為 1.05×10-3g·mm-3。使用拌線法加速邊界層轉(zhuǎn)捩,在距平板前緣100 mm處緊貼直徑為4 mm的拌線,以產(chǎn)生湍流邊界層。測(cè)量系統(tǒng)用波長(zhǎng)為 450 nm 的連續(xù)激光半導(dǎo)體激光器作光源,經(jīng)過柱透鏡輸出厚度為1 mm的片光,使用CCD相機(jī)配合片光進(jìn)行圖像記錄。整個(gè)測(cè)量系統(tǒng)安裝在一個(gè)沿流向平行布置的水平導(dǎo)軌上,移動(dòng)速率在 0~110 mm·s-1內(nèi)可調(diào)。測(cè)量系統(tǒng)布置和坐標(biāo)系定義如圖2所示,其中,x、y、z分別表示流向、法向和展向,原點(diǎn)為平板前緣的中心。

    圖2 測(cè)量系統(tǒng)和坐標(biāo)系定義Fig.2 Definition of coordinate system and measurement system

    本文實(shí)驗(yàn)中自由來流速度U∞=66 mm·s-1,相機(jī)平移速度Uc=0.83U∞,流場(chǎng)測(cè)量范圍從x=850 mm到x=1650 mm,兩個(gè)流向位置對(duì)應(yīng)的Reynolds數(shù)分別為 Reθ=θU∞/ν=429和 733,其中,θ代表動(dòng)量厚度,ν表示運(yùn)動(dòng)黏性系數(shù)。CCD相機(jī)單幀曝光時(shí)間為200 ms,在整個(gè)測(cè)量范圍內(nèi)連續(xù)曝光拍攝,采用0.14倍鏡頭,視野范圍為51 mm×39 mm,相機(jī)分辨率為1600 pixel×1200 pixel,對(duì)應(yīng)的空間分辨率為32 μm。

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

    2.1 湍流邊界層平均速度分布

    實(shí)驗(yàn)中首先采用單幀單曝光圖像法(SFSE)測(cè)量了距平板前緣850 mm與1650 mm處沿法向的邊界層平均流向速度分布。以相機(jī)側(cè)視的方式進(jìn)行流場(chǎng)測(cè)量,激光片光位于被測(cè)位置展向中心平面上,平行于來流垂直平板入射,CCD相機(jī)光軸與片光垂直進(jìn)行拍攝記錄。將測(cè)得的瞬時(shí)速度場(chǎng)進(jìn)行統(tǒng)計(jì)分析,并根據(jù)公式y(tǒng)+=yuτ/ν,u+=u/uτ將平均速度和壁面距離用摩擦速度 uτ量綱 1化,其中摩擦速度 uτ根據(jù)距離壁面1 mm內(nèi)的速度梯度獲得。得到的邊界層中量綱1平均速度分布曲線如圖3所示,該量綱1速度分布與Spalding[30]速度分布公式基本符合

    式中,K=0.4,B=5.5。從圖中可以看出兩個(gè)流向位置的平均速度分布呈現(xiàn)自相似,說明湍流邊界層在測(cè)量區(qū)域內(nèi)已經(jīng)充分發(fā)展。

    圖3 850 mm與1650 mm處邊界層平均速度分布Fig.3 Dimensionless velocity profile at positions of 850 mm and 1650 mm in boundary layer

    2.2 高低速條帶與發(fā)夾渦

    在湍流邊界層中,發(fā)夾渦和高低速條帶是重要的相干結(jié)構(gòu),其運(yùn)動(dòng)和發(fā)展過程是湍流猝發(fā)和阻力產(chǎn)生的重要原因,因此,了解其相互作用關(guān)系不僅是湍流相干結(jié)構(gòu)運(yùn)動(dòng)機(jī)理研究的核心問題之一,而且也有助于控制湍流流動(dòng)和減少壁面摩擦阻力。

    圖4和圖5是在距離壁面6 mm的流向-展向平面內(nèi)得到的相干結(jié)構(gòu)運(yùn)動(dòng)過程圖像及相干結(jié)構(gòu)關(guān)系,其中,x和z坐標(biāo)分別代表流向和展向,流體從圖片的右方流向圖片的左方,參考坐標(biāo)系即測(cè)量系統(tǒng)平移速度Uc=0.83U∞,方向自右向左。

    圖4 相干結(jié)構(gòu)的運(yùn)動(dòng)過程Fig.4 Movement process of coherent structure

    圖5 相干結(jié)構(gòu)關(guān)系Fig.5 Relations of coherent structure

    觀察連續(xù)拍攝的圖像,可以從示蹤粒子的軌跡得到在流場(chǎng)中各部分流體和旋渦相對(duì)于相機(jī)一定的速度。由于相機(jī)以旋渦的遷移速度移動(dòng),流速與相機(jī)速度相同的示蹤粒子在圖像中的軌跡基本上是個(gè)亮點(diǎn),高于遷移速度的示蹤粒子的軌跡在連續(xù)圖像中向下游移動(dòng),低于遷移速度的示蹤粒子的軌跡在連續(xù)圖像中向上游移動(dòng)。由此可以判斷流體及旋渦相對(duì)于遷移速度的運(yùn)動(dòng)快慢,從而識(shí)別高速條帶和低速條帶。

    以圖4(a)為例,此時(shí)記t1=0 s,條帶的流動(dòng)方向及旋渦的旋轉(zhuǎn)方向通過前后兩幀圖像之間的對(duì)比得出,用箭頭標(biāo)記于圖中。測(cè)量系統(tǒng)平移速度與旋渦的遷移速度相同,在整個(gè)200 ms曝光時(shí)間中,盡管旋渦已移動(dòng)了10.96 mm距離,但旋渦的圖像仍保持很清晰,表明旋渦相對(duì)相機(jī)的位置在該測(cè)量時(shí)間中基本沒有發(fā)生變化。圖中的低速條帶則相對(duì)于旋渦向上游流動(dòng),而高速條帶相對(duì)于旋渦向下游流動(dòng),高低速條帶沿展向交錯(cuò)排列。

    兩個(gè)低速條帶兩側(cè)分布著兩對(duì)反向旋轉(zhuǎn)的渦結(jié)構(gòu),逆時(shí)針旋轉(zhuǎn)的渦結(jié)構(gòu)分布在低速條帶的右側(cè),而順時(shí)針旋轉(zhuǎn)的渦結(jié)構(gòu)分布在低速條帶的左側(cè)。該兩對(duì)旋渦從三維空間考慮,應(yīng)該是穿過片光平面的發(fā)夾渦渦腿,旋渦的橢圓形狀是水平片光與傾斜的發(fā)夾渦渦腿相交的結(jié)果。A1、A2渦結(jié)構(gòu)的轉(zhuǎn)速從圖中可以判斷為約0.45 r·s-1。

    圖中還可以觀察到在低速條帶和高速條帶交界處存在著鞍點(diǎn)(stagnation point,Stag.Pnt)[7],其周圍形成強(qiáng)烈的剪切層,表現(xiàn)出很強(qiáng)的上拋(Q2)、下掃(Q4)特性,圖5給出了相干結(jié)構(gòu)關(guān)系。在經(jīng)過t1=0.5 s的發(fā)展后,如圖4(b)所示,由于上游高速條帶的沖擊,B1、B2渦結(jié)構(gòu)破裂進(jìn)而合并,而A1、A2渦結(jié)構(gòu)相對(duì)穩(wěn)定性較高,未發(fā)生合并,A1沿展向產(chǎn)生變形。在t1=1.9 s時(shí),如圖4(c)所示,可以看到在 A1的下游出現(xiàn)明顯的鞍點(diǎn),B渦結(jié)構(gòu)耗散在高速條帶中,使流動(dòng)發(fā)生彎曲。當(dāng)時(shí)間t1=2.9 s時(shí),如圖4(d)所示,由于兩側(cè)高速條帶的沖擊與擠壓,A1、A2渦結(jié)構(gòu)的相對(duì)位置沿流向位移,中間的低速條帶向上游延伸, B渦結(jié)構(gòu)在高速條帶中消亡。

    上述過程說明,不同穩(wěn)定性的渦結(jié)構(gòu)在高速條帶的沖擊下,會(huì)產(chǎn)生不同的結(jié)構(gòu)變化,穩(wěn)定性較低的渦結(jié)構(gòu)會(huì)發(fā)生合并與耗散,而穩(wěn)定性較高的渦結(jié)構(gòu)會(huì)發(fā)生位移與變形,整個(gè)運(yùn)動(dòng)過程中鞍點(diǎn)及附近強(qiáng)剪切層始終存在,說明其與相干結(jié)構(gòu)的發(fā)展演變密切相關(guān)。

    圖6記錄了高速條帶沖擊低速條帶產(chǎn)生發(fā)夾渦的過程。在圖6(a)中存在鞍點(diǎn),鞍點(diǎn)的上游是高速條帶,下游是低速條帶,高速條帶與低速條帶碰撞后分流成兩部分,包圍在低速條帶兩側(cè),在這個(gè)時(shí)刻記t2=0 s。在經(jīng)過t2=0.8 s的發(fā)展后,如圖6(b)所示,在低速條帶與左右兩側(cè)高速條帶之間形成強(qiáng)烈的剪切層,低速條帶開始破碎。在t2=1.3 s時(shí),如圖6(c)所示,在低速條帶兩側(cè)形成了完整的旋渦對(duì)-發(fā)夾渦形成。在t2=2.3 s時(shí),如圖6(d)所示,在高速條帶的推動(dòng)下,該發(fā)夾渦的運(yùn)動(dòng)速度加快,移出了相機(jī)視野,可以推測(cè)其逐漸遠(yuǎn)離壁面向下游運(yùn)動(dòng)。上述過程表明,高速條帶的沖擊可以導(dǎo)致低速條帶破碎,進(jìn)而形成發(fā)夾渦,同時(shí)鞍點(diǎn)及周圍強(qiáng)剪切層存在于整個(gè)發(fā)展過程中。

    3 結(jié) 論

    本文針對(duì)湍流邊界層相干結(jié)構(gòu)的測(cè)量提出了具有拉格朗日性質(zhì)的運(yùn)動(dòng)單幀長(zhǎng)曝光圖像法(MSFLE),并對(duì)充分發(fā)展的湍流邊界層流向-展向平面流場(chǎng)進(jìn)行了實(shí)驗(yàn)研究,測(cè)量平面距離平板 6 mm,流向測(cè)量范圍x從850 mm到1650 mm,對(duì)應(yīng)的Reynolds數(shù)分別為Reθ= 429和733。實(shí)驗(yàn)觀測(cè)了湍流邊界層流向-展向平面內(nèi)相干結(jié)構(gòu)的時(shí)空發(fā)展過程以及低速條帶破碎生成發(fā)夾渦的過程,并對(duì)其演化特征進(jìn)行了分析,得出以下結(jié)論。

    (1)基于單幀單曝光圖像法(SFSE)發(fā)展而來的運(yùn)動(dòng)單幀長(zhǎng)曝光圖像法(MSFLE)是具有拉格朗日性質(zhì)的流場(chǎng)測(cè)量方法,不需要旋渦判別準(zhǔn)則和伽利略速度分解,可以直觀地顯示出湍流邊界層內(nèi)的相干結(jié)構(gòu)形態(tài),能夠很好地實(shí)現(xiàn)對(duì)湍流邊界層流動(dòng)結(jié)構(gòu)的測(cè)量。

    圖6 發(fā)夾渦的形成過程Fig.6 Formation process of hairpin vortex

    (2)實(shí)驗(yàn)得到了湍流邊界層相干結(jié)構(gòu)包括高低速條帶、發(fā)夾渦和鞍點(diǎn)在流向-展向平面內(nèi)隨時(shí)間與空間的發(fā)展演變過程,表明利用MSFLE方法開展對(duì)湍流邊界層相干結(jié)構(gòu)的時(shí)空演化測(cè)量是一種有效的研究手段

    (3)高低速條帶之間的相互作用可以導(dǎo)致發(fā)夾渦的位移、合并和耗散,對(duì)發(fā)夾渦的發(fā)展演變有重要影響。高速條帶的沖擊可以引起低速條帶的破碎而生成發(fā)夾渦。鞍點(diǎn)及周圍剪切層的存在是湍流邊界層相干結(jié)構(gòu)運(yùn)動(dòng)演化的關(guān)鍵,為湍流邊界層的研究提供了新思路。

    符 號(hào) 說 明

    Reθ——基于動(dòng)量損失厚度的Reynolds數(shù)

    Uc——測(cè)量系統(tǒng)的平移速度,m·s-1

    U∞——自由來流速度,m·s-1

    u ——平均流向速度,m·s-1

    u+——量綱1平均流向速度

    uτ——摩擦速度,m·s-1

    y ——距離平板壁面的法向高度,m

    y+——量綱1壁面法向高度

    θ ——邊界層動(dòng)量損失厚度,m

    ν ——水的運(yùn)動(dòng)黏性系數(shù),m2·s-1

    [1] 連祺祥. 湍流邊界層相干結(jié)構(gòu)的實(shí)驗(yàn)研究[J]. 力學(xué)進(jìn)展, 2006, 3:373-388.LIAN Q X. Experimental studies of coherent structures in turbulent boundary layers[J]. Advances in Mechanics, 2006, 3: 373-388.

    [2] KLINE S J, REYNOLDS W C, SCHRAUB F A, et al. The structure of turbulent boundary layers[J]. Journal of Fluid Mechanics, 1967,30(4): 741-773.

    [3] SMITS A J, MCKEON B J, MARUSIC I. High-Reynolds number wall turbulence[J]. Annual Review of Fluid Mechanics, 2011, 43(1):353-375.

    [4] SMITS C R, METZLER S P. The characteristics of low-speed streaks in the near-wall region of a turbulent boundary layer[J]. Journal of Fluid Mechanics, 1983, 129: 27-54.

    [5] THEODORSEN T. Mechanism of turbulence[C]// Proceedings of the Second Midwestern Conference on Fluid Mechanics. Columbus:Ohio State University, 1952: 1-18.

    [6] LEE J, LEE J H, CHOI J I, et al. Spatial organization of large- and very-large-scale motions in a turbulent channel flow[J]. Journal of Fluid Mechanics, 2014, 749: 818-840

    [7] ADRIAN R J, MEINHART C D, TOMKINS C D. Vortex organization in the outer region of the turbulent boundary layer[J].Journal of Fluid Mechanics, 2000, 422(422): 1-54.

    [8] ADRIAN R J. Hairpin vortex organization in wall turbulence[J].Physics of Fluids, 2007, 19(4): 041301.

    [9] ZHOU J, ADRIAN R J, BALACHANDAR S, et al. Mechanisms for generating coherent packets of hairpin vortices in channel flow[J].Physics of Fluids, 2015, 387(10): 353-396.

    [10] KIM K, SUNG H J, ADRIAN R J. Effects of background noise on generating coherent packets of hairpin vortices[J]. Physics of Fluids,2008, 20(10): 041301.

    [11] BALAKUMAR B J, ADRIAN R J. Large- and very-large-scale motions in channel and boundary-layer flows[J]. Philosophical Transactions Mathematical Physical & Engineering Sciences, 2007,365(1852): 665-681.

    [12] ADRIAN R J, OUDHEUSDEN B W V, SCARANO F.Three-dimensional vortex organization in a high-Reynolds-number supersonic turbulent boundary layer[J]. Journal of Fluid Mechanics,2010, 644(4): 35-60.

    [13] CHEN Q, ADRIAN R J, ZHONG Q, et al. Experimental study on the role of spanwise vorticity and vortex filaments in the outer region of open-channel flow[J]. Journal of Hydraulic Research, 2014, 52(4):476-489.

    [14] TOMKINS C D, ADRIAN R J. Energetic spanwise modes in the logarithmic layer of a turbulent boundary layer[J]. Journal of Fluid Mechanics, 2005, 545: 141-162.

    [15] KIM K C, ADRIAN R J. Very large-scale motion in the outer layer[J].Physics of Fluids, 1999, 11(2): 417-422.

    [16] MCKEON B J, SHARMA A S. A critical-layer framework for turbulent pipe flow[J]. Journal of Fluid Mechanics, 2010, 658(9):336-382.

    [17] TOMKINS C D, ADRIAN R J. Spanwise structure and scale growth in turbulent boundary layers[J]. Journal of Fluid Mechanics, 2003,490(490): 37-74.

    [18] GANAPATHISUBRAMANI B, LONGMIRE E K, MARUSIC I.Characteristics of vortex packets in turbulent boundary layers[J].Journal of Fluid Mechanics, 2003, 478(478): 35-46.

    [19] 王洪平, 高琪, 王晉軍. 基于層析PIV的湍流邊界層渦結(jié)構(gòu)統(tǒng)計(jì)研究[J]. 中國(guó)科學(xué): 物理學(xué) 力學(xué) 天文學(xué), 2015, 45: 124707.WANG H P, GAO Q, WANG J J. The statistical study of vortex structure in turbulent boundary layer flow based on Tomographic PIV[J]. Sci. Sin-Phys. Mech. Astron., 2015, 45: 124707.

    [20] 潘翀, 王晉軍, 張草. 湍流邊界層 Lagrangian擬序結(jié)構(gòu)的辨識(shí)[J].中國(guó)科學(xué): 物理學(xué) 力學(xué) 天文學(xué), 2009, (4): 627-636.PAN C, WANG J J, ZHANG C. Identification of Lagrangian coherent structure of turbulent boundary layer[J]. Sci.Sin-Phys.Mech.Astron.,2009, (4): 627-636.

    [21] JIMENEZ J. Cascades in wall-bounded turbulence[J]. Annual Review of Fluid Mechanics, 2011, 44(1): 27-45.

    [22] DEL ALAMO J C, JIMENEZ J, ZANDONADE P, et al. Self-similar vortex clusters in the turbulent logarithmic region[J]. Journal of Fluid Mechanics, 2006, 561(561): 329-358.

    [23] JIMENEZ J. How linear is wall-bounded turbulence?[J]. Physics of Fluids, 2013, 25(11): 0814.

    [24] JIMENEZ J. Near-wall turbulence[J]. Physics of Fluids, 2013, 25(10):97-120.

    [25] LOZANODURAN A, JIMENEZ J, et al. Time-resolved evolution of coherent structures in turbulent channels: characterization of eddies and cascades[J]. Journal of Fluid Mechanics, 2015, 759:432-471.

    [26] 張晶晶. 基于單幀單曝光圖像法的多相流速度場(chǎng)和粒度分布測(cè)量研究[D]. 上海: 上海理工大學(xué), 2011.ZHANG J J. Study on measuring velocity field and granularity distribution of multi-phase flow based on single-frame single-exposure imaging[D]. Shanghai: University of Shanghai for Science and Technology, 2011.

    [27] 張晶晶, 范學(xué)良, 蔡小舒. 單幀單曝光圖像法測(cè)量氣固兩相流速度場(chǎng)[J]. 工程熱物理學(xué)報(bào), 2012, 33(1): 79-82.ZHANG J J, FAN X L, CAI X S. Measurement of flow field of gas-solid two-phase flow using single frame single exposure imaging[J]. Journal of Engineering Thermophysics, 2012, 33(1):79-82.

    [28] 朱贈(zèng)好, 周騖, 蔡小舒. 射流卷吸微米級(jí)結(jié)構(gòu)的圖像可視化研究[J]. 工程熱物理學(xué)報(bào), 2014, (6): 1123-1126.ZHU Z H, ZHOU W, CAI X S. Image visualization of micro-structures in the entrainment of jet flow[J]. Journal of Engineering Thermophysics, 2014, (6): 1123-1126.

    [29] 陳晶麗, 李琛, 蔡小舒, 等. 流動(dòng)多參數(shù)場(chǎng)的單幀圖像法測(cè)量方法研究[J]. 實(shí)驗(yàn)流體力學(xué), 2015, 29(6): 67-73.CHEN J L, LI C, CAI X S,et al.Research on the single frame imaging method for measuring multi-parameter fields in flow field[J].Journal of Experiments in Fluid Mechanics, 2015, 29(6): 67-73.

    [30] SPALDING D B. A single formula for the “l(fā)aw of the wall”[J].Journal of Applied Mechanics, 1961, 28(3): 455-458.

    date:2017-04-25.

    Prof. CAI Xiaoshu, usst_caixs@163.com

    supported by the National Natural Science Foundation of China (51327803, 51576130).

    Experimental measurement of coherent structures of turbulent boundary layer by motion single frame imaging method

    QU Ningning1,2, CAI Xiaoshu1,2, ZHOU Wu1,2, LIU Chaoqun3, ZHOU Lei1,2
    (1Institute of Particle and Two-phase Flow Measurement,University of Shanghai for Science and Technology,Shanghai200093,China;2Shanghai Key Laboratory of Multiphase Flow and Heat Transfer for Power Engineering in University of Shanghai for Science and Technology,Shanghai200093,China;3University of Texas at Arlington,USA)

    Based on the single frame and single exposure imaging (SFSE) method for flow field measurement, the motion single frame and long exposure imaging (MSFLE) method is proposed for measuring the coherent structure of turbulence in the plant turbulent boundary layer. The method is a Lagrangian property measurement technique. The camera moves at substantially the same speed as the migration speed of the coherent structure, and the trajectory of the particle is recorded by the method of long exposure. Experiments were carried out with this method to measure the development of the coherent structure in the turbulent boundary layer. The evolution of the coherent structure with time and space were intuitively recorded. The motion characteristics of the coherent structure are analyzed. The experimental results show that the interaction between the high and low speed streaks can lead to the displacement, merging and dissipation of the hairpin vortex. The existence of stagnation point and strong shear layer is the key to the evolution of coherent structure.

    turbulence; boundary layer; coherent structure; flow; image method; measurement

    O 357.5

    A

    0438—1157(2017)11—4088—07

    10.11949/j.issn.0438-1157.20170469

    2017-04-25收到初稿,2017-08-18收到修改稿。

    聯(lián)系人:蔡小舒。

    曲寧寧(1995—),男,碩士研究生。

    國(guó)家自然科學(xué)基金項(xiàng)目(51327803,51576130)。

    猜你喜歡
    測(cè)量結(jié)構(gòu)
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    把握四個(gè)“三” 測(cè)量變簡(jiǎn)單
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    滑動(dòng)摩擦力的測(cè)量和計(jì)算
    滑動(dòng)摩擦力的測(cè)量與計(jì)算
    測(cè)量的樂趣
    論《日出》的結(jié)構(gòu)
    測(cè)量
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
    亚洲天堂国产精品一区在线| 国产亚洲91精品色在线| 亚洲国产最新在线播放| 亚洲精品国产av蜜桃| 午夜免费男女啪啪视频观看| 肉色欧美久久久久久久蜜桃 | 亚洲国产精品国产精品| 成人亚洲欧美一区二区av| 亚洲天堂av无毛| 在线观看国产h片| 国产淫片久久久久久久久| 观看免费一级毛片| 精品国产一区二区三区久久久樱花 | 日韩av免费高清视频| 欧美日韩国产mv在线观看视频 | 国产毛片a区久久久久| 国产黄频视频在线观看| 午夜免费鲁丝| 免费看日本二区| 小蜜桃在线观看免费完整版高清| 亚洲av日韩在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩成人伦理影院| 街头女战士在线观看网站| 久久久久久国产a免费观看| 亚洲最大成人手机在线| 男人添女人高潮全过程视频| 免费观看av网站的网址| 我的女老师完整版在线观看| 嫩草影院入口| 免费黄频网站在线观看国产| 成人亚洲欧美一区二区av| 成人国产麻豆网| 少妇裸体淫交视频免费看高清| 国产精品国产三级国产专区5o| 80岁老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 欧美另类一区| 性色av一级| 国产精品不卡视频一区二区| 男女啪啪激烈高潮av片| 一级二级三级毛片免费看| 久久精品夜色国产| 久久久久久久久久成人| 国产精品无大码| 久久综合国产亚洲精品| 成人亚洲精品av一区二区| 中文资源天堂在线| av.在线天堂| 精品熟女少妇av免费看| 亚洲精品第二区| 亚洲欧美日韩另类电影网站 | 九九久久精品国产亚洲av麻豆| 精品一区二区免费观看| 亚洲自偷自拍三级| 国产黄a三级三级三级人| 日韩强制内射视频| 黄色怎么调成土黄色| 午夜亚洲福利在线播放| 日日摸夜夜添夜夜爱| 观看美女的网站| 一区二区三区乱码不卡18| 午夜激情福利司机影院| av在线亚洲专区| 日本与韩国留学比较| 偷拍熟女少妇极品色| 人妻系列 视频| 免费看不卡的av| freevideosex欧美| 亚洲色图综合在线观看| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区| 各种免费的搞黄视频| 亚洲四区av| 天堂网av新在线| 两个人的视频大全免费| 欧美一区二区亚洲| 蜜桃久久精品国产亚洲av| 三级经典国产精品| 国产精品国产三级国产av玫瑰| 男插女下体视频免费在线播放| 国产 一区 欧美 日韩| 全区人妻精品视频| 少妇的逼好多水| 人妻少妇偷人精品九色| 人妻一区二区av| 国产精品无大码| 激情 狠狠 欧美| 日本黄色片子视频| 亚洲av二区三区四区| 国产成人91sexporn| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 老司机影院毛片| 一本色道久久久久久精品综合| 最后的刺客免费高清国语| tube8黄色片| 国产精品99久久久久久久久| 纵有疾风起免费观看全集完整版| 最近最新中文字幕免费大全7| 欧美另类一区| 亚洲国产成人一精品久久久| 亚洲av二区三区四区| 99久久中文字幕三级久久日本| 亚洲av免费高清在线观看| 免费人成在线观看视频色| 日韩一本色道免费dvd| 舔av片在线| 干丝袜人妻中文字幕| 久久亚洲国产成人精品v| 男的添女的下面高潮视频| 欧美日韩一区二区视频在线观看视频在线 | 色视频在线一区二区三区| 婷婷色综合www| 晚上一个人看的免费电影| 欧美一级a爱片免费观看看| kizo精华| 白带黄色成豆腐渣| 日本免费在线观看一区| 日韩欧美精品v在线| 爱豆传媒免费全集在线观看| 欧美性猛交╳xxx乱大交人| 一区二区三区免费毛片| 国产人妻一区二区三区在| 女人被狂操c到高潮| av黄色大香蕉| 小蜜桃在线观看免费完整版高清| 久久久久久久久久人人人人人人| 中文欧美无线码| 91在线精品国自产拍蜜月| 欧美成人午夜免费资源| 久久久亚洲精品成人影院| 深爱激情五月婷婷| 免费黄网站久久成人精品| 肉色欧美久久久久久久蜜桃 | 在线播放无遮挡| 黄片wwwwww| 能在线免费看毛片的网站| 日韩电影二区| 久久精品久久久久久噜噜老黄| 一本久久精品| 秋霞伦理黄片| 亚洲欧洲国产日韩| 国产精品99久久99久久久不卡 | 国产成人a∨麻豆精品| 国产 一区 欧美 日韩| 22中文网久久字幕| 国产免费视频播放在线视频| 国产日韩欧美在线精品| av免费观看日本| 成人无遮挡网站| 大又大粗又爽又黄少妇毛片口| 亚洲av国产av综合av卡| 18禁裸乳无遮挡免费网站照片| 成年女人在线观看亚洲视频 | av在线app专区| 1000部很黄的大片| 舔av片在线| 性色av一级| 黄色日韩在线| 青春草亚洲视频在线观看| 成人亚洲精品av一区二区| 免费观看av网站的网址| 国产黄片视频在线免费观看| 最近最新中文字幕大全电影3| av网站免费在线观看视频| 欧美日韩视频高清一区二区三区二| 国产在线一区二区三区精| 成人毛片60女人毛片免费| 国产视频首页在线观看| 国产v大片淫在线免费观看| 久久久久久久久久成人| 伊人久久国产一区二区| 亚洲成人久久爱视频| 国产老妇伦熟女老妇高清| 男人狂女人下面高潮的视频| 日韩一区二区视频免费看| 成人免费观看视频高清| 亚洲精品亚洲一区二区| 日日撸夜夜添| 免费观看性生交大片5| 婷婷色综合www| 在线亚洲精品国产二区图片欧美 | 国产成人精品福利久久| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 欧美性感艳星| 国产人妻一区二区三区在| 日韩一区二区视频免费看| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 九草在线视频观看| 狂野欧美激情性xxxx在线观看| 狠狠精品人妻久久久久久综合| 日本-黄色视频高清免费观看| 久久ye,这里只有精品| 精品99又大又爽又粗少妇毛片| 人妻夜夜爽99麻豆av| 超碰av人人做人人爽久久| 涩涩av久久男人的天堂| 亚洲精品亚洲一区二区| 亚洲最大成人手机在线| 成人亚洲精品一区在线观看 | 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 亚洲精品乱久久久久久| 免费看av在线观看网站| 寂寞人妻少妇视频99o| 超碰av人人做人人爽久久| 男女那种视频在线观看| 国产一区二区在线观看日韩| 丝袜脚勾引网站| 国产欧美亚洲国产| 舔av片在线| 国产在线一区二区三区精| 亚洲一区二区三区欧美精品 | 天天躁夜夜躁狠狠久久av| 国产免费视频播放在线视频| 亚洲av中文字字幕乱码综合| 少妇人妻 视频| 亚洲欧美一区二区三区国产| 国产男女超爽视频在线观看| 欧美性猛交╳xxx乱大交人| 免费黄频网站在线观看国产| 国产真实伦视频高清在线观看| 高清欧美精品videossex| 特级一级黄色大片| 久久久久久国产a免费观看| 久久精品国产亚洲av天美| 男女边摸边吃奶| 18禁在线无遮挡免费观看视频| 午夜免费男女啪啪视频观看| 丰满人妻一区二区三区视频av| 亚洲成人一二三区av| 三级国产精品欧美在线观看| 卡戴珊不雅视频在线播放| 亚洲丝袜综合中文字幕| 99久久中文字幕三级久久日本| 热re99久久精品国产66热6| 制服丝袜香蕉在线| 丰满少妇做爰视频| 九草在线视频观看| 成人特级av手机在线观看| 新久久久久国产一级毛片| av在线观看视频网站免费| 日韩免费高清中文字幕av| 亚洲性久久影院| 一级片'在线观看视频| 亚洲人成网站在线播| 亚洲精品亚洲一区二区| 在线免费十八禁| 国产白丝娇喘喷水9色精品| 国产精品爽爽va在线观看网站| 日日啪夜夜爽| 少妇丰满av| 99热6这里只有精品| 精品国产露脸久久av麻豆| 久久久久久久精品精品| 亚洲精品一二三| 亚州av有码| 午夜精品国产一区二区电影 | 亚洲最大成人手机在线| 欧美区成人在线视频| 观看美女的网站| 99久久精品一区二区三区| 国产精品久久久久久精品电影| 国产精品国产三级国产专区5o| 日本黄大片高清| 一二三四中文在线观看免费高清| 免费av观看视频| 五月玫瑰六月丁香| 国产色爽女视频免费观看| 精品国产乱码久久久久久小说| 97精品久久久久久久久久精品| 免费观看在线日韩| 97热精品久久久久久| 麻豆成人av视频| 一边亲一边摸免费视频| 国产一区有黄有色的免费视频| 日韩中字成人| 亚洲精品中文字幕在线视频 | 午夜免费男女啪啪视频观看| 777米奇影视久久| 日本三级黄在线观看| 国产真实伦视频高清在线观看| 午夜爱爱视频在线播放| 日本黄大片高清| 欧美xxxx黑人xx丫x性爽| 天天一区二区日本电影三级| 国产精品爽爽va在线观看网站| 天堂网av新在线| 少妇人妻久久综合中文| 熟妇人妻不卡中文字幕| 国产一区二区在线观看日韩| 性色av一级| 日韩人妻高清精品专区| 亚洲欧美中文字幕日韩二区| 国产av码专区亚洲av| 一二三四中文在线观看免费高清| 色5月婷婷丁香| 麻豆国产97在线/欧美| 国产亚洲av片在线观看秒播厂| 美女主播在线视频| 欧美日韩亚洲高清精品| 亚洲精品,欧美精品| 熟女人妻精品中文字幕| 日韩欧美 国产精品| 国产成人a∨麻豆精品| 免费看光身美女| 噜噜噜噜噜久久久久久91| 国产探花在线观看一区二区| 亚州av有码| 极品少妇高潮喷水抽搐| 精品一区二区免费观看| av黄色大香蕉| 免费电影在线观看免费观看| 高清毛片免费看| 少妇人妻一区二区三区视频| 秋霞在线观看毛片| 日日啪夜夜撸| 久久久久性生活片| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 天堂中文最新版在线下载 | 日韩,欧美,国产一区二区三区| 国产一区亚洲一区在线观看| 久久久久性生活片| 日韩,欧美,国产一区二区三区| 日韩精品有码人妻一区| 精品久久久久久久末码| 午夜精品一区二区三区免费看| 日韩av在线免费看完整版不卡| av国产精品久久久久影院| 国产在线一区二区三区精| 制服丝袜香蕉在线| 超碰av人人做人人爽久久| 亚洲精品一区蜜桃| 身体一侧抽搐| 欧美性猛交╳xxx乱大交人| 亚洲精品一区蜜桃| 日韩av不卡免费在线播放| 久久97久久精品| 在线观看一区二区三区激情| 国产亚洲5aaaaa淫片| 国精品久久久久久国模美| tube8黄色片| 黄片无遮挡物在线观看| 久久人人爽人人片av| 日韩av在线免费看完整版不卡| 黄色视频在线播放观看不卡| 国内少妇人妻偷人精品xxx网站| 男人爽女人下面视频在线观看| 制服丝袜香蕉在线| 在线免费十八禁| 永久免费av网站大全| 99久国产av精品国产电影| 高清午夜精品一区二区三区| 婷婷色综合www| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 国产爱豆传媒在线观看| 久久精品夜色国产| 2021天堂中文幕一二区在线观| 韩国av在线不卡| 久久久午夜欧美精品| 亚洲精品第二区| 久久久午夜欧美精品| 免费在线观看成人毛片| 亚洲成人精品中文字幕电影| 干丝袜人妻中文字幕| 国产一区二区三区av在线| 欧美3d第一页| 亚洲第一区二区三区不卡| 亚洲精品国产av成人精品| 你懂的网址亚洲精品在线观看| 亚洲精品国产av成人精品| 肉色欧美久久久久久久蜜桃 | 免费观看av网站的网址| 精品一区二区三区视频在线| 久久久久久九九精品二区国产| 乱系列少妇在线播放| 国产永久视频网站| 久久人人爽人人片av| 久久99精品国语久久久| 一级av片app| 哪个播放器可以免费观看大片| 免费看日本二区| eeuss影院久久| 免费黄网站久久成人精品| 日本wwww免费看| 极品少妇高潮喷水抽搐| 免费不卡的大黄色大毛片视频在线观看| 精品人妻一区二区三区麻豆| 日韩三级伦理在线观看| 可以在线观看毛片的网站| 精品久久国产蜜桃| 欧美老熟妇乱子伦牲交| av国产精品久久久久影院| 五月开心婷婷网| 好男人视频免费观看在线| 99视频精品全部免费 在线| 免费观看的影片在线观看| 国产成人精品久久久久久| 国产老妇伦熟女老妇高清| 22中文网久久字幕| 成年免费大片在线观看| 看非洲黑人一级黄片| 久久精品国产亚洲av天美| 日日撸夜夜添| 精品视频人人做人人爽| 欧美性猛交╳xxx乱大交人| 亚洲在久久综合| 亚洲色图综合在线观看| 美女脱内裤让男人舔精品视频| 亚洲最大成人手机在线| 亚洲精品乱码久久久久久按摩| 国产视频内射| 免费播放大片免费观看视频在线观看| 国产 一区 欧美 日韩| 大片电影免费在线观看免费| kizo精华| 亚洲国产精品成人综合色| 黄色日韩在线| 亚洲精品日本国产第一区| 久久精品人妻少妇| 久热久热在线精品观看| 日韩制服骚丝袜av| 国产免费又黄又爽又色| a级毛片免费高清观看在线播放| 91精品一卡2卡3卡4卡| 免费观看性生交大片5| 网址你懂的国产日韩在线| 亚洲va在线va天堂va国产| 两个人的视频大全免费| 天堂俺去俺来也www色官网| 免费观看a级毛片全部| 高清毛片免费看| 爱豆传媒免费全集在线观看| 国产在线男女| 久久97久久精品| 欧美日韩综合久久久久久| 毛片女人毛片| 汤姆久久久久久久影院中文字幕| 亚洲av不卡在线观看| 99久久精品热视频| 日韩电影二区| 综合色丁香网| 精品国产一区二区三区久久久樱花 | 超碰av人人做人人爽久久| av免费观看日本| 国产精品久久久久久久久免| 一级毛片 在线播放| 亚洲美女搞黄在线观看| 久久精品国产自在天天线| 如何舔出高潮| 国产欧美日韩一区二区三区在线 | 亚洲怡红院男人天堂| 晚上一个人看的免费电影| av国产免费在线观看| 韩国高清视频一区二区三区| 亚洲av国产av综合av卡| 最近2019中文字幕mv第一页| 人妻夜夜爽99麻豆av| 国产淫片久久久久久久久| 成人毛片60女人毛片免费| 久久久色成人| 一级毛片黄色毛片免费观看视频| 六月丁香七月| 三级经典国产精品| 免费观看无遮挡的男女| 日韩国内少妇激情av| 色视频在线一区二区三区| 天天一区二区日本电影三级| 亚洲,一卡二卡三卡| 岛国毛片在线播放| 日日摸夜夜添夜夜爱| 亚洲欧美日韩东京热| 国产综合精华液| 99久久精品一区二区三区| 亚洲av福利一区| 亚洲三级黄色毛片| 18禁在线播放成人免费| 精品国产一区二区三区久久久樱花 | 亚洲av电影在线观看一区二区三区 | 国产精品爽爽va在线观看网站| 久久久久国产网址| 一本一本综合久久| 深夜a级毛片| 中文在线观看免费www的网站| 欧美国产精品一级二级三级 | 少妇的逼水好多| 欧美zozozo另类| 日本黄大片高清| 啦啦啦啦在线视频资源| 最近最新中文字幕免费大全7| 亚洲人成网站在线观看播放| 免费观看无遮挡的男女| 永久免费av网站大全| 亚洲精品456在线播放app| 99久久中文字幕三级久久日本| 另类亚洲欧美激情| 大片免费播放器 马上看| 中文天堂在线官网| 久久久久久久国产电影| 国产黄a三级三级三级人| 精品少妇黑人巨大在线播放| 日韩免费高清中文字幕av| 亚洲精品久久久久久婷婷小说| 国产久久久一区二区三区| 男女边吃奶边做爰视频| 亚洲色图av天堂| 午夜爱爱视频在线播放| 男女下面进入的视频免费午夜| 联通29元200g的流量卡| 国精品久久久久久国模美| 三级国产精品欧美在线观看| 久热久热在线精品观看| 又粗又硬又长又爽又黄的视频| 欧美+日韩+精品| 一级毛片我不卡| 22中文网久久字幕| 熟妇人妻不卡中文字幕| 97在线视频观看| 国产黄片美女视频| 免费黄色在线免费观看| 搡女人真爽免费视频火全软件| av专区在线播放| 色播亚洲综合网| 国产高潮美女av| 香蕉精品网在线| xxx大片免费视频| 日韩国内少妇激情av| 国产精品嫩草影院av在线观看| 亚洲欧美精品专区久久| 国产 一区 欧美 日韩| 一级毛片久久久久久久久女| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 在线观看三级黄色| 岛国毛片在线播放| 国产黄色视频一区二区在线观看| 99久久精品热视频| 欧美日本视频| 一本久久精品| 亚洲欧美日韩无卡精品| 精品一区二区三卡| 卡戴珊不雅视频在线播放| 久久久亚洲精品成人影院| 直男gayav资源| a级毛片免费高清观看在线播放| av福利片在线观看| 国产精品一二三区在线看| 97超视频在线观看视频| 亚洲欧美日韩东京热| 久久久久精品性色| av国产免费在线观看| 成人欧美大片| 91精品伊人久久大香线蕉| 亚洲一区二区三区欧美精品 | 亚洲最大成人av| 亚洲成人av在线免费| 性色av一级| 亚洲国产欧美在线一区| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 菩萨蛮人人尽说江南好唐韦庄| kizo精华| 国产亚洲av嫩草精品影院| 亚洲欧美日韩无卡精品| 免费人成在线观看视频色| 性插视频无遮挡在线免费观看| 天天躁夜夜躁狠狠久久av| 精品亚洲乱码少妇综合久久| 欧美性猛交╳xxx乱大交人| 欧美xxⅹ黑人| 中文字幕免费在线视频6| 97在线视频观看| 美女脱内裤让男人舔精品视频| 免费看av在线观看网站| 九色成人免费人妻av| 黄色配什么色好看| 少妇人妻精品综合一区二区| 尾随美女入室| 国产精品熟女久久久久浪| 日韩伦理黄色片| 又黄又爽又刺激的免费视频.| 男女国产视频网站| 免费大片黄手机在线观看| 日韩 亚洲 欧美在线| 国产精品一及| 久久鲁丝午夜福利片| 在线 av 中文字幕| 亚洲国产av新网站| 国产毛片在线视频| 亚洲国产精品成人久久小说| 精品一区二区免费观看| 内射极品少妇av片p| 丝袜脚勾引网站| 日日撸夜夜添| 久久鲁丝午夜福利片| 免费黄色在线免费观看| 国产精品人妻久久久影院| 黄色一级大片看看| 久久人人爽人人片av| 国产男人的电影天堂91| 伦理电影大哥的女人| 99热全是精品| 久久人人爽人人爽人人片va| 久久鲁丝午夜福利片| 亚洲欧美日韩另类电影网站 | 国产精品女同一区二区软件| 大话2 男鬼变身卡| 久久精品国产亚洲av天美|