余德磊,王 聰,何超杰

(哈爾濱工業(yè)大學 航天學院,哈爾濱 150001)

并聯(lián)入水是指兩個或多個運動體沿空間平行線同時或在極短時間間隔內(nèi)穿越自由水面進入水中的過程.在實際作戰(zhàn)中，往往需要在短時間內(nèi)連續(xù)發(fā)射多枚魚雷和射彈，對敵方單位造成有效打擊，該過程涉及到了典型的并聯(lián)入水問題。

國內(nèi)外對單個運動體入水開展了大量的研究，包括空泡形態(tài)特性、受載特性和彈道特性.Logvinovich[1]基于能量守恒原理得到了空泡獨立膨脹原理，為后來的空泡形態(tài)發(fā)展研究提供了理論基礎；May等[2]開展大量入水實驗工作，對鋼球入水初期的流動特點與入水載荷系數(shù)進行詳細研究；Shi等[3]開展了大量子彈高速垂直入水的試驗研究，得到了入水噴濺、回射流等非定常流動特性，同時發(fā)現(xiàn)子彈入水后的彈道偏移與入水深度有關(guān)；Aristoff等[4]進行了系列疏水球入水試驗，得到了不同入水速度下的空腔形態(tài)，并基于空泡壓力平衡理論和勢流理論成功預測了不同入水條件下的空腔演化；Guo等[5]開展進行高速水平入水試驗，研究了不同頭型、不同入水速度運動體的入水空泡形態(tài)與運動特性，并基于Rayleigh-Besant方程預測了入水過程中空泡的最大直徑；張偉等[6]建立了平頭運動體入水的空泡形態(tài)和彈道預測模型；王柏秋等[7]基于動網(wǎng)格技術(shù)進行超空泡射彈研究，得到了超空泡射彈的阻力系數(shù)變化規(guī)律；馬慶鵬等[8]基于VOF多相流模型研究了錐頭柱體高速入水過程，結(jié)果表明運動體入水速度越大航行體頭部壓力越高，隨著入水深度增加空化泡內(nèi)空化現(xiàn)象明顯；何乾坤等[9-15]對采用數(shù)值方法，對超空泡射彈的尾拍特性進行了系統(tǒng)研究。

目前國內(nèi)外對并聯(lián)入水的研究較少。王志東等[16]通過數(shù)值模擬研究了并列航行體的空泡形態(tài)及減阻特性；路麗睿等[17]開展了回轉(zhuǎn)體低速并聯(lián)入水過程的運動特性試驗，研究了入水速度對空泡及運動特性的影響，結(jié)果表明兩回轉(zhuǎn)體的平均偏轉(zhuǎn)角速度隨著入水速度的增大而增大；盧佳興等[18]開展了雙圓柱體低速并聯(lián)入水的試驗研究，結(jié)果表明入水空泡整體呈現(xiàn)良好的鏡面對稱特征,而圓柱體內(nèi)外側(cè)空泡存在明顯的非對稱性；宋武超等[19]基于勢流理論提出了回轉(zhuǎn)體低速并聯(lián)入水過程空泡形態(tài)發(fā)展的預測方法。綜上所述，對入水問題的研究多針對單體入水，涉及高速并聯(lián)入水的研究非常少，且未考慮空化現(xiàn)象對運動體并聯(lián)入水流體和運動的影響。

本文采用數(shù)值研究方法，將數(shù)值計算結(jié)果與文獻[6]中空泡半徑的預測公式進行對比，驗證了本文數(shù)值模擬方法的有效性，在此基礎上開展不同入水速度、不同初始凈距和不同橫流速度對回轉(zhuǎn)體并聯(lián)入水過程的數(shù)值模擬，研究了上述參數(shù)對并聯(lián)回轉(zhuǎn)體空泡特征尺寸、側(cè)向及偏航運動影響規(guī)律。

1 數(shù)值計算方法

1.1 控制方程及其求解

本文基于有限體積法對雷諾時均的納維-斯托克斯方程進行離散，引入realizablek-ε湍流模型[20]、VOF多相流模型和Schnerr and Sauer空化模型[21]來描述湍流流動、各相界面及空化現(xiàn)象，并利用重疊網(wǎng)格技術(shù)處理運動邊界，對并聯(lián)入水問題進行數(shù)值計算。

混合介質(zhì)的連續(xù)性方程為

(1)

其中，i=1, 2, 3.下同.動量方程為

式中:ρm=αlρl+αgρg+αvρv，μm=αlμl+αgμg+αvμv,其中αl、αg和αv分別為水、空氣和水蒸氣的體積分數(shù)，ρl、ρg和ρv為三相的密度，μl、μg和μv為三相的動力黏度;μt=ρmCμk2/ε為湍流黏性系數(shù)，其中Cμ為經(jīng)驗常數(shù)，k為湍動能，ε為湍動耗散率;ui、uj為速度分量;xi、xj為位移分量;Cμ由下式確定：

(2)

由于realizablek-ε模型適用于大雷諾數(shù)的流動，因此湍流模型采用realizablek-ε模型，其湍動能和湍動耗散率的輸運方程為:

Gk+Gb-ρε-YM+Sk

(3)

(4)

本文采用Schnerr and Sauer空化模型描述空化現(xiàn)象。水蒸氣相輸運方程為

(5)

式中:RB=1×10-6m為氣核半徑，αnuc=5×10-4為不可凝結(jié)氣體體積分數(shù)，p為遠場壓力，pv為飽和蒸氣壓，F(xiàn)vap=50,Fcond=0.001。

1.2 方法有效性驗證

文獻[6]中給出了同實驗結(jié)果符合良好的空泡形態(tài)模型，其表達式為

(6)

式中:R為空泡半徑，R0為頭部半徑，z為位移，z0為初始位移，σ0、σ分別為初始空化數(shù)和空化數(shù)，Cd=C0(1+σ),0.82≤C0≤0.83，N為經(jīng)驗系數(shù)，取為2。

為驗證計算方法的有效性，采用流體仿真軟件STARCCM+12.0對頭部半徑為5 mm，長細比λ=6的鋁質(zhì)平頭柱體以98.7 m/s入水過程進行數(shù)值計算.將計算結(jié)果和式(6)進行對比，結(jié)果如圖1所示?？梢钥闯鰞烧咻^吻合，說明本文數(shù)值計算方法是有效的。

圖1 空泡形態(tài)對比

利用本文問題的對稱性，僅取半邊模型進行計算。兩回轉(zhuǎn)體之間的重疊區(qū)域與背景區(qū)域尺寸如圖2所示。

圖2 重疊網(wǎng)格區(qū)域與背景域幾何示意

本文選取了3種網(wǎng)格，3種網(wǎng)格的參數(shù)見表1。

表1 不同密度和數(shù)量的網(wǎng)格

設置相同的時間步長1×10-6s進行計算.數(shù)值計算結(jié)果表明3種網(wǎng)格的阻力系數(shù)計算結(jié)果差異不大；而網(wǎng)格1與網(wǎng)格2和網(wǎng)格3在縱向速度和角速度的計算上較大偏差。本文為減少計算成本，綜合比較后決定選取網(wǎng)格2進行后續(xù)計算。

2 數(shù)值計算結(jié)果與分析

2.1 入水速度對并聯(lián)入水空泡及運動特性影響

本文中的參考截面為對稱平面，根據(jù)對稱平面的特征提取可以獲得空泡外側(cè)極徑rw、內(nèi)側(cè)極徑rn、空泡長度L、限制長度L1、L′1和噴濺高度h等空泡特征，如圖3所示。其中d為并聯(lián)凈距,D為回轉(zhuǎn)體直徑。由于對稱性，如非特別說明，本文僅考察左側(cè)回轉(zhuǎn)體得空泡于運動特性。

圖3 參考截面空泡形態(tài)參數(shù)示意

為研究入水初速度對回轉(zhuǎn)體并聯(lián)入水過程的影響，本文取初始凈距為0.4D，入水初速度分別為69.67、79.64和89.48 m/s(即初始空化數(shù)分別為0.040、0.031、0.024)的3種情況進行數(shù)值計算，并比較入水初速度對流場和運動特性的影響。

圖4、5分別為不同入水速度下外側(cè)空泡形態(tài)和限制長度演化過程。從圖中可以發(fā)現(xiàn)，較大的入水初始速度為回轉(zhuǎn)體頭部的液體提供了較大的初始動能，相應的液體徑向運動的初始速度較大，同一量綱一的時刻空泡的外側(cè)極徑和限制長度都隨入水初速度增大而增大。

圖4 入水初速度對外側(cè)空泡形態(tài)的影響

圖5 入水初速度對限制長度的影響

圖6、7分別給出了不同入水速度下的回轉(zhuǎn)體頭部壓力及側(cè)向壓差分布。可以看到入水速度對回轉(zhuǎn)體頭部的壓力分布影響較大，隨著速度增大，壓力峰值增大。對比圖7的側(cè)向壓差發(fā)現(xiàn)，較小的入水速度時回轉(zhuǎn)體初期分布特征明顯不同，其主要原因為空泡內(nèi)部的水蒸氣分布差異較大，此時存在明顯的回射流作用.在入水一段時間后，側(cè)向壓差的峰值隨著入水速度的增大而增大。

圖6 入水速度對頭部壓力分布的影響

圖7 入水速度對側(cè)向壓差分布的影響

圖8為不同入水速度下回轉(zhuǎn)體的側(cè)向位移與偏航角變化規(guī)律.從圖中可以看到速度增大，回轉(zhuǎn)體的側(cè)向位移和偏航角均呈增大趨勢。該趨勢與較大的頭部壓力密切相關(guān).入水初速度越大，頭部壓力造成的偏航力矩和附加側(cè)力越大，從而對側(cè)向和偏航運動的促進作用越強。另外，隨著空泡內(nèi)側(cè)的水蒸氣含量增大，回轉(zhuǎn)體受到的側(cè)向壓差增加，進一步促進了側(cè)向運動。

圖8 入水速度對回轉(zhuǎn)體運動特性的影響

2.2 初始凈距對并聯(lián)入水空泡及運動特性影響

為說明初始凈距對無橫流情況下回轉(zhuǎn)體并聯(lián)入水過程的影響，本文取入水初速度為98.7 m/s(即初始空化數(shù)為0.02)，初始凈距分別為0.2D、0.4D、0.6D和1.0D這4種情況進行數(shù)值計算，比較其對流場和運動特性的影響。

圖9、10分別為不同初始凈距下回轉(zhuǎn)體的空泡形態(tài)和限制長度演化過程?；讵毩⑴蛎浽恚鈧?cè)空泡的形態(tài)基本相同.內(nèi)側(cè)空泡的限制長度L′1差異明顯。從圖10中可以看到限制長度L′1隨初始凈距增大而增大，這是由于較小的初始間距造成頭部的速度駐點靠近內(nèi)側(cè)，回轉(zhuǎn)體內(nèi)側(cè)流體的分離速度較大，劇烈的空化作用促使內(nèi)側(cè)空泡在較小區(qū)域內(nèi)快速融合，限制長度反而越小。

圖9 初始凈距對外側(cè)空泡形態(tài)的影響

圖10 初始凈距對限制長度的影響

圖11、12分別為初始凈距對回轉(zhuǎn)體頭部壓力及側(cè)向壓差分布的影響。可以發(fā)現(xiàn)初始凈距對回轉(zhuǎn)體頭部壓力峰值影響不大，但入水初期較小間距的頭部壓力分布不對稱性明顯，且壓力峰值越靠近內(nèi)側(cè).通過側(cè)向壓差分布發(fā)現(xiàn)，在入水初期回轉(zhuǎn)體的側(cè)向壓力存在短暫振蕩，這與劇烈的空化過程有關(guān)。一段時間后初始間距越小側(cè)向壓差峰值越靠近回轉(zhuǎn)體頭部，這與上文中的較小限制長度相吻合。

圖11 初始凈距對頭部壓力分布的影響

圖12 初始凈距對側(cè)向壓差分布的影響

圖13為初始凈距對側(cè)向運動和偏航運動的影響規(guī)律。從圖中可以發(fā)現(xiàn)，回轉(zhuǎn)體的側(cè)向位移和偏航角隨初始凈距增大呈現(xiàn)先增大后減小的趨勢。這是由于初始凈距較小時側(cè)向壓差顯著地抑制了頭部壓力對側(cè)向運動的促進作用，從而起到限制側(cè)向運動的作用。當初始凈距足夠大時，隨著初始凈距的減小，頭部壓力造成的附加側(cè)力和偏航力矩越大，此時初始凈距越小，側(cè)向位移和偏航角越大；當初始凈距小于某一臨界值時，速度駐點幾乎不再向內(nèi)側(cè)偏移，但內(nèi)側(cè)靠近頭部由于水蒸氣的存在而使得內(nèi)側(cè)靠近頭部壓力較大，于是回轉(zhuǎn)體會受到正方向力矩，使偏航力矩減小，進而使附加側(cè)力減小。

圖13 初始凈距對回轉(zhuǎn)體運動特性的影響

2.3 橫流速度對并聯(lián)入水空泡及運動特性影響

為研究橫流對回轉(zhuǎn)體并聯(lián)入水影響，本文基于相同的入水條件(v0=98.7m/s；d/D=0.6)，進行橫流速度分別為5 m/s和20 m/s的入水數(shù)值仿真。

圖14為有橫流情況下單體入水和并聯(lián)入水的參考截面處空泡形態(tài)對比。可以看出，橫流作用下單回轉(zhuǎn)體的迎流和背流側(cè)空泡形態(tài)差異明顯，且隨著橫流速度增大，空泡的徑向尺寸差異增大；在并聯(lián)入水狀態(tài)下，迎流與背流回轉(zhuǎn)體的外側(cè)空泡徑向尺寸差異與單回轉(zhuǎn)體基本相同，由于空泡開口較大，較小的壓差作用導致迎流回轉(zhuǎn)體外側(cè)空泡較晚發(fā)生回卷與閉合，因此長度略長。背流回轉(zhuǎn)體外側(cè)空泡獨立膨脹，但迎流方向的回轉(zhuǎn)體的阻礙與兩回轉(zhuǎn)體中間液體的能量轉(zhuǎn)化，導致其空泡極徑減?。挥髋c背流回轉(zhuǎn)體的內(nèi)側(cè)空泡主要受到徑向間距限制，橫流作用下迎流回轉(zhuǎn)體內(nèi)側(cè)空泡略大于背流回轉(zhuǎn)體，最大極徑接近d+D/2；隨著橫流速度增大，迎流與背流回轉(zhuǎn)體的外側(cè)空泡尺寸與單回轉(zhuǎn)體的差異增大，并聯(lián)運動體的內(nèi)側(cè)空泡尺寸差異增大。

圖14 有橫流情況下單體和并聯(lián)入水過程參考截面空泡形態(tài)對比

圖15給出了兩種橫流速度下的水蒸氣相分布?？梢钥吹捷^大的橫流速度導致空泡不對稱增強，相應的內(nèi)部水蒸氣分布更靠近橫流方向。

圖15 有橫流情況下并聯(lián)入水蒸氣相分布

圖16為橫流對并聯(lián)入水運動體的側(cè)向運動和偏航運動的影響。當橫流速度較小時，空泡內(nèi)水蒸氣集中在兩回轉(zhuǎn)體中間區(qū)域，因此較大的側(cè)向壓差導致迎流回轉(zhuǎn)體向迎流方向運動，這與單獨入水時相反；當橫流速度較大時，橫流對回轉(zhuǎn)體的沖擊作用占主導地位，運動方向與有橫流情況下單回轉(zhuǎn)體的運動方向相同，但由于側(cè)向壓差的差異，單體入水的運動較大。對于背流回轉(zhuǎn)體，在頭部繞流產(chǎn)生的側(cè)向力作用下向外側(cè)移動，且隨著橫流速度增大而增大。

圖16 橫流速度對回轉(zhuǎn)體運動特性的影響

對于偏航運動，較小的橫流速度形成了較為均勻的頭部壓力分布，因此迎流回轉(zhuǎn)體偏航角基本為零。背流回轉(zhuǎn)體頭部壓力分布在橫流作用下更加不對稱，且稀疏的水蒸氣分布導致該處側(cè)向壓差較大，形成負方向的偏航力矩，在相反方向的力矩作用下出現(xiàn)兩回轉(zhuǎn)體頭部靠近而尾部遠離的運動狀態(tài)，這與無橫流情況下并聯(lián)入水過程中回轉(zhuǎn)體的相對姿態(tài)相似；橫流速度較大時，迎流回轉(zhuǎn)體尾部受到橫流的側(cè)向沖擊力遠大于側(cè)向壓差力，背流回轉(zhuǎn)體下部受到內(nèi)側(cè)流體的較大沖擊作用，在上述機理下兩回轉(zhuǎn)體分別產(chǎn)生負/正方向力矩，形成了兩回轉(zhuǎn)體頭部遠離而尾部靠近的運動狀態(tài)，這與無橫流情況下并聯(lián)入水過程中回轉(zhuǎn)體的相對姿態(tài)相反。

3 結(jié) 論

1)高速并聯(lián)入水運動時，兩運動體之間的空化程度與入水初始速度正相關(guān)，且同一量綱一的時刻空泡的外側(cè)極徑和限制長度增大，回轉(zhuǎn)體量綱一的側(cè)向位移和偏航角越大。

2)隨著初始凈距減小，同一量綱一的時刻空泡的限制長度越小，回轉(zhuǎn)體的量綱一的側(cè)向位移和偏航角先增大后減小。

3)小橫流作用下，迎流與背流回轉(zhuǎn)體的外側(cè)空泡徑向尺寸與單體基本相同，而長度略大，隨著橫流速度增大，并聯(lián)入水狀態(tài)的空泡尺寸差異同單體相比增大；內(nèi)側(cè)空泡的徑向尺寸差異較大，其中迎流回轉(zhuǎn)體的內(nèi)側(cè)空泡極徑較大，且隨著橫流速度增大極徑的最大值始終維持在d+D/2左右。

4)橫流速度較小時，并聯(lián)入水回轉(zhuǎn)體頭部靠近而尾部遠離；橫流速度較大時，兩回轉(zhuǎn)體頭部遠離而尾部靠近。