• <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ì)及分析
    老熟女久久久| 欧美激情国产日韩精品一区| 你懂的网址亚洲精品在线观看| 国产精品偷伦视频观看了| 一边亲一边摸免费视频| 中文字幕精品免费在线观看视频 | 亚洲经典国产精华液单| 成人黄色视频免费在线看| 男女边吃奶边做爰视频| av不卡在线播放| a级毛片黄视频| 国内精品宾馆在线| 一本久久精品| 亚洲四区av| 性色av一级| 内地一区二区视频在线| 国产亚洲精品第一综合不卡 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲天堂av无毛| 亚洲美女视频黄频| 夫妻午夜视频| 97在线视频观看| 国产视频首页在线观看| 黑人猛操日本美女一级片| 亚洲欧洲精品一区二区精品久久久 | 久久狼人影院| 日日啪夜夜爽| 国产亚洲av片在线观看秒播厂| 国产高清三级在线| av视频免费观看在线观看| 亚洲国产成人一精品久久久| 我要看黄色一级片免费的| 99久国产av精品国产电影| 国产精品女同一区二区软件| 欧美另类一区| 边亲边吃奶的免费视频| 国产亚洲欧美精品永久| 国产成人aa在线观看| 国产xxxxx性猛交| 中文字幕人妻丝袜制服| 十八禁网站网址无遮挡| 国产一区亚洲一区在线观看| 久久久久久伊人网av| 亚洲成国产人片在线观看| 久久久久人妻精品一区果冻| 两个人看的免费小视频| 亚洲 欧美一区二区三区| 国产免费视频播放在线视频| 国产免费视频播放在线视频| 性色avwww在线观看| 少妇 在线观看| 久久精品国产亚洲av天美| 午夜福利乱码中文字幕| 色5月婷婷丁香| 国产毛片在线视频| 高清视频免费观看一区二区| 18在线观看网站| 欧美另类一区| 久热久热在线精品观看| 国产成人精品在线电影| 亚洲成国产人片在线观看| 最后的刺客免费高清国语| 超碰97精品在线观看| 国产伦理片在线播放av一区| 人人妻人人澡人人看| 免费观看无遮挡的男女| 欧美成人精品欧美一级黄| 精品人妻熟女毛片av久久网站| 国产精品偷伦视频观看了| 久久热在线av| 99热全是精品| 九草在线视频观看| 亚洲av.av天堂| videossex国产| 伊人久久国产一区二区| 最近最新中文字幕大全免费视频 | av在线app专区| 欧美最新免费一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 欧美精品亚洲一区二区| 欧美人与善性xxx| 国产一区二区三区av在线| xxxhd国产人妻xxx| 久久精品国产亚洲av天美| 国产一区二区激情短视频 | 亚洲成人手机| 日本欧美视频一区| 丝袜人妻中文字幕| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 不卡视频在线观看欧美| 国产一区二区在线观看av| 咕卡用的链子| 精品视频人人做人人爽| 久久99蜜桃精品久久| 青青草视频在线视频观看| 妹子高潮喷水视频| 黑丝袜美女国产一区| 国产成人91sexporn| 大片免费播放器 马上看| 乱人伦中国视频| 99国产综合亚洲精品| av免费在线看不卡| 国产 精品1| 另类亚洲欧美激情| 日韩中字成人| 亚洲高清免费不卡视频| 成人综合一区亚洲| 看十八女毛片水多多多| 精品熟女少妇av免费看| 91午夜精品亚洲一区二区三区| 一级,二级,三级黄色视频| 日韩中字成人| 亚洲欧洲日产国产| 观看美女的网站| 春色校园在线视频观看| 99热这里只有是精品在线观看| 亚洲美女搞黄在线观看| 国产色婷婷99| 捣出白浆h1v1| 满18在线观看网站| 国产 一区精品| 18禁裸乳无遮挡动漫免费视频| 欧美性感艳星| 国产老妇伦熟女老妇高清| 少妇 在线观看| 久久久欧美国产精品| 日韩一本色道免费dvd| 啦啦啦在线观看免费高清www| 女性被躁到高潮视频| 亚洲国产看品久久| 汤姆久久久久久久影院中文字幕| 亚洲,一卡二卡三卡| 狠狠精品人妻久久久久久综合| 久久鲁丝午夜福利片| 69精品国产乱码久久久| 黄色毛片三级朝国网站| 天天影视国产精品| 国产黄色免费在线视频| 久久青草综合色| 飞空精品影院首页| 亚洲成人av在线免费| 少妇人妻久久综合中文| 中文乱码字字幕精品一区二区三区| 国产免费又黄又爽又色| 亚洲精品久久久久久婷婷小说| 搡老乐熟女国产| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验| 日本午夜av视频| 日日摸夜夜添夜夜爱| 久久久久久久国产电影| 插逼视频在线观看| 亚洲内射少妇av| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 国产成人免费观看mmmm| 国产精品99久久99久久久不卡 | 99久久综合免费| 久久精品国产a三级三级三级| 插逼视频在线观看| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频| 亚洲精品久久成人aⅴ小说| 久久精品久久久久久久性| 久久久欧美国产精品| 亚洲欧美一区二区三区黑人 | 大香蕉久久网| 亚洲精品乱码久久久久久按摩| 日韩制服骚丝袜av| 久久久久人妻精品一区果冻| 成年人午夜在线观看视频| 精品人妻一区二区三区麻豆| 国产精品.久久久| 香蕉国产在线看| 国产视频首页在线观看| 国产 精品1| 精品一区二区三卡| av.在线天堂| 国产欧美另类精品又又久久亚洲欧美| 又黄又粗又硬又大视频| 精品人妻一区二区三区麻豆| 夫妻性生交免费视频一级片| 国产精品国产三级专区第一集| 中文精品一卡2卡3卡4更新| 国产老妇伦熟女老妇高清| 亚洲,欧美,日韩| 精品一区二区免费观看| 欧美3d第一页| 交换朋友夫妻互换小说| av卡一久久| 日韩视频在线欧美| 日本wwww免费看| 男人爽女人下面视频在线观看| 99国产精品免费福利视频| 色婷婷av一区二区三区视频| 久久久精品区二区三区| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 精品视频人人做人人爽| 国产av精品麻豆| 久久精品夜色国产| 高清在线视频一区二区三区| 亚洲精品一二三| 国产成人精品无人区| 亚洲成国产人片在线观看| 午夜福利影视在线免费观看| 99久久人妻综合| 夫妻性生交免费视频一级片| 精品国产一区二区久久| 久久久久久久精品精品| 少妇精品久久久久久久| 香蕉国产在线看| 大香蕉97超碰在线| 国产探花极品一区二区| 欧美国产精品一级二级三级| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 99久久精品国产国产毛片| 久久亚洲国产成人精品v| 国产 一区精品| 黄色怎么调成土黄色| 老女人水多毛片| 国产精品一国产av| 水蜜桃什么品种好| 亚洲欧美一区二区三区黑人 | a级毛片黄视频| 五月开心婷婷网| 国国产精品蜜臀av免费| 免费少妇av软件| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 好男人视频免费观看在线| 国产精品不卡视频一区二区| 免费久久久久久久精品成人欧美视频 | 免费看不卡的av| 精品一区二区免费观看| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| 国产在视频线精品| 大香蕉久久网| 精品少妇久久久久久888优播| 久久国产精品男人的天堂亚洲 | 国产色婷婷99| 亚洲一级一片aⅴ在线观看| 国产69精品久久久久777片| 国产成人精品久久久久久| 啦啦啦中文免费视频观看日本| 午夜福利乱码中文字幕| 在线观看三级黄色| 精品亚洲成国产av| 国产极品粉嫩免费观看在线| 天堂俺去俺来也www色官网| 亚洲av日韩在线播放| 国产成人a∨麻豆精品| 国产又色又爽无遮挡免| 一级片'在线观看视频| 狠狠精品人妻久久久久久综合| 午夜日本视频在线| 国产精品久久久久久精品电影小说| videosex国产| 啦啦啦视频在线资源免费观看| 日韩欧美一区视频在线观看| av免费在线看不卡| 亚洲精品一区蜜桃| 亚洲美女视频黄频| 久久久久久人人人人人| 成人黄色视频免费在线看| 香蕉精品网在线| 中文字幕人妻丝袜制服| videossex国产| 十八禁网站网址无遮挡| 黄片播放在线免费| 日韩成人伦理影院| 久久ye,这里只有精品| av在线老鸭窝| 亚洲欧洲国产日韩| 欧美日韩亚洲高清精品| 国产 一区精品| 少妇人妻精品综合一区二区| 交换朋友夫妻互换小说| 五月天丁香电影| 成人亚洲欧美一区二区av| 青青草视频在线视频观看| 久久精品aⅴ一区二区三区四区 | 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 一级毛片电影观看| 国产av一区二区精品久久| 色吧在线观看| 亚洲精品一区蜜桃| 精品久久久久久电影网| 天天躁夜夜躁狠狠躁躁| 国产极品粉嫩免费观看在线| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 在线免费观看不下载黄p国产| 日本黄大片高清| 王馨瑶露胸无遮挡在线观看| 久久久国产精品麻豆| 综合色丁香网| 青春草视频在线免费观看| 欧美xxⅹ黑人| 欧美激情国产日韩精品一区| 寂寞人妻少妇视频99o| 婷婷色综合大香蕉| 午夜福利视频精品| 黄色配什么色好看| 日本爱情动作片www.在线观看| 乱人伦中国视频| 日韩在线高清观看一区二区三区| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 成人国产麻豆网| 热re99久久精品国产66热6| a级毛片在线看网站| 老司机亚洲免费影院| 亚洲人成77777在线视频| 国精品久久久久久国模美| 欧美激情 高清一区二区三区| 中文精品一卡2卡3卡4更新| 97在线人人人人妻| 亚洲av.av天堂| 男人舔女人的私密视频| 看非洲黑人一级黄片| 另类精品久久| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 男女边摸边吃奶| 成年人午夜在线观看视频| 肉色欧美久久久久久久蜜桃| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 性色avwww在线观看| 色94色欧美一区二区| 国精品久久久久久国模美| 久久久久久久久久成人| 在线观看美女被高潮喷水网站| 18+在线观看网站| 少妇的逼好多水| 夜夜爽夜夜爽视频| 有码 亚洲区| 免费少妇av软件| 少妇的丰满在线观看| av黄色大香蕉| 女人精品久久久久毛片| 又黄又爽又刺激的免费视频.| av片东京热男人的天堂| 国产精品国产三级国产专区5o| 免费av不卡在线播放| 亚洲一级一片aⅴ在线观看| 人妻系列 视频| 成年人免费黄色播放视频| 十分钟在线观看高清视频www| 97精品久久久久久久久久精品| 成年人免费黄色播放视频| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 精品国产露脸久久av麻豆| 欧美精品av麻豆av| 午夜精品国产一区二区电影| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 免费少妇av软件| 亚洲熟女精品中文字幕| 亚洲欧美清纯卡通| 国产成人午夜福利电影在线观看| 日韩制服骚丝袜av| 在线免费观看不下载黄p国产| 国产白丝娇喘喷水9色精品| 国产男女内射视频| 丝袜在线中文字幕| 日本与韩国留学比较| 久久久精品区二区三区| 日本色播在线视频| 久久久久久久国产电影| 色婷婷av一区二区三区视频| 又黄又爽又刺激的免费视频.| 免费不卡的大黄色大毛片视频在线观看| 老女人水多毛片| 免费不卡的大黄色大毛片视频在线观看| 老女人水多毛片| 一本色道久久久久久精品综合| 日本91视频免费播放| 美女国产视频在线观看| tube8黄色片| 各种免费的搞黄视频| 国产亚洲精品第一综合不卡 | 成年av动漫网址| 免费在线观看黄色视频的| 在线亚洲精品国产二区图片欧美| 日韩制服骚丝袜av| 一级爰片在线观看| 免费av中文字幕在线| 国产成人午夜福利电影在线观看| 天堂8中文在线网| 久久久久久久久久成人| 卡戴珊不雅视频在线播放| 亚洲av电影在线进入| 一级毛片我不卡| 免费观看在线日韩| 水蜜桃什么品种好| 91国产中文字幕| 免费看av在线观看网站| 26uuu在线亚洲综合色| 国产精品蜜桃在线观看| 日韩av不卡免费在线播放| 久久韩国三级中文字幕| 午夜免费男女啪啪视频观看| 亚洲国产精品专区欧美| 一级黄片播放器| 国产男女超爽视频在线观看| 国产视频首页在线观看| 在线观看人妻少妇| 亚洲精品久久久久久婷婷小说| 美国免费a级毛片| 欧美+日韩+精品| 日韩中文字幕视频在线看片| 亚洲精品456在线播放app| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩卡通动漫| 少妇被粗大的猛进出69影院 | 2021少妇久久久久久久久久久| 深夜精品福利| 亚洲精品色激情综合| 欧美成人精品欧美一级黄| 午夜日本视频在线| 一二三四在线观看免费中文在 | 老司机影院成人| 成人毛片a级毛片在线播放| 最近2019中文字幕mv第一页| 国产黄频视频在线观看| 欧美日韩亚洲高清精品| 我的女老师完整版在线观看| 秋霞伦理黄片| 天美传媒精品一区二区| 99精国产麻豆久久婷婷| 国产老妇伦熟女老妇高清| 男女边摸边吃奶| 另类精品久久| 国产成人91sexporn| 国产免费又黄又爽又色| 青春草亚洲视频在线观看| 日韩,欧美,国产一区二区三区| 午夜福利视频精品| 在线看a的网站| 国产黄色视频一区二区在线观看| 亚洲天堂av无毛| 久久精品国产a三级三级三级| 黑人高潮一二区| h视频一区二区三区| 亚洲欧洲日产国产| 十八禁网站网址无遮挡| 美国免费a级毛片| 男女免费视频国产| 最近手机中文字幕大全| 九草在线视频观看| av不卡在线播放| 亚洲av中文av极速乱| 香蕉精品网在线| 女人精品久久久久毛片| 赤兔流量卡办理| 色视频在线一区二区三区| 成人亚洲精品一区在线观看| 亚洲国产精品成人久久小说| 欧美日韩视频精品一区| av黄色大香蕉| 精品一区在线观看国产| 伦理电影大哥的女人| 国产亚洲精品久久久com| 国产精品 国内视频| 亚洲精品456在线播放app| 日本欧美国产在线视频| 免费观看在线日韩| 免费黄色在线免费观看| 69精品国产乱码久久久| 亚洲欧美色中文字幕在线| 国产视频首页在线观看| 日本黄色日本黄色录像| 免费大片18禁| 亚洲,欧美精品.| 777米奇影视久久| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 日韩成人伦理影院| 亚洲国产最新在线播放| 99热这里只有是精品在线观看| 看非洲黑人一级黄片| 国产精品麻豆人妻色哟哟久久| 久久ye,这里只有精品| 女性生殖器流出的白浆| 国产精品久久久久久精品电影小说| 免费高清在线观看视频在线观看| 老熟女久久久| 97在线视频观看| 国产成人精品福利久久| 中文乱码字字幕精品一区二区三区| av播播在线观看一区| 精品人妻一区二区三区麻豆| 亚洲精品国产色婷婷电影| 成人免费观看视频高清| 在线 av 中文字幕| 国产精品无大码| 久久久国产欧美日韩av| 欧美精品av麻豆av| 国产成人精品福利久久| 曰老女人黄片| 婷婷色综合大香蕉| 亚洲国产看品久久| 久久影院123| 免费大片黄手机在线观看| 在线观看免费高清a一片| 在线天堂中文资源库| 美女xxoo啪啪120秒动态图| 高清视频免费观看一区二区| 久久热在线av| 免费在线观看完整版高清| 久久综合国产亚洲精品| 午夜福利,免费看| 校园人妻丝袜中文字幕| 久久久久国产网址| 国产1区2区3区精品| 亚洲国产日韩一区二区| 卡戴珊不雅视频在线播放| 久久国产精品大桥未久av| 国语对白做爰xxxⅹ性视频网站| 欧美日韩综合久久久久久| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 五月天丁香电影| 精品人妻偷拍中文字幕| 97人妻天天添夜夜摸| freevideosex欧美| 另类精品久久| 99香蕉大伊视频| 九九爱精品视频在线观看| 亚洲国产av新网站| 丝瓜视频免费看黄片| 精品久久蜜臀av无| 日本黄色日本黄色录像| 曰老女人黄片| 人人妻人人澡人人看| 成人午夜精彩视频在线观看| 18禁在线无遮挡免费观看视频| 久久免费观看电影| 伊人久久国产一区二区| 两个人看的免费小视频| 亚洲欧美精品自产自拍| 国产成人欧美| 国产男女内射视频| 90打野战视频偷拍视频| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 亚洲国产精品999| 国产亚洲精品久久久com| 亚洲熟女精品中文字幕| 免费人成在线观看视频色| 久久韩国三级中文字幕| 这个男人来自地球电影免费观看 | 女的被弄到高潮叫床怎么办| 亚洲av福利一区| 亚洲中文av在线| 侵犯人妻中文字幕一二三四区| 成人国产麻豆网| 国产精品人妻久久久久久| 国产在视频线精品| 免费观看a级毛片全部| 久久狼人影院| 大片免费播放器 马上看| 日产精品乱码卡一卡2卡三| 男女午夜视频在线观看 | 99视频精品全部免费 在线| 天天躁夜夜躁狠狠躁躁| 一级毛片 在线播放| 亚洲久久久国产精品| 男人爽女人下面视频在线观看| 人人妻人人澡人人看| 亚洲精品一二三| 2021少妇久久久久久久久久久| 欧美日韩精品成人综合77777| 久久久精品免费免费高清| 日韩欧美精品免费久久| 男人操女人黄网站| 久热这里只有精品99| 最后的刺客免费高清国语| 欧美精品一区二区免费开放| 久久精品熟女亚洲av麻豆精品| av卡一久久| 大香蕉久久成人网| 欧美97在线视频| 美女内射精品一级片tv| av播播在线观看一区| 国产淫语在线视频| av国产久精品久网站免费入址| 亚洲美女黄色视频免费看| 久久精品熟女亚洲av麻豆精品| 好男人视频免费观看在线| 如日韩欧美国产精品一区二区三区| 国产xxxxx性猛交| 嫩草影院入口| 国产一区二区三区av在线| 一区二区日韩欧美中文字幕 | 国产一区二区激情短视频 | 亚洲成av片中文字幕在线观看 | 午夜精品国产一区二区电影| 中文字幕最新亚洲高清| 国产免费福利视频在线观看| 日韩制服丝袜自拍偷拍| 精品第一国产精品| 交换朋友夫妻互换小说| 涩涩av久久男人的天堂| 亚洲欧美日韩卡通动漫| 欧美日韩成人在线一区二区| 亚洲第一av免费看|