脈沖噴動(dòng)微波真空干燥器內(nèi)噴動(dòng)流化動(dòng)力特性的數(shù)值研究

李 靜,金光遠(yuǎn),2，*,張 慜,2,崔政偉,仲建勛

(1.江南大學(xué)食品先進(jìn)制造裝備技術(shù)江蘇省重點(diǎn)實(shí)驗(yàn)室,江蘇無(wú)錫 214122;2.江南大學(xué)食品科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇無(wú)錫 214122)

脈沖噴動(dòng)微波真空干燥器內(nèi)噴動(dòng)流化動(dòng)力特性的數(shù)值研究

李 靜1,金光遠(yuǎn)1,2，*,張 慜1,2,崔政偉1,仲建勛1

(1.江南大學(xué)食品先進(jìn)制造裝備技術(shù)江蘇省重點(diǎn)實(shí)驗(yàn)室,江蘇無(wú)錫 214122;2.江南大學(xué)食品科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇無(wú)錫 214122)

為探索脈沖噴動(dòng)微波真空干燥系統(tǒng)(PSMVD)中氣固噴動(dòng)流化床的動(dòng)力特性,采用基于無(wú)粘性兩相流理論,使用計(jì)算流體力學(xué)歐拉-歐拉方法(CFD)對(duì)二維噴動(dòng)流化床中噴動(dòng)流化特性進(jìn)行了數(shù)值研究,對(duì)噴動(dòng)流化床的動(dòng)力特性和顆粒運(yùn)動(dòng)分布特征進(jìn)行了揭示,這些因素影響顆粒空間位置分布和其微波加熱的均勻性。結(jié)果表明,采用歐拉-歐拉雙流體模型,成功獲得脈沖噴動(dòng)微波真空干燥腔內(nèi)的二維噴動(dòng)流化床的噴動(dòng)流化動(dòng)力特性;噴動(dòng)流化過(guò)程在特定氣流噴動(dòng)速度下有一個(gè)初始階段和準(zhǔn)穩(wěn)態(tài)流化階段;噴動(dòng)速度是控制顆粒在噴動(dòng)床中顆粒位置的重要因素,存在一個(gè)影響流化模式轉(zhuǎn)變的關(guān)鍵噴動(dòng)速度;射流穿透深度和噴射速度之間存在線性關(guān)系。

真空微波干燥器,噴動(dòng)流化床,氣固兩相流,數(shù)值模擬

近年來(lái),微波真空干燥(MVD)作為一種潛在的獲得優(yōu)質(zhì)食品產(chǎn)品的干燥方法在食品工業(yè)進(jìn)行了初步的應(yīng)用,包括干果制品、蔬菜、谷物[1-2]。目前由于干燥腔內(nèi)電磁場(chǎng)的空間分布不均勻?qū)е聹囟葓?chǎng)分布不均勻,產(chǎn)品出現(xiàn)冷熱斑點(diǎn),這個(gè)問(wèn)題阻礙了這一技術(shù)的廣泛應(yīng)用。同時(shí)微波干燥不均勻的溫度分布也產(chǎn)生了食品微生物安全問(wèn)題[3-4]。微波真空干燥加熱均勻性受多種因素影響,例如真空腔結(jié)構(gòu)、物料空間位置、微波干燥腔內(nèi)的空間電磁場(chǎng)強(qiáng)度[1,5]。一些研究顯示,微波能量均布可以通過(guò)電磁場(chǎng)平均技術(shù)、機(jī)械運(yùn)動(dòng)方式或者氣動(dòng)攪拌來(lái)獲得[6-8]。

流態(tài)化過(guò)程為干燥床中的顆粒提供了一種氣體攪拌。由于顆粒表面不斷生成的新邊界層,有利于傳熱傳質(zhì)的進(jìn)行[9-10]。由于粗糙的食物顆粒如切塊或切片材料,尤其是當(dāng)其水分含量較高,表面較粘時(shí)難以通過(guò)傳統(tǒng)的流化床流化,脈沖噴動(dòng)床可用于流化顆粒較粗,應(yīng)用于微波真空干燥被認(rèn)為是解決加熱不均勻的問(wèn)題的一種有效手段。雖然一些研究人員已經(jīng)采用脈沖噴動(dòng)微波真空干燥(PSMVD)這個(gè)新技術(shù)對(duì)食品物料的干燥特性進(jìn)行實(shí)驗(yàn)研究[11-13],對(duì)于流化機(jī)理仍需要進(jìn)一步理論研究,以了解在脈沖噴動(dòng)微波真空干燥技術(shù)的干燥機(jī)理。

計(jì)算流體力學(xué)(CFD)方法被認(rèn)為是一個(gè)有效的工具,以獲得復(fù)雜的氣固兩相流的流體動(dòng)力學(xué)特性[14]。兩種不同的方法常用于模擬流化床或噴動(dòng)床的兩相流動(dòng),即歐拉-歐拉(EE)和歐拉-拉格朗日(EL)方法[15-17]。在本研究中,對(duì)于在脈動(dòng)噴動(dòng)微波干燥系統(tǒng)中簡(jiǎn)化二維噴動(dòng)流化床的噴動(dòng)流化過(guò)程,采用歐拉-歐拉雙流體模型方法進(jìn)行探究。對(duì)流化噴動(dòng)床的動(dòng)力特性和不同噴動(dòng)氣流速度下顆粒運(yùn)動(dòng)分布特征進(jìn)行探索,有助于明確PSMVD系統(tǒng)干燥的控制因素,改善顆粒空間微波加熱的均勻性。

1 材料與方法

1.1 材料與儀器

以來(lái)自中國(guó)寧波海通集團(tuán)農(nóng)場(chǎng)里的新鮮萵苣作為研究對(duì)象。萵苣洗凈去皮,測(cè)量獲得萵苣的物理屬性,采用FA1104電子天平測(cè)量質(zhì)量,采用E5061矢量網(wǎng)絡(luò)分析儀測(cè)量介電特性,水分含量采用國(guó)標(biāo)烘箱方法,所有測(cè)量過(guò)程進(jìn)行三次,取平均值[12-13],本次研究中萵苣質(zhì)量為200.0g。

圖1 脈沖噴動(dòng)微波真空干燥系統(tǒng)Fig.1 Schematic diagram of a laboratory pulsed spouted microwave-vacuum drying system注：1.進(jìn)料閥,2.板閥,3.微波加熱腔,4.微波發(fā)生器,5.循環(huán)水冷裝置,6.干燥腔,7.和13壓力表,8.氣固分離器,9.氣體濃縮器,10.真空泵,11.氣流電磁閥,12.氣流調(diào)節(jié)閥,14.氣源,15.控制面板,16.水管,17.氣體分布器,18.噴動(dòng)管,19.硅橡膠塞,20.物料,21.干燥室,22.固定裝置,23.光纖溫度傳感器。

實(shí)驗(yàn)室PSMVD實(shí)驗(yàn)設(shè)備有七個(gè)基本部分組成,其原理圖如圖1所示。干燥箱內(nèi)的壓力可以調(diào)節(jié)在3.5～100kPa的范圍。包括:(1)圓柱多模微波腔,配有四個(gè)2450MHz微波發(fā)生器;(2)一個(gè)圓形管道真空干燥箱;(3)一個(gè)脈沖噴動(dòng)系統(tǒng),配備了一套可調(diào)氣體流動(dòng)和分配單元以及一套1cm3/min空氣處理單元;(4)供熱系統(tǒng)。每個(gè)磁控管的輸出功率是在0.1到1.0kW由GPA-1.8kW微波功率控制器所控制;(5)水負(fù)載系統(tǒng),被用來(lái)防止磁控管過(guò)熱;(6)真空系統(tǒng)配有冷卻器和水循環(huán)真空泵,泵輸送率為1cm3/min。

1.2 實(shí)驗(yàn)方法

1.2.1 計(jì)算模型 本文采用歐拉-歐拉無(wú)粘性雙流體理論研究真空壓力條件下在噴動(dòng)流化床中氣固兩相流的動(dòng)力學(xué)特性。努森數(shù)是稀薄氣體的一個(gè)無(wú)量綱參數(shù)。由于脈沖噴動(dòng)微波真空干燥系統(tǒng)中,真空努森數(shù)(絕對(duì)壓力,7000Pa)遠(yuǎn)遠(yuǎn)低于1,符合流體連續(xù)介質(zhì)假設(shè),可以由納維爾斯托克斯方程(N-S方程)求解真空流動(dòng)[18]。在CFD(ANSYS FLUENT 13.0,America)軟件中對(duì)PSMVD設(shè)備中的氣固兩相流模擬。其中萵苣顆粒假設(shè)為球形,直徑為5mm。

在目前的二維流化床(2D噴動(dòng)床),噴動(dòng)流動(dòng)被簡(jiǎn)化為一個(gè)二維流,如圖2所示。干燥室的寬度是40.0mm,腔高度H為 516.0mm;空氣分布噴口的寬度L是8.0mm,空氣分配器角度是45°。空氣和萵苣顆粒分別被看做氣相和顆粒相。

圖2 二維噴動(dòng)流化床計(jì)算模型Fig.2 Computational models of air distribution structure

1.2.2 數(shù)值方法 氣固兩相流場(chǎng)采用ANSYS FLUENT求解器,通過(guò)求解雷諾數(shù)平均值的N-S方程來(lái)計(jì)算獲得。采用體積分?jǐn)?shù)的隱式方案,層流粘性模型,壓力速度耦合SIMPLEC算法,最小二乘法為基礎(chǔ)的一階迎風(fēng)空間離散,和一階隱式瞬態(tài)格式。噴射空氣和萵苣顆粒分別被設(shè)置為主相和次相。空氣和顆粒之間的相互作用被認(rèn)為是顆粒流動(dòng)的阻力和碰撞。采用結(jié)構(gòu)網(wǎng)格策略,利用ANSYS ICEM進(jìn)行網(wǎng)格剖分。網(wǎng)格包含370000個(gè)單元,單元的最小尺寸為0.01mm。進(jìn)行網(wǎng)格無(wú)關(guān)性計(jì)算,保證當(dāng)前網(wǎng)格尺寸的計(jì)算結(jié)果不依賴于網(wǎng)格疏密。

本文模擬邊界條件定義如下:入口被指定為速度入口,并且指定為主相速度;出口被指定為壓力出口,重力和顆粒體積分?jǐn)?shù)被指定。瞬時(shí)時(shí)間步長(zhǎng)設(shè)置0.001秒。為了考察噴動(dòng)速度對(duì)流化特征的影響,選用若干氣流噴動(dòng)速度工況進(jìn)行模擬(氣流噴動(dòng)速度 uspout=3.5、7、14、28、56m/s)。

圖3 最小噴動(dòng)流化速度下顆粒體積分布Fig.3 Particle volume fraction at the minimum spout fluidization velocity注:uspout=3.5m/s,t=0.2～1.0s。

2 結(jié)果與討論

2.1 最小流化速度下噴動(dòng)流化特征

當(dāng)前歐拉模型下萵苣顆粒發(fā)生流化時(shí)的最小流化速度為3.5m/s,其顆粒相體積分?jǐn)?shù)和速度如圖3和圖4所示。由圖3和圖4可以看出,萵苣顆粒在氣流噴動(dòng)1.0s后保持傳統(tǒng)的固定床狀態(tài)。在噴動(dòng)流化床整體流態(tài)化過(guò)程包括初始階段和穩(wěn)定階段。初始階段是射流穿透床層的一個(gè)瞬態(tài)發(fā)展過(guò)程,其中,所述射流最初出現(xiàn),然后無(wú)干擾的周期性地分裂成氣泡的第二階段,是一個(gè)穩(wěn)定的循環(huán)過(guò)程,其中顆粒循環(huán)在一個(gè)固定的模式。在離噴射入口高度為29.0mm處形成懸浮顆粒界面,大多數(shù)顆粒占據(jù)空氣分布器底部,很少的顆粒能夠逃離這個(gè)界面,此時(shí)最大顆粒速度達(dá)0.1m/s。

圖4 最小噴動(dòng)流化速度下顆粒速度分布 Fig.4 Particle velocity at the minimum spout fluidization velocity注:m/s,uspout=3.5m/s,t=1.0～3.0s。

2.2 流動(dòng)穩(wěn)定性

隨著噴口處氣體流速的增加,噴動(dòng)流化床的氣固兩相流體動(dòng)力學(xué)特性開始改變。微波干燥腔內(nèi)兩相流動(dòng)穩(wěn)定性如圖5所示,給出了在不同氣體噴動(dòng)速度下,位于初始靜止時(shí)顆粒表面中間點(diǎn)處的靜壓力脈動(dòng)特征和萵苣顆粒的速度脈動(dòng)特征;整個(gè)噴動(dòng)流化有一個(gè)明確的約1.0秒的不穩(wěn)定初始過(guò)程,在這種初始波動(dòng)階段后,當(dāng)氣體噴射速度不高于21m/s時(shí)趨于穩(wěn)定。 隨噴動(dòng)速度增加,這種不穩(wěn)定狀態(tài)的持續(xù)時(shí)間增加,且靜壓值趨于降低,而靜壓波動(dòng)幅值趨于增大。在相同位置萵苣顆粒速度如圖6所示,顆粒速度在初始階段差異較大,當(dāng)噴動(dòng)速度不高于21m/s時(shí)穩(wěn)定。 當(dāng)噴動(dòng)速度達(dá)到28m/s,顆粒速度顯著增加,特點(diǎn)是周期性波動(dòng),并且表現(xiàn)為振動(dòng)的流態(tài)化。上述分析表明,噴動(dòng)速度是影響噴動(dòng)流化特征的一個(gè)重要因素,存在一個(gè)關(guān)鍵噴動(dòng)速度影響噴動(dòng)流化模式轉(zhuǎn)變;整個(gè)流化過(guò)程可分為兩個(gè)階段:初始階段和準(zhǔn)穩(wěn)態(tài)流化階段。

圖5 不同噴動(dòng)速度下干燥腔內(nèi)靜壓脈動(dòng)Fig.5 Fluctuation of static pressure at different spout gas velocity

圖6 不同噴動(dòng)速度下干燥腔內(nèi)粒子脈動(dòng)特征Fig.6 Fluctuation of particle velocity at different spout gas velocity

2.3 顆粒流的發(fā)展過(guò)程

當(dāng)噴口氣體速度達(dá)28m/s,考慮到顆粒流的波動(dòng)特性,分別討論在噴動(dòng)流速為14m/s和28m/s的情形來(lái)明確影響流動(dòng)模式轉(zhuǎn)變的關(guān)鍵噴動(dòng)速度。

圖7 粒子體積分?jǐn)?shù)分布Fig.7 Particle volume fraction注:uspout=14m/s,t=0.1～1.0s。

圖9 粒子體積分?jǐn)?shù)分布Fig.9 Particle volume fraction distribution注:uspout=28m/s,t=0.1～1.0s。

圖7和圖8顯示了噴動(dòng)速度為14m/s,噴動(dòng)時(shí)間為0.1～3.0s內(nèi)的瞬時(shí)顆粒流動(dòng)特征,顆粒流動(dòng)初始階段t=0.1～1.0s秒的粒子分布情況,如圖7所示,進(jìn)入準(zhǔn)穩(wěn)態(tài)流化階段,粒子運(yùn)動(dòng)特征如圖8所示。與噴動(dòng)速度3.5m/s相比,初始階段射流產(chǎn)生了更多、更大的氣泡,導(dǎo)致該區(qū)域更頻繁的顆粒運(yùn)動(dòng)。隨著射流發(fā)展,大多數(shù)顆粒在噴動(dòng)氣流攜帶作用下獲得提升,并形成一個(gè)球形分布區(qū)域。但仍有一小部分顆粒分布在空氣分配器壁的附近,顆粒積聚發(fā)生在空氣分配器底部。噴動(dòng)時(shí)間間隔0.5s從1.0s到3.0s,顆粒速度和最大提升高度分布結(jié)果如圖8。當(dāng)顆粒流動(dòng)模式處于一個(gè)穩(wěn)定的狀態(tài),在高度45mm噴口上方形成顆粒懸浮界面。最高的顆粒速度位于球形分布區(qū)域,這表明在此工況下噴動(dòng)流化床有兩條顆粒循環(huán)軌跡。少數(shù)顆粒從懸浮表面逸出,顆粒速度的最大值達(dá)0.3m/s(標(biāo)尺未給出)。

如圖9和圖10所示是噴動(dòng)速度在28m/s的瞬時(shí)顆粒流動(dòng)模式。瞬時(shí)流動(dòng)模式在1.0s的初始階段時(shí),帶有更多更大氣泡的射流迅速出現(xiàn),穿透顆粒床層,并與大多數(shù)顆粒顯著升高了。隨噴動(dòng)射流發(fā)展,積累在空氣分配器器壁的顆粒被攜帶到上部區(qū)域。在這種情況下,噴動(dòng)射流引起大多數(shù)顆粒在橫向和軸向方向的運(yùn)動(dòng)。顆粒懸浮高度位于噴口上方150到250mm,顆粒流動(dòng)具有一定的波動(dòng)特征,噴口處顆粒流動(dòng)是不穩(wěn)定的(圖10)。 與噴動(dòng)速度為14m/s的顆粒流相比,顆粒速度的最大值增大至1.6m/s(標(biāo)尺未給出),更高速度顆粒的分布位置移至20～100mm高度區(qū)域。存在一個(gè)含有較低速度顆粒的狹窄區(qū)域,這可能是由主噴口氣體流動(dòng)引起的。

在不同的噴口速度流動(dòng)模式的特征表明,噴動(dòng)氣流速度和顆粒狀態(tài)密切相關(guān),存在一個(gè)關(guān)鍵噴動(dòng)速度,影響顆粒流動(dòng)模式的轉(zhuǎn)變。

圖8 隨噴動(dòng)時(shí)間發(fā)展的粒子速度分布(uspout=14m/s)Fig.8 Particle velocity distribution at different spout time(uspout=14m/s)

圖10 隨噴動(dòng)時(shí)間發(fā)展的速度分布(uspout=28m/s)Fig.10 Particle velocity distribution at different spout time(uspout=28m/s)

2.4 噴射氣流的穿透特性

在噴動(dòng)流化床中,無(wú)論氣體以氣泡的形式離開噴口,還是以脈動(dòng)噴射或連續(xù)噴射的形式離開,噴射區(qū)域存在大量的氣體和固體混合[19]。射流穿透深度來(lái)表征噴口氣流的關(guān)鍵參數(shù),它也隨著噴動(dòng)流化床中噴射氣體的速度提高而提高。用來(lái)預(yù)測(cè)射流穿透深度的不同模型之間的差異主要是由于射流邊界的不同定義[20]。在本文中,射流穿透深度定義為孔隙率為0.8的等值線[21]。噴動(dòng)氣體速度對(duì)穿透深度和壓力降影響如圖11所示。在低于28m/s的噴動(dòng)氣流速度下穿透深度增加緩慢,直到56m/s才迅速增加。在射流穿透深度和氣體噴動(dòng)速度之間存在一個(gè)近似線性關(guān)系。

床層壓降是表征真空流化床干燥器顆粒循環(huán)模式的另一個(gè)重要的參數(shù),在準(zhǔn)穩(wěn)態(tài)階段,定義氣固混合物進(jìn)口和出口之間壓差為壓降(ΔP)。圖11所示床壓降表明,當(dāng)管口的氣體速度從14m/s增加到21m/s,該壓降從負(fù)壓變?yōu)檎龎骸＿@種現(xiàn)象表明粒顆粒循環(huán)模式的一個(gè)重要變化,噴口的氣流作用于顆粒的提升力克服了顆粒的重力,顆粒開始被從空氣分配器底部提升到真空腔室的上部區(qū)域。因此,從14m/s至28m/s的速度范圍是一個(gè)影響顆粒循環(huán)模式的關(guān)鍵噴動(dòng)速度區(qū)域。

圖11 噴動(dòng)流化穿透深度和壓降Fig.11 Penetration depth and Pressure drop for a 2D vacuum cavity

3 結(jié)論

本文成功采用歐拉-歐拉雙流體模型,獲得脈沖噴動(dòng)微波真空干燥腔內(nèi)的二維噴動(dòng)流化床的噴動(dòng)流化動(dòng)力特性。基于模擬結(jié)果,分析噴動(dòng)流化床的動(dòng)態(tài)特性和在不同噴動(dòng)速度工況下噴動(dòng)氣流對(duì)顆粒運(yùn)動(dòng)模型的影響規(guī)律。在一定噴動(dòng)速度下,流化過(guò)程可分為兩個(gè)階段:初始階段和準(zhǔn)穩(wěn)態(tài)流化階段。噴口速度是控制顆粒在噴動(dòng)流化床中位置的重要因素。存在一個(gè)關(guān)鍵噴動(dòng)速度區(qū)間,在此范圍內(nèi),顆粒流動(dòng)模式由穩(wěn)態(tài)轉(zhuǎn)化為非穩(wěn)態(tài)流動(dòng)。射流穿透深度與氣體噴動(dòng)速度之間存在一個(gè)近似的線性關(guān)系。在準(zhǔn)穩(wěn)態(tài)階段,隨著噴口速度增加顆粒循環(huán)模式明顯改變,粒子開始被提升到真空干燥器上部區(qū)域。

[1]Zhang M,Jiang H,Lim RX. Recent developments in microwave-assisted drying of vegetables,fruits,and aquatic products-Drying kinetics and quality considerations[J]. Drying Technology,2010,28(11):1307-1316.

[2]Li ZY,Wang RF,Kudra T. Uniformity issue in microwave drying[J]. Drying Technology,2011,29(6):652-660.

[3]Vadivambal R,Jayas DS. Non-uniform temperature distribution during microwave heating of food materials-A review[J]. Food and Bioprocess Technology,2010,3:161-171.

[4]Jangam SV. An overview of recent developments and some R&D challenges related to drying of foods[J]. Drying Technology,2011,29(12):1343-1357.

[5]Zhang M,Tang J,Mujumdar AS,et al. Trends in microwave-related drying of fruits and vegetables[J]. Trends in Food Science & Technology,2006,17(10):524-534.

[6]Torringa HM,Dijk EJV,Bartels PV. Microwave puffing of vegetables:modeling and measurements[J]. In:Proceedings of 31st microwave power symposium,1996,16-19.

[7]Balakrishnan M,Raghavan GSV,Sreenarayanan VV,et al. Batch drying kinetics of cardamom in a two-dimensional spouted bed[J]. Drying Technology,2011,29(11):1283-1290.

[8]Feng H,Tang J. Microwave finish drying of diced apples in a spouted bed[J]. Journal of Food Science,1998:63(4),679-683.

[9]Feng H,Tang J,Cavalieri RP,Plumb OA. Heat and mass transport in microwave drying of porous materials in a spouted bed[J]. AIChE Journal,2001,47(7):1499-1512.

[10]Jambhale AS,Barbadekar BV. Microwave drying system with high-tech phase controller:a modified applicator[J]. World Acad Sci Eng Technol,2008,46:1-5.

[11]Wang Y,Zhang M,Mujumdar AS,et al. Study of drying uniformity in pulse-spouted microwave-vacuum drying of stem lettuce slices with regard to product quality[J]. Drying Technology,2013,31(1):91-101.

[12]Wang Y,Zhang M,Mujumdar AS,et al. Microwave-assisted pulse-spouted bed freeze-drying of stem lettuce slices-effect on product quality[J]. Food and Bioprocess Tech. 2013,6:3530-3543.

[13]Wang Y,Zhang M,Mujumdar AS,et al. Experimental investigation and mechanism analysis on microwave freeze drying of stem lettuce cubes in a circular conduit[J]. Drying Technology,2012,30(11-12):1377-1386.

[14]Rahimi M R,Azizi N,Hosseini SH,et al. CFD study of hydrodynamics behavior of a vibrating fluidized bedusing kinetic-frictional stress model of granular flow[J]. Korean Journal of Chem. Eng.,2013,30(3):761-770.

[15]Pai M,Subramaniam SA. Comprehensive probability density function formalism for multiphase flows[J]. Journal of Fluid Mechanics,2009,628:181-228.

[16]Bouillard J,Lyczkowski R,Gidaspow D. Porosity distributions in a fluidized bed with an immersed obstacle[J]. AIChE Journal,1989,35:908-922.

[17]Pritchett J,Blake T,Garg S. A numerical model of gas fluidized beds[C]. AIChE Symp.,1978,176:134-148.

[18]Shen C. Rarefied gas dynamics:fundamentals,simulations and micro flows[B]. Publisher:Springer-Verlag Berlin and Heidelberg GmbH & Co. K. 2005.

[19]Hong R,Guo Q,Luo G,et al. On the jet penetration height in fluidized beds with two vertical jets[J]. Powder Technology,2003,133(1):216-227.

[20]Musmarra D. Influence of particle size and density on the jet penetration length in gas fluidized beds[J]. Industrial & Engineering Chemistry Research,2000,39(7):2612-2617.

[21]Gidaspow D,Ettehadleh B. Fluidization in two-dimensional beds with a jet hydrodynamic modeling[J]. Industrial & Engineering Chemistry Fundamentals,1983,22(2):193-201.

Numerical study on spout fluidization dynamics in pulsed-spout microwave-vacuum dryer

LI Jing1,JIN Guang-yuan1,2，*,ZHANG Min1,2,CUI Zheng-wei1,ZHONG Jian-xun1

(1.Jiangsu Province Key Laboratory of Advanced Food Manufacturing Equipment and Technology,Jiangnan University,Wuxi 214122,China;2. China State Key Laboratory of Food Science and Technology,Jiangnan University,Wuxi 214122,China)

The dynamic characteristics of gas-particle spout fluidized bed in a microwave-vacuum drying system(PSMVD)were investigated. The spout fluidization process in a Pseudo-2-D spouting bed was monitored by computational fluid dynamics(CFD)using the inviscid Two-fluid theory based on Euler-Euler method. The dynamic characteristics of the spout fluidized bed and the effect of spout elevation on the particle movement were revealed,which could be used to improve the uniformity of particle mixing and microwave heating. The results revealed that mathematical model succeeded to demonstrate the dynamic characteristics of gas-particle spout fluidized bed in a microwave-vacuum drying system;the spout fluidization process includes isolated,merged and transitional jets and the fluidization at specific spout gas velocity had a start-up stage and a quasi-steady fluidization stage. The spout velocity was an important factor controlling particle status in the spout fluidized bed and a critical velocity was identified for effect transition of the flow pattern. There was an approximately linear correlation between the jet penetration depth and the spout velocity.

microwave-vacuum dryer;spout fluidized bed;two phase flow;numerical simulation

2014-05-05

李靜(1972-),碩士,講師,主要從事食品加工裝備和無(wú)損檢測(cè)問(wèn)題研究。

*通訊作者:金光遠(yuǎn)(1979-)，女，博士，講師，研究方向：食品加工裝備。

江蘇省博士后基金項(xiàng)目,江蘇省青年基金項(xiàng)目(BK20130150)。

TS255-35

A

1002-0306(2015)11-0079-05

10.13386/j.issn1002-0306.2015.11.008