馮擁軍，顧鑫鑫，郭月姣，董才道，馮國增

(1.江蘇科技大學(xué) 資產(chǎn)經(jīng)營有限公司，江蘇 鎮(zhèn)江 212000；2.江蘇科技大學(xué) 能源與動(dòng)力學(xué)院，江蘇 鎮(zhèn)江 212000)

0 引言

半潛船主要用于裝運(yùn)超大、超重且不可分割的特殊物體，也可用來運(yùn)輸大型海上石油鉆井平臺(tái)、大型船舶、預(yù)制橋梁構(gòu)件等。半潛船功能多樣，可以實(shí)現(xiàn)上浮下潛，具有先進(jìn)的動(dòng)力定位系統(tǒng)。隨著深海資源的不斷開發(fā)和國際間合作項(xiàng)目的日益增多，半潛船的需求量必然增大。

船舶機(jī)艙作為船舶主動(dòng)力驅(qū)動(dòng)裝置放置的空間，機(jī)艙通風(fēng)狀況的好壞對(duì)機(jī)艙內(nèi)設(shè)備能否正常運(yùn)行影響重大。船舶機(jī)艙通風(fēng)的目的就是建立并維持機(jī)艙內(nèi)適宜的環(huán)境條件(溫度、濕度、空氣流速、清潔度等),以保證柴油發(fā)電機(jī)、鍋爐燃燒所必須的空氣量，同時(shí)也要保證機(jī)艙內(nèi)部良好的工作環(huán)境，改善機(jī)艙工作人員正常工作所需的環(huán)境。曾宏強(qiáng)等將全新風(fēng)系統(tǒng)、循環(huán)冷卻加新風(fēng)系統(tǒng)及射流通風(fēng)3種方式進(jìn)行對(duì)比，分析各自的優(yōu)缺點(diǎn)；楊衛(wèi)國提出一種應(yīng)用CFD技術(shù)進(jìn)行艙內(nèi)大空間氣流組織分布預(yù)測(cè)以判定機(jī)艙大空間通風(fēng)系統(tǒng)完善程度的研究方法；郭昂等運(yùn)用CFD技術(shù)對(duì)海監(jiān)船機(jī)艙通風(fēng)系統(tǒng)進(jìn)行數(shù)值模擬，并對(duì)送回風(fēng)口形式進(jìn)行優(yōu)化，改善了機(jī)艙內(nèi)的氣流組織分布；趙楠等采用CFD仿真方法和通風(fēng)評(píng)價(jià)指標(biāo)，分析了機(jī)艙常規(guī)通風(fēng)時(shí)的氣流溫度場，并在此基礎(chǔ)上引入射流通風(fēng)，改善機(jī)艙的通風(fēng)狀況；CHEN等運(yùn)用CFD數(shù)值模擬方法研究機(jī)艙內(nèi)不同通風(fēng)方案下的風(fēng)速、溫度和濕度分布，選擇合適的通風(fēng)方案。

機(jī)艙內(nèi)的氣流組織是由風(fēng)口位置及風(fēng)量大小決定的，風(fēng)口高度的合理性會(huì)影響艙內(nèi)氣體的流動(dòng)性，減少渦旋的產(chǎn)生。本文以80 000噸級(jí)半潛船機(jī)艙為研究對(duì)象，應(yīng)用Airpak軟件對(duì)不同風(fēng)口布置高度的機(jī)艙通風(fēng)進(jìn)行數(shù)值模擬，找到合適的通風(fēng)管道布置高度，以維護(hù)機(jī)器設(shè)備的穩(wěn)定運(yùn)行，保障工作人員正常的工作環(huán)境。

1 機(jī)艙結(jié)構(gòu)及物理模型

1.1 機(jī)艙結(jié)構(gòu)

80 000噸級(jí)半潛船設(shè)有前后2個(gè)機(jī)艙，機(jī)艙尺寸為48 m×22.5 m×9.2 m(長×寬×高)，2個(gè)機(jī)艙通風(fēng)相互獨(dú)立。每個(gè)機(jī)艙分上下層，內(nèi)設(shè)2道縱壁，將每個(gè)機(jī)艙分為6塊區(qū)域。機(jī)艙內(nèi)有6臺(tái)發(fā)電機(jī)組、2臺(tái)鍋爐及電氣設(shè)備等。機(jī)艙通風(fēng)系統(tǒng)設(shè)計(jì)需充分考慮結(jié)構(gòu)的影響，確保機(jī)艙送風(fēng)及排風(fēng)順暢，使發(fā)電機(jī)組、鍋爐得到足夠的燃燒空氣，設(shè)備散發(fā)的熱量被充分帶走。

1.2 物理模型

根據(jù)機(jī)艙的實(shí)際結(jié)構(gòu)尺寸，將其簡化為48 m×22.5 m×9.2 m(長×寬×高)的長方體結(jié)構(gòu)，機(jī)艙內(nèi)的主要熱源是柴油機(jī)和鍋爐，而電氣設(shè)備等散熱相對(duì)于柴油機(jī)和鍋爐較少，可以忽略不計(jì)。根據(jù)機(jī)艙內(nèi)設(shè)備的實(shí)際結(jié)構(gòu)尺寸，在不影響模擬結(jié)果的前提下，將艙內(nèi)設(shè)備進(jìn)行簡化，艙內(nèi)設(shè)備及通風(fēng)系統(tǒng)的數(shù)量和結(jié)構(gòu)尺寸見表1。采用Airpak軟件建立機(jī)艙、艙內(nèi)設(shè)備及通風(fēng)系統(tǒng)的三維模型，見圖1。

表1 設(shè)備數(shù)量及尺寸

1—6 500 kW發(fā)動(dòng)機(jī)；2—4 200 kW發(fā)動(dòng)機(jī)；

2 數(shù)學(xué)模型

本文假設(shè)艙室內(nèi)的空氣流動(dòng)處于常溫、低速、不可壓的狀態(tài)，并符合理想氣體狀態(tài)方程和Boussinesq假設(shè)：

(1)密度變化除了對(duì)動(dòng)量守恒方程浮升力項(xiàng)有影響外，對(duì)其他項(xiàng)的影響可忽略，室內(nèi)空氣物性參數(shù)可視為常數(shù)。

(2)密度差與溫度差成比例，即：

(-)=(-)

(1)

綜合引入的湍流動(dòng)能方程和湍流動(dòng)能耗散率方程，得到的模擬室內(nèi)氣流組織的控制方程組如下：

(1)連續(xù)性方程

(2)

(2)動(dòng)量方程

(-)

(3)

(3)能量方程

(4)

(4)方程

(5)

(5)方程

(6)

式(1)～式(6)中：為參考密度；為流體密度；為容積膨脹系數(shù)，1/，取值為3.67×10；為流體溫度；為參考溫度；、為坐標(biāo)位置；、分別為速度矢量在、方向上的分量，=1,2,3；為湍流粘性系數(shù)；為粘性系數(shù)；為湍流有效粘性系數(shù)；為時(shí)間；為重力加速度；為流體壓力；為湍動(dòng)能；為湍流動(dòng)能耗散率；、、、、、為常數(shù)，經(jīng)驗(yàn)常數(shù)取值可參考表2。

表2 常用經(jīng)驗(yàn)常數(shù)

3 模擬計(jì)算

根據(jù)船舶設(shè)計(jì)規(guī)范，機(jī)艙的通風(fēng)量為柴油機(jī)鍋爐燃燒所需空氣量及排除機(jī)艙內(nèi)設(shè)備散熱所需空氣量之和。通過計(jì)算得到機(jī)艙的總通風(fēng)量為248 400 m/h，每個(gè)送風(fēng)口風(fēng)量為2.875 m/s，送風(fēng)風(fēng)速為2.995 m/s，送風(fēng)溫度為35 ℃。排風(fēng)量與送風(fēng)量相等，則每個(gè)排風(fēng)口風(fēng)量為11.5 m/s。6 500 kW發(fā)電機(jī)的散熱量約為170 kW，4 200 kW發(fā)電機(jī)的散熱量約為150 kW，1 900 kW發(fā)電機(jī)的散熱量約為100 kW，鍋爐的散熱量約為40 kW。機(jī)艙墻壁為絕熱壁面；排煙風(fēng)道由于保溫層隔熱作用，表面近似設(shè)置為恒溫60 ℃。

模型建立完成后，利用Airpak軟件對(duì)其進(jìn)行網(wǎng)格劃分和獨(dú)立性驗(yàn)證。監(jiān)測(cè)點(diǎn)溫度隨網(wǎng)格數(shù)量的變化見圖2。最終選擇設(shè)置各向最大網(wǎng)格尺寸，網(wǎng)格數(shù)量為282 046。入口設(shè)為等速度入口邊界，出口設(shè)為流量出口邊界條件，采用了壓力-速度耦合求解器、SIMPLE算法和二階迎風(fēng)離散方法，通過前期數(shù)值模型對(duì)比選擇了標(biāo)準(zhǔn)-湍流模型。

圖2 網(wǎng)格獨(dú)立性驗(yàn)證

4 模擬結(jié)果與分析

選取典型截面的溫度場、速度場，對(duì)3種不同風(fēng)口高度的機(jī)艙通風(fēng)系統(tǒng)進(jìn)行分析比較。

4.1 溫度場

設(shè)備表面溫度分布見圖3。從圖3可以看出，設(shè)備表面存在著明顯的溫度梯度，由下到上逐漸遞減，可見機(jī)艙通風(fēng)系統(tǒng)對(duì)艙內(nèi)設(shè)備表面有著明顯的散熱效果。但是，設(shè)備下表面的溫度依然較高，設(shè)備上表面溫度雖低，但也高于100 ℃，因此機(jī)艙通風(fēng)系統(tǒng)只能排除部分設(shè)備產(chǎn)生的熱量，剩余熱量會(huì)使艙內(nèi)溫度升高。根據(jù)設(shè)備表面溫度分布，如果不是必要情況，工作人員應(yīng)盡可能遠(yuǎn)離艙內(nèi)工作設(shè)備，避免造成不必要的危險(xiǎn)。圖3為不同送風(fēng)高度對(duì)機(jī)艙內(nèi)設(shè)備表面溫度分布的影響。從圖3可以看出，3種送風(fēng)高度下，設(shè)備表面溫度分布相差甚微，說明不同送風(fēng)高度對(duì)設(shè)備表面溫度分布影響不大。

圖3 設(shè)備表面溫度分布

圖4為截面=1.73 m處不同送風(fēng)高度的溫度分布，即工作人員正常工作時(shí)人體呼吸區(qū)所在截面的溫度分布。在該截面上，溫度分布都相對(duì)均勻，平均溫度40 ℃左右，低于柴油機(jī)驅(qū)動(dòng)船舶機(jī)艙要求設(shè)計(jì)的最高溫度50 ℃，人員在該環(huán)境下可以正常工作，不會(huì)出現(xiàn)明顯不適感。

圖4 Z=1.73 m處機(jī)艙平面溫度分布

風(fēng)口布置高度為6.5 m時(shí)，溫度分布均勻性明顯好于其他兩種風(fēng)口布置高度。出現(xiàn)這種狀況的原因是風(fēng)口高度為7.0 m時(shí)，此時(shí)距離截面=1.73 m較遠(yuǎn)，所在截面的溫度分布受設(shè)備散熱的影響較大；而由于排煙風(fēng)道的中心線高度為5.0 m，排煙風(fēng)管半徑為0.8 m，因此排煙風(fēng)口和高度為6.0 m送風(fēng)口距離過近，部分新風(fēng)剛從送風(fēng)口送入機(jī)艙，隨即又通過排煙風(fēng)道排出，造成了能源的浪費(fèi)，能源利用率降低，設(shè)備正常工作的散熱排出效率下降，不能有效改善機(jī)艙內(nèi)的通風(fēng)狀況。

4.2 速度場

圖5是截面=35 m處3種不同送風(fēng)高度的送風(fēng)口速度矢量圖。從圖中可以明顯看出，設(shè)備之間、設(shè)備與墻壁之間都產(chǎn)生了渦流。由于設(shè)備間的布置距離大小的關(guān)系，渦流數(shù)量和大小也相應(yīng)受影響。渦流的存在增強(qiáng)了設(shè)備間空氣的擾動(dòng)，有利于設(shè)備表面的散熱。

圖5 X=35 m處機(jī)艙平面速度分布

當(dāng)送風(fēng)風(fēng)口高度為7.0 m時(shí)，空氣形成的渦流主要分布在機(jī)艙的中部，設(shè)備間的擾動(dòng)相對(duì)于其他兩種送風(fēng)風(fēng)口高度的空氣擾動(dòng)較弱，設(shè)備的散熱能力也相應(yīng)較弱；當(dāng)送風(fēng)風(fēng)口高度為6.0 m和6.5 m時(shí)，渦流主要分布在機(jī)艙中下部,但是送風(fēng)風(fēng)口高度為6.5 m時(shí)，渦流區(qū)平均速度為0.5 m/s左右，設(shè)備間渦流的擾動(dòng)明顯更加劇烈，設(shè)備表面的散熱能力顯著，設(shè)備在此狀態(tài)下可以穩(wěn)定運(yùn)行，因此船上工作人員在此環(huán)境下工作，不適性會(huì)有所降低。

3種工況的模擬結(jié)果見表3。由表3可知，當(dāng)送風(fēng)口布置高度從7.0 m下降至6.0 m時(shí)，設(shè)備表面平均空氣速度逐漸降低，說明設(shè)備表面換氣效率有所提高；送風(fēng)口平均送風(fēng)速度和排風(fēng)口平均排氣速度相差不大，說明風(fēng)口布置的高度對(duì)送、排風(fēng)風(fēng)口速度影響不大。

表3 3種工況的模擬結(jié)果

5 結(jié)論

(1)應(yīng)用Airpak軟件對(duì)80 000噸級(jí)半潛船機(jī)艙通風(fēng)系統(tǒng)進(jìn)行數(shù)值模擬，通過對(duì)溫度場、速度場進(jìn)行分析，模擬結(jié)果可以對(duì)船舶機(jī)艙通風(fēng)系統(tǒng)的設(shè)計(jì)提供參考。

(2)模擬結(jié)果表明：送風(fēng)風(fēng)口高度布置為6.5 m時(shí)，機(jī)艙內(nèi)的速度、溫度分布都較為合理，機(jī)艙內(nèi)設(shè)備可以穩(wěn)定運(yùn)行，人員正常工作無明顯不適感。

(3)雖然在人員工作區(qū)域溫度分布對(duì)人員正常工作無明顯不適感，速度分布也趨于合理，但是設(shè)備表面溫度依然較高，長時(shí)間工作會(huì)使周圍環(huán)境迅速升溫，可以考慮在設(shè)備工作區(qū)域添加輔助散熱設(shè)備，加速設(shè)備散熱，保障設(shè)備穩(wěn)定運(yùn)行。