余旻豐 彭旭東 孟祥鎧 梁楊楊

摘要：沖洗是強化機械密封換熱的主要措施之一，但在一些特殊場合下沖洗量的大小往往是給定的，因此，需要采取一些強制換熱措施來改善機械密封的運行環(huán)境。通過在密封環(huán)的外周面開設(shè)織構(gòu)，針對沖洗量一定的情況，基于SST k-ω湍流模型，采用Ω方法分析了不同轉(zhuǎn)速下織構(gòu)深徑比對端面溫度、外周面局部努塞爾數(shù)Nu和織構(gòu)區(qū)域流場的影響，對比研究了動環(huán)織構(gòu)和靜環(huán)織構(gòu)的換熱機理。研究結(jié)果表明：在相同工況、沖洗量和織構(gòu)幾何參數(shù)條件下，若動靜環(huán)外周面單獨開設(shè)織構(gòu)，則動環(huán)的換熱效果更佳。在相同轉(zhuǎn)速下，動環(huán)外周面開設(shè)織構(gòu)時，減小深徑比會使織構(gòu)內(nèi)部換熱效果差的區(qū)域增大，換熱效果減弱，但對于靜環(huán)織構(gòu)，減小深徑比會使織構(gòu)流體流動下游側(cè)換熱效果差的區(qū)域減小，換熱效果增強；隨著轉(zhuǎn)速的增大，較小深徑比動環(huán)織構(gòu)內(nèi)部出現(xiàn)了新的差換熱效果區(qū)域，單位面積換熱強度顯著下降。因此，為獲取較優(yōu)的換熱效果，在低轉(zhuǎn)速下宜選取小深徑比織構(gòu)，而在高轉(zhuǎn)速下應(yīng)選取大深徑比織構(gòu)。

關(guān)鍵詞：表面織構(gòu); 外周面; 接觸式機械密封; 強化換熱

中圖分類號：TH117.1

DOI：10.3969/j.issn.1004132X.2023.11.002

Research on Heat Transfer Enhancement Mechanism of Contact Mechanical

Seals with Textured Circumference Surfaces

YU Minfeng PENG Xudong MENG Xiangkai LIANG Yangyang

School of Mechanical Engineering，Zhejiang University of Technology，Hangzhou，310014

Abstract： Flushing was one of the main measures for heat transfer enhancement of mechanical seals. But in some special occasions， the flow rate of flushing was given， so some forced heat transfer measures were needed to improve the operating environment of mechanical seals. Through processing texture on the circumference surfaces of the sealing rings and aiming at the condition of a given amount of flushing， the influences of speed and depth to diameter ratio on interface temperature， local Nu on circumference surfaces and flow field in textured regions were analyzed by SST k-ω turbulence model and Ω method， heat transfer mechanism of texture processed on rotor or stator was compared and analyzed. The results show that under the same working conditions， flow rate of flushing and texture geometric parameters， the heat transfer of rotor is better if the texture is set on the circumference surfaces of rotor and stator separately. Under the same speed， when texture is processed on the circumference surface of rotor， reduction of the aspect ratio may increase the areas with poor heat transfer inside the texture and weakened heat transfer. But when texture is processed on the circumference surfaces of stator， the conclusion is completely opposite. As the speed increases， a new area of poor heat transfer appears in a small aspect ratio rotor texture， the heat transfer intensity per unit area is decreased significantly. Therefore， in order to obtain better heat transfer effectiveness， the small aspect ratio texture may be selected under low speed， and the large aspect ratio texture may be selected under high speed.

Key words： texture； circumference surface； contact mechanical seal； heat transfer enhancement

0 引言

機械密封是旋轉(zhuǎn)機器設(shè)備的重要零部件，廣泛應(yīng)用于泵、壓縮機、釜等設(shè)備的軸端密封［1-3］。通常，要求機械密封運行時端面間處于合適的流體潤滑狀態(tài)，以確保泄漏量小甚至零泄漏且端面磨損低，但是，在很多情況下，機械密封泄漏超標(biāo)的前期預(yù)兆或原因并非表現(xiàn)為端面磨損，而是端面過熱或溫升過高或摩擦扭矩過大，特別是針對航天渦輪泵等設(shè)備用機械密封其沖洗量無法調(diào)節(jié)的特殊使用場合，后一種原因更為常見［4-8］。因此，如何在沖洗量給定的情況下，提高機械密封的傳熱特別是換熱效果，改善機械密封的運行環(huán)境已成為航天渦輪泵的難點問題。近期，密封環(huán)外周面開設(shè)織構(gòu)的方法由于兼顧了較好的強化換熱效果和較小的攪拌損失而受到關(guān)注［9-11］。

關(guān)于織構(gòu)表面換熱機制的研究最初是以平板表面為研究對象。AFANASYEV等［12］在平板表面設(shè)置了半球形織構(gòu)，實驗測量了不同深度、不同半徑、不同排布織構(gòu)內(nèi)部的溫度分布和速度分布，確定了半球形織構(gòu)的強化換熱效果，這一結(jié)果被LIN等［13］通過數(shù)值模擬得到了驗證。CHYU等［14］同樣采用平板表面進行實驗，通過測量半球形和水滴形織構(gòu)的內(nèi)部溫度分布，發(fā)現(xiàn)盡管二者的局部溫度分布存在微小差異，但計算所得整體努塞爾數(shù)Nu基本相等。

MAHMOOD等［15］使用兩塊相隔一定間距的矩形平板組成的狹窄通道來模擬微通道換熱器，并在通道內(nèi)壁面加工交錯分布的規(guī)則半球形織構(gòu)，綜合采用粒子圖像測速（particle image velocimetry， PIV） 技術(shù)和壁面測溫技術(shù)，獲得了半球形織構(gòu)內(nèi)部流場和局部努塞爾數(shù)，分析了Re=600～11 000時量綱一高度H/D（其中H為平板間距，D為半球形織構(gòu)直徑）、入口溫度對通道內(nèi)壁面半球形織構(gòu)換熱效果的影響［16］，研究結(jié)果表明，隨著H/D的減小和入口溫度的降低，織構(gòu)的換熱效果增強。

在上述基礎(chǔ)上，多位學(xué)者［17-20］分別對圓柱形織構(gòu)、傾斜側(cè)壁織構(gòu)、六邊形織構(gòu)和橢圓形織構(gòu)進行研究，對比了這些織構(gòu)與半球形織構(gòu)的換熱效果。LIGRANI等［21］在通道的上下壁面設(shè)置對稱的凸起與織構(gòu)，考慮錯位因素的影響，對流場、局部努塞爾數(shù)和流動摩阻進行了測試，結(jié)果表明凸起表面強化換熱的效果更佳，但阻力增大了2～2.7倍。LIU等［22］在半球形織構(gòu)通道的內(nèi)壁面加入半球形凸起，并考慮了凸起和織構(gòu)的幾何和排布參數(shù)對換熱效果的影響，同樣發(fā)現(xiàn)凸起能有效提高下游織構(gòu)的強化換熱效果，但流動阻力也相應(yīng)增大了。XIE等［23］將方形直通道改為方形U形通道，以模擬葉片尖端流道，并在通道內(nèi)壁半球形織構(gòu)之間加工了圓柱形凸起，研究結(jié)果表明，不同部位的局部努塞爾數(shù)提高了3.2%～31.5%，流動阻力同時也增大了2.3% ～16.2%。

NIAN等［9］率先提出在機械密封外圓周表面加工織構(gòu)并進行了相關(guān)的實驗，結(jié)果證實這一技術(shù)在中低轉(zhuǎn)速和中低壓情況下能降低端面溫度10%左右；隨后他們又研究了中低速條件下外圓周圓柱形織構(gòu)的排布和幾何參數(shù)對換熱強化的影響［10］，發(fā)現(xiàn)小尺寸、高密度的織構(gòu)換熱效果更佳。周宇坤等［11］針對多種織構(gòu)的換熱效果開展了初步研究，結(jié)果表明旋轉(zhuǎn)角為90°的等邊三角形織構(gòu)具有最佳強化換熱效果，矩形織構(gòu)強化換熱效果最差。

為了更高效地對流場進行量化分析，研究者們采用了Ω渦識別方法。截至目前，對渦的分析方法已經(jīng)歷了三代［24-25］：①基于渦量的渦識別方法，即直接使用渦量定義法識別渦，認(rèn)為渦量的大小即為當(dāng)?shù)匦D(zhuǎn)運動的強度，但已經(jīng)被證明存在大誤差，且高渦量區(qū)域和實際的渦結(jié)構(gòu)關(guān)聯(lián)性很低；②第二代渦識別方法，其基本方法是設(shè)置一個判定參數(shù)和一個閾值，當(dāng)計算出的判定參數(shù)數(shù)值大于閾值時，認(rèn)為當(dāng)?shù)赜袦u，但閾值的選擇將極大程度地影響識別結(jié)果且判定參數(shù)物理意義不明；③第三代渦識別方法，包括Ω渦識別方法和Liutex向量法，其優(yōu)點是將參數(shù)的閾值歸一化至0.52，使渦中心位置和強度的判定更為精準(zhǔn)。

上述工作表明，在眾多換熱手段中，織構(gòu)能兼顧強換熱效果和低摩阻損失，同時無需占用額外空間，在空間受限的航天渦輪泵密封腔體內(nèi)是一種理想的強化換熱措施。但是，截止目前，絕大多數(shù)研究還是以無限大平板和微通道換熱器為研究對象，這與傳統(tǒng)密封腔內(nèi)部的流動與傳熱的情況存在明顯不同。此外，織構(gòu)應(yīng)用于機械密封的研究多以端面開設(shè)織構(gòu)為主［26-28］，而在沖洗量給定的條件下，高速機械密封外圓周面織構(gòu)的換熱機理鮮有研究。因此，本文采用SST湍流模型對比分析了密封動環(huán)或靜環(huán)外周面織構(gòu)的內(nèi)部流場和局部努塞爾數(shù)，使用Ω方法對密封腔內(nèi)部和織構(gòu)局部的漩渦進行了識別，討論了不同深徑比外周面織構(gòu)動環(huán)和靜環(huán)的換熱效果及其變化，研究了轉(zhuǎn)速對動環(huán)織構(gòu)換熱效果的影響，揭示了織構(gòu)的換熱機理。

1 數(shù)值模型

1.1 幾何模型

圖1為一種典型的高速渦輪泵用機械密封的結(jié)構(gòu)示意圖，其中，機械密封動環(huán)或靜環(huán)的外周面加工有織構(gòu)。為提高計算效率，考慮到動環(huán)和靜環(huán)沿圓周方向的對稱性，截取環(huán)整體的1/18作為計算域，建立包含動環(huán)、靜環(huán)、密封腔內(nèi)流體的計算模型，相關(guān)幾何參數(shù)如圖1、圖2和表1所示。

1.2 流動傳熱控制方程

本文重點分析研究在不同織構(gòu)深徑比和轉(zhuǎn)速條件下機械密封端面溫度變化及密封腔和織構(gòu)內(nèi)部的流動狀態(tài)。為便于分析，可作如下假設(shè)：流體在接觸面上無滑移；忽略密封端面的表面粗糙度與波度、密封環(huán)的熱力變形和端面泄漏；不考慮介質(zhì)物性變化及熱輻射；由于接觸式機械密封泄漏量極小，故忽略端面處介質(zhì)泄漏對溫度場的影響。考慮到密封環(huán)在腔體內(nèi)旋轉(zhuǎn)，腔內(nèi)流體容易處于湍流狀態(tài)，本文采用商業(yè)軟件ANSYS FLUENT 2020R2內(nèi)的SST k-ω湍流模型（下文簡稱SST模型）來計算腔內(nèi)流體的壓力場、速度場，并結(jié)合能量方程來模擬計算傳熱及溫度場。已有的研究［9，11］使用了RNG k-ε湍流模型配合壁面函數(shù)作為流動傳熱控制方程，但這種研究方法對壁面量綱一參數(shù)y+的依賴性較大［29-31］。其中，量綱一參數(shù)y+的定義以及RNG k-ε模型中y+應(yīng)滿足的取值范圍如下［31］：

式中，cμ為經(jīng)驗常數(shù)，取cμ=0.09；k為湍動能；ν為流體運動黏度；y為流體域內(nèi)第一個節(jié)點與固體壁面間的距離。

本文涉及織構(gòu)內(nèi)部速度場分析，若使用RNG k-ε模型配合壁面函數(shù)方法，則由于第一層網(wǎng)格不能過密，織構(gòu)內(nèi)小漩渦將會難以分辨，而SST模型通過混合函數(shù)F的不同取值，分別在近壁面區(qū)域和湍流旺盛區(qū)域使用k-ω模型和k-ε模型，這樣就能同時確保壁面換熱和織構(gòu)內(nèi)部流動的計算精度，因此，本文使用SST模型，該模型y+僅需要滿足近壁區(qū)湍流模型（即k-ω模型）的要求，確保y+≤1即可［32-34］。

FLUENT中SST k-ω模型的穩(wěn)態(tài)不可壓縮流表達式如下：

式中，ρ為介質(zhì)密度；k為湍動能；ω為湍流耗散率；Γk、Γω分別為k和ω的有效擴散項；Gk、Gω分別為k和ω的方程；Yk、Yω分別為k和ω的發(fā)散項；Dω為正交發(fā)散項；Sk、Sω為用戶自定義源項；xi、xj為各坐標(biāo)方向的張量表示，下標(biāo)i、j為啞標(biāo)；ui為液相速度矢量的分量。

式（2）和式（3）中各項除Sk與Sω以外均由FLUENT默認(rèn)給出，本文不使用源項。其他計算中所涉及的模型常數(shù)的取值見表2，其中φ1、φ2為參數(shù)φ在計算過程中使用的中間值，φ=φ1F+φ2 （1-F），F(xiàn)為混合函數(shù)，F(xiàn)在近壁區(qū)取1，在湍流旺盛區(qū)取0；α為邊界層內(nèi)部剪切應(yīng)力與湍動能之比；β為k方程的耗散項經(jīng)驗常數(shù);σk、σω分別為k方程和ω方程中的擴散項經(jīng)驗常數(shù)。

1.3 渦分析方法

Ω方法定義了一個Ω：

其中，e是一個任意取值很小的正數(shù)，A是渦量中非旋轉(zhuǎn)（變形）的部分，B是渦量中旋轉(zhuǎn)部分，V表示對速度矢量進行三維哈密頓運算。這一方法推薦選取Ω≥0.52的等值面來展示渦，這一推薦閾值意味著渦量的旋轉(zhuǎn)部分超過了變形耗散部分，且這一數(shù)值適用于大多數(shù)經(jīng)驗案例［24-25］。本文使用Tecplot的自定義方程計算密封腔內(nèi)各部位流場的Ω，并對內(nèi)部渦分布進行了分析。

圖3為主軸轉(zhuǎn)速為5400 r/min時，分別使用Ω方法和渦量分析法獲得的密封腔內(nèi)部軸對稱平面上的漩渦強度云圖，Ω和渦量的數(shù)值越大表明漩渦強度越高?？梢钥闯?，由于動環(huán)的旋轉(zhuǎn)，動環(huán)外圓周近壁區(qū)的速度梯度遠(yuǎn)大于其余區(qū)域的速度梯度，因此采用渦量分析法得到在動環(huán)外圓周近壁區(qū)出現(xiàn)了更大的渦量數(shù)值，而Ω方法則識別出了密封腔內(nèi)部存在兩

個旋轉(zhuǎn)方向相反的主要漩渦（下文簡稱主渦）。對于圖1所示的幾何結(jié)構(gòu)，密封腔可以近似視為一種內(nèi)壁面旋轉(zhuǎn)的有限長同心圓環(huán)腔體內(nèi)部的泰勒庫艾特流［35］，顯而易見，Ω方法識別出的主渦形態(tài)符合泰勒渦的形態(tài)，而渦量分析法則完全無法識別，因此，Ω方法相比于渦量分析法更適合本文設(shè)計研究對象的流動傳熱分析。

1.4 邊界條件

根據(jù)實際使用中高速渦輪泵機械密封的安裝特點，由圖1可知，密封腔內(nèi)流體的入口設(shè)置于動環(huán)背部，出口設(shè)置于靜環(huán)外側(cè)。常用沖洗流量一般為8～25 L/min，這里為安全起見，選取下限值8 L/min作為沖洗流量的計算用值，因此計算得到入口流速約為0.1 m/s。端面的摩擦熱參考NIAN等［9］提出的方法進行計算并加載于圖2所示的WR和WS面上，本文所使用摩擦熱公式中的關(guān)鍵參數(shù)為端面摩擦因數(shù)f，由摩擦狀態(tài)決定，軟硬配對的接觸式機械密封啟動時端面摩擦因數(shù)較大，穩(wěn)定運行后端面摩擦因數(shù)將會減小，一般取值范圍在0.1～0.2之間，極少數(shù)能達到0.25［36］，這一范圍適用于混合摩擦、邊界摩擦及干摩擦，又考慮到高轉(zhuǎn)速會導(dǎo)致摩擦因數(shù)增大，本文選取f=0.15，這既不影響中低轉(zhuǎn)速時的產(chǎn)熱計算，也考慮了高轉(zhuǎn)速對端面產(chǎn)熱模型的影響。本研究認(rèn)為動、靜環(huán)內(nèi)徑處與大氣之間存在熱對流，根據(jù)經(jīng)驗公式［37］可計算出對流傳熱系數(shù)分別為72 W/（m2·K）和0.1 W/（m2·K），最后，設(shè)置WR與WS為耦合壁面， 這樣FLUENT就會根據(jù)產(chǎn)熱量、溫度場及材料物性自動將端面摩擦熱分配至動環(huán)與靜環(huán)。具體邊界條件列于表3，表中T為壁面溫度，v為入口流速，p為出口壓力，q為壁面熱通量，h為對流傳熱系數(shù)，n為主軸轉(zhuǎn)速。各邊界編號見圖2，密封環(huán)材料和密封介質(zhì)的物性參數(shù)如表4所示。

1.5 參數(shù)定義

1.6 網(wǎng)格無關(guān)性與模型正確性驗證

采用ICEM軟件對所截取密封環(huán)和腔內(nèi)流體進行網(wǎng)格劃分，環(huán)與流體均采用完全六面體網(wǎng)格劃分，對網(wǎng)格進行規(guī)則加密。圖4a為網(wǎng)格局部加密總示意圖，圖4b為動、靜環(huán)外周加密示意圖，圖4c為端面加密示意圖，圖4d為織構(gòu)底面加密示意圖，圖4e為織構(gòu)側(cè)壁面加密示意圖，局部加密規(guī)則為：①確保壁面y+≤1；②邊界層至少鋪15層網(wǎng)格；③邊界層內(nèi)部網(wǎng)格尺寸增長倍率為1.2～1.4，以接近1.2為宜。

網(wǎng)格無關(guān)性驗證和模型正確性驗證如圖5所示，其中Z為本研究模型中密封環(huán)軸向坐標(biāo)，z為文獻［41］所給的密封環(huán)軸向坐標(biāo)，零值為端面所在平面處軸向坐標(biāo)。圖5a所示為網(wǎng)格無關(guān)性驗證結(jié)果，選取網(wǎng)格數(shù)N分別為5×105、1×106、2×106的網(wǎng)格進行網(wǎng)格無關(guān)性驗證，以端面徑向溫度分布為衡量標(biāo)準(zhǔn)。

如圖5a所示，當(dāng)網(wǎng)格數(shù)增加至1×106時，再增大一倍網(wǎng)格尺寸，端面溫度曲線各點數(shù)值變化小于0.5 ℃，相對誤差在0.7%以內(nèi)，綜合考慮計算時間和計算精度，選取1×106數(shù)量網(wǎng)格作為計算網(wǎng)格。圖5b所示為模型正確性驗證，可以看出，在對比文獻［41］所給參數(shù)條件下，使用本文模型計算所得的三個不同軸向位置的徑向溫度分布與對比文獻實驗值的變化規(guī)律均吻合較好，各點溫度偏差均小于0.5 ℃，相對誤差均在1.5%以內(nèi)，優(yōu)于對比文獻的數(shù)值模擬結(jié)果。究其原因，本文主要改進了所使用的湍流模型，RNG k-ε湍流模型要求壁面y+較大，而SST模型僅需加密邊界層確保y+≤1即可較好地符合實驗數(shù)據(jù)，在流動結(jié)構(gòu)簡單時，RNG k-ε湍流模型能有效減少計算所需資源，但本文涉及跨尺度流動問題，使用SST模型效果更佳。

2 結(jié)果分析

2.1 動、靜環(huán)織構(gòu)作用效果規(guī)律

圖6為不同深徑比情況下動、靜環(huán)織構(gòu)對端面最高溫度的降低幅度曲線圖。圖6a與圖6b表明，與靜環(huán)織構(gòu)相比，動環(huán)織構(gòu)強化換熱效果遠(yuǎn)勝，降溫幅度Δt平均高出5倍左右。除此之外，對不同深徑比γ下的動、靜環(huán)織構(gòu)進行分析發(fā)現(xiàn)，Δt隨γ的變化曲線存在轉(zhuǎn)折點。動環(huán)曲線的轉(zhuǎn)折點在γ=0.13附近，γ≥0.13時，改變γ會對Δt產(chǎn)生顯著影響，γ減小0.09，Δt增大了100%（5 ℃）；γ<0.13時，改變γ對Δt的影響較小，γ減小0.05，Δt僅增大了10%（1 ℃）。靜環(huán)也存在γ-Δt曲線突變點（對應(yīng)γ=0.1），但與動環(huán)不同，當(dāng)γ>0.1時，織構(gòu)換熱效果變化不顯著，γ減小0.3，Δt可以認(rèn)為不變（增大了0.7 ℃，33%）；當(dāng)γ<0.1時，γ減小0.02，Δt就已經(jīng)達到了之前的總和（0.7 ℃），但總體而言，靜環(huán)織構(gòu)的換熱效果差，降溫幅度小，不推薦在此工況下使用。由于動環(huán)換熱效果遠(yuǎn)勝于靜環(huán)換熱效果，進一步擴大轉(zhuǎn)速范圍，討論轉(zhuǎn)速n為1800～12 600 r/min工況下不同深徑比動環(huán)織構(gòu)的降溫幅度變化規(guī)律，如圖6c所示。五種織構(gòu)的Δt/tnmax均隨轉(zhuǎn)速提高而減小，但大深徑比織構(gòu)在經(jīng)歷初期的驟降之后，在中高轉(zhuǎn)速下Δt/tnmax的減小幾乎可以忽略不計，而小深徑比織構(gòu)Δt/tnmax近似呈直線下降。

2.2 動、靜環(huán)織構(gòu)強化換熱機理差異性分析

圖7為密封腔內(nèi)部軸對稱平面的漩渦強度云圖，圖7a所示為主軸轉(zhuǎn)速n=1800 r/min時的內(nèi)部漩渦強度，圖7b所示為主軸轉(zhuǎn)速n=12 600 r/min時的內(nèi)部漩渦強度，主渦集中于動環(huán)和靜環(huán)外部，近似于泰勒庫艾特流［35］。兩個主渦的產(chǎn)生機理各不相同，靠近軸末壓蓋（圖7中Z軸正方向末端）的漩渦（渦1）是由壓力出口的回流撞擊軸末壓蓋和靜環(huán)外周形成的。動環(huán)外側(cè)漩渦（渦2）的形成原因則較為復(fù)雜，渦1壓迫了入口流體的流動范圍，迫使入口流體僅能從渦1的右側(cè)穿過，如此形成Y軸正向速度，但出口流道由于回流占據(jù)了部分范圍而變得較為狹窄，因此部分流體將會沿Z軸負(fù)向回流并撞擊頂部壁面，又因為入口流的進入，渦2的Z向范圍也受到了限制。在不同轉(zhuǎn)速下均觀察到了這一現(xiàn)象，但當(dāng)轉(zhuǎn)速提高時，渦心的漩渦強度增大，從0.8 （n=1800 r/min）增大至0.95 （n=12 600 r/min）；對設(shè)置織構(gòu)的機械密封進行同樣的分析可知，兩個主渦的強度和位置幾乎不受織構(gòu)的影響。

進行換熱機理差異性分析時需要對織構(gòu)內(nèi)和織構(gòu)周邊的換熱效率分布和速度場等進行分析，由于本文織構(gòu)深度較淺（0.2 mm），側(cè)壁面面積過小，換熱能力可以忽略不計，因此本文將主要針對織構(gòu)底面和織構(gòu)周邊的外周壁面的換熱效果進行研究。選取γ=0.13作為分析對象，圖8為n=1800 r/min、γ=0.13時靜環(huán)和動環(huán)外周面換熱效率分布云圖，其中ηj0、ηd0 分別為n=1800 r/min時的靜環(huán)外周面換熱效率和動環(huán)外周面換熱效率?？棙?gòu)上游和織構(gòu)下游以圓心為界，主流進入側(cè)為織構(gòu)上游，主流離開側(cè)為織構(gòu)下游，織構(gòu)上游和織構(gòu)下游包括織構(gòu)外部附近的外圓周壁面，特別定義織構(gòu)下游較長的顯著影響區(qū)為尾跡區(qū)。

在上文分析的基礎(chǔ)上，選取了圖9所示的切片位置作進一步分析。圖9a所示為僅在動環(huán)設(shè)置織構(gòu)時的切片位置，圖9b所示為僅在靜環(huán)設(shè)置織構(gòu)時的切片位置。無論密封環(huán)外周面是否設(shè)置織構(gòu)，在圖1和圖2所示的幾何及邊界條件下均具有周向軸對稱性，因此，選取中軸面即X=0截面進行切片；Z=-1.5 mm和Z=4 mm位置分別為動環(huán)織構(gòu)和靜環(huán)織構(gòu)最佳換熱效果的中心平面；圖9中①和②的虛線圓所示范圍是綜合了周向軸對稱性和外周面換熱效果選取的位置。

圖10為Z=-1.5 mm橫截面流線圖，其中vm為流體流速，通過三維速度矢量求模獲得。與半球形織構(gòu)類似［39-41］，圓柱形織構(gòu)也可分為流動分離區(qū)、流動循環(huán)區(qū)、沖擊區(qū)和流動重附著區(qū)。入口側(cè)的流動分離區(qū)與半球形織構(gòu)類似，出現(xiàn)在幾何突變點附近，此區(qū)域換熱效率低，換熱效率值接近于1；流動循環(huán)區(qū)在動環(huán)中表現(xiàn)為規(guī)則漩渦，在靜環(huán)中表現(xiàn)為雜亂流動，但總體均沿幾何結(jié)構(gòu)向下游流動，且此部分流體與織構(gòu)底面不發(fā)生接觸，對織構(gòu)底面換熱效率無增強作用；沖擊區(qū)成因與半球形織構(gòu)不同，是織構(gòu)入口位置循環(huán)區(qū)之外的流體在進入織構(gòu)后形成垂直流動并沖擊織構(gòu)底面形成的，這一流動破壞了原本近壁面區(qū)域的熱邊界層，使織構(gòu)入口附近的換熱效率值接近于10；動環(huán)織構(gòu)內(nèi)大部區(qū)域的流體流動均以近似層流形式進行，相應(yīng)熱邊界層厚度也開始增大，使動環(huán)織構(gòu)底面換熱效率值規(guī)則地減??；最后，流體在出口附近出現(xiàn)了二次流動分離和離開織構(gòu)后的流動再附著，二次流動分離產(chǎn)生了圖8中尾跡區(qū)中間的低換熱區(qū)域，而流動再附著形成了高換熱長尾跡區(qū)。

圖11為Z=4 mm橫截面流線圖，靜環(huán)織構(gòu)內(nèi)部流體流動速度極低，同時內(nèi)部流動狀態(tài)也與動環(huán)織構(gòu)內(nèi)有較大的不同。入口分離區(qū)仍然位于幾何突變區(qū)附近，但循環(huán)區(qū)范圍增大，且深度范圍擴大至接近h0，同時在下游位置出現(xiàn)了二次的循環(huán)區(qū)，2個循環(huán)區(qū)不僅壓縮了沖擊區(qū)范圍，還增大了局部熱邊界層厚度，因此靜環(huán)內(nèi)部換熱效率值極小。最后一個區(qū)域是下游附近的流動再附著區(qū)，這一區(qū)域從下游循環(huán)區(qū)的上部開始，一直延續(xù)至下游較遠(yuǎn)的位置，對應(yīng)于圖8中靜環(huán)織構(gòu)下游ηj0>1.5的區(qū)域。

圖12為織構(gòu)內(nèi)部三維流線圖，為了能更直觀地顯示流線全貌，對固體壁面進行了透明處理。圖12b表明靜環(huán)織構(gòu)內(nèi)部流體進入織構(gòu)后會分為兩股，分別以紅、藍(lán)流線表示，紅色流線沿逆時針方向流動后緊貼壁面流出，藍(lán)色流線在織構(gòu)內(nèi)部形成較大范圍的再循環(huán)流動，隨后順時針流出，藍(lán)色流線部分的流體在流出織構(gòu)時出現(xiàn)了與動環(huán)相同的噴射流，這是圖8中靜環(huán)織構(gòu)下游僅有一塊再附著區(qū)且高換熱區(qū)范圍小的主要原因。圖12a中動環(huán)織構(gòu)入口流分為了兩層，下層貼壁流體以紅色流線表示，這部分流體產(chǎn)生了圖10中的入口流動分離區(qū)、循環(huán)區(qū)和部分沖擊區(qū)，隨后流體沿著側(cè)壁以漩渦的形式向兩側(cè)流動，離開織構(gòu)后在外部形成兩塊再附著區(qū)并進入下一織構(gòu)，是織構(gòu)下游長條狀尾跡區(qū)的成因；上層流體以藍(lán)色流線表示，這部分流體進入后，形成了圖10中部分沖擊區(qū)和下游附近的近似層流區(qū)，最后在下游織構(gòu)側(cè)壁產(chǎn)生分離區(qū)和噴射流，分離區(qū)和噴射流導(dǎo)致尾跡區(qū)中間出現(xiàn)低換熱區(qū)域。

2.3 深徑比對動、靜環(huán)織構(gòu)強化換熱的影響分析

2.1節(jié)已經(jīng)說明，深徑比對織構(gòu)強化換熱機理存在影響，本節(jié)將截取圖9中①和②虛線圓周邊區(qū)域，進一步分別分析動、靜環(huán)織構(gòu)γ發(fā)生變化時換熱效果的變化規(guī)律。圖13為各深徑比動環(huán)外圓周換熱效率分布云圖，可以看出，當(dāng)γ減小時，織構(gòu)底部高換熱效果的相對面積開始減小，表明織構(gòu)的整體換熱效率有隨γ減小而降低的趨勢；當(dāng)γ<0.13時，織構(gòu)的尾跡區(qū)開始受到下一排織構(gòu)的幾何尺度限制，尾跡區(qū)范圍減小，尾跡區(qū)換熱效果變差。這解釋了為何織構(gòu)γ<0.13時，端面最高溫度的降溫幅度隨γ減小而增大的趨勢開始急劇減緩。

圖14為各深徑比靜環(huán)外圓周換熱效率分布云圖。與動環(huán)相反，當(dāng)γ減小時，靜環(huán)織構(gòu)內(nèi)部換熱效果增強，且在本文研究的深徑比范圍內(nèi)，靜環(huán)織構(gòu)的主要換熱位置都是流體離開織構(gòu)后形成的尾跡區(qū)，改變深徑比僅僅改變了尾跡區(qū)的范圍和角度。

圖15為γ=0.08時動環(huán)織構(gòu)內(nèi)部主流的橫截面流線圖，可以看出，深徑比減小時，主要流場變化不大，入口的循環(huán)區(qū)和沖擊區(qū)也仍然存在，這主要是因為深度方向的幾何尺寸沒有發(fā)生變化，導(dǎo)致在入口附近小空間內(nèi)，深度方向的幾何突變幾乎沒有改變流場，但直徑尺寸的大小影響了入口循環(huán)流動形成的漩渦尺寸。當(dāng)深徑比過小時，入口處循環(huán)流動形成的漩渦將在深度方向增大，使得垂直流減弱，沖擊區(qū)范圍減小，高換熱區(qū)域縮??；其次，減小的深徑比使得織構(gòu)底面總面積增大，導(dǎo)致高換熱區(qū)域的相對面積減小，減弱了織構(gòu)的換熱能力。

圖16為γ=0.08時靜環(huán)織構(gòu)內(nèi)部主流的橫截面流線圖。對于靜環(huán)織構(gòu)，γ變小后入口位置循環(huán)區(qū)幾乎沒有發(fā)生變化，但右下角的二次循環(huán)區(qū)占比顯著減小了，這變相增大了沖擊區(qū)的面積占比。循環(huán)區(qū)是低換熱區(qū)，沖擊區(qū)是高換熱區(qū)，二者的相對面積變化使得靜環(huán)換熱效率隨深徑比減小而提高。此外，小γ時靜環(huán)織構(gòu)內(nèi)部流速仍然接近于0，因此雖然靜環(huán)內(nèi)部換熱效率值有所增大，仍然比動環(huán)內(nèi)部換熱效率值小得多。

雖然小深徑比動環(huán)織構(gòu)和大深徑比動環(huán)織構(gòu)在主流方向橫截面的流線相差不大，但上述兩種動環(huán)織構(gòu)的內(nèi)部垂直主流方向流場則出現(xiàn)了較大的不同。圖17為不同深徑比動環(huán)織構(gòu)X=0截面局部流線圖，可以看出，兩種情況下織構(gòu)內(nèi)部流體均存在左右兩側(cè)同時進入織構(gòu)的情況， γ=0.2時左側(cè)入口流體占據(jù)主體，右側(cè)入口流體僅在進入織構(gòu)后的小范圍內(nèi)生成漩渦；γ=0.08時，左右兩側(cè)進入織構(gòu)的流體在織構(gòu)中心區(qū)產(chǎn)生了碰撞，并且在碰撞界面附近形成了兩個強烈的漩渦，這些漩渦外層的流體又在織構(gòu)底面附近呈現(xiàn)出近似層流的流態(tài)，降低了小深徑比織構(gòu)中心及下游位置的換熱效率，但對于靜環(huán)織構(gòu)，并無這一現(xiàn)象。

2.4 主軸轉(zhuǎn)速對動環(huán)織構(gòu)強化換熱的影響分析

選取γ=0.1和γ=0.2分別代表大深徑比和小深徑比兩種情況進行進一步的分析。圖18為γ=0.1和γ=0.2時n分別為1800 r/min和12 600 r/min工況下的動環(huán)織構(gòu)換熱效率云圖，可以看出，在轉(zhuǎn)速提高后，換熱效率值減小，驗證了前文所獲得的規(guī)律，即在高速情況下，織構(gòu)的換熱效果變差。此外，轉(zhuǎn)速提高后，尾跡區(qū)換熱效率值也隨之減小。

3 結(jié)論

（1）織構(gòu)強化換熱機理是外周面織構(gòu)的設(shè)置促成了垂直于壁面的流動，該流動通過破壞流體繞流的熱邊界層，增大了局部努塞爾數(shù)，增強了換熱效果。

（2）靜環(huán)外周面織構(gòu)內(nèi)部及其周邊流速均低于動環(huán)外周面織構(gòu)對應(yīng)區(qū)域的流速，熱邊界層較厚，因此前者換熱效果遠(yuǎn)不如后者。

（3）當(dāng)深徑比減小時，入口處循環(huán)渦旋將會增大，入口垂直流受到循環(huán)渦旋影響，占據(jù)范圍減小，使得沖擊區(qū)范圍被壓縮，局部努塞爾數(shù)減小。

（4）轉(zhuǎn)速增大時，流體在織構(gòu)內(nèi)部平行于壁面方向流速的增幅遠(yuǎn)大于垂直于壁面方向流速的增幅，流體流動方向與壁面間夾角減小，導(dǎo)致織構(gòu)整體換熱效果變差。此外，對于小深徑比織構(gòu)，將會在下游形成新的循環(huán)區(qū)，使局部努塞爾數(shù)急劇減小。

參考文獻：

［1］ LUAN Z G， KHONSARI M M. Analysis of Conjugate Heat Transfer and Turbulent Flow in Mechanical Seals［J］. Tribology International， 2009， 42（5）：762-769.

［2］ TAKAMI M R， GERDROODBARY M B， GANJI D D. Thermal Analysis of Mechanical Face Seal Using Analytical Approach［J］. Thermal Science and Engineering Progress， 2018， 5：60-68.

［3］ MOSAVAT M， MORADI M， TAKAMI M R， et al. Heat Transfer Study of Mechanical Face Seal and Fin by Analytical Method［J］. Engineering Science and Technology—an International Journal， 2018， 21（3）：380-388.

［4］ LUAN Z G， KHONSARI M M. Heat Transfer Correlations for Laminar Flows within a Mechanical Seal Chamber［J］. Tribology International， 2009， 42（5）：770-778.

［5］ 顧永泉， 吳宗祥， 王麗娟. 機械密封的端面溫度［J］. 流體工程， 1985， 4：3-10.

GU Yongquan， WU Zongxiang， WANG Lijuan. End Face Temperature of Mechanical Seal［J］. Fluid Machinery， 1985， 4：3-10.

［6］ 張淑敏， 胡麗國， 孟祥鎧. 超高速燃?xì)鉁u輪泵機械密封的分析與研究［J］. 流體機械， 2010， 40（10）：23-27.

ZHANG Shumin， HU Liguo， MENG Xiangkai. Analysis of a Mechanical Seal for Superspeed Gas Turbopump［J］. Fluid Machinery， 2010， 40（10）：23-27.

［7］ 彭旭東， 金杰， 李定， 等. 高速渦輪泵機械密封端面溫度變化規(guī)律研究［J］. 摩擦學(xué)學(xué)報， 2019， 39（3）：313-318.

PENG Xudong， JIN Jie， LI Ding， et al. Analysis of Face Temperature in Mechanical Seal Applied in the High Speed Turbopump［J］. Tribology， 2019， 39（3）：313-318.

［8］ KENNEDY F E， KARPE S A. Thermocracking of a Mechanical Face Seal［J］. Wear， 1982， 79（1）：21-36.

［9］ NIAN X， KHONSARI M M. Thermal Performance of Mechanical Seals with Textured Side-wall［J］. Tribology International， 2012， 45（1）：1-7.

［10］ NIAN X， KHONSARI M M. Improving Thermal Performance of Mechanical Seals with Surface Texturing［J］. Proceedings of the Institution of Mechanical Engineers，Part J：Journal of Engineering Tribology， 2015， 229：350-361.

［11］ 周宇坤， 彭旭東， 趙文靜， 等. 機械密封動環(huán)外周表面織構(gòu)換熱機理及結(jié)構(gòu)優(yōu)化［J］. 摩擦學(xué)學(xué)報， 2020， 40（4）：538-550.

ZHOU Yukun， PENG Xudong， ZHAO Wenjing， et al. Study on Heat Transfer Mechanism and Optimization of Circumferential Texture of Mechanical Seal［J］. Tribology， 2020， 40（4）：538-550.

［12］ AFANASYEV V N， CHUDNOVSKY Y P， LEONTIEV A I， et al. Turbulent Flow Friction and Heat Transfer Characteristics for Spherical Cavities on a Flat Plate［J］. Experimental Thermal and Fluid Science， 1993， 7（1）：1-8.

［13］ LIN Y L， SHIH T I P， CHYU M K. Computations of Flow and Heat Transfer in a Channel with Rows of Hemispherical Cavities［C］∥Proceedings of the ASME 1999 International Gas Turbine and Aeroengine Congress and Exhibition. Indianapolis：ASME， 1999：99-GT-263.

［14］ CHYU M K， YU Y， DING H， et al. Concavity Enhanced Heat Transfer in an Internal Cooling Passage［C］∥Proceedings of the ASME 1997 International Gas Turbine and Aeroengine Congress and Exhibition. Orlando：ASME， 1997：97-GT-437.

［15］ MAHMOOD G I， HILL M L， NELSON D L， et al. Local Heat Transfer and Flow Structure on and above a Dimpled Surface in a Channel［J］. Journal of Turbomachinery， 2001， 123（1）：115-123.

［16］ MAHMOOD G I， LIGRANI P M. Heat Transfer in a Dimpled Channel：Combined Influences of Aspect Ratio， Temperature Ratio， Reynolds Number， and Flow Structure［J］. International Journal of Heat and Mass Transfer， 2002， 45（10）：2011-2020.

［17］ MOON S W， LAU S C. Turbulent Heat Transfer Measurements on a Wall with Concave and Cylindrical Dimples in a Square Channel［C］∥Proceedings of the ASME Turbo Expo 2002：Power for Land， Sea， and Air. Amsterdam：ASME， 2002：459-467.

［18］ DYACHENKO A Y， TEREKHOV V I， YARYGINA N I. Turbulent Flow Past a Transverse Cavity with Inclined Side Walls［J］. Journal of Applied Mechanics and Technical Physics， 2007， 48：486-491.

［19］ BUTT U， EGBERS C. Flow Structure due to Hexagonal Cavities and Bumps on a Plate Surface［J］. Thermophysics and Aeromechanics， 2016， 23：839-847.

［20］ TURNOW J， KORNEV N， ISAEV S， et al. Vortex Mechanism of Heat Transfer Enhancement in a Channel with Spherical and Oval Dimples［J］. Heat Mass Transfer， 2011， 47：301-313.

［21］ LIGRANI P M， MAHMOOD G I， HARRISON J L， et al. Flow Structure and Local Nusselt Number Variations in a Channel with Dimples and Protrusions on Opposite Walls［J］. International Journal of Heat and Mass Transfer， 2001， 44（23）：4413-4425.

［22］ LIU Jian， SONG Yidan. Numerical Modeling Flow and Heat Transfer in Dimpled Cooling Channels with Secondary Hemispherical Protrusions［J］. Energy， 2015， 79：1-19.

［23］ XIE Yonghui， SHI Dongbo， SHEN Zhongyang. Experimental and Numerical Investigation of Heat Transfer and Friction Performance for Turbine Blade Tip Cap with Combined Pin-Fin-Dimple/Protrusion Structure［J］. International Journal of Heat and Mass Transfer， 2017， 104：1120-1134.

［24］ LIU Chaoqun， WANG Yiqian， YANG Yong， et al. New Omega Vortex Identification Method［J］. Physics， Mechanics & Astronomy， 2016， 59（8）：684711.

［25］ YU Y F， SHRESTHA P， ALVAREZ O， et al. Investigation of Correlation between Vorticity， Q， λci， λ2， Δ and Liutex［J］. Computers & Fluids， 2021， 225：104977.

［26］ 陳文杰， 孟祥鎧， 王玉明， 等. 機械密封織構(gòu)化表面粗糙度效應(yīng)的有限元模型與摩擦學(xué)特性分析［J］. 摩擦學(xué)學(xué)報， 2019， 39（5）：523-530.

CHENG Wenjie， MENG Xiangkai， WANG Yu-ming， et al. Finite Element Model and Tribological Properties Analysis of Textured Mechanical Seals Considering Roughness Effect［J］. Tribology， 2019， 39（5）：523-530.

［27］ 程健梁. 機械密封復(fù)合織構(gòu)化端面密封性能與形貌演化［D］. 杭州：浙江工業(yè)大學(xué)， 2020.

CHENG Jianliang. Sealing Performance and Morphology Evolution of Composite Textured End Face of Mechanical Seal［D］. Hangzhou：Zhejiang University of Technology， 2020.

［28］ ADJEMOUT M， BRUNETIRE N， BOUYER J. Numerical Analysis of the Texture Effect on the Hydrodynamic Performance of a Mechanical Seal［J］. Surface Topography：Metrology and Properties， 2016， 4（1）：014002.

［29］ TEKRIWAL P. Heat Transfer Predictions with Extended k-ε Turbulence Model in Radial Cooling Ducts Rotating in Orthogonal Mode［J］. Journal of Heat Transfer， 1994， 116：369-380.

［30］ WANG L B，WANG Q W， HE Y L， et al. Experimental and Numerical Study of Developing Turbulent Flow and Heat Transfer in Convergent/Divergent Square Ducts［J］. Heat and Mass Transfer， 38：399-408.

［31］ 陶文銓. 數(shù)值傳熱學(xué)［M］. 2版.西安：西安交通大學(xué)出版社， 2001：353-362.

TAO Wenquan. Numerical Heat Transfer［M］. 2nd ed. Xian：Xian Jiaotong University Press， 2001：353-362.

［32］ MENTER F R. Review of the Shear-stress Transport Turbulence Model Experience from an Industrial Perspective［J］. International Journal of Computational Fluid Dynamic， 2009， 23（4）：305-316.

［33］ 王震. 基于鈍體風(fēng)壓數(shù)據(jù)驅(qū)動的SST k-ω湍流模型參數(shù)修正［D］. 哈爾濱：哈爾濱工業(yè)大學(xué)， 2020.

WANG Zhen. Parameter Modification of SST k-ω Turbulence Model Based on Bluff Body Wind Pressure Data［D］. Harbin：Harbin Institute of Technology， 2020.

［34］ STEPHEN K R. Coherent Motions in the Turbulent Boundary Layer［J］. Annual Review of Fluid Mechanics， 1991， 23：601-639.

［35］ 丁仕成. 溝槽結(jié)構(gòu)參數(shù)對湍流Taylor-Couette流動及其換熱特性影響的研究［D］. 鎮(zhèn)江：江蘇大學(xué)， 2018.

DING Shicheng. Research of Slit Structure Parameters Effect on the Turbulent Taylor-Couette Vortex Flow and Heat Transfer Characteristics［D］. Zhenjiang：Jiangsu University， 2018.

［36］ 顧永泉. 機械密封實用技術(shù)［M］. 北京：機械工業(yè)出版社， 2001：89-90.

GU Yongquan. Practical Technology of Mechanical Seal［M］. Beijing：Machinery Industry Press， 2001：89-90.

［37］ 彭旭東， 謝友柏， 顧永泉. 機械密封端面溫度的確定［J］. 化工機械， 1996， 23（6）：333-366.

PENG Xudong， XIE Youbai， GU Yongquan. Determination of the End Face Temperature of Mechanical Seal［J］. Chemical Engineering & Machinery， 1996， 23（6）：333-366.

［38］ GRIFFITH T S， AL-HADHRAMI L， HAN J. Heat Transfer in Rotating Rectangular Cooling Channels （AR=4） with Dimples［J］. Journal of Turbomachinery， 2003， 125（3）：555-563.

［39］ NOURIN F N， BLUM B L， AMANO R S. Evaluation of Heat Transfer Enhancement on Rotational Gas Turbine Blade Internal Cooling Channel with Dimpled Surface［J］. Journal of Energy Resources Technology， 2022， 144（11）：112105.

［40］ NOURIN F N， AMANO R S. Experimental Study on Flow Behavior and Heat Transfer Enhancement with Distinct Dimpled Gas Turbine Blade Internal Cooling Channel［J］. Journal of Energy Resources Technology， 2022， 144（7）：072101.

［41］ PARVIZ M， NORI A O， ROBERT L P， et al. Experimental and Computational Investigation of Flow and Thermal Behavior of a Mechanical Seal［J］. Tribology Transactions， 1999， 42（4）：731-738.

［42］ 張冠敏. 復(fù)合波紋板式換熱器強化傳熱機理及傳熱特性研究［D］. 濟南：山東大學(xué)， 2006.

ZHANG Guanmin. Research on Heat Transfer Enhancement Mechanism and Heat Transfer Characteristic of Compound Corrugation Plate Heat Exchanger［D］. Jinan：Shandong University， 2006.