摘要

基于comsolmultiphysics對(duì)電除塵器的二維模型進(jìn)行了數(shù)值模擬，建立了以靜電場(chǎng)、湍流場(chǎng)和流體粒子為追蹤模塊的耦合物理場(chǎng)，在有ｅｈｄ流的情況下分析了收塵極版形狀（ｐ型板、ｂｅ型板、ｃｗ型板和m型板）、風(fēng)速、電壓、粒徑和溫度對(duì)靜電除塵器內(nèi)靜電場(chǎng)、湍流場(chǎng)和除塵效率的影響。結(jié)果表明：收塵極板形狀對(duì)電勢(shì)分布有影響，ｂｅ板電勢(shì)不均勻程度最高；不同極板的除塵效率與放電極電壓和顆粒物粒徑正相關(guān)，與一次風(fēng)速和溫度負(fù)相關(guān)；在不同工況下，ｂｅ板除塵效率最佳；增強(qiáng)離子風(fēng)有助于提高除塵效率。

0 引言

根據(jù)世界能源統(tǒng)計(jì)年鑒，在2019年我國(guó)一次能源的消耗量達(dá)到了141.7ｅｊ，可以占到全球總消耗量的24.3％，是世界第一。而煤炭是我國(guó)的第一大能源，2019年我國(guó)對(duì)于煤炭的消耗在我國(guó)能源消費(fèi)結(jié)構(gòu)中可以占到57.7％。而火力發(fā)電是我國(guó)主要的發(fā)電方式，據(jù)國(guó)家統(tǒng)計(jì)局?jǐn)?shù)據(jù)，2019年我國(guó)的總發(fā)電量為75034.28億千瓦時(shí)，火力發(fā)電量達(dá)到52201.48億千瓦時(shí)，而火力發(fā)電量占我國(guó)總發(fā)電量的比例為69.6％。靜電除塵器是我國(guó)燃煤電廠廣泛使用的燃煤煙氣處理設(shè)施，占比能夠達(dá)到90％～95％，它具有能耗低、煙氣處理量大、效率高等優(yōu)點(diǎn)。

在近年來(lái)人們對(duì)于靜電除塵器的流場(chǎng)和提高除塵效率的因素進(jìn)行了研究，而在靜電除塵器中會(huì)產(chǎn)生離子風(fēng)，它會(huì)對(duì)電除塵器中流場(chǎng)和顆粒物運(yùn)動(dòng)規(guī)律產(chǎn)生一定的影響。離子風(fēng)（ionicwind）也稱為電風(fēng)或電暈風(fēng)（coronawind），它是指在電暈放電的過(guò)程中由電暈放電產(chǎn)生的離子推動(dòng)離子周圍的氣流而產(chǎn)生的一種現(xiàn)象。對(duì)于離子風(fēng)的研究最早可以追溯到十八世紀(jì)，1709年，hauksbee發(fā)現(xiàn)把帶電的圓管靠近臉就會(huì)感受到輕微的風(fēng)；1922年，deutsch在不考慮離子風(fēng)引起粒子驅(qū)進(jìn)速度變化的情況下，提出了著名的多依奇公式，這個(gè)公式表明：電除塵效率與粒子驅(qū)進(jìn)速度和收塵極板的面積成正比，與氣流量成反比。而在近年來(lái)國(guó)內(nèi)外在離子風(fēng)方面取得了很多成果，shimamoto在2004年發(fā)現(xiàn)電除塵器中離子風(fēng)增強(qiáng)了湍流強(qiáng)度，而且離子風(fēng)在氣流速度低的時(shí)候作用明顯；ｃｈｕｎ等人使用simplest算法研究了線板ｅｓｐ中ｋ－?湍流模型的二維電流體動(dòng)力學(xué)流動(dòng)（ｅｈｄ）模式；ｆｅｎｇ等采用有限元與有限體積混合方法（ｆｅm－ｆｖm）分析了電除塵器中渦流對(duì)ｅｈｄ流的影響，并且驗(yàn)證了幾種經(jīng)典的ｅｈｄ模型。

近年來(lái)國(guó)內(nèi)外很多學(xué)者利用數(shù)值模擬軟件對(duì)電除塵器進(jìn)行了數(shù)值模擬，王曉華對(duì)不同直徑的電暈線進(jìn)行了數(shù)值模擬發(fā)現(xiàn)隨著電暈線直徑的減小收塵板處的電荷密度和電場(chǎng)強(qiáng)度增大，電暈線處的空間電荷密度增大；沈恒對(duì)星形線進(jìn)行了數(shù)值模擬發(fā)現(xiàn)星形線電除塵通道內(nèi)二次流動(dòng)對(duì)流動(dòng)形態(tài)和粒子濃度分布存在顯著作用；裴藝凱利用數(shù)值模擬發(fā)現(xiàn)主流速度的增大減弱了ｅｈｄ流對(duì)電除塵器內(nèi)流場(chǎng)的干擾；張哲當(dāng)電場(chǎng)斷面風(fēng)速越小、粒徑越大、施加電壓越大時(shí)，除塵效率越高，對(duì)粒徑為6.9μm的粉塵最大捕集效率可以達(dá)到99.4％；高韜深對(duì)濕式電除塵器進(jìn)行了數(shù)值模擬，發(fā)現(xiàn)與干式狀態(tài)相對(duì)比，得出以下結(jié)論：濕式工作狀態(tài)下粒子的分段濃度減少率、沿程沉積效率、除塵總效率的變化趨勢(shì)與干式時(shí)相同，但數(shù)值上更大，除塵效果更好；劉凱對(duì)濕式電除塵器進(jìn)行了數(shù)值模擬，發(fā)現(xiàn)添加氣流均布板有利于電除塵器內(nèi)部流場(chǎng)分布的優(yōu)化。但是國(guó)內(nèi)外學(xué)者們對(duì)收塵極板的研究很少。

本文利用comsolmultiphysics數(shù)值軟件模擬了不同板型電除塵器模型在不同風(fēng)速、不同電壓、不同粒徑和不同溫度下的電勢(shì)分布、流場(chǎng)分布和除塵效率，旨在探究板型對(duì)離子風(fēng)和電除塵器性能的影響因素，優(yōu)化電除塵器的配置。

1 數(shù)值模型

高壓靜電除塵過(guò)程是電場(chǎng)、流場(chǎng)、顆粒荷電及運(yùn)動(dòng)場(chǎng)的多物理耦合過(guò)程。

1.1 靜電場(chǎng)模型

電場(chǎng)模型所涉及到的參數(shù)主要包括電勢(shì)與空間電荷密度等。其中電場(chǎng)分布由泊松方程和電流連續(xù)方程共同控制，已在文獻(xiàn)有了充分介紹。

1.2 氣流流動(dòng)模型

本文用ｋ－ε模型來(lái)描述湍流場(chǎng)，在有限元分析軟件comsolmultiphysics中進(jìn)行求解。其中ｎavier－ｓtokes方程中的質(zhì)量守恒方程和動(dòng)量守恒方程求解流速ｕ和壓力ｐ，標(biāo)準(zhǔn)ｋ－ε湍流模型求解湍動(dòng)能和耗散率。

1.3 粒子荷電模型

粒子的荷電方式有兩種分為電場(chǎng)荷電和擴(kuò)散荷電對(duì)于直徑大于1μm的顆粒，電場(chǎng)荷電通常大于擴(kuò)散荷電。粒子荷電量與粒子直徑、電場(chǎng)強(qiáng)度等因素有關(guān)，粒子荷電方程由文獻(xiàn)給出。

1.4 幾何模型

本文建立了四種收塵極板板型的除塵器并對(duì)其進(jìn)行了二維數(shù)值計(jì)算，其幾何尺寸如圖1所示。圖1（ａ）、（ｂ）、（ｃ）、（ｄ）分別為平行板收塵極板

圖1 幾何模型

（后續(xù)簡(jiǎn)稱ｐ型板）、ｂｅ型板、m型板、ｃｗ型板幾種電除塵器幾何模型。三種模型尺寸均為500mm×100mm，電極線采用半徑2mm的圓柱形電極，均為三線式布置。三根電極線在圖中坐標(biāo)（ｘ，ｙ）分別為：（83.3mm，0）、（250mm，0）、（416.6mm，0）。相較于平行板，其它類型的極板的板面場(chǎng)強(qiáng)分布和板面電流分布比較均勻，能夠在流場(chǎng)內(nèi)使氣流緊縮和氣流擴(kuò)大，進(jìn)而影響離子風(fēng)渦旋對(duì)流場(chǎng)擾動(dòng)的強(qiáng)弱，同時(shí)其它極板內(nèi)的凹槽較ｐ型板更容易捕集一些沿氣流方向運(yùn)動(dòng)的粉塵，從而提高電除塵器的除塵效率。本文數(shù)值計(jì)算選取的的參數(shù)值為：空氣粘度μ＝1.789×10－5ｋｇ/（m·ｓ）；空氣中離子遷移率ｋ＝2.2×10－4m2/（ｖ·ｓ）；離子擴(kuò)散系數(shù)ｄ＝5.3×10－5m2/ｓ。

2 模擬結(jié)果

不同收塵極板與電極線構(gòu)成了多種形式極配，收塵極板形狀對(duì)離子風(fēng)渦旋和除塵器內(nèi)流場(chǎng)具有不同的作用，進(jìn)而對(duì)顆粒物的捕獲有著不同的影響。本章建立了二維線板電除塵器的數(shù)值模型，引入ｂｅ型板、ｃｗ型板和m型板與ｐ型板對(duì)比，詳細(xì)研究了在三種板型下除塵器內(nèi)部的靜電場(chǎng)分布和離子風(fēng)湍流場(chǎng)及顆粒物荷電的規(guī)律，進(jìn)一步探尋在電流體動(dòng)力作用下顆粒運(yùn)動(dòng)及捕獲過(guò)程。

2.1 電勢(shì)分布

圖2電壓50kv時(shí)，三種板型電除塵器內(nèi)的電勢(shì)分布，通過(guò)在從電極線位置開(kāi)始ｙ方向ｙ＝0、10、20、30、40mm提取數(shù)據(jù)點(diǎn)，繪制分布圖。在圖2中可以看出，除塵器內(nèi)部電勢(shì)分布不均勻，但三種板型下的等勢(shì)線均以電極線為中心呈橢圓狀向外擴(kuò)散，電勢(shì)向收塵極板方向變化速度更大。在入口到出口過(guò)程中電勢(shì)分布在電極線處達(dá)到峰值，電極線與電極線之間有一低谷。由于收塵極板形狀不同，內(nèi)部電勢(shì)分布略有差異；由圖2可以看出電勢(shì)波動(dòng)變化最為明顯的是ｂｅ型板，然后是ｃｗ型板、m型板和ｐ型板。因?yàn)椴痪鶆虻碾妶?chǎng)分布更有利于除塵效率的提高，而由圖4可以對(duì)不同極板不均勻程度作比較，可以看出ｂｅ板＞ｃｗ型板＞m型板＞ｐ型板。電勢(shì)分布有差異說(shuō)明了收塵極板形狀對(duì)電勢(shì)分布有影響，但隨著趨近電極線，不同收塵極板板型對(duì)電勢(shì)分布的影響減弱。

圖2 電勢(shì)分布

2.2 湍流場(chǎng)分布

靜電除塵器的離子風(fēng)作用于流場(chǎng)產(chǎn)生更復(fù)雜的湍流，會(huì)影響到荷電顆粒的捕集，離子風(fēng)對(duì)提高顆粒收集效率方面起著重要的作用。以ｐ型板電除塵器為例進(jìn)行分析。由圖3（ａ）可以看到電力線是以電極線為中心指向收塵極板并垂直于電勢(shì)等值線（電勢(shì)等值線以電極為中心呈橢圓狀向外延展）。電力線與電力線之間并不會(huì)相交，所以可以看到在相鄰電極線之間存在著一由電力線之間相互隔離出來(lái)的隱形壁面，離子風(fēng)的產(chǎn)生原因是離子與空氣分子碰撞的動(dòng)量傳遞，進(jìn)而帶動(dòng)空氣流動(dòng)。把空氣中的分子或離子、單個(gè)細(xì)顆粒假設(shè)為一個(gè)小球體，在無(wú)主氣流流動(dòng)時(shí)，流體中與電極線極性相同的帶電離子或荷電細(xì)顆粒物，在電極線附近受到ｘ方向和垂直指向收塵極板ｙ方向的電場(chǎng)力ｆｘ和ｆｙ，不同的力協(xié)同作用下產(chǎn)生位移進(jìn)而帶動(dòng)空氣分子運(yùn)動(dòng)。

在電極線前ｘ方向球體運(yùn)動(dòng)與主氣流流動(dòng)方向相反，在電極線后ｘ方向球體運(yùn)動(dòng)與主氣流流動(dòng)方向相同，而ｙ方向受到的力使其撞擊到收塵極板面，由于空氣流體不可壓縮，最終在ｆｘ和ｆｙ和壁面共同作用下在電極線前（面向入口）產(chǎn)生逆時(shí)針渦旋，電極線后（面向出口）產(chǎn)生順時(shí)針渦旋。隨著主氣流速度不斷的增大，流動(dòng)氣體相對(duì)運(yùn)動(dòng)阻力ｆｄ會(huì)逐漸抵消部分與氣流方向相反的力，在電極線后的ｆｘ、ｆｙ與ｆｄ形成向上向前的合力，最終導(dǎo)致離子風(fēng)渦旋的逐漸減弱。這也說(shuō)明主氣流風(fēng)速?zèng)Q定了離子風(fēng)減弱速度，更高的風(fēng)速有更強(qiáng)的流體阻力，渦旋減弱越快。

圖3（ｂ）繪制了假設(shè)球點(diǎn)的受力運(yùn)動(dòng)軌跡與離子風(fēng)渦旋形成與消散示意圖，圖3（ｃ）為在電壓20kv，入口風(fēng)速0.05m/ｓ時(shí)ｐ型板電除塵器內(nèi)湍流流線圖，并標(biāo)出了渦旋流線方向，以旋轉(zhuǎn)箭頭表示。線板式靜電除塵器中一定數(shù)量的收塵極板相互之間排列組成除塵通道，板型的變化會(huì)對(duì)通道離子風(fēng)渦旋產(chǎn)生影響。

圖4是ｐ型板、ｂｅ型板、m型板、ｃｗ型板幾種板型下電除塵器在電壓為50kv，粒子直徑為2.5μm，風(fēng)速為0、0.2和0.5m/ｓ流速下的湍流場(chǎng)的流線圖。由圖4可以看出，在風(fēng)速為0時(shí)，大部分除塵器流場(chǎng)內(nèi)以除塵器中心線為中心分布著12個(gè)渦旋，每個(gè)電極線周圍分布著4個(gè)渦旋，而且上下分布的渦旋流線方向相反，最中心電極產(chǎn)生的離子風(fēng)對(duì)兩邊電極有一定的抑制作用；而ｂｅ型板則只有11個(gè)渦旋，m型板和ｃｗ板的離子風(fēng)旋渦分布不規(guī)則，這說(shuō)明除塵器板型對(duì)于除塵器內(nèi)流場(chǎng)分布有著很大的影響，同時(shí)電勢(shì)分布的不均勻?qū)Τ龎m器湍流場(chǎng)也造成了影響。隨著風(fēng)速的增大，當(dāng)主氣流風(fēng)速為0.2m/ｓ時(shí)，由于流體相對(duì)運(yùn)動(dòng)阻力增大，渦旋與主氣流相互作用逐漸消失，一部分渦旋消失。

圖3 離子風(fēng)產(chǎn)生原理

由圖4可以看出，這五個(gè)除塵器模型的湍流場(chǎng)流線分布比較類似，由于主氣流風(fēng)速的影響，流場(chǎng)中靠近進(jìn)口處的電極線右側(cè)的兩個(gè)渦旋被主氣流沖散，而且其它渦旋也都因?yàn)橹鳉饬鞯淖饔眯螤詈兔娣e發(fā)生了改變。當(dāng)主氣流風(fēng)速達(dá)到0.5m/ｓ時(shí)，主氣流在流場(chǎng)中起主要作用，所有板型的除塵器模型流場(chǎng)中的渦旋都被主氣流沖散。

圖5為相同條件下，四種板型在ｙ＝0mm處沿ｘ方向的速度分布對(duì)比。由圖5可以看出，在ｙ＝0mm，主氣流速度為0.2m/ｓ時(shí)沿氣流方向從入口到出口過(guò)程中四種板型除塵器內(nèi)速度變化幅度規(guī)律基本相似，速度變化均呈波浪狀曲線分布，在兩根電極線之間存在速度低谷。而各極板速度的峰值都在極線處，四種極板在電極線附近最大值分布的規(guī)律是：第一根＞第二根＞第三根?？梢钥闯鏊俣却笮‰Sｘ軸波動(dòng)變化大小的排序是：ｂｅ型板＞ｃｗ型板＞m型板＞ｐ型板，另外三種板型波動(dòng)大于ｐ型板的原因是，它們的形狀不規(guī)則，有些地方通道窄，有些地方通道寬，這樣更有利于顆粒物的捕集。由圖中第二根和第三根極線之間的波動(dòng)和第三根后的波動(dòng)可以看出，這兩個(gè)波動(dòng)的最高點(diǎn)都是處于上下兩個(gè)離子風(fēng)旋渦的交匯處，波動(dòng)大小總體上ｂｅ板＞ｃｗ型板＞m型板＞ｐ型板，這也代表了離子風(fēng)的強(qiáng)弱，波動(dòng)大的板型離子風(fēng)效應(yīng)更強(qiáng)，主氣流對(duì)其的影響更弱。

圖4 不同極板在不同風(fēng)速下的湍流場(chǎng)流線圖

圖5 沿ｘ軸速度分布對(duì)比

2.3 除塵效率

除塵效率是判斷一個(gè)電除塵器性能優(yōu)劣的主要標(biāo)準(zhǔn)，電除塵器的顆粒物捕集效率受到很多因素的影響。圖6是利用數(shù)值軟件ｃｏmｓｏｌmｕｌ?ｔｉｐｈｙｓｉｃｓ和控制變量法模擬的ｐ型板、ｂｅ型板、ｃｗ型板和m型板在不同一次風(fēng)速、不同電壓、不同粒徑和不同溫度下的除塵效率的對(duì)比。圖6表明不同極板隨著一次風(fēng)速的降低、放電極電壓的升高、顆粒物粒徑的增大和溫度的降低，它們的除塵效率隨之升高，同時(shí)它們的增長(zhǎng)率也是隨之升高。不同工況下不同極板的除塵效率作比較：ｂｅ型板＞ｃｗ型板＞m型板＞ｐ型板，根據(jù)前文總結(jié)，這是因?yàn)椴痪鶆虻碾妶?chǎng)分布、速度大小隨ｘ軸波動(dòng)變化大以及較強(qiáng)的離子風(fēng)渦旋更有利于粒子的捕集。

3 結(jié)論

本文對(duì)三種板型的靜電除塵器進(jìn)行數(shù)值計(jì)算，對(duì)不同板型條件下除塵器內(nèi)離子風(fēng)產(chǎn)生的湍流場(chǎng)和顆粒物的荷電與捕獲過(guò)程及其顆粒運(yùn)動(dòng)軌跡進(jìn)行了研究，并給出了通過(guò)改善離子風(fēng)提高對(duì)顆粒捕獲的理論分析，主要結(jié)論如下：收塵極板板型的研究一直以來(lái)受到國(guó)內(nèi)外學(xué)者的青睞，研究者發(fā)現(xiàn)通過(guò)極板形狀的改變可以提高電除塵器的除塵效率。所以本文分析比較了ｐ型版、ｂｅ型板、ｃｗ型板和m型板在有ｅｈｄ流情況下的電勢(shì)分布、流場(chǎng)分布和除塵效率，來(lái)探究板型對(duì)除塵器性能的影響。由不同板型除塵器的電勢(shì)分布可以看出，收塵極板形狀對(duì)電勢(shì)分布有影響，電勢(shì)不均勻程度：ｂｅ型板＞ｃｗ型板＞m型板＞ｐ型板，且隨著趨近電極線，不同極板板型對(duì)電勢(shì)分布的影響減弱；而從湍流場(chǎng)可以看出，流場(chǎng)內(nèi)的主氣流決定了離子風(fēng)的減弱速度，更高的流速有更強(qiáng)的流體阻力，渦旋減弱的更快，同時(shí)當(dāng)主氣流風(fēng)速為0時(shí)，最中心電極產(chǎn)生的離子風(fēng)對(duì)兩邊電極有一定的抑制作用，當(dāng)主氣流速度為0.2m/ｓ時(shí)，速度大小隨著ｘ軸波動(dòng)變化大小的排序是：ｂｅ板＞ｃｗ型板＞m型板＞ｐ型板，不均勻的電場(chǎng)會(huì)產(chǎn)生更強(qiáng)的離子風(fēng)，波動(dòng)大的離子風(fēng)效應(yīng)更強(qiáng)，主氣流對(duì)其的影響更弱。

最后對(duì)不同板型除塵器的除塵效率進(jìn)行了比較，發(fā)現(xiàn)不同極板隨著一次風(fēng)速的降低，放電極電壓的升高、顆粒物粒徑的增大和溫度的降低其除塵效率隨之升高，而在不同工況下除塵效率：ｂｅ板＞ｃｗ型板＞m型板＞ｐ型板。從不同板型電勢(shì)不均勻程度、湍流場(chǎng)速度分布、沿ｘ軸速度大小和除塵效率的比較可以發(fā)現(xiàn)，不均勻的電勢(shì)分布會(huì)產(chǎn)生更強(qiáng)的離子風(fēng)，而更強(qiáng)的離子風(fēng)會(huì)使除塵器的除塵效率提高。

圖6 不同極板在不同一次風(fēng)速、不同電壓、不同粒徑和不同溫度下的除塵效率a