• <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级片在线免费高清观看视频| 色在线成人网| 久久久水蜜桃国产精品网| 母亲3免费完整高清在线观看| 无人区码免费观看不卡| 国产成人av激情在线播放| 在线十欧美十亚洲十日本专区| 999精品在线视频| 久久国产乱子伦精品免费另类| a级片在线免费高清观看视频| 亚洲人成77777在线视频| 国产精品成人在线| 久久久久久久久免费视频了| 亚洲欧美精品综合一区二区三区| 波多野结衣av一区二区av| 99在线人妻在线中文字幕 | 亚洲欧美一区二区三区久久| 99久久人妻综合| 亚洲精品久久成人aⅴ小说| 免费黄频网站在线观看国产| 国产欧美日韩一区二区精品| 天堂动漫精品| 欧美+亚洲+日韩+国产| 一区二区三区激情视频| 黄网站色视频无遮挡免费观看| 国产视频一区二区在线看| 亚洲欧美日韩另类电影网站| 精品福利观看| 无限看片的www在线观看| av电影中文网址| 搡老熟女国产l中国老女人| 久久青草综合色| 免费日韩欧美在线观看| 久热这里只有精品99| 国产成人精品无人区| 女人精品久久久久毛片| 热re99久久国产66热| 午夜精品国产一区二区电影| 国产av精品麻豆| 国产精品秋霞免费鲁丝片| 丁香六月欧美| 女人被躁到高潮嗷嗷叫费观| 亚洲人成77777在线视频| 午夜老司机福利片| 午夜福利影视在线免费观看| 日韩精品免费视频一区二区三区| 日韩大码丰满熟妇| 一级毛片女人18水好多| 妹子高潮喷水视频| 亚洲精品美女久久久久99蜜臀| 久久九九热精品免费| 最新美女视频免费是黄的| 亚洲精品一卡2卡三卡4卡5卡| 国产淫语在线视频| 高清毛片免费观看视频网站 | 欧美一级毛片孕妇| 九色亚洲精品在线播放| 国产精品久久久人人做人人爽| 久久亚洲真实| 多毛熟女@视频| 水蜜桃什么品种好| 国产精品久久久人人做人人爽| 久久精品91无色码中文字幕| 我的亚洲天堂| 99国产极品粉嫩在线观看| 国产精品一区二区在线观看99| 欧美日韩一级在线毛片| 日日摸夜夜添夜夜添小说| 欧美 日韩 精品 国产| 精品久久久久久久久久免费视频 | 色播在线永久视频| 久久国产精品男人的天堂亚洲| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 午夜91福利影院| 男女午夜视频在线观看| 在线看a的网站| www.精华液| 美女国产高潮福利片在线看| 国产av一区二区精品久久| 日韩人妻精品一区2区三区| √禁漫天堂资源中文www| 午夜福利影视在线免费观看| 亚洲一区二区三区不卡视频| av网站免费在线观看视频| 亚洲av日韩在线播放| 久久久久久久久久久久大奶| 亚洲av日韩在线播放| 美女高潮到喷水免费观看| 日韩欧美一区视频在线观看| 91成人精品电影| 搡老岳熟女国产| 欧美日韩亚洲高清精品| 身体一侧抽搐| 久久久久久久国产电影| 极品少妇高潮喷水抽搐| 日本精品一区二区三区蜜桃| av超薄肉色丝袜交足视频| 国产熟女午夜一区二区三区| 亚洲色图 男人天堂 中文字幕| 国产日韩欧美亚洲二区| 国产在线一区二区三区精| 国产又色又爽无遮挡免费看| 国产真人三级小视频在线观看| 久久久久精品人妻al黑| 亚洲精品中文字幕一二三四区| 国产极品粉嫩免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人一区二区三| 中文亚洲av片在线观看爽 | 免费人成视频x8x8入口观看| 男人操女人黄网站| 国产精品免费大片| 91九色精品人成在线观看| 无遮挡黄片免费观看| 亚洲国产精品一区二区三区在线| 高潮久久久久久久久久久不卡| 女人被躁到高潮嗷嗷叫费观| 国产有黄有色有爽视频| 黑人猛操日本美女一级片| 亚洲成人免费电影在线观看| 在线天堂中文资源库| 色婷婷av一区二区三区视频| 久热爱精品视频在线9| 亚洲av美国av| 夜夜夜夜夜久久久久| 色老头精品视频在线观看| 久久久精品免费免费高清| 亚洲人成电影观看| 一本综合久久免费| 高清视频免费观看一区二区| 18禁黄网站禁片午夜丰满| 在线观看舔阴道视频| 极品少妇高潮喷水抽搐| 91精品三级在线观看| 免费少妇av软件| 亚洲一区中文字幕在线| ponron亚洲| 亚洲成人免费av在线播放| 黄色毛片三级朝国网站| 夫妻午夜视频| 日韩精品免费视频一区二区三区| 一级a爱片免费观看的视频| 丝袜美腿诱惑在线| 久久久久久亚洲精品国产蜜桃av| 精品卡一卡二卡四卡免费| 12—13女人毛片做爰片一| 男女床上黄色一级片免费看| 69av精品久久久久久| 亚洲国产欧美一区二区综合| 日韩大码丰满熟妇| 欧美日韩视频精品一区| 免费在线观看黄色视频的| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧洲精品一区二区精品久久久| 国产免费男女视频| 色在线成人网| 午夜福利在线免费观看网站| 免费人成视频x8x8入口观看| 精品熟女少妇八av免费久了| 日韩欧美在线二视频 | 亚洲一区中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 国产不卡av网站在线观看| 交换朋友夫妻互换小说| 亚洲第一欧美日韩一区二区三区| 亚洲中文日韩欧美视频| 亚洲精品国产区一区二| 亚洲视频免费观看视频| 中亚洲国语对白在线视频| 午夜精品在线福利| 亚洲avbb在线观看| 99久久99久久久精品蜜桃| 9191精品国产免费久久| 搡老熟女国产l中国老女人| 亚洲一区中文字幕在线| 国产亚洲精品第一综合不卡| 99精品欧美一区二区三区四区| 精品福利永久在线观看| aaaaa片日本免费| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 国产成人精品久久二区二区91| 免费在线观看影片大全网站| 99热只有精品国产| 国产不卡av网站在线观看| 亚洲av熟女| 亚洲人成77777在线视频| 久久人妻熟女aⅴ| 一区在线观看完整版| 国产成人啪精品午夜网站| 大片电影免费在线观看免费| 啪啪无遮挡十八禁网站| 中亚洲国语对白在线视频| 两性夫妻黄色片| 亚洲免费av在线视频| 亚洲精品国产色婷婷电影| 久久ye,这里只有精品| 丝袜在线中文字幕| 精品福利永久在线观看| 成人国语在线视频| 精品第一国产精品| 国产片内射在线| 久久精品国产99精品国产亚洲性色 | 久久精品国产99精品国产亚洲性色 | 日韩人妻精品一区2区三区| 色94色欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 一级片免费观看大全| 亚洲国产欧美一区二区综合| 亚洲成a人片在线一区二区| 亚洲国产欧美网| 久久久久久人人人人人| 一级a爱片免费观看的视频| 久久香蕉精品热| 久久精品国产a三级三级三级| 久热这里只有精品99| videosex国产| 亚洲七黄色美女视频| 亚洲色图 男人天堂 中文字幕| 精品一区二区三区av网在线观看| 午夜福利欧美成人| 欧美精品一区二区免费开放| 人人澡人人妻人| 麻豆成人av在线观看| 国产精品久久久久久精品古装| 欧美成人午夜精品| 999久久久精品免费观看国产| 成人国产一区最新在线观看| 欧美av亚洲av综合av国产av| 啦啦啦 在线观看视频| 在线观看免费午夜福利视频| 黄色怎么调成土黄色| 亚洲精品一卡2卡三卡4卡5卡| 黑人巨大精品欧美一区二区mp4| 亚洲国产毛片av蜜桃av| 国产免费现黄频在线看| 国产精品美女特级片免费视频播放器 | 免费看十八禁软件| 大香蕉久久成人网| 久久精品成人免费网站| 久99久视频精品免费| 亚洲国产中文字幕在线视频| 国产99久久九九免费精品| 亚洲成人国产一区在线观看| 人人妻人人澡人人爽人人夜夜| 交换朋友夫妻互换小说| tocl精华| 免费久久久久久久精品成人欧美视频| 在线国产一区二区在线| 超碰97精品在线观看| 精品国产美女av久久久久小说| 少妇裸体淫交视频免费看高清 | 女性被躁到高潮视频| a级毛片黄视频| 午夜精品国产一区二区电影| 国产又爽黄色视频| 亚洲欧美一区二区三区黑人| 精品午夜福利视频在线观看一区| 精品一区二区三区av网在线观看| 国产成人欧美| 成在线人永久免费视频| 国产精品电影一区二区三区 | 国产精品成人在线| 欧美精品av麻豆av| 亚洲色图av天堂| 国产一区有黄有色的免费视频| 他把我摸到了高潮在线观看| 一本一本久久a久久精品综合妖精| 国产成人精品无人区| 久久久精品国产亚洲av高清涩受| 久久久久久久久免费视频了| 国产精品乱码一区二三区的特点 | 好看av亚洲va欧美ⅴa在| 婷婷精品国产亚洲av在线 | 亚洲一区高清亚洲精品| 色尼玛亚洲综合影院| 国产在视频线精品| 久久精品国产a三级三级三级| 国产一区二区三区在线臀色熟女 | 交换朋友夫妻互换小说| 在线观看免费高清a一片| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 我的亚洲天堂| 午夜福利影视在线免费观看| 国产精品免费视频内射| 999久久久精品免费观看国产| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 91麻豆av在线| 美国免费a级毛片| 日韩欧美国产一区二区入口| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区mp4| 美女高潮到喷水免费观看| 精品国产一区二区久久| 黄色毛片三级朝国网站| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看 | 免费女性裸体啪啪无遮挡网站| 国产欧美日韩一区二区精品| 国产精品一区二区免费欧美| 美女国产高潮福利片在线看| 欧美黑人精品巨大| 18禁黄网站禁片午夜丰满| 亚洲精品成人av观看孕妇| 亚洲精品国产精品久久久不卡| x7x7x7水蜜桃| 亚洲av片天天在线观看| 天堂俺去俺来也www色官网| 婷婷精品国产亚洲av在线 | 久久久国产成人精品二区 | 国产片内射在线| 又黄又粗又硬又大视频| 欧美老熟妇乱子伦牲交| www.精华液| 成人国产一区最新在线观看| 国产亚洲精品久久久久5区| 国内久久婷婷六月综合欲色啪| 欧美日韩av久久| 国产精品香港三级国产av潘金莲| 国产有黄有色有爽视频| 久久这里只有精品19| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 亚洲九九香蕉| 操美女的视频在线观看| 欧美不卡视频在线免费观看 | 久久香蕉激情| 精品第一国产精品| 老司机靠b影院| 99精品欧美一区二区三区四区| 18禁黄网站禁片午夜丰满| 黄色毛片三级朝国网站| 日韩制服丝袜自拍偷拍| 国产成人精品无人区| av国产精品久久久久影院| 村上凉子中文字幕在线| 国产精品.久久久| 久久香蕉激情| 夜夜躁狠狠躁天天躁| 精品国产超薄肉色丝袜足j| 天堂√8在线中文| 久久天躁狠狠躁夜夜2o2o| 狂野欧美激情性xxxx| 欧美日韩瑟瑟在线播放| 黄片小视频在线播放| 欧美日韩一级在线毛片| 在线观看免费高清a一片| 亚洲国产精品合色在线| 极品少妇高潮喷水抽搐| 午夜福利在线免费观看网站| 久久久国产成人精品二区 | 欧美最黄视频在线播放免费 | 久热爱精品视频在线9| 国产精品香港三级国产av潘金莲| 免费高清在线观看日韩| 日韩一卡2卡3卡4卡2021年| av国产精品久久久久影院| 精品人妻1区二区| 最近最新免费中文字幕在线| 亚洲精品国产一区二区精华液| 曰老女人黄片| 国产一区二区激情短视频| 亚洲色图综合在线观看| 国产三级黄色录像| 搡老岳熟女国产| 久久久精品免费免费高清| 久久久久久久午夜电影 | 最近最新中文字幕大全免费视频| 99热国产这里只有精品6| 亚洲欧美一区二区三区久久| 亚洲片人在线观看| 免费观看a级毛片全部| 精品久久久精品久久久| 久久久精品免费免费高清| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 亚洲成人免费电影在线观看| 啦啦啦 在线观看视频| 女人被狂操c到高潮| 丝袜美足系列| 欧美日韩黄片免| 深夜精品福利| 久久午夜亚洲精品久久| 人妻丰满熟妇av一区二区三区 | 免费看a级黄色片| 男女下面插进去视频免费观看| 在线观看免费高清a一片| 亚洲国产看品久久| 午夜精品在线福利| 亚洲中文日韩欧美视频| 精品人妻在线不人妻| 麻豆国产av国片精品| 一本大道久久a久久精品| 精品电影一区二区在线| 午夜福利视频在线观看免费| 精品福利永久在线观看| 亚洲精品国产色婷婷电影| 每晚都被弄得嗷嗷叫到高潮| 黑丝袜美女国产一区| 很黄的视频免费| 亚洲国产精品合色在线| 啦啦啦 在线观看视频| 高清毛片免费观看视频网站 | 久久精品aⅴ一区二区三区四区| 夜夜爽天天搞| 免费在线观看日本一区| 在线视频色国产色| 男人操女人黄网站| 99久久99久久久精品蜜桃| 午夜福利一区二区在线看| 久久久久久亚洲精品国产蜜桃av| 国产精品免费视频内射| 女人久久www免费人成看片| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 亚洲精品乱久久久久久| netflix在线观看网站| 一级a爱片免费观看的视频| ponron亚洲| 日本黄色视频三级网站网址 | 十分钟在线观看高清视频www| 777久久人妻少妇嫩草av网站| 在线天堂中文资源库| 99久久综合精品五月天人人| 国产成人欧美| 99国产精品一区二区蜜桃av | 精品人妻熟女毛片av久久网站| 人妻久久中文字幕网| av天堂在线播放| 欧美色视频一区免费| 一本综合久久免费| 日韩欧美一区二区三区在线观看 | 国产精品久久电影中文字幕 | 一进一出好大好爽视频| 国产成人欧美在线观看 | 亚洲五月天丁香| 精品国产国语对白av| 男女床上黄色一级片免费看| 啦啦啦视频在线资源免费观看| 91成年电影在线观看| 精品国内亚洲2022精品成人 | 天天添夜夜摸| 欧美黄色淫秽网站| 女性被躁到高潮视频| 日本精品一区二区三区蜜桃| 久久久国产成人精品二区 | 变态另类成人亚洲欧美熟女 | 精品国产一区二区三区久久久樱花| 在线观看免费午夜福利视频| 午夜福利影视在线免费观看| 两个人免费观看高清视频| 两性夫妻黄色片| 激情在线观看视频在线高清 | 成年女人毛片免费观看观看9 | 老司机影院毛片| 人妻丰满熟妇av一区二区三区 | 99国产精品一区二区蜜桃av | 国产精品免费一区二区三区在线 | 18禁裸乳无遮挡动漫免费视频| 一二三四在线观看免费中文在| 亚洲九九香蕉| 精品高清国产在线一区| 久久精品国产清高在天天线| av网站在线播放免费| 国产在线观看jvid| 在线观看日韩欧美| 一级黄色大片毛片| 亚洲黑人精品在线| 一区二区三区国产精品乱码| 久久天堂一区二区三区四区| 99久久国产精品久久久| 欧美 亚洲 国产 日韩一| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区免费开放| 9热在线视频观看99| 叶爱在线成人免费视频播放| 免费在线观看日本一区| 80岁老熟妇乱子伦牲交| 精品久久蜜臀av无| 一进一出抽搐gif免费好疼 | 免费在线观看黄色视频的| 麻豆乱淫一区二区| 免费人成视频x8x8入口观看| 欧美国产精品一级二级三级| 成人av一区二区三区在线看| 交换朋友夫妻互换小说| 欧美一级毛片孕妇| 国产精品免费视频内射| 久久人妻福利社区极品人妻图片| 久久国产亚洲av麻豆专区| 国产蜜桃级精品一区二区三区 | 欧美精品亚洲一区二区| 久久香蕉国产精品| 久久热在线av| 国产欧美日韩一区二区三| 欧美成人免费av一区二区三区 | 交换朋友夫妻互换小说| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品二区激情视频| 国产区一区二久久| 手机成人av网站| 香蕉丝袜av| 91在线观看av| 亚洲久久久国产精品| 91九色精品人成在线观看| 露出奶头的视频| 精品国内亚洲2022精品成人 | 露出奶头的视频| 国产精品一区二区在线不卡| 国产成人系列免费观看| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片 | 日韩免费高清中文字幕av| 久久精品人人爽人人爽视色| 国产aⅴ精品一区二区三区波| 国产真人三级小视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美视频一区| 国产亚洲欧美98| 国产日韩欧美亚洲二区| 久久这里只有精品19| 999久久久国产精品视频| 亚洲aⅴ乱码一区二区在线播放 | 国产精品二区激情视频| 高清av免费在线| 免费少妇av软件| 国产亚洲欧美98| 韩国精品一区二区三区| 精品久久蜜臀av无| 亚洲国产看品久久| 亚洲精品国产一区二区精华液| 午夜福利乱码中文字幕| 女人精品久久久久毛片| 精品国产超薄肉色丝袜足j| 免费女性裸体啪啪无遮挡网站| 99国产精品99久久久久| 天堂动漫精品| 欧美大码av| 国产蜜桃级精品一区二区三区 | 国产男靠女视频免费网站| 人人妻人人澡人人看| 女人被躁到高潮嗷嗷叫费观| 午夜精品久久久久久毛片777| 欧美+亚洲+日韩+国产| 午夜老司机福利片| 精品亚洲成a人片在线观看| 精品久久蜜臀av无| 99re在线观看精品视频| 成人国语在线视频| 男人舔女人的私密视频| 青草久久国产| 国产亚洲精品久久久久久毛片 | 中文字幕高清在线视频| 久久精品人人爽人人爽视色| 午夜日韩欧美国产| 色94色欧美一区二区| 日韩免费高清中文字幕av| 男女之事视频高清在线观看| 久热爱精品视频在线9| 麻豆国产av国片精品| 乱人伦中国视频| 亚洲一区高清亚洲精品| 波多野结衣av一区二区av| 大码成人一级视频| 在线观看一区二区三区激情| 免费在线观看黄色视频的| 国产欧美日韩一区二区三区在线| 免费观看a级毛片全部| 好男人电影高清在线观看| 人人妻人人澡人人看| 国产精品综合久久久久久久免费 | 午夜日韩欧美国产| 国产在视频线精品| 日韩欧美三级三区| 国产高清videossex| 在线永久观看黄色视频| 久久久久精品国产欧美久久久| 亚洲精品一二三| 亚洲情色 制服丝袜| 日本精品一区二区三区蜜桃| 亚洲精品久久午夜乱码| 老司机影院毛片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av第一区精品v没综合| 下体分泌物呈黄色| a在线观看视频网站| 久久久久久久国产电影| aaaaa片日本免费| av在线播放免费不卡| 两个人看的免费小视频| 深夜精品福利| 黑人欧美特级aaaaaa片| 捣出白浆h1v1| 一级毛片高清免费大全| videosex国产| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 亚洲午夜精品一区,二区,三区| 日韩人妻精品一区2区三区| 国产亚洲av高清不卡| av福利片在线| 午夜福利视频在线观看免费| 熟女少妇亚洲综合色aaa.| 国产精品免费一区二区三区在线 | 国产精品久久久久久精品古装| www.精华液| 精品国产一区二区三区久久久樱花| 黄色 视频免费看| √禁漫天堂资源中文www| 天堂动漫精品|