彭昕，王亦飛，位宗瑤，陳福田，于廣鎖

?

含細長顆粒的洗滌冷卻室內(nèi)的多相分布特性

彭昕，王亦飛，位宗瑤，陳福田，于廣鎖

（華東理工大學潔凈煤技術(shù)研究所，煤氣化及能源化工教育部重點實驗室，上海 200237）

采用改進的直接取樣法，在按幾何尺寸縮小的工業(yè)氣化爐洗滌冷卻室冷模裝置內(nèi)，同時測量不同操作條件下的軸徑向局部固含率和氣含率，對細長顆粒在洗滌冷卻室內(nèi)的多相分布特性進行研究。結(jié)果表明：以下降管出口截面為界，洗滌冷卻室可分為上部氣液固混合區(qū)和下部固液流動區(qū)，其中上部區(qū)域由下降管出口區(qū)、破泡板作用區(qū)和氣墊層區(qū)組成，下部區(qū)域由氣相湍動作用區(qū)、回流區(qū)和二次流動區(qū)組成；在顆粒阻礙效應(yīng)減速沉降和團聚效應(yīng)加速沉降的共同作用下，軸向固含率呈現(xiàn)波動分布；環(huán)隙氣速、固相體積分數(shù)和長徑比的增加均會增強床層的湍動，促進氣體的徑向擴散；操作條件的改變使顆粒的漂移速度發(fā)生改變，徑向固含率分布出現(xiàn)波動；在氣相擾動和回流作用下，二次流動區(qū)呈現(xiàn)環(huán)狀流動，流體和顆粒的“壁面效應(yīng)”使該區(qū)域的固含率呈現(xiàn)中心高邊壁低的特點。

細長顆粒；洗滌冷卻室；多相流；相分布；軸向；徑向；湍流

引 言

洗滌冷卻室作為高溫合成氣的洗滌冷卻裝置，在激冷式氣流床煤氣化技術(shù)中發(fā)揮著重要作用。多噴嘴對置式水煤漿氣化爐洗滌冷卻室采用噴淋和環(huán)隙鼓泡相結(jié)合的復合型鼓泡床[1]。經(jīng)激冷降溫的氣體與固態(tài)渣沿下降管進入洗滌冷卻室液池，在穿越液池的過程中完成氣體的洗滌凈化和氣液固三相分離的工作。洗滌冷卻室內(nèi)部流體流動復雜，涉及氣-液-固三相體系的流動、分離及熱質(zhì)傳遞等諸多復雜的科學問題。

有關(guān)洗滌冷卻室的實驗研究工作多集中在氣液分離[2-5]、降膜流動[6-10]以及內(nèi)構(gòu)件流體力學[11-13]等方面，但對洗滌冷卻室液池內(nèi)的固相流動特性的研究卻鮮有涉及。賀必云[14]以煤渣和石英砂為模擬物對洗滌冷卻室內(nèi)的固相濃度分布進行了實驗研究，將固體流動劃為3個區(qū)域：沉降區(qū)、擴散區(qū)和懸浮區(qū)。吳晅[15]則利用數(shù)值模擬對前者的實驗結(jié)果進行了驗證，數(shù)值計算結(jié)果與實驗結(jié)果的總體趨勢一致；在洗滌冷卻室內(nèi)，顆粒濃度在軸向上呈波動形式分布，由下而上呈下降趨勢。吳晅等[16]對氣固兩相流穿越液池過程中的顆粒在三相體系中的微觀運動過程及分布規(guī)律進行了研究，與Euler-Lagrange框架下的數(shù)值模擬結(jié)果基本一致。付碧華等[17]采用高清攝像儀和壓差傳感器對洗滌冷卻管出口處射流深度和液池內(nèi)氣液界面波動特性進行研究發(fā)現(xiàn)，隨著表觀氣速的增大，射流深度呈指數(shù)式增大；洗滌冷卻管出口處液面波動對床層內(nèi)氣液兩相環(huán)流脈動的影響較為顯著。龔曉波等[18]通過氣液兩相的數(shù)值模擬，量化分析了出口動量與最大下降深度和夾層內(nèi)側(cè)出口氣含率的關(guān)系。

上述學者采用球狀顆粒進行了相關(guān)實驗和模擬工作，然而洗滌冷卻室液池內(nèi)的灰渣顆粒多為非球形顆粒。與球形顆粒相比，非球形顆粒在流體中運動時的阻力系數(shù)存在差異[19]，而阻力系數(shù)的差異將導致顆粒運動特征的改變。非球形顆粒一般具備各向異性的性質(zhì)，其中最為典型的是具有一定長徑比的細長顆粒。Gore等[20]發(fā)現(xiàn)，兩相流中顆粒的粒徑對流體的湍動強度有一定影響，小顆粒會削弱湍流強度而大顆粒則會使湍流強度增強。林建忠等[21]對含柱狀固粒的兩相流模擬發(fā)現(xiàn)，粒子對流場的湍動特性起著抑制作用。而在氣化爐運行過程中發(fā)現(xiàn)，具有一定的長徑比的纖維狀顆粒會在氣化爐洗滌冷卻室的破泡板上橋接、沉積，造成洗滌冷卻室內(nèi)液位不穩(wěn)定，合成氣帶灰、帶液等問題，嚴重影響氣化爐的正常運行。因此，對纖維狀細長顆粒在洗滌冷卻室內(nèi)的流體力學特性研究具有深刻的理論意義和實際應(yīng)用價值。

本實驗在按幾何尺寸縮小的工業(yè)氣化爐洗滌冷卻室冷模裝置內(nèi)進行。采用改進的直接取樣法，同時測量不同操作條件下的軸徑向局部固含率s和氣含率g，分析其變化規(guī)律；研究細長顆粒在洗滌冷卻室液池內(nèi)的分布特性和其對流體湍動的影響，為含細長顆粒的洗滌冷卻室內(nèi)的流體力學特性研究提供借鑒和指導。

1 實驗裝置、材料和方法

實驗裝置如圖1所示，洗滌冷卻室由有機玻璃制成，液池內(nèi)徑為200 mm，下降管內(nèi)徑為70 mm，外徑為80 mm。實驗中采用空氣-水-尼龍短纖作為氣液固三相，短纖由直徑分別為0.1 mm和0.15 mm的PA6長絲采用特殊工具裁剪為3 mm和5 mm而成。空氣由羅茨鼓風機送至洗滌冷卻室氣體頂部入口，流量由轉(zhuǎn)子流量計控制。水由水泵從水槽打出，將尼龍短纖固體顆粒與水混合后由進水口注入液池內(nèi)，并采用下部鼓泡的方式預混。距出水口143、193、243、293、383、423、473、523、618 mm處的位置，由下至上布置9個取樣口。

1—fan；2,12—valve；3—flow-meter；4—air inlet；5—air outlet； 6—downcomer；7—bubble-break plate；8—sample tap；9—water inlet；10—water pump；11—bubble tube；13—water outlet

參照文獻[22-24]設(shè)計了圖2所示的取樣裝置，取樣管和取樣器進口內(nèi)徑均為8 mm，取樣口錐形角角度為45°。取樣時，拉動手柄至取樣管出口處，同時抽動帶標尺的取樣器至固定刻度。通過取樣器上的刻度直接讀取氣相所占體積，得到局部氣含率g（體積比）。將取樣器內(nèi)的固體過濾、烘干、稱重并換算得到局部固含率s（體積比）。

1—handle；2—O-rings；3—sampling tube inlet；4—sampling tube outlet；5—sampler inlet；6—sampler

在初始靜態(tài)液位540 mm的條件下，分別測量不同環(huán)隙氣速g（0.29、0.44、0.59 m?s-1）、固相體積分數(shù)s（0.0094%、0.0281%、0.044%）和長徑比r（20、30、50）時軸徑向局部固含率s和氣含率g。

2 實驗結(jié)果與討論

2.1 不同環(huán)隙氣速g下的軸徑向多相分布

2.1.1 軸向局部固含率和氣含率分布 圖3是不同環(huán)隙氣速下洗滌冷卻室內(nèi)軸向固含率和氣含率分布。從圖3可以看出，以下降管出口截面為分界線，洗滌冷卻室可以分為兩部分區(qū)域：以氣液兩相湍流為主體的上部氣液固混合區(qū)和受氣相影響的下部固液流動區(qū)。在環(huán)隙氣速較低的情況下，上下兩區(qū)域的固含率分布較為均勻；隨著環(huán)隙氣速的增加，軸向固含率呈現(xiàn)梯度分布。由于氣相的攜帶作用導致上部固含率略有增加。如圖4所示，下部固液區(qū)主要以垂直于軸向方向的環(huán)狀流動為主，環(huán)隙氣速的增加勢必導致下部固液流動區(qū)的回流增強。下降管出口處為負浮力的逆壓力梯度突擴流動，在主流和回流流動的影響下會引發(fā)洗滌冷卻室底部的二次流動[18]。二次流動強度增大，既增大了顆粒間的相互碰撞概率，又加速了顆粒的橫向漂移和沉降，因此在固液混合區(qū)上部出現(xiàn)了低氣速固含率大于高氣速固含率的現(xiàn)象，而隨著軸向高度的減小，二次流動的影響減弱，固含率的增加主要以沉降為主，且環(huán)隙氣速越大沉降越快。

如圖3所示，隨著環(huán)隙氣速的增大，氣含率呈現(xiàn)先增大后減小的趨勢，在破泡板附近的氣含率變化尤其劇烈。383 mm采樣點處，環(huán)隙氣速增加，局部氣含率略有增大。這是由于下降管表觀氣速的增大使下降管突擴射流深度增大，進而導致下降管出口處氣含率增大。如圖5所示（白色曲線為氣體的輪廓線，黑色直線為下降管外壁），在負浮力和逆壓力梯度作用下，氣體在下降管外壁聚集并沿邊緣上升，形成羽狀流動，下降管出口處形成的氣泡較大，測量具有一定誤差，氣含率的變化并不明顯。

而在423 mm測量點之后，氣含率的變化則較為明顯，一方面由于破泡板的作用，大氣泡破碎成大量小氣泡；另一方面氣體沿下降管外壁反折向上流動的同時存在徑向擴散。因此，氣含率的變化在破泡板截面形成了轉(zhuǎn)折點。環(huán)隙氣速增大，位于破泡板下部和下降管出口上部區(qū)域，由于氣泡破碎和回流的作用氣含率大幅增加。而破泡板上部區(qū)域氣含率的減小則是由于氣速的增大增加了氣體對液體的擾動，導致該部分的氣墊層高度增加，氣體帶液量增大，出現(xiàn)氣含率減小，含液量增大的現(xiàn)象。

2.1.2 徑向局部固含率和氣含率分布 圖6是g= 0.29 m?s-1時的不同軸向高度的徑向固含率和氣含率分布。對液池中心到近壁面的徑向距離做了量綱1處理。

由圖6可以看出，383 mm采樣點位于下降管出口正下方，此位置氣相劇烈湍動，固含率較低；293 mm采樣點處于氣相體積和流體流速分布不均勻引起的回流區(qū)，呈現(xiàn)下降管外壁附近固含率高，近壁面固含率低的變化；243、193、143 mm采樣點則處于由主流和回流引發(fā)的二次流動區(qū)，該處固含率總體基本呈現(xiàn)中心高邊壁低的趨勢。流體在近壁處的流速大于床層內(nèi)部，壁面附近的空隙率大于床層內(nèi)部，流動阻力較小，該區(qū)域呈環(huán)狀流動，形成流體的“壁面效應(yīng)”。此外，細長顆粒是一類具有主軸的細長體，根據(jù)細長體理論，當流體運動方向和細長體顆粒的主軸不在同一方向上時會引起細長體的轉(zhuǎn)動，這種轉(zhuǎn)動有使細長體遠離壁面的趨勢，形成細長體的“壁面效應(yīng)”[25-28]。兩種“壁面效應(yīng)”的協(xié)同作用導致243 mm處邊壁處的固含率低于中心處固含率。而隨著軸向距離減小，環(huán)狀流動對顆粒的作用減弱，細長顆粒是一類具有主軸的顆粒，其在沉降過程中存在橫向的漂移，這種漂移是脈動的、不規(guī)律的，因此導致193和143 mm兩個采樣點處的徑向固含率變化的規(guī)律性較差，但總體趨勢基本呈現(xiàn)出中心高邊壁低。

從圖6中徑向氣含率隨軸向高度的變化可知，下降管出口處，由于流道突擴導致氣流在床層內(nèi)分布不均，氣體主要沿著下降管外壁向上流動；破泡板截面處，由于破泡作用和對氣流的反折作用，邊壁處氣含率顯著增大，但氣流還是主要集中在下降管外壁面一側(cè)；如圖7所示，氣墊層區(qū)處于湍流發(fā)展的充分段，該區(qū)域主要是氣泡的聚并和破碎，徑向氣含率略有波動，變化較小。

因此，可將固液流動區(qū)分為3個區(qū)域：氣相湍動作用區(qū)、回流區(qū)和二次流動區(qū)。氣液固混合區(qū)則分為下降管出口區(qū)、破泡板作用區(qū)和氣墊層區(qū)。

圖8為不同環(huán)隙氣速下不同區(qū)域的徑向固含率分布。從圖中可看出，氣相湍動作用區(qū)和回流區(qū)處于負浮力淺層突擴射流區(qū)，其徑向固含率變化波動較大。在氣相作用區(qū)，環(huán)隙氣速較小時，突擴射流段較短，氣相主要沿軸向擴散，中心區(qū)域的固含率較高；隨著環(huán)隙氣速的增大，下降管突擴射流深度增大，中心區(qū)域固含率急劇減小。在較強的氣相湍動作用下，氣相徑向擴散范圍增大，氣相的攜帶使得固相從中心向壁面擴散，導致邊壁處附近的固含率較高；當氣速增大到0.59 m?s-1時，氣相引起的液相回流的范圍和強度增大，表現(xiàn)出壁面附近的固含率大幅減小，近下降管區(qū)域的固含率顯著增加。在回流區(qū)，氣速較低時，液相回流強度較小，徑向固含率分布差異較??；隨著環(huán)隙氣速增大，回流范圍和強度增大，表現(xiàn)出中心和近下降管壁面區(qū)域固含率減小的趨勢；環(huán)隙氣速繼續(xù)增大時，氣相引起的液相回流更強，壁面附近的固含率大幅減小，近下降管區(qū)域的固含率則顯著增加。根據(jù)徑向固含率變化，回流區(qū)徑向范圍從下降管中心至反應(yīng)器壁面。在243 mm采樣點處，則出現(xiàn)了完全不一樣的固含率分布，從下降管中心至壁面處逐漸減小。環(huán)隙氣速增大同樣導致近壁面處固含率減小，但與前兩塊區(qū)域不同，流體和細長顆粒的“壁面效應(yīng)”起主要作用。環(huán)隙氣速的增大，必然導致二次流動區(qū)的環(huán)狀流動增強，近壁面處的流體流速也隨之增大。流體的剪切作用使得纖維狀顆粒在壁面附近的周期性轉(zhuǎn)動幅度加大，更難達到穩(wěn)定取向狀態(tài)，更易于產(chǎn)生“pole vaulting”現(xiàn)象，因此環(huán)隙氣速的增大對近壁面處固含率的減小更為明顯。

不同環(huán)隙氣速下不同區(qū)域的徑向氣含率分布如圖9所示。由下降管出口區(qū)至氣墊層區(qū)，徑向氣含率的差異逐漸明顯。在下降管出口區(qū)，由于負浮力突擴射流的特性，僅有/= 0.4和0.5兩點的徑向氣含率具有較為明顯的差異。可以看出，隨著環(huán)隙氣速增大，氣含率顯著升高，而/= 0.4處差異較為不明顯的原因在于該位置的氣泡直徑較大，取樣具有一定誤差。在破泡板作用區(qū)，環(huán)隙氣速增加時，氣相的徑向擴散略有增大，這是由氣相折返和液相回流共同作用的結(jié)果。在氣墊層區(qū)，隨著環(huán)隙氣速增大，床層的湍流強度增強，氣含率呈較為明顯的中心峰分布[29]。

2.2 不同固相體積分數(shù)s下的軸徑向多相分布

2.2.1 軸向局部固含率和氣含率分布 圖10為g= 0.29 m?s-1條件下，不同固相體積分數(shù)下的固含率和氣含率軸向分布。細長體顆粒群的溶液按照濃度不同分為3種：稀釋溶液、半稀釋溶液和濃溶液[30-31]。在本實驗中< (pp)2，屬于稀釋溶液，其中為細長體顆粒的體積分數(shù)，p為直徑，p為長度。在稀釋溶液中，任一個細長體顆粒與其周圍顆粒間的距離均大于p，顆粒間的作用幾乎可以忽略。但考慮到流體中纖維之間的動力學相互作用，即使?jié)舛缺容^低，顆粒也會受到鄰近顆粒的阻礙效應(yīng)[31]。

由圖10可看出，固相體積分數(shù)增加對上部氣液固混合區(qū)的固含率變化基本無影響，此區(qū)域氣液擾動劇烈，固含率分布的差異不明顯。而下部固液混合區(qū)由于阻礙效應(yīng)的影響，固含率的變化則較為明顯。濃度較低時，阻礙效應(yīng)不明顯，固相分布較均勻；隨著濃度增大，顆粒間相互碰撞的概率增大，纖維的交織效應(yīng)促使纖維團聚，顆粒間的絮凝現(xiàn)象加重。纖維聚集后沉降速度增大，因此在293、243、193 mm測量點間的固含率變化較大。顆粒間的阻礙效應(yīng)會隨著濃度的增大而增大，最終導致顆粒的沉降速度減小，因此193和143 mm間的固含率變化明顯減小。這與Kumar等[32]以及Herzhaft等[33]的發(fā)現(xiàn)基本一致。

氣含率則表現(xiàn)出下降管出口區(qū)和破泡板作用區(qū)的氣含率隨固相濃度增大而增大，氣墊層區(qū)氣含率隨固相濃度增大而減小的特點，與環(huán)隙氣速增大時的變化一致，這說明細長顆粒濃度的增大使氣相的湍動增強。根據(jù)Gore等[20]的理論，粒子的大小對兩相流流場的湍流強度有著較大影響。以

為界，p/e< 0.1的粒子會使流場湍流強度降低，p/e> 0.1的粒子則會使流場湍流強度增強。根據(jù)經(jīng)驗公式式(2)～式(4)[34-35]可得式(5)，計算可得在本實驗條件下p/e= 0.00200.0029 < 0.1，屬于小的柱狀粒子，但與上述理論相反，細長顆粒的增加反而使流場的湍動強度增強。其原因可能在于，首先與Gore等[20]的兩相流實驗條件不同，本實驗在含有內(nèi)構(gòu)件的氣液固三相流中進行，流場結(jié)構(gòu)復雜，湍流強度的變化受多種因素共同影響。其次，與Gore等[20]采用的球形顆粒不同，本實驗采用的細長顆粒具有各向異性，細長顆粒上的渦脫落會增加渦量的擴散，增大含能渦的湍動，進而增強氣相的湍動強度。且顆粒體積分數(shù)越大，渦脫落越明顯，對氣相的擾動越大，流場越易失穩(wěn)。此外，濃度增大導致顆粒間的相互作用增強，對顆粒的取向有很大影響，由顆粒取向隨機分布導致的剪切正應(yīng)力的擾動則對流場的失穩(wěn)具有一定的促進作用。

= 0.07(2)

(4)

(5)

式中，為湍動能，m2?s-2；為湍動能耗散率，m2?s-3；為湍流長度尺度，m；為水力直徑，可用式(6)計算，m；C為經(jīng)驗常數(shù)，取0.09；p為顆粒直徑，m；e為含能尺度，m。

=-out(6)

式中，為洗滌冷卻室內(nèi)徑，m；out為下降管外徑，m。

2.2.2 徑向局部固含率和氣含率分布 不同固相體積分數(shù)下不同區(qū)域的徑向固含率分布如圖11所示。隨著濃度增大，氣相湍動作用區(qū)的近壁面固相濃度有增高的趨勢。細長顆粒在運動時存在平動和轉(zhuǎn)動，均會受周圍顆粒的影響。既有其他顆粒通過流體施加的水動力長程影響，又有顆粒間的短程影響，本實驗顆粒間的短程影響可以忽略。顆粒間的相互作用對于漂移速度有(1) 量級的影響，極易導致顆粒分布的不均勻性[36]，該區(qū)域徑向固含率的變化可能是由于濃度增大影響顆粒間長程作用，導致顆粒橫向漂移速度的改變。同時，隨著細長顆粒體積濃度的增大，顆粒間的阻礙效應(yīng)增強，顆粒間的絮凝現(xiàn)象變得嚴重，在氣相突擴射流的卷吸作用下使得近下降管壁面的固含率逐漸增大。濃度由0.0094% 增大到0.028% 時，回流區(qū)的徑向固含率在/= 0.4處出現(xiàn)先減小后增大的轉(zhuǎn)折，濃度繼續(xù)增大時該趨勢保持不變，這應(yīng)該是液相回流和顆粒間相互作用共同影響該處顆粒橫向漂移的結(jié)果。同樣，顆粒漂移速度的變化導致了二次流動區(qū)中心處和近壁面處的固含率差異隨濃度增大而增大。

圖12是不同區(qū)域的徑向氣含率隨固相體積分數(shù)的變化。隨著濃度增高，在下降管出口區(qū)，徑向氣含率在/= 0.4和0.5處差異較大；而在破泡板作用區(qū)和氣墊層區(qū)，整個床層的徑向氣含率分布較為一致。對比氣含率的軸相分布可知，隨著軸向高度的增加，固相體積分數(shù)增加時，不同位置的氣含率分布趨于一致，導致473和618 mm處的徑向氣含率分布差異不明顯。細長顆粒的體積分數(shù)的增加增強了下降管出口處氣體的湍動，使上部床層氣體徑向擴散強度增大，徑向差異較為明顯。

2.3 不同長徑比r下的軸徑向多相分布

2.3.1 軸向局部固含率和氣含率分布 在g= 0.29 m?s-1，s= 0.0094% 的條件下，圖13為不同長徑比下固含率和氣含率軸向分布。其中r= 20和r= 30的細長顆粒具有相同的3 mm長度，r= 30和r= 50則具有相同的0.1 mm直徑，三者單個顆粒的體積為30<50<20。由于氣體更易攜帶細小粒子，因此上部氣液固混合區(qū)的固含率呈現(xiàn)r= 30最大，r= 50次之，r= 20最小的情況，下部液固混合區(qū)的固含率變化則比較復雜。

細長顆粒在流體中沉降時強烈依賴取向分布，當顆粒主軸方向與重力方向相同時，沉降最快，與重力方向垂直時，沉降最慢。隨著長徑比增大，最大沉降速度與最小沉降速度之差增加。Herzhaft等[37]對大長徑比顆粒的沉降研究發(fā)現(xiàn)，顆粒沉降時會結(jié)塊，平均沉降速度比單個垂直顆粒時的情況更大，沉降時具有非常大的各向異性的脈動速度。因此，在二次流動區(qū)上部，r= 50顆粒的固含率最小，隨著軸向高度減小，固含率逐漸增大到最大?？紤]到顆粒本身的慣性，r= 20的顆粒沉降較r= 30快，在二次流動區(qū)的固含率較低。而在回流區(qū)，隨著長徑比增加固含率減小，這可能是顆粒取向變化引起沉降速度的變化所致。

由圖13可知，氣含率表現(xiàn)出下降管出口區(qū)和破泡板作用區(qū)的氣含率隨長徑比增大而增大，氣墊層區(qū)氣含率隨長徑比增大而減小的特點，這說明細長顆粒長徑比的增大使氣相的湍動增強。造成這種現(xiàn)象的原因可能是由于長徑比增加，顆粒的長程水動力作用增大，增強了湍流中的大尺度結(jié)構(gòu)，增大了含能渦的湍動，進而提高了對流場的擾動。

2.3.2 徑向局部固含率和氣含率分布 圖14是不同長徑比下不同區(qū)域的徑向固含率分布。在氣相湍動作用下，隨著長徑比增大，邊壁附近的固含率逐漸增高，這說明顆粒在氣相的作用下，隨長徑比的增加，顆粒更易于向邊壁處橫向漂移。對比3個長徑比不同的單個顆粒的體積可看出，氣相更易于攜帶r= 30和50的顆粒，因此二者的固含率較高。而r= 30顆粒的固含率更大則是由于該顆粒的慣性小，對流體的跟隨性好，顆粒更易發(fā)生聚集。在回流區(qū)，隨著長徑比的增大/= 0.4處的固含率逐漸減小，與近壁面附近的固含率差距逐漸縮小，對比氣相湍動作用區(qū)的徑向固含率變化可知，這主要是由于液相回流所引起的。由于顆粒自身的慣性影響，r= 20的顆粒對流體的跟隨性較差，同r= 30和50的顆粒的固含率變化出現(xiàn)了較大差異。在二次流動區(qū)，大長徑比的顆粒的“壁面效應(yīng)”更加明顯，中心與近壁面的固含率差異較大，且由于沉降時的漂移速度脈動，徑向固含率的變化波動較大。

圖15是不同長徑比下不同區(qū)域的徑向氣含率分布。由圖中可看出，隨著長徑比加大，氣墊層區(qū)的氣含率在徑向變化上尤其明顯。長徑比的增大同樣增強了氣體在床層中的擴散，這種流場的失穩(wěn)可能與顆粒取向隨機分布而導致的剪切正應(yīng)力的擾動和顆粒間長程水動力作用的增加有關(guān)。

3 結(jié) 論

（1）根據(jù)洗滌冷卻室內(nèi)的軸向相分布變化，可將洗滌冷卻室分為兩個區(qū)域：上部氣液固混合區(qū)和下部固液流動區(qū)。上部氣液固混合區(qū)可分為：下降管出口區(qū)、破泡板作用區(qū)和氣墊層區(qū)；下部固液流動區(qū)可分為：氣相湍動作用區(qū)、回流區(qū)和二次流動區(qū)。

（2）細長顆粒的沉降與顆粒取向分布密切相關(guān)。細長顆粒的阻礙效應(yīng)會減小顆粒的沉降速度，而顆粒的團聚則會加大平均沉降速度，軸向固含率呈現(xiàn)波動分布。

（3）環(huán)隙氣速、固相體積分數(shù)和長徑比的增大均會增強床層的湍動程度，促進氣體在床層中的徑向擴散。

（4）顆粒間的長程水動力作用影響細長顆粒的橫向漂移速度，進而影響顆粒在床層中的徑向分布；在氣相湍動作用區(qū)和回流區(qū)，由于氣相和液相回流的共同作用，徑向固含率的波動較大，在床層中的分布不均勻。

（5）二次流動區(qū)為環(huán)狀流動，流體和細長顆粒的“壁面效應(yīng)”顯著影響二次流動區(qū)的徑向固含率分布，呈現(xiàn)中心區(qū)域高壁面附近低的特點。

符 號 說 明

ar——長徑比 c——細長顆粒在溶液中的體積分數(shù)，% cs——固相體積分數(shù)，% D——洗滌冷卻室內(nèi)徑，m dout——下降管外徑，m dp——細長顆粒直徑，m k——湍動能，m2?s-2 L——水力直徑，m l——湍流長度尺度，m le——含能尺度，m lp——細長顆粒長度，m ug——環(huán)隙氣速，m?s-1 ε——湍動能耗散率，m2?s-3 εg——氣含率，% εs——固含率，% 下角標 g——氣相 s——固相 p——顆粒相

References

[1] 王亦飛, 于廣鎖, 龔欣, 等. 高溫氣體洗滌冷卻裝置: 103232863A[P]. 2014-06-18. WANG Y F, YU G S, GONG X,. High temperature gas cooling and washing equipment: 103232863A[P]. 2014-06-18.

[2] 盧瑞華, 王亦飛, 蘇宜豐, 等. 復合鼓泡床洗滌冷卻室中液滴夾帶統(tǒng)計模型[J]. 化學反應(yīng)工程與工藝, 2006, 22(1): 37-42. LU R H, WANG Y F, SU Y F,. Droplets statistics model for water entrainment in the separation zone of the compound bubbling chamber[J]. Chemical Reaction Engineering and Technology, 2006, 22(1): 37-42.

[3] 王亦飛, 陳意心, 劉霞, 等. 新型洗滌冷卻室內(nèi)的氣體帶液問題[J]. 化學反應(yīng)工程與工藝, 2008, 24(1): 24-28. WANG Y F, CHEN Y X, LIU X,. Entraining liquid problem in a new type of scrubbing-cooling chamber[J]. Chemical Reaction Engineering and Technology, 2008, 24(1): 24-28.

[4] 林嵐, 王亦飛, 郭強強, 等. 液相性質(zhì)對新型洗滌冷卻室內(nèi)液滴夾帶的影響[J]. 過程工程學報, 2011, 11(5): 742-746. LIN L, WANG Y F, GUO Q Q,. Effects of liquid properties on liquid entrainment in a new scrubbing-cooling chamber[J]. The Chinese Journal of Process Engineering, 2011, 11(5): 742-746.

[5] 郭強強, 王亦飛, 潘飛, 等. 復合型洗滌冷卻室內(nèi)液滴夾帶規(guī)律研究[J]. 中國電機工程學報, 2014, 34(11): 1769-1778. GUO Q Q, WANG Y F, PAN F,. Research on droplet entrainment in the composite quench chamber[J]. Proceedings of the CSEE, 2014, 34(11): 1769-1778.

[6] 王亦飛, 吳宏濤, 賀必云, 等. 新型旋流降膜式洗滌冷卻環(huán)的開發(fā)與研究[J]. 華東理工大學學報(自然科學版), 2005, 31(6): 752-755. WANG Y F, WU H T, HE B Y,. Development of a new rotary falling film quench ring for scrubbing and cooling hot syngas[J]. Journal of East China University of Science and Technology (Natural Science Edition), 2005, 31(6): 752-755.

[7] 王靈萍, 王亦飛, 郭強強, 等. 洗滌冷卻管內(nèi)垂直降膜流動特性[J]. 化工學報, 2013, 64(6): 1959-1968. WANG L P, WANG Y F, GUO Q Q,. Flow characteristics of vertical falling film in scrubbing-cooling pipe[J]. CIESC Journal, 2013, 64(6): 1959-1968.

[8] 顏留成, 王亦飛, 王靈萍, 等. 洗滌冷卻環(huán)內(nèi)液體冷態(tài)流動行為數(shù)值模擬[J]. 化工學報, 2014, 65(8): 2927-2933. YAN L C, WANG Y F, WANG L P,. Numerical simulation of cold state liquid flow in scrubbing-cooling ring[J]. CIESC Journal, 2014, 65(8): 2927-2933.

[9] YAN L C, WANG Y F, WU Z W,. Research of vertical falling film behavior in scrubbing-cooling tube[J]. Chem. Eng. Res. Des., 2017, 117: 627-636.

[10] 王晶, 王亦飛, 顏留成, 等. 管內(nèi)垂直下降液膜速度與厚度分布特性[J]. 化工學報, 2016, 67(6): 2239-2245. WANG J, WANG Y F, YAN L C,. Characteristics of velocity and thickness distribution of liquid film in vertical falling tube[J]. CIESC Journal, 2016, 67(6): 2239-2245.

[11] 柴雪琴, 王亦飛, 尹柯, 等. 流體誘導下洗滌冷卻室內(nèi)部構(gòu)件的振動特性[J]. 華東理工大學學報(自然科學版), 2010, 36(3): 346-350. CHAI X Q, WANG Y F, YIN K,. Vibrational property of internal component in scrubbing-cooling chamber with flow induced[J]. Journal of East China University of Science and Technology (Natural Science Edition), 2010, 36(3): 346-350.

[12] 陳程, 徐峰, 袁竹林, 等. 破泡結(jié)構(gòu)對洗滌室內(nèi)氣泡特性的影響研究[J]. 熱力發(fā)電, 2011, 40(7): 24-27. CHEN C, XU F, YUAN Z L,. Study on influence of bubble-breaking structure upon the gas bubble properties in the scrubbing chamber[J]. Thermal Power Generation, 2011, 40(7): 24-27.

[13] 付碧華. 洗滌冷卻室內(nèi)部結(jié)構(gòu)優(yōu)化及熱質(zhì)傳遞過程的數(shù)值模擬[D]. 上海: 華東理工大學, 2012. FU B H. Structure optimization of WSCC and numerical simulation of heat and mass transfer process[D]. Shanghai: East China University of Science and Technology, 2012.

[14] 賀必云. 新型洗滌冷卻室內(nèi)多相流動特征的研究[D]. 上海: 華東理工大學, 2005. HE B Y. Study on multiphase flow characteristics of the new type of scrubbing-cooling chamber[D]. Shanghai: East China University of Science and Technology, 2005.

[15] 吳晅. 洗滌冷卻室內(nèi)氣液固三相流動及熱質(zhì)傳遞規(guī)律研究[D]. 南京: 東南大學, 2009. WU X. Study on gas-liquid-solid flow and heat & mass transfer in scrubbing-cooling chamber[D]. Nanjing: Southeast University, 2009.

[16] 吳晅, 梁盼龍, 王麗芳, 等. 氣固兩相流穿越液池過程顆粒運動及分布特性[J]. 化工學報, 2015, 66(3): 905-914. WU X, LIANG P L, WANG L F,. Distribution properties and movement of particles in gas-solids flow passing through liquid bath[J]. CIESC Journal, 2015, 66(3): 905-914.

[17] 付碧華, 王亦飛, 林嵐, 等. 洗滌冷卻管出口處的射流深度和界面波動特性[J]. 化工學報, 2011, 62(7): 1817-1823. FU B H, WANG Y F, LIN L,. Penetration depth and characteristics of interface fluctuation at quenching pipe outlet[J]. CIESC Journal, 2011, 62(7): 1817-1823.

[18] 龔曉波, 顧兆林, 林高平, 等. 水煤漿氣化爐激冷流程中氣液兩相負浮力流動的數(shù)值模擬[J]. 化工學報, 2003, 54(7): 930-935. GONG X B, GU Z L, LIN G P,. Numerical simulation of gas-liquid two-phase flow with reverse buoyancy in quench chamber of coal gasifier and its application[J]. Journal of Chemical Industry and Engineering（China）, 2003, 54(7): 930-935.

[19] HUANG L Z, GAO X, LIN J Z. Cylindrical particulate internal flows: a review[J]. Fron. Mech. Eng., 2012, 7(4): 385-393.

[20] GORE R A, CROWE C T. Effect of particle size on modulating turbulent intensity[J]. Int. J. Multiphase Flow, 1989, 15(2): 279-285.

[21] 林建忠, 林江, 石興. 兩相流中柱狀固粒對流體湍動特性影響的研究[J]. 應(yīng)用數(shù)學和力學, 2002, 23(5): 483-488. LIN J Z, LIN J, SHI X. Research on the effect of cylinder particles on the turbulent properties in particulate flows[J]. Applied Mathematics and Mechanics, 2002, 23(5): 483-488.

[22] NASR-EL-DIN H A, MAC TAGGART R S, MASLIYAH J H. Local solids concentration measurement in a slurry mixing tank[J]. Chem. Eng. Sci., 1996, 51(8): 1209-1220.

[23] GANDHI B, PRAKASH A, BERGOUGNOU M A. Hydrodynamic behavior of slurry bubble column at high solids concentrations[J]. Powder Technol., 1999, 103(2): 80-94.

[24] 曹俊雅, 馬斌, 李向陽, 等. 測定氣-液-固攪拌槽內(nèi)相含率分布取樣法的改進[J]. 黃金科學技術(shù), 2015, 23(4): 80-85. CAO J Y, MA B, LI X Y,. The improved sampling method for the measurement of phase holdup distributions in gas-liquid-solid stirred tank[J]. Gold Science and Technology, 2015, 23(4): 80-85.

[25] CARLSSON A, LUNDELL F, S?DERBERG D. The wall effect on the orientation of fibres in a shear flow[J]. Ann. Trans. Nordic Rheol. Soc., 2006, 14: 83-89.

[26] HOLM R, S?DERBERG D. Shear influence on fibre orientation[J]. Rheol. Acta, 2007, 46(5): 721-729.

[27] KU X K, LIN J Z. Effect of two bounding walls on the rotational motion of a fiber in the simple shear flow[J]. Fiber. Polym., 2009, 10(3): 302-309.

[28] OZOLINS A, STRAUTINS U. Simple models for wall effect in fiber suspension flows[J]. Math. Model. Anal., 2014, 19(1): 75-84.

[29] OHNUKI A, AKIMOTO H. Experimental study on transition of flow pattern and phase distribution in upward air-water two-phase flow along a large vertical pipe[J]. Int. J. Multiphase Flow, 2000, 26(3):367-386.

[30] DINH S M. On the rheology of concentrated fiber suspensions[D]. Cambridge: Massachusetts Institute of Technology, 1981.

[31] FOLGAR F, TUCKER C L. Orientation behavior of fibers in concentrated suspensions[J]. J. Reinf. Plast. Comp., 1984, 3(2): 98-119.

[32] KUMAR P, RAMARAO B V. Enhancement of the sedimentation rates of fibrous suspensions[J]. Chem. Eng. Commun., 1991, 108(1): 381-401.

[33] HERZHAFT B, GUAZZELLI é. Experimental study of the sedimentation of dilute and semi-dilute suspensions of fibres[J]. J. Fluid Mech., 1999, 384(4): 133-158.

[34] 唐家鵬. FLUENT14.0超級學習手冊[M]. 北京: 人民郵電出版社, 2013: 184-185. TANG J P. FLUENT14.0 Super Learning Manual[M]. Beijing: Posts & Telecom Press, 2013: 184-185.

[35] 張兆順, 崔桂香, 許春曉. 湍流理論與模擬[M]. 北京: 清華大學出版社, 2005: 79. ZHANG Z S, CUI G X, XU C X. Theory and Modeling of Turbulence[M]. Beijing: Tsinghua University Press, 2005:79.

[36] MACKAPLOW M B, SHAQFEH E S G. A numerical study of the sedimentation of fibre suspensions[J]. J. Fluid Mech., 1998, 376: 149-182.

[37] HERZHAFT B, GUAZZELLI é, MACKAPLOW M B,. Experimental investigation of the sedimentation of a dilute fiber suspension[J]. Phys. Rev. Lett., 1996, 77(2): 290-293.

Slender particle-containing multi-phase distribution characteristics in scrubbing-cooling chamber

PENG Xin, WANG Yifei, WEI Zongyao, CHEN Futian, YU Guangsuo

(Key Laboratory of Coal Gasification and Energy Chemical Engineering of Ministry of Education, Institute of Clean Coal Technology, East China University of Science and Technology, Shanghai 200237, China)

In order to study multi-phase distribution characteristics of slender particles in scrubbing-cooling chamber, an improved direct sampling method was used to simultaneously measure locally axial and radial distribution of solid concentration and gas holdup in a cold model apparatus at various operating conditions. The cold model apparatus was scaled down according to geometric dimensions of scrubbing-cooling chamber of industrial gasification system. The results showed that scrubbing-cooling chamber could be divided into two zones of the upper gas-liquid-solid mixing zone and the lower solid-liquid flowing zone with boundary at cross-section of the downcomer outlet. The gas-liquid-solid mixing zone was consisted of downcomer outlet, bubble-break plate and gas cap regions, while the solid-liquid flow zone was consisted of gas phase turbulence, recirculation, and secondary flow regions. Due to compounded effects of particle hindrance (slowed down sedimentation) and agglomeration (accelerated sedimentation), the axial distribution of solid concentration exhibited wavy fluctuations. The increase of annular gas velocity, solid volume concentration and aspect ratio of length over diameter enhanced bed turbulence and promoted radial gas dissipation. Change of operating conditions altered particle drifting velocity and resulted in fluctuation of radial solid concentration distribution. Under the effects of gas turbulence and recirculation, the secondary flow region showed circular flow where the “wall effect” of fluids and particles forced solid concentration higher at center but lower near wall.

slender particles; scrubbing-cooling chamber; multi-phase flow; phase distribution; axial; radial; turbulence flow

10.11949/j.issn.0438-1157.20170480

TQ 021

A

0438—1157（2017）09—3368—12

2017-04-27收到初稿，2017-06-19收到修改稿。

王亦飛。

彭昕（1989—），男，博士研究生。

2017-04-27.

Prof.WANG Yifei, wangyf@ecust.edu.cn