應(yīng)用STAR-CCM+的旋流噴嘴內(nèi)部三維流場數(shù)值模擬與分析

陳 娜， 吳 敏， 邱 華， 葛明橋

(生態(tài)紡織教育部重點(diǎn)實(shí)驗(yàn)室(江南大學(xué))， 江蘇 無錫 214122)

隨著紡織科技的不斷進(jìn)步，氣流在紡織生產(chǎn)中的應(yīng)用日益增多。除了運(yùn)用氣流來改變成紗機(jī)制的新技術(shù)外，還可借助氣流來改善紗線質(zhì)量，如集聚紡[1]和旋流噴嘴環(huán)錠紡技術(shù)。旋流噴嘴環(huán)錠紡是在環(huán)錠紡上利用旋流噴嘴降低紗線表面毛羽，改善紗線質(zhì)量的技術(shù)。褚結(jié)[2]、牟俊玲[3]、張英姿[4]等主要研究了旋流噴嘴紗道與氣道直徑、長度和紗道形狀等結(jié)構(gòu)參數(shù)以及紡紗工藝參數(shù)對減羽效果的影響。本文借助計(jì)算流體力學(xué)軟件 STAR-CCM+對旋流噴嘴紗道三維流場進(jìn)行數(shù)值模擬，在前人設(shè)計(jì)的旋流噴嘴的基礎(chǔ)上通過改變氣道與紗道相交方向和相切點(diǎn)的位置，探討其對噴嘴紗道內(nèi)部氣流分布的影響，并選擇出較好的噴嘴結(jié)構(gòu)模型，為今后研究其他參數(shù)對旋流噴嘴減羽效果的影響奠定了基礎(chǔ)。STAR-CCM+是由CD-adapco Group公司開發(fā)的新一代通用計(jì)算流體力學(xué)(CFD)分析軟件[6-8]。該軟件具有一體化集成環(huán)境、高重復(fù)性和適用性等優(yōu)點(diǎn)，可通過一體化用戶界面展示從CAD建模、表面網(wǎng)格準(zhǔn)備、體網(wǎng)格生成、模型設(shè)定、計(jì)算求解一直到后處理分析的整個(gè)模擬過程，使用方便，高效。

1 模型的建立

1.1 旋流噴嘴簡介

旋流噴嘴是為減少環(huán)錠紡紗線毛羽而發(fā)明的組件[2-5]，這種減羽方法是直接在普通環(huán)錠細(xì)紗機(jī)的前羅拉和導(dǎo)紗鉤之間裝一個(gè)旋流噴嘴，噴嘴由氣道和紗道組成，利用噴嘴內(nèi)部特殊結(jié)構(gòu)，可產(chǎn)生高速旋轉(zhuǎn)氣流，從而有效減少紗線表面毛羽，改善紗線性能。圖1示出環(huán)錠紡中使用旋流噴嘴的工藝流程。粗紗經(jīng)過牽伸系統(tǒng)送入到前羅拉，隨后紗線通過旋流噴嘴，表面毛羽減少，最后細(xì)紗在導(dǎo)紗鉤的作用下，經(jīng)過鋼絲圈和鋼領(lǐng)，被卷繞在筒管上，旋流噴嘴的安裝傾角應(yīng)盡量與前羅拉傾角一致。

圖1 環(huán)錠紡中使用旋流噴嘴的工藝流程Fig.1 Ring spinning process of using a swirl nozzle

1.2 物理模型建立

噴嘴物理模型如圖2所示。設(shè)計(jì)的2種噴嘴模型分別是A和B。2個(gè)模型中圓柱①為氣道，圓柱②主要是紗線通過的管道，稱為紗道，氣道與紗道壁相切并呈一定角度，分別為α1和α2，滿足紗道中氣流旋轉(zhuǎn)方向與Z捻紗加捻方向一致。A和B模型中紗道與氣道相切點(diǎn)沿X軸負(fù)方向長度均為7 mm，A噴嘴氣道與紗道相交方向向下，而B相交方向向上。A和B噴嘴的物理模型尺寸與實(shí)際尺寸一致，紗道直徑都為2 mm，氣道直徑為1.4 mm，紗道長度為21 mm。物理模型中忽略了紗線的存在，對噴嘴的結(jié)構(gòu)做了適當(dāng)簡化。

圖2 旋流噴嘴物理模型Fig.2 Physical model of swirl nozzles

2 旋流噴嘴內(nèi)部三維氣流模擬

2.1 網(wǎng)格劃分及物理模型選擇

由于Trimmed網(wǎng)格計(jì)算精度比四面體網(wǎng)格高，而收斂速度又比多面體網(wǎng)格快[7]，故旋流噴嘴網(wǎng)格劃分選擇Trimmed網(wǎng)格，網(wǎng)格基本尺寸為0.1 mm。噴嘴內(nèi)氣流總的流速在330 m/s左右，將氣流流動(dòng)計(jì)算界定為三維定?？蓧嚎s湍流流動(dòng)，噴嘴計(jì)算選擇分離求解器做定常計(jì)算，根據(jù)文獻(xiàn)[3]設(shè)定計(jì)算環(huán)境為2個(gè)大氣壓，不考慮重力的影響及熱交換，湍流采用壁面函數(shù)法。離散格式為二階迎風(fēng)方式。

2.2 邊界條件設(shè)置

在STAR-CCM+中首先要對每一個(gè)邊界進(jìn)行逐一定義，入口邊界為氣道口，出口邊界為紗道2個(gè)進(jìn)出口。由于紗道內(nèi)紗線的質(zhì)量較小，對壓縮流體的運(yùn)動(dòng)影響很小，可以忽略，因此在數(shù)值模擬中，假設(shè)紗道中只有壓縮空氣。此外，流場在氣流進(jìn)入后很短的時(shí)間內(nèi)即可完成，因此可認(rèn)為紗道壁是絕熱的。

2.2.1入口邊界

大氣壓強(qiáng)為0.101 3 MPa，大氣溫度為293 K，儲(chǔ)氣罐內(nèi)空氣的壓強(qiáng)為P，P=0.2 MPa(表壓)，由于(0.101 3+P)/0.101 3>0.528，可知調(diào)壓閥前后壓差比較大，氣流在旋流噴嘴氣道的入口處發(fā)生阻塞，故氣體能流通的流量只與入口壓力有關(guān)(即跟氣體的壓縮程度有關(guān))[9]，因此，氣室中氣體壓力P0=(表壓+1)×1.013×105(Pa)，經(jīng)計(jì)算，入口處氣流壓強(qiáng)P0為0.159 MP；根據(jù)理想氣體狀態(tài)方程ρ0RT0=P0(ρ0為氣室中氣體的密度，kg/m3，氣體常數(shù)R=287.1 N·m/(kg·K))，計(jì)算得到入口處氣流溫度T0為244.069 K。

2.2.2出口邊界和固體壁面

由于旋流噴嘴出口與大氣直接相連，噴嘴出口壓強(qiáng)近似等于外界大氣壓強(qiáng)，即設(shè)置為一個(gè)標(biāo)準(zhǔn)大氣壓，溫度也為大氣溫度。固體壁面處壓縮流體的運(yùn)動(dòng)受到較大阻礙，其速度越來越小，因此近似認(rèn)為壁面處氣流速度在各個(gè)方向均為零，壁面邊界采用無滑移條件。

3 計(jì)算結(jié)果與分析

3.1 管道相交方向?qū)喌纼?nèi)部氣流的影響

3.1.1噴嘴紗道內(nèi)部氣流分布特點(diǎn)

圖3示出2種噴嘴紗道內(nèi)部氣流沿徑向的速度矢量圖。紗道總長為21 mm,氣流分布主要在氣流入口以上6 mm到氣流入口下12 mm段有明顯特征分布，故將紗道以2 mm為單位，模型A和B劃分成10個(gè)平面，規(guī)定噴嘴入口平面為0。2個(gè)模型紗道中的氣流均是順時(shí)針旋轉(zhuǎn)，氣流入口以下，氣流是向下順時(shí)針旋轉(zhuǎn)的，氣流入口以上，則是向上順時(shí)針旋轉(zhuǎn)。A模型由于噴嘴傾角向下傾斜，大部分氣流向下運(yùn)動(dòng)，而氣流入口以上旋轉(zhuǎn)氣流較少；B模型傾角向上傾斜，相反大部分氣流向氣流入口上方運(yùn)動(dòng)，而下方氣流較少。圖4示出2個(gè)噴嘴紗道內(nèi)部氣流沿軸向的速度矢量圖。A和B噴嘴中旋轉(zhuǎn)氣流在噴嘴軸向形成不同程度的回流區(qū)，回流區(qū)基本在氣流入口附近，其中A模型在氣流入口上下均存在回流區(qū)，而B相對較少，只存在于氣流入口的下方區(qū)域。由于回流區(qū)氣流比較紊亂，會(huì)使紗線表面毛羽相互糾纏，形成棉結(jié)造成紗線條干不勻，應(yīng)盡可能避免。

圖3 噴嘴模型A和B內(nèi)部氣流沿徑向的速度矢量圖Fig.3 Velocity vector diagram along radial profiles in nozzle model A and B

圖4 噴嘴模型A和B內(nèi)部氣流沿軸向的速度矢量圖Fig.4 Velocity vector diagram along axial profiles in nozzle model A and B

由于紡紗加捻三角區(qū)的存在，使纖維尾端未能全部加捻到紗體內(nèi)，從而形成較多的毛羽，其中48.6%的毛羽分布如圖5中A所示，定義該毛羽為順向毛羽，而29.6%的毛羽分布方向與順向毛羽相反，定義為逆向毛羽，此外還有18.4%的雙向毛羽和3.4%的其他毛羽[10]。由此看出，當(dāng)紗道中的氣流軸向運(yùn)動(dòng)與紗線運(yùn)動(dòng)方向相反時(shí)，可使紗線表面順向毛羽更貼附于紗體，有助于減少較多的該向毛羽；相反，氣流與紗線運(yùn)動(dòng)方向相同時(shí)，能夠減少少量的逆向毛羽，但是不能使順向毛羽順著氣流方向前進(jìn)，從而被反向吹散，形成新形態(tài)的毛羽，有些則直接被吹出紗體，造成纖維的散失，因此要獲得較好的減羽效果，紗道中應(yīng)形成較多的與紗線軸向運(yùn)動(dòng)方向相反的旋轉(zhuǎn)氣流，由于該氣流主要分布在氣流入口以上，故紗道該部分的長度要適中。

圖5 氣流對紗線表面毛羽作用過程Fig.5 Effects of air flow on the surface of the yarn hairiness

3.1.2噴嘴紗道內(nèi)部沿軸向氣流三向速度分布

從圖5可看出，氣流進(jìn)入紗道后對紗體的作用力可分解為軸向、徑向和切向。切向氣流主要對紗線表面毛羽起包纏作用，切速增大有利于對毛羽的包纏；徑向和軸向氣流主要可以使紗線表面浮游纖維易于伏貼在紗線表面，從而輔助切向氣流使毛羽很好地包纏在紗體上[11]，軸速和徑速增大，有利于紗線表面毛羽的貼伏和彎曲，但過大會(huì)使紗線表面毛羽發(fā)生折疊、糾纏，不利于減羽作用。

圖6示出A和B噴嘴模型的三向速度分量。圖6(a)示出噴嘴A和B的切向氣流分量，為了使紗線表面毛羽盡可能多地包纏到紗體上，在氣流入口以上(X負(fù)半軸)，具有較大的切向速度，有利于對伸出紗體纖維端的良好包纏；氣流入口以下(X正半軸)，切速較小且相對穩(wěn)定，可以保持纖維良好的包纏狀態(tài)。從圖6(b)可看出，噴嘴B的徑向速度與切向速度具有相似的變化趨勢，這也有助于徑向速度配合切向氣流實(shí)現(xiàn)對紗線表面的良好包纏。圖6(c)示出噴嘴的軸向速度分量，結(jié)合圖4，在氣流入口以上，軸速為負(fù)，而在氣流入口以下，軸速為正。A噴嘴軸速在氣流入口附近速度與正常流速方向相反且較大，所以在該區(qū)域產(chǎn)生較大的回流區(qū)，并且在氣流入口以下軸速較大，不利于減少順向毛羽；而B噴嘴軸速在氣流入口以下較小且比較穩(wěn)定，在氣流入口附近也存在回流區(qū)，但比較小。結(jié)合以上幾點(diǎn)，由于B噴嘴在紗道中形成較多與紗線軸向運(yùn)動(dòng)方向相反的氣流，紗道中三向速度分布也比較合理且紗道中回流區(qū)較少，因此選擇B噴嘴模型來探討氣道與紗道相切點(diǎn)位置對噴嘴紗道內(nèi)部氣流分布規(guī)律的影響。

圖6 A和B噴嘴模型的三向速度分量Fig.6 Three velocity components in A and B nozzle model. (a) Tangential velocity component; (b) Radial velocity component; (c) Axial velocity component

3.2 管道相切點(diǎn)位置對紗道內(nèi)部氣流影響

通過對A、B模型中氣流分布規(guī)律的分析，結(jié)合氣流對毛羽的作用機(jī)制確定B噴嘴為合理模型，為研究噴嘴位置對噴嘴紗道內(nèi)部氣流的影響，設(shè)計(jì)了C和D噴嘴。C和D 噴嘴的其他參數(shù)與B噴嘴一致，只是C噴嘴中氣道與紗道相切點(diǎn)沿X軸負(fù)向長度為10.5 mm，而D噴嘴為14 mm。

3.2.1噴嘴紗道內(nèi)部氣流分布特點(diǎn)

圖7、8分別為B、C、D 3個(gè)噴嘴紗道內(nèi)部氣流沿徑向和軸向的速度矢量圖，規(guī)定氣流入口平面為X=0平面。3個(gè)噴嘴在X=0平面以上，氣流運(yùn)動(dòng)方向一致，只是隨著氣流入口以上長度的增加，紗道中向上運(yùn)動(dòng)的氣流增加，這正滿足了減羽的第一個(gè)要素，形成較多的與紗線沿軸向運(yùn)動(dòng)方向相反的旋轉(zhuǎn)氣流，從該點(diǎn)來看D噴嘴比較合理；同時(shí)，當(dāng)噴嘴的負(fù)半軸較短時(shí)，在紗道入口形成的氣流流速較大且紗道入口離紡紗三角區(qū)較近，易造成該區(qū)域纖維的散失，使紗線細(xì)度達(dá)不到要求；此外，3個(gè)模型均在噴嘴下方有少量的回流區(qū)，B和C由于氣流入口以下區(qū)域較長，使紗道出口處氣流分布較少，所以氣流運(yùn)動(dòng)方向較紊亂，而D噴嘴氣流在該區(qū)域分布相對來說比較穩(wěn)定。

圖7 3個(gè)噴嘴模型紗道內(nèi)部氣流沿徑向的速度矢量圖Fig.7 Velocity vector diagram along radial profile in three nozzle models

圖8 3個(gè)噴嘴模型紗道內(nèi)部氣流沿軸向的速度矢量圖Fig.8 Velocity vector diagram along axial profile in three nozzle models

3.2.2噴嘴紗道內(nèi)部沿軸向氣流三向速度分布

圖9示出B、C、D噴嘴紗道中氣流的三向速度分量。

圖9 B、C、D噴嘴模型的三向速度Fig.9 Three velocity components in B, C and D nozzle models. (a) Tangential velocity component； (b) Radial velocity component；(c) Axial velocity component

由圖9可看出，3個(gè)噴嘴的三向速度分布規(guī)律基本一致，但D噴嘴的切向速度和徑向速度在氣流入口以上比其余2個(gè)噴嘴速度大，而在氣流入口以下則最小，這正滿足了主要減羽區(qū)對速度的要求；從3個(gè)噴嘴紗道軸向速度分布規(guī)律中可以看出，D噴嘴在紗道入口附近軸速為正(與氣流正常運(yùn)動(dòng)方向相反)，說明該區(qū)域有回流產(chǎn)生，主要是由于噴嘴負(fù)半軸較長，在紗道入口附近氣流速度比較小，導(dǎo)致外部氣流流入噴嘴紗道中，3個(gè)噴嘴在X正半軸氣流入口附近(X=0附近)均有少量回流產(chǎn)生，而D噴嘴回流區(qū)域最小且在該區(qū)域軸速最小。故由紗道中氣流的三向速度分布規(guī)律驗(yàn)證了D噴嘴模型相對比較合理。

4 結(jié) 論

使用計(jì)算流體軟件STAR-CCM+對設(shè)計(jì)的旋流噴嘴紗道內(nèi)部氣流流場進(jìn)行了三維數(shù)值模擬，得到紗道內(nèi)氣流的速度矢量圖和三向速度分布規(guī)律。

1)模型A和B是2個(gè)氣道與紗道相交方向不同的噴嘴，從減羽效果來看，氣流運(yùn)動(dòng)方向與紗線軸向運(yùn)動(dòng)相反時(shí)有助于減少紗線表面較多的順向毛羽，同時(shí)在氣流入口以下區(qū)域，較小的流速和回流區(qū)對減羽效果有利，因此從以上來看噴嘴B模型比A模型合理。

2)紗道內(nèi)三向流速分布直接影響毛羽在紗體表面的貼附程度和對毛羽的包纏效果，當(dāng)氣流入口離紗線入口為14 mm時(shí)，噴嘴的切向速度和徑向速度在氣流入口以上較大，而在氣流入口以下則較小，且在紗道中的回流區(qū)也比較小，這正滿足了可減少大量順向毛羽的要求。

3)噴嘴紗道中氣流的三維數(shù)值模擬沒有考慮紗線的存在，模擬結(jié)果為噴嘴結(jié)構(gòu)設(shè)計(jì)提供了一定的理論基礎(chǔ)，但最終的減羽效果還需進(jìn)一步的實(shí)驗(yàn)驗(yàn)證。

FZXB

[1] 鄒專勇,竺韻德,俞建勇,等. 鉆孔皮圈集聚紡紗系統(tǒng)三維流場的模擬與分析[J].紡織學(xué)報(bào),2009,30 (4): 111-116.

ZOU Zhuanyong, ZHU Yunde, YU Jianyong, et al. Simulation and analysis of three dimensional flow field in punched apron compact spinning[J]. Journal of Textile Research,2009,30( 4): 111-16.

[2] 褚結(jié),葛明橋.旋流器降低環(huán)錠紡紗線毛羽的研究[J].紡織學(xué)報(bào),2007,28(12): 34-37.

ZHU Jie, GE Mingqiao. Research into reducing the Hairiness of ring spun yarn by a swirler[J]. Journal of Textile Research,2007,28(12): 34-37.

[3] 牟俊玲,邱華,葛明橋.旋流器對環(huán)錠紡紗性能的影響[J]. 紡織學(xué)報(bào),2009,30(11): 43-47.

MOU Junling, QIU Hua, GE Mingqiao. The influence of cyclone setting on yarn properties[J]. Journal of Textile Research,2009,30(11): 43-47.

[4] 張英姿,邱華,吳敏.旋流器紗道截面形狀對環(huán)錠紡紗線質(zhì)量的影響[J].紡織學(xué)報(bào)，2011,32(7): 137-141.

ZHANG Yingzi, QIU Hua, WU Min. Effect of yarn trajectory crosssectional shape of cyclone on yarn properties[J]. Journal of Textile Research,2011,32(7): 137-141.

[5] WANG X, MIAO M. Reducing yarn hairiness with an air-jet attachment during winding[J]. Textile Research Journal, 1997, 67(7): 481-485.

[6] 方坤.計(jì)算流體力學(xué)的幾種常用軟件[J].煤炭技術(shù)，2006, 25(12): 124-125.

FANG Kun. Optimization design of sound soot blower system in the boiler of electric power plant[J]. Coal Technology,2006, 25(12): 124-125.

[7] 李明, 李明高.STAR-CCM+與流場計(jì)算[M]. 北京: 機(jī)械工業(yè)出版社,2011: 18-22.

LI Ming, LI Minggao. STAR-CCM+ and Airflow Simu-lation[M]. Beijing: Machinery Industry Press,2011: 18-22.

[9] 牟俊玲,邱華,葛明橋.旋流器的氣道流場對成紗質(zhì)量的影響[J]. 山東紡織科技, 2009(1): 36-39.

MOU Junling, QIU Hua, GE Mingqiao. Effect of airway flow field of cyclone on yarn properties[J]. Shandong Textile Technology, 2009(1): 36-39.

[10] 高曉平.降低筒紗毛羽的理論與應(yīng)用研究[D].天津: 天津工業(yè)大學(xué),2004：2.

GAO Xiaoping.The theory and application research of reducing hairiness[D].Tianjing: Tianjin University of Technology, 2004：2.

[11] RENGASAMY R S, PATANAIK A, ANANDJIWALA R D. Simulation of airflow in nozzle-ring spinning using computational fluid dynamics [J]. Textile Research Journal, 2008, 78(5): 412-420.