劉厚根,郭鵬濤
(中南大學 機電工程學院,長沙 410083)
機械增壓器是現(xiàn)代汽車發(fā)動機的重點研究方向之一,具有結構簡單、減少廢氣排放、瞬態(tài)響應快等優(yōu)點[1]。機械增壓器不僅應用于汽車行業(yè),而且已廣泛應用于通用飛機領域,對航空發(fā)動機的性能有顯著的提升。目前羅茨式機械增壓器已經(jīng)得到廣泛的應用,其工作原理與羅茨鼓風機相同,國內(nèi)學者對羅茨式機械增壓器進行了研究。李旺[2]等介紹了伊頓公司改進型羅茨機械增壓技術;劉厚根[3~6]等對四葉羅茨機械增壓器進行了性能實驗研究、轉子型線的設計和改進以及基于CFD的增壓器消減脈動方法研究;葉帥奇[7]等利用流體力學方法研究轉子型線對四葉羅茨機械增壓器內(nèi)部流場的影響;趙厚繼[8]等研究三葉轉子增壓器轉子系統(tǒng)的動態(tài)特性。
轉子的葉數(shù)是影響增壓器性能的因素之一,但是在現(xiàn)有文獻中,并沒有轉子葉數(shù)對羅茨機械增壓器性能影響的研究,因此本文借助動網(wǎng)格技術和數(shù)值模擬的方法,對常用的、大小相同的三葉和四葉羅茨機械增壓器進行研究,對結果進行分析,為羅茨機械增壓器數(shù)值模擬研究、轉子優(yōu)化設計和改進提供參考。如圖1為羅茨機械增壓器的主要結構。
圖1 羅茨機械增壓器主要結構
本文主要研究的對象是3.0L羅茨機械增壓器,三葉和四葉轉子羅茨機械增壓器轉子基本尺寸如表1所示。
表1 增壓器轉子基本結構尺寸
增壓器的實際流體域是復雜的三維,根據(jù)增壓器的工作原理和過程,忽略轉子與殼體的軸向間隙對增壓器內(nèi)泄露的影響,將流體域簡畫為二維,節(jié)省計算量。保證對容積效率影響最小的前提下,合理的簡畫出氣口和進氣口的形狀和大小,簡畫后的模型如圖2所示。為方便劃分網(wǎng)格和減少仿真中途出錯幾率,流體域中存在細小間隙,三葉和四葉流體域網(wǎng)格達9萬個,網(wǎng)格最小0.05mm網(wǎng)格質量達到動網(wǎng)格計算要求。劃分流體域中的網(wǎng)格采用三角形非結構化網(wǎng)格,在解決離散化繞流問題和網(wǎng)格再生問題有不可比擬的優(yōu)越性。
圖2 增壓器簡化流體域模型
由于機械增壓器進排氣呈現(xiàn)周期性變化,計算域中的網(wǎng)格隨時間發(fā)生變形和位移,只有動網(wǎng)格能實現(xiàn)這種狀況下的動態(tài)模擬[9]。隨著左、右轉子相互旋轉,轉子邊界位置時刻發(fā)生變化,網(wǎng)格會被重新劃分和重整化處理,以滿足要求。邊界的變形和運動用邊界函數(shù)定義,轉子在旋轉過程中不發(fā)生變形或變形量極小,故將動網(wǎng)格運動方式定義為剛體運動,轉子邊界設置為動網(wǎng)格邊界[10]。
左轉子邊界函數(shù)方程:
((rotation_left 3 point)
(time 0 0.1 1)
(omega_z w w w))
右轉子邊界函數(shù)方程:
((rotation_right 3 point)
(time 0 0.1 1)
(omega_z -w -w -w))
其中,omega_z表示轉子旋轉軸為z軸,w是轉子的角速度大小。
采用Fluent模擬二維狀態(tài)下機械增壓器的工作過程,需要對模型的邊界條件和初始條件設置[11]。將劃分好網(wǎng)格的三葉與四葉增壓器模型導入Fluent中進行瞬態(tài)分析,進氣口為壓力入口邊界,壓力值為標準大氣壓101Kpa;出氣口為壓力出口邊界,壓力值隨壓比不同而調整設定。殼體及左、右轉子邊界均設置為wall。殼體邊界設為熱對流方式,傳熱系數(shù)30W/(m2·K),溫度300K,壁面厚度7mm。
Fluent中提供的RNG/k-湍流模型適用于該仿真分析過程,此湍流模型使仿真過程更接近實際情況[3]。參考值(Reference Values)中將網(wǎng)格的深度(Depth)與轉子實際長度L設為一樣,使數(shù)值模擬結果更接近真實。求解器設置采用SIMPLE算法,有利于加快迭代步的收斂速度。
增壓器進排氣過程呈現(xiàn)周期性變化的,因此通過瞬態(tài)模擬轉子旋轉一周過程中流場變化研究增壓器的工作情況。在Fluent的求解設置中監(jiān)測轉子旋轉過程中增壓器排氣口的流量變化,研究增壓器排氣口流量變化,等步長記錄每次增壓器排氣口流量值數(shù)據(jù),經(jīng)數(shù)據(jù)處理獲得轉子旋轉過程中增壓器排氣口平均流量Qe。增壓器排氣口平均流量的計算:首先計算某轉速下,轉子旋轉一周的時間T;保證迭代計算順利進行,中途求解不出現(xiàn)負網(wǎng)格,經(jīng)過多次仿真驗證,設定迭代時間步長T/1500能滿足迭代步長的要求;最后每四步記錄一次該迭代時刻的排氣口流量大小,并將所得數(shù)據(jù)求平均數(shù),求得轉子旋轉一周過程中增壓器排氣口平均流量Qe,即作為增壓器容積效率計算中的排氣口流量Qs。
如圖3所示,轉子轉速3000r/min,進出口壓比1.3,轉子旋轉一周(0.02s)時,某一時刻增壓器轉子流場壓力分布情況。單個三葉轉子在旋轉一周的過程中,基元容積至少有60°是封閉旋轉過程;單個四葉轉子在旋轉的過程中,基元容積至少90°是封閉旋轉過程。一對三葉轉子在旋轉一周的過程中存在一個基元容積,一對四葉轉子在旋轉一周過程中至少存在兩個基元容積。由于進、排氣口壓差作用,氣體通過轉子間隙由排氣腔向進氣腔流動,發(fā)生內(nèi)泄露,影響增壓器的容積效率。從壓力分布圖看出,四葉轉子排氣腔到進氣腔之間形成較多的多級壓降,壓降的存在可以對內(nèi)泄露流動具有更大的阻礙作用。
圖3 壓力分布圖
圖4為Fluent中監(jiān)測的轉速為3500r/min,壓比1.5時,三葉和四葉轉子增壓器轉子旋轉一周過程中排氣口流量Q的變化情況。三葉轉子增壓器在啟動經(jīng)過約T/6后,排氣口流量Q隨時間規(guī)律性變化,進入穩(wěn)定狀態(tài),旋轉一周過程中出現(xiàn)6次諧波變化;四葉轉子增壓器在啟動經(jīng)過約T/8后排氣口流量Q隨時間規(guī)律性變化,進入穩(wěn)定狀態(tài),旋轉一周過程中出現(xiàn)8次諧波變化,這是左右轉子交替作用的結果。增壓器不同時刻的排氣口流量Q在平均值附近上下波動。
圖4 排氣口流量隨時間變化圖
利用流量不均勻度來描述排氣口流量脈動情況:
式(1)中:Qmax為瞬時流量最大值,Qmin為瞬時流量最小值,Qave為平均流量。
由于數(shù)值模擬的起始階段排氣口流量波動不穩(wěn)定,流量不均勻度計算從排氣口流量隨時間變化穩(wěn)定后,轉子旋轉一周過程中取值計算。根據(jù)公式(1)知,壓比1.5,轉速3500r/min時,三葉轉子的流量不均系數(shù)為8.26,四葉轉子的流量不均系數(shù)為3.10,四葉轉子增壓器排氣口流量脈動高于三葉轉子增壓器。如圖5所示,不同轉速和壓比下,三葉與四葉轉子增壓器排氣口流量不均勻度。由圖5知,隨著壓比的增大,增壓器排氣口流量不均度逐漸增大。三葉轉子增壓器和四葉轉子增壓器的壓比為1.0,1.1時,兩者排氣口流量不均度相差較??;三葉轉子增壓器和四葉轉子增壓器在壓比為1.3,1.5時,三葉轉子增壓器的排氣口流量不均度大,在低速時脈動更明顯。主要原因是,相同轉速和壓比下,四葉轉子增壓器的排氣口均壓時間短于三葉轉子增壓器,排氣口流量脈動幅度小,脈動小,氣流噪聲小。
圖5 不同轉速和壓比下排氣口流量不均勻度
不考慮壓縮和泄露等實際因素影響,增壓器在單位時間內(nèi)輸送的氣體容積,稱為理論流量[12],即:
式(1)中:Qth為理論流量;λ為轉子面積利用系數(shù);D為轉子外徑;n為轉子轉速。
機械增壓器的容積效率是排氣口流量Qs與理論流量Qth之比,以η表示。數(shù)值模擬轉子轉速在2000~5000r/min,壓比為1.0~1.5情況下,三葉和四葉轉子增壓器排氣口流量Qs,獲得三葉和四葉轉子增壓器在不同轉速和壓比情況下的容積效率η,如圖6所示。
圖6 增壓器數(shù)值模擬容積效率
從圖可見,壓比大于1,轉速2000~5000r/min時,四葉轉子增壓器的容積效率高于三葉轉子增壓器,在2000~3000r/min時,四葉轉子增壓器的容積效率高于三葉轉子增壓器10~15%,可知,四葉轉子增壓器低速增壓性能優(yōu)于三葉轉子增壓器。主要原因是轉子葉數(shù)增多,轉子的面積利用系數(shù)增大,理論排氣口流量變大,轉子在旋轉過程中,對氣體內(nèi)泄露阻礙作用越大,效率隨之增大。但隨著轉子轉速的增加,三葉和四葉轉子增壓器的容積效率趨于相同。
本實驗分別針對自主設計的大小相同的一臺四葉增壓器樣機和一臺的三葉增壓器開展性能試驗,性能試驗參考一般羅茨鼓風機性能試驗方法標準[13]。機械增壓器性能試驗原理圖,如圖7所示,機械增壓器由變頻電機通過多楔帶直接驅動,通過調節(jié)變頻器電機轉速改變機械增壓器的轉速。在排氣口安裝一個壓力調節(jié)閥來控制排氣壓力。排氣口管道的末端加裝消聲器,模擬發(fā)動機及其排氣系統(tǒng)的阻抗消聲作用。電機和機械增壓器之間安裝轉速轉矩傳感器,記錄增壓器的轉速,扭矩數(shù)據(jù)。機械增壓器與壓力調節(jié)閥之間安裝壓力和溫度傳感器,記錄增壓器排氣口的壓力和溫度數(shù)據(jù)。
圖7 試驗系統(tǒng)原理
增壓器性能試驗主要測試轉速在2000~5000r/min,壓比在1.0,1.3情況下,兩臺增壓器的排氣口流量,經(jīng)處理計算,獲得兩臺增壓器的性能試驗容積效率,如圖8所示。性能試驗獲的增壓器容積效率變化趨勢與數(shù)值模擬獲得增壓器容積效率變化趨勢相同,四葉轉子增壓器低速增壓性能更好,證明數(shù)值模擬的可行性。
圖8 增壓器性能試驗容積效率
1)壓比1.0,1.1時,三葉轉子和四葉轉子增壓器排氣流量不均度基本一致;壓比1.3,1.5時,三葉轉子增壓器的排氣口流量不均度大于四葉轉子增壓器,并隨轉速的增加兩者差距呈縮小趨勢。轉子葉數(shù)越多增壓器每次均壓時間越短,排氣流量的脈動幅度越小,氣流噪聲小。
2)四葉轉子增壓器在低速階段,容積效率高于三葉轉子增壓器較多,說明四葉轉子低速增壓性能優(yōu)于三葉轉子增壓器,四葉轉子增壓器對于改善發(fā)動機的低速扭矩特性優(yōu)勢明顯,可為三葉轉子與四葉轉子增壓器與發(fā)動機性能匹配提供依據(jù)。
3)通過機械增壓器性能試驗,驗證了數(shù)值模擬的可行性,為以后增壓器轉子優(yōu)化設計制造提供一種方法。