潛艇熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性

張旭升，郭 亮，胡日查，常雯娟

?

潛艇熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性

張旭升，郭 亮，胡日查，常雯娟

（中國科學(xué)院長春光學(xué)精密機(jī)械與物理研究所，吉林 長春 130033）

選取1/72龍鯊Ⅱ號(hào)核潛艇為研究對象，采用有限體積法建立外流域的三維計(jì)算模型，運(yùn)用動(dòng)參考系實(shí)現(xiàn)尾部螺旋槳高速旋轉(zhuǎn)的仿真。基于該模型，評(píng)估下潛深度、航行速度、熱尾流噴速、熱尾流溫度等因素對熱尾流浮升擴(kuò)散規(guī)律及海表溫度特性的影響。研究表明：1）隨著下潛深度和航行速度的增大，海表溫度特性的后向延遲距離增大、海表溫差減小；2）隨著熱尾流噴速和熱尾流溫度的增大，海表溫度特性的后向延遲距離減小、海表溫差增大；3）增大高溫?zé)嵛擦鳠崛萘炕驕p小低溫海水熱容量均能促使海表溫度特性趨于明顯。本文結(jié)論可為某型號(hào)潛艇的紅外探測與反演定位提供針對性的參考與借鑒。

熱尾流；浮升擴(kuò)散規(guī)律；海表溫度特性；有限體積法；紅外探測

0 引言

在現(xiàn)代化海戰(zhàn)中，潛艇承擔(dān)著隱蔽、偵查、突襲敵方海上作戰(zhàn)單元的重要角色，是具有“航母獵手”稱號(hào)的關(guān)鍵戰(zhàn)略性武器，因此世界各國已著重開展?jié)撏c反潛技術(shù)的研究。依據(jù)探測原理不同，反潛方式可分為聲納探測[1-2]、可見光探測[3-4]、紅外探測[5-7]、電磁探測[4,7]等。隨著機(jī)械制造精度與消磁降噪技術(shù)的進(jìn)步，以聲納探測為主、非聲探測為輔的聯(lián)合式探測將成為反潛技術(shù)發(fā)展的必然趨勢[5]。在非聲探測領(lǐng)域，紅外探測具有全天候工作、抗干擾能力強(qiáng)、快速實(shí)時(shí)成像、易于配合使用等特點(diǎn)[6-7]，是未來反潛技術(shù)的重要研究方向。

張健等[8-9]搭建了拖行水下航行器的熱尾流溫度測量系統(tǒng)，獲得了溫度分層水池中熱尾流的溫度分布規(guī)律及主次影響因素。吳猛猛等[10-11]基于上述熱尾流溫度測試系統(tǒng)，通過實(shí)驗(yàn)仿真相結(jié)合的方式研究了橢圓體熱尾流浮升-摻混-衰減全過程的浮升擴(kuò)散規(guī)律和水面溫度特征。陳伯義等[12]探討了海水的垂直溫度分布和有無氣泡等因素對潛艇熱尾流的浮升擴(kuò)散軌跡和水面冷熱特征的影響。Chen等[13-14]開展了水下航行體冷/熱尾流浮升軌跡及表面溫差信號(hào)的研究，實(shí)現(xiàn)了仿真計(jì)算結(jié)果與紅外相機(jī)實(shí)測數(shù)據(jù)的基本一致。張昊春等[15]建立了簡化水下航行體熱尾流的二維計(jì)算模型，分析了熱尾流的浮升軌跡和溫度衰減等規(guī)律。戴天奇等[16]基于動(dòng)網(wǎng)格技術(shù)建立了潛艇熱尾流的二維計(jì)算模型。該方法能夠更加全面、準(zhǔn)確地捕獲熱尾流流動(dòng)-傳熱過程的瞬變特征。

不難發(fā)現(xiàn)，上述各文獻(xiàn)在熱尾流浮升擴(kuò)散規(guī)律及海表溫度特性的研究過程中存在如下簡化：1）受計(jì)算資源與計(jì)算效率的限制，采用二維仿真計(jì)算近似代替三維；2）采用橢圓體等簡化幾何模型代替流線型潛艇，忽略艇體輪廓對外流場的影響；3）未考慮潛艇尾部螺旋槳的高速旋轉(zhuǎn)，造成外流場湍流擾動(dòng)項(xiàng)不完備。鑒于上述原因，本文選取1/72龍鯊Ⅱ號(hào)核潛艇為研究對象，采用有限體積法建立外流域的三維計(jì)算模型，運(yùn)用動(dòng)參考系實(shí)現(xiàn)尾部螺旋槳高速旋轉(zhuǎn)的仿真?；谠撃Ｐ?，評(píng)估下潛深度、航行速度、熱尾流噴速、熱尾流溫度等因素對熱尾流浮升擴(kuò)散規(guī)律及海表溫度特性的影響，為某型號(hào)潛艇的紅外探測與反演定位提供針對性的參考與借鑒。

1 物理模型及邊界條件

1.1 物理模型

綜合對比國內(nèi)外幾種典型核潛艇縮比模型的幾何結(jié)構(gòu)特征和動(dòng)力性能參數(shù)，選取自航式龍鯊Ⅱ號(hào)核潛艇縮比模型作為仿真分析和實(shí)驗(yàn)研究的對象。該模型主要由艇體、指揮臺(tái)圍殼及側(cè)翼、尾翼、螺旋槳及環(huán)狀導(dǎo)管、熱尾流出口管等結(jié)構(gòu)組成，如圖1所示。長度＝866mm，最大直徑max＝110mm，熱尾流出口管內(nèi)徑＝7mm，與原模型的尺寸比例為1:72。通過控制側(cè)翼、尾翼姿態(tài)和螺旋槳轉(zhuǎn)速與方向，實(shí)現(xiàn)潛艇模型航行方向與速度的調(diào)節(jié)。

依據(jù)潛艇縮比模型的幾何特征尺寸和具體實(shí)驗(yàn)條件，確定外流場計(jì)算域的幾何尺寸。其中，潛艇前端流域長度2，潛艇后端流域長度9，計(jì)算域總長度12；潛艇兩側(cè)流域?qū)挾?.75m，計(jì)算域總寬度1.5m；基準(zhǔn)下潛深度下，潛艇中軸線距自由表面高度0.5m，計(jì)算域總高度1.0m。即：外流場計(jì)算域的長×寬×高分別為10.392m×1.5m×1.0m，如圖2所示。

基于外流場計(jì)算域的對稱性，以1/2計(jì)算域?yàn)檠芯繉ο筮M(jìn)行計(jì)算流體動(dòng)力學(xué)分析是合理、準(zhǔn)確且高效的。同時(shí)，將螺旋槳周圍的圓柱形流域單獨(dú)分離出來，由于其結(jié)構(gòu)相對復(fù)雜且與整體流域尺寸相差較大，故采用較為強(qiáng)健的Meshing模塊劃分螺旋槳計(jì)算域的三維非結(jié)構(gòu)網(wǎng)格，如圖3所示，再結(jié)合動(dòng)參考系即可實(shí)現(xiàn)尾部高速旋轉(zhuǎn)螺旋槳強(qiáng)湍流擾動(dòng)的仿真。

1.2 邊界條件

邊界條件是獲得計(jì)算域流場變量唯一解的必要條件。1/2計(jì)算域外表面的邊界條件如圖4所示，包括：速度入口、出口流動(dòng)等進(jìn)出口條件，靜止壁面、平移壁面、旋轉(zhuǎn)壁面、對稱面、交界面等壁面條件，潛艇計(jì)算域、螺旋槳計(jì)算域等單元分區(qū)。在動(dòng)參考系中設(shè)置旋轉(zhuǎn)軸、旋轉(zhuǎn)方向和絕對轉(zhuǎn)速等，以仿真潛艇尾部高速旋轉(zhuǎn)螺旋槳引起的強(qiáng)湍流擾動(dòng)。文中采用SIMPLE算法和RNG-湍流模型求解粘性流體的三維Navier-Stokes方程。

圖1 龍鯊Ⅱ號(hào)核潛艇縮比模型二維視圖

圖2 潛艇縮比模型的外流場計(jì)算域

圖3 螺旋槳計(jì)算域表面的網(wǎng)格劃分

2 分析與討論

結(jié)合龍鯊Ⅱ號(hào)核潛艇縮比模型的動(dòng)力性能參數(shù)，規(guī)劃了分析熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性的計(jì)算工況，如表1所示。其中，涉及參數(shù)包括：下潛深度、航行速度s、熱尾流噴速h和熱尾流溫度h，標(biāo)注星號(hào)（*）的參數(shù)為基準(zhǔn)參數(shù)，所有基準(zhǔn)參數(shù)共同組成基準(zhǔn)工況。

2.1 下潛深度

在基準(zhǔn)參數(shù)條件下，分析下潛深度對熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性的影響。圖5中給出1/2計(jì)算域?qū)ΨQ截面的熱尾流溫度分布，顯而易見，隨著下潛深度的增大，潛艇上側(cè)低溫海水的熱容量增大，造成熱尾流浮升擴(kuò)散的中心高溫區(qū)域縮小，進(jìn)而嚴(yán)重削弱海洋表面熱尾流的紅外輻射特性。圖6中給出海洋表面的熱尾流溫度分布，在基準(zhǔn)工況下，熱尾流呈現(xiàn)海表溫度特性的后向延遲距離為6.52m、海表溫差為8mK；而當(dāng)＝0.3m時(shí)，海表溫度特性的后向延遲距離僅為4.38m、比基準(zhǔn)工況超前32.8%，海表溫差達(dá)到20mK。即：隨著下潛深度的增大，海表溫度特性的后向延遲距離增大、海表溫差減小，不利于實(shí)現(xiàn)敵方潛艇的紅外探測與反演定位。

2.2 航行速度

在基準(zhǔn)參數(shù)條件下，分析航行速度s對熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性的影響。分析圖7和圖8得出，隨著航行速度s和螺旋槳轉(zhuǎn)速r的增大，熱尾流浮升擴(kuò)散的中心高溫區(qū)域逐漸縮小，主要由于尾部螺旋槳的高速旋轉(zhuǎn)促進(jìn)了高溫?zé)嵛擦髋c臨近低溫海水的摻混和區(qū)域溫度均勻化。因此，需要考慮螺旋槳高速旋轉(zhuǎn)的強(qiáng)湍流擾動(dòng)項(xiàng)對熱尾流特性的影響。同時(shí)，海表溫度特性的后向延遲距離略有增大、海表溫差減小。當(dāng)s＝2.0m/s時(shí)，海表溫度特性的后向延遲距離增大為6.60m、比基準(zhǔn)工況延后1.2%，海表溫差減小為6mK。

圖4 潛艇縮比模型1/2計(jì)算域的邊界條件

表1 仿真分析的計(jì)算工況

圖5 下潛深度對熱尾流浮升擴(kuò)散規(guī)律的影響

圖6 下潛深度對熱尾流海表溫度特性的影響

圖7 航行速度對熱尾流浮升擴(kuò)散規(guī)律的影響

圖8 航行速度對熱尾流海表溫度特性的影響

2.3 熱尾流噴速

在基準(zhǔn)參數(shù)條件下，分析熱尾流噴速h對熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性的影響。分析圖9和圖10得出，隨著熱尾流噴速h的增大，熱尾流浮升擴(kuò)散的中心高溫區(qū)域急劇增大，海表溫度特性的后向延遲距離減小、海表溫差增大。在h＝1.3m/s時(shí)，熱尾流呈現(xiàn)海表溫度特性的后向延遲距離為5.79m，比基準(zhǔn)工況超前11.2%，海表溫差增大為12mK。

2.4 熱尾流溫度

在基準(zhǔn)參數(shù)條件下，分析熱尾流溫度h對熱尾流的浮升擴(kuò)散規(guī)律及海表溫度特性的影響。分析圖11和圖12所示，隨著熱尾流溫度h的增大，熱尾流浮升擴(kuò)散的中心高溫區(qū)域增大，海表溫度特性的后向延遲距離減小、海表溫差增大。在h＝90℃時(shí)，熱尾流呈現(xiàn)海表溫度特性的后向延遲距離為6.26m、比基準(zhǔn)工況超前4.0%，海表溫差增大為10mK。

圖9 熱尾流噴速對熱尾流浮升擴(kuò)散規(guī)律的影響

圖10 熱尾流噴速對熱尾流海表溫度特性的影響

圖11 熱尾流溫度對熱尾流浮升擴(kuò)散規(guī)律的影響

圖12 熱尾流溫度對熱尾流海表溫度特性的影響

綜上所述，下潛深度和航行速度s對熱尾流特性的影響趨勢相同，均增加了潛艇外流場低溫海水的熱容量，但下潛深度的影響效果更為強(qiáng)烈；熱尾流噴速h和熱尾流溫度h對熱尾流特性的影響趨勢相同、影響效果相當(dāng)，均增加了潛艇排放高溫?zé)嵛擦鞯臒崛萘俊?/p>

3 結(jié)論

本文采用有限體積法建立1/72龍鯊Ⅱ號(hào)核潛艇外流域的三維計(jì)算模型，運(yùn)用動(dòng)參考系實(shí)現(xiàn)尾部螺旋槳高速旋轉(zhuǎn)的仿真?；谠撃Ｐ?，評(píng)估下潛深度、航行速度、熱尾流噴速、熱尾流溫度等因素對熱尾流浮升擴(kuò)散規(guī)律及海表溫度特性的影響，所得結(jié)論如下：

1）下潛深度和航行速度s對熱尾流特性的影響趨勢相同，但下潛深度的影響效果更為強(qiáng)烈。隨著某一參數(shù)的增大，熱尾流浮升擴(kuò)散的中心高溫區(qū)域縮小，海表溫度特性的后向延遲距離增大、海表溫差減小。

2）熱尾流噴速h和熱尾流溫度h對熱尾流特性的影響趨勢相同、影響效果相當(dāng)。隨著某一參數(shù)的增大，熱尾流浮升擴(kuò)散的中心高溫區(qū)域增大，海表溫度特性的后向延遲距離減小、海表溫差增大。

3）增大高溫?zé)嵛擦鳠崛萘炕驕p小低溫海水熱容量均能促使海表溫度特性趨于明顯，更有利于實(shí)現(xiàn)敵方潛艇的紅外探測與反演定位。

[1] 張士成, 楊楨, 楊立, 等. 水下航行體自由表面波浪尾跡紅外特征及探測[J]. 紅外與激光工程, 2012, 41(10): 2615-2620.

ZHANG Shicheng, YANG Zhen, YANG Li, et al. Infrared characterization and detection of free surface wave wake of underwater vehicle[J]., 2012, 41(10): 2615-2620.

[2] 王江安, 韋文濤, 李曉龍, 等. 水下航行體熱尾流特性的技術(shù)融合[J]. 艦船科學(xué)技術(shù), 2011, 33(12): 9-13.

WANG Jiangan, WEI Wentao, LI Xiaolong, et al. Research on characteristics of underwater vehicles’ thermal wake by technology convergence[J]., 2011, 33(12): 9-13.

[3] FU Huiping, WAN Pengcheng. Numerical simulation on ship bubbly wake[J]., 2011, 10: 413-418.

[4] 朱海榮, 朱海, 劉金濤, 等. 水下航行器光學(xué)隱蔽深度測量系統(tǒng)[J]. 光學(xué)精密工程, 2015, 23(10): 2778-2784.

ZHU Hairong, ZHU Hai, LIU Jintao, et al. Measurement system of optical concealment depth of underwater vehicle[J]., 2015, 23(10): 2778-2784.

[5] 王建勛, 孫國倉, 鄧海華, 等. 熱尾流紅外探測作用距離工程估算方法[J]. 艦船科學(xué)技術(shù), 2015, 37(4): 107-111.

WANG Jianxun, SUN Guocang, DENG Haihua, et al. Engineering evaluation method for detection distance of submarine thermal wake by IR[J]., 2015, 37(4): 107-111.

[6] 蔣志忠, 楊日杰, 李德鑫, 等. 機(jī)載紅外探測儀檢查搜潛建模與仿真[J]. 紅外與激光工程, 2011, 40(3): 390-396.

JIANG Zhizhong, YANG Rijie, LI Dexin, et al. Modeling and simulation for checking search submarine by airborne infrared detectors[J]., 2011, 40(3): 390-396.

[7] 王建華, 趙浩淞. 機(jī)載紅外探潛系統(tǒng)綜述[J]. 激光與紅外, 2013, 43(6): 599-603.

WANG Jianhua, ZHAO Haosong. Survey of airborned infrared system for submarine detection[J]., 2013, 43(6): 599-603.

[8] 張健, 楊立, 袁江濤, 等. 水下航行器熱尾流試驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué), 2008, 22(3): 7-13.

ZHANG Jian, YANG Li, YUAN Jiantao, et al. Experimental research of underwater vehicle thermal wakes[J]., 2008, 22(3): 7-13.

[9] 張健, 陳翾, 楊立, 等. 水下航行器排放冷卻水溫度分布特性研究[J]. 船舶力學(xué), 2009, 13(4): 533-539.

ZHANG Jian, CHEN Xuan, YANG Li, et al. Study of temperature characteristic of cooling water discharged by underwater vehicle[J]., 2009, 13(4): 533-539.

[10] 吳猛猛, 陳伯義, 張修峰, 等. 溫度分層海洋中水下航行體引起水面特征參數(shù)變化規(guī)律的研究[J]. 紅外技術(shù), 2010, 32(4): 242-246.

WU Mengmeng, CHEN Boyi, ZHANG Xiufeng, et al. Study on the varied laws of surface feature parameters caused by a going body underwater in the temperature stratification ocean[J]., 2010, 32(4): 242-246.

[11] 吳猛猛, 龔文超, 袁寶吉. 水下航行體熱尾流水面溫度特征的研究[J]. 紅外, 2011, 32(9): 29-34.

WU Mengmeng, GONG Wenchao, YUAN Baoji. Study of surface temperature features of thermal wake caused by underwater vehicle[J]., 2011, 32(9): 29-34.

[12] 陳伯義, 雷渡民, 吳猛猛, 等. 水下運(yùn)動(dòng)體尾流的浮升規(guī)律及其水面冷熱特征的數(shù)值模擬[J]. 紅外與激光工程, 2012, 41(5): 1140-1146.

CHEN Boyi, LEI Dumin, WU Mengmeng, et al. Numerical simulation on the buoyant laws and the cold-thermal surface features of the wake behind a moving underwater body[J]., 2012, 41(5): 1140-1146.

[13] CHEN Shengtao, LIU Huanying, QI Yi. The thermal character of the underwater heat exhausting source[J]., 2010, 19(5): 445-450.

[14] CHEN Shengtao, ZHONG Jingjun, SUN Peng. Numerical simulation and experimental study of the submarine’s cold wake temperature character[J]., 2014, 23(3): 253-258.

[15] 張昊春, 吉宇, 馬銳, 等. 水下航行體熱尾流浮升特性研究[J]. 艦船科學(xué)技術(shù), 2015, 37(7): 24-28.

ZHANG Haochun, JI Yu, MA Rui, et al. Buoyant characteristics of thermal wakes discharged by underwater vehicles[J]., 2015, 37(7): 24-28.

[16] 戴天奇, 姚世衛(wèi), 魏志國. 基于動(dòng)網(wǎng)格技術(shù)的潛艇熱尾流浮升規(guī)律研究[J]. 艦船科學(xué)技術(shù), 2015, 37(5): 86-89.

DAI Tianqi, YAO Shiwei, WEI Zhiguo. Numerical simulation of thermal wake buoyant law based on dynamic meshing technique[J]., 2015, 37(5): 86-89.

Buoyancy Diffusion Law and Ocean-surface Temperature Characteristic ofSubmarine Thermal Wake

ZHANG Xusheng，GUO Liang，HU Richa，CHANG Wenjuan

(,,,130033,)

1/72 LSⅡ nuclear submarine is selected as the research object, three-dimensional numerical model of external flow-field is established by using the finite volume method, and high-speed rotation of tail-propeller is simulated by motion frame. Based on the model, factors which affect the buoyancy diffusion law and ocean-surface temperature characteristic, such as diving depth, sailing velocity, and thermal wake velocity/temperature are evaluated. Research shows that 1) with the diving depth and sailing velocity increasing, backward delay distance increases and temperature difference decreases towards ocean-surface temperature characteristic; 2) with the thermal wake velocity/temperature increasing, backward delay distance decreases and temperature difference increases towards ocean-surface temperature characteristic; 3) increasing the heat capacity of high-temperature thermal wake or decreasing the heat capacity of low-temperate seawater makes the ocean-surface temperature characteristic more obvious. The conclusion could provide the targeted reference for infrared detection and inverse positioning of certain type submarine.

thermal wake，buoyancy diffusion law，ocean-surface temperature characteristic，finite volume method，infrared detection

U661.1，TN215

A

1001-8891(2016)08-0678-05

2016-03-11；

2016-06-05.

張旭升（1988-），男，工學(xué)碩士，助理研究員，主要從事紅外熱輻射傳輸、航天器熱控技術(shù)等方面研究。E-mail：zxs1933@126.com。

中國科學(xué)院長春光學(xué)精密機(jī)械與物理研究所創(chuàng)新基金（Y4CX1SS144）。