基于停留時間分布的纏繞管內(nèi)二次流研究

黃正梁，王超，郭燕妮，楊遙，孫婧元，王靖岱，陽永榮

（1浙江大學(xué)化學(xué)工程與生物工程學(xué)院，浙江杭州310027；2浙江大學(xué)化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗(yàn)室，浙江杭州310027；3浙江省化工高效制造技術(shù)重點(diǎn)實(shí)驗(yàn)室，浙江杭州310027）

引 言

纏繞管換熱器廣泛應(yīng)用于石油化工、空氣分離、食品加工等領(lǐng)域[1?6]。與普通列管式換熱器相比，纏繞管換熱器特殊的緊湊結(jié)構(gòu)提供了更大的傳熱面積[7?10]，纏繞管內(nèi)流體在離心力的作用下形成的二次流增加了管內(nèi)流體的混亂程度，使傳熱速率顯著提高，管內(nèi)污垢熱阻更小[11?14]。因此，掌握纏繞管內(nèi)二次流的形成機(jī)制及其對管內(nèi)結(jié)垢的影響規(guī)律，不僅是纏繞管換熱器的設(shè)計基礎(chǔ)，也一直是該領(lǐng)域的研究熱點(diǎn)[15?16]。

纏繞管內(nèi)二次流的研究方法包括墨水染色直接觀察[17]、壓降建模間接表征[18]和流體力學(xué)模擬計算[19]等。染色法在早期使用較多，可直接觀察透明纏繞管中的二次流。研究者采用染色法發(fā)現(xiàn)了纏繞管中層流與湍流的轉(zhuǎn)變行為與直管不同，纏繞直徑和管徑對纏繞管中流動具有顯著影響[20?22]。染色法只能用于定性分析，不能提供定量的數(shù)據(jù)。壓降可以反映纏繞管中的流動狀態(tài)[23?24]，有學(xué)者考察了纏繞直徑、螺距、纏繞角度等結(jié)構(gòu)參數(shù)對纏繞管壓降的影響規(guī)律，建立了纏繞管壓降的經(jīng)驗(yàn)預(yù)測模型，通過與直管壓降預(yù)測模型的對比，間接表征二次流的影響[24?28]。壓降建模間接表征法不能揭示二次流的形成原因及其對壓降的影響機(jī)制。近年來，一些研究者通過CFD 模擬研究纏繞管內(nèi)的流動狀態(tài)、二次流及其對傳熱的影響，證實(shí)了纏繞管和直管內(nèi)的流動規(guī)律不同[29?35]。然而，由于缺乏纏繞管中二次流的檢測方法和數(shù)據(jù)，不能對二次流的模擬結(jié)果直接進(jìn)行驗(yàn)證[31,36]。

彎管中流體的流速分布不均、分子擴(kuò)散和湍流擴(kuò)散等因素，使得流體經(jīng)過彎管時存在液體停留時間分布[37]。液體停留時間分布可以反映彎管中流體偏離理想流動的程度。有學(xué)者針對微流控設(shè)備，在纏繞管入口階躍式注入示蹤劑，測量出口液體停留時間分布，研究結(jié)果表明，相比其他流動形式，二次流可以減少纏繞管中的流體動力擴(kuò)散[38]。也有學(xué)者通過脈沖式注入示蹤劑的方法測量螺旋盤管和直管的液體停留時間分布，證明了螺旋盤管可以得到更好的停留時間分布[39]。受此啟發(fā)，本文提出采用電導(dǎo)率法測量纏繞管內(nèi)流體停留時間分布，利用無量綱方差表征二次流強(qiáng)度的方法。采用該方法研究了纏繞直徑、纏繞角度、纏繞管管徑等結(jié)構(gòu)參數(shù)對纏繞管內(nèi)二次流的影響規(guī)律。

1 實(shí)驗(yàn)裝置及方法

如圖1所示，實(shí)驗(yàn)裝置由水槽、離心泵、流量計、纏繞管、電導(dǎo)率儀以及計算機(jī)等構(gòu)成。如圖2所示，實(shí)驗(yàn)中采用了7個長度均為4.6 m 的不銹鋼纏繞管。纏繞管的結(jié)構(gòu)參數(shù)見表1，實(shí)驗(yàn)介質(zhì)為水，經(jīng)離心泵增壓后從水箱進(jìn)入纏繞管內(nèi)，其流量由轉(zhuǎn)子流量計和閥門控制，在纏繞管出口檢測電導(dǎo)率隨時間的變化。纏 繞 管 中 的 流 速 變 化 范 圍 為0.2～1.6 m·s?1，Reynolds數(shù)的變化范圍為2160～39600。

圖2 纏繞管結(jié)構(gòu)示意圖Fig.2 Geometric parameters of a helical coil

表1 纏繞管的結(jié)構(gòu)參數(shù)Table 1 Structural parameters of helical coils

以飽和KCl 溶液作為示蹤劑，采用脈沖示蹤法和電導(dǎo)率儀測量纏繞管內(nèi)的液體停留時間分布曲線。實(shí)驗(yàn)步驟如下：（1）調(diào)節(jié)纏繞管內(nèi)液體流量至指定流量并穩(wěn)定10 min；（2）在t=0 時刻，將一定體積（5 ml）的飽和KCl 溶液通過注入裝置加到纏繞管；（3）在纏繞管出口使用電導(dǎo)率儀實(shí)時記錄電導(dǎo)率值。隨著KCl 溶液不斷流出，纏繞管出口處液體的電導(dǎo)率k(t)先上升，達(dá)到最大值kmax后開始下降，當(dāng)所有的KCl 溶液從纏繞管出口流出時，k(t)趨于穩(wěn)定并接近初始值，記為kmin；（4）改變實(shí)驗(yàn)條件，繼續(xù)進(jìn)行實(shí)驗(yàn)。

為了消除背景流體（水）的電導(dǎo)率對實(shí)驗(yàn)結(jié)果的影響，對k(t)進(jìn)行歸一化處理，得到無量綱的電導(dǎo)率k′(t)。

由k′(t)可計算得到相應(yīng)的停留時間分布密度函數(shù)E(t)：

對E(t)進(jìn)行分析得到平均停留時間-t：

2 實(shí)驗(yàn)結(jié)果與討論

2.1 直管中σ2隨Re的變化

無量綱方差可以描述管內(nèi)流體的混亂程度，在理想的平推流反應(yīng)器中，σ2=0，說明管截面上的速度幾乎完全相等；在全混流反應(yīng)器中，σ2=1，說明管截面上的速度無序分布。因此，本文測量不同Re下直管（d=15 mm,長度l=4.6 m）中的液體停留時間分布，作為對比例。利用式（1）～式（5）計算σ2，σ2與液體Re 的關(guān)系如圖3 所示，隨著Re 的增加，σ2逐漸降低并趨向于平穩(wěn)。在本文的實(shí)驗(yàn)條件下，直管中Re的變化范圍是3000～20000，表明管中流體均處于湍流狀態(tài)，傳質(zhì)受到對流和湍流擴(kuò)散的雙重影響，導(dǎo)致直管出口截面的σ2隨Re 的增加而降低并逐漸趨向于0。這與文獻(xiàn)所示的直管中流體流動的基本規(guī)律[23]一致，證實(shí)了實(shí)驗(yàn)裝置和方法的可靠性。

圖3 直管出口σ2隨Re的變化Fig.3 Variation of σ2 with Re in the straight pipe(d=15 mm,l=4.6 m)

2.2 不同結(jié)構(gòu)纏繞管中σ2隨Re的變化

利用上述實(shí)驗(yàn)方法，對不同Re下具有不同結(jié)構(gòu)參數(shù)的纏繞管測量其RTD 曲線，并計算纏繞管出口處的無量綱方差σ2。在相同纏繞角度（α）和管徑（d）下，不同纏繞直徑（D）的纏繞管出口σ2隨Re 的變化如圖4 所示。隨著Re 的增大，σ2呈先減小再增加的趨勢，這與直管有較大的差異。當(dāng)Re小于臨界Reynolds 數(shù)ReS（圖4 中曲線極小值點(diǎn)）時，σ2隨著Re的增大而減小，這與直管中的結(jié)果一致，表明在Re較小時，纏繞管內(nèi)的流體仍主要受與直管中類似的湍流效應(yīng)影響。隨著Re的增加，σ2呈下降趨勢并逐漸趨向于0，此時湍流效應(yīng)的影響逐漸減弱，管內(nèi)流體流動逐漸趨于理想平推流。當(dāng)Re 大于ReS時，管內(nèi)流體在離心力的作用下產(chǎn)生較強(qiáng)的二次流，導(dǎo)致管截面速度的徑向差異增大，管內(nèi)流動逐漸遠(yuǎn)離理想平推流，從湍流作用區(qū)轉(zhuǎn)變至二次流作用區(qū)。在二次流作用區(qū)，隨著Re 的繼續(xù)增大，管內(nèi)流體的離心力增大，管內(nèi)二次流逐漸增強(qiáng)，對應(yīng)的σ2逐漸增大。這與Greenspan[17]的研究結(jié)果一致，然而文獻(xiàn)中并沒有發(fā)現(xiàn)在足夠低的Re 下存在與直管中一致的湍流效應(yīng)。

圖4 不同纏繞直徑的纏繞管出口σ2隨Re的變化Fig.4 Variation of σ2 with Re in helical coils with different coiling diameters(α=10°,d=15 mm,l=4.6 m)

從圖4還可以看出，當(dāng)纏繞管的管徑相同時，在Re 較低的區(qū)域，湍流效應(yīng)對流體流動的影響差異較小，因此相同Re下測得的不同纏繞管出口σ2沒有明顯差異。在Re較高的二次流作用區(qū)，管內(nèi)流體微元受到的離心力與纏繞直徑呈反比，即纏繞直徑越小，流體受離心力的影響越大，二次流強(qiáng)度越大，因此σ2越大，ReS越小。其中，纏繞直徑為219、273 和325 mm 的 纏 繞 管 的ReS分 別 為10500、12000 和12000。這與Jayakumar 等[19]模擬得到的纏繞直徑對纏繞管傳熱系數(shù)的影響規(guī)律一致。他們研究發(fā)現(xiàn)二次流會促使管內(nèi)流體強(qiáng)烈再循環(huán)和混合，使得對流傳熱系數(shù)增加；同時，隨著纏繞直徑的增加，二次流的影響逐漸減弱，對流傳熱系數(shù)減小?？梢?，纏繞直徑越小，二次流對纏繞管傳熱的影響越大。

圖5 不同纏繞角度的纏繞管出口σ2隨Re的變化Fig.5 Variation of σ2 with Re in helical coils with different coiling angles(D=325 mm,d=15 mm,l=4.6 m)

在相同纏繞直徑（D）和管徑（d）下，不同纏繞角度（α）的纏繞管出口σ2隨Re 的變化如圖5 所示。與圖4 類似，不同纏繞角度的纏繞管出口σ2隨Re 的增加也呈現(xiàn)先減小再增加的變化規(guī)律。在二次流作用區(qū)，不同纏繞管中的二次流強(qiáng)度在相同Re下較為接近，即出口σ2幾乎相同，這與Jayakumar 等[19]的傳熱模擬結(jié)果具有一致性。他們研究了纏繞角度對Nusselt 數(shù)的影響，發(fā)現(xiàn)當(dāng)纏繞角度在8°～10°的范圍內(nèi)時，局部Nu 幾乎相同，此時纏繞角度不會影響二次流的強(qiáng)度?？梢姡p繞角度對纏繞管傳熱效率的影響較小，可以忽略。

在相同纏繞直徑（D）和纏繞角度（α）下，不同管徑（d）的纏繞管出口σ2隨Re 的變化如圖6 所示。當(dāng)纏繞管的管徑為11 mm 時，σ2隨Re 的增加先減小后增大再趨于平穩(wěn)；當(dāng)管徑為15 mm和22 mm時，σ2隨Re 的增加先減小后增大，與圖4 和圖5 所示的規(guī)律相同。當(dāng)纏繞管的管徑為11、15 和22 mm 時，對應(yīng)的ReS分別為5500、12000 和22000?？梢姡p繞管的管徑越小，ReS越小。這與Mirgolbabaei 等[34]的傳熱模擬結(jié)果一致，他們研究發(fā)現(xiàn)當(dāng)管徑增大時，纏繞管的換熱效率降低，壓降也隨之減小。當(dāng)纏繞管的纏繞直徑和角度相同時，單位流體微元在徑向上受到的離心力大小相同，即流體微元的徑向加速度相同，在相同的流動時間內(nèi)其在管道徑向截面上的位移也相同。因此，纏繞管的管徑越小，徑向上流體相對位移越大，二次流作用越強(qiáng)。但是，流體在徑向上的遷移距離受到管徑的限制。在小直徑的纏繞管中，當(dāng)流體徑向遷移距離達(dá)到極限時，纏繞管出口處的σ2不再隨Re變化，對應(yīng)的二次流強(qiáng)度達(dá)到極限，此時管內(nèi)流體流動從二次流作用區(qū)轉(zhuǎn)變?yōu)槎瘟鳂O限區(qū)。由于實(shí)驗(yàn)條件的限制，本文僅在管徑11 mm的纏繞管中發(fā)現(xiàn)了二次流極限區(qū)。

圖6 不同管徑的纏繞管出口σ2隨Re的變化Fig.6 Variation of σ2 with Re in helical coils with different tube diameters(D=325 mm,α=10°,l=4.6 m)

2.3 纏繞管中流體流動的轉(zhuǎn)變規(guī)律

不同結(jié)構(gòu)參數(shù)的纏繞管中臨界Reynolds數(shù)如表2 所示，其中ReS是從湍流作用區(qū)轉(zhuǎn)變?yōu)槎瘟髯饔脜^(qū)的臨界Reynolds 數(shù)，ReL是從二次流作用區(qū)轉(zhuǎn)變?yōu)槎瘟鳂O限區(qū)的臨界Reynolds 數(shù)。由表可知，ReS隨著纏繞直徑（D）和纏繞管管徑（d）的減小而減小，纏繞角度（α）對ReS沒有影響。根據(jù)上述實(shí)驗(yàn)結(jié)果，構(gòu)建了如圖7 所示的纏繞管內(nèi)流動特性示意圖。在不同的操作條件下，纏繞管內(nèi)會出現(xiàn)湍流作用區(qū)、二次流作用區(qū)和二次流極限區(qū)三種不同的流型。其中，Ⅰ區(qū)域?yàn)橥牧髯饔脜^(qū)，纏繞管出口σ2隨Re的變化規(guī)律與直管類似，纏繞管出口σ2隨Re 的增大而降低；Ⅱ區(qū)域?yàn)槎瘟髯饔脜^(qū)，纏繞管出口σ2隨Re 的增大而增大；Ⅲ區(qū)域?yàn)槎瘟鳂O限區(qū)，纏繞管出口σ2隨Re 的增大不再明顯變化。當(dāng)纏繞管內(nèi)的徑向流動未達(dá)到極限遷移距離時，隨著Re 的增大，纏繞管中流體流動特征會從湍流作用區(qū)轉(zhuǎn)變?yōu)槎瘟髯饔脜^(qū)；當(dāng)纏繞管內(nèi)的徑向流動達(dá)到極限遷移距離時，纏繞管中的流體流動特征會從湍流作用區(qū)先轉(zhuǎn)變?yōu)槎瘟髯饔脜^(qū)，再轉(zhuǎn)變?yōu)槎瘟鳂O限區(qū)。

3 結(jié) 論

（1）提出了通過測量纏繞管內(nèi)液體的停留時間分布，利用無量綱方差σ2定量表征纏繞管內(nèi)二次流強(qiáng)度的方法。

（2）考察了纏繞直徑、纏繞角度、管徑對纏繞管內(nèi)二次流的影響規(guī)律。當(dāng)纏繞管內(nèi)的徑向流動未達(dá)到極限遷移距離時，纏繞管出口σ2隨著Re的增大先減小后增大，分別對應(yīng)湍流作用區(qū)和二次流作用區(qū)，轉(zhuǎn)變點(diǎn)為臨界Reynolds 數(shù)ReS；當(dāng)纏繞管內(nèi)的徑向流動達(dá)到極限遷移距離時，纏繞管出口σ2隨著Re的增大先減小后增大再趨于平穩(wěn)，分別對應(yīng)湍流作用區(qū)、二次流作用區(qū)和二次流極限區(qū)，轉(zhuǎn)變點(diǎn)分別為臨界Reynolds數(shù)ReS和ReL。

表2 不同纏繞管的臨界ReS和ReLTable 2 ReS and ReL in different helical coils

圖7 纏繞管中流體流動特性示意圖Fig.7 Schematic diagram of the fluid flow characteristics in helical coils

（3）從湍流作用區(qū)轉(zhuǎn)變?yōu)槎瘟髯饔脜^(qū)的臨界Reynolds 數(shù)ReS隨纏繞直徑和纏繞管管徑的減小而減小，纏繞角度對ReS的影響較小。