• <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
    艦船尾跡濁度計算方法
    国产极品粉嫩免费观看在线| 一区二区三区国产精品乱码| 无人区码免费观看不卡| 18美女黄网站色大片免费观看| 国产成人精品无人区| 精品久久久久久久毛片微露脸| 天天添夜夜摸| 国产精品自产拍在线观看55亚洲| 99精国产麻豆久久婷婷| 日韩av在线大香蕉| 大型av网站在线播放| 久久香蕉精品热| 热99国产精品久久久久久7| 午夜日韩欧美国产| 国产精品国产高清国产av| 亚洲欧美一区二区三区久久| 久久精品亚洲熟妇少妇任你| 国产精品一区二区精品视频观看| 久久香蕉精品热| 亚洲成人精品中文字幕电影 | 亚洲欧美日韩另类电影网站| 亚洲精品在线美女| 国产成人欧美| 亚洲第一av免费看| 国产精品国产av在线观看| 日日干狠狠操夜夜爽| 亚洲人成伊人成综合网2020| 成人永久免费在线观看视频| 美女大奶头视频| 大陆偷拍与自拍| 日本一区二区免费在线视频| 精品高清国产在线一区| 免费看a级黄色片| 久久精品人人爽人人爽视色| 国产精品成人在线| 国产99久久九九免费精品| 亚洲成人免费电影在线观看| 一区福利在线观看| 久久中文字幕人妻熟女| av有码第一页| 国产成人一区二区三区免费视频网站| 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看 | 无限看片的www在线观看| 久久久久亚洲av毛片大全| 99在线人妻在线中文字幕| 国产精品偷伦视频观看了| 男女高潮啪啪啪动态图| 久久国产亚洲av麻豆专区| 国产aⅴ精品一区二区三区波| 日本免费a在线| 免费高清在线观看日韩| 露出奶头的视频| 久久99一区二区三区| 91国产中文字幕| 国产麻豆69| 丰满人妻熟妇乱又伦精品不卡| 国产精品乱码一区二三区的特点 | 深夜精品福利| 黄色片一级片一级黄色片| 怎么达到女性高潮| 精品电影一区二区在线| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自拍偷在线| 香蕉丝袜av| 夜夜夜夜夜久久久久| 男男h啪啪无遮挡| 999久久久国产精品视频| 老汉色av国产亚洲站长工具| 麻豆国产av国片精品| 看免费av毛片| 一边摸一边做爽爽视频免费| 亚洲熟妇熟女久久| 老司机午夜十八禁免费视频| 亚洲一区二区三区不卡视频| 日韩精品免费视频一区二区三区| 午夜精品在线福利| 亚洲色图综合在线观看| 91成人精品电影| 免费女性裸体啪啪无遮挡网站| aaaaa片日本免费| 国产亚洲精品综合一区在线观看 | 夫妻午夜视频| 999精品在线视频| 99在线人妻在线中文字幕| 午夜精品久久久久久毛片777| 国产精品野战在线观看 | 久久草成人影院| 超碰成人久久| 80岁老熟妇乱子伦牲交| 这个男人来自地球电影免费观看| 中文字幕av电影在线播放| 伊人久久大香线蕉亚洲五| 亚洲av日韩精品久久久久久密| 热99re8久久精品国产| 亚洲av片天天在线观看| 免费看a级黄色片| 中文字幕人妻丝袜制服| 免费av中文字幕在线| 午夜两性在线视频| av有码第一页| 成熟少妇高潮喷水视频| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 日本wwww免费看| 在线观看舔阴道视频| 亚洲美女黄片视频| 69av精品久久久久久| 国产精品永久免费网站| 国产99白浆流出| 欧美激情极品国产一区二区三区| 成人国语在线视频| 久久国产精品影院| 琪琪午夜伦伦电影理论片6080| 99久久国产精品久久久| 日韩免费高清中文字幕av| 国产一区在线观看成人免费| 亚洲一区二区三区欧美精品| 制服人妻中文乱码| 国产精品国产高清国产av| 亚洲精品在线观看二区| 成人av一区二区三区在线看| 亚洲激情在线av| 国产精品自产拍在线观看55亚洲| 久久性视频一级片| 中文字幕人妻熟女乱码| 精品久久久久久久久久免费视频 | 欧美亚洲日本最大视频资源| av天堂在线播放| 天堂√8在线中文| 嫩草影院精品99| 久久这里只有精品19| 人人澡人人妻人| 午夜激情av网站| 欧美成人午夜精品| 99re在线观看精品视频| 91成人精品电影| 一进一出好大好爽视频| 啦啦啦 在线观看视频| 免费看十八禁软件| 免费在线观看黄色视频的| 精品久久蜜臀av无| 热re99久久国产66热| 淫妇啪啪啪对白视频| 亚洲av片天天在线观看| 国产成+人综合+亚洲专区| 欧美日韩亚洲高清精品| 国产精品久久久人人做人人爽| 欧美日韩中文字幕国产精品一区二区三区 | 国产av又大| 久久久国产成人免费| 国产亚洲欧美98| 欧美另类亚洲清纯唯美| av电影中文网址| 又黄又爽又免费观看的视频| 亚洲 国产 在线| 老司机福利观看| 日日夜夜操网爽| 久热爱精品视频在线9| 国产又色又爽无遮挡免费看| 18禁国产床啪视频网站| 中文字幕另类日韩欧美亚洲嫩草| 日日干狠狠操夜夜爽| 一进一出抽搐动态| 啦啦啦在线免费观看视频4| 亚洲精品在线美女| 美女午夜性视频免费| 极品教师在线免费播放| 久久久久国内视频| 女生性感内裤真人,穿戴方法视频| 精品福利永久在线观看| 黄色片一级片一级黄色片| 人人妻人人添人人爽欧美一区卜| 国产亚洲av高清不卡| 国产精品国产av在线观看| 一级片免费观看大全| 国产成人啪精品午夜网站| 色婷婷久久久亚洲欧美| 成人三级黄色视频| 黄色视频,在线免费观看| 国产成人av教育| 国产成人免费无遮挡视频| 黄色 视频免费看| 天堂√8在线中文| 欧美av亚洲av综合av国产av| 国产精品综合久久久久久久免费 | 精品高清国产在线一区| 悠悠久久av| 国产91精品成人一区二区三区| 国产亚洲精品久久久久5区| 亚洲成a人片在线一区二区| 欧美乱色亚洲激情| 日韩成人在线观看一区二区三区| 国产精品九九99| 久久精品成人免费网站| 精品久久久久久电影网| 国产有黄有色有爽视频| 99精品久久久久人妻精品| 伦理电影免费视频| 国产精品免费视频内射| 80岁老熟妇乱子伦牲交| 中国美女看黄片| 正在播放国产对白刺激| 精品少妇一区二区三区视频日本电影| 一区二区日韩欧美中文字幕| 视频在线观看一区二区三区| 久久精品亚洲精品国产色婷小说| 在线天堂中文资源库| 十分钟在线观看高清视频www| 国产精品一区二区精品视频观看| videosex国产| 狂野欧美激情性xxxx| 国产精品98久久久久久宅男小说| 免费在线观看完整版高清| 中亚洲国语对白在线视频| 看免费av毛片| 成人国产一区最新在线观看| 51午夜福利影视在线观看| 国产在线精品亚洲第一网站| 涩涩av久久男人的天堂| 91国产中文字幕| 丝袜美足系列| www.999成人在线观看| 99热国产这里只有精品6| 亚洲av五月六月丁香网| 国产精品99久久99久久久不卡| 在线观看免费视频网站a站| 中文欧美无线码| 国产伦一二天堂av在线观看| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| 国产精品电影一区二区三区| 巨乳人妻的诱惑在线观看| 欧美人与性动交α欧美精品济南到| av国产精品久久久久影院| 男女下面进入的视频免费午夜 | 欧美日韩瑟瑟在线播放| 国产精品二区激情视频| 成人三级做爰电影| 夜夜看夜夜爽夜夜摸 | 亚洲欧美一区二区三区黑人| 天天躁夜夜躁狠狠躁躁| 国产不卡一卡二| av网站免费在线观看视频| 成人国语在线视频| 午夜日韩欧美国产| 中文字幕色久视频| www日本在线高清视频| 大型黄色视频在线免费观看| 午夜免费观看网址| 亚洲片人在线观看| 亚洲精品粉嫩美女一区| 51午夜福利影视在线观看| x7x7x7水蜜桃| 精品久久久久久成人av| 久久国产精品影院| 日韩 欧美 亚洲 中文字幕| 又紧又爽又黄一区二区| 午夜两性在线视频| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 国产成人一区二区三区免费视频网站| 国产成人av激情在线播放| 欧美激情久久久久久爽电影 | 老熟妇乱子伦视频在线观看| 久久久久亚洲av毛片大全| 午夜亚洲福利在线播放| 成熟少妇高潮喷水视频| 久久人人97超碰香蕉20202| 久久精品成人免费网站| 日韩国内少妇激情av| 19禁男女啪啪无遮挡网站| 中文字幕色久视频| 亚洲全国av大片| 大型av网站在线播放| netflix在线观看网站| 免费av毛片视频| 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 日韩欧美一区二区三区在线观看| 十分钟在线观看高清视频www| 悠悠久久av| 欧美激情高清一区二区三区| 亚洲av成人一区二区三| 国产蜜桃级精品一区二区三区| 狂野欧美激情性xxxx| 激情视频va一区二区三区| cao死你这个sao货| 丝袜美足系列| 免费高清视频大片| 伦理电影免费视频| 无遮挡黄片免费观看| av有码第一页| 女性生殖器流出的白浆| 男女下面进入的视频免费午夜 | 美女国产高潮福利片在线看| 国产男靠女视频免费网站| 久久精品人人爽人人爽视色| 国产精华一区二区三区| 国产精品野战在线观看 | 天堂影院成人在线观看| www.自偷自拍.com| 波多野结衣一区麻豆| 国产成+人综合+亚洲专区| 久久久久久大精品| 久久人人爽av亚洲精品天堂| 欧美午夜高清在线| 欧美一级毛片孕妇| 国产精品久久久av美女十八| bbb黄色大片| 日韩中文字幕欧美一区二区| 在线看a的网站| 亚洲一区二区三区不卡视频| 成人免费观看视频高清| 午夜免费鲁丝| 久久精品国产亚洲av香蕉五月| 自拍欧美九色日韩亚洲蝌蚪91| 999久久久国产精品视频| 如日韩欧美国产精品一区二区三区| 久久精品国产综合久久久| 国产精品一区二区三区四区久久 | 精品国产一区二区三区四区第35| 91国产中文字幕| 亚洲精品成人av观看孕妇| 欧美日本中文国产一区发布| 正在播放国产对白刺激| 亚洲专区国产一区二区| 99精品欧美一区二区三区四区| 国产高清国产精品国产三级| 日韩免费av在线播放| 国产精品成人在线| 国产精品1区2区在线观看.| а√天堂www在线а√下载| 正在播放国产对白刺激| 黄片小视频在线播放| 婷婷丁香在线五月| av欧美777| 亚洲欧美激情综合另类| 色在线成人网| 女警被强在线播放| 韩国精品一区二区三区| 最近最新免费中文字幕在线| 色老头精品视频在线观看| 女生性感内裤真人,穿戴方法视频| 人妻丰满熟妇av一区二区三区| av欧美777| 一级片'在线观看视频| 不卡一级毛片| 欧美成狂野欧美在线观看| 日韩精品中文字幕看吧| 中文字幕另类日韩欧美亚洲嫩草| 99精品在免费线老司机午夜| 亚洲av成人av| 国产av一区在线观看免费| 国产av一区二区精品久久| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 香蕉国产在线看| 国产99白浆流出| av超薄肉色丝袜交足视频| 老熟妇乱子伦视频在线观看| 性色av乱码一区二区三区2| 国产伦一二天堂av在线观看| 国产成人精品在线电影| 欧美日韩国产mv在线观看视频| 在线免费观看的www视频| 国产高清videossex| 99精品在免费线老司机午夜| 久久精品人人爽人人爽视色| 亚洲午夜理论影院| 国产精品av久久久久免费| 免费观看精品视频网站| 欧美日本中文国产一区发布| 亚洲在线自拍视频| 亚洲精品粉嫩美女一区| 国产成人av激情在线播放| 亚洲成av片中文字幕在线观看| 欧美中文综合在线视频| 另类亚洲欧美激情| 日本三级黄在线观看| 婷婷六月久久综合丁香| 国产精品亚洲av一区麻豆| 亚洲成人免费av在线播放| 久久午夜亚洲精品久久| 好男人电影高清在线观看| 如日韩欧美国产精品一区二区三区| 亚洲一区二区三区色噜噜 | 搡老乐熟女国产| 国产av在哪里看| a在线观看视频网站| 黄色a级毛片大全视频| 999久久久国产精品视频| 国产伦人伦偷精品视频| 国产av一区在线观看免费| 久久久久亚洲av毛片大全| 99国产精品一区二区蜜桃av| 可以免费在线观看a视频的电影网站| 久久精品亚洲熟妇少妇任你| 午夜福利欧美成人| 日韩av在线大香蕉| 我的亚洲天堂| 成人国产一区最新在线观看| 久久久久久久久久久久大奶| 黑人欧美特级aaaaaa片| 韩国精品一区二区三区| 一边摸一边抽搐一进一小说| 久久久久久久久免费视频了| www.自偷自拍.com| 免费看十八禁软件| 亚洲av成人一区二区三| 午夜福利,免费看| 性欧美人与动物交配| 国产精品国产高清国产av| 国产精品爽爽va在线观看网站 | 国内久久婷婷六月综合欲色啪| 自线自在国产av| 交换朋友夫妻互换小说| 国产高清videossex| 超碰97精品在线观看| 亚洲性夜色夜夜综合| 91老司机精品| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 首页视频小说图片口味搜索| 精品国产亚洲在线| 成人免费观看视频高清| 看免费av毛片| 国产主播在线观看一区二区| 大码成人一级视频| 亚洲美女黄片视频| 亚洲欧美激情综合另类| 在线观看日韩欧美| 在线观看www视频免费| 久久久国产成人免费| 999久久久精品免费观看国产| 免费在线观看黄色视频的| 悠悠久久av| 国产成人av教育| 成人国产一区最新在线观看| 国产麻豆69| 亚洲国产精品999在线| 成人国语在线视频| 精品人妻1区二区| 日韩视频一区二区在线观看| 久久性视频一级片| 中文字幕最新亚洲高清| 久久婷婷成人综合色麻豆| 亚洲熟妇熟女久久| 99精品久久久久人妻精品| 中文字幕av电影在线播放| 亚洲国产精品999在线| 国产成人影院久久av| 美女福利国产在线| 99久久久亚洲精品蜜臀av| 一边摸一边抽搐一进一小说| 男男h啪啪无遮挡| 日韩欧美一区视频在线观看| 日韩有码中文字幕| 国产免费男女视频| 在线观看一区二区三区激情| 精品免费久久久久久久清纯| 在线av久久热| 一本大道久久a久久精品| 村上凉子中文字幕在线| 精品乱码久久久久久99久播| 久久香蕉激情| 亚洲欧美一区二区三区久久| 黄网站色视频无遮挡免费观看| 国产蜜桃级精品一区二区三区| 后天国语完整版免费观看| 最新美女视频免费是黄的| 精品电影一区二区在线| 啪啪无遮挡十八禁网站| 一级毛片女人18水好多| 天天影视国产精品| 母亲3免费完整高清在线观看| 丰满人妻熟妇乱又伦精品不卡| 又黄又粗又硬又大视频| 在线观看一区二区三区激情| 国产精品久久视频播放| 伦理电影免费视频| 男人舔女人的私密视频| 日本a在线网址| 50天的宝宝边吃奶边哭怎么回事| 国产又爽黄色视频| 两性夫妻黄色片| 国产高清国产精品国产三级| 亚洲国产精品合色在线| 巨乳人妻的诱惑在线观看| 男女下面插进去视频免费观看| 亚洲欧洲精品一区二区精品久久久| av天堂在线播放| 久久精品成人免费网站| 精品卡一卡二卡四卡免费| 久久亚洲精品不卡| 免费不卡黄色视频| 精品国产亚洲在线| 天堂中文最新版在线下载| 国产日韩一区二区三区精品不卡| 日本黄色视频三级网站网址| 国产av精品麻豆| 国内久久婷婷六月综合欲色啪| 自线自在国产av| 亚洲欧美一区二区三区黑人| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩精品亚洲av| 国产区一区二久久| 亚洲 欧美 日韩 在线 免费| 大陆偷拍与自拍| 最新美女视频免费是黄的| 天天躁夜夜躁狠狠躁躁| 免费观看人在逋| 两个人看的免费小视频| 一二三四在线观看免费中文在| 国产午夜精品久久久久久| 成年女人毛片免费观看观看9| 97人妻天天添夜夜摸| 色综合站精品国产| 两性夫妻黄色片| 国产乱人伦免费视频| 99在线人妻在线中文字幕| 国产成人av激情在线播放| 男人舔女人下体高潮全视频| 两人在一起打扑克的视频| 欧美亚洲日本最大视频资源| 交换朋友夫妻互换小说| 国产成人免费无遮挡视频| 91成人精品电影| 免费搜索国产男女视频| 人成视频在线观看免费观看| www.自偷自拍.com| 黄片小视频在线播放| 色综合欧美亚洲国产小说| 亚洲精品久久成人aⅴ小说| 成人国语在线视频| 亚洲国产看品久久| 两个人看的免费小视频| 宅男免费午夜| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产精品久久久不卡| 国产熟女xx| 日本免费a在线| 757午夜福利合集在线观看| 嫩草影院精品99| 热re99久久国产66热| 俄罗斯特黄特色一大片| 国产一区二区三区视频了| 国产精华一区二区三区| 国产男靠女视频免费网站| 亚洲欧美激情综合另类| 在线视频色国产色| 十分钟在线观看高清视频www| 免费av毛片视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久av网站| 亚洲va日本ⅴa欧美va伊人久久| 天天添夜夜摸| 午夜日韩欧美国产| 最近最新免费中文字幕在线| 欧美日韩黄片免| 国产成人av教育| 天堂俺去俺来也www色官网| 淫秽高清视频在线观看| 久久久国产成人精品二区 | 国产精品一区二区免费欧美| 亚洲精品国产区一区二| 亚洲激情在线av| 99久久国产精品久久久| 亚洲va日本ⅴa欧美va伊人久久| av在线播放免费不卡| 精品熟女少妇八av免费久了| 国产一区二区三区综合在线观看| 91国产中文字幕| 色综合欧美亚洲国产小说| 欧美亚洲日本最大视频资源| 亚洲av第一区精品v没综合| 国产日韩一区二区三区精品不卡| 涩涩av久久男人的天堂| 欧美成狂野欧美在线观看| 午夜激情av网站| 久久久精品欧美日韩精品| 两个人看的免费小视频| 欧美在线黄色| 操美女的视频在线观看| 久久婷婷成人综合色麻豆| 免费搜索国产男女视频| 91在线观看av| 国产黄a三级三级三级人| 成人影院久久| 男人舔女人下体高潮全视频| 亚洲人成伊人成综合网2020| 天堂俺去俺来也www色官网| 久久精品91无色码中文字幕| 精品国产超薄肉色丝袜足j| 亚洲人成电影免费在线| 国产主播在线观看一区二区| 人人澡人人妻人| 国产成人精品在线电影| 亚洲中文字幕日韩| 国产成人一区二区三区免费视频网站| 搡老岳熟女国产| 亚洲色图综合在线观看| 欧美激情高清一区二区三区| 国产亚洲精品久久久久久毛片| 国产区一区二久久| 国产蜜桃级精品一区二区三区| 免费观看精品视频网站| 中文字幕另类日韩欧美亚洲嫩草| 妹子高潮喷水视频| 国产精品秋霞免费鲁丝片| 在线观看一区二区三区激情| 国产99久久九九免费精品| 性色av乱码一区二区三区2|