馬光飛，方 勇，吳燕明，程文韜，李 桃，李 超

?

不同進口高度排水管道用渦流裝置截流特性模擬與試驗

馬光飛，方 勇※，吳燕明，程文韜，李 桃，李 超

（1.水利部杭州機械設(shè)計研究所（水利部產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所），杭州 310024；2. 水利機械及其再制造技術(shù)浙江省工程實驗室，杭州 310024）

為研究不同進口高度的渦流截流裝置的截流特性，選用瞬態(tài)不可壓縮流動的N-S方程和重整化群 RNG-湍流模型，采用歐拉模型中的VOF（volume of fluid）法對渦流截流裝置進行了數(shù)值模擬，分析了進口高度對渦流截流裝置截流特性、截流效率的影響，并且研究了渦流裝置內(nèi)部流場分布特性，同時對渦流截流裝置的截流特性進行了試驗驗證。結(jié)果表明，進口高度越小，截流能力越強；渦流截流裝置的截流效率，隨著高度增大先減小后增大，轉(zhuǎn)折點在3/5~4/5之間；當(dāng)進口寬度固定時，為了保證不發(fā)生溢流或者堵塞，且能夠?qū)崿F(xiàn)較好截流，渦流截流裝置的進口高度應(yīng)選擇4/5~5/5之間，此時，渦流裝置的進口截面面積大于出口截面面積。渦流截流裝置的截流主要是因為在裝置內(nèi)部形成了局部低壓，靜壓力轉(zhuǎn)換為流體運動的動壓力，使流體產(chǎn)生高速旋流，發(fā)生了旋流截流。該研究可為城市污水截流裝置提供新的可選設(shè)計方法，降低污水溢流或者內(nèi)澇發(fā)生率。

渦流；水利設(shè)備；進口高度；瞬態(tài)數(shù)值模擬；截流特性；內(nèi)部流動特性；試驗驗證

0 引 言

隨著社會的進步，工業(yè)企業(yè)排放的污染性氣體、汽車尾氣等會形成大氣干沉降，累積在地表上，沉積在管道中，隨著降雨，這些污染物會隨著雨水的沖刷進入水體對水質(zhì)產(chǎn)生污染[1-3]，使得城市雨水管道污染問題變得更加嚴(yán)重。因此，如何控制城市強降雨雨水污染問題及加大雨水利用率已成為國內(nèi)外學(xué)者廣泛關(guān)注的研究熱點[4-7]。國內(nèi)外普遍采用的處理方式是：規(guī)劃、設(shè)計和實施截流裝置及其系統(tǒng)，如截流堰、限流封板、止回閥、控制閘等傳統(tǒng)截流裝置，但是上述截流裝置都存在著不同程度的缺陷[8-10]：如截流量調(diào)節(jié)不精確，截流需要改變原有截流井結(jié)構(gòu)，截流管中污水量及充滿度的動態(tài)變化需要實時觀察、截污流量小，且水力學(xué)設(shè)計使其實際截流實施困難等問題突出。

目前，研究水利裝置及其他流體流動現(xiàn)象的方法主要有試驗研究、理論分析和各種數(shù)值模擬軟件計算分析。焦愛萍等[11]用計算流體動力學(xué)（computational fluid dynamics，CFD）的方法，利用Fluent計算軟件，研究了不同流量比射流對水墊塘的流動特征、動水墊效應(yīng)和水墊塘底板時均沖擊動壓分布以及消能效果的影響。崔寶玲等[12]分析起旋器入射角度對旋進旋渦流量計性能的影響；趙偉國等[13]基于遺傳算法、人工神經(jīng)網(wǎng)絡(luò)和三維CFD流場數(shù)值模擬相結(jié)合的方法得到了最優(yōu)葉輪模型。程正飛等[14]采用耦合VOF法的三維非穩(wěn)態(tài)水氣兩相流滲流模型，實現(xiàn)了對碾壓混凝土壩復(fù)雜自由滲流場的數(shù)值模擬。季磊磊等[15]基于ANSYS CFX軟件分析了流量工況對混流泵內(nèi)壓力分布流線分布、湍動能的影響，探討了不同流量工況下不同監(jiān)測點的壓力脈動時域和頻域響應(yīng)。吳波等[16]基于CFD數(shù)值仿真結(jié)合單因素試驗和正交試驗，研究螺旋角、槽深、槽數(shù)、槽寬比和槽長壩長比對螺旋槽干氣密封性能的影響。張勤昭[17]以試驗結(jié)果為基礎(chǔ)驗證湍流模型，選擇SST模型作為湍流模型，建立獲得三偏心蝶閥詳細流場的數(shù)值模型，分析開度對三偏心蝶閥全流阻系數(shù)的影響。張義龍等[18]利用CFD數(shù)值模擬技術(shù)進行數(shù)值模擬，分析不同布水形式與流速的關(guān)系以及對傳質(zhì)效果的影響。朱榮生[19]等應(yīng)用ANSYS-CFX 軟件，對離心泵內(nèi)部發(fā)生空化時進行定常和非定常計算，采用泵內(nèi)部不同流道的速度場、壓力場、氣體體積分?jǐn)?shù)以及壓力脈動情況的計算結(jié)果分析流道內(nèi)部瞬態(tài)非線性的流動變化規(guī)律。針對類似的數(shù)值模擬和試驗研究方法還有很多，如文獻[20-24]，均采用數(shù)值模擬的方法分析了研究對象的內(nèi)部流動特性。綜上所述，采用數(shù)值模擬的方法分析水力學(xué)問題是可行的。

本文采用FLUENT計算軟件，選用瞬態(tài)數(shù)值模擬N-S方程和RNG-紊流模型和VOF（volume of fluid）法，對渦流截流裝置的內(nèi)部流場進行數(shù)值模擬，同時進行了試驗驗證。詳細研究了不同進口高度渦流截流裝置的截流特性、內(nèi)部流場的靜壓力場、速度場分布，揭示其截流機理，為城市污水截流裝置提供新的可選設(shè)計，降低污水溢流或者內(nèi)澇的發(fā)生。

1 模型建立及計算方法

1.1 物理模型

渦流截流裝置主要應(yīng)用于排水系統(tǒng)，其工作原理為當(dāng)流體流過渦流截流裝置時，在渦流截流裝置中心產(chǎn)生偏心的螺旋氣柱（進而產(chǎn)生局部低壓，甚至負(fù)壓），改變出水口過流斷面，從而限制出水流量，產(chǎn)生了截流效應(yīng)。

本研究的渦流截流裝置，主要包括前后蓋板，蝸殼型腔體，進出口，其模型示意圖如圖1a、1b所示。為了充分研究渦流截流裝置的截流特性和內(nèi)部流動特性，將其整體布施在上游檢查井內(nèi)，在進口寬度固定的前提下對多組進口高度下的模型進行研究。模型基本參數(shù)為進口高度=230 mm，蝸殼結(jié)構(gòu)直徑=460 mm，寬度=150 mm。研究對象模型參數(shù)見表1。

注：L為渦流裝置的寬度，mm；h為渦流裝置進口高度，mm；D為基圓直徑，mm。

1.2 數(shù)學(xué)模型和計算方法

研究對象物理模型的數(shù)值模擬方法，主要采用流體動力學(xué)計算軟件FLUENT軟件[25-28]進行計算，由于渦流截流裝置主要應(yīng)用于集水井中，其內(nèi)部流動并不是簡單的單相流動，而是多相組分流動。為了對其的截流特性和內(nèi)部流動特性進行精確預(yù)測，減少試驗成本，在本文的數(shù)值計算中將其視為氣液兩相流動。也正是由于氣相的存在，隨著集水井中水位的高低變化，在渦流截流裝置出口中心位置形成不同直徑大小的截流氣柱，才使得研究對象擁有良好的截流特性而不需要再提供外部動力能源，即可達到排水管網(wǎng)系統(tǒng)中的流量調(diào)度，減緩下游集水井中的溢流發(fā)生率。

基于流體質(zhì)量守恒方程、動量守恒方程、流體體積函數(shù)模型（VOF model）的多相流動方程及組分湍流方程（RNG-）進行三維湍流計算。

1）質(zhì)量守恒方程

質(zhì)量守恒方程又稱連續(xù)性方程，可表述為單位時間內(nèi)流體微元體中質(zhì)量的增加，等于同一時間間隔內(nèi)流入該微元體的凈質(zhì)量

連續(xù)性方程的適用范圍沒有限制，無論是可壓縮的或者不可壓縮的流體，黏性或者無粘性流體，瞬態(tài)或者穩(wěn)態(tài)流動均適用。

2）動量守恒方程

動量守恒方程是任何流動系統(tǒng)都必須滿足的基本定律。動量守恒定律實際上是牛頓第二定律，可表述為：任何控制微元中流體動量對時間的變化率等于外界作用在微元上各種力之和，用數(shù)學(xué)式表示為

3）流體體積函數(shù)方程（VOF model）

在VOF模型中，認(rèn)為水和氣共有相同的速度和壓力場，因而對水氣兩相流可以像單相流那樣采用一組方程來描述流場。該方法采用1個流體體積分?jǐn)?shù)來描述自由表面的各種變化，在控制體內(nèi)對第相流體的容積分?jǐn)?shù)規(guī)定為：α=0表示控制體內(nèi)無相流體，α=1表示控制體內(nèi)完全充滿相流體，0<α<1表示控制體內(nèi)部相流體的體積分?jǐn)?shù)，對所有流體的容積分?jǐn)?shù)總和為1，即∑α=1。

采用VOF方法追蹤流體自由表面對第相流體的控制微分方程為

流體自由界面的跟蹤通過求解該方程來完成，其具體位置采用幾何重建格式來確定。它采用分段線性近似的方法來表示自由水面線。在每個單元中，水氣交界面是具有不變斜率的斜線段，并用此線性分界面形狀來計算通過單元面上的流體通量。根據(jù)每個單元的容積分?jǐn)?shù)值及其偏微分，線性的水氣交界面相對于每個部分充滿的單元中心的位置就可以計算出來，從而確定其具體位置。

4）組分湍流方程（RNG-）

RNG-模型來源于嚴(yán)格的統(tǒng)計技術(shù)。它和標(biāo)準(zhǔn)-模型很相似，但進行了改進，對于具有旋流產(chǎn)生的結(jié)構(gòu)預(yù)測精度更加準(zhǔn)確，其方程式如下

湍動能方程

湍動能耗散率方程

在上述方程（1）－（5）中，為流體密度，kg/m3；為時間，s；=1，2，3，即x=、、，u=、、，、、分別為笛卡爾坐標(biāo)系中的三坐標(biāo)的速度分量，即、、方向上速度矢量的分量；為求和下標(biāo)；方程中通用模型常數(shù)，C=1.44，C2=1.92，σ=1.0，σ=1.3。

G為平均速度引起的湍動能產(chǎn)生項

方程中通用模型常數(shù)C=0.09。

2 網(wǎng)格劃分與邊界條件

由于研究對象具有多種結(jié)構(gòu)，即不同進口高度值，且其主要應(yīng)用于截流井（或集水井）中進行截流，且不同進口高度值，流道結(jié)構(gòu)網(wǎng)格劃分原理相同。因此，這里只給出全值安裝于截流井中的單開口流道網(wǎng)格劃分結(jié)構(gòu)。由于其中存在蝸殼型結(jié)構(gòu)，選擇非結(jié)構(gòu)化網(wǎng)格進行劃分，同時給出上下游延長管道（分別為進出口的5倍直徑和出口10倍直徑）的網(wǎng)格，這樣在計算中湍流流動可得到充分發(fā)展，見圖2。

注：D1為流體進口直徑，D1=300 mm；D2為流體出進口直徑，D2=200 mm。

邊界條件：流體進口采用速度進口，0.7 m/s（該流速為管道設(shè)計最小流速[29]），進口全部為介質(zhì)水，密度為998.2 kg/m3，粘度為1.005×10-3Pa·s；流體出口采用壓力出口，直接通大氣；空氣密度為1.205 kg/m3，粘度為1.81×10-5Pa·s；單開口入口處采用交界面interface進行定義，數(shù)值計算時間步長0.01 s，計算收斂殘差為10-6，速度壓力耦合采用PISO算法。

3 結(jié)果與分析

本文研究了渦流截流裝置的截流特性，使其有自適應(yīng)截流效果，通過截流裝置使集水井的出口流量隨著集水井中的水位發(fā)生變化，有效控制集水井的出口流量。在遇到強降雨時，使上游高地勢的水流流動減慢，降低下游集水井中的匯流量，使其在還未溢流前，擁有更多的離開時間，有效地為城市污水截流裝置提供新的可選設(shè)計。

3.1 渦流截流裝置的截流特性

圖3為不同進口高度時，渦流截流裝置的截流特性曲線。渦流裝置的截面為×的矩形截面，由于寬度=150 mm和初始進流流速=0.7 m/s，均為一定值，因此文中以渦流裝置進口高度為模型變化參數(shù)。其中橫坐標(biāo)為集水井出口的質(zhì)量流量，kg/s；縱坐標(biāo)為集水井中的水位，m。從圖中可以看出，1）在所有進口高度下，渦流截流裝置出口流量與水位變化關(guān)系均產(chǎn)生了S型轉(zhuǎn)折，起到了截流作用。2）進口高度越小，發(fā)生截流時的出口流量越小，最終的截流效果越明顯，即11>12>13>14>15。3）在進口高度1/5、2/5、3/5、4/5時，S型特性曲線的幅度相對較小，而在5/5時，S型特性曲線的幅度相對較大，這主要是因為在進口高度較小時，渦流截流裝置進口截面的面積小于出口截面面積，進口起到了孔板截流的效果，孔板截流與渦流截流裝置截流量的疊加，使得在1/5、2/5、3/5、4/5時截流效果特別明顯；因此，此時集水井中的水位也較高。4）進口高度不同，出現(xiàn)S轉(zhuǎn)折時的集水井水位并不相同，這說明渦流截流裝置的截流特性不僅與集水井中的水位有關(guān)，還可能與渦流截流裝置內(nèi)部的水位高低有關(guān)。

注：Q11-15為渦流裝置的S拐點結(jié)束后的出口質(zhì)量流量；H11-15為該出口流量對應(yīng)下的集水井中的水位。Q21-25為渦流裝置的S拐點前的出口質(zhì)量流量；H21-25為該出口流量對應(yīng)下的集水井中的水位。

3.2 渦流截流裝置的截流效率

在流體中沒有形成螺旋氣柱時，渦流截流裝置沒有出現(xiàn)截流特性，當(dāng)流體出現(xiàn)了S型的出口流量曲線時，才會有截流特性。因此，本文的截流效率以出現(xiàn)S型的出口流量曲線為基準(zhǔn)，是在出現(xiàn)S拐點前的流量與出現(xiàn)S拐點前的流量減去S拐點結(jié)束后的流量之比，截流效率公式為

在方程（7）中，為截流效率，%，m為出現(xiàn)S拐點前的出口質(zhì)量流量，kg/s。n為S拐點結(jié)束后的出口質(zhì)量流量，kg/s。

圖4為不同進口高度時，渦流截流裝置的截流效率曲線。

圖4 渦流截流裝置的截流效率

從圖4中可以看出，1）進口高度為1/5時，截流效率最大，可以產(chǎn)生30%以上的截流量，但是該進口高度下，渦流截流裝置進口截面積過小，一旦有較大雜質(zhì)進入雨水管道，發(fā)生堵塞的可能性最大。同時，由于截流量過大，使得集水井中的水位上升很快，在較小出口流量下就會使集水井發(fā)生溢流，如果在集水井上沒有安裝引流管或者溢流堰之類的裝置，將會使上游地區(qū)很快發(fā)生溢流或者內(nèi)澇等現(xiàn)象，無法達到改善城市雨水截流、降低內(nèi)澇發(fā)生的目的。2）隨著進口高度的增加，截流效率先減小后增大，轉(zhuǎn)折點在3/5～4/5之間，但是，從截流不發(fā)生堵塞的前提出發(fā)，4/5～5/5之間應(yīng)為最優(yōu)設(shè)計，該區(qū)間也是渦流截流裝置進口截面面積大于出口截面面積。3）結(jié)合圖3，我們可以看出，為了保證在集水井中流量較小時不發(fā)生截流，降低集水井本身的溢流問題以及大的雜質(zhì)造成堵塞的問題，又能提高截流效率，進口高度5/5，為最優(yōu)的進口高度選擇。

3.3 渦流截流裝置的內(nèi)部流動特性

通過上述分析，我們知道截流特性曲線均為S型，因此，其內(nèi)部流動形態(tài)差別不是很大，且進口高度5/5，為最優(yōu)進口高度，本文僅對進口高度5/5時，數(shù)值模擬形成S型截流曲線后，渦流截流裝置內(nèi)部=0，=300，數(shù)值計算時間=55 s，此時的中心截面進行分析，用來解釋渦流截流裝置截流特性的內(nèi)部機理。

3.3.1 渦流截流裝置=0截面上的流動特性

圖5a、b、c分別為=0時，渦流截流裝置內(nèi)部流場的壓力等值線云圖、速度等值線云圖、速度流線云圖。

1）從圖5a中可以看出，在集水井中的靜壓力分布均勻，從下到上逐步減小。但是經(jīng)過渦流截流裝置的作用，壓力分布不再均勻，在渦流截流裝置內(nèi)部形成了低壓區(qū)，甚至出現(xiàn)負(fù)壓，且低壓區(qū)向管道下游延伸。

2）從圖5b中可以看出，在湍流充分發(fā)展后，經(jīng)過渦流截流裝置的作用，流體的速度增加，這說明渦流截流裝置使得流體的靜壓力轉(zhuǎn)換成了動壓力。增加流體速度，有利于管道中的懸浮顆粒物質(zhì)加速排除。

3）從圖5c中可以看出，在整個流場中，存在大小不同的渦流，主要存在于集水井的底部，氣液兩相交界面，渦流截流裝置內(nèi)部上側(cè)，渦流截流裝置外部上側(cè)的集水井中，這說明，渦流截流裝置的存在，使得流體獲得動壓力，從而造成了障礙物阻流。

圖5 z=0時，渦流截流裝置內(nèi)部流場分布

3.3.2 渦流截流裝置=300截面上的流動特性

圖6a、b、c分別為=300時，渦流截流裝置內(nèi)部流場的壓力等值線云圖、速度等值線云圖、速度流線云圖。1）從圖6a中可以看出，在渦流截流裝置內(nèi)部，靜壓力分布從渦流裝置外部向渦流裝置中心遞減。2）從圖6b中可以看出，在渦流截流裝置內(nèi)部，速度分布從渦流裝置外部向渦流裝置中心遞增。3）從圖6c中可以看出，在渦流截流裝置內(nèi)部，形成了高速旋流。上述現(xiàn)象說明，渦流截流裝置之所以具有S型截流特性，是因為在渦流截流裝置內(nèi)部產(chǎn)生了低壓甚至負(fù)壓，靜壓力降低使得動壓力增加，流體流速增加，在渦流截流裝置該截面上產(chǎn)生的圍繞低壓帶的高速旋流，使得流體在裝置內(nèi)部旋轉(zhuǎn)流出，從而產(chǎn)生了旋流截流。

圖6 x=300時，渦流截流裝置內(nèi)部流場分布

3.4 試驗驗證

為了驗證數(shù)值模擬的準(zhǔn)確性，本文進行了實驗室試驗。通過對本文的數(shù)值模擬結(jié)果進行分析，我們發(fā)現(xiàn)盡管進口截面積小，渦流裝置的截流量大，但是從客觀實際出發(fā)，這樣很容易發(fā)生堵塞。

因此，本文數(shù)值模擬的最優(yōu)結(jié)構(gòu)為進口截面面積大于出口截面面積的渦流截流裝置，即進口高度4/5~5/5之間，因此僅對進口高度5/5的渦流截流裝置進行試驗驗證。本試驗在水利部杭州機械設(shè)計研究所（水利部產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所）水利機械及其再制造技術(shù)浙江省工程實驗室進行。試驗流程圖如圖7所示。

圖7 渦流截流裝置試驗流程圖

圖8為渦流截流裝置，進口高度為5/5時，數(shù)值模擬和試驗對比圖。從圖中可以看出，數(shù)值模擬具有較高的可信度，可以用來分析渦流截流裝置的截流特性和內(nèi)部流場。

圖9為實驗室試驗圖片，其中a、b分別給出了試驗裝置的實物模型和試驗裝置的實物模型及渦流裝置試驗效果圖。從圖中可以看出，渦流截流裝置內(nèi)部形成了螺旋氣柱，且氣柱的延伸段很長，說明本文數(shù)值模擬準(zhǔn)確。

圖8 數(shù)值模擬驗證

圖9 實驗室試驗圖片

4 結(jié) 論

本文采用FLUENT流體動力學(xué)數(shù)值計算軟件，模擬了不同進口高度的渦流截流裝置內(nèi)的復(fù)雜流動和截流特性之間的關(guān)系，同時，進行了截流特性試驗驗證，給出了不同進口高度對渦流截流裝置的截流特性、截流效率的影響規(guī)律，并重點分析了最優(yōu)截流就夠設(shè)計的內(nèi)部流動特性，揭示渦流截流裝置的截流機制。通過分析可以得到以下結(jié)論：

1）渦流截流裝置之所以具有S型截流特性，是因為在渦流截流裝置內(nèi)部產(chǎn)生了低壓甚至負(fù)壓，靜壓力降低使得動壓力增加，流體流速增加，在渦流截流裝置內(nèi)部截面上產(chǎn)生了圍繞低壓帶的高速旋流，使得流體在裝置內(nèi)部旋轉(zhuǎn)流出，從而產(chǎn)生了旋流截流。

2）在所有進口高度下，渦流截流裝置均產(chǎn)生了S型轉(zhuǎn)折，起到了截流作用。進口高度越小，發(fā)生截流時的出口流量越小，最終的截流效果越明顯。進口高度不同，出現(xiàn)S轉(zhuǎn)折時的集水井水位并不相同，渦流截流裝置的截流特性不僅與集水井中的水位有關(guān)，還與渦流截流裝置內(nèi)部的水位高低有關(guān)。隨著進口高度的增加，截流效率先減小后增大，轉(zhuǎn)折點在3/5～4/5之間，但是為了保證在低水位時候不進行截流，且進口截面積越小，越容易發(fā)生堵塞，因此，當(dāng)渦流裝置進口寬度給定時，渦流截流裝置進口截面面積大于出口截面面積為最優(yōu)結(jié)構(gòu)設(shè)計，即進口高度應(yīng)為4/5～5/5之間。

3）通過試驗驗證表明，本文研究的數(shù)值模擬具有較高的可信度，可以為城市污水截流裝置提供新的可選設(shè)計。同時，通過觀察試驗現(xiàn)象，確實是因為在渦流裝置內(nèi)部產(chǎn)生螺旋型旋流氣柱，改變了出水口過流斷面，從而限制出水流量，使管道出口質(zhì)量流量減小，產(chǎn)生了截流效應(yīng)。

[1] Wijesiri B, Egodawatta P, Mcgree J, et al. Influence of pollutant build-up on variability in wash-off from urban road surfaces[J]. Science of the Total Environment, 2015, 527: 344－350.

[2] Flanagan K, Branchu P, Ramier D, et al. Evaluation of the relative roles of a vegetative filter strip and a biofiltration swale in a treatment train for road runoff[J]. Water Science & Technology, 2017, 75(4): 987.

[3] 宋貞. 低影響開發(fā)模式下的城市分流制雨水系統(tǒng)設(shè)計研究[D]. 重慶：重慶大學(xué)，2014.

Song Zhen. Study on Design of Urban Diversion System Rainwater System Under Low Impact Development [D]. Chongqing: Chongqing University, 2014. (in Chinese with English abstract)

[4] Kyoung J A, Yun F L, Song H, et al. Multi-purpose rainwater harvesting for water resource recovery and the cooling effect[J]. Water Research, 2015, 86: 116－121.

[5] Rwehumbiza F, Hatibu N, Machibya M. Land characteristics, run-off and potential for rainwater harvesting in semi-arid areas of Tanzania[J]. Tanzania Journal of Agricultural Sciences, 1999, 2(2): 141－149.

[6] Sendanayake S. Rainwater Harvesting For Urban Living[M]. South Asian: South Asian Institute of Technology and Medicine, 2016.

[7] 靳云輝，吳學(xué)昌，石小飛. 城市截污工程中常用的截流限流設(shè)施[J]. 科技資訊，2008(5)：87－89.

Jin Yunhui, Wu Xuechang, Shi Xiaofei. The commonly used closure limiting facilities information of[J]. Science and Technology City Sewage Interception Project, 2008(5): 87－89. (in Chinese with English abstract)

[8] Kaiser M. Dezentrale regenwasserbewirtschaftung als baustein einer nachhaltigen siedlungsentwicklung[J]. Raumforschung und Raumordnung, 2006, 64(2): 126－134.

[9] 李衛(wèi)群，程瑞，方韜. 雨水管截污裝置設(shè)計探討[J]. 中國市政工程，2014(2)：29－31.

Li Weiqun, Cheng Rui, Fang Tao. Discussion on design of sewage interception device for storm sewer[J]. Chinese municipal engineering rainwater sewage pipe, 2014(2): 29－31. (in Chinese with English abstract)

[10] Santos M S, Alves A, Madeira L M. Chemical and photochemical degradation of polybrominated diphenyl ethers in liquid systems-A review[J]. Water Research, 2016, 88: 39－59.

[11] 焦愛萍，劉沛清，劉憲亮. 基于計算流體動力學(xué)的不同流量比兩股射流水墊塘的消能特性研究[J]. 水科學(xué)進展，2008，19(2)：238－244.

Jiao Aiping, Liu Peiqing, Liu Xianliang. Study on energy dissipation of two-jet with varying discharge ratio on a plunge pool based on computational fluid dynamics[J]. Advances in Water Science, 2008, 19(2): 238－244. (in Chinese with English abstract)

[12] 崔寶玲，呂子強，陳德勝，等. 起旋器入射角度對旋進旋渦流量計性能的影響[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(2)：53－58.

Cui Baoling, Lü Ziqaing, Chen Desheng, et al. Influence of incident angle of swirler on performance of swirl meter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(2): 53－58. (in Chinese with English abstract)

[13] 趙偉國，盛建萍，楊軍虎，等. 基于CFD的離心泵優(yōu)化設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(21)：125－131.

Zhao Weiguo, Sheng Jianping, Yang Junhu, et al. Optimization design and experiment of centrifugal pump based on CFD[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(21): 125－131. (in Chinese with English abstract)

[14] 程正飛，王曉玲，呂鵬，等. 基于VOF法的碾壓混凝土壩自由滲流場數(shù)值模擬[J]. 水利學(xué)報，2015，46(5)：558－566.

Cheng Zhengfei, Wang Xiaoling, Lv Peng, et al. Numerical simulation of free seepage field of RCC dam based on VOF method[J]. Journal of hydraulic engineering, 2015, 46(5): 558－566. (in Chinese with English abstract)

[15] 季磊磊，李偉，施衛(wèi)東，等. 導(dǎo)葉式混流泵內(nèi)部非定常流動特性數(shù)值模擬[J]. 農(nóng)業(yè)機械學(xué)報，2016(S1)：155－162.

Ji Leilei, Li Wei, Shi Weidong, et al. Numerical simulation of unsteady flow characteristics in mixed-flow pump with guide vanes[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016(S1): 155－162. (in Chinese with English abstract)

[16] 吳波，陳志，李建明，等. 基于 CFD 正交試驗的螺旋槽干氣密封性能仿真研究[J]. 流體機械，2014, 42(1)：11—16.

Wu Bo, Chen Zhi, Li Jianming, et al. Numerical simulation research on the sealing performance of spiral groove dry gas seal by orthogonal experiments based on CFD method[J]. Fluid Machinery, 2014, 42(1): 11－16. (in Chinese with English abstract)

[17] 張勤昭，劉福生，王宏. 三偏心蝶閥的流場和阻力特性研究[J]. 流體機械，2013, 41(11)：1－5.

Zhang Qinzhao, Liu Fushneg, Wang Hong. Study of flow field and resistance characteristics of triple eccentric butterfly valve[J]. Fluid Machinery, 2013, 41(11): 1－5. (in Chinese with English abstract)

[18] 張義龍，翟計紅，馬敏杰. 基于CFD的廢水厭氧處理系統(tǒng)布水形式數(shù)值模擬[J]. 鐵道工程學(xué)報，2014，31(10)：98－102.

Zhang Yilong, Zhai Jihong, Ma Minjie. Numerical simulation of anaerobic wastewater treatment system distribution based on CFD[J]. Journal of Railway Engineering Society, 2014, 31(10): 98－102. (in Chinese with English abstract)

[19] 朱榮生，盧永剛，王秀禮，等. 基于非線性空化的離心泵內(nèi)部流動特性分析[J]. 原子能科學(xué)技術(shù)，2016，50(7)：1216－1223. Zhu Rongsheng, Lu Yonggang, Wang Xiuli, et al. Analysis of internal flow characteristics of centrifugal pump based on nonlinear cavitation[J]. Atomic Energy Science and Technology, 2016, 50(7): 1216－1223. （in Chinese with English abstract）

[20] Lopes R J G, Quintaferreira R M. Assessment of CFD Euler-Euler method for trickle-bed reactor modelling in the catalytic wet oxidation of phenolic wastewaters[J]. Chemical Engineering Journal, 2010, 160(1): 293－301.

[21] 田濟揚，白丹，于福亮，等. 基于 Fluent 軟件的滴灌雙向流流道灌水器水力性能數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(20)：65－71.

Tian Jiyang, Bai Dan, Yu Fuliang, et al. Numerical simulation of hydraulic performance on bidirectional flow channel of drip irrigation emitter using Fluent[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(20): 65—71. (in Chinese with English abstract)

[22] 黃炎，趙滿全. 基于數(shù)值模擬與風(fēng)洞試驗的旋風(fēng)分離式集沙儀優(yōu)化設(shè)計[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(16)：50－56.

Huang Yan, Zhao Manquan. Optimization design of performance test of cyclone separator sand sampler based on numerical simulation and wind erosion tunnel experiment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(16): 50－56. (in Chinese with English abstract)

[23] 劉新陽，羅金耀，高傳昌，等. 滴灌用水力旋流器中顆粒分離的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(2)：7－11.

Liu Xinyang, Luo Jinyao, Gao Chuanchang. Numerical simulation of particle separation in hydrocyclone for drip irrigation system.[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(2): 7－11. (in Chinese with English abstract)

[24] 王福軍. 計算流體動力學(xué)分析[M]. 北京：清華大學(xué)出版社, 2004.

[25] 于勇. FLUENT入門與進階教程[M]. 北京：北京理工大學(xué)出版社，2008.

[26] 朱紅鈞. FLUENT 15.0流場分析實戰(zhàn)指南[M]. 北京：人民郵電出版社，2015.

[27] 王東瑋，王三反，李樂卓，等. 水處理反應(yīng)器的計算流體力學(xué)[M]. 北京：中國建筑工業(yè)出版社，2016.

[28] 宋鵬云，焦鳳，朱孝欽. 過程流體力學(xué)[M]. 北京：化學(xué)工業(yè)出版社，2016.

[29] 中華人民共和國住房和城鄉(xiāng)建設(shè)部. 室外排水設(shè)計規(guī)范(2014年版)：GB50014-2006 [S]. 北京：中國計劃出版社，2014.

馬光飛，方 勇，吳燕明，程文韜，李 桃，李 超.不同進口高度排水管道用渦流裝置截流特性模擬與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(22)：46－52. doi：10.11975/j.issn.1002-6819.2017.22.006 http://www.tcsae.org

Ma Guangfei, Fang Yong, Wu Yanming, Cheng Wentao, Li Tao, Li Chao. Closure characteristics simulation and experiment of eddy current device with different imported height in drainage pipeline[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(22): 46－52. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.22.006 http://www.tcsae.org

Closure characteristics simulation and experiment of eddy current device with different imported height in drainage pipeline

Ma Guangfei, Fang Yong※, Wu Yanming,Cheng Wentao,Li Tao, Li Chao

(1.310024; 2.310024)

In order to provide a new alternative design scheme for the interception treatment of urban sewerage, and reduce the effluent overflow in the catchment or the incidence of water logging in the pipe network, the closure characteristics of the eddy current closure device, the static pressure field distribution, the velocity field distribution and the velocity streamline distribution of the internal flow field at different inlet height are investigated. In this article, the transient N-S equations of incompressible flow, the renormalization group RNG k-ε turbulence model and the VOF (volume of fluid) method of Euler model are selected to carry out three-dimensional numerical simulation on the eddy closure device. Under the change of inlet height, the influences of the eddy current on the closure properties and closure efficiency of closure device are analyzed, as well as the internal flow characteristics of device. At the same time, the closure characteristics of the eddy current interception device have been experimentally verified. The experimental results are in good agreement with that of the numerical simulation. Based on the analysis of the riveting characteristics and the internal flow field characteristics of the eddy current device, the results show that under all the inlet height, the eddy current closure device produces the S-type turning and plays the role of closure. The smaller the inlet height, the stronger the intercepting capacity of the eddy current closure device; but the possibility of overflow and blockage is also bigger. At the different the inlet height, when the S-type turning occurs, the water level in the catchment is not the same, indicating that the cutoff characteristics of the eddy current closure device not only are related to the water level in the water collecting well, but also may be related to the water level inside the eddy current closure device. When the inlet height is small, the area of the inlet cross section of the eddy current closure device is smaller than that of the outlet cross section, and the inlet function is the effect of the orifice interception. The cutoff is the superposition of the interception of the orifice and the cutoff of the device. With the increase of the inlet height, the cutoff efficiency of the eddy current closure device decreases at first and then increases, and the increase range is between 4/5 and 5/5 of the inlet height, indicating that to ensure no overflow or clogging and achieve good closure, the inlet area of eddy current closure device should be greater than the export cross-sectional area. In the whole flow field, there are different vortices of different sizes, which are mainly present in the bottom of the trap, the gas-liquid interface, and the upper part of the internal and external eddy current closure device, and the existence of the eddy current closure device makes the fluid obtain dynamic pressure, resulting in obstacles obstruction. In the eddy current closure device, the high-speed swirl forms. The eddy closure device has S-type shutoff characteristics because of the low pressure or even negative pressure inside the eddy current closure device. With the decrease of the static pressure, the dynamic pressure increases, and the fluid flow rate increases; on the transverse cross section of the eddy current closure device, a high-speed swirl is produced around the low pressure zone, so the fluid flows out rotationally from the swirl device, resulting in the closure of cyclone flow. Through the experimental verification, it can be explained that the high speed vortex generated by fluid is a spiral gas column, and has a performance of long extension, which is beneficial to realize cyclone cleaning.

eddy currents; water conservancy equipment; inlet height; transient numerical simulation; closure characteristics; internal flow characteristics; experimental verification

10.11975/j.issn.1002-6819.2017.22.006

TD744; TU992

A

1002-6819(2017)-22-0046-07

2017-06-15

2017-10-16

浙江省公益技術(shù)應(yīng)用研究項目（2017C33010）；水利部公益性行業(yè)科研專項經(jīng)費項目（201501026）

馬光飛，工程師。主要從事水利機械優(yōu)化設(shè)計及流體機械方面的研究。Email：maguangfei@126.com

方 勇，高級工程師。主要從事水利機械方面的研究。Email：fangyongshuilibu@126.com