基于Fluent的T型三通內(nèi)部流場數(shù)值模擬

2017-01-11 by:CAE仿真在線來源:互聯(lián)網(wǎng)

鐘振亞

山東豪邁化工技術

引言

三通管在工業(yè)、農(nóng)業(yè)諸如油田內(nèi)部集輸管網(wǎng)等很多方面使用非常廣泛, 主要用于改變流體方向,可用于輸水管路,輸油管路及各種液體化工材料輸送管路。因此對三通結構流動形式的研究可以為管道的設計提供一定的依據(jù)。采用Fluent對三通結構進行模擬可以降低費用,節(jié)約成本,具有比價好的經(jīng)濟效益,同時也可以減小試驗儀器的不穩(wěn)定度于試驗結果的影響。

1 數(shù)學模型的建立

對于不可壓縮粘性流體,其運動可用N-S方程標示,控制方程為:

能量守恒方程如下:

湍流動能運輸方程為:

湍流動能耗散輸運方程為:

2 問題描述及網(wǎng)格劃分

2.1

問題描述

密度為860Kg/m3 ,粘度為1.02Kg/m?s,比熱容為2100J/Kg?K 的原油以不同的速度和溫度分別從inlet1和inlet2進入三通管,混合后從outlet流出。Inlet1流入的的原油的溫度為300k,速度為3m/s。inlet2流入的原油的溫度為450k,速度為 4m/s,兩股原油在三通管內(nèi)發(fā)生劇烈的能量交換和能量損失, 混合后由outlet流出

2.2

模型的建立及網(wǎng)格劃分

計算模型如圖1所示。圖2為模型的網(wǎng)格劃分圖。在Solidworks中利用其三維建模功能,建立如圖 1所示尺寸的模型,并利用Ansys Fluent 的mesh 功能進行網(wǎng)格劃分,網(wǎng)格劃分方法選擇Hex Dominant , 網(wǎng)格大小設置為3,太小的網(wǎng)格在劃分過程中會因為計算機內(nèi)存不足而退出。網(wǎng)格劃分的過于稀疏又無法達到指定的計算精度。經(jīng)過反復試驗本文取網(wǎng)格大小為3。具體劃分結果如圖3所示。

圖3 (網(wǎng)格劃分)

2.3

邊界條件的設置

邊界條件的設置如表1所示:

3 模型的求解及結果分析

在fluent中對求解器進行設置,采用三維穩(wěn)態(tài)隱式求解方法。粘度的處理采用Standard k-epsilon turbulence model。設置收斂進度為0.0001,進行600次迭代后,計算結果如圖

計算結果如圖4~圖8所示。從圖4中可以看出,兩股流體在三通內(nèi)混合后壓力較混合前有較大幅度的降低,并且在豎管出口的右側形成一個低壓區(qū)。從圖4,圖5可以看出管由于粘性的影響,在管道壁面上的流體質點速度為零。由于管內(nèi)的流動是紊流,所以管中的速度分布不是通常的“拋物線形狀”,其速度分布(除管壁)比較均勻平坦。但在三通管的匯合處及其后的管道中與混合之前大致相同,由于處于湍流狀態(tài)混合后流速分布仍然未呈現(xiàn)“拋物線形狀”的速度分布,在混合后出現(xiàn)了一個低速區(qū),基本上與低壓區(qū)重合。這個區(qū)域內(nèi)的流速與管內(nèi)主流體流動方向不同。由此可以看出在混合后側靠近避免的位置出現(xiàn)了漩渦流動。從圖8可以看出,在小管流體噴出后在三通交匯處的下部形成兩個低壓區(qū),這是由于小管流體的射流作用造成的。此外,不同的邊界層網(wǎng)格劃分對計算結果也會有影響,網(wǎng)格的寬度及層數(shù)都應適合,太疏則不能反映流動狀態(tài)的變化,太密則會增大計算量。

4 結語

在計算過程中采用Segregated隱式解,標準K-Ε雙方程模型來處理湍流計算,經(jīng)結果驗證,此模型相對于其他湍流模型能獲得更合理,更精確的計算結果。本文通過對三通管道三維湍流流動計算,正確計算出三通管內(nèi)流體混合后流場和溫度場的分布情況。對于研究油品在通過三通時的流動狀況和管道的設計積累了經(jīng)驗,提供了依據(jù)。

參考文獻

[1]邱秀云.水力學[M], 烏魯木齊:新疆點子出版社, 2008.

[2]王福軍. 計算流體動力學分析. CFD 軟件原理與應用[M], 北京:清華大學出版社 2004.

[3]韓占中,王靜. Fluent- 流體工程仿真計算實力與應用[M],北京:北京理工大學出版社2010.

開放分享：優(yōu)質有限元技術文章,助你自學成才