陳 帥, 周勇軍, 林偉振, 史亦航, 陶蘭蘭

改進(jìn)型三斜葉-Rushton組合槳槽內(nèi)流場(chǎng)特性研究

陳 帥1, 周勇軍1, 林偉振1, 史亦航1, 陶蘭蘭2

(1. 南京工業(yè)大學(xué) 機(jī)械與動(dòng)力工程學(xué)院, 江蘇 南京 211816;2. 江蘇省特種設(shè)備安全監(jiān)督檢驗(yàn)研究院國(guó)家化工設(shè)備質(zhì)量監(jiān)督檢驗(yàn)中心, 江蘇 蘇州 215600)

為將改進(jìn)型三斜葉-Rushton組合槳高效應(yīng)用于實(shí)際工程，通過(guò)計(jì)算流體動(dòng)力學(xué)(CFD)對(duì)裝配改進(jìn)型三斜葉-Rushton組合槳攪拌槽內(nèi)的流場(chǎng)進(jìn)行了模擬研究，與三斜葉-Rushton組合槳進(jìn)行了對(duì)比；并通過(guò)粒子圖像測(cè)速(PIV)技術(shù)研究該組合槳槽內(nèi)攪拌轉(zhuǎn)速、槳葉間距2及離底距離1的變化對(duì)流場(chǎng)產(chǎn)生的影響。結(jié)果表明當(dāng)轉(zhuǎn)速為100 r×min-1時(shí)能有效改善槽底和上槳葉區(qū)域的混合效果，槽內(nèi)速度分布相對(duì)均勻，實(shí)驗(yàn)結(jié)果和模擬結(jié)果吻合效果較好；模擬結(jié)果表明改進(jìn)型三斜葉槳的泵吸效果較三斜葉槳更強(qiáng)，能有效縮小攪拌槽底部和內(nèi)壁附近的低速區(qū)域，槽內(nèi)混合效果更好；實(shí)驗(yàn)結(jié)果表明當(dāng)槳葉間距為0.25時(shí)，槽內(nèi)速度分布更合理；當(dāng)離底距離為0.3時(shí)，上槳葉區(qū)的較低速區(qū)域進(jìn)一步縮小。

攪拌槽；改進(jìn)型三斜葉-Rushton組合槳；流場(chǎng)；計(jì)算流體動(dòng)力學(xué)；粒子圖像測(cè)速；混合

1 前 言

進(jìn)行攪拌反應(yīng)時(shí)，釜內(nèi)的物料都伴隨著動(dòng)量、質(zhì)量、熱量的傳遞，攪拌器作為攪拌反應(yīng)過(guò)程中的核心部件直接影響具有攪拌流程產(chǎn)品的效率和質(zhì)量，攪拌器的研究對(duì)于降低生產(chǎn)成本、節(jié)約能耗十分重要[1]，因此研究特定攪拌器的攪拌釜內(nèi)流場(chǎng)特性對(duì)于工程應(yīng)用具有重要意義[2-4]。

攪拌釜內(nèi)流場(chǎng)的研究常用計(jì)算流體動(dòng)力學(xué)(computational fluid dynamics，CFD)模擬得到不同工況下的流場(chǎng)特性[5-8]，通過(guò)粒子圖像測(cè)速(particle image velocimetry，PIV)技術(shù)完成對(duì)模擬結(jié)果的驗(yàn)證[9-12]。為了進(jìn)一步實(shí)現(xiàn)對(duì)流場(chǎng)時(shí)間和空間的精確解析，體三維速度測(cè)量(volumetric three-component velocimetry measurements，V3V)技術(shù)[13-14]也開始應(yīng)用到流場(chǎng)研究當(dāng)中。謝明輝等[15]將不同濃度下的羧甲基纖維素鈉(CMC)溶液代替菌絲發(fā)酵液，對(duì)三寬葉翼型軸流槳及其組合槳進(jìn)行了功耗和體積氧傳質(zhì)效率研究，發(fā)現(xiàn)裝有HY+2KCXu組合槳(上兩層三寬葉上翻式翼型KCXu槳，底層六葉半圓管圓盤渦輪HY槳)的攪拌罐內(nèi)的菌球相對(duì)于其他組合槳更容易形成。在裝配有上提式和下壓式軸流式槳及其組合槳的攪拌槽中，朱姝等[16]采用雙電導(dǎo)探針?lè)ê蛿?shù)值解析重構(gòu)法對(duì)攪拌槽內(nèi)徑向和軸向的氣泡分布進(jìn)行了測(cè)定，發(fā)現(xiàn)上提式的平均氣泡尺寸較下壓式大13%。楊娟等[17]通過(guò)攪拌混合實(shí)驗(yàn)比較了三斜葉槳及其組合槳在非牛頓流體攪拌槽內(nèi)的功率準(zhǔn)數(shù)，發(fā)現(xiàn)在相同單位功耗下，三斜葉槳與Rushton組合槳的混合性能相對(duì)于其他單層槳的混合效率較高。周勇軍等[18]利用PIV技術(shù)對(duì)二斜葉框式組合槳攪拌槽內(nèi)的流動(dòng)特性進(jìn)行了研究，發(fā)現(xiàn)隨著槳葉離底距離的增大，并不利于兩槳葉之間的混合。吳高杰等[19]采用實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方式，研究了槳葉安裝間距、介質(zhì)黏度等因素對(duì)攪拌過(guò)程中氧傳遞的不同影響，發(fā)現(xiàn)渦輪斜葉翼型組合槳使得攪拌釜內(nèi)的氣體分布更加均勻。宮磊等[20]通過(guò)數(shù)值模擬的方式將六斜葉槳和渦輪槳進(jìn)行組合，研究了4種攪拌組合槳分別對(duì)發(fā)酵罐內(nèi)的流場(chǎng)和攪拌混合時(shí)間產(chǎn)生的影響，得到了最佳的組合槳葉為六斜葉和渦輪槳的組合槳。劉寶慶等[21]通過(guò)數(shù)值模擬的方法研究了中高黏度流體下攪拌釜內(nèi)的流場(chǎng)特性、混合特性和功率特性，發(fā)現(xiàn)六斜葉-六直葉組合槳的功耗較低，并且混合效果較好。周勇軍等[22]基于CFD對(duì)改進(jìn)的Intermig組合槳進(jìn)行數(shù)值模擬研究，并用PIV實(shí)驗(yàn)驗(yàn)證發(fā)現(xiàn)槳葉間距和離底距對(duì)槳葉周圍的局部流場(chǎng)有明顯的影響。

文獻(xiàn)發(fā)現(xiàn)相關(guān)傳統(tǒng)斜葉組合槳具備較好的攪拌效果，改進(jìn)型三斜葉槳是在傳統(tǒng)三斜葉槳基礎(chǔ)上進(jìn)行優(yōu)化得到的新型槳，目前國(guó)內(nèi)外對(duì)改進(jìn)型三斜葉-Rushton組合槳攪拌槽內(nèi)流場(chǎng)的研究還未見報(bào)道。本研究利用CFD模擬和PIV技術(shù)對(duì)改進(jìn)型三斜葉-Rushton組合槳槽內(nèi)流場(chǎng)特性進(jìn)行研究，分析不同轉(zhuǎn)速、槳葉間距及離底距離下改進(jìn)型三斜葉-Rushton組合槳攪拌槽內(nèi)流場(chǎng)特性，并將該組合槳與三斜葉-Rushton組合槳槽內(nèi)流場(chǎng)進(jìn)行對(duì)比分析。以期研究結(jié)果為改進(jìn)型三斜葉-Rushton組合槳的工程應(yīng)用提供參考。

2 實(shí)驗(yàn)和數(shù)值模擬方法

2.1 實(shí)驗(yàn)裝置及方法

攪拌槽及改進(jìn)型三斜葉槳結(jié)構(gòu)如圖1所示，其中攪拌槽為圓柱形透明玻璃槽和橢圓形封頭組成，改進(jìn)型三斜葉槳的特點(diǎn)是葉片部分彎折。攪拌槽壁厚1=10 mm，液位高度=605 mm，改進(jìn)型三斜葉槳和Rushton槳的直徑=212 mm，葉片厚度2=2 mm，折葉彎折角度=160°，主要尺寸如表1所示。

圖1 攪拌槽及改進(jìn)型三斜葉槳結(jié)構(gòu)尺寸圖

表1 主要結(jié)構(gòu)參數(shù)

本實(shí)驗(yàn)利用流場(chǎng)微觀性能測(cè)試平臺(tái)對(duì)改進(jìn)型三斜葉-Rushton組合槳攪拌槽內(nèi)的流場(chǎng)特性進(jìn)行研究，實(shí)驗(yàn)平臺(tái)由攪拌控制實(shí)驗(yàn)裝置和粒子成像測(cè)速裝置組成。實(shí)驗(yàn)開始時(shí)，把攪拌槽放置在儲(chǔ)水方槽中，在攪拌槽中注入甘油和水溶液并攪拌均勻，通過(guò)黏度計(jì)測(cè)得甘油-水溶液黏度為28.3 mPa×s，通過(guò)質(zhì)量除以體積計(jì)算得到質(zhì)量濃度為1 172 kg×m-3，再向儲(chǔ)水方槽中注入清水，使方槽的液位與攪拌槽液位保持齊平，避免折射現(xiàn)象。打開雙脈沖激光器發(fā)射2條激光，在攪拌槽內(nèi)形成2塊明亮的示蹤粒子區(qū)域，相機(jī)連續(xù)拍攝兩幀粒子圖像，從這2張圖像中可以得到在0.001 s內(nèi)示蹤粒子的位移，從而計(jì)算示蹤粒子的速度，本實(shí)驗(yàn)取200組照片計(jì)算得到整個(gè)流場(chǎng)內(nèi)的平均速度分布。

2.2 建模與網(wǎng)格劃分

通過(guò)三維建模軟件分別對(duì)攪拌槽、改進(jìn)型三斜葉-Rushton組合槳及三斜葉-Rushton組合槳進(jìn)行模型構(gòu)建，攪拌槽模型高度為液位高度，完成之后導(dǎo)入ANSYS分析軟件的Workbench模塊當(dāng)中，劃分流體域后再進(jìn)行流體域的布爾運(yùn)算。由于研究的組合槳附近的結(jié)構(gòu)較為復(fù)雜，對(duì)于結(jié)構(gòu)化網(wǎng)格的劃分相對(duì)困難，所以研究的幾何模型采用非結(jié)構(gòu)化四面體網(wǎng)格。為了保證計(jì)算域的精確性，對(duì)上下層槳葉和攪拌軸附近區(qū)域進(jìn)行網(wǎng)格加密處理，劃分網(wǎng)格結(jié)構(gòu)模型如圖2所示。

圖2 計(jì)算域網(wǎng)格示意圖

2.3 網(wǎng)格無(wú)關(guān)性驗(yàn)證

從Hosseini等[23]的研究中可以發(fā)現(xiàn)，當(dāng)網(wǎng)格數(shù)量增加，某一方向上速度矢量的變化量低于3% 時(shí)，則認(rèn)為網(wǎng)格數(shù)量對(duì)數(shù)值模擬產(chǎn)生的結(jié)果無(wú)影響。本研究根據(jù)一定的網(wǎng)格數(shù)量，計(jì)算功率準(zhǔn)數(shù)的相對(duì)偏差來(lái)判定模擬結(jié)果與網(wǎng)格數(shù)量的無(wú)關(guān)性，表2為劃分網(wǎng)格數(shù)與功率準(zhǔn)數(shù)的關(guān)系，觀察到網(wǎng)格數(shù)量為9.8′105時(shí)，功率準(zhǔn)數(shù)的相對(duì)偏差為1.774%，小于3%，此時(shí)認(rèn)為9.8′105的網(wǎng)格數(shù)量為進(jìn)行數(shù)值模擬的網(wǎng)格數(shù)。

表2 劃分網(wǎng)格數(shù)與功率準(zhǔn)數(shù)的關(guān)系

2.4 模擬方法

本研究采用大渦模擬(large eddy simulation, LES)方法進(jìn)行模擬，流場(chǎng)中渦的大小尺度通過(guò)濾波函數(shù)來(lái)劃分，大尺度的渦通過(guò)Navier-Stokes方程計(jì)算，小尺度的渦通過(guò)建立亞格子尺度模型來(lái)模擬。濾波函數(shù)將連續(xù)性方程和動(dòng)量方程進(jìn)行空間濾波后，得到大渦模擬流體運(yùn)動(dòng)的控制方程[24]：

通常標(biāo)準(zhǔn)Smagorinsky-Lily模型中的Smagorinsky常數(shù)S會(huì)對(duì)湍動(dòng)能耗散率產(chǎn)生一定的影響，故采用動(dòng)態(tài)的Smagorinsky常數(shù)S[25]，該模型中局部流體應(yīng)變率與其子應(yīng)力呈反比，具體表達(dá)形式如下：

3 數(shù)值模擬的結(jié)果分析

3.1 改進(jìn)型三斜葉-Rushton組合槳槽內(nèi)速度場(chǎng)分析

3.1.1 垂直剖面速度場(chǎng)分析

對(duì)靜區(qū)域的流體參數(shù)設(shè)置與實(shí)驗(yàn)的甘油-水溶液一致，在不同轉(zhuǎn)速下的攪拌槽中心處垂直剖面速度云圖分布如圖3所示，攪拌槳葉在運(yùn)動(dòng)過(guò)程中形成的射流對(duì)槽底的流體運(yùn)動(dòng)影響較小，轉(zhuǎn)速的提升使攪拌槽底的流體速度在不斷增大，槽內(nèi)流體速度峰值也逐步提高。當(dāng)轉(zhuǎn)速=80 r×min-1時(shí)計(jì)算功率為2.199 W，改進(jìn)型三斜葉槳附近形成了較為明顯的射流，但對(duì)下層Rushton槳形成的射流影響較小，且整個(gè)槽內(nèi)的流體速度較低，混合效果不佳。圖3(b)中=100 r×min-1時(shí)計(jì)算功率為4.246 W，槽底和內(nèi)壁附近的流體速度增大，上下層槳葉附近的流體區(qū)域的湍動(dòng)程度增強(qiáng)，攪拌槽內(nèi)的速度分布得到改善。當(dāng)轉(zhuǎn)速提升為120和140 r×min-1時(shí)，槽底的速度分布在改善，但對(duì)攪拌槽內(nèi)的槳葉區(qū)域速度分布的影響較小，而計(jì)算功率分別為6.883和10.513 W，功耗較大，不利于降低能耗。故最佳轉(zhuǎn)速為100 r×min-1。

圖3 垂直剖面速度云圖

3.1.2 水平剖面速度場(chǎng)分析

不同轉(zhuǎn)速下軸向高度/=0.53處的水平剖面速度云圖如圖4所示，圖中顯示轉(zhuǎn)速提升使相同軸向高度下的云圖分布的速度出現(xiàn)了不同趨勢(shì)的增大。流體沿改進(jìn)型三斜葉槳的折邊流動(dòng)，阻力有所減弱，槳的下半部分帶動(dòng)流體的運(yùn)動(dòng)增強(qiáng)湍動(dòng)能。圖4(a)中由于轉(zhuǎn)速較低，在攪拌槽內(nèi)壁處流體受到攪拌槳的作用較小，速度變化不明顯。當(dāng)轉(zhuǎn)速提升到100 r×min-1時(shí)流體受到較高轉(zhuǎn)速的作用，改進(jìn)型三斜葉槳的外排能力也在增強(qiáng)，流體沿圓周的切向方向運(yùn)動(dòng)，此時(shí)攪拌槽內(nèi)壁的低速區(qū)域的范圍也在減小。隨著轉(zhuǎn)速提升，如圖4(c)、(d)所示，流體速度峰值逐步提高，但對(duì)槳葉附近湍動(dòng)能較高區(qū)域的范圍影響較小。不利于節(jié)約能耗。

圖4 水平剖面速度云圖

3.2 組合槳槽內(nèi)速度場(chǎng)對(duì)比分析

3.2.1 速度流線云圖對(duì)比分析

當(dāng)=100 r×min-1時(shí)，2種不同組合槳的速度流線云圖如圖5所示。圖5(a)中三斜葉-Rushton組合槳中三斜葉槳對(duì)Rushton槳的作用相對(duì)較弱，未能在Rushton槳形成射流的上側(cè)形成渦旋。觀察圖5(b)中速度流型云圖的變化，與圖5(a)相比Rushton槳形成射流的方向?yàn)樾毕蛳路?，在射流上?cè)形成渦旋。由于改進(jìn)型三斜葉槳折邊上側(cè)的流體沿折邊運(yùn)動(dòng)，下側(cè)的流體受到折邊刮擋的作用，所以其泵吸作用相對(duì)于三斜葉的泵吸作用更強(qiáng)，軸向?qū)α餮h(huán)性能更好。改進(jìn)型三斜葉槳對(duì)Rushton槳形成的射流具有更強(qiáng)影響，縮小了攪拌槽底部和內(nèi)壁附近的低速區(qū)域，且射流的上側(cè)有渦旋形成，增強(qiáng)了整體循環(huán)性能。當(dāng)=100 r×min-1時(shí)，三斜葉-Rushton組合槳的計(jì)算功率為4.424 W，而改進(jìn)型三斜葉-Rushton組合槳為4.246 W，能耗更低，這是由于三斜葉槳沒有折邊的存在，葉片所受槳葉周邊流體的阻力較大。

圖5 2種組合槳的速度流線云圖

3.2.2 無(wú)因次化速度對(duì)比分析

圖6為=100 r×min-1時(shí)徑向位置/=0.3處的無(wú)因次化速度分布。圖6(a)中2種組合槳的上槳葉區(qū)出現(xiàn)了相反方向的速度分布，主要是因?yàn)楦倪M(jìn)型三斜葉槳的泵吸作用較強(qiáng)，增強(qiáng)了槳葉附近流體區(qū)域的湍動(dòng)能，在上槳葉區(qū)出現(xiàn)了渦旋，而在三斜葉槳-Rushton組合槳槽內(nèi)的上槳葉區(qū)則沒有出現(xiàn)渦旋。2種不同的組合槳形成的射流及兩槳葉間的徑向無(wú)因次化速度的趨勢(shì)大致相同，但在Rushton槳附近的/tip值改進(jìn)型三斜葉-Rushton組合槳大于三斜葉-Rushton組合槳。觀察到圖6(b)中2種組合槳的軸向無(wú)因次化速度分布的趨勢(shì)相一致，與圖6(a)徑向無(wú)因次化速度分布不同的是三斜葉-Rushton組合槳兩槳葉間的/tip值大于改進(jìn)型三斜葉-Rushton組合槳，但這對(duì)于整體的速度分布的影響較小。

圖6 無(wú)因次化速度分布圖

4 實(shí)驗(yàn)的結(jié)果及討論

4.1 轉(zhuǎn)速對(duì)流場(chǎng)的影響

實(shí)驗(yàn)工況與模擬相同的4種不同轉(zhuǎn)速下的速度云圖如圖7所示。當(dāng)=80 r×min-1時(shí)，實(shí)驗(yàn)功率為2.502 W，流場(chǎng)中的速度峰值偏小，在0.18

圖7 4種不同轉(zhuǎn)速的速度云圖

4.2 槳葉間距對(duì)流場(chǎng)特性的影響

在甘油-水溶液中=100 r×min-1和底槳離底距離1=0.3的條件下，3種不同槳葉間距攪拌槽內(nèi)流場(chǎng)的速度云圖分布如圖8所示。圖8(a)中2=0.22時(shí)，改進(jìn)型三斜葉槳的前端射流導(dǎo)致兩槳間流場(chǎng)區(qū)域在槳葉附近區(qū)域出現(xiàn)了較大的速度，但由于槳葉間距較小使槳葉對(duì)槽內(nèi)上部分區(qū)域影響弱，所以在/=0.19～0.50，/=0.65～1.0出現(xiàn)了較大范圍的低速區(qū)域，不利于攪拌槽內(nèi)上槳葉區(qū)流體的混合。如圖8(b)所示槳葉間距2=0.25時(shí)，兩槳葉區(qū)流體速度的分布較為均勻，上槳葉區(qū)較低速度區(qū)域減小，相比槳葉間距2=0.22時(shí)混合效率提升。圖8(c)顯示當(dāng)2=0.29時(shí)上下槳葉對(duì)兩槳葉間的流體作用有所減弱，致使兩槳葉間出現(xiàn)低速區(qū)域，槽內(nèi)流體的混合效果降低。故最佳槳葉間距2為0.25。

圖8 3種不同槳葉間距的速度云圖

4.3 離底距離對(duì)流場(chǎng)特性的影響

在甘油-水溶液中=100 r×min-1和兩槳間距離2=0.25的條件下，3種不同離底距離對(duì)流場(chǎng)產(chǎn)生影響的速度云圖分布如圖9所示。當(dāng)離底距離1=0.23時(shí)，上槳葉區(qū)的流體的攪拌效果較差，出現(xiàn)大面積的低速區(qū)域。在下槳葉區(qū)，Rushton槳浸沒在槽底，流體的流動(dòng)主要在槽底區(qū)域。當(dāng)1=0.26時(shí)上槳葉區(qū)的低速區(qū)域相對(duì)離底距離1=0.23有所縮小，由0.56

圖9 3種不同離底距的速度云圖

5 結(jié) 論

本研究分別通過(guò)CFD模擬和PIV實(shí)驗(yàn)對(duì)裝配改進(jìn)型三斜葉-Rushton組合槳攪拌槽內(nèi)的流場(chǎng)進(jìn)行了研究，得到的主要結(jié)論如下：

(1) 轉(zhuǎn)速的改變影響槽內(nèi)流場(chǎng)的速度分布，當(dāng)轉(zhuǎn)速為80 r×min-1時(shí)，在攪拌槽內(nèi)流場(chǎng)的速度較小，在兩槳葉間呈現(xiàn)低速不均勻的現(xiàn)象。當(dāng)轉(zhuǎn)速為100 r×min-1時(shí)，兩槳葉之間的速度分布相對(duì)均勻，提高了攪拌槽底部和內(nèi)壁附近的流體速度，槽內(nèi)整體混合效果更好。轉(zhuǎn)速的提高使上槳葉區(qū)的較低速度區(qū)域在縮小，但槳葉區(qū)域速度分布沒有進(jìn)一步改善，不利于降低功耗。故=100 r×min-1為最佳轉(zhuǎn)速，實(shí)驗(yàn)結(jié)果和模擬結(jié)果吻合效果較好。

(2) 通過(guò)CFD模擬發(fā)現(xiàn)改進(jìn)型三斜葉槳的折邊具有刮擋作用，使其泵吸作用較三斜葉槳更強(qiáng)，外排液體呈現(xiàn)較強(qiáng)的軸流特性使軸向?qū)α餮h(huán)性能更強(qiáng)，更有利于兩槳葉之間液體的混合。改進(jìn)型三斜葉-Rushton組合槳相比三斜葉-Rushton組合槳能有效縮小攪拌槽底部和內(nèi)壁附近的低速區(qū)域，增強(qiáng)整體循環(huán)性能，且功耗更低。

(3) 通過(guò)PIV實(shí)驗(yàn)對(duì)比研究的工況，當(dāng)2=0.25時(shí)為最佳槳葉間距，上槳葉區(qū)較低速度區(qū)域有所減小，兩槳葉間的流體運(yùn)動(dòng)逐漸趨于平穩(wěn)，槽內(nèi)速度分布更合理，混合效果更好；當(dāng)1=0.3時(shí)為最佳離底距離，能縮小上槳葉區(qū)的較低速區(qū)域，整個(gè)槽內(nèi)的混合能力增強(qiáng)。

b0?折葉片的寬度，mm?解析應(yīng)變速率張量，無(wú)量綱 b1?折邊的寬度，mm?解析應(yīng)變速率張量的范數(shù) C1?底槳離底距離，mmt?時(shí)間，s C2?兩槳間距離，mmt1?攪拌槽壁厚，mm CS?亞格子常數(shù)，無(wú)量綱t2?槳葉葉片厚度，mm D?槳葉直徑，mm?速度矢量，m×s-1 d?攪拌槽內(nèi)徑，mmu/Vtip?徑向無(wú)因次化速度，無(wú)量綱 d1?渦與壁面的距離，mmV?計(jì)算網(wǎng)格單元體積，mm3 ?質(zhì)量力，Nv/Vtip?軸向無(wú)因次化速度，無(wú)量綱 H?攪拌槽高度，mmx?方向，無(wú)量綱 h?攪拌槽內(nèi)液位高度，mz/h?軸向位置，無(wú)量綱 L1?折邊的長(zhǎng)度，mmθ?折葉彎折角度，(°) L2?折葉片的長(zhǎng)度，mmк?卡門常數(shù)，無(wú)量綱 LS?亞格子尺度混合長(zhǎng)度，mm?亞格子正應(yīng)力，N NP?功率準(zhǔn)數(shù)，無(wú)量綱?亞格子應(yīng)力，N n?槳葉轉(zhuǎn)速，r×min-1?亞格子渦粘系數(shù)，無(wú)量綱 ?壓力脈動(dòng)值，N下標(biāo) r/d?徑向位置，無(wú)量綱i , j?坐標(biāo)軸方向的分量

[1] 陳甘棠. 聚合反應(yīng)工程基礎(chǔ) [M]. 北京: 中國(guó)石化出版社, 1991.

CHEN G T. Engineering foundation of polymerization reaction [M]. Beijing: China Petrochemical Press, 1991.

[2] 周丹, 張勇, 蔣軍成, 等. 苯乙烯聚合反應(yīng)熱失控及其緊急抑制的數(shù)值模擬 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2016, 30(3): 618-625.

ZHOU D, ZHANG Y, JIANG J C,. Numerical simulation on runaway reaction and inhibition of styrene thermal polymerization [J]. Journal of Chemical Engineering of Chinese Universities, 2016, 30(3): 618-625.

[3] XU C Z, WANG J J, GU X P,. CFD modeling of styrene polymerization in a CSTR [J]. Chemical Engineering Research and Design, 2017, 125: 46-56.

[4] 丁緒淮, 周理. 液體攪拌 [M]. 北京: 化學(xué)工業(yè)出版社, 1983.

DING X H, ZHOU L. Liquid mixing [M]. Beijing: Chemical Industry Press, 1983.

[5] 楊鋒苓, 周慎杰, 王貴超, 等. 非標(biāo)準(zhǔn)擋板攪拌槽內(nèi)湍流流場(chǎng)的數(shù)值模擬 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2012, 26(6): 952-958.

YANG F L, ZHOU S J, WANG G C,. Study on the hydrodynamics of a stirred tank equipped with non-standard baffles [J]. Journal of Chemical Engineering of Chinese Universities, 2012, 26(6): 952-958.

[6] 陳濤, 吳大轉(zhuǎn), 杜紅霞, 等. 低密度固體顆粒與液體攪拌混合過(guò)程的數(shù)值模擬研究 [J]. 工程熱物理學(xué)報(bào), 2010,31(2): 271-274.

CHEN T, WU D Z, DU H X,. Numerical simulation on low-density solid-liquid mixing process [J]. Journal of Engineering Thermophysics, 2010, 31(2): 271-274.

[7] 李志鵬, 高正明. 渦輪槳攪拌槽內(nèi)流動(dòng)特性的大渦模擬 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2007, 21(4): 592-597.

LI Z P, GAO Z M. Large eddy simulation of flow field in a rushton impeller stirred tank [J]. Journal of Chemical Engineering of Chinese Universities, 2007, 21(4): 592-597.

[8] TAMBURINI A, GAGLIANO G, MICALE G,. Direct numerical simulations of creeping to early turbulent flow in unbaffled and baffled stirred tanks [J]. Chemical Engineering Science, 2018, 192: 161-175.

[9] 高娜, 包雨云, 高正明. 多層槳攪拌槽內(nèi)氣-液兩相局部氣含率研究 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2011, 25(1): 11-17.

GAO N, BAO Y Y, GAO Z M. Local void fraction in a sparged reactor with a multi-impeller agitator [J]. Journal of Chemical Engineering of Chinese Universities, 2011, 25(1): 11-17.

[10] LI G, GAO Z, LI Z,. Particle-resolved PIV experiments of solid-liquid mixing in a turbulent stirred tank [J]. AIChE Journal, 2018, 64(1): 389-402.

[11] 馮俊杰, 孫冰, 姜杰, 等. 射流場(chǎng)中單個(gè)氣泡破裂的CFD模擬與PIV實(shí)驗(yàn)研究 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2018, 32(2): 286-294.

FENG J J, SUN B, JIANG J,. CFD simulation and PIV investigation on single bubble breakup in turbulent jet flow [J]. Journal of Chemical Engineering of Chinese Universities, 2018, 32(2): 286-294.

[12] ZHOU Y J, WANG L Z, HE H,. Mixing process in a tank stirred with improved double intermig impellers [J]. Journal of Chemical Engineering of Japan, 2019, 52(9): 719-729.

[13] 鮑蘇洋, 周勇軍, 王璐璐, 等. 渦輪槳攪拌槽內(nèi)流場(chǎng)特性的V3V實(shí)驗(yàn) [J]. 化工學(xué)報(bào), 2016, 67(11): 4580-4586.

BAO S Y, ZHOU Y J, WANG L L,. V3V study on flow field characteristics in a stirred vessel with Rushton turbine impeller [J]. CIESC Journal, 2016, 67(11): 4580-4586.

[14] 施乃進(jìn), 周勇軍, 鮑蘇洋, 等. 渦輪槳攪拌槽內(nèi)湍流特性的V3V實(shí)驗(yàn)及大渦模擬 [J]. 化工學(xué)報(bào), 2017, 68(11): 4069-4078.

SHI N J, ZHOU Y J, BAO S Y,. V3V study and large eddy simulation of turbulence characteristics in a stirred vessel with Rushton turbine impeller [J]. CIESC Journal, 2017, 68(11): 4069-4078.

[15] 謝明輝, 魯洪中, 唐文俊, 等. 多層攪拌式生物反應(yīng)器內(nèi)冷模和熱模的傳質(zhì)特性 [J]. 化工學(xué)報(bào), 2016, 67(S1): 224-231.

XIE M H, LU H Z, TANG W J,. Mass transfer performance for multiple-impeller bioreactor in cold model and fermentation experiments [J]. CIESC Journal, 2016, 67(S1): 224-231.

[16] 朱姝, 包雨云, 陳雷, 等. 用電導(dǎo)探針測(cè)定氣-液多層槳攪拌槽內(nèi)氣泡尺寸分布 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2011, 25(6): 977-984.

ZHU S, BAO Y Y, CHEN L,. Bubble size distributions measurement in a gas-liquid multi-impeller stirred tank by using dual-conductivity probe [J]. Journal of Chemical Engineering of Chinese Universities, 2011, 25(6): 977-984.

[17] 楊娟, 張慶華, 楊超, 等. 不同組合槳攪拌槽內(nèi)非牛頓流體的微觀混合特性 [J]. 過(guò)程工程學(xué)報(bào), 2019, 19(5): 865-871.

YANG J, ZHANG Q H, YANG C,. Micro-mixing characteristics of non-Newtonian fluid in a stirred tank agitated with different impellers [J]. The Chinese Journal of Process Engineering, 2019, 19(5): 865-871.

[18] 周勇軍, 袁名岳, 孫存旭. 改進(jìn)型框式組合槳攪拌槽內(nèi)流場(chǎng)特性 [J]. 化工進(jìn)展, 2019, 38(12): 5306-5313.

ZHOU Y J, YUAN M Y, SUN C X. Investigating on flow field in stirred tank equipped with improved frame type combined impellers [J]. Chemical Industry and Engineering Progress, 2019, 38(12): 5306-5313.

[19] 吳高杰, 陳劍佩, 黃建科, 等. 攪拌對(duì)多黏類芽孢桿菌發(fā)酵氧傳遞過(guò)程的影響及CFD模擬 [J]. 華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 35(3): 339-345.

WU G J, CHEN J P, HUANG J K,. Effect of agitation on oxygen transfer in paenibacillus polymyxa HY96-2 culture broth and CFD simulation [J]. Journal of East China University of Science and Technology ( Natural Science Edition), 2009, 35(3): 339-345.

[20] 宮磊, 周麗, 崔文璟, 等. 攪拌槳組合數(shù)值模擬優(yōu)化及在谷氨酰胺轉(zhuǎn)胺酶發(fā)酵中的應(yīng)用 [J]. 化工進(jìn)展, 2014, 33(10): 2570-2575.

GONG L, ZHOU L, CUI W J,. Numerical simulation of different impellers and its effect on transglutaminase fermentation process [J]. Chemical Industry and Engineering Progress, 2014, 33(10): 2570-2575.

[21] 劉寶慶, 張義堃, 劉景亮, 等. 新型同心雙軸攪拌器功率與混合特性的數(shù)值模擬 [J]. 化工學(xué)報(bào), 2013, 64(4): 1135-1144.

LIU B Q, ZHANG Y K, LIU J L,. Numerical simulation of power consumption and mixing characteristic in stirred vessel with novel coaxial mixer [J]. CIESC Journal, 2013, 64(4): 1135-1144.

[22] ZHOU Y J, LIN W Z, YUAN M Y,. Investigation on the flow field and mixing efficiency of a stirred tank equipped with improved intermig impellers [J]. International Journal of Chemical Reactor Engineering, 2019, 17(11): 1-16.

[23] HOSSEINI S, PATEL D, EIN-MOZAFFARI F,. Study of solid?liquid mixing in agitated tanks through computational fluid dynamics modeling [J]. Industrial & Engineering Chemistry Research, 2010, 49(9): 4426-4435.

[24] 張兆順, 崔桂香, 許春曉, 等. 湍流理論與模擬 [M]. 北京: 清華大學(xué)出版社, 2017.

ZHANG Z X, CUI G X, XU C X,. Turbulence theory and simulation [M]. Beijing: Tsinghua University Press, 2017.

[25] MALIK S, Lévêquea E, BOUAIFIBC M,. Shear improved Smagorinsky model for large eddy simulation of flow in a stirred tank with a Rushton disk turbine [J]. Chemical Engineering Research & Design, 2016, 108(S1): 69-80.

Flow field characteristics in a stirred tank with improved three-pitched blade-Rushton combined impellers

CHEN Shuai1, ZHOU Yong-jun1, LIN Wei-zhen1, SHI Yi-hang1, TAO Lan-lan2

(1. College of Mechanical and Power Engineering, Nanjing Tech University, Nanjing 211816, China;2. National Quality Supervision & Inspection Center of Chemical Equipment, Special Equipment Safety Supervision Inspection Institute of Jiangsu Province, Suzhou 215600, China)

In order to efficiently apply an improved three-pitched blade-Rushton combined impeller in practical engineering, the flow field in a stirred tank equipped with the improved three-pitched blade-Rushton combined impeller was simulated by computational fluid dynamics (CFD), and the flow field of the improved three-pitched blade-Rushton combined impeller was compared with that of the three-pitched blade-Rushton combined impeller. In addition, the particle image velocimetry (PIV) technology was used to study the effects of rotational speed, impeller pitch2and installation height of impellers1on the flow field of the improved three-pitched blade-Rushton combined impeller in stirred tanks. The result shows that when the rotational speed=100 r×min-1, the mixing effect of the bottom and upper blade areas in the stirred tank can be effectively improved, and the velocity distribution of the flow field in the stirred tank is relatively uniform. The experimental results are in good agreement with the simulation results. Numerical simulation results demonstrate that the improved three-pitched blade-Rushton combined impeller has stronger pumping effects than the three-pitched blade-Rushton combined impeller, and it can effectively reduce the area of low-velocity near the bottom and inner wall in the stirred tank, which results in better mixing effects. The experimental results show that the velocity distribution is more reasonable when2=0.25, and the lower velocity area of the upper blade area is further reduced when the installation height of impellers1=0.3.

stirred tank; improved three-pitched blade-Rushton combined impeller; flow field; computational fluid dynamics(CFD); particle image velocimetry(PIV); mixing

TQ 027.2

A

10.3969/j.issn.1003-9015.2021.03.009

1003-9015(2021)03-0455-09

2020-07-28；

2020-10-20。

國(guó)家自然科學(xué)基金(51775262)。

陳帥(1996-)，男，湖北赤壁人，南京工業(yè)大學(xué)碩士生。

周勇軍，E-mail：zhouyj@njtech.edu.cn