不同通風形式下的CO2氣體泄漏擴散數值模擬

王文芳

（安徽理工大學，安徽淮南 232001）

1 概述

CO2是一種無色無味的氣體，密度大于空氣，按照氣體密度劃分原則屬于重氣類氣體，當CO2泄漏后，會與空氣混合形成重氣云，由于重力沉降作用會在地面大量積聚，隨后會因濃度差向周圍空氣中擴散，產生濃度梯度現象。高濃度的CO2氣體對人的生命安全會產生十分隱蔽的危害，雖然這種氣體本身無毒性，但吸入過量會造成人體CO2中毒，原因是高濃度的CO2氣體環境常伴隨O2含量的降低，形成缺氧環境。當人突然進入高濃度CO2氣體環境中，會出現頭暈、心悸、呼吸困難的癥狀，還有可能直接導致人暈厥甚至死亡。因此研究CO2泄漏擴散規律意義重大，不僅可以為CCUS（碳捕捉、利用與封存）技術中提供數據分析與理論支持，保證CO2的管道輸運安全，更重要的是能夠保障輸運人員的生命安全。

在CO2泄漏擴散的規律研究過程中，主要有兩種方法；

1）通過搭建實驗平臺進行泄漏實驗；

2）通過仿真軟件進行泄漏實驗。

相對于搭建實驗平臺的研究方法而言，通過仿真軟件進行模擬實驗的方法成本低，可重復性高，但對于物理模型與計算方法的選擇需要十分準確，一旦出現錯誤，預測的結果就會與實際相差巨大。目前FLUENT 軟件在氣體泄漏擴散模型中應用已十分廣泛，FLUENT 軟件涵蓋各種物理建模功能，可對工業應用中的流動、湍流、熱交換和各類反應進行建模，可模擬從可壓縮流動到不可壓縮流動等許多流體，很多學者采用FLUENT 軟件進行氣體擴散實驗，預測結果與真實數據十分吻合，因此本文采用FLUENT 軟件來對不同通風形式下的CO2擴散規律進行研究，再結合實驗進行物理模型的有效性驗證。

2 ICEM前處理

ICEM 是FLUENT 的網格劃分軟件，作為一款專業的CFD 前處理軟件，主要用于創建幾何模型和進行網格劃分。ICEM 軟件包括兩大類網格劃分方法，結構化網格劃分與非結構化網格劃分方法，結構化網格質量高，但需要手動分塊劃分，需大量人工操作，不適用于復雜模型。

1）非結構化網格可以自動生成，操作簡單。在網格質量要求不高的模型中可以使用。幾何建模的第一步是通過創建點、線、面，在ICEM 中建立一個大小為3m（x）*1.8m（y）*2.5m（z）的三維長方體幾何模型，在x=0m 平面上分割出一個泄漏孔，泄漏孔為一圓形孔，選?。?，0.88，0.865）為圓心，直徑大小為1mm。然后再在y=2.5m 平面上分割出一個自然通風口，自然通風口為一矩形口，大小為0.7m*1.0m，泄漏口命名為inlet，自然通風口命名為outlet。

2）進行網格劃分，首先將總體網格尺寸設置為0.1m，并在氣體流動速度大的區域進行網格加密，局部加密網格尺寸大小設置為0.001m，采用非結構化網格自動生成方法，網格生成方式為四面體為主的混合網格。然后進行網格光順處理，確保網格精度達到要求。

3）選擇FLUENT 求解器輸出網格即可。CO2幾何模型如下圖1所示，氣源從左側泄漏孔（inlet）進入，從右側通風口（outlet）流出。

圖1 幾何模型

3 FLUENT建模

3.1 材料參數設置及邊界條件

在ICEM 中生成的mesh 文件可直接在FLUENT 中打開，通常第一步是要進行網格導入與質量檢查，檢查導入模型的單位是否有誤，塊是否丟失等，然后將求解器設置為壓力基求解器，再設置瞬態模型，由于CO2為重氣，因此需要在z軸方向設置重力加速度。幾何模型內部氣體設置為空氣與CO2的混合氣體，混合氣體密度類型為不可壓縮理想氣體，由于該擴散過程本身不存在化學變化，因此只需要打開組分輸運方程即可。將inlet 類型設置為Mass Flow Rate，初始流量為0.000 5kg/s，inlet 處CO2所占質量分數設置為1。ouelet 類型設置為pressure outlet，初始表壓設置為0，在outlet 處CO2所占質量分數設置為0.000 5，壁面類型設置為標準壁面，近壁面無滑移，操作環境設置為標準大氣壓環境。湍流模型采用雙方程k-ε模型中的標準k-ε方程，打開全浮力效果選項。

3.2 計算求解及初始化設置

采用SIMPLE 算法進行求解計算，其余選項設置為默認，初始化位置從全局開始，在模型內部進行局部初始化，設置CO2質量分數為0.000 5，步長設置為1s，每一步迭代10 次，時間步數設置為70，點擊caculate 便可進行計算。通過改變outlet 類型，即可得到兩種不同通風形式的幾何模型，通過設置四個濃度監測點可觀察CO2氣體在兩種不同通風狀態下的空間濃度變化情況，四個監測點位置及個數如表1所示。

表1 四個監測點坐標

4 實驗與仿真

4.1 不同通風形式下的監測點濃度對比

為簡化模擬過程，作出如下幾點假設：

（1）泄漏氣體不發生相變；

（2）泄漏氣體與環境無能量交換，無溫度變化；

（3）泄漏氣體不發生化學反應。

如圖2，為存在自然通風窗口的空間與封閉空間同一位置的四個監測點濃度變化對比，模擬時長為670s?？梢钥闯?，泄漏時間小于200s 時，四個監測點的濃度變化是完全一致的，200s 之后，自然通風口開始發揮作用，大量的氣體開始從自然通風口向外擴散，空間內部的濃度場分布開始發生變化，相比較封閉空間，帶有自然通風窗口的空間內部濃度全部明顯下降，上方兩監測點（圖2b 和圖2d）比下方兩監測點（圖2a 和圖2c）出現虛實線濃度差的時間更早，這是因為自然通風口的位置設置較高所導致，由此可見重氣的典型重力沉降作用比浮力擴散作用在擴散前期更加明顯，當無外界風干擾，在泄漏前期的濃度場分布是以重力為主要支配動力。

圖2 四個監測點濃度對比

4.2 仿真可行性驗證

由于仿真所得預測數據的真實性最終要靠實驗來驗證，因此，參考仿真模型參數進行了CO2氣體泄漏實驗，泄漏實驗系統由與幾何模型同尺寸的實驗艙、CO2氣源及管線、數據采集及存儲系統等設備組成。實驗系統結構如圖3所示；實驗通過鋼氣瓶釋放氣源，緩沖罐進行穩壓，流量計進行流量設置，再經泄漏管道，從泄放口泄漏至實驗艙。實驗艙中布置有四個型號為COZIR-WX-20%的CO2濃度傳感器，幾何位置與仿真監測點位置保持一致。利用組態王開發監控平臺，通過連接數據庫即可實現CO2濃度實時監控。

圖3 實驗系統結構圖

圖4為帶有自然通風口的實驗與仿真的四個監測點濃度曲線圖，由仿真與實驗結果對比可知，四個監測點仿真數據與傳感器所測濃度走向大致吻合，因此，在數值模擬中所采用的物理模型符合要求，可以充分驗證該物理模型的正確性。

圖4 四個監測點濃度對比

4.3 泄漏擴散云圖分析

圖5a 和圖5c 為帶有自然通風窗口的不同時刻y=0.88m 平面濃度云圖；圖5b 和圖5d 為封閉空間的不同時刻y=0.88m 平面濃度云圖。由圖5a 和圖5c 可以看出，CO2氣體泄漏后首先向平面下方積聚，隨著艙內濃度的升高，CO2氣體向平面上方不斷擴散。整個泄漏過程中濃度在該平面上的分布是不均勻的。在180s時，圖5a和圖5b中均可觀察到泄漏口上方氣體已產生明顯的濃度梯度，但由于泄漏速度較大，泄漏出來的氣體射程較遠，導致大量氣體積聚在實驗艙右側，平面內右側濃度明顯高于左側。在第360s時，由圖5d可以看出封閉空間仍是右側濃度高于左側濃度，而由圖5c可見360s時自然通風口處濃度較平面其他位置明顯更低，此時空間內右側濃度已明顯低于左側濃度，說明該位置附近氣體已從通風口大量逸出，在整個擴散過程中，平面上的濃度最高點始終集中在泄漏口附近。

圖5 不同時刻y=0.88m平面上的濃度分布云圖對比

5 結論

1）經實驗證明，本文所建立的仿真模型是有效的。

2）當氣體發生泄漏時，泄漏初始時帶有通風窗口的空間與封閉空間監測點濃度變化十分接近，超過200s 后監測點濃度開始出現差別，帶有自然通風口的空間在實驗后期整體濃度更低，但總體濃度仍然保持上升趨勢。

3）當氣體泄漏速度較大時，受限空間內的濃度分布是不均勻的，濃度分層現象只出現在局部位置，不同通風形式左右兩側濃度分布隨時間變化規律也是不同的。