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

    離心壓氣機管式擴壓器研究進展及評述

    2017-11-20 01:19:55韓戈陽誠武李紫良趙勝豐盧新根
    航空學報 2017年9期
    關鍵詞:壓器管式壓氣機

    韓戈, 陽誠武, 李紫良,2, 趙勝豐,2, 盧新根,2,*

    1.中國科學院工程熱物理研究所 輕型動力重點實驗室, 北京 100190 2.中國科學院大學, 北京 100190

    離心壓氣機管式擴壓器研究進展及評述

    韓戈1, 陽誠武1, 李紫良1,2, 趙勝豐1,2, 盧新根1,2,*

    1.中國科學院工程熱物理研究所 輕型動力重點實驗室, 北京 100190 2.中國科學院大學, 北京 100190

    緊湊高效的擴壓器設計非常具有挑戰(zhàn)性,成為制約高壓比離心壓氣機應用于工程實際的主要技術障礙。管式擴壓器是解決高壓比離心壓氣機擴壓器設計難題的有效手段。目前管式擴壓器已經(jīng)在國外先進中小型航空發(fā)動機中得到了應用,有效地提升了離心壓氣機性能。本文從管式擴壓器設計參數(shù)對壓氣機性能的影響以及管式擴壓器內(nèi)部復雜流動機理研究兩個方面對管式擴壓器的國內(nèi)外研究進展進行了回顧,討論了管式擴壓器內(nèi)部流動機理及其對離心壓氣機性能的影響,并指出了管式擴壓器研究的發(fā)展趨勢。

    離心壓氣機; 管式擴壓器; 壓氣機性能; 內(nèi)部流動; 非定常流動

    由于具有單級壓比高、零部件少、結(jié)構緊湊、可靠性高和抗外物撞擊能力強等一系列優(yōu)勢,離心壓氣機在中小型航空發(fā)動機中得到了廣泛的應用?,F(xiàn)代高性能航空發(fā)動機對高推(功)重比的需求使壓氣機的壓比不斷提高。然而高壓比使得離心葉輪出口氣流很不均勻,同時高壓比使得徑向擴壓器進口馬赫數(shù)增加,甚至出現(xiàn)超聲,再加上擴壓器入口與離心葉輪出口之間距離很小,兩者之間存在強烈的非定常相互作用,進一步惡化了葉輪與擴壓器之間的匹配;另一方面,飛行平臺對發(fā)動機外徑的限制使得擴壓器更加緊湊,導致擴壓器負荷進一步提高,離心壓氣機性能惡化,使得擴壓器成為了限制高壓比離心壓氣機性能提高的瓶頸,緊湊高效擴壓器的設計非常具有挑戰(zhàn)性。

    管式擴壓器特殊的燕尾型前緣形狀能夠較好地適應離心葉輪出口超聲速不均勻氣流,同時可以通過直接將氣流轉(zhuǎn)到軸向的轉(zhuǎn)彎方式來縮小擴壓器外徑,因此采用管式擴壓器既能夠提高高壓比離心壓氣機的性能,又能夠一定程度上減小壓氣機外徑,是解決高壓比離心壓氣機擴壓器設計難題的有效手段。目前管式擴壓器已經(jīng)在國外先進中小型航空發(fā)動機中得到了應用[1-2]。但是由于保密原因,公開發(fā)表的相關文獻很少。本文從管式擴壓器概述、管式擴壓器對壓氣機性能的影響以及管式擴壓器內(nèi)部流動等3個方面對管式擴壓器的研究進展進行介紹。

    1 管式擴壓器概述

    管式擴壓器由沿周向均布的一系列擴壓管道構成。一般情況下,管式擴壓器的單個擴壓管道由燕尾型前緣、喉部及緊接其后的擴壓通道構成。燕尾型前緣(見圖1及圖2)由管式擴壓器相鄰的2個通道在進口部分相貫形成,這種燕尾型前緣對于高馬赫數(shù)的不均勻氣流適應性較好。雖然管式擴壓器具有較復雜的前緣結(jié)構,然而從氣流通道來看,管式擴壓器與傳統(tǒng)的葉片擴壓器和楔形擴壓器相比卻具有相似的結(jié)構(見圖1[3-4]),它們都存在無葉區(qū)、半無葉區(qū)、擴壓器喉部和擴壓通道。管式擴壓器與葉片擴壓器和楔形擴壓器相比唯一區(qū)別在于管式擴壓器進口部分有一個獨特的區(qū)域,這一區(qū)域由相鄰管式擴壓器通道相貫形成的燕尾型前緣組成,位于擴壓器的無葉區(qū)與半無葉區(qū)之間,該區(qū)域常稱之為偽無葉區(qū)。在管式擴壓器前緣之后圓柱形通道一般會向下游延伸一段,這一段圓柱形通道的截面積一直保持不變,即管式擴壓器的幾何喉部。因此管式擴壓器的幾何喉部是三維的,具有一定的長度,而并非像其他類型的擴壓器那樣幾何喉部為一個截面。在管式擴壓器喉部之后,氣流通過一段錐形擴壓通道減速擴壓。

    在錐形擴壓通道之后,管式擴壓器除了減速擴壓以外還需進一步將氣流轉(zhuǎn)到軸向。傳統(tǒng)的管式擴壓器中此部分擴壓通道采用的是魚尾型通道,如圖2(a)[5-6]所示。這種擴壓通道的顯著特點是將徑向擴壓器和軸向擴壓器設計成一體,也被稱為“ram’s horn”管式擴壓器。這些魚尾型通道被逐根插入到錐形擴壓通道中并用螺栓固定以實現(xiàn)其與管式擴壓器基座的連接。同時也存在另一種采用管式擴壓器的轉(zhuǎn)彎形式,即在管式擴壓器徑向通道后連接90° 轉(zhuǎn)彎通道,而后再連接軸向擴壓器。這樣氣流在管式擴壓器擴壓通道之后的90° 轉(zhuǎn)彎中混合,待氣流轉(zhuǎn)到軸向后再經(jīng)軸向擴壓器進一步擴壓整流,最終排出壓氣機,如圖2(b)[7]所示。

    圖1 管式擴壓器Fig.1 Pipe diffuser

    Vrana[8]于20世紀60年代首先提出了管式擴壓器的概念,并引起了離心壓氣機設計和研究者的很大興趣。起初加拿大聯(lián)合飛機公司(1975年改稱加拿大普惠公司,Pratt and Whitney Cananda,以下簡稱P&WC公司)[9]以及美國國家航空航天局(National Aeronautics and Space Administration,NASA)[10]對管式擴壓器進行了性能試驗,發(fā)現(xiàn)管式擴壓器是解決高亞聲和超聲速擴壓器設計的有效手段。而后,P&WC公司將管式擴壓器應用于其發(fā)動機型號中(見圖3)[11],例如PT6渦軸發(fā)動機以及JT15D渦扇發(fā)動機。

    圖2 2種管式擴壓器轉(zhuǎn)彎型式Fig.2 Two turning configurations of pipe diffuser

    雖然如此,對于管式擴壓器的設計以及內(nèi)部流動機理并不清楚,因此P&WC公司隨后在20世紀70年代針對管式擴壓器設計參數(shù)對壓氣機性能的影響進行了試驗研究[12]。為了進一步探明管式擴壓器內(nèi)部流動機理,P&WC公司從20世紀80年代開始到2004年借助于Carlton大學的低速大尺度管式擴壓器風洞試驗臺對魚尾型管式擴壓器進行了風洞試驗,測得了管式擴壓器內(nèi)部流場[13-15]。同時通用電器(General Electric,GE)公司也進行了類似的風洞試驗研究[16]。2000年以后,阿爾斯通燃機(Alstom Gas Turbines)公司[17-18]以及本田公司[19]也都加入了管式擴壓器研究的行列。其中阿爾斯通燃機公司在英國Cranfield大學的支持下對管式擴壓器幾何參數(shù)對壓氣機性能的影響進行了研究。2008年以后,研究者們開始著手對管式擴壓器內(nèi)部流場進行測量,德國亞琛理工大學[7,20-24]以及加拿大的幾所大學[5-6,25-26]分別應用先進流場測量手段為GE公司和P&WC公司測得了管式擴壓器內(nèi)部的流場。最近管式擴壓器內(nèi)部的非定常流動以及離心葉輪和管式擴壓器前緣的非定常相互作用成為了研究的熱點,川崎重工[27]、P&WC公司[28-30]以及GE公司[31]均對管式擴壓器內(nèi)的非定常流動進行了研究。

    目前P&WC公司的PW600系列渦扇發(fā)動機、JT-15D小型渦扇發(fā)動機和PT6系列渦軸發(fā)動機的離心壓氣機均應用了管式擴壓器,并且在流動機理研究的基礎上不斷改進管式擴壓器的氣動及結(jié)構設計(見圖4)[1]。而包括GE公司在內(nèi)的其他航空發(fā)動機及燃機公司也都在積極開展管式擴壓器研究和設計工作。

    圖3 P&WC公司早期的管式擴壓器[11]Fig.3 Early pipe diffuser of P&WC[11]

    圖4 管式擴壓器結(jié)構[1]Fig.4 Structure of pipe diffuser[1]

    2 管式擴壓器對壓氣機性能的影響

    常見的擴壓器有:無葉擴壓器、低稠度擴壓器、葉片式擴壓器以及楔形擴壓器等。一般認為無葉擴壓器具有較為寬廣的裕度但損失較大;葉片式擴壓器和楔形擴壓器效率較高但裕度相對較?。坏统矶葦U壓器則介于二者之間。因此無葉擴壓器和低稠度擴壓器多用于對效率要求不苛刻但是需要有寬廣裕度的過程工程中,例如制冷壓縮機和渦輪增壓器。而葉片擴壓器及楔形擴壓器在高壓比離心壓氣機中應用較多[32-35]。相比于楔形擴壓器,管式擴壓器具有更高的性能,圖5為P&WC公司的PT6發(fā)動機離心壓氣機分別采用葉片式擴壓器和管式擴壓器時的級效率對比,其中橫坐標為壓氣機級壓比,縱坐標為效率變化。可以看出采用管式擴壓器后,壓氣機在不同工作點的性能都有了提高。

    早期對于管式擴壓器的研究集中在探究其對壓氣機性能的影響。P&WC公司的Kenny[9]于1969年最早證實了管式擴壓器在高性能高壓比離心壓氣機上應用的可行性。Kenny通過對比壓比5∶1和6∶1的2個離心壓氣機分別采用管式擴壓器和其他類型擴壓器的性能,發(fā)現(xiàn)采用管式擴壓器時離心壓氣機的壓比和效率比采用其他類型擴壓器時均有較大提高,詳見表1。通過分析還發(fā)現(xiàn)管式擴壓器的喉部堵塞較小且總壓損失較低。

    1973年,美國NASA的Lewis研究中心的Klassen[10]對一個壓比為1.9的低壓比離心葉輪與不同擴壓器相匹配時的性能進行了對比。其中原始擴壓器為葉片擴壓器,同時采用具有錐形和喇叭形2種不同擴張段的管式擴壓器。實驗結(jié)果發(fā)現(xiàn)原擴壓器具有最高的效率和靜壓恢復系數(shù),穩(wěn)定工作范圍也最寬廣。而通過增大2種管式擴壓器的喉部面積能夠提高擴壓器的性能。由此可見當管式擴壓器應用于壓比較低的離心壓氣機時,并沒有優(yōu)勢。

    管式擴壓器形狀復雜,影響其性能的參數(shù)較多,因此在對管式擴壓器有了初步研究之后,部分學者對管式擴壓器的設計參數(shù)對壓氣機性能的影響進行了研究。1977年,Reeves[12]發(fā)表了3組高壓比離心壓氣機帶不同截面形狀的管式徑向擴壓器(擴壓器未轉(zhuǎn)到軸向,氣流徑向排出)的實驗結(jié)果。實驗結(jié)果顯示方形進口截面的管式擴壓器具有較大的穩(wěn)定工作范圍,但是類橢圓形進口截面的管式擴壓器性能較高。Izumi和Hirose[36]認為設計時使管式擴壓器通道中心線的相切圓直徑小于葉輪外徑能夠減小管式擴壓器前緣與葉輪出口的距離,減小甚至消除葉輪出口氣流角與管式擴壓器前緣幾何角的不匹配,進一步提高管式擴壓器的靜壓恢復系數(shù)。英國Cranfield大學Bennett等[17-18]在國家工程與物理科學研究委員會(Engineering and Physical Sciences Research Council)和阿爾斯通燃氣輪機公司的支持下,對管式擴壓器喉部面積及擴壓器通道數(shù)目等關鍵參數(shù)對壓氣機性能的影響進行了實驗研究。發(fā)現(xiàn)當喉部面積的大小使得壓氣機設計流量與擴壓器堵塞流量之比為0.98時壓氣機性能最好,而截面形狀的變化對于壓氣機性能影響不大。但是這一研究將管式擴壓器出口設置在了喉部之后,所以在實驗過程中管式擴壓器沒有喉部之后的擴壓通道,存在一定的局限性。2014年,Han等[3]研究了管式徑向擴壓器進口與葉輪出口半徑比、喉部長度、喉部面積以及擴張角等幾何參數(shù)對壓氣機性能的影響,認為管式擴壓器進口與葉輪出口半徑比在1.03左右較好,喉部長度與喉部直徑之比為0~0.5較好,喉部面積存在最佳值,而擴張角不宜大于6°。2016年,Han等[37]以采用魚尾型管式擴壓器的高壓比單級離心壓氣機為研究對象,對魚尾型管式擴壓器的面積分布規(guī)律以及中心線形狀等幾何參數(shù)對壓氣機性能和擴壓器內(nèi)部流動的影響進行了數(shù)值研究,得到了最佳面積分布規(guī)律以及壓氣機性能隨二維中心線形狀變化時的圖譜。

    圖5 葉片式擴壓器與管式擴壓器效率對比[2]Fig.5 Comparison of stage efficiency of vaned diffuser and that of pipe diffuser[2]

    表1 不同擴壓器效率比較Table 1 Comparison of efficiency of different diffusers

    總壓比擴壓器類型絕熱效率/%5∶1平板擴壓器73葉片擴壓器75管式擴壓器81.86∶1葉片擴壓器73管式擴壓器81.3

    此外,F(xiàn)ilipenco等[38]在不同來流條件下對管式擴壓器進行了實驗研究,發(fā)現(xiàn)管式擴壓器的靜壓恢復系數(shù)主要受擴壓器進口氣流角影響。但擴壓器進口軸向的不均勻性以及擴壓器進口馬赫數(shù)基本上對擴壓器的靜壓恢復系數(shù)不產(chǎn)生影響。而Han等[4]對同一壓氣機分別采用管式擴壓器和楔形擴壓器時在不同工況點的性能進行了對比,發(fā)現(xiàn)采用管式擴壓器時在設計轉(zhuǎn)速壓氣機性能要高于采用楔形擴壓器的壓氣機性能,然而在低轉(zhuǎn)速時采用管式擴壓器的壓氣機性能反而不如采用楔形擴壓器的壓氣機性能。

    國內(nèi)對于管式擴壓器性能的研究相對較少。2011年,王毅等[39]以德國“Radiver”高壓比離心壓氣機為研究對象,應用數(shù)值方法對管式徑向擴壓器喉部面積變化對壓氣機特性的影響進行了研究。同樣發(fā)現(xiàn),當喉部面積使設計流量與擴壓器堵塞流量之比為0.98時,壓氣機獲得最佳性能。2012年,王博和嚴明[40]通過數(shù)值方法分析了進口形狀和壓力面局部變形對壓氣機性能的影響,發(fā)現(xiàn)類橢圓形進口管式擴壓器的壓比、效率和裕度均比楔形擴壓器高;圓形進口管式擴壓器具有比類橢圓形進口的管式擴壓器更高的壓比和效率,但是裕度大大減小;壓力面局部變形有利于提高管式擴壓器的性能。2014年,韓戈等[41]通過數(shù)值方法得到了管式擴壓器的喉部長度和擴張角2個參數(shù)對壓氣機性能的影響規(guī)律,認為2個參數(shù)的取值均不宜太大。

    從發(fā)表的文獻來看,國內(nèi)外學者研究了管式擴壓器關鍵參數(shù)對壓氣機性能的影響,得到了部分參數(shù)對壓氣機性能的影響規(guī)律。然而研究中部分學者對管式擴壓器做了簡化,因此無法反映管式擴壓器內(nèi)部的真實流動情況,具有一定的局限性。

    3 管式擴壓器內(nèi)部流動研究

    3.1 低速風洞研究

    圖6 Carleton大學的風洞試驗臺[42]Fig.6 Wind tunnel of Carleton University[42]

    圖7 不同進口條件下擴壓器出口流場[13] Fig.7 Outlet flow field of diffuser with different inlet conditions[13]

    針對魚尾型管式擴壓器展開的風洞試驗研究,測得了魚尾型管式擴壓器內(nèi)的低速流場,加深了對管式擴壓器內(nèi)部流動的理解。然而所采用的研究對象為簡化后的魚尾型管式擴壓器,既不包含燕尾型前緣結(jié)構,又簡化了魚尾型擴壓通道的幾何形狀,因此無法反映管式擴壓器內(nèi)部的真實流動情況,同樣具有一定的局限性。

    3.2 內(nèi)部流動研究

    隨著流動測量手段的發(fā)展,部分學者開始對真實管式擴壓器內(nèi)部流動進行測量。2008年以來,德國亞琛理工大學噴氣推進和葉輪機械研究所在GE公司航空發(fā)動機項目支持下,持續(xù)多年開展了管式擴壓器研究。Zachau和Buescher[20]針對一個采用管式擴壓器的離心壓氣機展開研究,應用粒子圖像測速法(Particle Image Velocimetry,PIV)測得了管式擴壓器內(nèi)部部分截面的流場(實驗臺如圖8(a)所示,流場如圖9(a)所示,圖中Cnorm為無量綱速度)。研究發(fā)現(xiàn):管式擴壓器前半段在壓力面發(fā)生流動分離,而擴壓器通道后30%中擴壓已經(jīng)不明顯,因此推測適當縮短出口段長度應該對性能不會有太大影響。2009年,Zachau和Niehuis[21]繼續(xù)對管式擴壓器進行了更細致的實驗研究,應用PIV測量技術,分析了離心葉輪葉頂間隙、離心葉輪出口抽氣以及葉輪和擴壓器軸向相對位置3個參數(shù)對壓氣機流場和性能的影響。分析發(fā)現(xiàn)增大離心葉輪葉頂間隙會使壓氣機性能下降;減小離心葉輪出口抽氣會使壓氣機特性向小流量方向移動,性能有所下降;而改變離心葉輪與擴壓器的軸向相對位置對于壓氣機性能基本沒有影響。Kunte等[22]在Zachau和Niehuis的研究基礎之上進一步對離心葉輪葉頂間隙和離心葉輪出口抽氣對壓氣機性能和擴壓器流場的影響進行了研究。發(fā)現(xiàn)增大離心葉輪葉頂間隙會使擴壓器內(nèi)射流向吸力面及葉頂移動,增大流動分離。減少離心葉輪葉頂抽氣會使擴壓器通道內(nèi)邊界層增厚,導致擴壓器壓力面的流動分離提前發(fā)生,降低管式擴壓器性能。Kunte等[7]在研究團隊之前成果基礎上,將管式擴壓器的后30%切掉,通過PIV以及油流測量擴壓器內(nèi)流場,發(fā)現(xiàn)切掉后30%后,管式擴壓器壓力面的流動分離得到了抑制,壓氣機效率不但沒有下降反而使絕熱效率提高了0.3%,但是擴壓器出口氣流角發(fā)生了改變。因此重新設計了串列軸向擴壓器代替了原來的軸向擴壓器,并進一步對串列軸向擴壓器兩排葉片的周向相對位置對壓氣機性能的影響進行了研究,發(fā)現(xiàn)當后排葉片靠近前排葉片的壓力面時,壓氣機性能最優(yōu)。通過采用串列軸向擴壓器以及切掉管式擴壓器的后30%使壓氣機性能提高了1.3%,驗證了Zachau和Buescher的猜想[20]。Wilkosz等[23]在Zachau及Kunte研究的基礎之上進一步為切掉后30%的管式擴壓器設計了1個單排葉片的軸向擴壓器,葉片前緣在轉(zhuǎn)彎通道之內(nèi)且長度與Kunte設計的串列軸向擴壓器相當。數(shù)值計算發(fā)現(xiàn)該軸向擴壓器較好地組織了內(nèi)部流動,性能略高于Kunte等[7]設計的串列擴壓器。2015年,Schmidt等[24]對原有的具有串列葉片的軸向擴壓器的前面1排葉片加長,直接與管式徑向擴壓器相接。通過實驗對比了具有單排葉片的軸向擴壓器[23]、具有串列葉片的軸向擴壓器[7]和具有加長前排葉片的串列軸向擴壓器的性能,發(fā)現(xiàn)兩種串列擴壓器的壓比和效率相差不大,靜壓恢復系數(shù)略高,但是加長第1排葉片后串列擴壓器的裕度比另外2種軸向擴壓器小2%。

    圖8 管式擴壓器實驗裝置Fig.8 Test rigs of pipe diffuser

    可以看出德國亞琛理工大學對管式擴壓器的研究很大程度上集中在葉輪參數(shù)變化對擴壓器的影響以及管式擴壓器之后的軸向擴壓器的研究,對于管式擴壓器本身涉及的內(nèi)容并不多。

    圖9 管式擴壓器內(nèi)部流場Fig.9 Internal flow field of pipe diffuser

    也有部分學者對P&WC公司采用的魚尾型管式擴壓器內(nèi)部流動進行了研究。2004年,Roberts和Steed[5]針對某型發(fā)動機軸流-離心組合壓氣機的后面級(串列離心葉輪+管式擴壓器),分別采用k-ε和SST湍流模型進行了設計轉(zhuǎn)速下定常數(shù)值模擬,簡單分析了葉輪和擴壓器內(nèi)流場,并將特性與實驗數(shù)據(jù)比對。研究發(fā)現(xiàn)SST湍流模型具有更好的吻合度。Bourgeois等[6,25]對一采用串列葉輪和魚尾型管式擴壓器的離心壓氣機進行了數(shù)值及實驗研究。采用k-ε、SST以及RSM-SSG 3種湍流模型對該壓氣機進行了數(shù)值模擬,并用激光多普勒測速法(Laser Doppler Velocimetry,LDV)測出壓氣機內(nèi)部幾個截面上的流場(實驗臺如圖8(b)所示,出口流場見圖9(b),圖中Cx為軸向速度,U2為離心葉輪輪緣速度)。通過分別對比計算以及實驗的總體性能和壓氣機進口、轉(zhuǎn)靜子之間和擴壓器出口3個截面上的流場,發(fā)現(xiàn)數(shù)值結(jié)果與實驗結(jié)果符合較好,相比于其它2種湍流模型,SST模型能夠更好地反映壓氣機的總體性能以及內(nèi)部流動狀況。2016年,Ali等[26]應用曲率修正后的SST湍流模型(SST-CC)對Bourgeois在2009年計算的采用串列葉輪和魚尾型管式擴壓器的離心壓氣機重新進行了計算并且與采用SST計算模型的計算結(jié)果進行了對比。研究發(fā)現(xiàn)SST-CC湍流模型對壓氣機的效率和內(nèi)部流場的預測與SST模型相差不大,但是對于壓氣機壓比的預測要好于SST湍流模型。

    雖然P&WC公司已經(jīng)將管式擴壓器應用到了航空發(fā)動機上,但是目前能夠看到的相關文獻并不多??梢钥吹竭@些對于魚尾型管式擴壓器的研究雖然測得了魚尾型管式擴壓器內(nèi)部流場,但是其研究重點放在了數(shù)值算法和湍流模型上,對于管式擴壓器的內(nèi)部流動機理關注不夠。

    國內(nèi)對于管式擴壓器內(nèi)部流動的研究很少。胡駿等[43]認為管式擴壓器效率高的主要原因有3點:① 管式擴壓器前緣的銳邊能夠起到減弱激波強度的作用,適合超聲速來流;② 管式擴壓器前緣銳邊能夠把不均勻來流進行調(diào)勻,從而也可以減少損失;③ 對于小流量高增壓比的離心壓氣機來講,擴壓器的軸向?qū)挾群苄。~片擴壓器在角區(qū)二次流損失嚴重,而在管式擴壓器的圓管通道中,二次流損失大為減少。

    3.3 非定常研究

    管式擴壓器特殊的燕尾型前緣與離心葉輪出口距離比其他類型的擴壓器都要近,存在著較強的非定常相互作用。因此對采用管式擴壓器的離心壓氣機進行非定常數(shù)值研究,以探究管式擴壓器特殊前緣結(jié)構與葉輪的非定常相互作用機理十分有必要。Gould等[28]對采用管式擴壓器時離心葉輪葉片上的負荷變化較為關注。通過應用二維及三維非定常數(shù)值方法得到了主葉片以及分流葉片上的非定常負荷。計算結(jié)果顯示,2種方法得到的流場和特性相差不大,在設計點得到的結(jié)果與實驗結(jié)果較為接近,而在其他點偏差較大。離心葉輪葉片非定常峰值負荷為定常負荷的3.5倍,大流量工況葉片尾緣的非定常負荷在向葉片前緣傳播時衰減較快。Ding[29]進一步對Roberts和Steed[5]以及Bourgeois等[6,25]所研究的壓氣機進行了非定常數(shù)值模擬研究,并與實驗數(shù)據(jù)進行了詳盡比對。發(fā)現(xiàn)通過采用非定常數(shù)值模擬可以進一步提高壓氣機性能的預測精度。Grates等[30]對一采用管式擴壓器的跨聲速壓氣機進行了非定常數(shù)值研究。數(shù)值計算得到的時間平均特性與實驗結(jié)果符合良好。分析發(fā)現(xiàn)擴壓器通道內(nèi)主要存在2個渦:一個是由于負攻角在擴壓器前緣產(chǎn)生的通道渦,靠近擴壓器葉頂,渦的強度隨著時間變化而變化;另一個渦是由于流體經(jīng)過擴壓器的燕尾型前緣發(fā)生回流而產(chǎn)生的,靠近擴壓器葉根。2個渦相互作用使得高能流體向壓力面移動,低能流體從壁面向通道中心遷移(見圖10),對于抑制氣流分離有重要作用。2014年,Sugimoto等[27]對帶有管式擴壓器的高壓比單級離心壓氣機進行了實驗研究,并且對其進行了定常和非定常數(shù)值計算,其中非定常數(shù)值結(jié)果通過Harmonic方法獲得。對比數(shù)值和實驗結(jié)果發(fā)現(xiàn)Harmonic方法非定常計算的準確性較高,時均特性與實驗結(jié)果吻合較好。同時認為管式擴壓器前緣的特殊結(jié)構使得其進口幾何角在展向不斷變化,正好與葉輪出口葉頂?shù)臍饬鹘窍噙m應,改善了擴壓器進口葉頂區(qū)域的流動(見圖11,圖中α2為葉輪出口氣流角,α3V為擴壓器進口幾何角,①、②、③分別對應于近失速點、設計點和堵塞點)。同年,Wilkosz等[31]對管式擴壓器進行了非定常數(shù)值及實驗研究,發(fā)現(xiàn)非定常損失占到壓氣機設計點損失的2.4%,而非定常效應增加的損失主要來源于葉頂間隙泄漏流而非擴壓器。

    圖10 管式擴壓器前緣對渦[30]Fig.10 Leading edge vortices of pipe diffuser[30]

    圖11 管式擴壓器進口氣流角與幾何角分布[27]Fig.11 Distributions of flow angle and metal angle at pipe diffuser inlet[27]

    壓氣機中的非定常流動十分復雜,管式擴壓器與離心葉輪之間的非定常相互作用更是近幾年研究的熱點,雖然以上對于采用管式擴壓器的壓氣機非定常研究取得了一些積極研究成果,但是對于管式擴壓器的非定常研究,特別是離心葉輪與管式擴壓器的非定常相互作用還需進一步研究。

    4 管式擴壓器研究的發(fā)展趨勢

    對于管式擴壓器的研究已持續(xù)了半個世紀,國內(nèi)外學者利用實驗測試以及數(shù)值模擬等多種手段對管式擴壓器內(nèi)部流動進行了深入的研究,取得了不少積極成果。國內(nèi)外研究者一方面開展部分管式擴壓器關鍵參數(shù)對壓氣機性能的影響,得到部分參數(shù)對壓氣機性能的影響規(guī)律,為管式擴壓器的設計提供了參考;另一方面,利用實驗和數(shù)值模擬手段,開展管式擴壓器內(nèi)部復雜流動機理研究,初步揭示了管式擴壓器提高壓氣機性能的流動機理。

    然而相比于管式擴壓器在中小型航空發(fā)動機壓氣機研制中所取得的巨大效益而言,關于管式擴壓器內(nèi)部流動機理的認識以及關于離心葉輪與管式擴壓器之間非定常耦合研究,顯然滯后得多。盡管就管式擴壓器幾何設計、效果評定和內(nèi)部流動機理認識等方面,積累了一定的實驗數(shù)據(jù)和經(jīng)驗準則,盡管在某些局部問題上取得了突破性進展,但就目前管式擴壓器研究現(xiàn)狀來看,仍有一些問題困擾著設計者們。如:

    1) 雖然國內(nèi)外學者均對管式擴壓器設計參數(shù)對壓氣機性能的影響進行了研究,但是研究過于分散,缺乏完整系統(tǒng)的研究,與傳統(tǒng)擴壓器相比,尚沒有建立不同幾何參數(shù)對管式擴壓器性能影響的圖譜,管式擴壓器設計中幾何參數(shù)的選取仍然存在很大的盲目性。

    2) 應用管式擴壓器時,不同公司選用了不同的轉(zhuǎn)彎型式,例如P&WC公司偏向于采用魚尾型管式擴壓器,而GE公司則偏向于采用管式徑向擴壓器+軸向擴壓器的轉(zhuǎn)彎布局形式,而不同轉(zhuǎn)彎型式的優(yōu)劣尚無對比。因此,如何選擇管式擴壓器的轉(zhuǎn)彎型式也是一個亟待解決的問題。

    3) 管式擴壓器應用于高壓比離心壓氣機有利于提高壓氣機的性能,但是其內(nèi)部流動機理還不明確,尚沒有指導管式擴壓器設計的理論體系。

    因此,對于管式擴壓器依然需要進一步研究,今后研究的主要內(nèi)容應該包括以下幾個方面:

    1) 開展具有完整前緣的管式擴壓器風洞實驗,應用氣動探針及流動顯示技術測量管式擴壓器內(nèi)部包括前緣對渦以及魚尾型通道內(nèi)對渦在內(nèi)的渦系的產(chǎn)生、發(fā)展和相互作用機理。

    2) 系統(tǒng)開展級環(huán)境下管式擴壓器前緣幾何參數(shù)對離心壓氣機性能影響的實驗研究,獲得擴壓器實驗數(shù)據(jù),總結(jié)獲得不同來流條件下擴壓器典型性能圖譜,減小管式擴壓器設計的盲目性。

    3) 兼顧管式擴壓器氣動性能和加工難度,探索魚尾型管式擴壓器和管式徑向擴壓器+軸向擴壓器等不同轉(zhuǎn)彎布局型式的優(yōu)劣,為管式擴壓器工程設計提供指導。

    4) 利用高精度實驗和數(shù)值模擬手段探討不同關鍵幾何參數(shù)沿流程分布規(guī)律對管式擴壓器性能以及內(nèi)部二次流動的影響。

    5) 應用先進的穩(wěn)態(tài)和動態(tài)測試手段更加詳細地測量管式擴壓器燕尾型前緣與離心葉輪出口氣流的定常以及非定常耦合作用,得到管式擴壓器內(nèi)部前緣對渦以及通道內(nèi)渦系的產(chǎn)生和發(fā)展過程,進一步結(jié)合定常及非定常模擬手段,對擴壓器通道內(nèi)部流動進行詳細分析,揭示管式擴壓器提高壓氣機性能的內(nèi)部流動機理。

    6) 通過試驗及數(shù)值手段研究設計及非設計工況下管式擴壓器內(nèi)部流動特點以及其對壓氣機性能的影響。

    7) 總結(jié)管式擴壓器性能以及內(nèi)部流動研究結(jié)果,提煉管式擴壓器流動損失模型,發(fā)展管式擴壓器相關設計方法及設計準則,為緊湊式高性能管式擴壓器設計奠定理論基礎。

    [1] MOUSTAPHA H. Small gas turbine technology: Evolution and challenges: AIAA-2003-2559[R]. Reston, VA: AIAA, 2003.

    [2] WRONG C B. Turbine engine design: AIAA-1981-0915[R]. Reston, VA: AIAA, 1981.

    [3] HAN G, LU X G, ZHAO S F, et al. Parametric studies of pipe diffuser on performance of a highly loaded centrifugal compressor[J]. Journal of Engineering for Gas Turbines and Power, 2014, 136(12): 122604-1-122604-10.

    [4] HAN G, LU X G, ZHANG Y F, et al. Study of a highly-Loaded centrifugal compressor with pipe diffuser at design and off-design operating conditions: GT2015-43426[R]. New York: ASME, 2015.

    [5] ROBERTS D A, STEED R. A comparison of steady-state centrifugal stage CFD analysis to experimental rig data: Proceedings of the ANSYS user’s conference 10-13[R]. Pittsburgh, PA: ANSYS, 2004.

    [6] BOURGEOIS J A, ROBERT J M, ROBERTS D A, et al. Experimental and numerical investigation of an aero-engine centrifugal compressor: GT2009-59808[R]. New York: ASME, 2009.

    [7] KUNTE R, JESCHKE P, SMYTHE C. Experimental investigation of a truncated pipe diffuser with a tandem deswirler in a centrifugal compressor stage[J]. Journal of Turbomachinery, 2013, 135(3): 031019-1-031019-10.

    [8] VRANA J C. Diffuser for centrifugal compressor: US Patent, 3333762[P]. 1967-08-01.

    [9] KENNY D P. A novel low-cost diffuser for high-performance centrifugal compressor[J]. Journal of Engineering for Gas Turbines and Power, 1969, 91(1): 37-47.

    [10] KLASSEN H A. Performance of a low-pressure-ratio centrifugal compressor with four diffuser designs: NASA TN-D7237[R]. Washington, D.C.: NASA, 1973.

    [11] GROH F G, WOOD G M, KULP R S, et al. Evaluation of a high hub/tip ratio centrifugal compressor[J]. Journal of Basic Engineering, 1970(9): 419-428.

    [12] REEVES G. Design and performance of selected pipe-type diffusers: 77-GT-104[R]. New York: ASME, 1977.

    [13] YARAS M I. Effects of inlet conditions on the flow in a fishtail curved diffuser with strong curvature[J]. Journal of Fluids Engineering, 1996, 118(2): 772-778.

    [14] YARAS M I. Flow measurements in a fishtail diffuser with strong curvature[J]. Journal of Fluids Engineering, 1999, 121(1): 410-417.

    [15] YARAS M I, ORSI P. Effects of periodic inflow unsteadiness on the time-averaged velocity field and pressure recovery of a diffusing bend with strong curvature[J]. Journal of Fluids Engineering, 2004, 126(3): 229-237.

    [16] BLAIR L W, RUSSO C J. Compact diffusers for centrifugal compressors[J]. Journal of Aircraft, 1981, 19(1): 46-51.

    [17] BENNETT I, TOURLIDAKIS A, ELDER R L. Detailed measurements within a selection of pipe diffusers for centrifugal compressors: 88-GT-92[R]. New York: ASME, 1988.

    [18] BENNETT I, TOURLIDAKIS A, ELDER R L. The design and analysis of pipe diffusers for centrifugal compressors[J]. Proceedings of the Institution of Mechanical Engineers Part A: Journal of Power and Energy, 2000, 214: 87-96.

    [19] SHUICHI H, MINEYASU O, HIDEHIKO N, et al. Axial diffuser for a centrifugal compressor: US Patent, 0022586 A1[P]. 2009-01-22.

    [20] ZACHAU U, BUESCHER C. Experimental investigation of a centrifugal compressor stage with focus on the flow in the pipe diffuser supported by particle image velocimetry (PIV) measurements: GT2008-51538[R]. New York: ASME, 2008.

    [21] ZACHAU U, NIEHUIS R. Experiential investigation of the flow in the pipe diffuser of a centrifugal compressor stage under selected parameter variations: GT2009-59320[R]. New York: ASME, 2009.

    [22] KUNTE R, SCHWARZ P, WILKOSZ B, et al. Experimental and numerical investigation of tip clearance and bleed effects in a centrifugal compressor stage with pipe diffuser[J]. Journal of Turbomachinery, 2013, 135(1): 011005-1-011005-12.

    [23] WILKOSZ B, SCHMIDT J, GUENTHER C, et al. Numerical and experimental comparison of a tandem and a single vane deswirler used in an aero engine centrifugal compressor[J]. Journal of Turbomachinery, 2014, 136(4): 041005-1-041005-10.

    [24] SCHMIDT J, SCHWARZ P, WILKOSZ B, et al. Detailed performance analysis of a centrifugal compressor stage with pipe diffuser and immerged tandem deswirler: GT2015-43484[R]. New York: ASME, 2015.

    [25] BOURGEOIS J A, ROBERT J M, SAVORY E, et al. Assessment of turbulence model predictions for an aero-engine centrifugal compressor[J]. Journal of Turbomachinery, 2011, 133(1): 011025-1-011025-15.

    [26] ALI S, ELLIOTT K J, SAVORT E, et al. Investigation of the performance of turbulence models with respect to high flow curvature in centrifugal compressors[J]. Journal of Fluids Engineering, 2016, 138(5): 051101-1-051101-10.

    [27] SUGIMOTO T, KAWANISHI T, KUMAMARU H, et al. Performance investigation into supersonic diffuser for a high pressure centrifugal compressor: GT2014-25104[R]. New York: ASME, 2014.

    [28] GOULD K A, TAN C S, MACRORIE M. Characterization of unsteady impeller-blade loading in a centrifugal compressor with a discrete-passage diffuser: GT2007-28002[R]. New York: ASME, 2007.

    [29] DING M Y. CFD analysis of off-design centrifugal compressor operation and performance[D]. Toronto: University of Toronto, 2005.

    [30] GRATES D R, JESCHKE P, NIEHUIS R. Numerical investigation of the unsteady flow inside a centrifugal compressor stage with pipe diffuser[J]. Journal of Turbomachinery, 2014, 136(3): 031012-1-031012-14.

    [31] WILKOSZ B, ZIMMERMANN M, SCHWARZ P, et al. Numerical investigation of the unsteady interaction within a close-coupled centrifugal compressor used in an aero engine[J]. Journal of Turbomachinery, 2014, 136(4): 041006-1-041006-10.

    [32] JAPIKSE D. Centrifugal compressor design and performance[M]. Wilder VT: Concepts ETI, Inc., 1996.

    [33] 西安交通大學透平壓縮機教研室. 離心壓縮機原理[M]. 北京: 機械工業(yè)出版社, 1980.

    TURBINE AND COMPRESSOR LAB OF XI’AN JIAOTONG UNIVERSITY. Fundamentals of centrifugal compressor[M]. Beijing: China Machine Press, 1980 (in Chinese).

    [34] 李燕生, 陸桂林. 向心透平與離心壓氣機[M]. 北京: 機械工業(yè)出版社, 1987.

    LI Y S, LU G L. Radial inflow turbine and centrifugal compressor[M]. Beijing: China Machine Press, 1987 (in Chinese).

    [35] VAN DEN BRAEMBUSSCHE R A. Centrifugal compressors analysis & design: Course note 192[R]. Karman Institute for Fluid Dynamics, 2012.

    [36] IZUMI M, HIROSE E. Diffuser arrangement for centrifugal compressors: US Patent, 6537028[P]. 2003-05-25.

    [37] HAN G, LU X G, ZHANG Y F, et al. Study of geometric parameter influence on fishtail pipe diffuser performance: GT2016-57404[R]. New York: ASME, 2016.

    [38] FILIPENCO V G, DENIZ S, JOHNSTON J M, et al. Effects of inlet flow field conditions on the performance of centrifugal compressor diffusers: Part 1-discrete-passage diffuser[J]. Journal of Turbomachinery, 2000, 122(1): 1-10.

    [39] 王毅, 趙勝豐, 盧新根, 等. 高負荷離心壓氣機管式擴壓器特點及機理分析[J]. 航空動力學報, 2011, 26(3): 649-655.

    WANG Y, ZHAO S F, LU X G, et al. Analysis of characteristic and mechanism of pipe diffuser for a highly loaded centrifugal compressor[J]. Journal of Aerospace Power, 2011, 26(3): 649-655 (in Chinese).

    [40] 王博, 嚴明. 離心壓氣機管式擴壓器設計與計算分析[J]. 航空動力學報, 2012,27(6): 1303-1311.

    WANG B, YAN M. Design and numerical simulation of pipe diffuser in centrifugal compressor[J]. Journal of Aerospace Power, 2012,27(6): 1303-1311 (in Chinese).

    [41] 韓戈, 盧新根, 趙勝豐, 等. 喉部長度和擴張角對離心壓氣機管式擴壓器影響研究[J]. 推進技術, 2014, 35(12): 1607-1614.

    HAN G, LU X G, ZHAO S F, et al. Effects of diffuser throat length and divergence angle on performance of centrifugal compressor with pipe diffuser[J]. Journal of Propulsion Technology, 2014, 35(12): 1607-1614 (in Chinese).

    [42] SAROCH M F. Contributions to the study of turbo machinery aerodynamics: Part Ⅰ: Design of a fish-tail diffuser test sections; PartⅡ: Computations of the effects of AVDR on transonic turbine cascades[D]. Carleton: Carleton University, 1996.

    [43] 胡駿, 吳鐵鷹, 曹人靖. 航空葉片機原理[M]. 北京: 國防工業(yè)出版社, 2006.

    HU J, WU T Y, CAO R J. Fundamentals of aviation turbomachinery[M]. Beijing: National Defense Industry Press, 2006 (in Chinese).

    (責任編輯: 鮑亞平,李世秋)

    *Corresponding author. E-mail: xingenlu@hotmail.com

    A review of studies on pipe diffuser of centrifugal compressor

    HAN Ge1, YANG Chengwu1, LI Ziliang1,2, ZHAO Shengfeng1,2, LU Xin’gen1,2,*

    1.KeyLaboratoryofLight-dutyGas-turbine,InstituteofEngineeringThermophysics,ChineseAcademyofSciences,Beijing100190,China2.UniversityofChineseAcademyofSciences,Beijing100190,China

    It is a great challenge to design a diffuser with both high performance and compact size. The compact and high performance design of diffusers has become the main technical limitation for the application of diffusers on high pressure ratio centrifugal compressors. Using pipe diffusers are recognized as an effective method to solve the tough diffuser design problem for high pressure ratio centrifugal compressors. Pipe diffusers have been used for advanced overseas small and medium jet engines. The performance of the centrifugal compressor has been improved evidently by using pipe diffusers. This review mainly introduces the research development of pipe diffusers home and abroad in two perspectives: first, the influences of design parameters of the pipe diffuser on the stage performance of the centrifugal compressor; second, the complex mechanism of the internal flow in pipe diffusers. The influence of the mechanism of internal flow in pipe diffuser on the performance of the centrifugal compressor is discussed, and some recommendations are made for the development of future research on pipe diffusers.

    centrifugal compressor; pipe diffuser; compressor performance; internal flow; unsteady flow

    2016-11-17; Revised: 2017-02-26; Accepted: 2017-06-15; Published online: 2017-06-26 09:07

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

    s: National Natural Science Foundation of China (51606187); Defense Industry Technology Development Program (B1320133025)

    V231.3

    A

    1000-6893(2017)09-520949-12

    2016-11-17; 退修日期: 2017-02-26; 錄用日期: 2017-06-15; 網(wǎng)絡出版時間: 2017-06-26 09:07

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

    國家自然科學基金 (51606187); 國防基礎科研項目 (B1320133025)

    *通訊作者.E-mail: xingenlu@hotmail.com

    韓戈, 陽誠武, 李紫良, 等. 離心壓氣機管式擴壓器研究進展及評述[J]. 航空學報, 2017, 38(9): 520949. HAN G, YANG C W, LI Z L, et al. A review of studies on pipe diffuser of centrifugal compressor[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(9): 520949.

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

    10.7527/S1000-6893.2017.620949

    猜你喜歡
    壓器管式壓氣機
    管式太陽能集熱器的顆粒換熱模擬
    不同類型擴壓器下離心式壓縮機級性能的數(shù)值模擬與試驗研究
    壓縮機技術(2022年2期)2022-05-19 05:02:48
    軸流壓氣機效率評定方法
    重型燃氣輪機壓氣機第一級轉(zhuǎn)子葉片斷裂分析
    壓氣機緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    軸向及徑向擴壓器對小型高速離心風機氣動性能的影響研究
    風機技術(2019年4期)2019-06-24 05:42:18
    管式空氣空預器泄漏分析及改進
    不同類型擴壓器對離心式壓縮機結(jié)構振動與輻射噪聲影響的試驗研究*
    風機技術(2016年6期)2016-04-11 06:22:50
    平推流管式連續(xù)反應器合成耐熱ABS樹脂的研究
    中國塑料(2015年7期)2015-10-14 01:02:39
    旋入式檢驗級測壓器鑒選方法研究
    計量學報(2014年4期)2014-06-07 10:02:36
    别揉我奶头 嗯啊视频| 丝袜喷水一区| 午夜精品一区二区三区免费看| 岛国毛片在线播放| 久久欧美精品欧美久久欧美| 舔av片在线| 国产亚洲av嫩草精品影院| 国产精品久久久久久久久免| 天堂网av新在线| 婷婷精品国产亚洲av| 男人舔奶头视频| 久久精品影院6| 亚洲va在线va天堂va国产| 九九爱精品视频在线观看| 国产免费男女视频| 国产极品天堂在线| 男女边吃奶边做爰视频| 国产视频首页在线观看| 亚洲成人久久爱视频| 干丝袜人妻中文字幕| 十八禁国产超污无遮挡网站| 99视频精品全部免费 在线| 欧美一区二区精品小视频在线| 国产伦理片在线播放av一区 | 精品久久久久久久久av| 精华霜和精华液先用哪个| 18禁黄网站禁片免费观看直播| 18禁黄网站禁片免费观看直播| 在线免费观看的www视频| 欧美一区二区国产精品久久精品| 国产黄片视频在线免费观看| 性色avwww在线观看| 亚洲欧美中文字幕日韩二区| 成人鲁丝片一二三区免费| 亚洲国产精品久久男人天堂| 国产精品电影一区二区三区| 国产伦在线观看视频一区| 国产熟女欧美一区二区| 男的添女的下面高潮视频| 菩萨蛮人人尽说江南好唐韦庄 | 99riav亚洲国产免费| 欧洲精品卡2卡3卡4卡5卡区| 99视频精品全部免费 在线| 日韩欧美一区二区三区在线观看| 国产精华一区二区三区| 国产av麻豆久久久久久久| 国产免费一级a男人的天堂| 亚州av有码| 18禁裸乳无遮挡免费网站照片| 国产视频内射| 国产精品伦人一区二区| 亚洲国产欧美在线一区| 99热全是精品| 能在线免费观看的黄片| 一本久久精品| 丝袜美腿在线中文| 99九九线精品视频在线观看视频| 夜夜爽天天搞| 能在线免费观看的黄片| 成年av动漫网址| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av涩爱 | 一级毛片电影观看 | 亚洲内射少妇av| 校园人妻丝袜中文字幕| 1024手机看黄色片| 男人狂女人下面高潮的视频| 青春草视频在线免费观看| 天堂av国产一区二区熟女人妻| 欧美精品一区二区大全| av在线天堂中文字幕| 人妻系列 视频| 国产精品人妻久久久久久| 中文字幕av成人在线电影| 国产一区二区激情短视频| 欧美xxxx性猛交bbbb| 成人午夜精彩视频在线观看| 大型黄色视频在线免费观看| 成人漫画全彩无遮挡| 国产一级毛片在线| 国产老妇女一区| 国产精品一及| 一进一出抽搐动态| 日韩在线高清观看一区二区三区| 一区福利在线观看| 午夜老司机福利剧场| 久久久a久久爽久久v久久| 国产成人精品久久久久久| 亚洲av成人av| 成年免费大片在线观看| h日本视频在线播放| 久久精品国产亚洲av香蕉五月| 插阴视频在线观看视频| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 国产精品三级大全| 91精品国产九色| 免费无遮挡裸体视频| 亚洲综合色惰| 可以在线观看的亚洲视频| 在线国产一区二区在线| 日韩中字成人| 日本一二三区视频观看| kizo精华| 99久久精品热视频| 国产成人aa在线观看| 久久亚洲国产成人精品v| 中文欧美无线码| 亚洲av第一区精品v没综合| 亚洲久久久久久中文字幕| 99国产精品一区二区蜜桃av| 午夜福利成人在线免费观看| 你懂的网址亚洲精品在线观看 | 国产片特级美女逼逼视频| 国产av麻豆久久久久久久| av在线观看视频网站免费| 综合色丁香网| 亚洲美女搞黄在线观看| 国产精品一区二区性色av| 综合色av麻豆| 亚洲美女搞黄在线观看| 夜夜看夜夜爽夜夜摸| 亚洲熟妇中文字幕五十中出| 啦啦啦啦在线视频资源| 少妇高潮的动态图| 亚洲成人精品中文字幕电影| 国模一区二区三区四区视频| 丝袜喷水一区| 久久人人精品亚洲av| 国产黄片视频在线免费观看| 成人美女网站在线观看视频| av天堂在线播放| 久久综合国产亚洲精品| 久久人妻av系列| 久久热精品热| 一个人免费在线观看电影| 亚洲美女搞黄在线观看| 亚洲av中文av极速乱| 狠狠狠狠99中文字幕| 国产av一区在线观看免费| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 国产老妇伦熟女老妇高清| 国产黄片视频在线免费观看| 精品一区二区三区人妻视频| 人妻系列 视频| 国产精品永久免费网站| 欧美一级a爱片免费观看看| 国产白丝娇喘喷水9色精品| 此物有八面人人有两片| 亚洲av免费在线观看| 久久午夜福利片| 日本黄色片子视频| 看片在线看免费视频| 日本-黄色视频高清免费观看| 成人国产麻豆网| 观看美女的网站| 免费在线观看成人毛片| 欧美最黄视频在线播放免费| 精品熟女少妇av免费看| 久久久久久久久大av| 18禁裸乳无遮挡免费网站照片| 国产av一区在线观看免费| 中文资源天堂在线| 我的女老师完整版在线观看| 秋霞在线观看毛片| 亚洲国产高清在线一区二区三| 丰满乱子伦码专区| 国内久久婷婷六月综合欲色啪| 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 真实男女啪啪啪动态图| 美女大奶头视频| 青青草视频在线视频观看| 麻豆国产av国片精品| 国产日韩欧美在线精品| 日本成人三级电影网站| 国产蜜桃级精品一区二区三区| 欧美激情在线99| 日本黄色片子视频| 婷婷精品国产亚洲av| 精品人妻偷拍中文字幕| 亚州av有码| 中国美女看黄片| 蜜桃久久精品国产亚洲av| 国产免费男女视频| av天堂在线播放| 色综合站精品国产| 国产人妻一区二区三区在| 国产精品麻豆人妻色哟哟久久 | 嫩草影院入口| 亚洲不卡免费看| av福利片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av一区综合| 久久九九热精品免费| 亚洲不卡免费看| 两性午夜刺激爽爽歪歪视频在线观看| 日产精品乱码卡一卡2卡三| 国国产精品蜜臀av免费| 简卡轻食公司| 精品国内亚洲2022精品成人| 久久精品91蜜桃| 爱豆传媒免费全集在线观看| 六月丁香七月| 蜜桃亚洲精品一区二区三区| 欧美+日韩+精品| 神马国产精品三级电影在线观看| 久久久久久大精品| 日韩一区二区三区影片| 在线播放无遮挡| 午夜福利在线在线| 一级毛片我不卡| 18+在线观看网站| 亚洲综合色惰| 亚洲av电影不卡..在线观看| 人妻夜夜爽99麻豆av| 一级毛片久久久久久久久女| 国产大屁股一区二区在线视频| av天堂在线播放| 国产成人影院久久av| 波多野结衣巨乳人妻| 久久精品夜夜夜夜夜久久蜜豆| 乱系列少妇在线播放| 国产毛片a区久久久久| 欧美精品国产亚洲| 亚洲美女搞黄在线观看| 深夜a级毛片| 高清毛片免费观看视频网站| 免费无遮挡裸体视频| 亚洲精品国产成人久久av| 99久久人妻综合| 亚洲精品日韩在线中文字幕 | av在线蜜桃| 午夜亚洲福利在线播放| 中文字幕av在线有码专区| 亚洲欧美成人综合另类久久久 | 免费av观看视频| 久久久久久久久大av| 国产亚洲欧美98| 国产一级毛片七仙女欲春2| 日本免费a在线| 国内精品久久久久精免费| 欧美色欧美亚洲另类二区| 欧美日韩综合久久久久久| 听说在线观看完整版免费高清| 级片在线观看| 日韩精品有码人妻一区| 国产69精品久久久久777片| 成人特级黄色片久久久久久久| 免费看美女性在线毛片视频| 国产伦精品一区二区三区四那| 黄片无遮挡物在线观看| 能在线免费观看的黄片| 天堂中文最新版在线下载 | 免费人成在线观看视频色| 欧美bdsm另类| 国产亚洲欧美98| 成人毛片60女人毛片免费| 一个人观看的视频www高清免费观看| 国产午夜精品一二区理论片| 少妇的逼好多水| 99热这里只有精品一区| 亚洲精品日韩在线中文字幕 | 免费看a级黄色片| 亚洲四区av| 欧美不卡视频在线免费观看| 少妇裸体淫交视频免费看高清| 啦啦啦韩国在线观看视频| 91av网一区二区| 国产精品美女特级片免费视频播放器| 欧美人与善性xxx| 国产高清激情床上av| 日韩欧美精品免费久久| 国产高清三级在线| 精品一区二区三区人妻视频| 青春草亚洲视频在线观看| 一级av片app| 网址你懂的国产日韩在线| 国产精品av视频在线免费观看| 日本欧美国产在线视频| 亚洲欧美成人精品一区二区| 亚洲欧美精品专区久久| 一级二级三级毛片免费看| 麻豆精品久久久久久蜜桃| 麻豆成人av视频| 国产不卡一卡二| 欧美成人a在线观看| 一本精品99久久精品77| 午夜福利在线观看免费完整高清在 | 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 伊人久久精品亚洲午夜| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 亚洲国产精品国产精品| 国产极品精品免费视频能看的| 超碰av人人做人人爽久久| 国产又黄又爽又无遮挡在线| 在线天堂最新版资源| 国国产精品蜜臀av免费| 精品久久久噜噜| 99久久精品热视频| 国产91av在线免费观看| 亚洲自拍偷在线| 最近的中文字幕免费完整| 日本黄色视频三级网站网址| 观看美女的网站| 久久精品国产亚洲av天美| 少妇熟女aⅴ在线视频| 欧美色视频一区免费| 欧美人与善性xxx| 天堂√8在线中文| 国产美女午夜福利| 国产毛片a区久久久久| 亚洲国产欧美在线一区| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 色综合站精品国产| 欧美极品一区二区三区四区| 夜夜看夜夜爽夜夜摸| 男女边吃奶边做爰视频| 国产精品一二三区在线看| 欧美不卡视频在线免费观看| 国产精品三级大全| 久久草成人影院| 欧美三级亚洲精品| 18禁在线播放成人免费| 国产av一区在线观看免费| 91午夜精品亚洲一区二区三区| 丝袜美腿在线中文| 亚洲丝袜综合中文字幕| 日韩中字成人| 老司机影院成人| 色5月婷婷丁香| 我的女老师完整版在线观看| 人妻系列 视频| 三级男女做爰猛烈吃奶摸视频| 九色成人免费人妻av| 女的被弄到高潮叫床怎么办| 身体一侧抽搐| 青青草视频在线视频观看| 男女那种视频在线观看| 麻豆一二三区av精品| 亚洲欧美日韩无卡精品| 日韩亚洲欧美综合| 在线播放无遮挡| 99久久中文字幕三级久久日本| a级毛色黄片| 变态另类成人亚洲欧美熟女| 欧美性猛交╳xxx乱大交人| 亚洲一区高清亚洲精品| av天堂中文字幕网| 国国产精品蜜臀av免费| 别揉我奶头 嗯啊视频| 日韩欧美 国产精品| 国产高清有码在线观看视频| 亚洲丝袜综合中文字幕| 欧美又色又爽又黄视频| 国产精品久久视频播放| 久久亚洲国产成人精品v| 丰满人妻一区二区三区视频av| 桃色一区二区三区在线观看| 精品不卡国产一区二区三区| a级毛片a级免费在线| 蜜桃久久精品国产亚洲av| 欧美精品一区二区大全| 日本av手机在线免费观看| 真实男女啪啪啪动态图| 久久99精品国语久久久| 亚洲精品影视一区二区三区av| 精品国内亚洲2022精品成人| 日本与韩国留学比较| 午夜福利在线观看吧| 日韩一区二区视频免费看| 国内揄拍国产精品人妻在线| 狂野欧美白嫩少妇大欣赏| 伊人久久精品亚洲午夜| 五月玫瑰六月丁香| 嘟嘟电影网在线观看| 亚洲一区二区三区色噜噜| 日本av手机在线免费观看| 男女那种视频在线观看| 亚洲成人精品中文字幕电影| 欧美成人一区二区免费高清观看| 国产成人午夜福利电影在线观看| 九九久久精品国产亚洲av麻豆| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 精品久久国产蜜桃| 内地一区二区视频在线| 12—13女人毛片做爰片一| 少妇猛男粗大的猛烈进出视频 | 亚洲国产精品国产精品| 天天躁夜夜躁狠狠久久av| 99国产精品一区二区蜜桃av| 亚洲在线观看片| 国产一区二区三区av在线 | 欧美日韩一区二区视频在线观看视频在线 | 熟妇人妻久久中文字幕3abv| 亚洲综合色惰| 简卡轻食公司| 伦理电影大哥的女人| 国产成人精品婷婷| 夜夜夜夜夜久久久久| 99视频精品全部免费 在线| 国产一区二区亚洲精品在线观看| a级毛色黄片| 欧美三级亚洲精品| 校园人妻丝袜中文字幕| 观看免费一级毛片| av天堂在线播放| 99久久精品国产国产毛片| 午夜免费男女啪啪视频观看| 日本与韩国留学比较| 国产精品精品国产色婷婷| 高清在线视频一区二区三区 | 美女 人体艺术 gogo| 欧美日韩一区二区视频在线观看视频在线 | 成人毛片60女人毛片免费| 精品免费久久久久久久清纯| 精品国内亚洲2022精品成人| 看非洲黑人一级黄片| 99在线视频只有这里精品首页| 成人高潮视频无遮挡免费网站| 少妇的逼好多水| 中文字幕人妻熟人妻熟丝袜美| 色噜噜av男人的天堂激情| 人妻久久中文字幕网| 黄片无遮挡物在线观看| 亚洲成人久久爱视频| 欧美日韩一区二区视频在线观看视频在线 | 日韩人妻高清精品专区| 国模一区二区三区四区视频| 日韩国内少妇激情av| 亚洲国产欧美人成| 十八禁国产超污无遮挡网站| eeuss影院久久| 人妻系列 视频| 国产爱豆传媒在线观看| 一本久久中文字幕| www.av在线官网国产| 免费大片18禁| 国产av麻豆久久久久久久| 午夜a级毛片| 婷婷色av中文字幕| 亚洲成人久久性| 亚洲va在线va天堂va国产| 99久久九九国产精品国产免费| a级毛片a级免费在线| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 国产单亲对白刺激| 久久精品人妻少妇| 性色avwww在线观看| 成年版毛片免费区| а√天堂www在线а√下载| 在线国产一区二区在线| 天堂√8在线中文| 乱码一卡2卡4卡精品| 草草在线视频免费看| 午夜精品国产一区二区电影 | 欧美精品一区二区大全| 美女cb高潮喷水在线观看| 黄色一级大片看看| or卡值多少钱| 国内精品久久久久精免费| 99热只有精品国产| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 亚洲欧美精品自产自拍| 97热精品久久久久久| 欧美激情久久久久久爽电影| 国产黄色小视频在线观看| 国产成人午夜福利电影在线观看| 床上黄色一级片| 联通29元200g的流量卡| 99热只有精品国产| 国产乱人偷精品视频| 一区福利在线观看| 国产v大片淫在线免费观看| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 99久久久亚洲精品蜜臀av| 春色校园在线视频观看| 一个人看的www免费观看视频| 一区二区三区四区激情视频 | 深爱激情五月婷婷| 国产午夜精品一二区理论片| 国产亚洲91精品色在线| 国产在线精品亚洲第一网站| 日韩制服骚丝袜av| 中文字幕免费在线视频6| 91av网一区二区| 成人二区视频| 免费看a级黄色片| 久久精品国产亚洲网站| 欧美+亚洲+日韩+国产| 乱系列少妇在线播放| 国内精品美女久久久久久| 美女cb高潮喷水在线观看| 精品不卡国产一区二区三区| 午夜福利在线观看吧| 婷婷亚洲欧美| 欧美在线一区亚洲| 少妇熟女aⅴ在线视频| 一夜夜www| 日韩,欧美,国产一区二区三区 | 欧美zozozo另类| 成人永久免费在线观看视频| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 最近视频中文字幕2019在线8| 久久久久久久久中文| 淫秽高清视频在线观看| 美女被艹到高潮喷水动态| 亚洲欧洲日产国产| 99热精品在线国产| 九色成人免费人妻av| 免费人成在线观看视频色| 有码 亚洲区| 99久久精品一区二区三区| 国产精品三级大全| 99热全是精品| 蜜臀久久99精品久久宅男| 美女 人体艺术 gogo| 又粗又爽又猛毛片免费看| 日韩 亚洲 欧美在线| 22中文网久久字幕| 国产黄色视频一区二区在线观看 | 黄片无遮挡物在线观看| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 国产精品蜜桃在线观看 | 99国产精品一区二区蜜桃av| 大香蕉久久网| 久久韩国三级中文字幕| 免费无遮挡裸体视频| 天天躁日日操中文字幕| 老熟妇乱子伦视频在线观看| 欧美成人a在线观看| 99视频精品全部免费 在线| 亚洲精品国产av成人精品| 天堂中文最新版在线下载 | 日本一二三区视频观看| av.在线天堂| 国产精华一区二区三区| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 全区人妻精品视频| 18禁在线无遮挡免费观看视频| 成人亚洲欧美一区二区av| 一区二区三区免费毛片| 人人妻人人澡人人爽人人夜夜 | 国产午夜精品论理片| 黄色视频,在线免费观看| 国产伦在线观看视频一区| 日韩成人av中文字幕在线观看| 国产精品乱码一区二三区的特点| 岛国在线免费视频观看| 婷婷亚洲欧美| 一区二区三区免费毛片| 插逼视频在线观看| av在线播放精品| .国产精品久久| 美女cb高潮喷水在线观看| 亚洲电影在线观看av| 乱系列少妇在线播放| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久精品电影小说 | 亚洲精品久久国产高清桃花| 亚洲中文字幕日韩| 国产精品一区www在线观看| av天堂在线播放| 久久草成人影院| 菩萨蛮人人尽说江南好唐韦庄 | 非洲黑人性xxxx精品又粗又长| 久久99蜜桃精品久久| 嘟嘟电影网在线观看| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久v下载方式| 男女做爰动态图高潮gif福利片| 又爽又黄a免费视频| 国产一区亚洲一区在线观看| 免费一级毛片在线播放高清视频| 国产精品精品国产色婷婷| 岛国毛片在线播放| 一本久久中文字幕| 91久久精品电影网| 又爽又黄无遮挡网站| 精品一区二区免费观看| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女| 欧美日韩国产亚洲二区| 亚洲国产精品成人综合色| 国产精品.久久久| 欧美成人一区二区免费高清观看| 国产精品久久久久久久电影| 亚洲欧美中文字幕日韩二区| 日本五十路高清| 人妻久久中文字幕网| 国产高清激情床上av| 十八禁国产超污无遮挡网站| 色综合站精品国产| 不卡视频在线观看欧美| 国产精品久久视频播放| 1000部很黄的大片| 久久久久性生活片| 午夜a级毛片| 午夜精品一区二区三区免费看| 国产一区二区三区在线臀色熟女| 欧美激情久久久久久爽电影| 男插女下体视频免费在线播放| 日本免费一区二区三区高清不卡| 欧美3d第一页|