螺旋管內局部二次流強度計算及演變

前言：本站為你精心整理了螺旋管內局部二次流強度計算及演變范文，希望能為你的創作提供參考價值，我們的客服老師可以幫助你提供個性化的參考范文，歡迎咨詢。

摘要:為了掌握螺旋管內二次流強度的變化規律,對螺旋管內局部二次流強度及演變規律進行了研究,提出了無量綱數H作為局部二次流強度指標。通過數值模擬得到H的體積平均值Hm,與描述整體二次流強度的迪恩數進行對比。結果顯示,在不同的雷諾數、曲率、螺距下,Hm與迪恩數的相關系數約為1,且呈現較好的正比關系,說明H能夠準確計量局部二次流強度。同時,探究了幾何參數對螺旋管內的二次流演變規律的影響。依據二次流的強度,管內流動可分為發展段、過渡段和穩定段:二次流在發展段迅速增強;在過渡段出現最大值;在穩定段呈小幅度周期性波動。雷諾數增大會顯著增大二次流強度,但是不影響二次流在3個階段出現的位置。曲率的增大會略微增大二次流的強度,且延長了二次流達到穩定段的路程。螺距對二次流的強度和演變規律影響非常小。

關鍵詞:迪恩渦;迪恩數;局部二次流強度;螺旋管;有限體積法

二次流,是指一定的主流速度下和一定的幾何邊界條件下作曲線運動的黏性流體所產生的一種有規律的伴隨運動[1]。螺旋管中的流體由于受到離心力、科氏力、浮力等作用會產生的一對呈反向渦旋的二次流———迪恩渦[2]。迪恩渦有著獨特的流場結構:一方面迪恩渦可在避免湍流的情況下加強管內流體的擾動,使管內流體的橫向對流作用顯著增強[3-6];另一方面在迪恩渦作用下,靠近管壁的流體對管壁的剪切作用可以減薄邊界層厚度、防止管壁結垢[7]。因此,螺旋管被廣泛運用于石油、化工、制冷、空分等領域。威廉?雷金納德?迪恩在小曲率條件下得到了圓形截面彎管內流動的攝動解,并且提出了判斷管道內二次流的強度的準則數迪恩數(De)[2]。迪恩數代表的是離心力與黏性力的比值,同時也衡量著二次流的強度[8]。螺旋管作為強化傳熱的手段使用已久,以往對螺旋管的優化和二次開發通常是通過對比整體傳熱和流動特點來進行的。也有學者用了多種方法衡量彎管內的二次流強度,大體上可以分為4種。第1種是通過觀察流線圖來判定局部二次流強度,這種方法只能做定性的判斷,并且帶有主觀性。第2種是通過計算速度在垂直于主流方向的截面上的分量大小[9]。但是,在很多情況下垂直于主流方向的橫向運動并非由二次流引起,如在直管中的湍流流體也有垂直于主流方向的橫向運動,所以單純地計量垂直于主流方向的橫向運動速度不足以說明二次流的強度。第3種是通過計算各節點的流體速度與主流速度的比值來闡述螺旋管內二次流特性[5,10]。但螺旋管壁面附近流體速度受邊界層的影響,與主流速度的比值總是趨于0,因此該方法無法對邊界層流體的流動變化特性進行表達,也無法分析螺旋管結構對邊界層處流體的影響。第4種是通過計算局部渦量或縱向渦及所衍生出來的無量綱數來衡量二次流的強度[11]。但是彎管近壁面處速度變化劇烈,存在著不是因二次流而產生的渦量,且通常近壁面處的渦量非常大,這給觀察和計量主流區的二次流強度帶來了干擾。雖然衡量局部二次流強度的方法眾多,但是這些方法都存在著局限性,這就導致了工程人員難以精準地調控二次流強度。本文提出了一個描述局部二次流強度的無量綱指標H,驗證各種條件下的H的體積平均值與迪恩
數是否為線性關系,并探究了螺旋管內流體的二次流演變規律,為螺旋管內二次流的調控及強化管的開發提供了參考指標。

1理論分析

圖1為螺旋管截面流線圖,在靠近上壁面ab和下壁面ef的流體受黏滯力作用速度相對較小,而中心處的流速較大。因流速大的微粒離心慣性大,所以沿著中心線dc區域流線的壓強增量將大于ab和ef線的壓強增量。因此c處的壓強將大于a和e處的壓強,流體自c處流向a和e處,而d處的流體也將靠離心慣性不斷地向c處流去[12]。這樣,d處的壓強將小于b處和f處的壓強,二次流變在徑向平面內發生,其旋轉方向如圖1所示。這兩個旋轉運動與主流相結合,便產生了雙螺旋流動。綜合截面上的速度分量的特性與主流速度對二次流強度的影響,本文提出使用螺旋強度he來計算局部二次流強度,其計算方法如下:式中:Ω為局部渦量;U為速度矢量;渦量與速度矢量的點積,代表了其與速度的協同程度,取絕對值是因為無論二次流在截面上的流動是順時針還是逆時針,對傳熱的影響是一樣的。由于二次流強度受速度、密度、特征長度、黏度影響,所以,he為速度、密度、特征長度、黏度的函數,表示為he=f(U,ρ,d,μ)。(3)式中5個物理量的量綱均由3個基本量綱(時間量綱、長度量綱、質量量綱)組成,故可以組成2個無量綱數。選定特征長度、密度、黏度為基本物理量。在此只分析包含he的無量綱數H。設:H=he?daρbμc。(4)由量綱和諧原理得dimH=M0.5+a-3b-cT-1-cLb+c。(5)即a=1.5,b=1,c=-1。所以,代表著局部二次流強度的無量綱數H表為描述截面上的平均二次流強度,使用截面局部二次流強度H的面積平均值Ha,其計算式為∑dsi。(7)式中:i代表第i個網格節點;ds為節點在截面上的面積。由局部的二次流強度的無量綱數,通過對流域的體積積分并除以總體積,即可得到描述全場的二次流強度,即H的體積平均值:Hm=∑HidVi∑dVi。(8)為驗證本文所提出的局部二次流強度計算方法的準確性,以螺旋管內流體為對象,驗證不同參數下的Hm與De的相關性。兩者的相關性系數計算式為式中:Hm和De分別為所有樣本的Hm和De的平均值;j代表第j個樣本。

2數值模擬驗證

2.1幾何模型

螺旋管的幾何參數如圖2所示。d為螺旋管的內徑,D為螺旋管的纏繞直徑,P為螺距。L為螺旋管長度,在本文中所有模型的螺旋管長度保持一致。

2.2邊界條件

入口為速度入口,出口為壓力出口,壁面為靜止、無滑移,恒壁溫,不計壁厚。

2.3求解方法

本文數值模擬使用的求解器是FLUENT2019。模型為層流穩態模型。速度壓力耦合計算采用SIMPLE算法。同時設置收斂殘差標準為1×10-6,對出入口處質量流量以及壓力差進行監控,最終確定模擬達到收斂。

2.4網格無關性驗證

計算域的網格劃分方式采用六面體結構化網格,并對近壁面的網格進行了局部加密。網格無關性驗證選取6套網格尺寸不同的方案,以網格最小的方案為標準,計算其他5套方案的努塞爾數與壓降的相對誤差,結果如圖3所示。綜合考慮準確性與效率,選取第3套網格劃分方案(2190617)。2.5數值模擬方法驗證為驗證數值模擬方案的可靠性,本文以水為介質,在d為8mm、D為110mm、P為50mm,雷諾數為1000~5000的工況下將模擬值與文獻[16]中的公式(見式(18),適用范圍為500<De<5000)的計算結果做對比,結果如圖4所示。模擬值與實驗擬合公式得出的結果誤差在8%以內,說明數值模擬方案可行。

3結果與討論

3.1不同參數下局部二次流強度計算方法驗證

采用d為8mm、D為110mm、P為50mm的螺旋管結構,驗證不同雷諾數下De與Hm的相關性,結果如圖5所示。隨著雷諾數的增大,流體的主流方向速度增大,受到的離心力也增大,從而在徑向上的分速度增大,導致了二次流強度增加。在不同的雷諾數下,De與Hm都隨Re的增大而增大。進一步分析De與Hm的數值特征,兩者幾乎成正比關系,Hm/De在2.71~2.80。說明在不同雷諾數下,H的體積平均值能夠很好地描述整體二次流的強度,用H做為局部二次流強度指標是合理的。若螺距無窮大,此時螺旋管近似為直管,二次流不會發生,De數為零。但是在實際工程應用中,螺距的變化范圍有限,對二次流的影響相對很小。為驗證不同螺距下De數與Hm的關系,選取d為8mm、D為110mm、P為20~140mm(無量綱螺距λ為0.0579~0.4051)的螺旋管結構,雷諾數恒定為4000,結果如圖6所示。因螺距的增加導致了螺旋管整體的彎曲程度減小,流體受到的離心力下降,De與Hm都隨無量綱螺距的減小而減小,但是其減少的幅度很小。Hm/De在2.74~2.82。為驗證不同曲率下De數與Hm的關系,選取d為8~12mm(曲率為0.0727~0.1091)、D為110mm,P為50mm的螺旋管結構,雷諾數恒定為4000,結果如圖7所示,De與Hm都隨曲率的增加而增加,Hm/De在2.68~2.77。若局部二次流強度指標H能夠精準描述局部二次流強度,則依據公式(9),理論上H的體積平均值Hm與De數的相關系數應為1。但是由于舍入誤差的存在,所有樣本點的Hm與De數的相關性系數r為0.9949,十分接近于1,且Hm/De都在2.75左右,說明Hm和De數存在著較好的正比關系,同時也說明H作為螺旋管內局部二次強度指標是可行的。

3.2二次流演變規律

以水為工質,選取d為8mm、D為110mm、P為50mm、雷諾數恒定為4000的螺旋管結構為研究對象。圖8展示的是螺旋管中各個截面上的Hm。其中S表示的是截面中心沿螺旋中心線距入口的沿程距離。從圖8可以看出,流體依據二次流的強度變化趨勢,其流動狀態可以分為3個階段:發展段(入口~35mm)、過渡段(35mm~200mm)、穩定段(200mm~出口)。在發展段,流體從入口進入螺旋管后,在離心力的作用下迅速發展出二次流,二次流強度迅速增大,并在路程S=35mm時達到極大值。然后流體進入了過渡段,二次流強度在路程S=160mm處達到最大值,之后二次流強度略有下降。在穩定段,流體在截面上的平均二次流強度在2850上下呈小幅度規律性波動。圖9所示為各個截面上的局部二次流強度云圖。二次流強度沿著水平線呈對稱分布,二次流強度最大區域分布在上下壁面附近,管道中心區域二次流強度相對很小,并且局部二次流強度最大值在發展階段和過渡階段逐漸從內側向外側轉移(各云圖的右側為內側,左側為外側)。

3.3不同參數對二次流演變規律的影響

在不同的幾何參數或者工況下,螺旋管中的二次流強度和演變規律是不同的。同一根螺旋管在不同雷諾數下,流道各截面上的H的面積平均值分布如圖10所示。隨雷諾數的增大,螺旋管中的二次流在發展段增大得更快。不同雷諾數下,螺旋管中的流體都經歷3個階段,且3個階段都發生相同的位置,說明雷諾數只影響二次流的強度,不影響二次流的演變規律。圖11所示為不同曲率下螺旋管中流體的面積平均二次流強度隨沿程距離的變化圖。各管中的流體都經歷3個階段,但是強度和3個階段發生的位置都不同。具體來說,隨著曲率的增大,二次流的強度會略微增大,并且到達穩定階段所需的路程更長。圖12所示為不同螺距的螺旋管中流體的面積平均二次流強度隨沿程距離的變化圖。5條曲線基本重合,說明螺距對二次流的強度和演變規律影響非常小,隨螺距的增大,二次流強度略微減小,3個二次流階段發生的位置基本一致。

4結論

(1)無量綱螺旋強度H的體積平均值與迪恩數De相關系數約為1,H/De在2.75左右,H能夠準確地描述螺旋管內的局部二次流強度。(2)螺旋管內流體從垂直于入口流入到流出,依據二次流的強度可以劃分為3個階段,分別為發展段、過渡段和穩定段。二次流在發展段迅速增強,在過渡段出現最大值,在穩定段呈周期性波動。(3)雷諾數的增大會顯著增大二次流強度,但是不影響二次流3個階段的出現位置。曲率的增大會略微增大二次流的強度,且延長二次流達到穩定階段的路程。螺距對二次流的強度和發展規律影響非常小。

作者:吳金星 徐耀 李松歌 肖嘉邦 馬宇翔 單位:鄭州大學 機械與動力工程學院