權(quán)曉波，孔德才，李 巖

(北京宇航系統(tǒng)工程研究所，100076北京，cxz0813@sohu.com)

水下航行體具備隱蔽、機(jī)動、對敵攻擊具有突然性等優(yōu)勢，在未來戰(zhàn)爭中將發(fā)揮重要作用［1］.水下航行體的出水過程是一個穿越不同介質(zhì)界面，環(huán)境介質(zhì)密度和運(yùn)動速度等參數(shù)劇烈變化的非定常、非線性氣水兩相運(yùn)動過程，將引起流體動力或運(yùn)動載荷的大幅變化，影響水下航行體出水姿態(tài)、運(yùn)動穩(wěn)定性以及在空中正常飛行，甚至導(dǎo)致發(fā)射失敗，是水下航行體水下發(fā)射的關(guān)鍵技術(shù)［2］.文獻(xiàn)［3］指出，水下航行體水下發(fā)射動力學(xué)與空中發(fā)射飛行器動力學(xué)由于側(cè)向載荷分布極為不同而有很大區(qū)別.此外，水下航行體出水時的運(yùn)動規(guī)律和載荷變化規(guī)律還要受到真實(shí)海況下水面波浪情況的影響，不但與波浪傳播方向、浪級有關(guān)，還與其出水相位有關(guān)［2］.劉曜和馬震宇［4］分析了波浪力對水下航行體出水參數(shù)的影響，文獻(xiàn)［5］研究表明，水下航行體出水俯仰角在順/逆浪時的變化范圍在±3°;田兵等［6］的仿真結(jié)果表明，水下航行體在不同海況下及不同的波浪方向角和初相位角下所受到的波浪力和波浪力矩是不同的，即波浪對水下航行體的影響具有隨機(jī)性.

本文首先基于波浪理論得到了二階Stokes波的速度表達(dá)式，并通過給定入射邊界速度分布形式的方法實(shí)現(xiàn)波浪的數(shù)值模擬;然后利用動網(wǎng)格技術(shù)、UDF技術(shù)等，以五級海情為例，探討了水下航行體模型(鈍頭圓柱體)出水時波浪對其三維流場、力學(xué)特性等的影響.

1 波浪模擬研究

海洋近水面處情況極為復(fù)雜，除了受自由液面影響外，還存在風(fēng)、波浪、潮流、涌等多種自然現(xiàn)象，其中波浪對水下航行體出水過程影響較大，若忽略該因素甚至可能導(dǎo)致水下航行體發(fā)射失敗.基于對水波方程的不同展開形式，波浪理論大致可分為線性波和非線性波2類.波面振幅較小時，線性波理論與實(shí)際波浪情況較接近，但是在波面振幅較大時，線性微幅波理論會帶來較大誤差.在實(shí)際海況中，波浪的波高相對于波長(或相對于水深)一般是有限的，在這種有限振幅波中，波動的自由水面引起的非線性影響必須加以考慮，即自由表面的運(yùn)動條件和動力學(xué)條件是非線性的，因此對于實(shí)際海洋中的波浪進(jìn)行研究需要基于非線性波理論.

1847年Stokes將攝動法用于非線性波的研究，把有限振幅的非線性波視為無限多個線性波的Fourier組合，得到了有限水深條件下的二階解和無限水深下的三階解，式(1)即為Stokes波的二階解形式.

波面方程:

由于Stokes波理論適用于水深相對波長較大的深水情況［7］，用于文中模擬水下航行體發(fā)射出水過程的波浪，具體做法是根據(jù)二階Stokes波理論推導(dǎo)出如式(2)和式(3)所示的波面運(yùn)動速度表達(dá)式.

水平速度:

垂直速度:

利用UDF技術(shù)將其作為入射邊界條件耦合到計(jì)算方程中，從而實(shí)現(xiàn)非定常波浪模擬.式中H為波(浪)高，θ為相位角，L為波長，k為波數(shù)，d為不計(jì)及波面的水深，s為計(jì)及波面的水深，T為周期.

為檢驗(yàn)上述造波方法的可行性和準(zhǔn)確性，在二維條件下進(jìn)行了驗(yàn)證［8～9］.氣、水交界面通過VOF模型進(jìn)行模擬［10］，所模擬的二維波浪參數(shù)為:波高為1 m，波長為9 m，周期為2.4 s，水深為5 m.考慮重力和粘性作用.計(jì)算域如圖1(a)所示，為20 m×10 m，采用結(jié)構(gòu)化網(wǎng)格劃分.計(jì)算域上下邊界為固體壁面，參考源項(xiàng)造波的方法［11］，左邊界給定速度入口，右邊界為自由出流條件［12］.計(jì)算過程中采用的湍流模型為標(biāo)準(zhǔn)k-ε模型［13］.

圖1 t=13.6 s時數(shù)值模擬波形及其與理論波形對比曲線

圖1(a)給出了t=13.6 s時的波形圖.從圖中可以看到計(jì)算域中包含有2個波形，這與9 m波長、20 m計(jì)算域這一計(jì)算條件是一致的.圖1 (b)為根據(jù)式(1)得到的Stokes波二階理論解與波浪數(shù)值模擬結(jié)果的對比曲線，可以看出二者波形較為吻合，在4.3 m和13.5 m處分別形成波峰，波峰之間距離與造波時給定的9 m波長相一致.Stokes波二階理論解是在流體無粘、無旋假設(shè)前提下通過攝動法得到的，而波浪數(shù)值仿真中則考慮了流體的粘性作用，通過N-S方程求解得到，因此，數(shù)值波高較理論波高有一定的偏差，且隨著距入射邊界的距離增大，這種衰減越明顯.此外，從圖1(b)還可看出，波峰較為尖突，波谷較為平坦，并且波峰偏離初始平衡位置的距離要大于波谷偏離平衡位置的距離，這也是符合二階Stokes波理論波形的.無旋運(yùn)動的有限振幅Stokes波的質(zhì)點(diǎn)運(yùn)動是不封閉的，在運(yùn)動軌跡上部處于波峰的質(zhì)點(diǎn)比位于波谷的軌跡下部質(zhì)點(diǎn)運(yùn)動得快一些，如圖2(a)所示，一部分質(zhì)點(diǎn)隨著波浪運(yùn)動被輸運(yùn)前進(jìn)了，即Stokes波不是純粹的振動波，而帶有輸移波的特征.圖2(b)給出的計(jì)算結(jié)果也體現(xiàn)出了這一特點(diǎn).綜合以上分析，利用給定計(jì)算域入口速度分布造波的方法模擬出的波浪具有與理論相符且具有精度較高的特點(diǎn)，可以應(yīng)用于工程研究.

圖2 質(zhì)點(diǎn)運(yùn)動軌跡對比

2 波浪對水下航行體出水過程的影響

在上述工作的基礎(chǔ)上，本文開展了波浪對水下航行體出水過程中的三維流場結(jié)構(gòu)、力學(xué)特性、運(yùn)動軌跡等的影響數(shù)值仿真研究.水下航行體的運(yùn)動過程通過求解流體力學(xué)N-S方程得到，根據(jù)二階Stokes波理論得到速度式，利用UDF技術(shù)將其作為入射邊界條件，模擬五級海情下的波浪，具體參數(shù)為:水深 45 m，波高 2.65 m，周期3.58 s，波長20 m.水下航行體模型(鈍頭圓柱體)為水下垂直發(fā)射，發(fā)射水深為30 m，彈射速度為20 m/s，推力為105N.

計(jì)算域網(wǎng)格劃分如圖3所示，為了更為準(zhǔn)確的模擬波面以及波浪與航行體的相互作用，網(wǎng)格在水面附近進(jìn)行了加密處理.圖3中左邊界的入口涉及水氣兩相，通過UDF設(shè)定波面以下為水相，波面以上為氣相，右邊界為自由出流條件，上下邊界為固體壁面條件，前后邊界為對稱面條件，水下航行體給定固體壁面條件.湍流模型仍采用標(biāo)準(zhǔn)k-ε模型.從相圖(圖4)中可以看出，本文模擬的水下航行體模型在出水時處于波谷附近.為了研究波浪對航行體出水過程的影響，還模擬了水面為靜水面時的航行體出水過程，其他條件與有波浪時相同.圖中x坐標(biāo)為偏航方向，y坐標(biāo)代表軸向，z坐標(biāo)為俯仰方向，以下各圖中坐標(biāo)采用相同的定義.

圖3 計(jì)算域網(wǎng)格劃分

圖4 t=1.13 s時相圖

圖5為有波浪時計(jì)算域內(nèi)不同時刻的流線圖，由圖5可知，在水下航行體距水面較遠(yuǎn)時，航行體周圍流場結(jié)構(gòu)基本成軸對稱形，即波浪引起的水質(zhì)點(diǎn)運(yùn)動主要集中在近水面的區(qū)域內(nèi)，在距水面較遠(yuǎn)處波浪對航行體運(yùn)動的影響基本可以忽略.隨航行體向水面運(yùn)動，受波浪引起旋渦的誘導(dǎo)作用，航行體兩側(cè)的附體旋渦結(jié)構(gòu)差別較大，其左側(cè)旋渦變小且逐漸脫離頭部向下運(yùn)動，而右側(cè)旋渦較大且仍存在于肩部，這2個不對稱的旋渦使得水下航行體受到1個繞z軸的逆時針方向的力矩影響，使其逐漸向左側(cè)偏轉(zhuǎn).此外，波浪作用導(dǎo)致的航行體兩側(cè)旋渦結(jié)構(gòu)的差別還與水下航行體出水相位有關(guān)，不同相位出水時的旋渦結(jié)構(gòu)將會發(fā)生變化，這需要進(jìn)一步的研究.

圖5 不同時刻流線圖

從圖6所示的有無波浪時水下航行體運(yùn)動軌跡曲線對比圖中可看出，無波浪靜水面發(fā)射時的航行體的運(yùn)動軌跡近似一條直線，而有波浪存在時其運(yùn)動軌跡沿x軸負(fù)向發(fā)生了偏轉(zhuǎn)，水下航行體尾部完全出水后的偏移量達(dá)到1 m左右，這與水下航行體出水時位于波谷處有關(guān).此外，由于流動的三維性，水下航行體在偏航面上也出現(xiàn)了輕微的偏移.

圖6 水下航行體運(yùn)動軌跡曲線

圖7為水下航行體質(zhì)心沿坐標(biāo)方向的速度變化曲線，可以看出，波浪作用對水下航行體x向的速度分量影響最大，無波浪時x向速度分量基本為0，而波浪存在時，隨水下航行體接近水面，x向速度分量出現(xiàn)負(fù)值且逐漸增加，這與圖5所示的流場結(jié)構(gòu)變化相符.水下航行體y向速度分量在有無波浪時變化趨勢基本一致，但波浪存在時該值較大.這主要是因?yàn)樗潞叫畜w距自由水面的距離隨波浪起伏運(yùn)動而發(fā)生變化，且近水面的區(qū)域內(nèi)流體具有一定的速度，這2個因素導(dǎo)致y向速度分量在數(shù)值上出現(xiàn)差別.在z方向，有無波浪時的速度變化均可視為三維擾動，數(shù)值不大且偏向正方向或偏向負(fù)方向都有隨機(jī)性.

圖7 3個方向速度變化

3 結(jié)論

1)本文基于波浪理論得到了二階Stokes波的波面方程及其運(yùn)動速度表達(dá)式，利用UDF技術(shù)將其作為入射邊界條件耦合到計(jì)算方程中，實(shí)現(xiàn)了非定常波浪模擬，且與理論值吻合良好，精度較高.

2)數(shù)值模擬了五級海情條件下，水下航行體模型(鈍頭圓柱體)水下垂直發(fā)射時波浪對其三維流場、力學(xué)特性等的影響.

3)距水面較遠(yuǎn)處波浪對水下航行體運(yùn)動的影響基本可以忽略;近水面時，受波浪引起的旋渦的誘導(dǎo)作用，航行體兩側(cè)的附體旋渦結(jié)構(gòu)、所處位置等差別較大，使得水下航行體受到繞z軸的逆時針方向的力矩影響，其向左側(cè)偏轉(zhuǎn)，尾部完全出水后的偏移量達(dá)到1 m左右;由于水下航行體距自由水面的距離隨波浪起伏運(yùn)動而發(fā)生變化，且近水面的區(qū)域內(nèi)流體具有一定的速度，水下航行體y向速度分量比無波浪時稍大.

［1］賀小艷.水下氣體射流數(shù)值研究［D］.北京:北京空氣動力研究所，2001.

［2］殷崇一.潛射導(dǎo)彈發(fā)射與出水載荷研究［D］.西安:西北工業(yè)大學(xué)，2004.

［3］BURGDORFO O.Hydrodynamics of unsteady underwater launched missile with trailing cavities and crossflow drag［C］//Fluid Dynamics，Plasma Dynamics，and Lasers Conference.Honolulu:AIAA，1987:AIAA Paper 87-1319.

［4］劉曜，馬震宇.潛載導(dǎo)彈垂直發(fā)射水彈道和分離彈道研究［J］.船舶工程，2005，27(3):6-9.

［5］劉曜.波浪對運(yùn)載器出水姿態(tài)角的影響［J］.艦船科學(xué)技術(shù)，2005，27(3):32-34.

［6］田兵，徐新琦，肖鵬飛.波浪對潛地彈道導(dǎo)彈出水運(yùn)動的影響［J］.艦船科學(xué)技術(shù)，2009，31(4):81-84.

［7］BURCE J，WEST J，THOMSON A，et al.Statistical mechanics of ocean waves［J］.J Hydronautics，1975，9 (1):25-31.

［8］LI Benxia，YU Xiping.A 2-D Numerical Irregular Wave Tank and Its Verification［J］.Journal of hydrodynamics，2005，17(2):222～227.

［9］BRADFORD S F.Numerical simulation of surf zone dynamics［J］.Journal of Waterway Port Coastal and Ocean Engineering-ASCE，2000，126(1):1-13.

［10］ZONG Zhi，DONG Guohai.A second order volume of fluid(VOF)scheme for numerical simulation of 2-D breaking wave［J］.Journal of Marine Science and Application，2007，6(2):1-5.

［11］BRORSEN M，LARSEN J.Source generation of nonlinear gravity waves with boundary integral eqution method［J］.Coastal Engineering，1987，11:93-113.

［12］OJHYAMA T.Development of a numerical wave tank for analysis of nonlinear and irregular wave field［J］. Fluid Dynamics Research，1991(8):231-251.

［13］QI Peng，HOU Yijun.Numerical wave flume study on wave motion around submerged plates［J］.China Ocean Engineering，2003，17(3):397-406.