999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

燃?xì)廨啓C(jī)排氣擴(kuò)壓器研究現(xiàn)狀

2021-07-26 09:39:16邱彬付經(jīng)倫
發(fā)電技術(shù) 2021年4期
關(guān)鍵詞:優(yōu)化影響

邱彬,付經(jīng)倫

燃?xì)廨啓C(jī)排氣擴(kuò)壓器研究現(xiàn)狀

邱彬1,3,4,5,付經(jīng)倫1,2,3,4,5*

(1.中國(guó)科學(xué)院工程熱物理研究所先進(jìn)燃?xì)廨啓C(jī)實(shí)驗(yàn)室,北京市 海淀區(qū) 100190;2.中國(guó)科學(xué)院工程熱物理研究所南京未來能源系統(tǒng)研究院,江蘇省 南京市 211135;3.中國(guó)科學(xué)院大學(xué),北京市 海淀區(qū) 100049;4.中國(guó)科學(xué)院先進(jìn)能源動(dòng)力重點(diǎn)實(shí)驗(yàn)室,北京市 海淀區(qū) 100190;5.中國(guó)科學(xué)院輕型動(dòng)力創(chuàng)新研究院,北京市 海淀區(qū) 100190)

燃?xì)廨啓C(jī)是關(guān)系國(guó)家安全和國(guó)民經(jīng)濟(jì)發(fā)展的核心裝備,排氣擴(kuò)壓器氣動(dòng)性能對(duì)燃?xì)廨啓C(jī)整機(jī)效率有重要影響。排氣擴(kuò)壓器與上游流場(chǎng)的相互作用及其復(fù)雜的幾何結(jié)構(gòu),導(dǎo)致了擴(kuò)壓器內(nèi)流場(chǎng)易分離、多損失的流動(dòng)特征。為探明擴(kuò)壓器的流動(dòng)特征及其與透平的耦合機(jī)理,綜述了國(guó)內(nèi)外擴(kuò)壓器氣動(dòng)性能和優(yōu)化設(shè)計(jì)的研究現(xiàn)狀,討論了擴(kuò)壓器設(shè)計(jì)與優(yōu)化方面未來的研究方向,未來的研究需要充分考慮透平和擴(kuò)壓器的相互作用以及將參數(shù)化等方法引入擴(kuò)壓器優(yōu)化設(shè)計(jì)問題以提高工作效率。

燃?xì)廨啓C(jī);排氣擴(kuò)壓器;性能研究

0 引言

重型燃?xì)廨啓C(jī)是一種功率密度大、啟動(dòng)速度快、噪音低頻分量低、清潔先進(jìn)的動(dòng)力裝置[1],是關(guān)系國(guó)家安全和國(guó)民經(jīng)濟(jì)發(fā)展的核心裝備,是一個(gè)國(guó)家科技和工業(yè)整體實(shí)力的重要標(biāo)志之一。

高溫氣流在透平中膨脹做功后進(jìn)入排氣擴(kuò)壓器,排氣擴(kuò)壓器用于引導(dǎo)氣流同時(shí)回收末級(jí)透平的余速動(dòng)能。現(xiàn)代燃?xì)廨啓C(jī)排氣余速動(dòng)能占整個(gè)透平功的10%[2]。針對(duì)某H級(jí)燃?xì)廨啓C(jī),經(jīng)熱力計(jì)算,透平排氣損失降低1000Pa使燃機(jī)整機(jī)循環(huán)效率增加0.14,輸出功率增加1.4MW。可見,排氣擴(kuò)壓器性能對(duì)整機(jī)效率和輸出功率的影響不容小覷。

為了提高排氣擴(kuò)壓器氣動(dòng)性能,需要清楚認(rèn)識(shí)其內(nèi)部的流動(dòng)機(jī)理,明確擴(kuò)壓器研究和設(shè)計(jì)的發(fā)展方向。本文通過對(duì)燃?xì)廨啓C(jī)排氣擴(kuò)壓器發(fā)展歷程、研究現(xiàn)狀進(jìn)行梳理和分析,為擴(kuò)壓器的進(jìn)一步研究提供參考依據(jù)。

1 燃?xì)廨啓C(jī)擴(kuò)壓器發(fā)展歷程

重型燃?xì)鈹U(kuò)壓器有徑向和軸向2種形式[3],其主要排氣結(jié)構(gòu)和參數(shù)如表1所示。三大燃?xì)廨啓C(jī)自F級(jí)燃?xì)廨啓C(jī)逐漸改為冷端輸出軸向排氣,排氣擴(kuò)壓器直接與余熱鍋爐進(jìn)氣導(dǎo)管相聯(lián),排氣溫度和排氣量逐漸提高,軸向尺寸變化不大。

表1 世界主要重型燃?xì)廨啓C(jī)排氣結(jié)構(gòu)及排氣參數(shù)

以三菱F級(jí)燃?xì)廨啓C(jī)為例,如圖1(a)[4]所示,排氣擴(kuò)壓通道是由內(nèi)、外錐體形成的一個(gè)流通面積逐漸擴(kuò)大的流道,通過降低流速,形成盡量低的背壓和排氣損失。其中,內(nèi)錐體用于保護(hù)軸承箱,外錐體用于防止排氣缸體過熱。內(nèi)外錐體間通過支撐件支撐,支撐件可為管道系統(tǒng)和儀器儀表提供通道,支撐件外有支撐套,冷卻空氣通過支撐套冷卻支撐件。三菱F級(jí)燃?xì)廨啓C(jī)的排氣通道沿軸向分為2段,前排氣通道內(nèi)布置6個(gè)周向等距的支撐件,后排氣通道布置2個(gè)中心支撐件,如圖1(b)所示。前排氣通道段內(nèi)錐體半徑不變,外錐體半徑沿軸向增大;后排氣通道段內(nèi)錐體半徑沿軸向減小,外錐體半徑不變,前后段外錐體由膨脹節(jié)連接以抵消前后段熱膨脹差異。J級(jí)燃?xì)廨啓C(jī)前后段排氣通道內(nèi)采用相似的支撐結(jié)構(gòu),排氣通道內(nèi)錐體直徑不變,外錐體為錐形,如 圖1(c)[5]所示。

圖1 M701J型燃?xì)廨啓C(jī)示意圖

GE公司首臺(tái)E級(jí)燃機(jī)采用熱端輸出徑向排氣的方式,透平后接環(huán)形擴(kuò)壓器和轉(zhuǎn)折90°的蝸殼。環(huán)形擴(kuò)壓器內(nèi)徑不變,外環(huán)型線為具有擴(kuò)張角的直線,擴(kuò)壓器內(nèi)外錐體間有支撐結(jié)構(gòu),出口布置多個(gè)導(dǎo)流板[6]。F級(jí)燃機(jī)改為冷端輸出軸向排氣形式,內(nèi)錐體為等徑圓柱,外錐體為半徑逐漸增大的錐形。

西門子燃機(jī)擴(kuò)壓器的設(shè)計(jì)一脈相承,從E級(jí)到J級(jí)都是采用環(huán)形擴(kuò)壓器與錐形擴(kuò)壓器組合的軸向排氣方式,如圖2所示。環(huán)形擴(kuò)壓器內(nèi)有周向排布的支撐件,2000E排氣支撐件6個(gè), 4000F支撐件減少為5個(gè),支撐件軸向弦長(zhǎng)增大。環(huán)形擴(kuò)壓段后連接錐形擴(kuò)壓段,錐形擴(kuò)壓段無內(nèi)錐體。

圖2 西門子各級(jí)燃?xì)廨啓C(jī)

2 擴(kuò)壓器氣動(dòng)性能

2.1 擴(kuò)壓器氣動(dòng)性能衡量

擴(kuò)壓器氣動(dòng)性能通常采用總壓損失系數(shù)tpl和靜壓恢復(fù)系數(shù)pr[7]來衡量,tpl和pr定義如下:

tpl越大表示擴(kuò)壓器內(nèi)的流動(dòng)損失越大,pr越大表示擴(kuò)壓器擴(kuò)壓能力越好,擴(kuò)壓器幾何結(jié)構(gòu)和上游流場(chǎng)是影響擴(kuò)壓器氣動(dòng)性能的主要因素。

2.2 幾何結(jié)構(gòu)對(duì)擴(kuò)壓器性能的影響

擴(kuò)壓器內(nèi)外錐體的擴(kuò)張角、軸向長(zhǎng)度、支撐套的外部輪廓結(jié)構(gòu)等[8-9]幾何參數(shù)決定著擴(kuò)壓器的通流結(jié)構(gòu)。許多學(xué)者進(jìn)行了大量的實(shí)驗(yàn)和數(shù)值研究,對(duì)比分析這些幾何參數(shù)對(duì)擴(kuò)壓通道氣動(dòng)性能的影響。

2.2.1 不考慮支撐時(shí)的影響分析

不考慮支撐時(shí),排氣擴(kuò)壓器的基本形式類似擴(kuò)張噴管,這種結(jié)構(gòu)形式廣泛應(yīng)用于各種情況,被深入研究。相同進(jìn)口面積和擴(kuò)張角時(shí),隨著軸向長(zhǎng)度增加,擴(kuò)壓器出口面積增大,靜壓恢復(fù)增加[10]。王廣[11]采用擴(kuò)散度把擴(kuò)張角、軸向長(zhǎng)度2個(gè)參數(shù)統(tǒng)一考慮,擴(kuò)散度的定義為

式中:¢為擴(kuò)壓器軸向長(zhǎng)度0與擴(kuò)壓器入口環(huán)面外徑0的比值;為擴(kuò)壓器外錐面2條與中軸線共面的母線的夾角。在相同擴(kuò)散度時(shí),當(dāng)擴(kuò)壓器內(nèi)未出現(xiàn)流動(dòng)分離時(shí),靜壓恢復(fù)系數(shù)相同;當(dāng)擴(kuò)壓器內(nèi)出現(xiàn)流動(dòng)分離后,軸向長(zhǎng)度越大,擴(kuò)張角越小,流動(dòng)分離程度越小,靜壓恢復(fù)系數(shù)越大。

2.2.2 考慮支撐時(shí)的影響分析

擴(kuò)壓器內(nèi)支撐件的存在明顯增加擴(kuò)壓通道內(nèi)流動(dòng)的復(fù)雜性。Stefano等[9,12]通過實(shí)驗(yàn)對(duì)比有、無支撐的流動(dòng),發(fā)現(xiàn)無支撐擴(kuò)壓器靜壓恢復(fù)接近理想狀態(tài),加入支撐后靜壓恢復(fù)有所下降,這是因?yàn)橹螘?huì)產(chǎn)生尾跡,導(dǎo)致端壁附近發(fā)生流動(dòng)分離。Prakash等[13]得到相同的結(jié)論,但同時(shí)發(fā)現(xiàn)氣流壓力在支撐前緣附近降低,在后緣及近下游位置處則快速提高,如圖3所示。支撐套的輪廓形狀僅對(duì)支撐前緣至1.5倍弦長(zhǎng)下游范圍內(nèi)的壓力產(chǎn)生不同的影響,其中錐形較基線型的影響較小。Fric等[14]發(fā)現(xiàn)錐形支撐套輪廓對(duì)尾跡渦脫落和聲波響應(yīng)有解耦作用,能有效地抑止支撐尾跡渦脫落。

圖3 帶有不同支撐結(jié)構(gòu)擴(kuò)壓器靜壓恢復(fù)

支撐與軸向、徑向的傾角以及支撐的軸向位置都會(huì)影響擴(kuò)壓器內(nèi)的流動(dòng)。董雨軒等[15]對(duì)支撐的徑向傾斜角度和輪廓型線進(jìn)行了數(shù)值研究,發(fā)現(xiàn)支撐型線形狀相同時(shí),帶有徑向傾斜支撐的擴(kuò)壓器比徑向垂直支撐的擴(kuò)壓器通流面積變化更加平緩,總壓損失系數(shù)下降7%~20%。Pradeep等[16]對(duì)支撐與擴(kuò)壓器外錐體相對(duì)位置進(jìn)行研究,該擴(kuò)壓器由環(huán)形擴(kuò)壓段和錐形擴(kuò)壓段組成,環(huán)形擴(kuò)壓段的外環(huán)端壁型線為2段多項(xiàng)式曲線,如圖4[16]所示,圖中為軸向位置,為進(jìn)口直徑。研究結(jié)果見圖5,圖中表示總壓損失系數(shù),可以看出,2段曲線交點(diǎn)位于支撐最大厚度處,靜壓恢復(fù)提高22%,總壓損失降低了約43%;同時(shí)減小擴(kuò)壓器初始擴(kuò)張角,可延緩流動(dòng)分離,改善擴(kuò)壓器性能和出口流動(dòng)的均勻性。

圖4 原始擴(kuò)壓器和2種改進(jìn)環(huán)形擴(kuò)壓器示意圖

圖5 3種擴(kuò)壓器的靜壓恢復(fù)

2.3 上游來流條件對(duì)擴(kuò)壓器性能的影響

機(jī)組實(shí)際運(yùn)行中,透平和擴(kuò)壓器流場(chǎng)之間存在一定的相互作用。隨著燃機(jī)發(fā)展,流量更大,排氣溫度更高,結(jié)構(gòu)相對(duì)更加緊湊,各組件之間氣動(dòng)影響愈加強(qiáng)烈,排氣段內(nèi)仍可檢測(cè)到上游葉片尾跡流場(chǎng)[17]。上游透平葉頂泄漏流、透平出口的非均勻壓力分布、湍流度、旋流角等因素會(huì)影響擴(kuò)壓器內(nèi)的流場(chǎng)和氣動(dòng)性能[18-19]。

2.3.1 葉頂泄漏流對(duì)擴(kuò)壓器內(nèi)的流場(chǎng)和氣動(dòng)性能的影響

葉頂泄漏流強(qiáng)度與邊界層穩(wěn)定之間有直接關(guān)聯(lián)[20],但末級(jí)透平葉頂泄漏流對(duì)擴(kuò)壓器的影響目前尚有爭(zhēng)論。Willinger等[21]研究了末級(jí)透平葉頂間隙與排氣擴(kuò)壓器之間的氣動(dòng)相互作用,指出葉頂間隙從0增大到1.5%弦長(zhǎng)時(shí),葉尖泄漏損失增大,透平輸出功率減小;當(dāng)葉頂間隙大于2.5%弦長(zhǎng)時(shí),葉尖泄漏損失不再增長(zhǎng)。擴(kuò)壓器進(jìn)口機(jī)匣附近環(huán)形壁面射流強(qiáng)度隨間隙的增大而增大。當(dāng)葉頂間隙增大時(shí),擴(kuò)壓器壓力恢復(fù)只有微小增加,對(duì)葉頂間隙依賴性較弱。吳飛等[22]通過數(shù)值方法研究發(fā)現(xiàn),隨著透平動(dòng)葉頂部間隙增大,動(dòng)葉頂部氣流角增大,對(duì)擴(kuò)壓器性能不利;但更多高能氣體進(jìn)入擴(kuò)壓器,有利于提升擴(kuò)壓器性能。當(dāng)葉頂間隙相對(duì)于葉高值從0%增大至0.5%時(shí),高能氣流帶來有利影響占主導(dǎo),擴(kuò)壓器性能逐漸提升;葉頂間隙相對(duì)值在0.5%~1.5%范圍內(nèi)時(shí),氣流角引起的不利影響占主導(dǎo),擴(kuò)壓器性能隨葉頂間隙增大逐漸下降;最優(yōu)間隙值約為0.6%。Zimmermann等[23]和Farokhi[24]指出,增大葉頂間隙產(chǎn)生的葉頂高能射流對(duì)擴(kuò)壓器性能的提升并不足以補(bǔ)償透平中的損失。Mihailowitsch等[25]通過研究發(fā)現(xiàn),在部分工況下,動(dòng)葉出口旋流使支板處發(fā)生分離,增大葉頂間隙會(huì)增大分離;在設(shè)計(jì)和過載工況下,透平增加的損失被擴(kuò)壓器減少的損失補(bǔ)償,所有間隙結(jié)構(gòu)下靜葉進(jìn)口到擴(kuò)壓器出口的系統(tǒng)效率都得以提高;設(shè)計(jì)工況時(shí)小間隙結(jié)構(gòu)最大限度提高效率,過載時(shí)中等間隙結(jié)構(gòu)較理想。Babu等[26]發(fā)現(xiàn)當(dāng)擴(kuò)壓器在環(huán)形擴(kuò)壓段內(nèi)出現(xiàn)分離(因擴(kuò)張角較大)時(shí),葉尖泄漏流動(dòng)的增益才明顯。

Farokhi[25]給出透平功率、擴(kuò)壓器壓力恢復(fù)系數(shù)和葉頂泄漏流之間的無量綱數(shù)學(xué)關(guān)系為

2.3.2 透平出口的非均勻壓力分布的影響

氣流從末級(jí)透平流出時(shí),氣動(dòng)參數(shù)沿徑向不均勻。文獻(xiàn)[27]通過實(shí)驗(yàn)和數(shù)值計(jì)算研究了擴(kuò)壓器入口壓力分布的影響,發(fā)現(xiàn)非均勻壓力分布使擴(kuò)壓器獲得更好的擴(kuò)壓性能。Hirschmann等[28]通過研究發(fā)現(xiàn),葉尖處壓力較大將導(dǎo)致機(jī)匣處發(fā)生流動(dòng)分離,葉根處壓力較大會(huì)引起輪轂處發(fā)生流動(dòng)分離;進(jìn)口總壓徑向分布3%的微小變化將使擴(kuò)壓器內(nèi)具有完全不同的流動(dòng)結(jié)構(gòu)。David等[29]通過非定常計(jì)算發(fā)現(xiàn),葉片尾流和葉柵通道內(nèi)二次流共同作用下環(huán)形擴(kuò)壓器中流體會(huì)重新附著。David 等[30]認(rèn)為周向非均勻的進(jìn)氣條件對(duì)擴(kuò)壓器性能有很大影響,指出動(dòng)葉出口和擴(kuò)壓器進(jìn)口之間的最佳軸向距離為動(dòng)葉弦長(zhǎng)的0.431倍。

2.3.3 湍流度的影響

通常認(rèn)為透平出口湍流度提高有利于改善擴(kuò)壓器性能。Opilat等[31]通過模型實(shí)驗(yàn)發(fā)現(xiàn)來流湍流度提高可以增大擴(kuò)壓器靜壓恢復(fù)能力。Stevens等[32]對(duì)自然發(fā)展?fàn)顟B(tài)和高湍流強(qiáng)度2種進(jìn)口狀態(tài)進(jìn)行對(duì)比,指出擴(kuò)壓器入口湍流強(qiáng)度增大對(duì)流動(dòng)穩(wěn)定性和壓力恢復(fù)的增益可達(dá)20%,而總壓損失僅略有增加。Desideri等[33]和Vassiliev等[34]研究證實(shí)了透平出口湍流有利于抑止流動(dòng)分離發(fā)生,減少流動(dòng)分離區(qū)面積,從而獲得更高的壓力恢復(fù),縮短擴(kuò)壓器長(zhǎng)度。

2.3.4 旋流角的影響

上游透平將導(dǎo)致擴(kuò)壓器入口流動(dòng)方向發(fā)生改變,與流動(dòng)速度相比,速度方向?qū)ο掠螖U(kuò)壓的影響更為明顯。Vassiliev等[34]的研究顯示,擴(kuò)壓器進(jìn)口旋流角對(duì)擴(kuò)壓器壓力恢復(fù)能力的影響程度明顯強(qiáng)于進(jìn)口馬赫數(shù)的影響,進(jìn)口旋流一般有利于提高擴(kuò)壓器性能。Kumar等[17]指出進(jìn)口旋流可提高擴(kuò)壓器靜壓恢復(fù),擴(kuò)壓器軸向長(zhǎng)度越短,旋流的影響越明顯。Song等[35]的無支板擴(kuò)壓器實(shí)驗(yàn)表明,壓力恢復(fù)隨旋流增加而增加,并在旋流角是14°時(shí)達(dá)到峰值,之后迅速下降。Opilat等[31]研究認(rèn)為5°~10°的旋流角對(duì)壓力恢復(fù)有積極的影響,在-5°旋流角時(shí),壓力恢復(fù)值最低。徐倩楠等[36]基于改進(jìn)的流線曲率法發(fā)展了一種環(huán)形擴(kuò)壓器性能預(yù)測(cè)方法,發(fā)現(xiàn)旋流角在0°~17°時(shí),壓力恢復(fù)隨旋流角增大而增大;在17°~25°時(shí),性能下降。徐自榮等[37]對(duì)負(fù)旋流角進(jìn)行研究,發(fā)現(xiàn)旋流角從0°變化至-32°,總壓損失下降4%,在-20°時(shí)下降迅速;靜壓恢復(fù)系數(shù)隨著旋流角增大呈現(xiàn)先上升后降低的變化趨勢(shì),旋流角在-16°時(shí),壓力恢復(fù)系數(shù)達(dá)到最大。

旋流角的影響與擴(kuò)壓器結(jié)構(gòu)相關(guān),對(duì)于不同結(jié)構(gòu)的擴(kuò)壓器,旋流角對(duì)擴(kuò)壓器性能影響不同。Feldcamp等[38]通過實(shí)驗(yàn)發(fā)現(xiàn)在旋流的影響下,支撐形狀與性能有較強(qiáng)的相關(guān)性,支撐數(shù)量?jī)H為次要相關(guān)。無旋流時(shí),支撐對(duì)總壓損失影響較小,增大旋流支撐與性能的相關(guān)性增大。Olaf等[39]在不同旋流角度下,對(duì)不同擴(kuò)張角和支撐的擴(kuò)壓器進(jìn)行研究,發(fā)現(xiàn)擴(kuò)張角為15°的環(huán)形擴(kuò)壓器內(nèi)不存在流動(dòng)分離,旋流角對(duì)擴(kuò)壓性能影響不大;但旋流角對(duì)擴(kuò)張角為20°的環(huán)形擴(kuò)壓器影響明顯,擴(kuò)壓器的靜壓恢復(fù)隨旋流角增大而減小。

擴(kuò)壓器內(nèi)支撐的存在會(huì)改變旋流方向,進(jìn)而改變進(jìn)口旋流對(duì)擴(kuò)壓器性能的影響規(guī)律,同時(shí)不同來流條件下最優(yōu)的支撐型面結(jié)構(gòu)也不同。例如,Stefano等[9,12]通過研究發(fā)現(xiàn),支撐產(chǎn)生的尾跡和末級(jí)透平出口旋流之間的相互作用導(dǎo)致邊界層在輪轂處加厚,支撐下游誘發(fā)流動(dòng)分離,使擴(kuò)壓器性能下降。Flelge等[40]指出旋流角高于10°時(shí),圓柱型支撐比對(duì)稱異型支撐的擴(kuò)壓器性能更好。Pierrasch等[41]同樣發(fā)現(xiàn)在±45°到±22°范圍的旋流角下,無支撐的擴(kuò)壓器比有支撐的擴(kuò)壓器具有更高的壓力恢復(fù)能力;旋流角大于18°時(shí),旋流角對(duì)無支撐擴(kuò)壓器性能無影響;大于23°時(shí),圓柱和異形支撐擴(kuò)壓器壓力恢復(fù)均下降,隨著旋流角的增加,異形支撐擴(kuò)壓器的損失比圓柱形支撐擴(kuò)壓器增加得快。研究表明:出口旋流數(shù)與支撐剖面形狀有較強(qiáng)的相關(guān)性;中等程度的旋流對(duì)支撐的影響是相對(duì)良性的;更大的旋流會(huì)導(dǎo)致流體與支撐分離,并增加總壓損失。

2.3.5 不同工況的影響

燃機(jī)運(yùn)行中,透平與排氣擴(kuò)壓器內(nèi)流場(chǎng)會(huì)隨著工況變化而改變。Sultanian等[42]在3種負(fù)載工況下進(jìn)行了實(shí)驗(yàn)和數(shù)值模擬研究,結(jié)果顯示:全速中載和全速滿載工況下擴(kuò)壓器進(jìn)口到支撐末端段的損失比全速空載下?lián)p失明顯減小;支撐末端到擴(kuò)壓器出口段的損失較全速空載時(shí)略大;在全速滿載工況時(shí),整個(gè)擴(kuò)壓器內(nèi)的損失最小,靜壓恢復(fù)最高。Mimic等[43]首次建立擴(kuò)壓器壓力恢復(fù)與末級(jí)轉(zhuǎn)子整體參數(shù),即載荷系數(shù)、流量系數(shù)和折合轉(zhuǎn)速之間的關(guān)系,發(fā)現(xiàn)載荷系數(shù)、折合轉(zhuǎn)速的增加,流量系數(shù)的減少有利于提高擴(kuò)壓器的壓力恢復(fù)能力。

3 擴(kuò)壓器改進(jìn)和優(yōu)化設(shè)計(jì)

為了提高性能,以往研究人員對(duì)擴(kuò)壓器子午通道參數(shù)和支撐外輪廓進(jìn)行改進(jìn)。Schaefer等[44]在總壓進(jìn)口條件下對(duì)擴(kuò)壓器子午面形狀和支撐參數(shù)進(jìn)行改進(jìn),分別如圖6、7所示。圖6中LE是支撐前緣標(biāo)志點(diǎn),TE是支撐后緣標(biāo)志點(diǎn)。工況2和工況3通道面積的變化梯度較緩TE下游通道內(nèi)流動(dòng)相對(duì)較快,將減少邊界層厚度和流動(dòng)分離,增大壓力恢復(fù)能力。在工況3中對(duì)支撐進(jìn)行改進(jìn),使支撐下游流動(dòng)軸向速度更大,使得壓力恢復(fù)進(jìn)一步提高。

Vassiliev等[45]為適應(yīng)更大的流量和功率輸出,對(duì)GT26排氣擴(kuò)壓器進(jìn)行改進(jìn)。結(jié)果如圖8所示,可以看出,新的支撐形狀明顯降低負(fù)入口氣流角工況下的流動(dòng)損失。

圖6 擴(kuò)壓通道子午線對(duì)比

圖7 原始支撐和優(yōu)化后的支撐對(duì)比

圖8 改進(jìn)設(shè)計(jì)前后擴(kuò)壓器壓力分布圖

隨算法和仿真軟件的發(fā)展,各種優(yōu)化算法被應(yīng)用到擴(kuò)壓器設(shè)計(jì)中。Cerantola等[46]采用5個(gè)自由變量定義內(nèi)錐體輪廓,2個(gè)變量定義外錐體輪廓,采用遺傳算法對(duì)這7個(gè)參數(shù)進(jìn)行優(yōu)化,利用FLUENT對(duì)優(yōu)化結(jié)果進(jìn)行驗(yàn)證。優(yōu)化結(jié)果如圖9所示,圖中為進(jìn)口半徑,o為進(jìn)口直徑,為軸向位置,最優(yōu)的擴(kuò)壓器長(zhǎng)度為外錐體直徑的0.74倍,中心輪轂為擴(kuò)張角約14°的拋物面。由此可見,聚攏式外錐體結(jié)構(gòu)可獲得更均勻的速度分布和更低的總壓損失。

王鎮(zhèn)宇等[47]用三次貝塞爾曲線描述擴(kuò)壓器傾斜壁面,利用控制點(diǎn)改變端壁型線。采用三次正交試驗(yàn)對(duì)控制點(diǎn)進(jìn)行優(yōu)化,優(yōu)化后損失系數(shù)下降10.69%,靜壓恢復(fù)系數(shù)提高1.96%,優(yōu)化后的型線及流場(chǎng)如圖10所示。

圖9 優(yōu)化后的擴(kuò)壓器壓力分布圖

圖10 優(yōu)化模型和初始模型流線圖的對(duì)比

有學(xué)者嘗試對(duì)透平和擴(kuò)壓器子午輪廓一體優(yōu)化。Christian等[48]采用遺傳算法對(duì)葉片出口角和擴(kuò)壓器子午通道幾何參數(shù)進(jìn)行優(yōu)化,選擇末級(jí)透平的效率作為優(yōu)化目標(biāo)函數(shù),結(jié)果表明:優(yōu)化后壓力恢復(fù)能力明顯提高,流動(dòng)損失顯著減少,但透平葉片出口損失基本沒有變化;進(jìn)一步優(yōu)化擴(kuò)壓器后,葉片出口損失增加,表明擴(kuò)壓器和透平末級(jí)密切相關(guān)。林智榮等[49]采用iSIGHT優(yōu)化軟件對(duì)某燃機(jī)末級(jí)透平葉片葉型和擴(kuò)壓器外形進(jìn)行聯(lián)合優(yōu)化。優(yōu)化變量為靜、動(dòng)葉葉根、葉頂相對(duì)扭轉(zhuǎn)角shstrhrt,擴(kuò)壓器外錐體3段直線的擴(kuò)張角123。優(yōu)化目標(biāo)為二者的等熵效率和擴(kuò)壓器長(zhǎng)度。CFD計(jì)算結(jié)果表明:優(yōu)化后的末級(jí)透平和擴(kuò)壓器的氣動(dòng)效率提高2.8%,擴(kuò)壓器長(zhǎng)度縮短30.1%。

4 結(jié)論與展望

重型燃機(jī)排氣溫度在500~660 ℃,排氣擴(kuò)壓器殼體因受熱發(fā)生變形,會(huì)產(chǎn)生局部溫度過高、最大變形量超出允許值等情況,從而導(dǎo)致排氣擴(kuò)壓器損壞。而排氣擴(kuò)壓器殼體溫度分布與擴(kuò)壓通道內(nèi)的流動(dòng)密切相關(guān),目前的研究多將二者分開研究,建議未來采用多場(chǎng)耦合的方法,對(duì)排氣擴(kuò)壓器的流動(dòng)、傳熱進(jìn)行綜合分析。

1)針對(duì)擴(kuò)壓器內(nèi)流動(dòng)的研究有以下特點(diǎn):現(xiàn)有對(duì)擴(kuò)壓器影響因素的研究多是針對(duì)單個(gè)因素開展的,且關(guān)注點(diǎn)主要在上游流場(chǎng)對(duì)下游擴(kuò)壓器流場(chǎng)的影響;多個(gè)影響因素對(duì)擴(kuò)壓器氣動(dòng)傳熱性能的綜合影響,擴(kuò)壓器內(nèi)部流場(chǎng)改變對(duì)上游透平葉片排的氣動(dòng)性能和流場(chǎng)分布的影響,以及透平末級(jí)和擴(kuò)壓器之間流體的相互作用的機(jī)理有待于深入研究。

2)在以往的分析中大部分都是基于定常的結(jié)果開展的,而實(shí)際的流動(dòng)是非定常的。隨著數(shù)值計(jì)算和測(cè)量技術(shù)的發(fā)展,需考慮透平與擴(kuò)壓器耦合流場(chǎng)的非定常效應(yīng),研究二者流場(chǎng)內(nèi)的二次流、渦旋等隨時(shí)間的發(fā)展效應(yīng)。

3)在擴(kuò)壓器的優(yōu)化設(shè)計(jì)問題中利用參數(shù)化方法和優(yōu)化算法可以明顯提高工作效率。其中,梯度類方法中的伴隨優(yōu)化方法,因其計(jì)算量與設(shè)計(jì)變量數(shù)目呈現(xiàn)弱相關(guān)性,可以高效地求解優(yōu)化所需的梯度信息。將伴隨方法引入內(nèi)流非定常計(jì)算可成為下一步研究的重點(diǎn)。

總之,為了設(shè)計(jì)出高性能、結(jié)構(gòu)緊湊和可靠的擴(kuò)壓器,必須把透平和擴(kuò)壓器作為一個(gè)整體,從流動(dòng)傳熱和強(qiáng)度等多學(xué)科角度系統(tǒng)考慮。

[1] SARABANAMUTTOO H I H,ROGERS G F C,COHEN H,et al.Gas turbine theory[M].Seventh Edition.United Kingdom:Pearson Education Limited,2017:185-186.

[2] 鄭露霞,張士杰,王波,等.GE 公司重型燃?xì)廨啓C(jī)透平冷氣量和燃?xì)獬鯗赝茰y(cè)[J].中國(guó)電機(jī)工程學(xué)報(bào),2019,39(23):6934-6943.

ZHENG L X,ZHANG S J,WANG B,et al.The cooling air and firing temperature estimation for GE’s heavy duty gas turbines[J].Proceedings of the CSEE,2019,39(23):6934-6943.

[3] 付經(jīng)倫,周嗣京,劉建軍.軸流透平與排汽系統(tǒng)間流場(chǎng)相互作用研究[J].工程熱物理學(xué)報(bào),2008,29(4):567-572.

FU J L,ZHOU S J,LIU J J.Investigation of interactions between axial turbine and exhaust hood [J].Journal of Engineering Thermophysics,2008,29(4):567-572.

[4] 中國(guó)東方汽輪機(jī)有限公司.日本三菱重工株式會(huì)社.東方–三菱M701F型燃?xì)廨啓C(jī)說明書[Z].成都:中國(guó)東方汽輪機(jī)有限公司,2008.

China Dongfang Turbine Co.,Ltd.,Mitsubishi Heavy Industries Corporation of Japan.Orient-Mitsubishi M701F Gas Turbine Manual[Z].Chengdu:China Dongfang Turbine Co.,Ltd.,2018.

[5] 三菱電力株式會(huì)社.M701J系列燃?xì)廨啓C(jī)[EB/OL].[2021-04-01].https://power.mhi.com/jp/sitemap.html.

Mitsubishi Electric Power Corporation:M701J series gas turbine[EB/OL].[2021-04-01].https://power.mhi.com/jp/sitemap.html.

[6] 燃機(jī)技術(shù):GE9F燃?xì)廨啓C(jī)結(jié)構(gòu)介紹[EB/OL].(2017-11-14)[2021-04-01].https://www.sohu.com/a/204410160_752692.

Gas turbine technology: GE9F gas turbine structure is introduced[EB/OL].(2017-11-14)[2021-04-01].https://www.sohu.com/a/204410160_ 752692.

[7] MIMIC D.Turbine-diffuser interaction [D].Hannover: Institutionelles Repositorium der Leibniz Universit?t Hannover,2021.

[8] THAYER E B.Evaluation of curved wall annular diffusers[J].ASME Paper,1971,71:26-35.

[9] STEFANO U,UMBERTO D.Experimental performance analysis of an annular diffuser with and without struts[J].Experimental Thermal and Fluid Science,2000,22:183-195.

[10] RUNSTADLER P W,DOLAN F X,DEAN R C.Diffuser data book[M].Creare Inc.,1975:170-186.

[11] 王廣.汽輪機(jī)軸向排氣缸擴(kuò)壓器內(nèi)流動(dòng)的數(shù)值分析與模態(tài)分析[D].上海:上海交通大學(xué),2018.

WANG G.Numerical analysis and mode decomposition of flow in axial exhaust hood diffuser of steam turbine[D].Shanghai:Shanghai Jiao Tong University,2018.

[12] STEFANO U,UMBERTO D.Flow development and turbulence length scales within an annular gas turbine exhaust diffuser[J].Experimental Thermal and Fluid Science,2000,22(1/2):55-70.

[13] PRAKASH R,SUDHAKAR P,MAHALAKSHMI N V.An experimental analysis of flow through annular diffuser with and without struts[C]//Aachen ASME 2006 Internal Combustion Engine Division Spring Technical Conference,2006:87-92.

[14] FRIC T F,VILLARREAL R,AUER R O,et al.Vortex shedding from struts in an annular exhaust diffuser [J].American Society of Mechanical Engineers,1998,120(1):186-192.

[15] 董雨軒,李志剛,李軍.支撐板型線和徑向傾斜設(shè)計(jì)對(duì)燃?xì)廨啓C(jī)排氣擴(kuò)壓器氣動(dòng)性能的影響[J].推進(jìn)技術(shù),2021, 42(6): 1245-1255.

DONG Y X,LI Z G,LI J.Effects of airfoil and radial tilt design of struts on aerodynamic performance of gas turbine exhaust diffuser[J].Journal of Propulsion Technology,2021,42(6): 1245-1255.

[16] PRADEEP A M,ROY B,VAIBHA V,et al.Study of gas turbine exhaust diffuser performance and its enhancement by shape modifications[C]//ASME Turbo Expo:Power for Land,Sea,& Air,2010.

[17] KUMAR D S,KUMAR K L,Effect of swirl on pressure recovery in annular diffusers[J].Journal of Mechanical Engineering Science,1980,22:305-313.

[18] BROWN K,GUILLOT S,NG W,et al.Experimental investigation of gas turbine axial diffuser performance: part I:parametric analysis of influential variables[C]// ASME 2020 Turbo Expo Conference (TE20),2020.

[19] VOLKMER S,HIRSCHMANN A,CASEY M,et al.The impact of a tip leakage jet on flow separation in axial gas turbine diffusers[C]//9th European Conference on Turbomachinery,F(xiàn)luid Dynamics and Thermodynamics,2011.

[20] ZIERER T.Experimental investigation of the flow in diffuser behind an axial flow compressor[J].ASME Journal of Turbomachinery,1995,117:231-239.

[21] WILLINGER R,HASELBACHER H.The role of rotor tip clearance on the aerodynamic interaction of a last gas turbine stage and exhaust diffuser[C]//ASME International Gas Turbine & Aeroengine Congress,1998.

[22] 吳飛,戴斌,艾松等.透平末級(jí)動(dòng)葉頂部間隙變化對(duì)排氣擴(kuò)壓器氣動(dòng)性能的影響[J].東方汽輪機(jī),2017(4):24-29.

WU F,DAI B,AI S.Influence of last stage rotor tip clearance' variation on aerodynamic performance of exhaust diffuser[J].Dongfang Turbine,2017(4):24-29.

[23] ZIMMERMANN C,STETTER H.Experimental determination of the flow field in the tip region of a LP-steam turbine[C]//ASME 1993 International Gas Turbine and AeroengineCongress and Exposition,1993.

[24] Farokhi S.A trade-off study of rotor tip clearance flow in a turbine/exhaust diffuser system[C]//ASME International Gas Turbine Conference & Exhibition,1987.

[25] MIHAILOWITSCH M,SCHATZ M,VOGT D M.Numerical investigations of an axial exhaust diffuser coupling the last stage of a generic gas turbine[J].Journal of Engineering for Gas Turbines and Power,2019,141(3):1025-1033.

[26] BABU M,BHATIA D,SHUKLA R K,et al.Effect of turbine tip leakage flows on exhaust diffuser performance[C]//ASME Turbo Expo: Turbine Technical Conference & Exposition,2011.

[27] VASSILIEV V,IRMISCH S,F(xiàn)LORJANCIC S.CFD analysis of industrial gas turbine exhaust diffusers[C]// ASME Turbo Expo: Power for Land, Sea, & Air. 2002.

[28] HIRSCHMANN A,VOLKMER S.The influence of the total pressure profile on the performance of axial gas turbine diffusers[C]//ASME Turbo Expo: Power for Land, Sea, & Air,2010.

[29] DAVID K,HORST S.Effect of wakes and secondary flow on re-attachment of turbine exit annular diffuser flow[J].Journal of Turbomachinery,2009,131:1-12.

[30] DAVID K,ALEXANDER W,WERNER V.Impact of gas turbine outflow on diffuser performance-a numerical study[C]//Proceedings of ASME Turbo Expo 2004,2014.

[31] OPILAT V,SEUME J R.The effect of the operating conditions of the last turbine stage on the performance of an axial exhaust diffuser[C]//ASME 2011 Turbo Expo:Turbine Technical Conference and Exposition,2011:695-704.

[32] STEVEN S J,WILLIAMS Q J.The influence of inlet conditions on the performance of annular diffusers[J].Journal of Fluids Engineering,1980,102:357-363.

[33] DESIDERI U,MANFRIDA G.Flow and turbulence survey for a model of gas turbine exhaust diffuser [C]// ASME International Gas Turbine & Aeroengine Congress & Exposition,1995.

[34] VASSILIEV V,IRMISCH S,CLARIDEG M,et al.Experimental and numerical investigation of the impact of swirl on the performance of industrial gas turbines exhaust diffusers[C]//ASME Turbo Expo, Collocated with the International Joint Power Generation Conference,2003.

[35] SONG X,GUILLOT S,NG W F,et al.An experimental investigation of the performance impact of swirl on a turbine exhaust diffuser/collector for a series of diffuser strut geometries[J].Journal of Engineering for Gas Turbines and Power,2016,138(9):1-8.

[36] 徐倩楠,吳虎.基于改進(jìn)流線曲率法的環(huán)形擴(kuò)壓器數(shù)值分析[J].航空計(jì)算技術(shù),2012,42(4):28-31.

XU Q N,WU H.Numerical analysis of annular diffusers for turbomachinery[J].Aeronautical Computing Technique,2012,42(4):28-31.

[37] 徐自榮,袁浩,初鵬,等.入口旋流角對(duì)排氣擴(kuò)壓段氣動(dòng)性能影響的數(shù)值研究[J].燃?xì)廨啓C(jī)技術(shù),2018,31(1):30-34.

XU Z R,YUAN H,CHU P,et al.Numerical investigation of the impact of inlet swirl on the performance of gas turbines exhaust diffusers[J].Gas Turbine Technology,2018,31(1):30-34.

[38] FELDCAMP G K ,BIRK A M.Strut losses in a diverging annular diffuser with swirling flow[C]//ASME Turbo Expo 2006:Power for Land,Sea and Air,2006.

[39] OLAF S,JOERG R S.Effects of rotating blade wakes on separation and pressure recovery in turbine exhaust diffusers[C]//ASME Turbo Expo 2008:Power for Land,Sea,and Air,2008.

[40] FLELGE H U,RLESS W,SEUME J.Swirl and tip leakage flow interaction with struts in axial diffusers [C]//Proceedings of ASME Turbo Expo,2002.

[41] PIETRASCH R Z,SEUME J R.Interaction between struts and swirl flow in gas turbine exhaust diffusers[J].Journal of Thermal Science,2005,14(4):314-320.

[42] SULTANIAN B K,SHINICHLRO N.Experimental and 3D CFD investigation in a gas turbine exhaust system[C]//ASME 1998 International Gas Turbine and Aeroengine Congress and Exhibition, Stockholm, Sweden,1998:98-111.

[43] MIMIC D,DRECHSEL B,HERBST F.Correlation between pressure recovery of highly loaded annular diffusers and integral stage design parameters[C]// Proceedings of ASME Turbo Expo 2017:Turbomachinery Technical Conference and Exposition.Charlotte,USA,2017.

[44] SCHAEFER P,HOFMANN W H,GIEB P.Multiobjective optimization for duct and strut design of an annular exhaust diffuser[C]//ASME Turbo Expo 2012:Turbine Technical Conference and Exposition.Copenhagen,Denmark.2012.

[45] VASSILIEV V,ROTHBRUST M,IRMISCH S.Refitting of exhaust diffuser of industrial gas turbine[C]//ASME Turbo Expo:Power for Land,Sea,and Air,2008:979-987.

[46] CERANTOLA D J,BIRK A M.Numerically optimizing an annular diffuser using a genetic algorithm with three objectives[C]//ASME Turbo Expo:Power for Land,Sea,and Air.2012:1033-1042.

[47] 王鎮(zhèn)宇,錢佳,戴韌.基于CFD和正交試驗(yàn)方法的非對(duì)稱平面擴(kuò)壓器的流動(dòng)優(yōu)化[J].動(dòng)力工程學(xué)報(bào),2010,30(5):352-356.

WANG Z Y,QIAN J,DAI R.Flow simulation and optimization of the asymmetric plane diffuser based on the CFD and orthogonal test method[J].Journal of Chinese Society of Power Engineering,2010,30(5):352-356.

[48] CHRISTIAN M,HEINRICH S.Optimization strategy for a coupled design of the last stage and the successive diffuser in a low pressure steam turbine[J].Journal of Turbomachinery,2013(135):011013-1-7.

[49] 林智榮,石坂浩一,袁新.燃?xì)馔钙侥┘?jí)葉片及擴(kuò)壓器的聯(lián)合優(yōu)化[J].工程熱物理學(xué)報(bào),2005,26(1):47-50.

LIN Z R,ISHIZAKA K,YUAN X.Combined optimization of last stage and diffuser for gas turbine [J].Journal of Engineering Thermophysics,2005,26(1):47-50.

Research Status of Gas Turbine Exhaust Diffuser

QIU Bin1,3,4,5, FU Jinglun1,2,3,4,5*

(1. Advanced Gas Turbine Laboratory, Institute of Engineering Thermophysics, Chinese Academy of Sciences, Haidian District, Beijing 100190, China; 2. Nanjing Institute of Future Energy System, Institute of Engineering Thermophysics, Chinese Academy of Sciences, Nanjing 211135, Jiangsu Province, China; 3. University of Chinese Academy of Sciences, Haidian District, Beijing 100049, China; 4. Key Laboratory of Advanced Energy and Power, Institute of Engineering Thermophysics, Chinese Academy of Sciences, Haidian District, Beijing 100190, China;5. Innovation Academy for Light-duty Gas Turbine, Chinese Academy of Sciences, Haidian District, Beijing 100190, China)

Gas turbine is the important equipment used in security and economic field of the country. The aerodynamic performance of exhaust diffuser has an important influence on the thermodynamic efficiency of gas turbine. The flow interaction between the exhaust diffuser and the upstream turbine and the complex geometrical configurations leads to the strong vortex pattern and high flow losses in the diffuser which reduces the performance of the upstream turbine and exhaust diffuser. In order to investigate the flow characteristics of diffusers and the coupling mechanism between the diffuser and the turbine, the research status of aerodynamic performance and optimal design of diffusers at home and abroad was reviewed. The prospect of research on the design and optimization of the diffuser was discussed. The interaction between the turbine and the diffuser should be fully considered in future research. Parameterization and other methods were introduced into the optimal design of the diffuser to improve the working efficiency.

gas turbine; exhaust diffuser; performance research

2021-04-25。

10.12096/j.2096-4528.pgt.21032

TK 05

國(guó)家自然科學(xué)基金項(xiàng)目(51776201)。

Project Supported by National Natural Science Foundation of China (51776201).

(責(zé)任編輯 辛培裕)

猜你喜歡
優(yōu)化影響
超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
是什么影響了滑動(dòng)摩擦力的大小
民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
哪些顧慮影響擔(dān)當(dāng)?
由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
沒錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
媽媽寶寶(2017年3期)2017-02-21 01:22:28
擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
基于Simulink的跟蹤干擾對(duì)跳頻通信的影響
主站蜘蛛池模板: 国产女人喷水视频| 亚洲国产欧洲精品路线久久| 99re在线视频观看| 国产白浆视频| 性欧美精品xxxx| 国产成人精品18| 制服丝袜一区| 亚洲视频免费在线| www.精品国产| 在线看免费无码av天堂的| 日本高清有码人妻| AV不卡无码免费一区二区三区| 国产青榴视频在线观看网站| 国产欧美日韩精品综合在线| 亚洲成人网在线观看| 五月婷婷欧美| 一区二区三区国产| 欧美不卡二区| 久久99国产综合精品1| 在线网站18禁| 日韩高清中文字幕| 综合社区亚洲熟妇p| 国产剧情伊人| 毛片大全免费观看| 久久鸭综合久久国产| 国产永久免费视频m3u8| 成人在线天堂| 日韩欧美亚洲国产成人综合| 亚洲色图欧美视频| 亚洲av无码专区久久蜜芽| 亚洲欧美成人综合| 亚洲AV永久无码精品古装片| 美女毛片在线| 亚洲精品视频免费| a网站在线观看| 国产剧情国内精品原创| 亚洲天堂免费| 丝袜高跟美脚国产1区| 无码久看视频| 国产不卡在线看| 五月天综合婷婷| 热九九精品| 亚洲精品免费网站| 亚洲高清无码久久久| 精品亚洲麻豆1区2区3区| 亚洲日韩日本中文在线| 亚洲福利一区二区三区| 国产一区二区免费播放| 国产乱码精品一区二区三区中文 | 喷潮白浆直流在线播放| 国产视频一二三区| 免费不卡在线观看av| 波多野结衣一区二区三区88| 青青操视频在线| 欧美精品一区二区三区中文字幕| 日韩色图区| 被公侵犯人妻少妇一区二区三区| 色成人综合| 国产精品爽爽va在线无码观看| 中文字幕啪啪| 日本午夜网站| 在线观看欧美国产| 国产无码在线调教| 久久这里只有精品免费| 国产成人无码久久久久毛片| 亚洲欧洲综合| 亚洲 欧美 中文 AⅤ在线视频| 日韩天堂在线观看| 久久综合九色综合97婷婷| 热九九精品| 欧美成人午夜影院| 91国内外精品自在线播放| 亚洲欧美国产五月天综合| 国产a v无码专区亚洲av| 国产黄色片在线看| 久热99这里只有精品视频6| 亚洲av综合网| 久久伊人色| 国产在线视频欧美亚综合| 国产成人成人一区二区| 99视频精品全国免费品| 日本一区高清|