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

    基于合成射流的二維后臺(tái)階分離流主動(dòng)控制

    2016-11-15 06:01:14李斌斌姚勇顧蘊(yùn)松程克明
    航空學(xué)報(bào) 2016年6期
    關(guān)鍵詞:附點(diǎn)射流臺(tái)階

    李斌斌, 姚勇, 顧蘊(yùn)松, 程克明

    1.西南科技大學(xué) 土木工程與建筑學(xué)院, 綿陽(yáng) 621010 2.南京航空航天大學(xué) 航空宇航學(xué)院, 南京 210016

    ?

    基于合成射流的二維后臺(tái)階分離流主動(dòng)控制

    李斌斌1,*, 姚勇1, 顧蘊(yùn)松2, 程克明2

    1.西南科技大學(xué) 土木工程與建筑學(xué)院, 綿陽(yáng)621010 2.南京航空航天大學(xué) 航空宇航學(xué)院, 南京210016

    作為一種新的流動(dòng)控制激勵(lì)器,合成射流技術(shù)在流動(dòng)分離控制、降低壓力脈動(dòng)和抑制噪聲等方面具有廣闊的應(yīng)用前景。實(shí)驗(yàn)利用合成射流主動(dòng)控制技術(shù)對(duì)二維后臺(tái)階湍流分離再附流動(dòng)控制進(jìn)行了研究,應(yīng)用表面測(cè)壓、粒子圖像測(cè)速(PIV)和熱線等多種實(shí)驗(yàn)測(cè)試技術(shù)對(duì)后臺(tái)階表面壓力分布、流場(chǎng)結(jié)構(gòu)以及剪切層特性進(jìn)行了測(cè)試。結(jié)果表明,在臺(tái)階前緣施加合成射流可有效減小回流區(qū)范圍和降低再附長(zhǎng)度,當(dāng)合成射流的動(dòng)量系數(shù)為0.301×10-3時(shí),可使再附點(diǎn)長(zhǎng)度減小25%。合成射流控制使得沿臺(tái)階下游的湍動(dòng)能和雷諾應(yīng)力增強(qiáng),提高了臺(tái)階下游流場(chǎng)的混合效率。熱線動(dòng)態(tài)結(jié)果表明頻率是后臺(tái)階分離流動(dòng)控制的關(guān)鍵參數(shù),當(dāng)頻率為260 Hz、激勵(lì)頻率與剪切層渦脫落頻率之比為1.32、激勵(lì)頻率等同于旋渦脫落頻率時(shí),合成射流控制效果最好,僅需消耗較小的能量即可實(shí)現(xiàn)流動(dòng)控制的目的。

    合成射流; 后臺(tái)階; 流動(dòng)控制; 流動(dòng)分離; 再附點(diǎn)長(zhǎng)度

    后臺(tái)階(Backward Facing Step, BFS)流場(chǎng)是典型的邊界層分離后再附之流場(chǎng),亦稱為突張室流場(chǎng)[1-2]。雖然結(jié)構(gòu)簡(jiǎn)單,但因流場(chǎng)性質(zhì)穩(wěn)定、結(jié)構(gòu)堅(jiān)固,流場(chǎng)內(nèi)流體的混合效率極佳,故常被應(yīng)用于駐焰器、燃燒室、擴(kuò)張器、機(jī)翼與建筑物周圍風(fēng)場(chǎng)等工程中[3-5]。后臺(tái)階流動(dòng)的分離將會(huì)導(dǎo)致一些不利的影響,如高速旋渦的形成、壓力損失、脈動(dòng)增大以及噪聲等[6]。因此,有必要采取流動(dòng)控制措施來(lái)抑制其分離流動(dòng)。

    后臺(tái)階的流動(dòng)分離控制可分為被動(dòng)控制和主動(dòng)控制[7]。被動(dòng)控制是指沒(méi)有外部能量的注入,如在流場(chǎng)中加入各種形狀的擾流片或渦流發(fā)生器等。利用渦流發(fā)生器在流場(chǎng)下游形成流向渦,提高流場(chǎng)的混合效率,從而降低分離區(qū)范圍??刂品椒ê?jiǎn)單易行,但存在無(wú)法根據(jù)實(shí)際工況進(jìn)行調(diào)節(jié)的缺陷。Park等[8]利用三角形擾流片形成的流向渦,使分離區(qū)減小20%以上。僅增加了流場(chǎng)部分區(qū)域的動(dòng)量交換,對(duì)流場(chǎng)中大多數(shù)區(qū)域的動(dòng)量交換降低。

    主動(dòng)控制是在流場(chǎng)中直接施加適當(dāng)?shù)臄_動(dòng)模式,并與流動(dòng)的內(nèi)在模式相耦合來(lái)實(shí)現(xiàn)對(duì)流動(dòng)的控制。在以往的研究中,吹吸氣控制被廣泛應(yīng)用于后臺(tái)階的分離流動(dòng)控制。Sano等[9]對(duì)臺(tái)階階腳處平行主流和臺(tái)階處垂直主流的后臺(tái)階分離流動(dòng)控制進(jìn)行了實(shí)驗(yàn)研究,分析了吸氣流量系數(shù)對(duì)再附點(diǎn)長(zhǎng)度和流場(chǎng)結(jié)構(gòu)的影響。陳國(guó)定和明曉[10]對(duì)不同吹吸氣位置的后臺(tái)階分離流動(dòng)控制進(jìn)行了研究,結(jié)果表明,該方法能夠有效減小再附點(diǎn)長(zhǎng)度,控制再附區(qū)附近的摩擦應(yīng)力分布。Chun和Sung[11]發(fā)現(xiàn)利用振蕩射流產(chǎn)生的周期性尾渦可以控制后臺(tái)階分離再附流動(dòng),當(dāng)上游尾渦脫落頻率St為0.2時(shí),可使再附點(diǎn)長(zhǎng)度減小10%。Dejoan和Leschziner[12]利用周期性振蕩射流對(duì)后臺(tái)階分離流動(dòng)進(jìn)行了控制研究,分析了振蕩射流幅值和頻率參數(shù)變化對(duì)后臺(tái)階流場(chǎng)結(jié)構(gòu)的影響。

    從20世紀(jì)90年代中期起,以合成射流為主導(dǎo)的零質(zhì)量射流成為主動(dòng)流動(dòng)控制研究的熱點(diǎn)[13]。與常規(guī)射流相比,合成射流具有無(wú)需氣源供應(yīng)系統(tǒng)、結(jié)構(gòu)簡(jiǎn)單、響應(yīng)快、工作頻帶寬和零質(zhì)量流率等特點(diǎn),在分離流控制、推力矢量、前體渦控制、氣動(dòng)舵面控制、旋翼流動(dòng)控制以及無(wú)人機(jī)流動(dòng)控制等方面具有廣泛的應(yīng)用[14-17]。在流動(dòng)分離控制方面,合成射流具有控制邊界層流動(dòng)分離[18-19]、推遲翼型分離、延遲失速,從而提高升阻比、改善翼型氣動(dòng)特性的效果[20-23]。Donovan等數(shù)值模擬了合成射流與定常射流對(duì)翼型分離流動(dòng)控制的效果對(duì)比,驗(yàn)證了合成射流技術(shù)在推遲邊界層分離、改善翼型氣動(dòng)特性方面具有巨大的潛能[24]。合成射流對(duì)流動(dòng)分離的自由剪切層控制研究相對(duì)較少,本文利用在臺(tái)階前緣形成的合成射流微擾動(dòng)對(duì)后臺(tái)階流動(dòng)分離形成的剪切層進(jìn)行控制研究,結(jié)合表面測(cè)壓、粒子圖像測(cè)速(PIV)和熱線動(dòng)態(tài)測(cè)試技術(shù),重點(diǎn)研究了合成射流對(duì)后臺(tái)階再附點(diǎn)長(zhǎng)度和回流區(qū)流場(chǎng)結(jié)構(gòu)的影響特性,揭示了合成射流射流能量和擾動(dòng)頻率對(duì)剪切層的控制機(jī)制。

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

    1.1射流風(fēng)洞和后臺(tái)階模型

    合成射流后臺(tái)階分離流動(dòng)控制實(shí)驗(yàn)在南京航空航天大學(xué)的二元直流式低速射流風(fēng)洞中進(jìn)行,該射流風(fēng)洞實(shí)驗(yàn)段為開口形式,出口截面為二元矩形,展向與橫向?qū)挾鹊某叽鐬?50 mm×50 mm。實(shí)驗(yàn)中主射流速度U∞=15 m/s,以臺(tái)階高度為參考的實(shí)驗(yàn)雷諾數(shù)Reh=20 548,測(cè)得射流出口中心湍流度約為3‰,射流風(fēng)洞和后臺(tái)階實(shí)驗(yàn)裝置如圖1所示。

    實(shí)驗(yàn)所用后臺(tái)階模型為典型的背向二維階梯結(jié)構(gòu),為方便PIV流場(chǎng)測(cè)試,整個(gè)模型采用透明有機(jī)玻璃材料加工。后臺(tái)階模型入口截面尺寸為250 mm×50 mm,高度H1=50 mm,與射流風(fēng)洞出口尺寸相一致。突張室高度H2=70 mm,與二元射流風(fēng)洞試驗(yàn)段出口尺寸相一致。后臺(tái)階突張室高度為70 mm,出口截面尺寸為250 mm×70 mm,臺(tái)階高度h為20 mm,模型的寬高比AR為12.5,突張比ER為1.4。

    1.2合成射流激勵(lì)器

    圖2為實(shí)驗(yàn)設(shè)計(jì)的用于后臺(tái)階分離流動(dòng)控制的合成射流激勵(lì)器。主體結(jié)構(gòu)主要由激勵(lì)器振動(dòng)腔、陣列式射流出口和揚(yáng)聲器振動(dòng)膜組成。激勵(lì)器腔體尺寸為200 mm×200 mm×80 mm,射流出口由13個(gè)直徑d=2 mm的圓孔射流出口組成,相鄰射流出口間距為10 mm,合成射流與主流方向垂直,低音揚(yáng)聲器振動(dòng)膜的最大輸出功率為40 W。

    1.3測(cè)試技術(shù)和方法

    1) 表面測(cè)壓技術(shù)

    測(cè)壓數(shù)據(jù)采集系統(tǒng)為美國(guó)NI公司的測(cè)壓系統(tǒng),主要由:PXI-6284數(shù)據(jù)采集卡,信號(hào)輸入/輸出端子板,Labview軟件和數(shù)據(jù)采集電腦組成。單通道分辨率為18 bits,最大采樣率為625 kS/s,測(cè)壓傳感器量程為0.15 psi (1 psi=6.895 kPa)。沿臺(tái)階底部中心位置下游0.5h位置,每隔0.5h臺(tái)階高度依次布置有40個(gè)直徑為l.2 mm的靜壓孔,用于實(shí)現(xiàn)對(duì)臺(tái)階下表面沿程壓力分布的測(cè)量。

    2) PIV流場(chǎng)測(cè)試技術(shù)

    采用美國(guó)TSI公司PIV測(cè)試系統(tǒng),主要由:雙脈沖Nd∶YAG激光器、互相關(guān)CCD相機(jī)、同步器、Insight軟件和示蹤粒子組成。雙脈沖激光頻率為15 Hz,測(cè)試精度為2%。圖3給出了PIV流場(chǎng)測(cè)試裝置布局圖。

    PIV測(cè)試時(shí)激光片光位置對(duì)準(zhǔn)臺(tái)階中心位置且與射流出口垂直,選用了香燃燒產(chǎn)生的煙粒子作為示蹤粒子,煙粒子濃度滿足測(cè)量要求。

    3)熱線動(dòng)態(tài)測(cè)試技術(shù)

    熱線測(cè)試系統(tǒng)為南京航空航天大學(xué)自主研制的四路熱線風(fēng)速儀,主要由:?jiǎn)谓z熱線探針(直徑為5 μm)和日本小野公司的CF920動(dòng)態(tài)信號(hào)分析儀組成。熱線探針的頻響為10 kHz,實(shí)驗(yàn)中定義fs為臺(tái)階分離剪切層的渦脫落頻率。

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

    2.1后臺(tái)階的流動(dòng)特性

    依照Bradshaw和Wong[25]的實(shí)驗(yàn)結(jié)果,當(dāng)臺(tái)階的寬高比AR>10時(shí),則可將其流場(chǎng)視為二維流動(dòng)特性,本文臺(tái)階試驗(yàn)?zāi)P偷膶捀弑華R=12.5(大于10),為了對(duì)其二維流場(chǎng)進(jìn)行驗(yàn)證,利用七孔探針測(cè)試技術(shù)對(duì)主射流沿臺(tái)階不同占位的x-y截面的x方向的主流速度Vx分布進(jìn)行了測(cè)試,測(cè)試結(jié)果見圖4。從圖中可以看出,主流沿x-y剖面的七孔探針時(shí)均速度場(chǎng)分布具有一致性,說(shuō)明臺(tái)階內(nèi)部流場(chǎng)是滿足二維流動(dòng)特性的,有利于進(jìn)行合成射流的分離流動(dòng)控制研究。

    圖5給出了臺(tái)階下壁面壓力系數(shù)Cps隨實(shí)驗(yàn)雷諾數(shù)的變化特性。隨雷諾數(shù)的變化下壁面壓力系數(shù)呈現(xiàn)如下3個(gè)典型的分布特征:① 在010的區(qū)域,表面壓力系數(shù)隨雷諾數(shù)基本無(wú)變化,說(shuō)明在該區(qū)域臺(tái)階下壁面的流動(dòng)呈現(xiàn)為再附流動(dòng)。

    圖4 七孔探針時(shí)均速度場(chǎng)測(cè)試結(jié)果Fig.4 Results of seven hole probe average velocity field

    結(jié)合圖6中的后臺(tái)階時(shí)均速度和渦量ω圖可以看出:① 在Ⅰ區(qū)08范圍內(nèi),主流流線貼附于臺(tái)階下壁面,流動(dòng)已經(jīng)發(fā)生再附,該區(qū)域內(nèi)下壁面表面壓力系數(shù)基本無(wú)變化。

    圖5 壁面壓力系數(shù)分布Fig.5 Distribution of wall pressure coefficients

    圖6 時(shí)均速度和渦量分布Fig.6 Distribution of time-average velocity and vorticity 

    2.2合成射流對(duì)再附點(diǎn)長(zhǎng)度的控制特性

    再附點(diǎn)長(zhǎng)度是衡量后臺(tái)階回流區(qū)流動(dòng)的一個(gè)重要特征參數(shù),因此本文以控制電壓U=5 V為例,重點(diǎn)研究了合成射流控制時(shí)再附點(diǎn)長(zhǎng)度隨激勵(lì)頻率的變化特性。

    圖7(a)給出了利用自制的壓力監(jiān)測(cè)探針和壓力傳感器測(cè)得的無(wú)量綱再附點(diǎn)長(zhǎng)度XL/h隨激勵(lì)頻率的變化曲線。由圖可知,無(wú)量綱再附點(diǎn)長(zhǎng)度XL/h隨控制頻率的變化較顯著,當(dāng)激勵(lì)頻率f/fs=0.71和1.32時(shí),再附點(diǎn)長(zhǎng)度XL/h達(dá)到最小,與無(wú)控制相比降低了25%。結(jié)合圖7(b)射流速度V與頻率的關(guān)系可知,當(dāng)f/fs=0.71(f=140 Hz)時(shí),合成射流速度最大,射流能量最強(qiáng);而當(dāng)f/fs=1.32(f=260 Hz)時(shí),合成射流能量較弱。關(guān)于其內(nèi)在機(jī)制將在熱線測(cè)試部分加以解釋。

    圖7 再附點(diǎn)長(zhǎng)度隨控制頻率的變化(U=5 V)Fig.7 Changes of reattachment length with forcingfrequency (U=5 V)

    本文自制的壓力監(jiān)測(cè)探針是由一對(duì)直徑?1.2 mm的不銹鋼探針管設(shè)計(jì)加工成L型,頭部采用焊錫密封成橢球形,在L型頭部正對(duì)主流前后位置開設(shè)有一個(gè)直徑?0.5 mm的前后總壓孔,構(gòu)成對(duì)臺(tái)階下壁面底層進(jìn)行監(jiān)測(cè)的壓差探針。

    圖8給出了控制電壓U=5 V、控制頻率f/fs=0.71,1.32,利用PIV測(cè)得的合成射流控制時(shí),沿臺(tái)階下游流場(chǎng)的速度V、渦量ω、湍動(dòng)能(TKE)和雷諾應(yīng)力-u'v'的分布云圖。

    圖8 合成射流控制PIV時(shí)均流場(chǎng)結(jié)果Fig.8 PIV time-average flow field results under synthetic jet control

    由速度流線圖可以看出,在x/h<3區(qū)域,流線基本與臺(tái)階位置平行,施加合成射流控制對(duì)該區(qū)域流動(dòng)的影響不顯著。在表面壓力系數(shù)受雷諾數(shù)影響較顯著的38時(shí),隨著流場(chǎng)中小渦結(jié)構(gòu)的配對(duì)、合并過(guò)程的結(jié)束,渦量逐漸減弱。施加合成射流控制后,可以看出在x/h<3剪切層的發(fā)展區(qū)域內(nèi)旋渦的強(qiáng)度減弱,流場(chǎng)中旋渦配對(duì)、合并過(guò)程的長(zhǎng)度減小。

    從湍動(dòng)能和雷諾應(yīng)力分布可以看出,在3

    2.3激勵(lì)頻率對(duì)剪切層的控制特性

    對(duì)于出口位置和大小固定的合成射流激勵(lì)器,合成射流的控制效果主要取決于激勵(lì)頻率和合成射流注入的射流動(dòng)量。那么,合成射流的控制機(jī)制是取決于激勵(lì)頻率還是射流能量,利用熱線對(duì)剪切層特性隨激勵(lì)頻率的變化進(jìn)行了監(jiān)測(cè)。

    圖9(a)給出了實(shí)驗(yàn)狀態(tài)下熱線測(cè)得的分離剪切層波形圖和功率譜結(jié)果,熱線測(cè)試位置位于臺(tái)階下游(x=1.63h,y=1.04h)位置。從圖中可以看出,分離剪切層的渦脫落頻率fs=197.5 Hz,說(shuō)明主流沿臺(tái)階邊緣位置發(fā)生流動(dòng)分離后,是以高速運(yùn)動(dòng)的類平面自由剪切層形態(tài)存在的。

    從圖9(b)~圖9(g)可以看出,當(dāng)施加控制電壓U=5 V時(shí),隨激勵(lì)頻率的增加,合成射流的射流能量逐漸增強(qiáng),在激勵(lì)頻率f/fs=0.71時(shí)射流能量達(dá)到最強(qiáng),合成射流的擾動(dòng)控制作用也逐漸增強(qiáng)。當(dāng)激勵(lì)頻率f/fs=0.71時(shí),從對(duì)應(yīng)的波形圖可以看出此時(shí)剪切層的能量較強(qiáng),但功率譜中顯示此時(shí)合成射流控制僅是使得激勵(lì)頻率下的尖峰能量增強(qiáng),而剪切層的特性并沒(méi)有發(fā)生變化。說(shuō)明剪切層對(duì)合成射流的擾動(dòng)不敏感,剪切層的能量增強(qiáng)主要是由合成射流的射流能量引起的。

    圖9 合成射流控制時(shí)熱線實(shí)驗(yàn)結(jié)果Fig.9 Hot wire experimental results under synthetic jet control 

    隨控制頻率的繼續(xù)增加,當(dāng)頻率f/fs=1.32時(shí),從波形圖可以看出,盡管此時(shí)的射流能量較弱,但此時(shí)剪切層的能量卻要強(qiáng)于f/fs=0.71的共振頻率。從對(duì)應(yīng)的功率譜可以看出,剪切層的特性發(fā)生了顯著變化,合成射流使得位于1/2倍頻下的剪切層能量增強(qiáng)。說(shuō)明在該激勵(lì)頻率下,合成射流擾動(dòng)與分離剪切層的相互作用最強(qiáng),兩種流體的頻率相當(dāng),剪切層對(duì)合成射流的擾動(dòng)起到了放大作用。較高或較低的激勵(lì)頻率,雖然合成射流的射流能量較強(qiáng),但兩種擾動(dòng)的相互干擾作用不敏感,因此不具有顯著的控制效果,在進(jìn)行流動(dòng)控制時(shí)必定會(huì)消耗更多的射流能量。

    為了進(jìn)一步加以驗(yàn)證,圖10給出了沿臺(tái)階中心底層位置,熱線測(cè)得的底層速度Vd和均方根RMS)值沿壁面下游的變化特性。從圖10(a)中可以看出,在臺(tái)階下游6h~8h范圍內(nèi)的合成射流控制影響恢復(fù)區(qū),與無(wú)控制相比,施加合成射流控制后速度得到恢復(fù),在頻率f/fs=1.32時(shí)速度變化最顯著。

    圖10 底層速度和均方根(RSM)值沿下游變化Fig.10 Downstream changes of bottom velocity and root mean square (RSM)

    從圖10(b)中的均方根值沿臺(tái)階下游的變化可知,施加合成射流控制均方根最大值位置向臺(tái)階位置移動(dòng)。當(dāng)頻率f/fs=1.32時(shí),均方根值的峰值位置由無(wú)控制時(shí)的x/h=8移至x/h=6位置,降低了25%,與自制的壓差探針監(jiān)測(cè)結(jié)果相一致。

    2.4不同強(qiáng)度合成射流的控制特性

    熱線實(shí)驗(yàn)結(jié)果表明,當(dāng)激勵(lì)器頻率等同于剪切層頻率時(shí),合成射流控制效果最好。為了分析不同射流強(qiáng)度下的控制效果,這里選取控制頻率等同于渦脫落頻率(f/fs=1.32),引入了合成射流動(dòng)量系數(shù)Cμ的概念,即

    式中:Vsj為激勵(lì)器整個(gè)工作周期內(nèi)的平均速度;S為臺(tái)階流道入口截面面積;Ssj為合成射流出口截面面積;u(t)為激勵(lì)器出口瞬時(shí)空間平均速度;T為振動(dòng)膜的振動(dòng)周期。

    圖11給出了射流動(dòng)量系數(shù)Cμ變化時(shí),合成射流對(duì)表面壓力系數(shù)分布和再附點(diǎn)長(zhǎng)度的控制效果。由圖可知,隨射流動(dòng)量系數(shù)Cμ的增加,在受雷諾數(shù)影響的3

    圖12給出了激勵(lì)頻率f/fs=1.32,熱線測(cè)得的合成射流控制時(shí)剪切層特性隨射流動(dòng)量系數(shù)Cμ變化的波形圖和功率譜信號(hào)。從熱線測(cè)試結(jié)果中可以看出,隨射流動(dòng)量系數(shù)Cμ的增加,合成射流的控制效果增強(qiáng)。功率譜結(jié)果顯示,當(dāng)最大射流動(dòng)量系數(shù)大于0.109×10-3時(shí),合成射流才對(duì)后臺(tái)階分離剪切層具有控制效果,且隨射流動(dòng)量系數(shù)的增加,合成射流對(duì)后臺(tái)階分離剪切層的擾動(dòng)不斷增強(qiáng)。

    圖11 表面壓力系數(shù)和再附點(diǎn)長(zhǎng)度隨射流動(dòng)量系數(shù)的變化Fig.11 Variation of surface pressure coefficients and reattachment length with momentum coefficient of synthetic jet

    圖12 射流動(dòng)量系數(shù)變化時(shí)熱線測(cè)試結(jié)果Fig.12 Hot wire results when synthetic jet momentum coefficient changes

    3 結(jié) 論

    1) 綜合實(shí)驗(yàn)結(jié)果表明,利用合成射流主動(dòng)控制技術(shù)可有效控制后臺(tái)階分離流動(dòng),使得再附點(diǎn)長(zhǎng)度降低,回流區(qū)范圍減小。

    2) 合成射流后臺(tái)階分離流動(dòng)控制的PIV結(jié)果表明,合成射流控制使得沿臺(tái)階下游的湍動(dòng)能和雷諾應(yīng)力增強(qiáng),提高了臺(tái)階下游流場(chǎng)的混合效率和動(dòng)量摻混,使得再附點(diǎn)長(zhǎng)度降低。

    3) 熱線動(dòng)態(tài)結(jié)果表明頻率是后臺(tái)階分離流動(dòng)控制的關(guān)鍵參數(shù)。當(dāng)激勵(lì)頻率等同于旋渦脫落頻率時(shí),合成射流的控制效果最好。此時(shí)剪切層對(duì)合成射流的擾動(dòng)具有放大作用,僅需消耗較小的能量即可實(shí)現(xiàn)流動(dòng)控制的目的。對(duì)于較高或較低的射流頻率,合成射流的控制效果均較弱,且消耗的射流能量也較多。

    4) 隨合成射流動(dòng)量系數(shù)的增加,合成射流對(duì)后臺(tái)階分離剪切層的擾動(dòng)不斷增強(qiáng),臺(tái)階下游分離區(qū)范圍逐漸減小,再附點(diǎn)長(zhǎng)度降低。當(dāng)激勵(lì)頻率f/fs=1.32時(shí),合成射流的動(dòng)量系數(shù)僅需消耗0.301×10-3,即可使再附點(diǎn)長(zhǎng)度減小25%。

    [1]EATON J K, JOHNSTOR J P. A review of research on subsonic turbulent flow reattachment[J]. AIAA Journal, 1981, 19(9): 1093-1099.

    [2]FREDERODE W R, KRADSHAW J T. Control of coherent structure of in reattaching laminar and turbulent shear layer[J]. AIAA Journal, 1986, 24(12): 1956-1963.

    [3]ARMALY B F, DURST F, PEREIRA J C F, et al. Expreimental and theoretical investigation of backward facing step flow[J]. Journal of Fluid Mechanics, 1983, 127(1): 473-496.

    [4]蘭世隆, 王晉軍. 后向臺(tái)階層流邊界層特性研究[J]. 北京航空航天大學(xué)學(xué)報(bào), 1996, 22(5): 581-584.

    LAN S L, WANG J J. The characteristics of separated shear layer a backward facing step lamnar flow[J]. Journal of Beijing University of Aeronautics and Astronautics, 1996, 22(5): 581-584 (in Chinese).

    [5]URUBA V, JONAS P, MAZUR O. Control of a channel flow behind a backward facing step by suction/blowing[J]. International Journal of Heat and Fluid Flow, 2007, 28(4): 665-672.

    [6]CHEN Y T, NIE J H, ARMALY B F, et al. Turbulent separated convection flow adjacent to backward facing step effects of step height[J]. International Journal of Heat and Mass Transfer, 2006, 49(9): 3670-3680.

    [7]KRAL L D. Active flow control technology: ASME No. FEDSM 2001-18196[R]. New York: ASME, 2001.

    [8]PARK H, JEON W P, CHOI H, et al. Mixing enhancement behind a backward facing step using tabs[J]. Physics of Fluids, 2007, 19(10): 103-105.

    [9]SANO M, FUKAZAWA K, SAKURABA K. Control of turbulent channel flow over a backward facing step by suction or injection[J]. Heat Transfer—Asian Research, 2004, 33(8): 490-504.

    [10]陳國(guó)定, 明曉. 后臺(tái)階流動(dòng)控制[C]//第十四屆研究生學(xué)術(shù)會(huì)議論文集. 南京: 南京航空航天大學(xué), 2012.

    CHEN G D, MING X. Research on the backward facing step flow control[C]//The 14th Graduate Academic Conference Proceedings. Nanjing: Nanjing University of Aeronautics and Astronautics, 2012 (in Chinese).

    [11]CHUN K B, SUNG H J. Control of turbulent separated flow over a backward facing step by local forcing[J]. Experiments in Fluids, 1996, 21(6): 417-426.

    [12]DEJOAN A, LESCHZINER M A. Large eddy simulation of periodically perturbed separated flow over a backward facing step[J]. International Journal of Heat and Fluid Flow, 2004, 25(4): 581-592.

    [13]SMITH B L, GLEZER A. The formation and evolution of synthetic jets[J]. Physics of Fluids, 1998, 10(9): 2281-2297.

    [14]CHEN F, BEELER Y, BRYANT R. Development of synthetic jet actuators for active flow control at NASA Langley: AIAA-2000-2405[R]. Reston: AIAA, 2000.

    [15]肖中云, 顧蘊(yùn)松, 江雄. 一種基于引射效應(yīng)的流體推力矢量新技術(shù)[J]. 航空學(xué)報(bào), 2012, 33(11): 1967-1974.

    XIAO Z Y, GU Y S, JIANG X. A new fluidic thrust vectoring technique based on ejecting mixing effects[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33 (11): 1967-1974 (in Chinese).

    [16]李斌斌, 姜裕標(biāo), 顧蘊(yùn)松, 等. 合成射流大攻角非對(duì)稱渦控制的試驗(yàn)研究[J]. 航空學(xué)報(bào), 2015, 36(3): 764-771.

    LI B B, JIANG Y B, GU Y S, et al. Experimental study of asymmetric vortex control at high angle of attack with synthetic jet[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(3): 764-771 (in Chinese).

    [17]CIURYLA M, LIU Y, FARNSWORTH J, et al. Flow control and flight control on a Cessna 182 model[J]. Journal of Aircraft, 2007, 44(2): 642-653.

    [18]DANDOIS J, GARNIER E. Numerical simulation of active separation control by a synthetic jet[J]. Journal of Fluid Mechanics, 2007, 574(5): 25-58.

    [19]陳占軍, 王晉軍. 合成射流改善S形進(jìn)氣道流場(chǎng)特性的研究[J]. 中國(guó)科學(xué): 技術(shù)科學(xué), 2012, 55(9): 2578-2584.

    CHEN Z J, WANG J J. Numerical investigation on sythetic jet flow control inside an S inlet duct [J]. Science China Technological Sciences, 2012, 55(9): 2578-2584 (in Chinese).

    [20]YOU D, MOIN P. Active control of flow separation over an airfoil using synthetic jets[J]. Journal of Fluids and Structures, 2008, 24(8): 1349-1357.

    [21]王林, 羅振兵, 夏智勛, 等. 合成雙射流控制翼型分離流動(dòng)的數(shù)值研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2012, 30(3): 353-357.

    WANG L, LUO Z B, XIA Z X, et al. Numerical simulation of separated flow control on an airfoil using dual sythetic jets[J]. Acta Aerodynamica Sinica, 2012, 30(3): 353-357 (in Chinese).

    [22]韓忠華, 喬志德, 宋文萍. 零質(zhì)量射流推遲翼型失速的數(shù)值模擬[J]. 航空學(xué)報(bào), 2007, 28 (5): 1040-1046.

    HAN Z H, QIAO Z D, SONG W P. Numerical simulation of active flow control to airfoil stall using local sythetic jet[J]. Acta Aeronautica et Astronautica Sinica, 2007, 28(5): 1040-1046 (in Chinese).

    [23]KIM S H, KIM C. Separation control on NACA23012 using synthetic jet[J]. Aerospace Science and Technoogy, 2009, 13(4): 172-182.

    [24]DONOVAN F, LINDA D K, CARY W. Active flow control applied to an airfoil: AIAA-1998-0210[R]. Reston: AIAA, 1998.

    [25]BRADSHAW P, WONG F Y F. The reattachment and relaxation of turbulent shear layer[J]. Journal of Fluid Mechanics, 1972, 52(1): 113-135.

    李斌斌男, 博士, 講師。主要研究方向: 流動(dòng)測(cè)試與流體流動(dòng)控制。

    Tel: 0816-2461214

    E-mail: libinbin-8@163.com

    顧蘊(yùn)松男, 教授。主要研究方向:實(shí)驗(yàn)空氣動(dòng)力學(xué), 流體流動(dòng)測(cè)試與流動(dòng)控制。

    Tel: 025-84896361

    E-mail: yunsongggu@nuaa.edu.cn

    Active control of 2D backward facing step separated flow based on synthetic jet

    LI Binbin1,*, YAO Yong1, GU Yunsong2, CHENG Keming2

    1. School of Civil Engineering and Architecture, Southwest University of Science and Technology,Mianyang621010, China 2. College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing210016, China

    As a new type of flow control actuator, synthetic jet has a potentially broad application in the fields of flow separation control, pressure pulsation reduction and noise suppression. Experimental investigation on 2D backward facing step turbulent separated and reattachment flow control with synthetic jet arrays is conducted, in which the surface pressure distribution of backward facing step, the field structure and the prominent features of shear layer are tested with many experimental devices such as pressure transducers, particle image velocimetry (PIV) and hot wire anemometer. The results show that the perturbation of synthetic jet which is formed at the upper edge of the step can effectively reduce the recirculation zone and reattachment length; when the synthetic jet momentum coefficient is 0.301×10-3, the non-dimensional length of reattachment decreases by 25% at most. Synthetic jet control increases the turbulent kinetic energy and Reynolds stress along the downstream steps and enhances the mixing efficiency of the flow field along the downstream steps. The hot wire test results show that frequency is the key parameter of backward facing step flow separation control; when the disturbance frequency is 260 Hz, the ratio of disturbance frequency to shear layer vortex shedding frequency is 1.32 and the forcing frequency is equivalent to the vortex shedding frequency, the effect of synthetic jet control is the best and the flow control can be achieved only with low consumption of energy.

    synthetic jet; backward facing step; flow control; flow separation; reattachment length

    2015-11-02; Revised: 2015-12-14; Accepted: 2016-01-08; Published online: 2016-01-1413:25

    State Key Laboratory Foundation of Aerodynamics (JBKY14010201)

    . Tel.: 0816-2461214E-mail: libinbin-8@163.com

    2015-11-02; 退修日期: 2015-12-14; 錄用日期: 2016-01-08;

    時(shí)間: 2016-01-1413:25

    www.cnki.net/kcms/detail/11.1929.V.20160114.1325.004.html

    空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室基金 (JBKY14010201)

    .Tel.: 0816-2461214E-mail: libinbin-8@163.com

    10. 7527/S1000-6893.2016.0014

    V211.7

    A

    1000-6893(2016)06-1753-10

    引用格式: 李斌斌, 姚勇, 顧蘊(yùn)松, 等. 基于合成射流的二維后臺(tái)階分離流主動(dòng)控制[J]. 航空學(xué)報(bào), 2016, 37(6): 1753-1762. LI B B, YAO Y, GU Y S, et al. Active control of 2D backward facing step separated flow based on synthetic jet[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(6): 1753-1762.

    http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

    URL: www.cnki.net/kcms/detail/11.1929.V.20160114.1325.004.html

    猜你喜歡
    附點(diǎn)射流臺(tái)階
    慣性釋放在白車身關(guān)鍵接附點(diǎn)剛度分析中的應(yīng)用
    車身接附點(diǎn)動(dòng)剛度后處理方法對(duì)比
    深海逃逸艙射流注水均壓過(guò)程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    基于逆子結(jié)構(gòu)分析方法的車門振動(dòng)噪聲分析
    論節(jié)奏中特殊時(shí)值的簡(jiǎn)快教學(xué)方法
    走在除法的臺(tái)階上
    臺(tái)階
    77級(jí)臺(tái)階
    射流齒形噴嘴射流流場(chǎng)與氣動(dòng)聲學(xué)分析
    日本免费a在线| 看片在线看免费视频| 亚洲av五月六月丁香网| 亚洲成人精品中文字幕电影| 成年女人永久免费观看视频| 久久性视频一级片| 久久国产精品人妻蜜桃| 国产99白浆流出| 一进一出好大好爽视频| 国模一区二区三区四区视频 | 精品国产亚洲在线| 亚洲精品色激情综合| 亚洲国产精品成人综合色| 欧美丝袜亚洲另类 | 一个人免费在线观看电影 | 国产主播在线观看一区二区| 亚洲国产欧美一区二区综合| 色av中文字幕| 国语自产精品视频在线第100页| 亚洲成人中文字幕在线播放| 啪啪无遮挡十八禁网站| 特大巨黑吊av在线直播| 老司机在亚洲福利影院| 国产高清激情床上av| 日日干狠狠操夜夜爽| 国产97色在线日韩免费| 亚洲欧美日韩东京热| 精品人妻1区二区| 亚洲av日韩精品久久久久久密| av女优亚洲男人天堂 | 免费看美女性在线毛片视频| 国产成人欧美在线观看| 国产毛片a区久久久久| 少妇人妻一区二区三区视频| 亚洲欧美一区二区三区黑人| 国产精品 欧美亚洲| 激情在线观看视频在线高清| 91麻豆精品激情在线观看国产| 亚洲五月天丁香| 亚洲国产中文字幕在线视频| 日本成人三级电影网站| 国产精品一区二区三区四区久久| 精品熟女少妇八av免费久了| 国产亚洲精品综合一区在线观看| 97人妻精品一区二区三区麻豆| 天天添夜夜摸| 在线免费观看不下载黄p国产 | 国产高清有码在线观看视频| 亚洲欧美精品综合久久99| 国产成+人综合+亚洲专区| 18禁观看日本| 久久中文字幕一级| 日韩欧美在线二视频| h日本视频在线播放| avwww免费| 亚洲国产精品成人综合色| 精品99又大又爽又粗少妇毛片 | 一本一本综合久久| 99国产精品99久久久久| 国产三级中文精品| 一个人免费在线观看电影 | 99热只有精品国产| 可以在线观看的亚洲视频| 国产激情欧美一区二区| av中文乱码字幕在线| 欧美性猛交╳xxx乱大交人| 两人在一起打扑克的视频| 国产熟女xx| 日本五十路高清| 欧美性猛交╳xxx乱大交人| 国产69精品久久久久777片 | 国产美女午夜福利| 又黄又粗又硬又大视频| 亚洲成人中文字幕在线播放| 18禁裸乳无遮挡免费网站照片| 中文在线观看免费www的网站| 一二三四在线观看免费中文在| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 高潮久久久久久久久久久不卡| 麻豆久久精品国产亚洲av| 国产欧美日韩精品亚洲av| 真人一进一出gif抽搐免费| 日韩大尺度精品在线看网址| 神马国产精品三级电影在线观看| 日本免费a在线| 久久久久国产一级毛片高清牌| 不卡av一区二区三区| 亚洲片人在线观看| 亚洲九九香蕉| 熟女人妻精品中文字幕| 国产精品一区二区免费欧美| 搡老熟女国产l中国老女人| 免费av不卡在线播放| 亚洲成人久久性| 男女下面进入的视频免费午夜| 国产亚洲欧美在线一区二区| 国产精华一区二区三区| aaaaa片日本免费| 国产精品永久免费网站| 99久久99久久久精品蜜桃| 国产精品乱码一区二三区的特点| 成年女人看的毛片在线观看| 18禁黄网站禁片免费观看直播| 国产av不卡久久| 桃红色精品国产亚洲av| 国产av麻豆久久久久久久| 精品电影一区二区在线| 国产精品女同一区二区软件 | 久久久久久久久久黄片| 可以在线观看的亚洲视频| 欧美高清成人免费视频www| 亚洲乱码一区二区免费版| 一个人免费在线观看电影 | 99在线人妻在线中文字幕| 亚洲成人久久爱视频| 亚洲国产欧美网| 日韩成人在线观看一区二区三区| 最近最新中文字幕大全免费视频| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 一进一出抽搐gif免费好疼| 黄色女人牲交| 黄色片一级片一级黄色片| 国产伦一二天堂av在线观看| 在线免费观看的www视频| av女优亚洲男人天堂 | 久久天躁狠狠躁夜夜2o2o| 国产 一区 欧美 日韩| 老汉色∧v一级毛片| 90打野战视频偷拍视频| 成人18禁在线播放| 又粗又爽又猛毛片免费看| 性色av乱码一区二区三区2| 国产成人啪精品午夜网站| 全区人妻精品视频| 天天添夜夜摸| 久久热在线av| 又粗又爽又猛毛片免费看| 51午夜福利影视在线观看| 国产成人啪精品午夜网站| 国产午夜福利久久久久久| 免费高清视频大片| 中文字幕精品亚洲无线码一区| 99久久精品一区二区三区| 三级国产精品欧美在线观看 | 天堂网av新在线| 全区人妻精品视频| 国产成人系列免费观看| 日韩三级视频一区二区三区| 免费搜索国产男女视频| 国产97色在线日韩免费| 久久久久久久久中文| 日韩中文字幕欧美一区二区| 在线观看一区二区三区| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| av黄色大香蕉| 不卡av一区二区三区| 此物有八面人人有两片| 日韩欧美三级三区| 久久精品91无色码中文字幕| 国产v大片淫在线免费观看| 99国产精品99久久久久| 国产不卡一卡二| 俺也久久电影网| 亚洲av中文字字幕乱码综合| 亚洲片人在线观看| 三级毛片av免费| 亚洲最大成人中文| 国产一区二区在线av高清观看| 两人在一起打扑克的视频| 白带黄色成豆腐渣| 淫秽高清视频在线观看| 高清在线国产一区| 久久人人精品亚洲av| 国产精品久久久久久人妻精品电影| 久久国产精品人妻蜜桃| 欧美在线黄色| 精品一区二区三区视频在线 | 18禁黄网站禁片免费观看直播| 午夜免费激情av| 中文字幕av在线有码专区| www日本黄色视频网| 给我免费播放毛片高清在线观看| 天天一区二区日本电影三级| 亚洲精品国产精品久久久不卡| 久久久久国产一级毛片高清牌| 动漫黄色视频在线观看| 岛国在线免费视频观看| 看片在线看免费视频| 五月伊人婷婷丁香| 成年人黄色毛片网站| 久久精品国产综合久久久| 精品国产三级普通话版| x7x7x7水蜜桃| 1000部很黄的大片| 日本 欧美在线| 可以在线观看毛片的网站| 美女大奶头视频| 午夜福利视频1000在线观看| 手机成人av网站| 免费高清视频大片| 色尼玛亚洲综合影院| 不卡一级毛片| 国产精品日韩av在线免费观看| 好男人在线观看高清免费视频| 高潮久久久久久久久久久不卡| 两个人的视频大全免费| 亚洲中文av在线| 成年免费大片在线观看| 免费看光身美女| 国产av在哪里看| 国产精品自产拍在线观看55亚洲| 禁无遮挡网站| 久久久久精品国产欧美久久久| 免费无遮挡裸体视频| 一区福利在线观看| 欧美日韩乱码在线| 欧美丝袜亚洲另类 | 黄色片一级片一级黄色片| 欧美最黄视频在线播放免费| 久久久久久人人人人人| 成年女人永久免费观看视频| 又爽又黄无遮挡网站| 久久久久国产精品人妻aⅴ院| 精品国产三级普通话版| 男插女下体视频免费在线播放| 中文字幕最新亚洲高清| 亚洲国产高清在线一区二区三| 宅男免费午夜| 久久九九热精品免费| 最新美女视频免费是黄的| 国产亚洲av高清不卡| 一进一出好大好爽视频| 禁无遮挡网站| 99热精品在线国产| 伊人久久大香线蕉亚洲五| 综合色av麻豆| 欧美丝袜亚洲另类 | 波多野结衣高清无吗| 久久久精品欧美日韩精品| 亚洲精品456在线播放app | 国产一区二区三区视频了| 精品乱码久久久久久99久播| www.精华液| 一进一出抽搐gif免费好疼| 偷拍熟女少妇极品色| 99热6这里只有精品| 又黄又粗又硬又大视频| 日本与韩国留学比较| 亚洲av第一区精品v没综合| 国产v大片淫在线免费观看| 免费无遮挡裸体视频| 久久精品国产清高在天天线| 青草久久国产| 夜夜躁狠狠躁天天躁| 成人av一区二区三区在线看| 久久久久久久午夜电影| 精品无人区乱码1区二区| 欧美日本视频| 五月伊人婷婷丁香| 男人的好看免费观看在线视频| 极品教师在线免费播放| 欧美黄色淫秽网站| 婷婷亚洲欧美| 岛国视频午夜一区免费看| 男人和女人高潮做爰伦理| 国产欧美日韩精品一区二区| 精品一区二区三区视频在线观看免费| 又黄又爽又免费观看的视频| 美女午夜性视频免费| 高清毛片免费观看视频网站| 成人欧美大片| 91av网一区二区| 午夜福利在线观看吧| 窝窝影院91人妻| 亚洲成人中文字幕在线播放| 一边摸一边抽搐一进一小说| 免费在线观看日本一区| 国内少妇人妻偷人精品xxx网站 | 色尼玛亚洲综合影院| 亚洲中文字幕一区二区三区有码在线看 | 18禁黄网站禁片免费观看直播| a级毛片a级免费在线| 美女cb高潮喷水在线观看 | 午夜福利欧美成人| 午夜福利免费观看在线| 99久久国产精品久久久| 91在线观看av| 久久精品国产清高在天天线| 三级毛片av免费| 少妇熟女aⅴ在线视频| 精品国产三级普通话版| 日韩高清综合在线| 国产欧美日韩精品一区二区| 国产三级在线视频| 日韩av在线大香蕉| 精品久久久久久久久久免费视频| 国产成人系列免费观看| 国产成人影院久久av| 丁香六月欧美| 国产真人三级小视频在线观看| 日本一二三区视频观看| 久久久久久久久中文| 男人舔女人下体高潮全视频| 黄片大片在线免费观看| 看黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 天天躁日日操中文字幕| 国产欧美日韩精品一区二区| 巨乳人妻的诱惑在线观看| 非洲黑人性xxxx精品又粗又长| 1000部很黄的大片| 国产亚洲精品久久久com| 成年版毛片免费区| 母亲3免费完整高清在线观看| 成人永久免费在线观看视频| 俄罗斯特黄特色一大片| 黑人操中国人逼视频| 亚洲国产欧美人成| 看片在线看免费视频| 长腿黑丝高跟| 可以在线观看毛片的网站| 99精品欧美一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| xxxwww97欧美| 亚洲av成人不卡在线观看播放网| 国产成人av教育| 身体一侧抽搐| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 亚洲欧美日韩高清专用| 日本黄大片高清| 男插女下体视频免费在线播放| 亚洲色图av天堂| а√天堂www在线а√下载| 好看av亚洲va欧美ⅴa在| 亚洲无线在线观看| netflix在线观看网站| 久久精品综合一区二区三区| 麻豆一二三区av精品| 国产成人福利小说| 国产激情偷乱视频一区二区| 国产一区二区在线av高清观看| 免费观看精品视频网站| 亚洲国产欧美人成| 国产亚洲精品综合一区在线观看| 极品教师在线免费播放| 超碰成人久久| 亚洲精品中文字幕一二三四区| 欧美日韩国产亚洲二区| 亚洲av第一区精品v没综合| 成人永久免费在线观看视频| а√天堂www在线а√下载| 老司机深夜福利视频在线观看| 国产精品女同一区二区软件 | 五月伊人婷婷丁香| 这个男人来自地球电影免费观看| 小蜜桃在线观看免费完整版高清| 搡老岳熟女国产| 青草久久国产| 亚洲自偷自拍图片 自拍| 亚洲av免费在线观看| 美女高潮的动态| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸| 久久精品国产99精品国产亚洲性色| 日韩中文字幕欧美一区二区| 制服丝袜大香蕉在线| 啦啦啦观看免费观看视频高清| 日韩欧美精品v在线| av女优亚洲男人天堂 | 黄频高清免费视频| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 后天国语完整版免费观看| 在线十欧美十亚洲十日本专区| 90打野战视频偷拍视频| 成人欧美大片| 波多野结衣高清作品| 久久人妻av系列| 国产精品av视频在线免费观看| 欧美色欧美亚洲另类二区| 国产精品亚洲一级av第二区| 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 国产成人av激情在线播放| 99精品在免费线老司机午夜| 国产高清激情床上av| 97超级碰碰碰精品色视频在线观看| 精品一区二区三区视频在线观看免费| 亚洲av熟女| 特级一级黄色大片| 女生性感内裤真人,穿戴方法视频| 高清在线国产一区| 久久国产乱子伦精品免费另类| a在线观看视频网站| 一本综合久久免费| 在线免费观看的www视频| 国产v大片淫在线免费观看| 黄色 视频免费看| 久久亚洲真实| 国产成人精品久久二区二区免费| 国产蜜桃级精品一区二区三区| 国内精品久久久久久久电影| 一区二区三区激情视频| 日本a在线网址| 免费一级毛片在线播放高清视频| 男人舔女人下体高潮全视频| 最近在线观看免费完整版| 成年版毛片免费区| 亚洲中文av在线| 观看免费一级毛片| 国产又色又爽无遮挡免费看| 精品一区二区三区四区五区乱码| 色哟哟哟哟哟哟| 天天躁日日操中文字幕| 熟女电影av网| 精品一区二区三区视频在线 | 久久精品国产亚洲av香蕉五月| 中文在线观看免费www的网站| 88av欧美| 欧美高清成人免费视频www| 在线观看美女被高潮喷水网站 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲av美国av| 人人妻人人看人人澡| 少妇裸体淫交视频免费看高清| 日韩欧美在线二视频| 99re在线观看精品视频| 制服人妻中文乱码| 欧美一区二区国产精品久久精品| 超碰成人久久| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 国产成人一区二区三区免费视频网站| 91av网一区二区| 亚洲第一电影网av| 99在线人妻在线中文字幕| 欧美高清成人免费视频www| 51午夜福利影视在线观看| 亚洲精品一卡2卡三卡4卡5卡| 成年人黄色毛片网站| 成人国产一区最新在线观看| 老司机福利观看| 1024香蕉在线观看| 看免费av毛片| 欧美性猛交黑人性爽| 免费在线观看视频国产中文字幕亚洲| cao死你这个sao货| 亚洲一区二区三区色噜噜| 午夜激情福利司机影院| 男人舔女人下体高潮全视频| 夜夜躁狠狠躁天天躁| or卡值多少钱| 国语自产精品视频在线第100页| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| www.精华液| 中文字幕人妻丝袜一区二区| 好看av亚洲va欧美ⅴa在| 99久久无色码亚洲精品果冻| 999久久久精品免费观看国产| 久久这里只有精品19| 啦啦啦免费观看视频1| 制服人妻中文乱码| 精品久久久久久久毛片微露脸| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| 男人舔奶头视频| 欧美极品一区二区三区四区| 亚洲av美国av| 香蕉久久夜色| 99riav亚洲国产免费| 男人舔女人的私密视频| 欧美极品一区二区三区四区| 一二三四社区在线视频社区8| 国产亚洲精品av在线| 91久久精品国产一区二区成人 | 亚洲自拍偷在线| 法律面前人人平等表现在哪些方面| 中文字幕人妻丝袜一区二区| 久久精品亚洲精品国产色婷小说| 99riav亚洲国产免费| 欧美+亚洲+日韩+国产| 人妻久久中文字幕网| 国产成年人精品一区二区| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久久久久久久久| 19禁男女啪啪无遮挡网站| 好看av亚洲va欧美ⅴa在| 久久天躁狠狠躁夜夜2o2o| 午夜福利高清视频| 欧美另类亚洲清纯唯美| 哪里可以看免费的av片| 欧美精品啪啪一区二区三区| 真人一进一出gif抽搐免费| 日韩成人在线观看一区二区三区| 久久国产精品人妻蜜桃| 偷拍熟女少妇极品色| 日本成人三级电影网站| 精品无人区乱码1区二区| 国产高清激情床上av| 成人欧美大片| 免费人成视频x8x8入口观看| 老汉色∧v一级毛片| 又爽又黄无遮挡网站| 日韩欧美三级三区| 成人永久免费在线观看视频| 男插女下体视频免费在线播放| 国产av在哪里看| 久久久久精品国产欧美久久久| 视频区欧美日本亚洲| 曰老女人黄片| 可以在线观看的亚洲视频| 国产不卡一卡二| 99热精品在线国产| 99国产精品一区二区蜜桃av| 午夜福利高清视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩中文字幕国产精品一区二区三区| 免费无遮挡裸体视频| 精品熟女少妇八av免费久了| 男人舔女人下体高潮全视频| 丁香六月欧美| 中亚洲国语对白在线视频| 97碰自拍视频| 久久久久精品国产欧美久久久| a级毛片a级免费在线| 久久国产精品影院| 亚洲成a人片在线一区二区| 一进一出好大好爽视频| 久久久久免费精品人妻一区二区| 桃红色精品国产亚洲av| 日韩欧美精品v在线| 亚洲av成人精品一区久久| 偷拍熟女少妇极品色| 欧美三级亚洲精品| 成年女人毛片免费观看观看9| 97超视频在线观看视频| 一个人看视频在线观看www免费 | cao死你这个sao货| 亚洲专区中文字幕在线| 国内久久婷婷六月综合欲色啪| 午夜福利免费观看在线| 在线观看免费视频日本深夜| 免费在线观看成人毛片| 久久天堂一区二区三区四区| 国产伦在线观看视频一区| 香蕉国产在线看| 欧美在线黄色| 国产激情偷乱视频一区二区| 麻豆av在线久日| 久久久久亚洲av毛片大全| 国产又黄又爽又无遮挡在线| 国产精品av久久久久免费| 中文亚洲av片在线观看爽| 性欧美人与动物交配| 长腿黑丝高跟| 网址你懂的国产日韩在线| 热99re8久久精品国产| 欧美午夜高清在线| 欧美高清成人免费视频www| 嫁个100分男人电影在线观看| 在线视频色国产色| 免费观看精品视频网站| 一个人观看的视频www高清免费观看 | 久久久成人免费电影| 99久久99久久久精品蜜桃| 嫁个100分男人电影在线观看| 亚洲国产精品999在线| 啦啦啦韩国在线观看视频| 天天躁日日操中文字幕| 人人妻人人澡欧美一区二区| 久久久精品大字幕| 国产精品久久久久久人妻精品电影| 嫩草影视91久久| 观看免费一级毛片| av视频在线观看入口| 成人国产一区最新在线观看| 国产又黄又爽又无遮挡在线| 亚洲欧美一区二区三区黑人| 亚洲 国产 在线| 最近最新中文字幕大全电影3| 亚洲欧美一区二区三区黑人| 亚洲 国产 在线| 最新美女视频免费是黄的| 一夜夜www| 久久精品人妻少妇| 久久午夜综合久久蜜桃| 狠狠狠狠99中文字幕| 日韩精品青青久久久久久| 国产高清三级在线| 免费人成视频x8x8入口观看| 精品人妻1区二区| 怎么达到女性高潮| 久久精品91无色码中文字幕| 少妇人妻一区二区三区视频| 99久久精品一区二区三区| 十八禁人妻一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲av片天天在线观看| 一本一本综合久久| 国产伦精品一区二区三区视频9 | 全区人妻精品视频| 成熟少妇高潮喷水视频| 嫩草影院精品99| 国产99白浆流出| 国产精品永久免费网站| 91av网站免费观看| 伊人久久大香线蕉亚洲五| 国产精品一区二区免费欧美| 夜夜爽天天搞|