本文以P360 型液環泵噴射器作為研究對象，分別采用結構和非結構計算網格，應用RNG k-ε湍流模型對噴射器的超音速混合氣體流場進行了數值模擬。計算得到了兩種網格條件下噴射器內的壓力、馬赫數及溫度等參數的分布規律。通過與實測結果對比表明，采用結構網格計算得到的液環泵噴射器性能與實測結果較為吻合，也能較好地捕捉到噴射器內激波膨脹波交替循環的流場結構。

　　在液環泵真空系統中，噴射器經常作為液環泵的上游設備，起著提高吸入真空度和增大吸氣量的作用，還可有效改善液環泵腔內低壓造成的汽蝕、振動和噪聲等問題。如圖1 所示，噴射器出口與液環泵入口相接，當液環泵啟動時，噴嘴進氣口與排氣口之間形成壓差；工作氣體流經收縮噴嘴后形成高速射流，吸氣室的壓力低于引射壓力使引射氣體進入吸氣室內；兩股氣流在吸氣室內混合后進入液環泵，經增壓后排出泵外，即完成了整個系統的吸氣、排氣過程。

圖1 液環泵—噴射器工作原理

　　有關噴射器的研究工作目前大多采用理論分析、數值計算的方法。例如于文艷等采用FLUENT 軟件對蒸汽噴射泵器進行了數值模擬，分析了噴射器工作壓力、引射壓力等參數變化對噴射系數的影響。單勇等通過對波瓣噴管引射- 混合器的流場進行了數值計算，得到了流向渦和正交渦產生、衰減的規律。姚亞夫等采用ANSYS 軟件對Laval 噴管內的混合氣流場進行了有限元數值模擬，得到了較直觀的混合氣體流動軌跡。黃思等模擬計算了液環泵噴射器內的氣體流動，探討了結構參數對噴射器內部流場及性能的影響。

　　由于噴射器是復雜的三維結構，上述研究在計算中均采用了非結構的網格，在計算精度方面可能會造成一定的偏差。本文借助ANSYS-CFX 軟件，分別采用結構網格和非結構網格對常用的液環泵噴射器進行數值模擬并通過實測數據的驗證，為液環泵噴射器結構的優化設計提供理論依據。

1、計算模型及網格劃分

　　1.1、計算域和網格生成

　　選取常用的P360 型液環泵噴射器作為研究對象。如圖2 所示，該噴射器計算域由進氣段、噴嘴、引射吸入室、混合室和擴壓室組成。借助ICEM 網格生成軟件，對噴射器計算域分別構造了六面體的結構網格和四面體的非結構網格，網格單元總數均為80 萬左右，如圖3 所示。

圖2 噴射器計算區域

圖3 噴射器計算網格(中心剖面)

　　1.2、計算模型及邊界條件

　　液環泵噴射器內的氣流視為理想氣體的絕熱過程，選取RNG k-ε 湍流模型進行湍流計算。工作氣體入口、引射氣體入口以及噴射管出口邊界條件均采用靜壓值和溫度值，為配合實驗工況，選取以下的數值范圍作為邊界條件：引射壓力P 引=2000～16000Pa(絕對壓力，下同)，出口壓力(背壓)P 出=10500～17000Pa，工作壓力均采用一個標準大氣壓，P 工=101325Pa，環境溫度為28℃。

2、計算結果分析

　　2.1、馬赫數分布

　　根據上述的噴射器工況條件分別進行了模擬計算。作為示例，圖4 和圖5 分別給出了噴射器計算域中心截面的馬赫數分布云圖和軸心線的馬赫數分布曲線(P引=6000Pa，P 出=12400Pa)。由圖可見，使用兩種網格計算得到的噴射器射流核心區流速均超過音速，即M>1。由以往的研究結果得知，當超音速噴嘴出口壓力低于環境壓力時，會在噴嘴出口產生斜激波，經過混合層反射再轉為膨脹波，如此交替循環形成網格模式直至衰減。

　　射流核心區的壓力、流速、溫度等參數也相應地產生周期性波動。采用結構網格計算得到的馬赫數等值線云圖(圖4a)在射流核心區可以看出有明顯的網格模式，軸心線的馬赫數曲線沿著流動方向也有明顯的周期性波動(圖5)。而采用非結構網格計算得到的馬赫數分布則看不出有類似的網格模式(圖4b)和波動規律(圖5)，說明采用結構網格計算能較好地捕捉到噴射器內的激波膨脹波交替循環的流場結構。

 

圖4 噴射器中心截面的馬赫數分布云圖

(P 引=6000Pa，P 出=12400Pa)

圖5 馬赫數沿噴射器軸心線的分布曲線

(P 引=6000Pa，P 出=12400Pa)

　　2.2、壓力分布

　　圖6 和圖7 分別給出噴射器內中心截面的壓力分布云圖和軸心線上的壓力分布曲線(P 引=6000Pa，P 出=12400Pa)。由圖可見，使用結構網格計算得到的壓力分布在射流核心區以引射壓力值為基準上下波動(圖7)；氣體進入混合室后壓力值才穩定持續上升，最終達到噴射器出口壓力。使用非結構網格計算則捕捉不到壓力分布在射流核心區類似的波動規律。

圖6 噴射器中心截面的壓力分布云圖

(P 引=6000Pa，P 出=12400Pa)

圖7 壓力沿噴射器軸心線的分布曲線

(P 引=6000Pa，P 出=12400Pa)

　　2.3、溫度分布

　　圖8 和圖9 分別給出噴射器內中心截面的溫度分布云圖和軸心線上的溫度分布曲線(P 引=6000Pa，P 出=12400Pa)。由圖可見，使用結構網格計算得到的射流核心區的溫度分布在一段范圍內處于不穩定的低溫狀態(圖8a，圖9)。氣體進入擴壓室后溫度值迅速上升。而采用非結構網格計算得到的射流核心區僅有很小的低溫區域，隨后氣溫沿著流動方向平緩回升(圖8b，圖9)。

圖8 噴射器中心截面的溫度分布云圖

(P 引=6000Pa，P 出=12400Pa)

圖9 溫度沿噴射器軸心線的分布曲線

(P 引=6000 Pa，P 出=12400 Pa)

3、實驗驗證

　　液環泵噴射器的性能測試在廣東省佛山水泵廠有限公司的水泵測試站完成。該測試站采用多功能參數測量儀的自動測試系統，整套系統達到國家B 級精度水平。液環泵真空系統由P360型噴射器與2BE13563 型液環泵組合，液環泵轉速：590r/min，供水量：26m3/h，室溫：28℃。實驗方法是通過閥門調節不同的引射壓力和液環泵入口壓力，測出相應的吸入氣量值。圖10 是采用兩種網格計算的吸入氣量與實測值的對比結果。

圖10 模擬計算與實測結果的對比

　　由圖可見，采用結構網格計算得到的吸入氣量與實測結果吻合較好。

4、結論

　　本文借助ANSYS-CFX 流動軟件，分別采用結構網格和非結構網格對常用的液環泵噴射器氣體流場進行了數值模擬，通過與實測結果的對比驗證，得到以下結論。

　　(1)設計工況下液環泵噴射器內的氣流屬超音速流動，在射流核心區將產生斜激波和膨脹波交替循環的網格模式；射流核心區的壓力、馬赫數、溫度等參數沿流動方向相應地呈周期性震蕩并逐漸衰減。

　　(2)兩種網格計算得到的噴射器氣體流場對比結果表明，采用結構網格能較準確地模擬計算液環泵噴射器內流場和外部性能，能較好地捕捉到噴射器內激波膨脹波交替循環的流場結構。