鼓泡床模擬概要

1、稠密氣固兩相流動(dòng)模擬鼓泡床流動(dòng)計(jì)算討論2012/4/18目錄1. 歐拉-歐拉方法簡(jiǎn)介2. 計(jì)算條件3. 計(jì)算結(jié)果 3.1 有中心射流條件下鼓泡床計(jì)算結(jié)果 3.2 無(wú)中心射流條件下鼓泡床計(jì)算結(jié)果 3.3 定常流動(dòng)計(jì)算結(jié)果4. 分析總結(jié) 歐拉-歐拉方法簡(jiǎn)介其基本假設(shè)是：1.流場(chǎng)中的每一位置，顆粒相與流體相共存并且相互滲透，且每一相具有各自的速度、溫度。2.每一顆粒相在空間中具有連續(xù)的速度溫度和容積分?jǐn)?shù)的分布 3.顆粒相除與流體有質(zhì)量和能量的相互作用外，還具有自身的湍流脈動(dòng)4.對(duì)稠密顆粒懸浮體，顆粒碰撞會(huì)引起附加的顆粒粘性、擴(kuò)散和熱傳導(dǎo)。 歐拉-歐拉方法（Eulerian Method）是將顆?？闯?/p>

2、是Euler坐標(biāo)系中的擬流體，具備與流體相似的動(dòng)力學(xué)特性，且與流體之間存在大滑移和顆粒擴(kuò)散，并能充分考慮顆粒質(zhì)量、動(dòng)量及能量的湍流擴(kuò)散，在不考慮相間質(zhì)量交換的情況下，通過(guò)求解兩相N-S方程進(jìn)行計(jì)算。計(jì)算條件 2D鼓泡床計(jì)算幾何模型實(shí)測(cè)的顆粒體積分?jǐn)?shù)（t=0.4 s，右圖為距中心線4mm處的顆粒體積分?jǐn)?shù)軸向分布） 邊界條件：無(wú)滑移； 氣相按層流計(jì)算； 所有變量的迭代殘差10-4，保證得到收斂解； 網(wǎng)格：結(jié)構(gòu)化網(wǎng)格；數(shù)量和粗細(xì)自定，在參數(shù)變化劇烈處加密； 初場(chǎng)設(shè)定：用patch（補(bǔ)丁）命令預(yù)先給定氣固兩相的速度值（包括軸向和徑向分量）。yx 多相流模型(Multiphase Model): Eul

3、erian 第一相 (Primary-phase): air 第二相 (Secondary-phase): glassbead參數(shù)英文名稱選取結(jié)果粘性 Granular Viscositysyamlal-obrien體積粘性 Granular Bulk Viscosityun-et al摩擦粘性 Frictional Viscosityschaeffer內(nèi)摩擦角 Angle of Internal Friction30deg摩擦壓力 Frictional Pressurebased ktgf摩擦模數(shù) Frictional Modulusderived徑向分布 Radial Distributi

4、onlun et al彈性模量Elasticity Modulusderived堆積密度極限 Packing Limit0.63顆粒擬流體模型參數(shù)選取計(jì)算條件計(jì)算結(jié)果有中心射流條件下鼓泡床計(jì)算結(jié)果 (計(jì)算時(shí)間1.0s)0s0.1s0.2s0.3s0.4s0.5s0.6s0.7s計(jì)算結(jié)果有中心射流條件下鼓泡床計(jì)算結(jié)果0.8s0.9s1.0s 有中心射流時(shí)，可以從計(jì)算結(jié)果中清晰地觀察到氣泡的生成、長(zhǎng)大以及上升至床面并破裂的過(guò)程，而且氣泡的尾渦也十分明顯。床面由清晰可見(jiàn)變?yōu)閯×易儎?dòng)，氣泡破裂過(guò)程中發(fā)生揚(yáng)析和夾帶。 在1s時(shí)間內(nèi)，先后有34個(gè)氣泡生成，氣泡生成的頻率約為4Hz。 由于模擬計(jì)算的理想化，

5、計(jì)算結(jié)果呈現(xiàn)高度對(duì)稱性，這與實(shí)際結(jié)果是不相符的。實(shí)際過(guò)程中任何微小的擾動(dòng)都可能消除這種對(duì)稱性。計(jì)算結(jié)果有中心射流條件下鼓泡床計(jì)算結(jié)果 選取x=0.04m，0.14m和0.28m三個(gè)特征線，考察其靜壓分布隨時(shí)間的變化。t=0.2st=0.4st=0.6st=0.8s計(jì)算結(jié)果有中心射流條件下鼓泡床計(jì)算結(jié)果 因?yàn)椴捎脡毫Τ隹谶吔鐥l件，因此出口壓力一致。在鼓泡過(guò)程中，床層壓差隨之變化。值得注意的是，當(dāng)氣泡剛剛生成時(shí)，氣泡附近的壓力大于氣泡遠(yuǎn)處的壓力，顆粒有向外擴(kuò)散的趨勢(shì)；而當(dāng)氣泡已經(jīng)生成時(shí)，氣泡附近的壓力小于氣泡遠(yuǎn)處的壓力，顆粒有向內(nèi)填充的趨勢(shì)。這一結(jié)論已通過(guò)其他時(shí)間的結(jié)果得到驗(yàn)證。 計(jì)算得到的顆粒濃

6、度分布與實(shí)驗(yàn)值相差極大。我覺(jué)得實(shí)驗(yàn)值存在問(wèn)題，因?yàn)樵趖=0.4s時(shí)，剛剛形成一個(gè)大的氣泡，顆粒濃度怎么會(huì)在氣泡的位置會(huì)有一個(gè)很大的提升？0.00.10.20.30.40.50.60.00.20.40.60.81.0 Glassbead Volume Fraction Kuipers, 1990 CalculationY Position (m)無(wú)中心射流條件下鼓泡床計(jì)算結(jié)果計(jì)算結(jié)果0s0.2s0.4s0.6s0.8s1.0sug=1m/sug=2m/sug=5m/s無(wú)中心射流條件下鼓泡床計(jì)算結(jié)果計(jì)算結(jié)果 選取x=0.14m的特征線，考察不同平均氣速下顆粒濃度分布隨時(shí)間的變化。0.00.20.4

7、0.60.81.00.00.10.20.30.40.50.60.7 Glassbead Volume FractionY Position (m) ug=1m/s ug=2m/s ug=5m/s0.00.20.40.60.81.00.00.10.20.30.40.50.60.7 Glassbead Volume FractionY Position (m) ug=1m/s ug=2m/s ug=5m/s0.00.20.40.60.81.00.00.10.20.30.40.50.60.7 Glassbead Volume FractionY Position (m) ug=1m/s ug=2m/

8、s ug=5m/s0.00.20.40.60.81.00.00.10.20.30.40.50.60.7 Glassbead Volume FractionY Position (m) ug=1m/s ug=2m/s ug=5m/st=0.2st=0.4st=0.6st=0.8s無(wú)中心射流條件下鼓泡床計(jì)算結(jié)果計(jì)算結(jié)果 對(duì)于同一平均氣速，隨著時(shí)間的增加，床面逐漸升高，顆?；旌霞觿。w粒濃度從單峰分布轉(zhuǎn)變?yōu)槎喾宸植肌?隨著平均氣速的提高，床面升高加速， 對(duì)于氣速為2m/s和5m/s的工況，大量的顆粒被吹出床體，尤其是后者，顆粒幾乎被吹空。 與有中心射流條件下的工況類似，由于模擬計(jì)算的理想化，計(jì)算結(jié)果

9、呈現(xiàn)高度對(duì)稱性，這與實(shí)際結(jié)果是不相符的。0.00.20.40.60.81.00.00.10.20.30.40.50.60.7 Glassbead Volume FractionY Position (m) ug=1m/s ug=2m/s ug=5m/st=1.0s定常流動(dòng)計(jì)算結(jié)果計(jì)算結(jié)果 定常流動(dòng)條件下，計(jì)算中無(wú)法實(shí)現(xiàn)收斂，殘差始終維持在10-1上下。不同迭代步數(shù)下，得到的流場(chǎng)結(jié)果也不相同。這說(shuō)明該條件下鼓泡床的氣固兩相流是無(wú)法收斂到一個(gè)定常的結(jié)果的。定常流動(dòng)計(jì)算部分結(jié)果云圖分析總結(jié)1.利用Eulerian兩相流計(jì)算模型，對(duì)鼓泡床內(nèi)氣固流動(dòng)進(jìn)行了初步的模擬，分別得到了有中心射流和無(wú)中心射流工況下的床內(nèi)氣固分布、壓力分