不同含水率對(duì)油水體系中液滴分布預(yù)測(cè)模型的影響

王玥 鞏孟 李博

摘 要：油水分散系中，求解分散相液滴的分布有十分重要的研究意義。目前一般通過求解群體平衡方程（PBE）來計(jì)算分散相的分布，常用兩區(qū)域模型（葉輪區(qū)域和環(huán)流區(qū)域）進(jìn)行求解，但在求解過程中發(fā)現(xiàn)隨著油包水分散體系中含水率的升高，預(yù)測(cè)結(jié)果偏離實(shí)驗(yàn)結(jié)果。本文通過采用計(jì)算流體力學(xué)（CFD）的方法，用商業(yè)軟件FLUENT模擬發(fā)現(xiàn)，油包水分散系中隨著含水率的升高高能量耗散區(qū)（葉輪區(qū)）在減小，水包油分散系中隨著含水率的升高高能量耗散區(qū)（葉輪區(qū)）在增大。

關(guān) 鍵 詞：油水分散體系；群體平衡方程；計(jì)算流體力學(xué)

中圖分類號(hào)：TQ 028 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1671-0460（2017）07-1402-03

The Impact of Different Water Contents of Oil and

Water System on Predicting Model

WANG Yue1， GONG Meng2， LI Bo3

（1. Tianjin Fire Research Institute， Tianjin 300381， China； 2. Trans-Asia Gas Pipeline Co.， Ltd.， Beijing 100007， China；

3. Hualu Engineering&Technology Co.，Ltd.， Shaanxi Xian 7100651， China）

Abstract： In oil and water system， it has a very important research significance to solve the distribution of dispersed phase droplet. In general， Population Balance Equation （PBE） is used to calculate the distribution of droplet， two-region model （turbine region and circulation region） is used to solve the equation. But， with increasing of the water content of oil-in-water system， predicted results deviate farther from the experimental results. In this paper， based on computational fluid dynamics （CFD）， the commercial software FLUENT was used to simulate the oil and water system. The results show that， with increasing of water content， the high energy dissipation region decrease in water-in-oil system， but increase in oil-in-water system.

Key words： Oil and water system； PBE； CFD

1 背景介紹

目前，油水分散系主要分為兩類，即油包水分散系和水包油分散系。在工業(yè)生產(chǎn)中，對(duì)分散系的形態(tài)和分布進(jìn)行監(jiān)測(cè)控制具有十分重要的意義。油水分散體系是一種以油、水為介質(zhì)的分散體系，這種分散體系要能穩(wěn)定存在，就必須有一定的外部能量來克服形成分散顆粒的界面能，因此，油水分散體系趨向于小的顆粒絮凝聚結(jié)以至于出現(xiàn)相界面分層是一種自發(fā)過程，而顆粒分散并穩(wěn)定存在是一種非自發(fā)過程。一般來說，高速攪拌和添加乳化劑均可得到這種分散體系，通常稱這種分散體系為乳狀液，而攪拌罐的流場(chǎng)內(nèi)很容易形成這種乳狀液。攪拌罐內(nèi)顆粒的粒度和分布特征，一直是諸多學(xué)者研究的重點(diǎn)內(nèi)容。對(duì)于沒有化學(xué)反應(yīng)、熱傳導(dǎo)和質(zhì)量傳遞的系統(tǒng)內(nèi)，體系內(nèi)的液滴在穩(wěn)定的湍流流場(chǎng)內(nèi)是不均勻分布的，并且這些液滴是橢圓形的，其直徑的變化僅與液滴的破裂和聚并相關(guān)。群體平衡方程（Population Balance Equation，以下簡(jiǎn)稱PBE）的一般形式是由Tsouris和Tavlarides[1]提出的。

（1）

其中n（v，t）是在體積為v的液滴t時(shí)刻在系統(tǒng)內(nèi)的數(shù)量，DB和BB代表體積為v的液滴因破裂而消失和出現(xiàn)的概率，相反，DC和BC代表體積為v的液滴因聚并而消失和出現(xiàn)的概率。Coulaloglou和Tavlarides[2]首先提出用PBE模型來預(yù)測(cè)液液分散體系中的不同直徑的液滴分布。Prince和Blanch[3]提出他們的模型，他們認(rèn)為液滴和湍流漩渦的碰撞使得液滴破裂?？紤]到攪拌罐內(nèi)非理想環(huán)境，Alopaeus et al[4]提出一個(gè)多分塊的管內(nèi)模型，此模型將液滴PBE運(yùn)用到濃分散相中。Hu et al提出兩區(qū)域模型（葉輪區(qū)域和環(huán)流區(qū)域）來預(yù)測(cè)攪拌體系內(nèi)的反相問題。中國(guó)石油大學(xué)（北京）王瑋[5]和程為等提出利用兩區(qū)域模型（葉輪區(qū)域占總體積的20%，環(huán)流區(qū)域占80%）來預(yù)測(cè)液滴直徑分布，其模型認(rèn)為輸入功率大部分耗散在葉輪區(qū)即高能量耗散區(qū)，而在環(huán)流區(qū)能量耗散較少為低能量耗散區(qū)，并先后發(fā)表葉輪轉(zhuǎn)速和含水率對(duì)液滴直徑分布的影響，當(dāng)預(yù)測(cè)結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比時(shí)發(fā)現(xiàn)，隨著油包水型乳狀液中含水率的升高，預(yù)測(cè)出的沙德半徑高出實(shí)驗(yàn)數(shù)據(jù)的部分越來越大，當(dāng)擴(kuò)大葉輪區(qū)域體積時(shí)，預(yù)測(cè)結(jié)果接近實(shí)驗(yàn)結(jié)果。最近一些年，有很多人在研究油水分散系的攪拌和混合問題中，用到了計(jì)算流體力學(xué)（CFD）的軟件，如FLUENT、CFX等。計(jì)算流體力學(xué)可以預(yù)測(cè)流體動(dòng)力和湍流動(dòng)能，以及對(duì)裝置內(nèi)流體的局部特點(diǎn)有很完整的解決方法。粒子總體平衡模型（PBM）可以預(yù)測(cè)液滴的運(yùn)動(dòng)、破裂和聚結(jié)，以及分散相液滴的直徑分布。但這個(gè)方法在單獨(dú)運(yùn)用時(shí)，有其限制條件。本文就是利用FLUENT商業(yè)軟件對(duì)攪拌罐進(jìn)行模擬，從而研究不同分散相含量對(duì)于體系流場(chǎng)的影響，主要關(guān)注其對(duì)整個(gè)體系的能量耗散率的影響。

2 建立模型

2.1 實(shí)驗(yàn)條件介紹

本實(shí)驗(yàn)中采用的油為Exxon Mobil出售的Exxsol D80溶劑油（以下簡(jiǎn)稱D80）。Exxsol D80是一種脫芳烴、脫硫的溶劑油，無色無味透明液體，不易揮發(fā)。水采用的是經(jīng)ELGA純水儀處理后的去離子水（電阻率保持在15 MΩ）。油相和水相的主要物性見表1。

本實(shí)驗(yàn)中采用的攪拌設(shè)備為IKA RW 20型攪拌器，配四葉槳式攪拌器，攪拌槳的葉輪直徑為50 mm，葉輪的厚度為1 mm，葉片傾斜角度為45°。攪拌轉(zhuǎn)速范圍為0～2 000 r/min。

2.2 建模過程

首先，根據(jù)實(shí)驗(yàn)中所用到的攪拌器和燒杯的尺寸在SOLIDWORKS軟件中建立幾何模型，模型的剖面圖如下圖1所示。

然后，利用ICEM CFD軟件在建立好的幾何模型中劃分網(wǎng)格，網(wǎng)格采用非結(jié)構(gòu)四面體網(wǎng)格，最大尺寸為2 mm，葉輪附近的環(huán)形區(qū)域?yàn)榧用軈^(qū)，最大尺寸為0.5 mm，最終共劃分了104萬個(gè)四面體網(wǎng)格，結(jié)果如下圖2所示。

3 FLUENT模擬計(jì)算及結(jié)果

將劃分好的網(wǎng)格在FLUENT 14.0中進(jìn)行模擬計(jì)算，流體條件按照實(shí)驗(yàn)條件設(shè)置，湍流模型采用標(biāo)準(zhǔn) 模型（壁面采用標(biāo)準(zhǔn)壁面函數(shù)），用多參考系模型來描述葉輪的旋轉(zhuǎn)運(yùn)動(dòng)，整個(gè)體系壁面為無滑移條件，油水體系采用混合模型進(jìn)行計(jì)算，分別模擬不同含水率的油水乳狀液（水的體積分?jǐn)?shù)為10%，20%，30%），固定旋轉(zhuǎn)速度為700 r/min，主要分析體系內(nèi)能量耗散率的分布情況，曲線主要是取Y=0，X>0截面上距離葉輪邊緣0、1、2、3…15 mm處的能量耗散率的值沿Z軸的變化情況，具體結(jié)果圖3-5所示。

為了可以直觀的看出不同含水率的體系內(nèi)在能量耗散率的不同，并體現(xiàn)在云圖中，將油品的粘度擴(kuò)大10倍繼續(xù)上述模擬就可以得出能量耗散率云圖，如圖6所示。

同樣的方式可以得出含油體積分?jǐn)?shù)為10%，20%，30%情況下的能量耗散率云圖，如圖7所示。

4 結(jié)果分析

通過對(duì)比粘度未被放大的情況下的能量耗散率曲線發(fā)現(xiàn)，隨著含水率的增加葉輪區(qū)域的能量耗散率在不斷減小，之后對(duì)比粘度放大的情況下的能量耗散率云圖可以發(fā)現(xiàn)，隨著油包水型乳狀液中含水體積分?jǐn)?shù)的增加，高耗散率區(qū)域即葉輪區(qū)域在不斷縮小，而低能量耗散區(qū)即環(huán)流區(qū)域在不斷增大；與此相反，隨著水包油型乳狀液中含油體積分?jǐn)?shù)的增加，高耗散率區(qū)域即葉輪區(qū)域在不斷擴(kuò)大，而低能量耗散區(qū)即環(huán)流區(qū)域在不斷減小。

故得出結(jié)論：（1）整個(gè)體系內(nèi)能量耗散率不是均勻分布的，即使在高能量耗散區(qū)內(nèi)也是有差別的。（2）區(qū)域模型計(jì)算PBE時(shí)，葉輪區(qū)域和環(huán)流區(qū)域的比例對(duì)于不同油水比例的體系是不固定的，其中水的比例升高時(shí)應(yīng)減小葉輪區(qū)域，油的比例增加時(shí)要增加葉輪區(qū)域。（3）在計(jì)算油水乳狀液中液滴分布時(shí)，需要將PBE與流體流場(chǎng)的計(jì)算相互偶合起來才能更加接近實(shí)驗(yàn)數(shù)據(jù)。

參考文獻(xiàn)：

[1]Tsouris， C.， Tavlarides， L.L. Breakup and coalescence models for drops in turbulent dispersions [J]. A. I. Ch. E. Journal， 1994， 40 （3）： 395-406.

[2]Coulaloglou， C.A.， Tavlarides， L.L. Description of interaction processes in agitated liquid-liquid dispersions [J].Chemical Engineering Science， 1977， 32：1289-1297.

[3]Prince， M.， Blanch， H.W. Bubble coalescence and breakup in air-sparged bubble columns[J]. A. I. Ch. E. Journal， 1990， 36 （10）： 1485-1499.

[4]Alopaeus V， Koskinen J. Simulation of the population balances for liquid-liquid systems in a non-ideal stirred tank[J]. II Chem.Sci.Eng.， 2002， 57： 1815-1825.

[5] 王瑋.油水混合液物性及流動(dòng)規(guī)律特性研究（D）.北京：中國(guó)石油大學(xué)，2009.