基于CFD的匯流三通管阻力特性與流場可視化分析

尚星航 季日臣

摘要：三通管是輸配水系統(tǒng)中常用的形式。在同一邊界條件下選用規(guī)格為30°、45°、60°、90°夾角三通管，采用不同匯流速度比，結(jié)合CFD數(shù)值模擬軟件進行分析計算，得到三通管道輸水時局部水頭損失系數(shù)變化規(guī)律。結(jié)果表明：局部阻力系數(shù)在雷諾數(shù)大于3×105之后基本保持穩(wěn)定；在同一流速比下直管-主管流向的局部水頭損失主要是由流速梯度變化和直管渦流運動造成的；當角度為90°時，引發(fā)的水流回流撞擊比較強烈，實際應用中應盡量減少三通管夾角的大小以減少能量損失，提高過流能力。

關(guān)鍵詞：匯流三通管；局部阻力系數(shù)；分流比；數(shù)值模擬

中圖分類號：TV134文獻標志碼：A

基金項目：國家自然基金資助項目（項目編號：51669010）。

0引言

近年來在急速凸顯的節(jié)水灌溉技術(shù)中，滴灌、噴灌等新型技術(shù)節(jié)水效果明顯，尤其在甘肅、新疆、寧夏等地，其節(jié)水效率比普通灌提高近50%。目前，隨著三通管在多項灌溉技術(shù)中的應用增多，三通管已經(jīng)逐漸成為輸水系統(tǒng)中一種常見的輸水設(shè)施。多年來，國內(nèi)外學者對三通管進行了廣泛研究：陳江林等[1]通過試驗研究和數(shù)值模擬對T型三通管水流的流動特性進行分析，得到了不同分流比、管徑比對水頭損失系數(shù)的影響；許虎等[2]通過對不同出口夾角下的斜直三通管進行數(shù)值模擬得到不同入口雷諾數(shù)、斜直三通管主管與支管間夾角、局部水頭損失系數(shù)三者之間相關(guān)關(guān)系；石喜[3]等通過對不同PVC三通管進行了試驗與數(shù)值模擬研究，揭示了內(nèi)部流場分布特性及局部阻力系數(shù)與雷諾數(shù)、管徑比的變化關(guān)系；李玲等[4]對三岔管內(nèi)水流的流動進行數(shù)值模擬與實驗研究，分析初步得出岔管內(nèi)能量損失的一些特性。

國內(nèi)學者對分岔管道還進行了大量的研究[5～10]，但國內(nèi)研究主要集中在分流三通管的阻力、流場特性和在不同流態(tài)下的沿程水頭損失或局部水頭損失以及阻力系數(shù)的變化，而對匯流三通管研究頗少，而且顯而易見僅利用試驗與理論來分析此類問題，會變得繁瑣而困難。本論述采用計算流體動力學技術(shù)（CFD），運用標準k-ε紊流模型，進行不同角度、流速下三通管的阻力特性和流場可視化分析，研究引起阻力、流場特征變化及局部水頭損失的原因。

1模型與方法

1.1湍流數(shù)學模型建立

Fluent提供的湍流模型包括：單方程模型、雙方程模型及雷諾應力模型和大渦模擬。標準模型需要求解湍動能及其耗散率方程。湍動能輸運方程是通過精確的方程推導得到，但耗散率方程是通過物理推理，數(shù)學上模擬相似原形方程得到的。該模型假設(shè)流動為完全湍流，分子粘性的影響可以忽略。因此，標準模型只適合完全湍流的流動過程模擬。

1.2Gambit建立幾何模型及網(wǎng)格劃分

因模型內(nèi)部結(jié)構(gòu)復雜，為考慮各部分網(wǎng)格之間的耦合問題，Gambit的自動功能無法完成，所以本論述對模型進行了處理，將模型切割分成幾個部分進行網(wǎng)格劃分，然后將各部分的網(wǎng)格連成一個完整的模型，對流道突變的地方進行更密的網(wǎng)格劃分。

研究建立模型以及網(wǎng)格劃分、邊壁處理都在Fluent的前處理軟件Gambit軟件進行。三種規(guī)格的三通管均采用混合網(wǎng)格，在三通管分岔處采用網(wǎng)格局部加密，進行處理效果較好的Hex/Wedge與Cooper網(wǎng)格，其余部位劃分采用Tet/Hybrid與TGrid結(jié)合。30°、45°、60°、90°三通管網(wǎng)格數(shù)分別為4.18×105、4.26×105、4.45×105、4.61×105，網(wǎng)格劃分如圖1所示。

2三通管計算與阻力特性分析

2.1三維模型局部阻力系數(shù)計算方法

根據(jù)三通管豎直管為供水進水口主管的工況，分別進行30°、45°、60°、90°T型三通管試驗，探討不同匯流流速比下三通管局部阻力系數(shù)ξ1、ξ2隨雷諾數(shù)的變化趨勢，計算公式見表1所列；

式中：

Δp1——主管三通管截面；

Δp2——側(cè)管三通管截面；

vo——主管三通管匯流后水流流出速度；

S0——主管三通管截面；

S2——側(cè)管三通管截面；

Qo——主管流出三通管流量；

Q2——側(cè)管流入三通管流量。

2.2阻力系數(shù)計算成果及阻力特性分析

圖2、圖3分別給出了不同匯流速度比下90°、60°、45°、30°三通管局部阻力系數(shù)ξ1、ξ2隨雷諾數(shù)的變化趨勢。從圖中可以看出，四種不同角度的三通管中：90°三通管的阻力系數(shù)隨雷諾數(shù)的增大而呈減小趨勢；60°、45°、30°三通管的阻力系數(shù)隨雷諾數(shù)的增大而增大。本實驗中雷諾數(shù)最小值為2.3×105，水流此時的流動狀態(tài)處于紊流。粘性流體處于紊流狀態(tài)時，當雷諾數(shù)大于一定程度后，基本不影響流動狀態(tài)和流速分布，此時流動處于阻力平方區(qū)，阻力系數(shù)基本保持穩(wěn)定，由結(jié)果看出當雷諾數(shù)Re≥3×105時，四種三通管的局部阻力系數(shù)變化幅度不大，幾乎處于平穩(wěn)狀態(tài)?？梢哉J為已經(jīng)進入阻力平方區(qū)。

圖4、圖5給出了同一匯流速度比下，四種不同角度的壓力云圖和流線圖，可以看出：水流由直管進入主管時，經(jīng)過轉(zhuǎn)彎處形成渦流和環(huán)流較為明顯，管內(nèi)流態(tài)變得極不穩(wěn)定，由于水流方向垂直，相互碰撞，能看到明顯的旋渦運動。當三通管分岔口角度為90°時，尺度極大，引發(fā)的水流回流撞擊也比較強烈；當三通管分岔口角度為60°時，流速梯度減弱，水流之間的碰撞速率降低，旋渦運動變?。划斎ü芊植砜诮嵌葹?5°和30°時水流流速平順，流速分布比較均勻，未出現(xiàn)較大的流速梯度，也未形成渦流。通過分析可得：直管-主管流向的局部水頭損失主要是由流速梯度變化和直管渦流運動造成；隨著角度的減小，流速梯度減弱，管內(nèi)水流流速趨于平順，渦流運動變小。

4結(jié)論

（1）在不同匯流速度比下90°、60°、45°、30°三通管局部阻力系數(shù)ξ1、ξ2隨雷諾數(shù)的增大而趨于穩(wěn)定，在雷諾數(shù)大于3×105時認為局部阻力系數(shù)與雷諾數(shù)無關(guān)，具體如表2、表3、表4、表5、表6所列。

（2）三通管夾角為90°時，引發(fā)的水流回流撞擊比較強烈，造成較大的能量損失，隨著角度的減小，流速梯度減弱，管內(nèi)水流流速趨于平順，渦流運動變小，在實際應用中應減少三通管夾角的大小以減少能量損失，提高過流能力。

（3）直管-主管流向的局部水頭損失主要是由流速梯度變化和直管渦流運動造成。

參考文獻：

[1]陳江林，呂宏興，石喜，等.T型三通管水力特性的數(shù)值模擬與試驗研究[J].農(nóng)業(yè)工程學報，2012，28（5）：73-77.

[2]許虎，吳文勇，王振華，等.基于CFD的斜三通管水力特性分析及流場計算[J].排灌機械工程學報，2020，38（11）：1138-1144.

[3]李玲，李玉梁，黃繼湯，等.三岔管內(nèi)水流流動的數(shù)值模擬與實驗研究[J].水利學報，2001（3）：49-53.

[4]石喜，呂宏興，朱德蘭，等.PVC三通管水流阻力與流動特征分析[J].農(nóng)業(yè)機械學報，2013，44（1）：73-79，89.

[5]華紹曾，楊學寧.實用流體阻力手冊[M].北京：國防工業(yè)出版社，1985.

[6]李玲，陸豪，陳嘉范.抽水蓄能電站尾水岔管水流運動及阻力特性試驗研究[J].水力發(fā)電學報，2008，27（3）：101-104，109.

[7]李玉梁，李玲，陳嘉范，等.抽水蓄能電站對稱岔管的流動阻力特性[J].清華大學學報（自然科學版），2003，43（2）：270-272，280.

[8]楊校禮，高季章，劉之平.三岔管水流數(shù)值模擬研究[J].水利水電技術(shù)，2005，36（1）：48-50.

[9]楊校禮，高季章，劉之平.有加強肋板的三岔管水流數(shù)值模擬及水頭損失研究[J].水力發(fā)電，2004，30（5）：18-20，27.

[10]李玲，李玉梁，黃繼湯，等.三岔管內(nèi)水流流動的數(shù)值模擬與實驗研究[J].水利學報，2001（3）：49-53.