三維絲網(wǎng)網(wǎng)胞特征參數(shù)的數(shù)值模擬

王曉靜 樊俊曉 王 菁

（天津大學(xué)a.化工學(xué)院；b.環(huán)境科學(xué)與工程學(xué)院）

隨著我國社會的發(fā)展，工業(yè)廢氣的排放量日益增加，而未經(jīng)處理的廢氣直接排入大氣不僅會嚴(yán)重威脅人類的身體健康，同時較差的空氣質(zhì)量也會影響人類的心理健康［1］。因此，加強(qiáng)工業(yè)廢氣治理，落實環(huán)保理念，提高工業(yè)廢氣防治技術(shù)，實現(xiàn)可持續(xù)發(fā)展已刻不容緩［2］。

有色煙羽中含多種污染物，其中可溶性顆粒物和鹽類是導(dǎo)致霧霾天氣的原因之一，還可能含有可凝結(jié)顆粒物，這種顆粒物停留時間長、擴(kuò)散距離遠(yuǎn)、毒性強(qiáng)、危害大［3］。 揮發(fā)性有機(jī)物（VOCs）是工業(yè)廢氣的主要產(chǎn)物之一，也是導(dǎo)致霧霾天氣的主要因素。 目前，人們主要從兩個方面解決大氣污染問題：一是提高工藝水平，升級環(huán)保設(shè)施；二是生產(chǎn)高質(zhì)量產(chǎn)品，以減少污染廢氣的排放［4］。主要的工藝脫白方法有3種，包含：直接加熱法、先降溫再加熱法和降溫混風(fēng)法［5］。 然而，這3種脫白方法對于大風(fēng)量脫白處理而言成本較高。 目前常見的VOCs治理技術(shù)，傳統(tǒng)的主要有直接燃燒法和催化燃燒法，新興的主要有活性炭纖維治理技術(shù)、生物治理技術(shù)及蓄熱式熱力氧化技術(shù)等［6］。直接燃燒法對操作有嚴(yán)格的要求，如果焚燒溫度不達(dá)標(biāo)，容易產(chǎn)生其他污染物；催化燃燒法使用成本高，投資量大；活性炭纖維治理技術(shù)的使用環(huán)境如果過于潮濕會影響處理效果；生物治理技術(shù)尚未確立明確的使用方法；蓄熱式熱力氧化技術(shù)可能會產(chǎn)生碳氧化合物，造成大氣二次污染。

裝備是化工行業(yè)組織生產(chǎn)的基礎(chǔ)，提高化工設(shè)備性能，對提高化工生產(chǎn)過程和化工技術(shù)與產(chǎn)業(yè)整體水平至關(guān)重要［7］。 三維網(wǎng)胞形式的立體絲網(wǎng)是一種立體編制的復(fù)合結(jié)構(gòu)絲網(wǎng)，可通過捕集微小液滴并使之凝聚，達(dá)到脫白的效果，進(jìn)而在絲網(wǎng)上燒結(jié)附著催化劑，實現(xiàn)催化分解有害污染物的目的。

筆者主要研究不同絲網(wǎng)參數(shù)對水滴的捕集效率和對流體阻力的影響，為立體絲網(wǎng)工業(yè)化應(yīng)用提供參考依據(jù)。

1 模型構(gòu)建及網(wǎng)格劃分

三維絲網(wǎng)是V形結(jié)構(gòu)在經(jīng)、 緯方向上垂直交叉構(gòu)成的。 利用NX12.0建立三維絲網(wǎng)模型，通過參數(shù)驅(qū)動的功能， 輸入不同絲網(wǎng)直徑D、 目數(shù)M（絲網(wǎng)水平投影，每英寸含有的絲網(wǎng)個數(shù)）和層數(shù)N，自動生成不同的模型，從而大幅減少模型的構(gòu)造時間。 由于三維絲網(wǎng)具有對稱性，因此筆者選擇四分之一結(jié)構(gòu)進(jìn)行建模，得到的三維絲網(wǎng)模型如圖1所示。

對氣體流過絲網(wǎng)的區(qū)域建立流場區(qū)域，對進(jìn)、出口進(jìn)行一定的延長，以達(dá)到穩(wěn)定氣流的目的。 對絲網(wǎng)模型均采用非結(jié)構(gòu)化網(wǎng)格劃分，對絲網(wǎng)接觸部分進(jìn)行局部加密處理（圖2）。

圖1 三維絲網(wǎng)模型

圖2 三維絲網(wǎng)模型網(wǎng)格劃分

2 模擬方法

2.1 控制方程

流體流動需要滿足物理守恒定律， 包括質(zhì)量守恒定律、動量守恒定律和能量守恒定律，其控制方程分別為連續(xù)方程、 運(yùn)動方程和能量方程［8］。

2.1.1 連續(xù)方程

連續(xù)方程是質(zhì)量守恒定律以數(shù)學(xué)表達(dá)式的形式在運(yùn)動流體中的運(yùn)用，其表達(dá)式為：

式中 t——時間；

ū——速度；

ρ——流體密度。

對于不可壓縮定長流動，且密度沒有發(fā)生變化的流體，其連續(xù)方程為：

由于本模擬不涉及傳熱過程，幾乎沒有溫度變化，因此不涉及能量方程。

2.2 湍流模型選擇

Fluent提供的湍流運(yùn)動計算方法主要有雷諾時均法和大渦模擬法［9］。 使用雷諾時均法計算湍流運(yùn)動的模型主要有Spalart-Allmaras模型、 標(biāo)準(zhǔn)k-ε模型、RNG k-ε模型及Realizable k-ε模型等。其中Spalart-Allmaras模型的穩(wěn)定性較差， 尚無成功的案例；標(biāo)準(zhǔn)k-ε模型是由Launder和Spalding提出的，主要用于高雷諾數(shù)的湍流模型計算，穩(wěn)定性高，但對于非均勻湍流的計算，有較大的誤差［10］；RNG k-ε模型是對標(biāo)準(zhǔn)k-ε模型的優(yōu)化，強(qiáng)于標(biāo)準(zhǔn)k-ε模型；Realizable k-ε模型使用的是新的湍流粘度公式， 在雷諾應(yīng)力上和真實湍流保持一致，在復(fù)雜的湍流模型中更出色。 由于本模型是一種多孔介質(zhì)模型，湍流流動比較復(fù)雜，因此筆者選用Realizable k-ε模型。

2.3 除霧機(jī)理

除霧機(jī)理［11］主要有以下3種：

a. 慣性碰撞。 液滴保持原來的運(yùn)動狀態(tài)，在氣體繞過絲網(wǎng)時，由于液滴的慣性會撞擊在絲網(wǎng)表面，從而實現(xiàn)捕集。

b. 布朗擴(kuò)散。 當(dāng)氣流較低或液滴粒徑較小時，液滴做不規(guī)則運(yùn)動從而碰撞到絲網(wǎng)表面。

c. 直接攔截。 粒徑較小的或者氣速較小的液滴隨著氣流繞過絲網(wǎng)，當(dāng)液滴的半徑大于液滴圓心到絲網(wǎng)的距離時將被捕集。

本課題選用的水滴直徑為10μm，慣性碰撞是主要的影響因素，布朗擴(kuò)散和直接攔截影響非常小，可忽略不計［12］。 當(dāng)水滴與絲網(wǎng)表面發(fā)生碰撞，即認(rèn)為是捕集到了水滴。

壓降△p的計算式為：

2.4 邊界條件設(shè)置

進(jìn)口處設(shè)置為速度入口，入口速度5m/s，DPM設(shè)置為液滴入口；出口處設(shè)置為壓力出口，靜壓為0，DPM設(shè)置為逃逸；相鄰兩個矩形面設(shè)置為壁面，DPM設(shè)置為反射； 另兩個相鄰矩形面設(shè)置為對稱面；絲網(wǎng)表面設(shè)置為捕集。 本課題選用的氣體為空氣、液體為水。

3 模擬結(jié)果與分析

3.1 絲網(wǎng)直徑對壓降和捕集效率的影響

絲網(wǎng)直徑的變化會對壓降和捕集效率產(chǎn)生較大的影響。 對3種系列（D×31.8×10-D×31.8×10、D×25.4×10-D×25.4×10和D×21.2×10-D×21.2×10）不同直徑的絲網(wǎng)進(jìn)行模擬研究，數(shù)據(jù)繪制成曲線如圖3、4所示。

圖3 絲網(wǎng)直徑與壓降的關(guān)系曲線

圖4 絲網(wǎng)直徑與捕集效率的關(guān)系曲線

由圖3、4可以看出，隨著絲網(wǎng)直徑的增加，壓降呈現(xiàn)線性增加的趨勢，這是因為隨著絲網(wǎng)直徑的增加，V形邊上相鄰兩個絲網(wǎng)之間的氣流速度變大，絲網(wǎng)表面附近的流體湍動能增加，湍流強(qiáng)度增加， 進(jìn)而導(dǎo)致對空氣產(chǎn)生的阻力有所增加。同時，由于絲網(wǎng)直徑的增加，提高了對水滴的有效攔截面積利用率，使得水滴的捕集效率得到了提高。 當(dāng)D≥1/M時，絲網(wǎng)對水滴的捕集效率幾乎不發(fā)生變化， 但是壓降繼續(xù)呈現(xiàn)線性增加的趨勢，這是由于隨著絲網(wǎng)直徑的增加，湍流強(qiáng)度還會逐漸增加，但是對水滴的有效攔截面積利用率卻不再發(fā)生變化。對于D×31.8×10-D×31.8×10系列絲網(wǎng)，當(dāng)絲網(wǎng)直徑D=1.5mm時，壓降發(fā)生了劇烈的變化，達(dá)到了83.3Pa。 這是因為倒數(shù)第2層相鄰的兩個絲網(wǎng)距離非常近，間隙非常小，導(dǎo)致對空氣的粘性阻力很大，進(jìn)而造成整個模型的壓降迅速增加。

當(dāng)D≥1/M時， 雖然三維絲網(wǎng)的捕集效率變化已不大，但壓降還在繼續(xù)變化，因此筆者對D=1/M的三維絲網(wǎng)進(jìn)行進(jìn)一步模擬分析。 對于D×25.4×10-D×25.4×10系列絲網(wǎng)，D=1/M時直徑與壓降和捕集 效 率 的 關(guān) 系 曲 線 如 圖5、6 所 示。 由 圖5、6可以看出，當(dāng)D=1/M時，隨著絲網(wǎng)直徑的增加，壓降大幅變小，捕集效率整體呈下降趨勢但變化不大。 說明水滴的有效攔截面積利用率不變時，三維絲網(wǎng)對水滴的捕集效率也不會產(chǎn)生大的變化。

圖5 D=1/M時直徑與壓降的關(guān)系曲線

圖6 D=1/M時直徑與捕集效率的關(guān)系曲線

3.2 絲網(wǎng)目數(shù)對壓降和捕集效率的影響

絲網(wǎng)目數(shù)的變化也會對壓降和水滴捕集效率產(chǎn)生影響。 對3種系列 （0.8×M×10-0.8×M×10、1.0×M×10-1.0×M×10和1.2×M×10-1.2×M×10）不同目數(shù)的絲網(wǎng)進(jìn)行模擬研究， 數(shù)據(jù)繪制成如圖7、8所示曲線。

圖7 絲網(wǎng)目數(shù)與壓降的關(guān)系曲線

由圖7、8可以看出，隨著絲網(wǎng)目數(shù)的增加，壓降也是呈現(xiàn)增加的趨勢。當(dāng)1/M＞D時，水滴捕集效率隨著目數(shù)的增加而增加；當(dāng)1/M≤D時，目數(shù)的變化對水滴捕集效率的影響較小。 這是因為隨著目數(shù)的增加，絲網(wǎng)的密集度增加，孔隙率減少，導(dǎo)致壓降增加，而水滴的有效攔截面積利用率先增加然后不再發(fā)生變化。

圖8 絲網(wǎng)目數(shù)與捕集效率的關(guān)系曲線

3.3 絲網(wǎng)層數(shù)對壓降和捕集效率的影響

絲網(wǎng)層數(shù)的變化也會對壓降和水滴捕集效率產(chǎn)生影響。對3種系列（0.8×31.8×N-0.8×31.8×N、1.0×25.4×N-1.0×25.4×N和1.2×21.2×N-1.2×21.2×N）不同層數(shù)的絲網(wǎng)進(jìn)行模擬研究，數(shù)據(jù)繪制成如圖9、10所示曲線。 可以看出，隨著絲網(wǎng)層數(shù)的增加，壓降逐漸減小，而捕集效率變化不大，主要在98%～100%之間波動。這是因為隨著層數(shù)的增加，絲網(wǎng)的孔隙率增大，流過相鄰網(wǎng)孔的氣流干擾性降低，有助于氣流的穩(wěn)定，利于湍流的擴(kuò)散，從而使湍流強(qiáng)度降低，整體壓降降低；而絲網(wǎng)層數(shù)的增加并沒有改變對水滴的有效攔截面積利用率，同時由于模擬產(chǎn)生的系統(tǒng)誤差，所以捕集效率僅產(chǎn)生小幅的波動性變化。

圖9 絲網(wǎng)層數(shù)與壓降的關(guān)系曲線

圖10 絲網(wǎng)層數(shù)與捕集效率的關(guān)系曲線

4 結(jié)論

4.1 隨著絲網(wǎng)直徑的增加，壓降增加。 當(dāng)D＜1/M時， 水滴的捕集效率隨著絲網(wǎng)直徑的增加而增加；當(dāng)D≥1/M時，絲網(wǎng)直徑的變化對水滴的捕集效率影響不大。 當(dāng)D=1/M時， 隨著絲網(wǎng)直徑的增加，壓降減小，捕集效率變化范圍較小且整體呈下降趨勢。

4.2 隨著絲網(wǎng)目數(shù)的增加， 壓降呈現(xiàn)增加的趨勢。當(dāng)1/M＞D時，水滴的捕集效率隨著絲網(wǎng)目數(shù)的增加而增加；當(dāng)1/M≤D時，絲網(wǎng)目數(shù)的變化對水滴的捕集效率影響不大。

4.3 隨著絲網(wǎng)層數(shù)的增加，孔隙率增加，壓降減小， 但是沒有改變水滴的有效攔截面積利用率， 因此水滴的捕集效率沒有發(fā)生顯著的變化。