基于格子Boltzmann 方法的動壓氣體軸承黏性熱耗散研究*

鹿翔宇 許 波 陳振乾

（東南大學(xué)能源與環(huán)境學(xué)院 江蘇南京 210096）

動壓氣體軸承因其轉(zhuǎn)速高、 壽命長、 低摩擦、 無污染等優(yōu)點(diǎn)， 被廣泛應(yīng)用于透平膨脹機(jī)、 渦輪壓氣機(jī)、 飛機(jī)發(fā)動機(jī)、 制冷機(jī)等高速旋轉(zhuǎn)渦輪機(jī)械中， 在航空航天、 能源與動力等領(lǐng)域發(fā)揮著十分重要的作用［1］。 在動壓氣體軸承的運(yùn)轉(zhuǎn)過程中， 雖然轉(zhuǎn)子和箔片不會直接接觸， 但是軸承間隙內(nèi)流體的高速剪切流動依然會產(chǎn)生大量的熱量， 容易引起溫度的分布不均和局部過熱， 進(jìn)而燒蝕箔片。 此外， 軸承的零部件也會因溫度的升高而導(dǎo)致剛度和阻尼發(fā)生變化， 產(chǎn)生熱變形， 對軸承的性能造成一定的影響［2］。 因此， 研究動壓氣體軸承的氣動熱和熱特性規(guī)律， 以及相關(guān)的冷卻方法具有重要意義。

國內(nèi)外學(xué)者對于軸承熱特性的研究主要集中在近20 年。 SALEHI 等［3-4］使用基于Couette 近似的簡化能量方程和可壓縮雷諾方程獲得了箔片氣體軸承中的氣膜溫度， 并采用有限差分法和迭代法對簡化表達(dá)式進(jìn)行求解。 PENG 和KHONSARI［5-6］考慮空氣的可壓縮性和黏溫特性以及軸承表面的柔度， 建立了用于預(yù)測箔片氣體軸承三維溫度場的熱流體力學(xué)模型， 并研究了冷卻流量的影響， 發(fā)現(xiàn)膜的溫升隨轉(zhuǎn)速的增加呈指數(shù)增長。 LEE 等［7-8］提出了一個(gè)包含箔片結(jié)構(gòu)和轉(zhuǎn)子的軸承傳熱模型， 計(jì)算了潤滑氣膜、 頂箔、 波箔、 軸承套和轉(zhuǎn)子的溫度， 并采用廣義雷諾方程和三維能量方程預(yù)測了氣膜溫度。 SIM 和KIM［9］通過軸的二維軸對稱熱傳導(dǎo)、 氣膜中的三維能量傳輸、 箔片熱阻和軸承殼體中的三維熱傳導(dǎo)模型， 建立了箔片氣體軸承的熱流體動力學(xué)模型， 該模型還結(jié)合了箔片接觸熱阻和入口流混合的分析模型， 以提高模型精度。 SAMANTA和KHONSARI［10］分析了箔片氣體軸承的熱彈性不穩(wěn)定問題， 推導(dǎo)出一個(gè)封閉的解析解來預(yù)測在波箔中發(fā)生熱彈性不穩(wěn)定的臨界速度， 并驗(yàn)證了波數(shù)和膜厚對失穩(wěn)臨界速度的顯著影響。 ZDZIEBKO 和MARTOWICZ［11］開發(fā)了箔片氣體軸承的有限元模型， 分析表明， 波箔與頂箔和軸承套接觸次數(shù)最多的區(qū)域， 也就是頂箔的熱量傳遞到軸承襯套和其他軸承的部件最集中的區(qū)域。

以上對于箔片氣體軸承熱特性的研究， 主要采用的是雷諾方程結(jié)合能量守恒方程的宏觀方法。 然而軸承間隙的尺寸一般非常小， 通常在幾十微米的級別［12］， 間隙內(nèi)的流場已經(jīng)處于滑移區(qū)， 甚至過渡區(qū)，會出現(xiàn)一些微觀效應(yīng)， 例如速度滑移、 溫度躍變、 黏性熱耗散等， 此時(shí)采用一些微觀方法將更加合適。 分子動力學(xué)（Molecular Dynamics， MD） 模擬和直接模擬蒙特卡羅（Direct Simulation Monte-Carlo， DSMC）方法等基于粒子的方法在微流場模擬方面取得了一定的進(jìn)展［13-14］。 然而， 對于大多數(shù)實(shí)際應(yīng)用來說， 它們的計(jì)算量過于龐大。 介觀的格子玻爾茲曼方法（Lattice Boltzmann Method， LBM）［15］自從被提出以來， 已經(jīng)被廣泛應(yīng)用于微尺度流動的研究， 其計(jì)算量與宏觀方法相當(dāng)， 比較適合用于文中的研究。 與經(jīng)典雷諾潤滑方程相比， LBM 可以捕捉更加微觀的現(xiàn)象，邊界條件實(shí)施簡單。 并且當(dāng)軸承間隙內(nèi)氣流處于過渡區(qū)時(shí)， 此時(shí)不滿足連續(xù)性假設(shè)， 宏觀雷諾方程不再適用， LBM 將成為十分有潛力的研究方法。 然而鮮有學(xué)者將LBM 用于動壓氣體軸承的黏性熱耗散研究。

目前， 格子Boltzmann 方法應(yīng)用在熱流動研究的主要有3 種模型： 多速度 （Multi-speed， MS） 模型［16-17］， 雙分布函數(shù)（Double-Distribution Function，DDF） 模型［18-23］和混合模型［24］。 MS 模型是等溫LB模型的直接延伸， 其利用密度分布函數(shù)的高階速度矩來描述能量方程。 DDF 模型利用了2 種不同的分布函數(shù)， 一種用于流場， 另一種用于溫度場。 在混合方法中， 流動模擬采用LB 法， 而溫度場的求解則采用傳統(tǒng)的數(shù)值方法， 如有限差分法。 1998 年， HE 等［18］基于DDF 模型引入內(nèi)能分布函數(shù)， 第一次提出了包含黏性熱耗散和壓縮功的熱格子Boltzmann 模型， 隨后有學(xué)者在其基礎(chǔ)上進(jìn)行了一定的改進(jìn)［20，25］。 2007年， GUO 等［21］在單一速度分布函數(shù)的基礎(chǔ)上定義了代表總能量的分布函數(shù)， 基于DDF 模型提出了總能形式的含有黏性熱耗散和壓縮功的熱LBE 模型。 該方法在演化過程中不涉及復(fù)雜的導(dǎo)數(shù)項(xiàng)， 計(jì)算簡單，數(shù)值穩(wěn)定性高， 文中將基于該模型研究動壓氣體軸承的熱特性。

LBM 雖然在研究熱流動方面已經(jīng)發(fā)展了20 多年， 但是目前相關(guān)研究依然是針對一些比較簡單的物理模型， 像軸承間隙這種彎曲變截面的高速剪切微通道鮮有研究。 因此本文作者將LB 方法應(yīng)用到該領(lǐng)域， 基于GUO 等［21］提出的總能形式的DDF 離散速度模型， 得到了可以應(yīng)用在貼體網(wǎng)格中的有限差分格式的顯式LB 方法， 并引入Langmuir 滑移邊界， 實(shí)現(xiàn)對軸承間隙微通道的模擬， 并研究了不同溫差大小、 偏心率、 轉(zhuǎn)速以及溫度階躍對于黏性熱耗散分布規(guī)律的影響， 從更加微觀的層面探究了軸承間隙黏性熱耗散的機(jī)制。

1 數(shù)值模型

1.1 有限差分DDF-LBE 模型

文中使用的包含黏性熱耗散和壓縮功的總能形式雙分布函數(shù)離散速度模型［21］， 其基于Hermit 多項(xiàng)式展開， 并保留至二階項(xiàng)。 它由密度分布函數(shù)和總能分布函數(shù)構(gòu)成， 其中總能包括內(nèi)能和機(jī)械能， 2 個(gè)分布函數(shù)的演化方程為

式中：fi和hi分別是i方向的密度分布函數(shù)和總能分布函數(shù)；和是對應(yīng)的平衡態(tài)分布函數(shù)；ci是粒子的離散速度；τhf＝1／τh-1／τf， 其中τf和τh分別為動能和內(nèi)能松弛時(shí)間；Zi＝ci·u-u2／2，u是宏觀速度；Fi和qi是與外力有關(guān)的2 個(gè)項(xiàng)。

對于二維九速（D2Q9） 模型， 2 個(gè)平衡態(tài)分布函數(shù)為

式中：R為氣體常數(shù)；T0為參考溫度；ρ為宏觀密度，p0＝ρRT0； 權(quán)重系數(shù)w0＝4／9，w1＝w2＝w3＝w4＝1／9，w5＝w6＝w7＝w8＝1／36， 離散速度ci取值為

2 個(gè)外力項(xiàng)Fi和qi的表達(dá)式為

式中：a為質(zhì)量力， 總能，cv為定容比熱容。

為穩(wěn)定可靠地求解離散方程（1）、 （2）， 對其進(jìn)行有限差分格式的離散， 在時(shí)間［tn，tn＋1］ 內(nèi)進(jìn)行求積， 并對演化方程中的碰撞項(xiàng)做如下處理：

式中：tn＋1＝tn＋Δt； 碰撞算子，，Si＝-τfZiΩf／τhf＋qi。

為了能夠顯式求解上式， 可以引入2 個(gè)新的分布函數(shù)，

式中：ωf＝Δt／τf，ωh＝Δt／τh。

將新的分布函數(shù)代入式（8）、 （9）， 就能得到可應(yīng)用于貼體網(wǎng)格中的有限差分形式的玻爾茲曼方程：

動能和內(nèi)能的松弛時(shí)間有如下關(guān)系：

式中：μ為動力黏度；p是壓力；κ為熱傳導(dǎo)系數(shù)；cp是定壓比熱容。

于是普朗特?cái)?shù)為Pr＝μcp／κ＝τf／τh。

為了能夠?qū)⑸鲜龇匠虘?yīng)用到曲線坐標(biāo)系ξ中， 以方程（12） 為例， 對其空間梯度項(xiàng)進(jìn)行如下離散［26］：

?fi／?ξβ的離散理論上可以采用中心差分格式或者二階迎風(fēng)插值格式進(jìn)行求解：

式中： Δξβ是ξβ方向的網(wǎng)格間距。

中心差分格式數(shù)值耗散少， 但不穩(wěn)定， 而迎風(fēng)插值格式雖相對穩(wěn)定， 但是卻有較強(qiáng)的數(shù)值耗散。 為了解決這些問題， 可以將上述2 種方法組成一種混合的差分格式：

其中， 0≤?≤1 是調(diào)節(jié)中心差分格式和迎風(fēng)插值格式權(quán)重的系數(shù)。 根據(jù)文獻(xiàn)的建議， 在后續(xù)研究中無特殊說明時(shí)取0.1。

為了讓LB 方法可以模擬微尺度流動現(xiàn)象， 需要處理2 個(gè)重要問題， 第一個(gè)就是解決松弛時(shí)間τf和流動Kn數(shù)（Kn＝λ／L，λ為分子平均自由程，L為特征長度） 之間的關(guān)系， 通過結(jié)合動力學(xué)理論和LBM的運(yùn)動黏度表達(dá)式， 可得［27］：

式中：H為通道的高度。

由于軸承間隙內(nèi)的流場壓力和溫度分布不均， 因此對上式進(jìn)行如下修正：

式中：ρ0為參考密度；T0為參考溫度。

1.2 邊界條件

將LB 方法應(yīng)用到熱微流動的另一個(gè)重要問題就是如何處理速度滑移和溫度階躍邊界條件。 文中采用了Langmuir 滑移模型［28］， 它反映了氣體分子和壁面之間的界面相互作用， 可以表示為

式中： 下標(biāo)w 代表壁面處的局部值； 下標(biāo)g 代表靠近壁面處的局部值；αv和αT是與氣體種類和壁面材料屬性有關(guān)的調(diào)節(jié)系數(shù)。

為了方便確定這2 個(gè)調(diào)節(jié)系數(shù)， 可以將Langmuir滑移模型與Maxwell 滑移模型進(jìn)行聯(lián)立

式中：σv為動量調(diào)節(jié)系數(shù)；σT為熱調(diào)節(jié)系數(shù)；σv和σT在文中沒特殊說明的情況下均取1；γ是比熱容比。

于是

式中：N＝1／δx表示網(wǎng)格數(shù)量。

將Langmuir 模型的處理方法應(yīng)用到非平衡外推邊界中， 最終可以得到適用于當(dāng)前模型的速度滑移和溫度階躍的邊界格式：

式中：xw是邊界處的節(jié)點(diǎn)；xg是與邊界相鄰的節(jié)點(diǎn)。

1.3 求解步驟

數(shù)值計(jì)算程序的大致流程如圖1 所示。

圖1 計(jì)算流程Fig.1 Numerical calculation procedure

具體的實(shí)現(xiàn)步驟如下：

（1） 通過求解橢圓偏微分方程組生成貼體網(wǎng)格，并計(jì)算轉(zhuǎn)置速度。

（2） 根據(jù)物理問題， 初始化網(wǎng)格節(jié)點(diǎn)的密度、速度和溫度， 得到初始狀態(tài)的密度分布函數(shù)和總能分布函數(shù)及其平衡態(tài)。

（3） 根據(jù)Kn數(shù)確定τf， 進(jìn)一步得到ωf， 同理根據(jù)Pr數(shù)確定τh， 進(jìn)一步得到ωh。

（4） 密度分布函數(shù)的求解： 根據(jù)式（14） 計(jì)算密度和速度， 由結(jié)合式（10） 可得fni， 隨后根據(jù)式（19） 計(jì)算fni的有限差分， 最后根據(jù)式（12） 計(jì)算得到。

（5） 總能分布函數(shù)的求解： 根據(jù)式（14） 計(jì)算溫度， 由結(jié)合式 （11） 可得hni， 隨后根據(jù)式（19）計(jì)算hni的有限差分， 最后根據(jù)式（13） 計(jì)算得到。

（6） 判斷速度和溫度是否滿足收斂條件， 如果不滿足則根據(jù)溫度和密度調(diào)整Kn數(shù)， 并重復(fù)步驟（3） — （6）， 如果滿足則輸出結(jié)果并停止計(jì)算。

2 物理模型

簡化后的動壓氣體軸承的示意圖如圖2 所示。 軸頸和軸承套的半徑分別為r1和r2（r2＝r1＋C0， 其中C0為平均間隙）， 最大間隙為hmax， 最小間隙為hmin。軸頸以速度ω順時(shí)針旋轉(zhuǎn)，e為偏心距， 則偏心率為ε＝e／C0， 軸承的Kn數(shù)定義為Kn＝λ／hmin。 軸頸表面的溫度為T1， 軸承內(nèi)表面的溫度為T2， 兩壁面的溫差ΔT（＝T1-T2） 可以采用無量綱數(shù)Eckert 數(shù)來表征，其定義為

圖2 動壓氣體徑向軸承簡化示意Fig.2 Schematic of gas-lubricated journal bearing

式中：U為軸頸表面的線速度。

軸承的具體參數(shù)和運(yùn)行工況范圍見表1。

表1 軸承參數(shù)和運(yùn)行工況Table 1 Bearing parameters and operating conditions

3 模型與可靠性驗(yàn)證

以圓柱熱Couette 流為例， 對文中模型進(jìn)行了驗(yàn)證。 其幾何參數(shù)和貼體網(wǎng)格劃分如圖3 所示。 內(nèi)外圓柱半徑比Rβ＝r1／r2＝3／5， 內(nèi)圓柱以線速度U順時(shí)針旋轉(zhuǎn)， 溫度為T1， 外圓柱靜止， 溫度為T2。 考慮黏性熱耗散， 圓周方向和半徑方向的網(wǎng)格數(shù)分別為NX＝188，NY＝20。 圓柱熱Couette 流的量綱一溫度T*的解析解表達(dá)式［29］為

圖3 圓柱熱Couette 流示意Fig.3 Schematic of thermal cylindrical Couette flow

式中：r*＝r／r1， 為量綱一半徑。

文中在Pr＝0.7 的情況下， 分別模擬得到了Ec數(shù)為1、 5、 10 的3 種徑向溫度分布， 結(jié)果見圖4。 可以發(fā)現(xiàn)， 在相同網(wǎng)格分辨率下， 數(shù)值模擬結(jié)果與理論解吻合良好， 最大相對誤差不超過0.8%， 證明模型精度較高。

圖4 不同Ec 數(shù)下的溫度分布Fig.4 Temperature distribution at different Ec number

在后續(xù)的模擬中， 為了排除網(wǎng)格數(shù)量對于軸承間隙黏性熱耗散計(jì)算結(jié)果的影響， 還需進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證。 文中對軸承半徑為8 mm， 間隙尺寸為25 μm，偏心率為0， 轉(zhuǎn)速為5×104r／min 的動壓氣體軸承間隙熱流動進(jìn)行模擬，Ec數(shù)固定為10。 間隙流域圓周方向和半徑方向的網(wǎng)格數(shù)分別為10 053×5、 20 106×10、 30 159×15、 40 212×20、 50 265×25。 5 種網(wǎng)格分辨率下模擬得到的間隙流場中最高的量綱一氣膜溫度如圖5 所示。 可以看出， 隨著網(wǎng)格數(shù)量的增加， 黏性熱耗散產(chǎn)生的最高溫度也隨之升高， 但是升高的速率逐漸減小， 說明增加網(wǎng)格數(shù)對于提升模型精度的作用逐漸降低。 綜合考慮模型精度和計(jì)算資源的消耗， 文中最終選取的網(wǎng)格分辨率為30 159×15。

圖5 網(wǎng)格獨(dú)立性驗(yàn)證Fig.5 Grid independence verification

4 結(jié)果與討論

4.1 不同Ec 數(shù)的影響

軸承在運(yùn)轉(zhuǎn)過程中， 由于散熱條件的不同， 一般軸頸表面和軸承套內(nèi)表面的溫度是存在差異的。 為了更好地探究黏性熱耗散的變化規(guī)律， 假設(shè)軸頸和軸承套兩側(cè)為定溫壁面， 且存在一定的溫差， 并用Ec數(shù)來表示溫差大小。 軸承轉(zhuǎn)速為5×104r／min， 偏心率保持0.6， 圖6 所示為不考慮兩壁面處溫度階躍時(shí)，不同Ec數(shù)下間隙內(nèi)最高氣膜溫度處徑向量綱一溫度T*（＝（T-T2）／（T1-T2）） 的分布。 可以看出， 在Ec數(shù)較小時(shí)， 兩側(cè)溫差較大， 間隙溫度分布幾乎呈線性分布， 此時(shí)是以熱傳導(dǎo)為主， 黏性耗散的影響基本可以忽略。 而當(dāng)Ec數(shù)較大時(shí)， 黏性熱耗散的作用明顯， 溫度呈拋物線分布， 流場中最高溫度出現(xiàn)在靠近軸頸側(cè)， 且流場中大部分的溫度都已經(jīng)高于T1。 總體來看， 隨著Ec數(shù)的增加， 也就是溫差的減小， 黏性耗散作用增強(qiáng)， 氣膜溫度的最大值是逐漸遞增的。

圖6 不考慮溫度階躍時(shí)不同Ec 數(shù)下最高氣膜溫度處徑向量綱一溫度分布Fig.6 Radial dimensionless temperature distribution at the highest gas film temperature for different Ec number without considering temperature jump

圖7 所示為考慮了溫度階躍情況下的不同Ec數(shù)時(shí)的溫度分布。 當(dāng)Ec＝1 時(shí)， 軸頸側(cè)的溫度小于T1，溫度階躍為正值； 而當(dāng)Ec＝2 時(shí)， 軸頸側(cè)的溫度大于T1， 此時(shí)溫度階躍變成了負(fù)值。 隨著Ec數(shù)的增加，兩側(cè)的溫度階躍均呈增大的趨勢， 說明溫度階躍的影響將變得更加顯著。 因此， 軸頸和軸承套兩側(cè)溫差的增大， 將更有利于削弱黏性熱耗散的影響。

圖7 考慮溫度階躍時(shí)不同Ec 數(shù)下最高氣膜溫度處徑向量綱一溫度分布Fig.7 Radial dimensionless temperature distribution at the highest gas film temperature for different Ec number considering temperature jump

4.2 不同偏心率的影響

圖8 所示為轉(zhuǎn)速5×104r／min， 軸承平均間隙為25 μm，Ec數(shù)為10， 不考慮溫度階躍時(shí)不同偏心率下間隙內(nèi)最高氣膜溫度處徑向量綱一溫度分布。 可以看出， 隨著偏心率的增加， 間隙中的氣膜最高溫度也逐漸升高。 這是因?yàn)椋?當(dāng)偏心率變大時(shí)， 軸承間隙尺寸在圓周方向的變化幅度增加， 意味著氣膜的擠壓作用增強(qiáng)， 黏性熱耗散的量也隨之增大。 此外， 出現(xiàn)最大氣膜溫度處的位置隨著偏心率的增加逐漸由靠近軸頸側(cè)向流場中間偏移， 表明間隙內(nèi)的流動逐漸接近于平板Couette 流。

圖8 不考慮溫度階躍時(shí)不同偏心率下最高氣膜溫度處徑向量綱一溫度分布Fig.8 Radial dimensionless temperature distribution at the highest gas film temperature for different eccentricity ratios without considering temperature jump

圖9 展示了考慮溫度階躍時(shí)不同偏心率下最高氣膜溫度處徑向量綱一溫度分布。 隨著偏心率的增加，軸頸和軸承套兩側(cè)的溫度階躍均增大， 并且軸頸側(cè)的溫度階躍量要小于軸承套側(cè)的溫度階躍量， 原因是軸頸側(cè)的溫度梯度要小于軸承套側(cè)的溫度梯度。 由于溫度階躍的發(fā)生， 當(dāng)偏心率為0.8 時(shí)， 軸承套側(cè)的溫度已十分接近軸頸側(cè)的溫度， 流場中的溫度基本都在T1之上， 因此在大偏心率情況下， 溫度階躍的影響不容忽略。

圖9 考慮溫度階躍時(shí)不同偏心率下最高氣膜溫度處徑向量綱一溫度分布Fig.9 Radial dimensionless temperature distribution at the highest gas film temperature for different eccentricity ratios considering temperature jump

為了更好地看出溫度階躍對于黏性熱耗散的影響， 圖10 給出了不同偏心率下氣膜溫度最大值的變化。 可以發(fā)現(xiàn)， 最高氣膜溫度的值隨著偏心率的增加而線性遞增， 并且考慮溫度階躍時(shí)， 溫度增加的速率更大， 再次印證了大偏心率情況下， 溫度階躍效應(yīng)的影響顯著。

圖10 不同偏心率下的最高氣膜溫度變化Fig.10 Variation of maximum gas film temperature at different eccentricity ratios

4.3 不同轉(zhuǎn)速的影響

圖11 所示為偏心率ε＝0.6， 軸承平均間隙為25 μm，Ec數(shù)為10， 不考慮溫度階躍時(shí)不同轉(zhuǎn)速下間隙內(nèi)最高氣膜溫度處徑向量綱一溫度分布。 可以看出，當(dāng)轉(zhuǎn)速從3×104r／min 增加到1.1×105r／min 時(shí)， 間隙內(nèi)氣膜溫度的最高值增加明顯。 這是因?yàn)椋?隨著軸承轉(zhuǎn)速的提升， 間隙內(nèi)流體的速度梯度增加， 流體剪切力增大， 黏性熱耗散的作用增強(qiáng)。 此外， 雖然間隙的幾何尺寸未發(fā)生變化， 但是隨著轉(zhuǎn)速的變化， 間隙內(nèi)出現(xiàn)最大氣膜溫度的位置也在隨之改變， 可能是因?yàn)殡x心力的變化導(dǎo)致的流體流動和傳熱特性的改變。

圖11 不考慮溫度階躍時(shí)不同轉(zhuǎn)速下最高氣膜溫度處徑向量綱一溫度分布Fig.11 Radial dimensionless temperature distribution at the highest gas film temperature at different rotational speeds without considering temperature jump

圖12 所示為考慮溫度階躍時(shí)不同轉(zhuǎn)速下最高氣膜溫度處徑向量綱一溫度分布。 可以看出， 隨著轉(zhuǎn)速的增加， 軸頸和軸承套兩側(cè)的溫度階躍均明顯增大。 當(dāng)轉(zhuǎn)速為7×104r／min 時(shí)， 軸承套側(cè)的溫度已十分接近T1， 而當(dāng)轉(zhuǎn)速為9×104和1.1×105r／min 時(shí)，軸承套側(cè)的溫度將大于軸頸側(cè)的溫度T1， 此時(shí)間隙內(nèi)整個(gè)流場的溫度都高于T1， 這是由于溫度階躍和黏性熱耗散共同作用的結(jié)果。

圖12 考慮溫度階躍時(shí)不同轉(zhuǎn)速下最高氣膜溫度處徑向量綱一溫度分布Fig.12 Radial dimensionless temperature distribution at the highest gas film temperature at different rotational speeds considering temperature jump

圖13 給出了不同轉(zhuǎn)速下間隙氣膜最高溫度的變化趨勢。 氣膜溫度最大值隨著轉(zhuǎn)速的增加而增加， 并且升高的速率逐漸增大。 當(dāng)考慮溫度階躍時(shí)， 最高氣膜溫度的增加速度和不考慮溫度階躍時(shí)比較接近， 說明不考慮溫度階躍對于不同轉(zhuǎn)速下黏性熱耗散計(jì)算帶來的低估影響大致一樣。

圖13 不同轉(zhuǎn)速下的最高氣膜溫度變化Fig.13 Variation of maximum gas film temperature at different rotational speeds

5 結(jié)論

（1） 軸頸和軸承套兩側(cè)溫差較大時(shí)， 間隙溫度幾乎呈線性分布， 以熱傳導(dǎo)為主； 當(dāng)軸頸和軸承套兩側(cè)溫度較小時(shí)， 間隙溫度呈拋物線分布， 此時(shí)黏性熱耗散作用明顯。

（2） 隨著偏心率的增加， 最大氣膜溫度線性遞增， 軸頸和軸承套兩側(cè)的溫度階躍也隨之增加， 且偏心率越大， 溫度階躍對于黏性熱耗散量計(jì)算的影響越大。

（3） 隨著轉(zhuǎn)子轉(zhuǎn)速的增加， 最大氣膜溫度升高，變化的速率逐漸增大， 軸頸和軸承套兩側(cè)的溫度階躍增大， 但溫度階躍對于不同轉(zhuǎn)速下黏性熱耗散量計(jì)算的影響程度相差不大。