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

    低雷諾數(shù)下橫流-射流中剪切渦的試驗研究

    2017-11-22 09:58:22張保雷上官燕琴王嫻陳剛李躍明
    航空學(xué)報 2017年7期
    關(guān)鍵詞:渦的背風(fēng)雷諾數(shù)

    張保雷,上官燕琴,王嫻,*,陳剛,李躍明

    1.西安交通大學(xué) 航天航空學(xué)院 機械結(jié)構(gòu)強度與振動國家重點實驗室,西安 710049

    2.西安交通大學(xué) 航天航空學(xué)院 陜西省先進飛行器服役環(huán)境與控制重點實驗室,西安 710049

    低雷諾數(shù)下橫流-射流中剪切渦的試驗研究

    張保雷1,2,上官燕琴1,2,王嫻1,2,*,陳剛1,2,李躍明1,2

    1.西安交通大學(xué) 航天航空學(xué)院 機械結(jié)構(gòu)強度與振動國家重點實驗室,西安 710049

    2.西安交通大學(xué) 航天航空學(xué)院 陜西省先進飛行器服役環(huán)境與控制重點實驗室,西安 710049

    為深入分析橫流-射流(JICF)的流動特性及其中的復(fù)雜渦系結(jié)構(gòu),從流動機理上研究燃機葉片氣膜冷卻,揭示高溫燃氣流與冷卻流的摻混機理,本文對橫向流中單孔射流所形成的剪切渦開展了試驗研究。主要研究了速度比、雷諾數(shù)及射流角對JICF所形成剪切渦的影響。結(jié)果表明:速度比、雷諾數(shù)以及射流角會改變主流與射流之間的摻混程度,從而改變射流軌跡的曲率、高度及垂向滲透能力,最終改變剪切渦的特性;迎風(fēng)渦與背風(fēng)渦分別是由射流邊界層渦與主流邊界層渦形成的,當主流邊界層渦強度大于射流邊界層渦時,背風(fēng)渦是流場的主導(dǎo)結(jié)構(gòu),反之,迎風(fēng)渦將成為流場的主要渦系結(jié)構(gòu)。

    橫流-射流;速度比;雷諾數(shù);射流角度;剪切渦

    當流體從縫隙或者孔中以一定角度射出,且與外界環(huán)境中不同流向的流體相互作用時,便產(chǎn)生了橫向流中的射流(Jet In Cross-Flow,JICF)問題,如圖1所示。JICF在工程中具有重要應(yīng)用價值,如航空發(fā)動機和地面燃機的透平葉片采用的氣膜冷卻,冷卻氣體從氣膜孔射出,與高溫主流相互作用,在氣膜孔附近及下游會產(chǎn)生非常復(fù)雜的三維非定常流動結(jié)構(gòu)[1]。氣膜冷卻效果難以準確預(yù)測的固有因素是葉片表面的復(fù)雜流場[2],其優(yōu)化設(shè)計一直是工業(yè)界的難題之一。對高溫燃氣流與冷卻流的摻混機理研究有賴于對JICF的深入研究。又如,鍋爐燃燒室中,控制燃燒必須要調(diào)節(jié)燃油和空氣的摻混,最終也歸結(jié)為JICF問題。而諸多國防尖端技術(shù),如沖壓發(fā)動機的性能分析,火箭推進器的設(shè)計以及垂直起降飛行器(V/STOL)懸停時的氣動特性分析,這些問題的優(yōu)化、解決都依賴于對射流和主流間相互作用及摻混規(guī)律的深入了解。JICF作為流體力學(xué)中一種基本、典型的流動,涉及空氣動力學(xué)、環(huán)境學(xué)、燃燒學(xué)、熱力學(xué)等。同時,流動過程中大尺度擬序結(jié)構(gòu)的形成、發(fā)展、卷入、破裂等動力學(xué)演化過程對擬序結(jié)構(gòu)理論、渦動力學(xué)和非線性動力學(xué)的研究具有重要理論意義。圖2為Fric和Roshko[3]提出的JICF典型渦結(jié)構(gòu)概念圖。

    目前關(guān)于JICF的研究十分廣泛。典型如:Gutmark等[4]利用二維PIV(Particle Image Velocimetry)技術(shù)試驗研究了吹風(fēng)比為3時,橫流-射流中單孔、串聯(lián)雙孔射流的動力學(xué)性質(zhì),研究的流動特性包括射流軌跡、滲透深度、射流前后緣的擴展及尺寸大小。Salewski等[5]利用大渦模擬(LES)及試驗方法研究了橫流-射流中的擬序結(jié)構(gòu)和摻混特性,對比了圓形、橢圓形、方形射流孔的平均速度場及湍流統(tǒng)計量。Kelso等[6]通過流動可視化技術(shù)與熱線風(fēng)速儀研究了速度比為2~6橫向流中圓孔射流的流動結(jié)構(gòu),包括馬蹄渦、剪切渦和尾跡渦,并發(fā)現(xiàn)了在不同速度比下,尾跡渦在不同雷諾數(shù)下其不穩(wěn)定性機理也不相同。Lim等[7]試驗研究了JICF的流動結(jié)構(gòu),分析了渦環(huán)的變形及折疊,指出渦環(huán)的變形沒有完全顯示出實際的流動。Mahesh[8]總結(jié)了有關(guān)JICF的相關(guān)研究,并明確指出在低吹風(fēng)比和低雷諾數(shù)下JICF的摻混規(guī)律與高吹風(fēng)比、高雷諾數(shù)下大不相同。Sau和Mahesh[9]采用直接數(shù)值模擬(DNS)的方法對JICF中渦環(huán)的動力學(xué)特性及輸運、摻混特性進行了研究。關(guān)暉和吳錘結(jié)[10]利用LES對JICF進行了數(shù)值研究,發(fā)現(xiàn)JICF中占據(jù)主導(dǎo)地位的是渦環(huán)的拉伸、扭曲。郭婷婷等[11]研究了射流角度對射流軌跡的彎曲程度及湍動能的影響,但其并未對渦系結(jié)構(gòu)進行研究。本課題組基于GPU并行技術(shù),利用混合熱格子法對平板單孔射流、氣膜冷卻冷熱流摻混機理進行了大規(guī)模直接數(shù)值模擬與大渦模擬研究[12-16],使用近2億計算網(wǎng)格,捕捉到了精細的渦系結(jié)構(gòu),分析了渦系間的相互作用及其對氣膜冷卻效率的影響。雖然JICF的相關(guān)研究十分廣泛[17-21],但對于JICF中典型渦系結(jié)構(gòu)(剪切渦、馬蹄渦、尾跡渦及反對稱渦),大多數(shù)研究者只關(guān)注主導(dǎo)橫向主流與射流摻混過程中的反對稱渦(CVP)[22-25],對剪切渦的研究較少。然而,剪切渦也在摻混過程中承擔(dān)了重要角色,其特性對射流的滲透能力有重要影響,同時,也會影響反對稱渦的生成。與此同時,目前對于低雷諾數(shù)下的JICF相關(guān)研究也較少。Gopalan等[26]做了速度比為0.5~2.5的研究,發(fā)現(xiàn)以速度比2為界限,產(chǎn)生了不同的流場特征,其雷諾數(shù)高達1.9×104。Cambonie等[27]做了速度比為0.5~3.0下的試驗研究,雷諾數(shù)范圍為220~660,但是該文獻主要研究了流場中的反向旋轉(zhuǎn)渦對,并沒有對剪切渦的形成、發(fā)展進行相應(yīng)的研究分析。

    鑒于低雷諾數(shù)下JICF蘊含著豐富的物理信息,且目前此方面已有較多數(shù)值模擬結(jié)果,相關(guān)試驗研究很少。因此,本文設(shè)計、搭建試驗平臺,進行了低雷諾數(shù)下的平板單孔JICF試驗,并重點分析了不同流動參數(shù)對剪切渦特性的影響。試驗清晰地捕捉到了渦系結(jié)構(gòu)的發(fā)展過程,結(jié)果表明速度比、雷諾數(shù)以及射流角會改變主流與射流之間的摻混程度,速度比對流動結(jié)構(gòu)影響很大,雷諾數(shù)與射流角對剪切渦的發(fā)展特性及射流軌跡的高度、滲透能力有影響,最終改變剪切渦的特性。這些數(shù)據(jù)對于數(shù)值模擬結(jié)果的驗證,以及揭示高溫燃氣流與冷卻流的摻混機理有重要支撐作用。

    1 試驗裝置

    本試驗在低速風(fēng)洞中進行,該風(fēng)洞分為穩(wěn)定段、收縮段、試驗段和出口段。在穩(wěn)定段前部布置有蜂窩罩和尼龍濾網(wǎng),以消除環(huán)境中可能的較大湍渦。主流經(jīng)過收縮段后進入風(fēng)洞試驗段,如圖3所示。整套試驗裝置如圖4所示。圖5為試驗段示意圖,其截面尺寸(Lz×Ly)為460 mm×450 mm,長度(Lx)為1 000 mm,射流孔孔徑D為7 mm。

    在風(fēng)洞的尾部安裝有離心風(fēng)機,用以抽取空氣,該離心風(fēng)機與頻率控制臺相連,從而可以方便地調(diào)節(jié)主流速度,主流的最大速度為30 m/s。射流則由空氣壓縮機提供,空氣經(jīng)由流量調(diào)節(jié)閥穩(wěn)定流量后進入煙霧發(fā)生器,摻有少量煙霧的氣體經(jīng)過流量計,最終由集氣腔上方的射流孔進入試驗段。為獲取更為清晰的流動圖像,本試驗采用連續(xù)性片光源打亮呈現(xiàn)為一定白色的射流。片光源最大發(fā)射功率為5 W,波長為532 nm,片光源全角為60°。本試驗拍攝使用的高速攝像機為PHANTOM CAMERA V1210,最高拍攝速度可達7.1×105frame/s。試驗中流體速度較低,對流動圖像進行記錄時,相鄰圖像間的時間間隔為0.28 ms。利用相應(yīng)的Pcc軟件來實現(xiàn)對相機的控制、觸發(fā)、參數(shù)調(diào)整和數(shù)據(jù)處理。

    試驗中,速度比的計算公式為

    式中:Uj與U∞分別為射流速度與主流在實驗腔中的速度。射流速度可由流量計的讀數(shù)計算得到,主流速度直接由熱線風(fēng)速儀測得?;谏淞骺卓讖紻、主流速度的雷諾數(shù)公式為

    式中:ρ為主流空氣密度;μ為空氣的動力黏度;v為空氣的運動黏度。以射流孔圓心為坐標原點,相應(yīng)的坐標系如圖5所示,所獲取的試驗結(jié)果為展向中性面,即:z=0時,x Oy截面內(nèi)的流動數(shù)據(jù)。本試驗工況分別為:① 在主流速度(U∞=0.43 m/s,Re∞=203)不變的前提下,通過調(diào)節(jié)射流流量改變射流速度,實現(xiàn)速度比R分別為0.40、0.70、0.88、1.05、1.23、1.39;② 在保證速度比不變的情況下,利用頻率控制臺將離心風(fēng)機轉(zhuǎn)速調(diào)整為2、3、4、5 Hz,對應(yīng)的主流速度分別為0.25、0.43、0.65、0.86 m/s,從而使得雷諾數(shù)Re∞分別為118、203、307、406;③ 射流角度θ分別為60°、90°,通過更換試驗件實現(xiàn)。

    2 結(jié)果分析

    2.1 典型工況下剪切渦的形成與發(fā)展機制

    圖6為Re∞=203,θ=90°,R=0.88時剪切渦的形成發(fā)展圖像。從圖6可以看出,射流從射流孔射出后,射流前緣部分直接面對主流沖擊,前緣處的渦結(jié)構(gòu)會受到主流邊界層渦的較大影響,而射流后緣受到主流的沖擊較小。在主流的壓迫剪切作用下,射流產(chǎn)生了沿主流方向的加速度及速度,射流軌跡發(fā)生了彎曲,由于射流與主流間強烈的剪切作用,產(chǎn)生了開爾文-赫姆霍茲不穩(wěn)定性[3,6,28],最終形成了剪切渦。與此同時,在射流孔后緣附近形成了一個低壓區(qū),在負壓強梯度的驅(qū)使下,該區(qū)域出現(xiàn)了逆流現(xiàn)象。Fric與Roshko[3]認為這個低壓區(qū)的出現(xiàn)有利于射流后緣附近平板邊界層的分離,最終在一定條件下促使尾跡渦的形成。射流前緣的剪切渦為迎風(fēng)渦(Leading-Edge Vortex,LEV),其旋轉(zhuǎn)方向為逆時針,根據(jù)展向渦量的定義可知,渦量為正;射流后緣的剪切渦為背風(fēng)渦(Trailing-Edge Vortex,TEV),旋轉(zhuǎn)方向為順時針,渦量為負。往下游輸運過程中,迎風(fēng)渦被拉伸、卷吸、分離,背風(fēng)渦則逐漸匯聚,負渦量占據(jù)主導(dǎo),迎風(fēng)渦和背風(fēng)渦在向下游的運動中逐漸融合、卷入,逆時針旋轉(zhuǎn)的迎風(fēng)渦逐漸變?nèi)?順時針旋轉(zhuǎn)的背風(fēng)渦逐漸成為流場的主導(dǎo)結(jié)構(gòu),正渦量逐漸減小。在射流向下游輸運的過程中,在射流孔下游約為5.5倍孔徑處,由于低壓區(qū)的出現(xiàn),平板邊界層發(fā)生分離,剪切渦被拉伸并發(fā)生較大變形。由于此工況的速度比(R=0.88)較小,低壓區(qū)內(nèi)的逆壓力梯度較小,不利于尾跡渦的生成,因此在圖6中無法觀察到清晰的尾跡渦。Kelso等[6]的試驗研究表明,當速度比大于4時,可捕捉到清晰穩(wěn)定的尾跡渦。流動結(jié)構(gòu)繼續(xù)向下游發(fā)展,由于主流與射流的劇烈摻混造成了大量的動能損失,剪切渦逐漸分解、破裂為大量小尺度的渦系結(jié)構(gòu)。

    射流主體前后緣壓強差的存在加劇了主流對射流的壓迫、剪切作用,最終使得射流軌跡發(fā)生彎曲并形成剪切渦。剪切渦產(chǎn)生于射流與主流的初始摻混階段,在射流的兩側(cè)面上其強度最大,剪切渦的出現(xiàn)削弱了射流主體的動能,進一步增大了射流的彎曲程度。簡言之,x Oy平面內(nèi),射流孔附近的流場主要特征為射流的彎曲、變形,周期性卷起的剪切渦,以及剪切渦的拉伸、變形、重組。

    2.2 速度比對剪切渦的影響

    圖7為Re∞=203,θ=90°時,不同速度比情況下的剪切渦系圖像。從圖7的渦系結(jié)構(gòu)可以看出,不同速度比下,射流軌跡高度不同,射流在主流中的垂向滲透能力不同。本文研究的6種速度比下(圖7(a)~圖7(f)),射流的高度分別約為3D、5.5D、6D、6.3D、7.2D、8D。隨著速度比的增加,射流動量增加,射流對主流的垂向滲透能力增強,最終導(dǎo)致射流軌跡高度增大,這與文獻[29-30]得到的結(jié)論一致。此外,射流軌跡發(fā)生彎曲的位置也逐漸遠離射流孔,在6種速度比下,射流主體發(fā)生彎曲的位置高度分別約為0.8D、1D、1.4D、1.8D、2.5D、3D。隨著速度比的增大,射流軌跡的曲率減小,且在較小速度比時,射流首先發(fā)生彎曲,繼而產(chǎn)生剪切渦(圖7(a)~圖7(c));而在較大速度比時,剪切渦在射流發(fā)生彎曲之前產(chǎn)生(圖7(d)~圖7(f))。與此同時,速度比對流動渦系結(jié)構(gòu)影響很大。由于主流的橫向剪切作用,射流前緣形成了逆時針旋轉(zhuǎn)的迎風(fēng)渦。低速度比時(圖7(a)、圖7(b)),射流被主流強烈擠壓,射流與主流相互作用的結(jié)果是迎風(fēng)渦的渦強減弱并消失,背風(fēng)渦成為主要的渦系結(jié)構(gòu)。隨著速度比的增大(圖7(c)~圖7(f)),當射流離開射流孔后,射流發(fā)生變形與彎曲,同時交替出現(xiàn)逆時針旋轉(zhuǎn)的迎風(fēng)渦與順時針旋轉(zhuǎn)的背風(fēng)渦。值得注意的是,當速度比較大時(圖7(e)、圖7(f)),隨著主流與射流的摻混加劇,逆時針旋轉(zhuǎn)的迎風(fēng)渦成為了流場的主導(dǎo)結(jié)構(gòu)。該試驗觀察到的現(xiàn)象與課題組之前的數(shù)值研究結(jié)果[14]一致。在小速度比時,由于射流的動能較小,順時針旋轉(zhuǎn)的主流剪切層渦會抑制逆時針旋轉(zhuǎn)的射流剪切層渦的形成及發(fā)展,并在剪切作用的驅(qū)使下形成了順時針旋轉(zhuǎn)的背風(fēng)渦,最終在流場中占據(jù)主導(dǎo)地位。速度比增大,射流的動能增大,逆時針旋轉(zhuǎn)的射流剪切層渦的強度增大,從而抵抗了主流剪切層渦的抑制作用,并形成了逆時針旋轉(zhuǎn)的迎風(fēng)渦,當速度比增大到一定程度時,迎風(fēng)渦會成為主要的流動結(jié)構(gòu)。在上述研究的6種速度比工況下,均是在射流前緣首先發(fā)生凹陷,繼而逐漸引起射流后緣的失穩(wěn)。除R=0.40(圖7(a))時,其本身射流高度太小外,失穩(wěn)位置隨速度比的增大也逐漸靠近射流孔,圖7(b)~圖7(f)中,失穩(wěn)凹陷處距射流孔高度分別約為2.5D、2D、1.5D、1D、0.5D。在速度比很小時(圖7(a)),迎風(fēng)渦與背風(fēng)渦之間的干擾小,剪切渦的運動形態(tài)很穩(wěn)定。隨著速度比的增大,渦結(jié)構(gòu)之間干擾增強,剪切渦的形態(tài)逐漸紊亂。從圖7展示的結(jié)果可以看出,當速度比大于1(圖7(d)~圖7(f))時,剪切渦的連續(xù)性出現(xiàn)較大減弱,與小速度比(圖7(a)~圖7(c))的工況相比,在速度比大于1的工況下,剪切渦的形態(tài)更加紊亂、不規(guī)則,而且,剪切渦在距離射流孔更近的流向位置就開始發(fā)生變形與耗散。這是因為,射流初始速度大,使得剪切渦還沒來得及向下游輸運,上游相鄰的一個甚至多個剪切渦已經(jīng)趕上,多個剪切渦發(fā)生融合、卷入,使得剪切渦尺寸更大,渦系的形態(tài)更加紊亂,并在向下游的運動中不斷脫落形成了垂直并延伸至平板表面的尾跡渦[3]。

    最后,從圖7還可以發(fā)現(xiàn),迎風(fēng)渦和背風(fēng)渦的產(chǎn)生雖有時間差,但總是相伴而行,在向下游發(fā)展過程中逐漸融合,且下游剪切渦的渦核部分與上游相鄰的背風(fēng)渦相連,在自身沿主流方向的速度分量和主流的推動下,向下游流動。射流后下方產(chǎn)生的逆流區(qū),流動復(fù)雜,靠近射流的逆流區(qū)邊界部分總體呈現(xiàn)出與主流反向的運動,同時加劇主流與射流之間的摻混。在逆流區(qū)之前,剪切渦的運動較為穩(wěn)定,強度也較大,而當流過此區(qū)域時,其形態(tài)快速改變,渦運動變的紊亂,這說明逆流區(qū)在加劇主流與射流之間摻混程度的同時也干擾了剪切渦的發(fā)展。

    2.3 雷諾數(shù)對剪切渦的影響

    圖8為R=0.62,θ=90°時,不同雷諾數(shù)下的剪切渦系圖。圖8的結(jié)果顯示,隨著雷諾數(shù)的增大,主流對射流的壓迫作用增強,射流軌跡高度逐漸降低,射流的滲透能力下降,射流軌跡發(fā)生彎曲的位置離射流孔越近,曲率也越大。與此同時,隨著雷諾數(shù)的增大,主流對射流的輸運作用增強,剪切渦結(jié)構(gòu)越快發(fā)生變形與耗散(圖8(d))。而且,隨著雷諾數(shù)的增大,剪切渦的生成位置也逐漸由射流前緣向射流后緣發(fā)展,即初始剪切渦由逆時針旋轉(zhuǎn)的迎風(fēng)渦逐漸變?yōu)轫槙r針旋轉(zhuǎn)的背風(fēng)渦,且渦強度逐漸增大。這是因為順時針旋轉(zhuǎn)的主流邊界層渦的強度隨著雷諾數(shù)的增大而增大,抑制了逆時針旋轉(zhuǎn)的射流邊界層渦的形成與發(fā)展,最終使得順時針旋轉(zhuǎn)的背風(fēng)渦成為流場的主導(dǎo)渦系結(jié)構(gòu)。在沿主流方向流動過程中,逆時針旋轉(zhuǎn)的迎風(fēng)渦與順時針旋轉(zhuǎn)的背風(fēng)渦交替出現(xiàn)并連接形成鏈狀結(jié)構(gòu)。繼續(xù)往下游流動,渦強更大的背風(fēng)渦將上游渦強較弱的迎風(fēng)渦吞噬并發(fā)生變形、拉伸,最后耗散。在本文研究的雷諾數(shù)范圍內(nèi)可以看出,隨雷諾數(shù)的增大,剪切渦的變形越明顯。

    2.4 射流角對剪切渦的影響

    圖9為Re∞=203,R=1.10時,不同射流角度下的剪切渦系圖像。從圖9中可以發(fā)現(xiàn),隨著射流角度的增大,射流軌跡發(fā)生彎曲的位置離射流孔中心的流向距離減小,彎曲位置的曲率增大。這是因為隨著角度的增大,射流速度沿主流方向的速度分量減小,主流對射流的壓迫、加速作用越明顯。與此同時,射流對主流的垂向滲透能力增強。當射流角為60°時(圖9(a)),射流一經(jīng)流出便在距離射流孔約為0.8D的位置處形成迎風(fēng)渦。而當射流角增大至90°時

    (圖9(b)),在距離平板的垂向高度約為1.5D的位置處才發(fā)生凹陷繼而形成剪切渦。這個現(xiàn)象說明了在大傾斜角的工況下,射流的垂向滲透能力更強。

    3 結(jié) 論

    為深入揭示氣膜冷卻中高溫燃氣流與冷卻流的摻混機理,本文設(shè)計、搭建試驗臺,研究了JICF中剪切渦的形成、發(fā)展機制,結(jié)果發(fā)現(xiàn):

    1)射流孔附近區(qū)域流場的主要特征是射流彎曲、變形及剪切渦周期性卷起、拉伸、變形、重組。

    2)速度比、雷諾數(shù)以及射流角會改變主流與射流之間的摻混程度及射流在主流中的垂向滲透能力,從而改變射流軌跡的曲率與高度,最終改變剪切渦的特性。

    3)迎風(fēng)渦與背風(fēng)渦分別是由逆時針旋轉(zhuǎn)的射流邊界層渦與順時針旋轉(zhuǎn)的主流邊界層渦形成的。當主流邊界層渦強度大于射流邊界層渦時,背風(fēng)渦是流場的主導(dǎo)結(jié)構(gòu),反之,迎風(fēng)渦將成為流場的主要渦系結(jié)構(gòu)。

    本文對JICF中的剪切渦進行了試驗研究,得到了JICF中剪切渦的定性流動信息,今后將進行定量等更為深入的研究。此外,展向yOz平面內(nèi)的反向旋轉(zhuǎn)渦對(CVP),是JICF中的另一主要渦系結(jié)構(gòu),其演化機理與發(fā)展趨勢對JICF的諸多應(yīng)用,如氣膜冷卻效率、沖擊射流冷卻效率等有重要影響。今后將對CVP的演化過程與機理進行深入的試驗與數(shù)值模擬研究,為揭示JICF的復(fù)雜流動機理提供更多基礎(chǔ)。

    [1] TYAGI M,ACHARYA S.Large eddy simulation of film cooling flow from an inclined cylindrical jet[J].Journal of Turbomachinery,2003,125(4):734-742.

    [2] BOGARD D G,THOLE K A.Gas turbine film cooling[J].Journal of Propulsion&Power,2006,22(2):249-270.

    [3] FRIC T F,ROSHKO A.Vortical structure in the wake of a transverse jet[J].Journal of Fluid Mechanics,1994,279:1-47.

    [4] GUTMARK E J,IBRAHIM I M,MURUGAPPAN S.Dynamics of single and twin circular jets in cross flow[J].Experiments in Fluids,2011,50(3):653-663.

    [5] SALEWSKI M,STANKOVIC D,FUCHS L.Mixing in circular and non-circular jets in crossflow[J].Flow,Turbulence and Combustion,2008,80(2):255-283.

    [6] KELSO R M,LIM T T,PERRY A E.An experimental study of round jets in cross-flow[J].Journal of Fluid Mechanics,1996,306:111-144.

    [7] LIM T T,NEW T H,LUO S C.On the development of large-scale structures of a jet normal to a cross flow[J].Physics of Fluids,2001,13(3):770-775.

    [8] MAHESH K.The interaction of jets with crossflow[J].Annual Review of Fluid Mechanics,2013,45:379-407.

    [9] SAU R,MAHESH K.Dynamics and mixing of vortex rings in crossflow[J].Journal of Fluid Mechanics,2008,604:389-409.

    [10] 關(guān)暉,吳錘結(jié).湍流橫向射流的大渦模擬及其渦結(jié)構(gòu)特性[J].中國科學(xué),2006,36(6):662-677.GUAN H,WU C J.Large eddy simulation of turbulent transverse jet and its vortex structure[J].Science in China,2006,36(6):662-677(in Chinese).

    [11] 郭婷婷,徐忠,李少華.2種角度橫向紊動射流的實驗分析[J].西安交通大學(xué)學(xué)報,2003,37(11):1207-1210.GUO T T,XU Z,LI S H.Experimental study on turbulent jets injected obliquely into a crossflow[J].Journal of Xi'an Jiaotong University,2003,37(11):1207-1210(in Chinese).

    [12] LEI J,WANG X,XIE G,et al.Turbulent flow field analysis of a jet in cross flow by DNS[J].Journal of Engineering Thermophysics,2015,24(3):259-269.

    [13] LEI J,WANG X,XIE G.High performance computation of a jet in crossflow by lattice Boltzmann based parallel direct numerical simulation[J].Mathematical Problems in Engineering,2015,2015:1-11.

    [14] SHANGGUAN Y,WANG X,LI Y.Large-scaled simulation on the coherent vortex evolution of a jet in a crossflow based on lattice Boltzmann method[J].Thermal Science,2015,19(3):977-988.

    [15] 上官燕琴,王嫻,李躍明.基于格子Boltzmann方法的平板射流大渦模擬[J].計算物理,2015,32(6):669-676.SHANGGUAN Y Q,WANG X,LI Y M.High-performance numerical simulation of jet in cross-flow based on lattice Boltzmann method[J].Chinese Journal of Computational Physics,2015,32(6):669-676(in Chinese).

    [16] SHANGGUAN Y Q,WANG X,LI Y M.Investigation on the mixing mechanism of single-jet film cooling with various blowing ratios based on hybrid thermal lattice Boltzmann method[J].International Journal of Heat&Mass Transfer,2016,97:880-890.

    [17] SHAHZAD K,FLECK B A,WILSON D J.Small scale modeling of vertical surface jets in cross-flow:Reynolds number and downwash effects[J].Journal of Fluids Engineering,2007,129(3):311-318.

    [18] FEARN R,WESTON R P.Vorticity associated with a jet in a cross flow[J].AIAA Journal,2012,12(12):1666-1671.

    [19] COUSSEMENT A,GICQUEL O,DEGREZ G.Large eddy simulation of a pulsed jet in cross-flow[J].Journal of Fluid Mechanics,2012,695:1-34.

    [20] RUIZ A M,LACAZE G,OEFELEIN J C.Flow topologies and turbulence scales in a jet-in-cross-flow[J].Physics of Fluids,2015,27(4):531-553.

    [21] 朱呈祥,尤延鋮.橫向氣流中非牛頓液體射流直接數(shù)值模擬[J].航空學(xué)報,2016,37(9):2659-2668.ZHU C X,YOU Y C.Direct numerical simulation of non-Newtonian liquid jet in crossflow[J].Acta Aeronautica et Astronautica Sinica,2016,37(9):2659-2668(in Chinese).

    [22] MUPPIDIS,MAHESH K.Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet[J].Physics of Fluids,2006,18(8):085103.

    [23] CORTELEZZI L,KARAGOZIAN A R.On the formation of the counter-rotating vortex pair in transverse jets[J].Journal of Fluid Mechanics,2001,446:347-373.

    [24] DONNADIEU C,ORTIZ S,CHOMAZ J M,et al.Three-dimensional instabilities and transient growth of a counter-rotating vortex pair[J].Physics of Fluids,2009,21(9):94-102.

    [25] SCHLEGEL F,WEE D,MARZOUK Y M,et al.Contributions of the wall boundary layer to the formation of the counter-rotating vortex pair in transverse jets[J].Journal of Fluid Mechanics,2011,676:461-490.

    [26] GOPALAN S,ABRAHAM B M,KATZ J.The structure of a jet in cross flow at low velocity ratios[J].Physics of Fluids,2004,16(6):2067-2087.

    [27] CAMBONIE T,GAUTIER N,AIDER J L.Experimental study of counter-rotating vortex pair trajectories induced by a round jet in cross-flow at low velocity ratios[J].Experiments in Fluids,2013,54(3):1-13.

    [28] YUAN L L,STREET R L,FERZIGER J H.Large-eddy simulations of a round jet in crossflow[J].Journal of Fluid Mechanics,1999,379:71-104.

    [29] NEW T H,LIM T T,LUO S C.Effects of jet velocity profiles on a round jet in cross-flow[J].Experiments in Fluids,2006,40(6):859-875.

    [30] YUAN L L,STREET R L.Trajectory and entrainment of a round jet in crossflow[J].Physics of Fluids,1998,10(9):2323-2335.

    Experimental investigation on shear vortex of jet in cross-flow at low Reynolds number

    ZHANG Baolei1,2,SHANGGUAN Yanqin1,2,WANG Xian1,2,*,CHEN Gang1,2,Ll Yueming1,2

    1.State Key Laboratory for Strength and Vibration of Mechanical Structures,School of Aerospace,Xi'an Jiaotong University,Xi'an 710049,China
    2.Shaanxi Key Laboratory of Environment and Control for Flight Vehicle,School of Aerospace,Xi'an Jiaotong University,Xi'an 710049,China

    The flow characteristic and vortical structures associated with the jet in cross-flow(JlCF)is studied to reveal the mixing mechanism between coolant jet and hot cross-flow in film cooling of turbine blades.Experimental investigations on shear vortex of single-jet in cross-flow are carried out in this paper.The effects of velocity ratio,Reynolds number and inclined angle on the characteristics of shear vortex are studied mainly.Results show that the velocity ratio,Reynolds number and inclined angle have a great influence on the mixing between jet and cross-flow,and therefore affect the jet trajectory and jet penetration into cross-flow which results in various characteristics of shear vortex.We also find that the leading-edge shear vortex and trailing-edge shear vortex are formed by boundary layer vortex of jet flow and cross-flow,respectively.The trailing-edge shear vortex becomes the main flow structure when the cross-flow boundary layer vortex is stronger than jet boundary layer vortex.On the contrary,the leading-edge vortex becomes dominant.

    jet in cross-flow;velocity ratio;Reynolds number;inclined angle;shear vortex

    2016-10-04;Revised:2016-10-27;Accepted:2016-11-16;Published online:2016-12-14 16:29

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

    s:National Basic Research Program of China(2013CB035702);National Natural Science Foundation of China(11302165)

    O358

    A

    1000-6893(2017)07-120831-09

    10.7527/S1000-6893.2016.0305

    2016-10-04;退修日期:2016-10-27;錄用日期:2016-11-16;網(wǎng)絡(luò)出版時間:2016-12-14 16:29

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

    國家“973”計劃 (2013CB035702);國家自然科學(xué)基金 (11302165)

    *通訊作者.E-mail:wangxian@mail.xjtu.edu.cn

    張保雷,上官燕琴,王嫻,等.低雷諾數(shù)下橫流-射流中剪切渦的試驗研究[J].航空學(xué)報,2017,38(7):120831.ZHANG B L,SHANGGUANY Q,WANG X,et al.Experimental investigation on shear vortex ofjet in cross-flow at low Reynolds number[J].Acta Aeronautica et Astronautica Sinica,2017,38(7):120831.

    (責(zé)任編輯:王嬌)

    *Corresponding author.E-mail:wangxian@mail.xjtu.edu.cn

    猜你喜歡
    渦的背風(fēng)雷諾數(shù)
    完整
    長江文藝(2023年5期)2023-05-19 02:12:21
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運動研究
    海洋通報(2020年5期)2021-01-14 09:26:52
    溝槽對湍流邊界層中展向渦影響的實驗研究
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    開縫圓柱縫隙傾斜角對脫落渦的影響
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    三角柱體對機翼根部馬蹄渦的影響
    民機高速風(fēng)洞試驗的阻力雷諾數(shù)效應(yīng)修正
    The coupling characteristics of supersonic dual inlets for missile①
    熟女少妇亚洲综合色aaa.| 亚洲av美国av| 在线看a的网站| 国产成人一区二区三区免费视频网站| 脱女人内裤的视频| 97精品久久久久久久久久精品| 亚洲欧美色中文字幕在线| 成年人免费黄色播放视频| 超色免费av| 国产精品欧美亚洲77777| 亚洲专区字幕在线| 亚洲国产欧美日韩在线播放| 国产亚洲精品第一综合不卡| 性色av一级| 国产片内射在线| 国产精品一二三区在线看| 亚洲 欧美一区二区三区| 手机成人av网站| 精品国内亚洲2022精品成人 | 色综合欧美亚洲国产小说| 50天的宝宝边吃奶边哭怎么回事| 建设人人有责人人尽责人人享有的| 日本91视频免费播放| 在线精品无人区一区二区三| 亚洲黑人精品在线| 国产亚洲欧美精品永久| 99久久综合免费| 丰满饥渴人妻一区二区三| 久久精品熟女亚洲av麻豆精品| 99精国产麻豆久久婷婷| 99热网站在线观看| 精品久久久精品久久久| 丝袜喷水一区| 51午夜福利影视在线观看| 国产欧美日韩综合在线一区二区| 国产1区2区3区精品| 热re99久久精品国产66热6| 久久久精品免费免费高清| 黄频高清免费视频| 91麻豆av在线| 亚洲av男天堂| av电影中文网址| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 久久久精品94久久精品| 侵犯人妻中文字幕一二三四区| 美女午夜性视频免费| 国内毛片毛片毛片毛片毛片| 夜夜夜夜夜久久久久| 国产成人a∨麻豆精品| 免费在线观看影片大全网站| 亚洲男人天堂网一区| 满18在线观看网站| 亚洲国产中文字幕在线视频| 久久国产精品影院| 少妇被粗大的猛进出69影院| 亚洲精品日韩在线中文字幕| 成人三级做爰电影| 搡老熟女国产l中国老女人| 久久久久久久大尺度免费视频| 久久亚洲国产成人精品v| 这个男人来自地球电影免费观看| 我的亚洲天堂| 三上悠亚av全集在线观看| 极品少妇高潮喷水抽搐| 一级a爱视频在线免费观看| 成人国语在线视频| 久久人妻福利社区极品人妻图片| 亚洲精品一区蜜桃| 久热爱精品视频在线9| 啦啦啦在线免费观看视频4| 国产男女内射视频| 爱豆传媒免费全集在线观看| 免费在线观看黄色视频的| 欧美成狂野欧美在线观看| 国产黄色免费在线视频| 国产精品久久久久久人妻精品电影 | 精品一区二区三区四区五区乱码| 人人妻人人添人人爽欧美一区卜| 999久久久国产精品视频| www.999成人在线观看| 国产精品偷伦视频观看了| 国产高清videossex| 黄色视频在线播放观看不卡| 国产91精品成人一区二区三区 | 久久久精品94久久精品| 免费在线观看影片大全网站| 成在线人永久免费视频| 国产男人的电影天堂91| 亚洲三区欧美一区| 国产在线免费精品| 十分钟在线观看高清视频www| 一本久久精品| 午夜免费观看性视频| 母亲3免费完整高清在线观看| 精品欧美一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡 | 老司机午夜福利在线观看视频 | 国产精品久久久久久精品古装| 国产精品二区激情视频| 99久久国产精品久久久| 色老头精品视频在线观看| 精品视频人人做人人爽| 两性夫妻黄色片| 亚洲 欧美一区二区三区| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 久久这里只有精品19| 久久精品国产综合久久久| 在线观看人妻少妇| 最近最新免费中文字幕在线| 精品少妇久久久久久888优播| 悠悠久久av| 久久人妻福利社区极品人妻图片| 高清黄色对白视频在线免费看| 国产激情久久老熟女| 手机成人av网站| 他把我摸到了高潮在线观看 | 男女高潮啪啪啪动态图| 美女视频免费永久观看网站| 91av网站免费观看| 高清在线国产一区| 亚洲国产成人一精品久久久| 天堂中文最新版在线下载| 亚洲精品国产一区二区精华液| 欧美日韩精品网址| 亚洲成av片中文字幕在线观看| 久久香蕉激情| 不卡av一区二区三区| 国产精品香港三级国产av潘金莲| 国产色视频综合| 一级毛片精品| 日韩视频在线欧美| 成年人黄色毛片网站| 日日爽夜夜爽网站| 国产精品麻豆人妻色哟哟久久| 大型av网站在线播放| 久久ye,这里只有精品| 亚洲精品中文字幕在线视频| 国产伦人伦偷精品视频| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 久久精品人人爽人人爽视色| 亚洲国产中文字幕在线视频| 国产欧美日韩一区二区三区在线| 纵有疾风起免费观看全集完整版| 久久青草综合色| 韩国精品一区二区三区| 中亚洲国语对白在线视频| 人人澡人人妻人| 亚洲国产欧美日韩在线播放| 国产免费福利视频在线观看| 黄片播放在线免费| 日本一区二区免费在线视频| 精品高清国产在线一区| 日韩电影二区| 精品熟女少妇八av免费久了| 国产精品一区二区精品视频观看| 国产亚洲欧美在线一区二区| av一本久久久久| 亚洲国产欧美日韩在线播放| 亚洲黑人精品在线| 丝袜美足系列| 国产欧美日韩一区二区三 | 岛国在线观看网站| 在线天堂中文资源库| 国产97色在线日韩免费| 超碰97精品在线观看| 岛国毛片在线播放| 欧美少妇被猛烈插入视频| 999精品在线视频| 国产精品自产拍在线观看55亚洲 | 两人在一起打扑克的视频| 国产精品一区二区精品视频观看| 欧美日韩中文字幕国产精品一区二区三区 | 嫁个100分男人电影在线观看| 亚洲色图综合在线观看| av欧美777| 男人操女人黄网站| 美女国产高潮福利片在线看| 国产91精品成人一区二区三区 | 女人久久www免费人成看片| 伊人亚洲综合成人网| 成人国语在线视频| 91精品国产国语对白视频| 国产亚洲精品第一综合不卡| 欧美精品一区二区免费开放| 成年人午夜在线观看视频| 深夜精品福利| 中文精品一卡2卡3卡4更新| 一区二区三区激情视频| 亚洲精品中文字幕一二三四区 | 999久久久精品免费观看国产| 国产精品九九99| av网站在线播放免费| 婷婷成人精品国产| 啪啪无遮挡十八禁网站| 丰满少妇做爰视频| 精品人妻熟女毛片av久久网站| 十分钟在线观看高清视频www| 麻豆乱淫一区二区| 成年美女黄网站色视频大全免费| 少妇猛男粗大的猛烈进出视频| 国产淫语在线视频| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 国产三级黄色录像| 亚洲av电影在线观看一区二区三区| 女人高潮潮喷娇喘18禁视频| 成年人黄色毛片网站| 我要看黄色一级片免费的| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美在线一区| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 亚洲中文日韩欧美视频| 久久av网站| 最新在线观看一区二区三区| 午夜视频精品福利| 又大又爽又粗| 高清欧美精品videossex| 91精品三级在线观看| 国产区一区二久久| 久久久久久久国产电影| a在线观看视频网站| 99九九在线精品视频| 国产一卡二卡三卡精品| 美女扒开内裤让男人捅视频| 日本五十路高清| 又大又爽又粗| 美女国产高潮福利片在线看| 久久人妻福利社区极品人妻图片| 真人做人爱边吃奶动态| 在线精品无人区一区二区三| 国产人伦9x9x在线观看| 亚洲国产看品久久| 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 国产伦人伦偷精品视频| 精品国产国语对白av| 黄片小视频在线播放| 在线av久久热| 悠悠久久av| 最近最新免费中文字幕在线| 色老头精品视频在线观看| 亚洲中文av在线| 一级毛片电影观看| 日本欧美视频一区| 操美女的视频在线观看| 中文字幕制服av| 免费看十八禁软件| 岛国在线观看网站| 成在线人永久免费视频| 一区二区三区四区激情视频| 宅男免费午夜| 精品一区在线观看国产| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 91字幕亚洲| 免费人妻精品一区二区三区视频| 国产精品一区二区精品视频观看| 高清黄色对白视频在线免费看| 99国产精品99久久久久| 啦啦啦 在线观看视频| 久久久精品94久久精品| av欧美777| 国产区一区二久久| 十八禁高潮呻吟视频| 18在线观看网站| 欧美性长视频在线观看| 美女脱内裤让男人舔精品视频| 十分钟在线观看高清视频www| 97人妻天天添夜夜摸| 成人免费观看视频高清| 自拍欧美九色日韩亚洲蝌蚪91| 97人妻天天添夜夜摸| 亚洲九九香蕉| 亚洲av成人一区二区三| 国产亚洲欧美在线一区二区| 午夜免费鲁丝| 亚洲av美国av| 中文字幕最新亚洲高清| 人妻 亚洲 视频| 亚洲免费av在线视频| 成人av一区二区三区在线看 | 成年av动漫网址| 亚洲精品久久久久久婷婷小说| 国产精品久久久久成人av| 亚洲精品粉嫩美女一区| 久久久久久久久免费视频了| 久久久欧美国产精品| 成年女人毛片免费观看观看9 | 搡老熟女国产l中国老女人| 成年人免费黄色播放视频| av一本久久久久| 久久狼人影院| 久久久久国产一级毛片高清牌| av线在线观看网站| 精品国产一区二区久久| 999久久久国产精品视频| 久久性视频一级片| 秋霞在线观看毛片| 免费看十八禁软件| 一级黄色大片毛片| 日韩人妻精品一区2区三区| av不卡在线播放| 亚洲性夜色夜夜综合| 国产99久久九九免费精品| 黄片小视频在线播放| 欧美中文综合在线视频| av欧美777| 男人添女人高潮全过程视频| 亚洲成国产人片在线观看| 男人舔女人的私密视频| 精品国产乱码久久久久久小说| 亚洲精品在线美女| 久久精品亚洲av国产电影网| 妹子高潮喷水视频| 黄片大片在线免费观看| 激情视频va一区二区三区| 色精品久久人妻99蜜桃| 亚洲欧美日韩另类电影网站| 午夜福利在线观看吧| 欧美性长视频在线观看| 国产成人精品久久二区二区91| 黑人操中国人逼视频| 国产精品九九99| 色视频在线一区二区三区| kizo精华| 国产日韩一区二区三区精品不卡| 人妻一区二区av| svipshipincom国产片| 欧美在线黄色| 大片电影免费在线观看免费| av欧美777| 老司机福利观看| 免费黄频网站在线观看国产| 久久久久视频综合| 日韩 欧美 亚洲 中文字幕| 中文字幕色久视频| 国产成人系列免费观看| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 老司机午夜福利在线观看视频 | 亚洲,欧美精品.| √禁漫天堂资源中文www| 国产真人三级小视频在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲少妇的诱惑av| 桃花免费在线播放| 国产精品麻豆人妻色哟哟久久| 最新在线观看一区二区三区| 国产在线观看jvid| av又黄又爽大尺度在线免费看| av在线app专区| 性少妇av在线| 在线观看人妻少妇| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 爱豆传媒免费全集在线观看| 日本91视频免费播放| 99久久国产精品久久久| 91字幕亚洲| 波多野结衣一区麻豆| 日韩中文字幕欧美一区二区| 性色av乱码一区二区三区2| av天堂久久9| 激情视频va一区二区三区| 国产精品九九99| 狂野欧美激情性xxxx| 亚洲av片天天在线观看| 日韩欧美国产一区二区入口| 亚洲国产精品一区二区三区在线| 亚洲国产精品一区三区| 亚洲国产成人一精品久久久| 黄频高清免费视频| 9色porny在线观看| 丁香六月天网| 又大又爽又粗| 99久久99久久久精品蜜桃| 欧美黑人精品巨大| 欧美精品一区二区免费开放| 色播在线永久视频| 人妻 亚洲 视频| 18禁黄网站禁片午夜丰满| 老司机影院成人| 人人妻人人爽人人添夜夜欢视频| 亚洲av男天堂| 亚洲欧美激情在线| 色94色欧美一区二区| 亚洲第一欧美日韩一区二区三区 | 女人精品久久久久毛片| 亚洲成人免费电影在线观看| 日韩熟女老妇一区二区性免费视频| www.熟女人妻精品国产| 激情视频va一区二区三区| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 日韩欧美国产一区二区入口| 精品福利永久在线观看| 久久国产亚洲av麻豆专区| 少妇被粗大的猛进出69影院| 国产熟女午夜一区二区三区| 国产亚洲精品一区二区www | 欧美精品一区二区免费开放| 亚洲欧美清纯卡通| 成年动漫av网址| 欧美日韩一级在线毛片| 又黄又粗又硬又大视频| 热99久久久久精品小说推荐| 日韩制服骚丝袜av| 下体分泌物呈黄色| 在线精品无人区一区二区三| 最新在线观看一区二区三区| 亚洲精品第二区| 久久女婷五月综合色啪小说| 欧美一级毛片孕妇| 久久久久久免费高清国产稀缺| 满18在线观看网站| 午夜激情av网站| 天堂中文最新版在线下载| 我要看黄色一级片免费的| 狠狠婷婷综合久久久久久88av| 黑人操中国人逼视频| 国产男女超爽视频在线观看| 9191精品国产免费久久| 少妇裸体淫交视频免费看高清 | 精品人妻一区二区三区麻豆| 高清av免费在线| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 免费日韩欧美在线观看| 亚洲精品日韩在线中文字幕| 中文字幕制服av| 两人在一起打扑克的视频| 美女午夜性视频免费| 国产人伦9x9x在线观看| 国产三级黄色录像| 久久精品亚洲熟妇少妇任你| 国产成人影院久久av| 亚洲九九香蕉| 在线十欧美十亚洲十日本专区| 亚洲av欧美aⅴ国产| 一区二区三区激情视频| 国产成人一区二区三区免费视频网站| 在线观看人妻少妇| 精品一区在线观看国产| 国产精品九九99| 亚洲精品国产av成人精品| 亚洲欧美精品自产自拍| 又紧又爽又黄一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产野战对白在线观看| 国产成人啪精品午夜网站| 免费在线观看完整版高清| 免费在线观看日本一区| 精品国产一区二区三区久久久樱花| 国产精品熟女久久久久浪| 亚洲av国产av综合av卡| 久久人人97超碰香蕉20202| 一本—道久久a久久精品蜜桃钙片| 久久亚洲精品不卡| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 色婷婷久久久亚洲欧美| 国产精品自产拍在线观看55亚洲 | 日韩视频一区二区在线观看| 91精品国产国语对白视频| 日韩,欧美,国产一区二区三区| 丰满迷人的少妇在线观看| 亚洲成人国产一区在线观看| 国产亚洲欧美精品永久| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| 国产在线一区二区三区精| 国产亚洲欧美精品永久| 十八禁网站免费在线| 欧美乱码精品一区二区三区| 1024视频免费在线观看| 天天躁夜夜躁狠狠躁躁| 一级片'在线观看视频| 亚洲av美国av| av电影中文网址| 1024视频免费在线观看| 伊人久久大香线蕉亚洲五| 国产男人的电影天堂91| 国产精品二区激情视频| 亚洲精品国产av成人精品| 亚洲精品日韩在线中文字幕| 黑人巨大精品欧美一区二区蜜桃| 熟女少妇亚洲综合色aaa.| 蜜桃在线观看..| 中文字幕高清在线视频| 国产精品 国内视频| 日韩一区二区三区影片| 久久久精品国产亚洲av高清涩受| 午夜免费观看性视频| 人人妻人人添人人爽欧美一区卜| 肉色欧美久久久久久久蜜桃| 一区二区三区激情视频| 国产区一区二久久| 91麻豆精品激情在线观看国产 | 久久久久久久久久久久大奶| 中文字幕制服av| 亚洲综合色网址| 十八禁网站免费在线| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女 | 岛国在线观看网站| 天天添夜夜摸| 真人做人爱边吃奶动态| 国产在线免费精品| 精品国内亚洲2022精品成人 | 丰满人妻熟妇乱又伦精品不卡| 欧美97在线视频| 99热网站在线观看| 日本五十路高清| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 国产色视频综合| 黑人巨大精品欧美一区二区mp4| 国产成人精品久久二区二区91| 一边摸一边抽搐一进一出视频| 一区在线观看完整版| 亚洲成人免费av在线播放| 国产日韩一区二区三区精品不卡| 叶爱在线成人免费视频播放| 国产视频一区二区在线看| 日本vs欧美在线观看视频| 午夜福利一区二区在线看| 在线观看免费高清a一片| av网站免费在线观看视频| 纯流量卡能插随身wifi吗| 熟女少妇亚洲综合色aaa.| 亚洲精品粉嫩美女一区| 下体分泌物呈黄色| 黄色片一级片一级黄色片| 日韩电影二区| 91精品伊人久久大香线蕉| 中亚洲国语对白在线视频| www.自偷自拍.com| 欧美精品av麻豆av| 欧美日韩亚洲综合一区二区三区_| 建设人人有责人人尽责人人享有的| 免费不卡黄色视频| 久久天堂一区二区三区四区| 一区二区三区乱码不卡18| 国产精品偷伦视频观看了| 波多野结衣av一区二区av| 青青草视频在线视频观看| 欧美xxⅹ黑人| 每晚都被弄得嗷嗷叫到高潮| 777米奇影视久久| 亚洲精品美女久久久久99蜜臀| 岛国在线观看网站| bbb黄色大片| 午夜免费观看性视频| 午夜日韩欧美国产| 亚洲精品自拍成人| 91成人精品电影| 欧美久久黑人一区二区| 中文字幕人妻丝袜制服| 999精品在线视频| 亚洲精品国产一区二区精华液| 日韩熟女老妇一区二区性免费视频| 久久久久国产精品人妻一区二区| 一区二区av电影网| 久久天堂一区二区三区四区| 又黄又粗又硬又大视频| 午夜精品久久久久久毛片777| 日韩电影二区| 亚洲av片天天在线观看| 国产精品成人在线| 男人添女人高潮全过程视频| 亚洲中文日韩欧美视频| 一区二区三区精品91| 无遮挡黄片免费观看| 精品国内亚洲2022精品成人 | 国产一卡二卡三卡精品| 国产精品久久久久成人av| av福利片在线| 国产精品久久久人人做人人爽| 9191精品国产免费久久| 欧美成狂野欧美在线观看| 精品视频人人做人人爽| 亚洲av成人一区二区三| 欧美成狂野欧美在线观看| 日韩电影二区| 亚洲欧美一区二区三区久久| 欧美精品一区二区免费开放| 精品视频人人做人人爽| 9191精品国产免费久久| 亚洲熟女精品中文字幕| 久热这里只有精品99| 美国免费a级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 80岁老熟妇乱子伦牲交| 无限看片的www在线观看| 精品国产乱码久久久久久小说| 欧美午夜高清在线| 老司机亚洲免费影院| 国产极品粉嫩免费观看在线| 亚洲成人国产一区在线观看| 久久久久网色| 操美女的视频在线观看| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合一区二区三区| 久久久久久亚洲精品国产蜜桃av| 国产欧美亚洲国产| 少妇 在线观看| 一本色道久久久久久精品综合| 欧美日本中文国产一区发布| 亚洲五月色婷婷综合|