基于CFD-DEM的高聚物顆粒流化干燥特性研究

王曉靜 于 雪

（天津大學(xué)化工學(xué)院）

流化床干燥器因其高效的傳質(zhì)傳熱效率、均勻的床層溫度分布及靈活的操作條件等優(yōu)勢，廣泛地應(yīng)用于物料的干燥環(huán)節(jié)。 了解物料在流化干燥過程中的干燥特性對于提高干燥效率、改善干燥系統(tǒng)設(shè)計和優(yōu)化干燥工藝有重要意義。

顆粒的流化干燥機理十分復(fù)雜，實驗手段往往難以獲得氣固兩相在干燥過程中的詳細(xì)信息。而隨著計算機技術(shù)和數(shù)值模擬的發(fā)展，人們建立了各種數(shù)值模型用于模擬氣固兩相系統(tǒng)， 其中CFD-DEM（計算流體力學(xué)-離散單元法）模型所需經(jīng)驗參數(shù)較少，可以從顆粒尺度上揭示氣固兩相的運動和傳遞特性，成為氣固兩相系統(tǒng)的研究熱點之一［1～5］。

目前有關(guān)氣固兩相傳遞特性的數(shù)值計算，大多集中在熱量傳遞的研究，這對于不存在相變反應(yīng)的過程是適用的， 而對于物料的干燥過程，必須同時考慮熱量和質(zhì)量的傳遞。 為此，筆者在上述研究的基礎(chǔ)上，引入傳質(zhì)模型，結(jié)合CFD-DEM探討某高聚物顆粒在振動流化床中的失水規(guī)律，以期為流化床干燥器的工業(yè)應(yīng)用提供理論依據(jù)。

1 數(shù)值模型

1.1 流體相模型

流體相模型主要包括連續(xù)性方程、 動量方程、能量方程和組分輸運方程，分別為：

需要說明的是，以上能量方程中的壁面和流體的對流換熱通量Qf，w在本研究中不予考慮。

1.2 顆粒相模型

顆粒相模型包括顆粒運動方程和單顆粒能量方程（本顆粒的初始含水率較低，故不考慮干燥過程中顆粒質(zhì)量和粒徑的變化）：

1.3 耦合模型

耦合模型包括Gidaspow曳力模型［8］、Li和Mason顆粒-流體對流傳熱模型［9］、Batchelor和O′Brien顆粒-顆粒導(dǎo)熱模型［10］、水分汽化傳熱模型和傳質(zhì)模型，為節(jié)省篇幅，僅描述水分汽化吸熱模型和傳質(zhì)模型，其余模型具體可參見文獻(xiàn)［11］。

水分汽化吸熱模型為：

單顆粒在無限流體中的傳質(zhì)系數(shù)Shi為：

2 模擬對象及計算參數(shù)

研究對象為三維矩形振動流化床，結(jié)構(gòu)如圖1所示，兩側(cè)面為周期性邊界。

圖1 流化床結(jié)構(gòu)示意圖

計算開始前，先在床體上方隨機生成一定數(shù)量的條狀高聚物顆粒，顆粒在重力作用下自由下落，形成一定高度的堆積層。 模擬計算時，流化床以固定振動參數(shù)振動，從而使顆粒能夠均勻向前運動；床層底部通入熱氣流，通過改變熱氣流的溫度和流速， 探究顆粒在不同條件下的干燥過程。

3 計算結(jié)果與分析

3.1 不同氣體流速對干燥特性的影響

圖2為氣體溫度90 ℃和床層高度10 mm時，不同氣體流速下物料顆粒的干燥特性曲線（其中，水分比為某一時刻顆粒的含水率與初始含水率之比）。

圖2 氣體速度對物料干燥特性的影響

由圖2可以看出，顆粒物料一共經(jīng)歷了預(yù)熱、恒速干燥和降速干燥3個階段， 其中恒速干燥階段所占時間最長，顆粒在恒速階段的干燥速率隨氣速的增大而增大，當(dāng)氣體流速由1.8 m/s增至3.8 m/s時， 物料最大干燥速率由0.51 g/（m2·s）均勻提高至0.76 g/（m·s），干燥終了時刻物料的含水率隨氣體流速的增加而降低， 當(dāng)氣速分別為3.3、3.8 m/s時，顆粒最終含水率低于物料的臨界含水率，顆粒進入降速干燥階段。 其原因主要為，當(dāng)氣速增加時，氣流湍動程度增加，對流換熱系數(shù)增大， 單位時間內(nèi)從氣體傳遞至顆粒的熱量增加，水分蒸發(fā)速率隨之增大。 其次，氣速增加使得顆粒表面的邊界層厚度減小，降低了水蒸氣向熱氣流中擴散的阻力，提高了傳質(zhì)效率。 此外，當(dāng)熱空氣濕度一定時，氣速越大，單位時間內(nèi)流經(jīng)顆粒的氣體流量越大，且?guī)ё叩乃魵饩驮蕉啵瑥亩S持了顆粒表面與熱氣流主體間的水蒸氣濃度差，有利于干燥過程的進行。

3.2 不同氣體溫度對干燥特性的影響

圖3為氣體流速1.8 m/s和床層高度10 mm時，不同氣體溫度下物料顆粒的干燥特性曲線。 由圖3可看出，隨著氣體溫度的提高，顆粒預(yù)熱時間縮短，干燥速率增大， 達(dá)到相同干燥要求時所需干燥時間減少。 當(dāng)氣體溫度由95 ℃增加至110 ℃時，顆粒最大干燥速率由0.82 g/（m2·s）均勻地增加至1.74 g/（m2·s）。 其原因主要是氣體溫度越高，顆粒和氣體間的傳熱溫差越大，傳熱效率越高。 其次，同一氣體被加熱至更高的溫度時， 其相對濕度減小， 因而顆粒表面和熱流體間的水蒸氣濃度差增大，促進了水分的遷移過程，從而提高了干燥效率。

圖3 氣體溫度對物料干燥特性的影響

由圖3b可看出，隨著氣體溫度的增大，干燥速率的增長率先增大到一定值后減小。 在預(yù)熱階段，顆粒從熱氣流中吸收熱量加熱，耗費了一定能量，因此干燥速率較低，不同溫度下顆粒的干燥速率變化不大；隨著顆粒溫度的上升，用于水分蒸發(fā)的氣體熱量增多， 干燥速率不斷增大，干燥速率的增長率也逐漸增大，當(dāng)顆粒溫度升至水分蒸發(fā)溫度時， 氣體熱量全部用于水分蒸發(fā)，因而顆粒干燥效率達(dá)到最大值并保持不變，此階段干燥速率的增長率為一定值；當(dāng)顆粒表面的非結(jié)合水全部蒸發(fā)完全時， 顆粒進入降速干燥階段，干燥過程由內(nèi)部水分向顆粒表面的遷移速率控制，因而顆粒的干燥速率降低，干燥速率的增長率也隨之降低。

3.3 不同粒徑對干燥特性的影響

圖4為氣體速度1.8 m/s、氣體溫度95 ℃和床層高度25 mm時， 同形狀的顆粒在粒徑5～8 mm范圍內(nèi)的干燥特性。 由圖4可看出，顆粒粒徑越小，加熱越快，水分蒸發(fā)速率越高，隨著粒徑的增大，其恒速階段的干燥速率差距逐漸減小，顆粒粒徑為5、6、7、8 mm的最大干燥速率分別為0.89、0.78、0.69、0.63 g/（m2·s），在相同干燥時間內(nèi)，顆粒最終含水率隨粒徑的增大而增大，即干燥時間隨粒徑增大而延長。 其原因主要是形狀相同的顆粒，粒徑越小時，比表面積越大，即傳熱面積越大，單位時間內(nèi)傳遞的熱量增加，水分蒸發(fā)加快。 其次，粒徑越小，顆粒內(nèi)部的自由水能夠及時傳遞至顆粒表面， 維持表面的潤濕狀態(tài)，從而使顆粒保持較高的干燥速率，縮短干燥時間。

圖4 顆粒粒徑對物料干燥特性的影響

圖5展示了不同時刻床層內(nèi)顆粒含水率的變化情況。 由圖5可看出：不同粒徑顆粒均勻分布在床層中，從總體上看，床層含水率從頂層向底層遞減； 對于同一粒徑顆粒， 越接近床層底部，含水率越低；對于不同粒徑顆粒，顆粒粒徑越小，含水率越低，與其所處床層位置無關(guān)，由此說明， 顆粒粒徑對干燥速率的影響比床層位置更為顯著。

圖5 不同時刻顆粒含水率變化

4 干燥動力學(xué)模型

以氣體流度1.8 m/s、氣體溫度95 ℃和床層高度10 mm的干燥工況為例，選用5個常見的薄層干燥動力學(xué)模型［12，13］對模擬數(shù)據(jù)進行擬合，擬合結(jié)果見表1。 擬合后，采用相關(guān)系數(shù)R2、卡方值χ2和均方根誤差（RMSE）作為數(shù)學(xué)評價指標(biāo)，其中，R2越接近于1、χ2和RMSE越低， 說明該模型的擬合效果越好。

表1 不同干燥動力學(xué)模型的擬合參數(shù)

由表1可知， 修正Page模型的擬合效果最好，擬合得到的模型參數(shù)a=1.0170、k=0.0021 和n=1.8663，評價指標(biāo)為：相關(guān)系數(shù)R2=0.9981、卡方值χ2=1.3115×10-4、均方根誤差（RMSE）為0.011 5。

5 結(jié)論

5.1 顆粒干燥過程主要發(fā)生在恒速干燥階段，隨著氣體的流速和溫度提高，氣固兩相間的傳熱和傳質(zhì)效率均提升，顆粒干燥速率增大，達(dá)到相同干燥要求時所需干燥時間減少。

5.2 顆粒粒徑對干燥速率影響顯著，隨著粒徑的增大，顆粒干燥速率降低，干燥時間延長。 對于同一粒徑顆粒，越接近床層底部，含水率越低；對于不同粒徑顆粒，顆粒粒徑越小，含水率越低，與所處床層位置無關(guān)，可以確定顆粒粒徑對干燥速率的影響比床層位置的影響更加顯著。

5.3 選用常見的薄層干燥動力學(xué)模型對模擬數(shù)據(jù)進行擬合，發(fā)現(xiàn)MR=aexp（-ktn）模型能夠很好地描述顆粒的失水規(guī)律，其評價指標(biāo)分別為相關(guān)系數(shù)R2=0.9981、卡方值χ2=1.3115×10-4、均方根誤差（RMSE）為0.011 5。