朱俊，周政霖，劉作華,，鄭雄攀，劉仁龍，陶長(zhǎng)元，王運(yùn)東

（1重慶大學(xué)化學(xué)化工學(xué)院，重慶 400044；2清華大學(xué)化學(xué)工程系，北京 100084）

剛?cè)峤M合攪拌槳強(qiáng)化流體混合的流固耦合行為

朱俊1，周政霖1，劉作華1,2，鄭雄攀1，劉仁龍1，陶長(zhǎng)元1，王運(yùn)東2

（1重慶大學(xué)化學(xué)化工學(xué)院，重慶 400044；2清華大學(xué)化學(xué)工程系，北京 100084）

傳統(tǒng)剛性攪拌槳通過(guò)對(duì)流體的剪切作用實(shí)現(xiàn)能量的傳遞，而剛?cè)峤M合攪拌槳可通過(guò)其多體運(yùn)動(dòng)行為強(qiáng)化能量傳遞。基于攪拌槳槳葉與流體之間的耦合運(yùn)動(dòng)作用，結(jié)合ANSYS Workbench仿真平臺(tái)，采用雙向流固耦合方法，模擬計(jì)算了剛性攪拌槳與剛?cè)峤M合攪拌槳槳葉的等效應(yīng)力和總變形量，研究了流場(chǎng)的宏觀結(jié)構(gòu)；并通過(guò)測(cè)定混合時(shí)間和計(jì)算攪拌槳功耗對(duì)比分析了兩種不同攪拌體系的混合行為。結(jié)果表明：剛?cè)峤M合攪拌槳使體系的混合時(shí)間縮短了近32%，攪拌槳功耗下降了7%，其槳葉尖端的變形量是剛性攪拌槳的105倍，其應(yīng)力比剛性攪拌槳增加了83%；與剛性攪拌槳相比，剛?cè)峤M合攪拌槳在流固耦合作用下對(duì)流體的作用力更大，能夠更好地傳遞能量，增強(qiáng)流體運(yùn)動(dòng)，強(qiáng)化流體混合。

混合；流體；傳遞；剛?cè)峤M合攪拌槳；ANSYS Workbench；流固耦合

引 言

攪拌槳作為攪拌混合設(shè)備中的重要部件之一，提供了攪拌過(guò)程中的所需能量和適宜流動(dòng)狀態(tài)，對(duì)它進(jìn)行合理設(shè)計(jì)是流體實(shí)現(xiàn)高效、節(jié)能混合的重要途徑[1-2]。研究表明，攪拌過(guò)程中能量的耗散主要發(fā)生在槳葉外緣和槳葉后的尾渦處，約占攪拌槳輸入能量的70%。以剛性材質(zhì)為主的傳統(tǒng)攪拌槳通過(guò)對(duì)流體的剪切與聚并作用傳遞能量，導(dǎo)致混合能耗高、效率低；而剛?cè)狁詈闲问降臄嚢铇ㄟ^(guò)其柔性端與周圍流體相互作用，從“渦街”中汲取能量，提高流體混合效率[3-5]。事實(shí)上，攪拌槳結(jié)構(gòu)與流體混合效率有直接關(guān)系，柔性槳葉與流體之間存在力的耦合作用，柔性槳葉在流體載荷的作用下產(chǎn)生變形或運(yùn)動(dòng)，柔性槳葉的變形或運(yùn)動(dòng)又反過(guò)來(lái)影響流場(chǎng)，引起流場(chǎng)載荷分布和大小的變化，從而影響流體混合行為[6-8]。因此，流固耦合分析方法可用于指導(dǎo)攪拌槽內(nèi)剛?cè)峤M合攪拌槳設(shè)計(jì)與研究流體混合行為。

目前，流固耦合分析方法廣泛用于土木、航天航空、船舶等工程領(lǐng)域[9-10]。Munch等[11]、Gao等[12]運(yùn)用雙向流固耦合方法模擬分析了人體主動(dòng)脈內(nèi)血液的流動(dòng)情況。Young[13-14]采用面元法和有限元法建立了復(fù)合材料螺旋槳的流固耦合方程，研究了螺旋槳的振動(dòng)頻率和振型問(wèn)題，開(kāi)發(fā)了一種敞水性能優(yōu)良的螺旋槳。Bucchignani等[15]通過(guò)耦合流體域和結(jié)構(gòu)域兩個(gè)不同離散化模型的方法，研究了攪拌槽內(nèi)不可壓縮流體的混合行為，發(fā)現(xiàn)攪拌槳槳葉與流體之間存在力的耦合作用。結(jié)合流體動(dòng)力學(xué)及固體動(dòng)力學(xué)的計(jì)算方法，Karry等[16]研究了錨式攪拌槳的流固耦合行為，發(fā)現(xiàn)水壓力影響槳葉變形。周思柱等[17]對(duì)混砂車攪拌體系進(jìn)行了流固耦合分析，獲得了離心力、流體壓力等載荷對(duì)葉輪的影響。劉作華等[18-19]研發(fā)了組合式的剛?cè)峤M合攪拌槳，并進(jìn)行了能效分析及混沌特性研究，但剛?cè)峤M合攪拌槳流固耦合運(yùn)動(dòng)行為的研究還未見(jiàn)報(bào)道。目前，人們大多采用單向流固耦合方法研究攪拌體系中的流固耦合問(wèn)題，該方法不利于分析變形大的槳葉對(duì)流場(chǎng)結(jié)構(gòu)的影響。

本文結(jié)合ANSYS Workbench平臺(tái)，采用雙向流固耦合方法對(duì)剛?cè)峤M合攪拌槳強(qiáng)化流體混合進(jìn)行了研究，其中考察了槳葉總變形量、等效應(yīng)力及流場(chǎng)宏觀結(jié)構(gòu)，并用剛-柔-流耦合運(yùn)動(dòng)行為闡釋了強(qiáng)化流體混合的原理。同時(shí)，通過(guò)測(cè)定混合時(shí)間和計(jì)算攪拌槳功耗對(duì)比分析了兩種體系的混合行為，為剛?cè)峤M合攪拌槳強(qiáng)化流體混合提供實(shí)驗(yàn)理論依據(jù)。

1 實(shí)驗(yàn)裝置與數(shù)值計(jì)算

1.1 實(shí)驗(yàn)裝置

本實(shí)驗(yàn)裝置與劉作華等[20]在剛?cè)峤M合攪拌槳與剛性槳調(diào)控流場(chǎng)結(jié)構(gòu)的對(duì)比研究中所采用的設(shè)備一致。攪拌槽為無(wú)擋板平底圓柱形有機(jī)玻璃攪拌槽（圖1），槽內(nèi)徑T為200 mm，液面高度H為140 mm。單層槳配置，攪拌槳離底距離C為70 mm，類型為六斜葉圓盤剛?cè)峤M合攪拌槳[稱組合槳（RF-RDT）]，直徑D為150 mm，傾角θ為45°，短槳葉葉寬m1為15 mm、葉長(zhǎng)l1為20 mm、材質(zhì)為不銹鋼，長(zhǎng)槳葉葉寬m2為8 mm、葉長(zhǎng)l2為8 mm、材質(zhì)為硅膠。對(duì)應(yīng)攪拌槳為六斜葉圓盤剛性攪拌槳[簡(jiǎn)稱剛性槳（RDT）]，其構(gòu)型及大小與組合槳相同，差異在于其短槳葉和長(zhǎng)槳葉的材質(zhì)都為不銹鋼。攪拌槳結(jié)構(gòu)如圖2所示。

圖1 攪拌槽結(jié)構(gòu)Fig.1 Schematic illustration for stirred tank

實(shí)驗(yàn)工作介質(zhì)為水，25℃時(shí)，經(jīng)測(cè)定，水的密度為997.05 kg·m?3,黏度為0.8937 mPa·s。整個(gè)實(shí)驗(yàn)在相同轉(zhuǎn)速N為120 r·min?1下進(jìn)行。

1.2 數(shù)值計(jì)算

攪拌槽內(nèi)流體流動(dòng)的基本守恒定律包括質(zhì)量守恒定律、動(dòng)量守恒定律。這些守恒方程的數(shù)學(xué)描述便是控制方程[16]。對(duì)于一般的牛頓流體，守恒定律通過(guò)如下控制方程描述

圖2 攪拌槳結(jié)構(gòu)Fig.2 Schematic illustration for impellers

固體攪拌槳槳葉的動(dòng)力學(xué)守恒方程可以由牛頓第二定律導(dǎo)出

在流固耦合交界面處，流固耦合方程應(yīng)滿足流體與攪拌槳槳葉應(yīng)力（τ）、位移（d）等變量的守恒

數(shù)值計(jì)算首先利用Gambit軟件對(duì)攪拌槽和攪拌槳建立三維參數(shù)化模型，并劃分網(wǎng)格和設(shè)置邊界條件（圖3）。然后在ANSYS Workbench平臺(tái)中建立流固耦合計(jì)算體系，將攪拌槳模型導(dǎo)入Finite Element Modeler模塊進(jìn)行有限元化處理。接著，將經(jīng)過(guò)有限元化后的攪拌槳模型導(dǎo)入Transient Structural模塊設(shè)置固體面，剛性材質(zhì)不銹鋼密度為7850 kg·m?3、彈性模量為2×105MPa、泊松比為0.3，柔性材質(zhì)硅膠密度為1200 kg·m?3、彈性模量為2.14 MPa、泊松比為0.48；同時(shí)，將攪拌槽模型導(dǎo)入Fluent模塊設(shè)置流體面，結(jié)合多重參考系MRF模型將其劃分為靜區(qū)域和動(dòng)區(qū)域，并采用彈簧光順和局部網(wǎng)格重構(gòu)來(lái)防止網(wǎng)格在計(jì)算時(shí)因拉伸和彎曲導(dǎo)致的結(jié)果不收斂。最后將兩個(gè)模型一并導(dǎo)入System Coupling模塊進(jìn)行求解計(jì)算。此外，雙向流固耦合對(duì)流體域和固體域同時(shí)求解，時(shí)間步長(zhǎng)須保持一致，設(shè)置為0.001 s。整個(gè)模擬計(jì)算過(guò)程在內(nèi)存為32 G的戴爾T5600工作站上進(jìn)行。

圖3 網(wǎng)格劃分示意圖Fig.3 Schematic illustration for meshing

本文計(jì)算重點(diǎn)關(guān)注的是攪拌槽內(nèi)流體的混合行為，而混合行為與流體的速度場(chǎng)分布密切相關(guān)，故選取槽內(nèi)X=75 mm、Y=0 mm、Z從?70 mm到70 mm的直線line A上的流體速度分布來(lái)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證。以剛性槳攪拌槽瞬態(tài)數(shù)值模擬為例，在計(jì)算10 s后發(fā)現(xiàn)，當(dāng)網(wǎng)格數(shù)從約26萬(wàn)個(gè)變至約70萬(wàn)個(gè)時(shí)，3條速度曲線隨網(wǎng)格數(shù)增加基本重合（圖4），此時(shí)的網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果影響很小，可認(rèn)為已達(dá)到網(wǎng)格無(wú)關(guān)，這與劉作華等[20]在剛?cè)峤M合攪拌槳與剛性槳調(diào)控流場(chǎng)結(jié)構(gòu)的對(duì)比研究中所采用的方法基本一致，具有一定的有效性。因此，在本文計(jì)算中，剛性槳攪拌槽網(wǎng)格數(shù)可為399929個(gè)；同理，組合槳攪拌槽網(wǎng)格數(shù)可為387810個(gè)。

圖4 流體合速度分布Fig.4 Velocity magnitude profile for fluid

2 結(jié)果與討論

2.1 模擬分析

2.1.1 攪拌槳槳葉剛度對(duì)比 圖5表明，攪拌槳槳葉上最容易發(fā)生變形的位置位于槳葉尖端，從槳葉尖端到槳葉根部的變形量依次減小，位于圓盤附近的槳葉基本不發(fā)生變形。計(jì)算結(jié)果顯示，組合槳槳葉尖端的變形量為5.89 mm，而剛性槳?jiǎng)t為1.02×10?5mm，兩者相比，組合槳槳葉尖端的變形量是剛性槳的105倍。槳葉變形是由離心力和流體壓力共同引起的，介于兩種攪拌槳的轉(zhuǎn)速相同，槳葉受到的離心力亦相同。與剛性槳相比，組合槳槳葉受到流體壓力更大。介于槳葉與流體之間存在力的耦合作用，使具有高彈性的組合槳對(duì)流體的作用力更強(qiáng)。

圖5 槳葉總變形云圖Fig.5 Total deformation contour for impellers

2.1.2 攪拌槳槳葉強(qiáng)度對(duì)比 圖6反映槳葉與流體作用的過(guò)程中，因槳葉受力不均勻，槳葉根部與圓盤接觸地方的應(yīng)力非常集中，這些地方被稱為局部應(yīng)力集中區(qū)域，此處槳葉的安全系數(shù)比較低，攪拌槳可能發(fā)生機(jī)械失效。由圖6所示，組合槳槳葉應(yīng)力集中區(qū)域的最大應(yīng)力為3.19 kPa，而剛性槳?jiǎng)t為7.58 kPa。與剛性槳相比，組合槳槳葉應(yīng)力集中區(qū)域的最大應(yīng)力下降了58%，即采用組合槳可以緩解槳葉上的疲勞程度，增加攪拌槳的使用年限。槳葉尖端處出現(xiàn)的應(yīng)力主要是由槳葉工作時(shí)復(fù)雜的流體動(dòng)力環(huán)境加上葉輪離心力的共同作用引起的，其中流體動(dòng)力作用占主體地位。計(jì)算結(jié)果表明，組合槳槳葉尖端處的應(yīng)力為1.85 Pa，而剛性槳?jiǎng)t為1.01 Pa，兩者相比，組合槳槳葉尖端處的應(yīng)力上升了83%。由于攪拌槳槳葉與流體之間存在力的耦合效應(yīng)，采用組合槳可提高其對(duì)流體的作用力。

圖6 槳葉等效應(yīng)力云圖Fig.6 Equivalent stress contour for impellers

圖7 體系合速度分布云圖Fig.7 Integrated velocity contour for systems(7 s)

2.1.3 體系流場(chǎng)宏觀結(jié)構(gòu)對(duì)比 流場(chǎng)結(jié)構(gòu)的運(yùn)移和演化可有效地反映流體內(nèi)部傳遞行為，它受攪拌轉(zhuǎn)速等多種因素的共同作用，具有復(fù)雜的非線性現(xiàn)象和時(shí)空混沌行為[20]。為研究方便，僅考慮轉(zhuǎn)速恒定時(shí)的流場(chǎng)結(jié)構(gòu)。圖7和圖8為組合槳和剛性槳在不同時(shí)間下的合速度分布。當(dāng)計(jì)算時(shí)間為7或10 s時(shí)，槽內(nèi)底部區(qū)域流體在剛性槳作用下的流速為0～0.05 m·s?1，底部“死區(qū)”范圍較廣，最底部區(qū)域流體的流速接近0，不利于流體充分混合；而在組合槳作用下，槽內(nèi)底部區(qū)域流體的流速為0.1～0.3 m·s?1，最底部區(qū)域流體的流速為0.05～0.1 m·s?1，這是因?yàn)槠錁~尖端的柔性葉片在局部三維空間上存在“擾動(dòng)”和“波動(dòng)”行為，它能夠促進(jìn)能量從槳葉尖端傳遞給流體，使流體流速加強(qiáng)。與剛性槳相比，組合槳通過(guò)其柔性尖端與周圍流體的相互作用提高了槽內(nèi)底部區(qū)域流體的流速，減小了攪拌過(guò)程中流體存在的“死區(qū)”現(xiàn)象（即圖7或圖8的紅框區(qū)域），可提高流體的混合效率。

分析圖9可知，剛性槳槳葉上下區(qū)域流體的最大軸向速度為0.16 m·s?1；而組合槳為0.24 m·s?1，且流體的軸向速度分布更加均勻（即圖9的紅框區(qū)域）。兩者相比，組合槳強(qiáng)化了流體的軸向流，使更多的流體從槽內(nèi)底部運(yùn)動(dòng)到液面頂部，能夠促進(jìn)全槽流體的充分混合。攪拌槽內(nèi)流體在運(yùn)動(dòng)過(guò)程中容易形成使流體整體運(yùn)動(dòng)的“柱狀回流”，不利于能量耗散。從圖10可看出，剛性槳體系中，槽壁流體切向速度為0.6 m·s?1，圓盤流體切向速度為0.3 m·s?1；而在組合槳體系中，槽壁流體切向速度為0.35～0.45 m·s?1，圓盤流體切向速度為0.2～0.3 m·s?1。兩者相比，組合槳體系的槽壁和圓盤附近流體的切向速度分布更加均勻（即圖10的藍(lán)框和紅框區(qū)域），流體運(yùn)動(dòng)中的“柱狀回流”現(xiàn)象減少，流體混合加強(qiáng)。結(jié)合實(shí)驗(yàn)現(xiàn)象，分析其原因是組合槳槳葉柔性尖端的“波動(dòng)”和“擾動(dòng)”現(xiàn)象能夠加強(qiáng)能量在槳葉尖端和流體內(nèi)部的傳遞，使得全槽流體可獲得更加有效的能量，實(shí)現(xiàn)能量的充分利用。

圖8 體系合速度分布云圖Fig.8 Integrated velocity contour for systems(10 s)

圖9 體系軸向速度分布云圖Fig.9 Axial velocity contour for systems

2.2 混合實(shí)驗(yàn)

攪拌槳功耗及混合時(shí)間是描述混合性能的兩個(gè)重要指標(biāo)，是衡量攪拌反應(yīng)器內(nèi)部混合行為的重要因素，也是評(píng)價(jià)攪拌槳設(shè)計(jì)優(yōu)劣及攪拌反應(yīng)器混合效率的重要參數(shù)[21]。

本文采用軸上扭矩法[18]計(jì)算攪拌槳功耗。首先由扭矩傳感器測(cè)定扭矩M，由電動(dòng)機(jī)得到轉(zhuǎn)速N，再由式(5)～式(7)計(jì)算功耗

計(jì)算發(fā)現(xiàn)，在120 r·min?1下，組合槳功耗為2.65 kW·m?3，而剛性槳?jiǎng)t為2.85 kW·m?3，兩者相比，組合槳功耗下降了7%。這表明組合槳因其柔性尖端與周圍流體的相互作用，能夠有效地將槳葉尖端的能量傳遞給流體，進(jìn)而提高流體的混合性能。這一結(jié)果與劉作華等[18]在柔性槳強(qiáng)化高黏度流體混合的能效分析中所得到的結(jié)論相同。

圖10 體系切向速度分布云圖Fig.10 Tangential velocity contour for systems

圖11 體系混合時(shí)間（120 r·min?1）Fig.11 Mixing time for systems(120 r·min?1)

一種混合時(shí)間的研究方法是向槽內(nèi)加入一種與液體發(fā)生有顏色變化的示蹤劑，通過(guò)觀測(cè)或攝像記錄顏色變化的整個(gè)過(guò)程[22-23]。該方法現(xiàn)象明顯、操作簡(jiǎn)單，能直觀反映攪拌槽中流場(chǎng)的混沌區(qū)和規(guī)則區(qū)的演變情況，廣泛應(yīng)用于確定流體混合效果的實(shí)驗(yàn)中[22]。結(jié)合碘液脫色法[19]，本文對(duì)剛性槳和組合槳在120 r·min?1時(shí)體系的混合時(shí)間進(jìn)行對(duì)比分析，并用相機(jī)記錄槽內(nèi)流體混合的整個(gè)演變過(guò)程。

如圖11所示，對(duì)于剛性槳和組合槳體系來(lái)說(shuō)，流場(chǎng)的規(guī)則區(qū)主要存在于槳葉的上下兩側(cè)，槳葉的下側(cè)位置（靠近槽底）是阻礙流體混合的主要區(qū)域。組合槳通過(guò)自身的多體運(yùn)動(dòng)及形變對(duì)流場(chǎng)的規(guī)則區(qū)進(jìn)行“擾動(dòng)”，這些擾動(dòng)使規(guī)則區(qū)的大小和位置不斷調(diào)整，可增大流場(chǎng)的混沌混合區(qū)[19]。對(duì)比發(fā)現(xiàn)，在tm=10 s時(shí)，組合槳體系中槳葉上側(cè)的規(guī)則區(qū)已逐漸消失，而剛性槳體系的槳葉上側(cè)規(guī)則區(qū)依然存在；在tm=15 s時(shí)，組合槳體系流體已經(jīng)混合均勻，而剛性槳體系在tm=22 s時(shí)，流體才混合均勻，組合槳使體系的混合時(shí)間縮短了近32%。

由數(shù)值計(jì)算結(jié)果可知，剛性槳體系的能量主要集中在槳葉尖端，槽內(nèi)流體的軸向流速很小，特別是最底部區(qū)域流體的流速接近0；而組合槳通過(guò)其柔性尖端對(duì)周圍流體的“擾動(dòng)”和“波動(dòng)”，可將能量從槳葉尖端擴(kuò)散至全槽，使槽內(nèi)流體軸向流動(dòng)加強(qiáng)且切向速度分布得到明顯改善，可提高流體的混合效率。從混合實(shí)驗(yàn)結(jié)果可得，組合槳因其柔性尖端具有高彈性，在力的耦合作用下，增加了攪拌槳對(duì)流體的作用力，流體流動(dòng)加強(qiáng)，使體系的混合時(shí)間縮短了32%；且組合槳可以強(qiáng)化能量傳遞，其攪拌槳功耗降低了7%。這表明模擬結(jié)果和實(shí)驗(yàn)現(xiàn)象基本一致，二者均能反映組合槳更有利于流體的高效混合。

3 結(jié) 論

剛?cè)峤M合槳槳葉尖端的變形量是剛性槳的105倍，其應(yīng)力比剛性槳增加了83%，它可通過(guò)剛-柔-流耦合運(yùn)動(dòng)及多體運(yùn)動(dòng)行為傳遞能量，增加槳葉與流體之間的相互作用力，進(jìn)而使體系的混合時(shí)間縮短近32%，攪拌槳功耗下降7%，可提高流體的混合效率。

符 號(hào) 說(shuō) 明

a——加速度矢量

C——攪拌槳安裝高度，mm

D——攪拌槳直徑，mm

d——位移，m

F——體積力矢量

H——液面高度，mm

l1——槳葉剛性部分長(zhǎng)度，mm

l2——槳葉柔性部分長(zhǎng)度，mm

M——扭矩，N·m

m1——長(zhǎng)槳葉寬度，mm

m2——短槳葉寬度，mm

N——轉(zhuǎn)速，r·min?1

n——矩陣矢量

P——輸入功率，W

PV——功率耗散，kW·m?3

T——攪拌槽直徑，mm

t——時(shí)間，s

tm——混合時(shí)間，s

V——有效體積，m3

ν——速度矢量

δ——剪切力張量

θ——槳葉傾角，（°）

ρ——密度，kg·m?3

σ——柯西應(yīng)力張量

τ——應(yīng)力，Pa

下角標(biāo)

f——流體

s——固體

[1] Jin Zhaowen (靳兆文), Pan Jiazhen (潘家禎). New mixing impeller research [J].Chemical Equipment and Technology(化工設(shè)備技術(shù)), 2004, 25(4): 10-13.

[2] Liu Baoqing (劉寶慶), Zhang Yikun (張義堃), Liu Jingliang (劉景亮), Qin Fulei (秦福磊), Lin Xinghua (林興華), Jin Zhijiang (金志江). Numerical simulation of power consumption and mixing characteristic in stirred vessel with novel coaxial mixer [J].CIESC Journal(化工學(xué)報(bào)), 2013, 64(4): 1135-1144.

[3] Wiggins S. Coherent structures and chaotic advection in three dimensions coherent structures and chaotic advection in three dimensions [J].Journal of Fluid Mechanics, 2010, 654(10): 1-4.

[4] Shekhar C, Nishino K, Yamane Y, Huang J. Stereo-PIV measurement of turbulence characteristics in a flow mixer [J].Journal of Visualization, 2012, 15(4): 293-308.

[5] Liu Fangfang (劉芳芳), Yang Canjun (楊燦軍), Su Qi (蘇琦), Wang Donghai (王東海), Zhang Yizhong (張譯中). Simulation analysis and experimental research on the movements of biomimetic fin [J].Journal of Mechanical Engineering(機(jī)械工程學(xué)報(bào)), 2010, 46(19): 24-29.

[6] Zhang Chengbin (張程賓), Zhao Muwen (趙沐雯), Chen Yongping (陳永平), Shi Mingheng (施明恒). Effects of fluid density on velocity slip in nanochannels [J].CIESC Journal(化工學(xué)報(bào)), 2012, 63(S1): 12-16.

[7] Lee J, Lee S. Fluid-structure interaction for the propulsive velocity of a flapping flexible plate at low Reynolds number [J].Computers and Fluids, 2013, 71(1): 348-374.

[8] Gomes J P, Yight S, Lienhart H, Schafer M. Experimental and numerical study on a laminar fluid-structure interaction reference test case [J].Journal of Fluids and Structures, 2011, 27(1): 43-61.

[9] Antoci C, Gallati M, Sibilla S. Numerical simulation of fluid-structure interaction by SPH [J].Computers and Structures, 2007, 85(1): 879-890.

[10] Martino E S, Guadagni G, Fumero A, Ballerini G, Spirto R, Biglioli P, Redaelli A. Fluid-structure interaction within realastic three-dimensional models of the aneurysmatic aorta as a guidance to assess the risk of rupture of the aneurysm [J].Med. Eng. Phys., 2001, 23(9): 647-655.

[11] Munch C, Ausoni P, Braun O, Farhat M, Avellan F. Fluid-structure coupling for an oscillating hydrofoil [J].Journal of Fluids and Structures, 2010, 26(6): 1018-1033.

[12] Gao F, Guo Z, Sakamoto M, Matsuzawa T. Fluid-structure interaction within a layered aortic arch model [J].J. Biol. Phys., 2006, 32(5): 435-454.

[13] Young Y L. Fluid-structure interaction analysis of flexible composite marine propellers [J].Journal of Fluids and Structures, 2008, 24(1): 799-818.

[14] Young Y L. Time-dependent hydroelastic analysis of cavitating propulsors [J].Journal of Fluids and Structures, 2003, 23(2): 269-295.

[15] Bucchignani E, Stella F, Paglia F. A partition method for the solution of a coupled liquid-structure interaction problem [J].Applied Numerical Mathematics, 2004, 51(4): 463-475.

[16] Karry S, Driss Z, Kchaou H, Abid S. Numerical simulation of fluid-structure interaction in a stirred vessel equipped with an anchor impeller [J].Journal of Mechanical Science and Technology, 2011, 25(7): 1749-1760.

[17] Huang Tiancheng (黃天成), Wang Deguo (王德國(guó)), Zhou Sizhu (周思柱), Wu Hanchuan (吳漢川), Zhu Keqiang (祝克強(qiáng)). Modal research of fluid-solid coupling for the blender truck mixing impeller [J].Journal of Southwest Petroleum University:Science & Technology Edition(西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版), 2012, 34(1): 165-170.

[18] Liu Zuohua (劉作華), Zeng Qiqin (曾啟琴), Wang Yundong (王運(yùn)東), Liu Renlong (劉仁龍), Tao Changyuan (陶長(zhǎng)元). Energy efficiency analysis for high-viscosity fluid mixing enhanced by flexible impeller [J].CIESC Journal(化工學(xué)報(bào)), 2013, 64(10): 3620-3625.

[19] Liu Zuohua (劉作華), Chen Chao (陳超), Liu Renlong (劉仁龍), Tao Changyuan (陶長(zhǎng)元), Wang Yundong (王運(yùn)東). Chaotic mixing enhanced by rigid-flexible impeller in stirred vessel [J].CIESC Journal(化工學(xué)報(bào)), 2014, 65(1): 61-70.

[20] Liu Zuohua (劉作華), Zeng Qiqin (曾啟琴), Yang Xianyan (楊鮮艷), Liu Renlong (劉仁龍), Wang Yundong (王運(yùn)東), Tao Changyuan (陶長(zhǎng)元). Comparative study on flow field structure controlling between rigid-flexible combination impeller and rigid impeller [J].CIESC Journal(化工學(xué)報(bào)), 2014, 65(6): 2078-2084.

[21] Huang Nannan(黃男男). Analysis on flow field and structural optimization of fermentation stirred tank [D]. Wuxi: Jiangnan University, 2009.

[22] Shi Litian (施力田), Gao Zhengming (高正明), Min Jian (閔健). Large eddy simulation for mixing time in stirred tank with dual Rushton turbines [J].CIESC Journal(化工學(xué)報(bào)), 2010, 61(7): 1747-1752.

[23] Zhou Guozhong (周國(guó)忠), Wang Yingchen (王英琛), Shi Litian (施力田). CFD study of mixing process in stirred tank [J].CIESC Journal(化工學(xué)報(bào)), 2003, 54(7): 886-890.

Fluid-structure interaction in liquid mixing intensified by flexible-rigid impeller

ZHU Jun1, ZHOU Zhenglin1, LIU Zuohua1,2, ZHENG Xiongpan1, LIU Renlong1, TAO Changyuan1, WANG Yundong2
(1School of Chemistry and Chemical Engineering,Chongqing University,Chongqing400044,China;2Department of Chemical Engineering,Tsinghua University,Beijing100084,China)

Traditional rigid impeller transfers energy by shearing action, while flexible-rigid impeller can intensify energy transfer by multiple-body movement. Based on the interaction between impeller and fluid, the equivalent stress and total deformation are computationally simulated for flexible-rigid and rigid impellers. Macroscopic flow structure is obtained by two-way fluid-structure interaction technique with simulation platform ANSYS Workbench, and the mixing effect in two stirred systems are discussed with measured mixing time and calculated impeller power dissipation. Compared with rigid impeller, the mixing time of flexible-rigid impeller system is decreased by 32%, its power dissipation is declined by 7%, its total deformation of blade tip is 105times larger, and its equivalent stress of blade tip is 83% greater, so that flexible-rigid impeller exerts greater force to fluid with the fluid-structure interaction, which contributes to energy transmission, liquid flow and mixing intensification.

mixing; fluid; transfer; flexible-rigid impeller; ANSYS Workbench; fluid-structure interaction

Prof. LIU Zuohua, liuzuohua@cqu.edu.cn

10.11949/j.issn.0438-1157.20141571

TQ 027.2

：A

：0438—1157（2015）10—3849—08

2014-10-22收到初稿，2015-04-20收到修改稿。

聯(lián)系人：劉作華。

：朱俊（1982—），男，博士研究生，高級(jí)工程師。

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2012CBA01203）；清華大學(xué)化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放課題（SKL-ChE-12A02）；重慶市自然科學(xué)基金重點(diǎn)項(xiàng)目（CSTC2012JJB0006）。

Received date: 2014-10-22.

Foundation item: supported by the National Basic Research Program of China (2012CBA01203), the State Key Laboratory of Chemical Engineering (SKL-ChE-12A02) and the Key Natural Science Fund of Chongqing (CSTC2012JJB0006).