呂留根 黃 晨 陳劍昌 李 麗 白天宇

(上海理工大學環(huán)境與建筑學院 上海 200093)

基于實驗的不同風量下輻射時間因子修正

呂留根 黃 晨 陳劍昌 李 麗 白天宇

(上海理工大學環(huán)境與建筑學院 上海 200093)

利用輻射時間因子，可方便計算輻射得熱到輻射負荷的轉(zhuǎn)化過程。在變風量工況下，實測環(huán)境室24個時刻的逐時得熱和負荷，采用對流輻射分離方法計算24個時刻逐時輻射得熱和輻射負荷，整理出相應(yīng)的傳遞函數(shù)，通過構(gòu)造系數(shù)矩陣，計算得到24項輻射時間因子。研究結(jié)果表明：提出的輻射時間因子計算方法合理可行；輻射時間因子第一項值隨送風風量降低而減小，前五項值之和大于0.75，對于由ASHRAE軟件計算的輻射時間因子，對前五項進行修正，即可用于不同空氣循環(huán)類型房間輻射負荷的計算。

輻射得熱；輻射負荷；輻射時間因子；修正方法

輻射時間序列法[1-2]是ASHRAE TC4.1(設(shè)計負荷計算委員會)推薦的新的空調(diào)冷負荷計算方法。該方法將房間得熱劃分為對流和輻射兩部分，對流部分直接成為冷負荷；而輻射部分采用24項輻射時間因子將其轉(zhuǎn)化為空調(diào)負荷(稱為輻射負荷)，意義明確，避免了傳遞函數(shù)法計算輻射負荷的迭代計算，應(yīng)用簡便[3-4]。應(yīng)用輻射時間序列方法計算輻射負荷，首先要確定圍護結(jié)構(gòu)的輻射時間因子，輻射時間因子一般由輻射傳遞函數(shù)系數(shù)通過矩陣運算得到[5-6]。輻射傳遞函數(shù)的計算有兩種方法：1)由熱平衡方法計算得到，ASHRAE以此原理已建立典型區(qū)域的房間輻射傳遞函數(shù)系數(shù)數(shù)據(jù)庫[7-9]，并編制了輻射時間因子的計算軟件[5,10]，但該方法計算的輻射時間因子僅考慮了圍護結(jié)構(gòu)的蓄熱特性，并未體現(xiàn)送風風量變化的影響。2)實測逐時輻射得熱和輻射負荷，通過對實測數(shù)據(jù)進行數(shù)理分析得到，單寄平[11]通過實驗研究了房間太陽輻射得熱和冷負荷的傳遞函數(shù)；連之偉等[12]采用實驗研究了熱源得熱和冷負荷的房間傳遞函數(shù)，但并未區(qū)分熱源得熱的對流部分和輻射部分，且兩者均未計算圍護結(jié)構(gòu)的輻射時間因子。

本文提出了一種由實驗得到輻射時間因子的計算方法，并針對ASHRAE提供軟件計算的輻射時間因子未考慮送風風量變化的不足，研究了輻射時間因子隨送風風量改變的變化規(guī)律，并給出了考慮風量的修正方法。

1 輻射時間因子的計算

1.1 輻射時間因子與傳遞函數(shù)系數(shù)關(guān)系

輻射時間序列法在計算輻射負荷時，τ時刻的輻射負荷與τ時刻的輻射得熱及之前諸時刻的輻射得熱有關(guān)，采用相應(yīng)時刻的輻射時間因子來反應(yīng)對τ時刻的輻射負荷的影響。在時間輻射序列法中，逐時輻射負荷計算式為[5]：

qτ=r1Qτ+r2Qτ-Δτ+r3Qτ-2Δτ+…+r24Qτ-23Δτ

(1)

式中：qτ為τ時刻輻射負荷，W；Qτ-nτ為n時刻前的得熱量的輻射部分，W；rn為第n個輻射時間因子。應(yīng)用輻射時間序列方法計算輻射負荷，式(1)寫成矩陣形式為：

(2)

該矩陣可簡寫為：

q=RQ

(3)

式中：q為逐時輻射負荷列向量；Q為逐時輻射得熱列向量；R為輻射時間因子矩陣。

以房間輻射得熱Qτ為擾量，輻射負荷qτ為響應(yīng)，利用房間輻射傳遞函數(shù)，則輻射負荷為[13-15]：

qτ=ν0Qτ+ν1Qτ-1-ω1qτ-1

(4)

式中：νi、ωi為輻射傳遞函數(shù)系數(shù)，式(4)寫成矩陣形式為：

(5)

將式(5)矩陣形式進行簡化，并進行矩陣變換為：

(6)

式中：C1為ν系數(shù)矩陣；C2為ω系數(shù)矩陣。由式(3)和式(6)可知，采用傳遞函數(shù)法和輻射時間序列法，其輻射負荷計算結(jié)果基本相等，為此可得輻射時間因子為：

(7)

1.2 傳遞函數(shù)系數(shù)的實驗求解方法

建筑內(nèi)壁面是一種導熱、對流和輻射共存在復(fù)合換熱過程。采用對壁面進行對流和輻射的分離方法，則內(nèi)壁面i存在：

Qiλ+Qid+QiR=0

(8)

式中：Qiλ為導熱傳熱量，W；Qid為對流換熱量，W；QiR為輻射換熱量，W。輻射換熱量QiR可由有效輻射模型計算，即由各壁面溫度、角系數(shù)和壁面發(fā)射率聯(lián)立矩陣求得壁面有效輻射，然后再列矩陣求解壁面輻射換熱量如下[16]：

(9)

(10)

式中：εi為壁面發(fā)射率；Ai為壁面面積，m2；σ為黑體輻射常數(shù)，5.67×10-8W/(m2·K4)；θi為壁面溫度，K；Ji為壁面有效輻射，W/m2；Xij為角系數(shù)。實驗中用電熱膜散熱量模擬壁面導熱傳熱，其傳熱量可實測得到；壁面輻射換熱量由式(9)和式(10)計算；然后可由式(8)計算壁面對流換熱量。實驗中，房間輻射得熱為除加熱面以外的各個壁面的凈輻射得熱量之和，也等于加熱面的凈輻射失熱量。輻射負荷是指除加熱面以外其他壁面通過對流輻射分離方法分離出來的對流熱量之和。

由上述方法可得環(huán)境室24個時刻逐時輻射得熱和輻射負荷。利用最小二乘法，令傳遞函數(shù)系數(shù)滿足式(11)，并使M值最小，從而求出相應(yīng)的輻射傳遞函數(shù)系數(shù)值[11-12]。

(11)

2 實驗系統(tǒng)與參數(shù)測量

2.1 實驗室概況與測試方案

實驗系統(tǒng)由環(huán)境室、冷熱源、空氣處理系統(tǒng)、電氣控制與測試系統(tǒng)等部分組成。環(huán)境室尺寸為2.5 m×2.1 m×2.5 m，采用100 mm厚聚氨酯保溫庫板制作，詳見圖1。該環(huán)境室位于周邊環(huán)境溫度可控的恒溫房間內(nèi)，環(huán)境室外間見圖2。

實驗時溫度測點布置見圖3。實驗中共布有溫度測點27個，其中北墻電熱膜加熱，布有溫度測點7個；東墻、西墻、北墻、屋頂、地面各布有溫度測點2個。送回風風口各有溫度測點1個；室內(nèi)共8個測點測量空氣溫度，布有4根測線，每根測線有2個測點，距地面距離分別為0.8 m和1.6 m。熱流密度測點布置見圖4。實驗中布有熱流測點9個，其中北墻4個，其余5個壁面各1個。

2.2 測試儀器與測試工況

實驗中，空氣溫度、送回風溫度和壁面溫度通過溫度傳感器測量，壁面熱流密度則通過熱流傳感器測量，兩者通過建通采集儀(JTNT-C)采集，送回風風量采用風量罩測量，實驗中使用的測試儀器見表1。實驗前對各測試儀器進行了實驗室標定。為模擬周期擾量下非穩(wěn)態(tài)實驗工況，實驗中壁面電熱膜采用周期加熱方式，采用正弦擾量的均值和峰值見表2。

3 實驗結(jié)果與分析

3.1 傳遞函數(shù)系數(shù)的求解

圖5～圖7反映的是三個工況下輻射得熱和輻射負荷的關(guān)系，由圖可知：隨著時間的增加，輻射得熱和輻射負荷均是先增后減，輻射負荷相對于輻射得熱存在峰值上的衰減和時間上的延遲現(xiàn)象。具體來看，三種工況的峰值延遲時間均為2個時刻。且隨著送風風量減少，衰減倍數(shù)逐漸增加，工況1衰減倍數(shù)為1.20，工況2為1.24，工況3為1.28。分析其原因，這是因為在相同等級熱流密度下，送風風量越小，壁面處風速減小，壁面對流換熱減少，轉(zhuǎn)化成為室內(nèi)冷負荷的比例減少，導致衰減倍數(shù)的增加。

對工況1、工況2和工況3，將逐時輻射得熱和輻射負荷數(shù)據(jù)整理為傳遞函數(shù)三項式，見表3。由傳遞函數(shù)系數(shù)需滿足ν0+ν1-ω1=1，將求和誤差均分到傳遞函數(shù)三個系數(shù)，三個工況的系數(shù)值進行調(diào)整后使其滿足系數(shù)和為1的條件。傳遞函數(shù)ν0值表示某時刻的輻射得熱轉(zhuǎn)化為當時刻冷負荷的比例，從表3可知，隨著送風風量的減少，ν0值有所減小，這同樣是因為壁面處風速降低導致對流換熱能力減弱引起的。

說明[14-15]：中等風量(M)——通過地板、墻壁或天花板散流器送風；大風量(H)——由誘導器或風機盤管使室內(nèi)空氣循環(huán)；甚大風量(VH)——為達到室內(nèi)環(huán)境均勻而采用的超大風量循環(huán)。

表4中空氣循環(huán)類型和對應(yīng)對流換熱系數(shù)值摘自文獻[14]。本文實驗工況1、工況2和工況3對應(yīng)的壁面對流換熱系數(shù)如表4所示，實驗值是通過公式(8)利用實測獲得的空氣溫度與壁面溫度計算獲得。由表4可知對應(yīng)的空氣循環(huán)類型分別接近甚大風量、大風量和中等風量[14]，且隨風量減少，對流換熱系數(shù)顯著降低，從而說明其輻射負荷即輻射時間因子也將隨風量減少而降低。ASHRAE提供的傳遞函數(shù)數(shù)據(jù)庫中內(nèi)表面對流換熱系數(shù)壁面為4.68 W/(m2·K)，地板為4.37 W/(m2·K)[5]，可知對應(yīng)的空氣循環(huán)類型基本為大風量[14-15]。

3.2 輻射時間因子的修正

根據(jù)輻射得熱和輻射負荷之間的傳遞函數(shù)式(4)，構(gòu)造如式(7)的系數(shù)矩陣，通過矩陣運算可得相應(yīng)的輻射時間因子。輻射時間因子表示某時刻的輻射得熱轉(zhuǎn)化為當時刻及以后逐時刻輻射負荷的比例，第一項r1值即表示轉(zhuǎn)化為當時刻冷負荷的比例，其意義與傳遞函數(shù)系數(shù)ν0值相同。工況1、工況2和工況3中輻射時間因子的r1值均在0.5左右，第二項r2值已降到0.1以下，前五項值大于0.75，可知某時刻的冷負荷主要為前5個時刻的輻射得熱的貢獻。三個工況的輻射時間因子前五項值見圖8。

由圖8可知：三個工況的輻射時間因子均隨著時間的增加而減小，r1值與表3中ν0值也基本相同。同時，在1時刻，隨著送風風量減小，r1值有所減小，原因與傳遞函數(shù)ν0值變化規(guī)律相同。以工況2為基準，工況1和工況3的r1值的相對誤差絕對值在4%左右。在計算輻射時間因子時，假設(shè)區(qū)域是絕熱的，也就是所有表面的輻射得熱最終都轉(zhuǎn)化為冷負荷，即24項輻射時間因子之和為1[1]。因此，考慮將工況1的r1值減小4%，相應(yīng)將后面r2～r5各增加r1值減小量的25%(即r1值的1%)；同樣，將工況3的r1值增加4%，相應(yīng)將后面r2～r5各減小r1值增加量的25%(即r1值的1%)，調(diào)整后三個工況的輻射時間因子已基本相同。

在輻射時間序列方法中，采用輻射時間因子計算輻射得熱到輻射負荷的轉(zhuǎn)化過程。目前，房間輻射時間因子由ASHRAE提供的軟件計算得到，但該方法未考慮風量變化的影響。由表4可知，本文實驗中工況2基本對應(yīng)空氣循環(huán)類型為大風量，即ASHRAE提供的輻射時間因子軟件的計算條件。以上述調(diào)整方法的逆思路作為輻射時間因子隨風量變化的修正方法：以ASHRAE提供的軟件計算的輻射時間因子為基礎(chǔ)，如果房間送風為甚大風量，將計算的輻射時間因子的r1值增加4%，并對r2～r5各減小r1值增加量的25%(即r1值的1%)。如果房間送風為中等風量，將計算的輻射時間因子的r1值減小4%，并對r2～r5各增加r1值減小量的25%(即r1值的1%)。

4 結(jié)論

1)提出了一種通過實驗計算輻射時間因子的方法，采用對流輻射分離方法計算逐時輻射得熱和輻射負荷，利用最小二乘法得到輻射傳遞函數(shù)，再通過矩陣運算得到輻射時間因子，研究結(jié)果表明該方法合理可行。

2)隨著送風風量減小，傳遞函數(shù)系數(shù)ν0值略有減小，輻射時間因子第一項r1值也略有減小，表明某時刻的輻射得熱轉(zhuǎn)化為當時刻輻射負荷的比例有所下降。

3)目前房間輻射時間因子由ASHRAE提供的軟件計算得到，未能體現(xiàn)送風風量變化的影響。

本文研究了輻射時間因子隨送風風量改變的變化規(guī)律，并給出了修正方法，為輻射時間因子的準確應(yīng)用提供了依據(jù)。

本文受滬江基金項目(D14003)資助。(The project was supported by the Hujiang Foundation of China(No. D14003).)

[1] Spitler J D, Fisher D E, Pedersen C O. The radiant time series cooling load calculation procedure[J]. ASHRAE Transactions, 1997, 103(2): 503-515.

[2] ASHRAE. ASHRAE handbook-fundamentals[M]. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineer Inc, 2009.

[3] Spitler J D, Fisher D E. On the relationship between the radiant time series and transfer function methods for design cooling load calculations[J]. International Journal of HAVC & R Research, 1999, 5(2): 123-136.

[4] Mc Quiston F C, Parker J D, Spitler J D. Heating ventilating and air conditioning analysis and design[M]. 6th ed. San Francisco: John Wiley & Sons Inc, 2004: 231-279.

[5] I P SENG I U. Experiment validation of the radiant time series method for cooling load calculations[D]. Macau: University of Macau, 1999: 55-64.

[6] Spitler J D, Nigusse B A. Refinements and improvements to the radiant time series method[J]. ASHRAE Transactions, 2010, 116(2): 542-549.

[7] Sowell E F. Classification of 200640 parametric zones for cooling load calculation[J]. ASHRAE Transactions, 1988, 94(2): 754-777.

[8] Sowell E F. Cross-check and modification of the DOE-2 program for calculation of zone weighting factors[J]. ASHRAE Transactions, 1988, 94(2): 737-753.

[9] Sowell E F. Load calculation for 200640 zones[J]. ASHRAE Transactions, 1988, 94(2): 716-736.

[10] Calvin l. PRF/RFT Generator[EB/OL]. (2013-08-04) [2014-04-20]. http://www.hvac.okstate.edu/resources.html.

[11] 單寄平. 空調(diào)負荷實用計算法[M]. 北京: 中國建筑工業(yè)出版社, 1986: 78-88.

[12] 連之偉, 馬仁民, 李海林. 下送風室內(nèi)熱源條件下房間傳遞函數(shù)的實驗研究[J]. 西安冶金建筑學院學報, 1992, 24(2): 123-129. (Lian Zhiwei, Ma Renmin, Li Hailin. The experimental study of the room transfer functions formed by heat source inside the room with upward displacement[J]. J. Xi′an Inst. of Metall. & Cons. Eng., 1992, 24(2): 123-129.)

[13] 彥啟森, 趙慶珠. 建筑熱過程[M]. 北京：中國建筑工業(yè)出版社, 1986: 115-129.

[14] 陳沛霖, 曹叔維, 郭建雄. 空氣調(diào)節(jié)負荷計算理論與方法[M]. 上海: 同濟大學出版社, 1987: 166-182.

[15] 曹叔維. 房間熱過程和空調(diào)負荷[M]. 上海: 上海科學技術(shù)文獻出版社, 1991: 243-268.

[16] 章熙民, 任澤沛, 梅飛鳴. 傳熱學[M]. 北京: 中國建筑工業(yè)出版社, 2001: 218-255.

About the corresponding author

Huang Chen, female, professor, Ph. D., supervisor, School of Environment and Architecture, University of Shanghai for Science and Technology, +86 21-55273409, E-mail: hcyhyywj@163.com. Research fields: thermal environment in large space building, building energy conservation, cleaning air-conditioning.

Improvement to the Radiant Time Factors under Different Air Volume Based on Experiment

Lü Liugen Huang Chen Chen Jianchang Li Li Bai Tianyu

(School of Environment and Architecture, University of Shanghai for Science and Technology, Shanghai, 200093, China)

The radiant time factors are used to calculate radiant cooling load in the radiant time series method. Based on measurement data of hourly heat gain and cooling load collected within 24 hours, radiant heat gain and radiant cooling load can be calculated by separating convection with radiation. Transfer function is obtained and a 24 term radiant time factors can be calculated by transfer function coefficients matrix. By comparing the results with experiment, it shows that the calculation method of radiant time factors is feasible. The first value of radiant time factors decreases as supply air volume decreases and the sum of the first five factors is more than 0.75. For the radiant time factors calculated by ASHRAE software, the improvement method can be used to calculate radiant cooling load under different air circulation after correcting the first five factors.

radiant heat gain; radiant cooling load; radiant time factors; improvement method

0253- 4339(2015) 01- 0059- 06

10.3969/j.issn.0253- 4339.2015.01.059

國家自然科學基金(51278302、51108263)資助項目。(The project was supported by the National Natural Science Foundation of China (No. 51278302 & No. 51108263).)

2014年4月29日

TU113.3； TU831.2； O211.61

A

黃晨，女，教授，博士生導師，上海理工大學環(huán)境與建筑學院，(021)55273409，E-mail: hcyhyywj@163.com。研究方向：大空間建筑室內(nèi)熱環(huán)境，建筑節(jié)能，凈化空調(diào)技術(shù)。