袁 聰，朱麗莎，杜尊令，李 鶴,，宋錦春

(1.肇慶學(xué)院 機械與汽車工程學(xué)院，廣東 肇慶 526061； 2.東北大學(xué) 機械工程與自動化學(xué)院，遼寧 沈陽 110819)

引言

在液壓傳動技術(shù)中，溢流閥、順序閥等壓力控制元件普遍采用錐閥結(jié)構(gòu)。在系統(tǒng)運行期間，錐閥需要長時間承受較高的壓差，容易在開口區(qū)域誘發(fā)空化，形成空化射流。由于水的黏度遠小于礦物質(zhì)油，水壓錐閥內(nèi)的空化現(xiàn)象普遍更加嚴重，并對產(chǎn)品性能存在重要的影響。而當前，水壓技術(shù)在采礦、清洗、建筑等領(lǐng)域的應(yīng)用愈加廣泛[1-3]。通過分析水壓錐閥內(nèi)部空化射流的動態(tài)特性，探索空化分布規(guī)律及作用機制，能夠為產(chǎn)品性能優(yōu)化提供充足基礎(chǔ)及科學(xué)依據(jù)。

早期關(guān)于水壓錐閥的研究主要以實驗手段為主。OSHIMA S等[4]通過半切模型機對不同結(jié)構(gòu)的水壓錐閥開展了實驗研究，揭示了多種壓差條件下的空化分布情況。LIU Y等[5]測量了兩級式水壓錐閥的流動及空化特性，分析了不同工況條件下的水壓錐閥流量-壓差變化趨勢。LIU Y等[6]又對比了不同閥座結(jié)構(gòu)水壓錐閥的空化特性，證實閥座結(jié)構(gòu)對水壓錐閥的空化抑制能力有重要的影響。隨著計算機技術(shù)的進步及流場模擬算法的發(fā)展，近年來的研究大量采用CFD方法對錐閥空化射流進行流場模擬。NIE S等[7]分析了兩級式水壓錐閥的壓力分布，并針對空化抑制能力提出了結(jié)構(gòu)參數(shù)的優(yōu)化方案。LIANG J等[8]通過數(shù)值模擬方法分析了入口壓力脈動對水壓錐閥流量及空化特性的影響。HAN M等[9]通過流場仿真研究了水壓錐閥在空化狀態(tài)下的液流力變化趨勢。張增猛等[10]基于CFD數(shù)值研究對水壓球形錐閥進行了穩(wěn)態(tài)液流力分析。謝海波等[11]對不同形態(tài)錐閥的液流力開展了數(shù)值研究。

上述的研究工作從不同的角度揭示了水壓錐閥的空化特性。較低的流體黏度使水壓錐閥內(nèi)部射流的層流-湍流轉(zhuǎn)捩異常迅猛，與空化瞬態(tài)演變密切相關(guān)，但由于極短的時間和長度尺度，難以開展針對性的分析，因此該方面的研究成果仍較為欠缺。本研究擬采用CFD方法對水壓錐閥的空化射流動態(tài)行為開展三維動態(tài)仿真，根據(jù)計算結(jié)果開展流場特征結(jié)構(gòu)分析，探索空化分布情況及其潛在機理，研究不同空化結(jié)構(gòu)的相互關(guān)聯(lián)，揭示空化射流的周期性演變過程。

1 數(shù)值方法

有研究表明[12]，對于高速射流，需要在數(shù)值方法中加入液體的可壓縮性才能準確預(yù)測微小尺寸的漩渦結(jié)構(gòu)。本研究采用VOF算法處理空化兩相流，流相輸運方程為：

(1)

式中，αl，αv—— 液相及氣相的體積分數(shù)

ρl，ρv—— 液相及氣相的密度

U—— 氣液混相的平均速度

Ur—— 兩相間的相對滑移速率

(2)

其中，cα為氣液界面壓縮系數(shù)，為保證數(shù)值穩(wěn)定普遍取值低于1，本計算中取0.7。

速度方程為：

(3)

其中，σκ▽α為表面張力項，κ為相間的交界面曲率；p為流場壓力；σ為流體的張力系數(shù)。

采用Schnerr空化模型計算相間質(zhì)量轉(zhuǎn)換速率：

(4)

式中，pv—— 水在常溫下的汽化壓力

r，Rb—— 空化模型常系數(shù)

該數(shù)值方法的推導(dǎo)及求解可參考文獻[12]。在OpenFOAM開源平臺上通過PIMPLE循環(huán)對上述模型進行數(shù)值求解，時間相關(guān)項使用Backward離散格式，流相輸運方程使用MULES算法，將速度方程進行半離散處理構(gòu)建泊松壓力方程，并使用預(yù)測步疊加修正步進行速度壓力解耦。為了捕捉精細的漩渦結(jié)構(gòu)，使用圖1所示的網(wǎng)格模型進行三維計算，兩側(cè)面使用對稱邊界條件，入口及出口選用非反射性壓力條件，其余設(shè)置為無滑移壁面。

圖1 倒角式水壓錐閥結(jié)構(gòu)示意圖及計算網(wǎng)格

2 計算結(jié)果及分析

1) 空化結(jié)構(gòu)分布及作用機制

圖2展示了不同工況下的空化結(jié)構(gòu)分布情況。本研究涉及到的數(shù)值計算中，入口邊界的全壓力設(shè)置為5.0 MPa，通過調(diào)整出口壓力獲得不同空化強度的工況。當壓差Δp=2.5 MPa時，空化結(jié)構(gòu)首先出現(xiàn)在倒角閥座流道內(nèi)部，當壓差Δp=3.0 MPa時，倒角閥座流道出口下游處亦開始出現(xiàn)小尺寸空化結(jié)構(gòu)，與實驗拍攝照片較為吻合。另外照片中可以觀測到閥芯后沿有輕微的空化，此處的空化結(jié)構(gòu)也被數(shù)值方法預(yù)測到(可參見圖3)。當壓差Δp為4.0，4.4 MPa時，空化的強度明顯增強，并呈現(xiàn)如下3個分布規(guī)律：倒角閥座流道內(nèi)部的空化延伸至出口下游處；在流道出口的下游區(qū)域形成尺寸較大的汽泡群結(jié)構(gòu)；該汽泡群結(jié)構(gòu)可越過閥芯后沿，同時在下游處的汽泡結(jié)構(gòu)明顯增大。在壓差Δp=4.4 MPa的工況條件下，數(shù)值結(jié)果顯示的后沿下游空化結(jié)構(gòu)分布與后沿上游空化區(qū)域相比較為稀疏，并預(yù)測有大尺寸汽泡結(jié)構(gòu)產(chǎn)生，其覆蓋區(qū)域為4.6 mm×3.2 mm；實驗照片中，在該區(qū)域呈現(xiàn)出的相似的空化分布特征，后沿下游的空化范圍約為4.4 mm×3.4 mm。實驗照片中，流道出口下游至閥芯后沿中間區(qū)域的空化分布極為密集，表現(xiàn)為類似連續(xù)性的空化結(jié)構(gòu)；而數(shù)值計算結(jié)果盡管成功預(yù)測該區(qū)域密集的空化分布，但該處的空化結(jié)構(gòu)具有明顯的非連續(xù)性特征。這種差異極大可能是與空化結(jié)構(gòu)具有較高的運動速度及劇烈的演變有關(guān)，而反射型的高速拍攝需要有相對長的曝光時間，另外不同方位角的空化結(jié)構(gòu)亦能反射光線到照片中從而造成光亮的背景，難以在照片中識別該區(qū)域的空化結(jié)構(gòu)細微形態(tài)。

圖2 不同壓差條件下的空化形態(tài)(空化結(jié)構(gòu)通過αl=70%等值面識別)

圖3為空化射流的瞬時流場特征結(jié)構(gòu)，從左到右分別為漩渦結(jié)構(gòu)分布、空化結(jié)構(gòu)分布、空化-漩渦層疊圖及速度云圖，各運行工況與圖2一致。圖3a～圖3d為不同工況的數(shù)值模擬結(jié)果，層疊圖中，深色空化結(jié)構(gòu)通過αl=70%等值面識別，淺色漩渦結(jié)構(gòu)通過Q=1×1010等值面識別。在倒角閥座流道入口處，從速度云圖可以觀測到分離流現(xiàn)象，緊隨入口的閥座表面區(qū)域存在低壓區(qū)，從而誘發(fā)附著型空化。在倒角閥座出口外的區(qū)域，產(chǎn)生大尺寸的漩渦結(jié)構(gòu)，漩渦中心的低壓區(qū)域誘發(fā)空化，形成漩渦空化。各工況的空化-漩渦層疊圖表明，此處的空化結(jié)構(gòu)與漩渦中心吻合，同時大尺寸的漩渦結(jié)構(gòu)主要分布在射流勢核的自由剪切層一側(cè)，該區(qū)域的漩渦空化結(jié)構(gòu)亦相應(yīng)地集中在自由剪切層，此種分布規(guī)律與圖2a的照片呈現(xiàn)的空化結(jié)構(gòu)較為一致。從三維視角下的漩渦分布可以清晰判別出，該區(qū)域的漩渦空化并非是連續(xù)的空腔結(jié)構(gòu)。同時，閥芯表面有薄層型的漩渦產(chǎn)生，在壓差Δp=4.4 MPa的工況條件下射流勢核的壁面剪切層處形成具有薄層形態(tài)的漩渦空化，并且該空化結(jié)構(gòu)尺寸較小，射流勢核在錐形閥芯的發(fā)散趨勢使壁面上的漩渦結(jié)構(gòu)逐漸削弱消散，薄層形態(tài)的空化結(jié)構(gòu)在下游區(qū)域逐漸消失。在閥芯后沿處，射流勢核亦有明顯的分離流現(xiàn)象，并在此處同樣形成附著型空化。

圖3 三維視角下的空化結(jié)構(gòu)、漩渦結(jié)構(gòu)及截面速度云圖

2) 空化結(jié)構(gòu)的動態(tài)演變過程

圖4為總氣相體積Vvapor隨時間t的變化趨勢，主要由數(shù)個大幅度波動組成，但其重復(fù)性趨勢并不穩(wěn)定，表現(xiàn)為夾雜有時間跨度較短的波動，原因在于整體的空化結(jié)構(gòu)演變涉及多種漩渦空化與附著型空化的耦合，與單純的附著型空化或噴射流漩渦空化的周期性空化行為相比，牽涉到多個作用機制。圖5展示了某個大幅度波動內(nèi)的空化結(jié)構(gòu)具體演變過程。圖5a～圖5j為不同時刻的空化結(jié)構(gòu)，各時刻對應(yīng)于圖5k中的黑色方形點，空化結(jié)構(gòu)通過等值面αl=90%識別，Δp=4.4 MPa。在t1～t3時刻，倒角閥座流道的附著型空化脫落并在出口下游處形成大尺寸氣泡群結(jié)構(gòu)；t4～t6時刻，閥座流道空化產(chǎn)生二次脫落，形成的漩渦空化在t7時刻與下游的氣泡群結(jié)構(gòu)合并，同時氣泡群結(jié)構(gòu)到達閥芯后沿，與該處的附著型空化結(jié)合，在這兩種機制的共同作用下氣相體積在t7時刻到達周期內(nèi)最高點；之后的t8～t10時刻，空化結(jié)構(gòu)逐漸潰滅，氣相體積不斷減少。從圖4的漩渦結(jié)構(gòu)分布可以得知，閥座流道后半部的漩渦已具備明顯的三維特性，脫落空化在出口下游形成的漩渦空化結(jié)構(gòu)差異明顯，而整體的空化演變又主要圍繞漩渦空化開展，因此全局空化體積的周期性變化趨勢并不穩(wěn)定。具體的空化動態(tài)過程模擬結(jié)果有待進一步研究。

圖4 氣相體積的時間變化 (Δp=4.4 MPa)

圖5 單個波動周期內(nèi)的空化結(jié)構(gòu)演變過程

3) 漩渦空化與附著型空化的耦合

閥座流道內(nèi)的附著型空化的脫落行為伴隨著三維漩渦結(jié)構(gòu)的產(chǎn)生，并且脫落的空化基本上位于漩渦中心區(qū)域，表明漩渦空化已經(jīng)在后半部流道內(nèi)形成。但是，由于該處的漩渦具有較高的三維特性，流道內(nèi)的漩渦空化已經(jīng)不具備周向?qū)ΨQ性，分布并無明顯規(guī)律性。出口下游處，漩渦結(jié)構(gòu)的尺寸迅速增長，內(nèi)部包裹的空化結(jié)構(gòu)亦隨之快速膨脹，周向相鄰的空化結(jié)構(gòu)可形成相互連接，發(fā)展成環(huán)形的空化鏈，如圖6所示，連續(xù)圖片的時間間隔為2×10-5s，深色的空化結(jié)構(gòu)通過等值面αl=70%識別，淺色的漩渦結(jié)構(gòu)通過等值面Q=1×1010識別，Δp=4.4 MPa。應(yīng)該明確指出，該空化鏈是由不同方位角的脫落空化通過周向連接而成，夾雜著較明顯的三維漩渦結(jié)構(gòu)，因此并不具備周向?qū)ΨQ性。

圖6 倒角閥座流道后部脫落空化形成的鏈狀空化結(jié)構(gòu)

氣泡群在經(jīng)過閥芯后沿處時，與該區(qū)域的附著型空化相遇。一方面，促使該處的附著型空化發(fā)生脫落，因此閥芯后沿處的附著空化的獨立演變行為被強行打斷。另一方面，附著型空化與不斷逼近的漩渦空化發(fā)生合并，原本正緩慢潰滅的漩渦空化將由于附著空化的介入再次迎來增長。經(jīng)過二次增長的漩渦空化在閥芯后沿的下游區(qū)域形成大尺寸的氣泡結(jié)構(gòu)，呈現(xiàn)稀疏分布特征并覆蓋一個較大的區(qū)域，如圖7所示，連續(xù)圖片的時間間隔為2×10-5s，空化結(jié)構(gòu)通過等值面αl=70%識別，Δp=4.4 MPa。漩渦空化與后沿處附著空化的耦合合理地揭示了實驗照片中該區(qū)域空化分布特征的形成過程。

圖7 漩渦空化與閥芯后沿附著空化的相互作用過程

3 結(jié)論

研究對倒角形水壓錐閥開展流場仿真，捕捉了關(guān)鍵的流場特征，并取得了以下結(jié)論：

(1) 倒角形閥座水壓錐閥內(nèi)部的空化主要集中在3個區(qū)域。倒角閥座流道入口處及閥芯后沿存在分離流現(xiàn)象，并在緊鄰的下游區(qū)域形成附著型空化；而倒角閥座流道出口至閥芯后沿中間區(qū)域，在射流勢核的自由剪切層產(chǎn)生大尺寸漩渦結(jié)構(gòu)，并誘發(fā)漩渦空化；

(2) 由于不同類型的空化之間存在耦合，因此空化結(jié)構(gòu)的全局演變的周期特性受到影響；

(3) 倒角閥座流道后部的脫落空化伴有三維漩渦，并在自由剪切層處尺寸增大并形成空化鏈結(jié)構(gòu)。而漩渦空化與閥芯后沿的附著空化經(jīng)過相互作用，發(fā)展成大尺寸的氣泡。