圓形舭部附近氣泡與自由面非線性耦合研究

田昭麗， 劉云龍， 高利峰

(1.哈爾濱工程大學(xué)，哈爾濱　150001; 2 武漢第二船舶設(shè)計研究所，武漢　430205;3. 上海外高橋造船海洋工程設(shè)計有限公司，上海　200137)

?

圓形舭部附近氣泡與自由面非線性耦合研究

田昭麗1,2， 劉云龍1， 高利峰3

(1.哈爾濱工程大學(xué)，哈爾濱150001; 2 武漢第二船舶設(shè)計研究所，武漢430205;3. 上海外高橋造船海洋工程設(shè)計有限公司，上海200137)

摘要：用氣泡動力學(xué)邊界元模型分析圓形舭部附近氣泡與自由面非線性耦合作用?；诓豢蓧嚎s勢流理論建立水下爆炸氣泡動力學(xué)數(shù)值模型；針對自由面與壁面交界點須同時滿足壁面邊界條件、自由面邊界條件的過約束問題，采用節(jié)點拆分方法結(jié)合交界處約束條件對傳統(tǒng)數(shù)值模型進行改進；用所建數(shù)值模型對已有研究成果進行拓展，分析攻角、爆距對氣泡與自由面動態(tài)特性影響規(guī)律，初步討論浮力影響，以期為艦船抗沖擊設(shè)計、研究提供參考。

關(guān)鍵詞：水下爆炸；氣泡動力學(xué)；圓形舭部；自由面；邊界元

水下爆炸氣泡在不同邊界附近運動特性復(fù)雜。對單一邊界附近氣泡運動已有深入研究：如氣泡在自由面附近運動時膨脹階段被自由面吸引，而收縮階段被擊退，并在自由面產(chǎn)生水冢[1-4]；在壁面附近運動時特性則相反[5-7]。同時研究壁面、自由面對附近氣泡運動影響較少。Klaseboer等[8]用鏡像法模擬線性化自由面、水面船舶附近氣泡運動，但只適用水較深、自由面影響較弱情況。初文華等[9]用電火花氣泡實驗對此問題進行分析，發(fā)現(xiàn)自由面與壁面交界處產(chǎn)生明顯的非線性耦合現(xiàn)象。劉云龍等[10]用鏡像法模擬無限大平壁面作用進行分析、研究，但無法推廣到復(fù)雜結(jié)構(gòu)。通常氣泡與艦船尺度相近甚至更大，須考慮舷側(cè)形狀影響。為此，本文用邊界元數(shù)值模型對考慮舭部形狀的復(fù)雜邊界附近氣泡與自由面相互耦合特性進行分析。

本文在此基礎(chǔ)上將軸對稱模型推廣到三維問題分析中，并導(dǎo)出考慮結(jié)構(gòu)運動速度的交界點處理方法。結(jié)合氣泡動力學(xué)模型對文獻[10]的研究進行拓展，分析圓形舭部附近氣泡與自由面耦合動態(tài)特性，以期為艦船抗爆設(shè)計及氣泡與復(fù)雜邊界耦合研究提供參考。

1理論與數(shù)值模型

1.1氣泡動力學(xué)基本方程

以氣泡初始中心為原點建立直角坐標系，y軸指向壁面，見圖1。圖中，h為氣泡中心初始深度；θ為其相對舷側(cè)攻角；d為兩者間初始距離。

據(jù)以往氣泡動力學(xué)相關(guān)研究，水下爆炸周圍流場可假設(shè)為不可壓縮的理想流體[13-15]。據(jù)勢流理論，由自由面與圓舭壁面所圍流場中存在速度勢φ，并滿足邊界積分方程，即

(1)

式中：x，y為場點、源點空間坐標矢量；R=x-y；ζ為立體角；S為流場邊界；n為S指向流場內(nèi)部的單位法向量。

圖1　傾斜壁面附近氣泡與自由面相互作用示意圖Fig.1 Sketch of the interaction between the bubble and the free surface near round bilge

自由面、氣泡表面動力學(xué)邊界條件滿足伯努利方程，即

(2)

式中：P為壓力；g為重力加速度；z為垂向坐標；P為大氣壓，自由面取Patm，氣泡表面取氣泡內(nèi)部壓力Pb，對絕熱假設(shè)的氣泡內(nèi)部氣體，可表示為[16-17]

課程教學(xué)的“整體設(shè)計”是我們按照先進職教觀念提出的一個教學(xué)改革新概念。[1]本文立足于職業(yè)院校室內(nèi)設(shè)計技術(shù)專業(yè)教學(xué)，結(jié)合《室內(nèi)效果圖設(shè)計表現(xiàn)》課程對其進行教學(xué)整體設(shè)計。新的課程教學(xué)設(shè)計體現(xiàn)了現(xiàn)代職業(yè)教育的先進教學(xué)理念，突出能力目標，以實際崗位中的項目及具體工作任務(wù)為導(dǎo)向，以能力為本位，學(xué)生為主體，對知識理論實踐一體化，教學(xué)做一體化進行精心設(shè)計，將“工學(xué)結(jié)合”、“校企合作”等職業(yè)教育原則融合在課程教學(xué)中，充分體現(xiàn)職業(yè)院校教學(xué)特點，符合國家對職業(yè)院校人才培養(yǎng)模式的要求。

(3)

式中：V，V0為氣泡當前、初始體積；P0為氣泡初始壓力；γ為氣泡內(nèi)不可冷凝氣體熱容比，對TNT爆轟產(chǎn)物取1.25；Pc為可冷凝氣體飽和蒸汽壓。

1.2三維自由面與壁面交界處理

將邊界離散為單元及節(jié)點，邊界積分方程式(1)可改寫為

GV=HΦ

(4)

式中：G，H為式(1)右端第一項積分與計入立體角項第二項積分影響矩陣；V，Φ為速度勢法向?qū)?shù)及速度勢本身列向量。

因直接求解式(4)存在過約束問題，而自由面與壁面交界的節(jié)點速度勢與法向速度均已知，為此，將所有邊界分為4部分，即剛性壁面Sw、自由面Sf、氣泡表面Sb?？紤]自由面與氣泡表面邊界條件的相似性，合并兩者為Ss。設(shè)自由面、壁面交界節(jié)點為N，由于該處滿足自由面及壁面邊界條件，須進行特殊處理。本文將節(jié)點N拆分為自由面部分Fi及壁面部分Wi，見圖2,并將其從Sw及Ss中分離，則式(4)可表示成分塊形式，即

(5)

圖2　自由面與壁面交界處理示意圖Fig.2 Sketch of the treatment for the intersection between the free surface and body surface

由于Fi及Wi節(jié)點位置完全重合，故ΦFi=ΦWi，且矩陣G及H中Fi與Wi對應(yīng)的行相同，刪去Wi對應(yīng)行，并將已知項與未知項分別變換到等號兩側(cè)后得

(6)

由上式知，影響矩陣中已將交界點對應(yīng)的元素分解為自由面及壁面影響部分分開考慮。當Vw=0即壁面靜止時，上式即退化為文獻[12]的表達式，表明該數(shù)值模型的正確性。

2結(jié)果與討論

2.1無浮力情況

2.1.1基本現(xiàn)象

以γw=1.0,γf=1.0,θ=30°工況為例分析圓舭舷側(cè)附近氣泡與自由面動態(tài)特性。其中強度參數(shù)κ=100，浮力參數(shù)β=0。不同時刻氣泡及自由面形態(tài)見圖3。

圖3　γw=1.0,γf=1.0,θ=30°,圓舭附近氣泡動態(tài)特性Fig.3 Bubble and free surface dynamics near round bilge with γw=1.0,γf=1.0,θ=30°

由圖3看出，水下爆炸氣泡產(chǎn)生后迅速膨脹，自由面在氣泡推動下產(chǎn)生水冢，并在舷側(cè)壁面阻礙下其左上部被擠壓。在t= 0.86時達最大體積，隨后進入坍塌。在坍塌過程中因自由面與氣泡間高壓區(qū)作用，水冢高度進一步增加。壁面對氣泡的Bjerknes力具有水平方向分量，氣泡右上部產(chǎn)生斜射流。隨氣泡坍塌，在自由面排斥下氣泡被自由面擊退，距離增大將自由面作用削弱，從而使射流逐漸向左即壁面吸引方向傾斜。此現(xiàn)象也可解釋為推動射流發(fā)展的高壓區(qū)向右下方移動所致。

2.1.2角度影響

為獲得不同工況參數(shù)對氣泡、自由面動態(tài)特性影響規(guī)律，固定其它參數(shù)，通過改變目標參數(shù)進行對比分析。選γw= 1.0作為典型爆距，對比分析θ=30°,45°,60°,90°。以最大體積及氣泡被射流穿透時為典型時刻，不同攻角的氣泡及自由面形態(tài)見圖4。由圖4看出，對自由面而言，隨攻角增大氣泡對自由面影響逐漸減小。攻角為30°時自由面在氣泡膨脹階段產(chǎn)生明顯水冢，在氣泡坍塌階段迅速收縮，并向遠離舷側(cè)方向傾斜；攻角為90°時自由面僅產(chǎn)生輕微隆起，并在坍塌階段很快消失。對氣泡而言，隨攻角增大自由面影響逐漸減弱，壁面影響逐漸增強。氣泡形態(tài)主要表現(xiàn)在射流形態(tài)上，隨攻角增大氣泡射流形成越晚射流方向越偏向壁面。

圖4　不同攻角下舷側(cè)附近氣泡與自由面相互作用 Fig.4 Comparison between bubble shapes at typical moment with different attack angle

各工況下氣泡體積隨時間變化曲線見圖5。由圖5看出，隨攻角增大氣泡最大半徑及周期均逐漸增大。因攻角較大時氣泡主要受壁面影響，而壁面會對附近氣泡產(chǎn)生延遲作用即延長氣泡周期。而自由面則對附近氣泡作用相反。另外，由于氣泡最大體積時刻在其擾動下產(chǎn)生的自由面水冢含部分動能，因此氣泡不能膨脹到自由場中達到的最大體積，導(dǎo)致自由面附近氣泡最大半徑偏小。

圖5　不同攻角下氣泡體積時歷曲線Fig.5 Comparison between the bubble volume with different attack angle

2.1.3距離影響

分析爆距對氣泡及自由運動影響規(guī)律。選θ= 45°為典型攻角，選爆距γw=0.6,0.8,1.0,1.5進行對比分析。以最大體積及氣泡被射流穿透時為典型時刻，不同爆距下氣泡及自由面形態(tài)見圖6。由圖6看出，θ= 45°時氣泡最大體積，隨爆距增加壁面、自由面對氣泡影響減小，γw=1.5時氣泡基本呈球形，受邊界作用不明顯。而自由面水冢高度隨爆距增大迅速減小。氣泡進入坍塌階段時爆距越小，自由面誘導(dǎo)的射流發(fā)生越早，氣泡被射流穿透時間也越早。然而，由射流方向逐漸向豎直向下變化知，隨爆距增大自由面對氣泡影響的衰減慢于壁面影響。

2.2有浮力情況

實際中艦船遭水下爆炸產(chǎn)生的大尺度氣泡會受到巨大浮力作用。以θ=45°，ε=100，β2=0.25為例分析有浮力時圓形舭部附近氣泡與自由面的相互作用。不同時刻氣泡與自由面形態(tài)見圖7。該圖為典型時刻氣泡與自由面形態(tài)發(fā)展的計算結(jié)果。由圖7看出，氣泡膨脹至最大體積階段時除壁面阻擋作用外，自由面及浮力均未對氣泡發(fā)展產(chǎn)生明顯影響。氣泡進入坍塌階段在水下爆炸大浮力作用下，在壁面、自由面影響尚未明顯發(fā)展前浮力已使氣泡產(chǎn)生向上的射流，且方向基本垂直向上，未受到結(jié)構(gòu)的明顯作用。

圖6　不同爆距下氣泡與自由面形態(tài)對比 Fig.6 Comparison between bubble shapes at typical moment with different initial distance

為深入分析該過程的力學(xué)機理，在各典型時刻計算流場壓力并繪制云圖見圖8。由圖8看出，大浮力時氣泡收縮階段先在氣泡正下部產(chǎn)生明顯高壓區(qū)，在該高壓區(qū)作用下氣泡迅速產(chǎn)生向上的射流。隨氣泡進一步發(fā)展壁面作用逐漸積累，使高壓區(qū)緩慢向右側(cè)移動。盡管氣泡被射流穿透時射流仍豎直向上，但可預(yù)見穿透的氣泡運動中右側(cè)小部分氣泡被撕裂后，左側(cè)氣泡將在高壓區(qū)作用下產(chǎn)生朝向結(jié)構(gòu)的射流，造成二次沖擊損傷。

通過改變浮力參數(shù)分析其對氣泡、自由面運動影響規(guī)律。浮力參數(shù)取δ2=0.15,0.2,0.25,0.3。典型時刻氣泡、自由面形態(tài)見圖9。由圖9可知，不同浮力參數(shù)下自由面形態(tài)區(qū)別較小。氣泡進入坍塌階段后隨浮力參數(shù)增加正下方射流發(fā)展的越早，氣泡被射流穿透時刻發(fā)展越充分。浮力參數(shù)δ2大于0.25時氣泡在浮力主導(dǎo)下運動，自由面、壁面對氣泡影響均不顯著。

圖7　不同時刻氣泡與自由面形態(tài)Fig.7 The bubble and free surface dynamics

圖8　有浮力典型時刻流場壓力云圖 Fig.8 Contour of pressure and the vector of velocity of the fluid field at different time considering buoyancy

不同浮力參數(shù)下氣泡中心在垂向的位移曲線見圖10。由圖10可知，隨浮力參數(shù)增大氣泡向上的位移逐漸增加。由于氣泡產(chǎn)生初期體積較小，受浮力作用較小，因此該階段并未由于浮力作用產(chǎn)生顯著上浮。氣泡膨脹到一定尺寸時浮力作用逐漸顯著，氣泡開始加速上浮。而氣泡上浮靠近自由面時受到自由面向下的排斥作用，在氣泡運動后期浮力參數(shù)越大上浮速度越小。由此知，考慮浮力作用后氣泡在多種載荷聯(lián)合作用下運動，隨同一參數(shù)變化規(guī)律也更復(fù)雜。

圖9　浮力δ2=0.15,0.2,0.25,0.3時氣泡與自由面形態(tài)對比Fig.9 Comparison between bubble shapes at typical moment withdifferent initial distance

圖10不同浮力參數(shù)下氣泡垂向位移曲線Fig.10 Comparison between the bubble volume with different attack angle

3結(jié)論

基于不可壓縮勢流理論建立氣泡與復(fù)雜邊界的耦合動力學(xué)模型，針對流固交界處邊界條件特殊性，通過對交界節(jié)點拆分去除過約束問題。基于該數(shù)值模型計算、分析艦船圓形舭部附近水下爆炸氣泡與自由面的耦合特性，結(jié)論如下：

(1)在舷側(cè)與自由面聯(lián)合作用下, 氣泡坍塌階段流場高壓區(qū)移動是導(dǎo)致射流從豎直向下向垂直于壁面偏轉(zhuǎn)的主要原因。

(2)在研究范圍內(nèi)，隨攻角增大自由面影響逐漸減弱，而壁面影響逐漸增強；隨爆距增大自由面對氣泡影響的衰減慢于壁面。

(3)浮力參數(shù)δ2>0.25后氣泡運動由浮力主導(dǎo)；氣泡被射流穿透后殘余部分會在壁面吸引作用下產(chǎn)生新射流，對結(jié)構(gòu)產(chǎn)生二次損傷。

參 考 文 獻

[1] Blake J R, Gibson D C. Growth and collapse of a vapour cavity near a free surface[J]. Journal of Fluid Mechanics, 1981,111: 123-140.

[2] Chahine G L, Hsiao C T. 3DYNAFS? a three-dimensional free surface and bubble dynamics code bubble dynamics Version 4.1 [R].Dynaflow, INC, User Manual,2006.

[3] Pearson A, Cox E,Blake J R,et al. Bubble interactions near a free surface[J]. Engineering Analysis With Boundary Elements, 2004,28(4): 295-313.

[4] Wang Q X, Yeo K S,Khoo B C,et al. Strong interaction between a buoyancy bubble and a free surface[J]. Theoretical and Computational Fluid Dynamics, 1996,8(1): 73-88.

[5] Blake J R, Gibson D C. Cavitation bubbles near boundaries[J]. Annual Review of Fluid Mechanics, 1987, 19: 99-123.

[6] Wang Q X. The evolution of a gas bubble near an inclined wall[J]. Theoretical and Computational Fluid Dynamics, 1998, 12(1):29-51.

[7] Klaseboer E, Hung K C, Wang C, et al. Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient/rigid structure[J]. Journal of Fluid Mechanics, 2005,53(7): 387-413.

[8] Klaseboer E, Khoo B C, Hung K C. Dynamics of an oscillating bubble near a floating structure[J]. Journal of Fluids and Structures, 2005,10(2): 1-10.

[9] 初文華,張阿漫,王詩平. 壁面與自由液面聯(lián)合作用下氣泡動態(tài)特性實驗研究[J]. 振動與沖擊,2013,32(13):112-117.

CHU Wen-hua, ZHANG A-man,WANG Shi-ping. Experimental study on bubble pulse features under combined action of wall and free surface[J]. Journal of Vibration and Shock,2013,32(13):112-117.

[10] 劉云龍,汪玉,張阿漫. 有傾角的豎直壁面附近氣泡與自由面相互作用研究[J]. 物理學(xué)報,2013,62(21):267-276.

LIU Yun-long, WANG Yu, ZHANG A-man. Interaction between bubble and free surface nearvertical wall with inclination[J]. Acta Physica Sinica,2013,62(21):267-276.

[11] Wang Q X. Unstructured MEL modeling of nonlinear unsteady ship waves[J]. Journal of Computational Physics,2005,210(1):368-385.

[12] Ni B Y, Zhang A M, Wu G X. Numerical simulation of motion and deformation of ring bubble along body surface[J]. Applied Mathematics and Mechanics,2013,34(12):1495-1512.

[13] Cole R H. Underwater explosion[M]. Princeton USA: Princeton University Press,1948.

[14] Zhang A M, et al. Numerical simulation of column charge underwater explosion based on SPH and BEM combination[J]. Computers & Fluids, 2013,71(30): 169-178.

[15] Wang C, Khoo B C, Yeo K S. Elastic mesh technique for 3D BIM simulation with an application to underwater explosion bubbles[J]. Computers & Fluids, 2003,32(9): 1195-1212.

[16] Turangan C K, Ong G P,Klaseboer E,et al. Experimental and numerical study of transient bubble-elastic membrane interaction[J]. Journal of Applied Physics,2006,100(5):054910.

[17] Best J P, Kucera A. A numerical investigation of non-spherical rebounding bubbles[J]. Journal of Fluid Mechanics, 1992,245: 137-191.

Interaction between the bubbles and the free surface near the round bilge

TIANZhao-li1,2,LIUYun-long1,GAOLi-feng3

(1. Harbin Engineering University, Harbin 150001, China;2. Wuhan 2nd Ship Research and Design Institute, Wuhan 430205, China;3. Shanghai Waigaoqiao Shipbuilding & Offshore Engineering Design Co., Ltd., Shanghai 200137, China)

Abstract:The nonlinear interaction between the underwater explosion bubbles and the free surface near the round bilge was studied by using the bubble dynamics model based on the Boundary Element Method. Firstly, a numerical underwater explosion bubble dynamics model with the ideal and incompressible assumption was developed. Then, for the over-constrained problem caused by the intersection between the rigid wall and the free surface, the traditional numerical model was improved by using the double-node method to split each intersection node up to 2 parts while considering the continuity conditions. With the present model, the studies on the interaction between the bubbles and complex boundaries were extended. The influences of the attack angle, the initial distance of the charge, and the buoyancy of the dynamics of the bubbles and the free surface were analyzed, respectively. This work can serve as a reference to the anti-shock research and design of the warships.

Key words:underwater explosion; bubble dynamics; round bilge; free surface; Boundary Element Method

中圖分類號:O352

文獻標志碼：A

DOI:10.13465/j.cnki.jvs.2016.06.023

通信作者劉云龍 男，副教授，1988年5月生

收稿日期：2015-01-13修改稿收到日期：2015-03-31

第一作者 田昭麗 女，工程師，1987年12月生

E-mail：yunlong_liu@hrbeu.edu.cn