周秀英 白軍 姜梅 王超琪 顧巖城 吉華

(西藏農(nóng)牧學(xué)院水利土木工程學(xué)院，西藏 林芝 860000)

引言

西藏是我國的關(guān)鍵區(qū)域，對國家安全、環(huán)境保護和清潔能源供應(yīng)都有著重要作用[1]?！笆奈濉币?guī)劃計劃加快推動雅江下游、瀾滄江上游等多流域水電開發(fā)，預(yù)計到2025年末，西藏地區(qū)的水電建成和在建裝機容量將超過15GW[2]。由于西藏地區(qū)的平均海拔超過4km，平均氣壓低于65.25kPa，目前已建電站的運行情況調(diào)研結(jié)果顯示，此種特殊的地理環(huán)境導(dǎo)致了該地區(qū)水力機組的空化問題突出，對機組的功效和壽命造成了嚴重的影響[3]。羅紅英等[4]通過數(shù)值模擬，從理論上證明了高海拔地區(qū)空化規(guī)模更大、空化周期更長，水輪機的空蝕現(xiàn)象更加嚴重。因此，總結(jié)空化的發(fā)展歷程與研究現(xiàn)狀，深入研究空化現(xiàn)象，總結(jié)空化研究方法對于西藏地區(qū)的國家重點水利水電工程(如雅江下游)等具有重要意義。

空化會對水力機械等工程設(shè)備帶來極大的侵蝕，所以一直以來空化都是水動力學(xué)的研究熱點與難點，其本質(zhì)為相變，一般可簡述為“當(dāng)局部壓力降低，液體中固有的微小氣核發(fā)生膨脹以及潰滅的過程”，因此研究流體中氣泡空化行為的空泡動力學(xué)則是空化研究的基礎(chǔ)。

空泡動力學(xué)最初的建立可以追溯到1917年由Rayleigh提出的Rayleigh方程[4]，其研究了理想情況球形氣泡的運動過程。后續(xù)Plesset等[4]逐步考慮了表面張力、黏性等因素對Rayleigh方程進行了修正。雖然空泡的理論方程由于空化本身的強非線性等原因，計算結(jié)果與實際情況具有較大差異，無法直接運用于工程實際，但理論的發(fā)展催生了目前2種主流研究方法：實驗觀測和數(shù)值模擬。根據(jù)理論方程選擇合適參數(shù)建造的水洞等空化實驗平臺搭配Particle Image Velocimetry(PIV)觀測系統(tǒng)為空化認知提供了基礎(chǔ)。但實驗手段亦存在著成本過高、周期過長、可獲取數(shù)據(jù)有限等缺陷，局限性過高。數(shù)值模擬在理論方程的基礎(chǔ)上搭配高性能計算機，目前已成為了一種極為重要的空化現(xiàn)象研究手段。數(shù)值模擬是一種基于選定控制方程，對水力模型進行離散后，迭代求解的方法，相較于實驗手段，數(shù)值模擬有著成本低、時效快、易于調(diào)整工況參數(shù)等優(yōu)勢。綜上所述，本文對氣泡空化的數(shù)值模擬方法進行綜述，介紹了空化研究的發(fā)展歷程及空化基本原理；對空泡數(shù)值模擬方法現(xiàn)狀進行了總結(jié)與綜述，并探討了各方法的優(yōu)缺點；對空泡研究進行了總結(jié)與展望。

1 研究背景與意義

1.1 空泡動力學(xué)研究背景

1753年，科學(xué)家L Euler曾提出“當(dāng)水管中的某處壓強降至負值時，水會從管壁分離，形成間隙從而產(chǎn)生真空空間”的觀點。然而，這個所謂的真空空間實際上并非完全真空，而是由液體汽化形成的空穴。對空化現(xiàn)象的研究真正始于19世紀后期，由于船舶行業(yè)的興盛，隨之而來的是人們對于船舶螺旋槳等部位大面積剝蝕及運行效率下降的疑問。1917年Lord Rayleigh提出了著名的Rayleigh方程[4]，其形式如式(1)所示。

(1)

式中，R為氣泡半徑，是關(guān)于時間t的函數(shù)；pv為泡內(nèi)壓強；p∞為環(huán)境壓強。

該方程為空泡動力學(xué)奠定了基礎(chǔ)，但由于其忽略了液體的黏性、氣泡的表面張力、非冷凝氣體(NCG)、熱效應(yīng)等影響項，因此對于空化實際情況的分析并非完全合理。

1951年，Plesset在Rayleigh的基礎(chǔ)上，考慮了氣泡的表面張力及液體黏性，給出了經(jīng)典的Rayleigh-Plesset方程[4]，如式(2)所示。

(2)

式中，σ為表面張力，是關(guān)于時間t的函數(shù)；μ為黏性系數(shù)。

1952年Gilmore考慮了空泡潰滅后期產(chǎn)生的超聲速射流及流體的可壓縮性[4]，將Rayleigh方程進一步修正，如式(3)所示。

(3)

式中，c為聲速；H為焓值。

1980年Keller在Rayleigh-Plesset方程的基礎(chǔ)上考慮了超聲場中空泡的振蕩以及流體的壓縮性、黏性等影響因素，提出了Keller-Miksis方程[4]，如式(4)所示。

(4)

式中，M為馬赫數(shù)；pl為泡外壓強。

2023年，我國哈爾濱工程學(xué)院的張阿漫教授提出了氣泡統(tǒng)一方程[4]，形式如式(5)所示。

(5)

式中，H為泡內(nèi)外焓差。方程左側(cè)為氣泡遷移及環(huán)境耦合力、右側(cè)為體積加速度。該方程同時考慮邊界效應(yīng)、泡群、流場環(huán)境、遷移、可壓縮性、粘性、表面張力等因素，將不同源、不同尺度、不同環(huán)境下的氣泡動力學(xué)特性統(tǒng)一起來，并且統(tǒng)一了式(1)～(4)經(jīng)典方程，是空泡動力學(xué)中里程碑式的工作。

1.2 研究意義

目前空泡理論研究表明，空化是力學(xué)作用、化學(xué)腐蝕、電化學(xué)腐蝕、熱力學(xué)作用[5]等多種機制共同作用的結(jié)果，其中力學(xué)作用被普遍認為是空蝕破壞的主要原因，認為空化泡在潰滅和再生長的過程中會壓縮空化泡周圍的流體介質(zhì)，產(chǎn)生強度極大的壓力波，從而使壁面材料發(fā)生變形[6]。化學(xué)腐蝕理論認為，空化和腐蝕作用往往同步發(fā)生，互相促進?？栈瘯r的力學(xué)作用清除了材料表面的氧化膜，使得腐蝕作用繼續(xù)進行，而腐蝕也能使空化的力學(xué)作用更為集中[7]。此外，電化學(xué)作用能夠使材料內(nèi)部及其周圍形成電偶電池[8]，熱學(xué)作用則能在局部產(chǎn)生上萬度的高溫，這些因素在特定條件下都能對表面材料產(chǎn)生破壞作用[9]。

針對西藏地區(qū)的水利水電工程而言，在水力空化中，受高海拔復(fù)雜且難以控制與預(yù)測的環(huán)境因素的影響，空泡伴隨著液體流動經(jīng)歷初生、發(fā)展和潰滅的過程[10]。潰滅振蕩時產(chǎn)生局部高速微射流與壓力沖擊波，使水力機械、泄洪設(shè)施、冷卻系統(tǒng)等水力設(shè)施局部溫度過高，從而使其過流部件表面產(chǎn)生空化空蝕破壞，帶來巨大的經(jīng)濟損失，對水利設(shè)施造成運行效率降低、使用壽命縮短等問題，甚至可能導(dǎo)致嚴重事故[10]。因此，深入了解空化泡的動力學(xué)規(guī)律、研究空泡現(xiàn)象具有重要的工程價值和學(xué)術(shù)意義。

2 空化泡動力學(xué)研究方法及國內(nèi)外研究現(xiàn)狀

2.1 研究歷程

由于理論研究較大的局限性，目前關(guān)于空化的研究方法主要包含實驗法與數(shù)值模擬2種。2種主流方法的優(yōu)缺點對比如表1所示。

表1 實驗法與數(shù)值模擬法優(yōu)缺點對比

國內(nèi)外學(xué)者已對多相流系統(tǒng)中氣泡的運動特性開展了一系列實驗研究。然而，受氣泡運動過程的非線性、復(fù)雜多變且難以控制的環(huán)境、液態(tài)金屬的透明程度等因素的影響，難以準確分析氣泡的運動規(guī)律，難以在光學(xué)測量技術(shù)下獲取準確的氣泡運動參數(shù)[11]。另外，由于實驗法缺點的局限性，因此許多研究人員利用模擬軟件解決實驗方法中遇到的困難。數(shù)值模擬法具有突出的優(yōu)勢[12]，受到了廣大研究人員的青睞。最初的空泡潰滅過程模擬主要采用有限差分法和MAC方法，然而，這2種方法在空泡模擬方面存在一些限制，求解過程復(fù)雜且精度不高[13]。為了克服這些問題，研究者逐漸轉(zhuǎn)向使用流體體積法(Volume Of Fluid，VOF)和邊界元法(Boundary Element Method，BEM)進行數(shù)值模擬。在過去，閉源商業(yè)軟件Fluent被廣泛應(yīng)用于空泡模擬，然而，隨著開源軟件的發(fā)展，OpenFOAM和OpenFVM等開源軟件逐漸嶄露頭角[13]。21世紀初，格子波爾茲曼法(Lattice Boltzmann Method，LBM)亦成為空化領(lǐng)域一種重要數(shù)值模擬方法。數(shù)值計算的流程見圖1。

圖1 數(shù)值計算流程

2.2 研究現(xiàn)狀

目前，研究者常常采用實驗、理論和模擬相結(jié)合的方式進行研究。在數(shù)值模擬方面，常用的方法包括流體體積法、水平集法(Level Set，LS)、界面追蹤法(Finite Volume Method，F(xiàn)TM)、邊界元法、格子玻爾茲曼法以及有限體積法(Finite Volume Method，F(xiàn)VM)。

由于空化的復(fù)雜性，現(xiàn)有的理論模型在實際工程應(yīng)用中很難準確預(yù)測空化破壞，因此計算流體力學(xué)(CFD)為其提供了有效的方法。上述數(shù)值模擬方法的優(yōu)缺點對比如表2所示。

表2 數(shù)值模擬方法的優(yōu)缺點對比

2.2.1 流體體積法研究現(xiàn)狀

VOF法是一種基于歐拉網(wǎng)格的界面追蹤模擬方法，將氣液兩相視為均相流動，用標(biāo)量體積分數(shù)來區(qū)分2種流體[12]。通過求解質(zhì)量守恒輸運方程，此方法可模擬交界面拓撲結(jié)構(gòu)的變化，使用特殊平流格式(如幾何平流)求解以最小化數(shù)值擴散，其界面捕捉方程：

(6)

Sagar等[14]使用流體體積法研究了考慮相變的近壁區(qū)域空化氣泡的破裂過程，計算結(jié)果與實驗測量一致。Li等[15]使用VOF法研究了旋流中單個氣泡的動力學(xué)。Nguyen等[16]使用基于幾何VOF的數(shù)值模擬方法，研究了平壁、斜壁和自由表面附近氣泡潰滅過程中的射流和壓力載荷等動力學(xué)特性。Annaland等[17]提出了一種三維流體體積方法，可以處理較大的密度粘度比和表面張力系數(shù)。

通過使用Fluent軟件中VOF模型以及非穩(wěn)態(tài)方法，李疆等[18]對近水平壁面的空泡潰滅過程進行了數(shù)值研究。張凌新等[19]通過Rayleigh方程將球形氣泡隨時間變化的解析解成功求出，使用VOF法模擬了單個純氣體泡。夏冬生等[20]使用VOF計算了空泡潰滅(初始半徑=0.1mm)過程。胡影影等[21]通過VOF中Youngs法追蹤了空泡界面，直接求解納維爾-斯托克斯(Navier-Stokes，N-S)方程，對距離固壁不同位置的空泡進行了數(shù)值計算。劉蘭等[22]研究了無邊蒸汽空泡的潰滅過程，結(jié)果表明，蒸汽空泡的最小半徑遠小于氣體空泡，且從最大半徑到最小半徑所需的時間比氣體空泡更長。

2.2.2 水平集法研究現(xiàn)狀

LS法是一種用距離函數(shù)表示交界面的方法，可以在不重構(gòu)界面的情況下高精度地確定曲率和方向。但當(dāng)界面發(fā)生斷裂或合并時，不能滿足質(zhì)量守恒[23]。該方法的目標(biāo)是減少數(shù)值擴散，將界面定義為與界面距離函數(shù)的零水平集，隨局部流體速度變化而變化。其控制方程形式如下：

(7)

(8)

式中，φ為定義的距離函數(shù)；d為點到界面的歐幾里得距離；Ω+表示氣相；Ω-表示液相；Г表示氣液交界面。

Hu等[23]采用LS法校正守恒損失，模擬二維空泡坍縮沖擊波過程；Mark Sussman[24]將LS法和VOF法應(yīng)用于探究不可壓縮蒸汽泡的生長與潰滅。Lauer等[25]研究了沖擊波作用下并排放置的3個空泡，采用守恒的明銳界面模型和LS法，精確跟蹤了空化泡壁和自由液面的流動。Can和Prosperetti[26]基于改進的LS法研究了蒸汽泡的動力學(xué)。Huang等[27]采用非均勻網(wǎng)格的LS法模擬了空泡在剛性壁面附近的生長、潰滅和回彈過程。

2.2.3 邊界元法研究現(xiàn)狀

BEM，又稱為邊界積分法(Boundary Integral Method，BIM)，是一種常用的研究剛性壁面附近空泡振蕩動力學(xué)的方法，假設(shè)流體為無旋、無粘、不可壓縮，流動由拉普拉斯方程控制，速度場可表示為速度勢的梯度，被廣泛應(yīng)用于研究剛性邊界附近的空泡振蕩動力學(xué)行為。其控制方程形式如下：

(9)

式中，G為三維格林函數(shù)；s為氣泡與自由面表面；φ為速度勢。

BEM方法模擬的空化泡動態(tài)特性大致分為2個主要階段，即射流刺穿空化泡前和射流刺穿空化泡后。然而，對于射流刺穿空化泡后的復(fù)雜模擬需要特殊處理和復(fù)雜的輔助函數(shù)。Blak等[28]使用軸對稱構(gòu)型和三維BEM法計算了近剛性壁面和自由表面的氣泡坍塌過程，闡明了影響潰滅腔內(nèi)液體射流的形成和方向以及附近剛性邊界上產(chǎn)生壓力的一些因素。Aziz等[29]研究了偏心距對氣泡形狀、射流形成等的影響，以及波向和水平間距對氣泡動力學(xué)的影響。Wang等[30]基于經(jīng)典的BEM法和渦環(huán)模型研究了氣泡的分裂現(xiàn)象，建立了多渦環(huán)模型，研究了2個環(huán)形氣泡之間的相互作用，得到了與實驗結(jié)果吻合較好的數(shù)值計算結(jié)果。王起棣等[31]采用插值法對網(wǎng)格進行細分，模擬近壁空化特性(考慮表面張力的影響)，結(jié)果表明近壁空泡潰滅過程受表面張力的顯著影響。

2.2.4 格子玻爾茲曼法研究現(xiàn)狀

近年來，LBM成為了模擬多相流的一種新方法。其基于格子氣自動機的發(fā)展，可用來求解Boltzmann方程的離散化方法，其控制方程形式如下：

fi(x+eiδt，t+δt)=Ωi(x，t)+fi(x，t)

(10)

(11)

(12)

式中，x為格點坐標(biāo)；t為時間；ei為粒子離散速度集；fi為粒子分布函數(shù)；δt為離散時間步長；Ωi為碰撞算子。式(6)代表粒子的碰撞，式(7)代表粒子的遷移。

Shan等[32]使用LBM方法對氣泡的潰滅進行了模擬，結(jié)果與理論和實驗吻合較好，并從二維壓力場演化過程出發(fā)研究了氣泡破裂階段不同部位的動力學(xué)特性，討論了第2次塌陷在剛性邊界損傷中的作用，論證了2次塌陷之間的阻礙作用。袁曉龍等[33]采用LBM模型模擬近壁空泡潰滅的過程，以此分析了空泡潰滅過程受不同汽/液黏滯系數(shù)對的影響，不同條件下空泡潰滅的最大射流速度的變化以及最大壓強的變化。胡玉等[34]用LBM法模擬靜水中雙氣泡的上升過程，研究發(fā)現(xiàn)，當(dāng)氣泡直徑相同時，上層氣泡的變形程度與單個氣泡相似，而下層氣泡則明顯受到前一氣泡尾跡的變形。

2.2.5 有限體積法及界面追蹤法研究現(xiàn)狀

FVM及FTM相似，都是采用一組標(biāo)記點，并通過點來擬合出一組線條，用于標(biāo)記整個氣泡界面的方法。該方法的核心在于如何標(biāo)記點的位置，不同的學(xué)者采用不同的方法。但無論如何選擇方法，為了提高計算精度，都需要提高標(biāo)記點的數(shù)目，因此該方法的計算量較大，且無法模擬氣泡破裂等行為。

Popinet等[35]基于FTM研究了粘性對固體邊界附近氣泡潰滅射流形成的影響。Liu等[36]在FVM法的基礎(chǔ)上，利用固定網(wǎng)格上的交錯網(wǎng)格直接求解無粘不可壓縮假設(shè)的N-S方程，并用FTM模擬了固體附近環(huán)形氣泡的動力學(xué)，分析了不同壁面距離下水層的變化過程、飛濺流和徑向流的發(fā)展過程、環(huán)形氣泡的分裂現(xiàn)象以及固體壁面中心壓力的變化趨勢。Johnsen等[37]使用高階精確的激波和FTM模擬了平面剛性表面附近和自由場中氣泡的潰滅，通過比較潰滅時間、氣泡位移、界面速度和表面壓力與驅(qū)動壓力比和氣泡與壁面距離的函數(shù)，與現(xiàn)有的理論和實驗結(jié)果展現(xiàn)出良好的一致性。雷杰等[38]用FTM方法模擬同軸雙氣泡，發(fā)現(xiàn)雙氣泡的上升速度大于單氣泡的上升速度，合并后的雙氣泡的速度與直徑相同的單氣泡的速度相同，下層氣泡的速度隨著氣泡間距的減小而增大。

3 結(jié)論與展望

目前應(yīng)將提高氣泡測量方法的時間分辨率與空間分辨率、擴大氣泡測量范圍、增強氣泡在復(fù)雜流場環(huán)境中的適應(yīng)性作為重點。傳統(tǒng)的測量方法只能獲得氣泡的二維信息，而不能獲得氣泡的三維空間位置。目前數(shù)值模擬僅考慮了重力、表面張力、浮力等因素，對浮選過程中湍流特性與渦旋尺度等微觀方面的研究較少。部分模型的預(yù)測結(jié)果與實驗結(jié)果不吻合，難以準確計算浮選設(shè)備的內(nèi)力場，這些問題有待進一步研究。

未來的研究可以更加關(guān)注氣液固三相流體的探索；進一步深入理解氣泡形成和消失的機制；引入更精細的數(shù)值模擬方法并結(jié)合實驗驗證，以提高對氣泡行為的理解；實現(xiàn)氣泡的三維測量；研究氣泡的形態(tài)、運動特性以及其與固體表面的相互作用；探索氣泡在化工、生物醫(yī)學(xué)和能源等不同領(lǐng)域中的應(yīng)用潛力；改進研究方法和技術(shù)，提高研究的準確性和可靠性。