陳　杭，孫　澤，宋興福，于建國

（華東理工大學(xué)國家鹽湖資源綜合利用工程技術(shù)研究中心，上海 200237）

?

基于均齡理論高效計算萃取塔軸向混合分布

陳杭，孫澤，宋興福，于建國

（華東理工大學(xué)國家鹽湖資源綜合利用工程技術(shù)研究中心，上海 200237）

在驗證了CFD單相流場模擬的基礎(chǔ)上，采用均齡理論計算了中試轉(zhuǎn)盤塔內(nèi)的軸向混合分布，并將計算結(jié)果和理論平均停留時間以及組分輸運模型計算值進行對比。結(jié)果表明：均齡理論能準確預(yù)測轉(zhuǎn)盤塔內(nèi)的軸向混合信息，且其計算時間只需數(shù)十秒，遠小于傳統(tǒng)組分輸運模型所需的兩周時間，具有低計算量的特點；同時均齡理論克服了傳統(tǒng)組分輸運模型無法模擬軸向混合空間分布的缺陷，為萃取塔內(nèi)部結(jié)構(gòu)優(yōu)化提供了更多信息，是一種高效的模擬方法。后續(xù)均齡理論模擬結(jié)果的分析預(yù)示著轉(zhuǎn)盤塔內(nèi)的流動近似呈現(xiàn)出級內(nèi)全混、級間平推的特點，符合萃取操作的需求；而相對于轉(zhuǎn)盤間良好的混合作用，靜環(huán)間存在明顯的流動死區(qū)，造成一定的非理想性，其結(jié)構(gòu)有待于進一步的優(yōu)化。

萃取；轉(zhuǎn)盤塔；計算流體力學(xué)；均齡理論；數(shù)值模擬；軸向混合；

DOI：10.11949/j.issn.0438-1157.20151476

引　言

轉(zhuǎn)盤塔是一種重要的萃取設(shè)備，因其結(jié)構(gòu)簡單，處理量大，安裝維修方便的優(yōu)點而廣泛應(yīng)用于化工分離過程［1］。相比較于其他攪拌型萃取塔，轉(zhuǎn)盤塔的能耗較低，但其萃取效率也明顯偏低［2］，這與塔內(nèi)軸向混合有著密切的聯(lián)系［3］。萃取塔內(nèi)的攪拌混合在保證兩相充分接觸的同時也會引起返混、前混、溝流等各種非理想性流動，進而降低傳質(zhì)推動力及萃取效率［4］。大型工業(yè)塔中，通常有60%以上的塔高專門用于補償軸向返混所引起的傳質(zhì)效率的降低［5］。因此，轉(zhuǎn)盤塔內(nèi)軸向混合的研究對于萃取塔的設(shè)計和強化具有重要意義［6］。

目前，關(guān)于軸向混合的研究主要基于軸向混合模型，即采用單一的軸向混合系數(shù)來描述流動的非理想性［7］。一般在獲得停留時間分布后，可利用軸向混合模型近似解計算得到軸向混合系數(shù)，進而評估塔內(nèi)軸向返混程度。該方法已成功應(yīng)用于轉(zhuǎn)盤塔內(nèi)軸向混合系數(shù)的計算［8-9］，且其中停留時間分布曲線的獲得也從傳統(tǒng)的示蹤法實驗測量［10-12］發(fā)展到更為高效的計算流體力學(xué)（CFD）模擬［13-16］。但依據(jù)停留時間分布計算得到的軸向混合系數(shù)只體現(xiàn)了進料點與監(jiān)測點之間軸向混合的整體信息，無法反映塔內(nèi)軸向混合的具體分布情況，因而，其結(jié)果可用于整塔設(shè)計及校核計算［17］，但對塔內(nèi)內(nèi)構(gòu)件的開發(fā)和改進并無針對性的指導(dǎo)意義。另外，相比于實驗測量，CFD模擬確實提高了軸向混合的評估效率，但停留時間分布的模擬意味著長時間的非穩(wěn)態(tài)計算，故其應(yīng)用仍受到高額計算成本的限制，尤其是對于大規(guī)模塔器而言。

均齡理論［18-19］的提出為軸向混合的計算提供了一種新的途徑。其基本思想是從軸向混合模型出發(fā)，推導(dǎo)得出停留時間分布各階矩量的控制方程，并利用CFD求解進而直接計算得到軸向混合結(jié)果。由于各階矩量是停留時間分布隨時間積分處理后的結(jié)果，其應(yīng)用將控制方程內(nèi)原有的非穩(wěn)態(tài)項轉(zhuǎn)化為穩(wěn)態(tài)源項，因此，均齡理論避免了停留時間分布的直接模擬，在提高了計算效率的同時獲得了軸向混合的分布信息。簡而言之，均齡理論通過犧牲意義較小的時間信息得到了更多、更有用的空間分布信息，從而克服了傳統(tǒng)組分輸運模型的缺點。

目前，均齡理論已成功應(yīng)用于攪拌槽反應(yīng)器［20-22］及靜態(tài)混合器［23］的混合研究，但未有應(yīng)用于萃取塔的相關(guān)報道，且有關(guān)萃取塔內(nèi)軸向混合的研究大都基于傳統(tǒng)的停留時間分布計算方法，缺乏塔內(nèi)軸向混合空間分布的探索，另外，不同于傳統(tǒng)的二維或三維模型，轉(zhuǎn)盤塔的模擬往往采用二維旋轉(zhuǎn)軸對稱假設(shè)［24］，在此簡化條件下均齡理論是否仍適用亦未可知，為此，本文利用均齡理論研究了中試轉(zhuǎn)盤塔內(nèi)的軸向混合。

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

1.1幾何及物理模型設(shè)置

圖1　中試轉(zhuǎn)盤塔幾何模型及網(wǎng)格方案Fig.1　Geometry and grid scheme of pilot-scale RDC（a） column geometry； （b） geometrical sizes； （c） grid layout

中試轉(zhuǎn)盤塔塔徑為150 mm，轉(zhuǎn)軸直徑54 mm，模型采用二維旋轉(zhuǎn)軸對稱假設(shè)，其具體幾何及內(nèi)部尺寸如圖1所示。計算域采用四邊形網(wǎng)格進行離散，根據(jù)前期計算結(jié)果［15］，網(wǎng)格無關(guān)性尺寸選定為0.5 mm。模擬的流體介質(zhì)為水，其塔頂入口和塔底出口均采用速度邊界條件，除特殊說明外，其流量為200 L·h-1，轉(zhuǎn)盤轉(zhuǎn)速為150 r·min-1，所有壁面采用無滑移邊界條件。數(shù)值離散采用QUICK格式，模型計算收斂標準為10-4（組分輸運模型為10-7），時間步長為0.001 s，以保證Courant數(shù)小于1。

1.2CFD單相流模型

CFD流場模擬是組分輸運模型及均齡理論計算的基礎(chǔ)，為此建立了轉(zhuǎn)盤塔的CFD單相流模型，其控制方程主要包括連續(xù)性方程和動量守恒方程

式中，ρ為流體密度，t為時間，u為流體速度，p為壓力，τ為切應(yīng)力，g為重力加速度。

Drumm等［25］對不同湍流模型的模擬結(jié)果進行了比較，并推薦采用雷諾應(yīng)力模型（RSM）計算轉(zhuǎn)盤塔單相流場，故本文直接采用RSM封閉時均化后的動量方程，其具體模型可見文獻［26］。

1.3組分輸運模型

在CFD流場的基礎(chǔ)上進一步建立組分輸運模型，作為均齡理論的對比驗證。其控制方程為

式中，Y為組分濃度，D為擴散系數(shù)，Dm為分子擴散系數(shù)，μt為湍流黏度，Sct為Schmidt數(shù)。

在CFD流場穩(wěn)定后，采用0.5 s時間間隔內(nèi)純組分進樣的方式模擬脈沖信號，同時，在圖1中點P1、P2、P3、P4、P5、P6及出口處開始監(jiān)測組分濃度，從而得到各監(jiān)測位置離散化的停留時間分布曲線。萃取塔軸向混合以量綱1 Peclet數(shù)（Pe）表示，其計算［27］如下

式中，ta表示平均停留時間，σ表示停留時間分布的方差，u表示特征速度，L表示特征長度。

1.4均齡理論

Liu等［18］以一階矩量為例給出了均齡理論的推導(dǎo)過程，本文則列出了其通用形式。

首先，式（3）兩邊同乘以tn并積分可得

根據(jù)分部積分，式（8）左邊第1項有

在無限長時間域上，塔內(nèi)濃度終將為0，故有

式（10）代入式（8）則有

停留時間分布的n階矩量（Mn）定義為

則根據(jù)式（11）可得均齡理論的控制方程為

由矩量定義可知，停留時間分布零階矩為量綱1單元1，其一階矩為平均停留時間，求解得到一階和二階矩量值，即可利用式（6）和式（7）計算Peclet數(shù)?？刂品匠淌剑?3）在壁面和出口采用零通量邊界，而在入口處，組分剛進入計算域，因此有

所有模型均采用商業(yè)軟件Fluent 14.0求解，均齡理論則通過其中的自定義標量方程（UDS）和自定義函數(shù)（UDF）功能實現(xiàn)。

2　實驗結(jié)果與討論

2.1轉(zhuǎn)盤塔單相流場模擬結(jié)果

無論是組分輸運模型還是均齡理論，都是基于CFD穩(wěn)定流場進行計算的，因此，有必要對轉(zhuǎn)盤塔單相流場模擬結(jié)果進行驗證和討論。

Drumm等［25］利用粒子圖像測速儀（PIV）測量了中試轉(zhuǎn)盤塔內(nèi)沿直線L1（圖1）的速度分布，本文模擬結(jié)果與其報道實驗數(shù)據(jù)之間的對比如圖2所示。結(jié)果顯示模擬和實驗值有著較好的吻合，說明模型能合理預(yù)測轉(zhuǎn)盤塔的單相流場，為后續(xù)組分輸運模型和均齡理論計算提供了基礎(chǔ)保障。

圖3為轉(zhuǎn)盤塔內(nèi)單相流場的速度云圖及矢量圖，其中圖3（a）表示塔內(nèi)流體總速度，圖3（b）則表示徑向和切向合成的平面速度。從圖中可以看出，轉(zhuǎn)盤塔內(nèi)流體速度從轉(zhuǎn)軸至塔壁逐漸減小，轉(zhuǎn)盤之間的流速明顯大于靜環(huán)之間的速度，這主要歸因于轉(zhuǎn)盤轉(zhuǎn)動引起的顯著切向運動。其次，圖3（b）顯示轉(zhuǎn)盤及靜環(huán)間存在明顯的渦流結(jié)構(gòu)，這是攪拌所引起的離心力和重力共同作用的結(jié)果，也是圖2中速度分布近似呈現(xiàn)出M形狀的主要原因。由于渦流結(jié)構(gòu)和塔體幾何結(jié)構(gòu)存在一定差異，可以看出靜環(huán)角落處存在明顯的流動死區(qū)。

圖2　轉(zhuǎn)盤塔單相流場速度分布驗證Fig.2　Validation of CFD simulated velocity distribution for RDC

2.2均齡理論的模型驗證

為充分說明均齡理論模型的準確性，本文采用兩種方式進行了驗證，一是與理論平均停留時間計算值的對比，二是與傳統(tǒng)組分輸運模型計算結(jié)果的對比。

對于轉(zhuǎn)盤塔而言，在幾何與流量恒定的條件下，組分從進入塔體到離開塔體的理論平均停留時間（ta）可由式（15）計算

式中，V為轉(zhuǎn)盤塔有效體積，Q為轉(zhuǎn)盤塔操作流量，且該式已得到理論和實驗雙重驗證［28-29］。

通過均齡理論計算了中試轉(zhuǎn)盤塔在不同流量條件下出口處的平均停留時間，并與式（15）的理論計算值進行對比，其結(jié)果如圖4所示?？梢园l(fā)現(xiàn)計算值和理論值有著非常好的吻合，模型計算誤差均小于0.1%，有效說明了均齡理論的合理性和準確性。

圖4　均齡理論計算值與理論平均停留時間的對比Fig.4　Comparison of mean residence time data calculated by mean age theory with theoretical values

圖5　不同位置模擬得到的停留時間分布曲線Fig.5　Simulated residence time distribution curves at different positions

與此同時，均齡理論計算結(jié)果也與傳統(tǒng)組分輸運模型進行了對比。采用組分輸運模型，得到了點P1、P2、P3、P4、P5、P6（圖1）及出口處的停留時間分布曲線，其模擬結(jié)果如圖5所示，根據(jù)分布曲線所得的各位置處平均停留時間計算結(jié)果則如表1所示。通過對比發(fā)現(xiàn)，均齡理論模擬值和組分輸運模型計算值有著很好的吻合，兩者之間相差不超過0.5%，這進一步說明了均齡理論的準確性，并確保模型能用于轉(zhuǎn)盤塔內(nèi)軸向混合的考察。

表1　R均齡理論和組分輸運模型計算結(jié)果對比Table 1　Comparison of mean age theory and species transport model

根據(jù)上述結(jié)果，均齡理論在模擬精度上和傳統(tǒng)組分輸運模型大致相當，但其計算量得到了大幅度的下降。在惠普工作站（Xeon W3520，16GB內(nèi)存）上模擬圖5所示的停留時間分布曲線需要兩周時間，而利用均齡理論，在同樣計算條件下耗時不到1 min，其主要原因在于均齡理論將傳統(tǒng)組分輸運模型的非穩(wěn)態(tài)模擬轉(zhuǎn)化為了穩(wěn)態(tài)計算。因此，均齡理論在獲得軸向混合空間分布信息的同時大大提高了計算效率，是一種高效的計算方法。

2.3轉(zhuǎn)盤塔內(nèi)的軸向混合

均齡理論計算結(jié)果提供了轉(zhuǎn)盤塔內(nèi)有關(guān)于軸向混合空間分布的詳細信息。圖6顯示了中試轉(zhuǎn)盤塔內(nèi)平均停留時間的分布情況，為更清晰地體現(xiàn)塔內(nèi)平均停留時間分布狀況，采用了三維作圖方案。從圖中分布可以看出，在塔節(jié)間，平均停留時間從塔頂?shù)剿字饾u增加，而在塔節(jié)內(nèi)，平均停留時間則大致相當。近似來說，整個萃取塔呈現(xiàn)出塔節(jié)內(nèi)全混流、塔節(jié)間平推流的特點，這對于萃取操作而言是有利的，塔節(jié)內(nèi)的混合使得兩相充分接觸傳質(zhì)，而塔節(jié)間的平推又抑制了軸向混合所帶來的推動力的降低。

相比于轉(zhuǎn)盤之間的結(jié)果，靜環(huán)之間平均停留時間的分布存在明顯的不均勻，角落和中心處的平均停留時間值明顯偏高，這主要歸因于流動死區(qū)的存在。從2.1節(jié)關(guān)于流場的討論可知，靜環(huán)角落處的流動死區(qū)主要歸因于流場結(jié)構(gòu)和塔節(jié)幾何不一致所導(dǎo)致的低流速；而靜環(huán)隔室中心的流動死區(qū)則由循環(huán)渦流結(jié)構(gòu)引起，當流體被卷入靜環(huán)間的渦流時，循環(huán)作用會導(dǎo)致其擁有更大的停留時間，從這個角度講，流動死區(qū)不一定都是低流速引起的，而相比于流場分析，均齡理論在流動死區(qū)識別方面更具優(yōu)勢［18］。

圖6　轉(zhuǎn)盤塔內(nèi)平均停留時間的分布Fig.6　Distribution of mean residence time in RDC

圖7　轉(zhuǎn)盤塔內(nèi)軸向混合系數(shù)及Peclet數(shù)分布Fig.7　Distributions of axial mixing coefficient and Peclet number in RDC

圖7為轉(zhuǎn)盤塔內(nèi)軸向混合系數(shù)和Peclet數(shù)的分布云圖。從圖7（a）可以看到，轉(zhuǎn)盤間的軸向混合系數(shù)明顯大于靜環(huán)間的結(jié)果，這預(yù)示著轉(zhuǎn)盤間的隔室具有更為明顯的攪拌作用，主要用于流體混合，而靜環(huán)間的隔室混合較弱，主要起分相作用，這也是轉(zhuǎn)盤之間平均停留時間分布更為均勻的原因。根據(jù)均齡理論一階矩M1和二階矩M2計算得到的塔內(nèi)Peclet數(shù)分布如圖7（b）所示，從圖中可以看到，靜環(huán)間的Peclet數(shù)要大于對應(yīng)轉(zhuǎn)盤間的結(jié)果，且其分布更不均勻。由定義可知，Peclet數(shù)是對流和擴散作用相對強弱的綜合描述，靜環(huán)間Peclet數(shù)較大主要是因為靜環(huán)間軸向混合系數(shù)小，湍流擴散作用弱，另一方面則是由于靜環(huán)間渦流結(jié)構(gòu)所導(dǎo)致的強對流。實際上，靜環(huán)間渦流結(jié)構(gòu)的本質(zhì)可認為是重力和離心力雙重作用導(dǎo)致的流量分配不均勻性，進而形成了強對流與流動死區(qū)并存的流場結(jié)構(gòu)，由于靜環(huán)間擴散作用相對較弱，故這種流量分配的不均勻性最終亦體現(xiàn)在了Peclet數(shù)的分布上。因而， Peclet數(shù)的分布進一步說明萃取塔內(nèi)轉(zhuǎn)盤間的混合作用較好，而靜環(huán)間流動存在明顯的非理想性，有待于結(jié)構(gòu)上的優(yōu)化。

3　結(jié)　論

本文將均齡理論應(yīng)用于轉(zhuǎn)盤塔內(nèi)軸向混合的研究，在CFD模型驗證的基礎(chǔ)上對塔內(nèi)軸向混合的計算結(jié)果進行了討論，主要結(jié)論如下。

（1）流場與均齡理論計算的驗證結(jié)果表明基于CFD模擬的均齡理論能有效應(yīng)用于萃取塔內(nèi)軸向混合的評估，且均齡理論通過矩量定義巧妙地將非穩(wěn)態(tài)項轉(zhuǎn)化為穩(wěn)態(tài)源項，其計算時間不到1 min，遠遠小于傳統(tǒng)組分輸運模型的兩周，具有低計算成本的優(yōu)勢。

（2）傳統(tǒng)組分輸運模型對于塔內(nèi)軸向混合的計算主要是基于停留時間分布的模擬，這意味著其計算結(jié)果只能反映塔內(nèi)軸向混合的平均信息。相比之下，均齡理論能提供塔內(nèi)軸向混合的空間分布信息，為萃取塔內(nèi)部結(jié)構(gòu)的設(shè)計優(yōu)化提供更多依據(jù)，具有信息量大的特點。

（3）均齡理論計算結(jié)果顯示轉(zhuǎn)盤塔內(nèi)的流動近似呈現(xiàn)出級內(nèi)全混、級間平推的特點，符合萃取塔的操作需求。但相比于轉(zhuǎn)盤間充分的混合作用，靜環(huán)間的隔室存在兩種類型的流動死區(qū)，造成了明顯的非理想性，其結(jié)構(gòu)有待于進一步的優(yōu)化。

上述結(jié)果說明均齡理論能高效應(yīng)用于轉(zhuǎn)盤塔的軸向混合研究，下一步工作應(yīng)在現(xiàn)有基礎(chǔ)上將模型進一步推廣至多相流，以滿足實際工業(yè)的設(shè)計優(yōu)化需求。

符號說明

D，Dm——分別為總擴散系數(shù)及分子擴散系數(shù)，m2·s-1

g ——重力加速度，m·s-2

L ——特征長度，m

Mn——停留時間分布n階矩，sn

Pe——Peclet數(shù)

p——壓力，Pa

Q——流量，m3·s-1

Sct——湍流Schmidt數(shù)

t——時間，s

ta——平均停留時間，s

u——流體速度，m·s-1

u——特征速度，m·s-1

V——有效塔體積，m3

Y——組分濃度（質(zhì)量分數(shù)）

μt——湍流黏度，kg·m-1·s-1

ρ——密度，kg·m-3

σ——停留時間分布方差

τ——切應(yīng)力，kg·m-1·s-2

References

［1］ MOHANTY S. Modeling of liquid-liquid extraction column： a review［J］. Rev. Chem.， 2000， 16（3）： 199-248. DOI： 10.1515/REVCE. 2000.16.3.199.

［2］ TORAB-MOSTAEDI M， SAFDARI J， MOOSAVIAN M A， et al. Stage efficiency of Hanson mixer-settler extraction column ［J］. Chem. Eng. Process.， 2009， 48（1）： 224-228. DOI： 10.1016/j.cep. 2008.03.010.

［3］ 王運東， 費維揚， 劉小秦， 等. 新型轉(zhuǎn)盤萃取塔研究開發(fā)與工業(yè)應(yīng)用 ［J］. 化學(xué)工程， 2008， 36（4）： 1-5. DOI： 10.3969/j.issn.1005-9954. 2008.04.001. WANG Y D， FEI W Y， LIU X Q， et al. Research， development and industrial application of new rotating disc contactor ［J］. Chem. Eng.（China）， 2008， 36（4）： 1-5. DOI： 10.3969/j.issn.1005-9954. 2008.04.001.

［4］ 費維揚. 萃取塔設(shè)備研究和應(yīng)用的若干新進展 ［J］. 化工學(xué)報，2013， 64（1）： 44-51. DOI： 10.3969/j.issn.0438-1157.2013.01.007. FEI W Y. Progresses of study and application on extraction columns［J］. CIESC Journal， 2013， 64（1）： 44-51. DOI： 10.3969/j.issn. 0438-1157.2013.01.007.

［5］ LI H B， TANG X J， LUO G S， et al. Development of two-point dynamic method for evaluating extraction columns ［J］. Can. J. Chem. Eng.， 2004， 82（3）： 471-477. DOI： 10.1002/cjce.5450820307.

［6］ FEI W Y， WANG Y D， WAN Y K. Physical modelling and numerical simulation of velocity fields in rotating disc contactor via CFD simulation and LDV measurement ［J］. Chem. Eng. J.， 2000， 78（2）：131-139. DOI： 10.1016/S1385-8947（00）00156-X.

［7］ BART H J. Reactive extraction in stirred columns—a review ［J］. Chem. Eng. Technol.， 2003， 26（7）： 723-731. DOI： 10.1002/ceat. 200306102.

［8］ ZHANG S H， NI X D， SU Y F. Hydrodynamics， axial mixing and mass transfer in rotating disk contactors ［J］. Can. J. Chem. Eng.， 1981，59（5）： 573-583. DOI： 10.1002/cjce.5450590503.

［9］ ZHU J W， ZHANG S H， ZHOU X K， et al. Hydrodynamics， axial mixing and mass transfer in open turbine rotating disc contactor（OTRDC） ［J］. Chem. Eng. Technol.， 1991， 14（3）： 167-177. DOI：10.1002/ceat.270140304.

［10］ STEINER L， KUMAR A， HARTLAND S. Determination andcorrelation of axial-mixing parameters in an agitated liquid-liquid extraction column ［J］. Can. J. Chem. Eng.， 1988， 66（2）： 241-247. DOI：10.1002/cjce.5450660208.

［11］ MORíS M A， DíEZ F V， COCA J. Hydrodynamics of a rotating disc contactor ［J］. Sep. Purif. Technol.， 1997， 11（2）： 79-92. DOI：10.1016/S1383-5866（96）01003-9.

［12］ MURUGESAN T， REGUPATHI I. A correlation for dispersed phase axial mixing coefficients in rotating disc contactors ［J］. J. Chem. Technol. Biotechnol.， 2006， 81（4）： 721-726. DOI： 10.1002/jctb.1389.

［13］ DESHMUKH S S， SATHE M J， JOSHI J B. Residence time distribution and flow patterns in the single-phase annular region of annular centrifugal extractor ［J］. Ind. Eng. Chem. Res.， 2009， 48（1）：37-46. DOI： 10.1021/ie800231d.

［14］ CHARTON S， DUHAMET J， BORDA G， et al. Axial dispersion in pulsed disk and doughnut columns： a unified law ［J］. Chem. Eng. Sci.，2012， 75： 468-477. DOI： 10.1016/j.ces.2012.04.011.

［15］ CHEN H， SUN Z， SONG X F， et al. Key parameter prediction and validation for a pilot-scale rotating disk contactor by CFD-PBM simulation ［J］. Ind. Eng. Chem. Res.， 2014， 53（51）： 20013-20023. DOI： 10.1021/ie503115w.

［16］ 毛在砂， 楊超. 化學(xué)反應(yīng)器宏觀混合研究展望 ［J］. 化工學(xué)報，2015， 66（8）： 2795-2804. DOI： 10.11949/j.issn.0438-1157.20150628. MAO Z S， YANG C. Perspective to study on macro-mixing in chemical reactors ［J］. CIESC Journal， 2015， 66（8）： 2795-2804. DOI：10.11949/j.issn.0438-1157.20150628.

［17］ NETO A P， MANSUR M B. Transient modeling of zinc extraction with D2EHPA in a Kühni column ［J］. Chem. Eng. Res. Des.， 2013，91（12）： 2323-2332. DOI： 10.1016/j.cherd.2013.05.014.

［18］ LIU M， TILTON J N. Spatial distributions of mean age and higher moments in steady continuous flows ［J］. AIChE Journal， 2010，56（10）： 2561-2572. DOI： 10.1002/aic.12151.

［19］ LIU M. A method for computing the degree of mixing in steady continuous flow systems ［J］. Chem. Eng. Sci.， 2011， 66（13）：3045-3048. DOI：10.1016/j.ces.2011.03.049.

［20］ LIU M. Quantitative characterisation of mixing in stirred tank reactors with mean age distribution ［J］. Can. J. Chem. Eng.， 2011，89（5）： 1018-1028. DOI： 10.1002/cjce.20563.

［21］ LIU M. Prediction of tracer concentration and mixing in CFSTRs with mean age distribution ［J］. Ind. Eng. Chem. Res.， 2011， 50（9）：5838-5851. DOI： 10.1021/ie2002395.

［22］ LIU M. Age distribution and the degree of mixing in continuous flow stirred tank reactors ［J］. Chem. Eng. Sci.， 2012， 69（1）： 382-393. DOI：10.1016/j.ces.2011.10.062.

［23］ LIU M. Age distribution in the Kenics static micromixer with convection and diffusion ［J］. Ind. Eng. Chem. Res.， 2012， 51（20）：7081-7094. DOI： 10.1021/ie200716v.

［24］ CHEN H， SUN Z， SONG X F， et al. A pseudo-3D model with 3D accuracy and 2D cost for the CFD-PBM simulation of a pilot-scale rotating disc contactor ［J］. Chem. Eng. Sci.， 2016， 139： 27-40. DOI：10.1016/j.ces.2015.09.023.

［25］ DRUMM C， BART H J. Hydrodynamics in a RDC extractor： single and two-phase PIV measurements and CFD simulations ［J］. Chem. Eng. Technol.， 2006， 29（11）： 1297-1302. DOI： 10.1002/ceat. 200600212.

［26］ GIBSON M M. Ground effects on pressure fluctuations in the atmospheric boundary layer ［J］. J. Fluid Mech.， 1978， 86（3）： 491-511. DOI： 10.1017/S0022112078001251.

［27］ LEVENSPIEL O. Chemical Reaction Engineering［M］. 3rd ed. New York： John Wiley & Sons， Inc.， 1999： 293-310.

［28］ GAO H， FU Y. Mean residence time of Markov processes for particle transport in fluidized bed reactors ［J］. J. Math. Chem.， 2011， 49（2）：444-456. DOI： 10.1007/s10910-010-9751-x.

［29］ GAO J， WALSH G C， BIGIO D， et al. Mean residence time analysis for twin screw extruders ［J］. Polym. Eng. Sci.， 2000， 40（1）： 227-237. DOI： 10.1002/pen.11155.

Efficient estimation of axial mixing distribution for extraction in rotating disc contactor based on simulation using mean age theory

CHEN Hang， SUN Ze， SONG Xingfu， YU Jianguo
（National Engineering Research Center for Integrated Utilization of Salt Lake Resource， East China University of Science and Technology， Shanghai 200237， China）

In this study， on basis of verification of single-phase flow field simulation by computational fluid dynamics （CFD）， the mean age theory was employed to estimate axial mixing distribution in a pilot-scale rotating disc contactor （RDC）. In order to validate the mean age theory， two results of theoretical mean residence time and the simulated values by species transport model were compared. It showed that the mean age theory is capable to relatively accurately predict the axial mixing distribution in RDC within ten seconds of computional time， far shorter than two weeks by traditional species transport model， hence possesses the advantage of low computational cost. Moreover， the mean age theory is also capable to provide the spatial distribution information of axial mixing distribution， while the traditional species transport model does not. Hence， the mean age theory can provide more information for structural optimization， and is an efficient method to calculate the axial mixing distribution for extraction. Subsequent analysis on hydrodynamics in RDC showed that it is appropriately mixed flow in compartments and plug flow between compartments， which is favorable for extraction operation. Compared with the well mixing between rotors， the flow field between stators has apparent dead flow zones，which are mainly attributed to two factors of low flow velocity and cyclic vortex structure. The dead zones cause acertain flow non-idealities and are disadvantage for extraction efficiency. Hence， the structure of the stators in RDC needs to be further improved.

date： 2015-09-21.

Prof. SUN Ze， zsun@ecust.edu.cn； Prof. YU Jianguo， jgyu@ecust.edu.cn

supported by the National Natural Science Foundation of China （21206038）.

extraction；rotating disc contactor；CFD；mean age theory；numerical simulation；axial mixing

TQ 021.1

A

0438—1157（2016）05—1694—07

2015-09-21收到初稿，2016-01-26收到修改稿。

聯(lián)系人：孫澤，于建國。第一作者：陳杭（1989—），男，博士研究生。

國家自然科學(xué)基金項目（21206038）。