鄔勇奇 溫正慧

(華中科技大學)

流體溫度變化對反應器中RTD測定的影響分析

鄔勇奇**溫正慧

(華中科技大學)

采用脈沖法測量了攪拌反應器中的流體停留時間分布(RTD)，討論了流體溫度變化對停留時間分布的統(tǒng)計特征值及流動模型參數(shù)的影響。當實驗過程中流體溫度升高(或變化)，將導致電導率基線發(fā)生漂移。實驗數(shù)據(jù)處理發(fā)現(xiàn)，這種基線漂移會加大停留時間分布的計算偏差，進一步影響模型參數(shù)的計算正確性和反應器設計。消除基線漂移的影響對獲得正確的計算結果具有積極作用。

反應器 RTD 溫度變化 基線漂移

反應器內進行的過程可分為化學反應過程和傳遞過程。其中，傳遞過程與設備尺寸、結構、流體物性及流動狀況等密切相關。通過冷模試驗所獲得的傳遞過程規(guī)律可用于建立反應器數(shù)學模型。RTD測定就是最重要的實驗內容之一。工業(yè)上，連續(xù)反應器內的流動與混合狀態(tài)通常是各不相同的，反應的結果也受其影響。RTD測定是表征反應器內流動狀態(tài)的有效途徑，對于反應器運行的工況分析、設計與放大等具有重要意義。

有文獻針對設備內停留時間分布的研究進行了大量報道[1～8]。筆者通過RTD相關實驗和數(shù)據(jù)的處理與分析，考察了溫度變化對反應器中流體停留時間分布計算結果(如數(shù)學特征值及模型參數(shù)等)的影響，從而為正確評價反應器的運行狀態(tài)提供參考。

1 實驗裝置及流程

實驗裝置由反應器、水泵、水箱、電導電極、閥門、電導儀及轉子流量計等組成。其中大反應釜內徑150mm、內空高約160mm；小釜內徑110mm、內空高約140mm。主要儀器型號：玻璃轉子流量計，天津市自動化儀表十四廠LZB型，量程0～60L/h；數(shù)字電導率儀，天津市盛邦科技有限公司DDS-302；微型離心泵，陽江市新力工業(yè)有限公司粵華牌不銹鋼微型離心泵，量程1.2～4.8m3/h。此外還有轉速調節(jié)器及精密溫度計等，實驗裝置如圖1所示。

圖1 RTD測定實驗裝置示意圖1——π形管出口；2——攪拌電機； 3、4——示蹤劑加入口； 5——電導儀； 6——三通閥； 7——電導電極； 8——轉子流量計；9——水箱； 10——水泵

實驗中所用流體為水，示蹤劑為飽和KCl溶液。水從水箱由泵經(jīng)流量計計量后送入釜中。KCl示蹤劑用注射器從反應器頂部示蹤劑入口快速加入反應器中。實驗過程中的流體溫度變化由精密溫度計測量，流體電導率的變化由電導率儀測得，轉速由電機調節(jié)。實驗可分為單釜(大釜)與多釜(小釜)串聯(lián)測定。由于水泵提供的量程大于實驗所需流量，故多余流量經(jīng)泵出口旁路回流至水箱中。

2 結果與討論

2.1流體溫度變化對其電導率數(shù)值的影響

在停留時間分布測定實驗過程中，由于部分實驗時間較長(與設備容積及流體流量有關)，或其他原因導致流體溫度發(fā)生變化。這對RTD的測定和數(shù)據(jù)的處理會產生一定程度的影響。要消除溫度變化對計算結果的影響，首先要掌握流體電導率與溫度變化的關系曲線?？紤]到實驗環(huán)境是在12～18℃溫度范圍，本實驗對14～17℃的水流體(載體)電導率隨溫度變化的情況進行了測定，結果見表1(不同的水質測量數(shù)據(jù)會有所不同)。

表1 水流體電導率隨溫度變化

從表1可以看出，流體溫度每變化1℃，其電導率將產生(6.7～7.2)×10-3ms/cm的變化。本實驗是通過測量示蹤劑KCl溶液的電導率來獲得反應器的停留時間分布的。載流體水本身存在一個電導率值，也即實驗開始時流體電導率的初值。實驗過程中由于泵體做功等原因使得流體溫度升高，當采用脈沖法測量時，實驗結束時的電導率終值會一直大于實驗初值，導致電導率無法回歸原點。這就造成了實驗基線漂移(抬升)，并影響最終計算結果。

2.2實驗數(shù)據(jù)處理方法(原則)

由于實驗中使用的水流體并非去離子水，本身有一定的電導率，且隨實驗中溫度的變化而變化。設實驗中含示蹤劑水流體的電導率儀表讀數(shù)為μ1(t)，水流體本身的電導率值為μ0(t)，則示蹤劑本身的電導率值μ(t)為：

μ(t)=μ1(t)-μ0(t)

(1)

式中t——流體在反應器中的停留時間，s。

實驗中流體溫度不斷升高，水流體本身的電導率μ0(t)的計算式為：

μ0(t)=μ0(0)+k1(T-T0)

(2)

式中k1——溫度變化導致水流體電導率變化的斜率，k1=(6.7～7.2)×10-3，可根據(jù)實驗確立(本實驗取范圍內的均值)；

T——實驗測量的流體溫度；

T0——實驗流體初始溫度；

μ0(0)——水流體初始電導率，即實驗初值(原點)。

溫度未校正(不考慮基線漂移)的情況下，示蹤劑電導率μ(t)的計算公式為：

μ(t)=μ1(t)-μ0(0)

(3)

溫度已校正(已考慮基線漂移)的情況下，示蹤劑電導率μ(t)的計算公式為：

μ(t)=μ1(t)-μ0(t)=μ1(t)-[μ0(0)+k1(T-T0)]

(4)

不同時間t下的溫度T可通過精密溫度計實時測量。假設示蹤劑濃度為c(t)，實驗范圍內示蹤劑濃度近似正比于其電導率μ(t)：

c(t)=k·μ(t)

(5)

其中，k為比例系數(shù)。停留時間分布密度函數(shù)E(t)計算式為：

(6)

流體的平均停留時間和方差計算式為：

(7)

(8)

從以上公式可以看出，示蹤劑的電導率-時間曲線相當于示蹤劑濃度-時間曲線，該曲線可以比較全面地反映容器內流體的RTD分布情況。

2.3流體溫度變化未校正的RTD計算及結果

(9)

實驗條件下反應器的空間時間τ可按反應器有效體積V與進口體積流量Q之比計算:

τ=V/Q

(10)

幾組條件下的實驗數(shù)據(jù)計算結果見表2。

表2 不同條件下RTD數(shù)據(jù)計算結果(溫度未校正)

注：編號中的S表示單釜實驗，下同。

2.4流體溫度變化校正后的RTD計算及結果

按公式(4)修正溫度變化引起的基線漂移影響，再按式(5)～(10)對實驗數(shù)據(jù)進行處理后，所得計算結果列于表3。

表3 不同條件下RTD數(shù)據(jù)計算結果(溫度已校正)

2.5對比分析

比較分析實驗條件下的數(shù)據(jù)發(fā)現(xiàn)，溫度修正與否對RTD計算中的相關曲線影響很大。圖2是計算平均停留時間時t×E(t)-t曲線對比(實驗序號3S)。

圖2 t×E(t)-t曲線對比圖

在進行RTD方差的計算時，這種曲線的差異表現(xiàn)得更加明顯。在轉速80r/min、流量16L/h實驗條件下，所得的t2×E(t)-t曲線如圖3所示。

圖3 t2×E(t)-t曲線對比圖

圖3中兩條曲線的差異在實驗后半期十分顯著。其原因在于未進行基線漂移修正的電導率數(shù)據(jù)無法回歸原點，E(t)殘存值大，致使曲線發(fā)生翹尾現(xiàn)象(時間越長越嚴重)。這樣，模型參數(shù)的計算結果也將失真。

實際上，影響RTD實驗數(shù)據(jù)測定的因素較多。除流體溫度外，還與攪拌轉速(高速、中速、低速)、流體流量、反應器尺寸與結構以及實驗終止點的正確判定等一系列因素有關，可進一步深入探討。

3 結束語

在進行停留時間分布實驗測定時，如果流體的溫度在實驗過程中發(fā)生變化，這將導致基于電導率測量的實驗基線發(fā)生相應的漂移。在此基礎上的停留時間分布函數(shù)的計算(包括平均停留時間與方差等)將產生一定程度的偏差，從而進一步影響模型參數(shù)的計算與反應器設計。對此電導率基線漂移進行修正，實驗條件下所得平均停留時間計算結果與空間時間吻合，模型參數(shù)N值約為1.1～1.2，比較接近全混流模型。

[1] 李紹芬. 反應工程[M].北京：化學工業(yè)出版社,2013：135～150.

[2] Choi B S, Wan B, Philyaw S, et al. Residence Time Distributions in a Stirred Tank: Comparison of CFD Predictions with Experiment[J]. Industrial & Engineering Chemistry Research, 2004, 43(20): 6548～6556.

[3] 紀利俊,劉海峰,王輔臣，等. 撞擊流反應器停留時間分布[J].華東理工大學學報:自然科學版, 2006, 32(1): 24～27.

[4] Levenspiel O. Chemical Reaction Engineering[M]. New York: John Wiley & Sons,1999:260～338.

[5] 王勤輝,駱仲泱,倪明江，等. 循環(huán)流化床內顆粒停留時間分布[J]. 燃燒科學與技術, 2001, 7(4): 221～225.

[6] 陳國南,李廣贊,王嘉俊，等. 泰勒反應器中流體流動及停留時間分布研究[J].化學工程, 2005,33(6):22～25.

[7] Raghuraman J, Varma Y B G. A Model for Residence Time Distribution in Multistage Systems with Cross-flow Between Active and Dead Regions [J]. Chemical Engineering Science, 1973, 28(2): 585～591.

[8] Patwardhan,Ashwin W. Prediction of Residence Time Distribution of Stirred Reactors[J].Industrial & Engineering Chemistry Research, 2001, 40(24): 5686～5695.

AnalysisofFluidTemperatureInfluenceonRTDinAgitationReactor

WU Yong-qi, WEN Zheng-hui

(HuazhongUniversityofScienceandTechnology,Wuhan430074,China)

The pulse method was adopted to measure fluid’s residence time distribution (RTD) in agitation reactor. The temperature change’s influences on RTD’s statistical characteristic value and the fluid model parameters were discussed to show that, the temperature’s rise or variation would result in conductivity baseline drift which can increase the error in calculating RTD and it further influences both model parameters and reactor design. Removing the baseline drift’s influence on the calculation results has positive role.

reactor, RTD, temperature variation, baseline drift

**鄔勇奇，男，1964年1月生，副教授。湖北省武漢市，430074。

TQ052

A

0254-6094(2015)02-0184-04

2014-07-30，

2014-08-21)