吳興祥 劉冬 原光輝 張亮

通訊作者：劉冬，1997年5月，男，漢族，安徽池州人，就讀于上海海事大學(xué)，碩士。研究方向：電氣工程。

基金項目：港口機械與海工裝備智能裝配關(guān)鍵技術(shù)研究及其示范應(yīng)用（18040501600）;臨港新片區(qū)高新產(chǎn)業(yè)和科技創(chuàng)新專項項目（SH-LG-GK-2020-37）。

摘 要：通過對揮發(fā)性有機化合物（VOCs）收集處理裝置三維建模，并運用fluent流體數(shù)值分析方法，對VOCs催化設(shè)備內(nèi)流場進行數(shù)值分析，催化床用多孔介質(zhì)等效替換，其結(jié)果表明，其內(nèi)置流場速度分布與入口布置，流域形狀，多孔介質(zhì)層數(shù)直接相關(guān)。合理的入口布置與流域形狀有利流場流速的均勻性，多孔介質(zhì)的層數(shù)越多其流速均勻性越好。

關(guān)鍵詞：VOCs;fluent;流域形狀;多孔介質(zhì);流速均勻性;

一、引言

油漆間工作時會產(chǎn)生大量的揮發(fā)性有害氣體，會不同程度影響當(dāng)?shù)丨h(huán)境和人們的生活環(huán)境水平，對于有害氣體的收集和處理顯得愈發(fā)重要，同時新的公約標(biāo)準對工業(yè)區(qū)有害氣體的排放有進一步要求，針對這些要求，公司研發(fā)生產(chǎn)了一套VOCs（volatile organic compounds）收集處理裝置。為了更可靠更方便更有科學(xué)依據(jù)說明此套設(shè)備的性能，現(xiàn)對其展開研究與分析。

而催化設(shè)備在整套設(shè)備中是關(guān)鍵設(shè)備之一，催化設(shè)備工作性能的好壞直接影響整套設(shè)備的性能。催化設(shè)備中催化劑[1，2]為貴金屬整體式蜂窩型模式，模擬時采用多孔介質(zhì)形式來引入此模型。VOCs與空氣混合氣體從裝置進口流入，經(jīng)過蜂窩型催化劑，再經(jīng)過熱交換器流出。其冷卻氣體進出口就是裝置從油漆間收集后的混合氣體，在熱交換器處先進行預(yù)加熱，同時對催化裝置出口起到冷卻的作用。而催化劑的催化效果與處理混合氣體的含VOCs氣體的均勻程度，以及流經(jīng)催化劑區(qū)域的流速分布情況有關(guān)，不同流速會影響催化劑的催化效率[3]，同時流速分布不均勻，會導(dǎo)致整體式催化劑的利用效率下降。

為了研究其催化裝置的流場特性[4，5，6]，需要對其展開分析與研究，而采用fluent數(shù)值分析的方法，會更加直觀了解整個裝置流場分布和工作特性。

二、模型

為了準確地反應(yīng)實際裝置工作時的流場特性，需要對裝置流場區(qū)域進行三維建模，根據(jù)其裝置實際情況來搭建三維裝置流域[7]。同時為了減少相應(yīng)的無關(guān)計算量，對其局部等小細節(jié)對裝置流場分析影響很小的局部區(qū)域進行簡化、優(yōu)化。同時對催化區(qū)域通過引入多孔介質(zhì)模型[8]來等效替代其在流場中對流體流動特性的效果。根據(jù)實際裝置建立以下相關(guān)的流道區(qū)域模型，如圖1所示。

如圖1（a）所示為催化床的邊界情況，包括混合氣體進口、混合氣體出口、冷卻氣體進出口和多孔介質(zhì)模擬蜂窩型催化劑。

如圖1（b）所示為催化床的尺寸參數(shù)，尺寸信息如上圖所示，總高2337 mm，最大寬度860 mm，催化劑區(qū)域為255 mm。下部為熱交換器，冷混合氣體從換熱板間隔空隙中流過，催化燃燒后氣體從換熱板內(nèi)部區(qū)域流過，再從混合氣體出口排出。

如圖1所示，均直接定義催化床的內(nèi)部流體區(qū)域模型，外部殼體不顯示在流體模型場景中。

三、Fluent邊界條件確立

整個系統(tǒng)的流量與抽氣機的風(fēng)量直接相關(guān)，假設(shè)在理想狀況下的系統(tǒng)流量完全由抽氣機的風(fēng)量決定，設(shè)定整體系統(tǒng)的入口風(fēng)量：

Q風(fēng)機的風(fēng)量，Sin入口截面面積。將5000 m3/h，以及半徑r = 0.075 m，代入式中，求得Sin = 0.0177 m2，解得υV = 78.6 m/s;見圖2。

根據(jù)實際測得的VOCs濃度，來定義混合氣體的各組分含量，在此實例中，選用甲苯來代表整個VOCs的組分，來模擬催化處理裝置流道的流動分布情況。測得室內(nèi)VOCs氣體濃度為200 mg/m3，把其換算為進口質(zhì)量分數(shù)。在標(biāo)準狀況下一立方米空氣質(zhì)量大約為1.29 kg/m3，其進口甲苯質(zhì)量分數(shù)為200 mg/1.29 kg = 0.0155%，見圖3。

四、網(wǎng)格無關(guān)性分析

在fluent數(shù)值分析中，通過把流域離散為眾多單元區(qū)域的形式來計算其流場特性，而數(shù)值離散的方法會使其解算結(jié)果與真實值存在誤差。所以離散單元的數(shù)目會直接影響求解精度，在理論上，網(wǎng)格數(shù)目越多其求解精度越高，但是過高的求解精度會導(dǎo)致計算量無限制的增大，影響求解速度和電腦的資源配置。在滿足一定精度的網(wǎng)格數(shù)目和優(yōu)良的求解速度之間取得一個平衡就至關(guān)重要。研究單元網(wǎng)格數(shù)與出口流速的無關(guān)性。

如表1所示，當(dāng)網(wǎng)格尺寸達到20 mm后時，其出口流速變化趨于穩(wěn)定。

五、流場分析

假設(shè)VOCs收集處理裝置混合氣體經(jīng)過緩沖，過濾吸附裝置到達催化裝置入口時，已經(jīng)混合均勻，且比較穩(wěn)定。在此文分析中著重針對催化裝置內(nèi)部流場速度分布進行分析。因為速度的分布情況直接，影響催化裝置的催化效果包括催化轉(zhuǎn)換效率，以及整體催化劑的利用率。同時為催化劑的布置以及VOCs催化氧化仿真模擬提供理論依據(jù)。

如圖4所示為催化裝置內(nèi)部整體流速分布情況。催化裝置進口流速很高達到80 m/s，且入口處速度波動比較劇烈，根據(jù)上述云圖觀察，顏色鮮明處，速度流速較大，此催化裝置多孔介質(zhì)上部氣體流速差異較大，局部流速明顯。由于進口的布置，會導(dǎo)致氣體流速分布不均勻。同時由于進口流速較高，導(dǎo)致其他低流速區(qū)域云圖顏色變化狀況并不明顯，但并不意味著流速分布均勻。下面對催化裝置催化區(qū)域進行局部分析。

如圖5與圖6所示，多孔區(qū)域由于入口布置的影響，遠離入口端的混合氣體流速比靠近入口段的混合氣體流速明顯較大。同時上層多孔區(qū)域局部流速速度差明顯高于下層多孔區(qū)域流速速度差。上層多孔區(qū)域局部流速最高流速大致為5.278 m/s，最低流速大約為1.3 m/s，已接近四倍的流速差異。而下層多孔區(qū)域局部流速最高流速大約為2.162 m/s，最低流速大約為1.5 m/s。由此可見，在多孔區(qū)域處流速分布不均，且局部差異較大。而下層多孔區(qū)域流速分布均勻性比上層多孔區(qū)域更加均勻。

為了觀察混合氣體流經(jīng)多孔介質(zhì)時其黏性阻力對速度分布的影響，在多孔區(qū)域Z軸方向上，每隔0.16米建立一個平面，從左到右總共5個平面用來檢測多孔區(qū)域上層混合氣體流速分布情況。

如圖7所示，下端的云圖顏色數(shù)目明顯小于上端的云圖數(shù)目，且下端云圖最高流速低于上端的最高流速，且顏色越深的區(qū)域，其流速等高線變化越明顯。說明在經(jīng)過多孔介質(zhì)時，由于存在粘性阻力，對流速區(qū)域均勻性存在積極的影響。同時流速越高時，其效果越明顯。

在多孔區(qū)域Z軸方向上，每隔0.16米建立一個平面，從左到右總共5個平面用來檢測多孔區(qū)域下層混合氣體流速分布情況。

如圖8所示，在上層多孔區(qū)域的流體阻力的影響下，流速等高線變化不大，同時其云圖顏色跨度也進一步下降，下層多孔介質(zhì)中的速度分布均勻性明顯上升。

為了觀察混合氣體流經(jīng)多孔介質(zhì)時速度局部分布的情況，在多孔區(qū)域Y軸方向上，每個多孔區(qū)域建立三個平面，從上到下總共6個平面用來檢測多孔區(qū)域混合氣體流速分布情況。

如圖9所示，受入口布置和流域形狀的影響，在上層多孔介質(zhì)區(qū)域，出現(xiàn)兩個局部速度較大的流域，其上層的流域均勻性較差，局部流域速度差異較大且變化明顯。

其下層多孔介質(zhì)也存在一個局部流速較大的區(qū)域，但相對于上層多孔介質(zhì)而言，其流速均勻性也得到了很大改善，其最大流速區(qū)域面積也大幅度增大，流速均勻性得到一定的保證。且速度變化范圍也大大降低。

六、結(jié)論

通過fluent數(shù)值分析的方法，對VOCs催化裝置進行仿真分析，并得到如下結(jié)果。

第一，在此VOCs收集處理裝置中，催化裝置的流速分布由于進口的布置情況，出現(xiàn)局部流速差異較大的情況，流速分布不均勻，遠離入口處流速比靠近入口處偏大。

第二，催化劑的多層布置，對流速均勻性起到積極的作用，由于其存在較大的粘性阻力，有利于增強流速的穩(wěn)定性和均勻性，改善催化劑區(qū)域的流域流場速度分布。

第三，同時流速分布受眾多因素的影響，包括入口分布，流域形狀，以及多孔介質(zhì)布置等諸多因素，故為了改善流場流速分布，可通過這些方面進一步改善流場分布。

參考文獻：

[1]殷珂，陳瑞洋，劉志明.錳基氧化物上甲苯催化氧化的研究進展[J].材料導(dǎo)報， 2020，34（23）：23051-23056.

[2]陳丹，田樹梅，石靜，沈華瑤.多孔氧化硅負載銀催化劑催化消除VOCs的研究進展[J].環(huán)境化學(xué)， 2020，39（11）：3145-3152.

[3]梁文俊，劉迪，任思達，李堅.基于Fluent的整體式催化反應(yīng)系統(tǒng)降解甲苯的數(shù)值模擬[J].環(huán)境工程學(xué)報， 2020，14（02）：457-464.

[4]向靖鋒，王立華，黃亞宇.基于Fluent的打葉機內(nèi)部流場數(shù)值模擬[J].農(nóng)業(yè)裝備與車輛工程， 2021，59（02）：83-86.

[5]張忱，張人會，李仁年，孟凡瑞.液環(huán)泵葉片軸向葉頂凹槽間隙流場及其性能[J/OL].排灌機械工程學(xué)報：1-6[2021-04-14].

[6]郭禹，秦傳新，章守宇.基于結(jié)構(gòu)差異性正方體魚礁的流場效應(yīng)研究[J/OL].漁業(yè)科學(xué)進展：1-11[2021-04-14].

[7]陳磊，譚均軍，石小濤，柯森繁，譚紅林，李鴻.組合隔板式魚道水力特性數(shù)值模擬研究[J/OL].長江科學(xué)院院報：1-7[2021-04-14].

[8]呂兆華.泡沫型多孔介質(zhì)中非達西流動特性的研究[J].工程力學(xué)， 1998（02）：57-64.

中船海洋動力部件有限公司，上海 201306

上海海事大學(xué)，上海 201306