真空背壓下多截面細徑管內氣體流動特性分析

2015-03-21 任旭東 中國工程物理研究院總體工程研究所

  采用數值模擬方法研究了真空背壓下多截面細徑管的氣體流動問題,并結合理論和實驗結果對管內氣體流動特性進行了分析。研究結果表明:在初始時間段,容器壓力下降較快,細徑管內氣體處于流動壅塞狀態,壅塞狀態穩定后,流速下降緩慢,使得壅塞流動狀態可以持續很長時間;壅塞狀態下,容器壓力隨時間呈對數衰減關系,容器壓力越小,壅塞狀態持續時間越短。同時,壅塞狀態下截面突變區域出現了膨脹波和超聲速環。隨著容器壓力的進一步下降,導管入口的流速降低以及導管長度小于臨界管長,不足以使氣流加速到音速,產生壅塞條件不足,導管Ⅱ內的壅塞狀態率先結束,導管Ⅲ和導管Ⅰ內的壅塞現象依次消失。細徑管內氣體轉為層流狀態,容器壓力下降較慢。真空技術網(http://bjjyhsfdc.com/)認為數值模擬結果與實驗結果吻合較好,為多截面細徑管內氣體流動特性分析提供依據。

  在航空航天、國防軍事以及工業領域,細徑管路得到了廣泛的應用。因受空間限制,細徑管路系統中存在多處截面過渡區域,管徑變化多,形狀復雜,流場速度梯度較大,粘性耗散作用明顯。因此,如何提高細徑管路的流動效率是一個難點。氣體通過管道已經被廣泛研究,高壓或大氣環境下,李軍結合氣動充放氣系統的實際情況,提出一種考慮摩擦和傳熱的一維非定常流場方法。Jason通過對絕熱氣體流動和等溫氣體流動分析,研究了等直徑管路的天然氣壅塞現象,提出管道燃氣在兩種流動下最大流量的估算方法。劉慶堂研究了天然氣流動時變徑管道壅塞流動特性分析,引入臨界壓強比,建立了臨界壅塞的判定條件,建立了壅塞流動狀態下的流量計算公式。丁英濤采用實驗研究和數值計算相結合的手段研究了微型噴管內氣體的流動特性。低壓環境下,周曉對不同真空管道阻塞比條件下高速列車車體的空氣阻力進行了數值模擬,得出真空管道內阻塞比對列車空氣阻力的影響規律。劉加利深入系統地研究了真空管道高速列車氣動阻力特性和系統參數設計方法,研究了管道壓力、阻塞比和列車速度對列車氣動阻力的影響。然而通過細徑管從中低壓到真空環境中的流動研究甚少,目前已知僅有廖彬通過試驗的手段對帶細徑導管的容器放氣性能影響因素進行了研究,指出了容器的放氣速率與氣體種類、管道長度和截面尺寸、背壓等因素的關系。

  本文對真空背壓下多截面細徑管的氣體流動問題進行了數值模擬,并結合理論方法和實驗結果對多截面細徑管的氣體流動特性進行了分析。

1、數值模型與計算方法

  放氣過程中,放氣容器內充入0.6MPa氣體,通過不同幾何尺寸的細徑管排氣到真空室里。實際應用中,不同尺寸的細徑管通過閥門連接,使得流動情況比較復雜,進行完全數值模擬比較困難,因此對流動模型進行了簡化,模型中包括容器、三段不同內徑不同長度的細徑管、3m×3m 的模擬真空區域,如圖1所示。由于簡化后的模型結構以及氣體流動是軸對稱的,所以將數值模擬簡化為二維計算。

  數值模擬采用商業軟件FLUENT,基于有限體積法對圖1的模型進行數值模擬。計算域與網格劃分結果如圖2所示。計算初始條件為:容器內初始壓力為0.6MPa,細徑管與模擬真空區域的初始壓力為300Pa。計算湍流模型采用標準的k-ε 湍流模型并封閉方程組,考慮到細徑管內流場速度梯度較大,溫度和動量方程強烈耦合,對溫度和密度所用的亞松弛因子進行了亞松弛,壁面邊界均采用絕熱條件和無滑移邊界條件。

數值模型示意圖

圖1 數值模型示意圖

網格劃分示意圖

圖2 網格劃分示意圖

3、結論

  本文對真空背壓下多截面細徑管的氣體流動問題進行了數值模擬,并結合理論方法和實驗結果對多截面細徑管的氣體流動特性進行了分析。主要有以下結論:

  (1)在初始時間段,容器壓力下降較快,細徑管內氣體處于流動壅塞狀態,壅塞穩定后,細徑管內的流速下降緩慢,使得壅塞流動狀態可以持續很長時間;壅塞狀態下,容器壓力隨時間呈對數衰減關系,容器壓力越小,壅塞狀態持續時間越短,此外,壅塞狀態下截面突變區域出現了膨脹波和超聲速環。

  (2)導管入口的流速降低以及導管長度小于臨界管長,不足以使氣流加速到音速,產生壅塞條件不足,從而壅塞現象消失。

  (3)導管Ⅱ內的壅塞狀態率先結束,導管Ⅲ和導管Ⅰ內的壅塞現象依次消失。這與導管的幾何尺寸有關。

  (4)隨著容器壓力的進一步下降,細徑管內氣體處于層流狀態,容器壓力下降較慢,主要是因為氣流的流速受導管截面積和長度的限制而顯著降低。

  (5)數值模擬結果與實驗結果吻合較好,為多截面細徑管內氣體流動特性分析提供依據。