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

    側(cè)式進(jìn)/出水口偏流出流下的大渦模擬

    2024-04-04 16:12:25郭港歸劉亞坤魏杰曹澤張帝
    水科學(xué)進(jìn)展 2024年1期
    關(guān)鍵詞:出水口

    郭港歸 劉亞坤 魏杰 曹澤 張帝

    摘要:實(shí)際工程中受地形地質(zhì)條件影響,引水洞在水平上存在轉(zhuǎn)角,導(dǎo)致出流工況來流不均,進(jìn)/出水口內(nèi)部紊流特性的復(fù)雜程度顯著提升。利用大渦模擬對某帶水平彎段抽水蓄能電站側(cè)式進(jìn)/出水口進(jìn)行數(shù)值計(jì)算,其流速、概率密度分布與模型試驗(yàn)吻合較好。結(jié)果表明:偏流出流下進(jìn)/出水口各流道分流比分別為0.64、0.81、1.26和1.29,水平方向流速分布極為不均;垂直方向主流靠近中下部,垂向雷諾切應(yīng)力在擴(kuò)散段內(nèi)呈一正一負(fù)峰值分布,該現(xiàn)象主要由中上部的流動(dòng)分離和底部的壁面剪切造成;兩中間流道的回流區(qū)高度大于兩邊流道,導(dǎo)致中間流道的攔污柵更易受到反向流速影響;流動(dòng)分離使攔污柵附近存在三軸漩滾,靠近過水?dāng)嗝嬷猩喜亢偷撞慨a(chǎn)生2處能量集中的低頻脈動(dòng),且各流道在偏流條件下的紊動(dòng)強(qiáng)度相比均勻來流的紊動(dòng)強(qiáng)度分別提高11%、25%、29%、3%,不利的水流流態(tài)和較高的紊動(dòng)強(qiáng)度可能對攔污柵造成威脅。

    關(guān)鍵詞:側(cè)式進(jìn)/出水口;偏流出流;雷諾切應(yīng)力;流動(dòng)分離;低頻脈動(dòng)

    中圖分類號(hào):TV135??文獻(xiàn)標(biāo)志碼:A??文章編號(hào):1001-6791(2024)01-0132-13

    抽水蓄能電站是一種清潔、環(huán)保、可持續(xù)的儲(chǔ)能發(fā)電方式,目前在中國“雙碳”的戰(zhàn)略背景下迅速發(fā)展。進(jìn)/出水口是抽水蓄能電站引水系統(tǒng)中的重要組成部分,其幾何結(jié)構(gòu)存在水平、垂直雙向擴(kuò)散,造成2個(gè)方向上均發(fā)生流動(dòng)分離,三維流動(dòng)特性明顯,內(nèi)部流場結(jié)構(gòu)復(fù)雜、紊亂,需要重點(diǎn)關(guān)注。孫雙科等[1]對抽水蓄能電站進(jìn)/出水口進(jìn)行模型試驗(yàn),研究了不同斷面的流速分布,發(fā)現(xiàn)隧洞坡度是造成各斷面流速分布不佳的主要原因;王晨茜等[2]利用Realizable k-ε模型對進(jìn)/出水口進(jìn)行數(shù)值模擬,分析垂直擴(kuò)散角對擴(kuò)散段流動(dòng)分離的影響,發(fā)現(xiàn)當(dāng)垂直擴(kuò)散角逐漸減小,流動(dòng)分離減弱,該角度小于2°時(shí),無明顯流動(dòng)分離現(xiàn)象;高學(xué)平等[3]利用粒子圖像測速技術(shù)(PIV)對進(jìn)/出水口出流工況內(nèi)部流場進(jìn)行測量,發(fā)現(xiàn)擴(kuò)散段存在流動(dòng)分離,頂部存在反向流速;朱洪濤等[4]采用激光多普勒流速儀得到隧洞段、擴(kuò)散段、調(diào)整段和防渦梁段的紊動(dòng)強(qiáng)度,發(fā)現(xiàn)紊動(dòng)強(qiáng)度沿程先增大后減??;Ye等[5]采用數(shù)值模擬方法對進(jìn)/出水口出流工況進(jìn)行模擬,發(fā)現(xiàn)擴(kuò)散段頂部出現(xiàn)大范圍的回流區(qū),其主要原因是主流在橫向上擴(kuò)散。大部分文獻(xiàn)主要針對引水洞段為直線管道的進(jìn)/出水口水力特性研究[6-7],而實(shí)際工程中由于地形地質(zhì)的原因需要將引水洞設(shè)計(jì)為垂直或水平方向上偏轉(zhuǎn)的體型。張?zhí)m丁[8]對帶上彎段引水洞的響水澗抽水蓄能電站進(jìn)/出水口進(jìn)行試驗(yàn)研究,發(fā)現(xiàn)彎道效應(yīng)將水流推向邊壁;當(dāng)過渡段長達(dá)到5倍隧洞洞徑時(shí),流速分布符合直管流動(dòng)特性,彎道效應(yīng)消失。張從聯(lián)等[9]對惠州抽水蓄能電站進(jìn)行模型試驗(yàn)研究,發(fā)現(xiàn)進(jìn)/出水口與水平彎段之間直線洞段的長度需要大于40倍引水洞洞徑才能滿足各流道均勻分流的要求。實(shí)際上,水平彎段會(huì)使進(jìn)/出水口在出流工況下的流道流量分配極為不均、流量偏向一側(cè),同時(shí)考慮擴(kuò)散段造成水平和垂直方向上的流動(dòng)分離,極可能使流道內(nèi)部的流場更加紊亂,加劇水體質(zhì)點(diǎn)的脈動(dòng)。

    由于國內(nèi)缺乏對帶水平彎段的側(cè)式進(jìn)/出水口紊流特性研究,且尤為關(guān)注偏流出流情況下攔污柵段的流場結(jié)構(gòu)。本文以帶水平彎段的某抽水蓄能電站上庫側(cè)式進(jìn)/出水口為研究對象,利用大渦模擬方法對進(jìn)/出水口進(jìn)行模擬,分析水平及垂直方向的流速、雷諾切應(yīng)力分布和進(jìn)/出水口頂部回流區(qū)高度,并進(jìn)一步研究攔污柵段的流場結(jié)構(gòu);從紊動(dòng)能的角度分析進(jìn)/出水口水頭損失產(chǎn)生的源頭,研究偏流出流下各流道不同斷面的紊動(dòng)強(qiáng)度分布并與均勻來流進(jìn)行對比;對進(jìn)/出水口各流道不同斷面的流速測點(diǎn)做功率譜密度分析(PSD),研究流場渦旋結(jié)構(gòu)在時(shí)間尺度上的分布。

    1 研究方法

    1.1 幾何模型

    某抽水蓄能電站上水庫側(cè)式進(jìn)/出水口模型試驗(yàn)示意如圖1所示,采用三墩四流道布置,水平、垂向擴(kuò)散角分別為33.774°、6.12°,Uin為隧洞斷面的進(jìn)口流速。在引水洞段存在平面偏轉(zhuǎn)角,轉(zhuǎn)角為29.27°,洞徑為9 m,直線洞身段長約為洞徑的17倍。模型試驗(yàn)比尺采用1∶40,進(jìn)/出水口采用有機(jī)玻璃制作,4臺(tái)機(jī)滿發(fā)出流工況下流量為262 m3/s,各流道流速測點(diǎn)采用聲學(xué)多普勒三維點(diǎn)式流速儀量測。

    進(jìn)/出水口網(wǎng)格如圖2所示,采用混合網(wǎng)格計(jì)算,在擴(kuò)散段和攔污柵段局部加密,第一層邊界層為0.000 125 m,壁面處的y+(用于衡量近壁處湍流計(jì)算的精度是否足夠)約1.0~3.5;內(nèi)部水體為均勻的正方體網(wǎng)格,尺寸約0.004 m,邊界層網(wǎng)格和正方體網(wǎng)格之間采用六棱柱網(wǎng)格過渡,尺寸與正方體網(wǎng)格一致,總網(wǎng)格量為600萬。引水洞段采用結(jié)構(gòu)化網(wǎng)格,總網(wǎng)格量約130萬。

    采用大渦模擬的方法[10-12]對進(jìn)/出水口4臺(tái)機(jī)組滿發(fā)死水位出流工況進(jìn)行模擬,隧洞斷面采用流速進(jìn)口邊界,依據(jù)流量按平均流速給出,庫區(qū)為壓力出口邊界,依據(jù)庫水位按靜壓給出,庫區(qū)頂部和檢修門槽頂部為壓力進(jìn)口邊界。亞格子模型為Smagorinsky-Lilly[13-14],并以動(dòng)態(tài)計(jì)算[15]的方法調(diào)整Smagorinsky常數(shù),壓力耦合采用算子分裂的壓力隱式算法(PISO算法),時(shí)間離散為有界二階隱格式,空間離散采用二階迎風(fēng)格式,時(shí)間步長取0.005 s。

    1.2 量綱一化處理與監(jiān)測點(diǎn)布置

    坐標(biāo)、壓強(qiáng)、流速以及流場二階指標(biāo)都經(jīng)過量綱一化處理,從以往的研究[16-18]發(fā)現(xiàn)流速和長度分別以擴(kuò)散段進(jìn)口處時(shí)均流速(Ub)和寬度(H)進(jìn)行量綱一化處理,本文Ub和H分別為3.22 m/s和9 m。坐標(biāo)的量綱一化均以擴(kuò)散段進(jìn)口處的坐標(biāo)為基準(zhǔn),以z方向?yàn)槔M(jìn)/出水口各點(diǎn)的高程減去擴(kuò)散段進(jìn)口處底板的高程,再除以H,得到量綱一化的高程。

    為獲取水平、垂直方向的流速分布,沿水流方向取8個(gè)斷面,x/H分別為-1.33、-2.06、-3.00、-4.40、-5.06、-5.44、-5.60、-6.10,稱DM1—DM8。在每個(gè)流道中,沿8個(gè)斷面中垂線上取垂向監(jiān)測點(diǎn),每條垂線布置11個(gè)監(jiān)測點(diǎn),如圖3(a)所示,因此4個(gè)流道共32條監(jiān)測線,總計(jì)352個(gè)垂向監(jiān)測點(diǎn);同理,在z/H=0.5的位置處取水平剖面,對8個(gè)斷面沿水平方向取10個(gè)監(jiān)測點(diǎn),如圖3(b)所示,總計(jì)320個(gè)監(jiān)測點(diǎn)。最后,為了與模型試驗(yàn)流場進(jìn)行對比,在攔污柵斷面,即x/H=-5.6處,每個(gè)流道的左、中、右處分別布置5個(gè)監(jiān)測點(diǎn),4個(gè)流道總計(jì)60個(gè)監(jiān)測點(diǎn),用于驗(yàn)證數(shù)值模型結(jié)果的準(zhǔn)確性。以上共732個(gè)監(jiān)測點(diǎn),每個(gè)監(jiān)測點(diǎn)的采樣頻率為100 Hz。

    1.3 模型驗(yàn)證

    為了驗(yàn)證數(shù)值模型的有效性,對比了數(shù)值模型與物理模型的時(shí)均流場。表1為攔污柵斷面中垂線處的流速分布,數(shù)值模型和物理模型在攔污柵斷面時(shí)均流速分布較為接近,流道1—流道4的平均相對誤差分別為10.9%、11.6%、4.0%和10.6%。

    除了驗(yàn)證數(shù)值模型的時(shí)均流場外,還需要驗(yàn)證瞬時(shí)流場的準(zhǔn)確性。湍流的主要特征是不規(guī)則性,其流速、壓強(qiáng)的時(shí)間序列呈現(xiàn)不規(guī)則的振蕩運(yùn)動(dòng),屬于隨機(jī)過程,而在隨機(jī)過程中,概率密度包含隨機(jī)變量的全部統(tǒng)計(jì)信息[19]。因此,判斷數(shù)值模型和物理模型的瞬時(shí)場是否一致可以通過驗(yàn)證某點(diǎn)瞬時(shí)流速的概率密度實(shí)現(xiàn)。取每個(gè)流道中垂線上的第3個(gè)測點(diǎn)進(jìn)行概率密度驗(yàn)證,如圖4所示,橫坐標(biāo)u+為標(biāo)準(zhǔn)化后的瞬時(shí)流速,發(fā)現(xiàn)各測點(diǎn)的物理模型和數(shù)值模型的概率密度較為接近,且基本符合標(biāo)準(zhǔn)正態(tài)分布;數(shù)值模型和物理模型的概率密度在峰值和均值存在一定差異,其主要原因是試驗(yàn)采用的測量儀器可能干擾內(nèi)部流場,數(shù)值模型網(wǎng)格解析度不足造成脈動(dòng)流速與實(shí)際存在差異。數(shù)值模型各測點(diǎn)的偏態(tài)系數(shù)(Cs)為0.44~0.70,峰態(tài)系數(shù)(Ce)為2.70~3.34。綜上所述,本數(shù)值模型計(jì)算結(jié)果能夠反映模型試驗(yàn)的時(shí)均和瞬時(shí)流場。

    2 結(jié)果及分析

    2.1 偏流出流下的流場分布

    2.1.1 橫向紊流結(jié)構(gòu)

    根據(jù)進(jìn)/出水口的平面時(shí)均流速分布(圖5)可知,偏流出流下流量主要集中在流道3、流道4,各流道分流比分別為0.64、0.81、1.26和1.29,流量分配極為不均;而均勻來流時(shí)流道2、流道3流量較大,各流道分流比分別為0.90、1.09、1.12和0.87。各流道主流核心區(qū)均靠近分流墩內(nèi)壁,邊流道和中間流道分別在x/H為-1、-2時(shí)出現(xiàn)流動(dòng)分離。

    圖6為水平斷面各流道的時(shí)均流速分布,橫坐標(biāo)代表在平面上該點(diǎn)與壁面的量綱一距離(L為沿y方向不同截面處的流道寬度)。均勻來流的水平流場結(jié)構(gòu)與偏流來流下的結(jié)構(gòu)相似,因此僅展現(xiàn)偏流的結(jié)果。圖6中發(fā)現(xiàn),各流道DM1的主流集中、流速達(dá)到最大,到DM3流動(dòng)分離最明顯,主流偏向內(nèi)側(cè),在DM4之后流動(dòng)分離減弱,流速分布逐漸均勻。除偏流出流下流道3的流速在DM4—DM8呈雙峰分布以外,其他各流道斷面流速分布均呈單峰分布,僅在流速的數(shù)值上存在差異。

    雷諾切應(yīng)力是一個(gè)重要的紊流指標(biāo),代表了流體質(zhì)點(diǎn)的脈動(dòng)動(dòng)量通量,它將動(dòng)能從高流速區(qū)再分配到低流速區(qū)域,實(shí)現(xiàn)動(dòng)量交換[20],其公式如下:

    式中:ρ為水的密度;u′、v′和w′分別為x、y和z方向的脈動(dòng)流速。偏流出流和均勻出流下大部分流道的水平方向雷諾切應(yīng)力分布相似,只有偏流出流下流道3的雷諾切應(yīng)力分布存在差異,因此為避免重復(fù)僅展示偏流出流下流道3、流道4的雷諾切應(yīng)力分布。如圖7所示,DM1位于分流墩處,流道面積收縮,水流質(zhì)點(diǎn)受到擠壓,此時(shí)該斷面各點(diǎn)的雷諾切應(yīng)力幾乎為0,流層之間脈動(dòng)動(dòng)量交換極少;沿著水流方向,在擴(kuò)散角的影響下發(fā)生流動(dòng)分離,雷諾切應(yīng)力逐漸增大,在DM3達(dá)到最大;DM4之后的雷諾切應(yīng)力沿程逐漸降低。與其他流道雷諾切應(yīng)力呈單峰值相比,偏流出流下流道3的雷諾切應(yīng)力分布在DM4—DM8呈現(xiàn)出一正一負(fù)雙峰值分布,在文獻(xiàn)[12]中的有壓擴(kuò)散器也發(fā)現(xiàn)類似的現(xiàn)象。流道3的水平流場遠(yuǎn)比有壓擴(kuò)散器復(fù)雜,偏流對其流場結(jié)構(gòu)影響較大,此外還受到分流墩和垂直擴(kuò)散角的約束,造成該流道DM4之后兩側(cè)的流速高、中間低,從而在高流速和低流速區(qū)的相互剪切下導(dǎo)致流道3在y/L=0.2~0.8之間形成一正一負(fù)的雷諾切應(yīng)力分布。

    2.1.2 垂向紊流結(jié)構(gòu)

    圖8展示了流道3、流道4在偏流和均勻來流下的垂向時(shí)均流速分布,流道1、流道2的幾何結(jié)構(gòu)分別與流道3、流道4一致,計(jì)算結(jié)果接近不再重復(fù)展示。分流墩墩頭處水流流速達(dá)到最大,經(jīng)過垂向擴(kuò)散,在x/H=-1.1發(fā)生流動(dòng)分離,主流核心區(qū)靠近中下部,頂部為反向流速;沿著水流方向,主流區(qū)高度受回流的擠壓而逐漸變窄。

    進(jìn)/出水口頂板反向流速容易造成攔污柵處流態(tài)欠佳,值得重點(diǎn)關(guān)注。受垂向擴(kuò)散角造成流動(dòng)分離的影響,擴(kuò)散段頂部形成回流,通過回流系數(shù)(χ)定義,回流系數(shù)計(jì)算公式如下:

    式中:u(x,z)和w(x,z)分別為x和z方向的時(shí)均流速;α為垂向擴(kuò)散角。通過計(jì)算發(fā)現(xiàn)曲線以上的水體均為反方向流速,曲線以下的水體為正方向流速,該曲線可以視作流動(dòng)分離交界。如圖8所示,當(dāng)水流經(jīng)過垂向擴(kuò)散后發(fā)生流動(dòng)分離,各流道回流區(qū)的高度迅速增加,到達(dá)調(diào)整段后逐漸降低,兩邊流道在距離擴(kuò)散段末端約1倍進(jìn)口寬度的位置,即x/H=6.4,回流區(qū)高度降低為0,實(shí)現(xiàn)水流再附著。Tornblom等[18]通過對8.5°平面擴(kuò)散器進(jìn)行風(fēng)洞試驗(yàn),在距擴(kuò)散段末端約5倍進(jìn)口寬度流體重新附著,該附著現(xiàn)象與本文結(jié)果較為接近。兩中間流道的回流面積均遠(yuǎn)大于兩邊流道,幾乎占據(jù)擴(kuò)散段面積的1/3,因此位于流道2、流道3的攔污柵更容易受到反向流速的影響。

    圖9為各流道斷面的垂向流速分布,發(fā)現(xiàn)均勻來流和偏流出流下邊流道1、流道4的流速分布接近且中間流道2、流道3的流速分布接近。與邊流道相比,中間流道在頂部反向流速范圍大,受回流區(qū)擠壓造成中下部主流流速大、產(chǎn)生較大的逆壓梯度,各流道垂向流速分布與Tornblom等試驗(yàn)所測流速分布相似[18,21]。如圖10所示,均勻來流和偏離條件下的雷諾切應(yīng)力分布相似,因此僅展示偏流出流的垂向雷諾切應(yīng)力分布。在DM1處,水體質(zhì)點(diǎn)受分流墩的影響,各流層間的動(dòng)量交換較少,雷諾切應(yīng)力幾乎為0,此處與Tornblom等[18]試驗(yàn)中進(jìn)口的切應(yīng)力分布不同,其有壓擴(kuò)散器的上游為均勻槽道流銜接且紊流已充分發(fā)展,靠近兩側(cè)壁面的雷諾切應(yīng)力較大而中間較小。到DM3時(shí),流動(dòng)分離顯著,流層間相互混摻,在垂向上存在一正一負(fù)的峰值。DM3之后,其正向峰值逐漸降低直至為0而負(fù)向峰值逐漸減少,到擴(kuò)散段末端后僅存在單個(gè)峰值,切應(yīng)力分布特征與Tornblom試驗(yàn)結(jié)果基本一致。

    2.1.3 攔污柵處流場分布

    據(jù)現(xiàn)有文獻(xiàn)[22-23],攔污柵處的漩滾、不利流態(tài)是造成其發(fā)生破壞的主要誘因之一。實(shí)際上,沿水流方向攔污柵剛度大,而橫向的攔污柵腹板由于厚度薄更容易受到漩滾造成的疲勞破壞[6,24-25],因此攔污柵段的水流流態(tài)是重點(diǎn)關(guān)注的區(qū)域。如圖11所示,各流道在流動(dòng)分離交界處產(chǎn)生橫軸漩渦,靠近主流區(qū)域流速和流速梯度較大,流線嚴(yán)重彎曲。相比邊流道,中間流道的漩滾范圍更大,幾乎一半以上的攔污柵條處于橫軸漩滾中,進(jìn)一步說明中間流道容易產(chǎn)生疲勞破壞的風(fēng)險(xiǎn)。攔污柵橫斷面時(shí)均流速矢量圖如圖12所示,各流道在中部和底部均存在縱軸漩滾,前者由流動(dòng)分離造成,后者由主流與壁面的剪切產(chǎn)生。為避免該區(qū)域的三軸漩滾,應(yīng)嚴(yán)格控制垂向、橫向擴(kuò)散角,防止流動(dòng)分離帶來大面積的漩渦。

    2.2 偏流出流下的紊動(dòng)能和紊動(dòng)強(qiáng)度

    進(jìn)/出水口水頭損失關(guān)乎到水輪機(jī)的發(fā)電效益,是工程中關(guān)注的重點(diǎn)。前人[7]已研究各種幾何因素對進(jìn)/出水口水頭損失的影響,計(jì)算水頭損失的方法主要通過能量方程,而對進(jìn)/出水口內(nèi)部的水頭損失分布研究較少。紊動(dòng)能(K)代表脈動(dòng)流速的能量,紊動(dòng)能較大的區(qū)域紊流混摻劇烈,造成水流的能量損耗,反映水體質(zhì)點(diǎn)的水頭損失[6,18,26-27],其公式如下:

    垂向紊動(dòng)能分布如圖13(a)所示,紊動(dòng)能主要集中在流動(dòng)分離交界,交界以上流速較小而交界以下流速較大造成能量的再分配,同時(shí)產(chǎn)生水頭損失。從平面紊動(dòng)能分布圖13(b)發(fā)現(xiàn),當(dāng)x/H<-1.5時(shí),紊動(dòng)能主要是水平和垂直方向的流動(dòng)分離造成;而當(dāng)-1.5≤x/H≤-0.75時(shí),在分流墩壁面附近發(fā)現(xiàn)紊動(dòng)能較大的區(qū)域,造成這一現(xiàn)象的主要因素可能有以下2點(diǎn):其一是分流墩壁面的強(qiáng)剪切作用,其二是分流墩附近的流動(dòng)分離。

    紊動(dòng)強(qiáng)度(Tu)是影響攔污柵結(jié)構(gòu)安全性的重要因素,較大的紊動(dòng)強(qiáng)度可能會(huì)破壞攔污柵條、主次梁焊縫等結(jié)構(gòu)薄弱位置[6],其公式如下:

    圖14為各流道不同斷面的紊動(dòng)強(qiáng)度垂向分布,紊動(dòng)強(qiáng)度較大的區(qū)域分布在中上部和底部的流動(dòng)分離交界,而位于主流核心區(qū)的紊動(dòng)強(qiáng)度較小。DM1、DM2在主流區(qū)的紊動(dòng)強(qiáng)度幾乎為0,DM3之后主流區(qū)紊動(dòng)強(qiáng)度緩慢增加但始終小于中上部和底部的紊動(dòng)強(qiáng)度,該結(jié)果與Hajaali等[21]在橫縱向雙向擴(kuò)散器的風(fēng)洞試驗(yàn)中紊動(dòng)強(qiáng)度垂向分布一致。因此,在攔污柵斷面,流動(dòng)分離交界不僅紊動(dòng)強(qiáng)度較大,而且該區(qū)域存在三軸漩滾,流態(tài)欠佳,從而更容易造成攔污柵的破壞。

    圖15為對8個(gè)斷面上測點(diǎn)的紊動(dòng)強(qiáng)度取平均,得到8個(gè)斷面的平均紊動(dòng)強(qiáng)度。結(jié)果表明,各流道平均紊動(dòng)強(qiáng)度沿?cái)嗝娣植枷嗨?,從DM1—DM8呈現(xiàn)出先劇增后緩慢降低的趨勢,且流道1、流道2略大于流道3、流道4的紊動(dòng)強(qiáng)度。與均勻來流相比,偏流出流下各流道在DM7(攔污柵斷面)的紊動(dòng)強(qiáng)度均較大,從流道1—流道4的紊動(dòng)強(qiáng)度分別提高11%、25%、29%、3%。

    2.3 偏流出流下的頻譜分析

    為了研究進(jìn)/出水口的脈動(dòng)動(dòng)能在垂向上時(shí)間尺度的分布,對8個(gè)斷面上流道中監(jiān)測點(diǎn)的時(shí)間序列做PSD。結(jié)果如圖16所示,橫坐標(biāo)為與底板相隔的垂直距離,縱坐標(biāo)為頻率。各斷面的上下兩側(cè)產(chǎn)生能量集中的低頻脈動(dòng),主頻約0.02 Hz左右,與文獻(xiàn)[28]進(jìn)/出水口水流頻率一致。DM1的主頻位于z/H=0.9和z/H=0.1處,也是流動(dòng)分離交界,遠(yuǎn)離主頻區(qū)域的能量逐漸降低;到DM3時(shí),主頻幅值達(dá)到峰值,脈動(dòng)動(dòng)能最大,此時(shí)受頂部反向流速的影響,主頻的位置發(fā)生變化,底部主頻上升至z/H=0.25,頂部主頻下降至z/H=0.78;DM4之后,主頻位置基本不發(fā)生變化,其幅值逐漸降低。

    本試驗(yàn)攔污柵斷面水流的主頻約0.02 Hz,因此由紊流作用在柵條上的力為低頻脈動(dòng)荷載,王光綸等[29]利用試驗(yàn)發(fā)現(xiàn)柵條固有頻率為113 Hz,申永康等[30]通過對攔污柵進(jìn)行模態(tài)分析得到整個(gè)柵架的一階固有頻率約14.98 Hz,均遠(yuǎn)大于紊流脈動(dòng)壓強(qiáng)的頻率,因此,紊流脈動(dòng)不會(huì)對攔污柵產(chǎn)生共振破壞。但是,鑒于主頻區(qū)域的能量較大、紊動(dòng)強(qiáng)度高,攔污柵與橫梁、次梁相交的焊縫區(qū)域仍存在靜疲勞破壞的風(fēng)險(xiǎn)。

    3 結(jié)? 論

    本文利用大渦模擬方法對帶水平彎段引水洞的抽水蓄能電站側(cè)式進(jìn)/出水口進(jìn)行三維數(shù)值模擬,深入研究了其內(nèi)部流場規(guī)律。主要結(jié)論如下:

    (1) 偏流出流下,進(jìn)出水口各流道水平流速分布極不均勻,主流偏向一側(cè),造成該側(cè)流道的流速顯著增加。在垂直方向上,主流核心區(qū)靠近過水?dāng)嗝嬷邢虏?,擴(kuò)散段垂向雷諾切應(yīng)力呈一正一負(fù)峰值分布,兩峰值分別源自流動(dòng)分離和底部壁面剪切作用。中間流道回流區(qū)高度大于兩邊流道,幾乎占據(jù)擴(kuò)散段的1/3,因此,中間流道攔污柵更容易受到反向流速影響。

    (2) 由水平、垂直擴(kuò)散角引起的流動(dòng)分離,是造成攔污柵處三軸漩滾的主要原因,且中間流道漩滾區(qū)域面積大于兩邊流道;同時(shí)在流動(dòng)分離交界處的紊動(dòng)強(qiáng)度較為集中,與均勻來流相比,偏流出流下攔污柵斷面的紊動(dòng)強(qiáng)度有明顯的增加。已有工程實(shí)例表明不利的水流流態(tài)和較高的紊動(dòng)強(qiáng)度是造成攔污柵損害的主要誘因之一,因此應(yīng)重點(diǎn)關(guān)注處于偏流出流下流動(dòng)分離交界處的攔污柵條。

    (3) 進(jìn)/出水口的分流墩壁面和流動(dòng)分離交界處存在較大的水頭損失;且其內(nèi)部的流動(dòng)分離造成靠近斷面中上部和底部2處出現(xiàn)比較集中的低頻脈動(dòng),主頻約0.02 Hz,不會(huì)對攔污柵造成共振破壞,沿水流方向脈動(dòng)的能量先增大后降低。

    參考文獻(xiàn):

    [1]孫雙科,柳海濤,李振中,等.抽水蓄能電站側(cè)式進(jìn)/出水口攔污柵斷面的流速分布研究[J].水利學(xué)報(bào),2007,38(11):1329-1335.(SUN S K,LIU H T,LI Z Z,et al.Study on velocity distribution behind the trashrack in lateral intake/outlet of pumped storage power station[J].Journal of Hydraulic Engineering,2007,38(11):1329-1335.(in Chinese))

    [2]王晨茜,張晨,張翰,等.側(cè)式進(jìn)/出水口流動(dòng)分離現(xiàn)象研究[J].水力發(fā)電學(xué)報(bào),2017,36(11):73-81.(WANG C X,ZHANG C,ZHANG H,et al.Flow separation in side inlets/outlets of pumped storage power stations[J].Journal of Hydroelectric Engineering,2017,36(11):73-81.(in Chinese))

    [3]高學(xué)平,朱洪濤,劉殷竹,等.進(jìn)/出水口雙向流動(dòng)結(jié)構(gòu)流速分布與脈動(dòng)規(guī)律研究[J].水利學(xué)報(bào),2022,53(6):722-732.(GAO X P,ZHU H T,LIU Y Z,et al.Research on velocity distribution and fluctuation law of bidirectional flow structure of inlet/outlet[J].Journal of Hydraulic Engineering,2022,53(6):722-732.(in Chinese))

    [4]朱洪濤,高學(xué)平,劉殷竹.進(jìn)/出水口出流擴(kuò)散脈動(dòng)流速規(guī)律試驗(yàn)研究[J].水力發(fā)電學(xué)報(bào),2022,41(7):129-139.(ZHU H T,GAO X P,LIU Y Z.Experimental study of velocity fluctuations in inlet/outlet under outflow diffusion[J].Journal of Hydroelectric Engineering,2022,41(7):129-139.(in Chinese))

    [5]YE F,GAO X P.Numerical simulations of the hydraulic characteristics of side inlet/outlets[J].Journal of Hydrodynamics:Ser B,2011,23(1):48-54.

    [6]高學(xué)平,劉永朋,孫博聞,等.側(cè)式進(jìn)/出水口擴(kuò)散出流紊動(dòng)強(qiáng)度變化規(guī)律研究[J].水利水電技術(shù),2020,51(3):91-100.(GAO X P,LIU Y P,SUN B W,et al.Study on variation law of turbulence intensity of diffused outflow at lateral inlet/outlet[J].Water Resources and Hydropower Engineering,2020,51(3):91-100.(in Chinese))

    [7]徐準(zhǔn),吳時(shí)強(qiáng).抽水蓄能電站側(cè)式進(jìn)/出水口隔墩布置對水力特性的影響[J].水利水電科技進(jìn)展,2020,40(3):21-27,67.(XU Z,WU S Q.Influence of arrangement of division piers in lateral inlet and outlet of pumped storage plants on hydraulic characteristics[J].Advances in Science and Technology of Water Resources,2020,40(3):21-27,67.(in Chinese))

    [8]張?zhí)m丁.響水澗抽水蓄能電站上、下庫進(jìn)(出)水口分流特性研究[J].水利水電科技進(jìn)展,2010,30(6):48-52.(ZHANG L D.Flow-dividing characteristics of intakes/outlets in upper and lower reservoirs of Xiangshuijian pumped storage station[J].Advances in Science and Technology of Water Resources,2010,30(6):48-52.(in Chinese))

    [9]張從聯(lián),朱紅華,鐘偉強(qiáng),等.惠州抽水蓄能電站下庫進(jìn)出水口水工模型試驗(yàn)研究[J].人民珠江,2003,24(6):10-14.(ZHANG C L,ZHU H H,ZHONG W Q,et al.Hydraulic modelling of inlet/outlet of lower reservoir for Huizhou pump storage plant[J].Pearl River,2003,24(6):10-14.(in Chinese))

    [10]丁全林,汪德爟,王玲玲.基于LBGK方法的粘性流體晃動(dòng)問題大渦模擬[J].水科學(xué)進(jìn)展,2009,20(5):652-657.(DING Q L,WANG D G,WANG L L.LES of viscous liquid sloshing problem based on the LBGK method[J].Advances in Water Science,2009,20(5):652-657.(in Chinese))

    [11]王文全,張立翔,郭亞昆,等.彎曲槽道邊壁振動(dòng)情況下湍流特性的大渦模擬[J].水科學(xué)進(jìn)展,2008,19(5):618-623.(WANG W Q,ZHANG L X,GUO Y K,et al.LES of turbulent flow in a vibrating curved-wall channel[J].Advances in Water Science,2008,19(5):618-623.(in Chinese))

    [12]徐振山,陳永平,張長寬.波流環(huán)境中垂向圓管射流三維運(yùn)動(dòng)和稀釋過程模擬[J].水科學(xué)進(jìn)展,2017,28(2):265-275.(XU Z S,CHEN Y P,ZHANG C K.Study on three-dimensional movement and dilution of a vertical round jet under the wave-current effect[J].Advances in Water Science,2017,28(2):265-275.(in Chinese))

    [13]SMAGORINSKY J.General circulation experiments with the primitive equations:I:the basic experiment[J].Monthly Weather Review,1963,91(3):99-164.

    [14]徐振山,陳永平,張長寬,等.波流環(huán)境中射流的大渦數(shù)值模擬[J].水科學(xué)進(jìn)展,2014,25(2):239-245.(XU Z S,CHEN Y P,ZHANG C K,et al.Large eddy simulation of a jet in wave-current environment[J].Advances in Water Science,2014,25(2):239-245.(in Chinese))

    [15]LILLY D K.A proposed modification of the Germano subgrid‐scale closure method[J].Physics of Fluids A:Fluid Dynamics,1992,4(3):633-635.

    [16]HERBST A H,SCHLATTER P,HENNINGSON D S.Simulations of turbulent flow in a plane asymmetric diffuser[J].Flow,Turbulence and Combustion,2007,79(3):275-306.

    [17]OBI S,NIKAIDO H,MASUDA S.Reynolds number effect on the turbulent separating flow in an asymmetric plane diffuser[C]∥Proceedings of the ASME/JSMR Fluids Engineering Division Summer Meeting.New York:American Society of Mechanical Engineers,1999.

    [18]TRNBLOM O,LINDGREN B,JOHANSSON A V.The separating flow in a plane asymmetric diffuser with 8.5° opening angle:mean flow and turbulence statistics,temporal behaviour and flow structures[J].Journal of Fluid Mechanics,2009,636:337-370.

    [19]張兆順,崔桂香,許春曉.湍流理論與模擬[M].北京:清華大學(xué)出版社,2005.(ZHANG Z S,CUI G X,XU C X.Theory and modeling of turbulence[M].Beijing:Tsinghua University Press,2005.(in Chinese))

    [20]王文娥,廖偉,漆力健.寬窄相間河道水流紊動(dòng)特性試驗(yàn)研究[J].水科學(xué)進(jìn)展,2020,31(3):394-403.(WANG W E,LIAO W,QI L J.Experiment of turbulent characteristics of flow in wide-and-narrow channels[J].Advances in Water Science,2020,31(3):394-403.(in Chinese))

    [21]HAJAALI A,STOESSER T.Flow separation dynamics in three-dimensional asymmetric diffusers[J].Flow,Turbulence and Combustion,2022,108(4):973-999.

    [22]皮仙槎.抽水蓄能電站攔污柵的設(shè)計(jì)[J].水力發(fā)電學(xué)報(bào),1990,9(4):86-93.(PI X C.Design of trash rack in pumped storage power station[J].Journal of Hydroelectric Engineering,1990,9(4):86-93.(in Chinese))

    [23]XUE B,MAO Y P,ZHANG C H,et al.Nonlinear fluid-induced vibro-impact analysis on the fatigue failure pattern of a large-scale trashrack with a reduced-order model[J].Structures,2023,49:467-478.

    [24]潘欽.抽水蓄能電站攔污柵流激振動(dòng)研究[D].大連:大連理工大學(xué),2007.(PAN Q.Study on flow-induced vibration of trash racks of pumped-storage hydroelectric plants[D].Dalian:Dalian University of Technology,2007.(in Chinese))

    [25]鐘耀.基于試驗(yàn)的攔污柵的激振分析及體型優(yōu)化研究[D].大連:大連理工大學(xué),2008.(ZHONG Y.Research on flow-induced vibration and body type optimization of trashracks based on test[D].Dalian:Dalian University of Technology,2008.(in Chinese))

    [26]陳啟剛,鐘強(qiáng).明渠紊流中渦結(jié)構(gòu)的運(yùn)動(dòng)規(guī)律[J].水科學(xué)進(jìn)展,2017,28(4):579-587.(CHEN Q G,ZHONG Q.Experimental study on the movement of vortices in turbulent open-channel flows[J].Advances in Water Science,2017,28(4):579-587.(in Chinese))

    [27]廖偉,張維樂,王文娥,等.臂坡對堰槽組合設(shè)施紊流結(jié)構(gòu)影響試驗(yàn)研究[J].水科學(xué)進(jìn)展,2023,34(3):465-479.(LIAO W,ZHANG W L,WANG W E,et al.Experimental study on the effect of arm slope on turbulent structure of weir-flume combination[J].Advances in Water Science,2023,34(3):465-479.(in Chinese))

    [28]高學(xué)平,袁野,劉殷竹,等.攔污柵結(jié)構(gòu)對進(jìn)出水口水力特性影響試驗(yàn)研究[J].水力發(fā)電學(xué)報(bào),2023,42(2):74-86.(GAO X P,YUAN Y,LIU Y Z,et al.Experimental study on effect of trash rack structure on hydraulic characteristics of inlet and outlet[J].Journal of Hydroelectric Engineering,2023,42(2):74-86.(in Chinese))

    [29]王光綸,張文翠,李未顯.抽水蓄能電站攔污柵結(jié)構(gòu)振動(dòng)特性模型試驗(yàn)研究[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,41(2):114-118.(WANG G L,ZHANG W C,LI W X.Experimental study of the free vibration characteristics of the trashrack of the Tianhuangping pumped storage station[J].Journal of Tsinghua University (Science and Technology),2001,41(2):114-118.(in Chinese))

    [30]申永康,方寒梅,趙春龍,等.大型攔污柵結(jié)構(gòu)液固耦合流激振動(dòng)分析[J].振動(dòng)與沖擊,2014,33(21):137-141.(SHEN Y K,F(xiàn)ANG H M,ZHAO C L,et al.Flow-induced vibration of large trash-racks considering fluid-structure interaction[J].Journal of Vibration and Shock,2014,33(21):137-141.(in Chinese))

    Study on large eddy simulation of lateral inlet/outlet under bending outflow

    Abstract:Due to the influence of topography and geology in actual engineering projects,water tunnels often have horizontal bends.These bends result in non-uniform inflow conditions,complicating the turbulence characteristics at the inlet and outlet compared to straight-line tunnels.This study employs Large Eddy Simulation to analyze a real engineering problem.The calculated results are compared with experimental data,revealing close agreement in flow velocity and probability density distribution.The findings show that under outflow conditions,the flow-split ratios for each channel are 0.64,0.81,1.26,and 1.29,indicating a highly non-uniform horizontal velocity distribution.In the vertical direction,the main flow occurred primarily in the lower middle section,and the vertical Reynolds shear stress in the diffusion segment exhibits a distinctive pattern with alternating positive and negative peaks.This phenomenon is primarily attributed to flow separation in the upper middle area and wall shearing near the bottom.The recirculation region′s height in the two middle channels was found to be greater than that in the side channels,indicating that trash racks in the central channels are more susceptible to reverse flow velocities.Additionally,flow separation near the trash rack gives rise to a three-axis vortex,while two energy-concentrated pulsations occur in the upper middle section and near the bottom.Moreover,under bending flow conditions,turbulence intensities in different channels increase by 11%,25%,29%,and 3%,respectively,compared to those in uniform flow situations.Consequently,the unfavorable flow regime and high turbulence intensity associated with bending flow pose a threat to the trash rack.

    Key words:lateral inlet/outlet;bending flow;Reynolds shear stress;flow separation;low frequency pulsation

    猜你喜歡
    出水口
    低水頭長距離輸水隧洞進(jìn)出水口水力特性研究
    人民長江(2024年13期)2024-07-06 16:37:47
    民機(jī)飲用水系統(tǒng)供水管路流動(dòng)性能計(jì)算研究
    科技視界(2020年19期)2020-07-30 14:04:34
    輸水隧洞坡角對側(cè)式進(jìn)/出水口水力特性影響研究
    沒有水龍頭的洗手池
    陜西檔案(2019年2期)2019-05-07 06:07:42
    潛流式人工濕地對生活污水生物需氧量和化學(xué)需氧量的影響分析
    發(fā)動(dòng)機(jī)出水口金屬墊片的密封設(shè)計(jì)分析
    新型出水口保護(hù)體在小農(nóng)水工程中的應(yīng)用
    山東水利(2018年6期)2018-03-24 13:00:35
    江蘇揚(yáng)州寶應(yīng)明代劉堡減水閘
    東南文化(2016年6期)2016-12-29 07:27:15
    側(cè)式進(jìn)/出水口工程布置進(jìn)流工況新進(jìn)展
    適用于農(nóng)田低壓輸水管道的玻璃鋼出水口
    国产男女内射视频| 草草在线视频免费看| 久久热精品热| 两个人的视频大全免费| 在线看a的网站| 99国产综合亚洲精品| 大话2 男鬼变身卡| 亚洲精品日韩在线中文字幕| 两个人的视频大全免费| 国产一区二区三区av在线| av福利片在线| 亚洲精品乱久久久久久| 综合色丁香网| av电影中文网址| 国产 精品1| 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 国产亚洲精品第一综合不卡 | 男人爽女人下面视频在线观看| 最近的中文字幕免费完整| 青春草国产在线视频| 精品99又大又爽又粗少妇毛片| 亚洲av成人精品一区久久| 国产在视频线精品| 美女cb高潮喷水在线观看| 精品一区二区三区视频在线| 久久人人爽人人片av| 亚洲图色成人| 欧美老熟妇乱子伦牲交| 久久久久国产网址| 男女免费视频国产| 交换朋友夫妻互换小说| 国产成人精品无人区| av网站免费在线观看视频| 亚洲av中文av极速乱| 哪个播放器可以免费观看大片| 亚洲精品色激情综合| 人人澡人人妻人| 亚洲四区av| 亚洲精品成人av观看孕妇| 亚洲国产精品成人久久小说| 国产69精品久久久久777片| 涩涩av久久男人的天堂| 国产不卡av网站在线观看| 一区二区日韩欧美中文字幕 | 嘟嘟电影网在线观看| 亚洲欧美色中文字幕在线| 80岁老熟妇乱子伦牲交| 亚洲久久久国产精品| 人妻夜夜爽99麻豆av| 久久久亚洲精品成人影院| 亚洲精品,欧美精品| 最黄视频免费看| 久久久国产一区二区| 国产av国产精品国产| 一本一本综合久久| 亚洲精品第二区| 如何舔出高潮| 亚洲精品久久成人aⅴ小说 | 18禁裸乳无遮挡动漫免费视频| 国产日韩一区二区三区精品不卡 | av在线老鸭窝| 久久亚洲国产成人精品v| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 国产精品国产三级国产av玫瑰| 亚洲欧洲日产国产| 免费观看性生交大片5| 日本欧美国产在线视频| 999精品在线视频| 久久久久国产网址| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 日韩强制内射视频| 中文字幕免费在线视频6| 少妇高潮的动态图| 国产伦精品一区二区三区视频9| 丝袜美足系列| 久久精品国产鲁丝片午夜精品| 夫妻性生交免费视频一级片| 精品久久国产蜜桃| 成人二区视频| 国产一区二区三区综合在线观看 | 在线观看免费高清a一片| 九九爱精品视频在线观看| 波野结衣二区三区在线| 亚洲内射少妇av| 亚洲人与动物交配视频| 国产一区二区三区综合在线观看 | 午夜91福利影院| 国产色爽女视频免费观看| 日本-黄色视频高清免费观看| 中文精品一卡2卡3卡4更新| 亚洲国产色片| 欧美人与善性xxx| 久久精品国产自在天天线| 少妇猛男粗大的猛烈进出视频| 乱人伦中国视频| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 国产 精品1| 日韩成人伦理影院| 国产熟女午夜一区二区三区 | 满18在线观看网站| 精品一区在线观看国产| 内地一区二区视频在线| 高清午夜精品一区二区三区| 青春草亚洲视频在线观看| 性色av一级| 日韩,欧美,国产一区二区三区| 99热国产这里只有精品6| 久久午夜福利片| 亚洲国产毛片av蜜桃av| 免费高清在线观看视频在线观看| 国模一区二区三区四区视频| 两个人免费观看高清视频| 国产精品偷伦视频观看了| 国产高清有码在线观看视频| 久久精品国产亚洲av天美| 人体艺术视频欧美日本| 中文字幕最新亚洲高清| 新久久久久国产一级毛片| 亚洲国产精品一区二区三区在线| 日韩成人av中文字幕在线观看| 久久久久久久国产电影| 边亲边吃奶的免费视频| 亚洲欧美日韩卡通动漫| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 久久99精品国语久久久| 国产在视频线精品| 久久久欧美国产精品| 免费观看a级毛片全部| 亚洲婷婷狠狠爱综合网| av视频免费观看在线观看| 黄片无遮挡物在线观看| 中文字幕精品免费在线观看视频 | 汤姆久久久久久久影院中文字幕| 卡戴珊不雅视频在线播放| 国产亚洲精品第一综合不卡 | 亚洲欧洲精品一区二区精品久久久 | 日韩熟女老妇一区二区性免费视频| 亚洲av.av天堂| 日韩亚洲欧美综合| 欧美激情极品国产一区二区三区 | 精品少妇内射三级| 日韩欧美一区视频在线观看| 免费播放大片免费观看视频在线观看| 午夜激情av网站| 高清视频免费观看一区二区| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 大香蕉久久网| 97在线视频观看| 亚洲av综合色区一区| 久久久国产精品麻豆| 亚洲国产欧美在线一区| 女的被弄到高潮叫床怎么办| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 久久精品国产鲁丝片午夜精品| 人成视频在线观看免费观看| 国产黄片视频在线免费观看| 久久狼人影院| 少妇被粗大猛烈的视频| 国产成人aa在线观看| 人妻制服诱惑在线中文字幕| 又粗又硬又长又爽又黄的视频| 日韩一区二区三区影片| 全区人妻精品视频| 亚洲精品久久午夜乱码| 两个人免费观看高清视频| 精品久久久久久久久亚洲| 人妻人人澡人人爽人人| 一个人看视频在线观看www免费| 亚洲伊人久久精品综合| 日本黄色片子视频| 亚洲av日韩在线播放| 夫妻午夜视频| 国产爽快片一区二区三区| 大话2 男鬼变身卡| 美女内射精品一级片tv| 秋霞伦理黄片| 美女国产视频在线观看| 婷婷色av中文字幕| 久久久亚洲精品成人影院| 高清不卡的av网站| 国产午夜精品一二区理论片| 亚洲一级一片aⅴ在线观看| 99热6这里只有精品| 女性生殖器流出的白浆| 日韩强制内射视频| 黄色欧美视频在线观看| 亚洲欧洲国产日韩| 久久99热6这里只有精品| 黑丝袜美女国产一区| 欧美少妇被猛烈插入视频| 自线自在国产av| 久久免费观看电影| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看| 男女国产视频网站| tube8黄色片| 91国产中文字幕| 亚洲欧美清纯卡通| 亚洲性久久影院| 97超视频在线观看视频| 午夜精品国产一区二区电影| 18在线观看网站| 国模一区二区三区四区视频| 看免费成人av毛片| 国产精品三级大全| 少妇精品久久久久久久| 欧美日韩在线观看h| 亚洲欧美清纯卡通| 伊人久久精品亚洲午夜| 男女国产视频网站| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 最近2019中文字幕mv第一页| 久久人妻熟女aⅴ| 国产精品女同一区二区软件| 成人午夜精彩视频在线观看| 五月玫瑰六月丁香| 99久久人妻综合| 这个男人来自地球电影免费观看 | freevideosex欧美| 婷婷色综合www| 亚洲久久久国产精品| 亚洲av福利一区| 狠狠精品人妻久久久久久综合| a级毛片在线看网站| 王馨瑶露胸无遮挡在线观看| 不卡视频在线观看欧美| 青春草亚洲视频在线观看| 一级爰片在线观看| 亚洲国产av新网站| 国产精品蜜桃在线观看| 国产极品天堂在线| 大香蕉97超碰在线| 免费高清在线观看日韩| 一区二区三区四区激情视频| 一区二区日韩欧美中文字幕 | 丰满乱子伦码专区| 亚洲av成人精品一区久久| 久久精品国产自在天天线| 久久精品夜色国产| 亚洲av欧美aⅴ国产| 18禁在线无遮挡免费观看视频| 亚洲精品一二三| 亚洲一级一片aⅴ在线观看| 飞空精品影院首页| 91aial.com中文字幕在线观看| 亚洲四区av| 日韩熟女老妇一区二区性免费视频| 国产亚洲最大av| 精品亚洲成国产av| 草草在线视频免费看| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 精品一区二区三区视频在线| 欧美xxⅹ黑人| av电影中文网址| 九草在线视频观看| 久久99热这里只频精品6学生| h视频一区二区三区| 亚洲综合色惰| 久久久久久久久久久丰满| 欧美亚洲 丝袜 人妻 在线| 热99国产精品久久久久久7| av有码第一页| 高清不卡的av网站| 欧美日本中文国产一区发布| 成人漫画全彩无遮挡| 国产色婷婷99| 亚洲精品乱久久久久久| 午夜福利视频在线观看免费| 日日撸夜夜添| 十八禁高潮呻吟视频| 午夜日本视频在线| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 大香蕉97超碰在线| 亚洲第一区二区三区不卡| 毛片一级片免费看久久久久| 亚洲国产日韩一区二区| 在线观看美女被高潮喷水网站| 两个人免费观看高清视频| 一本色道久久久久久精品综合| 特大巨黑吊av在线直播| 日本午夜av视频| 久久人人爽人人片av| 交换朋友夫妻互换小说| 国产精品久久久久成人av| 国产成人精品在线电影| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 日韩av免费高清视频| 亚洲少妇的诱惑av| 亚洲国产精品专区欧美| 久久久久久久国产电影| 又大又黄又爽视频免费| 久久久久久久国产电影| 99热国产这里只有精品6| 日韩在线高清观看一区二区三区| 色哟哟·www| 一级爰片在线观看| 亚洲国产欧美在线一区| 国产色婷婷99| 高清av免费在线| 国产一区有黄有色的免费视频| 国产伦精品一区二区三区视频9| 91精品国产九色| 大香蕉久久成人网| 国产精品国产三级国产专区5o| 人妻夜夜爽99麻豆av| 夜夜骑夜夜射夜夜干| 人人妻人人澡人人爽人人夜夜| 天美传媒精品一区二区| 国产精品99久久久久久久久| 综合色丁香网| 国产一区二区三区综合在线观看 | 超碰97精品在线观看| 亚洲精品久久久久久婷婷小说| 熟女人妻精品中文字幕| 亚洲精品视频女| 男人爽女人下面视频在线观看| 久久人妻熟女aⅴ| av电影中文网址| 久久久久久久久久人人人人人人| 丝袜美足系列| 99精国产麻豆久久婷婷| 婷婷色综合www| 亚洲成人av在线免费| 中文字幕av电影在线播放| 国产乱来视频区| 99国产精品免费福利视频| 在线观看国产h片| 国产又色又爽无遮挡免| 国产高清国产精品国产三级| 在线播放无遮挡| 国国产精品蜜臀av免费| 欧美变态另类bdsm刘玥| 欧美少妇被猛烈插入视频| 秋霞伦理黄片| 日韩强制内射视频| 精品人妻在线不人妻| 九色成人免费人妻av| 成人综合一区亚洲| 一级a做视频免费观看| 成人综合一区亚洲| 午夜久久久在线观看| 大又大粗又爽又黄少妇毛片口| 大码成人一级视频| 亚洲精品久久成人aⅴ小说 | 国产亚洲精品第一综合不卡 | 我要看黄色一级片免费的| 色网站视频免费| 成人无遮挡网站| 免费av不卡在线播放| 亚洲成人手机| 国产高清不卡午夜福利| 婷婷色综合www| av.在线天堂| 美女cb高潮喷水在线观看| av.在线天堂| 国内精品宾馆在线| 边亲边吃奶的免费视频| 亚洲精品第二区| 亚洲精品久久久久久婷婷小说| 韩国av在线不卡| 国产黄色免费在线视频| a级毛色黄片| 国产精品偷伦视频观看了| 国产亚洲一区二区精品| 国产熟女午夜一区二区三区 | 亚洲国产欧美在线一区| 日韩欧美一区视频在线观看| 青春草亚洲视频在线观看| 尾随美女入室| 日本爱情动作片www.在线观看| 日韩,欧美,国产一区二区三区| 中文字幕av电影在线播放| 五月开心婷婷网| 91精品国产国语对白视频| 国产在线免费精品| 性色avwww在线观看| 日本免费在线观看一区| 精品一区二区三区视频在线| 亚洲av国产av综合av卡| 亚洲精品中文字幕在线视频| 久久精品久久久久久噜噜老黄| 亚洲人成网站在线播| 日本wwww免费看| 国产在线一区二区三区精| av免费在线看不卡| 一级毛片aaaaaa免费看小| 黑人猛操日本美女一级片| 久久久欧美国产精品| 亚洲国产精品一区三区| 国产在线免费精品| 亚洲欧美清纯卡通| 国产精品 国内视频| 久久精品国产亚洲网站| 男人操女人黄网站| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 久久久国产欧美日韩av| 91aial.com中文字幕在线观看| 热99国产精品久久久久久7| 国产伦理片在线播放av一区| 亚洲国产精品专区欧美| 久久久精品区二区三区| 高清黄色对白视频在线免费看| videossex国产| 国产探花极品一区二区| 久久久国产一区二区| 欧美成人午夜免费资源| 久久精品国产自在天天线| 热re99久久精品国产66热6| 日韩人妻高清精品专区| 成年人午夜在线观看视频| 国产男女超爽视频在线观看| 亚洲精品久久午夜乱码| 亚洲无线观看免费| 精品一区在线观看国产| 日韩视频在线欧美| 考比视频在线观看| 精品久久蜜臀av无| 亚洲图色成人| 最新中文字幕久久久久| 亚洲av综合色区一区| 汤姆久久久久久久影院中文字幕| 老司机亚洲免费影院| 欧美日韩综合久久久久久| 中国国产av一级| 免费黄频网站在线观看国产| 久久久久久伊人网av| 免费黄网站久久成人精品| av福利片在线| 亚洲成人一二三区av| 中文天堂在线官网| 国产高清有码在线观看视频| 中国美白少妇内射xxxbb| 又粗又硬又长又爽又黄的视频| 一级爰片在线观看| 少妇被粗大猛烈的视频| 久久人人爽人人片av| 人妻夜夜爽99麻豆av| 免费观看无遮挡的男女| 精品一区二区三卡| 3wmmmm亚洲av在线观看| 中文字幕精品免费在线观看视频 | 久久国产精品大桥未久av| av免费观看日本| 少妇熟女欧美另类| 美女大奶头黄色视频| 高清欧美精品videossex| 国产亚洲欧美精品永久| h视频一区二区三区| 国产片内射在线| 免费看光身美女| 亚洲高清免费不卡视频| 国产女主播在线喷水免费视频网站| 韩国av在线不卡| 91国产中文字幕| 日韩一区二区三区影片| 久久免费观看电影| 黄色视频在线播放观看不卡| 91午夜精品亚洲一区二区三区| 亚洲精品中文字幕在线视频| 色吧在线观看| 两个人免费观看高清视频| 亚洲精品久久午夜乱码| 久久人妻熟女aⅴ| 一边摸一边做爽爽视频免费| 少妇丰满av| 国产日韩欧美在线精品| 人人澡人人妻人| 亚洲成人av在线免费| 水蜜桃什么品种好| 26uuu在线亚洲综合色| 一级爰片在线观看| 人人澡人人妻人| 色网站视频免费| 水蜜桃什么品种好| 日韩伦理黄色片| 国产成人精品一,二区| 国产日韩一区二区三区精品不卡 | 久久综合国产亚洲精品| 人体艺术视频欧美日本| 韩国av在线不卡| 精品国产一区二区久久| 丝袜喷水一区| 国产精品女同一区二区软件| 国产av码专区亚洲av| 亚洲欧美一区二区三区国产| 国产国语露脸激情在线看| 国产免费又黄又爽又色| 中文字幕最新亚洲高清| 久久国产精品大桥未久av| 欧美精品亚洲一区二区| 26uuu在线亚洲综合色| 亚洲三级黄色毛片| av在线播放精品| 三级国产精品欧美在线观看| 蜜桃国产av成人99| 大码成人一级视频| 不卡视频在线观看欧美| 亚洲精品日本国产第一区| 欧美精品一区二区免费开放| 丝袜喷水一区| 在线观看美女被高潮喷水网站| 高清午夜精品一区二区三区| 久久久精品区二区三区| 97超视频在线观看视频| 伦精品一区二区三区| 国产在视频线精品| 欧美精品一区二区免费开放| 欧美精品高潮呻吟av久久| 亚洲国产精品一区三区| 一级毛片黄色毛片免费观看视频| 性色av一级| 18禁裸乳无遮挡动漫免费视频| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区| 99国产精品免费福利视频| av在线播放精品| 亚洲在久久综合| 天堂8中文在线网| 国产精品99久久久久久久久| 亚洲,一卡二卡三卡| 亚洲中文av在线| 美女内射精品一级片tv| 欧美xxⅹ黑人| 亚洲精品第二区| 亚洲天堂av无毛| 97在线人人人人妻| 亚洲综合色网址| kizo精华| 卡戴珊不雅视频在线播放| 狠狠婷婷综合久久久久久88av| 在现免费观看毛片| 亚洲精品久久久久久婷婷小说| 晚上一个人看的免费电影| 色婷婷av一区二区三区视频| 精品一区二区三卡| 国产精品一区二区在线不卡| 国产一区二区三区综合在线观看 | 麻豆精品久久久久久蜜桃| 看非洲黑人一级黄片| 久久久国产欧美日韩av| 亚洲欧洲日产国产| 国产探花极品一区二区| 一区二区三区乱码不卡18| 男的添女的下面高潮视频| 亚洲av综合色区一区| 亚洲综合色惰| 精品亚洲成国产av| 亚洲欧洲精品一区二区精品久久久 | 日韩视频在线欧美| av天堂久久9| 大香蕉久久网| av不卡在线播放| 亚洲av国产av综合av卡| 亚洲精品国产av蜜桃| 日韩中文字幕视频在线看片| 蜜桃在线观看..| 亚洲情色 制服丝袜| 亚洲精品美女久久av网站| 午夜福利影视在线免费观看| 美女大奶头黄色视频| 韩国av在线不卡| 久久免费观看电影| 一本一本综合久久| 国产精品 国内视频| 少妇 在线观看| a级毛片在线看网站| 国产黄色视频一区二区在线观看| 成人黄色视频免费在线看| 中文字幕制服av| 爱豆传媒免费全集在线观看| 考比视频在线观看| 国产视频首页在线观看| 日韩,欧美,国产一区二区三区| 午夜老司机福利剧场| 欧美精品一区二区大全| 久久99热6这里只有精品| 精品国产一区二区久久| 国产淫语在线视频| 人人澡人人妻人| 99国产综合亚洲精品| 最近的中文字幕免费完整| 久久精品国产自在天天线| 国产精品麻豆人妻色哟哟久久| 国产高清国产精品国产三级| 国产免费现黄频在线看| 国产精品熟女久久久久浪| 久久99蜜桃精品久久| av在线app专区| 亚洲精品国产色婷婷电影| 看十八女毛片水多多多| 成年人午夜在线观看视频| 日本与韩国留学比较| 极品少妇高潮喷水抽搐| 亚洲欧洲精品一区二区精品久久久 | 亚洲av在线观看美女高潮| 亚洲av福利一区| 午夜福利,免费看| 在线观看美女被高潮喷水网站| 国产成人精品福利久久|