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

    基于VOF的弧底梯形短喉道量水槽水力性能研究

    2017-03-21 02:04:21馬孝義景志芳甘學(xué)濤
    中國農(nóng)村水利水電 2017年12期
    關(guān)鍵詞:喉道測流水頭

    張 敏,馬孝義,景志芳,柳 燁,甘學(xué)濤

    (西北農(nóng)林科技大學(xué) 水利與建筑工程學(xué)院,陜西 楊凌 712100)

    水作為生命之源,與人類活動、社會發(fā)展緊密相連。目前我國的農(nóng)業(yè)灌溉水利用系數(shù)平均不足0.5,節(jié)水空間很大。隨著經(jīng)濟(jì)的發(fā)展及水資源的日益緊張,為了促進(jìn)農(nóng)業(yè)節(jié)水、實(shí)現(xiàn)水資源的可持續(xù)發(fā)展,加強(qiáng)灌區(qū)量水工作,對灌區(qū)用水進(jìn)行有效的控制和精確的量測變得十分重要[1,2]?;〉滋菪吻烙捎诮Y(jié)構(gòu)簡單,整體性強(qiáng),水流條件好,過流順暢,底部的弧形能夠有效減少渠道淤積,降低工程維護(hù)費(fèi)用,且抗凍脹性好,因此適宜在北方渠道大規(guī)模應(yīng)用。針對弧底梯形渠道量水設(shè)備不足的現(xiàn)狀[3,4],選定結(jié)構(gòu)簡單、經(jīng)濟(jì)適用、測量精度高的量水設(shè)施尤為重要。

    R.L.Parshall[5]在文丘里槽[6]基礎(chǔ)上研制出巴歇爾量水槽,它由短直喉道、上游收縮段和下游擴(kuò)散段三部分組成,該槽具有壅水高度低,量水精度高,測流幅度大,水頭損失小,不易淤積等特點(diǎn),是近年來使用較廣泛的短喉道槽,但其因結(jié)構(gòu)復(fù)雜,會造成施工困難。Blaisdel[7]對Parshall的測流數(shù)據(jù)和流量公式進(jìn)行深入分析,得出其測流相對誤差在 5%以內(nèi)。李杰[8]等對自由和淹沒出流情況下量水槽的水位流量關(guān)系及沖沙情況等水力性能進(jìn)行了研究。李香平[9]等將巴歇爾量水槽斜槽底變?yōu)樗?,并對測流范圍進(jìn)行重新標(biāo)定。以某灌區(qū)灌溉渠系改造工程為原型,筆者在巴歇爾量水槽基礎(chǔ)上進(jìn)行改造,提出了與之配套的弧底梯形短喉道式量水槽。對其在不同底坡、喉道長度及收縮比條件下的水力性能進(jìn)行研究,并對其影響因素進(jìn)行分析,以期為新型量水槽的推廣應(yīng)用提供一定的參考。

    1 量水槽體型設(shè)計(jì)

    1.1 量水槽結(jié)構(gòu)

    弧底梯形短喉道量水槽由進(jìn)口收縮段、喉道段、出口擴(kuò)散段三部分組成,其中進(jìn)口收縮段呈3∶1收縮,出口擴(kuò)散段呈6∶1擴(kuò)散。通過調(diào)整喉口斷面弧底半徑R、喉口邊坡系數(shù)m,使水流產(chǎn)生局部收縮,并在收縮段產(chǎn)生臨界流,由穩(wěn)定的水位流量單值函數(shù)關(guān)系得到過槽流量,達(dá)到測流目的。收縮比ε=A0/A,A0為量水槽喉口斷面面積,A為渠道斷面面積。短喉道量水槽由于增加了喉道段長度使水流流態(tài)變化得到緩沖,能夠較為平穩(wěn)的銜接,能更精確確定上游水深斷面位置,且由于水流內(nèi)部摻混作用的加強(qiáng),使得流速分布更趨于均勻化分布。量水槽結(jié)構(gòu)示意圖如圖1。

    圖中:W為量水槽喉寬;a為喉長;L+a為量水槽長度;B為渠道頂寬;h為渠槽水位。圖1 弧底梯形渠道短喉道量水槽結(jié)構(gòu)示意圖Fig.1 Structure of arc-based trapezoidal canal and short-throat flow measurement flume

    1.2 量水槽尺寸設(shè)計(jì)

    本文對收縮比ε=0.61,0.65,0.7的量水槽在不同喉道長度a=0,0.25,0.5,1 m下弗勞德數(shù)、水頭損失、水面線等的變化進(jìn)行了模擬分析。渠道參數(shù)及數(shù)值試驗(yàn)方案如表1所示。

    表1 量水槽尺寸設(shè)計(jì)方案Tab.1 Dimension design scheme of arc-based trapezoidal and measuring flume

    1.3 設(shè)計(jì)流量

    弧底梯形渠道設(shè)計(jì)流量Q為設(shè)計(jì)水位下的渠道流量,可按下式計(jì)算:

    (1)

    式中:Q為渠道設(shè)計(jì)流量;A為渠道斷面面積,R為底弧半徑;C為謝才系數(shù);i為渠道底坡。

    2 量水槽數(shù)學(xué)模型

    2.1 控制方程

    弧底梯形短喉道量水槽測流為不可壓縮黏性流體運(yùn)動,根據(jù)基本物理守恒定律,過槽流量為牛頓流體,則測流方程為連續(xù)性方程和Navier-Stokes方程[10]。

    連續(xù)性方程:

    (2)

    式中:ρ為流體密度,本文中取1 000 kg/m3;t為時間;ui為沿i(i為笛卡兒坐標(biāo)中的x、y或z)方向的速度分量;xi為沿i方向的坐標(biāo)分量。

    Navier-Stokes方程:

    (3)

    式中:ui、uj分別為沿i和j方向速度分量;xi、xj分別為沿i和j方向的坐標(biāo)分量;p為壓力;μ為分子黏性系數(shù);μ′i、μ′j為各脈動流速分量。

    本文研究的弧底梯形渠道短喉道量水槽過槽水流屬于具有高雷諾數(shù)的湍流,因此選擇RNGk-ε模型進(jìn)行湍流數(shù)值模擬,使控制方程組封閉的同時因考慮湍流漩渦而有效改善了精度。湍流脈動的動能k和耗散率ε如下:

    湍動能k方程:

    (4)

    其中:

    式中:μ為動力黏滯系數(shù),N·s/m2;μt為湍動黏度,Cμ=0.084 5;μeff=μ+μt,N·s/m2;Gk為平均速度引起的紊動能產(chǎn)生項(xiàng),αk=1.39。

    湍動耗散率ε方程:

    (5)

    2.2 三維模型建立、網(wǎng)格劃分與邊界條件設(shè)置

    本研究用Matlab進(jìn)行量水槽選型,Autocad實(shí)現(xiàn)三維建模,聯(lián)合Gambit完成網(wǎng)格劃分。渠道上下游邊界長度變化會導(dǎo)致網(wǎng)格數(shù)發(fā)生變化,從而對計(jì)算時長和計(jì)算精度產(chǎn)生影響。通過試驗(yàn)可得,選取的上游渠段過短,難以消除水面跌落對測流精度的影響[4]。理論上得知,上下游渠道長度越長,計(jì)算精度就越高,但隨著計(jì)算長度的增加,計(jì)算域內(nèi)的網(wǎng)格單元數(shù)會大幅增加,其計(jì)算時間也就越長,效率也越低。因此從實(shí)際情況出發(fā),在兼顧精度和計(jì)算時長的關(guān)系的前提下,分別采用5種不同上下游渠道長度弧底梯形短喉道量水槽對計(jì)算域進(jìn)行網(wǎng)格劃分[11]。實(shí)驗(yàn)結(jié)果見表2,綜合考慮,數(shù)值模擬選取渠道長250 m左右,以量水槽為零點(diǎn),上游取200 m,下游取50 m,能夠得到更加逼真的模擬渠道流態(tài),減小上下游邊界對計(jì)算域的影響。

    表2 上下游不同邊界長度對計(jì)算影響Tab.2 The influence of different boundary lengths on the calculation of the upper and lower reaches

    建模時,坐標(biāo)原點(diǎn)為渠道進(jìn)口斷面圓弧中點(diǎn),取順?biāo)鞣较驗(yàn)閤軸(順?biāo)鳛檎?,寬度方向?yàn)閥軸(河道左岸為正),高度方向?yàn)閦軸(向上為正)。在量水槽段采用非結(jié)構(gòu)網(wǎng)格,并進(jìn)行局部加密,其高適應(yīng)性與量水槽段復(fù)雜流動域結(jié)構(gòu)更為匹配,渠道段采用結(jié)構(gòu)網(wǎng)格,因渠道左右對稱,可取一半計(jì)算域,網(wǎng)格總數(shù)約8.0×105,量水槽尺寸見表1,邊界條件設(shè)置如下:進(jìn)口邊界由空氣進(jìn)口和水流進(jìn)口二部分組成,空氣進(jìn)口邊界條件設(shè)置Pressure Inlet,水流進(jìn)口設(shè)置為Velocity Inlet,出口設(shè)置為Pressure Outlet,渠道斷面中心面取對稱面Symmetry,渠道的整個底部以及邊壁選擇固壁邊界條件wall,并設(shè)置無滑移選項(xiàng),采用標(biāo)準(zhǔn)的壁面函數(shù)法處理,模型網(wǎng)格劃分如圖2。

    圖2 弧底梯形短喉道量水槽模型與網(wǎng)格劃分Fig.2 Model and mesh division of arc-based trapezoidal canal and short-throat flow measurement flume

    3 量水槽計(jì)算方法與模擬結(jié)果分析

    3.1 計(jì)算方法

    本文采用RNGk-ε湍流模型,它將非穩(wěn)態(tài)Navier-Stokes方程對一個平衡態(tài)作Gauss統(tǒng)計(jì)展開,并用對脈動頻譜的波數(shù)段作濾波的方法,從理論上導(dǎo)出高Re數(shù)k-ε模型,所得出的k-ε方程形式上同標(biāo)準(zhǔn)k-ε模型完全一樣,但不同的是5個系數(shù)之值不是根據(jù)實(shí)驗(yàn)數(shù)據(jù)而是由理論分析得出,因而在計(jì)算復(fù)雜的流動方面,理論上精度要高于標(biāo)準(zhǔn)的k-ε模型[12]。本研究通過VOF方法跟蹤自由表面,水流從入口流入渠槽直至出口, 通過對體積分?jǐn)?shù)的迭代求解,就能夠自動生成水氣的交界面,設(shè)置監(jiān)視器中Mass Flow Rate選項(xiàng),當(dāng)其值接近于0且基本保持不變,認(rèn)為計(jì)算收斂。采用CFD軟件之Fluent6.3對建立好的三維模型進(jìn)行仿真計(jì)算[13]。

    3.2 流態(tài)分析與水面線變化

    流態(tài)分析是研究弧底梯形短喉道量水槽水力特性的基本要素。為觀測不同喉道長度、不同流量下量水槽段水面線沿程變化,選擇對底弧半徑為260的弧底梯形短喉道量水槽進(jìn)行模擬研究,i=1/1 800、ε=0.61,四種體型喉長分別為0,0.25,0.5,1 m。圖3為不同喉道長度量水槽過槽水流的水面線變化情況,由圖3可見,量水槽上游水面水流均平穩(wěn),在量水槽進(jìn)口收縮段(L=200~202.7 m)水面線開始緩慢降落,在出口擴(kuò)散段前半段處水面線急劇下降,從喉口處開始有所不同,喉長a=0量水槽擴(kuò)散段水流較其他3個湍急,水面曲率大,水面下降位置比其他量水槽靠近量水槽喉口,其他量水槽因喉道長度較其長給水面跌落以緩沖力,使水面跌落比較平緩。水流流經(jīng)量水槽受到側(cè)收縮而產(chǎn)生臨界流,從而形成單一穩(wěn)定的流量關(guān)系,此時通過上游水深即可得出過槽流量。

    圖3 不同喉道長度下量水槽過槽水流水面線Fig.3 The surface line of the flume under different throat lengths

    圖4為i=1/1 400、ε=0.61、a=1時,不同流量下弧底梯形短喉道量水槽水面線變化情況。由圖可知,水面線變化均平緩,水深沿水流方向變化趨勢與圖3一致,且水面線隨流量增大而上移。

    圖4 不同流量下水面線變化Fig.4 Variation of surface line at different flow rate

    3.3 弗勞德數(shù)變化及臨界水深位置

    3.3.1 不同情況下槽前弗勞德數(shù)變化

    弗勞德數(shù)作為明渠水流流態(tài)判別的標(biāo)準(zhǔn),在量水槽測流過程中是十分重要的。由圖5可以看出,渠道上游水流平穩(wěn),F(xiàn)r小于0.5,滿足明渠測流規(guī)范要求。以ε=0.61,i=1/1 400為例,不同喉道長度下槽前弗勞德數(shù)變化見圖5,a=1,i=1/1 400時不同收縮比情況下的槽前弗勞德數(shù)變化見圖6。由圖可得,隨流量增大,弗勞德數(shù)呈先減小后增大趨勢,總的來看,弗勞德數(shù)與喉道長度為負(fù)相關(guān)關(guān)系,即喉長越長,弗勞德數(shù)越小。隨著收縮比減小,弗勞德數(shù)逐漸減小,同一收縮比情況下,弗勞德數(shù)呈先減小后增大的趨勢。

    圖5 不同喉道長度下弗勞德數(shù)變化Fig.5 Variation of Froude number at different throat lengths

    圖6 不同收縮比下弗勞德數(shù)變化Fig.6 Variation of Froude number at different shrinkage ratios

    圖7 不同喉道長度沿程弗勞德數(shù)變化Fig.7 Variation of Froude number along different throat length

    3.3.2 沿程弗勞德數(shù)變化及臨界水深位置

    佛勞德數(shù)Fr可以作為明渠水流流態(tài)判別的標(biāo)準(zhǔn)。Fr<1時,水流流態(tài)為緩流;Fr=1時,水流流態(tài)為臨界流;Fr>1時,水流流態(tài)為急流。ε=0.65,i=1/1600時,量水槽在不同喉道長度情況下,沿程弗勞德數(shù)的變化情況見圖7,由圖可得,水流在量水槽入口至喉口收縮斷面為緩流,F(xiàn)r<1;在量水槽收縮段內(nèi),隨著過流斷面減小,流速逐漸增大,水面線開始下降,佛勞德數(shù)不斷增大,水流流態(tài)仍為緩流;在量水槽出口擴(kuò)散段內(nèi),水面線急劇下降,斷面佛勞德數(shù)繼續(xù)增大,當(dāng)水面線下降到一定程度時形成臨界流,對應(yīng)的水深即為臨界水深,此時Fr=1;由于慣性作用,水面線繼續(xù)降低,流速不斷增大,其大小和方向均發(fā)生變化,水流流態(tài)變?yōu)榧绷?,此時Fr>1;在量水槽下游段,F(xiàn)r值達(dá)到最大??梢钥闯觯?條曲線Fr變化總體趨勢一致。

    為準(zhǔn)確確定Fr=1即臨界水深位置,將弗勞德數(shù)相關(guān)數(shù)據(jù)導(dǎo)入matlab中,對左圖虛線框部分進(jìn)行局部放大,用ginput進(jìn)行點(diǎn)的確定和坐標(biāo)的讀取。李杰[14]所得巴歇爾量水槽臨界水深斷面在喉道段前5 cm內(nèi),肖苡辀[10]所得田間便攜式短喉道量水槽臨界水深在喉道段內(nèi)。本文所選四種不同體型量水槽喉道末端x值分別為202.7,202.95,203.2,203.7,matlab導(dǎo)出四個交點(diǎn)從左至右分別為(203.601,1),(203.655 3,1),(203.785 3,1),(204.018 1,1),可以看出,隨喉道長度的增加,臨界水深位置后移,且臨界水深斷面發(fā)生在量水槽擴(kuò)散段內(nèi),約為喉道段后30 cm左右。

    3.4 水頭損失分析

    液體的黏滯性是產(chǎn)生水頭損失的根本原因?;〉滋菪味毯淼懒克圻^槽水流水頭損失包括兩部分,因與渠槽邊壁摩擦產(chǎn)生的沿程水頭損失和因側(cè)收縮存在,渠道漸變段形狀改變,水流結(jié)構(gòu)急劇調(diào)整而造成的局部水頭損失。對于喉口寬度小于22.86 m的短喉槽,水頭損失較小,可以用上、下游測量水頭差近似代替[15]??偹^:h=Z+P/γ+v2/2g,分別取弧底梯形短喉道量水槽收縮段首斷面和擴(kuò)散段末斷面為上、下游水頭控制斷面,上下游水頭差即為水頭損失hw。

    以ε=0.61,i=1/2 400為例,表3為弧底梯形短喉道量水槽在時不同喉道長度下水頭損失的變化情況,由表4可知,收縮比一定時,隨喉道長度增大,水頭損失呈增大趨勢。表4為喉長a=0.25,i=1/1 400情況下,不同收縮比所對應(yīng)水頭損失變化情況。由表可得,喉長一定時,隨收縮比增大,水頭損失呈減小趨勢,這是因?yàn)槭湛s比越小,使得側(cè)收縮越大邊界條件變化急劇,局部水頭損失變大;同一收縮比時,量水槽水頭損失隨著流量增大而減小,個別點(diǎn)小流量時誤差較大,因?yàn)樾×髁肯率湛s比較大流量的收縮比小,導(dǎo)致小流量相對大流量水頭損失要大。

    3.5 公式以及模擬誤差

    試驗(yàn)?zāi)M了3種弧底梯形短喉道量水槽不同流量工況下的水位流量關(guān)系,水位流量相關(guān)性較好,通過數(shù)據(jù)分析可知,流量與收縮比ε、上游水深h、喉道長度a有關(guān)。對流量數(shù)據(jù)進(jìn)行回歸分析,得到流量的經(jīng)驗(yàn)公式如下:

    a=1 mQ=1.086 5(εw)1.113 9h2.020 41

    a=0.5 mQ=1.276 7(εw)0.902 2h2.005 5

    a=0.25 mQ=1.348 5(εw)0.88h2.010 4

    (6)

    表3 不同喉道長度下水頭損失變化Tab.3 Variation of head loss at different throat lengths

    表4 不同收縮比下水頭損失變化Tab.4 Variation of head loss at different shrinkage ratio

    式中:Q為過流流量,m3/s;ε為量水槽收縮比;h為上游水深,m;W為喉口寬度,m;該公式測流范圍為0.408 6~5.589 6 m3/s。

    由公式(7)可以看出,流量受收縮比ε,喉道寬度W和上游水深h影響較大,且流量指數(shù)近似為2,可以嘗試重新構(gòu)造短喉槽自由出流流量公式,使過槽流量與上游水深的平方成正比,通過對不同喉道長度時的實(shí)驗(yàn)數(shù)據(jù)回歸分析可得形如下的統(tǒng)一流量計(jì)算公式:

    Q=A(εw)Bh2

    (7)

    式中:A=-0.353 8a+1.443 6,B=0.792 0ε0.329 6 a,a為喉道長度。

    圖8 渠道流量與計(jì)算流量對比結(jié)果Fig.8 Comparison results of channel flow and calculated flow

    用公式(7)流量公式所得計(jì)算流量與渠道流量進(jìn)行比較,圖8表明二者具有較好的一致性,相關(guān)性可達(dá)0.998 3。表5為部分工況下渠道流量、模擬流量、公式計(jì)算流量對比結(jié)果,最大誤差為6.98%,平均誤差為1.74%,滿足灌區(qū)測流要求。

    表5 部分工況下渠道流量、模擬流量與公式計(jì)算流量的對比結(jié)果Tab.5 Results of channel flow, simulated flow and formula calculation under partial working conditions

    4 結(jié) 語

    本研究采用 RNGk-ε湍流模型和VOF方法相耦合,基于Fluent6.3大型流體力學(xué)計(jì)算軟件對弧底梯形短喉道量水槽在上游渠長200 m,下游渠長50 m情況下進(jìn)行數(shù)值模擬研究,研究結(jié)果如下:

    (1)通過對9種弧底梯形短喉道量水槽和3種無喉道量水槽在不同工況下槽前和沿程弗勞德數(shù)變化情況進(jìn)行分析,確定臨界水深發(fā)生位置在量水槽擴(kuò)散段,喉道段下游30 cm內(nèi)。

    (2)通過對弧底梯形短喉道量水槽水頭損失分析得知,最大水頭損失為總水頭的5.98%,最小水頭損失為總水頭的1.09%,小于長喉道量水槽的13%,且比無喉道量水槽水流平穩(wěn),水面變化穩(wěn)定性好,該理論為修建灌區(qū)新型量水槽提供了新思路。

    (3)模擬流量與渠道流量誤差除個別小流量點(diǎn)外均小于5%,二者吻合程度較好,說明模擬具有一定的可靠性。所建立的自由出流狀態(tài)下弧底梯形短喉道量水槽的統(tǒng)一流量公式平均誤差為2.43%,精度較高,簡單可用,以期為后續(xù)其他量水槽測流公式的擬定提供參考。

    [1] 嚴(yán)曉軍. 我國灌區(qū)量水技術(shù)現(xiàn)狀及發(fā)展趨勢[J]. 中國科技信息,2010,(17):89-90.

    [2] 郝晶晶,馬孝義,王波雷,等. 灌區(qū)量水設(shè)備的研究應(yīng)用現(xiàn)狀與發(fā)展趨勢[J]. 中國農(nóng)村水利水電,2008,(4):39-41.

    [3] 鐘新銘,馬孝義,王文娥,許坤鵬,樊琨.弧底梯形渠道無喉道量水槽水力特性的數(shù)值模擬[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,(4):218-224.

    [4] 楊珮珮.弧底梯形無喉道量水槽的數(shù)值模擬與體型優(yōu)化研究[D]. 陜西楊凌:西北農(nóng)林科技大學(xué),2016.

    [5] Parshall RL. The Improved Venturi Flume[M]. USA: Colorado Experiment Station, Fort Collins, Coloado Agricultural College, 1928.

    [6] Cone VM. The Venturi flume [J]. Journal of Agricultural Research, 1917, 9(4):116-129.

    [7] Blaisdell F W. Results of Parshall Flume Tests[J]. Journal of Irrigation and Drainage Engineering, 1994,120(2):278-291.

    [8] 李 杰, 呂宏興, 胥維纖. 巴歇爾量水槽水流特性試驗(yàn)研究[J]. 節(jié)水灌溉,2009,(12): 45-47.

    [9] 李香平, 楊兆選, 張 濤. 改進(jìn)型帕歇爾水槽[J]. 天津大學(xué)學(xué)報(bào), 2004,(5): 419-422.

    [10] 肖苡辀,王文娥,胡笑濤. 基于FLOW_3D的田間便攜式短喉槽水力性能數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(3): 55-61.

    [11] 潘志寶. 機(jī)翼形量水槽水力特性試驗(yàn)與數(shù)值模擬研究[D]. 陜西楊凌:西北農(nóng)林科技大學(xué), 2009.

    [12] 張明亮,沈永明.RNGk-ε湍流模型在三維彎曲河流中的應(yīng)用[J].水力發(fā)電學(xué)報(bào),2007,(5):86-91.

    [13] 鐘新銘,馬孝義,楊珮珮. 弧底梯形渠道無喉道量水槽水位流量關(guān)系數(shù)值模擬[J]. 水土保持研究,2016,(2):340-345.

    [14] 李 杰.巴歇爾量水槽水力特性試驗(yàn)研究[D]. 陜西楊凌:西北農(nóng)林科技大學(xué),2010.

    [15] 王長德. 量水技術(shù)與設(shè)施[M]. 北京: 中國水利水電出版社,2005.

    猜你喜歡
    喉道測流水頭
    渠道斷面自動測流系統(tǒng)在位山灌區(qū)測水量水中的應(yīng)用
    水文測流技術(shù)方法與進(jìn)展分析
    石河子科技(2022年4期)2022-03-24 05:45:28
    玉龍水電站機(jī)組額定水頭選擇設(shè)計(jì)
    小水電(2021年6期)2021-12-15 02:00:06
    泵房排水工程中剩余水頭的分析探討
    曹店灌區(qū)渠首測流存在的問題及對策
    山東水利(2018年6期)2018-03-24 13:00:35
    M9在建設(shè)在線雷達(dá)測流設(shè)備選址中的應(yīng)用
    U型渠道無喉道量水槽流動規(guī)律數(shù)值模擬
    勝利油田致密砂巖油藏微觀孔隙結(jié)構(gòu)特征
    亞聲速二喉道流場不對稱現(xiàn)象研究
    溪洛渡水電站機(jī)組運(yùn)行水頭處理
    极品人妻少妇av视频| 欧美乱妇无乱码| 热99re8久久精品国产| 亚洲午夜精品一区,二区,三区| 超碰成人久久| 女人高潮潮喷娇喘18禁视频| 欧美日韩中文字幕国产精品一区二区三区 | 人妻 亚洲 视频| 国产精品久久久人人做人人爽| 欧美午夜高清在线| 搡老岳熟女国产| 丰满饥渴人妻一区二区三| 少妇粗大呻吟视频| 99国产极品粉嫩在线观看| 巨乳人妻的诱惑在线观看| 亚洲精品一二三| 老司机影院毛片| 99热国产这里只有精品6| 黄色丝袜av网址大全| 激情在线观看视频在线高清 | 午夜久久久在线观看| 超碰成人久久| av视频免费观看在线观看| 亚洲av日韩在线播放| 女人被狂操c到高潮| 亚洲av成人一区二区三| www.自偷自拍.com| 视频区欧美日本亚洲| 一区在线观看完整版| 搡老乐熟女国产| 亚洲专区国产一区二区| 午夜福利,免费看| 国产成人影院久久av| √禁漫天堂资源中文www| 欧美性长视频在线观看| 黄色女人牲交| 久久午夜亚洲精品久久| 男男h啪啪无遮挡| 精品电影一区二区在线| 精品国产一区二区三区久久久樱花| 热99re8久久精品国产| 黄色毛片三级朝国网站| 日本五十路高清| 久久国产精品大桥未久av| 无人区码免费观看不卡| 99久久99久久久精品蜜桃| 91字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av成人不卡在线观看播放网| 黑丝袜美女国产一区| 一进一出抽搐动态| 国产激情欧美一区二区| 国产精品综合久久久久久久免费 | ponron亚洲| 成年人免费黄色播放视频| 亚洲熟女精品中文字幕| 淫妇啪啪啪对白视频| 动漫黄色视频在线观看| 国产三级黄色录像| 日韩有码中文字幕| 丁香六月欧美| 女人高潮潮喷娇喘18禁视频| 亚洲av欧美aⅴ国产| 欧美另类亚洲清纯唯美| 精品国产超薄肉色丝袜足j| 欧美黄色片欧美黄色片| 免费在线观看黄色视频的| 一级毛片精品| 亚洲少妇的诱惑av| 免费女性裸体啪啪无遮挡网站| 在线av久久热| 欧美精品啪啪一区二区三区| 国产人伦9x9x在线观看| 中文字幕制服av| 久久精品人人爽人人爽视色| 十分钟在线观看高清视频www| 一本大道久久a久久精品| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 欧美一级毛片孕妇| 麻豆成人av在线观看| 亚洲,欧美精品.| 中文欧美无线码| 免费观看精品视频网站| 男女午夜视频在线观看| 最近最新免费中文字幕在线| 窝窝影院91人妻| 18禁观看日本| 亚洲欧美一区二区三区黑人| 91大片在线观看| 亚洲精品国产区一区二| 黄色视频,在线免费观看| 日本vs欧美在线观看视频| 日韩免费高清中文字幕av| 欧美av亚洲av综合av国产av| 亚洲av熟女| www.精华液| 亚洲欧美激情在线| 亚洲,欧美精品.| 成人亚洲精品一区在线观看| 欧美乱码精品一区二区三区| 久久国产精品影院| 欧美日韩亚洲高清精品| 巨乳人妻的诱惑在线观看| 国产一区二区激情短视频| 香蕉丝袜av| 亚洲人成电影观看| 99久久国产精品久久久| 欧美大码av| 成年人午夜在线观看视频| 成人国产一区最新在线观看| 国产91精品成人一区二区三区| 国产有黄有色有爽视频| 久9热在线精品视频| 国内久久婷婷六月综合欲色啪| 热re99久久国产66热| 亚洲熟女毛片儿| 国产xxxxx性猛交| 老汉色av国产亚洲站长工具| 亚洲一区二区三区欧美精品| 夜夜夜夜夜久久久久| 精品国产国语对白av| 免费少妇av软件| av超薄肉色丝袜交足视频| 最近最新中文字幕大全电影3 | 大香蕉久久网| 操美女的视频在线观看| 下体分泌物呈黄色| 亚洲精品粉嫩美女一区| 自拍欧美九色日韩亚洲蝌蚪91| 精品高清国产在线一区| 老熟妇仑乱视频hdxx| 国产成人免费无遮挡视频| 熟女少妇亚洲综合色aaa.| 色播在线永久视频| 建设人人有责人人尽责人人享有的| 日韩制服丝袜自拍偷拍| 午夜免费观看网址| 亚洲av欧美aⅴ国产| 久久久国产成人精品二区 | 国产片内射在线| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 久久香蕉国产精品| 如日韩欧美国产精品一区二区三区| 色尼玛亚洲综合影院| 亚洲欧美激情在线| 黄频高清免费视频| 淫妇啪啪啪对白视频| 亚洲精品在线观看二区| 在线观看午夜福利视频| av国产精品久久久久影院| 极品少妇高潮喷水抽搐| 大香蕉久久成人网| 欧美 亚洲 国产 日韩一| 日韩欧美一区视频在线观看| 无人区码免费观看不卡| 黑人巨大精品欧美一区二区mp4| 黄色视频不卡| 国产成人啪精品午夜网站| 国产成人欧美| av网站免费在线观看视频| 高潮久久久久久久久久久不卡| 午夜免费成人在线视频| 在线观看免费日韩欧美大片| 国产免费男女视频| 欧美日本中文国产一区发布| 亚洲精品国产一区二区精华液| 超碰97精品在线观看| 俄罗斯特黄特色一大片| 黄色a级毛片大全视频| 免费看a级黄色片| 欧美精品人与动牲交sv欧美| 一级毛片高清免费大全| 18禁国产床啪视频网站| 黄色片一级片一级黄色片| 色婷婷久久久亚洲欧美| 一级a爱片免费观看的视频| 18禁美女被吸乳视频| 国产精品亚洲av一区麻豆| 国产精品美女特级片免费视频播放器 | 新久久久久国产一级毛片| 99re6热这里在线精品视频| 母亲3免费完整高清在线观看| 精品久久蜜臀av无| 免费一级毛片在线播放高清视频 | 成熟少妇高潮喷水视频| 高清欧美精品videossex| 精品人妻1区二区| 一个人免费在线观看的高清视频| 欧美日韩一级在线毛片| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| 午夜福利乱码中文字幕| videosex国产| 久久精品成人免费网站| 亚洲成国产人片在线观看| netflix在线观看网站| 村上凉子中文字幕在线| 欧美精品一区二区免费开放| 亚洲一区中文字幕在线| 99国产综合亚洲精品| www.熟女人妻精品国产| 天天躁夜夜躁狠狠躁躁| 嫁个100分男人电影在线观看| 亚洲午夜理论影院| 亚洲熟妇中文字幕五十中出 | 好看av亚洲va欧美ⅴa在| 两个人免费观看高清视频| 日韩有码中文字幕| 在线观看免费视频日本深夜| 亚洲一码二码三码区别大吗| 欧美中文综合在线视频| 欧美日韩瑟瑟在线播放| 免费久久久久久久精品成人欧美视频| 韩国av一区二区三区四区| 国产精品乱码一区二三区的特点 | 久久亚洲真实| 在线播放国产精品三级| 丝袜在线中文字幕| 亚洲av日韩精品久久久久久密| 国产av又大| 一进一出抽搐动态| 黄色怎么调成土黄色| 在线观看免费视频日本深夜| 国产欧美日韩一区二区三| 久久久久国内视频| av一本久久久久| 午夜两性在线视频| 亚洲av美国av| 欧美日韩一级在线毛片| 成人精品一区二区免费| 三上悠亚av全集在线观看| 大香蕉久久网| 18禁黄网站禁片午夜丰满| 一区二区三区国产精品乱码| 丝瓜视频免费看黄片| 午夜亚洲福利在线播放| 大片电影免费在线观看免费| 免费看a级黄色片| 欧美乱妇无乱码| 又黄又爽又免费观看的视频| 欧美激情高清一区二区三区| 男男h啪啪无遮挡| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 国产成人欧美| 国产精品成人在线| 女同久久另类99精品国产91| 国产精品久久久久久人妻精品电影| 中文字幕av电影在线播放| 日韩欧美一区视频在线观看| 操美女的视频在线观看| www.自偷自拍.com| 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 国产1区2区3区精品| 国产av一区二区精品久久| 两人在一起打扑克的视频| 精品少妇一区二区三区视频日本电影| 91精品三级在线观看| 中国美女看黄片| 亚洲熟妇熟女久久| 这个男人来自地球电影免费观看| 久久久久国产一级毛片高清牌| 曰老女人黄片| 91麻豆精品激情在线观看国产 | 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| 国产高清视频在线播放一区| 美女高潮到喷水免费观看| 亚洲成a人片在线一区二区| 欧美激情久久久久久爽电影 | 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看免费视频日本深夜| 亚洲av第一区精品v没综合| 少妇猛男粗大的猛烈进出视频| 国产91精品成人一区二区三区| 精品午夜福利视频在线观看一区| 在线视频色国产色| 最新在线观看一区二区三区| 极品少妇高潮喷水抽搐| 亚洲欧美激情综合另类| 91字幕亚洲| 两性夫妻黄色片| 美女国产高潮福利片在线看| 久久人妻熟女aⅴ| 在线国产一区二区在线| 女性生殖器流出的白浆| 欧美最黄视频在线播放免费 | 精品久久久久久久毛片微露脸| 精品久久久久久久久久免费视频 | 美女高潮到喷水免费观看| 久久人妻av系列| 久久久精品免费免费高清| 久久久久久亚洲精品国产蜜桃av| 男女午夜视频在线观看| 免费不卡黄色视频| 午夜精品国产一区二区电影| 一级作爱视频免费观看| 精品国产亚洲在线| 一级,二级,三级黄色视频| 国产精品乱码一区二三区的特点 | 亚洲精品粉嫩美女一区| 两个人看的免费小视频| 亚洲第一青青草原| 精品第一国产精品| 99热只有精品国产| 国产日韩一区二区三区精品不卡| 欧美老熟妇乱子伦牲交| 精品第一国产精品| 亚洲av美国av| 国产成人av激情在线播放| 少妇被粗大的猛进出69影院| 黑人猛操日本美女一级片| 亚洲欧美色中文字幕在线| 精品国产一区二区久久| 男人舔女人的私密视频| bbb黄色大片| 久久婷婷成人综合色麻豆| 又大又爽又粗| 中文亚洲av片在线观看爽 | 三上悠亚av全集在线观看| 两性夫妻黄色片| 亚洲熟女毛片儿| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| 久久亚洲真实| 侵犯人妻中文字幕一二三四区| 交换朋友夫妻互换小说| 美女国产高潮福利片在线看| 母亲3免费完整高清在线观看| 在线十欧美十亚洲十日本专区| 成在线人永久免费视频| 国产在视频线精品| 久久人妻熟女aⅴ| 国产成+人综合+亚洲专区| 大陆偷拍与自拍| 午夜日韩欧美国产| 一区二区三区国产精品乱码| 少妇被粗大的猛进出69影院| 国产精品一区二区在线不卡| 欧美激情 高清一区二区三区| 黑人猛操日本美女一级片| avwww免费| 一级片'在线观看视频| 国产精品久久久av美女十八| 欧美一级毛片孕妇| av在线播放免费不卡| 18禁裸乳无遮挡动漫免费视频| 国产精品影院久久| av片东京热男人的天堂| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 伊人久久大香线蕉亚洲五| 国产欧美日韩综合在线一区二区| 好男人电影高清在线观看| 欧美日韩一级在线毛片| 日韩欧美国产一区二区入口| 69精品国产乱码久久久| 国产欧美日韩一区二区三区在线| 久久久精品国产亚洲av高清涩受| 涩涩av久久男人的天堂| 狂野欧美激情性xxxx| 日韩成人在线观看一区二区三区| 亚洲综合色网址| 亚洲欧美色中文字幕在线| 少妇的丰满在线观看| 无遮挡黄片免费观看| 女人被狂操c到高潮| 国产亚洲精品第一综合不卡| 婷婷丁香在线五月| 黄色女人牲交| 人人澡人人妻人| 日本精品一区二区三区蜜桃| 免费观看精品视频网站| 90打野战视频偷拍视频| 国产成+人综合+亚洲专区| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久成人av| 男女之事视频高清在线观看| 又黄又爽又免费观看的视频| 国产黄色免费在线视频| 亚洲精品在线美女| 日本一区二区免费在线视频| 日本五十路高清| 国产成人免费观看mmmm| 天天影视国产精品| 男女高潮啪啪啪动态图| 欧美激情 高清一区二区三区| 91老司机精品| 美女高潮到喷水免费观看| 中文字幕高清在线视频| 免费女性裸体啪啪无遮挡网站| 欧美精品啪啪一区二区三区| 中国美女看黄片| 日本黄色视频三级网站网址 | 日韩成人在线观看一区二区三区| 国产精品久久电影中文字幕 | 婷婷丁香在线五月| 欧美午夜高清在线| 精品国内亚洲2022精品成人 | 久久午夜综合久久蜜桃| 日本黄色日本黄色录像| 王馨瑶露胸无遮挡在线观看| 男女之事视频高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 18禁裸乳无遮挡免费网站照片 | 人成视频在线观看免费观看| 亚洲三区欧美一区| 欧美日韩一级在线毛片| 国产亚洲精品第一综合不卡| 亚洲avbb在线观看| 亚洲男人天堂网一区| 久久青草综合色| 国产一区二区激情短视频| 91在线观看av| 天天操日日干夜夜撸| 女同久久另类99精品国产91| 69精品国产乱码久久久| 午夜影院日韩av| av网站在线播放免费| 少妇被粗大的猛进出69影院| xxxhd国产人妻xxx| 国产av精品麻豆| 欧美精品啪啪一区二区三区| 丝袜人妻中文字幕| 悠悠久久av| 国产精品秋霞免费鲁丝片| 午夜福利视频在线观看免费| 国产精品 欧美亚洲| 国产一卡二卡三卡精品| 黄色怎么调成土黄色| 90打野战视频偷拍视频| 50天的宝宝边吃奶边哭怎么回事| 欧美激情 高清一区二区三区| www.精华液| 丁香欧美五月| 国产精品成人在线| 天天躁日日躁夜夜躁夜夜| 精品亚洲成国产av| 亚洲五月婷婷丁香| 9热在线视频观看99| 午夜福利在线免费观看网站| 我的亚洲天堂| 国产欧美日韩一区二区三| 热re99久久国产66热| 国产一区二区三区综合在线观看| 制服人妻中文乱码| 亚洲,欧美精品.| 国产亚洲欧美精品永久| 99re在线观看精品视频| 狂野欧美激情性xxxx| 亚洲第一av免费看| 狂野欧美激情性xxxx| 亚洲第一av免费看| 欧美成人免费av一区二区三区 | 18在线观看网站| 国产精品1区2区在线观看. | 国产人伦9x9x在线观看| 精品国产亚洲在线| 免费女性裸体啪啪无遮挡网站| 国产1区2区3区精品| 两个人免费观看高清视频| 久久久国产一区二区| 精品久久久久久,| 日本vs欧美在线观看视频| 好男人电影高清在线观看| 制服诱惑二区| 日本黄色视频三级网站网址 | 九色亚洲精品在线播放| 99国产精品一区二区蜜桃av | 大陆偷拍与自拍| 丝袜美足系列| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 日韩欧美三级三区| 在线播放国产精品三级| 一区二区日韩欧美中文字幕| 久久久久视频综合| 国产成人系列免费观看| 夜夜夜夜夜久久久久| 女人高潮潮喷娇喘18禁视频| 免费不卡黄色视频| 国产日韩欧美亚洲二区| av一本久久久久| 大香蕉久久网| 国产精品国产高清国产av | 18禁观看日本| 熟女少妇亚洲综合色aaa.| 一本一本久久a久久精品综合妖精| 精品高清国产在线一区| 欧美精品人与动牲交sv欧美| 一进一出抽搐gif免费好疼 | 国产乱人伦免费视频| 成人手机av| 99久久人妻综合| 欧美中文综合在线视频| 久久香蕉国产精品| 亚洲午夜精品一区,二区,三区| 欧美日韩视频精品一区| 亚洲成人免费av在线播放| 国产人伦9x9x在线观看| 国产男女超爽视频在线观看| 亚洲五月色婷婷综合| 国产成人av教育| 一区福利在线观看| 男女午夜视频在线观看| 国产精品一区二区在线不卡| 黑丝袜美女国产一区| 少妇 在线观看| 三级毛片av免费| 国产精品香港三级国产av潘金莲| av国产精品久久久久影院| 黄频高清免费视频| 身体一侧抽搐| 老司机午夜福利在线观看视频| 一本一本久久a久久精品综合妖精| 欧美乱码精品一区二区三区| 99riav亚洲国产免费| 午夜免费成人在线视频| 制服人妻中文乱码| 久久人妻福利社区极品人妻图片| 日本黄色视频三级网站网址 | 亚洲视频免费观看视频| 亚洲av片天天在线观看| 久久精品91无色码中文字幕| 老熟女久久久| 久久ye,这里只有精品| 午夜福利在线免费观看网站| 精品少妇一区二区三区视频日本电影| 极品人妻少妇av视频| 热99久久久久精品小说推荐| 91九色精品人成在线观看| 十八禁人妻一区二区| 精品亚洲成a人片在线观看| 成人影院久久| 91av网站免费观看| 亚洲精品在线美女| 欧美一级毛片孕妇| 黄色视频不卡| 国产精品99久久99久久久不卡| 黄色片一级片一级黄色片| 一区二区日韩欧美中文字幕| 亚洲成人手机| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩综合在线一区二区| 亚洲精品久久成人aⅴ小说| 在线观看免费高清a一片| a级毛片在线看网站| 成人国产一区最新在线观看| 狂野欧美激情性xxxx| 欧美黑人精品巨大| av天堂久久9| 两个人免费观看高清视频| 久久精品人人爽人人爽视色| 少妇猛男粗大的猛烈进出视频| 久久精品人人爽人人爽视色| 91在线观看av| 欧美日韩乱码在线| 欧美日韩瑟瑟在线播放| 男女下面插进去视频免费观看| 免费在线观看黄色视频的| 成在线人永久免费视频| 亚洲欧洲精品一区二区精品久久久| 黄片小视频在线播放| xxxhd国产人妻xxx| 老司机影院毛片| 色婷婷久久久亚洲欧美| 一级片免费观看大全| 又紧又爽又黄一区二区| 久久精品亚洲av国产电影网| 纯流量卡能插随身wifi吗| 亚洲成a人片在线一区二区| 国产真人三级小视频在线观看| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 日日摸夜夜添夜夜添小说| 大香蕉久久网| 国产精品久久久久成人av| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 日韩有码中文字幕| 在线视频色国产色| 男女高潮啪啪啪动态图| 十八禁网站免费在线| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 女同久久另类99精品国产91| 亚洲人成电影免费在线| 波多野结衣av一区二区av| 法律面前人人平等表现在哪些方面| 国产一卡二卡三卡精品| 国产成人免费无遮挡视频| 午夜亚洲福利在线播放| 亚洲欧美精品综合一区二区三区| 久久精品亚洲av国产电影网| 嫩草影视91久久| 中出人妻视频一区二区| 18禁国产床啪视频网站| 国产成人精品久久二区二区免费| 男女免费视频国产| 激情视频va一区二区三区| 热99re8久久精品国产| 搡老岳熟女国产| 久久人人爽av亚洲精品天堂| 久久国产精品大桥未久av| 无限看片的www在线观看| 十八禁人妻一区二区| av免费在线观看网站| 成人黄色视频免费在线看| 精品亚洲成国产av| 99久久综合精品五月天人人| 女人爽到高潮嗷嗷叫在线视频|