求解顆粒凝并與反彈的DSMC方法及應(yīng)用

吳志豪， 凡鳳仙， 胡曉紅

(上海理工大學(xué)能源與動(dòng)力工程學(xué)院， 上海200093)

氣固兩相流廣泛存在于能源、環(huán)境、化工等領(lǐng)域,如煤粉的氣力輸送[1]、 固體燃料顆粒在燃燒室內(nèi)的流動(dòng)與燃燒[2]、 細(xì)顆粒物在除塵器內(nèi)的流動(dòng)與分離[3]、 大氣顆粒物的擴(kuò)散和演化[4]、 流化床內(nèi)顆粒物料的流化[5-6]等。近年來(lái)，隨著計(jì)算機(jī)硬件的更新?lián)Q代和數(shù)值模擬技術(shù)的迅速發(fā)展，計(jì)算機(jī)數(shù)值模擬已成為研究氣固兩相流過(guò)程和機(jī)理的重要手段。建立可靠的描述顆粒運(yùn)動(dòng)、 碰撞、 凝并及反彈等物理過(guò)程的顆粒動(dòng)力學(xué)模型，對(duì)氣固兩相流的高精度模擬顯得非常重要和迫切。離散元方法(discrete element method，DEM)是顆粒動(dòng)力學(xué)行為建模和模擬的最直接的方法，該方法基于牛頓第二定律，根據(jù)顆粒與顆粒之間的相互作用力和顆粒與流體之間的耦合作用來(lái)求解每個(gè)顆粒的速度與運(yùn)動(dòng)軌跡，通過(guò)跟蹤所有顆粒的運(yùn)動(dòng)軌跡來(lái)判斷顆粒間碰撞的發(fā)生。 盡管DEM對(duì)顆粒動(dòng)力學(xué)過(guò)程的描述最為接近實(shí)際， 但其在處理涉及大量顆粒的氣固兩相流時(shí)， 尤其是當(dāng)固相為粉體顆粒時(shí)，計(jì)算成本極高。直接模擬蒙特卡洛(direct simulation Monte Carlo, DSMC)方法是對(duì)氣固兩相流中顆粒動(dòng)力學(xué)過(guò)程進(jìn)行建模和模擬的另一種方法。DSMC方法在硬球模型的基礎(chǔ)上結(jié)合了判斷顆粒碰撞的蒙特卡洛方法，在稀相氣固兩相流數(shù)值模擬中得到廣泛應(yīng)用。該方法立足于氣體分子運(yùn)動(dòng)論，采用取樣顆粒代替真實(shí)顆粒，為每個(gè)取樣顆粒賦予一定的權(quán)重，使其代表一定數(shù)目的真實(shí)顆粒，跟蹤取樣顆粒在流場(chǎng)中的受力與運(yùn)動(dòng)，通過(guò)概率判斷顆粒碰撞的發(fā)生，從而有效地節(jié)省了計(jì)算時(shí)間與計(jì)算機(jī)內(nèi)存，適宜于對(duì)大量顆粒構(gòu)成的稀疏氣固兩相流中顆粒動(dòng)力學(xué)模擬。

DSMC方法在氣固兩相流模擬的實(shí)際應(yīng)用中往往將顆粒碰撞的結(jié)果處理為單一事件，即所有碰撞都導(dǎo)致顆粒反彈或都導(dǎo)致顆粒凝并。例如，在顆粒碰撞導(dǎo)致顆粒反彈的假設(shè)下，Tsuji等[7]對(duì)循環(huán)流化床上升管內(nèi)的團(tuán)簇分布進(jìn)行了數(shù)值預(yù)測(cè)，F(xiàn)an等[8]研究了駐波聲場(chǎng)中顆粒的碰撞頻率，F(xiàn)ede等[9]探究了湍流中的顆粒碰撞。在相關(guān)顆粒碰撞-反彈的研究中，往往對(duì)所有取樣顆粒采用相同的數(shù)目權(quán)重，這種等權(quán)重算法便于應(yīng)用，但當(dāng)研究對(duì)象為寬分布多分散顆粒體系或?yàn)楹猩倭刻囟ǔ煞值幕旌项w粒體系時(shí)，等權(quán)重算法將可能導(dǎo)致代表少量顆粒的取樣顆粒數(shù)目不足，甚至缺失，從而造成統(tǒng)計(jì)噪聲增加，計(jì)算精度降低。為了解決這一問(wèn)題，He等[10]研究出一種處理顆粒碰撞-反彈行為的異權(quán)重DSMC方法。DSMC方法在黏性顆粒凝并過(guò)程模擬上得到廣泛應(yīng)用，目前已發(fā)展出常體積蒙特卡洛(constant-volume Monte Carlo，CVMC)方法[11]、 常數(shù)目蒙特卡洛(constant-number Monte Carlo，CNMC)方法[12]、 多重蒙特卡洛(multi-Monte Carlo，MMC)方法[13]。其中，CVMC方法和CNMC方法均采用等權(quán)重算法。CVMC方法保持計(jì)算區(qū)域不變，隨著顆粒凝并的進(jìn)行，取樣顆粒數(shù)目不斷減少，使得計(jì)算精度不斷降低，因此，為了保證計(jì)算精度，往往采用大量取樣顆粒，導(dǎo)致計(jì)算成本很高。CNMC方法通過(guò)增大計(jì)算區(qū)域體積以保持取樣顆粒數(shù)目不變，但該方法難以考慮顆粒的空間擴(kuò)散和計(jì)算區(qū)域的邊界條件，使其應(yīng)用范圍受到限制。針對(duì)這2種方法的不足，Zhao等[13]提出了一種異權(quán)重算法，即MMC方法，該方法在不改變計(jì)算區(qū)域體積的同時(shí)，通過(guò)調(diào)整取樣顆粒的數(shù)目權(quán)重，以保持計(jì)算過(guò)程中取樣顆粒數(shù)目不變。MMC方法應(yīng)用于亞微米顆粒的布朗凝并和聲凝并建模和模擬，展現(xiàn)出對(duì)亞微米顆粒凝并行為優(yōu)良的數(shù)值預(yù)測(cè)能力[13-16]。

盡管DSMC方法在氣固兩相流數(shù)值模擬中發(fā)揮著重要作用，并已得到廣泛應(yīng)用，然而該方法在對(duì)碰撞結(jié)果的處理上存在一定的局限性。分析顆粒碰撞的物理過(guò)程可以發(fā)現(xiàn)，大粒徑和高碰撞速度將會(huì)導(dǎo)致顆粒反彈，而小粒徑和低碰撞速度條件下顆粒碰撞的結(jié)果則傾向于凝并[17]。當(dāng)對(duì)粒徑分布寬廣、 氣相場(chǎng)復(fù)雜的氣固兩相流進(jìn)行建模與模擬時(shí)，如果仍對(duì)顆粒碰撞的后果進(jìn)行單一化處理(全部為反彈或凝并)，將導(dǎo)致數(shù)值預(yù)測(cè)結(jié)果與實(shí)際有較大的偏差。為解決這一問(wèn)題，本文中對(duì)已有研究進(jìn)行繼承和發(fā)展，構(gòu)建了一種能夠求解顆粒碰撞后凝并與反彈共存的問(wèn)題的DSMC方法，并基于該方法對(duì)寬粒徑雙峰分布的燃煤飛灰聲凝并效果進(jìn)行預(yù)測(cè)，獲得了與實(shí)驗(yàn)吻合良好的結(jié)果，驗(yàn)證了模型的可靠性。

1 顆粒凝并與反彈共存的DSMC方法

1.1 顆粒碰撞概率

DSMC方法認(rèn)為顆粒碰撞為隨機(jī)的二元碰撞，并利用概率的方法來(lái)判斷顆粒碰撞的發(fā)生。通常對(duì)計(jì)算區(qū)域劃分網(wǎng)格，認(rèn)為只有處于同一網(wǎng)格內(nèi)的顆粒才可能發(fā)生碰撞，則取樣顆粒i與取樣顆粒j所代表的真實(shí)顆粒發(fā)生碰撞的概率[16]可寫(xiě)為

Pij=βijwjΔtcol/Vci，

(1)

式中：Pij為碰撞概率；βij為碰撞核函數(shù)， m3/s；wj為取樣顆粒j的數(shù)目權(quán)重； Δtcol為求解顆粒間碰撞的時(shí)間間隔， s；Vci為取樣顆粒i所在網(wǎng)格的體積， m3。

取樣顆粒i和所有顆粒的碰撞概率Pi可寫(xiě)為

(2)

式中，N為取樣顆粒i所在網(wǎng)格內(nèi)取樣顆?？倲?shù)目。

1.2 顆粒碰撞的判斷

采用修正的南布方法[7]判斷顆粒碰撞事件的發(fā)生。在Pij<1/N的前提下，利用計(jì)算機(jī)產(chǎn)生的[0, 1)范圍均勻分布的隨機(jī)數(shù)R，選擇同一網(wǎng)格內(nèi)的取樣顆粒j作為候選被碰顆粒：

j=int[RN]+1，

(3)

式中，int[RN]表示RN的整數(shù)部分。

如果滿足關(guān)系式

R>j/N-Pij，

(4)

則認(rèn)為取樣顆粒i與j之間發(fā)生碰撞。

1.3 顆粒碰撞后果的確定

顆粒碰撞后果由顆粒碰撞后的動(dòng)能和顆粒間的黏附力共同決定。對(duì)于不帶電的固體干顆粒，范德華力為顆粒間唯一的黏附力?；谀芰科胶怅P(guān)系式，可得臨界碰撞速度vc的表達(dá)式[17]為

(5)

式中：A為Hamaker常數(shù)， J；e為恢復(fù)系數(shù)；z0為最小接觸距離， m；pc為極限接觸壓力， Pa；ρ為顆粒密度， kg/m3；ds和dl分別為較小和較大顆粒的直徑，m。

當(dāng)碰撞顆粒之間的法向相對(duì)速度不大于臨界速度，即滿足式(6)時(shí)，碰撞后果為凝并；反之為反彈。

(6)

式中：v為顆粒速度， m/s；n為由顆粒i指向顆粒j的法向單位向量； 上標(biāo)0表示碰前值。

1.4 顆粒凝并的處理

若取樣顆粒i和j發(fā)生凝并，則意味著它們代表的真實(shí)顆粒發(fā)生凝并。分析顆粒凝并的物理過(guò)程可知，真實(shí)顆粒發(fā)生碰撞-凝并的次數(shù)為min(wi,wj)，形成了min(wi,wj)個(gè)顆粒團(tuán)聚體，同時(shí)數(shù)目權(quán)重較大的取樣顆粒所代表的真實(shí)顆粒中有|wi-wj|個(gè)未參與碰撞-凝并過(guò)程。根據(jù)質(zhì)量和動(dòng)量守恒定律，可得：

(7)

(8)

式中：m為顆粒質(zhì)量， kg；Vagg為團(tuán)聚體體積， m3；vagg為團(tuán)聚體的速度。

對(duì)顆粒凝并的處理方案為保持取樣顆粒數(shù)目不變，調(diào)整取樣顆粒的數(shù)目權(quán)重、體積和速度以體現(xiàn)碰撞-凝并的發(fā)生[13-15]，具體如下：

當(dāng)wi=wj=2l(l=1, 2, 3, …)時(shí)，

(9)

當(dāng)wi=wj=2l-1(l=1, 2, 3, …)時(shí)，

(10)

當(dāng)wi>wj時(shí)，

(11)

當(dāng)wi

(12)

式中,上標(biāo)new表示碰撞后。

1.5 顆粒反彈的處理

對(duì)于每次真實(shí)顆粒碰撞-反彈過(guò)程，顆粒反彈速度可根據(jù)動(dòng)量方程計(jì)算。若不考慮顆粒轉(zhuǎn)動(dòng)，并認(rèn)為顆?；七^(guò)程所受的摩擦力遵守庫(kù)侖摩擦定律，且顆粒一旦停止滑移，不會(huì)再次發(fā)生滑移，顆粒反彈速度[7]為：

(13)

(14)

(15)

(16)

(17)

(18)

分析碰撞-反彈過(guò)程可知，對(duì)于等權(quán)重的取樣顆粒的碰撞-反彈，僅需按式(13)—(16)更新顆粒速度，以表征顆粒反彈的結(jié)果。對(duì)于不等權(quán)重的兩取樣顆粒的碰撞-反彈，采用“分離-合并”方法處理碰撞-反彈過(guò)程[10]。 以wi>wj時(shí)為例，真實(shí)顆粒碰撞-反彈次數(shù)為wj，即取樣顆粒i所代表的wj個(gè)真實(shí)顆粒與取樣顆粒j所代表的所有真實(shí)顆粒發(fā)生碰撞并產(chǎn)生反彈；而取樣顆粒i所代表的剩余的(wi-wj)個(gè)真實(shí)顆粒未參與到碰撞-反彈過(guò)程中。基于此，首先將取樣顆粒i分離為子取樣顆粒i*和i**，其數(shù)目權(quán)重分別為wj和(wi-wj)；接著，按照等權(quán)重取樣顆粒碰撞-反彈處理方法，獲得i*與j的反彈速度；最后，將碰撞-反彈后的子取樣顆粒i*和不參與碰撞-反彈過(guò)程的子取樣顆粒i**進(jìn)行合并，并將合并后的顆粒重新存儲(chǔ)為取樣顆粒i，其數(shù)目權(quán)重恢復(fù)為wi，其速度可由動(dòng)量守恒定律計(jì)算：

(19)

式中： 上標(biāo)mer表示合并后；下標(biāo)i*和i**分別表示顆粒i*和i**。

在合并過(guò)程中，不可避免地產(chǎn)生了能量損失，即

(20)

式中, ΔEi為能量損失。

為保持合并前后系統(tǒng)的能量不變，將能量損失補(bǔ)償?shù)絽⑴c后續(xù)碰撞的某一取樣顆粒s上，根據(jù)能量平衡方程，補(bǔ)償后顆粒s的速度為

(21)

式中, 上標(biāo)com表示補(bǔ)償后。

為在能量補(bǔ)償對(duì)系統(tǒng)擾動(dòng)盡量小的前提下能夠快速搜尋到顆粒s，采用的搜索規(guī)則為

(22)

2 顆粒凝并與反彈共存的DSMC方法在聲凝并中的應(yīng)用

聲凝并是指利用高強(qiáng)度外加聲場(chǎng)作用促使顆粒發(fā)生相互作用，進(jìn)而碰撞接觸，凝并為粒徑較大顆粒的現(xiàn)象。聲凝并作為富有工業(yè)應(yīng)用潛力的顆粒預(yù)處理技術(shù)，為燃燒源、工業(yè)源細(xì)顆粒物的高效脫除提供了一種重要技術(shù)途徑。聲凝并過(guò)程中涉及的顆粒既有亞微米尺度黏附顆粒，又有微米尺度慣性顆粒，適用于本文中建立的顆粒凝并與反彈共存的DSMC方法的驗(yàn)證。

2.1 取樣顆粒運(yùn)動(dòng)方程

在聲場(chǎng)中，顆粒除發(fā)生碰撞外，還將受到曳力、重力、浮力作用發(fā)生運(yùn)動(dòng)，顆粒i的運(yùn)動(dòng)方程[18-19]可寫(xiě)為：

(23)

式中：μg為氣體的動(dòng)力黏度， Pa·s；u為氣體速度， m/s；ρg為氣體密度， kg/m3；g為重力加速度， m/s2；Cine為慣性修正系數(shù)；CCun為肯寧漢滑移修正系數(shù)。

(24)

CCun=1+Kn[1.257+0.400exp(-1.100/Kn)]，

(25)

式中:Kn為克努森數(shù);Re為雷諾數(shù)。 其表達(dá)式分別為

(26)

(27)

式中,λg為氣體分子平均自由程，m。

水平平面行波聲場(chǎng)的波動(dòng)方程為

ux=uampsin[ω(t-x/c)]，

(28)

式中：ux為聲波引起的氣體介質(zhì)振動(dòng)速度，m/s；uamp為聲波速度振幅，m/s；ω為聲波的角頻率， rad/s；t為時(shí)間， s；x為距聲源的距離， m；c為聲速，m/s。

通常采用聲壓級(jí)L(單位dB)描述聲場(chǎng)，聲波速度振幅可根據(jù)聲壓級(jí)計(jì)算，

(29)

式中,Pr為參考聲壓，Pr=2×10-5Pa。

2.2 碰撞核函數(shù)

碰撞核函數(shù)是碰撞概率求解中的一個(gè)重要參數(shù)(見(jiàn)式(1))， 用于描述顆粒相互作用機(jī)理對(duì)顆粒碰撞的作用效果， 定義為單位時(shí)間內(nèi)顆粒因某種機(jī)理而能夠與其它顆粒發(fā)生碰撞的體積。 聲凝并中的碰撞核函數(shù)包括同向相互作用、 重力沉降作用、 聲尾流效應(yīng)、 互輻射壓力效應(yīng)以及布朗擴(kuò)散效應(yīng)核函數(shù)。

顆粒i和j之間的同向相互作用和重力沉降作用碰撞核函數(shù)βij，Orth-Grav可基于顆粒i和j運(yùn)動(dòng)方程的解vi和vj[16]給出，

(30)

聲尾流效應(yīng)碰撞核函數(shù)βij, AWE[20]可表示為:

(31)

式中：r為顆粒i與j的中心距， m；li和lj分別為顆粒i和j的滑移系數(shù)。

li的表達(dá)式為

(32)

其中，

(33)

(34)

(35)

(36)

互輻射壓力效應(yīng)碰撞核函數(shù)βij,MRPE的表達(dá)式[21]為

(37)

式中：gijr和gjir為顆粒i與顆粒j之間的流體力學(xué)相互作用函數(shù)，具體見(jiàn)文獻(xiàn)[21]。

布朗擴(kuò)散效應(yīng)碰撞核函數(shù)βij,Brow可表示[22]為

(38)

式中：T為溫度， K；kB為玻爾茲曼常數(shù)，kB= 1.38×10-23J/K。

總的碰撞核函數(shù)βij可寫(xiě)為

βij=βij,Orth-Grav+βij,AWE+βij,MRPE+βij,Brow

。

(39)

2.3 聲凝并數(shù)值模擬條件和流程

2.3.1 時(shí)間步長(zhǎng)

求解顆粒運(yùn)動(dòng)方程的時(shí)間步長(zhǎng)Δt[18]應(yīng)滿足

(40)

(41)

在求解顆粒碰撞的時(shí)間步長(zhǎng)Δtcol的選擇上，考慮到修正的南布方法的要求Pij<1/N，則有

(42)

基于上述選擇原則，為加速計(jì)算過(guò)程，本文中采用Δtcol=100Δt，即顆粒運(yùn)動(dòng)每求解100步，顆粒碰撞進(jìn)行搜索和求解一遍。

2.3.2 初始條件和邊界條件

沿聲波傳播方向取一個(gè)波長(zhǎng)λ范圍即λ×0.01 m×0.01 m的三維區(qū)域?yàn)橛?jì)算區(qū)域，并將其沿x向劃分成160個(gè)網(wǎng)格。依據(jù)每個(gè)粒徑區(qū)間內(nèi)顆粒的數(shù)目濃度，對(duì)取樣顆粒的初始數(shù)目權(quán)重進(jìn)行分配。將取樣顆粒隨機(jī)均勻地添加到計(jì)算區(qū)域，并為每個(gè)顆粒分配一個(gè)隨機(jī)的微小初始速度。采用循環(huán)邊界條件處理流出流進(jìn)計(jì)算區(qū)域的顆粒。

2.3.3 計(jì)算流程

數(shù)值計(jì)算流程如下：

1)初始化取樣顆粒場(chǎng)?；诔跏剂椒植籍a(chǎn)生取樣顆粒，為每個(gè)顆粒賦予一定的數(shù)目權(quán)重、位置和速度。

2)求解取樣顆粒的運(yùn)動(dòng)。采用四階龍格-庫(kù)塔算法，基于2.1節(jié)給出的顆粒運(yùn)動(dòng)方程，計(jì)算每一個(gè)取樣顆粒經(jīng)過(guò)一個(gè)時(shí)間步長(zhǎng)Δt后的速度和位置，并基于循環(huán)邊界條件對(duì)流出計(jì)算區(qū)域的顆粒進(jìn)行處理。該步驟重復(fù)100次。

3)判斷顆粒的碰撞并對(duì)碰撞后果進(jìn)行處理。按照第1部分結(jié)合第2.2節(jié)的碰撞核函數(shù)，對(duì)顆粒的碰撞及碰撞后的凝并與反彈進(jìn)行處理。

4)重復(fù)第2)、3)步的計(jì)算，直到達(dá)到設(shè)定的聲凝并時(shí)間。

2.4 數(shù)值模擬與實(shí)驗(yàn)結(jié)果的對(duì)比驗(yàn)證

為驗(yàn)證模型的準(zhǔn)確性，采用和Zhang等[23]的實(shí)驗(yàn)報(bào)道一致的燃煤飛灰顆粒粒徑分布和物性參數(shù)、氣相參數(shù)、聲場(chǎng)條件和聲波作用時(shí)間，開(kāi)展數(shù)值模擬。數(shù)值模擬參數(shù)見(jiàn)表1。

表1 數(shù)值模擬參數(shù)

圖1給出了數(shù)值模擬得到的聲凝并后顆粒粒徑分布與實(shí)驗(yàn)結(jié)果的對(duì)比情況。其中，數(shù)值模擬結(jié)果包括2種情況：1)將顆粒碰撞結(jié)果全部視為凝并；2)將顆粒碰撞結(jié)果視為凝并和反彈共存的情況。由圖可見(jiàn)，將顆粒碰撞結(jié)果全部視為凝并時(shí)，數(shù)值模擬得到的聲凝并后顆粒粒徑分布與實(shí)驗(yàn)差異顯著；當(dāng)利用臨界碰撞速度將顆粒碰撞結(jié)果考慮為凝并和反彈共存時(shí)，所得結(jié)果和實(shí)驗(yàn)吻合良好。這表明，對(duì)于粒徑涉及亞微米至數(shù)微米的飛灰顆粒，聲凝并過(guò)程中顆粒碰撞后的反彈不能忽略，此時(shí)若將顆粒碰撞后果全部視為凝并，將導(dǎo)致數(shù)值模擬的預(yù)測(cè)結(jié)果優(yōu)于實(shí)際情況，從而很大程度上高估了聲凝并的效果。

a)粒徑分布b)模擬結(jié)果的局部放大圖1 數(shù)值模擬得到的聲凝并后顆粒粒徑分布與實(shí)驗(yàn)結(jié)果的對(duì)比Fig.1 Comparison of particle size distribution after acoustic agglomeration obtained from numerical simulation with experimental data

由圖1中數(shù)值模擬結(jié)果可知，聲凝并后粒徑在10 μm以下顆粒的數(shù)量濃度降低，同時(shí)出現(xiàn)了粒徑在10 μm以上的顆粒，這些新顆粒是由聲凝并形成的團(tuán)聚體。在外加聲場(chǎng)作用下，粒徑在2.5 μm以下的細(xì)顆粒因慣性很小而幾乎被聲波完全夾帶，使得細(xì)顆粒間相對(duì)運(yùn)動(dòng)很弱，碰撞概率很低，凝并效果很差；粒徑在2.5 μm以上的粗顆粒，慣性逐漸增強(qiáng)，而越來(lái)越難以被夾帶，使得粗顆粒間相對(duì)運(yùn)動(dòng)較明顯，碰撞概率大，凝并效果也較好[16,19]。同時(shí)，細(xì)顆粒和粗顆粒之間夾帶程度的不同導(dǎo)致的顯著的顆粒間相對(duì)運(yùn)動(dòng)，也使得細(xì)顆粒易于與粗顆粒發(fā)生碰撞，從而細(xì)顆粒數(shù)目減小，因此，聲凝并過(guò)程中生成的10 μm以上的大顆粒往往是由捕獲了細(xì)顆粒的粗顆粒之間發(fā)生凝并或粗顆粒間凝并后再捕獲細(xì)顆粒形成的團(tuán)聚體。在今后的研究中，可進(jìn)一步挖掘DSMC方法的優(yōu)勢(shì)，在數(shù)值模擬過(guò)程中追蹤顆粒團(tuán)聚體的構(gòu)成信息，以直觀展示顆粒團(tuán)聚體的構(gòu)成和演變特性。此外，需要說(shuō)明的是：實(shí)驗(yàn)結(jié)果未能給出10 μm以上顆粒的粒徑分布，這是由于實(shí)驗(yàn)中采用靜電低壓沖擊器(ELPI)測(cè)量聲凝并前后顆粒的粒徑分布，而ELPI的測(cè)量范圍是0.007～10 μm[23]。

為了進(jìn)一步說(shuō)明數(shù)值模擬結(jié)果的可靠性，圖2給出了聲凝并前后顆粒的累積質(zhì)量濃度分布。由圖可見(jiàn)，數(shù)值模擬前后顆粒的總質(zhì)量濃度相等，這表明數(shù)值模擬中顆粒質(zhì)量是守恒的。綜合圖1和圖2可知，本文中建立的求解顆粒凝并與反彈的DSMC方法用于聲凝并過(guò)程模擬能夠準(zhǔn)確預(yù)測(cè)出顆粒粒徑分布的演變。該模型為準(zhǔn)確、高效模擬包含寬粒徑分布顆粒的稀疏氣固兩相流行為提供了重要基礎(chǔ)。

圖2 聲凝并前后顆粒的累積質(zhì)量濃度分布Fig.2 Cumulative particle mass concentration distribution before and after acoustic agglomeration

3 結(jié)論

DSMC方法作為一種高效、經(jīng)濟(jì)的數(shù)值模擬方法，在稀疏氣固兩相流數(shù)值模擬研究中得到了廣泛應(yīng)用，然而，已有相關(guān)研究往往認(rèn)為顆粒碰撞的結(jié)果全部是顆粒反彈，或者全部是顆粒凝并，使得DSMC方法在預(yù)測(cè)復(fù)雜氣相場(chǎng)或?qū)捔椒植碱w粒條件下的氣固兩相流時(shí)，具有明顯的局限性。為解決已有研究中存在的問(wèn)題，本文中做了以下工作。

1)基于顆粒碰撞后發(fā)生凝并和反彈的臨界速度，以及異權(quán)重DSMC方法的基本思想，發(fā)展了顆粒凝并與反彈共存的DSMC方法。

2)基于所發(fā)展的方法，建立了顆粒聲凝并模型，提出了相應(yīng)的數(shù)值模擬方法。

3)對(duì)聲凝并效果進(jìn)行了數(shù)值模擬，發(fā)現(xiàn)對(duì)顆粒碰撞結(jié)果進(jìn)行單一凝并處理時(shí)，數(shù)值模擬得到的聲凝并后顆粒粒徑分布與實(shí)驗(yàn)結(jié)果存在顯著差異；而將顆粒凝并后果處理為凝并與反彈共存時(shí)，則數(shù)值模擬與實(shí)驗(yàn)結(jié)果吻合良好。