差速流化床密相區內流動特性的數值模擬

宋興飛， 王慶功， 別如山， 紀曉瑜， 陳 佩， 劉茜茜

（哈爾濱工業大學 能源科學與工程學院，哈爾濱150001）

符號說明：

差速流化床采用多床面的結構布置，主、副床的非均勻布風形式使流化床底部密相區形成強烈的床料內循環，可實現不同粒徑燃料的分床燃燒.差速流化床具有物料混合強烈、燃料適應性廣、燃燒效果好、床溫穩定和負荷調節范圍大等優點，且能夠有效控制大氣污染物的生成及排放，尤其在燃用劣質燃料和生物質時具有獨特的優勢，應用前景光明［1－2］.

與常規循環流化床［3］相比，差速流化床內的流動過程比較復雜，且伴有強烈的內混特點，在運行過程中，主、副床風速，顆粒粒徑，主、副床高度差和副床布風板傾角等操作條件對床內的流動特性有很大影響［4－5］.筆者對差速流化床密相區內的流動過程進行數值模擬，分析氣固兩相流體特性，研究主、副床風速對流動特性的影響，為差速流化床的設計、運行和優化提供必要的參考依據.

1 氣固兩相流動模型

模擬對象為12MW差速流化床底部密相區二維流場，模擬區域的幾何結構見圖1，底部中間為主床布風板，兩側為副床布風板，采用四邊形網格，對近壁面處的網格進行加密，總網格數為34萬.

圖1 差速流化床密相區的幾何結構（單位：mm）Fig.1 Geometric structure of dense－phase zone in a high－low CFB boiler（unit：mm）

1.1 數值模型

氣固兩相流動滿足的質量和動量守恒方程［6－8］如下：

質量守恒方程

動量守恒方程

1.2 計算所需參數、邊界條件及算法

煤灰顆粒相屬性如下：形狀假設是球形，單一粒徑為0.6mm，密度為2 000kg／m3，填充值為0.63，顆粒碰撞恢復系數為0.3，初始床高為1 000mm.空氣氣相屬性如下：壓力為1.013 25×105Pa，動力黏度為1.785×10－5kg／（m·s）.

邊界條件如下：主、副床風入口均為速度入口，出口為壓力出口，采用無滑移壁面.

采用Simple算法［9］，時間步長取1×10－4s，計算時長為60s，后30s作為時均值計算樣本.

2 計算結果分析

2.1 差速流化床密相區內的氣固流動特性

根據實際運行工況，主床風速u1＝4.0m／s，副床風速u2＝2.0m／s，主、副床的高度差△h＝500 mm，副床布風板傾角ω＝5°，通過數值模擬得到密相區內的流場.不同高度H（距主床布風板距離）處的時均流動特性見圖2.

由圖2（a）可知，H＝0.8m處副床的時均顆粒體積分數較高，且越靠近邊壁時均顆粒體積分數越高；H＝1.5m處副床的時均顆粒體積分數稍微高于主床，整體來看此處兩床的時均顆粒體積分數大致相同；H＝2.2m處，除靠近邊壁處外，副床的時均顆粒體積分數略低于主床.總體來看，各高度處主床的時均顆粒體積分數均較低，越靠近邊壁時均顆粒體積分數越高，邊壁效應越明顯.由圖2（b）和圖2（c）可以看出，H＝0.8m 處密相區左半面（W ＝－2.8～0m）的時均顆粒橫向速度為正，右半面（W＝0～2.8m）的時均顆粒橫向速度為負（時均顆粒橫向速度由左至右為正，由右至左為負），大部分顆粒（靠近邊壁顆粒除外）的時均顆?？v向速度為正（時均顆?？v向速度由下至上為正，由上至下為負），說明該位置處的顆粒不斷地沿副床布風板回流到主床，且處于上升過程.H＝1.5m處時均顆粒橫向速度和縱向速度有正有負，且時均顆粒橫向速度的絕對值較小，時均顆?？v向速度的絕對值較大，顆粒主要進行縱向運動，內循環中心在該高度處.H＝2.2 m處密相區左半面的時均顆粒橫向速度為負，右半面的時均顆粒橫向速度為正，顆粒不斷地從主床涌向副床，該高度為床層沸騰高度.

由圖2（d）可知，H＝2.2m處主床壓力高于副床，形成了頂部從主床到副床的壓力差，該壓力差使顆粒涌向副床；H＝0.8m處副床壓力高于主床，形成了底部從副床到主床的壓力差，該壓力差使顆粒從副床回流到主床，主、副床間的壓力差為密相區顆粒的內循環流動提供了動力.差速流化床進行冷態流化試驗時，保持其他條件不變，在主床不流化時副床的空氣流量明顯大于主床流化時的情況，這是由于主、副床同時流化時，兩床頂部間存在的壓力差使得大量顆粒涌向副床，增大了副床阻力，從而導致副床空氣流量減小.

圖2 不同高度處顆粒或混合物的時均流動特性Fig.2 Time－averaged flow characteristics of particles or airparticle mixture at different heights

綜上所述，差速流化床密相區內顆粒流動過程可描述為：主床中顆粒在主床高速風的作用下向上運動至床面（H＝2.2m），由于主、副床間存在壓力差，使得主床顆粒翻騰到副床上，然后向下運動至副床底部（H＝0.8m），導致副床底部的壓力高于主床，顆粒重新返回到主床，形成了差速流化床密相區內特有的顆粒內循環，內循環中心在H＝1.5m處.

2.2 主床風速對密相區流動特性的影響

保持副床風速u2＝1.5m／s不變，當主床風速u1分別為3.0m／s、3.5m／s、4.0m／s和4.5m／s時，得到不同工況下的氣固流動過程.H＝1m處主、副床上時均顆粒體積分數隨主床風速的變化見圖3.由圖3可知，隨著主床風速的增大，主床時均顆粒體積分數降低，而副床時均顆粒體積分數呈升高的趨勢；當主床風速大于3.5m／s時，主床時均顆粒體積分數較低，副床時均顆粒體積分數較高，表明回流過程更強烈、內循環特性更明顯.

圖3 主、副床上時均顆粒體積分數隨主床風速的變化Fig.3 Time－averaged volume fraction of particles in primary and secondary bed under different primary bed velocities

氣泡通過床層會導致床內壓力降增大，瞬時的高壓力降說明床內出現氣泡，壓力降的變化反映了床內氣固的流動狀態，壓力降變化越大，床內氣固兩相流動越強烈，顆?；旌显骄鶆?壓力降方差可通過瞬時壓力降計算得到，壓力降方差的大小表明床內顆粒和氣體流動的程度，壓力降方差越大，氣固兩相混合效果越好［10］.圖4給出了主、副床壓力降方差隨主床風速的變化.由圖4可知，隨著主床風速的增大，主、副床壓力降方差均逐漸增大，顆粒和氣體流動越強烈，氣固混合效果越好，內循環越明顯.當u1為3.5～4.5m／s時，3種工況下的顆?；亓魉俣认嘟颐黠@大于u1＝3.0m／s時的回流速度（見圖5）.因此，在所研究的范圍內，保證床內顆粒和氣體維持良好的內循環過程及混合效果的最佳主床風速應為3.5～4.5m／s.

圖4 不同主床風速下的壓力降方差Fig.4 Standard deviation of pressure drop under different primary bed velocities

圖5 不同主床風速下顆粒的回流速度Fig.5 Reflux velocity of particles under different primary bed velocities

2.3 副床風速對密相區流動特性的影響

圖6 不同副床風速下顆粒的回流速度Fig.6 Reflux velocity of particles under different secondary bed velocities

保持主床風速u1＝4.0m／s，當副床風速u2分別為1.5m／s、2.0m／s和2.5m／s時，得到不同工況下的氣固流動過程，顆?；亓魉俣鹊淖兓妶D6.由圖6可知，主、副床風速差越大，回流速度越大.圖7給出了不同副床風速下的壓力降方差.由圖7可知，隨著副床風速的增大，壓力降方差先增大后減小，其最大值發生在副床風速為2.0m／s的工況下，此時壓力降方差最大，表明床內顆粒和氣體流動最強烈，氣固混合效果最好.

圖7 不同副床風速下的壓力降方差Fig.7 Standard deviation of pressure drop under different secondary bed velocities

3 結 論

（1）差速流化床密相區內的顆粒流動過程可描述為：密相區內的顆粒相在高、低速風作用下開始流化，上升至床面，由于主、副床間存在壓力差，使得顆粒從主床表面翻騰到副床表面，并導致副床壓力升高，顆粒下降至副床底部，主、副床底部的壓力差迫使副床上的顆?；亓鞯街鞔?，形成差速流化床密相區內特有的顆粒內循環過程.

（2）保持副床風速u2＝1.5m／s不變，隨著主床風速的增大，主床的時均顆粒體積分數降低，副床的時均顆粒體積分數升高，氣固混合更為激烈.在所研究的范圍內，主床風速最佳值應為3.5～4.5m／s.

（3）保持主床風速u1＝4.0m／s不變，隨著副床風速的增大，回流速度減小，壓力降方差先增大后減小，當u2＝2.0m／s時壓力降方差達到最大值，此時氣固流動最強烈，混合效果最好.

［1］劉洪鵬，肖劍波，李惟毅，等.65t／h高低差速循環流化床流動特性模擬［J］.化工進展，2013，32（2）：290－294，345.LIU Hongpeng，XIAO Jianbo，LI Weiyi，et al.Numerical simulation of flow characteristics in a 65t／h high－low bed CFB［J］.Chemical Industry and Engineering Progress，2013，32（2）：290－294，345.

［2］陳玉村，陳晗霞，余更孫.淺談65t／h燃低劣油頁巖差速流化床鍋爐的設計［J］.工業鍋爐，2008（5）：15－17，25.CHEN Yucun，CHEN Hanxia，YU Gengsun.Brief talk to the design of 65t／h oil－shale－fired circulating fluidized bed boiler with high－low bed［J］.Industrial Boiler，2008（5）：15－17，25.

［3］駱仲泱，何宏舟，王勤輝，等.循環流化床鍋爐技術的現狀及發展前景［J］.動力工程，2004，24（6）：761－767.LUO Zhongyang，HE Hongzhou，WANG Qinhui，et al.Status quo－technology of circulating fluidized bed boiler and its prospects of development［J］.Power Engineering，2004，24（6）：761－767.

［4］秦兵，陳晗霞，鄧勇.高低差速流化床鍋爐內循環流體動力特性及技術特點［J］.工業鍋爐，2007（3）：22－25.QIN Bing，CHEN Hanxia，DENG Yong.The internal circulation dynamic flow characteristic and technical characteristic of high－low bed boiler［J］.Industrial Boiler，2007（3）：22－25.

［5］桂北芳.高低差速流化床鍋爐在設計運行中關鍵技術問題探討［J］.工業鍋爐，2006（3）：18－21.GUI Beifang.The study of problems that occur during the design and operation on high low bed boiler［J］.Industrial Boiler，2006（3）：18－21.

［6］GIDASPOW D.Multiphase flow and fluidization：continuum and kinetic theory descriptions［M］.California，USA：Academic Press Inc，1994.

［7］KOCH D L，SANGANI A S.Particle pressure and marginal stability limits for a homogeneous monodisperse gas－fluidized bed：kinetic theory and numerical simulations［J］.Journal of Fluid Mechanics，1999，400：229－263.

［8］WEN C Y，YU Y H.A generalized method for prediction the minimum fluidized velocity［J］.AICHE Journal，1966，12（3）：610－612.

［9］陶文銓.數值傳熱學［M］.2版.西安：西安交通大學出版社，2001.

［10］趙云華，曾令艷，陸慧林，等.管式布風裝置對流化床焚燒爐內流化特性的影響［J］.過程工程學報，2007，7（3）：451－456.ZHAO Yunhua，ZENG Lingyan，LU Huilin，et al.Fluidizing behavior in a fluidized bed incinerator with tubular distributor［J］.The Chinese Journal of Process Engineering，2007，7（3）：451－456.