斜三通的CFD數值模擬及傳熱影響因素研究

朱海榮 盧紅亮 崔海亭 鄧小葉

摘要：為研究斜三通管內流體的流動與傳熱過程，應用CFD軟件Fluent分別對管間夾角θ為20°～ 90°、流速比λ為2～3.5以及支管位置l為500～1 000 mm的斜三通內流體的流動與傳熱過程進行了數值模擬，并通過正交試驗設計對影響三通壁面平均換熱系數的各參數進行了敏感性分析。結果表明，隨著θ和λ的增大，三通出口處流體最大速度逐漸增大，其中λ的影響更為顯著;壁面平均換熱系數隨θ和λ的增大而增大，隨l的減小而增大;支管位置l主要影響壁面平均換熱系數，對最大速度的影響可忽略不計。按照影響壁面平均換熱系數作用的強弱，支管位置l >流速比λ>管間夾角θ，即支管越靠近主管進口、流速比越大、管間夾角越大，則壁面換熱效果越好。因此，支管位置是影響斜三通壁面傳熱效果的首要因素。

關鍵詞：節能技術;斜三通;流動與傳熱;管間夾角;流速比;支管位置

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

Abstract：In order to study the flow and heat transfer process of the fluid in an obliquetee， the CFD simulation of flow and heat transfer processes of the fluid in oblique tee with the included angle of 20°～90°， the velocity ratio of 2～3.5 and the branch position of 500～1 000 mm is carried out based on Fluent software， and the sensitivity analysis of various parameters influencing the average heat transfer coefficient of wall surface is conducted by orthogonal experimental design. The results show that as the included angle or velocity ratio increases， the maximum velocity of fluid at the outlet of oblique tee increases gradually， and the influence of velocity ratio is more significant. The average heat transfer coefficient of wall surface also increases with the increase of included angle and velocity ratio; the branch position mainly affects the average heat transfer coefficient of wall surface， and its influence on the velocity is negligible. The average heat transfer coefficient of wall surface is increasing with the decrease of the branch position. According to the different effects of various parameters on average heat transfer coefficient of wall surface， branch position l>velocity ratio λ>included angle θ， which means that the closer the branch pipe is to the inlet of main pipe， the greater the velocity ratio， and the greater the included angle， the better the heat transfer effect of wall surface will be. Therefore， the branch position is the primary factor to influence the heat transfer effect of oblique tee.

Keywords：energy conservation technology; ?oblique tee; flow and heat transfer; ?included angle; velocity ratio; branch position

三通可用于水、油及各種液體化工材料的輸送管路，廣泛應用于醫學、水利（節水灌溉、給水排水）、能源動力（石油、天然氣、核工業）、建筑等工程領域。三通作為主要的分流部件，當流體流經主管與支管的交匯處時，會在主管道截面上產生二次流動并在支管內產生渦旋結構，不但會導致流體總壓和能量的損失，還會降低傳熱、傳質的效率[1-2]。因此，研究三通在不同結構形式下的流動和傳熱特性以及影響其流動傳熱性能的主要因素是非常有必要的。

有學者已經對T型三通管內流體流動、傳熱特性、結構設計等方面進行了一些研究。曹海兵等[3]和侯聰等[4]運用數值模擬的方法分析了等直徑的 T 型、Y 型和圓弧型三通管內流體的流動狀況和阻力性能。鞏啟濤等[5-6]通過對三通管內流體速度場和壓力場的分析，研究了流速比λ、雷諾數Re、管徑比d及管間夾角θ對大管徑三通局部阻力的影響。楊康等[7]利用Fluent軟件分析了不同的流體參數和管路結構下三通管內流體的速度場、壓力場、溫度場和湍動能等。韓方軍等[8]通過對不同管徑直角三通管的分析得出局部水頭損失和靜壓力的變化。王海軍等[9]研究了不同流速比下T型三通管內的內壁面傳熱系數以及內部溫度波動。魏顯達等[10]通過對三通的速度、溫度、壓力進行數值模擬對管道進行結構優化和選材優化。陳江林等[11]和盧冬華等[12]通過數值模擬與實驗研究分析了三通水頭損失的機理及管內流動狀態。谷雨等[13]用標準k-ε模型對三通進行三維模擬，分析了管內流體壓力和溫度的變化。石喜等[14]對PVC 三通管進行了實驗與數值模擬研究，分析影響其局部阻力的主要因素。卜琳等[15]、龔旭等[16]數值模擬了斜三通管內流體的流動結構、壓力分布與溫度分布情況，分析了斜三通幾何結構、流動條件和流體溫度變化對管內流體流動結構、溫度分布及熱應力的影響。孫鑫[17]通過對T型三通和圓弧型三通內流體的湍流流動進行數值模擬和實驗研究，分析其內部的流場分布，初步研究了三通對大氣出流情況下的水頭損失率。康朝陽[18]分析了在三通管內流體發生瞬態流動時，流體的速度場、壓力場的變化。SELVAM等[19]和GAO等[20]通過數值模擬分析驗證了流速比與射流深入主流的程度成正比，在不同區域流體與壁面的換熱程度也不同，在射流管與主流管的交匯處溫度波動最大。ANDREOPOULOS等[21-22]對當射流垂直入射至主流體中且不同流速比時射流與主流的平均速度、紊動能及平均溫度分布情況進行了研究;SHERIF等[23]對當水洞中熱水垂直向上注入冷流中且速比為1，4和7時的平均溫度及波動溫度進行了研究。

現有文獻從不同方面研究了三通管內流體流動時溫度、壓力的分布情況以及阻力和湍流強度的變化規律，但是還缺乏對影響三通管內流體流動傳熱性能因素的顯著性分析，而顯著性分析的結果恰恰可以有效指導三通的結構優化和工程應用。本文應用Fluent軟件數值模擬了斜三通在不同結構參數（如支管位置l、管間夾角θ）和運行工況（如流速比）下的流體流動狀況、速度場與平均換熱系數的分布情況，并得到了各參數對三通壁面平均換熱系數的影響規律。

1模型的建立

1.1物理模型

平均換熱系數是指流體與管壁面的平均對流換熱系數，不考慮三通管的管體內外徑和壁厚以及其內部的熱傳導。本文主要研究兩股流速、溫度不同的流體，分別從主管、支管進入斜三通，在三通管內混合流動并進行傳熱，最終從主管出口流出的過程。斜三通物理模型如圖1所示，主管和支管的管徑均為150 mm，主管總管長為1 500 mm，支管管長為500 mm。

1.2網格模型

采用分塊劃分網格方法，以六面體網格為基本單元，在三通主管和支管交匯處進行網格加密。近壁面第1個網格單元中心至壁面的距離滿足壁面函數法的要求。如圖2所示，當管間夾角θ=60°，l=750 mm時，三通共劃分了 157 666個網格。進一步進行網格無關性驗證，計算結果無明顯改變，表明此時網格滿足計算要求且較為合理。

2控制方程及邊界條件

2.1控制方程

3結果與討論

3.2壁面平均換熱系數的影響因素分析

圖6為λ=2，l=750 mm時，壁面平均換熱系數隨θ的變化規律。如圖6所示，支管流體的溫度較低，θ較小時支管流體與主管流體的混合不明顯，兩股流體各自平穩地向出口流動，支管壁面的換熱系數較大，換熱系數最大值出現在支管與主管的交匯處;隨著θ的增大，支管流體進入主管后被迫改變方向并與主管流體混合碰撞，產生渦流，兩股流體在主管內混合換熱，但是主要集中在主管軸線附近，因此壁面平均換熱系數略有下降;隨著θ的進一步增大，兩股流體產生了強烈的回旋渦流，在管道交匯處的湍流混合加劇，從而導致換熱系數的大幅度提高;θ=90°時，兩股流體垂直相遇，碰撞最為劇烈，引起強烈的能量和動量的交換，此時壁面換熱系數達到最大值。同時，由于流動方向的改變而產生的漩渦也引起了水頭損失，在θ=90°時支管產生的水頭損失最大。

圖7所示為θ=45°，λ=2時，壁面平均換熱系數隨l的變化規律。無論支管位置距離主管進口遠近，兩股流體交會時的流動狀態幾乎一樣，不同支管位置的壁面平均換熱系數最大值基本相同。然而，支管離主管進口越近，支管流體越早與主管流體混合。因此，相同時間內，混合后流體與壁面的對流換熱量會隨著l的減小而增大。

圖8為斜三通壁面平均換熱系數變化曲線圖。在相同λ下，隨著θ的增大，三通壁面平均換熱系數先減小后增大，但整體呈上升趨勢。在相同θ下，λ越大，壁面換熱系數越大。當λ較小時，支管流體匯入主管比較緩慢，兩股流體無法形成強烈的混合，近似保持層流狀態向出口流動，因此壁面換熱系數較小;隨著λ的增大，支管流體高速匯入主管，兩股流體強烈混合，在交匯處產生了強烈的回旋渦流，兩流體的湍流混合加劇，加速了流體間以及流體與壁面間的接觸換熱，從而使壁面平均換熱系數提高。

3.3正交試驗分析

對不同參數（支管位置l、管間夾角θ、流速比λ）下三通管內流體的傳熱過程進行正交試驗分析，正交試驗方案如表4所示。

表4所列的9組實驗結果中，以第3組實驗得出的壁面平均換熱系數最大，其水平組合為A3B3C3，分別是各影響因素中數值最大的組合。由此可以看出，本實驗各因素組合中的最優組合為A3B3C3。

表5為各參數的敏感性分析結果，表中k1，k2，k3為K1，K2，K3的平均值;R為極差，即每個因素下k的最大值與最小值的差值。通過R值的大小可以看出本實驗中各參數的影響作用存在顯著性順序，其主次關系為A>C>B，因此影響三通管壁面傳熱效果的因素依次為支管位置l、支管與主管的流速比λ、管間夾角θ。支管越靠近主管進口、流速比越大、管間夾角越大，則三通壁面的換熱效果越好。

4結論

利用Fluent軟件數值模擬了三通管內流體的流動傳熱過程，獲得了不同管間夾角、支管位置以及流速比對三通管內流體流動與傳熱的影響，并運用正交試驗設計分析了各參數對壁面換熱效果的影響規律，得出以下結論。

1）隨著λ和θ的增大，三通出口處流體的最大速度都呈現增大的趨勢，且λ對流體最大速度的影響要高于θ;而支管位置l對兩股流體交匯時的流動狀態影響較小。可見出口流體最大速度主要取決于支管流體的流速。

2）三通壁面的平均換熱系數隨θ的增大先略有減小后增大，隨λ的增大持續增大，隨l的減小而逐漸增大，說明θ和λ的增大均能使兩股流體的混合更為充分，加劇了能量的交換和熱量的傳遞，使壁面平均換熱系數提高;而l的減小意味著支管流體提前匯入主管，相同時間內增大了流體與壁面間的換熱量，從而使整個壁面的平均換熱系數增大。

3）通過正交試驗分析影響三通壁面平均換熱系數的因素，按照影響壁面傳熱效果的程度依次為支管位置l>流速比λ>管間夾角θ。即支管越靠近主管進口、流速比越大、管間夾角越大，壁面的換熱效果越好。因此，支管位置是影響斜三通壁面傳熱效果的首要因素。

4）筆者對全面認識三通管內流體的流動和傳熱過程提供了參考，然而三通管內流體在交匯時伴隨著能量損失和阻力性能的變化，因此在最小能量損失和阻力下研究流體輸送和流體與管壁之間的對流換熱是未來的一個研究方向。

參考文獻/References：

[1]魏志，王玉璋. 閥體后90°圓形彎管內部流場的數值模擬[J]. 動力工程學報，2013，33（8）：595-599.WEI Zhi，WANG Yuzhang. Numerical simulation on inner flow in a 90° circular-sectioned bend with fore-endvalve[J].Journal of Chinese Society of Power Engineering，2013，33（8）：595-599.

[2]李雄，劉偉軍，唐飄.排氣三通管道分散流流動特性的數值模擬及分析[J]. 河北科技大學學報，2014，35（3）：272-278.LI Xiong，LIU Weijun，TANG Piao.Numerical simulation and analysis on the dispersed flow haracteristics in tee branch pipe of exhaust system[J]. Journal of Hebei University of Science and Technology，2014，35（3）：272-278.

[3]曹海兵，高莉麗. 不同結構三通管內流體流動特性數值模擬[J]. 廣東化工，2016，43（13）：72-73.CAO Haibing，GAO Lili. Numerical simulation on flow characteristics of different three-way pipe[J].Guangdong Chemical Industry，2016，43（13）：72-73.

[4]侯聰，王合旭，張方駒.不同結構三通內冷熱水混流特性的數值仿真[J]. 計算機輔助工程，2017，26（1）：7-12.HOU Cong，WANG Hexu，ZHANG Fangju. Numerical simulation on characteristics of hot and cold fluids mixed in tee junctions of different structures[J].Computer Aided Engineering，2017，26（1）：7-12.

[5]鞏啟濤，楊俊紅，韓奎，等. 大管徑T型三通數值模擬及局部阻力特性分析[J]. 動力工程學報，2016，36（9）：753-758.GONG Qitao，YANG Junhong，HAN Kui，et al. Characteristic analysis on the flow and local resistance in large pipe tees [J].Journal of Chinese Society of Power Engineering，2016，36（9）：753-758.

[6]鞏啟濤. 大管徑三通局部阻力特性數值模擬及多熱源集中供熱運行分析[D]. 天津：天津大學，2015.GONG Qitao. Numerical Simulation of Large Pipe Tees Local Resistance Characteristics and Analysis of Central Heating System with Multi-heat Sources Operation[D]. Tianjin：Tianjin University，2015.

[7]楊康，劉吉普，馬雯波. 基于FLUENT軟件的T型三通管湍流數值模擬[J]. 化工裝備技術，2008，29（4）：33-36.

[8]韓方軍，孫鑫，張原. T型三通管內部流場數值模擬與結構優化[J].新疆水利，2010，4：1-3.HAN Fangjun，SUN Xin，ZHANG Yuan.Numerical simulation and structural optimization of internal flow field in T-type tee[J].Xinjiang ?Water ?Resources，2010，4：1-3.

[9]王海軍，盧冬華，羅毓珊，等. 壓水堆主系統T型三通管傳熱實驗研究[J]. 核動力工程，2002，23（3）：17-21.WANG Haijun，LU Donghua，LUO Yushan，et al. Experimental study on heat transfer characteristic of T-type tee pipe of PWR main system[J].Nuclear Power Engineering，2002，23（3）：17-21.

[10]魏顯達，王為民，徐建普. 基于 ?Fluent ?的三通管數值模擬及分析[J]. 當代化工，2011，40（2）：165-167.WEI Xianda，WANG Weimin，XU Jianpu. Numerical simulation and analysis of fluid in three-way connection pipe based on fluent software[J]. Contemporary Chemical Industry，2011，40（2）：165-167.

[11]陳江林，呂宏興，石喜，等. T型三通管水力特性的數值模擬與試驗研究[J]. 農業工程學報，2012，28（5）：73-77.CHEN Jianglin，LV Hongxing， SHI Xi，et al. Numerical simulation and experimental study on hydrodynamic characteristics of T-types[J]. Transactions of the Chinese Society of Agricultural Engineering，2012，28（5）：73-77.

[12]盧冬華，村松壽晴. T型三通管內熱分層流動 3D 數值模擬[J].核動力工程，2005，26（4）： 332-334.LU Donghua，Toshiharu Muramatsu. 3-D numerical simulation of stripping and stratification in mixing T-junction[J]. Nuclear Power Engineering，2005，26（4）： 332-334.

[13]谷雨，李明達.基于CFD 的三維三通管道湍流數值模擬研究[J].科協論壇，2013，12（5）：234-235.

[14]石喜，呂宏興，朱德蘭，等.PVC 三通管水流阻力與流動特征分析[J].農業工程學報，2013，44（1）：73-89.SHI Xi，LV Hongxing，ZHU Delan，et al. Flow resistance and characteristics of PVC tee pipes[J].Transactions of the Chinese Society of Agricultural Engineering，2013，44（1）：73-89.

[15]卜琳，李會雄，王海軍，等.斜向三通管內射流混合特性的三維數值模擬[J].核動力工程，2008，29（5）：124-128.PU Lin，LI Huixiong，WANG Haijun，et al. Numerical simulation of 3D inclined jet to cross flow[J].Nuclear Power Engineering，2008，29（5）：124-128.

[16]龔旭，劉德明，吳昆.建筑排水三通入流角度對流態影響的數值模擬[J].福建建設科技，2015，5：83-85.GONG Xu，LIU Deming，WU Kun.Numerical simulation of the flow state affected by inflow angle in the tee coupling tube[J].Fujian Construction Science and Technology，2015，5：83-85.

[17]孫鑫. T型三通和圓弧型三通流場的試驗研究及數值模擬[D]. 烏魯木齊：新疆農業大學，2009.SUN Xin. Physical Experiment and Numerical Simulation of "T and arc" Shaped Three-way Pipe's Flow Field[D]. Urumchi： Xinjiang Agricultural University，2009.

[18]康朝陽.主蒸汽三通對流傳熱、熱應力及壽命的計算分析[D].上海：上海發電設備成套設計研究院，2015.KANG Chaoyang. Calculation and Analysis of the Heat Transfer Thermal Stress and Life of the Main Stream Tee[D]. Shanghai：Shanghai Power Equipment Research Institute，2015.

[19]SELVAM P K，KULENOVIC R，LAURIEN E. Experimental and numerical analyses on the effect of increasing inflow temperatures on the flow mixing behavior in a T-junction[J]. International Journal of Heat and Fluid Flow，2016，61：323-342.

[20]GAO K，WANG P，LU T，et al. Experimental investigation and numerical simulation for weakening the thermal fluctuations in a T-junction [J]. Annals of Nuclear Energy，2015，78：180-187.

[21]ANDREOPOULOS J，RODI W. Experimental investigation of jets in cross flow[J]. J Fluid Mech， 1984，138（1）：93-137.

[22]ANDREOPOULOS J. Heat transfer measurements in a heated jet-pipe flow issuing into a cold cross stream[J]. Phys Fluids，1983，26（11）：3201-3210.

[23]SHERIF S A，PLETCHER R H. Measurements of the thermal characteristics of heated turbulent jets in cross flow[J].Journal of Heat Transfer，1989，111（11）：897-903.