蒸汽在含有不可溶核和可溶無機(jī)鹽的細(xì)顆粒物表面的核化特性

蒸汽在含有不可溶核和可溶無機(jī)鹽的細(xì)顆粒物表面的核化特性發(fā)展生物質(zhì)發(fā)電技術(shù)是減少溫室氣體排放的重要手段之一，對(duì)我國推動(dòng)實(shí)現(xiàn)“碳達(dá)峰·碳中和”目標(biāo)具有重要意義。然而，生物質(zhì)直接燃燒和生物質(zhì)與煤混合燃燒發(fā)電過程排放出大量細(xì)顆粒物。這些細(xì)顆粒物表面通常含有一定量的可溶無機(jī)鹽，其中氯化鈉（NaCl）、氯化鈣（CaCl2）和氯化鉀（KCl）占比較多[1]。由于細(xì)顆粒物粒徑微小，目前工業(yè)中應(yīng)用的常規(guī)除塵設(shè)備，如靜電除塵器、袋式除塵器等對(duì)其捕集效果仍較差。鑒于排放到大氣中的細(xì)顆粒物對(duì)人群健康和大氣環(huán)境帶來嚴(yán)重危害，在“雙碳”背景下，針對(duì)含有可溶無機(jī)鹽的細(xì)顆粒物排放控制技術(shù)開展研究非常重要且迫切。對(duì)細(xì)顆粒物進(jìn)行預(yù)處理，使其粒徑增大，以利于常規(guī)除塵設(shè)備效率的提高，是細(xì)顆粒物排放控制的重要技術(shù)途徑。目前，已發(fā)展出聲凝并[2]、電凝并[3]和蒸汽凝結(jié)[4-6]等顆粒預(yù)處理技術(shù)。其中，蒸汽凝結(jié)預(yù)處理技術(shù)與煙氣濕法、半干法脫酸工藝相結(jié)合富有工業(yè)應(yīng)用的潛力，引起了研究者的持續(xù)關(guān)注[7-9]。

異質(zhì)核化是蒸汽在顆粒表面凝結(jié)的必經(jīng)過程，由于大氣氣溶膠顆粒的異質(zhì)核化對(duì)于成云、降水以及霧霾的形成都有著重要影響，異質(zhì)核化也是大氣環(huán)境領(lǐng)域關(guān)注的重點(diǎn)[10]。現(xiàn)有描述過飽和蒸汽在細(xì)顆粒物表面異質(zhì)核化的理論主要有經(jīng)典異質(zhì)核化理論[11]、分子動(dòng)力學(xué)理論[12]、密度函數(shù)理論[13]等。其中，經(jīng)典異質(zhì)核化理論利用經(jīng)典熱力學(xué)方法，基于晶核形成自由能表達(dá)式，求解異質(zhì)核化特性參數(shù)，應(yīng)用該理論研究異質(zhì)核化具有計(jì)算簡單、計(jì)算量小的特點(diǎn)，在晶核尺度較大或蒸汽飽和度較低時(shí)能夠獲得可靠的預(yù)測結(jié)果[5,14-15]。鑒于蒸汽凝結(jié)顆粒預(yù)處理技術(shù)中蒸汽飽和度通常在1.1~1.5，經(jīng)典異質(zhì)核化理論適用于描述該情況下的異質(zhì)核化。經(jīng)典異質(zhì)核化理論[16]首次提出以來，不斷得到發(fā)展，主要體現(xiàn)在：將理論推廣到含有可溶成分[17-18]、表面粗糙[5]、表面不均勻[19]的顆粒，引入了線張力的影響[20]，考慮了晶核增長機(jī)制[21]等。然而，縱觀已有理論分析和數(shù)值模擬研究，對(duì)象集中在不可溶球形顆粒，針對(duì)生物質(zhì)直燃與混燃電站排放的包含可溶無機(jī)鹽的細(xì)顆粒物的研究存在明顯不足。盡管Gorbunov等[17]對(duì)蒸汽在含有球形不可溶核和NaCl的氣溶膠顆粒表面異質(zhì)核化的晶核形成自由能與平衡蒸汽飽和度進(jìn)行了數(shù)值預(yù)測，但由于研究對(duì)象為大氣環(huán)境，研究中采用的蒸汽飽和度接近于1，無法為蒸汽凝結(jié)預(yù)處理技術(shù)脫除細(xì)顆粒物的工程應(yīng)用提供參考。凡鳳仙等[22]在接觸角為80°條件下對(duì)表面含有微量NaCl的細(xì)顆粒物的異質(zhì)核化特性進(jìn)行了數(shù)值預(yù)測，采用的蒸汽飽和度為1.8，顯得過高；Fan等[18]結(jié)合生活垃圾焚燒顆粒的特性，利用數(shù)值模擬方法研究了飽和度為1.45條件下，蒸汽在同時(shí)含有不可溶球形核和NaCl、CaCl2兩種可溶無機(jī)鹽（二者含量接近）的顆粒表面的異質(zhì)核化特性，但缺乏對(duì)有、無可溶無機(jī)鹽時(shí)成核特性的對(duì)比研究，也未對(duì)不同成分可溶無機(jī)鹽的影響進(jìn)行探究。此外，這些研究在計(jì)算成核速率時(shí)引入了經(jīng)驗(yàn)參數(shù)成核動(dòng)能常數(shù)，使得計(jì)算結(jié)果具有不確定性。

針對(duì)已有研究存在的問題，本文綜合考慮晶核生長的直接沉積機(jī)制和表面擴(kuò)散機(jī)制，建立了改進(jìn)的蒸汽在含有球形不可溶核和可溶無機(jī)鹽的細(xì)顆粒物表面的異質(zhì)核化模型，基于該模型對(duì)蒸汽在不同種類細(xì)顆粒物（包括球形不可溶顆粒、表面含有不同可溶無機(jī)鹽的顆粒）表面異質(zhì)核化的晶核形成自由能、成核機(jī)制和臨界飽和度進(jìn)行對(duì)比分析。研究結(jié)果對(duì)于認(rèn)識(shí)和掌握蒸汽在復(fù)雜成分細(xì)顆粒物表面的異質(zhì)核化特性具有重要價(jià)值，對(duì)于利用蒸汽凝結(jié)顆粒預(yù)處理技術(shù)實(shí)現(xiàn)生物質(zhì)電站細(xì)顆粒物的超低排放具有重要意義。

1模型與算法

1.1異質(zhì)核化模型1.1.1幾何參數(shù)的計(jì)算經(jīng)典異質(zhì)核化理論認(rèn)為，蒸汽分子在固體顆粒表面聚集形成球冠形液體晶核，當(dāng)晶核生長到某一臨界尺寸成為臨界晶核后，蒸汽依托臨界晶核持續(xù)凝結(jié)，固體顆粒長大為含塵液滴[11,16]。圖1給出了本文所涉及的球形不可溶顆粒、表面含有可溶無機(jī)鹽的顆粒，以及蒸汽在顆粒表面異質(zhì)核化形成的球冠形晶核的示意圖。圖中，點(diǎn)O為球形不可溶顆?；蚯蛐尾豢扇芎说闹行模琑p為其半徑，O′為氣-液界面所在球面的中心，re為晶核半徑，s為點(diǎn)O到O′的距離，θ為晶核與顆粒之間的接觸角，φ為液-固界面的中心半角，ψ為氣-液界面的中心半角。

圖1

圖1顆粒及蒸汽在顆粒表面異質(zhì)核化形成晶核示意圖

Fig.1Schematicdiagramofparticlesandembryoformedbyheterogeneousnucleationofwatervaporontheparticlesurface

根據(jù)圖1的幾何關(guān)系，可得[16]

Lgls=2πresinψ(1)Sgl=2πre2(1-cosψ)(2)Sls=2πRp2(1-cosφ)(3)V=13πre3(2-3cosψ+cos3ψ)-13πRp3(2-3cosφ+cos3φ)(4)

式中，Lgls為氣-液-固三相接觸線的周長，m；Sgl、Sls分別為氣-液、液-固界面面積，m2；V為晶核體積，m3；ψ與φ的計(jì)算式分別為

ψ=arccosRpcosθ-res(5)φ=arccosRp-recosθs(6)

其中

s=Rp2+re2-2Rprecosθ(7)

由式(2)~式(4)可得[15]

dSgldr=2πRp2X1+1g-cosθg+X-cosθg3(1-Xcosθ)(8)dSlsdr=2πRpX2g1-X-cosθg2(9)dVdr=πRp22X2-(1-cos2θ)X2g5+21-XcosθgX2+1-cos2θg31-(cosθ-X)(1-Xcosθ)Xg2(10)

其中

X=Rpre(11)g=1+X2-2Xcosθ(12)

1.1.2晶核形成自由能蒸汽在不可溶顆粒表面異質(zhì)核化的晶核形成自由能ΔG可寫為[16]

ΔG=-kTlnS+σglSgl+(σls-σgs)Sls(13)

式中，k為Boltzmann常數(shù)，k=1.38×10-23J/K；T為溫度，K；S為蒸汽飽和度，S=pv/ps；pv為氣相蒸汽壓，Pa；ps為溫度為T時(shí)的飽和蒸氣壓，Pa；σgl、σls、σgs分別為氣-液、液-固、氣-固界面自由能，N/m。

根據(jù)Young方程，可得

σgs-σls=σglcosθ(14)

對(duì)于包含球形不可溶核和可溶無機(jī)鹽的顆粒，晶核生長過程伴隨著無機(jī)鹽的溶解。隨著晶核中可溶無機(jī)鹽的摩爾分?jǐn)?shù)的變化，晶核中水的化學(xué)勢和氣-液界面自由能也將發(fā)生變化，因而需要采用積分的方法計(jì)算晶核形成自由能。恒溫、恒壓下，晶核形成自由能ΔG與晶核半徑re的函數(shù)關(guān)系可寫為[17]

ΔG(re)=∫0reRgTlnawSdnwdrdr+∫0reσgldSgldrdr+∫0re(σls-σgs)dSlsdrdr(15)

式中，Rg為通用氣體常數(shù)，J/(mol·K)；aw為晶核中水的活度；nw為晶核中水的物質(zhì)的量；r為積分變量，m。

dnw/dr可寫為

dnwdr=dV/drdV/dnw(16)

根據(jù)晶核質(zhì)量和密度，其體積V又可寫為

V=nwMw+nsMsρs(17)

式中，Mw為水的摩爾質(zhì)量，kg/mol；ns為晶核中可溶無機(jī)鹽的物質(zhì)的量，mol；Ms為可溶無機(jī)鹽的摩爾質(zhì)量，kg/mol；ρs為晶核密度，kg/m3。

在晶核生長的起始階段，晶核為飽和溶液，其密度不隨晶核中水的摩爾分?jǐn)?shù)發(fā)生變化；隨著晶核的生長，晶核轉(zhuǎn)變?yōu)椴伙柡腿芤?，晶核密度隨晶核中水摩爾分?jǐn)?shù)的變化而變化?；谑?17)，當(dāng)可溶無機(jī)鹽未完全溶解時(shí)

dVdnw=Mw(Ss+100)100ρs0(18)

式中，Ss為可溶無機(jī)鹽的溶解度，g/100g；ρs0為飽和無機(jī)鹽溶液的密度，kg/m3。

當(dāng)可溶無機(jī)鹽完全溶解后

dVdnw=Mwρs(nw)-(nwMw-nw0Mw+ms0)dρs(nw)dnw1ρs(nw)(19)

式中，nw0為可溶無機(jī)鹽剛好完全溶解時(shí)晶核中水的物質(zhì)的量，mol；ms0為可溶無機(jī)鹽剛好完全溶解時(shí)晶核的質(zhì)量，kg；ρs(nw)表示ρs是nw的函數(shù)，即

ρs(nw)=(nw-nw0)Mw+ms0Vs0+(nw-nw0)Mw/ρw(20)

式中，ρw為水的密度，kg/m3；Vs0為可溶無機(jī)鹽剛好完全溶解時(shí)晶核的體積，m3。

1.1.3成核速率和臨界飽和度Fletcher經(jīng)典異質(zhì)核化理論模型[16]中給出的成核速率的表達(dá)式為

J=4πRp2Kexp-ΔG*kT(21)

式中，J為成核速率，s-1；ΔG*為臨界晶核形成自由能，即ΔG的極大值，J；K為成核動(dòng)力學(xué)系數(shù)，m-2·s-1。K的值取決于成核條件，具有一定的不確定度，數(shù)量級(jí)為1028~1031m-2·s-1，一般取K=1029m-2·s-1。

為避免采用不確定的成核動(dòng)力學(xué)系數(shù)K，依據(jù)晶核生長的兩種基本機(jī)制，即蒸汽分子在晶核表面的沉積機(jī)制與水分子在氣-液-固三相邊界處的表面擴(kuò)散機(jī)制，對(duì)K進(jìn)行計(jì)算，從而改進(jìn)的成核動(dòng)力學(xué)系數(shù)K可寫為[21]

K=Z(βdep+βdiff)n(22)

式中，Z為Zeldovich因子；βdep和βdiff分別為由于直接沉積機(jī)制和表面擴(kuò)散機(jī)制引起的水分子添加速率，s-1；n為顆粒表面吸附的水分子的數(shù)密度，m-2。

其中

Z=Vwπre*2σglkT×12+(1-Xcosθ)2-4Xcosθ-(cos2θ-3)X2(1-2Xcosθ+X2)3(23)βdep=pv2πmwkTSgl(24)βdiff=δυnexp-ΔGdiffkTLgls(25)n=pvυ2πmwkTexpΔGdeskT(26)

式中，Vw為單個(gè)水分子的體積，m3；mw為單個(gè)水分子的質(zhì)量，kg；δ為水分子的平均跳躍距離，m；υ為水分子的振動(dòng)頻率，Hz；ΔGdiff為單個(gè)水分子的表面擴(kuò)散能，J；ΔGdes為單個(gè)水分子的沉積能，J。

臨界飽和度是指成核速率J=1s-1時(shí)的蒸汽飽和度，是表征異質(zhì)核化性能的重要參數(shù)之一。經(jīng)典異質(zhì)核化理論認(rèn)為，當(dāng)蒸汽飽和度達(dá)到臨界飽和度時(shí)，蒸汽將能夠以顆粒表面為依托發(fā)生顯著的核化凝結(jié)現(xiàn)象。

1.1.4活度和氣-液界面自由能對(duì)于含有不可溶核和可溶無機(jī)鹽的顆粒，晶核形成自由能的計(jì)算需要建立在無機(jī)鹽溶液物化性質(zhì)參數(shù)——溶液中水的活度、氣-液界面自由能已知的基礎(chǔ)之上。晶核中可溶無機(jī)鹽組分和含量會(huì)對(duì)氣-液、液-固界面自由能帶來影響，但不影響氣-固界面自由能和接觸角。在氣-液界面自由能和接觸角已知的情況下，基于Young方程可獲得氣-固、液-固界面自由能之差，使得晶核形成自由能的計(jì)算成為可能。

單一組分無機(jī)鹽溶液中水的活度可寫為[23]

aw=1+∑i=14ai(100xi)i(27)

式中，xi為溶液中溶質(zhì)的質(zhì)量分?jǐn)?shù)；ai為常數(shù)，由溶質(zhì)確定。

單一組分無機(jī)鹽溶液的氣-液界面自由能可寫為[23]

σgl=σw+Bbi(28)

式中，σw為水的表面自由能，N/m；bi為溶液的濃度，mol/kg；B為常數(shù)，N·kg/(m·mol)。

氣溶膠無機(jī)-有機(jī)混合物官能團(tuán)活度系數(shù)（aerosolinorganic-organicmixturesfunctionalgroupsactivitycoefficient，AIOMFAC）模型[24]和擴(kuò)展的氣溶膠無機(jī)物模型（extendedaerosolinorganicsmodel，E-AIM）是計(jì)算氣溶膠系統(tǒng)中溶液物化性質(zhì)的重要模型[25]，前者給出了活度系數(shù)與溫度和溶液組分的函數(shù)關(guān)系，可以計(jì)算無機(jī)鹽溶液中水的活度；后者是氣溶膠系統(tǒng)中氣體/液體/固體分配的群體模型，可以計(jì)算無機(jī)鹽溶液的氣-液界面自由能和有機(jī)酸溶液中水的活度。本文利用AIOMFAC模型計(jì)算水的活度，利用E-AIM模型計(jì)算氣-液界面自由能。在此基礎(chǔ)上，通過對(duì)計(jì)算數(shù)據(jù)分別按照式(27)、式(28)進(jìn)行曲線擬合，獲得可溶無機(jī)鹽溶液中水的活度和溶液的氣-液界面自由能計(jì)算式中的系數(shù)值。圖2給出了常溫常壓下(T=298.15K、p=101325Pa)可溶無機(jī)鹽溶液中水的活度和溶液的氣-液界面自由能。圖中，離散的點(diǎn)表示模型計(jì)算結(jié)果，連續(xù)的線表示擬合曲線，所得的系數(shù)值見表1。

圖2

圖2常溫常壓下無機(jī)鹽溶液中水的活度和溶液的氣-液界面自由能

Fig.2Wateractivityandfreeenergyofgas-liquidinterfaceofinorganicsaltsolutionatroomtemperatureandpressure

表1常溫常壓下無機(jī)鹽溶液中水的活度和溶液的氣-液界面自由能計(jì)算系數(shù)

Table1Coefficientsforcalculatingwateractivityandsurfacetensionofinorganicsaltsolutionatroomtemperatureandpressure

參數(shù)符號(hào)NaClKClCaCl2活度a1-0.5286-0.3818-0.5507

a2-1.0298-0.77881.3166

a3-1.73300-17.9406

a40021.6473氣-液界面自由能B/(N·kg/(m·mol))1.6200×10-31.4923×10-33.4372×10-3新窗口打開|下載CSV

1.2計(jì)算方法和計(jì)算條件根據(jù)上述理論模型，建立蒸汽在顆粒表面異質(zhì)核化特性數(shù)值模擬平臺(tái)。在已知顆粒特性和氣相溫度、壓力、蒸汽飽和度條件下，對(duì)于不可溶顆粒，可根據(jù)式(13)計(jì)算任意晶核半徑下的晶核形成自由能；對(duì)于包含球形不可溶核和可溶無機(jī)鹽的顆粒，將晶核半徑由0以0.01nm的均勻間隔進(jìn)行離散，通過數(shù)值積分，即可獲得晶核形成自由能隨晶核半徑的變化關(guān)系?；诖说玫降木Ш诵纬勺杂赡艿臉O大值，即為臨界晶核形成自由能，其對(duì)應(yīng)的晶核半徑即為臨界晶核半徑，進(jìn)而可計(jì)算成核速率。成核臨界飽和度的計(jì)算，采用逐次逼近的方法，具體過程為：(1)將蒸汽飽和度S在1~3以均勻間隔0.01進(jìn)行離散，生成一系列飽和度值，記為Si=1+0.01i（i=0,1,2,…,100）；(2)計(jì)算每個(gè)Si條件下的Ji，若|Ji-1|0.01，則臨界成核飽和度Scr=Si；否則，搜尋使得Ji1、Ji+11的飽和度Si和Si+1，并在兩個(gè)飽和度之間以均勻間隔(Si+1-Si)/100進(jìn)行離散，生成一系列飽和度值；(3)重復(fù)步驟(2)直至獲得成核臨界飽和度。數(shù)值模擬采用的計(jì)算參數(shù)見表2，其中水分子特性參數(shù)來自文獻(xiàn)[26]。此外，受燃料特性和燃燒條件等的影響，生物質(zhì)燃燒過程排放細(xì)顆粒物中可溶無機(jī)鹽含量存在顯著差異；就單個(gè)細(xì)顆粒物而言，其可溶無機(jī)鹽質(zhì)量分?jǐn)?shù)fs可在0~100%之間取值。本文取fs=20%，事實(shí)上，由于對(duì)異質(zhì)核化特性起決定作用的臨界晶核尺寸僅為數(shù)納米，顆粒表面上可溶無機(jī)鹽含量fs≥0.01%時(shí)，晶核為飽和溶液，此條件下異質(zhì)核化特性不受可溶無機(jī)鹽含量的影響，使得本文結(jié)果適用于fs≥0.01%的情況。

表2數(shù)值模擬采用的計(jì)算參數(shù)

Table2Computationalparametersusedinnumericalsimulation

參數(shù)符號(hào)單位數(shù)值溫度TK298.15壓力pPa101325顆粒中不可溶核的密度ρpkg/m32400可溶無機(jī)鹽質(zhì)量分?jǐn)?shù)fs%20單個(gè)水分子的質(zhì)量mwkg2.99×10-26單個(gè)水分子的體積Vwm32.99×10-29水分子的平均跳躍距離δm3.2×10-10水分子的振動(dòng)頻率υHz1×1013單個(gè)水分子的表面擴(kuò)散能ΔGdiffJ2.9×10-21單個(gè)水分子的吸附能ΔGdesJ2.9×10-20新窗口打開|下載CSV

2結(jié)果與討論

課題組前期的研究中，將基于經(jīng)典異質(zhì)核化理論的蒸汽在不可溶球形顆粒表面的成核臨界飽和度數(shù)值模擬結(jié)果與文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)相結(jié)合，驗(yàn)證了模型和算法的可靠性[5]。本文將在前期研究的基礎(chǔ)上，考察包含不可溶球形核和可溶無機(jī)鹽的顆粒的核化特性，分析其與球形不可溶顆粒核化特性的差異，以及不同可溶組分對(duì)核化特性的影響規(guī)律。然而，由于目前尚缺少蒸汽在包含球形不可溶核和可溶無機(jī)鹽的顆粒表面異質(zhì)核化特性的實(shí)驗(yàn)報(bào)道，給相應(yīng)顆粒異質(zhì)核化數(shù)值模擬結(jié)果的實(shí)驗(yàn)驗(yàn)證帶來困難。盡管最近東南大學(xué)張軍教授團(tuán)隊(duì)[27]借助環(huán)境掃描電子顯微鏡實(shí)現(xiàn)了細(xì)顆粒物異質(zhì)核化過程的可視化，但研究對(duì)象均為不可溶顆粒。今后的研究中可制備包含球形不可溶核和可溶無機(jī)鹽的顆粒，并對(duì)其核化特性開展實(shí)驗(yàn)，為復(fù)雜顆粒異質(zhì)核化理論模型的驗(yàn)證提供數(shù)據(jù)支持，而本文模擬結(jié)果可以為實(shí)驗(yàn)研究中參數(shù)選擇提供先期指導(dǎo)。

2.1晶核形成自由能

圖3給出了S=1.2、Rp=0.1μm條件下得到的無量綱晶核形成自由能ΔG/(kT)隨晶核半徑re的變化關(guān)系。其中，圖3(a)~(c)分別為接觸角θ=30°、40°和50°下的結(jié)果。由圖可見，晶核形成自由能均呈現(xiàn)出隨著晶核半徑的增加而先增大后減小的趨勢，晶核形成自由能的極大值為臨界晶核形成自由能，其對(duì)應(yīng)的晶核半徑為臨界晶核半徑。在過飽和蒸汽與細(xì)顆粒物組成的系統(tǒng)中，水分子在顆粒表面聚集形成晶核，只有越過臨界晶核形成自由能、超過臨界晶核半徑的晶核才能持續(xù)長大，成為液滴。在Rp不變的情況下，隨著接觸角的增加，臨界晶核形成自由能增大，這是因?yàn)榻佑|角越大，顆粒的潤濕性能越差，越不利于顆粒充當(dāng)蒸汽相變的凝結(jié)核。由圖還可見，在接觸角一定的情況下，臨界晶核形成自由能的大小與顆粒表面是否含有可溶無機(jī)鹽以及可溶無機(jī)鹽的種類密切相關(guān)。其中，不可溶顆粒的臨界晶核形成自由能最大，含CaCl2顆粒的臨界晶核形成自由能最小。從物理化學(xué)機(jī)制上看，晶核形成自由能包括兩部分：一是體積項(xiàng)，即水分子由蒸汽態(tài)轉(zhuǎn)變?yōu)橐簯B(tài)，汽液兩相化學(xué)勢的差異，造成系統(tǒng)自由能降低；二是表面項(xiàng)，即液體晶核在顆粒表面生成，氣-液和液-固界面取代了原來的氣-固界面，使得系統(tǒng)自由能增加。與不可溶顆粒相比，含有可溶無機(jī)鹽的顆粒表面形成的晶核過程中由于可溶組分的溶解，氣液兩相化學(xué)勢之差與氣-液界面自由能均增大，從而體積項(xiàng)降低，表面項(xiàng)增加?？紤]到體積項(xiàng)與re的立方成正比，表面項(xiàng)與re的平方成正比，隨著re的增加，體積項(xiàng)的降低更為顯著，使得臨界晶核半徑與臨界晶核形成自由能均更小。此外，由于水的活度通過影響化學(xué)勢而影響體積項(xiàng)，其對(duì)晶核形成自由能的影響更為顯著。在模擬條件下，晶核均為飽和溶液，結(jié)合圖2可知，對(duì)于三種可溶無機(jī)鹽，飽和溶液中水的活度由小到大依次為CaCl2NaClKCl。綜上所述，晶核形成自由能按由小到大的順序排列為CaCl2NaClKCl。

圖3

圖3無量綱晶核形成自由能隨晶核半徑的變化關(guān)系

Fig.3Dimensionlessfreeenergyofembryoformationversusembryoradius

晶核生長過程中，可溶無機(jī)鹽的溶解不僅改變了晶核中水的活度和晶核的氣-液界面自由能，還由于溶解度和密度的不同而導(dǎo)致晶核體積隨晶核半徑的變化規(guī)律不同，使得含不同可溶無機(jī)鹽的顆粒表面上的晶核形成自由能隨晶核半徑的演變規(guī)律不同，從而引起臨界晶核形成自由能和臨界晶核半徑差異顯著。圖3(a)示出了不可溶顆粒的無量綱臨界晶核形成自由能ΔG*/(kT)和臨界晶核半徑re*。當(dāng)rere*時(shí)，由于系統(tǒng)具有自發(fā)地由非平衡態(tài)向平衡態(tài)轉(zhuǎn)變的特性，隨著ΔG→0，re→0，晶核將消失；而一旦晶核半徑達(dá)到re*，只要再捕集一個(gè)蒸汽分子，即可越過臨界晶核形成自由能，之后隨著ΔG→0，re單調(diào)增大，這表明，顆粒能夠充當(dāng)凝結(jié)核，隨后顆粒由于持續(xù)的凝結(jié)增長演變?yōu)榇罅降暮瑝m液滴。臨界晶核形成自由能和臨界晶核半徑的存在也表明蒸汽凝結(jié)過程中必須經(jīng)過異質(zhì)核化過程。一般說來，臨界晶核形成自由能和臨界晶核半徑越小，顆粒的核化性能越好，蒸汽越容易在顆粒表面發(fā)生異質(zhì)核化凝結(jié)。因此，四種顆粒相比，不可溶顆粒最難以發(fā)生異質(zhì)核化凝結(jié)，含KCl顆粒次之，含NaCl顆粒再次之，含CaCl2顆粒最容易發(fā)生異質(zhì)核化凝結(jié)。

2.2晶核生長機(jī)制

欲深入理解晶核生長過程，有必要探討兩種晶核生長機(jī)制對(duì)晶核生長的貢獻(xiàn)。由式(24)、式(25)可知，晶核生長的直接沉積機(jī)制和表面擴(kuò)散機(jī)制與晶核的幾何參數(shù)——?dú)?液界面面積及氣-液-固三相接觸線的周長密切相關(guān)，而晶核的幾何參數(shù)取決于不可溶顆?；虿豢扇芎说陌霃絉p與接觸角θ。在S=1.2條件下開展數(shù)值模擬，圖4給出了直接沉積和表面擴(kuò)散機(jī)制引起的水分子添加速率之比(βdiff/βdep)的模擬結(jié)果。其中，圖4(a)為θ=40°時(shí)βdiff/βdep隨Rp及晶核半徑re的變化情況，圖4(b)為θ=40°時(shí)臨界晶核半徑下βdiff/βdep隨Rp的變化情況，圖4(c)為Rp=0.1μm時(shí)臨界晶核半徑下βdiff/βdep隨θ的變化情況。

圖4

圖4直接沉積和表面擴(kuò)散機(jī)制引起的分子添加速率之比

Fig.4Ratioofmoleculeadditionratesduetodirectdepositionandsurfacediffusion

由圖4(a)可以看出，在模擬條件下，表面擴(kuò)散機(jī)制在晶核生長過程中起著決定性作用，對(duì)于給定的Rp，隨著晶核的增長，蒸汽分子在顆粒表面的直接沉積機(jī)制的貢獻(xiàn)趨于增大，然而對(duì)于re10nm的晶核，水分子的表面擴(kuò)散對(duì)晶核生長的貢獻(xiàn)仍顯著超過蒸汽分子的直接沉積的貢獻(xiàn)。此外，由式(24)、式(25)分析可得，兩種生長機(jī)制引起的分子添加速率均不受顆粒中可溶無機(jī)鹽的影響，這表明圖4(a)中的結(jié)果可以表征本文涉及的不可溶球形顆粒和3種包含可溶無機(jī)鹽和球形不可溶核的顆粒。

然而，顆粒組分的不同會(huì)導(dǎo)致臨界晶核半徑re*不同，進(jìn)而影響到表面擴(kuò)散和直接沉積兩種成核機(jī)制對(duì)臨界晶核增長的貢獻(xiàn)。由圖4(b)、(c)可以看出，顆粒含有可溶無機(jī)鹽時(shí)，臨界晶核半徑下βdiff/βdep顯著增加，即表面擴(kuò)散機(jī)制對(duì)臨界晶核生長的貢獻(xiàn)更大，特別是含有CaCl2時(shí)，βdiff/βdep與不可溶顆粒相比提高了5倍以上。這是因?yàn)?，?dāng)晶核中含有可溶無機(jī)鹽時(shí)，臨界晶核半徑大幅減?。▓D3），表面擴(kuò)散機(jī)制引起的分子添加速率βdiff取決于三相接觸線的周長，其與晶核半徑成正比，表面沉積機(jī)制引起的分子添加速率βdep取決于氣-液界面面積，其與晶核半徑的平方成正比，因此βdiff/βdep與晶核半徑成反比，從而臨界晶核半徑越小，臨界晶核半徑下βdiff/βdep就越大。臨界晶核條件下，對(duì)于不同組分顆粒，βdiff/βdep均隨Rp呈現(xiàn)出先略有增加，而后保持不變的變化規(guī)律，同時(shí)βdiff/βdep隨θ的增大而單調(diào)下降，這同樣可以通過臨界晶核尺寸的變化加以解釋。郭陽等[15]對(duì)氨法脫硫系統(tǒng)排放細(xì)顆粒物異質(zhì)核化特性的研究中發(fā)現(xiàn)，含有可溶組分(NH4)2SO4與(NH4)HSO4時(shí)，臨界晶核半徑隨著不可溶核半徑的增大而減小，但減小量迅速下降；而已有異質(zhì)核化數(shù)值模擬研究表明，接觸角越大，臨界晶核半徑越大[5,19]。因此，臨界晶核條件下，βdiff/βdep呈現(xiàn)出如圖4的變化規(guī)律。

2.3成核臨界飽和度

圖5給出了成核臨界飽和度隨不可溶顆粒半徑或不可溶核半徑以及接觸角的變化關(guān)系。成核臨界飽和度與成核速率密切相關(guān)，在成核臨界飽和度下，單位時(shí)間（1s）內(nèi)顆粒表面形成一個(gè)臨界晶核，即成核速率為1s-1；并且在氣相參數(shù)以及Rp相同時(shí)，成核速率越高，意味著顆粒異質(zhì)核化性能更好，顆粒的成核臨界飽和度往往更低。圖5(a)的數(shù)值模擬中保持θ=40°不變，圖5(b)的數(shù)值模擬中采用Rp=0.1μm。由圖5(a)可知，不可溶顆粒與3種含有可溶無機(jī)鹽的顆粒，成核臨界飽和度均隨Rp的增大而減小，但成核臨界飽和度在數(shù)值上差異較大。由圖5(a)還可知，在Rp=0.1μm條件下，對(duì)于不可溶顆粒，Scr=1.388；對(duì)于含有KCl和NaCl的顆粒，Scr分別為1.339、1.209；對(duì)于含CaCl2顆粒，Scr=1.004。這表明不可溶顆粒表面附著有可溶無機(jī)鹽時(shí)，成核臨界飽和度下降，特別是附著有CaCl2時(shí)，蒸汽只要稍微超過飽和狀態(tài)，就能夠在顆粒表面發(fā)生異質(zhì)核化。由圖5(b)可知，接觸角對(duì)成核臨界飽和度有重要影響，接觸角越大，成核臨界飽和度越大，這是由于接觸角大，意味著顆粒潤濕性差，水分子更難在顆粒表面聚集，造成顆粒難以充當(dāng)凝結(jié)核，表現(xiàn)為成核速率降低，成核臨界飽和度增高，即在更高的飽和度條件，蒸汽才能夠在顆粒表面發(fā)生異質(zhì)核化凝結(jié)。由圖還可以看出，當(dāng)接觸角較小（30°~50°）時(shí)，不可溶顆粒的臨界飽和度最大；當(dāng)接觸角較大（≥50°）時(shí)，含KCl和NaCl顆粒的成核臨界飽和度逐漸超過不可溶顆粒的成核臨界飽和度。例如，當(dāng)θ=30°時(shí)，不可溶顆粒與含KCl、NaCl、CaCl2顆粒的成核臨界飽和度分別為1.218、1.156、1.039、1.002，然而當(dāng)θ=70°，4種組分顆粒的成核臨界飽和度依次為2.257、2.519、2.315、1.165。含有可溶無機(jī)鹽的顆粒的成核臨界飽和度高于不可溶顆粒的臨界飽和度的原因可解釋為：接觸角增大，引起臨界晶核形成自由能和臨界晶核半徑增大，受可溶無機(jī)鹽溶解、晶核中水的活度、晶核氣-液界面自由能的影響，含可溶無機(jī)鹽的顆粒與不可溶顆粒之間、含不同可溶無機(jī)鹽的顆粒之間的臨界晶核形成自由能和臨界晶核半徑的增大程度不同，導(dǎo)致不同種類顆粒的成核速率有不同程度的減小，從而顆粒的成核臨界飽和度有不同程度的增加，產(chǎn)生了圖5(b)所示的結(jié)果。圖5中結(jié)果還表明，當(dāng)接觸角較大時(shí)，除含CaCl2顆粒外，其余3種顆粒的成核臨界飽和度均較高，在實(shí)際工程應(yīng)用中，可通過添加潤濕劑以降低接觸角，改善顆粒的潤濕性，從而在較低飽和度條件下實(shí)現(xiàn)顆粒的異質(zhì)核化凝結(jié)增長，以降低能耗和水耗。此外，生物質(zhì)燃燒過程排放細(xì)顆粒物有相當(dāng)一部分是分形團(tuán)聚體顆粒，已有數(shù)值模擬和實(shí)驗(yàn)研究均表明，分形結(jié)構(gòu)降低了顆粒接觸角，從而有利于顆粒的異質(zhì)核化，與球形顆粒相比，分形顆粒成核臨界飽和度有所下降[18-19]。

圖5

圖5成核臨界飽和度

Fig.5Criticalsaturationfornucleation

3結(jié)論

建立了改進(jìn)的蒸汽在不可溶球形顆粒與包含球形不可溶核