航空燃油離心泵內流特性研究

昝宇晗，孫晉

（中國航發西安動力控制科技有限公司，陜西 西安 710077）

在航空飛機發動機中，燃油泵屬于關鍵部件，可以為發動機提供燃油或增壓。從類型看，航空燃油泵包括了離心泵、柱塞泵、齒輪泵等。相比之下，航空燃油離心泵的比較優勢相對較大，具體表現在轉速高、尺寸小、揚程高、燃油溫升少、調整流量范圍大、安全可靠性強、質量功率比小、使用壽命長等方面。所以，在我國航空技術快速發展的過程中，對于航空燃油離心泵的使用頻率較高，對于它的研究相對較多。

1 航空燃油離心泵內流動研究方法

航空燃油離心泵內流動屬于三維湍流流動，復雜化程度較高，通常在研究過程中，無論采用理論分析方法，還是實驗研究方法均存在一定的難度。從現階段的對航空燃油離心泵內流動情況的相關研究情況看，在實際研究中為了提高研究效果與精準性，一般會選擇數值模擬方法。例如，在離心泵處于不同工況時的內流壓力場的動靜干涉影響研究、離心泵蝸殼內部渦流量場研究中，均可以采用SA模型、DDES湍流模型、玻爾茲曼方法、雷諾平均方程等進行數值模擬與分析等。本次選擇Fluent軟件結合某型燃油離心泵對不同流量、不同轉速條件下的離心泵內流特性展開分析。

2 航空燃油離心泵內流特性數值模擬分析

2.1 數學模型

2.1.1 幾何模型

航空燃油離心泵幾何模型構成要素包括了出口、入口、葉輪、蝸殼、輸入管道等，其中的參數指標如下：（1）入口直徑：80mm；（2）出口直徑：50mm；（3）葉輪直徑：170mm?；贔luent軟件的離心泵幾何模型，見圖1。

圖1 航空燃油離心泵幾何模型示意

2.1.2 網格生成

目前，在航空燃油離心泵內部非定常流動方面的數值模擬中，通常以動網格技術、滑移動網格技術為主。比較而言，后一種技術在應用時，無需重建網格、計算量相對較少，而且應用中相對穩定，所以，通常研究中一般以后一種技術應用為主進行網格生成。本次研究中按照網格生成需求，采用了滑移動網格技術。具體應用步驟如下：（1）確定獨立計算域，包括葉輪、蝸殼、入口管道區域。（2）生成四面體非結構網絡，主要是將作為研究對象的航空燃油離心泵相關數據輸入后形成具有三維特征的四面體非結構網格。（3）設置滑動網格、確定運動方式，具體是將計算域作為動域并將其設置為滑動網格，選擇繞葉軸軸葉定軸轉動作為確定的轉動方式。

2.1.3 流動控制方程

航空燃油離心泵內部流動屬于非定常不可壓縮黏性流動，主要特征為湍流。因此，在本次研究中選擇的流動控制議程為雷諾平均方向，假定密度為ρ、壓力為P、速度分量為ui、雷諾應力為公式為

2.1.4 湍流模型

在應用雷諾平均議程時，需要應用湍流模型對雷諾應力項進行封閉求解處理，所以，在確定流動控制方程時，本次應用了k?ε湍流模型，其中k表示湍動能，ε表示耗散率。假定動力黏度為u、湍動黏度為ut、平均速度梯度引起的湍動能產生項為Gk、湍動能與耗散率對應的普朗數為kσ=1.0、εσ=1.3，經驗系數C1ε=1.44、C2ε=1.92，那么可以得到輸運方程為

2.1.5 邊界條件

在本次設計的模型中，航空燃油離心泵的邊界條件包括入口邊界（壓力入口，壓力值=0.1MPa）、出口邊界（體積流量出口）、移動壁面（葉片，相對靜止對象為滑移網格區域）、其他壁面（靜止無滑移壁面），將interface設置為計算域交界面。

2.1.6 計算方法

本次研究中選擇的軟件為Fluent，在求解器方面主要以壓力基求解器為準，空間離散應用一階迎風格式，時間離散應用一階隱式格式，迭代方法選擇SIMPLE算法。相關參數指標如下：（1）迭代時間步長選擇1.25×10-4s；（2）流體區域材料：為航空煤油，密度=780kg?m-3；動力黏度 =2.5×10-3kg?（m?s）-3。

2.2 計算結果分析

2.2.1 不同流量下的計算結果分析

首先，將離心泵的轉速設置為每分鐘300r，然后通過對體積流量（出口邊界）的改變，對其工況進行計算，獲得相關物理量分布情況。本次研究中設置的流量共計4個：30m3/h、40m3/h、50m3/h、m3/h。其次，將葉輪作為動參考系，然后對相對速度rν分布（不同物理量分布情況見下圖2）。計算結果分析如下。

圖2 轉速300r/min時不同流量下的葉輪截面相對速度分布示意圖

（1）在不同流量下的葉輪截面相對速度分布情況表明，在流量逐次遞增時，葉輪中流體相對速度rν隨之增加。從原因看，主要是在這種條件下，葉輪輸出體積流量與穿過相對速度穿過葉輪出口時的通量保持一致。從關系方面看，可以表述為：相對速度大小與流量大小呈正相關關系。

（2）不同流量下的蝸殼截面絕對速度分布情況表明，在流量逐次遞增時，葉輪邊緣出流絕對速度ν逐漸遞減。從原因方面看，主要是葉輪邊緣速度合成過程中的絕對速度在右旋轉時的矢量合成結果關系所致。具體而言，確定葉輪邊緣上面的某點后，對葉輪旋轉時的情況進行分析，會發現其中出現了牽連速度矢量eν，假定小流量下的相對速度為νr1、絕對速度為1ν，大流量下的相對速度為νr2、絕對速度為2ν，此時，牽涉速度方向與葉輪周線相切，相對速度方向與葉稍處于近似相切狀態，當轉速相對較高時，牽涉速度大于相對速度，且是后者的2倍到5倍。因此，在這種矢量合成情況下，流量增大、相對速度增加、絕對速度矢量相對減小。從關系方面看，可以表述為：隨著流量增大葉輪邊緣速度相對降低。

（3）不同流量下蝸殼葉輪截面靜壓分布情況表明，從中心到邊緣，泵內靜壓處于逐漸升高狀態，而且在流量增大時，蝸殼內靜壓處于降低狀態。從原因看，泵內靜壓升高主要由葉輪做功造成，而蝸殼內靜壓降低則由葉輪邊緣出流絕對速度降低所致。流量增大后，蝸殼出口的動壓實際上處于較大狀態，此時轉化為壓能的動能比較少，靜壓自然處于較低狀態。

除此之外，葉根附近呈現出局部低壓現象。例如，當流量為30m3/h時，葉根附近靜壓分布情況比較復雜，流體繞過葉根的過程中，流道的狀態變化方式是先收縮再擴張，因此流速在急劇提高的情況下發生了驟然下降，在此過程中出現了部分壓能轉變為動能的情況，所以造成了葉根附近呈現出局部根壓的現象。因此，在設計結構時，可以通過對流動情況的改變，合理的使低壓區面積縮小，達到對氣蝕現象的控制。不同轉速下的葉輪截面相對速度分布。

2.2.2 不同轉速下的計算結果分析

首先，結合通常的數值模擬與實驗分析經驗，將流量設置為50m3/h，然后，設置轉速共計4個：2000r/min、3000r/min、4000r/min、5000r/min。其次，通過對滑移網格區域轉角速度的改變，可以獲得不同轉速工況下的相對速度rν、絕對速度ν、靜壓力曲線等相關物理量的分布情況（見圖3）。分析如下：

圖3 不同轉速下的壓力流量曲線分布情況示意

首先，不同轉速下的葉輪截面相對速度分布情況表明，轉速提高時葉根附近的存在局部高速區域逐漸擴大的現象，可是葉輪附近大部分區域相對速度的變化不是很大，其次，不同轉速下的蝸殼截面絕對速度分布情況表明，轉速提高時，受到牽連速度的同步增大影響，蝸殼葉輪邊緣附近的速度隨之增高。最后，轉速提高時泵內靜壓隨之增大，葉輪邊緣與蝸殼中的靜壓增大情況比較明顯。從原因看，除了轉速提高時民生的科氏力增大影響外，還包括葉輪邊緣出流靜壓的升高、動能的增加。所以，在流量確定的情況下，動能隨轉速增加、動壓相同的情況下動能轉化為使泵內靜壓升高的壓能。

當進口壓力不變的情況下，對不同轉速下的壓力流量進行工況計算并制作出上圖中的離心泵壓力流量曲線。此時可以看出轉速確定的條件下，流量增加泵輸出壓力有明顯降低的現象。同時，流量確定的情況下轉速提高泵輸出壓力增大。另外，轉速降低離心泵則不能在大流量下保持正常工況。

3 結語

總之，本次選擇Fluent軟件對不同流量、轉速條件下的離心泵內流特性進行了數值模擬分析，結果顯示，增強對幾何形狀的設計、重視低壓區設計中葉輪載荷與氣餒影響，均有利于促進離心泵產品研發設計中的性能提升。同時，在應用過程中應該注重出口壓力與轉速限定條件下的運行流量變化，并在充分認識到轉速低時大流量運行困難的應用特點后按照轉速合理地控制其運行流量，確保離泵的應用性能始終處于良好狀態。