设计参考暖通空调HV&AC2012年第42卷第5期
41
辐射地板供冷量及表面温度分布的简化计算方法*
清华大学张伦☆刘晓华高志宏江亿
摘要辐射地板的供冷量、表面温度分布的均匀性和表面最低温度是研究者以及工程人员重点关注的三个问题。通过求解导热方程,得出辐射地板热阻的解析表达式,提出一种辐射地板供冷量、表面温度分布的简化计算方法。在实际应用中,可通过公式便捷地计算出辐射地板供冷性能相关的主要参数。
关键词地板辐射供冷量表面温度简化方法
Simplifiedcalculationmethodforfhecoolingcapacityand
surfacelemperaturedistributionoffloorpaneI
ByZhongLU几★。LIuXiaohuo.GooZhihongondJiongyf
Abstracttemperature
Coolingcapacity,the
uniformityofsurfacetemperaturedistributionandthelowestsurface
engineersmainlyfocuson.Obtainstheheatresistance
a
are
threeaspectswhichresearchersand
analyrticalSOlutionbysolvingtheheatconductionequation,andputsforward
simplifiedmethodfor
calculatingthecoolingcapacityandsurfacetemperaturedistribution.Inactualapplicationof、floorpanel
cooling,the
keyparameters
call
becalculatedeasilybytheequations.
method
’
Keywords
floorpanel,coolingcapacity,surfacetemperature,simplified
★TsinghuaUniversity,Beijing,China
O引言Holopainen等人利用非均匀网格,模拟了辐射地板的导热[71;Koschenz等人建立并通过实验验证了辐射地板与空气换热的模型18]。3)针对实际系统,研究辐射地板的应用,其中Olesen等人研究了辐射地板应用的可行性和限制性[9],计算了辐射地板的换热系数及供冷量H妇;Song等人研究了辐射地板与置换通风相结合的系统[111;卢军等人通过实验研究了辐射地板供冷系统[12-13]。
工程应用中,在给定供回水温度和室内参数条件下,尚没有较简便的方法计算实际辐射地板的供冷量和表面温度分布。数值模拟方法由于计算量大,难以在工程中简化应用;采用实验测试结果的方法,对于不同结构材料的辐射地板的适用性较
☆张伦,男,1987年1月生,在读博士研究生
地板辐射供冷是一种利用建筑物内部地面进行供冷的系统,地板通过对流换热冷却周围空气,同时通过辐射与四周的围护结构进行换热。在辐射地板的实际应用中,辐射地板的供冷量、表面温度分布的均匀性和表面最低温度是研究者以及工程人员重点关注的3个问题。目前对于辐射地板的研究主要集中在以下几个方面:1)通过理论分析求解辐射地板的导热问题,其中Koschenz等人得出了单层辐射地板二维导热的解析解[13;DeMonte求解了一维多层固体的导热问题[z3;Lu等人求解了一维多层固体非稳态的导热问题[31;Beck等人研究了部分加热时的二维固体导热问题[41;Laouadi建立了辐射地板的半解析模型r53;Weber等人建立了辐射地板的RC网络模型[6]。2)通过数值模拟和实验分析辐射地板的性能,其中
*铁道部科技研究开发计划项目:铁路客站采暖空调技术综合
利用及关键技术研究(编号:2010G009一G)
’100084清华大学建筑技术科学系
(010)62773772
E-mail:zhanglun09@mails.tsinghua.edu.cn收稿日期:2010—07—05一次修回:2010—09—13二次修回:2012-04一05
42
暖通空调HV&AC靴年第旺誊第5期设计参考
差。本文通过对辐射地板导热过程的理论分析,得到辐射地板供冷量和表面最低、最高温度的表达式,可方便计算辐射地板的性能。1辐射地板供冷■计算
辐射地板的结构如图1所示,辐射地板的供冷量在不同的工作环境中差别较大。但辐射地板的热阻只与辐射地板的结构与材料相关.因此,通过求解导热方程,得到辐射地板热阻的表达式.再结合辐射地板工作环境的参数。例如室内空气温度、壁面温度、对流换热表面传热系数和辐射换热表面传热系数,即可求解出辐射地板的表面平均温度和供冷量。一般辐射地板在沿供回水管方向上温度变化比较缓慢,所以可以将三维的导热问题简化为二维的导热问题。在求解辐射地板供冷量时,核心的参数是辐射地板表面的平均温度,因此,可以将二维的导热问题通过热阻转化为一维问题求解。而在求解表面温度分布时,核心的参数是辐射地板表面的最高温度和最低温度。因此。必须考虑求解
二维的导热问题。
考刮,y刮lM
耋=吣=士专
(3)
(4)
距离,mlL为管同距,m}k为供回水平均温
式中q为辐射地板的供冷量,w/订。
通过求解导热方程,可以得到热阻的表达式:
R一2赤
‘7)
’
,
r=峙+哿+薹GC。s)r㈦
cc。=趾鲑蠹孑
B…i+2ns。lgdt一2声qq)一e铷I
Bi一2ns’。
u=陆1tTdt]1
(9)
Bi=华
(11)
式(8)~(11)中8为水管直径,m;Bi为毕渥数。
b自一十■
圈l辐射地摄的结构示膏田
1
1均一介质辐射地板的热阻
均一介质辐射地板(如图Ib所示)的导热方程
因此,从辐射地扳供回水管管壁到室内空气的热阻等于辐射地板供凰水管管壁到辐射地板表面
满足‘13如下关系;
蟊zt-t--ay2a'zt=。甜。够。
(1)…(2)
均一介质矩形固体的导热满足以下条件(其中
^考=讹叫・y=o
2042(5)
张伦.等:辐射地板供冷量度表面温度分布的筒化计算方法43
孚+器一o
扯2’劬2”
㈨,…7
(1s)(16)(17)
R.=R0+Rl+R2+Rd+÷
(23)
A・苗=h(top
t岫-0
如仅考虑辐射地板自身.那么辐射地板供回水管管壁到辐射地板表面的热阻为
R=Ro+RI+R:+Rd
t一,(z),,=Lm
(24)
箬:0.z:0,z:L∞
固体的热阻:
通过热阻,可以推导出辐射地板表面平均温度的表达式以及供冷量的表达式:
“
通过求解温度分布,可以得到均一介质的矩形
R:!!兰!12二!!!!墨12
‘
f,-址铲R=址PA警(25)
.^
+1
…7
口
q=三LiI=li±
^
^●
(26)
弘圳)如一舭础汕:拿(18)
^1
式中L。为均一介质厚度,m;舢为均一介质导热系数,W/(m・K)。
1
3多层介质辐射地板的热阻
如图2所示。组成辐射地扳的各个部分的热阻分别为
一
Rl=警
1.4辐射地板供冷量
通过模拟和实验两种方法验证辐射地板的热阻、表面平均温度和供冷量的表达式。模拟的方法可以得到辐射地板的温度分布,主要验证了辐射地
板每层介质的热阻;通过实验可以测得辐射地板表
面的温度和热流,主要验证了表面平均温度和供冷量的取值。
模拟和实验采用相同的参数设定口‘]:豆石混凝土层的导热系数为1.84W/(m・K)。水泥砂浆的导热系数为0.93w/(m・K),大理石地面的导热系数为3.83w/(m・K),PERT塑料管的导热系数为0.4W/(m・K)。实际辐射板结构如图3所示。对流和辐射混合抉热的表面传热系数h=
6.5~7w/(m2・K),其中辐射换热表面传热系数为5.5w/(m・K),冷表面在下的对流换热表面传
R。=丽L[-n去+争鲁+善掣]c・。,
(20)
R2=警
(21)
Rd一譬刍
(22’
热系数为1~1.5W/(m2・K)E”].对应的外界温度为28℃。图4为模拟的辐射地板内部温度分布
和表面温度分布。
式(21),(22)中工咕为另一种均一介质厚度,mI儿为另一种均一介质导热系数,W/(m・K);丸为供回水管管壁导热系数,W(m・K);Ld为供回水管
壁厚度。m。
■一
l^m十
鼍严
Lf
自_m
■■]■■
I
I
V
.
.
.1醐
昌}_■■J
盈2●甜地饭撼盥示l圈
懒鲵坳黝獬礁
田3蜜际辅射地拉耸掏
实验测量了6种工况的供回水温度、表面平均温度和热流,与用式(25),(26)计算得到的相同供回水温度下的表面平均温度和热流进行比较.详见
由热阻的叠加原理,得蓟对于多层介质的辐射地板,辐射地板供回水管管壁到辐射地板表面的热阻:
44
震
^内%-‘廑
目《^-
暖通空调HV&Ac2012年第让奄第5期设计参考
表1,误差均在10%以内。
在相同供回水温度下,通过模拟比较了辐射地板每一层的热阻、表面平均温度和热流。除供回水管管壁热阻外,其余误差均在lo%以内,见
表2。
2辐射地板裹面温度分布
在辐射地板的实际工程应用中.其表面温度分布的均匀性,以及表面的最低温度是研究者和工程应用人员关注的问题。辐射表面的最低温度必须高于辐射地板表面空气的露点温度,否则将会出现结露,因此辐射地板表面的最低温度成为辐射地板供冷量的主要限制。在此前提之下,得到辐射地板表面的最低温度就显得尤为重要。2.1均一介质辐射地板的表面温度分布
考虑到辐射地板表面的最高温度‰与最低温度£,之差一定小于辐射地板的供回水温度差。所以可以定义辐射地板的衰减系数s,s为一个小于1的常数,且s不随供回水温度变化而变化。
S:纽二立
£h—t。
(Z7)
b^■i■&
圈4横拟结果l供水■鹰15℃.四木汪麈18_e
裹1寰面平均温度、热流和热阻比较
供承温度l。/℃
回水温度thl℃
128217.63152320.635304
128213.0217.6415.33
—2.28—390
2仉7351.17
l&0217.6415.332n215n47
--2.51—137
13.20l7.6015.4520.7851.51
1320
17.6315
23
17.69ls.4520.2849.91
—241-3.11
供回水平均温度lm/℃衰面平均沮度t./'C热流∥(w/种)
20.1650
97
垫里堡垒堡:三!圣!!:!!!
实验
垒!!:
工况4
=!:!!
误差/%
垒!!!
实验
坠!!!
工况5解析
=!:些
谋差/%
垒!!!
实验
坠!!!
工况6
=!:::
误差/%
解析解析
通过求解均一介质的二维导热微分方程,可以得到衰减系数的表达式。
S:£!Q!盟=!f皇:盟:
t(0,0)一t(L,O)
(28)
20q2(5)
张伦,等:辐射地板供冷量及表面温度分布的简化计算方法
表2各层介质热阻比较
解析结果模拟结果误差/%
45
那么表面的最低温度和最高温度分别为
1
供水温度£。/℃回水温度th/℃
15.021.0
15.021.00.0470.02650.006
9.153.468.33
t。。一t。一去s(£h—t。)
厶1
(30)
豆石混凝土层热阻Ro/(m2・"c/w)0.051
t
水泥砂浆层热阻R1/(m2‘'C/W)0.027大理石层热阻Rz/(m2・C/w)0.006
供回水管管壁热阻Rj(mz.℃/w)0.012
s}1一t。+去S(“一£。)
‘
(31)
0.014—15.00
0
表3给出了不同工况下的衰减系数,表4给出了辐射地板表面温度分布的模拟值和由式(30)和(31)计算得到的表面最低温度与表面最高温度。
表3不同工况下衰减系数比较
供水温度/'C
O
水与管壁换热热阻R。/(m2・'c/w)0.001
总热阻R/(m2・℃/W)0.097
21.67表面平均温度i。/'C
0.0015—10.000.09421.68
3.30一O.31
垫迹吐!旦£垡2
.
:.
!!:!!i!:i!:二!:i!
分析式(28)中各项随m增加的变化趋势,在工程误差允许范围内,可将S简化为5。:
S。=
!
回水温度/'C
00
模拟S
0O
So
0
0
f土+11盟U]岛。1/2
婴查里堕!!
21.021.0
0.14940.1321
(29)
¨¨¨O
15.0
扒博MO
15.5
O
m4m3m7
警
l
0.1496
.
堡查望壅!兰
模拟解析
15.015.0
室!壅亘婆室坌查鲨茎室堕墨塑兰室亘望薹!!室亘兰塑塑鏖!!
0.900.79
21.1421.25
查亘量壹塑壅!!耋堕墨堡塑星!!
21.5821.65
20.6820.85
堡薹!兰
2-2
实际的辐射地板往往由多层介质构成,同时必须要考虑供回水管对表面温度分布的影响。
7嚣雾鬻巽誓囊爸夏雨物甫同时
.=【7一丛笔妄%≥蔓鲁警(33)
L■_十乞j十h。,
u纠
吼一丽Ate∞qh=瓦而Ath/x/kth
k十
!:!!!:!!.!:墼
出c1_出“
q(34)
L53’
&c—os—oc,础n2ts一“
对于多层的辐射地板,可以得到等效的导热系‘数:
然后通过管壁的热阻计算管壁上的温差变化:
&cw—g。Z“-gL_.d,Athw=qh等L_d
最后可以得到管壁对于供回水温差的减少值:
△fw一△£cw一△£hw
(36)
(37)
如果考虑供回水管的影响,供回水管的热容可以忽略,但其热阻不可忽略,因此必须计算供回水管管壁上的温差变化。在工程计算中,首先可以分别得到供水管与回水管的供冷量:
;
通过式(32)~(37)的推导,可以分别计算在未考虑管壁影响和考虑管壁影响下,多层介质辐射地板的温度分布,表5,6分别列出了两种工况的计算
绡米。
.麦!塞亘墨室坌查丝篁!垡坐里堕垄!!翌!旦垄墨壁垫!!旦2
表面平均温度/
供冷量/
53.89j2.6750.2249.67一1.10
供回水温差/
1110.9
衰减系数s
0.10700.13480.08270.1348
表面温差/
0.110.140.080.12
表面最高温度/
19.8720.0320.2820.360.39
表面最低温度/
19.7619.8920.1920.240.25
未考虑管壁修正模拟
解析
19.8219.9620.2420.30
考虑管壁修正
模拟解析
误差/%0.30
表6表面温度分布比较(供水温度为15℃。回水温度为21℃)
表面平均温度/
供冷量/供回水温差/衰减系数s
0.10710.13480.08180.1348
表面温差/7表面最高温度/・表面最低温度/
!
未考虑管壁修正模拟
解析
考虑管壁修正
模拟解析误差/%
21.2921.4021.6021.660.28
!!!翌竺
41.0740.1138.4437.84—1_56
旦
6665.3
里
0.540.800.490.71
里
21.6021.802l‘8422.020.82
竺
20.9621.0021.3521.31一O.19
46
暖通空调HV&AC2012年第42卷第5期设计参考
3结论
cooling
431
modelforbuildingenergysimulationsoftware
andEnvironment,2004,39(4):421—
3.1辐射地板热阻几乎不随辐射地板工作的环境变化,其解析表达式可参见式(23)和(24)。在已知辐射地板供回水温度的情况下,辐射地板的表面平均温度和供冷量分别由式(25)和(26)计算得到,比较计算值与实验值,误差在10%以内。误差产生的原因主要是接触热阻的影响和材料参数的误差。3.2辐射地板表面温度分布的均匀性问题以及表面最低温度,可通过式(30)和(31)得到。实际辐射地板必须考虑多层介质和供回水管管壁的影响,简化的计算方法通过式(32)~(37)表述。
3.3采用本文提出的方法,可在工程应用中简便计算和辐射地板相关的重要参数,而不需要数值计算传热方程,大幅度缩短了计算时间,同时保证了计算的精度。
参考文献:[1]Koscheaz
M,Lehmann
B
ThermoaktiveZurich:
[J].Building
f6]WeberT,J6hannessonG.AnoptimizedRC-network
forthermallyBuildingand
activated
building
components[J].
Environment,2005,40(1):1—14
R,Tuomaala
P,Piippo
J.Uneven
[7]Holopainen
griddingofthermalnodalnetworksinfloorheating
simulations[J].Energy
(10):1107—1114
andBuildings,2007,39
[8]Koschenz
M,Dorerv.Interactionof
core
an
air
system
and
withconcrete
conditioning[J].Energy
Buildings,1999,30(2):139—145
r9]Olesen
B
w.Possibilitiesandlimitationsofradiant
floorcooling[G3∥ASHRAETrans,1997,103(1):
42-48
[10]Olesen
B
W,Michel
E.Heatexchangecoefficient
betweenfloorsurfaceandspacebyfloorcoolingtheory
or
a
Bauteilsysteme
tabs[M].
EMPA
question
of
definition[G]∥ASHRAE
S,et
a1.Performance
EnergiesystemeHaustechnik,2000
Trans,2000,106(1):684—694
in
a
[2]De
MonteF.Transientheatconduction
one-
[11]song
with
D,Kim
a
T,Song
dimensionalcomposite
slab.A‘natural’analytic
JournalofHeatandMass
evaluationof
radiantfloorcoolingsystemintegrated
approach[J].International
dehumidifiedventilation口].AppliedThermal
Transfer,2000,43(19):3607—3619
Engineering,2008,28(11/12):1299—1311
inthe
[3]Lux,TervolaP.Transientheat
composite
slatranalyrical
conduction
[12]卢军,陈金花,高殿策,等。住宅用地板辐射供冷系统
method[J].Journalof
的实验分析口].暖通空调,2007,37(12):20—24[13]张瑜,胡松涛,王东.冷辐射地板系统的实验研究[J].
能源技术,2004,25(2):64—68
PhysicsA:MathematicalandGeneral,2005,38(1):
81-96
[4]BeckJV,Hail-SheikhA,AmosD
E,eta1.
[14]魏立峰,戴立生,杨强,等.基于地板辐射盘管试块散
冷测试的研究[c]∥全国暖通空调制冷2010学术年
会论文集,2010
Verificationsolutionforpartialheatingofrectangular
solids[J].InternationalJOurnal
[5j
of
ofHeatandMass
Transfer,2004,47(19/20):4243—4255
Laouadi八Development
a
[15]章熙民,任泽霈,梅飞鸣.传热学[M].4版.北京:中
heating
radiant
and
国建筑工业出版社,2001
(上接第40页)
参考文献:
[13刘大龙,刘加平,杨柳,等.建筑气候区域性研究口].
暖通空调,2009,39(5):93—96
[5]陆耀庆.暖通空调设计指南[Ml北京:中国建筑工
业出版社,1996:5—10
[6]黄季宜,韩维平.昆明新机场航站楼热工性能优化及
空调系统新风利用分析[J].暖通空调,2009,39(5):
7—10
7
・
[2]付祥钊,陈敏.关于民用建筑通风室内外计算参数的
思考[J].暖通空调,2011,41(5):26—31
[3]王伟,罗会龙,高亚锋,等.昆明地区夏季居住建筑室
内热环境测试及热舒适调查[J].暖通空调,2011,41
(6):93—95
[7]王昕,黄晨,曹伟武.t自然通风在大空间建筑空调系
统下室内热环境中应用的理论分析及实测研究EJ].
暖通空调,2009,39(5):62—66
[4]陆耀庆.实用供热空调设计手册[M].2版.北京:中
国建筑工业出版社,2008:1479—1481
IS]刘晓华,江亿,等.温湿度独立控制空调系统[M].北
京:中国建筑工业出版社,2006
辐射地板供冷量及表面温度分布的简化计算方法
作者:作者单位:刊名:英文刊名:年,卷(期):
张伦, 刘晓华, 高志宏, 江亿, Zhang Lun, Liu Xiaohua, Gao Zhihong, Jiang Yi清华大学 100084
暖通空调
Heating Ventilating & Air Conditioning2012,42(5)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_ntkt201205007.aspx
设计参考暖通空调HV&AC2012年第42卷第5期
41
辐射地板供冷量及表面温度分布的简化计算方法*
清华大学张伦☆刘晓华高志宏江亿
摘要辐射地板的供冷量、表面温度分布的均匀性和表面最低温度是研究者以及工程人员重点关注的三个问题。通过求解导热方程,得出辐射地板热阻的解析表达式,提出一种辐射地板供冷量、表面温度分布的简化计算方法。在实际应用中,可通过公式便捷地计算出辐射地板供冷性能相关的主要参数。
关键词地板辐射供冷量表面温度简化方法
Simplifiedcalculationmethodforfhecoolingcapacityand
surfacelemperaturedistributionoffloorpaneI
ByZhongLU几★。LIuXiaohuo.GooZhihongondJiongyf
Abstracttemperature
Coolingcapacity,the
uniformityofsurfacetemperaturedistributionandthelowestsurface
engineersmainlyfocuson.Obtainstheheatresistance
a
are
threeaspectswhichresearchersand
analyrticalSOlutionbysolvingtheheatconductionequation,andputsforward
simplifiedmethodfor
calculatingthecoolingcapacityandsurfacetemperaturedistribution.Inactualapplicationof、floorpanel
cooling,the
keyparameters
call
becalculatedeasilybytheequations.
method
’
Keywords
floorpanel,coolingcapacity,surfacetemperature,simplified
★TsinghuaUniversity,Beijing,China
O引言Holopainen等人利用非均匀网格,模拟了辐射地板的导热[71;Koschenz等人建立并通过实验验证了辐射地板与空气换热的模型18]。3)针对实际系统,研究辐射地板的应用,其中Olesen等人研究了辐射地板应用的可行性和限制性[9],计算了辐射地板的换热系数及供冷量H妇;Song等人研究了辐射地板与置换通风相结合的系统[111;卢军等人通过实验研究了辐射地板供冷系统[12-13]。
工程应用中,在给定供回水温度和室内参数条件下,尚没有较简便的方法计算实际辐射地板的供冷量和表面温度分布。数值模拟方法由于计算量大,难以在工程中简化应用;采用实验测试结果的方法,对于不同结构材料的辐射地板的适用性较
☆张伦,男,1987年1月生,在读博士研究生
地板辐射供冷是一种利用建筑物内部地面进行供冷的系统,地板通过对流换热冷却周围空气,同时通过辐射与四周的围护结构进行换热。在辐射地板的实际应用中,辐射地板的供冷量、表面温度分布的均匀性和表面最低温度是研究者以及工程人员重点关注的3个问题。目前对于辐射地板的研究主要集中在以下几个方面:1)通过理论分析求解辐射地板的导热问题,其中Koschenz等人得出了单层辐射地板二维导热的解析解[13;DeMonte求解了一维多层固体的导热问题[z3;Lu等人求解了一维多层固体非稳态的导热问题[31;Beck等人研究了部分加热时的二维固体导热问题[41;Laouadi建立了辐射地板的半解析模型r53;Weber等人建立了辐射地板的RC网络模型[6]。2)通过数值模拟和实验分析辐射地板的性能,其中
*铁道部科技研究开发计划项目:铁路客站采暖空调技术综合
利用及关键技术研究(编号:2010G009一G)
’100084清华大学建筑技术科学系
(010)62773772
E-mail:zhanglun09@mails.tsinghua.edu.cn收稿日期:2010—07—05一次修回:2010—09—13二次修回:2012-04一05
42
暖通空调HV&AC靴年第旺誊第5期设计参考
差。本文通过对辐射地板导热过程的理论分析,得到辐射地板供冷量和表面最低、最高温度的表达式,可方便计算辐射地板的性能。1辐射地板供冷■计算
辐射地板的结构如图1所示,辐射地板的供冷量在不同的工作环境中差别较大。但辐射地板的热阻只与辐射地板的结构与材料相关.因此,通过求解导热方程,得到辐射地板热阻的表达式.再结合辐射地板工作环境的参数。例如室内空气温度、壁面温度、对流换热表面传热系数和辐射换热表面传热系数,即可求解出辐射地板的表面平均温度和供冷量。一般辐射地板在沿供回水管方向上温度变化比较缓慢,所以可以将三维的导热问题简化为二维的导热问题。在求解辐射地板供冷量时,核心的参数是辐射地板表面的平均温度,因此,可以将二维的导热问题通过热阻转化为一维问题求解。而在求解表面温度分布时,核心的参数是辐射地板表面的最高温度和最低温度。因此。必须考虑求解
二维的导热问题。
考刮,y刮lM
耋=吣=士专
(3)
(4)
距离,mlL为管同距,m}k为供回水平均温
式中q为辐射地板的供冷量,w/订。
通过求解导热方程,可以得到热阻的表达式:
R一2赤
‘7)
’
,
r=峙+哿+薹GC。s)r㈦
cc。=趾鲑蠹孑
B…i+2ns。lgdt一2声qq)一e铷I
Bi一2ns’。
u=陆1tTdt]1
(9)
Bi=华
(11)
式(8)~(11)中8为水管直径,m;Bi为毕渥数。
b自一十■
圈l辐射地摄的结构示膏田
1
1均一介质辐射地板的热阻
均一介质辐射地板(如图Ib所示)的导热方程
因此,从辐射地扳供回水管管壁到室内空气的热阻等于辐射地板供凰水管管壁到辐射地板表面
满足‘13如下关系;
蟊zt-t--ay2a'zt=。甜。够。
(1)…(2)
均一介质矩形固体的导热满足以下条件(其中
^考=讹叫・y=o
2042(5)
张伦.等:辐射地板供冷量度表面温度分布的筒化计算方法43
孚+器一o
扯2’劬2”
㈨,…7
(1s)(16)(17)
R.=R0+Rl+R2+Rd+÷
(23)
A・苗=h(top
t岫-0
如仅考虑辐射地板自身.那么辐射地板供回水管管壁到辐射地板表面的热阻为
R=Ro+RI+R:+Rd
t一,(z),,=Lm
(24)
箬:0.z:0,z:L∞
固体的热阻:
通过热阻,可以推导出辐射地板表面平均温度的表达式以及供冷量的表达式:
“
通过求解温度分布,可以得到均一介质的矩形
R:!!兰!12二!!!!墨12
‘
f,-址铲R=址PA警(25)
.^
+1
…7
口
q=三LiI=li±
^
^●
(26)
弘圳)如一舭础汕:拿(18)
^1
式中L。为均一介质厚度,m;舢为均一介质导热系数,W/(m・K)。
1
3多层介质辐射地板的热阻
如图2所示。组成辐射地扳的各个部分的热阻分别为
一
Rl=警
1.4辐射地板供冷量
通过模拟和实验两种方法验证辐射地板的热阻、表面平均温度和供冷量的表达式。模拟的方法可以得到辐射地板的温度分布,主要验证了辐射地
板每层介质的热阻;通过实验可以测得辐射地板表
面的温度和热流,主要验证了表面平均温度和供冷量的取值。
模拟和实验采用相同的参数设定口‘]:豆石混凝土层的导热系数为1.84W/(m・K)。水泥砂浆的导热系数为0.93w/(m・K),大理石地面的导热系数为3.83w/(m・K),PERT塑料管的导热系数为0.4W/(m・K)。实际辐射板结构如图3所示。对流和辐射混合抉热的表面传热系数h=
6.5~7w/(m2・K),其中辐射换热表面传热系数为5.5w/(m・K),冷表面在下的对流换热表面传
R。=丽L[-n去+争鲁+善掣]c・。,
(20)
R2=警
(21)
Rd一譬刍
(22’
热系数为1~1.5W/(m2・K)E”].对应的外界温度为28℃。图4为模拟的辐射地板内部温度分布
和表面温度分布。
式(21),(22)中工咕为另一种均一介质厚度,mI儿为另一种均一介质导热系数,W/(m・K);丸为供回水管管壁导热系数,W(m・K);Ld为供回水管
壁厚度。m。
■一
l^m十
鼍严
Lf
自_m
■■]■■
I
I
V
.
.
.1醐
昌}_■■J
盈2●甜地饭撼盥示l圈
懒鲵坳黝獬礁
田3蜜际辅射地拉耸掏
实验测量了6种工况的供回水温度、表面平均温度和热流,与用式(25),(26)计算得到的相同供回水温度下的表面平均温度和热流进行比较.详见
由热阻的叠加原理,得蓟对于多层介质的辐射地板,辐射地板供回水管管壁到辐射地板表面的热阻:
44
震
^内%-‘廑
目《^-
暖通空调HV&Ac2012年第让奄第5期设计参考
表1,误差均在10%以内。
在相同供回水温度下,通过模拟比较了辐射地板每一层的热阻、表面平均温度和热流。除供回水管管壁热阻外,其余误差均在lo%以内,见
表2。
2辐射地板裹面温度分布
在辐射地板的实际工程应用中.其表面温度分布的均匀性,以及表面的最低温度是研究者和工程应用人员关注的问题。辐射表面的最低温度必须高于辐射地板表面空气的露点温度,否则将会出现结露,因此辐射地板表面的最低温度成为辐射地板供冷量的主要限制。在此前提之下,得到辐射地板表面的最低温度就显得尤为重要。2.1均一介质辐射地板的表面温度分布
考虑到辐射地板表面的最高温度‰与最低温度£,之差一定小于辐射地板的供回水温度差。所以可以定义辐射地板的衰减系数s,s为一个小于1的常数,且s不随供回水温度变化而变化。
S:纽二立
£h—t。
(Z7)
b^■i■&
圈4横拟结果l供水■鹰15℃.四木汪麈18_e
裹1寰面平均温度、热流和热阻比较
供承温度l。/℃
回水温度thl℃
128217.63152320.635304
128213.0217.6415.33
—2.28—390
2仉7351.17
l&0217.6415.332n215n47
--2.51—137
13.20l7.6015.4520.7851.51
1320
17.6315
23
17.69ls.4520.2849.91
—241-3.11
供回水平均温度lm/℃衰面平均沮度t./'C热流∥(w/种)
20.1650
97
垫里堡垒堡:三!圣!!:!!!
实验
垒!!:
工况4
=!:!!
误差/%
垒!!!
实验
坠!!!
工况5解析
=!:些
谋差/%
垒!!!
实验
坠!!!
工况6
=!:::
误差/%
解析解析
通过求解均一介质的二维导热微分方程,可以得到衰减系数的表达式。
S:£!Q!盟=!f皇:盟:
t(0,0)一t(L,O)
(28)
20q2(5)
张伦,等:辐射地板供冷量及表面温度分布的简化计算方法
表2各层介质热阻比较
解析结果模拟结果误差/%
45
那么表面的最低温度和最高温度分别为
1
供水温度£。/℃回水温度th/℃
15.021.0
15.021.00.0470.02650.006
9.153.468.33
t。。一t。一去s(£h—t。)
厶1
(30)
豆石混凝土层热阻Ro/(m2・"c/w)0.051
t
水泥砂浆层热阻R1/(m2‘'C/W)0.027大理石层热阻Rz/(m2・C/w)0.006
供回水管管壁热阻Rj(mz.℃/w)0.012
s}1一t。+去S(“一£。)
‘
(31)
0.014—15.00
0
表3给出了不同工况下的衰减系数,表4给出了辐射地板表面温度分布的模拟值和由式(30)和(31)计算得到的表面最低温度与表面最高温度。
表3不同工况下衰减系数比较
供水温度/'C
O
水与管壁换热热阻R。/(m2・'c/w)0.001
总热阻R/(m2・℃/W)0.097
21.67表面平均温度i。/'C
0.0015—10.000.09421.68
3.30一O.31
垫迹吐!旦£垡2
.
:.
!!:!!i!:i!:二!:i!
分析式(28)中各项随m增加的变化趋势,在工程误差允许范围内,可将S简化为5。:
S。=
!
回水温度/'C
00
模拟S
0O
So
0
0
f土+11盟U]岛。1/2
婴查里堕!!
21.021.0
0.14940.1321
(29)
¨¨¨O
15.0
扒博MO
15.5
O
m4m3m7
警
l
0.1496
.
堡查望壅!兰
模拟解析
15.015.0
室!壅亘婆室坌查鲨茎室堕墨塑兰室亘望薹!!室亘兰塑塑鏖!!
0.900.79
21.1421.25
查亘量壹塑壅!!耋堕墨堡塑星!!
21.5821.65
20.6820.85
堡薹!兰
2-2
实际的辐射地板往往由多层介质构成,同时必须要考虑供回水管对表面温度分布的影响。
7嚣雾鬻巽誓囊爸夏雨物甫同时
.=【7一丛笔妄%≥蔓鲁警(33)
L■_十乞j十h。,
u纠
吼一丽Ate∞qh=瓦而Ath/x/kth
k十
!:!!!:!!.!:墼
出c1_出“
q(34)
L53’
&c—os—oc,础n2ts一“
对于多层的辐射地板,可以得到等效的导热系‘数:
然后通过管壁的热阻计算管壁上的温差变化:
&cw—g。Z“-gL_.d,Athw=qh等L_d
最后可以得到管壁对于供回水温差的减少值:
△fw一△£cw一△£hw
(36)
(37)
如果考虑供回水管的影响,供回水管的热容可以忽略,但其热阻不可忽略,因此必须计算供回水管管壁上的温差变化。在工程计算中,首先可以分别得到供水管与回水管的供冷量:
;
通过式(32)~(37)的推导,可以分别计算在未考虑管壁影响和考虑管壁影响下,多层介质辐射地板的温度分布,表5,6分别列出了两种工况的计算
绡米。
.麦!塞亘墨室坌查丝篁!垡坐里堕垄!!翌!旦垄墨壁垫!!旦2
表面平均温度/
供冷量/
53.89j2.6750.2249.67一1.10
供回水温差/
1110.9
衰减系数s
0.10700.13480.08270.1348
表面温差/
0.110.140.080.12
表面最高温度/
19.8720.0320.2820.360.39
表面最低温度/
19.7619.8920.1920.240.25
未考虑管壁修正模拟
解析
19.8219.9620.2420.30
考虑管壁修正
模拟解析
误差/%0.30
表6表面温度分布比较(供水温度为15℃。回水温度为21℃)
表面平均温度/
供冷量/供回水温差/衰减系数s
0.10710.13480.08180.1348
表面温差/7表面最高温度/・表面最低温度/
!
未考虑管壁修正模拟
解析
考虑管壁修正
模拟解析误差/%
21.2921.4021.6021.660.28
!!!翌竺
41.0740.1138.4437.84—1_56
旦
6665.3
里
0.540.800.490.71
里
21.6021.802l‘8422.020.82
竺
20.9621.0021.3521.31一O.19
46
暖通空调HV&AC2012年第42卷第5期设计参考
3结论
cooling
431
modelforbuildingenergysimulationsoftware
andEnvironment,2004,39(4):421—
3.1辐射地板热阻几乎不随辐射地板工作的环境变化,其解析表达式可参见式(23)和(24)。在已知辐射地板供回水温度的情况下,辐射地板的表面平均温度和供冷量分别由式(25)和(26)计算得到,比较计算值与实验值,误差在10%以内。误差产生的原因主要是接触热阻的影响和材料参数的误差。3.2辐射地板表面温度分布的均匀性问题以及表面最低温度,可通过式(30)和(31)得到。实际辐射地板必须考虑多层介质和供回水管管壁的影响,简化的计算方法通过式(32)~(37)表述。
3.3采用本文提出的方法,可在工程应用中简便计算和辐射地板相关的重要参数,而不需要数值计算传热方程,大幅度缩短了计算时间,同时保证了计算的精度。
参考文献:[1]Koscheaz
M,Lehmann
B
ThermoaktiveZurich:
[J].Building
f6]WeberT,J6hannessonG.AnoptimizedRC-network
forthermallyBuildingand
activated
building
components[J].
Environment,2005,40(1):1—14
R,Tuomaala
P,Piippo
J.Uneven
[7]Holopainen
griddingofthermalnodalnetworksinfloorheating
simulations[J].Energy
(10):1107—1114
andBuildings,2007,39
[8]Koschenz
M,Dorerv.Interactionof
core
an
air
system
and
withconcrete
conditioning[J].Energy
Buildings,1999,30(2):139—145
r9]Olesen
B
w.Possibilitiesandlimitationsofradiant
floorcooling[G3∥ASHRAETrans,1997,103(1):
42-48
[10]Olesen
B
W,Michel
E.Heatexchangecoefficient
betweenfloorsurfaceandspacebyfloorcoolingtheory
or
a
Bauteilsysteme
tabs[M].
EMPA
question
of
definition[G]∥ASHRAE
S,et
a1.Performance
EnergiesystemeHaustechnik,2000
Trans,2000,106(1):684—694
in
a
[2]De
MonteF.Transientheatconduction
one-
[11]song
with
D,Kim
a
T,Song
dimensionalcomposite
slab.A‘natural’analytic
JournalofHeatandMass
evaluationof
radiantfloorcoolingsystemintegrated
approach[J].International
dehumidifiedventilation口].AppliedThermal
Transfer,2000,43(19):3607—3619
Engineering,2008,28(11/12):1299—1311
inthe
[3]Lux,TervolaP.Transientheat
composite
slatranalyrical
conduction
[12]卢军,陈金花,高殿策,等。住宅用地板辐射供冷系统
method[J].Journalof
的实验分析口].暖通空调,2007,37(12):20—24[13]张瑜,胡松涛,王东.冷辐射地板系统的实验研究[J].
能源技术,2004,25(2):64—68
PhysicsA:MathematicalandGeneral,2005,38(1):
81-96
[4]BeckJV,Hail-SheikhA,AmosD
E,eta1.
[14]魏立峰,戴立生,杨强,等.基于地板辐射盘管试块散
冷测试的研究[c]∥全国暖通空调制冷2010学术年
会论文集,2010
Verificationsolutionforpartialheatingofrectangular
solids[J].InternationalJOurnal
[5j
of
ofHeatandMass
Transfer,2004,47(19/20):4243—4255
Laouadi八Development
a
[15]章熙民,任泽霈,梅飞鸣.传热学[M].4版.北京:中
heating
radiant
and
国建筑工业出版社,2001
(上接第40页)
参考文献:
[13刘大龙,刘加平,杨柳,等.建筑气候区域性研究口].
暖通空调,2009,39(5):93—96
[5]陆耀庆.暖通空调设计指南[Ml北京:中国建筑工
业出版社,1996:5—10
[6]黄季宜,韩维平.昆明新机场航站楼热工性能优化及
空调系统新风利用分析[J].暖通空调,2009,39(5):
7—10
7
・
[2]付祥钊,陈敏.关于民用建筑通风室内外计算参数的
思考[J].暖通空调,2011,41(5):26—31
[3]王伟,罗会龙,高亚锋,等.昆明地区夏季居住建筑室
内热环境测试及热舒适调查[J].暖通空调,2011,41
(6):93—95
[7]王昕,黄晨,曹伟武.t自然通风在大空间建筑空调系
统下室内热环境中应用的理论分析及实测研究EJ].
暖通空调,2009,39(5):62—66
[4]陆耀庆.实用供热空调设计手册[M].2版.北京:中
国建筑工业出版社,2008:1479—1481
IS]刘晓华,江亿,等.温湿度独立控制空调系统[M].北
京:中国建筑工业出版社,2006
辐射地板供冷量及表面温度分布的简化计算方法
作者:作者单位:刊名:英文刊名:年,卷(期):
张伦, 刘晓华, 高志宏, 江亿, Zhang Lun, Liu Xiaohua, Gao Zhihong, Jiang Yi清华大学 100084
暖通空调
Heating Ventilating & Air Conditioning2012,42(5)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_ntkt201205007.aspx