【裝備理論與裝備技術(shù)】

非穩(wěn)態(tài)超空泡流動的數(shù)值模擬

趙怡，劉平安，苗成林

(哈爾濱工程大學(xué) 航天工程系，哈爾濱150001)

摘要：采用理論分析和數(shù)值模擬相結(jié)合的方法，對水下超高速運動體自然空化形成超空泡形態(tài)問題展開研究。分析了水下超高速運動體在不同速度下超空泡形態(tài)的變化規(guī)律，比較不同速度下空泡脫落計算結(jié)果的異同。通過數(shù)值模擬，驗證了LES數(shù)值模擬與試驗結(jié)果的關(guān)系，以及評估非穩(wěn)態(tài)空泡流的計算模擬能力。

關(guān)鍵詞：多相流模型；超空泡；LES模型

收稿日期：2014-08-29

基金項目：中央高?；究蒲谢?HEUCFD1404)

作者簡介：趙怡(1990—)，女，碩士研究生，主要從事固體火箭發(fā)動機燃燒、流動及內(nèi)彈道研究。

doi:10.11809/scbgxb2015.01.018

中圖分類號：TV131.3+2

文章編號：1006-0707(2015)01-0064-04

本文引用格式：趙怡，劉平安，苗成林.非穩(wěn)態(tài)超空泡流動的數(shù)值模擬[J].四川兵工學(xué)報，2015(1):64-67.

Citationformat:ZHAOYi,LIUPing-an,MIAOCheng-lin.NumericalSimulationofUnsteadyCavitatingFlow[J].JournalofSichuanOrdnance,2015(1):64-67.

NumericalSimulationofUnsteadyCavitatingFlow

ZHAOYi,LIUPing-an,MIAOCheng-lin

(CollegeofAerospaceandCivilEngineering,HarbinEngineeringUniversity,Harbin150001,China)

Abstract:The supercavity form which is formed by the natural cavitation of the ultra-high-speed underwater moving body was researched by using the method of the combination of theoretical analysis and numerical simulation, which aims to study the motion state of the supercavity form and to compare the similarities and differences of the vacuoles off calculations at different speeds. The relationship between LES simulation and experimental results was verified by using the method of numerical simulations to assess the capacity of numerical simulation of unsteady cavitating flow.

Keywords:multiphaseflowmodel;supercavitation;LESmodel;

發(fā)展水中高速兵器的關(guān)鍵在于大幅減少兵器的運動阻力，水中運動阻力可以分為壓差阻力與黏性阻力，為減小壓差阻力可以優(yōu)化運動體外形[1]。因黏性阻力占總阻力的大部分，所以減小黏性阻力至關(guān)重要，其方法有很多，如在運動介質(zhì)中增加高分子溶劑、吸出或吹出邊界層、采用柔性壁面、利用減阻涂層、安裝順流肋條、調(diào)節(jié)壁面溫度等，但這些方法十分有限，目前水下減租最有效的方法就是采用超空泡技術(shù)[2]。

對水下運動體來說，人工通氣超空泡運動的減阻，在數(shù)量上要低于自然汽化超空泡運動的減阻。隨著超空泡技術(shù)的不斷發(fā)展和完善，可大大提高潛水艇、艦只等的航行速度[3]。我國對模擬空泡多相流模型已做了許多研究，探討了軸對稱體超空泡流動問題[4]。

1數(shù)值模擬

空泡流動主要為湍流運動[5]，模擬湍流運動采用的是LES模型[6]，壓力梯度采用PRESTO!格式離散，動量方程的差分格式采用中心差分格式，壓力-速度耦合采用PISO算法。研究基于Mixture多相流模型，并結(jié)合輸運方程，用簡單幾何外形的水下航行體模型，數(shù)值模擬非穩(wěn)態(tài)超空泡流動現(xiàn)象。

2非穩(wěn)態(tài)超空泡數(shù)值研究

2.1主要內(nèi)容

采用Mixture兩相流模型模擬軸對稱運動體在不同速度下的超空泡流。目的是分析在不同速度的情況下，超空泡形態(tài)的變化和超空泡內(nèi)部參數(shù)的變化規(guī)律，進而分析空泡脫落現(xiàn)象。空化泡是液體流動在一定壓力和一定溫度下發(fā)生的汽化或氣化現(xiàn)象，前者稱為含汽型空化泡，后者稱為含氣型空化泡，并統(tǒng)稱為空化泡，常簡稱為空化或空泡。空化是水流或液流在低壓處突然發(fā)生的空泡(汽化或氣化空穴)現(xiàn)象。

空化數(shù)

(1)

式中：p∞為參考流體的絕對壓力；pv為液體的汽化壓強；U∞為參考流體的運動速度；ρ為液體密度。

流體物性參數(shù)：水的汽化壓強pv=3 540Pa，參考流體的絕對壓力p∞=101 325Pa，液體密度ρ=1 000kg/m3。運動體在不同速度下通過式(1)可計算出空化數(shù)，通常，自然超空泡流的產(chǎn)生要求自然汽化空化數(shù)小于0.1。

2.2建模

本文主要研究的是細錐形航行體的超空泡流動，因航行體為軸對稱運動體，所以把模型簡化為二維模型。彈體長l=1.2m，直徑d=0.1m，采用如圖1所示矩形結(jié)構(gòu)化網(wǎng)格對流場計算域進行劃分，并按需求對彈體進行局部網(wǎng)格加密，情況如圖2，外流場長為12m。

圖1　運動體整體尺寸

圖2　模型局部網(wǎng)格劃分情況

2.3數(shù)值計算

速度入口參數(shù)設(shè)置如表1所示。

表1　速度入口參數(shù)設(shè)置

2.4數(shù)據(jù)處理

2.4.1開始階段

1) 算例1

圖3為運動體速度為80m/s時，空泡形成初期，不同時刻的外流場密度云圖。觀察圖片發(fā)現(xiàn)，空泡從形成開始逐漸變長，隨時間的增加達到一定長度后開始脫落，脫落從空泡尾部較快完成。

2) 算例2

圖4展示的是運動體速度為100m/s時，空泡形成初期，不同時刻的外流場密度云圖。觀察圖片發(fā)現(xiàn)，與上一算例相同空泡也出現(xiàn)拖尾現(xiàn)象，脫落從空泡尾部較快完成第一次脫落。但是緊接著出現(xiàn)了二次脫落，使空泡又一次變短。

圖4　空泡脫落過程云圖

3) 算例3

經(jīng)觀察發(fā)現(xiàn)，速度為130m/s時的空泡脫落現(xiàn)象比1、2算例更明顯(見圖5)。仔細觀察可以發(fā)現(xiàn)圖5(c)也發(fā)生了二次脫落，且脫落時間更短才使得脫落現(xiàn)象明顯。本次模擬與實驗[9]結(jié)果吻合較好，說明大渦模擬的方法有效。

圖5　空泡脫落過程云圖

把圖3(e)、圖4(e)和圖5(e)進行比較，可以發(fā)現(xiàn)速度越大，脫落后的空泡長度越短，即在空泡上發(fā)生斷裂的部位越靠前。在這3個算例中，速度越大脫落得越快，實驗中加速時尾部空泡脫落實照如圖6所示。

圖6　實驗中加速時尾部空泡脫落實照

2.4.2穩(wěn)定階段

以上3個算例在接下來的一段時間內(nèi)均無明顯的脫落現(xiàn)象，雖然空泡的長度始終不穩(wěn)定，但變化不大，可近似穩(wěn)定狀態(tài)。在彈身后方軸線上和非軸線上分別取兩點1、2，計算中監(jiān)測兩點的壓力和速度(見圖7～圖10)。

圖7　點1軸向速度-時間圖

圖8　點1空化數(shù)-時間圖

圖9　點2軸向速度-時間圖

圖10　點2空化數(shù)-時間圖

經(jīng)過觀察發(fā)現(xiàn)，水和空氣在流體黏性、壁面效應(yīng)等因素的作用下不斷地進行無規(guī)則運動，各處的壓力和速度隨著時間而不斷地發(fā)生改變，速度和壓力沒有一定的變化規(guī)律，且波動較大。

3結(jié)論

采用CFD軟件FLUENT，基于Mixture多相流模型，湍流輸運方程和各相的連續(xù)性方程，數(shù)值模擬了超空泡的復(fù)雜流動現(xiàn)象。

本文就非穩(wěn)態(tài)超空泡的流動情況進行模擬，用大渦模擬的方法，通過不斷地修改嘗試， 1.2m的彈體模型的數(shù)值模擬計算結(jié)果與實驗數(shù)據(jù)吻合較好。經(jīng)比較發(fā)現(xiàn)，空泡脫落現(xiàn)象只明顯地發(fā)生在初始階段，且速度越大脫落的時間越短，也越明顯。脫落之后，空泡又重新生長，之后保持較穩(wěn)定的狀態(tài)，沒有周期的規(guī)律可循。若進一步降低彈體速度可能不會形成超空泡。3個算例的空化數(shù)均小于0.03，完全可以形成超空泡，且在形成過程中初始階段都有空泡脫落現(xiàn)象。而在實際水下發(fā)射魚雷，周圍環(huán)境復(fù)雜，稍有一點擾動對彈體就有很大的影響，使水下航行體不斷地加速減速。所以實際上空泡的脫落現(xiàn)象時刻存在。這一論題的研究也十分重要。

參考文獻：

[1]陳兢.新概念武器——超空泡水下高速武器[J].飛航導(dǎo)彈,2004(10):34-37.

[2]王獻孚.空化泡和超空化泡流動理論及應(yīng)用[M].北京：國防工業(yè)出版，2009：1-195.

[3]趙衛(wèi).超空泡高速魚雷技術(shù)綜合分析[D].哈爾濱：哈爾濱工程大學(xué)，2005.

[4]楊洪瀾.非定常超空泡繞流研究[D].哈爾濱：哈爾濱工業(yè)大學(xué),2007.

[5]褚學(xué)森,王志.自然空化流動數(shù)值模擬中參數(shù)取值影響的研究[J].船舶力學(xué),2007,11(1):32-39.

[6]RodiW.ComparisonofLESandRANScalculationsoftheflowaroundbluffbodies[J].JournalofWindEngineeringandIndustrialAerodynamics,1997,69-71:55-75.

[7]KolmogorovAN.ThelocalStructureofTurbulenceinIncompressibleViscousFluidforveryLargeReynoldsNumbers[J].Dokl.Akad.Nauk.SSSR,1941,30(4):299-303.

[8]ShinBR,IkohagiT.Anumericalstudyofunsteadycavitatingflows[C]//3rdInternationalSymposiumonCavitation.Grenoble,France,1998.

[9]熊永亮.水下高速航行體超空泡減阻特性研究[D].哈爾濱：哈爾濱工程大學(xué)，2008.

[10]郭小飛，劉平安， 時文，等.水下超高速運動體人工通氣流動過程研究[J].四川兵工學(xué)報，2014(6):6-8.

(責任編輯周江川)