姜 輝，曾 驥，黃朝春 ，劉 瑋，張?zhí)焯?，牟星?/p>

（上海海事大學(xué) 海洋工程技術(shù)研究中心，上海 201306）

0 引言

極地旅游的發(fā)展帶動(dòng)了極地小型郵輪的需求，但在惡劣的極地環(huán)境下，油箱內(nèi)的燃油粘度增加，甚至結(jié)冰，可能導(dǎo)致設(shè)備無(wú)法正常工作，目前絕大部分的船舶上采用的艙柜加熱形式都是蒸汽盤(pán)管加熱，這種加熱形式使用方便、技術(shù)成熟、熱效率較高，也存在燃油存儲(chǔ)艙表面散熱損失大、盤(pán)管破損修復(fù)困難、表面易積碳等缺點(diǎn)[1]。如何通過(guò)合理的系統(tǒng)設(shè)計(jì)，在極地郵輪的營(yíng)運(yùn)過(guò)程中將上述不利影響降至最低，以及采取更為有效的方法，減少燃油艙加熱過(guò)程中的熱量損失，對(duì)極地郵輪安全性方面的保障具有重要意義。

影響蒸汽盤(pán)管加熱效率的因素很多，單因素的對(duì)于加熱效率存在一般性規(guī)律，而在綜合各因素對(duì)于加熱效率存在非線性關(guān)系，因此，需要通過(guò)數(shù)值模擬研究通過(guò)非線性擬合找到設(shè)計(jì)參數(shù)和計(jì)算目標(biāo)之間的關(guān)系，為極地船舶的應(yīng)用提供了指導(dǎo)和設(shè)計(jì)依據(jù)。

目前，針對(duì)燃油艙加熱的研究較少，主要是針對(duì)運(yùn)輸船的原油貨艙的加熱，與本文研究類似，可以通用。如黎俊杰[2]對(duì)105000 DWT油輪NO.4貨艙的加熱過(guò)程的流場(chǎng)進(jìn)行數(shù)值模擬計(jì)算，確定油輪貨艙加熱盤(pán)管加熱效率的影響因素。Magazinovi?等[3]對(duì) 65000 DWT油輪NO.3油艙的貨油加熱過(guò)程進(jìn)行熱損失研究。Suhara等[4]研究33000 DWT 油輪的No.5貨艙加熱過(guò)程中的熱損失。Akagis等[5]采用傳熱學(xué)經(jīng)驗(yàn)公式及研究方法計(jì)算得到蒸汽加熱盤(pán)管的導(dǎo)熱效率和艙室艙壁的熱損失情況。本文以極地郵輪的T064號(hào)燃油艙柜為研究對(duì)象，對(duì)于燃油艙采用的蒸汽盤(pán)管防凍措施進(jìn)行仿真計(jì)算，分析防凍效果的單因素影響，符合參考文獻(xiàn)[6]的一般規(guī)律，再結(jié)合優(yōu)化算法，得出最后優(yōu)化的方案設(shè)計(jì)，提高了加熱效率，具有實(shí)際的應(yīng)用價(jià)值。

1 研究對(duì)象及方案

1.1 研究對(duì)象

以極地小型郵輪的T064號(hào)燃油艙（尺寸約為7 m×12 m×2.6 m）為研究對(duì)象，對(duì)蒸汽盤(pán)管的加熱效率進(jìn)行分析，艙柜大小為7 m×12 m×2.6 m，該模型外部為厚度20 mm的艙柜框架，內(nèi)部為燃油，在燃油艙內(nèi)部靠近船底500 mm處布置有蛇形蒸汽盤(pán)管（見(jiàn)圖1），盤(pán)管內(nèi)徑為10 mm，長(zhǎng)度共計(jì)約2000 mm，在簡(jiǎn)化模型時(shí)，直接將盤(pán)管在燃油區(qū)域進(jìn)行布爾減運(yùn)算。

圖1 極地小型郵輪燃油艙模型

1.2 研究方案

通過(guò)對(duì)燃油艙蒸汽盤(pán)管的加熱模擬分析，得出溫度與時(shí)間的變化關(guān)系，再研究單因素對(duì)于加熱效果的一般規(guī)律與參考文獻(xiàn)的一般規(guī)律是否一致，進(jìn)而綜合所有因素進(jìn)行優(yōu)化計(jì)算，得到最優(yōu)的蒸汽盤(pán)管結(jié)構(gòu)參數(shù)，研究方案流程如圖2所示。

圖2 研究方案流程圖

2 溫度場(chǎng)模擬計(jì)算

溫度場(chǎng)模擬計(jì)算需要假設(shè)以下3種條件：

1）所有溫度場(chǎng)計(jì)算都是基于盤(pán)管溫度已經(jīng)達(dá)到穩(wěn)態(tài)，忽略蒸汽和盤(pán)管間的加熱過(guò)程；

2）僅考慮燃油滿艙的狀態(tài)，不考慮蒸發(fā)和使用等影響燃油體積的因素，即同時(shí)忽略船舶晃蕩的影響；

3）艙壁與外界海水及空氣都是全接觸，不考慮環(huán)境和外部自由液面的影響。

實(shí)驗(yàn)需測(cè)取同種激勵(lì)不同電場(chǎng)條件下經(jīng)減振器減振后的加速度衰減的數(shù)據(jù)。實(shí)驗(yàn)采用型號(hào)為JZ-5激振器給減振器提供正弦力F，采用DH5299動(dòng)態(tài)信號(hào)采集器和GF-20型號(hào)的功率放大器對(duì)激振器的振幅A和頻率ω進(jìn)行控制以及采集各傳感器的數(shù)據(jù)，加速度傳感器型號(hào)是YJ9A的壓電傳感器；壓電力傳感器的型號(hào)為YFF-1。按照上述的實(shí)驗(yàn)器材以及實(shí)驗(yàn)想法設(shè)計(jì)出實(shí)驗(yàn)系統(tǒng)示意圖3，其中虛線框中的直流穩(wěn)壓電源和高壓放大器組成了直流穩(wěn)壓高壓電源，為后續(xù)對(duì)比實(shí)驗(yàn)提供1000V的穩(wěn)定高壓。

2.1 網(wǎng)格劃分

本文網(wǎng)格大小選用200 mm×200 mm，網(wǎng)格劃分后的燃油艙柜溫度場(chǎng)模型，網(wǎng)格質(zhì)量分布如圖3所示?？梢钥吹?，大部分網(wǎng)格質(zhì)量接近1，且平均值為0.72。

圖3 溫度場(chǎng)網(wǎng)格質(zhì)量分布

2.2 熱邊界條件

燃油艙在蒸汽盤(pán)管加熱的同時(shí)也伴隨著向外界環(huán)境散熱的過(guò)程，這種有熱量逸散的加熱過(guò)程被稱為非穩(wěn)態(tài)傳熱過(guò)程。在對(duì)燃油艙做有限元溫度場(chǎng)計(jì)算時(shí)，根據(jù)熱力學(xué)計(jì)算的3種邊界條件，應(yīng)用至本文，可以得到3種類型的壁面邊界條件，分別是：

1）規(guī)定燃油艙壁面的溫度值，使其成為恒溫度壁面條件；

2）規(guī)定燃油艙壁面的熱流密度值，使其成為恒定熱流壁面，當(dāng)改值為0時(shí)，壁面被認(rèn)為絕熱；

3）規(guī)定燃油艙壁面的對(duì)流換熱情況，這種情況認(rèn)為整個(gè)熱力場(chǎng)中的固液體具有對(duì)流換熱的現(xiàn)象，需要提前設(shè)置目標(biāo)溫度和對(duì)流換熱系數(shù)。

選擇第3種邊界條件，根據(jù)國(guó)家標(biāo)準(zhǔn)《油艙蒸汽加熱系統(tǒng)計(jì)算方法》（CB/T 3373-2013）確定蒸汽盤(pán)管的傳熱系數(shù)、散熱面（壁面）傳熱系數(shù)、油溫和環(huán)境溫度。

熱邊界條件由內(nèi)到外分別是燃油與艙壁鋼板的邊界、艙壁與壓載艙海水的邊界和艙壁與臨艙空氣的邊界，導(dǎo)熱系數(shù)分別為116.3 W/m2·℃， 17.445 W/m2·℃和4.652 W/m2·℃，同時(shí)把燃油和艙壁的交界面定義為固液分界面，燃油艙蒸汽和盤(pán)管的導(dǎo)熱系數(shù)為116.3 W/m2·℃。計(jì)算忽略了蒸汽盤(pán)管和蒸汽間的接觸面，以蒸汽溫度100℃作為蒸汽管道和燃油的邊界條件。

2.3 計(jì)算結(jié)果及分析

如圖4和圖5所示，艙柜與海水相鄰的左右2個(gè)面上的溫度顯然小于其他4個(gè)面，且4個(gè)頂角的溫度處于最低的52.692℃。在艙柜下表面由于更加靠近盤(pán)管的位置，高溫分布呈現(xiàn)和盤(pán)管布置方向相似的圖形。

圖4 艙柜上頂面穩(wěn)態(tài)溫度分布云圖

圖5 艙柜下底面穩(wěn)態(tài)溫度分布云圖

圖6 燃油上部分穩(wěn)態(tài)溫度分布云圖

圖7 燃油下部分穩(wěn)態(tài)溫度分布云圖

可知，盤(pán)管在加熱燃油時(shí)，會(huì)出現(xiàn)低溫區(qū)域死角的現(xiàn)象，若燃油柜的臨艙是壓載艙，海水艙等注水艙室，則更需要注意這些近水艙壁的溫度分布，提前做好保溫措施，防止燃油艙柜內(nèi)部熱量逸散過(guò)多。

在瞬態(tài)溫度場(chǎng)計(jì)算模塊下，可以看到艙柜內(nèi)溫度隨時(shí)間變化逐漸升高，根據(jù)曲線呈凸函數(shù)的變化規(guī)律可知，加熱的效率隨時(shí)間逐漸下降，在加熱約8000 s后，艙柜內(nèi)的最小溫度值基本保持不變，即在8000 s左右，艙柜已基本完成加熱，溫度恒定在50℃左右。

3 蒸汽盤(pán)管單因素布置方案及效果

3.1 蒸汽盤(pán)管的布置間距

在以蒸汽盤(pán)管間的布置間距為因素進(jìn)行計(jì)算時(shí)，保持盤(pán)管管徑10 mm和總長(zhǎng)度3.5 m不變，且對(duì)蛇形盤(pán)管做一定的簡(jiǎn)化，簡(jiǎn)化為5根平行放置的蒸汽盤(pán)管，此處提取了不同盤(pán)管間距D的燃油區(qū)域最低溫度，如圖8所示。

圖8 不同間距布置下的燃油區(qū)域溫度分布

可知：1）隨著盤(pán)管布置間距的增大，燃油區(qū)域的最低溫度升高，主要是因?yàn)槿加蛢?nèi)部的熱源越來(lái)越分散，高溫區(qū)域逐漸擴(kuò)大，因此燃油的邊角區(qū)域溫度上升；2）隨著布置間距的增加，相鄰兩盤(pán)管間的高溫區(qū)域溫度逐漸減低；3）雖然模型做了一些簡(jiǎn)化，但溫度的基本分布規(guī)律與蛇形盤(pán)管一致，都呈現(xiàn)高溫區(qū)域集中在相鄰盤(pán)管間，低溫區(qū)域分散在艙室內(nèi)的各個(gè)頂角上。

3.2 蒸汽盤(pán)管的長(zhǎng)度

在以蒸汽盤(pán)管的總長(zhǎng)度為因素進(jìn)行計(jì)算時(shí)，保持盤(pán)管管徑10 mm和布置間距1800 mm不變，且同3.1節(jié)一樣，對(duì)蛇形盤(pán)管做同樣的簡(jiǎn)化。此處提取了盤(pán)管不同總長(zhǎng)度L時(shí)燃油區(qū)域最低溫度，如圖9所示。

圖9 不同總長(zhǎng)度時(shí)燃油區(qū)域的最低溫度變化云圖

可知，當(dāng)盤(pán)管總長(zhǎng)度減小時(shí)，燃油區(qū)域的最低溫度也隨之降低，但相對(duì)于布置間距對(duì)溫度的影響，顯然總長(zhǎng)度對(duì)其影響更小一些。

3.3 蒸汽盤(pán)管的管徑

在以蒸汽盤(pán)管管徑為因素進(jìn)行計(jì)算時(shí)，保持布置間距和總長(zhǎng)度不變，分別計(jì)算不同管徑的加熱效率和效果，不同管徑時(shí)燃油區(qū)域的最低溫度，如圖10所示。

圖10 不同管徑時(shí)燃油區(qū)域的最低溫度

可以看出，在不同管徑的條件下，燃油區(qū)域在加熱至最終狀態(tài)時(shí)，溫度都在50℃～55℃變化，且有管徑越大，最低溫度越大的變化趨勢(shì)。從高溫區(qū)域的分布情況來(lái)看，顯然隨著管徑的變大，高溫區(qū)域的體積也逐漸變大，尤其是相鄰管路之間的高溫區(qū)；對(duì)于低溫區(qū)而言，則是面積逐漸變小，且往蒸汽管末端角落逐漸收縮，最終低溫區(qū)域集中在該處。

另外，從實(shí)船的角度來(lái)說(shuō)，管徑越大，則需要的蒸汽量就越大，對(duì)于船用鍋爐的工作符合要求就更高。因此，如要加快加熱速率，管徑并不是越大越好，而是要將它與其他參數(shù)一起考慮計(jì)算才能得到最優(yōu)值。

4 基于響應(yīng)面優(yōu)化法的盤(pán)管布置方案分析

4.1 BBD響應(yīng)面優(yōu)化法應(yīng)用

根據(jù)BBD三因素設(shè)計(jì)方案原理，軟件對(duì)3個(gè)因素分別進(jìn)行自適應(yīng)地設(shè)計(jì)安排，管徑從160 ～240 mm變化，單根盤(pán)管長(zhǎng)度為5600 ～8400 mm變化，盤(pán)管間距為1280～1920 mm變化。

4.2 燃油區(qū)域最低溫度響應(yīng)

經(jīng)過(guò)計(jì)算，任意2組參數(shù)對(duì)目標(biāo)參數(shù)的響應(yīng)面如圖11～圖13所示?？梢钥吹?，3種因素對(duì)燃油區(qū)域的最低溫度影響變化規(guī)律基本與第2部分所分析的結(jié)果一致，即間距越長(zhǎng)，總長(zhǎng)度越大，管徑越小，最低溫度越小的變化規(guī)律。

圖11 總長(zhǎng)度和管徑與最低溫度的響應(yīng)面

圖12 總長(zhǎng)度和盤(pán)管間距與最低溫度的響應(yīng)面

圖13 管徑和盤(pán)管間距與最低溫度的響應(yīng)面

根據(jù)計(jì)算結(jié)果顯示，當(dāng)管徑為200 mm，總長(zhǎng)為3.5 m，布置間距為1600 mm時(shí)，最低溫度達(dá)到最優(yōu)值50.1℃。

5 結(jié)語(yǔ)

本文以極地小型郵輪的T046號(hào)燃油艙柜為研究對(duì)象，對(duì)加熱盤(pán)管參數(shù)對(duì)加熱效率的影響——做了模擬實(shí)驗(yàn)和定量分析，單從數(shù)值模擬結(jié)果的表面來(lái)看，得出間距越長(zhǎng)，總長(zhǎng)度越大，管徑越小，最低溫度越小的變化規(guī)律；根據(jù)單因素實(shí)驗(yàn)得到的結(jié)論，結(jié)合響應(yīng)面分析理論，針對(duì)蒸汽盤(pán)管的布置間距、蒸汽盤(pán)管的長(zhǎng)度和蒸汽盤(pán)管的管徑進(jìn)行實(shí)驗(yàn)設(shè)計(jì)并進(jìn)行仿真模擬計(jì)算，根據(jù)計(jì)算結(jié)果建立數(shù)學(xué)模型，得到了加熱效率最高的貨艙加熱盤(pán)管結(jié)構(gòu)參數(shù)，達(dá)到了縮短加熱時(shí)間，提高加熱效率的目的。