二氧化碳致裂器相變流場有限元分析

徐坤鴻 胡軍科 牛奇斌 趙存

摘 要：本文主要研究二氧化碳致裂器相變流場問題。在現有致裂器結構的基礎上，借鑒前人對液化氣體儲罐的研究，利用Fluent流體仿真軟件，選用標準k-ε湍流模型和VOF多相流模型，通過用戶自定義函數來定義相變，建立有限元仿真模型。通過仿真得到致裂器內壓力與氣相體積分數的時間響應曲線，發現升壓速率隨時間逐漸升高，相變速率隨時間逐漸降低。

關鍵詞：二氧化碳致裂器；相變；Fluent

中圖分類號：TP391.9 文獻標識碼：A 文章編號：1003-5168（2018）08-0139-04

Finite Element Analysis of the Phase Transition

in Carbon Dioxide Fracturing Device Flow Field

XU Kunhong1 HU Junke1，2 NIU Qibin1 ZHAO Cun1

（1.College of Mechanical & Electrical Engineering， Central South University，Changsha Hunan 410083；2.State Key Laboratory of High Performance Complex Manufacturing， Central South University， Changsha Hunan 410083）

Abstract： This paper mainly studied the phase change flow field of CO2 cracking device. On the basis of the existing structure，reference to the research of liquefied gas storage tank， the finite element model was established by adopting standard turbulent model and VOF （volume of fluid） model，and user-defined function was compiled to define phase interaction. Time response curve about pressure and volume fraction of the gas phase was generated. It is found that the rising pressure rate increases with time， and the phase change rate decreases with time.

Keywords： carbon dioxide fracturing device；phase transition；Fluent

近年來，囿于炸藥爆破安全性低、噪聲和污染等問題，一種非炸藥爆破技術——二氧化碳致裂技術取得了極大發展。二氧化碳致裂技術是利用液態的二氧化碳在有限的空間內吸熱汽化后體積劇烈膨脹產生高壓，使煤層、巖石或混凝土破裂的技術[1]。這種技術在爆破過程中無明火產生，安全性相對較高，可廣泛用于煤礦開采、采石場等工程項目[2，3]。

目前，關于致裂器的研究大多集中于結構的改進方面，并沒有關于致裂器內部相變機理的研究。筆者在借鑒前人對噴火環境下液化氣儲罐熱響應行為數值模擬[4]及液化氣罐受熱引爆機理分析[5]，對二氧化碳致裂器相變流場機理進行研究，以期對行業內致裂器的結構改進提供理論指導。

1 仿真模型

本文選用Fluent流場分析軟件進行致裂器相變模擬，選用標準k-ε湍流模型，采用有限體積法離散瞬態控制方程，求解方式采用非耦合隱式方程。為了提高求解的精度，對流項和壓力項均采用二階迎風格式。為加快非穩態過程的求解速度，采用PISO算法。多相流模型選用Fluent自帶的VOF（Volume of Fluid）模型，對相界面進行重構。由于軟件本身并不帶有計算相變的模型，需要用戶自行加入所需的UDF（User-defined Function）程序段。

1.1 幾何模型

二氧化碳致裂器的作用機理是利用高壓氣體做功達到致裂效果。其工作原理是：將液態二氧化碳充入密閉的裝有發熱裝置的高強度金屬筒體內，當激活發熱裝置時，金屬筒體內會迅速形成高壓氣體。當氣體壓力上升達到筒體的強度極限時，筒體會沿軸向刻槽處破裂泄放高壓氣體，作用于巖石孔壁，產生的高壓應力達到致裂破壞的效果。致裂器的結構簡圖如圖1所示。

1.2 材料定義

本模型的初始相為液態CO2，第二相為氣態CO2，其物性參數如表1所示。

定義相間關系蒸發與冷凝，蒸發速率取0.5，臨界溫度304.2K。

1.3 網格模型

利用Workbench前處理軟件Design Modeler對模型的流體域進行抽取，因為模型的對稱性，取其一半，此時的流體域為原來致裂器流體部分的1/4，由于結構的規則性，適合采用結構化網格劃分。采用Fluent的前處理軟件ANSYS ICEM進行網格劃分，采用O型網格劃分，定義全局網格最大尺寸2.5mm，模型網格總數約為22萬，整體網格質量達到0.7以上，如圖2所示，符合計算要求。本文對流體域施加的是體積熱源，在Fluent中按照加熱管的尺寸對流體域進行分割，分離出體積熱源流體域，如圖3所示。

1.4 自定義UDF

由于本文涉及的是液態CO2和氣態CO2之間的相變傳質傳熱，液態CO2受熱發生相變，這里只列出蒸發過程相變方程。

蒸發過程（[TLiq>TS]），液體源相、蒸汽源相、能量源相分別如式：

[SM=-0.1·αLiq·ρLiq·TLiq-TSTS] （1）

[SM=0.1·αLiq·ρLiq·TLiq-TSTS] （2）

[SM=-0.1·αLiq·ρLiq·TLiq-TSTS·ΔHevap] （3）

式中：[αLiq]為液相的體積分數；[ρLiq]為液相的密度，kg/m3； [ΔHevap]為蒸發潛熱，kJ/kg。

1.5 初始條件及邊界條件

致裂器筒體內充滿液態CO2，充裝壓力為6MPa，溫度為295K，CO2的臨界溫度為304.2K，臨界壓力為7.38MPa，處于臨界狀態的密度為463.9kg/m3。根據筒體內液態CO2的量和其汽化潛熱值，確定加熱裝置的能量值為324.36kJ。為了確保計算的收斂性，定義0.000 01s為計算時間步長，每個時間步長下迭代30次，對數值模型進行10 000歩瞬態數值計算，以對致裂器相變過程進行流場仿真。

2 仿真結果分析

圖4為0～100ms下致裂器內壁面的壓力-時間響應曲線和致裂器筒體內氣相的體積分數。

從曲線上可以看出，致裂器內壁面上的壓力由初始狀態下的6MPa逐漸升高，升壓速率緩慢變大。致裂器筒體內氣相的體積分數在0～10ms內與時間基本呈現線性關系，汽化速率很大；在10～40ms內，氣相的體積分數與時間呈現的是非線性關系，汽化速率減小；在40～100ms內，氣相的體積分數與時間呈現的是線性關系。現在從壓力場和相變程度對致裂器筒體內相變過程的流場進行具體分析。

由于致裂器是用于工程爆破的，而致裂器筒體內的相變膨脹增壓最終是要作用于致裂器筒體的內壁使其破裂，所以這里對壁面的壓力場進行分析。圖4中壓力-時間響應曲線表明，致裂器筒體內的升壓速率是在緩慢增大的，這是筒體內相變程度和高溫溫度場傳導擴散耦合作用的結果。

圖5所示為不同時刻致裂器筒體內壁的壓力分布云圖。從圖中可以明顯看出：壓力變化劇烈，從初始的6MPa，10ms時的10.85MPa到100ms的79.42MPa，且4個時刻下壁面壓力分布基本上都是均勻的。

致裂器筒體的壓力能是由分子的不規則運動產生的。從分子運動理論可知，流體的壓強是大量分子頻繁碰撞壁面而產生的，單個分子對壁面的碰撞時間極短，作用是不連續的，但大量分子頻繁碰撞壁面，對物面產生的壓力具有各向同性，對壁面的作用力是持續的、均勻的，因此致裂器筒體內壁面的壓力在任意時刻都是均勻的。

圖6顯示的是致裂器筒體內流場1～100ms對稱面上氣相的體積分數云圖。

初始狀態時，加熱管所在流體域受到施加的體積熱源，所在流體域的溫度瞬間升高，該區域所在的流體瞬間發生汽化。如圖6（a）所示，當t=1ms時，某些流體單元的氣相體積分數達到0.22，而氣相基本上都分布于加熱管所在的流體域，遠離加熱管的液相區域沒有發生相變。

1～10ms時，兩相流層以加熱管流體域為中心急速向外擴張。同時，已經發生相變的加熱管流體域內的相變程度加劇。由于加熱管體積熱源的施加，整個流體域的過熱度持續增大，兩相流層開始迅速雙向擴展，分別擠壓氣相空間和液相空間。結合圖4可以看出，兩相流層的迅速擴張狀態一直持續到10ms才稍稍減緩。

10～40ms時間段內，由于流體域內壓強急劇增大，加之液化氣體汽化消耗的能量增多，因此兩相流層的擴張速度降低。但是，某些單元的氣相體積分數已經達到了1.0，這說明在某些單元內已經完全相變。到40ms時，兩相流層擴展至容器壁面，而遠離加熱管流體域的液相區域仍然沒有發生相變。

40～100ms時間段內，流場形態相對比較穩定，兩相流層保持著均勻速度擴張，致裂器筒體內介質基本處于混亂無序的沸騰狀態，壓力繼續增大，兩相流層基本覆蓋了加熱管等長的流體域，遠離加熱管所在的流體域仍然沒有發生相變。

100ms時，雖然兩相流層里很大一部分單元內的體積分數達到了1.0，但根據圖4得知，整個流體域的氣相體積分數為0.36。可以預見的是，在100ms之后，兩相流層內的相變會繼續進行，但由于壓力的繼續增大，相變速率會繼續降低，整個流體域的氣相體積分數不會得到較大增長，而遠離加熱管所在的流體域仍然也不會發生相變。

3 結論

利用Fluent流體仿真軟件建立了CO2致裂器仿真模型，詳述了建模的方法與過程，從壓力和相變的時間響應曲線、壓力場與相變流場的變化過程3個方面，闡述了CO2致裂器相變的機理與過程，并得出以下結論。

①致裂器筒體內壓分布均勻，隨著時間的推移，致裂器內壓的升壓速率逐漸變大，100ms時的壓力達到了79.4MPa。

②氣相的體積分數反映了筒體內CO2的相變程度，相變速率隨時間逐漸降低，因為高壓給相變的進行增加了難度，100ms時，加熱管等長的流體域都發生相變。

③由云圖可知，壓力與相變是兩個互相制約的因素；相變的程度與加熱管的結構有關。

參考文獻：

[1]張小康，丁亞明.非炸藥爆破破巖展望[J].廣東化工，2017（15）：132-133.

[2]王軍，肖永勝.用二氧化碳爆破技術開采某石灰石礦的大理石材[J].現代礦業，2015（6）：15-17.

[3]王大賽.二氧化碳爆破在掘進工作面的應用[J].化工管理，2016（15）：101-102.

[4]俞昌銘，單彥廣，肖金生，等.液化氣儲罐受熱引爆機理分析[J].北京科技大學學報，2013（04）：522-530.

[5]邢志祥.火災環境下液化氣儲罐熱響應動力過程的研究[D].南京：南京工業大學，2004.