吳小舟,趙加寧,王灃浩

(1.西安交通大學(xué) 人居環(huán)境與建筑工程學(xué)院，陜西 西安　710049；2.哈爾濱工業(yè)大學(xué) 市政環(huán)境工程學(xué)院，黑龍江 哈爾濱　150090)

?

基于形狀因子的輻射地板傳熱量計算等效熱阻模型*1

吳小舟1,2?,趙加寧2,王灃浩1

(1.西安交通大學(xué) 人居環(huán)境與建筑工程學(xué)院，陜西 西安710049；2.哈爾濱工業(yè)大學(xué) 市政環(huán)境工程學(xué)院，黑龍江 哈爾濱150090)

對輻射地板傳熱過程進行分析，提出了基于形狀因子的輻射地板傳熱量計算等效熱阻模型.為了驗證簡化模型的合理性，采用數(shù)值模擬計算結(jié)果進行驗證，并與標(biāo)準(zhǔn)手冊中采用的簡化模型計算結(jié)果進行對比分析.結(jié)果表明，當(dāng)管間距變化范圍為50～250mm,換熱管上端填充層厚度變化范圍為25～65mm及熱水平均溫度變化范圍為25～45 ℃時，ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型計算結(jié)果與數(shù)值模擬計算結(jié)果相對誤差均較大，最大分別為20.2%，30.4%和22.8%，而本文提出的基于形狀因子的等效熱阻模型計算結(jié)果與數(shù)值模擬計算結(jié)果相對誤差較小，最大不超過3%；當(dāng)管間距變化范圍為50～200mm,換熱管上端填充層厚度變化范圍為15～55mm及冷水平均溫度變化范圍為10～20 ℃時，ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型計算結(jié)果與數(shù)值模擬計算結(jié)果相對誤差也均較大，最大分別為80.1%，17.7%和16.8%，而本文提出的基于形狀因子的等效熱阻模型計算結(jié)果與數(shù)值模擬計算結(jié)果相對誤差較小，最大不超過2%.

輻射地板；傳熱量計算；等效熱阻模型；形狀因子

輻射地板供暖供冷系統(tǒng)不僅具有較高熱舒適性，而且能夠顯著地降低建筑運行能耗，在公共建筑及住宅中得到廣泛的應(yīng)用[1-3].輻射地板傳熱量計算模型是輻射地板供暖供冷系統(tǒng)設(shè)計和運行控制的關(guān)鍵，因此國內(nèi)外學(xué)者對輻射地板傳熱量計算模型進行了大量的研究，并提出了很多模型.根據(jù)輻射地板傳熱過程描述的差異，這些模型可以分為3類：解析模型[4-6]、數(shù)值模擬模型[7-9]和簡化模型[10-12].

解析模型結(jié)構(gòu)相對復(fù)雜，不便于工程師及設(shè)計師理解應(yīng)用.而數(shù)值模擬模型僅能得到分散的結(jié)果，由于輻射地板組成部分及其特性參數(shù)可選擇的組合較多，導(dǎo)致數(shù)值模擬工作量相當(dāng)大，很難包括所有的組合.因此，結(jié)構(gòu)相對簡單及計算時間相對較少的簡化模型被廣泛應(yīng)用于輻射地板設(shè)計手冊及標(biāo)準(zhǔn)中，如ISO標(biāo)準(zhǔn)中的冪函數(shù)計算模型[13]、ASHRAE手冊中的平面肋片模型[14]及地暖設(shè)計手冊中的等效熱阻模型[15]等.但這些簡化模型均存在一定問題，計算精度較低或者使用不方便.如ISO標(biāo)準(zhǔn)冪函數(shù)計算模型需要查表二次計算，過程比較繁瑣；ASHRAE手冊平面肋片模型忽略了地板垂直方向的傳熱，只考慮水平方向的傳熱，可能會引起較大的計算誤差；地暖設(shè)計手冊等效熱阻模型形式過于簡單，計算缺乏理論依據(jù).

本文通過對輻射地板傳熱過程進行分析，提出計算精度較高且結(jié)構(gòu)相對簡單的輻射地板傳熱量計算簡化模型，并與數(shù)值模擬模型及標(biāo)準(zhǔn)手冊中采用的模型計算結(jié)果進行對比分析，為輻射地板供暖供冷系統(tǒng)的設(shè)計及運行控制提供理論方法.

1　輻射地板供暖供冷系統(tǒng)傳熱量計算模型

輻射地板供暖供冷系統(tǒng)是將換熱管埋設(shè)于建筑圍護結(jié)構(gòu)(樓板)中的供暖供冷系統(tǒng)，其結(jié)構(gòu)原理如圖1所示.

圖1　輻射地板結(jié)構(gòu)原理圖Fig.1　Schematic diagram of floor heatingsystem with embedded tube

如圖1所示，輻射地板主要由地面層、填充層、換熱管、保溫層和結(jié)構(gòu)層組成.其中λco為地面層的導(dǎo)熱系數(shù)；λscr為填充層的導(dǎo)熱系數(shù)；λp為換熱管的導(dǎo)熱系數(shù)；λins為保溫層的導(dǎo)熱系數(shù)；λcon為結(jié)構(gòu)層的導(dǎo)熱系數(shù)；δco為地面層的厚度；δp為換熱管的厚度；δins為保溫層的厚度；δcon為結(jié)構(gòu)層的厚度；M為換熱管間距；Do為換熱管外徑；Di為換熱管內(nèi)徑；d1為換熱管上端的填充層厚度；d2為換熱管下端的填充層厚度.

1.1基于形狀因子的等效熱阻模型

由圖1可知，輻射地板傳熱過程為：熱(冷)水以對流換熱方式把熱(冷)量傳遞給管內(nèi)壁，然后管內(nèi)壁以導(dǎo)熱方式把熱(冷)量傳遞到管外壁，之后管外壁以導(dǎo)熱方式通過填充層、保溫層或結(jié)構(gòu)層向地板表面?zhèn)鳠幔詈蟮匕迳匣蛳卤砻嬉詫α髋c輻射傳熱方式把熱(冷)量傳遞到室內(nèi).對于輻射地板供暖供冷系統(tǒng)，其傳熱量的計算難點主要在于計算管外壁面與填充層外壁面之間的導(dǎo)熱量.這主要是由于輻射地板中與管外壁接觸的填充層部分形狀比較復(fù)雜，溫度及熱流密度在二個坐標(biāo)方向上都是變化的，如此采用分析法求解計算通常相當(dāng)麻煩和困難.

隨著建筑節(jié)能技術(shù)的不斷發(fā)展，輻射供暖供冷系統(tǒng)熱(冷)媒溫度與室內(nèi)溫度相差較小，此時管外壁面與填充層外壁面表面溫度可視為均勻分布的，即兩者均可作為等溫界面.如此管外壁面與填充層外壁面之間的導(dǎo)熱量計算可以采用基于形狀因子的半無限大物體與其平行的圓管排之間(如圖2所示)導(dǎo)熱計算公式[16]，其形狀因子如式(1)所示.

圖2　無限大平板中的圓排管Fig.2　Round pipe in the infinity plane

(1)

基于上述分析，本文提出了基于形狀因子的輻射地板傳熱量計算等效熱阻模型，其地板上表面?zhèn)鳠崃坑嬎闶饺缡?2)所示.

(2)

式中: qf為地板上表面?zhèn)鳠崃?，W/m2；tw為管內(nèi)平均熱水溫度，℃；top為室內(nèi)操作溫度，℃；供暖時為20 ℃，供冷時為26 ℃；R為管內(nèi)熱水到地板上表面室內(nèi)空氣之間總的傳熱熱阻，(m2K)/W，忽略接觸熱阻如式(3)所示：

(3)

式中：αf為地板上表面綜合換熱系數(shù)(W/(m2K))，對于輻射供冷熱系統(tǒng)，一般等于 10.8W/(m2K) ，對于輻射供冷系統(tǒng)，一般等于 6.5W/(m2K)[13]；αw為換熱管內(nèi)表面對流換熱系數(shù)(W/(m2K))；Req為基于形狀因子的地板填充層等效熱阻((m2K)/W)，如式(4)所示：

(4)

式中：l為軸向管長度，m.

輻射地板填充層中的圓排管如圖3所示，其填充層等效熱阻可視為圖2中所述熱阻的2倍，因此其形狀因子可如式(5)所示.

圖3　輻射地板中的圓排管Fig.3　Round pipe in the radiant heating plane

(5)

將式(1)及式(5)代入式(4)，可得輻射地板填充層的等效熱阻，如式(6)所示：

(6)

假設(shè)地板上、下表面?zhèn)鳠釡夭钕嗟?，已知地板上表面?zhèn)鳠崃亢蛡鳠釤嶙钑r，可得到地板下表面的傳熱量，如式(7)所示：

(7)

式中：qc為地板下表面?zhèn)鳠崃?，W/m2；Ru為換熱管到地板上表面室內(nèi)之間總的傳熱熱阻，m2K/W，如式(8)所示；Rd為換熱管到地板下表面室內(nèi)之間總的傳熱熱阻，m2K/W，如式(9)所示：

(8)

(9)

式中：αc為地板下表面綜合換熱系數(shù)，W/(m2k)；對于輻射供熱系統(tǒng)，一般等于 6.5W/(m2K) ，對于輻射供冷系統(tǒng)，一般等于10.8W/(m2K)[13].

1.2數(shù)值模擬計算模型

上節(jié)通過對輻射地板傳熱過程進行分析，提出了基于形狀因子的輻射地板傳熱量計算等效熱阻模型.為了驗證簡化模型的合理性，本文將采用與數(shù)值模擬計算結(jié)果進行對比的方式.這主要是由于輻射地板傳熱量的數(shù)值模擬計算結(jié)果可信度較高，相關(guān)計算結(jié)果已被輻射供暖供冷標(biāo)準(zhǔn)采納[13].

1．2．1輻射地板物理模型

相對管間地板溫差傳熱，忽略地板軸向的溫差傳熱能得到足夠的計算精度，輻射地板三維導(dǎo)熱問題一般可以簡化為二維導(dǎo)熱問題，如圖4所示.

圖4　二維地板簡化模型Fig.4　Simplified two-dimension model of radiant floor

1．2．2輻射地板數(shù)學(xué)模型

假設(shè)各層材料為各向同性，互相緊密接觸，忽略接觸熱阻，則描述該計算單元的二維穩(wěn)態(tài)導(dǎo)熱微分

方程式如式(10)所示：

(10)

1．2．3定解條件

1)單元對稱邊界

對每一個計算單元對稱溫度場邊界有：

(11)

2)地板上表面邊界條件

地板上表面以對流換熱和輻射換熱兩種方式向房間傳熱，將對流換熱系數(shù)和輻射換熱系數(shù)綜合成一個換熱系數(shù)，如式(12)所示：

(12)

3)地板下表面邊界條件

地板下表面也是以對流換熱和輻射換熱兩種方式向房間傳熱，將對流換熱系數(shù)和輻射換熱系數(shù)綜合成一個換熱系數(shù)，如式(13)所示：

(13)

二維穩(wěn)態(tài)導(dǎo)熱微分方程組邊界條件復(fù)雜，分析求解困難，一般都采用數(shù)值求解.HEAT2軟件采用有限差分法求解二維穩(wěn)態(tài)導(dǎo)熱問題，具有計算精度較高、耗時較短等優(yōu)點，目前已被輻射供暖供冷國際標(biāo)準(zhǔn)ISO11855采用[13].因此，本文采用HEAT2軟件對上述導(dǎo)熱微分方程進行求解計算[17].

2　結(jié)果與分析

根據(jù)國際標(biāo)準(zhǔn)ISO11855[13]，典型輻射地板供暖供冷系統(tǒng)結(jié)構(gòu)及傳熱參數(shù)如表1所示.

為了將本文提出的基于形狀因子的等效熱阻模型與數(shù)值模擬計算模型及目前標(biāo)準(zhǔn)規(guī)范及手冊中采用的模型計算結(jié)果進行對比分析，以下將分別采用數(shù)值模擬計算模型、ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型計算不同工況下輻射地板傳熱量.

2.1輻射供暖傳熱量計算結(jié)果與分析

當(dāng)管間距變化范圍為50～250mm,換熱管上端填充層厚度變化范圍為25～65mm及熱水平均溫度變化范圍為25～45 ℃時，不同輻射地板供暖傳熱量各計算方法的計算結(jié)果如表2～表4所示.

表2　不同管間距下輻射地板傳熱量計算結(jié)果(d1=45 mm，tw=35 ℃)Tab.2　Calculation results of different heat transfer calculation models(d1=45 mm，tw=35 ℃)

注ε為輻射地板傳熱量計算結(jié)果與數(shù)值模擬計算結(jié)果的計算相對誤差.

表3　不同換熱管上端填充層厚度下輻射地板傳熱量計算結(jié)果(M=150 mm,tw=35 ℃)Tab.3　Calculation results of different heat transfer calculation models(M=150 mm，tw=35 ℃)

表4　不同熱水平均溫度下輻射地板傳熱量計算結(jié)果(M=150 mm，d1=45 mm)Tab.4　Calculation results of different heat transfer calculation models(M=150 mm，d1=45 mm)

由表2知，當(dāng)換熱管上端填充層厚度為 45mm，熱水平均溫度為35 ℃及管間距變化范圍為50～250mm時，管間距的變化對ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差基本都超過5%，最大達(dá)到30%.而管間距的變化對基于形狀因子的熱阻模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過3%.

由表3可知，當(dāng)管間距為150mm，熱水平均溫度為35 ℃及換熱管上端填充層厚度變化范圍為25～65mm時，換熱管上端填充層厚度的變化對ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差都超過5%，最大達(dá)到20%.而換熱管上端填充層厚度的變化對基于形狀因子的熱阻模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過3%.

由表4可知，當(dāng)管間距為150mm、換熱管上端填充層厚度為45mm及熱水平均溫度變化范圍為25～45 ℃時，熱水平均溫度的變化 對ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差都超過5%，最大達(dá)到17%.而換熱管上端填充層厚度的變化對基于形狀因子的熱阻模型及ISO標(biāo)準(zhǔn)冪函數(shù)計算模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過3%.

2.2輻射地板供冷傳熱量計算結(jié)果與分析

當(dāng)管間距變化范圍為50～200mm，換熱管上端填充層厚度變化范圍為15～55mm及冷水平均溫度變化范圍為10～20 ℃時,不同輻射地板供冷傳熱量各計算方法的計算結(jié)果如表5～表7所示.

表5　不同管間距下輻射地板傳熱量計算結(jié)果(d1=35 mm，tw=15 ℃)Tab.5　Calculation results of different heat transfer calculation models(d1=35 mm，tw=15 ℃)

表6　不同換熱管上端填充層厚度下輻射地板傳熱量計算結(jié)果對比(M=100 mm，tw=15 ℃)Tab.6　Calculation results of different heat transfer calculation models(M=100 mm，tw=15 ℃)

表7　不同冷水平均溫度下輻射地板傳熱量計算結(jié)果(M=100 mm，d1=35 mm)Tab.7　Calculation results of different heat transfer calculation models(M=100 mm，d1=35 mm)

由表5知，當(dāng)換熱管上端填充層厚度為35mm，冷水平均溫度為15 ℃及管間距變化范圍為50～200mm時，管間距的變化對ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差基本都超過5%，最大達(dá)到50%.而管間距的變化對基于形狀因子的熱阻模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過2%.

由表6可知，當(dāng)管間距為100mm，冷水平均溫度為15 ℃及換熱管上端填充層厚度變化范圍為15～55mm時，換熱管上端填充層厚度的變化對ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差都超過5%，最大達(dá)到80%.而換熱管上端填充層厚度的變化對基于形狀因子的熱阻模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過2%.

由表7可知，當(dāng)管間距等于100mm，換熱管上端填充層厚度等于35mm及冷水平均溫度變化范圍為10～20 ℃時，冷水平均溫度的變化對ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型結(jié)果影響較大，這些方法的計算結(jié)果與數(shù)值模擬計算結(jié)果的相對誤差都超過5%，最大達(dá)到44%.而換熱管上端填充層厚度的變化對基于形狀因子的等效熱阻模型及ISO標(biāo)準(zhǔn)冪函數(shù)計算模型計算結(jié)果影響較小，與數(shù)值模擬計算結(jié)果的相對誤差最大不超過2%.

對比表2～4及表5～7可知，本文提出的基于形狀因子的等效熱阻模型在供冷工況下的計算誤差明顯小于供熱工況.這主要是因為供冷工況相比供熱工況，供水平均溫度更接近于室內(nèi)溫度，導(dǎo)致管外壁面與填充層外壁面表面溫度分布相對更均勻，與模型假設(shè)一致，因此模型計算誤差相對較小.

根據(jù)本文的計算結(jié)果可得到本文提出的基于形狀因子的等效熱阻模型與其他簡化模型的對比結(jié)果，如表8所示.

表8　不同簡化模型對比Tab.8　Comparison of different heat transfercalculation models

由表8可知，本文提出的基于形狀因子的等效熱阻模型計算精度較高，使用范圍較廣.一般而言，換熱管上端填充層厚度均大于水管的半徑，因此本文提出的基于形狀因子的等效熱阻模型可用于常規(guī)的輻射地板供暖供冷系統(tǒng)設(shè)計及運行調(diào)節(jié).

本文提出的基于形狀因子的等效熱阻模型使用過程與其他簡化模型基本一致.已知房間冷熱負(fù)荷及供回水溫度條件下，通過調(diào)整管間距及換熱管上端填充層厚度等使設(shè)計的輻射地板供熱供冷量等于或略大于房間冷熱負(fù)荷，此時的管間距及換熱管上端填充層厚度等即為輻射地板設(shè)計參數(shù).

3　結(jié)　論

本文通過對輻射地板傳熱過程進行分析，提出了合理的基于形狀因子的輻射地板傳熱量計算等效熱阻模型，并與標(biāo)準(zhǔn)手冊中采用的模型計算結(jié)果進行對比分析，得到以下結(jié)論：

1)ISO標(biāo)準(zhǔn)冪函數(shù)計算模型、ASHRAE手冊平面肋片模型及地暖設(shè)計手冊等效熱阻模型的計算精度均較低；

2)本文提出的基于形狀因子的等效熱阻模型不僅計算精度較高，而且形式相對簡單.

[1]RICHARDW,CHAPMANK.Radiantheatingandcoolinghandbook[M].NewYork:McGrawHillProfessional, 2002:12-15．

[2]JANB,OLESENBW,KIMKW.Lowtemperatureheatingandhightemperaturecooling[M].Brussets:REHVA, 2009：24-28．

[3]ROBERTB,OLESENBW,KIMKW.Historyofradiantheating&coolingsystems[J].ASHRAEJournal, 2010, 52(1):26-31.

[4]MARKUSK,LEHMANNB.Thermoaktivebauteilsystemetabs[M]．Zurich:EMPAEnergiesysteme/Haustechnik,2000：122-131.(InGerman)

[5]LUX,TERVOLAP,Transientheatconductioninthecompositeslab-analyticalmethod[J]．JournalofPhysicsA:MathematicalandGeneral, 2005, 38 (1) :81-86.

[6]ABDELAZIZL．Developmentofaradiantheatingandcoolingmodelforbuildingenergysimulationsoftware[J]．BuildingandEnvironment, 2004, 39 (4): 421-431.

[7]李安邦,徐新華.內(nèi)嵌管式輻射地板的頻域熱特性分析[J].湖南大學(xué)學(xué)報：自然科學(xué)版,2015,42(1):115-119．

LIAn-bang,XUXin-hua.Studyofthefrequencythermalcharacteristicsofpipe-embeddedradiantfloors[J]．JournalofHunanUniversity:NaturalSciences, 2015,42(1):115-119.(InChinese)

[8]SILVANAFL,FILIPPNC,LESINOG.Transientsimulationofastoragefloorwithaheating/coolingparallelpipesystem[J].BuildingSimulation, 2010, 3 (2) :105-115.

[9]JINXin,ZHANGXiao-song.Numericalsimulationofradiantfloorcoolingsystem:theeffectsofthermalresistanceofpipeandwatervelocityontheperformance[J]．BuildingandEnvironment, 2010,45 (11) :2545-2552.

[10]LIUYan-feng,WANGDeng-jia,LIUJia-ping.Studyonheattransferprocessforin-slabheatingfloor[J].BuildingandEnvironment, 2012, 54(10): 77-85.

[11]ZHANGLun,LIUXiao-hua,JIANGYi.Simplifiedcalculationforcooling/heatingcapacity,surfacetemperaturedistributionofradiantfloor[J].EnergyandBuildings, 2012, 55(12):397-404．

[12]LI,Qing-qing.Simplifiedthermalcalculationmethodforfloorstructureinradiantfloorcoolingsystem[J].EnergyandBuildings, 2014, 74(7): 182-190．

[13]ISO/DIS11855-2Buildingenvironmentdesign-design,constructionandoperationofradiantheatingandcoolingsystems-part2:determinationofthedesignheatingandcoolingcapacity[S].Brussels:InternationalOrganizationforStandardization,2010:7-9．

[14]HANDBOOKA.HVACsystemsandequipment[M].Atlanta,GA:AmericanSocietyofHeating,RefrigeratingandAirConditioningEngineers, 2011:11-13.

[15]魯翠.地板采暖設(shè)計施工手冊[M].北京：中國電力出版社，2009:16-19.

LUCui.Designhandbookforfloorheatingconstruction[M].Beijing:ChinaElectricPowerPress, 2009:16-19．(InChinese)

[16]錢濱江,伍貽文,常加芳,等.簡明傳熱學(xué)手冊[M].北京:高等教育出版社,1983:22-24.

QIANBin-jiang,WUYi-wen,CHANGJia-fang, et.al.Conciseheatconductionhandbook[M].Beijing:HigherEducationPress,1983:22-24.(InChinese)

[17]BLOMBERGT.Heat2-AheattransferPC-program:manualforheat2[D].Lund:DepartmentofBuildingPhysics,University,Lund, 1991: 22-25.

An Equivalent Thermal Resistance Model to Calculate the Heat Transfer of Radiant Floor Heating and Cooling Systems Using Shape Factor

WU Xiao-zhou1,2?, ZHAO Jia-ning2, WANG Feng-hao1

(1.School of Human Settlements and Civil Engineering, Xi'an Jiaotong Univ, Xi’an, Shaanxi710049, China;2.School of Municipal and Environmental Engineering, Harbin Institute of Technology, Harbin，Heilongjiang150090, China)

Anequivalentthermalresistancemodelfortheheattransfercalculationofradiantfloorbasedonshapefactorwasdevelopedinthisstudy.Theproposedmodelwasverifiedbynumericalsimulation,andcomparedwiththedatafromtheuniversalsinglepowerfunctionofISOstandard,thefinmodelinASHAREhandbook,andtheequivalentthermalresistancemodelindesignhandbook.TheheattransfersobtainedbythenumericalsimulationdisagreedsignificantlywiththeexistingmethodsincludingtheuniversalsinglepowerfunctionofISOstandard,thefinmodelinASHAREhandbook,andtheequivalentthermalresistancemodelindesignhandbook.Themaximumerrorrateswere20.2 %, 30.4 %,and22.8 %,respectively,whenthetubespacerangesfrom50to250mm,thethicknessoffilllayerabovepiperangesfrom25to65mm,andtheaveragehotwatertemperaturerangesfrom25to45 ℃.Ontheotherhand,themaximumdifferenceoftheheattransferpredictionsbetweentheproposedmodelandthesimulationsoftwarewaslessthan3%.Whenthetubespacerangesfrom50to200mm,thethicknessoffilllayerabovepiperangesfrom15to55mm,andtheaveragecoldwatertemperaturerangesfrom10to20 ℃,themaximumerrorratesoftheheattransferpredictionsbetweenthenumericalsimulationandtheexistingmethodswere80.1 %, 17.7 %,and16.8 %.Ontheotherhand,theheattransferpredictedbytheproposedmodelshowedlessthan2 %ofdifferencefromthatofthenumericalsimulation.

radiantfloor;heattransfercalculation;equivalentthermalresistancemodel;shapefactor

1674-2974(2016)05-0137-07

2015-05-08

國家自然科學(xué)基金資助項目(51408482)，National Natural Science Foundation of China(51408482); 西安市科技計劃項目(CXY1514(3))

吳小舟(1984-)，男，江西九江人，西安交通大學(xué)副教授，博士

?通訊聯(lián)系人，E-mail:fonen519@mail.xjtu.edu.cn

TU833.1

A