葉思施，唐巧，喬軍帥，王運東（清華大學化學工程系，化學工程聯合國家重點實驗室，北京 100084）

P507-煤油體系物性測量及其在澄清槽內的CFD模擬

葉思施，唐巧，喬軍帥，王運東
（清華大學化學工程系，化學工程聯合國家重點實驗室，北京 100084）

摘要：體系的界面張力對稀土萃取分離過程的澄清分相有著重要的影響。采用懸滴技術分別測量了不同P507-煤油體系與水以及不同稀土酸溶液在常溫常壓下的界面張力，采用旋轉黏度計對不同P507-煤油體系進行黏度測量，初步考察了有機相內稀土負載量對黏度的影響。討論了煤油稀釋劑及氨水皂化對P507-煤油體系與水以及萃取不同稀土濃度酸溶液界面張力及黏度的影響，發(fā)現煤油體積增加，萃取體系與水和鹽酸的界面張力均呈上升趨勢，且在相同P507-煤油體積比下，皂化度增加會使得界面張力急劇減小。固定萃取劑濃度時，有機相內稀土負載量增加在一定皂化度范圍內對界面張力影響波動不大，然而稀土濃度過低或者皂化度過高都會導致不同程度的乳化現象，且不同稀土元素以及元素數量對界面張力影響并不非常明顯。此外，對不同P507-煤油體系在澄清室內的分相性能進行了模擬比較與討論，油相體積分數分布表明界面張力越小，黏度越大，油水兩相澄清越困難。

關鍵詞：P507-煤油體系；界面張力；黏度；澄清室；流體力學模擬

2015-07-17收到初稿，2015-11-20收到修改稿。

聯系人：王運東。第一作者：葉思施（1990—），女，博士研究生。

Received date: 2015-07-17.

引　言

稀土金屬是元素周期表中鑭系15種元素及鈧、釔共17種元素的總稱，因其獨特優(yōu)異的性質，使之具有諸多其他元素所不具備的光、電、磁、熱等物理和化學特性，被視為重要的戰(zhàn)略資源?；谄涮匦裕⊥两饘?、稀土合金以及由此制成的稀土新材料，在電子、信息、通信、汽車、能源、航天航空等領域都獲得了廣泛應用。值得一提的是，稀土金屬是我國的優(yōu)勢資源，其資源儲量、生產量、消費量和出口量均占世界第一[1-5]。

稀土元素因各元素核外電子結構類似，化學性質相似，所以單一稀土元素分離困難。溶劑萃取技術是實現稀土元素分離的主要方法，其基本原理是利用萃取劑對不同稀土元素萃取能力的差異而實現元素分離。目前，工業(yè)上常用的萃取劑主要有酸性磷類萃取劑[1]，如P507、P204、P538等。本文所選用的P507萃取劑，因其密度較大（接近于水）、黏度大，不能直接用于萃取，為了改進萃取過程的水力學特性，通常都要添加密度小、黏度低且能與萃取劑混溶的稀釋劑，常用煤油作為稀釋劑[6]。而且，P507在非極性溶劑中常以二聚體的形式存在，用于金屬萃取時萃取率不高，所以常需要預先進行皂化。皂化時有機磷酸二聚體被打開，形成一種陰離子表面活性劑，可大大提高萃取率，但過度皂化會引起水相乳化[7]。

混合澄清槽是稀土分離過程中應用最為廣泛的逐級接觸式萃取設備，其擁有諸多優(yōu)點，如操作適應性強，萃取級數增減容易，設備加工容易、操作簡單，設備費用少和對廠房高度要求低等[8-11]。進行萃取的兩相在混合槽內經攪拌混合后，必須在相應的澄清槽內進行澄清分相，才能完成一定的萃取操作或實現連續(xù)的混合-澄清過程。針對澄清室，結合現工業(yè)生產情況，其存在以下幾點不足：①分相效率低，混合室內的過度攪拌造成分散相的小液滴、萃取劑雜質較多，乳化現象普遍，澄清分相困難；②試劑流失嚴重，積液量大，一次性投資成本高；③現工業(yè)所用混合澄清槽，設備的混合與澄清部分體積比約為1:3，可見適當減小澄清室體積可大大改善成本投入[12-14]。

計算流體力學是20世紀50年代以來隨著計算機的發(fā)展而產生的一個介于數學、流體力學和計算機之間的交叉學科，可通過計算機和數值方法求解流體力學的控制方程，對流體力學問題進行模擬和分析。近年來，許多學者采用CFD模擬的方法對單元操作過程，如混合、燃燒、分離等進行了模擬，探究實驗中難以觀測到的實驗現象（如體積分數分布、速度流場、壓降等），并得到了與實驗結果較為吻合的結果[15-19]。2011年，Roohollah等[14]對工業(yè)銅萃取設備的澄清部分進行了研究，采用CFD模擬出澄清室內流動特性等，為工業(yè)銅萃取設備澄清室提供了理論指導依據。2012年，Mohsen等[13]采用CFD模擬方法考察了擋板對澄清室分相性能的影響?？傊?，采用CFD模擬對混合澄清槽澄清室進行流動特性的模擬研究尚不完善，有待深入。

在稀土萃取分離過程中，多相流體系的界面張力對混合和澄清分相過程都有著重要的影響，這是因為在這兩個過程中都涉及液滴的形成、聚并以及破碎，而界面張力和黏度直接參與影響液滴的尺寸分布[20]，不同體系物性不同，則其在混合澄清槽內的混合和澄清性能也大不相同。

綜上所述，關于P507-煤油體系中，稀釋劑煤油及皂化對萃取劑物性的影響，尤其是對流體流動有直接影響的界面張力及黏度等物性的基礎數據測量及研究尚沒有報道。采用CFD對澄清槽內分相過程進行模擬的研究也大多基于簡單體系，而與實際體系的聯系較少。基于此，本文采用懸滴技術分別測量了不同P507-煤油體系與水，未皂化P507-煤油體系與鹽酸在常溫常壓下的界面張力，采用旋轉黏度計對不同P507-煤油體系進行黏度測量。討論了煤油稀釋劑及氨水皂化對P507-煤油體系界面張力及黏度的影響。此外，對不同P507-煤油體系在澄清室內的分相性能進行了模擬比較與討論。

1　實驗材料和方法

1.1實驗裝置

本實驗中采用懸滴法測量液-液兩相界面張力，所采用的視頻光學接觸角測量儀型號為OCAH200（德國Dataphysics公司）。

懸滴法測量表面張力原理[21]：借助于現代的完全數字、計算機化，由攝像機抓取一懸滴的圖像，并將整個圖像數字化。數字化后的圖像由計算機進行圖像處理，測定輪廓的坐標，通過將坐標擬合到描述懸滴輪廓的Bashforth-Adams方程式，可得到毛細管常數。已知界面兩相密度差以及重力加速度情況下，便可計算出界面張力值。

萃取體系黏度的測量則是采用數顯黏度計（型號LVDV-Ⅱ+P，Brookfield Engineering Laboratories股份有限公司）。

旋轉法測黏度原理[22]：儀器由同步電機以穩(wěn)定的速度旋轉，連接刻度圓盤，通過游絲和轉軸帶動轉子旋轉，若轉子未受到液體的阻力，則游絲、指針與刻度盤同速旋轉，指針在刻度盤上指出的讀數為“0”。反之，若轉子受到液體的黏滯阻力，則游絲產生扭矩，與黏滯阻力抗衡最后達到平衡，此時在刻度盤上指示一定的讀數（即游絲的扭轉角），將讀數乘上特定系數即可獲得所測液體黏度。

EDTA滴定方法：該方法適用于稀土氧化物中95%～100%稀土總量的測定。本文中稀土酸溶液濃度均由該法確定。準確稱取待測稀土酸溶液1～5 ml 于25 ml錐形瓶中，加入10 ml六亞甲基四銨緩沖溶液，加入3滴二甲酚橙指示劑，用EDTA（乙二胺四乙酸二鈉）標準溶液滴定至溶液由紫紅變?yōu)榱咙S色，即為滴定終點。稀土溶液濃度可由式（1）計算。

式中，CRE為稀土濃度，mol·L?1；CEDTA為EDTA濃度，mol·L?1，此處為0.001007 mol·L?1；VEDTA為滴定過程中所消耗EDTA體積，ml；VRE為測量時稀土酸溶液的取樣體積，ml。

1.2實驗試劑

實驗所用試劑有工業(yè)級P507、磺化煤油（物性參數見表1）、濃氨水、濃鹽酸、去離子水、乙醇、EDTA標準溶液、六亞甲基四銨、二甲酚橙。其中，濃氨水、濃鹽酸、乙醇、六亞甲基四銨、二甲酚橙均為分析純，EDTA標準液為滴定分析純，購置于北京世紀科博有限公司。

圖1　混合澄清槽的結構Fig.1　Structure of mixer-settler

表1　磺化煤油參數Table 1　Parameters of sulfonated kerosene

有機相的制備：工業(yè)級P507與磺化煤油按體積比進行混合攪拌，制備兩種P507濃度（其一為P507與煤油體積比為3:2，其二為配制P507濃度為0.072 mol·L?1）待用，待混合均勻后分別加入濃氨水進行皂化，可制備未皂化及皂化度分別為10%、20%、30%、40%的有機相。

水相的制備：采用去離子水與濃鹽酸配制pH 為3的鹽酸溶液[2-3]，稀土酸溶液采用稀土氯化鹽（帶結晶水）與pH為3的鹽酸溶液配制。

由于測量實際體系時，要求萃取反應達到平衡狀態(tài)，將所有體系油水兩相比均固定為1:1，常溫常壓下在恒溫振蕩箱中完全反應2 h。

1.3模擬方法

模擬研究所用幾何模型是基于實驗所用混合澄清槽結構（圖1）所建，其三維立體模型如圖2所示。

圖2　CFD模擬的幾何模型Fig.2　Geometrical model for settler

CFD模擬方法包括對一系列控制方程的求解，該研究采用ANSYS軟件FLUENT6.3進行。本文中幾何模型采用混合網格，為提高網格質量，在出口處采用四面體網格，同時澄清室主體采用六面體網格，為更準確描述入口處流體湍動情況，其網格尺寸更為密集，網格總數約為421萬個。

假設油水體系為不可壓縮流體，涉及兩相流動，采用歐拉-歐拉多相流模型，湍流模型采用標準k-ε方程，液滴尺寸采用統(tǒng)一粒徑。求解時，獲得單相收斂結果后加入油相，求解分散相控制方程進而獲得兩相收斂結果。針對計算域內采用一級迎風離散化方法，壓力與體積耦合。本模擬模型中采用速度入口，以壓力出口為邊界條件，殘差達到10?4，且所監(jiān)測參數均達到穩(wěn)定時認為模擬收斂。

2　實驗結果與討論

本實驗中所涉及的非實際萃取體系的密度見表2。

表2　常溫常壓下不同體系的密度Table 2　Density of different P507-kerosene system at room temperature and atmospheric pressure

2.1煤油與皂化對界面張力的影響

P507屬于酸性磷酸酯，未皂化的P507在稀釋劑煤油中以二聚體的形式存在，為了增加有機相的萃取容量，可采用濃氨水進行皂化，皂化后的萃取劑外觀仍顯現出清亮透明，不析出水來，是一個油包水的分散體系，非常穩(wěn)定。萃取劑皂化度越高越易達到飽和，因為皂化度高的萃取劑中二聚分子環(huán)被打開，有效萃取劑增加。因本文中所涉及體系相互溶解度很低或不溶，故在測量時未考慮密度因溶解而產生的極微小變化，直接采用純液體的密度來進行計算測量。

2.1.1煤油與皂化對P507-煤油體系與水界面張力的影響由圖3可以發(fā)現，隨著煤油體積分數的增加，P507-煤油體系與水的界面張力皆呈上升趨勢，且煤油體積分數分別為40%和60%時，皂化萃取劑體系的界面張力受煤油體積影響較為平穩(wěn)。

同時，圖4也驗證了不同皂化度下40%和60%煤油體系的界面張力曲線幾乎重合。圖4中給出了不同皂化度下，P507-煤油體系與水的界面張力變化，不難發(fā)現，隨著皂化度的升高，萃取體系界面張力均急劇下降。

圖3　不同煤油體積分數下P507-煤油體系與水的界面張力Fig.3　Interfacial tensions of P507/kerosene-water systems with different volume fraction of kerosene

圖4　不同皂化度下P507-煤油體系與水的界面張力Fig.4　Interfacial tensions of P507/kerosene-water systems with different saponification rate

2.1.2煤油對未皂化P507-煤油體系與鹽酸界面張力的影響在實際皂化P507萃取稀土過程中，稀土溶液一般配成酸性，萃取過程中會發(fā)生氫離子與稀土離子的交換，銨根離子也會參與稀土離子的交換過程，測量界面張力會因油水兩相發(fā)生反應，導致界面模糊（圖5），兩相密度變化大，故暫時不能直接測量皂化P507-煤油體系與鹽酸的界面張力。

圖5　皂化P507-煤油體系在鹽酸中的液滴Fig.5　Droplet in saponified P507/kerosene-hydrochloric acid system

由圖6可以發(fā)現，隨著煤油體積分數的增加，未皂化P507-煤油體系與鹽酸的界面張力呈上升趨勢，與萃取體系和水界面張力的規(guī)律一致。相同條件下，萃取體系與水的界面張力大于與鹽酸的界面張力，是因為鹽酸溶液中H+濃度高于去離子水，和萃取體系相互作用力相對較大而導致界面張力降低。

圖6　不同煤油體積分數下未皂化P507-煤油體系與鹽酸和水的界面張力Fig.6　Interfacial tensions of unsaponified P507/kerosene-hydrochloric acid systems with different volume fraction of kerosene

2.2P507-煤油體系與稀土酸溶液的界面張力測量

雖然由于酸性溶液會與皂化P507-煤油體系發(fā)生反應導致界面波動無法測量，但考慮到澄清室內分相澄清之前體系處于混合狀態(tài)，故假設混合均勻且達到萃取平衡狀態(tài)，此時可測量酸性溶液與P507-煤油體系的界面張力，同時也可以考察水相中負載稀土量、皂化度、不同稀土種類以及單稀土與多稀土對P507-煤油體系與稀土酸溶液界面張力的影響。為使得考察對象更具代表性，特選取鑭、釤、鉺分別代表輕、中、重稀土金屬。本實驗中所涉及的實際萃取體系的密度見表3，部分體系在萃取過程中乳化嚴重，無法測量。

2.2.1水相中稀土負載量及皂化度對P507-煤油體系與稀土酸溶液界面張力的影響圖7給出了無機相中負載不同濃度稀土時，稀土負載量對P507-煤油體系與酸溶液的界面張力，在該圖中也考察了皂化度對界面張力的影響。由表3以及圖7可以發(fā)現，當稀土濃度較低時，皂化度超過一定數值會使得體系乳化嚴重，導致界面張力無法測量，而當稀土濃度較高時，較高的皂化度也不一定會引起乳化現象。固定皂化度，考察稀土濃度對界面張力的影響可以發(fā)現，隨著濃度的增大，界面張力在小范圍內不斷增加，然而增幅不大。固定稀土濃度，改變皂化度，可以發(fā)現在一定范圍內增加皂化度會使得界面張力相應增加（然而增幅較小），而當超過一定數值，則會發(fā)生突變，這有可能是因為物性之間息息相關，超過一定量的皂化會使得萃取劑的體系物性變化較大，進而影響界面張力。

表3　常溫常壓下不同實際體系萃取平衡后的密度Table 3　Density of different organic and aqueous phases after extraction equilibrium at room temperature and atmospheric pressure

2.2.2輕中重稀土及稀土元素數量對P507-煤油體系與稀土酸溶液界面張力的影響稀土元素分類可分為輕、中、重3種，本文分別以鑭、釤、鉺作為代表。同時也對單元素與三元素體系的界面張力進行了區(qū)別和比較。

圖7　不同濃度稀土溶液及皂化度下P507-煤油體系與酸溶液的界面張力Fig.7　Interfacial tensions of P507/kerosene-hydrochloric acid systems with different concentration of La and saponification rate

圖8中給出了不同稀土種類及稀土數量下，P507-煤油體系與酸溶液的界面張力。由于釤、鉺及三元素溶液濃度較低，故采用低濃度鑭溶液來進行比較，然而由于濃度較低時，鑭溶液乳化嚴重無法進行測量，故同時加入了0.00761 mol·L?1鑭溶液進行比較。由圖8可以發(fā)現，在相同皂化度下，稀土種類對低濃度稀土溶液的界面張力影響不大，稀土元素因各元素核外電子結構類似，化學性質相似，這一點也體現在了界面張力上。正因為如此，單元素與三元素的界面張力區(qū)別也并不非常明顯，相較輕稀土金屬而言，中重稀土金屬在低濃度較高皂化度時并不發(fā)生明顯乳化現象。

圖8　不同稀土種類及稀土數量下，P507-煤油體系與酸溶液的界面張力Fig.8　Interfacial tensions of P507/kerosene-hydrochloric acid systems with different kinds and number of rare earths

2.3煤油與皂化對體系黏度的影響

本實驗中所采用的工業(yè)級P507的黏度約為30.8 mPa·s。表4中給出了純P507在不同皂化度下的黏度范圍。

表4　不同皂化度下，P507的黏度范圍Table 4　Viscosity range of P507 with different saponification rate

圖9給出了不同煤油體積比下，P507-煤油體系的黏度。由圖9可以發(fā)現，煤油黏度較小，增加煤油的體積分數可以使萃取體系的黏度急劇下降。

圖9　不同煤油體積分數下P507-煤油體系的黏度Fig.9　Viscosity of P507/kerosene systems with different volume fraction of kerosene

結合表4與圖10可以發(fā)現，當P507體積分數大于20%時，隨著皂化度的升高，黏度升高，且當P507體積比較大時，黏度增加的幅度較大。這是因為在P507-煤油體系中，濃氨水皂化影響的主要是P507。當P507體積分數為20%時，其黏度基本不受皂化度影響。

圖10　不同皂化度下P507-煤油體系的黏度Fig.10　Viscosity of P507/kerosene systems with different saponification rate

2.4有機相稀土負載量對體系黏度的影響

對萃取平衡后的體系，尤其是有機相的黏度進行了考察。當萃取完全后有機相內會負載一定量的稀土，表5分別給出了不同稀土濃度酸溶液經皂化20%的P507-煤油萃取體系（P507濃度為0.072 mol·L?1，油水兩相體積比為1:1）萃取后，有機相內稀土負載量對黏度的影響，同時也考察了稀土元素種類對黏度的影響。

表5　不同稀土種類及負載量對有機相黏度的影響Table 5　Viscosity of P507-kerosen with different rare earths capacity and species

圖11　不同P507體積分數下的澄清室內油相體積分數分布Fig.11　Contours of volume fraction of organic phase with different volume fraction of P507

由表5可以發(fā)現，在文中所給出的稀土濃度范圍內，有機相稀土負載量以及稀土種類對體系黏度影響不大，相較而言，煤油與皂化對體系黏度影響更大。

2.5澄清室內不同體系的分相性能

界面張力是分離過程中最重要的基礎數據之一[20]，是描述相對作用力的重要影響因素，界面張力也是分子間作用力的一種體現。黏度則是流體黏滯性的一種量度，流體的黏滯性作用，使得流體在流動時受到摩擦阻力和壓差阻力，是流體流動時其內部摩擦現象的一種表示，黏度越大則表示內摩擦力越大。本文針對部分上述P507-煤油體系進行了CFD模擬，固定入口油相體積分數為0.075，圖11給出了不同P507體積比下，澄清室內油相體積分數的分布情況。

由圖11可知，增加煤油體積分數，油相黏度下降，油水兩相界面張力增加，分散帶由寬分布變窄，且增加至一定量后對油相體積分數分布影響趨于穩(wěn)定。

觀察圖12中不同P507體積比下的水相速度矢量圖可以發(fā)現，增加煤油體積，黏滯力減小，水相返流明顯，且各種情況下皆在油相中有不同程度的夾帶，針對其返流現象，應適當減少澄清過程中的湍動，或將湍動限制在有限空間內。

圖12　不同P507體積比下澄清室內水相速度矢量圖Fig.12　Velocity vector diagram of aqueous phase with different volume fraction of P507

圖13　不同皂化度下澄清室內油相體積分數分布（P507體積比為60%）Fig.13　Contours of volume fraction of organic phase with different saponification rate

固定P507-煤油體積比為60%，考察不同皂化率下澄清室內的油相體積分數分布，由圖13發(fā)現，固定P507-煤油體積比后，皂化度不同所引起的界面張力與黏度的變化對油相體積分數分布并沒有明顯的影響。

圖14　不同皂化度下澄清室內水相速度矢量圖（P507體積比為60%）Fig. 14　Velocity vector diagram of organic phase with different saponification rate

事實上，實際體系的流場必然會在不同程度上根據不同操作條件的平衡狀況而有差異。然而，從澄清室內的水相速度矢量圖觀察到了和油相體積分數分布類似的結果，即皂化度對體系物性的改變對澄清室內的流場沒有非常大的影響。這可能是因為在該模擬體系中采取的水相為純水，忽略了反應平衡對流場的作用。圖中可以明顯觀察到水相層存有兩個旋渦，這是因為來流碰到底端壁面形成返流，返流與來流相遇形成大小兩個旋渦。

為了能夠更加直觀地表現出體系物性對分相性能的影響，圖15、圖16分別給出了不同P507體積比以及皂化度下，澄清室水相出口處油相夾帶量。

圖15　不同P507體積比下澄清室水相出口處油相夾帶量Fig.15　Entrainment in water outlet with different volume fraction of P507

圖16　不同皂化度下澄清室水相出口處油相夾帶量Fig.16　Entrainment in water outlet with different saponification rate

當澄清室中體系為純P507時，水相出口處油相體積分數為0.00455，夾帶量約占入口輸入的6%，而當采用煤油作為稀釋劑改善萃取劑體系物性后，夾帶現象明顯得到改善。由圖15可以觀察到，煤油體積含量越高，水相出口處夾帶量越低，分相效果越好。然而體系的選擇主要由實際操作確定，并不能因為其分相效果好就不權衡其萃取效果等，即分相過程應與混合過程相輔相成。

固定P507體積分數為60%，改變皂化度可以發(fā)現水相出口處夾帶量均處于同一量級，且隨著皂化度增加，夾帶量增加，這亦與皂化度增加使得黏度增加，界面張力減小，分相阻力增大相呼應。

總體而言，真實體系的分相性能與簡單體系的分相性能有明顯差別，模擬真實體系需要在簡單體系的基礎上循序漸進，然而在CFD模擬中采用真實體系進行研究勢在必行。

3　結　論

本文對稀土溶劑萃取常用P507-煤油體系進行了基礎數據測量，給出了P507-煤油體系與水以及不同稀土酸溶液在常溫常壓下的界面張力，同時測定了各個體系的黏度，考察了煤油與皂化度對P507-煤油體系與水以及萃取不同稀土濃度酸溶液界面張力以及黏度的影響，指出隨著煤油體積增加，界面張力增加，黏度急劇下降。固定P507與煤油體積比，增加皂化度會使得界面張力急劇減小。有機相內稀土負載量在一定范圍內對界面張力和黏度的影響不大，且由于稀土性質相似，不同稀土元素及元素數量對界面張力影響并不明顯。由澄清室內不同體系的CFD模擬可以發(fā)現，實際體系與簡單體系的分相性能差距大，僅研究簡單體系不足以描述實際工業(yè)分離過程。此外，本文并未對P507-煤油體系與稀土水溶液的界面張力進行測量，針對這部分工作本課題組將會結合實際體系的實驗流場測定一同開展。

References

[1] 徐光憲, 袁承業(yè). 稀土的溶劑萃取[M]. 北京: 科學出版社, 1987. XU G X, YUAN C Y. Solvent Extraction of Rare Earths [M]. Beijing: Science Press, 1987.

[2] 夏潔, 任鐘旗, 王運東. P507鹽酸體系萃取分離鈰和鐠[J]. 中國科技論文, 2012, 7(12): 907-912. XIA J, REN Z Q, WANG Y D. Solvent extraction separation of cerium and praseodymium from hydrochloric solution with P507 [J]. China Sciencepaper, 2012, 7(12): 907-912.

[3] 劉晶, 王運東. 皂化P507-HCl-煤油體系萃取鐠釹[J]. 化工學報, 2014, 65(1): 264-270. LIU J, WANG Y D. Extraction of praseodymium(Ⅲ) and neodymium(Ⅲ) in saponified P507-HCl-kerosene system [J]. CIESC Journal, 2014, 65(1): 264-270.

[4] 庾晉, 子蔭. 稀土材料的重要應用[J]. 有色設備, 2002, (6): 17-20. YU J, ZI Y. The important application of rare earth materials [J]. Nonferrous Equipment, 2002, (6): 17-20.

[5] 程建忠, 車麗萍. 中國稀土資源開采現狀及發(fā)展趨勢[J]. 稀土, 2010, 31(2): 65-69. CHENG J Z, CHE L P. The exploitation status and growing trend of rare earths in China [J]. Chinese Rare Earths, 2010, 31(2): 65-69.

[6] SILEM A, 聶國林. 未磺化煤油稀釋劑對從溶液中萃取鈾和可共萃取雜質的影響(Ⅱ): 工業(yè)磷酸[J]. 濕法冶金, 1990, 4: 48-51. SILEM A, NIE G L. The influence of non-sulfonated kerosene diluent to uranium and coextraction impurities by solvent extraction (Ⅱ): Industrial phosphoric acid [J]. Hydrometallurgy, 1990, 4: 48-51.

[7] 吳文遠, 張豐云, 尹少華, 等. 皂化P507體系有機相負載硅對稀土La萃取的影響[J]. 東北大學學報（自然科學版）, 2013, 34(5): 679-682. WU W Y, ZHANG F Y, YIN S H, et al. Effect of saponified P507 loaded with silicon on the extraction of La(Ⅲ) [J]. Journal of Northeastern University (Natural Science), 2013, 34(5): 679-682.

[8] WANG S H, ZHANG T A, ZHAO Q Y, et al. Experimental study on aqueous phase entrainment in a mixer-settler with double stirring mode [J]. China Petroleum Processing and Petrochemical Technology, 2013, 15(2): 59-62.

[9] 劉作華, 唐巧, 王運東, 等.剛柔組合攪拌槳增強混合澄清槽內流體宏觀不穩(wěn)定性[J]. 化工學報, 2014, 65(1): 78-86. LIU Z H, TANG Q, WANG Y D, et al. Enhancement of macro-instability in mixer-settler with rigid-flexible impeller [J]. CIESC Journal, 2014, 65(1): 78-86.

[10] 趙秋月, 張廷安, 劉燕, 等.攪拌對箱式混合澄清槽流動性能的影響[J]. 東北大學學報(自然科學版), 2012, 33(4): 559-562. ZHAO Q Y, ZHANG T A, LIU Y, et al. Effect of stirring on flow characteristic of box-type mixing/clarfying tank [J]. Journal of Northeastern University (Natural Science), 2012, 33(4): 559-562.

[11] 劉海旺, 張麗華. 稀土分離用混合澄清器研究[J]. 稀土, 2000, 21(6): 58-61. LIU H W, ZHANG L H. Study on mixer-settler for rare earth separation [J]. Chinese Rare Earths, 2000, 21(6): 58-61.

[12] GRAEME L, KRISHNA M, WILLIAM Y. Assessment of the flow pattern in a solvent extraction settler [C]// 9th International Conference on CFD in the Minerals and Process Industries. CSIRO, Melbourne, Australia, 2012.

[13] MOHSEN O S, MAZAHERY A, ALIZADEH M, et al. Computational fluid dynamics (CFD) simulation of effect of baffles on separation in mixer settler [J]. International Journal of Mining Science and Technology, 2012, 22(5): 703-706.

[14] ROOHOLLAH S, MOHEBBI A, SARRAFI A, et al. CFD simulation and optimization of the settler of an industrial copper solvent extraction plant: a case study [J]. Hydrometallurgy, 2011, 106(3/4): 148-158.

[15] WEI Z J, XU S M, YUAN Y J, et al. CFD simulation of hydrodynamics characteristics in stirred reactors equipped standard Rushton or 45°-upward PBT impeller [J]. Chinese Journal of Chemical Engineering, 2003, 11(4): 467-471.

[16] ZHANG L H, GAO G H, SUI H, et al. CFD simulation and experimental validation of fluid flow in pre-distributor [J]. Chinese Journal of Chemical Engineering, 2011, 19(5): 815-820.

[17] SINGH K K, MAHAJANI S M, SHENOY K T, et al. CFD modeling of pilot-scale pump-mixer: single-phase head and power characteristics [J]. Chemical Engineering Science, 2007, 62(5): 1308-1322.

[18] SRILATHA C, MUNDADA T P, PATWARDHAN A W. Scale-up of pump-mix mixers using CFD [J]. Chemical Engineering Research and Design, 2010, 88(1): 10-22.

[19] CHEN J B, LIU C J, YUAN X G, et al. CFD simulation of flow and mass transfer in structured packing distillation columns [J]. Chinese Journal of Chemical Engineering, 2009, 17(3): 381-388.

[20] WANG K, Lü Y C, XU J H, et al. Determination of dynamic interfacial tension and its effect on droplet formation in the T-shaped microdispersion process [J]. Langmuir, 2009, 25(4): 2153-2158.

[21] 宋瑛, 田宜靈, 肖衍繁, 等. 二元液液系統(tǒng)界面張力[J]. 化工學報, 1999, 50(5): 620-628. SONG Y, TIAN Y L, XIAO Y F, et al. Interfacial tension of binary liquid-liquid systems [J]. Journal of Chemical Industry and Engineering (China), 1999, 50(5): 620-628.

[22] 劉大紅, 趙慶, 李民. 旋轉黏度計理論公式的探討[J]. 石油鉆采工藝, 1994, 16(3): 40-42. LIU D H, ZHAO Q, LI M. The theoretical formula for rotational viscometer [J]. Oil Drilling Process, 1994, 16(3): 40-42.

DOI：10.11949/j.issn.0438-1157.20151142

中圖分類號：TQ 028.8

文獻標志碼：A

文章編號：0438—1157（2016）02—0458—11

基金項目：國家重點基礎研究發(fā)展計劃項目（2012CBA01203）。

Corresponding author:WANG Yundong, wangyd@tsinghua.edu.cn supported by the National Basic Research Program of China (2012CBA01203).

Physical properties measurements and CFD simulations in settler of different P507-kerosene systems

YE Sishi, TANG Qiao, QIAO Junshuai, WANG Yundong
(State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China)

Abstract:Interfacial tensions play an important role on settling in rare earth element separation process. The interfacial tensions of binary liquid-liquid systems (P507-kerosene / water) and (P507-kerosene / hydrochloric acid which contain rare earths) were measured by the pendant drop method. Viscosity of P507- kerosene with different proportions was obtained by rotational viscometer. The influence of the volume of kerosene, the saponification rate, the concentration of rare earths, the kinds of rare earths and number of rare earths on the interfacial tension and viscosity have been discussed that increase the volume of kerosene would raise the interfacial tensions and lower the viscosity, and the increase of saponification rate brought the decrease of interfacial tension and increase of viscosity. The varied concentration of rare earths did not show great impact on the interfacial tension within a certain degree of saponification rate. Too low concentration and too high saponification rate would raise emulsification, which would hinder the measurements. Different kinds and number of rare earths presented little difference on interfacial tensions for their similar properties. CFD simulations in settler with different P507-kerosene systems have been conducted and the results showed that the phase separation was tougher withlower interfacial tension and higher viscosity.

Key words:P507-kerosene systems; interfacial tension; viscosity; settler; CFD simulation