• <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| 午夜福利视频精品| 一边摸一边做爽爽视频免费| 国产xxxxx性猛交| 久久久国产精品麻豆| 叶爱在线成人免费视频播放| 国产精品成人在线| 最近最新中文字幕免费大全7| 少妇的丰满在线观看| 国产精品免费大片| 黄片无遮挡物在线观看| 免费女性裸体啪啪无遮挡网站| 老鸭窝网址在线观看| 看非洲黑人一级黄片| 成人三级做爰电影| 中文字幕高清在线视频| 欧美亚洲日本最大视频资源| 男女无遮挡免费网站观看| 97在线人人人人妻| 午夜91福利影院| 中文字幕av电影在线播放| 毛片一级片免费看久久久久| 下体分泌物呈黄色| 亚洲av在线观看美女高潮| 不卡视频在线观看欧美| 色精品久久人妻99蜜桃| 日韩伦理黄色片| 欧美亚洲 丝袜 人妻 在线| 国产成人啪精品午夜网站| av卡一久久| 一二三四在线观看免费中文在| 伊人久久大香线蕉亚洲五| 少妇被粗大的猛进出69影院| 免费女性裸体啪啪无遮挡网站| 国产精品成人在线| 欧美日韩一级在线毛片| 只有这里有精品99| 亚洲精品aⅴ在线观看| 满18在线观看网站| 极品人妻少妇av视频| 欧美久久黑人一区二区| 亚洲欧美一区二区三区国产| 日韩 欧美 亚洲 中文字幕| 啦啦啦 在线观看视频| 欧美在线一区亚洲| 各种免费的搞黄视频| 综合色丁香网| 韩国av在线不卡| 久久久国产欧美日韩av| 成人亚洲精品一区在线观看| 欧美日韩视频精品一区| 精品一区在线观看国产| 99久国产av精品国产电影| 飞空精品影院首页| 熟女av电影| 在线看a的网站| 欧美日韩一级在线毛片| xxxhd国产人妻xxx| 九色亚洲精品在线播放| 国产在线一区二区三区精| 精品亚洲乱码少妇综合久久| 免费少妇av软件| 精品少妇黑人巨大在线播放| 久久99精品国语久久久| 男女国产视频网站| 在线 av 中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清不卡午夜福利| 别揉我奶头~嗯~啊~动态视频 | 久久久久国产一级毛片高清牌| 亚洲成色77777| 午夜免费男女啪啪视频观看| 国产精品.久久久| 婷婷色麻豆天堂久久| 狂野欧美激情性xxxx| 国精品久久久久久国模美| 在线 av 中文字幕| 国产免费又黄又爽又色| av.在线天堂| 亚洲美女黄色视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 中文天堂在线官网| 免费不卡黄色视频| 在线观看三级黄色| 国产精品一区二区在线不卡| 丝袜脚勾引网站| 999久久久国产精品视频| 蜜桃在线观看..| 亚洲国产精品一区二区三区在线| 操美女的视频在线观看| 国产免费一区二区三区四区乱码| 性色av一级| 日韩视频在线欧美| 亚洲av日韩在线播放| 亚洲av成人精品一二三区| 亚洲情色 制服丝袜| 又黄又粗又硬又大视频| 97人妻天天添夜夜摸| 另类精品久久| 国产99久久九九免费精品| svipshipincom国产片| 精品国产一区二区久久| 国产黄频视频在线观看| 亚洲av国产av综合av卡| 欧美人与性动交α欧美精品济南到| 精品一区二区三区四区五区乱码 | 精品国产一区二区三区久久久樱花| 免费在线观看完整版高清| 久久精品熟女亚洲av麻豆精品| 亚洲成人一二三区av| 亚洲四区av| 亚洲欧美一区二区三区久久| 91精品伊人久久大香线蕉| 丰满少妇做爰视频| 久久午夜综合久久蜜桃| 成人18禁高潮啪啪吃奶动态图| 一二三四中文在线观看免费高清| 视频在线观看一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品免费福利视频| 国产免费又黄又爽又色| 激情视频va一区二区三区| 日日摸夜夜添夜夜爱| 女性被躁到高潮视频| 999精品在线视频| 亚洲av在线观看美女高潮| 亚洲熟女毛片儿| 我要看黄色一级片免费的| 天天躁日日躁夜夜躁夜夜| 黑人猛操日本美女一级片| 中文字幕人妻丝袜制服| 久久久国产精品麻豆| 女人久久www免费人成看片| 无遮挡黄片免费观看| netflix在线观看网站| 久久综合国产亚洲精品| 日日爽夜夜爽网站| 日日爽夜夜爽网站| 国产高清国产精品国产三级| 亚洲精品久久午夜乱码| 久久青草综合色| 久久99热这里只频精品6学生| 在线观看国产h片| 热re99久久国产66热| 午夜福利视频精品| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美清纯卡通| 日韩欧美精品免费久久| 亚洲欧美一区二区三区久久| 亚洲欧美一区二区三区国产| 啦啦啦啦在线视频资源| 亚洲精品国产av蜜桃| 国产无遮挡羞羞视频在线观看| 丝袜美足系列| 99久久精品国产亚洲精品| 1024视频免费在线观看| 香蕉丝袜av| 日韩精品免费视频一区二区三区| 亚洲精品中文字幕在线视频| 国产无遮挡羞羞视频在线观看| 亚洲三区欧美一区| 国产精品久久久久成人av| 国产亚洲av高清不卡| 建设人人有责人人尽责人人享有的| 99久国产av精品国产电影| 亚洲 欧美一区二区三区| 亚洲欧洲日产国产| 中文精品一卡2卡3卡4更新| 人妻人人澡人人爽人人| 亚洲欧美色中文字幕在线| 亚洲一码二码三码区别大吗| 丰满乱子伦码专区| 免费黄色在线免费观看| 99re6热这里在线精品视频| 一区在线观看完整版| 男人操女人黄网站| 国产男女内射视频| 嫩草影视91久久| a级片在线免费高清观看视频| 精品国产乱码久久久久久男人| 十分钟在线观看高清视频www| 女人爽到高潮嗷嗷叫在线视频| 婷婷色av中文字幕| 性高湖久久久久久久久免费观看| 一区二区av电影网| 久久影院123| 999久久久国产精品视频| 91精品三级在线观看| 婷婷色综合大香蕉| 男男h啪啪无遮挡| 999久久久国产精品视频| 成年人免费黄色播放视频| 下体分泌物呈黄色| 黑丝袜美女国产一区| 久久性视频一级片| 丰满乱子伦码专区| 我的亚洲天堂| 国产精品熟女久久久久浪| 久久人人爽人人片av| 下体分泌物呈黄色| 九九爱精品视频在线观看| 国产成人精品在线电影| 精品久久久久久电影网| 一边摸一边抽搐一进一出视频| 一区福利在线观看| 如日韩欧美国产精品一区二区三区| 日韩成人av中文字幕在线观看| 如何舔出高潮| 欧美激情高清一区二区三区 | 久久久久精品性色| 国产精品国产av在线观看| 汤姆久久久久久久影院中文字幕| www.精华液| 最近的中文字幕免费完整| 亚洲男人天堂网一区| 亚洲天堂av无毛| 日日摸夜夜添夜夜爱| 日日啪夜夜爽| 亚洲一码二码三码区别大吗| 色视频在线一区二区三区| 亚洲国产欧美日韩在线播放| 国产精品一区二区在线观看99| 色婷婷av一区二区三区视频| 9191精品国产免费久久| 蜜桃在线观看..| 国产免费福利视频在线观看| 色婷婷av一区二区三区视频| 日韩一区二区三区影片| 亚洲国产毛片av蜜桃av| 国产女主播在线喷水免费视频网站| 国产亚洲最大av| 国产免费一区二区三区四区乱码| 国产在线视频一区二区| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合www| 久久狼人影院| 精品国产乱码久久久久久小说| 热99国产精品久久久久久7| 欧美黑人精品巨大| av免费观看日本| 热re99久久国产66热| 国产一区有黄有色的免费视频| 日韩人妻精品一区2区三区| 国产精品女同一区二区软件| 十八禁高潮呻吟视频| 天天操日日干夜夜撸| 成人漫画全彩无遮挡| 999久久久国产精品视频| 免费看不卡的av| 人成视频在线观看免费观看| 精品国产一区二区久久| 久久久久网色| 国产男女超爽视频在线观看| 中文字幕最新亚洲高清| 丁香六月天网| av网站免费在线观看视频| 天堂8中文在线网| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 国产成人a∨麻豆精品| 亚洲成人免费av在线播放| 永久免费av网站大全| 91老司机精品| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 欧美日韩精品网址| 在线精品无人区一区二区三| 国产精品亚洲av一区麻豆 | 日本欧美视频一区| 亚洲av中文av极速乱| 免费观看性生交大片5| 亚洲成人免费av在线播放| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 在线观看一区二区三区激情| 中文欧美无线码| 亚洲国产欧美网| 日韩一区二区视频免费看| 国产精品成人在线| 在现免费观看毛片| 老汉色av国产亚洲站长工具| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 看十八女毛片水多多多| 老鸭窝网址在线观看| 大香蕉久久成人网| 国产深夜福利视频在线观看| 亚洲成人免费av在线播放| 欧美老熟妇乱子伦牲交| 国产精品一国产av| 毛片一级片免费看久久久久| 天天影视国产精品| 中国三级夫妇交换| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 九九爱精品视频在线观看| 两性夫妻黄色片| 少妇人妻精品综合一区二区| 男女免费视频国产| 一本大道久久a久久精品| 九草在线视频观看| 日本爱情动作片www.在线观看| 午夜福利一区二区在线看| 大码成人一级视频| 国产免费又黄又爽又色| 亚洲综合色网址| 亚洲av在线观看美女高潮| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 久久人人爽人人片av| 精品福利永久在线观看| 免费日韩欧美在线观看| 国产伦人伦偷精品视频| av不卡在线播放| 午夜福利在线免费观看网站| 久久久久久久精品精品| 久久久久精品人妻al黑| 亚洲色图综合在线观看| 国产在线免费精品| 91精品伊人久久大香线蕉| 亚洲国产日韩一区二区| 亚洲 欧美一区二区三区| 成人午夜精彩视频在线观看| 日韩不卡一区二区三区视频在线| 久久久久网色| 亚洲欧美精品自产自拍| 看免费av毛片| 国产精品香港三级国产av潘金莲 | 日韩av不卡免费在线播放| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 伊人亚洲综合成人网| 啦啦啦视频在线资源免费观看| 日韩一区二区三区影片| 天天影视国产精品| av在线老鸭窝| 热99久久久久精品小说推荐| 国产免费福利视频在线观看| 亚洲成色77777| 麻豆av在线久日| 一级片'在线观看视频| 久久99热这里只频精品6学生| 久热爱精品视频在线9| 成年av动漫网址| 久久久精品免费免费高清| 精品亚洲成a人片在线观看| 99国产精品免费福利视频| 欧美日韩一区二区视频在线观看视频在线| 十八禁人妻一区二区| bbb黄色大片| 日韩人妻精品一区2区三区| 人人妻人人澡人人看| 啦啦啦视频在线资源免费观看| 在线观看www视频免费| 成年av动漫网址| 在线精品无人区一区二区三| 熟女av电影| 亚洲欧美精品自产自拍| 国产精品亚洲av一区麻豆 | 久久精品久久久久久噜噜老黄| 制服诱惑二区| 中国三级夫妇交换| 看免费av毛片| 丝袜在线中文字幕| 黄色一级大片看看| 91aial.com中文字幕在线观看| 国产亚洲欧美精品永久| 色视频在线一区二区三区| 亚洲一区中文字幕在线| 99国产精品免费福利视频| 天天躁夜夜躁狠狠久久av| xxx大片免费视频| 久久久久久久精品精品| 日本欧美国产在线视频| 男人爽女人下面视频在线观看| 少妇 在线观看| 麻豆乱淫一区二区| av不卡在线播放| 欧美日韩精品网址| 久久久久久人妻| 成年动漫av网址| 欧美日韩亚洲综合一区二区三区_| 国产精品欧美亚洲77777| 九九爱精品视频在线观看| 色精品久久人妻99蜜桃| 精品少妇内射三级| 国产成人午夜福利电影在线观看| 一级毛片 在线播放| 老司机影院成人| 天天影视国产精品| 女性生殖器流出的白浆| 免费观看a级毛片全部| 中文字幕色久视频| 一边亲一边摸免费视频| 校园人妻丝袜中文字幕| 亚洲欧美一区二区三区久久| 建设人人有责人人尽责人人享有的| av线在线观看网站| 亚洲中文av在线| 亚洲熟女精品中文字幕| 精品一区二区三区av网在线观看 | 在线天堂最新版资源| 少妇被粗大猛烈的视频| www.精华液| 国产成人免费无遮挡视频| 国产精品久久久久成人av| 少妇人妻精品综合一区二区| 久久毛片免费看一区二区三区| 香蕉丝袜av| 看免费成人av毛片| 婷婷色麻豆天堂久久| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 亚洲国产看品久久| 狠狠婷婷综合久久久久久88av| 精品一区二区三区四区五区乱码 | 欧美亚洲日本最大视频资源| 深夜精品福利| 99精国产麻豆久久婷婷| 国产精品亚洲av一区麻豆 | 高清av免费在线| 精品少妇久久久久久888优播| 激情视频va一区二区三区| 80岁老熟妇乱子伦牲交| 女性生殖器流出的白浆| av在线老鸭窝| kizo精华| 亚洲精品美女久久av网站| 青草久久国产| 国产高清国产精品国产三级| 熟女少妇亚洲综合色aaa.| 中文字幕av电影在线播放| 女人被躁到高潮嗷嗷叫费观| 精品一品国产午夜福利视频| 亚洲精品aⅴ在线观看| 精品国产一区二区久久| av国产精品久久久久影院| av线在线观看网站| 国产高清不卡午夜福利| 久久久国产一区二区| 欧美日韩精品网址| 国产成人精品久久久久久| av一本久久久久| 亚洲一码二码三码区别大吗| 亚洲欧洲国产日韩| 中文字幕最新亚洲高清| 精品人妻一区二区三区麻豆| 一级,二级,三级黄色视频| 色播在线永久视频| 国产成人午夜福利电影在线观看| 中文字幕色久视频| 女人被躁到高潮嗷嗷叫费观| 麻豆av在线久日| 免费在线观看视频国产中文字幕亚洲 | 欧美97在线视频| 精品人妻熟女毛片av久久网站| 男女国产视频网站| 国产精品香港三级国产av潘金莲 | 色综合欧美亚洲国产小说| 亚洲av成人不卡在线观看播放网 | 十八禁高潮呻吟视频| 中文字幕人妻熟女乱码| 久久人人爽人人片av| 大话2 男鬼变身卡| 久久久国产一区二区| 国产野战对白在线观看| 日本wwww免费看| 女的被弄到高潮叫床怎么办| 亚洲精品视频女| 国产老妇伦熟女老妇高清| 国产免费现黄频在线看| 国产亚洲最大av| 欧美97在线视频| 亚洲一级一片aⅴ在线观看| 十八禁高潮呻吟视频| 亚洲欧美成人综合另类久久久| 欧美日韩亚洲高清精品| 超色免费av| 精品少妇一区二区三区视频日本电影 | 久久人妻熟女aⅴ| 夜夜骑夜夜射夜夜干| 国产有黄有色有爽视频| 久久韩国三级中文字幕| 男女床上黄色一级片免费看| www.精华液| 成年动漫av网址| 精品午夜福利在线看| 狠狠精品人妻久久久久久综合| 日韩av不卡免费在线播放| 老司机在亚洲福利影院| 男女免费视频国产| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 9191精品国产免费久久| 色94色欧美一区二区| av卡一久久| 国产精品一区二区在线观看99| 国产免费视频播放在线视频| 国产在线视频一区二区| 99久久综合免费| 国产福利在线免费观看视频| 丝袜在线中文字幕| 在线观看三级黄色| 精品人妻熟女毛片av久久网站| 亚洲精品国产av成人精品| 日韩av在线免费看完整版不卡| 国产精品蜜桃在线观看| av网站在线播放免费| 不卡视频在线观看欧美| 日本欧美国产在线视频| 波多野结衣一区麻豆| 亚洲国产毛片av蜜桃av| 精品一区二区三区av网在线观看 | 国产日韩一区二区三区精品不卡| 亚洲成色77777| 最近2019中文字幕mv第一页| 亚洲av成人精品一二三区| 欧美黑人精品巨大| 91精品国产国语对白视频| 在线看a的网站| 美女视频免费永久观看网站| 国产成人欧美在线观看 | 大香蕉久久成人网| 国产亚洲av片在线观看秒播厂| 欧美中文综合在线视频| 看十八女毛片水多多多| 最近的中文字幕免费完整| 精品国产一区二区久久| 国产午夜精品一二区理论片| 在线观看免费高清a一片| 欧美xxⅹ黑人| 欧美 日韩 精品 国产| 成人三级做爰电影| 国产成人精品福利久久| 一边摸一边抽搐一进一出视频| 亚洲av国产av综合av卡| 青春草国产在线视频| 一区二区三区精品91| 亚洲三区欧美一区| 国产精品国产三级专区第一集| 亚洲精品第二区| 亚洲成人av在线免费| 无限看片的www在线观看| 大话2 男鬼变身卡| 丝袜脚勾引网站| 丝袜人妻中文字幕| 国产日韩欧美视频二区| 只有这里有精品99| 黄色 视频免费看| 久久这里只有精品19| 日韩欧美一区视频在线观看| 国产成人免费无遮挡视频| 黄频高清免费视频| 青春草国产在线视频| 爱豆传媒免费全集在线观看| 老汉色av国产亚洲站长工具| 一边摸一边抽搐一进一出视频| 久久99精品国语久久久| 一级爰片在线观看| 校园人妻丝袜中文字幕| 交换朋友夫妻互换小说| 老熟女久久久| 国产日韩一区二区三区精品不卡| 亚洲第一区二区三区不卡| 悠悠久久av| 2018国产大陆天天弄谢| 90打野战视频偷拍视频| 在线免费观看不下载黄p国产| bbb黄色大片| 亚洲精华国产精华液的使用体验| 最近手机中文字幕大全| 观看av在线不卡| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av高清一级| 久久久久久久大尺度免费视频| 日韩av免费高清视频| 午夜福利免费观看在线| 少妇人妻精品综合一区二区| 亚洲欧美成人综合另类久久久| kizo精华| 日本欧美国产在线视频| 国产欧美日韩综合在线一区二区| 久久毛片免费看一区二区三区| 午夜福利影视在线免费观看| 哪个播放器可以免费观看大片| 日韩电影二区| 91精品伊人久久大香线蕉| 亚洲精品久久久久久婷婷小说| 亚洲国产成人一精品久久久| 亚洲av男天堂| 两个人免费观看高清视频| 日韩成人av中文字幕在线观看| 九草在线视频观看| 亚洲欧美清纯卡通| 母亲3免费完整高清在线观看| 亚洲欧美激情在线| 日韩av不卡免费在线播放| 热re99久久精品国产66热6| 在线观看免费午夜福利视频| 亚洲美女搞黄在线观看| 国产一区亚洲一区在线观看| 久久久欧美国产精品| 亚洲色图综合在线观看| 国产精品av久久久久免费| 777久久人妻少妇嫩草av网站| 欧美精品一区二区免费开放| 亚洲伊人色综图| 两个人看的免费小视频| 人成视频在线观看免费观看| 亚洲少妇的诱惑av| av在线播放精品| 汤姆久久久久久久影院中文字幕| 久久国产精品男人的天堂亚洲|