朱 杰，潘 騰，王麗軍，成有為，李 希

浙江大學(xué)化學(xué)工程與生物工程學(xué)院，浙江 杭州 310027

隨著連續(xù)化、自動化和智能化生產(chǎn)觀念的深入，連續(xù)流反應(yīng)技術(shù)在精細化工和藥物化工領(lǐng)域快速發(fā)展[1-2]。作為一種典型的連續(xù)流反應(yīng)器，微反應(yīng)器憑借自身高效的傳熱和傳質(zhì)性能在有機合成反應(yīng)中得到廣泛應(yīng)用[3-4]。然而，微尺度在帶來優(yōu)越的傳遞性能的同時也限制了流動通量，微反應(yīng)器的大規(guī)模工業(yè)應(yīng)用仍有許多挑戰(zhàn)[5-6]。

微反應(yīng)器的放大策略主要有三種[7]：一是平行數(shù)目放大，該方法的優(yōu)點是可以保持傳遞性能不變，難點在于如何實現(xiàn)成千上萬臺設(shè)備的流量均布；二是串聯(lián)數(shù)目放大，該方法無須流量分配器，可以應(yīng)用于較長反應(yīng)停留時間的反應(yīng)以及多步反應(yīng)，但是能夠?qū)崿F(xiàn)的通量有限，鮮有超過100 mL/min處理量的報道；三是幾何尺寸放大，該方法可以實現(xiàn)通量的數(shù)量級增長，降低通道堵塞風(fēng)險，難點在于保留微尺度的傳遞性能?？祵幑净趲缀纬叽绶糯蟮牟呗栽O(shè)計了一種先進流反應(yīng)器（Advanced Flow Reactor，AFR）[8-10]，該反應(yīng)器的特征幾何尺寸為寬1 mm、高1.1 mm，流動通量達到10～200 mL/min，氣-液相界面積為160～1 300 m-1[11]，液-液相界面積為1 000～10 000 m-1[12]，僅需七個AFR組合就可實現(xiàn)每分鐘公斤級的原料處理量，展現(xiàn)出了一定的工業(yè)應(yīng)用潛力。Fang等[13]自行設(shè)計了一種可以連續(xù)合成季銨鹽的反應(yīng)器，該反應(yīng)器由內(nèi)徑為9 mm、總長為200 m的鈦材蛇形圓管構(gòu)成，管內(nèi)容積為15.7 L，在管內(nèi)插入組合式靜態(tài)混合器以增強傳熱和傳質(zhì)效果，結(jié)果表明該反應(yīng)器可以在幾分鐘內(nèi)實現(xiàn)相當(dāng)于間歇反應(yīng)器10 h的生產(chǎn)能力，年生產(chǎn)能力超過1 000 t。

本研究基于擴大通道尺寸并內(nèi)置三維螺旋構(gòu)件的思路設(shè)計了一種高通量連續(xù)流反應(yīng)器（High-throughput Continuous-flow Reactor，HTCR），該反應(yīng)器在保持較優(yōu)的混合與傳熱性能的前提下，通量比常見微反應(yīng)器增加1～3個數(shù)量級，為微反應(yīng)器的幾何放大提供了一條新的思路。本工作首先介紹HTCR的幾何構(gòu)型，再采用計算流體力學(xué)（CFD）的方法對反應(yīng)器基本單元的對流傳熱過程進行模擬，分析其流動與傳熱規(guī)律，隨后采用單因素分析的方法考察雷諾數(shù)、普朗特數(shù)、螺旋片長徑比和扭轉(zhuǎn)角對傳熱的影響，最后通過數(shù)值擬合得到努塞爾數(shù)與流動摩擦系數(shù)的計算關(guān)聯(lián)式。

1 反應(yīng)器幾何構(gòu)型設(shè)計

流體的旋轉(zhuǎn)、拉伸及剪切，對強化流體傳熱、混合和傳質(zhì)有重要作用。借鑒SK靜態(tài)混合器的結(jié)構(gòu)，提出了一種新的連續(xù)流反應(yīng)器型式，其基本單元為內(nèi)置SK螺旋片的六角直管。如圖1所示，六角管橫截面為正六邊形，SK螺旋片為左旋、右旋交替出現(xiàn)的周期性結(jié)構(gòu)，相鄰元件頭尾交叉角90°，將六角管與SK螺旋片組合在一起，即構(gòu)成HTCR的基本單元。

圖1 內(nèi)置SK螺旋片的六角管示意Fig.1 Schematic diagram of hexagonal tube embedded SK helical plates

將基本單元沿六個外壁面法向陣列再切割形成反應(yīng)管束，如圖2所示。將反應(yīng)管束軸向的兩個端面套上固定板，在板上按一定的規(guī)則開槽使每根六角管依次連通，連通后的槽口布置如圖3所示。將槽口表面封閉后再將z方向兩個端面用矩形板封閉，并在第一根六角管上開若干進料口，在最后一根六角管上開一個出料口，此時反應(yīng)器的流體通道構(gòu)建完畢，反應(yīng)器內(nèi)物料流向如圖4所示。

圖2 反應(yīng)管束示意Fig.2 Schematic diagram of reaction tubes bunch

反應(yīng)管束中六角管外的區(qū)域為換熱介質(zhì)區(qū)域，在x方向上兩個端面分別設(shè)置換熱介質(zhì)進、出口，在進口與反應(yīng)管束間設(shè)置分布器，各部件的裝配關(guān)系如圖5所示。

圖5 零部件裝配Fig.5 Assembly diagram of parts

反應(yīng)器運行時，反應(yīng)物分別通過不同的進料口進入第一根六角管，在螺旋片的作用下進行快速混合與高效傳熱，然后沿立體蛇形路徑流至最后一根六角管再從出料口流出。換熱介質(zhì)通過管線進入分布器，從分布器孔流入反應(yīng)管束殼程，通過間壁傳熱與管內(nèi)反應(yīng)物進行熱交換，兩流體錯流換熱。

2 數(shù)學(xué)模型

2.1 幾何模型

抽取流體通道的最小重復(fù)單元作為幾何模型，如圖6所示，l為一個螺片的長度，橫截面正六邊形邊長為5 mm，內(nèi)置n組螺旋片。

圖6 幾何模型Fig.6 Geometric model

2.2 控制方程

對于不可壓縮牛頓流體的三維穩(wěn)態(tài)流動，忽略質(zhì)量力，以張量形式表示的控制方程如下：

連續(xù)性方程

動量方程

能量方程

由于模擬的雷諾數(shù)（Re）在湍流區(qū)（對于SK螺旋結(jié)構(gòu)，Re＞150即為湍流區(qū)[14]）。根據(jù)以往對于三維剪切流的研究，剪切應(yīng)力傳輸模型（SST）k-ω湍流模型可以更準確地描述流體行為[15]，故本研究選用該模型進行計算。

2.3 網(wǎng)格劃分及無關(guān)性分析

采用軟件Star ccm+的網(wǎng)格劃分功能，貼近壁面設(shè)置5層棱柱層網(wǎng)格以解析邊界層，保證算例中壁面y+值為0.1～1.2，網(wǎng)格劃分情況如圖7所示。

圖7 網(wǎng)格劃分進口段(a)和SK段(b)Fig.7 Mesh of inlet section(a) and SK structure section(b)

網(wǎng)格無關(guān)性分析的邊界條件設(shè)置：進口速度為0.3 m/s，進口溫度為60 ℃；出口為壓力邊界，其值為0 Pa；幾何壁面為無滑移邊界條件，外壁面溫度設(shè)為等溫80 ℃，螺旋片壁面設(shè)置為絕熱，流動介質(zhì)設(shè)置為水。由疏到密選取了163 659，262 874，418102，566 706，866 350和14 490 616六個網(wǎng)格總數(shù)目進行計算，以模擬對象的努塞爾數(shù)（Nu）和摩擦系數(shù)（f）作為判斷依據(jù)，結(jié)果如圖8所示。結(jié)果顯示，f在整個考察范圍內(nèi)均保持穩(wěn)定，當(dāng)網(wǎng)格數(shù)大于56.6萬時，Nu幾乎不隨網(wǎng)格數(shù)的增長而變化，此時網(wǎng)格尺寸已經(jīng)可以消除其本身對計算結(jié)果的影響，本文案例均采用該尺寸進行計算。

圖8 網(wǎng)格無關(guān)性分析Fig.8 Grid independence analysis

為了排除進出口對計算結(jié)果的干擾，進一步考察螺旋片組數(shù)（n）對流動和傳熱的影響，模擬結(jié)果如圖9所示。當(dāng)n≥6時，Nu和f趨于平緩，說明此時螺旋片組數(shù)對于傳熱和流動阻力已經(jīng)幾乎沒有影響，故本模擬計算的n均以6為基準。

圖9 螺旋片組數(shù)對Nu和f的影響Fig.9 Influence of group number of spiral pieces on Nu and f

3 結(jié)果與討論

3.1 速度場與二次流

圖10展示了模擬對象的流線圖?？梢钥闯隽黧w在管內(nèi)呈現(xiàn)出不斷分裂合并的流動狀況，與此同時，流體在徑向上還有大量的穿插和交互，這些徑向上的二次流對于強化傳熱具有重要意義。

圖10 模擬對象的流線Fig.10 Streamline of simulation object

為了細致分析流體在徑向上的流動，考察了第七個螺旋片對應(yīng)區(qū)域的橫截面速度場的變化情況，結(jié)果如圖11所示。z為模型橫截面距入口的長度，則z/l可便捷地表示橫截面所處位置，比如z/l為49/8就代表第七個螺旋片軸向1/8長度位置的橫截面，黑色箭頭矢量表示速度在橫截面上的投影；顏色等值線表示軸向速度大小。

圖11 第七個螺旋片不同橫截面的速度場Fig.11 Velocity field of different cross sections of the 7th helical plate

可以發(fā)現(xiàn)每個截面上都存在徑向二次流，這是因為流體在螺旋前進的過程中，軸心附近的流速高于管壁附近的流速，不同流速的流體微團的離心力不同，徑向上產(chǎn)生了逆壓梯度，引起流體從軸心向螺旋外側(cè)管壁的流動，同時還會引起流體沿管壁從螺旋外側(cè)向內(nèi)側(cè)的附加流動，形成一對反向的雙旋流渦，這對旋渦又被稱為Dean渦[16]。圖11中螺旋片兩側(cè)各有一對Dean渦，圖11（a）和圖11（b）中每一對都是兩個同向渦組成，這是因為相鄰的兩個螺旋片頭尾呈90°交叉布置，第六個螺旋片產(chǎn)生的兩對Dean渦進入第七個螺旋片時，四個渦被重組，旋向相同的兩個渦進入了第七個螺旋片的同一側(cè)。由于同向的兩個渦的動能會相互抵消，可以觀察到：從圖11（a）到圖11（b），渦的強度在逐漸變?nèi)?；圖11（c）時，其中一個渦消失；圖11（d）時，渦消失的區(qū)域形成了一個與消失前反向的渦；再從圖11（e）到圖11（h），兩對反向渦不斷發(fā)展、平衡，如此再進入第八個螺旋片。整個流體區(qū)域的二次流都表現(xiàn)出這種周期往復(fù)的特性。

3.2 溫度場及傳熱過程強化

在二次流的影響下，反應(yīng)器的傳熱過程也得到了強化，圖12所示為第七個螺旋片不同橫截面的溫度場情況。顏色等值線表示溫度場，線狀的紋理表示速度矢量在橫截面投影的線積分卷積。

圖12 第七個螺旋片不同橫截面的溫度場Fig.12 Temperature field of different cross sections of the 7th helical plate

由圖12可知，熱邊界層很薄，整個橫截面的絕大部分區(qū)域溫度是均一的，這是由于在二次流的作用下，流動主體區(qū)的冷流體不斷侵擾邊界層的熱流體，邊界層的熱流體不斷被卷入流動主體區(qū)的冷流體，這個過程中不僅流體微團的位置發(fā)生了交錯，流體微團的溫度也在快速交換和平衡，形成了大面積溫度均一的區(qū)域。

為了更加直觀地展現(xiàn)SK螺旋結(jié)構(gòu)強化傳熱的能力，以相同邊界條件下無構(gòu)件空管的溫度場為對照，三個截面的中軸線溫度分布如圖13所示，橫坐標(biāo)χ為橫截面徑向上的無量綱長度?？梢钥闯鲭S著軸向距離的增大，SK結(jié)構(gòu)對應(yīng)的主流區(qū)溫度快速上升并逐漸趨于均一，而空管對應(yīng)的主流區(qū)溫度僅緩慢上升，并且徑向分布很不均勻。在該條件下，SK螺旋結(jié)構(gòu)的傳熱系數(shù)為空管的4.2倍。

圖13 不同橫截面中軸線的溫度分布Fig.13 Temperature distribution of the central axis with different cross sections

3.3 場協(xié)同效應(yīng)分析

將對流傳熱比作有熱源的熱傳導(dǎo)，傳熱速率不僅與流體物性和流速有關(guān)，還與速度場與溫度梯度的夾角有關(guān)。當(dāng)流體被加熱時，應(yīng)使夾角盡量接近0°，當(dāng)流體被冷卻時應(yīng)使夾角接近180°，這一思想被稱為“場協(xié)同原理”[17-18]。

通過式（4）計算不同截面的協(xié)同角，結(jié)果如圖14所示?？梢钥闯鲇休^大面積區(qū)域的協(xié)同角偏離90°，這說明SK螺旋結(jié)構(gòu)產(chǎn)生的二次流對于強化傳熱起著重要的作用。觀察圖14（a）～圖14（c）紅色與藍色區(qū)域的演變，可以發(fā)現(xiàn)綠色區(qū)域占比逐漸增多，這表明流過第六個螺旋片的流體二次流在第七個螺旋片的反旋向?qū)Я飨轮饾u減弱。在第七個螺旋片的后半段，即圖14（d）～圖14（h）所示的截面位置，圖中紅藍色區(qū)域的占比明顯增多，這是由于流體在第七個螺旋片的主導(dǎo)下逐漸形成新的Dean渦，二次流再次發(fā)展增強。

圖14 不同截面的協(xié)同角分布Fig.14 Synergy angle distribution of different cross sections

3.4 雷諾數(shù)與普朗特數(shù)對傳熱的影響

雷諾數(shù)（Re）和普朗特數(shù)（Pr）是影響對流傳熱的重要因素，通過單因素分析的方式考察了兩者對傳熱過程的影響。為了綜合考慮傳熱與流動阻力兩方面的結(jié)果，引入性能評價指標(biāo)（PEC）來量化等泵功約束條件下強化傳熱的程度[19]，計算方法見式（5）。

選取等水力學(xué)直徑的中空圓管作為參照結(jié)構(gòu)，不同Re時的Nu和f值如圖15所示，對PEC的影響如圖16所示?？梢钥闯?，在其他因素保持不變時，Nu隨著Re的增大而增大，f隨著Re的增大而減小，PEC呈現(xiàn)先增大再減小的規(guī)律。在考察范圍內(nèi)，PEC數(shù)值始終大于1，說明SK螺旋結(jié)構(gòu)實現(xiàn)了強化傳熱。

圖15 Re對Nu和f的影響Fig.15 The effect of Re on Nu and f

圖16 Re對PEC的影響Fig.16 The effect of Re on PEC

通過擬合得到：

Pr對Nu和f的影響如圖17所示?？梢钥吹?，Nu隨Pr的增大而增大，f數(shù)值上與Pr無關(guān)。擬合得到Nu與Pr的關(guān)系見式（8），Pr項指數(shù)為0.406，與空管在完全湍流區(qū)的Nu實驗關(guān)聯(lián)式的Pr項指數(shù)0.4基本一致。

圖17 Pr對Nu和f的影響Fig.17 The effect of Pr on Nu and f

3.5 螺旋片長徑比與扭轉(zhuǎn)角對傳熱的影響

幾何結(jié)構(gòu)也是影響傳熱過程的因素之一，對于SK螺旋結(jié)構(gòu)，特征幾何參數(shù)為螺旋片長徑比（Ar）和扭轉(zhuǎn)角（φ）。不同Ar對Nu和f的影響如圖18所示，對PEC的影響如圖19所示。

圖18 Ar對Nu和f的影響Fig.18 The effect of Ar on Nu and f

圖19 Ar對PEC的影響Fig.19 The effect of Ar on PEC

隨著Ar的增大，Nu和f都呈現(xiàn)減小的趨勢，這是由于Ar增大會使得單位軸向長度螺旋片旋轉(zhuǎn)角度的下降，二次流的強度減弱。PEC則隨著Ar的增大緩慢減小，當(dāng)Ar＞2.4時減小速度變快，因此在實際應(yīng)用中Ar設(shè)計值為1.7～2.4較合適。通過擬合得到如下的關(guān)聯(lián)式：

不同φ對Nu和f的影響如圖20所示，對PEC的影響如圖21所示，其中φ0為180°?？梢钥闯鯪u和f隨著φ/φ0的增大而增大，PEC隨φ/φ0的增大而顯著增大，這說明增大扭轉(zhuǎn)角是一種強化SK螺旋結(jié)構(gòu)傳熱的理想方式，在實際應(yīng)用中選擇180°扭轉(zhuǎn)角的螺旋片是合適的。

圖20 φ對Nu和f的影響Fig.20 The effect of φ on Nu and f

圖21 φ對PEC的影響Fig.21 The effect of φ on PEC

擬合得到Nu和f與無因次扭轉(zhuǎn)角的關(guān)系式分別為：

3.6 計算關(guān)聯(lián)式

取3.4和3.5的所有數(shù)據(jù)進行擬合，得到SK螺旋結(jié)構(gòu)的Nu及f的總計算關(guān)聯(lián)式：

式中：866＜Re＜4 330，1＜Pr＜214，1.7＜Ar＜2.9，0.5＜φ/φ0＜1。

圖22和圖23為總關(guān)聯(lián)式計算得到的值與模擬值的比較。

圖22 Nu計算值與擬合值的對比Fig.22 Comparison of calculated Nu and fitted Nu

圖23 f計算值與擬合值的對比Fig.23 Comparison of calculated f and fitted f

由圖22和圖23可知，計算值和擬合值的最大偏差都在10%以內(nèi)。Nu計算關(guān)聯(lián)式的相關(guān)系數(shù)R2為0.998 1，顯著性檢驗中F值為23 089，F(xiàn)0.95=4.07，F(xiàn)10F0.95；f計算關(guān)聯(lián)式的相關(guān)系數(shù)為0.999 7，顯著性檢驗中F值為111 353，F(xiàn)0.95為4.16，F(xiàn)10F0.95。兩個關(guān)聯(lián)式都有著較好的符合性和預(yù)測性。

4 結(jié) 論

a）基于擴大通道尺寸并內(nèi)置三維螺旋內(nèi)構(gòu)件的思路設(shè)計了一種高通量連續(xù)流反應(yīng)器，通量可達1000 mL/min，比傳統(tǒng)連續(xù)流微反應(yīng)器高1～2個數(shù)量級；

b）在周期性螺旋片的作用下，流動徑向上產(chǎn)生大量二次流，這些漩渦顯著了強化傳熱過程；

c）熱邊界層很薄，徑向上溫度近乎均一，SK結(jié)構(gòu)的傳熱系數(shù)為空管的4.2倍；

d）雷諾數(shù)（Re）對性能評價指標(biāo)PEC影響較大，合適的Re為2 000～3 000；

e）螺旋片長徑比（Ar）和扭轉(zhuǎn)角（φ）對PEC影響較大，合適的Ar為1.7～2.4，φ為180°；

f）反應(yīng)器內(nèi)努塞爾數(shù)（Nu）和摩擦系數(shù)（f）計算關(guān)聯(lián)式如下：

式中：866＜Re＜4 330，1＜Pr＜214，1.7＜Ar＜2.9，0.5＜φ/φ0＜1。

符號說明

Ar—— 螺旋片長徑比vi—— 速度張量，m/s

cp—— 流體比熱容，J/(kg·K)z—— 模型橫截面距入口的長度，mm

f—— 流動摩擦系數(shù)θ—— 協(xié)同角，°

l—— 螺片的長度，mmλ—— 流體熱導(dǎo)率，W/(m·K)

Nu—— 努塞爾數(shù)μ—— 流體黏度，Pa·s

Pr—— 普朗特數(shù)ρ—— 流體密度，kg/m3

Re—— 雷諾數(shù)φ—— 扭轉(zhuǎn)角，°

T—— 流體溫度，Kχ—— 橫截面徑向上的無量綱長度