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

    偏心漸縮管和彎管對(duì)泵站進(jìn)水管道流態(tài)的影響研究

    2025-02-28 00:00:00李曉超高志鍇謝敏萍李曉鶴羅毅唐華林林曉恒李君
    關(guān)鍵詞:數(shù)值模擬泵站

    摘要:為深入研究不同偏心漸縮管長(zhǎng)度和彎管角度對(duì)開(kāi)敞式泵站進(jìn)水管道水力特性的影響,探索進(jìn)水管道水力性能參數(shù)隨彎管角度和偏心漸縮管長(zhǎng)度變化的規(guī)律,以黃河下游引黃提灌泵站為例,基于Fluent仿真軟件,對(duì)14種不同彎管角度和偏心漸縮管長(zhǎng)度方案進(jìn)行了數(shù)值模擬,并通過(guò)評(píng)價(jià)函數(shù)對(duì)典型斷面的進(jìn)水管道流場(chǎng)進(jìn)行了流態(tài)分析。對(duì)比分析方案1—5發(fā)現(xiàn):偏心漸縮管長(zhǎng)度的改變對(duì)進(jìn)水管道水力損失及出口斷面處的水流速度分布均勻度和速度加權(quán)平均角的影響較大,且偏心漸縮管長(zhǎng)度分別為600、700、800 mm時(shí)進(jìn)水管道水流的整體水力特性較好。對(duì)方案6—14的計(jì)算分析發(fā)現(xiàn):彎管角度的變化主要對(duì)彎管段水流流態(tài)產(chǎn)生影響,對(duì)其他部位的水流擾動(dòng)不大。進(jìn)一步結(jié)合評(píng)價(jià)函數(shù)圖,得出結(jié)論:①?gòu)澒芙嵌仍酱?,水流過(guò)渡越平穩(wěn),水力損失越小;彎管角度越小,則水力損失越大。②偏心漸縮管長(zhǎng)度越長(zhǎng),出口斷面水流速度分布均勻度和速度加權(quán)平均角的值越大,整流效果越好,但水力損失也越大;偏心漸縮管長(zhǎng)度越短,整流效果越差。

    關(guān)鍵詞:泵站;彎管;偏心漸縮管;進(jìn)水流態(tài);數(shù)值模擬

    中圖分類號(hào):TV765文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):2096-6792(2025)01-0103-12

    泵站進(jìn)水管道是連接泵站進(jìn)水池和水泵之間的一段過(guò)渡性流道,其水力流態(tài)對(duì)水泵性能產(chǎn)生直接影響,應(yīng)盡量減少進(jìn)水流道附件及其長(zhǎng)度,且布置平順,保證嚴(yán)密不漏氣,以便使水流以最小的水力損失均勻且平順地流入葉輪。一般來(lái)講,泵站進(jìn)水管道包括直管段、彎管段和漸變管段。其中,漸變管即管徑逐漸發(fā)生改變的管道,由于工程中常受其他條件限制,為提高管道水流水力性能,經(jīng)常將漸變管用作進(jìn)水管道穿墻、進(jìn)水管道與水泵進(jìn)口等管徑變化處的連接過(guò)渡,漸變管根據(jù)管道中心是否變化可分為偏心漸變管和同心漸變管,又可分為漸縮管和漸擴(kuò)管兩種。由于進(jìn)水管道中的彎管和漸變管會(huì)改變水流方向,對(duì)水流影響較大,故彎管和漸變管水力性能的優(yōu)劣對(duì)水流流態(tài)具有重要影響。

    在泵站工程目前的研究當(dāng)中,由于計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)[1-3]技術(shù)的迅猛發(fā)展以及CFD軟件計(jì)算分析方便快捷等特點(diǎn),CFD軟件已被廣泛應(yīng)用于泵站工程的模擬研究中[4-7]。如:張建偉等[8]通過(guò)CFD軟件,對(duì)泵站進(jìn)水管道岔管的不同接入角度進(jìn)行了流態(tài)和模態(tài)的計(jì)算與分析,發(fā)現(xiàn)岔管接入角度越小,水流對(duì)管道的沖擊越小。陳義春等[9]通過(guò)CFD相關(guān)軟件,針對(duì)泵站機(jī)組進(jìn)水管段經(jīng)常抽不上水和進(jìn)水管水流含氣量高的問(wèn)題,設(shè)計(jì)并計(jì)算分析了原始及優(yōu)化方案的機(jī)組工作狀況,發(fā)現(xiàn)優(yōu)化后的機(jī)組進(jìn)水管道水力性能有了顯著提升。高傳昌等[10]針對(duì)泵站進(jìn)水管路布置引發(fā)的流態(tài)惡化、泵裝置效率低下等問(wèn)題,借助Fluent軟件對(duì)優(yōu)化前后的進(jìn)水管路布設(shè)方案進(jìn)行了數(shù)值模擬,發(fā)現(xiàn)優(yōu)化方案效果顯著。董亮等[11]對(duì)比分析了不同湍流模型下90°彎管的流場(chǎng)分布,發(fā)現(xiàn)標(biāo)準(zhǔn)k-ε模型、Realizable k-ε模型、標(biāo)準(zhǔn)k-ω模型在管道直線段的流態(tài)模擬效果良好,RNG k-ε模型除了在彎管段θ=60°時(shí)的模擬存在偏差外,其余狀態(tài)下的模擬效果良好,SST k-ω模型則除了在彎管段θ=30°時(shí)的模擬存在偏差外,其余狀態(tài)下的模擬效果良好,LES(Large Eddy Simulation)模型則在θ=0°~90°下的流動(dòng)模擬效果均較好。夏學(xué)坤等[12]借助CFD軟件,對(duì)不同偏心異徑管連接進(jìn)水管道的3種優(yōu)化方案進(jìn)行了計(jì)算分析,發(fā)現(xiàn)3種方案下的管路水力性能均能滿足實(shí)際要求。李曉超等[13]對(duì)優(yōu)化后的泵站進(jìn)水管路進(jìn)行了數(shù)值模擬,給出了進(jìn)水管徑推薦值750 mm。樊新建等[14]探索了泵站正向進(jìn)水池?cái)U(kuò)散角分別為20°、25°、30°、35°、40°、46°時(shí)對(duì)泵站進(jìn)水池流場(chǎng)結(jié)構(gòu)的影響,發(fā)現(xiàn):大擴(kuò)散角對(duì)水流約束小,產(chǎn)生水流脫壁現(xiàn)象,回流擠壓主流,流態(tài)變差;而小擴(kuò)散角對(duì)水流約束強(qiáng),但水流不易擴(kuò)散,也使得水泵進(jìn)水條件惡化。黃先北等[15]研究了懸空高對(duì)進(jìn)水流道水力特性的影響,發(fā)現(xiàn)吸氣渦的強(qiáng)度基本隨著懸空高的增加而減小。郁片紅等[16]通過(guò)Fluent軟件對(duì)進(jìn)水流道存在的主流集中、回流及偏流等問(wèn)題進(jìn)行了研究分析。施偉等[17]針對(duì)豎井貫流式泵裝置的進(jìn)出水流道,進(jìn)行了性能影響的關(guān)鍵參數(shù)標(biāo)準(zhǔn)化研究。車曉紅等[18]對(duì)泵站肘形進(jìn)水流道和虹吸式出水流道的上升角和下降角進(jìn)行了設(shè)計(jì)分析。吳鵬飛等[19]通過(guò)試驗(yàn)與數(shù)值模擬相結(jié)合,分析了泵站進(jìn)水池的表面吸氣渦的機(jī)理并提出了相應(yīng)解決措施。張德勝等[20]通過(guò)數(shù)值模擬與PIV(Particle Image Velocimetry)試驗(yàn)結(jié)果對(duì)比,探討了泵站進(jìn)水池吸入口的漩渦結(jié)構(gòu)及其發(fā)展過(guò)程和相關(guān)特性??椎轮t等[21]分析探討了水平式進(jìn)水管的懸空高對(duì)泵站進(jìn)水池進(jìn)水特性的影響。鄭云浩等[22]通過(guò)極差分析法對(duì)模擬結(jié)果進(jìn)行分析,探討了簸箕形進(jìn)水流道的型線優(yōu)化。高傳昌等[23]通過(guò)數(shù)值模擬和試驗(yàn)分析了泵站進(jìn)水池超低水位下不同消渦組合方案的整流效果。初長(zhǎng)虹等[24]分析了5種不同進(jìn)水口淹沒(méi)深度對(duì)流場(chǎng)漩渦的影響。高學(xué)平等[25]對(duì)抽蓄電站進(jìn)/出水流道不同反坡坡比的水力性能進(jìn)行了分析探討。張晨等[26]通過(guò)SPH(Smoothed Particle Hydrodynamics)方法對(duì)豎井式進(jìn)/出水流道進(jìn)行模擬,并通過(guò)試驗(yàn)驗(yàn)證了其可行性。

    總體而言,目前專家學(xué)者關(guān)于泵站的研究分析大多集中于進(jìn)水前池和出水流道,關(guān)于進(jìn)水管道的研究較少,鮮有進(jìn)水管道中彎管和漸變管水流流態(tài)的研究。但進(jìn)水管道出口直接與水泵進(jìn)口相連,其水力性能又對(duì)水泵產(chǎn)生重要影響,故十分有必要開(kāi)展彎管和漸變管對(duì)泵站進(jìn)水管道水流流態(tài)的影響研究。

    為深入探討彎管和漸變管對(duì)泵站進(jìn)水管道的水力特性影響,結(jié)合引黃提灌工程實(shí)際,共設(shè)置了14種不同彎管角度和偏心漸縮管長(zhǎng)度的計(jì)算方案,通過(guò)數(shù)值模擬研究不同方案下的進(jìn)水管道水力性能,以期進(jìn)一步減少管道水力損失、探索管路布設(shè)最優(yōu)方式、提高泵站效益。

    1工程概況

    該引黃提灌泵站位于黃河下游,主要保障灌區(qū)當(dāng)?shù)氐霓r(nóng)業(yè)用水、工業(yè)用水和生活用水等,在當(dāng)?shù)厮W(wǎng)結(jié)構(gòu)中具有重要地位,為當(dāng)?shù)厣钌a(chǎn)和社會(huì)經(jīng)濟(jì)發(fā)展作出了重要貢獻(xiàn)。該泵站進(jìn)水池為開(kāi)敞式的,共裝設(shè)24臺(tái)離心泵,2臺(tái)機(jī)組共用一進(jìn)水池,4臺(tái)機(jī)組共用一引水池,進(jìn)水池水深3.4 m,布置示意圖如圖1(a)所示。其中進(jìn)水管采用垂直式布置,進(jìn)水管出口與水泵進(jìn)口相連,進(jìn)水管道包括喇叭管、進(jìn)口直管、彎管、出口直管、偏心漸縮管5部分,各部分長(zhǎng)度等詳細(xì)參數(shù)如圖1(b)所示。為深入分析進(jìn)水管道內(nèi)偏心漸縮管長(zhǎng)度和彎管角度變化對(duì)進(jìn)水管道水流水力特性的影響,結(jié)合計(jì)算資源,僅對(duì)1臺(tái)機(jī)組在運(yùn)行工況下的進(jìn)水管道水力特性進(jìn)行研究。

    2數(shù)值計(jì)算

    2.1計(jì)算方案根據(jù)泵站設(shè)計(jì)標(biāo)準(zhǔn)[27-29]可知,進(jìn)水管道流速宜在1.5~2.0 m/s范圍內(nèi)。針對(duì)進(jìn)水管道而言,喇叭管主要起整流作用,使得水流盡可能均勻、平順地流入進(jìn)水管道;直管段內(nèi)水流基本保持原有流動(dòng)慣性,沿管道方向流動(dòng),在實(shí)際工程中,為防止進(jìn)水管道頂部積氣,水平段會(huì)設(shè)計(jì)成傾斜向上;彎管段改變了水流流向,影響管道水流流態(tài);偏心漸縮管一方面引導(dǎo)水流以較小的水力損失流入水泵,一方面對(duì)管道水流進(jìn)行整流,使得流出進(jìn)水管道的水流更加均勻、平順地流入水泵,其水流流態(tài)直接決定進(jìn)入水泵的水流流態(tài),進(jìn)而影響水泵裝置的效率。因此,為探索偏心漸縮管長(zhǎng)度和彎管角度變化對(duì)進(jìn)水管道水力特性的影響,筆者針對(duì)偏心漸縮管不同長(zhǎng)度、彎管不同角度分別設(shè)置對(duì)應(yīng)方案,其中彎管角度為其兩端面所成夾角,如圖2所示。進(jìn)水管直徑、懸空高、后壁距等其他參數(shù)保持不變,各方案在設(shè)計(jì)流量0.88 m3/s下進(jìn)行數(shù)值模擬,參數(shù)設(shè)置詳見(jiàn)表1。

    2.2模型與網(wǎng)格劃分進(jìn)水流道計(jì)算模型包括前池、進(jìn)水池和進(jìn)水管,通過(guò)NX.UG軟件對(duì)計(jì)算模型按1∶1比例進(jìn)行建模,采用大型流體力學(xué)計(jì)算軟件ANSYS Fluent自帶的Meshing模塊對(duì)模型進(jìn)行網(wǎng)格劃分,剖分網(wǎng)格采用混合網(wǎng)格格式(內(nèi)部為六面體結(jié)構(gòu)網(wǎng)格,邊壁采用多面體結(jié)構(gòu)網(wǎng)格),為滿足工程要求及計(jì)算精度,對(duì)流速梯度變化大的喇叭口、彎管、偏心漸縮管、邊壁處進(jìn)行局部加密。以方案4為例,經(jīng)網(wǎng)格無(wú)關(guān)性分析(圖3(a)),確定網(wǎng)格數(shù)量為430萬(wàn)個(gè)。模型如圖3(b)所示。

    2.3控制方程與邊界條件數(shù)值模擬計(jì)算采用連續(xù)方程、雷諾時(shí)均N-S方程[30]和兩方程模型中能更好模擬進(jìn)水池表面渦及附壁渦的Realizable k-ε模型[31-34],采用SIMPLEC算法進(jìn)行求解。計(jì)算域的進(jìn)口斷面為前池距進(jìn)水流道足夠遠(yuǎn)處,該處可認(rèn)為流速分布均勻,為速度進(jìn)口。流道出口斷面即計(jì)算域出口斷面,設(shè)置為自由出流。由于水面穩(wěn)定,將其設(shè)為對(duì)稱面條件,其他固體邊壁均為壁面,采用壁面無(wú)滑移條件。

    2.4分析截面與評(píng)價(jià)函數(shù)為更好地分析進(jìn)水管道內(nèi)部流態(tài),探索相關(guān)變化規(guī)律,現(xiàn)選取如下特征截面進(jìn)行分析:在Z方向上選擇進(jìn)水流道中心剖面Z0;在X方向上選擇彎管出口斷面X0、偏心漸縮管進(jìn)口斷面X1和偏心漸縮管出口斷面X2;在Y方向上選擇喇叭口斷面Y0和彎管進(jìn)口斷面Y1。共計(jì)6個(gè)截面進(jìn)行流線、流速分析,截面示意圖如圖4所示。

    為更好地評(píng)判各方案下進(jìn)水流道的水力性能,以進(jìn)水管道水力損失、進(jìn)水管道出口斷面的水流速度分布均勻度和速度加權(quán)平均角為評(píng)價(jià)函數(shù)[35], 水流速度分布均勻度和速度加權(quán)平均角分別按式(1)和式(2)計(jì)算。

    Vu=1-1ua∑(uai-ua)2m×100%,(1)

    θ=∑uai90°-arctanutiuai∑uai。(2)

    式中:Vu為水流速度分布均勻度;uai和ua分別為進(jìn)水管道出口斷面各單元的軸向流速和平均軸向流速;m為進(jìn)水流道出口斷面的單元總數(shù),所謂單元即網(wǎng)格劃分的單元網(wǎng)格,通過(guò)該式利用各單元流速得出該面的Vu;θ為水流速度加權(quán)平均角;uai和uti分別為進(jìn)水管道出口斷面各單元的軸向流速和橫向流速。

    流速分布均勻度和速度加權(quán)平均角兩者越接近理想值,認(rèn)為水流性能越好,其中流速分布均勻度和速度加權(quán)平均角的理想值分別為100%和90°。

    水力損失hw是指運(yùn)動(dòng)過(guò)程中單位質(zhì)量液體的機(jī)械能損失,計(jì)算公式如下:

    hw=E2-E1=P2ρg-P1ρg+Z2-Z1+V222g-V122g。(3)

    式中:hw為水力損失;E2、E1分別為進(jìn)水管道進(jìn)口、出口的總能量;P2、P1分別為進(jìn)水管道進(jìn)口斷面和出口斷面處的平均動(dòng)水壓強(qiáng),Pa;V2、V1分別為進(jìn)水管道進(jìn)口斷面和出口斷面處的平均流速,m/s;ρ為流體密度,kg/m3;g為重力加速度,10 m/s2。

    3結(jié)果分析

    3.1相同彎管角度、不同偏心漸縮管長(zhǎng)度方案1—5的截面分析結(jié)果為90°彎管下、偏心漸縮管的長(zhǎng)度分別為400、500、600、700、800 mm時(shí)的運(yùn)算結(jié)果。

    圖5為進(jìn)水管道中心剖面的流速云圖。由圖5可知:在各方案下,水流經(jīng)喇叭管整流后比較平穩(wěn)地流入直管段;彎管內(nèi)側(cè)水流流速大于外側(cè)水流流速,符合流體力學(xué)規(guī)律(半徑越大,壓力越大,流速越低;反之亦反);水流從彎管流出后,保持原有運(yùn)動(dòng)慣性,向出口直管段的偏上部分流動(dòng),故直管段下部產(chǎn)生低流速區(qū),由于直管段不改變水流方向,低流速區(qū)呈增大趨勢(shì),直至出口偏心漸縮管處;水流受偏心漸縮管整流,流速增大,且偏心漸縮管的長(zhǎng)度影響明顯,漸縮管越長(zhǎng),整流效果越好,出口水流流速越來(lái)越接近真實(shí)流速,速度分布也越來(lái)越均勻。

    圖6為喇叭管進(jìn)口剖面的流速云圖。由圖6可知:各方案下喇叭管口處的流速分布均呈環(huán)狀,較為均勻,流速值大小也較接近,流線分布類似,整體無(wú)明顯差別,偏心漸縮管的長(zhǎng)度變化對(duì)喇叭管處水流影響不大。

    圖7為彎管進(jìn)口剖面圖。圖中各方案均出現(xiàn)高流速區(qū)和低流速區(qū),高流速區(qū)靠近管道內(nèi)側(cè),呈類圓狀分布,且隨漸縮管長(zhǎng)度的增加,彎管局部的高流速區(qū)域略有增大趨勢(shì)(偏心漸縮管越長(zhǎng),水流過(guò)渡越平緩,軸向速度損失越?。?外側(cè)則為低流速區(qū),呈月牙狀分布。各方案的整體流速分布類似,管道中心處均有漩渦產(chǎn)生。

    圖8為彎管出口剖面流速云圖。觀察圖8各方案的流線分布可知:各方案下均存在兩個(gè)位置對(duì)稱漩渦,分布在彎管左右兩側(cè),且右側(cè)漩渦強(qiáng)度大于左側(cè)漩渦強(qiáng)度,兩個(gè)漩渦強(qiáng)度有隨著偏心漸縮管長(zhǎng)度的增加而不斷減弱的趨勢(shì),這是因?yàn)槠臐u縮管的增長(zhǎng)使得流速降低,漩渦強(qiáng)度亦隨之減弱;各方案下的流線分布大致相似,除方案1外,其他方案的流速分布較類似;方案1中的低流速區(qū)略大,隨著偏心漸縮管長(zhǎng)度的增加而增加;彎管出口斷面的低流速區(qū)不斷減小,流速分布愈加均勻。

    圖9為偏心漸縮管進(jìn)口剖面流速云圖。由圖9可知:偏心漸縮管處管道斷面偏心收縮,水流從管道下部低流速區(qū)向管道上部出口流動(dòng),各方案下流線分布一致;各方案下水流流速分布相似,高流速區(qū)位于管道上側(cè),低流速區(qū)在下側(cè),隨著漸縮管長(zhǎng)度的增長(zhǎng),漸縮管進(jìn)口斷面的低流速區(qū)不斷變小,斷面水流流速分布趨均勻。

    圖10為偏心漸縮管出口剖面圖。由圖10可知:各方案下水流繼續(xù)保持向上的流動(dòng)趨勢(shì),流線分布一致;速度分布大致相同,近似成環(huán)狀分布;但高流速區(qū)面積隨著偏心漸縮管長(zhǎng)度的增加不斷變大,慢慢鋪滿整個(gè)斷面,這是因?yàn)槠臐u縮管越長(zhǎng),越利于管道水流的調(diào)整,使得流速整體分布更加均勻。

    結(jié)合評(píng)價(jià)函數(shù)的計(jì)算公式,方案1—5即偏心漸縮管長(zhǎng)度分別為400、500、600、700、800 mm的評(píng)價(jià)函數(shù)結(jié)果分布如圖11所示。由圖11可知:

    1)隨著偏心漸縮管長(zhǎng)度的增加,沿程損失增大,進(jìn)水管道水力損失變大,同時(shí)偏心漸縮管長(zhǎng)度變長(zhǎng),水流整流區(qū)域長(zhǎng)度越長(zhǎng),由于整流區(qū)域內(nèi)的流體流態(tài)較平直管段流體流態(tài)紊亂,故其整流區(qū)域內(nèi)部流動(dòng)損失有所增加,使得水力損失進(jìn)一步增大,但方案整體的水力損失差距不大,最大水力損失差距為0.35 cm。

    2)速度分布均勻度均在90%以上,且隨著偏心漸縮管長(zhǎng)度的增加不斷增大,但增長(zhǎng)幅度漸緩。其中:方案1速度分布均勻度最小,為90.132%;方案5的最大,為94.126%??梢?jiàn),偏心漸縮管的變化對(duì)出口斷面流速分布均勻度影響較大。

    3)各方案出口斷面的速度加權(quán)平均角在75°以上,隨著偏心漸縮管長(zhǎng)度的增加,其出口斷面速度加權(quán)平均角明顯增大,但其增長(zhǎng)速率逐漸減緩。方案1出口斷面的速度加權(quán)平均角最小,為75.998°;方案5的最大,為81.618°??梢?jiàn),偏心漸縮管的變化對(duì)出口斷面速度加權(quán)平均角影響較大。

    綜上所述,偏心漸縮管長(zhǎng)度的變化對(duì)進(jìn)水管道水流流態(tài)影響較大。偏心漸縮管越長(zhǎng),速度分布均勻度、出口斷面速度加權(quán)平均角越大,即整流效果越好,但同時(shí)水力損失在不斷增大。方案3、4、5下的綜合水流流態(tài)更好,水流漩渦強(qiáng)度較低,且水流速度分布均勻度、速度加權(quán)平均角較大,水力損失差距不是很大。因此,重點(diǎn)進(jìn)行偏心漸縮管長(zhǎng)度分別為600、700、800 mm時(shí),不同角度彎管下進(jìn)水管道的水流流態(tài)研究。

    3.2相同偏心漸縮管長(zhǎng)度、不同彎管角度偏心漸縮管長(zhǎng)度分別為600、700、800 mm時(shí),不同彎管角度方案下的流場(chǎng)分布類似,限于文章篇幅,以700 mm為例,分析彎管不同角度對(duì)進(jìn)水流道水力特性的影響。圖12—17中,從左至右方案4、方案9—11的截面分析結(jié)果分別為彎管角度取90°、60°、45°、35°時(shí)的運(yùn)算結(jié)果。

    圖12為流道中心斷面圖。觀察圖12可知:各方案下直管段內(nèi)水流流態(tài)差別不大;彎管處管道內(nèi)、外側(cè)形成高、低流速區(qū),且彎管角度越小,外側(cè)的低流速區(qū)面積越大、內(nèi)側(cè)的高流速區(qū)面積越小;直管段下部出現(xiàn)低流速區(qū),一直延伸至偏心漸縮管處;僅彎管段內(nèi)水流流態(tài)差別較大,其他區(qū)域的流速分布整體一致。

    圖13為喇叭管進(jìn)口剖面圖。由圖13可知:各方案喇叭口處的流速分布均呈環(huán)狀,且斷面內(nèi)大部分流速均為1.34 m/s左右;水流從四周向管道軸線方向流動(dòng),各方案的流線分布類似,故彎管角度的變化對(duì)喇叭管進(jìn)口處水流流態(tài)無(wú)明顯影響。

    圖14為彎管進(jìn)口剖面圖。由圖14可知:各方案斷面速度均大致呈月牙狀分布,管道內(nèi)側(cè)為高流速區(qū),外側(cè)為低流速區(qū);隨著彎管角度的減小,高流速區(qū)變小,低流速區(qū)變大,流速分層更加明顯,速度分布更不均勻。這是因?yàn)殡S著彎管角度的不斷減小,彎管邊壁過(guò)渡連接處越來(lái)越不平穩(wěn),邊壁產(chǎn)生漩渦,水流流態(tài)變差,流速分布越來(lái)越不均勻;觀察流線可知,該斷面中出現(xiàn)明顯漩渦,隨彎管角度的不斷增加,漩渦不斷從管道邊壁向彎管中心移動(dòng),強(qiáng)度不斷增加。這是因?yàn)?0°大角度彎管下,水流在直管、彎管連接的過(guò)渡段較為平穩(wěn),邊壁不易產(chǎn)生漩渦,僅管道中心水流受流速梯度影響,有漩渦產(chǎn)生。當(dāng)彎管角度減小后,過(guò)渡段不再平穩(wěn),水流流動(dòng)產(chǎn)生水力損失,局部損失增大,管道邊壁連接處水流流態(tài)變差,引發(fā)漩渦,但彎管角度減小,水流流速減小,漩渦強(qiáng)度變?nèi)酢?/p>

    圖15為彎管出口剖面圖。觀察可知:各方案下彎管出口斷面的流速分布與彎管進(jìn)口斷面分布類似,管道角度越小,流速區(qū)分層越明顯,低流速區(qū)位于管道外側(cè),高流速區(qū)位于管道內(nèi)側(cè);彎管出口斷面出現(xiàn)兩個(gè)漩渦,且隨著彎管角度的不斷減小,逐漸由管道邊壁向管道中下部靠攏,強(qiáng)度有所減弱;漩渦亦是由于過(guò)渡水流引起,彎管角度大的水力損失小,但較高的水流流速致使水流速度環(huán)量較大,從而引發(fā)較強(qiáng)的漩渦;彎管角度的降低,使得水流過(guò)渡更加不平穩(wěn),使水流流態(tài)進(jìn)一步紊亂,但因其水流流速的降低使漩渦強(qiáng)度不斷減弱。

    圖16為偏心漸縮管進(jìn)口剖面圖。由圖16可知:各方案偏心漸縮管進(jìn)口斷面流速分布相同,斷面下半部分為低流速區(qū),高流速區(qū)位于管道斷面上部;隨著偏心漸縮管的收縮作用,水流從管道下部低流速區(qū)向管道上部出口流動(dòng),流線分布整體一致,且隨著彎管角度的增加,漸縮管進(jìn)口斷面的流線分布更加平順均勻。

    圖17為偏心漸縮管出口剖面圖。由圖17可知:各方案偏心漸縮管出口斷面的流速分布相同,呈圓環(huán)狀,且管道內(nèi)大部分為高流速區(qū),流速值均無(wú)較大差別;流線分布相同,依然保持原有運(yùn)動(dòng)趨勢(shì),向管道上部流動(dòng),無(wú)漩渦產(chǎn)生。

    圖18為不同偏心漸縮管方案下的評(píng)價(jià)函數(shù)結(jié)果分布圖。由圖18可知:①不同偏心漸縮管長(zhǎng)度方案下,彎管角度的點(diǎn)線圖分布規(guī)律基本一致,出口斷面的速度分布均勻度在90°和60°彎管方案下較大,在45°彎管方案下劇烈下降,在30°彎管方案下又有回升,整體趨勢(shì)是隨著彎管角度的減小,先降低后增加。這是因?yàn)樗俣确植季鶆蚨戎饕茌S面速度和平均軸面速度影響,彎管角度較大情況下,其水力損失小,速度分布均勻度就高,但角度較小時(shí),整體損失增加,水流流速降低,再經(jīng)偏心漸縮管整流,其水流軸向流速及其平均軸向流速差值有所減小,出口斷面速度分布均勻度明顯提高。②速度加權(quán)平均角在90°和60°彎管方案下較大,但60°彎管方案下速度加權(quán)平均角比90°彎管下的大,后隨彎管角度減小,在45°和30°彎管下其值較小,整體變化趨勢(shì)是隨著彎管角度的減小,先略有增加而后減小。這是因?yàn)殡S著彎管角度的增加,整體水流過(guò)渡更加順暢,其水力損失小,整體的速度分布更趨均勻,但在較大的彎管角度下,由于其流速較高,漩渦強(qiáng)度較大,其水流產(chǎn)生的速度環(huán)量也大,則水流的軸向速度分量損失偏大,使得速度加權(quán)平均角反有略降。

    速度分布均勻度的最大差值是0.7%,速度加權(quán)平均角的最大差值為0.3°,可見(jiàn)彎管角度的變化對(duì)出口斷面的速度分布均勻度和速度加權(quán)平均角影響不顯著;水力損失隨著彎管角度的減小而增大,且增大趨勢(shì)越發(fā)明顯,彎管角度為30°時(shí)水力損失均為最大,最大水力損失差約為0.6 cm。對(duì)比分析可知,彎管角度變化對(duì)管道損失的影響較顯著。綜上:彎管角度為60°和90°時(shí),進(jìn)水流道的水流流態(tài)較好。

    圖19為不同彎管角度方案下的評(píng)價(jià)函數(shù)結(jié)果分布圖。由圖19可知:不同彎管角度下的評(píng)價(jià)函數(shù)折線分布類似,各方案下的偏心漸縮管長(zhǎng)度越大,相應(yīng)的速度分布均勻度和速度加權(quán)平均角越大,其管道水力損失也增大。這是因?yàn)橄嗤瑥澒芙嵌认?,偏心漸縮管長(zhǎng)度越長(zhǎng),沿程損失越大。偏心漸縮管內(nèi)紊亂的水流流態(tài)會(huì)使得流動(dòng)損失進(jìn)一步增加,但其整流效果隨著偏心漸縮管長(zhǎng)度的增加會(huì)越來(lái)越好。結(jié)合上述分析,推薦偏心漸縮管的長(zhǎng)度為700 mm。

    4結(jié)論與展望

    1)偏心漸縮管長(zhǎng)度的變化對(duì)進(jìn)水管道前面的過(guò)流部件內(nèi)水流流態(tài)的影響不大,其長(zhǎng)度變化主要對(duì)出口斷面水流的速度分布均勻度和加權(quán)平均角有較大影響,而水力損失的變化值差距不大。

    2)偏心漸縮管長(zhǎng)度越長(zhǎng),整流效果越好,出口斷面水流分布越均勻,但管道水力損失也相應(yīng)增大,其推薦值以700 mm為優(yōu)。

    3)彎管角度的變化主要影響彎管段內(nèi)的水流流態(tài),其角度越小,水力損失越大,低流速區(qū)面積也越大,流速不斷降低,漩渦強(qiáng)度也越小,推薦大彎管角度60°~90°為宜。

    4)彎管角度的變化對(duì)管道水力損失的影響較為顯著,出口流速分布均勻度及加權(quán)平均角雖也有變化,但波動(dòng)值較小。

    本文僅研究了特定條件下不同彎管角度和偏心漸縮管長(zhǎng)度下的進(jìn)水管道水流流態(tài)變化規(guī)律,進(jìn)水管直徑、懸空高、后壁距等關(guān)鍵參數(shù)改變后的進(jìn)水管道水流流態(tài)變化規(guī)律有待進(jìn)一步研究。

    參考文獻(xiàn)

    [1]HIRSCH C.Numerical computation of internal and external flows:fundamentals of computational fluid dynamics[M].England:Butterworth-Heinemann,2007.

    [2]VERSTEEG H K.An introduction to computational fluid dynamics the finite volume method[M].2nd edition.England:Prentice Hall,2007.

    [3]UDAYKUMAR H S.The finite volume method in computaional fluid dynamics:an advanced introduction with OpenFOAM and Matlab[M].Germany:Springer International Publishing,2016.

    [4]羅興锜,吳大轉(zhuǎn).泵技術(shù)進(jìn)展與發(fā)展趨勢(shì)[J].水力發(fā)電學(xué)報(bào),2020,39(6):1-17.

    [5]羅興锜,朱國(guó)俊,馮建軍.水輪機(jī)技術(shù)進(jìn)展與發(fā)展趨勢(shì)[J].水力發(fā)電學(xué)報(bào), 2020,39(8):1-18.

    [6]羅興锜,戈振國(guó),朱國(guó)俊,等.沖擊式水輪機(jī)技術(shù)進(jìn)展與發(fā)展趨勢(shì)[J].水力發(fā)電學(xué)報(bào),2023,42(2):116-134.

    [7]高傳昌,解克宇,黃丹,等.不同水位對(duì)泵站進(jìn)水池流態(tài)影響的數(shù)值模擬[J].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,36(6):10-14.

    [8]張建偉,華薇薇,王立彬.基于流態(tài)-模態(tài)聯(lián)合分析的泵站管道結(jié)構(gòu)優(yōu)化[J].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,40(6):89-96.

    [9]陳義春,張慧,戴盛,等.小型泵站吸水井進(jìn)水管位置數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2017,35(11):962-967.

    [10]高傳昌,高余鑫,孫龍?jiān)?基于標(biāo)準(zhǔn)k-ε紊流模型的泵站進(jìn)水管路水力特性研究[J].灌溉排水學(xué)報(bào),2021,40(5):106-113.

    [11]董亮,劉厚林,代翠,等.不同湍流模型在90°彎管數(shù)值模擬中的應(yīng)用[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,40(12):18-22.

    [12]夏學(xué)坤,楊祖強(qiáng),羅希.CFD方法在泵站進(jìn)水管更新改造設(shè)計(jì)中的運(yùn)用[J].水電與新能源,2014(8):25-27,34.

    [13]李曉超,謝敏萍,李君,等.開(kāi)敞進(jìn)水池泵站垂直進(jìn)水管道的水力特性研究[J].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,43(6):33-40.

    [14]樊新建,董春海,王之君,等.泵站正向進(jìn)水前池?cái)U(kuò)散角對(duì)池內(nèi)流場(chǎng)結(jié)構(gòu)的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2023,39(1):92-99.

    [15]黃先北,方濤,郭嬙,等.開(kāi)敞式進(jìn)水池喇叭口懸空高度對(duì)吸氣渦的影響[J].工程科學(xué)與技術(shù),2023,55(5):161-168.

    [16]郁片紅,李君菡,付小莉.大型雨水泵站進(jìn)水流道水力特性研究[J].中國(guó)農(nóng)村水利水電,2020(10):126-130.

    [17]施偉,劉軍,陸林廣,等.前置豎井式貫流泵裝置進(jìn)出水流道水力設(shè)計(jì)標(biāo)準(zhǔn)化[J].南水北調(diào)與水利科技(中英文),2020,18(5):115-122.

    [18]車曉紅,張帝,成立.立式混流泵站進(jìn)出水流道水力優(yōu)化[J].南水北調(diào)與水利科技(中英文),2020,18(5):144-150.

    [19]吳鵬飛,郭志偉,錢(qián)忠東,等.進(jìn)水池表面吸氣渦數(shù)值模擬與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2018,49(2):120-125.

    [20]張德勝,李普熙,趙睿杰,等.泵站進(jìn)水池吸入口渦旋結(jié)構(gòu)及湍流特性的大渦模擬[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(10):134-141.

    [21]孔德謙,韓靈杰,李見(jiàn)陽(yáng).水平吸水管懸空高對(duì)泵站進(jìn)水池進(jìn)水特性影響[J].人民長(zhǎng)江,2017,48(5):87-90,100.

    [22]鄭云浩,李彥軍,胡新益,等.基于正交試驗(yàn)的泵站簸箕形進(jìn)水流道型線優(yōu)化[J].人民長(zhǎng)江,2021,52(3):123-126.

    [23]高傳昌,曾新樂(lè),解克宇,等.泵站進(jìn)水池超低水位下組合整流方案與驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2017,33(23):101-108.

    [24]初長(zhǎng)虹,劉超,孫玉民,等.泵站進(jìn)水口淹沒(méi)深度對(duì)漩渦的影響[J].南水北調(diào)與水利科技,2019,17(6):178-186.

    [25]高學(xué)平, 陳思宇,朱洪濤,等.反坡段對(duì)側(cè)式進(jìn)/出水口水力特性影響研究[J].水力發(fā)電學(xué)報(bào),2021,40(5):87-98.

    [26]張晨,李天國(guó), 高學(xué)平.豎井式進(jìn)/出水口工程尺度SPH方法模擬研究[J].水力發(fā)電學(xué)報(bào),2022,41(11):34-45.

    [27]中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.泵站設(shè)計(jì)標(biāo)準(zhǔn):GB 50265—2022[S].北京:中國(guó)計(jì)劃出版社,2022.

    [28]KIM C G,KIM B H,BANG B H,et al.Experimental and CFD analysis for prediction of vortex and swirl angle in the pump sump station model[J].IOP Conference Series:Materials Science and Engineering,2015,72(4):42044.

    [29]Hydraulic Institute.American national standard for rotodynamic pumps for pump instake design:ANSI/HI19.8—2012[S].New Jersey:Hydraulic Institute,2012.

    [30]PARK Y K,DEY M K,CHOI Y H.Numerical visualization of air intake induced by free surface vortex[J].Journal of Thermal Science,2017,26(6):540-544.

    [31]王福軍,唐學(xué)林,陳鑫,等.泵站內(nèi)部流動(dòng)分析方法研究進(jìn)展[J].水利學(xué)報(bào),2018,49(1):47-61,71.

    [32]王福軍.流體機(jī)械旋轉(zhuǎn)湍流計(jì)算模型研究進(jìn)展[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(2):1-14.

    [33]REBOLLO T C,LEWANDOWSKI R.Mathematical and numerical foundations of turbulence models and applications[M].New York:Springer,2014.

    [34]CHIPONGO K,KHIADANI M,SOOKHAK L K.Comparison and verification of turbulence Reynoldsaveraged Navier-Stokes closures to model spatially varied flows[J].Scientific Reports,2020,10(1):19059.

    [35]陸林廣.泵站進(jìn)水流道設(shè)計(jì)理論的新進(jìn)展[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,29(1):40-45.

    Influence of Eccentric Tapered Pipe and Elbow Pipe on the

    Flow Pattern of Inlet Pipe in Pumping Station

    LI Xiaochao1,2, GAO Zhikai3, XIE Minping3, LI Xiaohe4, LUO Yi2,

    TANG Hualin2, LIN Xiaoheng5, LI Jun3

    (1.China Institute of Water Resources and Hydropower Research, Beijing 100038, China;

    2.China Water Hydropower Development Co., Ltd., Chengdu 610000, China;

    3.North China University of Water Resources and Electric Power, Zhengzhou 450045, China;

    4.South China Normal University, Guangzhou 510631, China;

    5.The Fourth Construction Co., Ltd. of China Construction Eighth Engineering Bureau, Qingdao 266100, China)

    Abstract:

    To investigate the effects of different eccentric tapered pipe lengths and elbow angles on the hydraulic characteristics of the inlet pipe in open pumping stations, this study explores the variation of hydraulic performance parameters with changes in elbow angles and eccentric tapered pipe lengths. Using the Yellow River pumping station in the lower reaches of the Yellow River as a case study, numerical simulations of 14 schemes were conducted based on Fluent software. The flow field of the inlet pipe at typical sections was analyzed using evaluation functions. Comparative analysis of schemes 1-5 reveals that changes in the length of the eccentric tapered pipe significantly affect the hydraulic loss of the inlet pipe, the uniformity of flow velocity distribution, and the velocity-weighted average angle at the outlet section. Optimal hydraulic performance is observed when the eccentric tapered pipe lengths are 600 mm, 700 mm, and 800 mm. The calculation and analysis of scheme 6-14 show that the change of elbow angle mainly affects the flow pattern of elbow section, and has little disturbance to the flow of other parts. Further evaluation shows that larger elbow angles result in smoother water flow transitions and lower hydraulic losses, whereas smaller elbow angles increase hydraulic losses, whereas smaller elbow angles increase hydraulic losses. Additionally, longer eccentric tapered pipe lengths improve the uniformity of flow velocity distribution and the velocity-weighted average angle at the outlet section, enhancing rectification effects, but at the cost of increased hydraulic losses.

    Keywords:

    pumping station; elbow pipe; eccentric tapered pipe; inflow flow pattern; numerical simulation

    (編輯:?jiǎn)檀淦剑?/p>

    收稿日期:2023-07-04

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51909094)。

    第一作者:

    李曉超(1997—),男,博士研究生,從事水利水電工程方面的研究。E-mail:hnscyxlxc@163.com。

    通信作者:李君(1979—),男,副教授,碩導(dǎo),博士,從事水利工程方面的研究。E-mail:uniquelijun@ncwu.edu.cn。

    猜你喜歡
    數(shù)值模擬泵站
    張家邊涌泵站建設(shè)難點(diǎn)及技術(shù)創(chuàng)新實(shí)踐
    2016年河南省己建成泵站數(shù)量
    張家灣煤礦巷道無(wú)支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開(kāi)切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設(shè)計(jì)與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    全省已建成泵站數(shù)量
    河南省2014年已建成泵站數(shù)量
    日本在线视频免费播放| 男人舔奶头视频| 午夜a级毛片| 色综合站精品国产| 国产中年淑女户外野战色| 久久人妻av系列| 黄色一级大片看看| 一a级毛片在线观看| 嫩草影院精品99| 精品久久久久久久末码| 国产大屁股一区二区在线视频| 伦理电影大哥的女人| 又黄又爽又免费观看的视频| 黄片wwwwww| 狂野欧美激情性xxxx在线观看| 一本久久中文字幕| 久久久久久久亚洲中文字幕| 亚洲国产高清在线一区二区三| 久久精品国产亚洲av香蕉五月| 日日干狠狠操夜夜爽| 欧美丝袜亚洲另类| 国产免费男女视频| 男人舔女人下体高潮全视频| 少妇猛男粗大的猛烈进出视频 | 人人妻人人看人人澡| 成人av一区二区三区在线看| 变态另类丝袜制服| 99久久精品国产国产毛片| 欧美高清性xxxxhd video| 桃色一区二区三区在线观看| 男女下面进入的视频免费午夜| 亚洲一区二区三区色噜噜| 一级黄色大片毛片| 变态另类成人亚洲欧美熟女| 欧美区成人在线视频| 亚洲经典国产精华液单| 久久精品91蜜桃| 成人毛片a级毛片在线播放| 亚洲激情五月婷婷啪啪| 亚洲国产精品sss在线观看| 搡女人真爽免费视频火全软件 | 国产av不卡久久| 91久久精品国产一区二区三区| 亚洲av第一区精品v没综合| 高清毛片免费观看视频网站| 亚洲国产色片| 亚洲成a人片在线一区二区| av专区在线播放| 99久久久亚洲精品蜜臀av| 日韩三级伦理在线观看| 国产美女午夜福利| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩无卡精品| 在线观看美女被高潮喷水网站| 夜夜夜夜夜久久久久| 91在线观看av| 美女免费视频网站| videossex国产| 午夜福利成人在线免费观看| 精品福利观看| 亚州av有码| 免费观看的影片在线观看| 国产高清三级在线| 中文字幕av成人在线电影| 国产一区二区在线av高清观看| 国产人妻一区二区三区在| 成人一区二区视频在线观看| 波多野结衣高清无吗| 99国产极品粉嫩在线观看| a级毛片免费高清观看在线播放| 欧美xxxx性猛交bbbb| 成人精品一区二区免费| 精品人妻视频免费看| 在线观看美女被高潮喷水网站| 国产精品伦人一区二区| 国产老妇女一区| 亚洲av.av天堂| 男插女下体视频免费在线播放| 丝袜喷水一区| 日韩人妻高清精品专区| 日韩欧美一区二区三区在线观看| 十八禁网站免费在线| 深夜精品福利| 国产国拍精品亚洲av在线观看| 综合色av麻豆| 亚洲av熟女| 日本欧美国产在线视频| 在线免费观看不下载黄p国产| 国产亚洲欧美98| 观看免费一级毛片| 我要看日韩黄色一级片| 久久午夜福利片| 精品免费久久久久久久清纯| 日日干狠狠操夜夜爽| 狂野欧美白嫩少妇大欣赏| 久久久欧美国产精品| 97碰自拍视频| 欧美zozozo另类| 国产91av在线免费观看| 麻豆久久精品国产亚洲av| 精品久久国产蜜桃| 女的被弄到高潮叫床怎么办| 日本黄色视频三级网站网址| 久久久精品欧美日韩精品| 成人二区视频| 网址你懂的国产日韩在线| 一进一出好大好爽视频| 国产麻豆成人av免费视频| 欧美日本视频| 国产高清激情床上av| 99riav亚洲国产免费| 国内精品久久久久精免费| 国产精品综合久久久久久久免费| 亚洲av熟女| 中国美女看黄片| 51国产日韩欧美| 国产在视频线在精品| 精品一区二区三区视频在线观看免费| 国产成人福利小说| 夜夜夜夜夜久久久久| 久久久久国产精品人妻aⅴ院| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡免费网站照片| 国产综合懂色| av在线观看视频网站免费| а√天堂www在线а√下载| 亚洲图色成人| 人妻夜夜爽99麻豆av| 中文字幕熟女人妻在线| 亚洲第一电影网av| 国产免费一级a男人的天堂| 欧美性猛交黑人性爽| 亚洲一区高清亚洲精品| 成人特级黄色片久久久久久久| 国产精品一区二区三区四区久久| 国产69精品久久久久777片| 看免费成人av毛片| 天堂√8在线中文| 久久人人爽人人片av| 久久久久国产网址| 亚洲熟妇熟女久久| 精品日产1卡2卡| 日本成人三级电影网站| 有码 亚洲区| 免费看av在线观看网站| 变态另类成人亚洲欧美熟女| 国产男人的电影天堂91| 黄色日韩在线| or卡值多少钱| 丝袜喷水一区| 亚洲av中文av极速乱| 免费观看人在逋| 99热全是精品| 亚洲欧美成人精品一区二区| 又黄又爽又免费观看的视频| 不卡一级毛片| 欧美日韩在线观看h| 亚洲无线在线观看| 国产亚洲精品久久久com| 亚洲久久久久久中文字幕| 欧美日本视频| 久久精品国产亚洲av涩爱 | 男人舔女人下体高潮全视频| 偷拍熟女少妇极品色| 美女高潮的动态| 99热精品在线国产| 色哟哟哟哟哟哟| 日韩欧美三级三区| 亚洲婷婷狠狠爱综合网| 老司机福利观看| 国产精品久久久久久久电影| 丰满乱子伦码专区| 久久精品国产鲁丝片午夜精品| 久久久久久九九精品二区国产| 日本a在线网址| 观看免费一级毛片| 国产一区二区在线观看日韩| 午夜精品在线福利| 又黄又爽又免费观看的视频| a级毛片a级免费在线| 欧美极品一区二区三区四区| 嫩草影视91久久| 亚洲国产精品成人综合色| 成人美女网站在线观看视频| 国产真实伦视频高清在线观看| 亚洲婷婷狠狠爱综合网| 亚洲天堂国产精品一区在线| 免费观看精品视频网站| 久久韩国三级中文字幕| 一级黄片播放器| 亚洲精品456在线播放app| 看黄色毛片网站| 亚洲电影在线观看av| 亚洲欧美成人综合另类久久久 | av在线亚洲专区| 国内少妇人妻偷人精品xxx网站| 99久久中文字幕三级久久日本| 日韩av在线大香蕉| 免费看av在线观看网站| 亚洲成av人片在线播放无| 99久久久亚洲精品蜜臀av| 国产精品国产三级国产av玫瑰| 六月丁香七月| 一卡2卡三卡四卡精品乱码亚洲| 欧美+日韩+精品| 国产三级在线视频| 日日干狠狠操夜夜爽| 亚洲精品乱码久久久v下载方式| 在线观看66精品国产| 国产精品三级大全| 日本在线视频免费播放| 深夜a级毛片| 美女cb高潮喷水在线观看| 2021天堂中文幕一二区在线观| av在线播放精品| 久久精品国产鲁丝片午夜精品| 欧美区成人在线视频| 亚洲自拍偷在线| 国国产精品蜜臀av免费| 亚洲va在线va天堂va国产| 女的被弄到高潮叫床怎么办| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 深夜精品福利| 美女cb高潮喷水在线观看| 观看免费一级毛片| 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 丝袜美腿在线中文| 老司机福利观看| 日韩精品中文字幕看吧| av国产免费在线观看| 天美传媒精品一区二区| 国产单亲对白刺激| 亚洲精品亚洲一区二区| 99在线视频只有这里精品首页| 少妇高潮的动态图| 免费观看人在逋| 青春草视频在线免费观看| 看非洲黑人一级黄片| 成人性生交大片免费视频hd| 欧美日韩一区二区视频在线观看视频在线 | 美女大奶头视频| 麻豆国产av国片精品| 可以在线观看毛片的网站| 欧美一区二区精品小视频在线| 国产精品av视频在线免费观看| 一个人免费在线观看电影| 国产麻豆成人av免费视频| 日韩欧美在线乱码| 国内揄拍国产精品人妻在线| 啦啦啦啦在线视频资源| 日韩精品青青久久久久久| 久久久久国内视频| 亚洲色图av天堂| 久久久精品欧美日韩精品| 亚洲在线自拍视频| 久久久久久伊人网av| 国产单亲对白刺激| 2021天堂中文幕一二区在线观| 久久午夜亚洲精品久久| 国产精品久久久久久久电影| 美女内射精品一级片tv| 日韩精品中文字幕看吧| 日韩欧美免费精品| 在线播放无遮挡| 波多野结衣高清作品| 高清毛片免费看| 午夜福利18| 国产真实乱freesex| 成人亚洲欧美一区二区av| 欧美日韩乱码在线| 国产午夜精品论理片| 99在线人妻在线中文字幕| 亚洲美女搞黄在线观看 | 九九爱精品视频在线观看| 国产三级在线视频| 国产日本99.免费观看| 免费观看在线日韩| 成人无遮挡网站| 日韩,欧美,国产一区二区三区 | 亚洲av美国av| 一本精品99久久精品77| 国产精品av视频在线免费观看| 欧美丝袜亚洲另类| 日韩精品青青久久久久久| 国内精品宾馆在线| 国产熟女欧美一区二区| 六月丁香七月| 国语自产精品视频在线第100页| 18禁在线播放成人免费| 蜜桃久久精品国产亚洲av| 精品熟女少妇av免费看| 男女下面进入的视频免费午夜| 日韩av在线大香蕉| 亚洲人成网站高清观看| 午夜福利在线观看免费完整高清在 | 日日撸夜夜添| 直男gayav资源| 日日摸夜夜添夜夜爱| 亚洲人成网站高清观看| 国产蜜桃级精品一区二区三区| 久久草成人影院| 亚洲精品一区av在线观看| 久99久视频精品免费| 国产精品久久久久久久电影| 一区二区三区四区激情视频 | 村上凉子中文字幕在线| 日本撒尿小便嘘嘘汇集6| 12—13女人毛片做爰片一| 午夜免费激情av| 国产白丝娇喘喷水9色精品| 免费观看人在逋| 久久久欧美国产精品| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 久久久成人免费电影| 99久国产av精品| 日产精品乱码卡一卡2卡三| 简卡轻食公司| ponron亚洲| 国产精品嫩草影院av在线观看| 亚洲专区国产一区二区| 久久午夜福利片| 99久久九九国产精品国产免费| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 精品国产三级普通话版| 内地一区二区视频在线| 亚洲美女视频黄频| 国产久久久一区二区三区| 搡老岳熟女国产| 成人鲁丝片一二三区免费| 神马国产精品三级电影在线观看| 久久人人爽人人片av| 国产乱人视频| 在线免费观看不下载黄p国产| 精华霜和精华液先用哪个| 亚洲成a人片在线一区二区| 亚洲精品亚洲一区二区| aaaaa片日本免费| 免费在线观看影片大全网站| 亚洲国产精品成人久久小说 | 亚洲av成人av| 淫妇啪啪啪对白视频| 在线免费观看的www视频| 国产乱人偷精品视频| 婷婷六月久久综合丁香| 中文字幕人妻熟人妻熟丝袜美| 搡女人真爽免费视频火全软件 | 国产精品福利在线免费观看| 免费看av在线观看网站| 亚洲精品成人久久久久久| 最新在线观看一区二区三区| 美女大奶头视频| 又爽又黄无遮挡网站| 村上凉子中文字幕在线| 99久国产av精品| 一级a爱片免费观看的视频| 又爽又黄无遮挡网站| 午夜影院日韩av| 成人毛片a级毛片在线播放| 一级毛片我不卡| 国产高清不卡午夜福利| 午夜免费激情av| 国产精品久久久久久久电影| 老司机影院成人| 国产一级毛片七仙女欲春2| 免费观看的影片在线观看| 久久精品91蜜桃| 99久久精品国产国产毛片| 亚洲电影在线观看av| 春色校园在线视频观看| 三级国产精品欧美在线观看| 欧美在线一区亚洲| 麻豆一二三区av精品| 国产精品一区二区三区四区久久| 联通29元200g的流量卡| 九九久久精品国产亚洲av麻豆| 狂野欧美白嫩少妇大欣赏| 嫩草影视91久久| 99热全是精品| 日本色播在线视频| 亚洲美女搞黄在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色| 丰满乱子伦码专区| 成人永久免费在线观看视频| 色哟哟·www| 中文字幕久久专区| 久久精品综合一区二区三区| 久久亚洲国产成人精品v| 亚洲性夜色夜夜综合| 色5月婷婷丁香| 男女之事视频高清在线观看| 日韩精品中文字幕看吧| 一级毛片我不卡| 欧美三级亚洲精品| 国产成人aa在线观看| 午夜福利高清视频| aaaaa片日本免费| 日韩欧美精品v在线| 欧美中文日本在线观看视频| 日韩精品有码人妻一区| 国产综合懂色| 麻豆国产97在线/欧美| 麻豆一二三区av精品| 最近视频中文字幕2019在线8| 日韩成人伦理影院| 国产高清有码在线观看视频| 在线观看66精品国产| 美女cb高潮喷水在线观看| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 国产三级在线视频| 免费黄网站久久成人精品| 亚洲综合色惰| 国产探花在线观看一区二区| 中文字幕熟女人妻在线| 久久鲁丝午夜福利片| 色播亚洲综合网| av天堂中文字幕网| 国产片特级美女逼逼视频| 久久精品人妻少妇| 老司机午夜福利在线观看视频| 男人舔女人下体高潮全视频| 最近最新中文字幕大全电影3| 五月玫瑰六月丁香| 在线国产一区二区在线| 91在线精品国自产拍蜜月| 亚洲av二区三区四区| 亚洲欧美清纯卡通| 大型黄色视频在线免费观看| 女人被狂操c到高潮| 午夜a级毛片| 丝袜喷水一区| a级毛片免费高清观看在线播放| 国产成人一区二区在线| 美女内射精品一级片tv| 1000部很黄的大片| 狂野欧美白嫩少妇大欣赏| 日韩欧美精品v在线| 高清毛片免费看| 黑人高潮一二区| a级一级毛片免费在线观看| 18禁在线播放成人免费| 一个人看的www免费观看视频| 免费观看的影片在线观看| 天天躁夜夜躁狠狠久久av| 成年免费大片在线观看| 草草在线视频免费看| 99热这里只有是精品50| 我要看日韩黄色一级片| 国产黄片美女视频| 亚洲天堂国产精品一区在线| 色综合亚洲欧美另类图片| 人妻丰满熟妇av一区二区三区| 亚洲无线在线观看| 伦理电影大哥的女人| 禁无遮挡网站| 国产av在哪里看| 午夜精品国产一区二区电影 | 我要看日韩黄色一级片| 亚洲色图av天堂| 亚洲四区av| 亚洲激情五月婷婷啪啪| 午夜福利成人在线免费观看| 国产一区亚洲一区在线观看| 91在线观看av| 久久天躁狠狠躁夜夜2o2o| 亚洲一区二区三区色噜噜| 韩国av在线不卡| 最新中文字幕久久久久| 中文字幕久久专区| 午夜亚洲福利在线播放| 九九在线视频观看精品| 99热6这里只有精品| 一a级毛片在线观看| 精品一区二区三区av网在线观看| av在线天堂中文字幕| 国语自产精品视频在线第100页| 久久久色成人| 国产一区二区三区av在线 | 欧美日韩一区二区视频在线观看视频在线 | 97热精品久久久久久| 精品一区二区三区人妻视频| 成人鲁丝片一二三区免费| 亚洲综合色惰| 在线国产一区二区在线| 亚洲五月天丁香| 99视频精品全部免费 在线| 插阴视频在线观看视频| 日韩精品中文字幕看吧| 亚洲最大成人中文| 偷拍熟女少妇极品色| 日本一本二区三区精品| 少妇人妻一区二区三区视频| 欧美中文日本在线观看视频| 亚洲精品在线观看二区| 97超级碰碰碰精品色视频在线观看| 22中文网久久字幕| 欧美在线一区亚洲| 欧美最新免费一区二区三区| 男女视频在线观看网站免费| 亚洲专区国产一区二区| 亚洲欧美清纯卡通| 成年女人看的毛片在线观看| 高清毛片免费观看视频网站| 高清午夜精品一区二区三区 | 97热精品久久久久久| 国产精品亚洲一级av第二区| 久久精品国产自在天天线| 精品日产1卡2卡| 久久精品国产自在天天线| 久久久国产成人精品二区| 亚洲成人久久性| 亚洲国产精品成人久久小说 | 成人国产麻豆网| 欧洲精品卡2卡3卡4卡5卡区| 搡女人真爽免费视频火全软件 | 国产免费男女视频| 亚洲国产精品成人综合色| 国产精品亚洲一级av第二区| 最近在线观看免费完整版| 久久精品91蜜桃| 国产欧美日韩一区二区精品| 国产三级中文精品| 久久久国产成人精品二区| 久久精品国产自在天天线| 国产久久久一区二区三区| 国产不卡一卡二| 久久精品国产亚洲av天美| 美女黄网站色视频| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件 | 白带黄色成豆腐渣| 12—13女人毛片做爰片一| 插阴视频在线观看视频| 亚洲国产高清在线一区二区三| 一级黄色大片毛片| 欧美性感艳星| 欧美高清成人免费视频www| 18禁黄网站禁片免费观看直播| 免费在线观看成人毛片| 亚洲激情五月婷婷啪啪| 99久久无色码亚洲精品果冻| 亚洲,欧美,日韩| 男女之事视频高清在线观看| 女同久久另类99精品国产91| 99久国产av精品国产电影| 内射极品少妇av片p| 激情 狠狠 欧美| 久久鲁丝午夜福利片| 一级毛片电影观看 | 国产一区二区在线观看日韩| 欧美国产日韩亚洲一区| 99久久精品热视频| 真实男女啪啪啪动态图| 国产一区二区激情短视频| 国产伦一二天堂av在线观看| 一个人观看的视频www高清免费观看| 精品熟女少妇av免费看| aaaaa片日本免费| 热99在线观看视频| 亚洲天堂国产精品一区在线| 亚洲内射少妇av| 国产亚洲精品av在线| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 亚洲人与动物交配视频| 国产男靠女视频免费网站| 中国国产av一级| 亚洲成av人片在线播放无| 亚洲七黄色美女视频| 九九在线视频观看精品| 国产成人aa在线观看| 无遮挡黄片免费观看| 亚洲国产精品成人久久小说 | 欧美日本亚洲视频在线播放| 国产v大片淫在线免费观看| 免费av毛片视频| 天堂av国产一区二区熟女人妻| 日本免费a在线| 少妇高潮的动态图| 99久久精品国产国产毛片| 国产爱豆传媒在线观看| 少妇熟女欧美另类| av国产免费在线观看| 日本三级黄在线观看| 久久精品国产亚洲av香蕉五月| 日本一本二区三区精品| 麻豆国产97在线/欧美| 欧美人与善性xxx| 香蕉av资源在线| 欧美xxxx性猛交bbbb| 国产不卡一卡二| 国产色婷婷99| 亚洲18禁久久av| 亚洲国产日韩欧美精品在线观看| 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站 | 欧美xxxx黑人xx丫x性爽| 久久精品人妻少妇| 成熟少妇高潮喷水视频| 国产精品国产三级国产av玫瑰| 波多野结衣巨乳人妻| 69人妻影院| 亚洲中文字幕日韩| 最近2019中文字幕mv第一页| 国产av在哪里看| 欧美国产日韩亚洲一区| 天堂动漫精品| 日韩欧美在线乱码| 99国产精品一区二区蜜桃av|