蓄冷材料應(yīng)用于艦船空調(diào)的數(shù)值仿真

李 鵬，吳贊俠

● （1.海軍駐上海江南造船（集團(tuán)）有限責(zé)任公司軍事代表室，上海 201913；2.舟山海軍監(jiān)修室，舟山 316000）

蓄冷材料應(yīng)用于艦船空調(diào)的數(shù)值仿真

李 鵬1，吳贊俠2

● （1.海軍駐上海江南造船（集團(tuán)）有限責(zé)任公司軍事代表室，上海 201913；2.舟山海軍監(jiān)修室，舟山 316000）

為了研究蓄冷材料在艦船空調(diào)系統(tǒng)中的應(yīng)用，分別對(duì)冰和某型高溫相變蓄冷球兩種材料的蓄釋冷過程用有限元分析方法進(jìn)行了模擬仿真，找出蓄冷量和蓄釋冷時(shí)間的變化規(guī)律，并分析比較兩種材料的蓄釋冷的性能，結(jié)果表明：采用高溫相變材料，其蓄冷周期短，制冷機(jī)組相對(duì)較節(jié)能，可滿足水面艦船蓄冷空調(diào)的要求。

蓄冷空調(diào)；冰蓄冷；高溫相變蓄冷材料；有限元法；數(shù)值模擬

0 引言

艦船使用蓄冷空調(diào)對(duì)于艦船節(jié)能和提高續(xù)航力意義重大，文獻(xiàn)[1][2]提出了在潛艇上使用蓄冷空調(diào)的方案，根據(jù)其研究結(jié)果，如在水面艦船應(yīng)用蓄冷空調(diào)將會(huì)提高艦船的戰(zhàn)技性能，水面艦船對(duì)設(shè)備重量和體積有著特定的要求，而蓄冷材料對(duì)蓄冷空調(diào)的性能有重要的影響，因此選擇何種蓄冷材料是較為重要的。

蓄冷空調(diào)目前主要為水蓄冷和冰蓄冷，對(duì)其他蓄冷材料如：共晶鹽蓄冷、氣體水合物蓄冷等國(guó)內(nèi)外均有研究。冰蓄冷在工程中應(yīng)用較為廣泛，優(yōu)點(diǎn)是其相變潛熱高，技術(shù)相對(duì)較為成熟，缺點(diǎn)是由于蒸發(fā)溫度低（-10℃～-8℃），使得空調(diào)的制冷系數(shù)COP和制冷量大幅度降低。而高溫相變材料蓄冷系統(tǒng)的相變溫度相對(duì)較高，在空調(diào)工況下，能較好地和蒸發(fā)溫度匹配，原有制冷機(jī)組及空調(diào)系統(tǒng)基本不需要改動(dòng)，而且冷損失小、制冷機(jī)COP值大，該蓄冷材料可滿足艦船蓄冷空調(diào)的要求。系統(tǒng)中，蓄冷材料常用的封裝形式是蓄冷球，下面對(duì)某種高溫相變材料[3]蓄冷球的蓄釋冷過程有限元分析方法進(jìn)行了模擬仿真，分析其蓄釋冷性能，研究在艦船上的適應(yīng)性。

1 建立數(shù)學(xué)模型

計(jì)算中采用直徑為50mm的蓄冷球，蓄冷時(shí)，低溫的冷媒水流過蓄冷球表面，吸收球內(nèi)的熱量，使得蓄冷材料由球壁向球內(nèi)逐漸凝固，直至整個(gè)蓄冷材料完全發(fā)生相變。釋冷時(shí)，空調(diào)流回的冷媒水流過蓄冷球表面吸收蓄冷球的熱量，蓄冷材料由球壁向球內(nèi)逐漸融化，直至完全發(fā)生相變。因傳熱過程較為復(fù)雜，同時(shí)存在導(dǎo)熱和對(duì)流換熱，各項(xiàng)物性參數(shù)會(huì)發(fā)生變化，為方便分析計(jì)算，考慮主要因素，所以對(duì)蓄冷材料及傳熱作如下假設(shè)：

1）蓄冷材料固液兩相物理性能不隨溫度變化，即比熱、導(dǎo)熱系數(shù)、密度均為常數(shù)。

2）熱量的傳遞以導(dǎo)熱為主，忽略對(duì)流換熱。

3）相變從外層到內(nèi)層發(fā)生，相界面在相變過程中溫度恒定。

根據(jù)上述假設(shè)，數(shù)學(xué)模型為：

初始條件：T（r , 0）=T0（0 ≤r≤R）

式中，T為蓄冷球的溫度，℃；t為蓄釋冷時(shí)間，s；Cps、Cp1為蓄冷材料的固、液相比熱，J／(kg·℃)；λs、λ1為蓄冷材料的固、液相導(dǎo)熱系數(shù)，w／(m·℃)；ρs、ρ1為蓄冷材料的固、液相密度，kg／m3；R為蓄冷球的半徑，m；L為蓄冷材料的固、液相變潛熱，J／kg；Tf為冷媒水的溫度，℃。

2 單元模型的建立和求解

采用ANSYS建立單元模型，基于有限元分析方法求解，并對(duì)結(jié)果進(jìn)行分析比較。

2.1 建立單元模型和參數(shù)的設(shè)置

單元模型見如圖1，其以1/4圓扇面作為分析模型，采用PLANE55單元類型，采用自由網(wǎng)格劃分，精度選擇為1，軸對(duì)稱方式求解。

圖1 單元模型

兩種蓄冷材料的物理性能參數(shù)如表1所示。

表1 兩種蓄冷材料參數(shù)

表2 兩種蓄冷材料溫度變化的焓值

2.2 加載求解

選用ANSYS的二維熱分析單元PLANE5進(jìn)行非線性的瞬態(tài)熱計(jì)算和分析，便捷條件和加載見表3。

表3 兩種蓄冷材料的初始值及載荷

選取第一類邊界條件，采用分析類型為“Transient”，對(duì)幾何模型施加溫度載荷并設(shè)定載荷選項(xiàng)，冰蓄釋冷和高溫相變材料蓄釋冷模擬時(shí)間分別為8000s和1000s；自動(dòng)時(shí)間步長(zhǎng)選項(xiàng)設(shè)定為on，時(shí)間步長(zhǎng)為20s，最小步長(zhǎng)、最大步長(zhǎng)分別為10s、30s，選用lines search，加載求解穩(wěn)定收斂。

2.3 數(shù)值仿真的結(jié)果

計(jì)算完成后進(jìn)行后處理，蓄冷過程的相變末期時(shí)刻云圖如圖2所示。蓄冷過程中兩種蓄冷材料球心處溫度和時(shí)間曲線見圖3。釋冷過程的相變末期時(shí)刻云圖如圖4所示。釋冷過程兩種蓄冷材料球心處溫度和時(shí)間曲線見圖5。

圖2 相變末期時(shí)刻云圖

圖3 兩種蓄冷材料球心處溫度和時(shí)間曲線

圖4 相變末期時(shí)刻云圖

圖5 兩種蓄冷材料球心處溫度和時(shí)間曲線

3 數(shù)值仿真的結(jié)果和分析

考慮艦船要求，蓄冷球體積相同更具有實(shí)際意義，在相同溫差條件下對(duì)兩種蓄冷材料進(jìn)行計(jì)算，用單個(gè)球蓄冷量進(jìn)行比較，仿真計(jì)算的結(jié)果及分析如下。

3.1 蓄冷過程的結(jié)果和分析

兩種蓄冷材料單個(gè)球蓄冷量、時(shí)間(蓄冷球從初始溫度降到低于相變溫度1℃的時(shí)間)、速度對(duì)比見表4。

表4 蓄冷過程性能對(duì)比

從表4可知，在相同體積條件下，冰蓄冷和高溫相變材料蓄冷相比，冰蓄冷的蓄冷量大，蓄冷時(shí)間長(zhǎng)，蓄冷速度低，這主要是因?yàn)楸罾錆摕岣?，但是?dǎo)熱率較低?？梢娫谛罾鋾r(shí)間相同的情況下，高溫相變材料的蓄冷量相對(duì)要更大。

3.2 釋冷過程的結(jié)果和分析

兩種蓄冷材料單個(gè)球釋冷量、時(shí)間(從初始溫度到高于相變溫度1℃的時(shí)間)、速度對(duì)比見表5。

表5 釋冷性能對(duì)比

從表5可知，在相同體積條件下，冰蓄冷和高溫相變蓄相比，冰蓄冷得冷釋冷量高，釋冷時(shí)間要長(zhǎng)，所以冰蓄冷的釋冷性好。但制冷機(jī)組在蓄冷工況時(shí)，蒸發(fā)溫度越低，制冷系數(shù)COP和制冷量越低，一般每下降1℃，COP下降3%左右，而冰蓄冷要比高溫相變材料低8℃左右，因此冰蓄冷相對(duì)耗電多近30%，制冷量低約一半，這將不利于艦船的節(jié)能和提高續(xù)航力。

4 結(jié)束語

1）綜上所述，與冰蓄冷相比，采用高溫相變材料，可提高制冷機(jī)組的制冷系數(shù)，同時(shí)由于和空調(diào)工況所需蒸發(fā)溫度相近，可使用單工況制冷機(jī)組，減少了設(shè)備的復(fù)雜性，有利于提高艦船的續(xù)航力。

2）盡管高溫蓄冷材料的蓄冷量小，但其蓄冷時(shí)間短，可以采用加快循環(huán)周期的方法彌補(bǔ)蓄冷量不足，以滿足艦艇空調(diào)的要求。

3）如何提高高溫相變蓄冷材料的潛熱以及延長(zhǎng)釋冷時(shí)間，達(dá)到進(jìn)一步節(jié)能，是有待解決的關(guān)鍵技術(shù)。

[1] 曾強(qiáng)洪, 吳鋼. 蓄冷技術(shù)應(yīng)用于常規(guī)潛艇的方案設(shè)計(jì)[J]. 制冷空調(diào), 2006, 34(11): 53-56.

[2] 吳鋼, 曾強(qiáng)洪. 利用蓄冷空調(diào)提高潛艇水下續(xù)航力研究[J]. 船海工程, 2006(6).

[3] 羅智特, 杜雁霞, 賈代勇. 相變材料蓄冷球蓄冷過程的數(shù)值模擬[J]. 制冷空調(diào), 2004(2).

Numerical Simulation of Cooling Storage Material Using in Warship’s Air Conditioning

LI Peng1, WU Zan-xia2
(1. Naval Military Representative Department of Jiangnan Shipyard (group) Co., Ltd, Shanghai 201913, China;2. Naval Military Representative Department of Zhongshan Ship Repairing, Zhoushan 316000, China)

For the study of application of cooling storage material in air conditioning system of warship, an ANSYS program is used for numerical simulation of ice storage ball and a kind of high temperature phase-change cooling storage ball in cooling storage and melting processes. The cooling storage capacity, cooling storage and melting duration of the two materials are obtained in cooling storage and melting processes. The cooling storage and melting performances of the two materials are compared. The result of the research shows that the high temperature phase-change material has shorter cooling storage cycle and better chillers performance in energy saving, and the material satisfies the request of warship cooling storage air conditioning。

cooling storage air-conditioning; ice cooling storage; high temperature phase-change cooling storage material; Finite Element Method (FEM); numerical simulation

U664

A

李鵬（1968-），男，工程師、碩士。研究方向：:輪機(jī)工程。