• <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在线直播| 能在线免费看毛片的网站| 久久久久网色| 自拍欧美九色日韩亚洲蝌蚪91 | 午夜激情久久久久久久| 久久av网站| 久久久久久久久久久丰满| 免费在线观看成人毛片| 毛片一级片免费看久久久久| 一级a做视频免费观看| 女人精品久久久久毛片| 欧美日韩精品成人综合77777| 男的添女的下面高潮视频| 美女脱内裤让男人舔精品视频| 国产高清不卡午夜福利| 国产一区有黄有色的免费视频| 免费在线观看成人毛片| 国产精品一区www在线观看| av国产精品久久久久影院| 我的女老师完整版在线观看| 国产一区二区三区av在线| 人体艺术视频欧美日本| 国产在线男女| 欧美日韩亚洲高清精品| 熟女人妻精品中文字幕| 美女cb高潮喷水在线观看| 久久久国产一区二区| 亚洲熟女精品中文字幕| 久久热精品热| 五月开心婷婷网| 亚洲丝袜综合中文字幕| 下体分泌物呈黄色| 国语对白做爰xxxⅹ性视频网站| 日日撸夜夜添| 午夜福利网站1000一区二区三区| 久久亚洲国产成人精品v| 亚洲国产成人一精品久久久| 国产美女午夜福利| 亚洲国产最新在线播放| av国产精品久久久久影院| 噜噜噜噜噜久久久久久91| 一级a做视频免费观看| 国产高清国产精品国产三级| 亚洲精品色激情综合| 99热国产这里只有精品6| 亚洲真实伦在线观看| 亚洲国产精品一区二区三区在线| av有码第一页| 人人妻人人看人人澡| 国产日韩欧美亚洲二区| 欧美3d第一页| 免费观看的影片在线观看| 中国三级夫妇交换| 伦理电影大哥的女人| 日韩av不卡免费在线播放| 中文在线观看免费www的网站| 在线亚洲精品国产二区图片欧美 | 国产亚洲午夜精品一区二区久久| 又大又黄又爽视频免费| 寂寞人妻少妇视频99o| 亚洲无线观看免费| 涩涩av久久男人的天堂| 亚洲不卡免费看| 在线观看www视频免费| 亚洲婷婷狠狠爱综合网| 久久午夜福利片| 日韩中文字幕视频在线看片| 久久99精品国语久久久| 天天操日日干夜夜撸| 永久网站在线| 日本黄色日本黄色录像| 亚洲无线观看免费| 亚洲国产精品999| 麻豆成人av视频| 女人久久www免费人成看片| 亚洲中文av在线| 熟女电影av网| 亚洲欧美精品自产自拍| 亚洲久久久国产精品| 久热这里只有精品99| 97在线人人人人妻| 在线亚洲精品国产二区图片欧美 | 人人澡人人妻人| a级毛片免费高清观看在线播放| 97超视频在线观看视频| 久久精品国产自在天天线| 亚洲精品久久久久久婷婷小说| 2022亚洲国产成人精品| 亚洲欧美精品专区久久| 自拍偷自拍亚洲精品老妇| 另类精品久久| av一本久久久久| 激情五月婷婷亚洲| 激情五月婷婷亚洲| 多毛熟女@视频| 国产精品免费大片| 亚洲国产欧美在线一区| 色婷婷av一区二区三区视频| 亚洲综合色惰| 国产精品偷伦视频观看了| 一区在线观看完整版| 只有这里有精品99| 欧美日韩精品成人综合77777| 久久精品久久久久久噜噜老黄| 哪个播放器可以免费观看大片| 美女国产视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 美女国产视频在线观看| 一级毛片电影观看| 久久久精品94久久精品| 校园人妻丝袜中文字幕| 免费观看性生交大片5| 卡戴珊不雅视频在线播放| 一级毛片我不卡| 亚洲四区av| 国产淫片久久久久久久久| 日本欧美国产在线视频| 中文乱码字字幕精品一区二区三区| 一本一本综合久久| 成年人午夜在线观看视频| 亚洲四区av| 男人爽女人下面视频在线观看| 一本一本综合久久| 国产精品久久久久久精品电影小说| 欧美精品亚洲一区二区| 麻豆乱淫一区二区| 男人狂女人下面高潮的视频| 晚上一个人看的免费电影| 亚洲不卡免费看| 中文在线观看免费www的网站| 高清在线视频一区二区三区| 高清在线视频一区二区三区| 国产免费又黄又爽又色| 精品一区二区三卡| 久久久欧美国产精品| 国产av一区二区精品久久| 蜜桃在线观看..| 亚洲在久久综合| 大香蕉久久网| 九草在线视频观看| 成人免费观看视频高清| 女性生殖器流出的白浆| 99精国产麻豆久久婷婷| 不卡视频在线观看欧美| 女性生殖器流出的白浆| 色5月婷婷丁香| 97超碰精品成人国产| 免费看不卡的av| 国产伦精品一区二区三区四那| 男女无遮挡免费网站观看| 新久久久久国产一级毛片| 黑人巨大精品欧美一区二区蜜桃 | 看免费成人av毛片| 伊人亚洲综合成人网| 成人亚洲欧美一区二区av| 日本黄大片高清| 亚洲一区二区三区欧美精品| av免费观看日本| 国内揄拍国产精品人妻在线| 成人国产麻豆网| 国产成人91sexporn| 在线观看三级黄色| 国产精品一区二区在线不卡| 亚洲色图综合在线观看| 男人和女人高潮做爰伦理| 欧美日韩av久久| 啦啦啦啦在线视频资源| tube8黄色片| 久久久欧美国产精品| 国产欧美亚洲国产| 久久久久久久久久人人人人人人| 欧美日韩综合久久久久久| 22中文网久久字幕| 久久精品国产自在天天线| 国产欧美另类精品又又久久亚洲欧美| 男女边吃奶边做爰视频| 国产高清有码在线观看视频| av不卡在线播放| 欧美精品亚洲一区二区| 免费黄频网站在线观看国产| 国产又色又爽无遮挡免| 黄色一级大片看看| av.在线天堂| 不卡视频在线观看欧美| 少妇被粗大的猛进出69影院 | 亚洲不卡免费看| 观看免费一级毛片| 男女无遮挡免费网站观看| 亚洲精品一二三| 欧美激情极品国产一区二区三区 | 99九九在线精品视频 | 99精国产麻豆久久婷婷| 国产精品国产三级专区第一集| 国产国拍精品亚洲av在线观看| 校园人妻丝袜中文字幕| 春色校园在线视频观看| 天堂俺去俺来也www色官网| 美女视频免费永久观看网站| 免费少妇av软件| 看免费成人av毛片| 国产免费一级a男人的天堂| 国内少妇人妻偷人精品xxx网站| 久久久久精品久久久久真实原创| 观看美女的网站| 欧美日韩亚洲高清精品| 国产在线免费精品| 99精国产麻豆久久婷婷| 欧美亚洲 丝袜 人妻 在线| 亚洲av.av天堂| 国产精品女同一区二区软件| 老司机亚洲免费影院| 三上悠亚av全集在线观看 | 国产免费一级a男人的天堂| 国产精品国产三级专区第一集| 十八禁网站网址无遮挡 | 自拍偷自拍亚洲精品老妇| 国产欧美日韩综合在线一区二区 | 日日摸夜夜添夜夜爱| 熟女人妻精品中文字幕| 久久久国产欧美日韩av| 亚洲性久久影院| 制服丝袜香蕉在线| 一本色道久久久久久精品综合| 成人无遮挡网站| 日韩一区二区三区影片| 国产色爽女视频免费观看| www.av在线官网国产| 亚洲精品久久久久久婷婷小说| 国产欧美另类精品又又久久亚洲欧美| 国产av国产精品国产| 国产精品久久久久久精品电影小说| 国产精品久久久久久久久免| 免费少妇av软件| 三级经典国产精品| 亚洲av国产av综合av卡| 国产一区亚洲一区在线观看| 校园人妻丝袜中文字幕| av女优亚洲男人天堂| 亚洲av综合色区一区| 精品视频人人做人人爽| 国产高清不卡午夜福利| 亚洲美女黄色视频免费看| 一边亲一边摸免费视频| 亚洲国产最新在线播放| 婷婷色av中文字幕| 80岁老熟妇乱子伦牲交| 免费观看av网站的网址| 亚洲熟女精品中文字幕| 亚洲国产精品一区三区| 99热全是精品| 亚洲三级黄色毛片| 欧美区成人在线视频| 日本与韩国留学比较| 亚洲欧洲日产国产| 在线 av 中文字幕| 黄色欧美视频在线观看| 成人免费观看视频高清| 国产午夜精品久久久久久一区二区三区| 日韩大片免费观看网站| 一区二区三区四区激情视频| 国产视频内射| 国产色爽女视频免费观看| av天堂中文字幕网| 久久精品国产自在天天线| 久久久久人妻精品一区果冻| 国产淫片久久久久久久久| 在线亚洲精品国产二区图片欧美 | 纯流量卡能插随身wifi吗| 老女人水多毛片| 不卡视频在线观看欧美| 欧美亚洲 丝袜 人妻 在线| 极品人妻少妇av视频| a级毛片免费高清观看在线播放| 国产精品熟女久久久久浪| 又黄又爽又刺激的免费视频.| 伦理电影大哥的女人| 精华霜和精华液先用哪个| 久久久久人妻精品一区果冻| 国产淫片久久久久久久久| 日本91视频免费播放| 欧美区成人在线视频| 国产成人免费观看mmmm| 91在线精品国自产拍蜜月| 免费黄色在线免费观看| 成人综合一区亚洲| 国产男女超爽视频在线观看| 亚洲欧美精品专区久久| 亚洲av福利一区| 国产成人a∨麻豆精品| a级毛色黄片| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区三区四区免费观看| 22中文网久久字幕| 一本—道久久a久久精品蜜桃钙片| 日日爽夜夜爽网站| 免费在线观看成人毛片| 亚洲国产毛片av蜜桃av| 国产在线免费精品| 亚洲精品,欧美精品| 少妇的逼好多水| 日韩免费高清中文字幕av| 高清毛片免费看| 色网站视频免费| 少妇丰满av| 最近中文字幕高清免费大全6| 国产乱来视频区| 99久久精品一区二区三区| 高清黄色对白视频在线免费看 | 欧美丝袜亚洲另类| 高清视频免费观看一区二区| 日韩伦理黄色片| 在线观看免费高清a一片| 亚洲av不卡在线观看| 亚洲在久久综合| 亚洲欧美精品专区久久| 国产精品女同一区二区软件| 在线观看www视频免费| 岛国毛片在线播放| 亚洲欧洲国产日韩| 99九九在线精品视频 | 亚洲欧美精品专区久久| 日本爱情动作片www.在线观看| 国产精品伦人一区二区| 99国产精品免费福利视频| 欧美日韩视频精品一区| 国精品久久久久久国模美| 亚洲精品乱码久久久v下载方式| 国产男女超爽视频在线观看| 亚洲欧美成人综合另类久久久| 久久国产亚洲av麻豆专区| 亚洲国产成人一精品久久久| 亚洲中文av在线| 国内揄拍国产精品人妻在线| 国产精品欧美亚洲77777| 国精品久久久久久国模美| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 日韩大片免费观看网站| 乱系列少妇在线播放| 欧美人与善性xxx| 中文欧美无线码| av国产精品久久久久影院| 丝袜在线中文字幕| 亚洲性久久影院| 国产精品偷伦视频观看了| 免费观看无遮挡的男女| 国产在线一区二区三区精| 全区人妻精品视频| 免费久久久久久久精品成人欧美视频 | xxx大片免费视频| 好男人视频免费观看在线| 熟妇人妻不卡中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品专区欧美| 日韩亚洲欧美综合| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 久久精品久久精品一区二区三区| 免费黄频网站在线观看国产| 制服丝袜香蕉在线| 久久久国产一区二区| 国产美女午夜福利| 午夜视频国产福利| 亚洲精品视频女| 我要看日韩黄色一级片| 免费看日本二区| 亚洲av二区三区四区| 国产精品人妻久久久影院| 99久久精品国产国产毛片| av专区在线播放| 男人舔奶头视频| 国产免费视频播放在线视频| 久久精品国产鲁丝片午夜精品| 日本免费在线观看一区| 国内揄拍国产精品人妻在线| 美女内射精品一级片tv| 日本91视频免费播放| 成人毛片a级毛片在线播放| 久久久久人妻精品一区果冻| 香蕉精品网在线| 欧美日韩综合久久久久久| 国产成人精品婷婷| 丰满饥渴人妻一区二区三| 尾随美女入室| 哪个播放器可以免费观看大片| 日日摸夜夜添夜夜爱| 伊人久久国产一区二区| 久久人妻熟女aⅴ| 国产男女内射视频| 建设人人有责人人尽责人人享有的| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄| 国产淫片久久久久久久久| 中国美白少妇内射xxxbb| 日日啪夜夜撸| 亚洲欧洲精品一区二区精品久久久 | 国产免费一区二区三区四区乱码| 观看免费一级毛片| 亚洲av国产av综合av卡| 午夜91福利影院| 免费看不卡的av| 精华霜和精华液先用哪个| 亚洲高清免费不卡视频| 精品一区二区三卡| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩卡通动漫| 极品少妇高潮喷水抽搐| 在线看a的网站| 久久影院123| 18禁动态无遮挡网站| 熟女电影av网| 精品久久国产蜜桃| 国产亚洲一区二区精品| 日本免费在线观看一区| 日本爱情动作片www.在线观看| 日韩精品免费视频一区二区三区 | 你懂的网址亚洲精品在线观看| 少妇的逼好多水| 欧美丝袜亚洲另类| 亚洲国产精品一区二区三区在线| 伦理电影免费视频| 久久精品国产a三级三级三级| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 麻豆精品久久久久久蜜桃| 国产成人精品福利久久| 高清av免费在线| 亚洲av在线观看美女高潮| 美女大奶头黄色视频| 91久久精品国产一区二区成人| 观看免费一级毛片| 国产一区二区三区av在线| 久久国内精品自在自线图片| 人人妻人人看人人澡| 91aial.com中文字幕在线观看| 乱系列少妇在线播放| 亚洲国产精品国产精品| 免费观看av网站的网址| 在线观看国产h片| 成年女人在线观看亚洲视频| 菩萨蛮人人尽说江南好唐韦庄| 特大巨黑吊av在线直播| 中文欧美无线码| 狂野欧美白嫩少妇大欣赏| 国产女主播在线喷水免费视频网站| av天堂中文字幕网| 国产一区二区在线观看日韩| 在现免费观看毛片| 欧美人与善性xxx| 99热这里只有是精品在线观看| 国产男女超爽视频在线观看| 我要看日韩黄色一级片| 欧美精品高潮呻吟av久久| 亚洲国产av新网站| 久久久久久人妻| 久久免费观看电影| 最近手机中文字幕大全| 国产黄片美女视频| av国产精品久久久久影院| 久热这里只有精品99| 亚洲伊人久久精品综合| 麻豆精品久久久久久蜜桃| 欧美xxⅹ黑人| 天美传媒精品一区二区| 久久毛片免费看一区二区三区| 色哟哟·www| 夜夜爽夜夜爽视频| 日本vs欧美在线观看视频 | 亚洲av日韩在线播放| 久久6这里有精品| 极品少妇高潮喷水抽搐| 亚洲av欧美aⅴ国产| av不卡在线播放| 夜夜看夜夜爽夜夜摸| 国产在视频线精品| 久久久久久久久久久免费av| 大片电影免费在线观看免费| 精品亚洲成a人片在线观看| 视频区图区小说| 亚洲美女搞黄在线观看| 日本午夜av视频| 少妇熟女欧美另类| 免费观看av网站的网址| 日本爱情动作片www.在线观看| 中文在线观看免费www的网站| 日韩中字成人| 最近的中文字幕免费完整| 欧美xxⅹ黑人| 午夜福利影视在线免费观看| 亚洲经典国产精华液单| 亚洲欧美日韩卡通动漫| 亚洲在久久综合| av国产精品久久久久影院| 亚洲欧美日韩东京热| 精品少妇久久久久久888优播| 色5月婷婷丁香| 亚洲国产精品一区三区| 赤兔流量卡办理| 男女国产视频网站| 国产色婷婷99| 最近最新中文字幕免费大全7| 欧美成人午夜免费资源| av福利片在线| 香蕉精品网在线| 日本色播在线视频| 久久99热6这里只有精品| 最新的欧美精品一区二区| 婷婷色综合www| 亚洲久久久国产精品| 久久综合国产亚洲精品| 人人妻人人看人人澡| 午夜久久久在线观看| 国产在线一区二区三区精| 少妇被粗大猛烈的视频| 狂野欧美激情性bbbbbb| 肉色欧美久久久久久久蜜桃| 美女内射精品一级片tv| 精品人妻熟女毛片av久久网站| 亚洲精品久久午夜乱码| 精品一区二区三区视频在线| 久久精品国产亚洲av涩爱| 色吧在线观看| 十八禁高潮呻吟视频 | 美女内射精品一级片tv| 黄色一级大片看看| 国产高清国产精品国产三级| 成年女人在线观看亚洲视频| 亚州av有码| 国产日韩欧美亚洲二区| 日韩伦理黄色片| 超碰97精品在线观看| 伊人亚洲综合成人网| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 国模一区二区三区四区视频| .国产精品久久| 国产日韩一区二区三区精品不卡 | 亚洲四区av| 欧美日韩一区二区视频在线观看视频在线| 91精品伊人久久大香线蕉| 久久久久久久久久成人| 91久久精品国产一区二区三区| 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 亚洲不卡免费看| 看非洲黑人一级黄片| 成人二区视频| 免费人妻精品一区二区三区视频| 国产成人aa在线观看| 91久久精品电影网| 国产淫片久久久久久久久| 国产白丝娇喘喷水9色精品| 在线观看免费高清a一片| 国产在线视频一区二区| 校园人妻丝袜中文字幕| 国产毛片在线视频| 99九九线精品视频在线观看视频| 中文资源天堂在线| 色视频在线一区二区三区| 国产一区二区三区综合在线观看 | 久久久久网色| 人体艺术视频欧美日本| 卡戴珊不雅视频在线播放| 日韩精品免费视频一区二区三区 | 在线观看免费高清a一片| 亚洲精品一二三| 在线观看国产h片| 日日摸夜夜添夜夜添av毛片| 激情五月婷婷亚洲| 91精品一卡2卡3卡4卡| 久久午夜综合久久蜜桃| 国产精品不卡视频一区二区| 搡老乐熟女国产| 中国美白少妇内射xxxbb| 久久午夜综合久久蜜桃| 18+在线观看网站| av天堂久久9| 国产伦在线观看视频一区| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| av卡一久久| 国产免费福利视频在线观看| 少妇熟女欧美另类| 中文欧美无线码| 丝瓜视频免费看黄片| 最近最新中文字幕免费大全7| av线在线观看网站| 日韩大片免费观看网站| 国产成人精品一,二区| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 日韩av免费高清视频| 少妇被粗大的猛进出69影院 | 欧美日韩国产mv在线观看视频| 美女脱内裤让男人舔精品视频| 王馨瑶露胸无遮挡在线观看| 亚洲av二区三区四区| 国语对白做爰xxxⅹ性视频网站| 永久网站在线| 国产成人aa在线观看| 午夜日本视频在线| 国产精品一二三区在线看| 亚洲美女黄色视频免费看| 九九爱精品视频在线观看| 国内揄拍国产精品人妻在线| 十八禁高潮呻吟视频 | 熟妇人妻不卡中文字幕| 精品久久久久久电影网|