《武汉工程大学学报》  2010年09期 31-34   出版日期:2010-09-30   ISSN:1674-2869   CN:42-1779/TQ
沥青混凝土路面融雪热流密度的分析


0引言道路交通事故在天气恶劣条件下会增加.据资料统计,2008年11月全国有19省(区、市)雨、雪、雾等恶劣天气交通事故导致死亡人数同比上升,山东省2009年11月阴、雨、雾、雪等恶劣天气条件下共发生道路交通事故、死亡人数分别占总数的16.8%、19.5%.在冬季下雪时,路面积雪、结冰,导致车辆行驶困难,极易发生交通事故.据统计,路面湿润时发生的事故是干燥路面的2倍,降雪时是干燥路面的5倍,结冰时是干燥路面的8倍[1].2010年1月9日晚22时许,受路面结冰影响,汉十高速公路湖北枣阳段张凹水库大桥处,发生4起道路交通事故,12台车辆相继连环追尾相撞,造成7人死亡、14人轻伤.可见道路融雪化冰对冬季道路交通条件的改善起到关键作用.目前国内外道路融雪化冰新技术是采用热力学法,包括导电混凝土、发热电缆及地源热泵.其中地源热泵融雪化冰技术利用太阳能蓄热,可节约能源[2],同时环保、高效,发展前景巨大.美国、日本、北欧等国家已经有一批典型的地源热泵道路、桥梁融雪示范工程,在国内已有学者开展了相关研究,胡文举[3]进行了沥青混凝土桥面融雪系统的基础研究,王华军[4]对流体加热水泥混凝土道路融雪传热介质特性进行了研究,但目前还没有对沥青混凝土路面融雪化冰的研究.本文以武汉地区2010年1月5日的气象条件为依据,分析京珠高速湖北段沥青混凝土路面融雪随机传热机理,并确定武汉地区沥青混凝土路面融雪热流密度.1分析条件1.1气象条件分析计算采用武汉2010年1月5日历史降雪记录.降雪时段是2:00 PM~10:00 PM,持续8 h;积雪厚度为12.8 mm,瞬态降雪强度为1.6 mm/h;气温-3~2 ℃;偏北风3~4级;相对湿度60%~95%.
1.2路面条件[5]湖北省京珠高速公路全部采用沥青路面,结构如下:a.4 cm SUPERPAVE12.5或AK16A抗滑表层;b.中面层为6 cm厚的AC201型沥青混凝土;c.底层为6 cm厚的AC20S沥青混凝土.1.3材料热物性参数由于摊铺沥青混合料温度不低于110~130 ℃,埋管采用耐高温高密度聚乙烯管,导热系数λp=0.39 W/(m·K),内径D1=22 mm,外径D2=25 mm.沥青混凝土λA=1.0 W/(m·K).2路面随机传热分析融雪化冰随机传热机理如图1所示.在沥青混凝土面层中埋耐高温高密度聚乙烯管,管内流体温度大于路面温度,于是沥青混凝土内部存在温差,热量由埋管传至路面(埋管下表面做绝热处理).路面上的积雪吸收热量,温度升高直至融化成水并蒸发,路面无积雪及水.这一过程是随机非稳态的,融雪所需的热量主要包括:升高雪的温度所需显热qs、融化雪所需潜热qm、对流换热qc、长波辐射换热qr、水蒸发所需热量qe;热源为地源热泵提供热量q0及吸收太阳辐射热qI.依据能量平衡,有q0+qI=qs+qm+Ar(qc+qr+qe)(1)图1融雪化冰随机传热机理
Fig.1Random heat transfer mechanism
based on snow melting则单位面积融雪化冰所需热流密度为q0=qs+qm+Ar(qc+qr+qe)(2)qs=ρsnowS(CwT0-CsnowTa(3)qm=ρsnowShf(4)qc=hc(Tf-Ta)(5)qr=εσ(T4f-T4sky)(6)qe=hmHfρw(100-φ)(7)q1=αI(8)
第9期屠艳平,等:沥青混凝土路面融雪热流密度的分析
武汉工程大学学报第32卷
其中:Ar为融雪等级;ρsnow为冰雪的密度,取917 kg/m3;S为降雪量,mm/h;Cw为水的比热,Csnow为冰雪的比热,按文献[4]中公式计算为2.05 kJ/(kg·K);T0为雪融化成水的温度,℃;Ta为空气温度,℃;hf为冰雪的融解热,取334 kJ/(kg);hc为对流换热系数,W/(m2·K),hc=5.678a+b(v/0.304)n,a=1.09,b=0.23,n=1;Tf为水膜温度,℃;ε为路面辐射率,沥青混凝土取0.96;σ取5.67×10-8 W/(m2·K4);Tsky为天空辐射温度,Tsky=0.522 T1.5db,Tdb为空气干球温度;hm为水蒸发时的传质速率,与风速有关;Hf为蒸发潜热,2 502 kJ/(kg);φ为空气相对湿度;α为太阳能集热效率,取0.6;I为大地表面太阳辐射量,W/m2.根据上述气象条件,可分别求出qs、qm、qc、qr、qe、qI,计算出q0.T0=0 ℃,Ta=-3 ℃,Tf=0 ℃=273 K,Tsky=2.71 ℃=270.29 K,qs=2.50 W/m2,qm=136 W/m2,hc=23.37 W/(m2·K),qc=70.11 W/m2,qr=5.99 W/m2,qe=40.3 W/m2 [6]武汉地区全年累计太阳辐射量为4 497 MJ/m2,I=142.60 W/m2,qI=85.56 W/m2,Ar取1,则q0=169.34 W/m2由于武汉地区冬天最低气温为-5 ℃,与北方相比较,气温较高,升高雪温所需显热少,同时辐射散热也很小,起主要作用的是融雪潜热、对流换热和蒸发热.融雪潜热与降雪量成线性关系,对流换热主要与风速和空气温度有关,蒸发热主要与风速有关.以上求出的是给定气象条件下的最大热负荷,同时考虑利用太阳能以节约资源,提高经济效益.3埋管传热分析埋管埋入沥青混凝土中,首先埋管内流体与管壁进行对流换热,然后管壁内传热,最后与沥青混凝土进行热交换(假设热量不向下传递).引入热阻的概念,管内平均温度为tw的乙二醇水溶液(质量分数为16%,流量为12.5 L/s),以定热流q给沥青混凝土进行传热.tw=tb=qRB(9)tb-t0=qRA(10)tw-t0=q(RB+RA)(11)其中tw、tb、t0分别为管内流体平均温度、孔壁温度和路面温度,℃;RB为埋管热阻,RA为埋管上沥青混凝土面层热阻,分别按式(12)、(16)计算.RB=Rp+Rf(12)其中Rp为管壁热阻、Rf管内流体换热热阻,分别按式(13)、(14)计算.Rp=ln(D2/D1)/(2πλp)(13)Rp=0.052 m2·K/W
Rf=1/(πD1hf)(14)其中hf为埋管内流体对流换热系数,W/(m2·K)hf=0.023 Re0.8P0.4rλf/D1(15)其中雷诺数Re=10 724,Pr为流体普朗特数=57.65,λf为管内流体导热系数,取0.558 W/(m·K),则hf=194 W/(m2·K)Rf=0.075 m2·K/W
RA=δ/λA(16)δ为埋管深度,取40 mm.RA=0.04 m2·K/Wtw-t0=q×(0.052+0.075+0.04),q取q0=169.34 W/m2,当t0取0 ℃时,tw=28.3 ℃.从以上分析计算可知,管径、管内流体热物性参数、流体流量、埋管深度对地源热泵融雪化冰效果有很大影响,如果设计不合理,会影响到系统的可靠性、经济性.4不同气候条件下的热流密度上述分析计算只是给定气象条件下的热流密度,不同天气条件下热流密度是不相同的.由以上分析知热流密度与降雪量、空气温度、风速和环境相对湿度有关,而这些因素孰主孰次从以上的计算还无法推段出,因此很有必要进行不同气候条件下融雪热流密度的计算.4.1不同降雪量时的热流密度空气温度为Ta=-3 ℃,风速为V=4 m/s,环境相当湿度为φ=60%,融雪等级Ar=1,太阳辐射qI=85.56 W/m2保持不变,降雪量变化时的热流密度见表1.表1不同降雪量时的热流密度W/m2
Table 1Heat flux of different snowfall W/m2
S/(mm/h)qsqmqcqrqeq011.568570.115.9940.3117.423.1217070.115.9940.3203.9634.6825570.115.9940.3290.524.2不同风速时的热流密度降雪量S=2 mm/h,空气温度Ta=-3 ℃,环境相对湿度φ=60%,融雪等级Ar=1,太阳辐射 qI=85.56 W/m2保持不变,风速变化时的热流密度见表2.表2不同风速时的热流密度W/m2
Table 2Heat flux of different wind speeds  W/m2
V/(m/s)qsqmqcqrqeq033.1217057.245.9930.5181.2943.1217070.115.9940.3203.9653.1217083.015.9945.7222.264.3不同空气温度时的热流密度降雪量S=2 mm/h,风速为V=4 m/s,环境相对湿度φ=60%,融雪等级Ar=1,太阳辐射qI=85.56 W/m2保持不变,空气温度变化时的热流密度见表3.表3不同室外空气温度时的热流密度W/m2
Table 3Heat flux of different outdoor air temperatures W/m2
Ta/℃qsqmqcqrqeq0-33.1217070.115.9940.3203.96-44.1617093.4820.5440.3248.91-55.20170116.927.4640.3280.294.4不同环境相对湿度时的热流密度空气温度为Ta=-3 ℃,降雪量S=2 mm/h,风速为V=4 m/s,融雪等级Ar=1,太阳辐射qI=85.56 W/m2保持不变,环境相对湿度发生变化时的热流密度见表4.表4不同相对湿度时的热流密度W/m2
Table 4The heat flux of different relative humidity W/m2
φqsqmqcqrqeq0503.1217070.115.9950.4213.76603.1217070.115.9940.3203.96703.1217070.115.9930.3193.96803.1217070.115.9920.2182.86903.1217070.115.9910.1179.764.5极端气象条件根据武汉地区近十年的气象数据,确定极端气象条件是S=2 mm/h,Ta=-5 ℃,V=5 m/s ,太阳辐射qI=85.56 W/m2不变.在极端气象条件下不同融雪等级的热流密度见表5.
表5极端气象条件不同融雪等级时的热流密度W/m2
Table 5Heat flux of different snowmelt level of extreme weather condition W/m2
Arqsqmqcqrqeq00.75.2170138.3527.4657.6246.030.85.2170138.3527.4657.6268.370.95.2170138.3527.4657.6290.711.05.2170138.3527.4657.6313.05由表1可知,降雪量从1 mm/h增大到2 mm/h、从2 mm/h增大到3 mm/h热流密度均增加86.56 W/m2,且最大热流密度出现在其他气象条件不变降雪量为3 mm/h情况下.由表2可知风速从3 m/s增大到4 m/s热流密度增加22.67 W/m2、从4 m/s增大到5 m/s热流密度增加18.30 W/m2.由表3可知室外温度从-3 ℃降到-4 ℃热流密度增加44.95 W/m2,从-4 ℃降到-5 ℃热流密度均增加31.33 W/m2,其他气象条件不变室外温度在-5 ℃情况下热流密度居第二位.由表4可知相对湿度每增加10%,热流密度约增加10 W/m2.可见武汉地区降雪量及室外温度是影响融雪化冰热流密度的关键因素.主要是因为雪的融解热很大,且根据式(5)在水膜温度为0 ℃情况下温度越低对流换热也越大.风速是主要因素,相对湿度是次要因素,是因为风速对换热及辐射换热有较大影响,而模型中相对湿度只对蒸发热有影响.由表5可知在极端气候条件下不同的融雪等级,热流密度密度增加不大,但其取值对系统的可靠性及经济性影响极大.5结语地源热泵融雪化冰高效、环保,在国外已经有应用,而在我国还处于起步阶段.通过对沥青混凝土路面融雪随机传热机理的研究可确定不同气象条件下的融雪热流密度,通过对埋管的传热分析,埋管的参数选择对融雪化冰的效果有很大影响.本文的研究为地源热泵融雪化冰系统的设计祭奠了基础.