基于ANSYS CFX的PVA纖維束攪拌分散作用機理研究

曹源文，周 華，肖 偉，楊 雪，夏柱林，鄭南翔

(1. 重慶交通大學 機電與車輛工程學院，重慶 400074； 2. 濟祁高速公路(碭山段)項目辦公室，安徽 宿州 235300；3. 長安大學 公路工程學院，陜西 西安 710064)

0 引 言

現(xiàn)代工程水泥混凝土中通過摻入聚乙烯醇纖維(PVA)來改善水泥混凝土的抗裂性能和水泥基復合材料的變形及破壞特性[1]。但在PVA纖維摻入水泥混凝土進行拌和之前，需將PVA纖維束分散成單絲狀并形成均勻的三維網(wǎng)狀結(jié)構(gòu)，使分散后的PVA纖維與水泥混凝土材料拌和均勻，發(fā)揮出其對材料的改善作用[2-3]，否則成束的纖維滲入水泥混凝土中也難以明顯提高其材料及工作性能。

針對某款PVA纖維束攪拌分散設(shè)備開展研究，通過建立該設(shè)備的數(shù)值模擬模型，設(shè)定參數(shù)并求解后分析仿真結(jié)果，觀察并分析不同攪拌轉(zhuǎn)速和攪拌槳類型下PVA纖維束的分散性能，得出攪拌設(shè)備對PVA纖維束的攪拌作用機理，確定能使分散效果呈現(xiàn)最佳狀態(tài)的轉(zhuǎn)速和攪拌槳類型，為進一步研究PVA纖維束分散技術(shù)提供可靠依據(jù)。

1 攪拌設(shè)備模型

1.1 攪拌設(shè)備基本結(jié)構(gòu)及工作原理

采用的PVA纖維束攪拌設(shè)備為一款立式精密電動攪拌設(shè)備。設(shè)備主要由8個部分組成：直流電機、調(diào)速器、攪拌軸、聯(lián)軸器、攪拌槳、攪拌槽、支架和底座。攪拌設(shè)備結(jié)構(gòu)如圖1。攪拌槳采用四直葉槳、四斜葉槳及其組合槳。其中，放置于攪拌槽中的攪拌物質(zhì)為PVA纖維束和HEC分散劑水溶液的混合物，所采用的纖維長度為24 mm、直徑為12 μm、拉伸強度為1 600 MPa、楊氏模量為35 GPa[4]。實際的攪拌分散過程極其復雜，主要是在對流混合、擴散混合和剪切混合這3類攪拌運動的共同作用下達到PVA纖維束的分散和均勻[5-6]。不同攪拌設(shè)備的主要流場形態(tài)不同，可分為徑向流、軸向流和切向流，其中徑向流和軸向流能對混合和剪切起到關(guān)鍵作用，促進攪拌分散，而切向流則應(yīng)該加以約束[7-8]。

圖1 攪拌設(shè)備結(jié)構(gòu)Fig. 1 Mixing equipment structure

1.2 攪拌槳幾何模型建立

采用四直葉槳、四斜葉槳及其組合槳進行研究，攪拌槳葉輪直徑均為70 mm，組合槳中上下層間距25 mm。在CATIA建立模型并通過布爾操作運算得到計算域的幾何模型，如圖2。

圖2 攪拌槳三維實體模型Fig. 2 Three dimensional solid model of stirring blade

2 流場數(shù)學模型建立

2.1 湍流數(shù)學模型

選取的RNGk-ε湍流模型是在標準k-ε模型基礎(chǔ)上的修正方程[9]。

標準k-ε模型湍流動能k方程為：

(1)

標準k-ε模型湍流動能耗散率ε方程為：

(2)

式中：σk、Cε、Cε1、Cε2為常數(shù)；μt為湍流黏度，其值與湍流動能和湍流動能耗散率有關(guān)；Pk是黏性力和浮力的湍流產(chǎn)物。μt和Pk的方程為：

(3)

Pk=μt

(4)

RNGk-ε模型湍流生成及湍流耗散方程同標準k-ε模型是一樣的，但是它的常數(shù)系數(shù)Cε替換為CεRNG，計算精度提高且適用范圍更廣[10]。

RNGk-ε模型湍流動能耗散方程為：

(5)

2.2 求解設(shè)定

采用八叉樹法生成四面體網(wǎng)格，面網(wǎng)格采用默認設(shè)置的按照面輪廓線生成四邊形和小部分三角形網(wǎng)格。在CFX-Pre中設(shè)置PVA纖維束的密度ρ1=1.3 g/cm3，動力粘度μ1=22 Pa·s,以及溶解度為15%的HEC分散劑混合水溶液的密度ρmHEC=0.95 g/cm3。設(shè)定計算域的靜止域和旋轉(zhuǎn)域參數(shù)，在靜止域上設(shè)置進出口邊界、壁面邊界，旋轉(zhuǎn)域上設(shè)置壁面邊界，在旋轉(zhuǎn)域和靜止域的交界面設(shè)置interface對[11]。求解控制參數(shù)設(shè)置為：求解格式為高階求解模式，設(shè)置最大迭代步數(shù)200步，選擇自動時間尺度，收斂判據(jù)標準選擇均方根殘差值(RMS)，殘差值為10-4。

在靜止域壁面設(shè)定時，使用默認的無滑移壁面No Slip Wall；在旋轉(zhuǎn)域的壁面設(shè)定時，設(shè)定為旋轉(zhuǎn)壁面，旋轉(zhuǎn)角速度為0。對于無特別指定的外部區(qū)域，系統(tǒng)都自動指定為無滑移、絕熱的壁面條件，無特別指出的內(nèi)部區(qū)域?qū)⒈缓雎訹12]。壁面條件設(shè)置如圖3。

圖3 靜止域和旋轉(zhuǎn)域壁面設(shè)置Fig. 3 Wall settings in static and rotating domains

由于CFX-Pre的材料庫中沒有PVA纖維束和分散劑水溶液，因此在計算域設(shè)定之前，在庫設(shè)定的材料設(shè)定中自定義生成PVA和dispersant solution兩種材料。材料屬性定義界面如圖4。

圖4 材料屬性定義界面Fig. 4 Definition interface of material properties

3 攪拌作用機理分析

3.1 不同轉(zhuǎn)速的作用機理分析

選取四直葉攪拌槳，添加HEC分散劑。攪拌槳轉(zhuǎn)速范圍為500 ～800 r/min，以50 r/min為間隔，得到不同攪拌轉(zhuǎn)速下流場的速度和湍流動能分布。

3.1.1 PVA纖維速度云圖分布

在速度梯度范圍為0～2.021 m/s，轉(zhuǎn)速范圍為500～800 r/min下的速度云圖如圖5。

圖5 不同攪拌轉(zhuǎn)速的速度云圖Fig. 5 Velocity contour of different stirring speed

由速度云圖可以看出：

1)攪拌槽內(nèi)PVA纖維與分散劑水溶液混合物作強制對流運動，且隨著轉(zhuǎn)速的增加，由葉輪旋轉(zhuǎn)的機械能轉(zhuǎn)化的混合物對流循環(huán)動能增加，剪切作用增強，功耗也隨之增加。葉端線速度最快，剪切力最強，沿葉輪半徑方向的出口流速最高，對纖維的破壞力最大，整體的徑向循環(huán)液流使PVA纖維束得到分散與混合。

2)轉(zhuǎn)速為500 r/min時，槽中靜止域內(nèi)基本為切向流主導的周向水平低速環(huán)流，葉端剪切力小，液流速度約為1.2 m/s。同時，葉輪間存在靠近中心的單環(huán)渦流和靠近槽壁的回流，其中槽壁回流速度較渦流略高，形成擴散混合。另外，由于PVA纖維的纏敷性，單環(huán)渦流易使纖維在攪拌軸上纏繞結(jié)團，極不利于分散。當轉(zhuǎn)速提高到550 r/min和600 r/min時，液流中心盲區(qū)(流體速度很小或趨于零的區(qū)域，過大則對分散不利)較500 r/min時有所減小，槽壁回流區(qū)增大，從而在葉端外靠近槽壁處形成低速渦流區(qū)，對分散不利。旋轉(zhuǎn)域與靜止域的各個部分速度差異較大，分布不均勻，影響纖維的分散與混合。

3)轉(zhuǎn)速為650 r/min時，葉端剪切作用增強，液流速度達到1.8 m/s，中心盲區(qū)進一步減小，葉片間單環(huán)渦流消失，槽壁回流融入對流混合運動，葉端外低速渦流區(qū)幾乎消失，速度分布均勻程度高，靜止域與旋轉(zhuǎn)域形成了宏觀整體對流循環(huán)的混合運動，促進了纖維的剪切與混合。當轉(zhuǎn)速繼續(xù)加大時，雖然葉端剪切作用增強，但槽壁回流間的低速渦流區(qū)在增大，相對速度差異較650 r/min時更加明顯。過高轉(zhuǎn)速使得整體對流循環(huán)運動開始分解，不利于達到分散均勻的效果。同時，中心盲區(qū)幾乎無變化且設(shè)備功耗增加。

3.1.2 不同轉(zhuǎn)速對湍流動能的影響

圖6是PVA纖維水溶液的湍流動能梯度范圍為1.505～20.4 m2/s2、轉(zhuǎn)速范圍為500～800 r/min的湍流動能云圖。

圖6 不同攪拌轉(zhuǎn)速的湍流動能云圖Fig. 6 Turbulent kinetic energy contour of different stirring speed

由圖6可以看出：

1)隨著轉(zhuǎn)速的增加，由機械能轉(zhuǎn)化的湍流動能增加，液流運動越來越劇烈且紊亂，對流混合對PVA纖維的混合作用增強。同時，葉端附近湍流擴散混合和剪切混合越來越劇烈，范圍逐漸擴大。

2)轉(zhuǎn)速為500 r/min時，葉端處湍流動能很小，約為9 m2/s2，靜止域葉輪間幾乎全部為低湍能區(qū)，混合效果差。同時，該區(qū)域的周向水平切向流與旋轉(zhuǎn)域內(nèi)周向低湍能區(qū)的剪切與混合作用弱，整體上沒有形成較為明顯的徑向流和對流混合。當轉(zhuǎn)速提高到550 r/min和600 r/min時，中心盲區(qū)及槽壁低湍能區(qū)有所減小，但葉端處湍流動能仍然較小，葉端外低速渦流區(qū)仍然較大，影響分散。

3)轉(zhuǎn)速為650 r/min和700 r/min時，葉端湍流動能約為19 m2/s2，剪切力增大且擴散混合增強，形成了比較劇烈的徑向流和整體對流混合。同時，中心盲區(qū)及槽壁低湍能區(qū)進一步減小，旋轉(zhuǎn)域與大部分靜止域的湍流動能梯度減小，動能分布較均勻。加大轉(zhuǎn)速后，靜止域的低湍能區(qū)大小及分布幾乎無變化，旋轉(zhuǎn)域內(nèi)則有少許增大，但湍流動能梯度擴大，影響混合效果。同時，對流混合開始分解且功耗上升。

綜上所述，低轉(zhuǎn)速時剪切與混合作用較差，隨著轉(zhuǎn)速的增大，攪拌槽內(nèi)PVA纖維混合液速度和湍流動能都隨之增大。但是在轉(zhuǎn)速大于約650 r/min后，速度和湍流動能的增大梯度都有所減緩，總體循環(huán)與分散效果變差，同時攪拌設(shè)備功耗變大。在滿足PVA纖維束分散效果的前提下，考慮攪拌設(shè)備的經(jīng)濟節(jié)能性，認為轉(zhuǎn)速650 r/min是PVA纖維束攪拌設(shè)備的最佳攪拌轉(zhuǎn)速。

3.2 不同攪拌槳的作用機理分析

選取轉(zhuǎn)速為650 r/min的3種不同攪拌槳進行CFD分析，通過不同攪拌槳的速度云圖和湍流動能云圖來研究其流場特性，從而分析不同攪拌槳的攪拌分散作用機理。

3.2.1 PVA纖維速度云圖分布

圖7是速度梯度范圍為0～2.021 m/s、轉(zhuǎn)速為650 r/min下3種攪拌槳的PVA纖維縱截面速度云圖。

圖7 不同攪拌槳的速度云圖Fig. 7 Velocity contour of different stirring blades

由速度云圖可以看出：

1)混合溶液在3種攪拌槳的旋轉(zhuǎn)攪拌下，均以大范圍的對流混合為主，但各自的速度大小及分布不同，產(chǎn)生的流場差異較大。三者均在葉端處產(chǎn)生最高速度，即剪切力最強，以葉端為中心的向外液流表現(xiàn)出不同的混合運動及流場形態(tài)。

2)同一轉(zhuǎn)速下，四直葉槳葉端處速度相對較小，剪切力較弱。四直葉槳屬于徑流型攪拌設(shè)備，主要流場形態(tài)為徑向流，包含以槳葉水平面為中心分界面的上下兩個徑向循環(huán)回路。下循環(huán)回路自葉端射出，沿槽壁向下再徑向流回槳葉下方；上循環(huán)回路則對稱相反。由于受重力作用影響，下方回路流速較高。中心盲區(qū)較小，主要分布在攪拌槳中心的上下區(qū)域，總體速度分布較均勻。四斜葉槳由于其強烈的軸向推動作用，在葉端產(chǎn)生相對較高的軸向速度，但其為垂直分流，剪切力較小。中心盲區(qū)分布于攪拌槽上部和攪拌槳中心以下區(qū)域，范圍很大。槽底旋轉(zhuǎn)域外分布大范圍的軸向混合液流，但其速度較低。總體來說，剪切作用弱，混合效果一般。

3)組合槳盲區(qū)較四直葉槳變化不大，較四斜葉槳大幅度減少。由于上方四斜葉產(chǎn)生的軸向循環(huán)流減少了混合液下部液流的壓力，下方四直葉槳葉端附近的液流速度較單四直葉更大，上下循環(huán)回路液流速度更高、范圍更廣，剪切及擴散混合作用增強。四斜葉槳旋轉(zhuǎn)域和靜止域由于受到四直葉槳向上射流的疊加影響，對流循環(huán)強度較單四斜葉更大?？傮w來說，組合槳的剪切與混合作用較前兩者更好。

3.2.2 不同攪拌槳的湍流動能分布

圖8為縱截面內(nèi)不同攪拌槳湍流動能云圖，注意湍流動能梯度不同?？梢钥闯觯?/p>

1)3種攪拌槳均在葉端附近產(chǎn)生最大湍流動能，組合槳能夠產(chǎn)生相對范圍較大的最高湍能，四直葉槳次之且范圍較小，而四斜葉槳湍能最小。3者在以對流混合為主的液流作用下，湍能分布也各不相同。

2)四直葉槳的高速徑向流自葉端射出后，遇槽壁后形成向上和向下兩個徑向回流，葉端處形成擴散混合，向外則是湍能梯度較小的對流混合，中心盲區(qū)與旋轉(zhuǎn)域內(nèi)徑向回流的湍能相近，整體分布較均勻，PVA纖維能夠得到一定的剪切和分散混合。四斜葉槳的曲面構(gòu)造對液流具有良好的推進作用并形成軸向流，當液流回到葉輪下方處，槳葉對其向上推送并產(chǎn)生輕微的剪切作用，因此，四斜葉槳排液量大但剪切力小。

3)組合槳中，由于四直葉槳離底間隙較小，徑向方向上槳葉的向下射流與槽底撞擊速度大，該區(qū)域湍流劇烈并帶動周圍的纖維運動，葉端附近產(chǎn)生較大范圍的高湍能區(qū)，并以葉端為中心向外呈均勻梯度減少，湍能較單四直葉槳更大且范圍更廣。四斜葉槳下方受四直葉槳向上射流的影響，該區(qū)域較單四斜葉槳擴散混合增強，同時由于軸向流作用其槳葉上方形成范圍很大的對流混合，總體對流循環(huán)的湍能分布最佳。

圖8 不同攪拌槳的湍流動能云圖Fig. 8 Turbulent kinetic energy contour of different stirring blades

綜上所述，通過分析不同攪拌槳的PVA纖維速度云圖及湍流動能云圖，得出組合槳對PVA纖維的剪切及總體循環(huán)混合效果最好。研究的3種攪拌槳對PVA纖維攪拌分散效果排序為：組合槳>四直葉槳>四斜葉槳。

4 結(jié) 語

以攪拌槳轉(zhuǎn)速和類型為研究對象，應(yīng)用ANSYS CFX對攪拌分散設(shè)備進行了數(shù)值模擬，根據(jù)流場速度云圖和湍流動能云圖對PVA纖維束攪拌作用機理進行分析，結(jié)論如下：

1)為使PVA纖維束進一步分散，在攪拌分散中應(yīng)盡量加強剪切與混合作用，達到纖維束的均勻分散。攪拌槳轉(zhuǎn)速影響PVA 纖維速度和湍流動能；四直葉槳和四斜葉槳能夠產(chǎn)生不同的纖維混合溶液流場形態(tài)，組合槳則會產(chǎn)生更加復雜的混合運動。

2)攪拌轉(zhuǎn)速越大，PVA纖維速度和湍流動能越大。但當轉(zhuǎn)速增加到約650 r/min之后，整個計算域液流循環(huán)效果變差且速度與湍能差異擴大，總體對流循環(huán)及分散效果開始不利于纖維的剪切與混合。綜合考慮，認為轉(zhuǎn)速650 r/min是PVA纖維束攪拌設(shè)備的最佳攪拌轉(zhuǎn)速。

3)四直葉槳屬于徑流型攪拌設(shè)備，具有較好的剪切作用；四斜葉槳則為軸向型，能夠產(chǎn)生良好的整體循環(huán)流動；組合槳在二者的交互作用下，總體對流循環(huán)與分散效果更好，能夠產(chǎn)生更強的剪切與混合作用，可以更有效地促進PVA纖維束的攪拌分散。