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

    粗糙度對(duì)高/低雷諾數(shù)跨聲壓氣機(jī)性能的影響

    2017-06-05 15:00:53楊榮菲楊小平向宏輝
    燃?xì)鉁u輪試驗(yàn)與研究 2017年2期
    關(guān)鍵詞:葉柵雷諾數(shù)激波

    楊榮菲,黃 進(jìn),楊小平,向宏輝,葛 寧

    (1.南京航空航天大學(xué)能源與動(dòng)力學(xué)院江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,南京210016;2.中國(guó)航發(fā)四川燃?xì)鉁u輪研究院航空發(fā)動(dòng)機(jī)高空模擬技術(shù)重點(diǎn)實(shí)驗(yàn)室,四川江油621703)

    粗糙度對(duì)高/低雷諾數(shù)跨聲壓氣機(jī)性能的影響

    楊榮菲1,黃 進(jìn)1,楊小平2,向宏輝2,葛 寧1

    (1.南京航空航天大學(xué)能源與動(dòng)力學(xué)院江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,南京210016;2.中國(guó)航發(fā)四川燃?xì)鉁u輪研究院航空發(fā)動(dòng)機(jī)高空模擬技術(shù)重點(diǎn)實(shí)驗(yàn)室,四川江油621703)

    以跨聲壓氣機(jī)Stage 35為研究對(duì)象,針對(duì)地面、20 km高空兩種雷諾數(shù)工況,數(shù)值研究了轉(zhuǎn)子壓力面、吸力面、整個(gè)葉片分別為光滑及5μm、20μm、45μm粗糙度時(shí)壓氣機(jī)的性能變化。結(jié)果表明:吸力面粗糙度較壓力面粗糙度對(duì)壓氣機(jī)性能的影響更大;粗糙度對(duì)低雷諾數(shù)壓氣機(jī)性能的影響小于高雷諾數(shù)壓氣機(jī);相較于粗糙度總是惡化高雷諾數(shù)壓氣機(jī)性能,在低雷諾數(shù)工況,小幅值粗糙度能改善壓氣機(jī)性能,而大幅值粗糙度惡化壓氣機(jī)性能。當(dāng)粗糙度為5μm時(shí),壓氣機(jī)峰值效率最大增量為0.79%。

    航空發(fā)動(dòng)機(jī);跨聲壓氣機(jī);表面粗糙度;低雷諾數(shù);高雷諾數(shù);氣動(dòng)性能;數(shù)值模擬

    1 引言

    航空發(fā)動(dòng)機(jī)長(zhǎng)期運(yùn)行過(guò)程中,腐蝕、磨損及污垢堆積引起葉片表面粗糙度大幅增加,壓氣機(jī)氣動(dòng)性能衰減。而飛行高度增加使壓氣機(jī)進(jìn)口雷諾數(shù)下降,也會(huì)惡化壓氣機(jī)性能。因此,有必要研究不同雷諾數(shù)下粗糙度對(duì)壓氣機(jī)氣動(dòng)性能的影響。

    粗糙度通常指加工或污染表面的幾何不平度,用中心線平均粗糙度Ra值表示,而流體力學(xué)粗糙度一般用當(dāng)量砂礫粗糙度ks描述。Bons[1]總結(jié)獲得ks=(2.0~10.0)Ra;Koch等[2]指出壓氣機(jī)污染粗糙度滿(mǎn)足關(guān)系式ks=6.2Ra,本文將采用此式計(jì)算。

    Back等[3]實(shí)驗(yàn)發(fā)現(xiàn),ks=12~850 μm的平面葉柵落后角及損失隨粗糙度的增加而增加。Gbadebo等[4]在靜子吸力面覆蓋Ra=25 μm砂紙,實(shí)驗(yàn)發(fā)現(xiàn)粗糙度位于前緣至吸力面峰值速度之間時(shí)壓氣機(jī)性能顯著降低,而位于峰值速度下游時(shí)其影響可以忽略。Chen等[5]在低速高負(fù)荷平面葉柵實(shí)驗(yàn)臺(tái)上,研究了葉片覆蓋不同面積及粗糙度砂紙(Ra=150、44 μm)的工況,發(fā)現(xiàn)粗糙度使分離流再附,大多數(shù)實(shí)驗(yàn)工況葉柵損失降低,與文獻(xiàn)[3-4]的結(jié)果矛盾,這是因?yàn)榇颂幋植诙雀笄颐黠@改變了葉型幾何所致。為獲得不同雷諾數(shù)下粗糙度對(duì)葉柵性能的影響,Leipold等[6]實(shí)驗(yàn)研究了Re=(3~10)×105工況下的高速高載荷光滑葉柵(Ra=0.70 μm)及粗糙葉柵(Ra=11.15 μm),發(fā)現(xiàn)葉柵損失隨雷諾數(shù)的增加而增加,粗糙葉柵較光滑葉柵損失更大,且隨著雷諾數(shù)的增加,粗糙葉柵損失在低雷諾數(shù)下增加較慢而在高雷諾數(shù)下增加劇烈,葉柵損失主要來(lái)源于粗糙度引起的層流分離泡減小及湍流分離增加。Back等[7]實(shí)驗(yàn)觀察了Re= (3.0~6.4)×105、Ra=0.38~2.89 μm,以及粗糙度位于壓力面、前緣、吸力面自前緣開(kāi)始至不同弦長(zhǎng)位置處的平面葉柵,發(fā)現(xiàn)壓力面及葉片前緣粗糙度對(duì)葉柵性能無(wú)影響,吸力面粗糙面積小于20%時(shí),葉柵損失基本不受雷諾數(shù)影響,吸力面粗糙面積大于20%時(shí),隨著粗糙面積的增加,損失隨雷諾數(shù)增加的速率變大。更進(jìn)一步,Ju等[8]采用LDV技術(shù),測(cè)量了平面葉柵吸力面前緣至2%弦長(zhǎng)處ks=100 μm時(shí)吸力面邊界層速度及葉柵損失。當(dāng)Re=(2.1~3.8)×105時(shí),前緣分離泡尺寸及湍流邊界層再附后的湍流水平降低,葉柵損失降低;而當(dāng)Re=(4.5~6.4)×105時(shí),前緣粗糙度降低了軸向正壓力梯度,使湍流邊界層提前分離,葉柵損失急劇增加。

    針對(duì)跨聲速壓氣機(jī),國(guó)內(nèi)外也開(kāi)展了粗糙度影響的研究。Suder等[9]在Rotor 37轉(zhuǎn)子葉片表面分別增加0.025 mm厚的光滑涂層(Ra=0.254~0.508 μm)及粗糙涂層(Ra=2.54~3.18 μm),實(shí)驗(yàn)發(fā)現(xiàn)粗糙度增加使壓氣機(jī)設(shè)計(jì)點(diǎn)效率、壓比分別下降了6%和9%,而葉片厚度增加導(dǎo)致壓氣機(jī)性能衰減僅為粗糙葉片引起壓氣機(jī)性能衰減值的一半,壓氣機(jī)性能衰減來(lái)源于激波與增厚邊界層作用引起的流動(dòng)堵塞增加和擴(kuò)散能力下降。陳紹文等[10]指出,增加壁面粗糙度對(duì)壓氣機(jī)穩(wěn)定工作范圍影響不大,而增加葉片厚度減小了穩(wěn)定工作范圍。以Stage 37為數(shù)值研究對(duì)象,Morini等[11]發(fā)現(xiàn)隨著轉(zhuǎn)子表面Ra從0.16 μm增至12.1 μm,壓氣機(jī)流量、壓比、效率下降,效率損失主要來(lái)源于流動(dòng)堵塞增加,粗糙度自身產(chǎn)生的損失對(duì)效率影響很小。Morini等[12]還發(fā)現(xiàn),轉(zhuǎn)子粗糙度對(duì)壓氣機(jī)性能的影響大于靜子,吸力面粗糙度的影響大于壓力面。Aldi等[13]數(shù)值觀察了轉(zhuǎn)子Ra從0.81 μm到6.45 μm展向線性變化工況,發(fā)現(xiàn)粗糙度分布改變了壓氣機(jī)葉片載荷分布。Chen等對(duì)Stage 35數(shù)值計(jì)算發(fā)現(xiàn),當(dāng)轉(zhuǎn)子Ra從8.06 μm增至16.12 μm時(shí),壓氣機(jī)性能衰減加劇[14];當(dāng)粗糙度靠近轉(zhuǎn)子端部時(shí)性能惡化程度增加,葉片前緣附近粗糙度較尾緣附近粗糙度對(duì)性能影響大,某些位置的粗糙度能一定程度改善壓氣機(jī)性能[15]。上述粗糙度研究都是針對(duì)地面高雷諾數(shù)工況,而粗糙度對(duì)低雷諾數(shù)跨聲壓氣機(jī)性能影響的研究很少。Chen等[14]數(shù)值觀察了轉(zhuǎn)子表面Ra=16.12 μm時(shí)跨聲壓氣機(jī)的高空性能,相較于高雷諾數(shù)工況,粗糙度引起低雷諾數(shù)壓氣機(jī)性能衰減較少,但作者并未給出粗糙度影響低雷諾數(shù)壓氣機(jī)性能的機(jī)理。同時(shí),不同尺寸、位置的粗糙度對(duì)低雷諾數(shù)壓氣機(jī)性能的影響及機(jī)理,是否等同于粗糙度對(duì)高雷諾數(shù)壓氣機(jī)的影響,此方面研究也未見(jiàn)報(bào)道。

    對(duì)此,本文研究粗糙度大小、位置對(duì)其高/低雷諾數(shù)跨聲壓氣機(jī)氣動(dòng)性能的影響,以期為受污染跨聲壓氣機(jī)的性能預(yù)估提供一定支撐。

    2 數(shù)值計(jì)算方法

    以NASA Stage 35壓氣機(jī)作為研究對(duì)象。其設(shè)計(jì)轉(zhuǎn)速17 188 r/min,動(dòng)/靜葉片數(shù)及展弦比分別為36/46和1.19/1.26,轉(zhuǎn)子葉尖間隙為0.2%葉高,詳細(xì)幾何參數(shù)見(jiàn)文獻(xiàn)[16]。

    計(jì)算網(wǎng)格采用Numeca/Autogrid5軟件自動(dòng)生成,葉片表面及機(jī)匣/輪轂處固壁加密,第一層網(wǎng)格距壁面30 μm,對(duì)應(yīng)高雷諾數(shù)下y+≈2.53,總網(wǎng)格數(shù)58萬(wàn)。

    數(shù)值計(jì)算采用Numeca/Fine軟件,選用真實(shí)大氣、低雷諾數(shù)Spalart-Allmaras湍流模型。當(dāng)進(jìn)口湍流度高達(dá)5%時(shí),壓氣機(jī)內(nèi)流動(dòng)近似全湍流。Sonoda等[17]基于全湍流假設(shè),數(shù)值研究了一臺(tái)低雷諾數(shù)小型渦扇發(fā)動(dòng)機(jī)中的跨聲壓氣機(jī),數(shù)值結(jié)果與實(shí)驗(yàn)符合很好。因此,本文給定進(jìn)口湍流度為5%,不考慮轉(zhuǎn)捩問(wèn)題。

    對(duì)于粗糙壁面流動(dòng),如果粗糙度特征高度相較邊界層外層高度足夠小,則粗糙度對(duì)流動(dòng)的影響被限制在近壁粘性子層中,此時(shí)粗糙壁面邊界層外層流動(dòng)與光滑壁面外層流動(dòng)相似[18]。在此假設(shè)下,粗糙度對(duì)流動(dòng)的影響可通過(guò)簡(jiǎn)單修改壁面函數(shù)得到[19]:

    進(jìn)口按照實(shí)驗(yàn)條件給定總溫、總壓、軸向進(jìn)氣,固壁絕熱無(wú)滑移,單通道采用周期性邊界,轉(zhuǎn)靜交界面采用摻混平面,出口截面給定靜壓。改變出口靜壓獲得壓氣機(jī)從堵塞到失速的不同工作點(diǎn)。

    圖1給出Stage 35在100%設(shè)計(jì)轉(zhuǎn)速地面工況、光滑葉片的數(shù)值計(jì)算與實(shí)驗(yàn)特性對(duì)比。相較實(shí)驗(yàn)值,計(jì)算總壓比略微偏小、等熵效率符合較好。實(shí)驗(yàn)與計(jì)算的堵點(diǎn)流量分別為20.95 kg/s、20.77 kg/s,相對(duì)誤差0.86%;實(shí)驗(yàn)與計(jì)算的最大效率分別為84.5%、84.2%,相對(duì)誤差0.35%。由于最大誤差小于3%,且本文針對(duì)粗糙度開(kāi)展對(duì)比性研究,不同方案的網(wǎng)格及數(shù)值方法相同,方案對(duì)比過(guò)程中數(shù)值計(jì)算誤差可以消除,故認(rèn)為上述數(shù)值方法可行。

    3 粗糙度研究方案

    以Stage 35設(shè)計(jì)點(diǎn)平均進(jìn)氣速度、轉(zhuǎn)子中徑弦長(zhǎng)為特征參數(shù),計(jì)算海平面及高空壓氣機(jī)進(jìn)口雷諾數(shù)分別為1.088×106、0.860×105,分別對(duì)應(yīng)高、低雷諾數(shù)工況。

    選取ks分別為5 μm、20 μm、45 μm。Roberts等[20]指出,轉(zhuǎn)子葉片經(jīng)工業(yè)標(biāo)準(zhǔn)拋光后ks=3.10~3.70 μm,精拋光后ks=0.775 μm。Syverud等[21]對(duì)GE J88-13發(fā)動(dòng)機(jī)中的壓氣機(jī)進(jìn)行鹽水吞入實(shí)驗(yàn),發(fā)現(xiàn)葉片表面深度污染后粗糙度尺寸ks/c=(0.06~1.10)×10-3,以Stage 35轉(zhuǎn)子葉中弦長(zhǎng)為特征參數(shù),相當(dāng)于轉(zhuǎn)子表面粗糙度ks=3.30~61.16 μm。Melino等[22]通過(guò)數(shù)值計(jì)算及理論推導(dǎo),獲得了壓氣機(jī)工作時(shí)間與葉片表面粗糙度的關(guān)系,表明壓氣機(jī)工作1 650 h后葉片表面ks=30 μm。因此,本文選取的5 μm可看作葉片表面低污染,而20 μm、45 μm分別對(duì)應(yīng)葉片表面中度污染和深度污染。

    表1 不同雷諾數(shù)下的數(shù)值計(jì)算方案Table 1 Calculation cases at high/low Reynolds number

    文獻(xiàn)[12]認(rèn)為,粗糙轉(zhuǎn)子對(duì)壓氣機(jī)氣動(dòng)性能的影響大于粗糙靜子。故本文僅在轉(zhuǎn)子吸力面、壓力面及整個(gè)葉片覆蓋粗糙度。表1給出了粗糙度大小、位置及雷諾數(shù)改變所對(duì)應(yīng)的20種數(shù)值計(jì)算方案。表中:High-Re、Low-Re分別代表發(fā)動(dòng)機(jī)海平面和高空工況,PS、SS、All分別表示在壓力面、吸力面、整個(gè)葉片疊加粗糙度,5 μm等數(shù)值表征粗糙度水平。

    4 粗糙度對(duì)高雷諾數(shù)壓氣機(jī)性能的影響

    圖2為轉(zhuǎn)子表面粗糙及光滑時(shí)高雷諾數(shù)壓氣機(jī)特性計(jì)算結(jié)果,橫坐標(biāo)為無(wú)量綱流量??梢?jiàn),粗糙度大小及位置對(duì)壓氣機(jī)性能影響顯著。相較于葉片表面光滑情況,當(dāng)粗糙度位于壓力面時(shí),壓氣機(jī)的壓比基本不變,效率隨粗糙度增加而小幅下降;當(dāng)粗糙度位于吸力面時(shí),壓氣機(jī)壓比、效率隨粗糙度增加而急劇下降;吸力面粗糙度對(duì)壓氣機(jī)性能的影響比壓力面粗糙度的大,此結(jié)論與文獻(xiàn)[12]的結(jié)果一致。當(dāng)整個(gè)葉片粗糙時(shí),壓氣機(jī)性能為壓力面與吸力面粗糙度影響的綜合,并接近吸力面粗糙度下的壓氣機(jī)特性線。因此,下文僅針對(duì)壓力面/吸力面粗糙度工況,分析粗糙度對(duì)壓氣機(jī)性能影響的流動(dòng)機(jī)理。

    當(dāng)粗糙度位于壓力面時(shí),從圖3(a)中壓氣機(jī)峰值效率點(diǎn)50%葉高壁面靜壓分布曲線可看出,不同粗糙度下葉片載荷近似為常數(shù),導(dǎo)致圖2(a)中壓力面粗糙度對(duì)壓比影響很小。圖4對(duì)比了葉片光滑及壓力面粗糙度45 μm時(shí)壓氣機(jī)峰值效率點(diǎn)轉(zhuǎn)子出口截面相對(duì)總壓,圖中黑線為轉(zhuǎn)子尾緣。可見(jiàn),壓力面粗糙度明顯增加了壓力面尾跡寬度,這與文獻(xiàn)[7]中壓力面粗糙度通過(guò)加強(qiáng)湍流摻混來(lái)增加葉柵尾跡寬度的現(xiàn)象一致。尾跡寬度增加使轉(zhuǎn)子出口摻混損失增加,壓氣機(jī)效率小幅下降(圖2(b))。

    當(dāng)粗糙度位于吸力面時(shí),由圖3(b)可看出,隨著粗糙度的增加,壓力曲線所包圍的有效面積減小,對(duì)應(yīng)葉片載荷減小、壓氣機(jī)壓比降低(圖2(a));葉片前緣附近載荷逐漸增加、葉片后半部分載荷逐漸降低,葉片總載荷降低主要來(lái)源于葉片后半部分載荷的減少,這與文獻(xiàn)[12]獲得的粗糙葉片載荷降低主要來(lái)源于轉(zhuǎn)子前緣載荷減少的結(jié)論相矛盾,其原因?yàn)槲墨I(xiàn)[12]針對(duì)遠(yuǎn)離峰值效率點(diǎn)的等壓比工作點(diǎn)分析,與此處流動(dòng)并不相似。

    圖3(b)中葉片最大負(fù)荷點(diǎn)、吸力面激波位置隨粗糙度增加而逐漸前移,葉片載荷前移使葉片后半部分流動(dòng)在逆壓力梯度下更易分離。從圖5中壓氣機(jī)峰值效率點(diǎn)吸力面極限流線圖也可看出,相較光滑葉片,粗糙葉片輪轂分離點(diǎn)前移、分離范圍增大,導(dǎo)致相應(yīng)的分離損失增大。同時(shí),葉片前緣負(fù)荷增加使轉(zhuǎn)子前緣附近泄漏流增大,發(fā)展至尾緣附近葉尖泄漏渦面積增大,這從圖5中葉尖附近分離區(qū)增大也可看出,從而增加了葉尖泄漏損失。另一方面,葉片吸力面激波位置前移使激波-邊界層作用提前,導(dǎo)致粗糙葉片較光滑葉片尾緣邊界層厚度增加、損失增加。因此,吸力面粗糙度引起的三維分離損失、葉尖泄漏損失、激波-邊界層作用下的邊界層損失的增大,是造成高雷諾數(shù)壓氣機(jī)效率大幅下降的主要原因。

    5 粗糙度對(duì)低雷諾數(shù)壓氣機(jī)性能的影響

    圖6為葉片光滑及粗糙時(shí)低雷諾數(shù)跨聲壓氣機(jī)特性線的數(shù)值結(jié)果。與圖2高雷諾數(shù)情況對(duì)比可知,低雷諾數(shù)下粗糙度對(duì)壓氣機(jī)性能影響較小。粗糙葉片對(duì)壓氣機(jī)性能的影響,為壓力面/吸力面粗糙度影響的綜合。當(dāng)粗糙度位于壓力面時(shí),壓氣機(jī)壓比基本不隨粗糙度變化,壓氣機(jī)效率在低粗糙度下稍微上升,且效率改善量隨粗糙度增加而逐漸減??;當(dāng)粗糙度位于吸力面時(shí),壓氣機(jī)壓比及效率隨粗糙度的增加先增加后降低,并最終惡化壓氣機(jī)性能。此現(xiàn)象與文獻(xiàn)[14]中粗糙度惡化低雷諾數(shù)壓氣機(jī)性能的結(jié)論相矛盾,這是因?yàn)槲墨I(xiàn)[14]中的粗糙度為100 μm,遠(yuǎn)大于本文中的最大粗糙度,而本文計(jì)算方案中使壓氣機(jī)性能改善的粗糙度約為5 μm,大幅值粗糙度45 μm也使低雷諾數(shù)壓氣機(jī)性能惡化。

    圖7給出了壓力面/吸力面粗糙度時(shí)壓氣機(jī)峰值效率點(diǎn)葉中表面壓力分布。首先對(duì)比葉片光滑時(shí)高低雷諾數(shù)轉(zhuǎn)子葉中截面壓力分布(圖3及圖7):高雷諾數(shù)工況時(shí)葉片最大負(fù)荷位于50%弦長(zhǎng)附近,吸力面激波位于45%弦長(zhǎng)處、激波下游為湍流分離邊界層;而低雷諾數(shù)工況時(shí)葉片最大負(fù)荷位于葉片前緣,吸力面激波前移至25%軸向弦長(zhǎng)處,激波下游先出現(xiàn)分離泡然后發(fā)生湍流分離。這些流動(dòng)特征使低雷諾數(shù)壓氣機(jī)性能較高雷諾數(shù)工況急劇下降,峰值效率下降5.92%。因此,低雷諾數(shù)與高雷諾數(shù)壓氣機(jī)葉片表面壓力分布不同,造成粗糙度影響高/低雷諾數(shù)壓氣機(jī)性能的機(jī)理不同。

    當(dāng)粗糙度位于壓力面時(shí),圖7(a)中不同粗糙度下葉片總載荷近似不變,對(duì)應(yīng)壓氣機(jī)壓比對(duì)壓力面粗糙度不敏感(圖6(a))。圖8對(duì)比了不同粗糙度下壓氣機(jī)轉(zhuǎn)子出口截面相對(duì)總壓??梢?jiàn),相較于光滑葉片,當(dāng)壓力面粗糙度從5 μm增至45 μm時(shí),壓力面尾跡先變薄后增厚,與圖6(b)中壓氣機(jī)效率先小幅上升后小幅下降的趨勢(shì)相同。轉(zhuǎn)子尾跡寬度變化的原因分析如下:圖7(a)中葉片前半部分壓力面的逆壓力梯度使光滑葉片壓力面出現(xiàn)分離泡,而壓力面粗糙度增強(qiáng)了流動(dòng)抵抗逆壓力梯度的能力、使得分離泡減小,有利于壓力面?zhèn)任槽E變窄;另外,葉片后半部分粗糙度加強(qiáng)湍流摻混使壓力面尾跡變寬,這兩個(gè)因素共同影響尾跡寬度,進(jìn)而影響壓氣機(jī)效率。

    當(dāng)粗糙度位于吸力面時(shí),相較光滑葉片,圖7(b)中葉片平均載荷隨粗糙度的增加先增加后降低,與圖6(a)中壓比的變化特征相同。

    從圖7(b)可看出,相較于光滑葉片,小尺寸粗糙度如5 μm時(shí),吸力面激波位置后移,波后分離泡長(zhǎng)度減小;大尺寸粗糙度如20 μm時(shí),吸力面激波位置前移,波后分離泡長(zhǎng)度增加,且隨粗糙度增加變化較慢。由于跨聲轉(zhuǎn)子吸力面激波位置影響端壁三維分離區(qū)大小,激波位置后移時(shí)三維分離區(qū)減小,反之增大。圖9中吸力面表面極限流線也證明了此現(xiàn)象,故激波引起的三維分離損失是影響壓氣機(jī)效率變化的一個(gè)原因。另外,吸力面流動(dòng)在低雷諾數(shù)下易分離,小幅粗糙度有利于抵抗流動(dòng)分離,使葉片吸力面前半部分邊界層較薄(圖10(b)),激波-邊界層作用后邊界層厚度較光滑葉片減小。因此,吸力面小幅值粗糙度下壓氣機(jī)效率大幅增加歸因于三維分離區(qū)減小,吸力面邊界層變薄、分離泡縮短。而當(dāng)吸力面覆蓋較大尺寸粗糙度時(shí),一方面大尺寸粗糙度增加了邊界層湍流水平使吸力面激波上游的邊界層較厚,激波位置前移;另一方面當(dāng)邊界層厚度達(dá)到一定水平后邊界層厚度對(duì)粗糙度敏感性降低,并隨粗糙度增加緩慢增加,導(dǎo)致激波前移量較小,從而使激波相關(guān)的三維分離、激波-邊界層作用損失隨著粗糙度緩慢增加。故大幅值粗糙度下,壓氣機(jī)效率隨著粗糙度的增加而緩慢下降。

    為量化粗糙度大小、位置對(duì)高/低雷諾數(shù)跨聲壓氣機(jī)氣動(dòng)性能的影響,圖11給出了不同粗糙度下壓氣機(jī)峰值效率點(diǎn)效率、壓比變化量。其中,壓氣機(jī)效率、壓比變化量都以光滑葉片高雷諾數(shù)工況下的峰值效率工況為基準(zhǔn)進(jìn)行計(jì)算,效率變化量為絕對(duì)差值,壓比變化量為相對(duì)差值??梢?jiàn),高雷諾數(shù)工況下,粗糙度對(duì)壓氣機(jī)性能的影響不可忽略。當(dāng)整個(gè)葉片進(jìn)入深度污染、表面粗糙度達(dá)45 μm時(shí),壓氣機(jī)效率、壓比分別下降5.98%和6.15%。即使在少量污染、表面粗糙度為5 μm時(shí),壓氣機(jī)最大效率下降量仍高達(dá)1.79%、壓比下降1.41%。這意味著為保證跨聲壓氣機(jī)低空性能,必須在發(fā)動(dòng)機(jī)連續(xù)工作一段時(shí)間后對(duì)壓氣機(jī)進(jìn)行清洗。與高雷諾數(shù)工況不同,深度污染、大幅值粗糙度對(duì)低雷諾數(shù)壓氣機(jī)性能的影響很小,且存在臨界粗糙度,當(dāng)葉片表面粗糙度小于此臨界值(本文為20~45 μm)時(shí)壓氣機(jī)效率上升,反之壓氣機(jī)效率下降。在光滑葉片與臨界粗糙度之間,還存在著一個(gè)最佳表面粗糙度,使壓氣機(jī)效率最高。在本文低雷諾數(shù)算例中,最佳表面粗糙度為5 μm,對(duì)應(yīng)峰值效率最大增量為0.79%、壓比最大增加量為1.41%,這意味著低雷諾數(shù)工況葉片表面并不是越光滑越好。因此,為保證壓氣機(jī)在不同雷諾數(shù)下都有高的氣動(dòng)性能,需在整個(gè)葉片表面覆蓋5 μm以下的加工粗糙度。

    6 結(jié)論

    通過(guò)數(shù)值計(jì)算分析了粗糙度大小、位置對(duì)高/低雷諾數(shù)跨聲壓氣機(jī)氣動(dòng)性能的影響,得到如下結(jié)論:

    (1) 相較于壓力面粗糙度,吸力面粗糙度對(duì)壓氣機(jī)氣動(dòng)性能的影響更大。

    (2) 高雷諾數(shù)工況下壓氣機(jī)性能隨粗糙度的增加而顯著降低;低雷諾數(shù)工況下存在臨界粗糙度水平,粗糙度大于此值時(shí)壓氣機(jī)性能衰減,粗糙度小于此值時(shí)壓氣機(jī)性能改善。

    (3) 粗糙度對(duì)壓氣機(jī)壓比的影響來(lái)源于葉片載荷變化。壓力面粗糙度對(duì)壓氣機(jī)載荷分布無(wú)影響,而吸力面粗糙度對(duì)載荷分布影響顯著。

    (4) 壓力面粗糙度通過(guò)改變尾跡寬度來(lái)影響壓氣機(jī)效率,高雷諾數(shù)工況下尾跡寬度變化與壓力面湍流摻混相關(guān),而低雷諾數(shù)工況下尾跡寬度變化還受壓力面分離泡影響。吸力面粗糙度通過(guò)改變高雷諾數(shù)工況葉片前半部分負(fù)荷、低雷諾數(shù)工況葉片前半部分邊界層厚度來(lái)影響激波位置,進(jìn)而影響與激波位置相關(guān)的三維分離損失、激波-邊界層作用損失,同時(shí)葉片負(fù)荷變化還改變了葉尖泄漏損失,最終改變壓氣機(jī)效率。

    [1]Bons J P.A review of surface roughness effects in gas tur?bines[J]. Journal of Turbomachinery, 2010, 132:021004-1—16.

    [2]Koch C C,Smith L H.Loss sources and magnitudes in axi?al flow compressors[J].ASME J.Eng.Gas Turbines Power,1976,98(2):411—424.

    [3]Back S C,Sohn J H,Song S J.Impact of surface roughness on compressor cascade performance[J].Journal of Fluids Engineering,2010,132:064502-1—6.

    [4]Gbadebo S A,Hynes T P,Cumpsty N A.Influence of sur?face roughness on three-dimensional separation in axial compressor[R].ASME GT2004-53619,2004.

    [5]Chen S W,Xu H,Wang S T,et al.Experimental research of surface roughness effects on highly-loaded compressor cascade aerodynamics[J].Journal of Thermal Science,2014,23(4):307—314.

    [6]Leipold R,Boese M,F(xiàn)ottner L.The influence of technical surface roughness caused by precision forging on the flow around a highly loaded compressor cascade[J].Journal of Turbomachinery,2000,122:416—424.

    [7]Back S C,Hobson G V,Song S J,et al.Effects of Reynolds number and surface roughness magnitude and location on compressor cascade performance[J].Journal of Turboma?chinery,2012,134:051013-1—6.

    [8]Ju H I,Ju H S,Hobson G V,et al.Effect of leading edge roughness and Reynolds number on compressor profile loss [R].ASME GT2013-95487,2013.

    [9]Suder K L,Chima R V,Strazisar A J,et al.The Effect of adding roughness and thickness to a transonic axial com?pressor rotor[J].Journal of Turbomachinery,1995,117(4):491—505.

    [10]陳紹文,張 辰,石 慧,等,軸流壓氣機(jī)內(nèi)污垢沉積影響的數(shù)值研究[J].推進(jìn)技術(shù),2012,33(3):377—383.

    [11]Morini M,Pinelli M,Spina P R,et al.Computational fluid dynamics Simulation of fouling on axial compressor stages [J].Journal of Engineering for Gas Turbines and Power,2010,132:072401-1—10.

    [12]Morini M,Pinelli M,Spina P R,et al.Numerical analysis of the effects of non-uniform surface roughness on com?pressor stage performance[J].Journal of Engineering for Gas Turbines and Power,2011,133:072402-1—8.

    [13]Aldi N,Morini M,Pinelli M,et al.Performance evaluation of nonuniform ly fouled axial compressor stages by means of computational fluid dynamic analyses[J].Journal of Tur?bomachinery,2014,136:021016-1—11.

    [14]Chen S W,Zhang C,Shi H,et al.Study on the impact of fouling on axial compressor stage[R]. ASME GT2012-68041,2012.

    [15]Chen S W,Sun S J,Xu H,et al.Influence of local surface roughness of a rotor blade on performance of an axial com?pressor stage[R].ASME GT2013-94816,2013.

    [16]Reid L,Moore R D.Performance of single-stage axialflow transonic compressor with rotor and stator aspect ratios of 1.19 and 1.26,respectively,and with design pressure ratio of 1.82[R].NASA-TP-1338,1978.

    [17]Sonoda T,Arima T,Endicott G,et al.The effect of Reyn?olds number on a transonic swept fan OGV in a small tur?bofan engine[R].ASME GT2012-69956,2012.

    [18]Bons J P,Christensen K T.A comparison of real and simu?lated surface roughness characterizations[R]. AIAA 2007-3997,2007.

    [19]ANSYS CFX-Solver Theory Guide[M].Release 13.0.AN?SYS,Inc,2010:129—131.

    [20]Roberts W B,Thorp S A,Prahst P S,et al.The effect of ul?trapolish on a transonic axial rotor[J].Journal of Turboma?chinery,2013,135:011001-1—6.

    [21]Syverud E,Brekke O,Bakken L E.Axial compressor dete?rioration caused by saltwater ingestion[J].Journal of Tur?bomachinery,2007,129:119—127.

    [22]Melino F,Morini M,Peretto A,et al.Compressor fouling modeling:relationship between computational roughness and gas turbine operation time[J].Journal of Engineering for Gas Turbine and Power,2012,134:052401-1—8.

    Effect of roughness on high/low Reynolds num ber transonic com p ressor perform ance

    YANG Rong-fei1,HUANG Jin1,YANG Xiao-ping2,XIANG Hong-hui2,GE Ning1
    (1.Jiangsu Province Key Laboratory of Aerospace Power Systems,College of Energy and Power Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China;2.Key Laboratory on Aero-engine Altitude Simulation Technology,AECC Sichuan Gas Turbine Establishment,Jiangyou 621703,China)

    The NASA Stage 35 was used as the case study geometry.Compressor performance curves and flow structure were obtained through numerical simulation with smooth rotor and roughened rotor at two dif?ferent Reynolds number.The roughened rotor contained cases of various roughness of 5 μm,20 μm and 45 μm imposed on suction surface,pressure surface and entire blade respectively.And the high Reynolds num?ber condition was correspond to sea level while the low Reynolds number condition was obtained at high al?titude of 20 km.The results show that suction surface roughness has much more influence on compressor performance than pressure surface roughness.And the compressor performance changes greater at high Reynolds number than that at low Reynolds number under the influence of surface roughness.At high Reyn?olds number,the compressor performance usually decreases with increasing roughness.However,at low Reynolds number,the compressor performance can improve with small surface roughness while it worsens at high surface roughness.When the compressor rotor is roughened with 5 μm,the maximum increment of compressor peak efficiency reaches 0.79%.

    aero-engine;transonic compressor;surface roughness;low Reynolds number;high Reynolds number;aerodynamic performance;numerical simulation

    V231.3 文獻(xiàn)標(biāo)志碼:A

    :1672-2620(2017)02-0037-08

    2016-10-20;

    :2017-04-09

    中國(guó)航空科學(xué)基金(20141052012)

    楊榮菲(1982-),女,湖北武漢人,講師,博士,主要從事葉輪機(jī)械氣動(dòng)熱力學(xué)研究。

    猜你喜歡
    葉柵雷諾數(shù)激波
    亞聲速壓氣機(jī)平面葉柵及其改型的吹風(fēng)試驗(yàn)
    一種基于聚類(lèi)分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計(jì)研究
    民機(jī)高速風(fēng)洞試驗(yàn)的阻力雷諾數(shù)效應(yīng)修正
    超音速透平葉柵造型設(shè)計(jì)及分析
    一进一出抽搐gif免费好疼| 老熟妇乱子伦视频在线观看| 高潮久久久久久久久久久不卡| 又黄又爽又免费观看的视频| 亚洲精品日韩av片在线观看 | 99久久综合精品五月天人人| 熟妇人妻久久中文字幕3abv| 精华霜和精华液先用哪个| 亚洲无线在线观看| 成年人黄色毛片网站| 日韩欧美在线二视频| 内射极品少妇av片p| 国产精品 国内视频| 亚洲精华国产精华精| 最新在线观看一区二区三区| 欧美乱色亚洲激情| 免费看a级黄色片| 亚洲精品国产精品久久久不卡| 桃红色精品国产亚洲av| 国产成人av激情在线播放| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片| 日本免费a在线| 国内精品久久久久久久电影| 91在线观看av| 在线视频色国产色| 一个人免费在线观看电影| 国产中年淑女户外野战色| 免费在线观看日本一区| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| av女优亚洲男人天堂| 两个人看的免费小视频| 人妻夜夜爽99麻豆av| 亚洲av免费在线观看| 最近最新免费中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 日韩精品中文字幕看吧| 久久婷婷人人爽人人干人人爱| 最近最新中文字幕大全免费视频| 亚洲第一欧美日韩一区二区三区| 国产乱人伦免费视频| av国产免费在线观看| 我要搜黄色片| 熟女少妇亚洲综合色aaa.| 国产探花极品一区二区| 久久精品夜夜夜夜夜久久蜜豆| 成人午夜高清在线视频| 欧美另类亚洲清纯唯美| 在线观看免费视频日本深夜| 99精品欧美一区二区三区四区| 日本一二三区视频观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本黄色视频三级网站网址| 哪里可以看免费的av片| 国内精品久久久久久久电影| 欧美日韩瑟瑟在线播放| 国产免费av片在线观看野外av| 国产精品久久久久久人妻精品电影| 中出人妻视频一区二区| 波野结衣二区三区在线 | 很黄的视频免费| 国产激情偷乱视频一区二区| 女人高潮潮喷娇喘18禁视频| 欧美激情在线99| 一卡2卡三卡四卡精品乱码亚洲| 好看av亚洲va欧美ⅴa在| 亚洲精品一区av在线观看| 嫩草影院精品99| 99久久精品国产亚洲精品| 亚洲成av人片在线播放无| 欧美一区二区精品小视频在线| 亚洲av成人av| 亚洲精华国产精华精| 精品99又大又爽又粗少妇毛片 | 婷婷丁香在线五月| 国内精品一区二区在线观看| 国产欧美日韩一区二区三| 狂野欧美激情性xxxx| 看黄色毛片网站| 国产三级中文精品| 亚洲av第一区精品v没综合| 网址你懂的国产日韩在线| 成人特级黄色片久久久久久久| 欧美bdsm另类| 最近视频中文字幕2019在线8| 最后的刺客免费高清国语| 天堂影院成人在线观看| ponron亚洲| 啦啦啦观看免费观看视频高清| 亚洲欧美精品综合久久99| 国产精品一区二区三区四区久久| 99精品欧美一区二区三区四区| 中文字幕熟女人妻在线| 十八禁人妻一区二区| 婷婷丁香在线五月| 丁香六月欧美| 三级男女做爰猛烈吃奶摸视频| 又粗又爽又猛毛片免费看| 精品久久久久久,| avwww免费| 免费看a级黄色片| 老熟妇乱子伦视频在线观看| 日韩中文字幕欧美一区二区| 亚洲成人久久爱视频| 国产日本99.免费观看| 欧美色视频一区免费| 在线观看日韩欧美| 老鸭窝网址在线观看| 18禁裸乳无遮挡免费网站照片| 女警被强在线播放| www.熟女人妻精品国产| 美女被艹到高潮喷水动态| 欧美不卡视频在线免费观看| 国产视频一区二区在线看| 在线免费观看的www视频| 欧美乱妇无乱码| 成人一区二区视频在线观看| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩卡通动漫| 久久久精品大字幕| 最好的美女福利视频网| 中国美女看黄片| 怎么达到女性高潮| 亚洲国产精品999在线| 久久久久久久精品吃奶| 国产精品一区二区免费欧美| 在线观看av片永久免费下载| 欧美性感艳星| 久久久久精品国产欧美久久久| 欧美三级亚洲精品| 国产伦人伦偷精品视频| 国产 一区 欧美 日韩| 亚洲中文字幕一区二区三区有码在线看| 99在线人妻在线中文字幕| 少妇人妻精品综合一区二区 | 午夜日韩欧美国产| 国产免费av片在线观看野外av| 成年人黄色毛片网站| 99热精品在线国产| 精华霜和精华液先用哪个| 一区二区三区免费毛片| 脱女人内裤的视频| eeuss影院久久| 国产精品三级大全| 老司机福利观看| 中文亚洲av片在线观看爽| 国产av不卡久久| x7x7x7水蜜桃| 久久国产精品影院| 99国产综合亚洲精品| 内地一区二区视频在线| 久久久久久九九精品二区国产| 亚洲性夜色夜夜综合| 精品日产1卡2卡| 久久精品91蜜桃| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 舔av片在线| 久久精品国产99精品国产亚洲性色| 婷婷亚洲欧美| 成年女人永久免费观看视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩中文字幕欧美一区二区| 国产伦人伦偷精品视频| 最新美女视频免费是黄的| 好男人电影高清在线观看| 波多野结衣高清作品| 在线免费观看不下载黄p国产 | 看黄色毛片网站| 3wmmmm亚洲av在线观看| 亚洲成人久久爱视频| 亚洲欧美日韩东京热| 宅男免费午夜| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 丰满乱子伦码专区| 亚洲第一电影网av| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| 岛国视频午夜一区免费看| 国内精品久久久久久久电影| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av电影不卡..在线观看| 又紧又爽又黄一区二区| 听说在线观看完整版免费高清| 色尼玛亚洲综合影院| 日本免费一区二区三区高清不卡| 国产精品 国内视频| 免费一级毛片在线播放高清视频| 免费av毛片视频| 精品久久久久久久毛片微露脸| 91麻豆av在线| 国产免费一级a男人的天堂| 一夜夜www| 黑人欧美特级aaaaaa片| 久久精品国产自在天天线| 亚洲成av人片免费观看| 婷婷亚洲欧美| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| 女人高潮潮喷娇喘18禁视频| 久久久国产成人免费| 婷婷亚洲欧美| 亚洲狠狠婷婷综合久久图片| 亚洲在线观看片| 无遮挡黄片免费观看| 久久久久九九精品影院| 日韩国内少妇激情av| 亚洲国产欧洲综合997久久,| 午夜福利成人在线免费观看| 久久久久亚洲av毛片大全| 嫩草影院入口| 午夜视频国产福利| 国产精品免费一区二区三区在线| 色综合欧美亚洲国产小说| 国内精品美女久久久久久| a级毛片a级免费在线| 国产精品永久免费网站| 欧美av亚洲av综合av国产av| 在线免费观看不下载黄p国产 | 在线十欧美十亚洲十日本专区| 99精品欧美一区二区三区四区| 久久精品国产99精品国产亚洲性色| 色哟哟哟哟哟哟| 国产高清有码在线观看视频| 日韩欧美在线乱码| 国产av不卡久久| 欧美另类亚洲清纯唯美| 国产又黄又爽又无遮挡在线| 狂野欧美激情性xxxx| 亚洲 欧美 日韩 在线 免费| 免费搜索国产男女视频| 麻豆一二三区av精品| 色综合亚洲欧美另类图片| 很黄的视频免费| 国产精品久久久久久久电影 | 精品不卡国产一区二区三区| 免费看美女性在线毛片视频| 国产精品野战在线观看| 国产精品久久久久久精品电影| 亚洲成a人片在线一区二区| 中文亚洲av片在线观看爽| 男女之事视频高清在线观看| 日韩欧美在线乱码| 久久香蕉国产精品| 欧美国产日韩亚洲一区| 丁香欧美五月| 黄色片一级片一级黄色片| 又紧又爽又黄一区二区| 国产精品久久久久久人妻精品电影| 人人妻人人看人人澡| 精品国产美女av久久久久小说| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| а√天堂www在线а√下载| 国产精品野战在线观看| 亚洲av成人av| 精品久久久久久久久久免费视频| 欧美午夜高清在线| 日本在线视频免费播放| www日本黄色视频网| 亚洲人成电影免费在线| 精品一区二区三区视频在线观看免费| 欧美xxxx黑人xx丫x性爽| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 精品久久久久久久久久免费视频| 久久精品国产清高在天天线| 久久久国产成人精品二区| 国产免费一级a男人的天堂| 精品久久久久久久末码| 欧美激情久久久久久爽电影| 超碰av人人做人人爽久久 | 国产精品久久久久久精品电影| 首页视频小说图片口味搜索| 此物有八面人人有两片| 午夜精品久久久久久毛片777| 一个人看的www免费观看视频| av中文乱码字幕在线| 毛片女人毛片| 黄色日韩在线| 国产精品 国内视频| 熟妇人妻久久中文字幕3abv| 国内毛片毛片毛片毛片毛片| 国产熟女xx| 女生性感内裤真人,穿戴方法视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜免费观看网址| 亚洲av电影不卡..在线观看| 国产午夜精品久久久久久一区二区三区 | www.熟女人妻精品国产| 日本 av在线| 欧美丝袜亚洲另类 | 精品国产亚洲在线| 午夜免费激情av| 白带黄色成豆腐渣| 免费人成视频x8x8入口观看| 亚洲欧美日韩东京热| 亚洲性夜色夜夜综合| 国产欧美日韩精品一区二区| 天堂√8在线中文| 亚洲天堂国产精品一区在线| 亚洲激情在线av| 久久草成人影院| 99国产极品粉嫩在线观看| 一级作爱视频免费观看| 三级国产精品欧美在线观看| 一本一本综合久久| 亚洲专区国产一区二区| 国产精品三级大全| 91麻豆av在线| 日韩国内少妇激情av| 国产精品av视频在线免费观看| 成人午夜高清在线视频| 日本在线视频免费播放| 两个人的视频大全免费| 舔av片在线| 蜜桃久久精品国产亚洲av| 久久精品91蜜桃| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 国产毛片a区久久久久| 丰满的人妻完整版| 亚洲欧美日韩高清专用| 欧美高清成人免费视频www| 一a级毛片在线观看| 久久久久国内视频| av专区在线播放| 亚洲精品日韩av片在线观看 | 国产又黄又爽又无遮挡在线| 国产爱豆传媒在线观看| 一夜夜www| 亚洲欧美日韩高清在线视频| 亚洲成人中文字幕在线播放| 90打野战视频偷拍视频| 高潮久久久久久久久久久不卡| 国产精品,欧美在线| 亚洲国产精品久久男人天堂| 色综合亚洲欧美另类图片| 内地一区二区视频在线| 九九久久精品国产亚洲av麻豆| 91麻豆av在线| 精华霜和精华液先用哪个| 性欧美人与动物交配| 国产三级在线视频| 又爽又黄无遮挡网站| 亚洲乱码一区二区免费版| 最新中文字幕久久久久| 俄罗斯特黄特色一大片| 国产av在哪里看| 国产黄色小视频在线观看| 亚洲天堂国产精品一区在线| 日韩免费av在线播放| а√天堂www在线а√下载| 男女那种视频在线观看| 老汉色∧v一级毛片| 中文字幕熟女人妻在线| 国产av一区在线观看免费| 999久久久精品免费观看国产| 亚洲成人久久性| 亚洲av一区综合| 九色国产91popny在线| 精品人妻偷拍中文字幕| 久久6这里有精品| 男女视频在线观看网站免费| 丁香六月欧美| 成年女人永久免费观看视频| 一级毛片高清免费大全| 99在线人妻在线中文字幕| 精品人妻偷拍中文字幕| 国产精品美女特级片免费视频播放器| 嫩草影视91久久| 欧美激情在线99| 国内精品一区二区在线观看| 女人十人毛片免费观看3o分钟| 久久久久国产精品人妻aⅴ院| 亚洲欧美日韩卡通动漫| 变态另类成人亚洲欧美熟女| 亚洲最大成人手机在线| 国产爱豆传媒在线观看| aaaaa片日本免费| 国产精品影院久久| 亚洲精品456在线播放app | 免费无遮挡裸体视频| 69av精品久久久久久| 欧美不卡视频在线免费观看| 午夜老司机福利剧场| 成人特级黄色片久久久久久久| 毛片女人毛片| 成人av一区二区三区在线看| 久久久久性生活片| 午夜福利18| 欧美中文日本在线观看视频| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 国产三级黄色录像| netflix在线观看网站| 18禁在线播放成人免费| 国产精品乱码一区二三区的特点| 最近最新免费中文字幕在线| 精品一区二区三区视频在线观看免费| 欧美成狂野欧美在线观看| 午夜福利在线观看免费完整高清在 | 亚洲欧美日韩高清专用| 97碰自拍视频| 国产v大片淫在线免费观看| 欧美bdsm另类| 观看免费一级毛片| 国产亚洲欧美在线一区二区| 国产精品亚洲美女久久久| 欧美成人性av电影在线观看| 国产精品98久久久久久宅男小说| 给我免费播放毛片高清在线观看| 欧美丝袜亚洲另类 | 女同久久另类99精品国产91| 亚洲av免费高清在线观看| 国产精品综合久久久久久久免费| x7x7x7水蜜桃| 欧美bdsm另类| 国产高清三级在线| 亚洲无线在线观看| 亚洲最大成人中文| 俺也久久电影网| 精品一区二区三区人妻视频| 在线看三级毛片| 97人妻精品一区二区三区麻豆| 蜜桃亚洲精品一区二区三区| 亚洲美女视频黄频| 天堂网av新在线| 久久久国产精品麻豆| 国产精品亚洲美女久久久| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 日本与韩国留学比较| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 国产激情偷乱视频一区二区| 精品免费久久久久久久清纯| 小说图片视频综合网站| 在线观看免费午夜福利视频| 少妇人妻一区二区三区视频| svipshipincom国产片| 狠狠狠狠99中文字幕| 偷拍熟女少妇极品色| 又黄又爽又免费观看的视频| www.www免费av| 国产高潮美女av| 两人在一起打扑克的视频| 国产真人三级小视频在线观看| 亚洲不卡免费看| 欧美黑人欧美精品刺激| 国产亚洲精品综合一区在线观看| 色播亚洲综合网| 动漫黄色视频在线观看| 免费在线观看成人毛片| 国产不卡一卡二| 精品人妻一区二区三区麻豆 | 中文在线观看免费www的网站| 长腿黑丝高跟| 久久久久免费精品人妻一区二区| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在 | 欧美3d第一页| 亚洲激情在线av| 国产av麻豆久久久久久久| 欧美日韩乱码在线| 亚洲成人免费电影在线观看| 51国产日韩欧美| av国产免费在线观看| 国产又黄又爽又无遮挡在线| 90打野战视频偷拍视频| 亚洲av五月六月丁香网| 亚洲av免费高清在线观看| 午夜福利在线在线| svipshipincom国产片| 久久精品影院6| 亚洲一区二区三区不卡视频| 亚洲精品久久国产高清桃花| 嫩草影视91久久| 成熟少妇高潮喷水视频| 十八禁网站免费在线| 欧美日韩国产亚洲二区| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 亚洲av电影不卡..在线观看| 老汉色av国产亚洲站长工具| 国产探花极品一区二区| 午夜福利视频1000在线观看| 男女午夜视频在线观看| 国产高清videossex| 久久精品综合一区二区三区| 婷婷亚洲欧美| 在线观看日韩欧美| 国产不卡一卡二| 人妻丰满熟妇av一区二区三区| 成人特级av手机在线观看| 成人国产综合亚洲| 别揉我奶头~嗯~啊~动态视频| 99国产极品粉嫩在线观看| 国内精品久久久久久久电影| 2021天堂中文幕一二区在线观| 亚洲人成网站在线播| 久久精品国产清高在天天线| 国产高清三级在线| 亚洲精品在线观看二区| 午夜激情欧美在线| 在线观看av片永久免费下载| 可以在线观看的亚洲视频| 亚洲色图av天堂| 成人午夜高清在线视频| 成人高潮视频无遮挡免费网站| 国产精品爽爽va在线观看网站| 嫩草影院入口| 看片在线看免费视频| 久久精品影院6| www.熟女人妻精品国产| av国产免费在线观看| 久久草成人影院| 啦啦啦免费观看视频1| 人妻久久中文字幕网| www日本黄色视频网| 亚洲av电影不卡..在线观看| avwww免费| 日本在线视频免费播放| 色综合亚洲欧美另类图片| 午夜久久久久精精品| 99热只有精品国产| 人妻丰满熟妇av一区二区三区| 国产中年淑女户外野战色| 国产爱豆传媒在线观看| 午夜激情欧美在线| 免费看a级黄色片| 国产成人欧美在线观看| 亚洲国产精品sss在线观看| 午夜福利成人在线免费观看| 国产主播在线观看一区二区| 成人av在线播放网站| 日韩欧美精品v在线| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 亚洲人成网站在线播放欧美日韩| 精品无人区乱码1区二区| 国产精品久久久久久精品电影| 免费av观看视频| 欧美xxxx黑人xx丫x性爽| 人人妻,人人澡人人爽秒播| 中文字幕人妻熟人妻熟丝袜美 | 精品国内亚洲2022精品成人| 丝袜美腿在线中文| 97碰自拍视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本免费a在线| 999久久久精品免费观看国产| 国产成人aa在线观看| 亚洲av成人不卡在线观看播放网| 国产69精品久久久久777片| 小说图片视频综合网站| 久久久久久久久久黄片| 两个人看的免费小视频| 淫秽高清视频在线观看| 国产探花极品一区二区| 国产蜜桃级精品一区二区三区| 高潮久久久久久久久久久不卡| 国产一区二区亚洲精品在线观看| 午夜日韩欧美国产| 69人妻影院| 又紧又爽又黄一区二区| 国产亚洲精品久久久com| 三级国产精品欧美在线观看| 日本黄色片子视频| 久久亚洲精品不卡| 久久久精品大字幕| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品合色在线| 嫁个100分男人电影在线观看| 一进一出抽搐gif免费好疼| av片东京热男人的天堂| 国产精品香港三级国产av潘金莲| 观看免费一级毛片| 丰满乱子伦码专区| 亚洲最大成人中文| 一本久久中文字幕| 91字幕亚洲| 国内精品久久久久久久电影| 18禁黄网站禁片免费观看直播| 熟女少妇亚洲综合色aaa.| 久久精品国产清高在天天线| 脱女人内裤的视频| 黄色视频,在线免费观看| 婷婷亚洲欧美| 黄片大片在线免费观看| 成人国产一区最新在线观看| 首页视频小说图片口味搜索| 国产免费av片在线观看野外av| 脱女人内裤的视频| 听说在线观看完整版免费高清| 婷婷亚洲欧美| 禁无遮挡网站| 精品人妻1区二区| 日韩欧美在线二视频| 亚洲国产欧洲综合997久久,| 国产亚洲精品综合一区在线观看| 亚洲av一区综合| 欧美最新免费一区二区三区 | 全区人妻精品视频| 国产成人aa在线观看| 国产成人系列免费观看| 天堂影院成人在线观看| 国产一区二区在线av高清观看| 又黄又粗又硬又大视频|