基于Workbench的高壓圓盤氣體軸承共軛傳熱研究

郭良斌，吳永良

(1.武漢科技大學(xué)冶金裝備及其控制教育部重點(diǎn)實(shí)驗(yàn)室，湖北武漢 430081；2.武漢科技大學(xué)機(jī)械傳動(dòng)與制造工程湖北省重點(diǎn)實(shí)驗(yàn)室，湖北武漢 430081)

隨著氣體軸承向高轉(zhuǎn)速、高供氣壓力方向發(fā)展，氣膜中熱量的產(chǎn)生和遷移愈加顯著，不僅導(dǎo)致軸承及相關(guān)部件的熱變形，而且還改變了軸承在實(shí)際工況下的工作間隙，進(jìn)而使軸承的實(shí)際特性偏離理想設(shè)計(jì)狀態(tài)[1]。因此，氣體軸承的熱穩(wěn)定性逐漸成為制約其工作性能正常發(fā)揮的一個(gè)關(guān)鍵因素[2-6]。郭良斌和徐凡[7]從二維氣體動(dòng)力學(xué)方程組出發(fā)，提出了一種雙對(duì)稱收縮段設(shè)計(jì)理論，采用該收縮段向靜壓圓盤氣體軸承供氣，可以使圓盤間隙出口處的氣體流速達(dá)到超音速，軸承承載力可隨供氣壓力的提高線性增加，供氣壓力可以達(dá)到2.03 MPa(20 atm)以上。供氣壓力的提高使得圓盤間隙中的氣體壓力和流速迅速增加，軸承流道內(nèi)氣體從低速發(fā)展到高亞音速、甚至超音速狀態(tài)，氣膜間隙內(nèi)部各部分之間、以及氣膜與圓盤固壁之間均存在明顯的溫度差，軸承圓盤內(nèi)部熱傳導(dǎo)與流道間隙內(nèi)高亞音速或超音速氣流的對(duì)流換熱過程通過圓盤壁面緊密耦合在一起，是一個(gè)典型的共軛傳熱問題。

共軛傳熱是指包含固體內(nèi)部熱傳導(dǎo)與流經(jīng)該固體表面的流體的自由、強(qiáng)制或混合對(duì)流之間的相互作用的一種復(fù)雜傳熱過程[8]。共軛傳熱現(xiàn)象在科學(xué)和工程領(lǐng)域中大量存在，如核反應(yīng)系統(tǒng)[9-12]，復(fù)合式氣膜沖擊冷卻[13-16]，高速飛行器[17-20]，以及油膜軸承等[21-22]。

目前共軛傳熱方法在工程領(lǐng)域有一定程度的應(yīng)用[10-11，23-24]，受到越來越多的關(guān)注[25]。高壓氣體軸承流道間隙內(nèi)氣流溫度變化對(duì)軸承圓盤的變形有直接的影響，因此研究軸承流道間隙內(nèi)部和軸承圓盤的溫度場(chǎng)，以及流固耦合界面上的熱流密度分布規(guī)律是十分必要的。

本文作者使用ANSYS Workbench 2020R2工作平臺(tái)的Fluid Flow(Fluent)模塊，對(duì)高壓圓盤氣體軸承進(jìn)行流體域與固體域的網(wǎng)格劃分，設(shè)置邊界條件并使用該模塊的Fluent程序進(jìn)行共軛傳熱數(shù)值模擬，得到了該軸承共軛傳熱條件下流場(chǎng)參數(shù)分布、流體域與固體域的溫度分布及流固耦合壁面上的熱流密度分布，并將其與非共軛傳熱恒溫壁面條件下流體域的流場(chǎng)參數(shù)計(jì)算結(jié)果進(jìn)行了對(duì)比，得到了高壓圓盤氣體軸承共軛傳熱的一些基本特性，為該類軸承的設(shè)計(jì)和制造提供了有益的指導(dǎo)。

1 軸承結(jié)構(gòu)與邊界命名

1.1 軸承結(jié)構(gòu)

高壓圓盤氣體軸承的結(jié)構(gòu)如圖1所示。高壓圓盤氣體軸承由上、下工作圓盤組件組成，兩圓盤設(shè)計(jì)工作間隙h1為0.2 mm。上、下工作圓盤組件關(guān)于氣膜對(duì)稱面14呈平面對(duì)稱，同時(shí)關(guān)于對(duì)稱軸線1呈軸對(duì)稱。以上工作圓盤組件為例，其由上穩(wěn)流管3、上密封墊圈5和上圓盤6通過螺紋連接緊固在一起，上穩(wěn)流管3外端與供氣管道的管接頭相連。上穩(wěn)流管3內(nèi)部加工有軸對(duì)稱收縮段2，目的是對(duì)供氣管道輸送進(jìn)來的氣體進(jìn)行整流，使氣體以均勻的速度進(jìn)入上供氣孔4。氣體從穩(wěn)流管3、10進(jìn)入，經(jīng)供氣孔4、11和上、下圓盤收縮段7、13及平行氣膜區(qū)15，最終從氣膜間隙出口流出至環(huán)境大氣。上、下圓盤收縮段7、13的作用是使軸向來流無分離地轉(zhuǎn)變?yōu)槠叫杏趫A盤壁面的徑向均勻氣流，逐漸加速后再進(jìn)入平行氣膜間隙。穩(wěn)流管3、10的入口直徑d2為20 mm，高度h3為70 mm；供氣孔4、11的直徑d1為10 mm，圓盤直徑D為120 mm，圓盤厚度h2為30 mm。

圖1 高壓圓盤氣體軸承結(jié)構(gòu)示意

文中采用Solidworks專業(yè)建模軟件對(duì)高壓圓盤氣體軸承進(jìn)行三維建模，建模時(shí)忽略穩(wěn)流管與工作圓盤之間的密封墊及穩(wěn)流管端面倒角，以及工作圓盤上、下端面的倒角，并在氣膜間隙出口建立封閉面，方便抽取流體域。簡(jiǎn)化后的高壓圓盤氣體軸承3D幾何模型如圖2所示。

圖2 高壓圓盤氣體軸承的簡(jiǎn)化三維結(jié)構(gòu)

1.2 軸承邊界命名

為了在保證計(jì)算精度的同時(shí)，盡量減少網(wǎng)格數(shù)量，以減少計(jì)算時(shí)間，文中采用氣體軸承三維模型的1/2作為計(jì)算域。在DesignModeler建模工具中導(dǎo)入Solidworks生成后綴名為x_t的模型文件，使用Fill功能對(duì)軸承內(nèi)流域進(jìn)行抽取，并利用布爾運(yùn)算功能將穩(wěn)流管和工作圓盤合并成一體。其流體域與固體域的邊界命名如圖3所示。

圖3 氣體軸承流體與固體計(jì)算域的邊界命名

2 網(wǎng)格劃分與邊界條件設(shè)置

2.1 網(wǎng)格劃分

在Fluid Flow(Fluent)模塊中，使用Mesh網(wǎng)格劃分工具，對(duì)高壓圓盤氣體軸承的流體域與固體域進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分，其最終劃分網(wǎng)格結(jié)果如圖4所示。

圖4 三維軸承模型的固體域網(wǎng)格與圓盤收縮段附近的流體域局部網(wǎng)格

流體域網(wǎng)格劃分方法選擇Multizone，采用尺寸函數(shù)類型為Curvature，網(wǎng)格平滑度為中等，曲率法向角度為3°。對(duì)流體域近壁區(qū)添加5層膨脹層，膨脹層設(shè)置方法選擇為Smooth Transition，增長(zhǎng)率設(shè)置為1.05，固體域網(wǎng)格劃分方法選擇Body sizing，且采用與流體域相同的單元尺寸。

2.2 單元基本尺寸的選擇

單元基本尺寸是控制模型網(wǎng)格數(shù)量的主要參數(shù)，基本尺寸越小，網(wǎng)格數(shù)量越大；單元基本尺寸的大小也會(huì)影響到網(wǎng)格質(zhì)量，進(jìn)而對(duì)計(jì)算結(jié)果產(chǎn)生影響。詳細(xì)分析不同單元基本尺寸算例的網(wǎng)格質(zhì)量，對(duì)選擇合適的單元尺寸作為最終分析算例的設(shè)置參數(shù)具有指導(dǎo)作用。單元基本尺寸從1 mm變化到0.42 mm時(shí)，5種算例對(duì)應(yīng)的網(wǎng)格數(shù)量和質(zhì)量參數(shù)如表1所示。

表1 不同單元基本尺寸對(duì)應(yīng)的網(wǎng)格數(shù)量和質(zhì)量

從表1可知，單元基本尺寸從1 mm變化到0.42 mm時(shí)，網(wǎng)格數(shù)量從100萬增大到340萬；網(wǎng)格的平均質(zhì)量、長(zhǎng)寬比、雅可比值均變化較小，但歪斜率和扭曲因子變化較大；歪斜度從0.118 49減小到0.061 04，扭曲因子從2.378 1×10-4減小到1.121 6×10-4，兩者都幾乎減小了1/2。因此，表1中算例4和算例5的劃分參數(shù)是較為合適的參數(shù)。圖5進(jìn)一步比較了5個(gè)算例在上圓盤流固耦合面(流體側(cè))上的熱流密度。

圖5 不同網(wǎng)格數(shù)量時(shí)耦合面的熱流密度

由圖5可知，在算例4和算例5中，算例4的上圓盤流固耦合面(流體側(cè))上的熱流密度計(jì)算值最大。對(duì)于具有相同導(dǎo)熱系數(shù)的同種材料，熱流密度越大，對(duì)應(yīng)的溫度梯度也越大，圓盤內(nèi)部溫度分布不均勻的程度也越大，進(jìn)而其熱變形也越嚴(yán)重。因此，采用算例4的參數(shù)進(jìn)行共軛傳熱求解，對(duì)后續(xù)的軸承圓盤變形分析、材料選擇和幾何參數(shù)設(shè)計(jì)而言，是保守的方案，得到的結(jié)果具有更大的安全裕量。故文中采用290萬網(wǎng)格的模型進(jìn)行共軛傳熱求解。

2.3 邊界條件設(shè)置

如圖3所示，將上進(jìn)氣口、下進(jìn)氣口均設(shè)置為壓力入口，設(shè)置參數(shù)保持一致，即：入口總壓為2.03 MPa，入口初始?jí)毫?.02 MPa，湍流強(qiáng)度為5%，水力直徑為20 mm，入口總溫為300.15 K；將氣膜出口設(shè)置為壓力出口，出口壓力為0.10 MPa，湍流強(qiáng)度為5%，水力直徑為120 mm，出口總溫為300.15 K；將流體域中的上、下圓盤流固耦合面(流體側(cè))與固體域中的上、下圓盤流固耦合面(固體側(cè))設(shè)置為固定無滑移耦合面；將流體計(jì)算域?qū)ΨQ面、上圓盤對(duì)稱面、下圓盤對(duì)稱面均設(shè)置為對(duì)稱面；將上穩(wěn)流管外側(cè)面、下穩(wěn)流管外側(cè)面、上圓盤外表面、下圓盤外表面設(shè)置為與空氣自然對(duì)流的壁面，對(duì)流換熱系數(shù)取10 W/(m2·K)，自由流溫度取300.15 K；考慮到穩(wěn)流管將與管接頭連接，故將上、下穩(wěn)流管端面設(shè)置為絕熱壁面。

設(shè)置固體材料為Steel，其屬性保持默認(rèn)；設(shè)置流體材料為Air，其密度采用理想氣體狀態(tài)方程，黏度采用三參數(shù)的薩瑟蘭公式，比熱容和導(dǎo)熱系數(shù)均采用溫度校正；設(shè)置操作壓力為0；采用SIMPLEC算法，其各項(xiàng)差分格式保持默認(rèn)。能量松弛因子取0.3，其他松弛因子保持默認(rèn)。

3 導(dǎo)熱方程與湍流模型

3.1 導(dǎo)熱方程

高壓圓盤氣體軸承工作時(shí)，一方面，流道內(nèi)的氣流與圓盤壁面之間存在溫差，會(huì)產(chǎn)生對(duì)流換熱；另一方面，由于軸承圓盤內(nèi)部存在著非均勻溫度場(chǎng)，固體圓盤內(nèi)部也會(huì)產(chǎn)生熱傳導(dǎo)。流體計(jì)算域的能量方程[26]為

(1)

高壓圓盤氣體軸承固體計(jì)算域的導(dǎo)熱方程如下：

(2)

3.2 湍流模型

文中使用Realizablek-ε湍流模型，它是帶旋流修正的k-ε模型，對(duì)比標(biāo)準(zhǔn)的k-ε模型，該模型增加了一個(gè)與湍流黏性有關(guān)的公式，為湍流耗散率增加了新的傳輸方程。其優(yōu)點(diǎn)之一是對(duì)平板和圓柱射流的發(fā)散比率有著更精確的預(yù)測(cè)。Realizablek-ε湍流模型的方程[26]如下：

Gb-ρε-YM

(3)

(4)

3.3 共軛傳熱問題的計(jì)算流程

Fluid Flow(Fluent)模塊中共軛傳熱問題的計(jì)算流程與單一流動(dòng)問題類似，區(qū)別在于計(jì)算域同時(shí)包括流體域和固體域，設(shè)置邊界條件時(shí)，既要設(shè)置流體域的邊界條件，也要設(shè)置固體域的邊界條件。計(jì)算時(shí)對(duì)流體域和固體域進(jìn)行整體迭代求解，得到耦合面的參數(shù)。計(jì)算流程如圖6所示。

圖6 共軛傳熱問題的計(jì)算流程

4 數(shù)值結(jié)果及分析

4.1 總體速度場(chǎng)和壓力場(chǎng)

圖7所示為高壓圓盤氣體軸承間隙流場(chǎng)的速度流線圖及圓盤收縮段附近的流場(chǎng)速度矢量圖。由圖7(a)可知，氣體從上、下穩(wěn)流管入口流入，流經(jīng)穩(wěn)流管收縮段、供氣孔、圓盤收縮段及平行氣膜區(qū)，最終在氣膜間隙出口處流出至環(huán)境大氣。穩(wěn)流管內(nèi)的氣流速度普遍較低，而平行氣膜區(qū)氣流速度普遍較高，氣流從進(jìn)入穩(wěn)流管開始逐漸加速，并在工作圓盤收縮段出口處達(dá)到極大值，而后在平行氣膜區(qū)速度有所降低，在氣膜間隙出口處，氣流速度再次加速達(dá)到最大值。盡管進(jìn)入上、下進(jìn)氣口的氣流都存在一定的湍流度，但在上、下穩(wěn)流管軸對(duì)稱收縮段的整流作用下，氣流形成了均勻的流線從上、下供氣孔進(jìn)入流道間隙。由圖7(b)可以直觀地觀察到，流道中心區(qū)域的氣流速度矢量分布及圓盤收縮段壁面邊界層的速度梯度變化，此處氣流經(jīng)供氣孔流入，在上、下圓盤收縮段的作用下，將軸向來流無分離地轉(zhuǎn)變?yōu)閺较蛄鲃?dòng)，并沿半徑R方向逐漸加速，最后從氣膜間隙出口流出。

圖7 氣體軸承間隙流場(chǎng)的速度流線圖及圓盤收縮段附近的速度矢量分布

使用Fluent軟件自帶的Iso-surface截線工具可以截取高壓圓盤氣體軸承的氣膜對(duì)稱線，并提取對(duì)稱線上的數(shù)據(jù)。如圖8所示為高壓圓盤氣體軸承氣膜對(duì)稱線上的馬赫數(shù)和靜壓分布。

圖8 氣膜對(duì)稱線上的馬赫數(shù)和靜壓分布

由圖8中的馬赫數(shù)曲線可知：氣流馬赫數(shù)在供氣孔區(qū)域a從接近于0逐漸增大至0.11，并在收縮段區(qū)域b快速增大至超音速1.13；而后又在平行氣膜區(qū)c先緩慢下降，平行氣膜區(qū)c的馬赫數(shù)整體維持在0.86左右的高亞音速狀態(tài)；最后在接近氣膜間隙出口時(shí)由于壓力遠(yuǎn)大于環(huán)境背壓，從而使馬赫數(shù)快速增加到超音速1.3。

由圖8中的靜壓曲線可知：氣流壓力在供氣孔區(qū)域a緩慢下降至1.95 MPa，而后在圓盤收縮段區(qū)域b快速下降至0.827 MPa，在到達(dá)平行氣膜區(qū)c后，氣流壓力經(jīng)過小幅上升后再次連續(xù)下降至氣膜出口處的0.349 MPa。因?yàn)槌羲倭鲃?dòng)的基本特性是下游的低壓擾動(dòng)不會(huì)影響上游的壓力分布，故外界環(huán)境的低壓擾動(dòng)不會(huì)影響軸承流道內(nèi)的壓力分布，所以在高壓圓盤氣體軸承流道內(nèi)的絕大部分區(qū)域，壓力能維持在0.592 MPa以上，較好地實(shí)現(xiàn)了高壓氣體潤(rùn)滑。

4.2 傳熱分析

4.2.1 溫度場(chǎng)分析

如圖9所示為軸承對(duì)稱面上的溫度云圖，軸承對(duì)稱面由流體計(jì)算域?qū)ΨQ面、上圓盤對(duì)稱面、下圓盤對(duì)稱面組合而成。

圖9 軸承對(duì)稱面上的溫度云圖

在來流總溫不變的條件下，馬赫數(shù)越大，溫度越低。結(jié)合圖8、9可知，供氣孔區(qū)域氣流速度較低，該區(qū)域的溫度與環(huán)境溫度相比變化很小，故上、下穩(wěn)流管的溫度變化不明顯；圓盤收縮段出口處及氣膜間隙出口處的馬赫數(shù)較高，相應(yīng)的溫度也較低；而在平行氣膜區(qū)絕大部分區(qū)域氣流仍然保持在0.86馬赫以上的高亞音速狀態(tài)，故整個(gè)軸承圓盤在鄰近平行氣膜區(qū)的部分溫度較低，而遠(yuǎn)離平行氣膜區(qū)的部分溫度較高。

圖10所示為與軸承對(duì)稱面相垂直的截面上的溫度等值線圖。

圖10 與軸承對(duì)稱面相垂直的截面上溫度等值線

由圖10可知，上、下圓盤的溫度以氣膜對(duì)稱面為對(duì)稱平面呈對(duì)稱式分布，圓盤收縮段出口處及氣膜間隙出口處的溫度最低，這兩處周圍的等溫線較為密集，溫度梯度較大，導(dǎo)熱熱流密度也較大。軸承圓盤內(nèi)的溫度梯度矢量以收縮段出口處及氣膜間隙出口處為中心，向圓盤外端面及穩(wěn)流管方向逐漸減小。熱流密度矢量也相應(yīng)減小，其方向與溫度梯度矢量的方向相反。

氣膜對(duì)稱線及流體域耦合壁面上的溫度分布如圖11所示。

圖11 氣膜對(duì)稱線及流體域耦合壁面上的溫度分布

由圖11可知，在氣膜對(duì)稱線上，氣流在供氣孔區(qū)域a的溫度變化并不明顯，而在圓盤收縮段區(qū)域b的前半部分變化較小，后半部分變化較大，收縮段出口處的溫度已陡降至241 K。而后在平行氣膜區(qū)c由于氣流的速度有所下降，其溫度有明顯上升，在平行氣膜區(qū)c的絕大部分區(qū)域，溫度都維持在260 K以上。當(dāng)氣流到達(dá)氣膜間隙出口處時(shí)，由于速度再次增大至超音速，其溫度驟降至235 K。圖11中流體域耦合壁面上的溫度變化趨勢(shì)與氣膜對(duì)稱線上氣流溫度的分布趨勢(shì)相類似，但耦合壁面在工作圓盤收縮段出口處及氣膜出口處的溫度明顯高于氣膜對(duì)稱線上相應(yīng)位置的氣流溫度，說明高亞音速特別是超音速狀態(tài)下圓盤壁面溫度邊界層中可以存在較高的溫度差。總的說來，共軛傳熱條件下，流固耦合壁面上的溫度分布仍然取決于間隙內(nèi)的氣流速度分布，即軸承初始的速度型設(shè)計(jì)。

4.2.2 熱流密度分析

由于耦合的流體域壁面與固體域壁面上的熱流密度大小相等方向相反，故僅分析流體域耦合壁面的熱流密度即可。如圖12所示為高壓圓盤氣體軸承流體域耦合壁面上的熱流密度分布。

圖12 流體域耦合壁面上的熱流密度分布

由圖12可知，圓盤收縮段區(qū)域b的流體域耦合壁面上熱流密度先小幅下降后再上升，然后再經(jīng)過平緩的下降后再明顯上升，并在收縮段區(qū)域b的出口處達(dá)到峰值，其值為73 467.1 W/m2，且收縮段區(qū)域 b的部分區(qū)域流體壁面熱流密度為負(fù)，說明在收縮段區(qū)域b內(nèi)的不同位置，同時(shí)存在流體吸熱和放熱的現(xiàn)象，這是由于軸承圓盤的熱傳導(dǎo)與間隙內(nèi)氣體對(duì)流換熱具有的復(fù)雜共軛作用機(jī)制造成的。平行氣膜區(qū)c的換熱達(dá)到熱充分發(fā)展?fàn)顟B(tài)，氣流吸熱，熱流密度為正。氣膜間隙出口處的流體域壁面熱流密度達(dá)到最大值135 966 W/m2。

4.3 共軛傳熱與非共軛傳熱的對(duì)比分析

為了研究高壓圓盤氣體軸承共軛傳熱解與非共軛傳熱解的差異性，以圖3(a)為研究模型，獲取高壓圓盤氣體軸承在恒溫壁面條件下的計(jì)算結(jié)果，此時(shí)整個(gè)計(jì)算域只包含單一的流體域，可看作是一類非共軛傳熱解。即先在Fluid Flow(Fluent)模塊的Mesh界面中，將單一流體域模型輸出成后綴名為msh的網(wǎng)格文件，然后單獨(dú)啟動(dòng)Fluent 2020R2軟件，導(dǎo)入前述的msh網(wǎng)格文件；將流體域上、下圓盤壁面設(shè)置為固定無滑移恒溫壁面，壁溫300.15 K，其余設(shè)置同2.3節(jié)。

將單一流體域的計(jì)算結(jié)果與共軛傳熱條件下的計(jì)算結(jié)果進(jìn)行對(duì)比，圖13所示為氣膜對(duì)稱線上的馬赫數(shù)、靜壓和靜溫對(duì)比。可知，共軛傳熱條件下高壓圓盤氣體軸承氣膜對(duì)稱線上的馬赫數(shù)、靜壓、靜溫和非共軛傳熱恒溫壁面條件下的變化趨勢(shì)是一致的，但在局部區(qū)域數(shù)值結(jié)果存在明顯差異。由圖13(a)可知，在供氣孔區(qū)域a和圓盤收縮段區(qū)域b的絕大部分區(qū)域，兩者僅存在微小差異，但從收縮段出口至平行氣膜區(qū)c的絕大部分區(qū)域，兩者存在明顯的差異；其中最大差異發(fā)生在圓盤收縮段出口處，兩者相差達(dá)到10.76%。

圖13 氣膜對(duì)稱線上的流場(chǎng)參數(shù)對(duì)比

由圖13(b)可知，共軛傳熱和非共軛傳熱恒溫壁面條件下軸承氣膜對(duì)稱線上的靜壓在供氣孔區(qū)域a、圓盤收縮段區(qū)域b和平行氣膜區(qū)c均存在較為明顯的差異，且最大差異發(fā)生在收縮段出口處，兩者相差15.78%。由圖13(c)可知，共軛傳熱和非共軛傳熱恒溫壁面條件下氣膜對(duì)稱線上的靜溫在供氣孔區(qū)域a和圓盤收縮段區(qū)域b的絕大部分區(qū)域存在微小差異，但在收縮段出口處及平行氣膜區(qū)c均存在較大差異，其中收縮段出口處溫度差為13.05 K，相對(duì)誤差為5.14%；最大差異發(fā)生在氣膜間隙出口處，溫度差為18.14 K，相對(duì)誤差為7.16%。綜上，造成2種條件下流場(chǎng)參數(shù)計(jì)算結(jié)果存在明顯差異的原因，在于2種情況下流固耦合壁面上的溫度分布不同，進(jìn)而使軸承間隙流道內(nèi)，尤其是平行氣膜區(qū)的馬赫數(shù)發(fā)生變化，進(jìn)而改變了間隙內(nèi)的壓力和溫度分布。

如圖14所示為高壓圓盤氣體軸承共軛傳熱與非共軛傳熱恒溫壁面條件下流體域耦合壁面上的熱流密度的對(duì)比。可知，2種情況下流體域耦合壁面上熱流密度的變化趨勢(shì)相類似，但圓盤收縮段b及平行氣膜區(qū)c的數(shù)值結(jié)果均存在較大差異，其中氣膜間隙出口處的熱流密度差異為60.06%，最大差異發(fā)生在收縮段出口處，兩者相差79.28%。非共軛傳熱恒溫壁面條件下，間隙內(nèi)的氣體只吸熱，流體域耦合壁面上的熱流密度均為正值。

圖14 流體域耦合壁面的熱流密度對(duì)比

共軛傳熱和非共軛傳熱恒溫壁面條件下，流固耦合壁面上熱流密度計(jì)算結(jié)果產(chǎn)生顯著變化的主要原因，在于恒溫壁面假設(shè)在很大程度上高估了圓盤壁面氣膜側(cè)的溫度。由于在流固耦合面附近，軸承圓盤的熱傳導(dǎo)與氣膜內(nèi)氣體的對(duì)流換熱發(fā)生了復(fù)雜的共軛作用機(jī)制，這與恒溫壁面簡(jiǎn)化下的傳熱機(jī)制有本質(zhì)的不同。軸承在實(shí)際工作中，流固耦合壁面上的溫度不可能保持恒溫，因此恒溫壁面假設(shè)使計(jì)算誤差偏大，共軛傳熱計(jì)算可以得到更為符合實(shí)際的結(jié)果[23]。

5 結(jié)論

(1)由于存在著復(fù)雜的共軛作用機(jī)制，高壓圓盤氣體軸承流固耦合壁面上的熱流密度及溫度分布曲線，均是隨R的增大發(fā)生劇烈變化的曲線；流固耦合壁面上既非恒熱流密度，也非恒壁溫，這是非共軛傳熱視角下難以獲得的結(jié)果；高壓圓盤氣體軸承的傳熱問題，是一個(gè)典型的共軛傳熱問題。

(2)共軛傳熱時(shí)，高壓圓盤氣體軸承流固耦合壁面上溫度分布曲線的變化趨勢(shì)，與氣膜對(duì)稱線上溫度變化趨勢(shì)類似，但二者在圓盤收縮段出口處和氣膜間隙出口處具體溫度值有明顯差異；因此，總的說來，流固耦合壁面上的溫度分布仍然取決于間隙內(nèi)氣流速度分布，即軸承初始的速度型設(shè)計(jì)。

(3)共軛傳熱與非共軛傳熱恒溫壁面條件下得到的氣膜對(duì)稱線上的馬赫數(shù)、靜壓、靜溫以及流固耦合壁面上的熱流密度、溫度均存在明顯差異。非共軛傳熱恒溫壁面條件下，間隙內(nèi)的氣體只吸熱，流體域耦合壁面上熱流密度均為正值；而共軛傳熱條件下流體域耦合壁面上熱流密度存在正負(fù)值，即間隙內(nèi)氣體在流動(dòng)過程中，吸熱和放熱同時(shí)存在，顯示出軸承圓盤的熱傳導(dǎo)與間隙內(nèi)氣體的對(duì)流換熱具有復(fù)雜的共軛作用機(jī)制。

(4)分析高壓圓盤氣體軸承的溫度場(chǎng)，推薦采用共軛傳熱方法進(jìn)行計(jì)算，可以得到更為符合實(shí)際的結(jié)果。