張昊春，吉　宇，馬　銳，李　垚，嚴利明，秦　江

(1．哈爾濱工業(yè)大學能源科學與工程學院，黑龍江哈爾濱150001; 2．中國艦船研究設(shè)計中心，湖北武漢430064; 3．哈爾濱工業(yè)大學復合材料研究所，黑龍江哈爾濱150001)

水下航行體熱尾流浮升特性研究

張昊春1，吉宇1，馬銳2，李垚3，嚴利明1，秦江1

(1．哈爾濱工業(yè)大學能源科學與工程學院，黑龍江哈爾濱150001; 2．中國艦船研究設(shè)計中心，湖北武漢430064; 3．哈爾濱工業(yè)大學復合材料研究所，黑龍江哈爾濱150001)

摘要:水下航行體排放的熱尾流的浮升規(guī)律及其水面溫度分布特征是紅外探測的基礎(chǔ)。本文從熱尾流形成過程出發(fā)，建立水下射流的二維計算模型。基于Fluent軟件的湍流數(shù)學模型，對溫度分層和均勻環(huán)境介質(zhì)中的熱尾流的浮升特性及水面溫度分布進行數(shù)值計算，并以VOF方法追蹤了溫度分層環(huán)境中的尾流界面，得到熱尾流的軌跡與溫度的衰減規(guī)律，對比分析溫度分層環(huán)境與均勻介質(zhì)環(huán)境中的異同。

關(guān)鍵詞:水下航行體;熱尾流;浮升規(guī)律;溫度分布;數(shù)值模擬

Buoyant characteristics of thermal wakes discharged by underwater vehicles

ZHANG Hao-chun1，JI Yu1，MA Rui2，LI Yao3，YAN Li-ming1，QIN Jiang1
(1．School of Energy Science and Engineering，Harbin Institute of Technology，Harbin 150001，China; 2．China Ship Development and Design Institute，Wuhan 430064，China; 3．Center for Composite Materials，Harbin Institute of Technology，Harbin 150001，China)

Abstract:The buoyant characteristics and temperature distribution of thermal wakes discharged by underwater vehicles are the basis of the detection with infrared technology．This paper started from the formation of thermal wakes，and established a two－dimensional model about the submerged buoyant jet．Using the turbulent flow model in Fluent software，the numerical simulation in temperature stratification case and homogeneous environmental media case were derived，and obtained the buoyant laws as well as temperature distribution．In addition，VOF method was adopted to trace the interface of jet in temperature stratification case; the trajectory of thermal wakes and the decay laws of its temperature were figured out as well．Then a comparison of in the temperature stratification case and in homogenous media case was carried out．

Key words:underwater vehicle; thermal wakes; buoyantcharacteristics; temperature distribution; numerical simulation

0　引言

水下航行體(如潛艇)由于其極強的隱蔽性與突襲性，成為現(xiàn)代海戰(zhàn)中的重要作戰(zhàn)手段。在發(fā)展?jié)撏Ъ夹g(shù)的同時，對反潛技術(shù)的研究同樣受到各國的高度重視。由于潛艇在水下航行時需要冷卻動力裝置，尤其是核潛艇，因此會向海洋中排放大量的溫熱海水而形成熱尾流，熱尾流浮升至海面會使海面水溫發(fā)生異常。國外有專家估測，1臺反應堆功率為190 MW的核潛艇每秒鐘釋放入海洋的熱能多

達1.89×108J，在速度為5 kn時，它排出的熱能使其身后的水溫升高0.2 K［1］。隨著遙感技術(shù)的發(fā)展，利用裝載在飛機或衛(wèi)星上的紅外探測裝置和高速攝像機等設(shè)備儀器可以發(fā)現(xiàn)水面上的熱尾流，從而進行判斷分析是否有潛艇。

為了解熱尾流的形成機理與浮升規(guī)律，從而為紅外探潛提供理論依據(jù)。國內(nèi)外已有許多學者針對此問題從不同的角度來進行研究。Garber等［2］在1945年對艦船熱尾流進行了開拓性的研究，并對其溫度分布特征進行了實船探測; Garrett等［3］研究了潛艇排放的冷卻水在相對穩(wěn)定環(huán)境運動條件下和分層流中的流跡和衰減; Peltzer等［4］利用機載紅外成像儀對Hayes號雙體船的尾流進行遙感測量以研究紊流尾流熱表面效應及少量有機物質(zhì)對熱表面效應的影響; Schwartz和Priest［5］從概率角度考慮了船尾熱尾流與粗糙海面的紅外對比度，對艦船尾流進行了紅外遙感研究;美國海軍研究實驗室(NRL)開發(fā)了一個艦船紊流尾流計算的通用程序TWAKE［6］，采用定常、三維、時間平均拋物線型N－S方程及k－ε封閉模式，利用有限元方法進行計算; Yang和Gu等［7－12］對水下棲息潛艇形成的熱尾流的浮升與衰減規(guī)律進行了大量的理論與實驗研究。此外，在射流研究方面，Larsen等［13］研究了橫流中的浮力射流特性并比較了積分模型與k－ε模型; K.Harby等［14］通過實驗研究了液態(tài)環(huán)境中淹沒的橫向氣體射流的流動特性; Masaki和Seiji［15］利用分布式注入模型研究了分層的熱水箱中的橫向射流的特性。

本文將水下航行體排放冷卻水與螺旋槳攪動行為的作用結(jié)果視為“熱射流”，分別建立溫度分層環(huán)境中與均勻介質(zhì)環(huán)境中的射流模型，利用有限體積法進行數(shù)值計算分析，得到了熱尾流的溫度分布特性以及浮升擴散過程中的衰減規(guī)律。同時，利用VOF (Volume of Fluid)方法追蹤“熱射流”的界面。

1　數(shù)理模型

1.1計算模型及參數(shù)

本文所采用的計算模型基于二維結(jié)構(gòu)，計算區(qū)域為50 m×10 m;水下航行物體模型位于計算區(qū)域的下部，距離底面2.5 m，距離左側(cè)速度入口18.5 m，整個計算區(qū)域中充滿海水，水面和水底均采用剛性壁面邊界條件。由于在實際情況下，所排放的冷卻水全被螺旋槳卷吸進入尾流的中心區(qū)域，因此熱尾流形成的初始階段因冷卻水而溫度升高的區(qū)域在尾流的中央，其尺寸一般由螺旋槳的直徑Dr決定，通常為潛艇直徑D的一半。因此本文的研究重心在于所形成的熱尾流在海水中的浮升擴散，而將螺旋槳與冷卻水的作用結(jié)果直接視為圓形“熱射流”，如圖1所示。

圖1　計算區(qū)域示意圖Fig.1　The schematic view of calculation area and its geometry

為了考察熱射流的浮升擴散特性，在本文的數(shù)值模擬中，分別建立了溫度分層和均勻介質(zhì)海水模型。考慮到具體的海洋情況，在溫度分層計算模型中，以交界面處為界，上層介質(zhì)為T1= 293.15 K的海水，下層介質(zhì)為T2=290 K的海水，相關(guān)物性均以實際數(shù)據(jù)為準。由于海洋的分層現(xiàn)象，在潛艇航行的水層，其相對于海水的速度應大于無潛艇航行的水層，因此設(shè)定上層水流速度為0.8 m/s，下層速度為1 m/s;出口為自由出流。相應地，在均勻海水介質(zhì)模型中，整個計算區(qū)域中的海水均為T = 293.15 K。

由于海水與水下物體壁面的摩擦、冷卻水的排放以及螺旋槳攪動等綜合因素，因此在本文的計算中，設(shè)定射流出口溫度為Tjet= 310 K，速度v=0.4 m/s。

本數(shù)值模型通過網(wǎng)格獨立性檢驗后采用混合網(wǎng)格進行計算，所選流動模型為非穩(wěn)態(tài)Realizable k－ε模型，并采用VOF方法跟蹤了溫度分層計算模型中的“熱射流”界面。

1.2流動傳熱控制方程

連續(xù)性方程:

動量方程:

能量方程:

其中

狀態(tài)方程:

第p相的體積分數(shù)方程［16］:

式中:ρ為流體密度; ui(i = 1，2)分別為x，y方向的速度分量; p為流體靜壓;μ為流體的動力粘度; cp為定壓比熱; T為溫度; Pr為流體的普朗特數(shù); gi為重力加速度分量;αp為第p相的體積分數(shù)。

湍流動能方程為:

湍流耗散率方程為:

2　計算結(jié)果與分析

水下航行體排放冷卻水時，由于冷卻水溫度較高，排出后則處于跟周圍的海水進行熱量、質(zhì)量和動量交換的過程?？紤]到溫度分層情況下，所得到的溫度場、速度場及相分布云圖如圖2～圖4所示。

圖2　溫度分層環(huán)境中的溫度場分布Fig.2　The temperature field in temperature stratification case

圖3　溫度分層環(huán)境中的速度場分布Fig.3　The velocity field in temperature stratification case

圖4　溫度分層環(huán)境中的相分布Fig.4　The phase distribution in temperature stratification case

由上述計算結(jié)果可以看出，在溫度分層環(huán)境中，冷卻水等形成的“熱射流”對整個計算區(qū)域的溫度場分布有較大的影響，在豎直方向上，原有的兩層溫度水域中出現(xiàn)了連續(xù)的溫度變化，而水下航行體的運動對速度場也產(chǎn)生了很大的影響，潛艇附近出現(xiàn)了明顯的繞流現(xiàn)象，進而造成上層水域的速度有一定范圍的增大。為了獲得熱尾流在擴散過程中的溫度衰減變化特性，對比艇后不同位置處海水溫度隨距離水面深度的變化曲線如圖5所示。

從圖5可看出，所排放的溫水對于水下航行體下部水域的影響較小，僅在熱尾流附近引起了較大的溫度變化，但在上部水域中形成了明顯的“冷特征”。此現(xiàn)象是由于熱射流在浮升擴散過程中，與外部水域存在較強烈的動量交換和質(zhì)量交換，導致下層溫度較低的水被翻至近水面處。同

時也可以看出，隨著艇后距離的加大，溫度分布曲線的峰值有所減小，即存在溫度衰減現(xiàn)象且在距離艇尾較近的區(qū)域，尾流溫度衰減較快，隨著與艇體距離的增大，熱尾流的溫度衰減越來越緩慢。

圖5　艇后不同位置處海水溫度深度變化情況Fig.5　Variation of temperature with depth in different distances away from submarine

同時，對于均勻介質(zhì)海水的計算模型，其溫度場分布與速度場分布如圖6和圖7所示。

圖6　溫度均勻海水中的溫度場分布Fig.6　The temperature field in uniform temperature case

圖7　溫度均勻海水中的速度場分布Fig.7　The velocity field in uniform temperature case

從圖中可看出，在均勻介質(zhì)環(huán)境中，艇體對于海水域的擾動也很大，速度場基本與溫度分層環(huán)境無太大差別，但是溫度場的分布卻存在很大的差異，除形成的“熱射流”附近有明顯的溫度變化之外，其余的區(qū)域基本上沒有受到影響。造成此現(xiàn)象的原因可能在于，均勻介質(zhì)環(huán)境中的密度差異較小以至于尾流不易于浮升至水面，從而很難在水面上形成較明顯的溫度變化。

將溫度分層與均勻介質(zhì)海水模型中所得到的水面溫度對比可以發(fā)現(xiàn)(見圖8)，前者呈現(xiàn)出明顯的“冷尾流”特征，且在所計算的區(qū)域內(nèi)隨著尾流的擴展，水面溫度越低，最大溫差可以達到0.35 K;而后者幾乎看不到水面上有任何的溫度變化，即尾流引起的水面冷熱特征很微弱。

圖8　溫度分層與均勻溫度海水模型中水面溫度變化Fig.8　The temperature on the water surface in the temperature stratification case and uniform temperature case

3　結(jié)語

本文針對于水下航行體在運動過程中所排放的冷卻水和螺旋槳攪動等行為對周圍的水域所造成的擾動，依據(jù)流體流動與傳熱特性，建立了二維的熱尾流浮升擴散模型，分析了熱尾流形成、浮升和衰減過程的特性，利用Fluent軟件進行了數(shù)值計算，得到主要結(jié)論如下:

1)相同運動條件下，溫度分層環(huán)境中的尾流更容易擴散浮升至水面，且在水面上形成明顯的“冷尾流”特征，最大溫差可達0.35 K;而在均勻介質(zhì)海水模型中，由于密度與溫度幾乎無差異，因此水下的“熱射流”難以浮升至水面上，水面上沒有顯著的溫差。

2)溫度分層環(huán)境中的熱尾流在浮升擴散的過程中，在豎直方向上形成了一定的溫度變化，且在越接近水下航行體的深度，溫度變化越劇烈，而在離水下物體較遠的深度處，溫度變化較緩慢。

3)隨著艇后距離的加大，峰值溫度在不斷的減小，即存在溫度衰減特性，且越靠近潛艇，尾流溫度衰減越快。這在一定程度上能給紅外探測提供相應的理論支持。

參考文獻:

［1］李蓉芳．探測潛艇的非聲吶方法［J］．科學，1998(7) :31．

LI Rong-fang．Non sonar method to detect submarine［J］．Science，1998(7) :31．

［2］GARBER D H，VRICK R J，CRYDEN J．Thermal wake detection［R］．Navy Radio andSoundLaboratory: AD－493809，1945．

［3］GARRETT W D，SMITH P M．Physical and chemical factors affecting the thermal IR imagery of ship wake［R］．Navy Radio and Sound Laboratory:AD－A144034，1984．

［4］PELTZER R D，GARRETT W D，SMITH P M．A remote sensing study of a surface ship wake［J］．International Journal of Remote Sensing，1978，8(5) :689－702．

［5］SCHWARTZ I B，PRIEST R G．Reflection driven ship wake contrasts in the infrared［M］，NRL－9144，1988．

［6］SWEAN T F．A User's/Programmers manual for TWAKE ［R］．NRL－MR－6201，1988．

［7］YANG Li，HUA Shun-fang，DU Xian-zhi．Investigation of buoyant trajectory of thermal wake behind a slender streamlined heated body of revolution in stratified fluid［J］．Journal of Hydrodynamics，Ser．B．1992，4(2) :59－63．

［8］ZHANG Jian，CHEN Xuan，YANG Li，et al．Study of temperature characteristic of cooling water discharged by underwater vehicle［J］．Journal of Ship Mechanics，2009，13(4) :533－539．

［9］WANG Jian-gan，GUO Yan，GU Jian-nong．The theoretical and trial study of thermal wake in the infrared detection of submarine［J］．Laser＆Infrared，2002，32(3) :159－162．

［10］顧建農(nóng)，鄭學嶺．艦船熱尾流研究述評［J］．魚雷技術(shù)，1999(2) :5－8．GU Jian-nong，ZHENG Xue-ling．Review on the thermal wake of ship［J］．Torpedo Technology，1999(2) :5－8．

［11］楊立．熱尾流浮升規(guī)律的理論與實驗研究［D］．武漢:海軍工程學院，1989． YANG Li．Theoretical and experimental studies on the buoyancy law of thermal wake［D］．Wuhan: Naval University of Engineering，1989．

［12］顧建農(nóng)，鄭學嶺，劉巨斌，等．熱分層環(huán)境下船模尾流溫度分布與擴展規(guī)律實驗研究［J］．大連理工大學學報，2001(11) :93－98．GU Jian-nong，ZHENG Xue-ling，LIU Ju-bin，et al．Experimental research on temperature distribution and spread law of a ship model wakes in thermal stratified flow ［J］．Journal of Dalian University of Technology，20011 (11) :93－98．

［13］LARSEN T，PETERSEN O，CHEN H B．Numerical experiment on computation on turbulent buoyant jets in flowing ambient［C］/ /Proceedings of the 8－th International Conference on Computational Methods in Water Resources．Berlin: Springer－Verlag，1990:339－334．

［14］HARBY K，CHIVA S，MUNOZ J．An experimental investigation on the characteristics of submerged horizontal gas jets in liquid ambient［J］．Experimental Thermal and Fluid Science，2014，53:26－29．

［15］TOYOSHIMA M，OKAWA S．An effect of horizontal buoyant jet on the temperature distribution inside a hot water storage tank［J］．International Journal of Heat and Fluid Flow，2013，44:403－413．

［16］陳伯義，雷渡民，吳猛猛，等．水下運動體尾流的浮升規(guī)律及其水面冷熱特征的數(shù)值模擬［J］．紅外與激光工程，2012，41(5) :1140－1146．CHEN Bo-yi，LEI Du-min，WU Meng-meng，et al．Numerical simulation on the buoyant laws and the coldthermal surface features of the wake behind a moving underwater body［J］．Infrared and Laser Engineering，2012，41(5) :1140－1146．

作者簡介:張昊春(1977－)，男，副教授，從事高速飛行器熱控技術(shù)、紅外隱身、工程熱力學和傳熱傳質(zhì)學研究。

基金項目:國家自然科學基金資助項目(91216123) ;國家自然科學基金創(chuàng)新研究群體資助項目(51121004)

收稿日期:2014－11－18;修回日期: 2014－12－08

文章編號:1672－7649(2015) 07－0024－05doi:10.3404/j.issn.1672－7649.2015.07.006

中圖分類號:U661．1

文獻標識碼:A