• <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)田低壓輸水管道的玻璃鋼出水口
    亚洲成人精品中文字幕电影 | 免费av中文字幕在线| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 欧美人与性动交α欧美精品济南到| 日韩av在线大香蕉| 三级毛片av免费| av网站在线播放免费| 国产一区二区三区综合在线观看| 午夜福利一区二区在线看| 老熟妇仑乱视频hdxx| 成人永久免费在线观看视频| 亚洲欧美日韩高清在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕一二三四区| 久久中文字幕人妻熟女| 国产人伦9x9x在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久国内视频| 免费看a级黄色片| 中文字幕精品免费在线观看视频| x7x7x7水蜜桃| 亚洲九九香蕉| 一a级毛片在线观看| 自线自在国产av| 久久人妻av系列| 中文欧美无线码| www.999成人在线观看| 男人舔女人的私密视频| 日韩免费高清中文字幕av| 午夜久久久在线观看| 高清毛片免费观看视频网站 | 国产精品爽爽va在线观看网站 | 久久香蕉国产精品| 亚洲少妇的诱惑av| 欧美精品啪啪一区二区三区| 熟女少妇亚洲综合色aaa.| 午夜免费激情av| 免费高清在线观看日韩| avwww免费| 热re99久久精品国产66热6| 岛国在线观看网站| 制服人妻中文乱码| 亚洲一卡2卡3卡4卡5卡精品中文| 国产免费现黄频在线看| 国产不卡一卡二| 一级a爱视频在线免费观看| 女警被强在线播放| 国产主播在线观看一区二区| 国产99白浆流出| 精品卡一卡二卡四卡免费| 91麻豆av在线| 制服诱惑二区| 久久人妻福利社区极品人妻图片| 国产av一区二区精品久久| 好男人电影高清在线观看| 岛国在线观看网站| 丰满迷人的少妇在线观看| 亚洲欧美激情综合另类| 国产精品 国内视频| 久久香蕉国产精品| 91成人精品电影| 天堂俺去俺来也www色官网| 91九色精品人成在线观看| 国产av一区在线观看免费| av福利片在线| 亚洲一码二码三码区别大吗| 国产成年人精品一区二区 | 亚洲人成电影免费在线| 国产成人精品久久二区二区91| 国产精品秋霞免费鲁丝片| 黄色视频,在线免费观看| 在线永久观看黄色视频| 纯流量卡能插随身wifi吗| 岛国在线观看网站| 免费一级毛片在线播放高清视频 | 窝窝影院91人妻| 国产精品亚洲av一区麻豆| 国产黄色免费在线视频| 久久久水蜜桃国产精品网| 男女下面进入的视频免费午夜 | 亚洲成av片中文字幕在线观看| 欧美av亚洲av综合av国产av| 日本黄色视频三级网站网址| 美国免费a级毛片| 精品一区二区三区av网在线观看| 十分钟在线观看高清视频www| 一区二区三区精品91| 少妇裸体淫交视频免费看高清 | 一边摸一边做爽爽视频免费| 亚洲欧美一区二区三区久久| 成人影院久久| 两人在一起打扑克的视频| 嫁个100分男人电影在线观看| 91av网站免费观看| 老汉色av国产亚洲站长工具| 国产在线观看jvid| 91字幕亚洲| 午夜两性在线视频| 又黄又爽又免费观看的视频| 免费在线观看日本一区| 夫妻午夜视频| av片东京热男人的天堂| 桃红色精品国产亚洲av| 久久人妻福利社区极品人妻图片| 国产午夜精品久久久久久| 在线观看一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲精品久久久久久毛片| 国产免费男女视频| 一级毛片女人18水好多| 午夜成年电影在线免费观看| 女人精品久久久久毛片| 男人操女人黄网站| 视频在线观看一区二区三区| 午夜亚洲福利在线播放| 琪琪午夜伦伦电影理论片6080| 欧美老熟妇乱子伦牲交| 日韩精品免费视频一区二区三区| 99久久人妻综合| 99久久久亚洲精品蜜臀av| 国产精品亚洲一级av第二区| 亚洲一区二区三区不卡视频| xxxhd国产人妻xxx| 国产在线精品亚洲第一网站| 校园春色视频在线观看| www国产在线视频色| 很黄的视频免费| 国产欧美日韩一区二区精品| 午夜福利在线观看吧| 天天添夜夜摸| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品自产拍在线观看55亚洲| 亚洲激情在线av| 久久精品国产亚洲av香蕉五月| 亚洲美女黄片视频| 不卡一级毛片| 国产亚洲av高清不卡| 黄色成人免费大全| 欧美日韩av久久| www日本在线高清视频| 国产亚洲精品综合一区在线观看 | 50天的宝宝边吃奶边哭怎么回事| 国产精品 欧美亚洲| 脱女人内裤的视频| 一级毛片精品| 女同久久另类99精品国产91| 亚洲成人精品中文字幕电影 | 国产成人av教育| 天堂影院成人在线观看| 国产亚洲精品一区二区www| 国产三级在线视频| 啦啦啦 在线观看视频| 午夜免费成人在线视频| 国产三级黄色录像| 国产国语露脸激情在线看| 人人妻人人添人人爽欧美一区卜| 欧美日韩黄片免| 麻豆久久精品国产亚洲av | 日韩大尺度精品在线看网址 | 亚洲欧美日韩无卡精品| 黄频高清免费视频| 亚洲人成伊人成综合网2020| 久久精品亚洲av国产电影网| 91国产中文字幕| 久久午夜亚洲精品久久| 丝袜人妻中文字幕| 18美女黄网站色大片免费观看| 十八禁人妻一区二区| 超色免费av| av网站在线播放免费| 成人18禁在线播放| 久久国产精品影院| 新久久久久国产一级毛片| 国产野战对白在线观看| 成人亚洲精品一区在线观看| 99久久国产精品久久久| 日日摸夜夜添夜夜添小说| 99久久久亚洲精品蜜臀av| 亚洲第一青青草原| 久久久久久久久久久久大奶| 成熟少妇高潮喷水视频| 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| www.999成人在线观看| 欧美日韩乱码在线| 中文字幕人妻熟女乱码| 一级毛片精品| 久久久久久免费高清国产稀缺| 1024香蕉在线观看| 国产成人精品久久二区二区91| 一进一出抽搐动态| 亚洲av成人一区二区三| 老汉色∧v一级毛片| 久久精品国产亚洲av香蕉五月| 99国产精品一区二区蜜桃av| 日本vs欧美在线观看视频| 国产激情欧美一区二区| 国产精品亚洲一级av第二区| av视频免费观看在线观看| 欧美日韩视频精品一区| 黑人欧美特级aaaaaa片| 一区二区三区激情视频| a级片在线免费高清观看视频| 久久久久久久精品吃奶| 中出人妻视频一区二区| 国产精品一区二区精品视频观看| 久久精品aⅴ一区二区三区四区| 黄片大片在线免费观看| 亚洲在线自拍视频| 色老头精品视频在线观看| 少妇被粗大的猛进出69影院| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| 日韩中文字幕欧美一区二区| 88av欧美| 亚洲人成电影免费在线| 欧美乱色亚洲激情| 男女下面进入的视频免费午夜 | svipshipincom国产片| 亚洲中文字幕日韩| 日韩三级视频一区二区三区| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3 | av天堂在线播放| 日韩视频一区二区在线观看| 久99久视频精品免费| 国产高清国产精品国产三级| 亚洲五月色婷婷综合| 亚洲片人在线观看| 高清毛片免费观看视频网站 | 自线自在国产av| 日韩欧美免费精品| 国产精品一区二区在线不卡| 国产精品综合久久久久久久免费 | 最好的美女福利视频网| 香蕉丝袜av| 十分钟在线观看高清视频www| 久久 成人 亚洲| 99国产极品粉嫩在线观看| 久久中文看片网| 在线观看免费日韩欧美大片| 亚洲精品国产色婷婷电影| 无限看片的www在线观看| 亚洲国产精品合色在线| 少妇粗大呻吟视频| 色在线成人网| 亚洲激情在线av| 99国产综合亚洲精品| 18美女黄网站色大片免费观看| 久久这里只有精品19| 亚洲精品国产色婷婷电影| 欧美人与性动交α欧美软件| 神马国产精品三级电影在线观看 | 悠悠久久av| 一a级毛片在线观看| 黄色视频不卡| 亚洲一区二区三区欧美精品| 免费不卡黄色视频| 精品久久久久久久久久免费视频 | 亚洲国产精品999在线| 免费高清视频大片| 大型av网站在线播放| 人妻久久中文字幕网| 又紧又爽又黄一区二区| 男女午夜视频在线观看| 免费在线观看影片大全网站| 动漫黄色视频在线观看| 国产精品国产av在线观看| 一边摸一边抽搐一进一小说| 亚洲 欧美一区二区三区| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 亚洲人成伊人成综合网2020| 国产97色在线日韩免费| 久久久久精品国产欧美久久久| 亚洲一区二区三区不卡视频| 老汉色∧v一级毛片| 精品人妻在线不人妻| 国产精品爽爽va在线观看网站 | 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久成人av| 欧美一区二区精品小视频在线| 国产精品国产av在线观看| 国产伦人伦偷精品视频| 亚洲成国产人片在线观看| 国产成人精品久久二区二区免费| 精品一区二区三区视频在线观看免费 | 国产亚洲精品久久久久久毛片| 深夜精品福利| 免费搜索国产男女视频| 黄色片一级片一级黄色片| 成熟少妇高潮喷水视频| 免费一级毛片在线播放高清视频 | 国产伦一二天堂av在线观看| 精品一区二区三区视频在线观看免费 | 成人av一区二区三区在线看| 欧美av亚洲av综合av国产av| 国产亚洲欧美在线一区二区| 成年女人毛片免费观看观看9| 亚洲av成人一区二区三| 国产一区在线观看成人免费| 国产亚洲精品久久久久5区| 亚洲一区高清亚洲精品| 黑人操中国人逼视频| 又紧又爽又黄一区二区| 亚洲成人精品中文字幕电影 | 国产成人精品久久二区二区91| 女警被强在线播放| 91大片在线观看| 757午夜福利合集在线观看| 97碰自拍视频| 天堂动漫精品| 亚洲av熟女| 人人妻人人澡人人看| 两个人看的免费小视频| 岛国在线观看网站| 亚洲精品在线观看二区| 国产精品一区二区三区四区久久 | 久久精品国产99精品国产亚洲性色 | 午夜a级毛片| 99国产综合亚洲精品| 久久久国产成人免费| svipshipincom国产片| 国产成+人综合+亚洲专区| 三上悠亚av全集在线观看| 宅男免费午夜| 人成视频在线观看免费观看| 97超级碰碰碰精品色视频在线观看| 淫妇啪啪啪对白视频| 亚洲国产精品合色在线| 国产av在哪里看| 在线天堂中文资源库| 国产黄色免费在线视频| 精品日产1卡2卡| 色播在线永久视频| 悠悠久久av| 欧美激情 高清一区二区三区| 免费人成视频x8x8入口观看| x7x7x7水蜜桃| 日日摸夜夜添夜夜添小说| 久久久久久免费高清国产稀缺| 精品一品国产午夜福利视频| 中文字幕人妻丝袜制服| x7x7x7水蜜桃| 日日摸夜夜添夜夜添小说| 国产99白浆流出| 国产91精品成人一区二区三区| 免费搜索国产男女视频| 久久久久国产一级毛片高清牌| 精品国产国语对白av| 999久久久国产精品视频| 亚洲精品美女久久久久99蜜臀| 日韩高清综合在线| 国产精品 欧美亚洲| 亚洲午夜精品一区,二区,三区| 精品久久久久久,| 90打野战视频偷拍视频| 亚洲人成电影免费在线| 91在线观看av| 一区二区日韩欧美中文字幕| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| 免费日韩欧美在线观看| av天堂在线播放| 国产精品爽爽va在线观看网站 | 超碰成人久久| 91麻豆精品激情在线观看国产 | 人妻丰满熟妇av一区二区三区| 一区二区三区激情视频| 欧美午夜高清在线| a级毛片在线看网站| 亚洲男人的天堂狠狠| 国产精品国产高清国产av| 久久国产乱子伦精品免费另类| 亚洲一区中文字幕在线| 亚洲精品国产精品久久久不卡| www.精华液| 久久人人精品亚洲av| 欧美日韩一级在线毛片| 亚洲熟女毛片儿| 狠狠狠狠99中文字幕| 国产精品野战在线观看 | 最近最新中文字幕大全电影3 | 欧美中文日本在线观看视频| 999精品在线视频| 日韩高清综合在线| 曰老女人黄片| 极品教师在线免费播放| 亚洲 欧美一区二区三区| 久久中文看片网| 午夜精品在线福利| 国产成人av激情在线播放| 欧美日韩黄片免| 欧美激情极品国产一区二区三区| 午夜精品在线福利| 啦啦啦在线免费观看视频4| 天天躁夜夜躁狠狠躁躁| 国产精品爽爽va在线观看网站 | 在线观看免费高清a一片| 神马国产精品三级电影在线观看 | 国产精品爽爽va在线观看网站 | 在线国产一区二区在线| 婷婷精品国产亚洲av在线| 国产成人精品在线电影| e午夜精品久久久久久久| 国产免费男女视频| 一a级毛片在线观看| 视频在线观看一区二区三区| 最新美女视频免费是黄的| 亚洲一区二区三区色噜噜 | 日韩高清综合在线| 国产成人av激情在线播放| 久久久久久人人人人人| 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到| 免费观看精品视频网站| 日本免费一区二区三区高清不卡 | 欧美黑人欧美精品刺激| 99精品久久久久人妻精品| 色播在线永久视频| 在线视频色国产色| 一级片免费观看大全| 级片在线观看| 亚洲成人久久性| 亚洲性夜色夜夜综合| 欧美成人免费av一区二区三区| 9热在线视频观看99| 1024香蕉在线观看| av福利片在线| 亚洲专区字幕在线| 婷婷六月久久综合丁香| 久久久久精品国产欧美久久久| 国产乱人伦免费视频| 免费在线观看完整版高清| av福利片在线| 欧美日韩国产mv在线观看视频| 国产又爽黄色视频| 一级片免费观看大全| 国产aⅴ精品一区二区三区波| 12—13女人毛片做爰片一| 淫秽高清视频在线观看| 大码成人一级视频| 久久香蕉精品热| 视频区欧美日本亚洲| 亚洲欧美精品综合久久99| 亚洲熟女毛片儿| videosex国产| 老司机午夜福利在线观看视频| 麻豆成人av在线观看| 欧美激情极品国产一区二区三区| 女人精品久久久久毛片| 亚洲精品一卡2卡三卡4卡5卡| 深夜精品福利| 国产麻豆69| 国产精品国产高清国产av| 男女床上黄色一级片免费看| 无人区码免费观看不卡| 免费在线观看视频国产中文字幕亚洲| 国产蜜桃级精品一区二区三区| tocl精华| 99国产精品一区二区三区| 久久天堂一区二区三区四区| 国产av又大| 国产亚洲精品一区二区www| 黄色女人牲交| 别揉我奶头~嗯~啊~动态视频| 男男h啪啪无遮挡| 成人三级黄色视频| 久久久久国内视频| 一级毛片精品| 搡老岳熟女国产| 男人操女人黄网站| 国产成+人综合+亚洲专区| 两个人免费观看高清视频| 欧美精品一区二区免费开放| 脱女人内裤的视频| 欧美av亚洲av综合av国产av| 又大又爽又粗| 国产熟女xx| 欧美激情 高清一区二区三区| 99久久精品国产亚洲精品| 一级a爱视频在线免费观看| 欧美激情 高清一区二区三区| 两人在一起打扑克的视频| 欧美日韩av久久| 欧美日韩中文字幕国产精品一区二区三区 | 欧美成人午夜精品| 久久久国产精品麻豆| 人妻久久中文字幕网| 动漫黄色视频在线观看| 91av网站免费观看| 国产精品一区二区在线不卡| 亚洲欧美激情综合另类| 亚洲全国av大片| 又大又爽又粗| 一本综合久久免费| 中文字幕人妻丝袜一区二区| 精品欧美一区二区三区在线| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 日韩欧美在线二视频| 婷婷精品国产亚洲av在线| 黄片大片在线免费观看| 国产免费男女视频| 99在线人妻在线中文字幕| 精品欧美一区二区三区在线| 成人av一区二区三区在线看| 高清黄色对白视频在线免费看| 久久 成人 亚洲| 1024香蕉在线观看| 可以免费在线观看a视频的电影网站| 精品国产超薄肉色丝袜足j| 亚洲av第一区精品v没综合| 级片在线观看| 久久人人精品亚洲av| 欧美乱妇无乱码| 色老头精品视频在线观看| 国产三级黄色录像| 国产人伦9x9x在线观看| 亚洲国产精品999在线| 久99久视频精品免费| 亚洲精品成人av观看孕妇| 91av网站免费观看| 怎么达到女性高潮| 80岁老熟妇乱子伦牲交| 亚洲一区高清亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 午夜福利在线观看吧| 亚洲欧美一区二区三区久久| 99热只有精品国产| 亚洲欧美精品综合久久99| 97碰自拍视频| 真人一进一出gif抽搐免费| 午夜激情av网站| 少妇裸体淫交视频免费看高清 | 91国产中文字幕| 国产乱人伦免费视频| 免费女性裸体啪啪无遮挡网站| 99国产精品99久久久久| 国产精品二区激情视频| 国产精品 欧美亚洲| 色在线成人网| 无人区码免费观看不卡| 不卡一级毛片| 国产精品 国内视频| 欧美国产精品va在线观看不卡| 国产精品爽爽va在线观看网站 | 亚洲中文字幕日韩| 久久久久久久久久久久大奶| 亚洲自偷自拍图片 自拍| 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 女人被躁到高潮嗷嗷叫费观| 男人的好看免费观看在线视频 | 99国产精品99久久久久| 91精品三级在线观看| 中文字幕最新亚洲高清| 一区福利在线观看| 精品欧美一区二区三区在线| 婷婷精品国产亚洲av在线| 99热国产这里只有精品6| 999精品在线视频| 国产精品偷伦视频观看了| av免费在线观看网站| 高清黄色对白视频在线免费看| 久久精品亚洲精品国产色婷小说| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品亚洲av| 亚洲情色 制服丝袜| 一区二区日韩欧美中文字幕| 一级片'在线观看视频| 成人三级黄色视频| 久久久国产成人精品二区 | 久久久国产成人免费| 在线观看www视频免费| 黄片大片在线免费观看| 国产视频一区二区在线看| 久久久久久久久久久久大奶| 久久久久国产精品人妻aⅴ院| 日韩欧美免费精品| 国产伦一二天堂av在线观看| 亚洲性夜色夜夜综合| 99久久99久久久精品蜜桃| 露出奶头的视频| 国产成人免费无遮挡视频| 88av欧美| 高潮久久久久久久久久久不卡| 国产成年人精品一区二区 | 欧美黑人精品巨大| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 精品午夜福利视频在线观看一区| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 久久欧美精品欧美久久欧美| 久久这里只有精品19| 精品人妻1区二区| 女人精品久久久久毛片| 亚洲美女黄片视频| 国产单亲对白刺激| 亚洲美女黄片视频| 久久草成人影院| 国产不卡一卡二| 多毛熟女@视频| 黄色女人牲交| 最近最新中文字幕大全电影3 | 国产黄色免费在线视频| 国产亚洲av高清不卡| 成人手机av| 欧美一区二区精品小视频在线| 国产国语露脸激情在线看| 国产精品影院久久| 久久久国产成人精品二区 |