李苗苗，徐 俊，王 新，徐 慶

（天津科技大學(xué) 機(jī)械工程學(xué)院；天津市輕工與食品工程機(jī)械裝備集成設(shè)計(jì)與在線監(jiān)控重點(diǎn)實(shí)驗(yàn)室，天津300222）

液滴在固體表面的流動(dòng)變化，歷來是工程領(lǐng)域的一個(gè)研究熱點(diǎn)，張建林等[1]就對(duì)液滴在固體表面的鋪展情況進(jìn)行了研究。流化床包覆是指芯材在某一溫度流化氣體的流化作用下處于流化狀態(tài)，霧化的壁材與芯材接觸將芯材包埋的技術(shù)，也是微膠囊常用的制備技術(shù)[2]。包覆的壁材既可以是親水性的糖類或蛋白類，又可以是親脂性材料。Litster等[3]認(rèn)為包覆的過程主要分為3個(gè)步驟：1）液滴的潤濕和成核；2）顆粒的固化和生長；3）顆粒間磨損破碎，該理論也被用來解釋流化床包覆過程。Link等[4]總結(jié)了影響包覆效果和包覆過程的影響因素，主要分為液體特性、固體特性、系統(tǒng)和操作變量。

在頂噴的流化床包覆過程中，主要包括氣體的流化作用、液滴的霧化作用、顆粒與液滴之間的碰撞以及液滴在顆粒表面的沉積鋪展等等。霧化液滴在流化床內(nèi)隨機(jī)運(yùn)動(dòng)進(jìn)行碰撞，碰撞過程中，如果2個(gè)表面潤濕的顆粒接觸，會(huì)形成液橋，顆粒在冷空氣的作用下，固化會(huì)導(dǎo)致液橋向固橋轉(zhuǎn)變，形成團(tuán)聚。顆粒的團(tuán)聚和破碎經(jīng)歷很多微觀過程，需要引進(jìn)群體平衡模型（population balance model，PBM）[5]計(jì)算顆粒的一些變量，如粒徑分布、潤濕顆粒的數(shù)量等。

PBM的模擬對(duì)象主要是多相流中副相的粒徑分布情況，如液滴、氣泡以及固態(tài)顆粒等，需要得到合適的粒徑分布規(guī)律[6-9]。計(jì)算粒徑分布主要有2種方法，一種是拉格朗日法，即離散相模型（discrete phase model，DPM），另一種是歐拉法，即PBM法。由于DPM主要模擬離散相，且僅限于稀疏相的模擬，群體平衡模型因?yàn)橄∈柘嗑梢阅M而被廣泛提出。

本文中基于PBM和DPM，對(duì)噴霧冷凍流化床包覆過程中的霧化液滴的產(chǎn)生、凝固、運(yùn)動(dòng)進(jìn)行了數(shù)值模擬，主要對(duì)霧化過程中的速度、溫度及顆粒的團(tuán)聚破碎規(guī)律進(jìn)行了研究。

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

在過去的幾十年時(shí)間里，計(jì)算流體動(dòng)力學(xué)（computational fluid dynamics，CFD）已被廣泛應(yīng)用于噴霧領(lǐng)域，比如霧滴的沉積特性研究[10]、霧化效果的研究[11]等，各種數(shù)值方法，例如流體體積法、DPM已經(jīng)用于預(yù)測(cè)噴嘴的霧化特性和預(yù)測(cè)液滴的軌跡[12]。DPM方法是基于拉格朗日計(jì)算方法[13]，液滴軌跡使用單獨(dú)的控制方程計(jì)算；而體積流體應(yīng)用歐拉法，把各相看作是相互聯(lián)系的相計(jì)算。模擬粒子過程主要有PBM和DPM兩大體系[14]。

1.1 氣相控制方程

計(jì)算流體力學(xué)需要了解流體的動(dòng)力特性，這對(duì)研究霧化過程中的流場(chǎng)分布有著很重要的意義，因此，為了得出流體的流速與壓力等動(dòng)力學(xué)特性，在計(jì)算中，將空氣作為連續(xù)相介質(zhì)處理，采用SIMPLE算法求解，方程的通用形式表示為

式中：Φ為通用變量；ΓΦ為廣義擴(kuò)散系數(shù)；SΦ為廣義源相。

1.2 群體平衡模型

霧化液滴在運(yùn)行的過程中會(huì)出現(xiàn)，聚并和破碎的行為，這在研究流化床包覆的過程中也是必須要考慮的，關(guān)于顆粒的聚并和破碎的研究，在實(shí)驗(yàn)中難以控制和實(shí)現(xiàn)，因此本研究中引入了群體平衡模型。群體平衡模型是一組描述顆粒聚并和破碎行為的偏微分方程[15-16]，反映顆粒成長和消亡的過程。其一般表示形式為

1.3 離散相模型

利用DPM可以對(duì)顆粒進(jìn)行追蹤，因?yàn)樗窃跉W拉-拉格朗日法的基礎(chǔ)上建立的，所以認(rèn)為液滴或顆粒是控制體內(nèi)動(dòng)量、質(zhì)量和能量的源項(xiàng)。

由流化床頂部霧化器生成液滴的直徑是一組離散的尺寸，并且每種尺寸的液滴都各有其運(yùn)動(dòng)軌跡，沿著該軌跡的液滴的流動(dòng)速度是固定的。在設(shè)定了空氣的溫度和速度后，各組尺寸液滴的控制方程就會(huì)進(jìn)行迭代，計(jì)算它們的速度、軌跡、尺寸和溫度，直到這些液滴到達(dá)噴霧流化床的出口或碰撞到壁面。

追蹤顆粒軌跡是通過計(jì)算力平衡方程實(shí)現(xiàn)的，力平衡方程為：

式中：v為連續(xù)相速度；up為顆粒速度；ρga為連續(xù)相密度；ρp為顆粒密度。

作用在顆粒上的力有重力、慣性力以及氣體曳力。

滑移雷諾數(shù)（Re）和曳力系數(shù)（CD）的計(jì)算公式為

式中：dp為顆粒直徑；α1、α2、α3為顆粒球形度系數(shù)。

2 模型的建立

本幾何模型對(duì)噴霧流化床的模型[17]進(jìn)行簡化，計(jì)算域如圖1所示。

圖1 幾何模型尺寸和網(wǎng)格Fig.1 Size of geometric model and mesh

本研究中模型為柱錐形流化床，計(jì)算域使用二維空間進(jìn)行模擬，幾何體相關(guān)尺寸參數(shù)及模型參數(shù)見表1。表中所列數(shù)據(jù)是參考了耿縣如等[17]的單液滴噴霧冷凍研究的有關(guān)數(shù)據(jù)，并進(jìn)行了相關(guān)修改。

3 計(jì)算結(jié)果及討論

3.1 液滴的溫度

圖2為液滴溫度隨時(shí)間的變化圖。

表1 邊界條件和參數(shù)Tab.1 Boundary conditions and simulation parameters

圖2 液滴溫度隨時(shí)間的變化Fig.2 Change of droplet temperature with time

由圖可知，液滴從入口噴出，在噴出時(shí)液滴溫度為300 K，由于環(huán)境溫度很低，為203 K，因此液滴在達(dá)到底部時(shí)已成凝固狀態(tài)。整個(gè)過程中，流化床內(nèi)溫度在X軸兩側(cè)呈對(duì)稱分布，經(jīng)過1.8 s到達(dá)錐體體部分，隨即底部溫度達(dá)到最大，繼而產(chǎn)生回流，之后流化床內(nèi)溫度又呈對(duì)稱分布。整個(gè)過程當(dāng)中，只有1.8 s到2.34 s底部溫度達(dá)到最大，到達(dá)5.85 s以后流化床內(nèi)溫度在X軸兩側(cè)又呈對(duì)稱分布，這是因?yàn)殪F化液滴在運(yùn)動(dòng)過程中，中心區(qū)域流速較大，所以會(huì)很快達(dá)到底部，而當(dāng)整個(gè)過程達(dá)到穩(wěn)定的過程中，霧化液滴與冷空氣接觸，溫度也逐漸下降，同時(shí)達(dá)到底部又產(chǎn)生回流，最終又呈現(xiàn)出兩邊的對(duì)稱分布。

本研究中以模型液滴入口為坐標(biāo)原點(diǎn)，軸向?yàn)閄方向，指向底部椎體部分的方向?yàn)閄軸正方向。結(jié)果顯示，以X軸為軸線，兩邊成對(duì)稱分布，因此以下圖像只采取了模型的一半進(jìn)行了研究。

圖3為流化床內(nèi)溫度隨Y軸的變化圖。如前所述，霧化液滴的溫度會(huì)對(duì)液滴的霧化產(chǎn)生影響，因?yàn)橐旱卧跍囟冗^低時(shí)會(huì)產(chǎn)生凝固現(xiàn)象。在流化床包覆過程中，一旦液滴涂敷在顆粒表面，是需要進(jìn)行冷卻的，因此本文也研究了流化床內(nèi)的溫度分布。

圖3 流化床內(nèi)溫度隨Y軸的變化Fig.3 Temperature change with axis Y in a fluidized bed

圖3說明了霧化液滴在不同軸向位置上的溫度從Y=0 m處到壁面處的溫度變化，冷氣入口接近液滴噴射入口。由于冷氣初始溫度和幾何模型壁面溫度都均為203 K，而液滴溫度為常溫300 K，因此液滴在與冷氣接觸后，溫度會(huì)很快下降。因?yàn)榱骰仓行牟课皇菤怏w的噴射區(qū)，所以中間部位的溫度很低，由中心向兩邊溫度逐漸升高，而在壁面處，由于液滴的回流，會(huì)使得液滴能量損失，溫度與壁面溫度接近，因此出現(xiàn)流化床內(nèi)接近壁面處溫度陡降，這與耿縣如等[17]所做關(guān)于液滴噴霧冷凍的溫度變化和分布規(guī)律相類似。

3.2 霧化液滴的運(yùn)動(dòng)

圖4為液滴在氣液兩相流的霧化特性圖。由圖4a可以看出，液滴從入口進(jìn)入流化床內(nèi)，液滴均具有一定初速度，因此液滴會(huì)很快到達(dá)流化床底部壁面，碰到壁面后，由于沖量作用，其方向發(fā)生了改變，這就會(huì)使得碰到流化床底部的液滴產(chǎn)生向上的回流，如圖4b所示。同時(shí)使得液滴的速度呈現(xiàn)出沿軸向逐漸減小，沿徑向中心高兩邊低的現(xiàn)象，如圖4c所示。這說明在液滴霧化的過程中，由于幾何模型結(jié)構(gòu)的作用，會(huì)產(chǎn)生一定的回流這將對(duì)液滴霧化產(chǎn)生影響。同時(shí)比較圖4c和圖4d也可以看出，PBM模型和DPM模型在計(jì)算液滴霧化過程中有著相同的速度場(chǎng)分布，這與Taghipour等[13]對(duì)氣固流化床的研究結(jié)論是一致的。同時(shí)PBM又可以用于液滴的粒徑分布研究，因此本文采用了PBM模型來研究液滴的粒徑變化和分布情況。

圖4液滴在氣液兩相流的霧化特性Fig.4 Atomization characteristics of droplets in gas-liquid two-phase flow

3.3 液滴的速度

圖5為液滴X方向速度隨Y坐標(biāo)的變化圖。

霧化液滴之間會(huì)出現(xiàn)團(tuán)聚和破碎的現(xiàn)象，這直接關(guān)乎液滴的霧化效果，同時(shí)也會(huì)對(duì)流化床包覆的研究產(chǎn)生影響。速度和溫度的分布對(duì)霧化液滴的團(tuán)聚和破碎有著很大影響。

圖5a、5b說明了PBM和DPM模型在對(duì)液滴霧化進(jìn)行計(jì)算過程中，在軸向上不同位置液滴軸向速度從Y=0處向壁面接近的變化。液滴的初速度為25 m/s，液滴噴出后，在幾何模型的中部，液滴有能量損失，因此圖中液滴在X=0 m處速度為23 m/s，并且越接近幾何體壁面，液滴的速度越小。在Y=0.21 m處液滴速度為零，從Y=0.21 m處到壁面即Y=0.4 m處，液滴速度方向改為向上，這是因?yàn)橐旱卧谂龅奖诿嬉院蟪霈F(xiàn)了回流現(xiàn)象，回流又會(huì)導(dǎo)致液滴之間的接觸和碰撞。這一現(xiàn)象將會(huì)使液滴的粒徑發(fā)生變化，同時(shí)由圖可知，2種模型的計(jì)算結(jié)果吻合，這也說明了由圖4c、d得出的可以用PBM來研究液滴的粒徑變化和分布情況的結(jié)論。

3.4 液滴的粒徑分布

本文中粒徑分布是指液滴在某種粒徑下的液滴數(shù)量，其具體分布情況如圖6所示。

圖6a、b、c表明了霧化液滴軸向上不同位置的粒徑分布，液滴的初始粒徑為1.26 μm。由圖可知，液滴在霧化過程中發(fā)生了團(tuán)聚和破碎現(xiàn)象，形成了粒度不一的分布。團(tuán)聚體數(shù)量最大的粒徑范圍為2.4～3 μm，最小范圍為0～0.5 μm。在軸向上，離噴口越遠(yuǎn)處，液滴的團(tuán)聚現(xiàn)象越明顯，即2.4～3 μm的液滴數(shù)量越多，而0～0.5 μm、0.5～0.79 μm的液滴數(shù)量比較少。這是因?yàn)樵诹骰驳闹行牟课粴怏w的溫度，速度較大，而隨著液滴向流化床邊緣不斷的運(yùn)動(dòng)，液滴出現(xiàn)了團(tuán)聚和破碎的現(xiàn)象如圖6d所示，ddroplet表示液滴的粒徑。由液滴的累積粒度曲線可以看出，液滴在噴口處和末端部位的數(shù)量是最多的，這主要是因?yàn)樵诹骰岔敳拷咏肟谔幰旱嗡俣容^大，而在接近底部椎體部分液滴的碰撞比較明顯，碰撞后形成的小顆粒數(shù)量較多。

圖5 液滴X方向速度隨Y坐標(biāo)的變化Fig.5 Change of X direction velocity of droplet with Y coordinate

圖6霧化液滴的粒徑分布Fig.6 Particle size distribution of atomized droplets

圖6b表示的是X=1.37 m直線上的液滴在徑向由中心向兩邊運(yùn)動(dòng)的過程中，1.8～2.4、0.6～1.2 μm粒徑的液滴數(shù)量的變化。由圖可知，由于初始粒徑為1.26 μm，在Y=0～0.2 m的過程中，液滴同時(shí)存在團(tuán)聚和破碎現(xiàn)象，因此才會(huì)出現(xiàn)了1.8～2.4、0.6～1.2 μm的液滴數(shù)量不斷增加的現(xiàn)象。在Y=0.2 m處以后流化床內(nèi)的溫度開始越來越接近壁面溫度即T=203 K，因此顆粒的溫度開始逐漸降低出現(xiàn)凝固現(xiàn)象，團(tuán)聚后難以破碎，2.5～3 μm的液滴數(shù)量分布比較多，所以1.8～2.4、0.6～1.2 μm的液滴，尤其是1.8～2.4 μm的數(shù)量在這一區(qū)域比較少（說明：霧化的液滴在實(shí)際過程中由于噴嘴的作用，液滴可以被分散成很細(xì)小的單個(gè)液滴，出現(xiàn)幾微米的液滴屬合理現(xiàn)象）。

綜上所述，由模擬結(jié)果可知，液滴在霧化過程中，由于流化床壁面和入口氣體的溫度都為203 K，與液滴初始溫度相差較大，因此接近壁面處流化床內(nèi)溫度會(huì)出現(xiàn)陡降。同時(shí)由于液滴的初速度為28 m/s，因此液滴在短時(shí)間內(nèi)到達(dá)底部椎體部分，產(chǎn)生回流，造成液滴自身的能量損失，而這將會(huì)使接近壁面處的液滴出現(xiàn)大量團(tuán)聚。

4 結(jié)論

本文中基于CFD中的群體平衡模型對(duì)流化床內(nèi)液滴的霧化行為進(jìn)行了二維數(shù)值模擬，并與離散相模型的計(jì)算結(jié)果進(jìn)行對(duì)比，結(jié)果表明，群體平衡模型和離散相模型可以得出相同的速度場(chǎng)分布，群體平衡模型可以進(jìn)一步得出更為詳細(xì)的粒徑分布，故群體平衡模型可以作為離散相模型的替代模型進(jìn)一步研究液滴的聚并和破碎行為。

流化床內(nèi)溫度呈現(xiàn)中間低、兩邊高的狀況，在液滴碰到錐形壁面時(shí)液滴產(chǎn)生回流，液滴在這一區(qū)域很容易團(tuán)聚。在中心區(qū)域液滴速度較大，破碎現(xiàn)象比較明顯。

（

）：

[1]張建林，劉邱祖，賴成，等.液滴在固體表面流動(dòng)特性的數(shù)值模擬及驗(yàn)證[J].中國粉體技術(shù)，2016（3）:19-21.

[2]NIENOWA W.Fluidised bed granulation and coating:applications to materials，agriculture and biotechnology[J].Chemical Engineering Communication，1995，139：233-253.

[3]LITSTER J，ENNIS B.The science and engineering of granulation processes[J].Particle Technology，2004，15.

[4]LINK K C，ERNST-ULRICH S.Fluidized bed spray granulation investigation of the coating process on a single sphere[J].Chemical Engineering and Processing 1997，36：443-457.

[5]MUBASHIR，HUSSAIN，JITENDRA，et al.A new framework for population balance modeling of spray fluidized bed agglomeration[J].Particuology，2015，19（2）:141-154.

[6]REYNOLDS G K.Modelling of pharmaceutical granule size reduction in a conical screen mill[J].Chemical Engineering Journal，2010，164（2/3）:383-392.

[7]GRIFFIN D W，MELLICHAMP D A，DOHERTY M F.Reducing the mean size of API crystals by continuous manufacturing with product classification and recycle[J].Chemical Engineering Science，2010，65（21）:5770-5780.

[8]MORTIER S T F C，GERNAEY K V，BEER T D，et al.Development of a population balance model of a pharmaceutical drying process and testing of solution methods[J].Computers&Chemical Engineering，2013，50（9）:39-53.

[9]BASHEER I A，HAJMEER M.Artificial neural networks:fundamentals，computing，design，and application[J].Journal of Microbiological Methods，2000，43（1）:3.

[10]孫國祥，汪小旵，丁為民.基于CFD離散相模型霧滴沉積特性的模擬分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（6）:13-19.

[11]梁欽，高貴軍，劉邱祖.壓力型霧化噴嘴射流噴霧氣-液兩相流數(shù)值模擬[J].中國粉體技術(shù)，2015，21（2）:5-9.

[12]DUANGKHAMCHAN W，RONSSEWF，DEPYPERE F，et al.CFD study of droplet atomisation using a binary nozzle in fluidised bed coating[J].Chemical Engineering Science，2012，68（1）:555-566.

[13]TAGHIPOUR F，ELLIS N，WONG C.Experimental and computational study of gas solid fluidized bed hydrodynamics[J].Chemical Engineering Science，2005，60（24）:6857-6867.

[14]BARRASSO D，TAMRAKAR A，RAMACHANDRAN R.A reduced order PBM ANN model of a multi-scale PBM DEM description of a wet granulation process[J].Chemical Engineering Science，2014，119（119）:319-329.

[15]侯廣信，陳建義.顆粒碰撞團(tuán)聚對(duì)旋風(fēng)分離器分離性能影響研究[J].石油機(jī)械，2008，36（5）:1-4.

[16]MCGRAWR.Description of aerosol dynamics by the quadrature method of moments[J].Aerosol Science&Technology，1997，27（2）:255-265.

[17]耿縣如，徐慶，李占勇，等.噴霧冷凍法單個(gè)液滴凍結(jié)過程模擬[J].化工進(jìn)展，2012，31（5）:981-986.