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

    上下游結(jié)構(gòu)對噴嘴調(diào)節(jié)給水泵汽輪機調(diào)節(jié)閥組流動特性的影響

    2016-04-18 06:35:41劉應(yīng)征
    動力工程學(xué)報 2016年2期
    關(guān)鍵詞:計算流體力學(xué)閥門

    王 鵬, 何 勁, 劉應(yīng)征

    (1. 上海交通大學(xué) 燃氣輪機研究院, 上海 200240;

    2. 上海交通大學(xué) 機械與動力工程學(xué)院, 葉輪機械研究所, 上海 200240;

    3. 上海汽輪機有限公司, 上海 200240)

    ?

    上下游結(jié)構(gòu)對噴嘴調(diào)節(jié)給水泵汽輪機調(diào)節(jié)閥組流動特性的影響

    王鵬1,2,何勁3,劉應(yīng)征1,2

    (1. 上海交通大學(xué) 燃氣輪機研究院, 上海 200240;

    2. 上海交通大學(xué) 機械與動力工程學(xué)院, 葉輪機械研究所, 上海 200240;

    3. 上海汽輪機有限公司, 上海 200240)

    摘要:針對某噴嘴調(diào)節(jié)給水泵汽輪機調(diào)節(jié)閥組在額定工況時的三維蒸汽流動,采用計算流體力學(xué)方法對其進行數(shù)值模擬分析,以獲取各調(diào)節(jié)閥均全開時閥組上下游結(jié)構(gòu)對閥組管道內(nèi)部穩(wěn)態(tài)蒸汽流場的影響.計算中采用剪切應(yīng)力輸運模型(SST)來封閉湍流模型控制方程組,并采用多重參考坐標系法(MRF)計算獲取調(diào)節(jié)級內(nèi)動葉排通道區(qū)域的三維蒸汽流場.重點分析了各閥門通道內(nèi)的壓力損失水平、流量分配情況和過熱蒸汽在調(diào)節(jié)閥組系統(tǒng)中的流動特性,以及多閥腔室、閥門出口下游擴壓管路和彎管等部件內(nèi)部的復(fù)雜三維流動特征.結(jié)果表明:調(diào)節(jié)閥流量主要由閥門下游噴嘴數(shù)量決定;調(diào)節(jié)閥距離閥組系統(tǒng)進口的距離對調(diào)節(jié)閥通流能力具有較大影響.

    關(guān)鍵詞:給水泵汽輪機; 閥門; 蒸汽流場; 計算流體力學(xué)

    給水泵汽輪機是電站熱力循環(huán)系統(tǒng)的主要部件之一,相對于電動泵而言,可以提高總體熱力循環(huán)的效率.因鍋爐給水的要求,給水泵常常需要調(diào)整工作點,給水泵汽輪機常常處于變工況運行狀態(tài),采用噴嘴調(diào)節(jié)可以得到較高的變工況效率.在給水泵汽輪機中,用于調(diào)節(jié)控制蒸汽流量的調(diào)節(jié)閥組系統(tǒng)常由多個緊湊相連的調(diào)節(jié)閥構(gòu)成,每個調(diào)節(jié)閥都對應(yīng)一部分調(diào)節(jié)級噴嘴,其連接管道結(jié)構(gòu)復(fù)雜且難以保證對稱分布,而調(diào)節(jié)級整周流量分配嚴重不均勻會影響轉(zhuǎn)子的長期安全運行.此外,即便在各調(diào)節(jié)閥均全開的額定運行工況下,調(diào)節(jié)閥組復(fù)雜管路的存在將不可避免地造成一定的蒸汽壓力損失.因此,在閥組全開情形下,充分了解調(diào)節(jié)閥組系統(tǒng)中各調(diào)節(jié)閥內(nèi)部的蒸汽壓力損失分布以及流量分配是否均衡等信息,對閥組氣動優(yōu)化設(shè)計和調(diào)節(jié)閥相對布置有著非常重要的指導(dǎo)意義.

    電站蒸汽閥門的流動特性在國內(nèi)外都得到了較為廣泛的重視.王煒哲等[1]對某1 000 MW超超臨界汽輪機主調(diào)閥的三維蒸汽流場進行了數(shù)值模擬,發(fā)現(xiàn)在調(diào)節(jié)閥喉部位置出現(xiàn)較強的分離渦,流動較為紊亂.屠珊等[2]對某汽輪機GX-1型調(diào)節(jié)閥內(nèi)的流動特性進行了實驗和數(shù)值模擬研究,發(fā)現(xiàn)非全開工況下閥門喉部流場很不均勻且不穩(wěn)定.高怡秋等[3]采用數(shù)值模擬得到了某型汽輪發(fā)電機進汽調(diào)節(jié)閥的閥門流量特性曲線.相曉偉等[4]針對某汽輪機調(diào)節(jié)閥在全工況運行參數(shù)下的數(shù)值計算表明,在閥門中等開度時,閥碟下方所出現(xiàn)的空穴區(qū)將降低閥門通流能力,并可能導(dǎo)致閥桿振動.Yonezawa等[5]對某發(fā)電廠主蒸汽閥門進行了三維瞬態(tài)蒸汽流場的數(shù)值分析和實驗研究,給出了閥門蒸汽流態(tài)與進出口壓比和開度之間的關(guān)系.此外,Clari等[6]通過數(shù)值模擬和實驗測量某汽輪機調(diào)節(jié)閥內(nèi)的三維流場,對閥門開度以及壓比變化所引起的閥門下游管道分離流動狀態(tài)進行了分類.對于多閥組成的閥組系統(tǒng),張寶等[7]在得到汽輪機主蒸汽流量后,綜合考慮單個調(diào)節(jié)閥的流量與噴嘴出口壓力的關(guān)系,計算出閥組系統(tǒng)各調(diào)節(jié)閥的蒸汽流量大小.但這一方法僅考慮了噴嘴對調(diào)節(jié)閥流量大小的影響.對于本文所研究的給水泵汽輪機閥組系統(tǒng)而言,各調(diào)節(jié)閥自身結(jié)構(gòu)及其連接管道不同,且閥門之間存在較大的相互干擾,這將對閥組的流量分配影響很大.且大多數(shù)閥門蒸汽流場的研究都僅涉及單個調(diào)節(jié)閥部件的蒸汽流場.

    筆者針對某給水泵汽輪機調(diào)節(jié)閥組,采用計算流體力學(xué)方法對額定工況時整個閥組管道內(nèi)部的氣動特性進行了系統(tǒng)地建模和數(shù)值模擬研究.計算中充分考慮了調(diào)節(jié)級噴嘴和動葉排區(qū)域?qū)﹂y門流場的影響,得到了調(diào)節(jié)閥組系統(tǒng)內(nèi)4個閥門通道的流量分配,并對各調(diào)節(jié)閥及其連接管道內(nèi)的蒸汽流場進行了對比及分析.

    1閥組計算模型

    圖1為國內(nèi)某給水泵汽輪機進配汽系統(tǒng)示意圖,其流場計算模型見圖2.從圖1可以看出,該閥組系統(tǒng)主要由1個主蒸汽閥和4個緊湊布置的調(diào)節(jié)閥組成,其中閥1、閥3和閥4的閥座喉部直徑為200 mm(DN200),閥2的閥座喉部直徑為160 mm(DN160).4個調(diào)節(jié)閥下游管道的結(jié)構(gòu)布置非常復(fù)雜,閥1~閥3的下游管道包括擴壓段、彎管段和收縮段,閥4的下游管道則包括較長的等截面管道段、2個90°彎管段及突擴段.閥1、閥3和閥4都對應(yīng)31個噴嘴通道,閥2對應(yīng)19個噴嘴通道,總共112個有效噴嘴通道,而調(diào)節(jié)級動葉排則由72個動葉通道組成.過熱蒸汽經(jīng)主蒸汽閥(含濾網(wǎng))后進入調(diào)節(jié)閥組腔室,沿各調(diào)節(jié)閥門通道流入到對應(yīng)的調(diào)節(jié)級噴嘴組,最后共同通過調(diào)節(jié)級動葉排從系統(tǒng)出口截面流出.以閥1為例,其閥門通道內(nèi)各部件名稱如圖2(b)所示.

    圖1 給水泵汽輪機閥組進配汽系統(tǒng)示意圖

    (a)調(diào)節(jié)閥組流場計算模型(b)閥1通道示意圖

    圖2調(diào)節(jié)閥組計算模型示意圖

    Fig.2Calculation model of the control valve system

    相對而言,主蒸汽閥在整個閥組系統(tǒng)中的壓力損失較低[8],且濾網(wǎng)整流使得主蒸汽閥擴壓管道出口處的蒸汽流動較為均勻[9].因此,為了降低數(shù)值計算成本,選取流場計算分析區(qū)域為4個調(diào)節(jié)閥流道、下游連接管路和調(diào)節(jié)級(噴嘴和動葉排),如圖1中實線方框區(qū)域所示.由于各段調(diào)節(jié)級噴嘴區(qū)所對應(yīng)的上游閥門通道非對稱分布,在計算中必須對調(diào)節(jié)級進行全周期三維建模.計算進口位于主蒸汽閥擴壓管與調(diào)節(jié)閥組連接段處,而出口位于調(diào)節(jié)級后壓力級前的環(huán)形腔室內(nèi),進出口位置處蒸汽參數(shù)分布相對比較均勻.

    2數(shù)值計算方法

    采用ANSYS ICEM CFD 14.5軟件對調(diào)節(jié)閥組(即調(diào)節(jié)閥組通道區(qū)域、調(diào)節(jié)級噴嘴區(qū)域以及調(diào)節(jié)級動葉排區(qū)域)計算模型進行非結(jié)構(gòu)網(wǎng)格劃分,并對閥頭、閥座、閥后管道以及調(diào)節(jié)級區(qū)域進行局部網(wǎng)格加密處理.首先對表1中的4套網(wǎng)格進行網(wǎng)格無關(guān)性測試,測試中以閥組系統(tǒng)出口靜壓為評判標準.測試結(jié)果表明,隨著網(wǎng)格數(shù)的增加,采用網(wǎng)格4計算得到的系統(tǒng)出口靜壓已經(jīng)很小,且與運行參數(shù)(0.889 MPa)相差不大,故最終計算采用網(wǎng)格4.網(wǎng)格4所對應(yīng)的各區(qū)域網(wǎng)格劃分結(jié)果如表2所示.整個區(qū)域網(wǎng)格數(shù)近1 900萬,近壁第一層網(wǎng)格厚度為0.1 mm,以保證y+值小于300.

    表1 網(wǎng)格無關(guān)性測試

    表2 計算區(qū)域網(wǎng)格劃分結(jié)果

    該給水泵汽輪機在額定工況下運行的主要參數(shù)為:閥組系統(tǒng)進口處蒸汽滯止壓力為1.04 MPa,滯止溫度為373.2 ℃,蒸汽質(zhì)量流量為55.74 kg/s,閥組系統(tǒng)出口處蒸汽靜壓為0.889 MPa,調(diào)節(jié)級動葉轉(zhuǎn)速為4 697 r/min.計算時采用進口滯止壓力和滯止溫度作為進口邊界條件,出口處則采用流量邊界條件,將出口處靜壓值作為計算結(jié)果準確性的評判依據(jù).所有固體壁面采用無滑移條件處理.

    調(diào)節(jié)閥組結(jié)構(gòu)內(nèi)蒸汽流場非常復(fù)雜,并存在較多的分離區(qū)和漩渦區(qū).較多文獻研究中,閥門穩(wěn)態(tài)流場計算[4,10]采用了標準κ-ε模型并通過驗證表明其計算閥門流場的可行性.剪切應(yīng)力輸運湍流模型(Shear Stress Transport, SST)[11]綜合考慮了κ-ε模型和κ-ω模型的優(yōu)點,即在近壁區(qū)域采用κ-ω模型求解,而在遠離壁面的自由剪切流動區(qū)域則自動采用κ-ε模型求解,避免了κ-ε模型求解近壁區(qū)域的不穩(wěn)定性和κ-ω模型對來流的敏感性.對于穩(wěn)態(tài)流場計算,SST模型整體求解精度較高,能夠滿足工程計算要求.惠偉等[12]通過對雙彎管蝶閥管路流場的數(shù)值計算和實驗對比,發(fā)現(xiàn)SST模型的準確性較高.因此,采用SST模型來封閉控制方程組.調(diào)節(jié)級中的動靜葉區(qū)域流動轉(zhuǎn)換通過多重參考坐標系(Multiple Reference Frame, MRF)[13]法來實現(xiàn).MRF法通過將旋轉(zhuǎn)區(qū)域的靜止坐標系轉(zhuǎn)化為相對旋轉(zhuǎn)坐標系,同時在動靜區(qū)域交界面采用連續(xù)性假設(shè),最終實現(xiàn)整個區(qū)域流場的定常求解.閥組通道的來流工質(zhì)為高溫高壓過熱蒸汽,故在模擬中采用IAPWS-IF97實際蒸汽模型[9]來計算工質(zhì)熱物性.

    3計算結(jié)果與分析

    流量系數(shù)定義為閥門通道實際蒸汽質(zhì)量流量與調(diào)節(jié)閥臨界質(zhì)量流量的比值,即:

    (1)

    式中:qm,ac為CFD計算得到的閥門通道實際蒸汽質(zhì)量流量;qm,cr為調(diào)節(jié)閥臨界質(zhì)量流量.

    (2)

    為便于對比不同閥門通道的壓損大小,壓損系數(shù)定義為每個調(diào)節(jié)閥通道的局部壓損與整個閥組系統(tǒng)進出口壓損的比值:

    (3)

    通流能力系數(shù)定義為流量系數(shù)與壓損系數(shù)的比值,以表征相同壓損下閥門的通流量大小,即:

    (4)

    表3為計算得到的調(diào)節(jié)閥組內(nèi)不同截面位置處的蒸汽質(zhì)量流量和滯止壓力,以及對應(yīng)各個調(diào)節(jié)閥通道的流量系數(shù)ξ、壓損系數(shù)Δp和通流能力系數(shù)φ.

    表3 調(diào)節(jié)閥組額定工況下的計算結(jié)果

    從表3可以看出,通過閥2的蒸汽質(zhì)量流量占系統(tǒng)總質(zhì)量流量的16.6%,與其所對應(yīng)的噴嘴數(shù)占總噴嘴數(shù)的比例16.96%(19/112)很相近;通過閥1、閥3和閥4的蒸汽質(zhì)量流量占系統(tǒng)總質(zhì)量流量的83.4%,與其所對應(yīng)噴嘴數(shù)占總噴嘴數(shù)的比例83.04%(93/112)十分接近.這說明調(diào)節(jié)閥的質(zhì)量流量分配主要由每個閥門所對應(yīng)的調(diào)節(jié)級噴嘴數(shù)決定.然而,同一規(guī)格閥門(閥1、閥3和閥4)的蒸汽質(zhì)量流量存在一定差別,即閥4通道的蒸汽質(zhì)量流量比閥3通道的高2.32%,而比閥1通道的高5.85%.從閥門通流能力系數(shù)來看,閥門通流能力的大小與閥門距離系統(tǒng)入口的距離有關(guān),距離越遠則通流能力越強.

    圖3為調(diào)節(jié)閥組通道區(qū)域內(nèi)的馬赫數(shù)分布圖,圖4為對應(yīng)圖3中虛線處各調(diào)節(jié)閥的閥座下游出口截面法向速度分布圖.綜合圖3和圖4可以看出,由于閥1閥座出口布置在下游管道入口的右側(cè),則在管道左側(cè)形成較大的低速漩渦區(qū),同時由于閥1通道的右側(cè)進氣速度高于左側(cè)進氣速度,導(dǎo)致在閥1下游管道的右下方也形成了低速漩渦區(qū),使得閥1通道主要流通面積變小,流動阻力變大.相對于閥1通道,閥3通道左右兩側(cè)進氣速度較為對稱,沒有在左下角形成低速漩渦區(qū),而只是由于結(jié)構(gòu)布置原因在管道左上方形成低速漩渦區(qū).對于閥4通道而言,閥門左右兩側(cè)的進氣速度大小最為均勻,蒸汽在管道內(nèi)得到充分發(fā)展,但在2個90°彎管處出現(xiàn)了明顯的流動分離現(xiàn)象.閥2的喉部直徑較小,且對應(yīng)的噴嘴數(shù)最少,故其蒸汽質(zhì)量流量最小.且閥2的閥座布置在下游管道入口截面的中間,在管道左右兩側(cè)形成比較對稱的低速漩渦區(qū),由于閥2通道的右側(cè)進氣速度高于左側(cè),通道內(nèi)蒸汽的主流方向向左偏轉(zhuǎn).

    圖3 調(diào)節(jié)閥組通道區(qū)域內(nèi)馬赫數(shù)的分布

    圖5為調(diào)節(jié)閥組管道區(qū)域內(nèi)的三維流線分布圖.從圖5(a)可以看出,由系統(tǒng)進口流入的過熱蒸汽在此側(cè)主要通過閥1、閥2和閥3流入到閥門下游管道,主流區(qū)域呈三角形態(tài)特征并位于腔室右側(cè)底部,從而在閥組腔室左上段形成較為紊亂的低速流動區(qū)域;從圖5(b)可以看出,由系統(tǒng)進口流入的過熱蒸汽在此側(cè)有較多部分流過閥2套筒后上升并繞閥4套筒旋轉(zhuǎn)下降,最后流入閥4下游管道.

    圖6為調(diào)節(jié)閥組腔室不同截面位置處的蒸汽流動狀態(tài)分布圖,其中負值表示蒸汽流向閥門下游通道,正值表示蒸汽流向調(diào)節(jié)閥組腔室頂部,箭頭表示蒸汽流動方向.截面a靠近調(diào)節(jié)閥組腔室頂部,整體法向速度較小,漩渦主要存在于不同閥門套筒之間,主要由蒸汽繞流導(dǎo)致.截面b位于調(diào)節(jié)閥組腔室中間,其中閥1套筒周圍蒸汽法向速度為負值,即蒸汽流向閥門下游通道,而在閥2左上方有一股速度較大的上升氣流存在,這股氣流也正是之前提到的流入到閥4通道的蒸汽氣流,而漩渦主要存在于截面上閥1和閥2套筒中間以及閥3和閥4套筒中間,這主要由蒸汽繞流以及流動碰撞摻混導(dǎo)致.截面c靠近調(diào)節(jié)閥組腔室底部,大部分氣流方向為負值,流向閥門下游通道,漩渦主要集中在閥4套筒周圍,主要由氣流碰撞、摻混和卷吸所致.

    (a)閥1(b)閥2(c)閥3(d)閥4

    圖4各調(diào)節(jié)閥閥座下游出口截面法向速度分布

    Fig.4Distribution of normal velocity at downstream outlet of each control valve

    (a)前視圖(b)后視圖

    圖5調(diào)節(jié)閥組管道區(qū)域的三維流線圖

    Fig.5Three-dimensional streamline pattern in different valve passages

    (a)截面a(b)截面b(c)截面c

    圖6調(diào)節(jié)閥組腔室不同截面處的流態(tài)圖

    Fig.6Flow patterns at different cross sections in

    multi-valve chamber

    圖7為沿蒸汽主流方向觀察到的閥4彎管段不同截面位置處的速度矢量分布圖,其中各箭頭長度一致以體現(xiàn)漩渦分布特征,標號①②③④對應(yīng)圖5中的虛線截面位置.從圖7可以看出,閥4管道的來流順時針旋轉(zhuǎn),并由于該來流的旋轉(zhuǎn)性,4個截面的蒸汽主流旋轉(zhuǎn)方向為順時針;蒸汽經(jīng)過第一個彎管后在截面②處形成類似于迪恩渦[14]的一對左右反向渦;在第2個彎管的進口截面③處,內(nèi)側(cè)法向速度較高,且內(nèi)側(cè)壁面受到蒸汽沖擊;蒸汽經(jīng)過第2個彎管后同樣在截面④上形成類似于迪恩渦的一對反向渦,并由于來流的旋轉(zhuǎn)性,該對反向渦位于截面上下兩側(cè).

    (a)截面①(b)截面②

    (c)截面③(d)截面④

    圖7閥4通道彎管段不同截面處的流態(tài)圖

    Fig.7Flow patterns at different cross sections in tube

    bends of valve 4

    4結(jié)論

    (1) 各調(diào)節(jié)閥均全開的狀態(tài)下,不同調(diào)節(jié)閥通道的通流量大小主要由各閥門所對應(yīng)的調(diào)節(jié)級噴嘴數(shù)決定,但閥門的緊湊布置、閥門在腔室的位置以及閥門下游管道布置等均會影響閥門通流量大小,導(dǎo)致同種規(guī)格閥門中,閥4通道的蒸汽質(zhì)量流量比閥3通道高2.32%,而比閥1通道高5.85%.

    (2) 從通流能力系數(shù)來看,閥門距離調(diào)節(jié)閥組系統(tǒng)進口越近,通流能力系數(shù)越小,說明了在相同壓損下各閥門通道的通流能力大小與其距離調(diào)節(jié)閥組系統(tǒng)進口的遠近直接相關(guān),這是由閥門下游管道入口蒸汽速度的不均勻性造成的:距離越近,入口蒸汽速度越不均勻,從而降低了閥門通道通流能力.

    (3) 調(diào)節(jié)閥組腔室內(nèi)由于存在4個調(diào)節(jié)閥以及其對應(yīng)的套筒、閥頭等結(jié)構(gòu),使得蒸汽在其內(nèi)部發(fā)生較多的碰撞、摻混、繞流和卷吸等現(xiàn)象,最終導(dǎo)致腔室內(nèi)部流動非常紊亂,內(nèi)部不同橫截面上流動特性均不同.

    (4) 由于閥1~閥3通道下游均包含了突擴段,導(dǎo)致下游管道內(nèi)形成了較多的漩渦低速區(qū),增大了蒸汽的流動阻力;而閥4通道由于其來流環(huán)繞閥4套筒旋轉(zhuǎn)下降,使得蒸汽在閥門下游管道內(nèi)順時針旋轉(zhuǎn)下降,且蒸汽流過2處彎管后出現(xiàn)了較大區(qū)域的流動分離,并均形成了類似于迪恩渦的一對反向渦.

    參考文獻:

    [1]王煒哲,施鎏鎏,柴思敏.1 000 MW超臨界汽輪機主調(diào)閥內(nèi)流動和噪聲計算分析[J].動力工程,2007,27(3):401-405.

    WANG Weizhe, SHI Liuliu,CHAI Simin.Numerical analysis of flow field and noise generation within the mainstop and governing valve of a 1 000 MW supercritical steam turbine[J].Journal of Power Engineering,2007,27(3):401-405.

    [2]屠珊,孫弼,毛靖儒.汽輪機GX-1型調(diào)節(jié)閥流動特性的試驗與數(shù)值研究[J].西安交通大學(xué)報,2003,37(11):1124-1127.

    TU Shan, SUN Bi,MAO Jingru.Experimental investigation and numerical simulation on flow characteristics in steam turbine control valve of GX-1 type[J].Journal of Xi'an Jiaotong University,2003,37(11):1124-1127.

    [3]高怡秋,周振東,張李偉.基于CFD的蒸汽調(diào)節(jié)閥流量特性研究之一[J].汽輪機技術(shù),2011,53(5):328-330.

    GAO Yiqiu,ZHOU Zhendong,ZHANG Liwei.Investigate flux character to steam adjustive valve through CFD method[J].Turbine Technology,2011,53(5):328-330.

    [4]相曉偉, 毛靖儒,孫弼,等.汽輪機調(diào)節(jié)閥全工況三維流場特性的數(shù)值研究[J].西安交通大學(xué)學(xué)報,2006,40(3):289-293.

    XIANG Xiaowei,MAO Jingru,SUN Bi,etal.Numerical investigation of flow characteristic of control valve of steam turbine in the entire range of operating mode[J].Journal of Xi'an Jiaotong University,2006,40(3):289-293.

    [5]YONEZAWA K, OGI K, TAKINO T,etal. Experiment and numerical investigation of flow induced vibration of steam control valve[C]//Proceedings of the ASME 2010 3rd Joint US-European Fluids Engineering Summer Meeting Collocated with 8th International Conference on Nanochannels, Microchannels, and Minichannels. Montereal, Canada: ASME, 2010.

    [6]CLARI M B, POLKLAS T, JOOS F. Three-dimensional flow separations in the diffuser of a steam turbine control valve[C]//ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition. Vancouver, Canada: ASME, 2011.

    [7]張寶, 顧正皓,童小忠.汽輪機調(diào)節(jié)閥中蒸汽流量的計算[J].發(fā)電設(shè)備,2012,26(6):389-391.

    ZHANG Bao,GU Zhenghao,TONG Xiaozhong.Calculation of mainsteam flow in various control valves of steam turbine[J].Power Equipment,2012,26(6):389-391.

    [8]鄔文睿,萬津津,王煒哲,等. 600 MW亞臨界機組主調(diào)閥蒸汽流動損失數(shù)值分析[J].水動力學(xué)研究與進展A輯,2008,23(5):501-506.

    WU Wenrui,WAN Jinjin,WANG Weizhe,etal. Numerical analysis on flowing loss of steam flow in a main stop and control valve of a 600 MW sub-critical steam turbine[J].Chinese Journal of Hydrodynamics, 2008,23(5):501-506.

    [9]DOMNICK C B, BENRA F-K, DOHMEN H J,etal. Numerical investigation on the time-variant flow field and dynamic forces acting in steam turbine inlet valves[J].J Eng Gas Turbines Power,2014,137(8):081601.

    [10]相曉偉,毛靖儒,孫弼.汽輪機調(diào)節(jié)閥通流及損失特性研究[J].西安交通大學(xué)學(xué)報,2006,40(7):762-766.

    XIANG Xiaowei,MAO Jingru,SUN Bi.Investigation on flowrate and total pressure loss of flow through steam-turbine control valve[J].Journal of Xi'an Jiaotong University,2006,40(7):762-766.

    [11]MENTER F R, KUNTZ M,LANGTRY R. Ten years of industrial experience with the SST turbulence model[J]. Turbulence,Heat and Mass Transfer4, 2003, 4: 625-632.

    [12]惠偉,劉應(yīng)征,王少飛.蝶閥閥體后雙彎管道流場的數(shù)值模擬[J].動力工程學(xué)報,2014,34(6):450-457.

    HUI Wei,LIU Yingzheng,WANG Shaofei.Numerical simulation on turbulent flow in a dual-elbow channel behind butterfly valve[J].Journal of Chinese Society of Power Engineering,2014,34(6):450-457.

    [13]何煒,馬靜,王東,等.多參考坐標系法和滑移網(wǎng)格法在汽車前端進氣數(shù)值模擬中的比較[J].計算機輔助工程, 2007,16(3):96-100.

    HE Wei,MA Jing,WANG Dong,etal.Comparing MRF method with sliding mesh method for automotive front end airflow simulation[J].Computer Aided Engineering,2007,16(3):96-100.

    [14]朱輝,湛含輝,李燦,等.迪恩渦運動規(guī)律的數(shù)值研究[J].湖南工業(yè)大學(xué)學(xué)報,2009,23(3):58-62.

    ZHU Hui,ZHAN Hanhui,LI Can,etal.Numerical investigation on the characteristics of dean vortices[J].Journal of Hunan University of Technology,2009,23(3):58-62.

    Influence of Upstream/Downstream Structures on Aerodynamics of Control Valve System in Feed-water Pump Steam Turbine

    WANGPeng1,2,HEJin3,LIUYingzheng1,2

    (1. Gas Turbine Research Institute, Shanghai Jiaotong University, Shanghai 200240, China;2. Institute of Turbomachinery, School of Mechanical Engineering, Shanghai Jiaotong University,Shanghai 200240, China; 3. Shanghai Turbine Co., Ltd., Shanghai 200240, China)

    Abstract:The three-dimensional steam flow in control valve system of a nozzle-governing feed-water pump steam turbine was extensively modeled at the rated condition by using computational fluid dynamics (CFD), so as to study the influence of upstream/downstream structures on internal steam flow field in the piping of control valve system when all the valves are open. In simulation, the shear stress transport (SST) model was used to solve the closure problem of classic governing equations, while the multiple reference frame (MRF) method was used to simulate the flow field around the rotor in the governing state, with focus on the pressure loss and flow distribution in different valve passages, and on the flow characteristics of superheated steam in the control valve system; meanwhile, an analysis was conducted to the complicated three-dimensional steam flow patterns in following components, such as in the multi-valve chamber, the diffuser downstream the valve exit and in the tube bends, etc. Results show that the flow capacity of control valve is mainly determined by the quantity of downstream nozzles, which is greatly influenced by the distance from the control valve to the inlet of the valve system.

    Key words:feed-water pump steam turbine; valve; steam flow field; CFD

    文章編號:1674-7607(2016)02-0107-06

    中圖分類號:TK261

    文獻標志碼:A學(xué)科分類號:470.30

    作者簡介:王鵬(1991-),男,河南駐馬店人,博士研究生,研究方向為動力機械氣動力學(xué).劉應(yīng)征(通信作者),男,教授,博士,電話(Tel.):021-34206719;E-mail:yzliu@sjtu.edu.cn.

    收稿日期:2015-02-03

    修訂日期:2015-05-20

    猜你喜歡
    計算流體力學(xué)閥門
    美嘉諾閥門(大連)有限公司
    VANESSA始終引領(lǐng)三偏心閥門的未來發(fā)展
    裝配式玻璃鋼閥門井的研發(fā)及應(yīng)用
    煤氣與熱力(2021年3期)2021-06-09 06:16:18
    核電站閥門緊急采購系統(tǒng)的構(gòu)建
    智富時代(2018年5期)2018-07-18 17:52:04
    《工程流體力學(xué)》教學(xué)方法探討
    基于預(yù)條件技術(shù)的風力機葉片計算方法研究
    風力發(fā)電機組機艙結(jié)構(gòu)散熱性能分析方法
    汽車前端冷凝器進風不均勻性分析
    民用飛機靜壓源位置誤差修正設(shè)計研究①
    科技資訊(2015年17期)2015-10-09 21:02:59
    省力閥門瓶蓋
    久久精品国产亚洲av天美| 成人二区视频| 简卡轻食公司| 乱系列少妇在线播放| 亚洲精品中文字幕在线视频 | 欧美日本中文国产一区发布| 最新的欧美精品一区二区| 欧美精品国产亚洲| 婷婷色综合www| 欧美丝袜亚洲另类| 人妻一区二区av| 91午夜精品亚洲一区二区三区| 日韩大片免费观看网站| 国产精品久久久久久av不卡| 久久久久久久久久久丰满| 亚洲在久久综合| 纵有疾风起免费观看全集完整版| 久久久久视频综合| 中文乱码字字幕精品一区二区三区| 又爽又黄a免费视频| 免费大片黄手机在线观看| 亚洲精品aⅴ在线观看| 桃花免费在线播放| 一级毛片我不卡| 自线自在国产av| 久久鲁丝午夜福利片| 国产精品福利在线免费观看| 我要看日韩黄色一级片| 国产69精品久久久久777片| 午夜免费男女啪啪视频观看| 国产午夜精品久久久久久一区二区三区| 亚洲欧洲日产国产| 少妇熟女欧美另类| 尾随美女入室| 国产黄色免费在线视频| 色视频在线一区二区三区| 国产精品一二三区在线看| 两个人免费观看高清视频 | 久久精品夜色国产| 久久亚洲国产成人精品v| 日日啪夜夜撸| 国国产精品蜜臀av免费| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 久久国产乱子免费精品| 精品少妇久久久久久888优播| 中文乱码字字幕精品一区二区三区| 精品99又大又爽又粗少妇毛片| 日韩一区二区视频免费看| 久久久国产精品麻豆| av福利片在线| 一级毛片我不卡| videos熟女内射| 老司机影院毛片| 日韩精品有码人妻一区| 午夜福利影视在线免费观看| videos熟女内射| 久久av网站| 亚洲经典国产精华液单| 天天操日日干夜夜撸| 又大又黄又爽视频免费| av又黄又爽大尺度在线免费看| 一级毛片久久久久久久久女| 成年av动漫网址| 女人久久www免费人成看片| 26uuu在线亚洲综合色| av在线播放精品| 女的被弄到高潮叫床怎么办| 人妻制服诱惑在线中文字幕| av福利片在线| 亚洲高清免费不卡视频| 成人综合一区亚洲| 一区二区三区精品91| 在线观看国产h片| 国产乱来视频区| 国产精品不卡视频一区二区| 午夜av观看不卡| 人人妻人人添人人爽欧美一区卜| 国产91av在线免费观看| 夜夜看夜夜爽夜夜摸| 97超视频在线观看视频| 女性被躁到高潮视频| 免费人成在线观看视频色| 男的添女的下面高潮视频| 高清av免费在线| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| 亚洲激情五月婷婷啪啪| 免费观看性生交大片5| 特大巨黑吊av在线直播| 精品国产一区二区三区久久久樱花| 99热网站在线观看| a 毛片基地| 国产精品国产av在线观看| 在线亚洲精品国产二区图片欧美 | 亚洲av欧美aⅴ国产| 国产亚洲午夜精品一区二区久久| 婷婷色av中文字幕| 欧美最新免费一区二区三区| 国产男女内射视频| 99热6这里只有精品| 中文资源天堂在线| 日韩中字成人| 男人和女人高潮做爰伦理| 男女无遮挡免费网站观看| 美女大奶头黄色视频| av一本久久久久| 99热这里只有是精品在线观看| 国产午夜精品一二区理论片| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 欧美97在线视频| 天天躁夜夜躁狠狠久久av| 亚洲国产精品一区三区| 晚上一个人看的免费电影| 伦理电影大哥的女人| 美女国产视频在线观看| 亚洲国产色片| 国产在线视频一区二区| 菩萨蛮人人尽说江南好唐韦庄| 最近手机中文字幕大全| 国产成人午夜福利电影在线观看| 26uuu在线亚洲综合色| 午夜老司机福利剧场| 夜夜爽夜夜爽视频| 国产 一区精品| 国产免费视频播放在线视频| 成人综合一区亚洲| 一区二区三区精品91| 日日爽夜夜爽网站| 夜夜爽夜夜爽视频| 亚洲欧美精品专区久久| 国产成人a∨麻豆精品| 国产色爽女视频免费观看| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图 | 在线观看人妻少妇| 成人毛片60女人毛片免费| 欧美另类一区| 熟妇人妻不卡中文字幕| 亚洲av欧美aⅴ国产| 美女xxoo啪啪120秒动态图| 女的被弄到高潮叫床怎么办| 纵有疾风起免费观看全集完整版| 亚洲欧美日韩卡通动漫| 亚洲第一av免费看| 深夜a级毛片| 热re99久久国产66热| 亚洲av成人精品一区久久| 国产成人精品无人区| av一本久久久久| 亚洲国产欧美日韩在线播放 | 狂野欧美激情性bbbbbb| 色5月婷婷丁香| 最新中文字幕久久久久| 日韩 亚洲 欧美在线| 99re6热这里在线精品视频| 成人国产麻豆网| 春色校园在线视频观看| 亚洲性久久影院| 少妇人妻久久综合中文| av女优亚洲男人天堂| 久久av网站| 九九久久精品国产亚洲av麻豆| 亚洲av免费高清在线观看| 国产午夜精品一二区理论片| 亚洲av.av天堂| 在线 av 中文字幕| 日本黄色日本黄色录像| 黑人高潮一二区| 亚洲精品第二区| 看免费成人av毛片| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 亚洲精品久久久久久婷婷小说| 少妇猛男粗大的猛烈进出视频| 少妇精品久久久久久久| 久久国产乱子免费精品| 亚洲欧美中文字幕日韩二区| 国产一区二区三区综合在线观看 | 欧美日韩精品成人综合77777| 春色校园在线视频观看| 九草在线视频观看| 成人特级av手机在线观看| 亚洲国产欧美日韩在线播放 | 99热这里只有是精品在线观看| 波野结衣二区三区在线| 久久影院123| 丝瓜视频免费看黄片| 天天操日日干夜夜撸| 免费观看性生交大片5| 欧美精品高潮呻吟av久久| 免费高清在线观看视频在线观看| 春色校园在线视频观看| 岛国毛片在线播放| 成人特级av手机在线观看| 免费黄网站久久成人精品| 在线看a的网站| 日韩成人av中文字幕在线观看| 91成人精品电影| 国产精品国产三级专区第一集| 国产精品久久久久久av不卡| 国产成人一区二区在线| 午夜久久久在线观看| 天堂俺去俺来也www色官网| 亚洲av二区三区四区| 免费av中文字幕在线| 国产爽快片一区二区三区| 中文天堂在线官网| 欧美日韩国产mv在线观看视频| 国产在线男女| 亚洲一级一片aⅴ在线观看| 国产精品国产三级国产专区5o| av天堂中文字幕网| 观看美女的网站| tube8黄色片| 自拍偷自拍亚洲精品老妇| av国产久精品久网站免费入址| 少妇的逼水好多| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩另类电影网站| 久久女婷五月综合色啪小说| 亚洲情色 制服丝袜| 精品亚洲成国产av| 3wmmmm亚洲av在线观看| 久久99热6这里只有精品| 欧美区成人在线视频| 国产亚洲最大av| 日日啪夜夜爽| 黄片无遮挡物在线观看| 国产精品一区www在线观看| www.av在线官网国产| 久久久国产一区二区| 99热这里只有精品一区| 黄片无遮挡物在线观看| 在现免费观看毛片| 精品国产露脸久久av麻豆| 青春草亚洲视频在线观看| 涩涩av久久男人的天堂| av福利片在线| 亚洲成人一二三区av| av.在线天堂| 啦啦啦中文免费视频观看日本| 亚洲激情五月婷婷啪啪| 99久国产av精品国产电影| 国产探花极品一区二区| a级毛片免费高清观看在线播放| av黄色大香蕉| 最黄视频免费看| 亚洲欧美一区二区三区国产| 全区人妻精品视频| 黄色怎么调成土黄色| 日韩视频在线欧美| 亚洲伊人久久精品综合| 欧美日韩亚洲高清精品| 久久99精品国语久久久| 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线观看99| 国产亚洲最大av| 国产高清三级在线| 色婷婷av一区二区三区视频| 久久精品夜色国产| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 18+在线观看网站| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频 | 亚洲精品国产av成人精品| 欧美+日韩+精品| 啦啦啦在线观看免费高清www| 纯流量卡能插随身wifi吗| 精品人妻偷拍中文字幕| 国产精品熟女久久久久浪| 看十八女毛片水多多多| 日韩成人av中文字幕在线观看| 日韩一区二区视频免费看| 久久久久久伊人网av| 男人狂女人下面高潮的视频| 国产精品福利在线免费观看| 久久久久国产精品人妻一区二区| 夜夜看夜夜爽夜夜摸| 国精品久久久久久国模美| 国产在线男女| 97在线视频观看| tube8黄色片| 成人国产麻豆网| 在线观看www视频免费| av.在线天堂| 国产成人a∨麻豆精品| 亚洲无线观看免费| 成人18禁高潮啪啪吃奶动态图 | 内射极品少妇av片p| 日韩av免费高清视频| 亚洲av在线观看美女高潮| 亚洲精品日韩av片在线观看| 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 亚洲人与动物交配视频| 精品人妻偷拍中文字幕| 大香蕉97超碰在线| 欧美人与善性xxx| 国产精品国产三级国产av玫瑰| 丝袜在线中文字幕| 制服丝袜香蕉在线| av在线播放精品| 日韩成人伦理影院| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 亚洲欧洲精品一区二区精品久久久 | 有码 亚洲区| 少妇精品久久久久久久| 97超视频在线观看视频| tube8黄色片| 少妇人妻久久综合中文| 十八禁高潮呻吟视频 | 国产精品麻豆人妻色哟哟久久| 看非洲黑人一级黄片| 三级国产精品片| 亚洲av欧美aⅴ国产| 美女cb高潮喷水在线观看| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 国产成人免费观看mmmm| 高清在线视频一区二区三区| 视频区图区小说| 亚洲国产欧美日韩在线播放 | av黄色大香蕉| 99久久中文字幕三级久久日本| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| a级毛片在线看网站| 男人爽女人下面视频在线观看| 国产av码专区亚洲av| 好男人视频免费观看在线| av有码第一页| 久久韩国三级中文字幕| 女性生殖器流出的白浆| 亚洲婷婷狠狠爱综合网| 美女脱内裤让男人舔精品视频| 免费观看无遮挡的男女| 久久久久久久久久久免费av| 99久久综合免费| 少妇高潮的动态图| 国产精品久久久久久av不卡| 久久国产精品大桥未久av | 肉色欧美久久久久久久蜜桃| 日韩免费高清中文字幕av| 国产精品一区二区在线观看99| 久久国产亚洲av麻豆专区| 99re6热这里在线精品视频| 国产日韩一区二区三区精品不卡 | 少妇高潮的动态图| 亚洲成人手机| 久久久久久久久久人人人人人人| 插阴视频在线观看视频| 亚洲电影在线观看av| 欧美日韩国产mv在线观看视频| 老司机影院毛片| 精品卡一卡二卡四卡免费| 国产精品一区二区三区四区免费观看| 一二三四中文在线观看免费高清| 免费av不卡在线播放| 少妇猛男粗大的猛烈进出视频| av有码第一页| 亚洲国产最新在线播放| 中文字幕人妻丝袜制服| 国产片特级美女逼逼视频| 成年美女黄网站色视频大全免费 | 国产色婷婷99| 九色成人免费人妻av| 一级毛片aaaaaa免费看小| 色吧在线观看| 青春草视频在线免费观看| 久久久国产一区二区| 日韩电影二区| 97在线视频观看| 少妇的逼水好多| 成人影院久久| 亚洲无线观看免费| 国精品久久久久久国模美| 国产成人精品福利久久| 男男h啪啪无遮挡| 黄色日韩在线| 中文字幕精品免费在线观看视频 | 五月开心婷婷网| 日本av手机在线免费观看| 精品人妻熟女毛片av久久网站| 少妇裸体淫交视频免费看高清| 亚洲性久久影院| 午夜免费鲁丝| 三级经典国产精品| 一级毛片我不卡| 免费观看av网站的网址| av女优亚洲男人天堂| 久久久久久久久久久久大奶| 麻豆精品久久久久久蜜桃| 在线观看免费日韩欧美大片 | 国产欧美亚洲国产| 亚洲国产最新在线播放| 久久99一区二区三区| 在线 av 中文字幕| 这个男人来自地球电影免费观看 | 亚洲,欧美,日韩| 最近中文字幕高清免费大全6| 大香蕉久久网| 成人无遮挡网站| 国产熟女欧美一区二区| 在线 av 中文字幕| 国产极品天堂在线| 男人狂女人下面高潮的视频| 精品国产一区二区久久| 亚洲av成人精品一区久久| 久久久国产一区二区| 伊人亚洲综合成人网| 熟女人妻精品中文字幕| 亚洲精品久久久久久婷婷小说| 欧美3d第一页| 亚洲美女视频黄频| 中文字幕人妻熟人妻熟丝袜美| 国产欧美另类精品又又久久亚洲欧美| 国产av一区二区精品久久| 青春草亚洲视频在线观看| 精品少妇黑人巨大在线播放| 午夜福利视频精品| 永久网站在线| 亚洲精品第二区| 如日韩欧美国产精品一区二区三区 | 成年人午夜在线观看视频| 欧美精品一区二区大全| 久久精品久久久久久久性| 最黄视频免费看| 伊人久久国产一区二区| 高清av免费在线| 蜜桃在线观看..| 亚洲精品久久午夜乱码| 久久国产精品大桥未久av | 国产午夜精品久久久久久一区二区三区| 特大巨黑吊av在线直播| av免费在线看不卡| 久久这里有精品视频免费| 搡老乐熟女国产| 伦精品一区二区三区| 卡戴珊不雅视频在线播放| 黑丝袜美女国产一区| 在现免费观看毛片| 一区二区三区精品91| 午夜精品国产一区二区电影| 色吧在线观看| 国产精品国产av在线观看| 日本爱情动作片www.在线观看| 亚洲av二区三区四区| 亚洲av免费高清在线观看| 91午夜精品亚洲一区二区三区| 亚洲成人一二三区av| 国产伦精品一区二区三区视频9| 国产一区有黄有色的免费视频| 黑人高潮一二区| 亚洲国产毛片av蜜桃av| 亚洲国产欧美在线一区| 亚洲精华国产精华液的使用体验| 久久 成人 亚洲| 精品熟女少妇av免费看| 99热6这里只有精品| 精品久久久久久久久亚洲| av黄色大香蕉| 丰满乱子伦码专区| av有码第一页| 成年人午夜在线观看视频| 一区二区三区精品91| 免费黄网站久久成人精品| 亚洲精品色激情综合| 成人午夜精彩视频在线观看| 亚洲精品久久午夜乱码| 肉色欧美久久久久久久蜜桃| 麻豆乱淫一区二区| 九九久久精品国产亚洲av麻豆| 亚洲三级黄色毛片| 少妇裸体淫交视频免费看高清| 老司机影院毛片| a 毛片基地| 亚洲熟女精品中文字幕| 久久久a久久爽久久v久久| 一级a做视频免费观看| 精品少妇黑人巨大在线播放| 亚洲欧美清纯卡通| 亚洲精品久久久久久婷婷小说| 久久久久久久国产电影| 一级片'在线观看视频| 久久毛片免费看一区二区三区| 黄色怎么调成土黄色| 国产有黄有色有爽视频| 午夜激情福利司机影院| 一区二区av电影网| 免费人成在线观看视频色| 午夜免费鲁丝| 日本-黄色视频高清免费观看| 日韩一区二区三区影片| 国产成人91sexporn| 成人午夜精彩视频在线观看| 91久久精品国产一区二区成人| 国产av码专区亚洲av| 少妇 在线观看| 国内揄拍国产精品人妻在线| 精品一区二区三卡| 国产黄片视频在线免费观看| 九九在线视频观看精品| 亚洲精品第二区| 搡女人真爽免费视频火全软件| 成年女人在线观看亚洲视频| av国产精品久久久久影院| 夜夜看夜夜爽夜夜摸| 国产乱来视频区| 日本vs欧美在线观看视频 | 国产成人免费无遮挡视频| 高清黄色对白视频在线免费看 | 丰满人妻一区二区三区视频av| 午夜免费男女啪啪视频观看| 王馨瑶露胸无遮挡在线观看| 亚洲一区二区三区欧美精品| 亚洲精品一区蜜桃| 97在线视频观看| 久久久国产精品麻豆| 日韩伦理黄色片| 亚洲国产毛片av蜜桃av| 在线观看美女被高潮喷水网站| 久久女婷五月综合色啪小说| 少妇人妻精品综合一区二区| 久久久久精品久久久久真实原创| 久久这里有精品视频免费| 国产成人免费无遮挡视频| 中文天堂在线官网| 蜜桃久久精品国产亚洲av| 国产精品三级大全| 亚洲经典国产精华液单| 亚洲精品第二区| 色婷婷av一区二区三区视频| 亚洲精华国产精华液的使用体验| 性色avwww在线观看| 久久99精品国语久久久| 久久久国产欧美日韩av| 欧美日韩视频高清一区二区三区二| 人体艺术视频欧美日本| 久久精品久久久久久久性| 97超视频在线观看视频| 欧美激情国产日韩精品一区| 亚洲av免费高清在线观看| 亚洲精品日韩av片在线观看| 青春草视频在线免费观看| 日日爽夜夜爽网站| 久久免费观看电影| 黑人巨大精品欧美一区二区蜜桃 | 街头女战士在线观看网站| 高清在线视频一区二区三区| 丰满迷人的少妇在线观看| 99视频精品全部免费 在线| 日本黄色片子视频| 777米奇影视久久| 看免费成人av毛片| 丝袜脚勾引网站| 搡女人真爽免费视频火全软件| 少妇的逼水好多| 日韩欧美精品免费久久| 国产精品成人在线| 国产精品麻豆人妻色哟哟久久| 嫩草影院新地址| 日韩欧美 国产精品| 亚洲精品乱码久久久v下载方式| 欧美另类一区| 插逼视频在线观看| 国产一级毛片在线| 国产精品三级大全| 大陆偷拍与自拍| 建设人人有责人人尽责人人享有的| 日韩三级伦理在线观看| 免费大片黄手机在线观看| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 麻豆成人av视频| 男人爽女人下面视频在线观看| 最新中文字幕久久久久| 国产精品人妻久久久久久| 中文字幕精品免费在线观看视频 | 婷婷色麻豆天堂久久| 91成人精品电影| 国产黄片视频在线免费观看| 精品人妻偷拍中文字幕| 中文在线观看免费www的网站| 亚洲情色 制服丝袜| 国产成人精品无人区| av线在线观看网站| 久久综合国产亚洲精品| 国产女主播在线喷水免费视频网站| 久久久久久久大尺度免费视频| 丰满乱子伦码专区| av.在线天堂| 亚洲精品国产av蜜桃| 久久久精品94久久精品| 人体艺术视频欧美日本| 亚洲国产欧美日韩在线播放 | 五月开心婷婷网| 国产精品国产三级国产专区5o| 日韩三级伦理在线观看| 伦理电影大哥的女人| 国产视频内射| 国产成人精品婷婷| av天堂中文字幕网| 国产一级毛片在线| 男女国产视频网站| 寂寞人妻少妇视频99o| 亚洲精品乱码久久久v下载方式| 一级二级三级毛片免费看| 日韩欧美 国产精品| videos熟女内射|