李偉光 李安邦 施紅梅 徐新華 謝軍龍

船舶圍壁空氣層自然對流對傳熱系數(shù)的影響分析

李偉光1李安邦2施紅梅1徐新華2謝軍龍2

1中國艦船研究設(shè)計(jì)中心
2華中科技大學(xué)建筑環(huán)境與設(shè)備工程系

本文采用數(shù)值模擬軟件對不同室內(nèi)外溫度條件下的船舶水平及垂直圍壁穩(wěn)態(tài)傳熱進(jìn)行模擬計(jì)算。結(jié)果表明水平圍壁在冬季和夏季的傳熱系數(shù)存在明顯的差異，而且冬季的傳熱系數(shù)要比夏季的傳熱系數(shù)大的多，這主要是水平空氣層在冬季的自然對流比夏季的自然對流強(qiáng)烈。垂直圍壁冬季和夏季傳熱系數(shù)的差異很小。結(jié)果也表明同一季節(jié)的室內(nèi)外溫差的變化對圍壁傳熱系數(shù)的影響均很小。因此在計(jì)算船舶空調(diào)負(fù)荷時(shí)，水平圍壁傳熱系數(shù)取值應(yīng)對冬季和夏季予以區(qū)分而垂直圍壁的傳熱系數(shù)選取不需區(qū)分冬季和夏季。

船舶圍壁自然對流數(shù)值計(jì)算傳熱系數(shù)

0引言

近幾十年來，隨著我國船舶行業(yè)發(fā)展，船舶行業(yè)總體能耗也隨之增大，在國家大力發(fā)展節(jié)能減排技術(shù)的政策背景下，研究船舶節(jié)能技術(shù)具有很大的現(xiàn)實(shí)意義。船舶空調(diào)系統(tǒng)是船舶主要的耗能設(shè)施之一，萬噸級(jí)以上遠(yuǎn)洋船舶空調(diào)系統(tǒng)能耗約占船舶總能耗的近1/3[1]。在船舶空調(diào)系統(tǒng)設(shè)計(jì)中，往往由于對艙室內(nèi)空調(diào)負(fù)荷的保守計(jì)算（甚至是估算）導(dǎo)致空調(diào)設(shè)備選型偏大，既加大了初投資也造成了運(yùn)行過程中能量的浪費(fèi)，因此獲得正確船舶艙室空調(diào)負(fù)荷是提高船舶空調(diào)系統(tǒng)經(jīng)濟(jì)性和實(shí)現(xiàn)空調(diào)系統(tǒng)節(jié)能的前提。

船舶艙室一般采用穩(wěn)態(tài)的方法計(jì)算空調(diào)冷熱負(fù)荷。船舶圍壁傳熱產(chǎn)生的空調(diào)負(fù)荷與船舶圍壁的傳熱系數(shù)成正比，其取值直接影響到艙室空調(diào)負(fù)荷計(jì)算結(jié)果。船舶圍壁結(jié)構(gòu)包括鋼板、梁（扶材）、隔熱保溫層、裝飾板，其中隔熱保溫層與裝飾板之間存在一定厚度的空氣層。船舶圍壁空氣層在溫差（密度差）的驅(qū)動(dòng)下會(huì)產(chǎn)生自然對流從而強(qiáng)化換熱，不同溫度條件下的自然對流換熱強(qiáng)度會(huì)有差異[2]。

目前船舶空調(diào)負(fù)荷計(jì)算時(shí)圍壁傳熱系數(shù)的選取一般參考文獻(xiàn)[3～4]。文獻(xiàn)中提供的船舶圍壁傳熱系數(shù)簡化算法以及傳熱系數(shù)參考值均沒有考慮船體圍壁空氣層自然對流對船舶圍壁傳熱系數(shù)的影響。本文采用數(shù)值模擬方法對船舶圍壁進(jìn)行穩(wěn)態(tài)傳熱模擬，對比分析在不同溫度條件（自然對流強(qiáng)度）下船舶圍壁傳熱系數(shù)的差異。

1船舶圍壁傳熱數(shù)值模擬方法

1.1物理模型以及邊界條件

船舶圍壁按幾何位置可以分為水平（甲板）圍壁和垂直圍壁兩種。文獻(xiàn)[4]中列舉了多種不同隔熱材料、尺寸的船舶圍壁隔熱結(jié)構(gòu)形式，本文摘取其中典型的一種有空氣層的隔熱結(jié)構(gòu)形式進(jìn)行分析研究。該種隔熱結(jié)構(gòu)組成及尺寸（單位：mm）見圖1，各部分材料的物性參數(shù)見表1。由于水平圍壁中空氣層主要的自然對流運(yùn)動(dòng)方向?yàn)榭諝鈱雍穸确较?，圍壁的長度對自然對流的影很小可以忽略，而垂直圍壁中空氣層自然對流的主要運(yùn)動(dòng)方向?yàn)榭諝鈱痈叨确较颍瑖诘母叨葘ψ匀粚α鞯挠绊戄^大須考慮[2]。因此選取結(jié)構(gòu)如圖1所示的750mm長度（水平圍壁的梁間距為750mm）的水平圍壁、3000mm高度的垂直圍壁（一般艙室的高度為2600mm到3000mm）作為船舶圍壁傳熱的物理模型。船舶圍壁內(nèi)外表面與空氣的對流換熱系數(shù)分別為8W/(m2·K)與80W/(m2·K)，夏季艙室內(nèi)空氣溫度為27℃，冬季艙室內(nèi)空氣溫度為22℃[4]，夏季艙室外空氣溫度取30～65℃（下文中將通過改變艙室外的空氣溫度來分析船舶圍壁的傳熱系數(shù)變化），冬季艙室外空氣溫度取-16～19℃。水平（垂直）圍壁物理模型的左右（上下）兩側(cè)壁面均視作絕熱面。FLUENT軟件是一款基于有限體積法的數(shù)值計(jì)算軟件，在求解流-固耦合傳熱問題方面已有較多的應(yīng)用[5～7]。有空氣層的船舶圍壁傳熱是典型的流-固耦合傳熱問題，本文采用FLUENT軟件對船舶圍壁進(jìn)行穩(wěn)態(tài)傳熱模擬。

圖1一種有空氣層的船舶圍壁隔熱結(jié)構(gòu)

表1隔熱結(jié)構(gòu)各部分材料的熱物性參數(shù)

1.2控制方程及數(shù)值計(jì)算模型

運(yùn)用質(zhì)量守恒定律、牛頓第二定律以及熱力學(xué)第一定律可以建立關(guān)于圍壁傳熱的一系列控制方程：固體區(qū)域的二維導(dǎo)熱微分方程，氣體區(qū)域的二維連續(xù)性方程、動(dòng)量方程及能量守恒方程。本文采用Boussinesq假設(shè)忽略流體中粘性耗散，除動(dòng)量方程中浮升力項(xiàng)考慮密度的變化外，其余各項(xiàng)中密度以及其它物性參數(shù)均看作用常數(shù)[8]。其中氣體區(qū)域x、y方向上的動(dòng)量方程分別如式（1）、（2），式（2）中ρβ(T-T0)g為基于Boussinesq假設(shè)的浮升力項(xiàng)，式（3）為空氣層自然對流的瑞利數(shù)。

式中：T為溫度變量；T0是參考溫度；u、v分別是速度在x、y方向上的分量；μ是空氣動(dòng)力粘性系數(shù)；λ是空氣導(dǎo)熱系數(shù)；cρ是空氣定壓比熱；ρ是空氣密度；β為空氣的熱膨脹系數(shù)；g為重力加速度；△T為空氣層上下壁面的溫差；L為空氣層的厚度；α為熱擴(kuò)散率。

在FLUENT中求解時(shí)，壓力一速度耦合方程采用SIMPLE算法，動(dòng)量、能量方程選擇二階迎風(fēng)格式。由式（3）計(jì)算出本文所研究的空氣層自然對流的瑞利數(shù)Ra均小于108，因此使用層流模型。改變船舶室內(nèi)外的空氣溫度（室內(nèi)空氣溫度恒定），通過FLUENT軟件模擬計(jì)算得出不同室內(nèi)外溫差下船舶圍壁的溫度分布、傳熱量及傳熱系數(shù)。

2結(jié)果比較分析

2.1水平圍壁

表2為水平圍壁在不同室內(nèi)外空氣溫差下的傳熱系數(shù)計(jì)算結(jié)果。在相同大小的室內(nèi)外溫差下水平圍壁在冬季的傳熱系數(shù)比夏季的傳熱系數(shù)大很多，前者甚至是后者的幾倍之多。夏季水平圍壁的空氣層冷面在下方而熱面在上方，空氣直接分層，自然對流很弱，空氣層的傳熱過程可近似地看做純導(dǎo)熱[2]。夏季水平圍壁的溫度分布云圖如圖2所示。冬天水平圍壁的冷面在上方而熱面在下方，下方靠近熱面的空氣在浮升力的作用下上升，上方靠近冷面的空氣被冷卻而下降，這樣就形成了在封閉的空氣層內(nèi)循環(huán)的自然對流，對流換熱大大地減少了空氣層部分的傳熱熱阻，因此水平圍壁在冬天的換熱系數(shù)一般要比夏天的換熱系數(shù)要大很多。冬季水平圍壁的溫度分布云圖如圖3所示。

表2不同室內(nèi)外空氣溫度差下的水平船舶圍壁傳熱系數(shù)

圖2夏季水平圍壁的溫度分布云圖（室內(nèi)外溫差18℃）

圖3冬季水平圍壁的溫度分布云圖（室內(nèi)外溫差-18℃）

表2還表明，改變冬夏季室內(nèi)外溫差大小對水平圍壁的傳熱系數(shù)的影響均很小。夏季水平圍內(nèi)空氣層的傳熱方式為導(dǎo)熱，改變溫差對圍壁內(nèi)空氣層的導(dǎo)熱熱阻的幾乎沒有影響，因此夏季室內(nèi)外溫差的變化對船舶圍壁的傳熱系數(shù)影響很小。冬季水平圍壁內(nèi)空氣層自然對流換熱較大程度地減小了空氣層熱阻，使得空氣層部分的熱阻相比隔熱材料等其它部分的熱阻小很多?？諝鈱由舷卤诿娴臏夭钫颊麄€(gè)室內(nèi)外溫差的比例很小，因此當(dāng)室內(nèi)外溫差有較大的變化時(shí)，空氣層上下壁面的溫差其實(shí)變化很小從而對流換熱的強(qiáng)度的改變也不明顯，因此冬天室內(nèi)外溫差的變化對船舶圍壁的傳熱系數(shù)的影響很小。

2.2垂直圍壁

表3為垂直船舶圍壁在不同室外溫差下的傳熱系數(shù)計(jì)算結(jié)果。對于垂直圍壁，在相同大小的室內(nèi)外溫差下冬季和夏季船舶圍壁的傳熱系數(shù)相差很小，且冬夏季室內(nèi)外溫差的大小改變對垂直圍壁傳熱系數(shù)的影響均很小。垂直圍壁空氣層兩側(cè)存在溫差，熱側(cè)的空氣因浮升力上升而冷側(cè)的空氣遇冷下降，這樣就在空氣層內(nèi)形成了自然對流。圖4與圖5分別為垂直圍壁在夏季及冬季的溫度分布云圖。垂直圍壁內(nèi)的空氣層在冬夏季均能形成自然對流，因此在大小相同的室內(nèi)外溫差下冬夏季垂直圍壁的傳熱系數(shù)相差很小?？諝鈱幼匀粚α鱾鳠岬臒嶙柘啾扔趪诟魺岵牧系绕渌糠值臒嶙栊『芏啵?dāng)室內(nèi)外溫差有較大的變化時(shí)，空氣層兩側(cè)壁之間的溫差變化實(shí)際很小。因此冬夏季室內(nèi)外溫差大小的變化對船舶垂直圍壁的傳熱系數(shù)的影響均很小。

表3不同室內(nèi)外空氣溫度差的垂直船舶圍壁傳熱系數(shù)

圖4夏季垂直圍壁溫度分布云圖（室內(nèi)外溫差為18℃）

圖5冬季垂直圍壁溫度分布云圖（室內(nèi)外溫差為-18℃）

3總結(jié)與討論

本文采用FLUENT軟件對船舶圍壁進(jìn)行了穩(wěn)態(tài)傳熱模擬計(jì)算，研究討論了不同的室內(nèi)外溫度條件下圍壁空氣層對圍壁傳熱系數(shù)的影響。由于水平圍壁的空氣層在夏季及冬季形成的空氣自然對流完全不一樣，水平外圍壁在冬季及夏季的傳熱系數(shù)差異很大。垂直圍壁的空氣層在夏季及冬季都能形成良好自然對流，在冬季及夏季的傳熱系數(shù)基本一樣。在冬夏季室內(nèi)外溫差大小的變化對水平圍壁和垂直圍壁的傳熱系數(shù)的影響均很小。因此在計(jì)算船舶艙室冬夏季空調(diào)負(fù)荷時(shí)，水平圍壁的傳熱系數(shù)在冬季和夏季應(yīng)取不同的值，而垂直圍壁的傳熱系數(shù)可以取相同的值。

[1]劉曉晨.船舶空調(diào)節(jié)能措施探討[J].青島遠(yuǎn)洋船員學(xué)院學(xué)報(bào), 2007,28(2):29-30

[2]章熙民,任澤霈,梅飛鳴.傳熱學(xué)(第五版)[M].北京:中國建筑工業(yè)出版社,2007

[3]《船舶設(shè)計(jì)實(shí)用手冊冷藏通風(fēng)》編寫組.船舶設(shè)計(jì)實(shí)用手冊冷藏通風(fēng)分冊[M].上海:國防工業(yè)出版社,1975

[4]中華人民共和國國家技術(shù)監(jiān)督局．船舶起居處所空氣調(diào)節(jié)與通風(fēng)設(shè)計(jì)參數(shù)和計(jì)算方法(GB/T 13409-92)[S].北京:中國標(biāo)準(zhǔn)出版社,1992

[5]A Gustavsen,B T Griffith,D Arasteh.Natural convection effects in three-dimensional window frames with internal cavities[J]. ASHRAE Transactions,2001,107(2):527-537

[6]張琪,王偉旭,黃蘇榮,等.高密度車用永磁電機(jī)流固耦合傳熱仿真分析[J].電機(jī)與控制應(yīng)用,2012,39(8):1-5

[7]劉硯文,李振海.入口氣流方向?qū)?zhǔn)逆流式全熱交換器性能的影響[J].建筑節(jié)能,2012,40(1):14-17

[8]陶文銓.數(shù)值傳熱學(xué)(第2版)[M].西安:西安交通大學(xué)出版社, 2005

Ana lys is of the Effe c ts of Na tura l Conve c tion in the Air La ye r of Ship Enve lope s on the He a t Tra ns fe r Coe ffic ie nt

LI Wei-guang1,LI An-bang2,SHI Hong-Mei1,XU Xin-hua2,XIE Jun-long2
1 China Ship Development and Design Center
2 Department of Building Environment and Energy Engineering,Huazhong University of Science and Technology

Numerical simulation software is used to simulate the steady heat transfer of ship envelopes under various indoor and outdoor temperature conditions.Results show that the heat transfer coefficient of horizontal ship envelops in winter season differs much from the one in summer season.The coefficient in winter season is much larger than that in summer season due to that the natural ventilation is much powerful in winter season.Results also show that the heat transfer coefficient of vertical ship envelops in summer season is almost the same to that in winter season due to similar natural convection.The simulation results also indicate that the heat transfer coefficient of ship envelops is almost the same in one season although the indoor and outdoor temperature difference changes.Therefore,when calculating the air conditioning cooling/heating load of ships,the heat transfer coefficient of horizontal ship envelops in winter and summer should be determined separately while the heat transfer coefficient of vertical ship envelops in winter and summer can be the same.

ship envelope,natural convection,numerical simulation,heat transfer coefficient

1003-0344（2014）05-038-4

2013-7-8

李偉光（1982～），男，碩士，工程師；湖北省武漢市中國艦船研究設(shè)計(jì)中心（430064）；E-mail:lwgabe@gmail.com

中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金（2013QN072）；教育部高校博士點(diǎn)基金（201201421 10078）