豎直降落式反應(yīng)器中顆粒流流動(dòng)傳熱的CFD-DEM模擬

葛少恒，陳德珍，袁?郡，尹麗潔，梅振飛

豎直降落式反應(yīng)器中顆粒流流動(dòng)傳熱的CFD-DEM模擬

葛少恒，陳德珍，袁?郡，尹麗潔，梅振飛

(同濟(jì)大學(xué)機(jī)械與能源工程學(xué)院熱能與環(huán)境工程研究所，上海 200092)

掌握豎直降落式熱解反應(yīng)器中污泥顆粒的流動(dòng)與傳熱規(guī)律是設(shè)計(jì)反應(yīng)器的關(guān)鍵．采用計(jì)算流體力學(xué)與離散元(CFD-DEM)耦合的方法模擬了豎直降落式反應(yīng)器中污泥顆粒的流動(dòng)與傳熱過(guò)程．首先，利用文獻(xiàn)中的實(shí)驗(yàn)結(jié)果與傳熱模擬結(jié)果比對(duì)，驗(yàn)證了傳熱模型的合理性．進(jìn)而考慮了熱氣流輸入、反應(yīng)器尺寸和污泥顆粒變化等對(duì)傳熱特性的影響．其中顆粒的質(zhì)量隨溫度的變化規(guī)律通過(guò)熱重實(shí)驗(yàn)測(cè)得，進(jìn)而通過(guò)擬合得到顆粒的粒徑與溫度的關(guān)系．計(jì)算結(jié)果表明，顆粒在反應(yīng)器內(nèi)的流動(dòng)接近于“活塞流”，氣相壓降較大(1329Pa/m)；顆粒與氣相之間的對(duì)流傳熱占據(jù)主導(dǎo)地位，污泥顆粒熱解傳熱過(guò)程顆粒形變的影響較大、不可忽略，反應(yīng)器管徑超過(guò)60mm以后，管徑增大會(huì)導(dǎo)致壁面?zhèn)鳠嵴急燃眲∠陆担擞?jì)算為污泥豎直降落式熱解反應(yīng)器的設(shè)計(jì)提供了依據(jù)．

CFD-DEM；傳熱；豎直降落式反應(yīng)器；污泥顆粒；顆粒形變

熱解技術(shù)在污泥的減量化、無(wú)害化、資源化方面有較大優(yōu)勢(shì)[1-2]，熱解工藝中熱解器的選擇至關(guān)重要，一方面影響能量的利用效率，另一方面對(duì)于某些特定的固體廢棄物，利用其理化特性，可以選擇最合適及經(jīng)濟(jì)的熱解器以滿足要求．目前已經(jīng)報(bào)道的廢棄物熱解反應(yīng)器有固定床熱解器、回轉(zhuǎn)窯熱解器、流化床熱解器等[3]．而豎直降落式反應(yīng)器借助重力作用自行下落，無(wú)需外力，易于密封，且傳熱面可以靈活布置．對(duì)污泥這種均勻物料比較適合[4]，但是對(duì)于其流動(dòng)傳熱的研究還非常不足．

一般熱解反應(yīng)器都是間接加熱，雖然保證了氧氣的隔絕，但傳熱效率較低．為了提升傳熱效果，可以利用少量的高溫?zé)煔庵苯訉?duì)污泥進(jìn)行加熱，以減少對(duì)受熱面的需求．且高溫?zé)煔膺€能促使焦油裂解；煙氣中水蒸氣能夠?qū)]發(fā)分中的焦油進(jìn)行重整提升熱解氣熱值[5]，此外水蒸氣對(duì)半焦氣化也有重要影響[6]．

污泥在豎直降落管式熱解反應(yīng)器內(nèi)的熱解過(guò)程是一個(gè)稠密氣固反應(yīng)系統(tǒng)．雙歐拉方法將顆粒作為擬流體，不能獲得顆粒層面的信息．歐拉-拉格朗日方法可以追蹤每一個(gè)顆粒的運(yùn)動(dòng)軌跡，在顆粒尺度上給出系統(tǒng)豐富的受力和運(yùn)動(dòng)信息[7]．計(jì)算流體力學(xué)與離散元(CFD-DEM)方法使用軟球模型計(jì)算顆粒間相互作用，因顆粒碰撞的搜索耗費(fèi)計(jì)算資源，再耦合傳熱、傳質(zhì)和化學(xué)反應(yīng)計(jì)算資源的消耗會(huì)非常巨大．在移動(dòng)床熱轉(zhuǎn)化過(guò)程的數(shù)值模擬方面，國(guó)內(nèi)外采用CFD-DEM方法的研究較少，使用CFD-DEM方法研究流化床的較多，但主要集中在二維和準(zhǔn)三維空間[8].

此外，熱解過(guò)程中顆粒的收縮對(duì)顆粒的傳熱和運(yùn)動(dòng)都有一定影響[9]．由于污泥顆粒熱解過(guò)程復(fù)雜，目前關(guān)于污泥顆粒熱解過(guò)程的模擬較少，多以單顆粒為主，劉秀如[10]研究單個(gè)污泥顆粒在流化床中的受熱升溫過(guò)程，分析了顆粒粒徑等對(duì)顆粒內(nèi)部各點(diǎn)升溫過(guò)程的影響規(guī)律．景亮晶[11]建立了生物質(zhì)單顆粒的熱解模型，研究了外界溫度、顆粒粒徑等對(duì)顆粒內(nèi)部的傳熱影響．這些研究均假設(shè)污泥顆粒為均勻球體，不考慮顆粒之間的碰撞和傳熱，也不考慮熱解過(guò)程中顆粒體積的變化．

本文基于CFD-DEM模型，考慮煙氣和污泥顆粒對(duì)流換熱、反應(yīng)器壁面和顆粒之間導(dǎo)熱、顆粒之間導(dǎo)熱以及顆粒在熱解過(guò)程的粒徑變化，對(duì)豎直降落式熱解器中污泥顆粒熱解過(guò)程進(jìn)行模擬，追蹤不同顆粒的運(yùn)動(dòng)軌跡，分析反應(yīng)器內(nèi)顆粒傳熱和溫度分布，為豎直降落式熱解反應(yīng)器的設(shè)計(jì)提供依據(jù)．

1?模型介紹

基于CFD-DEM方法，氣相用Navier-Stokes方法描述，顆粒相采用DEM方法跟蹤體系內(nèi)的每個(gè)顆粒[12]．污泥顆粒的熱解過(guò)程非常復(fù)雜，本文主要研究豎直降落式熱解反應(yīng)器內(nèi)部顆粒的運(yùn)動(dòng)和傳熱特性，對(duì)模型有如下假設(shè)：①污泥顆粒為球形、干顆粒；②不考慮污泥顆粒內(nèi)部的溫度梯度；③污泥顆粒熱解過(guò)程中，粒徑均勻減?。?/p>

1.1?氣相控制方程

對(duì)于密相系統(tǒng)內(nèi)的冷態(tài)和熱態(tài)氣固流動(dòng)模擬，絕大多數(shù)研究者將流體運(yùn)動(dòng)假設(shè)為層流[13]，且本文工況雷諾數(shù)較小，所以選擇層流模型．

氣相控制方程主要包括連續(xù)性方程、動(dòng)量方程及能量方程：

式中：g為氣體的比熱容；為氣體熱擴(kuò)散率；g,i為氣相與顆粒之間的熱流；g，wall為氣相與壁面間熱流．

1.2?固相控制方程

顆粒的運(yùn)動(dòng)采用離散單元法(DEM)描述[14]，顆粒在平動(dòng)和旋轉(zhuǎn)運(yùn)動(dòng)下的運(yùn)動(dòng)方程如下：

式中：、、i和分別是單個(gè)顆粒的質(zhì)量、慣性矩、線速度和角速度．g和g,i是重力和流體-顆粒相互作用力．彈性力e,ij和阻尼力d,ij用于表征顆粒間力．顆粒作用在顆粒上的扭矩包括：由切向力產(chǎn)生并引起顆粒旋轉(zhuǎn)的,ij，和由不對(duì)稱法向接觸力產(chǎn)生的,ij．當(dāng)法向力不通過(guò)粒子中心時(shí)，應(yīng)生成,ij．

流體-顆粒間曳力計(jì)算使用Gidaspow[15]模型．

顆粒的能量方程[14]為

式中：c,i和T分別為顆粒的比熱容和溫度；,j為顆粒和顆粒之間的導(dǎo)熱熱流；,g為顆粒與周圍氣相的對(duì)流換熱熱流；,rad顆粒與周圍環(huán)境的輻射換熱．

1.3?傳熱模型

在式(3)和(6)中包括了3種熱傳遞方式，即對(duì)流換熱(顆粒與氣相、氣相與壁面)、導(dǎo)熱(顆粒與顆粒、顆粒與壁面)和輻射換熱(顆粒與周圍環(huán)境)．煙氣的發(fā)射率[16]和溫度較低，因此忽略輻射換熱．

1.3.1?對(duì)流換熱

顆粒與氣相之間的對(duì)流換熱量, g為：

式中，和g, i為顆粒表面積和氣相溫度．

氣相與近壁面?zhèn)鳠崃縢,wall為

式中，g,wall和g,wall為傳熱系數(shù)和氣相與壁面接觸面積．

1.3.2?導(dǎo)熱換熱

顆粒之間導(dǎo)熱使用Chaudhuri等[17]提出的模型，顆粒之間導(dǎo)熱通過(guò)重疊面積傳遞，如圖1所示，相互碰撞的顆粒之間的傳熱計(jì)算公式如下：

其中，c為熱傳導(dǎo)系數(shù)．

其中，k、k是接觸顆粒，的導(dǎo)熱系數(shù)；F是顆粒間的接觸力；*是接觸顆粒、的幾何平均半徑；*是等效彈性模量．

圖1?顆粒導(dǎo)熱模型

1.4?顆粒粒徑變化模型

假定顆粒的密度不變，顆粒的粒徑隨著質(zhì)量的減小，通過(guò)熱重實(shí)驗(yàn)(WRT-3P 型熱分析儀，上海精密科學(xué)儀器有限公司)，測(cè)得顆粒的質(zhì)量變化與溫度的關(guān)系，擬合得到顆粒粒徑與溫度的關(guān)系．

由此得到：

式中，是污泥顆粒轉(zhuǎn)化率，%；是指前因子；是通用氣體常數(shù)；是顆粒溫度；是升溫速率；是活化能；0和new分別是顆粒初始粒徑和反應(yīng)過(guò)程中變化的粒徑；end是熱解實(shí)驗(yàn)終止時(shí)樣品的失重率，%．

圖2為質(zhì)量為7.4mg、粒徑為2.16mm的干污泥顆粒在加熱速率為7℃/min，反應(yīng)終溫為650℃的條件下的TG曲線．從圖中可以看出計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果較為吻合．

圖2?干燥污泥顆粒的TG曲線

2?模擬工況

模擬對(duì)象為豎直降落式熱解器，如圖3所示，長(zhǎng)170mm，直徑60mm，為節(jié)省計(jì)算資源，選取的長(zhǎng)度較短，顆粒從頂端落入，同時(shí)顆粒床層向下移動(dòng)，顆粒從底端流出．入口熱煙氣從頂端流入，流網(wǎng)格尺寸為2.5mm．顆粒、流體以及其他參數(shù)均列于表1．

圖3?豎直降落式熱解器示意

表1?模擬參數(shù)及工況

Tab.1?Simulation parameters and working conditions

3?結(jié)果與討論

3.1?模型驗(yàn)證

為驗(yàn)證采用的CFD-DEM及傳熱子模型的正確性，將模型計(jì)算的結(jié)果與Patil等[18]做的鼓泡床內(nèi)氣固傳熱實(shí)驗(yàn)結(jié)果進(jìn)行比對(duì)，驗(yàn)證本文傳熱模型在復(fù)雜密相系統(tǒng)內(nèi)的適用性．參數(shù)設(shè)定見(jiàn)表2．

表2?準(zhǔn)二維鼓泡床驗(yàn)證算例工況設(shè)置

Tab.2 Setting of quasi-two-dimensional bubbling bed verification case

圖4給出了Patil等[18]實(shí)驗(yàn)中所采用的準(zhǔn)二維鼓泡床幾何模型示意．溫度為363.15K的顆粒堆積在床層下部，隨后溫度為293.15K的氣體流化顆粒．入口設(shè)置為均勻速度，壁面為無(wú)滑移；出口設(shè)置為環(huán)境大氣壓．工況設(shè)置列于表1中．根據(jù)Patil等[18]的實(shí)驗(yàn)，對(duì)于速度的影響，研究1.20m/s、1.54m/s和1.71m/s 3組工況．

圖5比對(duì)不同工況下實(shí)驗(yàn)和模擬中鼓泡床內(nèi)平均顆粒溫度的演變(＝75.0g，p＝1.0mm)．如圖所示，增大流體速度會(huì)增大每個(gè)顆粒的對(duì)流傳熱速率，因此顆粒溫度降低較快．此外，較大的流體速度使得顆粒運(yùn)動(dòng)更為劇烈，顆粒碰撞頻率增加，顆粒-顆粒導(dǎo)熱以及顆粒-流體-顆粒導(dǎo)熱也相應(yīng)增加，帶走系統(tǒng)內(nèi)更多熱量，顆粒較快降溫．從圖5中可以看出模擬的顆粒平均溫度與實(shí)驗(yàn)值吻合情況良好，說(shuō)明本文中的傳熱模型適用于稠密堆積兩相傳熱．

圖4?準(zhǔn)二維鼓泡床幾何模型示意

Fig.4 Schematic of geometric model of quasi-two-dimensional bubbling bed

3.2?溫度分布及顆粒運(yùn)動(dòng)過(guò)程中溫度變化

圖6(a)為本文反應(yīng)器中顆粒的溫度分布，常溫顆粒落入反應(yīng)器中，與熱煙氣接觸，溫度迅速上升，同時(shí)隨著顆粒床向下移動(dòng)，移動(dòng)的過(guò)程繼續(xù)受到熱煙氣和壁面的加熱緩慢升溫．最終流出反應(yīng)器的顆粒溫度較為均勻，但貼近壁面處溫度稍高，中心處溫度稍低．圖6(b)為氣相溫度分布，氣相進(jìn)入反應(yīng)器后與顆粒換熱迅速降溫，同時(shí)與壁面換熱被壁面加熱，形成邊界溫度高、中心溫度較低的溫度分布．圖6(c)顯示了流體流動(dòng)過(guò)程的壓力變化，從進(jìn)入床層到流出床層，流體壓降為1329Pa/m，說(shuō)明這種形式的反應(yīng)器阻力相對(duì)較大，設(shè)計(jì)中長(zhǎng)度不應(yīng)過(guò)長(zhǎng)．

圖7為反應(yīng)器中顆粒的流動(dòng)狀態(tài)，初始階段，床層上顆粒流動(dòng)狀態(tài)基本上是平推流，在重力作用下顆粒平推向下流動(dòng)，移動(dòng)床表層的形狀基本保持不變；流動(dòng)一定高度后，床層中心顆粒與壁面處顆粒速度差越來(lái)越大，顆粒流動(dòng)逐漸向匯聚流轉(zhuǎn)變．

圖7?反應(yīng)器中顆粒不同時(shí)刻的流動(dòng)狀態(tài)

3.3?不同傳熱方式的影響

圖8分別給出了貼近壁面處與中心處顆粒的溫度及熱流變化．兩顆粒起始階段溫度均迅速上升，后壁面處顆粒緩慢上升至穩(wěn)定，中心處顆粒溫度穩(wěn)定后又有一定的上升，這是由于煙氣在初始階段加熱顆粒后降溫，后又被壁面加熱升溫，反過(guò)來(lái)又對(duì)顆粒進(jìn)行了傳熱，從圖6(b)可以看出．圖中還可以看到對(duì)于兩種顆粒，對(duì)流換熱均占據(jù)主導(dǎo)地位，但是對(duì)于壁面處顆粒，壁面導(dǎo)熱也不可忽略，且在54s后超過(guò)了對(duì)流換熱量，這主要是煙氣的溫度下降而壁面的溫度保持恒定造成的；對(duì)于中心處顆粒，包括對(duì)流換熱和顆粒間換熱，在相對(duì)短的時(shí)間內(nèi)，顆粒間換熱是可以忽略的．

圖8?不同位置處顆粒溫度及熱流變化

3.4?粒徑變化對(duì)傳熱的影響

圖9(a)考察了考慮顆粒在熱解過(guò)程中粒徑縮小對(duì)傳熱過(guò)程的影響，相對(duì)于假定粒徑不變，考慮了顆粒形變的過(guò)程顆粒換熱效果更好，出口溫度高出40℃，主要是因?yàn)閾Q熱過(guò)程中顆粒的粒徑隨溫度變化較大，從圖9(b)可以看出粒徑從最初的1.6mm，最小的減小到1.22mm，顆粒粒徑減小，床層單位體積的顆粒表面積相應(yīng)增加，增大了顆粒與煙氣的換熱面積，同時(shí)粒徑的減少，會(huì)增加顆粒與壁面以及顆粒之間的接觸，顆粒間導(dǎo)熱與接觸面積正相關(guān)，強(qiáng)化了導(dǎo)熱．

3.5?壁面溫度對(duì)傳熱的影響

圖10中不同壁面溫度情況下開(kāi)始時(shí)間顆粒溫升速度幾乎一致，這是因?yàn)樵陂_(kāi)始階段煙氣與顆粒溫度差巨大，對(duì)流換熱占據(jù)絕對(duì)主導(dǎo)地位，壁面溫度的影響可以忽略，但是隨著顆粒向下運(yùn)動(dòng)，壁面溫度的影響開(kāi)始顯現(xiàn)，高壁溫反應(yīng)器中的顆粒溫度維持相對(duì)較高的溫升速率，主要是由于氣體加熱顆粒溫度下降后，高壁溫壁面對(duì)氣體傳熱，氣體可以維持較高溫度，同時(shí)高壁溫壁面對(duì)顆粒傳熱量也大，所以壁面溫度影響主要體現(xiàn)在顆粒進(jìn)入反應(yīng)器一段時(shí)間以后．

圖9?考慮顆粒形變的溫度演變及粒徑分布

圖10?不同壁面溫度下顆粒溫度的演變

3.6?反應(yīng)器管徑對(duì)傳熱的影響

圖11(a)所示為不同管徑下顆粒溫度變化情況，60mm與80mm的管徑溫差較大，但是管徑繼續(xù)增大到100mm顆粒溫差差別很?。畧D11(b)給出了3種管徑下壁面導(dǎo)熱占總傳熱量的比例，60mm管徑導(dǎo)熱占比較大，同時(shí)加上壁面對(duì)流體的傳熱，所以小管徑情況下顆粒溫度相對(duì)高很多，當(dāng)管徑增大，壁面溫度的影響逐漸減小，管徑的影響存在一個(gè)臨界值，在本模擬工況下為60mm．

圖11?不同反應(yīng)器內(nèi)徑下顆粒溫度演變及壁面導(dǎo)熱占比

4?結(jié)?論

本文對(duì)豎直降落式反應(yīng)器中污泥顆粒的流動(dòng)與傳熱過(guò)程進(jìn)行了模擬，在考慮顆粒之間、顆粒與壁面、顆粒與熱煙氣、熱煙氣與壁面的傳熱基礎(chǔ)上，加入了顆粒粒徑變化的模型，研究了顆粒粒徑、反應(yīng)器尺寸等對(duì)豎直降落式反應(yīng)器傳熱過(guò)程的影響．

(1)驗(yàn)證了CFD-DEM模型以及傳熱子模型對(duì)豎直降落式反應(yīng)器中污泥顆粒的流動(dòng)與傳熱的適用性．

(2)反應(yīng)器中顆粒溫度呈邊界高、中心低的分布，流出反應(yīng)器的顆粒溫度較為均勻；流體在反應(yīng)器中壓降較大，為1329Pa/m．

(3)顆粒進(jìn)入反應(yīng)器后溫度先因強(qiáng)烈的對(duì)流換熱迅速上升，后緩慢上升，在顆粒升溫的過(guò)程中對(duì)流換熱占主導(dǎo)地位，壁面溫度主要是通過(guò)影響煙氣溫度進(jìn)而影響顆粒溫度分布的．

(4)污泥顆粒熱解傳熱過(guò)程顆粒形變的影響較大，不可忽略，反應(yīng)器管徑超過(guò)60mm以后，管徑增大會(huì)導(dǎo)致壁面導(dǎo)熱占比急劇下降．

［1］ Shao Jingai，Yan Rong，Chen Hanping，et al. Pyrolysis characteristics and kinetics of sewage sludge by thermogravimetry Fourier transform infrared analysis[J].，2008，22(1)：38-45.

［2］ 胡艷軍，吳亞男，高?濤，等. 污泥熱解中殘焦表面官能團(tuán)結(jié)構(gòu)演化特征[J]. 燃燒科學(xué)與技術(shù)，2018，24(2)：126-131.

Hu Yanjun，Wu Ya’nan，Gao Tao，et al. Characteristics of surface functional groups structure of residues during sewage sludge pyrolysis[J].，2018，24(2)：126-131(in Chinese).

［3］ 諸黃清. 污泥熱解過(guò)程中的傳熱特性及熱解殘焦種植研究[D]. 上海：同濟(jì)大學(xué)機(jī)械與能源工程學(xué)院，2016.

Zhu Huangqing. Study on Heat Transfer Characteristics of Sewage Sludge Pyrolysis Process and Biochar Applied to Cultivation[D]. Shanghai：School of Mechanical Engineering，Tongji University，2016(in Chinese).

［4］ Chen D Z，Yin L J，Wang H，et al. Pyrolysis technologies for municipal solid waste：A review[J].，2014，34：2466-2486.

［5］ Mei Zhenfei，Chen Dezhen，Zhang Jixuan，et al. Sewage sludge pyrolysis coupled with self-supplied steam reforming for high quality syngas production and the influence of initial moisture content[J].，2020，106：77-87.

［6］ Franco C，Pinto F，Gulyurtlu I，et al. The study of reactions influencing the biomass steam gasification process[J].， 2003，82(7)：835-842.

［7］ 謝?俊，鐘文琪，邵應(yīng)娟，等. 流化床半焦燃燒的DEM三維數(shù)值模擬[J]. 工程熱物理學(xué)報(bào)，2019，40(6)：1307-1312.

Xie Jun，Zhong Wenqi，Shao Yingjuan，et al. 3D CFD-DEM simulation of char combustion in fluidized bed[J].，2019，40(6)：1307-1312(in Chinese).

［8］ Ku X K，Tian L，Lovas T，et al. CFD-DEM simulation of biomass gasification with steam in a fluidized bed reactor[J].，2015，122：270-283.

［9］ Bruchmüller J，Wachem B G M，Gu S，et al. Modeling the thermochemical degradation of biomass inside a fast pyrolysis fluidized bed reactor[J].，2012，58(10)：3030-3042.

［10］ 劉秀如. 城市污水污泥熱解實(shí)驗(yàn)研究[D]. 北京：中國(guó)科學(xué)院研究生院(工程熱物理研究所)，2011.

Liu Xiuru. Experimental Study on Sewage Sludge Pyro-lysis[D]. Beijing：Graduate School of Chinese Acad-emy of Sciences(Institute of Engineering Thermophys-ics)，2011(in Chinese).

［11］ 景亮晶. 木材熱解過(guò)程中單顆粒內(nèi)部傳熱模型的建立與研究[D]. 北京：北京林業(yè)大學(xué)，2011.

Jing Liangjing. Heat Transfer Modeling and Study in a Single Wood Particle During the Pyrolysis Process[D]. Beijing：Beijing Forestry University，2011(in Chinese).

［12］ 劉?陽(yáng)，陸慧林，劉文鐵，等. 氣固流化床的離散顆粒運(yùn)動(dòng)-碰撞解耦模型與模擬[J]. 燃燒科學(xué)與技術(shù)，2003，9(6)：551-555.

Liu Yang，Lu Huilin，Liu Wentie，et al. Numerical simulations and discrete particle motion-collision decoupled model of gas-particle flow in bubbling fluidized beds[J].，2003，9(6)：551-555(in Chinese).

［13］ Müller C R，Holland D J，Sederman A J，et al. Granular temperature：Comparison of magnetic resonance measurements with discrete element model simulations [J].，2008，184(2)：241-253.

［14］ Cundall P A，Strack O D L. A discrete numerical model for granular assemblies[J].，1979，29(1)：47-65.

［15］ Chalermsinsuwan B，Gidaspow D，Piumsomboon，et al. Two- and three-dimensional CFD modeling of Geldart a particles in a thin bubbling fluidized bed：Comparison of turbulence and dispersion coefficients[J].，2011，171：301-313.

［16］ Yang W J，Zhou Z Y，Yu A B. Particle scale studies of heat transfer in a moving bed [J].，2015(281)：99-111.

［17］ Chaudhuri B，Muzzio F J，Tomassone M S. Modeling of heat transfer in granular flow in rotating vessels[J].，2006(61)：6348-6360.

［18］ Patil A V，Peters E A J F，Sutkar V S，et al. A study of heat transfer in fluidized beds using an integrated DIA/PIV/IR technique [J].，2015，259：90-106.

CFD-DEM Modeling of Heat Transfer Characteristics of Granular Flow in a Vertical Falling Reactor

Ge Shaoheng，ChenDezhen，YuanJun，YinLijie，Mei Zhenfei

(Thermal and Environmental Engineering Institute，School of Mechanical Engineering，Tongji University，Shanghai 200092，China)

The knowledge of the flow and heat transfercharacteristicsof sludge particles in a vertical falling reactor is the key to reactor design. The coupling of computational fluid dynamics and discrete element method(CFD-DEM)is carried out to investigate the flow and heat transfer characteristics of sludge particles in the vertical falling reactor. First，thesimulation results are compared with the experimental results from the literature，and the heat transfer model is verified. Then，the effects of hot air input，reactor size and sludge particle deformation on the heat transfer characteristics are comprehensively explored. The relationship between the mass change of particles and temperature is obtained by thermogravimetric experiment，and the function between diameter and temperature is fitted. The simulation results indicate that the flow of sludge particles approximates to a kind of plug flow. Additionally，the convective heat transfer between particles and gas phase plays a leading role and the pressure drop of gas phase is relatively high(1329Pa/m). Furthermore，the impact of the deformation of particles cannot be ignored. After the reactor tube diameter exceeds 60mm，an increase in pipe diameter will cause a sharp drop in the proportion of wall heat transfer. This simulation provides a basis for the design of thesludge vertical falling reactor.

CFD-DEM；heat transfer；vertical falling reactor；sludge particle；particle deformation

TK6

A

1006-8740(2022)02-0126-07

10.11715/rskxjs.R202202024

2021-04-13.

國(guó)家自然科學(xué)基金資助項(xiàng)目(51776141).

葛少恒（1989—??），男，碩士，助教，ericge@#edu.cn.

陳德珍，女，博士，教授，chendezhen@#edu.cn.

(責(zé)任編輯：梁?霞)