中空扭旋元件排列方式對管內(nèi)層流傳熱特性的影響

張傲宇,丁桂彬,王宗勇

(沈陽化工大學(xué) 能源與動力工程學(xué)院， 遼寧 沈陽 110142)

近年來，隨著人們環(huán)保意識的增強(qiáng)，提高能源利用率、節(jié)約能源越來越受人們關(guān)注.換熱器作為化工生產(chǎn)過程中一種必不可少的設(shè)備，減小換熱器體積、提高其換熱效率顯得尤為重要.為了達(dá)到換熱器結(jié)構(gòu)緊湊、換熱高效的目的，可以采取的最重要方法就是增加傳熱系數(shù)和減小壓降.傳熱系數(shù)的增加降低了熵產(chǎn)[1]；減小壓降可以提高能源的利用率.在圓管中插入扭帶是增強(qiáng)傳熱效果最常用的方法之一[2-5]，其主要原因是扭帶使管內(nèi)流體產(chǎn)生旋轉(zhuǎn)并引起二次流，促進(jìn)流體的徑向混合，在離心力影響下，管中心流體和壁面邊界層流體充分混合，形成一個溫度比較均勻、速度變化較為明顯的核心區(qū)域，同時在換熱壁面附近形成一個溫度和速度變化均比較顯著的邊界層區(qū)域，從而達(dá)到強(qiáng)化傳熱效果[6].

內(nèi)置扭帶的傳熱強(qiáng)化方法雖然能有效提高換熱器的傳熱效率，但是同時也增加了流動阻力，使換熱器能耗增加.Smith Eiamsa-ard等[7]對內(nèi)置帶中心桿的扭帶、不帶中心桿的扭帶和規(guī)則間隔扭帶3種渦流發(fā)生器進(jìn)行了強(qiáng)化傳熱研究，發(fā)現(xiàn)各種類型扭帶均可有效提高換熱效率；間隔扭帶在有效提升換熱效率的同時也能將壓降增長幅度控制在一個較小的范圍內(nèi).張曉嶼等[8]對內(nèi)置中空扭帶換熱管的傳熱與流動特性進(jìn)行了數(shù)值模擬研究，結(jié)果表明扭旋元件表面積的縮小能有效降低流體的阻力，并且能夠顯著提升換熱管的換熱綜合性能.此外，張曉嶼等[9]對在圓管內(nèi)間隔布置多個扭旋元件時的傳熱和流動性能進(jìn)行了數(shù)值模擬研究，研究結(jié)果表明：在充分發(fā)展的層流狀態(tài)下，間隔扭旋元件可以造成圓管傳熱性能的強(qiáng)化，并且流動阻力增長不大，從而提高換熱管的熱-水力學(xué)綜合性能.上述研究多基于元件連續(xù)排列方式，并未對元件交錯排列方式做出研究.王特等[10]對內(nèi)置多種錯位角的交替軸扭帶的管內(nèi)層流換熱特性進(jìn)行了數(shù)值模擬研究，結(jié)果顯示元件交錯排列不僅可以引起漩渦流，而且在元件交錯部分可以引導(dǎo)流體向兩側(cè)分流，引起流體擾動，促使流體更高效地混合，從而獲得更好的換熱效果.

郭劍等[11]對十字形扭帶的傳熱強(qiáng)化效果進(jìn)行了數(shù)值模擬，其研究結(jié)果表明若在受限空間的核心區(qū)添加擾流物，可使核心區(qū)的流體溫度趨于均勻，若擾流物所占空間比較小，不足以產(chǎn)生較大的阻力，則可以實(shí)現(xiàn)傳熱強(qiáng)化的同時阻力增加不大，從而提高換熱設(shè)備的綜合傳熱性能.Guo等[12]分別對內(nèi)置傳統(tǒng)扭帶熱管與中空扭帶熱管的傳熱性能進(jìn)行了數(shù)值模擬研究，結(jié)果表明與傳統(tǒng)扭帶相比，中空扭帶傳熱性能提高了7%～20%.Li等[13]通過數(shù)值模擬研究了一種新型窄邊中空扭帶的強(qiáng)化傳熱特性，與傳統(tǒng)扭帶相比，該新型窄邊中空扭帶的最佳整體傳熱性能提高了28.1%.當(dāng)扭帶數(shù)為4、雷諾數(shù)增加到600以上時，換熱管的整體換熱性能最佳.

縱觀現(xiàn)有文獻(xiàn)可知，中空扭旋元件與傳統(tǒng)扭帶相比具有良好的綜合傳熱性能，得到了中空扭旋元件結(jié)構(gòu)包括元件寬度、扭率等對傳熱性能的影響規(guī)律，但關(guān)于中空扭旋元件排列方式相關(guān)研究鮮有報道，這種現(xiàn)狀阻礙了該種類型換熱器的優(yōu)化設(shè)計(jì)及工業(yè)應(yīng)用.因此，筆者旨在通過探究中空扭旋元件連續(xù)和交錯兩種排列方式對管內(nèi)傳熱特性及流動特性的影響規(guī)律，期望為中空扭旋元件的進(jìn)一步開發(fā)利用提供理論基礎(chǔ)和技術(shù)支持.

1 物理模型及數(shù)值模擬方法

1.1 物理模型

筆者以內(nèi)置中空扭旋元件換熱管為研究對象，其中扭旋元件采用交錯和連續(xù)兩種排列方式.交錯排列中空扭旋元件由4組相同的扭旋元件組成，各組扭旋元件沿?fù)Q熱管軸向首尾相接排布，并同心內(nèi)置于管內(nèi)，相鄰各組扭旋元件錯開90°，每組扭旋元件包含2塊圍繞在同一圓柱面上的扭旋片，兩塊扭旋片成對稱布置，各扭旋片扭轉(zhuǎn)角度為180°.連續(xù)排列中空扭旋元件與交錯排列中空扭旋元件唯一不同之處在于相鄰各組元件交錯角為0°，其結(jié)構(gòu)相當(dāng)于將連續(xù)扭帶用圓筒截去中心部分.

兩種排列方式的換熱管結(jié)構(gòu)如圖1所示.換熱管結(jié)構(gòu)參數(shù)為：管道內(nèi)徑D=20 mm；管長L=200 mm；管道入口長度Li=40 mm；元件總長Le=120 mm；扭旋板厚度δ=1 mm；扭旋元件外徑do=19 mm，內(nèi)徑di=11 mm.

圖1 內(nèi)置交錯和連續(xù)中空扭旋元件的換熱管結(jié)構(gòu)示意圖

1.2 數(shù)值模擬方法

模擬過程以水為工質(zhì)，為貼近實(shí)際工況，采用初始溫度323 K，恒壁溫邊界條件(壁溫為363 K)，雷諾數(shù)Re=200～1800(Re=ρuD/μ).元件與流體接觸的邊界采用無滑移邊界條件，忽略自然對流及溫度對介質(zhì)性質(zhì)的影響.采用Fluent 17.0進(jìn)行CFD模擬，選用三維穩(wěn)態(tài)模型，壓力場和速度場的耦合采用Simplec算法，動量和能量方程均采用QUICK格式，連續(xù)性、動量和能量方程收斂條件均設(shè)定在10-6以下.

描述流體流動的封閉控制方程包括連續(xù)性方程、動量方程以及能量方程.

(1)

(2)

(3)

其中:ρ為密度，kg/m3；v為流體平均速度，m/s；μ為流體黏度，Pa·s;vi,vk為i，k(i，k=1,2,3)方向的速度分量,m/s;p為壓強(qiáng)，Pa；T為流體溫度，K；λ為導(dǎo)熱系數(shù)，W/(m·K);Cp為流體比熱容，J/(kg·K).

1.3 網(wǎng)格劃分及無關(guān)性檢驗(yàn)

筆者采用Solidworks建立幾何模型，利用hypermesh進(jìn)行網(wǎng)格劃分，采用四面體網(wǎng)格并在靠近壁面處進(jìn)行網(wǎng)格加密.網(wǎng)格質(zhì)量得到提升，縱橫比(aspect ratio)達(dá)5以下，雅可比值(Jacobian)達(dá)0.9以上，網(wǎng)格質(zhì)量較好，滿足模擬需要.網(wǎng)格局部加密如圖2所示.為消除網(wǎng)格尺寸對計(jì)算結(jié)果的影響，對含有連續(xù)中空扭旋元件結(jié)構(gòu)模型的換熱管在Re=200情況下進(jìn)行網(wǎng)格無關(guān)性檢驗(yàn).劃分5種尺寸網(wǎng)格，得到每種網(wǎng)格數(shù)量模型下的平均努賽爾數(shù)，如圖3所示.由圖3可以看出：當(dāng)換熱管的網(wǎng)格數(shù)為147萬左右時網(wǎng)格尺寸對傳熱系數(shù)的影響即可忽略不計(jì).

圖2 邊界層網(wǎng)格加密

圖3 換熱管網(wǎng)格獨(dú)立性驗(yàn)證

2 傳熱強(qiáng)化性能評價方法

2.1 綜合傳熱性能評價因子(PEC值)

換熱管的綜合傳熱性能由傳熱性能和阻力性能綜合體現(xiàn).傳熱性能由努賽爾數(shù)Nu表征;阻力性能由流動阻力系數(shù)f表征，其表達(dá)式為

(4)

式中：Δp為壓力降，Pa；ρ為流體密度，kg/m3；vb為流體平均速度，m/s.

綜合傳熱性能評價因子PEC值能夠考察換熱器整體特性的優(yōu)劣，其表達(dá)式為

PEC=(Nu/Nu0)/(f/f0)1/3.

(5)

式中:Nu0為對應(yīng)光管的努塞爾數(shù);f0為對應(yīng)光管的阻力系數(shù).

2.2 火積傳熱效率

Guo等[14]從導(dǎo)熱過程與導(dǎo)電過程的比擬出發(fā)，提出了表征物體傳熱能力的新物理量火積及火積耗散的概念，火積的提出使定義一個傳熱過程的效率成為現(xiàn)實(shí).火積耗散是物體傳遞熱量能力的體現(xiàn)，熱量傳遞是一個不可逆過程，在傳遞過程中部分火積將被耗散.

火積傳熱效率為傳熱性能的最直接的指標(biāo)，實(shí)現(xiàn)強(qiáng)化傳熱的最根本原因都?xì)w結(jié)為傳熱效率的提升.圖4表示了在恒壁溫條件下?lián)Q熱管內(nèi)流體的傳熱過程.整個過程的火積傳熱效率由下列方程求出：

(6)

式中：Q為壁面向流體傳遞的熱量，J/s；TW為壁面溫度，K；Ediss為火積耗散，J·K；Cv為定容比熱容，J/(kg·K)；m為質(zhì)量流量，kg/s；Tout為出口溫度，K；Tin為入口溫度，K.

Q=CvmTout-CvmTin.

(7)

(8)

整個過程的火積傳熱效率：

(9)

式中：η為火積傳熱效率；Ein為輸入火積，J·K；Eout為輸出火積，J·K.

圖4 恒壁溫邊界條件下管內(nèi)傳熱過程

3 數(shù)值模擬結(jié)果及分析

3.1 平均努塞爾數(shù)Nu

為了分析中空扭旋元件排列方式對傳熱性能的影響，筆者計(jì)算了Re=200～1800范圍內(nèi)共計(jì)9個雷諾數(shù)下的交錯排列、連續(xù)排列元件管以及光管Nu.Nu隨Re變化規(guī)律如圖5所示.由圖5可看出：3種結(jié)構(gòu)形式換熱管的Nu均隨著Re的增大而增大，且兩種內(nèi)置扭旋元件換熱管的努賽爾數(shù)明顯好于光管，相對于光管傳熱強(qiáng)化效果均隨著Re的增大而提高，在所分析的Re范圍內(nèi)連續(xù)和交錯扭旋元件換熱管的Nu分別為光管Nu的1.41～ 2.60倍和1.56～2.33倍.

當(dāng)Re<600時，交錯元件換熱管的Nu高于連續(xù)元件換熱管Nu.這是由于在較小的Re下，兩種排列形式元件區(qū)域范圍內(nèi)流體的平均流速相差微小，但是交錯元件換熱管由于相鄰各組元件存在交錯結(jié)構(gòu)，前一元件流出的流體會受到后一元件扭旋板的切割和分流作用，致使兩元件分界面處的近壁流體特別是傳熱邊界層內(nèi)流體產(chǎn)生擾動，傳熱邊界層減薄，傳熱系數(shù)提高；連續(xù)元件換熱管內(nèi)的近壁流體由于元件的連續(xù)結(jié)構(gòu)，不存在流體的局部擾動，因此傳熱系數(shù)偏低.兩種排列形式元件的努賽爾數(shù)在該Re范圍相差不大，說明2塊扭旋板的交錯分流引起的擾流作用影響有限，若增加扭旋板數(shù)量應(yīng)該能夠進(jìn)一步提高交錯元件的傳熱系數(shù).

當(dāng)600≤Re≤1800時，連續(xù)元件換熱管的Nu大于交錯元件換熱管Nu，并且兩者之間的差距隨著Re增大而增大.形成這種規(guī)律的原因在于：扭旋元件對流體既有旋流作用又有阻流作用，盡管交錯元件能夠?qū)α黧w產(chǎn)生局部擾動而促進(jìn)傳熱，但其對旋流的阻礙作用會隨著Re增大而明顯提高，導(dǎo)致扭旋元件區(qū)域流體旋流速度相對連續(xù)元件偏小.這種偏小程度會隨著Re的增大而增大，致使連續(xù)元件的旋流強(qiáng)化傳熱作用持續(xù)增強(qiáng)，明顯大于交錯元件的旋流強(qiáng)化作用，而交錯元件的擾流強(qiáng)化作用雖然隨著Re增大也會有所提高，但提高幅度極為有限，進(jìn)而導(dǎo)致連續(xù)元件的傳熱系數(shù)在該Re范圍內(nèi)要大于交錯元件，且隨著Re增加幅度在擴(kuò)大.

圖5 傳熱系數(shù)Nu變化曲線

3.2 阻力系數(shù)f

如圖6所示，在Re=200～1800時，連續(xù)元件換熱管的阻力系數(shù)f略大于交錯元件換熱管的阻力系數(shù)，但兩相之間相對偏差較小(<5%),特別是在Re<600范圍時，兩種排列形式換熱管的阻力系數(shù)幾乎完全相同.該Re范圍內(nèi)兩種元件排列方式換熱管的阻力系數(shù)分別為光管阻力系數(shù)的2.01～3.16倍、2.01～3.00倍.

連續(xù)元件換熱管的阻力系數(shù)略高于交錯元件換熱管的系數(shù)阻力，從直觀感覺上令人難以理解.為了便于理解這種流動現(xiàn)象，筆者繪制了Re=1800情況下兩種排列方式換熱管的橫截面速度云圖，如圖7所示.通過計(jì)算可得連續(xù)元件和交錯元件兩種換熱管在元件區(qū)域的平均軸向流速分別為0.88vm和0.94vm(vm為管內(nèi)平均軸向流速).由此可以看出：交錯元件區(qū)域范圍的流速并沒有因?yàn)樵诲e分流的原因致使流動阻力增大、流速下降，反而比連續(xù)元件流速有所提高.形成這種狀況的原因是交錯元件在分流的同時，在相鄰元件分界面處形成了徑向渦流，該渦流對流動邊界層產(chǎn)生了擾動，阻礙了其進(jìn)一步發(fā)展變厚，從而減小流動阻力，導(dǎo)致整個流動區(qū)域的壓力降與連續(xù)元件相比偏小，所以會形成交錯元件阻力系數(shù)小于連續(xù)元件阻力系數(shù)的狀況.

圖6 阻力系數(shù)f變化曲線

圖7 換熱管截面速度云圖

3.3 傳熱性能評價因子PEC值

傳熱性能評價因子PEC值能夠反應(yīng)傳熱與流體阻力的綜合影響，兩種內(nèi)置扭旋元件換熱管的PEC值變化曲線如圖8所示.由圖8可以看出：連續(xù)元件換熱管和交錯元件換熱管的PEC值均隨著Re的增大而增大且均大于1；在Re=1800時連續(xù)元件換熱管和交錯元件換熱管的PEC值取得最大值，分別為1.75和1.62，證明在空心軸上插入中空扭旋元件可以取得很好的換熱與流動綜合性能;在Re<600時，交錯元件換熱管的PEC值大于連續(xù)元件換熱管的PEC值,兩者的偏差程度隨著Re的增大而減小；當(dāng)600≤Re≤1800時，交錯元件換熱管的PEC值小于連續(xù)元件換熱管的PEC值，兩者的偏差程度隨著Re的增大而增大.兩種排列方式換熱管的PEC值隨Re的變化規(guī)律與傳熱系數(shù)Nu變化規(guī)律基本一致，這主要是由于兩種排列方式元件的摩擦系數(shù)差別不大導(dǎo)致的.

圖8 PEC值變化曲線

3.4 火積傳熱效率η

根據(jù)火積耗散極值原理作η-Re變化曲線，結(jié)果如圖9所示.

圖9 火積傳熱效率η變化曲線

由圖9可以看出：兩種換熱管的火積換熱效率η均隨著Re的增大而減小，連續(xù)元件換熱管η由90.54%降至92.17%，交錯元件換熱管η由90.48%降至92.20%；在所分析的Re范圍內(nèi)，連續(xù)元件換熱管的火積傳熱效率均高于交錯元件換熱管的火積傳熱效率，但差別不大；Re<980范圍內(nèi)，光管火積傳熱效率η高于兩種扭旋元件換熱管的火積傳熱效率，在Re>980范圍內(nèi)光管的火積傳熱效率則小于兩種扭旋元件換熱管的火積傳熱效率.由火積傳熱效率公式(9)可以看出，在流體入口溫度及管壁溫度恒定的情況下，傳熱效率僅取決于換熱管出口溫度.兩種排列方式元件換熱管的火積傳熱效率變化規(guī)律表明連續(xù)元件促進(jìn)了整個管內(nèi)流體的均布，促進(jìn)了徑向熱量傳遞，從火積傳熱效率角度應(yīng)優(yōu)先選用連續(xù)元件結(jié)構(gòu)，但從流動阻力角度考慮情況恰恰相反.

4 結(jié) 論

筆者對交錯和連續(xù)排列的中空扭旋元件換熱管傳熱性能進(jìn)行了數(shù)值模擬研究，并從傳熱系數(shù)、綜合傳熱性能評價因子和火積耗散極值原理3個方面對換熱管的傳熱性能進(jìn)行了評價.研究結(jié)論如下：

(1) 在Re=200～1800范圍內(nèi)，兩種排列方式的中空扭旋元件換熱管的Nu均隨著Re的增大而增大，且均高于光管的Nu；在較低Re(Re<600)范圍內(nèi)，交錯元件換熱管的Nu高于連續(xù)元件換熱管的Nu；在較高Re(600≤Re≤1800)范圍內(nèi)，連續(xù)中空扭旋元件換熱管的Nu較高.

(2) 在Re=200～1800范圍內(nèi)，交錯和連續(xù)元件排列方式換熱管的阻力系數(shù)f略高于交錯元件換熱管的阻力系數(shù)，但兩者相差不大，最大不超過5%.連續(xù)扭旋元件和交錯扭旋元件的阻力系數(shù)分別為光管阻力系數(shù)的2.01～3.16倍和2.01～3.00倍.

(3) 在Re=200～1800范圍內(nèi)，兩種排列方式元件換熱管的PEC值均高于1，在Re=1800情況下，連續(xù)元件和交錯元件換熱管的PEC值分別取得最大值，分別為1.75和1.62，說明兩種排列方式元件均能有效的提高換熱管的綜合傳熱性能.在低Re(Re<600)下，交錯元件換熱管的PEC值高于連續(xù)元件換熱管的PEC值，證明在低Re(Re<600)下，交錯元件換熱管的綜合傳熱性能優(yōu)于連續(xù)元件換熱管；在較高Re(600≤Re≤1800)范圍，連續(xù)扭旋元件換熱管的PEC值較高，證明在較高Re(600≤Re≤1800)下，連續(xù)元件換熱管的綜合傳熱性能更優(yōu).

(4) 根據(jù)火積散耗極值原理，在Re=200～1800范圍內(nèi)，連續(xù)扭旋元件換熱管的火積傳熱效率η基本均略高于交錯元件換熱管的火積傳熱效率，但兩種換熱管的火積傳熱效率差別不大.