黎 西，袁銳波*，錢俊兵，汪院林，袁安華，何邦華

1. 昆明理工大學(xué)機(jī)電工程學(xué)院，昆明市呈貢區(qū)景明南路727 號(hào) 650504

2. 云南中煙工業(yè)有限責(zé)任公司技術(shù)中心，昆明市五華區(qū)紅錦路367 號(hào) 650231

煙草加料是卷煙制造的重要工序之一，其中加料噴嘴霧化效果直接影響片煙和香精香料的混合均勻性[1-2]，而噴嘴霧化效果由霧化液滴粒徑分布進(jìn)行評(píng)價(jià)。目前，國(guó)內(nèi)外通常通過(guò)檢測(cè)片煙和香精香料混合均勻性來(lái)經(jīng)驗(yàn)式判斷霧化效果，而非通過(guò)系統(tǒng)性研究加料過(guò)程中料液霧化液滴粒徑分布情況進(jìn)行判斷[3-4]，使得加料工序的調(diào)控存在一定的局限性[5]。由于國(guó)內(nèi)外對(duì)噴嘴霧化過(guò)程的檢測(cè)與評(píng)價(jià)主要是采用圖像處理軟件和激光檢測(cè)噴霧測(cè)量系統(tǒng)相結(jié)合的方式進(jìn)行[6-7]，因此，結(jié)合噴嘴霧化液滴破碎過(guò)程和特點(diǎn)，考慮外混式空氣霧化噴嘴裝置內(nèi)液體流速快、易于調(diào)控氣壓和液體流量、適用高黏度液體分析等優(yōu)點(diǎn)[8-9]，采用高速相機(jī)背光源補(bǔ)償拍照，利用Halcon 軟件確定及選取特征霧化液滴，研究了霧化氣壓力、料液流量、料液溫度對(duì)霧化液滴粒徑分布的影響，評(píng)價(jià)了煙草加料效果，旨在為評(píng)價(jià)卷煙感官質(zhì)量影響因素提供科學(xué)理論驗(yàn)證和相關(guān)技術(shù)支持。

1 材料與方法

1.1 材料與儀器

甲醇，1,2-丙二醇（AR，昆明寶城化工有限公司）；片煙樣品在加料滾筒出口處連續(xù)隨機(jī)抽取，每次取樣100 g 左右，共抽取30 次[樣品取自紅云紅河煙草（集團(tuán)）有限責(zé)任公司制絲車間加料工序]；料液由紅云紅河煙草（集團(tuán)）有限責(zé)任公司提供。

FASTCAM Mini UX100 高速相機(jī)（日本Photron公司）；外混式空氣霧化噴嘴[紅云紅河煙草（集團(tuán)）有限責(zé)任公司提供]；KLB-CMFI-DN3 加料流量計(jì)（北京科力博奧儀表有限公司）；SWPC40 料液溫度傳感器（香港昌輝儀表有限公司）；CC.HBD/2.15/5.0 料液、料溫、霧化氣壓力供給系統(tǒng)（意大利COMAS-SRL 灌裝設(shè)備公司）。

1.2 試驗(yàn)裝置及測(cè)試原理

試驗(yàn)裝置示意圖見(jiàn)圖1。裝置測(cè)試原理為：①利用料液、料溫、霧化氣壓力供給系統(tǒng)進(jìn)行溫度、氣壓和流量的調(diào)節(jié)（如圖2 所示）；②利用高壓空氣的噴散作用，通過(guò)外混式空氣噴嘴以較高的速度將夾帶著料液的高壓空氣噴向加料空間；③利用兩臺(tái)分別安裝在試驗(yàn)臺(tái)軸向中部位置處和徑向中部位置處的高速相機(jī)（如圖3 所示）在背光光源的補(bǔ)償下捕捉高速運(yùn)動(dòng)以及快速破裂的料液霧化液滴，并自動(dòng)將圖片上傳至計(jì)算機(jī)；④采用德國(guó)MVtec 公司Halcon 圖像處理軟件處理和提取高速相機(jī)所拍攝的液滴圖，采用陰影法對(duì)圖像進(jìn)行處理，取得各個(gè)液滴的實(shí)測(cè)值[10]。

1.3 方法

1.3.1 霧化粒徑分布的檢測(cè)

圖1 霧化試驗(yàn)裝置示意圖Fig.1 Schematic diagram of atomization test device

圖2 料液、料溫、霧化氣壓力供給系統(tǒng)Fig.2 Supply systems for casing, material temperature, air pressure

圖3 高速相機(jī)徑向與軸向位置Fig.3 Radial and axial positions of high-speed camera

試驗(yàn)中采用空壓機(jī)為外混式空氣霧化噴嘴提供氣壓，霧化液體為片煙加料料液。試驗(yàn)臺(tái)為長(zhǎng)度2 400 mm，高度和寬度均為1 500 mm 的透明長(zhǎng)方體。由噴嘴噴射料液時(shí)，通過(guò)高速成像技術(shù)，利用分別設(shè)在試驗(yàn)臺(tái)軸向和徑向兩個(gè)位置的高速相機(jī)（圖1）并在背光光源的補(bǔ)償下，通過(guò)控制拍照間隔時(shí)間分別得到不同軸向和徑向位置處液滴霧化圖，最后自動(dòng)傳輸?shù)絇C 上保存及處理，獲得霧化液滴空間坐標(biāo)，進(jìn)而得到霧化液滴粒徑分布情況。同時(shí)采用Halcon 軟件對(duì)圖像進(jìn)行處理，選取了各所需位置特征霧化液滴情況[11-12]，如圖4 所示。在實(shí)際的噴霧過(guò)程中，液體霧化后會(huì)形成大小和形狀不同的液滴組，不同液滴組間尺寸可能相差幾十倍，并且噴霧場(chǎng)中的液滴分布也有較大差別，因此考慮雷諾數(shù)（Re）的影響，根據(jù)量綱分析法，得到SMD 與Re關(guān)系式，從而計(jì)算得到實(shí)驗(yàn)所得的粒徑值。

1.3.2 霧化粒徑區(qū)域劃分及均勻性統(tǒng)計(jì)量數(shù)學(xué)建模將透明試驗(yàn)臺(tái)軸向上分為10 個(gè)研究區(qū)域（見(jiàn)圖5），同時(shí)選取該試驗(yàn)臺(tái)徑向方向正方形截面的內(nèi)接圓作為徑向研究區(qū)域，每個(gè)軸向研究區(qū)域均對(duì)應(yīng)有1 個(gè)徑向研究區(qū)域。利用Halcon 軟件處理高速相機(jī)拍攝得到的霧化液滴像素圖，通過(guò)選取每個(gè)研究區(qū)域軸向和徑向上分布最均勻的霧化液滴區(qū)域作為特征分析區(qū)域，并將這兩個(gè)特征分析區(qū)域作圓約束處理，同時(shí)將該圓約束圖劃分為16等份，得到了基于不同顏色下不同粒徑的霧化液滴效果圖，如圖6 所示。據(jù)此，建立了用來(lái)描述霧化液滴在軸向、徑向區(qū)域分布情況的特征分析區(qū)域的分布函數(shù)?；陔x散型均勻分布法，構(gòu)造用于描述霧化液滴軸向和徑向的空間分布的統(tǒng)計(jì)量tm，見(jiàn)公式（1）。

圖4 Halcon 軟件特征液滴選取Fig.4 Characteristic droplet selection by Halcon software

圖5 軸向分析區(qū)域劃分示意圖Fig.5 Schematic diagram of axial division of testing area

圖6 區(qū)域粒徑分布及特征霧化液滴選取Fig.6 Particle size distribution and characteristic atomized droplet selection

式中：qc為每個(gè)等份內(nèi)的特征霧化液滴數(shù)量的期望值；Bci為每個(gè)等份內(nèi)期望含有的特征霧化液滴面積比；Bc為期望值；T1為液滴個(gè)數(shù)統(tǒng)計(jì)量；T2為液滴面積統(tǒng)計(jì)量；λ為權(quán)重；c 為軸向和徑向系數(shù)，當(dāng)c=α?xí)r為軸向，該軸向處統(tǒng)計(jì)量為tm=t1；當(dāng)c=β時(shí)為徑向，該徑向處統(tǒng)計(jì)量為tm= t2。

經(jīng)查閱相關(guān)論文資料[13]，判定公式（1）中λ=0.2。由公式（1）可得，t1和t2值越小，表示圖像中特征液滴分布的相對(duì)均勻性越好?；谝陨陷S向與徑向t1和t2，構(gòu)造霧化均勻性統(tǒng)計(jì)量T，見(jiàn)公式（2）。

經(jīng)過(guò)試驗(yàn)研究分析，判定公式（2）中的μ=0.3。

1.3.3 霧化粒徑仿真模擬

霧化粒徑由索太爾平均直徑（SMD）表征，在相同的工況下，SMD 值越小表示霧化均勻性越好。選取上述劃分的10 個(gè)區(qū)域中的其中1 個(gè)區(qū)域，通過(guò)分析計(jì)算高速相機(jī)拍攝時(shí)的時(shí)間參數(shù)，獲得該區(qū)域中的噴嘴霧化液滴分布情況。利用ANSYS-FLUENT 仿真軟件分析噴嘴液滴噴灑和霧化的過(guò)程，邊界條件設(shè)置為：入口壓力0.20 MPa，出口壓力0 MPa，料液溫度45 ℃，料液流量20 kg/h，霧化介質(zhì)壓力0.15 MPa。通過(guò)仿真研究，從區(qū)域4 開(kāi)始，霧化液滴總數(shù)開(kāi)始大量增加，因此將區(qū)域4 作為研究區(qū)域。仿真模擬料液霧化情況如圖7 所示。

1.3.4 霧化粒徑計(jì)算

索太爾平均直徑（SMD）是按假設(shè)的液滴群的體積表面積之比與真實(shí)的液滴群的體積表面積之比與真實(shí)的液滴群的體積表面積之比相同來(lái)進(jìn)行推算，根據(jù)量綱分析法，可以得到SMD 與雷諾數(shù)（Re）的關(guān)系[14]，見(jiàn)公式（3）。

圖7 料液流量20 kg/h、料液溫度45 ℃、霧化氣壓力0.2 MPa 時(shí)料液霧化仿真效果Fig.7 Simulation of casing atomization at casing flow rate of 20 kg/h,casing temperature of 45 ℃and air pressure of 0.2 MPa

式中：xi表示所劃分的10 個(gè)軸像區(qū)域（i=1,2,3，…，10）相關(guān)系數(shù)，工況不同，系數(shù)也不同。

式中：ρi表示該實(shí)驗(yàn)所用料液密度為1 037 kg/m3；Δp 表示噴嘴進(jìn)出口壓力差，MPa；ρa(bǔ)為氣體密度，kg/m3；Q1為料液出口流量，kg/h；T1為料液溫度，℃。

由公式（3）和（4）可以獲得區(qū)域4 于高速相機(jī)背光源捕捉拍照和Halcon 軟件特征液滴選取的試驗(yàn)噴嘴霧化液滴粒徑情況。

2 結(jié)果與分析

2.1 霧化液滴粒徑仿真模擬和試驗(yàn)計(jì)算結(jié)果

根據(jù)該煙草公司提供的實(shí)際工藝參數(shù)，分別設(shè)定3 個(gè)研究變量的參數(shù)為[15-17]：霧化氣壓力0.20、0.25、0.35、0.40 MPa；料液流量20、30、40、50 kg/h；料液溫度35、45、55、65 ℃。采用1.3.3 小節(jié)的方法獲得不同工況下區(qū)域4 霧化液滴粒徑仿真理論值，如表1 所示。采用1.3.4 小節(jié)的方法獲得不同工況下區(qū)域4 霧化液滴粒徑試驗(yàn)計(jì)算值，結(jié)果見(jiàn)表2。

2.2 不同工況對(duì)SMD 值的影響

由表1 仿真模擬得到的SMD 值和表2 試驗(yàn)計(jì)算得到的SMD 值，可分析在不同工況影響下區(qū)域4 霧化液滴粒徑參數(shù)的變化情況，如圖8 所示。

由圖8 可以看出：①在料液流量20 kg/h 時(shí)，試驗(yàn)中隨著料液溫度的增加，不同霧化氣壓力下SMD 值呈下降趨勢(shì)，但當(dāng)霧化氣壓力為0.30 MPa時(shí)，SMD 值直線增大；仿真中相同霧化氣壓力下SMD 值隨著料液溫度的增加整體呈減小趨勢(shì)，其中以霧化氣壓力0.20 MPa 時(shí)下降趨勢(shì)更加明顯。②在料液流量30 kg/h 時(shí)，試驗(yàn)和仿真模擬中不同霧化氣壓力下SMD 值整體均為下降趨勢(shì)；這是因?yàn)殡S著料液溫度的提高，料液黏度逐漸減小，破碎需要的動(dòng)能有所減小，所以粒徑逐漸減小[18]。③在料液流量40 kg/h 時(shí)，試驗(yàn)中不同霧化氣壓力下SMD 值趨于減小趨勢(shì)，但霧化氣壓力為0.30 MPa 時(shí)，SMD 值呈增大減小再增大的不規(guī)則變化；仿真中SDM 值均呈下降趨勢(shì)。④在料液流量50 kg/h 時(shí)，試驗(yàn)中不同霧化氣壓力下SMD 值隨著料液溫度的升高呈不規(guī)則的增加趨勢(shì)，這是因?yàn)殡S著料液流量的增加，液滴破碎所需能量增加，但霧化氣壓力提供的能量有限，所以料液流量增加到一定程度后霧化粒徑逐漸增大。其中以霧化氣壓力0.25、0.30 MPa 時(shí)變化最為明顯；在仿真中，SDM值均呈下降趨勢(shì)。

表1 不同工況下區(qū)域4 霧化液滴粒徑理論值Tab.1 Theoretical value of atomized droplet size in area No.4 under different working conditions

表2 不同工況下區(qū)域4 霧化液滴粒徑試驗(yàn)計(jì)算值Tab.2 Calculated value of atomized droplet size in area No.4 under different working conditions

由于SMD 值越小，霧化液滴粒徑越小，分布越均勻，卷煙生產(chǎn)后續(xù)工序（煙葉加料混合）物料的均勻性越好?；谏鲜鰧?shí)驗(yàn)結(jié)果，可得出工況條件為料液溫度35 ℃、料液流量30 kg/h、霧化氣壓力0.30 MPa 時(shí)，得到的SMD 值最好。

圖8 流體參數(shù)及溫度對(duì)區(qū)域4 液滴SMD 的影響Fig.8 Influences of fluid parameters and temperature on SMD in area No.4

圖8 中試驗(yàn)結(jié)果與仿真結(jié)果存在一定的誤差，這是因?yàn)椋孩俑咚傧鄼C(jī)所拍攝得到的霧化液滴粒徑均較大，采用陰影法求得的霧化粒徑也忽略了細(xì)小的霧化液滴。②試驗(yàn)采用手動(dòng)選取測(cè)量，然后進(jìn)行數(shù)學(xué)計(jì)算；而仿真則是直接得到結(jié)果。

2.3 不同區(qū)域內(nèi)霧化液滴粒徑分布

取最佳工況進(jìn)行舉例分析，通過(guò)高速相機(jī)拍攝和相應(yīng)軟件預(yù)處理，將特征霧化液滴的大小和位置數(shù)據(jù)導(dǎo)入Matlab 中編程計(jì)算，由公式（1）、（2）可以得到徑向處10 個(gè)區(qū)域不同工況下的均勻性，如表3 所示。

根據(jù)均勻分布理論，均勻性統(tǒng)計(jì)量越低表明霧化效果越好，霧化液滴體積越小。從表3 可以看出，所劃分的10 個(gè)區(qū)域霧化液滴粒徑均勻性統(tǒng)計(jì)量均逐漸降低，說(shuō)明霧化液滴粒徑分布也逐漸均勻。其中區(qū)域1 到區(qū)域3 霧化粒徑分布呈不規(guī)則變化，區(qū)域2 T 值最大，為0.206 79；但從區(qū)域4開(kāi)始，霧化效果逐漸增加，T 值由0.190 97 逐漸降到0.117 51。

利用高速相機(jī)拍攝得到的10 個(gè)區(qū)域徑向的霧化液滴像素圖，通過(guò)配套電腦系統(tǒng)性的處理，導(dǎo)出各徑向區(qū)域像素圖中霧化液滴的二維坐標(biāo)點(diǎn)，將其導(dǎo)入Matlab 中，處理得到10 個(gè)研究區(qū)域徑向霧化液滴分布圖，見(jiàn)圖9。由圖9 可以看出區(qū)域1、2、3 中霧化液滴處于聚合狀態(tài)，霧化液滴較多，說(shuō)明此階段霧化液滴并未開(kāi)始破碎；區(qū)域4 中霧化液滴開(kāi)始均勻分布整個(gè)空間；區(qū)域5、6、7 中霧化液滴數(shù)量開(kāi)始逐漸減少，說(shuō)明這3 個(gè)區(qū)域中液滴開(kāi)始破碎并分散；區(qū)域8、9、10 中拍攝得到的霧化液滴數(shù)量逐漸減少，說(shuō)明該階段液滴已經(jīng)大量破碎并分布于整個(gè)空間。圖9 比較完整地描述了霧化液滴破碎過(guò)程。

表3 不同區(qū)域霧化液滴徑向均勻統(tǒng)計(jì)量的變化Tab.3 Variation of radial uniform statistics of atomized droplets in different areas

圖9 不同徑向區(qū)域粒徑分布試驗(yàn)圖Fig.9 Particle size distribution in different radial areas

2.4 特征區(qū)域霧化粒徑均勻性分析

由于從區(qū)域4 開(kāi)始，霧化液滴開(kāi)始均勻分布整個(gè)空間，液滴大量破碎和均勻擴(kuò)散，該區(qū)域的霧化情況影響著后續(xù)整個(gè)實(shí)驗(yàn)過(guò)程中液滴的霧化效果。因此選取區(qū)域4 作為特征區(qū)域，通過(guò)研究該區(qū)域中不同工況下霧化液滴粒徑的均勻性，可以大體表征整個(gè)霧化過(guò)程中液滴粒徑分布均勻性。

圖10 流體參數(shù)及溫度對(duì)區(qū)域4 液滴霧化均勻性統(tǒng)計(jì)量的影響Fig.10 Influences of fluid parameters and temperature on atomization uniformity statistics in area No.4

區(qū)域4 中，流體參數(shù)及溫度對(duì)霧化均勻性統(tǒng)計(jì)量的影響如圖10 所示。由圖10 可以看出，隨著霧化氣壓力逐漸增加，軸向霧化液滴均勻性統(tǒng)計(jì)量t1直線下降，徑向霧化液滴均勻性統(tǒng)計(jì)量t2呈先增大后降低趨勢(shì)，液滴破碎狀態(tài)逐漸趨于小粒徑顆粒變化狀態(tài)，霧化效果T 逐漸變佳；其他條件不變的情況下，霧化氣壓力0.35 MPa 為最佳指標(biāo)。隨著料液流量逐漸增加，t1先減小后增大，而t2先增大后減小，總體霧化液滴均勻性統(tǒng)計(jì)量T 逐漸增加，霧化效果逐漸變差；其他條件不變情況下，料液流量50 kg/h 為最差工況。隨著料液溫度的增加，t1和t2及霧化液滴均勻性統(tǒng)計(jì)量先減小后增加，霧化效果先變佳后變差；最佳料液溫度為55 ℃。

3 結(jié)論

①流體參數(shù)固定不變的情況下，隨著料液溫度、霧化氣壓力的增加，SMD 值呈不規(guī)則變化，其中以料液溫度55 ℃、霧化氣壓力0.20 MPa 時(shí)SMD值最??；以料液溫度35 ℃、霧化氣壓力0.35 MPa時(shí)SMD 值最大。②隨著料液流量、霧化氣壓力的增加，SMD 值呈先減小后增大趨勢(shì)，其中以料液流量30 kg/h、霧化氣壓力0.3 MPa 時(shí)SMD 值最小。③隨著料液流量、料液溫度的增加，SMD 值逐漸減小，其中以料液溫度65 ℃、料液流量30 kg/h 時(shí)SMD 值最大；整體上實(shí)驗(yàn)和仿真的數(shù)據(jù)變化趨勢(shì)較為吻合，但由于實(shí)驗(yàn)器材的質(zhì)量，陰影法測(cè)量的不確定性以及仿真結(jié)果的差異性，產(chǎn)生了相對(duì)誤差。④利用離散型均勻分布法分析了不同區(qū)域內(nèi)噴嘴霧化效果，其中霧化液滴分布均勻性統(tǒng)計(jì)量以區(qū)域4 的差異最大，故選作進(jìn)行特征區(qū)域分析。在特征區(qū)域內(nèi)，噴嘴霧化液滴開(kāi)始大量破碎、擴(kuò)散；在其他變量一定下，隨著霧化氣壓力的增加T 值逐漸降低；隨著料液流增加T 值逐漸增加；隨著料液溫度增加T 值先減小后增加。T 值越小，霧化效果越好，反之越差。軸向和徑向均勻性統(tǒng)計(jì)量相對(duì)于霧化均勻性指標(biāo)誤差范圍在0.17%～0.46%之間。