呂 洪，蔡書娟，許 瀟，劉 湃，陳 專，馬建新

(同濟大學(xué)汽車學(xué)院，新能源汽車工程中心，上海 201804)

前言

隨著汽車工業(yè)的發(fā)展，能源危機和環(huán)境污染日益加劇。為解決這些問題，必須開發(fā)替代型新能源汽車。在各類新能源汽車中，燃料電池汽車(FCV)以工作溫度低、能量轉(zhuǎn)化效率高和環(huán)境友好等突出的優(yōu)越性，成為各國政府和各大公司的研究熱點［1］。

以質(zhì)子交換膜燃料電池(PEMFC)作為發(fā)動機的FCV都采用空氣作為氧化劑，因此城市空氣質(zhì)量對燃料電池發(fā)動機(FCE)的效率和壽命產(chǎn)生重要的影響。解決空氣污染物對燃料電池的影響變得更加緊迫。空氣中的有害氣體(SO2)對PEMFC陰極催化層具有毒化作用，影響其性能和壽命。研究表明，當(dāng)電池以500mA/cm2的電流密度放電，而空氣中SO2濃度達0.05×10-6(3.2mg/m3)時，電池性能會受到不可忽略的影響［2］。

目前通常使用兩種方法來消除空氣中有害氣體對燃料電池性能的影響。一種是開發(fā)新的催化劑，增強催化劑的抗毒性能;另一種是在空氣進入燃料電池之前，將空氣中的有害氣體過濾去除。然而空氣中有害氣體的種類繁多，開發(fā)一種可以抵抗多種有害氣體的新型陰極催化劑仍有較多問題需要攻克。而使用空氣過濾器，在空氣進入電池陰極之前將有害氣體從空氣中過濾去除是一種簡單而行之有效的辦法。

目前，國內(nèi)外對車用燃料電池發(fā)動機空氣過濾器的研究較少。文獻［3］中詳細闡述了設(shè)計燃料電池空氣過濾器結(jié)構(gòu)的方法，認為設(shè)計燃料電池空氣過濾器主要應(yīng)該考慮進氣屬性、濾芯體積、濾芯材料的吸附量和過濾器的壓降等因素。

對于酸性化學(xué)污染物，活性炭是一種比較好的吸附材料。文獻［4］中通過將活性炭浸漬于不同濃度的KOH，對活性炭進行化學(xué)改性，結(jié)果表明:改性劑增加了吸附劑的堿性，顯著提高了活性炭對SO2等氣體的吸附性能。

燃料電池空氣過濾器濾芯要求具有較小的壓降。文獻［5］中根據(jù)壓降和速度的試驗數(shù)據(jù)計算得出黏性阻力系數(shù)1/α和慣性阻力系數(shù)C2，并將其代入FLUENT軟件多孔介質(zhì)模型中，得出空氣通過過濾器的壓降特性，且實驗結(jié)果和FLUENT仿真結(jié)果基本一致，為燃料電池空氣過濾器結(jié)構(gòu)設(shè)計提供了重要的方法和依據(jù)。

文獻［6］中利用FLUENT和MATLAB軟件對多孔介質(zhì)下丙烷催化燃燒反應(yīng)進行了模擬，研究了形狀分別為球狀和柱狀的丙烷催化劑的不同粒徑和氣體溫度等因素對熱量交換、組分質(zhì)量輸運的影響。結(jié)果表明，F(xiàn)LUENT與MATLAB軟件計算結(jié)果比較吻合，對固定床反應(yīng)器的設(shè)計具有重要的指導(dǎo)意義。

本文中應(yīng)用FLUENT軟件對FCE空氣過濾器濾芯的化學(xué)吸附特性進行仿真研究，并通過實驗進行驗證，指導(dǎo)設(shè)計開發(fā)FCE空氣過濾器，提高其使用壽命。

1 計算模型

1.1 模型描述

主要針對FCE空氣過濾器的流體力學(xué)特性、壓力損失和濾芯化學(xué)吸附特性進行仿真。采用Auto-CAD建立的計算幾何模型如圖1所示。過濾器結(jié)構(gòu)包括上下兩個濾芯:上濾芯為紙質(zhì)或無紡布，用于去除空氣中的顆粒污染物;下濾芯為KOH改性的活性炭顆粒，用于去除空氣中的化學(xué)污染物。當(dāng)過濾器的SO2吸附曲線穿透時，即過濾器出口SO2濃度大于0時，定義為過濾器失效。

1.2 模型假設(shè)

(1)反應(yīng)氣體均為理想氣體，且在空氣過濾器中混合均勻;

(2)反應(yīng)氣體為不可壓縮的湍流流動;

(3)濾芯材料為各向同性的均相多孔介質(zhì);

(4)濾芯化學(xué)吸附反應(yīng)為有限速率體積反應(yīng);

(5)網(wǎng)格劃分后每個活性炭小單元的SO2吸附飽和量由實驗和分析得到;

(6)實驗中經(jīng)KOH改性后的活性炭化學(xué)吸附的SO2體積濃度(30%左右)遠大于物理吸附的體積濃度(1%)［7］，因此僅針對濾芯多孔介質(zhì)化學(xué)反應(yīng)區(qū)域進行模擬，不考慮空氣過濾器對SO2隨時間改變的物理吸附。

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

1.3.1 基本方程

由于將氣體視為不可壓縮流體的湍流流動，所以整個過程可以用質(zhì)量守恒方程、雷諾平均N-S方程、湍動能κ和湍動能耗散率ε的輸運方程等基本控制方程來描述。

質(zhì)量守恒方程為

式中:ui為沿i方向的速度分量;fi為沿i方向的質(zhì)量力;p為壓力;ρ為空氣密度;υ為運動黏性系數(shù);υt為渦黏系數(shù)。系數(shù)使用Launder和Spalding的推薦值。

1.3.2 組分方程

式中:左邊第1項為非穩(wěn)態(tài)項，第2項為對流項;右邊第1項為擴散項，第2項為擴散源項。其中，Ck、D、ε、Sk分別表示組分濃度、組分有效擴散系數(shù)、孔隙率和組分源項，下標(biāo)k為組分代碼。

1.3.3 化學(xué)反應(yīng)速率方程

式中:K為反應(yīng)速率常數(shù);K0為指數(shù)前因子(頻率因子);Eα為實驗活化能;Rg為氣體常數(shù);T為溫度。

1.3.4 約束條件方程

通過用戶自定義函數(shù)(UDF)編寫約束條件方程，每個小單元(cell)達到飽和后不再發(fā)生化學(xué)吸附反應(yīng)。

式中:mcSO2為每個小單元吸附SO2的質(zhì)量;qc為濾芯的單位質(zhì)量吸附量;mc為小單元質(zhì)量。

1.3.5 吸附性能-進氣屬性數(shù)學(xué)模型

本文中采用實驗和計算機模擬相結(jié)合的方法計算出平衡吸附量和進氣屬性的關(guān)系。首先通過希洛夫方程擬合實驗數(shù)據(jù)得出吸附穿透時間與床層堆積厚度的關(guān)系，然后利用Wheeler公式得出不同進氣屬性下的穿透時間，最后采用Langmuir吸附等溫公式建立平衡吸附量-空氣流量-初始濃度的三維模型［8］。通過如圖2所示的三維模型可以計算出在一定的SO2初始濃度和空氣流量范圍內(nèi)的平衡吸附量，將其作為編寫UDF約束方程的依據(jù)，用于FLUENT仿真計算。

1.4 模型主要尺寸和邊界條件

為減小運算量，對空氣過濾器進行適當(dāng)簡化，即只針對濾芯多孔介質(zhì)反應(yīng)區(qū)域進行仿真?？諝膺^濾器的邊界條件參考吸附實驗數(shù)據(jù)，如表1所示。

表1 模型邊界條件

2 結(jié)果及討論

2.1 過濾器仿真結(jié)果

圖3和圖4分別給出在過濾器入口處通入初始濃度為10×10-6和50×10-6的 SO2的仿真穿透曲線。從圖中可以看出，在SO2初始濃度為10×10-6時，過濾器的穿透時間為318h;當(dāng)初始濃度為50×10-6時，過濾器的穿透時間為64h。

為驗證該模型預(yù)測過濾器壽命的準(zhǔn)確性，采用實驗方法對該模型進行驗證。因為在低濃度時，過濾器失效時間較長，考慮到實驗的操作性，實驗分為2個階段。

第1階段空氣經(jīng)過過濾器前SO2的濃度為10×10-6，空氣入口速度為0.4m/s，此過程持續(xù)70h，空氣過濾器未失效;第2階段進口空氣的SO2濃度調(diào)整為50×10-6，空氣入口速度仍為0.4m/s，直至過濾器失效。仿真也分兩個階段進行，條件與實驗完全相同。圖5為實驗和仿真所得空氣過濾器濾芯吸附SO2近似的穿透曲線。由圖可見，仿真與實驗結(jié)果的區(qū)別在第2階段:實驗經(jīng)過52h后過濾器失效;而仿真經(jīng)過50h后過濾器失效。

表2給出實驗和仿真所得的穿透時間和穿透平均吸附量。穿透時間上，仿真相對實驗結(jié)果的誤差為-1.6%;而濾芯穿透平均吸附量的誤差為-3.1%。實驗和仿真數(shù)據(jù)基本相符，證明了UDF編寫的約束條件和所建模型與FLUENT仿真的可靠性。

表2 實驗和仿真濾芯穿透時間和穿透平均吸附量的對比

2.2 過濾器結(jié)構(gòu)分析

2.2.1 流速與壓力的分布

圖6和圖7分別為空氣過濾器流速和壓力分布圖?？諝膺^濾器流速最大值約為2.67m/s，楔形的結(jié)構(gòu)設(shè)計有利于氣體流速在顆?；钚蕴慷嗫捉橘|(zhì)主體區(qū)間分布的均勻性?？諝膺^濾器的壓力損失為811Pa，滿足FCE空氣過濾器的要求。

2.2.2 SO2濃度分布

圖8為FLUENT仿真所得空氣過濾器從濾芯開始吸附SO2到吸附飽和過程中不同時刻的SO2濃度分布圖。隨著帶有SO2的氣體進入吸附區(qū)域，濾芯中部吸附的SO2量明顯大于兩側(cè)，120h后出口處檢測出SO2，過濾器失效。顆?；钚蕴繛V芯在達到吸附飽和時，局部區(qū)域出現(xiàn)了“死區(qū)”、“短路”現(xiàn)象，因此活性炭吸附區(qū)域結(jié)構(gòu)需要改進。

2.3 過濾器結(jié)構(gòu)的改進

為提高濾芯的吸附效果，在整體式吸附層外側(cè)設(shè)置導(dǎo)流板，如圖9所示。設(shè)置導(dǎo)流板的作用在于強迫空氣分流經(jīng)過吸附劑后繼續(xù)向前流動。

圖10為結(jié)構(gòu)改進后的空氣過濾器吸附SO2的仿真穿透曲線。仿真條件同前，在過濾器中通入10×10-6的SO270h后，再將 SO2進氣濃度調(diào)整為50×10-6，過濾器經(jīng)過73h后失效(總的吸附時間為143h)。計算得出的濾芯穿透平均吸附量約為 44.84mg/g，比濾芯結(jié)構(gòu)改進前的吸附量提高了40.1%。

圖11為對結(jié)構(gòu)改進后的空氣過濾器進行仿真得到的濾芯接近吸附飽和時的SO2濃度分布圖。由圖可見，設(shè)置導(dǎo)流板使顆粒活性炭濾芯在達到吸附飽和時SO2的濃度分布較為均勻，同時，活性炭的利用率大大提高，“死區(qū)”明顯減小，過濾器使用壽命提高。

3 結(jié)論

通過UDF編寫約束條件方程，應(yīng)用FLUENT軟件，對燃料電池空氣過濾器流體力學(xué)和化學(xué)吸附性能進行仿真研究，并通過實驗驗證;其次，通過建立的仿真模型，對空氣過濾器進行結(jié)構(gòu)優(yōu)化，結(jié)果表明:

(1)仿真結(jié)果與實驗測試結(jié)果符合較好，表明建立的仿真模型可靠，可用于空氣過濾器的結(jié)構(gòu)設(shè)計;

(2)增設(shè)導(dǎo)流板后的空氣過濾器改善吸附效果明顯，提高了使用壽命。

［1］毛宗強.燃料電池［M］.北京:化學(xué)工業(yè)出版社，2005.

［2］楊代軍.大氣污染物對質(zhì)子交換膜燃料電池性能影響的研究［D］.上海:華東理工大學(xué)資源與環(huán)境學(xué)院，2006.

［3］Daniel M Kennedy，Donald R Cahela，Wenhua H Zhu，et al.Fuel Cell Cathode Air Filters:Methodologies for Design and Optimization［J］.Journal of Power Sources，2007，168(1):391-399.

［4］Ma X W，Yang D J，Zhou W，et al.Evaluation of Activated Carbon Adsorbent for Fuel Cell Cathode Air Filtration［J］.Journal of Power Sources，2008，175(1):383-389.

［5］陳專，呂洪，馬曉偉.FCE空氣過濾吸附材料壓降試驗與仿真研究［J］.計算機仿真，2009，26(5):296-298.

［6］Kolaczkowski S T，Chao R，Awdry S.Application of a CFD Code(Fluent)to Formulate Models of Catalytic Gas Phase Reactions in Porous Catalyst Pellets［J］.Chemical Engineering Research and Design，2007，85(11):1539-1552.

［7］馬曉偉.車用燃料電池空氣凈化器的研究［D］.上海:同濟大學(xué)，2007.

［8］劉湃.車用燃料電池空氣過濾器吸附材料化學(xué)吸附行為建模與仿真［D］.上海:同濟大學(xué)，2010.

［9］陳專.燃料電池發(fā)動機空氣過濾器的試驗和仿真研究［D］.上海:同濟大學(xué)，2009.