波度和錐度對(duì)液體潤(rùn)滑機(jī)械密封空化特性影響

李振濤，郝木明，楊文靜，韓婕，任寶杰

（中國(guó)石油大學(xué)（華東）密封技術(shù)研究所，山東 青島 266580）

?

波度和錐度對(duì)液體潤(rùn)滑機(jī)械密封空化特性影響

李振濤，郝木明，楊文靜，韓婕，任寶杰

（中國(guó)石油大學(xué)（華東）密封技術(shù)研究所，山東 青島 266580）

為進(jìn)一步探究機(jī)械密封液膜中空穴初生、演變等影響因素，基于質(zhì)量守恒空化模型，建立考慮微觀波度和錐度的液體潤(rùn)滑機(jī)械密封數(shù)學(xué)模型，采用有限控制體積法對(duì)控制方程進(jìn)行離散，綜合分析不同波幅、波數(shù)及錐度對(duì)空化特性的影響。結(jié)果表明：波度產(chǎn)生流體動(dòng)壓潤(rùn)滑效應(yīng)，錐度產(chǎn)生流體靜壓潤(rùn)滑效應(yīng)，兩者共同影響空穴的初生和演變，而又影響空穴沿周向和徑向的變化；空穴初生時(shí)，液膜破裂位置位于膜厚沿周向發(fā)散的低壓區(qū)，且在空穴演變時(shí)其沿某一半徑方向保持不變；量綱1波幅不超過(guò)0.5且波數(shù)不低于8或量綱1負(fù)錐度大于0.5的工況均可加快空化率，促進(jìn)空化發(fā)生；而相同波幅且波數(shù)小于8或錐度趨向于正錐度的工況均可減緩空化率，有效抑制空化發(fā)生。

機(jī)械密封；液體潤(rùn)滑；波度；錐度；空化；有限控制體積法

DOI：10.11949/j.issn.0438-1157.20151733

引　言

因加工精度、壓力變形或熱變形等因素導(dǎo)致產(chǎn)生的表面形貌（如表面粗糙度、周向波度及徑向錐度等）或各種型式激光加工微造型的客觀存在，液膜潤(rùn)滑機(jī)械密封端面間產(chǎn)生大量的微觀或宏觀的間隙發(fā)散區(qū)和收斂區(qū)；而端面間液膜在發(fā)散區(qū)易產(chǎn)生壓力下降，當(dāng)壓力降至低于密封介質(zhì)的飽和蒸汽壓時(shí)，液膜開(kāi)始發(fā)生汽化，并形成汽泡和空穴，被稱(chēng)為空化現(xiàn)象［1-2］?？栈F(xiàn)象對(duì)機(jī)械密封端面間液膜壓力的分布、液膜承載能力及密封穩(wěn)定性具有較大影響［3］。

Findley［4］于1968年最早提出機(jī)械密封中的空化問(wèn)題，隨后引起了國(guó)內(nèi)外學(xué)者的重視并進(jìn)行了多年的試驗(yàn)研究和理論研究。Nau［5］于1980年通過(guò)試驗(yàn)觀測(cè)到機(jī)械密封中的空化現(xiàn)象，并用分析法確定空化區(qū)的邊界；Khonsari等［6-7］對(duì)微造型端面密封進(jìn)行可視化觀察試驗(yàn)，觀察到空化發(fā)生位置并且分析認(rèn)為液膜厚度會(huì)影響空化初生，繼而在研究空化效應(yīng)對(duì)液膜承載力的研究中僅考慮了表面粗糙度和微造型結(jié)構(gòu)尺寸等影響因素；孟凡明等［8］研究了機(jī)械密封或軸承端面微造型的幅值、分布位置及長(zhǎng)寬比等結(jié)構(gòu)參數(shù)對(duì)空穴機(jī)理的影響；Dowson等［9］介紹了求解液膜潤(rùn)滑控制方程的空化邊界條件，即Sommerfeld、Half-Sommerfeld、Reynolds和JFO邊界；其中，Sommerfeld邊界為無(wú)空化邊界，Half-Sommerfeld邊界將計(jì)算結(jié)果中出現(xiàn)負(fù)壓部分置為零，適用于第一次求近似解，而Reynolds邊界應(yīng)用最廣但無(wú)法正確解釋液膜再形成的過(guò)程，不滿(mǎn)足質(zhì)量守恒；Jakobsson等［10-11］提出的JFO邊界很好地彌補(bǔ)了Reynolds邊界存在的缺陷，滿(mǎn)足質(zhì)量守恒。胡紀(jì)濱等［12］采用質(zhì)量守恒空化邊界條件，研究了不同工況和槽深對(duì)徑向直線(xiàn)槽端面密封的空化演變進(jìn)行了數(shù)值模擬；郝木明等［13］采用數(shù)值計(jì)算方法研究了空化效應(yīng)對(duì)螺旋槽液膜密封的穩(wěn)態(tài)特性的影響。總之，目前在機(jī)械密封空化研究領(lǐng)域，絕大多數(shù)主要集中于對(duì)該現(xiàn)象的試驗(yàn)觀察、端面微造型結(jié)構(gòu)參數(shù)對(duì)空穴影響及空穴對(duì)密封潤(rùn)滑性能影響的分析；但在誘發(fā)空化發(fā)生影響因素領(lǐng)域，尤其是密封環(huán)端面微觀形貌對(duì)其影響，未引起足夠重視，研究不夠深入。

鑒于上述存在的問(wèn)題，本文基于質(zhì)量守恒空化模型，建立考慮微觀周向波度和徑向錐度的液膜潤(rùn)滑機(jī)械密封數(shù)學(xué)模型，采用有限體積法對(duì)液膜潤(rùn)滑控制方程進(jìn)行離散并通過(guò)Gauss-Seidel松弛迭代方法進(jìn)行求解，以期獲得微觀波度和錐度對(duì)密封端面液膜中空穴初生、空穴區(qū)域演變及空化率等影響規(guī)律，為機(jī)械密封端面形貌參數(shù)優(yōu)化設(shè)計(jì)提供參考。

1　數(shù)學(xué)模型與數(shù)值求解方法

1.1數(shù)學(xué)模型

機(jī)械密封典型結(jié)構(gòu)如圖1所示，當(dāng)旋轉(zhuǎn)環(huán)和靜止環(huán)間形成的密封面處于純液體潤(rùn)滑狀態(tài)［14］時(shí)，彈簧壓力和密封介質(zhì)壓力形成的閉合力與密封面間的液膜產(chǎn)生的開(kāi)啟力相平衡。

圖1　機(jī)械密封結(jié)構(gòu)和潤(rùn)滑模型Fig.1　Schematics of mechanical seal and lubrication model

假設(shè)在配對(duì)的旋轉(zhuǎn)環(huán)和靜止環(huán)中，一個(gè)密封面為光滑平面，另一個(gè)密封面具有各向同性的周向波度和徑向錐度，密封面幾何模型如圖2所示。

對(duì)于圖2所示的密封面間液膜總厚度定義為

式中，密封面內(nèi)徑處最小膜厚作為參考膜厚（用符號(hào)refh表示），大小為1 μm。而錐度正切定義為

對(duì)式（1）進(jìn)行量綱1處理，得

圖2　密封面幾何模型Fig.2　Geometric model of sealing interface

圖2所示幾何模型結(jié)構(gòu)參數(shù)與操作工況參數(shù)及文獻(xiàn)［15］中相關(guān)參數(shù)分別如表1中Value 1和Value 2［15］所述，其中，本文計(jì)算的模擬介質(zhì)為22#液壓油。

表1　R密封面幾何結(jié)構(gòu)參數(shù)與工況參數(shù)Table 1　Sealing interface geometrical parameters and operational parameters

1.2控制方程

為確定密封端面間液膜控制方程及簡(jiǎn)化數(shù)值計(jì)算，做如下假設(shè)：

（1）液膜厚度為微米級(jí)，足夠小，適用于雷諾方程；

（2）密封面間流體為牛頓流體，工作狀態(tài)為穩(wěn)態(tài)層流且黏度不隨壓力變化；

（3）密封面間由兩相混合組成的薄膜分為全液膜區(qū)和空穴區(qū)，全液膜區(qū)流體為不可壓縮，空穴區(qū)壓力保持不變；

（4）密封流體的熱楔及配對(duì)密封環(huán)的熱變形忽略不計(jì)，且忽略除離心力外的慣性項(xiàng)的影響；

（5）該模型適用于等溫絕熱工況。

鑒于上述假設(shè)，液膜潤(rùn)滑控制方程在極坐標(biāo)下的表達(dá)式［16-17］為

對(duì)于機(jī)械密封空化，目前主要存在兩種空化模型：傳統(tǒng)的Reynolds模型和JFO模型，前者僅給出了液膜破裂位置的邊界條件，但無(wú)法正確描述液膜再形成位置的邊界條件，不滿(mǎn)足質(zhì)量守恒定律；而后者不僅給出了液膜破裂位置的邊界條件

同時(shí)，也給出了液膜再形成位置的邊界條件

該模型在數(shù)值求解過(guò)程中可實(shí)現(xiàn)跟蹤動(dòng)態(tài)液膜，能很好地滿(mǎn)足質(zhì)量守恒定律。因此本文采用滿(mǎn)足質(zhì)量守恒的JFO空化模型。

在全液膜區(qū)，液膜密度cρ保持不變，則式（3）簡(jiǎn)化為

在空穴區(qū)，液膜壓力為空化壓力cp，滿(mǎn)足，則式（3）簡(jiǎn)化為

為實(shí)現(xiàn)JFO空化模型在數(shù)值計(jì)算過(guò)程中自動(dòng)捕捉液膜破裂和再形成邊界及將全液膜區(qū)和空穴區(qū)的控制方程統(tǒng)一化，引入通用變量φ和空化數(shù)F，兩參數(shù)滿(mǎn)足

具體定義為

全液膜區(qū)：

空穴區(qū)：

對(duì)式（3）進(jìn)行量綱1化，可得到同時(shí)描述全液膜區(qū)和空穴區(qū)壓力分布的通用控制方程式中，量綱1參量為

1.3邊界條件設(shè)置

圖1所示機(jī)械密封配對(duì)的密封環(huán)關(guān)于中心對(duì)稱(chēng)，數(shù)值計(jì)算時(shí)選取一個(gè)波度為計(jì)算周期。量綱1化之前式（3）的初始邊界條件為

周期性邊界條件為）

1.4數(shù)值求解方法

1.4.1計(jì)算域網(wǎng)格劃分圖3所示計(jì)算域局部網(wǎng)格，圖中剖面線(xiàn)為以節(jié)點(diǎn)P為中心的控制體積。

圖3　極坐標(biāo)下控制體積示意圖Fig.3　Control volume in polar coordinates

為使模擬計(jì)算結(jié)果更接近于實(shí)際，計(jì)算域網(wǎng)格劃分的數(shù)目應(yīng)遵守隨著網(wǎng)格數(shù)目的增加計(jì)算結(jié)果不再有顯著變化的原則。因此，在滿(mǎn)足網(wǎng)格密度的條件下，應(yīng)當(dāng)盡量減少計(jì)算量，提高計(jì)算收斂的穩(wěn)定性和收斂速度。由于密封面間液膜中發(fā)生空穴的位置及區(qū)域未確定，控制體積網(wǎng)格劃分采用結(jié)構(gòu)化網(wǎng)格。

1.4.2控制方程離散化為獲得迭代代數(shù)方程，采用有限控制體積法對(duì)控制方程進(jìn)行離散［20］。首先，在控制體積內(nèi)對(duì)式（11）進(jìn)行積分，表達(dá)式為

在控制體積中，沿we→，角度為θΔ；沿sn→，距離為rΔ，則式（15）的離散式整理為

式（16）等號(hào)左邊的壓差流項(xiàng)采用中心差分格式，而等號(hào)右邊的剪切流項(xiàng)和離心慣性項(xiàng)采用迎風(fēng)差分格式［21］。則中心節(jié)點(diǎn)P處φ值的代數(shù)方程表達(dá)式如

式中，

1.4.3代數(shù)方程迭代方法式（17）中的PS項(xiàng)含有φ和F的二次項(xiàng)，致使方程具有非線(xiàn)性［22］；因此，在迭代過(guò)程中，將φ和F的上一次迭代值代替當(dāng)前的迭代值。為進(jìn)一步提高迭代速度并保證良好的收斂性，采用Gauss-Seidel迭代方法進(jìn)行求解，并引入松弛因子ξ和η，對(duì)每一步計(jì)算的φ和F進(jìn)行迭代更新，迭代更新方程為

式中，

代數(shù)方程迭代收斂判別準(zhǔn)則為

2　計(jì)算結(jié)果與分析

2.1程序算法驗(yàn)證與算例

以文獻(xiàn)［15］中的結(jié)構(gòu)參數(shù)和工況參數(shù)（表1）為例，采用本文數(shù)值計(jì)算方法進(jìn)行計(jì)算，將計(jì)算得到的半個(gè)波度周期內(nèi)膜壓分布（表2）與文獻(xiàn)［15］中的膜壓分布數(shù)據(jù)表（表3）進(jìn)行對(duì)比。

表2　R采用本文數(shù)值計(jì)算方法求解的計(jì)算周期內(nèi)膜壓分布數(shù)據(jù)Table 2　Film pressure distribution datasheet calculated by this paper's numerical method/MPa

表3　R文獻(xiàn)［15］中的膜壓分布數(shù)據(jù)Table 3　Film pressure distribution datasheet of Ref. ［15］/MPa

對(duì)比表2和表3可知，采用本文數(shù)值計(jì)算方法得到的計(jì)算域壓力值、空化率與文獻(xiàn)［15］中的計(jì)算結(jié)果誤差分別為2%～4%、2.3%，空化區(qū)域邊界沿半徑方向向外徑側(cè)偏移一個(gè)網(wǎng)格距離。由于本文算法精度較高，計(jì)算誤差在允許范圍內(nèi)，可認(rèn)為兩者計(jì)算結(jié)果十分吻合，從而證明本文算法的準(zhǔn)確性，可很好地開(kāi)展本文的后續(xù)計(jì)算。

以表1中的Values 1參數(shù)為算例，并根據(jù)表4所述的計(jì)算域網(wǎng)格無(wú)關(guān)性檢測(cè)結(jié)果，確定在單周期計(jì)算域內(nèi)沿徑向和周向網(wǎng)格劃分為50×50，計(jì)算得到的單周期內(nèi)膜壓分布如圖4所示。

表4　R網(wǎng)格劃分?jǐn)?shù)目對(duì)液膜空化率的影響Table 4　Effects of grid numbers on liquid film cavitation rate

圖4　算例量綱1膜壓分布Fig.4　Dimensionless film pressure distribution of example

由圖4可知，在液膜收斂區(qū)域、發(fā)散區(qū)域分別形成高壓區(qū)和低壓區(qū)。液膜經(jīng)當(dāng)前高壓區(qū)進(jìn)入低壓區(qū)后，由于局部壓力過(guò)低，液膜無(wú)法承受致使開(kāi)始破裂；經(jīng)低壓區(qū)后進(jìn)入下一個(gè)計(jì)算周期的高壓區(qū)，隨著壓力的升高，破裂液膜逐漸聚集再形成，因此在液膜破裂位置和再形成位置的區(qū)間形成空穴區(qū)，即圖4中顏色較深位置。

2.2對(duì)空穴初生及區(qū)域演變影響

2.2.1波幅的影響波幅是波度尺寸表征參數(shù)之一。取量綱1波幅A分別為0.05、0.10、0.65和1.50，量綱1錐度表征量B分別為1和0.4時(shí)，量綱1半處沿周向的液膜壓力及液膜空穴區(qū)分布如圖5、圖6和圖7所示。

圖5　不同A時(shí)，B=1、=1.07處的量綱1膜壓分布Fig.5　Effects of parameter A on dimensionless film pressure distribution at B=1=1.07

圖6　不同A時(shí)，B=0.4、=1.07處的量綱1膜壓分布Fig.6　Effects of parameter A on dimensionless film pressuredistribution at B=0.4=1.07

由圖5可知：不考慮錐度時(shí)即B=1，因波度的存在，液膜沿周向收斂間隙產(chǎn)生高壓區(qū)，沿發(fā)散間隙產(chǎn)生低壓區(qū)。隨著波幅的增加，流體動(dòng)壓潤(rùn)滑效應(yīng)逐漸增強(qiáng)，高壓區(qū)最高壓力顯著增加，而低壓區(qū)最低壓力明顯降低。當(dāng)A=1.5時(shí)，由于低壓區(qū)最低壓力仍高于空化壓力（量綱1空化壓力為0），未產(chǎn)生空穴。因此，液膜壓力發(fā)生顯著變化，主要是波度產(chǎn)生流體動(dòng)壓潤(rùn)滑效應(yīng)作用的結(jié)果。

考慮錐度后（量綱1負(fù)錐度值為0.6），如圖6所示，端面液膜產(chǎn)生流體靜壓潤(rùn)滑效應(yīng)，在Sommerfeld無(wú)空化邊界下，高壓區(qū)和低壓區(qū)的最高/最低壓力的數(shù)值大小相比不考慮錐度時(shí)均增加。因采用JFO質(zhì)量守恒空化邊界，當(dāng)A低于0.05時(shí)未產(chǎn)生空穴，增加至0.05～0.1區(qū)間的某值時(shí)即將發(fā)生空化；當(dāng)A為0.1時(shí)低壓區(qū)局部壓力已低于空化壓力，產(chǎn)生空穴。

圖7　B=0.4、不同A時(shí)液膜空穴分布（空白處為空穴區(qū)）Fig.7　Effects of parameter A on cavitation region distribution at B=0.4 （the blank is cavitation region）

隨著波幅的增加，流體動(dòng)壓潤(rùn)滑效應(yīng)逐漸增強(qiáng)，空穴邊界沿徑向向兩側(cè)擴(kuò)展，而沿周向單側(cè)擴(kuò)展，空化區(qū)域增大，如圖7所示。當(dāng)A達(dá)到0.65時(shí)，空化率達(dá)到最大值；繼續(xù)增大至1.5，液膜破裂位置幾乎不變，位于周向第31個(gè)節(jié)點(diǎn)所在半徑線(xiàn)上；但因流體動(dòng)壓潤(rùn)滑效應(yīng)的減弱，空化區(qū)域緩慢降低，液膜再形成位置沿周向向破裂位置收縮。

2.2.2波數(shù)的影響波數(shù)也是波度尺寸表征參數(shù)之一。取波數(shù)K分別為4、12、24和48，量綱1波幅A和量綱1錐度表征量B分別為0.7和0.4時(shí)，量綱1半徑=1.07處沿周向的液膜壓力及液膜空穴區(qū)分布如圖8和圖9所示。

由8圖可知：高壓區(qū)最高壓力基本位于同一位置未發(fā)生偏移，在Sommerfeld無(wú)空化邊界下，低壓區(qū)最低壓力分布位置相類(lèi)似。因采用JFO質(zhì)量守恒空化邊界，當(dāng)K=4時(shí)，低壓區(qū)已初生空穴。

隨著波數(shù)的增加，流體動(dòng)壓潤(rùn)滑效應(yīng)增強(qiáng)，液膜破裂位置形成所謂“負(fù)壓”的能力也隨之增強(qiáng)。由圖9可知，液膜開(kāi)始破裂位置基本不變且沿周向第31個(gè)節(jié)點(diǎn)所在半徑方向向兩側(cè)延伸，液膜再形成位置沿著周向單側(cè)外移，空化區(qū)域雖逐漸擴(kuò)大，但形狀基本保持不變。當(dāng)K達(dá)到24時(shí)，液膜空化率最大；隨著波數(shù)的進(jìn)一步增加，流體動(dòng)壓潤(rùn)滑效應(yīng)非增卻降，致使空化率降低；當(dāng)波數(shù)增加一倍即K=48時(shí)，空化區(qū)域降低約2.5%，變化較??；但空化區(qū)域形狀呈現(xiàn)狹長(zhǎng)變化趨勢(shì)。

圖8　不同K時(shí)，A=0.7、B=0.4、=1.07處的膜壓分布Fig.8　Effects of parameter K on dimensionless film pressure distribution at A=0.7， B=0.4，=1.07

圖9　A=0.7、B=0.4、不同K時(shí)液膜空穴分布Fig.9　Effects of parameter K on cavitation region distribution at A=0.7， B=0.4 （the blank is cavitation region）

2.2.3錐度的影響取量綱1錐度表征量B分別為0.20、0.55、0.75和1.10，量綱1波度A分別為0和0.7時(shí)，量綱1半徑=1.07處沿周向的液膜壓力及液膜空穴區(qū)分布如圖10、圖11和圖12所示。

由圖10可知，不考慮波度、量綱1錐度為負(fù)錐度時(shí)，密封面間隙沿徑向發(fā)散，進(jìn)入間隙內(nèi)的液體在最小膜厚處發(fā)生節(jié)流膨脹，壓力急劇降低，液膜壓力呈凹形拋物線(xiàn)變化；量綱1錐度為正錐度時(shí)，密封面間隙沿徑向收斂，進(jìn)入間隙的液體受阻，壓力緩慢降低，液膜壓力呈凸形拋物線(xiàn)變化。當(dāng)錐度由負(fù)錐度增加至正錐度，液膜壓力由凹拋物線(xiàn)漸變至凸拋物線(xiàn)，主要是流體靜壓潤(rùn)滑效應(yīng)作用的結(jié)果。

圖10　不同B時(shí)，A=0、=1.07處的量綱1膜壓分布Fig.10　Effects of parameter B on dimensionless film pressure distribution at A=0，=1.07

圖11　不同B時(shí)，A=0.7、=1.07處量綱1膜壓分布Fig.11　Effects of parameter B on dimensionless film pressure distribution at A=0.7=1.07

考慮波度后，隨著錐度由負(fù)錐度向正錐度增加，逐漸增強(qiáng)的流體靜壓潤(rùn)滑效應(yīng)較波度產(chǎn)生的流體動(dòng)壓潤(rùn)滑效應(yīng)占主導(dǎo)作用，促進(jìn)低壓區(qū)壓力呈增大趨勢(shì)，液膜破裂位置形成所謂“負(fù)壓”的能力卻隨之減弱，空化區(qū)域逐漸較小。同樣，在空化區(qū)域演化過(guò)程中，其形狀基本保持不變。

2.3對(duì)空化率影響

繪制密封面間隙液膜中空化區(qū)域所占的比率即空化率隨周向波度和徑向錐度的變化曲線(xiàn)，分別如圖13、圖14及圖15所示。

圖12　不同B時(shí)液膜空穴分布Fig.12　Effects of parameter B on cavitation region distribution（the blank is cavitation region）

圖13　A=0～2， 不同B時(shí)空化率的變化Fig.13　Cavitation rate variation at A=0—2 and different B

由圖13和圖14可知：當(dāng)量綱1波幅A為低波幅0～0.5，量綱1錐度（B-1）＜0時(shí)，空化率迅速增加，尤其是（B-1）≤-0.5時(shí)，顯著促進(jìn)了空穴的初生及空化區(qū)域演變；當(dāng)量綱1錐度趨向于正錐度時(shí)，空化率明顯降低至零，進(jìn)而抑制了空化的初生及空化區(qū)域演變。當(dāng)A＞0.5后，空化率呈緩慢下降趨勢(shì)。

由圖15可知：當(dāng)量綱1波幅A為低波幅0～0.5，波數(shù)為高頻波即K≥8時(shí)，空化率迅速增加，亦顯著促進(jìn)了空穴的初生及空化區(qū)域演變；當(dāng)波數(shù)為低頻波即K＜8，尤其是K=4時(shí)，未發(fā)生空化，有效抑制了空穴的初生及空化區(qū)域演變。

圖14　B=0～2， 不同A時(shí)空化率的變化Fig.14　Cavitation rate variation at B=0—2 and different A

圖15　A=0～2，不同K時(shí)空化率的變化Fig.15　Cavitation rate variation at A=0—2 and different K

3　結(jié)　論

（1）通過(guò)與文獻(xiàn)［15］計(jì)算結(jié)果對(duì)比，兩者十分吻合，本文采用有限控制體積法離散統(tǒng)一描述全液膜區(qū)和空穴區(qū)壓力分布的通用控制方程的方法是可行的。

（2）液膜潤(rùn)滑機(jī)械密封端面空穴初生及空化區(qū)域演變是密封面微觀波度和錐度共同作用的結(jié)果，而空穴初生即液膜開(kāi)始破裂位置位于膜厚沿周向發(fā)散的低壓區(qū)，且在空化區(qū)域演變時(shí)其沿某一半徑方向保持不變，空化區(qū)域形狀亦基本保持不變。

（3）在低波幅、高頻波數(shù)或負(fù)錐度較大時(shí)，液膜空化率迅速增加，促進(jìn)了空化發(fā)生；在低波幅、低頻波數(shù)或錐度由負(fù)錐度趨向于正錐度時(shí)，液膜空化率明顯降低甚至降為零，有效抑制了空化發(fā)生。

（4）本文考慮了端面微觀波度和錐度對(duì)空穴初生及其演變等影響，但未提及微觀粗糙度和角偏差等各因素影響；因此，后續(xù)將進(jìn)一步開(kāi)展相關(guān)研究。

符號(hào)說(shuō)明

A——量綱1波度波幅

Am——波度波幅，m

B——量綱1錐度表征量

F——空化數(shù)

h——密封面間隙總膜厚，m

href——參考膜厚，m

hid，min——密封面內(nèi)徑處最小膜厚，m

hod，min——密封面外經(jīng)處最小膜厚，m

K——波度波數(shù)

Pf——量綱1液膜壓力

p——密封面間隙液膜壓力

pi——密封面內(nèi)徑處壓力，Pa

po——密封面外徑處壓力，Pa

pc——空化壓力，Pa

Re?——密封面間隙流體Reynolds數(shù)

R——量綱1外半徑

r——密封面半徑，m

ri——密封面內(nèi)半徑，m

ro——密封面外半徑，m

Vn——外法線(xiàn)方向速度，m·s-1

α——密封面錐度，rad

φ——量綱1通用變量

γ——密封面間隙流體過(guò)程系數(shù)

η——空化數(shù)迭代松弛因子

μ——密封面間隙流體動(dòng)力黏度，Pa·s

θ——密封面圓周弧度，rad

θi——單周期波度初始弧度，rad

ρ——空化區(qū)氣液混相平均密度，kg·m-3

ρc——密封面間隙流體密度，kg·m-3

ω——密封面旋轉(zhuǎn)角速度，rad·s-1

ξ——量綱1壓力迭代松弛因子

下角標(biāo)

E——控制容積東側(cè)相鄰節(jié)點(diǎn)

e——控制容積東側(cè)邊界

i——密封面圓周方向網(wǎng)格節(jié)點(diǎn)下標(biāo)

imax——密封面圓周方向網(wǎng)格節(jié)點(diǎn)總數(shù)

j——密封面半徑方向網(wǎng)格節(jié)點(diǎn)下標(biāo)

jmax——密封面半徑方向網(wǎng)格節(jié)點(diǎn)總數(shù)

N——控制容積北側(cè)相鄰節(jié)點(diǎn)

n——控制容積北側(cè)邊界

new ——當(dāng)前迭代過(guò)程變量值

old ——上一次迭代過(guò)程變量值

P ——控制容積中心節(jié)點(diǎn)

S ——控制容積南側(cè)相鄰節(jié)點(diǎn)

s ——控制容積南側(cè)邊界

W ——控制容積西側(cè)相鄰節(jié)點(diǎn)

w ——控制容積西側(cè)邊界

References

［1］ 王濤， 黃偉峰， 王玉明. 機(jī)械密封液膜汽化問(wèn)題研究現(xiàn)狀及進(jìn)展［J］. 化工學(xué)報(bào)， 2012， 63 （11）： 3375-3382. DOI： 10.3969/j.issn.0438-1157.2012.11.001. WANG T， HUANG W F， WANG Y M. Research and progress of mechanical seals operating with vaporization transition ［J］. CIESC Journal， 2012， 63 （11）： 3375-3382. DOI： 10.3969/j.issn.0438-1157.2012.11.001.

［2］陳匯龍， 吳強(qiáng)波， 左木子， 等. 機(jī)械密封端面液膜空化的研究進(jìn)展［J］. 排灌機(jī)械工程學(xué)報(bào)， 2015， 33 （2）： 138-144. DOI： 10.3969/ j.issn.1674-8530.14.0085. CHEN H L， WU Q B， ZUO M Z， et al. Overview on liquid film cavitation in mechanical seal faces ［J］. Journal of Drainage and Irrigation Machinery Engineering， 2015， 33 （2）： 138-144. DOI：10.3969/j.issn.1674-8530.14.0085.

［3］唐飛翔， 孟祥凱， 李紀(jì)云， 等. 基于質(zhì)量守恒的Laserface液體潤(rùn)滑機(jī)械密封數(shù)值分析 ［J］. 化工學(xué)報(bào)， 2013， 64 （10）： 3694-3700. DOI： 10.3969/j.issn.0438-1157.2013.10.029. TANG F X， MENG X K， LI J Y， et al. Numerical analysis of Laser-Face liquid mechanical seal based on mass conservation ［J］. CIESC Journal， 2013， 64 （10）： 3694-3700. DOI： 10.3969/j.issn.0438-1157.2013.10.029.

［4］FINDLAY J A. Cavitation in mechanical face seals ［J］. Journal of Lubrication Technology， 1968， 90 （2）： 356-364.

［5］NAU B S. Observations and analysis of mechanical seal film characteristics ［J］. Journal of Lubrication Technology. 1980， 102 （3）：341-347.

［6］QIU Y， KHONSARI M M. On the prediction of cavitation in dimples using a mass-conservative algorithm ［J］. Journal of Tribology， 2009，131 （4）： 41702-41711.

［7］QIU Y， KHONSARI M M. Experimental investigation of tribological performance of laser textured stainless steel rings ［J］. Tribology International， 2011， 44 （5）： 635-644.

［8］孟凡明， 王戰(zhàn)江. 微造型表面上氣穴發(fā)生機(jī)理初步研究 ［J］. 摩擦學(xué)學(xué)報(bào)， 2011， 31 （2）： 137-143. DOI： 10.16078/j.tribology.2011.02.014. MENG F M， WANG Z J. Preliminary study on mechanism for occurrence of cavitation on a textured surface ［J］. Tribology， 2011， 31（2）： 137-143. DOI： 10.16078/j.tribology.2011.02.014.

［9］DOWSON D， TAYLOR C M. Cavitation in bearings， annual reviews［J］. Fluid Mechanics， 1979， （11）： 35-66.

［10］JAKOBSSON B， FLOBERG L. The finite journal bearing，considering vaporization ［J］. Wear， 1958， 2 （2）： 85-88.

［11］OLSSON K O. Cavitation in dynamically loaded bearings ［J］. Wear，1967， 55 （2）： 295-304.

［12］ 胡紀(jì)濱， 劉丁華， 魏超. 徑向直線(xiàn)槽端面密封空化特性數(shù)值模擬 ［J］.摩擦學(xué)學(xué)報(bào)， 2011， 31 （6）： 551-556. DOI： 10.16078/j.tribology.2011.06.007. HU J B， LIU D H， WEI C. Numerical simulation for cavitation of groove face seals ［J］. Tribology， 2011， 31 （6）： 551-556. DOI：10.16078/j.tribology.2011.06.007.

［13］郝木明， 莊媛， 章大海， 等. 考慮空化效應(yīng)的螺旋槽液膜密封特性數(shù)值研究 ［J］. 中國(guó)石油大學(xué)學(xué)報(bào) （自然科學(xué)版）， 2015， 39 （3）：132-137 . DOI： 10.3969/j.issn.1673-5005.2015.03.018. HAO M M， ZHUANG Y， ZHANG D H， et al. Numerical study on sealing performance of spiral groove liquid film seal considering effects of cavitation ［J］. Journal of China University of Petroleum（Natural Science）， 2015， 39 （3）： 132-137. DOI： 10.3969/j.issn.1673-5005.2015.03.018.

［14］LUBBINGE H. On the lubrication of mechanical face seals ［D］. Enschede： University of Twente， 1999.

［15］顧永泉. 機(jī)械密封的空化和空化邊界確定 ［J］. 流體機(jī)械， 1998， 26（12）： 16-20. GU Y Q. Determination of cavitation and cavitation boundary of mechanical seal ［J］. Fluid Machinery， 1998， 26 （12）： 16-20.

［16］PINKUS O， LUND J W. Centrifugal effects in thrust bearings and seals under laminar conditions ［J］. Transactions of the ASME， 1981，103： 126-136.

［17］YU T H， SADEGHI F. Groove effects on thrust washer lubrication ［J］. Journal of Tribology， 2001， 123 （2）： 295-304.

［18］ 李強(qiáng)， 鄭水英， 劉淑蓮. 計(jì)入JFO邊界條件的滑動(dòng)軸承性能分析［J］. 機(jī)械強(qiáng)度， 2010， 32 （2）： 270-274. LI Q， ZHENG S Y， LIU S L. Analysis of the performance of journal bearings with JFO boundary condition ［J］. Journal of Mechanical Strength， 2010， 32 （2）： 270-274.

［19］PAYVAR P， SALANT R F. A computational method for cavitation in a wavy mechanical seal ［J］. Journal of Tribology， 1992， 114： 199-204.

［20］PATANKAR S V. Numerical Heat Transfer and Fluid Flow ［M］. London： Taylor & Francis， 1980： 72-73.

［21］XIONG S W， WANG J Q. Steady-state hydrodynamic lubrication modeled with the Payvar-Salant mass conservation model ［J］. Journal of Tribology， 2012， 134： 031703-1-031703-16.

［22］趙一民， 苑士華， 胡紀(jì)濱， 等. 基于質(zhì)量守恒邊界條件的螺旋槽旋轉(zhuǎn)密封性能分析 ［J］. 機(jī)械工程學(xué)報(bào)， 2014， 50 （22）： 142-149. DOI：10.3901/JME.2014.22.142. ZHAO Y M， YUAN S H， HU J B， et al. Performance analysis of spiral-groove rotary seals considering mass conserving boundary condition ［J］. Journal of Mechanical Engineering， 2014， 50 （22）：142-149. DOI： 10.3901/JME.2014.22.142.

Effects of waviness and taper on cavitation characteristics of liquid lubricated mechanical seals

LI Zhentao， HAO Muming， YANG Wenjing， HAN Jie， REN Baojie
（Institute of Sealing Technology， China University of Petroleum， Qingdao 266580， Shandong， China）

In order to further probe into the influence factors of cavitation start and evolution， the mathematical model of liquid film lubricated mechanical seals considering the micro circumferential waviness and radial taper was established based on the mass-conservative cavitation model. Liquid film governing equation was discretized by the finite control volume method and was solved by the Gauss-Seidel relaxation iterative algorithm. The effects of different waviness amplitudes， waviness numbers and tapers on the cavitation characteristics were analyzed. The results indicated that the start and evolution of the cavitation were affacted by the combined action of the hydrodynamic lubrication effect produced by the circumferential waviness and the hydrostatic lubrication effect which was greatly influenced by radial taper. To a large extent， changing of the cavitation area along the circumferential direction was influenced by the waviness， while that along the radial direction depended on the taper. The rupture position of liquid was located where film thickness diverged in the depression region at the start stage of cavitation， while the subsequent rupture direction remained unchanged at the evolution stage of cavitation. The cavitation rate was dramatically accelerated at the situations of dimensionless waviness amplitude not beyond 0.5 and waviness number greater than or equal 8， or dimensionless negative taper value greater than 0.5， which promoted the occurrence of cavitation. But under the same waviness amplitude and waviness number less than 8or the taper tending to be positive， the cavitation rate was retarded， which effectively inhibited the occurrence of cavitation.

date： 2015-11-18.

Prof. HAO Muming， haomm@upc.edu.cn

supported by the National Natural Science Foundation of China （51375497）， and the Shandong Special Projects of Independent Innovation and Achievement Transformation （2014ZZCX10102-4）..

mechanical seals； liquid lubricated； waviness； taper； cavitation； finite control volume method

TH 117.2

A

0438—1157（2016）05—2005—10

2015-11-18收到初稿，2016-01-29收到修改稿。

聯(lián)系人：郝木明。第一作者：李振濤（1983—），男，博士研究生。

國(guó)家自然科學(xué)基金項(xiàng)目（51375497）；山東省自主創(chuàng)新及成果轉(zhuǎn)化專(zhuān)項(xiàng)（2014ZZCX10102-4）。