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

    基于非線性諧波法的離心壓氣機(jī)動靜干涉流動特性研究

    2025-02-28 00:00:00朱偉張巖郭敏付雪青曹曉琳劉欣源
    車用發(fā)動機(jī) 2025年1期
    關(guān)鍵詞:尾跡

    摘要: 針對高負(fù)荷、小尺寸離心壓氣機(jī)葉輪與擴(kuò)壓器之間存在強(qiáng)烈的非定常擾動問題,采用Numeca軟件對帶葉片擴(kuò)壓器的跨聲速離心壓氣機(jī)進(jìn)行數(shù)值分析,基于非線性諧波法開展動靜干涉非定常流動特性研究,探討動靜干涉非定常流動機(jī)理。通過試驗結(jié)果驗證數(shù)值設(shè)置的合理性,對比定常計算與非線性諧波法計算結(jié)果,探討壓氣機(jī)內(nèi)部流場差異,選取典型工況點(diǎn)對壓氣機(jī)葉輪和有葉擴(kuò)壓器的內(nèi)部流動及動靜干涉非定常流動進(jìn)行詳細(xì)分析。研究結(jié)果表明:非線性諧波法與定常計算結(jié)果在近失速點(diǎn)附近性能差別較大;葉輪與擴(kuò)壓器之間均存在顯著的相互干涉作用,包括葉輪與擴(kuò)壓器之間的相互勢干擾,以及葉輪尾跡對擴(kuò)壓器流動產(chǎn)生的顯著非定常影響;隨著工況惡化,葉輪出口周向速度變大,無葉區(qū)內(nèi)出現(xiàn)高速環(huán)流甚至封閉水環(huán)的現(xiàn)象,對葉輪和擴(kuò)壓器通道流動產(chǎn)生嚴(yán)重堵塞。

    關(guān)鍵詞: 離心壓氣機(jī);非線性諧波;動靜干涉;尾跡;高速環(huán)流

    DOI: 10.3969/j.issn.1001 2222.2025.01.009

    中圖分類號:TK422.5" 文獻(xiàn)標(biāo)志碼: B" 文章編號: 1001 2222(2025)01 0060 07

    離心壓氣機(jī)因其緊湊的結(jié)構(gòu)和較高的增壓比,已成為渦輪增壓器、渦輪發(fā)動機(jī)和無人機(jī)等關(guān)鍵設(shè)備的核心組成部分。當(dāng)離心壓氣機(jī)壓比大于3時,隨著輪緣速度的提高,其內(nèi)部流動為跨聲速流動,流場受到激波影響,內(nèi)部流動更為復(fù)雜,葉片尾跡在向下游傳播過程中會產(chǎn)生強(qiáng)烈摻混,形成復(fù)雜的動靜干涉非定?,F(xiàn)象。尤其是在高負(fù)荷、小尺寸的離心壓氣機(jī)中,葉輪 擴(kuò)壓器非定常相互作用更加劇烈,動靜葉之間的不匹配是高負(fù)荷離心壓氣機(jī)性能低下的常見原因1 5。因此,為提高離心壓氣機(jī)性能,掌握葉輪 擴(kuò)壓器之間的動靜干涉非定常特性具有重要的意義。

    壓氣機(jī)動靜干涉研究的方式主要是試驗研究和數(shù)值計算,從20世紀(jì)90年代開始,有學(xué)者對離心壓氣機(jī)動靜干涉現(xiàn)象進(jìn)行試驗研究。K. U. ZIEGLER等4 5試驗分析了動靜干涉對葉輪與擴(kuò)壓器流動的影響,發(fā)現(xiàn)葉輪與擴(kuò)壓器徑向距離越小,葉輪出口尾跡區(qū)在擴(kuò)壓器葉片壓力面的作用越明顯。B. WILKOSZ等6和R. KUNTE等7通過試驗發(fā)現(xiàn),動靜干涉對葉輪性能的影響較小,對擴(kuò)壓器性能影響較大。K. GALLIER等8通過試驗發(fā)現(xiàn),葉輪與擴(kuò)壓器之間無葉區(qū)的流場不均勻,擴(kuò)壓器葉片前緣沖角變化較大,引起擴(kuò)壓器葉片載荷的劇烈變化。基于工程實際需要,數(shù)值計算研究比試驗研究更加豐富,且主要采用非定常數(shù)值模擬。S. ANISH等9數(shù)值計算發(fā)現(xiàn)葉輪和擴(kuò)壓器之間的非穩(wěn)定流動具有周期性。W. N. DAWES[10發(fā)現(xiàn)影響擴(kuò)壓器不穩(wěn)定的主要因素是葉片前緣處氣流角的周向變化。劉立軍等11發(fā)現(xiàn)擴(kuò)壓器葉片前緣存在一個低速滯止區(qū),不均勻流場周向分布造成了葉輪尾跡損失增大。

    由于常規(guī)非定常模擬計算量過于龐大,而定常模擬所使用的摻混面方法未將上下游的非定常相互作用考慮在內(nèi),其存在一定模擬誤差12。因此近些年針對如何降低獲得收斂解時間以及時間精確法計算量,學(xué)術(shù)界開展了相關(guān)研究13 14。其中非線性諧波法(NLH)已在商業(yè)軟件Numeca Fine/Turbo 中實現(xiàn),并在葉輪機(jī)械算例上獲得成功的試驗驗證15 16。L. WANG等17采用NLH有效模擬了對轉(zhuǎn)壓氣機(jī)中的非定常流動,結(jié)果表明計算時間大大縮短。B. J. LIU等18采用定常、非定常和NLH方法,對比研究了跨聲速離心壓氣機(jī)的流動特性,結(jié)果顯示NLH具有計算簡單、對網(wǎng)格限制少等優(yōu)勢,相比于常規(guī)非定常模擬明顯節(jié)省了計算時間,同時結(jié)果精確度比定常模擬更高。但相比于其他類型壓氣機(jī), NLH 數(shù)值模擬離心壓氣機(jī)的研究成果相對較少,且不同壓氣機(jī)設(shè)計會導(dǎo)致動靜干涉影響存在差異。

    本研究以某帶葉片擴(kuò)壓器的跨聲速離心壓氣機(jī)為研究對象,結(jié)合臺架試驗結(jié)果,利用穩(wěn)態(tài)算法和NLH對其進(jìn)行數(shù)值研究,對比分析定常和NLH數(shù)值結(jié)果,研究其內(nèi)部流動特性,探討葉輪與有葉擴(kuò)壓器流道內(nèi)的三維非定常流態(tài)特性,研究動靜干涉非定常流動對葉片性能以及流動匹配的影響機(jī)理,為葉輪和有葉擴(kuò)壓器的設(shè)計提供參考。

    1 非線性諧波法

    非線性諧波法是由 L. HE等19在1998年提出,利用葉輪機(jī)械內(nèi)周期性流動特點(diǎn),將非定常流動分解為定常流動與線性擾動項,用諧波函數(shù)來逼近。由于動靜干涉是壓氣機(jī)流動的非定常因素,會周期性地影響葉輪和擴(kuò)壓器通道內(nèi)的流動,因此直接采用諧波函數(shù)進(jìn)行逼近,若周期性是所有擾動的特點(diǎn),那逼近時一律采用高階諧波疊加方法20。NLH可通過控制擾動源數(shù)目和傅里葉階數(shù)來控制求解精度,比非定常流場計算量低一至兩個數(shù)量級,一般2~3階諧波足夠滿足工程需求。NLH的基本假設(shè)是葉片通過頻率(BPF)引起的擾動最大,這些周期性出現(xiàn)的波動比湍流波動要大得多。守恒型變量在非定常流動中實際上就是周期性擾動與時均值的和,這些分量又可以分解成N次諧波:

    U(,t)=U()+∑U′(,t),

    U′(,t)=∑Nk=1(Uk()eIwkt+U-k()e-Iwkt)。

    式中:Uk和U-k為共軛復(fù)數(shù),其模為擾動幅值;t為時間;為位置矢量。

    在葉輪機(jī)械應(yīng)用中,第一次諧波是與葉片排葉片通過頻率相關(guān)的基頻,周期性擾動主要是由相鄰葉片排引起,這是來自上游尾跡效應(yīng)和上下游的潛在效應(yīng)。非定常擾動輸運(yùn)方程保留了基本非定常流動方程與時間平均方程之差的一階項,考慮周期擾動的前N次諧波分量,得到了擾動的每個復(fù)振幅輸運(yùn)方程。通過將一階線性化方程轉(zhuǎn)化為頻域,使諧波擾動方程與時間平均流動方程一樣與空間相關(guān),確定性應(yīng)力可以直接由這些諧波值計算,能夠解釋非定常對時均流的影響。由于將方程轉(zhuǎn)置到頻域,基于相移邊界條件,可只模擬一個葉片通道而不考慮葉片排數(shù)來模擬多級渦輪機(jī)械中的非定常流場。

    2 研究對象及數(shù)值設(shè)置

    2.1 研究對象

    離心壓氣機(jī)主要組成部件有進(jìn)氣道、離心葉輪和擴(kuò)壓器。圖1和圖2示出該離心壓氣機(jī)的三維模型和詳細(xì)的幾何參數(shù)。葉輪采用半開式結(jié)構(gòu),外徑為93 mm,包括9支主葉片及9支分流葉片;有葉擴(kuò)壓器有17支葉片,呈對稱布置,擴(kuò)壓器進(jìn)口直徑為112 mm,出口直徑為136 mm,葉片寬度為4.7 mm;從前緣至尾緣的葉頂間隙均設(shè)置為0.5 mm。

    2.2 數(shù)值設(shè)置

    由于NLH方法工作在頻域,僅需對葉輪與有葉擴(kuò)壓器進(jìn)行單通道計算,類似于定常計算。通過Navier Stokes 方程求解流體運(yùn)動,湍流模型采用S A模型,方程空間離散采用有限體積法的中心差分方法,采用多重網(wǎng)格技術(shù)加速收斂,CFL數(shù)取為3。葉輪和擴(kuò)壓器之間轉(zhuǎn)靜界面根據(jù)定常和NLH計算要求,選擇不同的混合面進(jìn)行數(shù)據(jù)傳遞,其中定常模擬采用混合平面法,NLH采用一維無反射交界面。每個擾動的諧波數(shù)設(shè)置為 3,每個葉排最大的擾動階數(shù)設(shè)為 1。葉片表面和兩邊端壁的固體壁面給定絕熱、無滑移壁面條件。

    網(wǎng)格質(zhì)量直接影響到數(shù)值解的計算精度,離心壓氣機(jī)流道特別是具有分流葉片的流道,幾何形狀復(fù)雜,加之對漩渦及激波結(jié)構(gòu)捕捉的高精度要求,物理域的網(wǎng)格劃分對壓氣機(jī)CFD結(jié)果的可信度有重大影響。因此,有必要定制出適用于高負(fù)荷壓氣機(jī)內(nèi)部流動計算的網(wǎng)格,圖3示出離心壓氣機(jī)葉輪和有葉擴(kuò)壓器的網(wǎng)格。采用Auto Grid模塊生成網(wǎng)格,葉片網(wǎng)格采用Hamp;I型拓?fù)浣Y(jié)構(gòu),B2B截面設(shè)置網(wǎng)格周期性匹配來保證周期性交界面能夠準(zhǔn)確地傳遞激波信息,葉頂間隙采用蝶型式網(wǎng)格。葉輪網(wǎng)格數(shù)為60萬,有葉擴(kuò)壓器網(wǎng)格數(shù)為40萬,可滿足計算精度要求。

    3 試驗驗證

    3.1 試驗系統(tǒng)

    試驗研究是研究離心壓氣機(jī)的重要方法,也是最接近實際運(yùn)行條件的研究方法。本試驗研究依托于中國北方發(fā)動機(jī)研究所車用動力系統(tǒng)全國重點(diǎn)實驗室的離心壓氣機(jī)試驗臺開展。壓氣機(jī)性能測試試驗臺包括如下子系統(tǒng):試驗臺控制系統(tǒng)、試驗臺測試系統(tǒng)、試驗臺燃燒室燃油供給系統(tǒng)、手動試驗備份系統(tǒng)等。試驗平臺原理21如圖4所示。

    3.2 試驗驗證

    離心壓氣機(jī)葉輪和擴(kuò)壓器試驗樣件如圖5所示,壓氣機(jī)性能仿真值與試驗值的對比如圖6所示。從流量 壓比和流量 效率特性看出,3種不同轉(zhuǎn)速工況下,壓比和效率隨流量的變化趨勢相同,相比于高轉(zhuǎn)速,低轉(zhuǎn)速條件下的特性線重合度更好。對比發(fā)現(xiàn),低轉(zhuǎn)速時CFD最大效率值比試驗結(jié)果高2%左右,壓比值基本重合;高轉(zhuǎn)速時CFD最大效率值比試驗結(jié)果低1%左右,最大壓比值偏差較小。原因可能是除了數(shù)值模擬存在的計算精度問題外,在壓氣機(jī)性能試驗中也存在測量截面和測量點(diǎn)的參數(shù)不均勻和波動現(xiàn)象,導(dǎo)致計算結(jié)果與實際情況有所差別。綜合來看,在各個轉(zhuǎn)速下,數(shù)值計算結(jié)果與試驗結(jié)果變化趨勢基本一致,驗證了數(shù)值計算的可靠性,可對壓氣機(jī)性能展開進(jìn)一步的研究。

    4 動靜干涉非定常流動特性分析

    4.1 穩(wěn)態(tài)與非線性諧波法計算結(jié)果對比

    以設(shè)計轉(zhuǎn)速95 000 r/min為例,選取近最高效率點(diǎn)進(jìn)行分析,定常及NLH計算殘差收斂曲線如圖7所示。NLH收斂曲線與定常曲線趨勢相似,流量偏高,但受模型精度影響,收斂殘差與定常曲線相差3個量級。

    對比分析多個轉(zhuǎn)速條件下定常及NLH壓氣機(jī)特性線,結(jié)果如圖8所示。由圖8可看出,在最高效率點(diǎn)附近,NLH預(yù)測的效率比定常模擬預(yù)測的效率偏高,趨勢基本一致,在背壓較高及靠近失速點(diǎn)時由于大范圍的分離流動導(dǎo)致二者產(chǎn)生偏差。

    以95 000 r/min轉(zhuǎn)速為例,對比穩(wěn)態(tài)和NLH計算的近最高效率點(diǎn)流場(如圖9所示)。由圖9可以看到,定常計算與NLH計算的流線結(jié)果基本一致,但對于葉輪出口處的流動分離,NLH比定常計算的范圍小,這是其計算效率略高于定常計算結(jié)果的主要原因。

    4.2 葉輪與擴(kuò)壓器非定常流動特性分析

    基于近失速點(diǎn)的非定常效應(yīng)更為強(qiáng)烈,也更具有代表性,因此以設(shè)計轉(zhuǎn)速95 000 r/min為例,選取近失速點(diǎn)工況進(jìn)行分析。圖10示出無葉區(qū)葉中截面通道在一個物理周期內(nèi)的靜壓分布。由圖10可看出,靜壓具有很好的周期性,表明NLH計算已經(jīng)收斂。由于該離心壓氣機(jī)葉輪長葉片、分流葉片和擴(kuò)壓器葉片數(shù)量分別為9,9,17,葉片數(shù)比例接近1∶1∶2,所以,一個周期內(nèi)葉輪通道要經(jīng)歷2個“小周期”,靜壓脈動幅值較大,約2%,表明葉輪與擴(kuò)壓器之間的相互作用比較強(qiáng)烈。

    由于NLH模擬濾除了高于三次諧波的諧波分量,因此不存在高頻分量,這表明該方法僅考慮周期流場中占主導(dǎo)地位的低頻分量。圖11所示分別為葉輪和擴(kuò)壓器中間葉高葉片表面靜壓幅值,圖中給出幅值的一次諧波、二次諧波和三次諧波,可見壓力幅值以一次諧波為主,說明低頻分量占周期流場的主導(dǎo)地位。動葉壓力面壓力幅值基本上大于吸力面

    壓力幅值,原因是通道內(nèi)激波阻止下游勢場向上游傳播,而且靜葉前緣比動葉尾緣靜壓攝動范圍和幅度要大得多,表明上游動葉尾跡對下游靜葉影響遠(yuǎn)大于下游靜葉勢場對上游葉輪的影響。

    圖12示出不同時刻葉輪氣流平均密度沿流向分布的對比。由圖12可以看到,不同時刻下葉輪前90%弦長的區(qū)域密度分布一樣,而后10%弦長的密度分布出現(xiàn)差別,越靠近葉輪后部差別越大。因此擴(kuò)壓器勢干擾的影響范圍達(dá)到了10%葉輪流向弦長。圖13示出在一個周期T內(nèi)不同時刻(0T,T/3和2T/3)葉輪葉片表面的壓力分布。由圖13可以看出,擴(kuò)壓器的勢干擾對葉輪后部負(fù)荷影響較大,且主要影響葉輪壓力面,因為擴(kuò)壓器相對于葉輪的轉(zhuǎn)動方向是從葉輪的壓力面指向葉輪的吸力面。

    綜上分析可知,擴(kuò)壓器對葉輪的勢影響比較大,從葉輪出口向上游的影響范圍達(dá)10%葉輪弦長。而且擴(kuò)壓器的勢干擾使得葉輪出口靜壓發(fā)生周期性變化,引起葉輪后部負(fù)荷和損失等也產(chǎn)生周期性的變化。

    圖14示出擴(kuò)壓器進(jìn)口葉中截面通道周期內(nèi)氣流角的時變曲線。由圖14可看出,進(jìn)口氣流角在不同時刻差別很大,最大相差約2°。

    圖15示出不同時刻擴(kuò)壓器葉中截面葉片表面的靜壓分布。全弦長范圍內(nèi)壓力面和吸力面靜壓均產(chǎn)生波動,表明葉輪對擴(kuò)壓器的影響范圍達(dá)整個弦長,葉輪對于擴(kuò)壓器的影響不僅僅是勢干擾,而且葉輪出口尾跡還受到了擴(kuò)壓器的切割并在通道中遷移,且尾跡流和附面層產(chǎn)生相互作用。

    圖 16示出不同時刻中間葉高截面的熵值分布。由圖16知,葉輪出口尾跡流動造成擴(kuò)壓器進(jìn)口損失分布不均勻。在0T時刻,擴(kuò)壓器通道前緣已經(jīng)出現(xiàn)了主流與尾跡;在T/3時刻尾跡向下移動進(jìn)入擴(kuò)壓器葉片3/4位置;在2T/3時刻,下一個葉輪通道的尾跡已到達(dá)擴(kuò)壓器葉片前緣,尾跡和主流氣流在擴(kuò)壓器內(nèi)相互摻混。因此,擴(kuò)壓器入口因上游來流不均而受到較大影響,主流和尾跡交替形成周期性非定常擾動。而且靜葉對尾跡進(jìn)行切割,尾跡與主流相互摻混,因此在尾跡內(nèi)部存在強(qiáng)烈的剪切作用,導(dǎo)致尾跡內(nèi)部產(chǎn)生嚴(yán)重的氣動損失。

    綜上分析可得,葉輪對擴(kuò)壓器的非定常影響主要是通過周期性地改變擴(kuò)壓器的進(jìn)口氣流角產(chǎn)生,而影響的大小取決于葉輪所引起擴(kuò)壓器進(jìn)口氣流角的變化幅度。葉輪對擴(kuò)壓器的非定常影響更大,其作用機(jī)制主要是尾跡干擾作用,而不僅僅是勢干擾。

    為直觀分析葉輪與擴(kuò)壓器之間動靜干涉流動狀態(tài),選取高轉(zhuǎn)速(95 000 r/min)和低轉(zhuǎn)速(50 000 r/min)兩種轉(zhuǎn)速下近最高效率點(diǎn)中間截面葉輪和擴(kuò)壓器的相對馬赫數(shù)云圖進(jìn)行分析(見圖17和圖18),選擇葉輪一個旋轉(zhuǎn)周期內(nèi)的非定常流態(tài)變化。

    在高轉(zhuǎn)速時,葉輪區(qū)域流速分布較為流暢,在長葉片壓力面近出口處存在激波現(xiàn)象,從葉輪出水邊至整個擴(kuò)壓器通道內(nèi),流動狀態(tài)變得復(fù)雜;擴(kuò)壓器壓力面存在因高速旋轉(zhuǎn)而產(chǎn)生的周期性高速環(huán)流。無葉區(qū)內(nèi)的周向速度為非對稱分布,且隨著時間的變化而變化。葉輪區(qū)域分流葉片兩邊通道內(nèi)的漩渦流動狀態(tài)各不相同,出口形成較強(qiáng)尾跡流。在0T時刻,無葉區(qū)高速環(huán)流速度由葉輪出口傳遞到擴(kuò)壓器壓力面,速度梯度方向一致;隨著轉(zhuǎn)輪的旋轉(zhuǎn),在T/3時刻,由于葉輪葉片轉(zhuǎn)到無葉高速環(huán)流區(qū),在高速旋轉(zhuǎn)作用下帶走流體,形成更復(fù)雜的流動;轉(zhuǎn)動到2T/3時刻,無葉高速環(huán)流區(qū)范圍有一定減小。

    在低轉(zhuǎn)速時,流道內(nèi)流態(tài)更為復(fù)雜,表現(xiàn)在無葉區(qū)和擴(kuò)壓器壓力面存在因高速旋轉(zhuǎn)而產(chǎn)生的周期性高速環(huán)流,嚴(yán)重堵塞通道。隨著時間變化,高速環(huán)流范圍沒有減少,如在T/3時,葉片旋轉(zhuǎn)到高速環(huán)流區(qū),流態(tài)仍沒有改變,表明葉片旋轉(zhuǎn)對無葉區(qū)內(nèi)周向速度的改變不能對無葉區(qū)流動產(chǎn)生明顯的影響。

    圖19和圖20示出在高轉(zhuǎn)速和低轉(zhuǎn)速下,近失速點(diǎn)中間截面葉輪和擴(kuò)壓器相對馬赫數(shù)云圖,選擇葉輪一個旋轉(zhuǎn)周期內(nèi)的非定常流態(tài)變化。在高轉(zhuǎn)速時,葉輪區(qū)域內(nèi)存在不同尺寸的渦流,從葉輪出水邊至整個有葉擴(kuò)壓器通道內(nèi),流動狀態(tài)較為復(fù)雜;擴(kuò)壓器壓力面靠近前緣存在高速流動,隨著時間變化無葉區(qū)內(nèi)高速區(qū)范圍先變大后變小。在0T時刻,無葉區(qū)高速環(huán)流速度梯度由葉輪出口傳遞到擴(kuò)壓器壓力面靠近前緣區(qū)域;隨著轉(zhuǎn)輪旋轉(zhuǎn),在T/3時刻,由于葉輪葉片轉(zhuǎn)到無葉區(qū)高速環(huán)流區(qū),在高速旋轉(zhuǎn)作用下帶走流體,堵塞擴(kuò)壓器通道;轉(zhuǎn)動到2T/3時刻,無葉區(qū)高速環(huán)流區(qū)范圍有一定減小。

    在低轉(zhuǎn)速時,流道內(nèi)的流態(tài)急劇惡化,表現(xiàn)在無葉區(qū)存在因高速旋轉(zhuǎn)而產(chǎn)生的封閉“水環(huán)”,通道完全被堵塞,隨著時間變化無葉區(qū)的流場流態(tài)沒有變化,阻礙了流體的流進(jìn)或流出,在葉輪和有葉擴(kuò)壓器形成多個大尺度的漩渦流。

    5 結(jié)論

    a) 壓氣機(jī)特性線結(jié)果顯示,NLH與定常模擬總體性能基本一致,在近失速點(diǎn)由于非周期性復(fù)雜流動導(dǎo)致較大偏差;

    b) 擴(kuò)壓器的勢干擾使得葉輪出口的靜壓發(fā)生周期性變化,由于擴(kuò)壓器的壓力波流動方向是從葉輪的壓力面指向吸力面,導(dǎo)致壓力面附近的流動影響較大;

    c) 葉輪對擴(kuò)壓器的非定常的影響主要是通過周期地改變擴(kuò)壓器的進(jìn)口氣流角產(chǎn)生,而且葉輪尾跡流由于強(qiáng)烈的剪切作用,使得流動損失加劇,導(dǎo)致高熵增;

    d) 高轉(zhuǎn)速時葉輪與擴(kuò)壓器之間的動靜干涉作用明顯,無葉區(qū)產(chǎn)生高速環(huán)流;低轉(zhuǎn)速時流態(tài)急劇惡化,表現(xiàn)在無葉區(qū)因高速旋轉(zhuǎn)而產(chǎn)生封閉“水環(huán)”現(xiàn)象,對葉輪和擴(kuò)壓器通道流動產(chǎn)生嚴(yán)重堵塞。

    參考文獻(xiàn):

    [1] DRING R P.Turbine rotor stator interaction[J].ASME Journal of Engineering for Power,1982,104:729 742.

    [2] 趙鵬飛.離心壓縮機(jī)擴(kuò)壓器氣動損失機(jī)理及動靜干涉現(xiàn)象研究[D].大連:大連理工大學(xué),2017.

    [3] ZIEGLER K U,GALLUS H E,NIEHUIS R.A study on impeller diffuser interaction part Ⅰ:Influence on the performance[J].Journal of Turbomachinery,2003,125(1):173 182.

    [4] ZIEGLER K U,GALLUS H E,NIEHUIS R.A study on impeller diffuser interaction part Ⅱ:Detailed flow analysis[J].Journal" of Turbomachinery,2003,125(1):183 192.

    [5] TAMAKI H.A study on matching between centrifugal compressor impeller and low solidity diffuser and its extension to vaneless diffuser[J].Journal of Engineering for Gas Turbines and Power,2018,141(4):1 42.

    [6] WILKOSZ B,SCHMIDT J,GUENTHER C,et al.Numerical and experimental comparison of a tandem and single vane deswirler used in an aero engine centrifugal compressor[J].Journal of" Turbomachinery,2014,136:041005.

    [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:031019.

    [8] GALLIER K,LAWLESS P B,F(xiàn)LEETER S.Particle imagevelocimetry characterization of high speed centrifugal compressor impeller diffuser interaction[J].Journal of Propulsion and Power,2010,26(4):784 789.

    [9] ANISH S,SITARAM N,KIM H D.A numerical study of the unsteady interaction effects on diffuser performance in a centrifugal compressor[J].Journal of Turbomachinery,2014,136:011012.

    [10] DAWES W N.Unsteady interaction of a centrifugal impeller with its vaned diffuser:Flow analysis[J].ASME Journal of Turbomachinery,1995,117(2):213 222.

    [11] 劉立軍,徐忠,張瑋.離心式壓氣機(jī)模型級內(nèi)非定常流動的數(shù)值試驗[J].航空動力學(xué)報,2002,17(1):58 64.

    [12] LEI L B W.Investigation of the unsteady flow in a counter rotating compressor using the nonlinear harmonic" method[J].International Journal of Energy Science,2012,2(5):182 188.

    [13] GAO F,MA W,ZAMBONINI G,et al.Large eddy simulation of 3 D corner separation in a linear compressor cascade[J].Physics of Fluids,2016,27(8):354 363.

    [14] 劉海霞.基于非線性諧波法的離心壓氣機(jī)性能研究[D].哈爾濱:哈爾濱工程大學(xué),2019.

    [15] VILMIN S,LORRAIN E,HIRSCH C,et al.Unsteady flow modeling across the rotor/stator interface using the nonlinear harmonic method[C].ASME Paper GT2006 90210.

    [16] VILMIN S,LORRAIN E,HIRSCH C.Application of a nonlinear harmonic method to the simulation of clocking effects[C].ASME Paper GT2009 59475.

    [17] WANG L,LIU B,YANG X D,et al.Unsteady simulation of a counter rotating research compressor using the nonlinear harmonic method[J].Proc IMechE,Part G:J Aerospace Engineering,2012(3):415 424.

    [18] LIU B J,ZHANG B,LIU Y W.Numerical investigations of impeller diffuser interactions in a transonic centrifugal compressor stage using nonlinear harmonic method[J].Proceedings of the Institution of Mechanical Engineers,Part A:Journal of Power and Energy,2014,228(8):862 877.

    [19] HE L,NING W.Efficient approach for analysis of unsteady viscous flows in turbomachines[J].AIAA Journal,1998,36(11):2005 2012.

    [20] 藥曉江,董景新,尚捷,等.非線性諧波法在葉輪機(jī)械非定常計算中的應(yīng)用[J].推進(jìn)技術(shù),2016,37(4):632 639.

    [21] 佟鼎,田紅艷,劉欣源,等.離心壓氣機(jī)非對稱有葉擴(kuò)壓器拓穩(wěn)流動方法研究[J].車用發(fā)動機(jī),2018(6):17 22.

    Study on Dynamic and Static Interaction Flow Characteristics of Centrifugal Compressor Based on Nonlinear Harmonic Method

    ZHU Wei1,ZHANG Yan1,GUO Min2,F(xiàn)U Xueqing1,CAO Xiaolin1,LIU Xinyuan1

    (1.China North Engine Research Institute(Tianjin),Tianjin 300406,China;2.Shanxi Diesel Engine Industry Co.,Ltd.,Datong 037036,China)

    Abstract: Aiming at the problem of strong unsteady disturbance between impeller and diffuser of high load and small size centrifugal compressor, numerical analysis of transonic centrifugal compressor with vane diffuser was carried out by Numeca software. Based on the nonlinear harmonic method, the unsteady flow characteristics of impeller diffuser interaction were studied, and the unsteady flow theory was discussed. The rationality of numerical setting was verified by the test results. The difference of compressor internal flow field was discussed by comparing the calculation results of steady and nonlinear harmonic method. The unsteady flow of compressor impeller and the vane diffuser were analyzed in detail by selecting typical working conditions. The results show that the performance of nonlinear harmonic method is different from that of steady calculation near the stall point. There is significant interaction between impeller and diffuser, including potential interference between impeller and diffuser and significant unsteady influence of impeller wake on diffuser flow. With the deterioration of working condition, the outlet circumferential velocity of impeller increases, and a new phenomenon of high speed circulation and even closed water ring appears in the bladeless zone, which severely blocks the flow of impeller and diffuser.

    Key words: centrifugal compressor;nonlinear harmonic method;dynamic and static interaction;wake;high speed circulation

    [編輯: 袁曉燕]

    猜你喜歡
    尾跡
    艦船尾跡光學(xué)信號異常特征貝葉斯識別方法
    浮式風(fēng)力機(jī)運(yùn)動形式對尾跡大尺度運(yùn)動的影響
    基于本征正交分解的水平軸風(fēng)力機(jī)非定常尾跡特性分析
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    兩種典型來流條件下風(fēng)力機(jī)尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    基于EEMD-Hilbert譜的渦街流量計尾跡振蕩特性
    SAR圖像艦船尾跡檢測及其討論
    粗糙海面上艦船Kelvin尾跡的電磁散射仿真
    基于改進(jìn)Hough變換的ASAR圖像船只尾跡檢測算法研究
    上海航天(2016年6期)2016-02-15 08:53:51
    艦船尾跡濁度計算方法
    岛国视频午夜一区免费看| 他把我摸到了高潮在线观看| 久久天堂一区二区三区四区| 亚洲七黄色美女视频| 亚洲精品久久国产高清桃花| 十分钟在线观看高清视频www| 亚洲精品国产一区二区精华液| 国产亚洲精品av在线| 一进一出抽搐gif免费好疼| 成人av一区二区三区在线看| 男女午夜视频在线观看| 亚洲狠狠婷婷综合久久图片| 99香蕉大伊视频| 制服人妻中文乱码| 不卡一级毛片| 免费高清视频大片| 亚洲色图 男人天堂 中文字幕| 国产亚洲欧美98| 亚洲国产精品久久男人天堂| 亚洲欧美激情在线| 一进一出好大好爽视频| 嫁个100分男人电影在线观看| www.www免费av| 免费久久久久久久精品成人欧美视频| 国产精品98久久久久久宅男小说| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看| 淫妇啪啪啪对白视频| 久久香蕉精品热| av视频免费观看在线观看| 国产私拍福利视频在线观看| 香蕉国产在线看| 天堂影院成人在线观看| 欧美精品亚洲一区二区| av在线播放免费不卡| 亚洲av电影在线进入| 欧美日韩福利视频一区二区| 黄色毛片三级朝国网站| 99国产极品粉嫩在线观看| 国产三级黄色录像| 性少妇av在线| 好看av亚洲va欧美ⅴa在| 亚洲av成人不卡在线观看播放网| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品爽爽va在线观看网站 | 国产单亲对白刺激| 国产精品国产高清国产av| 在线视频色国产色| 黄频高清免费视频| 欧美中文综合在线视频| 久久精品91蜜桃| 欧美激情 高清一区二区三区| 69av精品久久久久久| 欧美黄色淫秽网站| 免费在线观看视频国产中文字幕亚洲| 久久亚洲精品不卡| 国产精华一区二区三区| 国产99久久九九免费精品| 亚洲人成网站在线播放欧美日韩| 一二三四社区在线视频社区8| 亚洲色图 男人天堂 中文字幕| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三区在线| 天天添夜夜摸| 成人国语在线视频| 成人欧美大片| 国产色视频综合| 亚洲久久久国产精品| 男女床上黄色一级片免费看| 美女高潮到喷水免费观看| 叶爱在线成人免费视频播放| 日韩大码丰满熟妇| 国内精品久久久久久久电影| 丰满人妻熟妇乱又伦精品不卡| 真人做人爱边吃奶动态| 手机成人av网站| 嫩草影视91久久| 在线永久观看黄色视频| 久久久国产精品麻豆| 电影成人av| 精品一区二区三区四区五区乱码| 中文字幕最新亚洲高清| 亚洲va日本ⅴa欧美va伊人久久| 色哟哟哟哟哟哟| 久久中文看片网| 成人手机av| 怎么达到女性高潮| 大型黄色视频在线免费观看| 日韩欧美免费精品| www.熟女人妻精品国产| 又紧又爽又黄一区二区| 欧美日本中文国产一区发布| 亚洲免费av在线视频| 亚洲va日本ⅴa欧美va伊人久久| 一a级毛片在线观看| 亚洲黑人精品在线| 看免费av毛片| 欧美日本亚洲视频在线播放| 99国产精品一区二区蜜桃av| 黑人欧美特级aaaaaa片| www.自偷自拍.com| 高清黄色对白视频在线免费看| 色综合欧美亚洲国产小说| 亚洲色图av天堂| 亚洲激情在线av| 亚洲成av人片免费观看| 99久久99久久久精品蜜桃| 男女午夜视频在线观看| 麻豆av在线久日| 9色porny在线观看| 人人妻人人澡欧美一区二区 | 热re99久久国产66热| www.自偷自拍.com| 啪啪无遮挡十八禁网站| 日韩大尺度精品在线看网址 | 成熟少妇高潮喷水视频| 精品国产超薄肉色丝袜足j| 啦啦啦 在线观看视频| 午夜老司机福利片| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 香蕉国产在线看| 中文字幕高清在线视频| 女生性感内裤真人,穿戴方法视频| a在线观看视频网站| 50天的宝宝边吃奶边哭怎么回事| 久久人人精品亚洲av| 欧美日韩乱码在线| 在线av久久热| 亚洲国产高清在线一区二区三 | 亚洲av第一区精品v没综合| 校园春色视频在线观看| 久久久久久免费高清国产稀缺| 亚洲精品国产精品久久久不卡| aaaaa片日本免费| 天天一区二区日本电影三级 | 啦啦啦观看免费观看视频高清 | 亚洲国产精品久久男人天堂| 亚洲专区中文字幕在线| 69av精品久久久久久| 久久国产精品人妻蜜桃| 女人高潮潮喷娇喘18禁视频| 欧美在线一区亚洲| 久久国产精品人妻蜜桃| 久久国产乱子伦精品免费另类| 久久久久久久精品吃奶| 国产99久久九九免费精品| 涩涩av久久男人的天堂| 丝袜人妻中文字幕| 国产亚洲精品综合一区在线观看 | 最好的美女福利视频网| 69av精品久久久久久| 精品第一国产精品| 免费看a级黄色片| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清 | 高清在线国产一区| 亚洲少妇的诱惑av| 国产成人影院久久av| 看免费av毛片| 日本 av在线| 欧美乱妇无乱码| 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 中文字幕人妻熟女乱码| 97人妻天天添夜夜摸| 天堂√8在线中文| 女警被强在线播放| 国产精品亚洲一级av第二区| 俄罗斯特黄特色一大片| 50天的宝宝边吃奶边哭怎么回事| 99re在线观看精品视频| 三级毛片av免费| 天天一区二区日本电影三级 | 在线观看舔阴道视频| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 亚洲精品在线美女| 一夜夜www| svipshipincom国产片| 午夜视频精品福利| 久久天躁狠狠躁夜夜2o2o| 亚洲精品av麻豆狂野| 国产在线观看jvid| av福利片在线| 亚洲五月色婷婷综合| 色综合站精品国产| 最近最新中文字幕大全电影3 | 亚洲片人在线观看| 精品无人区乱码1区二区| 亚洲av五月六月丁香网| 国产av精品麻豆| 伦理电影免费视频| av中文乱码字幕在线| 久久亚洲真实| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 99香蕉大伊视频| 色老头精品视频在线观看| 欧美日本中文国产一区发布| 欧洲精品卡2卡3卡4卡5卡区| 91成人精品电影| 18禁裸乳无遮挡免费网站照片 | 午夜激情av网站| 国产一区二区激情短视频| 一卡2卡三卡四卡精品乱码亚洲| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 丁香六月欧美| 岛国视频午夜一区免费看| 亚洲成国产人片在线观看| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 精品一区二区三区av网在线观看| 国产又色又爽无遮挡免费看| 男人操女人黄网站| 国语自产精品视频在线第100页| 国产黄a三级三级三级人| 亚洲精品国产一区二区精华液| 女人被躁到高潮嗷嗷叫费观| 这个男人来自地球电影免费观看| 欧美日本视频| 日韩欧美三级三区| 校园春色视频在线观看| 久久人人精品亚洲av| 成年人黄色毛片网站| 老司机深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| 熟女少妇亚洲综合色aaa.| 一区二区三区高清视频在线| 午夜免费激情av| 亚洲一区二区三区不卡视频| 欧美激情极品国产一区二区三区| 欧美久久黑人一区二区| 久久亚洲精品不卡| 在线av久久热| 99国产极品粉嫩在线观看| 久久性视频一级片| 国语自产精品视频在线第100页| 午夜福利免费观看在线| 国产精品一区二区在线不卡| 精品国产一区二区三区四区第35| 操美女的视频在线观看| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放 | 精品高清国产在线一区| 亚洲熟妇中文字幕五十中出| 久久国产精品人妻蜜桃| 一区福利在线观看| 91麻豆av在线| 欧美av亚洲av综合av国产av| 免费在线观看亚洲国产| 国产高清视频在线播放一区| 香蕉国产在线看| 免费无遮挡裸体视频| 精品久久久久久久毛片微露脸| 亚洲情色 制服丝袜| 啦啦啦免费观看视频1| 18禁黄网站禁片午夜丰满| 成年版毛片免费区| 亚洲七黄色美女视频| 一区二区三区国产精品乱码| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| 国产三级在线视频| 性少妇av在线| 日韩欧美国产一区二区入口| 欧美日韩福利视频一区二区| 亚洲电影在线观看av| 大型av网站在线播放| 成人精品一区二区免费| 亚洲五月婷婷丁香| 国产熟女xx| 91精品三级在线观看| 美女扒开内裤让男人捅视频| 日本五十路高清| 给我免费播放毛片高清在线观看| 一边摸一边做爽爽视频免费| 中亚洲国语对白在线视频| 亚洲一区二区三区不卡视频| 大陆偷拍与自拍| 国产精品电影一区二区三区| 久久精品影院6| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 国产单亲对白刺激| 精品久久久久久,| 三级毛片av免费| 国产一区二区三区视频了| 麻豆av在线久日| 女性被躁到高潮视频| 久久久久久大精品| 精品一区二区三区av网在线观看| 久久精品aⅴ一区二区三区四区| 国产区一区二久久| 国产三级在线视频| 最新在线观看一区二区三区| 日本五十路高清| 操出白浆在线播放| 午夜福利高清视频| 在线观看免费午夜福利视频| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 757午夜福利合集在线观看| 91成年电影在线观看| 看免费av毛片| 国产一区二区三区综合在线观看| 精品一区二区三区av网在线观看| 一进一出好大好爽视频| 日韩欧美国产在线观看| 欧美老熟妇乱子伦牲交| xxx96com| 免费在线观看完整版高清| 亚洲第一电影网av| 久久天堂一区二区三区四区| 中文字幕最新亚洲高清| 1024香蕉在线观看| 久久人妻av系列| 在线永久观看黄色视频| 国产人伦9x9x在线观看| 久久中文看片网| svipshipincom国产片| 国产xxxxx性猛交| 十分钟在线观看高清视频www| 久久欧美精品欧美久久欧美| aaaaa片日本免费| 自线自在国产av| 久久久国产成人免费| 久久久久久免费高清国产稀缺| 久久久久国内视频| 欧美成人午夜精品| 青草久久国产| 老汉色∧v一级毛片| 又大又爽又粗| 夜夜看夜夜爽夜夜摸| 丝袜人妻中文字幕| 日韩国内少妇激情av| 久久中文字幕人妻熟女| 香蕉久久夜色| or卡值多少钱| 亚洲国产日韩欧美精品在线观看 | 日韩高清综合在线| 久久久久久久久免费视频了| 久久久水蜜桃国产精品网| 成年人黄色毛片网站| 久久人人爽av亚洲精品天堂| 国产亚洲av高清不卡| 免费高清在线观看日韩| 成人亚洲精品av一区二区| 动漫黄色视频在线观看| 亚洲性夜色夜夜综合| 成人18禁在线播放| 欧美+亚洲+日韩+国产| 国产成人精品无人区| 久久人人爽av亚洲精品天堂| 女人高潮潮喷娇喘18禁视频| 久久精品国产亚洲av高清一级| 9色porny在线观看| 午夜影院日韩av| 一区二区三区精品91| 欧美日本视频| 国产av一区二区精品久久| 黄色片一级片一级黄色片| 老司机在亚洲福利影院| 国产午夜精品久久久久久| 亚洲精品国产一区二区精华液| 12—13女人毛片做爰片一| 日日干狠狠操夜夜爽| 精品卡一卡二卡四卡免费| 精品久久久久久久毛片微露脸| 亚洲av成人一区二区三| 日韩欧美免费精品| 激情在线观看视频在线高清| 免费看十八禁软件| 不卡一级毛片| 欧美日韩福利视频一区二区| av天堂在线播放| 精品无人区乱码1区二区| 亚洲伊人色综图| 性少妇av在线| 美女免费视频网站| 日本在线视频免费播放| 欧美黄色片欧美黄色片| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 精品欧美一区二区三区在线| 日韩国内少妇激情av| 成人三级做爰电影| 此物有八面人人有两片| 日本 欧美在线| 国产精品电影一区二区三区| 国产精品免费一区二区三区在线| 亚洲欧美激情综合另类| 黄网站色视频无遮挡免费观看| av视频在线观看入口| 国产99久久九九免费精品| 男人舔女人下体高潮全视频| 国产区一区二久久| 亚洲精品在线观看二区| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 在线天堂中文资源库| 久久伊人香网站| 午夜福利欧美成人| 高清在线国产一区| 亚洲av成人av| 99久久99久久久精品蜜桃| 久久人妻福利社区极品人妻图片| 曰老女人黄片| 视频在线观看一区二区三区| 久久婷婷成人综合色麻豆| 操出白浆在线播放| 精品久久久久久,| 精品日产1卡2卡| 久热这里只有精品99| 久久午夜综合久久蜜桃| 国产欧美日韩精品亚洲av| 亚洲色图综合在线观看| 亚洲中文av在线| netflix在线观看网站| 午夜亚洲福利在线播放| 一二三四在线观看免费中文在| 亚洲五月婷婷丁香| 91字幕亚洲| aaaaa片日本免费| 日本欧美视频一区| 久久影院123| 天堂动漫精品| 国产亚洲精品av在线| 亚洲七黄色美女视频| 亚洲欧美精品综合一区二区三区| 91国产中文字幕| 99久久久亚洲精品蜜臀av| 18禁黄网站禁片午夜丰满| 国产亚洲精品第一综合不卡| 日韩欧美在线二视频| 岛国视频午夜一区免费看| 91大片在线观看| av视频免费观看在线观看| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看 | 免费在线观看视频国产中文字幕亚洲| 无遮挡黄片免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 免费无遮挡裸体视频| av在线天堂中文字幕| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 精品国产一区二区三区四区第35| 国产单亲对白刺激| 婷婷精品国产亚洲av在线| 国产成人影院久久av| 九色亚洲精品在线播放| 男女之事视频高清在线观看| 精品少妇一区二区三区视频日本电影| 日韩有码中文字幕| 日韩成人在线观看一区二区三区| 久久伊人香网站| 亚洲av成人av| 一本综合久久免费| 一级片免费观看大全| 亚洲中文字幕一区二区三区有码在线看 | 欧美日本视频| 亚洲精华国产精华精| 国产成人欧美| 国产97色在线日韩免费| 少妇被粗大的猛进出69影院| 久久 成人 亚洲| 香蕉丝袜av| 777久久人妻少妇嫩草av网站| 大香蕉久久成人网| 黑人巨大精品欧美一区二区蜜桃| 美女扒开内裤让男人捅视频| 一二三四社区在线视频社区8| 咕卡用的链子| 久久久久国内视频| 此物有八面人人有两片| 午夜久久久在线观看| 中文字幕人妻丝袜一区二区| 国产伦一二天堂av在线观看| www.www免费av| 日本黄色视频三级网站网址| a在线观看视频网站| 亚洲欧美激情在线| 日韩大码丰满熟妇| 女同久久另类99精品国产91| 黄频高清免费视频| 精品国产国语对白av| 女生性感内裤真人,穿戴方法视频| 怎么达到女性高潮| 99国产精品免费福利视频| 久久影院123| 国产精品,欧美在线| 伦理电影免费视频| 欧美午夜高清在线| 亚洲第一电影网av| 黄色女人牲交| 午夜成年电影在线免费观看| 国产主播在线观看一区二区| 国产精品久久电影中文字幕| 99re在线观看精品视频| 人妻久久中文字幕网| 中文字幕人成人乱码亚洲影| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| 久久久久久亚洲精品国产蜜桃av| 最好的美女福利视频网| 女人精品久久久久毛片| 亚洲激情在线av| 在线国产一区二区在线| 麻豆av在线久日| 嫩草影视91久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文av在线| 欧美不卡视频在线免费观看 | 久久久久国内视频| 久久久久亚洲av毛片大全| 色尼玛亚洲综合影院| 琪琪午夜伦伦电影理论片6080| 色老头精品视频在线观看| 久久久国产成人免费| 国产真人三级小视频在线观看| 精品国产乱子伦一区二区三区| 男女之事视频高清在线观看| 身体一侧抽搐| 成年女人毛片免费观看观看9| 国产成年人精品一区二区| 无限看片的www在线观看| 青草久久国产| 亚洲国产精品999在线| 精品一区二区三区av网在线观看| 国产亚洲精品久久久久5区| 久久精品91蜜桃| 性欧美人与动物交配| 日本精品一区二区三区蜜桃| 欧美性长视频在线观看| 如日韩欧美国产精品一区二区三区| 村上凉子中文字幕在线| 精品一区二区三区av网在线观看| 亚洲欧美一区二区三区黑人| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 亚洲精品国产精品久久久不卡| 超碰成人久久| 亚洲五月天丁香| 亚洲 欧美 日韩 在线 免费| 男人操女人黄网站| 久久久久亚洲av毛片大全| 人妻久久中文字幕网| 人人妻人人爽人人添夜夜欢视频| 中文字幕av电影在线播放| 国产精品影院久久| 亚洲av第一区精品v没综合| 午夜福利影视在线免费观看| 大陆偷拍与自拍| 手机成人av网站| 大香蕉久久成人网| 每晚都被弄得嗷嗷叫到高潮| 久久久精品国产亚洲av高清涩受| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 欧美在线黄色| 淫妇啪啪啪对白视频| 国产av精品麻豆| 国产乱人伦免费视频| av网站免费在线观看视频| 亚洲欧美日韩另类电影网站| 99精品欧美一区二区三区四区| 成人免费观看视频高清| 国产人伦9x9x在线观看| 免费在线观看完整版高清| 亚洲国产中文字幕在线视频| 美女高潮喷水抽搐中文字幕| xxx96com| 亚洲一区二区三区色噜噜| 欧美日韩福利视频一区二区| 中文字幕人成人乱码亚洲影| 丰满人妻熟妇乱又伦精品不卡| 中文字幕人妻丝袜一区二区| 中文字幕久久专区| 极品教师在线免费播放| 久久香蕉精品热| 久9热在线精品视频| 免费少妇av软件| 免费观看精品视频网站| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 午夜老司机福利片| 日韩精品免费视频一区二区三区| 午夜免费鲁丝| 亚洲九九香蕉| 亚洲欧美激情在线| 国产97色在线日韩免费| 亚洲激情在线av| 亚洲av片天天在线观看| 免费看十八禁软件| 国产男靠女视频免费网站| 99久久国产精品久久久| 日韩精品青青久久久久久| 午夜老司机福利片| 欧美不卡视频在线免费观看 | 夜夜夜夜夜久久久久| 12—13女人毛片做爰片一| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 亚洲中文av在线| 免费一级毛片在线播放高清视频 | av天堂久久9| 精品乱码久久久久久99久播| 我的亚洲天堂| 黄色视频不卡|