多孔陶瓷催化劑載體的流動特性

張任平，顏柳菁

(景德鎮陶瓷學院 江西 景德鎮 333001)

多孔陶瓷催化劑載體的流動特性

張任平，顏柳菁

(景德鎮陶瓷學院 江西 景德鎮 333001)

建立了氣體在多孔陶瓷中流動的三維模型，模型中考慮氣體在多孔陶瓷中流動的粘性損失項和慣性損失項。數值模擬計算所需的粘性和慣性阻力系數通過實驗獲得。研究表明氣體在多孔陶瓷中流動壓力線性下降，速度在入口處急劇下降，隨之下降較為緩慢，在不同截面處的分布趨勢基本相似。

多孔陶瓷；流動；催化劑載體；尾氣凈化

0　引　言

多孔陶瓷是一種經高溫燒成、體內具有大量彼此相通并與材料表面也相貫通的孔道結構的陶瓷材料。根據成孔方法和孔隙結構，多孔陶瓷可分為三類：(1)粒狀陶瓷； (2)泡沫陶瓷；(3)蜂窩陶瓷。泡沫陶瓷的氣孔率一般為80-90%，蜂窩陶瓷的氣孔率一般為70%左右，粒狀陶瓷結體的氣孔率一般為30-50%。多孔陶瓷在金屬熔體過濾凈化技術、精過濾技術、催化劑載體、敏感元件、隔膜材料和降低噪聲等領域有廣泛應用。

國內外學者對于多孔陶瓷的研究已經做了大量的工作，如David等[1]準靜態離散單元法(DEM)用來模擬部分燒結多孔陶瓷的彈性性能；夏建國等[2]綜述了多孔陶瓷應用在汽車尾氣催化凈化器載體的研究進展，著重介紹了近年來為提高凈化效率而新開發的多孔陶瓷載體材料；王一鳴等[3]提出了一種新型的多孔太陽墻采暖房, 其核心部分多孔太陽墻是由多孔陶瓷構成， 在晴朗的冬季, 對多孔太陽墻采暖房的采暖性能進行了實驗研究。研究重點檢測了不同外界環境(主要是指太陽輻射照度和環境溫度) 下, 多孔太陽墻采暖房的采暖及儲熱性能；張志金等[4]以中間相瀝青添加質量分數為50%的Si粉制備的炭泡沫預制體為坯體，在高溫感應燒結爐中結合反應燒結工藝制備了SiC多孔陶瓷預制體。利用擠壓鑄造工藝制備了SiC多孔陶瓷增強鋁基復合材料。張志金等[5]以中間相瀝青添加55%(質量分數,下同)的Si粉混合物為原料，制備了含Si的炭泡沫模板；蔣兵等[6]綜述了多孔陶瓷孔徑及孔徑分布的常見測定方法，比較了各種測試方法的優缺點,認為電子顯微鏡圖像分析法是最直接有效的測定方法，并對多孔陶瓷的測試表征方法提出了展望。鮑遠通等[7]實驗研究多孔泡沫陶瓷的流動阻力特性,測定了不同空氣溫度和流速下泡沫陶瓷的流動阻力特性,獲得了泡沫陶瓷流動阻力與氣流溫度、流速及孔隙結構之間的關系。白鳳武等[8]建立了碳化硅泡沫陶瓷吸熱體在聚光輻射能流加熱條件下空氣與吸熱體骨架間的傳熱過程的數學物理模型，分析了吸熱體厚度、吸熱體孔徑大小對傳熱性能的影響，為優化設計該類吸熱器提供了理論依據；李德波等[9]對甲烷氣體在泡沫陶瓷內燃燒進行了數值模擬和實驗研究；劉建建等[10]對采用新型SiC泡沫陶瓷填料的濕化器在加壓條件下的濕化性能進行了實驗研究,分析了水氣比、進口水溫、操作壓力以及進口空氣溫度對濕化過程的影響；朱華清等[11]分析了泡沫陶瓷材料的吸聲性能缺陷及原因,研究有機修飾對泡沫陶瓷吸聲性能的影響。從以上的研究文獻可以看出，研究主要集中在多孔陶瓷的制備方法和物理性能，對于氣體在多孔陶瓷的流動特性也只有實驗研究的報道[7]，而有關氣體在多孔陶瓷的流動特性的數值模擬研究還非常缺乏。為此，本文以應用于汽車尾氣處理的催化劑載體多孔陶瓷為研究對象，基于計算流體動力學的基本理論數值模擬了氮氣在裝有多孔陶瓷的汽車尾氣處理裝置中的流動特性。

1　物理模型及網格劃分

以多孔陶瓷為催化劑載體的汽車尾氣處理裝置如圖1所示，由殼體、載體和催化劑三部分組成。在其中間部分裝有多孔陶瓷，在多孔陶瓷中附有催化劑， 因多孔陶瓷表面的催化劑對汽車發動機排出的廢氣中的污染物CO、HC和NOx等有害氣體進行氧化和還原，將其轉變為無害的CO2、 H2O和N2氣體排入環境中。多孔陶瓷區的長度為200 mm，寬度為70．7 mm，高度為101．2 mm。多孔陶瓷的氣孔率為38．6%，平均孔徑為45 μm，滲透率0．3679 μm2。計算模型的網格劃分如圖2所示。

2　理論模型

2.1控制方程

氮氣在多孔陶瓷流動受基本守恒定理的支配，過程滿足質量、動量守恒。計算湍流時采用k-εr湍流計算模型，相關的控制方程可以表示如下。

質量守恒方程(又稱連續性方程)：

該方程是質量守恒方程的一般形式，它適用于可壓流動和不可壓流動。Sm是自定義源項。

圖1　裝有多孔陶瓷的汽車尾氣處理裝置Fig.1 Automobile exhaust gas purifier with porous ceramic catalyst carrier

圖2　計算模型的網格示意圖Fig.2 The mesh diagram of the calculation model

在慣性(非加速)坐標系中i方向上的動量守恒方程為：

其中，p是靜壓，τij是應力張量，ρgi和Fi分別為i方向上的重力體積力和外部體積力。

多孔介質的動量方程具有附加的動量源項。源項由兩部分組成，一部分是粘性損失項(Darcy)；另一個是內部損失項：

其中，Si是i向(x, y, 或z向)動量源項，Dij和Cij是規定的矩陣，v 是速度大小。在多孔介質單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度方陣）成比例。

對于簡單的均勻多孔介質：

式中，α是滲透性系數，C2是慣性阻力系數，也就是將Dij和Cij矩陣均簡化為對角矩陣，他們對角上的元素分別為1/α和C2，其它元素都是0。

2.2慣性和粘性阻力系數的計算過程

(1)通過實驗測定不同進口速度時通過多孔介質的壓力降；

(2)多孔介質的厚度為△n，根據多孔介質內速度的壓降實驗數據，通過擬合得到通過多孔介質的壓降和進口速度的關系式；

(3)動量方程的簡化形式就是將壓降與源項進行相關，因此可得：

(4)通過式(5)和擬合關系式進行比較，可以計算得到粘性和慣性阻力系數。

2.3邊界條件和數值計算

裝有多孔陶瓷的汽車尾氣處理裝置進口速度為22．6m/s，紊流強度為10%，水力直徑為28．284 mm。出口為壓力出口邊界條件，表壓為0，回流紊流強度為5%，水力直徑為33．137 mm。中間部分為多孔陶瓷區，粘性和慣性阻力系數通過實驗數據計算而得。為減小擴散誤差, 采用二階迎風格式離散控制方程。基于控制容積有限差分法和求解壓力耦合方程的半隱式(SIMPLE)算法求解控制方程式，當相關變量殘差滿足要求，并且進出口質量守恒，即認為計算結果收斂。

3　結果分析與討論

圖3給出了多孔陶瓷中不同截面的速度分布。如圖所示，當氣體流過多孔陶瓷時，中心的速度較高，截面速度分布形成一個菱形的形狀，選取的三個截面的速度分布趨勢非常相似，可見氣體在多孔陶瓷中流動時速度變化不是特別大。

圖4給出了多孔陶瓷中不同截面的壓力分布。如圖所示，壓力在整個長度方向變化特別明顯，氣體流過多孔陶瓷后，壓力明顯下降，在長度方向的不同位置處，壓力分布趨勢也相差非常大，在多孔陶瓷的進口和中間部分，截面中間部分的壓力呈菱形分布，四周壓力較低；在多孔陶瓷的出口位置，截面的壓力呈條紋形分布。

圖3　多孔陶瓷中不同截面的速度分布Fig.3 The gas velocity distribution in different cross-sections in porous ceramic

圖4　多孔陶瓷中不同截面的壓力分布Fig.4　The gas flow pressure distribution in different cross-sections of porous ceramic

圖5 給出了汽車尾氣處理裝置截面y=35 mm位置速度矢量分布圖。如圖所示，進口速度最大，當氣體進入到多孔陶瓷段時，速度分布比較平均，整個過程速度在緩慢減小。圖6給出了氣體流過過孔陶瓷時中心線的速度分布。如圖所示，當氣體流過多孔陶瓷的起始段，速度急劇下降，隨之速度下降非常緩慢，直到氣體流出多孔陶瓷時，速度有些許回升。氣體流過多孔陶瓷時中心線的壓力變化如圖7所示，壓力變化呈線性分布，從入口到出口壓力線性下降。

圖5　截面y=35 mm速度矢量分布圖Fig.5 The velocity vector distribution at y=35mm in the cross-section

圖6　多孔陶瓷中心線速度軸向分布Fig.6 The axial distribution of centerline velocity for porous ceramic

圖7　多孔陶瓷中心線壓力軸向分布Fig.7 The axial distribution of centerline pressure for porous ceramic

4　結　論

本文以多孔陶瓷應用于汽車尾氣處理的催化劑載體，利用計算流體動力學基本理論數值研究了氮氣在裝有多孔陶瓷的汽車尾氣處理裝置的流動特性。研究表明：

(1)當氣體流過多孔陶瓷時，中心的速度較高，截面速度分布形成一個菱形的形狀，選取的三個截面的速度分布趨勢非常相似；氣體流過多孔陶瓷的起始段，速度急劇下降，隨之速度下降非常緩慢，直到氣體流出多孔陶瓷時，速度有些許回升；

(2) 壓力在整個長度方向變化特別明顯，氣體流過多孔陶瓷后，壓力明顯下降，在長度方向的不同位置處，壓力分布趨勢也相差非常大，在多孔陶瓷的進口和中間部分，截面中間部分的壓力呈菱形分布，四周壓力較低；在多孔陶瓷的出口位置，截面的壓力呈條紋形分布。

[1] JAUFFRES D, MARTIN C L, LICHTNER A, et al． Simulation of the elastic properties of porous ceramics with realistic microstructure [J]． Modelling and Simulation in Materials Science and Engineering, 2012, 20(4): 450-459．

[2] 夏建國, 王少洪, 侯朝霞, 等． 多孔陶瓷在汽車尾氣處理中的研究進展[J]． 兵器材料科學與工程, 2012, 35(4): 93-97．

[3] 王一鳴, 楊昆, 劉偉, 等． 一種新型的多孔太陽墻-采暖房[C]．//2011年中國工程熱物理學會傳熱傳質學學術會議論文集, 2011:1-6．

[4] 張志金, 王揚衛, 于曉東, 等． 三維網絡SiC多孔陶瓷增強鋁基復合材料的制備[J]． 稀有金屬材料與工程, 2009, 38(2): 499-501．

[5] 張志金, 于曉東, 王揚衛, 等． 三維網絡碳化硅多孔陶瓷的制備[J]． 航空材料學報, 2012, 32(4): 57-61． DOI:10．3969/ j．issn．1005-5053．2012．4．011．

[6] 蔣兵, 翟涵, 李正民, 等． 多孔陶瓷孔徑及其分布測定方法研究進展[J]． 硅酸鹽通報, 2012, 31(2): 311-315, 321．

[7] 鮑遠通, 劉聯勝, 王華, 等． 多孔泡沫陶瓷中的流動阻力研究[J]． 河北工業大學學報, 2012, 41(2): 60-64．

[8] 白鳳武, 王艷． 碳化硅泡沫陶瓷空氣吸熱器性能數值模擬[C]．//2011年中國工程熱物理學會傳熱傳質學學術會議論文集, 2011: 1-11．

[9] 李德波, 樊建人, 宋正昶, 等． 甲烷氣體在泡沫陶瓷中預混燃燒的數值模擬和實驗研究[C]． //中國工程熱物理學會2008多項流學術會議論文集, 2008．

[10] 劉建建, 徐震, 肖云漢, 等． 泡沫陶瓷填料濕化器加濕性能實驗研究[J]． 熱能動力工程, 2011, 26(1): 27-30．

[11] 朱華清, 梁華銀． 有機修飾對泡沫陶瓷吸聲性能影響的研究[J]． 陶瓷學報, 2010, 31(4): 623-626．

Flow Characteristics of Porous Ceramic Catalyst Carrier

ZHANG RenpingYAN Liuqing
(Jingdezhen Ceramic Institute, Jingdezhen 333001, Jiangxi, China)

A three-dimensional model is developed for predicting the gas flow in porous ceramic, where the viscous and inertial losses are accounted for. Viscous and inertial resistance coefficients are calculated through the experimental data, and the resultant data are used for the model calculation. The study indicates that the gas flow pressure in porous ceramic decreases linearly, the gas velocity demonstrates a sharp decline at the entrance, then declines more slowly; the gas velocity distribution shows a similar trend in different cross-sections.

porous ceramics; flow; catalyst carrier; exhaust gas purification

date: 2015-03-15.Revised date: 2015-03-18.

TQ174.75

A

1006-2874(2015)03-0017-05

10.13958/j.cnki.ztcg.2015.03.004

2015-03-15。

2015-03-18。

通信聯系人：張任平，男，博士。

Correspondent author:ZHANG Renping, male, Doctor.

E-mail:rpzhang@126.com