豬舍靜電除塵器細(xì)顆粒脫除效應(yīng)數(shù)值模擬及試驗(yàn)*

胡捷，文建萍，黎波，楊衛(wèi)平，肖志鋒

(1. 江西農(nóng)業(yè)大學(xué)工學(xué)院，南昌市，330045； 2. 江西省畜牧設(shè)施技術(shù)開發(fā)工程研究中心，南昌市，330045)

0 引言

豬場(chǎng)中的顆粒物可分為細(xì)顆粒物(PM2.5，空氣動(dòng)力學(xué)直徑分別≤2.5 μm)和可吸入顆粒物[1]，其中PM2.5是細(xì)菌、病毒傳播的載體，產(chǎn)生巨大的危害[2-3]。此外，細(xì)菌、病毒、NH3會(huì)吸附在豬舍PM2.5表面發(fā)生物理和化學(xué)反應(yīng)形成危害性更強(qiáng)的二次污染[4-7]。因此，豬舍顆粒污染物脫除關(guān)系到舍內(nèi)外人畜健康及周邊環(huán)境空氣質(zhì)量。

豬舍EP可對(duì)顆粒物進(jìn)行有效脫除[8-9]。吳新構(gòu)建了一套EP系統(tǒng)，在保育豬舍內(nèi)進(jìn)行試驗(yàn)，結(jié)果表明該系統(tǒng)能夠去除70%的粉塵；王樹華等[10]在保育豬舍內(nèi)使用3DDF-450型EP系統(tǒng)，與對(duì)照豬舍相比后發(fā)現(xiàn)該系統(tǒng)最大可降低72.59%的舍內(nèi)TSP濃度(P<0.01)。然而，上述試驗(yàn)結(jié)果同時(shí)表明EP對(duì)細(xì)顆粒PM2.5的脫除效率比大顆粒要低。EP對(duì)豬舍顆粒污染物的捕集能力與EP內(nèi)電場(chǎng)和流場(chǎng)特性密切相關(guān)[11]。但是，強(qiáng)電場(chǎng)的存在給EP內(nèi)含塵氣體的流動(dòng)和顆粒輸運(yùn)的直接測(cè)量帶來極大的困難[11]。在此情況下，數(shù)值模擬成為揭示豬舍EP內(nèi)輸運(yùn)行為，實(shí)現(xiàn)豬舍EP對(duì)細(xì)顆粒脫除工藝優(yōu)化的有效工具[12-13]。

本文針對(duì)豬舍細(xì)顆粒PM2.5脫除效率低的問題，重點(diǎn)開展顆粒在EP內(nèi)的輸運(yùn)行為研究。將采用多場(chǎng)耦合數(shù)值方法，分析豬舍線板型EP內(nèi)電場(chǎng)、流場(chǎng)和細(xì)顆粒運(yùn)動(dòng)的特性，通過可視化試驗(yàn)驗(yàn)證電暈線電壓、流場(chǎng)流速對(duì)細(xì)顆粒物捕集效率的控制規(guī)律，提高豬舍空氣質(zhì)量。

1 設(shè)計(jì)與方法

1.1 工作原理

靜電除塵器(Electrostatic precipitator, EP)主要由電暈線(Corona wire, CW)、捕塵極板(Collecting plate, CP)、氣溶膠通道、靜電控制組件等組成，該器件又稱為靜電空間電荷系統(tǒng)(Electrostatic Space Charge System， ESCS)。ESCS系統(tǒng)一般安裝在天花板等不易接觸的位置，對(duì)人畜無干擾、維護(hù)要求低，如圖1所示。在當(dāng)下集約化、封閉化的豬舍中，大氣中的含塵空氣通過進(jìn)風(fēng)系統(tǒng)進(jìn)入豬舍中，含塵氣流混合豬舍內(nèi)產(chǎn)生的顆粒污染物通過舍內(nèi)空氣流動(dòng)進(jìn)入ESCS系統(tǒng)入口，由電暈效應(yīng)使得顆粒帶電，顆粒物或粉塵在電場(chǎng)力和流體曳力的共同作用下向捕塵極板聚集，實(shí)現(xiàn)顆粒污染物的吸附。含塵氣流通過ESCS系統(tǒng)后不僅輸出干凈、綠色的空氣至豬舍內(nèi)部，同時(shí)通過排風(fēng)系統(tǒng)輸出至豬舍外部，提高了豬舍內(nèi)外的空氣質(zhì)量，形成綠色、低排放的豬養(yǎng)殖場(chǎng)空氣循環(huán)。

圖1 豬舍靜電除塵器工作原理示意圖

1.2 可視化豬舍顆粒脫除試驗(yàn)

顆粒污染物輸運(yùn)可視化試驗(yàn)平臺(tái)如圖2所示，該裝置包括：氣溶膠顆粒發(fā)生系統(tǒng)、風(fēng)機(jī)、氣溶膠流道、高壓電源、高速相機(jī)、顆粒物荷電量裝置、DP-02激光粒度儀、計(jì)算機(jī)、流速傳感器。

圖2 顆粒污染物輸運(yùn)可視化試驗(yàn)平臺(tái)

采用蚊香顆粒來研究PM2.5粒徑范圍及物性特征顆粒物的運(yùn)動(dòng)及捕集特性。通過氣溶膠顆粒發(fā)生系統(tǒng)產(chǎn)生蚊香顆粒的方法為：蚊香置于容積為50 L的不銹鋼儲(chǔ)煙罐中，通過壓縮空氣瓶提供燃燒所需要的標(biāo)準(zhǔn)空氣，空氣流量由ALICAT 21系列標(biāo)準(zhǔn)型數(shù)字質(zhì)量流量控制器來控制。試驗(yàn)中，需要根據(jù)可視化系統(tǒng)的拍攝要求選擇合適的氣流速度，通過調(diào)節(jié)送料機(jī)入口處的給氣壓力調(diào)節(jié)旋鈕，可以實(shí)現(xiàn)對(duì)氣流速的控制。

蚊香氣溶膠流道采用特氟龍材質(zhì)，內(nèi)部尺寸為50 mm×40 mm×1 000 mm，在距入口200 mm處和800 mm處前后開有石英玻璃觀察窗口，距入口100 mm 處水平插入直徑1 mm的不銹鋼針電極，電極上下流道內(nèi)壁面嵌有接地鋼極板，同時(shí)，距出口400 mm 處水平等距插入4根直徑1 mm的不銹鋼針電極，其間距為100 mm，電極上下流道內(nèi)壁面嵌有接地鋼極板，其作為電暈放電區(qū)。石英玻璃窗口處為可視化觀測(cè)界面。

高壓電源的型號(hào)為TRC2020P50-150，其輸出參數(shù)為：電壓0～±100 kV連續(xù)可調(diào)，電流0～3 mA連續(xù)可調(diào)，分別通過正、負(fù)直流電源實(shí)現(xiàn)對(duì)顆粒物的荷電。利用DP-02型激光粒度儀可實(shí)現(xiàn)對(duì)顆粒物氣溶膠粒徑分布變化和濃度的實(shí)時(shí)在線檢測(cè)，激光粒度儀與電腦連接，通過計(jì)算機(jī)中分析軟件可以得到輸出顆粒的濃度分布，進(jìn)而獲得捕塵極板上的顆粒濃度，最終獲得顆粒捕集效率。該試驗(yàn)驗(yàn)證裝置的除塵效率可由式(27)表征。

1.3 數(shù)學(xué)模型

1.3.1 電暈?zāi)Ｐ?/p>

1) 基本假設(shè)。(1)忽略電暈放電產(chǎn)生的磁場(chǎng)。(2)在靜電除塵器內(nèi)，電荷移動(dòng)速度約100 m/s，空氣的流速約為1m/s[14]，因此空氣流動(dòng)對(duì)電場(chǎng)分布的影響可以忽略。(3)電遷移率恒定。

2) 控制方程。電暈?zāi)Ｐ筒捎秒姾墒睾惴匠毯筒此煞匠糖蠼廨d流子的傳輸問題。載流子的傳輸包括電場(chǎng)漂移和對(duì)流。沒有源項(xiàng)，域方程

·J=0

(1)

J=zqγρqE+ρu

(2)

ε02V=-ρq

(3)

式中：J——電流密度；

zq——電荷數(shù)；

γ——電遷移率；

ρq——空間電荷數(shù)密度；

E——電場(chǎng)；

u——流體速度；

V——電勢(shì)；

ε0——真空介電常數(shù)。

可以修改這組方程以獲得控制方程

(4)

值得注意的是，域方程不包含任何與等離子體產(chǎn)生和維持有關(guān)的信息。所有等離子體物理場(chǎng)都縮減在內(nèi)電極的邊界條件中。

3) 邊界條件。電暈電極上電場(chǎng)的法向分量用作泊松方程的邊界條件

n·E=E0

(5)

泊松方程的其他邊界條件是捕塵極板上的V=0以及入口和出口處的零電荷。方程(4)的邊界條件涉及采用拉格朗日乘子求得電暈電極上的空間電荷密度ρq，從而驗(yàn)證施加的電勢(shì)V0。

V-V0=0

(6)

在模型中，電勢(shì)和電場(chǎng)都施加在電暈電極上。為了獲得預(yù)測(cè)結(jié)果，電暈處電場(chǎng)值必須足夠接近真實(shí)值。這里用的是匹克定律

(7)

式中：E0——擊穿電場(chǎng)；

δ——133.322 Pa、293.15 K下歸一化到氣體密度的氣體數(shù)密度；

ri——電暈電極的直徑。

1.3.2 湍流模型

1) 基本假設(shè)。(1)靜電除塵器內(nèi)空氣是不可壓縮牛頓流體。(2)含塵氣體經(jīng)過濕簾、過濾等一次除塵進(jìn)入豬舍內(nèi)。進(jìn)入靜電除塵器的氣體僅含有少量的大顆粒粉塵和小顆粒粉塵，粉塵濃度很低，因此顆粒的存在對(duì)流體運(yùn)動(dòng)的影響可以忽略[15-16]。

(8)

u′——速度脈動(dòng)；

(9)

式中：μT——湍流動(dòng)力黏度。

因此只需求出μT即可。通過引入兩個(gè)額外的傳遞方程和兩個(gè)變量(湍動(dòng)能k和湍動(dòng)能耗散率ε)可求出μT。k-ε模型

(10)

(11)

(12)

式中：k——湍動(dòng)能；

ε——湍動(dòng)能耗散率；

σk、σε——k、ε對(duì)應(yīng)的普朗克常數(shù)；

Cμ、Cε1、Cε2——經(jīng)驗(yàn)常數(shù)，其值如表1所示。

k-ε模型由式(10)～式(12)構(gòu)成，其中式(10)～式(12)分別為湍流黏度方程，湍動(dòng)能k輸運(yùn)方程及生成項(xiàng)，湍動(dòng)能耗散率ε輸運(yùn)方程。本文利用式(8)～式(12) 對(duì)豬舍內(nèi)流場(chǎng)進(jìn)行建模求解，以獲得舍內(nèi)的速度、壓力和渦量分布。

表1 經(jīng)驗(yàn)常數(shù)Tab. 1 Empirical constant

3) 邊界條件。(1)在壁面處采用無滑移邊界條件。(2)在出口處，應(yīng)用壓力出口條件。

1.3.3 顆粒追蹤模型

1) 基本假設(shè)。(1)忽略顆粒之間的相互作用。(2)顆粒在到達(dá)捕塵極板之后會(huì)被完全吸附，不會(huì)有顆粒反彈。

2) 控制方程。遵循牛頓第二定律，通過求解顆粒位置矢量分量的二階運(yùn)動(dòng)方程來計(jì)算顆粒位置，表達(dá)式

(13)

式中：q——顆粒位置；

v——顆粒速度；

mp——顆粒質(zhì)量；

Ft——施加在顆粒上的總力。

模型中，作用在顆粒上的力是曳力和電力。由于顆粒半徑變得非常小，因此曳力中需要包含稀薄效應(yīng)。這里，曳力FD用Cunningham-Millikan-Davis模型

(14)

式中：τp——顆粒速度響應(yīng)時(shí)間。

(15)

式中：ρp——顆粒的密度；

dp——粒徑；

CD——曳力系數(shù)；

Rer——相對(duì)雷諾數(shù)。

(16)

S是曳力校正系數(shù)，定義

(17)

式中：Kn——無量綱系數(shù)。

作用在顆粒上的電力Fe(SI單位：N)定義

Fe=eZE

(18)

式中：e——是元電荷；

Z——每個(gè)顆粒上的累積電荷數(shù)。

使用Lawless模型計(jì)算顆粒上累積的電荷，表達(dá)式

(19)

式中：τc——特性充電時(shí)間。

(20)

式中：kB——玻爾茲曼常數(shù)；

Ti——離子溫度。

Rf和Rd分別是由場(chǎng)和擴(kuò)散傳輸引起的無量綱充電速率，定義

(21)

(22)

式(21)和式(22)中，ve和vs分別表達(dá)

(23)

(24)

(25)

式中：εr,p——顆粒的相對(duì)介電常數(shù)；

fa——將擴(kuò)散充電速率和場(chǎng)充電速率相結(jié)合的函數(shù)。

(26)

捕集效率是衡量靜電除塵器性能的重要指標(biāo)。

(27)

式中：η——捕集效率，無量綱；

N1——靜電除塵器入口處顆粒數(shù)量；

N2——靜電除塵器出口處顆粒數(shù)量。

3) 邊界條件。電暈線：采用對(duì)稱邊界條件。

1.3.4 計(jì)算條件、計(jì)算方法和收斂標(biāo)準(zhǔn)

參照通用型3DDF-450、3DDF-300等靜電除塵器的幾何尺寸，模擬需要的器件幾何參數(shù)和物性參數(shù)分別如表2、表3所示。

表2 靜電除塵器幾何尺寸范圍Tab. 2 Dimension scope of electrostatic precipitator

表3 靜電除塵器物性參數(shù)Tab. 3 Physical parameters of EP

采用有限元法(Finite Element Method, FEM)對(duì)二維模型進(jìn)行網(wǎng)格劃分，局部網(wǎng)格如圖3所示。本文采用自由剖分三角形網(wǎng)格，網(wǎng)格大小為標(biāo)準(zhǔn)，最大單元尺寸3.5 mm，最小單元尺寸0.1 mm，最大單元生長率為1.13，曲率因子0.3，狹窄區(qū)域解析度1。對(duì)于捕塵極板附近及電極線附近設(shè)置邊界層并進(jìn)行網(wǎng)格加密，邊界層數(shù)為2，邊界層拉伸因子為1.2，厚度調(diào)節(jié)因子為5。整個(gè)計(jì)算域劃分為25 272個(gè)三角形網(wǎng)格。并利用MATLAB求解上述偏微分方程。電場(chǎng)計(jì)算收斂標(biāo)準(zhǔn)為

(28)

(29)

式中：K——迭代次數(shù)。

在流場(chǎng)計(jì)算過程中，采用SIMPLE算法處理壓力和速度的耦合。收斂標(biāo)準(zhǔn)為

(30)

(31)

圖3 局部網(wǎng)格

2 模型驗(yàn)證

2.1 電勢(shì)分布

圖4給出了文獻(xiàn)[17]的電勢(shì)測(cè)量試驗(yàn)點(diǎn)位置A-A。計(jì)算過程中，參照文獻(xiàn)中電勢(shì)探針測(cè)量X=228.6 mm 處電勢(shì)剖面的電勢(shì)分布，X表示以入口為起點(diǎn)的橫軸位置，Y表示X=228.6 mm位置上以捕塵極板為起點(diǎn)到電極的縱軸位置。結(jié)果表明靜電除塵器內(nèi)的電勢(shì)分布是不均勻的，電極線周圍電勢(shì)最高。將X=228.6 mm剖面下測(cè)量的電勢(shì)與試驗(yàn)值進(jìn)行了比較，如圖5所示，表明電勢(shì)的數(shù)值解和文獻(xiàn)[17]的試驗(yàn)數(shù)據(jù)值吻合良好，計(jì)算得到的電勢(shì)分布與試驗(yàn)結(jié)果在誤差范圍內(nèi)，說明數(shù)值模擬建立的電勢(shì)是正確的。

圖4 文獻(xiàn)[17]的靜電除塵器

圖5 電勢(shì)計(jì)算值與試驗(yàn)值對(duì)比

2.2 流場(chǎng)分布

靜電除塵器的流場(chǎng)計(jì)算采用文獻(xiàn)[19]的試驗(yàn)條件，數(shù)值模擬分別采用了RNGk-ε湍流模型。圖6表明，使用該模型得到的流場(chǎng)流體速度非常接近試驗(yàn)值，在56 mm處計(jì)算獲得的流體速度與試驗(yàn)值誤差最大，為12.4%。需要指出的是，10 mm，90 mm分別為靠近收塵極和電暈線處；由圖可知，該兩處的計(jì)算結(jié)果與試驗(yàn)結(jié)果的誤差較小。同時(shí)，文獻(xiàn)指出在捕塵極板處，存在明顯的回流區(qū)，而在電暈線處存在明顯的邊界層分離現(xiàn)象。因此，RNGk-ε湍流模型對(duì)靜電除塵器內(nèi)流場(chǎng)的計(jì)算結(jié)果比較準(zhǔn)確。

圖6 流速計(jì)算值與試驗(yàn)值對(duì)比

2.3 捕集效率計(jì)算

靜電除塵器的收塵率計(jì)算采用文獻(xiàn)[20]的試驗(yàn)條件，采用RNGk-ε模擬獲得的數(shù)值結(jié)果如圖7所示。在不同入口速度下，靜電除塵器收塵效率的預(yù)測(cè)值與試驗(yàn)值基本吻合。當(dāng)入口流速從0.15 m/s上升到0.35 m/s 時(shí)，顆粒捕集效率η有明顯下降，顆粒捕集效率從81.2%下降到64.5%。當(dāng)入口速度為0.25 m/s 時(shí)，模擬值與試驗(yàn)值相對(duì)誤差最大，其值為2.6%，該誤差同時(shí)證明了本文數(shù)值建模與求解的準(zhǔn)確性和合理性。

圖7 不同速度下收塵效率的計(jì)算值與試驗(yàn)值

3 結(jié)果與討論

3.1 靜電除塵器內(nèi)電勢(shì)和電荷密度分布

圖8和圖9展示了在V=30 kV下，二維靜電除塵器XOY面的電場(chǎng)分析特點(diǎn)。(1)電勢(shì)和電荷密度均以電暈線的連線為對(duì)稱軸呈對(duì)稱分布。(2)在相鄰電暈線連線的中點(diǎn)處，電勢(shì)和電荷密度出現(xiàn)極小值。(3)電勢(shì)和電荷密度以電暈線為中心呈環(huán)狀分布；遠(yuǎn)離電暈線時(shí)，電勢(shì)和電荷密度逐漸變?yōu)闄E圓分布。這是因?yàn)樵谶h(yuǎn)離電暈線的過程中，鄰近電暈線對(duì)電勢(shì)和電荷密度影響越來越明顯。(4)可視化闡明了，根據(jù)式(4)，靜電除塵器內(nèi)部的電場(chǎng)由施加在放電極上的電壓V以及空間電荷ρq共同作用形成。

圖8 二維線-板式靜電除塵器XOY面的電勢(shì)分布

圖9 二維線-板式靜電除塵器XOY面的電荷密度分布

3.2 不同直徑顆粒的運(yùn)動(dòng)軌跡和捕集效率

圖10展示了當(dāng)u=1 m/s、V=20 kV時(shí)，顆粒收集效率與顆粒直徑的關(guān)系。從捕塵效率角度具體分析，(1)當(dāng)0.01 μm

圖11顯示了4個(gè)不同直徑顆粒的顆粒軌跡和顆粒中累積的電荷。顆粒在左側(cè)釋放，并在流體流動(dòng)中向右側(cè)出口傳輸。顆粒沿著其軌跡逐漸變?yōu)閹щ姞顟B(tài)，產(chǎn)生的電力使其軌跡在壁方向上發(fā)生偏轉(zhuǎn)。顆粒半徑影響顆粒受到的曳力和電力的平衡，從而影響顆粒軌跡和靜電除塵器的收集效率。綜上可知，在顆粒尺寸大小兩種極端情況下，收集效率更高，即dp大于5 μm 和dp小于0.05 μm時(shí)η均超過了50%。這兩個(gè)極端表現(xiàn)是因?yàn)?，顆粒越大會(huì)獲得越大的電荷，顆粒越小則受到的曳力越小。在這兩個(gè)極端之間，曳力對(duì)顆粒軌跡的影響最大，同時(shí)，電場(chǎng)荷電導(dǎo)致的顆粒電遷移率增大，擴(kuò)散荷電導(dǎo)致的顆粒電遷移率減小，二者的綜合作用在這兩個(gè)極端之間粒徑最弱，因此0.05 μm

圖10 0.01～10 μm粒徑顆粒與捕塵效率之間的關(guān)系

(a) dp=0.01 μm

(b) dp=0.2 μm

(c) dp=2.0 μm

(d) dp=5.0 μm

3.3 流速對(duì)細(xì)顆粒捕集效率的影響

圖12展示了當(dāng)V=20 kV時(shí)，顆粒收集效率與流速大小的關(guān)系。本文探究了PM2.5主要粒徑(dp為0.2 μm 和2.0 μm)在0.2～2.0 m/s范圍內(nèi)的捕集效率變化特征，該流速變化為典型豬舍舍內(nèi)空氣流速變化范圍。由圖可知，當(dāng)流速減小時(shí)，顆粒捕集效率η有明顯上升，均成線性變化。2 μm徑粒顆粒捕集效率從40%上升到100%；0.2 μm徑粒顆粒捕集效率從8%上升到71%。該變化規(guī)律說明了靜電除塵器除塵過程中空氣流速是影響顆粒污染物脫除的一個(gè)重要因素，同時(shí)，基于上述u-η曲線可有效調(diào)控舍內(nèi)顆粒分布。

圖12 顆粒收集效率與流速大小的關(guān)系

3.4 電暈線電壓對(duì)細(xì)顆粒捕集效率的影響

圖13展示了當(dāng)u=1 m/s時(shí)，顆粒收集效率與電暈線電壓的關(guān)系。由圖可知，施加在電暈線上的電壓越高，靜電除塵器對(duì)于dp=0.2 μm和2.0 μm的脫除效果就越高，V-η曲線呈現(xiàn)線性變化，且試驗(yàn)結(jié)果與算計(jì)結(jié)果相一致。因?yàn)殡S著V升高，除塵器中的電場(chǎng)強(qiáng)度也隨之增大，PM2.5顆粒所受到的電場(chǎng)力隨之上升，細(xì)顆粒向捕塵板的偏移效果加強(qiáng)，從而η增大。在較低電壓區(qū)間內(nèi)V的升高對(duì)η的提高比較明顯；在較高電壓區(qū)間內(nèi)V的升高對(duì)除塵率也有提高，但對(duì)于dp=2 μm在電壓達(dá)到45 kV后提高效果不明顯。

圖13 顆粒收集效率與電壓的關(guān)系

試驗(yàn)結(jié)果表明，當(dāng)u=1 m/s,V=50 kV時(shí)，η可達(dá)100%。因此，上述V-η曲線可有效的調(diào)控舍內(nèi)顆粒分布，提高PM2.5的脫除效果。

4 結(jié)論

本文以豬舍靜電除塵器為載體，針對(duì)豬舍細(xì)顆粒PM2.5脫除效率低的問題，利用電場(chǎng)-流場(chǎng)耦合的數(shù)值模擬和試驗(yàn)手段重點(diǎn)分析了顆粒在靜電除塵器內(nèi)的傳輸行為。

1) 靜電除塵器內(nèi)電勢(shì)和電荷密度以電暈線為中心呈環(huán)狀分布。0.01～10 μm粒徑范圍內(nèi)，捕集效率η隨粒徑dp的增大，先減小，再保持不變，再增大，直至恒定不變。

2) 0.2～2.0 m/s流速范圍內(nèi)，當(dāng)流速u減小時(shí)，PM2.5顆粒的η有明顯上升，u-η曲線成線性變化。15～50 kV電暈線電壓范圍內(nèi)，當(dāng)電暈線電壓V增大時(shí)，PM2.5顆粒的η有明顯上升，V-η曲線成線性變化。因此，利用上述數(shù)值模型所獲得的方法及規(guī)律為豬舍內(nèi)細(xì)顆粒污染物的綠色高效脫除技術(shù)奠定了科學(xué)方法，增強(qiáng)了對(duì)豬舍細(xì)顆粒的脫除效果，解決了豬舍內(nèi)PM2.5難脫除的問題。